USGS ScienceSearch

SEARCH · USGS Science

Results for “Transactions of the American Fisheries Society”

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

A regional model comparison between MODPATH and MT3D of groundwater travel time distributions

Groundwater quality changes in wells and streams lag behind changes to land use due to groundwater travel times. Two contaminant transport methods were compared to assess differences in their simulated travel time distributions (TTDs) to streams and wells in the Wisconsin Central Sands. MODPATH simulates advective groundwater flow with particle tracking, while MT3D simulates age-mass using a finite difference solution without dispersion to allow for direct comparison of the two methods. MODPATH appropriately simulates groundwater TTDs from the water table to surface discharge but is subject to inaccuracies at weak-sink well cells due to the flow-model grid discretization and imprecise location of well discharge within well cells. MT3D better represents weak-sink well cells since it removes mass in proportion to the prescribed pumping rate, although travel time within well cells is neglected. Conversely, MT3D's treatment of surface water boundary cells is not as accurate as MODPATH because mass should be removed from the water table rather than the full cell volume. MT3D simulations of TTDs can also be confounded by the instantaneous vertical distribution of mass introduced throughout recharge cells instead of at the water table, which initiates mass along deeper flow paths. We evaluated 9 MODPATH and 13 MT3D implementations, generating differences in median travel times of up to 18 years. Both methods have strengths and weaknesses, with MT3D better representing weak-sink well cell behavior and MODPATH better representing surficial recharge and discharge. The effect of these characteristics on simulated TTDs, along with ideas for ameliorating method weaknesses, is discussed.

Wisconsin

Geochemistry of quartz and fluid inclusions associated with gold, tungsten, and antimony in the Stibnite-Yellow Pine district, Idaho

