Search USGS⌕ Search

SEARCH · Search USGS

Results for “Journal of Geophysical Research Solid Earth”

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

Tectonic controls on volcanism and associated hydrothermal activity in a sediment-dominated mid-ocean ridge; Escanaba Trough

Mid-ocean ridges, the Earth's most extensive volcanic system, exhibit unique characteristics in sediment-dominant environments. Thick sediment cover insulates the crust and channels fluid along pathways that can lead to the formation of distinct crustal alteration patterns, exceptionally large mineral deposits, and specialized chemosynthetic ecosystems. This study presents an interdisciplinary investigation into the tectonics of the Escanaba Trough, a heavily sedimented axial valley at the southern Gorda Ridge in the Northeast Pacific Ocean. A primary challenge in such environments is overcoming the masking effect of thick sediments on basement structures that control magmatic and hydrothermal activity. We address this by employing three-dimensional (3D) magnetic modeling of high-resolution near-seafloor magnetic data collected by an autonomous underwater vehicle (AUV). The 2022 surveys with AUV Sentry provided data for 3D magnetic susceptibility models, refining our understanding of the geometry of sub-sediment laccoliths/saucer-shaped sills and hydrothermal alteration. In conjunction with a new 1:100,000 scale lithostratigraphic map, we outline the tectonic controls on the emplacement of Escanaba Trough's three main volcanic centers, characterize the geometry of its spreading segments, and provide volumetric data on the distribution of sub-sediment volcanism in the southern Gorda Ridge.

Escanaba Trough, Pacific Ocean↗
USGS2025

Near-surface geophysics: Environmental applications

The field of geophysics encompasses a broad and diverse compilation of methodologies that employs principles of physics to characterize properties of earth materials within the subsurface. While geophysical methods have a long history in resource exploration and studies of Earth’s interior, the subdiscipline of “near-surface geophysics” has evolved in recent decades for examination of the shallow, near-surface environment for a range of purposes ranging from archaeological or forensic investigations to assessment of geologic, hydrologic, biologic, and geochemical properties and processes. “Environmental geophysics” are near-surface geophysical studies and methods that focus on understanding natural systems (e.g., watershed hydrology, groundwater–surface water connections, biophysical processes) as well as research pertaining to anthropogenic impacts and land management, (e.g., contamination and remediation, saltwater intrusion, agricultural practices). This field can be further subdivided into subdisciplines focused on specific topics and applications, such as water resources and hydrology (hydrogeophysics) or biologic and microbial processes (biogeophysics). Studies in environmental geophysics span a range of scales, from pore-scale laboratory tests to watershed-scale or regional field experiments. Methods vary by the nature of physics employed, the specific measurement acquired, and how that data is ultimately processed and analyzed to produce interpretable results. There exists further diversity in the acquisition logistics, geometry, and timing of data collection. Geophysical data can be collected in boreholes (one-dimensional, 1-D, vertical profiles), along survey lines (two-dimensional, 2-D, cross-sections), or in dense sensor arrays or gridded profiles (three-dimensional, 3-D, models). Regarding the temporal aspect, studies can conduct one-time geophysical surveys to obtain detailed imaging of subsurface structure or use timelapse and continuous monitoring to investigate variations in subsurface properties over time. The cumulation of all possible permutations of these factors (method, acquisition geometry, survey design, and target application) results in an immense diversity among environmental geophysical studies. Nevertheless, this field remains unified in the pursuit of understanding natural and human-impacted near-surface environments through geophysical investigations. Here we highlight some key references within environmental geophysics. Resources on geophysical theory, acquisition logistics, processing and inversion workflows, and example case studies are categorized into the most common geophysical classes within Geophysical Methods. Lastly, example references for the dominant types of applications in environmental geophysical studies are catalogued in Environmental Applications.

Book chapter↗
USGS2024

Imaging of seismic discontinuities using an adjoint method

