USGS Science⌕ Search

SEARCH · USGS Science

Results for “Earth and Space Science”

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 49Linked to original sources
USGS2024

Paleoseismology of the Sawtooth fault and implications for fault behavior in the epicentral region of the 2020 Mw 6.5 Stanley, Idaho, earthquake

The 2020 moment magnitude ( M w ⁠⁠ ) 6.5 Stanley, Idaho, earthquake raised questions about the history and extent of complex faulting in the northwestern Centennial Tectonic Belt (CTB) and its relation to the Sawtooth normal fault and Eocene Trans‐Challis fault system (TCFS). To explore faulting in this area, we excavated a paleoseismic trench across the Sawtooth fault along the western margin of the CTB, and compared an early Holocene (9.1 ± 2.1 ka, 1δ ⁠ ) rupture at the site with lacustrine paleoseismic data and fault mapping in the 2020 epicentral region. We find: (1) a history of partial to full rupture of the Sawtooth fault ( M w ⁠ 6.8–7.4), (2) that shorter ruptures ( M w ≤ ⁠⁠6.9 ) are likely along distributed and discontinuous faults in the epicentral region, (3) that this complex system that hosted the 2020 earthquake is not directly linked to the Sawtooth fault, (4) that the northeast‐trending TCFS likely plays a role in controlling fault length and rupture continuity for adjacent faults, and (5) that parts of the TCFS may facilitate displacement transfer between normal faults that accommodate crustal extension and rotation. Our results help unravel complex faulting in the CTB and imply that relict structures can help inform regional seismic hazard assessments.

Idaho↗
USGS2024

Hydroacoustic observations reveal drivers of mixing and salinization of a karst subterranean estuary during intense precipitation

Karst subterranean estuaries within globally ubiquitous carbonate aquifers are coastal groundwater ecosystems that provide an essential water resource for human populations. To understand the drivers of salinization within a coastal aquifer in the Yucatan Peninsula (Mexico), we employed hydroacoustics in flooded caves to observe how oceanic and atmospheric events facilitate mixing between the meteoric lens (fresh-brackish groundwater) and the saline groundwater on tidal and episodic timescales. Precipitation during Tropical Storm Carlotta increased the flow and salinity of the meteoric lens without evidence for vertical mixing across the halocline. We postulate that vertical migration of haloclines in the conduit relative to those within the rock matrix during precipitation creates lateral density gradients that drive mixing, and ultimately creates a brackish layer within the meteoric lens. These results provide a mechanistic explanation for vertical and lateral exchange in a coastal carbonate aquifer, which has implications for groundwater response to future climatic change.

Ox Bel Ha cave system, Yucatan peninsula↗
USGS2024

Short-term sediment dispersal on a large retreating coastal river delta via 234Th and 7Be sediment geochronology: The Mississippi River Delta Front

Many Mississippi River Delta studies have shown recent declines in fluvial sediment load from the river and associated land loss. In contrast, recent sedimentary processes on the subaqueous delta are less documented. To help address this knowledge gap, multicores were collected offshore from the three main river outlets at water depths of 25–280 m in June 2017 just after the peak river discharge period, with locations selected based on 2017 U.S. Geological Survey seabed mapping. The coring locations included the undisturbed upper foreset, mudflow lobes, gullies, and the undisturbed prodelta. Nine multicores were analyzed for Beryllium-7 activity, and four cores were analyzed for excess Thorium-234 activity via gamma spectrometry, granulometry and X-radiography. Our results indicate a general trend of declining 7 Be and 234 Th activities and inventories with increasing distance from sources and in deeper water. The core X-radiographs are graded from the predominantly physically stratified nearshore to the more bioturbated offshore, consistent with the sedimentation patterns. Sediment focusing assessed via the 7 Be and 234 Th sediment inventories shows preferential sedimentation in gully and lobe environments, whereas the upper foreset and prodelta focusing factors are relatively depleted. Overall, short-term sediment deposition from the main fluvial source remains active offshore from all three major river outlets, despite the overall declining river load.

Louisiana↗
USGS2024

