USGS ScienceSearch

SEARCH · USGS Science

Results for “Open-File Report”

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.

756 records · Page 33Linked to original sources

Per- and polyfluoroalkyl substances in waters associated with oil and gas development in the Denver Basin

Use of per- and polyfluoroalkyl substances (PFAS) in the petroleum industry could be a cause for concern due to the large volumes of produced water (PW) generated during oil and gas extraction, the reuse of these wastes in water-stressed regions, and adverse health outcomes related to PFAS exposures. However, PW PFAS characterization is nearly absent in the literature, and hydraulic fracturing (HF) chemical disclosures often omit the identities of additives as proprietary. Here we evaluate PFAS in PW samples from three petroleum wells in the Denver Basin during their first year of production. Total concentrations of targeted PFAS (Σ 40 PFAS) were < 35 ng/L in PW samples, with short-chain PFAS like perfluorobutanoic acid persisting throughout the sampled duration. Analysis of freshwater inputs for hydraulic fracturing (Σ 40 PFAS ~ 113 ng/L) and mixed fracture fluid (Σ 40 PFAS ~ 69 ng/L) indicated much of the targeted PFAS content was derived from the input water, and not from HF additives, however samples subjected to oxidation indicated the presence of PFAS precursors that would not be detected by targeted analysis. This study highlights that while PFAS content is low in the studied PWs, the potential for redistribution of PFAS in the environment may be a consideration for reuse applications.

Colorado, Nebraska, Wyoming

Surficial geologic map of the Owlshead Mountains 30' x 60' quadrangle, Inyo and San Bernardino Counties, California

The surficial geologic map of the Owlshead Mountains 30' x 60' quadrangle depicts the distribution and characteristics of surficial-deposit materials and neotectonic deformation for an area of approximately 5,000 square kilometers (km 2 ) located in the western Basin and Range Province of eastern California. The map represents a new compilation of the surficial geology that encompasses deposits within the late Pliocene to Quaternary. The map is based primarily on new mapping conducted between 2001 and 2009. Map compilation was supported by field observations distributed across the map area, combined with reference to several published and unpublished mapping sources that mostly emphasized neotectonic deformation. The surficial-deposit units included in the map follow a classification scheme that systematically denotes depositional process, relative age, and any secondary sedimentologic or morphologic characteristics. Identification, correlation, and age estimation of map units are based primarily on the relative degree of development of certain time-dependent characteristics such as surface morphology, including local dissection and surface preservation, surface clast modification, and degree of soil development; these characteristics are implicitly incorporated into unit designations. The map represents a detailed and regionally uniform synthesis of the late Neogene geology for this large area that provides a framework applicable to many interpretative studies, such as regional patterns of deposition and dissection; surface drainage development and evolution; and the distribution, style, and timing of neotectonic deformation.

California

Episodic sediment-discharge events in Cascade Springs, southern Black Hills, South Dakota

Cascade Springs is a group of artesian springs in the southern Black Hills, South Dakota, with collective flow of about 19.6 cubic feet per second. Beginning on February 28, 1992, a large discharge of red suspended sediment was observed from two of the six known discharge points. Similar events during 1906-07 and 1969 were documented by local residents and newspaper accounts. Mineralogic and grain-size analyses were performed to identify probable subsurface sources of the sediment. Geochemical modeling was performed to evaluate the geochemical evolution of water discharged from Cascade Springs. Interpretations of results provide a perspective on the role of artesian springs in the regional geohydrologic framework. X-ray diffraction mineralogic analyses of the clay fraction of the suspended sediment were compared to analyses of clay-fraction samples taken from nine geologic units at and stratigraphically below the spring-discharge points. Ongoing development of a subsurface breccia pipe(s) in the upper Minnelusa Formation and/or Opeche Shale was identified as a likely source of the suspended sediment; thus, exposed breccia pipes in lower Hell Canyon were examined. Upper Minnelusa Formation breccia pipes in lower Hell Canyon occur in clusters similar to the discrete discharge points of Cascade Springs. Grain-size analyses showed that breccia masses lack clay fractions and have coarser distributions than the wall rocks, which indicates that the red, fine-grained fractions have been carried out as suspended sediment. These findings support the hypothesis that many breccia pipes were formed as throats of abandoned artesian springs. Geochemical modeling was used to test whether geochemical evolution of ground water is consistent with this hypothesis. The evolution of water at Cascade Springs could not be suitably simulated using only upgradient water from the Minnelusa aquifer. A suitable model involved dissolution of anhydrite accompanied by dedolomitization in the upper Minnelusa Formation, which is caused by upward leakage of relatively fresh water from the Madison aquifer. The anhydrite dissolution and dedolomitization account for the net removal of minerals that would lead to breccia pipe formation by gravitational collapse. Breccia pipes in the lower Minnelusa Formation are uncommon; however, networks of interconnected breccia layers and breccia dikes are common. These networks, along with vertical fractures and faults, are likely pathways for transmitting upward leakage from the Madison aquifer. It is concluded that suspended sediment discharged at Cascade Springs probably results from episodic collapse brecciation that is caused by subsurface dissolution of anhydrite beds and cements of the upper Minnelusa Formation, accompanied by replacement of dolomite by calcite. It is further concluded that many breccia pipes probably are the throats of artesian springs that have been abandoned and exposed by erosion. The locations of artesian spring-discharge points probably have been shifting outwards from the center of the Black Hills uplift, essentially keeping pace with regional erosion over geologic time. Thus, artesian springflow probably is a factor in controlling water levels in the Madison and Minnelusa aquifers, with hydraulic head declining over geologic time, in response to development of new discharge points. Development of breccia pipes as throats of artesian springs would greatly enhance vertical hydraulic conductivity in the immediate vicinity of spring-discharge points. Horizontal hydraulic conductivity in the Minnelusa Formation also may be enhanced by dissolution processes related to upward leakage from the Madison aquifer. Potential processes could include dissolution resulting from leakage in the vicinity of breccia pipes that are abandoned spring throats, active spring discharge, development of subsurface breccias with no visible surface expression or spring discharge, as well as general areal leakage from the Madison aquifer into the Minnelusa Formation.

