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 613 records · Page 34Linked to original sources
USGS2026

Avoiding common pitfalls in the implementation of adaptive forest wildlife management

Forest management planning requires the evaluation of tradeoffs in proposed management actions to meet multiple objectives, including providing the habitat required to sustain wildlife populations. However, various species-specific life history strategies and ecological uncertainties complicate our understanding of the responses of wildlife to disturbances, such as wildfires and forest management. We believe that adaptive management processes, applied through the close collaboration of researchers and managers (e.g., coproduction), will be valuable in this effort. Although adaptive management programs have been touted in the past, we stress that many of these programs have fallen short of the requisite steps to complete the learning cycles of this process (i.e., single- and double-loop learning). Thus, we provide guidance on what we believe to be the missing components of these efforts in hopes to guide future coordination of local and regional adaptive management programs for forest wildlife.

General Technical Report↗
USGS2026

A simplified two-station approach for modeling metabolism in dam tailwaters subject to diel flow variation

Tailwaters are ubiquitous and highly managed ecosystems whose food webs often rely disproportionately on autochthonous energy. In situ continuous dissolved oxygen data are increasingly being used to estimate gross primary productivity and ecosystem respiration in rivers, but this approach is complicated in tailwaters, where upriver discontinuities (i.e., dams) violate commonly employed one-station approaches. In such cases, two-station metabolism models can be applied, although substantial diel variation in flow (a common outcome of hydropower production) requires more complex treatment of water parcel travel times. Here, we present a new two-station metabolism model that allows estimation of reach-scale gross primary productivity and ecosystem respiration in streams and rivers that experience within-day variation in flow. Our approach simplifies two-station variable flow model implementation compared to previous efforts. We apply our model to a 6-yr dissolved oxygen time series and use Bayesian inference to estimate daily gross primary productivity, ecosystem respiration, and gas exchange velocity ( k 600 ) for a ~12-km reach of the Colorado River downriver of Glen Canyon Dam. We compare our model's performance to a more mechanistically detailed and computationally intensive Eulerian dynamic flow model and also to a widely-used one-station model that uses assumptions of reach uniformity that are often strongly violated in tailwaters. These comparisons show that our metabolism estimates conform with output from the more detailed dynamic flow model and that the one-station approach deviates substantially from both two-station approaches. Our new stream metabolism model can help resolve a fundamental analytical impediment in tailwater ecology.

Arizona↗
USGS2026

Sedimentological reconstruction of hurricane storm surges from a carbonate reef coastline at Anegada, British Virgin Islands

Anegada is a small, low-lying carbonate island in the British Virgin Islands that is highly vulnerable to marine inundation from frequent hurricanes. A major challenge in forecasting the timing, intensity, and impacts of future hurricanes to islands such as Anegada is the limited short-term instrumental record that does not capture the full spectrum of inundation events from the frequent but low-magnitude to rare, but high-magnitude events. Washover deposits, which are transported onshore by hurricane-related storm surges and preserved within the onshore stratigraphy, can provide insight into the recurrence of these events. Here, we use grain size distributions, X-ray fluorescence (XRF) data, and an inverse sediment transport model to characterize a series of anomalous sand deposits from shallow hypersaline ponds to determine the variability in storm surges impacting Anegada over at least the last 70 years. The analyses were successfully applied to modern analogues of storms impacting Anegada's southern shoreline (Hurricanes Donna in 1960 at Category 3; Earl in 2010 at Category 4; and Irma in 2017 at Category 5). These deposits are identified by coarse-grained material and elevated marine geochemical signals (Ca, Sr, ln(Br/Ti)). Inverse modelled flow parameters (average flow depth = 0.89–1.46 m relative to the reef flat, average flow speed = 2.95–3.78 m/s) agree with observational data and eyewitness accounts, and the deposit associated with Hurricane Donna indicates the greatest impacts occurred along the southern shoreline (modelled flow depth = 1.24–1.46 m relative to the reef flat, flow speed = 3.49–3.78 m/s). Using the modern analogues as a basis for comparison, we infer the origin of an older sand deposit to be from Hurricane Dog in 1950, when it tracked by Anegada as a Category 3 storm. This methodology provides a transferable framework for interpreting storm recurrence in carbonate-dominated coastlines worldwide that could be used to inform hurricane hazard assessments and risk management efforts.

