USGS Science⌕ Search

SEARCH · USGS Science

Results for “Oil and Gas Journal”

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,658 records · Page 62Linked to original sources
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

Software to support remote sensing of river discharge based on critical flow theory

Water resource management requires accurate observations of streamflow but standard field methods for measuring river discharge ( Q ) are costly and can be hazardous for equipment and personnel. Remote sensing has become a viable alternative, but many image-based techniques require field data for calibration and depth and velocity can seldom be mapped with a single sensor. A new approach based on critical flow theory, in contrast, allows both of these attributes to be inferred from readily available image data. This technique only pertains to sites with standing waves, called undular hydraulic jumps (UHJs), but a recent investigation demonstrated its potential to provide accurate discharge estimates. This paper introduces software designed to facilitate Inferring Q from UHJs Identified in River Images (InQUIRI). The package includes modules for retrieving data from image servers, making the measurements of wavelength and width required to calculate discharge, inferring a representative wavelength from a profile digitized along a wave train, combining multiple estimates to obtain an ensemble median discharge, and assessing accuracy via comparison to gage records from the U.S. Geological Survey. By making these steps easier to implement, InQUIRI enables users to apply the workflow to a variety of UHJ-containing images. Accumulating more case studies, some successful and others less so, would help constrain the range of applicability of the critical flow approach and foster development of refined guidelines for selecting and measuring waves. The software described herein could play an important role in promoting informed use of this new technique for non-contact streamflow measurement.

Arizona, Colorado, New Mexico, Utah↗
USGS2026

Influence of new residential construction varying in housing density on bird species, human tolerance guilds, and communities

Human population growth and changing settlement patterns fuel the development of urban fringe lands worldwide, with implications for biodiversity. We conducted a 12-year study of birds in the fast-developing urban fringe lands of the central Puget Sound region, Washington, USA, to examine the effect of development configuration on birds. We hypothesized that lower-intensity conservation developments, compared to higher-intensity planned community developments, would benefit the overall bird community, as well as native forest birds (avoiders of human development) and avian generalist species (adapters to human development), but that higher-intensity planned community developments would benefit synanthropic species (exploiters of human development). We fit single-species and multi-species occupancy models to test these hypotheses. Consistent with our hypotheses, we found that a greater proportion of the overall community, avoiders, and adapters occupied lower-intensity conservation developments compared to higher-intensity planned community developments. However, we did not detect an effect of development type on the exploiter guild, and we found that species in the exploiter guild are variable in their response to the configuration of suburban developments. We also hypothesized that human tolerance guilds would be a useful predictor of individual species responses to development type. This hypothesis was somewhat supported: we found that, for avoiders, 87% of species in the guild had the same response to development type as the overall guild; for adapters, 63% had the same response as the overall guild, and for exploiters, only 44% had the same response as the overall guild. Our results indicate that the configuration of suburban developments can have a meaningful impact on bird communities, particularly on those species that are most sensitive to any level of development.

Washington↗
USGS2026

Estimation of parental abundance using hierarchical Bayesian modeling with data augmentation

Pedigree-based estimation methods leverage the fact that each offspring in a cohort is genotypically “marked” by its parents and represent a recent and promising toolset for estimating population dynamics. This includes pedigree accumulation estimators that model the “accumulation” of inferred unique parents within a given cohort to estimate parental abundance. Unlike close-kin mark-recapture approaches, which rely on intercohort comparisons, pedigree accumulation modeling can be completed solely using intracohort samples. This is particularly advantageous for semelparous species, where intercohort pairs are impossible and adult life stages can be difficult to sample without affecting their likelihood of successfully reproducing. Previous work has evaluated a range of estimators for such datasets, concluding that the non-parametric Chao estimator provides the most accurate and precise estimates for feasible levels of sampling effort. We used simulated data to evaluate an alternative estimator based on hierarchical modeling and data augmentation in a Bayesian framework. Results indicate that estimates from the hierarchical Bayesian estimator had comparable accuracy and better precision than both the previously tested Chao1 estimator and the improved iChao formulation across a range of sample sizes and sex ratios. Furthermore, the Bayesian estimator was far more robust to simulated errors in pedigree reconstruction, especially the presence of false negatives. Hierarchical Bayesian pedigree accumulation models can also provide additional insight into underlying reproductive ecology through their use of an explicit observation process, allowing for the incorporation or estimation of species- and population-specific reproductive dynamics. More broadly, the parametric nature of these models offers opportunities to efficiently pool information among datasets as well as to propagate uncertainty within more complex models.

