Journal of Spatial Analysis

---

Gully erosion is one of the advanced forms of soil erosion, which needs to be analyzed and identified in order to protect the soil. In this research, according to the complex system of factors influencing the creation of ditch erosion, 23 factors were analyzed in the two famous Dempster-Schiffer models and the entropy model, and using Google Earth images and field visits, 331 ditch points were identified, recorded, and a ditch distribution map was prepared. Spatial data of gully erosion distribution were divided into two random training (70%) and experimental (30%) groups. In this research, two indicators of tolerance coefficient and variance inflation factor were used to check the collinearity test, and as a result, two indicators of waterway density and relative humidity index were removed and 21 factors were used in the modeling process. The output results of the layers, weighting and classification and integration in two Dempster-Schiffer and entropy models are the extraction of the zoning map of the gully's erodibility sensitivity. and 30% of the calibration and validation of the models, the area under the ROC system performance characteristic curve and the area under the AUC diagram of the Dempster-Schiffer model with an explanation factor of 0.934 and the maximum entropy model with an explanation factor of 0.936, both models have an acceptable percentage of the area under the curve were that this issue shows the high performance of both models in the region. Among other results of statistical analysis, the prioritization of the impact of 21 factors in causing ditch erosion in the region was determined. The scientific results of the research can be promoted and taught, and from the practical point of view, the relevant executive body to control ditch erosion can take the necessary measures using the results of this research.
 

---

Developing roads and constructing new highway are urban policices contributing to solve transportation problems in cities. These projects often being passed through urban fabrics, so it is nessessery to buy and demolish buildings from their owners including individuals or governmental and public institutions to imply the projects. However, acquiring land is not an easy task and completing these projects may hit with long-term delays. This paper aimed to analyze the impacts of delaying in constructing Shoosh Highway in Tehran. The investigated impacts originate from land acquisition problems. The research data was obtained from many sources including documents and research reports, a survey, and interviews with Tehran Municipality managers. The One Sample T-Test in SPSS software was performed to analyze the data obtained form the survey. Findings indicate that the residents are often dissatisfied with the project because since the beginning of the project, social security decreased a lot and people are less likely to respect citizens' rights than before, recreational sites are often demolished, the value of residential buildings slowed down significantly, living costs incresed, and businesses were stagnant. In conclusion, if urban highways are not contributing to proper planning and site selection, they will disrupt the physical, social, and economic structures of urban neighborhoods and cause to many environmental problemes including air pollution. To avoid these adverse outcomes, it should be thought in advance about sufficient financial resources and possible practical methods to acquire land for projects. These consequences are reduced by studying and managing the risk of projects.
 

---

Abstract
Dust storm is a complex process that it was affected by relation between earth-atmophere system and point of veiw climatologist and meteorologist that they assessing atmospheric and climatic change, in general of world veiw, monitoring from dust cover is a need structures.
The western region of Iran is the study area. The data used in this study are divided into two categories: ground-based observations in 27 synoptic stations extracted from the Iran’s Meteorological Organization during period (1998-2010) and satellite MODIS images during the first to fourth days of July 2008. Finally, the aim has analyzed using Arc GIS and ENVI softwares and NDDI index.
According to results, interpolated map for the number of dusty days during the study period over the western half of Iran showed that extent of case study have not a equal system aspect quantity of occuring from dust phenomenon and how is it trend. The number of dust days increase from north toward south and sites located in northen proprotion of studied area have experienced a lower dust events. While, maximum hotspots are occuring over southwestern sites such as: Ahvaz, Ilam, Boushehr and Shiraz. Therefore, principle offspring of dust input has been out of country boundaries and arrived far way area. On based resultes obtined on satellite images using NDDI index also idicate that maximun of intense cover dust is observed over Fars, Ilam, Boushehr and Ahvaz provinces on the first, second, thrid and forth of July. But, the lowest rate of index situated in extent far such as: Eastern Azarbayjan, Western Azarbayjan provinces. Thus, parts located on the north of the study area experienced less dusty days and the maximum dust core was located in the southwestern (mostly ahvaz). The long-term result was consistent with the use of NDDI index and the daily average of NDDI index in the whole study area indicated the hotspot areas (Ilam, Ahvaz, Omidieh, Bushehr and Shiraz) during the first to fourth days July 2008. However, in the region has reduced the level of dust cover when a wet and cloudy synoptic system pass over the central and northwestern parts of the study area.

---

