Search

Trench‐parallel translation of the Central American Forearc (CAFA) is the result of strain partitioning along the Cocos and Caribbean (CA) convergent margin. Unlike the tectonics of northwestern Costa Rica and El Salvador, CAFA‐CA relative motion in Nicaragua is not accommodated on margin‐parallel fault systems. Rather, the northwest‐trending dextral shear is accommodated on margin‐normal sinistral strike‐slip faults that approximate the motion of a margin‐parallel fault (i.e., bookshelf faulting). We compare a new Global Positioning System interseismic horizontal velocity field to analytical and numerical models to show that the bookshelf faulting model can produce the observed velocity field and provide insight into the kinematics and configuration of the margin‐normal fault system. We find that a fault system with 20 km‐long parallel to sub‐parallel margin‐normal sinistral faults, spaced ∼5 km apart, locked from the surface to 5 km depth, and with interseismic slip deficits of 4 mm yr−1, can replicate the observed velocity field. These findings have implications for the region's seismic hazard where shallow moderate‐magnitude earthquakes will have reoccurrence intervals of ∼50 years. These findings are also important for volcanic hazard estimation and unrest forecasting because the margin‐normal faults are in the volcanic arc and magma‐tectonic interactions have been documented along the CAFA. Plain Language Summary: The Central American Forearc (CAFA) is a northwest‐moving tectonic block that is located between the Middle America Trench and the Central American Volcanic Arc (CAVA). Northwest‐oriented faults (i.e., faults parallel to the CAVA motion), have been found in Costa Rica, El Salvador, and Guatemala, representing the boundary between CAFA and CAVA. In eastern El Salvador, Nicaragua, and western Costa Rica, a different fault style has been identified where the faults are oriented perpendicular to the motion of the CAVA (bookshelf faulting). These faults are ancient structures that take up the motion between the CAFA and CAVA, which is easier than creating new faults parallel to the motion. We show that, in northwestern Nicaragua, the bookshelf faulting can produce the observed GPS data. We find that these faults should move at 4 mm/yr, are close to populated areas, and will produce M5 to M6 earthquakes every 50 years. Key Points: Bookshelf faulting accommodates northwest‐directed forearc motion on margin‐normal sinistral faultsStrain accumulation on margin‐normal sinistral faults agrees with GPS‐derived horizontal velocity field in northwestern NicaraguaShallow M > 5 upper‐plate earthquakes in the region will have reoccurrence intervals of ∼50 years [ABSTRACT FROM AUTHOR]

Many examples of exposed giant dike swarms can be found where lateral magma flow has exceeded hundreds of kilometers. We show that massive magma flow into dikes can be established with only modest overpressure in a magma body if a large enough pathway opens at its boundary and gradual buildup of high tensile stress has occurred along the dike pathway prior to the onset of diking. This explains rapid initial magma flow rates, modeled up to about 7400 cubic meters per second into a dike ~15-kilometers long, which propagated under the town of Grindavík, Southwest Iceland, in November 2023. Such high flow rates provide insight into the formation of major dikes and imply a serious hazard potential for high-flow rate intrusions that propagate to the surface and transition into eruptions.

Increased rates of deformation and seismicity are well-established precursors to volcanic eruptions, and their interpretation forms the basis for eruption warnings worldwide. Rates of ground displacement and the number of earthquakes escalate before many eruptions 1-3 , as magma forces its way towards the surface. However, the pre-eruptive patterns of deformation and seismicity vary widely. Here we show how an eruption beginning on 19 March 2021 at Fagradalsfjall, Iceland, was preceded by a period of tectonic stress release ending with a decline in deformation and seismicity over several days preceding the eruption onset. High rates of deformation and seismicity occurred from 24 February to mid-March in relation to gradual emplacement of an approximately 9-km-long magma-filled dyke, between the surface and 8 km depth (volume approximately 34 × 10 6  m 3 ), as well as the triggering of strike-slip earthquakes up to magnitude M W 5.64. As stored tectonic stress was systematically released, there was less lateral migration of magma and a reduction in both the deformation rates and seismicity. Weaker crust near the surface may also have contributed to reduced seismicity, as the depth of active magma emplacement progressively shallowed. This demonstrates that the interaction between volcanoes and tectonic stress as well as crustal layering need to be fully considered when forecasting eruptions., (© 2022. The Author(s).)

