Physical volcanology, sedimentology, stratigraphy and petrochemistry of the Berry Creek metavolcanics: an Archean calc-alkaline complex, Lake of the Woods, Ontario

The steeply dipping, isoclinally folded early Precambrian (Archean) Berry Creek Metavolcanic Complex comprises primary to resedimented pyroclastic, epiclastic and autoclastic deposits. Tephra erupted from central volcanic edifices was dumped by mass flow mechanisms into peripheral volcanosedimentary depressions. Sedimentation has been essentially contemporaneous with eruption and transport of tephra. The monolithic to heterolithic tuffaceous horizons are interpreted as subaerial to subaqueous pumice and ash flows, secondary debris flows, lahars, slump deposits and turbidites. Monolithic debris flows, derived from crumble breccia and dcme talus, formed during downslope collapse and subsequent gravity flowage. Heterolithic tuff, lahars and lava flow morphologies suggest at least temporary emergence of the edifice. Local collapse may have accompanied pyroclastic volcanism. The tephra, produced by hydromagmatic to magmatic eruptions, were rapidly transported, by primary and secondary mechanisms, to a shallow littoral to deep water subaqueous fan developed upon the subjacent mafic metavolcanic platform. Deposition resulted from traction, traction carpet, and suspension sedimentation from laminar to turbulent flows. Facies mapping revealed proximal (channel to overbank) to distal facies epiclastics (greywackes, argillite) intercalated with proximal vent to medial fan facies crystal rich ash flows, debris flows, bedded tuff and shallow water to deep water lava flows. Framework and matrix support debris flows exhibit a variety of subaqueous sedimentary structures, e.g., coarse tail grading, double grading, inverse to normal grading, graded stratified pebbly horizons, erosional channels. Pelitic to psammitic AE turbidites also contain primary stru~tures, e.g., flames, load casts, dewatering pipes. Despite low to intermediate pressure greenschist to amphibolite grade metamorphism and variably penetrative deformation, relicts of pumice fragments and shards were recognized as recrystallized quartzofeldspathic pseudomorphs. The mafic to felsic metavolcanics and metasediments contain blasts of hornblende, actinolite, garnet, pistacitic epidote, staurolite, albitic plagioclase, and rarely andalusite and cordierite. The mafic metavolcanics (Adams River Bay, Black River, Kenu Lake, Lobstick Bay, Snake Bay) display _holeiitic trends with komatiitic affinities. Chemical variations are consistent with high level fractionation of olivine, plagioclase, amphibole, and later magnetite from a parental komatiite. The intermediate to felsic (64-74% Si02) metavolcanics generally exhibit calc-alkaline trends. The compositional discontinuity, defined by major and trace element diversity, can be explained by a mechanism involving two different magma sources. Application of fractionation series models are inconsistent with the observed data. The tholeiitic basalts and basaltic andesites are probably derived by low pressure fractionation of a depleted (high degree of partial melting) mantle source. The depleted (low Y, Zr) calc-alkaline metavolcanics may be produced by partial melting of a geochemically evolved source, e.g., tonalitetrondhjemite, garnet amphibolite or hydrous basalt.

TECHNOLOGICAL ADVANCEMENT IN REMOTE SENSING AND INNOVATIVE APPLICATION TO VOLCANOLOGY