In order to plan, manage and exploit water and land resources, awareness of the spatial variability of resources, as well as understanding the response behavior of the watershed in order to model physical processes, has an identical significant role. Due to its location in arid and semi-arid areas, special climatic and geomorphological conditions, Qom-Roud basin is prone to flash floods. Due to the lack of hydrometric and topographical data with high accuracy in the basin, the use of hydraulic models does not lead to accurate results of the hydraulic characteristics of floods. In such a situation, the methods based on the geomorphological features of the basin can be advantageous. In this article, Variable flood stage method (VFS) method is used, which combines the hydraulic characteristics of the river with the geomorphic characteristics of the basin in order to estimate the water depth in the river caused by floods with different return periods. The water depth was investigated for different return periods of two, five, ten, twenty-five, fifty and hundred years. In each period, the highest water depth was in the parts near the outlet and the lowest water depth was in the upstream parts of the river. The research illustrations there is a direct relationship between the depth of water and the area of the sub-basin. The results of this research can be used for basins without hydrometric and topographic statistics with high accuracy in order to estimate the peak speed and flood depth.
 

---

Dust storms affect different regions of the globe countries and Iran. The dust storm events were considered as one of the climatic components in arid and semi-arid regions and it is called one of the most important environmental problems of these regions. Therefore, in this research, the activity class of sand dunes were investigated based on the Lancaster index in Iranshahr city in Sistan and Baluchistan province. The aim of this study is to investigate the climatic different conditions on the mobility of sand dunes in Iranshahr city. Meteorological data for synoptic station of Iranshahr were collected from the Iranian Meteorological Organization for 15 years (2003 to 2018). After examining the wind speed velocity in Iranshahr, and extracting the seasonal and annual Wind rose diagrams, Dust Storm Index (DSI) was calculated.  Finally, the effects of possible changes in climatic elements on the mobility of sand dunes were predicted. The results showed a significant correlation between the Lancaster index and the amount of annual rainfall, wind and drought index. The results of the sensitivity analysis also showed that if the frequency of erosive winds and potential evaporation and transpiration increases to 30%, the activity of sand dunes in Iranshahr station increases by 25 and 23%, respectively, on the other hand, with a 30% increase in rainfall, the activity of sand dunes decreases by 30%.
 

---

The zoning of phenomena in order to recognize and evaluate the current situation as a basis for decision-making is one of the of scientific management essentials of natural resources. In the management of desert areas to prevent paying exorbitant costs and making incorrect decisions, strong and well-reasoned methods are needed to understand the current situation and the development process. Therefore, the main goal of the current research is to provide a suitable model for zoning the spatial distribution of desertification risk by using Moora, Aras scoring techniques and the use of geographic information system in Yazd-Khizrabad sub-basin. At first, after determining the working units by geomorphological method, the importance of indicators in each unit was obtained based on the Delphi method. Shannon's entropy method was used for the relative importance of indicators. Then the decision matrix was formed and after weighting, the potential of desertification intensity was estimated by calculating the desirability coefficient from Moora and Aras scoring method. The results of the analysis of both models showed that the land units of Mountain Agricultural Grounds (MAG) and Plan Agricultural Grounds (PAG) have the highest desertification potential, which covers 7135 ha (9.11%) of the entire study area and Quantitative value of desertification potential for the whole region based on the all of indicators was placed in the middle class (IV). It is suggested that the obtained results and ranking should be taken into consideration in the plans for controlling and reducing the effects of desertification and rehabilitating degraded lands.
 

---

The growth of settlements and the increase of human activities in the floodplains, especially the banks of rivers and flood-prone places, have increased the amount of capital caused by this risk. Therefore, it is very important to determine the extent of the watershed in order to increase risk reduction planning, preparedness and response and reopening of this risk. The present study uses the common pattern of the machine and the classification of Sentinel 2 images to produce land cover maps, in order to construct sandy areas and determine land issues affected by the flood of March 2018 in Aqqla city. Also, in order to check and increase the accuracy of the algorithms, three software indices of vegetation cover (NDVI), water areas (MNDWI) and built-up land (NDBI) were used using images. The different sets of setting of each algorithm were evaluated by cross-validation method in order to determine their effect on the accuracy of the results and prevent the optimistic acquisition of spatial correlation from the training and test samples. The results show that the combination of different indices in order to increase the overall accuracy of the algorithms and to produce land cover maps, the forest algorithm is used with an accuracy of 83.08% due to the use of the collection method of higher accuracy and generalizability than compared to. Other algorithms of support vector machine and neural network with accuracy of 79.11% and 75.44% of attention respectively. After determining the most accurate algorithm, the map of flood zones was produced using the forest algorithm in two classes of irrigated and non-irrigated lands, and the overall accuracy of the algorithm in the most optimal models and by combining vegetation indices (MNDWI) was 93.40%. Then, with overlapping maps of land cover and flood plains, the surface of built-up land, agricultural land and green space covered by flood was 4.2008 and 41.0772 square kilometers, respectively.
 

