A CO2-rich gas trigger of explosive paroxysms at Stromboli basaltic volcano,Italy

In addition to rhythmic slug-driven Strombolian activity, Stromboli volcano occasionally produces discrete explosive paroxysms (2 per year on average for the most frequent ones) that constitute a major hazard and whose origin remains poorly elucidated. Partial extrusion of the volatile-rich feeding basalt as aphyric pumice during these events has led to consider their triggering by the fast ascent of primitive magma blobs from possibly great depth. Here I examine and discuss the alternative hypothesis that most of the paroxysms could be triggered and driven by the fast upraise of CO₂-rich gas pockets generated by bubble foam growth and collapse in the sub-volcano plumbing system. Data for the SO₂ and CO₂ crater plume emissions are used to show that Stromboli's feeding magma may originally contain as much as 2 wt.% of carbon dioxide and early coexists with an abundant CO₂-rich gas phase with high CO₂/SO₂ molar ratio (≥ 60 at 10 km depth below the vents, compared to ～ 7 in time-averaged crater emissions). Pressure-related modelling indicates that the time-averaged crater gas composition and output are well accounted for by closed system decompression of the basalt–gas mixture until the volcano–crust interface (～ 3 km depth), followed by open degassing and crystallization in the volcano conduits. However, both the low viscosity and high vesicularity of the basaltic magma permit bubble segregation and bubble foam growth at deep sill-like feeder discontinuities and at shallower physical boundaries (such as the volcano–crust interface) where the gas-rich aphyric basalt interacts with the unerupted crystal-rich and viscous magma drained back from the volcano conduits. Gas pressure build-up and bubble foam collapse at these boundaries will intermittently trigger the sudden upraise of CO₂-rich gas blobs that constitute the main driving force of the paroxysms. Deeper-sourced gas blobs, driving the most powerful explosions, will be the richest in CO₂ and have highest CO₂/SO₂ ratios. This mechanism is shown to account well for the dynamic, seismic and petrologic features of Stromboli's paroxysms and, hence, to provide a potential alternative interpretation for their genesis and their forecasting. Enhanced bubble foam leakage prior to a paroxysm, or foam emptying in several steps, should lead indeed to precursory upstream of CO₂-rich gas and increasing CO₂/SO₂ ratio in crater plume emissions. The recent detection of such signals prior to two explosions in December 2006 and March 2007 strongly supports this expectation and the model proposed in this study.     

Quantification of seep-related methane gas emissions at Tommeliten,North Sea

Tommeliten is a prominent methane seep area in the Central North Sea. Previous surveys revealed shallow gas-bearing sediments and methane gas ebullition into the water column. In this study, the in situ methane flux at Tommeliten is re-assessed and the potential methane transport to the atmosphere is discussed, with regards to the hydrographic setting and gas bubble modeling. We have compiled previous data, acquired new video and acoustic evidence of gas bubble release, and have measured the methane concentration, and its C-isotopic composition in the water column. Parametric subbottom sonar data reveal the three-dimensional extent of shallow gas and morphologic features relevant for gas migration. Five methane ebullition areas are identified and the main seepage area appears to be 21 times larger than previously estimated. Our video, hydroacoustic, subbottom, and chemical data suggest that ～1.5×10⁶ mol CH₄/yr (～26 tons CH₄/yr) of methane gas is being released from the seepage area of Tommeliten. Methane concentration profiles in the vicinity of the gas seeps show values of up to 268 nM (～100 times background) close to the seafloor. A decrease in δ¹³C-CH₄ values at 40 m water depth indicates an unknown additional biogenic methane source within the well oxygenated thermocline between 30 and 40 m water depth. Numerical modeling of the methane bubbles due to their migration and dissolution was performed to estimate the bubble-derived vertical methane transport, the fate of this methane in the water column, and finally the flux to the atmosphere. Modeling indicates that less than ～4% of the gas initially released at the seafloor is transported via bubbles into the mixed layer and, ultimately, to the atmosphere. However, because of the strong seasonality of mixing in the North Sea, this flux is expected to increase as mixing increases, and almost all of the methane released at the seafloor could be transferred into the atmosphere in the stormy fall and winter time.     

