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 1,207 records · Page 67Linked to original sources

Dog attacks on wild desert tortoises: A risk model

We retrospectively compiled a database of 6,727 live, wild Agassiz's desert tortoises ( Gopherus agassizii ) and evaluated them for clinical signs of trauma to shell and limbs at 50 sites in the Mojave and Colorado (western Sonoran) deserts of California, USA, spanning the years 1977–2006. Our objectives were to 1) identify tortoises with severe trauma to shell, limbs, and gular horns typically seen in attacks from dogs ( Canis familiaris ); 2) identify locations where severe injuries occurred; and 3) develop a risk model based on distances of tortoises from settlements, towns, or cities. Our models identified multiple variables of importance for tortoises with severe damage to shells, limbs, and gular horns: relative age and sex of tortoises, decades of occurrence, and location. Females and very old tortoises were more vulnerable to attacks than other tortoises. In the decades between the 1970s and 2000s, the risk of severe overall trauma to shell and limbs increased 4 times and to gular horns 16.5 times. Compared to previous decades, by the early 2000s the percent of tortoises with severe trauma increased exponentially the closer a tortoise site was to a settlement; the exponential increase began at approximately 12 km from a settlement. We suggest that the risks may be higher now because of the growth of human populations within the geographic range of the tortoise. The threats to tortoises from dogs are based on whether dogs are off-leash in the Mojave and Colorado deserts.

Arizona, California, Nevada, Utah↗

Spatial and temporal geochemical variations of lava flows and tephra deposits from the December 2020 to September 2024 eruptions of Kīlauea volcano

Kīlauea volcano underwent dramatic morphological changes in 2018. That year recorded the end of the 35-year-long eruption of Puʻuʻōʻō (1983–2018) and 10-year-long (2008–2018) Halemaʻumaʻu lava lake and emplacement of the ~4-month-long lower East Rift Zone lava flows that coincided with ~500 m of summit caldera collapse. Starting on December 20, 2020, eruptions resumed at Kīlauea’s summit. There were five summit eruptions between December 2020 and September 2023, which ranged in duration from more than a year to as short as a week. Following these summit eruptions, seismicity and deformation increased in the upper Southwest Rift Zone in 2024, culminating in a ~8.5-h-long eruption in this region on June 3, 2024. Increased seismicity and deformation then shifted to the upper and middle East Rift Zone and after several months culminated in an eruption just west of, and within, Nāpau Crater in the middle East Rift Zone from September 15 to 20, 2024. Despite vast morphological changes at Kīlauea’s summit, the geochemical compositions (i.e., whole rock and glass) that erupted from December 2020 to September 2023 are all remarkably similar to each other. Whole-rock compositions appear distinct from the preceding 2008–2018 Halemaʻumaʻu lava lake and phase 3 (i.e., summit or uprift-derived mafic lavas) of the 2018 lower East Rift Zone lava flows, although glass compositions appear to have more overlap with 2018 lower East Rift Zone glasses. The June 3, 2024, upper Southwest Rift Zone spatter and lava flows exhibit a dramatic enrichment in whole-rock MgO that is not recorded in glass, which reflects accumulation of olivine (e.g., antecrysts or xenocrysts) during dike emplacement, and is consistent with the abundance of olivine in the lava flows (5–10%). June 2024 Southwest Rift Zone whole-rock and glass compositions overlap with those erupted at the summit from December 2020 to September 2023, whereas some whole-rock trace (i.e., Sc, Sr, and Zr) and major elements (i.e., CaO) are suggestive of mixing with a magmatic component that had fractionated plagioclase and pyroxene and/or a new parental magma influencing the summit reservoir system. The September 15–20, 2024, eruption at Nāpau Crater in the middle East Rift Zone involved the most differentiated magma since eruptive activity resumed in December 2020, with its magma fractionating olivine + plagioclase + pyroxene. The September 15–20, 2024, composition resembles Puʻuʻōʻō lava flows that erupted in, or near, Nāpau Crater in 1983 (episode 1), 1997 (episode 54), and 2011 (episode 59), with episode 59 having a compositional cluster that is most similar to that of the September 2024 lava flows. The data presented and provided herein open new research perspectives for long-term analyses of geochemical variations following caldera collapse at Kīlauea volcano and facilitate comparisons with other basaltic caldera systems worldwide.

