Journal of Environmental Science Studies

Journal of Environmental Science Studies

Identifying Conservation Hotspots for the Medicinal Plant Teucrium orientale under the Interactive Effects of Climate Change and Human Pressures

Document Type : Original Article

Authors
1 Associate Professor, Department of Natural Engineering, Faculty of Natural Resources and Earth Sciences, Shahrekord University, Shahrekord, Iran.
2 MS graduate, Department of Natural Engineering, Faculty of Natural Resources and Earth Sciences, Shahrekord University, Shahrekord, Iran.
3 Isfahan Agriculture and Natural Resources Research and Education Center, Isfahan, Iran.
4 Associate Professor, Department of Environmental Engineering, Faculty of Natural Resources and Earth Sciences, Shahrekord University, Shahrekord, Iran.
10.22034/jess.2026.581830.2456
Abstract
Abstract
Given the high vulnerability of arid regions to climate change and human pressures, this study investigated the current and future habitat suitability of Teucrium orientale in Isfahan Province to provide a scientific basis for conservation planning in these environmentally stressed ecosystems. To achieve the objectives of the study, a comprehensive set of advanced species distribution modeling algorithms was implemented within an ensemble modeling framework. The models predicted species distribution using CHELSA bioclimatic variables in combination with topographic parameters. Future habitat changes were assessed using outputs from the GFDL-ESM4 and MRI-ESM2.0 climate models under the SSP126, SSP370, and SSP585 emission scenarios, while also incorporating human pressure indices into the analysis. The results demonstrated that climatic factors, particularly temperature and precipitation, exert the strongest influence on the occurrence and distribution of this species. Simulations of future conditions further indicated that, under most scenarios, the extent of suitable habitats for the species is expected to decline, although the spatial pattern of these changes varies among climate models across the province. Among different regions, the western and southern parts of the province, especially high-elevation areas, showed greater potential for maintaining favorable environmental conditions in the future and were therefore identified as potential climatic refugia. However, spatial analysis revealed that a considerable proportion of these areas overlaps with regions experiencing substantial human activity, which may threaten the long-term stability of the remaining suitable habitats.Accordingly, effective land-use management, reduction of anthropogenic pressures, and targeted conservation of climatically stable areas could play a crucial role in preserving both the species and its habitats. The findings of this study also highlight the necessity of simultaneously considering climatic and human-related factors in conservation planning in order to support more scientifically robust decision-making processes.

EXTENDED ABSTRACT

Introduction

Climate change significantly alters ecosystem structures and plant distribution patterns, forcing species to shift their geographic ranges in response to changing temperature and precipitation regimes. Rangeland ecosystems in arid and semi-arid regions are particularly susceptible to these shifts, as existing water limitations are often exacerbated by both climatic stress and anthropogenic factors like overgrazing and land-use change. Consequently, utilizing predictive modeling to assess future habitat suitability has become a crucial tool for identifying areas where species can survive, which is essential for developing effective conservation strategies. This specific study evaluates the current and projected habitat of Teucrium orientale in Isfahan Province to provide a scientific foundation for sustainable rangeland management and biodiversity preservation in environmentally stressed landscapes.

Materials and methods
In this study, 30 occurrence records of T. orientale were used to model its habitat suitability. Nineteen bioclimatic variables for three time periods (current conditions, 2050, and 2080) were obtained from the CHELSA database. These climatic data were derived from two general circulation models, GFDL ESM4 and MRI ESM2 0. Future projections were evaluated under three emission scenarios: SSP126, SSP370, and SSP585. All variables were prepared at a high spatial resolution to improve modeling accuracy. In addition to climatic variables, physiographic factors including elevation, slope, and aspect were incorporated. Pearson’s correlation analysis was used to identify and remove highly collinear variables. This step helped reduce redundancy and improve model performance. After variable selection, ten predictors (nine bioclimatic variables and one physiographic variable) were retained. These predictors were used to run species distribution models using eight different algorithms. The outputs of these algorithms were combined within an ensemble modeling framework. Finally, the human pressure index map was overlaid with the predicted optimal habitat map of T. orientale to assess potential human impacts.


