USGS Science⌕ Search

SEARCH · USGS Science

Results for “Frontiers in Veterinary Sciences”

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,605 records · Page 56Linked to original sources
USGS2026

Tephra from Kīlauea’s 2008–2018 lava lake eruption—Proximal deposits and dispersal characteristics

A network of ten buckets was established early in the 2008–2018 summit eruption at Kīlauea to collect proximal tephra ejected from the new, informally named the “Overlook crater”; the buckets were emptied on most days of the eruption thereafter. This report summarizes the results of more than 2,400 different sampling intervals (most 1–3 days long) during the eruption, focusing on the physical and dispersal characteristics of the tephra deposits. The network was within about 300 meters south of the vent to capture tephra dispersed by the dominant northeast trade wind. The juvenile tephra mainly reflected spattering at the southeast (SE) sink, a downwelling area in the southeastern part of the lava lake in the Overlook crater that remained in the same area throughout the eruption, with admixtures of solid rock and secondary minerals derived from the wall of the crater. The proportion of juvenile material to lithic material ranged widely early in the eruption but was generally greater than 90 percent for the last 6 years of the eruption as lake level rose and the crater walls decreased in height and became more stable. The accumulation rate of tephra at each bucket was strongly dependent on the location of the bucket and reflects the interplay between lava lake level and wind direction and speed. The mass per unit area (m/a) of collected tephra was a maximum of about 97 kilograms per square meter for the entire eruption, equivalent to a thickness of about 75 millimeters. Thirty-two explosive events with a volcanic explosivity index of −2 to −4 deposited much of the tephra in the network. Fifteen of these rock-fall-induced events occurred within 6 days of one another. Pele’s hair, a volcanic glass predominantly associated with quieter activity at the vent in the final half of the eruption, was dispersed more than 60 kilometers downwind from the lava lake and formed a nearly continuous deposit near the Overlook crater. This eruption was probably the most frequently sampled long-lasting eruption in history, but most of the deposits are ephemeral. The collected samples, although generally of small mass, are retained by the Hawaiian Volcano Observatory and are available for detailed study.

Hawaii↗
USGS2024

Perspectives on transportable array Alaska background noise levels

Background seismic noise fundamentally sets a lower bound on our ability to record signals arising from earthquakes. The background noise spectrum at a station is a combination of cultural noise, ocean-generated microseism noise, intrinsic instrument self-noise, and the sensitivity of the instrument to nonseismic noise sources. The USArray-Transportable Array Alaska deployed 195 stations across Alaska and parts of Canada (Yukon, British Columbia, and Northwest Territories). These stations were all installed using similar techniques and made use of instruments with similar self-noise levels. As such, this network provides an opportunity to look at how geographic location influences seismic background. Using these broadband stations, we report background noise levels from 0.2 to 75 s period in six discrete bands. By constructing “noise maps,” we depict both spatial and temporal changes in the background noise field. Using these maps, combined with targeted analysis, we infer sources and contributing factors to noise levels in these different period bands. These include cultural noise, the formation of sea ice, seasonal changes in permafrost and wave activity in the Gulf of Alaska, and magnetic field variability. We use this study as an opportunity to review several previous studies examining seismic noise in Arctic regions.

Book chapter↗
USGS1992

Nonpoint-source pollutant discharges of the three major tributaries to Reelfoot Lake, west Tennessee, October 1987-September 1989

