USGS Science⌕ Search

SEARCH · USGS Science

Results for “Earth Surface Processes and Landforms”

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,665 records · Page 24Linked to original sources
USGS2026

Simulation-based scenario ShakeMaps for large magnitude (MW6.5+) crustal earthquakes on the Seattle, Tacoma, and southern Whidbey Island faults, Washington, USA

Scenario ground‐motion maps based on empirical ground‐motion models (GMMs) provide a rapid and generally reliable means of estimating the amplitude and distribution of earthquake shaking. However, because GMMs are designed for broad applicability, they often rely on simplified representations of Earth structure, which can limit their accuracy in regions with complex source, path, and site effects. This can substantially impact the accuracy of predicted shaking in areas like western Washington State, where deep, interconnected basin structure exerts a strong influence on seismic‐wave propagation. In this study, we present a new suite of simulation‐based scenario ShakeMaps that characterize ground shaking from large‐magnitude ( ⁠ M W 6.5–7.5) earthquakes on the Seattle, Tacoma, and southern Whidbey Island faults. These maps are developed using results from recent 3D wave propagation simulations ( Stone et al. , 2022 , 2023 , 2025 ) that incorporate realistic rupture geometries, variable slip distributions, and a regional 3D seismic velocity model with shallow soils. Broadband ground motions are estimated by combining the low‐frequency (<1 Hz) deterministic seismograms from these studies with high‐frequency (1–10 Hz) stochastic seismograms. Simulated ground motions are corrected to account for the enforced minimum shear‐wave velocity and nonlinear site response. The resulting ShakeMaps represent median ground‐shaking estimates derived from multiple rupture scenarios with varying slip distributions and hypocenter locations for each fault. To extend ShakeMap coverage beyond the simulation domain (i.e., into eastern Washington, northern Oregon, and southwestern British Columbia), we scale GMM‐based ground‐motion estimates using amplification patterns observed in the simulations. These new ShakeMaps reveal the substantial influence of deep basin structure on shaking intensity, underscoring the importance of considering crustal structure complexity in regional hazard assessments for the Pacific Northwest.

Washington↗
USGS2025

An evaluation of the effects of different deicing salt application rates on three watersheds in Essex County, New York

The U.S. Geological Survey, in cooperation with the New York State Department of Transportation, evaluated the effects of different deicing salt application rates on surface water, groundwater, and highway runoff quality near State highways in northern New York. Three reaches of State highways were tested with different deicing treatments between October 2019 and November 2022: a salt-sand mixture (Treatment A), a salt mixture applied at a lower rate (Treatment B), and a control mixture consistent with typical deicing salt amounts and application rates. Data on pavement conditions and the quality of surface water, highway runoff, and groundwater were collected. Surface electromagnetic data were also collected. Surface-water and groundwater quality downgradient from the State highways were compared with water quality at upgradient locations. The percentage of snow or ice coverage was used to evaluate the effectiveness of the salt applications. This report provides an overview of the transport of deicing salt. The Treatment B watershed had deicing mixture applied more frequently than other highway reaches, which caused it to have the highest annual total chloride application. Despite differences in chloride application, flow-weighted mean chloride concentrations in highway runoff were comparable across treatments. Chloride concentrations were elevated in surface water and groundwater downgradient from highways relative to chloride concentrations upgradient from highways. A chloride mass balance, calculated for one treatment watershed, indicated that groundwater affected by legacy deicing practices may be contributing additional chloride to surface water. Spatial patterns from electromagnetic surveys show a shallow saline plume alongside the highway in that area. Differences in winter severity and pavement-surface conditions drove deicing salt applications in the treatment areas. This study found that several factors affect chloride loads in the watersheds, including variable winter conditions, adaptive snow and ice management, legacy management practices, and area-specific aquifer and groundwater conditions.

New York↗
USGS2026

Salinas Valley integrated hydrologic and reservoir operations models, Monterey and San Luis Obispo Counties, California

