Journal of Spatial Analysis

---

 
Evaluation and comparison of the accuracy of fault and seismic data in fractal analysis of northwest Zagros tectonic
Introduction
Complexity of natural processes especially tectonic processes that shape landscapes cannot be evaluated by classic geometry. In comparison with integer dimension of Euclidean space, fractal geometry can analyze features with non-integer dimension (Turcotte, 1977:121). Fractal behavior in such features shows self-similarity that this component is independent of the accuracy of investigation (Baas, 2002, 311). In fact, fractal dimension, is scale-invariant (Phillips, 2002, 144). Spatial variations of fractal parameters are an important factor in studying the tectonic state of regions. By determining the fractal dimension of Linear structures such as faults, it is possible to compare their geometry disorder (Suk moon et al, 1996:5). This parameter affects seismic behavior of fault because earthquake is an event related to faulting (Bachmanov, et al, 2012: 221) and Their concentration in an area indicates tectonic activity. In this research we performed fractal analysis using box counting method on fault and seismic data of northwest of Zagros about different scales of fault and different time periods of earthquake epicenters of two organizations with various detail to find and examine their fractal behavior by fractal dimension.
Methods
Data in this research can be divided to three clusters: 1. Fault lines of two scales of geology maps (1:100000 and 1:250000), 2. Earthquake epicenters of two periods of times prepared by two organizations (20 century data of Institute of Geophysics and 1900-2020 data of International Institute of Earthquake Engineering and Seismology) and 3. The second cluster with exert of Magnitude of completeness of earthquakes that show the minimum magnitude above which the data in the earthquake catalog is complete. Fractal analysis applied on these data by box counting method. To achieve this goal firstly, under study area divided to 6 boxes that contain main fault trends horizontally and vertically (A: folded Zagros in west of Kermanshah, B: faulted Zagros around Kermansha and east of kermansha, C: folded Zagros near mountain front fault, D: An area between faulted and folded Zagros near Khoramabad, E: Area around Balarud fault and F: An area between Balarud and mountain front fault to faulted Zagros). To calculate fractal dimension of fault lines and distribution of earthquake epicenters, box counting method suggested by Turcotte (1997) were applied by using Hausdorff dimension, which in two quantity of size (side length of grids) and number (number of grid boxes containing earthquake epicenter or fault) are used to calculate FD (total fractal dimension) value (Schuller et al, 2001: 3). Relation between reciprocal of side length (quantity of size) and number of boxes containing point and linear features (quantity of Number) was drawn Logarithmically as a linear regression in Excel that shows fractal dimension.
Result and discussion
Larger values of fractal dimension indicate greater geometric disorder (Sukmono et al., 1996: 5). Analysis of faults of two scales represent that faults geometry is fractal and the amount of FD for scale of 1:100000 compared with scale of 1:250,000 is larger but their result approximately is same. The FD values for both scales are locate between 1 and 2 that expresses development of the fractal community of faults has a linear trend. On the other hand, for earthquakes, increase in FD values shows that earthquakes are not clustered and are distributed homogeneously (Oncel & Wilson, 2002: 339) along a line in understudy area. Calculated number-size values for faults and earthquakes represent both partial and popular FD changes. Based on partial FD, two populations can be classified: (a) Background with FD larger than popular FD; (b) Threshold with FD lower than popular FD.
Conclusion
Fractal analysis of faults of two scales of geology maps reveals that the order of active areas with high FD values in both scales are same but due to different details of faults in geology maps of geology survey and oil company, in scale of 1:100000 area labeled B and in scales of 1:250000 area labeled A is the most tectonically active region, however, area labeled E in both scales has lowest value. The order of active areas based on FD values for earthquake epicenters of 1900-2021 data of geophysics institute do not support other results because area labeled C with low density of faults and earthquake epicenters is in the first order and area labeled A is on the contrary of it. However, FD results of 20 century earthquake epicenters with exert of magnitude of completeness are reliable and higher magnitude of earthquakes spatially recent Ezgeleh earthquake in area labeled A is its evidence.
Keywords: Fractal, Tectonic, Northwest Zagros, Fault, Earthquake
 

---

Investigation of Temperature and Precipitation Changes in the Seymarreh Basin by Using CMIP5 Series Climate Models
 
Abstract
Panel reports on climate change suggest that climate change around the world is most likely due to human factors. Temperature and precipitation are two important parameters in the climate of a region whose variations and fluctuations affect different areas such as agriculture, energy, tourism and so on. Seymareh basin is one of the most significant sub-basins of Karkheh. The purpose of this study is to predict the impact of climate change on precipitation and temperature of the Seymareh Basin in 2021-2040 period. These effects were analyzed at selected stations with uncertainties related to atmospheric general circulation models (GCMs) of CMIP5 models under two scenarios of RCP45 and RCP85 through LARS-WG statistical model. Then the uncertainties of the models and scenarios were investigated by comparing the monthly outputs of the models by the coefficients of determination coefficient (R²) in the forthcoming period (2021-2040) with the base period (1980–2010). The root mean square error (RMSE) calculations presented the best model and scenarios for generating future temperature and precipitation data.            
The Seymareh catchment is the largest and the main Karkheh sub-basin that covers parts of Kermanshah, Lorestan and Ilam provinces. The length of the largest river at the basin level to the site of the Seymareh Reservoir Dam is approximately 475 km, and the area of the basin is 26,700 km². Geographic coordinates of the basin are from 33° 16 ́ 03 ̋to 34°59 ́ 29 ̋north latitudes and 46°6 ́9 ̋to ̋ 5 ́ 0 ° 49 Eastern longitudes, minimum basin height 698 m at the dam outlet and its maximum height 3,638 m. It is on the western highlands of Borujerd.
The information used in this study was obtained from the Meteorological Organization of the country. For this study, three synoptic stations of Kermanshah, Hamadan and Khorramabad, which had the highest statistical records and had appropriate distribution at basin level, were used. These data included daily and monthly temperature and precipitation information, and sunshine hours.
The LARS-WG fine-scale exponential model was proposed by Rasko et al., Semnoff and Barrow (1981). We used daily data at stations under current and future weather conditions. In order to select the best GCM model from the models mentioned above, minimum temperature, maximum temperature, precipitation and sunshine data were entered daily in the base period (1980–2010) and data were generated for five models under two scenarios of RCP45 and RCP85 for the period 2040–2021. The data were generated in 100 random series and the mean of required variables (minimum temperature, maximum temperature and rainfall) were extracted monthly in the period 2021-2040. Then, root mean square error (RMSE) and determination coefficient (R²) were used to evaluate the performance of the models and compare the results.
To ensure the models' ability to generate data in the coming period, computational data from the model and observational data at the stations under study should have been compared. The capability of the LARS-WG model in modeling the minimum temperature, maximum temperature, and radiation at the stations under study was completely consistent with the observed data. The model's ability to exemplify rainfall was also acceptable, however the highest modeling error was related to March rainfall.
By comparing the observed and produced data including monthly average precipitation, minimum and maximum temperatures through five mentioned models with their indices, the best model and scenario for future fabrication were determined. The results of this comparison showed that among the available models, HADGEM2-ES model under RCP 4.5 scenario had the best result for precipitation and HADGEM2-ES under RCP 8.5 scenario predicted the best result for maximum temperature. Determining the best model, precipitation data, minimum temperature and maximum temperature produced in the selected models and scenarios were analyzed to investigate the climate change temperature and precipitation for the future period.
The results of this study indicated that due to the wide range of output variations of different models and scenarios, by not taking into account the uncertainties of the models and scenarios can have a great impact on the results of the studies. It was also found in this study that the LARS-WG exponential model was capable of modeling precipitation data and baseline temperature in the study area, so that the radiation data, minimum and maximum temperatures were completely consistent with the data.
The observations are consistent and the models' ability to predict rainfall is very good and acceptable manner. In investigating the uncertainties caused by atmospheric general circulation models and existing scenarios, the best model to predict precipitation in the study area is HADGEM2-ES model under RCP 8.5 scenario, the best model for temperature estimation model HADGEM2-ES under RCP scenario No. 4.5.
The overall results of this study revealed that the average precipitation in the basin will decrease by 4.5% on average, while the minimum temperature will be 1.5° C and the maximum temperature will be 2.17° C. The highest increase will be due to the warmer months of the year. Notable are the disruptions of rainfall distribution and the high temperatures will have significantly negative consequences than rainfall reduction.
 

