USGS ScienceSearch

SEARCH · USGS Science

Results for “General Information Product”

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.

831 records · Page 36Linked to original sources

Global variability of the composition and temperature at the 410-km discontinuity from receiver function analysis of dense arrays

Seismic boundaries caused by phase transitions between olivine polymorphs in Earth's mantle provide thermal and compositional markers that inform mantle dynamics. Seismic studies of the mantle transition zone often use either global averaging with sparse arrays or regional sampling from a single dense array. The intermediate approach of this study utilizes many densely spaced seismic arrays distributed around the globe. We systematically compute teleseismic P-to-S receiver functions for each seismic array and invert for the 1-D seismic velocity structure of the mantle transition zone beneath each array to facilitate a comparison between densely sampled regions. We stack 3,600 receiver functions on average at 67 arrays in total. The stack is used in a probabilistic inversion to estimate the mantle transition zone interface depths and velocities beneath each array. We focus on the 410-km discontinuity (410) because it is a prominent seismic interface that is clearly linked to a single mineral phase transition between olivine and wadsleyite. The depths and velocity contrasts of the 410 are mapped to temperatures and compositions using mineral physics constraints. The depth of the 410 ranges from ∼405–440 km, which is consistent with a ∼360 K temperature range in a dry mantle and a ∼260 K temperature range in a wet mantle (2 wt. % water). The Vs contrast across the 410 ranges from ∼2.5–8 %, which is consistent with ∼20–70 vol. % olivine composition in a dry mantle and ∼25–80 vol. % in a wet mantle. The bulk composition of the upper mantle near the 410-km discontinuity is typically considered to be well-mixed because there is no thermodynamic impediment to convection at the olivine to wadsleyite phase transition. However, the wide range of inferred olivine content from our study suggests that there are large lateral variations in the bulk composition of the upper mantle near the 410-km discontinuity.

Earth and Planetary Science Letters

Sources and risk factors for nitrate, pathogens, and fecal contamination of private wells in rural southwestern Wisconsin, USA

Household well water can be degraded by contaminants from the land's surface, but private well owners lack means to protect the source water from neighboring disturbances. Rural residents of southwestern Wisconsin, USA, rely on private well water, and the combination of land use and fractured carbonate bedrock makes groundwater vulnerable to contamination. To identify the extent, sources, and risk factors of private well contamination, randomly selected wells sampled during two-day periods in fall (n = 301) and spring (n = 529) were analyzed for nitrate and indicator bacteria, and a subset (n = 138) was sampled across four seasonal events for analysis of pathogens and microbial source tracking markers by quantitative polymerase chain reaction. Risk factors representing land use, hydrology, geology, and well construction were analyzed for associations with contamination in multivariable models. The importance of risk factors varied by contaminant, illustrating the multifaceted nature of rural groundwater quality. Nitrate contamination was associated with agricultural land use, and wells with casings that extended below a shale aquitard accessed less contaminated water than those drawing water from above it. Human fecal microbes were detected in 64 wells (46%), and rainfall was the key risk factor for contamination, indicating that wastewater from septic systems was available to contaminate wells when transport conditions were favorable. Manure microbes from cattle/ruminants and pigs were detected in 33 and 13 wells, respectively, and concentrations increased with the hectarage of cultivated land near wells. Pathogen genes for viruses, bacteria, and protozoa were detected in 66 wells (48%), including more detections of zoonotic than human-specific pathogens, and human Bacteroides , an indicator of wastewater, was an equivocal predictor of pathogen presence in private wells. Characterizing important elements of the setting, like geology, and identifying sources and risk factors for contaminants can inform landscape-level policies to protect groundwater quality.

Wisconsin

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

3D Converted wave reverse time migration imaging

