Assessment of soil carbon pool,carbon sequestration and soil CO<Subscript>2</Subscript> flux in unreclaimed and reclaimed coal mine spoils

Surface coal mining inevitably deforests the land, reduces carbon (C) pool and generates different land covers. To re-establish the ecosystem C pool, post-mining lands are often afforested with fast-growing trees. A field study was conducted in the 5-year-old unreclaimed dump and reclaimed coal mine dump to assess the changes in soil CO₂ flux and compared with the reference forest site. Changes in soil organic carbon (SOC) and total nitrogen stocks were estimated in post-mining land. Soil CO₂ flux was measured using close dynamic chamber method, and the influence of environmental variables on soil CO₂ flux was determined. Woody biomass C and SOC stocks of the reference forest site were threefold higher than that of 5-year-old reclaimed site. The mean soil CO₂ flux was highest in 5-year-old reclaimed dump (2.37 μmol CO₂ m^?2 s^?1) and lowest in unreclaimed dump (0.21 μmol CO₂ m^?2 s^?1). Soil CO₂ flux was highly influenced by environmental variables, where soil temperature positively influenced the soil CO₂ flux, while soil moisture, relative humidity and surface CO₂ concentration negatively influenced the soil CO₂ flux. Change in soil CO₂ flux under different land cover depends on plant and soil characteristics and environmental variables. The study concluded that assessment of soil CO₂ flux in post-mining land is important to estimate the potential of afforestation to combat increased emission of soil CO₂ at regional and global scale.     

Energy and CO<Subscript>2</Subscript> exchanges and influencing factors in spring wheat ecosystem along the Heihe River,northwestern China

Spring wheat (Triticum aestivum Linn.) is an important crop for food security in the desert-oasis farmland in the middle reaches of the Heihe River in northwestern China. We measured fluxes using eddy covariance and meteorological parameters to explore the energy fluxes and the relationship between CO₂ flux and climate change in this region during the wheat growing seasons in 2013 and 2014. The energy balance closures were 70.5% and 72.7% in the 2013 and 2014 growing season, respectively. The wheat ecosystem had distinct seasonal and diurnal dynamics of CO₂ fluxes with U-shaped curves. The accumulated net ecosystemic CO₂ exchanges (NEE) were -111.6 and -142.2 g C/m² in 2013 and 2014 growing season, respectively. The ecosystem generally acted as a CO₂ sink during the growing season but became a CO₂ source after the wheat harvest. A correlation analysis indicated that night-time CO₂ fluxes were exponentially dependent on air temperature and soil temperature at a depth of 5 cm but were not correlated with soil-water content, water-vapour pressure, or vapour-pressure deficit. CO₂ flux was not correlated with the meteorological parameters during daytime. However, irrigation and precipitation, may complicate the response of CO₂ fluxes to other meteorological parameters.     

Bayesian hierarchical models for soil CO<Subscript>2</Subscript> flux and leak detection at geologic sequestration sites

Proper characterizations of background soil CO₂ respiration rates are critical for interpreting CO₂ leakage monitoring results at geologic sequestration sites. In this paper, a method is developed for determining temperature-dependent critical values of soil CO₂ flux for preliminary leak detection inference. The method is illustrated using surface CO₂ flux measurements obtained from the AmeriFlux network fit with alternative models for the soil CO₂ flux versus soil temperature relationship. The models are fit first to determine pooled parameter estimates across the sites, then using a Bayesian hierarchical method to obtain both global and site-specific parameter estimates. Model comparisons are made using the deviance information criterion (DIC), which considers both goodness of fit and model complexity. The hierarchical models consistently outperform the corresponding pooled models, demonstrating the need for site-specific data and estimates when determining relationships for background soil respiration. A hierarchical model that relates the square root of the CO₂ flux to a quadratic function of soil temperature is found to provide the best fit for the AmeriFlux sites among the models tested. This model also yields effective prediction intervals, consistent with the upper envelope of the flux data across the modeled sites and temperature ranges. Calculation of upper prediction intervals using the proposed method can provide a basis for setting critical values in CO₂ leak detection monitoring at sequestration sites.     

The influence of degradation of the swamp and alpine meadows on CH<Subscript>4</Subscript> and CO<Subscript>2</Subscript> fluxes on the Qinghai-Tibetan Plateau

