Magma Migration Research Articles

The magmatic plumbing beneath the Wudalianchi volcanic region, northeast China, is recharged by magma reservoirs under Keluo volcano

The intraplate magmatic systems of the Wudalianchi potassic volcanic region in northeast China have historically exhibited eruption activity, and there is evidence that they retain the potential for future eruptions. However, the characteristics of this complex volcanic region remain debated and poorly understood due to the relatively low resolution of existing crustal tomography. Here we present a high-resolution shear-wave velocity model of the crust and uppermost mantle obtained by applying full-waveform ambient noise tomography to data recorded at 106 recently deployed broadband seismic stations. Our results reveal two vertical low-velocity anomalies discovered beneath the Nuominhe and Keluo volcanoes, which appear to be magma ascending pathways connecting a mantle source and crustal reservoirs. A horizontally interconnected low-velocity zone in the middle-lower crust indicates a magmatic migration pathway transporting potassium-rich melts directly among the potassic volcanic units. In addition, the Wudalianchi volcano is likely recharged by deeper and larger crustal magma reservoirs under the Keluo volcano. The unique geometries of the magmatic plumbing systems provide valuable seismic evidence for the volcanic eruption dynamics and potential volcanic hazards.

Generation and evolution of the Choiyoi granitic magmatism based on U-Pb zircon studies, Cordón del Portillo, Frontal Cordillera (Argentina)

Based on a relevant geochronological U-Pb zircon dataset (n = 47) from a sample (PBL-109) of the Cerro Punta Blanca pluton (CPB), which is part of a calc-alkaline suite, we corroborate the development of protracted magmatic activity with three major crystallization events for this Permian magmatism: 278 ± 1, 283 ± 2, and 289 ± 2 Ma, which outcrop in the Cordón del Portillo, Cordillera Frontal (CF) of Argentina. These ages can be assigned to the lower section of the magmatic record of the Choiyoi magmatism (ca., 290-265 Ma), while the age of 289 Ma represents the oldest known age for the Choiyoi, indicating the start of this magmatism during the Artinskiense. Considering these geochronological data, we postulate the presence of a deep mush reservoir where protracted magmatic activity permitted the prolonged crystallization of antecrysts (ca. 283–289 Ma). Migration of the parental magma from the mush reservoir zone occurred near the time of emplacement and culminated in the formation of an ephemeral magma chamber at shallow levels, where zircon autocrysts crystallized (ca. 278 Ma). Age spectra reported within individual samples support the idea of massive magma migration when conditions were favorable (e.g., thermally matured crust). In this view, the studied “older” Choiyoi magmatism represents a continuous magmatic event lasting 11 Ma and corresponds to a single magmatic episode rather than different periods of magmatic activity and subsequent emplacements. A later alkali-calcic magmatic event is recorded at 265 ± 4 Ma from a sample of the Cerro Bayo pluton (MH-0113), which could represent the end of the lower section of the Choiyoi magmatism. Whole-rock Sm-Nd, Rb-Sr, and Lu-Hf data in zircon, along with ages reported for the studied igneous and inherited zircon from the CPB; together with isotopic data and ages from the detrital zircon found in the Carboniferous accretionary complex of Chile, indicate that the source of the Permian parental magma in this region was a heterogeneous continental crust mainly formed by Devonian and Carboniferous rocks, related probably to a magmatic arc. However, some contribution from the Carboniferous accretionary complex of Chile to the parental magma should be consider.

Time‐Resolved Trigger Processes Leading to the Plinian Eruptions at Sakurajima Volcano, Japan

AbstractMafic magma recharge of crustal reservoirs and subsequent magma mixing has been considered a direct trigger of volcanic eruptions. However, although recharge frequently occurs in many active volcanoes, it rarely leads to an eruption immediately, making its role as a trigger ambiguous. Sakurajima volcano, Japan, has vigorously erupted three times since the 15th century following a common process; mixed magmas after recharge were once stored in a shallow, thick conduit before each eruption (conduit pre‐charge). We reconstructed the magma migration with a high time resolution by diffusion modeling on orthopyroxene and magnetite. Orthopyroxene phenocrysts recorded prolonged diffusive re‐equilibration timescales of years or more after recharge‐and‐mixing. Magnetite, which has the fastest elemental diffusivity among the phenocrysts examined, predominantly lacks zoning. This demonstrates that the mineral phase was re‐equilibrated with surrounding magma and homogenized via elemental diffusion after the final magmatic perturbation, implying the final repose of the shallow pre‐charged magma body for more than several tens of days. After this shallow stagnation period, the Plinian magmas began to ascend and reached the surface within 55 hr. Mass balance calculations show that crystallization‐driven vesiculation upon pre‐charge can produce overpressure sufficient to cause an eruption. The Sakurajima cases demonstrate the hierarchical timescales of trigger processes leading to the explosive eruptions.

