USGS ScienceSearch

SEARCH · USGS Science

Results for “Journal of Research of the U.S. Geological Survey”

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,697 records · Page 27Linked to original sources

Cruise summary—Samoa Basin abyssal mapping—Box coring leg, Pago Pago, Territory of American Samoa to Pago Pago, American Samoa, April 11– May 1, 2026

Expedition Summary U.S. Geological Survey scientists led a box coring effort to the Samoa Basin to characterize minerals and the surrounding abyssal sediments and fauna. Thirty-eight box cores were deployed between April 13, 2026, and April 28, 2026. Thirty-six box cores recovered sufficient material to determine nodule density, and 35 recovered sufficient material for subcores to be collected. The purpose of this Data Report is to provide a summary of samples collected, initial results that were obtained shipboard, and briefly mention planned future analyses from this expedition.

Data Report

Footprints of past mining in Alaska (USA) derived from high-resolution satellite imagery

Mapping the land area used for mining in the past is essential for guiding the remediation of affected landscapes and assessing the resource potential of related waste products. Despite significant recent progress delineating footprints of active and inactive mining globally, the known inventory of such mine lands remains incomplete. Here, I describe a new map dataset of footprints of land surface disturbance and waste at sites of past mining in Alaska (USA) based on visual interpretation of satellite imagery. This dataset maps 6–14 times the area of previous regional and global mine footprint maps in Alaska and is the first in the region to explicitly delineate mine waste landforms (e.g., tailings piles). The data are publicly available from the U.S. Geological Survey under a “no rights reserved” Creative Commons (CC0) license agreement.

Alaska

Water resources related to breccia pipe uranium mining in the Grand Canyon region

Introduction In the arid Grand Canyon region, water resources are limited to primarily the Colorado River and associated tributaries and to groundwater in the form of seeps and springs. Groundwater resources in the region supply water for human use and support diverse and rich ecosystems in the locations immediately surrounding the seeps and springs. Throughout the region, uranium resources occur and may interact with water resources in both mined and unmined uranium deposits. There is a need to better understand groundwater in the region and the effects from uranium mining in order to better manage the limited water resources in the area. This Fact Sheet summarizes results from U.S. Geological Survey studies that were conducted on this topic from 2012 to 2023.

Arizona

An inset groundwater-flow model to evaluate the effects of layering configuration on model calibration and assess managed aquifer recharge near Shellmound, Mississippi

The U.S. Geological Survey has developed a high-resolution inset groundwater-flow model in the Mississippi Delta as part of an interdisciplinary collaboration coordinated by the Mississippi Alluvial Plain project to provide a tool that stakeholders can use to support water-resource management decisions. Groundwater withdrawals from the Mississippi River Valley alluvial (MRVA) aquifer have been vital to support agricultural production in the region, but substantial groundwater-level declines near Shellmound, Mississippi, have caused concerns for long-term sustainability of the aquifer. To better understand the subsurface and try to mitigate the long-term groundwater-level declines, stakeholders have undertaken actions including a Groundwater Transfer and Injection Pilot (GTIP) project using a riverbank filtration-based managed aquifer recharge approach. The pilot project consisted of extracting groundwater near the Tallahatchie River and reinjecting it into the aquifer 3 kilometers west where water levels have substantially declined. A high-resolution airborne electromagnetic (AEM) survey was also completed to collect electrical resistivity data to support the GTIP project and the development of the groundwater model. The inset groundwater-flow model was developed to (1) integrate the AEM data into the optimal layering configuration of the MRVA aquifer that the available observation data can support through calibration, and (2) assess the potential effect of the GTIP project on the groundwater levels. The AEM data were processed into three different layering configurations leading to the development of model A (18 layers), model B (16 layers), and model C (8 layers), all at a 100- x 100-meter cell spatial resolution using the U.S. Geological Survey modular finite-difference flow model 6 code with Newton-Raphson formulation. The model development process integrated recent advances in modeling, such as the incorporation of AEM data, the use of outputs from the soil-water-balance (SWB) model, and the Aquaculture and Irrigation Water-Use Model, and was facilitated by robust automation using the open-source python packages Modflow-setup and SFRmaker. Using Parameter Estimation ++ Iterative Ensemble Smoother, the three numerical groundwater-flow models (models A, B, and C) were calibrated against a set of observations, which included aquifer groundwater levels, streamflows, stream stage, and aquifer transmissivity. Results indicate that the detailed representation of MRVA aquifer layers in model A produced the best calibrated model by history matching, and the integration of data representing surficial connectivity played a key role in improving groundwater recharge and enhancing the ability of the model to match groundwater levels in the cone of depression. A forecast model simulated the managed aquifer recharge approach, and the results indicated that, given average irrigation and recharge conditions (2010–15), the GTIP project has the potential to induce groundwater-level increases of as much as 3 meters around the injection site, but a sustained increase would require repetition in subsequent years of water transfer at 2022 rates or above.

