Next Article in Journal
Organic Materials Promote Rhododendron simsii Growth and Rhizosphere Soil Properties in a Lead–Zinc Mining Wasteland
Previous Article in Journal
Biological Control of Three Major Cucumber and Pepper Pests: Whiteflies, Thrips, and Spider Mites, in High Plastic Tunnels Using Two Local Phytoseiid Mites
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Plants
Volume 13
Issue 6
10.3390/plants13060890
plants-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Distribution Pattern of Suitable Areas and Corridor Identification of Endangered Ephedra Species in China

by
Huayong Zhang
1,2,*,
Jiangpeng Li
1,
Hengchao Zou
1,
Zhongyu Wang
1,
Xinyu Zhu
3,
Yihe Zhang
4 and
Zhao Liu
2
1
Research Center for Engineering Ecology and Nonlinear Science, North China Electric Power University, Beijing 102206, China
2
Theoretical Ecology and Engineering Ecology Research Group, School of Life Sciences, Shandong University, Qingdao 250100, China
3
Dalian Eco-Environmental Affairs Service Center, No. 58 Lianshan Street, Shahekou District, Dalian 116026, China
4
School of Engineering, RMIT University, P.O. Box 71, Bundoora, VIC 3083, Australia
*
Author to whom correspondence should be addressed.
Plants 2024, 13(6), 890; https://doi.org/10.3390/plants13060890
Submission received: 9 February 2024 / Revised: 12 March 2024 / Accepted: 18 March 2024 / Published: 20 March 2024
(This article belongs to the Section Plant Ecology)
Browse Figures
Versions Notes

Abstract

:
The suitable habitat of endangered Ephedra species has been severely threatened and affected by climate change and anthropogenic activities; however, their migration trends and restoration strategies are still relatively understudied. In this study, we utilized the MaxEnt model to simulate the suitable habitats of five endangered Ephedra species in China under current and future climate scenarios. Additionally, we identified significant ecological corridors by incorporating the minimum cumulative resistance (MCR) model. Under the current climate scenario, the suitable area of Ephedra equisetina Bunge, Ephedra intermedia Schrenk ex Mey, Ephedra sinica Stapf, and Ephedra monosperma Gmel ex Mey comprised 16% of the area in China, while Ephedra rhytidosperma Pachom comprised only 0.05%. The distribution patterns of these five Ephedra species were primarily influenced by altitude, salinity, temperature, and precipitation. Under future climate scenarios, the suitable areas of E. equisetina, E. intermedia, and E. sinica are projected to expand, while that of E. monosperma is expected to contract. Notably, E. rhytidosperma will lose its suitable area in the future. Our identified ecological corridors showed that the first-level corridors encompassed a wider geographical expanse, incorporating E. equisetina, E. intermedia, E. sinica, and E. monosperma, while that of E. rhytidosperma exhibited a shorter length and covered fewer geographical areas. Overall, our study provides novel insights into identifying priority protected areas and protection strategies targeting endangered Ephedra species.

1. Introduction

Environmental change profoundly impacts the structure and functioning of ecosystems across the biosphere, as well as the composition of biological communities and patterns of species distribution [1]. Climate change is pivotal in shaping plant species’ distribution patterns, which have been intricately linked to recurrent climatic fluctuations since the Fourth Great Ice Age [2,3]. Warming temperatures and altered rainfall regimes greatly affect species diversity and even increase the risk of plant extinction [4,5], posing significant challenges to the stability of China’s ecosystems and the conservation of endangered species [6,7]. Under the influence of global warming, plant distribution patterns are expected to undergo varying degrees of change, potentially leading to the extinction of 20–30% of species [8,9], particularly those that are endangered or have limited geographical ranges and are designated as national key protected plants [10]. Although higher temperatures may increase the genetic diversity of species and thus allow for the expansion of their ranges [11,12], species are losing their original habitats due to the overuse of resources by humans and changes in land use patterns [13]. Therefore, simulating the suitable distribution of species and identifying ecological corridors within the context of climate stress can serve as a theoretical foundation for formulating and implementing sustainable measures to conserve endangered species.
In recent years, the Species Distribution Model (SDM) has emerged as a pivotal approach used in hotspot studies for species distribution determination [14]. Among various species distribution models, the maximum entropy model (MaxEnt) is widely regarded as the most accurate and extensively employed model [15]. It offers several advantages, including its ability to yield reliable predictions with small sample sizes, flexible variable processing, and high precision in budget estimation [16]. With the advancement of MaxEnt, it has gained extensive application in safeguarding endangered flora, harnessing medicinal resources, and combating invasive species [17,18,19,20,21,22,23,24,25,26,27]. After predicting species-suitable areas, ecological corridors are good at connecting fragmented and isolated habitats [28,29]. The most commonly employed approach for studying migratory ecological corridors in species is integrating the minimum cumulative resistance (MCR) model with Geographic Information Systems (GISs) [30,31]. Notably, both the MaxEnt model and MCR model have been successfully utilized in conserving numerous endangered animal species [32,33], yielding significant outcomes. However, there remains an urgent need to investigate potential ecological corridors for plants.
Ephedra species are widely distributed in northern and northwestern China [13], with a total of 12 Ephedra species identified in the country [34]. Among them, Ephedra equisetina Bunge, Ephedra intermedia Schrenk ex Mey, Ephedra sinica Stapf, Ephedra monosperma Gmel ex Mey, and Ephedra rhytidosperma Pachom have been classified as endangered at varying levels of protection (see Section 4 for references). Notably, E. rhytidosperma is currently on the brink of extinction [35,36]. As a renowned sand fixation and traditional Chinese medicine species, it plays an irreplaceable role in maintaining ecological stability and treating diseases [37,38]. Therefore, it is imperative to comprehensively understand its distribution pattern of suitable areas and implement effective conservation measures for this genus. Despite extensive research on Ephedra species encompassing distribution surveys, physiological indicators, population structure, captive cultivation techniques, suitable area studies, and spatial distribution mapping [13,22,23,24,25,39,40,41,42,43,44], limited attention has been devoted to conservation design based on assessing potential habitat suitability. Henceforth, identifying appropriate suitable areas and dispersal ecological corridors for endangered Ephedra species under climate stress conditions are crucial to conserving and utilizing these valuable resources.
In this study, we employed the MaxEnt model to predict the suitable areas of five endangered Ephedra species (E. equisetina, E. intermedia, E. sinica, E. monosperma, and E. rhytidosperma). Utilizing these distribution areas as a basis, we utilized the MCR model to identify ecological corridors suitable for the dispersal of Ephedra species and optimize the habitat pattern. The specific purposes of this study are as follows: (1) determine the distribution of suitable areas for endangered Ephedra species under the current climatic conditions and the main environmental factors affecting their distribution; (2) predict and analyze the distribution of suitable areas and migration trends of endangered Ephedra species under changing climate scenarios in the future; and (3) calculate the minimum cumulative resistance of the ecological source area based on the MCR model to clarify the location of the ecological corridor. This study aims to provide a scientific foundation and practical guidance for the conservation of endangered Ephedra species and the preservation of biodiversity.

2. Results

2.1. Contemporary Distribution Patterns of Endangered Ephedra Species

The MaxEnt model achieved AUC values of 0.901, 0.890, 0.903, 0.890, and 0.998 for the five endangered species of Ephedra in China (Figure S1). Over 16% of the total suitable areas are accounted for by the four species (E. equisetina, E. intermedia, E. sinica, and E. monosperma), with a certain degree of overlap and concentration being exhibited in their distributions, primarily in the southwestern part of North China, the eastern part of Northwest China, and the northern part of Central China (Figure 1). In contrast, E. rhytidosperma has an exceedingly small total suitable area, comprising merely 0.05%, predominantly concentrated within Ningxia’s Helan Mountain region (Figure 1 and Table 1). The key environmental factors influencing the spatial distribution of the suitable areas of E. equisetina, E. intermedia, E. sinica, E. monosperma, and E. rhytidosperma were topsoil calcium carbonate (16%), elevation (19.8%), precipitation of the wettest month (19.1%), elevation (33.1%), and topsoil pH (38.6%), respectively (Table 2).

2.2. Future Distribution Patterns of Endangered Ephedra Species

