USGS Science⌕ Search

SEARCH · USGS Science

Results for “Journal of Wildlife Management”

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,661 records · Page 58Linked to original sources
USGS2025

Geochemistry and radiogenic isotopes constrain the mantle source region of the Mountain Pass Intrusive Suite, California

The Mountain Pass carbonatite stock is the largest rare earth element (REE) deposit and only active REE mine in the United States. The carbonatite intrusion and spatially associated alkaline silicate intrusions constitute the Mountain Pass Intrusive Suite, which is located within the Mojave Province in California. Both the carbonatite and the alkaline silicate rocks are enriched in large ion lithophile elements and light REEs, and less enriched to depleted in high field strength elements, indicating the mantle source region was metasomatically enriched in incompatible trace elements. The cause of this metasomatic mantle enrichment and the genetic relationship between the carbonatite and the alkaline silicate stocks are poorly understood. In this study, major and trace element geochemical data and isotopic (Rb-Sr, Sm-Nd, and Lu-Hf) data are presented to constrain genesis of the Mountain Pass Intrusive Suite, from mantle source region to the intrusion of the stocks. Our geochemical data are consistent with derivation of the alkaline silicate and carbonatite melts through partial melting from a shared mantle source region rather than through liquid immiscibility or fractional crystallization and separation of a carbothermal fluid. Although the Rb-Sr isotopic system in the Mountain Pass Intrusive Suite is disturbed at the whole-rock scale, the isotopic systems for whole-rock Sm-Nd (εNd i = ‐2.2 ± 0.8) and zircon Lu-Hf (εHf i = 0.1 ± 1.1) are robust and support mantle derivation of the magmas. Geochemical modeling using experimentally derived partition coefficients was used to identify possible causes of enrichment in incompatible elements through metasomatism in the mantle source region. Modeling of metasomatism by melts derived by partial melting of deeply subducted carbonated sediments approximates observed Mountain Pass Intrusive Suite trace element chemistry. Scattered εHf i in inherited zircon (2.8 ± 2.6) is consistent with derivation from an arc-related environment with substantial crustal contamination. Paleotectonic studies in the Mojave Province indicate that regional subduction preceded emplacement of the Mountain Pass Intrusive Suite by ∼300 Ma. Melting of the Mountain Pass source region may have been caused by post-collisional thermal relaxation and extension.

California↗
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↗
USGS2024

New 10Be-26Al isochron burial dating informs the Pliocene and Pleistocene evolution of the lower Colorado River, southwestern United States

Four new 10 Be- 26 Al isochron burial ages ranging from 4.4 to 2 Ma on ancestral Colorado River deposits in the lower Colorado River corridor (LCRC) help constrain the river’s evolution during the Pliocene and early Pleistocene. They help fill a gap between previous work that focused on older and younger deposits: Older dated deposits include the 5 Ma Bouse Formation, which records the integration of the Colorado River through a series of preexisting basins to the Gulf of California and the ca. 4.5 - 3.5 Ma Bullhead Alluvium, a 200 to 300 m thick aggregational package that immediately followed integration. The much younger, 100 - 70 ka, Chemehuevi Formation is another major aggradation package mapped throughout the LCRC. The new burial ages on the facies of Santa Fe Railway (4.37 ± 0.71 Ma), boulder conglomerate of Bat Cave Wash (2.12 ± 0.26 and 2.05 ± 0.31 Ma), and the Palo Verde alluvium (3.03 ± 0.26 Ma) partially fill in a 3.5 M.y. gap between the deposition of the Bullhead Alluvium and the Chemeheuvi Formation and document the timescales over which the Colorado River was able to remove the Bullhead aggradational package and initiate newer and smaller aggradational pulses.

Arizona, California, Nevada↗
USGS2026

Status and understanding of groundwater quality in the San Joaquin Valley Kern County subbasin domestic-supply aquifer study unit, 2022—California GAMA Priority Basin Project