We describe a newly developed method for recovering high-resolution images of seismic discontinuities, such as subducting slabs, in 3D. Our method makes use of converted S → P or P → S waves observed by dense arrays of seismometers to infer the locations and relative strengths of seismic discontinuities at depth in a target region. Observed direct and converted waves are backpropagated to their times of origin. The time-reversed wavefield is then separated into its constituent P and S components via the Helmholtz decomposition, and those separated wavefields are used to compute imaging functions that characterize the locations and relative strengths of seismic discontinuities. Imaging functions may be designed to use either S → P or P → S waves, so that users can target those arrivals expected to be most dominant in a given dataset. We have previously demonstrated the efficacy of our method in two dimensions, and we now present a 3D implementation of our technique which addresses the significant computational challenges posed by the size of volumetric wavefield data in three dimensions. Through a series of synthetic examples, we demonstrate that our method is capable of recovering the fine scale structure of a subducting slab given realistic station coverage and earthquake sources. We investigate optimal seismic station geometries for our technique and explore image interpretability in regions with poor data coverage. We find that linear station geometries yield more optimal, interpretable imaging functions than collections of small arrays can. We also show that our method can successfully recover bothS → P or P → S images when realistic shear earthquake sources are used, and we explore the additional computational challenges presented by the high frequency content of S waves. Our results demonstrate the potential for our technique to recover high-resolution information about subducting slabs in real-world regions, given that relatively sparse seismic arrays with only approximately 100 stations are capable of recovering interpretable imaging functions from just a few realistic earthquake sources for multiple discontinuities at significant depth in an area of approximately 400~sq~km.

Seismica

Cambrian and Ordovician stratigraphy, conodont biostratigraphy, and microfacies analysis to support 1:24,000-scale geologic mapping of the southern Lake Champlain valley, New York and Vermont

Introduction Geologic mapping in the southern Lake Champlain valley of New York and Vermont (fig. 1) has required evaluation of stratigraphic nomenclature used since the early 1900s. The paleogeography of the Cambrian and Ordovician adjacent to the Adirondack Highlands, the high-relief mountains of Proterozoic igneous and metamorphic rocks of northern New York, has impacted the lithostratigraphic succession, facies changes, and several hiatuses that occur within the package of clastic and carbonate rocks that demonstrates overall deepening consistent from the southern to northern Appalachians. The stratigraphic nomenclature used for the geologic mapping of the Cambrian and Ordovician strata of the study area is derived from evaluation of names put forth by previous workers from areas of northern New York and western Vermont. Descriptions of type localities or type sections and other criteria set forth by the North American Stratigraphic Code (North American Commission on Stratigraphic Nomenclature [NACSN], 2021) were consulted. The North American Stratigraphic Code states that lithostratigraphic units, such as the fundamental unit, the formation, be defined by lithic character and be mappable (NACSN, 2021, Article 24). However, some early workers (Walcott, 1912; Kay, 1937; Fisher and Hanson, 1951; Oxley and Kay, 1959) have defined stratigraphic units based on fossil content, which may or may not be facies dependent, and hence may not be mappable as lithostratigraphic units. Although preservation of older names takes priority in stratigraphic nomenclature (NACSN, 2021, Article 7), some do not have adequately defined stratotypes or are located a distance away with different lithologies than that in the study area. Biostratigraphy can be a helpful tool for understanding correlations and facies changes. Herein, conodont biostratigraphy is utilized to date and correlate separate lithologic units. Further, microfacies analyses of units through inspection of petrographic thin sections provide paleogeographic information that helps in understanding the genesis of geologic units.

New York, Vermont

Vulnerability assessment of groundwater influenced ecosystems in the Northeastern United States