An investigation of the concentration and loads of nitrogen, phosphorus, and suspended sediment in storm runoff to Reelfoot Lake, in western Tennessee, was conducted from October 1987 through September 1989. Concentrations of selected herbicides also were defined. Reelfoot Lake, with a surface area of about 15,500 acres, is the largest natural lake in Tennessee and an important recreation and fisheries resource. Previous studies showed that the lake is hypereutrophic, a condition caused by high concentrations of nutrients in water and sediments discharged from the three principal tributaries (South Reelfoot Creek, North Reelfoot Creek, and Running Slough) to the lake. Pesticides, including herbicides, have been detected in the lake's bottom sediments. Storm runoff contributed about 87percent of the total water discharge of the three main tributaries to Reelfoot Lake. South Reelfoot Creek contributed about 4.7 tons per acre per year of suspended sediment, while North Reelfoot Creek contributed about 1.9 tons per acre per year. Running Slough contributed only about 0.31 ton per acre per year of suspended sediment. Most of the suspended sediment was transported by storm runoff between October and March. About 80 percent of the annual streamflow of the three tributaries occurs during these months. The North Reelfoot Creek basin contributed 8.2 pounds per acre per year of total nitrogen and 2.4 pounds per acre per year of total phosphorus. South Reelfoot Creek basin contributed about 6.5 and 1.3 pounds per acre per year of total nitrogen and phosphorus, respectively, while Running Slough basin contributions were 3.4 and 0.86 pounds per acre per year, respectively. The differences in nutrient yields appear to result from more row-crop agriculture and the relatively steeply sloping agricultural land in the North Reelfoot Creek basin. Ninety-one percent of the total nitrogen load and 95 percent of the total phosphorus load in the three streams was transported by storm runoff. Significant differences in the mean concentrations of nutrients in runoff were defined between the active agricultural months (April through September) and the inactive months (October through March). Storm-runoff samples were analyzed for II selected triazine herbicides. Alachlor and atrazine were the most commonly detected herbicides. Thirty-two percent of the samples contained detectable levels of alachlor and 93 percent of the samples contained detectable levels of atrazine. Ninety percent of the samples collected during the active agricultural months contained detectable leveki of alachlor and all 29 samples contained detectable levels of atrazine. Sixteen samples exceeded lifetime health-advisory levels for atrazine in drinking water (3 micrograms per liter); two samples collected from the April IS, 1988, storm at North Reelfoot Creek and South Reelfoot Creek contained 42 and 57 micrograms per liter of atrazine, respectively. Concentrations of the other nine triazine herbicides were generally less than the level of detection (0.1 microgram per liter).

Tennessee↗
USGS2024

Framework for mapping liquefaction hazard–Targeted design ground motions

Liquefaction-induced ground failure poses substantial challenges to geotechnical earthquake engineering design. Current approaches for designing against liquefaction hazards, as specified in most seismic provisions, focus on estimating a liquefaction factor of safety ( 𝐹⁢𝑆𝐿 ) and typically characterize earthquake loading using design parameters based on probabilistic or deterministic ground motion levels. Because 𝐹⁢𝑆𝐿 is estimated deterministically, this basis of design neglects considerable uncertainties for estimating liquefaction triggering and its consequences and results in a lack of liquefaction-specific design criteria, particularly as structural design has advanced toward risk-targeted performance objectives. This study presents a framework for developing liquefaction-targeted design criteria based on a minimum acceptable return period of liquefaction, informed by probabilistic liquefaction hazard analysis (PLHA). PLHA quantifies annualized rates of liquefaction by considering contributions from (1) the full ground-motion probability space, and (2) uncertainties in liquefaction triggering using probabilistic models. PLHA is used in this study to characterize the current, effective return periods of 𝐹⁢𝑆𝐿 ( 𝑇𝑅,𝐹⁢𝑆 ) obtained from conventional liquefaction hazard analysis (CLHA) using uniform-hazard ground motions. 𝑇𝑅,𝐹⁢𝑆 is evaluated in a parametric study of nearly 100 sites throughout the conterminous United States. The results indicate large geographic variations in acceptable liquefaction hazard levels, with implied 𝑇𝑅,𝐹⁢𝑆 ranging between approximately 1,000 to 3,000 years. To address these inconsistencies without the computational demands of full PLHA, a framework is proposed for developing a liquefaction-targeted design peak ground acceleration, 𝑃⁢𝐺⁢𝐴𝐿 , for use in liquefaction models that result in consistent liquefaction design levels across all geographic locations. The mapped 𝑃⁢𝐺⁢𝐴𝐿 is shown to be somewhat sensitive to site-specific properties, and adjustment factors are developed and presented. The proposed 𝑃⁢𝐺⁢𝐴𝐿 mapping procedure produces 𝐹⁢𝑆𝐿 estimates that are consistent with those obtained from full PLHA at a target 𝑇𝑅,𝐹⁢𝑆 , providing a promising roadmap to incorporating PLHA concepts into current liquefaction design methods.

Journal of Geotechnical and Geoenvironmental Engin↗
USGS2025

Disparate groundwater responses to wildfire

