Lava Rheology Research Articles

Real-time, in situ viscosity mapping of active lava

Viscosity is a fundamental physical property that controls lava flow dynamics, runout distance, and velocity, which are critical factors in assessing and mitigating risks associated with effusive eruptions. Natural lava viscosity is driven by a dynamic interplay among melt, crystals, and bubbles in response to the emplacement conditions. These conditions are challenging to replicate in laboratory experiments, yet this remains the most common method for quantifying lava rheology. Few in situ viscosity measurements exist, but none of those constrains the spatial evolution of viscosity along an entire active lava flow field. Here, we present the first real-time, in situ viscosity map of active lava as measured in the field at Litli-Hrútur, Iceland. We precisely measured a lava viscosity increase of over two orders of magnitude, associated with a temperature decrease, crystallinity increase, and vesicularity decrease from near-vent to distal locations, crossing the pāhoehoe−‘a‘ā transition. Our data expand the limited database of three-phase lava viscosity, which is crucial for improvements and validation of the current numerical, experimental, and petrological approaches used to estimate lava viscosity. Further, this study showcases that field viscometry provides a rapid, accurate, and precise assessment of lava viscosity that can be implemented in eruptive response modeling of lava transport.

The viscous-brittle transition in flowing crystal-bearing volcanic dome lavas

The ascent and advance of volcanic dome lava is non-linear and viscoelastic. There exists a mismatch between current theoretical approaches to dome lava rheology, which are based on rheological laws for viscous suspensions, and empirical experimental approaches to convolved viscous-brittle deformation, which show mixed evidence for simultaneous lava flow and fracturing. The missing requirement is a unified framework for understanding the transition between micro-mechanical flow mechanisms that are dominantly viscous, and those that include micro-cracking in multiphase suspensions such as magmas. Here, we use high-temperature compression rheology with sample-scale acoustic emission analysis to constrain the conditions under which crystal-rich volcanic dome lava can flow by mixed viscous and brittle fracturing processes at small scales, leading to ‘crackling’ acoustic signals, even at moderate shear stresses extant in nature. Using multi-directional permeability measurements on large 60 mm diameter quenched samples of natural magmas, we show that this micro-cracking flow mechanism leads to permeability anisotropy, localizing outgassing into pathways that are off-axis relative to the direction of flow. Finally, we use a scaling approach and a database of published observations from real eruptions to upscale our findings, and show that bulk, apparently ductile flow of low-porosity dome magma is likely to involve a local mixed-mode of micro-cracking and viscous flow during the shallowest portions of ascent and during emplacement on the Earth's surface. The micro-cracking involved in lava advance divorces real crystal-bearing lava emplacement from most current rheology models based on a purely viscous micro-mechanism and shows that a revised solution for the rheology of mixed brittle-viscous flow is required. By re-examining published numerical models for dome emplacement, we demonstrate that the viscous-brittle transition can be intercepted in spatially heterogeneous zones within the dome core.

Estimates of Crystallinity Utilizing Differential Scanning Calorimetry: Application to the Kīlauea 2018 Lower East Rift Zone Eruption

Abstract Rocks produced by diverse processes, from condensation in space to impacts on planetary surfaces to volcanism, contain both crystals and amorphous material. Crystallinity provides information on the thermal history of the sample and is especially important in characterizing volcanic rocks and pyroclasts because lava rheology is profoundly influenced by the crystal content. Crystallinity is typically quantified via microscopy, using transmitted light or backscattered electrons. However, many samples present visibly ambiguous textures such as intimate intergrowth of crystal phases, and/or crystal sizes extending down to the nanometer scale. Here, we apply calorimetric methods involving heat capacity and enthalpy to assess the crystallinity of a series of volcanic samples. We tested three different approaches, using differential scanning calorimetry, on 30–40 mg aliquots of powdered basalts from the 2018 Kīlauea lower East Rift Zone. The first approach involves determining the magnitude of the increase in heat capacity at the glass transition, which can determine crystallinity to a 1σ precision of ±3%. The second approach is based on the enthalpy of fusion, which requires a longer more complex procedure with results that are typically more uncertain than for the heat capacity method, with a 1σ of ±6%. A final method utilizing differences in enthalpies calculated from the heat capacities required the most complex procedure and has the greatest uncertainty of ±18%. Preliminary results for lavas with microscopically determined crystallinities ranging from 11 to 98% indicate that crystallinity based on calorimetric data can be tens of percent higher than the average value identified using microscopy and petrographic analysis. Image-based methodologies applied to sections of samples reveal spatial heterogeneity and details in texture and crystallinity, whereas calorimetry-based methodologies capture the overall ‘bulk sample’ properties, unbiased by section effects or imaging resolution limits. These techniques are a powerful combination that can present complementary views of crystallinity.

