Predicting climate change impacts on the distribution of endemic fish Cyprinion muscatense in the Arabian Peninsula

Freshwater fishes are facing considerable threats in the Arabian Peninsula which is considered as a highly stressed region in the Middle East. It is predicted that northern Oman is likely to face decreasing rainfall and increasing temperature in coming decades. In this study, we focused on an endemic cyprinid fish Cyprinion muscatense, as a model to investigate impacts of climate change on the mountain fishes inhibiting in this arid region. This species is expected to be strongly affected by climate change because of its limited distribution range in a montane area surrounded by lowlands and sea, limiting the species in shift to other areas. We used an ensemble approach by considering two regressions-based species distribution modeling (SDM) algorithms: generalized linear models (GLM), and generalized additive models (GAM) to model the species habitat suitability and predict the impacts of climate change on the species habitat suitability. Based on the distribution models, the montane area located in northeastern Oman was identified as the most suitable habitat for this species. Our results indicate that, even under the minimum greenhouse gas emissions scenario (RCP 2.6), climate change will produce a high reduction in its potential future habitats. According to the results of percent contribution, elevation and annual minimum temperature were the most important variables in predicting the species suitable habitats. Results also showed that only a small percentage of suitable habitats for the species within boundaries of protected areas. Therefore, the impact of climate change on the species appears particularly alarming. Although our study was restricted to a single cyprinid freshwater species, decreases in potential habitats are likely predicted for other cyprinid fish species restricted to the mountains of this region, suggesting severe consideration is needed for aquatic systems in future conservation planning, especially for endemic freshwater fishes.

Mapping current and future risk of scorpion sting from a species with low medical concern, Mesobuthus phillipsii (Scorpiones: Buthidae) in Iran

Scorpion stings are one of the most important health challenges and high priority research topic in public health. In this study, we aimed to model habitat suitability of the Mesobuthus phillipsii (Pocock 1889), a species with low medical concern, under current and future climatic conditions in Iran. We also identified vulnerable populations to scorpion stings in the country. Scorpion sting risk modeling was done using an ensemble approach by considering two species distribution modeling methods: MaxEnt and Random Forest methods. Distribution modeling was performed using the sdm R package. The results showed that due to climate change in 2070, the high-risk areas will increase from 20,839 to 79,212 km2. Habitats with a moderate risk of scorpion stings will also increase from 139,347 to 222,833 km2. Consequently, the number of villages in high-risk areas of scorpion stings will increase from 2,870 to 7,017, while this number will increase from 12,759 to 20,104 in the case of medium-risk villages. The results of this study can be used for scorpion stings management in Iran. This study can be used as an example for similar studies on scorpions with high medical emergency.

Predicting current and future high-risk areas for vectors and reservoirs of cutaneous leishmaniasis in Iran

Climate change will affect the distribution of species in the future. To determine the vulnerable areas relating to CL in Iran, we applied two models, MaxEnt and RF, for the projection of the future distribution of the main vectors and reservoirs of CL. The results of the models were compared in terms of performance, species distribution maps, and the gain, loss, and stable areas. The models provided a reasonable estimate of species distribution. The results showed that the Northern and Southern counties of Iran, which currently do not have a high incidence of CL may witness new foci in the future. The Western, and Southwestern regions of the Country, which currently have high habitat suitability for the presence of some vectors and reservoirs, will probably significantly decrease in the future. Furthermore, the most stable areas are for T. indica and M. hurrianae in the future. So that, this species may remain a major reservoir in areas that are present under current conditions. With more local studies in the field of identifying vulnerable areas to CL, it can be suggested that the national CL control guidelines should be revised to include a section as a climate change adaptation plan.

Monitoring spatiotemporal impacts of changes in land surface temperature on near eastern fire salamander (Salamandra infraimmaculata) in the Middle East