Translithospheric magma plumbing system fossilized in the Emeishan large igneous province

Lighting up the magma plumbing system beneath a large igneous province (LIP) is challenging because the complex magma migration paths are often covered by flood basalts and sediments. Here, we present a three-dimensional seismic image of the Permian Emeishan LIP in Southwest China, constructed by joint inversion of Rayleigh wave dispersions and receiver functions. The results outline a cylindrical, high-velocity anomaly extending to ~135 km depth below the inner zone of this continental LIP. The geometry and magnitude of the high-velocity anomaly suggest that it represents culminated crystallized materials of primary magmas, thereby mirroring a magma plumbing system fossilized in the lithosphere. Furthermore, our geodynamic models illustrate that the nearly vertical plumbing system was controlled by slow plate motion during the magma emplacement. The plume head beneath a nearly static plate has higher thermal buoyancy and thus promotes more intensive magma emplacement. This phenomenon may apply to other LIPs throughout Earth’s history

Tracing Lower Crustal Contamination in Continental Arc Magmas Using Sr–Nd–Hf Isotopes: A Combined In Situ and Bulk Rock Approach Applied to the Adamello Batholith

Abstract The incremental construction of plutons characterises magmatic activity in arc settings, where new continental crust is produced. This polyphasic growth entails interactions with one or more crustal components, which modulate the geochemical and isotopic compositions of the newly formed crust. However, the early stages of magmatism are not always preserved due to obliteration by later magmatic pulses. Spatial migration of magmatism during the construction of the Adamello batholith (Northern Italy) enables the examination of the early pulses of pluton formation, thus allowing a time-integrated study of the relative importance of crystallisation-differentiation and contamination in a continental arc setting. We conducted a detailed textural, major and trace element and Sr isotopic study of plagioclase from the first intrusive pulses of the Adamello batholith, combined with new major, trace element and Sr–Nd isotopic analyses of bulk rock samples across the entire Adamello batholith. We selected well-characterised samples with published CA-ID-TIMS 206Pb-238U ages and Hf isotopic composition for zircons. Strontium isotopes in plagioclase from the same samples were determined by laser ablation multi-collector ICP-MS. The tonalitic samples in the early magmatic stages show elevated but constant Sr isotopic compositions despite large variations in anorthite contents (An90 to An13), indicating that crustal contamination occurred before significant differentiation. Invariant bulk-rock 87Sr/86Sr with variable SiO2 in all superunits of the Adamello batholith further supports contamination preceding significant melt differentiation. Contamination by lower crustal basement lithologies is due to the increasing thermal anomaly triggered by consecutive magmatic injections coupled with the heterogeneous and less restitic nature of the basement in the early stage of the magmatic system (i.e., before consumption of fusible components). In addition, we observe significant variability in crustal contamination proxies (e.g., 87Sr/86Srplag, 87Sr/86Srbulk, εNdbulk, εHfzircon) during the initial phases of magmatism. This variability likely reflects the uneven distribution of positive thermal anomalies in the lower crust during early magmatic stages as well as the diverse lithological and isotopic makeup of the lower crust. The processes identified in our case study are pertinent to continental arc magmatism, particularly where magmas interact with a metapelitic lower crust.

Optimal spatial distribution of seismic stations to detect magma migration using the seismic amplitude ratio analysis

