USGS Science⌕ Search

SEARCH · USGS Science

Results for “Behavioral Ecology and Sociobiology”

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,656 records · Page 71Linked to original sources
USGS2024

Environmental monitoring of groundwater, surface water, and soil at the Ammonium Perchlorate Rocket Motor Destruction Facility at the Letterkenny Army Depot, Chambersburg, Pennsylvania, 2021

Letterkenny Army Depot in Chambersburg, Pennsylvania, built an Ammonium Perchlorate Rocket Motor Destruction (ARMD) Facility in 2016 to centralize rocket motor destruction and contain all waste during the destruction process. The U.S. Geological Survey has collected environmental samples from groundwater, surface water, and soils at ARMD since 2016. During 2021, samples were collected from four groundwater wells in September, one surface-water site in October, and five soil sites in November near the facility. Samples were analyzed for nutrients, trace metals, major ions, total volatile organic compounds, and perchlorate. Perchlorate was not detected in any 2021 samples. Groundwater results showed no constituents exceeded any U.S. Environmental Protection Agency (EPA) maximum contaminant level (MCL). Dissolved arsenic (As) was detected in one well above the reporting detection level (RDL) of 3 micrograms per liter (μg/L) at 5.4 μg/L but below its MCL of 10 μg/L. Dissolved iron (Fe) was the only inorganic constituent measured above an EPA secondary maximum contaminant level (SMCL). All groundwater samples collected in 2021 exceeded the Fe SMCL of 300 μg/L, with concentrations ranging from 390 μg/L to 3,500 μg/L. Surface-water data collected during 2021 showed no measured constituents in the surface-water sample that exceeded any EPA MCL or SMCL. Soil samples collected from 2016 through 2021 showed all concentrations of As exceeded the EPA soil screening levels of 3 milligrams per kilogram (mg/kg) but did not exceed the Pennsylvania medium-specific concentrations for As of 61 mg/kg. Arsenic concentrations in 2021 ranged from 9.1 mg/kg to 12.9 mg/kg. The 2021 results for the ARMD Facility indicate no increases in concentrations of reported compounds compared to data from 2016 to 2020. The contained burn treatment facility for demilitarization of rocket motors during 2021 appears to have operated without elevating concentrations of target compounds compared to previous years.

Pennsylvania↗
USGS2024

U.S.-Mexico Borderland & vegetation community map

People on both sides of the United States-Mexico border need a high-resolution, binational vegetation community map that spans the entire United States-Mexico borderlands. Traditionally, mapping efforts in this region were impeded by complex logistics related to the international border, differing national needs and plans, and resource allocations and priorities. To address this need, scientists from the U.S. Geological Survey (USGS) Southwest Biological Science Center partnered with the Sonoran Joint Venture, the U.S. Fish and Wildlife Service (FWS) Migratory Bird Program, data engineers from the Department of Biosystems Engineering at the University of Arizona, and collaborators from the Wildlands Network, the Borderlands Program to produce the first prototype land cover map within the overlapping Mojave Desert, Sonoran Desert, and the North American Bird Conservation Initiative’s Bird Conservation Region 33 (BCR33) using Landsat satellite data . BCR33 is an area of high biodiversity, providing habitat for bird species of concern and other wildlife. The land cover map supports FWS recovery plan efforts related to conservation planning activities for many species, including Yellow-billed Cuckoo ( Coccyzus americanus ), Cactus Ferruginous Pygmy-Owl ( Glaucidium brasilianum cactorum ), Southwestern Willow Flycatcher ( Empidonax traillii extimus ), Yuma Ridgway’s Rail ( Rallus obsoletus yumanensis ), Bendire’s thrasher ( Toxostoma bendirei ), LeConte’s thrasher ( Toxostoma lecontei ), Masked Bobwhite ( Colinus virginianus ridgwayi ), jaguar ( Panthera onca ), and endangered plants such as Bartram’s stonecrop ( Graptopetalum bartramii ) and the Pima pineapple cactus ( Coryphantha robustispina ssp. robustispina ). In 2024, a Phase-II map for the full BCR33 region was completed, increasing the understanding of the binational nature of natural communities. The published map and associated paper can be found here .

