Journal of Spatial Analysis

---

---

---

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.
 

---

Aging is one of the most prominent indicators of demographic decline that most modern societies experience. At this stage of demographic decline, alongside a decrease and stabilization of mortality rates, birth rates also sharply decline. The development of technology and the mechanization of tasks, the improvement of quality of life and health-related indicators, individual-centered lifestyles, and increased economic inflation are significant factors in this issue. Iran is also among the countries on the verge of entering the stage of demographic decline. However, the intensity of this trend varies in different regions of the country. This article examines and analyzes the state of aging in the neighborhoods of the metropolis of Tehran. This research falls into the category of applied research in terms of purpose and is descriptive-analytical in terms of method. The research is based on the census data from 2016 and utilizes spatial statistical analyses. The positive values of Moran's autocorrelation analysis for each of the indices: aging (0.664), old-age dependency ratio (0.644), youth ratio (0.653), aging ratio (0.664), and aging index (0.665) in the neighborhoods of Tehran indicate a clustered pattern. This means that the issue of aging is more acute in some neighborhoods and areas of Tehran. Accordingly, the density of the elderly population is higher in most neighborhoods of the central and northern parts of the city. The final result shows that the distribution of the elderly space follows the logic of the social macro-ecology of Tehran and is relatively consistent with its natural-social topography. Furthermore, the spatial analysis of aging in the neighborhoods of this city shows that although all neighborhoods generally grapple with the issue of aging, planning and management should be based on the patterns and nature of the spatial distribution of this issue.
 

---

Deforestation is one of the primary challenges and environmental threats facing forest ecosystems, including the Hyrcanian forests, and occurs under the influence of various natural and anthropogenic drivers. This study aimed to model the probability of deforestation occurrence within the Loveh forest management district located in northern Iran. The dataset comprised 104 documented deforestation points and 14 explanatory variables, derived through spatial analysis using GIS and environmental, topographic, and anthropogenic data. To assess the relationships among variables and predict the likelihood of deforestation, two statistical models were employed: logistic regression and the Generalized Additive Model (GAM). The results revealed that the GAM outperformed the logistic regression model, achieving a higher Kappa coefficient (0.84) and Area Under the Curve (AUC) value (0.956), and providing a more realistic spatial distribution of deforestation risk. The most influential variables included distance from roads, slope, wind effect, and elevation. Based on the GAM output, approximately 20% of the study area was categorized as high and very high risk. These findings underscore the pivotal role of access infrastructure, human pressure, and climatic factors in accelerating deforestation processes. The results of this study can serve as a scientific basis for prioritizing conservation interventions, reassessing road development policies, and enhancing spatial planning for sustainable forest management in northern Iran.
 

---

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.

---

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.

---

Waste management is one of the main challenges faced by modern cities. Given the population growth and the increasing generation of waste, there is a growing need for innovative and intelligent methods in this field. Smart growth indicators can serve as tools to improve urban waste management. A waste management system comprises a set of activities aimed at organizing community waste through engineering and sanitary approaches. One of the most significant problems of coastal areas is the lack of proper waste management. Smart growth in waste management focuses on integrating technology and sustainable practices to optimize waste collection, reduce environmental impacts, and promote recycling. This study presents key indicators and trends related to smart waste management. The research employs a mixed-methods approach, combining quantitative and qualitative data via a descriptive survey. The study collected opinions from 20 experts in waste management and urban growth issues, as well as from randomly selected residents of Mahmoudabad city. Data analysis was conducted using grounded theory for qualitative data and structural equation modeling for quantitative data. The results indicate that the smart growth indicator of modern leadership, with a mean score of 4.6, and adequate infrastructure, with a mean score of 4.04, hold the highest average values among the smart growth indicators affecting waste management in the coastal city of Mahmoudabad.
 

---