The integration of remote monitoring techniques at different scales is of crucial importance for monitoring of volcanoes and assessment of the associated hazard. In this optic, technological advancement and collaboration between research groups also play a key role. Vhub is a community cyberinfrastructure platform designed for collaboration in volcanology research. Within the Vhub framework, this dissertation focuses on two research themes, both representing novel applications of remotely sensed data in volcanology: advancement in the acquisition of topographic data via active techniques and application of passive multi-spectral satellite data to monitoring of vegetated volcanoes. Measuring surface deformation is a critical issue in analogue modelling of Earth science phenomena. I present a novel application of the Microsoft Kinect sensor to measurement of vertical and horizontal displacements in analogue models. Specifically, I quantified vertical displacement in a scaled analogue model of Nisyros volcano, Greece, simulating magmatic deflation and inflation and related surface deformation, and included the horizontal component to reconstruct 3D models of pit crater formation. The detection of active faults around volcanoes is of importance for seismic and volcanic hazard assessment, but not a simple task to be achieved using analogue models. I present new evidence of neotectonic deformation along a north-south trending fault from the Mt Shasta debris avalanche deposit (DAD), northern California. The fault was identified on an airborne LiDAR campaign of part of the region interested by the DAD and then confirmed in the field. High resolution LiDAR can be utilized also for geomorphological assessment of DADs, and I describe a size-distance analysis to document geomorphological aspects of hummock in the Shasta DAD. Relating the remote observations of volcanic passive degassing to conditions and impacts on the ground provides an increased understanding of volcanic degassing and how satellite-based monitoring can be used to inform hazard management strategies in nearreal time. Combining a variety of satellite-based spectral time series I aim to perform the first space-based assessment of the impacts of sulfur dioxide emissions from Turrialba volcano, Costa Rica, on vegetation in the surrounding environment, and establish whether vegetation indices could be used more broadly to detect volcanic unrest.

Towards realistic simulations of lava dome growth using the level set method

The level set method has been implemented in a computational volcanology context. New techniques are presented to solve the advection equation and the reinitialisation equation. These techniques are based upon an algorithm developed in the finite difference context, but are modified to take advantage of the robustness of the finite element method. The resulting algorithm is tested on a well documented Rayleigh–Taylor instability benchmark [19], and on an axisymmetric problem where the analytical solution is known. Finally, the algorithm is applied to a basic study of lava dome growth.

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.

Aplicação de técnicas de geodesia espacial ao estudo dos sistemas vulcano-tectónicos e hidrotermais do segmento definido pelas ilhas Terceira, São Jorge e Graciosa

Tese de Doutoramento em Geologia, especialidade de Geodesia.

Reconstructing the architectural evolution of volcanic islands from combined K/AR, morphologic, tectonic, and magnetic data. The Faial Island exemple (Azores)

The morpho-structural evolution of oceanic islands results from competition between volcano growth and partial destruction by mass-wasting processes. We present here a multi-disciplinary study of the successive stages of development of Faial (Azores) during the last 1 Myr. Using high-resolution digital elevation model (DEM), and new K/Ar, tectonic, and magnetic data, we reconstruct the rapidly evolving topography at successive stages, in response to complex interactions between volcanic construction and mass wasting, including the development of a graben. We show that: (1) sub-aerial evolution of the island first involved the rapid growth of a large elongated volcano at ca. 0.85 Ma, followed by its partial destruction over half a million years; (2) beginning about 360 ka a new small edifice grew on the NE of the island, and was subsequently cut by normal faults responsible for initiation of the graben; (3) after an apparent pause of ca. 250 kyr, the large Central Volcano (CV) developed on the western side of the island at ca 120 ka, accumulating a thick pile of lava flows in less than 20 kyr, which were partly channelized within the graben; (4) the period between 120 ka and 40 ka is marked by widespread deformation at the island scale, including westward propagation of faulting and associated erosion of the graben walls, which produced sedimentary deposits; subsequent growth of the CV at 40 ka was then constrained within the graben, with lava flowing onto the sediments up to the eastern shore; (5) the island evolution during the Holocene involves basaltic volcanic activity along the main southern faults and pyroclastic eruptions associated with the formation of a caldera volcano-tectonic depression. We conclude that the whole evolution of Faial Island has been characterized by successive short volcanic pulses probably controlled by brief episodes of regional deformation. Each pulse has been separated by considerable periods of volcanic inactivity during which the Faial graben gradually developed. We propose that the volume loss associated with sudden magma extraction from a shallow reservoir in different episodes triggered incremental downward graben movement, as observed historically, when immediate vertical collapse of up to 2 m was observed along the western segments of the graben at the end of the Capelinhos eruptive crises (1957-58).