Borderland↗
USGS2023

Slow slip detectability in seafloor pressure records offshore Alaska

In subduction zones worldwide, seafloor pressure data are used to observe tectonic deformation, particularly from megathrust earthquakes and slow slip events (SSEs). However, such measurements are also sensitive to oceanographic circulation-generated pressures over a range of frequencies that conflate with tectonic signals of interest. Using seafloor pressure and temperature data from the Alaska Amphibious Community Seismic Experiment, and sea surface height data from satellite altimetry, we evaluate the efficacy of various seasonal and oceanographic pressure signal proxy corrections and conduct synthetic tests to determine their impact on the timing and amplitude prediction of ramp-like signals typical of SSEs. We find that subtracting out the first mode of the complex empirical orthogonal functions of the pressure records on either the shelf or slope yields signal root-mean-square error (RMS) reductions up to 73% or 80%, respectively. Additional correction with proxies that exploit the depth-dependent spatial coherence of pressure records provides cumulative variance reductions up to 83% and 93%, respectively. Our detectability tests show that the timing and amplitude of synthetic SSE-like ramps can be well constrained for ramp amplitudes ≥4 cm on the shelf and ≥2 cm on the slope, using a fully automated detector. The principal limits on detectability are residual abrupt changes in pressure that occur as part of the transition to and from summer to winter conditions but are not adequately characterized by our seasonal corrections, as well as the inability to properly account for instrumental drift, which is not readily separated from the seasonal signal.

Alaska↗
USGS2024

Limited evidence of late Quaternary tectonic surface deformation in the eastern Tennessee seismic zone, USA

The ~300-km-long eastern Tennessee seismic zone (ETSZ), USA, is the second-most seismically active region east of the Rocky Mountains. Seismicity generally occurs below the Paleozoic fold-and-thrust belt within the Mesoproterozoic basement, at depths of 5–26 km, and earthquake magnitudes during the instrumental record have been moment magnitude ( M w )≤4.8. Evidence of surface deformation may not exist or be difficult to detect because of the vegetated and soil-mantled landscape, landslides, locally steep topography, anthropogenic landscape modification, or long, irregular recurrence intervals between surface-rupturing earthquakes. Despite the deep seismicity, analog models indicate that accumulation of strike-slip or oblique-slip displacement at depth could be expected to propagate upward through the Paleozoic section, producing a detectable surficial signal of distributed faulting. To identify potential surface deformation, we interrogated the landscape at different spatial scales. We evaluated morphotectonic and channel metrics, such as channel sinuosity and catchment-scale hypsometry. Additionally, we mapped possible fault-related topographic features on 1-m lidar. Finally, we integrated our observations with available bedrock and Quaternary surficial mapping and subsurface geophysical data. At a regional scale, most morphotectonic and channel metrics have a strong lithologic control. Within smaller regions of similar lithology, we observe changes in landscape metrics like channel sinuosity and catchment-scale hypsometry that spatially correlate with new lineaments identified in this study and previously mapped east–west Cenozoic faults. These faults have apparent left-lateral offsets, are optimally oriented to slip in the current stress field, and match kinematics from recent focal mechanisms, but do not clearly preserve evidence of late Pleistocene or Holocene tectonic surface deformation. Most newly mapped lineaments might be explained by either tectonic or non-tectonic origins, such as fluvial or karst processes. We also re-evaluated a previously described paleoseismic site and interpret that the exposure does not record evidence of late Pleistocene faulting but instead is explained by fluvial stratigraphy.

Tennessee↗
USGS2025

Broadband stochastic simulation of earthquake ground motions with multiple strong phases with an application to the 2023 Kahramanmaraş, Turkey (Türkiye), earthquake