The area surrounding the Salinas Valley groundwater basin in Monterey and San Luis Obispo Counties of California is a highly productive agricultural area, contributes substantially to the local economy, and provides a substantial portion of vegetables and other agricultural commodities to the Nation. This region of California provides about half of the Nation’s lettuce, celery, broccoli, and spinach each year. Thus, this agricultural area provides substantial volumes of agricultural products not just for California but for the United States. Changes in population and increased agricultural development, which includes a shift toward more water-intensive crops, and climate variability, have put increasing demand on both surface-water and groundwater resources in the valley. This situation has resulted in water management challenges in the Salinas Valley that generally relate to the distribution of the water supply throughout the basin. Where and when the water is present in the surface and subsurface does not coincide with where and when the water is needed. Historically, to deal with the distribution issue, water has been used conjunctively in the valley. Conjunctive use is a water management strategy that coordinates surface-water and groundwater use to maximize water availability. Groundwater is used throughout the Salinas Valley to meet water demands when surface-water supplies are insufficient. The availability of surface water is constrained by climate. Precipitation and streamflow vary seasonally and year to year. Although there are two reservoirs in the Salinas Valley to capture and store water during wet periods, the only conveyance of reservoir water to coastal agricultural areas is the Salinas River. Increasing demand for groundwater and surface-water resources throughout the Salinas Valley has resulted in undesirable effects from unsustainable water use, such as surface-water depletion, groundwater-level declines, storage depletion in the principal aquifers, and seawater intrusion. To address these escalating issues, local communities, water management agencies, and groundwater sustainability agencies are evaluating how to sustainably manage both their surface-water and groundwater resources. To meet water demands and reduce the undesirable effects of unsustainable water use, continued conjunctive management of surface water and groundwater would ideally incorporate strategies to deal with increases in demand and climate variability. To evaluate the challenging water management issues in the Salinas Valley, the U.S. Geological Survey, Monterey County Water Resources Agency, and the Salinas Valley Basin Groundwater Sustainability Agency developed a comprehensive suite of models that represent the Salinas Valley hydrogeologic system called the Salinas Valley System Model. The geologic framework is known as the Salinas Valley Geologic Framework and was developed to characterize the subsurface using various topographic and geologic data sources, including information on hydrogeologic units, their surfaces and extents, geologic structures, lithology, and elevations from borehole data and cross sections, as well as details on faults and existing models. The surface-water model is called the Salinas Valley Watershed Model and simulates the Salinas River watershed. Monthly surface-water inflows into the integrated hydrologic model domain were simulated using the Salinas Valley Watershed Model. The historical model uses historical climate data, water and land use data, and reservoir releases to simulate agricultural operations, including landscape water demands, diversions, and reclaimed wastewater. The operational model adds an embedded reservoir operations framework to the simulation of the historical model that allows specified operational rules to simulate reservoir releases and changes in reservoir storage. The operational model assumes current reservoir operations and constant land use, which differs from historical conditions. Thus, the operational model is a hypothetical baseline model that can be used by local water managers to evaluate and quantify potential benefits of water supply projects. Together, the geologic framework, watershed, historical, and operational models form a tool that can be used to simulate irrigated agriculture and associated reservoir operations of the integrated hydrologic system of the Salinas Valley.

California↗
USGS2026

Recent scientific contributions by the U.S. Geological Survey in the San Francisco Bay and Sacramento–San Joaquin Delta Estuary

Introduction The San Francisco Bay and Sacramento–San Joaquin Delta Estuary (hereafter, Bay-Delta) is the largest estuary on the West Coast of the United States. The Bay-Delta covers more than 1,600 square miles and drains a watershed of more than 75,000 square miles, which is greater than 40 percent of California. The region surrounding the Bay-Delta is home to about 10 million people, and its habitats (fig. 1) support more than 800 plant and animal species. The waterways of the Bay-Delta are the central hub of California’s extensive freshwater delivery system, supplying water to more than 27 million Californians and 4 million acres of farmland in the Central Valley. The U.S. Geological Survey (USGS) is the primary science agency of the U.S. Department of the Interior, providing a broad range of Earth, water, biological, and mapping data and expertise to inform natural resource management across the country. This fact sheet focuses on research conducted by the USGS in the Bay-Delta region, mostly within the past 5 years. The fact sheet is organized across five major themes: water flow and water quality, fish and wildlife, wetland restoration, invasive species, and hazards. In each of these areas, the USGS works closely with Federal, State, and local agencies; academia; and non-governmental organizations to provide objective, evidence-based science. The data and knowledge gained from USGS research and monitoring are publicly available and are routinely used by agencies, including the U.S. Bureau of Reclamation, U.S. Army Corps of Engineers, U.S. Fish and Wildlife Service, California Department of Water Resources, California Department of Fish and Wildlife, California State Water Resources Control Board, Delta Stewardship Council, and many additional organizations and stakeholders to inform management of the Bay-Delta region and its resources.

California↗
USGS2024

The EnMAP spaceborne imaging spectroscopy mission: Initial scientific results two years after launch

