CO2 output and δ13C(CO2) from Mount Etna as indicators of degassing of shallow asthenosphere

 An estimated average CO₂ output from Etna's summit craters in the range of 13±3 Mt/a has recently been determined from the measured SO₂ output and measured CO₂/SO₂ molar ratios. To this amount the CO₂ output emitted diffusely from the soil (≈ 1 Mt/a) and the amount of CO₂ dissolved in Etna's aquifers (≈ 0.25 Mt/a) must be added. Data on the solubility of CO₂ in Etnean magmas at high temperature and pressure allow the volume of magma involved in the release of such an amount of this gas to be estimated. This volume of magma (≈ 0.7 km³/a) is approximately 20 times greater than the volume of magma erupted annually during the period 1971–1995. On the basis of C-isotopic data of CO₂ collected in the Etna area and of new hypotheses on the source of Mediterranean magmas, significant contributions of CO₂ from non-magmatic sources to the total output from Etna are unlikely. Such large outputs of CO₂ and also of SO₂ from Etna could be due to an anomalously shallow asthenosphere beneath the volcano that allows a continuous escape of gases toward the surface, even without migration of magma. Received: 7 August 1996 / Accepted: 9 November 1996     

The eruptive history of the Tequila volcanic field, western Mexico: ages, volumes, and relative proportions of lava types

The eruptive history of the Tequila volcanic field (1600 km²) in the western Trans-Mexican Volcanic Belt is based on ⁴⁰Ar/³⁹Ar chronology and volume estimates for eruptive units younger than 1 Ma. Ages are reported for 49 volcanic units, including Volcán Tequila (an andesitic stratovolcano) and peripheral domes, flows, and scoria cones. Volumes of volcanic units 1 Ma were obtained with the aid of field mapping, ortho aerial photographs, digital elevation models (DEMs), and ArcGIS software. Between 1120 and 200 kyrs ago, a bimodal distribution of rhyolite (~35 km³) and high-Ti basalt (~39 km³) dominated the volcanic field. Between 685 and 225 kyrs ago, less than 3 km³ of andesite and dacite erupted from more than 15 isolated vents; these lavas are crystal-poor and show little evidence of storage in an upper crustal chamber. Approximately 200 kyr ago, ~31 km³ of andesite erupted to form the stratocone of Volcán Tequila. The phenocryst assemblage of these lavas suggests storage within a chamber at ~2–3 km depth. After a hiatus of ~110 kyrs, ~15 km³ of andesite erupted along the W and SE flanks of Volcán Tequila at ~90 ka, most likely from a second, discrete magma chamber located at ~5–6 km depth. The youngest volcanic feature (~60 ka) is the small andesitic volcano Cerro Tomasillo (~2 km³). Over the last 1 Myr, a total of 128±22 km³ of lava erupted in the Tequila volcanic field, leading to an average eruption rate of ~0.13 km³/kyr. This volume erupted over ~1600 km², leading to an average lava accumulation rate of ~8 cm/kyr. The relative proportions of lava types are ~22–43% basalt, ~0.4–1% basaltic andesite, ~29–54% andesite, ~2–3% dacite, and ~18–40% rhyolite. On the basis of eruptive sequence, proportions of lava types, phenocryst assemblages, textures, and chemical composition, the lavas do not reflect the differentiation of a single (or only a few) parental liquids in a long-lived magma chamber. The rhyolites are geochemically diverse and were likely formed by episodic partial melting of upper crustal rocks in response to emplacement of basalts. There are no examples of mingled rhyolitic and basaltic magmas. Whatever mechanism is invoked to explain the generation of andesite at the Tequila volcanic field, it must be consistent with a dominantly bimodal distribution of high-Ti basalt and rhyolite for an 800 kyr interval beginning ~1 Ma, which abruptly switched to punctuated bursts of predominantly andesitic volcanism over the last 200 kyrs.Electronic Supplementary Material Supplementary material is available in the online version of this article at Editorial responsility: J. Donnelly-NolanThis revised version was published online in January 2005 with corrections to Tables 1 and 3.An erratum to this article can be found at     

Volcanic risk ranking for Auckland, New Zealand. I: Methodology and hazard investigation