Persistence and coexistence of many pond-breeding amphibians depend on seasonality. Temperature, as a seasonal climate component, affects numerous physical and biological processes of pond-breeding amphibians. Satellite-derived land surface temperature (LST) is the radiative skin temperature of the land surface, which has received less attention in spatiotemporal seasonal habitat monitoring. The present study aims to evaluate the increasing and decreasing effects of LST trends at two levels: (1) habitat suitability and connectivity; (2) individual population sites and their longitudinal distribution (with increasing longitude). Habitat suitability modeling was conducted based on an ensemble species distribution model (eSDM). Using electrical circuit theory, the connectivity of interior and intact habitat cores was investigated. An average seasonal LST was prepared separately for each season from 2003 to 2021 and entered into Mann-Kendall (MK) analysis to determine the spatiotemporal effects of LST changes using the Z-Score (ZMK) at two confidence levels of 95 and 99%. Based on the results, in winter, 28.12% and 70.70% of the suitable habitat were affected by an increasing trend of LST at 95% and 99% confidence levels, respectively. The highest spatial overlap of the decreasing trend of LST with the suitable habitat occurred in summer and was 6.4% at the 95% confidence level and 4.2% at the 99% confidence level. Considering population site at 95% confidence interval, the increasing trend of LST was calculated to be 20.2%, 9.5%, 4.2%, and 6.3% of localities in winter, spring, summer, and autumn, respectively. At the 99% confidence level, these percentages reduced to 8.5%, 3.1%, 1%, and 1%, respectively. During winter and summer, based on the results of the longitudinal trend, an increasing trend of LST was observed in sites. Localities of Hatay and Iica village in Turkey experienced seasonally asynchronous climate change regimes. The approach used in this study allowed us to create a link between the life cycle and seasonal changes on a micro-scale (breeding sites) and macro-scale (distribution and connectivity). Findings of this paper can be effectively used by conservation managers to preserve S. infraimmaculata's metapopulation.

Identifying high snakebite risk area under climate change for community education and antivenom distribution

Snakebite is one of the largest risks from wildlife, however little is known about venomous snake distribution, spatial variation in snakebite risk, potential changes in snakebite risk pattern due to climate change, and vulnerable human population. As a consequence, management and prevention of snakebite is hampered by this lack of information. Here we used habitat suitability modeling for 10 medically important venomous snakes to identify high snakebite risk area under climate change in Iran. We identified areas with high snakebite risk in Iran and showed that snakebite risk will increase in some parts of the country. Our results also revealed that mountainous areas (Zagros, Alborz, Kopet-Dagh mountains) will experience highest changes in species composition. We underline that in order to improve snakebite management, areas which were identified with high snakebite risk in Iran need to be prioritized for the distribution of antivenom medication and awareness rising programs among vulnerable human population.

Medically important snakes and snakebite envenoming in Iran

Snakebite is a common health condition in Iran with a diverse snake fauna, especially in tropical southern and mountainous western areas of the country with plethora of snake species. The list of medically important snakes, circumstances and effects of their bite, and necessary medical care require critical appraisal and should be updated regularly. This study aims to review and map the distributions of medically important snake species of Iran, re-evaluate their taxonomy, review their venomics, describe the clinical effects of envenoming, and discuss medical management and treatment, including the use of antivenom. Nearly 350 published articles and 26 textbooks with information on venomous and mildly venomous snake species and snakebites of Iran, were reviewed, many in Persian (Farsi) language, making them relatively inaccessible to an international readership. This has resulted in a revised updated list of Iran's medically important snake species, with taxonomic revisions of some, compilation of their morphological features, remapping of their geographical distributions, and description of species-specific clinical effects of envenoming. Moreover, the antivenom manufactured in Iran is discussed, together with treatment protocols that have been developed for the hospital management of envenomed patients.

Disentangling the impacts of climate and land cover changes on habitat suitability of common pheasant Phasianus colchicus along elevational gradients in Iran