Imaging spectroscopy has been a recognized and established remote sensing technology since the 1980s, mainly using airborne and field-based platforms to identify and quantify key bio- and geo-chemical surface and atmospheric compounds, based on characteristic spectral reflectance features in the visible-near infrared (VNIR) and short-wave infrared (SWIR). Spaceborne missions, a leap in technology, were sparse, starting with the CHRIS/PROBA and EO1/Hyperion missions in the early 2000s, and providing spectroscopy data with limited spectral coverage and/or low data quality in the SWIR. Since 2019, several countries and agencies have successfully launched a number of spaceborne imaging spectroscopy systems into orbit or deployed them on the International Space Station (ISS) such as DESIS, PRISMA, HISUI, GF-5, EnMAP and EMIT. Among these recent missions, the German Environmental Mapping and Analysis Program (EnMAP) stands for its long-term development, sophisticated design with on-board calibration, high data quality requirements, and extensive accompanying science program. EnMAP was launched in April 2022 and, following a successful commissioning phase, started its operational activities in November 2022. The EnMAP mission encompasses global coverage from 80° N to 80° S through on-demand data acquisitions. Data are free and open access with 30 m spatial resolution, a high spectral resolution with a spectral sampling distance of 6.5 nm and 10 nm in the VNIR and SWIR regions respectively, and a high signal-to-noise ratio. In this paper, we aim to present the mission's current status, coverage, science capabilities and performance two years after launch. We show the potential of EnMAP for space-based imaging spectroscopy to operate in various environments, including high and low light levels, dense forests, Antarctic glaciers, and arid agricultural areas. EnMAP enables various applications in fields such as agriculture and forestry, soil compositional, raw materials, and methane mapping, as well as water quality assessment, and snow and ice properties. The results show that EnMAP's performance exceeds the mission requirements, and highlights the significant potential for contribution to scientific exploitation in various geo- and biochemical sciences. EnMAP is also expected to serve as a key tool for the development and testing of data processing algorithms for upcoming global operational missions.

Remote Sensing of Environment↗
USGS2025

Concentration-discharge relations and transient metal loads reveal spatiotemporal variability in solute-generation mechanisms in a mine-affected watershed

Concentration-discharge (CQ) relations are commonly used to understand geochemical and hydrologic controls on the generation of solutes in watersheds. Despite the widespread application of CQ relations, this technique has been infrequently applied to acid mine drainage (AMD) sites, but the CQ framework may allow mechanistic understanding of remedial outcomes such as impoundment of water within underground mines. Results of CQ analyses and changes in metal loads in an AMD affected watershed in Colorado, USA indicate that dissolved loads increased at many individual locations following water impoundment within mine workings. Although increased loads were observed at most individual locations, these increases were offset by a large decrease in loading from the largest mine. A loading analysis that included data from an instream monitoring location showed a statistically significant decrease in Fe and Zn after bulkhead emplacement, indicating a net positive effect of bulkheads. Streams generally displayed dilution CQ patterns whereas mines and springs showed either flushing or chemostatic patterns prior to bulkheading, which transitioned to chemostatic patterns following bulkheading, indicating a transition from dynamic to equilibrium geochemical processes. Saturation indices for sulfide and secondary minerals indicated that mines and springs were near equilibrium for phases including schwertmannite, fluorite, and gypsum. Saturation indices vary through time for mines suggesting progressive leaching of sulfide minerals as the mass of available minerals in the mine workings decreases. Together, these diverse analyses provide an integrated understanding of the variability in solute generating processes in this watershed and may inform remediation plans for similarly affected sites by indicating the nature of mineralogic controls on water quality.

Colorado↗
USGS2025

Salinas Valley integrated hydrologic and reservoir operations models, Monterey and San Luis Obispo Counties, California

The area surrounding the Salinas Valley groundwater basin in Monterey and San Luis Obispo Counties of California is a highly productive agricultural area, contributes significantly to the local economy, and provides a substantial portion of vegetables and other agricultural commodities to the Nation. This region of California provides about half of the Nation’s lettuce, celery, broccoli, and spinach each year. Thus, this agricultural area provides significant volumes of agricultural products not just for California but the entire United States. Changes in population and increased agricultural development, which includes a shift toward more water-intensive crops, and climate variability, have put increasing demand on both surface water and groundwater resources in the valley. This has resulted in water management challenges in the Salinas Valley that are predominantly related to distribution of water supply throughout the basin. Where and when the water is present in the surface and subsurface does not coincide with where and when the water is needed. To deal with the distribution issue, historically water has been used conjunctively in the valley. Conjunctive use is a water management strategy that coordinates surface water and groundwater use to maximize water availability. Groundwater is used throughout the Salinas Valley to meet water demands when surface water supplies are insufficient. Availability of surface water is constrained by climate. Precipitation and streamflow vary seasonally and year to year. Although there are two reservoirs in the Salinas Valley to capture and store water during wet periods, the only conveyance of reservoir water to coastal agricultural areas is the Salinas River. Increasing demand on groundwater and surface water resources throughout the Salinas Valley has resulted in undesirable effects of unsustainable water use, such as surface water depletion, groundwater level declines, storage depletion in the principal aquifers, and seawater intrusion. To address these escalating issues, local communities, water management agencies, and groundwater sustainability agencies are evaluating how to sustainably manage both their surface water and groundwater resources. To meet water demands and reduce undesirable effects of unsustainable water use, continued conjunctive management of surface water and groundwater would ideally incorporate strategies to deal with increases in demand and a variable climate. To evaluate the challenging water management issues in the Salinas Valley, the U.S. Geological Survey, Monterey County Water Resource Agency, and the Salinas Valley Basin Groundwater Sustainability Agency developed a comprehensive suite of models that represent the Salinas Valley Hydrogeologic system called the Salinas Valley System Model. The Salinas Valley Geologic Framework was developed to characterize the subsurface using various topographic and geologic data sources, including information on hydrogeologic units, their surfaces and extents, geologic structures, lithology, and elevations from borehole data and cross sections, as well as details on faults and existing models. The Salinas Valley Watershed Model simulates the entire Salinas River watershed. Monthly surface water inflows into the integrated hydrologic model domain were simulated using the Salinas Valley Watershed Model. The historical model uses historical climate data, water and land use data, and reservoir releases to simulate agricultural operations, including landscape water demands, diversions, and reclaimed wastewater. The operational model adds an embedded reservoir operations framework to the simulation of the historical model that allows specified operational rules to simulate reservoir releases and changes in reservoir storage. The operational model assumes current reservoir operations and constant land use, which differs from historical conditions. Thus, the operational model is a hypothetical baseline model that can be used by local water managers to evaluate and quantify potential benefits of water supply projects. Together, the geologic framework, watershed, historical, and operational models form a tool that can be used to simulate irrigated agriculture and associated reservoir operations of the integrated hydrologic system of the Salinas Valley.