Anegada↗
USGS2026

Mercury biomagnification across food webs with varying non-native fish presence: Implications for native fish conservation in the upper Colorado River Basin

Environmental contaminants and non-native species introductions negatively affect aquatic ecosystem conservation. Mercury (Hg) accumulates within food webs where it can biomagnify to toxic concentrations, which can be affected by altered trophic relationships from non-native species. This study examined Hg concentrations (n samples = 655) and stable carbon and nitrogen isotope ( n = 746) compositions in native ( n = 313) and non-native ( n = 342) species among four river subbasins (Gunnison, Lower Green, Colorado-Dolores, and White-Yampa), in the Upper Colorado River Basin (UCRB) during 2016 and 2017 to inform the potential combined effects of Hg exposure and non-native species interactions on native fishes of conservation concern. We assessed Hg biomagnification at the fish assemblage level, compared concentrations to risk thresholds for fish and human health, and evaluated niche overlap among trophic levels across subbasins. The generalist-invertivores feeding guild had the lowest Hg concentrations (0.001–0.191 μg g −1 wet weight), and no individuals exceeded the U.S. Environmental Protection Agency aquatic life criterion (0.225 μg g −1 ). Piscivores had the highest concentrations (0.008–1.840 μg g −1 ), and 74% of individuals exceeded the criterion. Biomagnification across trophic levels was observed throughout the UCRB, yet varied by subbasin, with higher magnification factors in the Lower Green and White-Yampa subbasins which also had the highest proportion of non-native species. Stable isotope compositions revealed niche overlap among native and non-native species. Additionally, trophic position varied within species among the subbasins. The proliferation of non-native species has likely reshaped food webs and may have intensified the potential threat Hg poses to native fish recovery and conservation. This study documented increased Hg concentrations relative to past assessments, suggesting that continued monitoring could assist in evaluating trends in Hg accumulation. Particularly studies that focus on sources entering the food web, vectors for accumulation, and competitive interactions among native and non-native species.

Colorado, Utah, Wyoming↗
USGS2026

Growth of a threatened desert fish becomes resource dependent when released from temperature limitation in a regulated river

1. Growth of ectotherms typically positively correlates with increasing temperature towards an optimal temperature range, followed by a decline in growth when temperatures exceed the optimal temperature range for that organism. When temperature is within or close to the optimal range, food availability and other environmental factors can play important roles in modifying growth rates. River regulation often alters both environmental conditions (e.g., flow, turbidity and temperature) and aquatic food webs, modifying the quality and quantity of resources available to consumers. Here we aim to understand how food availability, water temperature and other environmental factors are associated with fish growth (defined here as a change in length) in a large, regulated desert river. 2. We measured growth of humpback chub (Gila cypha) at two river reaches located 240 km apart in the Colorado River within the Grand Canyon, Arizona, USA. We then fit Bayesian state-space models of growth that account for variable time at large to test environmental predictors such as water temperature, flow, turbidity (which increases during tributary flooding and can be a proxy for allochthonous inputs) and gross primary productivity (a proxy for autochthonous food) using an inclusion parameter approach. Both river reaches are located below Glen Canyon Dam, which regulates flow, alters seasonal water temperatures and traps sediment, nutrients, and organic carbon; however, merging tributaries and a greater distance from the dam lead to differences in environmental conditions. 3. In the upriver reach, where warm season water temperatures are below the optimal range for humpback chub most of the time, temperature was the primary driver of variation in growth according to inclusion parameter values and the strength of the standardised effect size. Turbidity was also included in the final model for the upper reach, but the inclusion parameter values and standardised effect size were smaller than for temperature. In the downriver, warmer reach, temperature still had the highest inclusion parameter value, but GPP had the highest standardised effect and nearly as high an inclusion parameter value. 4. We conclude that temperature is the primary limiting factor on growth of humpback chub in the Grand Canyon; however, as water temperatures increase, other factors, especially GPP, are associated with temporal variation in growth rates. Over the past two decades, declines in reservoir elevations and the volume of summer releases have increased downstream water temperatures. If warm temperatures persist into the future, other factors may have increasing roles in regulating interannual variability in chub growth. 5. Field studies of ectotherm growth often emphasise the role of water temperature; however, the role of resource availability is also increasingly recognised, especially in bioenergetics studies. Quantifying resource availability directly can be difficult, especially in large, remote rivers. Here, we show here that using proxies for resource availability (i.e., GPP as a proxy for autochthonous inputs) can yield useful insights regarding resource limitations with implications for management of federally listed fish species in a larger regulated river system.