South Dakota

Bathymetric and velocimetric surveys at highway bridges crossing the Missouri and Mississippi Rivers on the periphery of Missouri, June 13–22, 2022

Bathymetric and velocimetric data were collected by the U.S. Geological Survey, in cooperation with the Missouri Department of Transportation, near seven bridges at six highway crossings of the Missouri and Mississippi Rivers on the periphery of Missouri from June 13–22, 2022. A multibeam echosounder mapping system was used to obtain channel-bed elevations for river reaches about 1,640 feet longitudinally and generally extending laterally across the active channel from bank to bank during minor flood-flow conditions. These surveys provided channel geometry and hydraulic conditions at the time of the surveys and provided characteristics of scour holes that may be useful in developing or verifying predictive guidelines or equations for computing potential scour depth. These data also may be useful to the Missouri Department of Transportation as a minor flood-flow assessment of the bridges for stability and integrity issues with respect to bridge scour during floods. Bathymetric data were collected around every in-channel pier. Scour holes were present at most piers for which bathymetry could be obtained, except those on banks or surrounded by riprap. Occasionally, scour holes were minor and difficult to discern from nearby dunes and ripples. All bridge sites in this study were surveyed and documented in previous studies. Although partial exposure of substructural support elements was observed at several piers, at most sites the exposure most likely is minimal compared to the overall substructure that remains buried in bed material at these piers. The notable exceptions are piers 12 and 13 at structure L0135 on State Highway 51 at Chester, Illinois, where the bedrock material was fully exposed around the piers. The average difference between the bathymetric surfaces between 2022 and 2018 varied from 0.41 foot higher to 1.86 feet lower. Between 2022 and 2014, the average difference between the bathymetric surfaces varied from 1.02 feet higher to 4.69 feet lower. Only the two sites on the Missouri River and the Caruthersville site were surveyed in 2011; for those sites, the average difference between the bathymetric surfaces varied from 5.83 feet higher to 1.34 feet lower. The most substantial overall net gain of sediment in a reach was between 2011 and 2022 at structure A1700 near Caruthersville, Mo. (site 38). This result was expected because structure A1700 is downstream from the confluences of the Missouri and Ohio Rivers, and therefore subject to the largest streamflows, the largest streamflow fluctuations, and the most substantial sediment flux, as has historically been observed at this site. The presence of riprap blankets, pier size and nose shape, and alignment to flow had a substantial effect on the size of the scour hole observed for a given pier. Piers that were surrounded by riprap blankets had scour holes that were substantially smaller (to nonexistent) compared to piers at which no rock or riprap were present. New riprap blankets were surveyed at pier 3 of structure L0098 at Brownville, Nebraska, and at piers 15–18 of structure A1700 near Caruthersville, Mo., that effectively mitigated the scour holes historically observed at these piers. Narrow piers having round or sharp noses that were aligned with flow often had scour holes that were difficult to discern from nearby bed features, whereas piers having wide or blunt noses resulted in larger, deeper scour holes. Several of the structures had piers that were skewed to primary approach flow. Scour holes near these piers consistently displayed greater depth on the side of the pier with impinging flow and deposition on the leeward side of the pier.

Missouri

Rewilding the upper Klamath River Basin: Rapid recolonization of Chinook salmon following the world’s largest dam removal

Habitat fragmentation is a global driver of decline in biodiversity, yet evaluations of ecological response following re-connection remain rare. In 2024, four hydroelectric dams on the Klamath River were removed, reconnecting over 640 km of anadromous habitat and marking the largest dam removal project completed to date. We used a stationary high-resolution multi-beam imaging SONAR to quantify adult Chinook salmon ( Oncorhynchus tshawytscha ) passage into newly reopened reaches during the first two migratory seasons following dam removal. An estimated 7,742 (95% CI: 7,702–7,778) Chinook salmon migrated upstream in 2024 and 13,310 (95% CI: 12,876–13,733) in 2025, representing 18–19% of fish returning to the Klamath River Basin. Within this period, Chinook salmon recolonized 88% of their documented historical distribution, expanding into 345 km of reconnected habitats and reaching elevations of up to 1,250 m. Recolonization required fish to swim through dam removal sites, complex hydraulic corridors and engineered fishways, demonstrating the propensity for rapid recolonization at basin scale despite these challenges. These findings provide the first quantitative assessment of salmon response to the world’s largest dam removal and demonstrate how reconnection can catalyze large-scale recolonization in historically fragmented river networks.