Results and discussion
The results indicated that climatic variables, particularly temperature and precipitation, play the most important role in determining the occurrence and distribution of this species. The net rate of change in habitat extent relative to current conditions was estimated to range from −24.21% to −89.84%. Simulations of future conditions suggested that under most scenarios the area of suitable habitats for this species will decline, although the spatial pattern of these changes varies across the province depending on the climate model used. Among different regions, the western and southern parts of the province, especially high elevation areas, demonstrated a greater capacity to maintain favorable environmental conditions in the future and were therefore identified as potential habitat refugia. However, spatial analysis revealed that a considerable portion of these areas overlaps with regions experiencing substantial human activities, which may pose challenges for the persistence of the remaining suitable habitats.

Conclusion
Effective land use management is vital for mitigating climate change impacts on vulnerable species and rangelands. By reducing anthropogenic pressures and protecting climatically stable refugia, conservationists can enhance ecosystem resilience and ensure long-term biodiversity. Achieving successful outcomes requires integrating ecological modeling into land-use planning and fostering collaboration between policymakers and local communities. Ultimately, robust conservation strategies must account for both climatic and human-related drivers to ensure scientifically informed and sustainable management.

EXTENDED ABSTRACT

Introduction

Climate change significantly alters ecosystem structures and plant distribution patterns, forcing species to shift their geographic ranges in response to changing temperature and precipitation regimes. Rangeland ecosystems in arid and semi-arid regions are particularly susceptible to these shifts, as existing water limitations are often exacerbated by both climatic stress and anthropogenic factors like overgrazing and land-use change. Consequently, utilizing predictive modeling to assess future habitat suitability has become a crucial tool for identifying areas where species can survive, which is essential for developing effective conservation strategies. This specific study evaluates the current and projected habitat of Teucrium orientale in Isfahan Province to provide a scientific foundation for sustainable rangeland management and biodiversity preservation in environmentally stressed landscapes.

Materials and methods
In this study, 30 occurrence records of T. orientale were used to model its habitat suitability. Nineteen bioclimatic variables for three time periods (current conditions, 2050, and 2080) were obtained from the CHELSA database. These climatic data were derived from two general circulation models, GFDL ESM4 and MRI ESM2 0. Future projections were evaluated under three emission scenarios: SSP126, SSP370, and SSP585. All variables were prepared at a high spatial resolution to improve modeling accuracy. In addition to climatic variables, physiographic factors including elevation, slope, and aspect were incorporated. Pearson’s correlation analysis was used to identify and remove highly collinear variables. This step helped reduce redundancy and improve model performance. After variable selection, ten predictors (nine bioclimatic variables and one physiographic variable) were retained. These predictors were used to run species distribution models using eight different algorithms. The outputs of these algorithms were combined within an ensemble modeling framework. Finally, the human pressure index map was overlaid with the predicted optimal habitat map of T. orientale to assess potential human impacts.


Results and discussion
The results indicated that climatic variables, particularly temperature and precipitation, play the most important role in determining the occurrence and distribution of this species. The net rate of change in habitat extent relative to current conditions was estimated to range from −24.21% to −89.84%. Simulations of future conditions suggested that under most scenarios the area of suitable habitats for this species will decline, although the spatial pattern of these changes varies across the province depending on the climate model used. Among different regions, the western and southern parts of the province, especially high elevation areas, demonstrated a greater capacity to maintain favorable environmental conditions in the future and were therefore identified as potential habitat refugia. However, spatial analysis revealed that a considerable portion of these areas overlaps with regions experiencing substantial human activities, which may pose challenges for the persistence of the remaining suitable habitats.

