USGS ScienceSearch

SEARCH · USGS Science

Results for “International Journal of Critical Infrastructure Protection”

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

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

1,661 records · Page 29Linked to original sources

Dissolved arsenic concentrations in surface waters within the upper portions of the Klamath River Basin, Oregon and California

Arsenic toxicity is an environmental health problem. Levels of arsenic in surface waters at some locations in the Klamath River Basin in southern Oregon and northern California can exceed the U.S. Environmental Protection Agency (EPA) standard for drinking water. There are both anthropogenic and natural sources of arsenic. The Klamath River Basin consists primarily of volcanic deposits and contains an underground geothermal system with hot springs and warm water wells, all known natural sources of arsenic. Anthropogenic sources of arsenic are related to the agricultural use of herbicides, fungicides, and insecticides. Surface water arsenic levels can also be affected by fertilizer amendments, evaporative concentration, oxygen-level depletion, and various geochemical transformations that can increase arsenic mobilization. In this study by the U.S. Geological Survey and the Bureau of Reclamation, dissolved concentrations of arsenic, copper, and lead were measured in surface waters at 39 unique sites within the upper portions of the Klamath River Basin between 2018 and 2022. In every year, except 2022, sites were sampled four times between April and November. Surface-water arsenic concentrations varied up to four-orders of magnitude among sites. Median arsenic concentration was lowest at Cherry Creek (0.03 micrograms per liter [μg/L]) and highest at Wood Kimball Spring (36.7 μg/L), two sites located north of Upper Klamath Lake. The highest arsenic concentrations (17.4±4.9 μg/L, n =3) were found in drain sites (defined here as a waterbody returning used irrigation water) while the lowest arsenic concentrations were found in an artesian well (0.8 μg/L, n =1). The elevated arsenic concentrations of the drain sites suggest that arsenic might be concentrated or mobilized by agricultural activities, water re-use practices, and (or) by geochemical processes occurring around water stored in drains (that is, in the water column and across sediment water boundaries). A source of arsenic in drain water in the Klamath Strait Drain area includes water used for irrigation originating from Ady Canal. Other potential sources include groundwater, geothermal water, and local soils and sediments. Seasonal differences in surface-water arsenic concentrations were detected at 13 sites, 10 of which had higher arsenic concentrations in summer than in either spring or fall. The sites sampled around Upper Klamath Lake, the impounded rivers, one of the two canal sites, and 5 of the 14 river sites had higher surface-water arsenic concentrations in the summer than in either spring or fall. Surface-water arsenic concentrations from groundwater sources (that is, springs and in the artesian well) did not vary significantly among seasons (p-values greater than 0.1). Median surface-water concentrations of copper and lead ranged from 0.03 to 3.7 μg/L, and from 0.013 to 0.175 μg/L ( n =2–18), respectively. Dissolved concentrations of both metals were below acute toxicity endpoints reported by the EPA for freshwater animals. Surface-water arsenic concentrations varied independently from corresponding changes in surface-water lead or copper concentrations. However, arsenic concentrations measured in bed-sediment samples collected from a subset of sites located north of Upper Klamath Lake correlated strongly and significantly with the corresponding sedimentary lead concentrations ( p =0.015). Aqueous arsenic speciation measured in a subset of sites in 2019 and 2022 showed that all the arsenic existed as arsenic (V), the most oxidized arsenic species, and presumably, the least toxic. The highest proportions of arsenite (As(III)), the presumably most toxic arsenic species, relative to total arsenic concentrations were found at drain sites. Our assessment of dissolved arsenic concentrations in various surface-water bodies in the Upper Klamath River Basin reveals geographical areas of consistently low (below 2.1 μg/L), moderate (below 10 μg/L) and high (above 10 μg/L) surface-water arsenic concentrations. South of Upper Klamath Lake, surface-water arsenic concentrations were consistently higher than 20 μg/L at two drain sites located in an area of predominant agricultural land use with extensive water re-use practices. North of Upper Klamath Lake, surface-water arsenic concentrations greater than 20 μg/L were consistently measured at sites with limited nearby agricultural activities, suggesting a geogenic source. The consistently high arsenic levels from the Wood River at Jackson F. Kimball State Park, Fort Creek, and Crooked Creek, which are sites located at or near headwater spring sources, suggest a natural background source of arsenic. Water flowing downstream from this area could be a potential source of arsenic to Upper Klamath Lake and the Upper Klamath River.

California, Oregon

Overview of The SCEC/USGS Community Stress Drop Validation Study using the 2019 Ridgecrest earthquake sequence

We present initial findings from the ongoing Community Stress Drop Validation Study to compare spectral stress‐drop estimates for earthquakes in the 2019 Ridgecrest, California, sequence. This study uses a unified dataset to independently estimate earthquake source parameters through various methods. Stress drop, which denotes the change in average shear stress along a fault during earthquake rupture, is a critical parameter in earthquake science, impacting ground motion, rupture simulation, and source physics. Spectral stress drop is commonly derived by fitting the amplitude‐spectrum shape, but estimates can vary substantially across studies for individual earthquakes. Sponsored jointly by the U.S. Geological Survey and the Statewide (previously, Southern) California Earthquake Center our community study aims to elucidate sources of variability and uncertainty in earthquake spectral stress‐drop estimates through quantitative comparison of submitted results from independent analyses. The dataset includes nearly 13,000 earthquakes ranging from M 1 to 7 during a two‐week period of the 2019 Ridgecrest sequence, recorded within a 1° radius. In this article, we report on 56 unique submissions received from 20 different groups, detailing spectral corner frequencies (or source durations), moment magnitudes, and estimated spectral stress drops. Methods employed encompass spectral ratio analysis, spectral decomposition and inversion, finite‐fault modeling, ground‐motion‐based approaches, and combined methods. Initial analysis reveals significant scatter across submitted spectral stress drops spanning over six orders of magnitude. However, we can identify between‐method trends and offsets within the data to mitigate this variability. Averaging submissions for a prioritized subset of 56 events shows reduced variability of spectral stress drop, indicating overall consistency in recovered spectral stress‐drop values.

California

Hydrogeology, karst, and groundwater availability of Monroe County, West Virginia

