Volcanic Centers Research Articles

Strain Localization at Volcanoes Undergoing Extension: Investigation of Long‐Term Deformation at Krafla and Askja Volcanic Systems in North Iceland

AbstractVolcanoes in extensional environments may show gradual subsidence over decades during quiescent periods, due to various processes such as magma withdrawal, cooling, contraction, plate spreading and viscoelastic response. If significant rheological anomalies reside in volcano roots, due to the presence of magma and hot rock, they can influence the style of deformation. We use Finite Element Method (FEM) models to explore how strain localization due to extension can lead to volcano deflation. We apply rheological models comprising an elastic layer overlying a viscoelastic domain and include local up‐doming regions of low viscosity material beneath volcanic centers. The models reveal a localized subsidence above the rheological anomaly, influenced by the tectonic extension, and by the up‐doming volume and its viscosity. The models suggest that plate divergence may account for 4–5 mm/yr of observed subsidence at Krafla and Askja volcanic systems (KVS and AVS, respectively) in North Iceland.

Boron isotopes identify deep-slab serpentinite in the source of Aleutian arc magma

Seafloor lavas of the Western Aleutian arc have isotopically heavy boron (δ11B to +13.4‰) that is negatively correlated with B content (ppm). Endmember samples are primitive dacites and rhyodacites (δ11B > +10‰, SiO2 = 63%−70%, Mg# > 0.60) with adakitic trace-element and isotopic characteristics that require roles for residual garnet and rutile in their formation. The source of isotopically heavy B is likely serpentinite in the mantle section of the subducting plate, which dewaters into an inverted geothermal gradient and drives melting within the overlying volcanic section at depths where prior effects of seawater alteration were minimal. Most volcanic rocks from the Aleutian Island locations have 10−30 ppm B with an average δ11B of ∼+1.0‰ ± 1.3‰, reflecting a mixed source dominated by subducted sediment. A subset of island samples has B that is isotopically light (δ11B < −2.4‰) and at low concentrations (<11.0 ppm), which is typical of arc lavas globally from rear-arc settings where depth-to-slab is high, and where δ11B may be interpreted to reflect a source in dehydrated (isotopically light) altered oceanic crust. Mass balance modeling indicates that isotopically heavy B from deep-slab serpentinite is present in the Aleutian source arc-wide but is typically masked by sediment-derived B at volcanic centers outside of the westernmost segment of the arc.

High-precision U–Pb geochronology of the Lundy igneous complex: implications for North Atlantic volcanism and the far-field Paleocene–Eocene ash record

The Lundy Island (Bristol Channel, UK) granite is a felsic expression of the southernmost igneous centre of the North Atlantic Igneous Province that emplaced millions of cubic kilometres of magma during the Paleogene. The granite's distinctive S-type, peraluminous, two-mica ± garnet ± tourmaline composition has led to the hypothesis that eruptions from the Lundy volcanic centre may be the source of thick felsic ash layers within the early Eocene Fur Formation (Denmark) that act as key marker horizons for the onset and duration of the Paleocene–Eocene Thermal Maximum. This paper presents high-precision zircon U–Pb emplacement ages of 57.24 ± 0.11/0.12/0.13 Ma for the granite and 55.970 ± 0.021/0.030/0.070 Ma for a felsic ‘lundyite’ dyke. Trace and rare earth element patterns indicate close similarities between late-stage Lundy activity and ash layer −33 in Denmark. This ash layer was deposited during the Paleocene–Eocene Thermal Maximum carbon isotope excursion, suggesting that the Lundy volcanic centre is likely to be the source of this key ash horizon and that magmatism at Lundy likely continued into the early Eocene.

Insights into the North Patagonian Massif Lower Crust: Petrology and Microstructure of Granulite Xenoliths

Abstract The continental lower crust constitutes a key zone for understanding the mantle–crust magmatic and mechanical transfers, but its study is hampered by the paucity of lower crust samples. Here, we characterise the petrological, geochemical and petrophysical processes structuring the lower crust of the North Patagonian Massif (NPM; Argentina) using a suite of representative mafic granulite and websterite xenoliths. These xenoliths were entrained by alkaline lavas from five volcanic centres that erupted between the Oligocene and Pleistocene. Electron microprobe and Laser Ablation Inductively Coupled Plasma Mass Spectrometer (LA-ICPMS) were used to obtain in situ geochemical data on the minerals, while microstructural data were obtained by Electron BackScatter Diffraction (EBSD). Both granulites and websterites display a granoblastic texture and sometimes a weak inherited magmatic layering. Mafic granulite xenoliths show a plagioclase + clinopyroxene ± orthopyroxene assemblage commonly associated with spinel or titanomagnetite. Websterite xenoliths show an association of clinopyroxene + orthopyroxene + spinel, along with accessory plagioclase. Mafic granulites and websterites have SiO2 contents ranging from 44 to 53 wt %, while their Mg# varies from 53 to 79. Clinopyroxenes are characterised by weak convex upward chondrite-normalised Rare Earth Elements - REE patterns (Light-REE </<< Mid-REE > Heavy-REE) which are similar to clinopyroxene phenocrysts and megacrysts from intra-plate basalts. Calculated liquid in equilibrium with clinopyroxene have similar REE patterns to those found in Cenozoic basalts from the NPM, suggesting that the xenolith suite represents evidence for underplating processes, possibly related to one of the magmatic events that have occurred in the NPM since the Permo-Trias. Mafic granulites and websterites show a weak mineral shape preferred orientation and an associated weak Crystal Preferred Orientation (CPO) related to the magmatic layering. Recorded plastic deformation is associated with the activation of both (100)[001] and (001)[100] slip systems in clinopyroxene, (100)[001] in orthopyroxene and (010)[001] in plagioclase. However, the activation of slip systems is generally not correlated with CPO in granulites, suggesting that the lower crust underwent subsolidus equilibration and weak plastic deformation in an inactive tectonic context, thereby preserving an inherited magmatic layering. Two-pyroxene (Fe–Mg) thermometer and pseudosection calculations define P–T conditions of the main paragenesis at 760°C to 1120°C and 7.2 to 10.3 kbar, which allows to define the Cenozoic geotherm of the NPM crust at 30°C/km and to reconsider the petrologic Moho depth at ca. 40 km.

Melt embayments record multi-stage magma decompression histories

Over the last decade, the melt embayment has proven its merit as a robust petrological tool capable of recording magma decompression rates for explosive eruptions. However, the models developed and applied to extract this information from embayments have not accounted for the complexity and nonlinearity of magma flow in the conduit. We present Embayment Decompression in Two Stages (EDiTS): a numerical model for extracting magma decompression rates from measured volatile diffusion profiles preserved in crystal-hosted embayments, approximating magma acceleration using two constant-rate decompression paths. This model solves for three unknown parameters: initial (deeper) and final (shallower) decompression rates, as well as the pressure where a transition occurs. We successfully benchmark EDiTS against existing numerical diffusion models, and use controlled multi-stage decompression experiments on natural quartz-hosted embayments to test the ability of our model to recover known decompression paths. We find that EDiTS is able to closely approximate the known two-stage path in the mixed-volatile (H2O + CO2) experiment, while a constant-rate modeling approach is unable to simultaneously fit H2O and CO2 gradients. However, in the H2O-saturated experiment, there is no unique solution to the resulting gradient, with both constant-rate and two-stage models reproducing the measured profile, and EDiTS notably overestimating the known total ascent time by several hours. Using decompression experiments, we show that constant-rate models can provide misleadingly good fits to embayment H2O gradients produced by more complex decompression histories, and thus the measurement and modeling of multiple diffusing species, when available, can provide crucial constraints. We then apply EDiTS to re-evaluate mixed-volatile embayment datasets from explosive silicic arc and caldera-forming eruptions from five volcanic centers (Yellowstone, WY, USA; Bandelier, NM, USA; Long Valley, CA, USA; Taupo, NZ; Mount St. Helens, WA, USA), where all systems contain initial CO2. In contrast to the minutes to hours of total ascent time extracted from embayment volatile profiles using constant-rate models, our two-stage model resolves slower initial ascent times that span several to >10 h. Final ascent rates are 1–2 orders of magnitude faster than the initial extracted rates, in agreement with theoretical conduit flow model predictions. Application of EDiTS to embayments from the May 18th, 1980 eruption of Mount St. Helens results in an initial stage of ascent on the order of hours consistent with the timing of magma arrival at the surface from the seismically-inferred storage region (7–9 km) ~3.5 h after the initial blast, and a final stage of ascent (<1–5 min) in close agreement with time-integrated bubble number densities. Our combined numerical, experimental, and natural results suggest that, with the application of more advanced models, the melt embayment can provide a complete picture of magma decompression timescales from the deep conduit to the surface.

40 Ar/ 39 Ar chronostratigraphy and geochemistry of Harrats Ishara and Kura in the Hamd–Jizl Basin, Madinah, Saudi Arabia

An extensive system of mafic dykes, volcanic centers, and provincial flood basalts (harrat) have erupted within the Arabian flank of the Red Sea rift since the early Oligocene. The age, spatial distribution, and composition of these features correlates with distinct changes in rift flank structural architecture as the margin progressed through successive deformational phases. A widespread pre-rift dyke swarm in the proto-Red Sea region c. 25 Ma created thermomechanical crustal conditions favourable for diffuse continental rifting and inboard extensional strain markers like the Hamd–Jizl Basin (HJB) originated amagmatically c. 23 Ma. Six million years later, basalt flows from Harrats Ishara and Kura entered the HJB at several points along strike during a 4.5 million year-long eruptive phase. Two distinct volcanic units developed in Harrat Ishara: a predominantly alkali-olivine basalt with an age range of 18.0–14.3 Ma (Ishara A) and a mostly olivine transitional basalt with an age range of 16.6–12.1 Ma (Ishara B). The contemporaneous eruption of these units indicates that multiple magma sources operated simultaneously during the middle Miocene. Both units generally decrease in alkalinity and rare earth element abundances up section, signaling either ongoing melt fractionation processes or output from different magma chambers with time.

Pneumatolytic type of ore mineralization within the elbrus volcanic center (North caucasus, Russia)

Pneumatolytic type of ore mineralization within the elbrus volcanic center (North caucasus, Russia)

The Cenozoic Abanico rift system: Implications of increased southward extension in the southern central Andes, in Chile

The late Eocene-Early Miocene intra-arc Abanico extensional basin represents a major feature in central Chile (∼31°-42°S). Such basin system concentrated the magmatic activity along this Andean region and hosted locally more than 3.000 m of volcanic and volcaniclastic deposits (Abanico Formation) over a southward thinning crust (<35 km thick). South of 34°30′S, major changes occur in and out of the basin realm that modify the rather regular distribution of the geological features observed northwards, namely: 1. Southward termination of the exposures of the Farellones Formation, 2. Increased presence of pre-Abanico volcanic and plutonic exposures in the basin realm, 3. Southward width increase of the basin, 4. Abrupt end in the Maule region of the Cretaceous plutonic swaths exposed in the eastern Coastal Cordillera and westward bend of the swaths of Mesozoic units, 5. Increased presence of primitive volcanic rocks in the west-side of the basin, and 6. Presence, south of 36°S, of the Cura-Mallín Formation on the east-side of the basin. These changes are caused by the oblique orientation of the Abanico basin relative to the Mesozoic structural and paleogeographic trends and are enhanced by the increased southward opening of the basin. This opening would have occurred through a westward scissors-like rotation of the block of continental crust located west of the basin, which corresponds to the present-day Coastal Cordillera south of 34°30′S. The eastern border, which is represented by the El Diablo fault, remained essentially fixed. The westward rotation is supported by the north-northeast orientation of the Mesozoic plutonic and stratified swaths south of 34°30′S. Additionally, the swaths of Mesozoic rocks are obliquely interrupted by the similarly oriented western border of the basin. Previous analogue models designed to understand the opening and closure mechanisms of an extensional basin show that greater extension and subsidence in the basin occur next to the mobile border, which in the Abanico basin was located on its west side. This observation aligns with the increased southward presence of primitive volcanic rocks, like the Colbún Formation, in the region where crustal thinning would have been greater. Analogue models show that the increased southward extension allowed the development of wider depocenters. In the Maule region, an eastern depocenter near the El Diablo fault, hosted the Early to Middle Miocene deposits of the Cura-Mallín Formation. During contraction, deeply rooted reverse faults can involve basement rocks of the basin and bring them up to the surface as observed in the Abanico basin realm in the Tinguiririca and Maule regions. This study also highlights the continuity and importance of the El Diablo fault in the configuration of the basin. It is shown that, notwithstanding the existence of major changes south of 34°30′S, this fault represents a major structural element that remained unchanged along the Andean segment comprised between ∼32° and 36°30′S separating the Cenozoic deposits to the west from Mesozoic deposits to the east. This fault can be traced continuously from, at least, the Aconcagua valley (33°S) to the Maule region (36°30′S), and most probably extends further south. Apart from separating the Cenozoic Abanico basin deposits from the Mesozoic units involved in the east-vergent Andean fold-thrust belt, its trace forms an almost straight line that connects the Palomo, Tinguiririca, Planchón-Peteroa, Descabezados, San Pedro and Rezago volcanic centers and complexes, hot-springs, and hydrothermal alterations.

Disentangling the Roles of Subducted Volatile Contributions and Mantle Source Heterogeneity in the Production of Magmas Beneath the Washington Cascades

AbstractThe compositional diversity of primitive arc basalts has long inspired questions regarding the drivers of magmatism in subduction zones, including the roles of decompression melting, mantle heterogeneity, and the amount and composition of slab‐derived materials. This contribution presents the volatile (H2O, Cl, and S), major, and trace element compositions of melt inclusions from basaltic magmas erupted at three volcanic centers in the Washington Cascades: Mount St. Helens (two basaltic tephras, 2.0–1.7 ka), Indian Heaven Volcanic Field (two &lt;600 ka basaltic hyaloclastite tuffs), and Glacier Peak (late Pleistocene to Holocene basaltic tephra from Whitechuck and Indian Pass cones). Compositions corrected to be in equilibrium with mantle olivine display variability in Nb and trace element ratios indicative of mantle source variability that impressively spans nearly the entire range of arc magmas globally. All volcanic centers have magmas with H2O and Cl contributions from the downgoing plate that overlap with other Cascade Arc segments. Volatile abundances and trace element ratios support a model of melting of a highly variably mantle wedge driven by a subduction component of variably saline fluid and/or slab partial melt. Magmas from Glacier Peak in northern Washington have unusually high Th/Yb ratios that are similar to Lassen region basalts, indicating possible contributions of “subcreted” metasediments that geophysical data suggest are not present beneath central Oregon and southern Washington. This data set adds to the growing inventory of primitive magma volatile concentrations and provides insight into spatial distributions of mantle heterogeneity and the role of slab components in the petrogenesis of arc magmas.

Magma storage conditions beneath a peralkaline caldera in the Main Ethiopian Rift

The numerous volcanic centres in the Main Ethiopian Rift (MER) present significant but poorly understood hazards to local populations. The MER is also an important site to gain insights into tectonic processes as it captures the transition from continental rifting (to the south) to incipient seafloor spreading (to the north). Peralkaline magmas account for around 90% of the volcanic products found in the MER. Determining the conditions under which these magmas evolve is critical to understanding rift-related volcanism and its associated hazards. Corbetti Caldera has an extensive record of large-scale, predominantly aphyric, peralkaline rhyolite eruptions. However, little is known about the mafic magmas from which these highly differentiated melts have evolved. Here we present data from the only basaltic deposit found within the caldera, coupled with whole rock, glass and mineral analysis of the peralkaline products, to investigate magma storage conditions at Corbetti. We demonstrate that magma mixing played a role in the evolution of the basaltic magmas and use RhyoliteMELTS modelling to show Corbetti's peralkaline magmas likely evolved at pressures between 100 and 250 MPa, from a magma with an initial water content of 0.5–1 wt%, at or below the QFM buffer. Mineral hygrometry on the sparse crystal populations corroborates the RhyoliteMELTS modelling, suggesting that the basaltic magma had 0.1–1.2 ± 0.32 wt% H2O, and the peralkaline magmas an average of ∼5.5 ± 1.25 wt% H2O. These results also match melt inclusion data for Corbetti and other peralkaline systems. We also provide new 40Ar/39Ar ages for two eruptions, a pre-caldera rhyolitic lava flow (206.7 ± 0.9 ka) and a post-caldera peralkaline ignimbrite (160 ± 0.8 ka). These results add to our understanding of the history of Corbetti and the storage conditions of peralkaline magmas within a continental rift setting and highlight the hydrous nature of Corbetti's magmas and the role that H2O plays during explosive eruptions.

Formation temperature of microcrystalline dolomite in oil-rich black shale: Implications for hydrocarbon accumulation in a volcanic-lacustrine rift basin

Lacustrine microcrystalline dolostones in the ancient records have been commonly interpreted as of an evaporitic origin, on the basis of their numerous occurrences in modern evaporative shallow lakes. This interpretation obscures the forming mechanisms of dolostones formed in a deep lacustrine setting, where dolostones interbedded with shale become a focus of hydrocarbon source rocks in recent years. In this study, an alternative non-evaporitic origin of dolomite is interpreted for a Permian lacustrine laminated dolostone. The dolostones are thinly intercalated with tuffaceous shales and composed mainly of volumetrically abundant, low-ordered, subhedral to euhedral, microcrystalline dolomites with a variable amount of pyroclastics. Evidence for evaporation (e.g., the presence of halite, gypsum and shrinkage cracks) and typical methanogenic microfossils is absent. Instead, our dolostones are unevenly and locally distributed as thin beds or lenses and intercalated with tuffaceous shale and have a formation temperature ranging from 31.4 to 73.8 °C, as estimated from dolomite clumped oxygen? isotope thermometry. In addition, the dolostones have a87Sr/86Sr ratio similar to that of the mantle (0.705154 on average) and extremely negative δ18O (−2.9‰ to −24‰) and positive δ13C (5.1‰–10.8‰) values. A lateral trend of increasing formation temperature of the dolostone toward the volcanic center correlates with the trends of: 1) an increasing amount of pyroclastics, euhedral dolomite crystals, Fe, and TOC, and an increasingly positive δEu value, and 2) decreasing values of 87Sr/86Sr and δ18O. The correlation suggests that a local temperature increase and frequent elemental (Fe and Mg) cycling between the conditions of dolomite supersaturation and undersaturation, which are controlled by episodic volcanic-hydrothermal activities are the causes for dolomite nucleation and growth, resulting in the formation of massive dolostone during deposition and shallow burial. The process also affected organic matter accumulation and local transformation to immature or early mature hydrocarbon during deposition and early burial. The results establish a volcanism-controlled temperature range for primary dolomite nucleation and growth, as well as hydrocarbon generation in a deep lake environment during the period of deposition and early burial, and provide insights on shale oil exploration in petroliferous rift basins.