Peak streamflow trends in South Dakota and their relation to changes in climate, water years 1921–2020

Peak-flow (flood) frequency analysis is essential to water-resources management applications, including the design of critical infrastructure such as bridges and culverts, and floodplain mapping. Federal guidelines for performing peak-flow flood frequency analyses are presented in a U.S. Geological Survey Techniques and Methods Report known as Bulletin 17C. A basic assumption within Bulletin 17C, which documents the guidelines for determining annual peak streamflow frequency, is that, for basins without major hydrologic alterations (for example, regulation, diversion, and urbanization), statistical properties of the distribution of annual peak streamflows are stationary; that is, the mean, variance, and skew are constant through time. Nonstationarity is a statistical property of a peak-flow series such that the long-term (on the order of decades) distributional properties change one or more times either gradually or abruptly through time. Individual nonstationarities may be attributed to one source such as flow regulation, land-use change, or climate but are often the result of a combination of sources, making detection and attribution of nonstationarities challenging. In response to a growing concern regarding nonstationarity in peak streamflows in the region, the U.S. Geological Survey, in cooperation with the Departments of Transportation of Illinois, Iowa, Michigan, Minnesota, Missouri, South Dakota, and Wisconsin; the Montana Department of Natural Resources and Conservation; and the North Dakota Department of Water Resources, assessed the potential nonstationarity in peak streamflows in the north-central United States. This chapter characterizes the effects of natural hydroclimatic shifts and potential climate change on annual peak streamflows in the State of South Dakota. Annual peak and daily streamflow as well as model-simulated gridded climatic data were examined for temporal monotonic trends, change points, and other statistical properties indicative of changing climatic and environmental conditions. Changes in annual peak and daily flows were evaluated among 13, 35, and 81 qualifying U.S. Geological Survey streamgages for the 75-, 50-, and 30-year trend periods through water year 2020 (the period from October 1, 2019, to September 30, 2020) in South Dakota, respectively. No qualifying streamgages were in the 100-year trend period in the State. Statistical tests for autocorrelation (independent and identically distributed assumption), monotonic trends, and change points in the median and scale are analyzed to evaluate potential stationarity violations (nonstationarity) for performing at-site peak-flow flood-frequency analysis. The trends are reported using a likelihood approach as an alternative to simply reporting significant trends with an arbitrary p -value cutoff point. A distinct east-west spatial pattern of likely upward and downward monotonic trends and change points, respectively, was detected in 75- and 50-year trend periods, but an inconsistent spatial pattern was detected in the 30-year trend period. Additionally, change points in the median annual peak streamflows were detected in the late 1970s and early 1980s in the western part of the State, but in the east, the change point was more commonly detected in 1992–93. A similar east-west spatial pattern of likely upward and downward trends was detected in the annual peak-flow timing, the day of the year of the annal peak streamflow. In the western part of the State, the annual peak streamflows are arriving earlier, but in the east, the annual peak streamflows are arriving later. A peaks-over-threshold (POT) analysis where, on average, there are two events per year (POT2) and four events per year (POT4) was also used to evaluate changes in the frequency (count) of daily streamflows exceeding the threshold. Similar to detected changes in the annual peak streamflow, an east-west likely upward or downward change corresponding to an increase or decrease, respectively, in the frequency of daily streamflow greater than a POT2 and POT4 threshold was detected. A monthly water-balance model was used to evaluate hydroclimatic variation in annual and seasonal precipitation, snowfall, potential evapotranspiration, and soil moisture storage for all qualifying streamgages in the 75-, 50-, and 30-year trend periods. Detected trends in the annual hydroclimatic metrics for the 75- and 50-year trend periods indicate a spatially consistent statewide increase in precipitation, decrease in snowfall, increase in potential evapotranspiration, and increase in soil moisture storage. Furthermore, detected trends in seasonal precipitation in the 75- and 50-year trend periods highlight a pronounced change in precipitation in winter and later into the summer season, especially in the 50-year trend period in the eastern part of the State. Statewide increases in seasonal soil moisture storage were also detected, highlighting year-round increasing flood magnitudes, particularly in the eastern part of the State. Based on the results of these stationarity tests for the qualifying streamgages in South Dakota among the 75-, 50-, and 30-year trend periods, consistent temporal and spatial patterns of nonstationarity were detected among the 75- and 50-year trend periods. Furthermore, when nonstationarity is detected in daily streamflow, increased streamflow and volume (increasing frequency in POT), as well as potentially bridge scour, may have implications on culvert and highway design in the eastern part of South Dakota. Thus, when performing at-site peak-flow flood-frequency analyses in South Dakota, potential nonstationarities and alternative approaches are important considerations.