The Control of Lava Rheology on the Formation of Lunar Sinuous Rilles by Substrate Thermal Erosion: Topographic and Morphometric Relationships with Eruption Rates, Erosion Rates, Event Durations, and Erupted Magma Volumes

We describe a model of the fluid mechanics and thermodynamics of turbulent lava flows capable of thermally eroding sinuous rille channels on bodies without atmospheres. The model assumes Bingham plastic rheology for the lava and shows how the effects of radiant cooling and consequent crystallization change the rheology and control the point at which turbulence ceases. A correlation is found between magma volume eruption rate and the length and width of the eroded rille channel. Thus, simple measurements of rille length and width in images can provide reliable estimates of magma eruption rates. The model also predicts rille floor erosion rates, so that if rille depths are measured, eruption durations and erupted magma volumes can also be found. The model is applied in detail to six well-preserved lunar rilles and more generally to a published catalog of 214 lunar rilles. We find that rille-forming eruptions have magma volume eruption rates of a few times 104 m3 s−1, durations of up to 3 months, and erupted magma volumes up to ∼200 km3, consistent with theoretical predictions of basaltic magma ascent and eruption from deep mantle sources. The key requirement for rille formation, rather than mare lava flow deposit formation, is the turbulent eruption of a sufficiently large volume of low-viscosity magma at a sufficiently low eruption rate.

Magnetic Resonance Imaging of Multi‐Phase Lava Flow Analogs: Velocity and Rheology

AbstractThe rheology of lavas and magmas exerts a strong control on the dynamics and hazards posed by volcanic eruptions. Magmas and lavas are complex mixtures of silicate melt, suspended crystals, and gas bubbles. To improve the understanding of the dynamics and effective rheology of magmas and lavas, we performed dam‐break flow experiments using suspensions of silicone oil, sesame seeds, and N2O bubbles. Experiments were run inside a magnetic resonance imaging (MRI) scanner to provide imaging of the flow interior. We varied the volume fraction of sesame seeds between 0 and 0.48, and of bubbles between 0 and 0.21. MRI phase‐contrast velocimetry was used to measure liquid velocity. We fit an effective viscosity to the velocity data by approximating the stress using lubrication theory and the imaged shape of the free surface. In experiments with both particles and bubbles (three‐phase suspensions), we observed shear banding in which particle‐poor regions deform with a lower effective viscosity and dominate flow propagation speed. Our observations demonstrate the importance of considering variations in phase distributions within magmatic fluids and their implications on the dynamics of volcanic eruptions.

Experimental Constraints on the Rheology of Lavas From 2021 Cumbre Vieja Eruption (La Palma, Spain)

AbstractThe 2021 Tajogaite eruption of Cumbre Vieja (La Palma, Spain) was typified by the emission of low viscosity lavas that flowed at high velocities and inundated a large area. We experimentally investigated the rheological evolution of melt feeding the eruption through concentric cylinder viscometry to understand the exceptional flowing ability of these lavas and constrain its emplacement dynamics. We conducted a set of cooling deformation experiments at different cooling rates (from 0.1 to 10 °C/min), and isothermal deformation experiments at subliquidus dwell temperatures between 1225 and 1175°C. All experiments were conducted at a shear rate of 10 s−1. Results show that disequilibrium crystallization and its timescale fundamentally control the rheological evolution of the melt, resulting in different rheological response to deformation of the crystal‐bearing magmatic suspension. Integrating rheological data with field observations allows us to shed light on the mechanisms that govern the high flowability of these lavas.