• : Climate Change, Climate Scenarios, Uncertainty, LARS-WG, Seymareh.

 
 

---

Prioritization analysis of effective factors in non-participation of local societies in desertification projects (Case Study: Ain Khosh region, Ilam province)
 
Introduction
Participation in social affairs is a commitment and acceptance of individual and social responsibility that all human beings will have to accept. This commitment and responsibilities may take the form of definite and unlimited activities. By increasing the population and the complexity of the goals and efforts of the human community to advance economic, cultural, social and political goals, we inevitably need partnership and cooperation. Participation means using personal resources to participate in a collective action. The first step is to increase popular participation in desertification initiatives, identify and remove barriers to effective non-participation in project implementation. The main objective of this research is to prioritize the factors affecting the lack of public participation in desertification plans using the FUZZY-AHP method and the Friedman nonparametric test.
 
Materials and methods
The statistical population of the study consisted of households in Ein-e-Khosh village of Dehloran Ilam and experts of Ilam University and natural resources organization of Ilam province and Dehloran County with more than ten years’ experience in combating desertification issues. In this research, the indexes and sub-indicators related to library studies, questionaire from experts of the university, experts from the Natural Resources Department of Ilam province and Dehloran city, as well as referring to the region and interviews with the residents of the region were identified. Then, the questionnaire designed by the FAHP method evaluated by inconsistency rate and its validity and reliability by Likert scale, and finally tried to prioritizing them based on the following steps. First, the prioritization of the indicators was performed from the expert's point of view using Fuzzy Analytical Hierarchy Process (FAHP). Then, the prioritization of the indicators from the perspective of experts by application of Friedman test and finally, the priority of indices and sub-indicators by the local point of view with Friedman's test.
 
Results and discussion
The findings of the research showed that the ranking of indices using Friedman's nonparametric test is based on the average rating from the viewpoint of residents of the region, respectively, economic index, design-executive, educational-promotional, and social. Also, this prioritization from the perspective of experts using the FUZZY-AHP test is design-executive, economic, educational, promotional and social priority, respectively. Also, the results showed that in the total of 15 identified subcategories influencing the lack of public participation in combating desertification projects, from the perspective of experts, the sub-design of the design-executive entitled "Short-term, mid-term and long-term non-planning for participation" with the average rating of 11.68 was the most important and the sub-index, "Migration of youth in the countryside" with the average rating of 3.59, is the most insignificant sub-indicator. However, from the perspective of residents in the region, the underlying economic indexes "disregarding people's income as a direct incentive to implement combating desertification projects" with an average of 11.24, the most significant and sub-indicator of the design-executive "lack of full allocation of funds during the implementation of combating desertification projects" with average rating of 5.63 is the most significant sub-indicator, which indicates that economic indicators and design-executive, along with the sub-indexes, are the most important reasons for non-participation of people in combating desertification projects in the study area.
Due to the fact that the indicators and sub-indicators are identified based on the opinions of experts and locals in the study area, this has led to familiarizing the respondents with the research. In this research, the FAHP and Friedman test were used. According to the topic of the research in the field of public participation, the best tool for measuring the comprehensive statistical view of experts including experts and locals with regard to the study area is considered. In the reliability of the FAHP questionnaire, the responsiveness questionnaire has high reliability with regard to the multi-stage and multi-stage couples comparing method and the incompatibility rate test (mean 0.043 inconsistency rate). Cronbach's alpha test was used for Likert scale questionnaires. Results (Cronbach's alpha = 0.83) showed that the questionnaire had acceptable reliability. So, results is in consistence with other researchers' findings, including Saleh Pourjem et al. (2017).
 
Conclusion
According to the results obtained from prioritization, it has been shown that in the subject of participation, in spite of the difference between the views of experts and the people of the region, in some cases, the main priorities in the discussion of non-participation are almost similar; these results are consistent with the studies of previous researchers such as Saleh Pourjem et al. (2017).
It is suggested that the removal of obstacles to public participation in combating desertification projects be put on the agenda of trusteeship organizations and public participation in all stages of design, implementation and future protection in the combating desertification projects will be considered.
 
Keywords: People’s Participation, MCDM, FUZZY-AHP, Friedman Test
 

---

Risk modeling of plant species diversity and extinction in Sorkheh_hesar National Park
 
Zahra Mosaffaei¹, Ali Jahani^2*, ³MohammadAli ZareChahouki, ⁴Hamid GoshtasbMeygoni, ⁵Vahid Etemad
 
¹ Masters of Natural Resources Engineering, Environmental Sciences, College of Environment, Karaj
^*2Associate Professor, Department of Natural Environment and Biodiversity, College of Environment, Karaj.
³ Professor, Department of Restoration of arid and mountainous regions, University of Tehran, Karaj
⁴ Associate Professor, Department of Natural Environment and Biodiversity, College of Environment, Karaj
⁵ Associate Professor, Department of Forestry and Forest Economics, University of Tehran, Karaj
 
 
Abstract
Full identification of hazards and prioritizing them for non-harm to nature is one of the first steps in natural resource management. Therefore, introducing a comprehensive system of evaluation, understanding, and evaluation is essential for controlling hazards. This study aimed to model and predict environmental hazards following increased degradation in natural environments by ANN. Thus, 600 soil and vegetation samples were collected from inhomogeneous ecological units. Soil samples were prepared by strip transect method according to soil depth in four profiles (5, 10, 15, 20 cm). Vegetation samples were also collected using a minimum level method using 2 2 square plots according to the type, density, and distribution of vegetation. Sampling was done in two safe zones and other uses were modeled using ANN in MATLAB environment. The optimal model of multilayer perceptron with two hidden layers, sigmoid tangent function and 19 neurons per layer and coefficient of determination of 0.90. The results of sensitivity analysis showed that soil moisture content would be effective in decreasing biodiversity and flood risk as well as increasing the risk of extinction of endemic species in the region, and then the apparent and true gravity and soil porosity and distance from the road play a key role in the degradation of cover. Vegetation has increased flooding and extinction risk. Therefore, it is recommended that measures related to soil and vegetation restoration in this park be taken to reduce future damages as soon as possible.
 
