The Journal of Engineering Geology Vol. 30, No. 2, June, 2020, pp. 119-129 pISSN : 1226-5268 eISSN : 2287-7169 https://doi.org/10.9720/kseg.2020.2.119

RESEARCH ARTICLE

Assessment of the Effect of Sand Dam on Level: A Case Study in Chuncheon, South Korea

Bisrat Yifru1ㆍMin-Gyu Kim2ㆍSun Woo Chang3ㆍJeongwoo Lee4ㆍIl-Moon Chung5* 1Researcher, Department of Land, Water and Environment Research, KICT Student of Ph.D. Program, Department of Civil and Environmental Engineering, University of Science and Technology 2Research Specialist, Department of Land, Water and Environment Research, KICT 3Senior Researcher, Department of Land, Water and Environment Research, KICT Adjunct Professor, Department of Civil and Environmental Engineering, University of Science and Technology 4Research Fellow, Department of Land, Water and Environment Research, KICT 5Senior Research Fellow, Department of Land, Water and Environment Research, KICT Full-time Professor, Department of Civil and Environmental Engineering, University of Science and Technology

Abstract Sand dam is a successful water harvesting method in mountainous areas with ephemeral rivers. The success is dependent on several factors including material type, hydrogeology, slope, riverbed thickness, , and streamflow. In this study, the effect of a sand dam on the groundwater level in the Chuncheon area, South Korea was assessed using the MODFLOW model. Using the model, multiple scenarios were tested to understand the groundwater head before and after the construction of the sand dam. The effect of groundwater abstraction before and after sand dam construction and the sand material type were also assessed. The results show, the groundwater level increases substantially after the application of a sand dam. The comparison of model outputs, simulated groundwater head before and after sand dam application with and without pumping well, shows a clear difference in the head. The material type has also an effect on the groundwater head. As the conductivity of the material increases, the head showed a significant rise.

Keywords: sand dam, recharge, groundwater level, MODFLOW, Chuncheon

  

OPEN ACCESS In mountainous areas with ephemeral and intermittent rivers, water harvesting is *Corresponding author: Il-Moon Chung E-mail: [email protected] indispensable to secure ecological, agricultural, industrial, and domestic water supply demands in dry periods (Subagyono and Pawitan, 2008; Ertsen and Hut, 2009; Lasage et Received: 7 May, 2020 Revised: 21 May, 2020 al., 2015; De Trincheria et al., 2018). In this regard, sand storage dam is a successful Accepted: 21 May, 2020 water harvesting technology particularly in arid and semiarid environments (Quilis et Ⓒ 2020 The Korean Society of Engineering Geology al., 2009; Jadhav et al., 2012; De Trincheria et al., 2016). Sand dam harvests excess This is an Open Access article distributed under the terms of water, stores excess water mainly in the sand, to overcome drought periods in an the Creative Commons Attri- bution Non-Commercial License (http://creative- commons.org/licenses/by-nc/4.0/) which permits ephemeral river (Lasage et al., 2008, 2015; Quinn et al., 2019). Studies also show that unrestricted non-commercial use, distribution, and reproduction in any medium, provided the original sand dams help to overcome water scarcity during long-term climate changes (Aerts et work is properly cited. 120 ∙ Bisrat YifruㆍMin-Gyu KimㆍSun Woo ChangㆍJeongwoo LeeㆍIl-Moon Chung

al., 2007; Ryan and Elsner, 2016). Water harvesting using sand dams has several more advantages including minimizing evaporation, groundwater recharge, less pollution, and simplicity of construction (Baurne, 1984; Ertsen and Hut, 2009; Stern and Stern, 2011; Villani et al., 2018). In aquifers that have a full connection with the river, a sand dam may be considered as managed recharge (Quinn et al., 2019). Even though it is difficult to quantify, several researchers confirmed that the seepage from the riverbed and riverbank recharges the surrounding aquifer depending on the groundwater head (Hoogmoed, 2007; Hut et al., 2008; Quilis et al., 2009; Quinn et al., 2019). This phenomenon may last throughout the year. Consequently the groundwater level in the sand dam area increases (Yifru et al., 2018; Quinn et al., 2019). As the increases the natural vegetation-indigenous trees and riparian plants will get more water and live longer. The ecosystem restored and biodiversity has the potential to increase with the increase in groundwater level after the sand dam construction (Stern and Stern, 2011; Ryan and Elsner, 2016). Evaporation from sand is relatively low (Hellwig, 1973; Hut et al., 2008; Quinn et al., 2019). This makes the water available for a long time but it is dependent on the depth of the accumulated sand. Coarse gravel sand can be filled with water to 35 percent of its volume (Lasage et al., 2008; Maddrell and Neal, 2012). As the water level in the sand storage decrease to a few centimeters, the evaporation loss becomes significantly low (Baurne, 1984). And thus taking into consideration this concept in the study and construction of a sand dam is essential. The construction and the basic idea of water storage are relatively simple (Ertsen and Hut, 2009; Villani et al., 2018). The construction is mainly dependent on the locally available materials and it is economical (Lasage et al., 2008, 2015; De Trincheria et al., 2018). Moreover, the harvested water, using a sand dam, is clean because of the filtering effect of the sand material (Quinn et al., 2019). However, the successful application of a sand dam needs a good understanding of the hydrological processes in the dam area before and after the dam construction (Quilis et al., 2009). The socioeconomic aspect is also an important parameter, which needs to be investigated in the sand dam development process. Since the material accumulated on the upstream of the dam quality determines the ability of water-storing in it and transferring to the groundwater investigating all the possible issues is important. Silting is the most commonly reported problem which reduces the efficiency of a sand dam (De Trincheria et al., 2016). The efficiency and capacity can be improved by maximizing the accumulation of coarser materials and constructing a series large network of sand dams (De Trincheria et al., 2018). Selecting an area rich with several small upstream tributaries of the sand dam is also another way to maximize water storage and to get an adequate discharge at the dam site. Usually, more than one sand dam built in sequence to maximize the local storage capacity in a wider area or watershed scale. However, the total effect of the sand dam in the area may not be the aggregate sum of all the dams built. The result is site-specific and different for each sand dam. Therefore, it is important to evaluate each sand dam based on the specific site condition by including the most possible details in the numerical models (Quilis, 2007). If the sand dam site selection, height, downstream protection, and other hydrological changes are not studied well and acted accordingly the disadvantages may outweigh. Sand dam technology has a long-time history and implemented all over the world (Baurne, 1984; Hut et al., 2008; Yifru et al., 2018). The application is more common in arid and semi-arid areas. Recently it is widely used in Ethiopia               ∙ 121

