USGS Science⌕ Search

Geology topics

Browse USGS reports

Search indexed USGS publications on groundwater, aquifers, geologic maps, mineral resources and earthquakes. Explore source records by subject and place.

At least 991 records · Page 55Linked to original sources

Predicted genetic consequences of alternative population control strategies for North American plains bison in Yellowstone National Park

North American plains bison ( Bison bison bison ) in Yellowstone National Park are descendant from a population low of 23 indigenous individuals from the late 1800s, and 21 individuals introduced from outside herds in the early 1900s to strengthen the possibility of recovery. Within the park, the herd has rebounded and now averages around 5,000 animals. The herd is managed as a closed population within the park and adjacent areas in Montana, USA, which involves annual removals of animals to control numbers. We used genetic simulation to first isolate the effects of population size on retention of genetic diversity, and second to predict the genetic consequences of different population control strategies that vary in the number of bison removed annually, and their relatedness, while considering varying degrees of male dominance and environmental catastrophes. We found through simulation of a stable population that a herd size >2,500 had a high probability (>0.9) of retaining genetic diversity at levels compatible with long-term conservation (>90% of the initial genetic variation over 200 years), but a herd of 3,500 or more had a high likelihood (>0.9) of retaining levels of genetic variation (95%) suggested for maintaining evolutionary potential. When simulating the removals based on plausible management strategies available to the National Park Service, we found that maintaining a herd size of 3,500 or more, limiting removals to less than 40% of the population at a time, keeping a balanced sex ratio, and preferentially removing related individuals would likely retain more than 95% of existing genetic variation. With a herd size of 3,500 animals or more, population catastrophes and male dominance only influenced genetic retention in the most extreme cases, such as a catastrophic population decline of over 95%, a prolonged decline from 90% to 10% over 5 years, or extreme male dominance where 1% of males produced all offspring each year.

Idaho, Montana, Wyoming↗

Variable pH habitats could help prepare stone crabs for coastal acidification

Coastal acidification is being exacerbated by terrestrial organic inputs, especially after high precipitation events. Florida’s rainy season coincides with stone crab ( Menippe mercenaria) reproduction, and pH extremes could limit future harvests by reducing reproductive output. Populations that experience pH variability can serve as “natural laboratories” for estimating the tolerance of a species to coastal acidification. Here, we conducted a series of experiments to determine if ovigerous stone crabs conditioned in more variable pH habitats (seagrass) would result in faster embryonic development, greater hatching success, and higher larval survival relative to crabs conditioned in habitats with lower pH variability (sandy habitats). After field conditioning, crabs were transported to the laboratory and randomly acclimated to either a control pH (pH = ~ 7.90) or a reduced pH condition (pH = ~ 7.60) until larval release. The rate of embryo development was slower in the laboratory reduced pH treatment, however, there were no observable field condition effects on embryo development rate. Crabs conditioned in the more pH variable seagrass habitat did have greater hatching success and higher larval survival than crabs in less pH variable sandy habitats; however, larval survival was low across all treatments. These results suggest that the pH variability experienced in seagrass habitats during brooding may serve as a mechanism for stone crabs to acclimatize to extremes in seawater pH.

Marine Biology↗

Thematic accuracy assessment of the National Land Cover Database (NLCD2021) for the conterminous United States