Large volume effusive eruptions with relatively minor observed precursory signals are at odds with widely used models to interpret volcano deformation. Here we propose a new modelling framework that resolves this discrepancy by accounting for magma buoyancy, viscoelastic crustal properties, and sustained magma channels. At low magma accumulation rates, the stability of deep magma bodies is governed by the magma-host rock density contrast and the magma body thickness. During eruptions, inelastic processes including magma mush erosion and thermal effects, can form a sustained channel that supports magma flow, driven by the pressure difference between the magma body and surface vents. At failure onset, it may be difficult to forecast the final eruption volume; pressure in a magma body may drop well below the lithostatic load, create under-pressure and initiate a caldera collapse, despite only modest precursors.

Crust at many divergent plate boundaries forms primarily by the injection of vertical sheet-like dykes, some tens of kilometres long. Previous models of rifting events indicate either lateral dyke growth away from a feeding source, with propagation rates decreasing as the dyke lengthens, or magma flowing vertically into dykes from an underlying source, with the role of topography on the evolution of lateral dykes not clear. Here we show how a recent segmented dyke intrusion in the Bárðarbunga volcanic system grew laterally for more than 45 kilometres at a variable rate, with topography influencing the direction of propagation. Barriers at the ends of each segment were overcome by the build-up of pressure in the dyke end; then a new segment formed and dyke lengthening temporarily peaked. The dyke evolution, which occurred primarily over 14 days, was revealed by propagating seismicity, ground deformation mapped by Global Positioning System (GPS), interferometric analysis of satellite radar images (InSAR), and graben formation. The strike of the dyke segments varies from an initially radial direction away from the Bárðarbunga caldera, towards alignment with that expected from regional stress at the distal end. A model minimizing the combined strain and gravitational potential energy explains the propagation path. Dyke opening and seismicity focused at the most distal segment at any given time, and were simultaneous with magma source deflation and slow collapse at the Bárðarbunga caldera, accompanied by a series of magnitude M > 5 earthquakes. Dyke growth was slowed down by an effusive fissure eruption near the end of the dyke. Lateral dyke growth with segment barrier breaking by pressure build-up in the dyke distal end explains how focused upwelling of magma under central volcanoes is effectively redistributed over long distances to create new upper crust at divergent plate boundaries.

Gradual inflation of magma chambers often precedes eruptions at highly active volcanoes. During such eruptions, rapid deflation occurs as magma flows out and pressure is reduced. Less is known about the deformation style at moderately active volcanoes, such as Eyjafjallajökull, Iceland, where an explosive summit eruption of trachyandesite beginning on 14 April 2010 caused exceptional disruption to air traffic, closing airspace over much of Europe for days. This eruption was preceded by an effusive flank eruption of basalt from 20 March to 12 April 2010. The 2010 eruptions are the culmination of 18 years of intermittent volcanic unrest. Here we show that deformation associated with the eruptions was unusual because it did not relate to pressure changes within a single magma chamber. Deformation was rapid before the first eruption (>5 mm per day after 4 March), but negligible during it. Lack of distinct co-eruptive deflation indicates that the net volume of magma drained from shallow depth during this eruption was small; rather, magma flowed from considerable depth. Before the eruption, a ∼0.05 km(3) magmatic intrusion grew over a period of three months, in a temporally and spatially complex manner, as revealed by GPS (Global Positioning System) geodetic measurements and interferometric analysis of satellite radar images. The second eruption occurred within the ice-capped caldera of the volcano, with explosivity amplified by magma-ice interaction. Gradual contraction of a source, distinct from the pre-eruptive inflation sources, is evident from geodetic data. Eyjafjallajökull's behaviour can be attributed to its off-rift setting with a 'cold' subsurface structure and limited magma at shallow depth, as may be typical for moderately active volcanoes. Clear signs of volcanic unrest signals over years to weeks may indicate reawakening of such volcanoes, whereas immediate short-term eruption precursors may be subtle and difficult to detect.

Pressure influences both magma production and the failure of magma chambers. Changes in pressure interact with the local tectonic settings and can affect magmatic activity. Present-day reduction in ice load on subglacial volcanoes due to global warming is modifying pressure conditions in magmatic systems. The large pulse in volcanic production at the end of the last glaciation in Iceland suggests a link between unloading and volcanism, and models of that process can help to evaluate future scenarios. A viscoelastic model of glacio-isostatic adjustment that considers melt generation demonstrates how surface unloading may lead to a pulse in magmatic activity. Iceland's ice caps have been thinning since 1890 and glacial rebound at rates exceeding 20 mm yr(-1) is ongoing. Modelling predicts a significant amount of 'additional' magma generation under Iceland due to ice retreat. The unloading also influences stress conditions in shallow magma chambers, modifying their failure conditions in a manner that depends critically on ice retreat, the shape and depth of magma chambers as well as the compressibility of the magma. An annual cycle of land elevation in Iceland, due to seasonal variation of ice mass, indicates an annual modulation of failure conditions in subglacial magma chambers.