USGS Science⌕ Search

SEARCH · USGS Science

Results for “Scientific Investigations Report”

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.

923 records · Page 25Linked to original sources
USGS2025

Cancer risk and estimated lithium exposure in drinking groundwater in the US

Importance Lithium is a naturally occurring element in drinking water and is commonly used as a mood-stabilizing medication. Although clinical studies have reported associations between receiving lithium treatment and reduced cancer risk among patients with bipolar disorder, to our knowledge, the association between environmental lithium exposure and cancer risk has never been studied in the general population. Objectives To evaluate the association between exposure to lithium in drinking groundwater and cancer risk in the general population. Design, Setting, and Participants This cohort study included participants with electronic health record and residential address information but without cancer history at baseline from the All of Us Research Program between May 31, 2017, and June 30, 2022. Participants were followed up until February 15, 2023. Statistical analysis was performed from September 2023 through October 2024. Exposure Lithium concentration in groundwater, based on kriging interpolation of publicly available US Geological Survey data on lithium concentration for 4700 wells across the contiguous US between May 12, 1999, and November 6, 2018. Main Outcome and Measures The main outcome was cancer diagnosis or condition, obtained from electronic health records. Stratified Cox proportional hazards regression models were used to estimate the hazard ratios (HRs) and 95% CIs for risk of cancer overall and individual cancer types for increasing quintiles of the estimated lithium exposure in drinking groundwater, adjusting for socioeconomic, behavioral, and neighborhood-level variables. The analysis was further conducted in the western and eastern halves of the US and restricted to long-term residents living at their current address for at least 3 years. Results A total of 252 178 participants were included (median age, 52 years [IQR, 36-64 years]; 60.1% female). The median follow-up time was 3.6 years (IQR, 3.0-4.3 years), and 7573 incident cancer cases were identified. Higher estimated lithium exposure was consistently associated with reduced cancer risk. Compared with the first (lowest) quintile of lithium exposure, the HR for all cancers was 0.49 (95% CI, 0.31-0.78) for the fourth quintile and 0.29 (95% CI, 0.15-0.55) for the fifth quintile. These associations were found for all cancer types investigated in both females and males, among long-term residents, and in both western and eastern states. For example, for the fifth vs first quintile of lithium exposure for all cancers, the HR was 0.17 (95% CI, 0.07-0.42) in females and 0.13 (95% CI, 0.04-0.38) in males; for long-term residents, the HR was 0.32 (95% CI, 0.15-0.66) in females and 0.24 (95% CI, 0.11-0.52) in males; and the HR was 0.01 (95% CI, 0.00-0.09) in western states and 0.34 (95% CI, 0.21-0.57) in eastern states. Conclusions and Relevance In this cohort study of 252 178 participants, estimated lithium exposure in drinking groundwater was associated with reduced cancer risk. Given the sparse evidence and unknown mechanisms of this association, follow-up investigation is warranted.

contiguous United States↗
USGS2026

Patterns and drivers of cliff erosion in Big Sur, California, USA using repeat photogrammetry, 2017–2023

Seacliff erosion in steep terrain poses major risks to transportation and critical infrastructure. In Big Sur, California, USA, seacliff erosion threatens the sustainability of the central coast stretch of California State Route 1, a transportation corridor that is critical to the region's economy. Published cliff retreat rates for the region range from 1 to 40 cm yr −1 , highlighting that high-resolution, process-based studies could enhance understanding of the causes of spatial and temporal variability. We quantified cliff erosion and investigated its drivers along ∼13 km of the Big Sur coastline at week–month timescales during the late fall to early spring wet seasons between January 2017 and June 2023 by analyzing 3D point clouds developed from aerial imagery using four-dimensional structure-from-motion (4D SfM) photogrammetry techniques. We calculated cliff face retreat rates of 2.23 ± 3.06 cm yr −1 (mean ±1 σ ), an order of magnitude lower than long-term estimated rates for the region (which included large deep-seated landslides), but in line with short-term rates reported across California. Change detection imagery comparison, cliff profiles through time, and statistical analysis reveal a cyclical cliff evolution process in which erosion by wave action at the cliff base destabilizes the cliff and primes it for subsequent failure during precipitation events. Although more erosion by volume could be attributed to precipitation-induced increases in soil moisture (784 m 3 km −1 yr −1 ) compared with erosion attributed to wave power (282 m 3 km −1 yr −1 ), our observations underscore the coupled nature of these processes in driving cliff evolution, consistent with established theory and observations.

California↗
USGS2024

Constraining the earthquake recording threshold of intraslab earthquakes with turbidites in southcentral Alaska’s lakes and fjords

