Search

The conditions under which halogens partition in favor of an exsolved fluid relative to the coexisting melt are key for understanding many magmatic processes, including volcanic degassing, evolution of crustal melt bodies, and ore formation. We report new F, Cl, and Br fluid/melt partition coefficients for intermediate to silicic melts, for which F and Br data are particularly lacking; and for varying CO2-H2O contents to assess the effects of changing fluid composition (XH2O) on Br fluid/melt partitioning for the first time. The experiments were conducted at pressures 50–120 MPa, temperatures 800–1100 °C, and volatile compositions [molar XH2O = H2O/(H2O +CO2)] of 0.55 to 1, with redox conditions around the Nickel-Nickel Oxygen buffer (fO2 ≈ NNO). Experiments were not doped with Cl, Br, or F and were conducted on natural crystal-bearing volcanic products at conditions close to their respective pre-eruptive state. The experiments therefore provide realistic constraints on halogen partitioning at naturally occurring, brine-undersaturated conditions. Measurements of Br, Cl, and F were made by Secondary Ion Mass Spectrometry (SIMS) on 13 experimental glass products spanning andesite to rhyolitic compositions, together with their natural starting materials from Kelud volcano, Indonesia, and Quizapu volcano, Chile. Fluid compositions were constrained by mass balance. Average bulk halogen fluid/melt partition coefficients and standard deviations are: DClfluid/melt = 3.4 (±3.7 1 s.d.), DFfluid/melt = 1.7 (±1.7), and DBrfluid/melt = 7.1 (±6.4) for the Kelud starting material (bulk basaltic andesite), and DClfluid/melt = 11.1 (±3.5), DFfluid/melt = 0.8 (±0.8), and DBrfluid/melt = 31.3 (±20.9) for Quizapu starting material (bulk dacite). The large range in average partition coefficients is a product of changing XH2O, pressure and temperature. In agreement with studies on synthetic melts, our data show an exponential increase of halogen Dfluid/melt with increasing ionic radius, with partitioning behavior controlled by melt composition according to the nature of the complexes forming in the melt (e.g., SiF4, NaCl, KBr). The fundamental chemistry of the different halogens (differing ionic size and electronegativities) controls the way in which partitioning responds to changes in melt composition and other variables. Experimental results confirm that more Cl partitions into the fluid at higher bulk Cl contents, higher melt Na, higher fluid XH2O ratios, and lower temperatures. Bromine shows similar behavior, though it seems to be more sensitive to temperature and less sensitive to Na content and XH2O. In contrast, F partitioning into the fluid increases as the melt silica content decreases (from 72 to 56 wt% SiO2), which we attribute to the lower abundance of Si available to form F complexes in the melt. These new data provide more insights into the conditions and processes that control halogen degassing from magmas and may help to inform the collection and interpretation of melt inclusions and volcano gas data.

There is a temporal correlation between the peak activity of the North Atlantic Igneous Province (NAIP) and the Paleocene–Eocene Thermal Maximum (PETM), suggesting that the NAIP may have initiated and/or prolonged this extreme warming event. However, corroborating a causal relationship is hampered by a scarcity of expanded sedimentary records that contain both climatic and volcanic proxies. One locality hosting such a record is Fur Island in Denmark, where an expanded pre- to post-PETM succession containing hundreds of NAIP ash layers is exceptionally well preserved. We compiled a range of environmental proxies, including mercury (Hg) anomalies, paleotemperature proxies, and lithium (Li) and osmium (Os) isotopes, to trace NAIP activity, hydrological changes, weathering, and seawater connectivity across this interval. Volcanic proxies suggest that NAIP activity was elevated before the PETM and appears to have peaked during the body of the δ13C excursion, but decreased considerably during the PETM recovery. This suggests that the acme in NAIP activity, dominated by flood basalt volcanism and thermogenic degassing from contact metamorphism, was likely confined to just ~200 kyr (ca. 56.0–55.8 Ma). The hundreds of thick basaltic ashes in the post-PETM strata likely represent a change from effusive to explosive activity, rather than an increase in NAIP activity. Detrital δ7Li values and clay abundances suggest that volcanic ash production increased basaltic reactive surface area, likely enhancing silicate weathering and atmospheric carbon sequestration in the early Eocene. Signals in lipid biomarkers and Os isotopes, traditionally used to trace paleotemperature and weathering changes, are used here to track seaway connectivity. These proxies indicate that the North Sea was rapidly cut off from the North Atlantic in under 12 kyr during the PETM recovery due to NAIP thermal uplift. Our findings reinforce the hypothesis that the emplacement of the NAIP had a profound and complex impact on Paleocene–Eocene climate, both directly through volcanic and thermogenic degassing, and indirectly by driving regional uplift and changing seaway connectivity.