Under the future (2050s and 2090s) climate scenarios (SSP126, SSP370, and SSP585), an expanding trend is shown in the suitable areas of E. equisetina, E. intermedia, and E. sinica, and a shrinking trend is shown in the areas of E. monosperma, while E. rhytidosperma completely loses its suitable areas (Figure 2 and Table S1). Among them, E. equisetina and E. intermedia have more suitable areas, mainly distributed in the areas north of the Qinling Mountains and northern Xinjiang. However, E. sinica mainly grows in North China and Northeast China. Thus, the highly suitable areas for E. monosperma are reduced to the Gansu and Qinghai provinces.
The suitable areas of E. equisetina expand by more than 11% under future climate scenarios, reaching 51.12% in the SSP585 scenario (2090s) and expanding to Inner Mongolia and the three eastern provinces (Figure 3 and Table S1). E. intermedia’s suitable distribution area extends to limited areas in North China and northwest and southwest regions. According to the projections under the SSP585 scenario (2050s) and SSP126 scenario (2090s), the expansion of suitable areas for E. intermedia will not exceed 7%. However, alternative scenarios show its expansion surpassing 13%, with the most substantial increase observed in the SSP585 scenario (2090s), reaching up to 46.96%. The distribution of E. sinica will primarily expand northward. In the SSP585 scenario (2090s), its total suitable areas will reach its maximum value, whereas in the SSP126 scenario (2050s), the proportion of highly suitable areas is the highest, approaching 25%. The suitable areas for E. monosperma are projected to experience a reduction of over 5% under future climate scenarios, with a significant decrease of 12.42% anticipated under the SSP585 scenario (2050s). This decline is primarily observed across Qinghai, Tibet, Sichuan, Inner Mongolia, and Shanxi provinces.

2.3. Corridor Identification of Endangered Ephedra Species

The study area of E. equisetina in this research was situated in the eastern and northwestern regions of China. A total of nine ecological source sites were identified, along with the selection of 16 ecological corridors, including seven first-level corridors (Figure 4 and Table S2). These corridors primarily serve to connect various ecological source points in Northwest China. Notably, the longest ecological corridor spans from Jingyuan County in Gansu Province to Wenshui County in Shanxi Province, covering a remarkable distance of 625.475 km. Conversely, the shortest corridor stretches from Yongdeng County to Jingyuan County within Gansu Province, measuring a total length of 129.533 km. The research area of E. intermedia was primarily situated in the eastern part of the northwest region, encompassing a total of seven selected ecological source areas. Ultimately, thirteen optimal solutions were derived as diffusion ecological corridors, comprising six first-level corridors with a combined length of 1123.455 km and seven secondary corridors spanning a total distance of 1300.094 km. For E. sinica, nine ecological source areas were established within the North China research area. Among them, seven first-level diffusion corridors exclusively cover short distances, primarily situated within Qinghai province; additionally, eight secondary diffusion corridors exist, with the longest one extending from Uxin Banner in Inner Mongolia to Ar Horqin Banner in Inner Mongolia.
The study area of E. monosperma was primarily located in the Qinghai, Gansu, and Sichuan provinces. A total of eight ecological source sites were selected, resulting in the identification of 14 optimal solution diffusion ecological corridors. The cumulative length of these corridors amounted to 1959.849 km, with eight being classified as first-level corridors predominantly situated in Qinghai province. The study area of E. rhytidosperma was mainly concentrated in the northern part of Ningxia and the eastern part of Gansu province. Eight ecological source sites were established, leading to the identification of seven first-level corridors and six secondary corridors. These corridors had respective lengths totaling 389.031 km and 659.515 km, with their minimum cost distances primarily located in the southern part of Helan Mountain.

3. Discussion