The impact of vent geometry on the growth of lava domes

SUMMARY Thick lava flows that are a feature of many volcanic fields on the Earth and Venus vary from sheet-like to nearly perfect axisymmetric domes. Here, we investigate how these geometrical characteristics depend on the shape of the feeder vent. We study the gravitational spreading of viscous lava erupting from elliptical vents onto a flat surface using 3-D numerical models. The aspect ratio of the vent, defined to be the major to minor axes ratio, varies between 1 and 25. In the limit of an aspect ratio of one, the vent is circular and spreading is axisymmetric. In the limit where the ratio is large, the vent behaves as a fissure. The numerical models rely on an isoviscous lava rheology and a constant volumetric eruption rate. In all cases, the initial phase of the dome’s evolution is in a lava-discharge dominated regime such that spreading is insignificant and the height of the dome increases at a constant rate over the vent area. For vent aspect ratios greater than five, three successive regimes of spreading are identified: 2-D spreading in the direction perpendicular to the major axis of the vent, a transient phase such that the dome shape evolves towards that of a circular dome and a late axisymmetric spreading phase that does not depend on the vent shape. These regimes are delimited by the times required for the flow thickness above the vent to reach a given height and for the flow to spread axisymmetrically up to a length equal to the semi-major axis of the vent. Numerical results for the flow height and runout length tend towards the similarity solutions in the 2-D and axisymmetric regimes. Two main implications for highly viscous (rhyolitic) fissure eruptions can be drawn. First, the fissure length determines the flow regimes. The longer the vent fissure length, the longer the early lava discharge regime and 2-D spreading perpendicular to the length of the fissure. Second, the aspect ratio of fissure-fed lava flows can be used as an indicator of the fissure length and the duration of lava discharge. The ellipticity of some terrestrial fissure-fed flows provides evidence for viscous gravity-driven spreading terminated before the onset of the axisymmetric regime. On the other hand, the circular domes on Venus appear to be the result of fissure-fed eruptions sustained enough for the spreading to reach the axisymmetric regime. We propose relationships providing estimates of the fissure length and the duration of lava discharge based on fossil dome dimensions.

Rheology of three-phase suspensions determined via dam-break experiments.

Three-phase suspensions, of liquid that suspends dispersed solid particles and gas bubbles, are common in both natural and industrial settings. Their rheology is poorly constrained, particularly for high total suspended fractions (≳0.5). We use a dam-break consistometer to characterize the rheology of suspensions of (Newtonian) corn syrup, plastic particles and CO2 bubbles. The study is motivated by a desire to understand the rheology of magma and lava. Our experiments are scaled to the volcanic system: they are conducted in the non-Brownian, non-inertial regime; bubble capillary number is varied across unity; and bubble and particle fractions are 0 ≤ ϕ gas ≤ 0.82 and 0 ≤ ϕ solid ≤ 0.37, respectively. We measure flow-front velocity and invert for a Herschel-Bulkley rheology model as a function of , , and the capillary number. We find a stronger increase in relative viscosity with increasing in the low to intermediate capillary number regime than predicted by existing theory, and find both shear-thinning and shear-thickening effects, depending on the capillary number. We apply our model to the existing community code for lava flow emplacement, PyFLOWGO, and predict increased viscosity and decreased velocity compared with current rheological models, suggesting existing models may not adequately account for the role of bubbles in stiffening lavas.

The effect of bubbles on the rheology of basaltic lava flows: Insights from large-scale two-phase experiments

