USGS Science⌕ Search

SEARCH · USGS Science

Search USGS Science

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.

At least 523 records · Page 29Linked to original sources
USGS2026

Geologic framework, hydrostratigraphy, and ichnology of the Blanco, Payton, and Rough Hollow 7.5-minute quadrangles, Blanco, Comal, Hays, and Kendall Counties, Texas

During 2023, the U.S. Geological Survey, in cooperation with the Edwards Aquifer Authority, revised an initial characterization completed during 2015–16 of the geologic framework, hydrostratigraphy, and ichnology of the Edwards (Balcones Fault Zone) and Trinity aquifers in the Blanco, Payton, and Rough Hollow 7.5-minute quadrangles in Blanco, Comal, Hays, and Kendall Counties, Texas. The purpose of this report is to present the updated geologic framework, hydrostratigraphy, and ichnology of the Trinity and Edwards Groups in those quadrangles. Rocks exposed in the study area are of the Lower Cretaceous Trinity Group and lower part of the Fort Terrett Formation of the Lower Cretaceous Edwards Group. The faulting and fracturing in the study area are part of the Balcones Fault Zone, an extensional system of faults active during the Paleocene to middle Eocene that generally trends southwest to northeast in south-central Texas. Hydrostratigraphically, the rocks exposed in the study area contain a section of the Edwards (Balcones Fault Zone) aquifer, upper zone of the Trinity aquifer, and middle zone of the Trinity aquifer. In the study area, the only hydrostratigraphic units of the Edwards (Balcones Fault Zone) aquifer remaining are (from top to bottom) VII and VIII, which cap several hills. The mapped hydrostratigraphic units of the upper zone of the Trinity aquifer are (from top to bottom) the cavernous, Camp Bullis, upper evaporite, fossiliferous, and lower evaporite. The mapped hydrostratigraphic units of the middle zone of the Trinity aquifer are (from top to bottom) the Bulverde, Little Blanco, Twin Sisters, Doeppenschmidt, Herff Falls (where present), Rust, Honey Creek, Hensell, and Cow Creek hydrostratigraphic units.

Texas↗
USGS2026

Variability of discharge, nutrients, and Escherichia coli from tile drains in a small agricultural stream, a synoptic study in the School Branch watershed, Hendricks County, Indiana

The U.S. Geological Survey, in cooperation with the Indiana Department of Environmental Management, sampled a 1-mile stream reach of School Branch in Hendricks County, Indiana. Sampling consisted of 24 tile drain sites and surface-water sites both upstream and downstream from the tile drains. This study intended to (1) assess the variability of discharge, nutrients, and Escherichia coli ( E. coli ) among tiles within a stream reach in variable conditions; (2) define how the tiles affect the stream hydrology and water quality; and (3) determine potential management implications. Samples were collected on three dates—round 1 on June 12, 2018, after fertilizer was applied; round 2 on September 11, 2018, in fall harvest conditions; and round 3 on April 16, 2019, before fertilizer was applied—to encompass differing precipitation events, multiple seasons, and field conditions. Sites were sampled for water-quality parameters, dissolved and total nutrients (ammonia, nitrate plus nitrite, total nitrogen, orthophosphate, and total phosphorus), E. coli , and instantaneous discharge. Instantaneous nutrient loads were calculated from tiles and at the upstream and downstream surface-water sites. Nutrient concentrations and instantaneous loads varied depending on the tile and the sampling date. Tiles contributed between 26.5 and 66.0 percent of increased discharge between the surface-water sites. Tiles contributed between 34.7 and 87.4 percent and between 19.7 and 65.9 percent of the increased total nitrogen and phosphorus loads, respectively, for the downstream site relative to the upstream site. Most of the total nitrogen and total phosphorus from the tiles was in the dissolved forms (nitrate plus nitrite and orthophosphate). Surface-water samples had a nitrate plus nitrite to total nitrogen ratio of 71 to 93 percent, with an average of 82 percent; tile samples ranged from 44 to nearly 100 percent, with an average of 91 percent dissolved nitrogen. Surface-water samples had an orthophosphate to total phosphorus ratio of 30 to 65 percent, with an average of 53 percent; tile samples ranged from 32 to nearly 100 percent, with an average of 85 percent dissolved phosphorus. E. coli was detected in all but 5.79 percent of tile and surface-water samples. The E. coli results indicated the possibility of human or animal fecal contamination by way of septic systems leaking or directly linked to the tiles draining into the study area. Consequently, five microbial source tracking samples were collected from the upstream and downstream surface-water sites between August 3, 2021, and May 18, 2022. Human, canine, ruminant, and avian sources of fecal contamination were evaluated. The microbial source tracking samples indicated that human-specific sources are likely consistently contributing to elevated fecal contamination. Differences in discharge, nutrient concentrations, nutrient loads, and E. coli concentrations show that many factors affect how tiles alter the hydrology and water quality of receiving waters. The timing and severity of precipitation events, season, and the antecedent conditions of the field or watershed before sampling can affect the results. Tiles with the largest discharge rate had the most potential effect on the water quality and quantity of the stream. Targeting high-flowing tiles may affect efforts to improve overall water quality.

