High‐resolution magnetic data collected along the axis and the south flank of the Puna Ridge, the submarine extension of the East Rift Zone of Kilauea Volcano, Hawaii, aid in the interpretation of magmatic processes. Crustal magnetization varies along the crest of the Puna Ridge. Three prominent magnetization highs are located between the depths of 1500 and 2500 m. Each of the magnetization highs is associated with distinct volcanic features: a flow field on top of fissured terrain, a lava plateau, and a staircase of lava terraces. The crustal magnetization highs can be explained by the presence of these volcanic features, taking into account water depth and the size of the volcanic features and making simple assumptions about the rock magnetization. No measurements of rock magnetization have been made on Puna Ridge samples. We assume two values: 50 and 10 A/m, based on values from mid‐ocean ridge basalts and the estimated ages of Puna Ridge rocks [ Clague et al. , 1995]. The two shallowest magnetization highs are adjacent to a landslide scarp. Lava flows overtop the scarp suggesting that a dike has been emplaced and eruptions have occurred since the landslide formed. A possible scenario is that the collapse of the flank led to dike intrusion and eruption at this location because of a reduction in the horizontal stress within the edifice. Our observations suggest that flows are starting to fill in the gap caused by the landslide to rebuild this section of the north flank.
This dataset contains the results of various Ensemble Kalman Filter (EnKF) inversions in which synthetic GNSS and InSAR observations from an inflating magma system are assimilated into numerical models of rock deformation around a pressurized ellipsoidal magma reservoir. Each inversion uses a different variant of the EnKF, with changes to workflow meta-parameters such as the number of ensemble members or the particular update algorithm used. In particular, each filter variant is evaluated by comparing the final output model to the original synthetic model. The specific performance criteria used include (1) the root mean square error (RMSE) between the model predictions and the assimilated observations, as well as normalized misfit terms measuring the filter's ability to resolve (2) reservoir wall tensile stress, (3) easily-observable unique parameters such as reservoir position and aspect ratio, and (4) difficult-to-derive non-unique parameters such as the specific size and internal pressure of the reservoir. The assimilated data include two different scenarios, one in which inflation is caused by pressurization and another in which it is driven by a lateral reservoir expansion. Both datasets are tested with each EnKF variant. Finally, we include example matrices from within an EnKF update step to demonstrate inter-parameter correlations that develop during the assimilation and how they can be mitigated through randomization.
Abstract We utilize 3‐D temperature‐dependent viscoelastic finite element models to investigate the mechanical response of the host rock supporting large caldera‐size magma reservoirs (volumes >10 2 km 3 ) to local tectonic stresses. The mechanical stability of the host rock is used to determine the maximum predicted repose intervals and magma flux rates that systems may experience before successive eruption is triggered. Numerical results indicate that regional extension decreases the stability of the roof rock overlying a magma reservoir, thereby promoting early‐onset caldera collapse. Alternatively, moderate amounts of compression (≤10 mm/year) on relatively short timescales (<10 4 years) increases roof rock stability. In addition to quantifying the affect of tectonic stresses on reservoir stability, our models indicate that the process of rejuvenation and mechanical failure is likely to take place over short time periods of hundreds to thousands of years. These findings support the short preeruption melt accumulation timescales indicated by U series disequilibrium studies.
Data for the paper will be published in JGR: solid earth Key points:A rotating pressure source determined by InSAR data indicates how crystal mush is pressurizedDilatancy in a pre-existing weak zone may open channels for pore fluid to inject into and to trigger volcano-tectonic earthquakesA water table drop due to dilatation in a pre-existing fault zone may have drained the crater lake at Korovin and triggered the phreatic eruption
New data reveal details of the post-caldera history at the Earth's youngest resurgent supervolcano, Toba caldera in Sumatra. Resurgence after the caldera-forming ~74 ka Youngest Toba Tuff eruption uplifted the caldera floor as a resurgent dome, Samosir Island, capped with 100m of lake sediments. 14C age data from the uppermost datable sediments reveal that Samosir Island was submerged beneath lake level (~900m a.s.l) ~33.7 ky. Since then, Samosir experienced 700m of uplift as a tilted block dipping to the west. Using 14C ages and elevations of sediment along a transect of Samosir reveal that minimum uplift rates were ~4.9 cm/yr from ~33.7 to 22.5 ka, but diminished to ~0.7 cm/yr after 22.5ka. Thermo-mechanical models informed by these rates reveal that detumescence does not produce the uplift nor the uplift rates estimated for Samosir. However, models calculating the effect of volume change of the magma reservoir within a temperature-dependent viscoelastic host rock reveal that a single pulse of ~475 km3 of magma produces a better fit to the uplift data than a constant flux. Reproducing the uplift rates require more sophisticated models. Motivation for resurgent uplift of the caldera floor is rebound of remnant magma as the system re-established magmastatic and isostatic equilibrium after the caldera collapse. Previous assertions that the caldera floor was apparently at 400m a.s.l or lower requires that uplift must have initiated between sometime between 33.7 ka and 74 ka at a minimum average uplift rate of ~1.1 cm/ year. The change in uplift rate from pre-33.7 ka to immediately post-33.7 ka suggests a role for deep recharge augmenting rebound. Average minimum rates of resurgent uplift at Toba are at least an order of magnitude slower than net rates of "restlessness" at currently active calderas. This connotes a distinction between resurgence and "restlessness" controlled by different processes, scales of process, and controlling variables.
Abstract Understanding how shallow reservoirs store and redirect magma is critical for deciphering the relationship between surface and subsurface volcanic activity on the terrestrial planets. Complementing field, laboratory and remote sensing analyses, elastic models provide key insights into the mechanics of magma reservoir inflation and rupture, and hence into commonly observed volcanic phenomena including edifice growth, circumferential intrusion, radial dyke swarm emplacement and caldera formation. Based on finite element model results, the interplay between volcanic elements – such as magma reservoir geometry, host rock environment (with an emphasis on understanding how host rock pore pressure assumptions affect model predictions), mechanical layering, and edifice loading with and without flexure – dictates the overpressure required for rupture, the location and orientation of initial fracturing and intrusion, and the associated surface uplift. Model results are either insensitive to, or can readily incorporate, material and parameter variations characterizing different planetary environments, and they also compare favourably with predictions derived from rheologically complex, time-dependent formulations for a surprisingly diverse array of volcanic scenarios. These characteristics indicate that elastic models are a powerful and useful tool for exploring many fundamental questions in planetary volcanology.