Arizona↗
USGS2026

Flood-frequency estimates for Kentucky streamgages based on data through water year 2021 and results of updating the fundamental layers in Kentucky StreamStats

The U.S. Geological Survey, in cooperation with the Kentucky Transportation Cabinet, analyzed flood-frequency statistics for streamgages in Kentucky. Using annual peak-flow data through water year 2021, flood-frequency estimates were computed for 261 streamgages, including unregulated and regulated sites as well as sites with mixed regulation records. Methods followed those outlined in “Guidelines for Determining Flood Flow Frequency—Bulletin 17C” (U.S. Geological Survey Techniques and Methods 4–B5). These estimates included flows corresponding to annual exceedance probabilities of 50, 20, 10, 4, 2, 1, and 0.2 percent. Temporal trend analyses using the Mann-Kendall test indicated that 18 percent of unregulated streamgages with (1) at least 30 years of peak-flow record and (2) peak-flow record at least as recent as water year 2000 showed statistically significant trends, most of which were weak to moderate increases in peak flows. Concurrently, the fundamental geospatial datasets that support the Kentucky StreamStats application were updated by using high-resolution digital elevation models and hydrography datasets to derive flow direction, flow accumulation, and stream definition rasters. Comparisons of regression-based flood-frequency models using the old and new layers demonstrated consistent results, with a statewide root-mean-square error of 0.019, in the base-10 logarithm of cubic feet per second. Furthermore, to assess model performance, flood-frequency estimates made by using the updated layers and previously published regression-based models were compared to flood-frequency estimates newly computed by following Bulletin 17C. This analysis showed the models performed adequately for most Kentucky stream locations. The updated statistics and geospatial layers provide stakeholders with more accurate, current data for flood-risk assessment, infrastructure design, and water-resource management.

Kentucky↗
USGS2026

Mass mortality of avian migrants in New Mexico, USA, that coincided with an extreme weather event

Many birds are migratory, and this life history strategy allows for maximized access to seasonally abundant resources and favorable climates. However, migration exposes birds to threats and stressors, resulting in high mortality during migration. Anthropogenic landscape alterations and climate change have intensified threats, and mass mortality events linked to extreme weather are more common in recent decades. Documenting mass mortality is critical for predicting future occurrences and implementing effective conservation. Here, we describe a mortality event that occurred throughout New Mexico, USA in fall 2020. Carcasses began appearing in the region in mid-August, during a period of extreme heat and drought. Following an extreme cold weather event on September 8–9 th , the number of carcasses increased dramatically and expanded throughout the state. In total, we collected 628 carcasses comprising 58 species within Doña Ana and Otero Counties in New Mexico. Necropsy determined emaciation was the cause of death for 74.6% of carcasses. Live birds captured during the period of peak mortality (n = 223) were in similarly poor condition. This event provides a striking example of how multiple types of extreme stressors, in this case widespread drought and unseasonal cold, coincided with a mortality event, indicating a possible synergistic relationship between these factors and the mass mortality. Mortality events are likely to increase in frequency with intensifying climate change. Establishment of networks of biologists and researchers could improve our ability to identify and communicate developing mortality events, organize data collection, and improve understanding of the causes and consequences of mortality.

New Mexico↗
USGS2026

Seasonal estimates of nutrient loading, sources, and impacts on water availability in streams across the conterminous United States, 2000–2020

Excess riverine nutrients can limit water availability for humans and ecosystems. Nutrient monitoring on individual stream reaches provides important information about current and past conditions, yet most river reaches are unmonitored. Additionally, monitored data integrate effects of upstream sources, sinks, and pathways by which nutrients get into streams, making explicit planning for effective nutrient reduction strategies challenging. We developed and applied dynamic SPAtially Referenced Regression On Watershed attributes (SPARROW) models to quantify total nitrogen (TN) and total phosphorus (TP) delivery and source-specific loading in all stream reaches of the conterminous United States (CONUS) at a seasonal timestep from 2000 through 2020. Between one third and one half of non-point source nutrients lagged more than one season in their delivery to CONUS streams over the study period. The largest sources of nutrients that were applied and delivered to streams within a season were agricultural fertilizers, treated wastewater, atmospheric deposition, and manure. Increases in loads were widespread across the CONUS and predominantly associated with increased loading from lagged delivery of non-point source inputs and current season fertilizer application. Model results also identified extensive areas where predicted seasonal concentrations frequently exceeded the National Rivers and Streams Assessment “poor” criteria. Increased loading occurred in areas with high yields and delivery to downstream receiving waterbodies, particularly across the Midwest, raising the risk of negative impacts to downstream water availability. The modeling approach presented here provides a continental-scale framework for identifying when, where, and why nutrients are generated and how that may ultimately impact receiving waterbodies.

