Search USGS⌕ Search

SEARCH · Search USGS

Results for “Journal of Hydrology, Regional Studies”

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,664 records · Page 5Linked to original sources
USGS2024

Evaluating the impact of uncertainty in ground motion forecasts for post-earthquake impact modeling applications

The US Geological Survey’s (USGS) ShakeMap system provides a rapid characterization of strong ground shaking in areas directly affected by an earthquake. This study focuses on studying the aggregate effects of macroseismic shaking estimates from ShakeMap, expressed in terms of modified Mercalli intensity (MMI), when accounting for the uncertainty in forecasted ground motions. We use a Monte Carlo approach to generate numerous spatially correlated realizations of ground motions by utilizing a combination of circulant embedding and kriging techniques for efficiently handling the correlations. We then assessed the aggregate effects of shaking by looking at bin counts across these realizations. We demonstrate that the aggregate shaking regarding the mean macroseismic intensity estimates (from the ShakeMap output) is a biased representation of the aggregate shaking when shaking uncertainty is included. Incorporating shaking uncertainty can help to improve various downstream earthquake impact applications, such as the USGS Prompt Assessment of Global Earthquakes for Response (PAGER) overall earthquake fatality distribution or estimates of shaking-induced ground failure impacts from consequential earthquakes.

Earthquake Spectra Journal↗
USGS2026

Mineral chemistry perspective on remobilization of stored magma at Kamakai'a Hills, Southwest Rift Zone of Kilauea, Island of Hawai'i, USA

Differentiated magmas stored in the rift zones of Kīlauea have received more attention in recent years following eruption of andesite during the early phase of 2018 lower East Rift Zone activity. Despite this growing interest, some of the most voluminous eruptions of differentiated rift zone magmas remain poorly studied. One such eruption, and the most voluminous exposed differentiated flow field at Kīlauea, is the Kamakaiʻa Hills. This eruption took place in the Southwest Rift Zone of Kīlauea, a region that is hypothesized to contain a long-lived rift zone reservoir. The Kamakaiʻa Hills flow field encompasses >250 × 10 6 m 3 of basaltic andesite and basalt compositions with a mineral assemblage of orthopyroxene + clinopyroxene + plagioclase during its early ʻaʻā phase and clinopyroxene + plagioclase + olivine during its late pāhoehoe phase. To better understand storage conditions and magma accumulation, this study focuses on major, minor, and trace elements from the mineral assemblage present within the early ʻaʻā and late pāhoehoe phases. The diversity of clinopyroxene and plagioclase compositions within the early ʻaʻā and late pāhoehoe phases, as well as diverse compositions of plagioclase and orthopyroxene within the early ʻaʻā phase, suggest multiple magma bodies and limited pre-eruption magma mixing within the broader Kamakaiʻa Hills reservoir. Oscillatory zoning patterns (particularly in clinopyroxene) imply processes such as recharge events, magma mixing or mingling, or convection within a differentially cooling, chemically stratified reservoir over protracted time intervals, whereas only limited resorbed mineral textures indicate incomplete mixing of heat and chemically distinct magmas during the dike intrusion that triggered the eruption. Mineral-mineral and mineral-melt thermobarometry indicate predominantly shallow (≤2.5 km depth) crustal storage conditions of the cooled, differentiated magma (∼1100 °C and cooler for the basaltic andesites) to hotter temperatures for the basalts (all >1100 °C). Despite the known large standard errors estimated for mineral-melt and mineral-mineral barometry (10s to >100 MPa), the calculated pressures and depths broadly correspond with earthquake swarm depths beneath the Kamakaiʻa Hills, and drill core and fluid inclusion barometry storage depths of differentiated magmas within the lower East Rift Zone. The Kamakaiʻa Hills differentiated magmas have H 2 O contents (∼0.5 wt%, using plagioclase-melt hygrometry) equivalent to typical Kīlauea basalts. Our data and interpretations demonstrate a complex, long-lived rift zone storage system that consisted of multiple magma bodies and was mobilized into eruption through intrusion of a hotter and more primitive summit-derived (uprift) magma.

Hawaii↗
USGS2026

Geophysical architecture and geochronology of the Neoarchean Mentor anorthosite intrusive complex, northwestern Minnesota: Largest anorthosite complex of the Superior Province?