Deposits in the historic Stibnite-Yellow Pine district in east-central Idaho were mined for Ag, Au, Hg, Sb, and W. Ore is hosted in breccia and fault zones along NE-striking faults that transect the Atlanta lobe of the Late Cretaceous Idaho batholith and metamorphosed Neoproterozoic to lower Paleozoic rocks. Geochronological studies have delineated five superimposed and zoned hydrothermal events in the mining district which range in age from Late Cretaceous to the Eocene. They formed the following consecutive assemblages: (1) vein and disseminated pre-ore milky quartz (± calcite), muscovite, and brassy pyrite with minor molybdenite and base metal sulfides; (2) vein and disseminated quartz, calcite, sericite, auriferous-pyrite, and arsenopyrite; (3) potassium feldspar, quartz, dolomite, and calcite with vein and disseminated scheelite ± pyrite and arsenopyrite; (4) Au-Ag veins that contain adularia; and (5) breccia, vein, and disseminated quartz, calcite, and stibnite with minor Ag sulfosalts and distal zones of cinnabar and realgar at Fern and Hermes mines. In this study, the pressure, temperature, and composition of ore-forming fluids were measured in Au-, W-, and Sb-bearing samples collected from drill core and surface exposures. Using optical petrography, scanning electron microscopy cathodoluminescence (SEM-CL) and electron microprobe cathodoluminescence (EMPA-CL) imaging, as well as electron microprobe (EMP) and laser ablation-inductively coupled plasma-mass spectrometry trace element analyses, four major overprinting generations of hydrothermal quartz are defined. Fluid inclusion assemblages (FIAs) in each generation of quartz were characterized by microthermometry and laser Raman spectroscopy. Quartz and fluid inclusion paragenesis are as follows: (Qi) altered and overprinted Qi with no coeval FIAs that formed before or during event 1; (Q1) early veins as well as disseminated quartz with dull gray SEM-CL and H 2 O-NaCl-CO 2 FIAs that formed at minimum temperatures and pressures of 281° to 358°C and 1346 to 1789 bar during event 2 (Au); (Q2) vein and disseminated quartz, often euhedral, with H 2 O-NaCl ± CO 2 ± CH 4 FIAs that formed at minimum temperatures and pressures of 183° to 213°C and 72 to 602 bar during event 3 (W); and (Q3) euhedral and disseminated quartz with H 2 O-NaCl FIAs that formed at minimum temperature and pressures of 156° to 161°C and 5 to 6 bar during event 5 (Sb). Stibnite hosts H 2 O-NaCl FIAs trapped at minimum temperatures and pressures of 150° to 189°C and 5 to 12 bar that are similar to those of Q3-hosted FIAs. The quartz paragenetic history in the district follows a progressive temporal path from Tienriched quartz deposited by a relatively deep and higher temperature aqueous-carbonic fluid to a Li, Al, and Sb enriched quartz precipitated from a lower temperature aqueous fluid. Gases in fluid inclusion extracts from ore and gangue minerals were analyzed by mass spectrometry. Noble gas isotope analysis of fluid inclusion extracts revealed three signatures. Pre-ore Qi has R/R A (the isotopic ratio of the sample 3 He/ 4 He to atmospheric 3 He/ 4 He (1.384 × 10 −6 )) values that reflect the composition of the crustal melt-derived Atlanta lobe of the Idaho batholith. Quartz (Q1) and auriferous arsenian pyrite of hydrothermal event 2 have 20 Ne/ 4 He and R/R A values that plot along the air-crust mixing line, and these gases were likely introduced by metamorphic fluids derived from metasedimentary rocks. Stibnite with associated quartz (Q3 ± Q2) have 20 Ne/ 4 He and R/R A values that extend from the air-crust mixing line to the air-mantle mixing line, showing a shift from a crustal signature to an overprinting mantle signature. Late calcite has 20 Ne/ 4 He and R/R A values that extend from air toward the mantle end member, which suggests that magmatic volatiles condensed into meteoric groundwater late in the paragenesis. A mantle volatile source is also supported by high 20 Ne/ 22 Ne ratios. The Stibnite-Yellow Pine district is similar to other Au-W-Sb districts around the world in ore fluid chemistries, with late-stage stibnite being a product of cooling. The fluid pressure-temperature-composition (P-T-X) estimates and the crustal source of He during hydrothermal event 2 gold mineralization are characteristic of orogenic Au deposits, and the observed inclusion reequilibration textures indicate that trapped fluid inclusions were deformed during transitions between the brittle-ductile environment. Scheelite was deposited during event 3 with Q2 quartz and is overgrown by Q3 quartz and stibnite. Our interpretation of the quartz texture, chemistry, and fluid inclusion P-T-X data suggest that Au, W, and Sb mineralization occurred as hydrothermal fluids evolved from H 2 O-NaCl-CO 2 at >300°C and depths of ~4 to 6 km near the ductile to brittle transition to H 2 O-NaCl at ~150°C and epithermal depths as the Idaho batholith was exhumed. Geochronology indicates that from the Late Cretaceous to early Eocene, overprinting hydrothermal events formed deep orogenic Au, intermediate W, and shallow epithermal Au-Ag and Sb. Without comprehensive geochronologic constraints, our fluid inclusion data and textural observations could be interpreted as a single evolving hydrothermal system, as has been done in many other mining districts with similar mineralization.

Idaho

Constraining large magnitude event source and path effects using ground motion simulations

The purpose of this study is to use ground motion simulations to investigate ways in which source and path effects for large magnitude events can be represented in non-ergodic GMMs. While we initially developed computation techniques using CyberShake simulations, the range of magnitudes and source-site combinations is not adequate to replicate what is observed empirically. We therefore designed a new ground motion simulation study, which includes earthquakes with a large range of magnitudes distributed uniformly on a fault plane, and sites covering a large range of rupture distances and azimuths. After running a large suite of simulations (M4-M7), we then develop a non-ergodic GMM with the simulation data. We find that the within-site residuals are dominated by the radiation pattern, rupture directivity, and slip patterns. Next, we modify an existing rupture directivity model to fit and remove the observed radiation pattern and rupture directivity from the residuals. We also minimize the contributions of slip patterns by averaging the within-site residuals among multiple source realizations. Finally, after removing the source effects from the within-site residuals, we compare the path effects computed with different magnitude groups using two approaches. The first approach only considers the small events that have the same shortest path to a site with the large events, while the second approach considers all small events on the fault plane. The results indicate that the path effects of large events cannot be satisfactorily approximated with that of small events using either approach.