Ecology and Evolution↗
USGS2024

Changes in soil erosion caused by wildfire: A conceptual biogeographic model

Soil erosion rates after wildfire are strongly controlled by intrinsic properties such as topography, weather, climate, soil, and vegetation. These landscape and hydroclimatic properties are important in determining post-fire erosion rates; however, their influence on post-fire erosion and their interaction with the intensity of a wildfire remains uncertain. A key limitation in resolving this uncertainty is the lack of conceptual models and frameworks for organising data related to the geomorphic sensitivity of landscapes to wildfire. Our aim is to develop a framework for consolidating understanding of post-fire erosion in the context of hydroclimatic conditions which contribute to system states, for example soil and vegetation properties, and wildfire regime. The framework is developed around a simple conceptual model where the change in erosion due to wildfire is a product of change in runoff generation and sediment supply, which is strongly related to landscape net primary productivity (NPP). We hypothesised that geomorphic sensitivity to wildfire should vary as a unimodal humped relationship across a gradient of NPP, peaking at an intermediate level. To develop this framework and to test the hypothesis, we first review intrinsic soil and vegetation properties related to the supply and transport of sediment from burned and unburned hillslopes. Net primary productivity is systematically related to these intrinsic properties because it integrates many processes involved in soil and vegetation development. Empirical data indicate a trend in the change in surface runoff generation with NPP after wildfire, peaking at an NPP of approximately 15 Mg C ha −1 y −1 . A simple model of fuel availability and soil heating are correlated with a similar “humped” trend in sediment supply. These results are consistent with our conceptual model, which indicates that sediment supply and runoff contribute towards a distinct peak in wildfire effects on erosion at an intermediate level of NPP. We propose that landscapes of intermediate NPP typically have the highest quantity of fuel available to burn, which cause large changes to the soil surface properties. Landscapes at intermediate NPP also tend to produce intrinsic soil and vegetation properties that promote erosion after wildfire. The interplay between these short and long-term landscape characteristics is strongest at intermediate levels of NPP. Our proposed biogeographic model of geomorphic sensitivity to wildfire was supported by erosion data from burned hillslope and zero-order catchments studies from a range fire-prone landscapes in Australia and North America. Our proposed conceptual model will help identify areas most vulnerable to post-fire erosion changes.

Geomorphology↗
USGS1983

Base flow and ground water in upper Sweetwater Valley, Tennessee

The upper Sweetwater Valley area has a flow system with complex interaction between surface and ground water. A water budget study indicated that during dry years approximately three-fourths of the annual flow to Sweetwater Creek may be derived from ground-water sources. Hydrograph analysis showed seasonal variation of recharge to the ground-water flow system. Streamflow records were analyzed to estimate the frequency of low flow of Sweetwater Creek at river mile 16.7, and indicated the lowest average flow for 1 day in 20 years to be about 5.1 cubic feet per second. Two periods of base-flow measurements of Sweetwater Creek identified channel reaches with significant gains and losses of streamflow. Base flow measurements also showed interbasin transfer of water among sub-basins of the valley. Major flow surpluses were associated with areas in which the majority of flow originated at a spring. Topographically low areas adjacent to the main stem of Sweetwater Creek generally have surplus flow. Topographically higher areas generally have deficient surface outflow unless significant spring flow occurs in the basin. Ground-water recharge occurs by water draining into sinkholes, faults, and fractures. Ground-water flow is regionally diffused across formation strikes from the topographically low areas unless the water is exposed to highly permeable formations or impervious formations. Ground water infiltrates the highly permeable formations and flows along strike. Ground water encountering impervious formations may discharge at small springs at the contact, or may reroute along the contact if the upgradient rock is sufficiently permeable or has well developed secondary porosity. Groundwater discharges to streams at innumerable seeps and at a few large springs. Areas of ground-water flow up-gradient of large springs are hypothesized as likely areas of significant ground-water reservoirs.

