Geology Reports⌕ Search

SEARCH · Geology Reports

Results for “Geostandards and Geoanalytical Research”

Search indexed USGS publications on groundwater, aquifers, geologic maps, mineral resources and earthquakes. Explore source records by subject and place.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

1,661 records · Page 10Linked to original sources
USGS2026

Quantitative mineral resource assessment of lithium pegmatite deposits in the southern Appalachian orogen

The first quantitative mineral resource assessment for undiscovered lithium pegmatite deposits in the southern Appalachian region of the United States was conducted. Permissive tracts for lithium pegmatite deposits were delineated by integrating lithological, tectonic, geochemical, geophysical and mineral occurrence data. Lithium pegmatite prospectivity of the tracts was ranked with simplified mappable criteria, including proximity to Paleozoic felsic intrusions and major lithotectonic structures, stream sediment geochemical anomalies, and pegmatite occurrence data. The geospatial data and permissive tracts were used to estimate the number of undiscovered lithium pegmatite deposits. These estimates were integrated into probabilistic simulations along with a new global lithium pegmatite grade and tonnage dataset to quantify potential contained undiscovered lithium resources. An economic filter was applied to convert the probabilistic estimates of contained lithium into recoverable material. The identified lithium pegmatite resources for the Carolina Lithium and Kings Mountain deposits, North Carolina, contain 1589 thousand tons (kt) of Li 2 O. The median contained undiscovered resource for the southern Appalachian orogen was estimated to be 2240 kt Li 2 O. At 90% confidence, the region contains at least 130 kt Li 2 O, and 10,700 kt at 10% confidence. After applying economic filters, the median recoverable contained resource was 1430 kt Li 2 O, corresponding to approximately 201 years of current lithium imports for consumption in the United States. North and South Carolina are likely to contain most of these resources. Coarse data resolution and intra-state variations in the geological data contribute to uncertainty of undiscovered lithium pegmatite resources. Continued efforts to harmonize disparate geospatial datasets with updated or new information can improve the accuracy and precision of estimated undiscovered lithium pegmatite resources in the study area and at broader scales.

Alabama, Georgia, Maryland, North Carolina, South ↗
USGS2025

The GorDAS Distributed Acoustic Sensing experiment above the Cascadia locked zone and subducted Gorda Slab

The southernmost portion of the Cascadia Subduction zone in Northern California produces high rates of moderate and large earthquakes owing to subduction of the Gorda slab and deformation associated with the Mendocino Triple Junction. Distributed Acoustic Sensing (DAS) is rapidly advancing as a method for detecting earthquakes and imaging crustal structure. We have begun a long-term DAS monitoring experiment on buried telecom fiber in Arcata, California, with the goal of increasing the available recordings of moderate to large earthquakes as well as imaging seismogenic structures. We have recorded over a year's worth of data, including most aftershocks of the 2022 M w 6.4 Ferndale earthquake, though not the mainshock itself. The dataset includes numerous magnitude 3.5 and larger earthquakes including the 2023/01/01 M w 5.4 Rio Dell earthquake. Here we present initial results comparing an earthquake detection algorithm, run in real-time on the processing unit of the interrogator system, with both the ShakeAlert earthquake early warning system as well as a post-processed earthquake catalog developed with deep-learning phase-picker algorithms. The rapid onboard processing of the detector demonstrates the potential utility of DAS-based edge computing for earthquake early warning. We also verify the quality of the strain waveforms both in terms of peak amplitudes and waveform similarity using about five months of nodal seismometer data. These instruments were deployed roughly every 300 m along the ~15km long cable and validate large variations in peak strain over short distances that are seen in the DAS data. All data from time windows surrounding both the local and teleseismic earthquakes are publicly available, which will improve our understanding of both the performance of DAS systems in moderate earthquakes and earthquake hazards associated with the Gorda subduction zone.

California↗
USGS2026

A ground-motion model derived using a generalized mean rupture distance for large slab interface earthquakes