A new unsteady-state method of determining two-phase relative permeability illustrated by CO2-brine primary drainage in berea sandstone

This study presents a new unsteady-state method for measuring two-phase relative permeability by obtaining local values of the three key parameters (saturation, pressure drop, and phase flux) versus time during a displacement. These three parameters can be substituted to two-phase Darcy Buckingham equation to directly determine relative permeability. To obtain the first two, we use a medical X-ray Computed Tomography (CT) scanner to monitor saturation in time and space, and six differential pressure transducers to measure the overall pressure drop and the pressure drops of five individual sections (divided by four pressure taps on the core) continuously. At each scanning time, the local phase flux is obtained by spatially integrating the saturation profile and converting this to the flux using a fractional flow framework. One advantage of this local method over most previous methods is that the capillary end effect is experimentally avoided; this improvement is crucial for experiments using low viscosity fluids such as supercritical and gas phases. To illustrate the new method, we conduct five CO₂-brine primary drainage experiments in a 60.8 cm long and 116 mD Berea sandstone core at 20 °C and 1500 psi. In return, we obtain hundreds of unsteady-state CO₂ and brine relative permeability data points that are consistent with steady-state relative permeability data from the same experiments. Due to the large amount of relative permeability data obtained by the new unsteady-state method, the uncertainties of the exponents in the Corey-type fits decrease by up to 90% compared with the steady-state method.     

Himalayan metamorphic CO2 fluxes: Quantitative constraints from hydrothermal springs

Hot springs in the Marsyandi Valley, Nepal, vent CO₂ sourced from metamorphic fluids that mix with shallow groundwaters before degassing near the Earth's surface. The δ¹³C of spring waters ranges up to + 13‰, while that of the coexisting free gas phase is close to ? 4‰. Empirical and thermodynamic modelling of this isotopic fractionation suggests > 97 ± 1% CO₂ degassing. The calculated minimum total CO₂ degassing in the Marsyandi catchment is 5.4 × 10⁹ mol/yr from a Cl-based estimate of the spring water discharge to the Marsyandi River and the fraction of CO₂ degassed. Extrapolated to the whole of the Himalayas, this implies a probable minimum metamorphic CO₂ flux of 0.9 × 10¹² mol/yr, or ～ 13% of solid Earth CO₂ degassing. The calculated flux is a factor of three greater than the estimated CO₂ drawdown by silicate weathering in the Himalayas. Himalayan metamorphic degassing contributes a significant fraction of the global solid Earth CO₂ flux and implies that metamorphism may cause changes in long-term climate that oppose those resulting from the orogenic forcing of chemical weatherability.     

Migration of carrier and trace gases in the geosphere: an overview

The migration mechanisms of endogenous gases in the geosphere are defined in relation to the fluid-rock conditions and analyzed by basic transport equations. Upon examining the geological factors that influence the physical parameters in the equations in porous and fracture media, and considering the widespread high-permeability of deep subsurface rocks, in terms of fracture aperture, (orders of 10⁻² to 10¹ mm at depths of thousands meters, as suggested by recent crustal surveys) advection of carrier gases, in its several forms (gas-phase flow, water displacement by gas, gas slugs and bubbles) seems to represent a major migration process. Accordingly, in contrast with early views, the role of gas diffusion and water advection in the transport of endogenous gas to the Earth surface should be strongly minimized in many contexts. In a wide range of geological settings, carrier gases (CO₂, CH₄) may assume a dominant role in controlling transport and redistribution toward the Earth’s surface of trace gases (Rn, He). Bubble movement in fissured rocks seems to be an effective way of rapid (gas velocities in the order of 10⁰ to 10³ m per day) and long-distance gas migration. The evolution from bubble regimes to continuous phase flow and vice versa, as gas pressure and fracture width change, is the most suitable mechanism towards determining the surface geochemical processes of seismo-tectonic, environmental and geo-exploration relevance. The transport effectiveness of trace gases by a carrier gas has yet to be studied in quantitative terms. It is already clear, however, that further studies on the distribution and behavior of trace gases approaching the Earth’s surface may not be meaningful unless accompanied by carrier gas dynamics analyses.     

Mantle-derived carbon in Hercynian granites. Stable isotopes signatures and C/He associations in the thermomineral waters,N-Portugal