Syn and post-glacial volcanic activity in the western retroarc area of the andean southern volcanic zone: Overo volcano

Located in a poorly studied zone of the rear arc in the Transitional Southern Volcanic Zone, the Overo volcano presents a complete record of syn-glacial, inter-glacial, and post-glacial activity, marked by significant changes in eruptive mechanisms over time. Following pre-glacial units, an effusive syn-glacial basaltic phase was succeeded by an extensive inter-glacial diking that fed andesitic monogenetic fields. Subsequently, an eruptive pulse, possibly coupled with the destabilization of the main edifice, resulted in a block and ash deposit channeled into a glacial valley towards the southwest, later covered by extensive basaltic lava flows. Then, different types of pyroclastic flows and glacial deposits covered these earlier units. Finally, these younger glacial deposits were overlain by a localized postglacial pyroclastic deposit and widespread lava flows, which ultimately reconstructed the deeply eroded glacial morphology of the Overo volcano. This volcanic record reveals the young and complex activity of this poorly known western rear-arc volcanic center in the Transitional Southern Volcanic Zone of the Andes, whose evolution is more akin to the arc front volcanoes to the west than to their eastern within-plate counterparts. We propose that this magmatic evolution was influenced by the Pleistocene-Holocene Last Glacial Maximum activity (LGM), beginning with an initial phase of volcanic activity encapsulated within glacial conditions. Glacial retreat then triggered isostatic uplift and crustal relaxation, resulting in dike formation, flank instability, and syn-eruptive collapse of the volcanic edifice. Subsequently, a lesser glacial pulse partially eroded the volcano morphology and produced the currently preserved glacial deposits. Finally, a latest and limited post-glacial pyroclastic eruption occurred simultaneously with the formation of an extensive monogenetic basaltic field, shaping the final morphology of the volcano.

