SEARCH · USGS Science

Results for “Journal of Earthquake Engineering (JEE)”

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

SourceResearch connectionFrom yearThrough year

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

1,705 records · Page 26Linked to original sources

USGS2026

Repeated coral bleaching events have eroded reef accretion potential in the Florida Keys

Coral bleaching is now the primary driver of coral-reef degradation globally, threatening the persistence of coral-reef functions and the invaluable ecosystem services they support. The Florida Keys have been impacted by four severe, regional-scale coral-bleaching events, which occurred in 1997–1998, 2005, 2014–2015, and 2023. Our aim in this study was to quantify the effects of those thermal-stress events on the persistence of coral-reef structures. We leveraged long-term coral-reef monitoring data to quantify the impacts of coral bleaching on coral cover, carbonate production, and reef-accretion potential throughout the Florida Keys from 1996–2024. Our study indicates that coral bleaching drove significant declines of reef-building corals, which successively diminished regional-scale reef-accretion potential (from 1.65 ± 0.33 mm y -1 in 1996) and culminated with a regional-scale transition to net reef erosion (-0.12 ± 0.08 mm y -1 ) following the 2023 event. Coral bleaching was the primary driver of declining reef-building capacity, with the cumulative impacts of the four coral-bleaching events accounting for >75% of total declines since 1996. The 1997–1998 event alone drove ~50% of overall declines, largely due to losses of the shallow-water ecosystem engineer, Acropora palmata . More moderate impacts of thermal stress on regional reef-accretion potential in 2005 (6%) and 2014–2015 (9%) could lend some support to the hypothesis that the remaining corals were more thermally tolerant; however, the unprecedented thermal stress event in 2023 caused widespread coral-bleaching and bleaching-related mortality, and caused an additional 12% of the total decline in reef-accretion potential. Stony-coral tissue loss disease accounted for an additional 6% of the decline, with the remaining 19% unexplained. Coral bleaching has already dramatically reshaped the structure and function of coral reefs in the Florida Keys over the last three decades. With the high likelihood that thermal stress will accelerate moving forward, the persistence of Florida’s essential coral-reef functions and services is uncertain.

Florida↗

USGS2025

Colored shaded-relief bathymetry and acoustic backscatter of Lake Sammamish, Washington

Evidence of strong earthquakes (such as underwater landslides and associated deposits) may be recorded within the lacustrine sediments of Pacific Northwest lakes. The floor of Lake Sammamish, Wash., an approximately 11 kilometer (6.8 mile) long, 2 kilometer (1.2 mile) wide, and 35 meter (114.8 feet) deep lake located in a populated region just east of Seattle, was mapped by the U.S. Geological Survey in November of 2021 to search for evidence of past earthquakes. Mapping was conducted using a SWATHplus-M 234-kHz interferometric side-scan sonar system was pole-mounted on the U.S. Geological Survey research vessel Parke Snavely, and the system collected full-coverage bathymetric and acoustic backscatter data, which were processed to 2-meter spatial resolution. Two maps were created, a colored shaded-relief bathymetric map showing lake floor morphology (sheet 1), and an acoustic-backscatter map showing backscatter intensities (sheet 2). The results may then be utilized together to investigate past earthquake activity.

Washington↗

USGS2024

Colored shaded relief bathymetry and acoustic backscatter of Ozette Lake, Washington