Source–station distance is a central input to ground‐motion models (GMMs) for predicting seismic shaking. GMM development uses distance metrics including the Joyner–Boore distance, which is the shortest distance from an observation point to the surface projection of the earthquake rupture, and R rup the shortest distance to the rupture in three dimensions. Thompson and Baltay (2018) proposed the generalized mean rupture distance R p to address observed near‐fault ground‐motion saturation. R p accounts for the contribution to the shaking of all parts of the rupture and provides a simple method for incorporating spatially variable slip. They used R p to develop a GMM for shallow crustal earthquakes, assuming uniform slip. Here, we investigate the improvement offered by an R p ‐based GMM for large subduction interface earthquakes by recalibrating the path term for a published GMM ( Parker et al. , 2022 ) using R p derived from distributed slip models (DSMs). Inspection of within‐event and total residuals indicates that the recalibrated model fits the data at least as well as an alternative in which the path term was recalibrated using R rup and the same dataset. Incorporating slip information in its entirety or trimming the DSM geometry is preferable to assuming uniform moment release on a prescribed model fault that extends beyond the actual ruptured area. R p tuned to minimize model uncertainty is closer to the maximum distance to rupture for peak ground velocity (PGV) than acceleration, possibly reflecting the contribution to high‐frequency shaking of slip on local asperities. However, when the moment is concentrated far from stations, PGV is fit adequately with R p closer to R rup ⁠ . Our results suggest incorporating slip‐derived moment release through R p could improve GMMs for slab interface events, especially if the implementation of R p ‐based models is refined using a larger dataset of earthquakes with greater geographic diversity to account for regional ground‐motion variations.

Seismological Research Letters↗
USGS2025

Sources and risk factors for nitrate, pathogens, and fecal contamination of private wells in rural southwestern Wisconsin, USA

Household well water can be degraded by contaminants from the land's surface, but private well owners lack means to protect the source water from neighboring disturbances. Rural residents of southwestern Wisconsin, USA, rely on private well water, and the combination of land use and fractured carbonate bedrock makes groundwater vulnerable to contamination. To identify the extent, sources, and risk factors of private well contamination, randomly selected wells sampled during two-day periods in fall (n = 301) and spring (n = 529) were analyzed for nitrate and indicator bacteria, and a subset (n = 138) was sampled across four seasonal events for analysis of pathogens and microbial source tracking markers by quantitative polymerase chain reaction. Risk factors representing land use, hydrology, geology, and well construction were analyzed for associations with contamination in multivariable models. The importance of risk factors varied by contaminant, illustrating the multifaceted nature of rural groundwater quality. Nitrate contamination was associated with agricultural land use, and wells with casings that extended below a shale aquitard accessed less contaminated water than those drawing water from above it. Human fecal microbes were detected in 64 wells (46%), and rainfall was the key risk factor for contamination, indicating that wastewater from septic systems was available to contaminate wells when transport conditions were favorable. Manure microbes from cattle/ruminants and pigs were detected in 33 and 13 wells, respectively, and concentrations increased with the hectarage of cultivated land near wells. Pathogen genes for viruses, bacteria, and protozoa were detected in 66 wells (48%), including more detections of zoonotic than human-specific pathogens, and human Bacteroides , an indicator of wastewater, was an equivocal predictor of pathogen presence in private wells. Characterizing important elements of the setting, like geology, and identifying sources and risk factors for contaminants can inform landscape-level policies to protect groundwater quality.

Wisconsin↗
USGS2026

Ground-motion simulations for the 2024 Mw 4.8 Tewksbury, New Jersey, earthquake

Ground-motion simulations of notable earthquakes in the central and eastern United States are limited and typically assume one-dimensional (1D) Earth structure. In this study, we use a three-dimensional (3D) seismic velocity model to better constrain the depth and focal mechanism of the April 5th, 2024, moment magnitude 4.8 Tewksbury earthquake and investigate the spatial variability of earthquake ground motions and the effects of nearby sedimentary basins. We perform earthquake ground-motion simulations up to 0.5 Hz using the 3D spectral-element wave-propagation solver SPECFEM3D over a region 280-km wide by 260-km long by 77-km deep. Topography and subsurface geophysical structure are assigned using the U.S. Geological Survey National Crustal Model with a minimum shear-wave velocity of 200 m/s. We use earthquake time series from 13 broadband seismic stations in the region that have a uniform azimuthal distribution and epicentral distances ranging from 76 to 131 km to compare with synthetics and explore the effects of 1D versus 3D seismic structure on focal mechanism and depth solutions. Ground-motion intensity metrics are also presented relative to the NGA-East ground-motion models (GMMs) currently used in seismic hazard assessments for the region. We find that the 3D model, which reveals a wide spatial variability of period-dependent ground motions, yields better predictions of earthquake ground motions relative to the 1D model and the NGA-East ergodic ground-motion model, with 76 percent reduction of residual variance in observed ground motions averaged over 3-, 5-, 7-, and 10-second periods. Use of the 3D model to solve for a focal mechanism yields a shallower focal depth at 4 km and a shallower east-dipping focal plane relative to the U.S. Geological Survey regional moment tensor and Global Centroid Moment Tensor. Our study demonstrates that use of 3D seismic velocity models can improve estimates of earthquake focal mechanisms, ground motions, and seismic hazard.

