Kilauea Research Articles

Magma fizz: Tremor during the Kīlauea summit reservoir decompression

Typical eruptions at Kīlauea volcano involve the evacuation of magma from the summit and/or south caldera reservoirs towards the East or Southwest rift zones. The reservoir drainage provokes the summit deflation, and on extreme occasions, such as in 2018, the summit caldera collapse. Systematically, seismic tremor, often with a particular multichromatic spectral signature characterized by frequency gliding, accompanies summit deflation episodes. In 2018, this type of continuous tremor accompanied the steady subsidence stage, whereas discrete earthquakes dominated the collapse stage. In this work, we aim to understand the source mechanism of the syn-deflation tremor of 2018. To locate the seismic source, we develop a novel machine-learning-based algorithm as an alternative to the amplitude source location technique. We use a large high-resolution catalog to resolve a composite amplitude decay function. Under these conditions, our method outperforms the traditional technique. We locate the tremor source 1 km below the eastern perimeter of the Halema‘uma‘u crater, which coincides with the position of the summit magma reservoir, as determined in many other studies. Furthermore, we model the seismic source as pressure oscillations driven by gas porous flow at the roof of the reservoir. In this model, gas accumulates temporarily in many gas pockets between the magma and the roof. Our modeling shows that the gas flux is responsible for the tremor amplitude modulations, whereas the gas pocket thickness controls the frequency variations. Beyond a critical point of depressurization, the magma cannot contribute further to the tremor oscillations via decompression-driven degassing, nor support the roof above it, resulting in rock failure. This work advances our understanding of magma-degassing dynamics and tremor generation at Kilauea volcano, and provides novel seismological techniques for volcano seismology monitoring and research.

Geomorphological evidence for volcano-tectonic deformation along the unstable western flank of Cumbre Vieja Volcano (La Palma)

In 2021, La Palma's southern volcanic complex Cumbre Vieja erupted for its longest period in historic times. Although the geological record shows no evidence for a collapse of Cumbre Vieja, ground deformation studies and field observations suggest that its western flank is moving seawards, following the direction of previous collapses of the island. To better estimate the hazard of a potential flank collapse of Cumbre Vieja, it is important to identify the lateral extent and depth of the mobile sector. Here, we analyse the volcano-tectonic deformation along Cumbre Vieja's western flank, based on geomorphological analysis of combined topographic and new ship-born bathymetric data as well as the analysis of shallow seismicity records associated with the 2021 eruption. In our interpretation, the shoreline-crossing Puerto Naos Ridge results from tectonic uplift accompanying transpressional deformation along the northern boundary of Cumbre Vieja's moving flank, therefore decoupling a stable sector in the north from the mobile sector farther south. The proposed moving sector is consistent in scale with previous ground deformation studies and documented flank collapses of structurally similar volcanoes. We present a workflow for semi-automatically detecting boundaries of unstable volcanic flanks based on morphological changes captured in digital elevation data. The method correctly delineated the known boundaries of the unstable flanks of Mt. Etna and Kilauea volcanoes. The ability to constrain potential boundaries of unstable volcanic flanks should inform the planning of future geophysical and geodetic campaigns aiming to identify precursory signals of potential flank failures.

Correction: Density structure of Kilauea volcano: implications for magma storage and transport

Correction: Density structure of Kilauea volcano: implications for magma storage and transport

From the dark side of paradise: a new natural replication of cave planthopper evolution from Hawaiian lava tubes (Hemiptera: Fulgoromorpha: Cixiidae)

Abstract The Hawaiian Islands are known to harbour a rich and diverse fauna of troglobionts (obligate subterranean species). To date, 74 obligate cavernicolous arthropod species have been documented from across the main Hawaiian islands, the majority of which were from Hawaiʻi Island, and mostly from lava tubes of Kilauea volcano, the youngest volcano on the island. A recent bioinventory of the Kipuka Kanohina lava tube system on the south-western side of Mauna Loa volcano revealed the existence of previously unknown cave-adapted species. Among them is the first cave-adapted species of the planthopper genus Iolania, Iolania frankanstonei Hoch & Porter sp. nov. Morphological and molecular data suggest that the species is closely related to the epigean (i.e. surface-dwelling) species Iolania perkinsi, which occurs in surface environments on Hawaiʻi Island. Thus, parapatric speciation is assumed, further corroborating the assumption that adaptive shifts are the major evolutionary patterns underlying the evolution of troglobionts on young oceanic islands.