Indiana↗
USGS2026

Space-based monitoring of penguin diet links sea ice, food webs, and population change

Rapid environmental change is reshaping Antarctic ecosystems through shifts in temperature and sea ice dynamics, with implications for species such as the iconic Adélie penguin ( Pygoscelis adeliae ). Because Antarctic predators rely on a relatively small number of key prey resources in the region 1 whose life histories are intrinsically linked to sea ice, 2 , 3 these abiotic changes are expected to cascade through food webs. However, the scale and remoteness of Antarctica have limited efforts to link environmental change to ecological responses across the continent. Combining tools from imaging spectroscopy, stable isotope analysis, and hierarchical statistical modeling, we reconstructed Adélie penguin diet across the entirety of the species’ global range over a three-decade period (1984–2013) using satellite imagery from the Landsat program. This approach leveraged the distinct spectral properties (a generalized measure of color) of penguin guano as observed by satellites, 4 , 5 which vary according to penguin diet. 6 , 7 We found pronounced differences in diet across the continent and that year-to-year dietary shifts were strongly associated with sea ice dynamics. Higher sea ice corresponded to more fish-based diets, whereas lower sea ice was associated with a greater reliance on krill. Additionally, spatial differences in penguin diet were associated with long-term trends in penguin abundance, linking food web processes to large-scale population dynamics. This study represents the first use of satellite observations to capture trophic dynamics at continental and decadal scales and highlights how environmental change might restructure Antarctic food webs and influence the future of a key sentinel species.

Current Biology↗
USGS2026

Fatal leptospirosis in southern sea otters from Central California: Pathologic findings and detection of Leptospira interrogans