Hawaii↗

Optimizing cryo-focused pyrolysis GC/MS for tracing soil organic matter across diverse ecosystems

The cycling of organic matter in terrestrial soils and sediments is central to a range of biogeochemical processes that regulate nutrient cycling, crop productivity, trace gas emissions, and contaminant transport. Pyrolysis-gas chromatography/mass spectrometry (py-GC/MS) is a powerful tool for characterizing bulk soil organic matter (SOM) at the molecular level. In this study, we used a cryo-focused py-GC/MS system to analyze soil samples from seven diverse ecosystems: vernal pool, prairie pothole, temperate forest, tropical forest, tundra, wildfire-affected boreal forest, and grassland. We addressed a key bottleneck in molecular-level SOM characterization by developing an automated data analysis pipeline to optimize py-GC/MS and complementary evolved gas analysis/mass spectrometry (EGA/MS) methods, incorporating advanced tools for data deconvolution, developing a custom compound class library, and implementing fragmentation spectrum-based molecular networking for the first time. This improved workflow was applied to soil samples from all seven ecosystems, including multiple depths and density fractions. Our findings demonstrate that ecosystem type plays a dominant role in shaping compositional differences in SOM. We also identified trends in the source of SOM compounds (e.g., microbial vs. plant-derived) across soil depth and density fractions, which are critical for understanding persistence and turnover of SOM. Our molecular networking analysis indicated that although many compounds are widespread across ecosystems, others are restricted to specific environments, such as wetlands. This underscores the utility of molecular-level data in elucidating the complexity of SOM composition and the environmental drivers that shape it. Such molecular-level insights can deepen our knowledge of biogeochemical SOM cycles.

Environmental Science and Technology↗

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

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

Geophysical Research Letters↗

Depositional interplay between the Ancestral Rocky Mountains and Ouachita–Marathon–Sonora orogenies: Insights from provenance records in the late Paleozoic Marfa Basin, West Texas, U.S.A.

The Marfa Basin in West Texas is a late Palaeozoic synorogenic depocenter associated with regional deformation linked to the Ancestral Rocky Mountains (ARM) and Ouachita–Marathon–Sonora (OMS) orogenies in southwestern Laurentia. Basin strata range in age from Middle Pennsylvanian to the middle Permian and include the Cieneguita, Alta, Pinto Canyon, Rose Mine and Mina Grande Formations. Sandstone petrography and detrital zircon (DZ) U–Pb and (U–Th)/He double dating data from these strata reveal three tectonically driven sedimentation stages: syntectonic ARM deposition, progressive OMS foredeep deposition and an orogenic transition. The Cieneguita and lower part of the Alta Formations exhibit a Mesoproterozoic DZ age signature (~1318 and ~1076 Ma age peaks) and quartzo-feldspathic sandstone compositions sourced from the adjacent ARM-related Diablo Platform basement uplift in the Middle Pennsylvanian to earliest Permian. In contrast, the upper part of the Alta Formation, as well as the Pinto Canyon and Rose Mine Formations, have peri-Gondwanan DZ age signatures, with Mesoproterozoic (~1069–1036 Ma age peaks), Neoproterozoic–Cambrian (~700–490 Ma) and Palaeozoic (~490–300 Ma) age modes and litho-quartzose sandstone compositions derived from the OMS fold-and-thrust belt and orogenic hinterland during the early to middle Permian. The lower to middle parts of the Alta Formation have alternating DZ age signatures and sandstone compositions from both ARM and OMS sources, revealing that the transition in the sediment supply occurred during the middle Wolfcampian. This transition was not characterised by source mixing, but rather by sediment interfingering alternately sourced from the Diablo Platform uplift and the advancing OMS belt. These observations are confirmed by the DZ He ages, which reveal distinct cooling histories for both source terranes. These results document a switch from ARM- to OMS-related syntectonic deposition in southwestern Laurentia during the early Permian, demonstrating that ARM-driven deformation largely preceded the continental collision along the Marathon segment of the OMS orogen.