Offshore of the Pacific Northwest of the United States is the Cascadia Subduction Zone, a 1,000-kilometer-long tectonic boundary defined by a large fault, called a megathrust, that extends from the Mendocino Junction off northern California to the Nootka Fracture Zone off Vancouver Island, Canada (U.S. Geological Survey, 2023). The Juan de Fuca and Gorda oceanic plates to the west of this boundary subduct under the North America continental plate to the east. Several other smaller faults that cut through the North America plate crust also affect the region. Although their effects upon Ozette Lake are uncertain, geological evidence for past earthquakes, such as underwater landslides, may be found in Pacific Northwest lakes. Underwater landslides caused by past earthquakes should be well preserved in these relatively undisturbed lake environments. The floor of Ozette Lake, Washington, located along the Pacific coast of the United States, west of the Puget Sound region and about 140 kilometers east of the megathrust was mapped by the U.S. Geological Survey in July of 2019 to search for evidence of past earthquakes. Mapping was completed using a SWATHplus-M 234-kHz interferometric side scan sonar system pole-mounted on the U.S. Geological Survey research vessel San Lorenzo. The system collected full-coverage bathymetric and acoustic backscatter data that were processed to 2-meter spatial resolution (Dartnell and others, 2024). This two-map series displays the results of this mapping. A colored shaded relief bathymetry map (sheet 1) and an acoustic backscatter map (sheet 2) show the lake floor morphology and backscatter intensities, respectively, that can be analyzed for evidence of past earthquakes.

Washington↗

USGS2023

A comprehensive fault system inversion approach: Methods and application to NSHM23

We present updated inversion‐based fault‐system solutions for the 2023 update to the National Seismic Hazard Model (NSHM23), standardizing earthquake rate model calculations on crustal faults across the western United States. We build upon the inversion methodology used in the Third Uniform California Earthquake Rupture Forecast (UCERF3) to solve for time‐independent rates of earthquakes in an interconnected fault system. The updated model explicitly maps out a wide range of fault recurrence and segmentation behavior (epistemic uncertainty), more completely exploring the solution space of viable models beyond those of UCERF3. We also improve the simulated annealing implementation, greatly increasing computational efficiency (and thus inversion convergence), and introduce an adaptive constraint weight calculation algorithm that helps to mediate between competing constraints. Hazard calculations show that ingredient changes (especially fault and deformation models) are the primary driver of hazard changes between NSHM23 and UCERF3. Updates to the inversion methodology are also consequential near faults in which the slip rate in UCERF3 was poorly fit or was satisfied primarily using large multifault ruptures that are now restricted by explicit b ‐value and segmentation constraints.

Bulletin of the Seismological Society of America↗

USGS2024

Separating signals in elevation data improves supervised machine learning predictions for hydrothermal favorability

A recent study identified topography (land surface elevation above sea level) as an important input dataset (feature) for predicting the location of hydrothermal systems in the Great Basin in Nevada. Yet, topography is generally a result of more than one geological process and may consequently contain multiple distinct signals. For example, the geologic evolution of the Great Basin has produced both crustal thickening (i.e., regional-scale trends in elevation) and thinning via Basin and Range extensional faulting (i.e., valley-scale topographic relief). We postulate that these geologic processes may affect the occurrence of hydrothermal systems differently. Therefore, we separate the regional trend from the valley-scale signal in the Great Basin, and then use them separately to evaluate the importance of each as predictors for hydrothermal favorability. Our prior work applying supervised machine learning (ML) using the data from the Nevada Machine Learning Project demonstrated that employing a training strategy that randomly selects negative training sites produces better performing models for predicting hydrothermal favorability than a training strategy that uses expert-selected negatives. The models created using both training strategies exhibited a west-east geographic trend in the predictions for the favorability of hydrothermal resources. These models generally predicted higher favorability in western Nevada and lower favorability in eastern Nevada. This west-east trend in predicted favorability correlates with elevation across the Great Basin, which trends higher from west to east. By separating the original elevation feature into distinct features for elevation trend (i.e., regional-scale topography) and detrended elevation (i.e., valley-scale or local relative topography), we find that models using the separated topographic signals consistently outperform competing models that use the original elevation feature. Although western Nevada still exhibits higher favorability than eastern Nevada, using separated signals for regional elevation and local structure reduces the west-east prediction trend in the region and emphasizes structures associated with hydrothermal upflow. This work emphasizes how carefully engineering features to represent geological conditions relevant to hydrothermal systems allows ML algorithms to detect important patterns for predicting hydrothermal resource favorability and leads to better model performance.

