Sebastian Heimann

Gradual caldera collapse at Bárdarbunga volcano, Iceland, regulated by lateral magma outflow

Driven to collapse Volcanic eruptions occur frequently, but only rarely are they large enough to cause the top of the mountain to collapse and form a caldera. Gudmundsson et al. used a variety of geophysical tools to monitor the caldera formation that accompanied the 2014 Bárdarbunga volcanic eruption in Iceland. The volcanic edifice became unstable as magma from beneath Bárdarbunga spilled out into the nearby Holuhraun lava field. The timing of the gradual collapse revealed that it is the eruption that drives caldera formation and not the other way around. Science, this issue p. 262 Magma flow from under the Bárdarbunga volcano drove caldera collapse during the 2014 eruption. INTRODUCTION The Bárdarbunga caldera volcano in central Iceland collapsed from August 2014 to February 2015 during the largest eruption in Europe since 1784. An ice-filled subsidence bowl, 110 square kilometers (km2) in area and up to 65 meters (m) deep developed, while magma drained laterally for 48 km along a subterranean path and erupted as a major lava flow northeast of the volcano. Our data provide unprecedented insight into the workings of a collapsing caldera. RATIONALE Collapses of caldera volcanoes are, fortunately, not very frequent, because they are often associated with very large volcanic eruptions. On the other hand, the rarity of caldera collapses limits insight into this major geological hazard. Since the formation of Katmai caldera in 1912, during the 20th century’s largest eruption, only five caldera collapses are known to have occurred before that at Bárdarbunga. We used aircraft-based altimetry, satellite photogrammetry, radar interferometry, ground-based GPS, evolution of seismicity, radio-echo soundings of ice thickness, ice flow modeling, and geobarometry to describe and analyze the evolving subsidence geometry, its underlying cause, the amount of magma erupted, the geometry of the subsurface caldera ring faults, and the moment tensor solutions of the collapse-related earthquakes. RESULTS After initial lateral withdrawal of magma for some days though a magma-filled fracture propagating through Earth’s upper crust, preexisting ring faults under the volcano were reactivated over the period 20 to 24 August, marking the onset of collapse. On 31 August, the eruption started, and it terminated when the collapse stopped, having produced 1.5 km of basaltic lava. The subsidence of the caldera declined with time in a near-exponential manner, in phase with the lava flow rate. The volume of the subsidence bowl was about 1.8 km3. Using radio-echo soundings, we find that the subglacial bedrock surface after the collapse is down-sagged, with no indications of steep fault escarpments. Using geobarometry, we determined the depth of magma reservoir to be ~12 km, and modeling of geodetic observations gives a similar result. High-precision earthquake locations and moment tensor analysis of the remarkable magnitude M5 earthquake series are consistent with steeply dipping ring faults. Statistical analysis of seismicity reveals communication over tens of kilometers between the caldera and the dike. CONCLUSION We conclude that interaction between the pressure exerted by the subsiding reservoir roof and the physical properties of the subsurface flow path explain the gradual near-exponential decline of both the collapse rate and the intensity of the 180-day-long eruption. By combining our various data sets, we show that the onset of collapse was caused by outflow of magma from underneath the caldera when 12 to 20% of the total magma intruded and erupted had flowed from the magma reservoir. However, the continued subsidence was driven by a feedback between the pressure of the piston-like block overlying the reservoir and the 48-km-long magma outflow path. Our data provide better constraints on caldera mechanisms than previously available, demonstrating what caused the onset and how both the roof overburden and the flow path properties regulate the collapse. The Bárdarbunga caldera and the lateral magma flow path to the Holuhraun eruption site. (A) Aerial view of the ice-filled Bárdarbunga caldera on 24 October 2014, view from the north. (B) The effusive eruption in Holuhraun, about 40 km to the northeast of the caldera

Systematic Changes of Earthquake Rupture with Depth: A Case Study from the 2010Mw 8.8 Maule, Chile, Earthquake Aftershock Sequence

The very shallow part of subduction megathrusts occasionally hosts tsunami earthquakes, with unusually slow rupture propagation. The aftershock sequence of the 2010 Mw 8.8 Maule earthquake, offshore Chile, provides us with the opportunity to study systematic changes in source properties for smaller earthquakes within a single segment of a subduction zone. We invert amplitude spectra for double-couple moment tensors and centroid depths of 71 aftershocks of the Maule earthquake down to magnitudes Mw 4.0. In addition, we also derive average source durations.We find that shallower earthquakes tend to have longer normalized source durations on average, similar to the pattern observed previously for larger magnitude events. This depth dependence is observable for thrust and normal earthquakes. The normalized source durations of normal- faulting earthquakes are at the lower end of those for thrust earthquakes, probably because of the higher stress drops of intraplate earthquakes compared to interplate earthquakes. We suggest from the similarity of the depth dependence of normal and thrust events and between smaller and larger magnitude earthquakes that the depth-dependent variation of rigidity, rather than frictional conditional stability at the plate interface, is primarily responsible for the observed pattern. Tsunami earthquakes probably require both low rigidity and conditionally stable frictional conditions; the presence of longduration moderate-magnitude events is therefore a helpful but not sufficient indicator for areas at risk of tsunami earthquakes.