Texas↗

Invasion resistance varies by fuel break type in sagebrush ecosystems

Background Wildfire is an increasingly important driver of changes within sagebrush (Artemisia spp. L.) ecosystems of the western USA, often resulting in increased spread of exotic annual grasses, such as cheatgrass (Bromus tectorum L.), and subsequent losses of native vegetation and wildlife habitat. Fuel breaks— areas of land treated to reduce or redistribute fuel loads — are widely implemented to help prevent the spread of wildfires and provide areas to facilitate firefighting efforts. However, localized installation and maintenance of fuel breaks directly reduce or remove vegetation and may propagate the spread of exotic annual grasses into fuel break boundaries and surrounding areas, inadvertently weakening ecological resilience to disturbance. To investigate if exotic annual grass cover was associated with mowed or green strip fuel breaks across the sagebrush biome, we combined multiple data sources and methodologies. We used targeted field surveys and land-management agency monitoring data within a space-for-time substitution framework coupled with a progressive-change before-after control-impact (PC BACI) study design using historical remotely sensed vegetation cover data which allowed us to account for potential confounding effects of roads on annual grass cover. Results Models using both field collected and remotely sensed vegetation indices estimated increases in exotic annual grass cover over time following mowed fuel break installation, and higher exotic annual grass cover closer to mowed fuel breaks. These increases in exotic annual grass occurred within, at 500 m and at 1000 m from mowed fuel breaks. However, we found variable patterns of exotic annual grass after green strip fuel break installation depending on the data source. No increase in exotic annual grass were indicated by either analysis at distances greater than 500 m from green strip fuel breaks. However, our and field data analyses disagreed on the direction of the association of exotic annual grass cover and green strip fuel breaks. Conclusions Although fuel breaks are an important tool in managing wildland fire, our analysis underscores the importance of planting fire-resistant vegetation, rather than mowing alone, to reduce invasion by annual grasses within and around fuel breaks in sagebrush ecosystems. In addition, site characteristics that hinder the proliferation of exotic annual grasses could be evaluated when installing fuel breaks to minimize unintended effects of exotic annual grass on surrounding sagebrush habitat.

California, Idaho, Nevada, Oregon, Utah↗

The U.S. Geological Survey 2025 Puerto Rico and U.S. Virgin Islands time-independent earthquake rupture forecast

We present the 2025 U.S. Geological Survey Puerto Rico and U.S. Virgin Islands (PRVI) time‐independent earthquake rupture forecast (ERF), developed for the 2025 update to the National Seismic Hazard Model (NSHM) for PRVI. The updated ERF improves upon a prior model from 2003, including an expanded fault inventory with slip‐rate estimates, updated seismicity catalogs, and refined subduction zone geometries and deformation models. It applies the fault‐system inversion methodology to solve for rates of ruptures on modeled faults, adapted from the 2023 NSHM (NSHM23) for the western United States, including the first application of the inversion to model rates on a U.S. subduction interface. Off‐fault and intraslab seismicity are constrained by observed seismicity and use updated methods developed for NSHM23. Uncertainties in model components are substantial, and the ERF represents epistemic uncertainties through a comprehensive logic tree consisting of 1.7 billion logic‐tree branches combined across all sources.

Puerto Rico, U.S. Virgin Islands↗

Stress states on the eve of past earthquakes inform earthquake rupture through fault complexity along the San Andreas and San Jacinto faults