Conference Paper

Mineral commodity summaries 2025

Introduction Each mineral commodity chapter of the 2025 edition of the U.S. Geological Survey (USGS) Mineral Commodity Summaries (MCS) includes information on events, trends, and issues for each mineral commodity as well as discussions and tabular presentations on domestic industry structure, Government programs, tariffs, 5-year salient statistics, and world production, reserves, and resources. The MCS is the earliest comprehensive source of 2024 mineral production data for the world. More than 90 individual minerals and materials are covered by two-page synopses. Abbreviations and units of measure and definitions of selected terms used in the report are in Appendix A and Appendix B, respectively. Reserves and resources information is in Appendix C, which includes “Part A—Resource and Reserve Classification for Minerals” and “Part B—Sources of Reserves Data.” A directory of USGS minerals information country specialists and their responsibilities is in Appendix D. The USGS continually strives to improve the value of its publications to users. Constructive comments and suggestions by readers of the MCS 2025 are welcomed.

Mineral Commodity Summaries

Mineral commodity summaries 2026

Introduction Each mineral commodity chapter of the 2026 edition of the U.S. Geological Survey (USGS) Mineral Commodity Summaries (MCS) includes information on events, trends, and issues for each mineral commodity as well as discussions and tabular presentations on domestic industry structure, Government programs, tariffs, 5-year salient statistics, and world production, reserves, and resources. The MCS is the earliest comprehensive source of 2025 mineral production data for the world. More than 90 individual minerals and materials are covered by two-page synopses. Abbreviations and units of measure and definitions of selected terms used in the report are in Appendix A and Appendix B, respectively. Reserves and resources information is in Appendix C, which includes “Part A—Resource and Reserve Classification for Minerals” and “Part B—Sources of Reserves Data.” A directory of USGS minerals information country specialists and their responsibilities is in Appendix D. The USGS continually strives to improve the value of its publications to users. Constructive comments and suggestions by readers of the 2026 MCS are welcomed.

Mineral Commodity Summaries

Ungulate migrations of the Western United States, volume 4

Broadly distributed across the Western United States, ungulates (hooved mammals) play an important role in ecosystem function by affecting vegetation communities and forming the prey base for large carnivores. Additionally, ungulates provide economic benefits to regional communities through tourism and hunting and hold cultural significance for many Tribal communities. Many ungulates migrate seasonally between distinct summer and winter ranges to take advantage of spatially and temporally variable food sources and avoid threats such as predators and deep snow. Increasingly, these migrations are threatened by the growing human footprint and associated subdivisions, energy development, and increased traffic volume. Efforts to study ungulate populations and conserve their migrations received support in recent years from the U.S. Department of the Interior Secretarial Order No. 3362, which provided Federal support for enhancing habitat quality for ungulates across the Western States. In response to Secretarial Order No. 3362, the U.S. Geological Survey (USGS) established the Corridor Mapping Team, a collaboration among USGS and participating State and Federal wildlife management agencies and numerous Tribal Nations. Together, the Corridor Mapping Team maps ungulate migrations throughout the Western United States in the USGS “Ungulate Migrations of the Western United States” report series. This report (volume 4) details migrations and seasonal ranges from 31 new herds throughout nine Western States. Additionally, this report includes updates to two herds published in previous reports. Including this report, the report series has provided the mapped migrations and seasonal ranges of 182 unique herds and has provided a map-based inventory of the documented ungulate migrations across the Western United States for biologists, managers, policy makers, and conservation practitioners. This report also discusses how the mapping efforts associated with the Corridor Mapping Team can be used to guide management and policy regarding renewable energy development and ungulate disease, specifically chronic wasting disease, in the Western United States.