EarthArXiv↗
USGS2025

Nanometer-scale relationships between sedimentary organic matter molecular composition, fluorescence, cathodoluminescence, and reflectance: The importance of oxygen content at low thermal maturities

Molecular characterization of sedimentary organic matter (SOM), termed macerals, is a common goal when seeking to understand petroleum generation as well as other geologic processes in deep time. However, unambiguous measurement of discrete macerals is challenging due to the small size of organic particles in sedimentary rocks, the proximity of different organic matter types to one another, mineral-organic matter interactions, and maceral mixing that occurs during SOM isolation prior to ex situ analysis. The recent advent of infrared spectrometers capable of nanometer-scale resolution and the application of these technologies to geologic samples has enabled advances in rapid, in situ molecular characterization of SOM allowing for insights into paleoenvironmental processes, such as organic matter productivity and preservation, among others. Here we employ one such technology, optical photothermal infrared (OPTIR) spectroscopy, to map SOM functional group distributions at 500-nm resolution in a sample from the Lower Cretaceous Sunniland Limestone of the South Florida Basin. Examined fields of view include occurrences of amorphous organic matter (AOM), inertinite, micrinite, solid bitumen, telalginite, and vitrinite. OPTIR data from these macerals are compared against traditional organic petrographic data from the same organic grains including fluorescence intensity and white light reflectance as well as against cathodoluminescence response, an emerging organic petrographic approach. Maceral oxygen content (using carbonyl functional group abundance as a proxy) is observed to vary widely between maceral types but correlates strongly with fluorescence and cathodoluminescence intensity as well as against reflectance. These findings highlight the important role that oxygen content plays in determining the optical properties of SOM and further demonstrate the ability of OPTIR to discriminate subtle molecular differences between SOM types.

Organic Geochemistry↗
USGS2025

International data gaps at the Center for Engineering Strong Motion Data

The Center for Engineering Strong Motion Data (CESMD) is utilized by seismologists, engineers, and disaster management professionals in the US and has historically achieved and distributed waveforms from across the globe for significant earthquakes. The increased access to the waveforms via Web API (Application Programming Interface) offers a unique opportunity to provide the community complete datasets, sampling a variety of tectonic environments and geologic conditions, increasing the number of available ground motion records for use in ground motion models (GMMs) and improving the accuracy of earthquake engineering evaluations. The objective of this study is to programmatically identify gaps in global event data from the past decade and backfill missing data gaps at CESMD. We first compare the CESMD catalog with the Advanced National Seismic System (ANSS) Comprehensive Earthquake Catalog identifying regions and time periods where strong-motion data is limited or inadequate. To backfill datasets at CESMD for significant events, we pinpoint regions and time intervals that lack information, creating a list of events for which we’d like to obtain data. An important facet of this work is identifying the source of data and metadata across earthquake repositories around the world and integrating these data repositories into our current strong-motion data processing workflow. In parallel with these newly processed datasets, we are developing a script to produce data origination citations to include provenance and attribution information to associate with respective datasets at CESMD. We showcase our methodology for identifying and filling data gaps at CESMD using three case studies (the 2018 Anchorage Alaska earthquake sequence, seismicity associated with the 2018 Hawaiian Kilauea volcano eruption, and several earthquakes in Turkey) and then outline our strategy to apply our data gap backfilling methods on an international scale.

Conference Paper↗
USGS2026

Ground-motion simulations for the 2024 Mw 4.8 Tewksbury, New Jersey, earthquake