Mississippi

The Great Lakes Geologic Mapping Coalition—Working collaboratively to understand the geology of the Great Lakes Region

Introduction The Great Lakes Geologic Mapping Coalition (GLGMC), commonly referred to as the “Coalition,” is a partnership between the U.S. Geological Survey (USGS), the U.S. States of Illinois, Indiana, Michigan, Minnesota, New York, Ohio, Pennsylvania, and Wisconsin and the Canadian province of Ontario. The member States receive funding for geologic mapping work from the USGS National Cooperative Geologic Mapping Program (NCGMP), whereas Ontario participates as a nonfunded partner. The mission of the GLGMC is to produce three-dimensional (3D) geologic maps that depict unconsolidated sediments and near-surface bedrock in the Great Lakes region of North America. Geologic maps are the basis of most earth science investigations and help support resource exploration (energy, minerals, groundwater), natural hazard mitigation, infrastructure development, and land-use planning, all of which can be used to advance economic development and strengthen national security in the Great Lakes region. During the last few million years, the Great Lakes region has experienced repeated glacial advances and retreats, leaving behind extensive sediments, abundant natural resources, and widespread effects on the underlying bedrock geology (Swezey and others, 2022). Linked by shared histories of past glaciations, industrial agriculture, and legacy automotive, coal, steel, and manufacturing industries, the GLGMC member States collaborate to improve the understanding of the 3D distribution of the sediments overlying the region’s bedrock (fig. 1). Developing a comprehensive subsurface 3D framework of this glaciated terrain can provide earth science data to policymakers at all levels. These insights facilitate informed decisions on the exploration, use, and protection of vital resources, such as critical minerals, industrial materials, and aquifers, thereby supporting economic prosperity and the well-being of the citizens of this region. Since its inception in 1998, the Coalition has completed more than 100 geologic mapping projects across the Great Lakes region. Each project aims to deliver geologic maps, 3D datasets, and other information that improves understanding of the geology of the Great Lakes region, with an emphasis on economic and water resources. Key deliverables include 3D geologic maps and models typically portraying sediment thickness, often derived from top-of-bedrock and borehole data. These products are developed through a combination of fieldwork, subsurface modeling, and the collection and analysis of rock and sediment cores. To support Coalition goals, member States collaborate with scientists working on related STATEMAP, EDMAP, and FEDMAP projects. Coalition scientists also engage with Tribal Nations in the Great Lakes region to ensure that Tribal interests pertaining to Coalition work are addressed. Through this collaboration, the Coalition unites the efforts of State, Federal, and Tribal Nation stakeholders to advance geologic data production and enhance understanding of the geologic resources of the Great Lakes region.

Illinois, Indiana, Michigan, Minnesota, New York,

Methodology and technical input for the 2025 U.S. List of Critical Minerals—Assessing the potential effects of mineral commodity supply chain disruptions on the U.S. economy

The Secretary of the Interior, acting through the Director of the U.S. Geological Survey, is tasked by section 7002 (“Mineral Security”) of title VII (“Critical Minerals”) of the Energy Act of 2020 (Public Law 116–260, December 27, 2020, 116th Congress) with reviewing and revising the methodology used to evaluate mineral commodity supply risk and the U.S. List of Critical Minerals (LCM) no less than every 3 years. Following two previous LCM assessments, this analysis represents the latest technical input for evaluating each mineral commodity’s supply risk and determining their recommended status on the LCM. We evaluated mineral commodity supply risk using two criteria: (1) an economic effects assessment that quantified the potential effects of various trade disruption scenarios on the U.S. economy, and (2) an examination of whether the mineral commodity’s U.S. supply chain relied on a sole domestic producer that represented a single point of failure. For the first criterion, postdisruption equilibrium quantities and prices for each mineral commodity were calculated based on their price elasticities of supply and demand and the availability of excess production capacity for each yearlong foreign trade disruption scenario. Subsequently, a nonlinear optimization routine was used with detailed economic input-output tables to estimate the potential economic effects on the U.S. economy of over 1,200 scenarios for 84 mineral commodities. After accounting for the probability of each scenario’s occurrence, the overall results are presented in terms of changes in U.S. gross domestic product (GDP) by individual industry and the economy overall. The results, which ranged from a net decrease in U.S. GDP of nearly $4.5 billion to a net increase of $33 million, largely reflect U.S. import dependency and world production concentration. Using the Jenks natural breaks optimization method, a statistical classification technique, we categorized the mineral commodities into several classes based on this overall risk quantification. Mineral commodities with annualized probability-weighted net decreases in U.S. GDP greater than $2 million were recommended for inclusion on the LCM. If a mineral commodity did not meet the threshold for inclusion on the LCM under the first criterion, its domestic supply chain was examined under the second criterion, which recommended a mineral commodity for inclusion on the LCM if there was only a single domestic producer. Ultimately, the two criteria resulted in the recommendation of the addition of six mineral commodities (in descending risk order, potash, silicon, copper, silver, rhenium, and lead) to and the removal of two mineral commodities (arsenic and tellurium) from the LCM. By using an economic effects assessment, the results of this analysis provide a prioritization that can also be compared directly against other risk analyses and the cost of various risk mitigation strategies.

