تحلیل بیلان آب زیرزمینی برای ارزیابی پایداری آبخوان دشت ایذه با روش IDW و تیسن

نوع مقاله : مقاله پژوهشی

نویسندگان

1 دانش‌آموخته کارشناسی ارشد، دانشکده علوم زمین، دانشگاه شیراز، شیراز، ایران.

2 دانش‌آموخته کارشناسی ارشد، گروه مهندسی منابع آب، دانشکده مهندسی آب و محیط زیست، دانشگاه شهید چمران اهواز، اهواز، ایران.

چکیده

این پژوهش با هدف ارزیابی کمی پایداری آبخوان دشت ایذه، رویکردی تلفیقی را اتخاذ کرد که در آن روش‌های درون‌یابی وزن‌دهی معکوس فاصله (IDW) و چندضلعی‌های تیسن با محاسبات بیلان هیدروژئولوژیک ترکیب شدند. داده‌های ماهانه سطح آب ۱۰ چاه در سامانه اطلاعات جغرافیایی (GIS) تحلیل شده و یک هیدروگراف نماینده استخراج شد. اگرچه مقایسه دو روش نشان‌دهنده همبستگی بسیار بالایی میان آن‌ها است، اما روش IDW با میانگین خطای مطلق 2.8 سانتی‌متر، دقت عملیاتی برتری در بازسازی تغییرات مکانی سطح آب از خود نشان داد. نتایج محاسبات بیلان، وضعیت بحرانی آبخوان را آشکار ساخت: ورودی‌های مؤثر سالانه حدود 18.9 میلیون مترمکعب در مقابل خروجی‌های 20.2 میلیون مترمکعب است. این عدم تعادل منجر به کسری خالص سالانه حدود 1.3 میلیون مترمکعب، افت متوسط سطح آب و ایجاد بیلان هیدرودینامیکی منفی شده است. در نهایت، یافته‌ها حاکی از آن است که بهره‌برداری فعلی فراتر از توان تغذیه طبیعی آبخوان رفته و پایداری بلندمدت آن را تهدید می‌کند. با این حال، تلفیق این روش‌های تحلیلی ساده، ابزاری کاربردی برای پایش و مدیریت پایدار منابع آب زیرزمینی فراهم می‌آورد. تحقق تعادل پایدار در دشت ایذه مستلزم اجرای فوری راه‌بردهای کاهش برداشت، اقدامات تعادل‌بخشی و بازنگری در الگوی کشت است.
 

کلیدواژه‌ها

موضوعات


عنوان مقاله [English]

Groundwater balance analysis for assessing the sustainability of the Izeh plain aquifer using IDW and Thiessen polygon methods

نویسندگان [English]

  • Mohammad Bagher Eskandari 1
  • Niloofar Rahimi 2
1 Master's Degree Graduate, Faculty of Earth Sciences, Shiraz University, Shiraz, Iran.
2 Master's Degree Graduate, Department of Water Resources Engineering, Faculty of Water and Environmental Engineering, Shahid Chamran University of Ahvaz, Ahvaz, Iran.
چکیده [English]

Abstract
Objective: This research, which quantitatively assessed the groundwater sustainability of the Izeh Plain, adopted an integrated approach that combined Inverse Distance Weighting (IDW) and Thiessen polygon methods with hydrogeological water balance calculations.
Method: Monthly groundwater-level data from 10 wells were analyzed in a Geographic Information System (GIS), and a representative hydrograph was extracted.
Results: Although comparison of the two methods indicates a very high correlation (with a determination coefficient exceeding 99%) between them, the IDW method demonstrated superior operational accuracy in reconstructing the spatial variations of groundwater levels, with a Mean Absolute Error (MAE) of 2.8 cm. The results of the water balance calculations revealed a critical condition of the aquifer: annual effective recharge amounts to approximately 18.9 million cubic meters (MCM) against outflows of 20.2 MCM. This imbalance leads to a net annual deficit of approximately 1.3 MCM, a decrease in the average groundwater level, and the creation of a negative hydrodynamic balance. Finally, the findings suggest that current exploitation exceeds the natural recharge capacity of the aquifer, threatening its long-term sustainability.
Conclusion: However, the integration of these simple analytical methods provides a practical tool for monitoring and sustainable management of groundwater resources. Achieving a sustainable balance in the Izeh Plain requires the immediate implementation of strategies to reduce abstraction, balancing measures, and a revision of cropping patterns. 
 
 
English Extended Abstract
Introduction
Despite extensive research on groundwater resources in Iran and numerous studies on interpolation methods, no study has yet conducted an integrated evaluation of the Izeh plain aquifer's stability by combining simple, practical interpolation techniques (Thiessen and IDW) within a GIS framework alongside groundwater budget calculations. Previous research in this domain has either relied on complex numerical and artificial intelligence models or remained confined to separate analyses; consequently, the synergistic combination of these straightforward interpolation methods with balance calculations to simultaneously analyze spatio-temporal variations of water levels and derive an accurate representative hydrograph has been overlooked. This research addresses this gap by presenting a systematic integration of Thiessen and IDW methods in GIS with groundwater budgeting to comprehensively evaluate the stability of the Izeh plain aquifer.
 