Na–HCO₃–CO₂-rich thermomineral waters issue in the N of Portugal, within the Galicia-Trás-os-Montes region, linked to a major NNE-trending fault, the so-called Penacova-Régua-Verin megalineament. Along this tectonic structure different occurrences of CO₂-rich thermomineral waters are found: Chaves hot waters (67 °C) and also several cold (16.1 °C) CO₂-rich waters. The δ²H and δ¹⁸O values of the thermomineral waters are similar to those of the local meteoric waters. The chemical composition of both hot and cold mineral waters suggests that water–rock reactions are mainly controlled by the amount of dissolved CO₂ (g) rather than by the water temperature. Stable carbon isotope data indicate an external CO₂ inorganic origin for the gas. δ¹³C_CO2 values ranging between ? 7.2‰ and ? 5.1‰ are consistent with a two-component mixture between crustal and mantle-derived CO₂. Such an assumption is supported by the ³He/⁴He ratios measured in the gas phase, are between 0.89 and 2.68 times the atmospheric ratio (Ra). These ratios which are higher than that those expected for a pure crustal origin (≈ 0.02 Ra), indicating that 10 to 30% of the He has originated from the upper mantle. Release of deep-seated fluids having a mantle-derived component in a region without recent volcanic activity indicates that extensive neo-tectonic structures originating during the Alpine Orogeny are still active (i.e., the Chaves Depression).     

New insights into the origin and evolution of Lake Vida,McMurdo Dry Valleys,Antarctica — A noble gas study in ice and brines

Unlike other lakes in the McMurdo Dry Valleys, Antarctica, Lake Vida has a thick (~ 19 m) ice cover sealing a liquid brine body of unusually high salinity (~ 245 g/L) from the atmosphere. To constrain the conditions under which the atypical Lake Vida ice cover formed and evolved, 19 ice samples were collected down to a depth of ~ 14 m, together with three brine samples trapped in the ice at ~ 16 m for analysis of helium, neon, argon, krypton, and xenon concentrations. The broad pattern of noble gas concentrations for Lake Vida samples is fundamentally different from that of air saturated water (ASW) at 0 °C and an elevation of 340 m for salinities of 0 (ice) and 245 g/L (brine). Overall, ice samples are enriched in He and depleted in Ne with saturation relative to ASW averages of 1.38 and 0.82, respectively, and strongly depleted in Ar, Kr, and Xe with relative saturations of 0.10, 0.06, and 0.05, respectively. By contrast, brine samples are generally depleted in He and Ne (relative saturation averages of 0.33 and 0.27, respectively) but enriched in Ar, Kr, and Xe, with relative saturation averages of 1.45, 3.15, and 8.86, respectively. A three-phase freezing partitioning model generating brine, ice and bubble concentrations for all stable noble gases was tested and compared with our data. Measured brine values are best reproduced for a salinity value of 175 g/L, a pressure of 1.1 atm, and a bubble volume of 20 cm³ kg^?1. Sensitivity tests for ice + bubble samples show an ideal fit for bubble volumes of ~ 1–2 cm³ kg^?1. Our results show that the conditions under which ice and brine formed and evolved at Lake Vida are significantly different from other ice-covered lakes in the area. Our brine data suggest that Lake Vida may be transitioning from a wet to a dry-based lake, while the ice + bubble data suggest at least partial re-equilibration of residual liquid with the atmosphere as ice forms at the top of Lake Vida ice cover.     

Deep pooling of low degree melts and volatile fluxes at the 85°E segment of the Gakkel Ridge: Evidence from olivine-hosted melt inclusions and glasses