Shallow storage of the explosive Earthquake Flat Pyroclastics magma body, Okataina Volcanic Center, Taupo Volcanic Zone, New Zealand: evidence from phase-equilibria experiments

Rhyolitic tuffs range widely in their crystal contents from nearly aphyric to crystal-rich, and their crystal cargoes inform concepts of upper crustal magma reservoirs. The Earthquake Flat pyroclastics (Okataina Volcanic Center, Taupo Volcanic Zone, New Zealand) are 10 km3 of rhyolitic tuffs with abundant (~ 40 vol.%) plagioclase and quartz, minor biotite, hornblende, and orthopyroxene, and accessory Fe-Ti oxides, apatite, and zircon, set in high-silica rhyolitic glass. Major minerals form large, euhedral phenocrysts and abundant glomerocrysts with few disequilibrium textures excepting some faintly resorbed quartz. Plagioclase phenocrysts have thick rims of nearly constant composition near An30, and hornblende is weakly zoned or unzoned. The abundant and texturally complex mineral assemblage contrasts with the nearby (~ 25 km), nearly synchronous, but more voluminous and crystal-moderate rhyolite tuffs from Rotoiti caldera. New H2O-saturated phase-equilibria results on the erupted Earthquake Flat melt (glass) determine its co-saturation with the partial phenocryst assemblage of plagioclase, quartz, biotite, and Fe-Ti oxides at: 140 MPa, 755 ºC. These closely approximate the conditions of the pre-eruptive magma body assuming it was saturated with nearly pure H2O and at an fO2 of ~ Ni–NiO. Absence of hornblende and orthopyroxene from the synthesized assemblages may result from those minerals being in a peritectic reaction relation with melt to produce biotite, so they would not grow from the liquid used as starting material. Experimental results on Rotoiti rhyolite (Nicholls et al. 1992) show that the two bodies resided at similar pressures, temperatures, and fO2s. Lower crystal abundance of the Rotoiti tuffs may result from slight compositional differences. We interpret that the Earthquake Flat pyroclastics were sourced from the crystal-rich periphery of a mushy reservoir system with the Rotoiti occupying a more melt-rich central location. Uncertain is whether this was a single intrusion zoned continuously in crystallinity, or discrete adjacent intrusions, but our results illustrate and quantify complexities of magma storage across relatively short distances.