Conclusion
Effective land use management is vital for mitigating climate change impacts on vulnerable species and rangelands. By reducing anthropogenic pressures and protecting climatically stable refugia, conservationists can enhance ecosystem resilience and ensure long-term biodiversity. Achieving successful outcomes requires integrating ecological modeling into land-use planning and fostering collaboration between policymakers and local communities. Ultimately, robust conservation strategies must account for both climatic and human-related drivers to ensure scientifically informed and sustainable management.
Keywords

1.      Araújo, M.B., Anderson, R.P., Barbosa, A.M., Beale, C.M., Dormann, C.F., Early, R., Garcia, R.A., Guisan, A., Maiorano, L., Naimi, B., O’Hara, R.B., Zimmermann, N.E., Rahbek, C. 2019. Standards for distribution models in biodiversity assessments. Science Advances, 5(1), eaat4858. https://doi.org/10.1126/sciadv.aat4858
2.      Ardestani, E.G., Tarkesh, M., Bassiri, M., Vahabi, M.R. 2014. Potential habitat modeling for reintroduction of three native plant species in central Iran. Arid Land, 7: 381–390.
4.      Cahyaningsih, R., Phillips, J., Brehm, J.M., Gaisberger, H., Maxted, N. 2021. Climate change impact on medicinal plants in Indonesia. Global Ecology and Conservation, 30: e01752.
5.      Fick, S.E., Hijmans, R.J. 2017. WorldClim 2: new 1-km spatial resolution climate surfaces for global land areas. International Journal of Climatology, 37: 4302–4315.
7.      Ghehsareh Ardestani, E., Heidari Ghahfarrokhi, Z. 2021. Ensembpecies distribution modeling of Salvia hydrangea under future climate change scenarios in Central Zagros Mountains Iran. Glob. Ecol. Conserv., 26: e01488.
8.      Guisan, A., Thuiller, W. Zimmermann, N. 2017. Habitat suitability and distribution models: with applications in R. Cambridge, United Kingdom, 462p.
9.      Hao, T., Elith, J., Lahoz‑Monfort, J.J., Guillera‑Arroita, G. 2019. A review of evidence about use and performance of species distribution modelling ensembles like BIOMOD. Diversity and Distributions, 25(5): 839–852. https://doi.org/10.1111/ddi.12892
10.  Hijmans, R.J., Cameron, S.E., Parra, J.L., Jones, P.G., Jarvis, A. 2005. Very high resolution interpolated climate surfaces for global land areas. International Journal of Climatology: A Journal of the Royal Meteorological Society, 25(15): 1965-1978.
11.  Karami, S., Ejtehadi, H., Moazzeni, H., Vaezi, J. Behroozian, M. 2022. Minimal climate change impacts on the geographic distribution of Nepeta glomerulosa, medicinal species endemic to southwestern and central Asia. Scientific Reports, 12: 19893.
12.  Kunwar, R.M., Thapa-Magar, K.B., Subedi, S.C., Kutal, D.H., Baral, B., Joshi, N.R., Adhikari, B., Upadhyaya, K.S., Thapa-Magar, S., Ansari, A.S., Thapa, G.J. and Bhandari, A. R., 2023. Distribution of important medicinal plant species in Nepal under past, present, and future climatic conditions. Ecological Indicators. 146: 1-12, https://doi.org/ 10.1016/j.ecolind.2023.109879
13.  Lee, C.O., Arge, C.N., Odstril, D., Millward, G., Pizzo, V., Quinn, J.M., Henney, C.J. 2013. Ensemble modeling of CME Propagation. Solar Physics, 285(1-2): 349 368.
14.  Lin, C.T., Chiu, C.A. 2018. The relic Trochodendron aralioides Siebold & Zucc. (Trochodendraceae) in Taiwan: Ensemble distribution modeling and climate change impacts. Forests, 9(10): 1–17. https://doi.org/10.3390/f9100620
15.  Mohapatra, J., Singh, C.P., Hamid, M., Verma, A., Semwal, S.C., Gajmer, B., Khuroo, A.A., Kumar, A., Nautiyal, M.C., Sharma, N., Pandya, H.A. 2019. Modelling Betulautilis distribution in response to climate-warming scenarios in Hindu-Kush Himalaya using random forest. Biodiversity and Conservation, 28 (8): 2295-2317.
16.  Moraitis, M.L., Valavanis, V.D., Karakassis, I. 2019. Modelling the effects of climate change on the distribution of benthic indicator species in the Eastern Mediterranean Sea. Science of The Total Environment, 667:16–24.
17.  Mu H. Li, X. Wen Y. Huang J. Du P. SuW. Miao S. Geng M. 2022. A global record of annual terrestrial Human Footprint datasetfrom 2000 to 2018.Scientific Data, 9(1): 1-26.
18.  Noroozi, J., Körner, C., Doostmohammadi, M., Akbarinia, M., Khoshnazar, H., Mosadegh, H. 2019. Endemic diversity and distribution of the Iranian vascular flora across phytogeographical regions, biodiversity hotspots and areas of endemism. Scientific Reports, 9: 12991. https://doi.org/10.1038/s41598-019-49417-1
19.  Noroozi, J., Moser, D., Essl, F. 2016. Diversity, distribution, ecology and description rates of alpine endemic plant species from Iranian mountains. Alpine Botany, 126: 1–9.
20.  Rew, J., Cho, Y., Hwang, J., Kim, S., Lee, W. 2020. Habitat suitability estimation using a two‑stage ensemble approach. Remote Sensing, 12(9): 1475. https://doi.org/10.3390/rs12091475
22.  Salas, E.A.L., Valdez, R., Michel, S., Boykin, K.G. 2018. Habitat assessment of Marco Polo sheep (Ovis ammon polii) in Eastern Tajikistan: Modeling the effects of climate change. Ecology and Evolution, 8(10): 5124-5138.
23.  Tang, C.Q., Dong, Y.F., Herrando-Moraira, S., Matsui, T., Ohashi, H., He, L.Y., Nakao, K., Tanaka, N., Tomita, M., Li, X.S., Yan, H.Z. 2017. Potential effects of climate change on geographic distribution of the tertiary relict tree species Davidia involucrata in China. Scientific Reports, 7: 43822.
24.  Tarkesh, M., Jetschke, G. 2016. Investigation of current and future potential distribution of Astragalus gossypinus in Central Iran using species distribution modelling. Arabian Journal of Geosciences, 9(1): 1-11.
26.  Thuiller, W., Lafourcade, B., Engler, R. Araújo, M. B. 2009. BIOMOD - a platform for ensemble forecasting of species distributions. Ecography (Cop.), 32: 369–373.
27.  Thuiller, W., Lavorel, S., Araujo, M. B., Sykes, M. T., Prentice, I. C. 2005. Climate change threats to plant diversity in Europe. Proceedings of the National Academy of Sciences, 102(23): 8245-8250. https://doi.org/10.1073/pnas.0409902102
28.  Travis, J.M.J., Delgado, M., Bocedi, G., Baguette, M., Barton, K., Bonte, D., Boulangeat, I., Hodgson, J.A., Kubisch, A., Penteriani, V., Saastamoinen, M., Stevens, V.M., Bullock, J.M. 2013. Dispersal and species’ responses to climate change. Oikos, 122(11), 1532–1540. https://doi.org/10.1111/j.1600-0706.2013.00399.x
30.  Wang, R., Li, Q., He, S., Liu, Y., Wang, M., Huang, J. 2018. Modeling and mapping the current and future distribution of Pseudomonas syringae pv. actinidiae under climate change in China. PLOS ONE, 13(1): e0192153. https://doi.org/10.1371/journal.pone.0192153
31.  Wei, B., Wang, R., Hou, K., Wang, X., Wu, W. 2018. Predicting the current and future cultivation regions of Carthamus tinctorius L. using MaxEnt model under climate change in China. Global Ecology and Conservation, 16: e00477.
32.  Wisz, M.S., Pottier, J., Kissling, W.D., Pellissier,L., Lenoir, J., Damgaard, C.F., Dormann, C.F., Forchhammer, M.C., Grytnes, J.A., Guisan, A., Heikkinen, R.K., Høye, T.T., Kühn, I., Luoto, M., Maiorano, L., Nilsson, M.C., Normand, S., Öckinger, E., Schmidt, N. M., Svenning, J.C. 2013. The role of biotic interactions in shaping distributions and realised assemblages of species: implications for species distribution modelling. Biological Reviews, 88(1), 15–30. https://doi.org/10.1111/j.1469-185X.2012.00235.x
33.  Zhang, K., Yao, L., Meng, J., Tao, J. 2018. Maxent modeling for predicting the potential geographical distribution of two peony species under climate change. Science of the Total Environment, 634: 1326-1334.
34.  امیری، ح. 1387. شناسایی مواد تشکیل دهنده روغن اسانسی گیاه Teucrium orentale L. subsp. talori (Boiss.) Rech.f.. فصلنامه گیاهان دارویی، 7(4): 1-4.
35.  حیدری، ر.، خداقلی، م.، زارعان، م. 1394. پهنه‌بندی زیست اقلیم مریم نخودی شرقی (Teucrium orentale) در استان اصفهان با استفاده از روش‌های آماری چند متغیره و سیستم اطلاعات جغرافیایی (GIS). دوماهانامه علمی-پژوهشی تحقیقات گیاهان دارویی و معطر ایران، 31(4): 637-650.
36.  خداحامی، ق.، رضایی، م. 1393. پهنه‌بندی گیاه دارویی Teucrium polium L. در رویشگاه‌های مختلف حوزه آبخیز مهارلو استان فارس. فصلنامه اکوفیتوشیمی گیاهان دارویی 8(4): 27-35.
37.  زهیری، م.، صالحی سلمی، م.ر.، نگارش، ک.، الماسیه، ک. 1404. مدلسازی مناطق مطلوب رویشگاهی گونه Euphorbia hypericifolia در استان خوزستان. مطالعات علوم محیط زیست، 10(2): 10147-10156. https://doi.org/10.22034/jess.2025.470407.2283
38.  فخیمی، ا.، خداقلی، م. 1402. اثر تغییر اقلیم بر گستره رویشی گونه Bromus tomentellus Boiss در رویشگاه‌های مرتعی زاگرس مرکزی، استان چهارمحال و بختیاری. مطالعات علوم محیط زیست، 8(2): 6730-6740. https://doi.org/10.22034/jess.2023.380017.1949
39.  قلیچ نیا، ح. 1403. اثر تغییر اقلیم بر رویشگاه گونه Artemisia Aucheri Boisss. بر پایه مدل پیش‌بینی اقلیم در رویشگاه‌های مرتعی البرز، استان مازندران. مطالعات علوم محیط زیست. 9(4): 9719-9728. https://doi.org/10.22034/jess.2024.434233.2206
40.  قهساره اردستانی، ا.، یادگاری، ح.، اسدی، ا.، بهمنی، ب. 1404. تاثیر تغییر اقلیم بر زویشگاه‌های مطلوب زرین گیاه (Dracocephalum kotschyi) در ایران مرکزی. فصلنامه علوم محیطی نوین، 23(2): 477-490.
41.    مظفریان، و. (1396). فلور چهارمحال و بختیاری. مرکز تحقیقات و آموزش کشاورزی و منابع طبیعی استان چهار محال و بختیاری، شهرکرد، 895 ص.
42.  معتمدی، ج. خداقلی، م. خلیفه زاده، ر. 1401. پیش‌‌بینی گستره کنونی و گستره بالقوه آینده گونه Artemisia aucheri تحت دو مدل هشدار اقلیمی Rcp4.5 و Rcp8.5 در رویشگاه‌‌های مرتعی البرز جنوبی، استان قزوین. مطالعات علوم محیط زیست، 7(2): 5015-5023. https://doi.org/10.22034/jess.2022.331331.1732