Upper crustal structure of Madeira Island revealed from ambient noise tomography

We present the first image of the Madeira upper crustal structure, using ambient seismic noise tomography. 16 months of ambient noise, recorded in a dense network of 26 seismometers deployed across Madeira, allowed reconstructing Rayleigh wave Green's functions between receivers. Dispersion analysis was performed in the short period band from 1.0 to 4.0 s. Group velocity measurements were regionalized to obtain 20 tomographic images, with a lateral resolution of 2.0 km in central Madeira. Afterwards, the dispersion curves, extracted from each cell of the 2D group velocity maps, were inverted as a function of depth to obtain a 3D shear wave velocity model of the upper crust, from the surface to a depth of 2.0 km. The obtained 3D velocity model reveals features throughout the island that correlates well with surface geology and island evolution. (C) 2015 Elsevier B.V. All rights reserved.

The role of basalt replenishment in the generation of basaltic andesites of the ongoing activity at Arenal volcano, Costa Rica: evidence from clinopyroxene and spinel

The bulk composition of magma erupted from Volcan Arenal has remained nearly constant (SiO2 = 53.6-54.9 wt%; MgO = 5.0-4.5 wt%) during almost 30 years of continuous activity (1969-1996). None the less, clinopyroxene (cpx) phenocrysts and their spinel inclusions record a much more complex open-system evolution in which steady-state production of the erupted basaltic andesitic magma is linked to episodic injections of basalt into Arenal's magma conduit/reservoir system. High-resolution major element zoning profiles (electron microprobe) on a large number of phenocrysts (>14,000 analyses), tied to back-scattered electron (BSE) images, have been used to assess the compositional characteristics of the magmatic end members as well as the timing and dynamics of magma replenishment events. No two cpx phenocrysts have exactly the same zoning profile. The vast majority of our analyses record the crystallization of cpx (Cr2O3 < 0.12 wt%; Mg# = 65-79; Al/Ti = 2-7) from a liquid comparable to or more evolved than erupted magma compositions. However, half of all cpx grains are cored by high-Cr cpx (Cr2O3 = 0.2-0.72 wt%) or contain similar basaltic compositions as abrupt growth bands in phenocrysts with and without high-Cr cores; phenocrysts with high-Cr cpx occur throughout the ongoing activity. In a few cases, high-Cr cpx occurs very near the outer margin of the grain without an apparent growth hiatus, particularly in 1968/69 and 1992/93. The main conclusions are: (1) all basaltic andesitic lavas erupted at Arenal during the ongoing activity that began in July, 1968, are the products of magma mixing, (2) clinopyroxenes record multiple replenishment events of basaltic magma in contrast to the near constancy of erupted bulk compositions, (3) some phenocrysts preserve records of multiple interactions with basaltic magmas requiring magmatic processes to operate on time-scales shorter than residence times of some phenocrysts, (4) multiple occurrences of clinopyroxene with high-Cr rims suggest that basalt replenishment events have occurred with sub-decadal frequency and may predate eruption by months or less. From this we infer that Arenal volcano is underlain by a continuously active, small-volume magmatic reservoir maintained in quasi-steady state by basalt recharge over several decades. The monotony of erupting Arenal magmas implies that fractionation, recharge, ascent, and eruption are well balanced in order for magmas to be essentially uniform while containing phenocrysts with vastly different growth histories at the time of eruption.

Geochemical evidence for mixing of magmatic fluids with seawater, Nisyros hydrothermal system, Greece

