USGS Science⌕ Search

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 90Linked to original sources
USGS1990

Preconditioned Conjugate-Gradient 2 (PCG2), a computer program for solving ground-water flow equations

This report documents PCG2: a numerical code to be used with the U.S. Geological Survey modular three-dimensional, finite-difference, ground-water flow model. PCG2 uses the preconditioned conjugate-gradient method to solve the equations produced by the model for hydraulic head. Linear or nonlinear flow conditions may be simulated. PCG2 includes two reconditioning options: modified incomplete Cholesky preconditioning, which is efficient on scalar computers; and polynomial preconditioning, which requires less computer storage and, with modifications that depend on the computer used, is most efficient on vector computers. Convergence of the solver is determined using both head-change and residual criteria. Nonlinear problems are solved using Picard iterations. This documentation provides a description of the preconditioned conjugate gradient method and the two preconditioners, detailed instructions for linking PCG2 to the modular model, sample data inputs, a brief description of PCG2, and a FORTRAN listing.

Water-Resources Investigations Report↗
USGS2024

Potential climate and human water-use effects on water-quality trends in a semiarid, western U.S. watershed: Fountain Creek, Colorado, USA

Nutrients, total dissolved solids (TDS), and trace elements affect the suitability of water for human and natural needs. Here, trends in such water-quality constituents are analyzed for 1999–2022 for eight nested monitoring sites in the 24,000 km 2 Fountain Creek watershed in Colorado, USA, by using the weighted regressions on time, discharge, and season (WRTDS) methodology. Fountain Creek shares characteristics with other western U.S. watersheds: (1) an expanding but more water-efficient population, (2) a heavy reliance on imported water, (3) a semiarid climate trending towards warmer and drier conditions, and (4) shifts of water from agricultural to municipal uses. The WRTDS analysis found both upward and downward trends in the concentrations of nutrients that reflected possible shifts in effluent management, instream uptake, and water conservation by a watershed population that grew by about 40%. Selenium, other trace elements, and TDS can pose water-quality challenges downstream and their concentrations were found to have a downwards trend. Those trends could be driven by either a warming and drying of the local climate or decreased agricultural irrigation, as both would reduce recharge and subsequent mobilization from natural geologic sources via groundwater discharge. The patterns illustrate how changes in climate and water use may have affected water quality in Fountain Creek and demonstrate the patterns to look for in other western watersheds.

Colorado↗
USGS2026

Special issue “The 2024 M7.6 Noto Peninsula earthquake and seismic swarm”

The 2024 Mw7.5 (M JMA 7.6) Noto Peninsula earthquake struck the Noto Peninsula, central Japan on January 1, 2024, causing widespread damage. This Mw7.5 mainshock, together with the prolonged seismic swarm that persisted in the northeastern Noto Peninsula since the end of 2020, of which the largest event was the 2023 M JMA 6.5 earthquake, constitutes one of the most significant seismic sequences in Japan in recent decades. The sequence provides a rare opportunity to investigate the complex interplay among swarm activity, fluid migration, fault geometry and distribution, and the preparatory conditions required for a large inland crustal earthquake. Motivated by this scientific context, this special issue brings together studies examining the sequence from diverse disciplinary perspectives, incorporating seismic, geodetic, and tsunami-related observations and modelling.

Noto Peninsula↗
USGS2026

Cnidarian–algal partnerships structure bacterial communities during strobilation in Cassiopea xamachana