The MultiResolution Land Characteristics (MRLC) consortium’s National Land Cover Database (NLCD) was established as an operational, stakeholder-oriented land cover monitoring program with the release of the NLCD2001 database. Both land cover and land cover change data became available with the release of NLCD2006. Here, we report the land cover and land cover change accuracies for NLCD2021, documenting Level II and I accuracies for the 2016, 2019, and 2021 land cover datasets and 2016–2021 land cover change. An additional objective, not addressed in previous NLCD assessments, is the estimation of the total variance of area estimates, where total variance includes sampling variance and variability in reference label assignment. We estimate the total variance of the area estimates for the stable and change classes targeted by the strata used in the sampling design. For the 16 Level II classes, the overall accuracy (OA) was 73% ± 1% when agreement was defined using only the primary reference class label and 89% ± 0.7% when agreement was defined as a match between the map label and either the primary or alternate reference label (± standard error, SE). For the eight Level I classes, the corresponding OA estimates were 81% ± 0.9% and 89% ± 0.7%. Level II and Level I OA tended to be higher for 2016 land cover and lower and more similar for 2019 and 2021 land cover. Land cover change user and producer accuracies (UA and PA) tended to be <50%. When the alternate reference label was included in the definition of agreement, exceptions of high class-specific accuracies (≥75%) were more prevalent for PA than UA. The consistency of reference label interpretations was greater for no-change classes than for change classes. The estimates of the total variance incorporating interpreter variability were often smaller than the standard variance estimates, indicating the possibility that the total variance estimator is highly unstable. Further study is needed to improve the utility of this total variance estimator for practical applications.

conterminous United States↗

Beyond pro-vs-anti hatchery: A typology of fishing community perceptions of salmonid hatcheries in Oregon

Objectives: Pacific salmon and steelhead hatcheries are central to fisheries management along the West Coast of North America, yet they remain a source of persistent controversy. While often framed in binary terms of pro- or anti-hatchery, public perspectives are more nuanced. Methods: We conducted 18 in-depth interviews with members of Oregon’s non-tribal fishing communities, including recreational anglers, professional fishers, non-governmental organizations, hatchery workers, and policymakers, to explore perceptions of salmonid hatcheries. Results: We developed a typology of perceptions of salmonid hatcheries in Oregon, identifying five distinct belief-based types, which we label Advocates, Reformers, Balancers, Skeptics, and Restorers. These types differed in how they evaluated hatcheries’ social and ecological impacts, governance legitimacy, and management acceptability. While Advocates and Reformers emphasized hatcheries’ social, economic, cultural, and ecological benefits, Skeptics and Restorers prioritized minimizing ecological harm and questioned hatchery effectiveness. Balancers acknowledged both social benefits and ecological risks, supporting context-specific management. All groups articulated different meanings of the term “wild,” suggesting that different perspectives across typologies may exacerbate fundamental miscommunications. Conclusion: Our findings highlight the complexity of stakeholder views and the limitations of binary framings. We suggest that tailored communication and engagement strategies, informed by typologies, may foster more inclusive and productive dialogue in hatchery policy and management.

Oregon↗

Laboratory evaluation of phototaxis in cultured Cisco larvae

Objective Field observations in a Lake Ontario embayment suggested that there is a phototaxic response of larval coregonines that changes with development. Methods We used a controlled laboratory experiment with hatchery-reared larvae exposed to three different light regimes (24-h light, 24-h dark, or 12-h light followed by 12-h dark) to test for phototaxis in Cisco Coregonus artedi . Results Larval depth distributions were not uniform, and contrasts between vertical distributions within early and late portions of the time series when larvae were exposed to light showed evidence of phototaxis and transition to a preference for the bottom around 68 d posthatch. Larval behavior in dark treatments contrasted with those exposed to light and also suggested that other influences on depth distribution could be operating. Conclusions As observed in other species, larval Cisco phototaxis response would likely result in diel vertical migrations and may be an adaptation to changes in food distributions balanced with predation risk as the larvae grow and become stronger swimmers. This has implications for the success of Cisco restoration in the Great Lakes, particularly as climate change affects environmental conditions and phenology.

Transactions of the American Fisheries Society↗

Annotated bibliography of scientific research on new world screwworm ( Cochliomyia hominivorax ) myiasis in wildlife

The New World screwworm ( Cochliomyia hominivorax ; NWS) is a parasitic blowfly that lays its eggs in open wounds of live, warm-blooded animals including livestock, wildlife, and potentially humans. The larvae consume living animal tissue, and if untreated, the infestation can lead to death. Although NWS was eradicated in the United States in 1966, it has been moving northward from its endemic range in South America during the past decade and could seriously threaten the health of U.S. wildlife populations, making detection, treatment, and surveillance of the disease far more difficult across this multi-sector disease system. As the likelihood of NWS reintroduction to the United States increases, veterinarians, wildlife managers, and conservation specialists need to be informed and prepared to respond. The existing knowledge about NWS interactions with wildlife hosts is lacking, especially regarding North American species where the NWS has been eradicated for more than 50 years. To address this knowledge gap, we compiled an annotated bibliography that consolidates key information from the existing literature on NWS infestation in wild animals.

