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

39223258
71208
10.1038/s41598-024-71208-6
Article
Mapping climate suitability index for rainfed cultivation of medicinal plants by developing an AI-based probabilistic framework
http://orcid.org/0000-0002-9018-547X
Sadeghfam Sina s.sadeghfam@maragheh.ac.ir

1
Rahmani Mohammad Sina 1
Moazamnia Marjan 2
Morshedloo Mohammad Reza 3
1 https://ror.org/0037djy87 grid.449862.5 0000 0004 0518 4224 Department of Civil Engineering, Faculty of Engineering, University of Maragheh, P.O. Box 55136-553, Maragheh, East Azerbaijan Iran
2 https://ror.org/01app8660 grid.440821.b 0000 0004 0550 753X Department of Civil Engineering, Faculty of Engineering, University of Bonab, Bonab, East Azerbaijan Iran
3 https://ror.org/0037djy87 grid.449862.5 0000 0004 0518 4224 Department of Horticultural Science, Faculty of Agriculture, University of Maragheh, Maragheh, East Azerbaijan P.O. Box 55136-553, Iran
2 9 2024
2 9 2024
2024
14 2041314 5 2024
26 8 2024
© The Author(s) 2024
2024
https://creativecommons.org/licenses/by/4.0/ Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, 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 changes were made. 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/4.0/.
The Climate Suitability Index (CSI) can increase agricultural efficiency by identifying the high-potential areas for cultivation from the climate perspective. The present study develops a probabilistic framework to calculate CSI for rainfed cultivation of 12 medicinal plants from the climate perspective of precipitation and temperature. Unlike the ongoing frameworks based on expert judgments, this formulation decreases the inherent subjectivity by using two components: frequency analysis and Particle Swarm Optimization (PSO). In the first component, the precipitation and temperature layers were prepared by calculating the occurrence probability for each plant, and the obtained probabilities were spatially interpolated using geographical information system processes. In the second component, PSO quantifies CSI by classifying a study area into clusters using an unsupervised clustering technique. The formulation was implemented in the Lake Urmia basin, which was distressed by unsustainable water resources management. By identifying clusters with higher CSI values for each plant, the results provide deeper insights to optimize cultivation patterns in the basin. These insights can help managers and farmers increase yields, reduce costs, and improve profitability.

Keywords

Agro-climate
Frequency analysis
Lake Urmia
Particle swarm optimization
Rainfed agriculture
Subject terms

Agroecology
Hydrology
issue-copyright-statement© Springer Nature Limited 2024
==== Body
pmcIntroduction

Extending agricultural species with medicinal and commercial values can provide new opportunities for local communities to increase efficiency in agriculture1. The limitation in water and soil resources, pressure on natural ecosystems, climate change, drought, land use change, and increasing cultivation on slopes highlight the necessity for paying attention to sustainable agriculture2. Various strategies and approaches, including low-input, ecological, and conservation agriculture, have been developed to achieve sustainable agricultural systems3. Low-input agriculture aims to identify areas suitable for low-water-demand medicinal plants, replacing high-water-consumption crops.

Identifying areas with appropriate climate requirements and ecological environments is beneficial for cultivating high-quality medicinal plants and increasing the efficiency of agriculture. Soil, precipitation, and temperature are variables that affect the growth and quality of medicinal plants. While all these variables are important, soil characteristics change insignificantly over time, making external climate factors like precipitation and temperature more influential on the long-term growth dynamics of medicinal plants4. Different studies were undertaken to investigate the Climate Suitability Index (CSI), which reflects the the adaptivity degree of the climatic conditions in a study area for a specific type of crop (for example, see5). The literature review emphasizes indices known as CSI, which play a crucial role in shaping biodiversity. These indices have various applications in various fields, such as agriculture. However, the results are on a large scale and provide a general overview of agriculture without considering specific types of cultivation.

The common components for calculating cultivation potential are the incorporated data layers and overlaying techniques, summarized in Table 1 for a set of previous studies together with the plant name and study area. The table shows that different data layers are incorporated, which can be categorized into climatic variables (e.g., precipitation, temperature, and humidity), soil quality (e.g., pH, EC), and basin physiography (e.g., slope and altitude). The table also compared the studies from the point of view of overlaying techniques. The critical issue in this part is the relative importance of the incorporated data layers, referred to as the weights. The Analytical Hierarchical Process (AHP) is the most popular overlaying technique, which calculates the weights based on expert judgment using a pairwise comparison matrix. Other incorporated techniques are also relying on expert judgment. Notably, subjectivity can lead to inconsistencies in the outcomes, making the reproducibility of results difficult. Including subjective elements in modeling can also make it challenging to interpret the results, and it may not be transparent due to the influence of personal biases. Table 1 List of previous studies from different point of views.

Row	Data layers	Overlaying techniques	Plant	Study area	References	
1	Average/Min/Max annual temperature, Altitude, Soil covariates	Boolean and fuzzy logic	Rapeseed	Turkey	6	
2	Precipitation, temperature, slope, altitude, soil texture, slope direction	AHP	Astragalus hypsogeton bunge	Iran	6	
3	Land surface temperature, bareness of land, slope, surface elevation, humidity and rainfall	Equal weight overlap analysis	Jatropha curcas	Pakistan	7	
4	Temperature, precipitation, humidity, slope, altitude, evaporation, slope direction, pH	AHP	Saffron	Iran	8	
5	Soil parameters including, CaCO3, EC, PH, Slope, Texture	AHP	Rainfed wheat	Turkey	9	
6	Average temperature in the flowering period, average temperature in the winter season, average annual relative humidity, average annual rainfall, average daily sunny time, altitude, slope, land use capability	AHP	Citrus	Turkey	5	
7	Annual average temperature and precipitation, slope, slope direction, altitude, land use, soil texture	Weighted linear combination	Rainfed Chickpea and Irrigated Chickpea	Iran	10	
8	Average isothermal temperature, seasonal temperature,