Strong ground motion from intraslab earthquakes, which do not produce primary paleoseismic evidence, may initiate gravity-driven turbidity flows in subaqueous basins. The resulting deposits (turbidites) can provide a paleoseismic proxy if the conditions that initiate these flows are known. To better constrain the initiating conditions, we use two recent intraslab earthquakes in southcentral Alaska, the M w 7.1 30 November 2018 Anchorage earthquake and the M w 7.1 24 January 2016 Iniskin earthquake, as calibration events. Through a multilake investigation, we document the occurrence, or the absence, of earthquake-generated turbidity flows from these two events. Both earthquakes are recorded by centimeter-scale turbidites that can be differentiated from climatically generated deposits, as well as other seismic sources based on deposit thickness, sedimentological properties, and deposit age. We show that a Modified Mercalli Intensity (MMI) of ∼V–V1/2 is the minimum shaking intensity required to generate localized sediment remobilization from deltaic slopes, and an MMI of ∼V1/2 is required to produce a deposit of sufficient thickness that a seismic origin can be confidently assigned. The documentation of seismically generated deposits in quick succession (∼2 years) with diagnostic features highlights the utility of using recent earthquakes as calibration events to investigate the subaqueous response to strong ground motion.

Alaska↗
USGS2025

Lunar grid systems, coordinate systems, and map projections for the Artemis missions and lunar surface navigation

Foreward This document contains design specifications of a navigational standard for the Moon, including a Lunar Transverse Mercator system, a Lunar Polar Stereographic system, a Lunar Grid Reference System, and a unique coordinate structure, Artemis Condensed Coordinates, for Artemis mission navigation and lunar surface science. The National Aeronautics and Space Administration (NASA) Artemis campaign seeks to place humans on the Moon for the first time since the Apollo missions. Early Artemis missions are heavily focused on the lunar south pole, which promises to return valuable data on the Moon’s geologic record, amongst other mission objectives. Coordinate systems in use today for the lunar south pole provides crew members on the surface neither an efficient nor intuitive means to communicate their position and orientation. A novel grid coordinate system, the Lunar Grid Reference System, is proposed to address these concerns for use in real-time extravehicular activity operations on the lunar surface. The many stakeholders involved in the Artemis missions will need a common system to communicate position and orientation while astronauts are operating on the lunar surface. To that end, Artemis crew members will need that system to be efficient and intuitive to promote efficient extravehicular activity timelines and reduce confusion. In the context of this document, these characteristics are addressed on the design of lunar coordinate systems: Efficient.—The number of characters required to communicate a location within a desired precision level in both local and global contexts, and how many steps are required for a recipient or sender to interpret a location. Intuitive.—How well the system aligns with human perceptual abilities, and whether the system yields distances that have the same relationship to actual lunar surface distance in all directions from the point where a person is located. Technological systems are currently being investigated to supplement the crew members’ ability to locate and orient themselves and other assets on the lunar surface; however, it is unlikely that those systems will be fully operational for the first few landed missions. Even with future positional aids, crew members will still need an efficient and intuitive means to communicate position and orientation. In addition, if technological systems fail, the crews will require land navigation skills and have maps available, thus providing further motivation for a crew-centric coordinate system. The contents of this U.S. Geological Survey (USGS) document detail a comprehensive framework for standardizing lunar crewed surface navigation within NASA and outlines the protocols, methods, and designs necessary for achieving consistency and interoperability across relevant space mission teams and lunar surface navigators. Key components of this document include designs of map projections, projected coordinate reference systems (Lunar Transverse Mercator and Lunar Polar Stereographic systems), and a grid system (Lunar Grid Reference System and Artemis Condensed Coordinates) for the Moon. The work proposed in this document seeks to accomplish something similar to the National Geospatial-Intelligence Agency (NGA) document SIG 0012 (NGA, 2014a), but for using grid systems for the Moon. This report incorporates initial feedback and input from NASA’s Artemis Geospatial Data Team, NASA’s Flight Operations Directorate, National Geodetic Survey, USGS Astrogeology Science Center, and NGA and is intended to serve as a resource for all involved with the Artemis missions, as well as for engineers designing and operating lunar infrastructure.

Techniques and Methods↗
USGS1973

Geohydrology of the Cross-Florida Barge Canal area, with special reference to the Ocala vicinity