California, Oregon

New Hampshire and Landsat

At its widest point, a mere 80 miles separate the eastern and western borders of New Hampshire. Its northern and southern borders are just 175 miles apart. Even so, few States can boast as much rugged natural beauty per mile as the Nation’s fifth smallest. Nestled within New Hampshire are 93 State parks teeming with moose, Ursus americanus (Pallas, 1780; black bears), coyotes, beavers, river otters, and foxes. The largest section of White Mountain National Forest cuts across north-central New Hampshire, drawing visitors to its lakes, streams, mountain peaks, and hardwood forests. New Hampshire also is home to Lake Winnipesaukee, the State’s largest lake, notable for its floating post offices, the annual “ice-out” contest that sees residents vying to guess the date its surface ice dissipates, and its supporting role in films such as “On Golden Pond” and “What About Bob?” However, the scenic forests of New Hampshire face challenges in the form of invasive species such as Lymantria dispar (Linnaeus, 1758; spongy moth), Adelges piceae (balsam woolly adelgid), and Agrilus planipennis (emerald ash borer). In recent years, New Hampshire’s lakes and streams have seen more cyanobacterial blooms as well. The U.S. Geological Survey Landsat Program offers a consistent, reliable, and historically unmatched source of Earth observations that can aid in the mapping, monitoring, and management of New Hampshire’s land and water resources. Here are a few ways Landsat data have been used in the Granite State.

New Hampshire

Climate change and future water availability in the United States

The steady rise in global temperature as a result of human activity is causing changes in Earth’s water cycle. The balance of water stored within and moving between vapor, liquid, and frozen states in the water cycle is shifting, with consequences for water availability that include increases in drought, fire weather, flooding, and heavy precipitation, as well as cryosphere decline and sea-level rise. In this chapter of the U.S. Geological Survey Integrated Water Availability Assessment—2010–20, we provide an overview of climate-change observations and projections from Earth-system model simulations that relate to future water availability, from global and national climate assessments and from the published literature. Effects of climate change on primary water-cycle components are discussed in context of how global-scale hydroclimate drivers influence regional processes within the United States. Understanding the major climate drivers impacting the water cycle is crucial to predicting future changes in water availability and developing adaptation strategies to ensure human and ecosystem water supplies. First, we provide background information on the water cycle, the climate-model ensemble simulations developed to produce projections based on warming scenarios, and attribution and certainty levels. Tipping points, self-reinforcing feedbacks, cascading effects, and compound extremes are introduced. The framework of climatic impact drivers (CIDs) outlined in the Intergovernmental Panel on Climate Change Sixth Assessment Report (IPCC AR6) is used to show primary drivers of physical change to the water cycle and to understand and predict changes in future water availability. Specific climate-change related observations and projections are discussed for water cycle components of precipitation, evapotranspiration, soil moisture, streamflow, lakes and wetlands, ice and snow, and groundwater, as well as their implications for future water availability for humans and ecosystems. The chapter concludes with a synthesis discussion of three examples of complex regional-scale hydroclimate processes that influence water availability for populations in the United States, including (1) mountain and coastal precipitation, (2) aridification and drought, and (3) the influence of forest-cover change on terrestrial water-vapor recycling.

Professional Paper

Geologic map of Scoggins Dam, Henry Hagg Lake, and Scoggins Valley, Washington County, Oregon

New geologic mapping (Wells and others, 2020b) and geophysical mapping (Blakely and others, 2000; McPhee and others, 2014; Wells and others, 2020a) document kilometers of Cenozoic right-lateral offset along the Gales Creek Fault Zone, a major, northwest-striking fault zone forming the boundary between the Tualatin Valley and the Coast Range. The Bureau of Reclamation’s (Reclamation) Scoggins Dam (fig. 1), in the Coast Range foothills west of Forest Grove, Oregon, lies within the Gales Creek Fault Zone as mapped by Wells and others (2020a, 2020b; fig. 2). Active faults of the Gales Creek Fault Zone defined by paleoseismic trenching (Redwine and others, 2017, 2019b, Horst and others, 2018, 2019, 2021, and Wells and others, 2020a) are presently mapped as projecting through the existing dam. The Pacific Northwest Region of Reclamation requested assistance with geologic studies around Scoggins Dam to provide better understanding of fault locations and their activity, which are needed to design a modification of the dam (Maguire, 2019a, b). The scope of this project includes detailed geology of the existing Scoggins Dam site, Henry Hagg Lake, the reservoir behind the dam, and Scoggins Valley downstream of the existing dam, particularly around a potential new dam site, where Scoggins Creek cuts through a narrow gap formed by a resistant felsic tuff bed that crosses the valley.

