Next Article in Journal
Spatial Inequalities in Access to Safe Drinking Water in an Upper-Middle-Income Country: A Multi-Scale Analysis of Brazil
Next Article in Special Issue
Fish Beta Diversity Patterns across Environmental Gradients in 63 European Shallow Lakes: Effects of Turbidity, Nutrient Enrichment, and Exotic Species
Previous Article in Journal
Effects of the Eating Habits of Romanian Residents on the Water Footprint
Previous Article in Special Issue
Diet and Trophic Structure of the Fish Community in a Small Sub-Tropical Lake in Central Mexico
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Effects of Climate Change on the Habitat Suitability and Distribution of Endemic Freshwater Fish Species in Semi-Arid Central Anatolian Ecoregion in Türkiye

by
Mustafa Korkmaz
1,2,*,
Fatih Mangıt
3,
İlayda Dumlupınar
4,
Mehmet Arda Çolak
5,
Mustafa Berkay Akpınar
6,
Meltem Koru
2,
Juan Pablo Pacheco
7,8,
Arely Ramírez-García
9,
Gültekin Yılmaz
2,
Cihelio Alves Amorim
1,
İbrahim Kaan Özgencil
1,
Deniz İnnal
10,
Sedat Vahdet Yerli
3,
Korhan Özkan
2,11,
Zuhal Akyürek
5,6,11,
Meryem Beklioğlu
1,11 and
Erik Jeppesen
1,8,11,12,13,14
1
Department of Biological Sciences, Middle East Technical University, Ankara 06800, Türkiye
2
Institute of Marine Sciences, Middle East Technical University, Mersin 33731, Türkiye
3
Department of Biology, Hacettepe University, Ankara 06800, Türkiye
4
Graduate School of Science and Engineering, Hacettepe University, Ankara 06800, Türkiye
5
Department of Geodetic and Geographic Information Technologies, Middle East Technical University, Ankara 06800, Türkiye
6
Department of Civil Engineering, Middle East Technical University, Ankara 06800, Türkiye
7
Centro Universitario de la Regional del Este (CURE), Universidad de la República, Maldonado 11200, Uruguay
8
Department of Ecoscience, Aarhus University, 8000 Aarhus, Denmark
9
Laboratorio de Biología Acuática, Facultad de Biología, Universidad Michoacana de San Nicolás de Hidalgo, Morelia 58000, Mexico
10
Department of Biology, Burdur Mehmet Akif Ersoy University, Burdur 15030, Türkiye
11
Centre for Ecosystem Research and Implementation (EKOSAM), Middle East Technical University, Ankara 06800, Türkiye
12
Arctic Research Centre (ARC), Aarhus University, 8000 Aarhus, Denmark
13
Sino-Danish Centre for Education and Research (SDC), Beijing 101408, China
14
Institute for Ecological Research and Pollution Control of Plateau Lakes, School of Ecology and Environmental Sciences, Yunnan University, Kunming 650500, China
*
Author to whom correspondence should be addressed.
Water 2023, 15(8), 1619; https://doi.org/10.3390/w15081619
Submission received: 16 March 2023 / Revised: 5 April 2023 / Accepted: 18 April 2023 / Published: 21 April 2023

Abstract

:
Climate change is affecting freshwater ecosystems globally, particularly those in semi-arid and arid regions. The Central Anatolian Ecoregion (CAE) in Türkiye has a semi-arid climate and is home to numerous endemic fish species. We used species distribution modelling to elucidate the distribution of sixteen endemic fish species in CAE and predicted their potential distributions for 2041–2060 and 2081–2100 based on the CMIP6 climate model. Half of the species are predicted to experience a significant loss of climatically suitable areas. Anatolichthys fontinalis, Gobio gymnostethus, Gobio hettitorum, and Pseudophoxinus burduricus will face a complete loss of suitable areas by 2081–2100 under a high emissions climate scenario, whereas Cobitis bilseli, Egirdira nigra, Gobio intermedius, and Squalius anatolicus will experience a significant loss. The other eight species can potentially benefit from climate warming if all other stressors remain equal. Anthropogenic stressors, such as water abstraction for irrigation, pollution, invasive species introductions, and dam construction, are already putting endemic fish populations in CAE under extreme pressure. Climate change is expected to exacerbate these threats. Regular monitoring of freshwater ecosystems and fish fauna in the CAE and protecting the region from key anthropogenic stressors are recommended to successfully conserve these endemic freshwater fishes under climate change.

1. Introduction

Freshwater ecosystems are among the most endangered on Earth [1,2]. Despite only making up 2.3% of the Earth’s surface, they support high biodiversity, ecosystem functions, and services that underpin human well-being [3,4,5]. They are undergoing drastic changes globally due, e.g., to changes in land and water use, water abstraction, pollution, and colonisation by invasive species [4]. Climate change further promotes or aggravates these pressures as it increases water extraction, alters hydrological regimes, amplifies the effects of changes in land use, and affects the habitat and dispersion of species [3,6,7]. According to predictions, even slight warming scenarios pose severe threats to global biodiversity, from cellular to population levels [8]. Changing environmental conditions might favour specific species or induce changes in their distribution areas [9,10].
In the context of climate change, species are expected to exhibit one of three responses: They adapt to the altered climate, track favourable habitats, or become extinct [11]. Climate change is likely to have adverse effects on most species, independent of other underlying anthropogenic factors; however, some species may experience benefits when suitable habitat availability increases. Defining the winners and losers of climate change is emerging as a key subject in climate change biodiversity research [12,13]. Once the winners and losers have been identified, it may be possible to develop management plans targeting each species according to its unique response to climate change.
Freshwater fish comprise 40% of all fish species and 25% of all vertebrates on Earth [4,14], and they are among the most diverse taxonomic groups threatened by global changes [6,15,16]. Several studies have revealed that the phenology and distribution of freshwater fish species have already changed due to the climate change [17,18], and drastic declines in population size and their distribution range have occurred in recent decades reflecting an increase in various threats, including habitat loss and degradation, water abstraction, invasive species, overfishing, water pollution, and climate change [7,19,20]. To estimate the consequences of such habitat alterations for fish populations, it is necessary to understand how these threats affect the niche and spatial distribution of the species.
Species distribution modelling (SDM) combines the principles of ecology, biogeography, and statistical and machine learning techniques to predict the spatial distribution of species based on their environmental preferences. SDM provides insights into species distribution–environment relationships and can be used to estimate the bioclimatic niche of a species by correlating species occurrence or abundance records with climatic data [21,22]. SDM has been widely used in ecology, biogeography, conservation biology, and wildlife management as a tool to predict the potential distributions of a species using projected scenarios based on the likelihood of the existence of a targeted species in response to various environmental factors [22,23,24,25]. SDM has also proven to be useful in predicting how species may respond to changes in climate conditions [26]. It represents a practical approach for assessing the vulnerability of a species to climate change and determining its geographical distribution and fundamental niches [27]. Climate projections indicate that (semi) arid regions will be among those most affected by climate change and other anthropogenic impacts that add to habitat alteration. It is, therefore, essential to understand how climate change can affect the distribution of species, particularly those that are more sensitive, such as species endemic to these ecosystems.
There are 25 water basins in Türkiye with more than 420 fish species (including some that still need confirmation) [28,29]. According to the zoogeographical delineation conducted in the Freshwater Ecoregions of the World, Türkiye has 14 freshwater ecoregions [30], among which 3 are specific to Türkiye: the Central Anatolia, Northern Anatolia, and Lake Van ecoregions. Despite the diverse aquatic systems, Anatolia does not have abundant water resources. According to World Resources Institute data, Türkiye is ranked 32 on the water-stress countries list [31]. The Central Anatolia Ecoregion (CAE), which includes the Konya Closed, Burdur, and Akarçay River Basins, is expected to be significantly affected by climate change [32,33,34]. A significant portion of the CAE can be characterised as having a semi-arid climate [35,36] and a high level of endemism [37]. Geologically, CAE is the remnant of the Central Anatolia Lake System, which showed significant fluctuations through time [38,39,40] and was inhabited by fishes from central and western Europe, Mesopotamia, and central and western Asia [28,38]. Many endemic fish species in the region are already threatened [32,41,42], and climate change is likely to intensify these threats [43,44]. The current state of the endemic freshwater fish in the region is alarming, with 30 out of the 43 endemic freshwater fish taxa recognised as being threatened [29,32,33].
Here, we elucidate the present habitat suitability and forecast habitat changes for endemic freshwater fish in the Central Anatolia Ecoregion using SDM and various climate change scenarios. We aimed to determine how climate change affects endemic freshwater fish species in the region and to identify winners and losers in the species distribution under climate change scenarios.

2. Materials and Methods

2.1. Study Area