For imaging of seismic discontinuities at depth, reverse time migration (RTM) is a powerful method to apply to recordings of seismic events. It is especially powerful when an extensive receiver array, numerous seismic sources, or both, permit adequate reconstruction of incident and scattered wavefields at depth. Reconstructing either the incident or scattered wavefield at depth becomes less accurate when relatively few recordings of seismic events are available. Here we explore an inverse scattering approach to imaging discontinuities based on an adjoint method, employing sensitivity kernels (Frechet derivatives) that represent jumps in material properties across seismic-discontinuity surfaces. When combined with ray-based requirements on scattering geometry, it constitutes a powerful approach to determining the locations and amplitudes of the discontinuities, recovering only those properties that can be resolved by a spatially limited source and/or receiver distribution. This is illustrated by synthetic examples with local sources followed by a field example in a subduction zone setting

Washington↗
USGS2026

An entropic explanation for Gutenberg-Richter scaling

We develop a simple explanation for Gutenberg-Richter (G-R) size scaling of earthquakes on a single fault. We discretize the fault and consider all possible contiguous ruptures at that level of discretization. In this static model, we assume that slip scales with rupture length, and that the rupture rates at each point along the fault are consistent with an a priori long-term slip rate. These simple assumptions define an (under-determined) non-negative least-squares inverse problem. Each solution to this inverse problem is a set of earthquake rates that matches the slip-rate constraint. We use a Markov Chain Monte Carlo (MCMC) algorithm to uniformly sample the solution space assuming constant slip rates along the fault. At finer discretizations, deviations from G-R behavior decrease, which is consistent with an entropic pressure towards G-R solutions. When the fault is discretized into 10 or more segments, random solutions found by the MCMC algorithm have G-R size scaling, even though there are trivial solutions that, for example, have earthquakes of only one size. This is because there are simply far more solutions that have G-R scaling; as the problem size increases, the strong degeneracy of GR solutions results in other solutions becoming improbably rare. Also, the entropically favored G-R distribution has a b -value of approximately 1, which agrees with measured b -values in real earthquake catalogs.

JGR Solid Earth↗
USGS2025

Induced earthquakes are generally not tidally triggered in Oklahoma and Kansas

Human-induced earthquakes occur along critically stressed faults as injected wastewater simultaneously heightens fluid pressure and pushes faults to failure. We investigate the possibility that small stresses imposed by Earth tides could trigger earthquakes in the induced seismicity region of Oklahoma and Kansas from 2011 to 2018. We decluster a catalog consisting of ∼110,000 earthquakes using three methods (Reasenberg, nearest-neighbor distance, and phase-bin). We find no significant tidal earthquake triggering using Schuster's p -value test for the declustered catalogs as a whole. We search for localized triggering using discretized space-time cells and find ∼0–6% of cells have significant tidal triggering which is close to what is randomly expected (5%) and indicates there is an insignificant amount of tidal triggering for the full study region. One area that has significant p -values across multiple time windows, ∼2014–2016 is ∼15 km from a region of large wastewater injection volume. It is possible that localized tidal triggering occurs for this time and area because faults remain critically stressed and are particularly susceptible to slip under the small stress load from the semidiurnal tide. Possible explanations for the lack of tidal triggering in our broader study are that the pre-seismic stressing rate in the earthquake nucleation area is faster than the tidal stressing rate (∼3 kPa/day), faults are not close enough to critically stressed to be affected by tidal forcing, and that nucleation occurs over longer periods than the tides considered in this study (∼1, ∼14 days). Fluid injection could be the source of a higher pre-seismic stress rate.

Kansas, Oklahoma↗
USGS2026

A temporal look at the influence of topographic amplification on earthquake-triggered landslides in 3D seismic simulations