The chemical and isotopic compositions (deltaD(H2O), delta(18)O(H2O), delta(18)O(CO2), delta(13)C(CO2), delta(34)S, and He/N-2 and He/Ar ratios) of fumarolic gases from Nisyros, Greece, indicate that both arc-type magmatic water and local seawater feed the hydrothermal system. Isotopic composition of the deep fluid is estimated to be +4.9+/-0.5parts per thousand for delta(18)O and -11+/-5parts per thousand for deltaD corresponding to a magmatic water fraction of 0.7. Interpretation of the stable water isotopes was based on liquid-vapor separation conditions obtained through gas geothermometry. The H-2-Ar, H-2-N-2, and H-2-H2O geothermometers suggest reservoir temperatures of 345+/-15 degreesC, in agreement with temperatures measured in deep geothermal wells, whereas a vapor/liquid separation temperature of 260+/-30 degreesC is indicated by gas equilibria in the H2O-H-2-CO2-CO-CH4 system. The largest magmatic inputs seem to occur below the Stephanos-Polybotes Micros crater, whereas the marginal fumarolic areas of Phlegeton-Polybotes Megalos craters receive a smaller contribution of magmatic gases.

El depósito volcanogénico de Cu-Zn-Pb-Au El Cobre, Cuba Oriental: estructura y mineralogía

El depósito "El Cobre" encaja en series volcanosedimentarias del arco de islas Paleógeno de Cuba Oriental. Los tipos de mineralización presentes son estratiforme, estratoligada y filoniana. Los depósitos estratiformes (anhidrita, barita o criptomelana) se formaron por precipitación de fluidos hidrotermales en un fondo submarino. La mineralización filoniana constituye la zona de raíz de los depósitos estratiformes, presentando morfología de stockwork en sectores superficiales, pero hacia los sectores profundos disminuye la anchura de la zona mineralizada. En ella predominan pirita y calcopirita en asociación con cuarzo, excepto en algunos filones más superficiales en donde predomina esfalerita asociada a anhidrita. Las alteraciones asociadas son cloritización y silicificación en los sectoresmás profundos; en los más superficiales se dan sericitización y caolinitización. La mineralización estratoligada es contemporánea con la mineralización filoniana y se produce por reemplazamiento de diversos niveles favo r a bles de la serie: calizas en la parte inferior y el nivel exhalativo de sulfato en la parte superior (donde se producen las concentraciones de oro más elevadas del depósito). En base a los tipos de mineralización presentes y de acuerdo al carácter singenético de los depósitos con las rocas vulcanosedimentarias de la secuencia media del Grupo El Cobre, se propone un modelo de formación exhalativo vulcanogénico en un contexto geodinámico de arco de islas, similar al modelo Kuroko.

Mineral inventory of continuously erupting basaltic andesites at Arenal volcano, Costa Rica: implications for interpreting monotonous, crystal-rich, mafic arc stratigraphies

