Geology Reports⌕ Search

SEARCH · Geology Reports

Results for “Arizona Geological Survey Contributed Report”

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.

1,128 records · Page 26Linked to original sources
USGS2025

A comprehensive geologic framework of the National Crustal Model for seismic hazard studies in the conterminous United States

A three-dimensional (3D) geologic framework has been developed for the conterminous United States (U.S.) as part of the U.S. Geological Survey National Crustal Model to enhance seismic hazard modeling. The geologic framework is created from geologic maps and multiple subsurface geologic unit boundaries including the base of the Miocene, Cenozoic, Phanerozoic, and the Mohorovičić discontinuity. Modifications are made to surficial geologic maps to remove discontinuities across state and country borders. The subsurface distribution of rock type and age is extrapolated from the surface, seeded with subsurface geologic information, and constrained by a map of basement geology. The framework provides the basis for estimates of subsurface seismic velocity and density that is needed to improve estimates of earthquake ground shaking and seismic hazard. The present framework greatly expands and updates a previously published 3D geologic framework of the western part of the U.S. that was itself a first-of-its-kind digital 3D portrayal of the nation.

conterminous United States↗
USGS2026

Widespread landslide activity in an extreme wet season and implications for regional sediment management, eastern San Francisco Bay area, California

Watershed sediment production is expected to increase in a warmer future with more extreme rain, with cascading effects throughout drainage and sediment-transport networks. This study investigated landscape-scale sediment movement in the eastern San Francisco Bay area, California, USA, during the extreme 2016–2017 wet season that brought major rainfall, landslides, and flooding. Mapping 8,928 landslides across a 1,050-km 2 study area revealed new sediment yield of 510–956 t/km 2 , equivalent to denudation of 193–361 mm/ky. These results correspond closely to long-term denudation rates in the northern and central California Coast Ranges, indicating that mass wasting in very wet years dominates long-term sediment mobilization. However, due to long residence times in drainage networks, the 2017 landslides contributed at most ∼1%–2% of the estimated locally derived fluvial sediment transport to San Francisco Bay. Although the amount of sediment mobilized did not threaten municipal water supplies, small rangeland impoundments in this mixed-use landscape lost storage capacity to new sedimentation. Considering regional sediment supply and demand, even the exceptionally large sediment delivery in an extreme wet year cannot meet the need for sediment to accrete tidal wetlands in the bay. To keep pace with rising sea levels, this abnormally high terrestrial sediment input would need to occur in 50 of the next 75 years, an unlikely occurrence due to the prevalence of recent drought years. Shoreline protection and restoration in the bay would need additional sources of sediment, such as through management of dredged sediment through beneficial-reuse programs.

California↗
USGS2026

Quantifying methane emissions from a rich fen with uncrewed aircraft systems in boreal Alaska

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.

JGR Atmospheres↗
USGS2024

Pesticide occurrence in shallow groundwater in three regions of agricultural land use: Baldwin County, the Wiregrass region, and the Tennessee River valley region of Alabama, 2009–20

As part of a cooperative investigation between the U.S. Geological Survey and the Alabama Department of Agriculture and Industries, a network of 22 groundwater wells were sampled from 2014 through 2020 for about 230 pesticide and pesticide degradate compounds. Wells were located in three regions of intensive agricultural land use in Alabama: Baldwin County, the Wiregrass region, and the Tennessee River valley region. Metolachlor sulfonic acid, a degradate of the herbicide metolachlor, was the most frequently detected compound, occurring in about 70 percent of the samples. Three other compounds, metolachlor, atrazine, and 2-chloro-4-isopropylamino-6-amino-s-triazine, were also detected in more than half of the samples. Metolachlor and its degradates accounted for 33 of the 50 greatest compound concentrations study-wide, including the maximum pesticide concentration across all compounds (62,500 nanograms per liter). The frequency and magnitude of detections of many specific pesticide compounds varied among the three regions, but all detected pesticide concentrations were well below the U.S. Environmental Protection Agency maximum contaminant levels and applicable human health benchmarks. Sample results were combined with results of previous (2009–13) sampling to provide a continuous time-series of data for 2009–20. More than half of the 289 pesticide compounds analyzed during 2009–20 were not detected in any samples. Only four compounds were detected at great enough frequency throughout the 10 sampling years to evaluate patterns of change through time. Metolachlor and its degradate, metolachlor sulfonic acid, were frequently detected in all regions. Atrazine and its degradate, 2-chloro-4-isopropylamino-6-amino-s-triazine, were also detected in wells from all regions, but the variability and magnitude of concentrations were greatest in the Tennessee River valley region. No apparent temporal pattern in concentrations was found.

