Uncertainty analysis of CO<Subscript>2</Subscript> flux components in subtropical evergreen coniferous plantation

We present an uncertainty analysis of ecological process parameters and CO₂ flux components (R _eco, NEE and gross ecosystem exchange (GEE)) derived from 3 years’ continuous eddy covariance measurements of CO₂ fluxes at subtropical evergreen coniferous plantation, Qianyanzhou of ChinaFlux. Daily-differencing approach was used to analyze the random error of CO₂ fluxes measurements and bootstrapping method was used to quantify the uncertainties of three CO₂ flux components. In addition, we evaluated different models and optimization methods in influencing estimation of key parameters and CO₂ flux components. The results show that: (1) Random flux error more closely follows a double-exponential (Laplace), rather than a normal (Gaussian) distribution. (2) Different optimization methods result in different estimates of model parameters. Uncertainties of parameters estimated by the maximum likelihood estimation (MLE) are lower than those derived from ordinary least square method (OLS). (3) The differences between simulated R_eco, NEE and GEE derived from MLE and those derived from OLS are 12.18% (176 g C·m⁻²·a⁻¹), 34.33% (79 g C·m⁻²·a⁻¹) and 5.4% (92 g C·m⁻²·a⁻¹). However, for a given parameter optimization method, a temperature-dependent model (T_model) and the models derived from a temperature and water-dependent model (TW_model) are 1.31% (17.8 g C·m⁻²·a⁻¹), 2.1% (5.7 g C·m⁻²·a⁻¹), and 0.26% (4.3 g C·m⁻²·a⁻¹), respectively, which suggested that the optimization methods are more important than the ecological models in influencing uncertainty in estimated carbon fluxes. (4) The relative uncertainty of CO₂ flux derived from OLS is higher than that from MLE, and the uncertainty is related to timescale, that is, the larger the timescale, the smaller the uncertainty. The relative uncertainties of R_eco, NEE and GEE are 4%−8%, 7%−22% and 2%−4% respectively at annual timescale. Supported by the National Natural Science Foundation of China (Grant No. 30570347), Innovative Research International Partnership Project of the Chinese Academy of Sciences (Grant No. CXTD-Z2005-1) and National Basic Research Program of China (Grant No. 2002CB412502)     

Dynamics of carbon dioxide emissions,crystallization, and magma ascent: hypotheses,theory, and applications to volcano monitoring at Mount St. Helens

 Measurements of CO₂ fluxes from open-vent volcanos are rare, yet may offer special capabilities for monitoring volcanos and forecasting activity. The measured fluxes of CO₂ and SO₂ from Mount St. Helens decreased from July through November 1980, but the record includes variations of CO₂/SO₂ in the emitted gas and episodes of greatly increased fluxes of CO₂. We propose that the CO₂ flux variations reflect two gas components: (a) a component whose flux decreased in proportion to 1/ √t with a CO₂/SO₂ mass ratio of 1.7, and (b) a residual flux of CO₂ consisting of short-lived, large peaks with a CO₂/SO₂ mass ratio of 15. We propose two hypotheses: (a) the 1/ √t dependence was generated by crystallization in a deep magma body at rates governed by diffusion-limited heat transfer, and (b) the gas component with the higher CO₂/SO₂ was released from ascending magma, which replenished the same magma body. The separation of the total CO₂ flux into contributions from known processes permits quantitative inferences about the replenishment and crystallization rates of open-system magma bodies beneath volcanos. The flux separations obtained by using two gas sources with distinct CO₂/SO₂ ratios and a peak minus background approach to obtain the CO₂ contributions from an intermittent source and a continuously emitting source are similar. The flux separation results support the hypothesis that the second component was generated by episodic magma ascent and replenishment of the magma body. The diffusion-limited crystallization hypothesis is supported by the decay of minimum CO₂ and SO₂ fluxes with 1/ √t after 1 July 1980. We infer that the magma body at Mount St. Helens was replenished at an average rate (2.8×10⁶m³d^–1) which varied by less than 5% during July, August, and September 1980. The magma body volume (2.4–3.0 km³) in early 1982 was estimated by integrating a crystallization rate function inferred from CO₂ fluxes to maximum times (20±4 years) estimated from the increase of sample crystallinity with time. These new volcanic gas flux separation methods and the existence of relations among the CO₂ flux, crystallization rates, and magma body replenishment rates yield new information about the dynamics of an open-vent, replenished magma body. Received: 15 February 1995 / Accepted: 30 March 1996     