annual rainfall and seasonal rainfall, evaporation

	Maximum entropy model (Maxent)	Amorpha canescens,

Leucaena collinsi,

Cynodon dactylon,

Atriplex canescens (Pursh)

	Mexico	11	
9	Temperature, precipitation	Maxent	Ephedra sinica	China	12	
10	19 climatic and edaphic factors	Maxent	Cuphea aequipetala Cav	Mexico	13	
11	Four bioactive compounds (madecassoside, asiaticoside, madecassic acid, and asiatic acid)	–	Centella asiatica	India	14	
12	Temperature and precipitation	Entropy Catastrophe Scheme	Wheat	Iran	15	

Some overlaying techniques, such as the maximum entropy model, decrease the inherent subjectivity by modeling the geographic distribution of species with presence records and determining the distribution through the maximum entropy probability11. Also, the ensemble ecological niche model predicts the ecological position of a species by determining the most vital environmental elements14 and the weighted linear combination, which assigns a weight to each layer based on its importance5.

Unsupervised learning algorithms for cluster analysis employ unstructured data that are grouped based on similarities and patterns. Examples of these algorithms are the genetic algorithm16, density-based spatial clustering of applications with noise, and Particle Swarm Optimization (PSO)17. PSO was incorporated in previous studies as an optimizer model, e.g., mapping and protecting agricultural land under spatial constraints18 and determining soil nutrient management areas19. This study used PSO to combine precipitation and temperature layers to decrease the inherent subjectivity. PSO offers several advantages over other clustering techniques, including the ability to efficiently handle high-dimensional data, simplicity of implementation, and quick convergence to a solution20. Additionally, PSO is less likely to be trapped in local optima compared to other optimization techniques21. Notably, the study does not aim to evaluate the clustering or optimization techniques; further techniques are available in the literature22–24.

In addition to the economic benefits, the rainfed cultivation of medicinal plants can take a significant step toward restoring the distressed Lake Urmia. This lake, located in the northwest of Iran, experienced a state of distress regarding water resource unsustainability25,26, and it lost 96% of its surface and 91% of its volume between 1996 and 2016. One of the main factors that contributed to the drying up of Lake Urmia is the overexploitation of water resources due to the unsustainable development of agriculture, which the evidence is a 400% increase in the number of wells in the basin between 1984 and 2012 and a 60% increase in croplands development between 1987 and 202027. The agricultural sector has played a vital role, accounting for almost 90% of the total water consumption in the basin28,29. Therefore, the indiscriminate expansion of agriculture can cause serious environmental impacts on the lake.

The paper formulates a framework based on probability and Artificial Intelligence (AI) modeling to delineate CSI for 12 medicinal plants within Lake Urmia. Most of the previous studies are based on expert judgments, and the present framework decreases the inherent subjectivities, as outlined as follows: (i) conducting frequency analysis to calculate the occurrence probability of desired precipitation and temperature, unlike the previous studies, which deterministically use subjective rate or expected values in the geospatial analysis; and (ii) overlaying the precipitation and temperature layers using the PSO algorithm as an unsupervised clustering technique to decrease the inherent subjectivity related to scoring systems. The results offer different cultivation options for managers and local farmers, giving them the opportunity to switch from their current cultivation practices. On the other hand, the unsustainable management of water resources in the agricultural sector has caused a crisis in the Lake Urmia basin, making it critical to apply the developed formulation in the area.

Study area

The Lake Urmia basin, located in the northwest of Iran, has an approximate area of 51,876 km2, which occupies 3.2% of the total area of Iran. The basin is surrounded by mountains such as the southern slopes of Mount Sablan and the northern part of the Zagros Mountains, as well as the north, western, and southern slopes of Mount Sahand. Figure 1a shows the geographical location and altitude of the basin, varying from 1171 to 3700 m, considering the Digital Elevation Model (DEM) (SRTM-1 satellite in 2014). Figure 1b shows the land use map of the basin using Sentinel-2 satellite images from 2020 to 2022. The figure indicates that the study area was occupied by water bodies (1.8%), cultivated lands (28.7%), build-up lands (3.4%), salty lands (6.8%), and barren lands (59.3%).Fig. 1 Study area: (a) geological location; (b) land use; (c) annual average precipitation in 2000–2020; (d) daily average temperature in 2000–2020. Note The figure is produced by the authors using QGIS 3.10.

Figure 1c,d show the spatial distributions for the annual precipitation in 57 stations and daily temperature in 11 stations in 2000–2020 obtained by the IDW interpolation technique. Notably, the data were provided by the Iran Meteorological Organization. The precipitation in Fig. 1c varies from 192 to 890 mm, with an average of 390 mm in 2000–2020. The figure indicates that precipitation decreases from the southwest parts to the north of the basin. The daily average temperature in Fig. 1d varies from 9 to 15 °C, with an average value of 12 °C. Based on average precipitation and temperature values, the climate of the basin is cold and semi-arid, as per Emberger (1930)30.

Results

Precipitation distributions