Drought is one of the most significant natural hazards, characterized by complex spatiotemporal patterns. This study analyzes the structural and spatial characteristics of droughts in Iran across three temporal scales: annual, seasonal, and monthly. To achieve this, the intensity and extent of droughts were calculated using the RAI index and ERA5 monthly precipitation data over 42 years (1979–2021). Initially, the spatial distribution and directional trends of the drought centroid were examined, and its spatial variations over the years were analyzed. Additionally, the relationship between the location of the drought centroid and its extent was investigated. The results revealed that during the cold season, the drought centroid is primarily concentrated in central Iran, while in the warm season, it shifts toward the northwest, the Caspian Sea coast, and the southeastern regions of the country. The distribution pattern of droughts at all scales predominantly follows a northwest-to-southeast trajectory. Furthermore, shifts in the drought centroid toward the northeast, east, southeast, and south were observed to coincide with an increase in drought extent, whereas shifts toward the north, northwest, and west were associated with a reduction in drought extent. Overall, the findings of this study demonstrate a direct relationship between the location of the drought centroid and changes in drought extent, despite the fact that droughts in Iran lack consistent and predictable spatiotemporal patterns
 

---

subtropical high pressure (STHP) and Mediterranean cyclone are among the most important synoptic systems affecting Iran's climate. In this study, the effect of the high altitude location of the sthp on the Mediterranean gyres during the droughts and wetness of Iran during 1979 to 2020 was analyzed. In this regard, two datasets were used. Station data were used to identify drought and wetness periods, and ECMWF-ERA5 grid data was used to identify the location of high pressure in the subtropical region. The results showed that STHP with 3 anticyclone cells (ridge) affects the position of atmospheric waves affecting Iran's rainfall. The STHP system, especially the Arabian Subtropical anticyclone (ASA) and North Africa, play a more important role in the location of the cyclone affecting Iran's rainfall, so that widespread droughts with the expansion of the ASA to the west and its integration with the African anticyclone, the lack of expansion of the Mediterranean trough to the sea Redness and reduction of Sudan low and Mediterranean integration systems occur. With the eastward movement of the ASA over the Arabian Sea and the northern Indian Ocean, the Mediterranean trough deepens and the amount of waves and consequently the rainfall of the country increases. Therefore, the eastward expansion of the Arabian Peninsula and the strengthening of the North African Ridge provide the conditions for the expansion of the Mediterranean Sea. Whenever the ASA is located in its easternmost position on the Oman Sea and the Arabian Sea, it will lead to the advection of moisture for Iran through the access to the large areas of southern water and eventually rainfall. The main cause of the occurrence of drought and wetness in Iran is the spatial variations of atmospheric waves due to the spatial variations in the ASA.

---

Abstract
Objective:The Kabul Basin is one of the most vulnerable regions in Afghanistan due to the frequency of heavy rainfalls and devastating floods. This study aims to identify heavy rainfall events (above 20 mm) and analyze their synoptic mechanisms, focusing on their causes and patterns.
Methods:The study uses a descriptive-analytical approach based on daily rainfall data from 18 hydrometeorological stations in the Kabul Basin over the statistical period of 2008 to 2022. Heavy and flooding rainfall events were identified using the environmental-circulation method. Cluster analysis was conducted using Ward’s hierarchical clustering technique, and GrADS software was employed to extract and interpret synoptic maps.
Results:The analysis revealed three main synoptic circulation patterns responsible for heavy rainfalls in the basin. Three representative days were selected for detailed analysis: March 23, 2009 (31 mm rainfall at Qala-e-Malak), March 17, 2014 (59 mm at Bagh-e-Umumi), and February 5, 2017 (60 mm at Qala-e-Malak). These events were associated with Mediterranean troughs, cold Siberian air intrusions, and Indian anticyclone influence, which collectively intensified rainfall. The findings suggest that these systems can be monitored in advance for early warning.
Conclusions:Heavy and flooding rainfalls in the Kabul Basin are strongly influenced by specific synoptic systems and atmospheric interactions. Recognizing these patterns enables early detection of risk and can improve the efficiency of disaster preparedness, water resource management, and regional warning systems. This study provides valuable insight for reducing vulnerabilities and mitigating the impacts of extreme weather events in the region.