and Kenya dry riverbeds to increase the local storage capacity for small schemes and drinking purposes. The community gets water from the sand dam area in different ways. The most common water abstraction ways are using embedded pipe near the bottom of the sand dam wall on the downstream side, using pumping wells installed on the riverbank, and in the rural area of Africa scoop holes are also other ways (Borst and de Haas, 2006; Quinn et al., 2018; Villani et al., 2018). However, the studies on hydrogeology and the effect of a sand dam on the groundwater level are limited (Aerts et al., 2007). This work focuses on the assessment of the effect of the sand dam on the groundwater level in Chuncheon, South Korea. The groundwater head was evaluated using the MODFLOW model before and after the construction of a sand dam. The effect of groundwater abstraction and conductivity of sand material on the fluctuation of the groundwater table was also assessed.

    

 The location of the sand dam has a substantial effect on its effectiveness. To bring the desired hydrological impact the location should be selected cautiously (Ryan and Elsner, 2016). The selection of a sand dam site is practically similar to any surface water-storing dams. Keeping the length of the dam wall short, maximized storage in a minimum surface area, accessibility, suitable slope that maximize storage and accumulation of sand, good material type in the riverbed, social aspect, and water availability are the key issues need to be addressed in the selection process (Nissen-Petersen, 2000; Nissen-Patterson, 2006; Pauw et al., 2008). To minimize the length of the dam but get the maximum possible storage without submerging large area of land or valuable property the site should be a narrow valley or deep gorge. The upstream area should be large enough and ideal slope to store the maximum amount of sand. The availability of local construction materials and thickness of the riverbed are also important things to consider (Baurne, 1984). The current study area has a pronounced elevation difference, which varies from 350 m to 240 m above sea level. Fig. 1 shows the topography and model boundary of the study area. The elevation difference downstream and upstream of the dam site is easily noticeable. The slope difference along the river (from point A to point B) is also depicted in Fig. 1. Since the DEM is 30 × 30 m of horizontal resolution, it is not easy to see the gorge. However, the overlain contour is based on field data and it shows the elevation difference and the gorge clearly. A river, which runs from north to south, crosses the conceptualized model area. The river has an average discharge of 1,755.27 m3/day. The flow in this river is seasonal; it has almost zero flow rate in dry periods. This data does not include the low flow in the winter. The flow measurement is mostly ignored at the winter low flows. Fig. 2 shows the average monthly precipitation pattern of the study area. The precipitation in the area varies from zero to 750 mm per month. The area gets the majority of the precipitation from June to September. The maximum precipitation occurs in August. In the modeling, the precipitation is used to approximate the groundwater recharge. 122 ∙ Bisrat YifruㆍMin-Gyu KimㆍSun Woo ChangㆍJeongwoo LeeㆍIl-Moon Chung