We present new analyses of volatile, major, and trace elements for a suite of glasses and melt inclusions from the 85°E segment of the ultra-slow spreading Gakkel Ridge. Samples from this segment include limu o pele and glass shards, proposed to result from CO₂-driven explosive activity. The major element and volatile compositions of the melt inclusions are more variable and consistently more primitive than the glass data. CO₂ contents in the melt inclusions extend to higher values (167–1596 ppm) than in the co-existing glasses (187–227 ppm), indicating that the melt inclusions were trapped at greater depths. These melt inclusions record the highest CO₂ melt concentrations observed for a ridge environment. Based on a vapor saturation model, we estimate that the melt inclusions were trapped between seafloor depths (～ 4 km) and ～ 9 km below the seafloor. However, the glasses are all in equilibrium with their eruption depths, which is inconsistent with the rapid magma ascent rates expected for explosive activity. Melting conditions inferred from thermobarometry suggest relatively deep (25–40 km) and cold (1240°–1325 °C) melting conditions, consistent with a thermal structure calculated for the Gakkel Ridge. The water contents and trace element compositions of the melt inclusions and glasses are remarkably homogeneous; this is an unexpected result for ultra-slow spreading ridges, where magma mixing is generally thought to be less efficient based on the assumption that steady-state crustal magma chambers are absent in these environments. All melts can be described by a single liquid line of descent originating from a pooled melt composition that is consistent with the aggregate melt calculated from a geodynamic model for the Gakkel Ridge. These data suggest a model in which deep, low degree melts are efficiently pooled in the upper mantle (9–20 km depth), after which crystallization commences and continues during ascent and eruption. Based on our melting model and the assumption that CO₂ is perfectly incompatible, we show that the highest CO₂ concentrations of the melt inclusions (～ 1600 ppm) are consistent with the calculated CO₂ concentrations of primary undegassed melts. The highest measured CO₂/Nb ratio (443) of Gakkel Ridge melt inclusions predicts a mantle CO₂ content of 134 ppm and would result in a global ridge flux of 2.0 × 10¹² mol CO₂/yr.     

Multifluid flow in bedded porous media: laboratory experiments and numerical simulations

Understanding light nonaqueous-phase liquid (LNAPL) movement in heterogeneous vadose environments is important for effective remediation design. We investigated LNAPL movement near a sloping fine- over coarse-grained textural interface, forming a capillary barrier. LNAPL flow experiments were performed in a glass chamber (50 cm×60 cm×1.0 cm) using two silica sands (12/20 and 30/40 sieve sizes). Variable water saturations near the textural interface were generated by applying water uniformly to the sand surface at various flow rates. A model LNAPL (Soltrol® 220) was subsequently released at two locations at the sand surface. Visible light transmission was used to quantitatively determine water saturations prior to LNAPL release and to observe LNAPL flow paths. Numerical simulations were performed using the Subsurface Transport Over Multiple Phases (STOMP) simulator, employing two nonhysteretic relative permeability–saturation–pressure (k–S–P) models. LNAPL movement strongly depended on the water saturation in the fine-grained sand layer above the textural interface. In general, reasonable agreement was found between observed and predicted water saturations near the textural interface and LNAPL flow paths. Discrepancies between predictions based on the van Genuchten/Mualem (VGM) and Brooks–Corey/Burdine (BCB) k–S–P models existed in the migration speed of the simulated LNAPL plume and the LNAPL flow patterns at high water saturation above the textural interface. In both instances, predictions based on the BCB model agreed better with experimental observations than predictions based on the VGM model. The results confirm the critical role water saturation plays in determining LNAPL movement in heterogeneous vadose zone environments and that accurate prediction of LNAPL flow paths depends on the careful selection of an appropriate k–S–P model.     

Degassing of stratified magma by compositional convection

A new model is proposed for passive degassing from sub-volcanic magma chambers. The water content in stably stratified shallow magma chamber will be equated to its solubility at the upper boundary by convection. Water from a lower layer high in water content can enrich the contact zone of the upper layer and lead to further convective overturn of this boundary layer. A complete set of equations describing convection with bubble formation and dissolution is reduced to a simplified form by assuming a small bubble content. The development and pattern of flow driven by vesiculation is modeled numerically in a 2D magma chamber for relatively low Raleigh numbers (5×10⁵). Bubbles rising from the magma will collect near the roof in a layer of 8–10 vol% and then escape upward to fumaroles. The Stokes flux of bubbles escaping from an andesitic magma with viscosity 10⁴ P and a top surface of about 500×500 m corresponds with observed total magmatic water fluxes of 35 kg/s. Pressure within the chamber is buffered by elastic (and local visco-elastic) deformations in the solid rocks bounding the chamber to the range between ambient and close to lithostatic values. In a chamber closed to fresh magma inputs, the decrease in volume due to such gentle volatile escape lowers the reference pressure. Bubbles flux from the lower layer induced by variation of the saturation level around stratification boundary may be efficient mechanism for the water transport between layers.     