Stochastic ground motion simulation models are often less accurate at lower frequencies than at higher frequencies when fitting recorded data unless supplemented by a deterministic forward directivity velocity pulse model. Moreover, time-modulated stochastic models, which adjust ground motion amplitudes over time, typically use functions that fail to capture multiple strong-motion phases. The February 2023 Turkey (Türkiye) earthquake exhibited diverse recordings, including near-fault and far-field motions with pulse-like and non-pulse-like characteristics, along with single and multiple strong-motion phases. To better represent such a diverse set of recordings, this study enhances a fully non-stationary site-based stochastic model without combining it with a deterministic model. Improvements include a new band-pass filter with upper- and lower-frequency limits, which refines the representation of the low-frequency content. Moreover, a time-modulating function that can represent energy arrival in multiple strong phases is introduced. The reference model’s parameters are identified by fitting to the energy content, zero-level crossings, and cumulative counts of positive-minima and negative-maxima of a target accelerogram. This fitting procedure is modified to address the increased number of parameters. These improvements broaden the reference model’s applicability while preserving its simplicity, a key aspect appealing to engineering practitioners. The improved model’s applicability is demonstrated by simulating a dataset from the February 2023 Türkiye earthquake, and the accuracy is tested using a pulse-like Next Generation Attenuation Relationships for Western United States dataset. Validations are performed based on total energy, zero-level crossings, Fourier amplitude spectrum, elastic response spectra, and peak ground motion parameters. Validations are performed schematically in the time and frequency domains and quantitatively using goodness-of-fit scores, various validation-metrics errors, and inter-period correlations. Overall, the improved stochastic model can effectively simulate a set of diverse ground motion recordings, including near-fault pulse-like records, records with multiple strong phases, and far-field motions across a broad frequency range.

Kahramanmaraş↗
USGS2026

Evaluation of stream capture related to groundwater pumping, middle Humboldt River Basin, Nevada

Historical, future, and potential stream capture from groundwater pumping in the middle Humboldt River Basin (MHRB), Nevada, is estimated using a calibrated numerical groundwater flow model. The model was developed to estimate (1) stream capture, which is the change in flux between the groundwater system and the Humboldt River and tributaries, and (2) change in streamflow, which is the change in streamflow estimated for the Imlay gage on the Humboldt River (U.S. Geological Survey streamgage 10333000). Historical stream capture for water years (WYs) 1961–2015 is estimated using recorded and estimated groundwater pumping during that period. Future (predictive) stream capture was based on historical stresses (WYs 1961–2015) using a scenario that simulated non-mine pumping from WY 2015 at a uniform rate for 100 years into the future. Potential stream capture throughout the middle Humboldt River Basin from groundwater pumping during varying durations of time are presented in a series of capture maps. Maps also are presented that show the potential to capture from groundwater evapotranspiration, as well as the storage changes for pumping duration of 100 years. Estimates of historical stream capture from the mainstem Humboldt River during the early 1960s are less than 400 acre-feet per year (acre-ft/yr) when groundwater withdrawals and pumping rates were relatively small compared to more recent times. In the late 1980s and early 1990s, groundwater withdrawals increased and estimated historical stream capture also increased from about 4,000 acre-ft/yr in the late 1980s and early 1990s to as much as 18,800 acre-feet (acre-ft) in WY 1998. In WY 2015, estimated historical stream capture declined to about 13,000 acre-ft because of decreasing groundwater withdrawals and lower streamflow during the drought of WYs 2012–15, resulting in less stream water available for capture. Stream capture was estimated for 100 years into the future based on WY 2015 non-mine pumping rates and mine-dewatering activity through WY 2015. Stream capture is forecast to increase to about 23,000 acre-ft/yr, and streamflow in the Humboldt River could decrease by as much as 19,000 acre-ft/yr. Pumping for mine-dewatering and the associated discharge of that water affects streamflow in the Humboldt River at Imlay, Nevada (U.S. Geological Survey streamgage 10333000). Historically, from WYs 1991 to 2015, streamflow was greater at Imlay gage during active mine-dewatering from mine-water discharge operations and increased by as much as 105,000 acre-ft in WY 1998. The increase was attributed mostly to the discharge of groundwater from mine-related dewatering operations directly into the mainstem Humboldt River or its tributaries, with some of this increase associated with return flows from discharge to rapid infiltration basins. Results indicate that streamflow at Imlay gage is expected to decrease by as much as 1,600 acre-ft/yr 30 years after mine-related pumping and discharge are discontinued. The streamflow reductions at the Imlay gage are expected to then decrease to around 500 acre-ft/yr, 100 years after mine-related pumping and discharge are discontinued. Potential capture maps were produced for pumping durations of 10, 25, 50, and 100 years. Capture map results indicate that areas of greater potential stream capture occur adjacent to the Humboldt River and for upstream tributaries areas north of the Humboldt River.

