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     

The Dorn gold deposit in northern Victoria Land,Antarctica: Structure,hydrothermal alteration,and implications for the Gondwana Pacific margin

The Dorn gold deposit in northern Victoria Land, Antarctica is a fault related gold-only deposit and it represents the first described occurrence of gold mineralization in Paleozoic terranes of the Antarctic continent. The deposit is hosted in lower greenschist facies Middle Cambrian metavolcanic and metasedimentary rocks of the Bowers terrane. Gold-bearing veins are located in a brittle–ductile reverse high-strain zone, which was produced by transpressional deformation that overprints the structures of the Cambrian–Ordovician Ross Orogeny. The vein system is surrounded by a hydrothermal alteration zone that is as much as 300-m-wide, where the host rocks are partially to completely transformed into Fe–Mg carbonate-rich rocks with different degrees of replacement of the original mineralogy and texture.The type of host rock, the temperature estimates for mineralization from 290–320 °C, the quartz dominant vein system with sulfides and Fe-rich carbonates, and the controlling structures linked to a convergent margin tectonic setting together suggest that this mineralized vein system can be classified as an orogenic gold deposit. Close analogies are found with deposits of the Stawell zone in western Victoria, which is consistent with the correlation between units and hydrothermal events in northern Victoria Land and southeastern Australia.     

Multi-stage melt–rock interaction in the Mt. Maggiore (Corsica, France) ophiolitic peridotites: microstructural and geochemical evidence

Spinel and plagioclase peridotites from the Mt.Maggiore (Corsica, France) ophiolitic massif record a composite asthenosphere–lithosphere history of partial melting and subsequent multi-stage melt–rock interaction. Cpx-poor spinel lherzolites are consistent with mantle residues after low-degree fractional melting (F = 5–10%). Opx + spinel symplectites at the rims of orthopyroxene porphyroclasts indicate post-melting lithospheric cooling (T = 970–1,100°C); this was followed by formation of olivine embayments within pyroxene porphyroclasts by melt–rock interaction. Enrichment in modal olivine (up to 85 wt%) at constant bulk Mg values, and variable absolute REE contents (at constant LREE/HREE) indicate olivine precipitation and pyroxene dissolution during reactive porous melt flow. This stage occurred at spinel-facies depths, after incorporation of the peridotites in the thermal lithosphere. Plagioclase-enriched peridotites show melt impregnation microtextures, like opx + plag intergrowths replacing exsolved cpx porphyroclasts and interstitial gabbronoritic veinlets. This second melt–rock interaction stage caused systematic chemical changes in clinopyroxene (e.g. Ti, REE, Zr, Y increase), related to the concomitant effects of local melt–rock interaction at decreasing melt mass, and crystallization of small (<3%) trapped melt fractions. LREE depletion in minerals of the gabbronoritic veinlets indicates that the impregnating melts were more depleted than normal MORB. Preserved microtextural evidence of previous melt–rock interaction in the impregnated peridotites suggests that they were progressively uplifted in response to lithosphere extension and thinning. Migrating melts were likely produced by mantle upwelling and melting related to extension; they were modified from olivine-saturated to opx-saturated compositions, and caused different styles of melt–rock interaction (reactive spinel harzburgites, vs. impregnated plagioclase peridotites) depending on the lithospheric depths at which interaction occurred. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Volcanic history, geologic analysis and map of the Prometheus Patera region on Io

Data from Jupiter's moon Io returned by the Galileo spacecraft have been used to create a geologic map of Prometheus Patera, its associated flow field, and nearby features. We have identified the location of the vent that fed the Prometheus flow field during the Galileo epoch in the north-eastern portion of the main Prometheus flow field. This vent is the probable source of a small sulphur-rich plume. Previous studies suggested that the vent may be atop a tectonic fault but we find that the vent is offset from the putative fault. It is plausible that, in the past, magma exploited the fault to reach the surface at Prometheus Patera, but subsequent magma cooling in the conduit could have caused an obstruction preventing further eruptions from providing significant contributions to the Prometheus flow field. We also speculate on how a new Prometheus plumbing system may be fed by mafic magmas after melt stalls in magma reservoirs during its ascent through the lithosphere from the mantle.     