The quality of water accessed by domestic wells (here referred to as domestic groundwater resources) in the San Joaquin Valley Kern County subbasin (basin number 5-022.14) was assessed as part of the California Groundwater Ambient Monitoring and Assessment (GAMA) Program Priority Basin Project (GAMA-PBP), in cooperation with the California State Water Resources Control Board. Kern County is at the southern end of the San Joaquin Valley in California, and about 30,000 residents are estimated to use privately owned domestic wells for drinking water. Domestic wells typically draw from shallower parts of the aquifer system than public-supply wells and can be more vulnerable to effects from surface activities. Kern County is host to a highly productive agricultural industry, with Bakersfield as the main urban center. The Kern River runs through Bakersfield from the southern Sierra Nevada and intersects the Kern Water Bank, one of the largest groundwater banking operations in California, at the Kern River Intertie. The section of the Kern River running through the Kern Water Bank is dry most years. Kern County also encompasses some of the most productive oil and gas basins in California, with extensive underground and surface disposal of oil-field wastewater. This study was based on data collected from 33 sites sampled by the U.S. Geological Survey for the GAMA-PBP in 2022. To provide context for the water quality assessment, measured concentrations were compared to regulatory and non-regulatory health-based and aesthetic benchmarks. A grid-based method was used to estimate the proportions of the groundwater resources used for domestic-supply wells that have water-quality constituents below (low relative concentration), approaching (moderate relative concentration), or above (high relative concentration) benchmark concentrations. At least one measured constituent with a regulatory benchmark was categorized as having a high relative concentration in 72 percent of the aquifer area used for domestic groundwater resources. Inorganic constituents were detected at high concentrations in 45 percent of the domestic groundwater resources, and the constituents detected above regulatory benchmarks were arsenic, nitrate, and uranium. At least one organic constituent was detected at high concentrations in 41 percent of the domestic groundwater resources, and the constituents exceeding regulatory benchmarks were the fumigants 1,2,3-trichloropropane (1,2,3-TCP), 1,2-dibromo-3-chloropropane (dibromochloropropane [DBCP]), 1,2-dibromoethane (EDB), and the per-and polyfluoroalkyl substance (PFAS) perfluorooctanesulfonate. The disinfection by-product chloroform, the fumigant 1,2-dichloropropane, the herbicides atrazine and hexazinone, and the herbicide degradates 2-chloro-6-ethylamino-4-amino-s-triazine, 2-chloro-4,6-diamino-s-triazine, 4-hydroxychlorothalonil, and metolachlor sulfonic acid were detected in more than 10 percent of domestic groundwater resources, but concentrations did not exceed regulatory benchmarks. Land use, groundwater age (fraction of modern water and mean age), and geochemical environment (oxic or anoxic conditions, pH, alkalinity) were associated with the distribution of high relative concentrations of inorganic and organic constituents. Young, oxygenated water is recharged along the Kern River and adjacent recharge ponds, or as irrigation water in the agricultural areas. High concentrations of nitrate and volatile organic compounds occurred in the oxic water in urban and agricultural areas. The fumigants 1,2,3-TCP, DBCP, and EDB were reported throughout the agricultural areas, whereas chloroform, tetrachloroethene, and PFAS were associated with urban land use. High uranium concentrations were associated with young, modern groundwater in agricultural areas with low pH and high bicarbonate. Total dissolved solids increased with distance from the Kern River, as the contributions of fresh, oxic water decreased. High concentrations of arsenic were present in older anoxic or alkaline groundwater away from areas of recharge. Overall, groundwater age, redox conditions, and the source of recharge as a result of different land uses contribute to large aquifer-scale portions of domestic groundwater resources that exceed health-based benchmarks for nitrate, uranium, and fumigant concentrations.

California↗
USGS2025

Petrogenesis and mineralization potential of the Bradley Peak komatiitic basalts, Wyoming Province

The ca. 2.72 Ga Bradley Peak greenstone terrane in the Wyoming Province contains spinifex- and cumulate-textured komatiitic volcanic rocks that may host Ni-Cu-PGE sulfide mineralization. New whole rock geochemistry classifies these rocks as Al-undepleted komatiitic basalts derived from a parental melt with ~17 wt. % MgO. Isotopic data (Sm-Nd, Re-Os) and REE profiles suggest a plume source and possible assimilation of juvenile footwall volcanic rocks. Calculated sulfur concentrations at sulfide saturation and low-detection PGE geochemistry indicate these magmas did not saturate a sulfide liquid resulting in sulfide mineralization. Although the Bradley Peak komatiitic basalt flows are unlikely to host magmatic sulfide mineralization, they are likely the source of Au for regional orogenic gold deposits.

Wyoming↗
USGS2026

Patterns and drivers of cliff erosion in Big Sur, California, USA using repeat photogrammetry, 2017–2023