In this study, all species within the endangered Ephedra genus in China exhibited AUC values ranging from 0.85 to 1, indicating a high level of reliability in our predictions and excellent performance of the model’s simulation [25,26].
The findings demonstrate that E. equisetina, E. intermedia, E. sinica, and E. monosperma have extensive distribution areas suitable for their growth under current climate models in China with a total proportion exceeding 16%, making them widely distributed species primarily concentrated in Gansu, Yunnan, Shaanxi, Shanxi, Hebei Henan, and northwestern Inner Mongolia. In contrast, the suitable area for E. rhytidosperma is relatively small at only 0.05%, mainly concentrated in the Helan Mountain region of Ningxia as a narrowly distributed species. These results are consistent with records in “Flora of China” and the conclusions drawn by most scholars regarding Ephedra ‘s distribution [25,35,36,38,39,40,41,42,43]. The highly suitable areas for these five endangered Ephedra species are primarily concentrated in Central–Western China while being relatively scarce or absent in southern regions due to differences in the climatic conditions between north and south [45].
In the arid regions of Northwest China, natural environmental factors such as salinity, temperature, and precipitation play pivotal roles in shaping the potential geographic distribution and growth development of species [17,20,44,46]. The findings from this study revealed that elevation, salinity, and temperature significantly influenced the distribution patterns of E. equisetina, E. intermedia, and E. monosperma. Additionally, precipitation, salinity, and elevation were pivotal in determining the distribution ranges of E. sinica and E. rhytidosperma. Among them, the primary environmental factor influencing E. equisetina was topsoil calcium carbonate, while for E. rhytidosperma, it was topsoil pH (H2O), both of which are closely associated with the physiological characteristics of medicinal and desert plants [25,39,47,48]. It has been demonstrated that osmoregulation plays a pivotal role in enhancing the water retention capacity under arid conditions [49]. Osmoregulation facilitates the accumulation of both organic and inorganic solutes, thereby reducing osmotic pressure and promoting plant metabolism [40]. Alkaloids, as crucial organic osmoregulators, are indispensable for effective osmoregulation [50]. The primary environmental determinant for E. sinica is the precipitation of the wettest month. Previous research has demonstrated that “extreme” precipitation significantly influences the growth of E. sinica [13,42], which aligns with the findings of this study. For E. intermedia and E. monosperma, elevation is the dominant environmental factor due to its association with significant variations in water availability, temperature, and light intensity. These changes at different altitudes create a more diverse habitat conducive to the survival of rare and endangered plant species [51].
Under future climate scenarios, the suitable areas for three Ephedra species (E. equisetina, E. intermedia, and E. sinica) are projected to increase significantly. Particularly under the SSP585 scenario in the 2090s, there is a substantial expansion of total suitable areas by 51.12%, 46.96%, and 19.61%, respectively, for these three species. This phenomenon can be attributed to warmer temperatures resulting in higher minimum temperatures and reduced frost damage, which lowers shrub mortality rates while creating favorable conditions for their expansion and colonization [52,53,54,55]. Moderate temperature warming, along with increased precipitation, positively influence the dispersal and population growth of these species [55,56]. Studies indicate that surface temperatures in China will rise by 2.7–2.9 °C under the future global warming context, accompanied by an average annual increase in precipitation by 20%, particularly in the northern and northwestern regions [57,58]. Chinese scholars have also found that increased rainfall combined with rising temperatures will have a positive impact on the suitable habitats of these three Ephedra species, corroborating our study’s findings [25,39,40,41,42].
Climate change will result in the expansion of the suitable distribution areas for certain wide-ranging and invasive species while simultaneously reducing the size of suitable distribution areas for narrow-ranging and rare species [59]. Drastic climate change may even lead to the fragmentation of species distributions and habitat loss, thereby posing a greater threat to endangered species with limited natural ranges and exacerbating their risk of extinction [60]. Under projected future climate scenarios, E. monosperma is expected to undergo habitat reduction while E. rhytidosperma faces a complete loss of its habitat, highlighting the urgent need for conservation research. The decrease in suitable areas for E. monosperma under different climate scenarios can be attributed to the multifaceted impact of temperature, which exerts a more pronounced negative influence on vegetation than precipitation, particularly as temperatures continue to rise [60,61,62,63,64]. Previous studies have identified a negative correlation between rare and endangered plant species, such as Sinowilsonia henryi Hemsl, E. rhytidosperma, and Phoebe bournei Yen C Yang, with annual precipitation and average temperature. Moreover, future habitat predictions indicate that E. rhytidosperma is projected to experience a complete loss of its suitable distribution area, which corroborates the findings of this study [65]. Meanwhile, E. rhytidosperma exhibits a restricted distribution, confined solely to the flood fan area or shallow mountain foothills of Helan Mountain, characterized by a limited range and small population size in an unfavorable growth environment [35,36,65]. With recent land development and utilization activities in the eastern foothills of Helan Mountain, the habitat of E. rhytidosperma has suffered severe degradation, exacerbating its endangered status.
Previous research has demonstrated that, in response to future warming temperatures, most species are expected to undergo latitudinal shifts toward higher latitudes [19,20,66]. However, due to variations in their suitable capacities toward environmental changes, different species exhibit distinct migration patterns [67]. In the SSP370 scenario, suitable areas for E. equisetina, E. intermedia, E. sinica, and E. monosperma are projected to shift toward higher latitudes. Similarly, under the SSP585 scenario, suitable areas will also move poleward (with E. intermedia and E. monosperma shifting equatorward by the 2050s and 2090s, respectively). Conversely, only E. equisetina and E. sinica migrate to higher latitudes in the SSP126 scenario (2050s). These findings highlight that there are differences in direction and distance with respect to suitable areas among various species within the genus Ephedra under climate stress conditions. This confirms that the migratory abilities of species constitute one factor influencing their adaptability to future changes [66].
Ecological corridors play a crucial role in addressing wildlife and plant habitat fragmentation issues by breaking down biological isolation barriers, facilitating gene exchange among populations, mitigating biodiversity impacts caused by fragmented habitats, and further improving ecosystem services [13,64,68,69,70,71]. Our research findings indicate that habitats connected by ecological corridors are more effective in preserving local species than isolated patches [13,64,69,71]. Due to the dispersed growth characteristics of species in the Ephedra genus, establishing individual protected areas becomes challenging [72]. Therefore, the identification and establishment of ecological corridors can provide a more effective means of protecting Ephedra species. These ecological corridors not only consider the distance between points but also account for cumulative environmental resistance effects on Ephedra plant dispersal, enhancing their operational feasibility [73]. Moreover, this positive impact becomes stronger over time and plays a crucial role in the maintenance of local biodiversity, as well as the promotion of ecosystem balance and stability [74].
Notably, the selected five endangered Ephedra species exhibited largely unaffected ecological corridors by environmental factors such as topography and slope. At the same time, they are also situated at a considerable distance from residential areas [13]. Since Ephedra is a plant with a relatively weak dispersal ability, we prioritized short-range corridors in the selection process as important ecological corridors. However, this study predicted suitable areas for Ephedra under future climate stress conditions, and as a result, long-distance corridors were also chosen [13,72].
Based on the outcomes of assessing the adaptability of endangered Ephedra species to climate-change-induced stress, suitable areas for E. equisetina in Henan and Hebei will shrink in future climate scenarios but expand toward Inner Mongolia and the three eastern provinces. Consequently, over time, the service value of three ecological corridors extending from Etoke Banner in Inner Mongolia to Zichang City in Shaanxi Province, from Etoke Banner in Inner Mongolia to Wenshui County in Shanxi Province, and from Jingyuan County in Gansu Province to Wenshui County in Shanxi Province will progressively gain significance. Therefore, it is crucial to prioritize their construction and restoration [25]. In the future climate scenario, more suitable areas for E. intermedia will be concentrated in Gansu, Ningxia, and Qinghai but will also expand into Inner Mongolia and Shaanxi. Henceforth, it is imperative to enhance the restoration of ecological corridors in Inner Mongolia within the existing corridor design [24]. Under future climate change scenarios, the suitable areas and high-suitableness areas of E. sinica will shift northward, resulting in a decrease in suitable areas within Shanxi and Hebei provinces. Consequently, the value of ecological corridors in Shanxi and Hebei provinces will decrease in value under future climate conditions. Therefore, it is crucial to prioritize constructing ecological corridors within Inner Mongolia, specifically from Siziwang Banner to Ar Horqin Banner, from Wengniute Banner to Ar Horqin Banner, and from Shangyi County in Hebei Province to Wengniute Banner in Inner Mongolia [13,39]. Regarding E. monosperma under future climate conditions, its highly suitable areas are still concentrated at the junction of the Qinghai, Gansu, and Sichuan provinces. Corridor construction should be dominated by the first-level corridor and supplemented by the secondary corridor [24,64,68,75,76]. E. rhytidosperma has a small distribution area and low population size and is highly endangered. A single corridor restoration can no longer satisfy the protection and management needs of the Ephedra species. On this basis, as suggested for other threatened rare plant species with a very restricted range or very small populations [77], protection is fundamental. In our case, this can be carried out by taking measures such as establishing nature reserves for special protection and management, maintaining the integrity of the habitat, and avoiding further fragmentation of the habitat caused by human factors [35,36].

4. Materials and Methods

4.1. Data Sources and Data Processing