Open-File Report

A seismic nodal deployment to understand magmatic structure in the vicinity of the Pahala earthquake swarm

In summer-fall 2022, 80 three-component SmartSolo IGU-BD3C-5 nodal seismometers were deployed surrounding the Pāhala seismic swarm on the Island of Hawaiʻi, with the goal of improving seismicity catalogs, and seismic velocity images of the crust and upper mantle in this region. The Pāhala swarm, located south of Mauna Loa and Kīlauea, has been the site of a multiyear sustained swarm of seismicity at depths of ~ 25 – 40 km, with order of magnitude increases in rate in 2015, and then again in 2019. This seismicity is possibly related to the input of magma from the mantle plume below, which may then be subsequently transported to volcanic edifices. However, these processes remain enigmatic, in part due to a lack of precise earthquake locations and seismic velocity models in this region. Here we provide an overview of the deployment, an assessment of the quality of the collected data, and discuss the viability of the dataset for local earthquake relocation, tomography, and teleseismic receiver functions. Through comparisons with proximal permanent broadband and short period instruments, we find that the nodes produce high quality data, particularly at periods shorter than 5 s, although we find, document, and correct discrepancies with the gain and polarities of the instruments. We successfully record signals from teleseismic earthquakes, even at periods longer than 5 s (the corner of the flat response of the nodes). We also record local earthquakes, including details related to source characteristics. This indicates that the data is likely to prove useful for investigations using both local and teleseismic earthquake signals to better understand the connections between the deep and shallow magmatic systems of Hawaiʻi. While this deployment provides a snapshot in time, its success may provide a useful benchmark for future studies as the volcanic systems of Hawaiʻi continue to evolve in the future.

Hawaii

Inset groundwater-flow models for the Cache and Grand Prairie Critical Groundwater Areas, northeastern Arkansas

The water resources in the Mississippi alluvial plain, located in parts of Missouri, Kentucky, Tennessee, Mississippi, Louisiana, and Arkansas, supports a multibillion-dollar agricultural industry that relies heavily on pumping of groundwater for irrigation of crops and aquaculture. The primary source of groundwater for agricultural-related pumping is the Mississippi River Valley alluvial aquifer, which has declined in storage for decades; secondary groundwater sources include the middle Claiborne aquifer and Wilcox aquifer system. Two areas in northeastern Arkansas that lie within the Mississippi alluvial plain, part of the Cache and Grand Prairie regions, have been designated as Critical Groundwater Areas owing to decades of groundwater declines that resulted from past and current water use. The multidisciplinary Mississippi Alluvial Plain project, led by the U.S. Geological Survey, and funded by their Water Availability and Use Science Program, included objectives to develop numerical groundwater models in focus regions, including the part of the Cache and Grand Prairie regions of northeastern Arkansas. Two inset models were developed using the child model capabilities of MODFLOW 6, the U.S. Geological Survey’s Modular Hydrologic Model simulation software. Both models, called the Cache model and Grand Prairie model, simulated the groundwater system and surface-water/groundwater interactions for the Mississippi River Valley alluvial aquifer and underlying Tertiary-age aquifers and confining units to the Midway confining unit. Each model was spatially discretized into 500-meter x 500-meter orthogonal cells on a grid with 5-meter constant-thickness vertical layers that represented the Mississippi River Valley alluvial aquifer and increasing thickness layers for the aquifers and confining units below the alluvial aquifer. The Cache and Grand Prairie models were calibrated with the PEST++ iterative ensemble smoother Version 5 and employed high dimensional parameterization schemes of 13,740 and 30,436 parameters, respectively. The Cache mean absolute residual for groundwater-level observations within each model domain for the priority well was 1.58 meters. Grand Prairie mean absolute residuals for the alluvial aquifer and middle Claiborne aquifer groundwater-level observations were 2.71 and 10.78 meters, respectively. The groundwater budgets for the Cache and Grand Prairie models were characterized by substantial outflows to irrigation wells, which constituted about 52 and 54 percent of all outflows, with the primary source of water to those wells being releases from unconfined aquifer storage.

Arkansas

Updating the crustal fault model for the 2023 National Seismic Hazard Model for Alaska