The Cross-Florida Barge Canal route commences at Palatka on the St. Johns River, about 75 miles upstream from the Atlantic Ocean, and extends 110 miles southwestward across Peninsular Florida into deep water in the Gulf of Mexico near Yankeetown. The canal will be equipped with five locks, each 600 feet long and 84 feet wide, and the channel will be a minimum of 12 feet deep and 150 feet wide. From near Ocala northeastward, the canal channel will replace much of the natural channel of the Oklawaha River, and will be excavated into beds of the so-called shallow sand aquifer of Miocene age and younger, which overlies limestone of the Floridan aquifer. Westward from Ocala, most of the canal will be excavated below the potentiometric surface into limestone and dolomite of the Floridan aquifer. Water levels of Rodman, Eureka, and Inglis Pools will be controlled by dams and spillways with the limited exchange of water between the pools and the aquifers. The water levels in the Summit Pools will fluctuate with the natural changes in the ground-water level of the Floridan aquifer, although the stage of the pool will be controlled partly by the stage held in the Eureka Pool. A dynamic inflow-outflow relationship will exist between the Summit Pool and the Floridan aquifer. The Floridan aquifer in the canal area is 1,000 to 1,200 feet thick and consists of limestone and dolomite of middle Eocene Miocene age, including from older to younger, the Lake City, Avon Park, and Ocala limestones plus permeable sandy, dolomitic limestone in the lower part of the Hawthorn Formation. It is possible that most of the flow to the two major springs in the area occurs in the upper 100 feet or so of the aquifer in the Ocala Limestone. The aquifer is underlain by the Oldsmar limestone of early Eocene age and is overlain by sand, clayey sand, clay and shell beds of Miocene through Holocene age, in thickness from a few feet to 300 feet. The permeable beds overlying the Floridan aquifer constitute the shallow aquifer, while the poorly permeable ones act as confining beds where the Floridan aquifer is under artesian conditions. A north-south line drawn separating the head of Silver Springs on the west from the Oklawaha River on the east marks the approximate western limit of a continuous blanket of materials of Miocene-Pliocene(?) age covering the rocks of the Floridan aquifer. East of the line, much of the aquifer is under artesian conditions, particularly in the Oklawaha River valley, although in some areas east of the valley, direct recharge through thick permeable Miocene-Pliocene(?) sands occurs. West of the line, only scattered remnants of a once continuous Miocene-Pliocene(?) cover remain. Lack of the cover is a result of erosion on the crest and flank of the Ocala Uplift, a broad northwest-southeast trending anticlinal upwarp, the axis of which is crossed by the canal route in the Dunnellon area. Over most of this area the Floridan aquifer is unconfined and receives direct recharge through a cover of a few tens of feet of sand and clayey sand of Quaternary age. Tensional stresses during the structural evolution of the Ocala Uplift produced an intersecting system of fractures and normal faults in rocks of the Floridan aquifer. The fractures and faults are important controls for orientation of solution channels and, therefore, for development of ground-water circulation patterns. When the system surface streams, which once drained the Barge Canal area, eroded the poorly permeable Miocene-Pliocene(?) cover from the flanks of the Ocala Uplift, surface runoff was reduced and precipitation began to directly infiltrate the underlying limestones. Now only principal streams remain, such as the Oklawaha and Withlacoochee Rivers and a few short tributaries, while one of the most highly developed subsurface drainage systems in the world has evolved in cavernous limestone of the Floridan aquifer. Two of the larger freshwater springs in the world now discharge from the Floridan aquifer in the canal area. Silver Springs near Ocala discharges an average 531 mgd (million gallons per day) down the 4-mile long Silver River, which flows on poorly permeable beds to the Oklawaha River. Rainbow Springs near Dunnellon discharges on average 468 mgd from numerous orifices in the bed of the 5-mile-long Rainbow River, which flows into the Withlacoochee River. The heads of the springs have migrated to their present positions partly because of a tendency of ground-water levels to decline as permeability in the aquifer is increased due to removal of limestone by solution, and because of mechanical erosion of the limestone in the vicinity of the spring heads. Also, points of principal spring discharges have shifted in the past due to changes in ground-water levels in response to changes in sea level. The subsurface drainage system is continuing to evolve today, as evidenced in part by frequent occurrence of new sinkholes and by the presence of significant amounts of calcium bicarbonate in the spring waters. Rodman Pool, at the east end of the canal, is separated from the Floridan aquifer by poorly permeable materials. The pool's operating water level will be only a few feet above the potentiometric surface at the downstream end, and at or slightly below the potentiometric surface at the upstream end. Little exchange of water between the Rodman Pool and the Floridan aquifer is expected. Eureka Pool, just upstream from Rodman Pool, will also be separated from the Floridan aquifer by poorly permeable beds. However, the stage of the pool will be about 15 feet higher than the natural potentiometric surface at the pool's downstream end, and some seepage into the Floridan aquifer is anticipated through faults and leaky parts of the poorly permeable beds, with a consequent rise in ground-water levels in areas adjacent to the lower end of the pool. Possibilities for particulate contamination of the aquifer will tend to be minimized because of the filtering capacity of the materials through which water must pass to reach the aquifer, although the natural filter will not preclude movement into the aquifer of contaminants which might become dissolved in the pool waters. No significant interchange of water between the pool and the aquifer is expected at the upstream end of the Eureka Pool. Present construction plans indicate an operating stage for Eureka Pool which will range between 38 and 40 feet above mean sea level, although it is possible to dredge the pool deep enough to permit a range in stage of 36 to 40 feet. A backwater effect extending up Silver River from the Eureka Pool is expected to regulate the stage at the head of Silver Springs between 39 and 44 feet above mean sea level if the pool ranges between 36 and 40 feet. If Eureka Pool ranges only between 38 and 40 feet, the range of stage at the head of the springs should be about 41 to 44 feet above mean sea level. From Inglis Lock west, the canal will have direct connection with the Gulf, and canal stage will fluctuate with the Gulf tide. Since the canal stage will be slightly lower than the adjacent ground-water levels along much of the reach, there will be some ground-water inflow to the canal. No significant changes in the existing ground-water regime are expected in the vicinity of Inglis Pool, the first step up in the canal east of the Gulf. Existing ground-water and surface-water levels in the area will not change appreciable, and the natural stage and flow of Rainbow Springs, which will flow by way of Rainbow River into Inglis Pool, should not be affected by canal operations. A possible adverse effect of the canal on the Inglis Pool area could result if sea water is locked up from the Gulf through Inglis Lock. However, the high step of 25 feet at the lock, flushing action of continuous flow from Inglis Pool to the lower reaches of the Withlacoochee River, and use of possible preventive locking procedures should minimize the problem. The potential for adverse effects on the ground-water regime is greatest in the area of the Summit Pool. Through most of the length of the pool, the canal channel will be excavated into limestone of the Floridan aquifer to depths of 12 to 27 feet below the potentiometric surface. Changes that will take place in the ground-water flow system in the Silver Springs drainage area, once the canal is completed, were estimated by flow-net analysis. Variation in aquifer transmissivity was determined by calculating transmissivity in 25 different flow cells surrounding Silver Springs. Transmissivity in the 25 cells averages about 15,600,000 gpd/ft [2,090,000 ft2/day (feet squared)], but transmissivity in the six cells through which the Summit Pool passes ranges from 9,000,000 to 44,000,000 gpd/ft (1,210,000 to 5,900,000 ft2/day). Transmissivity was used to compute static stage of the Summit Pool under given ground-water level conditions. Had the canal existed in May 1968 and had the stage of Eureka Pool been held at 36 feet at the time, the static stage in Summit Pool would have been about 42.1 feet above mean sea level. Thus, a conceptual model of the changes in the potentiometric surface wrought by the finished canal was drawn, and zones of ground-water inflow and outflow were delineated. Most outflow from the Summit Pool to the aquifer should be limited to one 4-mile-long zone along the north side of the pool, about 5 miles south of Silver Springs. It is estimated that a water volume equivalent to about 8 percent of the daily flow of Silver Springs will enter the Summit Pool each day from the southern one-third of the Silver Springs drainage area. A like amount will reenter the aquifer at the main zone of outflow and move toward Silver Springs at an estimated average velocity of about 200 feet per day, if something close to the natural static stage of the pool is maintained by return pumpage of the lockage losses. At a velocity of 200 feet per day, water from the Summit Pool would discharge at Silver Springs about 140 days later. However, any estimate of velocity in the highly cavernous limestone aquifer in the area should be used with caution, because difficult to measure changes in porosity and thickness of major zones of flow may cause large variations in velocity. If all lockage losses are returned to the Summit Pool by pumping from Eureka Pool, no net loss from the Silver Springs drainage area, except for some evaporation from the water surface in the canal and possible leakage around locks, will result from canal operations. The zone of outflow from the Summit Pool to the aquifer will be in a natural potentiometric trough, and the zone of inflow will be in a potentiometric ridge area. The equilibrium water level in the pool will tend to be about 1 foot higher than the altitude of the lowest level in the pre-canal potentiometric trough, and about 2 feet lower than the highest level on the pre-canal potentiometric ridge. West of the Silver Springs drainage area just east of Dunnellon Lock, in the area of a local potentiometric high, the water level in the Summit Pool is expected to be about 15 feet below the natural potentiometric surface. In most areas 2 to 3 miles away from the Summit Pool, effects of the canal on the natural potentiometric surface should be slight. The stage of the Summit Pool, judging from the 36-year record for ground-water level changes and the anticipated indirect effect of the controlled stage in Eureka Pool, should have a maximum of about 10.5 feet with a maximum water level of about 51.5 feet above mean sea level and a minimum of about 41.0 feet above mean sea level. Of particular importance in the Summit Pool is an implementation of well planned construction and operational procedures designed to minimize risks of ground-water contamination.

