Search USGS⌕ Search

SEARCH · Search USGS

Results for “Scientific Drilling”

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.

370 records · Page 21Linked to original sources
USGS2015

Water quality of groundwater and stream base flow in the Marcellus Shale Gas Field of the Monongahela River Basin, West Virginia, 2011-12

The Marcellus Shale gas field underlies portions of New York, Pennsylvania, Ohio, Virginia, Maryland, Tennessee, and West Virginia. Development of hydraulic fracturing and horizontal drilling technology led to extensive development of gas from the Marcellus Shale beginning about 2007. The need to identify and monitor changes in water-quality conditions related to development of the Marcellus Shale gas field prompted the U.S. Geological Survey, in cooperation with the West Virginia Department of Environmental Protection, Division of Water and Waste Management, to document water quality for comparison with water quality in samples collected at a future date. The identification of change in water-quality conditions over time is more difficult if baseline water-quality conditions have not been documented. U.S. Geological Survey personnel sampled groundwater and surface water in West Virginia’s Monongahela River Basin during 2011–12. A groundwater survey, in which 39 wells and 2 springs were sampled, was conducted during June through September 2011. A base-flow survey was conducted during July through October 2012; 50 stream sites were sampled under base-flow conditions in this survey. Because additives to hydraulic fracturing fluids are variable and decrease in flowback water over a relatively short time, water-quality analyses for this study focused on documenting the water-quality characteristics typical of water from shallow aquifers; water derived from contact with the Marcellus Shale (flowback from hydraulic fracturing or formation water); and water with constituents from conventional oil and gas development, sewage effluent, and coal-mine drainage. All samples were analyzed for field properties (water temperature, pH, specific conductance, dissolved oxygen, and turbidity), major ions, trace elements, naturally occurring radioactive materials, and stable isotopes. In addition to documenting baseline water-quality conditions for an area of shale-gas development, these data were examined for patterns in water quality. Groundwater and base-flow survey data were compared to historical data from the Monongahela River Basin in West Virginia. Additionally, groundwater- and base-flow survey samples were grouped by Marcellus Shale gas production in the subbasin in which that sampling site was located. The comparisons of data collected as part of this study with historical data identified few differences. No significant difference was found in a comparison of groundwater survey data and historical data. Base-flow survey samples differed significantly from historical data for pH, chloride, and strontium, all of which had higher concentrations in the base-flow survey samples. Differences in pH are likely related to changes in mining regulation beginning in 1977. Concentrations of chloride and strontium elevated above background concentrations may be related to saline groundwater; saline water is within 300 feet of the land surface in parts of the study area. In the comparison of base-flow survey samples grouped by shale-gas-production setting, significant differences were found for fluoride and barium. Concentrations of fluoride and barium were higher in stream subbasins with active Marcellus Shale production than in subbasins not near active Marcellus Shale production. Elevated fluoride and barium are associated with deep brines. Generally, naturally occurring radioactive materials were not found in elevated concentrations in either groundwater or base-flow samples. Only 3 samples, 2 from the groundwater survey and one from the base-flow survey, exceeded the U.S. Environmental Protection Agency maximum contaminant level for radium isotopes of 5.0 picocurie per liter for either a single isotope or a combined value of radium-226 and radium-228. Stable isotope composition indicates broad similarity among surface water, shallow groundwater, and precipitation in the region. Neither shallow groundwater nor surface water showed a marked similarity with the deep brines associated with shale gas. In most of the groundwater survey samples, 38 of 41 samples, dissolved gas profiles were similar to those previously found in samples from shallow, domestic wells in the region. This study provides a baseline of water-quality conditions in the Monongahela River Basin in West Virginia during the early phases of development of the Marcellus Shale gas field. Although not all inclusive, the results of this study provide a set of reliable water-quality data against which future data sets can be compared and the effects of shale-gas development may be determined.

West Virginia↗
USGS2015

Lithology, hydrologic characteristics, and water quality of the Arkansas River Valley alluvial aquifer in the vicinity of Van Buren, Arkansas

A study to assess the potential of the Arkansas River Valley alluvial aquifer in the vicinity of Van Buren, Arkansas, as a viable source of public-supply water was conducted by the U.S. Geological Survey in cooperation with the Little Rock, District, U.S. Army Corps of Engineers. An important study component was to identify possible changes in hydrologic conditions following installation of James W. Trimble Lock and Dam 13 (December 1969) on the Arkansas River near the study area. Data were gathered for the study in regard to the lithology, hydrologic characteristics, and water quality of the aquifer. Lithologic information was obtained from drillers’ logs of wells drilled from 1957 through 1959. Water-quality samples were collected from 10 irrigation wells and analyzed for inorganic constituents and pesticides. To evaluate the potential viability of the alluvial aquifer in the Van Buren area, these data were compared to similar stratigraphic, lithologic, and groundwater-quality data from the Arkansas River Valley alluvial aquifer at Dardanelle, Ark., where the aquifer provides a proven, productive, sole-source of public-supply water. Drillers’ logs for 59 wells in the Van Buren study area revealed well depths ranging from 25 to 52 feet (ft), with a mean depth of 42 ft. The thickness of the lower sand/gravel interval serving as the water-producing zone ranged from 5 to 47 ft, with a mean thickness of 29 ft. The presence of gravel was noted in only 4 of 59 well logs available for review from the study area. Percent sand was calculated from well logs in the study area, and these sand percentages were overlain onto an orthophotograph map to examine the areal distribution of sand percentage in relation to geomorphologic features of the flood plain in the study area. The logs denoting the greatest percent sand tend to occur in areas near to the river and on the concave (point bar) side of abandoned channels, while the lower percent sand tends to occur on the convex (channel fill and backswamp deposits) side of the abandoned channels. Comparison of hydrographs from water levels collected between 1957 and 1972 to cumulative departure from mean monthly and mean annual precipitation showed overall good fit and explained the long-term decreasing water levels from the earliest period of record through October 1967, followed by a sharp rise in water levels concurrent with rises in cumulative departure from mean monthly and mean annual precipitation. Hydrographs for four wells ranging from 0.8 to 4.5 miles upstream from the dam and potentially affected by rising river stage were compared to graphs of river stage and cumulative departure from mean monthly precipitation. Water levels for these wells showed minimal discernible effect by rising river stage following dam completion. Periods of increased precipitation compared closely to increases in water level for all hydrographs, regardless of river stage, and periods of no precipitation resulted in declining water levels, although river stage continued to slowly rise during these same periods. The Arkansas River has greater salinity than local groundwater, providing a quantitative tracer for any groundwater recharge originating from the river. Comparison of predam and postdam groundwater-chloride concentrations showed no increase in chloride concentrations after dam installation, which is consistent with hydrologic data. These data suggest that the dominant source of groundwater recharge in the Arkansas River Valley alluvial aquifer is infiltration of precipitation through proximal, coarse channel deposits, with minimal influx of river water. Groundwater-quality data collected from 10 wells in the study area indicated a calcium-bicarbonate water type. No primary drinking-water standards were exceeded for any constituents, and iron and manganese were the only constituents exceeding secondary drinking-water regulations. Six of the 10 well-water samples were analyzed for the presence of pesticides, as row-crop agriculture is the dominant land use in the study area. Six herbicide compounds and one herbicide metabolite were detected at concentrations substantially below those of the Federal primary drinking-water standards and health advisories. The hydrologic and geochemical data gathered for this study provide a qualitative assessment of the potential of the Arkansas River Valley alluvial aquifer as a source of public water supply in the Van Buren area. Results indicate minimal influx of water from the Arkansas River, and recharge to the aquifer appears to be dominantly by infiltration of precipitation through overlying alluvium. If vertical wells are used as a source of public water supply, then several wells will have to be used in combination at relatively low pumping rates and placed in areas with a greater percent sand. Use of a horizontal well configuration near the river to increase production may depend on infiltration of river water to supplement water removed from storage, especially where areas of lower permeability sediments might be encountered within the surrounding alluvium. If a poor hydraulic connection exists between the river and the alluvium, as indicated by this study, then production will depend on ample precipitation and recharge throughout the year and groundwater storage sufficient to prevent declining water levels where pumping rates exceed recharge.