Open-File Report↗

Assessment of groundwater quantity and quality contributions to Lake Huron

Lake Huron, one of the five Great Lakes, borders the United States and Canada, with Michigan as the only U.S. State on its shoreline. Like other freshwater lakes, it faces water-quality challenges from nutrients and chemicals applied across its drainage basin. Although past studies focused on surface-water sources, groundwater contributions remain less understood. To address this gap, the U.S. Geological Survey, as part of the Cooperative Science and Monitoring Initiative, classified drainage basins to Lake Huron into eight hydrogeologic zones based on bedrock rock type and glacial sediment transmissivity. Utilizing existing data and empirical field data, we quantified groundwater discharge and identified areas of concern for loading of chloride and nitrate to Lake Huron. Groundwater contributions, including indirect and shoreline discharge, ranged from 5.8 to 11.5 inches annually, totaling 1.9 cubic miles and 0.09 cubic mile, respectively. Hydrogeologic zones with higher glacial sediment transmissivity yielded greater indirect groundwater discharge. Chloride levels above the U.S. Environmental Protection Agency’s 250-mg/L recommendation were mainly in the Saginaw lowlands, whereas nitrate above the 10-mg/L standard was rare—found in only 11 wells. Together, the analysis of where groundwater discharge is occurring in the Lake Huron Basin and the identification of areas with potential groundwater-quality concerns can help prioritize areas that are critical to protecting the long-term health of Lake Huron.

Michigan↗

Line transect distance sampling and genetic analyses reveal a small but genetically diverse coastal Gopher Tortoise (Gopherus polyphemus) population

Gopher Tortoises inhabit coastal systems, including barrier islands, across the southeastern U.S. St. Vincent National Wildlife Refuge is an uninhabited barrier island located off the coast of northwestern Florida. Although tortoises have been observed on the island, no information is available on the status of the population. We conducted a line transect distance sampling survey to evaluate the Gopher Tortoise population on St. Vincent Island. Additionally, we collected samples for genetic analyses from 11 individual tortoises captured opportunistically and via bucket traps on the island and 15 tortoises captured at nearby mainland sites. Surveys covered approximately 43% of the sampling frame and resulted in 55 burrows, 28 of which were occupied. The abundance estimate for the island was 52 tortoises (95% confidence interval [CI] = 27–100) and the density was relatively low at 0.071 tortoises/ha (95% CI = 0.037–0.136). Genetic analyses of two mtDNA markers identified a new haplotype unique to St. Vincent Island and another three haplotypes previously found across the southeastern U.S. The genetic composition of Gopher Tortoises on St. Vincent Island is representative of the entire southeastern U.S. but most closely aligns with tortoise populations east of the Apalachicola-Chattahoochee River system. Although the tortoise population on this island is small, the extent of seemingly appropriate habitat on the island and the genetic diversity of the population suggests the potential for growth with added management intervention.

Florida↗

Quantitative mineral resource assessment of lithium pegmatite deposits in the northern Appalachian orogen, USA