A revised stratigraphic model for the ∼ 1910–1835 Ma Tanami Group, the Northern Territory, Australia: Implications for exploration targeting

The Paleoproterozoic Tanami Group of the Granites—Tanami Orogen hosts a known endowment of > 20 Moz of gold mineralization across numerous deposits. Throughout the region, host lithologies impart a first-order control on the style, grade and endowment of mineralization. However, the genetic relationship and regional distribution of three regionally correlative formations (the Dead Bullock, Mount Charles and Stubbins formations) remains enigmatic. Here, we combine lithostratigraphic analysis, lithogeochemistry, detrital zircon geochronology and geophysical observations to propose a new evolutionary model for the Tanami Group. We find that the Tanami Group was deposited in a continental back-arc basin setting and developed in four stages. Package A marked the onset of the rifting of the continental Archean basement. Sedimentary rocks within Package A are dominantly felsic in composition and display a unimodal detrital zircon age population at ∼ 2530 Ma. Package B marks the onset of mafic volcanism; consequently, volcaniclastic and sedimentary rocks within Package B display a distinct mafic geochemical signature. Detritus shed from an Archean basement persisted throughout this period, and a unimodal detrital age component at ∼ 2530 Ma persists. Package C marked the waning stages of voluminous mafic volcanism across the Tanami Basin. Sedimentary volcaniclastic rocks within Package C display a mixed mafic and felsic composition. Detrital zircon age spectra record multiple peaks at ∼ 2500, 3000 and 3350 Ma, which may reflect erosion of older Archean basement domains or changes in paleo-drainage. Package D marked a significant change in the source of detritus entering the basin from dominantly Archean to Paleoproterozoic, with dominant detrital age peaks at ∼ 1860 Ma. Despite significant lithofacies variations, the presented model highlights multiple correlations between the Dead Bullock and Mount Charles formations. The Mount Charles Formation likely represents a volcanic source-proximal, high-energy depocenter. In contrast, the Dead Bullock Formation was deposited in a lower energy setting, distal to volcanic centers. These observations are important to mineral exploration. Notably, the Mount Charles Formation is distributed over a significantly greater area than previously thought, and the under-represented Hangingwall and Footwall sequences of this formation are interpreted to host significant gold mineralization at the Oberon and Groundrush deposits.