Deformation of “Villafranchian” lacustrine sediments in the Chisone Valley (Western Alps, Italy)

The Chisone Valley is located in the internal NW Alps, in the Pinerolese District, an area characterized by present low to medium seismicity. Fine-grained sediments (sand, silt and clay with interbedded gravel) crop out in the lower Chisone Valley: they were first interpreted as glaciolacustrine deposits, and then as a lacustrine infilling of the valley floor probably due to differential uplifting of the valley mouth. Review of this data, together with new field and palynological observations, lead us to refer the lacustrine deposits to approximately the Lower Pleistocene (Villafranchian). In many outcrops, the lacustrine deposits show strong soft-sediment deformation such as convolute laminations, water-escape structures and disrupted beds, some of them associated with folds and faults (cm to dm in size); only two sites show metric to decametric folds and faults trending E-W and N-S. Detailed structural analysis conducted along a recently exposed section (Rio Gran Dubbione site) shows several soft-sediment deformation features on the limbs of mesoscale folds. Because of their intimate structural association, the origin of these minor structures seems to be connected to synsedimentary activity on reverse and normal faults (m to dm in size) affecting the lacustrine deposits in the same locality. Soft-sediment deformation features can be interpreted as possible paleoseismites. If so, the present seismicity of the Pinerolese District, which is the major area of such activity in NW Italy, cannot be considered an isolated episode in the geological evolution of the region; even if there is no supporting evidence for continuous seismicity, the deformations in the lacustrine sediments of the Chisone Valley testify to Early Pleistocene seismic activity, probably related to the recent tectonic evolution of the internal side of the NW Alps.     

Toward a full multiscale approach to interpret potential fields

The way potential fields convey source information depends on the scale at which the field is analysed. In this sense a multiscale analysis is a useful method to study potential fields particularly when the main field contributions are caused by sources with different depths and extents. Our multiscale approach is built with a stable transformation, such as depth from extreme points. Its stability results from mixing, in a single operator, the wavenumber low‐pass behaviour of the upward continuation transformation of the field with the enhancement high‐pass properties of n‐order derivative transformations. So, the complex reciprocal interference of several field components may be efficiently faced at several scales of the analysis and the depth to the sources may be estimated together with the homogeneity degrees of the field. In order to estimate the source boundaries we use another multiscale method, the multiscale derivative analysis, which utilizes a generalized concept of horizontal derivative and produces a set of boundary maps at different scales. We show through synthetic examples and application to the gravity field of Southern Italy that this multiscale behaviour makes this technique quite different from other source boundary estimators. The main result obtained by integrating multiscale derivative analysis with depth from extreme points is the retrieval of rather effective information of the field sources (horizontal boundaries, depth, structural index). This interpretative approach has been used along a specific transect for the analysis of the Bouguer anomaly field of Southern Apennines. It was set at such scales, so to emphasize either regional or local features along the transect. Two different classes of sources were individuated. The first one includes a broad, deep source with lateral size of 45∼50 km, at a depth of 13 km and having a 0.5 structural index. The second class includes several narrower sources located at shallowest depths, ranging from 3–6 km, with lateral size not larger than 5 km and structural indexes ranging from 1–1.5. Within a large‐scale geological framework, these results could help to outline the mean structural features at crustal depths.     

Simulating the role of tides and sediment characteristics on tidal flat sorting and bedding dynamics