Cnidarian–algal (Symbiodiniaceae) symbioses rely on complex interactions among the cnidarian host, algal symbionts, and associated bacterial communities. In the upside-down jellyfish Cassiopea xamachana , the polyp-to-medusa transition (strobilation) requires the establishment of symbiosis with Symbiodiniaceae algal partners, yet bacterial community dynamics during this developmental process remain unknown. Here, we experimentally induced symbiosis in aposymbiotic polyps using four algal treatments: xenic Symbiodinium microadriaticum (native symbiont), xenic Breviolum minutum , antibiotic-treated B. minutum , and a photosynthetically impaired B. minutum mutant. We combined 16S rRNA gene sequencing with measurements of photosynthetic efficiency, asexual budding, and algal surface N-glycan profiles to characterize holobiont assembly during symbiosis onset and strobilation. Algal treatment structured bacterial communities in both algal cultures and polyp tissues. Our analyses identified a set of amplicon sequence variants that consistently distinguished strobilating polyps from non-strobilating aposymbiotic and mutant polyps, in addition to potential bacterial biomarkers associated with successful metamorphosis. Strobilation was associated with the enrichment of bacterial communities putatively involved in sulfur and nitrogen cycling, whereas non-strobilating aposymbiotic and mutant polyps were characterized by opportunistic bacteria and increased community variability. Together, these results reveal coordinated changes in algal physiology, surface glycan profiles, and bacterial community structure associated with successful strobilation in C. xamachana and support a model in which tripartite host–alga–bacteria interactions influence cnidarian life stage transitions.

ISME Communications↗
USGS2024

Groundwater quality near the Placerita Oil Field, California, 2018

Groundwater-quality data and potential fluid-migration pathways near the Placerita Oil Field in Los Angeles County, California, were examined by the U.S. Geological Survey to determine if oil-field fluids (water and gas from oil-producing and non-producing zones) have mixed with groundwater resources. Six of the 13 new groundwater samples collected for this study contained petroleum hydrocarbons, thermogenic gas, inorganic chemical signatures, and (or) isotopic values consistent with potential mixing with fluids from hydrocarbon-bearing formations. For historical groundwater samples, benzene was the most detected petroleum hydrocarbon. The historical groundwater samples with a benzene concentration greater than 0.5 micrograms per liter were from environmental monitoring wells at industrial or commercial facilities unrelated to oil and gas development that, in many cases, have identified soil or groundwater contamination and were not typically analyzed for other constituents that could provide additional lines of evidence for potential mixing with oil-field fluids. Methane was not detected in any of the 12 historical samples with a reported measurement. Reviewing historical data revealed factors that could potentially adversely affect groundwater quality in the study area. These factors include modified hydraulic gradients caused by large volumes of water extracted from the main production area and reinjected downgradient into nonproducing zones, well-barrier failures in wells constructed in the northern part of the oil field before the 1970s, well-barrier failures in produced-water disposal wells downgradient from the main production area, and naturally occurring hydrocarbons at shallow intervals. The groundwater samples most geochemically similar to samples from hydrocarbon-bearing formations were in areas where hydrocarbons are naturally occurring at shallow intervals and where oil development is at shallow depths. Additional data for hydraulic heads, water quality, and formation temperatures at multiple depths in areas with large injection volumes and well-integrity issues are needed to evaluate whether those factors have contributed to mixing between fluids from oil-producing or injection formations and groundwater resources.

California↗
USGS2025

Reservoir thermal energy storage pre-assessment for the United States

