Mapping Saline Groundwater in a Coastal System Robert Anders1, Danny A. Cronquist1, 2, Lindsey E. White1, 3, and Wesley R. Danskin1 1–U.S. Geological Survey, CaWSC, 4165 Spruance Road, Ste. 200, , CA, 92101 H53B-1524 2–San Diego State University, 5500 Campanile Drive, San Diego, CA 92182 3–Naval Facilities Engineering Command, SW, 1220 Pacific Hwy, San Diego, CA 92132 Major-ion Chemistry Stable Isotopes of H and O

Extensive data, including lithologic information, geophysical logs, and water- Major-ion composition of groundwater samples collected from 9 multiple-well monitoring sites located in coastal San The stable isotopes of hydrogen and oxygen 0 Seawater quality data have been collected in the San Diego, California area since 2001 to Diego County are presented here using Stiff diagrams. As seen in both the West and East cross sections, much of the can be used to identify the different sources of groundwater in the coastal southern California aquifer system is of poor quality, having high sodium and chloride concentrations. recharge throughout the region. These different Global Meteoric Water Line evaluate the suitability of the San Diego Formation and overlying alluvial depos- –10 SDMC and SDNB have particularly poor quality water, requiring scale increases of 50% and 200%, respectively, in order to depict sources are distinguishable by: (1) isotopic values its for use as a drinking-water supply. Lithologic information was compiled from the high concentrations of sodium and chloride at these sites. comprised of a mixture of groundwater and sea- 100 descriptions of drill cuttings collected during the installation of 9 multiple-well water: (2) lighter (more negative) isotopic values –20 Group 1 SDEP SDCC SDOR that are characteristic of recharge originating in monitoring sites to depths of more than 700 meters and was combined with avail- SDHF SDSW Sea the fractured crystalline rock more than 20 km east Seawater Mixing Line able geospatial data sets, surface geologic maps, and various literature references level and upgradient from the coastline; and (3) interme- –30 to wells and outcroppings to generate a three-dimensional geologic framework diate isotopic values that are characteristic of pre- Group 3 cipitation originating in the coastal plain as the model of the coastal San Diego area. Natural gamma-ray and electromagnetic in- –100 –40 duction (EM) logs were obtained prior to the installation of the monitoring wells source of recharge. δ DEUTERIUM, IN PER MIL Group 2 and during subsequent site visits. EM logs are sensitive to changes in lithology and –200 –50 water quality; because the lithology remains constant, repeated EM logs were used to show changes in water-quality due to natural recharge processes, seawater in- –300 –60 trusion, or other processes. Water samples were collected and analyzed for selected –8 –6 –4 –2 0 δOXYGEN-18, IN PERMIL major- and minor-ions, the stable isotopes of hydrogen, oxygen, and strontium, 405030 20 10 0 10 20 30 40 50 –400 milliequivalent per liter and the radiogenic isotope of carbon-14. Through the use of the geophysical logs, 40 30 20 10 0 10 20 30 40 milliequivalent per liter Na/Ca Molar Ratios

water-chemistry data, and cross sections derived from the 3D geologic framework ELEVATION ABOVE OR BELOW SEA LEVEL, IN METERS 40 30 20 10 0 10 20 30 40 –500 milliequivalent per liter N-S East 40 30 20 10 0 10 20 30 40 model, it was possible to map the groundwater quality in the coastal aquifer Stiff Diagrams milliequivalent per liter EXPLANATION The relationship between Na/Ca molar ratios Na+ + K+ Cl− 200 ++ − Ca HCO3 system. In addition, by using sodium-to-calcium molar ratios, in combination Quaternary Oligocene Cretaceous ++ −− Mg SO4 (Otay Formation) and Cl concentrations can be used to distin- –600 (SD Formation) Eocene with the isotopic data, the occurrence of saline groundwater as a result of seawa- 405030 20 10 0 10 20 30 40 50 guish between the processes responsible for milliequivalent per liter the occurrence of saline groundwater in a ter intrusion was distinguishable from groundwater in a previously-saline aquifer 100 150 which has been “freshened” by continental recharge. These ndings indicate mul- SDBP coastal aquifer system. For the San Diego area, SDNB the composition of fresh groundwater is repre- tiple groundwater ow paths exist in the coastal aquifer system. SDMC SDOT “Freshening” of a Sea sented by low Na/Ca molar ratios with low to inter- previously-saline level mediate Cl concentrations. In contrast, increasing aquifer Na/Ca molar ratios without a corresponding in- 100 Background –100 crease in Cl concentrations indicate "freshening" of a previously-saline aquifer by continental recharge. Na/Ca ratio