Seacliff erosion in steep terrain poses major risks to transportation and critical infrastructure. In Big Sur, California, USA, seacliff erosion threatens the sustainability of the central coast stretch of California State Route 1, a transportation corridor that is critical to the region's economy. Published cliff retreat rates for the region range from 1 to 40 cm yr −1 , highlighting that high-resolution, process-based studies could enhance understanding of the causes of spatial and temporal variability. We quantified cliff erosion and investigated its drivers along ∼13 km of the Big Sur coastline at week–month timescales during the late fall to early spring wet seasons between January 2017 and June 2023 by analyzing 3D point clouds developed from aerial imagery using four-dimensional structure-from-motion (4D SfM) photogrammetry techniques. We calculated cliff face retreat rates of 2.23 ± 3.06 cm yr −1 (mean ±1 σ ), an order of magnitude lower than long-term estimated rates for the region (which included large deep-seated landslides), but in line with short-term rates reported across California. Change detection imagery comparison, cliff profiles through time, and statistical analysis reveal a cyclical cliff evolution process in which erosion by wave action at the cliff base destabilizes the cliff and primes it for subsequent failure during precipitation events. Although more erosion by volume could be attributed to precipitation-induced increases in soil moisture (784 m 3 km −1 yr −1 ) compared with erosion attributed to wave power (282 m 3 km −1 yr −1 ), our observations underscore the coupled nature of these processes in driving cliff evolution, consistent with established theory and observations.

California↗
USGS2025

Using gridded seismicity to forecast the long-term spatial distribution of earthquakes for the 2025 Puerto Rico and U.S. Virgin Islands National Seismic Hazard Model

Gridded (or background) seismicity models are a critical component of probabilistic seismic hazard assessments, accounting for off‐fault and smaller‐magnitude earthquakes. They are typically developed by declustering and spatially smoothing an earthquake catalog to estimate a long‐term seismicity rate that can be used to forecast future earthquakes. Here, we present new gridded seismicity models for use in the 2025 National Seismic Hazard Model (NSHM) for Puerto Rico and the U.S. Virgin Islands (PRVI). The previous PRVI NSHM was released in 2003, and our new models incorporate updates to both data and methodology. We utilize an updated earthquake catalog based on improved Puerto Rico Seismic Network data with newly characterized completeness epochs. The catalog is divided into crustal, subduction interface, and intraslab seismicity using new methods and Slab2 subduction zone geometries. To forecast the long‐term spatial distribution of earthquakes, we use an updated methodology developed for the 2023 U.S. 50‐state NSHM, considering three declustering methods and two spatial smoothing methods based on 2D Gaussian kernels. To adapt it for the complex seismotectonics of the region, we also adopt probabilistic methods to account for events with unknown depths and uncertainties in tectonic classification, and develop a new method for spatial scaling to counteract the effects of spatial variability in network coverage while maintaining the use of smaller events. Finally, we test the performance of these spatial models in forecasting the location of M w ≥ 5earthquakes in the region. Our updated methodology improves the representation of epistemic uncertainty relative to the 2003 model, and our results demonstrate the effectiveness of the new measures we have introduced to address heterogeneities in network detection and systematically evaluate forecast performance.

Puerto Rico, U.S. Virgin Islands↗
USGS2026

Treatability study to evaluate bioremediation of trichloroethene at Site K, former Twin Cities Army Ammunition Plant, Arden Hills, Minnesota, 2020–22