Magma migrations frequently trigger seismic swarms, resulting in seismic events that overlap in time and hinder real-time phase picking commonly used for hypocenter location. Addressing this challenge, seismic amplitude ratio analysis (SARA) allows identification of seismic migrations in real-time by simply tracking the relative seismic amplitude between a pair of seismic stations.This paper aims to identify key statistical features of the seismic network array locations that improve their ability to detect seismic migrations using SARA. We evaluated the capability to detect the most frequently oriented magma migrations in over 100 volcanoes, using a criterion previously proposed to study vertical magma migrations in Piton de la Fournaise. Additionally, we investigate the influence of vent-station proximity on magma conduit coverage and identify the distance ratio that yields improved detection.Furthermore, we estimate the seismic network efficiency by calculating the detection capability volume per station. We then use the random forest regression algorithm to identify which statistical features of the seismic network location contribute more to the efficiency disparity among different volcanoes. Notably, our findings reveal that optimizing seismic network coverage entails maximizing the standard deviation of relative pair station distances, while maintaining a prescribed minimum separation distance between station pairs. Our results reveal important criteria that can be used to optimize seismic network location design.

Evolution of Devonian and Carboniferous pre-Andean arc magmatism in the Frontal Cordillera (27°-35°S), Argentina: Insights from U–Pb zircon and isotopic studies

We propose a conceptual model for the magmatic system that gave rise to the Devonian and Carboniferous subduction-related magmatism in the Frontal Cordillera of Argentina, (representing the pre-Andean margin of SW Gondwana), integrating previous geochronological studies with new geochronological data from granitoids and subvolcanic dikes, which include reported petrological, geochemical, and isotopic data. This conceptual magmatic system postulates an extended magmatic activity, similar to that reported in a previous study for the Devonian foreland magmatism located in the present-day Sierras Pampeanas of Argentina. The geochronological data play a central role in the model, and lead us to postulate the presence of a deep mush reservoir, where a prolonged magmatic activity, permitted the extended crystallization of Devonian (ca. 400 ± 3 and 414 ± 3 Ma) and Carboniferous antecrysts (334 ± 2 Ma, 341 ± 2, and 348 ± 2, where the age of 348 Ma maybe considered as derived from the source, i.e. xenocrysts). Migration of the parental magma from the mush reservoir zone occurred near the time of emplacement and culminated in the formation of an ephemeral magma chamber located in shallow levels, where Devonian and Carboniferous zircon autocrysts crystallized (382 ± 5 Ma and 325 ± 2 Ma, respectively). Age spectra reported within individual sample, favors the idea of a massive migration of magma when conditions were favorable (e.g., thermally matured crust). Notably, the Carboniferous geochronological data is strongly consistent with previous geochronological data recently reported for the granitic rocks of the Tabaquito batholith. Although this model supports a protracted magmatic system for the Devonian and Carboniferous subduction-related magmatism, previous and new isotopic data reveal a significant difference in the petrogenetic processes. While the Carboniferous arc magmatism shows significant juvenile material contribution, this is not the case for the Devonian magmatic arc, where isotopic data suggest an older continental lithosphere as the dominant source. These different sources are attributed to two contrasting geodynamic settings, with advancing and retreating oceanic slabs, respectively.

Uplift of the Tibetan Plateau driven by mantle delamination from the overriding plate

The geodynamic evolution of the Tibetan Plateau remains highly debated. Any model of its evolution must explain the plateau’s growth as constrained by palaeo-altitude studies, the spatio-temporal distribution of magmatic activity, and the lithospheric mantle removal inferred from seismic velocity anomalies in the underlying mantle. Several conflicting models have been proposed, but none of these explains the first-order topographic, magmatic and seismic features self-consistently. Here we propose and test numerically an evolutionary model of the plateau that involves gradual peeling of the lithospheric mantle from the overriding plate and consequent mantle and crustal melting and uplift. We show that this model successfully reproduces the successive surface uplift of the plateau to more than 4 km above sea level and is consistent with the observed migration of magmatism and geometry of the lithosphere–asthenosphere boundary resulting from subduction of the Indian plate and delamination of the mantle lithosphere of the Eurasian plate. These comparisons indicate that mantle delamination from the overriding plate is the driving force behind the uplift of the Tibetan Plateau and, potentially, orogenic plateaus more generally.

The structurally‐controlled polyphase Au mineralization in the southern Youjiang Basin: Insights from magnetotelluric imaging