Geochemical evidence of the renewal of volcanic activity inferred from CO<Subscript>2</Subscript> soil and SO<Subscript>2</Subscript> plume fluxes: the 2007 Stromboli eruption (Italy)

On 27 February 2007, a new eruption occurred on Stromboli which lasted until 2 April. It was characterized by effusive activity on the Sciara del Fuoco and by a paroxysmal event (15 March). This crisis represented an opportunity for us to refine the model that had been developed previously (2002–2003 eruption) and to improve our understanding of the relationship between the magmatic dynamics of the volcano and the geochemical variations in the fluids. In particular, the evaluation of the dynamic equilibrium between the volatiles (CO₂ and SO₂) released from the magma and the corresponding fluids discharged from the summit area allowed us to evaluate the level of criticality of the volcanic activity. One of the major accomplishments of this study is a 4-year database of summit soil CO₂ flux on the basis of which we define the thresholds (low–medium–high) for this parameter that are empirically based on the natural volcanological evolution of Stromboli. The SO₂ fluxes of the degassing plume and the CO₂ fluxes emitted from the soil at Pizzo Sopra la Fossa are also presented. It is noteworthy that geochemical signals of volcanic unrest have been clearly identified before, during and after the effusive activity. These signals were found almost simultaneously in the degassing plume (SO₂ flux) and in soil degassing (CO₂ flux) at the summit, although the two degassing processes are shown to be clearly different. The interpretation of the results will be useful for future volcanic surveillance at Stromboli.     

First <Superscript>13</Superscript>C/<Superscript>12</Superscript>C isotopic characterisation of volcanic plume CO<Subscript>2</Subscript>

We describe analytical details and uncertainty evaluation of a simple technique for the measurement of the carbon isotopic composition of CO₂ in volcanic plumes. Data collected at Solfatara and Vulcano, where plumes are fed by fumaroles which are accessible for direct sampling, were first used to validate the technique. For both volcanoes, the plume-derived carbon isotopic compositions are in good agreement with the fumarolic compositions, thus providing confidence on the method, and allowing its application at volcanoes where the volcanic component is inaccessible to direct sampling. As a notable example, we applied the same method to Mount Etna where we derived a δ¹³C of volcanic CO₂ between −0.9 ± 0.27‰ and −1.41 ± 0.27‰ (Bocca Nuova and Voragine craters). The comparison of our measurements to data reported in previous work highlights a temporal trend of systematic increase of δ¹³C values of Etna CO₂ from ~ −4‰, in the 1970’s and the 1980’s, to ~ −1‰ at the present time (2009). This shift toward more positive δ¹³C values matches a concurrent change in magma composition and an increase in the eruption frequency and energy. We discuss such variations in terms of two possible processes: magma carbonate assimilation and carbon isotopic fractionation due to magma degassing along the Etna plumbing system. Finally, our results highlight potential of systematic measurements of the carbon isotopic composition of the CO₂ emitted by volcanic plumes for a better understanding of volcanic processes and for improved surveillance of volcanic activity.     

CO<Subscript>2</Subscript> and He degassing at El Chichón volcano,Chiapas, Mexico: gas flux,origin and relationship with local and regional tectonics