CH₄ and CO₂ fluxes from a high-cold swamp meadow and an alpine meadow on the Qinghai-Tibetan Plateau, subject to different degrees of degradation, were measured over a 12-month period. Air temperature, soil temperature and moisture, and the depths of the water table and thawing-freezing layer were determined. For swamp meadows, the greater the degradation, the lesser the carbon efflux. CH₄ emissions at the nondegraded swamp meadow site were 1.09–3.5 and 2.5–11.27 times greater, and CO₂ emissions 1.08–1.69 and 1.41–4.43 times greater, respectively, than those from moderately and severely degraded sites. For alpine meadows, the greater the degradation, the greater the CH₄ consumption and CO₂ emissions. CH₄ consumption at the severely degraded alpine meadow site was 6.6–21 and 1.1–5.25 times greater, and CO₂ emissions 1.05–78.5 and 1.04–6.28 times greater, respectively, than those from the nondegraded and moderately degraded sites. The CH₄ and CO₂ fluxes at both sites were significantly correlated (R ² > 0.59, P < 0.05) with air temperature, soil temperature, and topsoil (0–5 cm depth) moisture, indicating these to be the main environmental factors affecting such fluxes.     

Analysis of the spatial and temporal changes in soil CO<Subscript>2</Subscript> flux in alpine meadow of Qilian Mountain

Field experiments on the CO₂ flux of alpine meadow soil in the Qilian Mountain were conducted along the elevation gradient during the growing season of 2004 and 2005. The soil CO₂ flux was measured using the Li-6400-09 soil respiration chamber attached to the Li-6400 portable photosynthesis system. The effects of water and heat and roots on the soil CO₂ flux were statistically analyzed. The results show that soil CO₂ flux along the elevation gradient gradually decreases. The soil CO₂ flux was low at night, with lowest value occurring between 0200 and 0600 hours, started to rise rapidly during 0700–0830 hours, and then descend during 1600–1830 hours. The peak CO₂ efflux appears during 1100–1600 hours. The diurnal average of soil CO₂ efflux was between 0.56 ± 0.32 and 2.53 ± 0.76 μmol m⁻² s⁻¹. Seasonally, soil CO₂ fluxes are relatively high in summer and autumn and low in spring and winter. The soil CO₂ efflux, from the highest to the lowest in the ranking order, occurred in July and August (4.736 μmol m⁻² s⁻¹), June and September, and May and October, respectively. The soil CO₂ efflux during the growing season is positively correlated with soil temperature, root biomass and soil water content.     

Diel Aquatic CO<Subscript>2</Subscript> System Dynamics of a Bermudian Mangrove Environment

Mangrove ecosystems play an important, but understudied, role in the cycling of carbon in tropical and subtropical coastal ocean environments. In the present study, we examined the diel dynamics of seawater carbon dioxide (CO₂) and dissolved oxygen (DO) for a mangrove-dominated marine ecosystem (Mangrove Bay) and an adjacent intracoastal waterway (Ferry Reach) on the island of Bermuda. Spatial and temporal trends in seawater carbonate chemistry and associated variables were assessed from direct measurements of dissolved inorganic carbon, total alkalinity, dissolved oxygen (DO), temperature, and salinity. Diel pCO₂ variability was interpolated across hourly wind speed measurements to determine variability in daily CO₂ fluxes for the month of October 2007 in Bermuda. From these observations, we estimated rates of net sea to air CO₂ exchange for these two coastal ecosystems at 59.8 ± 17.3 in Mangrove Bay and 5.5 ± 1.3 mmol m⁻² d⁻¹ in Ferry Reach. These results highlight the potential for large differences in carbonate system functioning and sea-air CO₂ flux in adjacent coastal environments. In addition, observation of large diel variability in CO₂ system parameters (e.g., mean pCO₂: 390–2,841 μatm; mean pH_T: 8.05–7.34) underscores the need for careful consideration of diel cycles in long-term sampling regimes and flux estimates.     

Human impact on the historical change of CO<Subscript>2</Subscript>degassing flux in River Changjiang