Arizona, California, Colorado, New Mexico, Nevada,

Magnetic storms and geoelectric hazards

Magnetic storms induce geoelectric fields at Earth's surface that can interfere with grounded long-line systems. The September 1859 storm disrupted global telegraph operations, the March 1989 storm caused a blackout in Canada and interfered with electric-power-transmission systems in the United States, and other storms have had related impacts. The geographic and temporal dependence of geoelectric fields are functions of both geomagnetic variation and local surface impedance, which differ considerably across different geological regions. These dependencies can be mapped across the contiguous United States by combining magnetotelluric impedance tensors with ground magnetometer time series. This review illustrates such mapping for the 1989 storm and shows that power-system interference was experienced where surface impedance is high, and when and where geoelectric fields were intense. Statistical analyses indicate that storms comparable to that of March 1989 occur roughly once every four solar cycles. Ongoing developments in numerical modeling and real-time monitoring are anticipated to enable prediction of geoelectric hazards. ▪ Magnetic storms can induced electric fields in the solid Earth that interfere with electric-power-transmission systems. ▪ Geoelectric hazards depend on the storm-time geomagnetic disturbance and the electrical conductivity structure of Earth. ▪ Historically, impacts on telecommunication and power-transmission systems in the United States have been concentrated in the East and Midwest. ▪ The future occurrence of a magnetic superstorm could cause widespread disruption of electric-power-transmission systems.

Annual Review of Earth and Planetary Sciences

Updating and recalibrating the integrated Santa Rosa Plain Hydrologic Model to assess stream depletion and to simulate future climate and management scenarios in Santa Rosa, Sonoma County, California

The Santa Rosa Plain Hydrologic Model (SRPHM) was developed and published in 2014 through a collaboration between the U.S. Geological Survey (USGS) and Sonoma Water to analyze the hydrologic system in the Santa Rosa Plain watershed, help meet the increasing demand for fresh water, and prepare for future uncertainties in water resources. The original model simulated hydrological conditions and water use from water years 1975 to 2010. Recently (2023), the USGS, in cooperation with Sonoma Water and the California State Water Resources Control Board, updated the SRPHM model to extend its simulation period to the end of the 2018 calendar year, incorporate new estimates of rural and agricultural water use, and use efficient input format for climate variables. The updated model was recalibrated, and evaluation of the new model calibration is included in this report. This report presents the results of comparing the hydraulic heads, streamflow, and groundwater budget simulated by the updated model with those generated by the original model and observed data. The main difference in the simulated budget between the original and updated SRPHM is the estimates of agricultural pumping, rural domestic pumping, and return flow generated from rural water use that was not simulated in the original model. The revised agricultural pumping is simulated using the agricultural package, which constrains pumping to available groundwater. The use of the agricultural package leads to a more realistic estimation of agricultural water use, with revised agricultural pumping being one-third less than that in the original model. The revised rural pumping is about half of the pumping in the original model because of using detailed parcel data to estimate population density in rural areas instead of coarse census tracts. Overall, average total inflows for water years 2006–10 simulated by the updated model were about 2 percent less than the original model, and the average total updated outflows were nearly 5 percent less than the original model. The updated model was then used to generate stream depletion maps, simulate climate change scenarios during 2019–99, and simulate water rights allocation using the Model for Decision Support in Integrated River Basin Management (MODSIM). The results from simulating eight future climate scenarios indicated either an increase in groundwater storage or no significant change in the next 80 years, along with an increase in recharge, an increase in actual evapotranspiration in six out of eight climate projections, and an increase in surface runoff. The increases in the simulated future groundwater storage, recharge, evapotranspiration, and runoff in most climate projections are mainly driven by the projected increase in precipitation in most of the future climate scenarios. The updated model also was used to test a pilot case study demonstrating water-resource allocation among different users with different water rights using the integrated MODSIM-Groundwater and Surface-Water Flow Model (GSFLOW) platform. The updated SRPHM serves as a valuable tool for analyzing historical and future hydrologic conditions in the Santa Rosa Plain watershed and preparing for future uncertainties.