Bridging Supervised and Unsupervised Learning to Build Volcano Seismicity Classifiers at Kilauea Volcano, Hawaii

Abstract Real-time classification of volcano seismicity could become a useful component in volcanic monitoring. Supervised learning provides a powerful means to achieve this but often requires a large amount of manually labeled data. Here, we build supervised learning models to discriminate volcano tectonic events (VTs), long-period events (LPs), and hybrid events in Kilauea by training with pseudolabels from unsupervised clustering. We test three different supervised models, and all of them achieve >93% accuracy. We apply the model ensemble to the six-day seismicity during the eruption in 2018 and show that they were mainly VTs (62%), in comparison with the dominance of LPs prior to the eruption (68%). The success of our method is aided by the accuracy of the majority of pseudolabels and the consistency of the three models’ performance. Using Shapley additive explanations, we show that the frequency contents at 1–4 Hz are the most important to differentiate volcano seismicity types. This work, together with our previous clustering analysis, provides an example of bridging unsupervised and supervised learning to construct potential real-time seismic classifiers from scratch.

Topography‐Incorporated Adjoint‐State Surface Wave Traveltime Tomography: Method and a Case Study in Hawaii

AbstractIn this study we recast surface wave traveltime tomography as an inverse problem constrained by an eikonal equation and solve it using the efficient adjoint‐state method. Specifically, recognizing that large topographic variations and high surface wave frequencies can make the topographic effect too significant to ignore, we employ an elliptically anisotropic eikonal equation to describe the traveltime fields of surface waves on undulated topography. The sensitivity kernel of the traveltime objective function with respect to shear wave velocity is derived using the adjoint‐state method. As a result, the newly developed method is inherently applicable to any study regions, whether with or without significant topographic variations. Hawaii is one of the most seismically and magmatically active regions. However, its significant topographic variations have made it less accurate to investigate using conventional surface wave traveltime tomography methods. To tackle this problem, we applied our new method to invert ambient noise Rayleigh wave phase traveltimes and construct a 3D shear wave velocity model. Our results reveal features that are consistent with geological structures and previous tomography results, including high velocities below Mauna Loa Volcano and Kilauea Volcano, and low velocities beneath the Hilina Fault Zone. Additionally, our model reveals a high‐velocity anomaly to the South of Hualalai's summit, which may be related to a buried rift zone. Our findings further demonstrate that including topography can lead to a correction of up to 0.8% in the shear wave velocity model of Hawaii, an island spanning approximately 100 km with volcanoes reaching elevations exceeding 4 km.

In-situ monitoring of 3He/4He in summit gases of Kilauea Volcano (Hawaii) prior to the 2020 eruption

We present He isotope (3He/4He) data from a fumarole and near-ground gases measured in-situ at the Sulfur Banks solfatara field at the summit of Kilauea Volcano, Hawaii. We used a field-deployable mass-spectrometer-based system: the Helium Isotope Monitor (HIM) previously described in McMurtry et al. (2019a, b). The in-situ instrument was deployed using solar power for the first time and results were ground-truthed against data determined using conventional gas analytical and noble gas mass spectrometry techniques. The HIM instrumentation, associated Vent Gas Purification System (VGPS), and related sampling equipment and strategy are described. Cloudy and rainy weather conditions hampered the deployment, which was reorganized to reduce power loads and resulted in less sampling than planned. Nevertheless, we obtained daily sampling of the volcanic vent gas. Results from the Old Well fumarole indicate a ~ 2 RA increase in 3He/4He on the day of the December 20th, 2020 eruption of nearby Halema‘uma‘u Crater, reaching 17.0 RA using the in-situ instrument and 16.0 ± 0.67 RA using conventional techniques. This finding suggests that a new 3He-enriched magma source is driving the current, ongoing eruption phase of Kilauea and, if so, confirms that the deep summit caldera fault system that hosts the Sulfur Banks field is connected to the Halema‘uma‘u Crater magmatic system. Overall, these findings illustrate how time-series helium isotope data, which are well established by ongoing discrete monitoring at low temporal resolution, can help forecast forthcoming eruptive events that may not be foreseen by other volcanic monitoring methods.