The CAE is a 1000 m high plateau consisting of three major endorheic water basins; Konya Closed, Burdur, and Akarçay Basins, as well as Lake Eğirdir, which is connected to the Antalya Basin and the Mediterranean Sea only historically through karst aquifers [45]. These basin systems have changed considerably in recent years [32,33].
Due to its geological history and geographical structure, Türkiye is known for its high endemism and genetic diversity. The CAE is among the highest contributors to this diversity, with 43 endemic fish species (Table S1 in the Supplementary Materials) [28,29,32,37,41,44]. The region’s biodiversity is still incomplete, as new species are regularly being discovered. Many endemic species are distributed in limited areas in central Anatolia’s aquatic systems that contain numerous important water bodies and watercourses (Figure 1) but are under serious threats.

2.2. Fish Presence Data

In this study, we focused on the endemic freshwater fish species in CAE using data from various data sources (Figure 1, Table S1 in the Supplementary Materials), including the Global Biodiversity Information Facility (GBIF) [46], FishBase [47], Eschmeyer’s Catalog of Fishes [48], literature [44,49,50,51,52,53], and data from field studies performed between 2012 and 2022 by some of the authors of this article. To reduce the effects of spatial sampling biases [54], the records for each species present were spatially filtered by reducing multiple records within 1 km to a single record with the “spThin” package [55] in R software (version 4.1.3) [56]. In the analysis, we assumed that fish could disperse freely to track suitable climatic conditions. We performed modelling on 16 of these species within CAE, while 15 species were excluded because their distribution extended beyond the CAE. Additionally, another nine species were excluded due to limited distribution data, and a further three species were excluded due to their extremely limited distribution range. Such exclusions were essential for maintaining the accuracy and reliability of our results.

2.3. Hydro-Climate Data

Long-term mean annual temperature and total annual precipitation data for the period 1970–2021 were obtained from the ERA5-Land product. ERA5-Land [57] is a climate reanalysis product generated by ECMWF, and it is the successor of the ERA5 [58] with a long temporal coverage from 1950 to the present (the ERA5-Land data released so far covers the period from 1981 to 2–3 months before the present). ERA5-Land provides hourly data on various atmospheric, land-surface, and sea-state parameters together with estimates of uncertainty with its enhanced 9 km spatial resolution. Surface and subsurface runoff values are used to observe the spatial distribution of runoff within the study area. In the ERA5-Land product, subsurface water fluxes are determined by Darcy’s law, and the surface runoff is obtained by the Hortonian runoff formulation [59].
Trend analysis was conducted by using generalized linear models (GLM) [60]. We modelled the temperature and precipitation response variables by defining year as the sole explanatory variable. The models were built by using the glmmTMB function from the “glmmTMB” package [61] in R software (version 4.1.3) [56]. We validated the models by following Zuur and Ienao [62] and checked for the residuals autocorrelation by using autocorrelation and partial autocorrelation plots.

2.4. Environmental Data

The present and future bioclimatic data were obtained from the WorldClim database version 2.1, having a spatial resolution of 30 arc-seconds [63]. For the 2041–2060 and 2081–2100 periods, we used data from the Coupled Model Intercomparison Project Phase 6—CMIP6 [64]. Downscaling and calibration were conducted with the WorldClim database version 2.1 as the baseline climate and based on the Model for Interdisciplinary Research on Climate Version 6 (MIROC6), a global climate model, and two Shared Socio-economic Pathways (SSP2-4.5 and SSP5-8.5) were used. The Shared Socio-economic Pathways (SSP) framework comprises a set of scenarios that outline possible trajectories of future socio-economic development and greenhouse gas emissions. Specifically, SSP2-4.5 envisions a future where greenhouse gas emissions peak around 2040 and subsequently decline. In contrast, SSP5-8.5 projects a future where greenhouse gas emissions continue to increase throughout the century [65]. Nineteen scenopoetic bioclimatic variables [66] were derived from monthly temperature and precipitation values. To ensure that the selected variables were not highly correlated, multicollinearity was assessed using variance inflation factor (VIF) testing, with a threshold of 10 for VIF and 80% for correlation using the “usdm” package (version 1.1-18; [67]) in shinyBIOMOD [68] and the “vifstep” and “vifcor” functions in R software (version 4.1.3) [56].

2.5. Species Distribution Modelling

Maxent v3.4.1 [69] software was used to model the endemic fish species and infer their bioclimatic suitability in the study area in both present and future scenarios. Maxent is a machine learning algorithm that predicts niches by correlating species presence with environmental variables, making it suitable for presence-only data [69]. It iteratively minimises the relative entropy between the probability densities at presence sites and the probability densities of the landscape, and it has been shown to be a useful tool for developing successful distribution models, particularly for small sample sizes [70,71].
The optimal variables and model settings were selected [70,72] in the WALLACE v1.0.6.3 platform [73,74], which is an R-scripted modern workflow. A user-specified study region was created from the presence records to select the background, and 1000 background points were used in the model. Candidate models were tested by combining five combinations of feature classes (linear; linear and quadratic; hinge; linear, quadratic, and hinge; linear, quadratic, hinge, and product) and five values of regularisation multiplier (1 to 5 in increments of 1).
The modelling process involved determining the model settings and developing a final model for each species using the entire dataset, including presence records and 1000 background points, projected onto present and future bioclimatic conditions for the study area. The validation of the models was based on the area under the receiver operating characteristic curve (AUC), which measures the ability of the model to distinguish between presence and absence points. AUC values below 0.5 indicate that the model’s performance is not better than random variation [75], while values above 0.7 indicate a good fit, and values above 0.9 indicate an excellent fit.
To identify extrapolation risks in model transfers, multivariate environmental similarity surface (MESS) was used [76,77]. Response curves were used to evaluate the impact of individual variables on model predictions, while a “jackknife” test was used to determine the relative contribution of each variable to the final model. To create maps of habitat suitability, “cloglog” output was used, with values ranging from 0 (low suitability) to 1 (high suitability). Model performance was evaluated using a partial ROC analysis with “ntbox” (NicheToolBox) software [78] (Proportion of omission = 0.05, Percentage of random points = 50, Bootstrap iterations = 1000). Maxent model outputs were converted into binary maps using “SDMtoolbox” software (version 2.5) [79], and calculation of the areas was conducted with the “biomod2” package in R software (version 4.1.3) [56].

3. Results

3.1. Model Performance

We found that the models for all 16 endemic species performed well based on the area under the receiver operating characteristic curve (AUC); all were above 0.85 (Table 1).

3.2. Importance of Climate Variables

The VIF analysis eliminated 13 of the 19 climate variables (Table S2 in the supplementary materials), and the remaining variables used as model input were annual mean temperature (Bio1), isothermality (Bio3), temperature seasonality (Bio4), mean temperature of wettest quarter (Bio8), annual precipitation (Bio12), and precipitation seasonality (Bio15).
Annual precipitation was the most critical influential factor as it determined the distribution of the selected fish species, with 7 of the 16 species exhibiting a solid association with specific annual precipitation. Furthermore, annual mean temperature, temperature seasonality, and precipitation seasonality for the remaining nine species emerged as critical predictors of their distributions with 43 endemic fish species.

3.3. Climatic Scenarios

In the scenario analysis, we found that five species (Anatolichthys fontinalis, Cobitis bilseli, Gobio gymnostethus, Gobio hettitorum, and Pseudophoxinus burduricus) are likely to experience a reduction in their predicted distribution range in both periods (2041–2060 and 2081–2100) and scenarios (SSP2-4.5 and SSP5-8.5) due to climate change (Figure 2; Table 2). Additionally, according to the predictions, Egirdira nigra is expected to undergo a reduction in the 2081–2100 period in both scenarios, and the 2041–2060 period in the SSP5-8.5. Gobio intermedius is expected to experience a range reduction in the 2081–2100 period under both scenarios. Furthermore, Squalius anatolicus is predicted to undergo a range reduction during the 2081–2100 period, specifically under the SSP5-8.5 scenario (Figure 3). Eight species (Anatolichthys anatoliae, Anatolichthys iconii, Chondrostoma beysehirense, Capoeta pestai, Oxynoemacheilus nasreddini, Pseudophoxinus anatolicus, Pseudophoxinus crassus, and Squalius recurvirostris; Table 2; Figures S1 and S2 in the Supplementary Materials), however, are predicted to have a potential distribution range increase.

3.4. Hydro-Climate Trends

The model we built for the trend in mean annual temperature yielded an effect size estimate of 0.056 (gamma GLM with log link, standard error = 0.001), and the one we built for the trend in total annual precipitation yielded an effect size estimate of −0.032 (gamma GLM with log link, standard error = 0.018).
During the 52-year study period (1970–2021), the increasing trend in mean annual temperature indicated by ERA5-Land was 0.041 °C/year, and the annual total precipitation decreased by 50 mm, from 465 mm to 416 mm (Figure 4). According to the ERA5-Land reanalysis dataset, total annual runoff in the study area showed a significant decreasing trend (p < 0.05) for May in all of the study areas and for August in the major part of the study area (Figure 5).

4. Discussion