Volcanic eruptions typically produce a number of hazards, and many regions are at risk from more than one volcano or volcanic field. So that detailed risk assessments can be carried out, it is necessary to rank potential volcanic hazards and events in terms of risk. As it is often difficult to make accurate predictions regarding the characteristics of future eruptions, a method for ranking hazards and events has been developed that does not rely on precise values. Risk is calculated individually for each hazard from each source as the product of likelihood, extent and effect, based on the parameters order of magnitude. So that multiple events and outcomes can be considered, risk is further multiplied by the relative probability of the event occurring (probability_e) and the relative importance of the outcome (importance_o). By adding the values obtained, total risk is calculated and a ranking can be carried out.This method was used to rank volcanic hazards and events that may impact the Auckland Region, New Zealand. Auckland is at risk from the Auckland volcanic field, Okataina volcanic centre, Taupo volcano, Tuhua volcano, Tongariro volcanic centre, and Mt. Taranaki volcano. Relative probabilities were determined for each event, with the highest given to Mt. Taranaki. Hazards considered were, for local events: tephra fall, scoria fall and ballistic impacts, lava flow, base surge and associated shock waves, tsunami, volcanic gases and acid rain, earthquakes and ground deformation, mudflows and mudfills, lightning and flooding; and for distal events: tephra fall, pyroclastic flows, poisonous gases and acid rain, mudflows and mudfills, climate variations and earthquakes. Hazards from each source were assigned values for likelihood, with the largest for tephra fall from all sources, earthquakes and ground deformation, lava flows, scoria fall and base surge for an Auckland eruption on land, and earthquakes and ground deformation from an Auckland eruption in the ocean. The largest values for extent were for tephra fall and climate variation from each of the distal centres. However, these parameters do not give a true indication of risk. In a companion paper the effect of each hazard is fully investigated and the risk ranking completed.     

Surface area, porosity and water adsorption properties of fine volcanic ash particles

Our understanding on how ash particles in volcanic plumes react with coexisting gases and aerosols is still rudimentary, despite the importance of these reactions in influencing the chemistry and dynamics of a plume. In this study, six samples of fine ash (<100 m) from different volcanoes were measured for their specific surface area, a_s, porosity and water adsorption properties with the aim to provide insights into the capacity of silicate ash particles to react with gases, including water vapour. To do so, we performed high-resolution nitrogen and water vapour adsorption/desorption experiments at 77 K and 303 K, respectively. The nitrogen data indicated a_s values in the range 1.1–2.1 m²/g, except in one case where a a_s of 10 m²/g was measured. This high value is attributed to incorporation of hydrothermal phases, such as clay minerals, in the ash surface composition. The data also revealed that the ash samples are essentially non-porous, or have a porosity dominated by macropores with widths >500 Å. All the specimens had similar pore size distributions, with a small peak centered around 50 Å. These findings suggest that fine ash particles have relatively undifferentiated surface textures, irrespective of the chemical composition and eruption type. Adsorption isotherms for water vapour revealed that the capacity of the ash samples for water adsorption is systematically larger than predicted from the nitrogen adsorption a_s values. Enhanced reactivity of the ash surface towards water may result from (i) hydration of bulk ash constituents; (ii) hydration of surface compounds; and/or (iii) hydroxylation of the surface of the ash. The later mechanism may lead to irreversible retention of water. Based on these experiments, we predict that volcanic ash is covered by a complete monolayer of water under ambient atmospheric conditions. In addition, capillary condensation within ash pores should allow for deposition of condensed water on to ash particles before water reaches saturation in the plume. The total mass of water vapour retained by 1 g of fine ash at 0.95 relative water vapour pressure is calculated to be ~10^–2 g. Some volcanic implications of this study are discussed.Editorial responsibility: J. Gilbert     

The Tambora eruption in 1815 provides a test on possible global climatic and chemical perturbations in the past

A coupled one-dimensional radiative-convective-photochemical diffusion model, which takes into account the influence of ocean inertia on global radiative perturbations is used to investigate the possible climatic and other atmospheric effects of a major volcanic eruption, thought to be similar in magnitude to that of the Tambora eruption, Indonesia, which took place in 1815. A volcanic cloud was introduced in the model stratosphere between 20–25 km and the global average peak aerosol optical thickness was assumed to be 0.25. Both the aerosol optical thickness and aerosol composition, which determine the optical properties, were allowed to vary in the model atmosphere during the life cycle of the volcanic cloud. The results indicate that the global average surface temperature decreases steadily from the date of eruption (7–12 April 1815) with maximum cooling of 1° K occurring in the spring of 1816. The calculations also show significant warming of the stratosphere, with temperature increasing up to 15° K at 25 km in less than six months after the date of eruption. The important effects of the Tambora eruption on stratospheric ozone and UV-B radiation at the surface are also mentioned.     

Evolution of a complex isolated dome system,Cerro Pizarro,central México