California

Monthly variability and possible sources of nitrate in ground water beneath mixed agricultural land use, Suwannee and Lafayette Counties, Florida

In an area of mixed agricultural land use in Suwannee and Lafayette Counties of northern Florida, water samples were collected monthly from 14 wells tapping the Upper Floridan aquifer during July 1998 through June 1999 to assess hydrologic and land-use factors affecting the variability in nitrate concentrations in ground water. Unusually high amounts of rainfall in September and October 1998 (43.5 centimeters total for both months) resulted in an increase in water levels in all wells in October 1998. This was followed by unusually low amounts of rainfall during November 1998 through May 1999, when rainfall was 40.7 centimeters below 30-year mean monthly values. The presence of karst features (sinkholes, springs, solution conduits) and the highly permeable sands that overlie the Upper Floridan aquifer provide for rapid movement of water containing elevated nitrate concentrations to the aquifer. Nitrate was the dominant form of nitrogen in ground water collected at all sites and nitrate concentrations ranged from less than 0.02 to 22 milligrams per liter (mg/L), as nitrogen. Water samples from most wells showed substantial monthly or seasonal fluctuations in nitrate concentrations. Generally, water samples from wells with nitrate concentrations higher than 10 mg/L showed the greatest amount of monthly fluctuation. For example, water samples from six of eight wells had monthly nitrate concentrations that varied by at least 5 mg/L during the study period. Water from most wells with lower nitrate concentrations (less than 6 mg/L) also showed large monthly fluctuations. For instance, nitrate concentrations in water from four sites showed monthly variations of more than 50 percent. Large fluctuations in nitrate concentrations likely result from seasonal agricultural practices (fertilizer application and animal waste spreading) at a particular site. For example, an increase in nitrate concentrations observed in water samples from seven sites in February or March 1999 most likely results from application of synthetic fertilizers during the late winter months. Lower nitrate concentrations were detected in water samples from five of eight wells sampled during high-flow conditions for the Suwannee River in March 1998 compared to low-flow conditions in November 1998. Evidence for reduction of nitrate due to denitrification reactions was observed at one site (AC-1), as indicated by elevated concentrations of nitrogen gas and a corresponding increase in nitrogen isotope (d15N-NO3) values with a decrease in nitrate concentrations. Denitrification is unlikely at other sites based on the presence of dissolved oxygen concentrations greater than 2 mg/L in ground water and no observed trend between nitrate concentrations and values d15N-NO3 values. Nitrate was the dominant nitrogen species in most monthly rainfall samples; however, ammonium concentrations were similar or greater than nitrate during November and December 1998. During February through May 1999, both nitrate and ammonium concentrations were substantially higher in monthly rainfall samples collected at the study area compared to mean monthly concentrations at the Bradford Forest site located east of the study area, which is part of the National Atmospheric Deposition Program/National Trends Network. Also, higher nitrogen deposition rates in the study area compared to those at Bradford Forest could indicate that substantial amounts of ammonia are volatilized from fertilizers and animal wastes, released to the atmosphere, and incorporated as nitrate and ammonium in rainfall deposited in the middle Suwannee River Basin. Ground-water samples from most sites had d15N-NO3 values that indicated a mixture of inorganic and organic sources of nitrogen, which corresponded to multiple land uses where both synthetic fertilizers and manure are used on fields near these sites. Distinct d15N-NO3 signatures, however, were observed at some sites. For example, water samples from areas of row-crop farming as the dominant land use had d15N-NO3 values less than 4 per mil, indicating an inorganic nitrogen source such as synthetic fertilizer. In contrast, d15N-NO3 values greater than 9 per mil were found in water samples from three sites where manure from dairy and or poultry operations is spread on fields throughout the year.

Florida