During 2007–2008, three CO₂ flux surveys were performed on El Chichón volcanic lake, Chiapas, Mexico, with an additional survey in April 2008 covering the entire crater floor (including the lake). The mean CO₂ flux calculated by sequential Gaussian simulation from the lake was 1,190 (March 2007), 730 (December 2007) and 1,134 g m⁻² day⁻¹ (April 2008) with total emission rates of 164 ± 9.5 (March 2007), 59 ± 2.5 (December 2007) and 109 ± 6.6 t day⁻¹ (April 2008). The mean CO₂ flux estimated from the entire crater floor area was 1,102 g m⁻² day⁻¹ for April 2008 with a total emission rate of 144 ± 5.9 t day⁻¹. Significant change in CO₂ flux was not detected during the period of survey, and the mapping of the CO₂ flux highlighted lineaments reflecting the main local and regional tectonic patterns. The ³He/⁴He ratio (as high as 8.1 R _A) for gases in the El Chichón crater is generally higher than those observed at the neighbouring Transmexican Volcanic Belt and the Central American Volcanic Arc. The CO₂/³He ratios for the high ³He/⁴He gases tend to have the MORB-like values (1.41 × 10⁹), and the CO₂/³He ratios for the lower ³He/⁴He gases fall within the range for the arc-type gases. The high ³He/⁴He ratios, the MORB-like CO₂/³He ratios for the high ³He/⁴He gases and high proportion of MORB-CO₂ (M = 25 ±15%) at El Chichón indicate a greater depth for the generation of magma when compared to typical arc volcanoes.     

Time-dependent CO<Subscript>2</Subscript> variations in Lake Albano associated with seismic activity

Lake Albano (Alban Hills volcanic complex, Central Italy) is located in a densely populated area near Rome. The deep lake waters have significant dissolved CO₂ concentrations, probably related to sub-lacustrine fluid discharges fed by a pressurized CO₂-rich reservoir. The analytical results of geochemical surveys carried out in 1989–2010 highlight the episodes of CO₂ removal from the lake. The total mass of dissolved CO₂ decreased from ∼5.8 × 10⁷ kg in 1989 to ∼0.5 × 10⁷ kg in 2010, following an exponential decreasing trend. Calculated values of both dissolved inorganic carbon and CO₂ concentrations along the vertical profile of the lake indicate that this decrease is caused by CO₂ release from the epilimnion, at depth <9 m, combined with (1) water circulation at depth <95 m and (2) CO₂ diffusion from the deeper lake layers. According to this model, Lake Albano was affected by a large CO₂ input that coincided with the last important seismic swarm at Alban Hills in 1989, suggesting an intimate relationship between the addition of deep-originated CO₂ to the lake and seismic activity. In the case of a CO₂ degassing event of an order of magnitude larger than the one that occurred in 1989, the deepest part of Lake Albano would become CO₂-saturated, resulting in conditions compatible with the occurrence of a gas outburst. These results reinforce the idea that a sudden CO₂ input into the lake may cause the release of a dense gas cloud, presently representing the major volcanic threat for this densely populated area.     

Soil organic carbon storage and soil CO<Subscript>2</Subscript> flux in the alpine meadow ecosystem

High-resolution sampling, measurements of organic carbon contents and ¹⁴C signatures of selected four soil profiles in the Haibei Station situated on the northeast Tibetan Plateau, and application of ¹⁴C tracing technology were conducted in an attempt to investigate the turnover times of soil organic carbon and the soil-CO₂ flux in the alpine meadow ecosystem. The results show that the organic carbon stored in the soils varies from 22.12×10⁴ kg C hm⁻² to 30.75×10⁴ kg C hm⁻² in the alpine meadow ecosystems, with an average of 26.86×10⁴ kg C hm⁻². Turnover times of organic carbon pools increase with depth from 45 a to 73 a in the surface soil horizon to hundreds of years or millennia or even longer at the deep soil horizons in the alpine meadow ecosystems. The soil-CO₂ flux ranges from 103.24 g C m⁻² a⁻¹ to 254.93 gC m⁻² a⁻¹, with an average of 191.23 g C m⁻² a⁻¹. The CO₂ efflux produced from microbial decomposition of organic matter varies from 73.3 g C m⁻² a⁻¹ to 181 g C m⁻² a⁻¹. More than 30% of total soil organic carbon resides in the active carbon pool and 72.8%281.23% of total CO₂ emitted from organic matter decomposition results from the topsoil horizon (from 0 cm to 10 cm) for the Kobresia meadow. Responding to global warming, the storage, volume of flow and fate of the soil organic carbon in the alpine meadow ecosystem of the Tibetan Plateau will be changed, which needs further research. Supported by the National Natural Science Foundation of China (Grant Nos. 40231015, 40471120 and 40473002) and the Guangdong Provincial Natural Science Foundation of China (Grant No. 06300102)     

