اقدامات مدیریت برداشت آب زیرزمینی با مدل ریاضی تفاضل محدود: مطالعه موردی آبخوان مرودشت استان فارس

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

نویسندگان

1 گروه مهندسی و مدیریت آب، دانشگاه تربیت مدرس، تهران، ایران

2 دانشیار، گروه مهندسی و مدیریت آب، دانشگاه تربیت مدرس، تهران، ایران

3 دانشکده مهندسی، دانشگاه بریتیش کلمبیای شمالی، پرینس جورج، کانادا

چکیده

در این مطالعه یک مدل ریاضی تفاضل محدود برای شبیه‌سازی جریان آب زیرزمینی در آبخوان یک محدوده مطالعاتی واقعی در استان فارس استفاده شد. در این روند، ابتدا در یک دوره شبیه‌سازی 15 ساله، تغییرات تراز آب زیرزمینی و بیلان آبخوان مورد تجزیه و تحلیل قرار گرفته و سپس مقادیر ضرایب آب برگشتی، ضریب نفوذ از بارش و ضریب تغذیه آب زیرزمینی واسنجی شد. براساس واسنجی مشخص شد که ضریب نفوذ از بارش برای کل محدوده، در حدود 23 درصد است اما در مقابل ضریب آب برگشتی در بخش‌های مختلف محدوده متفاوت است و مقدار بین 20 الی 50 درصد را به خود اختصاص داده است. سپس به جهت پیش‌بینی وضعیت منابع آب زیرزمینی برای یک دوره آینده، چندین اقدام مدیریتی مدنظر قرار گرفت. نتایج حاصل از شبیه‌سازی این اقدامات نشان داد که اقدام تغییر سیستم آبیاری از سنتی و سطحی به تحت فشار بیشترین تغییرات را نسبت به سایر اقدامات بر روی بیلان و تراز آبخوان در محدوده مطالعاتی داشته به طوری که در محل اعمال این اقدام، 17 متر افت تراز سطح آب زیرزمینی مشاهده شده و کسری مخزن 8/7 میلیون مترمکعب افزایش یافته است. دو اقدام تغییر الگوی کشت و کاهش برداشت از آبخوان تاثیر مثبتی را بر روی بیلان آبخوان و تراز آب زیرزمینی داشته‌اند، طوری که این دو اقدام به ترتیب سبب افزایش تراز آب زیرزمینی به میزان 4 و 7/11 متر در محل پیاده‌سازی شدند. برای مدیریت پایدار و استفاده بهینه از منابع آب بهتر است سیاست‌های تغییر الگوی کشت و کاهش برداشت از آبخوان به جای سیاست تغییر سیستم آبیاری مورد توجه قرار گیرد.

کلیدواژه‌ها

موضوعات


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

Groundwater withdrawal management measures using a finite difference mathematical model: A case study of Marvdasht aquifer in Fars province

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

  • Mohammad Ghaffari-Sola 1
  • Hamed Ketabchi 2
  • Davood Mahmoodzadeh 3
1 Department of Water Engineering and Management, Tarbiat Modares University, Tehran, Iran
2 Associate Professor, Department of Water Engineering and Management, Tarbiat Modares University, Tehran, Iran
3 School of Engineering, University of Northern British Columbia, British Columbia, Canada
چکیده [English]