Lithium demand is projected to increase more than 48 times by 2040 due to electric vehicle production and other energy storage needs. Most lithium production is outside of the USA, thereby increasing supply chain vulnerability. The combined end use importance and heightened supply risk of lithium make this lightest metallic element a critical commodity to the USA. To mitigate this supply risk, the US Geological Survey is actively assessing lithium deposits in the USA. Herein, we detail an assessment for lithium-mineralized pegmatites in the US northern Appalachian Mountains. Permissive tracts were generated by cross-referencing tectonic and geologic maps and mineral occurrence data with mappable criteria derived from generalized and region-specific lithium pegmatite ore deposit models; tracts were then ranked as having high, medium, or low permissibility. Available geophysical and geochemical data were found to be of minimal utility for this deposit type at the scale of the assessment. The number of undiscovered deposits were estimated and integrated into probabilistic simulations, which included an expanded and updated global grade and tonnage model of pegmatite-hosted lithium ore. The estimated total amount of undiscovered resources for the northern Appalachian Orogen has a median value of 1,410,000 metric tons of Li 2 O when considering moderate correlation across sub-regions. At a confidence level of 90%, a resource of at least 90,000 metric tons of Li 2 O remains undiscovered, and at a 10% confidence level, a resource of as much as 7,380,000 metric tons Li 2 O remains undiscovered. After applying an up-to-date economic filter to convert median contained lithium to recoverable material, a correlated total of 900,000 metric tons of Li 2 O may be economically extractable, equating to enough Li 2 O to provide the current annual US lithium supply deficit (presently obtained through net imports) for 127 years at 2025 rates of apparent consumption. This period of provision will inevitably shorten with projected increasing consumption rates, emphasizing that further research could be completed to better delineate regions of high lithium resource potential and support exploration and domestic production.

Connecticut, Delaware, Maine, Massachusetts, New H↗

The role of groundwater in contributing to surface water salinization in the Upper Colorado River Basin

Freshwater salinization impacts the availability of water for human use and ecosystem needs worldwide. It has been estimated that total dissolved solids (TDS) in the Colorado River Basin cause $350 million/year in damages and substantial resources are devoted to reducing TDS loading to streams. This study describes the development and application of coupled watershed models that enable TDS source tracking through the subsurface and across the landscape at a seasonal timestep for 35 years in the Upper Colorado River Basin. Results indicate that, on average, 75% of TDS loading to streams originates as baseflow, and 50% of loading is lagged in delivery by longer than one season. Snowmelt was identified as a dominant process controlling the transport of lagged TDS to streams. This approach informs when and where TDS mitigation efforts may be effective in a watershed that serves as a critical water supply for the southwestern United States.

Arizona, Colorado, New Mexico, Utah, Wyoming↗

A Bayesian hierarchical modeling approach for species diversity in ecology

Species diversity is the foundation of many ecological disciplines. This metric is often approximated using species richness and evenness, even though actual richness likely exceeds observations due to imperfect sampling methods. Estimating the “true” species richness, which includes identifying the number of missing species, has intrigued ecologists for decades. We adopted a parametric model that appeared in Fisher et al. (1943), which models the numbers of individuals from different species as random samples from a negative binomial distribution, and developed a Bayesian computational approach to directly estimate the distribution model parameters. The model parameters represent species abundance and evenness, and can be used to derive species richness. We evaluated our parametric approach using (1) a simulation study and (2) three historical data sets. Furthermore, we illustrated the hierarchical modeling approach to combine data from multiple parallel studies using a biannual fishery survey data set. Our parametric model formulation is computationally efficient, and the hierarchical structure facilitates embedding diversity estimation into broader application, such as assessing spatial and temporal trends in species diversity associated with environmental stressors. Additionally, because the two parameters of the negative binomial distribution model represent species abundance and evenness of a community, this parametric approach facilitates a deeper understanding of the ecological systems under study. The negative binomial distribution model works with a wide range of species frequency distribution types. As a result, our emphasis on a parametric model can help us characterize the structure of an ecosystem and provide a greater depth of ecologically meaningful information.

Ecological Informatics↗

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

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

JGR Solid Earth↗

A practical decision tool for marine bird mortality assessments