Groundwater-influenced ecosystems (GIEs) are increasingly vulnerable due to groundwater extraction, land-use practices, and climate change. These ecosystems receive groundwater inflow as a portion of their baseflow or water budget, which can maintain water levels, water temperature, and chemistry necessary to sustain the biodiversity that they support. In some systems (e.g., springs, seeps, fens), this connection with groundwater is central to the system’s integrity and persistence. Groundwater management decisions for human use often do not consider the ecological effects of those actions on GIEs. This disparity can be attributed, in part, to a lack of information regarding the physical relationships these systems have with the surrounding landscape and climate, which may influence the environmental conditions and associated biodiversity. We estimate the vulnerability of areas predicted to be highly suitable for the presence of GIEs based on watershed (U.S. Geological Survey Hydrologic Unit Code 12 watersheds: 24–100 km 2 ) and pixel (30 m × 30 m pixels) resolution in the Atlantic Highlands and Mixed Wood Plains EPA Level II Ecoregions in the northeastern United States. We represent vulnerability with variables describing adaptive capacity (topographic wetness index, hydric soil, physiographic diversity), exposure (climatic niche), and sensitivity (aquatic barriers, proportion urbanized or agriculture). Vulnerability scores indicate that ~26% of GIEs were within 30 m of areas with moderate vulnerability. Within these GIEs, climate exposure is an important contributor to vulnerability of 40% of the areas, followed by land use (19%, agriculture or urbanized). There are few areas predicted to be suitable for GIEs that are also predicted to be highly vulnerable, and of those, climate exposure is the most important contributor to their vulnerability. Persistence of GIEs in the northeastern United States may be challenged as changes in the amount and timing of precipitation and increasing air temperatures attributed to climate change affect the groundwater that sustains these systems.

Connecticut, Delaware, District of Columbia, Maine

Haunted Summerville: Ghostly lights or earthquake lights?

Among the colorful local lore in the Charleston, South Carolina, area, are a number of ghost stories, shared not only over campfires but also in published books. Among the most well-known of the stories is the tale of the Summerville Light. Local lore holds that a strange light sometimes seen in a remote area is a lantern carried by the ghost of a local woman who once waited hours for her husband, who turned out to have been decapitated earlier that day in a train accident (DePoppe, 2023). Extant sources suggest the ghost stories began to circulate in the 1950s to 1960s. So pervasive was the lore that (Old) Sheep Island Road became known among local residents as Light Road, with a local stretch of road known today as Old Light Road. Reviewing the location where the lights appear as well as the nature of accounts, I suggest that many if not all of the anecdotal observations can be most readily attributed to natural phenomena, including earthquake lights from earthquakes that were too small to be felt. Accounts of lights near Summerville cluster in proximity to the generally accepted epicenter of the 1886 Charleston, South Carolina, earthquake, where foreshocks to the 1886 mainshock were apparently concentrated, and within a few kilometers of three M3.5 – 4.4 earthquakes in 1959 and 1960.

South Carolina

Profiles of gamma-ray and magnetic data for aerial surveys over parts of the Western United States from longitude 108 to 126 degrees W. and from latitude 34 to 49 degrees N.

This CD-ROM contains images generated from geophysical data, software for displaying and analyzing the images and software for displaying and examining profile data from aerial surveys flown as part of the National Uranium Resource Evaluation (NURE) Program of the U.S. Department of Energy. The images included are of gamma-ray data (uranium, thorium, and potassium channels), Bouguer gravity data, isostatic residual gravity data, aeromagnetic anomalies, topography, and topography with bathymetry. This publication contains image data for the conterminous United States and profile data for the conterminous United States within the area longitude 108 to 126 degrees W. and latitude 34 to 49 degrees N. The profile data include apparent surface concentrations of potassium, uranium, and thorium, the residual magnetic field, and the height above the ground. The images on this CD-ROM include graytone and color images of each data set, color shaded-relief images of the potential-field and topographic data, and color composite images of the gamma-ray data. The image display and analysis software can register images with geographic and geologic overlays. The profile display software permits the user to view the profiles as well as obtain data listings and export ASCII versions of data for selected flight lines.

Data Series

A great tsunami earthquake component of the 1957 Aleutian Islands earthquake