The impact of water quality changes in River Changjiang (formally known as the Yangtze River) on dissolved CO₂ and silicate concentrations and seasonal carbon flux in the past several decades (1960s–2000) was evaluated, based on monitoring data from hydrographic gauge. It was found that dissolved CO₂ and silicate in Changjiang decreased dramatically during this decades, as opposed to a marked increase in nutrient (e.g. NO₃ ^-) concentrations. Our analyses revealed that dissolved CO₂ in Changjiang was over-saturated with the atmosphere CO₂, and its concentration had showed a declining trend since the 1960s, despite that fluvial DIC flux had maintained stable. Analysis results also suggested that the decrease in dissolved CO₂ concentration was attributed to changes on the riverine trophic level and river damming activities in the Changjiang drainage basin. Due to the economic innovation (e.g. agriculture and industry development) across the Changjiang watershed, fertilizers application and river regulations have significantly altered the original state of the river. Its ecosystem and hydrological condition have been evolving toward the "lacustrine/reservoir" autotrophic type prevailing with plankton. Accordingly, average CO₂ diffusing flux to the atmosphere from the river had been reduced by three-fourth from the 1960s to 1990s, with the flux value being down to 14.2 mol.m^-2.yr^-1 in the 1990s. For a rough estimate, approximately 15.3 Mt of carbon was degassed annually into the atmosphere from the entire Changjiang drainage basin in the 1990s.     

Carbon Dioxide Fluxes and Their Environmental Control in a Reclaimed Coastal Wetland in the Yangtze Estuary

Large areas of natural coastal wetlands have suffered severely from human-driven damages or conversions (e.g., land reclamations), but coastal carbon flux responses in reclaimed wetlands are largely unknown. The lack of knowledge of the environmental control mechanisms of carbon fluxes also limits the carbon budget management of reclaimed wetlands. The net ecosystem exchange (NEE) in a coastal wetland at Dongtan of Chongming Island in the Yangtze estuary was monitored throughout 2012 using the eddy covariance technique more than 14 years after this wetland was reclaimed using dykes to stop tidal flooding. The driving biophysical variables of NEE were also examined. The results showed that NEE displayed marked diurnal and seasonal variations. The monthly mean NEE showed that this ecosystem functioned as a CO₂ sink during 9 months of the year, with a maximum value in September (?101.2 g C m^?2) and a minimum value in November (?8.2 g C m^?2). The annual CO₂ balance of the reclaimed coastal wetland was ?558.4 g C m^?2 year^?1. The ratio of ecosystem respiration (ER) to gross primary production (GPP) was 0.57, which suggests that 57 % of the organic carbon assimilated by wetland plants was consumed by plant respiration and soil heterotrophic respiration. Stepwise multiple linear regressions suggested that temperature and photosynthetically active radiation (PAR) were the two dominant micrometeorological variables driving seasonal variations in NEE, while soil moisture (M _s) and soil salinity (PS_s) played minor roles. For the entire year, PAR and daytime NEE were significantly correlated, as well as temperature and nighttime NEE. These nonlinear relationships varied seasonally: the maximum ecosystem photosynthetic rate (A _max), apparent quantum yield (?), and Q ₁₀ reached their peak values during summer (17.09 μmol CO₂?m^?2 s^?1), autumn (0.13 μmol CO₂?μmol^?1 photon), and spring (2.16), respectively. Exceptionally high M _s or PS_s values indirectly restricted ecosystem CO₂ fixation capacity by reducing the PAR sensitivity of the NEE. The leaf area index (LAI) and live aboveground biomass (AGB_L) were significantly correlated with NEE during the growing season. Although the annual net CO₂ fixation rate of the coastal reclaimed wetland was distinctly lower than the unreclaimed coastal wetland in the same region, it was quite high relative to many inland freshwater wetlands and estuarine/coastal wetlands located at latitudes higher than this site. Thus, it is concluded that although the net CO₂ fixation capacity of the coastal wetland was reduced by land reclamation, it can still perform as an important CO₂ sink.     

Site characterization for CO<Subscript>2</Subscript> geologic storage and vice versa: the Frio brine pilot,Texas, USA as a case study

