Southwest Rift Research Articles

Rift Zone Architecture and Inflation‐Driven Seismicity of Mauna Loa Volcano

AbstractThe 2022 eruption at Mauna Loa, Hawai'i, marked the first extrusive activity from the volcano after 38 years of quiescence. The eruption was preceded by several years of seismic unrest in the vicinity of the volcano's summit. Characterizing the structure and dynamics of seismogenic features within Mauna Loa during this pre‐eruptive interval may provide insights into how pre‐ and co‐eruptive processes manifest seismically at the volcano. In particular, the extent to which seismicity may be used to forecast the location and timing of future eruptions is unclear. To address these questions, we construct a catalog of relocated seismicity on Mauna Loa spanning 2011–2023. Our earthquake locations image complex, sub‐kilometer‐scale seismogenic structures in the caldera and southwest rift zone. We additionally identify a set of streaks of seismicity in the volcano's northwest flank that are radially oriented about the summit. Using a rate‐and‐state friction model for earthquake occurrences, we demonstrate that the seismicity rate in this region can be modeled as a function of the stressing history caused by magma accumulation beneath the summit. Finally, we observe a mid‐2019 step change in the seismicity rate in the Ka'oiki region that may have altered the stress state of the northeast rift zone in the three years before the eruption. Our observations provide a framework for interpreting future seismic unrest at Mauna Loa.

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.

Tracking magma pathways and surface faulting in the Southwest Rift Zone and the Koaʻe fault system (Kīlauea volcano, Hawai ‘i) using photogrammetry and structural observations

Volcanic islands are often subject to flank instability, resulting from a combination of magmatic intrusions along rift zones and gravitational spreading causing extensional faulting at the surface. Here, we study the Koaʻe fault system (KFS), located south of the summit caldera of Kīlauea volcano in Hawaiʻi, one of the most active volcanoes on Earth, prone to active faulting, episodic dike intrusions, and flank instability. Two rift zones and the KFS are major structures controlling volcanic flank instability and magma propagation. Although several magmatic intrusions occurred over the KFS, the link between these faults, two nearby rift zones and the flank instability, is still poorly studied. To better characterize the KFS and its structural linkage with the surrounding fault and rift zones, we performed a detailed structural analysis of the extensional fault system, coupled with a helicopter photogrammetric survey, covering part of the south flank of Kīlauea. We generated a high-resolution DEM (~ 8 cm) and orthomosaic (~ 4 cm) to map the fracture field in detail. We also collected ~ 1000 ground structural measurements of extensional fractures during our three field missions (2019, 2022, and 2023). We observed many small, interconnected grabens, monoclines, rollover structures, and en-echelon fractures that were in part previously undocumented. We estimate the cumulative displacement rate across the KFS during the last 600 ~ 700 years and found a decrease toward the west of the horizontal component from 2 to 6 cm per year, consistent with GNSS data. Integrating morphology observations, fault mapping, and kinematic measurements, we propose a new kinematic model of the upper part of the Kīlauea’s south flank, suggesting a clockwise rotation and a translation of a triangular wedge. This wedge is bordered by the extensional structures (ERZ, SWRZ, and the KFS), largely influenced by gravitational spreading. These findings illustrate a structural linkage between the two rift zones and the KFS, the latter being episodically affected by dike intrusions.

Dynamics of magma mixing and magma mobilisation beneath Mauna Loa—insights from the 1950 AD Southwest Rift Zone eruption