Earthquakes are a primary trigger for landslides, often leading to catastrophic consequences. While numerous studies have explored the spatial distribution of earthquake-triggered landslides, understanding the interaction between seismic waves and topography remains a critical challenge. Topographic irregularities can cause seismic wave amplification, altering ground shaking, and can trigger landslides that are challenging for predictive models to anticipate. This study investigates the spatial and temporal evolution of topographically amplified landslides, focusing on coseismic landslides triggered by the Mw 7.5 mainshock of the 2018 Papua New Guinea earthquake and post-seismic landslides associated with its four aftershocks, each exceeding Mw 6.0. We employ low-frequency, three-dimensional numerical ground shaking simulations and data-driven multivariate analyses to examine how landslides evolved from the coseismic to post-seismic periods. Our findings reveal a spatial shift in landslide distribution, in which the mainshock triggered slope failures predominantly on steep hillslopes, and the aftershocks triggered landslides on gentler slopes, often near geologic boundaries. We attribute this transition partly to the earthquake legacy effect of the mainshock, where the mainshock caused weakening of these hillslopes, making them more prone to failure when aftershocks occur. Additionally, the concentration of failures along geologic contacts in the post-seismic phase suggests that site amplification, stemming from contrasts in subsurface materials, exerts a key influence on landslide occurrence. Although not explicitly captured in our current numerical simulations, this mechanism warrants further investigation for more accurate hazard modeling.

JGR Solid Earth↗
USGS2024

Repeated coseismic uplift of coastal lagoons above the Patton Bay Splay Fault System, Montague Island, Alaska, USA

Coseismic slip on the Patton Bay splay fault system during the 1964 M w 9.2 Great Alaska Earthquake contributed to local tsunami generation and vertically uplifted shorelines as much as 11 m on Montague Island in Prince William Sound (PWS). Sudden uplift of 3.7–4.3 m caused coastal lagoons along the island's northwestern coast to gradually drain. The resulting change in depositional environment from marine lagoon to freshwater muskeg created a sharp, laterally continuous stratigraphic contact between silt and overlying peat. Here, we characterize the geomorphology, sedimentology, and diatom ecology across the 1964 earthquake contact and three similar prehistoric contacts within the stratigraphy of the Hidden Lagoons locality. We find that the contacts signal instances of abrupt coastal uplift that, within error, overlap the timing of independently constrained megathrust earthquakes in PWS—1964 Common Era, 760–870 yr BP, 2500–2700 yr BP, and 4120–4500 yr BP. Changes in fossil diatom assemblages across the inferred prehistoric earthquake contacts reflect ecological shifts consistent with repeated draining of a lagoon system caused by >3 m of coseismic uplift. Our observations provide evidence for four instances of combined megathrust-splay fault ruptures that have occurred in the past ∼4,200 years in PWS. The possibility that 1964-style combined megathrust-splay fault ruptures may have repeated in the past warrants their consideration in future seismic and tsunami hazards assessments.

Alaska↗
USGS2025

A high-resolution 3-D P-wave velocity structure of the south-central Cascadia subduction zone from wide-angle shore-crossing seismic refraction data

This study addresses a significant gap in understanding the features of the south-central Cascadia subduction zone, a region characterized by complex geologic, tectonic, and seismic transitions both offshore and onshore. Unlike other segments along this margin, this area lacks a 3-D velocity model to delineate its structural and geological features on a fine scale. To address this void, we developed a high-resolution 3-D P-wave velocity model using active source seismic data from ship-borne seismic shots recorded on temporary and permanent onshore seismic stations and ocean-bottom seismometers. Our model shows velocity variations across the region with distinct velocity-depth profiles for the Siletz, Franciscan, and Klamath terranes in the overlying plate. We identified seaward dipping high-velocity static backstops associated with the Siletz and Klamath terranes, situated near the shoreline and further inland, respectively. Regions of reduced crustal velocity are associated with crustal faults. Moreover, there is significant along-strike depth variation in the subducting slab, which is about 4 km deeper near the thick, dense Siletz terrane and becomes shallower near the predominantly less-dense Franciscan terrane. This highlights a sudden tectonic and geologic transition at the southern boundary of the Siletz terrane. Our velocity model also indicates slightly increased hydration, though still minimal, in both the oceanic crust and the upper mantle of the subducting plate compared to other parts of the margin.

British Columbia, California, Oregon, Washington↗
USGS2025

False positives in the identification of dynamic earthquake triggering