Given the rise in anthropogenic, environmental, and disease events contributing to marine bird mortality, there is a critical need to improve the rigor of mortality assessments. Deficits in data collection and mortality estimation can hinder a manager’s ability to document the scale of events and assess population level impacts. Therefore, to inform decisions required during activities, such as conservation status assessments or harvest management, organizations may choose to incorporate mortality assessments into response plans. Resources, capacity, and assets to assess mortality vary across jurisdictions (federal, state, Indigenous, local, etc.), and clear guidance to support mortality assessments is often unavailable or not clearly addressed. Here, we present a decision support tool to help managers identify and evaluate survey options to assess bird mortality in a diverse array of scenarios. The objective of the decision tool is to improve data collection and availability, which will increase the ability to estimate mortality robustly, given situation-specific attributes and constraints. This decision tool is designed to guide the response when a mortality event is initially encountered and offers suggestions for assessment and reporting procedures in the absence of other guidance or to complement existing protocols. The decision tool is also meant to inform decision making for response determination and resource allocation. The tool facilitates examination of options for further assessment and monitoring, which users determine by examining questions pertaining to species prioritization, determination of mortality minimum spatial extent, and the potential magnitude of impacts on affected species. Finally, identification of appropriate survey methods that address imperfect detection when a complete census is not possible are determined by exploring location, spatial and temporal extent, and the type of species affected. Ultimately, this decision tool aims to facilitate and improve the standardization of mortality assessments, equipping managers with a practical resource to navigate the decision-making process for marine bird mortality estimation.

Ornithological Applications↗

Hydrogeology, groundwater salinity distributions, and assessment of the effect of oil-production activities on groundwater in the Midway Valley area, western Kern County, San Joaquin Valley, California

This study seeks to determine the effects of oil field produced water disposal operations and well mechanical integrity issues on groundwater quality in oil fields in the southwest San Joaquin Valley, California. Whereas previous studies used groundwater wells to study shallow aquifers outside the oil fields, this study demonstrates that future approaches may use oil well geophysical logs to map groundwater head gradients, create salinity profiles and document changes in salinity over time in oil field areas with sparse groundwater well data and at depths greater than 330 m. We also incorporate an analysis of well histories to determine potential effects of compromised wellbore seals on changes in aquifer quality that cannot be explained by water disposal practices. Water quality in the aquifers is naturally brackish across most of the area, with better quality groundwater occurring in the eastern part. Geophysical logs are used to determine salinity variations within aquifers including the depth at which TDS exceeds 10,000 mg/L. This depth ranges from 366 m in the northwest to approximately 1,500 m in the southeast. Oil well porosity logs are used to determine water table elevations. These logs indicate the water table slopes south-southeast, showing the predominant groundwater flow direction is from oil field disposal areas toward better quality groundwater east of the oil fields. Geophysical logs show formation resistivity near some disposal facilities has decreased over time, indicating the salinity of the aquifer has increased due to disposal of saline produced water in injection wells and ponds. Oil well history analysis suggests that increased salinity over time in water-saturated sand intervals >1.5 km from disposal facilities may be caused by mechanical failures and/or incomplete borehole seals in poorly constructed or abandoned wellbores prevalent throughout the study area—particularly wells drilled prior to 1930.

California↗

ShakeAlert Earthquake Early Warning System performance during the Mw 7.0 offshore Cape Mendocino earthquake

The 5 December 2024 M w 7.0 Offshore Cape Mendocino earthquake was a challenging test of the U.S. West Coast ShakeAlert earthquake early warning system due to its offshore epicenter and limited near‐source station coverage. We analyzed real‐time performance of all components of the ShakeAlert system, including the seismic algorithms (earthquake point‐source integrated code [EPIC] and Finite‐fault rupture Detector [FinDer]), the geodetic algorithm (Geodetic First Approximation of Size and Time–peak ground displacement [GFAST‐PGD]), and network telemetry during the event. EPIC created the first solution for this earthquake 15 s after origin time with an initial magnitude estimate of M 5.6 and location error of 10 km from the Advanced National Seismic System epicenter. An early spurious trigger from station CE.89101 fortuitously maintained location accuracy and, correspondingly, magnitude accuracy. FinDer contributed its first solution at 18 s with a location estimate closer to the seismic network and produced two distinct rupture geometries, leading to minor fluctuations in estimated intensity contours. GFAST‐PGD did not meet alerting thresholds but otherwise performed as expected. Network latencies were <2 s for most stations, supporting the rapid detection of this earthquake by the system. Roughly five million alerts were delivered to cell phone devices in California and Oregon during this event. This was also the first instance of a school district‐wide ShakeAlert‐powered system being activated. Comparisons to recorded seismograms demonstrate that the maximum warning times before potentially damaging shaking (intensity 6+) were in the range of 5–55 s. Although the ShakeAlert system provided accurate solutions and useful alert delivery, this earthquake raised awareness of potential issues within the system, including the need for improved offshore location estimates, a combination of solutions from ShakeAlert servers, and handling of spurious triggers.

