HR: 0800h
AN: V11B-1425 [Abstracts]
TI: Comparative investigation of soil CO$_{2}$ flux measurements and geostatistical estimation methods on
Masaya volcano, Nicaragua
AU: * Lewicki, J L
EM: jllewicki@lbl.gov
AF: Lawrence Berkeley National Lab, Earth Sciences Division, 1 Cyclotron Road, Berkeley, CA 94720
United States
AU: Bergfeld, D
EM: debergfel@usgs.gov
AF: U.S. Geological Survey, 345 Middlefield Road, Menlo Park, CA 94025
United States
AU: Cardellini, C
EM: geochem@unipg.it
AF: Universita di Perugia, Dipartimento di Scienze della Terra, Perugia, 00000
Italy
AU: Chiodini, G
EM: chiod@ov.ingv.it
AF: Osservatorio Vesuviano, Istituto Nazionale di Geofisica e Vulcanologia, via Diocleziano 328, Naples,
80124
Italy
AU: Granieri, D
EM: granieri@ov.ingv.it
AF: Osservatorio Vesuviano, Istituto Nazionale di Geofisica e Vulcanologia, via Diocleziano 328, Naples,
80124
Italy
AU: Varley, N
EM: eruption@softhome.net
AF: Universidad de Colima, Facultad de Ciencias, Colima, 28010
Mexico
AU: Werner, C
EM: c.werner@gns.cri.nz
AF: Institute of Geological and Nuclear Sciences, Wairakei Research Centre, Taupo, 2730
New Zealand
AB:
Measurements of soil CO$_{2}$ flux (F$_{CO2}$) and its natural spatial and temporal variability in volcanic and hydrothermal
environments are important for volcano monitoring, geothermal exploration, delineation of fault and fracture zones, and
estimation of the contribution of CO$_{2}$ from volcanic and hydrothermal sources to the global carbon cycle. However, the
choice of measurement and geostatistical methodologies may affect individual F$_{CO2}$ measurements and characterization of
their natural spatial and temporal variability, the total CO$_{2}$ emission rate estimated, and the ability to assess the
uncertainty associated with this estimate. We present a comparative study of F$_{CO2}$ measured by five research groups
(Groups 1-5) at the 2003 IAVCEI-CCVG Eighth Workshop on Volcanic Gases on Masaya volcano, Nicaragua. Groups 1-5 measured
F$_{CO2}$ using the accumulation chamber method at 5-m spacing within a 900~m$^{2}$ grid during a morning (AM) period. These
measurements were repeated by Groups 1-3 during an afternoon (PM) period. Measured F$_{CO2}$ ranged from 218 to
14,719~g~m$^{-2}$d^{-1}$. The variability of the five measurements made at each grid point ranged from $\pm$5 to 167%;
however, the arithmetic means of fluxes measured over the entire grid and associated total CO$_{2}$ emission rate estimates
varied between groups by only $\pm$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, this change is likely due
in large part to natural temporal variability of gas flow, rather than to measurement error. To estimate the mean and
associated CO$_{2}$ emission rate of one data set and to map the spatial F$_{CO2}$ 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$_{2}$ emission rates estimated using the different techniques only varied by
$\pm$1.1%, the F$_{CO2}$ maps showed important differences. We suggest that the sequential Gaussian simulation method
yields the most realistic representation of the spatial distribution of F$_{CO2}$, but a variety of geostatistical methods
are appropriate to estimate the total CO$_{2}$ emission rate from a study area, which is a primary goal in volcano monitoring
research. Acknowledgement: Part of this work was completed at Lawrence Berkeley National Laboratory, under U.S. Department
of Energy Contract No. DE-AC03765F00098.
DE: 8419 Eruption monitoring (7280)
DE: 8424 Hydrothermal systems (8135)
DE: 8494 Instruments and techniques
DE: 1875 Unsaturated zone
SC: Volcanology, Geochemistry, Petrology [V]
MN: 2004 AGU Fall Meeting