**Numerical Simulation of Groundwater Flow and Solute Transport in a Karst Aquifer with Conduits**

Bill X. Hu and Zexuan Xu

[34] Lam L, Fredlund D. G, Barbour S. L. Transient seepage model for saturated/unsaturated soil systems: a geotechnical engineering approach. Canadian Geotechnical Journal.

[35] Auvinet G, López-Acosta N. P. Rapid drawdown condition in submerged slopes. In: Mexican Society for Geotechnical Engineering (SMIG) and Institute of Engineering (UNAM), editors. 15 Presentations of Friends and Colleagues in the Tribute to Jesús Alberro. 1st ed. Mexico DF: Mexican Society for Geotechnical Engineering (SMIG);

[36] Fredlund M. Theory manual: a knowledge-based databases system for saturated/ unsaturated soil properties. 1st ed. Saskatoon, Saskatchewan, Canada: SoilVision

[37] Mendoza-Promotor J. A. Contribution to the analyses of water flow in unsaturated soils [thesis]. Mexico, D.F: School of Engineering and Architecture (IPN); 2014. 172 p. (in

[38] Brooks R. H, Corey A. T. Hydraulic properties of porous media. Hydrology Papers.

[39] Van Genuchten M. Th. A closed form-equation for predicting the hydraulic conduc‐ tivity of unsaturated soils. Journal of Soil Science Society of America. 1980;44:892–898.

[40] Fredlund D. G, Xing A. Equations for the soil-water characteristic curve. Canadian

[41] Aubertin M, Mbonimpa M, Bussière B, Chapuis R. P. A model to predict the water retention curve from basic geotechnical properties. Canadian Geotechnical Journal.

[42] Kovacs G. Seepage hydraulics (developments in water science). 1st ed. House of the Hungarian Academy of Sciences, Budapest, Hungary: Elsevier Scientific Publishing

[43] Sánchez-Solís M. A. Hydrogeomechanical behaviour of levees under transient-state flow conditions [thesis]. Mexico DF: Faculty of Engineering, National University of

1987;24:565–580.

114 Groundwater - Contaminant and Resource Management

2010. 167–189 pp. (in Spanish).

Geotechnical Journal. 1994;31(3):521–532.

Mexico (UNAM); 2013. 153 p. (in Spanish).

Systems Ltd; 2011. 217 p.

Spanish).

1964;3:1–36.

2003;40(6):1104–1122.

Company; 1981. 729 p.

Additional information is available at the end of the chapter

http://dx.doi.org/10.5772/63766

#### **Abstract**

Numerical simulation of groundwater flow and solute transport in a dual-permeabili‐ ty karst aquifer is a challenging issue, since groundwater flow in karst conduit network can become non-darcian even turbulent flow. The discrete-continuum model is a relatively new modeling method, which accounts for turbulent and laminar flow in karst aquifer. MODFLOW-CFP (Conduit Flow Process) is compared to the MODFLOW, a numerical code based on Darcy law, to evaluate the accuracy in a sub-regional scale karst aquifer. MODFLOW-CFP is more accurate than the MODFLOW when compar‐ ing the head simulation results with field measurements. After that, the CFPv2 and UMT3D numerical models are applied in the WKP to establish a sub-regional scale model to simulate chloride transport processes in the last four decades, and to predict contamination development. Numerical simulation results indicate sprayfields are the major chloride source in the study region. Conduit networks significantly control solute transport and contaminant distribution in the study region. Chloride transports through conduits rapidly and spread to several large contamination plumes in a short period. Chloride concentration started to increase in 1980s due to the operation of sparyfield. Solute transport simulation results by discrete-continuum models are more accurate because of the precise description of conduit network.

**Keywords:** karst aquifer, numerical modeling, discrete‐continuum numerical model, conduit flow process, Solute Transport