Table S1 presents the goodness-of-fit results for the selected PDFs based on the Kolmogorov–Smirnov (K–S) test. Figure 2 depicts the spatial distribution of the probability of occurrence for desired annual precipitation, interpolated by IDW and classified into five bands. Notably, the occurrence probability values were calculated based on the frequency analyses by the fitted PDFs given in Table S1, considering the desired precipitation in Table 2. According to the figure, Achillea millefmlant can be cultivated throughout the basin. Also, the cultivation of Allium ascalonicom, Matricaria chamomilla, Rosa damasena, and Silybum marianum in the north and northeast of the basin has an average potential, but the potential increases towards the south. For Capparis spinosa and Cuminum cyminum, the south of the basin is located in areas with the highest potential, and the north is located in relatively high areas. Echinacea purpurea, Carthamus tinctorius, Lavandula angustifolia, and Thymus daenensis have the highest potential in the southern parts, but the potential decreases towards the north of the basin. For Satureja hortensis, the high potential part was only observed in the southwest parts of the basin.Fig. 2 Probability occurrence values of desired precipitation. Note The figure is produced by the authors using QGIS 3.1.

Table 2 Investigated medicinal plants and desired precipitation.

Plant	Desired annual precipitation (mm)	Temperature	
Desired value (°C)	Growth period	
1	Achillea millefolium	120	18 to 26	Early April to late August	
2	Allium ascalonicom	250	9 to 13	Late October to early July	
3	Capparis spinose	200	− 8 to 50	Early March to early October	
4	Carthamus tinctorius	300	15 to 20	Late March to late July	
5	Cuminum cyminum	150	15 to 21	Late March to late July	
6	Echinacea purpurea	300	18 to 26	Early March to early October	
7	Lavandula angustifolia	300	18 to 25	Mid-April to late September	
8	Matericaria chamomilla	140	15 to 19	Mid-April to late August	
9	Rosa damasena	250	15 to 22	Early April to mid-July	
10	Satureja hortensis	450	20 to 25	Mid-April to mid-September	
11	Silybum marianum	240	18 to 20	Late March to mid-August	
12	Thymus daenensis	300	18 to 21	Early May to late September	

Temperature distributions

Figure 3 shows the spatial distribution of the occurrence probability of the desired temperature in the related growth period of each plant, which was interpolated by IDW and classified into five bands. Notably, a set of frequency analyses were conducted to calculate the occurrence probability using the fitted PDFs presented in Table S2, considering the desired temperature in Table 2. According to the figure, Capparis spinosa can be cultivated throughout the basin. Also, the cultivation of Satureja hortensis, Matricaria chamomilla, Echinacea purpurea, Achillea millefmlant, and Silybum marianum in the east of the basin has medium and weak potential. However, the potential increases towards the central and southwestern regions. Carthamus tinctorius and Lavandula angustifolia have the highest potential in the west and central regions, but the potential is low in the northeast. Thymus daenensis has the highest potential in the northeastern and southern areas, but the potential in the rest of the parts is medium and weak. For Rosa damasena, the central and western regions have the highest potential, and the potential decreases towards the northern and northeastern regions. The potential for Cuminum cyminum is quite high in the north and northwest and moderate from the central to the southeast. Allium ascalonicom has high potential in some parts of the west and central areas, but the potential is low in other areas.Fig. 3 Probability occurrence values of desired temperature. Note The figure is produced by the authors using QGIS 3.1.

CSI distributions by PSO

CSI values for the investigated medicinal plant were calculated using PSO to decrease the inherent subjectivity within expert judgment with the weights of precipitation and temperature layers. As an AI-based unsupervised clustering technique, PSO conducts cluster analysis to classify the study area into five bands by minimizing the objective function described in the method section. The PSO configuration includes a maximum iteration to complete the optimization process equal to 200, and the number of clusters is five. The constriction coefficients method was used to calculate the coefficients of inertia, in which two positive numbers were considered with a sum greater than four. Figure 4 depicts the objective function diagrams for each plant, which all converged, considering their tendency to the x-axis by increasing the iteration. After conducting cluster analysis, the precipitation and temperature layers were classified into five clusters using the PSO algorithm, and CSI clusters were distributed in Fig. 5 for each plant. Notably, CSI increases from Cluster 1 to 5, corresponding to the red-to-green color palette.Fig. 4 Objective function diagrams of PSO for clustering precipitation and temperature.

Fig. 5 CSI for the investigated medicinal plants. Note The figure is produced by the authors using QGIS 3.1.

Figure 5 indicates that the CSI for Achillea millefolium and Matricaria chamomilla is low in the eastern regions and becomes high in the central and western regions. The CSI for Allium ascalonicom is low in the northern parts but increases from the south to the central parts. The CSI for Capparis spinosa is low in the north half of the basin, but it increases from the central parts to the south of the basin. Carthamus tinctorius has the lowest CSI in the areas around Lake Urmia and the eastern parts, but the higher CSI clusters exist in the south and southwest parts. The CSI for Cuminum cyminum is high in northern and western regions but gradually decreases towards the central and southeastern regions. Echinacea purpurea can be cultivated in the central and southern regions, but cultivation is limited to the northern parts of the basin. The CSI for Lavandula angustifolia is low in the northeastern regions but moderate in the northern regions and increases towards the southern parts. The CSI for Rosa damasena varies from the low clusters in the northeastern regions to the central and western regions. Satureja hortensis has lower CSI in the northeastern and southeastern regions, and the CSI in northern regions changes from low and moderate to high clusters towards the south. The CSI for Silybum marianum is unsuitable in the central and northeastern areas, but it gradually increases from the southeast to the northwest. Thymus Daenensis has a lower CSI cluster in the central and northeastern regions; it increases in the northwest and southern regions.

Discussion