Florida↗
USGS2026

A temporal look at the influence of topographic amplification on earthquake-triggered landslides in 3D seismic simulations

Earthquakes are a primary trigger for landslides, often leading to catastrophic consequences. While numerous studies have explored the spatial distribution of earthquake-triggered landslides, understanding the interaction between seismic waves and topography remains a critical challenge. Topographic irregularities can cause seismic wave amplification, altering ground shaking, and can trigger landslides that are challenging for predictive models to anticipate. This study investigates the spatial and temporal evolution of topographically amplified landslides, focusing on coseismic landslides triggered by the Mw 7.5 mainshock of the 2018 Papua New Guinea earthquake and post-seismic landslides associated with its four aftershocks, each exceeding Mw 6.0. We employ low-frequency, three-dimensional numerical ground shaking simulations and data-driven multivariate analyses to examine how landslides evolved from the coseismic to post-seismic periods. Our findings reveal a spatial shift in landslide distribution, in which the mainshock triggered slope failures predominantly on steep hillslopes, and the aftershocks triggered landslides on gentler slopes, often near geologic boundaries. We attribute this transition partly to the earthquake legacy effect of the mainshock, where the mainshock caused weakening of these hillslopes, making them more prone to failure when aftershocks occur. Additionally, the concentration of failures along geologic contacts in the post-seismic phase suggests that site amplification, stemming from contrasts in subsurface materials, exerts a key influence on landslide occurrence. Although not explicitly captured in our current numerical simulations, this mechanism warrants further investigation for more accurate hazard modeling.