New Jersey↗
USGS2026

As above, so below? A framework for integrating long-term water quantity trends reveals divergent patterns in groundwater and low streamflow across the United States

Climate, land-use, and disturbance drive long-term global trends in groundwater levels and streamflow. At large scales, these trends are typically considered separately, despite the well-established concept that groundwater and surface water comprise a single resource. Joint trend assessment at national scales is challenging because it requires pairing and aggregating data from spatially disparate streamflow and groundwater monitoring sites for which no established framework exists. Here, we evaluate alternative approaches for integrating groundwater and streamflow data to enable joint trend analysis—a critical step toward understanding how water-budget components respond concurrently and interactively to environmental drivers. Mann–Kendall trends were computed for individual groundwater (annual mean depth) and streamflow (annual low of 7 d averages) sites across the U.S over 21- (2000–2020), 31- (1990–2020), and 41-year (1980–2020) periods. Regional Kendall trends were calculated using five spatially contiguous and noncontiguous regional classifications for aggregation based on subsurface (e.g. aquifer, geology) and surface (e.g. watershed, landscape) characteristics. Site-level results revealed contrasting trends, with tendencies toward increasing low flows (wetting) and increasing groundwater depths (drying). Agreement between streamflow and groundwater trends increased with regional aggregation and longer timeframes, though persistent skew toward streamflow wetting and groundwater drying remained. Results varied by region and trend period, with notable consistencies: unified drying in the West/Southwest and wetting in the Upper Midwest. Directional mismatches in long-term trends were prominent in the High Plains and Mississippi Alluvial Plain, whereas near-term mismatches were most evident in the Northwest. Aggregation by hydrologic landscape regions (HLR) yielded the greatest agreement between groundwater and streamflow trends. These findings indicate that coupled responses may represent combined influences of climate, relief, and geology, as captured by HLR, more strongly than geography or geology alone. Integrated water availability assessments may benefit from a multi-characteristic classification framework to treat groundwater and surface water as a unified resource.

Environmental Research: Water↗
USGS2025

Time-varying rates of organic and inorganic mass accumulation in southeast Louisiana marshes: Relationships to sea-level anomalies and tropical storms

Louisiana's coastal wetlands are complex systems that require a continuous input of organic and inorganic material to keep pace with relative sea-level rise. Coastal restoration projects such as sediment diversions are being implemented to mitigate land loss and increase availability of inorganic sediment to coastal wetlands, and marshes specifically rely on organic material to build soil volume and maintain surface elevation. Interannual-to-decadal sea-level anomalies such as hurricanes can affect marsh accretion, mineral deposition, and plant productivity. In this light, complex ecogeomorphic feedback controls whether a marsh surface is sustainable or eroded/drowns. This study performs some of the first differential vertical accretion rates (VARs) and organic and inorganic mass accumulation rates (MARs) over time in SE Louisiana marshes determined from the 210 Pb Constant Rate of Supply model, coupled with standard 137 Cs VARs. These accumulation rates over the past ∼100 years were measured from a total of six brackish and salt marsh locations in Barataria Basin near the proposed Mid-Barataria Sediment Diversion. They were then related to interannual sea-surface elevations at Grand Isle, Louisiana, over the last ∼60 years and recorded hurricane activity in the delta. Results show VARs range from 0.63 cm/y to 1.69 cm/y and total MARs range from 0.11 to 0.43 g/cm 2 /y. Temporally, VARs and MARs (total, inorganic, and organic) are characterized by gradual increases in rates with decreasing age along with episodic peaks in VARs and MARs. The findings of this study indicate that no relationship occurs between sea-level anomalies and VARs or organic and inorganic MARs; however, a strong relationship appears to occur between major hurricanes to VAR and MAR contributions. Furthermore, high water content (81 ± 8%) and organic-rich soils in the sediment cores highlight the significance of belowground biomass and associated pore volume in maintaining marsh elevation in the study area.

