Estimating high hydraulic conductivity locations through a 3D simulation of water flow in soil and a resistivity survey
Keisuke Inoue 1 3 Hiroomi Nakazato 1 Tomijiro Kubota 1 Koji Furue 2 Hiroshi Yoshisako 1 Michiaki Konno 1 Daisuke Shoda 11 Institute for Rural Engineering, National Agriculture and Food Research Organisation, 2-1-6, Kannondai, Tsukuba-shi, Ibaraki 305-8609, Japan.
2 Kagoshima Prefectural Institute for Agricultural Development, 2200, Ono, Kinpo-cho, Minamisatsuma-shi, Kagoshima 899-3401, Japan.
3 Corresponding author. Email: ksk@affrc.go.jp
Exploration Geophysics 49(3) 299-308 https://doi.org/10.1071/EG17054
Submitted: 10 April 2017 Accepted: 17 April 2017 Published: 30 June 2017
Abstract
In this study, we propose a method to estimate high hydraulic conductivity locations that uses 3D simulation of soil water flow and in-line resistivity survey data acquired during a groundwater recharge experiment, and we apply this method to numerical and field experiments. The high hydraulic conductivity locations are estimated from a combination of field-observed and simulated apparent resistivities using the following simple steps. (1) Assuming that high hydraulic conductivity zones exist in the first layer, simulations of saturated-unsaturated seepage are conducted for several possible water-flow models that have high hydraulic conductivity zones in different locations. (2) The simulated volumetric water contents are converted into bulk resistivities, which are used to produce apparent resistivity data through simulation of a resistivity survey. (3) The differences between the simulated apparent resistivities and the field-observed data are examined, and the best-fit hydraulic conductivity model is identified by minimising the above differences. In the numerical experiment, 3D inversion of the simulated resistivity survey provides an image of the preferential flow, although the infiltration locations are unclear. Comparing the field model with the possible models, the high hydraulic conductivity location in the field model corresponds to the high hydraulic conductivity location in the possible model with the minimum errors. In the field, an in-line resistivity survey was conducted during a groundwater recharge experiment on a pyroclastic plateau. The 3D inversion of the in-line resistivity survey data provides an image of the preferential flow. Comparing the field apparent resistivity data with the simulated apparent resistivity data, the high hydraulic conductivity location of the possible model that provides the minimum error corresponds to the recharge water range, whereas the hydraulic conductivity location of the possible model that gives the maximum errors corresponds to ranges with no recharge water. These results indicate that it is possible to estimate high hydraulic conductivity locations using 3D simulations of the soil water flow and a resistivity survey.
Key words: infiltration path estimation, in-line data, resistivity survey, water flow in soil, water flow simulation.
References
Amidu, S. A., and Dunbar, J. A., 2007, Geoelectric studies of seasonal wetting and drying of a Texas Vertisol: Vadose Zone Journal, 6, 511–523| Geoelectric studies of seasonal wetting and drying of a Texas Vertisol:Crossref | GoogleScholarGoogle Scholar |
Archie G. E. 1942
Batlle-Aguilar, J., Schneider, S., Pessel, M., Tucholka, P., Coquet, Y., and Vachier, P., 2009, Axisymetrical infiltration in soil imaged by noninvasive electrical resistivimetry: Soil Science Society of America Journal, 73, 510–520
| Axisymetrical infiltration in soil imaged by noninvasive electrical resistivimetry:Crossref | GoogleScholarGoogle Scholar | 1:CAS:528:DC%2BD1MXjt12rtLo%3D&md5=bd1f50f0539d4afc2566eabbac24d9a3CAS |
Binley, A., Shaw, B., and Henry-Poulter, S., 1996, Flow pathways in porous media: electrical resistance tomography and dye staining image verification: Measurement Science & Technology, 7, 384–390
| Flow pathways in porous media: electrical resistance tomography and dye staining image verification:Crossref | GoogleScholarGoogle Scholar | 1:CAS:528:DyaK28Xhs1GhsLs%3D&md5=195099d8df88756215471d3a3f233d64CAS |
Binley, A., Winship, P., West, L. J., Pokar, M., and Middleton, R., 2002, Seasonal variation of moisture content in unsaturated sandstone inferred from borehole radar and resistivity profiles: Journal of Hydrology, 267, 160–172
| Seasonal variation of moisture content in unsaturated sandstone inferred from borehole radar and resistivity profiles:Crossref | GoogleScholarGoogle Scholar |
Coscia, I., Greenhalgh, S. A., Linde, N., Doetsch, J., Marescot, L., Günther, T., Vogt, T., and Green, A. G., 2011, 3D crosshole ERT for aquifer characterization and monitoring of infiltrating river water: Geophysics, 76, G49–G59
| 3D crosshole ERT for aquifer characterization and monitoring of infiltrating river water:Crossref | GoogleScholarGoogle Scholar |
Daily, W., Ramirez, A., LaBrecque, D., and Nitao, J., 1992, Electrical resistivity tomography of vadose water movement: Water Resources Research, 28, 1429–1442
| Electrical resistivity tomography of vadose water movement:Crossref | GoogleScholarGoogle Scholar |
Daily, W., Ramirez, A., Binley, A., and LeBrecque, D., 2004, Electrical resistance tomography: The Leading Edge, 23, 438–442
| Electrical resistance tomography:Crossref | GoogleScholarGoogle Scholar |
Farzamian, M., Santos, F. A. M., and Khalil, M. A., 2015a, Application of EM38 and ERT methods in estimation of saturated hydraulic conductivity in unsaturated soil: Journal of Applied Geophysics, 112, 175–189
| Application of EM38 and ERT methods in estimation of saturated hydraulic conductivity in unsaturated soil:Crossref | GoogleScholarGoogle Scholar |
Farzamian, M., Santos, F. A. M., and Khalil, M. A., 2015b, Estimation of unsaturated hydraulic parameters in sandstone using electrical resistivity tomography under a water injection test: Journal of Applied Geophysics, 121, 71–83
| Estimation of unsaturated hydraulic parameters in sandstone using electrical resistivity tomography under a water injection test:Crossref | GoogleScholarGoogle Scholar |
Gharibi, M., and Bentley, L. R., 2005, Resolution of 3-D electrical resistivity images from inversions of 2-D orthogonal lines: Journal of Environmental & Engineering Geophysics, 10, 339–349
| Resolution of 3-D electrical resistivity images from inversions of 2-D orthogonal lines:Crossref | GoogleScholarGoogle Scholar |
Huisman, J. A., Rings, J., Vrugt, J. A., Sorg, J., and Vereecken, H., 2010, Hydraulic properties of a model dike from coupled Bayesian and multi-criteria hydrogeophysical inversion: Journal of Hydrology, 380, 62–73
| Hydraulic properties of a model dike from coupled Bayesian and multi-criteria hydrogeophysical inversion:Crossref | GoogleScholarGoogle Scholar |
Inoue, K., Nakazato, H., Kubota, T., Takeuchi, M., Sugimoto, Y., Kim, H. J., and Furue, K., 2016, Three-dimensional inversion of in-line resistivity data for monitoring a groundwater recharge experiment in a pyroclastic plateau: Exploration Geophysics, ,
| Three-dimensional inversion of in-line resistivity data for monitoring a groundwater recharge experiment in a pyroclastic plateau:Crossref | GoogleScholarGoogle Scholar |
Inoue, K., Nakazato, H., Takeuchi, M., Sugimoto, Y., Kim, H. J., Yoshisako, H., Konno, M., and Shoda, D., 2017, Investigation of the line arrangement of 2D resistivity surveys for 3D inversion: Exploration Geophysics, ,
| Investigation of the line arrangement of 2D resistivity surveys for 3D inversion:Crossref | GoogleScholarGoogle Scholar |
Jackson, P. D., Earl, S. J., and Reece, G. J., 2001, 3D resistivity inversion using 2D measurements of the electric field: Geophysical Prospecting, 49, 26–39
| 3D resistivity inversion using 2D measurements of the electric field:Crossref | GoogleScholarGoogle Scholar |
Koestel, J., Kemna, A., Javaux, M., Binley, A., and Vereecken, H., 2008, Quantitative imaging of solute transport in an unsaturated and undisturbed soil monolith with 3-D ERT and TDR: Water Resources Research, 44, W12411
| Quantitative imaging of solute transport in an unsaturated and undisturbed soil monolith with 3-D ERT and TDR:Crossref | GoogleScholarGoogle Scholar |
Kubota, T., Masumoto, T., Matsuda, S., and Furue, K., 2005, Hydrogeological features of Kasanohara plateau based on water quality and its circulation: Technical Report of the National Institute of Rural Engineering, Issue 203, 81–100 [in Japanese with English abstract].
LaBrecque, D. J., and Yang, X., 2001, Difference inversion of ERT data: a fast inversion method for 3-D in situ monitoring: Journal of Environmental & Engineering Geophysics, 6, 83–89
| Difference inversion of ERT data: a fast inversion method for 3-D in situ monitoring:Crossref | GoogleScholarGoogle Scholar |
Li, J., Li, X., Lv, N., Yang, Y., Xi, B., Li, M., Bai, S., and Liu, D., 2015, Quantitative assessment of groundwater pollution intensity on typical contaminated sites in China using grey relational analysis and numerical simulation: Environmental Earth Sciences, 74, 3955–3968
| Quantitative assessment of groundwater pollution intensity on typical contaminated sites in China using grey relational analysis and numerical simulation:Crossref | GoogleScholarGoogle Scholar |
Looms, M. C., Binley, A., Jensen, K. H., Nielsen, L., and Hansen, T. M., 2008, Identifying unsaturated hydraulic parameters using an integrated data fusion approach on cross-borehole geophysical data: Vadose Zone Journal, 7, 238–248
| Identifying unsaturated hydraulic parameters using an integrated data fusion approach on cross-borehole geophysical data:Crossref | GoogleScholarGoogle Scholar |
Nakano, M., Amemiya, Y., Muto, I., Shiozawa, S., and Nakamura, T., 1981, Physical and hydrologic properties of a Shirasu Hill — soil properties and water movement in a volcanic deposit (Shirasu) area (II): Transactions of The Japanese Society of Irrigation, Drainage and Reclamation Engineering, , Issue 93, 7–12
Naruo, H., and Kobayashi, T., 2002, Two large-scale earthquakes triggered by a 6.5ka BP eruption from Kikai caldera, southern Kyushu, Japan: The Quaternary Research (Daiyonki-Kenkyu), 41, 287–299
| Two large-scale earthquakes triggered by a 6.5ka BP eruption from Kikai caldera, southern Kyushu, Japan:Crossref | GoogleScholarGoogle Scholar |
Nimmo, J. R., Perkins, K. S., Schmidt, K. M., Miller, D. M., Stock, J. D., and Singha, K., 2009, Hydrologic characterization of desert soils with varying degrees of pedogenesis: 1. Field experiments evaluating plant-relevant soil water behaviour: Vadose Zone Journal, 8, 480–495
| Hydrologic characterization of desert soils with varying degrees of pedogenesis: 1. Field experiments evaluating plant-relevant soil water behaviour:Crossref | GoogleScholarGoogle Scholar |
Nishigaki, M. 2001, AC-UNSAF3D user’s manual [in Japanese].
Oldenborger, G. A., Knoll, M. D., Routh, P. S., and LaBrecque, D. J., 2007, Time-lapse ERT monitoring of an injection/withdrawal experiment in a shallow unconfined aquifer: Geophysics, 72, F177–F187
| Time-lapse ERT monitoring of an injection/withdrawal experiment in a shallow unconfined aquifer:Crossref | GoogleScholarGoogle Scholar |
Park, S., 1998, Fluid migration in the vadose zone from 3-D inversion of resistivity monitoring data: Geophysics, 63, 41–51
| Fluid migration in the vadose zone from 3-D inversion of resistivity monitoring data:Crossref | GoogleScholarGoogle Scholar |
Sasaki, Y., 1994, 3-D resistivity inversion using the finite-element method: Geophysics, 59, 1839–1848
| 3-D resistivity inversion using the finite-element method:Crossref | GoogleScholarGoogle Scholar |
Sugimoto, Y., 1995, Monitoring of electrolyte tracer in groundwater using resistivity tomography—numerical experiment: Proceedings of the 92nd SEGJ Conference, 57–62 [in Japanese].
Sugimoto, Y., Nakazato, H., Takeuchi, M., Kim, H. J., Inoue, K., Yamada, N., and Aono, T., 2004, Practical 3-D electrical resistivity survey method using measurements from a few 2-D survey lines: Proceedings of the 111th SEGJ Conference, 165–168 [in Japanese].
Suzuki, K., Oyama, T., Kawashima, F., Tsukada, T., and Jyomori, A., 2010, Monitoring of grout material injected under a reservoir using electrical and electromagnetic surveys: Exploration Geophysics, 41, 69–79
| Monitoring of grout material injected under a reservoir using electrical and electromagnetic surveys:Crossref | GoogleScholarGoogle Scholar |
Takakura, S., Yoshioka, M., Ishizawa, T., and Sakai, N., 2013, Geoelectrical monitoring of the slope of an embankment using a large-scale rainfall simulator: Proceedings of the 11th SEGJ International Symposium, 228–231.
Takeshita, Y., Konishi, K., Inoue, M., and Kohno, I., 1995, Studies on neural-network approach to the determination of unsaturated soil hydraulic properties by electric prospecting method: Proceedings of the 92nd SEGJ Conference, 240–243 [in Japanese].