Executive Summary Chlorinated solvents, including trichloroethene (TCE) and other chlorinated volatile organic compounds (cVOCs), are widespread contaminants that can be treated by bioremediation approaches that enhance anaerobic reductive dechlorination. Reductive dechlorination can be enhanced either through the addition of an electron donor (biostimulation) or the addition of a known dechlorinating culture (bioaugmentation) along with an electron donor. Although bioremediation has been applied at many TCE-contaminated groundwater sites, application in source zones at sites where residual dense nonaqueous phase liquid (DNAPL) is present is more limited. In this study, laboratory and field treatability tests were completed to evaluate the potential application of anaerobic bioremediation for a shallow groundwater plume containing TCE in a perched alluvial aquifer at Site K, former Twin Cities Army Ammunition Plant, Arden Hills, Minnesota, which was on the National Priorities List as the New Brighton/Arden Hills Superfund site until 2019. In addition to the presence of residual DNAPL at the site, temporal variability in groundwater flow directions and input of oxygenated recharge were possible complicating factors for the application of enhanced anaerobic biodegradation in the shallow plume. The Site K plume extends beneath the footprint of Building 103, which was demolished in 2006, and soil excavations to a maximum depth of 6 feet (ft) below ground surface in 2014 were known to leave some deeper contaminated soil in place in the TCE source area. Groundwater treatment at the site, formalized as part of the 1997 Record of Decision, has been in operation since 1986 and consists of an extraction trench at the downgradient edge of the plume to collect groundwater, which is then pumped to an on-site air stripper. Groundwater concentrations in the plume have been relatively stable since treatment began, indicating a continued source of TCE in the aquifer. The desire for a destructive remedy that would enhance the removal of cVOCs in the aquifer at Site K and shorten the remediation timeframe led the U.S. Army to request that the U.S. Geological Survey conduct a groundwater treatability study to assess bioremediation. This report describes the U.S. Geological Survey bioremediation treatability study conducted during 2020–22, including pre-design site characterization to assist in formulating the bioremediation approach, laboratory experiments to support the design of the field pilot test, and implementation and 1-year performance monitoring results for the pilot test. Pre-design site characterization included the collection of soil cores for cVOC analysis and lithologic descriptions and the re-installment of three wells to obtain hydrologic measurements and initial groundwater chemistry. Relatively flat head gradients were measured at the site, and substantial decreases in water-level elevations occurred from spring to summer (May–July 2021). Continuous water-level monitoring indicated a rapid response to precipitation. Groundwater flow velocities were consistently less than 0.5 foot per day, and the pilot bioremediation test was therefore designed with short lateral distances (about 5 ft) between injection and individual monitoring points. Soil analyses confirmed that high volatile organic compound contamination was left in place in the source area. The highest concentrations were near or in clay at the base of the perched aquifer. Concentrations of cVOCs measured in the replaced wells were consistent with historical data and had a maximum TCE concentration of 57,700 micrograms per liter (μg/L), indicative of nearby residual DNAPL based on the general rule of observed concentrations exceeding 1 percent of solubility. The primary TCE daughter product detected was 1,2-cis-dichloroethene (cisDCE), which indicated limited reductive dechlorination in the plume. Groundwater in both the source and downgradient areas was relatively reducing during the pre-design characterization, particularly in the source area where methane concentrations greater than 400 μg/L were measured. Initial laboratory tests conducted using native aquifer microorganisms from the three replacement wells showed that anaerobic TCE biodegradation rates were low when biostimulated with the addition of sodium lactate as an electron donor, also known as a carbon donor, and resulted in the production of only cisDCE. Addition of a known dechlorinating culture, WBC-2, however, resulted in rapid biodegradation and production of ethene, verifying complete reductive dechlorination of TCE. Microcosms constructed with aquifer soil collected from the site were used to evaluate other electron donors besides lactate to support reductive dechlorination by WBC-2, including corn syrup as an alternative fast-release compound and whey, soy-based vegetable oil, and 3-D Microemulsion (Regenesis, San Clemente, California) as slow-release compounds. First-order rate constants for total organic chlorine removal in these WBC-2 amended microcosms were greatest with either lactate or vegetable oil as the donor, ranging between 0.061 and 0.047 per day or corresponding half-lives of 11–15 days. Testing of commercial products in other WBC-2-bioaugmented microcosms led to selection for the field pilot test of an emulsified vegetable oil product that also contained some sodium lactate as a fast-release donor. Delaying the addition of WBC-2 relative to the donor in the microcosms resulted in the most rapid overall biodegradation rates. The selected design for the pilot test utilized three separate test plots, each about 30-ft wide and 60-ft long: plots GS1 and GS2 in the source area of the plume and plot GS3 in the downgradient area of the plume near the excavation trench. Each test plot had one injection well, one monitoring well upgradient from the injection point, and 12 surrounding monitoring wells in a grid to capture variable groundwater flow directions. Donor injections, which included a bromide tracer, were completed in October 2021, immediately following baseline sampling, and the WBC-2 culture was injected about 40 days later, between November 30 and December 2, 2021. Performance monitoring conducted until December 2022 included hydrologic measurements and analyses of cVOCs, redox-sensitive constituents, dissolved organic carbon, bromide, volatile fatty acids, compound-specific carbon isotopes, and microbial communities. The biogeochemical data collected during the pilot tests in the three treatment plots showed that enhanced, complete reductive dechlorination of cVOCs in the groundwater was achieved in the GS1 and GS3 plots. In contrast, evidence of distribution of the injected amendments and subsequent biodegradation was limited in GS2, which was in an area of more heterogeneous soil lithology and low water table elevations. The molar composition of volatile organic compounds in the GS1 and GS3 plots was dominated by ethene in wells that were reached by the injected amendments by the end of the monitoring period. In the GS1 and GS3 plots, similar patterns were observed of cVOC concentrations decreasing to near detection levels, or below, at some wells sampled in July and October 2022, whereas ethene became dominant and indicated sustained complete reductive dechlorination. Baseline cVOC concentrations were more than a factor of 10 higher in the groundwater in the GS1 plot than in GS3, but no apparent inhibition of complete dechlorination occurred. As expected from the initial pre-design site data and the laboratory experiments, enhanced dissolution of residual DNAPL coupled to biodegradation was evident in the GS1 plot, where a marked increase in dichloroethene (DCE) above the initial baseline and upgradient TCE and DCE concentrations occurred. DCE concentrations subsequently declined where DNAPL dissolution was evident, concurrent with production of vinyl chloride and then predominantly ethene. Thus, overall biodegradation rates outpaced the DNAPL dissolution and desorption and DCE production in the source area. This success in complete degradation to predominantly ethene was achieved even in areas where the DCE concentrations reached a maximum of about 30,000 μg/L. Compound specific isotope analysis of carbon in TCE, cisDCE, trans-1,2-dichloroethene, and vinyl chloride was conducted to provide another line of evidence of the occurrence and extent of anaerobic biodegradation. Along a flow path in each plot that was affected by the injected amendments, carbon isotopes in the TCE and daughter cVOCs in the groundwater became isotopically heavier, indicating biodegradation. Enhanced biodegradation rates calculated from the field tests in GS1 and GS3 showed half-lives of 36.9–75.3 days for DCE degradation and 9.48–38.5 days for ethene production. Notably, these ethene production rates calculated from the field tests are consistent with the results of WBC-2-bioaugmented microcosms amended with either lactate or vegetable oil, which had half-lives for total organic chlorine removal that ranged from 11 to 15 days. These rates indicated rapid enhanced biodegradation, which is promising for application of a full-scale bioremediation remedy. Ultimately, however, the mass of residual or sorbed TCE in the aquifer that remains accessible for dissolution and biodegradation would likely control the time required for a full-scale bioremediation effort to achieve performance goals for TCE and cisDCE specified in the Record of Decision for Site K. The field pilot tests showed that the relatively low hydraulic head gradients and temporal changes in groundwater flow directions in the shallow aquifer would add complexity to a full-scale bioremediation effort. The radius of influence (ROI) at GS1 and GS3 (16.3 ft and 12.7 ft, respectively) were close to the design ROI of 15 ft. The estimated ROI at GS2 was about four times the design ROI, but may be less reliable at this location owing to groundwater flow direction. In addition, the low temperatures following WBC-2 injection in late November to early December 2021, in combination with the low hydraulic head gradients, were probably major factors in the delay observed before the onset of enhanced biodegradation following injection of the culture. Additional test injections could be beneficial to optimize the timing of donor and culture injections with the variable temperatures and hydraulic head in the shallow aquifer.