South Dakota↗
USGS2025

Capturing the uncertainty of seismicity observations in earthquake rate estimates: Implications for probabilistic seismic hazard analysis and the USGS National Seismic Hazard Model

The rate of earthquakes in a region is a fundamental input to Probabilistic Seismic Hazard Analysis. We present a Monte Carlo method for computing that rate from seismicity catalogs while including a range of data and analysis uncertainties. This method is applied to regions for which the b value is assumed to be spatially invariant. Each region is broken down into epochs for which each epoch is estimated to have a uniform magnitude of completeness ( ⁠⁠ M c ). The distribution of earthquake rates for M ≥ M c is determined for each epoch by considering the Poisson likelihood of rates given the number of observed earthquakes with M ≥ M c ⁠ . We use a Monte Carlo process to include the uncertainty in b , ⁠ M c , and individual event magnitudes. The result for each epoch is the joint distribution of the Poisson rate of earthquakes with magnitudes larger than the minimum value used to calculate hazard ( ⁠⁠ M 1 ) and the Gutenberg–Richter b values, which control the extrapolation to other magnitudes. The rate for each region is either the duration‐weighted average over the epochs or, to better capture temporal variations, we also consider mixture models. The mixture models also provide an avenue to allow temporal variations in b values. To implement this joint distribution in a logic tree, we use the mean and 95% confidence branches, each of which is parameterized with an M ≥ M 1 rate and b value. We explore different ways of defining those branches, as well as non‐Gutenberg–Richter branches, and their impact on hazard estimates. The mean hazard, but not the fractiles, is robust with respect to these choices. To illustrate these new methods, we use synthetic data and catalogs from recent U.S. Geological Survey National Seismic Hazard Models for the Central and Eastern United States and for Puerto Rico and the U.S. Virgin Islands.

central and eastern United States, Puerto Rico, U.↗
USGS2026

Storage timescales and the crystal growth history recorded in the layered reservoir of the Unit 3 Kulanaokuaiki Tephra Member of the Uēkahuna Ash from Kīlauea volcano

The Kulanaokuaiki Tephra Unit 3 (900 C.E.) was a subplinian basaltic scoria eruption, the largest known to have occurred at Kīlauea’s summit. The initiation mechanism and cause(s) of this highly hazardous eruption remain poorly understood, particularly as this style of activity has not been observed in the modern period. We use olivine diffusion timescales from Fe-Mg and lithium zoning to assess magma storage timescales and evaluate evidence for mixing prior to eruption. Of the 200 crystals assessed, 55% are normally zoned for forsterite and 45% are non-zoned. Zoned crystals (> 1 mm) ( n = 96) have average compositions of Fo 89 cores and Fo 81 rims, while small crystals 0.5–1 mm ( n = 14) have Fo 85 cores and Fo 82 rims. Modeled Fo diffusion timescales range from < 1 to 12.6 years (median of 3.0 years). Lithium maps in 32 of the forsterite-zoned crystals reveal a variety of lithium zoning types, most of which we attribute to crystal growth. Many lithium profiles are coupled to phosphorus within the outer ~ 100 μm of the crystal ( n = 17), reflecting growth processes and precluding their use for diffusion studies. Overall, our Unit 3 olivine dataset displays no evidence of magmatic priming shortly prior to eruption, unlike prior observations for the basaltic subplinian eruption of the Keanakāko‘i Tephra Unit E (1650 C.E.). Thus, explosive eruption magmas may be stored for considerably longer than previously thought, and magmatic priming is not a necessary precursor for Kīlauea’s explosive eruptions. We hypothesize that gas-driven filter pressing from crystallization of microlites and macrocrysts in the reservoir slurry contributed to the explosivity of the Unit 3 eruption.

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