Arkansas↗
USGS2023

Arsenic, chromium, uranium, and vanadium in rock, alluvium, and groundwater, Mojave River and Morongo Areas, western Mojave Desert, southern California

Trace elements within groundwater that originate from aquifer materials and pose potential public-health hazards if consumed are known as geogenic contaminants. The geogenic contaminants arsenic, chromium, and vanadium can form negatively charged ions with oxygen known as oxyanions. Uranium complexes with bicarbonate and carbonate to form negatively charged ions having aqueous chemistry similar to oxyanions. The concentrations of arsenic, chromium, uranium, and vanadium in groundwater result from the combined effects of (1) geologic abundance within aquifer materials; (2) the fraction of these elements that have weathered from and sorbed to the surfaces of mineral grains and are potentially available to groundwater; and (3) the aqueous chemistry of dissolved oxyanions in groundwater during different redox conditions and pH, both of which are affected by hydrogeology, including the length of time groundwater has been in contact with aquifer materials. Concentrations of arsenic, chromium, uranium, and vanadium were measured in samples of (1) rock, surficial alluvium, and drill cuttings using portable (handheld) X-ray fluorescence (pXRF); (2) operationally defined fractions extractable from these materials; and (3) water from wells sampled between 2000 and 2018 within the 3,500 square mile Mojave River area and Morongo area of the western Mojave Desert, southern California. Regionally, rock and surficial alluvium in the Mojave River and Morongo areas are high in arsenic, low in chromium and uranium, and near the average bulk continental crust concentration for vanadium. Locally, high chromium concentrations are present in mafic rock within the San Gabriel Mountains; high uranium concentrations are present in felsic rock within the San Bernardino Mountains; and high arsenic, uranium, and vanadium concentrations are present in extrusive (volcanic) felsic rock within uplands surrounding groundwater basins along the Mojave River downstream from Barstow, California. Elemental assemblages identified using principal component analyses (PCA) of pXRF data were used to characterize felsic, mafic, and felsic volcanic source terranes in rock, surficial alluvium, and in geologic material penetrated by selected monitoring wells drilled between 1994 and 2018. Highly felsic alluvium associated with recent deposition from the Mojave River was identified along the 90-mile length of the floodplain aquifer along the river. The thickness of these highly felsic alluvial deposits ranged from 200 feet (ft) near Victorville and near Barstow to a thin veneer about 30 ft thick downstream from Victorville and downstream portions of the floodplain aquifer within the Mojave Valley. Groundwater in the Mojave River and Morongo areas was generally oxic and alkaline (pH≥7.5). Maximum concentrations of arsenic, hexavalent chromium [Cr(VI)], uranium, and vanadium in water from as many as 498 wells sampled between 2000 and 2018 were 360, 140, 1,470, and 690 micrograms per liter (μg/L), respectively. Water from 22 percent of sampled wells exceeded the U.S. Environmental Protection Agency (EPA) maximum contaminant level (MCL) for arsenic of 10 μg/L, with arsenic concentrations commonly exceeding the MCL in water from wells east of Barstow, deep wells in the Victorville fan, and in suboxic or reduced groundwater within the floodplain aquifer. Water from about 1 percent of sampled wells had Cr(VI) concentrations greater than the California MCL for total chromium of 50 μg/L, whereas 13 percent of sampled wells had Cr(VI) concentrations greater than the former California MCL of 10 μg/L. Hexavalent chromium concentrations were highest in water from wells in the Sheep Creek alluvial fan, eroded from mafic rock in the San Gabriel Mountains, although Cr(VI) concentrations greater than the former California MCL also were present elsewhere in the study area where mafic materials or older groundwater were present. Water from about 9 percent of sampled wells exceeded the EPA MCL for uranium of 30 μg/L, with concentrations exceeding the MCL commonly associated with irrigation return from agricultural land overlying the floodplain aquifer. Water from about 7 percent of sampled wells had vanadium concentrations greater than the California notification level of 50 μg/L; most of these wells were in the Victorville fan within the Mojave River area. In general, arsenic concentrations were higher in suboxic or reduced water; chromium concentrations were higher in oxic, alkaline (pH≥7.5) water; uranium concentrations were higher in circumneutral to slightly alkaline water (pH≤7.4); and vanadium concentrations were higher in highly alkaline (pH≥8.0) water, independent of redox status. Concentrations within geologic source terranes are not the sole factor controlling the concentrations of geogenic elements in groundwater. Differences in mineral weathering, pH-dependent sorption to surface-exchange sites on mineral grains, and aqueous geochemistry (especially redox status and pH) affect geogenic element concentrations in groundwater. Consequently, the relative abundances of arsenic, Cr(VI), uranium, and vanadium in groundwater differ from their relative abundances in the average bulk continental crust and their regional abundances in rock and surficial alluvium within groundwater basins of the western Mojave Desert. Processes that control the concentrations of arsenic, chromium, uranium, and vanadium in groundwater operate at the mineral-grain and aquifer scale. At the mineral-grain scale, sequential chemical extraction data show arsenic and uranium are more available to groundwater (under specific geochemical conditions) than chromium or vanadium, which largely are unavailable within unweathered mineral grains. Additionally, chromium and vanadium form few aqueous complexes and bind tightly with iron minerals within surface coatings on mineral grains making them less available to groundwater, whereas complexation with other dissolved ions enhances the solubility of uranium and, to a lesser extent, arsenic. Complexation also increases the valence (less negative charge) and increases the size of dissolved oxyanions and uranium complexes with bicarbonate and carbonate making them less readily sorbed to aquifer materials. At the aquifer scale, hydrogeology (including isolation of water in aquifers from surface sources of recharge, older groundwater age, and long contact times between groundwater and aquifer materials) combined with geochemical processes (such as silicate weathering) to produce alkaline groundwater. Desorption from sorption sites on the surfaces of mineral grains with increasing pH increases arsenic, chromium, and vanadium concentrations in water from wells and increases Cr(VI) concentrations as long as water remains oxic. Aqueous geochemistry and concentrations of geogenic contaminants also are affected by anthropogenic activities including (1) discharge of treated municipal wastewater, which may change the redox status of groundwater; (2) return from irrigated agriculture, which may alter the chemistry of groundwater and increase the solubility of trace elements such as uranium; and (3) groundwater pumping and subsequent water-level declines, which may change the source of water yielded by wells. The quality of water imported from northern California and infiltrated from ponds for groundwater recharge may be altered by naturally present trace elements, especially uranium in areas of agricultural land use or chromium within mafic alluvium.

