
==== Front
Sci Rep
Sci Rep
Scientific Reports
2045-2322
Nature Publishing Group UK London

39237616
71782
10.1038/s41598-024-71782-9
Article
Climate and land use/land cover changes increasing habitat overlap among endangered crested ibis and sympatric egret/heron species
Tuohetahong Yilamujiang 1
Lu Ruyue 1
Guo Ruiyan 1
Gan Feng 1
Zhao Fuyue 1
Ding Sheng 1
Jin Saisai 1
Cui Huifang 1
Niu Kesheng 2
Wang Chao 2
Duan Wenbing 2
Ye Xinping yexinping@snnu.edu.cn

134
Yu Xiaoping yuxp64@163.com

134
1 https://ror.org/0170z8493 grid.412498.2 0000 0004 1759 8395 College of Life Sciences, Shaanxi Normal University, Xi’an, 710119 China
2 Shaanxi Hanzhong Crested Ibis National Nature Reserve, Hanzhong, 723300 China
3 https://ror.org/0170z8493 grid.412498.2 0000 0004 1759 8395 Research Center for UAV Remote Sensing, Shaanxi Normal University, Xi’an, 710119 China
4 https://ror.org/0170z8493 grid.412498.2 0000 0004 1759 8395 Changqing Teaching & Research Base of Ecology, Shaanxi Normal University, Xi’an, 710119 China
5 9 2024
5 9 2024
2024
14 2073628 3 2024
30 8 2024
© The Author(s) 2024
2024
https://creativecommons.org/licenses/by-nc-nd/4.0/ Open Access This article is licensed under a Creative Commons Attribution-NonCommercial-NoDerivatives 4.0 International License, which permits any non-commercial use, sharing, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if you modified the licensed material. You do not have permission under this licence to share adapted material derived from this article or parts of it. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by-nc-nd/4.0/.
Climate and land use/land cover (LULC) changes have far-reaching effects on various biological processes in wildlife, particularly interspecific interactions. Unfortunately, interspecific interactions are often overlooked when assessing the impacts of environmental changes on endangered species. In this study, we examined niche similarities and habitat overlaps between wild Crested Ibis and sympatric Egret and Heron species (EHs) in Shaanxi, China, using Ecological niche models (ENMs). We aimed to forecast potential alterations in habitat overlaps due to climate and LULC changes. The results showed that although EHs possess a broader niche breadth compared to the Crested Ibis, they still share certain niche similarities, as indicated by Schoener’s D and Hellinger’s I values exceeding 0.5, respectively. Notably, despite varying degrees of habitat reduction, the shared habitat area of all six species expands with the changes in climate and LULC. We suggest that with the climate and LULC changes, the habitats of sympatric EHs are likely to suffer varying degrees of destruction, forcing them to seek refuge and migrate to the remaining wild Ibis habitat. This is primarily due to the effective conservation efforts in the Crested Ibis habitat in Yangxian County and neighboring areas. Consequently, due to the niche similarity, they will share and compete for limited habitat resources, including food and space. Therefore, we recommend that conservation efforts extend beyond protecting the Crested Ibis habitat. It is crucial to control human activities that contribute to LULC changes to safeguard the habitats of both Crested Ibis and other sympatric birds.

Keywords

Crested Ibis
Egrets
Herons
Ecological niche model
Niche
Interspecific competition
Subject terms

Ecology
Zoology
Environmental sciences
issue-copyright-statement© Springer Nature Limited 2024
==== Body
pmcIntroduction

Climate change and its impact on biodiversity have emerged as a growing global concern1–3. The rising temperatures, shifting precipitation patterns, and alterations in habitat suitability directly influence the biological responses of species4, as well as their behavioral patterns5, circadian and circannual activity rhythms6, reproduction7, and both intraspecific and interspecific interactions8,9. These climate-driven shifts can also alter resource availability and potentially intensify competition among species that share the same ecological niche10,11. Moreover, anthropogenic land use and land cover (LULC) changes, combined with the effects of climate change, have significant implications for the diversity and structure of local ecological assemblages12,13, species interactions11, and local and regional distributions14. At the local scale, land fragmentation can also negatively impact the richness and composition of waterbird species15.

Furthermore, climate-LULC change can lead to shifts in the geographic distribution of species' ecological niches, creating opportunities for new interactions, including competition, among previously coexisting or non-coexisting similar species. For example, the interplay of climate change and anthropogenic activities has significantly impacted the distribution patterns of the Andean Condor (Vultur gryphus) and the Black Vulture (Coragyps atratus) throughout South America, leading to habitat overlapping that facilitate cohabitation16. Stenseth et al. investigated the competitive interactions between the Blue Tits (Cyanistes caeruleus) and the Great Tits (Parus major) in Western Europe and found that climate change can generate equilibrium conditions with local differences in species composition and spatial structures17. In the central Indian region, characterized by a mosaic of human-made and natural landscapes, three partially sympatric owl species exhibit a positive correlation with the presence of agriculture. The conversion of agricultural land may result in spatial segregation among these owl species18. Climate change can also impact the spatial overlaps between carnivores, intensifying competition due to water scarcity in arid and semi-arid regions and decreasing rainfall19. Studies have demonstrated that future climate change will significantly influence the patterns of spatial overlap between the Giant Pandas (Ailuropoda melanoleuca) and coexisting mammalian species in Qinling Mountains. Under moderate and worst-case climate change scenarios, competition between species will intensify11.

The Crested Ibis (Nipponia nippon), once thought to be nationally extirpated20, has experienced a remarkable recovery in its wild population over the last four decades. According to unpublished data from the Shaanxi Hanzhong Crested Ibis National Nature Reserve, the current wild population stands at approximately 6000 individuals. Despite the dedicated efforts to protect its habitat, the population still faces the inherent risk of extinction during the recovery phase21. Currently, the Hanjiang River Basin in southern Shaanxi is the sole remaining habitat for the Crested Ibis22. Among the coexisting species that partially overlap in diet with the Crested Ibis include the Little Egret (Egretta garzetta), Grey Heron (Ardea cinerea), Chinese Pond Heron (Aedeola bacchus), Cattle Egret (Bubulcus ibis), and Black-crowned Night Heron (Nycticorax nycticorax)23,24. Researchers have used multi-temporal niche models to predict potential shifts in the ecological niche and distribution of the Crested Ibis under future climate and land use scenarios25. Furthermore, studies have found that the Crested Ibis, as a tactile forager, can form mixed-species groups with the Little Egrets to reduce predation risk and can tolerate some intraspecific competition through resource depletion, aligning with the resource exploitative competition hypothesis23. However, the impact of future climate and land use changes on the habitat distribution of the current wild population, particularly in the context of competition with closely related sympatric species, remains unknown.