Louisiana↗
USGS2024

Evidence for low effective stress within the crust of the subducted Gorda plate from the 2022 December Mw 6.4 Ferndale earthquake sequence

Stress levels on and adjacent to megathrust faults at seismogenic depths remain a key but difficult to constrain parameter for assessing seismic hazard in subduction zones. Although strong ground motions have been observed to be generated from distinct, high-stress regions on the downdip end of the megathrust rupture areas in many great earthquakes, we lack direct constraints on the stress level in the lower seismogenic portion of the Cascadia megathrust. On 2022 December 20, a M w 6.4 strike-slip earthquake occurred near Ferndale, California in southern Cascadia and likely ruptured the Gorda slab crust in the lower seismogenic portion, providing an opportunity to assess the stress level in this region. Here, we relocate the Ferndale mainshock and the first two weeks of aftershocks using a high-resolution 3-D velocity model and estimate rupture dimensions, directivity, and stress drop for several M w 4-5 aftershocks and recent earthquakes. The aftershocks define a strike-slip fault in the slab crust striking ENE, consistent with the mainshock focal mechanism. The orientation of this fault is about 45 ° off the ideally oriented fault plane given the stress state in the slab. The aftershock zone is extensive and broad in the forward direction of the mainshock rupture but still constrained within the volume of high Vp/Vs within the slab crust. Our stress drop estimates are generally lower for M w 4-5 earthquakes located within the slab crust compared to those a few km deeper in the slab mantle. Combined, our results support a relatively low effective stress level in the vicinity of the megathrust in the lower portion of the seismogenic zone in southern Cascadia, likely due to elevated fluid pressures. Consequently, the ground motion in the onshore region above this low-stress seismogenic portion in southern Cascadia may not be as intense as that observed during great earthquakes in other subduction zones.

California↗
USGS2026

Distributed faulting of the northern West Napa Fault Zone in Napa Valley, California

Mapped surface ruptures from the 24 August 2014 M w 6.0 South Napa earthquake in the Napa Valley, California, show a 2‐km‐wide zone of distributed faulting in the southern and central West Napa fault zone (WNFZ). In the northern WNFZ at Hendry Winery (HW), however, the mapped 2014 surface ruptures encompass an ∼100‐m‐wide zone, implying significant narrowing of the near‐surface fault zone to the north. We present a tomographic shear‐wave velocity ( ⁠⁠ V S ) model and guided‐wave data that indicate the northern WNFZ is at least 400‐m wide, with multiple near‐surface fault traces. Our V S model shows that the 2014 surface ruptures are underlain by discrete low‐velocity zones (LVZs), and coincident guided‐wave data show that the LVZs carry fault‐zone guided waves. If nearby (<500 m) mapped faults to the east of HW are part of the WNFZ, the entire WNFZ is more than 1 km wide in the northern Napa Valley. WNFZ guided waves travel up to 38% slower than S body waves, and low‐strain guided‐wave shaking is up to five times stronger than the associated body‐wave shaking. Our data suggest that guided waves, traveling along distributed faults, may result in an increased shaking hazard over a 1‐km‐wide area of the northern Napa Valley during future significant earthquakes. In places, the 2014 surface ruptures were difficult to find one year after the earthquake, and paleoseismic trenching showed only weak evidence for faulting, which may not have been identified in trenches if the locations of the 2014 surface ruptures had not been previously mapped ( Prentice et al. , 2015 ). Guided‐wave and V S tomography data, however, show strong evidence for faulting beneath the 2014 surface ruptures and at locations to the east. Although paleoseismic trenching is the gold standard for identifying near‐surface faulting, methods such as peak ground velocities of guided waves may better identify immature near‐surface fault traces.

California↗
USGS2024

Geologic input databases for the 2025 Puerto Rico – U.S. Virgin Islands National Seismic Hazard Model update: Crustal faults component

