Voluminous Storage and Rapid Magma Ascent Beneath La Palma Revealed by Seismic Tomography

Seismic tomography provides a window into magmatic plumbing systems; however, obtaining sufficient data for ‘real-time’ imaging is challenging. Until now, syn-eruptive tomography has not been successfully demonstrated. For the first time, we obtained high-resolution images of Earth's interior during an ongoing volcanic eruption. We used data from 11,349 earthquakes, most of which during La Palma eruption (19 September-13 December, 2021), to perform travel-time seismic tomography. We present high-precision earthquake relocations and 3D distributions of P and S-wave velocities highlighting the geometry of magma sources. We identified three distinct structures: (1) a shallow localised region (< 3 km) of hydrothermal alteration; (2) spatially extensive, consolidated, oceanic crust extending to ~10 km depth and; (3) a large (> 400 km3) sub-crustal magma-filled rock volume intrusion extending from ~7 to 25 km depth. Our results suggest that this large magma reservoir feeds the La Palma eruption continuously for almost three months. Prior to eruption onset, magma ascended from ~10 km depth to the surface in < 7 days. In the upper 3 km, melt migration is along the western contact between consolidated oceanic crust and altered hydrothermal material. Similar structural weaknesses along the eastern contact could potentially cause new eruptive centres in the future.

Transgressive nature and chilled margins of the Upper Zone in the western Bushveld Complex, South Africa

ABSTRACT The Upper Zone of the Bushveld Complex has long been known to have formed from a major influx of magma into the chamber that caused large-scale erosion of the chamber floor cumulates. The most dramatic manifestations of this process are two major gap areas (Northern and Southern) in the western Bushveld Complex in which the Upper Zone appears to have eroded away the underlying cumulates down to the very base of the Complex. However, due to almost complete lack of outcrops in the gap areas, no direct field observations have ever been reported to confirm the transgressive nature of the Upper Zone. Here, we present for the first time such observations from the Kameelhoek chromite mine located at the margin of the Northern Gap. In the open pit we have documented several transgressive depressions (up to 40 m in width) in the orthopyroxenite and chromitites of the Lower Critical Zone that are filled in with magnetite gabbro of the Upper Zone. The magnetite gabbro is chilled against the sidewalls of the depressions, forming glassy and fine-grained textured rocks with plagioclase laths arranged in radial clusters. Mineralogically and chemically, the magnetite gabbro correlates with cumulates from the lowermost part of the Upper Zone at its normal position in the complex. Three major points that have emerged from this study are: (1) the Critical Zone has been eroded away by magma that was parental to the Upper Zone, (2) this eroding magma was not the one that initiated formation of the Pyroxenite Marker, but rather the evolved melt that replenished the chamber at some later stage, and (3) the melt was phenocryst-free and likely derived from a deep-seated staging chamber. Our study thus supports a recent notion that even during the formation of the Upper Zone, the Bushveld chamber had still been operating as an open system that was replenished by melts from deeper magma sources.

Age and Geochemistry of Late Jurassic Mafic Volcanic Rocks in the Northwestern Erguna Block, Northeast China

The northwestern Erguna Block, where a wide range of volcanic rocks are present, provides one of the foremost locations to investigate Mesozoic Paleo-Pacific and Mongol-Okhotsk subduction. The identification and study of Late Jurassic mafic volcanic rocks in the Badaguan area of northwestern Erguna is of particular significance for the investigation of volcanic magma sources and their compositional evolution. Detailed petrological, geochemical, and zircon U-Pb dating suggests that the Late Jurassic mafic volcanic rocks formed at 157–161 Ma. Furthermore, the geochemical signatures of these mafic volcanic rocks indicate that they are calc-alkaline or transitional series with weak peraluminous characteristics. The rocks have a strong MgO, Al2O3, and total alkali content, and a SiO2 content of 53.55–63.68 wt %; they are enriched in Rb, Th, U, K, and light rare-earth elements (LREE), and depleted in high-field-strength elements (HFSE), similar to igneous rocks in subduction zones. These characteristics indicate that the Late Jurassic mafic volcanic rocks in the Badaguan area may be derived from the partial melting of the lithospheric mantle as it was metasomatized by subduction-related fluid and the possible incorporation of some subducting sediments. Subsequently, the fractional crystallization of Fe and Ti oxides occurred during magmatic evolution. Combined with the regional geological data, it is inferred that the studied mafic volcanic rocks were formed by lithospheric extension after the closure of the Mongol-Okhotsk Ocean.

Disposition of magmatic eruptions and fault distribution in northwestern Saudi Arabia using pseudo-depth slice magnetic anomaly