Careful site characterization is critical for successful geologic storage of carbon dioxide (CO₂) because of the many physical and chemical processes impacting CO₂ movement and containment under field conditions. Traditional site characterization techniques such as geological mapping, geophysical imaging, well logging, core analyses, and hydraulic well testing provide the basis for judging whether or not a site is suitable for CO₂ storage. However, only through the injection and monitoring of CO₂ itself can the coupling between buoyancy flow, geologic heterogeneity, and history-dependent multi-phase flow effects be observed and quantified. CO₂ injection and monitoring can therefore provide a valuable addition to the site-characterization process. Additionally, careful monitoring and verification of CO₂ plume development during the early stages of commercial operation should be performed to assess storage potential and demonstrate permanence. The Frio brine pilot, a research project located in Dayton, Texas (USA) is used as a case study to illustrate the concept of an iterative sequence in which traditional site characterization is used to prepare for CO₂ injection and then CO₂ injection itself is used to further site-characterization efforts, constrain geologic storage potential, and validate understanding of geochemical and hydrological processes. At the Frio brine pilot, in addition to traditional site-characterization techniques, CO₂ movement in the subsurface is monitored by sampling fluid at an observation well, running CO₂-saturation-sensitive well logs periodically in both injection and observation wells, imaging with crosswell seismic in the plane between the injection and observation wells, and obtaining vertical seismic profiles to monitor the CO₂ plume as it migrates beyond the immediate vicinity of the wells. Numerical modeling plays a central role in integrating geological, geophysical, and hydrological field observations.     

Potential risks to freshwater resources as a result of leakage from CO<Subscript>2</Subscript> geological storage: a batch-reaction experiment

In assessing the feasibility of widespread deployment of CO₂ geological storage, it is prudent to first assess potential consequences of an error or accident that could lead to CO₂ leakage into groundwater resources above a sequestration interval. Information about the sensitivity of the groundwater system to introduction of CO₂ is needed in order to design groundwater monitoring program. A laboratory-batch experiment was conducted to explore the range of CO₂ impact on groundwater quality of a spectrum of representative aquifers, in the Gulf Coast region, USA. Results show that CO₂ elevated concentrations of many cations within hours or days. Two types of cations were recognized according to their concentration trends. Type I cations—Ca, Mg, Si, K, Sr, Mn, Ba, Co, B, Zn—rapidly increased following initial CO₂ flux and reached stable concentrations before the end of the experiment. Type II cations—Fe, Al, Mo, U, V, As, Cr, Cs, Rb, Ni and Cu—increased at the start of CO₂ flux, but declined, in most cases, to levels lower than pre-CO₂ concentrations. Dissolution of dolomite and calcite caused the largest increase in concentrations for Ca, Mg, Mn, Ba and Sr. Cation release rates decreased linearly as pH increased during mineral buffering. Experiment results suggest that carbonate minerals are the dominant contributor of changes in groundwater quality. Risk assessments of potential degradation of groundwater and monitoring strategies should focus on these fast-reacting minerals. Mobilization risk of Type II cations, however, may be self-mitigated because adsorption occurs when pH rebounds.     

Multi-spectral imaging of vegetation for detecting CO<Subscript>2</Subscript> leaking from underground

Practical geologic CO₂ sequestration will require long-term monitoring for detection of possible leakage back into the atmosphere. One potential monitoring method is multi-spectral imaging of vegetation reflectance to detect leakage through CO₂-induced plant stress. A multi-spectral imaging system was used to simultaneously record green, red, and near-infrared (NIR) images with a real-time reflectance calibration from a 3-m tall platform, viewing vegetation near shallow subsurface CO₂ releases during summers 2007 and 2008 at the Zero Emissions Research and Technology field site in Bozeman, Montana. Regression analysis of the band reflectances and the Normalized Difference Vegetation Index with time shows significant correlation with distance from the CO₂ well, indicating the viability of this method to monitor for CO₂ leakage. The 2007 data show rapid plant vigor degradation at high CO₂ levels next to the well and slight nourishment at lower, but above-background CO₂ concentrations. Results from the second year also show that the stress response of vegetation is strongly linked to the CO₂ sink–source relationship and vegetation density. The data also show short-term effects of rain and hail. The real-time calibrated imaging system successfully obtained data in an autonomous mode during all sky and daytime illumination conditions.     

CO<Subscript>2</Subscript> flux variation and its contribution area in the debris-covered area of Koxkar Glacier,Mt. Tianshan in China