Alabama↗
USGS2026

A roadmap for identifying and interpreting physical processes and national water model prediction bias associated with baseflow index regimes across the contiguous United States

Understanding how groundwater–surface water interactions shape streamflow variability is critical for diagnosing low flow behavior and prediction bias in continental scale hydrologic models. We present a process informed framework that links observed baseflow (BF) dynamics, watershed attributes, and National Water Model (NWM) performance across the contiguous United States. Using daily observed streamflow from 797 reference quality streamgages, we developed monthly baseflow index (BFI) signatures using a streamgage specific, calibrated digital filter. Hierarchical clustering of these signatures identified seven distinct BFI regimes capturing regional and seasonal variability. We evaluated NWM v3.0 retrospective streamflow performance within each regime using multiple hydrograph and flow duration curve-based metrics. Model skill varied systematically across regimes: mixed flow systems were simulated most accurately, while predominantly BF dominated and quickflow dominated regimes exhibited substantially poorer performance. Across nearly all regimes, the NWM underestimated observed BFI magnitude and frequently failed to reproduce seasonal BF patterns, indicating systematic biases in simulated low flow contributions. To relate these regimes to potential process controls, we trained a Random Forest classifier using static watershed attributes and applied Shapley Additive Explanations to identify features most strongly associated with each regime. Results highlight regionally varying influences, including the dominant role of snow fraction and seasonal runoff timing in snow dominated basins and the importance of evapotranspiration and aridity in quickflow dominated systems. Collectively, these findings demonstrate how hydrologic signatures combined with interpretable machine learning can diagnose regime specific model biases and generate process-based hypotheses about limitations in large scale hydrologic prediction systems.

contiguous United States↗
USGS2025

Megafaunal community structure on ferromanganese and phosphorite hardgrounds in the Southern California Borderland

The Southern California Borderland (SCB) is a topographically complex region on the active continental margin that hosts varied hardground habitats, including ferromanganese (FeMn) crusts and phosphorites, marine minerals being considered for resource extraction. The SCB is influenced by seasonal upwelling and terrestrial inputs, and has a well-defined oxygen-minimum zone (OMZ). We analyzed megafaunal community composition, density, and diversity across 41 ROV video transects at 10 SCB hard substrate sites spanning depth (378–2765 m), temperature (1.79–7.96 °C), and oxygen (3.91–105.67 μM) gradients, with varying rock types (FeMn crusts, phosphorites, other rocks). We counted 32,426 individuals representing 146 unique taxa over a total area of 21,935 m 2 . Echinodermata contributed 48 % of the total abundance; Cnidaria 24 %; Porifera 13 %; Annelida 6 %; Arthropoda 6 %; Chordata 2 %; Mollusca, Foraminifera, and Hemichordata <1 %. Megafauna communities showed clear heterogeneity, with density, diversity, and community composition varying among transects within sites and across sites. Rock type recorded no influence on megafaunal density, but habitats that included FeMn crusts had higher diversity and distinct taxonomic composition. Location, depth, and oxygen were the most important variables explaining variance among megafaunal communities, with distinct communities formed at deeper depths and within the OMZ. This study provides baseline information about the megafauna inhabiting SCB hardgrounds, particularly those associated with minerals considered for their resource potential. It extends existing documentation of such communities to greater depths, reveals increased representation of sponges and corals below 1000 m, and highlights the variability within and among deep-sea megafaunal communities, offering context for informed policy decisions.

Southern California borderland↗
USGS2025