The great 1957 Aleutian Islands earthquake ruptured ∼1200 km of the plate boundary along the Aleutian subduction zone and produced a destructive tsunami across Hawaiʻi. Early seismic and tsunami analyses indicated that large megathrust fault slip was concentrated in the western Aleutian Islands, but tsunami waves generated by slip in the west cannot explain the large observed runup in Hawaiʻi far to the southeast. Recently mapped 1957 geologic deposits on eastern Aleutian Islands suggest occurrence of very large nearby slip. Jointly modeling tsunami runup along the eastern Aleutian and Hawaiian Islands together with tide gauge recordings across the Pacific resolves 12-26 m shallow slip along 600 km of the eastern Aleutian Islands in addition to modest, deeper western slip inferred from seismic records. The eastern near-trench slip results in an M W 8.3-8.6 tsunami earthquake component of the M W 8.6-8.8 rupture, comparable in size to the adjacent 1946 Aleutian tsunami earthquake to the east. The reexamination of the 1957 rupture confirms the tsunami hazards posed by the eastern Aleutian subduction zone to Hawaiʻi and lays the groundwork for investigation of large prehistoric earthquakes through modeling tsunami runup inferred from stratigraphic observations to constrain their rupture processes.

Alaska

Wind River subbasin restoration: Annual report of U.S. Geological Survey activities January 2024 through December 2024