This study utilized ensemble Ecological Niche Models (ENMs) to simulate the overlapping distribution patterns of the Crested Ibis and Egret and Heron species (EHs), and analyzed their responses to climate and LULC changes. We propose the following hypotheses: Firstly, visually foraging EHs have a broader ecological niche breadth compared to the Crested Ibis, resulting in spatial overlap in their suitable habitats. Secondly, in comparison to climate change, alterations in LULC patterns within the densely populated habitats of the Crested Ibis will primarily influence its future distribution. Ultimately, our findings aim to provide new insights into niche similarities that contribute to the coexistence and competitive interactions among EHs and the Crested Ibis in a protected landscape (Fig. 1).Fig. 1 Study area and occurrence points of the Crested Ibis and the EHs. The map was created using ArcGIS Pro (https://www.esri.com/en-us/arcgis/products/arcgis-pro/).

Results

Predicted LULC changes

The evaluation of the model performance yielded a Kappa coefficient of 0.83 and an OA of 0.88, demonstrating the effectiveness of the model in capturing actual variations in LULC within the study area. These findings strongly support our decision to select the most appropriate parameter sets for predicting LULC changes up to the year 2050. Based on the simulation, assuming the transition rate remains consistent with that observed in the previous decade (2010–2020), significant LULC changes are anticipated in the study area by 2050, particularly along the Hanjiang River basin. Comparing the projected LULC for 2050 with the actual distribution in 2020 (Fig. 2), the following alterations are anticipated in different land use categories: (1) Cropland: a negligible decrease of − 0.88%; (2) Forest: a slight increase of 0.67%; (3) Grassland: a substantial decline of − 10.43%; (4) Shrubland: a moderate decrease of − 3.68%; (5) Wetland: a significant reduction of − 27.27%; (6) Water: a minor decrease of − 4.00%; (7) Urban: a significant increase of 65.11%. Analyzing the percentage of each land use type within the total area, in 2020, Cropland, Forest, Grassland, Shrubland, Wetland, Water, and Urban accounted for 17.4674%, 76.7046%, 4.1552%, 0.2323%, 0.0094%, 0.5457%, and 0.7222% of the total area, respectively. By 2050, these proportions are projected to shift to 17.31422%, 77.2211%, 3.7217%, 0.2237%, 0.0068%, 0.5239%, and 1.1925%, respectively.Fig. 2 LULC maps of the study area, illustrating the observed LULC in 2020 and the simulated LULC in 2050. The map was created using ArcGIS Pro (https://www.esri.com/en-us/arcgis/products/arcgis-pro/).

ENM evaluation and key environmental factors

Figure 3a showcases the average evaluation metrics obtained from 10 test runs, which were conducted to assess the predictive performance of ecological niche modeling. The average overall accuracy, calculated from six metrics, is 0.80 for the Crested Ibis, 0.60 for the Little Egret, 0.61 for the Grey Heron, 0.69 for the Chinese Pond Heron, 0.76 for the Eastern Cattle Egret, and 0.73 for the Black-crowned Night Heron. These results indicate that the selected predictors accurately simulated the distribution of these species, suggesting their suitability for future projections.Fig. 3 (a) Evaluation metrics of each individual model and (b) Importance of environmental variables for each species’ ecological niche model (CI, Crested Ibis; LE, Little Egret; GH, Grey Heron; CPH, Chinese Pond Heron; ECE, Eastern Cattle Egret; BNH, Black-crowned Night Heron).

Habitat selection for Crested Ibis and EHs is influenced by various factors, each with distinct and significant impacts. In general, elevation, slope, and human disturbance exert a more pronounced influence on all of these species. Specifically for the Crested Ibis, several factors are found to be crucial in shaping its environmental niche. These factors include precipitation of the coldest quarter (Bio19), elevation, nightlight intensity, slope, and LULC patterns (Fig. 3b).

Quantification of niche similarities

There was a substantial geographic overlap of niches among the species, as indicated by the equivalence values (Schoener's D and Hellinger's I), which generally exceeded 0.5 (Fig. 4). Among the species pairs, the Grey Heron and the Chinese Pond Heron exhibited the highest niche overlap values, with D = 0.71 and I = 0.92. Notably, the Little Egret (D = 0.62 and I = 0.93) and the Eastern Cattle Egret (D = 0.57 and I = 0.94), two species closely resembling the Crested Ibis (the focal species for conservation efforts), demonstrated a significant degree of niche similarity with the Crested Ibis. Figure 5 presents kernel density plots that illustrate the density distribution of each species in relation to the key environmental variables. The Crested Ibis exhibits a limited distribution range and a low occurrence density in terms of important factors such as Bio19, elevation, slope, and human disturbance. Moreover, it appears to be overshadowed by the EHs, which indicates a narrow niche breadth for the Crested Ibis (Table 1).Fig. 4 Ecological niche similarity (D: Schoener's D; I: Hellinger's I) between species (CI, Crested Ibis; LE, Little Egret; GH, Grey Heron; CPH, Chinese Pond Heron; ECE, Eastern Cattle Egret; BNH, Black-crowned Night Heron).

Fig. 5 Kernel density distribution of the Crested Ibis and EHs with respect to key environmental factors: (a) Elevation; (b) Slope; (c) Nightlight; (d) Precipitation of coldest quarter (Bio19).

Table 1 The maximum training sensitivity plus specificity (MTSS) threshold and niche breadths of the Crested Ibis and EHs.

Species	Crested Ibis	Little Egret	Grey Heron	Chinese Pond Heron	Eastern Cattle Egret	Black-crowned Night Heron	
MTSS threshold	0.29	0.31	0.58	0.44	0.56	0.51	
Niche breadth	0.34	0.52	0.57	0.64	0.48	0.57	

Changes in suitable habitats under current and future conditions

The predicted potentially suitable habitats for the Crested Ibis are primarily located on the southern slope of the Qinling Mountains. These suitable habitats, defined with a MTSS threshold of ≥ 0.29, covers a total area of 12,083.50 km2, mainly concentrated in the counties of Yangxian, Chenggu, Hantai, Mianxian, Xixiang, Shiquan, Foping, and Ningshan (Fig. 6a). The five species of Egrets and Herons exhibit varying geographical distributions (Figs. 6 and 7). Among them, the Little Egret has the widest distribution range, covering an area of 15,814.40 km2, while the Eastern Cattle Egret has the narrowest distribution range, occupying only 1,943.82 km2 (Fig. 8b). With the exception of the Eastern Cattle Egret, the other EHs partially overlap in their distribution areas with the Crested Ibis. However, the Eastern Cattle Egret exhibits nearly complete overlap in its distribution area with the Crested Ibis. Notably, the suitability maps reveal that areas with high suitability for EHs (indicated by warm-colored zones) coincide with regions where wild Crested Ibis are densely distributed (Fig. 6).Fig. 6 Habitat suitability of (a) Crested Ibis, (b) Little Egret, (c) Grey Heron, (d) Chinese Pond Heron, (e) Eastern Cattle Egret, and (f) Black-crowned Night Heron. The lower limit of the color band is the Maximum Training Sensitivity plus Specificity (MTSS) threshold. The map was created using ArcGIS Pro (https://www.esri.com/en-us/arcgis/products/arcgis-pro/).

Fig. 7 Current and future potential habitat distribution of (a) Crested Ibis, (b) Little Egret, (c) Grey Heron, (d) Chinese Pond Heron, (e) Eastern Cattle Egret, and (f) Black-crowned Night Heron. C, current; CC, climate change scenario; CLC, climate and LULC change combined scenario. The map was created using ArcGIS Pro (https://www.esri.com/en-us/arcgis/products/arcgis-pro/).

Fig. 8 Species richness distribution map (a), suitable habitat areas for individual species (b) and multi-species overlap (c) under current condition and different environmental change scenarios by 2050. The numbers displayed under the color bands and histogram represent the overlap of the corresponding numbers of species (CI, Crested Ibis; LE, Little Egret; GH, Grey Heron; CPH, Chinese Pond Heron; ECE, Eastern Cattle Egret; BNH, Black-crowned Night Heron). The map was created using ArcGIS Pro (https://www.esri.com/en-us/arcgis/products/arcgis-pro/).

Future projections indicate that the suitable habitat range of the Crested Ibis is expected to potentially expand by 16.22% under the scenario of climate change alone. However, when considering both climate and LULC changes together, the suitable habitat range is projected to contract by 27.07%. The EHs are anticipated to be more significantly affected by both climate change and LULC changes. Specifically, under the climate change scenario, the habitat areas of the Little Egret, Grey Heron, Chinese Pond Heron, Eastern Cattle Egret, and Black-crowned Night Heron are projected to decrease by 32.35%, 12.43%, 9.38%, 16.42%, and increase by 10.41%, respectively. In the combined scenario of climate change and LULC changes, the projected habitat decreases for the Little Egret, Grey Heron, Chinese Pond Heron, and Black-crowned Night Heron are 62.16%, 25.83%, 55.98%, and 30.88% respectively, while the habitat of the Eastern Cattle Egret is projected to increase by 149.20% (Figs. 7, 8b).

Changes in overlaps (richness) under current and future conditions

Under the current climate and LULC conditions, the region with the highest species richness, encompassing both the Crested Ibis and five EHs (richness = 6), is located in an area of 785.70 km2, predominantly within the Hanjiang River basin (Fig. 8a). Additionally, areas with higher species richness, containing the Crested Ibis and four EHs (richness = 5), as well as the Crested Ibis and three EHs (richness = 4), cover 966.33 km2 and 1,715.22 km2, respectively (Fig. 8a, c).

Under the scenario of climate change alone, the overlap of suitable areas between the Crested Ibis and three, four, and five EHs is projected to decrease by 12.76%, increase by 6.11%, and decrease by 10.54%, respectively. However, when considering both climate change and LULC changes together, the overlap of suitable areas between the Crested Ibis and three, four, and five EHs is expected to decrease by 20.55%, decrease by 1.83%, and increase by 40.16%, respectively (Fig. 8a, c).

Discussion

Niche Competition Theory believes that when different species inhabit similar niches, they compete for limited resources such as space and nutrients26–29. The competitive exclusion principle, a fundamental concept in ecology, suggests that complete competitors are unable to coexist30. This exclusion can occur through exploitative competition for limited resources or apparent competition facilitated by shared natural enemies31. This subtlety can be leveraged to explore the relevance of niche competition theory to our study subjects, offering insights into the interactions and potential coexistence of species with partially overlapping niches. Drawing upon niche theory, the coexistence of sympatric Crested Ibis and EHs can be attributed to the spatial and temporal distribution of their diet composition and the segregation of their dietary niches32. For instance, the Crested Ibis employs tactile foraging, whereas EHs rely on visual foraging. The Crested Ibis derives benefits from associating mixed-species flocks that include visual foragers like Little Egrets, enhancing antipredator advantages through collective detection and early warning mechanisms23. In contrast, tactile foragers generally have lower night vision capabilities compared to visual nocturnal feeders33. Furthermore, tactile foragers tend to exhibit more limited habitat use and rely on short-term environmental factors that concentrate prey during periods of reduced food availability34. These overlaps and distinctions in their foraging behaviors mitigate competitive pressures, facilitating their coexistence. However, increasing spatial overlap in their distribution after environmental alterations may lead to heightened competitive interactions, as they share the consumption of aquatic creatures within the same space35–38.

Our results indicated that precipitation in the coldest quarter (Bio19), altitude, nightlight intensity, slope and LULC are important factors in the habitat choice of the Crested Ibis. This finding aligns with previous studies that have highlighted the significance of precipitation and anthropogenic activities in influencing the habitat distribution of the Crested Ibis25. The kernel density distributions of the first four important factors for the Crested Ibis displayed variations compared to those observed in EHs in terms of the ranges and densities of these factors. However, the distribution range and density of all EHs on these factors exceeded that of the Crested Ibis, with the Ibis’s curve almost encompassed by them. Furthermore, our research revealed that the niche breadth of the Crested Ibis is narrower compared to that of the EHs. This finding supports the assertion by Boulangeat et al. that species with a restricted niche breadth are at a greater risk of rarity and extinction39. This is attributed to a combination of factors, including habitat loss and environmental change, which can disproportionately affect specialist species40.

The Crested Ibis currently has a limited population confined to the counties along the banks of the Hanjiang River, sharing its habitat with the EHs. Our modeling predictions indicate that the geographical distribution range of the Crested Ibis is expected to expand in the future due to climate change, which aligns with a previous study by Yiwen Sun et al. that also predicted an increase in habitat range and abundance22. However, when considering the combined effects of climate change and LULC changes, the suitable habitat range for the Crested Ibis actually shows a decrease. This suggests that the species remains at significant risk due to human activities, highlighting its vulnerability. Previous studies have identified excessive hunting and agrochemical overuse as the primary threats to the wild Crested Ibis41,42. These studies have also reported a continued decline in the effective population size of the species over time. Our findings align with these previous studies and emphasize the importance of controlling human activities that contribute to LULC changes, such as the extensive expansion of agricultural land and impervious areas. Although the Crested Ibis does benefit from human intervention to some extent35,41, it is crucial to regulate the intensity of human pressure and the rates of LULC changes to ensure the species' population recovery. It is worth mentioning that changes in climate and LULC will affect sympatric EHs differently. While climate change is expected to have a minor negative effect on the coexisting EHs, the impact of LULC changes is expected to be significant. Based on the variable importance test, their habitat distributions are predominantly influenced by human activities. In our model, we considered a moderated Shared Socioeconomic Pathways scenario, in which these species may experience different levels of habitat loss as human activities intensify.

Our findings provide novel evidence supporting the notion that the Crested Ibis and EHs exhibit overlapping but slightly distinct environmental niches, indicating a potential competition for resources, particularly food and space. The EHs, known for their visual foraging behavior, possess a broader ecological niche compared to the Crested Ibis, enabling them to inhabit a wide range of habitats and exploit a variety of food sources. This adaptability and flexibility in their survival strategies make them more resilient to environmental changes43,44. In contrast, the Crested Ibis, as a tactile forager, has a narrower ecological niche and specialized requirements. This specialization makes the Crested Ibis more reliant on specific resources, such as larger prey like loaches, eels, crabs, crayfish, and aquatic insects. This differs from EHs, which predominantly focus on small fish. Moreover, the Crested Ibis necessitates specific environmental conditions, such as winter-flooded rice paddies and riverbanks, in contrast to EHs, which can thrive in diverse wetland habitats23,45,46. Consequently, it is less competitive and more vulnerable to environmental changes, potentially facing increased competition from other species with similar but broader niche requirements22,25,35.

The findings of our study indicate that alterations in climate and LULC have a substantial influence on species richness. Climate change is anticipated to lead to a reduction in areas with high species richness, while areas with low species richness are expected to increase. In central Iran, mirroring our findings, the habitats shared by carnivores are projected to decline due to climate change19. This implies that climate change may potentially reduce competition for resources and space within certain species' habitats. However, the combined effects of projected LULC changes and climate change are predicted to decrease the extent of regions characterized by low species richness. Conversely, regions with high species richness, including the Crested Ibis, are expected to expand in response to these changes. Specifically, the areas where the five EHs and the Crested Ibis coexist, considered diversity hotspots, are projected to increase due to the combined impacts of climate and LULC changes. The decrease in suitable habitats for individual species and the increased overlaps with the habitat of Crested Ibis indicate that EHs are more inclined to inhabit the protected Crested Ibis' habitat, if their habitats are threatened in the future, intensifying competition between them.

Our LULC simulation results suggest that urban land areas are expected to expand, while wetland and water areas will likely decrease in the study area. In accordance with the present results, Yu et al. (2022) demonstrated that the loss and deterioration of wetlands and forests, as well as the reduction of foraging and nesting habitats in historical Crested Ibis distribution areas, were consequences of agricultural intensification, expansion of impervious surfaces, and the clearing of natural vegetations25. It is quite possible that the ranges of the EHs now overlap with the historical habitat range of the Crested Ibis47. As a result of human activities leading to the destruction of their own habitats, the EHs are compelled to seek refuge in the remaining wild Ibis habitats around Yangxian County, where the natural environment is more effectively protected. This can lead to an increase in species richness and competition intensity in these areas. It's worth noting that some EHs have a global distribution, unlike the wild Crested Ibis, which are exclusively found on the southern slope of the Qinling Mountains. We acknowledge that the EHs within our study area might be constrained by altitude from crossing the northern slopes of the Qinling, yet they could potentially migrate southward in the future, which may alleviate the competitive pressures arising from habitat overlap. However, this aspect falls beyond the scope of our current research, necessitating further studies focusing on population dispersal and migration dynamics. The findings of this study suggest that potential habitat overlaps, along with certain ecological niche similarities, indicate that a significant level of habitat overlap could pose challenges to the sustained growth of the Crested Ibis population.

To comprehensively understand interspecific competition, it is essential to carefully examine various elements, including the impact of both biotic and abiotic factors, the scale of analysis, diurnal activity patterns, food abundance, and population abundance and dynamics48,49. Unfortunately, the scarcity of comprehensive data prevented us from fully considering the impact of these factors in our study. Instead, we relied solely on quantifying the spatial distribution of habitats, which inevitably introduced some inherent bias. However, despite this limitation, our understanding of distributional overlap can still provide valuable insights for restoration efforts by considering the competitive interactions among target species. Conservationists can leverage this knowledge to implement measures that prevent the dominance of a particular species and promote the coexistence of multiple species. This may involve habitat restoration, careful selection of suitable sites for reintroduction, controlled harvesting or translocation of dominant species, or promoting the colonization of less competitive species.

Methods

Target species and presence records

In our study, we selected five EHs, along with the Crested Ibis, based on their popularity, coexistence, and the availability of occurrence data. These EHs include the Little Egret, Grey Heron, Chinese Pond Heron, Eastern Cattle Egret, and Black-crowned Night Heron. Presence records were gathered through a combination of extensive field surveys and the utilization of additional resources. Initially, we obtained the annual reports from the Shaanxi Hanzhong Crested Ibis National Nature Reserve, which formed the basis of our dataset. To enhance this dataset with information on the distribution of the Crested Ibis beyond the protected area, we conducted extensive field observations, documenting both the Crested Ibis and EHs in Yangxian County and the surrounding 18 counties (32°02′–34°23′N, 106°11′–109°53′E) in Hanzhong and Ankang City, Shaanxi Province, China, spanning from 2017 to 2020 (Fig. 1). By compiling a comprehensive dataset, we amassed a total of 238 occurrence points for wild Crested Ibis and 93 points for EHs. Furthermore, we enriched the EHs records by incorporating observations recorded in the China Bird Report database (http://www.birdreport.cn/). In total, we gathered 263 georeferenced occurrence records, including 68 records for the Little Egret, 36 for the Grey Heron, 42 for the Chinese Pond Heron, 62 for the Eastern Cattle Egret, and 55 for the Black-crowned Night Heron (Fig. 1). To address spatial autocorrelation, we applied a filtering process where we randomly selected one occurrence point within each 300 × 300 m grid, corresponding to the maximum cell size of the environmental variables.

Environmental variables

To comprehensively examine the habitat preferences, environmental niches, and distribution dynamics of the Crested Ibis and the EHs coexisting in the study area, our model incorporated a wide range of environmental factors. These factors included climate conditions, topographical attributes, human impacts, and the availability of water and food resources. To incorporate plausible variations in future climate, we utilized climate projections for the years 2041–2060 (referred to as 2050) from the BCC-CSM2-MR, MIROC6, and CMCC-ESM2 Global Circulation Models (GCMs) featured in the sixth assessment report of the Intergovernmental Panel on Climate Change (IPCC6). We considered the Shared Socioeconomic Pathways 2–45 (SSP2-45) scenario for the year 2050, averaging three GCMs. The SSP2-45 scenario is characterized by a “middle of the road” development pathway, with moderate economic growth, a growing global population, and a moderate level of environmental protection. In this scenario, there is a gradual transition towards a more sustainable environmental condition, but it is not as rapid or comprehensive as in some other scenarios50.

To obtain the bioclimatic variables, we sourced data from WorldClim, which has a spatial resolution of 30 arc-seconds (http://www.worldclim.org). The LULC variables were obtained from GlobeLand30 (http://www.globeland30.org/), while the topographical variables were sourced from Geospatial Data Cloud (http://www.gscloud.cn). Lastly, the human influence variables were obtained from the Resource and Environmental Science and Data Center (https://www.resdc.cn/). We calculated the distance to variables using the Euclidean Distance tool in the Spatial Analyst in ArcGIS 10.8.

To reduce multicollinearity between the variables, we performed a subset selection of the variable set using Pearson's correlation coefficient. We selected the variables with a coefficient of |r|< 0.75, reducing dependence among selected variables51. Ultimately, we retained the ten most independent and ecologically significant environmental factors to build the niche models. These factors included Aspect, Precipitation of the coldest quarter (Bio19), Isothermality (Bio3), Temperature annual range (Bio7), Elevation, Farmland, LULC, Nightlight, Slope, and Distance to water (Table 2). Table 2 Environmental variables for Ecological Niche Models (ENMs) and the Patch-generating Land Use Simulation (PLUS) model: sources and applications.

Data type	Environmental variables (Unit)	Modelling usage	Source	
Climatic data	Bio3, Isothermality (-)	ENMs	https://www.worldclim.org/	
Bio19, Precipitation of coldest quarter (℃)	ENMs	
Bio7, Temperature annual range (℃)	ENMs	
Bio1, Annual mean temperature (℃)	PLUS	
Bio12, Annual precipitation (mm)	PLUS	
Natural Environment data	Elevation (m)	ENMs, PLUS	Geospatial Data Cloud

http://www.gscloud.cn

	
Aspect (°)	ENMs	
Slope (°)	ENMs	
Soil type (-)	PLUS	Resource and Environmental Science and Data Center https://www.resdc.cn/	
Distance to water (m)	ENMs, PLUS	
LULC (2000, 2010, 2020) (-)	ENMs, PLUS	GLOBELAND30 (http://www.globeland30.org/)	
Human interference and socioeconomic data	Population density (persons/km2)	PLUS	Resource and Environmental Science and Data Center (2015) https://www.resdc.cn/	
GDP (ten thousand Yuan/km2)	PLUS	
Nightlight (-)	ENMs, PLUS	
Farmland (m)	ENMs	
Distance to primary road (m)	PLUS	National Geomatic Center of China (2015)

http://www.ngcc.cn/ngcc/

	
Distance to the secondary road (m)	PLUS	
Distance to highway (m)	ENMs, PLUS	
Distance to railway (m)	PLUS	
Distance to prefectural center (m)	PLUS	
Distance to county center (m)	PLUS	
Distance to residential area (m)	PLUS	

LULC simulation

To understand the dynamics of land in the study area and forecast its future trends, we conducted a LULC simulation using an advanced model known as the Patch-generating Land Use Simulation (PLUS) model52. The PLUS model extends the well-recognized cellular automata (CA) model to capture intricate spatio-temporal interactions and the underlying rules governing land use changes53,54. It incorporates a novel rule-mining strategy and a patch-generating mechanism, significantly enhancing its capability to represent the non-linear nature of LULC changes. This model is particularly suitable for simulating complex geographical evolution processes and has been widely employed in studies related to land use simulation55–57. To apply the PLUS model in this study, we utilized two historical LULC maps derived from GlobeLand30 as the fundamental input for future simulations. Recognizing the significant impact of climatic conditions on vegetation growth, soil moisture levels, and water availability crucial for biological processes; the influence of socio-economic factors on land use for residential, commercial, or industrial purposes; and the effects of natural environmental factors on agricultural productivity and land suitability, we identified environmental drivers that could potentially influence LULC changes. Following the research methodology outlined by Liang et al.52, we integrated a total of 15 driving factors, consisting of ten socioeconomic, two climatic, and three natural environmental factors, believed to play a role in the evolution of land type transitions (Table 2). The PLUS model consists of two main modules: the Land Expansion Analysis Strategy (LEAS) and the CA Model based on Multi-type Random Patch Seeds (CARS). We established specific parameters for the PLUS model as follows: (1) LEAS module: consisting of 20 regression trees, a sampling rate of 0.05, and mTry set to 15; (2) CARS module: with a neighborhood size of 3, patch generation rate of 0.7, expansion coefficient of 0.5, seed percentage of 0.0001, and neighborhood weights allocated as follows: 0.103 (Cropland), 0.495 (Forest), 0.151 (Grassland), 0.056 (Shrubland), 0.012 (Wetland), 0.039 (Water), and 0.376 (Urban). The neighborhood weights are established by calculating the ratio of the expansion areas specific to particular land use category in relation to the overall land expansion. Subsequently, these weights are fine-tuned based on expert insights and knowledge.

To evaluate the accuracy of the model, we initially simulated the LULC for the year 2020 based on the development probability derived from the land type transition that occurred between 2000 and 2010. Subsequently, we employed the validation function within the PLUS model to compare the simulated LULC data for 2020 with the actual contemporaneous data obtained from GlobeLand30-202058. For this comparison, we utilized the Confusion Matrix calculation method, with a selected sampling rate of 0.05. To evaluate the model’s performance, we computed two key metrics: the overall accuracy (OA)52 and the Kappa coefficient59. Higher values of OA and the Kappa coefficient indicate greater accuracy, with a score above 0.8 generally considered to indicate statistically satisfactory model performance60,61. Finally, we projected the LULC for the year 2050 to conduct future ecological niche analysis.

Ecological niche modeling

In our study, we utilized an ensemble modeling approach to enhance the calibration of ENMs for both the Crested Ibis and sympatric EHs. To create robust ENMs, we combined three widely used algorithms: Maximum Entropy (MaxEnt), Generalized Linear Model (GLM), and Random Forest (RDF). This ensemble approach allows us to capitalize on the strengths of each algorithm and improve the overall predictive performance.

To train the ENMs, we randomly divided the available data into two sets using the bootstrap random partition method. We allocated 80% of the data for model generation, and the remaining 20% was reserved for assessing the predictive accuracy of each model. The modeling process was implemented using the R package ENMTML62, and the significance of environmental variables was assessed using "imp_var" function. To allocate pseudo-absence points within the background area, we used the "GEO_ENV_KM_CONST" method. This method combines environmental and geographical approaches with a k-means non-agglomerative clustering procedure to evenly disperse points across environmental space. By doing so, it can mitigate the effects of data bias by providing a corrective measure, highlighting regions where the species is less likely to occur62–64.

To evaluate the performance of each algorithm, we considered six different evaluation metrics: True Skill Statistic (TSS)65, Boyce66, Kappa59, Sorensen67, Jaccard67, and Area Under the Curve (AUC)68. Higher values for these metrics indicate greater accuracy and reliability of the individual ENMs. To investigate the potential distribution of the selected species, we constructed a final model by combining the results of all the algorithms using a weighted mean approach69. This approach calculates suitability values by considering the performance of the algorithms. The weight of each algorithm was determined using the TSS values, as illustrated in the following equation:1 Wj=rj∑j=1hrj

where Wj represents the weight assigned to the model result j, rj represents the TSS value associated with the model result j, and h represents the number of model results.

To examine the potential impacts of future climate and LULC changes on the habitat distribution pattern of the species, we projected ENMs under two different scenarios: “climate change only” and “combination of climate and LULC changes”. The post-modelling analysis was performed in three steps using ArcGIS10.8. In the initial step, we categorized habitats into a binary map by applying the Maximum Training Sensitivity plus Specificity (MTSS) logistic threshold70. Subsequently, we performed an overlay analysis to compare the distribution of suitable habitats in the present timeframe and the projected future scenarios for the 2050s. Finally, we used SDM toolbox 2.071 to estimate species richness, which represents the total number of species in a grid cell. We also calculated the areas encompassing various classes of richness that include the habitat of the Crested Ibis.

Quantification of niche similarities

To assess and compare the environmental niches of the Crested Ibis and EHs, we employed ENMTools 1.4.372, a valuable tool for automating ENM creation, calculating similarity measures, and conducting statistical comparisons of habitat distributions. To quantify the niche overlap between the species, we calculated two similarity metrics: Schoener’s metric D and Hellinger’s index I, derived from the ecological environmental niche space73,74. Schoener’s D considers the probability distributions of species occurrence across different regions or cells, measuring the overlap of niches based on species abundance in those locations. Hellinger's I, on the other hand, relies solely on probability distributions without assuming species abundance patterns. Both metrics range from 0 to 1, where 0 indicates no overlap or complete niche divergence, and 1 indicates complete overlap or identical niches. To assess the niche breadth for each species, we calculated the Levins’ inverse concentration metric72,75. Niche breadth quantifies the range of habitats occupied by a species. A value of 0 indicates that the species is restricted to a single grid cell, while a value of 1 suggests an equally suitable environmental distribution across all grid cells. Therefore, species with a wider environmental distribution have higher niche breadth values. Finally, to explore the critical environmental variables associated with the niche similarity between the Crested Ibis and EHs, we employed the Kernel density estimation method. If the Kernel density plots for each factor show significant overlap, it implies that the environmental preferences of the species share similar values or exhibit similar patterns in their distribution. This indicates a potential correlation or dependency between the factors for each species. Conversely, if the plots do not overlap much and are distinct, it suggests that the factors are independent of each other76.

Acknowledgements

We would like to extend our heartfelt gratitude for the financial support received from the National Natural Science Foundation of China (32270541). Additionally, we would like to express our special appreciation to the Shaanxi Natural Reserve and Wildlife Management Station for its invaluable administrative and financial assistance.

Author contributions

X.Y., X.Y., Y.T. conceived and designed the study; R.L., R.G., F.G., F.Z., S.D., S.J., H.C., .K.N, C.W., W.D. collected and managed the study data; Y.T. carried out the analyses; Y.T., R.L. prepared the draft of the manuscript; X.Y., X.Y. revised the draft; all authors read and approved the final manuscript.

Data availability

All the data generated/analyzed during the study are available with the corresponding author on reasonable request.

Competing interests

The authors declare no competing interests.

Publisher's note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
==== Refs
References

1. Newbold T Future effects of climate and land-use change on terrestrial vertebrate community diversity under different scenarios Proc. Royal Soc. B Biol. Sci. 2018 10.1098/rspb.2018.0792
Newbold, T. Future effects of climate and land-use change on terrestrial vertebrate community diversity under different scenarios. Proc. Royal Soc. B Biol. Sci.10.1098/rspb.2018.0792 (2018).10.1098/rspb.2018.0792
2. Torre I Puig-Montserrat X Díaz M Global change effects on Mediterranean small mammal population dynamics: Demography of Algerian mice (Mus spretus) along land use and climate gradients Sci. Total Environ. 2023 10.1016/j.scitotenv.2022.160875 38056636
Torre, I., Puig-Montserrat, X. & Díaz, M. Global change effects on Mediterranean small mammal population dynamics: Demography of Algerian mice (Mus spretus) along land use and climate gradients. Sci. Total Environ.10.1016/j.scitotenv.2022.160875 (2023).38056636 10.1016/j.scitotenv.2022.160875
3. Zahoor B Projected shifts in the distribution range of Asiatic black bear (Ursus thibetanus) in the Hindu Kush Himalaya due to climate change Eco. Inform. 2021 10.1016/j.ecoinf.2021.101312
Zahoor, B. et al. Projected shifts in the distribution range of Asiatic black bear (Ursus thibetanus) in the Hindu Kush Himalaya due to climate change. Eco. Inform.10.1016/j.ecoinf.2021.101312 (2021).10.1016/j.ecoinf.2021.101312
4. Wester, P., Mishra, A., Mukherji, A. & Shrestha, A. The Hindu Kush Himalaya Assessment: Mountains, Climate Change, Sustainability and People. (2019).
5. Muñoz A-R Márquez AL Real R An approach to consider behavioral plasticity as a source of uncertainty when forecasting species' response to climate change Ecol. Evol. 2015 5 2359 2373 10.1002/ece3.1519 26120426
Muñoz, A.-R., Márquez, A. L. & Real, R. An approach to consider behavioral plasticity as a source of uncertainty when forecasting species’ response to climate change. Ecol. Evol. 5, 2359–2373. 10.1002/ece3.1519 (2015).26120426 10.1002/ece3.1519
6. Kumar V Biological clocks and regulation of seasonal reproduction and migration in birds Physiol. Biochem. Zool. 2010 83 827 835 10.1086/652243 20604684
Kumar, V. et al. Biological clocks and regulation of seasonal reproduction and migration in birds. Physiol. Biochem. Zool. 83, 827–835. 10.1086/652243 (2010).20604684 10.1086/652243
7. Halupka L Halupka K The effect of climate change on the duration of avian breeding seasons: A meta-analysis Proc. Royal Soc. B Biol. Sci. 2017 284 20171710 10.1098/rspb.2017.1710
Halupka, L. & Halupka, K. The effect of climate change on the duration of avian breeding seasons: A meta-analysis. Proc. Royal Soc. B Biol. Sci. 284, 20171710. 10.1098/rspb.2017.1710 (2017).10.1098/rspb.2017.1710
8. Bista M Panthi S Weiskopf SR Habitat overlap between Asiatic black bear Ursus thibetanus and red panda Ailurus fulgens in Himalaya PLOS ONE 2018 13 e0203697 10.1371/journal.pone.0203697 30188937
Bista, M., Panthi, S. & Weiskopf, S. R. Habitat overlap between Asiatic black bear Ursus thibetanus and red panda Ailurus fulgens in Himalaya. PLOS ONE 13, e0203697. 10.1371/journal.pone.0203697 (2018).30188937 10.1371/journal.pone.0203697
9. Bocedi G Effects of local adaptation and interspecific competition on species’ responses to climate change Ann. New York Acad. Sci. 2013 1297 83 97 10.1111/nyas.12211 23905876
Bocedi, G. et al. Effects of local adaptation and interspecific competition on species’ responses to climate change. Ann. New York Acad. Sci. 1297, 83–97. 10.1111/nyas.12211 (2013).23905876 10.1111/nyas.12211
10. Rather TA Kumar S Khan JA Multi-scale habitat selection and impacts of climate change on the distribution of four sympatric meso-carnivores using random forest algorithm Ecol. Process. 2020 9 60 10.1186/s13717-020-00265-2
Rather, T. A., Kumar, S. & Khan, J. A. Multi-scale habitat selection and impacts of climate change on the distribution of four sympatric meso-carnivores using random forest algorithm. Ecol. Process. 9, 60. 10.1186/s13717-020-00265-2 (2020).10.1186/s13717-020-00265-2
11. Zhao N Zhang X Shan G Ye X Evaluating the effects of climate change on spatial aggregation of giant pandas and sympatric species in a mountainous landscape Animals 2021 11 3332 10.3390/ani11113332 34828063
Zhao, N., Zhang, X., Shan, G. & Ye, X. Evaluating the effects of climate change on spatial aggregation of giant pandas and sympatric species in a mountainous landscape. Animals 11, 3332 (2021).34828063 10.3390/ani11113332
12. Newbold T Global effects of land use on local terrestrial biodiversity Nature 2015 520 45 50 10.1038/nature14324 25832402
Newbold, T. et al. Global effects of land use on local terrestrial biodiversity. Nature 520, 45–50. 10.1038/nature14324 (2015).25832402 10.1038/nature14324
13. Dornelas M Assemblage time series reveal biodiversity change but not systematic loss Science 2014 344 296 299 10.1126/science.1248484 24744374
Dornelas, M. et al. Assemblage time series reveal biodiversity change but not systematic loss. Science 344, 296–299. 10.1126/science.1248484 (2014).24744374 10.1126/science.1248484
14. Gibbard S Caldeira K Bala G Phillips TJ Wickett M Climate effects of global land cover change Geophys. Res. Lett. 2005 10.1029/2005GL024550
Gibbard, S., Caldeira, K., Bala, G., Phillips, T. J. & Wickett, M. Climate effects of global land cover change. Geophys. Res. Lett.10.1029/2005GL024550 (2005).10.1029/2005GL024550
15. Demétrio Luis G Ângela Schmitz P Luis Fernando Carvalho P Leonardo M Spatial and temporal patterns of waterbird assemblages in fragmented wetlands of Southern Brazil Waterbirds 2005 28 261 272 10.1675/1524-4695(2005)028[0261:SATPOW]2.0.CO;2
Demétrio Luis, G., Ângela Schmitz, P., Luis Fernando Carvalho, P. & Leonardo, M. Spatial and temporal patterns of waterbird assemblages in fragmented wetlands of Southern Brazil. Waterbirds 28, 261–272. 10.1675/1524-4695(2005)028[0261:SATPOW]2.0.CO;2 (2005).10.1675/1524-4695(2005)028[0261:SATPOW]2.0.CO;2
16. Martínez-Meyer E Pérez-Torres J Rojas-Soto O Sheppard JK Sáenz-Jiménez F Effects of climate change and human influence in the distribution and range overlap between two widely distributed avian scavengers Bird Conserv. Int. 2021 31 77 95 10.1017/S0959270920000271
Martínez-Meyer, E., Pérez-Torres, J., Rojas-Soto, O., Sheppard, J. K. & Sáenz-Jiménez, F. Effects of climate change and human influence in the distribution and range overlap between two widely distributed avian scavengers. Bird Conserv. Int. 31, 77–95. 10.1017/S0959270920000271 (2021).10.1017/S0959270920000271
17. Stenseth NC Testing for effects of climate change on competitive relationships and coexistence between two bird species Proc. Royal Soc. B Biol. Sci. 2015 10.1098/rspb.2014.1958
Stenseth, N. C. et al. Testing for effects of climate change on competitive relationships and coexistence between two bird species. Proc. Royal Soc. B Biol. Sci.10.1098/rspb.2014.1958 (2015).10.1098/rspb.2014.1958
18. Khan Z Habitat determinants of species occupancy and niche partitioning among sympatric owlets: The paradoxical role of agricultural lands for the endangered Forest Owlet, Athene blewitti J. Field Ornithol. 2023 10.5751/jfo-00244-940201
Khan, Z. et al. Habitat determinants of species occupancy and niche partitioning among sympatric owlets: The paradoxical role of agricultural lands for the endangered Forest Owlet, Athene blewitti. J. Field Ornithol.10.5751/jfo-00244-940201 (2023).10.5751/jfo-00244-940201
19. Rezaei S Mohammadi A Shadloo S Ranaie M Wan HY Climate change induces habitat shifts and overlaps among carnivores in an arid and semi-arid ecosystem Ecol. Inf. 2023 77 102247 10.1016/j.ecoinf.2023.102247
Rezaei, S., Mohammadi, A., Shadloo, S., Ranaie, M. & Wan, H. Y. Climate change induces habitat shifts and overlaps among carnivores in an arid and semi-arid ecosystem. Ecol. Inf. 77, 102247. 10.1016/j.ecoinf.2023.102247 (2023).10.1016/j.ecoinf.2023.102247
20. Archibald GW Lantis SDH Lantis LR Munetchika I Endangered ibises Threskiornithinae: Their future in the wild and in captivity Int. Zoo Yearbook 1980 20 6 17 10.1111/j.1748-1090.1980.tb00936.x
Archibald, G. W., Lantis, S. D. H., Lantis, L. R. & Munetchika, I. Endangered ibises Threskiornithinae: Their future in the wild and in captivity. Int. Zoo Yearbook 20, 6–17. 10.1111/j.1748-1090.1980.tb00936.x (1980).10.1111/j.1748-1090.1980.tb00936.x
21. Zhang Y Can we reestablish a self-sustaining population? A case study on reintroduced Crested Ibis with population viability analysis Avian Res. 2021 12 14 10.1186/s40657-021-00250-z
Zhang, Y. et al. Can we reestablish a self-sustaining population? A case study on reintroduced Crested Ibis with population viability analysis. Avian Res. 12, 14. 10.1186/s40657-021-00250-z (2021).10.1186/s40657-021-00250-z
22. Sun Y Predicting and understanding spatio-temporal dynamics of species recovery: Implications for Asian crested ibis Nipponia nippon conservation in China Divers. Distrib. 2016 22 893 904 10.1111/ddi.12460
Sun, Y. et al. Predicting and understanding spatio-temporal dynamics of species recovery: Implications for Asian crested ibis Nipponia nippon conservation in China. Divers. Distrib. 22, 893–904. 10.1111/ddi.12460 (2016).10.1111/ddi.12460
23. Ye Y What makes a tactile forager join mixed-species flocks? A case study with the endangered Crested Ibis (Nipponia nippon) The Auk 2017 134 421 431 10.1642/AUK-16-191.1
Ye, Y. et al. What makes a tactile forager join mixed-species flocks? A case study with the endangered Crested Ibis (Nipponia nippon). The Auk 134, 421–431. 10.1642/AUK-16-191.1 (2017).10.1642/AUK-16-191.1
24. Wang M On the sustainability of a reintroduced Crested Ibis population in Qinling Mountains, Shaanxi Central China Restor. Ecol. 2017 25 261 268 10.1111/rec.12409
Wang, M. et al. On the sustainability of a reintroduced Crested Ibis population in Qinling Mountains, Shaanxi Central China. Restor. Ecol. 25, 261–268. 10.1111/rec.12409 (2017).10.1111/rec.12409
25. Yu F Linking the past and present to predict the distribution of Asian crested ibis (Nipponia nippon) under global changes Integr. Zool. 2022 17 1095 1105 10.1111/1749-4877.12581 34496146
Yu, F. et al. Linking the past and present to predict the distribution of Asian crested ibis (Nipponia nippon) under global changes. Integr. Zool. 17, 1095–1105. 10.1111/1749-4877.12581 (2022).34496146 10.1111/1749-4877.12581
26. Macarthur R Levins R The limiting similarity, convergence, and divergence of coexisting species Am. Nat. 1967 101 377 385 10.1086/282505
Macarthur, R. & Levins, R. The limiting similarity, convergence, and divergence of coexisting species. Am. Nat. 101, 377–385. 10.1086/282505 (1967).10.1086/282505
27. Broennimann O Measuring ecological niche overlap from occurrence and spatial environmental data Global Ecol. Biogeogr. 2012 21 481 497 10.1111/j.1466-8238.2011.00698.x
Broennimann, O. et al. Measuring ecological niche overlap from occurrence and spatial environmental data. Global Ecol. Biogeogr. 21, 481–497. 10.1111/j.1466-8238.2011.00698.x (2012).10.1111/j.1466-8238.2011.00698.x
28. May RM Arthur RHM Niche overlap as a function of environmental variability Proc. Nat. Acad. Sci. 1972 69 1109 1113 10.1073/pnas.69.5.1109 4504323
May, R. M. & Arthur, R. H. M. Niche overlap as a function of environmental variability. Proc. Nat. Acad. Sci. 69, 1109–1113. 10.1073/pnas.69.5.1109 (1972).4504323 10.1073/pnas.69.5.1109
29. Sahade R Climate change and glacier retreat drive shifts in an Antarctic benthic ecosystem Sci. Adv. 2015 1 e1500050 10.1126/sciadv.1500050 26702429
Sahade, R. et al. Climate change and glacier retreat drive shifts in an Antarctic benthic ecosystem. Sci. Adv. 1, e1500050. 10.1126/sciadv.1500050 (2015).26702429 10.1126/sciadv.1500050
30. Hardin G The competitive exclusion principle Science 1960 131 1292 1297 10.1126/science.131.3409.1292 14399717
Hardin, G. The competitive exclusion principle. Science 131, 1292–1297. 10.1126/science.131.3409.1292 (1960).14399717 10.1126/science.131.3409.1292
31. Johnson CA Bronstein JL Coexistence and competitive exclusion in mutualism Ecology 2019 100 e02708 10.1002/ecy.2708 30924140
Johnson, C. A. & Bronstein, J. L. Coexistence and competitive exclusion in mutualism. Ecology 100, e02708. 10.1002/ecy.2708 (2019).30924140 10.1002/ecy.2708
32. Pianka ER Niche overlap and diffuse competition Proc. National Acad. Sci. 1974 71 2141 2145 10.1073/pnas.71.5.2141
Pianka, E. R. Niche overlap and diffuse competition. Proc. National Acad. Sci. 71, 2141–2145. 10.1073/pnas.71.5.2141 (1974).10.1073/pnas.71.5.2141
33. Rojas LM McNeil R Cabana T Lachapelle P Behavioral, morphological and physiological correlates of diurnal and nocturnal vision in selected wading bird species Brain Behav. Evol. 1999 53 227 242 10.1159/000006596 10473901
Rojas, L. M., McNeil, R., Cabana, T. & Lachapelle, P. Behavioral, morphological and physiological correlates of diurnal and nocturnal vision in selected wading bird species. Brain Behav. Evol. 53, 227–242. 10.1159/000006596 (1999).10473901 10.1159/000006596
34. Beerens JM Gawlik DE Herring G Cook MI Dynamic habitat selection by two wading bird species with divergent foraging strategies in a seasonally fluctuating wetland The Auk 2011 128 651 662 10.1525/auk.2011.10165
Beerens, J. M., Gawlik, D. E., Herring, G. & Cook, M. I. Dynamic habitat selection by two wading bird species with divergent foraging strategies in a seasonally fluctuating wetland. The Auk 128, 651–662. 10.1525/auk.2011.10165 (2011).10.1525/auk.2011.10165
35. Ma, L. et al. Changes in the habitat preference of Crested Ibis (Nipponia nippon) during a period of rapid population increase. Animals 11 (2021).
36. Vardi R Kotler BP Altstein O Abramsky Z Social behaviour and foraging success of little egrets (Egretta garzetta) Behav. Process. 2021 183 104318 10.1016/j.beproc.2021.104318
Vardi, R., Kotler, B. P., Altstein, O. & Abramsky, Z. Social behaviour and foraging success of little egrets (Egretta garzetta). Behav. Process. 183, 104318. 10.1016/j.beproc.2021.104318 (2021).10.1016/j.beproc.2021.104318
37. Rani M Health risk assessment by toxic metals in little egrets (Egretta garzetta) and food chain contaminations Saudi J. Biol. Sci. 2022 29 366 370 10.1016/j.sjbs.2021.08.106 35002431
Rani, M. et al. Health risk assessment by toxic metals in little egrets (Egretta garzetta) and food chain contaminations. Saudi J. Biol. Sci. 29, 366–370. 10.1016/j.sjbs.2021.08.106 (2022).35002431 10.1016/j.sjbs.2021.08.106
38. Fasola M Changes in rice cultivation affect population trends of herons and egrets in Italy Global Ecol. Conserv. 2022 36 e02135 10.1016/j.gecco.2022.e02135
Fasola, M. et al. Changes in rice cultivation affect population trends of herons and egrets in Italy. Global Ecol. Conserv. 36, e02135. 10.1016/j.gecco.2022.e02135 (2022).10.1016/j.gecco.2022.e02135
39. Boulangeat I Lavergne S Van Es J Garraud L Thuiller W Niche breadth, rarity and ecological characteristics within a regional flora spanning large environmental gradients J. Biogeogr. 2012 39 204 214 10.1111/j.1365-2699.2011.02581.x
Boulangeat, I., Lavergne, S., Van Es, J., Garraud, L. & Thuiller, W. Niche breadth, rarity and ecological characteristics within a regional flora spanning large environmental gradients. J. Biogeogr. 39, 204–214. 10.1111/j.1365-2699.2011.02581.x (2012).10.1111/j.1365-2699.2011.02581.x
40. Slatyer RA Hirst M Sexton JP Niche breadth predicts geographical range size: A general ecological pattern Ecol. Lett. 2013 16 1104 1114 10.1111/ele.12140 23773417
Slatyer, R. A., Hirst, M. & Sexton, J. P. Niche breadth predicts geographical range size: A general ecological pattern. Ecol. Lett. 16, 1104–1114. 10.1111/ele.12140 (2013).23773417 10.1111/ele.12140
41. Li X Tian H Li D Why the crested ibis declined in the middle twentieth century Biodivers. Conserv. 2009 18 2165 2172 10.1007/s10531-009-9580-z
Li, X., Tian, H. & Li, D. Why the crested ibis declined in the middle twentieth century. Biodivers. Conserv. 18, 2165–2172. 10.1007/s10531-009-9580-z (2009).10.1007/s10531-009-9580-z
42. Feng S The genomic footprints of the fall and recovery of the Crested ibis Curr. Biol. 2019 29 340 349.e347 10.1016/j.cub.2018.12.008 30639104
Feng, S. et al. The genomic footprints of the fall and recovery of the Crested ibis. Curr. Biol. 29, 340-349.e347. 10.1016/j.cub.2018.12.008 (2019).30639104 10.1016/j.cub.2018.12.008
43. Qiao H Saupe EE Soberón J Peterson AT Myers CE Impacts of niche breadth and dispersal ability on macroevolutionary patterns Am. Nat. 2016 188 149 162 10.1086/687201 27420781
Qiao, H., Saupe, E. E., Soberón, J., Peterson, A. T. & Myers, C. E. Impacts of niche breadth and dispersal ability on macroevolutionary patterns. Am. Nat. 188, 149–162. 10.1086/687201 (2016).27420781 10.1086/687201
44. Batstone RT Carscadden KA Afkhami ME Frederickson ME Using niche breadth theory to explain generalization in mutualisms Ecology 2018 99 1039 1050 10.1002/ecy.2188 29453827
Batstone, R. T., Carscadden, K. A., Afkhami, M. E. & Frederickson, M. E. Using niche breadth theory to explain generalization in mutualisms. Ecology 99, 1039–1050. 10.1002/ecy.2188 (2018).29453827 10.1002/ecy.2188
45. Robin AB Nathan JD Mark IC Nestling diet of three sympatrically nesting wading bird species in the Florida Everglades Waterbirds 2012 35 154 159 10.1675/063.035.0116
Robin, A. B., Nathan, J. D. & Mark, I. C. Nestling diet of three sympatrically nesting wading bird species in the Florida Everglades. Waterbirds 35, 154–159. 10.1675/063.035.0116 (2012).10.1675/063.035.0116
46. Carlos M Food and niche overlap of the scarlet ibis and the yellow-crowned night heron in a tropical mangrove swamp Waterbirds 2004 27 1 8 10.1675/1524-4695(2004)027[0001:FANOOT]2.0.CO;2
Carlos, M. Food and niche overlap of the scarlet ibis and the yellow-crowned night heron in a tropical mangrove swamp. Waterbirds 27, 1–8. 10.1675/1524-4695(2004)027[0001:FANOOT]2.0.CO;2 (2004).10.1675/1524-4695(2004)027[0001:FANOOT]2.0.CO;2
47. Huang Y Tall trees drive the nest-site selection of wild Crested Ibis Nipponia nippon Bird Conserv. Int. 2022 32 486 497 10.1017/S0959270921000526
Huang, Y. et al. Tall trees drive the nest-site selection of wild Crested Ibis Nipponia nippon. Bird Conserv. Int. 32, 486–497. 10.1017/S0959270921000526 (2022).10.1017/S0959270921000526
48. Li Z Competitive interactions and coexistence of sympatric flagship carnivores in Asia Integr. Zool. n/a 2023 10.1111/1749-4877.12724
Li, Z. et al. Competitive interactions and coexistence of sympatric flagship carnivores in Asia. Integr. Zool. n/a10.1111/1749-4877.12724 (2023).10.1111/1749-4877.12724
49. Carter NH Linnell JDC Co-adaptation is key to coexisting with large carnivores Trends Ecol. Evol. 2016 31 575 578 10.1016/j.tree.2016.05.006 27377600
Carter, N. H. & Linnell, J. D. C. Co-adaptation is key to coexisting with large carnivores. Trends Ecol. Evol. 31, 575–578. 10.1016/j.tree.2016.05.006 (2016).27377600 10.1016/j.tree.2016.05.006
50. Kriegler E The need for and use of socio-economic scenarios for climate change analysis: A new approach based on shared socio-economic pathways Global Environ. Change 2012 22 807 822 10.1016/j.gloenvcha.2012.05.005
Kriegler, E. et al. The need for and use of socio-economic scenarios for climate change analysis: A new approach based on shared socio-economic pathways. Global Environ. Change 22, 807–822. 10.1016/j.gloenvcha.2012.05.005 (2012).10.1016/j.gloenvcha.2012.05.005
51. Sillero N Want to model a species niche? A step-by-step guideline on correlative ecological niche modelling Ecol. Model. 2021 456 109671 10.1016/j.ecolmodel.2021.109671
Sillero, N. et al. Want to model a species niche? A step-by-step guideline on correlative ecological niche modelling. Ecol. Model. 456, 109671. 10.1016/j.ecolmodel.2021.109671 (2021).10.1016/j.ecolmodel.2021.109671
52. Liang X Understanding the drivers of sustainable land expansion using a patch-generating land use simulation (PLUS) model: A case study in Wuhan, China Comput. Environ. Urban Syst. 2021 10.1016/j.compenvurbsys.2020.101569
Liang, X. et al. Understanding the drivers of sustainable land expansion using a patch-generating land use simulation (PLUS) model: A case study in Wuhan, China. Comput. Environ. Urban Syst.10.1016/j.compenvurbsys.2020.101569 (2021).10.1016/j.compenvurbsys.2020.101569
53. Zhai H Understanding spatio-temporal patterns of land use/land cover change under urbanization in Wuhan, China, 2000–2019 Remote Sens. 2021 10.3390/rs13163331
Zhai, H. et al. Understanding spatio-temporal patterns of land use/land cover change under urbanization in Wuhan, China, 2000–2019. Remote Sens.10.3390/rs13163331 (2021).10.3390/rs13163331
54. Li C Multi-scenario simulation of ecosystem service value for optimization of land use in the Sichuan-Yunnan ecological barrier China Ecol. Indic. 2021 10.1016/j.ecolind.2021.108328
Li, C. et al. Multi-scenario simulation of ecosystem service value for optimization of land use in the Sichuan-Yunnan ecological barrier China. Ecol. Indic.10.1016/j.ecolind.2021.108328 (2021).10.1016/j.ecolind.2021.108328
55. Sun S Mu L Feng R Wang L He J GAN-based LUCC prediction via the combination of prior city planning information and land-use probability IEEE J. Sel. Top. Appl. Earth Observ. Remote Sens. 2021 14 10189 10198 10.1109/jstars.2021.3106481
Sun, S., Mu, L., Feng, R., Wang, L. & He, J. GAN-based LUCC prediction via the combination of prior city planning information and land-use probability. IEEE J. Sel. Top. Appl. Earth Observ. Remote Sens. 14, 10189–10198. 10.1109/jstars.2021.3106481 (2021).10.1109/jstars.2021.3106481
56. Wang Z Dynamic simulation of land use change and assessment of carbon storage based on climate change scenarios at the city level: A case Ecol. Indic. 2022 10.1016/j.ecolind.2021.108499
Wang, Z. et al. Dynamic simulation of land use change and assessment of carbon storage based on climate change scenarios at the city level: A case. Ecol. Indic.10.1016/j.ecolind.2021.108499 (2022).10.1016/j.ecolind.2021.108499
57. Zhang S Land use/land cover prediction and analysis of the middle reaches of the Yangtze River under different scenarios Sci. Total Environ. 2022 833 155238 10.1016/j.scitotenv.2022.155238 35427604
Zhang, S. et al. Land use/land cover prediction and analysis of the middle reaches of the Yangtze River under different scenarios. Sci. Total Environ. 833, 155238. 10.1016/j.scitotenv.2022.155238 (2022).35427604 10.1016/j.scitotenv.2022.155238
58. Chen J Global land cover mapping at 30m resolution: A POK-based operational approach ISPRS J. Photogramm. Remote Sens. 2015 103 7 27 10.1016/j.isprsjprs.2014.09.002
Chen, J. et al. Global land cover mapping at 30m resolution: A POK-based operational approach. ISPRS J. Photogramm. Remote Sens. 103, 7–27. 10.1016/j.isprsjprs.2014.09.002 (2015).10.1016/j.isprsjprs.2014.09.002
59. Brennan RL Prediger DJ Coefficient kappa: Some uses, misuses, and alternatives Educ. Psychol. Meas. 1981 41 687 699 10.1177/001316448104100307
Brennan, R. L. & Prediger, D. J. Coefficient kappa: Some uses, misuses, and alternatives. Educ. Psychol. Meas. 41, 687–699. 10.1177/001316448104100307 (1981).10.1177/001316448104100307
60. Huang D Huang J Liu T Delimiting urban growth boundaries using the CLUE-S model with village administrative boundaries Land Use Policy 2019 82 422 435 10.1016/j.landusepol.2018.12.028
Huang, D., Huang, J. & Liu, T. Delimiting urban growth boundaries using the CLUE-S model with village administrative boundaries. Land Use Policy 82, 422–435. 10.1016/j.landusepol.2018.12.028 (2019).10.1016/j.landusepol.2018.12.028
61. Lin W Sun Y Nijhuis S Wang Z Scenario-based flood risk assessment for urbanizing deltas using future land-use simulation (FLUS): Guangzhou Metropolitan Area as a case study Sci. Total Environ. 2020 739 139899 10.1016/j.scitotenv.2020.139899 32540659
Lin, W., Sun, Y., Nijhuis, S. & Wang, Z. Scenario-based flood risk assessment for urbanizing deltas using future land-use simulation (FLUS): Guangzhou Metropolitan Area as a case study. Sci. Total Environ. 739, 139899. 10.1016/j.scitotenv.2020.139899 (2020).32540659 10.1016/j.scitotenv.2020.139899
62. Andrade AFA Velazco SJE De Marco Júnior P ENMTML: An R package for a straightforward construction of complex ecological niche models Environ. Model. Softw. 2020 125 104615 10.1016/j.envsoft.2019.104615
Andrade, A. F. A., Velazco, S. J. E. & De Marco Júnior, P. ENMTML: An R package for a straightforward construction of complex ecological niche models. Environ. Model. Softw. 125, 104615. 10.1016/j.envsoft.2019.104615 (2020).10.1016/j.envsoft.2019.104615
63. Senay SD Worner SP Ikeda T Novel three-step pseudo-absence selection technique for improved species distribution modelling PLOS ONE 2013 8 e71218 10.1371/journal.pone.0071218 23967167
Senay, S. D., Worner, S. P. & Ikeda, T. Novel three-step pseudo-absence selection technique for improved species distribution modelling. PLOS ONE 8, e71218. 10.1371/journal.pone.0071218 (2013).23967167 10.1371/journal.pone.0071218
64. Lobo JM Jiménez-Valverde A Hortal J The uncertain nature of absences and their importance in species distribution modelling Ecography 2010 33 103 114 10.1111/j.1600-0587.2009.06039.x
Lobo, J. M., Jiménez-Valverde, A. & Hortal, J. The uncertain nature of absences and their importance in species distribution modelling. Ecography 33, 103–114. 10.1111/j.1600-0587.2009.06039.x (2010).10.1111/j.1600-0587.2009.06039.x
65. Pineda E Lobo JM Assessing the accuracy of species distribution models to predict amphibian species richness patterns J. Animal Ecol. 2009 78 182 190 10.1111/j.1365-2656.2008.01471.x
Pineda, E. & Lobo, J. M. Assessing the accuracy of species distribution models to predict amphibian species richness patterns. J. Animal Ecol. 78, 182–190. 10.1111/j.1365-2656.2008.01471.x (2009).10.1111/j.1365-2656.2008.01471.x
66. Boyce MS Vernier PR Nielsen SE Schmiegelow FKA Evaluating resource selection functions Ecol. Model. 2002 157 281 300 10.1016/S0304-3800(02)00200-4
Boyce, M. S., Vernier, P. R., Nielsen, S. E. & Schmiegelow, F. K. A. Evaluating resource selection functions. Ecol. Model. 157, 281–300. 10.1016/S0304-3800(02)00200-4 (2002).10.1016/S0304-3800(02)00200-4
67. Leroy B Without quality presence–absence data, discrimination metrics such as TSS can be misleading measures of model performance J. Biogeogr. 2018 45 1994 2002 10.1111/jbi.13402
Leroy, B. et al. Without quality presence–absence data, discrimination metrics such as TSS can be misleading measures of model performance. J. Biogeogr. 45, 1994–2002. 10.1111/jbi.13402 (2018).10.1111/jbi.13402
68. Bell JF Fielding AH A review of methods for the assessment of prediction errors in conservation presence/absence models Environm. Conserv. 1997 24 38 49 10.1017/S0376892997000088
Bell, J. F. & Fielding, A. H. A review of methods for the assessment of prediction errors in conservation presence/absence models. Environm. Conserv. 24, 38–49. 10.1017/S0376892997000088 (1997).10.1017/S0376892997000088
69. Thuiller W Lafourcade B Engler R Araújo MB BIOMOD—a platform for ensemble forecasting of species distributions Ecography 2009 32 369 373 10.1111/j.1600-0587.2008.05742.x
Thuiller, W., Lafourcade, B., Engler, R. & Araújo, M. B. BIOMOD—a platform for ensemble forecasting of species distributions. Ecography 32, 369–373. 10.1111/j.1600-0587.2008.05742.x (2009).10.1111/j.1600-0587.2008.05742.x
70. Di Cola V ecospat: An R package to support spatial analyses and modeling of species niches and distributions Ecography 2017 40 774 787 10.1111/ecog.02671
Di Cola, V. et al. ecospat: An R package to support spatial analyses and modeling of species niches and distributions. Ecography 40, 774–787. 10.1111/ecog.02671 (2017).10.1111/ecog.02671
71. Brown JL Bennett JR French CM SDMtoolbox 2.0: The next generation Python-based GIS toolkit for landscape genetic, biogeographic and species distribution model analyses PeerJ 2017 5 E4095 10.7717/peerj.4095 29230356
Brown, J. L., Bennett, J. R. & French, C. M. SDMtoolbox 2.0: The next generation Python-based GIS toolkit for landscape genetic, biogeographic and species distribution model analyses. PeerJ 5, E4095. 10.7717/peerj.4095 (2017).29230356 10.7717/peerj.4095
72. Warren DL Glor RE Turelli M ENMTools: A toolbox for comparative studies of environmental niche models Ecography 2010 33 607 611 10.1111/j.1600-0587.2009.06142.x
Warren, D. L., Glor, R. E. & Turelli, M. ENMTools: A toolbox for comparative studies of environmental niche models. Ecography 33, 607–611. 10.1111/j.1600-0587.2009.06142.x (2010).10.1111/j.1600-0587.2009.06142.x
73. Schoener TW The Anolis lizards of Bimini: Resource partitioning in a complex fauna Ecology 1968 49 704 726 10.2307/1935534
Schoener, T. W. The Anolis lizards of Bimini: Resource partitioning in a complex fauna. Ecology 49, 704–726. 10.2307/1935534 (1968).10.2307/1935534
74. Warren DL Glor RE Turelli M Environmental niche equivalency versus conservatism: Quantitative approaches to niche evolution Evolution 2008 62 2868 2883 10.1111/j.1558-5646.2008.00482.x 18752605
Warren, D. L., Glor, R. E. & Turelli, M. Environmental niche equivalency versus conservatism: Quantitative approaches to niche evolution. Evolution 62, 2868–2883. 10.1111/j.1558-5646.2008.00482.x (2008).18752605 10.1111/j.1558-5646.2008.00482.x
75. Regmi S Niche breadth and overlap of spotted deer and domestic cattle with swamp deer in tropical region of Nepal Ecol. Process. 2022 11 22 10.1186/s13717-022-00367-z
Regmi, S. et al. Niche breadth and overlap of spotted deer and domestic cattle with swamp deer in tropical region of Nepal. Ecol. Process. 11, 22. 10.1186/s13717-022-00367-z (2022).10.1186/s13717-022-00367-z
76. Mao J-F Wang X-R Distinct Niche divergence characterizes the homoploid hybrid speciation of Pinus densata on the Tibetan Plateau Am. Nat. 2011 177 424 439 10.1086/658905 21460565
Mao, J.-F. & Wang, X.-R. Distinct Niche divergence characterizes the homoploid hybrid speciation of Pinus densata on the Tibetan Plateau. Am. Nat. 177, 424–439. 10.1086/658905 (2011).21460565 10.1086/658905