Minnesota↗
USGS2025

Estimating the probability of export restrictions to inform mineral criticality

As demand for advanced technologies rises, mineral commodities will increase in geopolitical importance. To assess risks associated with mineral commodity supply chain disruptions, governmental agencies and others have developed "criticality" assessments, with criticality described using the economic impact and probability of supply chain disruptions. In previous work, subjective supply risk indicators were developed to approximate this probability, typically combining several factors such as supply diversity and political stability of trading partners, where indicator weightings can substantially impact results. This work explicitly quantifies trade barrier probability using an ensemble of several machine learning classifiers, with probability estimates informed by exogenous variables such as prior trade barrier implementation and global export dominance. Major differences in the high-probability countries and commodities are observed across models, but the ensemble method highlights Indonesia, China, Tanzania, and the United States as particularly high risk. This approach enables a direct, quantitative, objective approach to assessing trade barrier probability, enhancing risk identification and prioritization for policymakers.

SSRN↗
USGS2026

Localization of spatiotemporally heterogeneous subsurface flows using autoencoder-based deep learning framework for time-lapse self-potential tomography

Self-potential (SP) monitoring has emerged as a valuable method for characterizing subsurface hydrogeological features and processes due to its sensitivity to fluid-induced electrokinetic effects. Despite advancements in SP inversion, challenges remain in imaging groundwater dynamics from SP activities due to complex hydrological settings and transient noise. In this study, a deep learning autoencoder (AE)-based framework is proposed for the spatiotemporal localization of subsurface fluid movement from time-lapse SP tomography. Temporal segments of time-lapse numerical inversions were first derived from long-term SP monitoring conducted from a floodplain site in Oak Ridge, Tennessee, known for active hyporheic exchange. Subsequently, AE models based on vision transformer (ViT), convolutional long short-term memory (ConvLSTM), convolutional neural network, and temporal convolutional network were individually trained and compared on the SP tomography segments for reconstruction performance. Finally, the reconstruction error over time serves as an anomaly score to identify moments of active SP variation, whereas spatial distributions of errors within these moments are analyzed to image and localize regions associated with anomalous subsurface fluid movement. The results demonstrate that ConvLSTM- and ViT-AE are most capable for the localization task with contrasting error distributions and consistent delineation of anomalies. Applying the method to both SP arrays parallel and perpendicular to the stream produced consistent anomaly zones near a fault or karst feature, validating the robustness and generalization of the approach. These results demonstrate the potential of the proposed framework as a scalable and interpretable tool for spatiotemporal analysis of subsurface flow dynamics in complex hydrogeological systems.