Field observations and logs from the Rose Hip trench exposure across a north-facing scarp within the Seattle Fault Zone, southern Bainbridge Island, Washington

The Seattle Fault Zone is an approximately 70-km-long, east-west-trending zone of south-dipping blind reverse faults within the Puget lowland region in Washington. Because of the proximity, the Seattle Fault Zone poses a significant earthquake hazard to the Puget sound and Seattle metropolitan regions. We present preliminary mapping and trench-site information from a paleoseismic investigation across a newly identified active fault scarp located within the hanging wall of the Seattle Fault Zone on southern Bainbridge Island, Washington. The trench exposed monoclinally folded Miocene bedrock, fractured and faulted glacial-related deposits, and laminated lacustrine deposits capped by slope-derived colluvium. The observations from this investigation record late Pleistocene to Holocene north-vergent folding and faulting along this new fault scarp.

Washington↗
USGS2025

Conceptual and numerical groundwater flow model of the Iowa River alluvial aquifer near Tama County, Iowa, 1980 through 2022

The Iowa River alluvial aquifer is an important source of water on the Meskwaki Settlement in Tama County, Iowa, which is land owned by the Sac & Fox Tribe of the Mississippi in Iowa (commonly known as the Meskwaki Nation). The U.S. Geological Survey constructed a groundwater flow model, including a conceptual and numerical model, of the Iowa River alluvial aquifer and underlying hydrogeologic units near the Meskwaki Settlement in Tama County, Iowa, for the period of January 1980–August 2022 to estimate the fraction of water pumped from the Iowa River alluvial aquifer by Meskwaki Settlement wells that is derived from streamflow depletion in the Iowa River and its tributaries. Streamflow depletion is a reduction in streamflow caused by groundwater pumping and includes the interception by groundwater production wells of water that otherwise would have been discharged to streams (called “captured groundwater discharge”) and induced infiltration of streamflow to the production wells. Calibrated model runs were performed with no simulated pumping and simulated pumping only at Meskwaki Settlement wells, and the change in simulated flow rates between the groundwater system and streams for the two model runs represents the amount of streamflow depletion in the Iowa River and tributary streams resulting from pumping at the Meskwaki Settlement wells. Streamflow depletion in the Iowa River and its tributaries as a percentage of simulated pumping at the Meskwaki Settlement wells was calculated by dividing this difference by the total simulated pumping rate for the Meskwaki Settlement wells. The model results demonstrate that the mean monthly streamflow depletion, including induced infiltration and captured discharge, in the Iowa River and its tributary streams as a percentage of mean monthly pumping at the Meskwaki Settlement wells was 97.4 percent and ranged from 65.4 to 112 percent. Of the total streamflow depletion, mean monthly induced recharge was 20.9 percent and ranged from 4.9 to 37.2 percent. Mean monthly captured discharge was 76.5 percent and ranged from 57.1 to 97.1 percent. These results indicate that most of the water pumped from the Meskwaki Settlement wells is the result of streamflow depletion, in the form of both induced infiltration and captured discharge.

Iowa↗
USGS2025

Mechanics and statistics of postseismic shaking

Analysis of two weeks of continuous post-seismic shaking after the 2019 M7.1 Ridgecrest, CA earthquake sequence using 4 nearby borehole seismometers reveals that continuous ground motions decay as Omori’s law in time and follow the Gutenberg-Richter distribution in logarithmic amplitude. The measured temporal decay in amplitudes agrees with predictions of the rate-and-state framework and indicates shaking amplitudes are proportional to the velocity of afterslip. Our ground motion-based statistical framework provides a basis to forecast shaking intensity in the minutes to hours after a large earthquake.

California↗
USGS2024

Near-real-time earthquake-induced fatality estimation using crowdsourced data and few-shot large-language models

