JOURNAL OF VULCANOLOGY AND GEOTHERMAL RESEARCH

Resumen biblográfico de revista Vulcanology and Geothermal

Using the level set method to model endogenous lava dome growth

Modeling volcanic phenomena is complicated by free-surfaces often supporting large rheological gradients. Analytical solutions and analogue models provide explanations for fundamental characteristics of lava flows. But more sophisticated models are needed, incorporating improved physics and rheology to capture realistic events. To advance our understanding of the flow dynamics of highly viscous lava in Peléean lava dome formation, axi-symmetrical Finite Element Method (FEM) models of generic endogenous dome growth have been developed. We use a novel technique, the level-set method, which tracks a moving interface, leaving the mesh unaltered. The model equations are formulated in an Eulerian framework. In this paper we test the quality of this technique in our numerical scheme by considering existing analytical and experimental models of lava dome growth which assume a constant Newtonian viscosity. We then compare our model against analytical solutions for real lava domes extruded on Soufrière, St. Vincent, W.I. in 1979 and Mount St. Helens, USA in October 1980 using an effective viscosity. The level-set method is found to be computationally light and robust enough to model the free-surface of a growing lava dome. Also, by modeling the extruded lava with a constant pressure head this naturally results in a drop in extrusion rate with increasing dome height, which can explain lava dome growth observables more appropriately than when using a fixed extrusion rate. From the modeling point of view, the level-set method will ultimately provide an opportunity to capture more of the physics while benefiting from the numerical robustness of regular grids.

Magma Flow Instabilities in a Volcanic Conduit: Implications for Long-Period Seismicity

Silicic volcanic eruptions are typically accompanied by repetitive Long-Period (LP) seismicity that originates from a small region of the upper conduit. These signals have the capability to advance eruption prediction, since they commonly precede a change in the eruption vigour. Shear bands forming along the conduit wall, where the shear stresses are highest, have been linked to providing the seismic trigger. However, existing computational models are unable to generate shear bands at the depths where the LP signals originate using simple magma strength models. Presented here is a model in which the magma strength is determined from a constitutive relationship dependent upon crystallinity and pressure. This results in a depth-dependent magma strength, analogous to planetary lithospheres. Hence, in shallow highly-crystalline regions a macroscopically discontinuous brittle type of deformation will prevail, whilst in deeper crystal-poor regions there will be a macroscopically continuous plastic deformation mechanism. This will result in a depth where the brittle-ductile transition occurs, and here shear bands disconnected from the free-surface may develop. We utilize the Finite Element Method and use axi-symmetric coordinates to model magma flow as a viscoplastic material, simulating quasi-static shear bands along the walls of a volcanic conduit. Model results constrained to the Soufrière Hills Volcano, Montserrat, show the generation of two types of shear bands: upper-conduit shear bands that form between the free-surface to a few 100 metres below it and discrete shear bands that form at the depths where LP seismicity is measured to occur corresponding to the brittle-ductile transition and the plastic shear region. It is beyond the limitation of the model to simulate a seismic event, although the modelled viscosity within the discrete shear bands suggests a failure and healing cycle time that supports the observed LP seismicity repeat times. However, due to the paucity of data and large parameter space available these results can only be considered to be qualitative rather than quantitative at this stage.

Yarrol terrane of the northern New England Fold Belt: forearc or backarc?