JGR Solid Earth↗
USGS2025

Assessing earthquake risks to lifeline infrastructure systems in the United States

The security and economic stability of the United States rely heavily on robust lifeline infrastructure systems and yet the risks to such systems are seldom quantified at the national scale. For example, while earthquake risks to buildings in the United States have been investigated at the national scale regularly, such risks to gas pipelines have rarely been investigated nationally. In this paper, we use examples from two critical infrastructure sectors to demonstrate (1) the nature of earthquake risks to lifeline infrastructure systems, (2) complexities involved in regional seismic risk assessments, and (3) how such risks change with time. We found that bridge risks can be underestimated by at least 64 % when viewed from repair costs instead of traffic demands and that regional risks can be underestimated by 19 % when spatial correlations of ground motion are ignored. Further, exceedance of traffic demand can be 50 times more likely to occur when viewed at the regional scale than when viewed at an individual bridge. Similarly, exceedance of repairs can be 180 times more likely to occur when viewed at the pipeline network level than at a segment-specific level. Finally, sensitivity analyses with the 2018 and 2023 USGS National Seismic Hazard Models indicate an increase in bridge risk of at least 24 % and an increase in exposed gas pipeline mileage of 43 %. The evolution of risks, complexities involved in assessments, and limited resources jointly underscore the need for more routine updates to nationwide seismic risk assessments of lifeline systems in the United States.

International Journal of Critical Infrastructure P↗
USGS2026

Regional variations in sea ice and primary productivity in the Bering Sea during Marine Isotope Stage 11

Marine Isotope Stage (MIS) 11 (424-374 ka) has long been an analog for Holocene climate, because it is the most recent interglacial period with similar orbital conditions. However, there is significant global and regional variability in the climate response to MIS 11 warmth. Here, we review sediment core records from across the Bering Sea to investigate changes in paleoceanographic conditions during Marine Isotope Stages 12-10. Sea ice was present over much of the Bering Sea during MIS 11, but today, none of the sites investigated are ever ice-covered. This suggests that sea ice regimes in the Bering Sea during MIS 11 were different to those of the Holocene. There are also regional differences in the response of sea ice to MIS 11 warming. At the Umnak Plateau, Southeastern Bering Sea, sea ice concentrations decline during deglaciation, but they remain high at the slope sites until Peak MIS 11. Sea ice re-advances over the Umnak Plateau during peak interglacial warmth, at the same time that it declines over the slope sites. Late MIS 11 is characterized by high concentrations of seasonal sea ice at the Umnak Plateau, whilst sea ice at the slope sites fluctuates between consolidated and unconsolidated ice cover. This east-west dichotomy may be explained by changes in the behavior of the Aleutian Low. Productivity increases dramatically during deglaciation due to increased upwelling and sea level rise bringing fresh nutrients into the oceans. This is characterized by increased diatom productivity, increased terrestrial carbon deposition, and laminations at all sites.

Alaska↗
USGS2025