California↗
USGS2025

Groundwater hydrology, groundwater and surface-water interactions, water quality, and groundwater-flow simulations for the Wet Mountain Valley alluvial aquifer, Custer and Fremont Counties, Colorado, 2017–19

In 2017, the U.S. Geological Survey, in cooperation with the Upper Arkansas Water Conservancy District, began a study to provide a comprehensive analysis of the Wet Mountain Valley alluvial aquifer, Custer and Fremont Counties, Colorado. The study included collection of data pertaining to groundwater hydrology, groundwater and surface-water interactions, and water quality in the alluvial aquifer. In addition to providing foundational information on the hydrology of the alluvial aquifer, a numerical groundwater-flow model was developed to estimate the potential effects of additional storage of groundwater in the alluvial aquifer. Groundwater-level elevation data from 30 wells were used to estimate groundwater-flow directions in the alluvial aquifer, which were generally from the southwest to northeast, away from the Sangre de Cristo Mountains and towards perennial streams in the center of the valley. Although some seasonal variation was apparent in groundwater-level elevation records, no statistically significant seasonal trends were indicated. Statistically significant long-term trends were indicated in groundwater-level elevation records for 8 of the 30 wells, and of these wells with statistically significant trends, all but 1 indicated a negative trend of groundwater-level elevations. Spatial evaluation of wells with statistically significant negative groundwater-level elevation trends showed many are in areas of denser well drilling for domestic or other uses, indicating increasing groundwater use could potentially be causing groundwater-level elevation declines. There were instances of wells with no statistically significant groundwater-level elevation trends also located in areas of greater density of well completions. Additional investigations may be necessary to more fully characterize the processes responsible for negative groundwater-level elevation trends. Streamflow gain or loss calculations were completed for low flow in 2017–19 and for high flow in 2018 in nine reaches of streams within the study area. Stream reaches of the upper Texas Creek, upper Grape Creek, upper-middle Grape Creek, and Taylor Creek displayed consistent streamflow loss in each period from 2017 to 2019. These stream reaches represent long-term sources of recharge to the alluvial aquifer. Streamflow gain or loss varies through time in other stream reaches (lower Texas Creek, lower-middle Grape Creek, lower Grape Creek below Westcliffe, and lower Grape Creek above DeWeese Reservoir). The temporally variable behavior indicates these stream reaches may be sources of groundwater recharge or areas of groundwater discharge, likely depending on temporal dynamics between the elevation of the water table and the stream. Water-quality samples were collected from 10 groundwater wells and 10 stream sites during September through November 2019. All groundwater and stream samples were analyzed for major and trace elements and stable isotopes of water. A subset of groundwater samples was also analyzed for the environmental tracers sulfur hexafluoride, tritium, and noble gases. Comparison of water-quality results to U.S. Environmental Protection Agency drinking water-quality standards indicated no constituents exceeded primary standards for human health. Spatial evaluation of water quality indicated the concentrations of various constituents are likely controlled by groundwater and surface-water interactions and by spatial variability in bedrock geology underlying the alluvial aquifer. Specifically, streams shown to gain from groundwater had water chemistry constituent compositions similar to groundwater, whereas streams exiting the Sangre de Cristo Mountains tended to have compositions consistent with snowmelt. Groundwater geochemistry appeared to be partially controlled by oxidation-reduction processes and by proximity to igneous rocks in the Wet Mountains. Environmental tracers used to estimate groundwater age indicated all sampled groundwater contained tracers representing modern recharge (approximately less than 65 years old) but mixing of premodern recharge (approximately more than 65 years old) also occurs. Spatial evaluation of environmental tracers indicated large faults may be conduits for upwelling of older groundwater. No trends were observed in groundwater age with well depth, indicating all sampled wells are located within the zone of active groundwater flow. The presence of modern groundwater in wells with statistically significant negative groundwater-level elevation trends indicates groundwater storage depletions may be partially offset by capture of modern recharge. Repeated sampling of groundwater age would be necessary, however, to determine if any trends in groundwater age exist, which may indicate changing groundwater recharge, storage, or discharge. Additional investigations could also consider quantifying groundwater age in deeper wells to more fully define the depth of active groundwater flow. A numerical groundwater-flow model was developed to estimate components of the water budget, simulate groundwater and surface-water interactions, and evaluate the potential effects of aquifer storage and recovery. Simulated groundwater-level elevations from the calibrated groundwater-flow model are similar to the observed pattern of groundwater-level elevations with higher elevations in the western part of the study area along the Sangre de Cristo Mountains. Simulated water-budget components indicate most of the recharge to the alluvial aquifer is derived from streamflow losses, which is consistent with observations of losing streams along the mountain front. The largest groundwater discharge component of the alluvial aquifer was to streams in the center of the valley, where observations of stream gain or loss indicated the predominance of gaining conditions. Comparison of groundwater and surface-water interactions between the calibrated groundwater-flow model for 2000-19 (the base-case model) and a simulation including additional recharge, representing potential aquifer storage and recovery operations, indicated the additional recharge distributed throughout the area had minimal effects on streamflow in the nearby Grape Creek. An analysis of subregional groundwater budgets showed approximately 54 percent of the additional recharge flowed back to nearby Grape Creek, and the other 46 percent was distributed laterally into adjacent cells in the alluvial aquifer. The comparison of simulations and subregional water budget show the additional recharge did not substantially alter groundwater-level elevations or basin wide groundwater storage. Although the analysis of additional recharge provided in the numerical groundwater-flow model considers only one of many possible recharge scenarios, the model provides a useful tool that could be modified for various scenarios to understand potential effects of managed aquifer recharge.