Keywords: Modeling, Artificial Neural Network, Environmental Hazards, National Park, Vegetation
 

---

Despite extensive studies on the relationship between land subsidence and groundwater level, less research were focused on the impacts of distance to pumping wells on variation of land subsidence area. This study presented the linkage between the ground surface displacement rate and groundwater pumping area and the associated geomorphic consequences. The land subsidence rate was extracted from Sentinel-1A images. Then, to evaluate the relationship between the ground surface displacement extent and distance from the pumping wells, 30 pumping wells were identified within the study area. Different buffers at specified distances (500, 700, 1,000, 1,300, 1,500, 2,000 m) were created around each well. To test the effect of the distance to the pumping wells on the spatial extent of critical and slight subsidence areas, average annual images of land subsidence were classified into two classes, including areas with a maximum subsidence rate and a minimum subsidence rate. Further, earth fissure identified by GPS were transformed to the land subsidence classification map. The results showed that there is a significant relationship between the distances to pumping wells and displacement extent. The spatial extent of areas with the maximum subsidence rates decreased as the distance from the pumping wells increased. By contrast, the spatial extent of areas occupied by the minimum subsidence rates increased with increasing the distance from the pumping wells. Also, the density distribution of the earth fissures increased in areas with the maximum subsidence rate.

---

Geography and the Paradigm of Sustainable Development
 
Extended Abstract
Geography and sustainable development
The relation between society and environment has gone through different phases. During the years before the World War II, the environmental determinism controlled this relation. However, after the 1950’s the anthropocentrism replaced the environmental determinism and humans began to overuse the nature in such a way that nature lost its sustainability and many hazards and crises occurred. These destructions were so intense and widespread that some researchers compared with the episodes of geologic time and named the era beginning from 1970’s the Anthropocene epoch. During this period, the planetary boundaries were crossed in some areas such climate change, nitrogen cycle and biodiversity. Climate change has created most of other hazards.
To overcome these problems in 1978 the Brandtland report   announced the sustainable development as not to spend resources more than the nature’s production capacity and not to pollute the nature more than it can assimilate. In other words, the nature should remain in its sustainable state so that the future generations can live with no limitations. The principles of the sustainable development were defined in the earth summit of Rio in 1992 such as social equity, economic viability, and environmental sustainability. These principles were broken down in 17 goals. The Rio conference asked all countries to achieve the sustainable development goals by 2030. 
Methodologically the sustainable development requires integrated multidisciplinary approach to investigate the complex system of human- environment in different temporal and spatial scales to achieve the social equity, economic viability, and environmental sustainability. For this reason, many disciplines such as natural resources, environmental sciences, ecology and geography have contributed to the field. Different data from natural resources, human needs and drivers and environmental changes are required to process in very complicated models. In addition to different variables, the hazards are very important component of the sustainable development research, which also requires multi-aspect complicated approach and models. Spatial dependency is another aspect of sustainable development as it differs from place to place in many characteristics. In brief, from the spatial perspective the sustainable development is an integrated multi-approach research about the human-environment system to establish the sustainability on the earth. All of the related fields should study the sustainable development in collaboration with each other. However, the geography due to its long history of studying the relation between human and environment and its spatial dependency is the best single scientific field which can afford studying the sustainable development. Since the earliest times geography has developed quantitative methods, spatial techniques such as geostatistics, and environmental ethics to conserve the nature and human prosperity. The multi approach and systematic works are the main characteristics of Geography. On the other hand, Geography’s vision of defining the location for human’s activities while saving the nature’s sustainability covers the sustainable development completely. Therefore, geography is the overarching field for the sustainable development and it is the main intention of geographers to research and implement the sustainable development to reduce the environmental hazards and develop the sustainable environment for all the human beings at present and in the future. Geography studies the sustainable development through three steps including spatial analysis, spatial planning, and spatial management. In addition, geographers should learn different skills such remote sensing, multivariate statistics and above all develop a common language between different branches of geography.
 
Keywords: geography, sustainable development, environmental ethics, human nature relationship, Anthropocene, planetary boundaries, sustainability.
 
 

---

Abastract
Introduction
Road accidents is the outcome of driver behavior, road condition, vehicle status, and environmental factors. Therefore, identification and assessment of effective parameters on road accidents can be considered as an appropriate way to reduce the accident events, driving violations and increase the road safety. Determining the effects of meteorological factors on the road accident events has gained more attention in recent years.
  The The main objective of this study was to investigate the relationship between the number of road accidents and the meteorological variables in the intercity road of Grmi-Ardabil in the Barzand route.
Methodology:
In this regard, the effects of climatic factors (including rainfall amount, the minimum absolute temperature, and the number of frost days) on the frequency of perilous events were analyzed. The data of accident events (in recent 4 years) were obtained from the trooper department of Ardabil Province along with the meteorological parameters of Germi station through a 11-year period. The statistical tests were performed using R programming software through statistical analysis.
Findings and Discussion:
The results showed that the majority of accidents were occurred in winter season which is in consistent with the frequency of frost days and also corresponded to the absolute minimum temperature. According to the results, the highest significant positive correlation at (R²= 0.43) was observed between the number of injured people and frost days. In addition, the relationship between the absolute minimum temperature and the number of were identified as significant negative correlation.
Conclusion:
As a concluding remark, the poor road conditions caused by climate element can be considered increasing the frequency of accident events. Accordingly, the proper strategies related to behavior change could be
considered in setting the rules and regulations to reduce the accidents and the number of injuries.

Keywords: Climatic hazards, Correlation analysis, Frost days, Minimum absolute temperature, Germi-Ardabil road


 

---

Introduction
Land use reflects the interactive characteristics of humans and the environment and describes how human exploitation works for one or more targets on the ground. Land use is usually defined on the basis of human use of the land, with an emphasis on the functional role of land in economic activities. Land use, which is associated with human activity, is undergoing change over time. Land use information and land cover are important for activities such as mapping and land management. Over time, land cover patterns and, consequently, land use change, and the human factor can play a major role in this process. Today, satellite-based measurements with geographic information systems are increasingly being used to identify and analyze land-use change and land cover. With regard to the problems of changes and transformations in the studied area, remote sensing can allow managers to categorize images and evaluate land use changes, in addition to saving time and costs, which allows planners to make plans based on changes, more resources are lost. To be prevented.