This study presents the disposition of magmatic eruptions with a fault distribution in northwestern Saudi Arabia, where intensive magma invades the lithosphere. Structural and magmatic features are traced at successive depths through high-resolution magnetic anomaly pseudo-depth slices. The total horizontal gradient technique is applied to pseudo-depth slice magnetic anomalies to enhance the linear trends of faults and related magmatic activity. A comprehensive cross-section constructed from the projection of gradient horizontal maxima relative to pseudo-slices allows the visualization of the vertical behavior of faults and magma sources. Three major fault systems were identified, primarily aligned in the N–S, NE–SW and NW–SE directions. They are characterized by increasing length and width with depth. The N-S fault system is a major non-planar deep system throughout the area, affected by the NW–SE and NE–SW deep discontinuities. The evolution of these discontinuities with depth successfully shows magma uprising zones represented by a circular horizontal gradient, which starts to appear at a depth of 4500 m with a vertical continuity to the surface. They are interpreted as possible locations of ascending magma chambers or vents. The disposition of these magma sources with fault distribution can show a close relationship between the fault systems and the magma eruptions. The interpreted magma vents appear where the NE-trending transform faults intersect the NW or N–S fault zone. These intersections may represent weak zones that act as vertical conduits through which magma discontinuously erupts into the overlying crust, forming major volcanic fields in the eastern Red Sea.

Isotopic evolution of prehistoric magma sources of Mt. Etna, Sicily: Insights from the Valle Del Bove

Mount Etna in NE Sicily occupies an unusual tectonic position in the convergence zone between the African and Eurasian plates, near the Quaternary subduction-related Aeolian arc and above the down-going Ionian oceanic slab. Magmatic evolution broadly involves a transition from an early tholeiitic phase (~ 500 ka) to the current alkaline phase. Most geochemical investigations have focussed on either historic (> 130-years old) or recent (< 130-years old) eruptions of Mt. Etna or on the ancient basal lavas (ca. 500 ka). In this study, we have analysed and modelled the petrogenesis of alkalic lavas from the southern wall of the Valle del Bove, which represent a time span of Mt. Etna’s prehistoric magmatic activity from ~ 85 to ~ 4 ka. They exhibit geochemical variations that distinguish them as six separate lithostratigraphic and volcanic units. Isotopic data (143Nd/144Nd = 0.51283–0.51291; 87Sr/86Sr = 0.70332–0.70363; 176Hf/177Hf = 0.28288–0.28298; 206Pb/204Pb = 19.76–20.03) indicate changes in the magma source during the ~ 80 kyr of activity that do not follow the previously observed temporal trend. The oldest analysed Valle del Bove unit (Salifizio-1) erupted basaltic trachyandesites with variations in 143Nd/144Nd and 87Sr/86Sr ratios indicating a magma source remarkably similar to that of recent Etna eruptions, while four of the five subsequent units have isotopic compositions resembling those of historic Etna magmas. All five magma batches are considered to be derived from melting of a mixture of spinel lherzolite and pyroxenite (± garnet). In contrast, the sixth unit, the main Piano Provenzana formation (~ 42–30 ka), includes the most evolved trachyandesitic lavas (58–62 wt% SiO2) and exhibits notably lower 176Hf/177Hf, 143Nd/144Nd, and 206Pb/204Pb ratios than the other prehistoric Valle del Bove units. This isotopic signature has not yet been observed in any other samples from Mt. Etna and we suggest that the parental melts of the trachyandesites were derived predominantly from ancient pyroxenite in the mantle source of Etna.

Combining High- and Low-Rate Geodetic Data Analysis for Unveiling Rapid Magma Transfer Feeding a Sequence of Violent Summit Paroxysms at Etna in Late 2015

We propose a multi-temporal-scale analysis of ground deformation data using both high-rate tilt and GNSS measurements and the DInSAR and daily GNSS solutions in order to investigate a sequence of four paroxysmal episodes of the Voragine crater occurring in December 2015 at Mt. Etna (Italy). The analysis aimed at inferring the magma sources feeding a sequence of very violent eruptions, in order to understand the dynamics and to image the shallow feeding system of the volcano that enabled such a rapid magma accumulation and discharge. The high-rate data allowed us to constrain the sources responsible for the fast and violent dynamics of each paroxysm, while the cumulated deformation measured by DInSAR and daily GNSS solutions, over a period of 12 days encompassing the entire eruptive sequence, also showed the deeper part of the source involved in the considered period, where magma was stored. We defined the dynamics and rates of the magma transfer, with a middle-depth storage of gas-rich magma that charges, more or less continuously, a shallower level where magma stops temporarily, accumulating pressure due to the gas exsolution. This machine-gun-like mechanism could represent a general conceptual model for similar events at Etna and at all volcanoes.

Secular variations of magma source compositions in the North Patagonian batholith from the Jurassic to Tertiary: Was mélange melting involved?