Nevada↗
USGS2024

Upper crustal seismic velocity structure of the Hayward fault zone, San Francisco Bay, California, USA: Results from the 2016 East Bay Seismic Experiment (EBSI-16)

We developed Vp, Vs, Vp / Vs ratio, and Poisson’s ratio models of the uppermost crust (<4 km depth) from the eastern San Francisco (SF) Bay (California, USA) to near the Calaveras fault along a 15-km-long, linear profile. Upper crustal velocities are highly variable beneath, west, and well east of the Hayward fault. We observe eight notable features, from west to east: (1) Near San Francisco Bay, there is an ~2-km-wide structure with high Vp / Vs ratios (up to 5) and Poisson’s ratios (up to 0.48) extending from the surface to the base of our model, which we suggest the structure is a near-vertical fault that lies along a straight-line projection between the Silver Creek fault to the south and the Point Richmond fault to the north. The structure may be part of an ~90-km-long fault along the eastern SF Bay. (2) The western East Bay Plain, the lower lying area between the bay and the hills, includes up to 800 m of low-velocity sediments ( Vp ~1600–3000 m/s, Vs ~500 m/s to ~1000 m/s), underlain by higher velocity basement rocks ( Vp ~3000–5800 m/s; Vs ~1000–1500 m/s). (3) Between ~1 km and 3 km east of the Bay shoreline, sediments thin in a series of steps (likely faults) toward the Hayward fault. (4) Between ~3 km west and ~1 km east of the Hayward fault (at the East Chabot fault) at depths greater than 1 km, basement Vp (up to 6000 m/s) and Vs (up to 2800 m/s) are high, and Vp / Vs ratios (<2) and Poisson’s ratios (<0.3) are low, suggesting crystalline rocks. Furthermore, a near-vertical zone of low Vp / Vs ratios and Poisson’s ratios is between near-surface traces of the Hayward and East Chabot faults, likely corresponding to the San Leandro Gabbro of Ponce et al. (2003) . (5) Eastward of the East Chabot fault in the upper 1.5 km, basement Vp (~3000 m/s to ~4200 m/s) and Vs (~1200–2000 m/s) are lower than those west of the fault. (6) In the eastern Hayward/Oakland Hills, there are zones of laterally varying, high- and low-velocity ( Vp ~2500–3000 m/s) Jurassic–Cretaceous and Tertiary sediments in the shallow subsurface that likely extend much deeper than imaged. (7) Seismic energy that propagates westward from sources east of the Hayward fault (HF) appear weaker than energy that propagates eastward from sources west of the HF, suggesting that the HF acts as a partial barrier to shallow seismic energy propagation into the more populated eastern SF Bay area. (8) Unlike many fault zones, it appears that the active trace of the Hayward fault (in our study area) is not cored by a prominent, low-velocity zone relative to rocks to the east and west of the active trace. However, the active trace does mark a prominent change from relatively higher velocities to the west and lower velocities to the east.

California↗
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↗
USGS2026

Revisiting the utility of regional-scale, high-quality geophysical data in mineral exploration - A case study featuring the Mammoth Magnetic Anomaly, Pinal County, Arizona