Tennessee↗
USGS2026

King Rails in southeastern Arkansas are majority resident

King Rails ( Rallus elegans ), a declining species of secretive marsh bird, have a mixed migratory strategy, with some individuals migrating and some being year-round residents. Determining the distribution of migratory versus resident King Rails is critical to guiding habitat management and assessing risks to populations. We examined the migratory status of adult King Rails at the Choctaw West Wildlife Management Area in southeastern Arkansas, the only known site in Arkansas with multiple breeding pairs. Of the 23 adult King Rails we tagged in 2022 and 2023, two birds migrated away from the local area, one in each year. One migrated 145 km, and the other 450 km. Our findings suggest that the line between the migratory and resident range of Arkansas King Rails could be north of our site.

Arkansas↗
USGS2026

Tringa flavipes (Lesser Yellowlegs) from separate breeding sites subdivides the Prairie Pothole Region in space and time during southbound migration

Some staging regions support multiple groups of the same migratory species, each of which may use the region differently. Characterizing the ways, in which separate groups use such regions can therefore help to identify vulnerabilities during this sensitive period of the annual cycle. The Prairie Pothole Region (PPR) is a massive wetland complex in the northern Great Plains of North America used by ∼11 million shorebirds during migration. The PPR has been heavily modified by agriculture and is experiencing varied effects of global climate change, threatening the health of the shorebirds that rely on it. Here, we used 6 seasons of southbound tracking data of Tringa flavipes (Lesser Yellowlegs)—a long-distance migratory shorebird species with an estimated population decline of 63% over the last 4 decades—from 9 sites across their breeding range to explore differences in migratory behavior within this important staging region. We found that 75% of tracked individuals used the region during southbound migration, and T. flavipes from different breeding sites detoured 110–875 km from their most direct migratory route to access the PPR. Individuals that arrived later stayed longer and made more stops within the region than those that arrived early. Individuals originating from different breeding sites also displayed spatial and temporal segregation within the region: T. flavipes from southwest and central Alaska relied heavily on the northwestern PPR, while those from Canada used the central and southeastern portions of the PPR. Finally, timing of use varied among groups, but the southeastern PPR became increasingly important over the course of the southbound migratory window, as other wetlands likely dried out. Our study highlights the portions of the PPR of critical importance to migrating T. flavipes and the diversity of ways, in which different groups from within the same species can use a single staging region.

Alaska, Manitoba, Northwest Territories, Ontario, ↗
USGS2024

Earthquake scenario development in conjunction with the 2023 USGS National Seismic Hazard Model

We present earthquake scenarios developed to accompany the release of the 2023 update to the US Geological Survey National Seismic Hazard Model (NSHM). Scenarios can serve a range of local and regional needs, from developing proactive-targeted mitigation strategies for minimizing impending risk to aiding emergency management planning. These deterministic scenarios can also be used to communicate seismic hazard and risk to audiences who are not well versed in methods, such as probabilistic seismic hazard analyses. Specifically, we discuss the scenarios developed, challenges, and lessons learned in the development process, and how this work aided the development of the 2023 NSHM itself. In total, 28 scenarios were developed for Hawaii, Utah, Alaska, and Virginia considering the 2023 NSHM science, past scenario efforts, and input from local experts and stakeholders. Finally, we investigate how NSHM modeling decisions can change estimated impacts to Utah and Hawaii in more detail showing, for example, that a shallower dip of the Wasatch fault under Salt Lake City can increase predicted ground-motion intensities and therefore estimated losses and deaths.

Earthquake Spectra↗
USGS2026

Bayesian belief network model to predict human-wildlife conflict in protected areas