Sulphur emissions to the stratosphere from explosive volcanic eruptions

 Two methods were used to quantify the flux of volcanic sulphur (as the equivalent mass of SO₂) to the stratosphere over different timescales during the Holocene. A combination of satellite-based measurements of sulphur yields from recent explosive volcanic eruptions with an appropriate rate of explosive volcanism for the past 200 years constrains the medium-term (∼10²years) flux of volcanic sulphur to the stratosphere to be ∼1 Mt a^–1, with lower and upper bounds of 0.3 and 3 Mt a^–1. The short-term (∼10- to 20-year) flux due to small magnitude (10¹⁰–10¹²kg) eruptions is of the order of 0.4 Mt a^–1. At any time the instantaneous levels of sulphur in the stratosphere are dominated by the most recent (0–3 years) volcanic events. The flux calculations do not attempt to address this very short timescale variability. Although there are significant errors associated with the raw sulphur emission data on which this analysis is based, the approach presented is general and may be readily modified as the quantity and quality of the data improve. Data from a Greenland ice core support these conclusions. Integration of the sulphate signals from presumed volcanic sources recorded in the GISP2 core provides a minimum estimate of the 10³–year volcanic SO₂ flux to the stratosphere of 0.5–1 Mt a^–1 over the past 9000 years. The short-term flux calculations do not account for the impact of rare, large events. The ice-core record does not fully account for the contribution from small, frequent events. Received: 27 September 1995 / Accepted: 13 December 1995     

Structure and CO<Subscript>2</Subscript> budget of Merapi volcano during inter-eruptive periods

Soil temperature and gas (CO₂ concentration and flux) have been investigated at Merapi volcano (Indonesia) during two inter-eruptive periods (2002 and 2007). Precise imaging of the summit crater and the spatial pattern of diffuse degassing along a gas traverse on the southern slope are interpreted in terms of summit structure and major caldera organization. The summit area is characterized by decreasing CO₂ concentrations with distance from the 1932 crater rim, down to atmospheric levels at the base of the terminal cone. Similar patterns are measured on any transect down the slopes of the cone. The spatial distribution of soil gas anomalies suggests that soil degassing is controlled by structures identified as concentric historical caldera rims (1932, 1872, and 1768), which have undergone severe hydrothermal self-sealing processes that dramatically lower the permeability and porosity of soils. Temperature and CO₂ flux measurements in soils near the dome display heterogeneous distributions which are consistent with a fracture network identified by previous geophysical studies. These data support the idea that the summit is made of isolated and mobile blocks, whose boundaries are either sealed by depositional processes or used as pathways for significant soil degassing. Within this context, self-sealing both prevents long-distance soil degassing and controls heat and volatile transfers near the dome. A rough estimate of the CO₂ output through soils near the dome is 200–230 t day⁻¹, i.e. 50% of the estimated total gas output from the volcano summit during these quiescent periods. On Merapi’s southern slope, a 2,500 m long CO₂ traverse shows low-amplitude anomalies that fit well with a recently observed electromagnetic anomaly, consistent with a faulted structure related to an ancient avalanche caldera rim. Sub-surface soil permeability is the key parameter that controls the transfer of heat and volatiles within the volcano, allowing its major tectonic architecture to be revealed by soil gas and soil temperature surveys.     

Diffuse CO<Subscript>2</Subscript> degassing at Vesuvio,Italy

