USGS Science⌕ Search

SEARCH · USGS Science

Results for “Geophysical Investigations Map”

Search indexed USGS publications on groundwater, aquifers, geologic maps, mineral resources and earthquakes. Explore source records by subject and place.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

1,201 records · Page 67Linked to original sources
USGS2000

Hydrogeology of the regional aquifer near Flagstaff, Arizona, 1994-97

Sandstones, siltstones, and limestones that are Pennsylvanian to Permian in age underlie the southern part of the Colorado Plateau near Flagstaff, Arizona, and contain a complex regional aquifer that has become increasingly important as a source of water for domestic, municipal, and recreational uses. Ground-water flow in the regional aquifer is poorly understood in this area because (1) depth of the aquifer limits exploratory drilling and testing and (2) the geologic structure increases the complexity of the aquifer characteristics and the ground-water flow system. Four methods were used to improve the understanding of the hydrogeology of the regional aquifer near Flagstaff. • Remote-sensing techniques and geologic mapping provided data to identify many structural features that indicate a more complex structural environment and history than previously realized. • Data from surface-geophysical techniques that included ground-penetrating radar, seismic reflection and seismic refraction, and square-array resistivity, verified that some of the geologic structures expressed at land surface propagate deep into the subsurface and through the principal water-bearing zones of the regional aquifer at near-vertical angles. • A well and spring inventory, borehole-geophysical methods, and well and aquifer tests provided additional information relating aquifer and ground-water flow characteristics to geologic structure. • Water-chemistry data, which included major ion, nutrient, trace-element, and radioactive and stable-isotope analyses, provided an independent means of verifying the hydrogeologic characteristics of the aquifer and were used to determine recharge and discharge areas, groundwater movement, and ground-water age. Ground-water recharge occurs throughout the area but is greatest at higher altitudes where precipitation is greater and in areas where heavily fractured rock units of the aquifer are exposed. The estimated annual average recharge to the regional aquifer in the study area is about 290,000 acre-feet. Ground water flows laterally and vertically through pore spaces in the rock and along faults and other fractures from high-altitude areas in the southern part of the study area to regional drains north of the study area along the Little Colorado and Colorado Rivers, and to drains south of the study area along Oak Creek and the Verde Valley. Ground-water discharge in these areas—about 400,000 acre-feet per year—exceeds the annual recharge to the aquifer in the Flagstaff area, but ground water from areas outside the study area contributes to this discharge as well. The saturated thickness of the regional aquifer averages about 1,200 feet, and the amount of water in storage could be as much as 4,800,000 acrefeet, or about 10 percent of the total volume of the aquifer. The quality of water in the regional aquifer in terms of dissolved-solids concentrations is good for most uses throughout the area. Dissolvedsolids concentrations generally are less than 500 milligrams per liter. Water in the regional aquifer is primarily a calcium magnesium bicarbonate type. In some areas near the Rio de Flag, the water has significant nitrate and chloride components, which indicate direct recharge in these areas from the Rio de Flag. Oxygen and deuterium data indicate a common recharge source for water in the aquifer and that some sites receive recharge from surface waters where evaporation has occurred. Estimated carbon-14 ages and tritium activities indicate ground-water ages from less than 200 years in the Lake Mary area to more than 5,000 years in the Wupatki area. The regional aquifer is heterogeneous and anisotrophic and has a complex ground-water flow system. The most productive water-bearing material tends to be fine- to medium-grained sandstones, and ground-water flow and potential well yields are related to geologic structure. Fracturing associated with structural deformation increases recharge locally and also increases the potential for high well yields. Surface-geophysical techniques provided information on the orientation of high-angle, deep-seated structure in the saturated zone. Borehole-geophysical data identified horizontal to near-horizontal fractures as significant components of the fracture-flow system not apparent in the surface-geophysical data. Structural features that strike northwest appear to be areas that have the greatest potential for high well yields. A north-northeastward-striking structure may be just as promising, but additional data are needed to verify this relation.

Arizona↗
USGS2009

Sand resources, regional geology, and coastal processes of the Chandeleur Islands Coastal System: An evaluation of the Breton National Wildlife Refuge

Breton National Wildlife Refuge, the Chandeleur Islands chain in Louisiana, provides habitat and nesting areas for wildlife and is an initial barrier protecting New Orleans from storms. The U.S. Geological Survey (USGS) in partnership with the University of New Orleans Pontchartrain Institute for Environmental Sciences undertook an intensive study that included (1) an analysis of island change based on historical maps and remotely sensed shoreline and topographic data; (2) a series of lidar surveys at 3- to 4-month intervals after Hurricane Katrina to determine barrier island recovery potential; (3) a discussion of sea level rise and effects on the islands; (4) an analysis of sea floor evolution and sediment dynamics in the refuge over the past 150 years; (5) an assessment of the local sediment transport and sediment resource availability based on the bathymetric and subbottom data; (6) a carefully selected core collection effort to groundtruth the geophysical data and more fully characterize the sediments composing the islands and surrounds; (7) an additional survey of the St. Bernard Shoals to assess their potential as a sand resource; and (8) a modeling study to numerically simulate the potential response of the islands to the low-intensity, intermediate, and extreme events likely to affect the refuge over the next 50 years. Results indicate that the islands have become fragmented and greatly diminished in subaerial extent over time: the southern islands retreating landward as they reorganize into subaerial features, the northern islands remaining in place. Breton Island, because maintenance of the Mississippi River-Gulf Outlet (MRGO) outer bar channel requires dredging, is deprived of sand sufficient to sustain itself. Regional sediment transport trends indicate that large storms are extremely effective in transporting sand and controlling the shoreline development and barrier island geometry. Sand is transported north and south from a divergent zone near Monkey Bayou at the southern end of the Chandeleur Islands. Numerical simulation of waves and sediment transport supports the geophysical results and indicates that vast areas of the lower shoreface are affected and are undergoing erosion during storm events, that there is little or no fair weather mechanism to rework material into the littoral system, and that as a result, there is a net loss of sediment from the system. Lidar surveys revealed that the island chain immediately after Hurricane Katrina lost about 84 percent of its area and about 92 percent of its prestorm volume. Marsh platforms that supported the islands’ sand prior to the storm were reduced in width by more than one-half. Repeated lidar surveys document that in places the shoreline has retreated about 100 m under the relatively low-energy waves since Hurricanes Katrina and Rita; however, this retreat is nonuniform. Recent high-resolution geophysical surveys of the sea floor and subsurface within 5–6 km of the Chandeleur Islands during 2006 and 2007 show that, in addition to the sand that is rebuilding portions of the island chain, a large volume of sand is contained in Hewes Point, in an extensive subtidal spit platform that has formed at the northern end of the Chandeleur Islands. Hewes Point appears to be the depositional terminus of the alongshore transport system. In the southern Chandeleurs, sand is being deposited in a broad tabular deposit near Breton Island called the southern offshore sand sheet. These two depocenters account for approximately 70 percent of the estimated sediment volume located in potential borrow sites. An additional large potential source of sand for restoration lies in the St. Bernard Shoals, which are estimated to contain approximately 200 × 10 6 m 3 of sand. Successful restoration planning for the Breton National Wildlife Refuge should mimic the natural processes of early stages of barrier island evolution including lateral transport to the flanks of the island chain from a centralized sand source that will ultimately enhance the ability of the islands to naturally build backbarrier marsh, dunes, and a continuous sandy shoreline. Barrier island sediment nourishment should be executed with the understanding that gulf shoreline erosion is inevitable but that island area can be maintained and enhanced during retreat (thus significantly prolonging the life of the island chain) with strategic sand placement.

Louisiana↗
USGS2003

Ground-water conditions and studies in Georgia, 2001