Carbon isotope composition of basalts from Loihi Seamount: Primordial or recycled carbon in the Hawaiian mantle plume?

We analyzed the carbon isotope composition of vesicle CO2, plus He isotopes and He and CO2 concentrations in the vesicle (vapor) and glass (melt) phases of 37 submarine basalts from the summit, north and south rifts, and east flank of Loihi Seamount. Tholeiites and transitional basalts lie in a narrow range of vesicle δ13C = −4.6 to −0.9‰, while alkali basalts range from −7.2 to −2.1‰. Calculated total (vesicle+glass) δ13C for the majority of the basalts ranges from −6 to −2‰ assuming the vapor-melt fractionation factor Δ (= δvapor - δmelt) is +2 to +4‰ as measured in basaltic systems. This relatively narrow range of δ13C resembles mantle source values deduced from gas-rich mid-ocean ridge basalts and basalts from Iceland, and for Kilauea volcano deduced from its fumarole gas. However, this similarity presents a conundrum because Loihi basalts have degassed >97% of their initial CO2 as deduced from CO2 - Ba systematics and crystal fractionation modeling.Loihi parental magma (MgO=18 wt.%) had initial CO2 concentrations of 0.6 to 1.9 wt.%. Most tholeiitic and transitional basalts appear to have followed a quasi closed-system degassing history. Correcting for this degassing indicates the median δ13C for Loihi undegassed parental magmas is −1.5‰ and the δ13C range is −4 to +1‰. Estimates of this δ13C range are only weakly dependent on the choice of Δ and initial [CO2] in closed-system degassing scenarios. The Loihi mantle plume source is therefore characterized by δ13C values that are higher than the range of −6 to −4‰ prevalent in mantle-derived basalts and many diamonds. This could be due to primordial carbon isotope heterogeneity in Earth's mantle, exchange of carbon at the core-mantle boundary between ultra-low velocity zone silicates and the metallic liquid outer core, or to the presence of a small fraction (<1% by mass) of surficial carbonate that was tectonically recycled to the Hawaiian plume source region. Currently, an origin from recycled carbonate has the most supporting evidence.

Unraveling the Magnetic Signal of Individual Grains in a Hawaiian Lava Using Micromagnetic Tomography

AbstractMicromagnetic Tomography (MMT) is a new technique that allows the determination of magnetic moments of individual grains in volcanic rocks. Current MMT studies either showed that it is possible to obtain magnetic moments of relatively small numbers of grains in ideal sample material or provided important theoretical advances in MMT inversion theory and/or its statistical framework. Here, we present a large‐scale application of MMT on a sample from the 1907‐flow from Hawaii's Kilauea volcano producing magnetic moments of 1,646 grains. We produced 261,305 magnetic moments in total for these 1,646 grains, an increase of three orders of magnitude compared to earlier studies to assess the robustness of the MMT results, and a major step toward the number of grains that is necessary for paleomagnetic applications of MMT. Furthermore, we show that the recently proposed signal strength ratio is a powerful tool to scrutinize and select MMT results. Despite this progress, still only relatively large iron‐oxide grains with diameters &gt;1.5–2 μm can be reliably resolved, impeding a reliable paleomagnetic interpretation. To determine the magnetic moments of smaller (&lt;1 μm) grains that may exhibit pseudo‐single domain behavior and are therefore better paleomagnetic recorders, the resolution of the microcomputed tomography and magnetic scans necessary for MMT must be improved. Therefore, it is necessary to reduce the sample size in future MMT studies. Nevertheless, our study is an important step toward making MMT a useful paleomagnetic and rock‐magnetic technique.