---

Natural biological crises are a constant threat to human societies; Therefore, it is considered important to be prepared to control and prevent such disasters. In dealing with such urban crises, addressing the category of non-pharmacological interventions against infectious diseases can be a way forward. Therefore, the aim of the current research is to find the main effective indicators in the epidemic and to explain the environmental components resulting from the interaction and correlation of the relevant measurable indicators at the locality level based on previous environmental studies.
 In this research, in order to obtain a community of opinion about effective indicators in epidemiological crises, a systematic search was first conducted using the keyword of pandemic resilience. For this purpose, several widely used databases such as Web of Science, Scopus and Elsevier were searched between 2013 and 2023. After checking the quality of the conducted researches, 42 indicators were selected. Further, in order to find out the underlying variables and identify the basic factors or criteria in order to explain the correlation pattern between the observed variables, the exploratory factor analysis method was used and SPSS version 26 software was used to analyze the data. Data from 118 Tabriz neighborhoods were used for exploratory factor analysis. After collecting the data, the process of change, standardization has been done to prepare the data and convert the raw data into percentage, growth rate, average and ratio. Then, the status of urban resilience components against epidemic diseases was calculated and normalized separately based on factor load for each neighborhood. Finally, the results were illustrated using Arc GIS software.
 The findings of the research indicate that 42 indicators affecting the spread of epidemic diseases at the level of localities can be re-categorized into 8 components of accumulation, diversity and design, social factors, density, economic factors, health infrastructure, environmental pollution and green spaces. Also, based on the findings, it can be said that the condition of most of the components in the neighborhoods of Tabriz city is at an average level.
 Based on the results, the effects of environmental factors on the transmission of Covid-19 are differentiated spatially. These components represent more than 82% of the changes in effective environmental factors.
 

---

Considering that there are different forms of mass movements in the ChamGardlan watershed, especially along communication routes, agricultural and residential areas, it is necessary to create refugee maps Therefore, it is impossible to examine the factors influencing its situation in order to prevent and control this phenomenon. Therefore, during the field visits, geographical, physiography, land use, vegetation cover, soil, climatology, geography and geomorphology maps were produced in the GIS environment. The method of this research has been accomplished base of distinction of the geomorphological units while using aerial photos and crossing  basis maps. Then, the effective factors on the occuration of mass movements were studied using logistic regression equations. So that, the factors such as slope, geological formation type, pedologic, climatic, etc. were taked into consideration as independant variables, and mass movements occurance frequency as function of mentioned factors.The result indicated that the effective factors related to frequency of land sliding happening in the area in arrangement, are slope, geological formation type and mass  movements  type (both kind and amount of salts in soil) and also landuse.
 

---

he spread of Covid-19 in the rural areas of the country has caused more dangers due to the common rural cultures, and the ignorance and lack of efficient management of this crisis in the villages has caused irreparable consequences for these areas. In such cases, the existence of social capital can be very vital in creating national consensus and successful policies to pass this critical stage. Leading research is applied in terms of purpose and based on descriptive-analytical nature. To determine the sample size of villagers, using Cochran's formula, from the total of 6903 households in sample villages, 362 households were calculated as sample households to complete the questionnaire. In order to investigate the effects of social capital on economic and social indicators that are effective in reducing the vulnerability of the outbreak of Covid-19 in the border villages of Sistan, a wide range of indicators was determined, and from the one-sample T-test and the analysis of variance of the regression model to measure the The evaluation of the effects of social capital on socio-economic indicators effective in reducing the vulnerability of the outbreak of Covid-19 in the border villages of Sistan was used. The results of the research showed that the higher the level of people's participation and their trust towards each other, the higher the level of responsibility and knowledge of people, it has a positive role and effect on social and economic indicators in order to reduce the vulnerability of the spread of the covid disease. has had 19
 

---