conterminous United States↗
USGS2026

Transcriptional changes in wild Yukon River Chinook Salmon associated with Ichthyophonus infections

Objective We compared differentially expressed genes in Chinook Salmon Oncorhynchus tshawytscha with three divergent Ichthyophonus statuses (undetected, subclinical infections, or clinical disease; n = 100) to investigate associated transcriptomic responses. Disease associated with the fish parasite Ichthyophonus sp. was first diagnosed in adult Chinook Salmon from the Yukon River in the late 1980s and has subsequently been implicated in premature host mortality. Methods Chinook Salmon tissue sample collections and Ichthyophonus infection data were leveraged from a multi-agency collaboration during summer 2022 at three locations along the main-stem Yukon River that spanned nearly 2,000 km of freshwater migration (lower, middle, and upper river). We sequenced the transcriptome and compared this to infection status based on routine diagnostic procedures. Results Among the 17,569 genes for which messenger RNA was detected, we identified a transcription signature in the skeletal muscle that was associated with Ichthyophonus infections and included 53 differentially expressed genes. The differentially expressed genes and their pathways included those known for involvement in immune functions, energy synthesis, cellular breakdown, and reproduction—all processes that are known to be influenced by senescence during spawning migrations. Conclusions Results demonstrate a clear transcriptional difference between diseased fish (clinical disease group) and those in which Ichthyophonus was undetected, including identifying candidate markers for infection in this population. These results provide a foundation for development of nonlethal biomarkers to evaluate potential Ichthyophonus infections in Chinook Salmon based on gene transcription, protein products, or gene variants (e.g., polymorphisms).

Journal of Aquatic Animal Health↗
USGS2026

Review of Bayesian dynamic occupancy models and the implications of incomplete sampling

One popular tool for modeling long-term wildlife monitoring data is the dynamic occupancy model. Dynamic occupancy models require detection/non-detection observations recorded during repeated visits to the same site across many seasons. However, sampling all sites repeatedly for multi-year studies is not always possible due to limited funding, sites being difficult to access, the addition of sites after monitoring has begun, or even global pandemics. We explore the ramifications of incomplete sampling of sites across seasons on the modeling of dynamic occupancy data. We review four popular choices for modeling dynamic occupancy data, then, with simulation, we explore the implications of incomplete sampling on the inferences gleaned from the different models. These models are also fit to a data set monitoring American Bullfrog breeding populations along the Yellowstone River floodplain in Montana, USA. Results from the bullfrog data analyses are consistent with the results from the simulation, and help support our modeling conclusions. Ideally, if there are no limitations on data collection, we suggest repeatedly visiting all sites in all seasons, however this is often difficult, or impossible for a variety of reasons. For data already collected, regardless of incomplete sampling, we suggest using the dynamic or auto-logistic models.

Montana↗
USGS2026

A sea ice entrapment event in the southern Chukchi Sea: Analysis and prediction

Amplified Arctic warming is reducing sea ice cover, which is driving an increase in geopolitical interest in the region as it offers the possibility of reducing shipping times between Asia, Europe, and eastern North America, at the risk of increased ice hazards. Here, we examine the case of the Norseman II research ship that was trapped by sea ice in the southern Chukchi Sea for 14 days in June 2024. This is the first study of its sort in this region. We show that anomalously thick and extensive sea ice was present north of the region prior to the event, and that strong northerly winds in early June advected this ice southward, trapping the ship. Later in June, southerly winds advected ice northwards away from the ship, helping to free it. We further show that the event was forecastable.

Alaska↗
USGS2026

Role of volcano-tectonic interactions during early-phase magma-assisted continental rifting