Dynamic earthquake triggering is commonly identified through the temporal correlation between increased seismicity rates and global earthquakes that are possible triggering events. However, correlation does not imply causation. False positives may occur when unrelated seismicity rate changes coincidently occur at around the time of candidate triggers. We investigate the expected false positive rate in Southern California with global M ≥ 6 earthquakes as candidate triggers. We compute the false positive rate by applying the statistical tests used by DeSalvio and Fan (2023), https://doi.org/10.1029/2023jb026487 to synthetic earthquake catalogs with no real dynamic triggering. We find a false positive rate of ∼3.5%–8.5% when realistic earthquake clustering is present, consistent with the 95% confidence typically used in seismology. However, when this false positive rate is applied to the tens of thousands of spatial-temporal windows in Southern California tested in DeSalvio and Fan (2023), https://doi.org/10.1029/2023jb026487 , thousands of false positives are expected. The expected false positive occurrence is large enough to explain the observed apparent triggering following 70% of large global earthquakes (DeSalvio & Fan, 2023, https://doi.org/10.1029/2023jb026487 ), without requiring any true dynamic triggering. Aside from the known triggering from the nearby El Mayor-Cucapah, Mexico, earthquake, the spatial and temporal characteristics of the reported triggering are indistinguishable from random false positives. This implies that best practice for dynamic triggering studies that depend on temporal correlation is to estimate the false positive rate and investigate whether the observed apparent triggering is distinguishable from the correlations that may occur by chance.

JGR Solid Earth↗
USGS2026

Fully contained laboratory earthquakes: The effect of asperity aspect ratio and free surfaces

Corner frequency ( f c ) and seismic moment ( M 0 ) are key parameters derived from seismic signals that are used to characterize earthquake stress drop, rupture area, and slip. These parameters are also affected by fault geometry and boundary conditions. However, the systematic study of these effects in laboratory settings has been challenging. This study presents laboratory earthquake experiments that examine how rupture dynamics are influenced by (a) the aspect ratio of rectangular PMMA velocity-weakening (VW) asperities surrounded by the Teflon velocity-strengthening (VS) patches, and (b) whether the sides of a VW asperity are confined with VS patches or are free surfaces. We found that increasing confinement by reducing free surfaces or increasing the VW asperity aspect ratio stabilizes fault slip, so that higher normal stress is required to transition from aseismic to seismic slip. Increased confinement and high aspect ratios also reduced M 0 and increased f c , both of which were determined from the radiated seismic waves. M 0 and f c were primarily controlled by the shorter dimension of the VW asperity. Analysis of high-frequency acoustic emission signals revealed that ruptures on high-aspect-ratio VW asperities propagated more unidirectionally, whereas ruptures on square VW asperities were more complex. Further, the high-aspect-ratio asperities were more likely to be eroded by surrounding VS regions while low-aspect-ratio asperities were more likely to rupture into the VS surroundings. These results demonstrate that both the confinement from surrounding stable areas and the geometry of the seismogenic patch can affect rupture nucleation, propagation, and seismic source characteristics.

JGR Solid Earth↗
USGS2025

To heal or not to heal?: 1. The effect of pore fluid pressure on the frictional healing behavior of lithologies in Oklahoma

The competition between fault healing (i.e., re-strengthening) and fault loading determines the timing and magnitude of fault failure within the seismic cycle. Repeating earthquakes can give observational estimates of fault healing rates, however, it is difficult to link laboratory studies of frictional healing and observed healing rates from repeating earthquakes in part because of uncertainty in lithology at depth. Due to well-constrained and relatively simple geology, earthquakes in Oklahoma can be linked to the granitic basement rock and to the Arbuckle Group, which is primarily composed of dolomite at earthquake depths. Here, we conduct friction experiments to measure healing rates of the two earthquake-bearing lithologies at confining pressures representative of earthquake depths and pore pressures ranging from 0% to 80% of the confining pressure. We measure frictional healing by executing slide-hold-slide tests with hold times ranging from 3 s to 3000 s. The friction experiments on the Troy Granite indicate that pore fluid pressure does not greatly affect healing rate. On the other hand, the dolomite of the Arbuckle Group exhibits decreased healing with increased pore fluid pressure, with weakening at the highest pore pressure. We hypothesize that this is due to an increase in dissolution of dolomite at high pore pressures/low effective normal stress. These healing rates are used in the companion paper to understand the moment-recurrence time behavior of repeating earthquakes in Prague, Oklahoma. This work has implications for possible enhanced dissolution and weakening behavior of the Arbuckle Group during wastewater injection activities.