Oregon

Bathymetric and velocimetric surveys at highway bridges crossing the Missouri River between Kansas City and St. Louis, Missouri, May 19–26, 2021

Bathymetric and velocimetric data were collected by the U.S. Geological Survey, in cooperation with the Missouri Department of Transportation, near nine bridges at eight highway crossings of the Missouri River between Kansas City and St. Louis, Missouri, from May 19 to 26, 2021. A multibeam echosounder mapping system was used to obtain channel-bed elevations for river reaches about 1,640 to 1,840 feet (ft) longitudinally and generally extending laterally across the active channel from bank to bank during low to moderate flood-flow conditions. These surveys provided channel geometry and hydraulic conditions at the time of the surveys and provided characteristics of scour holes that may be useful in developing or verifying predictive guidelines or equations for computing potential scour depth. These data also may be useful to the Missouri Department of Transportation as a low to moderate flood-flow assessment of the bridges for stability and integrity issues with respect to bridge scour during floods. Bathymetric data were collected around every in-channel pier. Scour holes were present at most piers for which bathymetry could be obtained, except those on banks or surrounded by riprap. Occasionally, scour holes were minor and difficult to discern from nearby dunes and ripples. All the bridge sites in this study were previously surveyed and documented in previous studies. Comparisons between bathymetric surfaces from the previous surveys and those of the current (2021) study do not indicate any consistent correlation between channel-bed elevations and streamflow conditions. The average difference between the bathymetric surfaces varied from 1.59 ft higher to 0.95 ft lower in 2021 than 2017, which corresponds to a gain of 100,200 cubic yards and a loss of 55,800 cubic yards, respectively. The average difference between the bathymetric surfaces varied from 2.74 ft higher to 3.05 ft lower in 2021 than 2013, which corresponds to a gain of 111,500 cubic yards and a loss of 169,200 cubic yards, respectively. The average difference between the bathymetric surfaces varied from 4.52 ft higher to 1.38 ft lower in 2021 than 2011, which corresponds to a gain of 221,100 cubic yards and a loss of 90,300 cubic yards, respectively. The most substantial overall net gain was 221,100 cubic yards between 2011 and 2021 at structures L0550 and A4497 at Jefferson City (site 20). The large net gain likely results from a combination of the mitigation of the scour holes near pier 4 of both bridges and the substantially lower flow in 2021 than in 2011. Alternatively, the most substantial overall net loss was 169,200 cubic yards between 2013 and 2021 at structure A6288 at Hermann (site 21), despite comparable streamflows. Pier size, nose shape, and skew to approach flow had a substantial effect on the size of the scour hole observed at a given pier. Larger and deeper scour holes were present at piers with wide or blunt noses caused by exposed footings or caissons. When a pier was skewed to primary approach flow, the scour hole was generally deeper and larger than at a similar pier without skew; furthermore, the shape of the scour hole near skewed piers in this study generally was longer and deeper on the side with impinging flow. At structure A6288 at Hermann (site 21), the scour hole near pier 5 was difficult to discern from nearby dunes and ripples, whereas the upstream edge of the footing was visible at pier 4, which likely contributes to the larger scour hole near that pier; the top of the footing may blunt the horseshoe vortex at pier 5, but the exposed front of the footing may exacerbate the vortex at pier 4.

Missouri

Global geologic map of Europa

Discovered by Galileo Galilei more than 400 years ago and imaged in detail by the Voyager 2 Galileo spacecraft, Jupiter’s icy moon Europa has been a source of intrigue. A range of science investigations indicate that it contains the key ingredients for habitability, notably energy, chemistry, and liquid water. Europa’s surface is geologically complex and, based on the dearth of impact craters, interpreted to be as young as ~60 Ma. The array of geologic features that characterize the surface include extensive ridged plains, regions of broad disruption termed chaotic terrain, long, quasi-linear ridges that span thousands of kilometers, and bands as much as 60 kilometers wide and that extend 100s of kilometers. These features, along with other geophysical measurements, indicate the presence of a global briny liquid water ocean beneath the ice shell. It was not until the arrival of the Galileo spacecraft in 1995 that the true nature and level of complexity of the surface was revealed. Although image data returned by Galileo provided insight into the structure of a variety of regions, the entire satellite has yet to be observed at a regional scale (less than 250 meters per pixel) and the detailed geologic nature of much of its surface remains a mystery. Establishing the global context of the distribution and timing of Europan geologic units forms a basis to understand regional and local scale processes, serves as a tool for the planning of future missions, and most of all is essential to gaining insight into the potential habitability of this icy world.

Scientific Investigations Map

Potentiometric surface maps and groundwater-level hydrographs for confined aquifers of the New Jersey Coastal Plain, 2018