We sampled juvenile wild Oncorhynchus mykiss (Steelhead Trout) in headwater streams of the Wind River, Washington, to characterize population attributes and investigate life-history metrics, particularly migratory patterns, and early life-stage survival. Additionally, this sampling provides an index of abundance of non-native Brook Trout Salvelinus confluentus . We used passive integrated transponder (PIT) tagging and a series of instream PIT-tag interrogation systems (PTISs) to track juveniles and adults. The Wind River subbasin is considered a wild Steelhead refuge by Washington Department of Fish and Wildlife (WDFW). No hatchery Steelhead Trout have been released in the Wind River subbasin since 1997, and hatchery adults are estimated at less than one percent of spawners in most years. Over 25 years of Steelhead Trout status and trend monitoring and research in the subbasin is contributing to understanding of population response to numerous restoration actions in the subbasin, including removal of Hemlock Dam from Trout Creek in 2009, which had an inadequate adult ladder and contributed to increased water temperatures, to near lethal levels for salmonids. Data from our study, and companion work by WDFW, are contributing to the Columbia Basin Fish and Wildlife Program (2009) Research, Monitoring, and Evaluation (RM&E) Strategy of Fish Population Status Monitoring. Specifically, this work addresses the sub-strategies of 1) Assessing the Status and Trends of Diversity of Natural Origin Fish Populations and Uncertainties Research regarding differing life histories of a wild Steelhead Trout population, 2) Assessing the Status and Trend of Adult Natural Origin Fish Populations, and 3) Monitoring and Evaluating the Effectiveness of Tributary Habitat Actions Relative to Environmental, Physical, or Biological Performance Objectives. During summer and fall 2024, we PIT-tagged 2,075 Steelhead Trout parr (age-0 and age-1), in the Trout Creek and upper Wind River watersheds. Age-0 and age-1 parr were at high abundance due to strong spawner numbers in 2023 and 2024 (Estimate of 932 in 2023; and a moderate number of spawners in 2024 though an estimate is not available at this time; 24-year median = 619; Charlie Cochran, WDFW Fish Biologist, personal commun., 2025). An additional 100 parr (age-0 and age-2 or older) were tagged to provide fish for estimating detection efficiencies at PTISs. Steelhead Trout parr were recaptured and detected through repeat headwater sampling, smolt trapping, instream PTISs and Columbia River mainstem PIT-tag detection. We maintained and upgraded six instream PTISs to detect PIT-tagged Steelhead Trout parr, smolts, and adults, providing data for population assessments, and life-cycle research. Detection data from PIT-tagged adult Steelhead Trout at our four primary PTISs allow assessment of adult escapement to tributary watersheds within the Wind River subbasin. Most adult Steelhead detected were tagged by WDFW at an adult trap in the lower Wind River and some as smolts at screw traps. Adult Steelhead Trout detection efficiency estimates at our primary PTIS in Trout Creek have been greater than 99 percent during eight of the past ten years and have exceeded 97% at our primary PTIS in the Wind River during nine of the past ten years. Adult escapement estimates to tributary watersheds are helping evaluate the efficacy of the 2009 removal of Hemlock Dam from rkm 2.0 of Trout Creek. Detections at the instream PTISs have shown trends of age-0 and age-1 Steelhead Trout parr emigration from natal areas during summer and fall, in addition to the expected movement of parr and smolts in spring. Our data suggest that most fish from a cohort that migrate downstream will do so at age-1 for additional rearing in the subbasin. We have estimated that from 15% to 56% of parr tagged as age-0 fish in headwater areas make downstream migrations at age-1 for additional rearing (median = 29%). We have estimated that from 1% to 27% of Steelhead Trout parr, tagged as age-1 fish, make downstream migrations during fall (median = 3.5%): this is especially pronounced in the upper Wind River portion of the watershed (median 10%). These findings raise questions about parr rearing habitat and whether migrations are density- or habitat-quality driven, and answers to such questions are long-term goals of this study as active and passive habitat restoration occurs. Repeat sampling at sites in the subbasin within and between years has enabled assessment of juvenile Steelhead Trout growth patterns. Growth rates (relative change in weight) of age-0 PIT-tagged parr during summer have been similar across the subbasin, though slightly lower in the Trout Creek watershed (Trout Creek sites, 0.57 - 1.25%/day, Wind River sites = 0.94 – 2.48%/day). Summer growth rates have been lower for age-1 parr in the Trout Creek watershed than the upper Wind River watershed (Trout Creek sites, = 0.03 - 0.14%/day, Wind River sites = 0.05 – 0.19%/day). Yearly relative growth was similar across the subbasin for both age-0 and age-1 tagged parr. During 2024 range of mean daily relative change in weight from the previous year for age-0 tagged fish was: Trout Creek sites = 0.96 – 1.41%/day, Wind River sites = 0.80 – 1.30 %/day; and for age1 tagged fish was: Trout Creek sites = 0.37 - 0.50 %/day, Wind River sites = no age-1 tagged recaptures. Non-native Brook Trout are present in the subbasin, chiefly the Trout Creek watershed. Mean percent-of-catch that is Brook Trout, at four sample sites in Trout Creek, has declined from the period 1998 – 2003 to the period 2011 – 2024 (range of decrease = 19 - 82%). Percent-of-catch (of 64 samples, 52 < 10%; 3 > 20%) and number of Brook Trout at the Trout Creek sites from 2011 through 2024, though variable, have generally declined. Habitat restoration evaluation efforts are critical to ensure efficient use of money and resources but require long-term monitoring of population metrics. Detailed viable salmonid population and life-history data, such as that provided by PIT-tagging and instream PTIS networks inform fisheries policy and management and enable assessment of long-term effects of habitat restoration actions such as the removal of Hemlock Dam on Trout Creek and proposed major instream habitat restoration in the upper Wind River.

Washington

Unscrambling the Proterozoic supercontinent record of northeastern Washington State, USA

The time interval from Supercontinent Nuna assembly in the late Paleoproterozoic to Supercontinent Rodinia breakup in the Neoproterozoic is considered by some geologists to comprise the “Boring Billion,” an interval possibly marked by a slowdown in plate tectonic processes. In northeastern Washington State, USA, similar to much of western Laurentia, early workers generally thought the tectonostratigraphic framework of this interval of geologic time consisted of two major sequences, the (ca. 1480–1380 Ma) Mesoproterozoic Belt Supergroup and unconformably overlying (<720 Ma) Neoproterozoic Windermere Supergroup. However, recent research indicates that strata considered by early workers as Belt Supergroup equivalents are actually younger, and a post-Belt, pre-Windermere record is present within the <1360 Ma Deer Trail Group and <760 Ma Buffalo Hump Formation. Thus, the northeastern Washington region perhaps comprises the most complete stratigraphic record of the “Boring Billion” time interval in the northwestern United States and holds important insights into global Proterozoic supercontinent tectonic processes. In light of these exciting developments, this field guide will address the early historic economic geology and original mapping of these Proterozoic sequences in the northeastern Washington region, and from that foundation explore more recent isotopic provenance data and their regional to global context. Finally, the guide will end with a discussion of remaining questions with a goal of stimulating interest in these relatively understudied, yet important, rocks.