At Vesuvio, a significant fraction of the rising hydrothermal–volcanic fluids is subjected to a condensation and separation process producing a CO₂–rich gas phase, mainly expulsed through soil diffuse degassing from well defined areas called diffuse degassing structures (DDS), and a liquid phase that flows towards the outer part of the volcanic cone. A large amount of thermal energy is associated with the steam condensation process and subsequent cooling of the liquid phase. The total amount of volcanic–hydrothermal CO₂ discharged through diffuse degassing has been computed through a sequential Gaussian simulation (sGs) approach based on several hundred accumulation chamber measurements and, at the time of the survey, amounted to 151 t d^–1. The steam associated with the CO₂ output, computed assuming that the original H₂O/CO₂ ratio of hydrothermal fluids is preserved in fumarolic effluents, is 553 t d^–1, and the energy produced by the steam condensation and cooling of the liquid phase is 1.47×10¹² J d^–1 (17 MW). The location of the CO₂ and temperature anomalies show that most of the gas is discharged from the inner part of the crater and suggests that crater morphology and local stratigraphy exert strong control on CO₂ degassing and subsurface steam condensation. The amounts of gas and energy released by Vesuvio are comparable to those released by other volcanic degassing areas of the world and their estimates, through periodic surveys of soil CO₂ flux, can constitute a useful tool to monitor volcanic activity.Editorial responsibility: H. Shinohara     

Diffuse emission of CO2 from the Fossa crater,Vulcano Island (Italy)

 Numerous measurements of CO₂ degassing from the soil, carried out with the accumulation chamber method, indicate that in the period April–July 1995 the upper part of the Fossa cone released a total output of 200 t d^–1 of CO₂, which corresponds to approximately 1000 t d^–1 of steam. These large amounts of fluids are of the same order of magnitude as those released by the high temperature fumarolic field located inside the crater. The spatial distribution of soil gas fluxes shows that the main structures releasing CO₂ are the inner slopes of the crater and a NW–SE line, located NE of the crater rim, which correspond to the main direction of Vulcano Island active faults. The comparison of the φCO₂ maps with the soil temperature distribution, derived from both direct measurements and airborne infrared images, indicates the occurrence of extensive condensation of fumarolic steam within the upper part of the Fossa cone, whose total amount is comparable to the rainfall budget. Part of the condensate which originates from this process contributes to the recharge of the phreatic aquifer of Porto Plain, modifying the chemical and isotopic composition of the groundwater. Received: 1 September 1995 / Accepted: 8 January 1996     

Soil CO2 flux baseline in an urban monogenetic volcanic field: the Auckland Volcanic Field, New Zealand

The Auckland Volcanic Field (AVF) is a dormant monogenetic basaltic field located in Auckland, New Zealand. Though soil gas CO₂ fluxes are routinely used to monitor volcanic regions, there have been no published studies of soil CO₂ flux or soil gas CO₂ concentrations in the AVF to date or many other monogenetic fields worldwide. We measured soil gas CO₂ fluxes and soil gas CO₂ concentrations in 2010 and 2012 in varying settings, seasons, and times of day to establish a baseline soil CO₂ flux and to determine the major sources of and controlling influences on Auckland's soil CO₂ flux. Soil CO₂ flux measurements varied from 0 to 203 g m^?2 day^?1, with an average of 27.1 g m^?2 day^?1. Higher fluxes were attributed to varying land use properties (e.g., landfill). Using a graphical statistical approach, two populations of CO₂ fluxes were identified. Isotope analyses of δ¹³CO₂ confirmed that the source of CO₂ in the AVF is biogenic with no volcanic component. These data may be used to assist with eruption forecasting in the event of precursory activity in the AVF, and highlight the importance of knowing land use history when assessing soil gas CO₂ fluxes in urban environments.     

Contrasting hydrothermal activity at Sierra Negra and Alcedo volcanoes, Galapagos Archipelago, Ecuador