The U.S. Geological Survey, in cooperation with the New Jersey Department of Environmental Protection (NJDEP), prepared potentiometric surface maps for 10 confined aquifers of the New Jersey Coastal Plain physiographic province based on water-level measurements collected during late 2018 and early 2019 from 951 wells in New Jersey and parts of Pennsylvania and Delaware. Maps were prepared for the confined Cohansey aquifer, Rio Grande water-bearing zone, Atlantic City 800-foot sand, Piney Point aquifer, Vincentown aquifer, Wenonah-Mount Laurel aquifer, Englishtown aquifer system, and the upper, middle, and lower aquifers of the Potomac-Raritan-Magothy aquifer system. Potentiometric surface maps indicate regional cones of depression in the following aquifers and the counties in which they are centered: Atlantic City 800-foot sand in Atlantic County, the Piney Point aquifer in Cumberland County, the Wenonah-Mount Laurel aquifer and Englishtown aquifer system in Monmouth and Ocean Counties, the Wenonah-Mount Laurel aquifer in Camden and Gloucester County, the upper aquifer of the Potomac-Raritan-Magothy aquifer system in Ocean County, and the upper, middle, and lower aquifers of the Potomac-Raritan-Magothy aquifer system in Camden County. Cones of depression with smaller areal extents were in the confined Cohansey aquifer, the Rio Grande water-bearing zone, the Atlantic City 800-foot sand centered in Cape May County, the Piney Point aquifer centered in Ocean County, the Wenonah-Mount Laurel aquifer in Salem and Burlington Counties, the Englishtown aquifer system in Camden County, the upper aquifer of the Potomac-Raritan-Magothy aquifer system in Monmouth County, and the middle aquifer of the Potomac-Raritan-Magothy aquifer system in Monmouth, Ocean, and Salem Counties. No cone of depression was interpreted in the Vincentown aquifer. Long-term hydrographs are presented for 75 wells spanning each of the 10 confined aquifers, and contain a mix of discrete water-level measurements and daily mean water levels based on continuously recorded 15-minute data. Changes of water levels during 2014–19, as indicated by the hydrographs, were compared with those of previous periods to assess any departures from historical data. During 2014–19, water levels were stable and fluctuated within similar ranges as previous periods in the following aquifers and locations: all wells in the confined Cohansey aquifer, the Rio Grande water-bearing zone, the Vincentown aquifer, the Englishtown aquifer system, the Piney Point aquifer wells in Burlington and Ocean Counties, six of eight wells in the Wenonah-Mount Laurel aquifer, all wells in the upper and lower aquifers of the Potomac-Raritan-Magothy aquifer system outside NJDEP Critical Areas, and all wells in the middle aquifer of the Potomac-Raritan-Magothy aquifer system except those within NJDEP Critical Area II. Increasing water levels in 2014–19, ongoing since historical periods, were indicated in the following aquifers and locations: Atlantic County wells in the Piney Point aquifer, all wells in the upper and lower aquifers of the Potomac-Raritan-Magothy aquifer system outside NJDEP Critical Areas, and all wells in the middle aquifer of the Potomac-Raritan-Magothy aquifer system within NJDEP Critical Area II. Water levels in the Atlantic City 800-foot sand also increased during 2014–19 in wells in Atlantic County and northern Cape May County closer to the center of the cone of depression in that aquifer, which is a response unique to this period and absent from previous periods. During 2013–19, continued decreasing water levels, ongoing since previous periods, were indicated by hydrographs of Atlantic City 800-foot sand wells in southern Cape May County, Piney Point aquifer wells in Cumberland County where the regional cone of depression is located, and two wells in the Wenonah-Mount Laurel aquifer—070478, which in 2014–19 departed from previous periods, and 330020, which continued a gradual decrease throughout its period of record.

New Jersey

Plan to coordinate post-earthquake investigations supported by the National Earthquake Hazards Reduction Program (NEHRP)

Introduction This report presents a plan supported by the National Earthquake Hazards Reduction Program (NEHRP) to coordinate domestic and international post-earthquake investigations (herein called “the Plan”). Post-earthquake scientific and engineering investigations are undertaken to capture critical information to understand the causes and impacts of the event, lessons from which can substantially improve the Nation’s resilience after future earthquakes. NEHRP is the Federal Government’s coordinated nationwide program to reduce risks to life and property from earthquakes. The Plan describes the activation and coordination of the four designated NEHRP Agencies in the Federal Government: Federal Emergency Management Agency, National Institute of Standards and Technology, U.S. National Science Foundation, and U.S. Geological Survey (USGS). The Plan also describes coordination between NEHRP Agencies and other organizations that may participate in pre-event and post-earthquake investigations, including non-NEHRP Federal agencies; State, regional, local, Tribal, and territorial agencies; domestic nongovernmental organizations; academic institutions and affiliated organizations; private companies; foreign governmental agencies and nongovernmental organizations; and international organizations. The Plan delineates the coordination of NEHRP post-earthquake scientific and engineering investigations to document the direct, indirect, and cascading physical and societal impacts from fault rupture and ground shaking hazards and from secondary hazards such as landslides, liquefaction, and tsunamis. In addition, the Plan identifies pre-event activities necessary to ensure that post-earthquake investigations are executed effectively. The USGS is the lead NEHRP Agency for activating and coordinating NEHRP post-earthquake investigations and for implementing this Plan. The USGS also leads coordination of the NEHRP Agencies in completing the pre-event activities identified in the Plan. The Plan has new information and supersedes USGS Circular 1242, “The Plan to Coordinate NEHRP Post-Earthquake Investigations,” which was published in 2003. This second-generation Plan was developed with the assistance of the Applied Technology Council (ATC) of Redwood City, California, under USGS contract 140G0121P0309, ATC-155. A 17-member Project Review Panel provided guidance on plan development, and input was solicited from subject matter experts representing key stakeholder groups and from participants at a public workshop.