Monroe County is in southeastern West Virginia, encompassing an area of 474 square miles. The area consists of karst and siliciclastic aquifers of Ordovician, Silurian, Devonian, and Mississippian age and is in parts of two physiographic provinces: the Valley and Ridge Province to the east of Peters Mountain, and the Appalachian Plateau Province to the west of Peters Mountain. This study was developed in response to inquiries from the Monroe County Commission requesting assessment of the water resources of the county to better understand the quantity of the county’s groundwater resources, for both current [2023] and future demand, and to provide information to support protection and management of the county’s valuable groundwater resources. Various products were developed for this study that provide knowledge with respect to water availability and contamination susceptibility of the karst aquifers within the county. U.S. Geological Survey (USGS) geologists conducted extensive geologic mapping in support of the project, producing (1) a countywide bedrock geologic map, (2) a countywide hydrogeologic map, and (3) a light detection and ranging (lidar)-derived countywide digital elevation model and associated sinkhole map. A significant part of this work was to map in detail the Greenbrier Group at the formation level, which prior to this study had only partially been completed. The report also includes (4) a description of the lithologic units identified as part of the geologic mapping process. U.S. Geological Survey hydrologists completed several additional products for the hydrology part of the effort, including development of (1) a countywide potentiometric surface (water-table) map, (2) a countywide base-flow stream assessment, (3) countywide water-budget estimates, (4) well log surveys for 15 wells to better understand subsurface controls on groundwater flow within the study area, (5) two groundwater tracer tests to better refine the groundwater divide from the northern and southern parts of the karst aquifer in Monroe County; and finally, based on all available data collected for the study including the potentiometric surface map, geologic map, current [2023] and legacy fluorometric groundwater tracer tests, and base-flow stream assessments, (6) groundwater-basin delineations were reassessed for principal groundwater basins within the Greenbrier aquifer. In Monroe County, four principal hydrogeologic settings produce large yields of water for residential, agricultural, and other uses. The most relied upon water-bearing zone with respect to current [2023] public water supply is from springs along Peters Mountain. These springs are derived from intervals of fractured sandstone and resultant alluvial deposits. Groundwater flows downslope through these permeable alluvial deposits and discharges at the contact with less permeable strata, such as the Reedsville Shale. The second most relied upon water-bearing zone in Monroe County is within the karstic Greenbrier Group aquifer, in which the basal Hillsdale Limestone overlies the less permeable Maccrady Shale. This geologic contact between the Hillsdale Limestone and Maccrady Shale is not only targeted as a source of water for agricultural supply but also is targeted as a source of water for residential supply. The third most relied upon water-bearing zone is composed of shallow perched aquifers within the Greenbrier Group. The discontinuous nature of these perched aquifers makes mapping their extent impossible, but they are related to permeable geologic strata, such as karstified limestones with solutionally enhanced permeability that overlies less permeable shale or chert bedrock. During geologic mapping of the county, several of these perched aquifers were documented in the Pickaway, Union, and Alderson Limestones. A fourth zone consists of springs from Ordovician carbonates at the base of Peters Mountain, which are influenced by sinking streams as well as upwelling along faults. In terms of water quantity, the most sustainable springs are those having deeper-sourced flows. Public supplies are a principal source of water used for residential and commercial supply in the region, accounting for 0.49 million gallons per day (Mgal/d) of fresh-water withdrawals (0.14 Mgal/d of groundwater and 0.35 Mgal/d of surface water) for residential and commercial use and serving 6,645 individuals (49.2 percent of the population). An estimated 6,861 people, (50.8 percent of the population) primarily rely on private wells or other unregulated sources, such as springs, and withdraw 0.55 Mgal/d of groundwater for their residential use. Public water supply in the region is primarily (71.4 percent) derived from springs and augmented by stream withdrawals (backup sources mainly during low-flow periods), with the remaining portion (28.6 percent) derived from groundwater withdrawals from wells. For rural residents, however, 100 percent of their withdrawals are derived from groundwater (wells or springs).

West Virginia

Isotopic evaluation of the National Water Model reveals missing agricultural irrigation contributions to streamflow across the western United States