493350

A 252.05 254.85 256.00 254.60 251.97

253.77 252.96 250.15 252.89 257.46 254.80 248.55 249.66 251.19 248.82 Model boundary cp8251.74 251.00 247.32

249.86 255.32 253.15 600 247.38 249.49 249.25 249.32 248.86 245.14 255.25

244.35 246.16 254.74 248.29 244.47 244.30 493300 252.65 253.54 244.15

246.38 246.02 243.90 252.71 245.73 249.13 243.38 Profile of the river 252.21

244.42 246.66 254.11 244.48 242.87 242.94 246.01 250.45 400 244.20 242.87242.96 247.52 250.29 245.57 243.79 251.06

240.30 245.12244.57 from A to B

245.09 242.77 240.46 244.10 242.51 239.95 242.49242.35 242.17cp7

246.26 240.08 249.49 238.61 241.91 245.25 247.27 243.16 240.45 238.43 251.04 241.23 239.76 237.34 239.11 237.78 200 245.96 240.64

241.96 244.07 237.09 240.08 240.15 236.53

238.76 236.77 238.76 239.80 242.55 Dam site 238.73 238.07238.22 236.41 238.06 243.71243.85 247.88 cp6 234.89 234.78 240.38 237.52 250.54 236.74 235.62 236.32 244.15 237.17 245.30 233.60 233.49 235.08 233.24 234.54 237.02 0 243.43 231.78 237.59 231.48231.85 233.18 233.56 236.54 240.50 238.25 232.99 231.37231.78 244.71 493250 238.37 231.20 232.05 240.98 0 200 400 600 800 1000 1200 1400 231.75 237.58 231.68 231.40231.45231.75231.77231.75 231.38 231.47 229.98 232.12 230.80 237.22 229.79 231.68 230.69 230.53 241.32 231.97 229.79 230.30 229.15 232.74 232.19 229.63 230.37 228.78 229.56 231.48 233.24 232.53 230.35230.15 230.15 231.72 230.97 230.38 240.65 228.79 228.96 231.22 228.70 235.76 229.64229.50 229.54 228.21 230.89 230.47 232.10 229.67 236.88 228.62 228.32 227.60 229.95 350 m 230.30 230.52 229.12 235.99

229.56

229.17 237.81 230.13 229.70 228.43 226.22 228.38 229.33 226.65 234.92

229.45 227.99 226.22 236.74

227.25 227.12 226.56 232.57 228.59 226.31 232.05 232.19 226.91 233.12 227.85 231.41 227.27

227.31227.32 227.36 cp3 233.53 230.89 227.77 227.97 225.45 cp4 227.36 234.19

232.89 231.22 227.62 226.49 228.21 234.88 E 232.47 232.08 233.97

231.62 232.94 234.89 235.95 230.43 493200 232.15 227.55 227.42 224.74 233.87 225.66 232.59 230.54 234.85 235.92 228.92

230.83 227.46 227.38 227.48 229.23 225.29

229.25

227.87 232.52 227.41 231.28

223.78 226.82

224.15 227.73 227.10 227.18227.17227.06 133

226.90

226.52E

223.90 226.18 223.45

223.04

222.95 223.35

224.97

222.52 227.64 222.27 224.63 224.82 240 m

223.77

223.73 223.95223.82 113

226.70

223.70 223.47

220.79

220.73

493150 223.01 223.09 B 221.61 221.45

220.81 220.56 E 220.57

220.01 220.08

219.49219.43 111 493100 218.74

218.74

