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.
Using 3D Elevation Program data, a modern illustration of each State's landscape showcases its elevation differences in vivid detail. Each State has been consistently adjusted and colorized according to its highest and lowest elevations, accentuating differences even in States with low relief. To further emphasize landforms, shaded relief has been added on top of the elevation colors.
Using 3D Elevation Program data, a modern illustration of each State's landscape showcases its elevation differences in vivid detail. Each State has been consistently adjusted and colorized according to its highest and lowest elevations, accentuating differences even in States with low relief. To further emphasize landforms, shaded relief has been added on top of the elevation colors.
Using 3D Elevation Program data, a modern illustration of each State's landscape showcases its elevation differences in vivid detail. Each State has been consistently adjusted and colorized according to its highest and lowest elevations, accentuating differences even in States with low relief. To further emphasize landforms, shaded relief has been added on top of the elevation colors.
The giant gartersnake ( Thamnophis gigas ) is a semi-aquatic species of snake endemic to marshes and sloughs in the Central Valley of California. Because the Central Valley has lost a substantial amount of wetland habitat, giant gartersnake populations are typically found in aquatic habitats associated with rice ( Oryza sativa ) agriculture. Therefore, changes to the amount and location of rice agriculture on the landscape could affect giant gartersnake populations. We present analyses of 9 years of data examining how the occurrence of giant gartersnakes is affected by rice idling. We sampled agricultural canals associated with active rice fields and fallowed fields at 119 sites from 2016 through 2024. We captured giant gartersnakes at 95 of 119 sites sampled during our study, totaling 5,269 captures of 2,687 individuals. We used these capture data, as well as data on the trapping effort, habitat, vegetation, prey availability, and amount of active rice growing at each site, to estimate occupancy and detection probabilities of giant gartersnakes. Mean daily detection probability at a site was 0.38 (95-percent credible interval 0.22–0.56). Given the mean daily detection probability, the probability of detecting giant gartersnakes at an occupied site at least once during a 21-day sampling period within a given year was greater than (>) 0.99 (>0.99 to much greater than [>>] 0.99). The spatial extent of active rice growing within a 2-kilometer buffer of a site was positively correlated with initial site occupancy and the probability of site colonization. Giant gartersnake population density (hereinafter, density) estimated using spatially explicit capture-recapture models, varied among traplines and years, with densities consistently highest in the Butte basin. The annual population growth rate also varied among site-clusters and years, with no temporal trend of growth or decline among sites. Annual apparent survival of giant gartersnakes varied slightly between years and among site clusters, with overall mean annual apparent survival of 0.51 (95-percent credible interval 0.34–0.71) for female snakes and 0.45 (0.30–0.63) for male snakes. The amount of rice in the surrounding landscape was positively correlated to giant gartersnake survival. Rice contributes positively to giant gartersnake occurrence, colonization, and survival, and maximizing the amount of rice agriculture within 4 kilometers of canals with giant gartersnake populations would likely minimize the effects of rice idling on these populations.
The U.S. Geological Survey, in cooperation with the Gulf Coast Ecosystem Restoration Council, has used MODFLOW 6 to develop a groundwater-flow model to simulate groundwater and surface-water interactions in the Mobile Bay, Alabama, area. The model results indicated that, on average, groundwater discharge near the coastline is equal to 2.5 percent of the surface water that flows into the bay. The model was also used to determine how changes in recharge, sea level, and groundwater pumping affect groundwater levels and discharge rates. The results indicate that more groundwater discharge occurred in the winter and spring when recharge was higher, sea level was lower, and groundwater pumping was lower than during the summer and fall. Additionally, the amount of emergent groundwater was closely related to sea level; when sea level was higher, there was more area with emergent groundwater. Furthermore, the depth of nonemergent groundwater was related to trends in recharge and pumping. During periods of increased recharge and reduced pumping, a significant portion of the model area exhibited a depth to the water table of less than 1 meter below the land surface. Conversely, during periods of decreased recharge and heightened pumping, much of the model area showed a depth to the water table ranging from 1 to 5 meters.
Detailed and accurate historical data are crucial to characterizing landscape change, whether it is caused by natural phenomena or of anthropogenic origin. Aerial imagery has been periodically collected by government and nongovernment agencies, for a broad range of purposes, since the early 20th century. The acquisition of aerial images involves using standardized equipment with large-format, carefully calibrated cameras and lens assemblies. Historical archives of aerial images preserve a record of the landscape and are invaluable sources of information that predate satellite and lidar. Structure from motion (SfM) aided photogrammetry provides an efficient means for creating digital three-dimensional data from the vast archives of historical aerial imagery. This report presents techniques, methods, best practices, and a standard workflow for creating point clouds, digital elevation models, and orthomosaic digital scans of overlapping (stereographic) vertical aerial photographs using photogrammetry software. A list of critical issues relevant to SfM-aided photogrammetry of historical aerial imagery is provided, along with potential contingencies and solutions. The workflow was tested in many geographic study areas and with historical imagery datasets of different scales and sources. Three case studies are included. By integrating known SfM methods, these procedures enable the consistent and optimal extraction of historical quantitative data from archives like the U.S. Geological Survey (USGS) Aerial Photo Single Frame archive in the USGS EarthExplorer data portal.
Subsurface water ice deposits on Mars are an important resource for potential future human exploration. They are also an indicator of the planet’s past climate. However, the distribution of subsurface water ice in Mars’s midlatitudes is uncertain because spacecraft imagery cannot directly observe subsurface ice in most cases. Various spacecraft remote sensing instruments are sensitive to subsurface water ice, including thermal imaging spectrometers, radar sounders, and neutron spectrometers. Geomorphic analyses of images can also implicate subsurface ice. Building upon the data products from the Mars Subsurface Water Ice Mapping project, we provide a probabilistic framework to jointly interpret existing data and estimate the likelihood of subsurface water ice in the Martian midlatitudes between 60 ∘ S and 60 ∘ N with uncertainty. Broadly, we find that near-surface ice is likely present poleward of ∼45 ∘ in both the northern and southern hemispheres. However, closer to the equator, existing remote sensing data cannot uniquely constrain the presence of subsurface water ice. Our probabilistic results provide a framework for quantifying the abundance of ice on Mars, and our uncertainty estimates allow future analysis and exploration to target regions of high uncertainty.
Biological nitrogen (N) fixation (BNF) provides the N needed to produce proteins and other biological building blocks, helping feed humanity and mitigate climate change. Due to its high energetic cost compared to other forms of N acquisition, biotic investment in BNF indicates N limitation. Globally gridded BNF flux data provide an opportunity to determine the energetic investment in BNF across ecosystems, which would help reconcile conflicting indicators of N limitation. Here, we use a new BNF synthesis to quantify the relative importance of BNF in different N-fixing niches, in different biomes, and across the globe by calculating the fraction of net primary productivity (NPP) invested in BNF and the fraction of plant N acquisition provided by BNF. Larger fractions of non-agricultural NPP were invested in BNF in less-productive, higher-latitude biomes. This pattern was driven by biocrusts and mosses. Similarly, non-agricultural symbiotic N-fixing plants invested relatively more of their own NPP in BNF in less-productive, higher-latitude biomes. This symbiotic plant pattern was driven by shrubs and herbs, overriding the opposite pattern in trees. Symbiotic plants also acquired a higher fraction of their N from BNF at higher latitudes and in less productive biomes. Investments in symbiotic BNF were 10× higher in agricultural (2.9% of NPP) versus natural (0.29%) biomes, providing 26% versus 3.1% of ecosystem-scale plant N acquisition. These results support the paradigm of strong N limitation at higher latitudes, help understand the rarity of N-fixing trees at higher latitudes, underscore the dominance of human activity, and inform terrestrial biosphere models.
Thawing of permafrost in northern latitudes is accelerating, potentially releasing substantial amounts of methane (CH 4 ) as forested permafrost plateaus transition into wetlands. This ecosystem shift alters the carbon exchange between the soil and atmosphere, influencing the permafrost-carbon feedback. Monitoring these changes may require measurement platforms operating across varied spatial and temporal scales. Recent advancements in small uncrewed aircraft systems (sUAS) enable high resolution CH 4 flux quantification in remote, complex terrains; however, comparisons with established methods such as eddy covariance flux towers remain limited. We used a hexacopter sUAS to quantify CH 4 emissions from the Alaska Peatland Experiment, a wetland within the Bonanza Creek Experimental Forest. Using an ensemble of methods to define the background CH 4 concentration, along with near surface emissions from soil chambers, helped constrain our flux estimates. The sUAS method yielded an average flux of 0.0077 ± 0.0019 mol s −1 CH 4 , within a factor of two concurrent tower-derived total source flux estimates (0.0036 ± 0.00042 mol s −1 CH 4 ). To assess spatial drivers of observed fluxes, we conducted a 2D footprint analysis and overlaid the results with high-resolution hyperspectral land cover classification, quantifying vegetative contributions within each footprint. This revealed higher fen representation in sUAS measurements (73.8%) than in tower footprints (58.8%), and lower tussock meadow representation (15.6% and 30.3%, respectively). These differences were consistent with known variation in vegetation-specific CH 4 emissions. Our results highlight that combining footprint modeling with land cover characterization can enhance interpretations of CH 4 fluxes and guide cross-platform comparisons.
Management decisions for novel wildlife threats are challenging to optimize due to deep ecological uncertainty, and multiple, often competing objectives. We present a transparent, transferable, multi-criteria decision framework that integrates quantitative forecasts of management outcomes with decision-maker value weights to identify robust strategies. Specifically, we evaluated options for responding to the imminent threat of Batrachochytrium salamandrivorans (Bsal) to North American amphibians. We worked with managers across the northeastern United States to characterize the full complexity of their multi-criteria decision contexts for addressing Bsal risk. We provide quantitative evidence that, across a range of value weights and landscape contexts, acting ahead of invasion is an optimal strategy for mitigating disease risk. We find clear evidence that proactive management is optimal, but our results also show that local context affects the optimal management action. Though demonstrated with Bsal, our framework generalizes to other biosecurity threats, providing decision support that unites ecological prediction, structured decision making, and stakeholder values to meet conservation goals.
Dispersal is a critical ecological and evolutionary process, yet the drivers of individual variation in emigration, particularly under human-modified conditions, remain relatively understudied. We investigate the determinants and patterns of emigration in a managed population of cooperative breeders using genetic samples from harvested individuals. We used a novel genetic approach to identify natal emigration through sibship analyses of DNA from harvested grey wolves ( Canis lupus ), a cooperative-breeding species, and evaluated how prey biomass, wolf density, human-caused mortality and sex influenced natal emigration across Idaho, USA, from 2017 to 2023. From our analysis of 182 sibling groups (representing packs), we evaluated natal emigration status for 224 individuals. Of these, 86 were classified as emigrants and 138 as non-dispersers. We assessed how covariates influenced the probability of emigration and evaluated sex differences in apparent dispersal distance. We found no sex differences in emigration propensity or apparent dispersal distance. Emigrants dispersed an average of 81.8 km (range: 30.0–338.9 km) from their inferred natal territories. Emigration was more likely from territories with higher prey biomass and lower wolf density, while human-caused mortality had no detectable effect. We found that ecological drivers, such as prey availability and population density, have a stronger effect on the propensity to emigrate than harvest intensity. These results highlight the importance of resource availability and population dynamics in shaping dispersal patterns in a cooperative breeder. Our study used a novel genetic method to assess natal emigration in a group-living carnivore and developed a framework broadly applicable to social systems where kin structures can be traced to dispersers. Our findings identify landscape-scale conditions that influence dispersal and can inform the management of wild populations.
This study provides a comprehensive assessment of decadal changes in the quality of groundwater used for public drinking-water supply at 444 monitoring sites across California during 2004–23. We assessed decadal step trends in groundwater quality for 145 water-quality constituents and geochemical indicators statewide and across geographic and land-use based network groups. We evaluated the statistical significance of directional changes (predominant increase or decrease of constituent concentrations) and the magnitude of those changes across all network groups. Uranium showed the most widespread directional and high-magnitude increases of all constituents with regulatory benchmarks statewide, particularly in the agriculture-dominated Central Valley as well as urban-and desert-dominated regions of Southern California. Fluoride and perchlorate showed the most widespread directional and high-magnitude decreases of all constituents with regulatory benchmarks statewide, which were also most pronounced in Southern California. Although arsenic and nitrate did not often register significant directional changes across network groups, they showed widespread, high-magnitude changes in both directions (increase and decrease) at levels often exceeding 10 percent of respective regulatory benchmarks statewide. Triazine herbicides (atrazine and simazine) and the gasoline oxygenate methyl tert -butyl ether showed significant directional decreases statewide, but not at levels considered to be of high magnitude compared to respective regulatory benchmarks. We observed significant directional and high-magnitude increases of total dissolved solids statewide, which were most pronounced in agricultural areas. Analysis of explanatory geochemical indicators indicated that prevalent statewide increases of alkalinity and calcium were the predominant components of the observed statewide increases in total dissolved solids by mass. Widespread increases in groundwater alkalinity and calcium across agricultural and urban areas may be related, in part, to warm-season irrigation and other anthropogenic factors that have shifted soil weathering dynamics over the long term. Increasing alkalinity concentrations were related to increasing uranium concentrations, particularly in areas with aquifer materials derived from granitic rocks. Conversely, increasing calcium concentrations were related to decreasing fluoride concentrations, particularly in areas where fluoride occurred naturally at elevated concentrations. Decrease of perchlorate, triazine herbicides, and methyl tert -butyl ether are likely related to decreased anthropogenic source inputs over time and natural attenuation in aquifers.
Introduction: Nonnative fishes in streams and inland lakes can substantially alter aquatic ecosystems by degrading habitat conditions, competing with native species, and restructuring food webs. Identifying invasion hotspots can inform proactive management and prevention of nonnative species establishment. Methods: In this study, we identified streams and lakes across the conterminous United States that are vulnerable to invasion by 18 nonnative stream and 28 nonnative lake fish species. Habitat suitability models were developed separately for stream and lake species across nine ecoregions using natural and anthropogenic predictors, and species occurrences were predicted by integrating model results with climate match scores based on comparisons between environmental conditions in species’ native and potential nonnative ranges. We classified areas with both high habitat suitability and strong climate similarity as having high invasion potential and generated hotspot maps by stacking species-level predictions. Results: Florida exhibits the highest invasion risk for both habitat types. Stream hotspots were concentrated near major metropolitan areas, whereas lake hotspots were more widely distributed across the Great Lakes basin and the northeastern U.S. Twelve species were predicted to invade both habitat types, with goldfish ( Carassius auratus ) and pirapitinga ( Piaractus brachypomus ) showing particularly high invasion potential. Discussion: These findings provide a comprehensive assessment of nonnative freshwater fish invasion risk and support targeted monitoring and management strategies.
Northern bobwhite ( Colinus virginianus ; bobwhite) populations have experienced an 85% decline across most of its range since the 1960s. The most drastic decreases have been in the southeastern United States where biologists attribute the decline to widespread habitat loss, including conversion of native grasslands to tall fescue ( Lolium arundinaceum ) and other exotic cool-season grass (CSG) pasture and hayfields. Including agricultural lands in conservation programs could improve habitat conditions on a regional scale. One working-lands conservation strategy involves the use of native warm-season grasses (NWSGs) rather than cool-season exotic grasses as economically viable cattle forage. To evaluate this management practice for creating bobwhite habitat, we conducted a field study on a bobwhite population in eastern Kentucky, USA, where grazed NWSG and burned NWSG fields were adjacent to exotic grazed and exotic hayed CSG fields. Between April 2019 and October 2022, we used radio telemetry to track bobwhites to evaluate resource selection at second- and third-order scales. Bobwhites used traditionally managed CSGs less than would be expected by chance at the second-order scale in both breeding and non-breeding seasons. Bobwhite use of grazed NWSG and burned NWSG was greater than would be expected by chance at the second-order scale in both breeding and non-breeding seasons. Similarly, at the third-order scale, bobwhites used CSGs less than would be expected by chance, whereas NWSGs and woody stems were used more than by chance alone. Our study suggests that bobwhites select NWSGs at multiple spatial scales despite broad-scale availability of managed exotic CSGs. Under a working-lands framework, the integration of NWSGs into working lands could create bobwhite habitat in the southeastern United States, especially if woody stems are present.
Artificial light at night (ALAN) poses a threat to ecosystems globally. It includes both direct and indirect light, or skyglow, which occurs when ALAN scatters in the atmosphere, extending beyond its original source. We analyzed ALAN trends in Lake Washington, WA, from 2014 to 2023 using Visible Infrared Imaging Radiometer Suite (VIIRS) nighttime light measurements, evaluated the relationship between in situ and satellite measurements, and modeled predator search volumes for a juvenile salmon predator, cutthroat trout ( Oncorhynchus clarki ), as a function of ambient light and turbidity conditions. Open water regions experienced significant increases in ALAN, while nearshore areas primarily showed no or negative trends, revealing the role of skyglow in shaping open water light environments. Using a visual foraging model for cutthroat trout, we found that juvenile salmon at the shallow southern pelagic site experienced light 28 times brighter, resulting in a 168% greater predation vulnerability than those at the northern site. In the Ship Canal, a narrow corridor for outmigrating salmon, predator search volumes were 249% higher than at the southern site. These contrasts in predation vulnerability demonstrate how local conditions influence predator–prey dynamics and provide critical insight for targeting mitigation of both nearshore and distant light sources.
Autonomous Recording Units (ARUs) are increasingly used in research to support passive acoustic monitoring, but deployment in remote or inaccessible locations remains logistically challenging. Traditional manual placement is labor-intensive, time-consuming, potentially hazardous, and often limited to habitat edges, creating sampling biases and restricting spatial coverage. To address placement limitations in wetland systems, we developed and field-tested a low-cost, lightweight floating platform paired with a drone-based deployment and retrieval system. The drone-deployable ARU platform, constructed from inexpensive, off-the-shelf materials (~US $21 per unit), weighed ~560 g and was mounted with an ARU (AudioMoth) at ~1.2 m above water. The platform was designed for stability, portability, and compliance with Federal Aviation Administration Part 107 regulations. Field trials were conducted across 50 operations in various types of historical rice-field impoundments along the South Carolina coast, which provide critical habitat for secretive marsh birds such as rails (e.g ., Rallus, Laterallus , and Porzana ) and bitterns (e.g ., Botaurus and Ixobrychus ). All 50 deployments and retrievals were successful, demonstrating the robustness of the method under diverse hydrologic and vegetative conditions. Deployment times were significantly shorter than retrieval times (median 6 vs. 9 min), with retrieval requiring greater precision for hook engagement. Operation times scaled predictably with distance: deployment increased by 0.93 min/100 m (R 2 = 0.94) and retrieval by 1.17 min/100 m (R 2 = 0.95). No platform damage or displacement was observed, as emergent vegetation was likely providing natural anchoring. Our approach offers an affordable and effective solution for expanding ARU coverage in large, inaccessible wetlands, reducing sampling bias and enhancing biodiversity monitoring. With continued advances in drone payload capacity, battery endurance, and beyond-visual-line-of-sight flight regulations, our workflow could be adapted for deploying multi-sensor platforms (e.g., ARUs, eDNA samplers, water and air quality sensors, plant sample collectors, insect traps, and trail cameras) to support integrated biodiversity monitoring in challenging landscapes worldwide.
Lead exposure has been documented in a multitude of bird and mammal species but less frequently in reptiles. Of the studies that have evaluated lead concentrations in reptiles, few have focused on snakes. We analyzed lead concentrations in liver and two types of bone from 24 gophersnakes ( Pituophis catenifer; Pc) and 12 Great Basin rattlesnakes ( Crotalus oreganus lutosus; Col) found dead along roads in the Morley Nelson Snake River Birds of Prey National Conservation Area in southwestern Idaho, U.S.A from 2018 to 2021. Because this area has been heavily impacted by both military and recreational activity, we also quantified lead concentrations in the soil near where the dead snakes were found. All soil samples had detectable lead (median: 12.53 µg/g; range: 0.07–16.77). Within individuals of both species, dry weight lead concentrations in liver generally were lower than those in bone (n samples < LOQ Liver: 9 of 13 ( Pc ), 9 of 10 ( Col ) (only a subset of individuals were in good enough condition to provide liver samples); Bone: 0 of 24 ( Pc ), 3 of 12 ( Col ). We did not detect any differences in lead concentrations between the species for either tissue type, nor between bone from the front (anterior) or back (posterior) of the snake skeleton. Additionally, we did not detect any correlation between lead concentration in soil collected at the roadkill sites and the tissue in the snakes. These data provide baseline information for lead exposure of two snake species in southwestern Idaho and demonstrate how lead is present in higher trophic levels within this ecosystem.
Subduction zones are important drivers of Earth's geochemical evolution; the presence of water within the subduction system has substantial implications for geohazards, including the style of seismicity and the volume of magma generated at arc volcanoes. We utilize estimates of pore fluid distributions in the outer-rise region of the Central American margin imaged using controlled-source electro-magnetic surveys and novel two-dimensional forward models simulating poro-elastic deformation of slab bending to constrain the background permeability structure of the Cocos Plate and within outer-rise fault damage zones. The model results suggest that a relatively low background permeability structure is required to match the observed porosity at depths greater than 3 km (less than ~10 −21 m 2 . Similarly, increasing the permeability within the damage zones of outer-rise faults by more than a factor of 10 also over-predicts observed porosity distributions at depths greater than 3 km. However, the simulations consistently underpredict pore fluid distributions between depths of 1–3 km, which may reflect nonlinear elastic or anelastic deformation processes not included in the simulations or inherent uncertainty in the conversion of resistivity measurements to porosity distributions. The best-fitting depth-permeability relationship is relatively low, which has important implications for the timing of pore fluid release down-dip of the trench, as it restricts shallow fluid release, likely leading to higher pore fluid pressures, more widespread hydraulic fracturing, and deeper release of pore fluid along the plate interface.