The Neoarchean Mentor anorthosite intrusive complex (MAIC) lies within the Wawa subprovince in northwestern Minnesota, in a region where the Wawa, Quetico, and Wabigoon subprovinces are juxtaposed in close proximity. Archean rocks are entirely concealed, and interpretations are developed from aeromagnetic, gravity, and borehole samples. The MAIC includes both anorthosite (dense, weakly magnetized) and oxide-rich gabbro (strongly magnetized) lithologies. Anorthosite is coarse-grained to megacrystic, intensely altered, and locally brecciated. Pervasive epidote alteration enhanced the density of the anorthosite via introduction of higher-density mineral assemblages, explaining why the MAIC produces a significant gravity high. Oxide-rich gabbro forms a border phase of the MAIC and has potential for vanadium, chromium, and titanium mineralization, and produces a strong aeromagnetic high. The MAIC is interpreted to extend over an area of 640 km 2 , making it the largest known anorthosite complex of the Superior Province, as measured by preserved areal extent. Modeling indicates the MAIC extends more than 6 km into the subsurface. A new Pb–Pb zircon age of 2737.2 ± 4.5 Ma is interpretated as the crystallization age of anorthosite within the MAIC, showing that the MAIC formed well before the ca. 2690 Ma Shebandowanian orogeny, and raising new questions about correlations with other parts of the Wawa subprovince. A low-density batholith, here informally called the Fertile batholith, is interpreted to intrude the southern part of the MAIC. A new Pb–Pb zircon age of 2701.1 ± 6 Ma is interpreted as the magmatic age of the Fertile batholith.

Minnesota↗
USGS2025

New constraints on location and timing of the Great Lakes tectonic zone, central Upper Peninsula, Michigan, USA

The Great Lakes tectonic zone (GLTZ) forms the boundary between the Wawa–Abitibi and Minnesota River Valley subprovinces within the Archean Superior Province. The GLTZ is concealed for all of its 1100 km length, except for a segment in the central Upper Peninsula of Michigan. There, it is exposed as a northwest-striking mylonite zone along a 11 km segment, extending to the onlap of Paleozoic rocks to the east. Farther east, its location has been unknown. Here, we use aeromagnetic and gravity data to develop interpretations of the expression of the GLTZ and to define its extent under cover. Aeromagnetic gradients over the mylonite zone are interpreted to be produced by structurally juxtaposed rocks with varying magnetizations. Gravity data show a regional gradient along the GLTZ, produced by the juxtaposition of a dense greenstone belt on the north against lower-density gneisses and granites on the south. The GLTZ is interpreted to extend ∼55 km under cover to the east. The GLTZ is terminated on the east by the buried eastern arm of the ca. 1100 Ma Midcontinent Rift. An undeformed granitic dike that cuts the mylonitic foliation produces a U–Pb apatite age of 2523 ± 33 Ma, implying no major post-Archean shearing occurred, and is at odds with previous interpretations of major Proterozoic reactivation. A granite intrusion in the Minnesota River Valley subprovince produces a Pb–Pb zircon age of 2606.9 ± 3.6/7.4 Ma. This suggests that magmatism related to the Sacred Heart orogeny, previously known in Minnesota, extended to Michigan.

Michigan↗
USGS1996

Technique for estimating magnitude and frequency of peak flows in Maryland

A convenient and reliable technique for estimating flood magnitudes is required for effective flood-plain management and for the efficient design of bridges, culverts, embankments, and flood-protection structures. Methods are presented for estimating peak-flow magnitudes of selected frequencies, ranging from 2 to 500 years, for all nontidal drainage basins in Maryland. The methods were developed by generalized least-squares regression techniques using data from 219 gaged basins in and near Maryland. The State is divided into five hydrologic regions: the Appalachian Plateaus and Allegheny Ridges region, the Blue Ridge and Great Valley region, the Piedmont region, the Western Coastal Plain region, and the Eastern Coastal Plain region. These regions correspond to the physiographic provinces of the State, with the exceptions that (1) the Coastal Plain Province is divided into two hydrologic regions, and (2) there is no distinct hydrologic region corresponding to the Valley and Ridge Province as it is divided into its constituent Allegheny Ridges and Great Valley subdivisions. Sets of equations for calculating peak discharges based on physical basin characteristics are provided for each of the regions. Based on the peak-flow equations, methods for estimating peak flows are presented for ungaged and gaged streams in Maryland. The methods and equations are supported by generalized least-squares analysis of basin and flood-frequency characteristics data from 219 drainage basins in and near Maryland. estimates for each of the five regions are calculated using combinations of the fol-lowing basin characteristics: drainage area, forest cover, basin relief, carbonate rock coverage, storage, and runoff-curve number. Drainage area contributes to the estimate in all five study regions. Carbonate rock coverage is used only in the Blue Ridge and Great Valley region. Storage and runoff-curve number are used solely in the Eastern Coastal Plain region. All other basin characteristics are used in two or more regions. Standard errors of estimate for the regression equations range from 19 to 31 percent in the Appalachian Plateaus and Allegheny Ridges region, 34 to 47 percent in the Blue Ridge and Great Valley region, 33 to 48 percent in the Piedmont region, 45 to 64 percent in the Western Coastal Plain region, and 36 to 42 percent in the Eastern Coastal Plain region.