---

Climate change stands out as one of the most pressing environmental challenges of the modern era, exerting profound impacts on aquatic ecosystems—particularly wetlands. This study investigates the influence of climate change on three wetlands in Chaldoran County, West Azerbaijan Province—Pir-Ahmadkandi, Naver, and Zavieh-ye Sofla—spanning the period from 1984 to 2023. To achieve this, climate data were obtained from the TerraClimate database and CMIP6 model outputs under four emission scenarios. Landsat and Sentinel-2 satellite imagery, along with JRC/GSW data, were processed to evaluate changes in wetland surface areas. Annual wetland extents were extracted and compared against climatic parameters (temperature, precipitation, actual evapotranspiration, and snow water equivalent) using time-series analysis, Pearson correlation, and multivariate regression. Additionally, the Delta Method was employed for downscaled climate data to project possible trends over the next 20 years.
The results indicate that rising temperatures and evapotranspiration constitute the primary drivers of wetland shrinkage. Pir-Ahmadkandi and Naver have lost over 27% and around 20% of their surface area, respectively, whereas Zavieh-ye Sofla exhibits an irregular, seasonal reduction due to human interventions and agricultural runoff. Projections suggest that wetland surfaces—especially in Pir-Ahmadkandi and Naver—will continue to decline, potentially exacerbating drought conditions, diminishing biodiversity, and reducing water quality. These findings underscore the necessity of implementing sustainable water resource policies, controlling evaporation, and incorporating human impact assessments into conservation measures. Moreover, harnessing advanced hydrological modeling techniques and integrating remote sensing data with machine learning approaches may offer more effective strategies for safeguarding these vital wetland ecosystems.
 

---

Objective: During a landscape, it is not facile to discriminate land parts that have dissimilar amounts and types of vegetation. Plant Ecological Units (PEUs) are known as management units and are a reflection of the management actions and natural disturbances in the region. This research aims to fuse different resolutions of satellite images to increase the PEUs classification accuracy.
Methods: For this purpose, the Marjan-Borujen watershed in Chaharmahal va Bakhtiari province was selected. After field monitoring and surveys, four dominant PEUs groups were identified in the study area. In this study, bands from the Landsat_8 satellite images with 30 m  spatial resolution (bands 7_2) and a 15 m panchromatic band (band 8) were used, as well as the Sentinel_2 satellite images including panchromatic bands (8, 4, 2, 3) with 10 m spatial resolution. First step, using the Landsat panchromatic band, the 30-m bands were upgraded to 15 m through the pen-sharpening process; so the 15 m  data set was prepared from the Landsat_8 satellite. Then, to increase the spatial resolution of the 15-meter data set to 10 m, the Sentinel_2 panchromatic bands were used. In this way, the Sentinel_2 panchromatic bands were geometrically matched with the Landsat_8 15 m data set, and the Co-Registration process was performed with the minimum RMSE(0.05). Finally,  two data sets (2 to 8 bands) of the Landsat_8 satellite images with 15 m and 10 m spatial resolution, the PEUs classification maps were prepared using the RF classification algorithm, and the maps' accuracy was displayed as an error matrix.
Results: The results show that increasing the spatial resolution significantly enhances the accuracy of PEUs classification maps. The 15 m set shows an overall classification map accuracy of 66%, while increasing the spatial resolution to 10 m enhances the overall accuracy to 82%. As well as, the error matrix results show that the classification map procured from the 10 m set, all four PEUs groups have improved the producer accuracy, user accuracy, and kappa agreement index. So, in this map, PEU 2 and PEU 3 have the highest kappa agreement coefficient (83 percent).
Conclusions: This study shows that using the Gram-Schmidt fusion algorithm and consequently increasing the spatial resolution of Landsat 8 images from 30 m to 10 m reduces mixed pixels and increases pure pixels, which in turn improves the quality of PEU classification maps.
 

---