Eruptions from Mauna Loa’s Southwest Rift Zone (SWRZ) pose a significant threat to nearby communities due to high eruption rates and steep slopes resulting in little time for evacuation. Despite the large body of research done on Mauna Loa, knowledge of the timing and duration of magma residence and transfer through its internal plumbing system is still poorly constrained. This study presents a first quantitative look at thermochemical conditions and timescales of potentially deep storage and disaggregation of magmatic mush during the run-up to the voluminous 1950 AD SWRZ eruption. Details of heterogeneous compositions and textures of the macrocryst and glomerocryst cargo in 1950 AD lavas suggest magma mixing and crystal recycling along the entire plumbing system. Furthermore, the crystal cargo contains evidence for the direct interaction between primitive, deeply stored magma and pockets of more evolved magma stored at shallow to intermediate depths. An enigmatic attribute of 1950 near-vent lava is the near-ubiquitous presence of subhedral, unreacted Mg-rich orthopyroxene phenocrysts (Mg#>80). Phase relations of Mauna Loa olivine-tholeiite indicate that orthopyroxene joins olivine as a primary phase at pressures higher than 0.6 GPa. Coexisting Mg-rich olivine and orthopyroxene and the occurrence of harzburgitic (olivine-orthopyroxene) glomerocrysts provide evidence for cognate crystallisation at near-Moho (~ 18 km) depths (Thornber and Trusdell 2008). Petrogenetically diverse populations of glomerocrysts and macrocrysts alongside evidence of multilevel magma storage indicate a network of ephemeral and possibly interconnected magma pockets from near-Moho depths to the upper/mid-crust. Fe-Mg diffusion chronometry applied to 1950 AD olivine populations implies rapid mobilisation and transport of large volumes of magma (376×106 m3) from near-Moho storage to the surface within less than 8 months, with little residence time (~ 2 weeks) in the shallow (3–5 km) plumbing system.

Chemistry and petrography of early 19th century basaltic andesites and basalts from the Kamakaiʻa Hills in the Southwest Rift Zone of Kīlauea volcano, Hawaiʻi

Kīlauea is a frequently active, open-system volcano on the Island of Hawaiʻi known for erupting olivine-dominated tholeiitic basalt compositions. On rare occasions it erupts more differentiated magmas (<1% of erupted volume), such as basaltic andesites and andesites, from its rift zones. These differentiated magmas offer an opportunity to understand better the petrology, magma storage, magma mixing, and eruptive triggers that occur in Kīlauea's rift zone reservoirs. This study focuses on an eruption from the Southwest Rift Zone of Kīlauea, which is dominantly basaltic andesite with subordinate basalt. This eruption originated at the Kamakaiʻa Hills during the early 19th century and has two eruptive phases: 1) an early ‘a‘ā phase that is primarily exposed in the eastern part of the flow field, with minor western lobes, and 2) a late pāhoehoe phase that makes up most of the western part of the flow field. The early ‘a‘ā phase covers at least 5.8 km2 with an erupted volume of ∼150 × 106 m3 and consists of uniform composition basaltic andesites with 3.72–4.15 wt% MgO over its ∼7 km flow length. The late pāhoehoe phase reached >10 km from its vent, covers an area of ∼7.1 km2, has a volume of ∼100 × 106 m3, and initially erupted basaltic andesite near its vent (4.50–5.64 wt% MgO extending to 3.8 km from vent) with channel and tube-fed basalt (6.21–12.38 wt% MgO sampled at >3.8 km from vent) emplaced during its waning stages. Most Kamakaiʻa Hills lavas are crystal-poor, containing ≤1.5% glomerocrysts and individual phenocrysts of plagioclase + clinopyroxene + FeTi oxides ± orthopyroxene, as well as olivine in lavas with >6 wt% MgO.Major-oxide and trace-element concentrations throughout the Kamakaiʻa Hills lavas demonstrate the involvement of three distinct magmatic processes. First, the basaltic andesites of the early ‘a‘ā phase are the products of fractionation of plagioclase + clinopyroxene + FeTi oxides ± orthopyroxene, indicative of magmas that have been stored in rift zone reservoirs for decades or longer. Second, the near-vent (within ∼400 m of vent) basaltic andesites of the late pāhoehoe phase yield chemical concentrations that indicate magma mixing with a more differentiated magma (of a similar evolved composition to basaltic andesites at ∼55–56 wt% SiO2 and ∼3.4–4.1 wt% MgO that erupted in the lower East Rift Zone in 2018). Third, the progressively more mafic magma (containing olivine + plagioclase + clinopyroxene) that continued to erupt throughout the waning stages of activity suggests an eruptive triggering process whereby an intruding summit or uprift reservoir basalt overpressurized and forced out the stored, differentiated magma of the Kamakaiʻa Hills rift zone reservoir.

Pre-eruptive timescales from the historical Hapaimamo eruption at Mauna Loa, Hawai‘i

The 240 yr B.P. Hapaimamo eruption occurred on the southwest rift zone (SWRZ) of Mauna Loa. A large tephra deposit and an extensive lava flow field were produced that extend from 1880 m above sea level down to the southwest coast. Detailed petrological study of olivine crystals from both deposit types are used to determine the timescales of crystal disaggregation from a cumulate mush zone. Over 70% of studied olivines in both, lava flow and tephra samples, exhibit normal zoning. Tephra and flow olivine core compositions are between Fo89–84 and Fo89–81. Olivine rims in the tephra have compositions of Fo86–81, whereas rims in the flow samples extend to lower values within the range Fo83–71. The remaining tephra olivines are unzoned and fall within the same core compositional range. Using FeMg interdiffusion chronometry, we calculate tephra olivine diffusion timescales from 13 to 1600 days (days to years), with over 70% of timescales <100 days, and lava flow olivine diffusion timescales from 25 to 1110 days with 65% of timescales >100 days. Zoning widths were defined as narrow (type I) or broad (type II) and together with the diversity of timescales indicate an ongoing process of crystal mush disaggregation into the erupted host melt.Comparison of tephra and flow diffusion timescales show consistently longer timescales across the lava flow population, as expected due to continued diffusive re-equilibration. Our results suggest that where textural evidence allows, and tephra timescales can be compared, olivines from lava flows can be used to provide temporal information for magmatic processes within the subsurface and the lava flow. Comparison of published elemental diffusion timescales between Mauna Loa and neighbouring Kīlauea suggest that magma mobilisation and transfer processes occur over similar timescales within the plumbing systems of both Mauna Loa and Kīlauea.

Magnetotelluric Investigations of the Kīlauea Volcano, Hawaii

AbstractIn 2002 and 2003 a collaborative effort was undertaken between Lawrence Berkeley National Laboratory, Sandia National Laboratories, the U.S. Geological Survey (USGS) Menlo Park, the USGS Hawaiian Volcano Observatory, and Electromagnetic Instruments Inc. to study the Kīlauea volcano in Hawaii using the magnetotelluric (MT) technique. The work was motivated by a desire to improve understanding of the magma reservoirs and conduits within Kīlauea and the East and Southwest Rift zones, which has implications for understanding Kīlauea's plumbing system. An improved understanding of the rift zones has implications in understanding large‐scale landslides that are generated in the Hilina Slump, which produce significant impacts on coastal communities. Up to eight stations operated simultaneously, with multiple remote reference sites, and data were processed using multi‐station robust processing techniques. In total, data were acquired at 70 sites over the Southwest and East rift zones. Good to excellent quality data were obtained even in the harshest conditions, such as those encountered on the fresh lava flows of the East Rift Zone, where electrical contact resistances are on the order of 100 kΩ. A three‐dimensional (3D) MT model study was done to guide interpretation of the observed MT measurements. Synthetic modeling demonstrates that conductive bodies in the upper 3 km can be spatially resolved where MT station sampling is good. Resistivity anomalies in the 3D inversions have a high degree of spatial correlation with previously published seismic velocity anomalies beneath Kīlauea. Melt fractions between 0.096 and 0.117 are calculated for the Kīlauea and Puʻuʻōʻō low resistivity anomalies, respectively.

Iron isotope systematics during igneous differentiation in lavas from Kīlauea and Mauna Loa, Hawai'i

The Big Island of Hawai'i is currently the most active site of magmatic activity in the Hawaiian chain of volcanoes that are associated with a mantle plume in the Pacific Ocean. Mantle source variability and magmatic processes have both been proposed to significantly affect lava chemistry on the Big Island, which is most prominently illustrated in a radiogenic isotope dichotomy of the so-called ‘Kea’ and ‘Loa’ trends. Here, we present stable Fe isotopes (δ57Fe), complemented by radiogenic 176Hf/177Hf, in lavas from the Puʻu ʻŌʻō vent on the East Rift Zone of Kīlauea (Kea trend) and from the Southwest Rift Zone from Mauna Loa (Loa trend). In Hf isotopic space, lavas range from εHf = +8 to +13, with Puʻu ʻŌʻō samples aligning with published values of lavas from the ‘Kea-trend’ in SrHf isotope space, whereas Mauna Loa basalts and picrites presented here coincide with values of the Loa-trend. The Puʻu ʻŌʻō lavas have average δ57Fe of +0.15 ‰ and an average of MgO ~ 7 wt%. Mauna Loa basalts and picrites show on average lighter δ57Fe of +0.06 ‰ but higher MgO ~ 16 wt%. Combined, in a global OIB comparison, the Hawaiian samples have among the lowest δ57Fe reported to date. Both Puʻu ʻŌʻō and Mauna Loa lavas lie on a trend of decreasing δ57Fe with increasing MgO that is ascribed to olivine fractionation and accumulation. Geochemical indices of partial melting (e.g., (Gd/Yb)PM, TiO2) do not align with Fe isotopic signatures, whether corrected for crystal fractionation or not, and can thus not explain the heavy Fe isotopic signatures observed here. Likewise, a lack of co-variations of radiogenic isotopes with Fe isotopes do not support source inheritance from enriched, crustal-like plume components as the dominant process to explain the Fe isotopic signatures. On a three‑iron isotope plot, it appears that kinetic in combination or in exchange with equilibrium fractionation of stable isotopes can account for Fe isotopes lighter than δ57Fe < +0.15 ‰. An inherently heavy Fe isotope signature (> 0.15 ‰) is required to explain some of the high values measured here for which no geochemical proxy can be identified. Assuming fast recharge into the plumbing systems feeding Puʻu ʻŌʻō and Mauna Loa lavas with no presumed homogenisation in subsurface magma chambers, we surmise that the Hawaiian samples are reflecting lavas that experienced a combination of kinetic and equilibrium processes, which ultimately lead to a seemingly random distribution of Fe isotopes. The exposure to different magmatic processes, e.g., re-charging events, and the missing covariation between MgO vs δ57Fe argue for an open magma chamber underneath Puʻu ʻŌʻō and Mauna Loa, masking processes that influence Fe isotopes, i.e., source inheritance or melting degree. From our new stable Fe isotope data, we show that lavas, which travel through a plumbing system, can undergo multiple fractionation processes.

Diffusion Timescales of Magmatic Processes in the Moinui Lava Eruption at Mauna Loa, Hawai`i, as Inferred from Bimodal Olivine Populations

Abstract The 2·1 ka Moinui lava flow field, erupted from the southwest rift zone of Mauna Loa, Hawai`i, exhibits striking textural and geochemical variations, that can be used to interpret magma processes pre-, syn- and post-eruption. From this lava flow, the duration of magma storage and storage conditions, the timescales over which magma is transported to the surface, and flow emplacement mechanisms at Mauna Loa are determined. Electron microprobe analysis (EMPA) and diffusion chronometry of olivine crystals identify two distinct crystal populations: a primitive, polyhedral olivine population with core compositions of Fo90–88 and a more evolved, platy olivine population with core compositions of Fo83–82. Fe–Mg diffusion modelling of these olivine populations gives distinct timescales for each population; platy olivines yield timescales of days up to a few weeks, while polyhedral olivines yield timescales of months to years. Despite the nature of a well-insulated pāhoehoe flow, meaning that post-emplacement diffusion continues for some time, a wealth of time information can be retrieved concerning pre-eruptive magmatic processes as well as the processes associated with the lava extrusion. The short timescales obtained from the platy olivine crystals and the observed equilibrium between its cores and ambient melt suggest late-stage nucleation and crystal growth in the shallow conduit and during lava emplacement. Conversely, the longer timescales and olivine-melt disequilibrium of the polyhedral olivine crystals suggests accumulation from a deeper source and subsequent transportation to shallow magma storage beneath the summit of Mauna Loa months, or even years before eruption. The chemical and textural details of the Moinui lava reflect the mode of flow emplacement and may have implications for the interpretation of the distribution of spinifex and cumulate olivine within komatiites; high-temperature, low-viscosity lavas, common in the Archean.

Effects of habitat loss and human influence on the survival of black-crowned cranes in Africa

This paper reviews the effects of habitat loss and human influence on the survival of Black-crowned cranes in Africa. Black-crowned crane is distributed through Sahel and Sudan Savannah regions of Africa, ranging from the Senegal basin and Guinea-Bissau drainage in West Africa to the western Ethiopian Highlands and south-west Rift Valley in East Africa. Habitat loss, transformation and degradation, and illegal capture and trade for pet industry, are the major threats of the species. Because of habitat loss, the population of Black-crowned cranes is declining across its home range and it is even disappearing in some countries in West Africa. It is predicted that the population decline will continue in the future due to habitat loss and trapping of cranes for domestication. The species is classified as Lower Risk from 1988 to 2000, Near Threatened from 2004 to 2008, and then Vulnerable from 2010 onwards. The rapid conversion of wetlands and intensive land use and different economic development with poor environmental protection threat the survival of Black-crowned cranes. In order to overcome the threats of Black-crowned cranes ( Balearica pavonina ), habitat conservation and community awareness are essential, not optional.

A retrospective look at the February 1993 east rift zone intrusion at Kīlauea volcano, Hawaii

The February 1993 dike intrusion in the East Rift Zone (ERZ) of Kīlauea Volcano, Hawai'i, was recognized from tilt and seismic data, but ground-based geodetic data were too sparse to constrain the characteristics of the intrusion. Analysis of Interferometric Synthetic Aperture Radar (InSAR) from the Japan Aerospace Exploration Agency (JAXA) JERS-1 satellite reveals a maximum of ~30 cm of line-of-sight (LOS) displacement occurring near Makaopuhi Crater in the middle ERZ of Kīlauea. We model this deformation signal as a subvertical dike using a 3D-Mixed Boundary Element Method (3D-MBEM) paired with a nonlinear inversion algorithm to find the best-fit model. The best-fit dike is located just to the west of Makaopuhi Crater striking N50°W, extends to within 100 m of the surface, is ~1.3 km in length by ~4.2 km in width along strike, and has a total volume of ~7.4 × 106 m3. In addition, a post-intrusion interferogram from JERS-1 spanning 1993–1997 was analyzed. Guided by previous results, our model for the 4-year period consists of opening of the deep rift zones by about 0.5 m at 3–8.5 km depth beneath the Southwest Rift Zone (SWRZ), ERZ and the summit. A sub-horizontal detachment fault is connected to the seaward side of the vertical dike-like source to mimic the décollement known to exist beneath the volcano. We classify the 1993 dike intrusion as a passive intrusion similar to those that occurred in 1997 and 1999. Passive intrusions lack precursory inflation at Kīlauea's summit, and the likely triggering mechanism is persistent deep rift opening combined with seaward motion of the south flank along the basal décollement. Passive intrusions make forecasting and hazard assessment difficult since they are not preceded by inflation nor by large increases in seismicity.

Focused seismicity triggered by flank instability on Kīlauea's Southwest Rift Zone

Swarms of earthquakes at the head of the Southwest Rift Zone on Kīlauea Volcano, Hawaiʻi, reveal an interaction of normal and strike-slip faulting associated with movement of Kīlauea's south flank. A relocated subset of earthquakes between January 2012 and August 2014 are highly focused in space and time at depths that are coincident with the south caldera magma reservoir beneath the southern margin of Kīlauea Caldera. Newly calculated focal mechanisms are dominantly dextral shear with a north-south preferred fault orientation. Two earthquakes within this focused area of seismicity have normal faulting mechanisms, indicating two mechanisms of failure in very close proximity (10's of meters to 100m). We suggest a model where opening along the Southwest Rift Zone caused by seaward motion of the south flank permits injection of magma and subsequent freezing of a plug, which then fails in a right-lateral strike-slip sense, consistent with the direction of movement of the south flank. The seismicity is concentrated in an area where a constriction occurs between a normal fault and the deeper magma transport system into the Southwest Rift Zone. Although in many ways the Southwest Rift Zone appears analogous to the more active East Rift Zone, the localization of the largest seismicity (>M2.5) within the swarms to a small volume necessitates a different model than has been proposed to explain the lineament outlined by earthquakes along the East Rift Zone.