Method
Local methods, such as inverse distance weighting (IDW) and radial basis functions (RBF), perform prediction based on a subset of adjacent points and are classified as accurate (Exact Interpolators) in terms of values, meaning that predicted values at the sample locations are equal to the observed values.
In the inverse distance weighting method (IDW), which is known as a deterministic approach in spatial interpolation, the weight of each sample point in determining the value of unknown points is calculated according to the inverse distance from that point to the expected position.
In the applied part of this study, which is focused on the Izeh plain, the IDW method was used for interpolation of groundwater and rainfall data.
For spatial analysis of point data and providing continuous maps of hydrogeological parameters (such as surface water level), deterministic spatial interpolation methods were used. Spatial-temporal analysis of the hydrogeological variable in the Izeh plain using the IDW interpolation method based on monthly data was performed.
 
Results
The infiltration rate was 714 mm at the Izeh station, and with an aquifer area of 313 square meters, the effective infiltration volume was calculated.Applying an effective diffusivity of 20%, based on previous studies and the dominant lithology of the area (silty sandstone), this volume is estimated to be approximately 18.9 m³ per year.
Side currents were assessed to estimate lateral flow, initially by determining the hydro field potential using Surfer software. Subsequently, the hydraulic gradient (i) was calculated based on the difference in water level at the boundaries of the aquifer. The effective cross-sectional area of flow was obtained by multiplying the distance between streamlines by the effective thickness of the aquifer. Using a hydraulic conductivity coefficient (k) of 0.22 m³/day, extracted from previous studies, the discharge of lateral inflow was estimated using Darcy’s equation and based on the dimensions of the entrance boundary (750 m), resulting in 69 m³/day and an approximate annual volume of 295 m³.
Similarly, the lateral outflow rate was calculated as 76 m³/day, and the annual volume is close to 824 m³, considering the dimensions of the output boundary (450 m). In the remaining portion of the budget, the annual withdrawal from the aquifer, based on official statistics from the regional water organization for 10 wells, was estimated to be 200 m³.
 
Conclusion
The results of this research, with an integrated and operational approach to assess the stability of the Izeh plain and comparison of simple interpolation methods in hydrogeological analysis. The results clearly showed that integration of spatial - temporal variations analysis using deterministic methods such as inverse distance weighting and thiessen polygons, provides an effective and efficient tool for monitoring and understanding the dynamic situation of aquifers. This approach allows the drawing of more accurate representation of the hydrograph without requiring complex and expensive numerical modeling.
Quantitative evaluation of different components of the balance revealed the warning of the studied aquifer and showed a negative balance and a continuous drop in the water level. These conditions confirm the hydrodynamic instability of the reservoir and the increased pressure induced by the operation beyond renewable capacity. In general, the findings of this study emphasize the necessity of adopting immediate and sustainable management strategies to control the discharge and restore the balance of this vital resource. The proposed integrated approach, as an extensible analytical framework, can be used to monitor and evaluate the initial evaluation of other aquifers in the restricted area.

کلیدواژه‌ها [English]

  • Hydrogeological Balance
  • Natural Aquifer Recharge
  • Interpolation
  • Water Resources Management