Climate and land cover change are critical drivers of avian species range shift. Thus, predicting avian species' response to the land and climate changes and identifying their future suitable habitats can help their conservation planning. The common pheasant (Phasianus colchicus) is a species of conservation concern in Iran and is included in the list of Iran's protected avian species. The species faces multiple threats such as habitat destruction, land cover change, and overhunting in the country. In this study, we model the potential impacts of future climate and land cover change on the habitat suitability of common pheasant (Phasianus colchicus) along elevational gradients in Mazandaran province in Iran. We used shared socioeconomic pathways (SSP) scenarios and the 2015-2020 trend to generate possible future land cover projections for 2050. As for climate change projections, we used representative concentration pathway (RCP) scenarios. Next, we applied current and future climate and land cover projections to investigate how habitat suitability of common pheasant will change between 2020 and 2050 using species distribution modeling (SDM). Our results show that the species has 6000 km² suitable habitat; however, between 900 and 1965 km² of its habitat may be reduced by 2050. Furthermore, we found that the severity of the effects of climate and land cover change varies at different altitudes. At low altitudes, the impact of changing land structure is superior. Instead, climate change has a critical role in habitat loss at higher altitudes and imposes a limiting role on the potential range shifts. Overall, this study demonstrates the vital role of land cover and climate change in better understanding the potential alterations in avian distribution.

Iranian wetland inventory map at a spatial resolution of 10 m using Sentinel-1 and Sentinel-2 data on the Google Earth Engine cloud computing platform

Detailed wetland inventories and information about the spatial arrangement and the extent of wetland types across the Earth's surface are crucially important for resource assessment and sustainable management. In addition, it is crucial to update these inventories due to the highly dynamic characteristics of the wetlands. Remote sensing technologies capturing high-resolution and multi-temporal views of landscapes are incredibly beneficial in wetland mapping compared to traditional methods. Taking advantage of the Google Earth Engine's computational power and multi-source earth observation data from Sentinel-1 multi-spectral sensor and Sentinel-2 radar, we generated a 10 m nationwide wetlands inventory map for Iran. The whole country is mapped using an object-based image processing framework, containing SNIC superpixel segmentation and a Random Forest classifier that was performed for four different ecological zones of Iran separately. Reference data was provided by different sources and through both field and office-based methods. Almost 70% of this data was used for the training stage and the other 30% for evaluation. The whole map overall accuracy was 96.39% and the producer's accuracy for wetland classes ranged from nearly 65 to 99%. It is estimated that 22,384 km² of Iran are covered with water bodies and wetland classes, and emergent and shrub-dominated are the most common wetland classes in Iran. Considering the water crisis that has been started in Iran, the resulting ever-demanding map of Iranian wetland sites offers remarkable information about wetland boundaries and spatial distribution of wetland species, and therefore it is helpful for both governmental and commercial sectors.

Global habitat suitability modeling reveals insufficient habitat protection for mangrove crabs

Mangrove crabs are important components of mangrove forests however their large scale habitat suitability and conservation received little attention. The Metopograpsus thukuhar/cannicci species complex is a mangrove dwelling species occurs in the Indo-Pacific mangrove forests. Since identifying the complex suitable habitat is critical for its conservation, we modeled global habitat suitability of the complex within marine biogeographic realms and estimated representation of the complex suitable habitats within marine protected areas. We found that the complex' largest and smallest suitable ranges are located in Central Indo-Pacific and Temperate Southern Africa realms, respectively. Only 12.5% of the complex suitable habitat is protected. The highest proportion of the complex' protected suitable habitat (22.9%) is located in Western Indo-Pacific realm while the lowest proportion of the complex' protected suitable habitat (1.38%) is located in Central Indo-Pacific realm. Suitable unprotected habitats of the complex identified in this study have high priority for conservation and should be included in marine protected areas to ensure species conservation. Our results show that species distribution models are practical tools to study marine species distribution across large spatial scales and help marine conservation planning.