Human-wildlife conflict (HWC) poses a pervasive global challenge, affecting livelihoods and threatening biodiversity. To better anticipate and mitigate HWC risk, we developed a large-scale predictive model using a Bayesian Belief Network (BBN). We surveyed 1,011 park rangers across 135 terrestrial protected areas in three Andean countries, documenting recent HWC incidents involving wildlife persecution or killing, livestock depredation, crop damage, or threats to human safety and property. We identified key drivers of HWC risk, including governance, wildlife acceptance, participation, and habitat quality. A sensitivity analysis revealed that enhancing governance and improving wildlife acceptance could reduce HWC risk by > 85%. The BBN model demonstrated scalability, effectively identifying strategies to reduce HWC risk at multiple scales, from individual protected areas to national networks. Our findings highlight the importance of strengthening governance, increasing wildlife acceptance, and enhancing community participation in conservation efforts. BBNs provide a flexible, cost-effective, and data-driven tool to guide protected areas and wildlife managers in monitoring, anticipating, and making informed decisions to mitigate conflict and promote coexistence.

Scientific Reports↗
USGS2024

Post-wildfire debris flows

Post-wildfire debris flows pose severe hazards to communities and infrastructure near and within recently burned mountainous terrain. Intense heat of wildfires changes the runoff characteristics of a watershed by combusting the vegetative canopy, litter, and duff, introducing ash into the soil and creating water repellant soils. Following wildfire, rainfall on bare ground is less able to infiltrate into the fire-altered soils and overland flow is less impeded by vegetation. Rainfall runoff in recently burned areas can erode hillslopes owing to the removal of soil binding organic matter near the soil surface by fire. In channels, loose, dry-ravel deposits composed of sand and gravel are readily entrained by concentrated runoff in channels. Entrainment of soil on hillslopes and in channels bulks up the sediment concentration of the rainfall runoff to generate debris flows capable of transporting boulders and large woody debris. Post-wildfire debris flows can be triggered by rainfall conditions that would typically produce little runoff during unburned conditions. The primary rainfall trigger for post-wildfire debris flows is high intensity rainfall during short duration convective rainstorms or periods of high rainfall intensity embedded within a long-duration frontal storm. Numerous observations of debris flows triggered by storms lasting less than an hour following periods of little to no rainfall indicate that antecedent rainfall is not a requirement for initiation of post-wildfire debris flows. Post-wildfire debris-flow hazard assessment entails estimating probability and magnitude of debris flows in the burned area, estimating debris-flow runout and intensity, and defining rainfall intensity-duration thresholds for debris-flow initiation. In the United States, probability and magnitude is estimated using empirically derived models largely based on data collected in southern California. The models provide maps to identify watersheds and drainage paths where post-wildfire hazards are most pronounced. Rainfall intensity-duration thresholds can be incorporated into flood hazard forecasting tools. Currently, work is underway to identify how to best implement debris-flow runout models in burned areas with efficiency and accuracy. Post-wildfire debris flows have been a long-recognized process in the Transverse Ranges of southern California; however, climate change is driving more frequent wildfires to burn more mountainous terrain throughout the western United States and worldwide. As a result, post-wildfire debris flows are becoming a more common threat in areas where they were once infrequent. As the threat of post-wildfire debris flow expands into new areas, evaluating the hazard becomes challenging because the degree to which wildfire increases debris-flow susceptibility varies from region to region. This chapter summarizes the knowledge to date for evaluating post-wildfire debris-flow susceptibility and hazard assessment. We summarize the characteristics of wildfire burn severity, topography, underlying soil and geology, and rainfall conditions that contribute to making a watershed most likely to produce post-wildfire debris flows. Methods for hazard assessment in the United States and other countries are summarized. We highlight knowledge gaps for how post-wildfire debris-flow susceptibility varies throughout the western United States and worldwide and identify research needs to improve hazard assessment methods in different geographies.

Book chapter↗
USGS2025

Central Valley Hydrologic Model version 2 (CVHM2): Decision support tool for groundwater and land subsidence management