The pre-eruptive volatile contents of recent basaltic and pantelleritic magmas at Pantelleria (Italy)

Pantelleria Island, located in the Sicily Channel Rift Zone (Italy), is the type locality for the peralkaline rhyolitic rocks called pantellerites. In the last 50 ka, after the large Green Tuff caldera-forming eruption, volcanic activity at Pantelleria has consisted of effusive and explosive eruptions mostly vented inside and along the rim of the caldera and producing silicic lava flows, lava domes and poorly dispersed pantelleritic pumice fall deposits. Basaltic cinder cones and lava flows are only present outside the caldera in the NW sector of the island. The most recent basaltic (Cuddie Rosse, ～ 20 ka) and pantelleritic (Cuddia Randazzo and Cuddia del Gallo, ～ 6 ka) pyroclastic products were sampled to investigate magmatic volatile contents through the study of melt inclusions.The melt inclusions in pyroxene and olivine phenocrysts of Cuddie Rosse scoriae have an alkali basalt composition. The dissolved volatiles comprise 0.9–1.6 wt.% H₂O, several hundred ppm of CO₂, 1600–2000 ppm of sulphur and 500–900 ppm of chlorine. The water–carbon dioxide couple gives a confining pressure ～ 2 kbar prior to the eruption. This result indicates that episodes of magma ponding and crystallization occurred in the upper crust prior to eruption. The melt inclusions in feldspar, fayalite and aenigmatite phenocrysts of Cuddia del Gallo and Cuddia Randazzo pumice have a pantelleritic composition (Agpaitic Indices 1.3–2.1), up to 4.4 wt.% H₂O, 8700 ppm Cl, 6000 ppm F, and CO₂ below the detection limit. Sulphur averaging 420 ppm has been measured in Cuddia Randazzo melt inclusions. These data indicate relatively high volatile contents for these low-energy Strombolian-type eruptions. Melt inclusions in Cuddia del Gallo pumice show the most evolved composition (Agpaitic Indices 2–2.1) and the highest volatile content, in agreement with fluid saturation conditions in the magma chamber prior to the eruption. This implies a confining pressure of ～ 1 kbar for the top of the pantelleritic reservoir. The composition of melt inclusions and mineralogical assemblage of Cuddia Randazzo pumice indicate that it has a lower evolutionary degree (Agpaitic Indices 1.3–1.8) and lower pre-eruptive Cl and H₂O contents than Cuddia del Gallo pumice. An increase in pressure due to the exsolution of volatiles in the upper part of the pantelleritic reservoir may have triggered the Cuddia del Gallo explosive eruption. Evidence of widespread pre-eruptive mingling between trachytes and pantellerites suggests that the intrusion of trachytic magma into the pantelleritic reservoir likely played a major role in destabilizing the magma system just prior to the Cuddia Randazzo event.     

Cave air ventilation and CO2 outgassing by radon-222 modeling: How fast do caves breathe?