The Youjiang Basin, located on the southwestern margin of the Yangtze Craton, hosts significant Upper Triassic and Lower Cretaceous gold (Au) deposits. Isotope studies suggest that the polyphase Au mineralization in the Youjiang Basin is derived from Indosinian and Yanshanian magmatic events, but the structural plumbing system beneath the basin remains poorly understood. We collected 32 magnetotelluric (MT) stations extending from the southern margin of the Youjiang Basin which hosts the Zhesang, Badu, and Jinya Au deposits at the western margin of the Jiangnan Orogenic (JNO) belt. The full impedance tensor of MT data is utilized to generate a high-precision lithospheric resistivity model by three-dimensional (3-D) inversion. Importantly, our MT data reveals several crustal conductors (<30 Ω·m) underlying major structures in the northwestern part of the study area. These crustal conductive anomalies spatially overlap with high-magnetic anomalies and polymetallic Au-Sn deposits. We interpret these anomalies as solid conductive materials such as graphite and/or sulfides related to the structurally controlled Mesozoic Au mineralization. Additionally, we interpret these conductive anomalies beneath the known Au deposits as the conduits of magmas and solidified paleo-magma chambers that were emplaced along deep-seated northwestern faults and likely documenting the genetic relationship between Au mineralization and Mesozoic magmatism. The synthesis of the electrical resistivity structure and metallogenic considerations is presented in a schematic sketch model illustrating that the deep-seated northwest-trending structures likely acted as a plumbing system controlling both the migration and emplacement of magmas and related oxidizing fluids that subsequently precipitated the polyphase Au deposits.

Evolution of Silurian to Devonian magmatism associated with the Acadian orogenic cycle in eastern and southern Newfoundland Appalachians: Evidence for a three-stage evolution characterized by episodic hinterland- and foreland-directed migration of granitoid magmatism

The migration and character of magmatism over time can provide important insights into the tectonic evolution of an orogen. We present evidence for three separate stages of compositionally distinct granitoid magmatism associated with the Acadian orogenic cycle in the eastern and southern Newfoundland Appalachians. The interpretations are based on new zircon U-Pb ages, geochemical data, and Sr-Nd-Hf-O isotopic data for 18 samples from 15 Silurian and Devonian granitoid plutons, combined with previously published data. The three stages outline hinterland- and foreland-directed migration trends and represent subduction (435−420 Ma), syncollision (415−405 Ma), and postcollision (395−370 Ma) settings in the Acadian orogenic cycle. The Silurian plutons (435−420 Ma) of the first stage consist mainly of quartz diorite, tonalite, granodiorite, monzogranite, and syenogranite, with high-K calc-alkaline and enriched Sr-Nd-Hf-O isotopic compositions (e.g., εNd[t] = −5 to −2; εHf[t] = −3 to −1; δ18O = +6‰ to +8‰). They are interpreted to record the subduction of oceanic lithosphere of the Acadian seaway that separated the leading edge of composite Laurentia, represented by the Gander margin, and Avalonia. Early Devonian plutons (415−405 Ma) of the second stage contain more voluminous monzogranite and syenogranite; they have calc-alkaline to high-K calc-alkaline features, adakite-like compositions, and more depleted Sr-Nd-Hf-O isotopic compositions (e.g., εNd[t] = −6 to 0; εHf[t] = +1 to +3; δ18O = +5‰ to +6‰). Plutons of this stage occur mostly to the northwest of the Silurian granitoids, indicating a regional-scale northwestward (hinterland-directed) migration of magmatism with a rate of &amp;gt;9 km/m.y. The migration is interpreted to have been related to the progressive shallow underthrusting of Avalonia beneath the Gander margin (composite Laurentia) at least as far as 90 km inboard. The Middle to Late Devonian plutons of the third stage (395−370 Ma) consist mainly of monzogranite, syenogranite, and alkali-feldspar granite, which are silica- and alkali-rich granites with large negative Eu anomalies. These rocks are concentrated along both sides of the Dover−Hermitage Bay fault zone, which represents the boundary between Avalonia and composite Laurentia, to the southeast of the Silurian and Early Devonian igneous rocks. This stage of magmatism represents a foreland-directed (retreating) migration. The Early Devonian and Middle to Late Devonian episodes of magmatism were separated by a gap between 405 Ma and 395 Ma and recorded an evolution from (high-K) calc-alkaline to alkaline compositions, ascribed to partial delamination of Avalonian lithospheric mantle in a postcollisional setting.

Magma intrusion and migration under the continental large igneous provinces