During the formation and development of glacial meltwater runoff, hydrochemical erosion is abundant, especially the hydrolysis of K/Na feldspar and carbonates, which can consume H⁺ in the water, promote the formation of bicarbonate by dissolving atmospheric CO₂, and affect the regional carbon cycle. From July 21, 2015, to July 18, 2017, the CO₂ concentration and flux were observed by the eddy covariance (EC) method in the relatively flat and open moraine cover area of Koxkar Glacier in western Mt. Tianshan, China. We found that: (1) atmospheric CO₂ fluxes ranged from ??408.95 to 81.58 mmol m^?2 day^?1 (average ? 58.68 mmol m^?2 day^?1), suggesting that the study area is a significant carbon sink, (2) the CO₂ flux footprint contribution areas were primarily within 150 m of the EC station, averaging total contribution rates of 93.30%, 91.39%, and 90.17% of the CO₂ flux in the snow accumulation, snow melting, and glacial melting periods, respectively. Therefore, the contribution areas with significant influences on CO₂ flux observed at EC stations were concentrated, demonstrating that grassland CO₂ flux around the glaciers had little effect at the EC stations, (3) in the predominant wind direction, under stable daytime atmospheric stratification, the measurement of CO₂ flux, as interpreted by the Agroscope Reckenholz Tanikon footprint tool, was 79.09% ± 1.84% in the contribution area. This was slightly more than seen at night, but significantly lower than the average under unstable atmospheric stratification across the three periods of interest (89%). The average distance of the farthest point of the flux footprint under steady state atmospheric conditions was 202.61?±?69.33 m, markedly greater than that under non-steady state conditions (68.55?±?10.34 m). This also indicates that the CO₂ flux observed using EC was affected primarily by hydrochemical erosion reactions in the glacier area, (4) a good negative correlation was found between net glacier exchange (NGE) of CO₂ and air temperature on precipitation-free days. Strong ice and snow ablation could promote hydrochemical reactions of soluble substances in the debris area and accelerated sinking of atmospheric CO₂. Precipitation events might reduce snow and ice melting, driven by reduced regional temperatures. However, a connection between NGE and precipitation, when less than 8.8 mm per day, was not obvious. When precipitation was greater than 8.8 mm per day, NGE decreased with increasing precipitation, (5) graphically, the slope of NGE, related to daily runoff, followed a trend: snow melting period?>?snow accumulation period?>?early glacial ablation period?>?late glacier ablation period?>?dramatic glacier ablation period. The slope was relatively large during snow melting, likely because of CO₂ sinking caused by water–rock interactions. The chemical reaction during elution in the snow layer might also promote atmospheric CO₂ drawdown. At the same time, the damping effect of snow cover and the almost-closed glacier hydrographic channel inhibited the formation of regional runoff, possibly providing sufficient time for the chemical reaction, thus promoting further CO₂ drawdown.     

Effects of carbon dioxide variations in the unsaturated zone on water chemistry in a glacial-outwash aquifer

The research site at Otis Air Base, Cape Cod, Massachusetts, has been developed for hydrogeological and geochemical studies of sewage-effluent contaminated ground water since 1982. Research of hydrologic properties, transport, and chemical and biological processes is ongoing, but the origin of background water chemistry has not been determined.The principal geochemical process giving rise to the observed background water chemistry is CO₂-controlled hydrolysis of Na feldspar. Geochemical modeling demonstrated that CO₂ sources could vary over the project area. Analyses of unsaturated zone gases showed variations in CO₂ which were dependent on land use and vegetative cover in the area of groundwater recharge. Measurements of CO₂ in unsaturated-zone gases showed that concentrations of total inorganic C in recharge water should range from about 0.035 to 1.0 mmoles/L in the vicinity of Otis Air Base. Flux of CO₂ from the unsaturated zone varied for 4 principal land uses, ranging from 86 gC/m²/yr for low vegetated areas to 1630 gC/m²/yr for a golf course. Carbon dioxide flux from woodlands was 220 gC/m²/yr, lower than reported fluxes of 500 to 600 gC/m²/yr for woodlands in a similar climate. Carbon dioxide flux from grassy areas was 540 gC/m²/yr, higher than reported fluxes of 230 to 490 gC/m²/yr for grasslands in a similar climate.     