The San Joaquin Valley (SJV) of California is one of the world’s most productive agricultural regions. Reliance on groundwater has led to some of the greatest rates of human-induced land subsidence in the world in the 20th century, as well as more recently. The United States Geological Survey (USGS) has recently developed an integrated surface–subsurface hydrologic model, the Central Valley Hydrologic Model 2 (CVHM2), that represents the major components of the hydrologic system of California’s Central Valley. In this study, CVHM2 was applied as a decision support tool while simulating various management strategies to mitigate the land subsidence caused by the extraction of groundwater. CVHM2 was extended through to 2073 and applied to simulate management scenarios in terms of three primary drivers and their impact on subsidence along the Delta–Mendota Canal (DMC), a critical piece of infrastructure in the western SJV. The drivers considered were agricultural water demands, managed aquifer recharge (MAR), and changes in future climate. The results show that future subsidence is most sensitive to water demands, second most sensitive to future changes in climate, and relatively insensitive to MAR when it is applied as a surface application in the western SJV. However, we demonstrate via proof-of-concept scenarios that the MAR is capable of arresting subsidence when implemented via injection below the Corcoran Clay Member of the Tulare Formation instead of as a surface application. We also examine the uncertainty that is the result of climate variability and how to use the tool to identify the most appropriate strategies to constrain future subsidence to acceptable levels.

California↗
USGS2026

Decoding the oxidative digestion mechanism for polystyrene nanoplastic detection in the Great Lakes using a customizable Raman spectral processing algorithm

Despite the concerns surging around nanoplastics (NPs) regarding their prevalence and bioavailability in freshwater systems, robust detection of NPs in complex environmental matrices is hindered by the lack of standardized sample pretreatment and a mechanistic understanding of oxidative digestion. Here, we systematically investigate the interaction between hydrogen peroxide (H 2 O 2 ) and polystyrene (PS) NPs during digestion in deionized (DI) water and four environmental matrices from in and around the Great Lakes Basin. To facilitate high-throughput analysis, we develop Pre_peak, a customizable Raman spectral processing algorithm that achieves >99% accuracy for both NP identification and interference rejection, allowing reliable NP quantification via pixel counting and systematic decoding of the oxidative digestion mechanisms. In DI water, varying H 2 O 2 doses from 0 to 30% has negligible effects on the recovery and Raman signal intensity of PS NPs over 24 hours of digestion. However, morphological changes and aggregation of PS NPs are observed when the H 2 O 2 dose exceeds 20%. Prolonged digestion further leads to progressive NP loss. In natural waters, the optimal dosage and digestion duration depend on matrix characteristics, including dissolved organic matter (DOM) and ion composition. This study provides mechanistic insights into NP–oxidant interactions and underscores the need for matrix-tailored digestion protocols to advance standardized NP detection in freshwater environments.

Great Lakes↗
USGS2023

Don’t Let Negatives Hold You Back: Accounting for Underlying Physics and Natural Distributions of Hydrothermal Systems When Selecting Negative Training Sites Leads to Better Machine Learning Predictions

Selecting negative training sites is an important challenge to resolve when utilizing machine learning (ML) for predicting hydrothermal resource favorability because ideal models would discriminate between hydrothermal systems (positives) and all types of locations without hydrothermal systems (negatives). The Nevada Machine Learning project (NVML) fit an artificial neural network to identify areas favorable for hydrothermal systems by selecting 62 negative sites where the research team had confidence that no hydrothermal resource exists. Herein, we compare the implications of the expert selection of negatives (i.e., the NVML strategy) with a random sample strategy, where it is assumed that areas outside the favorable structural ellipses defined by NVML are negative. Because hydrothermal systems are sparse, it is highly probable that, in the absence of a favorable geological structure, hydrothermal favorability is low. We compare three training strategies: 1) the positive and negative labeled examples from NVML; 2) the positive examples from NVML with randomly selected negatives in equal frequency as NVML; and 3) the positive examples from NVML with randomly selected negatives reflecting the expected natural distribution of hydrothermal systems relative to the total area. We apply these training strategies to the NVML feature data (input data) using two ML algorithms (XGBoost and logistic regression) to create six favorability maps for hydrothermal resources. When accounting for the expected natural distribution of hydrothermal systems, we find that XGBoost performs better than the NVML neural network and its negatives. Model validation was less reliable using F1 scores, a common performance metric, than comparing probability estimates at known positives, likely because of the extreme natural class imbalance and the lack of negatively labeled sites. This work demonstrates that expert selection of negatives for training in NVML likely imparted modeling bias. Accounting for the sparsity of hydrothermal systems and all the types of locations without hydrothermal systems allows us to create better models for predicting hydrothermal resource favorability.