Ground-motion simulations of notable earthquakes in the central and eastern United States are limited and typically assume one-dimensional (1D) Earth structure. In this study, we use a three-dimensional (3D) seismic velocity model to better constrain the depth and focal mechanism of the April 5th, 2024, moment magnitude 4.8 Tewksbury earthquake and investigate the spatial variability of earthquake ground motions and the effects of nearby sedimentary basins. We perform earthquake ground-motion simulations up to 0.5 Hz using the 3D spectral-element wave-propagation solver SPECFEM3D over a region 280-km wide by 260-km long by 77-km deep. Topography and subsurface geophysical structure are assigned using the U.S. Geological Survey National Crustal Model with a minimum shear-wave velocity of 200 m/s. We use earthquake time series from 13 broadband seismic stations in the region that have a uniform azimuthal distribution and epicentral distances ranging from 76 to 131 km to compare with synthetics and explore the effects of 1D versus 3D seismic structure on focal mechanism and depth solutions. Ground-motion intensity metrics are also presented relative to the NGA-East ground-motion models (GMMs) currently used in seismic hazard assessments for the region. We find that the 3D model, which reveals a wide spatial variability of period-dependent ground motions, yields better predictions of earthquake ground motions relative to the 1D model and the NGA-East ergodic ground-motion model, with 76 percent reduction of residual variance in observed ground motions averaged over 3-, 5-, 7-, and 10-second periods. Use of the 3D model to solve for a focal mechanism yields a shallower focal depth at 4 km and a shallower east-dipping focal plane relative to the U.S. Geological Survey regional moment tensor and Global Centroid Moment Tensor. Our study demonstrates that use of 3D seismic velocity models can improve estimates of earthquake focal mechanisms, ground motions, and seismic hazard.

New Jersey↗
USGS2024

New developments at the Center for Engineering Strong-Motion Data (CESMD)

The Center for Engineering Strong-Motion Data (CESMD), an internationally utilized joint center of the U.S. Geological Survey (USGS) and the California Geological Survey (CGS), provides a single access point for earthquake strong-motion records and station metadata from the CGS California Strong-Motion Instrumentation Program (CSMIP), the USGS National Strong-Motion Project (NSMP), the USGS Advanced National Seismic System, and other affiliates. The CESMD has been continuously improving its webtools to facilitate the access of strong-motion data and metadata for use in post-earthquake response and for scientific and engineering research applications. The Center provides raw and processed strong-motion data via the Engineering Data Center (EDC) and the Virtual Data Center (VDC) web portals. This paper focuses on the strong-motion products provided by the EDC where more than 48,000 records with peak ground accelerations greater than 0.1% g from over 2400 earthquakes are currently hosted. and on the ongoing efforts to develop data access tools and applications. The new developments and ongoing efforts in the EDC include: 1) enhancements to the CESMD webservices to facilitate access to station metadata, earthquake information, and strong motion records 2) new features to the interactive map interface, improving the visualization and access to earthquake, station, and record information, 3) efforts to develop a new web application tool for data format conversion from a number of data formats, 4) efforts to unify varying waveform data formats into a consistent format, 5) ongoing efforts to compile seismic station site geology, measured or inferred Vs30 values, shear-wave profiles, NEHRP site class, and available structural instrument deployment schematics, and 6) a special studies pages for research topic-specific ground motion datasets that offer uniform processing of records from a variety of sources.

Conference Paper↗
USGS1996

Analysis of nutrients in the surface waters of the Georgia–Florida Coastal Plain study unit, 1970–91