Estimating the evolving state of stress along active fault systems can provide insight into the conditions that generated past ground‐rupturing earthquakes and influenced their ability to propagate through areas of geometric complexity, such as fault branches and stepovers. We use quasi‐static forward numerical models that incorporate the 3D complex configuration of active faults in southern California to estimate shear tractions on the geometrically complex southern San Andreas and San Jacinto faults from 1000 to 1900 C.E. These tractions include interseismic accumulation of traction due to tectonic loading, viscoelastic relaxation of shear stress within the upper crust between earthquakes, and effects of other earthquakes on the fault network. We simulate ground‐rupturing earthquakes based on the along‐strike earthquake extents modeled by Scharer and Yule (2020) , assuming that stress drop is complete in each earthquake. We use Monte Carlo simulations to estimate uncertainty in evolving shear tractions due to uncertainties in earthquake timing and in upper‐crustal viscosity. Pre‐earthquake shear tractions typically do not exceed ∼2 MPa. Although ruptures with length <200 km have pre‐earthquake shear tractions that range from near zero to ∼1.75 MPa, these tractions are not less than ∼0.4 MPa for earthquakes with rupture length >200 km. Earthquakes with long (>200 km) ruptures occur only in the single‐stranded part of the system, whereas those with short (<125 km) rupture length and high pre‐earthquake shear traction occur near fault stepovers and branches. This suggests that high accumulated shear traction encourages longer rupture propagation, but may not be sufficient to overcome geometric complexities. This modeling approach informs our understanding of rupture propagation and provides estimates of fault shear tractions that are unavailable from direct measurements.

Callifornia↗

Summertime methane and carbon dioxide emission rates and associated variables from a national-scale survey of 146 reservoirs in the United States

Reservoirs are globally important sources of greenhouse gases, but the magnitude of their emissions is highly uncertain. Here we present data for 146 reservoirs from two surveys of reservoir methane and carbon dioxide emissions, one at the regional scale in the midwestern United States and one at the national scale in the conterminous United States, plus data from one reservoir in Washington and another in Puerto Rico. At all reservoirs, ebullitive and diffusive emissions and basic physiochemistry were measured at 15-70 locations during one 22 to 64-hour period during the summers of 2016-2023, with four reservoirs revisited a second time. Concomitant water chemistry measurements were also made at an index site. The dataset is comprised of two geospatial files and seven .csv files containing greenhouse gas emissions, water chemistry, morphology, and other relevant data. These data comprise the largest multi-reservoir emissions dataset ever assembled using consistent measurement methods.

conterminous United States, Puerto Rico↗

Understanding the links between functional performance and environment in cold desert ecosystems through the flagship biocrust forming moss Syntrichia caninervis

Background and aims Biocrust communities shape the soil surface in drylands, regulating important ecosystem properties. Despite this, few works address key functional attributes in these poikilohydric photosynthetic communities, including the length of metabolic activity in the field and its relationship with microclimate. We aim to disentangle the links between functional performance of a keystone biocrust species and environment in drylands, in order to better understand ecosystem functioning and services provided to the soil by these organisms. Methods We developed a 23 months monitoring experiment using very high resolution measurements of chlorophyll a fluorescence and microclimate in the field, comparing habitats with soils of different particle size composition and microclimatic conditions for the flagship biocrust moss Syntrichia caninervis on the Colorado Plateau, Utah, USA. Results While growing in soils of different composition did not impair the functional performance of the species, the sun/shade plasticity did affect key physiological attributes. Our results quantified a strategy based on avoidance of the most extreme environmental conditions, where metabolic activity only occur during 19.9% of the total time monitored. Relative humidity was the best microclimatic predictor of metabolic activity, and the annual length of metabolically active periods showed strong convergences with other drylands of the world. Conclusions These results collectively suggest a strong control of the environment over the ecosystem services provided by biocrusts. Our results include evidence for an improved understanding of the plasticity of S. caninervis to its variety of habitats, and its response to potential threats under different climate change scenarios.