The last National Seismic Hazard Model (NSHM) for Puerto Rico and the U.S. Virgin Islands (PRVI) was published in 2003. In advance of the 2025 PRVI NSHM update, we created three geologic input databases to summarize new onshore and offshore fault source information in the northern Caribbean region between 62°–70° W and 16°–21° N. These databases, of fault sections, fault‐zone polygons, and geologic estimates of fault activity (fault‐slip rate and earthquake recurrence intervals) at specific sites, document updates to fault parameters used in prior seismic hazard models in PRVI. Fault sources were reviewed from published studies since 2003, which document substantial changes to the understanding of fault location, geometry, or activity. New fault section sources were added for features that meet the criteria of (1) length ≥7 km, (2) unequivocal evidence of recurrent tectonic Quaternary activity, and (3) documentation that is publicly available in a peer‐reviewed source. In addition, we revised several broad areal sources, such as the Mona and Anegada extensional zones. The 2003 model included three fault sections and two fault‐zone polygons (areal sources). These databases include 35 fault sections, 6 fault‐zone polygons, and 51 earthquake geology sites. To characterize fault activity rates, slip‐rate bins were assigned based on landscape expression and paleoseismic trench observations for faults without published slip‐rate sites. Additional fault sources were evaluated but not included in these databases due to a lack of published information about fault location, geometry, or recurrent Quaternary activity. The PRVI NSHM 2025 geologic input databases describe crustal faulting; the geometries and coupling of Puerto Rico subduction zone and Muertos Trough models are considered in a separate database. Updates to the fault sections, fault‐zone polygons, and earthquake geology databases can help inform the location and recurrence rate of damaging earthquakes in the PRVI NSHM implementation.

Puerto Rico, U.S. Virgin Islands↗
USGS2024

Integrating Sr isotopes, microchemistry, and genetics to reconstruct Salmonidae species and life history

Recent approaches to fisheries research emphasize the importance of the coproduction of knowledge in building resilient and culturally mindful fisheries management frameworks. Despite widespread recognition of the need for Indigenous knowledge and historical reference points as baseline data, archaeological data are rarely included in conservation biology research designs. Here we propose a novel multiproxy method to learn from former fisheries stewards by generating archaeological data on past salmonid population parameters. We used a newly developed, high throughput qPCR (HT-qPCR) chip, originally designed for environmental DNA (eDNA), for species identification of archaeological salmonid vertebrae. We combine this with the laser ablation split-stream (LASS) approach to identify ocean-migration versus freshwater residency. We test this multidisciplinary approach using both contemporary and archaeological salmonid samples and new radiocarbon dates from the Tronsdal Site on the Skagit River, Washington State, USA. This is a useful approach for extracting information about Salmonidae species and life history diversity from archaeological remains to reconstruct historic baselines for several population parameters in anadromous species with long periods of freshwater residency. The approach outlined in this paper may be particularly useful for research investigating past fisheries dynamics, offering hundreds to thousands of years of temporal depth for modern fisheries management, harvest policies, restoration ecology, and conservation biology.

Idaho, Oregon, Washington↗
USGS2026

Surface rupture and slip distribution of the 2025 Mw7.7 Mandalay earthquake and updated length scaling of supershear earthquakes

The 2025 M w 7.7 Mandalay, Burma (Myanmar), earthquake ruptured 475 km of the central Sagaing fault and is the longest continental strike-slip rupture on record. The observed rupture length is 1.6–4.7 times the value expected (100–300 km) from existing length-magnitude scaling relations for strike-slip earthquakes. The earthquake resulted from shallow dextral faulting and ruptured bilaterally with supershear speeds south of the epicenter, rupturing close to three major cities in Myanmar and exposing over six million people to violent or extreme shaking. We report on the surface rupture character, length, and slip distribution based on sub-pixel correlation of Sentinel-2 (10 m) and Planet Dove (3 m) optical images and visual analysis of SkySat and WorldView (0.3–0.5 m) optical images. The earthquake had moderate surface slip (average = 3.3 m, maximum = 5.6 m, 25–75% range = 3.0–4.0 m), narrow deformation zone width (1–10 pixels in sub-pixel correlation and up to 190 meters for the detailed surface rupture mapping), and simple fault geometry (no stepovers or large changes in strike, 87% of the rupture that was mapped in detail is single-stranded). We attribute the extreme length of the Mandalay earthquake to supershear rupture speed, simple fault geometry, narrow down-dip width, and moderate surface slip. Based on a compilation of 25 supershear strike-slip earthquakes ( M w 6.5–8.6; 1979–2025), we find that the rupture length of supershear earthquakes does not fit empirical scaling relationships for strike-slip earthquakes that predict length from magnitude. A length-magnitude scaling relationship based on supershear earthquakes has a best fit of log 10 (surface rupture length) = 0.89 M w – 4.44, indicating that supershear earthquakes tend to be longer than their subshear counterparts for any given magnitude and thus may expose a greater population to shaking.