Over the past decade, the USGS has been conducting a regional assessment of the groundwa- Finally, Na/Ca ratios approaching seawater with in- Na/Ca MOLAR RATIO for seawater ter resources in the San Diego area. The regional assessment of groundwater resources in the San –200 creasing chloride concentrations indicate the 50 10 0 22468 4 6 810 Diego area was designed as an integrated set of five drainage-basin investigations to evaluate the milliequivalent per liter source of saline groundwater to these wells is the suitability of the San Diego Formation and overlying alluvial deposits for use as a drinking water result of seawater intrusion. Seawater intrusion supply. The San Diego Formation is composed of thinly bedded sandstone and conglomerate, which –300 originated as marine and non-marine sediment during the late Pliocene and early Pleistocene, ranges 10 in thickness from about 30 meters to more than 250 meters, and is overlain by unconsolidated Quater- 10 100 1,000 10,000 100,000 –400 CHLORIDE, IN MILLIGRAMS PER LITER nary alluvium of the Linda Vista and Bay Point Formations. An integral part of the regional assess- 40 30 2010 100 20 30 40 ment has been the installation of 9 multiple-well monitoring sites, each with as many as 6 monitoring milliequivalent per liter 1,200 600 200 0 200 600 1,200 140 100 60 200 20 60 100 140 30 20 10 0 10 20 30 wells, to depths of more than 700 meters. ELEVATION ABOVE OR BELOW SEA LEVEL, IN METERS N-S West milliequivalent per liter milliequivalent per liter milliequivalent per liter –500 Stiff Diagrams Changes in Geophysical Logs and Cl Concentrations EXPLANATION Na+ + K+ Cl− ++ − Ca HCO3 Additionally, the multiple-depth monitoring-well sites have been equipped with real-time, water-level ++ −− Mg SO4 recording equipment and the data is available via the project website. (For more information see- –600 Gamma and electromagnetic induction (EM) logs were obtained prior to the installation of the monitoring wells, and during http://ca.water.usgs.gov/sandiego) 0 1,000 2,000 3,000 4,000 5,000 6,000 7,000 8,000 9,000 10,000 11,000 12,000 13,000 14,000 15,000 subsequent site visits. A natural gamma tool is designed to measure the total intensity of gamma-ray emissions from the for- DISTANCE ALONG SECTION, IN METERS mation. EM logs are sensitive to changes in lithology and water quality; because the lithology remains constant repeated EM logs show changes in water-quality due to natural recharge processes, seawater intrusion, or other processes. Recent in- Uncorrected Groundwater Age creases in EM conductivity and Cl concentrations for two multiple-well monitoring sites (SDNB and SDMC) located along the 3D Geologic Framework Model coast indicate seawater intrusion is occurring within these zones. In contrast, decreases in EM conductivity and Cl concen- Carbon-14 was used to derive the uncorrected groundwater ages from calculations using the percent modern carbon trations indicate zones of “freshening” of a previously-saline aquifer by recently recharged water. Lithologic information was compiled from descriptions of drill cuttings collected at each measured in the various samples. The cross sections show a general trend which depicts groundwater in the Quaternary allu- CHANGE IN ELECTROMAGNETIC CONDUCTIVITY, IN MILLIMHOS PER METER borehole and from observations recorded during drilling. These observations provided direct, vium as less than 10,000 years before present, groundwater in the San Diego Formation as generally less than 20,000 years −200 0 200 400 600 −50 0 50 100 150 SDNB 0 0 SDMC before present, and groundwater in the Otay Formation and the Eocene layer as between 25,000 and 45,000 years before present. 320 reliable data for the geologic boundaries of the subsurface. Additional data from wells, borings and 60 100 100 seismic shot points provided age and depth information. These values were coupled with pre- The presence of modern groundwater in only the shallowest monitoring wells indicates most of the groundwater in the coastal −600 700 aquifer system was recharged many thousands of years ago. 200 200 existing GIS data sets, surface geologic maps, drilling and well logs, the aforementioned direct ob- Chloride 86 servations of lithology, and various literature references to wells and outcroppings to provide the 100 300 300 Seawater input data used by RockWorks modeling software to generate a three-dimensional geologic frame- SDCC intrusion SDHF SDEP SDOR 400 400 work model. Two north-south cross-sections were derived from the 3D geologic framework model, SDSW Sea 500 500 one located along the San Diego Bay (West) and the other one located along the coastal plain away level 9,125 13 6,168 3,847 600 600 from the bay (East). The East and West cross sections were used to map the groundwater quality in 9,983 1,942 21,120 −128 Chloride 700 700 the coastal southern California aquifer system. 8,188 −40 –100 13,616 Groundwater 800 “freshening” 800 22,520 8,695 14,062 900 900 12,528 –200 25,102 1,000 1,000

17,276 1,100 1,100 17,548 −100 27,231 24,346 –300 18,475 1,200 1,200 DEPTH, FEET IN BELOW LAND SURFACE DEPTH, FEET IN BELOW LAND SURFACE 1,300 1,300 36,029 26,385 36,701 1,400 1,400 −130 –400 500 32,055 1,500 1,500 −500 0 500 1,000 1,500 2,000 2,500 3,000 −200 −100 0 100 200 300 400 500 29,528 28,824

SDCP ELEVATION ABOVE OR BELOW SEA LEVEL, IN METERS CHANGE IN CHLORIDE CONCENTRATION, IN MILLIGRAMS PER LITER –500 N-S East SDHF Carbon-14 age, SDBP in years before present Quaternary Oligocene Cretaceous Significant Findings (Otay Formation) 44,856 Pliocene –600 (SD Formation) Eocene NC Not collected SDEP 1. A three-dimensional geologic framework model was generated and used to map the groundwater SDNB 100 SDBP quality in the coastal southern California aquifer system. SDSW SDMC SDNB SDMC SDOT Sea 4,852 2. Stiff diagrams are used to depict the water quality of groundwater samples collected from 9 level 2,172 NC 1,024 13,375 multiple-well monitoring sites located in coastal San Diego County. As seen in both the West and East 9,095 1,656 SDCC NC cross sections, much of the groundwater in the coastal southern California aquifer system is of poor 117 116 30’ 116 –100 quality. San Dieguito River SDOT SDOR 16,306 15 Drainage Basin 12,272 5 San Diego River 3. The uncorrected C-14 ages of groundwater in the coastal aquifer system generally correlates with 33 Drainage Basin –200 14,333 14,390 Sweetwater River depth and indicates that most of the groundwater was recharged many thousands of years ago. Drainage Basin 8 27,643 22,960 Otay River –300 4. The stable isotopes of hydrogen and oxygen indicate the different sources of recharge for the ground- Drainage Basin United States 19,202 Mexico water in the San Diego area are distinguishable with more locally recharged groundwater being found 32 30’ River at shallower depths and more distantly recharged groundwater found at greater depths. Drainage –400 Basin 30,202 29,154 20,828 30,363 5. The relationship between Na/Ca molar ratios and Cl concentrations can be used to distinguish be-

ELEVATION ABOVE OR BELOW SEA LEVEL, IN METERS N-S West CALIFORNIA –500 tween the presence of saline groundwater as a result of seawater intrusion and “freshening” of a Carbon-14 age, previously-saline aquifer by continental recharge. in years before present Area of map –600 6. Changes in EM conductivities provide evidence for the occurrence of seawater intrusion and “freshen- 0 1,000 2,000 3,000 4,000 5,000 6,000 7,000 8,000 9,000 10,000 11,000 12,000 13,000 14,000 15,000 DISTANCE ALONG SECTION, IN METERS ing” of a previously-saline aquifer by recently recharged water.