Hydrovolcanic Explosions at the Lava Ocean Entry of the 2018 Kilauea Eruption Recorded by Ocean-Bottom Seismometers

AbstractFrom the beginning of May 2018, the Kilauea Volcano on the island of Hawaii experienced its largest eruption in 200 yr followed by a period of unrest for months. Because hot molten lava entered the ocean from the ocean-entry point near the lower East Rift Zone, the lava–water interaction led to explosions. Some explosions were near the water surface and ejected fragments of lava, also known as lava bombs. In the early morning on 16 July 2018, one of those lava bombs, which was almost the size of a basketball, hit a sightseeing boat and injured 23 people. In this study, we analyzed the hydrophone data recorded from July to mid-September by ocean-bottom seismometers (OBSs) deployed offshore near the ocean entry point to identify and locate the hydroacoustic signals of the lava–water explosions. Acoustic signals of hydrovolcanic explosions are characterized by a short duration (less than a few seconds) and a broad frequency range (at least up to 100 Hz). To automate event detection, a short-term average versus long-term average method was applied to the complete dataset. Approximately 4300 events were detected and located near the coastline and further used to prepare a catalog. The distribution of the lava–water explosions is consistent with the pattern of the offshore lava delta formed during the 2018 eruption. Identifying such hydroacoustic signals recorded by OBSs may provide new avenues of research using various seismoacoustic events associated with volcanic eruptions.

Effects of Volcanic Aerosols on the Genesis of Tropical Cyclone Wukong (2018)

AbstractThe short‐range impacts of volcanic aerosols from the eruption of the Kilauea volcano on the genesis of tropical cyclone (TC) Wukong in 2018 are investigated by numerical simulations, as distinct from the studies on the aerosol climate effects. The aerosol−radiation and aerosol−cloud effects are isolated through sensitivity experiments. For the aerosol−radiation effect, the radiative cooling at the low levels and heating above by aerosols southeast of Wukong's periphery induce stable sinking and decrease latent heating in the environment far from the storm center. Therefore, the sea‐level pressure increases there, enhancing low‐level radial inflow to the storm, inner‐core convection and the release of latent heat, which is conducive to TC genesis. For the aerosol−cloud effect, volcanic aerosols in storm circulation can be converted into cloud condensation nuclei, invigorating inner‐core convection, and subsequently favoring TC genesis. However, the nonlinear effect resulting from the interaction between aerosol−radiation and aerosol−cloud effects is not conducive to TC genesis, because it weakens the aerosol−cloud effect through the deposition of aerosols and also limits the aerosol−radiation effect through the increasing of peripheral convection.

First Implementation of a Normalized Hotspot Index on Himawari-8 and GOES-R Data for the Active Volcanoes Monitoring: Results and Future Developments

The Advanced Himawari Imager (AHI) and Advanced Baseline Imager (ABI), respectively aboard Himawari-8 and GOES-R geostationary satellites, are two important instruments for the near-real time monitoring of active volcanoes in the Eastern Asia/Western Pacific region and the Pacific Ring of Fire. In this work, we use for the first time AHI and ABI data, at 10 min temporal resolution, to assess the behavior of a Normalized Hotspot Index (NHI) in presence of active lava flows/lakes, at Krakatau (Indonesia), Ambrym (Vanuatu) and Kilauea (HI, USA) volcanoes. Results show that the index, which is used operationally to map hot targets through the Multispectral Instrument (MSI) and the Operational Land Imager (OLI), is sensitive to high-temperature features even when short-wave infrared (SWIR) data at 2 km spatial resolution are analyzed. On the other hand, thresholds should be tailored to those data to better discriminate thermal anomalies from the background in daylight conditions. In this context, the multi-temporal analysis of NHI may enable an efficient identification of high-temperature targets without using fixed thresholds. This approach could be exported to SWIR data from the Flexible Combined Imager (FCI) instrument aboard the next Meteosat Third Generation (MTG) satellites.

Distributed Scatterer Processing Based on Binary Partition Trees with Multi-Baseline PolInSAR Data