Mandalay↗
USGS2026

Ten new insights in climate science 2025

Interdisciplinary understanding is vital for delivering sound climate policy advice. However, navigating the ever-growing and increasingly diverse scholarly literature on climate change is challenging for any individual researcher. This annual synthesis highlights and explains recent advances across a variety of fields of climate change research. This year, the 10 insights focus on: (1) the record-warmth of 2023/2024 and the elevated Earth energy imbalance; (2) acceleration of ocean warming and intensifying marine heatwaves; (3) northern land carbon sinks under strain; (4) reinforcing feedback between biodiversity loss and climate change; (5) accelerated depletion of groundwater; (6) global dengue incidence; (7) global income losses and labour productivity declines; (8) strategic scaling of CDR; (9) integrity challenges in carbon credit markets and emerging responses; and (10) effective policy mixes for emissions reductions. The insights have been written to be accessible to researchers from different fields, serving as entry-points to specific topics, as well as providing an overview of the evolving landscape of climate change research. In the final section, the insights are used to develop overarching policy-relevant messages. This paper provides the basis for a science-policy report that was shared with all Party delegations ahead of COP30 in Belém, Brazil.

Global Sustainability↗
USGS2025

2023 Earthquake Ground-Motion Workshop for the Central and Eastern United States, with a focus on the Gulf and Atlantic Coastal Plains—Agenda and abstracts

The U.S. Geological Survey held a virtual workshop December 7–8, 2023, to share research and ideas about earthquake ground motions in the Central and Eastern United States, with a focus on the Atlantic and Gulf Coastal Plains. The workshop was organized to learn about potential regionalization of ground-motion characteristics (source, path, and site), consider new explanatory variables for site response, and hear and discuss updates on ground-motion research on the Atlantic and Gulf Coastal Plains. The workshop was organized into a series of contributed presentations and three panel discussions held during 2 days. This report documents the agenda, contributed abstracts, and panel summaries.

Scientific Investigations Report↗
USGS2025

U.S. Geological Survey science strategy to address white-nose syndrome and bat health in 2025–2029

Since its discovery in 2006, the fungal disease known as white-nose syndrome (WNS) has killed millions of bats. Of the 47 bat species native to the conterminous United States, Alaska, Hawaii, and Canada, 12 have been affected by WNS, including 3 endangered species and 1 proposed endangered species. WNS has also been detected in 40 States and 9 Canadian Provinces. U.S. Geological Survey (USGS) scientists have been critical in identifying the causal fungus for WNS ( Pseudogymnoascus destructans [Pd]), characterizing the effects of WNS, and tracking the spread of Pd in many bat populations in North America. The mission of the USGS WNS and Bat Health Science Team is to deliver integrated science in order to build resiliency into free-ranging bat populations through more effective WNS management, build capacity for bat health science, and enhance bat health information sharing across USGS science centers and cooperative research units as well as with stakeholders. The USGS can play an important role in supporting regional and national capacity building by providing resources and guidance to local, State, and Tribal management entities and by providing tools to enhance disease management. The USGS Ecosystems Mission Area’s Biological Threats and Invasive Species Research Program is the lead Federal program for free-ranging wildlife disease research and surveillance. As of 2024, guided by the science priorities set by the WNS Steering Committee, USGS scientists are engaged in a nationwide response to WNS. This work is done in close coordination with our partners at the U.S. Fish and Wildlife Service, National Park Service, Bureau of Land Management, U.S. Forest Service of the U.S. Department of Agriculture, U.S. Department of Defense, as well as State and Tribal agencies. In addition to conducting WNS research, the USGS is mapping the spread of WNS and coordinating the North American Bat Monitoring Program (NABat) to understand how WNS and other stressors affect the status and trends of native bats across their range. The USGS is supporting the national WNS response through four science goals: (1) provide situational awareness on the health of bat populations; (2) conduct ecological studies of bats along the gradient of disease vulnerability; (3) contribute actionable science to enhance the resiliency of bat populations; and (4) implement an adaptive, holistic approach to bat health.