We present the crustal fault model for Alaska, based on geologic observations, as a primary input for the 2023 revision of the U.S. Geological Survey National Seismic Hazard Model. We update the 2013 Alaska Quaternary fault and fold database to produce a simplified model of 105 fault sections and four fault zone polygons with basic geologic parameters including slip sense and rate. Significant updates include the following: (1) a slip rate of ∼53 mm/year on the Queen Charlotte Fault indicating it accommodates all of the plate boundary motion; (2) quantified slip rates on megathrust splay faults in the southern Prince William Sound region and near Kodiak Island; (3) improved details of structures in the Chugach-St. Elias orogen; (4) revision of the Castle Mountain Fault from right-lateral slip to a predominantly reverse fault; (5) improved Interior Alaska tectonic models that clarify relationships between the Denali, Totschunda, and thrust faults on both sides of the Alaska Range; (6) identified large earthquake sources in the eastern Brooks Range; and (7) omission of the Chatham Strait section of the Denali Fault. The fault model underscores that the collision of the Yakutat microplate is the dominant driver of active crustal faulting in most of Alaska.

Alaska, Yukon

Investigation of land cover within wetland complexes at Dixie Meadows, Churchill County, Nevada, from October 2015 to January 2022

The U.S. Geological Survey investigated land cover at subannual time steps within six wetland areas in Dixie Valley, Churchill County, Nevada, from October 2015 to January 2022. As requested by the U.S. Fish and Wildlife Service, we used aerial photography and satellite remote sensing data to map surface water and other land cover types within the wetland complexes. We identified five land cover classes using the green normalized difference vegetation index (gNDVI) and its inverse relationship to the normalized difference water index (NDWI) within three U.S. Department of Agriculture National Agriculture Imagery Program aerial images (acquired in 2015, 2017, and 2019) and 110 European Space Agency Sentinel-2 satellite images (acquired 2015–2022). The relative wetness of soil conditions within each land cover class is estimated by comparison to previously published observations of relative conductivity measured by 79 field-based sensors within the wetlands from 2019 to 2021. We mapped the areal coverage of the five land cover classes for approximately 385 acres (1,559,000 square meters [m²]) comprising six individual wetland complexes as well as a larger 1,298- acre (5,254,000-m 2 ) area of interest inclusive of the wetland complexes and adjacent landscape. Land cover of open water (Class 5) primarily within ponds at one of the wetland complexes comprised 8,333 m 2 , on average, of the wetland complexes. Land cover of mixed shallow surface water, saturated soil, and vegetation (Class 4) comprised 111,723 m 2 on average of the wetland complexes. Land cover of dense green vegetation canopy cover (Class 3) that often (46 percent of observations) had underlying surface water or saturated soil conditions comprised 592,522 m 2 on average of the wetland complexes. The remaining areas of the wetland complexes not mapped as these three land cover types (Classes 2 and 1) had sparse vegetation or bare soil cover and commonly (greater than or equal to 67 percent of observations) had dry soil conditions. The investigation of land cover detailed in this report could inform future efforts to map land cover more precisely via higher resolution remote sensing or ground-based surveying or could be incorporated with other environmental monitoring data to characterize habitat and hydrology of the wetland complexes at Dixie Meadows.

Nevada

SURF: An automated method for building nonplanar 3D fault models from earthquake hypocenters

Accurately characterizing 3D fault geometry is vital for improving our understanding of earthquake behavior and informing the development of seismic hazard models. Despite their importance, subsurface fault structures tend to be poorly constrained because of limitations in observational data. Improvements to the seismic networks and earthquake detection algorithms have increased the precision and volume of earthquake catalogs, which help illuminate detailed subsurface fault structure and provide the most direct information available about fault geometries at depth. We present a Python package to automate generating 3D fault geometries directly from hypocentral seismicity patterns. This method begins with clustering events based on their spatial density, identifying coherent patterns. Nearby clusters are then merged based on the similarity of their orientations. We fit nonplanar surfaces using support vector regression to balance surface accuracy with minimal deviations from planarity. The fault models are output as quadrilateral meshes at user‐defined resolution. In the process of generating the 3D fault surfaces, we compute the spatial density of seismicity around the surface and the planarity as quantitative metrics of the model outputs. As a proof of concept, we apply this approach to the San Andreas–Calaveras fault junction region and the 2019 Ridgecrest earthquake sequence, both in California, which contain complex subparallel faults well defined at the Earth’s surface and abundant microseismicity. These case studies demonstrate the method’s ability to model complex fault structures, including long continuous fault surfaces, crossing faults, variably dipping segments, and subparallel faults. We test the method on both standard network catalogs and double‐difference relocated catalogs. We find that our seismicity‐based fault model results align with published 3D models that incorporate additional constraints and interpretations ( Plesch et al. , 2020 ; Aagaard and Hirakawa, 2021 ). This workflow provides a low‐user‐input solution for estimating fault geometries at depth from earthquake catalogs.