Plateaus and sinuous ridges as the fingerprints of lava flow inflation in the Eastern Tharsis Plains of Mars

The Tharsis Montes rift aprons are composed of outpourings of lava from chaotic terrains to the northeast and southwest flank of each volcano. Sinuous and branching channel networks that are present on the rift aprons suggest the possibility of fluvial processes in their development, or erosion by rapidly emplaced lavas, but the style of lava flow emplacement throughout rift apron development is not clearly understood. To better characterize the style of lava emplacement and role of fluvial processes in rift apron development, we conducted morphological mapping of the Pavonis Mons southwest rift apron and the eastern Tharsis plains using images from the High Resolution Imaging Science Experiment (HiRISE), Mars Orbiter Camera (MOC), Context Camera (CTX), Thermal Emission Imaging System (THEMIS), and High Resolution Stereo Camera (HRSC) along with the Mars Orbiter Laser Altimeter (MOLA) Precision Experiment Data Records (PEDRs) and gridded data. Our approach was to: (1) search for depositional fans at the slope break between the rift apron and adjacent low slope plains; (2) determine if there is evidence that previously formed deposits might have been buried by plains units; (3) characterize the Tharsis plains morphologies east of Pavonis Mons; and (4) assess their relationship to the rift apron units. We have not identified topographically significant depositional fans, nor did we observe evidence to suggest that plains units have buried older rift apron units. Flow features associated with the rift apron are observed to continue across the slope break onto the plains. In this area, the plains are composed of a variety of small fissures and low shield vents around which broad channel-fed and tube-fed flows have been identified. We also find broad, flat-topped plateaus and sinuous ridges mixed among the channels, tubes and vents. Flat-topped plateaus and sinuous ridges are morphologies that are analogous to those observed on the coastal plain of Hawai‘i, where lava flows have advanced from the volcano's several degree flank onto the nearly zero degree coastal plain. When local volumetric flow rates are low, flow fronts tend to spread laterally and often thicken via endogenous growth, or inflation, of the sheet-like flow units. If flow advance is restricted by existing topography into narrow pathways, inflation can be focused into sinuous, elongate ridges. The presence of plateaus and ridges—emplaced from the rift zones, across the plains to the east of Pavonis Mons—and a lack of fan-like features, or evidence for their burial, are consistent with rift apron lavas crossing a slope break with low local volumetric flow rates that led to inflation of sheet-like and tube-fed lava flows.