California↗

Stratigraphic revision of the Neogene rocks of Saipan Island, Commonwealth of the Northern Mariana Islands

Geologic mapping in 2006 and 2007, and calcareous nannofossil biostratigraphy combined with comparative analysis of other published data, have resulted in a reinterpretation and revision of the ages and stratigraphic relationships of the Neogene lithostratigraphic units exposed on Saipan Island. These units, chiefly limestone and admixed volcaniclastic rocks, form much of the bedrock of Saipan Island. Stratigraphic revisions include a redefinition of the Tagpochau Limestone, restricting it to rocks of the newly named I Etdot, Tipo Poli, and Marpi Members. Other rocks, formerly included in the Tagpochau Limestone (the Donni Sandstone and the Machegit Conglomerate Members, and an informal transitional facies), are herein assigned to the newly defined Gloria Formation. The Gloria Formation comprises three members: the Machegit Conglomerate Member, the new As Teo Member, and the Donni Sandstone Member. The Tagpochau Limestone contains larger foraminifera that indicate an early Miocene (Aquitanian through middle Burdigalian) age, or a numerical range of about 23 to 18.7 Ma. Overlying the Tagpochau Limestone, both the Fina-sisu and the Gloria Formations are recognized herein to be distinctly younger than the Tagpochau Limestone and not correlative with any parts of that formation as had been previously thought. Basaltic lava flows in the Fina-sisu have yielded radiometric ages that indicate a middle Miocene age, whereas planktonic foraminifera recovered from interbedded sedimentary rock indicate a latest early Miocene age (Langhian). Assemblages of planktonic foraminifera and calcareous nannofossils recovered from the Gloria Formation indicate a latest Miocene age (Messinian) and a numerical age range of about 5.8 to 5.5 Ma. Additional studies could help to constrain the ages of these units more precisely because none of the biostratigraphic samples were taken in the context of measured stratigraphic sections. On an island dependent on groundwater, the revised stratigraphic understanding can inform the development of groundwater models and strategies for identifying water-well sites and targeting drilling depths for potable water. The new stratigraphic understanding can inform future interpretations of the geologic structure of Saipan Island and the Neogene tectonic history of Saipan Island and the Mariana subduction zone.

Professional Paper↗

Evaluating snow depth measurements from ground-penetrating radar and airborne lidar in boreal forest and tundra environments during the NASA SnowEx 2023 campaign

Snow is a vital component of high-latitude terrestrial systems, but environmental factors (e.g., permafrost) and complex vegetation challenge the accurate measurement of key snowpack properties. We evaluated local-scale ground-penetrating radar (GPR) and large-scale airborne lidar retrievals of snow depth collected during the NASA SnowEx 2023 campaign in tundra and boreal forest environments in Alaska along 44 short (3–12 m) transects. Compared to in situ observations, we identified modest biases for GPR snow depths (bias <0.03 m in tundra, +0.06 m in boreal forests) and larger biases for lidar snow depths in the boreal forests (–0.16 m). At the Upper Kuparuk-Toolik tundra site, lidar snow depths exhibited a small bias (–0.02 m), whereas the bias was much larger at the Arctic Coastal Plain tundra site (+0.19 m). For most sites, biases were primarily related to sub-snow vegetation, tussocks, and seasonally dynamic ground. However, we identified vertical alignment issues with the Arctic Coastal Plain lidar snow depth dataset that likely contributed to the higher bias. The complex ground surface and sub-snow vegetation in these environments present a challenge to established snow depth measurement methods, which needs to be considered when evaluating novel remote sensing approaches.

Alaska↗