Estimation of Natural and Artificial Recharge of Shahreza Plain Groundwater in Isfahan Using CRD and Hantush Models
Total Page:16
File Type:pdf, Size:1020Kb
American Journal of Engineering and Applied Sciences Original Research Paper Estimation of Natural and Artificial Recharge of Shahreza Plain Groundwater in Isfahan Using CRD and Hantush Models 1Mohsen Nekooei, 1Jahangir Abedi Koupai, 1Saeid Eslamian, 2Vijay P. Singh and 3Kaveh Ostad-Ali-Askari 1Department of Water Engineering, College of Agriculture, Isfahan University of Technology, Isfahan, Iran 2Department of Biological and Agricultural Engineering and Zachry, Department of Civil Engineering, Texas A and M University, 321 Scoates Hall, 2117 TAMU, College Station, Texas, USA 3Department of Civil Engineering, Najafabad Branch, Islamic Azad University, Najafabad, Iran Article history Abstract: Cumulative Rainfall Departure model (CRD) is widely used for Received: 13-11-2018 estimation of ground natural recharge. This model is based on the Revised: 13-01-2020 groundwater balance and it requires random parameters. Moreover, an Accepted: 08-06-2020 increase in the static level under artificial recharge facility is simulated using diffusion equations. Since Shahreza Plain is located in an arid and Corresponding Author: Mohsen Nekooei semi-arid area and given that rainfall is the major recharge source in such Department of Water regions, in this study CRD model is applied in order to estimate natural Engineering, College of recharge of Shahreza Plain groundwater. The amount of groundwater basin Agriculture, Isfahan University fed by rainfall was estimated about 177.2 million cubic meters using the of Technology, Isfahan, Iran CRD model. Also, by investigating and comparing available artificial Email: [email protected] recharge plans in Shahreza Plain it was known that in Kohrouyeh Project due to high soil permeability the increase in water table depth was higher than other projects and it achieves the maximum level of the water table within less time. By investigating the course of groundwater quality change in upstream and downstream boreholes of artificial recharge plans and given dewatering periods, we found that artificial recharge plans are generally effective in improving groundwater quality in the direction of groundwater movement but these changes are not significant. Hence, it is demonstrated that groundwater change mostly depends on the groundwater drop. Keywords: Artificial Recharge, Shahreza Plain, Groundwater Quality, Cumulative Rainfall Departure Model, Hantush Model Introduction surface currents join groundwater resources known as an effective recharge and it is defined as the amount One of the effective ways to supply the water of water penetrated into lower layers of the unsaturated required by different sectors, particularly in arid and area and heightens water table (Barry et al ., 1999). semi-arid regions and also overcoming the crisis of In South Africa, simulation of water fluctuation in drought is to assist in increasing groundwater water table against rainfall time was conducted in reservoirs and preventing Qanat dry out, control of order to determine recharge value. The result of this floods and basins artificial recharge (Behnia, 1988). method is the same as that in CRD. This model is In the last few decades, artificial recharge methods for applicable to those regions where the only recharge restoring groundwater resources dried out due to source is rainfall. Therefore, concerning the regions excessive exploitation in arid regions has been applied where there is water flow between basins, it may not in regions where there are coarse grains and moderate be applicable (Gerdefaramarzi and Koupaei, 2011). slope alluvial fans (Minnig et al ., 2018). Khazaei et al . (2003) studied groundwater recharge Rainfall is the main source of groundwater in arid and semi-arid regions of the south of Iran and resources. Every year, a certain amount of rainfall or they proposed a model in which estimated water table © 2020 Mohsen Nekooei, Jahangir Abedi Koupai, Saeid Eslamian, Vijay P. Singh and Kaveh Ostad-Ali-Askari. This open access article is distributed under a Creative Commons Attribution (CC-BY) 3.0 license. Mohsen Nekooei et al., / American Journal of Engineering and Applied Sciences 2020, 13 (2): 283.295 Aida DOI: 10.3844/ajeassp.2020.283.295 3 is consistent to the water table measured from Qp = Pumping Outflow L /T 3 observed boreholes. Sensitivity analysis indicated that Qout = Aquifer Outflow L /T 3 when rainfall in these regions is replaced by the Qin = Aquifer Inflow L /T rainfall in moderate and highlands, recharge value is ∆h i = Water level change within certain time reduced by 42-87%. The results indicated that in order interval [L] to achieve an exact estimation of recharge it is A = Aquifer area [L2] necessary to divide the region into similar hydrologic S = aquifer storability coefficient areas (Rebhun and Schwarz, 1968). Artificial recharge includes water into a permeable Bredenkamp assumed a linear relationship between formation; this method is widely used for protecting monthly changes of water level and cumulative aquifers and utilization of stored water under the rainfall departure. This relationship may be expressed ground (Bouwer, 2002). as follows: Matson 1961 believes that high permeability and = + − soil depth are two factors that make the soil suitable CRDi CRD i−1 ( p i C ) (2) for artificial recharge. He believes that light soil lands are the best-qualified ones for artificial recharge. Where: th So many studies have been done about artificial CRD i = Cumulative Rainfall Departure for i month recharge principles, water penetration into soil and (mm.); [L] th formation of hydraulic mound under recharge CRD i-1 = Cumulative Rainfall Departure for–i1 facilities that have led to the proposition of theoretical (mm.); [L] 3 equations for expressing hydraulic mound behavior Pi = Rainfall (mm./month) L /T (Cook et al ., 1989). By reaching wetting front to the C = A threshold value where recharge does water table, the hydraulic mound develops, notoccur anymore development of hydraulic mound under the modified bed of rivers may be considered as 2D (Khazaei et al ., 2003). This equation indicates the relationship between In order to evaluate the impacts of artificial the two values of CRD in two consecutive periods. The recharge on the water table, many computer programs general equation used for converting pure recharge have been developed. One of them is Glover and from the CRD model is as follows: Hantush’s analytical solution for rectangular and circular ponds that predicts the value of the mound R Q− Q Q =++in out − p growth and shape (Behnke, 1969). hi h i −1 (3) S SA× SA × Groundwater recharge via rainfall is the most important recharge source in Southern Shahreza Area. where, represents pure recharge expressed as a The average annual rainfall in this region is 204 mm coefficient of rainfall. h is the frequency in ith month of which a small part is stored and much of it’s i and h -1 is the frequency in i-1th month calculated on the evaporates. The aim of the present study is the i meter. The values of pumping outflow and rainfall estimation of natural recharge value in Southern statistics are usually measured accurately. Therefore, if Shahreza Plain using CRD model and investigating the the inflow and storability coefficient is known, the only impact of artificial recharge plans on the quantity and unknown parameter in above equation will be R. Using quality of groundwater of this region (WMSD, 2007). CRD model water table fluctuations obtained from Materials and Methods Equation (3) is portrayed against measured water table measured from observed boreholes. The best recharge Cumulative Rainfall Departure Model estimation is achieved when the difference between measured and estimated water table values are CRD model was initially proposed by minimized that will be possible through shrinking the Bredenkamp et al ., 1995. The theory of this method is amount of root mean square error between measured and based on groundwater balance in an aquifer (Raeisi and calculated water table values. Afzal, 1997). This balance may be expressed as follows: = + − +∆×× Hantush’s Model RQQp out Q in hAS i (1) The equation that predicts the value of groundwater Where: level rise or drops in artificial recharging ponds proposed 3 R = Total Aquifer Recharge L /T by (Hantush, 1967). This model could yield a good 284 Mohsen Nekooei et al., / American Journal of Engineering and Applied Sciences 2020, 13 (2): 283.295 Aida DOI: 10.3844/ajeassp.2020.283.295 estimation of groundwater bed recharge value and it is where, ts represents the time of water infiltration stoppage. applied for rectangular recharge level (Fig. 1). Since at the time of groundwater level decline, the Prediction of hydraulic mound development of an area of water saturation is drained, f value must be open aquifer under rectangular recharge ponds is considered for drainage conditions that are usually less performed via the following equation: than real f value (Xu and van Tonder, 2001). Vt W L Shahreza Plain Study −=a + + hx, y , t H F xn , yn 4f 2 2 Shahreza County is located at 85th Km of the south 2 W L of Isfahan. Shahreza’s 2500 km basin originates from +F + xn , − yn 2 2 Gavkhouni basin chief branches and Shahreza County is WL WL located in its outlet. +−F xn , ++ ynF − xn , − yn 22 22 Climatic conditions of this region according to De 1 (4) Martonne and Coupon’s method is arid and it is arid- = =−0.5 × − 0.5 F()αβ, errorfunction. ∫ erf()() α t erf β t cold based on Emberger’s method. The average annual 0 rainfall of this region is 204 mm and its annual 4tT n = ^()− 0.5 evaporation rate is 2650 mm. Where the maximum value f of rainfall in July equals 410 mm and the minimum W value of rainfall in January equals 65 mm. The α = ± x n 2 maximum absolute temperature of this region is 42°C L and minimum temperature is-27°C. β = ± y n 2 The formations of the second geologic period have produced one of the most important water reservoirs of this basin.