Circular

Altitude of the potentiometric surface and depth to water in the Mississippi River Valley alluvial aquifer, spring 2022

Potentiometric-surface and depth-to-water maps for spring 2022 were created for the Mississippi River Valley alluvial aquifer (MRVA) using groundwater-altitude data from 1,136 wells completed in the MRVA and from the altitude of the top of the water surface in area rivers from 160 streamgages. The potentiometric-surface and depth-to-water maps for 2022 were created to support investigations to characterize the MRVA as part of the U.S. Geological Survey Water Availability and Use Science Program. Sufficient data were available to map the potentiometric surface and depth to water of the MRVA for spring 2022 for about 83 percent of the aquifer area. The potentiometric contours ranged from 0 to 340 feet (ft) above the North American Vertical Datum of 1988. The regional direction of groundwater gradient was generally to the south-southwest, except in areas of groundwater-altitude depressions, where the groundwater gradient was into the depression, and near rivers, where the groundwater gradient can be from aquifer to the river or from the river into the aquifer. There are large depressions in the potentiometric-surface map in the lower one-half of the Cache region and in much of the Grand Prairie and Delta regions. Depth to water in the MRVA, spring 2022, by well ranged from 5.00 ft above land surface to 145.66 ft below land surface.

Mississippi River Valley alluvial aquifer

Determination of per- and polyfluoroalkyl substances in water by direct injection of matrix-modified centrifuge supernatant and liquid chromatography/tandem mass spectrometry with isotope dilution

A direct-injection liquid chromatography/tandem mass spectrometry method was developed to determine 34 per- and polyfluoroalkyl substances (PFAS), including selected branched isomers, in centrifuge supernatant of matrix-modified (amended with approximately 50 percent methanol) water samples. The method has been validated in reagent water, surface water, groundwater, and wastewater effluent. Other water types (for example, drinking water, untreated wastewater, and landfill leachate) have been analyzed by the method but not systematically validated. Recovery of isotope-dilution standards, added to each sample, may be used to assess method performance in nonvalidated matrices on a sample-by-sample basis. Using this method, PFAS concentrations were determined in the range of 2–2,000 nanograms per liter in water samples. This range can be extended by diluting concentrated samples. At circumneutral pH, most compounds are present in the environment in their ionized form, and data are reported as such (for example, perfluorooctanoic acid is referred to as “perfluorooctanoate” [PFOA], perfluorooctane sulfonic acid is referred to as “perfluorooctane sulfonate” [PFOS]). Sample preparation procedures were designed without the use of filtration and with minimum sample handling steps to mitigate procedural losses of target compounds due to sorption to surfaces. Further, isotope-dilution quantification allowed for the correction of bias that may result from procedural losses, matrix-induced signal suppression or enhancement, and other factors. Validation experiments to characterize bias and variability, method detection level, and holding time were done in four distinct water matrices—reagent water, surface water, treated wastewater effluent, and groundwater—at multiple concentration levels. Mean PFAS recoveries met data quality objectives of bias and variability studies in all four validation matrices except for two compounds with low and variable recovery in the reagent water matrix only. Isotope-dilution standards, treated as surrogate compounds, were analyzed in more than 1,500 customer-submitted environmental samples with aggregate recovery of 102.5±6.5 percent (mean±standard deviation). Maximum holding times for all target compounds in the four validation matrices were 28 days for refrigerated samples and 90 days for frozen samples.

Techniques and Methods

Geophysical identification of potential groundwater – surface water interactions in the Bonita Peak Mining District, San Juan County, Colorado

