Non-Newtonian rheology of crystal-bearing magmas and implications for magma ascent dynamics

The eruptive dynamics of volcanic systems are largely controlled by the viscosity of deforming magma. Here we report the results of a series of high-temperature, high-pressure experiments at conditions relevant for volcanic conduits (250 MPa confining pressure and temperature between 500 °C and 900 °C) that were undertaken to investigate the rheology of magma with crystal fractions varying between 0.5 and 0.8 (50 to 80 wt.%) at different strain-rate conditions. The experiments demonstrate that the presence of crystals increases the relative viscosity (ratio between the viscosity of the mixture and the viscosity of the melt phase) of magmas and additionally induces a decrease of the relative viscosity with increasing strain-rate (shear thinning, non-Newtonian behavior). The experimental results, combined with existing data at low crystal fractions (0–0.3), were used to develop a semi-empirical parameterization that describes the variations of relative viscosity for crystal fractions between 0 and 0.8 and accounts for the complex non-Newtonian rheology of crystal-bearing magmas. The new parameterization, included into numerical models simulating the magma ascent dynamics, reveals that strain-rate-dependent rheology significantly modifies the dynamic behavior inside volcanic conduits, particularly affecting the magma fragmentation conditions

Multidisciplinary approach to constrain kinematics of fault zones at shallow depths: a case study from the Cameros–Demanda thrust (North Spain)

Determining transport direction in thrusts is one of the main issues to study deformation and to achieve reliable balanced cross-sections reconstruction of inverted basins and thrust systems. Anisotropy of Magnetic Susceptibility provides a tool to study deformation in fault rocks through the relationships between the magnetic susceptibility ellipsoid and deformational structures. The Cameros-Demanda Thrust (N Spain) shows a relatively simple history, with an average N-directed movement during the Cenozoic, and provides the possibility of determining the relationships between magnetic fabrics and transport directions in thrusts. The outcropping rocks are Mesozoic limestones in the hangingwall of the main thrust and Cenozoic conglomerates and Albian sandstones and coal in its footwall. Illite crystallinity and organic matter maturity indicate P-T conditions on the order of....The study of AMS in fault rocks (more than 400 samples distributed in 25 sites with fault gouge, breccia and microbreccia) in the Cameros thrust and its comparison with kinematic indicators (foliation, S/C structures and slickenside striations) indicates, in spite of the a priori simple relationships inferred from thrust geometry, a complex history of movements, changing from top-to-the-NW to top-to-the-NE along the history of Cenozoic thrusting. The transport direction is either oblique to the magnetic lineation and perpendicular to the strike of magnetic foliation, and can be checked with shear structures observed in thin sections and other kinematic indicators. The results obtained indicate that AMS can give clues about the transport direction in thrusts depending on the particular structures developed in each studied area

Preliminary data on the structure and potential of the Tocomar geothermal field (Puna plateau, Argentina).

AbstractThis study presents new stratigraphic, structural and hydrogeological data on the Tocomar geothermal volcanic area (Puna plateau, Central Andes, NW Argentina), together with preliminary geochemical and magnetotelluric data.The main geothermal reservoir is located within the fractured Pre-Palaeozoic–Ordovician units. The reservoir is recharged by meteoric waters. Geothermal fluids upwell where main regional structures intersect secondary structures associated with the development of the Tocomar basin. Preliminary data indicate a reservoir temperature of ∼ 200° C and a local geothermal gradient of ∼ 130° C/km associated with the Quaternary volcanic activity in the Tocomar area

Electrical conductivity during incipient melting in the oceanic low-velocity zone