The potential of rainfed cultivation for 12 medicinal plants in the Lake Urmia basin can be one of the promising solutions for the ecological and economic future of the basin and the lake distressed by unsustainable use of water resources. Previous studies investigated the relationship between the decrease in the Lake Urmia level and the agriculture growth27. If policymakers and managers adopt the agricultural patterns based on the present and similar studies dealing with rainfed agriculture, the restoration of Lake Urmia is possible since the agriculture sector consumes almost 90% of water. Even with this, the developed framework is not confined to specific case studies and can be adjusted to various basins and regions. Its application in other case studies is imperative to address its vulnerabilities and enhance its efficacy.

The study employed precipitation and temperature historical records, but using projected climate variables under climate change scenarios can provide a deeper insight into the climate suitability index. Utilizing projected climate variables allows policymakers to explore adaptive strategies to mitigate the impacts of climate change. While previous studies31 have investigated these impacts on indices concerning cultivation potential, the innovation in the suggested formwork sets it apart from the goals of the current study.

The investigated medicinal plants are compatible with the climate condition of the basin and, additionally, can generate higher economic benefits for local farmers. However, climate suitability is just one aspect of achieving the goals of this study. It is important to address other perspectives, such as the social impacts of the local farmers’ change in cultivation type. The change in cultivation type is closely tied to the livelihoods of local farmers, so it is crucial to ensure more comprehensive knowledge. Certain plants, such as Rosa damasena, may pose challenges to the local farmers, and they are difficult to cultivate due to their blades. Moreover, changing the cultivation type practiced for many years may not bring enough economic benefits in the initial stages of cultivation. Therefore, governments need to provide financial support for such programs.

Pouladi et al.32 discussed some socio-hydrological barriers to the restoration of Lake Urmia. They declared that farmers experiencing financial problems were unsympathetic to the drying up or restoration of Lake Urmia. Changing the cultivation pattern in low economic conditions towards medicinal plants requires financial support for farmers because it will take time to benefit the farmer. The successful implementation of such decisions requires the participation of farmers, which can be achieved by understanding and reflecting their concerns in the decision-making process.

In addition to social aspects, economic and environmental aspects play significant roles in developing a comprehensive index for evaluating plant cultivation. When assessing plants for suitability, economic considerations such as market demand and value can be crucial in determining which species are prioritized for further investigation. By understanding the economic significance of these plants, policymakers can make informed decisions about resource allocation and conservation efforts. From an environmental standpoint, it is essential to recognize the potential impact of introduced plants on native species and ecosystems. Some plant species may exhibit invasive behavior, outcompeting native flora and disrupting the ecological balance. Therefore, careful evaluation of a plant's potential invasiveness is essential to safeguarding the biodiversity and stability of an area. Another environmental aspect to consider in plant evaluation is the quality and texture of the soil. Soil characteristics can significantly influence plants' growth and survival, affecting their overall suitability for a particular environment. By considering soil conditions, researchers can better understand the constraints and opportunities for plant species in a given habitat, leading to a more thorough and holistic assessment.

The cultivation of medicinal plants in the Lake Urmia basin is limited, and even sparse data is not available, making it difficult to verify the results directly. However, the goodness of fit of the fitted Probability Density Functions (PDFs) based on the K–S test confirms the accuracy of the results. It should be noted that errors due to interpolation may affect the results, but we have tried to manage this by making maximum use of the available precipitation stations. The number of precipitation stations is relatively adequate and covers the study area. Although the number of temperature stations is relatively low, gradual changes in temperature compared to precipitation make it acceptable.

The CSI is affected by climate variables obtained by data from point stations, and the IDW technique spatially distributes the point results. This procedure neglects microclimates between the stations, which arise from elevation, urban development, and proximity to bodies of water. This is regarded as the other limitation of present studies. The main sources of microclimate in the basin are Sahand Mountain in the western part of the basin and the water crisis in Lake Urmia. Previously, the large size of Lake Urmia played an influential role in regulating microclimates. However, as it continues to shrink, evidence shows that the daily temperature range has widened, humidity levels have decreased, and the annual precipitation pattern has shifted33,34. These variations provide room for future studies.

Conclusion

The study developed a probabilistic framework for identifying the high-potential areas for rainfed cultivation of 12 medicinal plants based on precipitation and temperature. The study delineated the high-potential areas by calculating an index, referring to SCI. Two main components for SCI calculations include (i) climate variables, prepared by frequency analysis to calculate the occurrence probability of desired precipitation and temperature, and (ii) the combination of climate variables to calculate CSI, conducted by PSO as an unsupervised clustering technique. These components reduce the inherent subjectivity associated with expert judgment. We implemented this framework in the Lake Urmia basin, facing challenges from unsustainable water resources management. The study identified suitable areas for rainfed cultivation of each plant within the basin, which could help policymakers make informed decisions to enhance agriculture efficiency and preserve water resources. This breakthrough has the potential to transform the region's agriculture sector and set an example for other areas facing similar challenges. The limitations of the study provide fertile ground for future studies, including microclimate changes due to Lake Urmia dryness, the interpolation error due to microclimate changes between stations, the impact of climate change on CSI, and the incorporation of social, economic, and environmental aspects of the investigated plants within the framework.

Method

This study formulates a probabilistic framework for mapping CSI of 12 medicinal plants using climate variables of precipitation and temperature. Figure 6 shows the methodological flowchart of the framework, and the essential components outlined as follows: (i) data pre-processing; (ii) fitting appropriate Probability Distribution Functions (PDFs) to precipitation and temperature data at each station considering the desired precipitation and temperature for each plant (see Table 2); (iii) and conducting the frequency analysis for calculating the occurrence probability for the desired precipitation and temperature; (iv) spatially distributing the calculated occurrence probabilities of precipitation and temperature for each plant; and (iv) clustering precipitation and temperature layers for each plant by the PSO algorithm.Fig. 6 Methodological flowchart of the probabilistic framework for mapping CSI.