Subsidence as an environmental hazard is caused by various natural and human factors. The drastic changes in land use, the increase in the number of deep wells, and the effects of the subsidence phenomenon in Dehgolan plain show the need to investigate these influencing factors. In such a situation, adequate understanding of the degree of vulnerability and investigation of the influencing factors in that process provides the opportunity for planning and environmental preparation of the space in order to reduce vulnerability. In this research, first, the NDVI index of the plain was investigated with the help of 15 Sentinel-2 and Landsat 8 satellite images, and the best date was selected for the Sentinel-1 images. In this way, 8 Sentinel-1 satellite images were analyzed over a period of 8 years (2014-2021) and all the images were analyzed and processed in eight stages with the help of SNAP software. 3 Landsat 7 and 8 satellite images were used to investigate land use changes (2000-2021).By applying atmospheric and radiometric corrections and finally performing the supervised classification method using Arc GIS software, land use was extracted and its changes were checked. The interferometric results showed that the Dehgolan plain suffered a total of 480 mm of subsidence. So that 60 mm of subsidence has occurred in this plain every year. In the end, with the preparation of the map of land use changes, the classes of irrigated agricultural and residential lands increased by 6.98, 1.47 percent, and the uses of pasture, forest and rainfed lands were faced with a sharp decrease, so that irrigated lands increased by 8477 and residential by 672 hectares. Is. The results obtained from the analysis of the relationship between water use and subsidence showed that rapid subsidence occurs mainly in water and urban land use classes. This is a consequence of increasing water extraction for agriculture and drinking. Usually, the pattern of land use conversion with more human influences has increased the rate of subsidence.
 

---

Human health depends on living in a healthy environment. Various factors determine environmental health, which should be explored to be able to improve them. The purpose of this research is to analyze the factors affecting the health level of urban environment in district 20 of Tehran municipality and to find out their spatial inequality. This is a descriptive-analytical research based on documentary data of 6 indicators and survey data of 26 indicators. For the analysis, various tools of the GIS, especially the spatial analysis tool of Interpolation, are used. Quantitative analyzes are carried out by calculating statistical parameters in attribute tables. Findings show that the score of 2.29 out of 5 in the evaluation range of 1 to 5, indicates a poor situation of environmental health. Higher weakness belongs to the 6 indicators of the economic dimension including saving, sufficient income for buying cultural goods or going to recreational journeys, and residential home per capita. Access to urban parks and green spaces with scores of 4.55 and 3.43 respectively, show a superiority compared to other indicators. In the outskirts, the environmental health condition is much weaker than in the interior areas. Despite the vastness of worn-out fabric and informal constructions and abandoned spaces, what is more effective are the impacts of external factors, 3 sides of the district are limited to fringe lands, where the establishment of polluting industries, the passage of sewage canals and the replacement of foreign and the population involved in informal and illegal jobs are affected the health level of the urban environment.
 

---

In this research, by preparing a distribution map of landslide areas and assessing the potential of landslide vulnerability in the Sivan basin, it has been tried to identify the resources at risk in different zones of vulnerability, then the environmental stability of the basin against the potential Vulnerability due to landslides should be investigated.. For this purpose, the information layers of landslide areas, agricultural lands, forests, residential areas and roads were prepared in the GIS environment. Then, the prepared layers based on the model (Index Overlay) were weighted and classified using the Class Maps Multi method and using expert opinion.According to the results obtained from the research, 44 landslides were determined in the Sivan basin, 17 of which could be measured and investigated in the field, 10 of which were in mountainous areas and were inaccessible and 12 of which could not be reached at all. The arena was not identified and 5 cases have been leveled due to the change of use to gardens and facilities. In the following, the potential vulnerability map of Sivan basin landslides in five vulnerability classes: very low (19846 hectares), low (1793 hectares), medium (741 hectares), high (2089 hectares) and very high (871 hectares) was prepared. In terms of environmental stability, the most damages caused by landslides in the Sivan basin are respectively related to communication roads and forest areas, which are in the range of high and very high vulnerability, agricultural lands are in the range of damage. Medium and lowest acceptability is related to residential areas, which is in a very low range. In the investigation of the environmental stability of the basin, it has been determined that there is a potential for vulnerability in all the mentioned cases, but it will be more intense in the direction of communication roads.
 

---

---