Bubbles are known to have a dramatic effect on lava rheology, but their impact on lava flow run out distance has not been extensively investigated. To quantify bubble effects on basaltic lava flow emplacement, we conducted a series of large-scale two-phase (melt+bubbles) experiments at the Syracuse University Lava Project facility. In this study, nine basaltic lava flows with low to intermediate vesicle contents (approximately 3-24 vol%) were emplaced over gently sloping (4°-9°) dry sand surfaces. Their velocities were measured by particle image velocimetry, and their bulk viscosities were calculated using Jeffreys' equation. Comparison with melt-only viscosities (calculated through the GRD model) allowed the isolation of the effect of bubbles on the viscosity of the experimental flows. We find that within the investigated range of vesicularities, bubbles increase the relative viscosity of the experimental flows by up to an order of magnitude. The increase appears to be steeper for sheet flows than for lobate and toey flows. Our results, which quantitatively assess the effect of bubbles on the rheology of basaltic lava flows, will provide a basis for the refinement of models for multi-phase flows that are fundamental to accurately forecasting the emplacement of lava flows.

Remote Characterization of Dominant Wavelengths From Surface Folding on Lava Flows Using Lidar and Discrete Fourier Transform Analyses

AbstractSurface folding is common in lava flows of all compositions and is believed to be due to changes in viscosity and flow velocity between the cooling crust and the more fluid flow interior. However, our understanding of the relationship between surface folding and flow rheology is incomplete. In this study we analyze digital terrain models of eight lava flows ranging in composition from basaltic andesite to rhyolite using a discrete Fourier transform analysis to quantitatively determine dominant surface fold wavelengths. Our discrete Fourier transform analyses show that each lava flow has multiple fold generations and that dominant wavelengths are more closely related to calculated effective viscosity than to lava composition. At our Oregon sites, average dominant wavelengths generally increase with viscosity (r2=0.68), and the correlation improves (r2=0.87) when expanded by including previously measured fold wavelengths and viscosities from the global database. However, there are a few exceptions to this positive trend where a few lava flows have lower or higher than expected dominant fold wavelengths, which we infer are due to secondary factors such as differences in eruption conditions (eruption rate, temperature, etc.). Additionally, over a 5 order of magnitude range in viscosity, there is significant overlap between the ranges of fold wavelengths, particularly from 10 to 20m, for lavas from basaltic andesite to rhyolite, making it difficult to determine a numeric correlation between surface folds and lava rheology that would allow remote characterization of lava.

Far-field boundary conditions in channeled lava flow with viscous dissipation

Cooling and dynamics of lava flowing in a rectangular channel driven by the gravity force is numerically modeled. The purpose is to evaluate the thermal process as a function of time involving the liquid lava in contact with the solid boundary that flanks lava. Lava rheology is dependent on temperature and strain rate according to a power law function. The model couples dynamics and thermodynamics inside the lava channel and describes the thermal evolution of the solid boundary enclosing the channel. Numerical tests indicate that the solution of the thermo-dynamical problem is independent of the mesh. The boundary condition at the ground and at the levees is treated assuming a solid boundary around the lava flow across which lava can exchange heat by conduction. A far field thermal boundary condition allows to overcome the assumption of constant temperature or constant heat flow as boundary conditions, providing more realistic results. The effect of viscous heating is evaluated and discussed .

Multiple-generation folding and non-coaxial strain of lava crusts

Viscoelastic strain in lava flows is commonly expressed as gravity-driven buckling of the lava crust. This surface folding process creates the well-known ropy pāhoehoe texture of basaltic lavas and the ogives and surface ridges of more compositionally evolved lava flows. Previous work has shown that surface fold wavelengths are proportional to the viscosity contrast between the lava crust and core, and to the thickness of the crust. Thus, fold analysis can be an important tool for understanding lava flow rheology. We analyze fold wavelength patterns of solidified natural lava flows from the Myvatn lava fields (Iceland), Piton de La Fournaise (La Reunion), and experimental lava flows from the Syracuse University Lava Project. In each case, lava flows exhibited two dominant wavelengths, consistent with multiple generations of coaxial folding. The ratio of the two dominant wavelengths for basalt (Iceland, La Reunion) is ~ 5:1 whereas the wavelength ratio for basaltic andesite (Syracuse) is ~ 3:1, suggesting a compositional control on deformation, as proposed by previous studies. Video analysis of incrementally folded Syracuse lava crusts reveals significant non-coaxial strain, which violates the assumptions of plane strain used in crustal buckling models. These results show that interpreting lava rheology from finite strain requires careful consideration of complex three-dimensional strain fields. Despite these complexities, the correlation between fold wavelength ratios and lava flow composition persists and may provide important insight into flow characterization.