Pre-processing

The outlier test and the gap data estimation were conducted in the pre-processing stage. Outlier data were detected using the interquartile range method35. Also, gap data were estimated by fitting the regression equations between stations with gap data and the station having the highest correlation coefficient.

Precipitation layers

Precipitation in rainfed and irrigated cultivations plays a vital role in increasing the efficiency of agriculture. The medicinal plants were selected based on their adaptability in the cold and semi-dry climate of the Lake Urmia basin, as well as their optimal water requirement during the growth period. Table 2 presents the investigated plants and corresponding desired annual precipitation values and desired temperature corresponding to the growth periods, which are collected referring to propagation books published by the Ministry of Agriculture Jihad.

Temperature layers

Temperature is crucial for calculating the CSI of medicinal plants, and desired values are required during the growth period. Table 2 presents the growth period and the corresponding desired temperature refereeing propagation books published by the Ministry of Agriculture Jihad. Notably, the Capparis Spinosa has deep and extensive roots and can tolerate severe summers in areas with low precipitation and temperatures of 50 °C. This plant not only shows considerable resistance to water shortage and high temperatures, but it is also resistant to cold and can survive up to − 8 °C36.

Goodness-of-fit for selecting PDFs

The precipitation and temperature data were extracted at each station, and a broad set of PDFs were fitted to the extracted data. The goodness-of-fit test for the fitted PDFs was conducted in the study by the K–S test37. Then, the appropriate PDFs were accepted in this test with a significance level of α = 0.05.

Spatial interpolation

The study employed the Inverse Distance Weighted (IDW) technique to spatially interpolate the probability occurrence of desired precipitation and temperature at different stations. IDW uses the neighboring points by averaging the sample points located in the neighborhood of each unknown point. Also, IDW assumes that the effect of each point is proportional to the power of the inverse of its distance, so the effect of the known point decreases with increasing distance38.

Artificial intelligence model

As a AI-based algorithm, PSO was first introduced by Kennedy et al.39. PSO has the capability to optimize non-linear problems with immediate convergence rate and a relatively low calculation burden. The design of this algorithm is inspired by the movement of flocks and fish to optimize the path to food40. In PSO, each community member is called a particle, and the position of each particle is a potential solution to the optimization problem. The position of the particles changes according to their experience that of their neighbors during the execution of the algorithm, and thus, the particles float and move in the search space. As the particles move, each stores a memory of its best value and corresponding position along its path. The value of the particles is measured according to the criterion function related to the optimization problem. For the ith particle, the best value and its corresponding position are stored in the variables pbesti and pi=(pi1.pi2.⋯.piD), respectively. Notably, D shows the dimension of the search space. Also, the best particle position in all iterations and its corresponding position are stored in the same variable. The collection of n particles in a community is called population and is expressed as pop=(p1.p2.⋯.pn)37. The basis of the PSO method is based on the principles accelerating towards their pbesti and gbesti in each iteration, with the addition of the displacement vector in each particle, a new position in the search space for each particle is found. At the end of the algorithm, the position corresponding to gbesti is considered as the optimal solution of the problem. In each iteration, the position and velocity of the ith particle are updated according to Eqs. (1) and (2):1 viiter+1=wviditer+c1randOpid-xiditer+c2randOpgd-xiditer

2 xiditer+1=xiditer+Viditer+1

where Vi=(vi1.vi2.⋯.vid) is the velocity vector of the ith particle; and Xi=(xi1.xi2.⋯.xiD) is its position vector. In this regard, c1 and c2 are positive numbers called personal and global learning coefficient, respectively. The randΟ is a function that generates random numbers between zero and one. w is the inertia factor that controls the effect of the current speed on the current speed, and its allowed value in the standard PSO algorithm changes in the range of 0.4 to 0.9.

The study transfers the clustering problem into the optimization problem by defining the objective function and connecting it to the PSO algorithm. The objective function is defined by minimizing the distances of each data from its cluster center, i.e., within-cluster distance by Eq. (3):3 ObjFunction=min∑j=1k∑i=1ndxi.mj

where xi is a member of the total data of X=x1.x2.⋯.xn; mj is a member of the cluster centers as X=x1.x2.⋯.xn; and the d represents the Euclidian distance.

Supplementary Information

Supplementary Tables.

Supplementary Information

The online version contains supplementary material available at 10.1038/s41598-024-71208-6.

Author contributions

All authors contributed to the formulation of the framework. Data collection, literature review and preliminary analysis were performed by S.S., S.R. The study was supervised and drafted by S.S., M.M. and M.R.M. The vision and review of the drafted manuscript at all of its stages and interpretation of results were done by S.S. All authors read and approved the final manuscript.

Data availability