Conference Paper↗

USGS2025

Evidence for fluid pressurization of fault zones and persistent sensitivity to injection rate beneath the Raton Basin

Subsurface wastewater injection has increased the seismicity rate within the Raton Basin over more than two decades, with the basin-wide injection rate peaked between 2009-2015. To understand the evolution of injection-induced earthquakes, we systematically analyzed 2016-2024 broadband recordings with a machine-learning-based phase picker and constructed a catalog with 95,993 earthquakes (-1≤ M L ≤4.3). We then inverted for full centroid moment tensors (CMT) for 90 M L ≥ 2 events, with a special interest in constraining the non-double-couple components via probabilistic metrics. Both relocations and CMT solutions support basement-rooted normal faults, including graben and half-graben structures. Furthermore, we observe the non-double-couple components that imply elevated pore pressure in the fault zones. An earthquake cluster emerged in the north-central basin in 2023, preceded by ~1-yr of increased injection volume from wells within 15km. Despite a basin-wide decrease in the injection volume, we highlights the persistence of seismicity that remains to sensitive to injection rates within the Raton Basin.

Colorado, New Mexico↗

USGS2024

How does the onset of offset influence geologic slip rates?

Geologic slip rates are typically based on the displacement accrued by a geomorphic or stratigraphic feature and the age of the offset feature. Because slip rates are commonly calculated by dividing the displacement of a faulted marker by its age, they contain two open time intervals: the elapsed time between the age of an offset feature and the age of the earthquake that displaced the feature, and the time between the present‐day and the most recent earthquake. Here, we explore the influence of including unconstrained open intervals in geologic slip rate calculations. We test the degree to which these open intervals affect geologic slip rates and their uncertainties, and we find that their influence depends primarily on mean earthquake recurrence intervals (RIs). Slip rates on faults with longer RIs, such as the Wasatch fault, can be greatly influenced by an increase of up to 20% when accounting for open intervals. In contrast, slip rates on faults with shorter RIs, such as the San Andreas fault, are only slightly influenced by the assumption that slip rates calculated over open intervals approximate those calculated over closed intervals. Our analyses indicate that faults with moderate slip rates (∼0.2–5 mm/yr) are sensitive to both open interval effects themselves, as well as methods to quantify and account for these effects. We re‐evaluate how slip rates are calculated and defined in displacement–time space using published deformation records. We explore the utility of assigning a probability distribution to the initiation of offset of the oldest faulted feature and the timing of the most recent earthquake (MRE). We find that calculating geologic slip rates without using probability distributions that capture the timing of the MRE and the onset of offset of the oldest faulted feature, especially on slow‐to‐moderate slip rate faults, can lead to systematic underestimation of average geologic slip rates.

Seismological Research Letters↗

USGS2025

Local, regional, and distal recordings of seismic unrest at Tau Island volcano, American Samoa

A seismic swarm near Taʻū Island, a volcanic island in eastern American Samoa, occurred from July to October 2022. The earliest unrest was noted as felt shaking reports in late July, and instrumentation varied in the beginning of the sequence as the U.S. Geological Survey (USGS) Hawaiian Volcano Observatory responded by installing temporary and then permanent seismometers to monitor the activity. This network variability made it difficult to characterize the earliest seismicity and contextualize the entire sequence to discriminate between an underlying tectonic or volcanic source. Here, we present results analyzing hydroacoustic detections from an International Monitoring System hydrophone array near Wake Island, 4500 km northwest of Taʻū Island volcano. Using least-squares beamforming analysis, we create a catalog of T-wave detections from the direction of Taʻū Island to track the earthquakes, some of which were located by the USGS National EarthquakeInformation Center. Both the rate and hydroacoustic pressures, which we interpret as a proxy for earthquake size, gradually increased from late July to August, peaking on August 19 (rate) and August 24 (size), before decreasing to background in late September. Minutes-long bursts of tremor were also contemporaneously recorded as local network data became avail-able on August 20. Tremor activity continued throughout the rest of August, peaking on August 25, before ending in earlySeptember. These tremor bursts were band-limited to ~ 1–5 Hz and recorded as S waves at a regional station on the island of Upolu in Samoa, 250 km to the west of Taʻū Island. Our results do not constrain the tremor locations, but comparisons of earthquake and tremor reduced displacements recorded locally and regionally suggest a deeper tremor source. We interpret the increase in earthquake size and rate, together with the occurrence, characteristics, and relative depth of the tremor to be the result of magmatic activity beneath Taʻū Island volcano.