Storing thermal energy underground for later use in electricity production or direct-use heating/cooling is a promising, viable, and economical green energy option. Reservoir thermal energy storage (RTES) is one such option, which stores energy in underutilized permeable strata with low ambient groundwater flow rates and more geochemically evolved (e.g. brackish/saline) waters relative to overlying principal aquifer systems. The U.S. Geological Survey has begun assessing RTES potential nationally by focusing on five generalized geologic regions (Basin and Range, Coastal Plain, Illinois Basin, Michigan Basin, Pacific Northwest) across the United States. Hydrogeologic reservoir models are developed for the following eight metropolitan area cities within those regions to evaluate RTES performance across different climates and subsurface conditions: Albuquerque, New Mexico; Charleston, South Carolina; Chicago and Decatur, Illinois; Lansing, Michigan; Memphis, Tennessee; Phoenix, Arizona; and Portland, Oregon. Evaluated metrics include estimated required well spacing, thermal storage capacity, and thermal recovery efficiency through time. Also considered for each reservoir are potential complicating factors, including reservoir depth, thermally driven free convection, and groundwater salinity. This work focuses on direct-use cooling because the need for cooling modern office buildings greatly exceeds that for heating in most parts of the country (Falta and others, 2016); however, the evaluated metrics are also relevant to heating and electricity applications. Results indicate that favorable RTES conditions exist in each region, with the Coastal Plain and Basin and Range being especially favorable for thermal storage capacity, while the Pacific Northwest and Michigan Basin excel at energy recovery for the evaluated cooling application. The results underscore the utility of developing maps of thermal storage capacity, subsurface temperature models, and volumetric estimates of thermal storage capacity to serve as key RTES resource classification standards. Overall, this pre-assessment provides a basic understanding of RTES potential in several cities and geologic regions throughout the country and will aid ongoing thermal energy storage assessment efforts.

Arizona, Illinois, Michigan, New Mexico, Oregon, S↗
USGS2026

Ice sheet dynamics drive pronounced changes in the subsurface freshwater-saltwater interface

Saltwater is migrating into freshwater aquifers globally with water quality and biogeochemical implications, yet saltwater intrusion in glaciated regions is sparsely investigated. Field observations suggest that groundwater head in glaciated systems is influenced by ice sheet forcings and provides evidence that seawater infiltrated into offshore aquifers during past deglaciation events. To understand links between ice sheet dynamics, groundwater head, and saltwater intrusion, we use numerical models to explore the effects of deglaciation on nearshore head and salinity distributions. We find that ice sheet thinning diminishes groundwater head, and the resulting shift in subsurface pressure gradients drives rapid landward movement of the subsurface freshwater-saltwater interface up to 4.0 km or 1.3 m per m ice sheet loss. Results highlight an overlooked saltwater intrusion mechanism that aligns with field observations and affects glaciated coastlines undergoing ice sheet retreat, underscoring the need to consider this mechanism in studies of contemporary coastal water quality.

Geophysical Research Letters↗
USGS2026

The digital archivist: Automating legacy macroseismic data processing using large language models

Macroseismic data are a key resource to investigate shaking and damage from preinstrumental and early instrumental eras. However, data are often stored as inconsistently formatted reports describing observed shaking and damage, making manually parsing and interpreting accounts labor‐intensive. We introduce a novel workflow using Google’s Gemini 2.5 Pro large language model (LLM) to automate the extraction and structuring of macroseismic observations from summary reports. We apply this workflow to the 22 March 1957 M 5.3 Daly City, California, earthquake as a case study. We used Gemini to extract addresses, originally assigned modified Mercalli intensity values, and descriptions from each report. To address coordinate precision limits, addresses were geocoded via Google’s Geocoding application programming interface. This workflow yielded over 2300 geocoded intensity reports for the Daly City earthquake. We use the geocoded accounts, with the original report intensity assignments, to develop a shaking intensity map that in some respects rivals modern Did You Feel It? Maps. We also extract and present data for the 9 February 1971 M L 6.7 Sylmar, California, earthquake. Our results demonstrate the potential of LLMs for reliably extracting and analyzing large, unstructured macroseismic datasets. LLMs offer a scalable solution for rapidly digitizing macroseismic archives, enabling their broader use to constrain ground‐motion models in modern seismic hazard analysis and to improve our understanding of site effects in urban areas. The concepts explored here may also be applied to the handling of other legacy seismological and earth science data.

Seismological Research Letters↗
USGS2026

Corundum discovered by SuperCam and the Perseverance rover at Jezero crater, Mars