During the early phase of the Georgia-Florida National Water Quality Assessment study, existing information on nutrients was compiled and analyzed in order to evaluate the nutrient concentrations within the 61,545 square mile study unit. Evaluation of the nutrient concentrations collected at surface-water sites between October 1, 1970, and September 30, 1991, utilized the environmental characteristics of land resource provinces, land use, and nonpoint and point-source discharges within the study unit. Long-term trends were investigated to determine the temporal distribution of nutrient concentrations. In order to determine a level of concern for nutrient concentrations, the U.S. Environmental Protection Agency (USEPA) guidelines were used: (1) for nitrate concentrations, the maximum contaminant level in public-drinking water supplies (10 mg/L); (2) for ammonia concentrations, the chronic exposure of aquatic organisms to un-ionized ammonia (2.1 mg/L); (3) for total-phosphorus concentrations, the recommended concentration in flowing water to discourage excessive growth of aquatic plants (0.1 mg/L); and (4) for kjeldahl concentrations, however, no guidelines were available. For sites within the 10 major river basins, median nutrient concentrations were enerally below USEPA guidelines, except for total-phosphorus concentrations where 45 percent of the medians exceeded the guideline. The only median ammonia concentration that exceeded the guideline occurred at the Swift Creek site (3.4 mg/L), in the Suwannee River basin, perhaps due to wastewater discharges. For all sites within the Withlacoochee, Aucilla, and St. Marys River basins, median concentrations of nitrate, ammonia, and total phosphorus were below the USEPA guidelines. Nutrient data at each monitoring site within each major basin were aggregated for comparisons of median nutrient concentrations among major basins. The Ochlockonee and Hillsborough River basins had the highest median nutrient concentrations, the Aucilla River basin had the lowest. Median concentrations of nitrate and ammonia among all major basins were below USEPA guidelines. The median total-phosphorus concentrations for the following river basins exceeded the USEPA guideline: Hillsborough, St. Johns, Suwannee, Ochlockonee, Satilla, Altamaha, and Ogeechee. Although nutrient concentrations within the study unit were low, long-term increasing trends were found in all four nutrients. All 18 study-unit wide nitrate trends had increasing slopes ranging from less than 0.01 to 0.07 (mg/L)/yr. The range in slope for the 13 ammonia trends was -0.03 to 0.01 (mg/L)/yr with 6 increasing trends in the northern part of the study unit. Of the 17 total-phosphorus trends found in the study unit, 10 were found at sites where the median concentration exceeded the USEPA guideline. At these 10 sites, 4 sites had increasing trends with slopes ranging from less than 0.01 to 0.07 (mg/L)/yr, 5 sites had decreasing trends with slopes ranging from -0.01 to -0.24 (mg/L)/yr, and one site showed a seasonal concentration trend. Median nutrient concentrations were significantly different among the four land resource provinces: Southern Piedmont, Southern Coastal Plain, Coastal Flatwoods, and Central Florida Ridge. As a result, nutrient concentrations among basins with similar nutrient inputs but located within different land resource provinces are not expected to be the same due to differences in the combination of factors such as soil permeability, runoff rates, and stream channel slopes. This concept is an important consideration in designing a surface-water quality network within the study area. For the most part, the Coastal Flatwoods showed the lowest median nutrient concentrations and the Southern Coastal Plain had the highest median nutrient concentrations. Lower median nitrate concentrations in surface-water basins were associated with the forest/wetland land-use category and higher median concentrations of nitrate and ammonia with the urban category when land-use percentages were classified into four land-use categories (agriculture, forest/wetland, mixed, and urban). These results were reasonable based on expected high nutrient inputs from urban areas and low inputs from forested and wetland areas. However, the lack of association between high nutrient concentrations and the agricultural land-use category was not expected since high nutrient inputs are generally needed for agriculture production.

Florida, Georgia↗
USGS2026

Upper Triassic igneous rocks of the southern Kenai Peninsula, Alaska – Prelude to Early Jurassic subduction along the western Wrangellia composite terrane margin

New U–Pb zircon geochronology identifies a latest Triassic (ca 214–201 Ma) igneous suite of tuff, hypabyssal dikes, and a pluton on the southern Kenai Peninsula, Alaska. The igneous suite was emplaced within Upper Triassic sedimentary rocks along the southern margin of Western Wrangellia, the western-most fragment of the Wrangellia composite terrane. The igneous rocks range from mafic (50.6% SiO 2 ) to felsic (78.3% SiO 2 ), characteristically have less than 1.55% K 2 O, and generally have low trace element abundances. The tonalitic and trondhjemitic magmas were largely sourced in mafic-rich lower crust and incompletely assimilated quartz and other mineral xenocrysts are common. Fractionation involving plagioclase and amphibole is indicated for some magmas and composite intrusions and igneous xenoliths indicate magma mixing was possible. Paleozoic and Precambrian inherited zircons and initial 87 Sr/ 86 Sr (0.704103–0.705609) and 143 Nd/ 144 Nd (0.512396–0.512777) ratios indicate that the Western Wrangellia crustal sources are heterogeneous and contain sialic components. The latest Triassic magmatism reflects processes that preceded Early Jurassic subduction along the Wrangellia composite terrane and Pacific Ocean plate boundary. These processes involved heating and melting of mantle lithosphere and lower crust as mantle instabilities accompanied the breaking of the plate boundary linkages. The Late Triassic transition to subduction along the Wrangellia composite terrane margin coincided with the transition to subduction cessation in the Late Triassic arcs of the western Intermontane terranes of Canada. The shift to subduction along the outboard Wrangellia composite terrane margin marks the beginning of the Pacific Ocean–Cordillera plate interactions that came to dominate the tectonic evolution of the northern Cordillera from the Early Jurassic to today.

Alaska↗
USGS2025

Regional characterization of coal resources in the U.S. Gulf Coast