Challenges in Regional Moment Tensor Resolution and Interpretation

The seismic moment tensor is a representation of a seismic source, described mathematically by a symmetric tensor of second order.

Seismicity in the block mountains between Halle and Leipzig, Central Germany: centroid moment tensors, ground motion simulation, and felt intensities of two M ≈ 3 earthquakes in 2015 and 2017

On April 29, 2017 at 0:56 UTC (2:56 local time), an MW = 2.8 earthquake struck the metropolitan area between Leipzig and Halle, Germany, near the small town of Markranstädt. The earthquake was felt within 50 km from the epicenter and reached a local intensity of I0 = IV. Already in 2015 and only 15 km northwest of the epicenter, a MW = 3.2 earthquake struck the area with a similar large felt radius and I0 = IV. More than 1.1 million people live in the region, and the unusual occurrence of the two earthquakes led to public attention, because the tectonic activity is unclear and induced earthquakes have occurred in neighboring regions. Historical earthquakes south of Leipzig had estimated magnitudes up to MW ≈ 5 and coincide with NW-SE striking crustal basement faults. We use different seismological methods to analyze the two recent earthquakes and discuss them in the context of the known tectonic structures and historical seismicity. Novel stochastic full waveform simulation and inversion approaches are adapted for the application to weak, local earthquakes, to analyze mechanisms and ground motions and their relation to observed intensities. We find NW-SE striking normal faulting mechanisms for both earthquakes and centroid depths of 26 and 29 km. The earthquakes are located where faults with large vertical offsets of several hundred meters and Hercynian strike have developed since the Mesozoic. We use a stochastic full waveform simulation to explain the local peak ground velocities and calibrate the method to simulate intensities. Since the area is densely populated and has sensitive infrastructure, we simulate scenarios assuming that a 12-km long fault segment between the two recent earthquakes is ruptured and study the impact of rupture parameters on ground motions and expected damage.

An Active Seismic Zone in Intraplate West Iberia Inferred From High‐Resolution Geophysical Data

Intraplate Iberia is a region of slow lithopsheric deformation (<1 mm/yr) with significant historical earthquake activity. Recent high-quality instrumental data have shown that small-magnitude earthquakes collapse along clusters and lineaments, which however do not bear a clear relationship to geologically mapped active structures. In this article, we investigate the controls of these earthquake clusters. In particular, we study two of the identified clusters—the Arraiolos and the Évora seismic zones (ASZ and ESZ), located in the Western Ossa Morena Zone, southwest Iberia. The ASZ marks a sharp boundary between a seismically active region to its south and a more quiet region to its north. We revise historical earthquakes in order to clarify whether earthquake activity in the region is persistent. We use data from a local network to compute accurate epicenters, focal depth, focal mechanisms, and spatiotemporal clustering, thus characterizing ongoing small-scale fracturing. Finally, we analyze complementary data sets, including tomographic models, Global Navigation Satellite Systems data, magnetic anomalies, and gravity anomalies, in order to discuss the factors that control seismogenesis in the two seismic zones. Consistency between earthquake locations, focal mechanisms and Global Navigation Satellite Systems data suggests that the ASZ is an active right-lateral shear zone, which divides two blocks within the Western Ossa Morena Zone. The ESZ seems to localize microseismicity due to its granitic lithology. These results suggest that high-resolution geophysical data have the potential to reveal blocks with different seismogenic and rheological behaviors, which may be used to improve our understanding of fault systems and the assessment of earthquake hazard in slowly deforming regions. Plain Language Summary Mainland Portugal is a region of slow lithospheric deformation. This means that changes in Earth’s outmost layer—the lithosphere—occur at very low rates (<1 mm/yr). In such environments, faults producing earthquakes are not easy to identify at the Earth’s surface, both because their evidence can be gradually erased by wind and water or simply because they do not reach the surface. Recent studies have shown that small earthquakes in mainland Portugal group together delineating seismically more active regions. In this article we focus in two particular groups of earthquakes—the Arraiolos and the Évora seismic zones (ASZ and ESZ) and we investigate why they occur in these particular locations. We obtain precise maps of earthquake epicenters. When possible, we also analyze the direction of slip during the earthquake and the orientation of the fracture on which it occurred. We compare our results with other data sets, such as images of the Earth’s interior, that could give hints about the constitution of crust beneath the ASZ and the ESZ. Earthquakes epicenters show fault sections at depth in the ASZ. These faults separate two crustal blocks with distinct material properties. In the ESZ earthquakes are associated to contrasts in crustal materials.