Tennessee↗
USGS2025

Comparative assessment of a restored and natural wetland using 13C-DNA SIP reveals a higher potential for methane production in the restored wetland

Wetlands are the largest natural source of methane (CH 4 ), a potent greenhouse gas produced by methanogens. Methanogenesis rates are controlled by environmental factors such as redox potential, temperature, and carbon and electron acceptor availability and are presumably dependent on the composition of the active methanogen community. We collected intact soil cores from a restored and natural freshwater depressional wetland on Maryland’s Delmarva Peninsula (USA) to assess the effects of wetland restoration and redox shifts on microbial processes. Intact soil cores were incubated under either saturated (anoxic) or unsaturated (oxic) conditions and amended with 13 C-acetate for quantitative stable isotope probing (qSIP) of the 16S rRNA gene. Restored wetland cores supported a distinct community of methanogens compared to natural cores, and acetoclastic methanogens putatively identified in the genus Methanosarcina were among the most abundant taxa in restored anoxic and oxic cores. The active microbial communities in the restored wetland cores were also distinguished by the unique presence of facultatively anaerobic bacteria belonging to the orders Firmicutes and Bacteroidetes . In natural wetland incubations, methanogen populations were not among the most abundant taxa, and these communities were instead distinguished by the unique presence of aerobic bacteria in the phyla Acidobacteria , Actinobacteria , and class Alphaproteobacteria . Iron-reducing bacteria, in the genus Geobacter , were active across all redox conditions in both the restored and the natural cores, except the natural oxic–anoxic condition. These findings suggest an overall higher potential for methanogenesis in the restored wetland site compared to the natural wetland site, even when there is evidence of Fe reduction.

Maryland↗
USGS2024

Onset of aftershocks: Constraints on the Rate-and-State model

Aftershock rates typically decay with time t after the mainshock according to the Omori–Utsu law, R (t)=K(c+t) −p ⁠ , with parameters K , c , and p . The rate‐and‐state (RS) model, which is currently the most popular physics‐based seismicity model, also predicts an Omori–Utsu decay with p = 1 and a c ‐value that depends on the size of the coseismic stress change. Because the mainshock‐induced stresses strongly vary in space, the c ‐value should vary accordingly. Short‐time aftershock incompleteness (STAI) in earthquake catalogs has prevented a detailed test of this prediction so far, but the newly developed a ‐positive method for reconstructing the true earthquake rate now allows its testing. Using previously published slip models, we calculate the coseismic stress changes for the six largest mainshocks in Southern California in recent decades and estimate the maximum shear as a scalar proxy of the coseismic stress tensor. Aftershock rates reconstructed for events in different stress ranges show that the rates follow a power law with p = 1 independent of stress with no clear sign of a c ‐value. The onset of the power‐law decay is abrupt and more delayed in areas with smaller stress changes. The observations do not necessarily contradict the RS model, as STAI limits the resolution for early aftershocks, and the RS model can reproduce the observations for specific Aσ values. However, the observations lead to strong constraints, namely Aσ <10 kPa and a power‐law decay of the background rate with distance to the fault, with exponent 2.7.

Seismological Research Letters↗
USGS2025

Characterizing the scale of regional landslide triggering from storm hydrometeorology

Rainfall strongly affects landslide triggering; however, understanding how storm characteristics relate to the severity of landslides at the regional scale has thus far remained unclear, despite the societal benefits that would result from defining this relationship. As mapped landslide inventories typically cover a small region relative to a storm system, here we develop a dimensionless index for landslide-inducing rainfall, A * , based on extremes of modeled soil water relative to its local climatology. We calibrate A * using four landslide inventories, comprising over 11 000 individual landslides over four unique storm events, and find that a common threshold can be applied to estimate regional shallow-landslide-triggering potential across diverse climatic regimes in California (USA). We then use the spatial distribution of A * , along with topography, to calculate the landslide potential area (LPA) for nine landslide-inducing storm events over the past 20 years, and we test whether atmospheric metrics describing the strength of landfalling storms, such as integrated water vapor transport, correlate with the magnitude of hazardous landslide-inducing rainfall. We find that although the events with the largest LPA do occur during exceptional atmospheric river (AR) storms, the strength of landfalling atmospheric rivers does not scale neatly with landslide potential area, and even exceptionally strong ARs may yield minimal landslide impacts. Other factors, such as antecedent soil moisture driven by storm frequency and mesoscale precipitation features within storms, are instead more likely to dictate the patterns of landslide-generating rainfall throughout the state.