According to the “National Key Protected Wild Plants” (2021 edition) and the “China Rare and Endangered Plant Information System-Protection of Wild Plants in Provinces and Municipalities”, there are five rare and endangered species of Ephedra in China [78,79]: E. equisetina, E. intermedia, E. sinica, E. monosperma, and E. rhytidosperma. Distribution data for these species were obtained from the Chinese Virtual Herbarium (CVH: https://www.cvh.ac.cn (accessed on 6 April 2023)), the Global Biodiversity Information Facility (GBIF: https://www.gbif.org (accessed on 6 April 2023)), and the Species Diversity Data Platform (http://www.especies.cn (accessed on 6 April 2023)). To prevent overfitting caused by dense sample distribution, ENMTools (http://purl.oclc.org/enmtools (accessed on 7 April 2023)) were used to eliminate duplicate distribution points while retaining only one valid point within a 1 km grid [80,81].
The base map of China used in this study was acquired from the National Center for Basic Geographic Information (https://ngcc.cn/ngcc (accessed on 6 April 2023)). The contemporary climate data were obtained from WorldClim (https://worldclim.org (accessed on 6 April 2023)) with a spatial resolution of 1 km, encompassing 19 bioclimatic variables. The elevation factors were obtained from WorldClim (https://worldclim.org (accessed on 6 April 2023)) with a spatial resolution of 1 km, and the slope and aspect data were extracted by using the surface analysis tool of ArcGIS 10.8.2. The soil factor data were obtained from the World Soil Database (http://www.iiasa.ac.at/web/home/research/researchPrograms/water/HWSD.html (accessed on 6 April 2023)) constructed by the Food and Agriculture Organization of the United Nations and the International Institute for Applied Systems Analysis, Vienna, with a spatial resolution of 1 km. The annual China Land Cover Dataset (CLCD) derived from Landsat and produced by Yang et al. is based on the GEE platform [82], with a spatial resolution of 30 m. NDVI data were obtained from the Research Center for Resource and Environmental Sciences of the Chinese Academy of Sciences (https://www.resdc.cn (accessed on 6 April 2023)) with a resolution of 1 km and annual averages from 2010 to 2020.
The future climate data (19 bioclimatic variables) were obtained from WorldClim version 2.1 (https://worldclim.org (accessed on 6 April 2023)) with a spatial resolution of 1 km, including the base period (1970–2000) and future (2050s: 2041–2060 and 2090s: 2081–2100) climate data. Three future climate scenarios (SSP126, SSP370, and SSP585) representing low, medium, and high CO2 emission scenarios were selected. Due to the lack of data on future topography and soils under climate change, this study defaults to no change in the nation’s topography and soils over the projected time period [22].
The environmental variables utilized in the MaxEnt and MCR models were adjusted to a spatial resolution of 1 km by using ArcGIS 10.8.2. Three categories of environmental factors, namely climate, elevation, and soil (Table 3), were selected for inclusion in the MaxEnt (version 3.4.4) model. In order to mitigate the adverse impact of multicollinearity among similar environmental factors on model overfitting, which can compromise the accuracy and precision of the results, this study extracted 36 environmental factors and conducted a Pearson correlation analysis by using IBM SPSS 27 software. If a strong correlation (R2 > 0.8) was identified between similar factors, the factor with the lowest percentage contribution was eliminated [22,25,27]. The MCR model incorporates four environmental variables, namely elevation, slope, land use type, and NDVI, to determine the minimum resistance surface influencing the dispersal of the Ephedra species [13].

4.2. Suitable Habitat Evaluation

The screened environmental factors and the natural distribution point data of five rare and endangered Ephedra species were imported into the MaxEnt model, and a value of 75% of the distribution point data was randomly selected as a training set for building the model, with the remaining 25% used as a test set for verification [83]. The results were output after being repeated 10 times, the output format was Logistic, and the rest of the options were adopted in the model default settings. The accuracy of the model prediction results was examined by using the AUC, which is 0–1 [84]. It is generally believed that AUC < 0.6 indicates that the prediction results fail, 0.6 ≤ AUC ≤ 0.8 indicates that the prediction results are generally accurate, 0.8 < AUC ≤ 0.9 indicates that the prediction results are better, and AUC > 0.9 indicates that the prediction results are excellent [85,86].
Using tools for reclassification in ArcGIS, the maximum test sensitivity plus specificity threshold (MTSPS) was used as a boundary to divide the suitable habitat into four levels: unsuitable areas, p ≤ MTSPS; generally suitable areas, MTSPS < p ≤ 0.5; moderately suitable areas, 0.5 < p ≤ 0.75; and highly suitable areas, p > 0.75 [87]. The ArcGIS10.8.2 tool was used to simulate and analyze the changes in the suitable areas of five endangered Ephedra species in different periods. The expansion area, contraction area, and retention area can be obtained by overlapping the current and future grid maps to analyze the change in the suitable areas over time.

4.3. MCR Modeling for Corridor Identification

Utilizing the MCR model, this study identifies ecological corridors that are conducive to the dispersal of five endangered Ephedra species. The model elucidates the movement process of these species from their source to target locations by identifying the most efficient or least resource-intensive path required to overcome resistance. A lower level of resistance facilitates smoother ecological flow [68,70,76].
This study combines the data of highly suitable and moderately suitable areas for E. equisetina, E. intermedia, E. sinica, E. monosperma, and E. rhytidosperma generated by the MaxEnt model as the fundamental data source for running the MCR model [13,88]. The relatively independent regions in these base data are characterized and classified, and attributes are extracted to form individual ecological source sites. Subsequently, the rasterized evaluation results are converted into vector polygons and transformed into individual ecological source points by using feature-to-point tools. By utilizing the cost distance module in ArcGIS 10.8.2, the cost distances and cost backlinkages between each pair of ecological source points are calculated to simulate the minimum cumulative resistance paths formed among them as a basis for identifying ecological corridors [89]. Next, four environmental variables including elevation, slope, land use type, and NDVI were employed to determine the minimum resistance surface affecting Ephedra species dispersal. The resistance gradient of the Ephedra species was divided into six resistance classes (1–6), with the highest resistance class being 6 and the smallest resistance class being 1 [72]. A correlation matrix of the resistance values was established by analyzing the resistance values between factors by using the statistical analysis software yaahp (version 10.3) [75,90]. Finally, weights for each resistance factor were calculated [32].
The gravity model can be employed to quantitatively assess the magnitude of interaction between sources and targets, thereby identifying crucial corridors [91]. In this study, we employed the gravity model to calculate the interaction forces among various ecological source points and delineated the resistance value range by using a natural breakpoint method [92]. The significance of ecological corridors for five Ephedra species was classified into first-level and secondary corridors and general ecological corridors. General ecological corridors were not considered in this study due to their long distances and high resistance values [13,23,24].

5. Conclusions

In this study, we employed the MaxEnt model to predict the distribution patterns of five endangered Ephedra species (E. equisetina, E. intermedia, E. sinica, E. monosperma, and E. rhytidosperma) under current and future climate scenarios. Additionally, using the MCR model, we identified the minimum resistance corridors between ecological source sites while quantitatively analyzing the corridor classes by applying the gravity model. The following conclusions can be drawn:
(1)
Under the current climate scenario, four species, namely E. equisetina, E. intermedia, E. sinica, and E. monosperma, inhabit over 16% of their suitable area in China. Conversely, E. rhytidosperma occupies a comparatively smaller proportion of its suitable areas in China at only 0.05%. Altitude was the most critical factor limiting the growth of E. intermedia and E. monosperma, while salinity played a pivotal role in constraining the growth of E. equisetina and E. rhytidosperma. In contrast, the distribution of E. sinica exhibited a stronger dependence on precipitation factors.
(2)
Under the future climate scenario, the expansion of the suitable area of E. equisetina, E. intermedia, and E. sinica was maximized under the SSP585 scenario (2090s). The suitable area of E. monosperma shrank, with the greatest degree of shrinkage in the SSP585 scenario (2050s), by 12.42%. In contrast, E. rhytidosperma loses its suitable area under future climate scenarios.
(3)
By employing the MCR model and gravity model, we successfully identified 71 crucial ecological corridors for five Ephedra species, all of which are strategically located away from anthropogenic surfaces. These corridors play a pivotal role in ensuring the long-term survival of Ephedra species. Identifying these ecological corridors provides valuable insights into the conservation and management strategies for Ephedra species.
In conclusion, E. equisetina, E. intermedia, and E. sinica demonstrate promising potential for development under future climate scenarios, while E. monosperma and E. rhytidosperma face unfavorable prospects or even the risk of extinction. This study provides valuable scientific guidance for the conservation planning of these five Ephedra species by understanding their habitats and identifying potential ecological corridors to promote sustainable development and biodiversity conservation.

Supplementary Materials

The following supporting information can be downloaded at https://www.mdpi.com/article/10.3390/plants13060890/s1: Figure S1. AUC training values for E. equisetina (a), E. intermedia (b), E. sinica (c), E. monosperma (d), and E. rhytidosperma (f). Table S1. Spatial distribution of suitable areas changes of E. equisetina, E. intermedia, E. sinica, E. monosperma, and E. rhytidosperma under future climate scenarios. Table S2. Ecological source coordinates of endangered Ephedra species.

Author Contributions

Conceptualization, H.Z. (Huayong Zhang); methodology, H.Z. (Huayong Zhang); software, J.L.; validation, H.Z. (Huayong Zhang), H.Z. (Hengchao Zou), X.Z. and Y.Z.; formal analysis, J.L.; writing—original draft preparation, H.Z. (Huayong Zhang), J.L. and H.Z. (Hengchao Zou); writing—review and editing, H.Z. (Huayong Zhang), J.L., H.Z. (Hengchao Zou), Z.W., X.Z., Y.Z. and Z.L.; visualization, J.L.; supervision, H.Z. (Huayong Zhang) and H.Z. (Hengchao Zou); funding acquisition, H.Z. (Huayong Zhang). All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by National Water Pollution Control and Treatment Science and Technology Major Project (2017ZX07101) and the Discipline Construction Program of Huayong Zhang, Distinguished Professor of Shandong University, School of Life Sciences (61200082363001).

Data Availability Statement

All links to the input data are reported in the manuscript, and all the output data are available upon request to the authors.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Srivastava, V.; Griess, V.C.; Keena, M.A. Assessing the Potential Distribution of Asian Gypsy Moth in Canada: A Comparison of Two Methodological Approaches. Sci. Rep. 2020, 10, 22. [Google Scholar] [CrossRef] [PubMed]
  2. Box, E.O.; Fujiwara, K. A Comparative Look at Bioclimatic Zonation, Vegetation Types, Tree Taxa and Species Richness in Northeast Asia. Bot. Pacifica 2012, 1, 5–20. [Google Scholar] [CrossRef]
  3. Ge, X.J.; Hsu, T.W.; Hung, K.H.; Lin, C.J.; Huang, C.C.; Huang, C.C.; Chiang, Y.C.; Chiang, T.Y. Inferring Multiple Refugia and Phylogeographical Patterns in Pinus massoniana Based on Nucleotide Sequence Variation and DNA Fingerprinting. PLoS ONE 2012, 7, e43717. [Google Scholar] [CrossRef] [PubMed]
  4. Hoover, D.L.; Knapp, A.K.; Smith, M.D. Resistance and resilience of a grassland ecosystem to climate extremes. Ecology 2014, 95, 2646–2656. [Google Scholar] [CrossRef]
  5. Olsen, S.L.; Töpper, J.P.; Skarpaas, O.; Vandvik, V.; Klanderud, K. From facilitation to competition: Temperature-driven shift in dominant plant interactions affects population dynamics in seminatural grasslands. Glob. Chang. Biol. 2016, 22, 1915–1926. [Google Scholar] [CrossRef] [PubMed]
  6. Xue, T.T.; Gadagkar, S.R.; Albright, T.P.; Yang, X.D.; Li, J.; Xia, C.Y.; Wu, J.Y.; Yu, S.X. Prioritizing conservation of biodiversity in an alpine region: Distribution pattern and conservation status of seed plants in the Qinghai-Tibetan Plateau. Glob. Ecol. Conserv. 2021, 32, e01885. [Google Scholar] [CrossRef]
  7. Wang, A.L.; Li, W.W. Genetic diversity of Rheum tanguticum (Polygonaceae), an endangered species on Qinghai-Tibetan Plateau. Biochem. Syst. Ecol. 2016, 69, 132–137. [Google Scholar] [CrossRef]
  8. Thomas, C.D.; Cameron, A.; Green, R.E.; Bakkenes, M.; Beaumont, L.J.; Collingham, Y.C.; Erasmus, B.F.N.; de Siqueira, M.F.; Grainger, A.; Hannah, L.; et al. Extinction risk from climate change. Nature 2004, 427, 145–148. [Google Scholar] [CrossRef]
  9. Fang, J.Y.; Zhu, J.L.; Shi, Y. The responses of ecosystems to global warming. Chin. Sci. Bull. 2018, 63, 136–140. [Google Scholar] [CrossRef]
  10. Cai, Q.Q.; An, M.T.; Yu, J.H.; Tian, L. Suitable habitat prediction of Pinus kwangtungensis in China under climate change. J. Shaanxi Norm. Univ. 2023, 52, 90–102. [Google Scholar] [CrossRef]
  11. Bálint, M.; Domisch, S.; Engelhardt, C.H.M.; Haase, P.; Lehrian, S.; Sauer, J.; Theissinger, K.; Pauls, S.U.; Nowak, C. Cryptic biodiversity loss linked to global climate change. Nat. Clim. Chang. 2011, 1, 313–318. [Google Scholar] [CrossRef]
  12. He, Y.Z.; Huang, W.D.; Zhao, X.; Lv, P. Review on the impact of climate change on plant diversity. J. Desert Res. 2021, 41, 59–66. [Google Scholar] [CrossRef]
  13. Liu, M.X.; Li, L.; Wang, S.Y.; Xiao, S.R.; Mi, J.L. Forecasting the future suitable growth areas and constructing ecological corridors for the vulnerable species Ephedra sinica in China. J. Nat. Conserv. 2023, 73, 126401. [Google Scholar] [CrossRef]
  14. Tracy, J.L.; Kantola, T.; Baum, K.A.; Coulson, R.N. Modeling fall migration pathways and spatially identifying potential migratory hazards for the eastern monarch butterfly. Landsc. Ecol. 2019, 34, 443–458. [Google Scholar] [CrossRef]
  15. Tsiftsis, S.; Djordjevic, V.; Tsiripidis, I. Neottia cordata (Orchidaceae) at its southernmost distribution border in Europe: Threat status and effectiveness of Natura 2000 Network for its conservation. J. Nat. Conserv. 2019, 48, 27–35. [Google Scholar] [CrossRef]
  16. Xu, W.M.; Du, Q.L.; Yan, S.; Cao, Y.; Liu, X.; Guan, D.X.; Ma, L.Q. Geographical distribution of As-hyperaccumulator Pteris vittata in China: Environmental factors and climate changes. Sci. Total Environ. 2022, 803, 149864. [Google Scholar] [CrossRef] [PubMed]
  17. Su, Z.H.; Pan, B.R.; Zhuo, L.; Li, W.J.; Liu, H.L. Impact of future climate change on distribution pattern of Tamarix taklamakanensis and its conservation revelation. Arid. Zone Res. 2018, 35, 150–155. [Google Scholar] [CrossRef]
  18. Zhang, K.L.; Yao, L.J.; Meng, J.S.; Tao, J. Maxent modeling for predicting the potential geographical distribution of two peony species under climate change. Sci. Total Environ. 2018, 634, 1326–1334. [Google Scholar] [CrossRef]
  19. Duan, Y.Z.; Wang, H.T.; Wang, C.; Du, Z.Y. Potential distribution of endangered plant Helianthemum songaricum in China under climate change. J. Plant Resour. Environ. 2020, 2, 55–68. [Google Scholar] [CrossRef]
  20. Lu, K.; He, Y.M.; Mao, W.; Zy, D.U.; Wang, L.J.; Liu, G.M.; Feng, W.J.; Duan, Y.Z. Potential geographical distribution and changes of Artemisia ordosicain China under future climate change. J. Appl. Ecol. 2020, 31, 3758–3766. [Google Scholar] [CrossRef]
  21. Liu, Y.; Miao, C.; Wang, H.S. Influence of climate change on distribution of suitable areas of larix plantation in China. Acta Ecol. Sin. 2023, 43, 9686–9698. [Google Scholar] [CrossRef]
  22. Zhang, H.Y.; Zhao, B.J.; Huang, T.S.; Chen, H.; Yue, J.J.; Tian, Y.L. Responses of the Distribution Pattern of the Suitable Habitat of Juniperus tibetica Komarov to Climate Change on the Qinghai-Tibet Plateau. Forests 2023, 14, 434. [Google Scholar] [CrossRef]
  23. Ma, S.M.; Wei, B.; Li, X.C.; Ji, X.R. The impacts of climate change on the potential distribution of Haloxylon ammodendron. Chin. J. Ecol. 2017, 36, 1243. [Google Scholar] [CrossRef]
  24. Zhang, D.F.; Zhang, Q.; Guo, J.; Sun, C.Z.; Xie, C.X. Research on the global ecological suitability and characteristics of regions with Angelica sinensis based on the MaxEnt model. Acta Ecol. Sin. 2017, 37, 5111–5120. [Google Scholar] [CrossRef]
  25. Rong, W.W.; Huang, X.; Niu, P.X.; Liu, J.Q.; Yan, R.; Su, J.J.; Hu, S.C.; Chu, G.M. Potentially suitable areas for traditional Chinese medicinal material Ephedra equisetina based on MaxEnt model. Acta Ecol. Sin. 2023, 43, 8631–8646. [Google Scholar] [CrossRef]
  26. Zhang, H.; Sun, P.; Zou, H.; Ji, X.; Wang, Z.; Liu, Z. Adaptive Distribution and Vulnerability Assessment of Endangered Maple Species on the Tibetan Plateau. Forests 2024, 15, 491. [Google Scholar] [CrossRef]
  27. Li, M.; Zhang, Y.; Yang, Y.; Wang, T.; Wu, C.; Zhang, X. Prediction of Historical, Current, and Future Configuration of Tibetan Medicinal Herb Gymnadenia orchidis Based on the Optimized MaxEnt in the Qinghai–Tibet Plateau. Plants 2024, 13, 645. [Google Scholar] [CrossRef]
  28. Jantz, P.; Goetz, S.; Laporte, N. Carbon stock corridors to mitigate climate change and promote biodiversity in the tropics. Nat. Clim. Chang. 2014, 4, 138–142. [Google Scholar] [CrossRef]
  29. Fei, Y.Q.; Huang, A.M.; Luo, Y.; Zhao, L. Ecological restoration of abandoned mines within giant panda habitat based on Maxent model and circuit theory. J. Environ. Eng. Technol. 2024, 14, 1–11. [Google Scholar] [CrossRef]
  30. Ye, Y.P.; Wang, S.N. Estimating urban suitable ecological land based on the minimum cumulative resistance model: A case study in Nanjing, China. In Proceedings of the IOP Conference Series: Earth and Environmental Science, Macao, China, 16–19 July 2019. [Google Scholar]
  31. Li, Z.; Ding, Y.; Wang, Y.L.; Chen, J.; Wu, F.M. Construction of Ecological Security Pattern in Mountain Rocky Desertification Area Based on MCR Model: A Case Study of Nanchuan, Chongqing. J. Ecol. Rural Environ. 2020, 36, 1046–1054. [Google Scholar] [CrossRef]
  32. Su, X.K.; Liu, S.L.; Dong, S.K.; Zhang, Y.; Wu, X.Y.; Zhao, H.D.; Zhao, Z.Z.; Sha, W. Effects of potential mining activities on migration corridors of Chiru (Pantholops hodgsonii) in the Altun National Nature Reserve, China. J. Nat. Conserv. 2015, 28, 119–126. [Google Scholar] [CrossRef]
  33. Zhuge, H.J.; Lin, D.Q.; Li, X.W. Identification of ecological corridors for Tibetan antelope and assessment of their human disturbances in the alpine desert of Qinghai-Tibet Plateau. Chin. J. Appl. Ecol. 2015, 26, 2504–2510. [Google Scholar] [CrossRef]
  34. Chen, S.L.; Gilbert, M.G. Flora of China Editorial Committee (2006) Flora of China; Science Press and Missouri Botanical Garden Press: Beijing, China, 2006; Volume 22. [Google Scholar]
  35. Shi, B.; Dou, J.D.; Huang, W.; Li, X.W. Community characteristics of Ephedra rhytidosperma in Helan Mountain of Ningxia, China. Chin. J. Plant Ecol. 2022, 46, 362–367. [Google Scholar] [CrossRef]
  36. Zhao, L.; Yang, B.; Liu, W.D.; Ma, L.; Gao, J.H. Sex ratio and spatial pattern of Ephedra rhytidosperma in Helan Mountain of Ningxia. Acta Ecol. Sin. 2022, 42, 10297–10304. [Google Scholar] [CrossRef]
  37. Xiu, L.M.; Miura, A.B.; Yamamoto, K.; Kobayashi, T.; Song, Q.H.; Kitamura, H.; Cyong, J.C. Pancreatic islet regeneration by ephedrine in mice with streptozotocin-induced diabetes. Am. J. Chin. Med. 2001, 29, 493–500. [Google Scholar] [CrossRef] [PubMed]
  38. Tao, C.; Du, C.G. Medicinal herbs for the management of desertification. J. Agric. Biotechnol. Ecol. 2010, 3, 124–132. [Google Scholar]
  39. Ma, X.H.; Lu, Y.Y.; Huang, D.D.; Zhu, T.T.; Lv, P.L.; Jin, L. Ecology suitability study of Ephedra intermedia. China J. Chin. Mater. Medica 2017, 42, 2068–2071. [Google Scholar] [CrossRef]
  40. He, P.; Li, J.Y.; Li, Y.F.; Xu, N.; Gao, Y.; Guo, L.F.; Huo, T.T.; Peng, C.; Meng, F.Y. Habitat protection and planning for three Ephedra using the MaxEnt and Marxan models. Ecol. Indic. 2021, 133, 108399. [Google Scholar] [CrossRef]
  41. Nzei, J.M.; Ngarega, B.K.; Mwanzia, V.M.; Musili, P.M.; Wang, Q.F.; Chen, J.M. The past, current, and future distribution modeling of four water lilies (Nymphaea) in Africa indicates varying suitable habitats and distribution in climate change. Aquat. Bot. 2021, 173, 103416. [Google Scholar] [CrossRef]
  42. Zhang, K.; Liu, Z.Y.; Abdukeyum, N.; Ling, Y.B. Potential Geographical Distribution of Medicinal Plant Ephedra sinica Stapf under Climate Change. Forests 2022, 13, 2149. [Google Scholar] [CrossRef]
  43. Guo, L.; Gao, Y.; He, P.; He, Y.; Meng, F. Modeling for Predicting the Potential Geographical Distribution of Three Ephedra Herbs in China. Plants 2023, 12, 787. [Google Scholar] [CrossRef]
  44. Cianfaglione, K.; Crișan, F.; Gafta, D. Soil Quality Enhances Seed Germination Success in Ephedra major—A Pilot Experiment. Plants 2023, 12, 438. [Google Scholar] [CrossRef] [PubMed]
  45. Remya, K.; Ramachandran, A.; Jayakumar, S. Predicting the current and future suitable habitat distribution of Myristica dactyloides Gaertn. using MaxEnt model in the Eastern Ghats, India. Ecol. Eng. 2015, 82, 184–188. [Google Scholar] [CrossRef]
  46. Sun, J.; Wang, N.A.; Niu, Z. Effect of Soil Environment on Species Diversity of Desert Plant Communities. Plants 2023, 12, 3465. [Google Scholar] [CrossRef] [PubMed]
  47. Yu, J.; Cao, G.C.; Rong, Z.L.; Li, H.F. Prediction of potential distribution of Ophiocordyceps sinensis in Chain based on Maxent model. Ecol. Sci. 2023, 42, 202–212. [Google Scholar] [CrossRef]
  48. Chen, M.Y.; Zhang, L.L.; Cao, B.; Bai, C.K. Predictive suitable distributions and principal environment variables of six medicinal plants of Scutellaria based on MaxEnt and GIS. J. Tradit. Chin. Vet. Med. 2022, 41, 7–11. [Google Scholar] [CrossRef]
  49. Song, H.; Feng, G.; Tian, C.Y.; Zhang, F.S. Osmotic adjustment traits of Suaeda physophora, Haloxylon ammodendron and Haloxylon persicum in field or controlled conditions. Plant Sci. 2006, 170, 113–119. [Google Scholar] [CrossRef]
  50. Jiang, L.; Wu, D.; Li, W.; Liu, Y.; Li, E.; Li, X.; Yang, G.; He, X. Variations in Physiological and Biochemical Characteristics of Kalidium foliatum Leaves and Roots in Two Saline Habitats in Desert Region. Forests 2024, 15, 148. [Google Scholar] [CrossRef]
  51. Cantón, Y.; Del Barrio, G.; Solé-Benet, A.; Lázaro, R. Topographic controls on the spatial distribution of ground cover in the Tabernas badlands of SE Spain. Catena 2004, 55, 341–365. [Google Scholar] [CrossRef]
  52. Polley, H.W.; Tischler, C.R.; Jobnson, H.B. Elevated atmospheric CO2 magnifies intra-specific variation in seedling growth of honey mesquite: An assessment of relative growth rates. Rangel. Ecol. Manag. 2006, 59, 128–134. [Google Scholar] [CrossRef]
  53. Wang, Y.X.; Chen, X.J.; Lou, S.N.; Hu, A.; Ren, J.F. Woody-plant encroachment in grasslands: A review of mechanisms and aftereffects. Acta Pratacult. Sin. 2018, 27, 219–227. [Google Scholar] [CrossRef]
  54. Kang, H.M.; Li, H.H.; Xu, D.H.; Qin, P.; Liu, M.L. Effects of elevated CO2 and temperature on Caragana jubata (Alpineshrub) growth and antioxidant systems. Acta Ecol. Sin. 2020, 40, 367–376. [Google Scholar] [CrossRef]
  55. Hong, C.C.; Wang, B.Z.; Bai, Y.H.; Zhu, Y.J.; Zhang, K.B.; Yang, X.H. Spatial and temporal variations of potential suitable areas of caragana roborovskyi in China under climate change scenarios. Acta Bot. Boreali-Occident. Sin. 2023, 43, 0856–0866. [Google Scholar] [CrossRef]
  56. Schöb, C.; Kammer, P.M.; Choler, P.; Veit, H. Small-scale plant species distribution in snowbeds and its sensitivity to climate change. Plant Ecol. 2009, 200, 91–104. [Google Scholar] [CrossRef]
  57. Wang, C.H.; Zhang, S.N.; Zhang, F.M.; Li, K.C.; Yang, K. On the increase of precipitation in the Northwestern China under the global warming. Adv. Earth Sci. 2021, 36, 980–989. [Google Scholar] [CrossRef]
  58. Li, B.F.; Chen, Y.N.; Chen, Z.S.; Xiong, H.G.; Lian, L.S. Why does precipitation in northwest China show a significant increasing trend from 1960 to 2010? Atmos. Res. 2016, 167, 275–284. [Google Scholar] [CrossRef]
  59. Fan, X.W.; Miao, C.Y.; Duan, Q.Y.; Shen, C.W.; Wu, Y. The Performance of CMIP6 Versus CMIP5 in Simulating Temperature Extremes Over the Global Land Surface. J. Geophys. Res.-Atmos. 2020, 125, e2020JD033031. [Google Scholar] [CrossRef]
  60. Root, T.L.; Price, J.T.; Hall, K.R.; Schneider, S.H.; Rosenzweig, C.; Pounds, J.A. Fingerprints of global warming on wild animals and plants. Nature 2003, 421, 57–60. [Google Scholar] [CrossRef]
  61. Xu, M.L.; Xue, X. Analysis on the effects of climate warming on growth and phenology of alpine plants. J. Arid. Land Resour. Environ. 2013, 27, 137–141. [Google Scholar] [CrossRef]
  62. Wang, D.; You, Q.L.; Jiang, Z.H.; Wu, W.B.; Jiao, Y. Analysis of extreme temperature changes in China based on the homogeneity-adjusteddata. Plateau Meteorol. 2016, 35, 1352–1363. [Google Scholar] [CrossRef]
  63. Fan, X.; Pan, J.W.; He, S.T. Prediction of the potential distribution of Rosa roxburghii under the background of climate change based on MaxEnt model. Acta Bot. Boreali-Occident. Sin. 2021, 41, 159–167. [Google Scholar] [CrossRef]
  64. Li, Q.W.; Li, S.; Cao, M.C.; Xu, H.G. Habitat Suitability Evaluation and Corridor Design of Muntiacus crinifrons in Qianjiangyuan National Park. J. Ecol. Rural. Environ. 2021, 37, 778–785. [Google Scholar] [CrossRef]
  65. Li, Q.L. Community Characteristics and Conservation Strategies of Endangered Plant adiantum nelumboides X. C. Zhang in Three Gorges Reservoir Region; Southwest University: Chongqing, China, 2018. [Google Scholar]
  66. Zhao, Y.F.; Pan, B.R. Potential geographical distributions of Tugarinovia in China under climate change scenarios. Arid. Zone Res. 2023, 40, 949–957. [Google Scholar] [CrossRef]
  67. Zhai, X.Y.; Shen, Y.F.; Zhu, S.H.; Tu, Z.H. Potential impacts of climate change in future on the geographical distributions of relic Liriodendron chinense. J. Trop. Subtrop. Bot. 2021, 29, 151–161. [Google Scholar] [CrossRef]
  68. Li, F.; Ye, Y.P.; Song, B.W.; Wang, R.S. Evaluation of urban suitable ecological land based on the minimum cumulative resistance model: A case study from Changzhou, China. Ecol. Model. 2015, 318, 194–203. [Google Scholar] [CrossRef]
  69. Li, R.; Jiao, Y.M.; Liu, X.; Liu, Z.L.; Gao, X. Suitability evaluation and corridor design of habitats for Green Peafowl based on MaxEnt Model. Chin. J. Ecol. 2019, 38, 919–926. [Google Scholar] [CrossRef]
  70. Akbaba, B.; Bulut, S. Inventory of Large Mammal Species in the Ilgaz Mountains (ankr): A Major Ecological Corridor in Anatolia. Hitte J. Sci. Eng. 2020, 7, 73–80. [Google Scholar] [CrossRef]
  71. Kong, L.Q.; Xu, W.H.; Xiao, Y.; Pimm, S.L.; Shi, H.; Ouyang, Z.Y. Spatial models of giant pandas under current and future conditions reveal extinction risks. Nat. Ecol. Evol. 2021, 5, 1309. [Google Scholar] [CrossRef]
  72. Yang, W.; Tong, Y.F.; Ma, T.; Wu, G.L.; Li, S.P. Study on characters of N, P and K of aeolian soil and fertilization for Ephedra sinica growing base. Pratacult. Sci. 2008, 25, 19–25. [Google Scholar]
  73. Shi, N.N.; Han, Y.; Wang, Q.; Quan, Z.J.; Luo, Z.L. Construction and optimization of ecological network for protected areas in Qinghai province. Chin. J. Ecol. 2018, 37, 1910–1916. [Google Scholar] [CrossRef]
  74. Rushdi, A.M.A.; Hassan, A.K. Reliability of migration between habitat patches with heterogeneous ecological corridors. Ecol. Model. 2015, 304, 1–10. [Google Scholar] [CrossRef]
  75. Li, L.; Shi, Z.H.; Yin, W.; Zhu, D.; Ng, S.L.; Cai, C.F.; Lei, A.L. A fuzzy analytic hierarchy process (FAHP) approach to eco-environmental vulnerability assessment for the danjiangkou reservoir area, China. Ecol. Model. 2009, 220, 3439–3447. [Google Scholar] [CrossRef]
  76. Tang, Y.H.; Gao, C.; Wu, X.F. Urban Ecological Corridor Network Construction: An Integration of the Least Cost Path Model and the InVEST Model. ISPRS Int. Geo-Inf. 2020, 9, 33. [Google Scholar] [CrossRef]
  77. Perrino, E.V.; Mahmoud, Z.N.A.; Valerio, F.; Tomaselli, V.; Wagensommer, R.P.; Trani, A. Synecology of Lagoecia cuminoides L. in Italy and evaluation of functional compounds presence in its water or hydroalcoholic extracts. Sci. Rep. 2023, 13, 20906. [Google Scholar] [CrossRef] [PubMed]
  78. Lu, Z.K.; Qin, H.N.; Jin, X.H.; Zhang, Z.X.; Yang, Q.W.; Hong, D.Y.; Li, D.Z.; Li, K.F.; Yuan, L.C.; Zhou, Z.H. On the necessity, principle, and process of updating the List of National Key Protected Wild Plants. Biodivers. Sci. 2021, 29, 1577–1582. [Google Scholar] [CrossRef]
  79. Wang, S.; Xie, Y. China Species Red List; Higher Education Press: Beijing, China, 2004. [Google Scholar]
  80. Warren, D.L.; Beaumont, L.J.; Dinnage, R.; Baumgartner, J.B. New methods for measuring ENM breadth and overlap in environmental space. Ecography 2019, 42, 444–446. [Google Scholar] [CrossRef]
  81. Warren, D.L.; Glor, R.E.; Turelli, M. ENMTools: A toolbox for comparative studies of environmental niche models. Ecography 2010, 33, 607–611. [Google Scholar] [CrossRef]
  82. Yang, H.F.; Zhong, X.N.; Deng, S.Q.; Xu, H. Assessment of the impact of LUCC on NPP and its influencing factors in the Yangtze River basin, China. Catena 2021, 206, 105542. [Google Scholar] [CrossRef]
  83. Moreno, R.; Zamora, R.; Molina, J.R.; Vasquez, A.; Herrera, M.A. Predictive modeling of microhabitats for endemic birds in South Chilean temperate forests using Maximum entropy (Maxent). Ecol. Inform. 2011, 6, 364–370. [Google Scholar] [CrossRef]
  84. Ward, D.F. Modelling the potential geographic distribution of invasive ant species in New Zealand. Biol. Invasions 2007, 9, 723–735. [Google Scholar] [CrossRef]
  85. Yang, Z.B.; Bai, Y.; Alatalo, J.M.; Huang, Z.D.; Yang, F.; Pu, X.Y.; Wang, R.B.; Yang, W.; Guo, X.Y. Spatio-temporal variation in potential habitats for rare and endangered plants and habitat conservation based on the maximum entropy model. Sci. Total Environ. 2021, 784, 147080. [Google Scholar] [CrossRef] [PubMed]
  86. Zhang, H.; Yuan, H.; Zou, H.; Zhu, X.; Zhang, Y.; Wang, Z.; Liu, Z. Global Warming Drives Expansion of Endangered Spruce Forest on the Tibetan Plateau. Sustainability 2024, 16, 2164. [Google Scholar] [CrossRef]
  87. Chen, Y.G.; Le, X.G.; Chen, Y.H.; Chen, W.X.; Du, J.G.; Zhong, Q.L.; Cheng, D.L. Identification of the potential distribution area of Cunninghamia lanceolata in China under climate change based on the MaxEnt model. Chin. J. Appl. Ecol. 2022, 33, 1207–1214. [Google Scholar] [CrossRef]
  88. Van der Windt, H.J.; Swart, J.A.A. Ecological corridors, connecting science and politics: The case of the Green River in the Netherlands. J. Appl. Ecol. 2008, 45, 124–132. [Google Scholar] [CrossRef]
  89. Teng, Y.; Zhang, Z.; Zhang, S.L.; Yang, Y.X.; He, W.; Wang, N.; Zhang, Z.Y.; Bao, W.D. Habital suitabilility analysis and ecological corridor designs for red deer (Cervus elaphus) in the southern Greater Khingan Mountains. Acta Ecol. Sin. 2022, 42, 5990–6000. [Google Scholar] [CrossRef]
  90. Bastille-Rousseau, G.; Wittemyer, G. Characterizing the Landscape of Movement to Identify Critical Wildlife Habitat and Corridors; John Wiley & Sons, Ltd.: Hoboken, NJ, USA, 2021. [Google Scholar] [CrossRef]
  91. Brandao, M.M.; Vieira, F.D.; Nazareno, A.G.; de Carvalho, D. Genetic diversity of neotropical tree Myrcia splendens (Myrtaceae) in a fragment-corridor system in the Atlantic rainforest. Flora 2015, 216, 35–41. [Google Scholar] [CrossRef]
  92. Yao, C.Y.; He, Y.M.; Cheng, J.X.; Zhang, T.Y.; Pan, H.Y.; Ma, H.J. Evaluation of ecological security pattern and optimization suggestions in MinjiangRiver Basin based on MCR model and gravity model. Acta Ecol. Sin. 2023, 43, 7083–7096. [Google Scholar] [CrossRef]
Plants 13 00890 g001
Figure 1. Spatial distribution of suitable areas of E. equisetina, E. intermedia, E. sinica, E. monosperma, and E. rhytidosperma under the current climate.
Figure 1. Spatial distribution of suitable areas of E. equisetina, E. intermedia, E. sinica, E. monosperma, and E. rhytidosperma under the current climate.
Plants 13 00890 g001
Plants 13 00890 g002
Figure 2. The spatial distribution of suitable areas of E. equisetina, E. intermedia, E.sinica, E. monosperma, and E. rhytidosperma under future climate scenarios.
Figure 2. The spatial distribution of suitable areas of E. equisetina, E. intermedia, E.sinica, E. monosperma, and E. rhytidosperma under future climate scenarios.
Plants 13 00890 g002
Plants 13 00890 g003
Figure 3. Distribution of suitable areas expansion and contraction of E. equisetina, E. intermedia, E. sinica, and E. monosperma under future climate scenarios in comparison to the current status.
Figure 3. Distribution of suitable areas expansion and contraction of E. equisetina, E. intermedia, E. sinica, and E. monosperma under future climate scenarios in comparison to the current status.
Plants 13 00890 g003
Plants 13 00890 g004
Figure 4. Ecological corridor identification of E. equisetina, E. intermedia, E. sinica, E. monosperma, and E. rhytidosperma.
Figure 4. Ecological corridor identification of E. equisetina, E. intermedia, E. sinica, E. monosperma, and E. rhytidosperma.
Plants 13 00890 g004
Table 1. Suitable areas of E. equisetina, E. intermedia, E. sinica, E. monosperma, and E. rhytidosperma under the current climate.
Table 1. Suitable areas of E. equisetina, E. intermedia, E. sinica, E. monosperma, and E. rhytidosperma under the current climate.
SpeciesGenerally Suitable Areas
/×104 km2		Moderately Suitable Areas
/×104 km2		Highly Suitable Areas
/×104 km2		Total Suitable Areas
/×104 km2Percentage/%
E. equisetina107.9047.2916.40171.6017.9
E. intermedia134.5562.2425.30222.1016.1
E. sinica89.7750.8414.11154.7216.1
E. monosperma124.4334.1630.51161.6416.8
E. rhytidosperma1.792.050.644.480.05
Table 2. Percent contribution of main environmental factors.
Table 2. Percent contribution of main environmental factors.
SpeciesDriving FactorContribution
E. equisetinaT_CACO3: topsoil calcium carbonate16
Elevation14.5
Bio6: min temperature of coldest month11.7
Bio2: mean diurnal range10.4
Bio12: annual precipitation10
E. intermediaElevation19.8
T_BS: topsoil base saturation14.2
Bio15: precipitation seasonality10.6
Bio9: mean temperature of driest quarter9.5
T_CACO3: topsoil calcium carbonate7.8
E. sinicaBio13: precipitation of wettest month19.1
T_CACO3: topsoil calcium carbonate15.8
Bio19: precipitation of coldest quarter13
Elevation10.9
Bio15: precipitation seasonality10.9
E. monospermaElevation33.1
Bio2: mean diurnal range12.4
Slope8.9
Bio19: precipitation of coldest quarter8.8
Bio12: annual precipitation6.7
E. rhytidospermaT_PH_H2O: topsoil pH (H2O)38.6
Bio12: annual precipitation15.5
Bio11: mean temperature of coldest quarter13.3
Aspect10.2
T_CASO4: topsoil gypsum7.9
Table 3. Environment variables driving MaxEnt.
Table 3. Environment variables driving MaxEnt.
SymbolEnvironmental FactorsUnitSymbolEnvironmental FactorsUnit
Bio1Annual mean temperature°CBio19Precipitation of coldest quartermm
Bio2Mean diurnal range°CElevElevationm
Bio3Isothermality\AspectAspect°
Bio4Temperature seasonality\SlopeSlope°
Bio5Max temperature of warmest month°CT_PH_H2OTopsoil pH (H2O)−log (H+)
Bio6Min temperature of coldest month°CT_GRAVELTopsoil gravel content%
Bio7Temperature annual range°CT_SILTTopsoil silt fraction%
Bio8Mean temperature of wettest quarter°CT_CLAYTopsoil clay fraction%
Bio9Mean temperature of driest quarter°CT_SANDTopsoil sand fraction%
Bio10Mean temperature of warmest quarter°CT_OCTopsoil organic carbon%
Bio11Mean temperature of coldest quarter°CT_CEC_CLAYTopsoil CEC (clay)cmol/kg
Bio12Annual precipitationmmT_CEC_SOILTopsoil CEC (soil)cmol/kg
Bio13Precipitation of wettest monthmmT_BSTopsoil base saturation%
Bio14Precipitation of driest monthmmT_TEBTopsoil TEBcmol/kg
Bio15Precipitation seasonality\T_CACO3Topsoil calcium carbonate%
Bio16Precipitation of wettest quartermmT_CASO4Topsoil gypsum%
Bio17Precipitation of driest quartermmT_ESPTopsoil sodicity (ESP)%
Bio18Precipitation of warmest quartermmT_ECETopsoil salinity (Elco)dS/m
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

Zhang, H.; Li, J.; Zou, H.; Wang, Z.; Zhu, X.; Zhang, Y.; Liu, Z. Distribution Pattern of Suitable Areas and Corridor Identification of Endangered Ephedra Species in China. Plants 2024, 13, 890. https://doi.org/10.3390/plants13060890

AMA Style

Zhang H, Li J, Zou H, Wang Z, Zhu X, Zhang Y, Liu Z. Distribution Pattern of Suitable Areas and Corridor Identification of Endangered Ephedra Species in China. Plants. 2024; 13(6):890. https://doi.org/10.3390/plants13060890

Chicago/Turabian Style

Zhang, Huayong, Jiangpeng Li, Hengchao Zou, Zhongyu Wang, Xinyu Zhu, Yihe Zhang, and Zhao Liu. 2024. "Distribution Pattern of Suitable Areas and Corridor Identification of Endangered Ephedra Species in China" Plants 13, no. 6: 890. https://doi.org/10.3390/plants13060890

APA Style

Zhang, H., Li, J., Zou, H., Wang, Z., Zhu, X., Zhang, Y., & Liu, Z. (2024). Distribution Pattern of Suitable Areas and Corridor Identification of Endangered Ephedra Species in China. Plants, 13(6), 890. https://doi.org/10.3390/plants13060890

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