Post-wildfire investigations of groundwater response reveal a range of outcomes, varying from substantial increases to notable decreases in recharge and baseflow, with some studies indicating negligible or short-lived effects. This review assesses these varied responses within five critical categories: climate, vegetation, hydrogeology, fire characteristics, and the cryosphere, examining both short-term (within 2 years) and intermediate (2–10 years post-fire) effects. Despite considerable variability, some consistent patterns emerge. For instance, in hydroclimatic settings where water input and evaporative demand cycles are out of sync, post-wildfire groundwater responses tend to be positive (i.e., increased flux or storage), whereas under low fire severity conditions or in vegetation types that quickly recover, groundwater responses tend to be negative (i.e., decreased flux or storage). We synthesize relevant findings into a compendium of testable hypotheses aimed at explaining the spatiotemporal variability in observed post-wildfire groundwater responses. A recurring theme is the critical influence of the pre-wildfire groundwater regime on expected response and recovery. We identify opportunities for specific improvements in post-wildfire monitoring and modeling that would further advance capabilities to predict groundwater response. A key area for further research is understanding how wildfire effects on snow dynamics and other cryospheric processes translate to changes in groundwater.

WIREs Water↗
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↗
USGS2024

Subducting plate structure and megathrust morphology from deep seismic imaging linked to earthquake rupture segmentation at Cascadia

The origin of rupture segmentation along subduction zone megathrusts and linkages to the structural evolution of the subduction zone are poorly understood. Here, regional-scale seismic imaging of the Cascadia margin is used to characterize the megathrust spanning ~900 km from Vancouver Island to the California border, across the seismogenic zone to a few tens of kilometers from the coast. Discrete domains in lower plate geometry and sediment underthrusting are identified, not evident in prior regional plate models, which align with changes in lithology and structure of the upper plate and interpreted paleo-rupture patches. Strike-slip faults in the lower plate associated with oblique subduction mark boundaries between regions of distinct lower plate geometry. Their formation may be linked to changes in upper plate structure across long-lived upper plate faults. The Juan de Fuca plate is fragmenting within the seismogenic zone at Cascadia as the young plate bends beneath the heterogeneous upper plate resulting in structural domains that coincide with paleo-rupture segmentation.

Oregon, Washington↗
USGS2024

Deep syntectonic burial of the Anthracite belt, Eastern Pennsylvania

Fluid inclusion microthermometry and Raman spectroscopy of fluid inclusions in quartz veins from the Pennsylvanian rocks of the Anthracite belt, eastern Pennsylvania support a deep burial model of coalification in favor of focused orogenic hot fluid flow. High-temperature (250 to 255 °C) trapping of CH 4 ± CO 2 saturated aqueous fluids and CH 4 ± CO 2 inclusions indicate fluid trapping at depths of 11.5 to 13.4 km under a cover of Pennsylvanian to Permian(?) syntectonic load. In the folded rocks to the south of the Anthracite belt, CH 4 ± CO 2 fluid inclusions indicate a sediment load that was up to 16.3 km thick. Re -equilibrated aqueous fluid inclusions from veins in Silurian through Devonian rocks give the same range of trapping conditions but a wide range of fluid salinities suggesting that folding, fracturing, and meteoric recharge resulted in the intermixing of fluids from throughout the stratigraphic succession.

Pennsylvania↗
USGS2024

Shallow faulting and folding beneath south‐central Seattle, Washington State, from land‐based high‐resolution seismic‐reflection imaging

The geologic framework of the Seattle fault zone (SFZ) has been extensively studied, but the structure and fault strand locations in the central portion of the fault zone through the city of Seattle have remained controversial. Much of what is known about the SFZ has come from light detection and ranging (lidar)‐topographic surveys and paleoseismic investigations of fault scarps primarily west of Puget Sound, regional gravity and aeromagnetic modeling, and multiscale marine seismic imaging in waters both west and east of Seattle. We analyze ∼24 km of land‐based P ‐wave seismic‐reflection data that fill in a critical gap in our understanding of the SFZ beneath the urban areas of West Seattle, south‐central Seattle, and Mercer Island. These data image deformed strata in the upper 1 km, including upwarped Tertiary rock and younger sediments. Collectively, these data provide evidence for multiple Quaternary‐active thrust faults, back thrusts, and sub‐basins within the SFZ beneath the city of Seattle. The results indicate that multiple and potentially active back thrusts in the upper ∼500 m extend across the length of the SFZ and the entire urban corridor that may be analogous to those on Bainbridge Island west of Puget Sound.

Washington↗
USGS2026