Distributed scatterers (DSs) are necessary to increase point density in multi-temporal InSAR (MT-InSAR) monitoring. The identification of homogeneous pixels (HPs) is the first and key step for DS processing to overcome the low signal-to-noise ratio condition. Since multi-polarization data are good at describing geometrical structures and dielectric properties of ground objects, they have been applied for HP identification. However, polarimetric information is not enough for identifying areas with similar ground objects but different deformation. We propose a novel DS preprocessing algorithm based on polarimetric interferometric homogeneous pixel (PIHP) identification. Firstly, a novel Polarimetric InSAR (PolInSAR) similarity that combines polarimetric intensity, interferometric coherence, and phase is proposed, which is readily available in multi-baseline and multi-polarization data and flexible by controlling weighting factors. Secondly, based on the binary partition tree (BPT) framework, object-orientated multi-scale PIHP identification is achieved, which is suitable for complex deformation scenes. Tested with simulated quad-polarization data, our method shows improvement in phase quality and point density, especially in the deformed areas, compared with the traditional HP identification method based on the polarimetric homogeneity (PolHom) test and the method with ground object type map. Tested with 30 quad-polarization Radarsat-2 images over Kilauea Volcano, the point density of our method is three times higher than that of the PolHom test in vegetation areas. Our method is proven to be more sensitive and mechanically more advanced to homogeneous pixels identification than the traditional ones, which is helpful for phase optimization, spatial enlargement of monitoring points, and stability of the MT-InSAR algorithm.

Tides and Volcanoes: A Historical Perspective

In this paper, we examine the origins and the history of the hypothesis for an influence of tidal forces on volcanic activity. We believe that exploring this subject through a historical perspective may help geoscientists gain new insights in a field of research so closely connected with the contemporary scientific debate and often erroneously considered as a totally separated niche topic. The idea of an influence of the Moon and Sun on magmatic processes dates back to the Hellenistic world. However, it was only since the late 19thcentury, with the establishment of volcano observatories at Mt. Etna and Vesuvius allowing a systematic collection of observations with modern methods, that the “tidal controversy” opened one of the longest and most important debates in Earth Science. At the beginning of the 20thcentury, the controversy assumed a much more general significance, as the debate around the tidal influence on volcanism developed around the formulation of the first modern theories on the origins of volcanism, the structure of the Earth’s interior and the mechanisms for continental drift. During the same period, the first experimental evidence for the existence of the Earth tides by Hecker (Beobachtungen an Horizontalpendeln über die Deformation des Erdkörpers unter dem Einfluss von Sonne und MondVeröffentlichung des Königl, 1907, 32), and the Chamberlin–Moulton planetesimal hypothesis (proposed in 1905 by geologist Thomas Chrowder Chamberlin and astronomer Forest Ray Moulton) about the “tidal” origin of the Solar System, influenced and stimulated new researches on volcano-tides interactions, such as the first description of the “lava tide” at the Kilauea volcano by Thomas Augustus Jaggar in 1924. Surprisingly, this phase of gradual acceptance of the tidal hypothesis was followed by a period of lapse between 1930 to late 1960. A new era of stimulating and interesting speculations opened at the beginning of the seventies of the 20thcentury thanks to the discovery of the moonquakes revealed by the Apollo Lunar Surface Experiment Package. A few years later, in 1979, the intense volcanism on the Jupiter’s moon Io, discovered by the Voyager 1 mission, was explained by the tidal heating produced by the Io’s orbital eccentricity. In the last part of the paper, we discuss the major advances over the last decades and the new frontiers of this research topic, which traditionally bears on interdisciplinary contributions (e.g., from geosciences, physics, astronomy). We conclude that the present-day debate around the environmental crisis, characterized by a large collection of interconnected variables, stimulated a new field of research around the complex mechanisms of mutual interactions among orbital factors, Milankovitch Cycles, climate changes and volcanism.

A Compound Faulting Model for the 1975 Kalapana, Hawaii, Earthquake, Landslide, and Tsunami