When a damaging earthquake occurs, immediate information about casualties (e.g., fatalities and injuries) is critical for time-sensitive decision-making by emergency response and aid agencies in the first hours and days. Systems such as the Prompt Assessment of Global Earthquakes for Response (PAGER) by the U.S. Geological Survey (USGS) were developed to provide a forecast of such impacts within about 30 min of any significant earthquake globally. However, existing disaster-induced human loss estimation systems often rely on early casualty reports manually retrieved from global traditional media, which are labor-intensive, time-consuming, and have significant time latencies. Recent approaches use keyword matching and topic modeling to identify human casualty-relevant information from social media but tend to be error-prone when dealing with complex semantics in multi-lingual text data and parsing dynamically changing and conflicting human death and injury numbers shared by various unvetted sources in social media platforms. In this work, we introduce an end-to-end framework to significantly improve the timeliness and accuracy of global earthquake-induced human loss forecasting using multi-lingual, crowdsourced social media. Our framework integrates (i) a hierarchical casualty extraction model built upon large language models, prompt design, and few-shot learning to retrieve quantitative human loss claims from social media, (ii) a physical constraint-aware, dynamic-truth discovery model that discovers the truthful human loss from massive noisy and potentially conflicting human loss claims, and (iii) a Bayesian updating loss projection model that dynamically updates the final loss estimation using discovered truths. We test the framework in real-time on a series of global earthquake events in 2021 and 2022 and show that our framework effectively automates the retrieval of casualty information faster but with comparable accuracy to those now retrieved manually by the USGS. The code associated with this work is made available at: https://github.com/SusuXu-s-Lab/Hierarchical-Earthquake-Casualty-Information-Retrieval

International Journal of Disaster Risk Reduction↗
USGS1976

Uranium and trace elements in stream sediments as an exploration tool

More than 45 trace elements have been reported in anomalous amounts in the uranium ores. The specific suite of elements associated with any one uranium deposit varies according to deposit type and geologic province. The primary geochemical halo of uranium and associated trace elements in the host rock, together with secondary dispersion halo in soils and alluvium, offers a potential geochemical exploration target. Sediment from streams near low-grade uranium occurrences in arkosic sandstone of the Denver basin, Colorado, and sediment from streams in the igneous and metamorphic terrace near the Midnite mine of eastern Washington have been analyzed for acid-extractable U, Cu,Pb, Zn, Ni, Co, Fe, Mn, Mo, Cd, Cr, Ag, V, and Se, and for total Fe, Mn, Hg, As, and organic C. The -80 mesh fraction of the alluvium downstream from a small uraniferous limonite occurrence in the Denver basin shows anomalous concentrations of U, Pb, As, Ni, Mo, V, Zn, and Se. Anomalous concentrations of U persist farther downstream (as much as 1,000 m) than any of the trace elements. In another stream system southeast of Denver, above-background levels of U, Pb, As, Zn, V, Fe, and Mn occur near a uranium anomaly (defined by other methods) which has no surface mineralization. In the Midnite mine area, sediments of streams draining prospect areas showed anomalous concentrations of uranium (as much as 10 times background), but no accompanying anomalous trace-element concentrations, although analyses of Midnite mine ores show that a broad suite of trace elements are associated with uranium.

Open-File Report↗
USGS2025

Geologic map of the Guinevere Planitia quadrangle (V–30), Venus