256550 256600 256650 256700 256750 256800        

           

      The effect of a sand dam on the groundwater level was tested by modeling a small area near the proposed sand dam site using the MODFLOW model. The MODFLOW solves the three-dimensional flow of groundwater of constant density through porous earth material which can be described by Eq. (1) (Langevin et al., 2003).

            ±    (1)                 

where, Kxx, Kyy, Kzz are principal components of the hydraulic conductivity tensor in x, y, z special directions, W is the source or sink, Ss is specific storage (1/L), h is the hydraulic head (L), and t is time.               ∙ 123

The model is built as a single layer with a depth of 100 m and discretized into four mathematical layers to see the fluctuation of the groundwater level with different scenarios developed. The model has a 100 × 100-cell size which resulted in 21320 active cells. The river flow was modeled using a river package in MODFLOW. The river package was used to simulate the influence of the river water level on groundwater. The package needs river stage, riverbed bottom, and conductance, which are all difficult to measure. In this model, riverbed bottom was approximated from DEM, river stage from measured data, and the conductance is set by trial and error. The entrance and the exit of the river are assigned as specified head boundary (by parallel flow concept) and the rest of the model is considered as a no-flow boundary. Surface recharge is another important boundary condition. In this study, the recharge was approximated from the precipitation in the area and applied to the upper-most active cells. Nearly 20 percent of the precipitation was considered as surface recharge. In this study, six numerical models were built based on the current site information and some assumptions. The first two models are bases used to understand the groundwater condition before the construction of the sand dam and with no pumping wells. The second two models were built to evaluate the effect of the sand dam on the groundwater level with and without pumping. And the last two models were built to assess the effect of the material with and without pumping wells. All six models are based on the same grid size and surface recharge. To understand the groundwater head in the study area before the construction of sand dam two models were built. The first model represents the site condition with no sand dam and no pumping wells. It was built to have a good background on the natural hydrological condition and use it as a base of comparison (as a reference). On the second model, three pumping wells with a pumping capacity of 10, 20, and 30 m3/day were included using the well package of MODFLOW. Two of the wells with the higher pumping capacity are placed on the upstream area of the proposed sand dam site, the other one at the downstream.

     

The model results are presented in Figs. 3 to 9. The presented results are the head contours of each scenario developed including the base scenario. The discussion is focused on the comparison of the head contours with the base scenario.

    The first scenario-the simulated head contours before the sand dam construction with no pumping is presented in Fig. 3. As can be seen from the figure, the head varies from 295 m to 233 m. The head variation follows the topography and has moderate variation from upstream to downstream in general. The head around the sand dam site is 253.8 m. However, the head change is sudden from the dam site to downstream of it. The reason is mainly the elevation change. As the flow vectors show, the flow is from north to south parallel to the river flow and the river and aquifer interaction is at limited places with a river leakage of 602.4 m3/day. At the upstream and downstream section of the model, the flow is to the river, around the dam site the flow is from the river to the aquifer. 124 ∙ Bisrat YifruㆍMin-Gyu KimㆍSun Woo ChangㆍJeongwoo LeeㆍIl-Moon Chung

               

Fig. 4 shows the simulated groundwater potential of the second model-after the pumping wells are added. The head shows a substantial change; at the sand dam site, the head dropped to 245 m, which is eight meters less than the head before the pumping is applied. Especially downstream of the dam site the head dropped due to the pumping. The river linkage, in this case, is almost zero. From this, it can be inferred that groundwater potential is limited to abstract around 60 m3/day without affecting the groundwater level.

                             ∙ 125

   A sand dam height should be limited so that the water level and the maximum flood level will be not higher than the riverbank. The flood level is a serious case in the sand dam fixing process. If the flood level rises above the height of the riverbank after the construction of the sand dam, the area will be endangered from flood risk, erosion, and failure of the sand dam. The sand dam constructed in this study area has a height of three-meter at the wings and two-meter in the middle (spillway). The third, fourth, and fifth models were after the sand dam construction-considering the increase in the water level of the river in the upstream area. The riverbed around the dam was also elevated due to the accumulation of sand and considered in the model and therefore, the head increase considered as constant head cells. Fig. 5 shows the simulated groundwater head after the construction of a sand dam. Compared to the first case (Fig. 3) here the head shows significant change especially closer to the sand dam both in upstream and downstream of the dam. At the dam site, the head increased to 267.5 m. In the dam area and downstream of it, the model shows flooded cells. The river leakage is also reduced to 34.24 m3/day.

              