Joint analysis of geodetic and earthquake fault-plane solution data to constrain magmatic sources: A case study from Kīlauea Volcano

A joint analysis of geodetic and seismic datasets from Kīlauea Volcano during a period of magmatic unrest in 2006 demonstrates the effectiveness of this combination for testing and constraining models of magma dynamics for a complex, multi-source system. At the end of 2003, Kīlauea's summit began a four-year-long period of inflation due to a surge in magma supply to the volcano. In 2006, for the first time since 1982, Kīlauea's Southwest Rift Zone (SWRZ) also experienced inflation. To investigate the characteristics of active magma sources and the nature of their interactions with faults in the SWRZ during 2006, we integrate, through Coulomb stress modeling, contemporary geodetic data from InSAR and GPS with a new catalogue of double-couple fault-plane solutions for volcano-tectonic earthquakes. We define two periods of inflation during 2006 based on the rate of deformation measured in daily GPS data, spanning February to 15 March 2006 (Period 1) and 16 March to 30 September 2006 (Period 2). InSAR data for these two periods are inverted to determine the position, change in size, and shape of inflation sources in each period. Our new models are consistent with microseismic activity from each period. They suggest that, during Period 1, deformation in the SWRZ can be explained by pressurization of magma in a spherical reservoir beneath the south caldera, and that, during Period 2, magma was also aseismically intruded farther to the southwest into the SWRZ along a sub-horizontal plane. Our Coulomb stress analysis shows that the microseismicity recorded in the SWRZ is induced by overpressurization of the south caldera reservoir, and not by magma intrusion into the SWRZ. This study highlights the importance of a joint analysis of independent geophysical datasets to fully constrain the nature of magma accumulation.