This 1:5,000,000-scale geologic map of the Guinevere Planitia quadrangle divides the region into 15 geologic material units, defined using Magellan synthetic aperture radar (12.6-centimeter-wavelength radar system; 75 meters per pixel) datasets and including upland terrain units (2.4 percent of the surface area), plains materials units (59 percent), flow materials associated with named and unnamed eruptive centers (37.2 percent), small volcanic edifices, and impact crater materials (1.4 percent). Upland terrain units consist of tessera and lineated upland material, plains materials consist of Guinevere regional plains and Guinevere lineated and mottled plains, and flow materials consist of lobate flow material and plains-forming flow material. Specific lobate flows associated with Atanua Mons, Tuli Mons, Var Mons, and Uilata Fluctus are mapped separately. Other mapped units are impact crater material and small volcanic edifice. In addition to geologic units, we mapped linear features that show patterns of deformation or flow across the quadrangle. These consist of faults, wrinkle ridges, broad arches, channels, troughs, and flow direction indicators. The map region also contains several small volcanic features: shields, depressions, and craters. These, in combination with the plains, large volcanoes, and coronae, show the pervasive influence of volcanism across Venusian lowlands. The rims of nine identified impact features are delineated; large bright and dark haloes, which in some cases are associated with individual impact craters, are mapped as surficial mantling deposits. We documented spatial relationships using the stratigraphic and cross-cutting relationships of the quadrangle’s geologic units and features to provide a synthesis of the region’s geologic history. The upland terrain of the quadrangle indicates intense tectonic deformation and uplift. It is exposed as embayed remnants, typically within the plains, and represents the oldest geologic materials locally and across the region. Guinevere plains and the plains-forming flow unit appear to be assemblages of volcanic flows from multiple sources, including distinct coronae and corona-like structures. The temporal evolution of Guinevere lineated and mottled plains was likely protracted, with continued formation of small volcanic edifices over a long period. The morphologic and radar brightness characteristics of volcanoes in the region indicate their growth may have involved (1) multiple large-scale eruptive centers with recognizable spatial and temporal sequences, (2) extensive lava flow fields with a multitude of flows producing complex, overlapping patterns, and (3) numerous small volcanic edifices, including shields, domes, and cones. Although geologic patterns common to other regions of Venus are evident in the Guinevere Planitia quadrangle, local relative age relationships are inconsistent or unclear, preventing robust stratigraphic correlation. The mapping results do, however, indicate complicated local sequences of volcanic and tectonic activity.

Scientific Investigations Map↗
USGS2024

The 3D National Topography Model Call for Action—Part 1. The 3D Hydrography Program

The U.S. Geological Survey is initiating the 3D Hydrography Program (3DHP), the first systematic remapping of the Nation’s surface waters since the original 1:24,000-scale topographic mapping program was active from 1947 to 1992. Building on decades of experience maintaining the National Hydrography Dataset (NHD), the Watershed Boundary Dataset (WBD), and the NHDPlus High Resolution (NHDPlus HR), the 3DHP will completely refresh the Nation’s hydrography data and improve discovery and sharing of water-related data. The design of the 3DHP is based on the results of a study that estimated that the fully implemented program would have the potential to provide more than $1 billion in benefits to Federal, State, Tribal, Territorial, and local governments and to private and nonprofit organizations every year, in addition to myriad societal benefits. The 3DHP would directly support better decision making regarding water resources by providing more accurate, complete, and integrated information than is currently available. The 3DHP datasets will include a three-dimensional (3D) hydrography network generated from and integrated with elevation data from the 3D Elevation Program (3DEP) to better represent stream gradients and channel conditions, along with waterbodies, hydrologic units, hydrologically enhanced elevation and other surfaces, and more consistent and accurate attributes. The 3DHP datasets will inherit key attributes of the NHD, WBD, and NHDPlus HR, and they also will include new attributes and links to other data such as the U.S. Fish and Wildlife Service National Wetlands Inventory, groundwater data, and engineered hydrologic systems such as stormwater networks. The 3DHP will be designed to provide a set of open and interoperable web-based tools, maps, and data catalogs, creating a robust system for users to reference their information about water; the system elements are collectively referred to as the “infostructure.” The 3DHP and the infostructure can provide a foundational geospatial underpinning for the Internet of Water, a community-based effort to modernize tools and technologies to share water data. As proposed, the 3DHP would begin providing products and services to the public in 2024.

Circular↗
USGS2000

Use of a ground-penetrating radar system to detect pre- and post-flood scour at selected bridge sites in New Hampshire, 1996-98