Objective: Land use/land cover (LULC) changes, as one of the main anthropogenic drivers, significantly influence runoff patterns and intensify flood hazards. This study aims to assess the impact of land use changes on flood hazard zonation over the period 2015 to 2024 in the Samian watershed, located in Ardabil Province, Iran.
Methodology: Satellite imagery from Landsat 7, Landsat 8, and Sentinel-2 was utilized to extract land use maps for the years 2015 and 2024 using the Google Earth Engine platform. LULC classification was performed using the Classification and Regression Trees (CART) algorithm. Subsequently, the Modified Flash Flood Potential Index (MFFPI) model was applied by integrating key environmental layers, including slope, flow accumulation, land use, geology, curvature, and soil texture, within the ArcMap environment to generate flood hazard zonation maps.
Findings: The results indicated substantial LULC changes between 2015 and 2024, including an 18.47% increase in irrigated agricultural lands, a 9.38% increase in residential areas, and a 25.85% rise in sparse rangelands. In contrast, dry farming lands decreased by 25.21%, dense rangelands by 9.14%, and snow-covered areas by 98.61%. These changes have led to a notable expansion of high-risk flood zones. The LULC classification achieved a high overall accuracy and Kappa coefficient exceeding 0.98, indicating reliable results.
Conclusion: The expansion of impervious surfaces and reduction in natural vegetation cover have increased surface runoff and, consequently, the extent of high-risk flood-prone areas. The MFFPI model, by incorporating both environmental and anthropogenic factors, proved to be an effective tool for flood hazard prediction and management.
 

---

Objective: This study aims to analyze greenhouse gas variations across Iran and to identify the gases that exert the greatest influence on their overall dynamics. The findings enhance understanding of atmospheric pollution patterns and support the development of effective mitigation strategies. These results provide a scientific basis for climate-change mitigation planning in Iran. The study relies on satellite-based remote sensing datasets.
Methods: This study analyzes the temporal and spatial variations of major greenhouse gases including carbon monoxide, nitrogen dioxide, ozone, water vapor, and methane across Iran from 2019 to 2024. Sentinel-5P satellite data were extracted via the Google Earth Engine platform, and after filtering and removing low-quality observations, the data were standardized using the Z-Score method to enhance comparability and correlation analysis. Principal Component Analysis (PCA) was applied to reduce data dimensionality and identify dominant variation patterns. Temporal and spatial trends were then quantified using complementary statistical techniques.
Results:
Methane exhibited a consistent increasing trend from late 2021 through 2024 and accounted for the largest share of total variance (R² = 0.87), likely reflecting intensified anthropogenic activities and regional climatic shifts. CO, NO₂, and O₃ were mainly affected by seasonal fluctuations and nonlinear factors, and no clear long-term increasing or decreasing trends were observed. Water vapor showed a direct relationship with temperature variations, water sources, and atmospheric patterns, with its lowest concentrations recorded during the cold months and increases observed in the warm months. PCA analysis indicated that the first two principal components explained more than 70% of the total data variance, with CH₄, O₃, and NO₂ contributing the most to the overall variations.
Conclusions: The study results indicated that greenhouse gas variations in Iran are simultaneously influenced by natural factors and human activities. The combination of satellite data, statistical analysis, and PCA enabled a precise assessment of the temporal and spatial trends of greenhouse gases, providing valuable information for planning pollutant reduction and developing strategies to combat climate change.



 

---