Oklahoma↗
USGS2025

Three-dimensional seismic velocity models for the San Francisco Bay region, California from joint body-wave and surface-wave tomography validated by waveform simulation

A high-resolution seismic velocity model for both P and S waves is essential for the San Francisco Bay (SFB) region to accurately simulate earthquake ground motions and assess seismic hazards in this high-risk region. We present new three-dimensional P-wave and S-wave velocity (Vp and Vs) models of the SFB region developed using joint tomographic inversion with a rich set of earthquake and active-source body-wave data and ambient-noise surface-wave data. We designed a dynamic balance strategy to effectively balance the contributions of body-wave and surface-wave data during each iteration of the inversion. Owing to the complementary resolution of body-wave and surface-wave data, our new models are significantly improved compared to previous tomographic models. Major geologic features in the SFB region, including low-velocity basins, high-velocity Franciscan Complex, granite, and gabbro rocks, and across-fault velocity contrasts, are imaged in our Vp and Vs models. Compared to the geology-based U.S. Geological Survey velocity model, our new model has better defined structures for some important areas that are of great concern for seismic hazard assessment, such as the Sacramento-San Joaquin Delta, and the 2014 magnitude 6.0 South Napa earthquake area. We validated our model by simulating waveforms for a set of moderate-magnitude earthquakes. The good performance of our model in waveform simulation indicates its potential for improving ground motion simulation and seismic hazard assessment in the SFB region.

California↗
USGS2025

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

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

Geophysical Journal International↗
USGS2026

From collision to subduction: Thermal-kinematic inversions constrain plate boundary structure and dip-slip activity in southwestern New Zealand

Constraints on three-dimensional fault system geometry and kinematics are required to understand both long-term tectonics and contemporary seismic hazard. In many cases, these parameters must be inferred from a combination of different data or analyses. Here, we invert low-temperature thermochronological data to quantitatively parameterize both the crustal-scale geometry and late Cenozoic dip-slip activity of the plate boundary system in southwest New Zealand. We find that the modern transpressive plate boundary system has accommodated the bulk of Australian-Pacific Plate convergence since establishment in the middle Miocene-early Pliocene. Shortening has occurred via (a) dextral-reverse slip on the moderate angle, listric central Alpine Fault décollement, and (b) reverse-slip on either the sharply curved, anti-listric northern Puysegur subduction interface or low-angle offshore thrust fault system that splays from the high-angle southern Alpine Fault at depth. We resolve late Cenozoic shortening rates on these structures that are indistinguishable, and interpret these to reflect the transfer of convergent deformation outboard of the predominantly strike-slip southern Alpine Fault. Our work helps elucidate the transition from strain localized oblique continental collision to strain partitioned subduction. We also provide an estimate of the across-strike width of the central-southern Alpine Fault section boundary at seismogenic depths. This is the first quantitative constraint on the geometry of a structural feature inferred to have arrested several major prehistoric earthquakes. This study furthers understanding of fault structure and dip-slip kinematics along the greatest sources of seismic hazard in and offshore of southwest New Zealand.

South Island↗
USGS2026

The United States Magnetotelluric Array and the National Impedance Map