This report includes the results and interpretation of geophysical investigations to aid in the location, identification, and conceptual site model (CSM) development of the near surface groundwater and surface water hydrology of the Upper Cement Creek and California Gulch of the Animas Rivers within the Bonita Peak Mining Districts (BPMD). Mine influenced water (MIW) from numerous mine adits within the BPMD has altered the geochemistry and ecology of these ground and surface waters. Many studies of varying scales have investigated the geology, hydrogeology, geochemistry, and aqueous chemistry of these adits and impacted areas. Field investigations using non-invasive geophysical methods were conducted in the summer of 2019 to aid in identifying potential groundwater inputs to these streams and further the CSM of the near surface within and along these streams. Two field deployments were completed in August and September of 2019. In August 2019, Phase 1 was completed using multifrequency electromagnetic induction (EMI) and thermal infrared imaging (TIR) to measure the bulk electrical conductivity, magnetic susceptibility, and the temperature response at known and suspected groundwater seeps into the Upper Cement Creek and California Gulch stretches of these alpine streams. Results and interpretation from Phase 1 were used to co-locate stretches of these streams for further detailed fiber optic distributed temperature system (FODTS) monitoring of the temperature dynamics for submeter identification of potential groundwater seeps during Phase 2. Phase 2 included additional EMI measurements, FODTS measurements, and repeated some Phase 1 measurements for improved data quality. All these data were used for the combined qualitative interpretation identifying seep locations. The data from both Phase 1 and Phase 2 are available in the USGS data release at: https://www.sciencebase.gov/catalog/item/5f18a30082cef313ed84879b. These data include detailed near-surface electrical conductivity and magnetic susceptibility mapping, discrete TIR imaging, and several kilometers of FODTS data. These data were georeferenced and used to identify potential groundwater seeps into these streams. The results indicate a complex near surface geophysical and hydrogeologic regime. Seeps are classified and spatially correlated to the EMI, TIR, and geochemical sampling results. This report is a thorough compilation of the geophysical field efforts in the summer of 2019, the data processing, and interpretations of these data. These data alone should not be used for final interpretations of the complex hydrogeology within these alpine systems. Rather, these data should be combined with other known data from various investigations to further refine the CSM and guide decisions for future investigations, watershed protective measures, and other site decisions. Geographical information system (GIS) layers of these results are available to project managers for such integration and combined interpretations.

Colorado

Eruptive history within the vicinity of Al Madīnah in northern Harrat Rahat, Kingdom of Saudi Arabia

The northernmost part of the Harrat Rahat volcanic field contains early Pleistocene to Holocene mafic eruptive products within the vicinity of the city of Al Madīnah, Kingdom of Saudi Arabia. A detailed geologic investigation into the eruptive history of a 570 square kilometer (km 2 ) area covering Al Madīnah and the surrounding area has yielded 33 mapped Quaternary volcanic units consisting of lava flows, scoria cones, and shield volcanoes. These eruptive products consist of continental, intraplate alkalic and minor transitional basalts, hawaiites, and a single mugearite that were emplaced from at least 1,014±14 thousand years ago (ka) to a single Holocene eruption in 1256 C.E. Lava flows are generally 10 to 15 kilometers (km) long (but can reach 23 km long), 1 to 3 km wide, and at least 10 meters thick. Most of the mapped units erupted episodically between 400 and 340 ka and 180 and 100 ka. Despite small individual volumes (less than 1 cubic kilometers dense rock equivalent), each unit represents eruption of a distinct magma batch that was strongly influenced by clinopyroxene, olivine, and plagioclase fractionation. Some of these units are interpreted to have undergone magma mixing pre- and (or) syneruptively. Integrating eruption ages, geochemistry, and paleomagnetic data yields evidence that some eruptions were temporally and (or) spatially clustered. Aligned scoria cones and elongate vent edifices were constructed atop fissure vent systems that reflect the local stress field, which controls dike ascent through the middle and upper crust.

Professional Paper

Salinas Valley integrated hydrologic and reservoir operations models, Monterey and San Luis Obispo Counties, California

The area surrounding the Salinas Valley groundwater basin in Monterey and San Luis Obispo Counties of California is a highly productive agricultural area, contributes substantially to the local economy, and provides a substantial portion of vegetables and other agricultural commodities to the Nation. This region of California provides about half of the Nation’s lettuce, celery, broccoli, and spinach each year. Thus, this agricultural area provides substantial volumes of agricultural products not just for California but for the United States. Changes in population and increased agricultural development, which includes a shift toward more water-intensive crops, and climate variability, have put increasing demand on both surface-water and groundwater resources in the valley. This situation has resulted in water management challenges in the Salinas Valley that generally relate to the distribution of the water supply throughout the basin. Where and when the water is present in the surface and subsurface does not coincide with where and when the water is needed. Historically, to deal with the distribution issue, water has been used conjunctively in the valley. Conjunctive use is a water management strategy that coordinates surface-water and groundwater use to maximize water availability. Groundwater is used throughout the Salinas Valley to meet water demands when surface-water supplies are insufficient. The availability of surface water is constrained by climate. Precipitation and streamflow vary seasonally and year to year. Although there are two reservoirs in the Salinas Valley to capture and store water during wet periods, the only conveyance of reservoir water to coastal agricultural areas is the Salinas River. Increasing demand for groundwater and surface-water resources throughout the Salinas Valley has resulted in undesirable effects from unsustainable water use, such as surface-water depletion, groundwater-level declines, storage depletion in the principal aquifers, and seawater intrusion. To address these escalating issues, local communities, water management agencies, and groundwater sustainability agencies are evaluating how to sustainably manage both their surface-water and groundwater resources. To meet water demands and reduce the undesirable effects of unsustainable water use, continued conjunctive management of surface water and groundwater would ideally incorporate strategies to deal with increases in demand and climate variability. To evaluate the challenging water management issues in the Salinas Valley, the U.S. Geological Survey, Monterey County Water Resources Agency, and the Salinas Valley Basin Groundwater Sustainability Agency developed a comprehensive suite of models that represent the Salinas Valley hydrogeologic system called the Salinas Valley System Model. The geologic framework is known as the Salinas Valley Geologic Framework and was developed to characterize the subsurface using various topographic and geologic data sources, including information on hydrogeologic units, their surfaces and extents, geologic structures, lithology, and elevations from borehole data and cross sections, as well as details on faults and existing models. The surface-water model is called the Salinas Valley Watershed Model and simulates the Salinas River watershed. Monthly surface-water inflows into the integrated hydrologic model domain were simulated using the Salinas Valley Watershed Model. The historical model uses historical climate data, water and land use data, and reservoir releases to simulate agricultural operations, including landscape water demands, diversions, and reclaimed wastewater. The operational model adds an embedded reservoir operations framework to the simulation of the historical model that allows specified operational rules to simulate reservoir releases and changes in reservoir storage. The operational model assumes current reservoir operations and constant land use, which differs from historical conditions. Thus, the operational model is a hypothetical baseline model that can be used by local water managers to evaluate and quantify potential benefits of water supply projects. Together, the geologic framework, watershed, historical, and operational models form a tool that can be used to simulate irrigated agriculture and associated reservoir operations of the integrated hydrologic system of the Salinas Valley.