Failure criteria for porous dome rocks and lavas: a study of Mt. Unzen, Japan

Abstract. The strength and macroscopic deformation mode (brittle vs. ductile) of rocks is generally related to the porosity and pressure conditions, with occasional considerations of strain rate. At high temperature, molten rocks abide by Maxwell's viscoelasticity and their deformation mode is generally defined by strain rate or reciprocally by comparing the relaxation timescale of the material (for a given condition) to the observation timescale – a dimensionless ratio known as the Deborah (De) number. Volcanic materials are extremely heterogeneous, with variable concentrations of crystals, glass–melt, and vesicles (of different sizes), and a complete description of the conditions leading to flow or rupture as a function of temperature, stress and strain rate (or timescale of observation) eludes us. Here, we examined the conditions which lead to the macroscopic failure of variably vesicular (0.09–0.35), crystal-rich (∼ 75 vol %), pristine and altered dome rocks (at ambient temperature) and lavas (at 900 °C) from Mt. Unzen volcano, Japan. We found that the strength of the dome rocks decreases with porosity and is commonly independent of strain rate; when comparing pristine and altered rocks, we found that the precipitation of secondary mineral phases in the original pore space caused minor strengthening. The strength of the lavas (at 900 °C) also decreases with porosity. Importantly, the results demonstrate that these dome rocks are weaker at ambient temperatures than when heated and deformed at 900 °C (for a given strain rate resulting in brittle behaviour). Thermal stressing (by heating and cooling a rock up to 900 °C at a rate of 4 °C min−1, before testing its strength at ambient temperature) was found not to affect the strength of rocks.In the magmatic state (900 °C), the rheology of the dome lavas is strongly strain rate dependent. Under conditions of low experimental strain rate (≤ 10−4 s−1), ductile deformation dominated (i.e. the material sustained substantial, pervasive deformation) and displayed a non-Newtonian shear thinning behaviour. In this regime, the apparent viscosities of the dome lavas were found to be essentially equivalent, independent of vesicularity, likely due to the lack of pore pressurisation and efficient pore collapse during shear. At high experimental strain rates ( ≥ 10−4 s−1) the lavas displayed an increasingly brittle response (i.e. deformation resulted in failure along localised faults); we observed an increase in strength and a decrease in strain to failure as a function of strain rate. To constrain the conditions leading to failure of the lavas, we analysed and compared the critical Deborah number at failure (Dec) of these lavas to that of pure melt (Demelt = 10−3–10−2; Webb and Dingwell, 1990). We found that the presence of crystals decreases Dec to between 6.6×10−4 and 1×10−4. The vesicularity (φ), which dictates the strength of lavas, further controls Dec following a linear trend. We discuss the implications of these findings for the case of magma ascent and lava dome structural stability.

Equilibrium Viscosity and Disequilibrium Rheology of a high Magnesium Basalt from Piton De La Fournaise volcano, La Reunion, Indian Ocean, France

Lava flows are a common hazard at basaltic to intermediate volcanoes in general and have posed a significant threat to La Reunion Island specifically over the past centuries. In sustained flow units, the efficiency of lava transport away from the vent is dominated by cooling. For basaltic to intermediate lavas it is the ability of the lava to solidify during cooling which exerts a first-order control on the extent of flow. As a consequence, understanding the sub-liquidus rheology of lavas has become a key focus in lava flow research in the past decade. To date, a systematic understanding of lava rheology during emplacement conditions has been largely hampered by a lack of experimental data. Here we present new data on the rheological evolution of a crystalizing high Mg basalt from Piton de la Fournaise. Sub-liquidus experiments were performed at constant cooling rates of 0.5 to 5 K/min, to mimic conditions similar to those experienced by lava during flow on the surface and magma during dike and sill emplacement. The data show that the effective viscosity of the crystallizing suspension increases until reaching a specific sub-liquidus temperature, the so-called “rheological cut-off temperature” (T cutoff ), at which the lava effectively solidifies. This departure to high viscosity is a consequence of rapid crystallization and is found to be primarily controlled by the imposed cooling-rate. Based on our experimental data, we adapt the failure forecasting method (FFM), which is commonly used to describe the self-accelerating nature of seismic signals to forecast material failure, to predict the rheological cut-off temperature (T cutoff ). The data presented serve to expand the modest existing experimental database on non-equilibrium rheology of lavas as a step towards understanding the underlying process dynamics.