The United States Magnetotelluric Array (USMTArray) data set, collected in the years 2006–2024, consists of more than 1,700 long-period magnetotelluric stations covering the entirety of the contiguous United States on a quasi-regular 70 km grid. Funding across multiple federal agencies was critical to sustaining this effort to its completion. Important components of the project included active guidance and participation from the MT community, the open and timely availability of all data, and the application of consistent instrumentation and robust data processing. Together with parallel advancement in the development of publicly available three-dimensional (3D) inversion codes, the USMTArray has revitalized the US magnetotelluric community and increased the visibility of magnetotellurics within the Earth-science community. Taken as a whole, these data are visualized as the National Impedance Map, which, together with a 3D synthesis conductivity model of the nation, reveals the electrical architecture of the contiguous US. USMTArray data are used by researchers worldwide for fundamental and applied studies, including investigations of continental architecture and evolution, estimation of hazards to critical infrastructure due to geomagnetic storms, and assessment of the nation's undiscovered geothermal and mineral resources. We here review the history and development of the project, discuss the challenges and successes in its execution, present the National Impedance Map and synthesis conductivity model, and highlight the breadth of research stemming from this rich data set.

conterminous United States↗
USGS2025

To heal or not to heal?: 2. The moment-recurrence time behavior of repeating earthquakes in the 2011 Prague, Oklahoma aftershock sequence is consistent with laboratory healing rates

The timing and failure conditions of an earthquake are governed by the interplay between fault reloading and restrengthening. The moment-recurrence time behavior of repeating earthquakes can give observational estimates of fault healing rates; however, it is difficult to link these observed healing rates to laboratory studies of frictional healing in part because of uncertainty in lithology. Here, we study the 2011 Prague earthquake sequence, which includes repeating earthquakes in the Arbuckle group and the granitic basement, and compare them to laboratory experiments on samples of the Arbuckle and Troy granite (representative of the basement rock) (Okamoto et al., 2025, https://doi.org/10.1029/2024JB030573 ). We find three spatially distinct groups of repeating earthquakes with different moment-recurrence behavior: (a) constant moment-recurrence time in the Arbuckle group, (b) scattered moment-recurrence time at the intersection of the foreshock-mainshock fault in the granitic basement, and (c) moment-predictable behavior outside of the foreshock-mainshock fault intersection also in the granitic basement. Our observation of stagnant healing for repeating sequences in the Arbuckle group is consistent with laboratory observations of low healing rates for moderately high pore fluid pressures in Arbuckle samples. For the moment-predictable group, the source radius that is required in order to match healing rates is consistent with source radius estimations when taking into account reasonable attenuation of the P- pulse width. Overall, we observe diverse healing behaviors in the seismic families that are consistent with laboratory healing rates, providing seismic evidence that contact-scale frictional mechanisms are relevant to large-scale earthquake dynamics.

Oklahoma↗
USGS2026

Cajon Pass and the southern San Andreas Fault System: Earthquake cycle stress accumulation and present-day loading

With over a century since the last major rupture affecting the wider Los Angeles region, tectonic stress has steadily built along the southern San Andreas and San Jacinto fault systems, raising concerns of an imminent large earthquake. Cajon Pass, located at the junction of these faults, represents a critical site for potential through-going ruptures in Southern California. We constructed new 4D earthquake cycle simulations using a 1000-year paleoseismic rupture history of the San Andreas Fault System (SAFS) to assess spatial and temporal variations in stress. A semi-analytic Fourier transform model was used to compute stress from 3D dislocations in an elastic plate overlying a Maxwell viscoelastic half-space, assuming a complete coseismic reset of resolved shear stress on ruptured elements. Results show highest stress accumulation north of Cajon Pass (∼1.8 MPa/100 years) due to greater slip rates, and lower rates south of Cajon Pass (∼1.0–1.5 MPa/100 years). By 2025, Coulomb stress is estimated at 2.8 MPa on the Mojave South (MOS) segment, 1.8 MPa on the North San Bernardino (NSB1) segment and 3.6 MPa on the San Jacinto Bernardino (SJB) segment. Segments accumulate stresses with characteristic ranges of pre-event stress interpreted as failure thresholds: 1.2–2.7 MPa for MOS, 0.4–1.6 MPa for NSB1, and 1.2–2.9 MPa for SJB. When the stress disparity between segments SJB and MOS narrows, the faults appear to rupture jointly, suggesting that stress levels may control how Cajon Pass acts as an earthquake gate. These results may inform seismic hazard assessments by linking stress evolution to fault interactions.

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