Utah↗

Evaluating drivers of environmental change in a lake sediment core: Insights from spectroscopic metrics of water-extractable organic matter and stable carbon isotopes

Freshwater lakes play a critical role in the global carbon cycle by storing and transforming organic matter (OM) from both terrestrial and aquatic sources. Small lakes in northern temperate regions, despite their limited surface area, disproportionately influence regional carbon budgets. Buried sediments integrate OM inputs over time and archive ecosystem responses to natural and anthropogenic disturbances. However, the direction and magnitude of recent environmental changes on sediment carbon (C) dynamics remain poorly understood. A 23-cm core was collected from a small temperate lake in northeastern USA to evaluate sediment OM content and composition over timescales relevant to historical land-use change, damming, and recovery from acid deposition. Patterns in OM burial and source contributions were revealed via elemental and isotopic analyses of bulk OM and ultraviolet-visible spectrophotometry of water-extractable organic matter (WEOM). The optical metrics expanded observations of likely OM sources beyond the information gained from bulk carbon metrics (total carbon, δ 13 C). The aromaticity of WEOM increased downcore, which is consistent with a shift from increased terrestrial inputs during early logging and damming activity (pre ∼1920) to more microbial-derived OM in recent surficial sediments. Future applications of WEOM optical properties as complements to traditional geochemical metrics can enhance interpretations of lake ecosystem responses recorded in lake sediments to environmental perturbations in temperate lakes.

New York↗

A diatom-based quantitative sea-ice proxy for the Bering and Chukchi seas

Sea ice affects Earth's climate system on both regional and global scales. Its incorporation into climate can be used to achieve more accurate predictions of future climate. However, instrumental records of sea-ice concentration do not extend earlier than 1978. In an effort to extend this record, we constructed a proxy using the generalized additive model based on relative abundances of five easy-to-identify diatom species found in sediment samples across the Bering and Chukchi seas. Here we present the first quantitative diatom-based sea-ice proxy developed for Beringia. The developed proxy has been applied to two sediment cores in the Bering Sea ranging from 0 to 25.7 ka (HLY0204 51JPC) and 369 to 430 ka (IODP Exp 323 Site U1345) and one in the Chukchi Sea ranging from 2.7 to 10 ka (HLY0204 24JPC). The obtained reconstructions of sea-ice concentrations are similar, but not identical to previously published qualitative and nearby records based on other proxies. Because our results are quantitative, they can be incorporated into regional climate models. The proxy is publicly available as an R Shiny application (app) and can be applied to any diatom count from marine sediments in the region.

Bering Sea, Chukchi Sea↗

Laboratory measurements of rise velocity for individual, hydrate-free and hydrate-coated gas bubbles in water

Tracking methane transport becomes more complicated in the deep ocean where seafloor release of methane gas bubbles occurs at the high pressures and low temperatures conducive to hydrate formation on bubble surfaces. Gas hydrate formation can make the bubble surface rigid, preventing the dynamic interplay between bubble size, shape and rise velocity that gas-transport models commonly rely upon when using bubble size to predict bubble rise velocity. To better constrain gas-transport model predictions, we conducted controlled laboratory measurements of rise velocity, u z , for hydrate-free air, methane, and xenon bubbles and hydrate-coated xenon bubbles. Experimental results for u z were compared to predicted u z values from several published parameterizations used to study dissolution of gas bubbles rising in the ocean. For both hydrate-free and hydrate-coated gas bubbles, the McGinnis et al. (2006) parameterization provides the most accurate u z predictions.

Geochemistry, Geophysics, Geosystems↗

Sentinel-2 for chlorophyll-a water quality monitoring: A review of validation evidence and application potential