Mars is primarily composed of mafic mineral assemblages and their alteration products, but small, scattered rocks strewn across the landscape offer clues to greater petrological diversity. While traversing the Jezero crater rim, the Perseverance rover encountered several plagioclase-rich light-toned float rocks. SuperCam identified the distinctive signature of corundum (α-Al 2 O 3 ) in these rocks using time-resolved luminescence spectroscopy. Two strong peaks (692.7 and 694.1 nm) with millisecond lifetimes, and additional supporting lines, are consistent with Cr 3+ substitution for Al 3+ in corundum. Corundum forms in Al-rich, Si-depleted environments through magmatic or metamorphic processes. Given the rocks' small size, association with plagioclase, and location on the crater rim, we interpret the most plausible formation scenario to be impact induced metamorphism at the interface of a felsic and a mafic/ultramafic member with the likely action of fluids at some stage, although other possibilities are not excluded.

Geophysical Research Letters↗
USGS1999

Reconfigured-Channel Monitoring and Assessment Program

A U.S. Geological Survey program to monitor and assess channel reconfiguration activity is described. A data base available on the world wide web will enable land-management agencies and other interested parties to evaluate the long-term success of various channel reconfiguration projects. A demonstration project on the Lake Fork of the Gunnison River, Colorado, illustrates the program objectives and approach.

Colorado↗
USGS2025

Updating regional‐scale geospatial liquefaction models with locally available geotechnical data

We present a method to update the geospatial liquefaction model used by the U.S. Geological Survey’s near‐real‐time ground failure product with subsurface geotechnical data. The geospatial model estimates liquefaction probability from peak ground velocity (via ShakeMap) and geospatial susceptibility proxies. In many regions, additional information relevant to constraining liquefaction likelihood is also available, including surface geology maps and subsurface geotechnical measurements. There is currently no mechanism to use these data in the ground failure product liquefaction model, even though these data could provide more precise constraints on spatial variations in the lithologic character of the soil (surface geology) and direct measurements of the subsurface mechanical properties that affect liquefaction occurrence and severity (geotechnical measurements). In this study, we develop a method to integrate these data with the geospatial model and assess how these data can improve regional‐scale predictions. We develop a Bayesian updating framework and apply it to the 1989 magnitude 6.9 Loma Prieta, California, earthquake, for which mapped observations are available to evaluate performance. We constrain the Bayesian framework with 373 Northern California cone penetration tests and liquefaction susceptibility classes based on the mapped surface geology. This Bayesian model incorporates geotechnical information into the geospatial model and more accurately predicts liquefaction occurrences than the geospatial model, while sacrificing less accuracy in terms of predicting the absence of liquefaction than the geotechnical model. In future applications, this approach could be adapted to update other geospatial models using locally available subsurface data.

California↗
USGS2026

Key minerals in energy transmission

The electrical grid relies on minerals with unique properties to transform, transmit, and convert electricity. Transmission needs are growing—to support the expansion of data centers and link diverse energy resources into the grid. This expansion will require greater quantities of minerals. Although some minerals are produced in the United States, the United States relies on imports for others. This infographic shows the key minerals that support energy transmission, and how reliant the United States is on imports for these minerals.

General Information Product↗
USGS2026

Geochemical disequilibrium at the brittle-ductile transition

