ارزیابی آسیب پذیری آبخوان کارستی با استفاده از مدلcop و pi (مطالعه موردی: آبخوان بیستون و پرآو)

منابع آب کارستی قابل شرب بخصوص در مناطق خشک و نیمه خشک از جمله ایران به عنوان منابع آبی ارزشمند و راهبردی مطرح هستند. آلودگی آبخوان‌ های کارستی تبعات جبران ناپذیری در زمینه‌ های توسعه و بهداشت سکونتگاه ‌های شهری و روستایی دارد بنابراین حفاظت از آنها از اهمیت زیادی برخوردار است. در این پژوهش برای ارزیابی آسیب‌پذیری آبخوان بیستون و پرآو که به لحاظ کارستی توسعه یافته و دارای درز و شکاف و لندفرم‌های متنوع است؛ از مدل‌ های آسیب ‌پذیری cop و pi برای شناسایی پراکندگی پهنه‌ های در معرض خطر آلودگی استفاده شده است. در این پژوهش از فاکتور‌های جریان (c)، لایه های فوقانی(o) و بارندگی(p)در مدل copو دو فاکتور پوشش حفاظتی (p) و شرایط نفوذ (i) در مدلpi و داده‌های کیفی آب زیرزمینی برای ارزیابی دقت مدلسازی مکانی آسیب‌ پذیری استفاده شده است. نتایج پژوهش نشان می‌ دهد که بیشتر آبخوان کارستی بیستون و پرآو در پهنه با آسیب پذیری خیلی زیاد و زیاد پراکنده هستند. همچنین نتایج ارزیابی دقت مدل‌ ها نشان داد که هرچند مدل pi در شناسایی مناطق با آسیب پذیری زیاد مطابق با داده های کیفی آب زیرزمینی است اما مدل cop در مدلسازی آسیب پذیری آبخوان کارستی موفق تر بوده است به طوری که پهنه های با آسیب پذیری بالا در شمال غرب، شرق و جنوب شرق منطقه با پهنه های آلودگی بیشتر آب زیرزمینی منطبق است.

vulnerability assessment of karst aquifer using cop and pi model (case study: bisotun and paraw aquifers)

introduction: drinking karst water resources, especially in arid and semi-arid regions, like iran, are considered as valuable and strategic water resources. a sharp decrease in rainfall reduces the quality and quantity of karst water sources (christensen et al., 2007). on the other hand, urban and industrial development, which is accompanied by the increase in population growth, increases the risk of underground water pollution caused by the dumping of chemicals, waste and change of use (mcdonald et al., 2011). protection of karst aquifer is one of the most important measures in the management of karst water resources due to its vulnerability and high sensitivity to pollution (khoshakhlagh et al., 2014, afrasiabian, 2007). therefore, with the advancement of geographic information system technology, rapid progress was made in the ability to identify and model groundwater pollution, as well as the vulnerability of water sources from these pollutants (babiker et al., 2004, rahman, 2008). the pollution potential decreases from the center to the periphery (saffari et al., 2021).materials and methods: in this study to evaluate the vulnerability of bisotun and paraw aquifer which is karstically developed and has, crack and fissure and various landforms; cop and pi vulnerability models have been used to identify areas at risk of contamination. the cop model includes three main factors including concentration of flow (c), overlaying layers (o) and precipitation (p). factor c, which indicates surface features (sf), slope and vegetation (sv). it was obtained between 0.8-0.0 in 5 classes. from the overlap of the subfactores soil, layer index and lithology, the o factor map was prepared in three classes, including class 2 with low protection value, 2-4 with medium protection value and 4-8 with high protection value.  the p factor, which is the temporal distribution of precipitation along with the intensity and duration of precipitation, can show the ability of precipitation to transfer pollutants from the surface to the underground water. p factor was 0.8 in 2 layers in the northwest of the study area and 0.8-0.9 with low protection value. furthermore, top soil, precipitation, net recharge, fracture density, bedrock and lithology maps were used for the protective cover factor (p) in the pi model. the zoning of the p factor showed 2 classes such as very low and low most of the study area is in the low class. the infiltration condition factor (i) using the characteristics of the soil, the slope layer, and the land use in four layers showed high, aamedium, low, very low, which due to the high slope of the area of ​​the high layer has the highest dispersion, which causes the reduction of the protective cover.results and discussion: consequently, cop vulnerability map in 5 classes with very high vulnerability (0-0.5) equal to 38774.74 hectares (41.4%) and very low vulnerability (4-9-4) with 57.86 hectares (0.06%) of the largest and smallest area respectively. also, the pi vulnerability map of the combination of these two factors showed very high vulnerability with the largest area of ​​about 68,783 hectares and 72.9% scattered throughout the study area and the high vulnerability class with an area of ​​about 25,526 hectares and 27%.conclusion: the results of this research showed that the simulation performance of each cop and pi vulnerability model is closely related to the amount of pollution in the environment. it seems that the cop vulnerability model can better and more accurately showed the level of vulnerability in the karst aquifers of bisotun and paraw.