Large igneous provinces (LIPs) are often the surface expressions of mantle plume process. Intensive magma could occur in the continental crust when a mantle plume penetrates the Moho. How magmas migrate in the crust, forming LIPs, remains debated. In this study, we employ 2D thermo‐mechanical numerical modelling to study the dynamics of magma intrusion and migration within the continental crust affected by a mantle plume. Our results suggest that (1) lateral magma migration dominates crustal deformation, promoting the possible formation of metamorphic core complexes; (2) two distinct crustal deformation modes are recognized regarding the presence or absence of crustal break‐up, affected by whether a significant vertical magma migration occurs; (3) the crustal thickness, Moho temperature and the length of the preset weak crust–mantle decoupling zone are the three key parameters controlling crustal magma migration. This study highlights the importance of lateral magma migration in the lower crust and provides physical mechanisms for the interpretation of magma migration under LIPs.

A conceptual model for the evolution of a magmatic system based on U-Pb zircon studies on the devonian granites of Sierra de San Luis, Argentina

ABSTRACT It is widely acknowledged that thermal models clearly demonstrate that crustal magma bodies should solidify rapidly upon emplacement. Small plutons can cool below the solidus in thousands of years, while even large plutons require hundreds of thousands of years, but not more than a million years. However, recent U-Pb zircon geochronological data contradict these results, strongly suggesting that magmatic systems are often characterized by protracted events. Therefore, a conceptual framework that reconciles thermal models and the geochronology data is necessary. The Devonian foreland magmatism of the Sierra de San Luis is made up of two distinctive suites, the Monzonite suite (<65 wt.% SiO2) and the Granite suite (>65 wt.% SiO2), both emplaced at ca. 3.7 kbar. Classification of the studied Devonian granitoids is debatable because they have a hybrid I- to A-type granite signature. Based on a robust geochronological U-Pb zircon dataset we corroborate the development of a protracted magmatic activity with three major crystallization events for this Devonian magmatism: 391 ± 1, 384 ± 1, and 379 ± 2 Ma. Considering these geochronological data, we postulate the presence of a deep mush reservoir, where a prolonged magmatic activity, permitted the prolonged crystallization of antecrysts (ca. 395–384 Ma). Migration of the parental magma from the mush reservoir zone occurred near the time of emplacement and culminated in the formation of an ephemeral magma chamber located at shallow levels, where zircon autocrysts crystallized (ca. 379 Ma). Age spectra reported within individual samples support the idea of a massive migration of magma when conditions were favourable (e.g. thermally matured crust). Individual crystallization ages recorded by monazites hosted in two-mica granites are comparable to those obtained from zircons, supporting the presence of a long-lasting hot source. Additional geochronological data indicate that a later thermal event (ca. 353 Ma) could have partially affected some areas of the Devonian magmatic zircons promoting Pb loss, with subsequent partial resetting of the isotopic clock.

Tracing Magma Migration at Mt. Etna Volcano during 2006–2020, Coupling Remote Sensing of Crater Gas Emissions and Ground Measurement of Soil Gases

The geochemical monitoring of volcanic activity today relies largely on remote sensing, but the combination of this approach together with soil gas monitoring, using the appropriate parameters, is still not widely used. The main purpose of this study was to correlate data from crater gas emissions with flank emissions of soil gases at Mt. Etna volcano from June 2006 to December 2020. Crater SO2 fluxes were measured from fixed stations around the volcano using the DOAS technique and applying a modeled clear-sky spectrum. The SO2/HCl ratio in the crater plume was measured with the OP-FTIR technique from a transportable instrument, using the sun as an IR source. Soil CO2 efflux coupled with the 220Rn/222Rn activity ratio in soil gases (named SGDI) were measured at a fixed monitoring site on the east flank of Etna. All signals acquired were subject both to spectral analysis and to filtering of the periodic signals discovered. All filtered signals revealed changes that were nicely correlated both with other geophysical signals and with volcanic eruptions during the study period. Time lags between parameters were explained in terms of different modes of magma migration and storage inside the volcano before eruptions. A comprehensive dynamic degassing model is presented that allows for a better understanding of magma dynamics in an open-conduit volcano.

Interaction between magmatism and polygonal faults revealed by three-dimensional seismic data in the Zhongjiannan Basin, South China Sea