Creation of aerial camera system boresight installations in Bangor, Maine; Monroe, Louisiana; and Bassville Park, Florida

Aerial camera systems rely on boresights for defining internal and external metrics. The U.S. Fish and Wildlife Service’s Migratory Bird Surveys Branch required boresights to calibrate their camera systems for creating georeferenced imagery that can then be used in geographic information systems. The U.S. Geological Survey created high-accuracy boresight installations in Bangor, Maine, in July 2024; Monroe, Louisiana, in September 2024; and Bassville Park, Florida, in January 2025 to provide this metric-defining capacity to the U.S. Fish and Wildlife Service. This U.S. Geological Survey Data Report details the methods and results of the boresight installation process and provides camera operators with the data needed to use each of the installations.

Florida, Louisiana, Maine

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

Method for estimating water use and interbasin transfers of freshwater and wastewater in an urbanized basin

Techniques for management of drainage basins that use water budgets to balance available water resources with actual or anticipated water use require accurate and precise estimates of basin withdrawals, interbasin transfers of freshwater, unaccounted-for use, water use, consumptive use, inflow and infiltration, basin return flow, and interbasin transfers of wastewater. Frequently, interbasin transfers of freshwater and wastewater are not included in basin water budgets because they occur within public water-delivery and wastewater-collection systems. A new 10-step method was developed to improve estimates of inflow and infiltration and interbasin transfers using readily available statewide data. The accuracy and precision of water-use estimates determined by this method are improved through careful application of coefficients for small users and the use of metered values for large users. The method was developed and tested with data for the Ten Mile River Basin in southeastern Massachusetts. This report uses examples from the basin to illustrate each step of the method.

Massachusetts

Ground magnetic and VLF studies at Midnite uranium mine, Stevens County, Washington

Ground magnetic and E-mode VLF measurements were made at more than 200 stations at the Midnite uranium Mine, Stevens County, Washington. The equipment was found to be light and foolproof, and each set of data was collected in one day, exclusive of surveying and marking the station points. The maps of contoured magnetic and VLF data resembled geologic maps of the area, showing many common features, as well as a few unique features on each map. Each showed some evidence of a small normal fault cutting the region; taken together, the evidence for this fault was judged to be much stronger, so that further structural details were inferred. Because of soil cover and inadequate outcrops, the presence of this fault had hitherto been unsuspected. A target area for further uranium mineralization north of the fault was indicated by this quick, inexpensive geophysical survey. The drilling program which Midnite Mine management had underway at that time had already outlined a mineralized zone in the indicated area.

Washington

Reaction kinetics and accelerant effects of sulfides in early mature hydrocarbon generation using hydrous pyrolysis

Hydrocarbon generation in organic-rich sediments is influenced by the molecular organic composition and relative abundance of associated minerals. Certain mineral-derived elements act as catalysts and reaction intermediaries, facilitating early-stage hydrocarbon formation in potential source rocks. This study investigated the role of sulfur contributed from pyrite as an accelerant in thermal reaction, focusing on its effects on early maturation and consequent hydrocarbon generation from gilsonite (low-sulfur solid petroleum). Hydrous pyrolysis (HP) experiments were conducted on mixtures of gilsonite and pyrite in varying ratios (1:0.1, 1:0.5, 1:1, 1:2, and 1:10 w/w gilsonite:pyrite) at 320, 350, and 370 °C for 72 h. Untreated and thermally altered residues were analyzed using solid bitumen reflectance (BR o , %), total organic carbon (TOC) content, programmed temperature pyrolysis, scanning electron microscopy with energy-dispersive spectroscopy (SEM-EDS), and X-ray diffraction (XRD) to evaluate the potential accelerant role of pyritic sulfur in hydrocarbon formation. The results show HP residues at 320 and 350 °C with greater pyrite concentrations had higher BR o , while reflectance values were similar in the 370 °C residues, regardless of pyrite concentration, suggesting enhanced reaction at lower thermal conditions. Increasing pyrite content systematically decreased hydrogen index (HI) values while increasing the transformation ratio (TR) and production index (PI), indicating enhanced conversion of organic matter to hydrocarbons with increasing pyrite concentrations. Gas yields increased with pyrite addition, particularly at 350 °C, confirming secondary cracking effects. However, gas production stabilized or declined at higher pyrite loadings (1:10), suggesting alternative reaction pathways such as coke formation. Our data indicate the presence of pyrite lowers the activation energy for thermal cracking, shifting peak experimental hydrocarbon generation temperatures downward by 20–30 °C, with the most pronounced accelerant effects observed at moderate pyrite concentrations (1:0.5 and 1:1). The thermodynamic framework reveals that pyrite stability is influenced by experimental conditions, with pyrrhotite formation favored in the presence of gilsonite due to reduced oxygen fugacity. Pyrite transformation to pyrrhotite, as observed through XRD, SEM-EDS, and predicted by thermodynamic data, further supports the accelerant role of S, as pyrrhotite exhibits a higher hydrogen transfer potential, promoting early oil generation. These findings highlight the importance of pyrite in modulating hydrocarbon generation pathways in organic-rich systems.