Data are available from the corresponding author upon 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. Aghaloo K Sharifi A A GIS-based agroecological model for sustainable agricultural production in arid and semi-arid areas: The case of Kerman Province, Iran Curr. Res. Environ. Sustain. 2023 6 100230 10.1016/j.crsust.2023.100230
Aghaloo, K. & Sharifi, A. A GIS-based agroecological model for sustainable agricultural production in arid and semi-arid areas: The case of Kerman Province, Iran. Curr. Res. Environ. Sustain. 6, 100230 (2023).10.1016/j.crsust.2023.100230
2. Rocha J Quintela A Serpa D Keizer JJ Fabres S Water yield and biomass production for on a eucalypt-dominated Mediterranean catchment under different climate scenarios J. For. Res. 2023 34 5 1263 1278 10.1007/s11676-022-01590-2
Rocha, J., Quintela, A., Serpa, D., Keizer, J. J. & Fabres, S. Water yield and biomass production for on a eucalypt-dominated Mediterranean catchment under different climate scenarios. J. For. Res. 34(5), 1263–1278 (2023).10.1007/s11676-022-01590-2
3. Amiri Z Asgharipour MR Campbell DE Azizi K Kakolvand E Moghadam EH Conservation agriculture, a selective model based on emergy analysis for sustainable production of shallot as a medicinal-industrial plant J. Clean. Prod. 2021 292 126000 10.1016/j.jclepro.2021.126000
Amiri, Z. et al. Conservation agriculture, a selective model based on emergy analysis for sustainable production of shallot as a medicinal-industrial plant. J. Clean. Prod. 292, 126000 (2021).10.1016/j.jclepro.2021.126000
4. Tang X Yuan Y Wang L Chen S Liu X Zhang J Identifying prioritized planting areas for medicinal plant Thesium chinense Turcz. under climate change in China Ecol. Inform. 2021 66 101459 10.1016/j.ecoinf.2021.101459
Tang, X. et al. Identifying prioritized planting areas for medicinal plant Thesium chinense Turcz. under climate change in China. Ecol. Inform. 66, 101459 (2021).10.1016/j.ecoinf.2021.101459
5. Nabati J Nezami A Neamatollahi E Akbari M GIS-based agro-ecological zoning for crop suitability using fuzzy inference system in semi-arid regions Ecol. Ind. 2020 117 106646 10.1016/j.ecolind.2020.106646
Nabati, J., Nezami, A., Neamatollahi, E. & Akbari, M. GIS-based agro-ecological zoning for crop suitability using fuzzy inference system in semi-arid regions. Ecol. Ind. 117, 106646 (2020).10.1016/j.ecolind.2020.106646
6. Shokati B Feizizadeh B Sensitivity and uncertainty analysis of agro-ecological modeling for saffron plant cultivation using GIS spatial decision-making methods J. Environ. Plan. Manag. 2019 62 3 517 533 10.1080/09640568.2018.1427561
Shokati, B. & Feizizadeh, B. Sensitivity and uncertainty analysis of agro-ecological modeling for saffron plant cultivation using GIS spatial decision-making methods. J. Environ. Plan. Manag. 62(3), 517–533 (2019).10.1080/09640568.2018.1427561
7. Dedeoğlu M Dengiz O Generating of land suitability index for wheat with hybrid system aproach using AHP and GIS Comput. Electron. Agric. 2019 167 105062 10.1016/j.compag.2019.105062
Dedeoğlu, M. & Dengiz, O. Generating of land suitability index for wheat with hybrid system aproach using AHP and GIS. Comput. Electron. Agric. 167, 105062 (2019).10.1016/j.compag.2019.105062
8. Khan MJ Hasan SA GiS-based screening model of coastal city Karachi for plantation of biofuel source Sci. Rep. 2020 10 1 4666 10.1038/s41598-020-61052-9 32170096
Khan, M. J. & Hasan, S. A. GiS-based screening model of coastal city Karachi for plantation of biofuel source. Sci. Rep. 10(1), 4666 (2020).32170096 10.1038/s41598-020-61052-9
9. Tercan E Dereli MA Development of a land suitability model for citrus cultivation using GIS and multi-criteria assessment techniques in Antalya province of Turkey Ecol. Indic. 2020 117 106549 10.1016/j.ecolind.2020.106549
Tercan, E. & Dereli, M. A. Development of a land suitability model for citrus cultivation using GIS and multi-criteria assessment techniques in Antalya province of Turkey. Ecol. Indic. 117, 106549 (2020).10.1016/j.ecolind.2020.106549
10. Buendía-Espinoza JC Martínez-Ochoa EDC Díaz-Aguilar I Cahuich-Damián JE Zamora-Elizalde MC Identifying potential planting sites for three non-native plants to be used for soil rehabilitation in the Tula Watershed Forests 2022 13 2 270 10.3390/f13020270
Buendía-Espinoza, J. C., Martínez-Ochoa, E. D. C., Díaz-Aguilar, I., Cahuich-Damián, J. E. & Zamora-Elizalde, M. C. Identifying potential planting sites for three non-native plants to be used for soil rehabilitation in the Tula Watershed. Forests 13(2), 270 (2022).10.3390/f13020270
11. Zhang K Liu Z Abdukeyum N Ling Y Potential geographical distribution of medicinal plant ephedra sinica stapf under climate change Forests 2022 13 12 2149 10.3390/f13122149
Zhang, K., Liu, Z., Abdukeyum, N. & Ling, Y. Potential geographical distribution of medicinal plant ephedra sinica stapf under climate change. Forests 13(12), 2149 (2022).10.3390/f13122149
12. Garibay-Castro LR Gutiérrez-Yurrita PJ López-Laredo AR Hernández-Ruíz J Trejo-Espino JL Potential distribution and medicinal uses of the Mexican plant Cuphea aequipetala Cav. (Lythraceae) Diversity 2022 14 5 403 10.3390/d14050403
Garibay-Castro, L. R., Gutiérrez-Yurrita, P. J., López-Laredo, A. R., Hernández-Ruíz, J. & Trejo-Espino, J. L. Potential distribution and medicinal uses of the Mexican plant Cuphea aequipetala Cav. (Lythraceae). Diversity 14(5), 403 (2022).10.3390/d14050403
13. Tri CL Tuan TA Kiem TN Application of AHP-GIS model to assess the ecological suitability of Codonopsis javanica in Kon Plong District Kon Tum Province, Viet Nam J. Ecol. Eng. 2022 10.12911/22998993/150027
Tri, C. L., Tuan, T. A. & Kiem, T. N. Application of AHP-GIS model to assess the ecological suitability of Codonopsis javanica in Kon Plong District Kon Tum Province, Viet Nam. J. Ecol. Eng.10.12911/22998993/150027 (2022).10.12911/22998993/150027
14. Shaban M Ghehsareh Ardestani E Ebrahimi A Borhani M Climate change impacts on optimal habitat of Stachys inflata medicinal plant in central Iran Sci. Rep. 2023 13 1 6580 10.1038/s41598-023-33660-8 37085511
Shaban, M., Ghehsareh Ardestani, E., Ebrahimi, A. & Borhani, M. Climate change impacts on optimal habitat of Stachys inflata medicinal plant in central Iran. Sci. Rep. 13(1), 6580 (2023).37085511 10.1038/s41598-023-33660-8
15. Sadeghfam S Saghayeshi H Moazamnia M Mapping cultivation potential index of rainfed wheat from climatic perspective using Shannon entropy catastrophe scheme Comput. Electron. Agric. 2024 217 108616 10.1016/j.compag.2024.108616
Sadeghfam, S., Saghayeshi, H. & Moazamnia, M. Mapping cultivation potential index of rainfed wheat from climatic perspective using Shannon entropy catastrophe scheme. Comput. Electron. Agric. 217, 108616 (2024).10.1016/j.compag.2024.108616
16. Hruschka ER Ebecken NF A genetic algorithm for cluster analysis Intell. Data Anal. 2003 7 1 15 25 10.3233/IDA-2003-7103
Hruschka, E. R. & Ebecken, N. F. A genetic algorithm for cluster analysis. Intell. Data Anal. 7(1), 15–25 (2003).10.3233/IDA-2003-7103
17. Verma S Jyoti K A PSO optimized layered approach for parametric clustering on weather dataset Int. J. Mod. Eng. Res. (IJMER) 2013 3 1 504 508
Verma, S. & Jyoti, K. A PSO optimized layered approach for parametric clustering on weather dataset. Int. J. Mod. Eng. Res. (IJMER) 3(1), 504–508 (2013).
18. Zhang R Li J Du Q Ren F Basic farmland zoning and protection under spatial constraints with a particle swarm optimisation multiobjective decision model: A case study of Yicheng, China Environ. Plan. B Plan. Design 2015 42 6 1098 1123 10.1068/b130213p
Zhang, R., Li, J., Du, Q. & Ren, F. Basic farmland zoning and protection under spatial constraints with a particle swarm optimisation multiobjective decision model: A case study of Yicheng, China. Environ. Plan. B Plan. Design 42(6), 1098–1123 (2015).10.1068/b130213p
19. Fu Q Wang Z Jiang Q Delineating soil nutrient management zones based on fuzzy clustering optimized by PSO Math. Comput. Model. 2010 51 11–12 1299 1305 10.1016/j.mcm.2009.10.034
Fu, Q., Wang, Z. & Jiang, Q. Delineating soil nutrient management zones based on fuzzy clustering optimized by PSO. Math. Comput. Model. 51(11–12), 1299–1305 (2010).10.1016/j.mcm.2009.10.034
20. Yang XS Nature-Inspired Metaheuristic Algorithms 2010 Luniver Press
Yang, X. S. Nature-Inspired Metaheuristic Algorithms (Luniver Press, 2010).
21. Mendes R Kennedy J Neves J The fully informed particle swarm: Simpler, maybe better IEEE Trans. Evol. Comput. 2004 8 3 204 210 10.1109/TEVC.2004.826074
Mendes, R., Kennedy, J. & Neves, J. The fully informed particle swarm: Simpler, maybe better. IEEE Trans. Evol. Comput. 8(3), 204–210 (2004).10.1109/TEVC.2004.826074
22. Rahmati K Ashofteh PS Loáiciga HA Application of the grasshopper optimization algorithm (GOA) to the optimal operation of hydropower reservoir systems under climate change Water Resour. Manag. 2021 35 13 4325 4348 10.1007/s11269-021-02950-z
Rahmati, K., Ashofteh, P. S. & Loáiciga, H. A. Application of the grasshopper optimization algorithm (GOA) to the optimal operation of hydropower reservoir systems under climate change. Water Resour. Manag. 35(13), 4325–4348 (2021).10.1007/s11269-021-02950-z
23. Azadi F Ashofteh PS Shokri A Loáiciga HA Development of the FA-KNN hybrid algorithm and its application to reservoir operation Theor. Appl. Climatol. 2024 155 2 1261 1280 10.1007/s00704-023-04688-7
Azadi, F., Ashofteh, P. S., Shokri, A. & Loáiciga, H. A. Development of the FA-KNN hybrid algorithm and its application to reservoir operation. Theor. Appl. Climatol. 155(2), 1261–1280 (2024).10.1007/s00704-023-04688-7
24. Moghadam SH Ashofteh PS Loáiciga HA Optimal water allocation of surface and ground water resources under climate change with WEAP and IWOA modeling Water Resour. Manag. 2022 36 9 3181 3205 10.1007/s11269-022-03195-0
Moghadam, S. H., Ashofteh, P. S. & Loáiciga, H. A. Optimal water allocation of surface and ground water resources under climate change with WEAP and IWOA modeling. Water Resour. Manag. 36(9), 3181–3205 (2022).10.1007/s11269-022-03195-0
25. Sadeghfam S Mirahmadi R Khatibi R Mirabbasi R Nadiri AA Investigating meteorological/groundwater droughts by copula to study anthropogenic impacts Sci. Rep. 2022 12 8285 10.1038/s41598-022-11768-7 35585219
Sadeghfam, S., Mirahmadi, R., Khatibi, R., Mirabbasi, R. & Nadiri, A. A. Investigating meteorological/groundwater droughts by copula to study anthropogenic impacts. Sci. Rep. 12, 8285 (2022).35585219 10.1038/s41598-022-11768-7
26. Hooshyaripor F Sardari J Dehghani M Noori R A new concept of drought feeling against the meteorological drought Sci. Rep. 2022 12 1 16711 10.1038/s41598-022-21181-9 36202951
Hooshyaripor, F., Sardari, J., Dehghani, M. & Noori, R. A new concept of drought feeling against the meteorological drought. Sci. Rep. 12(1), 16711 (2022).36202951 10.1038/s41598-022-21181-9
27. Barideh R Nasimi F Investigating the changes in agricultural land use and actual evapotranspiration of the Urmia Lake basin based on FAO’s WaPOR database Agric. Water Manag. 2022 264 107509 10.1016/j.agwat.2022.107509
Barideh, R. & Nasimi, F. Investigating the changes in agricultural land use and actual evapotranspiration of the Urmia Lake basin based on FAO’s WaPOR database. Agric. Water Manag. 264, 107509 (2022).10.1016/j.agwat.2022.107509
28. Maghrebi M Iran’s agriculture in the anthropocene Earths Future 2020 8 9 e2020EF001547 10.1029/2020EF001547
Maghrebi, M. et al. Iran’s agriculture in the anthropocene. Earths Future 8(9), e2020EF001547 (2020).10.1029/2020EF001547
29. Rahmani J Danesh-Yazdi M Quantifying the impacts of agricultural alteration and climate change on the water cycle dynamics in a headwater catchment of Lake Urmia Basin Agric. Water Manag. 2022 270 107749 10.1016/j.agwat.2022.107749
Rahmani, J. & Danesh-Yazdi, M. Quantifying the impacts of agricultural alteration and climate change on the water cycle dynamics in a headwater catchment of Lake Urmia Basin. Agric. Water Manag. 270, 107749 (2022).10.1016/j.agwat.2022.107749
30. Emberger, L. La végétation de la région méditérranéenne. Essai d’une classification des groupements végétaux. Revue de Botanique n 503, 642–662 (1930).
31. Zhu Q Modeling soybean cultivation suitability in China and its future trends in climate change scenarios J. Environ. Manag. 2023 345 118934 10.1016/j.jenvman.2023.118934
Zhu, Q. et al. Modeling soybean cultivation suitability in China and its future trends in climate change scenarios. J. Environ. Manag. 345, 118934 (2023).10.1016/j.jenvman.2023.118934
32. Pouladi P Interconnected governance and social barriers impeding the restoration process of Lake Urmia J. Hydrol. 2021 598 126489 10.1016/j.jhydrol.2021.126489
Pouladi, P. et al. Interconnected governance and social barriers impeding the restoration process of Lake Urmia. J. Hydrol. 598, 126489 (2021).10.1016/j.jhydrol.2021.126489
33. Schmidt M Gonda R Transiskus S Environmental degradation at Lake Urmia (Iran): Exploring the causes and their impacts on rural livelihoods GeoJournal 2021 86 5 2149 2163 10.1007/s10708-020-10180-w
Schmidt, M., Gonda, R. & Transiskus, S. Environmental degradation at Lake Urmia (Iran): Exploring the causes and their impacts on rural livelihoods. GeoJournal 86(5), 2149–2163 (2021).10.1007/s10708-020-10180-w
34. Azizpour F Riahi V Tagiloo AA Karimi K Disaster management and rural vulnerability (case study Urmia County) Int. J. Sci. Res. 2015 6 2136 2140
Azizpour, F., Riahi, V., Tagiloo, A. A. & Karimi, K. Disaster management and rural vulnerability (case study Urmia County). Int. J. Sci. Res. 6, 2136–2140 (2015).
35. Data MC Salgado CM Azevedo C Proença H Vieira SM Noise versus outliers Second. Anal. Electron. Health Rec. 2016 10.1007/978-3-319-43742-2_14
Data, M. C., Salgado, C. M., Azevedo, C., Proença, H. & Vieira, S. M. Noise versus outliers. Second. Anal. Electron. Health Rec.10.1007/978-3-319-43742-2_14 (2016).10.1007/978-3-319-43742-2_14
36. Fakhri, M., Bakhshi, K. G., & Sadeghi, S. M. A survey on ecological character of Capparis spinosa L. Bushehr Province (2008).
37. Vrbik J Small-sample corrections to Kolmogorov-Smirnov test statistic Pioneer J. Theor. Appl. Stat. 2018 15 1–2 15 23
Vrbik, J. Small-sample corrections to Kolmogorov-Smirnov test statistic. Pioneer J. Theor. Appl. Stat. 15(1–2), 15–23 (2018).
38. Isaaks EH Srivastava RM Introduction to Applied Geostatistics 1989 Oxford University Press
Isaaks, E. H. & Srivastava, R. M. Introduction to Applied Geostatistics (Oxford University Press, 1989).
39. Kennedy, J., & Eberhart, R. Particle swarm optimization (PSO), in Proc. IEEE International Conference on Neural Networks, Perth, Australia, Vol. 4, No. 1, 1942–1948 (1995).
40. Chen ZS Zhu B He YL Yu LA A PSO based virtual sample generation method for small sample sets: Applications to regression datasets Eng. Appl. Artif. Intell. 2017 59 236 243 10.1016/j.engappai.2016.12.024
Chen, Z. S., Zhu, B., He, Y. L. & Yu, L. A. A PSO based virtual sample generation method for small sample sets: Applications to regression datasets. Eng. Appl. Artif. Intell. 59, 236–243 (2017).10.1016/j.engappai.2016.12.024