Water quality monitoring is integral to preserving the health of freshwater ecosystems, and satellite remote sensing has emerged as one monitoring method. Sentinel-2, in particular, has been valuable for water quality monitoring due to its 5-day global temporal revisit time and spatial resolution that ranges from 10 to 60 metres. Sentinel-2 can be used to measure and monitor chlorophyll-a, which historically has been used as an indicator of water quality, eutrophication and harmful algal blooms. Our goal was to review aquatic chlorophyll-a Sentinel-2 research to assess the types of validation evidence reported. Validation evidence is defined here as the set of information key to assessing algorithm performance, and include the spatial and temporal scales of satellite validation, reported in situ sampling method context information, demonstration of validation results through plots, and appropriate algorithm performance metrics. We highlight how the body of literature collectively contributes to advancing a national scale chlorophyll-a product that could support future resource management applications. Our review of 122 published studies indicated that much of the validation evidence corresponded to early stages, as defined by the NASA data maturity framework, due to a limited focus on individual lakes and limited detail on methodology for reproducibility. Prioritizing data accessibility for both in situ data and satellite workflows used in published studies; reporting methods with transparency and consistency; and using standard algorithm performance metrics could provide a consistent framework to support and enhance the utility of satellite inland water quality research. These three quality assurance mechanisms can promote effective evaluation of approaches for remote sensing of chlorophyll-a. Adopting these quality criteria could enable the integration of validation evidence from multiple studies, supporting more spatially and temporally representative products that would advance these approaches towards maturation for broader application.

International Journal of Remote Sensing↗

Modeling species co-occurrence effects to inform invasive barred owl management and recovery of the northern spotted owl

Robust estimation of wildlife populations represents a cornerstone of wildlife research and provides critical information to guide management, including identifying at-risk species, setting harvest rates, and evaluating predator and invasive species control programs. Efforts to enhance population estimation have long included influences one species may have on another, beginning with direct effects of predation on prey populations. More recently, researchers have incorporated co-occurrence effects, such as fear of a competitor, into Lotka–Volterra competition models to generate more robust wildlife population estimates. Here, we introduce two modified Lotka–Volterra competition models, which incorporate one- and two-way co-occurrence effects, to estimate populations of two competing species. Using the test case of northern spotted ( Strix occidentalis caurina ) and barred owl ( Strix varia ) populations in the Pacific Northwest region of the United States, we evaluate if these new co-occurrence models can generate more robust population estimates than previous models. We then evaluate if potential co-occurrence effects among barred and northern spotted owls are uni- or bidirectional. Lastly, we leverage the best-performing model to evaluate the degree to which a recently proposed barred owl culling program may help recover northern spotted owl populations. Our model results suggest that incorporating co-occurrence effects improves model fit compared to classical Lotka–Volterra competition models. We found strong evidence for unidirectional co-occurrence effects of barred owls on northern spotted owls, but not vice versa. Our simulations of barred owl culling suggest that barred owls would need to be culled from approximately 40% of all occupied barred owl territories each year to reverse ongoing northern spotted owl population declines.

Oregon↗

Estimation of parental abundance using hierarchical Bayesian modeling with data augmentation

Pedigree-based estimation methods leverage the fact that each offspring in a cohort is genotypically “marked” by its parents and represent a recent and promising toolset for estimating population dynamics. This includes pedigree accumulation estimators that model the “accumulation” of inferred unique parents within a given cohort to estimate parental abundance. Unlike close-kin mark-recapture approaches, which rely on intercohort comparisons, pedigree accumulation modeling can be completed solely using intracohort samples. This is particularly advantageous for semelparous species, where intercohort pairs are impossible and adult life stages can be difficult to sample without affecting their likelihood of successfully reproducing. Previous work has evaluated a range of estimators for such datasets, concluding that the non-parametric Chao estimator provides the most accurate and precise estimates for feasible levels of sampling effort. We used simulated data to evaluate an alternative estimator based on hierarchical modeling and data augmentation in a Bayesian framework. Results indicate that estimates from the hierarchical Bayesian estimator had comparable accuracy and better precision than both the previously tested Chao1 estimator and the improved iChao formulation across a range of sample sizes and sex ratios. Furthermore, the Bayesian estimator was far more robust to simulated errors in pedigree reconstruction, especially the presence of false negatives. Hierarchical Bayesian pedigree accumulation models can also provide additional insight into underlying reproductive ecology through their use of an explicit observation process, allowing for the incorporation or estimation of species- and population-specific reproductive dynamics. More broadly, the parametric nature of these models offers opportunities to efficiently pool information among datasets as well as to propagate uncertainty within more complex models.