Facies models for rocky shorelines and their application to transgressed basement highs in the North Sea

Rocky shorelines form where basement highs are eroded and flooded during marine transgressive events. Although the Mesozoic North Sea rift generated numerous platform margins and rotated fault blocks that acted as basement highs, rocky shoreline deposits have not previously been reported in this region. Rocky shoreline deposits are usually represented in the rock record by thin conglomerates overlying major unconformities. They are typically characterized by their ichnological aspects rather than by their depositional facies. This study used the sedimentological aspects of modern and Miocene rocky shorelines from Spain and Austria to create facies models, which were then applied to the recognition of rocky shorelines in the Mesozoic of the Central North Sea. Our results demonstrate that structureless, clast-supported, poorly to moderately sorted conglomerate–breccia deposits are associated with competent basement lithologies, which produce hard, resistant coastal cliffs around previously overlooked volcanic centres in the subsurface of the North Sea. The basement lithologies in most of the Central North Sea favoured the formation of softer coastal cliffs, with less resistant lithologies that did not generate or preserve gravel-sized particles. These softer coastal cliffs were mostly characterized by low-angle, unconformity-bounded sandstones and fine-grained deposits, precluding the preservation and recognition of Mesozoic rocky shores in much of the North Sea's stratigraphic record.

Formation and hydrothermal alteration of a volcanic center: Melt pooling and mass transfers at Langseth Ridge (Gakkel Ridge, Arctic Ocean)