Colorado↗
USGS2011

Borehole geophysical investigation of a formerly used defense site, Machiasport, Maine, 2003-2006

The U.S. Geological Survey, in cooperation with the U.S. Army Corps of Engineers, collected borehole geophysical logs in 18 boreholes and interpreted the data along with logs from 19 additional boreholes as part of an ongoing, collaborative investigation at three environmental restoration sites in Machiasport, Maine. These sites, located on hilltops overlooking the seacoast, formerly were used for military defense. At each of the sites, chlorinated solvents, used as part of defense-site operations, have contaminated the fractured-rock aquifer. Borehole geophysical techniques and hydraulic methods were used to characterize bedrock lithology, fractures, and hydraulic properties. In addition, each geophysical method was evaluated for effectiveness for site characterization and for potential application for further aquifer characterization and (or) evaluation of remediation efforts. Results of borehole geophysical logging indicate the subsurface is highly fractured, metavolcanic, intrusive, metasedimentary bedrock. Selected geophysical logs were cross-plotted to assess correlations between rock properties. These plots included combinations of gamma, acoustic reflectivity, electromagnetic induction conductivity, normal resistivity, and single-point resistance. The combined use of acoustic televiewer (ATV) imaging and natural gamma logs proved to be effective for delineating rock types. Each of the rock units in the study area could be mapped in the boreholes, on the basis of the gamma and ATV reflectivity signatures. The gamma and mean ATV reflectivity data were used along with the other geophysical logs for an integrated interpretation, yielding a determination of quartz monzonite, rhyolite, metasedimentary units, or diabase/gabbro rock types. The interpretation of rock types on the basis of the geophysical logs compared well to drilling logs and geologic mapping. These results may be helpful for refining the geologic framework at depth. A stereoplot of all fractures intersecting the boreholes indicates numerous fractures, a high proportion of steeply dipping fractures, and considerable variation in fracture orientation. Low-dip-angle fractures associated with unloading and exfoliation are also present, especially at a depth of less than 100 feet below the top of casing. These sub-horizontal fractures help to connect the steeply dipping fractures, making this a highly connected fracture network. The high variability in the fracture orientations also increases the connectivity of the fracture network. A preliminary comparison of all fracture data from all the boreholes suggests fracturing decreases with depth. Because all the boreholes were not drilled to the same depth, however, there is a clear sampling bias. Hence, the deepest boreholes are analyzed separately for fracture density. For the deepest boreholes in the study, the intensity of fracturing does not decline significantly with depth. It is possible the fractures observed in these boreholes become progressively tighter or closed with depth, but this is difficult to verify with the borehole methods used in this investigation. The fact that there are more sealed fractures at depth (observed in optical televiewer logs in some of the boreholes) may indicate less opening of the sealed fractures, less water moving through the rock, and less weathering of the fracture infilling minerals. Although the fracture orientation remained fairly constant with depth, differences in the fracture patterns for the three restoration sites indicate the orientation of fractures varies across the study area. The fractures in boreholes on Miller Mountain predominantly strike northwest-southeast, and to a lesser degree they strike northeast. The fractures on or near the summit of Howard Mountain strike predominantly east-west and dip north and south, and the fractures near the Transmitter Site strike northeast-southwest and dip northwest and southeast. The fracture populations for the boreholes on or near the summit of Howard Mountain show more variation than at the other two sites. This variation may be related to the proximity of the fault, which is northeast of the summit of Howard Mountain. In a side-by-side comparison of stereoplots from selected boreholes, there was no clear correspondence between fracture orientation and proximity to the fault. There is, however, a difference in the total populations of fractures for the boreholes on or near the summit of Howard Mountain and the boreholes near the Transmitter Site. Further to the southwest and further away from the fault, the fractures at the Transmitter Site predominantly strike northeast-southwest and northwest-southeast.Heat-pulse flowmeter (HPFM) logging was used to identify transmissive fractures and to estimate the hydraulic properties along the boreholes. Ambient downflow was measured in 13 boreholes and ambient upflow was measured in 9 boreholes. In nine other bedrock boreholes, the HPFM did not detect measurable vertical flow. The observed direction of vertical flow in the boreholes generally was consistent with the conceptual flow model of downward movement in recharge locations and upward flow in discharge locations or at breaks in the slope of land surface. Under low-rate pumping or injection rates [0.25 to 1 gallon per minute (gal/min)], one to three inflow zones were identified in each borehole. Two limitations of HPFM methods are (1) the HPFM can only identify zones within 1.5 to 2 orders of magnitude of the most transmissive zone in each borehole, and (2) the HPFM cannot detect flow rates less than 0.010 + or - 0.005 gal/min, which corresponds to a transmissivity of about 1 foot squared per day (ft2/d). Consequently, the HPFM is considered an effective tool for identifying the most transmissive fractures in a borehole, down to its detection level. Transmissivities below that cut-off must be measured with another method, such as packer testing or fluid-replacement logging. Where sufficient water-level and flowmeter data were available, HPFM results were numerically modeled. For each borehole model, the fracture location and measured flow rates were specified, and the head and transmissivity of each fracture zone were adjusted until a model fit was achieved with the interpreted ambient and stressed flow profiles. The transmissivities calculated by this method are similar to the results of an open-hole slug test; with the added information from the flowmeter, however, the head and transmissivity of discrete zones also can be determined. The discrete-interval transmissivities ranged from 0.16 to 330 ft2/d. The flowmeter-derived open-hole transmissivity, which is the combined total of each of the transmissive zones, ranged from 1 to 511 ft2/d. The whole-well open-hole transmissivity values determined with HPFM methods were compared to the results of open-hole hydraulic tests. Despite the fact that the flowmeter-derived transmissivities consistently were lower than the estimates derived from open-hole hydraulic tests alone, the correlation was very strong (with a coefficient of determination, R2, of 0.9866), indicating the HPFM method provides a reasonable estimate of transmissivities for the most transmissive fractures in the borehole. Geologic framework, fracture characterization, and estimates of hydraulic properties were interpreted together to characterize the fracture network. The data and interpretation presented in this report should provide information useful for site investigators as the conceptual site groundwater flow model is refined. Collectively, the results and the conceptual site model are important for evaluating remediation options and planning or implementing the design of a well field and borehole completions that will be adequate for monitoring flow, remediation efforts, groundwater levels, and (or) water quality. Similar kinds of borehole geophysical logging (specifically the borehole imaging, gamma, fluid logs, and HPFM) should be conducted in any newly installed boreholes and integrated with interpretations of any nearby boreholes. If boreholes are installed close to existing or other new boreholes, cross-hole flowmeter surveys may be appropriate and may help characterize the aquifer properties and connections between the boreholes.