Drought is one of the most important challenges faced by pastoralists around the world. This phenomenon can have significant negative effects on livestock health, production, and livelihoods. However, pastoralists can adapt to drought and reduce its negative effects by adopting various strategies. Semi-nomadic people in Darbandi, Kalat-Naderi County, have been facing drought since 2007 due to their livestock farming. Since livestock farming has profound impacts on the lifestyle and livelihoods of these communities, the present study examines their experience in facing drought and identifies their management strategies in these conditions. The research method is qualitative and the research tool is in-depth interviews with 20 semi-nomadic people in Darbandi, Kalat-Naderi. Sampling was purposeful and carried out until theoretical saturation was reached to ensure that a wide range of perspectives and experiences were collected. The data from the interviews were analyzed using a qualitative grounded theory approach to extract key patterns and concepts. According to the findings, the semi-nomadic Darbandi people of Kalat County have adopted various strategies in the face of drought, which are classified into four categories: rangeland and grazing management strategies, livestock nutrition management, water consumption management, and livelihood diversification. These results can be used as a basis for formulating better policies in the field of crisis management and rural development. Also, these results can be used for more effective planning to reduce the vulnerability of nomads to drought.

---

The occurrence of persistent rainfall, resulting from the integration of multi-scale cyclonic systems, is prone to producing heavy to severe precipitation. Therefore, it is highly significant due to its potential to cause disasters such as floods and landslides, as well as its importance in water resource management for agricultural purposes. In this study, persistent rainfall refers to rainfall events lasting at least three consecutive days with a total precipitation of more than 1 millimeter. The aim of this research is to identify and classify these types of rainfall for the western region of Iran over a 30-year statistical period (1993 to 2022) for the first time using the Lamb-Jenkinson method. In this method, the central coordinates of the study area are used as the reference point in the calculations. Sixteen additional points are also determined around the reference point. With the availability of instantaneous standard sea-level pressure data for these points, it becomes possible to calculate the values of geostrophic wind and vorticity. By comparing these two quantities, the types of weather patterns in the L-J method—which include four types: directional, cyclonic, hybrid, and undefined—are identified and categorized. Disregarding the undefined category, it was ultimately determined that the provinces of Kermanshah, Kurdistan, Hamedan, Lorestan, and Ilam each had 50, 50, 50, 40, and 39 occurrences of the directional state, respectively. Moreover, the frequencies of the cyclonic state for these provinces were 69, 94, 65, 66, and 38, respectively, with cyclonic rotation percentages of 100%, 98%, 97%, 95%, and 97%, respectively. As for the hybrid state, the frequencies obtained for each province were 49, 53, 43, 41, and 38, respectively.
 

---

The vulnerability of rural settlements to environmental hazards necessitates attention to local resilience. The aim of this study was to prioritize the factors affecting the physical-environmental and institutional resilience of Divshal rural district. The approach of the present study is descriptive-analytical and survey. Data collection was done in two forms: documentary and field (expert questionnaire). The statistical population of the study is 16 Rural managers of Divshal rural district. Using DEMATEL and ANP techniques, the obtained data were analyzed. According to the results of the ANP method, the sub-criteria of respecting the privacy of roads and preserving indigenous and cultural identity in new constructions have the greatest impact on the resilience of the study area. Also in the quality index of passages, sub-criteria of observing the privacy of passages; In the index of quality of access to services and infrastructure, quality of access to welfare and administrative services; In the shelter index, the existence of a strong public shelter in critical situations; In the index of quality of housing and constructions, preservation of indigenous and cultural identity in new constructions; In the environmental quality index, waste collection and disposal; In land use quality index, proper location of land uses; In the context of institutions, the existence of administrative organizations and institutions to help people; And in the Index of Institutional Relations, the cooperation of institutions in facilitating laws and giving credit to the people has the greatest effect on physical-environmental and institutional resilience.

---