The features of magmatism and polygonal faults are unveiled through a comprehensive analysis of high-resolution and reprocessed three-dimensional (3D) seismic data covering 1200 km2 in the Zhongjiannan Basin, South China Sea. These structures influence the hydrocarbons formation and migration, and the interaction between them is significant for hydrocarbon exploration. We identified magmatic bodies, including eight volcanoes, fifty sills, and a single igneous diapir, are mainly developed in three stages: 27.9–23 Ma (between T70 and T60), around 11.6 Ma (T40), and after 5.3 Ma (T30). The polygonal faults, characterized by linear, polygonal, and radial configurations, are extensively distributed within Miocene and Quaternary sediments. In the platform, polygonal faults manifest linear or polygonal attributes, while in the basin floor, they show a polygonal shape. Particularly, above volcanoes and the igneous diapir, the polygonal faults exhibit radial configurations. However, among volcanoes, few polygonal faults are discerned based on the coherence attribute slice. The controlling factors influencing the distribution of polygonal faults are the lithological disparity, temperature, and pressure. Radial fault formation is promoted by volcanoes and igneous diapir through differential compaction. The resultant vertical pathways formed by volcanoes and overlying radial faults play a crucial role in adjusting pressure differentials or water discharges. Consequently, these pathways disrupt the requisite conditions for the formation of polygonal faults among volcanoes. In the basin floor, magma intrusion occurs along the planes of polygonal faults, giving rise to inclined, spoon-shaped, and bowl-shaped sills. As weak zones, polygonal faults worked as connectors for sills at different depths. Hence, pre-existing polygonal faults are considered pathways facilitating magma migration from the lower Miocene to Pliocene strata, ultimately contributing to sill formation in shallow sediments. This study advances our understanding of the intricate interaction between polygonal faults and magmatism within sedimentary basins.

Tectonic and Crustal Processes Drive Multi-Million Year Arc Magma Evolution Leading up to Porphyry Copper Deposit Formation in Central Chile

Abstract Subduction zone magmatism is a major control of volcanism, the generation of modern continental crust and the formation of economically important porphyry Cu–(Mo–Au) deposits. Reading the magmatic record of individual arc segments and constraining the rates of magmatic changes are critical in order to fully understand and quantify the processes that drive magma evolution in subduction settings during arc growth. This study focuses on the San Francisco Batholith and the Rio Blanco-Los Bronces porphyry deposit cluster in central Chile, which provides an igneous rock record over ~13.5 Myr of arc evolution. We use whole-rock geochemistry, zircon geochronology and Hf isotope geochemistry to track changes in the crustal magmatic system of this arc segment during crustal thickening and porphyry Cu deposit formation. By combining the analytical dataset with Monte Carlo fractional crystallisation and assimilation fractional crystallisation modelling, we test a model for significant crustal involvement during magma evolution. Systematic and continuous increases in Dy/Yb, La/Yb, V/Sc and Sr/Y in the magmas over time indicate a transition in the main fractionation assemblage from plagioclase-dominated to amphibole-dominated that reflects deeper crystallisation and/or a higher meltwater content. Concomitant decreases in εHf and Th/La as well as increasing Ba/Th are best explained by assimilation of progressively deeper crustal lithologies from low (Chilenia) to high Ba/Th (Cuyania) basement terranes. Our study highlights that an increasingly hydrous magma and a deepening locus of crustal magma differentiation and assimilation, driven by crustal thickening contemporaneous with increased tectonic convergence and ingression of the aseismic Juan Fernandez ridge, can account for all investigated aspects of the multi-Myr magmatic evolution leading up to the formation of the Rio Blanco-Los Bronces porphyry Cu deposits. Our findings corroborate the importance of high-pressure differentiation of hydrous magma for the formation of Andean-style porphyry deposits. Once magmas favourable for porphyry Cu mineralisation were generated in the lower crust, multiple episodes of efficient magma migration into the upper crust fed several, discrete, shallow magmatic-hydrothermal systems over ~3.5 Myr to form the world’s largest known Cu resource at Rio Blanco-Los Bronces.

Early Cretaceous continental arc magmatism in the Wakhan Corridor, South Pamir: Mantle evolution and geodynamic processes during flat subduction of the Neo-Tethyan oceanic slab