California↗
USGS2024

Rainfall intensification amplifies exposure of American Southwest to conditions that trigger postfire debris flows

Short-duration, high-intensity rainfall can initiate deadly and destructive debris flows after wildfire. Methods to estimate the conditions that can trigger debris flows exist and guidance to determine how often those thresholds will be exceeded under the present climate are available. However, the limited spatiotemporal resolution of climate models has hampered efforts to characterize how rainfall intensification driven by global warming may affect debris-flow hazards. We use novel, dynamically downscaled (3.75-km), convection-permitting simulations of short-duration (15-min) rainfall to evaluate threshold exceedance for late 21st-century climate scenarios in the American Southwest. We observe significant increases in the frequency and magnitude of exceedances for regions dominated by cool- and warm-season rainfall. We also observe an increased frequency of exceedance in regions where postfire debris flows have not been documented, and communities are unaccustomed to the hazard. Our findings can inform planning efforts to increase resiliency to debris flows under a changing climate.

Arizona, California, Nevada, New Mexico, Utah↗
USGS2026

Apatite (U-Th)/He thermochronology from Marsh Creek anticline reconciles Cenozoic and Holocene strain patterns and elucidates the Cenozoic canning displacement zone of Arctic Alaska

The northeastern Brooks Range is a vast fold-thrust belt that records the northernmost expression of Cenozoic deformation in Alaska, likely with complex dynamic linkages to distant tectonic elements. The principal focus of this manuscript is the Marsh Creek anticline (MCA), which is a large and recently active structure within the broader northeastern Brooks Range. The MCA is a complex, polyphase structure mostly buried beneath the Arctic coastal plain. We present apatite (U-Th)/He (AHe) dates from six outcrop samples, all near the structural crest of the anticline, as well as supporting detrital zircon U/Pb (ZUPb) age spectra from five outcrop samples. The AHe sample transect spans ~3 km of structural relief. The four shallowest AHe samples are not reset. However, the two deepest samples exhibit low intrasample variability among replicates and strong date versus grain radius correlations, suggesting post-depositional burial heating to within the AHe partial retention zone. Inverse thermal history modeling of these data resolves a protracted phase of exhumational cooling that began around 33–28 Ma. Detrital ZUPb age spectra from Paleogene strata support an interpretation of diverse inherited thermal histories recorded by the AHe data, and they are also consistent with extant reconstructions of Paleogene sediment dispersal. Placed into regional context, the results constrain a phase of increasing structural relief in the northeastern Brooks Range fold-thrust belt that began by the Eocene. The relief developed in a zone characterized by diminishing structural relief to the west called the Canning displacement zone. Although this zone was active since the Eocene, several structures including the studied part of the MCA were rapidly exhumed in the Oligocene. Lastly, we argue for kinematic compatibility between Cenozoic contraction and active strike-slip in the northeastern Brooks Range, with the temporal change in strain style likely reflecting the increase in structural relief and topography.

Alaska↗
USGS2000

Ground-water quality atlas of Oakland County, Michigan

The U.S. Geological Survey (USGS), in cooperation with Oakland County Health Division (OCHD), collected 140 water samples from 38 wells in Oakland County during 1998 to better understand ground-water quality. OCHD had observed temporal variations in concentrations of various constituents, so two additional sets of samples were collected to evaluate potential short-term variability related to sample collection procedures and long-term seasonal variability. Replicate samples from 28 wells were analyzed in the Michigan Department of Environmental Quality (MDEQ) Drinking Water Laboratory to compare MDEQ&rsquo;s analytical results to those obtained from the USGS National Water Quality Laboratory. Several additional databases describing population, land use, water supply, soils, geology, and flows of ground water and surface water are presented in the first part of the report to assist in interpreting the water-quality data. Maps created from these databases are provided in the first portion of the report as an extension of the study-area description. The U.S. Environmental Protection Agency (USEPA) has established Maximum Contaminant Levels (MCL) and Secondary Maximum Contaminant Levels (SMCL) for which samples were analyzed in this study. Water from the 38 wells sampled by the USGS did not exceed the SMCL or MCL for sulfate, fluoride, or nitrite. However, water from 26 wells exceeded the SMCL for iron, water from 12 wells exceeded the SMCL for manganese, and water from 12 wells exceeded the SMCL for dissolved solids. Water from two wells exceeded the MCL for nitrate, although nitrate concentrations in water from most wells was below the detection limit. Water from seven wells exceeded the SMCL for chloride, and water from all wells contained detectable concentrations of chloride. Water from five wells exceeded the MCL for arsenic, and most of the wells sampled contained detectable concentrations of arsenic. These five wells were identified from previous MDEQ analyses to have elevated arsenic concentrations, and were sampled to obtain additional chemistry information. Replicate samples were collected from 26 of the 38 wells for analysis at the MDEQ Drinking-Water Laboratory to compare the results with the USGS National Water Quality Laboratory. The results of the replicate analyses indicate close agreement between the laboratories, with mean differences for nitrate, chloride, and arsenic of 0.10 milligrams per liter (mg/L) as nitrogen, 6.8 mg/L, and 0.0008 mg/L, respectively between the USGS and MDEQ analyses. Potential health effects associated with ingesting nitrate, chloride, and arsenic are provided with the water-quality data, along with references for further information.