Materials & Methods
In order to classify and detect the marginal land of the river, TM and OLI image images were selected for a specific month (August, August) for the years 1987 and 2017. The purpose of this study was to investigate the changes occurring in the studied area with an emphasis on agricultural lands. To do this, the images before processing in the ENVI software took radiometric, atmospheric and geometric corrections on them. After that, the main components of the river route were extracted. Five basic algorithms were used to classify the base pixel, but eCognition software was used to classify the object. Supervised classification identifies homogeneous regions with examples of land use and land cover, in which pixels are assigned in known information classes. Education is a process that determines the criteria for these patterns. Learning output is a set of spectral signatures of proposed classes. The first step in object-oriented classification is the segmentation of the image and the creation of distinct objects, consisting of homogeneous pixels. The main purpose of image segmentation is to combine pixels or small objects to create large image objects based on the spectral and spatial characteristics of the image. In order to evaluate the accuracy and compare the resulting maps, the overall accuracy and Kappa coefficient are used. When the sampling of pixels is done as a spectral or informational class pattern, the evaluation of the spectral reflection of classes and their resolution can also be done. An algorithm with the highest accuracy and accuracy will be the basis for the detection. Detection of changes, which leads to a two-way matrix and shows variations of the main types of land use in the study area, was carried out in this study. Pixel-based cross-tabulation analysis on pixels facilitates the determination of the conversion value from a specific user class to another user category and areas associated with these changes over the given time period.

Results & Discussion
The results showed that the object-oriented method is more accurate than the base pixel algorithms for providing user-defined maps. The amount of accuracy in the method based on object-oriented classification depends largely on choosing the appropriate parameters for classification, defining the rules, and applying the appropriate algorithm to obtain the degree of membership. The Kappa coefficient for each image is approximately 0.90. So these maps are the basis for the discovery of change. According to the results, the agricultural and residential lands have been increased and this increase has been accompanied by a decrease in rangelands. A general overview of this 30-year period shows that the arable and dry farming, respectively, increased by 2418.79 and 719.61 hectares and the rangelands had a decrease of 2848.86 hectares. However, the residential class and human effects show an increase of 428.88 hectares or a growth of 178.87%, which indicates the importance of agriculture in the studied area.

Conclusion
Identifying and discovering land cover changes can help planners and planners identify effective factors in land use change and land cover, and have a useful planning to control them. For this reason, maps are needed with precision and speed, and object-oriented processing methods make this possible with very high precision. The results of this study, in addition to proving the precision and efficiency of object-oriented processing in land cover estimation, between 1987 and 2017, have witnessed a decrease in the area of rangeland lands and, on the other hand, agricultural and residential lands, which is indicative of the overall trend Destruction in the area through the replacement of pastures by other uses such as rainfed farming.

Keywords: Land Use, Gamasiab, Object Oriented, Pixel Base, Kappa Coefficient
 

---

Introduction
Rural areas are more vulnerable to earthquake hazards than urban areas but the vulnerability of rural areas has always been neglected and few studies are worrying about it. Given the importance of villages and played the crucial role in socio-economic development and national security, providing adequate housing for villagers and addressing the problems in this area, in particular, providing them with security and relieving their vulnerability are of particular importance. One of the policies of Iran to reduce the risk of damage; improvement and rehabilitation of rural housing by the Housing Foundation of the Islamic Revolution Which has become one of the most important strategies in Iran due to the extent of natural disasters and their financial and financial losses and their expectations beyond ensuring security against accidents are also a continuation of rural life. neyriz Township is subjected to major and minor faults that the existence of these faults and the probability of earthquake causes vulnerability of the villages of the region. The housing estate of the Islamic Revolution of the Islamic Republic of Iran, from 2004 to 2013, has provided 5255 villagers of more than 20 households with facilities for the renovation and renovation of houses. The number of facilities was paid to 66 villages of more than 20 households in the city and supervision of the construction process was carried out. So what seems to be important is the activity that the Housing Foundation has had in its housing estate, its impact on the rural areas, and how much it has been able to achieve resistance and stability in rural housing; in addition to what degree, they have been able to influence their satisfaction. Therefore, the present study addresses the vulnerability of rural dwellings. In this regard, vulnerability is initially studied then the satisfaction of the villagers is examined finally, the share of each vulnerability level criterion is measured on the satisfaction of the villagers and appropriate solutions.

Data and Methodology
 The research methodology is based on its descriptive-analytical nature. Data gathering was conducted through surveying, library and field method. A small portion of the sample includes 230 households from 18 rural in the Neyriz Township. Reliability of the questionnaire was calculated using Cronbach Alpha (alpha = 0.79). In the qualitative section for the implementation of grounded theory, an interview was conducted with 40 villagers.

Results and Discussion
By studying the vulnerability of rural nursing homes in the Township of Neyriz, the villages of the studied villages are in an unfavorable position in terms of economic and social dimensions. Objective satisfaction indicators show that 51.8% of the villagers' homes were constructed responsive to concrete. 80.4% of the walls of the houses are made with bricks. The roofs of the houses are covered with 75% block and block. Of the studied rural households, 94.6% have personal housing and only 5.4% of the tenants. In the area of providing services in residential units, all the studied villages have water, electricity, telephone and 2.05 villagers are satisfied with the crop of agricultural products, parking lots, agricultural machinery and heating and cooling equipment for their housing.

Conclusion
 Findings showed that the vulnerability of rural housing is not only physical and environmental in nature, but also in social, economic and institutional-organizational dimensions. Meanwhile, vulnerability in physical and environmental dimensions in the study area is lower than other dimensions. Therefore, the existence of unstable housing in rural areas has led to a decline in their quality of life. This situation is strongly influenced by internal and external factors and forces. The lack of financial support, the traditional housing structure, poor design, poor monitoring and enforcement, social constraints, lack of building facilities, weaknesses in government support and regulatory policies, and government institutions are among the main problems of rural dwellings.
According to interviews with villagers, the following solutions can be made to reduce rural housing problems:
 - The costs of facilities and infrastructure are not at the expense of villagers. Therefore, the creation of continuous financial resources for councils and departments can be open.
 - Reducing the problems and obstacles facing applicants for loan use (through increasing credit, decreasing profits and raising the age), providing welfare services and reducing the total poverty of the rural community, granting loans or with benefits and installments Low for women-headed households; increasing the number of borrowers and creating rural people's interest and motivation for living in the countryside.
 - Preserving indigenous architecture, using indigenous materials, avoiding blind imitation of urban housing, etc., are unfortunately much neglected, and new rural houses have become homogeneous and adapted to the natural and physical environment of the countryside.
- In anti-poverty programs, the problems of villagers have been underestimated, which has led to their vulnerability. Therefore, investing in villages, creating complementary agricultural businesses, increasing production and, consequently, increasing rural incomes, can accelerate the growth and development of this sector.Energy saving is considered to be a problem with rural housing problems in terms of access to fossil fuels and mechanical equipment for heating and cooling buildings.

 Key words: Vulnerability, pathology, Earthquake, Neyriz Township

 

---

---

Designing a Volunteer Geographic Information-based service for rapid earth quake damages estimation


Introduction
The advent of Web 2.0 enables the users to interact and prepare free unlimited real time data. This advantage leads us to exploit Volunteer Geographic Information (VGI) for real time crisis management. Traditional estimation methods for earthquake damages are expensive and time consuming. In contrast, volunteer and web-based service are near real time with almost no cost services. the lack of accessible real time data collection services causes delayed-emergency responses for disasters like an earthquake. This drawback is critical when we encounter a problem like buried people with valuable seconds for emergency rescue operation.
The current research aims to design and implement a web-based volunteer data collection service for rapid estimation of earthquake damages and number of buried people.