We investigate the microtextural, microchemical, and isotopic effects of late-stage ductile deformation in quartzite mylonites and kyanite–muscovite–quartz veins from the Raft River shear zone (Utah). Quartz microstructures record pervasive disequilibrium, expressed by unannealed features including undulatory extinction, deformation lamellae, and poorly defined fabrics, typical of waning deformation in shear zones. Microchemical disequilibrium is best preserved in kyanite quartzites, where CL-zoned kyanite records repeated fractures, overgrowth, and mineral precipitation, and in muscovite from muscovite-poor quartzite mylonites that shows minor-element zoning consistent with syn-deformational overgrowth on detrital cores. In contrast, muscovite from muscovite-rich kyanite quartzites exhibits minimal chemical zoning. These microchemical variations correlate with 40 Ar/ 39 Ar age systematics. Chemically zoned muscovite preserves variable single-step ages, including ∼150 Ma ages, reflecting retention of detrital cores. In contrast, syndeformational muscovite consistently yields Miocene ages, indicating recrystallization and new growth below argon closure temperatures that reset inherited isotopic signatures. Similar trends are observed in quartzite mylonites, where increasing quartz recrystallization and stronger crystallographic preferred orientations occur toward deeper structural levels. Together, these observations indicate increasing retrograde deformation and recrystallization with depth in the Raft River shear zone and demonstrate that strain-driven recrystallization exerts a first-order control on muscovite ages. We suggest that apparent thermochronologic gradients in retrograde shear zones may reflect recrystallization gradients rather than temperature gradients. Where cooling limits the thermal driving force for recrystallization, isotopic relics are preserved, and local deformation and fluid availability control re-equilibration. Consequently, isotopic disequilibrium—particularly in the 40 Ar/ 39 Ar system—may be the rule rather than the exception of retrograde tectonic environments.

Idaho, Utah↗
USGS2024

Using stable oxygen isotope dual-inlet isotope-ratio mass spectrometry to elucidate uranium transport and mixed 230Th/U calcite formation ages at the seminal Devils Hole, Nevada, natural laboratory

Rationale Vein calcite in Devils Hole has been precipitating continuously in oxygen-isotope equilibrium at a constant temperature for over 500 000 years, providing an unmatched δ 18 O paleoclimate time series. A substantial issue is that coeval calcite (based on matching δ 18 O values) has uranium-series ages differing by 12 000 years. Methods An unparalleled high-accuracy δ 18 O chronology series from continuously submerged calcite was used to correct the published uranium-series ages of non-continuously formed calcite in two cores, cyclically exposed by water-table decline during glacial–interglacial transitions. This method relies on the premise that the δ 18 O values of coevally precipitated calcite are identical, allowing matching calcite δ 18 O values to establish formation ages. Results Exposed calcite can have apparent ages that are 12 000 years too young due to unrecognized uranium mobility and resulting mixed ages identified in over 50 mixed uranium-series ages from previous studies. Secondary uranium in fluids, sourced from the formation or dissolution of porous carbonate deposits (folia) with high uranium-238 ( 238 U) concentrations, has migrated up to 10 mm into vein calcite. Conclusions The continuously submerged Devils Hole δ 18 O chronology is not explained by orbital forcing. Rather, this chronology represents a regional climate record in the southern Great Basin of sea-surface-temperature (SST) variations off California, variations that preceded the last and penultimate deglaciations by 5000 to approximately 10 000 years. Temporal discrepancies between the continuously submerged Devils Hole chronology and other regional δ 18 O records (e.g., the Leviathan chronology) can be explained by unrecognized cryptic, pernicious uranium mobility, leading to model estimations that may be thousands of years younger than actual ages. Consequently, paleo-moisture availability, water-table, and groundwater recharge models based on these mixed uranium-series ages are too young by as much as 12 000 years. The potential for post-formation uranium addition in subaerial cores and speleothems underscores the need for caution in uranium-series dating, highlighting δ 18 O time-series comparisons as a method for identifying mixed ages.

Nevada↗
USGS2025

Megafaunal community structure on ferromanganese and phosphorite hardgrounds in the Southern California Borderland