Ecohydrological response of a forested headwater catchment to a flash drought in the Southeastern U.S.

Flash droughts differ from traditionally defined droughts in their rapidity of intensification and often associated high vapor-pressure deficit. These droughts can lead to declines in streamflow and water table depth and induce water stress to vegetation at a greater rate than droughts that manifest over longer periods. However, little is known regarding the response of forested environments to flash drought because most studies of impacts have been conducted in agricultural settings. In this study we investigated water-use patterns of riparian trees using sap flow methods and examined the role of groundwater as a source of moisture over three periods that were delimited by antecedent soil moisture conditions. For a longer-term perspective we also examine monthly streamflow over the 35-year record. We observed that trees at only one monitoring plot showed a decrease in water use relative to evaporative demand during a flash drought. Total reverse sap flow (flow toward the roots rather than the canopy) greatly increased during the flash drought period, suggesting the likely occurrence of hydraulic redistribution to the excessively dry soils. Over the drought period groundwater became a more dominant source of moisture for sustaining forest water use. Monthly mean streamflow during the flash drought approached levels observed in past multiyear droughts. This is the first study, to our knowledge, to specifically investigate the response of multiple water budget components to flash drought in a humid forest. As more studies are conducted, a better understanding of the range of expected responses are likely to emerge.

Georgia↗
USGS2025

Relationship of atmospheric nitrogen deposition to soil nitrogen cycling along an elevation gradient in the Colorado Front Range

Microbial processing of atmospheric nitrogen (N) deposition regulates the retention and mobilization of N in soils, with important implications for water quality. Understanding the links between N deposition, microbial communities, N transformations, and water quality is critical as N deposition shifts toward reduced N and remains persistently high in many regions. Here, we investigated these connections along an elevation transect in the Colorado Front Range. Although rates of N deposition and pools of extractable N increased down the elevation transect, soil microbial communities and N transformation rates did not follow clear elevational patterns. The subalpine microbial community was distinct, corresponding to a high C:N ratio and low pH, while the microbial communities at the lower elevation sites were all very similar. Net nitrification, mineralization, and nitrification potential rates were highest at the Plains (1,700 m) and Montane (2,527 m) sites, suggesting that these ecosystems mobilize N. In contrast, the net immobilization of N observed at the Foothills (1,978 m) and Subalpine (3,015 m) sites suggests that these ecosystems retain N deposition. The contrast in N transformation rates between the plains and foothills, both of which receive elevated N deposition, may be due to spatial heterogeneity not captured in this study and warrants further investigation. Stream N concentrations from the subalpine to the foothills were consistently low, indicating that these soils are currently able to process and retain N deposition, but this may be disrupted if drought, wildfire, or land-use change alter the ability of the soils to retain N.

Colorado↗
USGS2025

Physicochemical properties and bioreactivity of sub-10 µm geogenic particles: Comparison of volcanic ash and desert dust

Exposure to ambient particulate matter (PM) with an aerodynamic diameter of <10 μm (PM 10 ) is a well-established health hazard. There is increasing evidence that geogenic (Earth-derived) particles can induce adverse biological effects upon inhalation, though there is high variability in particle bioreactivity that is associated with particle source and physicochemical properties. In this study, we investigated physicochemical properties and biological reactivity of volcanic ash from the April 2021 eruption of La Soufrière volcano, St. Vincent, and two desert dust samples: a standardized test dust from Arizona and an aeolian Gobi Desert dust sampled in China. We determined particle size, morphology, mineralogy, surface texture and chemistry in sub-10 μm material to investigate associations between particle physicochemical properties and observed bioreactivity. We assessed cellular responses (cytotoxic and pro-inflammatory effects) to acute particle exposures (24 hr) in monocultures at the air-liquid interface using two types of cells of the human airways: BEAS-2B bronchial epithelial cells and A549 alveolar type II epithelial cells. In acellular assays, we also assessed particle oxidative potential and the presence of microorganisms. The results showed that volcanic ash and desert dust exhibit intrinsically different particle morphology, surface textures and chemistry, and variable mineralogical content. We found that Gobi Desert dust is more bioreactive than freshly erupted volcanic ash and Arizona test dust, which is possibly linked to the presence of microorganisms (bacteria) and/or nanoscale elongated silicate minerals (potentially clay such as illite or vermiculite) on particle surfaces.

GeoHealth↗
USGS2024

Effects of stochastically-simulated near-fault ground motions on soil liquefaction