Methodology
To investigate the capacity of VGI in rapid estimation of earthquake, a technical frame work based on the web technologies has been programmed and implemented. The designed service is comprised of server and client sides.
The client side is a two-side browser-based service includes volunteers (users) and managers pages. On the user page, volunteers have a web page to enter and fill in the blank forms and taking a photograph of the target building and compare it with pictures. They watch the sample pictures in different level of damages and compare their building with the samples and give a grade of the most similar picture with their building. This grading process leads the server to analysis and classify the incoming data and create the heatmaps for managers. On the managers page two online discrete heatmaps for the both earthquake damages and buried people are displayed. In fact, the heatmaps present the online and real time quantitative situation of the building damages and buried persons as hot spots. These hotspots have the first priority for giving emergency services. The manager page also exploits query tools to request different level of details and classes from the server side.
The server side is responsible for receiving, saving, spatial analysis and transmission of the requested result to the client side. This task is carried out by the exchange side. As the citizens are entered to the browser-based service and fill in the blank forms for building damages based on the mentioned guideline and report the buried people, These forms are transmitted to the server side and a geo-server performs spatial analysis including Heatmap, distance and clustering analysis. Then, a real time damage and buried people map are prepared and delivered to the client side. The server updates the created maps whenever a new data is submitted. By this, a real time damage and buried people maps are accessible for official managers to conduct a goal-oriented emergency operation and a preliminary earth quake damages on city building blocks.
After the technical frame work has been designed, it was tested in Oshanvieh city by 132 volunteers on the scene for an earthquake.

Results and discussion
To investigate the capability of volunteer geographic information for earth quake afterwards, the designed service mentioned in the methodology was utilized on Oshnavieh city. It was assumed that an earthquake has occurred. 132 volunteers participated for the data collection process. According to the crisis management organization experts, 102 reports of the total 132 reports are correct that shows the accuracy of 76.52 percent. Besides the building damage level based on the defined guideline, the citizens also select their vital needs after the earthquake.
  the most requested vital needs are warm stuffs, medicine, water and foods respectively. Unfortunately, the participation rate is ranged from some seconds after the earthquake to three days. This means that some citizens have filled and transmitted their data three days after the earthquake.
In the following, a comparison between the designed service and traditional earthquake damage estimation methods (in situ) was carried out. The result shows that field-based methods for a city like Oshnavieh need about 20 days. However, the designed volunteer-based service what is programmed and implemented in the current research does this job by 3 days.

Conclusion
As the results show, the proposed service designed in this research implements the damage estimation process 6.5 times faster than the governmental procedures. This proves the efficiency of the research achievements. Besides the velocity, traditional damage estimation methods are expensive compare to volunteer-based data collection and processing which are almost free, scalable and pervasive.

Keywords: Volunteer Geographic Information (VGI), earthquake damage estimation, heatmap, oshnavieh city.



 

---

Introduction
Geographical situation of Iran is a place for interacting many physical and human processes which lead to specific precipitation climatology in the country. The month to month variation of precipitation is one of  the features which the precipitation climatology may reflect due to tempo - spatial characteristics. In fact, monthly distribution of precipitation is one of precipitation normal features building up the climate structure. In order to recognize this fundamental characteristic three following questions have been raised:
1) Have the month to month distribution of precipitation changed over recent four decades?
2) How is the pattern of relationship of month to month distribution of precipitation and spatio - topographical variables?
3) Is it possible to find a spatial pattern for decadal changes of precipitation of month to month distribution?

Data and Methods
In order to find a responses for the abovementioned questions the distribution of month to month precipitation and its decadal changes was considered by adopting coefficients of variations (CV) for 46 years (1970-2016)  and using the third version of Asfazari dataset. The relationship of precipitation data and spatio-topographical variables calculated based on regression techniques. Moreover, the spatial pattern considered by using cluster analysis.  The CV calculated as follow:

here ،،  are ith raw's and jth column's CV, standard deviation, and monthly mean, respectively.
CV and its relationships with spatio-topographical variables were calculated in two temporal scale, for whole the under investigation period (1970-2016) and in decadal period for four decades (1977-1986, 1987-1996, 1997-2006, 2007-2016).
Discussion
 The results of current study proved that the month to month different in precipitation amounts have had spatial variations, whilst the temporal trends is not statistically significant. In addition, the minimum, maximum, and consequently, the range of values also the averages have not experienced significantly changes. However, the region experiencing the same values of precipitation illustrated oscillatory behavior. Accordingly, the decadal variations have happened in different areas. Although the there have been statistically significant relationships between monthly CV and spatio - topographical factors, the correlations were low. Based on cluster analysis, we found 5 regions according to CV and its anomalies in compares with normal CV for all under investigation period. These regions generally follow the latitudes from 32 N toward northern latitudes, whilst the region in the south of 32 N generally follow the longitude patterns.
Results
Precipitation is known a chiastic and complicated climate element. One of chiastic behaviors which precipitation shows in its different time - scale behavior is its month to month distribution among a given year. In current research the decadal variation of  above-mentioned behavior among recent four decades and the variation of its relationships and the spatio - topographical features , as parts of climate structure of the country, have investigated in details. 
Our finding illustrated that the month to month different in precipitation amounts have had tempo - spatial variations, whilst the temporal long - term trends is not statistically significant. Moreover, the values of minimum, maximum, and consequently, the range of month to month CV also the decadal averages have not experienced significantly changes over four under study decades. However, the region experiencing the same values of precipitation depicted oscillatory behavior. consequently, the decadal variations have happened in different areas. Although there have been statistically significant relationships between monthly CV and spatio - topographical variables, the correlations were not considerably high. Based on cluster analysis technique, we found 5 regions according to CV and its anomalies in compares to normal CV for all under study decades. These regions generally follow the latitudes from 32 N toward northern latitudes, whilst the region in the south of 32 N generally follow the longitude patterns.

KeyWords: Iran precipitation, Month to month changes in precipitation, Inter annual variation of precipitation, Precipitation anomaly, Spatial analysis of precipitation

 

---