Understanding sediment sorting and bedding dynamics has high value to unravelling the mechanisms underlying geomorphological, geological, ecological and environmental imprints of tidal wetlands and hence to predicting their future changes. Using the Nanhui tidal flat on the Changjiang (Yangtze) Delta, China, as a reference site, this study establishes a schematized morphodynamic model coupling flow, sediment dynamics and bed level change to explore the processes that govern sediment sorting and bedding phenomena. Model results indicate an overall agreement with field data in terms of tidal current velocities, suspended sediment concentrations (SSCs), deposition thicknesses and sedimentary structures. Depending on the variation of tidal current strength, sand-dominated layers (SDLs) and mud-dominated layers (MDLs) tend to form during spring and neap tides, respectively. Thinner tidal couplets are developed during daily scale flood–ebb variations. A larger tidal level variation during a spring–neap tidal cycle, associated with a stronger tidal current variation, favours the formation of SDLs and tidal couplets. A larger boundary sediment supply generally promotes the formation of tidal bedding, though the bedding detail is partially dependent on the SSC composition of different sediment types. Sediment properties, including for example grain size and settling velocity, are also found to influence sediment sorting and bedding characteristics. In particular, finer and coarser sediment respond differently to spring and neap tides. During neap tides, relatively small flow velocities favour the deposition of finer sediment, with limited coarser sediment being transported to the upper tidal flat because of the larger settling velocity. During spring tides, larger flow velocities transport more coarser sediment to the upper tidal flat, accounting for distinct lamination formation. Model results are qualitatively consistent with field observations, but the role of waves, biological processes and alongshore currents needs to be included in further studies to establish a more complete understanding.     

Spectroscopic Investigation on Hydrophobic and Hydrophilic Fractions of Dissolved Organic Matter Extracted from Soils at Different Salinities

Organic matter can be considered one of the most important indicators of the extent of soil desertification processes. Among the causes of desertification, salinization induced by different factors is raising the greatest concern in the Mediterranean area. In the present research, hydrophilic (HI) and hydrophobic (HO) fractions of dissolved organic matter (DOM) extracted from soils at different degrees of salinization have been investigated by means of spectroscopic techniques such as tridimensional fluorescence spectroscopy in the mode of emission excitation matrix (EEM) and Fourier transform infrared spectroscopy (FT‐IR). The FT‐IR spectra were distinctive in differentiating HI from HO fractions and each DOM fraction as a function of soil salinity. The EEM spectra of HO fractions exhibited a shift toward longer emission wavelengths and higher fluorescence intensity (FI) values as compared to that of the HI fractions. These results could be ascribed to the different molecular complexities of HI and HO fractions. Further, a marked quenching effect was observed in the FI of both the DOM fractions with increasing soil salinity, which allowed to obtain immediate information on the soil salinity degree by comparing the fluorescence intensity.     

Steady-state operation of Stromboli volcano,Italy: constraints on the feeding system

Stromboli volcano has been in continuous eruption for several thousand years without major changes in the geometry and feeding system. The thermal structure of its upper part is therefore expected to be close to steady state. In order to mantaim explosive activity, magma must release both gas and heat. It is shown that the thermal and gas budgets of the volcano lead to consistent conclusions. The thermal budget of the volcano is studied by means of a finite-element numerical model under the assumption of conduction heat transfer. It is found that the heat loss through the walls of an eruption conduit is weakly sensitive to the dimensions of underlying magma reservoirs and depends mostly on the radius and length of the conduit. In steady state, this heat loss must be balanced by the cooling of magma which flows through the system. For the magma flux of about 1 kg s^-1 corresponding to normal Strombolian activity, this requires that the conduits are a few meters wide and not deeper than a few hundred meters. This implies the existence of a magma chamber at shallow depth within the volcanic edifice. This conclusion is shown to be consistent with considerations on the thermal effects of degassing. In a Strombolian explosion, the mass ratio of gas to lava is very large, commonly exceeding two, which implies that the thermal evolution of the erupting mixture is dominated by that of the gas phase. The large energy loss due to decompression of the gas phase leads to decreased eruption temperatures. The fact that lava is molten upon eruption implies that the mixture does not rise from more than about 200 m depth. To sustain the magmatic and volcanic activity of Stromboli, a mass flux of magma of a few hundred kilograms per second must be supplied to the upper parts of the edifice. This represents either the rate of magma production from the mantle source feeding the volcano or the rate of magma overturn in the interior of a large chamber.