AbstractThe Kalapana, Hawaii, MW 7.7 earthquake on November 29, 1975 generated a local tsunami with at least 14.3 m runup on the southeast shore of Hawaii Island adjacent to Kilauea Volcano. This was the largest locally generated tsunami since the great 1868 Ka'u earthquake located along‐shore to the southwest. Well‐recorded tide gauge and runup observations provide a key benchmark for studies of statewide tsunami hazards from actively deforming southeast Hawaii Island. However, the source process of the earthquake remains controversial, with coastal landsliding and/or offshore normal or thrust faulting mechanisms having been proposed to reconcile features of seismic, geodetic, and tsunami observations. We utilize these diverse observations for the 1975 Kalapana earthquake to deduce a compound faulting model that accounts for the overall tsunamigenesis, involving both landslide block faulting along the shore and slip on the island basal décollement. Thrust slip of 4.5–8.0 m on the offshore décollement produces moderate near‐field runup but controls the far‐field tsunami. The slip distribution implies that residual strain energy was available for the May 4, 2018 MW 7.2 thrust earthquake during the Kilauea‐East Rift Zone eruption. Local faulting below land contributes to geodetic and seismic observations, but is non‐tsunamigenic and not considered. Slip of 4–10 m on landslide‐like faults, which extend from the Hilina Fault Zone scarp to offshore shallowly dipping faults reaching near the seafloor, triples the near‐field tsunami runup. This compound model clarifies the roles of the faulting components in assessing tsunami hazards for the Hawaiian Islands.

Subdivision of Seismicity Beneath the Summit Region of Kilauea Volcano: Implications for the Preparation Process of the 2018 Eruption

AbstractLong‐period (LP), hybrid, and volcano‐tectonic (VT) seismicity are important indicators for tracking the evolution of volcanic processes. Here, we propose an unsupervised learning method to classify 5,949 seismic events in Kilauea volcano, Hawai'i, during a 4‐month period before the collapse of Pu'u’ O'o on April 30, 2018. The LPs and hybrids exhibit three major episodes, which progressively intensified and had increasing shallow events toward the eruption. The most intense episode starting 3 weeks before eruption coincides with changes in near‐caldera deformation and lava lake elevation in Halema'uma'u, serving as possible immediate precursors. However, the first two episodes imply magma migration was already active months prior to the eruption. The spatiotemporal patterns of abundant hybrids reveal that they are associated with magma movement but mixed with shear‐failure or near‐surface resonance. Our results provide useful constraints on the magmatic processes in the preparation phase of the Kilauea eruption in 2018.

Incorporation of volcanic SO&amp;lt;sub&amp;gt;2&amp;lt;/sub&amp;gt; emissions in the Hemispheric CMAQ (H-CMAQ) version 5.2 modeling system and assessing their impacts on sulfate aerosol over the Northern Hemisphere

Abstract. The state-of-the-science Community Multiscale Air Quality (CMAQ) Modeling System has recently been extended for hemispheric-scale modeling applications (referred to as H-CMAQ). In this study, satellite-constrained estimation of the degassing SO2 emissions from 50 volcanoes over the Northern Hemisphere is incorporated into H-CMAQ, and their impact on tropospheric sulfate aerosol (SO42-) levels is assessed for 2010. The volcanic degassing improves predictions of observations from the Acid Deposition Monitoring Network in East Asia (EANET), the United States Clean Air Status and Trends Network (CASTNET), and the United States Integrated Monitoring of Protected Visual Environments (IMPROVE). Over Asia, the increased SO42- concentrations were seen to correspond to the locations of volcanoes, especially over Japan and Indonesia. Over the USA, the largest impacts that occurred over the central Pacific were caused by including the Hawaiian Kilauea volcano, while the impacts on the continental USA were limited to the western portion during summertime. The emissions of the Soufrière Hills volcano located on the island of Montserrat in the Caribbean Sea affected the southeastern USA during the winter season. The analysis at specific sites in Hawaii and Florida also confirmed improvements in regional performance for modeled SO42- by including volcanoes SO2 emissions. At the edge of the western USA, monthly averaged SO42- enhancements greater than 0.1 µg m−3 were noted within the boundary layer (defined as surface to 750 hPa) during June–September. Investigating the change on SO42- concentration throughout the free troposphere revealed that although the considered volcanic SO2 emissions occurred at or below the middle of free troposphere (500 hPa), compared to the simulation without the volcanic source, SO42- enhancements of more than 10 % were detected up to the top of the free troposphere (250 hPa). Our model simulations and comparisons with measurements across the Northern Hemisphere indicate that the degassing volcanic SO2 emissions are an important source and should be considered in air quality model simulations assessing background SO42- levels and their source attribution.