Regional aeromagnetic surveys passively measure the total magnetic intensity (TMI) and are a foundational tool used in mineral exploration (Airo, 2015). With the increased global demand and the number of critical mineral resources required for manufacturing high-tech devices, developing high-quality, regional-scale geophysical surveys could aid critical mineral exploration efforts and geologic mapping. In 2019, the U. S. Geological Survey launched the Earth Mapping Resources Initiative (Earth MRI) to modernize the geologic and geophysical mapping of regions that have the potential to contain critical mineral resources within the United States. In support of planning Earth MRI geophysical surveys, Drenth and Grauch (2019) defined five aeromagnetic data quality rankings (rank 1 through rank 5) applying them to the airborne geophysical survey inventory of the United States (Johnson et al., 2021). Rank 1 aeromagnetic surveys are of the highest quality, meeting modern standards and allowing best practices for qualitative and quantitative interpretation; whereas rank 5 aeromagnetic surveys are of the lowest quality, being useful only for qualitative interpretation of broad features. Through the Earth MRI effort, 48 high-quality, regional-scale rank 1 and 2 airborne magnetic and radiometric geophysical surveys have been planned, collected, or publicly release through May 2025 (U. S. Geological Survey, 2025). Here, a portion of a rank 1 Earth MRI aeromagnetic survey in southeast Arizona is presented and compared to a legacy rank 5 aeromagnetic survey over the Mammoth Magnetic Anomaly (MMA), demonstrating how modern, high-quality aeromagnetic data improves our view of crustal geology, aiding mineral exploration.

Arizona↗
USGS2025

lasertram: A Python library for time resolved analysis of laser ablation inductively coupled plasma mass spectrometry data

Laser ablation inductively coupled plasma mass spectrometry (LA-ICP-MS) data has a wide variety of uses in the geosciences for in-situ chemical analysis of complex natural materials. Improvements to instrument capabilities and operating software have drastically reduced the time required to generate large volumes of data relative to previous methodologies. Raw data from LA-ICP-MS, however, is in counts per unit time (typically counts per second), not elemental concentrations and converting these count ratesto concentrations requires additional processing. For complex materials where the ablated volume may contain a range of material compositions, a moderate amount of user input is also required if appropriate concentrations are to be accurately calculated. In geologic materials such as glasses and minerals that potentially have numerous heterogeneities (e.g., microlites or other inclusions) within them, this is typically determiningwhether the total ablation signal should be filtered to remove these heterogeneities. This necessitates that the LA-ICP-MS data processing pipeline is one that is not automated, but is also designed to enable rapid and efficient processing of large volumes of data. Here we introduce , a Python library for the time resolved analysis of LA-ICP-MS data. We outline its mathematical theory, code structure, and provide an example of how it can be used to provide the time resolved analysis necessitated by LA-ICP-MS data of complex geologic materials. Throughout the pipeline we show how metadata and data are incrementally added to the objects created such that virtually any aspect of an experiment may be interrogated and its quality assessed. We also show, that when combined with other Python libraries for building graphical user interfaces, it can be utilized outside of a pure scripting environment. can be found at https://doi.org/10.5066/P1DZUR3Z

Applied Computing and Geosciences↗
USGS1986

Water-quality appraisal of NASQAN stations below impoundments, eastern Tennessee