Textural, thermal, and topographic constraints on lava flow system structure: the December 2010 eruption of Piton de la Fournaise

In this study, we examine the channel-fed ‘a‘ā lava flow system that was emplaced during a very short (less than 15 h long) eruption at Piton de la Fournaise (La Reunion) in December 2010. The system had four branches, the longest of which was 1100 m long. Three branches were emplaced over a smooth-surfaced pāhoehoe flow field with a vertical relief of 1–2 m and did not undergo burial by subsequent events. The fourth branch erupted from the same vent as the 1957 eruption and re-used the pre-existing channels of that eruption. In the proximal–medial sections of the three systems that were unconfined, we identified channelized flow sections that were characterized by the presence of either a single channel or multiple braided channels. These fed short (30–260 m long) zones of dispersed flow in the distal sections. We subsequently investigated the role of lava rheology (as controlled by downflow variations in crystal and bubble content) and pre-existing topography in triggering the transitions between single-channel and braided channel flow sections. Crystal content was 10 to 70 vol% and vesicle content was 18 to 55 vol%; cooling rates over distance (derived from glass chemistry) were 11 °C/km to 27 °C/km. However, downflow textural and thermal evolution appeared to neither affect, nor be affected by, whether the channel was single or braided. Instead, the channel network architecture could be related to even modest underlying slope variations. Here, a slope increase resulted in channel confluence, and a slope decrease resulted in channel bifurcation. This process was reversible, in that downflow slope variation could drive the channel network architecture to switch back and forth between a single channel and multiple braided channels several times along its length. Dispersed flow is always present immediately behind the flow front, irrespective of underlying topography. Three previous studies of basaltic lava flows found that steeper slopes favored braided channels, the opposite of what was observed here. We suggest that the underlying substrate and lava type may exert a control on this behavior, but further studies remain necessary.

The viscosity of pāhoehoe lava: In situ syn-eruptive measurements from Kilauea, Hawaii

Viscosity is one of the most important physical properties controlling lava flow dynamics. Usually, viscosity is measured in the laboratory where key parameters can be controlled but can never reproduce the natural environment and original state of the lava in terms of crystal and bubble contents, dissolved volatiles, and oxygen fugacity. The most promising approach for quantifying the rheology of molten lava in its natural state is therefore to carry out direct field measurements by inserting a viscometer into the lava while it is flowing. Such in-situ syn-eruptive viscosity measurements are notoriously difficult to perform due to the lack of appropriate instrumentation and the difficulty of working on or near an active lava flow. In the field, rotational viscometer measurements are of particular value as they have the potential to measure the properties of the flow interior rather than an integration of the viscosity of the viscoelastic crust + flow interior. To our knowledge only one field rotational viscometer is available, but logistical constraints have meant that it has not been used for 20 yr. Here, we describe new viscosity measurements made using the refurbished version of this custom-built rotational viscometer, as performed on active pāhoehoe lobes from the 61G lava flow of Kilauea's Pu'u ‘Ō‘ō eruption in 2016. We successfully measured a viscosity of ∼380 Pa s at strain-rates between 1.6 and 5 s−1 and at 1144 °C. Additionally, synchronous lava sampling allowed us to provide detailed textural and chemical characterization of quenched samples. Application of current physico-chemical models based on this characterization (16 ± 4 vol.% crystals; 50 ± 6 vol.% vesicles), gave viscosity estimates that were approximately compatible with the measured values, highlighting the sensitivity of model-based viscosity estimates on the effect of deformable bubbles. Our measurements also agree on the range of viscosities in comparison to previous field experiments on Hawaiian lavas. Conversely, direct comparison with sub-liquidus rheological laboratory measurements on natural lavas was unsuccessful because recreating field conditions (in particular volatile and bubble content) is so far inaccessible in the laboratory. Our work shows the value of field rotational viscometry fully-integrated with sample characterization to quantify three-phase lava viscosity. Finally, this work suggests the need for the development of a more versatile instrument capable of recording precise measurements at low torque and low strain rate, and with synchronous temperature measurements.