Landsliding follows signatures of wildfire history and vegetative regrowth in a steep coastal shrubland

Five years after the deadly and destructive 9 January 2018 Montecito debris flows (Santa Barbara County, California, USA), an atmospheric river storm on 9 January 2023 triggered widespread landsliding that affected many of the same drainages in the Santa Ynez Mountains. Using high-resolution aerial imagery, we identified >10,000 landslides over an ∼160 km 2 area. Most of the landslides were shallow (<1 m in depth) translational debris slides that initiated on steep (∼40°), south-facing hillslopes, with the highest incidence of landsliding in a sandstone-dominated bedrock unit. The landslides mobilized into debris flows and delivered substantial quantities of sediment downstream, which contributed to costly infrastructure impairments. We detected order-of-magnitude differences in landslide density across the study area that could not be attributed to variations in geomorphology (topographic aspect and slope), geology (bedrock type), or hydrology (seasonal antecedent rainfall, peak hourly storm rainfall intensity, total storm rainfall), which are usually considered relevant factors for shallow landsliding. Rather, we found that vegetation regrowth following wildfire was likely a relevant factor associated with the highly variable landslide densities. Hillslopes with less and different types of vegetation regrowth after fire appear to have been more susceptible to shallow landslides. We identify a possible vegetation control on postfire landsliding, which highlights an opportunity for hypothesis testing using more advanced techniques to track the evolution of vegetation cover and vegetation type in steep shrubland environments following wildfire.

California↗
USGS2024

Groundwater flow model for the Des Moines River alluvial aquifer near Des Moines, Iowa

Des Moines Water Works (DMWW) is a regional municipal water utility that provides residential and commercial water resources to about 600,000 customers in Des Moines, Iowa, and surrounding municipalities in central Iowa. DMWW has identified a need for increased water supply and is exploring the potential for expanding groundwater production capabilities in the Des Moines River alluvial aquifer, where it operates two radial collector wells (RCWs). The U.S. Geological Survey, in cooperation with DMWW, completed a study of the Des Moines River alluvial aquifer and interactions of the RCWs with the aquifer; no previously published model has included the existing well locations, which is the focus of this model. A conceptual and numerical groundwater flow model have been developed to characterize the Des Moines River alluvial aquifer under existing conditions, to simulate water levels observed in the RCWs, and to provide publicly accessible hydrologic data and research that advance understanding of the regional hydrologic system and can potentially be used in the future to evaluate groundwater production scenarios. Model performance was assessed by comparing observed and simulated groundwater levels that included water level elevations, water level changes, water level inequality observations, surface water streamflow, and change in surface water volume from upstream to downstream. Water table elevation in the aquifer layers is on average slightly overestimated with average absolute value error less than 1.5 meters at both RCWs and less than 2.5 meters for all observation wells in the alluvial aquifer layers. The model also accurately simulated water tables greater than the RCW design minimum (a water level threshold at which RCW pumping is reduced) in all timesteps for which water level observation data existed. Water table elevation error was higher in other model layers that were not the focus of the study, and the model did not accurately match streamflow targets.

Iowa↗
USGS2024

Simulation of groundwater flow in the Long Island, New York regional aquifer system for pumping and recharge conditions from 1900 to 2019