Maryland↗
USGS2026

Seismic characteristics of the April 2024 eruption of Ruang Volcano, North Sulawesi, Indonesia

The two most recent eruptions of Ruang volcano, in 2002 and 2024, have been large (VEI 4) with very short precursory periods, 3 and 6 days, respectively. Background seismicity at Ruang volcano is generally low with 0–2 volcano-tectonic (VT) earthquakes per day. The 2024 eruptions were preceded by a notable increase in VT earthquakes from 11 to 15 April with a rapid escalation from 16 to 17 April in the form of increased VT swarms, and the appearance of low frequency earthquakes and short bursts of VT drumbeats shortly before the eruption. Similar VT seismic unrest was recorded in March 2015, May 2016, and April 2022, although none of these episodes culminated in an eruption. A comparison of seismic energy release between April 2022 and April 2024 shows distinct differences. In 2024, cumulative seismic energy reached its peak within 5 days, just before the onset of eruption on 17 April. In contrast, during the 2022-episode, energy release followed a more gradual increase over 8 days peaking on 16 April, then declining gradually without leading to eruption. Notably, low frequency (LF) and drumbeat earthquakes were absent in non-eruptive unrest episodes. The transition to proximal LFs in 2024 likely signals pressure fluctuations, intensified magma degassing, and conduit formation, acting as an important precursor to eruption. Additionally, the occurrence of short bursts of drumbeat seismicity likely reflects the onset of the shallow magma ascent and was an important short-term precursor to the eruption. Finally, we examine the potential influence of regional tectonic earthquakes on the escalation of precursory unrest in 2024.

North Sulawesi, Ruang volcano↗
USGS2025

Gas emissions from the Sulphur Bank Mercury Mine hydrothermal system, Clear Lake volcanic field, California

The Sulphur Bank Mercury Mine (SBMM) hydrothermal system offers insights into active degassing processes in the Clear Lake volcanic field (CLVF), a high-threat region based on its record of Holocene eruptions and proximity to populated areas. Here we present chemical and isotopic analyses of gas samples collected between 2015 and 2023, along with the first comprehensive CO 2 flux survey of the SBMM area conducted in 2023. Sampled gases are CO 2 - and CH 4 -rich (≥84 and 6 mol% in dry gas, respectively) with high mantle-derived helium contributions ( 3 He/ 4 He = 6.54–7.86 R C /R A ). Carbon isotopic compositions of CO 2 (δ 13 C = −10.0 to −9.5 ‰) and CH 4 (δ 13 C = −35.8 ‰) indicate mixed sources, with significant contributions from metamorphism of organic-rich Franciscan Complex rocks hosting the hydrothermal system. Modeling of gas compositions shows that scrubbing by interaction with air-saturated groundwater strongly influences observed compositional variability. From our CO₂ flux measurements, we estimate the deeply derived CO 2 emission rate from the SBMM hydrothermal area (0.2 km 2 ) at 240 t d −1 , comparable to many quiescently degassing volcanoes worldwide. We also provide a first-order estimate of CH 4 emissions at approximately 0.5 t d −1 . Our findings establish crucial baseline data for future volcanic monitoring efforts, enhancing detection capabilities for potential changes in this active hydrothermal system. This work contributes to the broader understanding of volatile contributions from volcanic and metamorphic sources to the global carbon budget, while highlighting the strong influence of bedrock geology on gas compositions in the CLVF.

California↗
USGS2025

VIPER site analysis