The effect of inflation on the morphology-derived rheological parameters of lava flows and its implications for interpreting remote sensing data - A case study on the 2014/2015 eruption at Holuhraun, Iceland

Morphology-derived lava flow rheology is a frequently used tool in volcanology and planetary science to determine rheological parameters and deduce the composition of lavas on terrestrial planets and their moons. These calculations are usually based on physical equations incorporating 1) lava flow driving forces: gravity, slope and flow-rate and 2) morphological data such as lava flow geometry: flow-width, -height or shape of the flow outline. All available methods assume that no geometrical changes occur after emplacement and that the measured flow geometry reflects the lava's apparent viscosity and/or yield strength during emplacement. It is however well-established from terrestrial examples that lava flows may inflate significantly after the cessation of flow advance. This inflation affects, in turn, the width-to-height ratio upon which the rheological estimates are based and thus must result in uncertainties in the determination of flow rheology, as the flow height is one of the key parameters in the morphology-based deduction of flow properties. Previous studies have recognized this issue but, to date, no assessment of the magnitude of this error has been presented. This is likely due to a lack of digital elevation models (DEMs) at sufficiently high spatial and temporal resolution.The 2014/15 Holuhraun eruption in central Iceland represents one of the best monitored large volume (1.5 km3) lava flow fields (85 km2) to date. An abundance of scientific field and remote sensing data were collected during its emplacement. Moreover, inflation plays a key role in the emplacement dynamics of the late stage of the lava field. Here, we use a time series of high resolution DEMs acquired by the TanDEM-X satellite mission prior, during and after the eruption to evaluate the error associated with the most common methods of deriving lava flow rheology from morphological parameters used in planetary science.We can distinguish two dominant processes as sources of error in the determination of lava flow rheology from morphology 1) wholesale inflation of lava channels and 2) post halting inflation of individual lava toes. These result in a 2.4- to 17 - fold overestimation of apparent viscosity and a 0.7- to 2.4 - fold overestimation of yield strength. When applied in planetary sciences, this overestimation in rheological parameters translates directly to an overestimation of the respective lavas silica content. We conclude that, although qualitatively informative, morphological analysis is insufficient to discern lava rheology and composition. Instead, in-situ analysis together with high resolution remote sensing data is needed to properly constrain the compositions involved in planetary volcanism.

Emplacement dynamics and timescale of a Holocene flow from the Cima Volcanic Field (CA): Insights from rheology and morphology

We present a rheological and morphological study of a Holocene lava flow emitted by a monogenetic cinder cone in the Cima Volcanic Field, eastern California. Our field observations focused on surface morphology, which transitions from smooth core extrusions near the vent to jagged ‘a’ā blocks over the majority of the flow, and on channel and levée dimensions. We collected airborne photogrammetry data and used it to generate a digital elevation model. From this, the total flow volume was estimated and surface roughness was quantified in terms of standard deviation of the real surface (5cm resolution) from the software-generated 1m-average plane. Sample textural analyses revealed that the near-vent portion of the flow is significantly more crystalline (ϕxtal=0.95±0.04) than the main flow body (ϕxtal=0.66±0.11). The rheology of Cima lavas was determined experimentally by concentric cylinder viscometry between 1550°C and 1160°C, including the first subliquidus rheology measurements for a continental intraplate trachybasaltic lava. The experimentally determined effective viscosity increases from 54Pa·s to 1361Pa·s during cooling from the liquidus (~1230°C) to 1160°C, where crystal fraction is 0.11. The lava viscosity over this range is still lower than most basaltic melts, due to the high alkali content of Cima lavas (~6wt% Na2O+K2O). Monte Carlo simulations were used to account for and propagate experimental uncertainties, and to determine which rheological model (Bingham, power law, or Herschel-Bulkley) provides the best-fit of the obtained rheological data. Results suggest that Bingham and Herschel-Bulkley models are statistically indistinguishable from each other, and that both fit the data better than a power law model. By combining field observations and experimental results, we reconstructed the eruption temperature and few days-long emplacement history of the Cima flow.