The National Stream Quality Accounting Network (NASQAN) is a network of stations at which systematic and continuing water-quality data are collected. Major objectives of this U.S. Geological Survey program are (1) to depict areal variability of streamflow and water-quality conditions nationwide on a year-by-year basis and (2) to detect long-term changes in streamflow and stream quality. Several NASQAN stations in East Tennessee are downstream from impoundments which have a significant effect on water quality. NASQAN data obtained from the Tennessee River below Watts Bar Dam and the Clinch River below Melton Hill Dam were compared to water-quality data from the basins upstream. The comparison indicates that NASQAN data obtained below impoundment may not be adequate to describe a composite picture of water quality in the accounting unit. Detention time of storage in the impoundments is believed to moderate the range of constituent values observed at the NASQAN stations. Data obtained upstream and downstream from Watts Bar Dam indicate that the water sampled at the NASQAN station comes from stratified layers of the impoundment and is not representative of an integrated sample of water from the impoundment. Values of total recoverable iron suggest that, because of adsorption to sediments in impoundments, some constituents are not accurately described by sampling below impoundments. Relations between water-quality constituents and flow at stations on the Clinch River and Tennessee River are not well defined due to regulation. Direct load computations for many constituents were therefore not possible, which diminished the utility of data from these NASQAN stations to account for quantity versus quality of the water. Load computations were only possible for ionic constituents through use of a continuous specific-conductance record as an intermediary. Compensation for the effects of discharge prior to application of the Seasonal Kendall test for trends could not be done and identification of trends in water-quality constituents caused by some process (source) change was not possible. Some water-quality trends indicated by data from the Clinch and Tennessee Rivers might reflect the decreasing trend in discharge during the 1972-82 water years. Thus the stations below Watts Bar Dam and below Melton Hill Dam do not adequately meet the NASQAN objective to detect and assess long-term changes in stream quality.

Tennessee↗
USGS1976

Neutron activation analysis in hydrology

Neutron activation analysis in hydrogeochemical and related investigations is increasingly practical because of developments in radiation detectors and computerized data processing. While neutron activation is the most widely used kind of activation analysis other kinds (for example, charged particle activation) are available for special applications. Neutron activation analysis is promising because approximately 70 percent of the elements in the periodic table are detected with satisfactory sensitivity, contamination from reagents is eliminated, simultaneous determination of more than 20 elements is often possible, technique is versatile, is applicable to small samples and is non-destructive in the instrumental form. Developments in theU.S. Geological Survey include a special sampling technique for water; a comprehensive technique using carriers precipitated with sulfide to concentrate desired trace elements; a determination of mercury, silver, and gold; determination of aluminum, vanadium and nine other elements; determinations of two groups of rare earth elements; and adaptations of the analytical methods to sediment and biota.

Open-File Report↗
USGS2026

Characterizing changes in postfire debris-flow hazard as burned areas recover

Emergency assessments of postfire debris-flow hazards that are performed by the U.S. Geological Survey (USGS) provide estimates of debris-flow likelihood and rainfall triggering conditions that are used for evaluating and managing runoff-generated debris-flow hazards in recently burned areas throughout the western United States. Although the immediate postfire period, within roughly one year after fire, is typically the most susceptible to runoff-generated debris flows, the hazard evolves in time and space as the burned area recovers. The recovery trajectory a given burned area will take depends on local climate and weather and can be difficult to predict. Some burned areas recover quickly, whereas others experience debris flows for multiple years after fire. As a result, extending our ability to update debris-flow likelihood estimates and rainfall thresholds based on observed recovery of the burned area would be beneficial. We present a method for multi-year runoff-generated debris-flow hazard assessment that leverages the USGS “M1” debris-flow likelihood model and integrates updated, satellite-derived, normalized burn ratio data to estimate vegetation recovery. We predict recovery-aware rainfall thresholds and validate them against a multi-year debris-flow hazard prediction and could be adapted for use with other debris-flow models that incorporate burn severity data.

Arizona, California, Colorado, New Mexico, Washing↗
USGS2026

Structural evolution of iron coordination in proteins across Earth’s oxygenation history

Protein metal-binding sites support essential biological functions shaped by protein fold, subunit interactions, and cofactor chemistry. Because these sites encode both biochemical function and environmental constraint, they offer a route to connect protein evolution with changes in Earth’s surface environment through time. Of particular interest is iron (Fe), the most widely used metal in biology and a cofactor central to both anaerobic and aerobic metabolism. Here, we systematically compare the immediate chemical environments of functionally essential Fe-binding sites in three-dimensional protein structures to test whether Fe coordination spheres differ across oxygen contexts. Using a curated dataset of experimentally determined structures, we identify a clear shift in the local chemistry of Fe-binding environments from anaerobic to aerobic proteins. Aerobic Fe sites are significantly more hydrophilic than anaerobic ones, and amino-acid composition analyses show reduced cysteine use in aerobic Fe-binding neighborhoods. These patterns suggest that as Earth’s surface environments became more oxygenated, proteins retained Fe as a core redox metal while reconfiguring local coordination chemistry in ways less vulnerable to oxidative damage. More broadly, this study introduces and applies the Coordination Sphere Analysis and Comparison (CSAC) workflow, an open and archived Python workflow for extracting local metal-binding environments from structure datasets, providing a framework for linking metalloprotein structure to evolutionary and geobiological transitions across Earth history.