Streambank erosion hazard analysis by BEHI method, case study: Sajadroud stream, Mazandaran province
Extended Abstract
Introduction
River bank erosion is a complex natural process and plays an important role in the dynamic equilibrium of the river. The amount of river bank erosion affects the river plan, cross section and dimensions of the river and is the main factor controlling channel migration and the evolution of the river planform. This research was conducted with two main objectives, which are: 1. River bank erosion susceptibility analysis in incised rivers in mountainous areas, 2. Comparison of river bank erodibility with two methods original BEHI and modified BEHI.
Method
In this study, Bank Erosion Hazard Index (BEHI) and the modified BEHI method along a part of Sajadrood Stream in Mazandaran province have been investigated. The original BEHI (Rosgen, 1996, 2001) evaluates the river bank erosion field measurements. In this method, several parameters are measured, including bank height, Bankfull height, bank angle, root depth, root density, surface protection, bank material structure and stratification. From this parameters, the bank angle and height, root depth can be measured, but indicators such as plant root density and surface protection are visually estimated as a percentage. The score of each index varies from 1 to 10. The total scores of all indicators are classified into 6 groups: very low, low, medium, high, very high and extreme.
Newton and Drenten (2015) Based on the modified BEHI proposed a protocol for estimating the river bank erosion. In this protocol, the lengths of similar riverbanks are first defined as a uniform section of bank. These uniform sections can be identified due to differences in bank slope, differences in bank material and a break in vegetation. The characteristics of uniform sections are then evaluated in a pre-screening questionnaire consisting of six questions. If the answer to two or more questions is "yes", there is a high probability of erosion and the BEHI evaluation will be performed. Otherwise, the measurement will not be taken because the erosion is low or very low.
Result and discussion
The study area was divided into four reaches and 36 sites were surveyed. These river reaches have a deep bed (incised), high slope, low sinuosity, low width to depth ratio and predominant cobble sediments in the bed. They are in type A3 according to the Rosgen River classification. The height of the banks of the stream is high and its average is 2.6 meters with a minimum and a maximum of 0.4 to 9.7 meters. The average height of bankfull was 0.44 meters and varied from 0.15 to 0.85 meters. Hence, the ratio of bank height to bankfull height has been high. The root depth of plants was low and their average was 0.34 meters. Root density of plants was also low and averaged 5.7% The stream bank angle is measured from 31 to 90 degrees and an average of 51degree. The percentage of surface protection varied between 25 and 65% and averaged 42%. From a total of 36 sites, 8% are in the medium group, 39% in the high erosion category, 45% in the very high erosion group and 8% in the extreme erosion category.
Among the various variables, the β coefficients of the surface protection index and the ratio of bank height to bankfull height were -0.62 and 0.51, respectively. To evaluate the modified BEHI method, all reaches were first examined according to the pre-screening table and more than two "yes" answers were confirmed for each reach. In the modified BEHI method, 26 out of 36 sites were in the erosion group. But in the original BEHI method, 16 sites are in a very high class. The total score of the studied sites was evaluated by two BEHI methods with Pearson correlation coefficients, which obtained a coefficient of 0.21 and shows a relatively low correlation.
Conclusion
Sajadrood stream has high and steep banks due to the incision created in the channel bed. The deposition of large boulders at the toe of the streambank has caused its protection and the flood currents of bankfull are not able to carry this piece of rock. Under these conditions, calculating the ratio of bank height to bankfull height cannot indicate the erodibility of the streambank in the bankfull stage. Nevertheless, the erodible potential of the bank for larger flood currents is confirmed. Comparison of the original BEHI with the modified BEHI showed that the modified method has an overestimate than the original BEHI.

Key word: river bank erosion, BEHI, Sajadroud, Mazandaran
 

---

Application of multivariate techniques in-line with spatial regionalization of AOD over Iran

Introduction
Models, satellites and terrestrial datasets have been used to detect and characterize aerosol. Nontheless, micoscale classification using remote sensing parameters considers as a deficiency. Thus, regionalizion and modeling aerosol without regard to political boundaries or a specific stations over Iran  demonstrates the spatial distribution of simple AOD structures.
Materials and methods
 Present study attempted to simulate and detect homogeneous areaes of aerosol in Iran using AOD (areosol optical depth) datast at 550 nm across Iran. Among the eigen techniques, principal component analysis (PCA) is the most applicable and controversial classification applied as multivariate analysis approach. In the line of the target, PCA, S-Mode separate the AOD subgroups with similar correlations. In the mode, m time series apply to each n station or grid points as a variable in the analysis, which is the territory of the region or geographical area. Mathematically, if the input data column in the Z matrix is applied as mathematical variables and the Z matrix has n points in the time series and m is the time step, then in the Zs decomposition has 3654×9985. In addition, the scree test and North's rule were used to cut-off the principal components and to select the number of appropriate special vectors to be kept.
Results and Discussion
For the study purpose, 85 percentaile of loadings were used to determine AOD areas over Iran. Using the method, the spatial patterns of Iran's aerosolshave been divided into six subregions, which are the major centers affected by the AOD. These major AOD hotspots affect by AOD extermes that are originated from aerosol surrounding sources. So that, the geographical location of sources areas have caused the northeastern atmosphere of Iran to be influenced by severe storms originating from the Karakum Desert. The same is correct concerning the East and Southeast regions. While, the intensification and transfer of aerosol from the Sistan plain to the south is increased AOD load over southeast Iran. Moreover, this study revealed a set points associated with distinguishing spatial differences between the west-northwest and southwest regions as well as central region that have not addressed in previous studies because of focus on ground-based observations. Also, the method illustrated that formation of the identified regions are a function of the volume, growth, and spread of aerosol particles resulting from the source regions in the Middle East. Finally, the classification techniques converting dynamic phenomenon such as aerosol into simpler structures presented a interpretable understanding of the geographical distribution of the phenomenon.
Conclusion
The present study identified the spatial patterns of AOD hotspots into six distinct regions including northeast, west-northwest, southeast, southwest, central and eastern Iran affected by the aerosol as well as major centers or high gradient areas. In addition, the present study not only supported by previous studies, but also it  make sense a regionalization that was neglected by former studies, whileseperated  the boundaries of the AOD areas without considering  provincial boundaries. Overall, the classification techniques, PCA, simplified a dynamic phenomenon such as aerosol into a simpler and illustrated geographical and interpretable understanding of the spatial distribution of the phenomenon.

Keywords: Aerosol Optical Depth (AOD), Multivariate Techniques, Regionalization, Iran
 

---