The complex rheology of megacryst-rich magmas: The case of the mugearitic “cicirara” lavas of Mt. Etna volcano

Multiphase magma rheological properties play a fundamental role on lava flow transport, emplacement and morphology. To date, however, the three-phase (melt+crystals+bubbles) rheology of natural magma remains relatively understudied. We present here a series of high-temperature experiments designed to investigate the multiphase rheology of a mugearitic megacryst-bearing lava from Mt. Etna. A peculiar textural feature of this magma is the abundance of cm-size plagioclase crystals (megacrysts) together with smaller size crystals (phenocrysts and microlites), yielding a very wide crystal size distribution. We combined different experimental techniques (rotational and compressional rheometry) to investigate the rheology of this natural lava under different degrees of partial melting at subliquidus conditions. Results indicate that natural megacryst-bearing mugearite magmas from Mt. Etna display a wide range of behaviors as a function of temperature (T=1000–1200°C) and crystal content (ϕX=0.2–0.7). In the investigated T range, the deformation mechanism of these magmas varies from mainly brittle (T<1050°C) to mainly ductile (T>1085°C). At T=1075°C, both ductile and brittle behavior have been observed. In the ductile regime, these magmas behave as non-Newtonian fluids (at least up to T=1100°C) showing marked apparent shear thinning behavior. The observed rheological behavior is due to a complex response related to a non-homogenous deformation of the natural sample (e.g. viscous and/or brittle shear localization), favored by the presence of bubbles. Consequently, the obtained flow parameters can be considered as representative of the bulk rheology of natural magmas, commonly characterized by similar non-homogeneous deformation styles. We applied the obtained data to discuss the flow and emplacement conditions of these peculiar lava flows. We demonstrated that at eruptive temperatures, the presence of a pre-eruptive crystal cargo and bubbles facilitates the achievement of critical crystal content during flow and cooling, ultimately controlling lava transport and emplacement. Flow conditions can be maintained in the presence of an efficient insulation and, importantly, of deformable vesicles promoting and enhancing shear localization.

Field and experimental constraints on the rheology of arc basaltic lavas: the January 2014 Eruption of Pacaya (Guatemala)

We estimated the rheology of an active basaltic lava flow in the field, and compared it with experimental measurements carried out in the laboratory. In the field we mapped, sampled, and recorded videos of the 2014 flow on the southern flank of Pacaya, Guatemala. Velocimetry data extracted from videos allowed us to determine that lava traveled at ∼2.8 m/s on the steep ∼45° slope 50 m from the vent, while 550 m further downflow it was moving at only ∼0.3 m/s on a ∼4° slope. Estimates of effective viscosity based on Jeffreys’ equation increased from ∼7600 Pa s near the vent to ∼28,000 Pa s downflow. In the laboratory, we measured the viscosity of a representative lava composition using a concentric cylinder viscometer, at five different temperatures between 1234 and 1199 °C, with crystallinity increasing from 0.1 to 40 vol%. The rheological data were best fit by power law equations, with the flow index decreasing as crystal fraction increased, and no detectable yield strength. Although field-based estimates are based on lava characterized by a lower temperature, higher crystal and bubble fraction, and with a more complex petrographic texture, field estimates and laboratory measurements are mutually consistent and both indicate shear-thinning behavior. The complementary field and laboratory data sets allowed us to isolate the effects of different factors in determining the rheological evolution of the 2014 Pacaya flows. We assess the contributions of cooling, crystallization, and changing ground slope to the 3.7-fold increase in effective viscosity observed in the field over 550 m, and conclude that decreasing slope is the single most important factor over that distance. It follows that the complex relations between slope, flow velocity, and non-Newtonian lava rheology need to be incorporated into models of lava flow emplacement.