Development of a hydrologic flow model of the Zuñi Salt Lake and surrounding areas, west-central New Mexico

The terminal Zuñi Salt Lake is in a maar in west-​central New Mexico and contains hypersaline water that has long been used by Native Americans for religious purposes and the collection of salt. The U.S. Geological Survey (USGS), in cooperation with the Zuni Tribe of the Zuni Reservation, New Mexico, and Bureau of Reclamation, developed the Zuñi Salt Lake hydrologic flow model to simulate the steady-state conditions of the lake and surrounding groundwater-flow system. The model is a tool that can be used to analyze the potential hydrologic responses to resource development in the region under different water-use scenarios, which can support negotiations that could help protect future use of the lake. The USGS modular hydrologic model, MODFLOW 6, was used to simulate the hydrologic system of the lake and surrounding areas. To define the layering of the model, a three-dimensional hydrogeologic framework was constructed, which included seven informal hydrogeologic units and the locations of faults that truncate and offset the units. Calibration of the model was performed using the iterative ensemble smoother included in a parameter estimation software suite (PEST++). The iterative ensemble smoother approach resulted in a posterior parameter ensemble of model realizations with parameter and simulated values that show similar success in matching observations. Groundwater flux to the Zuñi Salt Lake simulated by using the posterior parameter distribution ranged from 26 to 730 acre-feet per year with a mean, median, and standard deviation of 530, 570, and 170 acre-feet per year, respectively. The relative contribution of groundwater discharge to the lake was simulated using a particle-tracking model for MODFLOW (MODPATH) to be from, in order of decreasing contribution, the Mesaverde, Cenozoic, Mancos, and Dakota hydrogeologic units. Simulations of future conditions using the ensemble of models can provide an empirical estimate of forecast uncertainty without substantial additional analysis.

New Mexico↗
USGS2025

Cooling perspectives on the risk of pathogenic viruses from thawing permafrost

Climate change is inducing wide-scale permafrost thaw in the Arctic and subarctic, triggering concerns that long-dormant pathogens could reemerge from the thawing ground and initiate epidemics or pandemics. Viruses, as opposed to bacterial pathogens, garner particular interest because outbreaks cannot be controlled with antibiotics, though the effects can be mitigated by vaccines and newer antiviral drugs. To evaluate the potential hazards posed by viral pathogens emerging from thawing permafrost, we review information from a diverse range of disciplines. This includes efforts to recover infectious virus from human remains, studies on disease occurrence in polar animal populations, investigations into viral persistence and infectivity in permafrost, and assessments of human exposure to the enormous viral diversity present in the environment. Based on currently available knowledge, we conclude that the risk posed by viruses from thawing permafrost is no greater than viruses in other environments such as temperate soils and aquatic systems.

mSystems↗
USGS2025

Salinas Valley integrated hydrologic and reservoir operations models, Monterey and San Luis Obispo Counties, California