Abstract:
Clean air is a necessity for human well-being and health. Air pollution is a major threat to humans and other organisms and is considered as one of the environmental challenges. Today, with the increase in air pollution, the need to know more about the causes of its occurrence has been raised. The various consequences of air pollution have made air quality monitoring and control inevitable in all societies at the forefront of environmental issues. In recent years, air pollutants have caused serious risks to human health and the environment. One of these pollutants Tropospheric ozone is the cause of health and environmental problems, especially respiratory problems and lung dysfunction and asthma attacks. Other effects of tropospheric ozone can be reduced lung capacity, cough, chest pain, sore throat, condition Nausea, damage to plants (growth disorders and the effect on germination) and reduced tire life, hence, it is necessary to know and study the tropospheric ozone in large and industrial cities. Tropospheric ozone is a pollutant because it plays an effective role in converting primary pollutants into secondary pollutants. Therefore, the aim of this study was to investigate the trend of changes in tropospheric ozone concentration with meteorological parameters, ozone precursors (nitrogen dioxide and nitrogen oxides). In this study, data from Kermanshah synoptic station and air quality station of Kermanshah General Department of Environmental Protection (Ziba Park station) in a long-term period of 10 years (2007-2016) have been used. Also, in this study using sensor images The OMI satellite Aura was surveyed in February and July 2016. Tropospheric ozone is known as a pollutant in Kermanshah. Therefore, no systematic studies have been conducted on the recognition of tropospheric ozone and the relationship between tropospheric ozone and meteorological parameters in Kermanshah over a long period of time. Tropospheric ozone and its relationship with changes in nitrogen oxides, nitrogen dioxide and synoptic parameters in Kermanshah were studied and the correlation between tropospheric ozone concentration and meteorological parameters was studied by Pearson test and the relationship between them was studied by linear regression. Based on the results. The maximum concentration of ozone occurs in the afternoon between 14:00 and 17:00 and the maximum amounts of nitrogen oxides occur at night and in the early morning of the year. Also, the study of seasonal changes in ozone concentration showed that in warm seasons due to the conditions of tropospheric ozone formation, including the intensity of sunlight, temperature and time of radiation and the presence of pollutants including nitrogen oxides, the concentration of tropospheric ozone was much higher. Ozone concentrations are highest in June, July, August, and spring and summer. The results also showed that there is a direct relationship between solar radiation and ozone concentration. Simultaneously with increasing solar radiation, it increases the air temperature, which increases the photochemical activity and thus increases the ozone concentration. This can be seen in the warm months of the year (June, July and August). Wind speed is also directly related to the concentration of tropospheric ozone. As the wind speed increases, the reactants mix faster and the tropospheric ozone concentration increases. However, precipitation is inversely related to the concentration of tropospheric ozone, which decreases with the occurrence of precipitation in the months associated with the onset of precipitation, and in the dry months of the year, the concentration is increasing. Therefore, meteorological factors and parameters play an important role in tropospheric ozone changes. Which can be seen by linear regression and Pearson test. The results of the study of nitrogen dioxide and nitrogen oxides showed that the highest concentration of nitrogen dioxide and nitrogen oxides during the day is the opposite of the concentration of tropospheric ozone and the lowest concentration of ozone occurs in summer due to increased solar radiation, increased oxidation of di Nitric oxide and nitrogen oxides, and as a result increase the concentration of tropospheric ozone, in autumn and winter, this amount has an increasing trend. Therefore, the trend of changes in tropospheric ozone concentration is the opposite of the concentration of nitrogen dioxide and nitrogen oxides, which can be seen in the daily, monthly and seasonal sections, which linear regression and Pearson test show this important and OMI sensor images confirm this fact. In conclusion of this study, all parameters related to the concentration of pollutants along with meteorological parameters have been effective factors in the concentration of tropospheric ozone.
Keywords: air pollution, meteorological parameters, ozone tropospheric, NO₂, NO_X
 

---

In many areas, snow cover in the mountains is a major source of surface and underground water supply. Due to climate change and its effect on the time of melting ,it  is very important for environmental planning to predict the arrival time of water from snow melt to water consumption cycle. The purpose of this study is to investigate the volumetric changes and time distribution of snow flood flows in future by integrating remote sensing , GIS and climatic models.The studied area is the Upper Basin of Amir Kabir Dam, which is located on the southern slopes of Alborz Mountains. In this study, digital elevation maps (DEM) and GIS software were used to estimate parameters such as area, environment, main length, highest and lowest elevation points. In order to complete the snow cover data, MODIS products (MOD10A100) were extracted and the snow cover was extracted in the Upper Basin of Amir Kabir Dam. Next, runoff and snow melting models were simulated using SRM software. Calibration and validation of the model's acceptable performance were estimated. Then, in order to investigate the effects of climate change on the future of snowmelt runoff production in the basin of Amir Kabir Dam, the latest CMIP5 climatic models were used under four scenarios RCP2.6, RCP4.5, RCP6.0 and RCP8.5. A survey on the relationship between snow cover area , temperature and precipitation was used to predict snow cover in the future. The increase in temperature in the autumn and winter season has led to a reduction in the shape of precipitation in the form of snow, and as a result, the amount of snow storm is reduced. The results show that the amount of runoff in the autumn and winter increases due to increased rainfall in the form of rain, and it will be  increased late winter and spring due to the increase in the amount of water resulting from snow melting. The results of this study are based on the increase of snow melt as a result of increased runoff volume, reduction of snow reserves and maximum flow transmission to earlier than normal conditions due to early snow melting due to temperature rise. Generally, in the future, the average annual runoff will be decreased about 1.1 cubic meters per second, and the average annual melting share will be about 13.9%

---

Introduction
Climate change is the greatest price society is paying for decades of environmental neglect. The impact of global warming is most visible in the rising threat of climate-related natural disasters. Globally, meteorological disasters more than doubled, from an average of forty-five events a year to almost 120 events a year (Vinod, 2017). Climate change refers to changes in the distributional properties of climate characteristics like temperature and precipitation that persist across decades (Field et al., 2014). Because precipitation is related to temperature, scientists often focus on changes in global temperature as an indicator of climate change. Valipour et al. (2021) reported the mean of monthly the global mean surface temperature (GMST) anomalies in 2000–2019 is 0.54 C higher than that in 1961–1990. Many studies have been done on climate change in Iran. These studies have mostly studied the mean and extreme temperature trends (Alijani et al., 2011; Masoudian and Darand, 2012). In general, the results of previous studies showed that the statistics of mean, maximum and minimum air temperature in most parts of the Iranian plateau have increased in recent decades. Also, the increase of minimum temperature is greater than maximum temperature.
A review of the research background shows that we need to understand more about regional climate change in Iran. Therefore, present study performs the climate change of 14 extreme temperature indices using multivariate statistical methods at the regional scale.

Data and methodology
Historical climate observations including daily maximum and minimum temperature were obtained from the Iranian Meteorology Organization for the period 1968 to 2017 at 39 stations. In this paper, 14 extreme temperature indices defined by ETCCDI were analyzed. The indices are as follows: (1) Annual maxima of daily maximum temperature (TXx); (2) Annual maxima of daily minimum temperature (TNx); (3) Annual minima of daily maximum temperature (TXn); (4) Annual minima of daily minimum temperature (TNn); (5) Cold nights (TN10p); (6) Cold days (TX10p); (7) Warm night (TN90p); (8) Warm day (TX90p); (9) Frost days (FD); (10) Icing days (ID); (11) Summer day (SU); (12) Tropical nights (TR); (13) The warm spell duration index (WSDI) and (14) the cold spell duration index (CSDI). The extreme temperature indices were extracted using R software environment, RclimDex extension. The Mann–Kendall Test and Sen’s Slope Method was employed to assess the trends in 14 extreme temperature indices. To identify homogeneous groups of stations with similar annual thermal regimes, Principal Component analysis (PCA) and Clustering (CL) was applied. Pearson correlation coefficient was used to investigate the relationship between height and trend slope.

