Simulation of Regional Karst Aquifer System and Assessment of Groundwater Resources in Manatí-Vega Baja, Puerto Rico


The North Coast karst aquifer system of Puerto Rico, the most productive aquifer of the island, is a vital water source for drinking water and local ecosystems. High freshwater demands alter the coastal groundwater system that impacts both human populations and coastal ecosystems of the island. To predict how this system might respond to rainfall events and high pumping demands, we used the equivalent porous medium (EPM) technique to develop a three-dimensional groundwater flow model to estimate hydrogeological parameters and assess groundwater resources in the Manatí-Vega Baja karst aquifer. The approach is based on the hypothesis that the simplified EPM approach will reproduce groundwater hydrodynamics in this complex karst environment. The steady-state model was calibrated with trial and error and parameter estimation methods using an observed groundwater table of 1995 (r = 0.86, p < 0.0001, n = 39). The large-scale simulation suggested that groundwater flow roughly follows the elevation slope (i.e. south to north). Calibrated hydraulic conductivities range from 0.5 to 86 m/d, whereas the hydro-geologic data strongly suggest higher permeability in the middle karst section of the study area. The transient model adequately estimates the observed groundwater fluctuations in response to rainfall events from 1980 until 2014. The transient results indicate that the conceptual model accuracy is more acceptable with a mean error (ME) of -0.132 m, mean absolute error (MAE) of 0.542 m and root mean square (RMSE) error of 0.365 m. The results of water budget simulation show that the total recharge satisfies the total groundwater withdrawal rate in the past, but continuous closure of more contaminated wells causes groundwater levels to increase in the future. The results indicate that the assumption of applicability of EPM approach is sustained and supported by measured data in the study area. Taking future water demands into account, this model could be applied further to predict the changes of groundwater levels and mass balance under different exploitation scenarios.

Share and Cite:

Maihemuti, B. , Ghasemizadeh, R. , Yu, X. , Padilla, I. and Alshawabkeh, A. (2015) Simulation of Regional Karst Aquifer System and Assessment of Groundwater Resources in Manatí-Vega Baja, Puerto Rico. Journal of Water Resource and Protection, 7, 909-922. doi: 10.4236/jwarp.2015.712075.

Conflicts of Interest

The authors declare no conflicts of interest.