Spatiotemporal model of Kīlauea's summit magmatic system inferred from InSAR time series and geometry‐free time‐dependent source inversion

AbstractKīlauea volcano, Hawai‘i Island, has a complex magmatic system including summit reservoirs and rift zones. Kinematic models of the summit reservoir have so far been limited to first‐order analytical solutions with predetermined geometry. To explore the complex geometry and kinematics of the summit reservoir, we apply a multitrack wavelet‐based InSAR (interferometric synthetic aperture radar) algorithm and a novel geometry‐free time‐dependent modeling scheme. To map spatiotemporally distributed surface deformation signals over Kīlauea's summit, we process synthetic aperture radar data sets from two overlapping tracks of the Envisat satellite, including 100 images during the period 2003–2010. Following validation against Global Positioning System data, we invert the surface deformation time series to constrain the spatiotemporal evolution of the magmatic system without any prior knowledge of the source geometry. The optimum model is characterized by a spheroidal and a tube‐like zone of volume change beneath the summit and the southwest rift zone at 2–3 km depth, respectively. To reduce the model dimension, we apply a principal component analysis scheme, which allows for the identification of independent reservoirs. The first three PCs, explaining 99% (63.8%, 28.5%, and 6.6%, respectively) of the model, include six independent reservoirs with a complex interaction suggested by temporal analysis. The data and model presented here, in agreement with earlier studies, improve the understanding of Kīlauea's plumbing system through enhancing the knowledge of temporally variable magma supply, storage, and transport beneath the summit, and verify the link between summit magmatic activity, seismicity, and rift intrusions.