California

A newly identified creeping strand of the Concord fault, San Francisco Bay Area

The Concord fault constitutes a major branch of the Pacific–North America transform plate boundary in Northern California, bridging the strike‐slip Bartlett Springs ‐ Green Valley Fault system to the north with the Greenville and Calaveras Faults to the south. Like many faults in the San Francisco Bay Area its long‐term slip is partially accommodated by aseismic slip (creep). Although creep has been recognized and monitored on the northern half of the fault for decades, the precise location of the southern half of the fault and its slip rate—whether accommodated seismically or aseismically—has remained enigmatic. How slip transfers between the Concord and Greenville or Calaveras faults to the south remains an outstanding question. New field observations presented here indicate that the active trace of the fault south of downtown Concord is not where previously interpreted and is indeed actively creeping. We report observations of shallow creep continuing >7 km farther south along the Concord fault than previously reported, along a fault strand not previously recognized for most of its length. This is evident as right‐laterally deflected concrete curbs and sidewalk slabs on both sides of every street that crosses the fault at a high angle in southeast Concord and northeast Walnut Creek. We document the magnitude and location of these deflections to estimate accumulated right‐lateral aseismic slip expressed in engineered structures. Offsets of these piercing lines range from 8 to 18 cm, over widths varying from narrow breaks along centimeter‐scale concrete joints to 10‐m‐wide zones of deflection. Significantly, this active trace is ∼400 m west of where the Quaternary active trace has previously been inferred, placing it within—rather than bounding—the built area of suburban Concord. Slip along the fault has already caused infrastructure damage. These results revise our understanding of the southern Concord fault and help constrain its seismic potential.

California

Deterministic physics-based earthquake sequence simulators match empirical ground-motion models and enable extrapolation to data poor regimes: Application to multifault multimechanism ruptures

We use the deterministic earthquake simulator RSQSim to generate complex sequences of ruptures on fault systems used for hazard assessment. We show that the source motions combined with a wave propagation code create surface ground motions that fall within the range of epistemic uncertainties for the Next Generation Attenuation‐West2 set of empirical models. We show the model is well calibrated where there are good data constraints, and has good correspondence in regions with fewer data constraints. We show magnitude, distance, and mechanism dependence all arising naturally from the same underlying friction. The deterministic physics‐based approach provides an opportunity for better understanding the physical origins of ground motions. For example, we find that reduced stress drops in shallow layers relative to constant stress drop with depth lead to peak ground velocities in the near field that better match empirical models. The simulators may also provide better extrapolations into regimes that are poorly empirically constrained by data because physics, rather than surface shaking data parameterizations, is underlying the extrapolations. Having shown the model is credible, we apply it to a problem where observations are lacking. We examine the case of crustal faults above a shallow subduction interface seen to break coseismically in simulations of the New Zealand fault system. These types of events were left out of consideration in the most recent New Zealand national seismic hazard model due to the modeling complexity and lack of observational data to constrain ground‐motion models (GMMs). Here, we show that in the model, by breaking up the coseismic crustal and interface rupturing fault motions into two separate subevents, and then recombining the resulting ground‐motion measures in a square‐root‐of‐sum‐of‐squares incoherent manner, we reproduce well the ground‐motion measures from the full event rupture. This provides a new method for extrapolating GMMs to more complex multifault ruptures.

Seismological Research Letters

Streams, springs, and volcanic lakes for volcano monitoring