Leptospira interrogans serovar Pomona infections cause periodic outbreaks in California sea lions (CSLs; Zalophus californianus ) and sporadic deaths in phocids. However, the frequency of infection and associated health impacts remain uncharacterized in sympatric threatened southern sea otters (SSOs; Enhydra lutris nereis ), which serve as important sentinels of coastal health. Given the broad impacts of L. interrogans on other marine mammals, our objective was to screen selected SSOs for infection, determine whether leptospirosis contributes to SSO mortality, and describe leptospiral-associated lesions. A retrospective review (2005–2025) identified 19 candidate cases that underwent detailed review, including Leptospira immunohistochemistry (IHC), serology, and polymerase chain reaction (PCR), with a special focus on renal and hepatic lesions. Kidney samples were PCR-positive for 74% (14/19) of suspected cases. For eight of these, DNA sequence-based serogroup typing detected L. interrogans serogroup Pomona. Seven of the 14 PCR-positive leptospirosis cases were classified as fatal based on positive renal IHC and moderate to severe tubulointerstitial nephritis. All fatal cases had anti- L. interrogans serovar Pomona antibody titers ≥1:25,600. The remaining seven PCR-positive cases were considered nonfatal leptospirosis due to minimal and/or unrelated renal lesions and negative IHC. Nonfatal Leptospira -infected cases ranged from seronegative to low positive (1:400) for serovar Pomona. Antibody titers for Leptospira PCR-negative cases were negative. In fatal cases, gross renal changes were often inapparent or characterized by miliary white cortical foci. Renal histologic lesions included tubulointerstitial nephritis, acute tubular necrosis, and suppurative tubulitis with intratubular bacteria, along with positive IHC staining for leptospiral antigen in the lesions. Gross hepatic changes were also inapparent in fatal cases, and histologic lesions were rare, characterized in one animal by hepatocellular dissociation and in two sea otters by limited leptospiral antigen detection by IHC. Most Leptospira -infected sea otters (71%, 10/14) stranded during higher rainfall months in California, suggesting possible land-to-sea transmission from terrestrial hosts. Given these findings, and because L. interrogans serovar Pomona infections have been confirmed in sympatric CSLs and terrestrial mammals from adjacent watersheds, a focused investigation of potential marine and terrestrial disease transmission dynamics could provide new information to reduce SSO mortalities.

California↗
USGS2026

Ecological Benthic Units (EBUs): A new characterization of the global seafloor for ocean spatial planning and management

Effective management of deep-sea ecosystems and the high seas is hindered by the absence of a globally consistent framework for characterizing benthic habitats. Here we present the first global ecological classification of the seafloor, comprising 250 unique ecological benthic units (EBUs), distributed on the seafloor as nearly 700,000 EBU polygon occurrences, generated by intersecting a high-resolution geomorphic map with multivariate environmental seascapes. Using 17 million seafloor data points and 0.05° resolution biophysical datasets—including bottom temperature, dissolved oxygen, pH, carbon flux, sediment thickness, crustal age, and bottom currents—we identified 57 benthic regions across six major geomorphic groups (shelves, slopes, seamounts/guyots, spreading ridges, abyssal/hadal areas, plateaus). The resulting EBUs reveal previously unrecognized ecological gradients, quantify global patterns of benthic heterogeneity, and expose large-scale environmental vulnerabilities. Notably, we find that 95.6 million km 2 (26% of the ocean area) of abyssal seafloor lies below the carbonate compensation depth, that 4.16 million km 2 (1% of the ocean) of continental slopes intersect severe oxygen minima, and <1% of seamounts occur in seascapes most environmentally favorable to life. These insights provide a powerful basis for identifying rare habitat configurations, assessing exposure to climate-driven stressors, and prioritizing areas for high seas marine protected area planning, as well as a policy-relevant foundation for environmental impact assessment and biodiversity baseline proxies under the new United Nations High Seas Treaty.

Oceanography↗
USGS2026

Decision analysis in support of proactive planning for chronic wasting disease in Vermont, USA

Chronic wasting disease (CWD), a fatal, transmissible disease in white-tailed deer ( Odocoileus virginianus ) and related species, is spreading across North America but has not yet been detected in Vermont, United States (U.S.). The Vermont Department of Fish and Wildlife, along with partner agencies, wants to develop a proactive prevention and response plan in anticipation of the eventual detection of the disease. Between September 2023 and September 2025, staff from the U.S. Geological Survey and the University of Vermont facilitated a structured decision-making (SDM) process with seven State and Federal agencies that have jurisdiction over some aspect of CWD management in Vermont. The aim of this process was to generate and evaluate alternative response plans against a range of long-term objectives important to the agencies. To aid in the evaluation of the alternatives, we developed a linked set of models for white-tailed deer population and disease dynamics, hunter participation and health, forest health, economic consequences, and agricultural opportunities related to the actions being contemplated as part of the response plan. We evaluated over 256 different permutations of management actions and used multi-criteria decision analysis, a branch of decision analysis designed to help decision makers navigate tradeoffs among competing objectives, to summarize the performance of those alternative strategies against the desired outcomes. Proactive actions—those designed to slow the arrival of CWD to Vermont—were moderately effective, but the most important proactive action was surveillance to detect the disease early after arrival, which triggered response actions after detection. With the insights generated by the SDM process and the results of the analyses, the participating agencies were able to identify a preferred strategy and outline the elements of a proactive response plan. This report describes the SDM process, the technical details of the modeling work, and the results of the analyses. It is intended to serve as the technical basis for Vermont’s response plan.