The scarcity of historically recorded near-fault ground motions poses a challenge to systematically understanding the influence of near-fault effects on various types of seismic demands for engineering purposes. In particular, the current state of knowledge of the influence of ground-shaking intensity on soil liquefaction and its consequences does not specifically account for the effects of near-fault ground motion characteristics. In this study, the influence of near-fault ground motions on liquefaction triggering and lateral spreading are investigated using non-linear modeling of a hypothetical liquefiable soil column in the finite-element computational platform OpenSees subjected to simulated ground motion time series that represent strong earthquake shaking in the near field. The simulated ground motion time series and resulting datasets are based on a parametric stochastic model and are developed for a range of source and path parameters to represent a realistic variability of ground motion characteristics. Dependencies between ground motion intensity measures (IMs) and liquefaction demand parameters are investigated for near-fault pulse and nonpulse-like ground motion sets. Evolutionary IMs, such as cumulative absolute velocity (CAV) and the time-varying magnitude-adjusted peak ground acceleration (PGAM), are considered in developing liquefaction triggering probability density functions. Post-liquefaction triggering responses such as lateral spreading displacements are examined in relation to PGA M and CAV. The ground motion simulations are validated by comparing their liquefaction-capacity PGA M fragilities and post-triggering CAV vulnerability relationships to historical records from the 1994 Northridge earthquake in California, USA. Finally, a path forward for future studies that includes finding systematic differences in the IM-liquefaction demand relationships between near-fault and far-field stochastic ground motion sets is outlined.

Conference Paper↗
USGS2025

Relationship of basin structure and bedrock lithology to faulting in the 2019 Ridgecrest earthquake region, California, from gravity and aeromagnetic data

We investigate patterns of cumulative offsets on the faults that ruptured in 2019 and along the Garlock Fault in the Ridgecrest region, California using recently published gravity and aeromagnetic data. We also examine the relationship of basin structure and bedrock structure to the 2019 M7.1 Ridgecrest earthquake ruptures (Fig. 1A), which were primarily along a dextral northwest-striking fault system, and along a sinistral northeast-striking fault, which ruptured hours earlier with a M6.4 event.

California↗
USGS2026

Evaluating machine learning approaches to identify and predict oil and gas produced water lithium concentrations

Recently, the demand for battery-grade lithium has substantially increased, largely due to electrification of the transportation sector. The search for new lithium sources has turned to produced waters (frequently brines), a large-volume wastewater by-product of oil and gas extraction. Geochemical analysis indicates the presence of varying concentrations of lithium from produced water samples collected across the United States and represented in the U.S. Geological Survey’s National Produced Water Geochemical Database, as well as mixtures of Marcellus Shale produced water included in the Pennsylvania Department of Environmental Protection’s Oil and Gas Well Waste Reports. We first examined whether the geochemical signature of the lithium-bearing produced waters is sufficiently distinct so that machine learning (ML) can be used to correctly classify samples to the formation of origin. The produced water sample data used to assess classification accuracy were from the Marcellus Shale, Utica Shale and Point Pleasant Formation (Utica), and Smackover Formation oil and gas wells. Further, we evaluated the potential for ML to accurately classify Marcellus Shale produced water spatially (i.e., northeast versus southwest Pennsylvania). We then investigated whether ML algorithms applied to a suite of geochemical concentration data (i.e. Ba, Br, Cl, K, Mg, Sr) may be used to predict the lithium concentration of an unknown sample. Finally, we applied an estimated economic lithium grade cutoff of 150 milligrams per liter (mg/l) and assessed the utility of ML to predict whether a produced water sample would fall above or below the grade cutoff based on the suite of geochemical parameters. Four machine learning algorithms—Random Forest (RF), Gradient Boosting Trees (GBT), Extreme Boosting (XGBoost), and Deep Neural Networks (DNN) were assessed. This study successfully demonstrates that all four machine learning methods can precisely and accurately estimate lithium concentrations and geologic formation classification. The products of this study contribute to the growing body of knowledge aimed at expanding the lithium resource base within the United States.

Alabama, Arkansas, Florida, Georgia, Louisiana, Mi↗
USGS2026

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

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

Washington↗
USGS2026

The mineral chemistry networks of tin and tungsten reflect metallogenic eras of the Mesozoic