In this study involves a finite difference mathematical model was used to simulate the groundwater flow in an aquifer within a real study area in Fars Province. Firstly, a 15-year simulation period was conducted to analyze groundwater level fluctuations and the aquifer water balance. Subsequently, key parameters, including return flow coefficients, infiltration rate from precipitation, and groundwater recharge coefficient, were calibrated. Based on the calibration results, it was determined that the infiltration rate from precipitation across the entire study area was approximately 23%. However, the return flow coefficient varied across different parts of the region, ranging between 20% and 50%. To predict the future status of groundwater resources, several management measures were considered. The simulation results indicated that changing the irrigation system from traditional surface irrigation to a pressurized irrigation system had the most significant impact on the aquifer balance and groundwater levels compared to other scenarios. Specifically, in the area where this scenario was implemented, the groundwater level declined by 17 meters, and the aquifer storage deficit increased by 7.8 million cubic meters. In contrast, the crop pattern modification and groundwater withdrawal reduction scenarios had a positive effect on the aquifer balance and groundwater levels. These two scenarios led to an increase in groundwater levels by 4 meters and 11.7 meters, respectively, in the affected areas. Overall, for sustainable water management and optimal use of water resource, it is recommended that crop pattern modification and groundwater withdrawal reduction policies be prioritized over changing the irrigation system policies.

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

  • Marvdasht aquifer
  • Fars province
  • recharge estimation
  • finite difference method
  • groundwater modeling

EXTENDED ABSTRACT

 

Introduction

Groundwater resources are a vital source of freshwater worldwide, especially in arid and semi-arid regions. The over-extraction of groundwater can have significant consequences and pose a potential threat to human life and survival. Therefore, it is essential to manage the use of these resources in a way that ensures their sustainability and prevents damage to their balance. One effective approach for evaluating the sustainability and management of groundwater resources is through the estimation of groundwater balance and the analysis of water table fluctuations.This approach quantitatively determines the groundwater system's inflows, outflows, and water table fluctuations. Also, this approach provides managers/policymakers with the tools to make informed decisions and plan effectively for these resources. There are various methods for evaluating and estimating the groundwater balance. One such method is numerical modeling, which has been widely used in many studies as an effective and efficient tool for assessing the groundwater balance and water table fluctuations, demonstrating satisfactory performance. This study aimed to develop a numerical groundwater model to assess and estimate the changes in groundwater levels and their fluctuations and the aquifer balance in a real region within Fars Province.

Material & Method

This study evaluates and analyzes the water balance changes and groundwater table fluctuations of the Marvdasht aquifer in Fars Province using a finite-difference mathematical modeling approach. The model is designed as a two-dimensional, single-layer representation of the groundwater system. The simulation period spans 15 years, covering the hydrological years from 2004 to 2018. For the model, both calibration and validation processes were conducted to minimize the differences between simulated and observed water table levels at observation wells.Finally, after completing the calibration and validation processes, the model was employed to predict and evaluate groundwater table levels for the next 15 years.

Results

Data from 54 observation wells in the study area were used to detemine the aquifer’s hydrograph over a 15-year period. The results show that groundwater levels declined by an average of 20 meters during the simulation period. The model calibration was conducted using recharge parameters dependent on return flow and rainfall infiltration coefficients, hydraulic conductivity, transmissivity, inflow and outflow rates at the boundaries. The rainfall infiltration coefficient set to 23% for all grid cells and the return flow coefficient varying between 20% and 50% across different sections of the region. During the 15-year simulation period, 12 years were used for calibration and 3 years for validation, and the model was evaluated using RMSE, ME, and MAE metrics. The results of RMSE = 0.65 m, ME = -0.58 m, and MAE = 0.58 m were calculated for the calibration period and RMSE = 0.12 m, ME = 0.24 m, and MAE = 0.24 m were calculated for the validation period, indicating that the model successfully simulated and estimated groundwater levels. The simulation revealed a -46.7 million cubic meter change in aquifer storage, and of the management scenarios applied to predict the aquifer's future status, two scenarios had positive impacts, improving groundwater levels and water balance, while two scenarios had negative impacts, leading to further declines, with the maximum groundwater level decline recorded at 17 meters and the maximum rise at 11.7 meters across the study area. On a long-term average, the contributions to recharge were estimated as total recharge at 421.55 million cubic meters, rainfall infiltration at 62.75 million cubic meters, surface water infiltration at 145.54 million cubic meters, and return flow at 213.25 million cubic meters.

Conclusions

