Simulation of saltwater movement in the Upper Floridan aquifer in the Savannah, Georgia-Hilton Head Island, South Carolina, area, predevelopment-2004, and projected movement for 2000 pumping conditions

A digital model was developed to simulate ground-water flow and solute transport for the Upper Floridan aquifer in the Savannah, Georgia–Hilton Head Island, South Carolina, area. The model was used to (1) simulate trends of saltwater intrusion from predevelopment to the present day (1885–2004), (2) project these trends from the present day into the future, and (3) evaluate the relative influence of different assumptions regarding initial and boundary conditions and physical properties. The model is based on a regional, single-density ground-water flow model of coastal Georgia and adjacent parts of South Carolina and Florida.

Variable-density ground-water flow and solute transport were simulated using the U.S. Geological Survey finite-element, variable-density solute-transport simulator SUTRA, 1885–2004. The model comprises seven layers: the surficial aquifer system, the Brunswick aquifer system, the Upper Floridan aquifer, the Lower Floridan aquifer, and the intervening confining units.

The model was calibrated to September 1998 water levels, for single-density freshwater conditions, then refined using variable density and chloride concentration to give a reasonable match to the trend in the chloride distribution in the Upper Floridan aquifer inferred from field measurements of specific conductance made during 2000, 2002, 2003, and 2004. The model was modified to simulate solute transport by allowing saltwater to enter the system through localized areas near the northern end of Hilton Head Island, at Pinckney Island, and near the Colleton River, and was calibrated to match chloride concentrations inferred from field measurements of specific conductance. This simulation is called the "Base Case." Water-level residuals ranged from –5.3 to 23.4 feet for September 1998 conditions and single-density freshwater conditions. When chloride transport was simulated, water-level residuals ranged from –12.5 to 23.3 feet. The simulated chloride distribution captures the general trends in the field data. Chloride transport is sensitive to the permeabilities assigned to the confining units in the source areas and the porosity assigned to the Upper Floridan aquifer.

Results of the study indicate that

1. if present-day (year 2000) pumping conditions are maintained, plumes of saltwater in the Upper Floridan aquifer will continue to expand and move toward Savannah and across Hilton Head Island;
2. the rate of movement of the 250-mg/L (milligram per  liter) isochlor toward Savannah is between 144  feet per year and 190 feet per year and that the 250-mg/L isochlor could reach the pumping center at Savannah in 800 years;
3. if effective porosities are lower than those used in the model, as is likely, higher rates of solute transport would result; and
4. plumes may have occurred along the northern shore of Hilton Head Island before substantial development began in the mid-1960s, and lesser amounts of intrusion may have already occurred prior to the onset of pumping during 1885.

Model limitations include uncertainty in (1) field data, (2) the conceptual model, (3) the physical properties and representation of the hydrogeologic framework, and (4) uncertainty in the boundary and initial conditions. Results of simulations projected far into the future must be interpreted with caution because they are based on an assumed future pumping distribution and fixed boundary conditions, and because these conditions may differ substantially from those for which the model is calibrated.