Maine↗
USGS2013

Geology, water-quality, hydrology, and geomechanics of the Cuyama Valley groundwater basin, California, 2008--12

To assess the water resources of the Cuyama Valley groundwater basin in Santa Barbara County, California, a series of cooperative studies were undertaken by the U.S. Geological Survey and the Santa Barbara County Water Agency. Between 2008 and 2012, geologic, water-quality, hydrologic and geomechanical data were collected from selected sites throughout the Cuyama Valley groundwater basin. Geologic data were collected from three multiple-well groundwater monitoring sites and included lithologic descriptions of the drill cuttings, borehole geophysical logs, temperature logs, as well as bulk density and sonic velocity measurements of whole-core samples. Generalized lithologic characterization from the monitoring sites indicated the water-bearing units in the subsurface consist of unconsolidated to partly consolidated sand, gravel, silt, clay, and occasional cobbles within alluvial fan and stream deposits. Analysis of geophysical logs indicated alternating layers of finer- and coarser-grained material that range from less than 1 foot to more than 20 feet thick. On the basis of the geologic data collected, the principal water-bearing units beneath the monitoring-well sites were found to be composed of younger alluvium of Holocene age, older alluvium of Pleistocene age, and the Tertiary-Quaternary Morales Formation. At all three sites, the contact between the recent fill and younger alluvium is approximately 20 feet below land surface. Water-quality samples were collected from 12 monitoring wells, 27 domestic and supply wells, 2 springs, and 4 surface-water sites and were analyzed for a variety of constituents that differed by site, but, in general, included trace elements; nutrients; dissolved organic carbon; major and minor ions; silica; total dissolved solids; alkalinity; total arsenic and iron; arsenic, chromium, and iron species; and isotopic tracers, including the stable isotopes of hydrogen and oxygen, activities of tritium, and carbon-14 abundance. Of the 39 wells sampled, concentrations of total dissolved solids and sulfate from 38 and 37 well samples, respectively, were greater than the U.S. Environmental Protection Agency’s secondary maximum contaminant levels. Concentrations greater than the maximum contaminant levels for nitrate were observed in five wells and were observed for arsenic in four wells. Differences in the stable-isotopic values of hydrogen and oxygen among groundwater samples indicated that water does not move freely between different formations or between different zones within the Cuyama Valley. Variations in isotopic composition indicated that recharge is derived from several different sources. The age of the groundwater, expressed as time since recharge, was between 600 and 38,000 years before present. Detectable concentrations of tritium indicated that younger water, recharged since the early 1950s, is present in parts of the groundwater basin. Hydrologic data were collected from 12 monitoring wells, 56 domestic and supply wells, 3 surface-water sites, and 4 rainfall-gaging stations. Rainfall in the valley averaged about 8 inches annually, whereas the mountains to the south received between 12 and 19 inches. Stream discharge records showed seasonal variability in surface-water flows ranging from no-flow to over 1,500 cubic feet per second. During periods when inflow to the valley exceeds outflow, there is potential recharge from stream losses to the groundwater system Water-level records included manual quarterly depth-to-water measurements collected from 68 wells, time-series data collected from 20 of those wells, and historic water levels from 16 wells. Hydrographs of the manual measurements showed declining water levels in 16 wells, mostly in the South-Main zone, and rising water levels in 14 wells, mostly in the Southern Ventucopa Uplands. Time-series hydrographs showed daily, seasonal, and longer-term effects associated with local pumping. Water-level data from the multiple-well monitoring sites indicated seasonal fluctuations as great as 80 feet and water-level differences between aquifers as great as 40 feet during peak pumping season. Hydrographs from the multiple-well groundwater monitoring sites showed vertical hydraulic gradients were upward during the winter months and downward during the irrigation season. Historic hydrographs showed water-level declines in the Southern-Main, Western Basin, Caliente Northern-Main, and Southern Sierra Madre zone ranging from 1 to 7 feet per year. Hydrographs of wells in the Southern Ventucopa Uplands zone showed several years with marked increases in water levels that corresponded to increased precipitation in the Cuyama Valley. Investigation of hydraulic properties included hydraulic conductivity and transmissivity estimated from aquifer tests performed on 63 wells. Estimates of horizontal hydraulic conductivity ranged from about 1.5 to 28 feet per day and decreased with depth. The median estimated hydraulic conductivity for the older alluvium was about five times that estimated for the Morales Formation. Estimates of transmissivity ranged from 560 to 163,400 gallons per day per foot and decreased with depth. The median estimated transmissivity for the younger alluvium was about three times that estimated for the older alluvium. Geomechanical analysis included land-surface elevation changes at five continuously operating global positioning systems (GPS) and land-subsidence detection at five interferometric synthetic aperture radar (InSAR) reference points. Analysis of data collected from continuously operating GPS stations showed the mountains to the south and west moved upward about 1 millimeter (mm) annually, whereas the station in the center of the Southern-Main zone moved downward more than 7 mm annually, indicating subsidence. It is likely that this subsidence is inelastic (permanent) deformation and indicates reduced storage capacity in the aquifer sediments. Analysis of InSAR data showed local and regional changes that appeared to be dependent, in part, on the time span of the interferogram, seasonal variations in pumping, and tectonic uplift. Long-term InSAR time series showed a total maximum detected subsidence rate of approximately 12 mm per year at one location and approximately 8 mm per year at a second location, while short-term InSAR time series showed maximum subsidence of about 15 mm at one location and localized maximum uplift of about 10 mm at another location.

California↗
USGS2024

Hydrologic analysis of an earthen embankment dam in southern Westchester County, New York