In general, the rate and timing of calcite precipitation is in part affected by variations in cave air CO₂ concentrations. Knowledge of cave ventilation processes is required to quantify the effect variations in CO₂ concentrations have on speleothem deposition rates and thus paleoclimate records. In this study we use radon-222 (²²²Rn) as a proxy of ventilation to estimate CO₂ outgassing from the cave to the atmosphere, which can be used to infer relative speleothem deposition rates. Hollow Ridge Cave, a wild cave preserve in Marianna, Florida, is instrumented inside and out with multiple micro-meteorological sensor stations that record continuous physical and air chemistry time-series data. Our time series datasets indicate diurnal and seasonal variations in cave air ²²²Rn and CO₂ concentrations, punctuated by events that provide clues to ventilation and drip water degassing mechanisms. Average cave air ²²²Rn and CO₂ concentrations vary seasonally between winter (²²²Rn = 50 dpm L^? 1, where 1 dpm L^? 1 = 60 Bq m^? 3; CO₂ = 360 ppmv) and summer (²²²Rn = 1400 dpm L^? 1; CO₂ = 3900 ppmv). Large amplitude diurnal variations are observed during late summer and autumn (²²²Rn = 6 to 581 dpm L^? 1; CO₂ = 360 to 2500 ppmv).We employ a simple first-order ²²²Rn mass balance model to estimate cave air exchange rates with the outside atmosphere. Ventilation occurs via density driven flow and by winds across the entrances which create a ‘venturi’ effect. The most rapid ventilation occurs 25 m inside the cave near the entrance: 45 h^? 1 (1.33 min turnover time). Farther inside (175 m) exchange is slower and maximum ventilation rates are 3 h^? 1 (22 min turnover time). We estimate net CO₂ flux from the epikarst to the cave atmosphere using a CO₂ mass balance model tuned with the ²²²Rn model. Net CO₂ flux from the epikarst is highest in summer (72 mmol m^? 2 day^? 1) and lowest in late autumn and winter (12 mmol m^? 2 day^? 1). Modeled ventilation and net CO₂ fluxes are used to estimate net CO₂ outgassing from the cave to the atmosphere. Average net CO₂ outgassing is positive (net loss from the cave) and is highest in late summer and early autumn (about 4 mol h^? 1) and lowest in winter (about 0.5 mol h^? 1). Modeling of ventilation, net CO₂ flux from the epikarst, and CO₂ outgassing to the atmosphere from cave monitoring time-series can help better constrain paleoclimatic interpretations of speleothem geochemical records.     

Earthquakes produce carbon dioxide in crustal faults

Fourier transform infrared (FTIR) microanalysis of pseudotachylytes (i.e. friction-induced melts produced by seismic slip) from the Nojima fault (Japan) reveals that earthquakes almost instantaneously expel 99 wt.% of the wall rock CO₂ content. Carbon is exsolved because it is supersaturated in the friction melts. By extrapolation to a crustal-scale fault rupture, large events such as the M7.2 Kobe earthquake (1995) may yield a total production of 1.8 to 3.4 × 10³ tons CO₂ within a few seconds. This extraordinary release of CO₂ can cause a flash fluid pressure increase in the fault plane, and therefore enhance earthquake slip or trigger aftershocks; it may also explain the anomalous discharge of carbon monitored in nearby fault springs after large earthquakes. Because carbon saturation in silicate melts is pressure-dependent, FTIR can be used as a new tool to constrain the maximum depth of pseudotachylyte formation in exhumed faults.     

Partial pressure and air–sea CO2 flux in the Aegean Sea during February 2006

Data on the distribution of fCO₂ were obtained during a cruise in the Aegean Sea during February 2006. The fCO₂ of surface water (fCO₂^sw) was lower than the atmospheric fCO₂ (fCO₂^atm) throughout the area surveyed and ΔfCO₂ values varied from ?34 to ?61 μatm. The observed under-saturation suggests that surface waters in the Aegean represent a sink for atmospheric CO₂ during the winter of 2006. Higher fCO₂^sw values were recorded in the ‘less warm’ and ‘less saline’ shallow northernmost part of the Aegean Sea implying that the lower seawater temperature and salinity in this area play a crucial role in the spatial distribution of fCO₂^sw.A first estimate of the magnitude of the air–sea CO₂ exchange and the potential role of the Aegean Sea in the transfer of atmospheric CO₂ was also obtained. The air–sea CO₂ fluxes calculated using different gas transfer formulations showed that during February 2006, the Aegean Sea absorbs atmospheric CO₂ at a rate ranging from ?6.2 to ?11.8 mmol m^?2 d^?1 with the shipboard recorded wind speeds and at almost half rate (?3.5 to ?5.5 mmol m^?2 d^?1) with the monthly mean model-derived wind speed. Compared to recent observations from other temperate continental shelves during winter period, the Aegean Sea acts as a moderate to rather strong sink for atmospheric CO₂.Further investigations, including intensive spatial and temporal high-resolution observations, are necessary to elucidate the role of the Aegean Sea in the process of transfer of atmospheric CO₂ into the deep horizons of the Eastern Mediterranean.     

Multiphysics modeling of CO2 sequestration in a faulted saline formation in Italy