International audienceThe low-viscosity layer in the upper mantle, the asthenosphere, is a requirement for plate tectonics1. The seismic low velocities and the high electrical conductivities of the asthenosphere are attributed either to subsolidus, water-related defects in olivine minerals2, 3, 4 or to a few volume per cent of partial melt5, 6, 7, 8, but these two interpretations have two shortcomings. First, the amount of water stored in olivine is not expected to be higher than 50 parts per million owing to partitioning with other mantle phases9 (including pargasite amphibole at moderate temperatures10) and partial melting at high temperatures9. Second, elevated melt volume fractions are impeded by the temperatures prevailing in the asthenosphere, which are too low, and by the melt mobility, which is high and can lead to gravitational segregation11, 12. Here we determine the electrical conductivity of carbon-dioxide-rich and water-rich melts, typically produced at the onset of mantle melting. Electrical conductivity increases modestly with moderate amounts of water and carbon dioxide, but it increases drastically once the carbon dioxide content exceeds six weight per cent in the melt. Incipient melts, long-expected to prevail in the asthenosphere10, 13, 14, 15, can therefore produce high electrical conductivities there. Taking into account variable degrees of depletion of the mantle in water and carbon dioxide, and their effect on the petrology of incipient melting, we calculated conductivity profiles across the asthenosphere for various tectonic plate ages. Several electrical discontinuities are predicted and match geophysical observations in a consistent petrological and geochemical framework. In moderately aged plates (more than five million years old), incipient melts probably trigger both the seismic low velocities and the high electrical conductivities in the upper part of the asthenosphere, whereas in young plates4, where seamount volcanism occurs6, a higher degree of melting is expected

Unexpected large eruptions from buoyant magma bodies within viscoelastic crust

Large volume effusive eruptions with relatively minor observed precursory signals are at odds with widely used models to interpret volcano deformation. Here we propose a new modelling framework that resolves this discrepancy by accounting for magma buoyancy, viscoelastic crustal properties, and sustained magma channels. At low magma accumulation rates, the stability of deep magma bodies is governed by the magma-host rock density contrast and the magma body thickness. During eruptions, inelastic processes including magma mush erosion and thermal effects, can form a sustained channel that supports magma flow, driven by the pressure difference between the magma body and surface vents. At failure onset, it may be difficult to forecast the final eruption volume; pressure in a magma body may drop well below the lithostatic load, create under-pressure and initiate a caldera collapse, despite only modest precursors

Amphibole and apatite insights into the evolution and mass balance of Cl and S in magmas associated with porphyry copper deposits

Chlorine and sulfur are of paramount importance for supporting the transport and deposition of ore metals at magmatic–hydrothermal systems such as the Coroccohuayco Fe–Cu–Au porphyry–skarn deposit, Peru. Here, we used recent partitioning models to determine the Cl and S concentration of the melts from the Coroccohuayco magmatic suite using apatite and amphibole chemical analyses. The pre-mineralization gabbrodiorite complex hosts S-poor apatite, while the syn- and post-ore dacitic porphyries host S-rich apatite. Our apatite data on the Coroccohuayco magmatic suite are consistent with an increasing oxygen fugacity (from the gabbrodiorite complex to the porphyries) causing the dominant sulfur species to shift from S2− to S6+ at upper crustal pressure where the magmas were emplaced. We suggest that this change in sulfur speciation could have favored S degassing, rather than its sequestration in magmatic sulfides. Using available partitioning models for apatite from the porphyries, pre-degassing S melt concentration was 20–200 ppm. Estimates of absolute magmatic Cl concentrations using amphibole and apatite gave highly contrasting results. Cl melt concentrations obtained from apatite (0.60 wt% for the gabbrodiorite complex; 0.2–0.3 wt% for the porphyries) seems much more reasonable than those obtained from amphibole which are very low (0.37 wt% for the gabbrodiorite complex; 0.10 wt% for the porphyries). In turn, relative variations of the Cl melt concentrations obtained from amphibole during magma cooling are compatible with previous petrological constraints on the Coroccohuayco magmatic suite. This confirms that the gabbrodioritic magma was initially fluid undersaturated upon emplacement, and that magmatic fluid exsolution of the gabbrodiorite and the pluton rooting the porphyry stocks and dikes were emplaced and degassed at 100–200 MPa. Finally, mass balance constraints on S, Cu and Cl were used to estimate the minimum volume of magma required to form the Coroccohuayco deposit. These three estimates are remarkably consistent among each other (ca. 100 km3) and suggest that the Cl melt concentration is at least as critical as that of Cu and S to form an economic mineralization

Controls on explosive-effusive volcanic eruption styles

One of the biggest challenges in volcanic hazard assessment is to understand how and why eruptive style changes within the same eruptive period or even from one eruption to the next at a given volcano. This review evaluates the competing processes that lead to explosive and effusive eruptions of silicic magmas. Eruptive style depends on a set of feedbacks involving interrelated magmatic properties and processes. Foremost of these are magma viscosity, gas loss, and external properties such as conduit geometry. Ultimately, these parameters control the speed at which magmas ascend, decompress and outgas en route to the surface, and thus determine eruptive style and evolution