The Upper Devonian to Lower Carboniferous volcanosedimentary rocks of the Yarrol terrane of the northern New England Fold Belt have previously been ascribed to a forearc basin setting. New data presented here, however, suggest that the Yarrol terrane developed as a backarc basin during the Middle to early Late Devonian. Based on field studies, we recognise four regionally applicable strati graphic units: (i) a basal, ?Middle to Upper Devonian submarine mafic volcanic suite (Monal volcanic facies association); (ii) the lower Frasnian Lochenbar beds that locally unconformably overlie the Monal volcanic facies association: (iii) the Three Moon Conglomerate (Upper Devonian - Lower Carboniferous): and (iv) the Lower Carboniferous Rockhampton Group characterised by the presence of oolitic limestone. Stratigraphic and compositional differences suggest the Monal volcanic facies association post-dates Middle Devonian silicic-dominated magmatism that was coeval with gold-copper mineralisation at Mt Morgan. The Lochenbar beds, Three Moon Conglomerate and Rockhampton Group represent a near-continuous sedimentary record of volcanism that changed in composition and style from mafic effusive (Late Devonian) to silicic explosive volcanism (Early Carboniferous). Palaeocurrent data from the Three Moon Conglomerate and Rockhampton Group indicate dispersal of sediment to the west and northwest, and are inconsistent with derivation from a volcanic-are source situated to the west (Connors-Auburn Arch). Geochemical data show that the Monal volcanic facies association ranges from tholeiitic subalkaline basalts to calc-alkaline basaltic andesite. Trace and rare-earth element abundances are distinctly MORE-like (e.g, light rare earth element depletion), with only moderate enrichment of the large-ion lithophile elements in some units, and negative Nb anomalies, suggesting a subduction-related signature. Basalts of the Monal volcanic facies association are best described as transitional between calc-alkali basalts and N-MORB. The elevated high field strength element contents (e.g. Zr, Y, Ti) are higher than modern island-are basalts, but comparable to basalts that floor modern backarc basins. This geochemical study, coupled with stratigraphic relationships, suggest that the eruption of backarc basin basalts followed widespread Middle Devonian, extension-related silicic magmatism (e.g. Retreat Batholith, Mt Morgan), and floored the Yarrol terrane. The Monal volcanic facies association thus shows similarities in its tectonic environment to the Lower Permian successions (e.g. Rookwood Volcanics) of the northern New England Fold Belt. These mafic volcanic sequences are interpreted to record two backarc basin-forming periods (Middle - Late Devonian and Late Carboniferous - Early Permian) during the Late Palaeozoic history of the New England Orogen. Silicic-dominated explosive volcanism, occurring extensively across the northern New England Fold Belt in the Early Carboniferous (Varrol terrane, Campwyn Volcanics, Drummond and Burdekin Basins), reflects another period of crustal melting and extension, most likely related to the opening of the Drummond Basin.

Welding and rheomorphism of phonolitic fallout deposits from the Las Canadas caldera, Tenerife, Canary Islands

The Las Canadas caldera is a nested collapse caldera formed by the successive migration and collapse of shallow magmatic chambers. Among the pyroclastic products of this caldera are phonolitic fallout deposits that crop out in the caldera wall and on the extracaldera slopes. These deposits exhibit an uninterrupted facies gradation from nonwelded to lava-like and record continuous volcanic deposition. Densely welded and lava-like facies result from the extreme attenuation and complete homogenization of juvenile clasts that destroy original clast outlines and any evidence of fallout deposition. Agglutination contributes significantly to the final degree of flattening observed in the welded facies. After deposition, rheomorphic flowage occurs. Emplacement temperatures for one of the welding sequences are calculated from magmatic temperatures and a model of tephra cooling during fallout. Results are 486 degreesC for the nonwelded facies and 740 degreesC for the moderately welded facies. For the same welding sequence, a cooling time between 25 and 54 days is estimated from published experimental and computational data as the possible duration of welding and rheomorphism. Following deposition and agglutination, the lava-like pyroclastic facies had the rheological properties of viscous lavas and flowed down the outer slopes away from the caldera. Some lava-like masses detached from proximal areas to more distal regions. During deposition, the eruptive style evolved from Plinian fallout to fountain-fed spatter deposition. This evolution was accompanied by a decrease in explosive power and a lower height of the eruptive column, which produce higher emplacement temperatures and more effective heat retention of pyroclasts.

Heat pulse line-source method to determine thermal conductivity of consolidated rocks

An instrument designed to measure thermal conductivity of consolidated rocks, dry or saturated, using a transient method is presented. The instrument measures relative values of the thermal conductivity, and it needs calibration to obtain absolute values. The device can be used as heat pulse line source and as continuous heat line source. Two parameters to determine thermal conductivity are proposed: TMAX, in heat pulse line source, and SLOPE, in continuous heat line source. Its performance is better, and the operation simpler, in heat pulse line-source mode with a measuring time of 170 s and a reproducibility better than 2.5%. The sample preparation is very simple on both modes. The performance has been tested with a set of ten rocks with thermal conductivity values between 1.4 and 5.2 W m¿1 K¿1 which covers the usual range for consolidated rocks.

The Quaternary pyroclastic succession of southeast Tenerife, Canary Islands: Explosive eruptions, related caldera subsidence, and sector collapse