Continental remobilization is a crucial driver of metallogenesis and the formation of ore deposits. Some of the world’s largest mineral deposits of the economically valuable elements tin (Sn), tungsten (W), gold (Au), copper (Cu), lead (Pb), and zinc (Zn) formed during the Mesozoic Era. Additionally, the chemistry and distribution of the elements Sn and W have been investigated in previous studies to understand planetary formation and differentiation processes. These two elements are largely co-located during certain South China Mesozoic metallogenic events but are not co-located during other time periods in the same regions. Here, we investigated the mineral chemistry network similarities and dissimilarities of Sn and W to understand their mineral formation and distribution during the Mesozoic Era and throughout Earth history. Mineral chemistry network community detection analysis and electronegativity associations among mineral constituent elements of Sn minerals and W minerals indicate that the elements have similar chemistry among their oxide minerals. However, Sn forms a much wider range of minerals that also contain S compared to W, which occurs in a limited number of S-containing minerals. The divergent constituent element interactions among S-containing Sn minerals and W minerals reflect the redox sensitivity and importance of oxygen (O) fugacity in Sn mineral formation. Conversely, extensive W mineral deposits are known to form at both high and low O fugacities. The similarities and differences between the mineral chemistry networks of Sn and W reflect the mineral distribution of the two elements in the Sn-W mineralization event from 160 to 139 Ma vs. the Sn–uranium (U) mineralization event from 125 to 98 million years ago (Ma). The mineral chemistry and distribution of Mesozoic Sn and W deposits illustrate the contrasting importance of redox and O fugacity on the mineral formation of different elements, and the dynamic crustal evolution that took place during this period of Earth history.

Geosciences↗
USGS2025

Pre-eruptive characteristics of “suspect” silicic magmas in Carlin-type Au-forming systems

World-class Carlin-type Au deposits hosted in sedimentary rock were formed when profuse Eocene silicic magmatism swept across northern Nevada in response to arc migration. Carlin-type Au deposits formed along with porphyry/skarn Cu-Mo-W-Au deposits, epithermal Ag-Au deposits, and distal disseminated Ag-Au deposits. But unlike these other Au-bearing deposits that have clear associations with igneous intrusions, Carlin-type ore deposits appear to have formed distant from concealed plutons, and their origin remains controversial. Despite decades of abundant geophysical, geochronological, and geochemical studies suggesting the involvement of magmas, concrete evidence for magmatic involvement is still lacking. Consequently, the involvement of contemporaneous igneous systems remains inferred based on age, proximity, and variable isotopic, geochemical, and geophysical clues. A recent synthesis of deposit models postulates that Carlin-type Au deposits are intrusion-related, but that the causative magmas reside deeper (∼6–12 km) than in typical porphyry and peripheral systems (∼3–5 km), meaning that Carlin-type deposits are perhaps more distal expressions of igneous intrusions. We investigate a collection of “suspect” magmatic systems over a ∼7 m.y. timespan (∼41–34 Ma) that are contemporaneous with and near known Carlin-type ore deposits. We report results of a multifaceted array of in situ geochemical analyses (FTIR, EMP, SHRIMP-RG, LA-ICP-MS) of quartz-hosted melt inclusions, biotite, and quartz to better characterize the pre-eruptive characteristics of these magmas. We also report results of thermobarometry and thermodynamic phase equilibria modeling to help place constraints on magmatic reservoir depths and processes. Rather than a single “flavor” of silicic magma, we observe a surprisingly broad compositional spectrum of rhyolites, with one end of the spectrum exhibiting more arc-like (I-type) characteristics and the other end displaying more post-subduction, thick-crust extensional (A-type) characteristics. This broad compositional spectrum suggests a more complex picture of silicic crustal magmatism operating over a narrow span of time during slab rollback. Despite this spectrum, magmatic systems in this study are consistently ferroan and generally peraluminous, which we interpret as an expression of the relatively elevated geotherm at the time and incorporation of variable amounts of highly peraluminous metasedimentary crustal components. The silicic magma spectrum encompasses a range of mineralization associations, including subduction-related Cu-Mo-W-Au-Ag and post-subduction, thick-crust extensional rare-metal Mo-Sn-W-F-Be-Ag-Au, consistent with the prolific and diverse array of ore deposits that formed during this time. Carlin-type Au deposition appears to be associated with nearly the entire magmatic spectrum. This apparent indifference to silicic magma “flavor” would seem to imply that if magmas are involved in Carlin-type Au deposit genesis, they perhaps do not need to be compositionally specialized and/or possibly are only relevant as heat sources driving circulation to remobilize and redistribute metals.

Nevada↗
USGS2024

Distinguishing natural sources from anthropogenic events in seismic data

As seismic data are increasingly used to investigate a diverse range of subsurface phenomena beyond regular fast-rupturing earthquakes (Peng and Gomberg, 2010; Beroza and Ide, 2011), it is important to acknowledge that human-generated ground vibrations may be mistaken for naturally generated subsurface processes (Larose et al., 2015; Li et al., 2018). Correct discrimination of natural processes from anthropogenic noise is especially pressing given the trend in seismic detection research toward automated algorithms and machine learning methods (Yoon et al., 2015; Kong et al., 2019;Mousavi and Beroza, 2022) and the growth in seismic data collection in new environments such as urban and industry settings (e.g., Díaz et al.,2017).

Seismological Research Letters↗