In 2001, the New York City Department of Environmental Protection installed 25 wells on the southern embankment of the Hillview Reservoir in Westchester County in an unsuccessful attempt to locate the source of a large seep (seep A) that began flowing continuously in 1999. In 2005, the U.S. Geological Survey began a cooperative study with the NYCDEP to characterize the hydrology of the local groundwater system and identify potential sources of seep A and other seeps on the embankment. At least two groundwater-flow zones—one shallow and the other deep—overlie the bedrock at the Hillview Reservoir in southern Westchester County, New York. Analyses of slug tests of wells drilled into the southern embankment of the reservoir were used to determine the three-dimensional distribution of hydraulic conductivity of the embankment materials. The wells with the minimum and maximum hydraulic conductivity values are in the deep saturated zone on the southern embankment, where hydraulic conductivity ranges from 0.0012 to 2 feet per day. Hydraulic conductivity ranges from 0.0026 to 1 foot per day in the shallow saturated zone and from 0.021 to 0.27 foot per day in the toe of the embankment. A hydraulic conductivity of 0.016 foot per day was determined for one toe well partially screened in the crystalline-bedrock aquifer. In 2005, the U.S. Geological Survey began a cooperative study with New York City Department of Environmental Protection to characterize the local groundwater-flow system and identify potential sources of seeps on the southern embankment of the Hillview Reservoir in southern Westchester County, New York. Long-term hydrologic data indicated that water levels trended downward in 29 of 41 sites, including the reservoir basin that was monitored during the 12-year study period; data from a National Weather Service precipitation gage at Central Park indicated annual precipitation also trended downward during the same 12-year period. Of the seven wells in which water levels trended upward during the study, two of the wells are on the west side of the southern embankment, proximal to a major water supply conduit, whereas the five remaining wells are screened in the toe. These data indicate an increasing hydrostatic pressure within the deep system and the toe of the dam, which could result in future seeps on the southern embankment near these wells. Results of 11 suspended-sediment samples collected from seeps along the southern embankment at 234.1- and 221.6-feet elevation, and another drainage outflow point between 2007 and 2015 indicate a poor correlation between suspended-sediment concentration and discharge. From the flowing seep at 234.1 feet, suspended-sediment concentrations ranged from 1 milligram per liter at a flow of 2.6 gallons per minute (that is, 1 milligram per 0.26 gallons) during March 2008 to 16 milligrams per liter at 12 gallons per minute during July 2014. At about 12 gallons per minute discharge, suspended-sediment concentration from samples collected at that seep during different sampling events, ranged from 3 to 16 milligrams per liter. From the seep at 221.6 feet elevation, the suspended sediment concentration was 2 milligrams per liter at a discharge of 3.4 gallons per minute and 2 milligrams per liter at a discharge of 1.1 gallons per minute. Only one sample was collected at the drainage outflow point, for which the suspended sediment concentration was 2 milligrams per liter at a discharge of 2.4 gallons per minute. Anomalously high-water levels were recorded in deep-system wells between June 5, 2013, and January 14, 2014. The period for the increase and the decrease back to more typical water-level elevations occurred rapidly during a 13-hour period in each instance. The sudden and rapid changes, in addition to the spatial distribution of magnitude of water-level response indicate that leaky water infrastructure was the source of recharge to the affected wells. A major water supply conduit was drained for repairs between July 7 and 10, 2010. The seeps indicated an immediate response and a substantial hydraulic connection to the water supply conduit. Approximately 10.5 hours after the water supply conduit was drained, flow from a seep on the southern embankment decreased from about 20 gallons per minute to less than 1 gallon per minute. This seep is located at about the same elevation and within the vicinity of the water supply conduit. A travel-time of about 10.5 hours from the source to the seep at 234.1 feet elevation was estimated from the dewatering timeline. During the 3-month shutdown of the water supply conduit, the previously flowing seeps remained dry until precipitation resulted in discharge of about 0.7 gallon per minute at the higher elevation seep, indicating a minor contribution from precipitation to the total seepage discharge. Discharge from the seeps resumed almost immediately coincident with the refilling of the water supply conduit, supporting the hydraulic connection observations during the drainage stage. In addition, during the refilling of the water supply conduit on September 21, 2010, a new seep (I) was observed on the southern embankment. Discharge from this new seep remained relatively constant until it became inaccessible under construction stone from subsequent embankment repairs by the New York City Department of Environmental Protection. Precipitation after the refilling stage of the shutdown seemed to have induced a rise in water levels in the toe wells and an increase in discharge from the seep at 234.1 feet elevation. The post shutdown discharge was less than 12 gallons per minute, compared to a discharge of about 20 gallons per minute before the repairs. The lower discharge rate measured during the period of historically higher discharge rates for the fall season indicates that the repair of the major water supply conduit may have contributed to a reduced discharge from the seeps. There were no definitive responses to the shutdown in any of the wells near the major water supply conduit. The more transmissive deep system of the southern embankment near the major water supply conduit and its associated infrastructure seems to be the preferential flow path for leaking infrastructure. The wells screened in this system showed a response during the deep system anomaly and have some of the highest hydraulic conductivities of the tested wells. All the seeps are in the elevation range of the deep system from approximately the crystalline bedrock surface around 200 feet elevation to the contact between the deep and shallow saturated zones of the reservoir at about 250 feet elevation.

New York↗
USGS2020

Geochemical and mineralogical study of the Red Mountain porphyry copper-molybdenum deposit and vicinity, Santa Cruz County, Arizona