Cerro Pizarro is an isolated rhyolitic dome in the intermontane Serdán-Oriental basin, located in the eastern Trans-Mexican Volcanic Belt. Cerro Pizarro erupted ~1.1 km³ of magma at about 220 ka. Activity of Cerro Pizarro started with vent-clearing explosions at some depth; the resultant deposits contain clasts of local basement rocks, including Cretaceous limestone, ~0.46-Ma welded tuff, and basaltic lava. Subsequent explosive eruptions during earliest dome growth produced an alternating sequence of surge and fallout layers from an inferred small dome. As the dome grew both vertically and laterally, it developed an external glassy carapace due to rapid chilling. Instability of the dome during emplacement caused the partial gravitational collapse of its flanks producing various block-and-ash-flow deposits. After a brief period of repose, re-injection of magma caused formation of a cryptodome with pronounced deformation of the vitrophyric dome and the underlying units to orientations as steep as near vertical. This stage began apparently as a gas-poor eruption and no explosive phases accompanied the emplacement of the cryptodome. Soon after emplacement of the cryptodome, however, the western flank of the edifice catastrophically collapsed, causing a debris avalanche. A hiatus in eruptive activity was marked by erosion of the cone and emplacement of ignimbrite derived from a caldera to the north of Cerro Pizarro. The final growth of the dome growth produced its present shape; this growth was accompanied by multiple eruptions producing surge and fallout deposits that mantle the topography around Cerro Pizarro. The evolution of the Cerro Pizarro dome holds aspects in common with classic dome models and with larger stratovolcano systems. We suggest that models that predict a simple evolution for domes fail to account for possibilities in evolutionary paths. Specifically, the formation of a cryptodome in the early stages of dome formation may be far more common than generally recognized. Likewise, sector collapse of a dome, although apparently rare, is a potential hazard that must be recognized and for which planning must be done.Editorial responsibility: J. Gilbert     

Geometry and growth of sill complexes: insights using 3D seismic from the North Rockall Trough

Doleritic sill complexes, which are an important component of volcanic continental margins, can be imaged using 3D seismic reflection data. This allows unprecedented access to the complete 3D geometry of the bodies and an opportunity to test classic sill emplacement models. The doleritic sills associated with basaltic volcanism in the North Rockall Trough occur in two forms. Radially symmetrical sill complexes consist of a saucer-like inner sill at the base with an arcuate inclined sheet connecting it to a gently inclined, commonly ragged, outer rim. Bilaterally symmetrical sill complexes are sourced by magma diverted from a magma conduit feeding an overlying volcano. With an elongate, concave upwards, trough-like geometry bilaterally symmetrical sills climb away from the magma source from which they originate. Both sill complex types can appear as isolated bodies but commonly occur in close proximity and consequently merge, producing hybrid sill complexes. Radial sill complexes consist of a series of radiating primary flow units. With dimensions up to 3 km, each primary flow unit rises from the inner saucer and is fed by primary magma tube. Primary flow units contain secondary flow units with dimensions up to 2 km, each being fed by a secondary magma tube branching from the primary magma tube. Secondary flow units in turn are composed of 100-m scale tertiary flow units. A similar branching hierarchy of flow units can also be seen in bilaterally symmetrical sill complexes, with their internal architecture resembling an enlarged version of a primary flow unit from a radial sill complex. This branching flow pattern, as well as the interaction between flow units of varying orders, provides new insights into the origin of the structures commonly seen within sill complexes and the hybrid sill bodies produced by their merger. The data demonstrate that each radially symmetrical sill complex is independently fed from a source located beneath the centre of the inner saucer, grows by climbing from the centre outwards and that peripheral dyking from the upper surface is a common feature. These features suggest a laccolith emplacement style involving peripheral fracturing and dyking during inner saucer growth and thickening. The branching hierarchy of flow units within bilaterally symmetrical sill complexes is broadly similar to that of primary flow units within a radially symmetrical sill complex, suggesting that the general features of the laccolith emplacement model also apply.Editorial responsibility: J. Stix     

Doppler radar sounding of volcanic eruption dynamics at Mount Etna