American Samoa, Taʻū Island↗

USGS2024

Integration of rupture directivity models for the US National Seismic Hazard Model

Several rupture directivity models (DMs) have been developed in recent years to describe the near-source spatial variations in ground motion amplitudes related to propagation of rupture along the fault. We recently organized an effort towards incorporating these directivity effects into the USGS National Seismic Hazard Model (NSHM), by first evaluating the community's work and potential methods to implement directivity adjustments into probabilistic seismic hazard analysis (PSHA). Guided by this evaluation and comparison among the considered DMs, we selected an approach that can be readily implemented into the USGS hazard software, that provides an azimuthally varying adjustment to the median ground motion and its aleatory variability. This method allows assessment of the impact on hazard levels and provides a platform to test the DM amplification predictions using a generalized coordinate system, necessary for consistent calculation of source-to-site distance terms for complex ruptures. We give examples of the directivity-related impact on hazard, progressing from a simple, hypothetical rupture, to more complex fault systems, composed of multiple rupture segments and sources. The directivity adjustments were constrained to strike-slip faulting, where DMs have good agreement. We find that rupture directivity adjustments using a simple median and aleatory adjustment approach can impact hazard both from a site perspective and on a regional scale, increasing shaking off the end of the fault trace up to 30--40\% and potentially reducing it for sites along strike. Statewide hazard maps of California show that the change in shaking along major faults can be a factor to consider for assessing long-period (>ls) near-source effects within the USGS NSHM going forward, reaching up to 10--20\%. Finally, we suggest consideration of minimum parameter ranges and baseline requirements as future DMs are developed to minimize single approach adaptations, to enable more consistent application within both ground motion and hazard studies.

Earthquake Spectra↗

USGS2025

Mechanics and statistics of postseismic shaking

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

California↗

USGS2024

SCEC/USGS Community Stress Drop Validation Study: How spectral fitting approaches influence measured source parameters

Spectral source parameters used to estimate an earthquake’s stress drop (Δσ) can vary significantly across measurement approaches. The Statewide California Earthquake Center/U.S. Geological Survey Community Stress-Drop Validation Study was initiated to compare source parameter estimates, focusing initially on a dataset from the 2019 Ridgecrest earthquake sequence. As part of that validation effort, here we focus on one potential source of uncertainty: whether spectral fitting approaches alone, applied to a common set of spectra from the 2019 Ridgecrest sequence result in different source parameter estimates. By using a common set of benchmark spectra analyzed across a consistent frequency band of 1–40 Hz, we eliminate many sources of variability. A subgroup of validation study participants volunteered to estimate the low-frequency displacement (Ω0) and corner frequency ( f c ) by fitting a smooth function to benchmark displacement spectra. Participants used linear- or log-sampled spectra, assumed a Brune or Boatwright spectral model, and applied different misfit criteria. We compare 17 approaches used to estimate Ω0, f c , and Δσ for 54 earthquake spectra. Our results reveal that 35% of events have Δσ estimates within a factor of two, whereas others exhibit variations exceeding an order of magnitude. The variability in Ω0 and f c can largely be attributed to whether a spectrum is consistent with the smooth function of an idealized simple crack model. The trade-off between Ω0 and f c may be more pronounced when using linearly sampled spectra, as higher frequency spectral bumps control the fits. As expected, methods that assumed a Boatwright model tended to have lower Ω0 and somewhat higher f c compared to those assuming a Brune model, although resulting Δσ estimates are similar. When compared to the overall validation study results, the fitting approach alone may account for between 5% and 90% (25% on average) of the total variability in spectral Δσ.