Washington

PEST++IES how many iterations & realizations, finding the point of diminishing returns

PEST++IES (White 2018; White et al. 2020) is widely used in the groundwater modeling community for its ability to perform computationally efficient history matching and uncertainty analysis in a highly parameterized context. One primary advantage of using an iterative ensemble smoother is that the number of model runs required per iteration depends on the number of realizations in an ensemble, not the number of parameters in each realization. However, this raises the question: what is the optimal number of realizations and iterations to use for any one model before the point of diminishing returns? Using a modified version of the Freyberg model (Freyberg 1988; Hunt et al. 2020), different parameter and observation scenarios were evaluated for four iterations and ensembles of 10, 25, 50, 100, 250, 500, 1000, and 2000 realizations. To match observations, PEST++IES altered hydraulic conductivity ( k ), both globally across the model and locally at three different pilot point densities, as well as global recharge (via a single multiplier), global river conductance, and individual well flow rates. Risk-based well capture zone results (Fienen et al. 2022a) and estimated hydraulic conductivity fields from each scenario were quantitatively and qualitatively compared against the “truth” model and its outputs. Across the cases examined, ensemble sizes of 100 to 250 realizations and two PEST++IES iterations were generally sufficient to achieve good results.

Groundwater

Concentrations and distribution of manmade organic compounds in the Lake Tahoe Basin, Nevada and California, 1997-99