BioRxiv↗
USGS2026

Testing the efficacy of industrial mitigation measures for caribou in the Arctic

1. Mitigation measures are commonly employed to reduce the negative effects of industrial development on wildlife but are not often evaluated for their efficacy. For example, oil fields in the Arctic typically incorporate design features intended to increase permeability for migratory, barren-ground caribou ( Rangifer tarandus) , even though there is limited empirical evidence of the effectiveness of some of these features. 2. Given expected increases in energy development in the North American Arctic, we examined whether two mitigation measures commonly used for migratory caribou, elevating pipelines and separating roads and pipelines, were effective at increasing the probability caribou would cross infrastructure. 3. We conducted our investigation on adult female caribou in the Central Arctic Herd of Alaska during summer, analyzing movement data from telemetry collars (2015-2020) in conjunction with spatial data on oil field infrastructure. To evaluate whether caribou would cross infrastructure as a function of the mitigation measures, we employed a generalized additive modeling framework capable of detecting non-linear and threshold responses. 4. We found that caribou were more likely to cross a pipeline when the nearest pipeline was elevated (≥1.9-m), a result that supports current mitigation recommendations. We also found that caribou were more likely to cross both a road and pipeline when they were directly adjacent to one another, as opposed to being spatially separated, a result that contradicts recommended mitigation strategies. 5. Synthesis and applications . As new energy projects are designed and implemented in environments around the globe, it is important to ensure that mitigation efforts for wildlife are scientifically validated for their efficacy. In the North American Arctic, such efforts will be critical for minimizing the impacts of expanding industrial development on migratory caribou and on the human communities that rely on them for subsistence.

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

A multiproxy reconstruction of changes in sea ice and primary productivity at IODP Site U1339 (Umnak Plateau, Bering Sea) during Marine Isotope Stage 11

The extreme warmth of Marine Isotope Stage (MIS) 11 has been extensively studied in the terrestrial realm; however, less is known about its expression in the North Pacific. Here, we present a multiproxy record from Integrated Ocean Drilling Program (IODP) Site U1339 in the southeastern Bering Sea (Umnak Plateau). We use sedimentology, diatoms, and organic geochemistry to reconstruct sea ice and primary productivity from the end of MIS 12 to the beginning of MIS 10 (430–368 ka) at this site. In late MIS 12, the Umnak Plateau experienced extensive seasonal sea ice cover, with spring sea ice concentration near 100%. Based on diatoms and ice-rafted debris, sea ice declined to year-round ice-free conditions during deglaciation but readvanced during the interglacial climatic optimum. In contrast, sea ice has not been present at the Umnak Plateau throughout the Holocene. During MIS 11, sea ice at this site may have persisted due to a more easterly and/or weaker Aleutian Low. Laminated sediments indicate enhanced seasonal productivity during deglaciation. However, productivity does not remain high throughout the interglacial. The decline in productivity during MIS 11 occurs independently of fluctuations in sea ice and upwelling and may be due to increasing nutrient limitation.

Alaska↗
USGS2026

Estimating basal area change by tree size with Sentinel-2 imagery following four fires in California, USA