Circular↗
USGS2025

Identifying precursors and tracking pulses of magma ascent in multidisciplinary data during the 2018–2023 phreatomagmatic eruption at Semisopochnoi Island, Alaska

The 2018–2023 phreatomagmatic eruptions at Semisopochnoi Island, Alaska produced abundant long-period (LP) seismicity, harmonic and broadband tremor, and explosion signals over several well-monitored periods of eruption and quiescence. The corresponding dataset provides an excellent opportunity to investigate precursory and syn -eruptive geophysical signals of long-lived phreatomagmatic eruptions using multiparameter observations. We generated explosion and LP event catalogs through novel implementations of the REDPy ( Hotovec-Ellis, 2024 ) repeating event detector in mid-2021 following a network upgrade and the onset of a new phase of the eruption. The hundreds of detected explosions show a high degree of infrasound waveform similarity over more than a year, indicating a repeating source mechanism likely associated with explosive magma-water interaction. The seismic LP catalog shows that events began over a month prior to renewed explosive activity at the beginning of August 2021, and that lower frequency index (FI) LPs were generated in the week prior to the onset of explosions. We applied a recently developed machine learning tool (VOISS-Net, Tan et al., 2024 ) to catalog abundant broadband and harmonic seismic tremor recorded before and during the renewed explosive activity, along with LPs and explosions. The tremor catalogs complement the LP and explosion catalogs by filling out the seismic sequence with the dominant signal types. Together, these catalogs reveal a seismic sequence of renewed unrest that started with several weeks of LP events, followed by LPs with lower FI values and harmonic tremor in the days prior to explosive activity, and finally the onset of discrete explosions and broadband eruption tremor. We interpret this sequence as the ascent of a new pulse of magma that first interacted with the hydrothermal/groundwater system to produce LPs, followed by harmonic tremor, and that ultimately drove explosive magma-water interactions and periods of continuous ash emissions. The 2021 seismic sequence, in combination with long-term records of satellite SO 2 emissions, deformation from interferometric synthetic aperture radar (InSAR) analysis, ash sample analysis, infrasound, and volcano tectonic seismicity, allows us to interpret the entire 9-year period of unrest and eruption that began with an intrusion and earthquake swarm in 2014.

Alaska↗
USGS2023

Evaluation of models for estimating hydraulic conductivity in glacial aquifers from NMR logging

Nuclear magnetic resonance (NMR) logging is a promising method for estimating hydraulic conductivity ( K ). During the past ∼60 years, NMR logging has been used for petroleum applications, and different models have been developed for deriving estimates of permeability. These models involve calibration parameters whose values were determined through decades of research on sandstones and carbonates. We assessed the use of five models to derive estimates of K in glacial aquifers from NMR logging data acquired in two wells at each of two field sites in central Wisconsin, USA. Measurements of K , obtained with a direct push permeameter (DPP), K DPP , were used to obtain the calibration parameters in the Schlumberger-Doll Research, Seevers, Timur-Coates, Kozeny-Godefroy, and sum-of-echoes (SOE) models so as to predict K from the NMR data; and were also used to assess the ability of the models to predict K DPP . We obtained four well-scale calibration parameter values for each model using the NMR and DPP measurements in each well; and one study-scale parameter value for each model by using all data. The SOE model achieved an agreement with K DPP that matched or exceeded that of the other models. The Timur-Coates estimates of K were found to be substantially different from K DPP . Although the well-scale parameter values for the Schlumberger-Doll, Seevers, and SOE models were found to vary by less than a factor of 2, more research is needed to confirm their general applicability so that site-specific calibration is not required to obtain accurate estimates of K from NMR logging data.

Wisconsin↗
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↗