The present work describes the results of a modeling study addressing the geological sequestration of carbon dioxide (CO₂) in an offshore multi-compartment reservoir located in Italy. The study is part of a large scale project aimed at implementing carbon capture and storage (CCS) technology in a power plant in Italy within the framework of the European Energy Programme for Recovery (EEPR). The processes modeled include multiphase flow and geomechanical effects occurring in the storage formation and the sealing layers, along with near wellbore effects, fault/thrust reactivation and land surface stability, for a CO₂ injection rate of 1 × 10⁶ ton/a. Based on an accurate reproduction of the three-dimensional geological setting of the selected structure, two scenarios are discussed depending on a different distribution of the petrophysical properties of the formation used for injection, namely porosity and permeability. The numerical results help clarify the importance of: (i) facies models at the reservoir scale, properly conditioned on wellbore logs, in assessing the CO₂ storage capacity; (ii) coupled wellbore-reservoir flow in allocating injection fluxes among permeable levels; and (iii) geomechanical processes, especially shear failure, in constraining the sustainable pressure buildup of a faulted reservoir.     

Quantification of the gas mass emitted during single explosions on Stromboli with the SO2 imaging camera

We performed measurements using an SO₂ imaging camera of the SO₂ gas mass emitted during five discrete explosive events on Stromboli volcano on 3 October 2006. The SO₂ gas mass released during discrete explosions was 15–40 kg per explosion, producing 3–8% of the total daily SO₂ gas emission, demonstrating that in terms of gas flux Strombolian explosions are a second-order phenomenon compared with quiescent degassing. Using the typical gas composition measured with OP-FTIR allows us to determine the total gas mass released during an explosion as 360–960 kg with a volume of 1500–4100 m³ at 1 bar. At the probable source pressure of gas slug formation of 75 MPa this gas amount would occupy a volume equivalent to a sphere with a radius of 0.8–1 m, comparable with estimates of Stromboli's conduit geometry.     

Carbon remineralization in the Amazon–Guianas tropical mobile mudbelt: A sedimentary incinerator

The Amazon River spawns a vast mobile mudbelt extending ∼1600 km from the equator to the Orinoco delta. Deposits along the Amazon–Guianas coastline are characterized by some of the highest C_org remineralization rates reported for estuarine, deltaic, or shelf deposits, however, paradoxically, except where stabilized by mangroves or intertidal algal mats, they are usually suboxic and nonsulfidic. A combination of tides, wind-driven waves, and coastal currents forms massive fluid muds and mobile surface sediment layers ∼0.5–2 m thick which are dynamically refluxed and frequently reoxidized. Overall, the seabed functions as a periodically mixed batch reactor, efficiently remineralizing organic matter in a gigantic sedimentary incinerator of global importance. Amazon River material entering the head of this dynamic dispersal system carries an initial terrestrial sedimentary C_org loading of ∼ 0.7 mg C m⁻² particle surface area. Total C_org loading is lowered to ∼ 0.2 mg C m⁻² in the proximal delta topset, ∼60–70% of which remains of terrestrial origin. Loading decreases further to 0.12–0.14 mg C m⁻² (∼60% terrestrial) in mudbanks ∼600 km downdrift along French Guiana, values comparable to those found in the oligotrophic deepsea. DOC/ΣCO₂ ratios in pore waters of French Guiana mudbanks indicate that >90% of metabolized organic substrates are completely oxidized. Within the Amazon delta topset at the head of the dispersal system, both terrestrial and marine organic matter contribute substantially to early diagenetic remineralization, although reactive marine substrate dominates (∼60–70%). The conditional rate constant for terrestrial C_org in the delta topset is ∼0.2 a⁻¹. As sedimentary C_org is depleted during transit, marine sources become virtually the exclusive substrate for remineralization except very near the mangrove shoreline. The δ¹³C and Δ¹⁴C values of pore water ΣCO₂ in mudbanks demonstrate that the primary source of remineralized organic matter within ∼1 km of shore is a small quantity of bomb signature marine plankton (+80‰). Thus, fresh marine organic material is constantly entrained into mobile deposits and increasingly drives early diagenetic reactions along the transit path. Relatively refractory terrestrial C_org is lost more slowly but steadily during sedimentary refluxing and suboxic diagenesis. Amazon Fan deposits formed during low sea level stand largely bypassed this suboxic sedimentary incinerator and stored material with up to ∼3X the modern high stand inner shelf C_org load (Keil et al., 1997b. Proceedings of the Ocean Drilling Program, Scientific Results. Vol. 155. pp. 531–537). Sedimentary dynamics, including frequency and magnitude of remobilization, and the nature of dispersal systems are clearly key controls on diagenetic processes, biogeochemical cycling, and global C storage along the continental margins.     