Introduction: Assessing the vulnerability and pollution of the aquifer is necessary for the management, development and allocation of land use, quality monitoring, prevention and protection of groundwater pollution. The purpose of this research is to identify and analyze the qualitative vulnerability of the Mashhad plain aquifer in order to monitor and manage underground water resources and prevent its future pollution.
Methodology: Mashhad plain is located in the northeast of Iran between Binaloud and Hezarmasjed mountains and in the watershed of the Kasfroud river, and its area is 2527 square kilometers. In this research, the vulnerability of the Mashhad Plain aquifer was evaluated with DRASTIC and SI models, and ArcGIS was used to analyze the parameters and prepare the vulnerability map. DRASTIC model is one of the overlap and index methods. In this method, the seven measurable parameters for the hydrogeological system include the depth of the groundwater level(D), net recharge(R), aquifer environment(A), soil environment(S), topography(T), Impact of the unsaturated Zone(I) and hydraulic conductivity(C) is used. The ratings for the sub-layers of each criterion vary from one to ten depending on their impact on the vulnerability potential. In SI method, five parameters of groundwater depth(D), net recharge (R), aquifer lithology(A), topography(T) and landuse(LU) are used for aquifer vulnerability. After preparing the SI model layers and weighting each of the layer classes using the functions available in the ArcGIS, the sensitivity index is obtained from the weighted sum of the mentioned parameters.
Conclusion: Study area is divided into four zones with very low vulnerability(21.85%), low(32.09%), medium to low(31.05%) and medium to high vulnerability(14.59%). Also, based on the results of the SI model, the study area is divided into five areas with very low vulnerability(0.4%), low(24.63%), medium to low(23.98%), medium to high(18.71%) and high vulnerability(32.25%). In general, the vulnerability of the aquifer increases from the southeast to the northwest.For verification, statistical method and calculation of correlation coefficient between vulnerability maps and TDS layer was used in TerrSet software and the results showed that both DRASTIC and SI models have high accuracy in zoning the vulnerability of Mashhad plain aquifer, so that the correlation coefficient of vulnerability maps with index The quality of TDS in Drastic model is (0.996) and in SI model (0.995); Therefore, the results of the following research can be used in environmental assessments and analysis of various pollutions and can be used as a basis for management decisions.

---

Reducing the emissions caused by deforestation and forest degradation REDD is a strategy to moderate climate change, which is used to reduce the intensity of deforestation and greenhouse gas emissions in developing countries. In the last few decades, drastic changes in land use have caused a significant decrease in Hyrkan forests located in Mazandaran province. For this purpose, the aim of this study is to investigate the changes in land use and its prediction for the year 2050 using the Markov chain and the REDD project to reduce carbon dioxide emissions for the cities of Nowshahr and Chalus. Using the images of TM and ETM⁺ sensors of Landsat satellite, a land use map has been prepared in three time periods related to the years 1989, 2000 and 2021. Maximum likelihood method was used to classify images from supervised classification. From the error matrix, the Kappa coefficient in this evaluation was equal to 0.83 for 1989, 0.81 for 2000, and 0.92 for 2021. The results show that the forest cover decreases in 2050. In contrast, the area of ​​range land, city, barren land, agriculture and wetland will increase. Based on the goals of the REDD project, the amount of carbon dioxide emissions was calculated until 2050. If the REDD project is not implemented, a large area of ​​forest cover will be destroyed and a lot of carbon dioxide is released. The amount of carbon dioxide in the project area in 2021 is 49,681 tons and will reach 806,732 tons by 2051, and with the implementation of the REDD project in the region, this amount of gas can be increased to the equivalent of 402,321 tons. 404411 tons of carbon dioxide was prevented from entering the upper atmosphere of the earth. Examining changes using satellite images can help managers and planners to make more informed decisions.
 

---

Objective: Investigating the occurrence of land subsidence in the country and the extent to which rural settlements are exposed to the phenomenon of land subsidence.
Methods: The present study was based on library surveys and studies, field observations and impressions. Using information received from the National Mapping Organization, a map of the country's subsidence zones and the degree of risk of each zone, including very low, low, medium, high and very high risk, was drawn, and the aforementioned maps were compared with the location of the villages.
Results: According to the analysis of the available data, 302 villages are at very high risk, 768 villages are at high risk, 834 villages are in the medium risk zone, and 573 villages are in the low risk zone. In terms of percentage weight, about 4 percent of the country's villages are at medium to very high risk of subsidence, of which 1,904 villages are at medium to very high risk, and 573 villages are at low risk.
Conclusions: The highest provincial distribution of villages at risk of subsidence in the country with a very high degree is in the provinces of Alborz, Tehran, Khorasan Razavi, Qazvin, Kerman, Golestan and Hamedan, and the highest provincial distribution with a high degree is in the provinces of West Azerbaijan, Isfahan, Alborz, Tehran, Khorasan Razavi, Semnan, Qazvin, Kerman, Golestan, East Azerbaijan, Hamedan and Yazd. Also, the highest provincial distribution of villages at risk of medium-level subsidence is in the provinces of East and West Azerbaijan, Isfahan, Alborz, Tehran, Semnan, Qazvin, Kerman, Golestan, Mazandaran, Markazi, Hormozgan, Hamedan, and Yazd.