Result
All the extreme temperature intensity indices (TXx, TNx, TXn, and TNn) showed increasing trends during 1968 to 2017. The increasing trends of TXx, TNx, TXn, and TNn were 0.2, 0.3, 0.44, and 0.5 ° C per decade, respectively. These results indicated that the extreme warm events increased and the extreme cold events decreased. The average of the extreme temperature frequency indices over Iran showed that the frequency of warm night (TN90p) and warm day (TX90p) significantly increased with a rate of 6.9 and 4.2 day per decade, respectively. Also, the frequency of cold night (TN10p) and cold day (TX10p) significantly fell with a decrease rate of 3.8 and 3.8 day per decade, respectively. The frequency of warm nights (TN90p) was higher than that of warm days (TX90p). The result indicated that the trend of nighttime extremes were stronger than those for daytime extremes. The average of frost days (FD) and icing days (ID) indices over Iran showed decreasing trends during 1968 to 2017 with rates of 3 and 1.1 d per decade, respectively. While, the averaged of summer days (SU) and tropical days (TR) indices over Iran showed increasing trends with rates of 4.4 and 6.4 day per decade, respectively. The warm spell duration index (WSDI) indices showed a clear increase, with a rate of 2.1 per decade. In contrast the cold spell duration index (CSDI) showed a significant decrease, with a rate of 1.7 per decade. In general, the cold indices displayed decreasing trends, whereas the warm indices displayed increasing trends over most of Iran. Pearson correlation coefficient between height and Sen’s Slope was estimated to be equal to -0.62 (p < 0.01). In general, the results of this study showed that there is a negative correlation between the elevation factor and the Sen’s Slope of warm extreme indices. That is, as the altitude decreases, the Sen’s Slope increases. Therefore, the stations located in low altitude have experienced stronger increasing trends than in high altitude. The area of ​​Iran was classified into four clusters using PCA and CL methods. Cluster 1 has experienced the strongest increasing trends. The average height of cluster 1 is 535 meters. Approximately 38% of the studied stations were located in cluster 1. Cluster 2 showed a moderate heating trends. 33% of the stations were located in cluster 2. Most of the stations of cluster 2 are located in the northwest and west of Iran. Cluster 3 showed a weak increasing trends compared to clusters 1 and 2. The stations of cluster 3 did not show a special geographical concentration and were scattered in all parts of Iran. 18% of the studied stations are located in cluster 3. The stations of Cluster 4, have experienced weak decreasing trends, which was different from the other three clusters

Conclusion
In this study we analyzed the climate change of extreme temperature indices in Iran. The result showed that the frequency of warm nights, warm days, summer days and tropical days increased. Also, the frequency of cold nights, cold days, Frost days and icing days decreased. The warm spell duration index showed a clear increase. In contrast the cold spell duration index showed a significant decrease. In general, the extreme warm events increased and the extreme cold events decreased over most of Iran. There is a negative correlation between the elevation factor and the Sen’s Slope of extreme warm indices (R = -0.62). Therefore, the stations located in low altitude have experienced stronger increasing trends than in high altitude. The area of ​​Iran was classified into four clusters using PCA and CL methods. Cluster 1 has experienced the strongest increasing trends. The average height of cluster 1 is 535 meters. Therefore, the most heating have occurred in Low-lying areas of Iran. Cluster 2 and Cluster 3 showed a moderate and weak heating trends, respectively. The stations of Cluster 4, have not experienced clear trends.

Key words: climate change; Extreme temperature; clustering; Iran



 

---

Objective: “The main and overall objective of this research was to identify natural hazards threatening Saluk National Park and Protected Area. In this regard, based on the assumptions of this research, in addition to the main objective, three sub-objectives were also determined for it. The sub-objectives include identifying hazards arising from climatic, geomorphological, and hydrological conditions.
Method: This research was conducted using field and library methods and descriptive analytical methods. The aim of the field studies was to identify the plant and animal species of this park and protected area, which was done using interviews with environmentalists, people living in villages adjacent to the park, and the results of censuses conducted by the Environment Organization. In library studies, concepts and definitions related to the environment and natural hazards were extracted from books available in libraries and documents of the Environment Organization. Statistics from meteorological stations in the region were used to study the climate of the region.
Results: Initial studies suggested that climatic, geomorphological, and hydrological hazards may be the most important natural hazards of the park. Based on the studies, climatic hazards including drought, frost, and frost have greater impacts in the study area. Geomorphological hazards of the park resulting from external dynamic processes including landslides, landslides, and erosion are ranked second, and hydrological hazards, especially river floods in the eastern part of the park in the Ruin River valley and the southern alluvial fans of the park, are ranked third.
Conclusions:
The most important conclusion of this research is that the natural hazards of this park naturally act as an ecological moderator. And these hazards can have detrimental effects on this park when they interact with negative human interventions, but if they are aligned with positive human actions, their effects will not have a significant destructive impact on this park.





 

---

Tehran metropolitan with its large population, daily migrant workforce and many students, needs to planning and designing watch/warning system to reduce the climatic problems for human health.for this purpose, we need to study the climate accurately and Since the factors affecting the climate of warm and cold periods in Iran are different, in this study , the meteorological variables of Tehran warm period (May to September 2002) turned into 4 components in Temporal Synoptic Index (TSI) using PCA Method and using P-Array and Varimax rotation.By the scores of components for each day, the clustering method (in ward method) were used and, the warm days of the year was divided into two cluster named favorable and oppressive airmasses. The average maximum air temperature that is more effective in mortality, was 36.13 ° C. Days with temperatures above  34 ° C, less pressure, mild winds , dryness and more sunshine resulted in more adverse weather conditions, which resulted in a 34% increasing in mortality compare with favorable weather. The total number of deaths from cardiovascular disease during the study period was 154046 that about  67%of deaths have been simultaneous with oppressive airmass.The epidemiological study of mortality also confirms the results of previous research in this area and shows that the incidence of mortality is higher in older people as well as in men. It is clear that not all mortality can be attributed to the effects of climate, but results show that change in climatic conditions will affect on mortality and also for study the effect of climatic hazards on human health, it is better that we study the effect of all variables together on humans.

---

Objective: Risk assessment in tourism destinations plays an important role in maintaining tourist safety and environmental sustainability. Identifying potential threats such as natural disasters, overcrowding, cultural damage, and social hazards provides the basis for more accurate planning and more effective tourism management. The present study was conducted with the aim of assessing environmental risks in the tourist areas of Dezful County.
Methods: The statistical population included tourists from rural areas of Dezful County, which, according to the statistics of the Cultural Heritage Organization, are visited by more than 10 million people annually. The sample size was determined as 384 people using the Krejci and Morgan table, and sampling was carried out randomly. In this study, 33 villages in four districts of the county were examined. The data were analyzed using a questionnaire and confirmatory factor analysis.
Results: According to the results of the study, only the Choghamish district faces high political risks, and the other districts have a moderate threat in this area. In the area of financial and economic risks, the matter is different, and only the central district faces a moderate threat, and the other districts face a high level of threat. However, the most important result of this study was in the area of social risks, where all districts showed a low level of threat, which indicates appropriate security and comfort for tourists traveling to this city. The results of the study also show a moderate level of environmental risks in the Shahyoun, Sardasht, and Choghamish districts, and even in the central district, these risks are in the low priority range. The results of the confirmatory factor analysis of tourism risk reduction strategies in Dezful city were classified into information and educational factors, legal and security factors, precautionary factors, and infrastructure factors.
Conclusions: The present study showed that social risks in Dezful county are at a low level and the security of tourists is ensured. However, economic and environmental risks in some sectors require more attention and management. Implementing the proposed solutions can help improve the quality of the tourism experience and the sustainable development of the region.