Space weather in the thermospheric–ionospheric domain over the Brazilian region: Climatology of ionospheric plasma bubbles in the subequatorial and low-latitude region

The climatology of ionospheric plasma bubbles is studied here by means of a comparison of the frequency of occurrence of the spread-F/plasma bubble events over the South American region using the images from two OI 630 nm imager systems located at the subequatorial station São João do Cariri—CA (7.4 S, 36.5 W, 20 S dip) and the low-latitude station Cachoeira Paulista—CP (22.5 S, 45 W, 33 S dip) in Brazil during the years of 2004 and 2005. The results are discussed in the light of current theory and geomagnetic parameters of the two observation stations.     

Volatile and trace elements in basaltic glasses from Samoa: Implications for water distribution in the mantle

We report volatile (H₂O, CO₂, F, S, Cl) and trace element data for submarine alkalic basalt glasses from the three youngest Samoan volcanoes, Ta'u, Malumalu and Vailulu'u. Most samples are visibly sulfide saturated, so have likely lost some S during fractionation. Cl / K ratios (0.04–0.15) extend to higher values than pristine MORBs, but are suspected to be partly due to source differences since Cl / K roughly varies as a function of ⁸⁷Sr / ⁸⁶Sr. There are no resolvable differences in the relative enrichment of F among sources, and compatibility of F during mantle melting is established to be nearly identical to Nd. Shallow degassing has affected CO₂ in all samples, and H₂O only in the most shallowly erupted samples from Vailulu'u. Absolute water contents are high for Samoa (0.63–1.50 wt.%), but relative enrichment of water compared to equally incompatible trace elements (Ce, La) is low and falls entirely below normal MORB values. H₂O / Ce (58–157) and H₂O / La (120–350) correlate inversely with ⁸⁷Sr / ⁸⁶Sr compositions (0.7045–0.7089). This leads us to believe that, because of very fast diffusion of hydrogen in mantle minerals, recycled lithospheric material with high initial water and trace element content will lose water to the drier ambient mantle during storage within the inner Earth. The net result is the counter-intuitive appearance of greater dehydration with greater mantle enrichment. We expect that subducted slabs will experience a two-stage dehydration history, first within subduction zones and then in the ambient mantle during long-term convective mixing.     

UV camera measurements of fumarole field degassing (La Fossa crater,Vulcano Island)

The UV camera is becoming an important new tool in the armory of volcano geochemists to derive high time resolution SO₂ flux measurements. Furthermore, the high camera spatial resolution is particularly useful for exploring multiple-source SO₂ gas emissions, for instance the composite fumarolic systems topping most quiescent volcanoes. Here, we report on the first SO₂ flux measurements from individual fumaroles of the fumarolic field of La Fossa crater (Vulcano Island, Aeolian Island), which we performed using a UV camera in two field campaigns: in November 12, 2009 and February 4, 2010. We derived ~ 0.5 Hz SO₂ flux time-series finding fluxes from individual fumaroles, ranging from 2 to 8.7 t d^?1, with a total emission from the entire system of ~ 20 t d^?1 and ~ 13 t d^?1, in November 2009 and February 2010 respectively. These data were augmented with molar H₂S/SO₂, CO₂/SO₂ and H₂O/SO₂ ratios, measured using a portable MultiGAS analyzer, for the individual fumaroles. Using the SO₂ flux data in tandem with the molar ratios, we calculated the flux of volcanic species from individual fumaroles, and the crater as a whole: CO₂ (684 t d^?1 and 293 t d^?1), H₂S (8 t d^?1 and 7.5 t d^?1) and H₂O (580 t d^?1 and 225 t d^?1).