The U.S. Geological Survey has developed a transient, groundwater-flow model that simulates hydrologic conditions in the Long Island aquifer system as part of an ongoing (since 2016) multiyear, cooperative investigation with the New York State Department of Environmental Conservation. The goals of this investigation are to assist stakeholders and resource managers to evaluate the response of the hydrologic system to changes in future hydraulic stresses. Responses in the hydrologic system include changes in water levels in the hydrogeologic units; discharge to streams, coastal waters, and subsurface infrastructure; and the extent of saline groundwater in the aquifers. Hydraulic stresses include future water-supply management and changes in land use and infrastructure. The numerical model synthesizes a diverse set of physiographic, geologic, climatic, land-use, and historical population, water use, and infrastructure data to physically represent the Long Island aquifer system from land surface to bedrock and to simulate annual hydrologic conditions between 1900 and 2019. A three-dimensional hydrogeologic framework was developed from existing and recently collected borehole geologic and geophysical data collected as part of a companion drilling program. Water-transmitting properties of the principal aquifer sediments were defined in three dimensions from new and existing lithologic logs. The distribution of recharge from precipitation was estimated from landscape characteristics and climate data. Anthropogenic recharge from wastewater, leaky infrastructure, and storm runoff were estimated from population, infrastructure, and pumping data. Water-use data, including well locations, depths, and pumping rates, were obtained from historical sources and records and used to estimate pumping stresses continuously in time and space, at an annual average time scale. The data were incorporated into a three-dimensional numerical model using the U.S. Geological Survey finite difference modeling code MODFLOW 6; the model encompassed all of Long Island and surrounding surface waters and simulated historical hydrologic conditions from 1900 to 2019. The calibration process involved trial and error adjustments using prior knowledge to improve general fit to observations followed by an inverse calibration to update and optimize input parameters, using an iterative ensemble smoother algorithm implemented in PEST++ version 5.0. This resulted in a model that generally was in good agreement with observed, dynamically varying hydrologic conditions from 1900 to 2019. The calibrated model was used to develop two base-case models for scenario testing of future, hypothetical conditions where one represented average-annual conditions, and one represented average-seasonal conditions from 2010 to 2019. The model representing average-annual conditions was modified further to represent an alternate sea-level position of 6 feet above the North American Vertical Datum of 1988, and the model representing average-seasonal conditions was modified to represent the average seasonal effects of a 5-year drought imposed upon current hydrologic conditions. Recharge is the sole source of water to the aquifer system; groundwater discharges to coastal water and streams and is withdrawn by pumped wells. Model-estimated annual recharge ranged from about 11 inches in 1965 to 41 inches in 1983. On average, from 2010 to 2019, about 23 percent of water was pumped from wells, and about 47 and 27 percent discharged to coastal waters and streams, respectively; the remaining 4 percent was water that moved into storage in the aquifer matrix. Water levels on Long Island vary naturally during time in response to changes in recharge; the amount of variation is largest in the interior of the island, in areas with highest water table altitudes near groundwater divides and lowest near streams and the coastal waters. The total range of water table altitudes on Long Island between 1900 and 2019 ranged from near 0 to more than 70 feet in western parts of Long Island. The largest range in altitudes is in New York City and is associated with areas of large historical withdrawals between the 1920s and the late 1980s. Water table altitudes generally varied by less than 10 feet in eastern Suffolk County, where the aquifer is under more natural conditions. Saltwater intrusion is of great concern on Long Island, particularly in western Long Island where both the unconfined and confined parts of the aquifer system have been intruded in response to large-scale groundwater withdrawals; however, the volume of freshwater in the islandwide aquifer system only has changed by about 5 percent between 1900 and 2019. The decadal change in the freshwater volume was largest during the early and mid-20th century, corresponding to the largest historical pumping, but that volume change did not exceed 1 percent. The negligible change in freshwater volume suggests that saltwater intrusion as of 2019 was limited at an islandwide scale but continues to occur in local areas of Queens and Nassau Counties, adversely affecting current water supplies and limiting future water supplies for affected communities. The regional groundwater model developed for this investigation is a tool that can be used to help determine the viability of current and future water supplies at a regional scale and can be used to support development of additional models at finer scale to support more focused assessments of groundwater sustainability.

New York↗
USGS2026

Water use in Louisiana, 2020