The Southern California Borderland (SCB) is a topographically complex region on the active continental margin that hosts varied hardground habitats, including ferromanganese (FeMn) crusts and phosphorites, marine minerals being considered for resource extraction. The SCB is influenced by seasonal upwelling and terrestrial inputs, and has a well-defined oxygen-minimum zone (OMZ). We analyzed megafaunal community composition, density, and diversity across 41 ROV video transects at 10 SCB hard substrate sites spanning depth (378–2765 m), temperature (1.79–7.96 °C), and oxygen (3.91–105.67 μM) gradients, with varying rock types (FeMn crusts, phosphorites, other rocks). We counted 32,426 individuals representing 146 unique taxa over a total area of 21,935 m 2 . Echinodermata contributed 48 % of the total abundance; Cnidaria 24 %; Porifera 13 %; Annelida 6 %; Arthropoda 6 %; Chordata 2 %; Mollusca, Foraminifera, and Hemichordata <1 %. Megafauna communities showed clear heterogeneity, with density, diversity, and community composition varying among transects within sites and across sites. Rock type recorded no influence on megafaunal density, but habitats that included FeMn crusts had higher diversity and distinct taxonomic composition. Location, depth, and oxygen were the most important variables explaining variance among megafaunal communities, with distinct communities formed at deeper depths and within the OMZ. This study provides baseline information about the megafauna inhabiting SCB hardgrounds, particularly those associated with minerals considered for their resource potential. It extends existing documentation of such communities to greater depths, reveals increased representation of sponges and corals below 1000 m, and highlights the variability within and among deep-sea megafaunal communities, offering context for informed policy decisions.

Southern California borderland↗
USGS2026

Beaver dam analogs as nature-based solutions to mitigate snowpack loss in northern New Mexico

Reductions in snow have left many streams in northern New Mexico dry or with very low flows during the summer. Base flow, defined as the contribution of groundwater to streamflow, can sustain streamflow during the summer and during drought conditions. Beaver dam analogs can be used to “slow the flow,” or increase the infiltration of rain-based runoff to replenish shallow groundwater reservoirs and increase base flow and potentially increase summertime streamflow.

New Mexico↗
USGS2025

Climate-driven sulfate export in alpine watersheds may stimulate methylmercury production

Climate change is increasing sulfate export and changing wetland extent in mountain regions. These changes may increase microbially mediated production of the neurotoxic substance methylmercury due to enhanced sulfate metabolism in mountain environments. Here, we assess methylmercury concentrations and formation rates across high-elevation wetlands in the Colorado Rocky Mountains. We also investigate sulfate controls on methylmercury production within subalpine peatlands by amending soils with sulfate to mimic increased stream export of sulfate from the alpine zone and measuring methylmercury formation rates for different sulfate treatments. We found that subalpine peatlands have statistically significant higher methylmercury concentrations and formation rates compared to alpine, mineral-soil wetlands. Methylmercury production in subalpine peatlands also increased significantly ( p < 0.05) following sulfate additions; the highest rates occurred in sediments with intermediate extractable sulfate concentrations (∼0.60–1.4 mg sulfate g −1 dry soil). Our study is the first to identify soil sulfate-related thresholds for methylmercury production and sulfate-limitation of methylmercury production in subalpine peatlands. These findings highlight important linkages between climate-driven mineral weathering and mercury cycling in mountain regions globally.

Colorado↗
USGS1977

A multi-channel digital telemetry system for low frequency geophysical data

An inexpensive general purpose digital telemetry system for collection of low frequency geophysical data from U.S. Geological Survey instruments (eg. tilt, strain, gravity, creep, water level, radon, magnetic field, resistivity, telluric current, temperature, etc.) has been designed and built. This system provides data for a more general interactive data acquisition, retrieval and analysis system. The field stations are self-contained, battery operated and housed in weather proof containers. Each accepts up to 15 analog data inputs in the range of -5 to +5 volts. The dynamic range is 70db. The units transmit information as FSK (Frequency Shift Keyed) tones onto either a phone line or radio link with up to 150 transmitters sharing one line. The average power consumption is 0.06 nR watts where n is the 1 number of input channels transmitted and R is the sample rate in minutes -1 . The central receiver-recorder unit accepts and decodes the FSK tones and converts, formats and records the digital data together with time information and station identification on IBM combatible magnetic tape. The digital data are also converted and recorded in analog form for visual monitoring.

Open-File Report↗