2 and approximately 85% SO₂ of the total sulfur gas. Relative amounts of He, Ar, and N₂ show a distinct hot-spot signature ( ). The δ¹³C–CO₂ is approximately −3.6‰ and δ³⁴S_T is approximately +3.3‰. The δD/δ¹⁸O of fumarole H₂O indicates steam separation from local meteoric waters whose estimated minimum mean residence time from ³H analyses is ≤40 years. Fumarolic activity at Alcedo is controlled by a caldera-margin fault containing at least seven hydrothermal explosion craters, and by an intracaldera rhyolite vent. Two explosion craters which formed in 1993–1994 produce approximately 15 m³/s of steam, yet discharge temperatures are ≤97°C. Water content of the total gas is 95–97 mol.%, noncondensible gas is 92–98 mol.% CO₂, and sulfur gas is dominated by H₂S. Relative amounts of He, Ar, and N₂ show extensive mixing between hot spot and air or air-saturated meteoric water components but the average . The δ¹³C–CO₂ is approximately −3.5‰ and δ³⁴S_T is approximately −0.8‰. The δD/δ¹⁸O of fumarole steam indicates separation from a homogeneous reservoir that is enriched 3–5‰ in ¹⁸O compared with local meteoric water. ³H indicates that this reservoir water has a maximum mean residence time of approximately 400 years and empirical gas geothermometry indicates a reservoir temperature of 260–320°C. The intracaldera hydrothermal reservoir in Alcedo is probably capable of producing up to 150 MW; however, environmental concerns as well as lack of infrastructure and power users will limit the development of this resource. Received: 19 April 1999 / Accepted: 23 October 1999     

Seasonal variations and mechanism for environmental control of NEE of CO<Subscript>2</Subscript> concerning the <Emphasis Type="Italic">Potentilla fruticosa</Emphasis> in alpine shrub meadow of Qinghai-Tibet Plateau

The study by the eddy covariance technique in the alpine shrub meadow of the Qinghai-Tibet Plateau in 2003 and 2004 showed that the net ecosystem carbon dioxide exchange (NEE) exhibited noticeable diurnal and annual variations, with more distinct daily changes during the warmer seasons. The CO₂ emission of the shrub ecosystem culminated in April and September while the CO₂ absorption capacity reached a maximum in July and August. The absorbed carbon dioxide during the two consecutive years was 231.4 and 274.8 g CO₂·m^?2 respectively, yielding an average of 253.1 gCO₂·m^?2 per year: that accounts for a large proportion of absorbed CO₂ in the region. Obviously, the diurnal carbon flux was negatively related to temperature, radiation and other atmospheric factors. Still, minute discrepancies in kurtosis and duration of carbon emission/absorption were detected between 2003 and 2004. It was found that the CO₂ flux in the daytime was similarly affected by photosynthetic photon flux density in both years. Temperature appears to be the most important determinant of CO₂ flux: specifically, the high temperature during the plant growing season inhibits the carbon absorption capacity. One potential explanation is that soil respiration is enhanced under such condition. Analysis of biomass revealed that the annual net carbon fixed capacity of aboveground and belowground biomass was 544.0 in 2003 and 559.4 g C·m^?2 in 2004, which coincided with the NEE absorption capacity (63.1 g C·m^?2 in 2003 and 74.9 g C·m^?2 in 2004) in the corresponding plant growing season.     

Simulation of crop growth and energy and carbon dioxide fluxes at different time steps from hourly to daily

Understanding the exchange processes of energy and carbon dioxide (CO₂) in the soil–vegetation–atmosphere system is important for assessing the role of the terrestrial ecosystem in the global water and carbon cycle and in climate change. We present a soil–vegetation–atmosphere integrated model (ChinaAgrosys) for simulating energy, water and CO₂ fluxes, crop growth and development, with ample supply of nutrients and in the absence of pests, diseases and weed damage. Furthermore, we test the hypotheses of whether there is any significant difference between simulations over different time steps. CO₂, water and heat fluxes were estimated by the improving parameterization method of the coupled photosynthesis–stomatal conductance–transpiration model. Soil water evaporation and plant transpiration were calculated using a multilayer water and heat‐transfer model. Field experiments were conducted in the Yucheng Integrated Agricultural Experimental Station on the North China Plain. Daily weather and crop growth variables were observed during 1998–2001, and hourly weather variables and water and heat fluxes were measured using the eddy covariance method during 2002–2003. The results showed that the model could effectively simulate diurnal and seasonal changes of net radiation, sensible and latent heat flux, soil heat flux and CO₂ fluxes. The processes of evapotranspiration, soil temperature and leaf area index agree well with the measured values. Midday depression of canopy photosynthesis could be simulated by assessing the diurnal change in canopy water potential. Moreover, the comparisons of simulated daily evapotranspiration and net ecosystem exchange (NEE) under different time steps indicated that time steps used by a model affect the simulated results. There is no significant difference between simulated evapotranspiration using the model under different time steps. However, simulated NEE produces large differences in the response to different time steps. Therefore, the accurate calculation of average absorbed photosynthetic active radiation is important for the scaling of the model from hourly steps to daily steps in simulating energy and CO₂ flux exchanges between winter wheat and the atmosphere. Copyright © 2007 John Wiley & Sons, Ltd.     