Making it and breaking it in the Midwest: Continental assembly and rifting from modeling of EarthScope magnetotelluric data

A three-dimensional lithospheric-scale resistivity model of the North American mid-continent has been estimated based upon EarthScope magnetotelluric data. Details of the resistivity model are discussed in relation to lithospheric sutures, defined primarily from aeromagnetic and geochronologic data, which record the southward growth of the Laurentian margin in the Proterozoic. The resistivity signature of the 1.1Ga Mid-continent Rift System is examined in detail, in particular as relates to rift geometry, extent, and segmentation. An unrecognized expanse of (concealed) Proterozoic deltaic deposits in Kansas is identified and speculated to result from axial drainage along the southwest rift arm akin to the Rio Grande delta which drains multiple rift basins. A prominent conductor traces out Cambrian rifting in Arkansas, Missouri, Tennessee, and Kentucky; this linear conductor has not been imaged before and suggests that the Cambrian rift system may have been more extensive than previously thought. The highest conductivity within the mid-continent is imaged in Minnesota, Michigan, and Wisconsin where it is coincident with Paleoproterozoic metasedimentary rocks. The high conductivity is attributed to metallic sulfides, and in some cases, graphite. The former is a potential source of sulfur for multiple mineral deposits types, occurrences of which are found throughout the region. Finally, the imprint left within the mantle following the 1.1Ga rifting event is examined. Variations in lithospheric mantle conductivity are observed and are interpreted to reflect variations in water content (depleted versus metasomatized mantle) imprinted upon the mantle by the Keweenawan mantle plume.