Journal of Analytical and Applied Pyrolysis

Uncertainty quantification of geologic energy storage in depleted gas reservoirs using material balance equations embedded in a hierarchical errors-in-variables model

The storage potential for gas in geologic settings, such as depleted hydrocarbon reservoirs and solution-mined salt caverns, is becoming salient to future energy infrastructure planning. Technologies such as carbon capture, utilization, and storage, carbon dioxide-enhanced oil recovery, and natural gas and hydrogen storage help meet growing energy demands, reduce carbon emissions to meet climate goals, and provide energy security amid geopolitical uncertainties. Therefore, estimates of underground gas storage capacity could be useful for efficiently navigating the energy transitions. Material balance is a fundamental method in reservoir engineering for estimating original gas in place and potential storage capacity at the scale necessary for national assessments of subsurface pore space resources. However, the deterministic method cannot accommodate multiple data sources or quantify uncertainty in predictions. In this study, a method that embeds material balance equations within a hierarchical errors-in-variables model is proposed which allows the estimation of the distributions of reservoir properties needed for assessments. Uncertainties associated with these reservoir properties have traditionally been expert-elicited, whereas the uncertainty estimates from the proposed models are data-driven. Capacity and uncertainty estimates can be used in a probabilistic resource assessment, supplementing information traditionally used by assessors or even replacing this expert elicitation step when data are unavailable. Various regression models are compared in a case study of the Michigan Basin, a large contributor to the United States’ current natural gas storage capacity. In particular, errors-in-variables models help ameliorate regression dilution and can quantify uncertainty in predictions of pressure in addition to storage capacity. Overfitting is addressed by quantifying generalization error and model averaging in simple and stratified cross-validation against reported working gas capacity, representing the varying quality and quantity of available data. Incorporating a statistical framework into existing numerical methods in reservoir engineering can improve the quality of estimation, and in particular, this method brings rigor to uncertainty quantification as part of a larger effort by the U.S. Geological Survey to assess domestic energy gas storage resources in depleted hydrocarbon reservoirs.

Mathematical Geosciences

Aftershock forecasting

Aftershocks can compound the impacts of a major earthquake, disrupting recovery efforts and potentially further damaging weakened buildings and infrastructure. Forecasts of the probability of aftershocks can therefore aid decision-making during earthquake response and recovery. Several countries issue authoritative aftershock forecasts. Most aftershock forecasts are based on simple statistical models that were first developed in the 1980s and remain the best available models. We review these statistical models, and the wide-ranging research to advance aftershock forecasting through better statistical, physical, and machine learning methods. Physics-based forecasts based on mainshock stress changes can sometimes match the statistical models in testing, but don’t yet outperform them. Physical models are also hampered by unsolved problems such as the mechanics of dynamic triggering and the influence of background conditions. Initial work on machine learning forecasts shows promise, and new machine learning earthquake catalogs provide an opportunity to advance all types of aftershock forecasts.

