Geology Reports⌕ Search

Geology topics

Jefferson C. Chang

Publications and source records attributed to Jefferson C. Chang.

13 recordsLinked to original sources

Geophysical insights into Kīlauea’s eruptive reawakening after its 2018 eruption and caldera collapse

A little more than two years after the destructive 2018 lower East Rift Zone eruption and summit collapse ended, renewed volcanic activity began at Kīlauea volcano, Hawaiʻi. During 2020–2025, activity at the volcano included six eruptions in the summit caldera, one in the Southwest Rift Zone, and one in the middle East Rift Zone, as well as multiple non-eruptive intrusions in all three locations. We categorize the volcano’s activity post-2018 into four phases based on eruption locations and style, as well as patterns of magma storage and migration inferred from geophysical data. In some cases, these observations illuminated parts of the volcano’s magmatic system for the first time since the advent of modern geophysical monitoring—particularly deformation datasets that include tilt, continuous GNSS, and InSAR. Analyses of these geodetic data and complementary seismic, geological, and geochemical observations provide an opportunity to explore summit and rift zone connectivity, rift-zone magma storage conditions, the variable geophysical manifestations of eruption onsets, magma-fault interactions, and other long-standing questions associated with eruptions and intrusions at Kīlauea.

Hawaii↗

Insights from more than a year of ongoing episodic lava fountaining at Kilauea

Episodic lava fountaining is rare globally, with only three occurrences from Kīlauea (Hawaii) since 1823. Beginning in December 2024, an eruption within the summit crater, Halemaʻumaʻu, has produced 52 fountaining episodes to date, each lasting typically for hours and separated by days to weeks. Although multiple mechanisms have been proposed, fundamental questions remain about how magma reservoir pressure, conduit permeability, and volatile exsolution control fountaining onset. We present a multi-parameter record of episodic fountaining at Kīlauea using modern geochemical, geophysical, and geologic monitoring. These observations revealed systematic correlations between fountaining and summit inflation, enabling unprecedented quantitative forecasts. This eruption will provide new constraints on the timing, evolution, and processes driving episodic fountaining globally, and highlights Kīlauea’s role as a natural laboratory for studying volcanic processes.

Hawaii↗

Commentary on paper by Matoza et al. (2021): Catalog revision to a common depth datum