Volcanic centers are characteristic features of ultraslow-spreading mid-ocean ridges, the least-explored parts of the global ridge system. Volcanic centers can provide insights into deep magmatic and metamorphic processes at these ridges. Here, we present data from the largest volcanic center on the Gakkel Ridge, the Langseth Ridge, situated at 60–62°E. Langseth is ∼10 km wide, consisting of three peaks that rise to 585 m water depth, some 3–4 km above the surrounding seafloor. It strikes perpendicular to Gakkel's spreading direction and can be traced for ∼40 km, which translates to an age of ∼8 Myr. Seafloor imaging revealed abundant (pillow) basalt but also fissures and geologic faults across the Langseth Ridge. Basaltic rocks were sampled at all summits and diabase at the slope of the northern summit that dips into the rift valley.Our samples are of normal to depleted mid-ocean ridge basalt composition and exhibit a wide range of major and trace element contents, due to magmatic processes, accumulation of macrocrysts, and hydrothermal alteration. Radiogenic isotope ratios, most notably 143Nd/144Nd and 208Pb/206Pb, trend from typical rift valley compositions to isotopically enriched values with increasing distance to the rift valley. This trend may imply melt pooling from different sources, potentially representing a shift from shallow melting beneath the rift valley to deeper melting of enriched sources and higher degrees of melting underneath Langseth. Mineral compositions and plagioclase sieve textures imply prolonged storage of magma at depth prior to eruption. Hydrothermal alteration occurred over a range of conditions. Basalt from the summits is weakly altered at temperatures ≪100 °C, which likely occurred in situ at the summit sites. Diabase samples experienced chloritization and albitization and display epidote and quartz veins, which formed at >300 °C. These assemblages and temperatures are typical for lower crustal levels and imply uplift of the samples of >1 km. Diabase samples from the Afanasenkov Seamount, another volcanic center on the Gakkel Ridge that we investigated for comparison, were altered under comparable conditions.Our findings suggest a combined volcanic–tectonic origin of the studied volcanic centers, potentially implying that such complexes may generally form due to the interplay of magmatism and tectonics. Researching volcanic centers has the potential to further our understanding of both deep and shallow crustal processes at ultraslow-spreading ridges, providing further insights into the role of these centers as linkages between lithosphere and hydrosphere and the (deep) biosphere they sustain.