We needed to evaluate available orbital data of NASA’s Volatiles Investigating Polar Exploration Rover (VIPER) mission area in order to derive a variety of maps to help the science team identify scientifically interesting places for the rover to visit and to provide scientific context for our mission. Some of these maps also fulfilled engineering and mission design needs to enable safe and efficient landing and roving. We incorporated data from the Lunar Reconnaissance Orbiter Camera, the Lunar Orbital Laser Altimeter, the Mini-RF instrument, the Chandrayaan-2 Orbital High Resolution Camera, the Korean Pathfinder Lunar Orbiter’s Shadowcam, the Kaguya Spectral Profiler and Multiband Imager, and the Chandrayaan-1 Moon Mineralogy Mapper. We used a variety of techniques to build these maps, including stereogrammetry, shape-from-shading, ice stability depth and surface temperature calculations, and the horizon method for solar illumination and direct-to-Earth communications maps. Altogether, these maps allowed us to survey for boulders, evaluate features in permanently shadowed regions that VIPER might explore, provide mineralogic context for what VIPER’s instruments may learn, estimate the ages and radar properties of craters in the VIPER mission area, and evaluate the potential for gravity traverses with the rover. These data and techniques provided a rich set of information from which both the VIPER science team and engineering teams were able to draw in order to plan a safe landing and to plan a VIPER surface mission that will be both scientifically valuable and robust from an operational perspective.

The Planetary Science Journal↗
USGS2026

Detecting snow avalanche activity using infrasound: Hooker Valley, New Zealand

Snow avalanches pose considerable hazards to people and infrastructure in alpine environments. Traditional avalanche monitoring relies on meteorological data and visual observations, which can be limited in scope and timeliness. Infrasound offers a promising complementary monitoring tool by detecting the low-frequency sound waves generated by avalanches. Here, we present infrasound and camera observations during a 50-day field campaign in the Hooker Valley of Aoraki/Mount Cook National Park, New Zealand. Our study detected seven avalanches with the cameras, whereas the infrasound system identified only one of these events, which was the largest and occurred under conditions that likely favoured infrasound propagation. The infrasound system recorded numerous other events not captured by the cameras, indicating the benefit of further investigation to determine their sources. These findings highlight the potential of infrasound technology for detecting avalanches and providing broad spatial coverage, capturing events in areas not monitored by cameras, while also showcasing limitations in infrasound capabilities. The limited detection of smaller avalanches underscores the opportunity for further research to enhance detection capabilities and understand environmental influences such as snow cover and wind noise. Overall, this study emphasises the utility of multidisciplinary monitoring techniques to improve avalanche detection in alpine environments.

Hooker Valley↗
USGS1998

Peak-discharge frequency and potential extreme peak discharge for natural streams in the Brazos River basin, Texas

The 2-, 5-, 10-, 25-, 50-, and 100-year peak discharges were estimated for 186 streamflow-gaging stations with at least 8 years of data for natural streams in and near the Brazos River Basin, Texas. Multiple regression equations were developed to estimate peak-discharge frequency for the 2-, 5-, 10-, 25-, 50-, and 100-year recurrence intervals for each of three hydrologic regions that compose the Brazos River Basin. The equations for each region are a function of significant basin characteristics (explanatory variables). The significant explanatory variables among six that were tested are the contributing drainage area and stream slope for regions 1 and 2 and the contributing drainage area for region 3. For the three sets of equations, the coefficient of determination ranges from 0.59 to 0.93, and the standard error ranges from 0.184 to 0.391 log units. A larger coefficient of determination and a lower standard error generally are associated with the equations for hydrologic regions 2 and 3. Statistics from the regression analysis allow computation of the prediction interval associated with a given significance level for a peak-discharge frequency estimate. The regression equations can be used to estimate peak discharges for sites at, near, or away from sites with streamflow-gaging stations. The potential extreme peak-discharge curves as related to contributing drainage area were estimated for each of the three hydrologic regions from measured extreme peaks of record at 186 sites with streamflow-gaging stations and from measured extreme peaks at 37 sites without streamflow-gaging stations in and near the Brazos River Basin. The potential extreme peak-discharge curves generally are similar for hydrologic regions 1 and 2, and the curve for region 3 consistently is below the curves for regions 1 and 2, which indicates smaller peak discharges.

Texas↗
USGS2026

Assessment of extreme subsurface hydrologic conditions captured during atmospheric river storms in the San Francisco Bay area (California, USA) with applications to shallow landslide early warning