Objective: in recent decades, population growth, urbanization development, and change in land use have led flooding as one of the most destructive natural disasters in the world. Therefore, our goal is to identify flood areas and the synoptic patterns that lead to it, which are among the most important issues in preventing and reducing the effects of flooding and dealing with it.
Methods: In this study, in order to prepare a map of flooded areas, the extent of the floodwater that occurred in June (2024) in Ardabil province, were processed SAR radar images before and after the flood. Then, to identify synoptic patterns, daily maps of geopotential height at 500 hectopascals, sea level pressure at 1000 hectopascals, omega pressure at 500 hectopascals, and relative humidity at 700 hectopascals with a spatial resolution of 2.5 degrees in 2.5 degrees latitude were received and analyzed from the National Center for Environmental Prediction and the National Center for Atmospheric Research (NCEP/NCAR) of the United States.
Results: The flood area study indicated that in the studied province, Bilehsavar city with an area of 593 hectares, Parsabad city with 505 hectares, Meshkin-shahr with 245 hectares, and Germi city with 192 hectares were flooded due to the waterlog. The analysis of the flood zones also showed that the largest volume of flood entering Ardabil Province during the studied period was related to the northern cities of the province, where the provision of all moisture conditions and instability at the full depth of the troposphere layer led to the occurrence of heavy flood-causing rainfall in these areas.
Conclusions: The results of this study indicate that the use of radar data, due to its outstanding capabilities, is a useful tool in detecting and continuously monitoring of floods. Therefore, by detecting flood-prone areas and synoptic conditions that produce floods, executive managers can make the best decisions to deal with possible future floods.
 

---

This study investigates the spatio-temporal dynamics of evapotranspiration (ET) and its modulation by biophysical variables and land use/land cover (LULC) changes in the Karun River Basin, southwestern Iran, from 2000 to 2023. The basin, spanning 67,257 km² and characterized by diverse topography, experiences significant annual water loss (72% of 413 billion m³ national precipitation) due to ET, leading to salt and sediment accumulation. Data from MODIS products (MCD12Q1, MOD13A1, MCD43A3, MOD11A2, MOD16A3, CHIRPS) provided land cover, NDVI, albedo, LST, precipitation, and ET at 500-meter resolution, supplemented by Landsat imagery (30-meter resolution) for validation. Multiple regression and Geographically Weighted Regression (GWR) analyses revealed a 39.5% ET increase (31.48 to 43.92 mm/year), a 32.78% NDVI rise (0.18 to 0.239), and a 16.35% LST decrease (33.52°C to 28.05°C), correlated with a 6.90% agricultural decline (6,939,225 to 6,460,335 ha), a 6.94% rangeland increase (3,840,375 to 4,106,780 ha), and a 42.76% forest expansion (156,000 to 222,700 ha). GWR (AdjR² > 0.97, peak 0.9887 in 2010) identified spatial non-stationarity, with overprediction in mountainous northeast regions and underprediction in agricultural southwest plains, reflecting LULC influences. Landsat-derived false color composites and classifications (overall accuracy 85–90%, Kappa 0.85–0.90) validated a 2,477 km² forest loss to high-ET rangelands/agriculture, driving warm-season ET elevation. Results emphasize the need for integrated hydrological models incorporating irrigation data and high-resolution analyses to enhance sustainable water management in this water-stressed region.

---

                   

This study employs the Analytic Hierarchy Process (AHP) to assess the vulnerability and resilience of the urban water distribution network in the Feyzabad and Baharan districts of Sanandaj against the parameter of Peak Ground Velocity (PGV). The main objective is to identify the key factors influencing network vulnerability and to propose strategies for enhancing the resilience of this critical infrastructure. PGV values were derived based on data from 40 faults longer than 10 km within a 70 km radius of the city, using empirical attenuation relationships. Geological, geomorphological, soil type, and pipe diameter and material data were collected from reliable local sources.In the AHP model, the main criteria including PGV, geology, soil, pipe material, and pipe diameter were integrated with weights of 0.460, 0.112, 0.243, and 0.182, respectively, and vulnerability maps of the network were generated. Results showed that PGV values across the city range between 35 and 39 cm/s. In Feyzabad, lower PGV values combined with thick steel pipes and Quaternary alluvial soils resulted in 81% of the network falling into the low-vulnerability class and only 2.1% into the high-vulnerability class. Conversely, in Baharan, higher PGV values (39 cm/s), combined with small-diameter asbestos pipes and shale bedrock, placed 34% of the network in the very high-vulnerability class.