We predicted the impact of climate change on the distribution of 16 out of 43 endemic fish species in the CAE that were suitable for modelling using SDM. We found that these species exhibited distinct responses to climate change, with eight species predicted to lose their climatically suitable habitats and four being at risk of complete extinction, while eight may potentially expand their climatically suitable habitats. Our findings of loss of climatically suitable habitats align with prior studies assessing the potential effects of climate change on freshwater fish fauna using SDM. For instance, Stewart et al. [81] projected that, by 2080, 9 of the 13 native fish taxa in southwestern Australia would lose their suitable habitat. Similarly, Yousefi et al. [12] predicted that 15 endemic freshwater fish species in Iran would experience range reductions, with 5 endemic species facing significant losses. Frederico et al. [82] suggested that the Amazon basin may lose 2% of its freshwater fish fauna due to unsuitable climatic conditions, with a further 34% expected to be adversely affected. Our finding is consistent with a study by Yılmaz et al. [32], revealing that, of the 38 fish species occurring in the Konya Closed Basin (one of the basins in CAE), 61% are threatened and highly vulnerable to the effects of climate change and water use, including a decline in water level.
Specifically, the models predict that A. fontinalis, G. gymnostethus, G. hettitorum, and P. burduricus will experience a complete loss of climatically suitable area by 2081–2100 under a high emissions scenario (SSP5-8.5). This potential loss highlights the extreme consequences of climate change in semi-arid regions, considering that the current conservation status for these species are NE, CR, CR, and EN, respectively [80]. Additionally, C. bilseli, E. nigra, G. intermedius, and S. anatolicus are predicted to have a significant loss of climatically suitable areas.
Of the 22 Squalius species in Türkiye, 4 are distributed in CAE. Members of this genus prefer pelagic areas and can be found in streams and lakes. They are considered tolerant to environmental alterations [80]. The model predicts that S. anatolicus may potentially expand its distribution by 2041–2060 but lose climatically suitable habitat by 2081–2100, and S. cappadocicus, a critically endangered species according to IUCN [80], is expected to lose its climatically suitable areas by 2081–2100 under the high emissions scenario (SSP5-8.5). The Gobio genus consists of benthic and rheophilic species. They generally prefer streams, but some species can also be found in lakes. G. gymnostethus and G. hettitorum are predicted to lose their suitable habitats completely by the end of this century. However, climatically suitable areas for G. intermedius showed a slight increase by 2041–2060 in both scenarios (SSP2-4.5 and SSP5-8.5), followed by a significant decrease by 2081–2100 in both scenarios.
Anatolia is a vital diversification centre for the genus Pseudophoxinus [83], and 24 species appear in Anatolia, accounting for nearly 90% (except 3 species) of the known species, and P. burduricus is endemic to the Burdur Basin. Even though the current distribution model defines Beyşehir and Akarçay basins as climatically suitable areas, it is not possible for this species to naturally disperse to these basins. Therefore, even a predicted increase in suitable habitats can only be functional if human management intervention, such as assisted migrations, is anticipated [84]. Based on the SDM results for C. bilseli, the species will experience a reduction in its climatically suitable habitats in both scenarios and periods. C. bilseli primarily inhabits small rivers and streams in Beyşehir lake tributaries. Despite its limited distribution, the species faces multiple threats, such as widespread pollution, excessive water abstraction, and the introduction of non-native fish species [44,85]. The distribution areas of E. nigra are predicted to decrease in both scenarios by 2081–2100, while an increase in climatically suitable habitats is projected by 2081–2100 under the SSP2-4.5 scenario. E. nigra is facing threats from water abstraction, pollution, and the introduction of non-native fish species, and E. nigra was previously distributed throughout the Lake Eğirdir Basin, including the lake; however, the introduction of the predatory non-native pikeperch Sander lucioperca (Linnaeus 1758) caused a significant decline in population size as well as significant range reduction [44].
The SDM identified potentially climatically suitable distribution areas for four endangered (EN) species of the eight that are predicted to experience an expansion in climatically suitable habitats: Capoeta pestai, Chondrostoma beysehirense, Pseudophoxinus anatolicus, and Pseudophoxinus crassus. These are the predicted winners of climate change, and the combined impact of climate change and other stressors may have unforeseen consequences due to the model uncertainty and other limitations. However, multiple other stressors threaten these species, including invasive species, water abstraction, drought, pollution, and dam construction [41,44,80,86,87,88]. For C. pestai, introduction of the non-native predatory pikeperch is one of the primary concerns [89]. C. beysehirense, which is found in the Beyşehir basin, is also mainly threatened by the introduction of the non-native predatory pikeperch [90]. P. anatolicus is under threat from water abstraction for irrigation, pollution, overfishing, restricted water levels, and non-native pikeperch introduction [91], and two of the three drainage basins where it occurred were already drained for agriculture [44]. Finally, P. crassus is limited to Central Anatolia, specifically İnsuyu Creek and Lake Gökgöl, located west of the Lake Tuz basin, and it is threatened by excessive water abstraction and the construction of dams and weirs [80].
Estimating the potential climatic distribution area of freshwater fish is not as straightforward as for terrestrial elements such as birds, reptiles, and mammals, and there is an ongoing debate about using climatic variables as surrogates for instream conditions [92]. For example, air temperature and precipitation do not directly reflect water temperature and hydrology. Still, some studies have shown that the inclusion of instream parameters (unavailable in our case) does not necessarily improve model accuracy [22]. Representativeness of climatic variables is one thing, but the capacity to settle in a suitable area is quite another. The dispersal of freshwater fish almost entirely depends on suitable freshwater conditions and connectivity. In addition to temperature, environmental factors such as flow (dams, seasonality), habitat, and pollution (physicochemical properties) are important [72] as they may prevent species from reaching new areas [76,93,94]. Biotic interactions (local fauna, predation, population/colonisation density) also affect dispersal capacity [95].
SDM use bioclimatic variables and presence/absence data to predict species distribution. However, correlations between variables do not always imply causality, and the model accuracy has limitations, as it is not accounting for other influential factors such as species adaptations, habitat fragmentation, biotic interactions, and anthropogenic stressors. In addition, estimating species distribution with large-scale predictors, as in this study (climate data), is known to pose a risk of overestimating distributions [96,97]. However, due to the complexity of the distribution process, the requirement for detailed data [92] and the significance of climate for species distributions [98], SDM frequently builds on climate variables alone [99,100]. Furthermore, the large individual grid cells of all climate data sets might not be the best indicators of the occurrence of patches with small-scale (micro) climate conditions that are favourable for a species [101,102]. Notwithstanding the possible drawbacks and limitations, SDM has emerged as a crucial instrument in conservation management for forecasting the dispersion patterns of species. To enable a comparison of how variations in a crucial class of variables might affect model performance, we concentrated on climate predictors.
The ERA5-Land reanalysis datasets revealed a highly significant increasing trend for temperature and a decreasing trend for precipitation for the period 1970–2021. Thus, surface and subsurface runoff values estimated by ERA5-Land have a decreasing trend within the CAE in most of the months. The decreasing trend is significant in May and August around the locations of the studied species. According to trend analysis, an increase in evaporation and transpiration has resulted in water loss in the CAE basin after 2000, indicating a water deficit [32]. The primary driver of water loss in lakes, however, is unsustainable water use for agriculture in the Konya and Burdur Basins, which are the main basins of CAE. Changes in land use and irrigation have detrimental effects on semi-arid and arid regions worldwide, including the Konya Closed and Burdur Basins [32,33]. The uncontrolled use of water for irrigation in the basins, leading to water abstractions, poses one of the primary threats to endemic fish species, irrespective of the impact of climate change. Water abstraction effects directly impact the winners and losers identified in this study. Therefore, species identified as losers may actually lose their predicted climatically suitable distribution areas more rapidly than predicted by the SDM. Furthermore, the species identified as winners may not attain the predicted gains and may potentially become losers if the water management strategy is not changed in the direction of a more sustainable use of water. Thus, our predictions based on the SDM may only be used to provide a conservative prediction of the future for these endemic freshwater fish species.

5. Conclusions