Marine tetrapods occupy important roles in modern marine ecosystems and often gather in large aggregations driven by patchy prey distribution, 1 ,2 social or reproductive behaviors, 3 ,4 or oceanographic factors. 5 Here, we show that similar grouping behaviors evolved in an early marine tetrapod lineage, documented by dozens of specimens of the giant ichthyosaur Shonisaurus in the Luning Formation in West Union Canyon, Nevada, USA. 6 ,7 A concentration of at least seven skeletons closely preserved on a single bedding plane received the bulk of previous attention. However, many more specimens are preserved across ∼106 square meters and ∼200 stratigraphic meters of outcrop representing an estimated >105–6 years. Unlike other marine-tetrapod-rich deposits, this assemblage is essentially monotaxic; other vertebrate fossils are exceptionally scarce. Large individuals are disproportionately abundant, with the exception of multiple neonatal or embryonic specimens, indicating an unusual demographic composition apparently lacking intermediate-sized juveniles or subadults. Combined with geological evidence, our data suggest that dense aggregations of Shonisaurus inhabited this moderately deep, low-diversity, tropical marine environment for millennia during the latest Carnian Stage of the Late Triassic Period (237–227 Ma). Thus, philopatric grouping behavior in marine tetrapods, potentially linked to reproductive activity, has an antiquity of at least 230 million years.

Satellite instruments play a valuable role in detecting, monitoring and characterising emissions of ash and gas into the atmosphere during volcanic eruptions. Plumes of ash and sulfur dioxide (SO2) from the April 2021 eruption of La Soufrière volcano on St Vincent in the Eastern Caribbean were observed by a multiple satellite instruments. This study looks at these plumes with two satellite instruments: the Advanced Baseline Imager (ABI) on the Geostationary Operational Environmental Satellite (GOES), and the Infrared Atmospheric Sounding Interferometer (IASI) on the MetOp platforms. Using true and false colour images, and brightness temperature difference images produced from the ABI data, a minimum of 32 eruptive events were identified. The ABI images were used to determine the approximate start and end times and character of each event. In this way the eruption has been divided into four phases: (1) an initial explosive event, (2) a sustained event lasting over nine hours, (3) a pulsatory phase with 23 explosive events in a 54 hour period and (4) a waning sequence of explosive events. The IASI instrument was used to study the dispersion of SO2 from this eruption. The results showed a highly complex structure to the plume, in terms of the column amounts and height, which is likely linked to the multiple explosive events. The SO2 is shown to have largely been emitted between 13 and 19 km. This was primarily in the upper troposphere and around the height of the tropopause, but with some emission into the stratosphere. The SO2 was transported around the globe with parts of the plume reaching as far as 45° S and 45° N. The largest SO2 atmospheric burden measured with IASI was 0.31±0.09 Tg, recorded on the 13 April 2021 (descending orbits). The SO2 masses were converted into fluxes. The SO2 flux was shown to peak on 10 April and then shown to decrease over time. By summing the IASI SO2 flux results, it is estimated that a total of 0.57±0.44 Tg of SO2 was emitted to the atmosphere. However, due to the limitations associated with the retrieval this should be considered a minimum estimate of the total mass of SO2 emitted. An average e-folding time of 7.09±5.70 days was computed based on the IASI SO2 results: similar to other tropical eruptions of this magnitude. There are a number of similarities between the 1979 and 2021 eruptions at La Soufrière. For example, both eruptions consisted of a series of explosive events with varied heights including some emission into the stratosphere. The similarities between the 1979 and 2021 highlight the importance of studying these eruptions to be prepared for future activity.

Past large igneous province (LIP) emplacement is commonly associated with mantle plume upwelling and led to major carbon emissions. One of Earth’s largest past environmental perturbations, the Toarcian oceanic anoxic event (T-OAE; ~183 Ma), has been linked to Karoo-Ferrar LIP emplacement. However, the role of mantle plumes in controlling the onset and timing of LIP magmatism is poorly understood. Using global plate reconstruction models and Lower Toarcian sedimentary mercury (Hg) concentrations, we demonstrate (i) that the T-OAE occurred coevally with Karoo-Ferrar emplacement and (ii) that timing and duration of LIP emplacement was governed by reduced Pangean plate motion, associated with a reversal in plate movement direction. This new model mechanistically links Earth’s interior and surficial processes, and the mechanism is consistent with the timing of several of the largest LIP volcanic events throughout Earth history and, thus, the timing of many of Earth’s past global climate change and mass extinction events.

The latest Permian mass extinction (LPME) was triggered by magmatism of the Siberian Traps Large Igneous Province (STLIP), which left an extensive record of sedimentary Hg anomalies at Northern Hemisphere and tropical sites. Here, we present Hg records from terrestrial sites in southern Pangea, nearly antipodal to contemporaneous STLIP activity, providing insights into the global distribution of volcanogenic Hg during this event and its environmental processing. These profiles (two from Karoo Basin, South Africa; two from Sydney Basin, Australia) exhibit significant Hg enrichments within the uppermost Permian extinction interval as well as positive Δ

The caldera-forming Samalas eruption of the Samalas–Rinjani volcanic complex on Lombok, Indonesia, in 1257 CE ranks as one of the most explosive and sulphur-rich eruptions of the past thousand years. Along with other significant events (e.g. Tambora 1815, El Chichón 1982), Samalas 1257 forms a class of powerful eruptions of magmas of intermediate alkaline composition, which remain relatively understudied. This study aims to better understand the pre-eruptive magma reservoir conditions and volatile storage capacities of intermediate alkaline systems, using the 1257 Samalas eruption as a case study. Systematic partial-equilibrium experiments were run to constrain potential magma reservoir conditions (25–200 MPa and 850–1000°C) by examining changes in the chemistry and textures of experimental glass and mineral products. Natural trachydacite pumice from the 1257 Samalas eruption was used as the experimental charge for a series of water-saturated experiments run in cold-seal pressure vessels, buffered at an oxygen fugacity ƒO2 of NNO/NNO + 1 log units. Assessing the match between natural and experimental glass compositions and mineral phases and textures reveals that pre-eruptive magma reservoir conditions were between 875°C and 930°C and ca. 100–150 MPa (4.5 ± 1km depth). Breakdown of experimental amphibole at pressures below 75 MPa, and plagioclase instability at 950°C/100 MPa, provide the strongest constraints on phase stabilities that are consistent with the magma storage region. In the observed temperature and pressure range, the natural compositions of plagioclase, amphibole and orthopyroxene are replicated in experimentally precipitated rims. Clinopyroxene and high anorthite plagioclase (An>64) in the natural mineral assemblage could not be replicated in the experiments, implying that these compositions most likely did not grow in the shallow pre-eruptive reservoir of Mt Samalas. The pre-eruptive storage depths of trachydacite magma are significantly shallower than the postulated magma storage at the currently active basaltic-andesite magmatic system at Rinjani, and thus monitoring magma storage depth may be a useful indicator of changing magma composition in the future. Our findings show that highly explosive VEI 7 eruptions of intermediate alkaline magma can be fed from a relatively limited range of storage pressures (100–150 MPa), suggesting accumulation of magma in one place, rather than the rapid extraction of magma from a vertically extensive, transcrustal magma system.

