Abstract
In coastal lowland plains, increased water demand on a limited water resource has resulted in declining groundwater levels, land subsidence and saltwater encroachment. In southwestern Kyushu, Japan, a sinking of the land surface due to over pumping of groundwater has long been recognized as a problem in the Shiroishi lowland plain. In this paper, an integrated model was established for the Shiroishi site using the modular finite difference groundwater flow model, MODFLOW, by McDonald and Harbaugh (1988) and the modular three-dimensional finite difference groundwater solute transport model, MT3D, by Zheng (1990) to simulate groundwater flow hydraulics, land subsidence, and solute transport in the alluvial lowland plain. Firstly, problems associated with these groundwater resources were discussed and then the established model was applied. The simulated results show that subsidence rapidly occurs throughout the area with the central prone in the center part of the plain. Moreover, seawater intrusion would be expected along the coast if the current rates of groundwater exploitation continue. Sensitivity analysis indicates that certain hydrogeologic parameters such as an inelastic storage coefficient of soil layers significantly contribute effects to both the rate and magnitude of consolidation. Monitoring the present salinization process is useful in determining possible threats to fresh groundwater supplies in the near future. In addition, the integrated numerical model is capable of simulating the regional trend of potentiometric levels, land subsidence and salt concentration. The study also suggests that during years of reduced surface-water availability, reduction of demand, increase in irrigation efficiency and the utilization of water exported from nearby basins are thought to be necessary for future development of the region to alleviate the effects due to pumping.
Similar content being viewed by others
References
Bakker M, Kraemer SR, Lange WJ, Strack ODL (2000) Analytic element modeling of coastal aquifers. EPA/600/R-99/110, U.S. Environmental Protection Agency, Washington, DC
Goldberg ED (1963) Chemistry: the oceans as a chemical system. In: Hill HM (ed) Composition of a sea water, comparative and descriptive oceanography of the sea, vol 2. Wiley, New York
Guo W, Langevin CD (2002) User’s Guide to SEAWAT: a computer program for simulation of three-dimensional variable-density ground-water flow. Techniques of Water-Resources Investigations 6-A7, USGS, Reston, Virginia
Hachiya Y, Miura N, Sakai S (1988) Characteristics of land subsidence in Saga Plain. In: Miura M, Bergado DT (eds) Improvement of soft ground: design, analyses, and current researches. Asian center for soil improvement and geosynthetics, Asian institute of technology, Bangkok, Thailand
Holzbecher E (1998) Modeling density-driven flow in porous media, principles, numerics, software. Springer, Berlin Heidelberg New York, 286 pp
Larson JK, Basagaoglu H, Marino MA (2001) Prediction of optimal safe groundwater yield and land subsidence in the Los Banos-Kettleman City area, California, using a calibrated numberical simulation model. Journal of Hydrology, 242:79–102
Leake SA, Lilly MR (1997) Documentation of a computer program (FHB1) for assignment of transient specified-flow and specified-head boundaries in applications of the modular finite-difference ground-water flow model (MODFLOW), Open-File Report 97-571, U.S. Geological Survey, Reston, Virginia, 50 pp
Leake SA, Prudic DE (1991) Documentation of a computer program to simulate aquifer-system compaction using the modular finite-difference groundwater flow model. Techniques of Water Resources Investigation, US Geological Survey, Reston, Virginia
McDonald MG, Harbaugh AW (1988) A modular three-dimensional finite difference groundwater flow model. Techniques of Water Resources Investigations, US Geological Survey, Reston, Virginia
Miura N, Taesiri Y, Sakai A (1988) Land subsidence and its influence to geotechnical aspect in Saga plain. Proc. of the Int. Symposium on Shallow Sea and Lowland, October 1988, Saga Univ., Saga, Japan
Onitsuka K (1988) Mechanical properties of the very sensitive Ariake clay. Proc. of the Int. symposium on shallow sea and lowland, October 1988, Saga Univ., Saga, Japan
Sakai A (2001) Land subsidence due to seasonal pumping of groundwater in saga plain, Japan. Lowland Technol Int 3(1):24–39
Shimoyama S, Matsumoto N, Yumura H, Takemura K, Iwao Y, Miura N, Tohno I (1994) Quaternary geology of lowland along the north coast of Ariake bay, west Japan. Science Reports, vol 18, no. 2, Dept. of Earth and Planetary Science, Kyushu Univ., Japan, pp 103–129
Simmons CT, Narayan KA, Wooding RA (1999) On a test case for density-dependent groundwater flow and solute transport models: the salt lake problem. Water Resour Res 35(12):3607-3620
Sorey ML (1978) Numerical modeling of liquid geothermal systems. Prof. Paper 16044-D, U.S. Geological Survey, Reston, Virginia
Stuyfzand PJ (1986) A new hydro-chemical classification of watertypes: principles and application to the coastal dunes aquifer system of the Netherlands. Proceedings of the Ninth Salt Water Intrusion Meeting, May 1986, Delft University of Technology, Delft, the Netherlands, pp 641–655
Tanaka S (1990) Optimum pumping volume of groundwater in Saga and Shiroishi areas. Proc. of the Inter. seminar on geotechnical and water problems in lowland, November 1990, Saga Univ., Saga, Japan
Terzaghi K (1925) Principles of soil mechanics: IV settlement and consolidation of clay. Engineering News Record, McGraw-Hill, New York
Voss CI (1984) SUTRA–A finite element simulation for saturated-unsaturated, fluid density dependent groundwater flow with energy transport or chemically reactive single species solute transport. Water-Resources Investigations Report 84–4369, USGS, Reston, Virginia, 409 pp
Weast RC (1982) Handbook of chemistry and physics, 63rd edn. CRC, Boca Raton, Florida
Wilson AM, Gorelick S (1996) The effect of pulsed pumping on land subsidence in the Santa Clara valley, California. J Hydrol 174:375–396
Zheng Ch (1990) MT3D-A modular three-dimensional transport model for simulation of advection, dispersion and chemical reactions of contaminants in groundwater systems. Report to the U.S. Environmental Protection Agency, Ada, Oklahoma, 170 pp
Zhou G, Esaki T, Mori J (2003) GIS-based spatial and temporal prediction system development for regional land subsidence hazard mitigation. Environ Geol 44:665–678
Acknowledgements
The authors sincerely are thankful for the reviews by editing committee of the Journal of Environmental Geology, and for their valuable comments. This research was supported by the Ministry of Education, Science, Sport, and Culture (Monbusho), Japan.
Author information
Authors and Affiliations
Corresponding author
Rights and permissions
About this article
Cite this article
Don, N.C., Araki, H., Yamanishi, H. et al. Simulation of groundwater flow and environmental effects resulting from pumping. Env Geol 47, 361–374 (2005). https://doi.org/10.1007/s00254-004-1158-1
Received:
Accepted:
Published:
Issue Date:
DOI: https://doi.org/10.1007/s00254-004-1158-1