The Central Anatolian Ecoregion (CAE) is expected to experience increased temperature and decreased precipitation due to climate change, which is similar to other arid and semi-arid regions globally [103]. These changes and extensive land and water use in the region are likely to profoundly impact the flow regimes of freshwater systems, making the freshwater fauna highly vulnerable [18]. Our analyses identified 8 out of 16 analysed species in the CAE that will experience a loss of suitable habitats under climate change, making them highly vulnerable. By the end of this century, four of these species are projected to lose their climatically suitable habitats entirely. Moreover, human activities, including but not limited to the destruction of habitats, availability of surface water, and an increase in water usage, are expected to exacerbate the adverse effects of climate change in the upcoming decades. It is plausible that the unconsidered variables in our study could potentially influence the outcomes of this study. Consequently, our findings may underestimate the negative impacts of climate change.
The CAE endemic freshwater fish fauna is currently facing a high level of threat, with 65% of endemic species considered threatened or nearly threatened by the IUCN Red List, and two species are already extinct. Climate change is expected to exacerbate the effects of other anthropogenic stressors, such as water abstraction, irrigation, pollution, the introduction of invasive species, and dam construction, particularly in semi-arid and arid areas and, more specifically, in the CAE, where fish communities are already under extreme pressure. The results of this study indicate that protection of the CAE water resources is essential for preserving the endemic freshwater fish fauna in Türkiye.
In order to preserve freshwater fish populations, appropriate conservation measures need to be implemented. A policy framework is urgently needed to restrict the exploitation of water resources to sustainable limits while simultaneously promoting conservation efforts to reverse ecosystem degradation or maintain current conditions. In this context, a better understanding of the response of endemic fish species to the combined effect of changes in climate and other stressors is needed. Regular monitoring of all freshwater ecosystems in the CAE is recommended, as well as protection of the region from other anthropogenic stressors to effectively conserve freshwater fishes under climate change. Regular monitoring programs can detect the presence of new species or the absence of previously recorded species, track changes in species distribution, richness, and composition, and support conservation planning for freshwater fishes in the face of climate change.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/w15081619/s1, Figure S1: Current and future distribution of Anatolichthys anatoliae, Anatolichthys iconii, Chondrostoma beysehirense, and Capoeta pestai in the CAE. Colours show habitat suitability, and the value of suitability increases in red areas; Figure S2: Current and future distribution of Oxynoemacheilus nasreddini, Pseudophoxinus anatolicus, Pseudophoxinus crassus, and Squalius recurvirostris in the CAE. Colours show habitat suitability, and the value of suitability increases in red areas; Table S1: Endemic fish fauna of the Central Anatolian Ecoregion and their IUCN conservation status; Table S2: The environmental predictors used in the SDM with their variance inflation factors (VIF) after the exclusion of collinear variables. References [104,105,106,107,108,109,110,111] are cited in the supplementary materials.

Author Contributions

Conceptualization: M.K. (Mustafa Korkmaz), E.J., F.M., İ.D., S.V.Y., D.İ., K.Ö., M.B. and Z.A.; methodology, validation, formal analysis, and investigation: İ.D., M.K. (Mustafa Korkmaz), F.M., M.B.A., Z.A. and E.J.; data curation: M.K. (Mustafa Korkmaz), F.M., İ.D., D.İ., M.B.A. and Z.A.; writing—original draft preparation: M.K. (Mustafa Korkmaz), F.M., İ.D., M.A.Ç., M.B.A., M.K. (Meltem Koru), J.P.P., A.R.-G., G.Y., C.A.A., İ.K.Ö., D.İ., S.V.Y., K.Ö., Z.A., M.B. and E.J.; writing—review and editing: M.K. (Mustafa Korkmaz), F.M., İ.D., M.A.Ç., M.B.A., M.K. (Meltem Koru), J.P.P., A.R.-G., G.Y., C.A.A., İ.K.Ö., D.İ., S.V.Y., K.Ö., Z.A., M.B. and E.J.; visualization: İ.D., M.K. (Mustafa Korkmaz), F.M., M.K. (Meltem Koru), M.A.Ç. and İ.K.Ö.; funding acquisition: E.J. All authors have read and agreed to the published version of the manuscript.

Funding

The project and M.K. (Mustafa Korkmaz), M.A.Ç., G.Y., C.A.A., K.Ö., M.B., Z.A. and E.J. were supported by TÜBİTAK program BIDEB2232, project number 118C250.

Data Availability Statement

The data are available from the authors on reasonable request.

Acknowledgments