[1] Gurdak JJ, Hanson RT, McMahon PB, Bruce BW, McCray JE, Thyne GD, Reedy RC. Climate variability controls on unsaturated water and chemical movement, High Plains aquifer, USA. Vadose Zone Journal. 2007; 6(3): 533-547. https://doi.org/10.2136/vzj2006.0087.
[2] Thomas BF, Famiglietti JS, Landerer FW, Wiese DN, Molotch NP, Argus DF. GRACE groundwater drought index: Evaluation of California Central Valley groundwater drought. Remote Sensing of Environment. 2017; 198: 384-392. https://doi.org/10.1016/j.rse.2017.06.026.
[3] Minnig M, Moeck C, Radny D, Schirmer M. Impact of urbanization on groundwater recharge rates in Dübendorf, Switzerland. Journal of Hydrology. 2018; 563: 1135-1146. https://doi.org/10.1016/j.jhydrol. 2017.09.058.
[4] Di Curzio D, Di Giovanni A, Lidori R, Montopoli M, Rusi S. Comparing rain gauge and weather RaDAR data in the estimation of the pluviometric inflow from the Apennine Ridge to the Adriatic Coast (Abruzzo Region, Central Italy). Hydrology. 2022; 9(12): 225. https://doi.org/10.3390/hydrology9120225.
[5] De Vera A, Alfaro P, Terra R. Operational implementation of satellite-rain gauge data merging for hydrological modeling. Water. 2021; 13(4): 533. https://doi.org/10.3390/w13040533.
[6] Ahmad AY, Saleh IA, Balakrishnan P, Al-Ghouti MA. Comparison GIS-Based interpolation methods for mapping groundwater quality in the state of Qatar. Groundwater for Sustainable Development. 2021; 13: 100573. https://doi.org/10.1016/j.gsd.2021.100573.
[7] Feng D, Chen Y, Jiang P, Ni J. Research on optimizing rainfall interpolation methods for distributed hydrological models in sparsely networked rainfall stations of watershed. Water. 2025; 17(22); 3: 237. https://doi.org/10.3390/w17223237.
[8] Antonakos A, Lambrakis N. Spatial interpolation for the distribution of groundwater level in an area of complex geology using widely available GIS tools. Environmental Processes. 2021; 8(3): 993-1026. https://doi.org/10.1007/s40710-021-00529-9.
[9] Kisiki CP, Ayenew T, Mjemah IC. Estimation of groundwater recharge variability using a GIS-based distributed water balance model in Makutupora basin, Tanzania. Heliyon. 2023; 9(4). https://doi.org/10.1016/j.heliyon.2023.e15117.
[10] Rajabi K, Nikbakht Shahbazi A, Fathian H, Zohrabi N. Izeh Plain Subsidence Modeling Using MODFLOW Mathematical Code. Iran-Water Resources Research. 2021; 16(4): 112-126. https://doi.org/20.1001.1.17352347.1399.16.4.8.8 [In Persian].
[11] Nadiri AA, Yousefzadeh S. A comparison of the performance of artificial neural network, fuzzy logic and adaptive neuro-fuzzy inference systems models in the estimation of aquifer hydraulic conductivity. A case study: Maraghe-Bonab Aquifer. Hydrogeomorphology. 2017; 4(10): 21-40 [In Persian].
[12] Morris DA, Johnson AI. Summary of hydrologic and physical properties of rock and soil materials, as analyzed by the hydrologic laboratory of the US Geological Survey, 1948-60. US Government Printing Office; 1967.
[13] Mehrabinejad A, Kalantari N, Alijani F, Mohammadi H. Zoning of denitrification with emphasis on hydrochemistry in the urban aquifer of Izeh, Northeast Khuzestan. Journal of Water and Sustainable Development. 2025; 12(1): 46-63. https://doi.org/10.22067/jwsd.v12i1.2401-1401 [In Persian].
[14] Taghizadeh-Mehrjardi R, Minasny B, Sarmadian F, Malone BP. Digital mapping of soil salinity in Ardakan region, central Iran. Geoderma. 2014; 213: 15-28. https://doi.org/10.1016/j.geoderma.2013.07.020  [In Persian].
[15] Croley II TE, Hartmann HC. Resolving Thiessen polygons. Journal of Hydrology. 1985; 76(3-4): 363-379. https://doi.org/10.1016/0022-1694(85)90143-X.
[16] Mushtaha AM, Van Camp M, Walraevens K. Quantification of recharge and runoff from rainfall using new GIS tool: Example of the Gaza Strip Aquifer. Water. 2019; 11(1): 84. https://doi.org/10.3390/w11010084.
[17] Bahtiar R, Wijayanto Y, Budiman SA, Saputra TW. Perbedaan karakteristik sebaran spasial hujan di kabupaten jember menggunakan metode inverse distance weighted (IDW) dan poligon thiessen. Jember, Indonesia: Berkala Ilmiah Pertanian. 2022; 5(1): 1–5.
دوره 2، شماره 1
فروردین 1405
صفحه 1-12
  • تاریخ دریافت: 02 بهمن 1404
  • تاریخ بازنگری: 20 اسفند 1404
  • تاریخ پذیرش: 02 فروردین 1405
  • تاریخ اولین انتشار: 08 فروردین 1405
  • تاریخ انتشار: 08 فروردین 1405