On the same model setup, the pumping wells are added into the model to see the effect of the abstraction on the groundwater level (Fig. 6). Compared to case two (Fig. 4) the effect of pumping wells on the groundwater head is not significant. The comparison of Figs. 5 and 6 shows also the groundwater head is almost the same. The head at the dam stayed the same as case three (sand dam no pumping wells). This implies that the effect of the sand dam on the increase of groundwater potential is significant and we can abstract water without affecting the aquifer. Moreover, in the dam area, the river shows a much higher water level than before the construction of the sand dam. 126 ∙ Bisrat YifruㆍMin-Gyu KimㆍSun Woo ChangㆍJeongwoo LeeㆍIl-Moon Chung

         

The first four cases are compared graphically with the simulated head along the river (Fig. 7). The head difference before the sand dam and after the sand dam is substantial, especially around the dam site.

                         

Several trials were performed to evaluate the effect of the material types on the effectiveness of the sand dam on the groundwater head. The sand type (material type in general) as Figs. 8 and 9 show the effect of the increase of hydraulic               ∙ 127

conductivity of the sand accumulated upstream of the sand dam. As the conductivity increase, the groundwater head shows a significant change, especially in the downstream area. In this model, the conductivity increased from 0.256 to 20 m/day to see the effect clearly and the result is as shown in the figure. The downstream area is flooded due to the increase in the head.

                        

                     128 ∙ Bisrat YifruㆍMin-Gyu KimㆍSun Woo ChangㆍJeongwoo LeeㆍIl-Moon Chung

    

A sand dam has several advantages in the water harvesting process. In addition to storing water for later use, it gets the water filtered, protect from getting polluted, and minimize evaporation. Moreover, the environmental benefit of a sand dam is several folds. In this work, the influence of the sand dam on groundwater was tested using a numerical model-MODFLOW. Before the application of sand dam, the groundwater head in the area with pumping wells and without pumping wells was evaluated using models. Then the effect of the sand dam was evaluated both with and without pumping wells by comparing the head with the simulated head result before the sand dam application. The effect of material type accumulated upstream of the dam was also evaluated by varying the conductivity. As the results show, the groundwater level increases substantially after the application of a sand dam. The comparison of model outputs, simulated groundwater head before and after sand dam application with and without pumping well, shows a clear head difference. The material type has also an effect on the groundwater head. As the conductivity of the material increases, the head, especially at the downstream, showed a significant increase. After the construction of the sand dam, the groundwater level is not affected much. In this study, groundwater level after the construction of a sand dam with pumping wells and without pumping wells is the same. In this work, the hydraulic conductivity, river conductance, and other important flow parameters are obtained by trial and error. The material type accumulated upstream of the dam helps to increase the reliability of the model results. This preliminary work can be improved with further hydrogeological investigations and calibration of the simulation with measured groundwater head in the area.

   

This study was funded by the Korea Ministry of Environment (MOE) as a Demand Responsive Water Supply Service Program (146515).

    

Aerts, J., Lasage, R., Beets, W., de Moel, H., Mutiso, G., Mutiso, S., de Vries, A., 2007, Robustness of sand storage dams under climate change, Vadose Zone Journal, 6(3), 572-580. Baurne, G., 1984, “Trap-dams”: Artificial subsurface storage of water, Water International, 9(1), 2-9. Borst, L., de Haas, S.A., 2006, of sand storage dams: A case study in the Kiindu catchment, Kitui district, Kenya, Masters Thesis, VU University Amsterdam, 25-68. De Trincheria, J., Leal, W.F., Otterpohl, R., 2018, Towards a universal optimization of the performance of sand storage dams in arid and semi-arid areas by systematically minimizing vulnerability to siltation: A case study in Makueni, Kenya, International Journal of Sediment Research, 33(3), 221-233. De Trincheria, J., Wibbing, J., Leal Filho, W., Otterpohl, R., 2016, Practical recommendations to prevent, restore and rehabilitate silted-up sand storage dams in arid and semi-arid areas, Proceedings of the 7th Rural Water Supply Network Forum “Water for Everyone”, Abidjan, 1-14.               ∙ 129