We are grateful to Anne Mette Poulsen for linguistic improvements to the manuscript. We acknowledge the World Climate Research Programme, which, through its Working Group on Coupled Modelling, coordinated and promoted CMIP6. We thank the climate modelling groups for producing and making available their model output, the Earth System Grid Federation (ESGF) for archiving the data and providing access, and the multiple funding agencies who support CMIP6 and ESGF.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Nel, J.L.; Roux, D.J.; Abell, R.; Ashton, P.J.; Cowling, R.M.; Higgins, J.V.; Thieme, M.; Viers, J.H. Progress and Challenges in Freshwater Conservation Planning. Aquat. Conserv. 2009, 19, 474–485. [Google Scholar] [CrossRef]
  2. Döll, P.; Zhang, J. Impact of Climate Change on Freshwater Ecosystems: A Global-Scale Analysis of Ecologically Relevant River Flow Alterations. Hydrol. Earth Syst. Sci. 2010, 14, 783–799. [Google Scholar] [CrossRef]
  3. Reid, A.J.; Carlson, A.K.; Creed, I.F.; Eliason, E.J.; Gell, P.A.; Johnson, P.T.J.; Kidd, K.A.; MacCormack, T.J.; Olden, J.D.; Ormerod, S.J.; et al. Emerging Threats and Persistent Conservation Challenges for Freshwater Biodiversity. Biol. Rev. 2019, 94, 849–873. [Google Scholar] [CrossRef] [PubMed]
  4. Dudgeon, D.; Arthington, A.H.; Gessner, M.O.; Kawabata, Z.; Knowler, D.J.; Lévêque, C.; Naiman, R.J.; Prieur-Richard, A.H.; Soto, D.; Stiassny, M.L.J.; et al. Freshwater Biodiversity: Importance, Threats, Status and Conservation Challenges. Biol. Rev. Camb. Philos. 2006, 81, 163–182. [Google Scholar] [CrossRef] [PubMed]
  5. Dawson, M.N. Species Richness, Habitable Volume, and Species Densities in Freshwater, the Sea, and on Land. Front. Biogeogr. 2012, 4, 105–116. [Google Scholar] [CrossRef]
  6. Janse, J.H.; Kuiper, J.J.; Weijters, M.J.; Westerbeek, E.P.; Jeuken, M.H.J.L.; Bakkenes, M.; Alkemade, R.; Mooij, W.M.; Verhoeven, J.T.A. GLOBIO-Aquatic, a Global Model of Human Impact on the Biodiversity of Inland Aquatic Ecosystems. Environ. Sci. Policy 2015, 48, 99–114. [Google Scholar] [CrossRef]
  7. Jeppesen, E.; Brucet, S.; Naselli-Flores, L.; Papastergiadou, E.; Stefanidis, K.; Nõges, T.; Nõges, P.; Attayde, J.L.; Zohary, T.; Coppens, J.; et al. Ecological Impacts of Global Warming and Water Abstraction on Lakes and Reservoirs Due to Changes in Water Level and Related Changes in Salinity. Hydrobiologia 2015, 750, 201–227. [Google Scholar] [CrossRef]
  8. IPCC. IPCC. IPCC Summary for Policymakers. In Climate Change 2021: The Physical Science Basis. Contribution of Working Group I to the Sixth Assessment Report of the Intergovernmental Panel on Climate Change; Masson-Delmotte, V., Zhai, P., Pirani, A., Connors, S.L., Péan, C., Berger, S., Caud, N., Chen, Y., Goldfarb, L., Gomis, M.I., et al., Eds.; IPCC: Geneva, Switzerland, 2021; ISBN 978-92-9169-158-6. [Google Scholar]
  9. Munday, P.L.; Jarrold, M.D.; Nagelkerken, I. Ecological Effects of Elevated CO2 on Marine and Freshwater Fishes: From Individual to Community Effects. Fish Physiol. 2019, 37, 323–368. [Google Scholar] [CrossRef]
  10. Parmesan, C. Ecological and Evolutionary Responses to Recent Climate Change. Annu. Rev. Ecol. Evol. Syst. 2006, 37, 637–669. [Google Scholar] [CrossRef]
  11. Hannah, L. Climate Change Biology; Academic Press: New York, NY, USA, 2021. [Google Scholar]
  12. Yousefi, M.; Jouladeh-Roudbar, A.; Kafash, A. Using Endemic Freshwater Fishes as Proxies of Their Ecosystems to Identify High Priority Rivers for Conservation under Climate Change. Ecol. Indic. 2020, 112, 106137. [Google Scholar] [CrossRef]
  13. Muths, E.; Chambert, T.; Schmidt, B.R.; Miller, D.A.W.; Hossack, B.R.; Joly, P.; Grolet, O.; Green, D.M.; Pilliod, D.S.; Cheylan, M.; et al. Heterogeneous Responses of Temperate-Zone Amphibian Populations to Climate Change Complicates Conservation Planning. Sci. Rep. 2017, 7, 17102. [Google Scholar] [CrossRef]
  14. Tedesco, P.A.; Oberdorff, T.; Cornu, J.-F.; Beauchard, O.; Brosse, S.; Dürr, H.H.; Grenouillet, G.; Leprieur, F.; Tisseuil, C.; Zaiss, R.; et al. A Scenario for Impacts of Water Availability Loss Due to Climate Change on Riverine Fish Extinction Rates. J. Appl. Ecol. 2013, 50, 1105–1115. [Google Scholar] [CrossRef]
  15. Carpenter, S.R.; Stanley, E.H.; Vander Zanden, M.J. State of the World’s Freshwater Ecosystems: Physical, Chemical, and Biological Changes. Annu. Rev. Env. Resour. 2011, 36, 75–99. [Google Scholar] [CrossRef]
  16. Dudgeon, D. Asian River Fishes in the Anthropocene: Threats and Conservation Challenges in an Era of Rapid Environmental Change. J. Fish Biol. 2011, 79, 1487–1524. [Google Scholar] [CrossRef] [PubMed]
  17. Arthington, A.H.; Dulvy, N.K.; Gladstone, W.; Winfield, I.J. Fish Conservation in Freshwater and Marine Realms: Status, Threats and Management. Aquat. Conserv. 2016, 26, 838–857. [Google Scholar] [CrossRef]
  18. Comte, L.; Olden, J.D. Climatic Vulnerability of the World’s Freshwater and Marine Fishes. Nat. Clim. Chang. 2017, 7, 718–722. [Google Scholar] [CrossRef]
  19. Jeppesen, E.; Mehner, T.; Winfield, I.J.; Kangur, K.; Sarvala, J.; Gerdeaux, D.; Rask, M.; Malmquist, H.J.; Holmgren, K.; Volta, P.; et al. Impacts of Climate Warming on the Long-Term Dynamics of Key Fish Species in 24 European Lakes. Hydrobiologia 2012, 694, 1–39. [Google Scholar] [CrossRef]
  20. Barbarossa, V.; Bosmans, J.; Wanders, N.; King, H.; Bierkens, M.F.P.; Huijbregts, M.A.J.; Schipper, A.M. Threats of Global Warming to the World’s Freshwater Fishes. Nat. Commun. 2021, 12, 1701. [Google Scholar] [CrossRef]
  21. Pearson, R.G.; Dawson, T.P. Predicting the Impacts of Climate Change on the Distribution of Species: Are Bioclimate Envelope Models Useful? Glob. Ecol. Biogeogr. 2003, 12, 361–371. [Google Scholar] [CrossRef]
  22. Elith, J.; Leathwick, J.R. Species Distribution Models: Ecological Explanation and Prediction Across Space and Time. Annu. Rev. Ecol. Evol. Syst. 2009, 40, 677–697. [Google Scholar] [CrossRef]
  23. Cianfrani, C.; Broennimann, O.; Loy, A.; Guisan, A. More than Range Exposure: Global Otter Vulnerability to Climate Change. Biol. Conserv. 2018, 221, 103–113. [Google Scholar] [CrossRef]
  24. Marcer, A.; Sáez, L.; Molowny-Horas, R.; Pons, X.; Pino, J. Using Species Distribution Modelling to Disentangle Realised versus Potential Distributions for Rare Species Conservation. Biol. Conserv. 2013, 166, 221–230. [Google Scholar] [CrossRef]
  25. Gama, M.; Crespo, D.; Dolbeth, M.; Anastácio, P.M. Ensemble Forecasting of Corbicula Fluminea Worldwide Distribution: Projections of the Impact of Climate Change. Aquat. Conserv. 2017, 27, 675–684. [Google Scholar] [CrossRef]
  26. Beaumont, L.J.; Hughes, L.; Pitman, A.J. Why Is the Choice of Future Climate Scenarios for Species Distribution Modelling Important? Ecol. Lett. 2008, 11, 1135–1146. [Google Scholar] [CrossRef]
  27. Buisson, L.; Thuiller, W.; Lek, S.; Lim, P.; Grenouillet, G. Climate Change Hastens the Turnover of Stream Fish Assemblages. Glob. Chang. Biol. 2008, 14, 2232–2248. [Google Scholar] [CrossRef]
  28. Kuru, M.; Yerli, S.V.; Mangıt, F.; Ünlü, E.; Alp, A. Fish Biodiversity in Inland Waters of Turkey. J. Acad. Doc. Fish. Aquac. 2014, 1, 93–120. [Google Scholar]
  29. Çiçek, E.; Sungur, S.; Fricke, R. Freshwater Lampreys and Fishes of Turkey; A Revised and Updated Annotated Checklist 2020. Zootaxa 2020, 4809, 241–270. [Google Scholar] [CrossRef]
  30. Abell, R.; Thieme, M.L.; Revenga, C.; Bryer, M.; Kottelat, M.; Bogutskaya, N.G.; Coad, B.W.; Mandrak, N.; Balderas, S.C.; Bussing, W.; et al. Freshwater Ecoregions of the World: A New Map of Biogeographic Units for Freshwater Biodiversity Conservation. BioScience 2008, 58, 403–414. [Google Scholar] [CrossRef]
  31. WRI World Resources Institute Aqueduct. Available online: http://www.aqueduct.wri.org (accessed on 4 February 2023).
  32. Yılmaz, G.; Çolak, M.A.; Özgencil, İ.K.; Metin, M.; Korkmaz, M.; Ertuğrul, S.; Soyluer, M.; Bucak, T.; Tavşanoğlu, Ü.N.; Özkan, K.; et al. Decadal Changes in Size, Salinity, Waterbirds, and Fish in Lakes of the Konya Closed Basin, Turkey, Associated with Climate Change and Increasing Water Abstraction for Agriculture. Inland Waters 2021, 11, 538–555. [Google Scholar] [CrossRef]
  33. Çolak, M.A.; Öztaş, B.; Özgencil, İ.K.; Soyluer, M.; Korkmaz, M.; Ramírez-garcía, A.; Metin, M.; Yılmaz, G.; Ertuğrul, S.; Tavşanoğlu, Ü.N.; et al. Increased Water Abstraction and Climate Change Have Substantial Effect on Morphometry, Salinity, and Biotic Communities in Lakes: Examples from the Semi-Arid Burdur Basin (Turkey). Water 2022, 14, 1241. [Google Scholar] [CrossRef]
  34. Yagmur, N.; Bilgilioglu, B.B.; Dervisoglu, A.; Musaoglu, N.; Tanik, A. Long and Short-Term Assessment of Surface Area Changes in Saline and Freshwater Lakes via Remote Sensing. Water Environ. J. 2021, 35, 107–122. [Google Scholar] [CrossRef]
  35. Şenkul, Ç.; Ören, A.; Doğan, U.; Eastwood, W.J. Late Holocene Environmental Changes in the Vicinity of Kültepe (Kayseri), Central Anatolia, Turkey. Quatern. Int. 2018, 486, 107–115. [Google Scholar] [CrossRef]
  36. Türkeş, M.; Öztaş, T.; Tercan, E.; Erpul, G.; Karagöz, A.; Dengiz, O.; Doğan, O.; Şahin, K.; Avcıoğlu, B. Desertification Vulnerability and Risk Assessment for Turkey via an Analytical Hierarchy Process Model. Land Degrad. Dev. 2020, 31, 205–214. [Google Scholar] [CrossRef]
  37. Çiçek, E.; Fricke, R.; Sungur, S.; Eagderi, S. Endemic Freshwater Fishes of Turkey. FishTaxa 2018, 3, 1–39. [Google Scholar]
  38. Kosswig, C. Zoogeography of the Near East. Syst. Biol. 1955, 4, 49–96. [Google Scholar] [CrossRef]
  39. Meijers, M.J.M.; Brocard, G.Y.; Whitney, D.L.; Mulch, A. Paleoenvironmental Conditions and Drainage Evolution of the Central Anatolian Lake System (Turkey) during Late Miocene to Pliocene Surface Uplift. Geosphere 2020, 16, 490–509. [Google Scholar] [CrossRef]
  40. Kuzucuoğlu, C. Geomorphological Landscapes in the Konya Plain and Surroundings. In Landscapes and Landforms of Turkey; World Geomorphological Landscapes; Springer: Cham, Switzerland, 2019; pp. 353–368. [Google Scholar]
  41. Giannetto, D.; Innal, D. Status of Endemic Freshwater Fish Fauna Inhabiting Major Lakes of Turkey under the Threats of Climate Change and Anthropogenic Disturbances: A Review. Water 2021, 13, 1534. [Google Scholar] [CrossRef]
  42. Boll, T.; Levi, E.E.; Bezirci, G.; Özuluğ, M.; Tavşanoğlu, Ü.N.; Çakıroğlu, A.İ.; Özcan, S.; Brucet, S.; Jeppesen, E.; Beklioğlu, M. Fish Assemblage and Diversity in Lakes of Western and Central Turkey: Role of Geo-Climatic and Other Environmental Variables. Hydrobiologia 2016, 771, 31–44. [Google Scholar] [CrossRef]
  43. Morrongiello, J.R.; Beatty, S.J.; Bennett, J.C.; Crook, D.A.; Ikedife, D.N.E.N.; Kennard, M.J.; Kerezsy, A.; Lintermans, M.; McNeil, D.G.; Pusey, B.J.; et al. Climate Change and Its Implications for Australia’s Freshwater Fish. Mar. Freshw. Res. 2011, 62, 1082–1098. [Google Scholar] [CrossRef]
  44. Freyhof, J.; Bergner, L.; Ford, M. Threatened Freshwater Fishes of the Mediterranean Basin Biodiversity Hotspot: Distribution, Extinction Risk and the Impact of Hydropower. Euronatur and Riverwatch; Museum für Naturkunde Berlin (MfN)–Leibniz Institute for Evolution and Biodiversity Science: Berlin, Germany, 2020. [Google Scholar] [CrossRef]
  45. Kesici, E.; Kesici, C. The Effects of Interferences in Natural Structure of Lake Egirdir (Isparta) to Ecological Disposition of the Lake. Su Ürün. Derg. 2006, 23, 99–103. [Google Scholar]
  46. Lane, M.A.; Edwards, J.L. The Global Biodiversity Information Facility (GBIF). In Biodiversity Databases; CRC Press: Boca Raton, FL, USA, 2007; ISBN 978-0-429-15233-7. [Google Scholar]
  47. Froese, R.; Pauly, D. FishBase. World Wide Web Electronic Publication. Available online: www.fishbase.org (accessed on 5 March 2023).
  48. Fricke, R.; Eschmeyer, W.N.; van der Laan, R. Eschmeyer’s Catalog of Fishes: Genera, Species, References. Available online: http://researcharchive.calacademy.org/research/ichthyology/catalog/fishcatmain.asp (accessed on 5 March 2023).
  49. Freyhof, J.; Yoğurtçuoğlu, B. A Proposal for a New Generic Structure of the Killifish Family Aphaniidae, with the Description of Aphaniops teimorii (Teleostei: Cyprinodontiformes). Zootaxa 2020, 4810, 421–451. [Google Scholar] [CrossRef] [PubMed]
  50. Kelleci, M.; Seçer, B.; ÇiÇek, E.; Sungur, S. Aksaray İli (Türkiye) İhtiyofaunası. Acta Aquat. Turc. 2021, 17, 279–289. [Google Scholar] [CrossRef]
  51. Kuyumcu, M.; Aksu, İ.; Bektaş, Y. Genetic Analysis of Aphaniidae Hoedeman, 1949 (Teleostei: Cyprinodontiformes) Family in Anatolia. J. Anatol. Environ. Anim. Sci. 2021, 6, 627–634. [Google Scholar] [CrossRef]
  52. Yoğurtçuoğlu, B.; Kaya, C.; Freyhof, J. Oxynoemacheilus nasreddini, a New Nemacheilid Loach from Central Anatolia (Teleostei: Nemacheilidae). Zootaxa 2021, 4974, 135–150. [Google Scholar] [CrossRef]
  53. Mangit, F.; Korkmaz, M.; Sü, U.; Yerli, S.V. Actual Status of Eber Lake in Terms of Fish Community Structure. J. Limnol. Fish. Res. 2017, 3, 101–106. [Google Scholar] [CrossRef]
  54. Boria, R.A.; Olson, L.E.; Goodman, S.M.; Anderson, R.P. Spatial Filtering to Reduce Sampling Bias Can Improve the Performance of Ecological Niche Models. Ecol. Model. 2014, 275, 73–77. [Google Scholar] [CrossRef]
  55. Aiello-Lammens, M.E.; Boria, R.A.; Radosavljevic, A.; Vilela, B.; Anderson, R.P. SpThin: An R Package for Spatial Thinning of Species Occurrence Records for Use in Ecological Niche Models. Ecography 2015, 38, 541–545. [Google Scholar] [CrossRef]
  56. RCoreTeam R. R Foundation for Statistical Computing, Version 4.1.3; RCoreTeam: Vienna, Austria, 2022.
  57. Muñoz-Sabater, J.; Dutra, E.; Agustí-Panareda, A.; Albergel, C.; Arduini, G.; Balsamo, G.; Boussetta, S.; Choulga, M.; Harrigan, S.; Hersbach, H.; et al. ERA5-Land: A State-of-the-Art Global Reanalysis Dataset for Land Applications. Earth Syst. Sci. Data 2021, 13, 4349–4383. [Google Scholar] [CrossRef]
  58. Hersbach, H.; Bell, B.; Berrisford, P.; Hirahara, S.; Horányi, A.; Muñoz-Sabater, J.; Nicolas, J.; Peubey, C.; Radu, R.; Schepers, D.; et al. The ERA5 Global Reanalysis. Q. J. R. Meteor. Soc. 2020, 146, 1999–2049. [Google Scholar] [CrossRef]
  59. ECMWF IFS Documentation CY47R3—Part I: Observations. Available online: https://www.ecmwf.int/en/elibrary/81268-ifs-documentation-cy47r3-part-i-observations (accessed on 15 March 2023).
  60. Zuur, A.; Ieno, E.N.; Walker, N.; Saveliev, A.A.; Smith, G.M. Mixed Effects Models and Extensions in Ecology with R.; Springer Science & Business Media: Berlin/Heidelberg, Germany, 2009; ISBN 978-0-387-87458-6. [Google Scholar]
  61. Brooks, M.E.; Kristensen, K.; Van Benthem, K.J.; Magnusson, A.; Berg, C.W.; Nielsen, A.; Skaug, H.J.; Machler, M.; Bolker, B.M. GlmmTMB Balances Speed and Flexibility Among Packages for Zero-Inflated Generalized Linear Mixed Modelling. R J. 2017, 9, 378. [Google Scholar] [CrossRef]
  62. Zuur, A.F.; Ieno, E.N. A Protocol for Conducting and Presenting Results of Regression-Type Analyses. Methods Ecol. Evol. 2016, 7, 636–645. [Google Scholar] [CrossRef]
  63. Fick, S.E.; Hijmans, R.J. WorldClim 2: New 1-Km Spatial Resolution Climate Surfaces for Global Land Areas. Int. J. Climatol. 2017, 37, 4302–4315. [Google Scholar] [CrossRef]
  64. Eyring, V.; Bony, S.; Meehl, G.A.; Senior, C.A.; Stevens, B.; Stouffer, R.J.; Taylor, K.E. Overview of the Coupled Model Intercomparison Project Phase 6 (CMIP6) Experimental Design and Organization. Geosci. Model Dev. 2016, 9, 1937–1958. [Google Scholar] [CrossRef]
  65. Riahi, K.; Van Vuuren, D.P.; Kriegler, E.; Edmons, J.; O’neill, B.C.; Fujimori, S.; Tavoni, M. The Shared Socioeconomic Pathways and Their Energy, Land Use, and Greenhouse Gas Emissions Implications: An Overview. Glob. Environ. Chang. 2017, 42, 153–168. [Google Scholar] [CrossRef]
  66. Peterson, A.T.; Soberón, J.; Pearson, R.G.; Anderson, R.P.; Martínez-Meyer, E.; Nakamura, M.; Araújo, M.B. Ecological Niches and Geographic Distributions (MPB-49); Princeton University Press: Princeton, NJ, USA, 2011; ISBN 978-0-691-13688-2. [Google Scholar]
  67. Naimi, B.; Hamm, N.A.S.; Groen, T.A.; Skidmore, A.K.; Toxopeus, A.G. Where Is Positional Uncertainty a Problem for Species Distribution Modelling? Ecography 2014, 37, 191–203. [Google Scholar] [CrossRef]
  68. Ondo, I.; Thuiller, W.; Guéguen, M.; Pironon, S. A New R Application for Modelling Species Distribution. In Proceedings of the Species On The Move 2019, Kruger National Park, South Africa, 22–26 July 2019. [Google Scholar] [CrossRef]
  69. Phillips, S.J.; Anderson, R.P.; Dudík, M.; Schapire, R.E.; Blair, M.E. Opening the Black Box: An Open-Source Release of Maxent. Ecography 2017, 40, 887–893. [Google Scholar] [CrossRef]
  70. Elith, J.; Phillips, S.J.; Hastie, T.; Dudík, M.; Chee, Y.E.; Yates, C.J. A Statistical Explanation of MaxEnt for Ecologists. Divers. Distrib. 2011, 17, 43–57. [Google Scholar] [CrossRef]
  71. Pearson, R.G.; Raxworthy, C.J.; Nakamura, M.; Townsend Peterson, A. Predicting Species Distributions from Small Numbers of Occurrence Records: A Test Case Using Cryptic Geckos in Madagascar. J. Biogeogr. 2007, 34, 102–117. [Google Scholar] [CrossRef]
  72. Merow, C.; Smith, M.J.; Silander, J.A., Jr. A Practical Guide to MaxEnt for Modeling Species’ Distributions: What It Does, and Why Inputs and Settings Matter. Ecography 2013, 36, 1058–1069. [Google Scholar] [CrossRef]
  73. Kass, J.M.; Vilela, B.; Aiello-Lammens, M.E.; Muscarella, R.; Merow, C.; Anderson, R.P. Wallace: A Flexible Platform for Reproducible Modeling of Species Niches and Distributions Built for Community Expansion. Methods Ecol. Evol. 2018, 9, 1151–1156. [Google Scholar] [CrossRef]
  74. Muscarella, R.; Galante, P.J.; Soley-Guardia, M.; Boria, R.A.; Kass, J.M.; Uriarte, M.; Anderson, R.P. ENMeval: An R Package for Conducting Spatially Independent Evaluations and Estimating Optimal Model Complexity for Maxent Ecological Niche Models. Methods Ecol. Evol. 2014, 5, 1198–1205. [Google Scholar] [CrossRef]
  75. Phillips, S.J.; Anderson, R.P.; Schapire, R.E. Maximum Entropy Modeling of Species Geographic Distributions. Ecol. Model. 2006, 190, 231–259. [Google Scholar] [CrossRef]
  76. Elith, J.; Kearney, M.; Phillips, S. The Art of Modelling Range-Shifting Species. Methods Ecol. Evol. 2010, 1, 330–342. [Google Scholar] [CrossRef]
  77. Pearce, J.; Ferrier, S. Evaluating the Predictive Performance of Habitat Models Developed Using Logistic Regression. Ecol. Model. 2000, 133, 225–245. [Google Scholar] [CrossRef]
  78. Osorio-Olvera, L.; Lira-Noriega, A.; Soberón, J.; Peterson, A.T.; Falconi, M.; Contreras-Díaz, R.G.; Martínez-Meyer, E.; Barve, V.; Barve, N. Ntbox: An r Package with Graphical User Interface for Modelling and Evaluating Multidimensional Ecological Niches. Methods Ecol. Evol. 2020, 11, 1199–1206. [Google Scholar] [CrossRef]
  79. Brown, J.L. SDMtoolbox: A Python-Based GIS Toolkit for Landscape Genetic, Biogeographic and Species Distribution Model Analyses. Methods Ecol. Evol. 2014, 5, 694–700. [Google Scholar] [CrossRef]
  80. IUCN. The IUCN Red List of Threatened Species. Version 2022-2. Available online: https://www.iucnredlist.org (accessed on 5 March 2023).
  81. Stewart, B.A.; Ford, B.M.; Benson, J.A. Using Species Distribution Modelling to Identify ‘Coldspots’ for Conservation of Freshwater Fishes under a Changing Climate. Aquat. Conserv. 2022, 32, 576–590. [Google Scholar] [CrossRef]
  82. Frederico, R.G.; Dias, M.S.; Jézéquel, C.; Tedesco, P.A.; Hugueny, B.; Zuanon, J.; Torrente-Vilara, G.; Ortega, H.; Hidalgo, M.; Martens, K.; et al. The Representativeness of Protected Areas for Amazonian Fish Diversity under Climate Change. Aquat. Conserv: Mar. Freshw. Ecosyst. 2021, 31, 1158–1166. [Google Scholar] [CrossRef]
  83. Hrbek, T.; Stölting, K.N.; Bardakci, F.; Küçük, F.; Wildekamp, R.H.; Meyer, A. Plate Tectonics and Biogeographical Patterns of the Pseudophoxinus (Pisces: Cypriniformes) Species Complex of Central Anatolia, Turkey. Mol. Phylogenet. Evol. 2004, 32, 297–308. [Google Scholar] [CrossRef]
  84. Morueta-Holme, N.; Fløjgaard, C.; Svenning, J.-C. Climate Change Risks and Conservation Implications for a Threatened Small-Range Mammal Species. PLoS ONE 2010, 5, e10360. [Google Scholar] [CrossRef]
  85. Erk’akan, F.; Atalay-Ekmekçi, F.G.; Nalbant, T.T. A Review of the Genus Cobitis in Turkey (Pisces: Ostariophysi: Cobitidae). Hydrobiologia 1999, 403, 13–26. [Google Scholar] [CrossRef]
  86. Küçük, F. Extinct Endemic Fishes of Turkey: Alburnus akili (Gövce) and Pseudophoxinus handlirschi (Kavinne) (Pisces: Cyprinidae). Turk. J. Fish. Aquat. Sci. 2012, 12, 21. [Google Scholar] [CrossRef]
  87. Innal, D.; Erk’akan, F. Effects of Exotic and Translocated Fish Species in the Inland Waters of Turkey. Rev. Fish Biol. 2006, 16, 39–50. [Google Scholar] [CrossRef]
  88. Yoğurtçuoğlu, B.; Uyan, U.; Ekmekçi, F.G. The Influence of Environmental Instability on the Reproductive Strategy of the Critically Endangered Acıgöl Killifish (Aphanius transgrediens). J. Fish Biol. 2020, 97, 246–256. [Google Scholar] [CrossRef] [PubMed]
  89. Yerli, S.V.; Alp, A.; Yeğen, V.; Uysal, R.; Apaydın Yağcı, M.; Balık, İ. Evaluation of the Ecological and Economical Results of the Introduced Alien Fish Species in Lake Eğirdir, Turkey. Turk. J. Fish. Aquat. Sci. 2013, 13, 795–809. [Google Scholar] [CrossRef] [PubMed]
  90. İlhan, A. Threatened Fishes of the World: Chondrostoma beysehirense Bogutskaya, 1997 (Cyprinidae). Environ. Biol. Fish. 2009, 86, 483–484. [Google Scholar] [CrossRef]
  91. Özuluğ, M.; Oǧuz Öztürk, M. Threatened Fishes of the World: Pseudophoxinus anatolicus (Hanko 1924) (Cyprinidae), Central Anatolia, Turkey. Environ. Biol. Fish. 2008, 83, 183–184. [Google Scholar] [CrossRef]
  92. McGarvey, D.J.; Menon, M.; Woods, T.; Tassone, S.; Reese, J.; Vergamini, M.; Kellogg, E. On the Use of Climate Covariates in Aquatic Species Distribution Models: Are We at Risk of Throwing out the Baby with the Bath Water? Ecography 2018, 41, 695–712. [Google Scholar] [CrossRef]
  93. Pépino, M.; Rodríguez, M.A.; Magnan, P. Fish Dispersal in Fragmented Landscapes: A Modeling Framework for Quantifying the Permeability of Structural Barriers. Ecol. Appl. 2012, 22, 1435–1445. [Google Scholar] [CrossRef]
  94. Kates, D.; Dennis, C.; Noatch, M.R.; Suski, C.D. Responses of Native and Invasive Fishes to Carbon Dioxide: Potential for a Nonphysical Barrier to Fish Dispersal. Can. J. Fish. Aquat. Sci. 2012, 69, 1748–1759. [Google Scholar] [CrossRef]
  95. Maynard, D.S.; Crowther, T.W.; Bradford, M.A. Competitive Network Determines the Direction of the Diversity–Function Relationship. Proc. Natl. Acad. Sci. USA 2017, 114, 11464–11469. [Google Scholar] [CrossRef] [PubMed]
  96. Lake, P.S.; Bond, N.; Reich, P. Linking Ecological Theory with Stream Restoration. Freshw. Biol. 2007, 52, 597–615. [Google Scholar] [CrossRef]
  97. Schmidt, H.; Radinger, J.; Teschlade, D.; Stoll, S. The Role of Spatial Units in Modelling Freshwater Fish Distributions: Comparing a Subcatchment and River Network Approach Using MaxEnt. Ecol. Model. 2020, 418, 108937. [Google Scholar] [CrossRef]
  98. Hijmans, R.J.; Graham, C.H. The Ability of Climate Envelope Models to Predict the Effect of Climate Change on Species Distributions. Glob. Chang. Biol. 2006, 12, 2272–2281. [Google Scholar] [CrossRef]
  99. Logez, M.; Bady, P.; Pont, D. Modelling the Habitat Requirement of Riverine Fish Species at the European Scale: Sensitivity to Temperature and Precipitation and Associated Uncertainty. Ecol. Freshw. Fish 2012, 21, 266–282. [Google Scholar] [CrossRef]
  100. Conti, L.; Comte, L.; Hugueny, B.; Grenouillet, G. Drivers of Freshwater Fish Colonisations and Extirpations under Climate Change. Ecography 2015, 38, 510–519. [Google Scholar] [CrossRef]
  101. Scheffers, B.R.; Edwards, D.P.; Diesmos, A.; Williams, S.E.; Evans, T.A. Microhabitats Reduce Animal’s Exposure to Climate Extremes. Glob. Chang. Biol. 2014, 20, 495–503. [Google Scholar] [CrossRef]
  102. Suggitt, A.J.; Gillingham, P.K.; Hill, J.K.; Huntley, B.; Kunin, W.E.; Roy, D.B.; Thomas, C.D. Habitat Microclimates Drive Fine-Scale Variation in Extreme Temperatures. Oikos 2011, 120, 1–8. [Google Scholar] [CrossRef]
  103. MoEUCC. Climate Adaptation Strategy and Action Plan of Türkiye; Ministery of Environment, Urbanisation and Climate Change: Ankara, Turkey, 2011.
  104. Ergönül, M.B.; Breine, J.; Atasağun, S. Length-Weight Relationships of Two Threatened Species Endemic to Turkey: Gobio insuyanus Ladiges and Gobio microlepidotus Battalgil. Fish. Aquat. Life 2019, 27, 118–121. [Google Scholar] [CrossRef]
  105. Freyhof, J.; Özuluğ, M.; Saç, G. Neotype Designation of Aphanius iconii, First Reviser Action to Stabilise the Usage of A. fontinalis and A. meridionalis and Comments on the Family Group Names of Fishes Placed in Cyprinodontidae (Teleostei: Cyprinodontiformes). Zootaxa 2017, 4294, 573–585. [Google Scholar] [CrossRef]
  106. Yegen, V.; Sari, H.M. Endemik Eğirdir Yağ Balığı (Pseudophoxinus egridiri) Populasyonunun Büyüme Özellikleri. J. Limnol. Fish. Res. 2021, 7, 241–249. [Google Scholar] [CrossRef]
  107. Seçer, B.; Sungur, S.; Çiçek, E.; Ceylan, M. Freshwater Fish Fauna of Niğde. J. Limnol. Fish. Res. 2020. [Google Scholar] [CrossRef]
  108. Çiçek, E. Seminemacheilus dursunavsari, a New Nemachelid Species (Teleostei: Nemacheilidae) from Turkey. Iran. J. Ichthyol. 2020, 7, 68–77. [Google Scholar]
  109. Yoğurtçuoğlu, B.; Kaya, C.; Geiger, M.F.; Freyhof, J. Revision of the Genus Seminemacheilus, with the Description of Three New Species (Teleostei: Nemacheilidae). Zootaxa 2020, 4802, 477–501. [Google Scholar] [CrossRef]
  110. Ekmekçi, F.G. The Effects of High Salinity on the Production of Capoeta tinca in a Naturally Contaminated River. Turk. J. Zool. 2002, 26, 265–270. [Google Scholar]
  111. İlhan, A.; Balık, S.; Sarı, H.M. Orta ve Batı Anadolu Endemik Içsu Balıklarının Günümüzdeki Dağılımları ve Koruma Statüleri. J. Fish. Aquat. Sci. 2014, 29, 9–34. [Google Scholar] [CrossRef]