The U.S. Geological Survey (USGS) collects ground-water data and conducts studies to monitor hydrologic conditions, to better define ground-water resources, and address problems related to water supply and water quality. Data collected as part of ground-water studies include geologic, geophysical, hydraulic property, water level, and water quality. A ground-water-level network has been established throughout most of the State of Georgia, and ground-water-quality networks have been established in the cities of Albany, Savannah, and Brunswick and in Camden County, Georgia. Ground-water levels are monitored continuously in a network of wells completed in major aquifers of the State. This network includes 17 wells in the surficial aquifer, 12 wells in the upper and lower Brunswick aquifers, 73 wells in the Upper Floridan aquifer, 10 wells in the Lower Floridan aquifer and underlying units, 12 wells in the Claiborne aquifer, 1 well in the Gordon aquifer, 11 wells in the Clayton aquifer, 11 wells in the Cretaceous aquifer system, 2 wells in Paleozoic-rock aquifers, and 7 wells in crystalline-rock aquifers. In this report, data from these 156 wells were evaluated to determine whether mean-annual ground-water levels were within, below, or above the normal range during 2001, based on summary statistics for the period of record. Information from these summaries indicates that water levels during 2001 were below normal in almost all aquifers monitored, largely reflecting climatic effects from drought and pumping. In addition, water-level hydrographs for selected wells indicate that water levels have declined during the past 5 years (since 1997) in almost all aquifers monitored, with water levels in some wells falling below historical lows. In addition to continuous water-level data, periodic measurements taken in 52 wells in the Camden County-Charlton County area, and 65 wells in the city of Albany-Dougherty County area were used to construct potentiometric-surface maps for the Upper Floridan aquifer. Ground-water quality in the Upper Floridan aquifer is monitored in the cities of Albany, Savannah, and Brunswick and in Camden County; and monitored in the Lower Floridan aquifer in the Savannah and Brunswick areas. In the Albany area since 1998, nitrate concentrations in the Upper Floridan aquifer have increased in 4 of the 11 wells monitored, and in 1 well, concentrations were above the U.S. Environmental Protection Agency's (USEPA) 10 milligrams per liter (mg/L) drinking-water standard. In the Savannah area, chloride concentration in water from four wells in the Upper Floridan aquifer showed no appreciable change during 2001, remaining within the USEPA 250 mg/L drinking-water standard; in seven wells completed in the Lower Floridan aquifer and in underlying zones, the chloride concentration remained above the drinking-water standard, with one well showing an increase over previous years. In the Brunswick area, water samples from 66 wells completed in the Upper or Lower Floridan aquifers were collected during June 2001 and analyzed for chloride. A map showing chloride concentrations in the Upper Floridan aquifer during June 2001 indicates that concentrations remained above USEPA drinking-water standards across a 2-square-mile area. In the north Brunswick area, chloride concentrations in the Upper Floridan aquifer continued to increase, whereas in the south Brunswick area, concentrations continued to decrease. In the Camden County area, chloride concentrations in six wells completed in the Upper Floridan aquifer remained within drinking-water standards. With the exception of one well, concentrations remained the same and were below 40 mg/L. In one well, concentrations showed a sharp decline during 2001, but remained above 130 mg/L. Ongoing studies during 2001 include evaluation of agricultural chemicals in shallow ground water in southwestern Georgia; evaluation of saltwater intrusion and water-level and water-quality m

Georgia↗
USGS2004

Stratigraphy and tectonic history of the Tucson Basin, Pima County, Arizona, based on the Exxon state (32)-1 well