Effect of volcanic emissions on clouds during the 2008 and 2018 Kilauea degassing events

Abstract. Volcanic eruptions in otherwise clean environments are “natural experiments” wherein the effects of aerosol emissions on clouds and climate can be partitioned from meteorological variability and anthropogenic activities. In this work, we combined satellite retrievals, reanalysis products, and atmospheric modeling to analyze the mechanisms of aerosol–cloud interactions during two degassing events at the Kilauea volcano in 2008 and 2018. The eruptive nature of the 2008 and 2018 degassing events was distinct from long-term volcanic activity for Kilauea. Although previous studies assessed the modulation of cloud properties from the 2008 event, this is the first time such an analysis has been reported for the 2018 event and that multiple degassing events have been analyzed and compared at this location. Both events resulted in significant changes in cloud effective radius and cloud droplet number concentration that were decoupled from local meteorology and in line with an enhanced cloud albedo. However, it is likely that the effects of volcanic emissions on liquid water path and cloud fraction were largely offset by meteorological variability. Comparison of cloud anomalies between the two events suggested a threshold response of aerosol–cloud interactions to overcome meteorological effects, largely controlled by aerosol loading. In both events, the ingestion of aerosols within convective parcels enhanced the detrainment of condensate in the upper troposphere, resulting in deeper clouds than observed under pristine conditions. Accounting for ice nucleation on ash particles led to enhanced ice crystal concentrations at cirrus levels and a slight decrease in ice water content, improving the correlation of the model results with the satellite retrievals. Overall, aerosol loading, plume characteristics, and meteorology contributed to changes in cloud properties during the Kilauea degassing events.

Spatiotemporal Clustering of Seismicity During the 2018 Kilauea Volcanic Eruption

AbstractThe intense volcanic and seismic activity at Kilauea volcano, Hawaii in 2018 provides an opportunity to understand the characteristics and driving forces of volcano earthquake clusters. Nearest neighbor distance analysis of a high‐resolution earthquake catalog reveals two modes of clustering: M1 and M2. M1 clusters focus on events within short space and time windows, and are concentrated within isolated fault patches. Examination of three structures shows a gradual increase of source radius (assuming a constant stress drop) and spatial migration of seismicity, likely representing aseismic growth of peripheral faults. M2 clusters focus on temporal clustering of all seismicity within the caldera and most clusters represent complete seismicity cycles between major collapse events. Beginning in late June, the M2 cluster duration gradually increased. The event numbers show a strong linear relationship with the tilt rates on the transverse component, suggesting earthquake rates are controlled by the deformation rate within the caldera.

Three-dimensional strain descriptors at the Earth’s surface through 3D retrieved co-event displacement fields of differential interferometric synthetic aperture radar

In deformation analysis, invariant strain descriptors that are independent of the coordinate system are significant. An approach is proposed to determine the strain tensor and invariant strain quantities through the 3D retrieved co-event displacement of differential interferometric synthetic aperture radar (DInSAR). An optimal cell count is selected based on the proposed rule-of-thumb formula. This set of coherent cells with prescribed acceptable minimum value coherence is distributed on the domain through a quasi-random number generator. The recursive moving least squares method is applied to estimate the components of strain and differential rotation tensors. In practice, this method is not limited to the homogeneous strain field hypothesis. The objective here is to propose an approach to determine the local and regional strain patterns in SAR’s image resolution details. This approach is tested on both synthetic and real data. The real data consist of 3D displacements of 19 global navigation satellite system stations and SAR images of Kilauea Volcano eruption in July 2007. Results indicate that the areas with maximum shear strain and maximum dilatation are located exactly near the east rift zone of the volcano.