The National Water Model (NWM) provides critical analyses and projections of streamflow that support water management decisions. However, the NWM performs poorly in lower-elevation rivers of the western United States (US). The accuracy of the NWM depends on the fidelity of the model inputs and the representation and calibration of model processes and water sources. To evaluate the NWM performance in the western US, we compared observations of river water isotope ratios ( 18 O / "> 16 O and 2 H / "> 1 H expressed in δ notation) to NWM-flux-estimated (model) river reach isotope ratios. The modeled estimates were calculated from long-term (2000–2019) mean summer (June, July, and August) NWM hydrologic fluxes and gridded isotope ratios using a mass balance approach. The observational dataset comprised 4503 in-stream water isotope observations in 877 reaches across 5 basins. A simple regression between observed and modeled isotope ratios explained 57.9 % ( δ 18 O) and 67.1 % ( δ 2 H) of variance, although observations were 0.5 ‰ ( δ 18 O) and 4.8 ‰ ( δ 2 H) higher, on average, than mass balance estimates. The unexplained variance suggest that the NWM does not include all relevant water fluxes to rivers. To infer possible missing water fluxes, we evaluated patterns in observation–model differences using δ 18 O diff ( δ 18 O obs − δ 18 O mod ) and d diff ( δ 2 H diff - 8 ⋅ δ 18 O diff "> ). We detected evidence of evaporation in observations but not model estimates (negative d diff and positive δ 18 O diff ) at lower-elevation, higher-stream-order, arid sites. The catchment actual-evaporation-to-precipitation ratio, the fraction of streamflow estimated to be derived from agricultural irrigation, and whether a site was reservoir-affected were all significant predictors of d diff in a linear mixed-effects model, with up to 15.2 % of variance explained by fixed effects. This finding is supported by seasonal patterns, groundwater levels, and isotope ratios, and it suggests the importance of including irrigation return flows to rivers, especially in lower-elevation, higher-stream-order, arid rivers of the western US.

Hydrology and Earth Systems Science

U.S. Geological Survey science strategy to address white-nose syndrome and bat health in 2025–2029

Since its discovery in 2006, the fungal disease known as white-nose syndrome (WNS) has killed millions of bats. Of the 47 bat species native to the conterminous United States, Alaska, Hawaii, and Canada, 12 have been affected by WNS, including 3 endangered species and 1 proposed endangered species. WNS has also been detected in 40 States and 9 Canadian Provinces. U.S. Geological Survey (USGS) scientists have been critical in identifying the causal fungus for WNS ( Pseudogymnoascus destructans [Pd]), characterizing the effects of WNS, and tracking the spread of Pd in many bat populations in North America. The mission of the USGS WNS and Bat Health Science Team is to deliver integrated science in order to build resiliency into free-ranging bat populations through more effective WNS management, build capacity for bat health science, and enhance bat health information sharing across USGS science centers and cooperative research units as well as with stakeholders. The USGS can play an important role in supporting regional and national capacity building by providing resources and guidance to local, State, and Tribal management entities and by providing tools to enhance disease management. The USGS Ecosystems Mission Area’s Biological Threats and Invasive Species Research Program is the lead Federal program for free-ranging wildlife disease research and surveillance. As of 2024, guided by the science priorities set by the WNS Steering Committee, USGS scientists are engaged in a nationwide response to WNS. This work is done in close coordination with our partners at the U.S. Fish and Wildlife Service, National Park Service, Bureau of Land Management, U.S. Forest Service of the U.S. Department of Agriculture, U.S. Department of Defense, as well as State and Tribal agencies. In addition to conducting WNS research, the USGS is mapping the spread of WNS and coordinating the North American Bat Monitoring Program (NABat) to understand how WNS and other stressors affect the status and trends of native bats across their range. The USGS is supporting the national WNS response through four science goals: (1) provide situational awareness on the health of bat populations; (2) conduct ecological studies of bats along the gradient of disease vulnerability; (3) contribute actionable science to enhance the resiliency of bat populations; and (4) implement an adaptive, holistic approach to bat health.

Circular

Quantitative mineral resource assessment of lithium pegmatite deposits in the northern Appalachian orogen, USA

Lithium demand is projected to increase more than 48 times by 2040 due to electric vehicle production and other energy storage needs. Most lithium production is outside of the USA, thereby increasing supply chain vulnerability. The combined end use importance and heightened supply risk of lithium make this lightest metallic element a critical commodity to the USA. To mitigate this supply risk, the US Geological Survey is actively assessing lithium deposits in the USA. Herein, we detail an assessment for lithium-mineralized pegmatites in the US northern Appalachian Mountains. Permissive tracts were generated by cross-referencing tectonic and geologic maps and mineral occurrence data with mappable criteria derived from generalized and region-specific lithium pegmatite ore deposit models; tracts were then ranked as having high, medium, or low permissibility. Available geophysical and geochemical data were found to be of minimal utility for this deposit type at the scale of the assessment. The number of undiscovered deposits were estimated and integrated into probabilistic simulations, which included an expanded and updated global grade and tonnage model of pegmatite-hosted lithium ore. The estimated total amount of undiscovered resources for the northern Appalachian Orogen has a median value of 1,410,000 metric tons of Li 2 O when considering moderate correlation across sub-regions. At a confidence level of 90%, a resource of at least 90,000 metric tons of Li 2 O remains undiscovered, and at a 10% confidence level, a resource of as much as 7,380,000 metric tons Li 2 O remains undiscovered. After applying an up-to-date economic filter to convert median contained lithium to recoverable material, a correlated total of 900,000 metric tons of Li 2 O may be economically extractable, equating to enough Li 2 O to provide the current annual US lithium supply deficit (presently obtained through net imports) for 127 years at 2025 rates of apparent consumption. This period of provision will inevitably shorten with projected increasing consumption rates, emphasizing that further research could be completed to better delineate regions of high lithium resource potential and support exploration and domestic production.

Connecticut, Delaware, Maine, Massachusetts, New H

Geochemical geodatabase of sedimentary strata (coal, coal-adjacent rocks, tuffaceous oil shale, phosphate-rich rocks) and produced water in the Uinta region, Utah and Colorado

The Geochemical Geodatabase of Sedimentary Strata (Coal, Coal-adjacent Rocks, Tuffaceous Oil Shale, Phosphate-rich Rocks) and Produced Water in the Uinta Region, Utah and Colorado, consists of compiled datasets acquired as part of the Carbon Ore, Rare Earth, and Critical Mineral (CORE-CM) Uinta Region assessment funded by the U.S. Department of Energy (DEFE0032046, 2021–2024; Birgenheier et al., 2024). The CORE-CM assessment focused on providing comprehensive geological and geochemical characterization of current and prospective sedimentary-hosted resources including coal, oil shale, phosphatic limestone, and produced water from oil and gas targets present in eastern Utah and northwestern Colorado (Figure 1). This Data Series includes a geodatabase that consists of analytical geochemical data collected September 2021 through December 2024 via portable X-ray fluorescence (pXRF), and laboratory measured analyses produced by inductively coupled plasma mass spectrometry (ICP-MS) and inductively coupled plasma optical emission spectroscopy (ICP-OES). The coal-related geochemical data are derived primarily from the Cretaceous Blackhawk Formation and Ferron Sandstone of Utah, and the Mesaverde Group of Colorado. Additional non-coal resources assessed include oil shale-bearing strata of the Eocene upper Green River Formation (Utah and Colorado), phosphate-rich limestone of the Permian Park City Formation (Utah) and produced water from oil and gas-bearing strata of the Eocene Green River and Wasatch Formations (Uinta Basin) and the Pennsylvanian Paradox Formation (Paradox Basin) (Table 1). The CORE-CM assessment included a wide range of lithologies present in the coal, oil shale, and phosphate geologic resource systems whether or not the specific lithology has current economic value. Geochemical analyses of produced water from oil and gas wells focused on current and emerging hydrocarbon targets in the central Uinta Basin and northern Paradox Basin. A total of 13,092 geochemical analyses from these geologic systems is provided in the included geodatabase. A series of coal quality data (e.g., composition and maceral analyses) is also included in the database and was digitized from archived coal samples from the Utah Geological Survey (Appendix A).

Colorado, Utah

Age and colony variation in Adélie penguin metapopulation vital rates: Insights from a 25-year mark–recapture study

Understanding how vital rates vary with age, life-history stage, and among populations is fundamental for predicting the demographic consequences of environmental change, especially in longer-lived species with complex life histories. These species often exhibit delayed maturity and iteroparity, making nuanced demographic insights critical for assessing their long-term viability. This study investigated age- and colony-related variation in Adélie penguin vital rates including survival, recruitment, and breeding propensity. We used mark–recapture data collected over 25 years (1996–2020) from three Adélie penguin breeding colonies that differed in population sizes and trends but comprised a metapopulation located at capes Royds, Bird, and Crozier on Ross Island, Antarctica. We used multi-state models to estimate survival and detection rates relative to reproductive state and breeding colony and estimated transition probabilities reflecting movements between reproductive states and colonies. Apparent survival varied by reproductive state, colony, and age and averaged 0.80 (SD = 0.02) at Bird, and 0.72 (SD = 0.03) and 0.73 (SD = 0.03) at Crozier and Royds, respectively, for pre-breeders age 2–7 years with strong declines in pre-breeder survival after age 8. We observed less age-related variation in survival of breeders and non-breeders, but we observed differences between colonies with lower survival for breeders (0.72 to 0.80) compared to non-breeders (0.75 to 0.82). The average probability of surviving the first 2 years after fledging ranged from 0.43 (SD = 0.14) at Royds and Crozier (0.43, SD = 0.10) to 0.55 (SD = 0.15) at Bird. Movement between colonies was highest for pre-breeders (0.00%–12.00% depending on age and colony) and lowest for breeders (<0.20%). We observed the lowest age-related recruitment rates at Royds, with recruitment at Crozier almost twice as high, and intermediate at Bird. Breeding propensity was highest at Crozier and lowest at Bird. Colony-specific variation in vital rates likely contributed differently to population trajectories, suggesting that care must be taken to extrapolate vital rate estimates across colonies even within a metapopulation. These findings also highlight the importance of considering age, life-history stage, and geographic variation when assessing population vital rates in species with complex life histories.

Frontiers in Ecology and Environment

A global database of soil microbial phospholipid fatty acids and enzyme activities

Soil microbes drive ecosystem function and play a critical role in how ecosystems respond to global change. Research surrounding soil microbial communities has rapidly increased in recent decades, and substantial data relating to phospholipid fatty acids (PLFAs) and potential enzyme activity have been collected and analysed. However, studies have mostly been restricted to local and regional scales, and their accuracy and usefulness are limited by the extent of accessible data. Here we aim to improve data availability by collating a global database of soil PLFA and potential enzyme activity measurements from 12,258 georeferenced samples located across all continents, 5.1% of which have not previously been published. The database contains data relating to 113 PLFAs and 26 enzyme activities, and includes metadata such as sampling date, sample depth, and soil pH, total carbon, and total nitrogen. This database will help researchers in conducting both global- and local-scale studies to better understand soil microbial biomass and function.

Scientific Data

Ground deformation and gravity for volcano monitoring

Introduction When magma accumulates or migrates, it can cause pressurization and related ground deformation. Characterization of surface deformation provides important constraints on the potential for future volcanic activity, especially in combination with seismic activity, gas emissions, and other indicators. A wide variety of techniques and instrument types have been applied to the study of ground deformation at volcanoes (sidebar, p. 2; Dzurisin, 2000, 2003, 2007). Geodetic instruments include continuously recording Global Navigation Satellite System (GNSS; of which the United States’ Global Positioning System is one example) stations (fig. D1), borehole tiltmeters, and interferometric synthetic aperture radar (InSAR) measurements (from satellites, occupied and unoccupied aircraft systems, and ground-based sensors). Additional geodetic measurements like continuous- and survey-mode gravity (fig. D2) can contribute substantially to interpreting these data. Borehole strainmeters (see chapter K , this volume, by Hurwitz and Lowenstern, 2024) also have outstanding utility for monitoring deformation, although because of cost and permitting challenges, we do not include them as part of standard volcano monitoring networks for U.S. volcanoes. Still other techniques like light detection and ranging (lidar), structure from motion, and optical satellite data can be used to derive gross topographic changes, which can be used to map volcanic deposits, infer eruption rates, and gain insights into the source processes associated with eruptive activity (see chapter G , this volume, on tracking surface changes caused by volcanic activity; Orr and others, 2024). Experience has shown that no single geodetic monitoring technique is adequate to detect and track the entire range of ground-motion patterns that occur at volcanoes, primarily because of the temporal and spatial diversity of volcano deformation (fig. D3). Similarly, the magnitude of surface deformation varies widely. Geodetic monitoring strategies should therefore include multiple techniques and instrument types to cover a wide range of spatial and temporal scales. In identifying recommendations for geodetic instrumentation for volcano monitoring networks, we attempted to maximize the diversity of instrument types to measure the full range of deformation signals and minimize their expense and number; thus, we do not include several well-known deformation-monitoring techniques in our recommendations. Extensometers, for example, measure strains over distances of a few meters and have an excellent record of success in detecting changes in preeruptive localized ground motion across existing cracks, including at Mount St. Helens, Washington (Iwatsubo and others, 1992), and Piton de la Fournaise, Réunion Island (Peltier and others, 2006). Despite being relatively inexpensive, extensometers are best used primarily when localized ground displacements (for example, ground cracks) need to be tracked, and are not necessary at all volcanoes. In considering volcano deformation monitoring strategies, two complicating factors are deserving of special attention. First, not all deformation is driven by subsurface magmatic activity—for example, at many large stratovolcanoes (for example, Mount Rainier), flank collapses and landslides are significant geologic hazards (Reid and others, 2001) that may occur even in the absence of magmatic activity. Monitoring the stability of volcanoes is thus another critical application of geodetic monitoring networks to inform hazard assessment. One of the most famous examples of edifice instability is the large flank collapse that initiated the May 18, 1980, eruption of Mount St. Helens. Deformation monitoring had detected a bulge on the north flank of the mountain in April 1980 that was expanding by several meters per day (Lipman and others, 1981). Given that flank collapses can happen at any time during a period of volcanic unrest (or even outside a period of unrest), the capability to assess edifice stability is critical. Second, although volcanoes are commonly treated as idealized structures that erupt from single points, like centralvent stratovolcanoes, many are characterized by long rift zones from which eruptions may originate, and distributed volcanic fields are characterized by broadly spaced vents. For example, linear dikes are common at Kīlauea, Mauna Loa, and between Mount Shasta and Medicine Lake in California. At Kīlauea, one of these linear dikes emerged more than 40 kilometers (km) away from the summit of the volcano during the lower East Rift Zone eruption in 2018. Other volcanic fields, like Lassen volcanic center, California, or the San Francisco Volcanic Field, Arizona, have many small vents spread over a wide area. Although the instrumentation guidelines presented in this chapter remain phrased for central-vent volcanoes, they should be modified as needed in the context of the eruptive characteristics of each individual volcanic system. Spatial analysis of geodetic network coverage could help to ensure adequate instrumentation in areas where volcanism can occur over a broad area as opposed to a central vent. As an example, consider the adjacent volcanoes Mount Shasta and Medicine Lake. If station locations are chosen based only on the distance from the centers of the volcanoes, then any geodetic anomalies between the two volcanoes—an area of potential volcanism as indicated by the presence of volcanic features—may remain undetected by ground-based instrumentation. The spatial analysis is accomplished via a grid of pressure point sources (Mogi, 1958) evenly distributed across the map area, at a depth of 5 km in this example (fig. D4). Each source is inflated until predicted deformations exceed the GNSS white noise uncertainty estimates at one site (Langbein, 2017; Murray and Svarc, 2017). This volume of detectable magma provides a measure of the quality of the coverage (fig. D4). The results indicate that, as of 2022, there is a large area between Mount Shasta and Medicine Lake volcano with existing mapped dikes in which a substantial amount of magma could intrude without being detected geodetically. Applying this style of analysis to individual volcanic systems can provide a guide for designing network geometry given the expected locations of future eruptions.

Scientific Investigations Report

A time-space model of graphite mineral systems

Increasing demand for graphite in energy storage systems warrants review of graphite ore genesis in a mineral systems framework. Orogenic graphite encompasses the metamorphic and orogenic mineral systems that produce flake graphite and hydrothermal vein (lump and chip) graphite deposits, respectively. A common feature of orogenic graphite deposits is an association with upper amphibolite- to granulite-facies metasedimentary rocks in continent-continent or continent-island arc collisional orogens. Orogenic flake graphite deposits form primarily through graphitization of organic carbon during regional metamorphism, but strain localization and partial melting of pelitic protoliths are likely important processes for graphite grade and quality enrichment. Orogenic vein graphite deposits precipitate from hydrothermal fluids, possibly derived from metamorphism or anatexis at depth. Decarbonation reactions in mixed calcareous-carbonaceous metasedimentary sequences are the most likely carbon sources for the veins. In contrast, intrusion-related graphite includes magmatic-hydrothermal and metamorphic mineral systems that form primarily in continental arc settings via the interaction of magmas with carbonaceous sedimentary packages. Magmatic-hydrothermal flake graphite deposits are hosted in plutonic and volcanic rocks, and result from the exsolution of CO 2 -CH 4 -rich fluids from melts contaminated by such packages. Contact metamorphism of carbonaceous sedimentary rocks by plutons produces some microcrystalline (amorphous) graphite deposits, including many in China. Compilation of geologic data from known graphite deposits globally suggests that pulses of carbon deposition in the Paleoproterozoic, Mesoproterozoic, and Neoproterozoic provided source carbon material. Subsequent supercontinent orogenesis at ca. 2,100 to 1,700 Ma (Columbia), ca. 1,300 to 1,000 Ma (Rodinia), and ca. 650 to 500 Ma (Gondwana) resulted in the genesis of orogenic flake and vein graphite deposits, where favorable geologic components overlapped with organic ± carbonate carbon-rich strata. Cryogenian deposition of graphite protoliths and Cryogenian – Cambrian metamorphic mineralization account for nearly 75% of all known resources globally and coincide with profound carbon isotope excursions and climate variability, implying a link with the global carbon budget. Comparatively few graphite deposits are associated with Pangea-forming orogens, attributed to less exhumation and/or denudation. High-temperature metasedimentary belts containing organic carbon-rich protoliths are most favorable for hosting orogenic flake graphite deposits, whereas sequences that also contain carbonate protoliths are favorable for orogenic graphite veins. Continent-scale orogenic belts may host both deposit types along with vanadium deposits. Use of a time-space mineral systems framework for graphite deposits can improve exploration models needed to ensure future supply of this critical mineral and provide insights into Earth’s long-term carbon cycle.

Mineralium Deposita

The value of electronic tagging and tracking studies for understanding fish–hypoxia interactions

Objective Hypoxia (i.e., low dissolved oxygen availability) is a natural phenomenon but can also be induced by human activities (e.g., nutrient enrichment from runoff). Given that dissolved oxygen is essential for aquatic life, periods of hypoxia tend to have negative consequences (e.g., sublethal disturbances, mortality) for most fishes. Extensive laboratory research has documented hypoxia thresholds and physiological and behavioral outcomes for a variety of freshwater and marine fishes and there is also an extensive body of fieldwork assessing population-level responses (e.g., survival, distribution). Yet, studying how individual fish respond to hypoxia in the wild has proved challenging; electronic tagging and tracking tools (e.g., biotelemetry, biologging) have made it easier to study individual responses to hypoxia and complement other tools like hydroacoustics that tend to focus on population-level responses. Methods We review what has been learned from contemporary studies that employ electronic tagging and tracking tools to understand how fish respond to hypoxia in the wild. Topics explored include identifying and validating hypoxia thresholds, habitat compression, connectivity, mortality, physiological and bioenergetic consequences, and extreme weather conditions. We also consider what we have learned about evaluating various management plans for hypoxia and reflect on how electronic tags have also been used in aquaculture systems and in hybrid studies that combine laboratory and field research. We highlight fishes in Lake Erie as a research narrative to demonstrate how electronic tagging and tracking have markedly improved our understanding of a longstanding hypoxia issue. Results Our synthesis revealed that electronic tagging and tracking have provided critical information on how fish respond to hypoxia in the field, revealing complex trade-offs and compensatory mechanisms as well as cryptic hypoxia-induced mortality. Beyond just illuminating space use and mortality, tags equipped with various sensors are revealing how fish deal with hypoxia in real time in terms of physiology, bioenergetics, and behavior. Conclusions As electronic tagging and tracking methods experience further innovation and are increasingly applied to understand the effects of hypoxia on fish, we expect more unanticipated findings about the effects of hypoxia on fish in all aquatic ecosystems, which will strengthen our ability to manage and mitigate hypoxia. Combining tools and approaches (e.g., lab and field) is perhaps the best way to generate comprehensive understanding.

Transactions of the American Fisheries Society

Editorial: From cold seeps to hydrothermal vents: Geology, chemistry, microbiology, and ecology in marine and coastal environments

This Research Topic compiles contemporary studies on cold seeps, hydrothermal vents, mud volcanoes, and related seafloor features that are associated with focused fluid emissions and the transfer of carbon, other chemical species, and sometimes heat from the geosphere to the ocean. Because these features sometimes tap fluids and gas originating kilometers below the seafloor, they provide an important window into deep processes that are otherwise inaccessible to scientists. At the shallow portion of their journey, migrating fluids nearing the seafloor contribute to a range of unique biological, physical, and chemical processes within the sediments themselves and at the sediment-water interface. Seafloor fluid emissions play a critical role in global biogeochemical cycles, ocean chemistry, and possibly even climate change. Seafloor leakage points often emit hydrocarbon gases (especially methane and CO 2 ) and are sometimes the loci for deposition of seafloor minerals that have economic value. A burgeoning area of research focuses on natural products generated at these features, seeking compounds with potential pharmaceutical or other applications. Multidisciplinary studies have become routine for characterization of seafloor fluid emission sites, attesting to the inseparability of geologic, physical, chemical, and biological processes in these settings. It is increasingly common for researchers to combine in a single research cruise: subbottom imaging and seafloor mapping; porewater and water column geochemistry and gas sampling; sediment retrieval for lithologic, biostratigraphic, and solid phase analyses; and studies of benthic and subseafloor communities at the microbial to macrofaunal scales. This multidisciplinary approach has the advantage of ensuring the spatial and temporal coincidence of surveys and samples, an important factor at highly dynamic seafloor fluid emission sites. In addition, researchers often use remotely operated vehicles (ROVs), autonomous underwater vehicles (AUVs), or human-occupied vehicles (HOVs) to record video of the seafloor, compile photomosaics, collect targeted samples, and survey with high-resolution geophysical near-seafloor systems, providing a degree of detail about seafloor fluid emission sites that is unprecedented compared to most areas of the deep ocean. While rarer, long-term cabled observatories or shorter-term deployments of portable observatories are also used at some loci for seafloor fluid flux and are particularly helpful for capturing temporal variations at these dynamic features. Here we summarize the Research Topic’s contribution to multidisciplinary seafloor emission studies in the categories of cold seeps, mud volcanoes, and hydrothermal vents. Figure 1 shows the geographic distribution of the studies in this Research Topic and key features referred to in this Introduction.

Frontiers in Earth Science

Insights and strategic opportunities from the USGS 2024 Per- and Polyfluoroalkyl Substances (PFAS) Interagency Workshop

Introduction In 2021, the U.S. Geological Survey (USGS) published Circular 1490 titled, “Integrated Science for the Study of Perfluoroalkyl and Polyfluoroalkyl Substances (PFAS) in the Environment: A Strategic Science Vision for the U.S. Geological Survey” (Tokranov and others, 2021). Circular 1490 was created to be a resource for USGS scientists prioritizing and planning research related to per- and polyfluoroalkyl substances (PFAS) and to be a guide for developing partnerships with other scientists, State and Federal agencies, and stakeholders engaged in PFAS research and management and mitigation of the environmental and human-health effects of PFAS. This USGS PFAS Strategic Science Vision document was intended to be the foundation for a “living strategic vision,” periodically providing updates on the state of USGS PFAS research, emerging PFAS data gaps and needs, and progress on interagency and stakeholder PFAS partnerships and priorities. To meet this objective, the USGS planned to host an Interagency and Stakeholder PFAS Workshop every 2–3 years. During September 10–12, 2024, the USGS hosted the first Interagency and Stakeholder PFAS Workshop in Reston, Virginia. The Workshop brought together experts from other Federal agencies (U.S. Environmental Protection Agency, National Institute of Environmental Health Sciences, Food and Drug Administration, Department of Defense [Air Force, Army]), State agencies (Washington Fish and Wildlife, Virginia Department of Transportation), and academia (Harvard University, University of Maryland) to address key challenges relating to the measurement and modeling of PFAS and the implications for environmental health. Participants engaged in in-depth discussions centered around six pivotal topics related to PFAS: (1) sampling protocols, methods and interpretation; (2) environmental sources, source apportionment, and occurrence; (3) environmental fate and transport; (4) human and wildlife exposure routes and risk; (5) bioconcentration, bioaccumulation, and biomagnification; and (6) ecotoxicology and effects. Each topic had three breakout sessions. A recurrent theme of workshop discussions was how data on a nationwide scale for PFAS occurrence in various environmental matrices, including air, water, food crops, biota, soil, and streambed sediment could help to advance scientific understanding. Participants noted significant geospatial data gaps, particularly in the midwestern and southern United States and the Pacific Northwest. PFAS data collection tends to be more robust along the eastern seaboard and in California. Participants stressed how enhancing the integration of large and small datasets across various agencies could help to support national scale understanding of PFAS. To address these gaps, attendees suggested leveraging datasets from Federal entities like the USGS and the U.S. Department of Defense, State agencies, and municipal utility services to develop predictive contaminant detection and transport models. Improved coordination between water quality programs and USGS research could help to facilitate access to valuable data, leading to comprehensive databases that inform PFAS point (wastewater treatment plants and landfills) and nonpoint (runoff from land, atmospheric deposition, food packaging) sources, environmental transport mechanisms, environmental detection and concentrations, potential exposure routes, and health effects on different biota, including humans. A specific request was made to develop a map demarking the depth of modern (1953 or later) groundwater, which is susceptible to surface-derived anthropogenic (that is, human-made) contamination, based on tritium-age dating. Emphasis was placed on incorporation of hydrology, groundwater flow paths, groundwater–surface water interactions, and landscape factors in predictive statistical models as a step to improve contaminant source identification and tracking. Molecular fingerprinting approaches garnered attention as techniques to link specific PFAS mixtures detected in a sample to environmental sources and levels in biota (Dávila-Santiago and others, 2022). Integrating data from abiotic (that is, water, soil, and air) and biotic (that is, living organisms) systems identified as a research opportunity. For example, understanding the composition of soils and sediments, which include a mixture of mineral, plant, and animal components, could advance understanding of exposure pathways. The discussions highlighted opportunities to explore and understand the potential redistribution and biotic exposures of PFAS from biosolid and wastewater treatment plant effluent land application practices, in addition to atmospheric releases and discharges from landfill and wastewater treatment plants. Participants identified research gaps surrounding how these sources may contribute to contamination and may affect surrounding ecosystems, including a better definition of anthropogenic background concentrations. Moving forward, the collection of co-occurrence data was noted as a means to improve understanding of complex mixtures and to leverage companion modeling efforts focused on areas with high and low contamination levels to identify areas of concern and unaffected resources. Participants emphasized how centralized USGS databases and the establishment of sample-metadata archives can help to ensure that samples are preserved and accessible for future research. In conclusion, the workshop participants identified opportunities to bridge data gaps and improve measurement techniques, modeling frameworks, databases, and communication, to enhance the understanding of PFAS and their effects on environmental and human health. Upon completion of the workshop, participants indicated an interest in developing strategic data collection, modeling, and analytical approaches to address these challenges.

Open-File Report

Evaluation of best management practices at an edge-of-field site in the Eagle Creek watershed, Ohio, 2012–20

In 2010, the U.S. Geological Survey worked in partnership with the Great Lakes Restoration Initiative and the Natural Resources Conservation Service to identify farm fields in priority watersheds—watersheds critically important to Great Lake health. Intensive best management practices were implemented in five States to test the efficacy of best management practices, also referred to as agricultural conservation practices, to reduce sediment and nutrient runoff on agricultural fields. These fields were chosen as representatives of the priority watersheds because their farming practices and geographical conditions were common among farms in those watersheds. A farm site located near Findlay, Ohio, in the Eagle Creek watershed used two conservation practices, cover crop (Natural Resources Conservation Service practice standard 340) and variable rate technology (VRT; Natural Resources Conservation Service practice standard 590), which represent common practices in the Maumee priority watershed. This study monitored surface runoff and subsurface tile runoff at the Eagle Creek watershed site. The effects of cover crop and VRT nutrient application on sediment and nutrient runoff were assessed from October 2012 to September 2020. Cover crops were applied in fall of 2016, 2017, and 2018, coinciding with VRT applications. The following parameters were analyzed as part of the study: runoff discharge; peak discharge; and concentrations and loads of total phosphorus, particulate phosphorus, dissolved reactive phosphorus, and suspended sediment. High peak runoff events disproportionally affected nutrient concentrations and loads in runoff. A threshold of the 85th percentile of the peak discharge was chosen to assess the effect of cover crop and VRT nutrient application on nutrient mitigation at each gage. Runoff events below the 85th percentile were considered small, and runoff events over the 85th percentile were considered large. Concentrations from collected water-quality samples were used for the analysis; estimated concentrations were excluded from the statistical evaluation. All surface runoff parameters during small runoff events, except surface runoff volume, were significantly lower in the period with cover crops and VRT implementation. Contrarily, subsurface runoff parameters during small events were not statistically different between the periods with and without cover crops and VRT implementation. Surface and subsurface parameters during large events were also not statistically different between the periods with and without cover crops and VRT implementation. These results indicate cover crops and VRT may improve water quality during small runoff events, but additional best management practices that mitigate large runoff events may lead to greater water-quality improvement given the contribution of large events to overall losses.

Ohio

Assimilation of reduced carbon triggers platinum alloy saturation in mafic and ultramafic magmas

It is generally observed that magmatic sulfide ores have higher ratios of Pd/Pt than the mantle-like values of their parental magmas. This discrepancy has defied simple explanation because the partitioning behavior of both elements between sulfide and silicate liquids is very similar. Assimilation of sulfur- and carbon-rich country rocks by mafic and ultramafic magmas is considered a critical, if not essential, step in the formation of magmatic base metal sulfide deposits. Although there is general consensus that the assimilation of external sulfur and carbon promotes sulfide saturation, the effect of carbon assimilation on the solubilities of platinum-group elements in natural S-bearing silicate melt has been overlooked. In this study, we investigate the variations of platinum and palladium solubilities during assimilation of graphite and methane through thermodynamic modeling, in comparison with data from an array of highly distinctive magmatic sulfide ore systems representing ages from Archean to Paleozoic, melt compositions from komatiite to basalt, and magmatic settings including lavas, hypabyssal intrusions, plutonic continental arc roots, and plutonic layered intrusions, namely: Raglan, Norilsk-Talnakh, Lac des Iles, and the J-M Reef of the Stillwater Complex. We model assimilation-fractional crystallization processes to estimate the reduction of oxygen fugacity ( ⁠fO2⁠ ) of the melt due to incorporation of graphite and methane. The simulations show that although Pd remains highly soluble during the progressive assimilation of reduced carbon, Pt solubility decreases significantly as the silicate melt becomes increasingly reduced. With less than 8 % of sediment assimilation, Pt alloy may saturate and then deviate from sulfide-undersaturated silicate melts, concomitantly increasing the Pd/Pt value of the remaining melts of the Raglan and Norilsk-Talnakh systems. For the Lac des Iles and Stillwater systems, a higher extent of assimilation is needed to reach Pt saturation because of the relatively carbon-poor nature of the lower crustal rocks. The assimilation of methane volatiles is shown to be more effective than graphite assimilation, and it provides a pathway to Pt alloy fractionation in the absence of detectable amounts of bulk host-rock assimilation. High Pd/Pt values have been documented in many world-class magmatic sulfide deposits whose parental magmas have demonstrably experienced crustal contamination. Our model suggests that although anomalous Pd/Pt values may be explained by other mechanisms such as incongruent melting of preexisting sulfide or differences in the diffusivities of the metals within achieving equilibration, the assimilation of graphite or methane may play an important role in the global occurrence of magmatic sulfide ores with elevated Pd/Pt values.

Economic Geology

Constraining landslide frequency across the United States to inform county-level risk reduction

Informative landslide hazard estimates are needed to support landslide mitigation strategies to reduce landslide risk across the United States. Whereas existing national-scale landslide susceptibility products assess where landslides are likely to occur, they do not address how often , which is a critical element of landslide hazard and risk assessments. In particular, the U.S. Federal Emergency Management Agency's National Risk Index (NRI) requires landslide frequency estimates to inform expected annual loss estimates. We present county-level landslide frequency (landslides per area per year) estimates for the 50 US states. We applied Bayesian negative binomial regression to estimate both the expected (average) reported landslide frequency and full distribution of annual landslide counts for each county. We compared a suite of models that used combinations of landslide-susceptible area, probability of potentially triggering earthquakes, frequency of potentially triggering precipitation, and ecological region as predictors. We trained our models with landslide inventory data from counties with the most comprehensive records available nationwide and used zero-inflated negative binomial distributions as an incompleteness model to correct for temporal reporting gaps. We selected a preferred frequency model to inform the NRI based on information criteria and physically plausible parameter estimates. The model showed that average annual reported landslide frequencies vary by 5 orders of magnitude across US counties, ranging from 0.002 (0.00015–0.05) landslides 1000 km −2 yr −1 in Kusilvak Census Area, Alaska, to 29 (19–46) landslides 1000 km −2 yr −1 in Lake County, California, reflecting the country's strong variations in landslide susceptibility, earthquake probability, and other factors for which ecological region serves as a proxy. Counties with estimated frequencies in the top 20 % of all counties are predominately along the West Coast of the continental United States, in mountainous regions of the Pacific Northwest and Intermountain West, in locally steep or earthquake-prone regions of the Midwest and Southeast, along the Appalachians, in southern and southeastern Alaska, and on some Hawaiian islands. By examining the number of landslides predicted in 99th percentile years for each county, we identified that 26 % of US counties likely have potential for widespread landsliding with more than 10 landslides 1000 km −2 yr −1 , even when such large events have not been reported in the training data for that county. Overall, our results better represent the range of possible landslide frequencies and spatial variations than previous national-scale estimates reported in the NRI, and our approach can inform other risk-reduction and loss-mitigation efforts across the United States and globally.

Natural Hazards and Earth System Sciences

Assessing spatial variability of nutrients, phytoplankton, and related water-quality constituents in the California Sacramento–San Joaquin Delta at the landscape scale—2018 high resolution mapping surveys

Executive Summary This study examined the abundance and distribution of nutrients and phytoplankton in the tidal aquatic environments of the Sacramento–San Joaquin Delta (Delta) and Suisun Bay, comprising three spatial surveys conducted in May, July, and October of 2018 that used continuous underway high frequency sampling and measurements onboard a high-speed boat to characterize spatial variation across the extent of the Delta. The method used involves simultaneously collecting information about the concentration and spatial distribution of all major nutrient forms with analogous information about the major classes of phytoplankton and associated water-quality conditions. The results showed substantial variation across space and time, providing an unprecedented snapshot of the dynamic environmental processes that shape the ways nutrients interact with and affect aquatic habitats in the Delta. The purposes of this study were to improve our understanding of how hydrodynamics, landscape features, and aquatic primary productivity interact to drive nutrient cycling and transport in the Delta and to provide insights into the underlying processes most directly responsible for the conditions at the time of this study, and thus into the range of conditions that may be expected following the wide array of prospective future changes to the Delta. One major anticipated change at the time of this study was the planned upgrade to the Sacramento Regional Wastewater Treatment Plant, but the study also informs our understanding of potential effects from other changes to the Delta, such as those caused by other nutrient-management actions, flow actions, large-scale wetland restoration, drought, flood, levee failure, and changes to water management. Nutrient loading is the primary driver of nutrient concentrations in the Delta, but several other major drivers interact to shape their distribution and effects: geomorphology, hydrodynamics, landscape features, and aquatic productivity. Hydrodynamics affect timescales of transport and dilution of nutrient loads in the Delta. During transit through the system, channel geometry, tidal mixing, and water exports affect hydrodynamics in diverse ways that influence water-residence and transport times, thereby markedly affecting the range of times during which natural internal cycling can alter nutrient concentrations and forms. Channel geometry and location shape tidal energy and river currents into these observed dynamics. Interactions with Delta aquatic landscapes such as herbaceous tidal marsh, submerged aquatic vegetation, and large expanses of intertidal or subtidal sediments (all highly productive landscapes) exert demand on available nutrient supplies but can also simultaneously transform and generate nutrients. Finally, while phytoplankton require nutrients to sustain production and thus are a potential nutrient sink, the amount and form of nutrients also can influence the occurrence of harmful algal blooms (HABs) that adversely affect aquatic organisms as well as affect the occurrence of beneficial algal blooms that result in production of algae that are favorable for imperiled Delta pelagic aquatic food webs. The surveys revealed a complex mosaic of spatial variation, with nutrient concentrations varying from near zero to well above concentrations considered eutrophic; nutrient concentrations were more often related to the extent of hydrologic transport and mixing than to specific geographic locations or to specific landscape features. Similarly, the surveys identified phytoplankton abundance ranging from near detection to the level of large phytoplankton blooms, with large variation in phytoplankton community composition. Although the study occurred during a period of low bloom activity, phytoplankton productivity appeared to be the strongest potential sink for inorganic nutrients in the Delta, indicating that it is a larger control on nutrient concentrations and distribution than previously understood. Cycling and transformation within the water column only appeared to substantially lower total nutrient concentrations at the longest estimated transport timescales. Contrary to expectations, we did not observe substantial nutrient depletion near landscape-scale features such as open-water habitats, submerged aquatic vegetation beds, extensive wetlands, or exposed sediments, indicating that these habitat types did not act as major sinks for nutrients in the Delta during these surveys. These results indicated that nutrient reduction efforts may have the greatest effect on pelagic phytoplankton productivity in the more productive reaches of the Delta and estuary, but these effects are unlikely to be magnified by changes to nutrient loss within the Delta over conceivable changes in flow conditions, Delta water management actions, or large-scale wetland restoration activities. Nevertheless, local processes were shown to cause substantial loss, and thus integrating of nutrient effects with other indicators of aquatic habitat conditions will help inform planning future actions at specific sites. Finally, we note that the primary contribution of this study was intended to be the survey data themselves. Aside from the results highlighted in this report, the surveys are a benchmark against which future environmental change may be evaluated, including changes to nutrient management or water exports, drought, large-scale wetland restoration, and climate change. Further, although we highlight some of the main findings from the surveys in this report, the necessarily limited scope precludes examination of many topics for which these surveys may be highly informative. To facilitate the utility of these data to stakeholders, managers, and researchers, we have released the data online (Bergamaschi and others, 2020) and created an online data exploration portal ( https:​//ca.water​.usgs.gov/​bay-​delta/​2018-​delta-​wide-​mapping-​surveys.html ) where users may query the surveys in a variety of ways to test hypotheses, examine relationships, assess spatial trends, and download data. The data exploration portal is intended to be an immersive experience that allows users to gain greater understanding of the complex interactions that shape Delta aquatic environments. This report is intended as a companion to the portal, allowing the reader to challenge and further explore the highlighted findings. This study was a collaboration between the U.S. Geological Survey and the Delta Regional Monitoring Program, with additional funding provided from U.S. Geological Survey Cooperative Matching Funds Program.

California