California

Assessment and characterization of ephemeral stream channel stability and mechanisms affecting erosion in Grand Valley, western Colorado, 2018–21

The Grand Valley in western Colorado is in the semiarid Southwest United States. The north side of the Grand Valley has many ungaged ephemeral streams, which are of particular interest because (1) the underlying bedrock geology, Late Cretaceous Mancos Shale, is a sedimentary rock deposit identified as a major salinity contributor to the Colorado River and (2) despite infrequent streamflows of short duration, monsoon-derived floods in these ephemeral streams can carry substantial amounts of sediment downstream, affecting upstream and downstream banks and channel cross sections. The study area is of interest, because salinity, or the total dissolved solids concentration, in the Colorado River causes an estimated $300 million to $400 million per year in economic damages in the United States, and it is estimated 62 percent of the Upper Colorado River Basin’s total dissolved solid loads originate from geologic sources. In an effort to minimize salt contributions to the Colorado River from public lands administered by the Bureau of Land Management, a comprehensive salinity control approach is typically used to reduce nonpoint sources of salinity through land management techniques and practices. In 2018, the U.S. Geological Survey, in cooperation with the Bureau of Land Management, began an assessment of ephemeral streams located on the north side of the Grand Valley, western Colorado, to characterize stream channel stability and identify mechanisms affecting erosion. The U.S. Geological Survey developed a method for automatically extracting channel cross-section geometry from existing remotely sensed terrain models. Based on estimated flood stage and surrogate streamflows, hydraulic characteristics were calculated. Furthermore, the channel geometries and hydraulic characteristics were used to estimate channel stability using a statistical model. Cross-section stabilities were determined from a stream channel stability assessment for a subset of 1,406 visited (field observed) locations out of 13,415 cross sections, which were delineated from remotely sensed terrain models. The application of Manning’s resistance equation in combination with multiple logistic regression models demonstrated channel stability can be estimated with a 0.845 goodness of fit for a validation dataset when using a combination of drainage area, width-to-depth ratio, sinuosity, and shear stress as the explanatory variables. Stream channel stability was extrapolated for 13,415 unvisited (not field observed) cross sections using the multiple logistic regression model and defined explanatory variables. Mapping of the ephemeral streams and their associated stabilities may be used by the Bureau of Land Management to prioritize areas for remediation or changes in management strategies to reduce sediment and salinity loading to the Colorado River. The study found channel stability within the ephemeral streams to be spatially variable, longitudinally discontinuous, and dictated by changes in channel bed slope. The stable ephemeral streams were relatively wide and shallow and often had smaller drainage areas with less potential for producing shear stresses capable of overcoming channel adhesion. A change in channel bed slope can provide the means necessary to generate shear stresses appropriate to initiate erosion and a subsequent stability transition to incising channels. Channel widening happens when either or both banks of an incising channel reach a critical height for mass wasting, or when channel curvature causes higher sidewall stress. Regardless, widening channels can promote increases in sinuosity and subsequently reduce steep channel bed slopes. Consequently, stable and widening channels can have comparable bed slopes, making channel bed slope a poor explanatory variable to predict channel stability overall, despite its function to initiate channel instability. The results were based on a surrogate 0.10 annual exceedance probability (AEP; return period equal to the 10-year flood) interval streamflow, although it was recognized fluctuations in streamflow would also affect channel stability. Past and current changes within the study area affect streamflow; therefore, mechanisms affecting erosion include land use disturbances, soil compaction, loss of vegetation cover, drought, less frequent and more extreme precipitation, and fires—which all intensify the potential runoff and erosion within the study area.

Colorado