Onverwacht Group, Barberton Greenstone Belt, South Africa: 300 m.y. development of a Paleoarchean stagnant lid

The 3.55−3.26 Ga Onverwacht Group, Barberton Greenstone Belt, South Africa and Eswatini, records nearly 300 m.y. of Paleoarchean history dominated by basaltic, komatiitic, and felsic volcanism. It provides the oldest well-preserved record of crustal development and the nature and evolution of associated surface environments, ocean, atmosphere, and biosphere on early Earth. Sedimentary layers within this 10- to 12-km-thick sequence represent a range of proximal to distal sedimentary environments relative to the centers of volcanism. Proximal deposits include coarse felsic breccias, conglomerates, and sandstones and mafic to komatiitic lapillistones that often show evidence for deposition in relatively shallow water. Distal deposits are composed of fine pyroclastic debris, chemical sediments, and biogenic materials deposited under subaqueous conditions during local volcanic quiescence. They show abundant current-produced features and are interpreted to have formed at water depths of a few hundred meters or less under the influence of tidal and/or ocean-circulation currents. Terrigenous clastic sediments formed by the weathering and erosion of older rocks are essentially absent. The Onverwacht Group was deposited under marine conditions on what appears to have been a water world with little evidence of large land areas and no evidence of active tectonism. It is interpreted to represent the upper part of a Paleoarchean stagnant lid overlying but decoupled from an active mantle. The Barberton Greenstone Belt and other Paleoarchean terranes offer a view of this lid over the past 300 m.y. of its development and during the early stages of its fragmentation and disruption.

Petrogenetic evolution of the Jalapa del Marqués pluton: Miocene arc magmatism in southern Mexico and its tectonic implications