Ground-penetrating radar was used to measure the depth and extent of existing and infilled scour holes and previous scour surfaces at seven bridges in New Hampshire from April 1996 to November 1998. Ground-penetrating-radar survey techniques initially were used by the U.S. Geological Survey to study streambed scour at 30 bridges. Sixteen of the 30 bridges were re-surveyed where floods exceeded a 2-year recurrence interval. A 300-megahertz signal was used in the ground-penetrating radar system that penetrated through depths as great as 20 feet of water and as great as 32 feet of streambed materials. Existing scour-hole dimensions, infilled thickness, previous scour surfaces, and streambed materials were detected using ground-penetrating radar. Depths to riprap materials and pier footings were identified and verified with bridge plans. Post data-collection-processing techniques were applied to assist in the interpretation of the data, and the processed data were displayed and printed as line plots. Processing included distance normalization, migration, and filtering but processing was kept to a minimum and some interference from multiple reflections was left in the record. Of the 16 post-flood bridges, 22 ground-penetrating-radar cross sections at 7 bridges were compared and presented in this report. Existing scour holes were detected during 1996 (pre-flood) data collection in nine cross sections where scour depths ranged from 1 to 3 feet. New scour holes were detected during 1998 (post-flood) data collection in four cross sections where scour depths were as great as 4 feet deep. Infilled scour holes were detected in seven cross sections, where depths of infilling ranged from less than 1 to 4 feet. Depth of infilling by means of steel rod and hammer was difficult to verify in the field because of cobble and boulder streambeds or deep water. Previous scour surfaces in streambed materials were identified in 15 cross sections and the depths to these surfaces ranged from 1 to 10 feet below the streambed. Riprap materials or pier footings were identified in all cross sections. Calculated record depths generally agree with bridge plans. Pier footings were exposed at two bridges and steel pile was exposed at one bridge. Exposures were verified by field observations.

New Hampshire↗
USGS2026

Metallurgical coal—Deposits, production, resources, market dynamics, and supply chain risks

Plain Language Summary Metallurgical coal (met coal; consumed to produce coke for steelmaking) must meet specific chemical and physical specifications. In 2023, the conterminous United States produced 66 million short tons (mst) of met coal, consumed 15.85 mst domestically, exported 51.1 mst, and imported 0.7 mst. Most met coal was produced in the Appalachian Basin, but there are also resources that meet the specifications for met coal in the Western United States.

conterminous United States↗
USGS2026

Utilizing hydrophones to detect streambed mobilization in the Wild and Scenic reach of the Rio Chama

This paper describes a study that deployed underwater acoustic sensors (hydrophones) to detect streambed mobilization in the Wild and Scenic reach of the Rio Chama, New Mexico. Fine sediment accumulation in the gravel‑cobble bed affects brown trout spawning habitat, and understanding incipient motion thresholds is essential for planning effective environmental flushing flows released from El Vado Dam. Two hydrophone stations were operated during 2021–2022, capturing audio data during a high‑flow pulse in late 2021. Acoustic processing techniques, including spectral analysis and threshold‑based impact detection, were used to identify sediment‑generated noise associated with coarse sediment transport. Results indicate detectable sediment mobilization during the high‑flow pulse, demonstrating hydrophones’ potential as a passive monitoring tool for environmental flow planning and sediment‑management decision support.

Conference Paper↗
USGS2025

Subduction zone earthquake catalog separation tool: Implementation in the USGS 2025 Puerto Rico and U.S. Virgin Islands National Seismic Hazard Model

The U.S. Geological Survey (USGS) periodically releases updates to National Seismic Hazard Model (NSHM) for the United States and its territories leveraging current scientific knowledge and methodologies to guide public policy, building codes, and risk assessments regarding potential ground shaking due to earthquakes that may result in infrastructure damage. In subduction zones, there is a need to separate the earthquake catalog into tectonic regimes to create specific seismicity models for which the most appropriate ground‐motion models are then applied. Here, we describe newly developed methods and software, called CatSep, that classifies subduction zone events into three primary tectonic regimes: crustal, interface, and intraslab. This method incorporates information about the location of the earthquake relative to the subducting slab, the depth of the Mohorovičić discontinuity, and the earthquake’s moment tensor. Applying this method is a first step in the NSHM workflow for regions covering U.S. subduction zones. Results using this subduction zone earthquake catalog separation tool for the 2025 Puerto Rico and U.S. Virgin Islands NSHM earthquake catalog are presented and analyzed.

Puerto Rico, U.S. Virgin Islands↗