Modeling climate change impacts on the distribution of an endangered brown bear population in its critical habitat in Iran

Climate change is one of the major challenges to the current conservation of biodiversity. Here, by using the brown bear, Ursus arctos, in the southernmost limit of its global distribution as a model species, we assessed the impact of climate change on the species distribution in western Iran. The mountainous forests of Iran are inhabited by small and isolated populations of brown bears that are prone to extinction in the near future. We modeled the potential impact of climate change on brown bear distribution and habitat connectivity by the years 2050 and 2070 under four representative concentration pathways (RCPs) of two general circulation models (GCMs): BCC-CSM1-1 and MRI-CGCM3. Our projections revealed that the current species' range, which encompasses 6749.8 km² (40.8%) of the landscape, will decline by 10% (2050: RCP2.6, MRI-CGCM3) to 45% (2070: RCP8.5, BCC-CSM1-1). About 1850 km² (27.4%) of the current range is covered by a network of conservation (CAs) and no-hunting (NHAs) areas which are predicted to decline by 0.64% (2050: RCP2.6, MRI-CGCM3) to 15.56% (2070: RCP8.5, BCC-CSM1-1) due to climate change. The loss of suitable habitats falling within the network of CAs and NHAs is a conservation challenge for brown bears because it may lead to bears moving outside the CAs and NHAs and result in subsequent increases in the levels of bear-human conflict. Thus, re-evaluation of the network of CAs and NHAs, establishing more protected areas in suitable landscapes, and conserving vital linkages between habitat patches under future climate change scenarios are crucial strategies to conserve and manage endangered populations of the brown bear.

Genetic Diversity Assessment of Iranian Kentucky Bluegrass Accessions: I. ISSR Markers and Their Association with Habitat Suitability Within and Between Different Ecoregions

Poa pratensis L. is a commonly used cool-season turfgrass and endemic to Iran. This research was carried out to examine the genetic diversity of this plant within and between ecoregions of Iran and the impact of climatic variables and elevation on the distribution of its genotypes, as well as habitat suitability modeling. We used fifty accessions collected from six ecoregions (West, South, North, North-West and North-East) for genetic diversity assessment using 20 ISSR marker primers. The prospective ecoregions for Kentucky bluegrass production were projected using habitat suitability modeling, which took into account important environmental parameters, such as annual mean temperature, annual mean rainfall, and elevation. According to the UPMGA dendrogram, the accessions were divided into two major types and four subclasses. The genetic distance between the North and North-east accessions, as well as the Center accessions, was greater than that of the other genotypes. Center accessions had the greatest levels of polymorphism, effective number of alleles, Shannon index, and Nei's genetic diversity. The FR method was used to create the habitat suitability map based on environmental factors. Rainfall had the largest influence on the genotype distribution of P. pratensis L. The findings of this study can be used as raw materials in future breeding programs to improve and generate new cultivars with superior characteristics. It can also assist programs in identifying rare cultivars as well as preserving and developing native P. pratensis L. genotypes.

Biogeography of rodents in Iran: species richness, elevational distribution and their environmental correlates

Rodent biogeographic studies are disproportionately scarce in Iran, however, they are an ideal system to understand drivers of biodiversity distributions in the country. The aims of the present research are to determine (i) the pattern of rodent richness across the country, (ii) quantify their elevational distribution patterns, and (iii) explore the underlying mechanisms. To reach these goals, an updated species list was compiled based on the latest taxonomic revisions, published until December 2021. We mapped all 76 rodent species distributions to develop the first map of rodent richness in Iran. We furthermore investigated their elevational distribution patterns in the following four geographic regions based on 100 m intervals; the Zagros Mountains, north of Alborz Mountains, south of Alborz and Kopet-Dagh mountains, and central and east mountains. North-east of Iran, Zagros Mountains, Alborz Mountains, and northwestern of Iran showed the highest richness and were identified as biodiversity hotspots of rodents in the country. This study highlights the importance of past climate change as the key driver of rodent richness in Iran. We showed that rodents’ elevational distribution patterns differ among geographic regions. Areas and elevational zones with the highest species richness should be prioritized for the conservation planning of rodents in Iran.