An increase in soil pore water pressure is the typical trigger for the majority of landslides caused by rainfall. Atmospheric river storms, common to the west coast of North America during the winter season, can deliver landslide triggering rainfall resulting in severe impacts to coastal communities. Using a network of hydrologic monitoring stations situated within landslide-prone terrain in the San Francisco Bay area of California (USA), we assess the meteorologic conditions and resulting hydrologic and landslide response resulting from eight consecutive storm events that caused thousands of shallow landslides during the winter of 2022–2023. We find disparate hydrological responses and resultant degrees of observed landsliding ranging from < 1 landslide/km2 to 18 landslides/km2 at the monitoring sites that reflect the interplay and differences between rainfall delivery, subsurface hydrological characteristics, and geotechnical properties at each site. Antecedent soil moisture from both early season rainfall and the first storm in the sequence played a critical role in setting up some hillslopes for failure. Subsequent storms then generated elevated pore water pressures for several hours with associated landsliding. However, we find that the occurrence of widespread landsliding required not only sufficient pore water pressure magnitude in susceptible hillslopes, but also full and prolonged effective soil saturation throughout hillslope profiles. Landslides may still occur at lower values and durations of effective saturation but are likely to be less extensive regionally. We present these findings within the context of research directions and improvements to landslide early warning systems first suggested by researchers 40 years ago.

California↗
USGS2024

Earthquake cycle mechanics during caldera collapse: Simulating the 2018 Kīlauea eruption

In multiple observed caldera-forming eruptions, the rock overlying a draining magma reservoir dropped downward along ring faults in sequences of discrete collapse earthquakes. These sequences are analogous to tectonic earthquake cycles and provide opportunities to examine fault mechanics and collapse eruption dynamics over multiple events. Collapse earthquake cycles have been studied with zero-dimensional slider-block models, but these do not account for the complicated interplay between fluid and elastic dynamics or for factors such as the heterogeneous fault properties and non-vertical ring fault geometries often inferred at volcanoes. We present two-dimensional axisymmetric mafic piston-like collapse earthquake cycle models that include rate-and-state friction, fully-dynamic elasticity, and compressible viscous fluid magma flow. We demonstrate that collapse earthquake intervals and magnitudes are highly sensitive to inertial effects, evolving stress fields, fault geometry, and depth-varying fault friction. Given the consistent earthquake cycles observed in most eruptions, this suggests that ring faults can quickly stabilize and often become nearly vertical at depth. We use the well-monitored 2018 collapse sequence at Kı̄lauea as a case study. Our model can produce many features of Kı̄lauea seismic and geodetic observations, except for a significant amount of interseismic slip, which cannot be readily explained with simple rate-and-state friction parameterizations.

Hawaii↗
USGS2026

Synthesizing a twelve-year sediment trap time series of planktic foraminiferal flux in the Gulf of America (Mexico)

Sediment trap time series provide powerful frameworks for testing hypotheses about planktic foraminiferal assemblage composition, seasonality, and geochemical responses to environmental variability, all of which are central to improving paleoceanographic reconstructions. We present high resolution foraminiferal assemblage data from a long-running sediment trap (2008–2020) in the northern Gulf of America (Mexico). This study summarizes the species composition, seasonality, and size distribution of the fifteen most abundant species of planktic foraminifera, which account for 98% of total flux. Foraminiferal flux peaks in winter and reaches a minimum in summer, following the seasonal pattern of primary production in the northern Gulf. Winter assemblages are dominated by non-spinose taxa, whereas spinose taxa prevail during summer. Across nearly all species, average monthly test size covaries with temperature, independent of seasonal flux trends. Notably, Trilobatus sacculifer and Neogloboquadrina dutertrei show a significant decline in relative abundance and flux during 2017–2020, nearly disappearing from the assemblage.

Journal of Foraminiferal Research↗
USGS2025

National population exposure and evacuation potential in the United States to earthquake-generated tsunami threats

Previous efforts to characterize tsunami threats to people have focused primarily on individual scenarios in specific areas but have not recognized multiple scenarios across an entire country. This study addresses this gap by quantifying population exposure and evacuation potential in the United States to 102 earthquake-related, tsunami-hazard zones, including 92 local scenarios, 8 distant scenarios, and 2 probabilistic products. Geospatial path-distance modeling quantified evacuation potential and the influence of departure delays. We focused on residents to support other national, multi-hazard risk analyses. Millions of residents are in distant-tsunami zones, and hundreds of thousands of residents are in local-tsunami zones. In 41 scenarios, there is at least one resident that may have insufficient time to evacuate before wave arrival. Tens of thousands of residents may have insufficient time to evacuate from local tsunamis that impact the U.S. Pacific Northwest or Puerto Rican coastlines. The largest improvements in evacuation potential may come from reducing departure delays in some areas but may involve vertical-evacuation structures or changing land use in other areas.

International Journal of Disaster Risk Reduction↗