The area surrounding the Salinas Valley groundwater basin in Monterey and San Luis Obispo Counties of California is a highly productive agricultural area, contributes significantly to the local economy, and provides a substantial portion of vegetables and other agricultural commodities to the Nation. This region of California provides about half of the Nation’s lettuce, celery, broccoli, and spinach each year. Thus, this agricultural area provides significant volumes of agricultural products not just for California but the entire United States. Changes in population and increased agricultural development, which includes a shift toward more water-intensive crops, and climate variability, have put increasing demand on both surface water and groundwater resources in the valley. This has resulted in water management challenges in the Salinas Valley that are predominantly related to distribution of water supply throughout the basin. Where and when the water is present in the surface and subsurface does not coincide with where and when the water is needed. To deal with the distribution issue, historically water has been used conjunctively in the valley. Conjunctive use is a water management strategy that coordinates surface water and groundwater use to maximize water availability. Groundwater is used throughout the Salinas Valley to meet water demands when surface water supplies are insufficient. Availability of surface water is constrained by climate. Precipitation and streamflow vary seasonally and year to year. Although there are two reservoirs in the Salinas Valley to capture and store water during wet periods, the only conveyance of reservoir water to coastal agricultural areas is the Salinas River. Increasing demand on groundwater and surface water resources throughout the Salinas Valley has resulted in undesirable effects of unsustainable water use, such as surface water depletion, groundwater level declines, storage depletion in the principal aquifers, and seawater intrusion. To address these escalating issues, local communities, water management agencies, and groundwater sustainability agencies are evaluating how to sustainably manage both their surface water and groundwater resources. To meet water demands and reduce undesirable effects of unsustainable water use, continued conjunctive management of surface water and groundwater would ideally incorporate strategies to deal with increases in demand and a variable climate. To evaluate the challenging water management issues in the Salinas Valley, the U.S. Geological Survey, Monterey County Water Resource Agency, and the Salinas Valley Basin Groundwater Sustainability Agency developed a comprehensive suite of models that represent the Salinas Valley Hydrogeologic system called the Salinas Valley System Model. The Salinas Valley Geologic Framework was developed to characterize the subsurface using various topographic and geologic data sources, including information on hydrogeologic units, their surfaces and extents, geologic structures, lithology, and elevations from borehole data and cross sections, as well as details on faults and existing models. The Salinas Valley Watershed Model simulates the entire Salinas River watershed. Monthly surface water inflows into the integrated hydrologic model domain were simulated using the Salinas Valley Watershed Model. The historical model uses historical climate data, water and land use data, and reservoir releases to simulate agricultural operations, including landscape water demands, diversions, and reclaimed wastewater. The operational model adds an embedded reservoir operations framework to the simulation of the historical model that allows specified operational rules to simulate reservoir releases and changes in reservoir storage. The operational model assumes current reservoir operations and constant land use, which differs from historical conditions. Thus, the operational model is a hypothetical baseline model that can be used by local water managers to evaluate and quantify potential benefits of water supply projects. Together, the geologic framework, watershed, historical, and operational models form a tool that can be used to simulate irrigated agriculture and associated reservoir operations of the integrated hydrologic system of the Salinas Valley.

EarthArXiv↗
USGS2026

Top Elevation of Glacial Till and Thickness of the Big Sioux Aquifer Delineated From Electrical Resistivity Tomography Surveys Near Sioux Falls, South Dakota, 2022 and 2025

The City of Sioux Falls, South Dakota, requested the U.S. Geological Survey perform electrical resistivity surveys on three parcels of land north of the city. Electrical resistivity data were collected along a total of 22 transects during March 14–18, 2022, and November 17–21, 2025. Results from electrical resistivity surveys were used to delineate the top of glacial till deposits for the purpose of characterizing the Big Sioux aquifer near the city. Delineating geologic contacts provides important information on groundwater storage, flow dynamics, well design and placement, contaminant transport, groundwater–surface-water interactions, and regional water modeling. The top elevation of glacial till and the thickness of the Big Sioux aquifer varied among the three survey areas. The interpreted top elevation of glacial till in the North survey area decreases from east to west toward a slough, with elevations ranging from 1,403 to 1,418 feet (ft). The estimated thickness of the Big Sioux aquifer in the North survey area increased from east to west, with thicknesses ranging from 23 to 38 ft. The top elevation of glacial till in the Well 72 survey area generally decreases from northwest to southeast. Top elevations of the glacial till in the Well 72 survey area ranged from 1,400 to 1,409 ft along the southern end of transect W72_2. The estimated thickness of the Big Sioux aquifer in the Well 72 survey area was greatest along a southeast to northwest trending channel, with thicknesses ranging from 28 to 40 ft. The top elevation of glacial till in the Nose survey area generally decreases west toward the Big Sioux River. Top elevations of the glacial till in the Nose survey area ranged from 1,362 to 1,395 ft. The estimated thickness of the Big Sioux aquifer in the Nose survey area ranged from 33 to 70 ft.

South Dakota↗
USGS2023

Explosive trachyte eruptions from the Al Efairia volcanic center in northern Harrat Rahat, Kingdom of Saudi Arabia