Michigan↗
USGS2025

False positives in the identification of dynamic earthquake triggering

Dynamic earthquake triggering is commonly identified through the temporal correlation between increased seismicity rates and global earthquakes that are possible triggering events. However, correlation does not imply causation. False positives may occur when unrelated seismicity rate changes coincidently occur at around the time of candidate triggers. We investigate the expected false positive rate in Southern California with global M ≥ 6 earthquakes as candidate triggers. We compute the false positive rate by applying the statistical tests used by DeSalvio and Fan (2023), https://doi.org/10.1029/2023jb026487 to synthetic earthquake catalogs with no real dynamic triggering. We find a false positive rate of ∼3.5%–8.5% when realistic earthquake clustering is present, consistent with the 95% confidence typically used in seismology. However, when this false positive rate is applied to the tens of thousands of spatial-temporal windows in Southern California tested in DeSalvio and Fan (2023), https://doi.org/10.1029/2023jb026487 , thousands of false positives are expected. The expected false positive occurrence is large enough to explain the observed apparent triggering following 70% of large global earthquakes (DeSalvio & Fan, 2023, https://doi.org/10.1029/2023jb026487 ), without requiring any true dynamic triggering. Aside from the known triggering from the nearby El Mayor-Cucapah, Mexico, earthquake, the spatial and temporal characteristics of the reported triggering are indistinguishable from random false positives. This implies that best practice for dynamic triggering studies that depend on temporal correlation is to estimate the false positive rate and investigate whether the observed apparent triggering is distinguishable from the correlations that may occur by chance.

JGR Solid Earth↗
USGS2024

Layered intrusions in the Precambrian: Observations and perspectives

Layered intrusions are plutonic bodies of cumulates that form by the crystallization of mantle-derived melts. These intrusions are characterized by igneous layering distinguishable by shifts in mineralogy, texture, or composition. Layered intrusions have been fundamental to our understanding of igneous petrology; however, it is their status as important repositories of critical metals – such as platinum-group elements, chromium, and vanadium – that has predominantly driven associated research in recent decades. Many layered intrusions were emplaced during the Precambrian, predominantly at the margins of ancient cratons during intervals of supercontinent accretion and destruction. It appears that large, layered intrusions require rigid crust to ensure their preservation, and their geometry and layering is primarily controlled by the nature of melt emplacement. Layered intrusions are best investigated by integrating observations from various length-scales. At the macroscale, intrusion geometries can be discerned, and their presence understood in the context of the regional geology. At the mesoscale, the layering of an intrusion may be characterized, intrusion-host rock contact relationships studied, and the nature of stratiform mineral occurrences described. At the microscale, the mineralogy and texture of cumulate rocks and any mineralization are elucidated, particularly when novel microtextural and mineral chemical datasets are integrated. For example, here we demonstrate how mesoscale observations and microscale datasets can be combined to understand the petrogenesis of the perplexing snowball oiks outcrop located in the Upper Banded Series of the Stillwater Complex. Our data suggest that the orthopyroxene oikocrysts did not form in their present location, but rather formed in a dynamic magma chamber where crystals were transported either by convective currents or within crystal-rich slurries. Critical metals may be transported to the level of a nascent intrusion as dissolved components in the melt. Alternatively, ore minerals are entrained from elsewhere in a plumbing system, potentially facilitated by volatile-rich phases. There are many ore-forming processes propounded by researchers to occur at the level of emplacement; however, each must address the arrival of the ore mineral, its concentration of metals, and its accumulation into orebodies. In this contribution, several of these processes are described as well as our perspectives on the future of layered intrusion research.

Precambrian Research↗