We study volcanic-tectonic interactions during the early phase of rifting in the Natron Rift, a magmatic segment of the East African Rift System. We use the open-source crustal deformation modeling code, PyLith, to model deflation during the 2007–2008 explosive eruptions and inflation during 2016–2023. We compare predictions of surface displacements resulting from magma reservoir inflation from 3D models with and without topography and homogeneous and heterogeneous material properties against Global Navigation Satellite System displacements from 2016 to 2023. We prefer the model with heterogeneous material properties and topography for fitting the observed displacements. Including topography improves the fit to GNSS data and modulates stress changes from shallow magma sources. Our numerical modeling of Coulomb stress changes suggests that deflation of a magma reservoir under the active Ol Doinyo Lengai volcano could trigger oblique normal slip on the Natron Fault with Coulomb stress increases exceeding ∼0.1 MPa on the volcano-adjacent fault segment. Fault slip on the Natron Fault is less likely during inflation of the magma reservoir because normal compressive stress increases on the fault. We also demonstrate that Coulomb stress changes are sensitive to the inclusion of topography in our model, particularly because the depth of the magma reservoir is similar to the amplitude of the topographic relief.

Ol Doinyo Lenga↗
USGS2026

Geochemistry of the 2022 Mauna Loa eruption: A comparison with earlier historical summit reservoir eruptions, with implications for magma supply and recharge

On November 28th, 2022, following a record historical repose period of 38 years, Mauna Loa erupted about 145 × 10 6 m 3 of lava and tephra over a 15-day period. The eruption was confined to the summit caldera region and the upper Northeast Rift Zone and is remarkably homogeneous in composition in both time and space. In these respects, it is typical of prior shallow summit reservoir magma bodies, recently estimated to be at a depth of around 1–2 km beneath the caldera. In contrast with these earlier magma bodies, which typically contain 6.7–7.1% MgO and are perched at the low-MgO end of olivine-control trends, the 2022 lava and tephra are more evolved with 6.24 + / − 0.03% MgO. This implies a temperature difference of around 11 °C with the prior 1984 magma. The simplest explanation is that over 38 years, cooling and crystallization of the remaining 1984 magma body has significantly exceeded magma recharge, giving rise to the evolved 2022 magma. The problem with this model is that we know from a variety of geophysical observations that in those 38 years, Mauna Loa has been erratically inflating, with heightened periods since around 2000 attributed to magma recharge. To reconcile these differences, we suggest instead that the 1984 magma cooled and crystallized much more extensively, from 1166 °C to around 1106 °C, co-crystallizing plagioclase, clinopyroxene, pigeonite and subsequently enstatite instead of pigeonite. At this point, the residual 1984 magma would have an MgO content around 4.2% and been about 50% solidified. Subsequent recharge and mixing by dominantly reservoir magmas, derived from a deeper 3–4-km intermediate magma reservoir, eventually produced the 2022 magma.

Hawaii↗
USGS2026

Loss and transformation of coastal wetlands due to global change in the conterminous United States: Past, present, and future

Coastal wetlands are being transformed by global change, impacting the ecological and societal benefits provided by these ecosystems. Synthesizing knowledge of historical and expected future transformations in coastal wetlands can help inform forward-looking planning and stewardship efforts. Here, we review anticipated future ecological transformations in coastal wetlands of the conterminous United States and contrast them with past transformations. We examine trends at the national scale but focus in detail on the following six wetland-rich regions: (1) Everglades, (2) Mississippi River delta, (3) Chesapeake Bay, (4) San Francisco Estuary, (5) Puget Sound, and (6) Great Lakes. This review is among the first to integrate and compare coastal wetland transformations in the Great Lakes with their oceanic counterparts. Between 1850 and 1970, wetland losses across the nation were high because coastal wetlands were viewed as flooded wastelands needing to be drained, filled, or excavated to make room for other land uses. Over the last 50 years, growing public and scientific appreciation of the ecological and societal benefits provided by coastal wetlands has fostered tremendous gains in wetland restoration and legal protection. Looking to the future, climate change and accelerated sea-level rise have become critical threats to coastal wetlands and are emerging as major drivers of ecological loss and transformation. This review synthesizes regional knowledge regarding past, present, and future changes to help coastal scientists, environmental managers, and the public better anticipate and prepare for future coastal wetland transformations due to climate change, accelerated sea-level rise, land-use change, and other aspects of global change.