Harrat Rahat is an alkali basalt, continental, intraplate volcanic field located within the central-western part of the Kingdom of Saudi Arabia. The northern quarter of Harrat Rahat contains evolved volcanic products that achieve trachyte compositions (>60 weight percent SiO 2 ). Within the Al Efairia volcanic center, pyroclastic-flow and -surge deposits that reflect explosive trachyte volcanism (and minor exposed lava domes that reflect effusive trachyte eruptions) sit at the surface as the youngest expression of volcanic activity within this part of the Harrat Rahat volcanic field. Five trachyte deposits emplaced explosively have been identified within the Al Efairia volcanic center based on geologic mapping, petrography, geochemistry, and paleomagnetism. These units are the trachytes of Um Rgaibah, Gura 5, Gura 4, Al Efairia, and Al Qayf, in descending stratigraphic order. Here, we present 14 40 Ar/ 39 Ar analyses from four of these units, which yield eruption ages of 4.2±5.2 thousand years (ka) for the trachyte of Um Rgaibah, 79.7±1.6 ka for the trachyte of Gura 5, 84.3±1.6 ka for the trachyte of Gura 4, and 88.0±1.8 ka for the trachyte of Al Efairia. The eruption age of the trachyte of Al Qayf has been constrained to between 410.3±3.4 and 418.8±1.9 ka using paleomagnetic correlations and 40 Ar/ 39 Ar ages from overlying and underlying intermediate composition lava flows. Most of these trachytes have distinct geochemical compositions, petrographic characteristics, and directions of remanent magnetization. The exceptions are for the trachytes of Gura 4 and Gura 5, which overlap in their geochemical, petrographic, paleomagnetic, and geochronologic affinities. Based on these similarities, we interpret the trachytes of Gura 4 and Gura 5 to have erupted during a closely spaced (a few decades) time interval from the same magma batch but from craters that are >2 kilometers (km) apart. The eruption of the trachyte of Al Efairia at 88.0±1.8 ka is the result of a different magma batch that erupted a few thousand years prior to the trachytes of Gura 4 and Gura 5. The Al Efairia volcanic center is remarkably different from the Matan volcanic center located ~10 km to the north, which has also erupted young (<150 ka) trachytes. The Matan volcanic center has been shown to produce trachyte compositions only after eruption of basalt followed by intermediate lava flows, whereas only trachyte compositions have erupted within the Al Efairia volcanic center over this same time interval.

Al Efairia volcanic center, Harrat Rahat↗
USGS2025

3D Converted wave reverse time migration imaging

We describe a newly developed method for recovering high-resolution images of seismic discontinuities, such as subducting slabs, in 3D. Our method makes use of converted S → P or P → S waves observed by dense arrays of seismometers to infer the locations and relative strengths of seismic discontinuities at depth in a target region. Observed direct and converted waves are backpropagated to their times of origin. The time-reversed wavefield is then separated into its constituent P and S components via the Helmholtz decomposition, and those separated wavefields are used to compute imaging functions that characterize the locations and relative strengths of seismic discontinuities. Imaging functions may be designed to use either S → P or P → S waves, so that users can target those arrivals expected to be most dominant in a given dataset. We have previously demonstrated the efficacy of our method in two dimensions, and we now present a 3D implementation of our technique which addresses the significant computational challenges posed by the size of volumetric wavefield data in three dimensions. Through a series of synthetic examples, we demonstrate that our method is capable of recovering the fine scale structure of a subducting slab given realistic station coverage and earthquake sources. We investigate optimal seismic station geometries for our technique and explore image interpretability in regions with poor data coverage. We find that linear station geometries yield more optimal, interpretable imaging functions than collections of small arrays can. We also show that our method can successfully recover bothS → P or P → S images when realistic shear earthquake sources are used, and we explore the additional computational challenges presented by the high frequency content of S waves. Our results demonstrate the potential for our technique to recover high-resolution information about subducting slabs in real-world regions, given that relatively sparse seismic arrays with only approximately 100 stations are capable of recovering interpretable imaging functions from just a few realistic earthquake sources for multiple discontinuities at significant depth in an area of approximately 400~sq~km.

Seismica↗
USGS2025

Using subducting plate motion to constrain Cascadia slab geometry and interface strength

Subduction zones are home to multiple geohazards driven by the evolution of the regional tectonics, including earthquakes, volcanic eruptions and landslides. Past evolution builds the present-day structure of the margin, while the present-day configuration of the system determines the state-of-stress in which individual hazardous events manifest. Regional simulations of subduction zones provide a tool to synthesize the tectonic history of a region and investigate how geologic features lead to variations in the state of stress across the subduction system. However, it is challenging to design regional models that provide a force-balance that is consistent with the large-scale motion of surrounding tectonic plates while also not over-constraining the solution. Here, we present new models for the Cascadia subduction zone that meet these criteria and demonstrate how the motion of the subducting Juan de Fuca plate can be used to determine the along-strike variations in the viscous (long-term) coupling across the plate boundary. All successful models require lower viscous coupling in the northern section of the trench compared to the central and southern sections. However, due to uncertainties in the geometry of the Cascadia slab, we find that there is a trade-off between along-strike variation in viscous coupling and slab shape. Better constraints on the slab shape, and/or use of other observations are needed to resolve this trade-off. The approach presented here provides a framework for further exploring how geologic features in the overriding plate and the properties of the plate boundary region affect the state-of-stress across this and other subduction zones.