Annual Review of Earth and Planetary Sciences

Determination of baseline water-quality conditions using continuous specific conductance and discrete measurements of total dissolved solids and selected major ions, Arkansas River Basin, Colorado

From 2022 through 2024, the U.S. Geological Survey (USGS) in cooperation with the Regional Resources Planning Group (Southeastern Colorado Water Conservancy District, Colorado Springs Utilities, Pueblo Water, Lower Arkansas Valley Water Conservancy District, Aurora Water, and Upper Arkansas Water Conservancy District) established a basinwide water-quality monitoring network in the Arkansas River Basin in Colorado. The study leveraged ongoing continuous measurements of streamflow and specific conductance (SC) at 13 sites along the main stem Arkansas River and Fountain Creek from the USGS 07108600 Arkansas River at Granite, Colo., (Ark Granite) site to the USGS 07137500 Arkansas River near Coolidge, Kansas, (Ark Coolidge) site. Baseline water-quality conditions were established using data from October 1989 through September 2022 to help identify changes in water quality that could result from changes in land and water use. Total dissolved solids (TDS) concentration was the main indicator of general water quality in this study because it can be accurately estimated from SC and is sensitive to changes in streamflow, water operations, and source contributions. Median TDS concentrations in Arkansas River samples (October 1989 through December 2023) ranged from 75 milligrams per liter (mg/L) at Ark Granite near the headwaters to 3,660 mg/L at Ark Coolidge near the Colorado–Kansas State line. Distinct seasonal differences in median streamflow and TDS were apparent during the winter and summer months. Seasonal Kendall tests showed no significant trends in monthly mean streamflow (October 1995 through September 2022) at Arkansas River sites upstream from Pueblo Reservoir. Upward trends in TDS concentrations and loads occurred only at Ark Granite. Downward trends in streamflow, TDS concentrations, and TDS loads occurred at Arkansas River sites from USGS 07109500 Arkansas River near Avondale, Colo., to Ark Coolidge. Downward trends in streamflow and SC (surrogate for TDS) were more prevalent during winter base flow months and may indicate changes in groundwater quantity, and possibly groundwater quality, potentially affecting the river water-quality condition.

Colorado

The 3D Elevation Program—Supporting Rhode Island’s economy

Introduction High-resolution elevation data are critical to applications of landscape modeling and planning, both of which have a significant effect on Rhode Island’s economy. In these and other enterprises, program managers, while aiming to strike a balance between accuracy and cost, strive to obtain the best available elevation data to help them address a range of issues. Programs focused on climate change, environmental management, transportation design and asset management, aviation navigation and safety, riverine ecosystem management, wildlife habitat characterization and management, shellfish aquaculture, and the management and mapping of forests, parks and recreation areas, soils, wetlands, and impervious surfaces are also among the critical applications that meet the State’s management needs and depend on light detection and ranging (lidar) data that provide a highly detailed three-dimensional (3D) model of the Earth’s surface and aboveground features. The 3D Elevation Program (3DEP) is managed by the U.S. Geological Survey (USGS) in partnership with Federal, State, Tribal, U.S. territorial, and local agencies to acquire consistent lidar coverage at quality level 2 or better to meet the many needs of the Nation and Rhode Island. The status of available and in-progress 3DEP baseline lidar data in Rhode Island is shown in figure 1. 3DEP baseline lidar data include quality level 2 or better, 1-meter or better digital elevation models, and lidar point clouds, and must meet the Lidar Base Specification version 1.2 ( https://www.usgs.gov/3dep/lidarspec ) or newer requirements. The National Enhanced Elevation Assessment identified user requirements and conservatively estimated that availability of lidar data would result in at least $178,560 in new benefits annually to the State. The top 10 Rhode Island business uses for 3D elevation data, which are based on the estimated annual conservative benefits of 3DEP, are shown in table 2.

Rhode Island