The U.S. Geological Survey (USGS), in cooperation with the Louisiana Department of Transportation and Development, collected water-withdrawal and water-use data from a 2020 inventory of water withdrawals in Louisiana. In 2020, approximately 8,700 million gallons per day (Mgal/d) of water was withdrawn from groundwater and surface-water sources in Louisiana, which represented a 0.22-percent decrease from 2015. Total groundwater withdrawals were about 1,900 Mgal/d, an increase of 7.1 percent from 2015, and total surface-water withdrawals were about 6,800 Mgal/d, a decrease of 2.1 percent from 2015 to 2020. Total water withdrawals, in million gallons per day, in 2020 for the various categories of use were as follows: public supply, 720; industry, 2,100; power generation, 4,100; rural domestic, 39; livestock, 7.0; rice irrigation, 930; general irrigation, 250; and aquaculture, 590. From 2015 to 2020, Louisiana’s total withdrawals for public supply increased by 1.4 percent, industry decreased by 2.3 percent, power generation decreased by 4.9 percent, rural domestic decreased by 1.2 percent, livestock increased by 11 percent, rice irrigation increased by 13 percent, general irrigation increased by 12 percent, and aquaculture increased by 20 percent. About 51 percent (approximately 960 Mgal/d) of all groundwater withdrawn was from the Chicot aquifer system and 24 percent (approximately 450 Mgal/d) was withdrawn from the Mississippi River alluvial aquifer. Since 2015, withdrawals from the Chicot aquifer system increased by 13 percent, and withdrawals from the Mississippi River alluvial aquifer increased by 18 percent. About 72 percent (4,900 Mgal/d) of all surface water withdrawn was from the Mississippi River main stem. This value represents a 1.1-percent decrease in withdrawals from 2015 to 2020. All water-withdrawal and water-use data presented in this report should be considered estimates. Because of rounding, totals and percentages presented in the tables, figures, and text in the report may differ slightly from totals or percentages calculated individually.

Louisiana↗
USGS2026

Evaluation of stream capture related to groundwater pumping, middle Humboldt River Basin, Nevada

Historical, future, and potential stream capture from groundwater pumping in the middle Humboldt River Basin (MHRB), Nevada, is estimated using a calibrated numerical groundwater flow model. The model was developed to estimate (1) stream capture, which is the change in flux between the groundwater system and the Humboldt River and tributaries, and (2) change in streamflow, which is the change in streamflow estimated for the Imlay gage on the Humboldt River (U.S. Geological Survey streamgage 10333000). Historical stream capture for water years (WYs) 1961–2015 is estimated using recorded and estimated groundwater pumping during that period. Future (predictive) stream capture was based on historical stresses (WYs 1961–2015) using a scenario that simulated non-mine pumping from WY 2015 at a uniform rate for 100 years into the future. Potential stream capture throughout the middle Humboldt River Basin from groundwater pumping during varying durations of time are presented in a series of capture maps. Maps also are presented that show the potential to capture from groundwater evapotranspiration, as well as the storage changes for pumping duration of 100 years. Estimates of historical stream capture from the mainstem Humboldt River during the early 1960s are less than 400 acre-feet per year (acre-ft/yr) when groundwater withdrawals and pumping rates were relatively small compared to more recent times. In the late 1980s and early 1990s, groundwater withdrawals increased and estimated historical stream capture also increased from about 4,000 acre-ft/yr in the late 1980s and early 1990s to as much as 18,800 acre-feet (acre-ft) in WY 1998. In WY 2015, estimated historical stream capture declined to about 13,000 acre-ft because of decreasing groundwater withdrawals and lower streamflow during the drought of WYs 2012–15, resulting in less stream water available for capture. Stream capture was estimated for 100 years into the future based on WY 2015 non-mine pumping rates and mine-dewatering activity through WY 2015. Stream capture is forecast to increase to about 23,000 acre-ft/yr, and streamflow in the Humboldt River could decrease by as much as 19,000 acre-ft/yr. Pumping for mine-dewatering and the associated discharge of that water affects streamflow in the Humboldt River at Imlay, Nevada (U.S. Geological Survey streamgage 10333000). Historically, from WYs 1991 to 2015, streamflow was greater at Imlay gage during active mine-dewatering from mine-water discharge operations and increased by as much as 105,000 acre-ft in WY 1998. The increase was attributed mostly to the discharge of groundwater from mine-related dewatering operations directly into the mainstem Humboldt River or its tributaries, with some of this increase associated with return flows from discharge to rapid infiltration basins. Results indicate that streamflow at Imlay gage is expected to decrease by as much as 1,600 acre-ft/yr 30 years after mine-related pumping and discharge are discontinued. The streamflow reductions at the Imlay gage are expected to then decrease to around 500 acre-ft/yr, 100 years after mine-related pumping and discharge are discontinued. Potential capture maps were produced for pumping durations of 10, 25, 50, and 100 years. Capture map results indicate that areas of greater potential stream capture occur adjacent to the Humboldt River and for upstream tributaries areas north of the Humboldt River.