[1] Regli, C., Rauber, M. and Huggenberger, P. (2003) Analysis of Aquifer Heterogeneity within a Well Capture Zone, Comparison of Model Data with Field Experiments: A Case Study from the River Wiese, Switzerland. Aquatic Sciences, 65, 111-128.
[2] Angelini, P. and Dragoni, W. (1997) The Problem of Modeling Limestone Springs: The Case of Bagnara (North Apennines, Italy. Ground Water, 35, 612-618.
[3] Kaufmann, G. and Braun, J. (2000) Karst Aquifer Evolution in Fractured, Porous Rocks. Water Resources Research, 36, 1381-1391.
[4] Ford, D.C. and Williams, P.W. (2007) Karst Hydrogeology and Geomorphology. Wiley, Chichester, 576.
[5] Király, L. (2003) Karstification and Groundwater Flow. Speleogenesis and Evolution of Karst Aquifers, 1, 1-26.
[6] Hu, C.H., Hao, Y.H., Yeh, T.C.J., Pang, B. and Wu, Z.N. (2008) Simulation of Spring Flows from a Karst Aquifer with an Artificial Neural Network. Hydrological Processes, 22, 596-604.
[7] Bakalowicz, M. (2005) Karst Groundwater: A Challenge for New Resources. Hydrogeology Journal, 13, 148-160.
[8] Faulkner, J., Hu, B.X., Kish, S. and Hua, F. (2009) Laboratory Analog and Numerical Study of Groundwater Flow and Solute Transport in a Karst Aquifer with Conduit and Matrix Domains. Journal of Contaminant Hydrology, 110, 34-44.
[9] Greene, E.A. (1992) Hydrologic Properties of the Madison Aquifer System in the Western Rapid City Area, South Dakota. USGS Water Resources Investigation, 56, 93-108.
[10] Teutsch, G. (1993) An Extended Double-Porosity Concept as a Practical Modeling Approach for a Karstified Terraine. Hydrogeological Processes in Karst Terraines. Proceedings of the Antalya Symposium and Field Seminar, International Association of Hydrological Sciences, Wallingford, October 1990, 281-292.
[11] Keeler, R.R. and Zhang, Y.K. (1997) Modeling of Groundwater Flow in a Fractured-Karst Aquifer in the Big Springs Basin, Iowa. GSA Abs with Programs, 29, 25.
[12] Larocque, M., Banton, O., Ackerer, P. and Razack, M. (1999) Determining Karst Transmissivities with Inverse Modeling and an Equivalent Porous Media. Ground Water, 37, 897-903.
[13] Bauer, S., Liedl, R. and Sauter, M. (2003) Modeling of Karst Aquifer Genesis: Influence of Exchange Flow. Water Resources Research, 39, 1285.
[14] Abusaada, M. and Sauter, M. (2012) Studying the Flow Dynamics of a Karst Aquifer System with an Equivalent Porous Medium Model. Groundwater, 51, 641-650.
[15] Scanlon, B.R., Mace, R.E., Barret, M.E. and Smith, B. (2003) Can We Simulate Regional Groundwater Flow in a Karst System Using Equivalent Porous Media Models? Case Study, Barton Springs Edwards Aquifer, USA. Journal of Hydrology, 276, 137-158.
[16] Pankow, J.F., Johnson, R.L., Hewetson, J.P. and Cherry, J.A. (1986) An Evaluation of Contaminant Migration Patterns at Two Waste Disposal Sites on Fractured Porous Media in Terms of the Equivalent Porous Medium (EPM) Model. Journal of Contaminant Hydrology, 1, 65-76.
[17] Neuman, S.P. (1987) Stochastic Continuum Representation of Fractured Rock Permeability as an Alternative to the REV and Fracture Network Concepts. In: Custodio, E., Gurgui, A. and Lobo-Ferreira, J.P., Eds., NATO Advanced Workshop on Advances in Analytical and Numerical Groundwater Flow and Quality Modelling, NATO ASI Series, Series C: Mathematical and Physical Sciences, Volume 224, Reidel Publications, Dordrecht, 331-362.
[18] Huntoon, P.W. (1994) Is It Appropriate to Apply Porous Media Groundwater Circulation Models to Karstic Aquifers? In: El-Kadi, A.I., Ed., Groundwater Models for Resources Analysis and Management, CRC/Lewis Publishers, Boca Raton, 339-358.
[19] Lugo, A.E. and Helmer, E. (2004) Emerging Forests on Abandoned Land: Puerto Rico’s New Forests. Forest Ecology and Management, 190, 145-161.
[20] Monroe, W.H. (1976) The Karst Landforms of Puerto Rico. US Geological Survey Professional Paper 899, 68.
[21] Troester, J.W. (1992) The Northern Karst Belt of Puerto Rico: A Humid Tropical Karst. In: Back, W., Herman, J.S., and Paloc, H., Eds., Hydrogeology of Selected Karst Regions, Verlag Heinz Heise, Hannover, 475-486.
[22] Lugo, A.E., Castro, L.M., Vale, A., et al. (2001) Puerto Rican Karst—A Vital Resource. US Department of Agriculture, Forest Service, General Technical Report WO-65, 100.
[23] Monroe, W.H. (1980) Geology of the Middle Tertiary Formations of Puerto Rico. US Geological Survey Professional Paper 953, 93. (1 pl., Scale 1:250,000)
[24] Renken, R.A., Ward, W.C., Gill, I.P., Gómez, G.F., Rodríguez, J.M., et al. (2002) Geology and Hydrogeology of the Caribbean Islands Aquifer System of the Commonwealth of Puerto Rico and the US Virgin Islands. US Geological Survey Professional Paper 1419, 1-139.
[25] Yu, X., Ghasemizadeh, R., Padilla, I., Irizarry, C., Kaeli, D. and Alshawabkeh, A. (2015) Spatiotemporal Changes of CVOC Concentrations in Karst Aquifers: Analysis of Three Decades of Data from Puerto Rico. Science of The Total Environment, 511, 1-10.
[26] Heisel, J.E., González, J.R. and Cruz, C. (1983) Analog Model Analysis of the North Coast Limestone Aquifer, Puerto Rico. US Geological Survey Open-File Report 82-52, 49.
[27] Gomez-Gomez, F. and Torres-Sierra, H. (1988) Hydrology and Effects of Development on the Water-Table Aquifer in the Vega Alta Quadrangle, Puerto Rico. US Geological Survey Water-Resources Investigations Report 87-4105, 54.
[28] Torres-González, S., Planert, M. and Rodríguez, J.M. (1996) Hydrogeology and Simulation of Ground-Water Flow in the Upper Aquifer of the Río Camuy to Río Grande de Manatí Area, Puerto Rico. US Geological Survey Water-Resources Investigations Report 95-4286, 102.
[29] Bennett, G.D. and Giusti, E.V. (1972) Ground Water in the Tortuguero Area, Puerto Rico, as Related to Proposed Harbor Construction. US Geological Survey Water-Resources Bulletin 10, 25.
[30] Quinones-Marquez, F. and Fuste, L.A. (1978) Limnology of Laguna Tortuguero, Puerto Rico. US Geological Survey Water-Resources Investigations Report 77-122, 86.
[31] McDonald, M.G. and Harbaugh, A.W. (1988) A Modular Three-Dimensional Finite-Difference Groundwater Flow Model. US Geological Survey, Techniques of Water-Resources Investigations, Book 6, Chapter A1, USGS, Reston, 586.
[32] Bennett, G.D. and Giusti, E.V. (1976) Ground Water in the Tortuguero Area, Puerto Rico, as Related to Proposed Harbor Construction. US Geological Survey Water-Resources Bulletin 10, 25.
[33] Doherty, J., Brebber, L. and Whyte, P. (1994) PEST: Model-Independent Parameter Estimation. Watermark Computing, Brisbane.
[34] Weiss, M. and Gvirtzman, H. (2007) Estimating Ground Water Recharge Using Flow Models of Perched Karstic Aquifers. Ground Water, 45, 761-773.
[35] Gregory, S.C. (2001) Simulation of Flow in the Upper North Coast Limestone Aquifer, Manatí-Vega Baja Area, Puerto Rico. Water-Resources Investigations Report, San Juan, Puerto Rico, 4266.
[36] Ghasemizadeh, R., Hellweger, F., Butscher, C., Padilla, I., Vesper, D., Field, M. and Alshawabkeh, A. (2012) Review: Groundwater Flow and Transport Modeling of Karst Aquifers, with Particular Reference to the North Coast Limestone Aquifer System of Puerto Rico. Hydrogeology Journal, 20, 1441-1461.
[37] Moujabber, M.E., Samra, B.B., Darwish, T. and Atallah, T. (2006) Comparison of Different Indicators for Groundwater Contamination by Seawater Intrusion on the Lebanese Coast. Water Resources Management, 20, 161-180.
[38] El Moujabber, M., Atallah, T., Darwish, T. and Bou Samra, B. (2004) Monitoring of Groundwater Salination by Seawater Intrusion on the Lebanese Coast. Lebanese Science Journal, 5, 21-36.

Copyright © 2023 by authors and Scientific Research Publishing Inc.

Creative Commons License

This work and the related PDF file are licensed under a Creative Commons Attribution 4.0 International License.