Modeling the effects of climate change on the potential distribution of the rangeland species Gymnocarpus decander Forssk (case study: Arid region of southeastern Iran)

The phenomenon of climate change is the biggest environmental challenge in the world. Climate is a determinant factor in species distribution, and climate change will affect the species' abilities to occupy geographic regions. In this study which was conducted in May of 2019, spatio-temporal changes in potential habitats of Gymnocarpus decander were assessed using the MRI-CGCM3 climate change model for RCP2.6 and RCP8.5 scenarios for the near future (2041-2061) and far future (2061-2080) periods for this purpose, climatic variables of 24 synoptic stations across a case study, bio-climatic data and vegetation cover maps of G. decander were used. First, using the factor analysis process, the dimensions of the station-observed climatic variables were reduced to five factors with a total variance of 88.3%. Then, the region was divided into five homogeneous climatic regions using partitional clustering analysis. In this study by using the logistic regression modeling technique, the probability of the presence of the desired species for two groups of independent variables including climatic factors and bioclimatic variables in each of the groups was modeled. The results showed that the best models for determining the potential habitats of G. decander are logistic regression models in groups with independent bioclimatic variables. According to the results obtained from both scenarios, the habitats of G. decander species will decrease in the future. In the most optimistic case, about 8% of G. decander habitats will be lost by 2060 and about 12% by 2080. According to modeling results, currently, 48.2% total area of the region under study has a high potential for the presence of G. decander. Also, results indicate that region number 4 in this study with an altitude range of about 800-1250 m, 16 °C average temperature in the growing season and annual precipitation around 150-170 mm is the major habitat for G. decander. According to climate change under the RCP2.6 scenario, the area of potential habitats of G. decander will decrease to 40% in the near future and 36.4% in the far future; and according to climate change under the RCP8.5 scenario, the area of potential habitats of G. decander will decrease to 23.9% in the near future and 32.5% in the far future. In the far future, because of the increase in total precipitation, some of the lost potential habitats during the near future will be suitable again for G. decander. Due to its stability in harsh environmental conditions, G. decander appears as a type-forming species in a wide range of natural habitats in the study area and is therefore important in terms of soil protection and forage production.

The risk of forfeiting the ranges of reptiles under nonrandom and stochastic scenarios of moving climate conditions: a case study for 115 species in China

Revealing the hazard features of forfeiting areal ranges for nonidentical scenarios of shifting climatic conditions is pivotal for the conformation of reptiles to climatic warming. Taking 115 reptiles in China as an example, the indefiniteness and danger of shrinking geographical range for the reptiles under stochastic and nonrandom scenarios of moving climatic situations were inspected via exploiting the scenarios of shifting climatic status associated with the representative concentration pathways, Monte Carlo simulation, and the classifications scheme based on the fuzzy set. For non-stochastic states of altering climatic elements, the richness of 115 reptiles improved in certain sites of northeastern, and western China and dropped in several areas of northern, eastern, central China, and southeastern China: roughly 59-74 reptiles forfeiting less than 20% of their present ranges, roughly 25-34 reptiles narrowing less than 20-40% of their present areal ranges, and roughly 105-111 reptiles inhabited more than 80% of their overall areal ranges. For the random status of shifting climatic elements, the count of reptiles that forfeited the various extent of the present or entire areal ranges descended with raising the eventuality; with a possibility of over 0.6, the count of reptiles that minified less than 20%, 20-40%, 40-60%, 60-80% and over 80% of the present ranges was roughly 28-49, 5-10, 1-3, 0-1 and 13-18, separately; the count of reptiles that inhabited below 20%, 20-40%, 40-60%, 60-80% and more than 80% of the entire real ranges was roughly 0-1, 5-6, 1-5, 0-2 and 35-36, separately. About 30% of 115 reptiles would face disappearance danger in response to moving climate conditions in the absence of adaption steps, and the conformation measures were indispensable for the reptiles that shrunk their areas.