Nevada↗
USGS2026

Numerical modeling of Late Pleistocene to Holocene earthquake-induced progressive rock slope damage in the paraglacial Serpentine valley, Prince William Sound, Alaska

Landslides in deglaciating fjords pose a potential tsunami threat to nearby communities; however, processes contributing to long-term progressive rock damage and landslide conditioning remain poorly constrained in paraglacial settings. Here, we analyse the role of earthquake-induced rock mass damage over late Pleistocene to Holocene time scales, as a conditioning factor for modern landslides, using distinct element numerical modelling to assess spatial and temporal patterns of fracture propagation influenced by varying glacier thickness. Conceptualised numerical models were parameterised by in situ rock mass, glacial and topographic conditions in Serpentine valley, located in Prince William Sound, Alaska, where several large landslides are actively developing along the western valley wall. Results show that although rigid glacier buttressing reduces co-seismic rock mass damage, it does not suppress it completely, and damage occurs both above and below the glacier surface elevation. We demonstrate that topography, and especially steep slopes with topographic convexity, as well as preexisting damage of joint networks and faults inherited from tectonic and exhumation induced stresses, exert primary control on the location of new co-seismic damage. Simulations representing a simplified deglaciation sequence over the past 25 ka, with a series of 10 evenly spaced earthquakes, generated rock mass damage patterns that qualitatively match in situ landslide structural and kinematic observations at one instability in Serpentine valley. Our conceptual study helps clarify the role of repeated seismicity over glacial timescales as a long-term conditioning process for paraglacial rock slope failure, highlighting spatial and temporal patterns of progressive damage accumulation, with outcomes relevant for modern landslide hazard assessment.

Alaska↗
USGS2024

Perchlorate, metals, organic compounds, and lead isotopes in groundwater, surface water, shallow groundwater, and soil within and near the Middleton Municipal Airport–Morey Field (C29), Middleton, Wisconsin, 2022

The Middleton Municipal Airport–Morey Field (C29) is in the City of Middleton and adjacent to the towns of Middleton and Springfield, Wisconsin. Nearby homes in the towns rely on private drinking water wells, and residents are concerned about the potential contamination of groundwater and surface water by airport activities, including flights by small aircraft that use leaded aviation fuel and a fireworks display in July 2021. The U.S. Geological Survey, in cooperation with the Town of Middleton, completed a study in 2022 to characterize the occurrence and sources of perchlorate, metals (including lead), and organic compounds in samples of groundwater, surface water, shallow groundwater, and soils within and near the airport. Lead isotopes were also measured to determine sources of lead by comparing samples to environmental references. Magnitudes of concentrations from samples of water and soil collected in 2022, and their spatial patterns across site locations, indicate the fireworks display in July 2021 was a likely source of perchlorate and metals in the airport study area. The highest perchlorate concentration was measured in surface water at the southeastern corner of the airport near the fireworks launch site; the highest concentrations of fireworks-associated metals were measured in shallow groundwater near the same location. Fireworks were not the only possible source of perchlorate and metals in the airport study area because both were also detected upgradient and away from the fireworks launch site. Ratios of lead isotopes indicate that lead measured in water and soil within the airport study area was primarily sourced from background atmospheric lead deposition or Wisconsin galena lead ore. However, two groundwater samples (one upgradient and one downgradient from the airport; both with concentrations less than 1 microgram per liter) had isotopic signatures matching leaded aviation fuel sold at the airport.

Wisconsin↗
USGS2025

Geochemical and hydrological investigations of historical data collected at the Lee Acres Landfill and Giant Bloomfield Refinery, New Mexico, 1985–2020