The petrogenesis of continental arc magmas provides critical insights into thermal evolution and geodynamics of the continental lithosphere, crust-mantle interaction, and deep dynamic processes. In this study, we report new zircon U-Pb ages along with isotopic and elemental whole-rock geochemistry, mineral chemistry, and Hf-O isotope data for the Kalaqigu diorites and monzogranites of the Chinese Wakhan Corridor, South Pamir. Zircon U-Pb dating indicates that the Kalaqigu pluton was emplaced in the Early Cretaceous (ca. 108−106 Ma). The diorites are geochemically characterized by low SiO2 (51.9−54.5 wt%) and CaO (7.7−9.4 wt%) contents, but high MgO (5.3−8.3 wt%), Al2O3 (12.8−16.8 wt%), and TiO2 (0.6−1.1 wt%) contents as well as high Mg# (56−65) values. Thus, they are similar to high-Mg diorites: enriched in large ion lithophile elements (e.g., K, Sr, and Ba) and light rare earth elements, while depleted in high field strength elements (i.e., Nb, Ta, Zr, and Hf). Combined with negative εNd(t) (−6.9 to −14.0) and εHf(t) (−9.9 to −12.2), and high (87Sr/86Sr)i (0.7075−0.7086) ratios, these observations indicate that they originated from an enriched lithospheric mantle source. High δ18Ozrn (7.49‰−9.01‰) values, in conjunction with relatively high 207Pb/206Pb and 208Pb/206Pb ratios, suggest that the source was modified by subducted sediment-derived melts. Variable Cr contents (54−117 ppm) are likely controlled by minor fractionation of olivine and orthopyroxene. The monzogranites show high SiO2 contents (69.2−72.0 wt%), and low Rb/Sr (0.4−0.6), (K2O + Na2O)/CaO (2.6−4.8), and FeOT/MgO ratios (2.6−3.2). They contain diagnostic cordierite and show strongly peraluminous characteristics (A/CNK &amp;gt; 1.1) with high δ18Ozrn (7.82‰−8.85‰) values that are compatible with those of typical S-type granites. Their abundant inherited zircons, with age populations similar to those of detrital zircons from regional early Paleozoic metasedimentary rocks, indicate that they were derived from partial melting of ancient metasedimentary rocks. Phase equilibrium modeling is consistent with biotite-dehydration melting of metagreywacke, probably at ∼750 °C and ∼6.0 kbar, as indicated by the biotite chemistry. Based on regional geochronology, a south-to-north magmatic migration suggests that northward flat-slab subduction of the Neo-Tethyan oceanic slab played an important role in the generation of these widespread Early Cretaceous continental arc magmatic rocks. However, the granitoids were generated earlier than the mantle-derived mafic rocks, which suggests that crustal melting occurred during the early stage of subduction. The continuous flat-subduction resulted in partial melting of subducted sediments, which metasomatized the mantle wedge. Contemporaneous regional compression primarily occurred far north of the subduction zone (i.e., North and Central Pamir), inducing deformation as well as crustal shortening. With the flare-up of continental arc magmatism in South Pamir, crustal shortening moved southward. These processes, combined with the addition of voluminous, mantle-derived magmas, played an important role in crustal thickening in Pamir during the Early Cretaceous.

Pre-existing structural control on the recent Holuhraun eruptions along the Bárðarbunga spreading center, Iceland

The active rift zones in Iceland provide unique insight into the geodynamic processes of divergent plate boundaries. The geodynamics of Iceland are studied intensively, particularly, by geophysical methods sensitive to active and/or visible structures such as earthquake seismic and Synthetic Aperture Radar observations or aerial photographs. However, older and less active structures, that may exert a strong control on the presently active geodynamics, are often buried beneath recent volcanic or sedimentary deposits and are—due to their passive mode—overseen by the typical geophysical investigations. Aeromagnetic surveys provide spatial information about subsurface magnetization contrasts relating to both active and inactive structures. However, the aeromagnetic data in Iceland were collected in the 1970-80s and are relevant only to large-scale regional rift studies. With the availability of reliable drones and light-weight atomic scalar sensors, high-quality drone magnetic surveys can provide an unprecedented spatial resolution of both active and passive structures of rift systems as compared to conventional airborne surveys. Here, we present the results of a drone-towed magnetic scalar field and scalar gradiometry study of the north-northeast trending Bárðarbunga spreading center to the north of the Vatnajökull ice cap, Iceland. Our results provide new information about the structural complexity of rift zones with evidence of densely-spaced, conjugate and oblique faults throughout the area. Evidence is shown of a hitherto unknown and prominent east-northeast trending fault structure that coincides with the northern tip of the main eruption edifice of the 1797 and 2014-15 Holuhraun volcanic events. We suggest that this pre-existing structure controlled the locus of vertical magma migration during the two Holuhraun events.