The Red Mountain porphyry copper-molybdenum deposit (Cu-Mo deposit or PCD) is located in the northern part of the Patagonia Mountains, Santa Cruz County, Arizona. Extensive core drilling has delineated a large, deep-seated, structurally intact mineral system that extends from the present surface to depths of more than 1,765 meters. This system is hosted in a thick complex of predominantly felsic to andesitic volcanic rocks of the Cretaceous Period. This complex was intruded by scattered bodies of the Tertiary Period that are predominantly quartz monzonite porphyry; no major associated source intrusion has yet been found at depth. A total of 818 samples of core were analyzed for as many as 44 elements. The abundances and distributions at depth of at least 17 of these elements (silver [Ag], arsenic [As], gold [Au], boron [B], bismuth [Bi], copper [Cu], mercury [Hg], potassium [K], molybdenum [Mo], lead [Pb], sulfur [S], antimony [Sb], tin [Sn], tellurium [Te], thallium [Tl], tungsten [W], and zinc [Zn]) are related mostly to events that generated the Red Mountain system. Many of these same samples were also analyzed by X-ray diffraction for a suite of minerals. The multielement and mineralogical analyses of the core samples provide important information about the concentrations, associations, and distributions of select elements and minerals, including zoning patterns that may not be apparent from visual examination of core samples. The distributions of selected elements and minerals in these samples reveal an unusually complete mineral system that extends from a typical PCD with potassic alteration at depth to peripheral zones of phyllic and advanced argillic alteration as well as a copper-rich supergene enriched zone and the remnants of a leached cap. R-mode factor analysis was run with 34 elements for a set of samples from the deep part of the hypogene Cu-Mo deposit and another set from the part of the supergene zone with the highest copper enrichment. For the hypogene zone dataset, five factors are related to the PCD: (1) Ag, Cu, Mo, S, and Te; (2) As, B, Hg, and Sb; (3) Au and sodium (Na); (4) manganese (Mn), Pb, and Zn; and (5) K and Tl. For the supergene dataset, the deposit-related factors include (1) Cu, Mo, S, and Te; (2) Ag, As, Hg, Pb, Sb, and Tl; (3) Au and Na; and (4) K and rubidium (Rb). The changes in element associations between the two datasets indicate that some of these new associations are a result of formation of several suites of hypogene minerals in the deep part of the deposit and different hypogene mineral suites in the peripheral part of the deposit. Some changes may be because of the effects of supergene processes. Zones containing deposit-related elements and minerals common to many PCDs are present at Red Mountain. These zones include a crude, inverted cup-shaped shell containing anomalous copper accompanied by high concentrations of Ag, Au, K, Mo, total S, sulfate S, Sb, Te, and Tl, as well as local concentrations of As, B, Hg, Pb, and Zn. Hydrothermal minerals spatially associated with the deep hypogene Cu-Mo deposit include chalcopyrite, molybdenite, pyrite, plagioclase, orthoclase, biotite, magnetite, calcite, quartz, and anhydrite. Many of the hydrothermally deposited elements that are spatially related to the deposit are also concentrated in zones above the deep part of the deposit, including Ag, As, K, Pb, Sb, Te, Tl, and Zn. These elements are concentrated either (1) in generally wide, flat zones present in the upper part of the system or (2) in crudely arcuate peripheral zones found mainly in the middle part of the system and surrounding the deep part of the deposit. Near-surface, restricted hypogene anomalies are present for bismuth, mercury, tin, and tungsten. The upper part of the deposit has been subjected to supergene enrichment and weathering. Deposit-related elements that remain anomalous in this area include Ag, As, Au, B, Bi, cobalt (Co), Cu, Hg, Mo, Pb, S, Sb, Sn, Te, Tl, uranium (U), and W. These positive concentrations indicate that, with the exception of copper and possibly mercury and uranium, these elements had relatively low chemical mobilities in the supergene enrichment and later weathering environments at Red Mountain. Most may have been deposited during one or more hypogene events and then redistributed locally during later events. Zinc is the only deposit-related element that has clearly been depleted as a result of supergene and (or) weathering events. Minerals that are common in the unweathered upper part of the system include chalcocite, pyrite, quartz, sericite, alunite, and pyrophyllite, as well as less common covellite, enargite, tennantite, tourmaline, barite, anglesite, and other sulfide or sulfate minerals. Subsequent to formation of the Red Mountain Cu-Mo deposit and supergene enrichment, chemical weathering produced an area of pervasive hematite and other iron oxides in the near-surface part of the deposit to form a leached cap. These iron-rich minerals formed primarily as a result of the oxidation of pyrite. This event was accompanied by losses of cobalt, mercury, magnesium, and zinc, as well as destruction of sericite, plagioclase, pyrite, clay minerals, and pyrophyllite. A total of 122 rock samples, 119 soil samples, and samples of three plant species (57 mesquite, 108 oak, and 68 juniper) were collected over and around Red Mountain. For the rock and soil samples, the distributions of anomalous Ag, As, Bi, Cu, Fe, Mo, Pb, Sb, Te, and Tl best delineated the exposed part of the deposit. The highest concentrations of many of these elements are centered on one or both of two main areas with exposures of quartz monzonite porphyry. The high concentrations of arsenic in the deposit area (as much as 390 parts per million (ppm) in rock and 1,500 ppm in soil) and of lead (as much as 2,370 ppm in rock and 1,490 ppm in soil) are particularly noteworthy. The concentrations of various elements in the plant ash vary widely among the three species and are species dependent. Many of the deposit-related elements are either nonessential for plant growth or are considered toxic at certain concentration ranges. In spite of this, the distributions of potentially toxic Ag, As, Bi, Cd, Cu, Mo, Pb, Sb, selenium (Se), and Zn produce deposit-related anomalies for one or more of the three species. Vegetation sampling offered no advantage over rock or soil sampling as an exploration tool. From an environmental standpoint, however, the plant analyses provide baseline data for both essential and nonessential elements that might be useful, for example, for selecting native plant species for revegetating mine waste areas. The exposed part of the Red Mountain deposit has not been greatly disturbed as a result of mining and other activities. However, some of the rock, soil, and plant samples that were collected near the Harshaw Creek and Alum Gulch drainages, which are peripheral to Red Mountain, are also anomalous for various deposit-related elements. These anomalies are probably the result of dispersion of stream sediments contaminated with material from past mining.

Arizona↗
USGS2021

Geologic and hydrogeologic characteristics of the White River Formation, Lance Formation, and Fox Hills Sandstone, northern greater Denver Basin, southeastern Laramie County, Wyoming