Geothermal Resources Council Transactions↗
USGS2024

Uncertainty and spatial correlation in station measurements for mb magnitude estimation

The body‐wave magnitude (⁠⁠) is a long‐standing network‐averaged, amplitude‐based magnitude used to estimate the magnitude of seismic sources from teleseismic observations. The U.S. Geological Survey National Earthquake Information Center (NEIC) relies on in its global real‐time earthquake monitoring mission. Although waveform modeling‐based moment magnitudes are the modern standard to characterize earthquake size, is important because (1) in many cases, waveform modeling is not possible (e.g., low signal‐to‐noise events), (2) is applicable over a broad range of magnitudes, ∼M 4–7, and (3) there is a many decades‐long history of estimating magnitudes. We use the NEIC Preliminary Determination of Epicenters earthquake catalog to investigate the uncertainty in NEIC station measurements. We show that measurements are spatially correlated, which can bias event ⁠, and we describe an empirical relation between this spatial correlation and station‐to‐station distance. We further describe an approach to mitigate bias from the spatial correlation. Accounting for the spatial covariance of observations can change the event from −0.15 to 0.07 units (10th to 90th percentile) for smaller events (⁠⁠). These smaller events have the largest standard deviations ranging from 0.05 to 0.15 units (10th to 90th percentile).

The Seismic Record↗
USGS2024

Seismicity and anisotropic imaging reveal an active detachment beneath the northern Alaska Range foothills

North of the Denali Fault, the collision between the Yakutat block with North America is accommodated by a fold-thrust belt giving rise to the northern Alaska Range foothills. At the western end, the Kantishna Hills anticline hosts prominent microseismicity and surface deformation, interpreted as active folding of the Kantishna Hills anticline above a midcrustal detachment. We test for this detachment by using anisotropy-aware receiver functions to image fabric contrasts within the crust in context with seismicity. Seismic stations near the crest of the Kantishna Hills anticline and near its southern flank show a single strong contrast in dipping fabric at depths of 12–13 km near microseismicity clustering depths, consistent with a detachment plane beneath the fold. A minimum b -value at 10–13 km depth is consistent with seismicity on the detachment, compatible with the imaged anisotropic contrast, while off-fault seismicity is shallower and deeper with smaller magnitudes. South-dipping imbricate thrusts in schist characterize the northern Alaska Range foothills structure. This supports our interpretation of the observed anisotropy as reflecting SSW–SSE-dipping foliation above a detachment at ∼10–13 km depth that exploits existing crustal weaknesses along subtle fabric contrasts observed in the seismically quiescent region north of the actively deforming belt.

Alaska↗
USGS1991

Reconnaissance of the occurrence of agricultural chemicals in ground water in Haywood, Lake, Obion and Shelby Counties, Tennessee

Data on the occurrence of agricultural chemicals in ground wafer in Tennessee are sparse. The surficial alluvial aquifer is an important source of domestic water supply in West Tennessee, and potentially is subject fo contamination from the application of agricultural chemicals in the area. Nineteen shallow wells completed in the alluvial aquifers in areas of high density agricultural use were sampled in the winter and again in the summer of 1988 to ascertain the occurrence of agricultural chemical in ground water. Although no triazine herbicides or organophosphorus insecticides were detected in any of the wells sampled, elevated nitrite plus nitrate (as nitrogen) concentrations were detected. Results from the winter sampling period indicate a range of nitrite plus nitrate (as nitrogen) concentrations of less than 0.1 to 7.8 milligrams per liter with a median concentration of 2.6 milligrams per liter. Results from the summer sampling period indicate a range of nitrite plus nitrate (as nitrogen) concentrations of less than 0.1 to 8.9 milligrams per liter, median, 2.5 milligrams per liter. The highest concentrations occurred in the shallowest wells, and, in one instance, in a shallow well near a heavily irrigated field.

Tennessee↗
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↗