conterminous United States↗
USGS2026

Improved estimates of relative occurrence and abundance using opportunistic surveys and presence-only observations: A zero-inflated integrated species distribution model

Modeling tools for estimating and forecasting shifts in species distributions are becoming increasingly valuable for conservation planning and response. This is especially true for wild bird populations, which have been declining across habitats and regions. Species distribution models (SDM) represent a diverse set of tools with options for addressing various sources of bias. Complex spatial processes associated with rare or clustered species can be accounted for using zero-inflated SDMs, whereas biased survey data can be integrated with additional data sources to improve estimates. Each option addresses an important and common source of bias, but the two SDM frameworks have not been implemented together previously. We present a novel zero-inflated extension of an integrated SDM framework for a Poisson regression (ZI-iSDM), allowing for the estimation of independent occurrence and abundance processes by integrating opportunistic survey and presence-only data. We validated performance of this ZI-iSDM using simulated datasets under different degrees of species rarity and density on the landscape as well as sampling bias in opportunistic datasets. We additionally applied this model to real survey data for multiple wild bird species using publicly available observation data combined with open-access environmental information to describe habitat associations. We found that integrating presence-only data, such as banding or harvest events, can compensate for potential deficiencies in opportunistic surveys by expanding sampling to be more representative of available and used habitat. Additionally, models that first differentiated occurrence and abundance using a ZI term were better suited for approximating distributions of spatially clustered species.

Ecological Modeling↗
USGS2026

A practical framework for GT-seq panel optimization

Genotyping-in-thousands by sequencing (GT-seq) panels are powerful tools in ecological, evolutionary and conservation genomics, yet the optimization process critical for robust and reproducible genotyping remains poorly formalized. Here, we present an iterative workflow for GT-seq panel optimization that emphasizes systematic refinement, quality control and structured decision-making to improve panel performance across diverse populations and study contexts. We illustrate this framework through the development and optimization of a GT-seq panel for white-tailed deer, a widely distributed and ecologically important North American species. From an initial set of 1200 candidate SNPs selected from a commercial microarray (OVSNP60, containing 72,728 SNPs) and prioritized for high heterozygosity, primers were designed for 646 loci. The final optimized panel contains 508 high-performing markers retained after iterative removal of overamplifying primer pairs, adjustment of primer concentrations, PCR conditions and bioinformatic filtering. The overall proportion of SNPs with more than 70% genotype rate increased from 25.5% in the first optimization round to 87.8% in the final round. Consequently, the overall genotype rate increased from 39.4% to 84%. We also identify key quality-control checkpoints and practical criteria to guide panel refinement and ensure consistent performance. By prioritizing optimization as an integral component of GT-seq panel development, this work provides a reproducible framework for generating robust, high-throughput genotyping tools in non-model species and underscores the importance of iterative refinement to maximize data quality and utility.

Molecular Ecology Resources↗
USGS2026

Bee-ing similar: Low diversity, no population structure, and signals of adaptation in two North Dakota bumble bees

Habitat fragmentation and destruction are driving widespread declines in ecosystem function, species abundance, and genetic diversity across the globe. Among the affected taxa are bumble bees (genus Bombus ), which are essential pollinators. Bumble bee declines have been linked to anthropogenic pressures such as land-use change and climate change. To better understand how the environment is shaping bumble bee populations, we employed a landscape genetics approach to examine the genetic diversity, population structure, and potential local adaptation of two widespread species— Bombus ternarius and Bombus griseocollis —across North Dakota. From 2017–2020, 161 B. ternarius and 200 B. griseocollis bumble bees were sampled in ND across 13 and 17 sites, respectively. We found low levels of heterozygosity and inbreeding across all populations for both species, with no evidence of population structure or isolation by distance. Our results revealed signatures of potential local adaptation to climatic variables and land cover characteristics, suggesting that our species are adapted to environmental gradients across the state. Having similar results between both species across all analyses suggests that other wide-ranging bumble bees in the region may share these patterns. High connectivity can buffer short-term population losses, but low genetic diversity may constrain adaptive potential and reduce resilience to future environmental change and emerging threats. Our findings have broad applicability to other pollinators and taxa facing similar environmental pressures. Both species have populations that continually exchange genes creating widespread connectivity but are still locally shaped by adaptation.

North Dakota↗