Besides their common use in atmospheric studies, Doppler radars are promising tools for the active remote sensing of volcanic eruptions but were little applied to this field. We present the observations made with a mid-power UHF Doppler radar (Voldorad) during a 7-h Strombolian eruption at the SE crater of Mount Etna on 11–12 October 1998. Main characteristics of radar echoes are retrieved from analysis of Doppler spectra recorded in the two range gates on either side of the jet axis. From the geometry of the sounding, the contribution of uprising and falling ejecta to each Doppler spectrum can be discriminated. The temporal evolution of total power backscattered by uprising targets is quite similar to the temporal evolution of the volcanic tremor and closely reproduces the overall evolution of the eruption before, during and after its paroxysm. Moreover, during the sharp decrease of eruptive activity following the paroxysm, detailed analysis of video (from camera recording), radar and seismic measurements reveals that radar and video signals start to decrease simultaneously, approximately 2.5 min after the tremor decline. This delay is interpreted as the ascent time through a magma conduit of large gas slugs from a shallow source roughly estimated at about 500 m beneath the SE crater. Detailed analysis of eruptive processes has been also made with Voldorad operating in a high sampling rate mode. Signature of individual outburst is clearly identified on the half part of Doppler spectra corresponding to rising ejecta: temporal variations of the backscattered power exhibit quasi periodic undulations, whereas the maximum velocity measured on each spectrum displays a sharp peak at the onset of each outburst followed by a slow decay with time. Periodicity of power variations (between 3.8 and 5.5 s) is in agreement with the occurrence of explosions visually observed at the SE vent. Maximum vertical velocities of over 160 m s^–1 were measured during the paraoxysmal stage and the renewed activity. Finally, by using a simplified model simulating the radar echoes characteristics, we show that when Voldorad is operating in high sampling rate mode, the power and maximum velocity variations are directly related to the difference in size and velocity of particles crossing the antenna beam.Editorial responsibility: A. Woods     

Nonlinear Time Series Analysis of Volcanic Tremor Events Recorded at Sangay Volcano,Ecuador

The assumption that volcanic tremor may be generated by deterministic nonlinear source processes is now supported by a number of studies at different volcanoes worldwide that clearly demonstrate the low-dimensional nature of the phenomenon. We applied methods based on the theory of nonlinear dynamics to volcanic tremor events recorded at Sangay volcano, Ecuador in order to obtain more information regarding the physics of their source mechanism. The data were acquired during 21–26 April 1998 and were recorded using a sampling interval of 125 samples s^–1 by two broadband seismometers installed near the active vent of the volcano. In a previous study Johnson and Lees (2000) classified the signals into three groups: (1) short duration (<1 min) impulses generated by degassing explosions at the vent; (2) extended degassing chugging events with a duration 2–5 min containing well-defined integer overtones (1–5 Hz) and variable higher frequency content; (3) extended degassing events that contain significant energy above 5 Hz. We selected 12 events from groups 2 and 3 for our analysis that had a duration of at least 90 s and high signal-to-noise ratios. The phase space, which describes the evolution of the behavior of a nonlinear system, was reconstructed using the delay embedding theorem suggested by Takens. The delay time used for the reconstruction was chosen after examining the first zero crossing of the autocorrelation function and the first minimum of the Average Mutual Information (AMI) of the data. In most cases it was found that both methods yielded a delay time of 14–18 samples (0.112–0.144 s) for group 2 and 5 samples (0.04 s) for group 3 events. The sufficient embedding dimension was estimated using the false nearest neighbors method which had a value of 4 for events in group 2 and was in the range 5–7 for events in group 3. Based on these embedding parameters it was possible to calculate the correlation dimension of the resulting attractor, as well as the average divergence rate of nearby orbits given by the largest Lyapunov exponent. Events in group 2 exhibited lower values of both the correlation dimension (1.8–2.6) and largest Lyapunov exponent (0.013–0.022) in comparison with the events in group 3 where the values of these quantities were in the range 2.4–3.5 and 0.029–0.043, respectively. Theoretically, a nonlinear oscillation described by the equation ++g(x)=fcost can generate deterministic signals with characteristics similar to those observed in groups 2 and 3 as the values of the parameters ,,f, are drifting, causing instability of orbits in the phase space.     

松辽盆地徐家围子断陷火山岩气藏天然气成因类型和气源对比

松辽盆地徐家围子断陷发育营城组火山岩气藏.钻井统计、地球化学特征及天然气组份分析和碳同位素研究表明徐中和徐东地区天然气源岩有机质偏腐泥型,生气潜力更大.对主力源岩层沙河子组和次要源岩层营城组源岩分布有了新的认识.对营城组天然气的烃类气体和CO₂气的来源进行对比和探讨,认为烃类气受气源的控制,CO₂气藏受深大断裂的控制,与火山岩分布关系不大.连通地幔的气藏中,也没有发现无机烃类气体聚集.徐家围子断陷营城组天然气以烃类气占绝对优势,而且主要为甲烷,二氧化碳成藏并不普遍.天然气甲烷成因类型以煤型气为主,甲烷碳同位素变化范围较小.乙烷、丙烷和丁烷碳同位素均有较大的变化范围,而且具有较好的相关性,说明母质类型多样,乙烷、丙烷和丁烷具有相同的来源.含量大于50%的二氧化碳来源于地幔,属无机成因.含量小于50%的二氧化碳既有无机成因,也有有机成因.