USGS Science⌕ Search

SEARCH · USGS Science

Results for “Geostandards and Geoanalytical Research”

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,660 records · Page 54Linked to original sources
USGS1976

Transport and dispersion of fluorescent tracer particles for the dune-bed condition, Atrisco Feeder Canal near Bernalillo, New Mexico

A tracer technique in which mineral particles were coated with fluorescent dyes was used to study the rates of transport and dispersion of sediment particles`of various diameters and specific gravities for a dune-bed condition in an alluvial channel. The experiment was conducted in the Atrisco Feeder Canal near Bernalillo, N. Mex., between May 1 and July 14, 1967. A continuous point source of tracers, approximated by injections at 5- or 10-minute intervals, was maintained for 7 days so that the steady-dilution procedure could be used to calculate the transport rate of bed material. After termination of the injection process, the spatial-integration procedure was used to follow the movement of the tracers downstream and to calculate the transport rate. Samples of the bed material in transport and the accompanying tracers moving along the surface of the dune bed were obtained periodically throughout the study with the "dustpan" sampler especially designed for fluorescent tracer studies. In addition, the spatial distributions of the tracers in the dune bed were determined three times during the study by core sampling.

New Mexico↗
USGS2025

Reply to, “Comment on ‘The 1886 Charleston, South Carolina, earthquake: Relic railroad offset reveals rupture,’ by Roger Bilham and Susan E. Hough”

We welcome this opportunity to respond to Pratt et al. (2024) (hereinafter P24). Bilham and Hough (2023) proposed a “first-cut” elastic deformation model for the 1886 earthquake, a quantitative source model constrained by identified coseismic constraints. A key observation was the measurement of a lateral offset of a railroad line south of Summerville, leading to a model with predominately dextral slip and minor convergence, from which we concluded that active faulting had raised the Penholoway Marine Terrace >6 m since ∼770 ka. P24 questioned these constraints and proposed an alternative rupture model with predominantly reverse slip. This alternative model is neither consistent with coseismic constraints nor with other geophysical data. In a revised model presented here, we recognize that uplift of the Penholoway Terrace is confined to the eastern edge of the terrace, which we conclude results from active folding and tectonic transpression centered on the dextral fault that offset the railroad in 1886.

South Carolina↗
USGS2026

Altitude of the potentiometric surface and depth to water in the Mississippi River Valley alluvial aquifer, spring 2024

Potentiometric-surface and depth-to-water maps for spring 2024 were created for the Mississippi River Valley alluvial aquifer using groundwater-altitude data from 1,151 wells completed in the Mississippi River Valley alluvial aquifer and from the altitude of the top of the water surface in area rivers from 160 U.S. Geological Survey and U.S. Army Corps of Engineers streamgages. The potentiometric-surface and depth-to-water maps for 2024 were created to support investigations to characterize the Mississippi River Valley alluvial aquifer as part of the U.S. Geological Survey Water Availability and Use Science Program. Sufficient data were available to map the potentiometric surface and depth to water of the Mississippi River Valley alluvial aquifer for spring 2024 for about 83 percent of the aquifer area. The potentiometric contours ranged from 0 to 330 feet (ft) above the North American Vertical Datum of 1988. The regional direction of the groundwater gradient was generally to the south-southwest, except in areas of groundwater-altitude depressions, where the groundwater gradient direction was into the depression, and near rivers, where the groundwater gradient direction was either from the aquifer to the river or from the river into the aquifer. Groundwater depressions in the potentiometric-surface map are in the lower one-half of the Cache region and in most of the Grand Prairie and Delta regions. Depth to water by well in the Mississippi River Valley alluvial aquifer in spring 2024 ranged from 0.27 ft above land surface to 145.65 ft below land surface.

Arkansas, Illinois, Kentucky, Louisiana, Mississip↗
USGS2024

Will there be water? Climate change, housing needs, and future water demand in California

Climate change in California is expected to alter future water availability, impacting water supplies needed to support future housing growth and agriculture demand. In groundwater-dependent regions like California's Central Coast, new land-use related water demand and decreasing recharge is already stressing depleted groundwater basins. We developed a spatially explicit state-and-transition simulation model that integrates climate, land-use change, water demand, and groundwater gain-loss to examine the impact of future climate and land use change on groundwater balance and water demand in five counties along the Central Coast from 2010 to 2060. The model incorporated downscaled groundwater recharge projections based on a Warm/Wet and a Hot/Dry climate future from a spatially explicit hydrological process-based model. Two urbanization projections from a parcel-based, regional urban growth model representing 1) recent historical and 2) state-mandated housing growth projections were used as alternative spatial targets for future urban growth. Agricultural projections were based on recent historical trends from remote sensing data. Annual projected changes in groundwater balance were calculated as the difference between land-use related water demand, based on historical estimates, and climate-driven recharge plus agriculture return flows. Results indicate that future changes in climate-driven groundwater recharge, coupled with cumulative increases in agricultural water demand, result in overall declines in future groundwater balance, with a Hot/Dry future resulting in cumulative groundwater decline in all but Santa Cruz County. Cumulative declines by 2060 are especially prominent in San Luis Obispo (−2.9 to −5.1 Bm 3 ) and Monterey counties (−6.5 to −8.7 Bm 3 ), despite limited changes in agricultural water demand over the model period. These two counties show declining groundwater reserves in a Warm/Wet future as well, while San Benito and Santa Barbara County barely reach equilibrium. These results suggest future groundwater supplies may not be able to keep pace with regional demand and declining climate-driven recharge, resulting in a potential reduction in water security in the region. However, our county-scale projections showed new housing and associated water demand does not conflict with California's groundwater sustainability goals. Rather, future climate coupled with increasing agricultural groundwater demand may reduce water security in some counties, potentially limiting available groundwater supplies for new housing.

California↗
USGS2025

Groundwater dominates snowmelt runoff and controls streamflow efficiency in the western United States

Climate change in seasonally snow-covered mountain catchments is reducing water supply and decreasing streamflow predictability. Here, we use tritium age dating to show that contrary to the common assumption that snowmelt quickly contributes to runoff, streamflow during snowmelt in western US catchments is dominated by older groundwater. The average age of streamwater during snowmelt runoff (5.7 ± 4.3 years) was intermediate to the average age of groundwater (10.4 ± 4.5 years) and recent precipitation, indicating that 58% (±34%) of snowmelt runoff was derived from groundwater. Water ages, streamflow, and groundwater storage were mediated by bedrock geology: low-permeability hard rock/shale catchments exhibited younger ages, less storage, and more efficient streamflow generation than high-permeability sandstone/clastic catchments. Our results demonstrate that snowmelt runoff is the result of multiple prior years of climate mediated by groundwater storage. Including these interactions will be crucial for predicting water resources as climate and landscape changes accelerate.

western united States↗
USGS2025

Gravity and magnetic surveys of the Skaergaard intrusion, East Greenland

Aeromagnetic and gravity surveys of the Skaergaard intrusion in East Greenland were carried out in July–August 1971 as part of a grant to the University of Oregon Center for Volcanology to refine the models of crystallization and differentiation of the intrusion, specifically to test whether the intrusion is underlain by dense rocks of a reservoir 20 kilometers (km) thick (referred to as a “hidden zone”). The Skaergaard intrusion is a source of platinum group elements that are critical mineral resources for many technologies, and because no new data have been collected these legacy datasets remain a valuable asset. The total-intensity aeromagnetic survey was flown in early July 1971 with a proton precession magnetometer at a constant barometric altitude of 1.5 km (5,000 feet) with a nominal line spacing of 1 km. Two gravimeters were used to acquire 168 stations of which 86 were at known altitudes (mainly sea level) and 82 had altitudes measured by altimetry in late July–August 1971. Finally, a north-south ground vertical-intensity magnetic traverse was completed across the intrusion together with collection of oriented hand specimens. The hand specimens were measured for remnant magnetization and density, along with density measurements of more specimens collected by expedition geologists for other purposes. The intrusion is composed of layered gabbro with extensive crystal fractionation that is dense and strongly reversely polarized. After terrain correction and standard Bouguer gravity reduction, the gravity anomaly dataset was corrected for all rock above sea level using the density measurements of the various zones of the intrusion and the topographic and geologic maps (variable density Bouguer gravity reduction). A large regional gradient in the gravity anomaly data was removed using orthogonal polynomial fitting to the gridded data. The zonal volumes of rock below sea level were calculated from the dipping polygonal layer gravity model of the intrusion below sea level and combined with elliptic cross–section cylinders for the various zones above sea level to approximate the original zonal volumes of the intrusion. The residual gravity anomaly of 18–20 milligals (mGal) was only about half of the expected anomaly if a large hidden zone proposed from petrologic considerations were present, and both two-dimensional and three-dimensional models imply that the exposed series of intrusion zones explain the gravity anomaly by their down-dip extension below sea level together with a small hidden-zone volume. A three-dimensional model of the exposed rocks and their down-dip extension below sea level also can account for the aeromagnetic anomaly with little or no requirement for hidden-zone rock. The middle and upper zone units of the intrusion contain the most magnetite and account for most of the aeromagnetic anomaly.

Skaergaard intrusion↗
USGS2025

Ages and trace element fertility of porphyry-related mineralization in the Philipsburg polymetallic district, Montana, with a comparison to Butte

The Philipsburg mining district is a Mo-Cu porphyry system with associated Cordilleran polymetallic veins. Geochronology was employed to date the porphyry (~66 Ma, U/Pb in zircon) and molybdenite mineralization from the veins (~76 Ma, Re-Os). Age results suggest that the two-mineralization events model proposed by Lund et al. (2018) for the Butte district can be applied to the Philipsburg deposit. Furthermore, common fertility indicators from zircon trace elements were analyzed with variable success to characterize the existing Mo-Cu mineralization of the Philipsburg porphyry. This raises concerns about the benefits of this method and/or the importance of sampling, which could impact exploration for similar porphyry deposits.

Montana↗
USGS2024

Static and dynamic strain in the 1886 Charleston, South Carolina, earthquake

During the 1886 Mw 7.3 Charleston, South Carolina, earthquake, three railroads emanating from the city were exposed to severe shaking. Expansion joints in segmented railroad tracks are designed to allow railroad infrastructure to withstand a few parts in 10,000 of thermoelastic strain. We show that, in 1886, transient contractions exceeding this limiting value buckled rails, and transient extensions pulled rails apart. Calculated values for dynamic strain in the meizoseismal region are in reasonable agreement with those anticipated from the relation between strain and moment magnitude proposed by Barbour et al. (2021) and exceed estimated tectonic strain released by the earthquake by an order of magnitude. Almost all of the documented disturbances of railroad lines, including evidence for shortening of the rails, can thus be ascribed to the effects of dynamic strain changes, not static strain. Little or no damage to railroads was reported outside the estimated 10 −4 dynamic strain contour. The correspondence between 10 −3 and 2×10 −4 contours of dynamic strain and Mercalli intensity 9 and 8, anticipated from the dependence of each quantity on peak ground velocity, suggests it may be possible to use railroad damage to quantitatively estimate shaking intensity. At one location, near Rantowles, ≈20 km west of Charleston, a photograph of buckled track taken one day after the earthquake has been cited as evidence for shallow dextral slip and has long focused a search for a causal fault in this region. Photogrammetric analysis reveals that the buckle was caused by transient contraction of <10 cm with no dextral offset. Our results further weaken the evidence for faulting in the swamps and forests south of the Ashley River in 1886, hitherto motivated by the photograph and limited macroseismic evidence for high‐intensity shaking.

South Carolina↗
USGS2024

Improving crop-specific groundwater use estimation in the Mississippi Alluvial Plain: Implications for integrated remote sensing and machine learning approaches in data-scarce regions

Study region The Mississippi Alluvial Plain (MAP) in the United States (US). Study focus Understanding local-scale groundwater use, a critical component of the water budget, is necessary for implementing sustainable water management practices. The MAP is one of the most productive agricultural regions in the US and extracts more than 11 km 3 /year for irrigation activities. Consequently, groundwater-level declines in the MAP region pose a substantial challenge to water sustainability, and hence, we need reliable groundwater pumping monitoring solutions to manage this resource appropriately. New hydrological insights for the region We incorporate remote sensing datasets and machine learning to improve an existing lookup table-based model of groundwater use previously developed by the U.S. Geological Survey (USGS). Here, we employ Distributed Random Forests, an ensemble machine learning algorithm to predict annual and monthly groundwater use (2014–2020) throughout this region at 1-km resolution, using pumping data from existing flowmeters in the Mississippi Delta. Our model compares favorably with the existing USGS model, with higher R 2 (0.51 compared to 0.42 in the previous model), and lower root mean square error (RMSE) and mean absolute error (MAE)— 0.14 m and 0.09 m, respectively in our model, compared to 0.15 m and 0.1 m in the previous model. Therefore, this work advances our ability to predict groundwater use in regions with scarce or limited in-situ groundwater withdrawal data availability.

Journal of Hydrology Regional Studies↗
USGS2025

‘The fish that stop’: Drivers of historical decline for Pacific cod and implications for modern management in an era of rapidly changing climate

n the Gulf of Alaska, a series of marine heat waves depleted Pacific cod ( Gadus macrocephalus ) biomass to the lowest abundance ever recorded and led to the fishery’s closure in 2020. Although the fishery has been productive for decades, this collapse may have historical precedents. Traditional knowledge holders refer to cod as ‘the fish that stop’, and there is a suggested period of decline in the 1930s. Here we conduct a catch reconstruction of the early commercial fishery (1864–1950), confirming a rapid catch decline in the 1920s and 1930s. Next, we evaluate evidence for possible drivers. We document changes to demand and technology that contributed to declining catch. However, we also find both qualitative and quantitative evidence of depletion, suggesting catch declines were not driven entirely by social factors. Overfishing may have contributed to localized catch declines as evidenced by declining catch rates in heavily fished localities. We also find evidence for climate as a driver of regional decline, with the period of catch decline characterized by up to 2°C higher temperatures as compared to the earlier period of high fisheries production. Our analysis underscores the importance of understanding long-term drivers of fisheries productivity and the value of linking fisheries and climate histories.

Alaska↗
USGS2025

Assessing nonpoint-source uranium pollution in an irrigated stream-aquifer system

Uranium (U) in rocks and soils of arid and semi-arid environments can be mobilized by irrigation and fertilization, posing environmental and health risks. Elevated U, along with selenium (Se) and nitrate (NO 3 ) co-constituents, necessitates careful monitoring and management. We developed a distributed-parameter numerical model to assess U pollution in an irrigated stream-aquifer system, applying it to a 552 km 2 region in Colorado's Lower Arkansas River Valley (LARV) over 14 years. A MODFLOW model, describing groundwater and stream flow, was coupled with an RT3D-OTIS model to portray reactive U transport. Calibration using the PESTPP-iES iterative ensemble smoother (iES) software indicated good agreement with observed U concentrations. The model revealed substantial and variable U levels across the LARV, highlighting potential hotspots and possible contributing factors, such as geological composition of the bedrock and near-surface shale and aquifer sediments derived from them, irrigation practices, and riparian landscape. U levels exceed the chronic standard (85th percentile = 30 μg/L, set by the US Environmental Protection Agency), which is the permissible regulatory threshold, in groundwater across 44 % of the region and along the river by an average factor of 2.9. Simulated average U concentrations in the non-riparian aquifer and river are 124 μg/L and 60 μg/L, respectively, compared with 112 μg/L and 62 μg/L for measured values. The average 85th percentile U concentration is 222 μg/L in the aquifer and 82 μg/L in the river. Average simulated U mass loading to the river is 0.17 kg/day per km, compared to an estimated 0.23 kg/day per km. Findings provide a baseline for comparing future simulated outcomes of alternative best management practices (BMPs) for U pollution mitigation and offer a methodology applicable to other irrigated regions.

Colorado↗
USGS1990

Effects of storm-water runoff on local ground-water quality, Clarksville, Tennessee

Storm-related water-quality data were collected at a drainage-well site and at a spring site in Clarksville, Tennessee, to define the effects of storm-water runoff on the quality of ground water in the area. A dye-trace test verified the direct hydraulic connection between the drainage well and Mobley Spring. Samples of storm run off and spring flow were collected at these sites for nine storms during the period February to October 1988. Water samples were collected also from Mobley Spring and two other springs and two observation wells in the area during dry-weather conditions to assess the general quality of ground water in an urban karst terrain. Evaluation of the effect of storm-water runoff on the quality of local ground water is complicated by the presence of other sources of contaminants in the area Concentrations and load for most major constituents were much smaller in storm-water runoff at the drainage well than in the discharge of Mobley Spring, indicating that much of the chemical constituent load discharged from the spring comes from sources other than the drainage well. However, for some of the minor constituents associated with roadway runoff (arsenic, copper, lead, organic carbon, and oil and grease), the drainage well contributed relatively large amounts of these constituents to local ground water during storms. The close correlation between concentrations of total organic carbon and concentrations of most trace metals at the drainage-well and Mobley Spring sites indicates that these constituents are transported together. Many trace metals were flushed early during each runoff event. Mean storm loads for copper, lead, zinc, and four nutrient species (total nitrogen, ammonia nitrogen, total phosphorus, and orthophosphorus) in storm-water runoff at the drainage-well site were lower than mean storm load predicted from an existing regression model. The overprediction by the model may be a result of the small size of the drainage area relative to the range of drainage areas used in the development of the models, or to the below-normal amounts of rainfall during the period of sampling for this investigation. Loads& in storm-water runoff for 22 constituents were extrapolated from sampled storms to total loads for the period February to October 1988. Calculated loads for trace metals for the period ranged from 0.030pound.s for cadmium to 12pound.s for strontium. Loads of the primary nutrients ranged from 0.97pounds for nitrite as nitrogen to 34pounds of organic nitrogen. Storm-water quality at the drainage-well and Mobley Spring sites was compared to background water quality of the local aquifer; as characterized by dry-weather samples from three springs and two observation wells in the Clarksville area. Concentrations of total-recoverable cadmium, chromium, copper, lead, and nickel were higher in many stormwater samples from both the drainage-well and Mobley Spring sites than in samples from any other site. In addition, concentrations of total organic carbon, methylene blue active substances, and total-recoverable oil and grease were generally higher in storm-water samples from the drainage-well site than in any ground-water sample. Densities of fecal coliform and fecal streptococcus bacteria and concentrations of total recoverable iron, manganese, and methylene blue active substances in storm samples from the drainage-well site exceeded the maximum contaminant levels listed in Tennessee’s drinking-water standards (1988) by as much as 2,500 and 5,500 colonies per 100 milliliters, and 2.7, 0.29, and 0.05 milligrams per liter, respectively. Densities of fecal coliform and fecal streptococcus bacteria and concentrations of total-recoverable iron, manganese, and lead in storm samples from Mobley Spring exceeded the maximum contaminant levels by as much as 500 and 4,500 colonies per 100 milliliters, and 18.7,0.65, and 0.02 milligrams per liter, respectively. For iron, manganese, and bacteria, these undesirable levels are not necessarily attributable to storm-water recharge, because concentrations of these constituents also exceeded drinking-water standards in one or more of the dry-weather samples from selected springs and observation wells in the area.

Tennessee↗
USGS2025

Assimilation of reduced carbon triggers platinum alloy saturation in mafic and ultramafic magmas

It is generally observed that magmatic sulfide ores have higher ratios of Pd/Pt than the mantle-like values of their parental magmas. This discrepancy has defied simple explanation because the partitioning behavior of both elements between sulfide and silicate liquids is very similar. Assimilation of sulfur- and carbon-rich country rocks by mafic and ultramafic magmas is considered a critical, if not essential, step in the formation of magmatic base metal sulfide deposits. Although there is general consensus that the assimilation of external sulfur and carbon promotes sulfide saturation, the effect of carbon assimilation on the solubilities of platinum-group elements in natural S-bearing silicate melt has been overlooked. In this study, we investigate the variations of platinum and palladium solubilities during assimilation of graphite and methane through thermodynamic modeling, in comparison with data from an array of highly distinctive magmatic sulfide ore systems representing ages from Archean to Paleozoic, melt compositions from komatiite to basalt, and magmatic settings including lavas, hypabyssal intrusions, plutonic continental arc roots, and plutonic layered intrusions, namely: Raglan, Norilsk-Talnakh, Lac des Iles, and the J-M Reef of the Stillwater Complex. We model assimilation-fractional crystallization processes to estimate the reduction of oxygen fugacity ( ⁠fO2⁠ ) of the melt due to incorporation of graphite and methane. The simulations show that although Pd remains highly soluble during the progressive assimilation of reduced carbon, Pt solubility decreases significantly as the silicate melt becomes increasingly reduced. With less than 8 % of sediment assimilation, Pt alloy may saturate and then deviate from sulfide-undersaturated silicate melts, concomitantly increasing the Pd/Pt value of the remaining melts of the Raglan and Norilsk-Talnakh systems. For the Lac des Iles and Stillwater systems, a higher extent of assimilation is needed to reach Pt saturation because of the relatively carbon-poor nature of the lower crustal rocks. The assimilation of methane volatiles is shown to be more effective than graphite assimilation, and it provides a pathway to Pt alloy fractionation in the absence of detectable amounts of bulk host-rock assimilation. High Pd/Pt values have been documented in many world-class magmatic sulfide deposits whose parental magmas have demonstrably experienced crustal contamination. Our model suggests that although anomalous Pd/Pt values may be explained by other mechanisms such as incongruent melting of preexisting sulfide or differences in the diffusivities of the metals within achieving equilibration, the assimilation of graphite or methane may play an important role in the global occurrence of magmatic sulfide ores with elevated Pd/Pt values.

Economic Geology↗
USGS2025

Relating systematic molecular and textural properties of graptolite pyrolyzed via gold tube hydrous pyrolysis: Implications for thermal proxies in lower Paleozoic marine shales

A series of gold tube pyrolysis experiments (72 h, 300–550 °C, 50 MPa) conducted on a graptolite-rich lower Paleozoic marine shale generated pyrolysis residues for a comprehensive evaluation of the molecular and structural variability of three types of graptolite periderm. Organic petrology, Raman spectroscopy, and field emission scanning electron microscopy (FE-SEM) with energy dispersive spectroscopy (EDS) were combined to evaluate the thermal evolution process. The three types of graptolite periderm, namely granular, non-granular, and nodular graptolite, were analyzed by Raman spectroscopy wherein point measurements were obtained after the maceral was identified and the location verified by organic petrology. Distinct thermal evolution pathways among non-granular, granular, and nodular graptolite periderms were recorded. The evolution patterns of the Raman parameters, particularly D1 and G bands, highlight the differences in geochemical composition of the graptolite periderm types and the alteration of molecular structure with increasing thermal maturity. Raman parameters D1 (position of the D1 peak), G-FWHM (full width at half maximum of the G peak), and ratios D1-FWHM/G-FWHM (full width at half maximum of the D1 peak ratioed to G-FWHM) and A D1 /A G (ratio of D1 and G peak intensities) showed effectiveness in assessing thermal maturity. Bireflectance with increasing gold tube pyrolysis temperature followed a hierarchy: non-granular > granular > nodular, reflecting different molecular alignment intensities. Qualitative FE-SEM evaluation showed that fine-grained mineral inclusions (primarily Fe-sulfide as determined via EDS) were associated with the graptolite populations, with granular graptolite containing greater amounts of coarser-grained (e.g., ∼300–1400 nm) mineral inclusions relative to non-granular and nodular graptolite, which contain finer-grained (e.g., ∼100–200 nm) inclusions difficult to resolve with optical microscopy. These findings are investigated to highlight the mechanisms that drive organic matter evolution within graptolite during thermal maturation, as well as to explore some of the limitations of using spectroscopic parameters as thermal maturity proxies.

International Journal of Coal Geology↗
USGS2025

Global maps of critical mineral production in 2023

Introduction The global production of many mineral commodities, especially critical minerals, is concentrated in a few countries that have mineral resources and the infrastructure necessary to mine and process those resources. For this reason, the type and amount of mineral production differ by country. For example, many countries produce such metallic ores as gold and silver, whereas only a few countries produce magnesium, niobium, platinum-group metals, and rare earths. The concentration of mining and processing in certain countries necessitates the existence of a global supply chain. A mineral supply chain is the sequence of mining and processing of minerals and manufacturing of products. Mineral supply chains are global in scale, complex, and dynamic. Supply chain data can be used to understand how a country’s mineral resources and various economic, technical, and environmental factors affect the complexity of global supply chains. This fact sheet summarizes the world’s leading countries (those accounting for 5 percent or more of a commodity’s global production in 2023) for production of select mineral commodities (mainly critical minerals) in the mining and processing stages. These countries and the minerals they produce are synthesized on global maps to communicate the status of, and potential risk to, mineral commodity supply chains from geographic production concentration. Trade data from United Nations Statistics Division (2025) is used to support assessments of the observed production data.

Fact Sheet↗
USGS2025

Local, regional, and distal recordings of seismic unrest at Tau Island volcano, American Samoa

A seismic swarm near Taʻū Island, a volcanic island in eastern American Samoa, occurred from July to October 2022. The earliest unrest was noted as felt shaking reports in late July, and instrumentation varied in the beginning of the sequence as the U.S. Geological Survey (USGS) Hawaiian Volcano Observatory responded by installing temporary and then permanent seismometers to monitor the activity. This network variability made it difficult to characterize the earliest seismicity and contextualize the entire sequence to discriminate between an underlying tectonic or volcanic source. Here, we present results analyzing hydroacoustic detections from an International Monitoring System hydrophone array near Wake Island, 4500 km northwest of Taʻū Island volcano. Using least-squares beamforming analysis, we create a catalog of T-wave detections from the direction of Taʻū Island to track the earthquakes, some of which were located by the USGS National EarthquakeInformation Center. Both the rate and hydroacoustic pressures, which we interpret as a proxy for earthquake size, gradually increased from late July to August, peaking on August 19 (rate) and August 24 (size), before decreasing to background in late September. Minutes-long bursts of tremor were also contemporaneously recorded as local network data became avail-able on August 20. Tremor activity continued throughout the rest of August, peaking on August 25, before ending in earlySeptember. These tremor bursts were band-limited to ~ 1–5 Hz and recorded as S waves at a regional station on the island of Upolu in Samoa, 250 km to the west of Taʻū Island. Our results do not constrain the tremor locations, but comparisons of earthquake and tremor reduced displacements recorded locally and regionally suggest a deeper tremor source. We interpret the increase in earthquake size and rate, together with the occurrence, characteristics, and relative depth of the tremor to be the result of magmatic activity beneath Taʻū Island volcano.

American Samoa, Taʻū Island↗
USGS2025

Divide and conquer: Separating the two probabilities in seismic phase picking

There are two fundamental probabilities in the seismic phase picking process – the probability of the existence of a seismic phase (detection probability) and the probability of correctly identifying the phase arrival time (timing probability). The nearly ubiquitous approach in developing deep learning phase picking models is to use a kernel, such as a truncated Gaussian, to mask the labeled phase arrival time, and train a segmentation model. Once a model is trained, the times of the peaks in the output are taken as phase arrival times (picks) and the height of the peaks are taken as “probability” of the picks. Here, we show that this “probability” represents neither the detection nor the timing probabilty because this approach forces the output to follow the shape of the kernel. We introduce an approach using two models to estimate these two distinct probabilities. We use a binary classifier with a calibrated confidence to address the detection probability and a multi-class classifier to obtain a probability mass function to address the timing probability. This new approach makes the deep learning-based phase picking process more interpretable and gives us options to logically control seismic monitoring workflows.

Geophysical Journal International↗
USGS2025

Flood-Inundation Maps of the Current and Jacks Fork Rivers including the Ozark National Scenic Riverways, Southeast Missouri, 2023

Digital flood-inundation maps for a 131.8-mile reach of the Current River and a 44.6-mile reach of the Jacks Fork River, in southeast Missouri, were created by the U.S. Geological Survey (USGS) in cooperation with the Ozark Foothills Regional Planning Commission and the South Central Ozark Council of Governments. The maps also encompass the 134 miles of the Current and Jacks Fork Rivers within the Ozark National Scenic Riverways, which is the first national park area to protect a river system. The flood-inundation maps, which can be accessed through the USGS Flood Inundation Mapping Program website at https://water.usgs.gov/osw/flood_inundation/ , depict estimates of the areal extent and depth of flooding corresponding to selected water levels (hereafter referred to as “stages”) at eight reference USGS streamgages—five on the Current River (USGS station numbers 07064440, 07064533, 07066510, 07067000, and 07068000) and three on the Jacks Fork River (USGS station numbers 07065200, 07065495, and 07066000). Near-real-time stages at these streamgages may be obtained from the USGS National Water Information System at https://doi.org/10.5066/F7P55KJN or the National Weather Service National Water Prediction Service at http://water.noaa.gov/ , which also forecasts flood hydrographs at four of these sites (USGS station numbers 07067000, 07068000, 07065495, and 07066000). Flood profiles were computed for seven of the eight map reaches by means of two-dimensional hydraulic models and the remaining reach by a one-dimensional hydraulic model. The models were calibrated by using stage-streamflow relations or streamflow measurements at the USGS streamgages and from high-flow stage measurements from water-level loggers distributed throughout the reaches. The hydraulic models were used to compute water-surface profiles for flood stages at 1-foot intervals referenced to the streamgage datums. The profile stages ranged from the National Weather Service “action stage” or near bankfull, to a stage exceeding the highest recorded water level at each streamgage. The simulated water-surface profiles were then combined with a digital elevation model (derived from light detection and ranging data having a nonvegetated vertical accuracy of a maximum 10-centimeter root mean square error) to delineate the area flooded at each water level and the associated water depths. The availability of these maps, along with information regarding current stage from the USGS streamgage and forecasted high-flow stages from the National Weather Service, will provide emergency management personnel, resource managers, and residents with information that is critical for flood-response activities such as evacuations and road closures, as well as for postflood recovery efforts.

Missouri↗