Modelling the dynamics and thermodynamics of volcanic degassing

 The rates of passive degassing from volcanoes are investigated by modelling the convective overturn of dense degassed and less dense gas-rich magmas in a vertical conduit linking a shallow degassing zone with a deep magma chamber. Laboratory experiments are used to constrain our theoretical model of the overturn rate and to elaborate on the model of this process presented by Kazahaya et al. (1994). We also introduce the effects of a CO₂–saturated deep chamber and adiabatic cooling of ascending magma. We find that overturn occurs by concentric flow of the magmas along the conduit, although the details of the flow depend on the magmas' viscosity ratio. Where convective overturn limits the supply of gas-rich magma, then the gas emission rate is proportional to the flow rate of the overturning magmas (proportional to the density difference driving convection, the conduit radius to the fourth power, and inversely proportional to the degassed magma viscosity) and the mass fraction of water that is degassed. Efficient degassing enhances the density difference but increases the magma viscosity, and this dampens convection. Two degassing volcanoes were modelled. At Stromboli, assuming a 2 km deep, 30% crystalline basaltic chamber, containing 0.5 wt.% dissolved water, the ∼700 kg s^–1 magmatic water flux can be modelled with a 4–10 m radius conduit, degassing 20–100% of the available water and all of the 1 to 4 vol.% CO₂ chamber gas. At Mount St. Helens in June 1980, assuming a 7 km deep, 39% crystalline dacitic chamber, containing 4.6 wt.% dissolved water, the ∼500 kg s^–1 magmatic water flux can be modelled with a 22–60 m radius conduit, degassing ∼2–90% of the available water and all of the 0.1 to 3 vol.% CO₂ chamber gas. The range of these results is consistent with previous models and observations. Convection driven by degassing provides a plausible mechanism for transferring volatiles from deep magma chambers to the atmosphere, and it can explain the gas fluxes measured at many persistently active volcanoes. Received: 26 September 1997 / Accepted: 11 July 1998     

三峡水库澎溪河消落区土-气界面CO2和CH4通量初探

水库近岸湿地(消落区)温室气体(CO₂、CH₄)产汇是水库温室气体效应问题的重要组成部分.本文以三峡水库支流澎溪河的白家溪、养鹿两处大面积消落区为研究对象,于2010年6 9月水库低水位运行期间,对近岸消落区土-气界面CO₂、CH₄通量进行监测.白家溪消落区土-气界面CO₂通量均值为12.38±2.42 mmol/(m²·h);CH₄通量均值为0.0112±0.0064 mmol/(m²·h).养鹿消落区CO₂、CH₄通量均值分别为10.54±5.17、0.14±0.16 mmol/(m²·h).总体上,6 9月土-气界面CO₂通量呈增加趋势,而CH₄通量水平呈现显著的递减趋势.消落区土地出露后植被恢复,在一定程度上促进了土壤有机质含量的增加,使得6 9月CO₂释放通量的总体趋势有所增加.消落区退耕后,其甲烷氧化菌的活性得到恢复,加之在土地出露曝晒过程中土壤透气性增强,使得消落区土壤对大气中CH₄吸收氧化潜势增强.尽管如此,仍需进一步的研究以明晰消落区土-气界面CO₂、CH₄产汇的主要影响因素.     

Surge in sulphur and halogen degassing from Ambrym volcano,Vanuatu