Figure 1. Distribution records of the studied species and their locations in the Central Anatolia Ecoregion in Türkiye.
Figure 1. Distribution records of the studied species and their locations in the Central Anatolia Ecoregion in Türkiye.
Water 15 01619 g001
Figure 2. Current and future distribution of Anatolichthys anatoliae, Cobitis bilseli, Egirdira nigra, and Gobio gymnostethus in the CAE. Colours show habitat suitability, and the value of suitability increases in red areas.
Figure 2. Current and future distribution of Anatolichthys anatoliae, Cobitis bilseli, Egirdira nigra, and Gobio gymnostethus in the CAE. Colours show habitat suitability, and the value of suitability increases in red areas.
Water 15 01619 g002
Figure 3. Current and future distribution of Gobio hettitorum, Gobio intermedius, Pseudophoxinus burduricus, and Squalius anatolicus in the CAE. Colours show habitat suitability, and the value of suitability increases in red areas.
Figure 3. Current and future distribution of Gobio hettitorum, Gobio intermedius, Pseudophoxinus burduricus, and Squalius anatolicus in the CAE. Colours show habitat suitability, and the value of suitability increases in red areas.
Water 15 01619 g003
Figure 4. Changes in mean annual temperature (a) and annual total precipitation (b) for the period 1970–2021. Points: data points; solid lines: predicted values (from the generalized linear models—see Section 2.3); ribbons: 95% confidence intervals of the predicted values (from the generalized linear models).
Figure 4. Changes in mean annual temperature (a) and annual total precipitation (b) for the period 1970–2021. Points: data points; solid lines: predicted values (from the generalized linear models—see Section 2.3); ribbons: 95% confidence intervals of the predicted values (from the generalized linear models).
Water 15 01619 g004
Figure 5. Runoff trend in the Central Anatolia Ecoregion. Black: no trend, white: decreasing trend. (a) Runoff trend in May; (b) Runoff trend in August.
Figure 5. Runoff trend in the Central Anatolia Ecoregion. Black: no trend, white: decreasing trend. (a) Runoff trend in May; (b) Runoff trend in August.
Water 15 01619 g005
Table 1. Area under the receiver operating characteristic curve (AUC) values for the 16 endemic freshwater fish distribution models.
Table 1. Area under the receiver operating characteristic curve (AUC) values for the 16 endemic freshwater fish distribution models.
SpeciesAUC Values
Anatolichthys anatoliae (Leidenfrost 1912)0.88
Anatolichthys fontinalis (Akşiray, 1948)0.98
Anatolichthys iconii (Akşiray, 1948)0.98
Capoeta pestai (Pietschmann, 1933)0.85
Chondrostoma beysehirense Bogutskaya, 19970.93
Cobitis bilseli Battalgil, 19420.96
Egirdira nigra (Kosswig & Geldiay, 1952)0.97
Gobio gymnostethus Ladiges, 19600.97
Gobio hettitorum Ladiges, 19600.94
Gobio intermedius Battalgil, 19440.96
Oxynoemacheilus nasreddini Yoğurtçuoğlu, Kaya & Freyhof, 20210.94
Pseudophoxinus anatolicus (Hankó, 1925)0.88
Pseudophoxinus burduricus Küçük, Gülle, Güçlü, Çiftçi & Erdoğan, 20130.95
Pseudophoxinus crassus (Ladiges, 1960)0.89
Squalius anatolicus (Bogutskaya, 1997)0.91
Squalius recurvirostris Özuluğ & Freyhof, 20110.97
Table 2. Conservation status [80] of the endemic freshwater fishes of the CAE and changes in their distribution range for various future climate scenarios and periods (grey shaded areas indicate the predicted losses in species distribution range).
Table 2. Conservation status [80] of the endemic freshwater fishes of the CAE and changes in their distribution range for various future climate scenarios and periods (grey shaded areas indicate the predicted losses in species distribution range).
SpeciesConservation Status *Range Size Change (%)
SSP2-4.5
2041–2060
SSP2-4.5
2081–2100
SSP5-8.5
2041–2060
SSP5-8.5
2081–2100
Anatolichthys anatoliaeNT214216226239
Anatolichthys fontinalisNE−100−100−100−100
Anatolichthys iconiiNE1386152114681172
Capoeta pestaiCR150175162189
Chondrostoma beysehirenseEN13812310747
Cobitis bilseliEN−21−48−47−100
Egirdira nigraEN7.8−76−70−100
Gobio gymnostethusCR−100−100−100−100
Gobio hettitorumCR−100−100−100−100
Gobio intermediusEN6.1−3017−70
Oxynoemacheilus nasreddiniNE106919775
Pseudophoxinus anatolicusEN155164162172
Pseudophoxinus burduricusEN−97−100−100−100
Pseudophoxinus crassusEN118132130146
Squalius anatolicusLC251162177−26
Squalius recurvirostrisVU4922608.5
Note: * CR, critically endangered; EN, endangered; VU, vulnerable; NT, near threatened; LC, least concern; NE, not evaluated based on the IUCN Red List.
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Korkmaz, M.; Mangıt, F.; Dumlupınar, İ.; Çolak, M.A.; Akpınar, M.B.; Koru, M.; Pacheco, J.P.; Ramírez-García, A.; Yılmaz, G.; Amorim, C.A.; et al. Effects of Climate Change on the Habitat Suitability and Distribution of Endemic Freshwater Fish Species in Semi-Arid Central Anatolian Ecoregion in Türkiye. Water 2023, 15, 1619. https://doi.org/10.3390/w15081619

AMA Style

Korkmaz M, Mangıt F, Dumlupınar İ, Çolak MA, Akpınar MB, Koru M, Pacheco JP, Ramírez-García A, Yılmaz G, Amorim CA, et al. Effects of Climate Change on the Habitat Suitability and Distribution of Endemic Freshwater Fish Species in Semi-Arid Central Anatolian Ecoregion in Türkiye. Water. 2023; 15(8):1619. https://doi.org/10.3390/w15081619

Chicago/Turabian Style

Korkmaz, Mustafa, Fatih Mangıt, İlayda Dumlupınar, Mehmet Arda Çolak, Mustafa Berkay Akpınar, Meltem Koru, Juan Pablo Pacheco, Arely Ramírez-García, Gültekin Yılmaz, Cihelio Alves Amorim, and et al. 2023. "Effects of Climate Change on the Habitat Suitability and Distribution of Endemic Freshwater Fish Species in Semi-Arid Central Anatolian Ecoregion in Türkiye" Water 15, no. 8: 1619. https://doi.org/10.3390/w15081619

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop