Turmoil at Turrialba Volcano (Costa Rica): Degassing and eruptive processes inferred from high-frequency gas monitoring

3He/4He Signature of Magmatic Fluids from Telica (Nicaragua) and Baru (Panama) Volcanoes, Central American Volcanic Arc

Constraining the magmatic 3He/4He signature of fluids degassed from a magmatic system is crucial for making inferences on its mantle source. This is especially important in arc volcanism, where variations in the composition of the wedge potentially induced by slab sediment fluids must be distinguished from the effects of magma differentiation, degassing, and crustal contamination. The study of fluid inclusions (FIs) trapped in minerals of volcanic rocks is becoming an increasingly used methodology in geochemical studies that integrates the classical study of volcanic and geothermal fluids. Here, we report on the first noble gas (He, Ne, Ar) concentrations and isotopic ratios of FI in olivine (Ol) and pyroxene (Px) crystals separated from eruptive products of the Telica and Baru volcanoes, belonging to the Nicaraguan and Panamanian arc-segments of Central America Volcanic arc (CAVA). FIs from Telica yield air corrected 3He/4He (Rc/Ra) of 7.2–7.4 Ra in Ol and 6.1–7.3 in Px, while those from Baru give 7.1–8.0 Ra in Ol and 4.2–5.8 Ra in Px. After a data quality check and a comparison with previous 3He/4He measurements carried out on the same volcanoes and along CAVA, we constrained a magmatic Rc/Ra signature of 7.5 Ra for Telica and of 8.0 Ra for Baru, both within the MORB range (8 ± 1 Ra). These 3He/4He differences also reflect variations in the respective arc-segments, which cannot be explained by radiogenic 4He addition due to variable crust thickness, as the mantle beneath Nicaragua and Panama is at about 35 and 30 km, respectively. We instead highlight that the lowest 3He/4He signature observed in the Nicaraguan arc segment reflects a contamination of the underlying wedge by slab sediment fluids. Rc/Ra values up to 9.0 Ra are found at Pacaya volcano in Guatemala, where the crust is 45 km thick, while a 3He/4He signature of about 8.0 Ra was measured at Turrialba volcano in Costa Rica, which is similar to that of Baru, and reflects possible influence of slab melting, triggered by a change in subduction conditions and the contemporary subduction of the Galapagos hot-spot track below southern Costa Rica and western Panama.

Volcanic CO2 tracks the incubation period of basaltic paroxysms

The ordinarily benign activity of basaltic volcanoes is periodically interrupted by violent paroxysmal explosions ranging in size from Hawaiian to Plinian in the most extreme examples. These paroxysms often occur suddenly and with limited or no precursors, leaving their causal mechanisms still incompletely understood. Two such events took place in summer 2019 at Stromboli, a volcano otherwise known for its persistent mild open-vent activity, resulting in one fatality and damage to infrastructure. Here, we use a post hoc analysis and reinterpretation of volcanic gas compositions and fluxes acquired at Stromboli to show that the two paroxysms were preceded by detectable escalations in volcanic plume CO₂ degassing weeks to months beforehand. Our results demonstrate that volcanic gas CO₂ is a key driver of explosions and that the preparatory periods ahead of explosions in basaltic systems can be captured by precursory CO₂ leakage from deeply stored mafic magma.

Seven potential sources of arsenic pollution in Latin America and their environmental and health impacts

This review presents a holistic overview of the occurrence, mobilization, and pathways of arsenic (As) from predominantly geogenic sources into different near-surface environmental compartments, together with the respective reported or potential impacts on human health in Latin America. The main sources and pathways of As pollution in this region include: (i) volcanism and geothermalism: (a) volcanic rocks, fluids (e.g., gases) and ash, including large-scale transport of the latter through different mechanisms, (b) geothermal fluids and their exploitation; (ii) natural lixiviation and accelerated mobilization from (mostly sulfidic) metal ore deposits by mining and related activities; (iii) coal deposits and their exploitation; (iv) hydrocarbon reservoirs and co-produced water during exploitation; (v) solute and sediment transport through rivers to the sea; (vi) atmospheric As (dust and aerosol); and (vii) As exposure through geophagy and involuntary ingestion. The two most important and well-recognized sources and mechanisms for As release into the Latin American population's environments are: (i) volcanism and geothermalism, and (ii) strongly accelerated As release from geogenic sources by mining and related activities. Several new analyses from As-endemic areas of Latin America emphasize that As-related mortality and morbidity continue to rise even after decadal efforts towards lowering As exposure. Several public health regulatory institutions have classified As and its compounds as carcinogenic chemicals, as As uptake can affect several organ systems, viz. dermal, gastrointestinal, peptic, neurological, respiratory, reproductive, following exposure. Accordingly, ingesting large amounts of As can damage the stomach, kidneys, liver, heart, and nervous system; and, in severe cases, may cause death. Moreover, breathing air with high As levels can cause lung damage, shortness of breath, chest pain, and cough. Further, As compounds, being corrosive, can also cause skin lesions or damage eyes, and long-term exposure to As can lead to cancer development in several organs.

BVLOS UAS Operations in Highly-Turbulent Volcanic Plumes

Long-range, high-altitude Unoccupied Aerial System (UAS) operations now enable in-situ measurements of volcanic gas chemistry at globally-significant active volcanoes. However, the extreme environments encountered within volcanic plumes present significant challenges for both air frame development and in-flight control. As part of a multi-disciplinary field deployment in May 2019, we flew fixed wing UAS Beyond Visual Line of Sight (BVLOS) over Manam volcano, Papua New Guinea, to measure real-time gas concentrations within the volcanic plume. By integrating aerial gas measurements with ground- and satellite-based sensors, our aim was to collect data that would constrain the emission rate of environmentally-important volcanic gases, such as carbon dioxide, whilst providing critical insight into the state of the subsurface volcanic system. Here, we present a detailed analysis of three BVLOS flights into the plume of Manam volcano and discuss the challenges involved in operating in highly turbulent volcanic plumes. Specifically, we report a detailed description of the system, including ground and air components, and flight plans. We present logged flight data for two successful flights to evaluate the aircraft performance under the atmospheric conditions experienced during plume traverses. Further, by reconstructing the sequence of events that led to the failure of the third flight, we identify a number of lessons learned and propose appropriate recommendations to reduce risk in future flight operations.

Aerial strategies advance volcanic gas measurements at inaccessible, strongly degassing volcanoes

Volcanic emissions are a critical pathway in Earth's carbon cycle. Here, we show that aerial measurements of volcanic gases using unoccupied aerial systems (UAS) transform our ability to measure and monitor plumes remotely and to constrain global volatile fluxes from volcanoes. Combining multi-scale measurements from ground-based remote sensing, long-range aerial sampling, and satellites, we present comprehensive gas fluxes-3760 ± [600, 310] tons day^-1 CO₂ and 5150 ± [730, 340] tons day^-1 SO₂-for a strong yet previously uncharacterized volcanic emitter: Manam, Papua New Guinea. The CO₂/S_T ratio of 1.07 ± 0.06 suggests a modest slab sediment contribution to the sub-arc mantle. We find that aerial strategies reduce uncertainties associated with ground-based remote sensing of SO₂ flux and enable near-real-time measurements of plume chemistry and carbon isotope composition. Our data emphasize the need to account for time averaging of temporal variability in volcanic gas emissions in global flux estimates.

Application of hydrogeochemistry and isotopic characterization for the assessment of recharge in a volcanic aquifer in the eastern region of central Costa Rica

In the eastern region of central Costa Rica, land use in the sub-basins of the Maravilla-Chiz and Quebrada Honda rivers (47 km²) is dominated by agricultural and livestock production, while groundwater resources constitute the main drinking water supply. This study aimed to (a) evaluate the location of groundwater recharge areas and groundwater flow paths, and (b) provide a characterization of the hydrochemistry and possible anthropic impacts. Groundwater was collected from 20 sites during the dry and rainy seasons and analysed for major ions, water stable isotopes and ²²²Rn. Approximated recharge areas were estimated through a local altitudinal line based on isotopic compositions in springs. The hydrochemical and isotopic characterization of groundwater showed that the main recharge areas occur in the upper part of the basin, except for springs in the middle part of the basin probably due to a certain hydraulic disconnection from the upper part that facilitates local recharge processes. In the lower basin, groundwater exhibited greater transit times and longer flow paths. Low nitrate, chloride and sulphate concentrations found in groundwater indicate low leaching of fertilizers or urban wastewaters. Our results are focused to improve water resources and agricultural management plans in a dynamic tropical landscape.

Combined ground and aerial measurements resolve vent-specific gas fluxes from a multi-vent volcano

Volcanoes with multiple summit vents present a methodological challenge for determining vent-specific gas emissions. Here, using a novel approach combining multiple ultraviolet cameras with synchronous aerial measurements, we calculate vent-specific gas compositions and fluxes for Stromboli volcano. Emissions from vent areas are spatially heterogeneous in composition and emission rate, with the central vent area dominating passive emissions, despite exhibiting the least explosive behaviour. Vents exhibiting Strombolian explosions emit low to negligible passive fluxes and are CO₂-dominated, even during passive degassing. We propose a model for the conduit system based on contrasting rheological properties between vent areas. Our methodology has advantages for resolving contrasting outgassing dynamics given that measured bulk plume compositions are often intermediate between those of the distinct vent areas. We therefore emphasise the need for a vent-specific approach at multi-vent volcanoes and suggest that our approach could provide a transformative advance in volcano monitoring applications.

Combining multiple ultraviolet cameras with synchronous aerial measurements, the authors here present vent-specific gas compositions and fluxes for Stromboli volcano. The results show that gas compositions vary between different vents, mirroring differences in eruptive behavior.

Bubbles and Dust: Experimental Results of Dissolution Rates of Metal Salts and Glasses From Volcanic Ash Deposits in Terms of Surface Area, Chemistry, and Human Health Impacts

Explosive volcanic eruptions lead to ash deposition and subsequent leaching of contaminants into soils or surface water, impacting flora and fauna, including human health. This study determined the control of ash surface area and chemical composition on ash dissolution rates. Fresh, unhydrated ash samples from four contrasting volcanoes were analyzed in the laboratory. Column leachate tests were used to compare leaching rates over a range of basaltic to andesitic ashes as a function of time and surface area, to analyze the effects of ash deposition. It was found that surface area, measured both geometrically and by multipoint Brunauer-Emmett-Teller analysis, generally increases for a short time, gradually decreases, then increases over the rest of the leaching experiment, due to area to mass ratio fluctuations. After the column leachate tests, postleaching water analyses for elemental compositions were conducted by inductively coupled plasma-mass spectrometry and ion chromatography. Steady state dissolution rates initially decayed rapidly due to the smallest size fraction of ash (dust), which provides a large area of fresh leachable surfaces as well as the rapid dissolution of highly soluble metal salts. Some of the dissolved concentrations of elements relevant to human and ecosystem health such as F, Cd, Se, As, and Cr rose above World Health Organization (WHO) drinking water standards within an hour of experimental leaching. In nature, however, safe consumption standards are further dependent upon bioaccumulation and chronic exposure. As such, individual and recurring ash deposition events have applications to emergency response and preparedness in volcanic regions.

Subducting carbon

A hidden carbon cycle exists inside Earth. Every year, megatons of carbon disappear into subduction zones, affecting atmospheric carbon dioxide and oxygen over Earth's history. Here we discuss the processes that move carbon towards subduction zones and transform it into fluids, magmas, volcanic gases and diamonds. The carbon dioxide emitted from arc volcanoes is largely recycled from subducted microfossils, organic remains and carbonate precipitates. The type of carbon input and the efficiency with which carbon is remobilized in the subduction zone vary greatly around the globe, with every convergent margin providing a natural laboratory for tracing subducting carbon.

Forearc carbon sink reduces long-term volatile recycling into the mantle

Carbon and other volatiles in the form of gases, fluids or mineral phases are transported from Earth's surface into the mantle at convergent margins, where the oceanic crust subducts beneath the continental crust. The efficiency of this transfer has profound implications for the nature and scale of geochemical heterogeneities in Earth's deep mantle and shallow crustal reservoirs, as well as Earth's oxidation state. However, the proportions of volatiles released from the forearc and backarc are not well constrained compared to fluxes from the volcanic arc front. Here we use helium and carbon isotope data from deeply sourced springs along two cross-arc transects to show that about 91 per cent of carbon released from the slab and mantle beneath the Costa Rican forearc is sequestered within the crust by calcite deposition. Around an additional three per cent is incorporated into the biomass through microbial chemolithoautotrophy, whereby microbes assimilate inorganic carbon into biomass. We estimate that between 1.2 × 10⁸ and 1.3 × 10¹⁰ moles of carbon dioxide per year are released from the slab beneath the forearc, and thus up to about 19 per cent less carbon is being transferred into Earth's deep mantle than previously estimated.

CO2 flux emissions from the Earth's most actively degassing volcanoes, 2005-2015

The global carbon dioxide (CO₂) flux from subaerial volcanoes remains poorly quantified, limiting our understanding of the deep carbon cycle during geologic time and in modern Earth. Past attempts to extrapolate the global volcanic CO₂ flux have been biased by observations being available for a relatively small number of accessible volcanoes. Here, we propose that the strong, but yet unmeasured, CO₂ emissions from several remote degassing volcanoes worldwide can be predicted using regional/global relationships between the CO₂/S_T ratio of volcanic gases and whole-rock trace element compositions (e.g., Ba/La). From these globally linked gas/rock compositions, we predict the CO₂/S_T gas ratio of 34 top-degassing remote volcanoes with no available gas measurements. By scaling to volcanic SO₂ fluxes from a global catalogue, we estimate a cumulative “unmeasured” CO₂ output of 11.4 ± 1.1 Mt/yr (or 0.26 ± 0.02·10¹² mol/yr). In combination with the measured CO₂ output of 27.4 ± 3.6 Mt/yr (or 0.62 ± 0.08·10¹² mol/yr), our results constrain the time-averaged (2005–2015) cumulative CO₂ flux from the Earth’s 91 most actively degassing subaerial volcanoes at 38.7 ± 2.9 Mt/yr (or 0.88 ± 0.06·10¹² mol/yr).

Controls on explosive-effusive volcanic eruption styles

One of the biggest challenges in volcanic hazard assessment is to understand how and why eruptive style changes within the same eruptive period or even from one eruption to the next at a given volcano. This review evaluates the competing processes that lead to explosive and effusive eruptions of silicic magmas. Eruptive style depends on a set of feedback involving interrelated magmatic properties and processes. Foremost of these are magma viscosity, gas loss and external properties such as conduit geometry. Ultimately, these parameters control the speed at which magmas ascend, decompress and outgas en route to the surface, and thus determine eruptive style and evolution.

Eruptive styles at a single volcano may transition from explosive to effusive behaviour (or vice versa) at any given time. This review examines the underlying controls on eruptive styles such as magma viscosity, degassing and conduit geometry at volcanoes with silicic compositions.

Understanding and forecasting phreatic eruptions driven by magmatic degassing

This paper examines phreatic eruptions which are driven by inputs of magma and magmatic gas. We synthesize data from several significant phreatic systems, including two in Costa Rica (Turrialba and Poás) which are currently highly active and hazardous. We define two endmember types of phreatic eruptions, the first (type 1) in which a deeper hydrothermal system fed by magmatic gases is sealed and produces overpressure sufficient to drive explosive eruptions, and the second (type 2) where magmatic gases are supplied via open-vent degassing to a near-surface hydrothermal system, vaporizing liquid water which drives the phreatic eruptions. The surficial source of type 2 eruptions is characteristic, while the source depth of type 1 eruptions is commonly greater. Hence, type 1 eruptions tend to be more energetic than type 2 eruptions. The first type of eruption we term "phreato-vulcanian", and the second we term "phreato-surtseyan". Some systems (e.g., Ruapehu, Poás) can produce both type 1 and type 2 eruptions, and all systems can undergo sealing at various timescales. We examine a number of precursory signals which appear to be important in understanding and forecasting phreatic eruptions; these include very long period events, banded tremor, and gas ratios, in particular H₂S/SO₂ and CO₂/SO₂. We propose that if these datasets are carefully integrated during a monitoring program, it may be possible to accurately forecast phreatic eruptions.

Spaceborne detection of localized carbon dioxide sources

Spaceborne measurements by NASA's Orbiting Carbon Observatory-2 (OCO-2) at the kilometer scale reveal distinct structures of atmospheric carbon dioxide (CO₂) caused by known anthropogenic and natural point sources. OCO-2 transects across the Los Angeles megacity (USA) show that anthropogenic CO₂ enhancements peak over the urban core and decrease through suburban areas to rural background values more than ~100 kilometers away, varying seasonally from ~4.4 to 6.1 parts per million. A transect passing directly downwind of the persistent isolated natural CO₂ plume from Yasur volcano (Vanuatu) shows a narrow filament of enhanced CO₂ values (~3.4 parts per million), consistent with a CO₂ point source emitting 41.6 kilotons per day. These examples highlight the potential of the OCO-2 sensor, with its unprecedented resolution and sensitivity, to detect localized natural and anthropogenic CO₂ sources.