Ertsen, M., Hut, R., 2009, Two waterfalls do not hear each other. Sand-storage dams, science and sustainable development in Kenya, Physics and Chemistry of the Earth, Parts A/B/C, 34(1-2), 14-22. Hellwig, D.H.R., 1973, Evaporation of water from sand, 2: Diurnal variations, Journal of Hydrology, 18(2), 109-118. Hoogmoed, M., 2007, Analyses of impacts of a sand storage dam on groundwater flow and storage: Groundwater flow modeling in Kitui district, Kenya, Master Thesis, VU University Amsterdam, 20-28. Hut, R., Ertsen, M., Joeman, N., Vergeer, N., Winsemius, H., van de Giesen, N., 2008, Effects of sand storage dams on groundwater levels with examples from Kenya, Physics and Chemistry of the Earth, Parts A/B/C, 33(1-2), 56-66. Jadhav, M.V., Shaikh, A.P., Gite, B.E., Yadav, A.P., 2012, Sand dam reservoir - Need of semi arid areas, International Journal of Engineering Research and Applications, 2(6), 1690-1694. Langevin, C.D., Shoemaker, W.B., Guo, W., 2003, MODFLOW-2000, the U.S. geological survey modular ground-water model--documentation of the SEAWAT-2000 version with the variable-density flow process (VDF) and the integrated MT3DMS transport process (IMT), U.S. Geological Survey Open-File Report 2003-426, 43p. Lasage, R., Aerts, J., Mutiso, G.C.M., de Vries, A., 2008, Potential for community based adaptation to droughts: Sand dams in Kitui, Kenya, Physics and Chemistry of the Earth, Parts A/B/C, 33(1-2), 67-73. Lasage, R., Aerts, J.C.J.H., Verburg, P.H., Sileshi, A.S., 2015, The role of small scale sand dams in securing water supply under climate change in Ethiopia, Mitigation and Adaptation Strategies for Global Change, 20(2), 317-339. Maddrell, S., Neal, I., 2012, Sand dams: A practical guide, Excellent Development, 62-81. Nissen-Petersen, E., 2000, Water from sand rivers: A manual on site survey, design, construction and maintenance of seven types of water structures in riverbeds, RELMA, Technical Handbook Series 23, 67p. Nissen-Patterson, E., 2006, Water from dry riverbeds, ASAL Consultants Ltd., 1-19. Pauw, W.P., Mutiso, S., Mutiso, G., Manzi, H.K., Lasage, R., Aerts, J.C.J.H., 2008, An assessment of the social and economic effects of the Kitui sand dams, SASOL & Institute for Environmental Studies, 49-50. Quilis, R.O., 2007, Modelling sand storage dams systems in seasonal rivers in arid regions: Application to Kitui district (Kenya), Master Thesis, Delft University of Technology, 1-57. Quilis, R.O., Hoogmoed, M., Ertsen, M., Foppen, J.W., Hut, R., de Vries, A., 2009, Measuring and modeling hydrological processes of sand-storage dams on different spatial scales, Physics and Chemistry of the Earth, Parts A/B/C, 34(4-5), 289-298. Quinn, R., Avis, O., Decker, M., Parker, A., Cairncross, S., 2018, An assessment of the microbiological water quality of sand dams in southeastern Kenya, Water, 10(6), 708. Quinn, R., Rushton, K., Parker, A., 2019, An examination of the hydrological system of a sand dam during the dry season leading to water balances, Journal of Hydrology X, 4, 100035. Ryan, C., Elsner, P., 2016, The potential for sand dams to increase the adaptive capacity of East African drylands to climate change, Regional Environmental Change, 16, 2087-2096. Stern, J.H., Stern, A., 2011, Water harvesting through sand dams, ECHO Technical Note No. 70, http://edn.link/q9xrkd. Subagyono, K., Pawitan, H., 2008, Water harvesting techniques for sustainable water resources management in catchments area, Proceedings of the International Workshop on Integrated Watershed Management for Sustainable Water Use in a Humid Tropical Region, Tsukuba, 18-30. Villani, L., Castelli, G., Hagos, E.Y., Bresci, E., 2018, Water productivity analysis of sand dams irrigation farming in northern Ethiopia, Journal of Agriculture and Environment for International Development, 112(1), 139-160. Yifru, B., Kim, M.G., Chang, S.W., Lee, J., Chung, I.M., 2018, Numerical modeling of the effect of sand dam on groundwater flow, The Journal of Engineering Geology, 28(4), 529-540.