The Lee Acres Landfill and Giant Bloomfield Refinery are adjacent properties near the City of Farmington, New Mexico, each having undergone monitoring and remediation related to historical site activities. At the landfill, site cleanup has included the installation of a capillary barrier over former liquid waste lagoons and periodic monitoring of groundwater elevations and groundwater quality. At the refinery, remediation has focused on several petrochemical and crude oil release areas and included soil excavation, groundwater treatment, and regular monitoring of groundwater elevations and quality. Groundwater at both sites has higher concentrations of volatile organic compounds and trace metals than background aquifer concentrations. In 2022, the U.S. Geological Survey compiled the Lee Acres-Giant Bloomfield Refinery Database (LAGBRD), which contains publicly available groundwater-elevation data and organic and inorganic groundwater-quality data from both sites, spanning from 1985 to 2020. Data from the LAGBRD and precipitation data from other sources were used to better understand the cause of relatively high manganese concentrations observed in some groundwater wells at the site through comparison of groundwater chemistry to chemical end members, interpretation of spatial and temporal patterns in the groundwater chemistry, and interpretation of groundwater flow properties. In this study, elevated chloride concentrations in groundwater downgradient from the landfill have been attributed to landfill leachate based on the temporal and spatial variability of chloride concentrations and chloride-to-bromide ratios. Installation of a capillary barrier and surface-water runoff controls at the landfill in 2005 appears to have altered infiltration patterns at that site, resulting in a decrease in chloride at some wells but an increase in chloride and dissolved manganese at others. The timing and relation among groundwater elevation, chloride concentration, and manganese concentration suggest that leachate stored in the vadose zone provides a continued source of contamination to groundwater.

New Mexico↗
USGS2025

Statewide cumulative human health risk assessment of inorganics-contaminated groundwater wells, Montana, USA

Across the United States, rural residents rely on unregulated and generally unmonitored private wells for drinking water, which may pose serious health risks due to unrecognized contaminants. We assessed the nature, degree, and spatial distribution of cumulative health risks from inorganic contaminants in groundwater. Our analysis included nearly 84,000 data points from 6500+ wells, across 51 of Montana's 98 watersheds, using a public groundwater database. We compared a drinking water screening level cumulative risk assessment (CRA) for inorganics based on the U.S. Environmental Protection Agency (EPA) protective health thresholds (Maximum Contaminant Level Goals, Health Advisories [MCLG-HAs]) to a CRA based on EPA public supply enforceable Maximum Contaminant Levels (MCLs). Based on median concentrations of 19 inorganics (antimony, arsenic, barium, beryllium, boron, cadmium, chromium, copper, fluoride, manganese, molybdenum, nickel, nitrate, lead, selenium, strontium, thallium, uranium, zinc), 75% of watersheds had MCLG-HA-based cumulative risk values > 1.0; arsenic and uranium contributed the most risk, followed by strontium, fluoride, manganese and boron. Hence, this screening level (Tier I) CRA indicated widespread potential for unrecognized human health risk to private well users from inorganic contaminants considering both carcinogenic and non-carcinogenic risks. Sensitivity analysis showed that benchmarks applied (MCLG-HAs versus MCLs) exerted the largest control on results. Our findings identify priority regions for Tier 2 risk assessments to elucidate local sources and distributions of geogenic versus anthropomorphic contaminants. Our study is the first statewide assessment of cumulative health risk from groundwater that we are aware of, and results support increased statewide drinking water education and testing to reduce human health risks from contaminated private well water.

Montana↗
USGS2025

What 25+ years of "Did You Feel It" intensities tell us about shaking in California