Design of Protected Area by Tracking and Excluding the Effects of Climate and Landscape Change: A Case Study Using Neurergus derjugini

This study aimed to use the applications of Ensemble Species Distribution Modelling (eSDM), Geographical Information Systems (GISs), and Multi-Criteria Decision Analysis (MCDA) for the design of a protected area (PA) for the critically endangered yellow-spotted mountain newt, Neurergus derjugini, by tracking and excluding the effects of climate and landscape changes in western Iran and northeastern Iraq. Potential recent and future distributions (2050 and 2070) were reconstructed by eSDM using eight algorithms with MRI-CGCM3 and CCSM4 models. The GIS-based MCDA siting procedure was followed inside habitats with high eSDM suitability by eliminating the main roads, cities, high village density, dams, poor vegetation, low stream density, agricultural lands and high ridge density. Then, within the remaining relevant areas, 10 polygons were created as “nominations” for PAs (NPAs). Finally, for 10 different NPAs, the suitability score was ranked based on ratings and weights (analytical hierarchy process) of the number of newt localities, NPA connectivity, NPA shape, NPA habitat suitability in 2070, NPA size, genetic diversity, village density and distance to nearest PAs, cities, and main roads. This research could serve as a modern realistic approach for environmental management to plan conservation areas using a cost-effective and affordable technique.

Accessing habitat suitability and connectivity for the westernmost population of Asian black bear (Ursus thibetanus gedrosianus, Blanford, 1877) based on climate changes scenarios in Iran

Climate change, as an emerging phenomenon, has led to changes in the distribution, movement, and even risk of extinction of various wildlife species and this has raised concerns among conservation biologists. Different species have two options in the face of climate change, either to adopt or follow their climatic niche to new places through the connectivity of habitats. The modeling of interpatch landscape communications can serve as an effective decision support tool for wildlife managers. This study was conducted to assess the effects of climate change on the distribution and habitat connectivity of the endangered subspecies of Asian black bear (Ursus thibetanus gedrosianus) in the southern and southeastern Iran. The presence points of the species were collected in Provinces of Kerman, Hormozgan, and Sistan-Baluchestan. Habitat modeling was done by the Generalized Linear Model, and 3 machine learning models including Maximum Entropy, Back Propagation based artificial Neural Network, and Support Vector Machine. In order to achieve the ensemble model, the results of the mentioned models were merged based on the method of "accuracy rate as weight" derived from their validation. To construct pseudo-absence points for the use in the mentioned models, the Ensemble model of presence-only models was used. The modeling was performed using 15 habitat variables related to climatic, vegetation, topographic, and anthropogenic parameters. The three general circulation models of BCC-CSM1, CCSM4, and MRI-CGCM3 were selected under the two scenarios of RCP2.6 and RCP8.5 by 2070. To investigate the effect of climate change on the habitat connections, the protected areas of 3 provinces were considered as focal nodes and the connections between them were established based on electrical circuit theory and Pairwise method. The true skill statistic was employed to convert the continuous suitability layers to binary suitable/unsuitable range maps to assess the effectiveness of the protected areas in the coverage of suitable habitats for the species. Due to the high power of the stochastic forest model in determining the importance of variables, this method was used. The results showed that presence/absence models were successful in the implementation and well distinguished the points of presence and pseudo-absence from each other. Based on the random forests model, the variables of Precipitation of Driest Quarter, Precipitation of Coldest Quarter, and Temperature Annual Range have the greatest impact on the habitat suitability. Comparing the modeling findings to the realities of the species distribution range indicated that the suitable habitats are located in areas with high humidity and rainfall, which are mostly in the northern areas of Bandar Abbas, south of Kerman, and west and south of Sistan-Baluchestan. The area of suitable habitats, in the MRI-CGCM3 (189731 Km2) and CCSM4 (179007 Km2) models under the RCP2.6 scenario, is larger than the current distribution (174001 Km2). However, in terms of the performance of protected areas, the optimal coverage of the species by the boundary of the protected areas, under each of the RCP2.6 and RCP8.5 scenarios, is less than the present time. According to the electric circuit theory, connecting the populations in the protected areas of Sistan-Baluchestan province to those in the northern Hormozgan and the southern Kerman would be based on the crossing through the heights of Sistan-Baluchestan and Hormozgan provinces and the plains between these heights would be the movement pinch points under the current and future scenarios. Populations in the protected areas of Kerman have higher quality patch connections than that of the other two provinces. The areas such as Sang-e_Mes, Kouh_Shir, Zaryab, and Bahr_Aseman in Kerman Province and Kouhbaz and Geno in Hormozgan Province can provide suitable habitats for the species in the distribution models. The findings revealed that the conservation of the heights along with the caves inside them could be a protective priority to counteract the effects of climate change on the species.