Cascadia subduction zone↗
USGS2025

Groundwater budget for the surficial aquifer surrounding Lake Nokomis, Minneapolis, Minnesota

During prolonged periods of above-average precipitation, rising groundwater levels have the potential to cause damage to and interfere with underground infrastructure and building foundations. To understand the relations between precipitation and groundwater in the vicinity of Lake Nokomis, the U.S. Geological Survey, in collaboration with the University of Minnesota, quantified five components of the groundwater budget: groundwater recharge, change in surficial aquifer storage, surficial aquifer groundwater discharge to Lake Nokomis, groundwater evapotranspiration, and groundwater discharge to underlying bedrock aquifers. Field data, geologic records, and empirical calculation methods were used to quantify groundwater budget components for April 2023 through April 2024. Lake water budget data indicate that Lake Nokomis is a flowthrough system during periods with no outflow through the weir, with groundwater inputs equal to outputs. Roughly 40 percent of precipitation that fell in the study area was added to the surficial aquifer as recharge. Uncertainty in the vertical hydraulic conductivity resulted in wide-ranging estimates (spanning three orders of magnitude) of water discharging from the surficial aquifer to the underlying bedrock aquifer. Drought conditions persisted for the duration of this study and were not representative of the conditions that motivated this study. This study is a start towards understanding relations between precipitation, Lake Nokomis levels, and groundwater levels that could affect local underground infrastructure.

Minnesota↗
USGS2024

Geologic map of the northwest flank of Mauna Loa volcano, Island of Hawai‘i, Hawaii

Mauna Loa, the largest active volcano on Earth, has erupted 34 times since written descriptions became available in A.D. 1832. The most recent eruption of Mauna Loa occurred on November 27, 2022, after a 38 year hiatus; it lasted for 12 days. Some eruptions began with only brief seismic unrest, whereas others followed several months to a year of increased seismicity. Once underway, Mauna Loa’s eruptions can produce lava flows that may reach the sea in less than 24 hours, severing roads and utilities. For example, lava flows that erupted from the Southwest Rift Zone in 1950 advanced at an average rate of 9.3 kilometers per hour (5.8 miles per hour); all three lobes reached the ocean within ~24 hours. Near the eruptive vents, the flows likely traveled even faster. In terms of eruption frequency, pre-eruption warning, and rapid flow emplacement, Mauna Loa has great volcanic-hazard potential for the Island of Hawai‘i. Volcanic hazards on Mauna Loa can be anticipated, and risk substantially mitigated, by documenting its past activity to refine our knowledge of the hazards, and by alerting the public and local government officials of our findings and their implications for hazards assessments and risk. The map of the north and west flanks of Mauna Loa shows the distribution and relation of volcanic and surficial sedimentary deposits. It incorporates previously reported work published as generalized small-scale maps and a more detailed map. Within the mapped area, lava has flowed from three different source regions: the Northeast Rift Zone (22 percent), the summit (64 percent), and radial vents (14 percent). All three have different points of origin which, in turn, affect the flow characteristics and periodicity of activity. The map area includes the uppermost part of the NERZ and extends from the highest elevation––13,040 feet at the south end of the Kokoolau quadrangle, just below the summit caldera––to the sea northwest and west of the summit. Lava that erupts from the north and west flanks typically flows to the west, northwest, or north, depending on the vent location. Both morphologic lava flow types—‘a‘ā and pāhoehoe—are present. Pāhoehoe units tend to spread out or widen in low-slope regions, such as in the saddle regions between Mauna Loa and Mauna Kea or between Mauna Loa and Hualālai. In comparison, ʻaʻā flows generally produce narrower flow lobes that have higher relief. This map is the fifth in a series of five maps that will cover Mauna Loa volcano. NOTE: Map sheet 1 contains lines and type with overprint. This feature may be turned on or off in the Adobe Acrobat page display preferences.

Hawaii↗