A much-revised Quaternary stratigraphy is presented for ignimbrites and pumice fall deposits of the Bandas del Sur, in southern Tenerife. New Ar-41/Ar-39 data obtained for the Arico, Granadilla, Fasnia, Poris, La Caleta and Abrigo formations are presented, allowing correlation with previously dated offshore marine ashfall layers and volcaniclastic sediments. We also provide a minimum age of 287 +/- 7 ka for a major sector collapse event at the Gaimar valley. The Bandas del Sur succession includes more than seven widespread ignimbrite sheets that have similar characteristics, including widespread basal Plinian layers, predominantly phonolite composition, ignimbrites with similar extensive geographic distributions, thin condensed veneers with abundant diffuse bedding and complex lateral and vertical grading patterns, lateral gradations into localized massive facies within palaeo-wadis, and widespread lithic breccia layers that probably record caldera-forming eruptions. Each ignimbrite sheet records substantial bypassing of pyroclastic material into the ocean. The succession indicates that Las Canadas volcano underwent a series of major explosive eruptions, each starting with a Plinian phase followed by emplacement of ignimbrites and thin ash layers, some of coignimbrite origin. Several of the ignimbrite sheets are compositionally zoned and contain subordinate mafic pumices and banded pumices indicative of magma mingling immediately prior to eruption. Because passage of each pyroclastic density current was characterized by phases of non-deposition and erosion, the entire course of each eruption is incompletely recorded at any one location, accounting for some previously perceived differences between the units. Because each current passed into the ocean, estimating eruption volumes is virtually impossible. Nevertheless, the consistent widespread distributions and the presence of lithic breccias within most of the ignimbrite sheets suggest that at least seven caldera collapse eruptions are recorded in the Bandas del Sur succession and probably formed a complex, nested collapse structure. Detailed field relationships show that extensive ignimbrite sheets (e.g. the Arico, Poris and La Caleta formations) relate to previously unrecognized caldera collapse events. We envisage that the evolution of the nested Las Cahadas caldera is more complex than previously thought and involved a protracted history of successive ignimbrite-related caldera collapse events, and large sector collapse events, interspersed with edifice-building phases.

Volcanic setting of the Bajo de la Alumbrera porphyry Cu-Au deposit, Farallon Negro Volcanics, Northwest Argentina

The late Miocene Farallon Negro volcanics, comprising basaltic to rhyodacitic volcano-sedimentary rocks, host the Bajo de la Alumbrera porphyry copper-gold deposit in northwest Argentina. Early studies of the geology of the district have underpinned the general model for porphyry ore deposits where hydrothermal alteration and mineralization develop in and around porphyritic intrusions emplaced at shallow depths (2.5-3.5 km) into stratovolcanic assemblages. The Farallon Negro succession is dominated by thick sequences of volcano-sedimentary breccias, with lavas forming a minor component volumetrically. These volcaniclastic rocks conformably overlie crystalline basement-derived sedimentary rocks deposited in a developing foreland basin southeast of the Puna-Altiplano plateau. Within the Farallon Negro volcanics, volcanogenic accumulations evolved from early mafic to intermediate and silicic compositions. The younger and more silicic rocks are demonstrably coeval and comagmatic with the earliest group of mineralized porphyritic intrusions at Bajo de la Alumbrera. Our analysis of the volcanic stratigraphy and facies architecture of the Farallon Negro volcanics indicates that volcanic eruptions evolved from effusive to mixed effusive and explosive styles, as magma compositions changed to more intermediate and silicic compositions. Air early phase of mafic to intermediate voleanism was characterized by small synsedimentary intrusions with peperitic contacts, and lesser lava units scattered widely throughout the district, and interbedded with thick and extensive successions of coarse-grained sedimentary breccias. These sedimentary breccias formed from numerous debris- and hyperconcentrated flow events. A later phase of silicic volcanism included both effusive eruptions, forming several areally restricted lavas, and explosive eruptions, producing more widely dispersed (up to 5 kin) tuff units, some tip to 30-m thickness in proximal sections. Four key features of the volcanic stratigraphy suggest that the Farallon Negro volcanics need not simply record the construction of a large steep-sided polygenetic stratovolcano: (1) sheetlike, laterally continuous debris-flow and other coarse-grained sedimentary deposits are dominant, particularly in the lower sections; (2) mafic-intermediate composition lavas are volumetrically minor; (3) peperites are present throughout the sequence; and (4) fine-grained lacustrine sandstone-siltstone sequences occur in areas previously thought to be proximal to the summit region of the stratovolcano. Instead, the nature, distribution, and geometry of volcanic and volcaniclastic facies suggest that volcanism occurred as a relatively low relief, multiple-vent volcanic complex at the eastern edge of a broad, > 200-km-wide late Miocene volcanic belt and oil ail active foreland sedimentary basin to the Puna-Altiplano. Volcanism that occurred synchronously with the earliest stages of porphyry-related mineralization at Bajo de la Alumbrera apparently developed in an alluvial to ring plain setting that was distal to larger volcanic edifices.