Continuous in situ measurements of volcanic gases with a diode-laser-based spectrometer: CO<Subscript>2</Subscript> and H<Subscript>2</Subscript>O concentration and soil degassing at Vulcano (Aeolian islands: Italy)

We report on a continuous-measurement campaign carried out in Vulcano (Aeolian islands, Sicily), devoted to the simultaneous monitoring of CO₂ and H₂O concentrations. The measurements were performed with an absorption spectrometer based on a semiconductor laser source emitting around a 2-μm wavelength. The emitted radiation was selectively absorbed by two molecular ro-vibrational transitions specific of the investigated species. Data for CO₂ and H₂O concentrations, and CO₂ soil diffusive flux using an accumulation chamber configuration, were collected at several interesting sampling points on the island (Porto Levante beach- PLB, Fossa Grande Crater – FOG- and Valley of Palizzi, PAL). CO₂/H₂O values, measured on the ground, are very similar (around 0.019 (± 0.006)) and comparable to the previous discrete detected values of 0.213 (Fumarole F5-La Fossa crater rim) and 0.012 (Fumarole VFS – Baia Levante beach) obtaid during the 1977–1993 heating phase of the crater fumaroles.     

Dynamics of CO<Subscript>2</Subscript> fluxes and concentrations during a shallow subsurface CO<Subscript>2</Subscript> release

A field facility located in Bozeman, Montana provides the opportunity to test methods to detect, locate, and quantify potential CO₂ leakage from geologic storage sites. From 9 July to 7 August 2008, 0.3 t CO₂ day⁻¹ were injected from a 100-m long, ~2.5-m deep horizontal well. Repeated measurements of soil CO₂ fluxes on a grid characterized the spatio-temporal evolution of the surface leakage signal and quantified the surface leakage rate. Infrared CO₂ concentration sensors installed in the soil at 30 cm depth at 0–10 m from the well and at 4 cm above the ground at 0 and 5 m from the well recorded surface breakthrough of CO₂ leakage and migration of CO₂ leakage through the soil. Temporal variations in CO₂ concentrations were correlated with atmospheric and soil temperature, wind speed, atmospheric pressure, rainfall, and CO₂ injection rate.     

Human footprints on greenhouse gas fluxes in cryogenic ecosystems

Various human footprints on the flux of biogenic greenhouse gases from permafrost-affected soils in Arctic and boreal domains in Russia are considered. Tendencies of significant growth or suppression of soil CO₂ fluxes change across types of human impact. Overall, the human impacts increase the mean value and variance of local soil CO₂ flux. Human footprint on methane exchange between soil and atmosphere is mediated by drainage. However, all the types of human impact suppress the sources and increase sinks of methane to the land ecosystems. N₂O flux grew under the considered types of human impact. Based on the results, we suggest that human footprint on soil greenhouse gases fluxes is comparable to the effect of climate change at an annual to decadal timescales.     

Volatile (H<Subscript>2</Subscript>O,CO<Subscript>2</Subscript>, Cl,S) budget of the Central American subduction zone