Background Failure to account for tree size when estimating burn severity may not accurately capture post-fire tree mortality and post-fire forest structure. Aims We explored whether basal area mortality by tree size class could be determined from remotely-sensed burn severity indices based solely on Sentinel-2 satellite imagery. Methods We used data collected in four large California wildfires to model the relationship between proportional basal area mortality and burn severity indices derived from Sentinel-2 imagery for three tree diameter class thresholds: small (15 to 30 cm), medium (30 to 50 cm) and large (>50 cm). Key results Our models showed that for a given burn severity index value, the proportion of mortality was greater overall in smaller trees, and that the proportion of mortality in large trees changed more slowly than that of smaller trees with changing burn severity index values. Conclusions We found that models that accounted for tree size can more precisely estimate changes in forest size structure than a similar model that did not account for tree size. Implications Explicitly accounting for tree size can improve estimates of post-fire forest structure, including for large trees which make up the bulk of stand biomass and post-fire seed sources.

California↗
USGS2026

Hydrogeologic conceptualization of the Arroyos subarea and Northeast Air Force Research Laboratory subarea of the Air Force Research Laboratory, Edwards Air Force Base, Antelope Valley, California, 2017–21

Research and testing within the Air Force Research Laboratory at Edwards Air Force Base in the southwestern Mojave Desert has been conducted since the 1950s. This testing requires large quantities of groundwater for cooling and operations involving hazardous constituents that have entered the groundwater. To better understand the potential movement of contaminants through the groundwater system, the U.S. Geological Survey entered into a cooperative agreement with the U.S. Air Force Civil Engineer Center with the purpose of (1) gaining a more detailed understanding of the hydraulic connection between the unconsolidated basin-fill sediment and the underlying weathered and crystalline bedrock and (2) investigating faults and their potential effect on groundwater flow. This study was specific to the Arroyos subarea and Northeast Air Force Research Laboratory subarea of the Air Force Research Laboratory. The refined hydrogeologic characterization was derived by evaluating surface geologic maps, a regional gravity model, new subsurface geophysical surveys, lithology data from boreholes, and groundwater-level data. Refinement of the hydrogeologic characterization of the study area involved developing a map of basin-fill sediment thickness, estimating the geometry and depth to bedrock, determining the presence of weathered bedrock material, and identifying the location of faults and subsurface structures. The thickness of basin-fill sediment was derived from a regional gravity model. Basin-fill sediment ranges from less than 20 feet (ft) thick to about 150 ft thick and increases substantially toward a deep structural basin in the northwestern part of the study area. The top of bedrock estimated from the new surface geophysical surveys is generally deeper than the top of bedrock from the regional gravity model. These differences in interpreted bedrock tops indicate that there may be a zone of weathered bedrock in the subsurface, overlying crystalline bedrock. Additional subsurface data evaluating the geometry and thickness of weathered bedrock are necessary for a comprehensive evaluation of weathered bedrock in the study area. Previously unmapped extensions of the Leuhman Fault, Spring Fault Zone, and Boron Gate Fault were identified in the new surface geophysical surveys. Evaluation of groundwater flow between bedrock and basin-fill sediment, and the effects of faults on groundwater flow, was done using groundwater-level data from 2020. The groundwater table was in crystalline bedrock throughout most of the Arroyos and Northeast Air Force Research Laboratory subareas, with a few exceptions where it was in basin-fill sediment and weathered bedrock. Basin-fill sediment, and likely weathered bedrock, was partially saturated in parts of the study area where basin-fill sediment was as thick as 150 ft. Weathered bedrock material may be present near the Spring Fault Zone, and if so, was likely saturated in 2020. The groundwater-level data near and within Spring Fault Zone indicate that the Spring Fault Zone acts as a vertical-along-strike conduit for, and as an across-strike impediment to, groundwater flow. In addition, steep groundwater-level gradients in the southeastern part of the study area indicate that a subsurface structure affecting groundwater flow could be related to the Leuhman Fault. Sufficient groundwater-level data were not available to conclusively determine if the Arroyos, Rich, and Boron Gate Faults, and their inferred extensions from the geophysical surveys, affect the movement or direction of groundwater flow.