Volcanoes provide important contributions to atmospheric budgets of SO₂ and reactive halogens, which play significant roles in atmospheric oxidative capacity and radiation. However, the global source strengths of volcanic emissions remain poorly constrained. These uncertainties are highlighted here by the first measurements of gas emission rates from Ambrym volcano, Vanuatu. Our initial airborne ultraviolet spectroscopic measurements made in January 2005 indicate fluxes of 18–270 kg s^-1 of SO₂, and 62–110 g s^-1 of BrO, into the atmosphere, placing Ambrym amongst the largest known contemporary point sources of both these species on Earth. We also estimate high Cl and F fluxes of ~8–14 and ~27–50 kg s^-1, respectively, for this period. Further observations using both airborne and spaceborne remote sensing reveal a fluctuating SO₂ output between 2004 and 2008, with a surge in the first half of 2005, and underline the substantial contribution that a single passively degassing volcano can make to the atmospheric budget of sulfur and halogens.     

基于同步观测信息利用的高寒湖泊湍流通量数据插补方法

源区划分和质量过滤提高湖面涡动相关通量数据可靠性的同时,却降低了通量时间序列的连续性.为此,本文基于TensorFlow机器学习框架构建了一种超宽人工神经网络(ANN)模型.在选择输入ANN模型的特征变量信息时,我们采取了尽可能获取湍流输送过程中热力、动力学同步观测背景强迫信息的原则.通过ANN模型模拟通量的插补,本文实现了通量时间序列连续性的优化,插补后的羊卓雍错湖面通量数据的时间覆盖率从不足0.40提升至超过0.98.基于10次折叠交叉验证的ANN模型通量模拟性能检验则表明,各个检验组之间ANN模型的模拟误差波动较小,这显示出了较好的稳健性.具体地讲,感热通量、潜热通量和水汽通量原始观测平均值分别约为18.8 W/m~2、81.5 W/m~2和1.84 mmol/(s·m~2),10组交叉验证的插补感热通量、潜热通量和水汽通量平均绝对误差分别为5.4 W/m~2、15.7 W/m~2和0.35 mmol/(s·m~2).这表明本文所探索的ANN建模结构和同步观测变量筛选原则可更充分地利用观测点局地同步观测信息估算通量强度,有效地优化湍流通量数据的时间连续性,从而提升通量数据的可分析性.     

Comparative soil CO2 flux measurements and geostatistical estimation methods on Masaya volcano, Nicaragua

We present a comparative study of soil CO₂ flux () measured by five groups (Groups 1–5) at the IAVCEI-CCVG Eighth Workshop on Volcanic Gases on Masaya volcano, Nicaragua. Groups 1–5 measured using the accumulation chamber method at 5-m spacing within a 900 m² grid during a morning (AM) period. These measurements were repeated by Groups 1–3 during an afternoon (PM) period. Measured ranged from 218 to 14,719 g m⁻² day⁻¹. The variability of the five measurements made at each grid point ranged from ±5 to 167%. However, the arithmetic means of fluxes measured over the entire grid and associated total CO₂ emission rate estimates varied between groups by only ±22%. All three groups that made PM measurements reported an 8–19% increase in total emissions over the AM results. Based on a comparison of measurements made during AM and PM times, we argue that this change is due in large part to natural temporal variability of gas flow, rather than to measurement error. In order to estimate the mean and associated CO₂ emission rate of one data set and to map the spatial distribution, we compared six geostatistical methods: arithmetic and minimum variance unbiased estimator means of uninterpolated data, and arithmetic means of data interpolated by the multiquadric radial basis function, ordinary kriging, multi-Gaussian kriging, and sequential Gaussian simulation methods. While the total CO₂ emission rates estimated using the different techniques only varied by ±4.4%, the maps showed important differences. We suggest that the sequential Gaussian simulation method yields the most realistic representation of the spatial distribution of , but a variety of geostatistical methods are appropriate to estimate the total CO₂ emission rate from a study area, which is a primary goal in volcano monitoring research.Editorial responsibility: H Shinohara