Evaluating the Impact of Drought Stresses on Groundwater System in Bagh- Malek Plain by Discharge Pattern Changes

Document Type : Research Article

Authors

• R. Lalehzari ¹
• Saeid Boroomand Nasab ²
• Hadi Moazed ³
• A. Haghighi ¹

¹ Shahid Chamran University of Ahvaz

² Professor, Faculty of Water Sciences Engineering, Shahid Chamran University of Ahvaz, Iran

³ Shahid Chamran University of Ahvaz, Iran

https://doi.org/10.22067/jsw.v0i0.31087

Abstract

Introduction: Groundwater is the largest resource of water supplement and shortages of surface water supplies in drought conditions that requires an increase in groundwater discharge. Groundwater flow dependson the subsurface properties such as hydraulic gradient (water table gradient or head loss in artesian condition) and hydrodynamic coefficients. The flow treatment is analyzed with an accurate estimation of effective parameters in groundwater equation. This function is obtained using the continuous equation. Inlet and outlet flows of a cell are equal to storage amounts in the continuous equation. Analytical solution of this equation is complex, so numerical methods are developed including finite element and finite difference methods. For example, Feflow is a groundwater modeling tool that makesuse of finite element method (Reynolds and Marimuthu, 2007). Modflow as a finite difference three-dimensional model simulated underground flow under steady and unsteady conditions in anisotropic and non-homogeneous porous media. Modflow is designed to simulate aquifer systems in which saturated-flow conditions exist, Darcy’s Law applies, the density of groundwater is constant, and the principal directions of horizontal hydraulic conductivity or transmissivity do not vary within the system. In Modflow, an aquifer system is replaced by a discretized domain consisting of an array of nodes and the associated finite difference blocks. Groundwater modeling and water table prediction by this model have the acceptable results, because many different informations of water resource system are applied. Many people and organizations have contributed to the development of an effective groundwater monitoring system, as well as experimental and modeling studies (Lalehzari et al., 2013). The objective of this paper is investigation of hydraulic and physical conditions. So, a numerical model has to be developed by PMWIN software for Bagh-i Malek aquifer to calculate hydrodynamic coefficients and predict water table in the future.
Materials and Methods: Bagh-i Malek aquifer located in Khuzestan province is mainly recharged by inflow at the boundaries, precipitation, local rivers and return flows from domestic, industrial and agricultural sectors. The discharge from the aquifer is through water extraction from wells, springs, and qanats as well as groundwater outflow and evapotranspiration. In this study, conceptual model of Bagh-i Malek aquifer on the framework of finite difference numerical approach has been used in simulating groundwater flow treatment. Water table data of 8 piezometers was collected for the 10 year duration from 2002 to 2012. The study years are divided into 40 seasonal stress periods with daily time step. Hydraulic conductivity, specific yield and recharge were calibrated in these periods. Verification was made between the simulated and measured hydraulic heads in the next calibration year. To simulate the groundwater table elevation in this study area, the PMWIN model is used. Bagh-i Malek aquifer is considered as a single layered aquifer, and therefore only the horizontal hydraulic conductivity is estimated. Modflow was used to simulate both steady state and transient flow systems. In steady conditions it is assumed that the total of time simulation is a time period and it does not change inlet data in the modeling duration. In unsteady conditions,the duration of study is divided into some stress periods that data is changed in every period.
Results and Discussion: Estimation of hydraulic conductivity is the first step of calibration process at steady state conditions. The correct assignment of hydraulic conductivity has a main effect on other parameters accuracy. Hydraulic conductivity mapping indicated that the maximum values are in the Eastern North (6-7 m/day) of the aquifer. The twice calibrated parameter is specific yield in unsteady conditions. Specific yield mapping indicated that the values vary from 0.03 to 0.08 and are maximum in the Southern regions of the plain similar to hydraulic conductivity. The results confirm that the flow model has the tolerable simulation accuracy by variances of 3.1 and 3.84 in calibration and verification processes, respectively. The sensitivity of the flow model to decreasing the hydraulic conductivity is more than increasing it. 50 percentage declined into the hydraulic conductivity causes the increase of the variance from 3.1 of initial value to 44.
Conclusions: Mapping of calibrated hydraulic conductivity showed that the Eastern North of aquifer has the higher transmissivity and discharge capability in comparison to Southern parts. At last, the result show that the Bagh-i Malek aquifer model is sensitive to recharge, hydraulic conductivity and specific yield, respectively.

Keywords

• Calibration
• Finite difference
• Hydrodynamic coefficients
• simulation

20.1001.1.20084757.1394.29.3.3.0

References