After more than a decade of multidisciplinary studies of the Central American subduction zone mainly in the framework of two large research programmes, the US MARGINS program and the German Collaborative Research Center SFB 574, we here review and interpret the data pertinent to quantify the cycling of mineral-bound volatiles (H₂O, CO₂, Cl, S) through this subduction system. For input-flux calculations, we divide the Middle America Trench into four segments differing in convergence rate and slab lithological profiles, use the latest evidence for mantle serpentinization of the Cocos slab approaching the trench, and for the first time explicitly include subduction erosion of forearc basement. Resulting input fluxes are 40–62 (53) Tg/Ma/m H₂O, 7.8–11.4 (9.3) Tg/Ma/m CO₂, 1.3–1.9 (1.6) Tg/Ma/m Cl, and 1.3–2.1 (1.6) Tg/Ma/m S (bracketed are mean values for entire trench length). Output by cold seeps on the forearc amounts to 0.625–1.25 Tg/Ma/m H₂O partly derived from the slab sediments as determined by geochemical analyses of fluids and carbonates. The major volatile output occurs at the Central American volcanic arc that is divided into ten arc segments by dextral strike-slip tectonics. Based on volcanic edifice and widespread tephra volumes as well as calculated parental magma masses needed to form observed evolved compositions, we determine long-term (10⁵ years) average magma and K₂O fluxes for each of the ten segments as 32–242 (106) Tg/Ma/m magma and 0.28–2.91 (1.38) Tg/Ma/m K₂O (bracketed are mean values for entire Central American volcanic arc length). Volatile/K₂O concentration ratios derived from melt inclusion analyses and petrologic modelling then allow to calculate volatile fluxes as 1.02–14.3 (6.2) Tg/Ma/m H₂O, 0.02–0.45 (0.17) Tg/Ma/m CO₂, and 0.07–0.34 (0.22) Tg/Ma/m Cl. The same approach yields long-term sulfur fluxes of 0.12–1.08 (0.54) Tg/Ma/m while present-day open-vent SO₂-flux monitoring yields 0.06–2.37 (0.83) Tg/Ma/m S. Input–output comparisons show that the arc water fluxes only account for up to 40 % of the input even if we include an “invisible” plutonic component constrained by crustal growth. With 20–30 % of the H₂O input transferred into the deeper mantle as suggested by petrologic modeling, there remains a deficiency of, say, 30–40 % in the water budget. At least some of this water is transferred into two upper-plate regions of low seismic velocity and electrical resistivity whose sizes vary along arc: one region widely envelopes the melt ascent paths from slab top to arc and the other extends obliquely from the slab below the forearc to below the arc. Whether these reservoirs are transient or steady remains unknown.     

A stochastic inversion workflow for monitoring the distribution of CO<Subscript>2</Subscript> injected into deep saline aquifers

We present a two-step stochastic inversion approach for monitoring the distribution of CO₂ injected into deep saline aquifers for the typical scenario of one single injection well and a database comprising a common suite of well logs as well as time-lapse vertical seismic profiling (VSP) data. In the first step, we compute several sets of stochastic models of the elastic properties using conventional sequential Gaussian co-simulations (SGCS) representing the considered reservoir before CO₂ injection. All realizations within a set of models are then iteratively combined using a modified gradual deformation algorithm aiming at reducing the mismatch between the observed and simulated VSP data. In the second step, these optimal static models then serve as input for a history matching approach using the same modified gradual deformation algorithm for minimizing the mismatch between the observed and simulated VSP data following the injection of CO₂. At each gradual deformation step, the injection and migration of CO₂ is simulated and the corresponding seismic traces are computed and compared with the observed ones. The proposed stochastic inversion approach has been tested for a realistic, and arguably particularly challenging, synthetic case study mimicking the geological environment of a potential CO₂ injection site in the Cambrian-Ordivician sedimentary sequence of the St. Lawrence platform in Southern Québec. The results demonstrate that the proposed two-step reservoir characterization approach is capable of adequately resolving and monitoring the distribution of the injected CO₂. This finds its expression in optimized models of P- and S-wave velocities, density, and porosity, which, compared to conventional stochastic reservoir models, exhibit a significantly improved structural similarity with regard to the corresponding reference models. The proposed approach is therefore expected to allow for an optimal injection forecast by using a quantitative assimilation of all available data from the appraisal stage of a CO₂ injection site.     

Relationship between soil CO<Subscript>2</Subscript> flux and volcanic tremor at Mt. Etna: Implications for magma dynamics

Large variations of the CO₂ flux through the soil were observed between November 2002 and January 2006 at Mt. Etna volcano. In many cases, the CO₂ flux was strongly influenced by changes in air temperature and atmospheric pressure. A new filtering method was then developed to remove the atmospheric influences on soil CO₂ flux and, at the same time, to highlight the variations strictly related to volcanic activity. Successively, the CO₂ corrected data were quantitatively compared with the spectral amplitude of the volcanic tremor by cross correlation function, cross-wavelet spectrum and wavelet coherence. These analyses suggested that the soil CO₂ flux variations preceded those of volcanic tremor by about 50 days. Given that volcanic tremor is linked to the shallow (a few kilometer) magma dynamics and soil CO₂ flux related to the deeper (~12 km b.s.l.) magma dynamics, the “delayed similarity” between the CO₂ flux and the volcanic tremor amplitude was used to assess the average speed in the magma uprising into the crust, as about 170–260 m per day. Finally, the large amount of CO₂ released before the onset of the 2004–2005 eruption indicated a deep ingression of new magma, which might have triggered such an eruption.