Except for the first 2 years since July 29, 1968, Arenal volcano has continuously erupted compositionally monotonous and phenocryst-rich (similar to35%) basaltic andesites composed of plagioclase (plag), orthopyroxene (opx), clinopyroxene (cpx), spinel olivine. Detailed textural and compositional analyses of phenocrysts, mineral inclusions, and microlites reveal comparable complexities in any given sample and identify mineral components that require a minimum of four crystallization environments. We suggest three distinct crystallization environments crystallized low Mg# (<78) silicate phases from andesitic magma but at different physical conditions, such as variable pressure of crystallization and water conditions. The dominant environment, i.e., the one which accounts for the majority of minerals and overprinted all other assemblages near rims of phenocrysts, cocrystallized clinopyroxene (Mg# similar to71-78), orthopyroxene (Mg# similar to71-78), titanomagnetite and plagioclase (An(60) to An(85)). The second environment cocrystallized clinopyroxene (Mg# 71-78), olivine (<Fo(78)), titanomagnetite, and very high An (similar to90) plagioclase, while the third cocrystallized clinopyroxene (Mg# 71-78) with high (>7) Al/Ti and high (>4 wt.%) Al2O3, titanomagnetite with considerable Al2O3 (10-18 wt.%) and possibly olivine but appears to lack plagioclase. A fourth crystallization environment is characterized by clinopyroxene (e.g., Mg#=similar to78-85; Cr2O3=0.15-0.7 wt.%), Al-, Cr-rich spinel olivine (similar toFo(80)), and in some circumstances high-An (>80) plagioclase. This assemblage seems to record mafic inputs into the Arenal system and crystallization at high to low pressures. Single crystals cannot be completely classified as xenocrysts, antecrysts (cognate crystals), or phenocrysts, because they often contain different parts each representing a different crystallization environment and thus belong to different categories. Bulk compositions are mostly too mafic to have crystallized the bulk of ferromagnesian minerals and thus likely do not represent liquid compositions. On the other hand, they are the cumulative products of multiple mixing events assembling melts and minerals from a variety of sources. The driving force for this multistage mixing evolution to generate erupting basaltic andesites is thought to be the ascent of mafic magma from lower crustal levels to subvolcanic depths which at the same time may also go through compositional modification by fractionation and assimilation of country rocks. Thus, mafic magmas become basaltic andesite through mixing, fractionation and assimilation by the time they arrive at subvolcanic depths. We infer new increments of basaltic andesite are supplied nearly continuously to the subvolcanic reservoir concurrently to the current eruption and that these new increments are blended into the residing, subvolcanic magma. Thus, the compositional monotony is mostly the product of repetitious production of very similar basaltic andesite. Furthermore, we propose that this quasi-constant supply of small increments of magma is the fundamental cause for small-scale, decade-long continuous volcanic activity; that is, the current eruption of Arenal is flux-controlled by inputs of mantle magmas. (C) 2004 Elsevier B.V. All rights reserved.

Geochemistry of the thermal springs and fumaroles of Basse-Terre Island, Guadeloupe, Lesser Antilles

The purpose of this work was to study jointly the volcanic-hydrothermal system of the high-risk volcano La Soufriere, in the southern part of Basse-Terre, and the geothermal area of Bouillante, on its western coast, to derive an all-embracing and coherent conceptual geochemical model that provides the necessary basis for adequate volcanic surveillance and further geothermal exploration. The active andesitic dome of La Soufriere has erupted eight times since 1660, most recently in 1976-1977. All these historic eruptions have been phreatic. High-salinity, Na-CI geothermal liquids circulate in the Bouillante geothermal reservoir, at temperatures close to 250 degrees C. These Na-CI solutions rise toward the surface, undergo boiling and mixing with groundwater and/or seawater, and feed most Na-CI thermal springs in the central Bouillante area. The Na-Cl thermal springs are surrounded by Na-HCO3 thermal springs and by the Na-Cl thermal spring of Anse a la Barque (a groundwater slightly mixed with seawater), which are all heated through conductive transfer. The two main fumarolic fields of La Soufriere area discharge vapors formed through boiling of hydrothermal aqueous solutions at temperatures of 190-215 degrees C below the ``Ty'' fault area and close to 260 degrees C below the dome summit. The boiling liquid producing the vapors of the Ty fault area has SD and delta(18)O values relatively similar to those of the Na-CI liquids of the Bouillante geothermal reservoir, whereas the liquid originating the vapors of the summit fumaroles is strongly enriched in O-18, due to input of magmatic fluids from below. This process is also responsible for the paucity of CH;I in the fumaroles. The thermal features around La Soufriere dome include: (a) Ca-SO4 springs, produced through absorption of hydrothermal vapors in shallow groundwaters; (b) conductively heated, Ca-Na-HCO3 springs; and (c) two Ca-Na-Cl springs produced through mixing of shallow Ca-SO4 waters and deep Na-Cl hydrothermal liquids. The geographical distribution of the different thermal features of La Soufriere area indicates the presence of: (a) a central zone dominated by the ascent of steam, which either discharges at the surface in the fumarolic fields or is absorbed in shallow groundwaters; and (b) an outer zone, where the shallow groundwaters are heated through conduction or addition of Na-Cl liquids coming from hydrothermal aquifer(s).