1- Bazargan-Lari M.R., Kerachian R., and Mansoori A. 2009. A conflict-resolution model for the conjunctive use of surface and groundwater resources that considers water-quality issues: a case study. Environmental Management, 43:470–482.

2- Delavar M., Moghadasi M., and Morid S. 2011. A real-time model for optimal water allocation in irrigation systems during droughts. Journal of Irrigation and Drainage Engineering, ASCE. http://dx.doi.org/10.1061/(ASCE)IR.1943-4774.0000440.

3- Heidari S. 2012. Quantitative simulation of groundwater source of Qale-Tol plain in different scenarios for agricultural uses by PMWIN Model. M.Sc. Thesis of Irrigation and Drainage. Shahid Chamran University of Ahvaz. Pp:150. (in Persian with English abstract)

4- Huang J., Ridoutt B.G., Chang X., Zheng H., and Chen F. 2012. Cropping pattern modifications change water resource demands in the Beijing metropolitan area. Journal of Integrative Agriculture, 11(11):1914-1923.

5- Karamouz M., and Araghinejad S. 2008. Drought mitigation through long-term operation of reservoirs: case study. Journal of Irrigation and Drainage Engineering, 134: 471-478.

6- Karamouz M., Mohammad Rezapour Tabari M., and Kerachian R. 2007. Application of genetic algorithms and artificial neural networks in conjunctive use of surface and groundwater resources. Journal of Water International. 32(1):163-176.

7- Kresic N. 1997. Quantitative solutions in hydrogeology and ground water modeling. CRC Press LLC. Pp: 115.

8- Lalehzari R. 2008. The effect of wastewater recharge on nitrate distribution in Shahrekord aquifer using MT3D model. M.Sc. Thesis of Irrigation and Drainage, Shahrekord University Pp, 130. (in Persian with English abstract)

9- Lalehzari R. 2015. Optimal allocation of surface and groundwater resources to cropping pattern in Baghmalek plain by multi-objective planning based on non-dominated sorting algorithm. Ph.D. Thesis of Irrigation and Drainage, Shahid Chamran University of Ahvaz. Pp: 250. (in Persian with English abstract)

10- Lalehzari R., Tabatabaei S.H., and Kholghi M. 2013. Simulation of nitrate transport and wastewater seepage in groundwater flow system. International Journal of Environmental Science and Technology. 10:1367-1376. DOI 10.1007/s13762-013-0213-4.

11- Lalehzari R., and Tabatabaei S.H. 2014. Simulating the impact of subsurface dam construction on the change of nitrate distribution. Environmental Earth Science. Springer: DOI 10.1007/s12665-015-4362-2.

12- Mohammad Rezapour Tabari M., Maknoon R., and Ebadi T. 2008. Multi-objective optimal model for surface and groundwater conjunctive use management using SGAs and NSGA-II. Journal of Water and Wastewater, 1:1-12. (in Persian with English abstract)

13- Mohammad Rezapour Tabari M., Maknoon R., and Ebadi T. 2012. Development structure for optimal long-term planning in conjunctive use. Journal of Water and Wastewater. 4:56-69. (in Persian with English abstract)

14- Nejati Jahromi Z. 2009. Simulation of the groundwater sources in Aghili’s plain by finite differences mathematical model. M.Sc. Thesis of Hydrogeology, Shahid Chamran University of Ahvaz. Pp: 158. (in Persian with English abstract)

15- Pourseyyedi A. 2009. Groundwater model development of Jiroft plain. M.Sc. Thesis of Irrigation and Drainage, Shahid Chamran University of Ahvaz. Pp: 203. (in Persian with English abstract)

16- Reynolds D.A., and Marimuthu S. 2006. Deuterium composition and flow path analysis as additional calibration targets to calibrate groundwater flow simulation in a coastal wetlands system. Hydrogeology Journal, 15:515- 535.

17- Selle B., Minasny B., Bethune M., Thayalakumaran T., and Chandra S. 2011. Applicability of Richards' equation models to predict deep percolation under surface irrigation. Geoderma. 160:569-578.

18- Shirzadi S., and Sabouhi M. 2009. Application of multiobjective programming in surface and groundwater resources management in the Savojbolagh region. Agricultural Economics. 3(2): 83-98. (in Persian with English abstract)

19- Thorley M., and Callander P. 2005. Christhurch city groundwater model. Environment Canterbury report U05/53. Pp: 10.

20- Yang C.C., Chang L.C., Chen C.H., and Yeh M.S. 2008. Multi-objective planning for conjunctive use of surface and subsurface water using genetic algorithm and dynamics programming. Water Resource Management. DOI 10.1007/s11269-008-9281-5.