There is increasing interest in extracting critical minerals (CM), including rare earth elements (REE), from coals in the United States to address the overreliance on imported REE. The U.S. Gulf Coast and the Williston basins are the two major lignite-bearing basins within the country. Recent REE and CM studies of the lignite in these basins have indicated that the coals may be a viable source material for REE and CM extraction. To evaluate in-place coal as a potential source of REE and CM, the coal resources need to be quantified. This study presents the results of a regional analysis of the U.S. Gulf Coast lignite and bituminous coal resources that might be available as potential sources of REE and CM. The resource analysis used kriging methods to develop isopleth maps of cumulative coal thickness throughout the region using data from 31,181 drill holes and other data points. The estimated total coal resource in the Gulf Coast is about 83 billion metric tons in the upper 90 m (~ 300 ft) of the subsurface. Texas accounted for 40 percent (32 billion metric tons) of the total resource, followed by Mississippi (24 %, 20 billion metric tons), Louisiana (14 %, 12 billion metric tons), Tennessee (10 %, 8.5 billion metric tons), and Arkansas (6 %, 5.1 billion metric tons). The remaining states each accounted for less than 5 percent of the total resource. Georgia had the smallest resource estimated at 7 million metric tons. Here we report the first known state-wide lignite resource estimates for Georgia, Kentucky (820 million metric tons), and Missouri (1,800 million metric tons). A comparison of the results of this study with those of previous Gulf Coast and Williston Basin resource studies is difficult because each study used different data sources, assessment methodologies, overburden depths, and qualifying coal thicknesses. Coal-power electric generation has sharply decreased in past decades and mining of these coals for CM and REE could provide additional co-products such as activated carbon and other uses such as fertilizer (soil enhancer).

EarthArXiv↗
USGS2026

Indirect mineral import reliance and provenance

Mineral commodity supply chain analyses rely on international trade data reported by individual countries as quantities of a mineral commodity form imported from (or exported to) a partner. However, export quantities frequently exceed a country’s domestic production, or occur when no production data are reported, suggesting that the trade partner is merely an intermediary in a transshipment. These discrepancies can result in misleading conclusions regarding supply chain vulnerabilities and dependencies. We present a two-stage methodology to reconcile gaps between reported material sources and actual producers. First, we construct trade networks for specific mineral forms, treating production as a type of import to distinguish producing nations from entrepôts. By tracing flows through these networks, we attribute a target country’s imports to original producers via both direct (in a single trade link) and indirect (transferring through intermediaries) pathways. Second, these production-attributed flows are incorporated into multi-stage supply chains to determine the upstream provenance of feedstock for domestic refining and processing. This approach provides a more representative picture of trade reliance. For example, while the United States (U.S.) Geological Survey reports no imports of unwrought antimony metal from Russia in 2022 (U.S. Geological Survey (2025). Mineral Commodity Summaries 2025. 10.3133/mcs2025), our analysis reveals that over 16% of U.S. imports can be traced back to Russian mining through intermediate processing in countries such as China, India, and Vietnam. Additionally, our analysis of the aluminum supply chain shows that while the U.S. is reported as 52% net import reliant on aluminum materials in 2022, it is 100% reliant on foreign bauxite, 7% of which arrived indirectly. This unreported reliance, which is predominantly tied to bauxite mined in Brazil (43%) and Jamaica (28%), highlights our methods ability to capture the supply chain’s dependence on foreign feedstock that may be missing in single-stage trade data.

Mineral Economics↗
USGS2026

Modeling legacy nitrogen transport under instantaneous, steady-state, and transient groundwater flow conditions

In hydrologic settings where groundwater discharge contributes substantially to surface waters, legacy nitrogen in groundwater can confound surface water nitrogen loads estimated exclusively from current terrestrial sources. Additionally, legacy nitrogen in groundwater can contribute to lagged responses to nitrogen management efforts. Some methods of estimating groundwater contributions to surface water nitrogen loads account for legacy nitrogen, while others do not. The resulting differences are rarely quantified. We used a numerical modeling framework to compare three methods of estimating time-varying annual groundwater nitrogen loads to surface water receptors on eastern Long Island, New York. The instantaneous load method used steady-state contributing areas and includes no temporal groundwater lag. The second method used numerical simulations of nitrogen loads under steady-state flow, which captures groundwater transport lags but omits the annual variability in transient hydrologic stresses. The third method numerically simulated both transient groundwater flow and nitrogen transport to explicitly capture the effects of legacy nitrogen in groundwater. Depending on antecedent nitrogen and hydrologic conditions, historical nitrogen loads estimated from the numerical simulations were sometimes similar (<10% difference) and other times substantially different (±100%) from the instantaneous load estimates. Additionally, simulated future surface water nitrogen loads responded asymptotically over several decades following reductions in terrestrial nitrogen sources, further highlighting the effect of groundwater transport lag times. The comparison of the three methods, quantification of historical interannual variability, and prediction of lagged responses to nitrogen source reductions provide important context for decision makers using estimated groundwater nitrogen loads to help evaluate nitrogen management efficacy.

New York↗
USGS2025

Potentiometric surface maps and groundwater-level hydrographs for confined aquifers of the New Jersey Coastal Plain, 2018