In this study, a finite difference mathematical simulation model was employed to simulate groundwater level fluctuations and the water balance of an actual aquifer study area. After calibrating and validating the model, four management scenarios were applied to predict and evaluate the aquifer's status. Among the scenarios, transitioning irrigation systems from traditional surface irrigation methods to pressurized irrigation systems reduced groundwater levels by 17 meters in the scenario-affected region, while reducing groundwater extraction by 50% compared to current conditions led to an increase in groundwater levels by 11.7 meters across the study area. These findings highlight the significant impact of improved irrigation practices and controlled water extraction on the sustainability of aquifer resources.

Author Contributions

Conceptualization, H.K., D.M. and M.G.; methodology, H.K.; software, M.G.; validation, H.K. and D.M.; formal analysis, H.K., D.M.; investigation, M.G.; writing—original draft preparation, M.G.; writing—review and editing, M.G., H.K. and D.M.; visualization, M.G.; supervision, H.K., D.M.; project administration, H.K., D.M.;

Data Availability Statement

Not applicable

Acknowledgements

The authors would like to thank all participants of the present study.

Ethical considerations

The authors avoided data fabrication, falsification, plagiarism, and misconduct.

Azadegan, B., Hanafi, H., Varavipour, M., & Kardan Moghaddam, H. (2024). Numerical simulation of Taybad plain aquifer using MODFLOW and prediction of management scenarios. Journal of Water and Irrigation Management, 14 (2), 343-35 . (In Persian).
Ahmadi Badejani, A., MalekMahmoudi, M., Rajayee, T., & Javadi Rad, M.. (2020). Numerical Simulation for Present a Relationship for Optimize Exploitation Groundwater (Case Study: Damaneh Aquifer). Iranian Journal of Irrigation & Drainage, 14(4 ), 1362-1375. (In Persian).
Abraham, M. K., Priyadarshini, & Manikannan, K. (2021). Numerical Modeling for Groundwater Recharge (pp. 165–177).
Anderson, M. P., Woessner, W. W., & Hunt, R. J. (2015). Applied Groundwater Modeling. In Applied Groundwater Modeling.
Aredo, M. R., Lohani, T. K., & Mohammed, A. K. (2023). Numerical groundwater modelling under changing water abstraction in Weyib watershed , Ethiopia Numerical groundwater modelling under changing water abstraction in Weyib watershed ,. Cogent Engineering, 10(2).
Dousti Rezaie, M. , Zeinalzadeh, K. , Besharat, S. and Amirataee, B. (2022). Effects of Management and Climate Scenarios on Groundwater Level Changes: Case Numerical Modeling Study in Salmas Plain Aquifer. Iranian Journal of Irrigation & Drainage, 16(2), 280-293. (In Persian).
Ghadimi, S., & Ketabchi, H. (2019). Possibility of cooperative management in groundwater resources using an evolutionary hydro-economic simulation-optimization model. Journal of Hydrology, 578, 124094.
Ghadimi,S. & Ketabchi, H. (2019). Impact assessment of different management strategies implementation on the aquifer using numerical simulation (Case study: Namdan aquifer, Fars province, Iran). Journal of Water and Soil Conservation, 25(6), 1-23. (In Persian).
Gumuła-Kawęcka, A., Jaworska-Szulc, B., Szymkiewicz, A., Gorczewska-Langner, W., Pruszkowska-Caceres, M., Angulo-Jaramillo, R., & Šimůnek, J. (2022). Estimation of groundwater recharge in a shallow sandy aquifer using unsaturated zone modeling and water table fluctuation method. Journal of Hydrology, 605.
Kakavand, Y., & Mozaffari, J. (2023). Simulation of Groundwater Level in Nahavand Plain Using MODFLOW. Irrigation and Water Engineering, 13(3), 305–319.
Ketabchi, Hamed and Mahmoodzadeh, Davood. (2018). Critique of estimating groundwater resources balance in Iran: Strategies and suggestions.Tehran, Iran: Center of the Research Islamic Consultative Assembly. (In Persian).
Khodaiari, M. , Hessari, B. , Ahmadi, H. and Mohammadpour, M. (2023). Evaluation of artificial recharge project efficiency for groundwater recovery with mathematical modeling (A case study on Firuraq of Khoy Plain in West Azerbaijan province). Hydrogeology, 8(1), 170-186. (In Persian).
Lozano Hernández, B. L., Marín Celestino, A. E., Martínez Cruz, D. A., Ramos Leal, J. A., Hernández Pérez, E., García Pazos, J., & Almanza Tovar, O. G. (2024). A Systematic Review of the Current State of Numerical Groundwater Modeling in American Countries: Challenges and Future Research. Hydrology, 11(11).
Mahmoodzadeh, D., & Karamouz, M. (2022). A hydroeconomic simulation-optimization framework to assess the cooperative game theory in coastal groundwater management. Journal of Water Resources Planning and Management, 148(1).
ParsaSadr, H. , Mohammadzadeh, H. and Nassery, H. R. (2016). Numerical simulating of Sabzevar Roudab aquifer and checking of influences of constructing Sabzevar Roudab dam on it. Journal of Water and Soil Conservation, 23(1), 119-135. (In Persian).
Raheem, A., Ahmad, I., Arshad, A., Liu, J., Rehman, Z. U., Shafeeque, M., Rahman, M. M., Saifullah, M., & Iqbal, U. (2024). Numerical Modeling of Groundwater Dynamics and Management Strategies for the Sustainable Groundwater Development in Water-Scarce Agricultural Region of Punjab, Pakistan. Water (Switzerland), 16(1).
Rezaee, S., Javadi, S. and Kardan moghaddam, H. (2020). Assessment of groundwater resources management solutions by finite element in numerical simulation. Hydrogeology, 5(2), 32-42. (In Persian).
Raja, O. , Parsinejad, M. and Tajrishi, M. (2022). Simulation of Groundwater Balance Using Integrated Surface and Groundwater SWAT-MODFLOW-NWT Model (Case Study: Mahabad Plain). Water and Soil, 36(1), 31-52.
Shittu, M. T., Usman, M. A., & Ariyo, S. O. (2024). Mathematical Modeling of Transient and Steady-State Groundwater Flow in a Confined Aquifer. Journal of Applied Sciences and Environmental Management, 28(3), 741–755.
Sheikha BagemGhaleh, S. , Babazadeh, H. , Rezaei, H. and Sarai Tabrizi, M. (2023). Numerical Modeling and Trend Analysis of Mahabad Aquifer Quantitative Status. Water and Soil Management and Modelling, 3(2), 1-17. (In Persian).
Soltani, F., Javadi, S., Roozbahani, A., Massah Bavani, A. R., & Lotfi, S. (2022). Simulation of surface water-groundwater interaction using MODFLOW-OWHM (case study: Shazand plain). Iranian Journal of Ecohydrology, 9(1), 199–210.
Tao, Z. C., Cui, Z. N., Yu, J. Q., & Khayatnezhad, M. (2022). Finite Difference Modelings of Groundwater Flow for Constructing Artificial Recharge Structures. Iranian Journal of Science and Technology - Transactions of Civil Engineering, 46(2), 1503–1514.
Tripura, S. I. N. (2016). Mass Balance With Hydraulic Head Using Groundwater Modeling Estimating Groundwater Volumetric Mass Balance With Hydraulic Head Using Groundwater Modeling System in Tripura , India. 1(August), 19–29.
Usman, M., Qamar, M. U., Becker, R., Zaman, M., Conrad, C., & Salim, S. (2020). Numerical modelling and remote sensing based approaches for investigating groundwater dynamics under changing land-use and climate in the agricultural region of Pakistan. Journal of Hydrology, 581(April 2019), 124408.
Yifru, B. A., Lee, S., Bak, S., Bae, J. H., Shin, H., & Lim, K. J. (2024). Estimating exploitable groundwater for agricultural use under environmental flow constraints using an integrated SWAT-MODFLOW model. Agricultural Water Management, 303(July), 109024.