California↗
USGS2026

Groundwater source water assessment area delineation of public supply wells in the Treasure Valley, western Idaho

The Idaho Department of Environmental Quality (IDEQ) is responsible for assessing the source water for public water supply systems including the delineation of source water assessment areas for wells, which are defined as the subsurface areas contributing water to a well during a period of time. Part of the IDEQ assessment plan is to update source water assessments over time. The U.S. Geological Survey (USGS) developed a new groundwater-flow model for part of the western Snake River Plain aquifer in southwestern Idaho. This groundwater- flow model provided a new opportunity for IDEQ to update source water assessment for public water supply wells within the model boundary. The USGS, in cooperation with IDEQ, used the USGS Treasure Valley Groundwater Flow Model along with a MODPATH particle tracking post-processing program to reassess public-water supply source water assessment areas. The USGS created an automated script that requests a well location, depth, and pumpage rate, and then incorporates that information into the Treasure Valley Groundwater Flow Model combined with the MODPATH model to create 3-, 6-, and 10-year time of travel zones that are output to a geographic information system shapefile. The TVSWA_Delineator script provides a fast, transparent tool for IDEQ to reassess or create new public-water supply source water assessment area delineations based on new and updated information and is available as a Python script and a Windows executable with the compiled script and all necessary components. TVSWA_Delineator utilizes a quadtree approach to rediscretize the Treasure Valley Groundwater Flow Model grid around the pumping well location from 1-mile cell lengths telescoped down to 1/16 mile. This rediscretization increases the source water assessment area delineation resolution near the pumping well and minimizes the effect of weak sinks in the MODPATH model. The stress inputs for each month of the final 10 years (2006 to 2015) of the Treasure Valley Groundwater Flow model historical simulation were averaged for the particle tracking analysis. The MODPATH model used backwards tracking of particles released at each stress period to create the source water assessment area delineations. This report documents the creation of the TVSWA_Delineator Python script, including adaptations made to the Treasure Valley Groundwater Flow Model, creation of a companion MODPATH model, and post-processing results to provide source water assessment areas.

Idaho↗
USGS2026

Groundwater and surface-water interactions in the Waihe‘e-Kahalu‘u watershed, O‘ahu, Hawai‘i—Analysis of historical data and numerical groundwater-model simulations

The Waihe‘e-Kahaluʻu watershed and surrounding areas lie in a dike-impounded groundwater setting of the Ko‘olau Range, O‘ahu, Hawai‘i, where groundwater withdrawal from wells and tunnels have resulted in reductions in streamflow. Resource managers and surface-water users seek information that can be used to balance the needs of competing uses of groundwater and streamflow in the watershed. This study assesses the effects of four groundwater developments in the Waihe‘e-Kahalu‘u watershed—Kahaluu Tunnel, Waihee Tunnel, Kahaluu Well, and Waihee Incline Wells—on Kahalu‘u and Waihe‘e Streams. Analyses of historical streamflow, withdrawal, and rainfall data indicate that base flow (the component of streamflow that is supplied from groundwater discharge) in Kahalu‘u and Waihe‘e Streams decreased in response to withdrawals from wells and tunnels within and outside the watershed. Groundwater-model simulations indicate that shutdown of selected tunnels and wells in the Waihe‘e-Kahaluʻu watershed will result in increases of groundwater discharge, but not all of the increase will occur in streams within the watershed; increases will also occur in streams outside the basin, the ocean, and tunnels other than those that were shut down. The analysis of historical data and model simulations is consistent in showing that the effects of withdrawal changes on streams are not constrained by the topographic boundaries of watersheds. Stream base flows in the Waihe‘e-Kahalu‘u watershed are subject to a complex interaction of multiple withdrawals from wells and tunnels installed at various times within and outside the watershed. The correspondence between withdrawals and base-flow changes is consistent with the conceptual model and understanding of the connection between groundwater and streamflow in the dike-impounded groundwater setting of the Waihe‘e-Kahalu‘u watershed and the Ko‘olau Range.