Applying species distribution models in public health research by predicting snakebite risk using venomous snakes' habitat suitability as an indicating factor

Snakebite envenoming is an important public health problem in Iran, despite its risk not being quantified. This study aims to use venomous snakes' habitat suitability as an indicator of snakebite risk, to identify high-priority areas for snakebite management across the country. Thus, an ensemble approach using five distribution modelling methods: Generalized Boosted Models, Generalized Additive Models, Maximum Entropy Modelling, Generalized Linear Models, and Random Forest was applied to produce a spatial snakebite risk model for Iran. To achieve this, four venomous snakes' habitat suitability (Macrovipera lebetinus, Echis carinatus, Pseudocerastes persicus and Naja oxiana) were modelled and then multiplied. These medically important snakes are responsible for the most snakebite incidents in Iran. Multiplying habitat suitability models of the four snakes showed that the northeast of Iran (west of Khorasan-e-Razavi province) has the highest snakebite risk in the country. In addition, villages that were at risk of envenoming from the four snakes were identified. Results revealed that 51,112 villages are at risk of envenoming from M. lebetinus, 30,339 from E. carinatus, 51,657 from P. persicus and 12,124 from N. oxiana. Precipitation seasonality was identified as the most important variable influencing distribution of the P. persicus, E. carinatus and M. lebetinus in Iran. Precipitation of the driest quarter was the most important predictor of suitable habitats of the N. oxiana. Since climatic variables play an important role in shaping the distribution of the four venomous snakes in Iran, thus their distribution may alter with changing climate. This paper demonstrates application of species distribution modelling in public health research and identified potential snakebite risk areas in Iran by using venomous snakes' habitat suitability models as an indicating factor. Results of this study can be used in snakebite and human-snake conflict management in Iran. We recommend increasing public awareness of snakebite envenoming and education of local people in areas which identified with the highest snakebite risk.

Making Future Teachers More Aware of Issues Related to Sustainability: An Assessment of Best Practices

The aim of the present research was to uncover the opinions of future teachers about whether the formative knowledge of education for sustainability received in seminars during their practicums improved their competencies in sustainability and their values in educational teaching for their future classroom practice. The study entailed qualitative research with a thematic analysis. Semi-structured interviews comprising 14 questions, which were validated by 15 experts, were used to solicit the opinions of 52 future teachers with respect to the importance of sustainability and the development of intercultural competencies among their students. The information from the interviews was analyzed using thematic analysis on the basis of education for sustainable development (ESD). We conclude that the teachers’ practicum placements and training experiences were very positive because a lecturing pedagogical approach, taught via practical application in a school environment by an expert of recognized prestige, can be a useful resource to develop awareness of both sustainable development and its education, as well as to learn didactic strategies to apply ESD, which addresses aspects that are relevant in multicultural contexts, such as tolerance and empathy.