Arc-related plutonic rocks in southern Mexico occur along the Pacific coast in the Sierra Madre del Sur (SMS) igneous province. These rocks record the migration of the Caribbean–North America–Farallon/Cocos junction towards the east between the Late Cretaceous and Miocene. However, little is known about the petrogenesis of igneous rocks that were emplaced during the latest expressions of magmatism in the SMS (i.e., early–middle Miocene). Based on fieldwork, petrography, whole-rock and mineral geochemistry, thermobarometric calculations, and zircon UPb dating, we document the petrogenetic evolution of the Jalapa del Marqués pluton and adjoining intrusions, which are the youngest so far dated in the SMS and are located ∼150 km inland from the present-day trench. This subduction-related magmatic complex consists of granites, granodiorites, diorites, dacite porphyries, and intermediate-to-felsic dikes, whose geochemical features suggest early evolutionary processes driven by fractionation of amphibole ± clinopyroxene (moderate–high pressures) or plagioclase (low pressures). Pulses of evolved magmas were episodically injected into upper-crustal reservoirs (7–14 km) over ∼9 Myr, starting at 22.14 ± 0.65 Ma (2σ)—earlier than previously thought (∼15.7 Ma)—and producing subsequent pulses at 19.45 ± 0.36 Ma, from 16.36 ± 0.37 to 15.8 ± 0.34 Ma, and from 13.9 ± 0.5 to 13.3 ± 0.2 Ma. These ages correlate with the onset and early stages of volcanic centers in the Trans-Mexican Volcanic Belt (TMVB), some of which were located up to 500 km landward from the present-day trench. The large difference in the distance-to-trench between coeval magmatism in the Jalapa del Marqués pluton and the TMVB, indicates a variation in the subduction angle between two segments of the Farallon/Cocos Plate that were likely separated by a slab-tear. The eastward migration of this slab-tear (along with the triple junction) terminated the magmatic activity in the SMS by switching the locus of magmatism to the TMVB. This probably occurred shortly after the studied magmatic complex was assembled.

Petrology and geochemistry of Mount Lurus lavas: evidence for alkaline to subalkaline transition in the rear arc of Java, Sunda Arc

Abstract Mount Lurus, a small Quaternary volcanic center in the rear of Java Island, Sunda arc, has erupted both alkaline and subalkaline lavas in a short span of active period. Furthermore, this small volcanic edifice, together with its neighboring alkaline volcanic center of Ringgit-Beser, corresponds to a shallower slab depth (∼160 km) than Muria in Central Java (270 km). We report the mineralogical and geochemical characteristics of Lurus lavas and compare them to the Ringgit-Beser lavas. All Lurus lavas show typical characteristics of island arcs with elevated Ba, Sr, Pb, and K and negative anomalies of Nb and Ti. The older Lurus edifice (Old Lurus) produces alkaline lavas with high alkalinity, which are trachyandesite and phonotephrite in composition. The younger Lurus (Young Lurus) generates calc-alkaline lavas with basaltic andesite composition. Both groups have plagioclase, clinopyroxene, olivine, and opaque minerals as phenocrysts phase. Leucite is found in one sample of Old Lurus with the lowest SiO2 content (48.92 wt%) and highest alkali content (5.38 wt%). Old Lurus lavas tend to show higher contents of incompatible elements (Pb, Ba, Nd, and Rb) than Young Lurus. The mineralogical and compositional differences between these lavas suggest that these magmas are sourced from at least three magma chambers of distinct composition: highly potassic magma chamber (comparable to the evolved-K series of Ringgit-Beser) and hydrous-potassic magma which fed Old Lurus edifice; and calc-alkaline magma chamber which fed Young Lurus edifice.

Crustal Characterization of the Hengill Geothermal Fields: Insights From Isotropic and Anisotropic Seismic Noise Imaging Using a 500‐Node Array

AbstractThe Hengill volcano and its associated geothermal fields represent Iceland's most productive harnessed high‐temperature geothermal fields, where resources are fueled by cooling magmatic intrusions connected to three volcanic systems. The crustal structure in this area is highly heterogeneous and shaped by the intricate interplay between tectonic forces and magmatic/hydrothermal activities. This complexity makes detailed subsurface characterization challenging. In this study, we aim to push the current resolution limits using a 500‐node temporary seismic array and perform an isotropic and, for the first time, radially‐anisotropic velocity model of the area. The high‐resolution isotropic velocity model reveals the characteristic N30ºE fissure swarm that crosses the area within the top 500 m and outlines a deep‐seated low‐velocity body composed of cooling magmatic intrusions at 5 km depth. This deeper body is located near the eastern part of the three volcanic centers and connected to a shallower body at 2–3 km depth that strikes westward toward Hengill volcano. Additionally, our study discovered that non‐induced earthquakes deeper than 2 km align with velocity contrasts that reflect structural variability, indicating the potential to identify deep permeable pathways using dense array imaging. The anisotropic model indicates that the shallow crust of Hengill within the top 2 km is dominated by vertical fractures or cracks, likely attributed to overall divergent deformation from rifting in the study area. This characteristic is diminished at depths greater than 2–3 km, replaced by a layering pattern where the lava flows and/or subhorizontal intrusions become the primary factors influencing the observed anisotropy.