Hawaii↗
USGS2026

Thermal detectability of subsurface water ice on Mars: A comparative analysis for the Subsurface Water Ice Mapping (SWIM) Project

We have developed a new global map of shallow ground ice distribution, SWIM23, based on Mars Global Surveyor Thermal Emission Spectrometer data and made systematic comparisons between this new map and two similar, previously developed data products. We have explored the origin of differences between the three ice maps by detailing technical and procedural differences in their development, by making global pixel-by-pixel comparisons, and by carrying out a series of one-dimensional thermal simulations to explore fundamental physical limitations of thermal ice-detection techniques. These efforts and the production of a composite thermal ice-consistency map supported integration of multiple geophysical data products relevant to ice detection in the upper meter of the Martian regolith into the larger Mars Subsurface Water Ice Mapping project. Our work also highlights fundamental physical limitations to thermal ice detection as a technique, particularly the rapid fall-off in ice detection sensitivity at depths >30 cm, which produces maximum uncertainty in the presence and depth of ice within regions preferred for potential human landing sites. A future Mars orbiter mission designed to detect ice and support crewed landing site selection in the midlatitude region should give payload priority to an instrument capable of probing the 1–5 m depth range (i.e., a high-frequency radar), over a next-generation thermal spectrometer, which is unlikely to offer clarity on ice table depths or lateral continuity of the ice table in the locations of highest interest.

Planetary Science Journal↗
USGS2026

Holistic understanding of uncertainty for collaborative and proactive global change decision making

Global change is accelerating and pushing the planet's ecosystems beyond the range of historical observations, creating increasing uncertainty in future system conditions. Despite general agreement that proactive environmental action is warranted, environmental decision conversations often end by identifying additional data needed to reduce uncertainty before taking novel action. Given the inherent uncertainty in complex issues such as global change, quantitative data alone are likely insufficient to support proactive environmental action. Holistic understanding of uncertainty includes scientific quantification of uncertainty paired with emotional responses and transcendental grounding to help people work together toward proactive action in uncertain decision contexts. Holistic understanding arises from the four ways in which humans perceive the world, termed the Four Realms: Physical (e.g., how I observe), Mental (e.g., how I think), Emotional (e.g., how I feel), and Transcendental (e.g., how I connect to greater meaning or purpose). Environmental scientists and decision makers are generally trained in Physical and Mental Realm observation and analysis, but not in how to apply Emotional and Transcendental Realm understanding. Emotional and Transcendental processing occurs in scientists and decision makers whether it is acknowledged or not and contributes to different people interpreting the same information in different ways. Thus, when the role of Emotional and Transcendental Realms in an individual's interpretation process is not understood, it can derail conversations and perpetuate the status quo. Explicitly recognizing all Four Realms can bring people together across differences and inspire shared, novel decision making even in increasing uncertainty. To illustrate the benefits of holistic understanding, we share stories from our experiences in environmental decision contexts. Because accessing the Four Realms requires experiential and embodied techniques, while still relying on core scientific tenets of observation and analysis, we also present techniques for readers to learn to feel their own emotional understanding and connect to their own transcendental understanding. Holistic understanding can enhance data-driven decisions by recognizing that human responses to uncertainty inherently include interactions between emotions, thoughts, transcendental connections, and behavior. Ultimately, holistic understanding can help anchor data-driven decisions in intra- and interpersonal connections, inspiring action in the face of uncertainty.

Ecosphere↗
USGS2026

Digging into soil: Effects of soil texture on RT-QuIC performance for environmental prion surveillance.