This study was supported by the Leverhulme Trust grant 2016–21 (Nature and impacts of Middle Pleistocene volcanism in the Ethiopian Rift). KF was supported by the UK's Natural Environment Research Council (NERC) grant NE/L013932/1 (RiftVolc: The Past, Present and Future of Rift Volcanism in the Main Ethiopian Rift), a Boise Fund grant from the Department of Zoology, University of Oxford, and acknowledges Fonds de Recherche Scientifique – FNRS MIS grant F.4515.20. Tephra work on the Chew Bahir cores in the Cambridge Tephra Lab by AA, AP and CL was made possible by NERC grant NE/K014560/1. Ar–Ar dating was supported by grants NIGFSC IP-1683-1116 and IP-1680-1116. The iCRAG lab is supported by SFI 13/RC/2092. The Main Ethiopian Rift (MER) is renowned as a focus of investigations into human origins. It is also the site of many large volcanic calderas, whose eruptions have spanned the timeframe of speciation, cultural innovation, and dispersal of our species. Yet, despite their significance for dating human fossils and cultural materials, the timing and geochemical signatures of some of the largest eruptions have remained poorly constrained at best. Here, through a programme of field surveys, geochemical analysis and 40Ar/39Ar dating, we report the ages of MER ignimbrites and link them to widespread tephra layers found in sequences of archaeological and paleoenvironmental significance. We date major eruptions of Fentale (76 ± 18 ka), Shala (ca. 145–155 ka), Kone (184 ± 42 ka and ca. 200 ± 12 ka) and Gedemsa (251 ± 47 ka) volcanoes, and correlate a suite of regionally important tephra horizons. Geochemical analysis highlights the predominantly peralkaline rhyolitic melt compositions (7.5–12 wt% Na2O + K2O, 70–76 wt% SiO2) across the central MER and remarkable similarity in incompatible trace element ratios, limiting the correlation of deposits via glass composition alone. However, by integrating stratigraphic and geochronological evidence from proximal deposits, lake sediment cores and distal outcrops at archaeological sites, we have traced ash layers associated with the ca. 177 ka Corbetti, ca. 145–155 ka Shala and ca. 108 ka Bora-Baricha-Tullu-Moye eruptions across southern Ethiopia. In addition to strengthening the tephrochronological framework that supports paleoenvironmental and archaeological work in the region, our findings have wider implications for evaluating the hypothesis of a middle Pleistocene ‘ignimbrite flare-up’ in the MER, and for evaluating the impacts of these great eruptions on landscapes, hydrology, and human ecology. Publisher PDF

Gas and particulate matter (PM) emissions from Masaya volcano, Nicaragua, cause substantial regional volcanic air pollution (VAP). We evaluate the suitability of low-cost SO2 and PM sensors for a continuous air-quality network. The network was deployed for six months in five populated areas (4–16 km from crater). The SO2 sensors failed and recorded erroneous values on multiple occasions, likely due to corrosion, requiring significant maintenance commitment. The PM sensors were found to be robust but data required correction for humidity. SO2 measurements could not be used as stand-alone tools to detect occurrence of VAP episodes (VAPE), but an SO2/PM correlation reliably achieved this at near-field stations, as confirmed by meteorological forecasts and satellite imagery. Above-background PM concentrations reliably identified VAPE at both near-field and far-field stations. We suggest that a continuous network can be built from a combination of low-cost PM and SO2 sensors with a greater number of PM-only sensors.

Understanding the timing of critical changes in volcanic systems, such as the beginning and end of eruptive behaviour, is a key goal of volcanic monitoring. Traditional approaches to forecasting these changes have used models motivated by the underlying physics of eruption onset, which assume that geophysical precursors will consistently display similar patterns prior to transition in volcanic state. We present a machine learning classification approach for detecting significant changes in patterns of volcanic activity, potentially signalling transitions during the onset or end of volcanic activity, which does not require a model of the physical processes underlying critical changes. We apply novelty detection, where models are trained only on data prior to eruption, to the precursory unrest at Augustine Volcano, Alaska in 2005. This approach looks promising for geophysically-monitored volcanic systems which have been in repose for some time, as no eruptive data is required for model training. We compare novelty detection results with multi32 class classification, where models are trained on examples of both non-eruptive and eruptive data. We contextualise the results of these classification models using constraints from petrological, satellite and visual observations from the 2006 eruption of Augustine Volcano. The transition from non-eruptive to eruptive behaviour we identify in mid-November 2005 is in agreement with previous estimates of the initiation of dike intrusion prior to the 2006 eruption. We find that models which include multiple types of data (seismic, deformation and gas emissions) can better distinguish between non-eruptive and eruptive data than models formulated on single data types.