Rapid accumulation and ascent precedes caldera forming eruption of low viscosity magma

Mafic magma is commonly associated with effusive eruptions, however several mafic volcanoes throughout the globe have produced explosive eruptions. Here we present one such volcano – Colli Albani. Colli Albani is 20 km SE of Rome and produced seven large volume ignimbrites. Field observations, mineral chemistry, and Sr and Nd isotopes in clinopyroxene show that the high potassic, silica undersaturated and CO2\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$_{2}$$\\end{document}-rich magma typical of Colli Albani is produced by partial melting of a metasomatized mantle. Clinopyroxene based thermobarometry combined with thermal modelling, indicates rapid accumulation of magma into the shallow crust preceding the last caldera forming event (355 ka). The crystallization of high Mg# and high Cr2\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$_{2}$$\\end{document}O3\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$_{3}$$\\end{document} clinopyroxenes at low pressures and high temperatures indicates rapid magma ascent from the mantle. We suggest that a final rapid input of this deeply sourced magma destabilised the shallow and fast assembled magma reservoir and lead to the caldera forming event. Our findings have significant implications for the evaluation of the timescales of reactivation of Colli Albani and other similar long-quiescent calderas erupting low viscosity magmas, as rapid migration of magma to shallow reservoirs may result in short unrest periods prior to a large eruption.

Crustal structure of the southern Lower Yangtze region and its geological implications: a deep reflection seismic profile

Abstract The southern part of the Lower Yangtze Region, located in the area bounded by the Jiangnan and Dabie orogenic belts, is characterized by a complex fault structure and unique magmatism, and is one of the key areas in regional tectonic evolution and magmatism. In this study the authors used a 45-km-long deep reflection seismic profile to study the crustal structure of the southern Lower Yangtze Region. The results show that the crust in the study area is ∼31.5–33.6 km thick, with the Huaining Basin as the thinnest part. A detachment surface at a depth of around 12 km divides the crust into the upper and lower crust, act as a key tectonic decoupling layer. Two sets of thrust nappe faults with opposite tendencies have horizontally developed above the detachment surface, with the Huaining Basin as the centre. A typical wedge-shaped structure has developed below the Qianshan Basin, suggesting that two sets of thrust nappe faults may have formed in the same dynamic system. The reflection patterns of the lower crust on both sides of the profile exhibit a monoclinic feature with opposite tendencies, while the lower part of the Huaining Basin exhibits an arcuate upward-arch feature. A prominent ductile shear zone in the lower crust has developed on the west side of the Huaining Basin, which can provide a favourable channel for magma migration. The results of this study deepen our understanding of the deep structure of the Lower Yangtze Region, and provide important constraint data for research on dynamic mechanisms.

Magma migration below Tancítaro and Paricutin volcanoes revealed by seismology

SUMMARY In this study, we analyse the 2020 seismic swarm that lasted two months and occurred between the Tancítaro and the Paricutin volcanoes in the Michoacán Guanajuato Volcanic Field, Mexico. We developed a new method to automatically detect and locate about 100 000 earthquakes, enabling us to track the magma migration through narrow dykes. Additionally, we reveal the presence of two magma reservoirs from two seismic noise tomography results. The first reservoir is located from 8 to 20 km below sea level and beneath the Tancítaro volcano and probably corresponds to a complex network of dykes and sills. This crustal reservoir is fed by a mantle reservoir with a wide horizontal extension between 35 and 50 km below sea level. The seismic swarm initiated beneath the Tancítaro summit in the lower portion of the crustal magma reservoir. At this stage, the seismicity migration was mainly horizontal, which we interpret as its response to the higher normal stress caused by the gravitational load of Tancítaro. Once the magma was displaced laterally from beneath Tancítaro, magma migration became more vertical. The swarm reached the upper portion of the crustal magma reservoir but did not escape it. We also reveal the effect of a distant but strong tectonic earthquake on the seismic swarm. Before its occurrence, magma migration followed several paths; afterwards, it became more focused along a single path. Finally, after the swarm, we observed a second type of seismicity called post-swarm seismicity, with a lower earthquake rate but with higher magnitudes. The hypocentres were diffuse and horizontally centred on the previous swarm location. Furthermore, some earthquakes were aligned along shallow faults, generating a high seismic risk to the different Tancítaro nearby localities.