Discover Life↗
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↗
USGS2024

Thermo-hydrologic processes governing supra-permafrost talik dynamics in discontinuous permafrost near Umiujaq (Québec, Canada)

Widespread supra-permafrost talik formation is currently recognized as a critical mechanism that could accelerate permafrost thaw in the Arctic (e.g., Connon et al. 2018; Farquharson et al. 2022). However, the trajectory of permafrost dynamics following talik formation may prove difficult to predict. Physically-based cryohydrogeologic models provide a powerful tool for understanding processes and factors controlling talik dynamics and, ultimately, how permafrost will respond to climate change. Such models are typically used to represent multiple non-linear processes relevant for groundwater systems in cold regions, such as coupled heat and groundwater movement, including freeze-thaw dynamics and the effects on the surface energy balance and the subsurface thermal and hydraulic properties (Lamontagne-Hallé et al. 2020). Though cryohydrogeologic modeling advances have been made in simulating talik dynamics, few applications have been tested against robust long-term hydrometeorological and subsurface observations.

Quebec↗
USGS1976

The water quality of Livingston Reservoir on the Trinity River, Southeastern Texas

The concentrations of dissolved solids, chloride, and sulfate in Livingstone Reservoir on the Trinity River in southeastern Texas usually average less than 250 mg/1 (milligrams per litre), 40 mg/1, and 50 mg/1, respectively. The water is usually hard or moderately hard (61 to 180 mg/1 as calcium carbonate). The concentrations of principal dissolved constituents in the reservoir are usually maximum during summer and fall when evaporation is high and inflow is low. Thermal stratification of the reservoir usually begins in March and persists until September or October. Neither the seasonal variation of dissolved constituents in inflow to the reservoir nor thermal stratification has resulted in significant stratification of the principal dissolved constituents. However, thermal stratification has resulted in significant seasonal and areal variations of dissolved oxygen, which results in higher concentration of dissolved iron, dissolved manganese, total phosphorus, and total inorganic nitrogen.

Texas↗
USGS2025

A probabilistic assessment methodology for the evaluation of geologic energy storage capacity—Natural gas storage in depleted hydrocarbon reservoirs

The need for energy storage, particularly underground, where capacity and duration may far exceed battery storage technologies, is especially relevant given the increasing demands for reliable power alongside the development of intermittent renewable electricity sources. Geologic energy storage facilities already exist, and expanded use would enable storing gases such as methane and hydrogen. In 2018, a National Academies of Sciences, Engineering, and Medicine report, “Future Directions for the U.S. Geological Survey's Energy Resources Program,” recommended that the U.S. Geological Survey (USGS) prioritize assessing underground energy storage in geologic formations in the United States. The U.S. Geological Survey has since developed a methodology for assessing natural gas storage capacities in depleted hydrocarbon reservoirs on a national scale. The methodology introduced in this report prescribes three approaches for calculating gas storage capacity. This methodology relies on the availability of input data, including cumulative hydrocarbon production records, reservoir petrophysical properties, and reservoir pressure data. Assessment inputs can be obtained from public, State-level databases and propriety national-scale databases, although the use of analogs could be warranted for estimating input parameters. Probabilistic assessment results are aggregated to play, petroleum province, regional, and national scales. The steps defined in this report are demonstrated on the Michigan Basin Province, which includes the Mississippian Sandstone Gas Play and the Clinton Structural Play. This methodology could be used to systematically and consistently assess hydrocarbon plays and provinces for natural gas storage capacity across the United States.

Michigan↗