“When will the Big One happen?” is a question that people often have for earthquake scientists. But while waiting for the “Big One” to occur, people will usually experience frightening or damaging shaking from multiple relatively smaller‐magnitude earthquakes. Given this context, it raises the question: “Where does most of the damage come from?” Could smaller, yet more frequent, earthquakes account for the majority of reported impactful shaking? To explore this question, we consider reports of earthquake damage and felt shaking experiences from a catalog of community‐collected intensity values from the U.S. Geological Survey’s “Did You Feel It?” system. Comparing these intensities to expectations from a ground‐motion model, we find that earthquakes of magnitudes smaller than expected are responsible for most reported intensities of community decimal intensities (CDI) 4.5 and above (moderate and higher shaking intensity levels). (Here “expected value” is meant in its mathematical sense of the mean or equivalently median shaking intensity.) We also present a regional analysis of observed earthquake shaking for specific areas, detailing the maximum intensity experienced within a predetermined area. We identify several instances of M < 4.5 events that generated maximum intensities of CDI > 5 in regions in California surrounding the cities of Eureka, San Francisco, Los Angeles, and San Diego. Our results motivate the need to include smaller‐magnitude earthquakes in communications about earthquake hazard and risk reduction.

California↗
USGS2024

Mapping Arundo donax (Arundo cane) with multispectral imagery before, during, and after herbicide treatment along the Rio Grande in Webb County, Texas, 2020–21

Arundo donax , commonly called Arundo cane, giant reed, or Carrizo cane, is an invasive bamboo-like perennial grass common in riparian areas throughout the southwestern United States. In Texas, not only does it negatively affect riparian ecosystems, but it has also become a problem for border security because it reduces visibility along the Rio Grande. To address these problems, in 2015 the Texas State Soil and Water Conservation Board was authorized by the Texas State Legislature to develop a program to eradicate Arundo cane along the Rio Grande. In 2020, the Texas State Soil and Water Conservation Board applied imazapyr and glyphosate herbicides along a 19.3-kilometer reach of the Rio Grande, northwest of Laredo, Texas. The U.S. Geological Survey, in cooperation with the Texas State Soil and Water Conservation Board and the Webb Soil and Water Conservation District, used WorldView-3 Standard high-resolution satellite imagery to map Arundo cane extent along the reach before , during , and after the herbicide-treatment period on June 30, 2020, September 26, 2020, and May 7, 2021, respectively. A maximum likelihood supervised classification analysis was computed on the images to map the spatial extent and estimate the area covered by Arundo cane. The estimated area covered by Arundo cane in the before classification was 1,282,000 square meters, in the during classification was 1,064,000 square meters, and in the after classification was 1,108,000 square meters. The qualitative comparison of the three images shows that there was an overall decrease in vegetation classified as Arundo cane throughout the study area.

Texas↗
USGS2025

Evaluating Laramide orogenesis via flexural basin response in the San Juan basin, New Mexico and Colorado

A challenge in interpreting the location, timing, and magnitude of ancient orogenic events is that ongoing uplift and erosion in the hinterlands often destroys much of the primary record of these events. However, basin-thickness patterns in the sedimentary record can provide complimentary evidence of uplift via flexural effects. Here, we deploy well-log correlation, isochores, basin modeling, flexural modeling, and subcrop mapping to evaluate the Late Cretaceous to Paleogene basin response to Laramide tectonism in the San Juan basin. A wedge of upper Campanian to Maastrichtian sedimentary rock thickens from 200 to 800 meters from southeast to northwest in the basin. This pattern can be successfully simulated via flexural modeling if we infer early Laramide uplift along the northwest basin flank that produced a 0.8 km high topographic load. The Laramide unconformity bounds the top of this Upper Cretaceous sedimentary wedge and truncates progressively older strata to the east, further supporting a westward tilt of the basin. The onset of Campanian Laramide flexure may have also contributed to the profound transgression from the upper Menefee Formation to the Lewis Shale. The Paleocene isochore map displays an approximately symmetrical pattern, with thickening towards the center of the basin. This suggests the possibility of competing flexural loads. The base Eocene structure indicates an asymmetric deep on the northeast flank of the basin, providing flexural evidence of contemporaneous uplift/loading of the Nacimiento uplift and Archuleta arch; this has been modeled as ~2.1 km load height. Both Cretaceous and Paleocene sedimentary wedges are narrow, suggesting low flexural rigidity; modeled effective elastic thicknesses (EET) are 20-30 km, comparable to estimates of modern EET for the region.

Colorado, New Mexico↗