This study of Sr-Nd initial isotopic ratios of plutons from the North Patagonian batholith (Argentina and Chile) revealed that a secular evolution spanning 180 m.y., from the Jurassic to Neogene, can be established in terms of magma sources, which in turn are correlated with changes in the tectonic regime. The provenance and composition of end-member components in the source of magmas are represented by the Sr-Nd initial isotopic ratios (87Sr/86Sr and 143Nd/144Nd) of the plutonic rocks. Our results support the interpretation that source composition was determined by incorporation of varied crustal materials and trench sediments via subduction erosion and sediment subduction into a subduction channel mélange. Subsequent melting of subducted mélanges at mantle depths and eventual reaction with the ultramafic mantle are proposed as the main causes of batholith magma generation, which was favored during periods of fast convergence and high obliquity between the involved plates. We propose that a parental diorite (= andesite) precursor arrived at the lower arc crust, where it underwent fractionation to yield the silicic melts (granodiorites and granites) that formed the batholiths. The diorite precursor could have been in turn fractionated from a more mafic melt of basaltic andesite composition, which was formed within the mantle by complete reaction of the bulk mélanges and the peridotite. Our proposal follows model predictions on the formation of mélange diapirs that carry fertile subducted materials into hot regions of the suprasubduction mantle wedge, where mafic parental magmas of batholiths originate. This model not only accounts for the secular geochemical variations of Andean batholiths, but it also avoids a fundamental paradox of the classical basalt model: the absence of ultramafic cumulates in the lower arc crust and in the continental crust in general.

An Overview on the Classification and Tectonic Setting of Neoproterozoic Granites of the Nubian Shield, Eastern Desert, Egypt

Granites constitute the main rock components of the Earth’s continental crust, which suggested to be formed in variable geodynamics environments. The different types of granitic rocks, their compositional characteristics, tectonic settings and magma sources are outlined. Mineralogical classification of granites includes four rock types: tonalites, granodiorites, granite (monzogranite and syenogranites) and alkali-feldspar granites. Alphabetical classification subdivided granites into: I-type, S-type, A-type and M-type granites. Moreover, formation of granitic magmas requires distinctive geodynamic settings such as: volcanic arc granite (Cordilleran); collision-related granites (leucogranites); intra-plate and ocean ridge granites. The Eastern Desert of Egypt (ED) forms the northern part of Nubian Shield. Both older and younger granites are widely exposed in the ED. Old granites (OG) comprise tonalites and granodiorites of syn- to late-orogenic granitoid assemblages. They are calcalkaline, I-type, metaluminous and display island arc tectonic setting. Younger granites (YG) on the other hand, include granites, alkali-feldspar granites and minor granodiorites. They are of I- and A-type granites and of post-orogenic to anorogenic tectonic settings. The majority of the YG are alkaline, A-type granite and of within-plate tectonic setting (WPG). The A-type granites are subdivided into: A2-type postorogenic granites and A1-type anorogenic granites. Granite magma genesis involves: (a) fractional crystallization of mafic mantle-derived magmas; (b) anatexis or assimilation of old, upper crustal rocks (c) re - melting of juvenile mafic mantle – derived rocks underplating the continental crust. Generally, older I-type granitoids were interpreted to result from melting of mafic crust and dated at approximately 760–650 Ma, whereas younger granites suggested to be formed as a result of partial melting of a juvenile Neoproterozoic mantle source. Moreover, they formed from anatectic melts of various crustal sources that emplaced between 600 and 475 Ma.

Geomorphological and Spatial Characteristics of Underwater Volcanoes in the Easternmost Australian-Antarctic Ridge

Underwater volcanoes and their linear distribution on the flanks of mid-ocean ridges are common submarine topographic structures at intermediate- and fast-spreading systems, where sufficient melt supplies are often available. Such magma sources beneath the seafloor located within a few kilometers of the corresponding ridge-axis tend to concentrate toward the axis during the upwelling process and contribute to seafloor formation. As a result, seamounts on the flanks of the ridge axis are formed at a distance from the spreading axis and distributed asymmetrically about the axis. In this study, we examined three linearly aligned seamount chains on the flanks of the KR1 ridge, which is the easternmost and longest Australian-Antarctic Ridge (AAR) segment. The AAR is an intermediate-spreading rate system located between the Southeast Indian Ridge and Macquarie Triple Junction of the Australian-Antarctic-Pacific plates. By inspecting the high-resolution shipboard multi-beam bathymetric data newly acquired in the study area, we detected 20 individual seamounts. The volcanic lineament runs parallel to the spreading direction of the KR1 segment. The geomorphologic parameters of height, basal area, volume, and summit types of the identified seamounts were individually measured. We also investigated the spatial distribution of the seamounts along the KR1 segment, which exhibits large variations in axial morphology with depth along the ridge axis. Based on the geomorphology and spatial distribution, all the KR1 seamounts can be divided into two groups: the subset seamounts of volcanic chains distributed along the KR1 segment characterized by high elevation and large volume, and the small seamounts distributed mostly on the western KR1. The differences in the volumetric magnitude of volcanic eruptions on the seafloor and the distance from the given axis between these two groups indicate the presence of magma sources with different origins.