Ecology and Evolution↗

Potential impacts of groundwater pumping on stream temperature are greatest in streams with substantial cold groundwater inflows

Groundwater pumping-induced reductions in streamflow (known as ‘streamflow depletion’) have been documented worldwide, but potential impacts of streamflow depletion on stream temperature are not well understood. Here, we use two types of models to identify potential impacts of pumping on stream temperature across the conterminous United States (CONUS) to determine which aspects of a stream's annual thermograph (thermal signatures) can be used to monitor and manage streamflow depletion impacts on stream temperature. We used long-term streamflow and stream temperature data from 30 streamgages across CONUS and surrogate models of streamflow depletion to analyse potential stream temperature impacts at each site. We compared two different stream temperature modelling approaches: (i) a process-based energy balance model and (ii) statistical regression models based on air temperature and stream discharge. We calculated a suite of thermal signatures under depleted and non-depleted conditions for each stream and found that maximum annual 7-day temperature and annual temperature range are potentially the most sensitive to streamflow depletion, with potential changes of at least 2°C at > 70% of the sites when using the process-based model. We also found that the regression-based models predicted much less sensitivity of stream temperature to streamflow depletion than the process-based model. This work provides an initial evaluation and sensitivity analysis of the potential impacts of streamflow depletion on stream temperature. We demonstrate that stream temperature may be most sensitive to pumping in streams with a high proportion of flow sourced from relatively cold groundwater inputs, and that regression-based stream temperature models may underpredict stream temperature changes caused by streamflow depletion.

conterminous United States↗

Spatial behavior of socially isolated wild pigs (Sus scrofa) following sounder removal via trapping

BACKGROUND The rapid expansion of wild pig ( Sus scrofa ) populations across North America, coupled with increased concern over disease transmission, has intensified the need for effective management strategies. Trapping is widely regarded as the most effective control method; however, trapping events often fail to capture entire sounders. The behavioral responses of untrapped individuals following partial sounder removal remain poorly understood, potentially undermining management efficiency. We evaluated the spatiotemporal movement responses of wild pigs that became socially isolated following trapping events. RESULTS We deployed GPS collars on 18 female wild pigs from multiple sounders and quantified post-trapping movement patterns using distance from trap site, step length, persistence velocity, space covered and overlap area over a 30-day period. Movement responses were highly variable among individuals, but wild pigs travelled an average of 1.2 km from the trap, with a maximum observed distance of 6.37 km. Space-use overlap was limited, and range sizes remained relatively stable. Individuals originating from sounders with a greater proportion of females moved farther from the trap, whereas wild pigs in better body condition exhibited lower movement velocities. CONCLUSION Socially isolated wild pigs generally remained close to trap sites following partial sounder removal and rarely dispersed from the area. This behavioral pattern suggests a predictable post-trapping window during which untrapped individuals remain spatially accessible. These findings provide critical empirical support for adaptive trapping strategies, indicating that follow-up removal efforts can be effectively concentrated near original trap locations to improve management efficiency and reduce the risk of population persistence or disease spread. © 2026 The Author(s). Pest Management Science published by John Wiley & Sons Ltd on behalf of Society of Chemical Industry. This article has been contributed to by U.S. Government employees and their work is in the public domain in the USA.

Alabama↗