The U.S. Geological Survey, in cooperation with the Tahoe Regional Planning Agency and the Lahontan Regional Water-Quality Control Board, sampled Lake Tahoe, major tributary streams to Lake Tahoe, and several other lakes in the Lake Tahoe Basin for manmade organic compounds during 1997-99. Gasoline components were found in all samples collected from Lake Tahoe during the summer boating season. Methyl tert -butyl ether (MTBE), benzene, toluene, ethylbenzene, and xylenes (BTEX) were the commonly detected compounds in these samples. Most samples from tributary streams and lakes with no motorized boating had no detectable concentrations of gasoline components. Motorized boating activity appears to be directly linked in space and time to the occurrence of these gasoline components. Other sources of gasoline components to Lake Tahoe, such as the atmosphere, surface runoff, and subsurface flow, are minor compared to the input by motorized boating. Water sampled from Lake Tahoe during mid-winter, when motorized boating activity is low, had no MTBE and only one sample had any detectable BTEX compounds. Soluble pesticides rarely were detected in water samples from the Lake Tahoe Basin. The only detectable concentrations of these compounds were in samples from Blackwood and Taylor Creeks collected during spring runoff. Concentrations found in these samples were low, in the 1 to 4 nanograms per liter range. Organochlorine compounds were detected in samples collected from semipermeable membrane devices (SPMD's) collected from Lake Tahoe, tributary streams, and Upper Angora Lake. In Lake Tahoe, SPMD samples collected offshore from urbanized areas contained the largest number and highest concentrations of organochlorine compounds. The most commonly detected organochlorine compounds were cis- and trans-chlordane, p , p' -DDE, and hexachlorobenzene. In tributary streams, SPMD samples collected during spring runoff generally had higher combined concentrations of organochlorine compounds than those collected during baseflow conditions. Upper Angora Lake had the fewest number of organochlorine compounds detected of all lake samples. Dioxins and furans were not detected in SPMD samples from two sites in Lake Tahoe or from two tributary streams. The number of polycyclic aromatic hydrocarbon (PAH) compounds and their combined concentrations generally were higher in samples from Lake Tahoe than those from tributary streams. Areas of high-motorized boating activity at Lake Tahoe had the largest number and highest concentrations of PAH's. PAH compounds were detected in samples from SPMD's in four of six tributary streams during spring runoff, all tributary streams during baseflow conditions, and at all lake sites. The most commonly detected PAH's in tributary streams during spring runoff were phenanthrene, fluoranthene, pyrene, and chrysene, and during baseflow conditions were phenanthrene, 1-methylphenanthrene, diethylnaphthalene, and pyrene. Upper Truckee River, which has an urban area in its drainage basin, had the largest number and highest combined concentration of PAH's of all stream samples. Bottom-sediment from Lake Tahoe had detectable concentrations of p-cresol, a phenol, in all but one sample. A sample collected near Chambers Lodge contained phenol at an estimated concentration of 4 micrograms per kilogram (µg/kg). Bottom-sediment samples from tributary streams had no detectable concentrations of organochlorine or PAH compounds. Several compounds were detected in bottom sediment from Upper Angora Lake at high concentrations. These compounds and their concentrations were p , p' -DDD (10 µg/kg), p , p' -DDE (7.4 µg/kg), 2,6-dimethylnaphthalene (estimated at 190 µg/kg), pentachlorophenol (3,000 µg/kg), and p-cresol (4,400 µg/kg).

California, Nevada

Lithium from magma to mine in an early Yellowstone hotspot caldera

Renewable energy technologies rely on the extraction of metals not historically in high demand, such as lithium (Li), for which ore deposit models are incompletely understood. One of the world’s largest Li deposits is hosted in lake sediments of the 16.4 Ma McDermitt caldera, which formed during the early stages of Yellowstone hotspot volcanism in the western United States. Eruptive and posteruptive mobility of Li are major challenges in elucidating deposit formation. Melt inclusions preserved in quartz crystals provide a means to assess pre-eruptive magmatic Li contents. Concentrations of Li determined by ion microprobe for melt inclusions in a McDermitt rhyolite lava are 400−1350 ppm, compared to 20−70 ppm Li in matrix rhyolite glasses. Synthesis with melt inclusion data for eight additional calderas demonstrates a recurrence of Li-rich rhyolitic magmas (200−2000 ppm Li) in the western part of the Yellowstone hotspot track. However, unlike the multicyclic caldera complexes with overlapping fault networks that may have compromised Li retention, the McDermitt caldera remained a closed hydrologic system throughout its evolution. Modeling indicates 100 km 3 of resurgent magma could yield 25−150 Mt Li in a magmatic fluid and supports accumulation of Li-rich magmatic fluid in a closed intracaldera lake, followed by evaporative concentration and sequestration of Li within clay minerals to generate the McDermitt deposit.

Idaho, Nevada, Oregon, Wyoming

Development of liquefaction-targeted design basis in U.S. seismic provisions

Seismic design criteria for new buildings in the United States have historically focused on life safety performance objectives through minimizing the potential for structural collapse. Development of design criteria to meet this performance objective has evolved over time, leading to the current, risk-targeted maximum considered earthquake (MCE R ) design basis. Corresponding MCE R seismic design loads incorporate the full range of probabilistic ground motion hazard information available from U.S. Geological Survey National Seismic Hazard Models, as well as a representation of the uncertainty in structural collapse capacity. Despite these considerable advances in structural design criteria, guidelines for geotechnical hazards such as liquefaction-induced damage, which continue to be based on uniform ground-shaking hazard (i.e., MCE ground motions), are not well-connected with risk-based, collapse prevention performance objectives. This study, undertaken as part of Building Seismic Safety Council efforts to update recommended seismic provisions through the National Earthquake Hazards Reduction Program, presents probabilistic, liquefaction-targeted design objectives and corresponding ground motion parameters, with a focus on improving consistency in first-level screening criteria for liquefaction hazard assessment. This study identifies a potential maximum acceptable annualized probability of liquefaction triggering ( P L,max ) of 2.5% in 50 years, above which liquefaction consequences and potential mitigation measures would be considered. The P L,max -based design basis is calibrated to effective probabilistic liquefaction hazard levels obtained in practice using current MCE peak ground accelerations, and could be implemented via a new liquefaction-targeted maximum considered earthquake (MCE L ) peak ground acceleration ( PGA L ). The proposed PGA L is derived from probabilistic liquefaction hazard curves and includes uncertainties inherent to liquefaction hazard modeling. The potential design effects of the proposed PGA L are demonstrated via comparison with current MCE PGA levels, effects on first-level liquefaction screening analyses, and improvements in the consistency of P L,max levels across the United States, and for a wide range of site conditions.

conterminous United States

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

Simulation of groundwater flow in Wake County, North Carolina, 2000 through 2070

In 2019, the U.S. Geological Survey and Wake County Environmental Services began a collaborative study to evaluate groundwater resources and long-term groundwater availability in the county’s fractured-rock groundwater system. Wake County, in central North Carolina, is experiencing rapid population growth, associated land development, and changing water use. Hydrogeologic data including groundwater levels, aquifer testing, borehole fracture flow measurements, water-quality samples, and groundwater age-dating tracers were collected, along with findings from previous investigations, to help inform a conceptual model of the flow system used to develop a modular three-dimensional finite-difference groundwater-flow model (MODFLOW) for simulating historical and future groundwater conditions from 2000 to 2070. Hydraulic conductivity and transmissivity ranges were estimated from 17 slug tests and 21 borehole-flow measurements. Groundwater-quality analytical results from 19 sampling sites indicate that oxidation-reduction (redox) conditions varied within the regolith and bedrock and that minimal evaporation occurred before recharge entered the groundwater system. Age dating revealed mixtures of older and younger water, ranging from the 1940s to the 1990s—indicating variable flow pathways of recharge within permeable bedrock fracture zones. To simplify the complex fractured-rock groundwater system, two layers representing the regolith and the fractured bedrock were used in the MODFLOW model. Model calibration included parameter estimation and provided a reasonable fit to observed groundwater levels and estimated stream base flows. The model forecast scenarios incorporated future climate-model data for two emissions scenarios with land cover change projections to simulate potential impacts to future groundwater levels, recharge, and base flows. Recharge and base flow projections were largely within historical ranges, with no apparent long-term trends, but did indicate a slight downward shift in median values—likely, in part, because of differences in spatial resolution of input climate datasets. Seasonal patterns were consistent with historical data, with projections of possible increases in future winter recharge. Model limitations are discussed, and additional monitoring and model refinement needs are highlighted to support decision making for local groundwater management.

North Carolina

Revision of ModelMuse to support the use of PEST software with MODFLOW and SUTRA models

Executive Summary ModelMuse is a graphical user interface for several groundwater modeling programs. ModelMuse was updated to generate the input files for the parameter estimation software suite PEST. The software is used with MODFLOW or SUTRA models to run PEST-based parameter estimation and display the updated model inputs after parameter estimation. The PEST input files can also be used with the PEST++ version 5 software suite. Parameter estimation typically requires defining the parameters being adjusted during calibration and observations for assessing calibration quality. After a parameter is defined in ModelMuse, it can be applied to all or part of a model dataset. Pilot points—a parameterization device that facilitates higher levels of parameterization—can be used to assign spatially variable distributions of model inputs. Parameters can be applied to temporally varying features, such as boundary conditions, by either applying them to all the values in a series in one step or by applying separate parameters to individual members of a series. ModelMuse allows the definition of many observation types from various model output files. For MODFLOW 6 and SUTRA models, new options were added to ModelMuse to allow it to display the changed input after parameter estimation is complete. For MODFLOW–2005 and MODFLOW–NWT models, ModelMuse can import an entire model for visualization. An example illustrates the use of PEST with a MODFLOW 6 model in ModelMuse.

Techniques and Methods