The Tucson Basin is a relatively large late Cenozoic extensional basin developed in the upper plate of the Catalina detachment fault in the southern Basin and Range Province, southeastern Arizona. In 1972, Exxon Company, U.S.A., drilled an exploration well (Exxon State (32)-1) near the center of the Tucson Basin that penetrated 3,658 m (12,001 ft) of sedimentary and volcanic rocks above granitoid basement. Detailed study of cuttings and geophysical logs of the Exxon State well has led to revision of the previously reported subsurface stratigraphy for the basin and provided new insight into its depositional and tectonic history. There is evidence that detachment faulting and uplift of the adjacent Catalina core complex on the north have affected the subsurface geometry of the basin. The gravity anomaly map of the Tucson Basin indicates that the locations of subbasins along the north-trending axis of the main basin coincide with the intersection of this axis with west-southwest projections of synforms in the adjacent core complex. In other words, the subbasins overlie synforms and the ridges between subbasins overlie antiforms. The Exxon State well was drilled near the center of one of the subbasins. The Exxon well was drilled to a total depth of 3,827 m (12,556 ft), and penetrated the following stratigraphic section: Pleistocene(?) to middle(?) Miocene upper basin-fill sedimentary rocks (0-908 m [0-2,980 ft]) lower basin-fill sedimentary rocks (908-1,880 m [2,980-6,170 ft]) lower Miocene and upper Oligocene Pantano Formation (1,880-2,516 m [6,170-8,256 ft]) upper Oligocene to Paleocene(?) volcanic and sedimentary rocks (2,516-3,056 m [8,256-10,026 ft]) Lower Cretaceous to Upper Jurassic Bisbee Group (3,056-3,658 m [10,026-12,001 ft]) pre-Late Jurassic granitoid plutonic rock (3,658-3,827 m [12,001- 12,556 ft]). Stratigraphy and Tectonic History of the Tucson Basin, Pima County, Arizona, Based on the Exxon State (32)-1 Well The 1,880 m (6,170 ft) of basin-fill sedimentary rocks consist of alluvial-fan, alluvial-plain, and playa facies. The uppermost unit, a 341-m-thick (1,120-ft) lower Pleistocene and upper Pliocene alluvial-fan deposit (named the Cienega Creek fan in this study), is an important aquifer in the Tucson basin. The facies change at the base of the alluvial fan may prove to be recognizable in well data throughout much of the basin. The well data show that a sharp boundary at 908 m (2,980 ft) separates relatively unconsolidated and undeformed upper basin fill from denser, significantly faulted lower basin fill, indicating that there were two stages of basin filling in the Tucson basin as in other basins of the region. The two stages apparently occurred during times of differing tectonic style in the region. In the Tucson area the Pantano Formation, which contains an andesite flow dated at about 25 Ma, fills a syntectonic basin in the hanging wall of the Catalina detachment fault, reflecting middle Tertiary extension on the fault. The formation in the well is 636 m thick (2,086 ft) and consists of alluvial-fan, playa, and lacustrine sedimentary facies, a lava flow, and rock- avalanche deposits. Analysis of the geophysical logs indicates that a K-Ar date of 23.4 Ma reported previously for the Pantano interval of the well was obtained on selected cuttings collected from a rock-avalanche deposit near the base of the unit and, thus, does not date the Pantano Formation. The middle Tertiary volcanic and sedimentary rocks have an aggregate thickness of 540 m (1,770 ft). We obtained a new 40Ar/ 39Ar age of 26.91+0.18 Ma on biotite sampled at a depth of 2,584-2,609 m (8,478-8,560 ft) from a 169-m-thick (554-ft) silicic tuff in this interval. The volcanic rocks probably correlate with other middle Tertiary volcanic rocks of the area, and the sedimentary rocks may correlate with the Cloudburst and Mineta Formations exposed on the flanks of the San Pedro Basin to the northeast. The Bisbee Group in the Exxon well is 602 m (1,975 f

Scientific Investigations Report↗
USGS2014

Origins and delineation of saltwater intrusion in the Biscayne aquifer and changes in the distribution of saltwater in Miami-Dade County, Florida

Intrusion of saltwater into parts of the shallow karst Biscayne aquifer is a major concern for the 2.5 million residents of Miami-Dade County that rely on this aquifer as their primary drinking water supply. Saltwater intrusion of this aquifer began when the Everglades were drained to provide dry land for urban development and agriculture. The reduction in water levels caused by this drainage, combined with periodic droughts, allowed saltwater to flow inland along the base of the aquifer and to seep directly into the aquifer from the canals. The approximate inland extent of saltwater was last mapped in 1995. An examination of the inland extent of saltwater and the sources of saltwater in the aquifer was completed during 2008–2011 by using (1) all available salinity information, (2) time-series electromagnetic induction log datasets from 35 wells, (3) time-domain electromagnetic soundings collected at 79 locations, (4) a helicopter electromagnetic survey done during 2001 that was processed, calibrated, and published during the study, (5) cores and geophysical logs collected from 8 sites for stratigraphic analysis, (6) 8 new water-quality monitoring wells, and (7) analyses of 69 geochemical samples. The results of the study indicate that as of 2011 approximately 1,200 square kilometers (km 2 ) of the mainland part of the Biscayne aquifer were intruded by saltwater. The saltwater front was mapped farther inland than it was in 1995 in eight areas totaling about 24.1 km 2 . In many of these areas, analyses indicated that saltwater had encroached along the base of the aquifer. The saltwater front was mapped closer to the coast than it was in 1995 in four areas totaling approximately 6.2 km 2 . The changes in the mapped extent of saltwater resulted from improved spatial information, actual movement of the saltwater front, or a combination of both. Salinity monitoring in some of the canals in Miami-Dade County between 1988 and 2010 indicated influxes of saltwater, with maximum salinities ranging from 1.4 to 32 practical salinity units (PSU) upstream of the salinity control structures. Time-series electromagnetic induction log data from monitoring wells G–3601, G–3608, and G–3701, located adjacent to the Biscayne, Snapper Creek, and Black Creek Canals, respectively, and upstream of the salinity control structures, indicated shallow influxes of conductive water in the aquifer that likely resulted from leakage of brackish water or saltwater from these canals. The determination that saltwater influxes were recent is supported by the similarity in the oxygen and hydrogen stable isotope composition in samples from the Snapper Creek Canal, 1.6 kilometers (km) inland of a salinity control structure, and in samples from well G–3608, which is adjacent to the canal, as well as by the relative ages of the water sampled from well G–3608 and other wells open to the aquifer below the saltwater interface. Historical and recent salinity information from the Card Sound Road Canal, monitoring well FKS8 located adjacent to the canal, and the 2001 helicopter electromagnetic survey indicated that saltwater may occasionally leak from this canal as far inland as 15 km. This leakage may be prevented or reduced by a salinity control structure that was installed in May 2010. Saltwater also may have leaked from the Princeton Canal. Results of geochemical sampling and analysis indicate a close correspondence between droughts and saltwater intrusion. Tritium/helium-3 apparent (piston-flow) ages determined from samples of saltwater with chloride concentrations of about 1,000 milligrams per liter (mg/L) or greater generally corresponded to a period during which droughts were frequent. Comparison of average daily air temperatures in Miami, Florida, with estimates of recharge temperatures determined from the dissolved gas composition in water samples indicated that saltwater likely entered the aquifer in April or early May when water levels are typically at their lowest during the year. Conversely, most of the samples of freshwater with chloride concentrations less than about 1,000 mg/L indicate recharge temperatures corresponding to air temperatures in mid to late May when rainfall and water levels in the aquifer increase, and the piston-flow ages of these samples correspond to wet years. The piston-flow ages of freshwater samples generally were younger than ages of samples of saltwater. Saltwater samples that were depleted in boron, magnesium, potassium, sodium, and sulfate, and enriched in calcium relative to the concentrations theoretically produced by freshwater/seawater mixing, generally were found to be associated with areas where saltwater had recently intruded. The calcium to (bicarbonate + sulfate) molar ratios (Ca/(HCO 3 +SO 4 )) of these samples generally were greater than 1. Saltwater samples from some of the monitoring wells, however, indicated little or no enrichment or depletion of these ions relative to the theoretical freshwater/seawater mixing line, and the Ca/(HCO 3 +SO 4 ) molar ratios of these samples generally were less than 1. Results indicated that aquifer materials are approaching equilibrium with seawater at these well locations.

Florida↗
USGS2007

Relation of Lake-Floor Characteristics to the Distribution of Variable Leaf Water-Milfoil in Moultonborough Bay, Lake Winnipesaukee, New Hampshire, 2005

Geophysical, water, and sediment surveys were done to characterize the effects of surficial geology, water and sediment chemistry, and surficial-sediment composition on the distribution of variable leaf water-milfoil in Moultonborough Bay, Lake Winnipesaukee, New Hampshire. Geophysical surveys were conducted in a 180-square-kilometer area, and water-quality and sediment samples were collected from 24 sites in the survey area during July 2005. Swath-bathymetric data revealed that Moultonborough Bay ranged in depth from less than 1 meter (m) to about 15 m and contained three embayments. Seismic-reflection profiles revealed erosion of the underlying bedrock and subsequent deposition of glaciolacustrine and Holocene lacustrine sediments within the survey area. Sediment thickness ranged from 5 m along the shoreward margins to more than 15 m in the embayments. Data from sidescan sonar, surficial-sediment samples, bottom photographs, and video revealed three distinct lake-floor environments: rocky nearshore, mixed nearshore, and muddy basin. Rocky nearshore environments were found in shallow water (less than 5 m deep) and contained sediments ranging from coarse silt to very coarse sand. Mixed nearshore environments also were found in shallow water and contained sediments ranging from silt to coarse sand with different densities of aquatic vegetation. Muddy basin environments contained the finest-grained sediments, ranging from fine to medium silt, and were in the deepest waters of the bay. Acoustic Ground Discrimination Systems (AGDS) survey data revealed that 86 percent of the littoral zone (the area along the margins of the bay and islands that extends from 0 to 4.3 m in water depth) contained submerged aquatic vegetation (SAV) in varying densities: approximately 36 percent contained SAV bottom cover of 25 percent or less, 43 percent contained SAV bottom cover of more than 25 and less than 75 percent, and approximately 7 percent contained SAV bottom cover of more than 75 percent. SAV included variable leaf water-milfoil, native milfoil, bassweed, pipewort, and other species, which were predominantly found near shoreward margins and at depths ranging from less than 1 to 4 m. AGDS data were used in a Geographic Information System to generate an interpolated map that distinguished variable leaf water-milfoil from other SAV. Furthermore, these data were used to isolate areas susceptible to variable leaf water-milfoil growth. Approximately 21 percent of the littoral zone contained dense beds (more than 59 percent bottom cover) of variable leaf water-milfoil, and an additional 44 percent was determined to be susceptible to variable leaf water-milfoil infestation. Depths differed significantly between sites with variable leaf water-milfoil and sites with other SAV (p = 0.04). Variable leaf water-milfoil was found at depths that ranged from 1 to 4 m, and other SAV had a depth range of 1 to 2 m. Although variable leaf water-milfoil was observed at greater depths than other SAV, it was not observed below the photic zone. Analysis of constituent concentrations from the water column, interstitial pore water, and sediment showed little correlation with the presence of variable leaf water-milfoil, with two exceptions. Iron concentrations were significantly lower at variable leaf water-milfoil sites than at other sampling sites (p = 0.04). Similarly, the percentage of total organic carbon also was significantly lower at the variable leaf water-milfoil sites than at other sampling sites (p = 0.04). Surficial-sediment-grain size had the greatest correlation to the presence of variable leaf water-milfoil. Variable leaf water-milfoil was predominantly growing in areas of coarse sand (median grain-size 0.62 millimeters). Surficial-sediment-grain size was also correlated with total ammonia plus organic nitrogen (Rho = 0.47; p = 0.02) and with total phosphorus (Rho = 0.44; p = 0.05) concentrations in interstitial pore-water samples.

Scientific Investigations Report↗
USGS2001

Mars Global Surveyor Thermal Emission Spectrometer experiment: Investigation description and surface science results

The Thermal Emission Spectrometer (TES) investigation on Mars Global Surveyor (MGS) is aimed at determining (1) the composition of surface minerals, rocks, and ices; (2) the temperature and dynamics of the atmosphere; (3) the properties of the atmospheric aerosols and clouds; (4) the nature of the polar regions; and (5) the thermophysical properties of the surface materials. These objectives are met using an infrared (5.8- to 50-μm) interferometric spectrometer, along with broadband thermal (5.1- to 150-μm) and visible/near-IR (0.3- to 2.9-μm) radiometers. The MGS TES instrument weighs 14.47 kg, consumes 10.6 W when operating, and is 23.6×35.5×40.0 cm in size. The TES data are calibrated to a 1-σ precision of 2.5 −6 ×10 −8 W cm −2 sr −1 /cm −1 , 1.6×10 −6 W cm −2 sr −1 , and ∼0.5 K in the spectrometer, visible/near-IR bolometer, and IR bolometer, respectively. These instrument subsections are calibrated to an absolute accuracy of ∼4×10 −8 W cm −2 sr −1 /cm −1 (0.5 K at 280 K), 1–2%, and ∼1–2 K, respectively. Global mapping of surface mineralogy at a spatial resolution of 3 km has shown the following: (1) The mineralogic composition of dark regions varies from basaltic, primarily plagioclase feldspar and clinopyroxene, in the ancient, southern highlands to andesitic, dominated by plagioclase feldspar and volcanic glass, in the younger northern plains. (2) Aqueous mineralization has produced gray, crystalline hematite in limited regions under ambient or hydrothermal conditions; these deposits are interpreted to be in-place sedimentary rock formations and indicate that liquid water was stable near the surface for a long period of time. (3) There is no evidence for large-scale (tens of kilometers) occurrences of moderate-grained (>50-μm) carbonates exposed at the surface at a detection limit of ∼10%. (4) Unweathered volcanic minerals dominate the spectral properties of dark regions, and weathering products, such as clays, have not been observed anywhere above a detection limit of ∼10%; this lack of evidence for chemical weathering indicates a geologic history dominated by a cold, dry climate in which mechanical, rather than chemical, weathering was the significant form of erosion and sediment production. (5) There is no conclusive evidence for sulfate minerals at a detection limit of ∼15%. The polar region has been studied with the following major conclusions: (1) Condensed CO 2 has three distinct end-members, from fine-grained crystals to slab ice. (2) The growth and retreat of the polar caps observed by MGS is virtually the same as observed by Viking 12 Martian years ago. (3) Unique regions have been identified that appear to differ primarily in the grain size of CO 2 ; one south polar region appears to remain as black slab CO 2 ice throughout its sublimation. (4) Regional atmospheric dust is common in localized and regional dust storms around the margin and interior of the southern cap. Analysis of the thermophysical properties of the surface shows that (1) the spatial pattern of albedo has changed since Viking observations, (2) a unique cluster of surface materials with intermediate inertia and albedo occurs that is distinct from the previously identified low-inertia/bright and high-inertia/dark surfaces, and (3) localized patches of high-inertia material have been found in topographic lows and may have been formed by a unique set of aeolian, fluvial, or erosional processes or may be exposed bedrock.

Journal of Geophysical Research E: Planets↗
USGS2019

Potential interaction of groundwater and surface water including autonomous underwater vehicle reconnaissance at Nolin River Lake, Kentucky, 2016

The U.S. Geological Survey collaborated with the U.S. Army Corps of Engineers, Louisville District, on a synoptic study of water quality at Nolin River Lake during August 2016. The purpose of the study was to develop a better understanding of the potential for interaction between groundwater and surface water at Nolin River Lake, Kentucky. Groundwater can have properties that are measurably different from those in adjacent surface water, and inflows and outflows can be an important component of water quality and quantity. An improved understanding of potential interaction of groundwater and surface water at Nolin River Lake may be used to refine lake-management strategies. This study (1) compiled and interpreted existing information to characterize the hydrogeological setting and implications for potential interaction of groundwater and surface water in the Nolin River Lake watershed; (2) collected transects of onsite water-quality parameters using an autonomous underwater vehicle (AUV) in areas with potential for interaction of groundwater and surface water, including five sites on Nolin River Lake and one site on the Nolin River; and (3) collected discrete water-quality and phytoplankton community data at the same six sites. A review of existing hydrogeologic information did not indicate the presence of karst features adjacent to or beneath Nolin River Lake that would facilitate groundwater interaction with the reservoir. Observations leading to this conclusion include (1) limestone that is adjacent to the shoreline and perhaps beneath the lake, is overlain with siliciclastic rocks and fine-grained sediment that inhibits infiltration and development of karst features that encourage rapid groundwater flow; (2) the geologic deposits surrounding the reservoir are described as having limited or no potential for development of karst features, some exceptions may exist in tributary valleys; (3) very few karst features were mapped within 1 mile of the reservoir or in the area currently occupied by the reservoir; and (4) faults that intersect the reservoir but may not possess hydraulic properties that cause the faults to be conduits for groundwater flow. Groundwater interaction with reservoir tributaries is likely more common in areas of the watershed upstream from Nolin River Lake where karst hydrogeology is prevalent. Results of water-quality surveys using an AUV from August 15 to 19, 2016, did not identify areas of anomalous values that might indicate groundwater inflows through preferential flow zones. Spatial distributions of water-quality parameters were generally uniform within each constant-depth layer. The constant-depth layers were selected to be above, within, and below the thermocline and ranged from the water surface to 25 feet. Surveys near the bottom of the reservoir that might have been more sensitive to groundwater inflows were not done because presurvey data were not available to indicate locations of obstacles that could ensnare the AUV. Water-quality data collected with the AUV did identify water-quality anomalies where stream tributaries were discharging to the reservoir. The discrete water-quality samples indicated uniformity among the five reservoir sites. The riverine site that is immediately upstream from Nolin River Lake, however, had some unique water-quality characteristics relative to sites on the reservoir. The highest concentrations of nitrate plus nitrite as nitrogen (0.145 milligrams per liter [mg/L]), total phosphorous (0.07 mg/L), chlorophyll a (36.1 micrograms per liter), and pheophytin a (10.2 micrograms per liter) were measured at the Nolin River Lake riverine site (site 2NRR20034). The concentrations of nutrients and chlorophyll a at the riverine site did exceed the 25th percentile of median concentrations measured by the U.S. Environmental Protection Agency (EPA) at other lakes and reservoirs in EPA level IV ecoregion 71a. Concentrations of most nutrients and chlorophyll a at the five reservoir sites also exceeded the 25th percentile of median concentrations in EPA level IV ecoregion 72h. The exception was the concentrations of total phosphorus as phosphorus at the reservoir sites that were at or below the 25th percentile of median concentrations measured by EPA (0.03 mg/L). Concentrations of orthophosphate as phosphorus were less than the method detection limit of 0.004 mg/L at all sites. The phytoplankton community in Nolin River Lake was almost exclusively (greater than 90 percent of total phytoplankton abundance) cyanobacteria, also known as blue-green algae. A species of Cylindrospermopsis dominated the cyanobacterial community at the five reservoir sites, while Chroococcus microscopicus was most abundant at the riverine site. Cyanobacterial cell densities ranged from 10,000 to 198,067,460 cells per liter in five areas in the reservoir and from 4,800 to 73,751,253 cells per liter at the riverine site. Multiple potential sources of water to Nolin River Lake include direct precipitation, overland flow, interflow, groundwater, and surface water. Understanding the exact contribution of each of these components to the water budget at Nolin River Lake may help the U.S. Army Corps of Engineers manage the water quality, water quantity, and biological communities in the reservoir. Additional hydrogeologic and water-quality data that builds on the results of this study may refine the inferences of this study; for example, deeper AUV surveys that target the largest fault zones might further the understanding of the potential for groundwater flow through those features. A complete understanding of the reservoir hydrology, however, may require the use of scientific methods intended for water bodies as large as Nolin River Lake, such as aerial infrared photography and imagery; water mass, chemical, and isotopic balance studies; geophysical measurements; and numerical simulations.

Kentucky↗
USGS2013

Geohydrology, geochemistry, and groundwater simulation (1992-2011) and analysis of potential water-supply management options, 2010-60, of the Langford Basin, California

Groundwater withdrawals began in 1992 from the Langford Basin within the Fort Irwin National Training Center (NTC), California. From April 1992 to December 2010, approximately 12,300 acre-feet of water (averaging about 650 acre-feet per year) has been withdrawn from the basin and transported to the adjacent Irwin Basin. Since withdrawals began, water levels in the basin have declined by as much as 40 feet, and the quality of the groundwater withdrawn from the basin has deteriorated. The U.S. Geological Survey collected geohydrologic data from Langford Basin during 1992–2011 to determine the quantity and quality of groundwater available in the basin. Geophysical surveys, including gravity, seismic refraction, and time-domain electromagnetic induction surveys, were conducted to determine the depth and shape of the basin, to delineate depths to the Quaternary-Tertiary interface, and to map the depth to the water table and changes in water quality. Data were collected from existing wells and test holes, as well as 11 monitor wells that were installed at 5 sites as part of this study. Water-quality samples collected from wells in the basin were used to determine the groundwater chemistry within the basin and to delineate potential sources of poor-quality groundwater. Analysis of stable isotopes of oxygen and hydrogen in groundwater indicates that present-day precipitation is not a major source of recharge to the basin. Tritium and carbon-14 data indicate that most of the basin was recharged prior to 1952, and the groundwater in the basin has an apparent age of 12,500 to 30,000 years. Recharge to the basin, estimated to be less than 50 acre-feet per year, has not been sufficient to replenish the water that is being withdrawn from the basin. A numerical groundwater-flow model was developed for the Langford Basin to better understand the aquifer system used by the Fort Irwin NTC as part of its water supply, and to provide a tool to help manage groundwater resources at the NTC. Measured groundwater-level declines since the initiation of withdrawals (1992–2011) were used to calibrate the groundwater-flow model. The simulated recharge was about 46 acre-feet per year, including approximately 6 acre-feet per year of natural recharge derived from precipitation runoff and as much as 40 acre-feet per year of underflow from the Irwin Basin. Between April 1992 and December 2010, an average of about 650 acre-feet per year of water was withdrawn from the Langford Basin. Groundwater withdrawals in excess of natural recharge resulted in a net loss of 11,670 acre-feet of groundwater storage within the basin for the simulation period. The Fort Irwin NTC is considering various groundwater-management options to address the limited water resources in the Langford Basin. The calibrated Langford Basin groundwater-flow model was used to evaluate the hydrologic effects of four groundwater-withdrawal scenarios being considered by the Fort Irwin NTC over the next 50 years (January 2011 through December 2060). Continuation of the 2010 withdrawal rate in the three existing production wells will result in 70 feet of additional drawdown in the central part of the basin. Redistributing the 2010 withdrawal rate equally to the three existing wells and two proposed new wells in the northern and southern parts of the basin would result in about 10 feet less drawdown in the central part of the basin but about 100 feet of additional drawdown in the new well in the northern part of the basin and about 50 feet of additional drawdown in the new well in the southern part of the basin. Reducing the withdrawals from the three existing production wells in the central part of the basin from about 45,000 acre-feet to about 32,720 acre-feet would result in about 40 feet of additional drawdown in the central basin near the pumping wells, about 25 feet less than if withdrawals were not reduced. The combination of reducing and redistributing the cumulative withdrawals to the three existing and two proposed new wells results in about 40 feet of additional drawdown in the central and southern parts of the basin and about 70 feet in the northern part of the basin. These results show that reducing and redistributing the groundwater withdrawals would maintain the upper aquifer at greater than 50 percent of its predevelopment saturated thickness throughout the groundwater basin. The scenarios simulated for this study demonstrate how the calibrated model can be utilized to evaluate the hydrologic effects of different water-management strategies.

California↗
USGS2015

Simulation of groundwater flow and analysis of the effects of water-management options in the North Platte Natural Resources District, Nebraska

The North Platte Natural Resources District (NPNRD) has been actively collecting data and studying groundwater resources because of concerns about the future availability of the highly inter-connected surface-water and groundwater resources. This report, prepared by the U.S. Geological Survey in cooperation with the North Platte Natural Resources District, describes a groundwater-flow model of the North Platte River valley from Bridgeport, Nebraska, extending west to 6 miles into Wyoming. The model was built to improve the understanding of the interaction of surface-water and groundwater resources, and as an optimization tool, the model is able to analyze the effects of water-management options on the simulated stream base flow of the North Platte River. The groundwater system and related sources and sinks of water were simulated using a newton formulation of the U.S. Geological Survey modular three-dimensional groundwater model, referred to as MODFLOW–NWT, which provided an improved ability to solve nonlinear unconfined aquifer simulations with wetting and drying of cells. Using previously published aquifer-base-altitude contours in conjunction with newer test-hole and geophysical data, a new base-of-aquifer altitude map was generated because of the strong effect of the aquifer-base topography on groundwater-flow direction and magnitude. The largest inflow to groundwater is recharge originating from water leaking from canals, which is much larger than recharge originating from infiltration of precipitation. The largest component of groundwater discharge from the study area is to the North Platte River and its tributaries, with smaller amounts of discharge to evapotranspiration and groundwater withdrawals for irrigation. Recharge from infiltration of precipitation was estimated with a daily soil-water-balance model. Annual recharge from canal seepage was estimated using available records from the Bureau of Reclamation and then modified with canal-seepage potentials estimated using geophysical data. Groundwater withdrawals were estimated using land-cover data, precipitation data, and published crop water-use data. For fields irrigated with surface water and groundwater, surface-water deliveries were subtracted from the estimated net irrigation requirement, and groundwater withdrawal was assumed to be equal to any demand unmet by surface water. The groundwater-flow model was calibrated to measured groundwater levels and stream base flows estimated using the base-flow index method. The model was calibrated through automated adjustments using statistical techniques through parameter estimation using the parameter estimation suite of software (PEST). PEST was used to adjust 273 parameters, grouped as hydraulic conductivity of the aquifer, spatial multipliers to recharge, temporal multipliers to recharge, and two specific recharge parameters. Base flow of the North Platte River at Bridgeport, Nebraska, streamgage near the eastern, downstream end of the model was one of the primary calibration targets. Simulated base flow reasonably matched estimated base flow for this streamgage during 1950–2008, with an average difference of 15 percent. Overall, 1950–2008 simulated base flow followed the trend of the estimated base flow reasonably well, in cases with generally increasing or decreasing base flow from the start of the simulation to the end. Simulated base flow also matched estimated base flow reasonably well for most of the North Platte River tributaries with estimated base flow. Average simulated groundwater budgets during 1989–2008 were nearly three times larger for irrigation seasons than for non-irrigation seasons. The calibrated groundwater-flow model was used with the Groundwater-Management Process for the 2005 version of the U.S. Geological Survey modular three-dimensional groundwater model, MODFLOW–2005, to provide a tool for the NPNRD to better understand how water-management decisions could affect stream base flows of the North Platte River at Bridgeport, Nebr., streamgage in a future period from 2008 to 2019 under varying climatic conditions. The simulation-optimization model was constructed to analyze the maximum increase in simulated stream base flow that could be obtained with the minimum amount of reductions in groundwater withdrawals for irrigation. A second analysis extended the first to analyze the simulated base-flow benefit of groundwater withdrawals along with application of intentional recharge, that is, water from canals being released into rangeland areas with sandy soils. With optimized groundwater withdrawals and intentional recharge, the maximum simulated stream base flow was 15–23 cubic feet per second (ft 3 /s) greater than with no management at all, or 10–15 ft 3 /s larger than with managed groundwater withdrawals only. These results indicate not only the amount that simulated stream base flow can be increased by these management options, but also the locations where the management options provide the most or least benefit to the simulated stream base flow. For the analyses in this report, simulated base flow was best optimized by reductions in groundwater withdrawals north of the North Platte River and in the western half of the area. Intentional recharge sites selected by the optimization had a complex distribution but were more likely to be closer to the North Platte River or its tributaries. Future users of the simulation-optimization model will be able to modify the input files as to type, location, and timing of constraints, decision variables of groundwater withdrawals by zone, and other variables to explore other feasible management scenarios that may yield different increases in simulated future base flow of the North Platte River.

Nebraska↗
USGS2016

Hydrologic assessment of the shallow groundwater flow system beneath the Shinnecock Nation tribal lands, Suffolk County, New York

Defining the distribution and flow of shallow groundwater beneath the Shinnecock Nation tribal lands in Suffolk County, New York, is a crucial first step in identifying sources of potential contamination to the surficial aquifer and coastal ecosystems. The surficial or water table aquifer beneath the tribal lands is the primary source of potable water supply for at least 6 percent of the households on the tribal lands. Oyster fisheries and other marine ecosystems are critical to the livelihood of many residents living on the tribal lands, but are susceptible to contamination from groundwater entering the embayment from the surficial aquifer. Contamination of the surficial aquifer from flooding during intense coastal storms, nutrient loading from fertilizers, and septic effluent have been identified as potential sources of human and ecological health concerns on tribal lands. The U.S. Geological Survey (USGS) facilitated the installation of 17 water table wells on and adjacent to the tribal lands during March 2014. These wells were combined with other existing wells to create a 32-well water table monitoring network that was used to assess local hydrologic conditions. Survey-grade, global-navigation-satellite systems provided centimeter-level accuracy for positioning wellhead surveys. Water levels were measured by the USGS during May (spring) and November (fall) 2014 to evaluate seasonal effects on the water table. Water level measurements were made at high and low tide during May 2014 to identify potential effects on the water table caused by changes in tidal stage (tidal flux) in Shinnecock Bay. Water level contour maps indicate that the surficial aquifer is recharged by precipitation and upgradient groundwater flow that moves from the recharge zone located generally beneath Sunrise Highway, to the discharge zone beneath the tribal lands, and eventually discharges into the embayment, tidal creeks, and estuaries that bound the tribal lands to the east, south, and west. Water levels in many of the wells in the network fluctuated in response to precipitation, upgradient groundwater flow, and tidal flux in Shinnecock Bay. Water level altitudes ranged from 6.66 to 0.47 feet (ft) above the North American Vertical Datum of 1988 during the spring measurement period, and from 5.25 to -0.24 ft (NAVD 88) during fall 2014. Historically, annual and seasonal precipitation seem to indicate long-term water level trends in an index well located in the town of Southampton, correlates with changes in storage in the upper glacial aquifer, but does not necessarily indicate water level extremes in the shallow groundwater system. To place the study period in perspective, calendar year 2014 was the 32d wettest year on record, with precipitation for the year totaling 48.1 inches, a 2.6-percent increase from the annual average (46.9 inches per year), based on 81 years of complete record at the National Oceanographic and Atmospheric Administration, National Weather Service cooperative meteorological station at Bridgehampton, New York. Estimated recharge to the water table beneath the tribal lands from precipitation for 2014 is 25.4 inches. Tidal flux caused water levels in wells to fluctuate from 0.30 to -0.24 ft during May 2014. Water levels in wells located north of Old Fort Pond and beneath the southernmost extent of the tribal lands were most influenced by tidal flux. During June 2014, hydrographs indicate that tidal flux influenced water levels by 0.48 ft in a well located near the southernmost extent of the tribal lands approximately 0.3 miles north of Shinnecock Bay, and was zero at a well located approximately 0.5 miles south of Montauk Highway, and 0.4 miles west of Heady Creek, near the geographic center of the tribal lands. Tidal-influence delay time (time interval between peak high-tide stage and corresponding peak high-water level) ranged from 1.75 hours at the well located near the southernmost extent of the tribal lands, to more than 4 hours at a well located north of Old Fort Pond, near the northwestern part of the tribal lands. Estimated hydraulic-conductivity values derived from the results of specific-capacity tests that were completed at nine observation wells during March 2015 were used to calculate average linear velocity. Average linear velocity along conceptualized flow-path segments of the upper glacial aquifer located beneath the tribal lands was estimated using an assumed effective porosity value, and hydraulic-conductivity and hydraulic-head values that were interpolated from measured values. Groundwater travel times were estimated by dividing the length of the flow-path segment by the average linear velocity along the flow-path segment. Total estimated groundwater travel time along a conceptualized flow path, beginning near Sunrise Highway and terminating at Shinnecock Bay, is approximately 45 years using a porosity value of 30 percent. A surficial-silty unit was identified from approximately 0 to 10 ft below land surface at multiple locations beneath the tribal lands. The lithology of the surficial unit was verified by interpreted gamma log results obtained from select wells, and auger-rig drill cuttings from an observation well located near the geographic center of the tribal lands. The altitude of the unit varies with topography and was delineated along a cross section line that trends north-south along the approximate centerline (spine) of the tribal lands. The altitude of the hydrogeologic contact between the upper glacial and the Magothy aquifers generally decreases from northwest to southeast, occurs at a depth ranging from about 150 to 200 ft beneath the tribal lands, and was identified at two locations north of the tribal lands, near Sunrise Highway and Sebonac Road. Results of electrical geophysical surveys indicate that the depth to the freshwater/saltwater interface decreases from north to south with decreasing water level altitude, and the Magothy and upper glacial aquifers contain saltwater at varying depths along the north-south trending section. Results of the surveys also indicate that the Magothy aquifer beneath the tribal lands contains brackish and salty water and is not considered a source of potable water supply. In general, depth to the interface increases with increasing geographic distance from the coastline. Low water table altitudes can result in increased saltwater encroachment into the surficial aquifer beneath the tribal lands. This upward movement and shallow depth of the freshwater/saltwater interface can jeopardize water quality in wells that supply water for domestic use.

New York↗
USGS2006

Hydrogeologic framework refinement, ground-water flow and storage, water-chemistry analyses, and water-budget components of the Yuma area, southwestern Arizona and southeastern California

The ground-water and surface-water system in the Yuma area in southwestern Arizona and southeastern California is managed intensely to meet water-delivery requirements of customers in the United States, to manage high ground-water levels in the valleys, and to maintain treaty-mandated water-quality and quantity requirements of Mexico. The following components in this report, which were identified to be useful in the development of a ground-water management model, are: (1) refinement of the hydrogeologic framework; (2) updated water-level maps, general ground-water flow patterns, and an estimate of the amount of ground water stored in the mound under Yuma Mesa; (3) review and documentation of the ground-water budget calculated by the Bureau of Reclamation, U.S. Department of the Interior (Reclamation); and (4) water-chemistry characterization to identify the spatial distribution of water quality, information on sources and ages of ground water, and information about the productive-interval depths of the aquifer. A refined three-dimensional digital hydrogeologic framework model includes the following hydrogeologic units from bottom to top: (1) the effective hydrologic basement of the basin aquifer, which includes the Pliocene Bouse Formation, Tertiary volcanic and sedimentary rocks, and pre-Tertiary metamorphic and plutonic rocks; (2) undifferentiated lower units to represent the Pliocene transition zone and wedge zone; (3) coarse-gravel unit; (4) lower, middle, and upper basin fill to represent the upper, fine-grained zone between the top of the coarse-gravel unit and the land surface; and (5) clay A and clay B. Data for the refined model includes digital elevation models, borehole lithology data, geophysical data, and structural data to represent the geometry of the hydrogeologic units. The top surface of the coarse-gravel unit, defined by using borehole and geophysical data, varies similarly to terraces resulting from the down cutting of the Colorado River. Clay A is nearly the same as the previous conceptual hydrogeologic model definition (Olmsted and others, 1973), except for a minor westward extension from the city of Yuma. Clay B is extended to the southerly international boundary and increased in areal extent by about two-thirds of the original extent (Olmsted and others, 1973). The other hydrogeologic units generally are the same as in the previous conceptual hydrogeologic model. Before development, the Colorado and Gila Rivers were the sources of nearly all the ground water in the Yuma area through direct infiltration of water from river channels and annual overbank flooding. After construction of upstream reservoirs and clearing and irrigation of the floodplains, the rivers now act as drains for the ground water. Ground-water levels in most of the Yuma area are higher now than they were in predevelopment time. A general gradient of ground-water flow toward the natural discharge area south of the Yuma area still exists, but many other changes in flow are evident. Ground water in Yuma Valley once flowed away from the Colorado River, but now has a component of flow towards the river and Mexicali Valley. A ground-water mound has formed under Yuma Mesa from long-term surface-water irrigation; about 600,000 to 800,000 acre-ft of water are stored in the mound. Ground-water withdrawals adjacent to the southerly international boundary have resulted in water-level declines in that area. The reviewed and documented water budget includes the following components: (1) recharge in irrigated areas, (2) evapotranspiration by irrigated crops and phreatophytes, (3) ground-water return flow to the Colorado River, and (4) ground-water withdrawals (including those in Mexicali Valley). Recharge components were calculated by subtracting the amount of water used by crops from the amount of water delivered. Evapotranspiration rates were calculated on the basis of established methods, thus were appropriate for input to the ground-water flow model developed by the Bureau of Reclamation (William Greer, hydrologist, Bureau of Reclamation, written commun., 2005). Evapotranspiration by crops and phreatophytes were calculated by using crop coefficient methods and meteorological data. Other methods of calculating evapotranspiration rates by using combinations of satellite imagery and ground-based data could be used for higher spatial and temporal resolution. Ground-water return flow during years of low flow on the Colorado River (1972–82, 1987–92, and 1994–96) averaged 79,000 acre-ft per year. Ground-water withdrawal data for 1970–99 were similar to other estimates made by the U.S. Geological Survey for the Yuma area. New water-chemistry data were collected in 12 wells and 8 canals/drains to characterize spatial patterns in chemical constituents, determine isotopic ages of water, infer possible sources of ground water, and locate the vertical intervals of the aquifer that contribute most water to wells. Depth-dependent samples were collected at one of the wells (YM-10). A large quantity of water-quality data were compiled from Bureau of Reclamation and U.S. Geological Survey records and merged into the U.S. Geological Survey National Water Information System database. New samples were analyzed for major ions, nutrients, stable isotopes of oxygen and hydrogen, tritium ( 3 H), and carbon-14 ( 14 C) (along with C 13 /C 12 ratios). Light values of oxygen-18 ( 18 O) and deuterium ( 2 H, D) in well 242-2 indicate recharge from the Colorado River. Heavy water samples from wells 242-22, CADC, and Mesa del Sol indicate local recharge sources. Tritium data indicate there is young water in wells in the valleys and near the edge of Yuma Mesa, while older water is found far from the Colorado River. 14 C data indicate that water from wells near the southerly international boundary is at least several thousand years old.

Arizona, California↗
USGS2018

Santa Barbara and Foothill groundwater basins Geohydrology and optimal water resources management—Developed using density dependent solute transport and optimization models

Groundwater has been a part of the city of Santa Barbara’s water-supply portfolio since the 1800s; however, since the 1960s, the majority of the city’s water has come from local surface water, and the remainder has come from groundwater, State Water Project, recycled water, increased water conservation, and as needed, seawater desalination. Although groundwater from the Santa Barbara and Foothill groundwater basins only accounts for a small percentage of the long-term supply, it is an important source of supplemental water during times of surface-water shortages. During the late 1980s and early 1990s, production wells extracted additional groundwater to compensate for drought related water-delivery shortfalls from other sources; in response, water levels declined substantially in the Santa Barbara and Foothill groundwater basins (below sea level in the Santa Barbara groundwater basin). In coastal basins that have groundwater extraction near shore, seawater intrusion is often a problem. Seawater intrusion in the Santa Barbara groundwater basin is thought to be more limited than in other coastal basins because of an offshore fault that acts as a partial barrier to groundwater flow. During the late 1980s and early 1990s, seawater intrusion was observed in the Santa Barbara groundwater basin, as indicated by increased chloride concentrations at several monitoring wells that ranged from 200 ft to 1,300 ft from the ocean and as close as 2,900 ft to the nearest pumping well. This demonstrated that seawater can intrude into the Santa Barbara groundwater basin when groundwater levels fall below sea level near the coast. The city of Santa Barbara is interested in developing a better understanding of the sustainability of its groundwater supplies. In 2014, California adopted historic legislation to manage its groundwater: the Sustainable Groundwater Management Act (SGMA). The SGMA requires the development and implementation of “Groundwater Sustainability Plans” in 127 priority groundwater basins; although Santa Barbara was not a designated priority basin, the city is taking steps to achieve sustainability. Sustainability was defined in the SGMA in terms of avoiding undesirable results: significant and unreasonable groundwater-level declines, reduction in groundwater storage, seawater intrusion, water-quality degradation, land subsidence, and surface-water depletion. In this project, a cooperative study between the U.S. Geological Survey (USGS) and the city of Santa Barbara, sustainable yield is defined as the volume of groundwater that can be pumped from storage without causing water-level drawdowns and the associated increases in seawater intrusion (as indicated by increases in measured chloride concentrations) at selected wells. In order to estimate the sustainability of Santa Barbara’s groundwater basins, a three-dimensional density-dependent groundwater-flow and solute-transport model (the Santa Barbara Flow and Transport Model, or SBFTM) was developed on the basis of an existing groundwater-flow model. To simulate seawater intrusion to the Santa Barbara Basin under various management strategies, the SBFTM uses the USGS code SEAWAT to simulate salinity transport and variable-density flow. The completed SBFTM was coupled with a management optimization tool, in this case a multi-objective evolutionary algorithm, to determine optimal pumping strategies that maximize the sustainable yield and at the same time satisfy user-defined drawdown and chloride-concentration constraints. As part of this study, a three-dimensional hydrogeologic framework model was developed to quantify the extent and hydrogeologic characteristics of the Santa Barbara and Foothill groundwater basins and to help define the discretization and hydraulic properties used in the SBFTM. The development of the hydrogeologic framework model required the collection and reconciliation of geologic and geophysical data from existing maps, reports, and databases, along with geologic and hydrologic data from recently drilled wells. These data were integrated into a three-dimensional hydrogeologic framework model that defines the stratigraphy and geometry of the aquifer zones and the major geologic structures in the basin. The hydrogeologic framework model also quantifies the variation in sediment grain size within each aquifer zone as the percentage of coarse-grained sediment. Previous studies indicated that there are two principal water-producing zones in the Santa Barbara groundwater basin, the upper and lower producing zones; an additional thin, productive zone was identified as part of this study. This “middle producing zone” is not as areally extensive as the upper and lower producing zones and only exists in the coastal part of Storage Unit I. These producing zones are bounded at depth by less productive shallow, middle, and deep zones. Two versions of the SBFTM were constructed: an initial-condition model and a modern transient model. The initial-condition model is a long-term transient model that simulates flow and solute-transport conditions during a period with limited anthropogenic influences preceeding the modern transient model. The simulation-transient model simulates flow and transport conditions from 1929 through 2013; however, because of data availability, the focus of the model calibration was 1972–2013. The SBFTM was calibrated to measured groundwater levels and drawdown, as well as measured chloride concentrations and change in concentrations, using a combination of automated and trial-and-error parameter-estimation techniques. A sensitivity analysis indicated that, in general, the SBFTM was most sensitive to recharge- and pumping-distribution parameters, specifically those controlling the amount of small-catchment recharge and the distribution of water extraction by hydrogeologic layer for production wells. The model was also sensitive to parameters controlling stream-recharge rates, horizontal and vertical hydraulic conductivity, and porosity. From 1929 to 1971, most of the water entering the area represented by the SBFTM was from creek and small-catchment recharge, and the majority of water leaving the SBFTM area was from pumping, discharge to creeks, and drains. In addition, about 37 percent of the total pumpage came from a net reduction in groundwater storage. From 1972 to 2013, the amount of water entering and leaving the SBFTM was fairly similar as that from 1929 to 1971, except the reduction in pumpage added about 17,000 acre-ft of water to storage. During this later period, there were also times of storage loss. For example, during July 1990, a month when approximately 705 acre-ft of groundwater was pumped in the study area, the pumpage was much greater than all sources of recharge combined, and about 382 acre-ft of water was removed from groundwater storage. Simulated hydraulic heads replicated the observed data to an acceptable matching of the measured water-level, flow direction, and vertical gradients. Simulated hydrographs for selected wells were in good agreement with the measured data, with an average residual of -2.7 ft and a standard deviation of 14.5 ft, indicating that the simulated heads, on average, underestimated the observed water levels. An examination of the model fit indicated that most of the discrepancies were lower simulated heads at wells proximal to production well sites. The simulated chloride concentrations reasonably matched the rising limbs of the measured breakthrough curves in terms of timing and magnitude; however, the simulation overestimated the chloride concentrations on the falling limbs. The overestimation of low chloride concentrations was attributed to the model overestimating the advance of the chloride front during periods of heavy pumping and underestimating the retreat of the chloride front during periods of low pumping. These simulation errors would result in a conservative response by local water managers to seawater intrusion. The SBFTM was used to develop a collection of predictive simulations optimized to produce pumping schedules that maximize yield, subject to a set of constraints and competing objectives. The simulations were grouped as scenarios that differed in their time horizon, initial conditions for groundwater levels and chloride concentrations, as well as precipitation, which was incorporated into the model through simulated recharge. Overall, five scenarios were developed in a multi-objective framework to obtain optimal pumping rates for all of the wells managed by the city, while minimizing excessive drawdown and seawater intrusion. For the current study, complexities in the simulation model and the optimization formulation required additional considerations. Incorporating the solute-transport equations to simulate chloride transport added a highly nonlinear process that is solved iteratively in each time step of the groundwater-flow model. These nonlinearities, coupled with the highly refined grid in the current model, creates challenges for many traditional optimization methods. Therefore, an optimization method was needed that could address nonlinear relationships as well as a very large problem size. Lastly, the optimization problem was reformulated to include multiple objectives without requiring convergence to a single solution. This approach, guided by the city’s objectives, allowed the maximum extraction of information from the complex simulation. Borg, a multi-objective evolutionary algorithm, was chosen as the optimization algorithm for this study for several reasons: (1) it is very computationally efficient; (2) it can run in parallel; (3) it requires little user input; and (4) it can solve for multiple competing objectives. The first three points allow the algorithm to proceed toward the optimal solutions at the fastest possible rate. The fourth point is advantageous for large, complex optimization problems because it is difficult to formulate the optimization problem in a way that produces only one optimal solution. The problem formulation consisted of four competing objectives and a constraint set in accordance with the main concerns of the city. The objectives were maximizing total pumpage, minimizing seawater intrusion, minimizing total drawdown in production wells, and minimizing the maximum drawdown. The constraints were pump capacity, meeting drinking-water standards for chloride, maintaining a specified minimum flowrate to a groundwater treatment plant, and maintaining minimum water levels in pumping wells. The decision variables either were quarterly pumpage by well or total pumpage by basin. Five optimization scenarios were developed that allow the decision makers to evaluate a range of optimal solutions for a variety of water levels and chloride concentrations as well as potential future climatic conditions. Three scenarios (1, 2, and 5) were multi-objective optimization formulations that allowed for variations in management preferences and climatic conditions. The other two scenarios (3 and 4) were designed to examine the optimization results to answer specific questions. Scenario 1 described the best-case sustainable yield assuming a “full” basin (that is, high initial water levels) and typical climate conditions for 10 years. Scenario 2 also started with a “full” basin; however, this was followed by a 10-year drought. Scenario 3 determined if an “empty” basin (that is, low initial water levels) would recover to full conditions (1998 conditions) given climate assumptions and optimal pumping schedules from scenarios 1 and 2. Scenario 4 was designed to produce decision rules that can be used by water managers to help choose an optimal pumping schedule based on measured water-level or chloride data. Scenario 5 identified future pumping schedules based on short-term climate variations during a 2-year management horizon. The results from scenarios 1 and 2 described the differences in maximum pumpage in the basin under typical and dry long-term climate projections, respectively. The scenario 1 results indicated the maximum 10-year pumpage of the basin was about 31,300 acre-ft under typical conditions and controlling simulated seawater intrusion and drawdowns. For scenario 2, less recharge over the 10-year dry climate produced a maximum pumpage estimate of 30,000 acre-ft to control seawater intrusion and drawdowns. The larger pumpage for scenario 1 resulted in more seawater intrusion, but less total drawdown, compared to that of scenario 2. Results for scenarios 3 and 4 showed the basin’s response to management actions combined with climate projections. Both scenarios used the optimal pumping schedules and the 10-year climates from scenarios 1 and 2. The scenario 3 results showed that under minimal pumping, the basin did not fully recover to 1998 water levels within 10 years under either climate scenario. The relatively larger recharge from the typical climate resulted in less drawdown at coastal monitoring wells after the 10-year recovery period than that from the dry climate. The location of the seawater intrusion front was not appreciably different between the scenarios, however. Scenario 4 used the optimal results from scenarios 1 and 2 to produce decision-rule curves that illustrated the pumpage for each basin, given measured levels of chloride concentration or drawdown. This allowed the use of additional measurements at monitoring wells to assess future management decisions on the basis of the sensitivity of observations of drawdown and seawater intrusion to various pumping rates. Scenario 5 allowed managers to investigate the effects of short-term climate variations on optimal pumping schedules. Three specific 2-year simulations were optimized: typical-to-dry (scenario 5A), dry-to-typical (scenario 5B), and dry-to-dry (scenario 5C). The most noteable result from scenario 5 was the overall reduction in optimal pumpage for most schedules in scenario 5C, when the climate is simulated as dry-to-dry. There are also many optimal pumping schedules that produced an overall increase in waterlevels over the two-year simulation period, regardless of climatic condition. Similar to scenario 2, the scenario 5C results represents conservative yield estimates under a minimal-precipitation climatic condition.

California↗