Chronic wasting disease (CWD) is a fatal neurodegenerative disease caused by infectious prions affecting wild and captive cervids. Transmission occurs directly between hosts or indirectly through contact with prion-contaminated environments. Soils, particularly those rich in clay, are hypothesized to enhance prion stability, retention, and bioavailability. Accurate detection of prions is therefore important for understanding environmental transmission risks. Real-time quaking-induced conversion (RT-QuIC) is a sensitive assay used to detect PrP CWD in tissue, excreta, and environmental materials. However, RT-QuIC performance across soil textures has not been evaluated. This study assessed RT-QuIC sensitivity and specificity using laboratory-prepared soils spiked with CWD-positive brain homogenate or water controls under a standardized extraction method. Conditional on the extraction method used, results suggest that RT-QuIC performance depends on soil texture, and thus, an optimal time-to-threshold (TTT) cutoff required to balance sensitivity and specificity will also vary with soil texture. RT-QuIC exhibited higher sensitivity and moderate specificity in soils with low clay (<20%) and moderate to high silt content, whereas high-clay soils (>20%) with low to moderate silt content (2%–60%) reduced both sensitivity and specificity, and required shorter TTT cutoffs. These findings highlight the importance of accounting for soil texture in environmental CWD surveillance.

Environmental Science and Technology↗
USGS2026

Mercury mobilization and export from the Greenland Ice Sheet using an ice-to-ocean approach

The Greenland Ice Sheet (GrIS) is a poorly constrained source of mercury (Hg) to Arctic ecosystems. We measured Hg concentrations and stable isotopes along an ice-to-ocean continuum to identify controls on GrIS Hg export. Early-season permafrost melt and rainfall produced high filtered total mercury (fTHg, ~17 pM) and monomethylmercury (MMHg, ~2 pM). As subglacial drainage evolved, particulate Hg doubled (from ~8 to 17 pM) and MMHg production remained elevated, indicating Hg mobilization from subglacial environments. Shifts in Hg stable isotope ratios and Δ 199 Hg mass balance show supraglacial sources contribute 20–48% of exported Hg, suggesting subglacial inputs dominate the seasonal Hg flux. Fjord waters were enriched in fTHg ( ~ 10 pM) and MMHg ( ~ 2 pM) relative to rivers, consistent with particulate Hg transformations and terrestrial Hg inputs. The estimated GrIS Hg yield ( ~ 23 mmol km −2 yr −1 ) is similar to that of Arctic rivers and will likely increase with climate-driven mass loss.

Communications Earth & Environment↗
USGS2026

An automatable method for developing preliminary stream crossing designs using lidar and hydraulic modeling in Massachusetts

Upgrading aging and undersized stream crossings (places where culverts and bridges intersect a stream) with replacements designed for aquatic organism passage can improve infrastructure resiliency, geomorphic stability, and maintenance costs. The goal of this work is to develop preliminary culvert designs that convey hydraulic design floods and meet the Massachusetts Stream Crossings Standards for select existing stream crossing sites across Massachusetts. We derive stream and infrastructure data from lidar and other geospatial datasets as input for one-dimensional hydraulic models. Culvert dimensions are iteratively increased until the structures can convey the hydraulic design flood flows with adequate freeboard. The approach was evaluated by comparing lidar-derived designs to designs produced with surveyed data at 30 stream crossings across Massachusetts. The automatable methodology was then applied to 343 stream crossing sites in the Housatonic River watershed in southwestern Massachusetts. For each stream crossing site, preliminary designs were developed to meet the Massachusetts Stream Crossing Standards and convey the 10-, 4-, 2-, and 1-percent annual exceedance probability flood flows. The U.S. Geological Survey StreamStats web application hosts the 3-sided box, 3-sided arch, and pipe culvert preliminary designs. Provided supporting information includes site information, habitat quality assessments, and modeled water surface elevations. The publicly available preliminary designs can be used by municipalities, engineers, and others to evaluate design options, prioritize replacements, and analyze ecological benefits. These designs also may help municipalities pursue grant opportunities and estimate replacement costs prior to investing in field assessments and engineering design.

EarthArXiv↗