California↗

USGS2025

Three-dimensional seismic velocity model for the Cascadia Subduction Zone with shallow soils and topography, version 1.7

The U.S. Geological Survey’s seismic velocity model for the Cascadia Subduction Zone provides P- and S-wave velocity ( V P and V S , respectively) information from 40.2° to 50.0° N. latitude and −129.0° to −121.0° W. longitude, and is used to support a variety of research topics, including three-dimensional (3D) earthquake simulations and seismic hazard assessment in the Pacific Northwest. This report describes an update to the previous version (v) 1.6 of the 3D seismic velocity model for the Cascadia Subduction Zone. This new model (herein referred to as v1.7) contains more detailed near-surface structure for improved earthquake ground motion modeling. Updated features include the addition of a new shallow soil velocity model in the top few hundred meters and the option of adding user-specified topography. Although v1.6 of the Cascadia seismic velocity model has a minimum V S of 600 meters per second (m/s), the new model (v1.7) has a minimum V S of approximately 40 m/s. Overall, this update will allow for more accurate ground motion estimates from 3D simulations of scenario earthquakes in the Cascadia Subduction Zone region.

British Columbia, California, Oregon, Washington↗

USGS2024

Assessment of the interconnection between Tampa Bay and the Floridan aquifer system: Historical groundwater data compilation and analysis, 1976–2022

The U.S. Geological Survey used existing data collected after the last major navigational channel modification in the mid-1980s to investigate groundwater levels and chloride concentrations in wells in the Floridan aquifer system and other aquifers beneath and near Tampa Bay. Tampa Bay is located on the west-central coast of Florida and provides access for commercial shipping. In 2021, the U.S. Army Corps of Engineers began to investigate alternatives to improve the efficiency of the deep-draft navigation channels within Tampa Bay, reduce costs, improve safety, and analyze the implications of modifying navigational channels. The Floridan aquifer system underlies Tampa Bay and is the primary source of public water supply in the region. Modifications to the channels have the potential to expose the Floridan aquifer system to the saltwater in Tampa Bay, with the potential to increase salinity in wells in the region. Other factors affecting the groundwater levels and the location of the freshwater/saltwater interface were also examined, including changes in sea level, groundwater extraction, and variations in climate. Groundwater levels and well-construction reports were used to identify whether different aquifer units are well-connected. Twelve wells had available data before and after the last major channel modifications, which took place in the 1980s, with six datasets of chloride concentration available in areas along the northern and eastern coastline of Tampa Bay, which is nearest to historical dredging activities. Of these six, Kendall’s τ and p -values indicated increasing trends in chloride concentration for three datasets (TR 11-2, TR 10-2, 51), no trend in chloride concentration for two datasets (TR 9-3, 50), and a decreasing trend in chloride concentration for one dataset (TR 9-1). The upward trends in chloride concentration observed for TR 10-2 and 51 are likely the result of changes in local groundwater withdrawals. Well TR 11-2 had a gradual increasing trend in chloride concentration, fresh groundwater throughout the period of record, and a 3- to 4-foot increase in hydraulic head during the period of record, possibly caused by the construction and control of the Tampa Bypass Canal, resulting in changes to the regional potentiometric surface.

Florida↗

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↗

USGS2026

The digital archivist: Automating legacy macroseismic data processing using large language models