We present a revision and update to the high-precision relocated seismicity catalog presented by Matoza et al. (2021, https://doi.org/10.1029/2020ea001253 ) for the Island of Hawai'i from 1986 to 2018. The starting catalog of hypocenters (input data), on which the study by Matoza et al. (2021, https://doi.org/10.1029/2020ea001253 ) was based, contained an inconsistent depth datum for events before and after 00:00 UT, 29 December 2017. Here we present a recomputed version of the catalog using a consistent reference depth. We corrected the starting catalog to a common depth datum (all events now use the model depth reference datum) and re-ran the entire workflow as described in the paper by Matoza et al. (2021, https://doi.org/10.1029/2020ea001253 ). This included pairing, cross-correlating, and relocating all seismic events again based on the updated starting catalog. We consider 347,446 events representing 32 years of seismicity on and around the island from 1986 to 2018. We now successfully relocate 299,966 (86%) events using ∼2.53 billion differential times ( P and S ) from ∼194 million similar-event pairs, derived from cross-correlations between ∼887 million event pairs total, a significant increase from our original analysis. The resolution of fine-scale seismicity features is improved and the median depth of shallow events (<5 km) under Kaluapele (Kīlauea summit caldera) in 2018 is shifted 926 m deeper as a result of the change. The interpretations and other major conclusions in the paper by Matoza et al. (2021, https://doi.org/10.1029/2020ea001253 ) are unchanged.

Earth and Space Science↗

A decade of precursory seismicity at Mauna Loa: 2013–2022

The 2022 eruption of Mauna Loa was preceded by nearly 10 years of seismic unrest. Using the earthquake catalog of the Hawaiian Volcano Observatory, we describe and explore many aspects of unrest between January 1, 2013, and the start of the eruption on November 28, 2022 (UTC). Earthquakes associated with unrest occurred above an inflating dike at shallow depths under and to the south of Moku‘āweoweo and in the cluster of seismicity to the west of Moku‘āweoweo associated with a buried rift zone. Earthquakes also occurred in a broad region at the base of the crust, likely associated with a zone of magma accumulation. Unrest is broken down into four phases of increasing intensity with time and into individual geographic regions. While activity occurring at shallow depths under the caldera often alternated with activity on the west flank in the first two phases of unrest, the third phase of unrest featured increased rates of earthquake activity in both areas acting concurrently. We also document a slow shift in seismicity across multiple years from the south of Moku‘āweoweo to an area under the caldera that likely represents a change in the location of the intrusion through time. Broadly speaking, unrest observed in advance of the 2022 eruption occurred in the same areas as the unrest preceding eruptions in 1975 and 1984. Furthermore, precursory unrest prior to eruptions in 1975, 1984 and 2022 lasted years and culminated in an intense phase of seismic activity lasting on the order of hours that represented the final dike propagation to the surface. Concurrent earthquake activity under Moku‘āweoweo and the west flank of Mauna Loa, paired with a shift in earthquakes under Moku‘āweoweo were the two most notable long-term precursors to an eruption; however, neither provided an actionable timeline of an impending eruption. The decade of precursory seismicity prior to the 2022 eruption is similar to long precursory phases prior to eruptions in 1975 and 1984 and underscores the challenges of informing on eruption hazards in a timely manner. By combining seismic activity with other geophysical fields such as geodesy and gravity, we present a conceptual model to aid in the interpretation of future unrest at Mauna Loa.

Hawaii↗

Remote single-station seismic monitoring of the July–October 2022 earthquake swarm at Tau volcano, American Samoa

From July to October 2022, a non-eruptive volcanic earthquake swarm occurred within ~15 km of Taʻū Island, located in eastern American Samoa. Felt reports from local residents were the only available information about the swarm when it started, as American Samoa lacked a seismic monitoring network. We developed a consistent single-station catalog for the entire swarm, using seismic data from the nearest station IU.AFI, ~250 km away. We applied the EQTransformer deep-learning model (Mousavi et al., 2020), automatically picking Pn and Sn arrivals on IU.AFI continuous data. We retained only events with Sn-Pn times of 22.5–25 seconds, consistent with the expected locations based on felt reports, then detected smaller swarm events with subsequent template-matching. This single-station catalog characterized the swarm’s onset and escalation to peak activity before a multi-agency field response team installed a local seismic network in mid-August 2022. This permanent seismic network captured the swarm’s decline. EQTransformer identified short S-P times on the first two locally deployed seismometers, both Raspberry Shake sensors, to constrain the swarm’s distance from Taʻū Island. Modern seismological processing methods, combined with basic observations such as felt reports, can quickly contribute useful information during an earthquake response in a poorly monitored region.

American Somoa↗

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↗

Triggering the 2022 eruption of Mauna Loa

Distinguishing periods of intermittent unrest from the run-up to eruption is a major challenge at volcanoes around the globe. Comparing multidisciplinary monitoring data with mineral chemistry that records the physical and spatio-temporal evolution of magmas fundamentally advances our ability to forecast eruptions. The recent eruption of Mauna Loa, Earth’s largest active volcano, provides a unique opportunity to differentiate unrest from run-up and improve forecasting of future eruptions. After decades of intermittent seismic and geodetic activity over 38 years of repose, Mauna Loa began erupting on 27 November 2022. Here we present a multidisciplinary synthesis that tracks the spatio-temporal evolution of precursory activity by integrating mineral and melt chemistry, fluid inclusion barometry, numerical modeling of mineral zoning, syn-eruptive gas plume measurements, the distribution and frequency of earthquake hypocenters, seismic velocity changes, and ground deformation. These diverse data indicate that the eruption occurred following a 2-month period of sustained magma intrusion from depths of 3–5 km up to 1–2 km beneath the summit caldera, providing a new model of the plumbing system at this very high threat volcano. Careful correlation of both the geochemistry and instrumental monitoring data improves our ability to distinguish unrest from the run-up to eruption by providing deeper understanding of the both the monitoring data and the magmatic system—an approach that could be applied at other volcanic systems worldwide.

Hawaii↗

Olivine diffusion constrains months-scale magma transport within Kīlauea volcano’s summit reservoir system prior to the 2020 eruption

The unprecedented 2018 summit collapse at Kīlauea and subsequent 2020–2021 eruption within the newly deepened Halema‘uma‘u Crater provide an unparalleled opportunity to understand how collapse events impact a volcano’s shallow reservoir system and magmatic processes. Glass and olivine from tephra ejected by lava fountains and several explosions on 20–21 December, within a few hours of the 2020 eruption onset, yield information about pre-eruptive magma storage and transport. The olivine population is bimodal with zoned and non-zoned phenocrysts. Normally zoned olivine crystals with core compositions around Fo 88 have 30–50 μm wide Fo 82 overgrowth rims that have skeletal textures. Two skeletal xenocrysts (cores Fo 74 and Fo 81 ) are also reversely zoned up to Fo 82 rims. The crystal cores have trace element records of at least two cycles of growth and dissolution prior to the formation of the overgrowth rims. These rims and a separate population of non-zoned Fo 82 crystals are in Fe–Mg equilibrium with their host glass (average MgO of 6.9 ± 0.4 wt% (1σ), Mg# [Mg / (Mg + Fe 2+ )] of 0.57), which suggests undercooling after intrusion of magma to shallow levels in the plumbing system. In the years prior to the 2018 collapse, non-zoned Fo 81 olivine and slightly lower MgO glasses (6.8 wt%) reflected continuous mixing and compositional buffering of magma recharge into several km 3 of stored magma in the Halema‘uma‘u reservoir (1–2 km depth). The 2020 olivine crystals lack evidence of an intrusion mixing with resident shallow magma, indicating that magma transport occurred in a disrupted system, and/or it may not have significantly mixed with stored magma remaining in the Halema‘uma‘u reservoir after the events of 2018. Diffusion modeling of Fe–Mg exchange in the zoned 2020 olivine crystals yield timescales that are mostly 60 days prior to the eruption or less, which aligns well with 22–24 October 2020 and subsequent seismic swarms at Nāmakani Paio ~ 5 km west of Kīlauea’s summit caldera. This correlation indicates that magma intruding beneath the summit (volume accommodation, recorded by the olivine crystals) was expressed by tectonic earthquakes along the Ka‘ōiki fault zone (stress accommodation). The absence of precursory SO 2 within minutes prior to eruption also indicates that the 2020 December magma may have risen from 1 to 2 km depth to the surface in as little as 10 min.

Hawaii↗

Neotectonic and paleoseismic analysis of the northwest extent of Holocene surface deformation along the Meers Fault, Oklahoma

TheMeers fault (Oklahoma) is one of fewseismogenic structures with evidence for Holocene surface rupture in the stable continental region of North America. The 37-kilometer-long southeast section of the full 54-kilometer-long Meers fault is interpreted to be Holocene active. The 17-kilometer-long northwest section is considered Quaternary active, but not Holocene active.We reevaluate surface expression and earthquake timing of the northwest Meers fault to improve seismic source characterization.We use airborne light detection and ranging and historical stereopaired aerial photos to evaluate the fault scarp and local faultzone geomorphology. In the northwest, complex surface deformation includes fault splays, subtle monoclinal warping, and a minor change in fault strike. We interpret that the alongstrike transition from surface faulting on the southeastMeers fault to surface folding on the northwest Meers fault occurs at the lithologic contact between Permian Post Oak conglomerate and Hennessey shale. We excavated a paleoseismic trench to evaluate the timing of surface-deforming earthquakes on the northwest section of the fault. The excavation revealed weathered Permian Hennessey shale and an ∼1–2-meter-thick veneer of Holocene alluvial deposits that were progressively deformed during two surface-folding earthquakes likely related to blind fault rupture beneath the site. Repeated onlapping to overlapping stratigraphic sequences and associated unconformities are intimately related to folding events along the monocline. OxCal paleoearthquake age modeling indicates that earthquakes occurred 4704–3109 yr B.P. and 5955–4744 yr B.P., and that part of the northwest section of the Meers fault is Holocene active. We find the Holocene-active section of the Meers fault should be lengthened 6.1 km to the northwest, to a total Holocene-active fault length of 43 km. Empirical scaling relationships between surface rupture length and magnitude reveal that the fault could generate an Mw 7.0 earthquake.

Oklahoma↗

The susceptibility of Oklahoma’s basement to seismic reactivation

Recent widespread seismicity in Oklahoma is attributed to the reactivation of pre-existing, critically stressed and seismically unstable faults due to decades of wastewater injection. However, the structure and properties of the reactivated faults remain concealed by the sedimentary cover. Here, we explore the major ingredients needed to induce earthquakes in Oklahoma by characterizing basement faults in the field, in seismic surveys and via rock-mechanics experiments. Outcrop and satellite mapping reveal widespread fault and fracture systems with trends that display a marked similarity to the trends of recent earthquake lineaments. Our three-dimensional seismic analyses show steeply dipping basement-rooted faults that penetrate the overlying sedimentary sequences, representing pathways for wastewater migration. Experimental stability analysis indicates that Oklahoma’s basement rocks become seismically unstable at conditions relevant to the dominant hypocentral depths of the recent earthquakes. These analyses demonstrate that the geometry, structure and mechanical stability of Oklahoma’s basement make it critically susceptible to seismic reactivation.

Oklahoma↗

The 2018 rift eruption and summit collapse of Kilauea Volcano

In 2018, Kīlauea Volcano experienced its largest lower East Rift Zone (LERZ) eruption and caldera collapse in at least 200 years. After collapse of the Pu'u 'Ō'ō vent on 30 April, magma propagated downrift. Eruptive fissures opened in the LERZ on 3 May, eventually extending ~6.8 km. A 4 May earthquake (M6.9) produced ~5 m of fault slip. Lava erupted at rates exceeding 100 m3/s, eventually covering 35.5 km2. The summit magma system partially drained, producing minor explosions and near-daily collapses releasing energy equivalent to M4.7-M5.4 earthquakes. Activity declined rapidly on 4 August. Summit collapse and lava flow volume estimates are roughly equivalent--about 0.8 km3. Careful historical observation and monitoring of Kīlauea enabled successful forecasting of hazardous events.

Hawaii↗

Evolution of wear and friction along experimental faults

We investigate the evolution of wear and friction along experimental faults composed of solid rock blocks. This evolution is analyzed through shear experiments along five rock types, and the experiments were conducted in a rotary apparatus at slip velocities of 0.002&ndash;0.97 m/s, slip distances from a few millimeters to tens of meters, and normal stress of 0.25&ndash;6.9 MPa. The wear and friction measurements and fault surface observations revealed three evolution phases: A) An initial stage (slip distances <50 mm) of wear by failure of isolated asperities associated with roughening of the fault surface; B) a running-in stage of slip distances of 1&ndash;3 m with intense wear-rate, failure of many asperities, and simultaneous reduction of the friction coefficient and wear-rate; and C) a steady-state stage that initiates when the fault surface is covered by a gouge layer, and during which both wear-rate and friction coefficient maintain quasi-constant, low levels. While these evolution stages are clearly recognizable for experimental faults made from bare rock blocks, our analysis suggests that natural faults &ldquo;bypass&rdquo; the first two stages and slip at gouge-controlled steady-state conditions.

Pure and Applied Geophysics↗

Rapid acceleration leads to rapid weakening in earthquake-like laboratory experiments

After nucleation, a large earthquake propagates as an expanding rupture front along a fault. This front activates countless fault patches that slip by consuming energy stored in Earth’s crust. We simulated the slip of a fault patch by rapidly loading an experimental fault with energy stored in a spinning flywheel. The spontaneous evolution of strength, acceleration, and velocity indicates that our experiments are proxies of fault-patch behavior during earthquakes of moment magnitude (M w ) = 4 to 8. We show that seismically determined earthquake parameters (e.g., displacement, velocity, magnitude, or fracture energy) can be used to estimate the intensity of the energy release during an earthquake. Our experiments further indicate that high acceleration imposed by the earthquake’s rupture front quickens dynamic weakening by intense wear of the fault zone.

Science↗