Karst systems have a high degree of heterogeneity and anisotropy, which makes them behave very differently from other aquifers. Slow seepage through the rock matrix and fast flow through conduits and fractures result in a high variation in spring response to precipitation events. Contaminant storage occurs in the rock matrix and epikarst, but contaminant transport occurs mostly along preferential pathways that are typically inaccessible locations, which makes modeling of karst systems challenging. Computer models for understanding and predicting hydraulics and contaminant transport in aquifers make assumptions about the distribution and hydraulic properties of geologic features that may not always apply to karst aquifers. This paper reviews the basic concepts, mathematical descriptions, and modeling approaches for karst systems. The North Coast Limestone aquifer system of Puerto Rico (USA) is introduced as a case study to illustrate and discuss the application of groundwater models in karst aquifer systems to evaluate aquifer contamination.

Karst aquifers, capable of storing and transmitting large amount of water, are the main source of drinking water in many regions worldwide. Their excessive permeability leads to an enhanced vulnerability to retain and spread the contamination accordingly. From sustainability perspective, the environmental, economic and social impacts of karst contamination on water resources management are gaining more attention. In this study, an overview of hydrogeological processes and concepts regarding groundwater flow and contaminant transport in karstic systems is presented, followed by a short discussion on surface water and groundwater interaction. Due to the complexity of karstic systems, different approaches have been developed by researchers for investigating and understanding hydrogeological processes and groundwater behavior in karst which are reviewed herein. Additionally, groundwater contamination issues and the most common and effective remediation techniques in karstic terrains are discussed. Lastly, modeling techniques and remote sensing methods, as beneficial and powerful tools for assessing groundwater flow and contaminant transport in karst terrains, are reviewed and evaluated. In each section, relevant research works conducted for Puerto Rico are discussed and some recommendations are presented to complement the ongoing hydrogeological investigations on this island.

We studied the spatial and temporal distribution patterns of Chlorinated Volatile Organic Compounds (CVOCs) in the karst aquifers in northern Puerto Rico (1982-2013). Seventeen CVOCs were widely detected across the study area, with the most detected and persistent contaminated CVOCs including trichloroethylene (TCE), tetrachloroethylene (PCE), carbon tetrachloride (CT), chloroform (TCM), and methylene chloride (DCM). Historically, 471 (76%) and 319 (52%) of the 615 sampling sites have CVOC concentrations above the detection limit and maximum contamination level (MCL), respectively. The spatiotemporal patterns of the CVOC concentrations showed two clusters of contaminated areas, one near the Superfund site “Upjohn” and another near “Vega Alta Public Supply Wells.” Despite a decreasing trend in concentrations, there is a general northward movement and spreading of contaminants even beyond the extent of known sources of the Superfund and landfill sites. Our analyses suggest that, besides the source conditions, karst characteristics (high heterogeneity, complex hydraulic and biochemical environment) are linked to the long-term spatiotemporal patterns of CVOCs in groundwater.

SUMMARY We studied the fractal scaling behavior of groundwater level fluctuation for various types of aquifers in Puerto Rico using the methods of (1) detrended fluctuation analysis (DFA) to examine the monofractality and (2) wavelet transform maximum modulus (WTMM) to analyze the multifractality. The DFA results show that fractals exist in groundwater fluctuations of all the aquifers with scaling patterns that are anti-persistent (1 < β < 1.5; 1.32 ± 0.12, 18 wells) or persistent (β > 1.5; 1.62 ± 0.07, 4 wells). The multi-fractal analysis confirmed the need to characterize these highly complex processes with multifractality, which originated from the stochastic distribution of the irregularly-shaped fluctuations. The singularity spectra of the fluctuation processes in each well were site specific. We found a general elevational effect with smaller fractal scaling coefficients in the shallower wells, except for the Northern Karst Aquifer Upper System. High spatial variability of fractal scaling of groundwater level fluctuations in the karst aquifer is due to the coupled effects of anthropogenic perturbations, precipitation, elevation and particularly the high heterogeneous hydrogeological conditions.