Macroseismic data are a key resource to investigate shaking and damage from preinstrumental and early instrumental eras. However, data are often stored as inconsistently formatted reports describing observed shaking and damage, making manually parsing and interpreting accounts labor‐intensive. We introduce a novel workflow using Google’s Gemini 2.5 Pro large language model (LLM) to automate the extraction and structuring of macroseismic observations from summary reports. We apply this workflow to the 22 March 1957 M 5.3 Daly City, California, earthquake as a case study. We used Gemini to extract addresses, originally assigned modified Mercalli intensity values, and descriptions from each report. To address coordinate precision limits, addresses were geocoded via Google’s Geocoding application programming interface. This workflow yielded over 2300 geocoded intensity reports for the Daly City earthquake. We use the geocoded accounts, with the original report intensity assignments, to develop a shaking intensity map that in some respects rivals modern Did You Feel It? Maps. We also extract and present data for the 9 February 1971 M L 6.7 Sylmar, California, earthquake. Our results demonstrate the potential of LLMs for reliably extracting and analyzing large, unstructured macroseismic datasets. LLMs offer a scalable solution for rapidly digitizing macroseismic archives, enabling their broader use to constrain ground‐motion models in modern seismic hazard analysis and to improve our understanding of site effects in urban areas. The concepts explored here may also be applied to the handling of other legacy seismological and earth science data.

Seismological Research Letters↗

USGS2026

Only very strong shaking can break a tree

A handful of credible accounts, both recent and historical, describe trees being broken during earthquakes. Apart from landslides and tsunamis, such accounts are uncommon, suggesting that trees break only when ground acceleration is extreme or when trees are especially vulnerable to breakage. Here, I consider the question: what ground acceleration is required to break a living tree? I present simplified first‐principles calculations, considering the properties of common trees, focusing on tall, thin trees that can be approximated as uniform vertical cantilevers. The results suggest that, possibly excepting tall trees with heavy canopies, standing trees can only be broken by shaking approaching or exceeding 1 g , as anticipated by the rarity of documented accounts of snapped trees. I briefly consider several well‐documented instances of tree damage during strong earthquakes, drawing inferences about local ground motions.

The Seismic Record↗

USGS2025

Paleoseismology and paleogeodesy using coral microatolls

Establishing the rupture extent and slip distribution of individual paleo-earthquakes is vital for assessing fault behavior including the persistence of rupture segmentation, recurrence patterns, and similarity of successive events, key issues in both fault mechanics and hazard assessment. Techniques with high temporal and geodetic precision as well as a wide distribution of study sites are necessary to investigate past earthquakes in such detail. Coral microatoll growth is one of the best types of geologic record for paleoseismology and paleogeodesy given these needs, as it provides long, continuous, widely distributed records of centimeter-scale vertical tectonic motion with potentially annual-level temporal precision. This chapter describes the process of interpreting microatoll growth records to obtain time series of relative sea level, tectonic vertical deformation fields, and finally slip and coupling parameters on an underlying fault interface.

Book chapter↗

USGS2024

Seismic velocity changes from repetitive seismicity at Mauna Loa prior to and during its 2022 eruption

Mauna Loa’s short-lived eruption from late November to early December 2022 marked the culmination of nearly a decade of elevated seismic activity and geodetic inflation. The volcano has been monitored by a network of permanent, short period and broadband seismometers. I used the continuous waveform data from that network starting in 2012 to generate a catalog of seismicity that enhances the US Geological Survey Hawaiian Volcano Observatory’s public seismic catalog with four times the number of earthquakes, which were then grouped by waveform similarity. Analysis of subtle delays in the timing of arrivals of scattered waves between pairs of earthquakes in this catalog yields a history of small changes in the shallow seismic velocity structure of the volcano. Seismic velocities have been shown at other volcanoes to change during unrest and eruption. My results show a decrease in seismic velocity centered on the summit beginning in September 2022, corresponding to the onset of a vigorous precursory swarm of seismic activity and shallow inflation. During the eruption itself, I observe large changes due likely to dike opening along the northeast rift zone and deflation of the summit reservoir. However, seismic velocity changes associated with non-volcanic sources such as ground shaking from large earthquakes and meteorological influences at seasonal and diurnal time scales are also observed, and these dominate the velocity changes prior to the eruption. Proper accounting of these effects will be a requirement for use in real-time monitoring, and this work serves as a starting point in that endeavor for Mauna Loa.

Hawaii↗