Local earthquake tomography with the inclusion of full topography and its application to Kīlauea volcano, Hawai'i

We develop a new three-dimensional local earthquake tomography algorithm with the inclusion of full topography (LETFT). We present both synthetic and real data tests based on the P- and S-wave arrival time data for Kīlauea volcano in Hawai'i. A total of 33,768 events with 515,711 P-picks and 272,217 S-picks recorded by 35 stations at the Hawaiian Volcano Observatory are used in these tests. The comparison between the new and traditional methods based on the synthetic test shows that our new algorithm significantly improves the accuracy of the velocity model, especially at shallow depths. In the real data application, the P- and S-wave velocity models of Kīlauea show several intriguing features. We observe discontinuous high Vp (>7.0km/s) and Vs (>3.9km/s) zones at 5–14km depth below Kīlauea caldera, its East Rift Zone (ERZ) and the Southwest Rift Zone, which may represent consolidated intrusive gabbro-ultramafic cumulates. At Kīlauea caldera, Vp and Vs decrease from ~3.9km/s and ~2.6km/s from the surface to ~3.7km/s and ~2.3km/s at 2km depth. We resolve a high Vp zone (>7.0km/s) at 5–14km depth and high Vs zone (>3.9km/s) at 5–11km depth. This high Vp and Vs zone extends to the north of the ERZ at 5–10km depth and to the upper ERZ at 8–12km depth. In the Hilina Fault System, there is a high Vp layer (~7.0km/s) at 4–6km depth and a low Vp body of ~5.7km/s at 6–11km depth. The high Vp layer could be associated with the intrusive ultramafic gabbro sills. The velocity contrast on the north and south sides of the Koa'e Fault System indicates that the intrusive activities mainly occur to the north of the fault. Our new LETFT method performs well in both the synthetic and real data tests and we expect that it will reveal more robust velocity structures in areas with larger topographic variations.

231Pa systematics in postglacial volcanic rocks from Iceland