In cooperation with the Wyoming State Engineer’s Office, the U.S. Geological Survey studied the geologic and hydrogeologic characteristics of Cenozoic and Upper Cretaceous strata at a location in southeastern Laramie County within the Wyoming part of the Cheyenne Basin, the northern subbasin of the greater Denver Basin. The study aimed to improve understanding of the aquifers/aquifer systems in these strata, motivated in part by declining groundwater levels and interest in exploring future groundwater supplies. Based on detailed geologic characterization using information obtained by drilling and coring a 960-foot-(ft) deep exploratory borehole, and comparisons with previously published descriptions, identified Cenozoic lithostratigraphic units included 40 ft of Quaternary older alluvial fan deposits consisting of an unconsolidated mixture of sand and gravel with lesser quantities of silt and clay in varying proportions and the underlying 407.3-ft-thick White River Formation of late Eocene-Oligocene age consisting largely of mudrocks with sparse thin beds of sandstone, muddy gravel, and conglomeratic mudrocks. Identified Upper Cretaceous lithostratigraphic units included the 351.6-ft-thick Lance Formation, consisting of terrestrial sedimentary rocks including mudrocks (muddy shale and silty and sandy shale, siltstone, claystone, and mudstone) interbedded with much smaller quantities of very fine- to medium-grained muddy and silty sandstone and coal; the 79.6-ft-thick Fox Hills Sandstone, consisting of a transitional marine sequence of muddy or silty sandstone present in five individual beds; and 86.7 ft of the upper transition member of the Pierre Shale, consisting largely of marine sedimentary rocks such as muddy shale. Beds of the upper and lower Fox Hills Sandstone were separated by tongues of the Lance Formation and upper transition member of the Pierre Shale, respectively. The White River hydrogeologic unit, consisting of the entire White River Formation or Group at the study site, did not contain any substantial secondary permeability features in the mudrocks that composed almost all the unit. A monitoring well (BR–1) was completed in the White River aquifer with the well screen open to the only coarse-grained unit (muddy sandstone) that had sufficient thickness and permeability to be considered as an aquifer. Sampling of the well for a broad suite of constituents indicated groundwater generally was of excellent quality except dissolved arsenic was detected at a concentration greater than the U.S. Environmental Protection Agency (EPA) Maximum Contaminant Level, and dissolved sodium was measured at a concentration greater than several EPA Drinking Water Advisory Levels (DWAs) for the constituent. Well development, well purging for groundwater sampling, and calculated aquifer properties indicated the sandstone aquifer screened by monitoring well BR–1 was not very productive. Analysis of the well water-level responses in BR–1 to atmospheric loading and Earth tides indicated the responses were consistent with a confined-aquifer response with wellbore-storage effects. Hydraulic properties estimated based on these responses yielded values of hydraulic conductivity ( K , 0.057 foot per day [ft/d]), specific storage ( Ss , 1.6×10 −6 per foot [ft −1 ]) and porosity ( n , 0.43). Water levels filtered to remove the effects of atmospheric loading and Earth tides indicated an upward trend (+1.13 foot per year [ft/yr]) during the period analyzed, September 5, 2014, to September 30, 2017. Lithologic characteristics of the Lance hydrogeologic unit, consisting of the entire Lance Formation at the study site, indicated a potential aquifer in a “sandy” interval in the upper part of the unit. Most of the Lance hydrogeologic unit below the “sandy” interval consisted of various low-permeability lithologies unlikely to yield substantial quantities of water. This lower part of the hydrogeologic unit likely functions as a confining unit separating the underlying Lance-Fox Hills aquifer. A geologic cross section constructed for this study indicated fine-grained sediments composed most of the Lance Formation/hydrogeologic unit not only at the study location, but also throughout southern Laramie County along the line of section and throughout the Wyoming and Colorado parts of the Cheyenne Basin. A monitoring well (LN–1) completed in a sandstone bed in the “sandy” interval of the Lance hydrogeologic unit produced a mean of about 23 gallons per minute (gal/min) during well development, indicating sandstone beds can form moderately productive confined subaquifers in this part of the hydrogeologic unit. Analysis of the well water-level responses in well LN–1 to atmospheric loading and Earth tides indicated the responses were consistent with a confined-aquifer response. Hydraulic properties estimated based on these responses yielded values for a lower bounding K of 0.60 ft/d, Ss of 1.6×10 −6 ft −1 , and n of 0.38. Water levels filtered to remove the effects of atmospheric loading and Earth tides indicated a downward trend (−0.86 ft/yr) during the period analyzed (November 8, 2014, to September 30, 2017). Analyses for a broad suite of constituents in samples from well LN–1 indicated groundwater quality generally was excellent, although dissolved sodium was measured at a concentration greater than two EPA DWA levels for the constituent. Because of the absence of any overlying or intertonguing sandstone beds belonging to the lower/basal part of the Lance Formation, the Lance-Fox Hills aquifer at the study site consisted only of the five sandstone beds of the Fox Hills Sandstone. The cross section constructed for this study illustrated how the Fox Hills Sandstone, and thus, most of the Lance-Fox Hills aquifer, consists of a series of sandstone bodies that overlap (shingle) upward to the east across southern Laramie County. These bodies collectively form a fairly continuous body of sandstone, thus potentially forming an areally extensive aquifer across southern Laramie County, and by extension, throughout most of the formation’s extent in the Wyoming part of the Cheyenne Basin, as is the case in the Colorado part of the basin. A monitoring well (FH–1) completed in part of the thickest sandstone bed of the Lance-Fox Hills aquifer was moderately to highly productive and easily produced 25 to 30 gal/min after development. Substantially larger water production rates likely could be obtained by penetrating the full thickness of this bed and by completing a well open to the other overlying and underlying sandstone beds of the aquifer. Analysis of the water-level responses in well FH–1 to atmospheric loading and Earth tides indicated the responses were consistent with a confined-aquifer response. Hydraulic properties computed based on these responses yielded values for a lower bounding estimate for K of 0.26 ft/d, for Ss of 1.0×10 −6 ft −1 , and for n of 0.41. Water levels filtered to remove the effects of atmospheric loading and Earth tides indicated a downward trend (−1.74 ft/yr) during the period analyzed, December 19, 2014, to September 30, 2017. Sampling of monitoring well FH–1 and two production wells completed in the Fox Hills Sandstone in other parts of Laramie County indicated groundwater quality generally is excellent, although pH exceeded a recommended EPA aesthetic drinking-water standard (Secondary Maximum Contaminant Level) in two of three sampled wells, total dissolved solids concentrations exceeded the Secondary Maximum Contaminant Level in one of the two sampled production wells, and dissolved sodium was measured in all three sampled wells at a concentration greater than two EPA DWA levels for the constituent. The Wyoming Class II agricultural (irrigation) sodium adsorption ratio standard of 8 was exceeded in all three sampled wells, indicating these waters are not suitable for irrigation use. Computed vertical hydraulic gradients indicated a strong potential for downward flow throughout the groundwater system at the study site, including from the low-yielding aquifer in the upper White River Formation/hydrogeologic unit (monitoring well BR–1) to the sandstone subaquifer in the Lance Formation/hydrogeologic unit (monitoring well LN–1), and from the Lance subaquifer (monitoring well LN–1) to the sandstone bed/aquifer that composes much of the Lance-Fox Hills aquifer thickness at the study site (monitoring well FH–1). However, large hydraulic-head differences between wells indicated high resistance to vertical flow attributable to the low vertical hydraulic conductivity of intervening strata, which consisted almost entirely of low-permeability mudrocks. The confined nature of the sandstone aquifers monitored by the various wells coupled with dissimilarities between groundwater-level fluctuations and trends in groundwater levels indicated downward flow through the intervening strata (primarily mudrocks in the various lithostratigraphic/hydrogeologic units) between the examined sets of wells likely was small.

Wyoming↗
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↗