The U.S. Geological Survey, in cooperation with the New Jersey Department of Environmental Protection (NJDEP), prepared potentiometric surface maps for 10 confined aquifers of the New Jersey Coastal Plain physiographic province based on water-level measurements collected during late 2018 and early 2019 from 951 wells in New Jersey and parts of Pennsylvania and Delaware. Maps were prepared for the confined Cohansey aquifer, Rio Grande water-bearing zone, Atlantic City 800-foot sand, Piney Point aquifer, Vincentown aquifer, Wenonah-Mount Laurel aquifer, Englishtown aquifer system, and the upper, middle, and lower aquifers of the Potomac-Raritan-Magothy aquifer system. Potentiometric surface maps indicate regional cones of depression in the following aquifers and the counties in which they are centered: Atlantic City 800-foot sand in Atlantic County, the Piney Point aquifer in Cumberland County, the Wenonah-Mount Laurel aquifer and Englishtown aquifer system in Monmouth and Ocean Counties, the Wenonah-Mount Laurel aquifer in Camden and Gloucester County, the upper aquifer of the Potomac-Raritan-Magothy aquifer system in Ocean County, and the upper, middle, and lower aquifers of the Potomac-Raritan-Magothy aquifer system in Camden County. Cones of depression with smaller areal extents were in the confined Cohansey aquifer, the Rio Grande water-bearing zone, the Atlantic City 800-foot sand centered in Cape May County, the Piney Point aquifer centered in Ocean County, the Wenonah-Mount Laurel aquifer in Salem and Burlington Counties, the Englishtown aquifer system in Camden County, the upper aquifer of the Potomac-Raritan-Magothy aquifer system in Monmouth County, and the middle aquifer of the Potomac-Raritan-Magothy aquifer system in Monmouth, Ocean, and Salem Counties. No cone of depression was interpreted in the Vincentown aquifer. Long-term hydrographs are presented for 75 wells spanning each of the 10 confined aquifers, and contain a mix of discrete water-level measurements and daily mean water levels based on continuously recorded 15-minute data. Changes of water levels during 2014–19, as indicated by the hydrographs, were compared with those of previous periods to assess any departures from historical data. During 2014–19, water levels were stable and fluctuated within similar ranges as previous periods in the following aquifers and locations: all wells in the confined Cohansey aquifer, the Rio Grande water-bearing zone, the Vincentown aquifer, the Englishtown aquifer system, the Piney Point aquifer wells in Burlington and Ocean Counties, six of eight wells in the Wenonah-Mount Laurel aquifer, all wells in the upper and lower aquifers of the Potomac-Raritan-Magothy aquifer system outside NJDEP Critical Areas, and all wells in the middle aquifer of the Potomac-Raritan-Magothy aquifer system except those within NJDEP Critical Area II. Increasing water levels in 2014–19, ongoing since historical periods, were indicated in the following aquifers and locations: Atlantic County wells in the Piney Point aquifer, all wells in the upper and lower aquifers of the Potomac-Raritan-Magothy aquifer system outside NJDEP Critical Areas, and all wells in the middle aquifer of the Potomac-Raritan-Magothy aquifer system within NJDEP Critical Area II. Water levels in the Atlantic City 800-foot sand also increased during 2014–19 in wells in Atlantic County and northern Cape May County closer to the center of the cone of depression in that aquifer, which is a response unique to this period and absent from previous periods. During 2013–19, continued decreasing water levels, ongoing since previous periods, were indicated by hydrographs of Atlantic City 800-foot sand wells in southern Cape May County, Piney Point aquifer wells in Cumberland County where the regional cone of depression is located, and two wells in the Wenonah-Mount Laurel aquifer—070478, which in 2014–19 departed from previous periods, and 330020, which continued a gradual decrease throughout its period of record.

New Jersey↗
USGS2025

Do bighorn sheep use desert tortoise burrow spoil piles as mineral licks in southern California?

In a previous study, bighorn sheep ( Ovis canadensis ) were periodically photographed by trail cameras when they visited desert tortoise ( Gopherus agassizii ) burrows in the San Bernardino Mountains near Palm Springs, California, USA. The authors suggested that bighorn sheep may utilize the excavated spoil material from burrows as a mineral lick. To test that supposition, we collected soil samples from the spoil piles and undisturbed surface soil up to 1 m away at two desert tortoise study sites with bighorn sheep at the aforementioned site and another nearby site. We hypothesized that sodium concentrations would be greater in spoil than in nearby surface soil samples. Common major cations were quantified in soil samples using inductively coupled plasma-mass spectrometry following microwave-assisted digestion to solubilize the fraction of elements that could become biologically or environmentally available. Mean sodium concentrations were greater in burrow spoil samples than in undisturbed soil samples at both sites, likely due to aqueous leaching of sodium from surface to subsurface soil horizons. Although we did not collect evidence of geophagy by bighorn sheep, the behavior is well known. However, trail camera photographs of bighorn sheep with their head next to the entrance of tortoise burrows suggests that they may periodically utilize sodium and other elements brought to the surface by desert tortoises that are otherwise unavailable to bighorn sheep on the surface.

California↗