Introduction Volcanic unrest can trigger appreciable change to surface waters such as streams, springs, and volcanic lakes. Magma degassing produces gases and soluble salts that are absorbed into groundwater that feeds streams and lakes. As magma ascends, the amount of heat and degassing will increase, and so will any related geochemical and thermal signal. Subsurface magma movement can cause pressurization that alters hydrostatic head and may induce groundwater discharge. Fluid-pressure changes have been linked to distal volcano-tectonic earthquakes (White and McCausland, 2016; Coulon and others, 2017) and phreatic eruptions (for example, Yamaoka and others, 2016). Clearly, changes in groundwater and surface waters are both indicators of unrest and clues to how and where magma is rising toward the surface. Where possible, it is prudent to incorporate real-time hydrologic data into multiparameter monitoring of restless volcanoes. Hydrologic dynamics can also be tracked by changes in groundwater levels that are commonly measured in shallow boreholes (see of Flinders, A.F., Lowenstern, J.B., Coombs, M.L., and Poland, M.P., eds., Recommended capabilities and instrumentation for volcano monitoring in the United States: U.S. Geological Survey Scientific Investigations Report 2024–5062–K, 5 p., https://doi.org/10.3133/sir20245062k. ">chapter K , this volume, on boreholes; Hurwitz and Lowenstern, 2024). Although inferred to be common, relatively few volcano-hydrology anomalies are well documented, and many are essentially anecdotal (Newhall and others, 2001), reflecting the fact that high-resolution time series remain rare. Extreme examples include the 2008 eruption of Nevado del Huila, Colombia, where relatively minor phreatomagmatic eruptions were accompanied by expulsion of as much as 300 million cubic meters of groundwater from fissures high on the volcano (Worni and others, 2011), generating large lahars. Substantial decreases in flow rate from springs about 8 kilometers from the summit of Mayon Volcano, Philippines, have been noted before most eruptions in the 20th century (Newhall and others, 2001). Stream monitoring at Redoubt Volcano in 2009 allowed Werner and others (2012) to recognize that groundwater was unable to absorb (or scrub) the high flux of volcanic gas and that a high CO 2 /SO 2 precursor signal had been evident for 5 months prior to the eruption. A key to better interpreting hydrologic anomalies—or even identifying them—is therefore obtaining adequate baseline data. Most hydrologic monitoring at U.S. volcanoes has been accomplished by intermittent sampling surveys with annual or less frequent sampling (for example, https://hotspringchem.wr.usgs.gov/index.php ). More frequent sampling, however, generally is needed to establish reliable baselines. A recent hydrologic and hydrothermal monitoring experiment at 25 sites and 10 of the 12 level 4 (very high threat) volcanoes in the U.S. portion of the Cascade Range demonstrated that there is sufficient temporal variability in hydrothermal fluxes, even during quiescent periods, that one-time measurements will commonly have limited interpretive value (Crankshaw and others, 2018). Thus, surveys are best augmented with data from streamgages (for example, Evans and others, 2004; Bergfeld and others, 2008). Streamflow (water discharge) data allow measured temperature and specific conductance to be converted to heat and solute mass fluxes, which could be insightful parameters for detecting anomalous activity (McCleskey and others, 2012). At the Yellowstone Caldera, long-term monitoring of river solutes has allowed calculation of the chloride flux, a proxy for heat discharge (Hurwitz and others, 2007; McCleskey and others, 2016) from the subsurface magma. This is readily accomplished because data from streamgages are continuously recorded and archived by the U.S. Geological Survey (USGS) National Water Information System (NWIS) (USGS, 2024). Similar studies on stratovolcanoes or shield volcanoes would be scientifically useful, and yet are logistically challenging, requiring streamgages on numerous radial drainages complemented by either frequent manual sampling or numerous deployments of equipment to measure water temperature and specific conductance as a proxy for water chemistry. Another challenge is that some volcanic areas, especially shield volcanoes, are characterized by near-surface porous rocks and soils, such that surface streams are rare and replaced by distant, dilute large-volume springs with only a trace of any original volcanically sourced water (Manga, 2001; Hurwitz and others, 2021). Volcanic lakes are worthy of special attention for monitoring efforts, as their temperature and composition can provide evidence of increased flux of volatile-rich fluids from below. Quantifying changes in volatile and heat release from magma can be simpler in lakes than for volcanoes with radial drainages and no major lakes. Moreover, volcanic lakes pose a range of hazards themselves, including phreatomagmatic eruptions, debris flows, flank collapse, tsunamis, and toxic gas release (Mastin and Witter, 2000; Delmelle and others, 2015; Manville, 2015; Rouwet and others, 2015)—hazards that have historically been responsible for substantial loss of life at many volcanoes worldwide (Manville, 2015). Catastrophic CO 2 release at Lake Nyos, Cameroon, in 1986 suffocated about 1,750 people and about 3,500 livestock and was probably triggered by a large landslide into the gas-saturated lake (Kling and others, 1987; Evans and others, 1993). Gas-charged springs in Soda Bay within Clear Lake (California) have caused almost a dozen deaths to bathers in the past hundred years (ABC News, 2000). A 2005 example of lake overturn and abundant gas release was documented at Mount Chiginagak in Alaska (Schaefer and others, 2008) but did not result in any human casualties. Although thermally stratified lakes, which promote trapping of exsolved magmatic gas, tend to develop in tropical regions, the phenomenon can also arise where salinity creates meromixis (a condition in which a lake does not mix completely), as occurs in Mono Lake, California (Jellison and Melack, 1993; Jellison and others, 1998). If magma erupts or flows into a lake, the interaction between hot magma and cold water can be explosive (Mastin and others, 2004; Zimanowski and others, 2015) and substantially expand the area affected by the eruption. Another hazard is the breaching of crater rims by landslides triggered by volcanic and (or) seismic activity. Under some circumstances, substantial volumes of water can be displaced, leading to large floods and lahars. Late Holocene lake flooding from Aniakchak Crater in the Alaska Peninsula (Waythomas, 2022) and from Paulina Lake in Newberry Crater, Oregon (Chitwood and Jensen, 2000), caused by the failure of outlet sills, testify to the substantial hazards at lake-filled calderas. Several volcanic systems in the United States host lakes known to receive heat and gas from underlying magma. These lakes vary widely in area, depth, and chemical composition. Lakes are present at level 4 volcanoes, including Crater Lake and Newberry Volcano in Oregon; Yellowstone Caldera in Wyoming; Long Valley Caldera, Clear Lake volcanic field, Medicine Lake, and Salton Buttes in California; and Aniakchak Crater, Mount Katmai, Fisher Caldera, Mount Okmok, and Kaguyak Crater, among others, in Alaska. A water lake was present in Halemaʻumaʻu, the crater of Kīlauea, Hawai‘i (fig. F1), from October 2019 to December 2020. Level 3 volcanoes with lakes include Mono Lake volcanic field (Calif.), Mount Bachelor (Ore.), Ukinrek Maars and Mount Chiginagak (Alaska), and Soda Lake (Nevada). In addition, there are lakes at many levels 1 and 2 volcanoes. In the United States, there are no strongly acidic lakes that receive abundant input of magmatic gas, such as those found at Mount Ruapehu (New Zealand), Ijen and Kelud (Indonesia), and Poás (Costa Rica). Nevertheless, many contain fluids that provide clues to magmatic processes below. Since publication of a previous report on recommended instrumentation for volcano monitoring (Moran and others, 2008), continuous hydrologic monitoring has become increasingly feasible. However, changes in water pressure, temperature, and chemistry remain, in general, poorly studied phenomena at volcanoes (Sparks, 2003; National Academies of Sciences, Engineering, and Medicine, 2017). Recent efforts by the USGS have included the temporary study of Cascade Range volcanoes, which included frequent (15 minute to hourly) temporal sampling of temperature, depth, and conductivity (Crankshaw and others, 2018; Ingebritsen and Evans, 2019). At Yellowstone Caldera, many streamgages have now added thermistors and specific conductance sensors, allowing estimation of time-dependent chloride flux as a proxy for variations in subsurface heat flux (McCleskey and others, 2012, 2016). Efforts to better understand lakes have also accelerated, with bathymetric mapping and sampling carried out at several locations in the United States. Especially thorough work was done at Yellowstone Lake thanks to the Hydrothermal Dynamics of Yellowstone Lake (HD-YLAKE, https://hdylake.org ) project, funded primarily by the National Science Foundation. In addition to geophysical surveys and recovery of cores and other samples, HD-YLAKE investigations included remotely operated vehicle (ROV) investigations of hydrothermal vents on the lake floor (fig. F2). Data collected by the ROV provided a better understanding of the thermal and chemical influx from lake-bottom hydrothermal systems (Sohn and others, 2017). In this chapter, we focus on detecting changes in the chemistry, temperature, discharge, or water levels of streams, springs, and lakes that can be caused by seismicity, volumetric strains, or increases in gas flux associated with ascending magma. There is unavoidable overlap with other chapters of this report. Samples of water and gas can also be obtained in boreholes ( of Flinders, A.F., Lowenstern, J.B., Coombs, M.L., and Poland, M.P., eds., Recommended capabilities and instrumentation for volcano monitoring in the United States: U.S. Geological Survey Scientific Investigations Report 2024–5062–K, 5 p., https://doi.org/10.3133/sir20245062k. ">chapter K , this volume; Hurwitz and Lowenstern, 2024), both shallow and deep. Gas monitoring ( of Flinders, A.F., Lowenstern, J.B., Coombs, M.L., and Poland, M.P., eds., Recommended capabilities and instrumentation for volcano monitoring in the United States: U.S. Geological Survey Scientific Investigations Report 2024–5062–E, 11 p., https://doi.org/10.3133/sir20245062e.">chapter E , this volume; Lewicki and others, 2024) relies in part on samples from springs and wells, particularly where measurable gas plumes are absent. Water acts as a trigger and lubricant for landslides and sediment-rich floods, and so hydrology has obvious relevance for lahar monitoring, as discussed in of Flinders, A.F., Lowenstern, J.B., Coombs, M.L., and Poland, M.P., eds., Recommended capabilities and instrumentation for volcano monitoring in the United States: U.S. Geological Survey Scientific Investigations Report 2024–5062–H, 6 p., https://doi.org/10.3133/sir20245062h. ">chapter H (this volume; Thelen and others, 2024). Shared situational awareness among scientists engaged in geophysical, gas, and hydrologic monitoring will improve overall understanding of the volcanic hazard.

Scientific Investigations Report

The GorDAS Distributed Acoustic Sensing experiment above the Cascadia locked zone and subducted Gorda Slab

The southernmost portion of the Cascadia Subduction zone in Northern California produces high rates of moderate and large earthquakes owing to subduction of the Gorda slab and deformation associated with the Mendocino Triple Junction. Distributed Acoustic Sensing (DAS) is rapidly advancing as a method for detecting earthquakes and imaging crustal structure. We have begun a long-term DAS monitoring experiment on buried telecom fiber in Arcata, California, with the goal of increasing the available recordings of moderate to large earthquakes as well as imaging seismogenic structures. We have recorded over a year's worth of data, including most aftershocks of the 2022 M w 6.4 Ferndale earthquake, though not the mainshock itself. The dataset includes numerous magnitude 3.5 and larger earthquakes including the 2023/01/01 M w 5.4 Rio Dell earthquake. Here we present initial results comparing an earthquake detection algorithm, run in real-time on the processing unit of the interrogator system, with both the ShakeAlert earthquake early warning system as well as a post-processed earthquake catalog developed with deep-learning phase-picker algorithms. The rapid onboard processing of the detector demonstrates the potential utility of DAS-based edge computing for earthquake early warning. We also verify the quality of the strain waveforms both in terms of peak amplitudes and waveform similarity using about five months of nodal seismometer data. These instruments were deployed roughly every 300 m along the ~15km long cable and validate large variations in peak strain over short distances that are seen in the DAS data. All data from time windows surrounding both the local and teleseismic earthquakes are publicly available, which will improve our understanding of both the performance of DAS systems in moderate earthquakes and earthquake hazards associated with the Gorda subduction zone.

California

A ground-motion model derived using a generalized mean rupture distance for large slab interface earthquakes

Source–station distance is a central input to ground‐motion models (GMMs) for predicting seismic shaking. GMM development uses distance metrics including the Joyner–Boore distance, which is the shortest distance from an observation point to the surface projection of the earthquake rupture, and R rup the shortest distance to the rupture in three dimensions. Thompson and Baltay (2018) proposed the generalized mean rupture distance R p to address observed near‐fault ground‐motion saturation. R p accounts for the contribution to the shaking of all parts of the rupture and provides a simple method for incorporating spatially variable slip. They used R p to develop a GMM for shallow crustal earthquakes, assuming uniform slip. Here, we investigate the improvement offered by an R p ‐based GMM for large subduction interface earthquakes by recalibrating the path term for a published GMM ( Parker et al. , 2022 ) using R p derived from distributed slip models (DSMs). Inspection of within‐event and total residuals indicates that the recalibrated model fits the data at least as well as an alternative in which the path term was recalibrated using R rup and the same dataset. Incorporating slip information in its entirety or trimming the DSM geometry is preferable to assuming uniform moment release on a prescribed model fault that extends beyond the actual ruptured area. R p tuned to minimize model uncertainty is closer to the maximum distance to rupture for peak ground velocity (PGV) than acceleration, possibly reflecting the contribution to high‐frequency shaking of slip on local asperities. However, when the moment is concentrated far from stations, PGV is fit adequately with R p closer to R rup ⁠ . Our results suggest incorporating slip‐derived moment release through R p could improve GMMs for slab interface events, especially if the implementation of R p ‐based models is refined using a larger dataset of earthquakes with greater geographic diversity to account for regional ground‐motion variations.

Seismological Research Letters

Analysis of the potential effects of Uzbekistan’s mineral endowment on the critical mineral supply of tungsten

Tungsten appears on the 2018 and 2022 U.S. Geological Survey critical mineral lists in part because of a very high global production concentration in China, which produces almost 83 percent of the world’s mined tungsten. Using known parameters and values from other tungsten mining operations, we created hypothetical scenarios in which three tungsten deposits in Uzbekistan are considered for development. Our results show that all three deposits are likely to be economically viable to develop under 2024 market conditions. If the three studied tungsten deposits were put into production, Uzbekistan could become the third-leading tungsten-producing country in the world and increase world output of tungsten by 2.7 percent. Putting these tungsten deposits in Uzbekistan into production could slightly reduce the tungsten global market concentration, therefore reducing the supply disruption potential for tungsten.

Open-File Report

Application of Hydrologic Simulation Program—FORTRAN (HSPF) as part of an integrated hydrologic model for the Salinas Valley, California

The U.S. Geological Survey (USGS), in cooperation with the Monterey County Water Resources Agency, conducted studies to help evaluate the surface-water and groundwater resources of the Salinas Valley study area, consisting of the entire Salinas River watershed and several smaller, adjacent coastal watersheds draining into Monterey Bay. The Salinas Valley study area is a highly productive agricultural region that depends on the coordinated use of surface water and groundwater to meet demand for irrigation and public water supply. To continue to meet these demands, a better understanding of the historical water balance and the effects of water-resource development on the long-term sustainability of water resources in the Salinas Valley study area is needed.

California