Several recent studies have highlighted the potential of combined 238U–230Th and 235U–231Pa systematics to constrain upwelling rates and the role of recycled mafic lithologies in mantle plume-derived basalts. Accordingly, we present measurements of the 231Pa concentrations from 26 mafic volcanic rocks from Iceland, including off-axis basalts from the Snaefellsnes Peninsula, to complement previously published 238U–230Th–226Ra data. 231Pa concentrations vary from 27 to 624fg/g and (231Pa/235U) ratios from 1.12 to 2.11 with the exception of one anomalous sample from the Southeast Rift which has a 231Pa deficit with (231Pa/235U)=0.86. An important new result is that basalts from the Southeast Rift and the Snaefellsnes Peninsula define a trend at relatively low (231Pa/235U) for a given (230Th/238U) ratio. Many of the remaining samples fall in or around the global field for ocean island basalts but those from the Mid-Iceland Belt and the Southwest Rift/Reykjanes Peninsula extend to higher (231Pa/235U) ratios at a given (230Th/238U), similar to mid-ocean ridge basalts. In principle, these lavas could result from melting of peridotite at lower pressures. However, there is no reason to suspect that the Mid-Iceland Belt and the Southwest Rift lavas reflect shallower melting than elsewhere in Iceland. In our preferred model, these lavas reflect melting of garnet peridotite whereas those from the Southeast Rift and the Snaefellsnes Peninsula contain a significant contribution (up to 20%) of melt from garnet pyroxenite. This is consistent with incompatible trace element and radiogenic isotope evidence for recycled oceanic crust in these lavas. There is increasing agreement that the displacement of ocean island basalts to lower (231Pa/235U) ratios at a given (230Th/238U), compared to mid-ocean ridge basalts, reflects the role of recycled mafic lithologies such as garnet pyroxenite as well as higher average pressures of melting. It now seems likely that this interpretation may have more general application to plume-affected MORB that exhibit the same U-series disequilibria characteristics.

The origin of Mauna Loa's Nīnole Hills: Evidence of rift zone reorganization

In order to identify the origin of Mauna Loa volcano's Nīnole Hills, Bouguer gravity was used to delineate density contrasts within the edifice. Our survey identified two residual anomalies beneath the Southwest Rift Zone (SWRZ) and the Nīnole Hills. The Nīnole Hills anomaly is elongated, striking northeast, and in inversions both anomalies merge at approximately −7 km above sea level. The positive anomaly, modeled as a rock volume of ~1200 km3 beneath the Nīnole Hills, is associated with old eruptive vents. Based on the geologic and geophysical data, we propose that the gravity anomaly under the Nīnole Hills records an early SWRZ orientation, now abandoned due to geologically rapid rift-zone reorganization. Catastrophic submarine landslides from Mauna Loa's western flank are the most likely cause for the concurrent abandonment of the Nīnole Hills section of the SWRZ. Rift zone reorganization induced by mass wasting is likely more common than currently recognized.

Seismic tomography of compressional wave attenuation structure for Kı̄lauea Volcano, Hawai‘i

AbstractWe present a frequency‐independent three‐dimensional (3‐D) compressional wave attenuation model (indicated by the reciprocal of quality factor Qp) for Kı̄lauea Volcano in Hawai‘i. We apply the simul2000 tomographic algorithm to the attenuation operator t* values for the inversion of Qp perturbations through a recent 3‐D seismic velocity model and earthquake location catalog. The t* values are measured from amplitude spectra of 26708 P wave arrivals of 1036 events recorded by 61 seismic stations at the Hawaiian Volcanology Observatory. The 3‐D Qp model has a uniform horizontal grid spacing of 3 km, and the vertical node intervals range between 2 and 10 km down to 35 km depth. In general, the resolved Qp values increase with depth, and there is a correlation between seismic activity and low‐Qp values. The area beneath the summit caldera is dominated by low‐Qp anomalies throughout the entire resolved depth range. The Southwest Rift Zone and the East Rift Zone exhibit very high Qp values at about 9 km depth, whereas the shallow depths are characterized with low‐Qp anomalies comparable with those in the summit area. The seismic zones and fault systems generally display relatively high Qp values relative to the summit. The newly developed Qp model provides an important complement to the existing velocity models for exploring the magmatic system and evaluating and interpreting intrinsic physical properties of the rocks in the study area.