Outline
Comptes Rendus

Rheological changes in melts and magmas induced by crystallization and strain rate
Comptes Rendus. Géoscience, Volume 354 (2022) no. S1, pp. 227-248.

Abstract

This review highlights the rheological and phase proportions variation induced by cooling events from superliquidus temperature (melt) to subliquidus temperatures. It provides a comprehensive view of the rheological response of magmatic systems undergoing dynamic cooling and shear deformation. The two main parameters which are of importance to model the rheological properties of such crystallizing systems and which are simultaneously poorly investigated so far are crystallization and strain rates. The response to relatively high deformation rates results in shear thinning behavior in partly crystallized systems under variable shear rate and it should be considered in magmatic processes. Due to the sluggish crystallization of SiO 2 -rich melts, data are mainly available for mafic systems, which does not allow a general reappraisal. An attempt to model available literature data for less evolved systems in dynamic scenarios and a comparison with MELTS algorithm approach (thermodynamic equilibrium conditions) is provided. Since there are difficulties in comparing experimental data gained using different methodologies, we focus mainly on data obtained with the concentric cylinder technique. This highlights the fact that a general experimental protocol is needed in order to compare and model viscosity data to predict the dynamic rheological evolution for volcanic rocks.

Supplementary Materials:
Supplementary material for this article is supplied as a separate file:

Metadata
Received:
Accepted:
Online First:
Published online:
DOI: 10.5802/crgeos.125
Keywords: Viscosity, Crystal nucleation and growth, Magma rheology, Cooling rate, Modelling

Francesco Vetere 1, 2, 3; Gianluca Iezzi 2, 4; Diego Perugini 5; Francois Holtz 3

1 Dipartimento di Scienze Fisiche, della Terra e dell’Ambiente (DSTA) Università di Siena 53100, Siena, Italy
2 Dipartimento di Ingegneria & Geologia (InGeO), Università G. D’Annunzio di Chieti-Pescara, Via dei Vestini 31, 66100, Chieti, Italy
3 Institut of Mineralogy, Leibniz Universität Hannover, Callinstrasse 3, 30167 Hannover, Germany
4 Istituto Nazionale di Geofisica e Vulcanologia, INGV, Via di Vigna Murata 605, 00143, Roma, Italy
5 Department of Physics and Geology, University of Perugia, Piazza Università 1, 06100 Perugia, Italy
License: CC-BY 4.0
Copyrights: The authors retain unrestricted copyrights and publishing rights
@article{CRGEOS_2022__354_S1_227_0,
     author = {Francesco Vetere and Gianluca Iezzi and Diego Perugini and Francois Holtz},
     title = {Rheological changes in melts and magmas induced by crystallization and strain rate},
     journal = {Comptes Rendus. G\'eoscience},
     pages = {227--248},
     publisher = {Acad\'emie des sciences, Paris},
     volume = {354},
     number = {S1},
     year = {2022},
     doi = {10.5802/crgeos.125},
     language = {en},
}
TY  - JOUR
AU  - Francesco Vetere
AU  - Gianluca Iezzi
AU  - Diego Perugini
AU  - Francois Holtz
TI  - Rheological changes in melts and magmas induced by crystallization and strain rate
JO  - Comptes Rendus. Géoscience
PY  - 2022
SP  - 227
EP  - 248
VL  - 354
IS  - S1
PB  - Académie des sciences, Paris
DO  - 10.5802/crgeos.125
LA  - en
ID  - CRGEOS_2022__354_S1_227_0
ER  - 
%0 Journal Article
%A Francesco Vetere
%A Gianluca Iezzi
%A Diego Perugini
%A Francois Holtz
%T Rheological changes in melts and magmas induced by crystallization and strain rate
%J Comptes Rendus. Géoscience
%D 2022
%P 227-248
%V 354
%N S1
%I Académie des sciences, Paris
%R 10.5802/crgeos.125
%G en
%F CRGEOS_2022__354_S1_227_0
Francesco Vetere; Gianluca Iezzi; Diego Perugini; Francois Holtz. Rheological changes in melts and magmas induced by crystallization and strain rate. Comptes Rendus. Géoscience, Volume 354 (2022) no. S1, pp. 227-248. doi : 10.5802/crgeos.125. https://comptes-rendus.academie-sciences.fr/geoscience/articles/10.5802/crgeos.125/

Version originale du texte intégral (Propose a translation )

1. Introduction

Rheological behavior of silicate melts and magmas are of high importance to understand how they form, travel, and emplace into and on the Earth’s surface. A special role in the magmatic evolution is exerted by the melt/magma undercooling degree(s); it is defined as the difference between the system temperature and the temperature at which a phase saturates in the given liquid [Kirkpatrick 1981; Giuliani et al. 2021]. From a thermodynamic point of View, the appearance of solid phases from a magma is driven by the undercooling degree. However, the nucleation and crystals growth of a mineral phase is also governed by kinetics processes [Lofgren 1980; Dowty 1980; Hammer 2008; Iezzi et al. 2008, 2009; Vetere et al. 2013, 2015; Giuliani et al. 2020a, b]. The overall variation in nucleation rate spans many orders of magnitude, while growth rates variation is more limited [Swanson 1977; Kirkpatrick 1981; Lasaga 1997; Cashman 1988; Hammer and Rutherford 2002; Zhang 2008; Iezzi et al. 2009; Mollo and Hammer 2017; Giuliani et al. 2020a]. For a similar bulk magmatic composition, a low nucleation rate produces few crystals that will became large, while a high nucleation rate will result in a large number of tiny crystals [Hammer 2008; Iezzi et al. 2009; Higgins 2006; Mollo and Hammer 2017]. Nucleation is favored in silicate liquids with a high chemical diffusion, i.e., poor in SiO2 and rich in H2O, whereas SiO2- and H2O-poor liquids are reluctant to nucleate and prone to vitrify [Iezzi et al. 2009; Vetere et al. 2015].

Nucleation and crystallization processes have profound implications for the rheology of magmas and lavas, which depends on the volumetric amount of solid phases as well as on their size and shape [Gonnermann and Manga 2007]. Hence, a clear knowledge of nucleation and crystallization processes driven by solidification, either by cooling or by degassing-induced depressurization, is the key to understand magma rheology and its evolution before turning into a solid rock.

The eruptive style of volcanoes is strongly affected by magma viscoelastic properties [Dingwell 1996; Gonnermann and Manga 2007; Cassidy et al. 2018, and reference therein] and, as shown by Petrelli et al. [2018], their associated timescales. Indeed, the PT–fO2 path(s) for a starting magmatic composition (X) governs its solidification behavior and rheological pathway [Lanzafame et al. 2013, 2017].

Crystallization under kinetic conditions imposed by 𝛥T∕𝛥t and/or 𝛥P∕𝛥t was experimentally studied and investigated mainly for mafic compositions so far [Arzilli et al. 2019; Polacci et al. 2018; Hammer 2008; Conte et al. 2006; Fiege et al. 2015, and reference therein], while highly evolved melts are harsher to experimental investigation (time consuming), especially in a dry context [Vetere et al. 2013, 2015]. Nowadays, it is recognized that disequilibrium conditions often occur in natural processes [Kolzenburg et al. 2018; Vetere et al. 2019, 2021, and reference therein]. Moreover, ambiguities persist in our understanding of the interplay between crystal nucleation and growth on magma rheological variation (bulk viscosity and density) as it evolves from liquidus temperature down to the glass transition temperature, Tg. The cooling rate and degassing rate induced by depressurization of magmatic batches influence textural, petrological, and geochemical features, as well as volcanic behavior, i.e., from mildly lava effusion to violent explosive paroxysm [Armienti 2008; Hammer 2008; Applegarth et al. 2013; Shea and Hammer 2013; Fiege et al. 2015; Mollo and Hammer 2017; Arzilli et al. 2019, and references therein].

Due to the changes in residual melt chemistry and solid phase growth (as melt solidifies below the liquidus temperature, TL), it is expected that properties such as the apparent viscosity and density will vary, with significant impact on the magma’s ability to move [Giuliani et al. 2021]. To have an overview of this natural scenario and because of the paucity of direct field data [Chevrel et al. 2019, 2018, 2015], experimental approach is of paramount importance to further understand how magma behaves as it moves into and on the Earth.

In this review, after a general discussion on the nucleation and growth of minerals in silicate melts and magmas and on the general and most used laws and models governing magma rheology, we will focus both on recent experimental data on natural/synthetic systems and numerical approach with the aim to examine equilibrium versus disequilibrium conditions on the solidification of mafic compositions and its impact on rheological behavior. Specifically, we will constrain and model the rheological changes induced by cooling and also by shear stresses. Results point to the complex interactions between crystallization processes and rheological behavior as the key for understanding how magmas behave in natural systems. We will also emphasize the possible future experimental approaches in order to shed new light on these important natural processes that could be applicable to solidification induced by depressurization for volatile-rich magmas.

2. The state of the art

Petrological processes and scenarios are interpreted and reconstructed mainly by studying crystal-chemical attributes of phases [Giuliani et al. 2021, 2022]. This approach commonly assumes thermodynamic equilibrium between magmatic liquid and their mineral/volatile phases. From thermodynamically derived laws, intensive variables of magmatic systems can be obtained by crystal-chemistry of minerals and coexisting liquids, such as to track evolution of magmas [e.g., Putirka 2008, and references therein]. However, we now have clear evidence of disequilibrium features in both intrusive and especially volcanic rocks [e.g., De Paolo and Getty 1996; Roselle et al. 1997; Müller et al. 2004; Nabelek 2007; Baker 2008; Mollo and Hammer 2017; Giuliani et al. 2021]. Hence, the scientific community is aware that compositional changes of minerals in volcanic rocks are not necessarily related to equilibrium conditions, but rather and frequently to a disequilibrium state [e.g., Baker and Grove 1985; Hammer 2006, 2008; Iezzi et al. 2008, 2011; Del Gaudio et al. 2010; Mollo et al. 2010, 2011a, b]. For a bulk chemical composition of a magmatic batch, the increase of undercooling degree and deformation rate (i.e., stretching and folding mechanisms) progressively induces and facilitates nucleation and crystallization processes [Swanson and Ottino 1990; Perugini et al. 2008; Rossi et al. 2019; Tripoli et al. 2019]; in turn, they strongly affect rheological behavior of magmas that eventually rise from depth and emplace on the Earth’s surface. Deformations occurring in dynamic magmas enhance the nucleation and crystal growth; the movement of a magma furnishes continuous and fresh “feeding ingredients” around early formed crystals [Vetere et al. 2020; Vetere and Holtz 2020] when compared to a static condition. Also, in order to better understand the rheology of complex systems, another important aspect is related to the delay in crystal nucleation and growth that became progressively higher at low undercooling degree, i.e., close to liquidus temperature [Rusiecka et al. 2020; Giuliani et al. 2020b].

Rheology is defined as the science studying the deformation of materials under well-defined conditions. As defined in Mezner [2014], the term rheology (i.e., the flow science) derives from the Greek word “rhein” that is equivalent to the word “flow.” When applied to volcanology, the rheology of magmas essentially studies melt or magma flows as they are subjected to different stress, in magma chambers, along the way to the surface or when flowing on the Earth’s surface or other planetary surface (Mercury or Mars, for example). Starting from the definition of magma, it also implies that the rheological behavior depends mainly on phases present in a system (fluid ± crystals ± exsolved and dissolved volatiles), which depend on a variety of parameters such as pressure, temperature, and fO2 for a well-defined chemical system. In addition, it is now clear that such a behavior can change if phases are subjected to certain stress [e.g., shear rate dependence; Vetere and Holtz 2020 and reference therein]. In the conventional description of rheological data, in order to parameterize any suspension, there is a need to clearly define the relationship between the imposed shear stress (𝜎 in Pa) and the shear strain-rate ( 𝛾 ̇ in s−1; i.e., deformation rate). This relationship can be used to characterize the flow behavior and related flow curve. The flow regime is usually defined as Newtonian when there is a simple proportionality between stress and strain:

𝜎 = 𝜂 𝛾 ̇ , (1)
where 𝜂 (in Pa⋅s) is the Newtonian viscosity (Figure 1).

Figure 1.

Newtonian and non-Newtonian behavior [Lenk 1967]. 𝜏y is the yield stress that is to be overcome in order to initiate flow; n is the flow index describing the degree of non-Newtonian behavior. For Newtonian materials 𝜏y = 0, n = 1, and the flow consistency K =𝜂; shear-thickening and shear-thinning materials are defined by n > 1 and n < 1, respectively.

However, in multiphase and complex systems, such as magma, deviation from Newtonian behavior is observed and shear stress is not simply proportional to strain rate and the following cases are common: (1) a minimum stress is required before viscous flow initiates (also called the Bingham regime) or (2) both the above behaviors are combined, which is expressed by the Herschel–Bulkley model [Herschel and Bulkley 1926]:

𝜎 = 𝜏 y + K 𝛾 ̇ n , (2)
where 𝜏y is the yield stress, K is the flow consistency, and n the flow index (details are shown in Figure 1) [Herschel and Bulkley 1926; Wildemuth and Williams 1984; Mueller et al. 2010; Mader et al. 2013; Moitra and Gonnermann 2015].

The onset of non-Newtonian rheology in silicate melts has been investigated by Webb and Dingwell [1990a, b] using fiber elongation techniques. The chosen melt compositions were a high silica rhyolite, an andesite, a tholeiitic basalt, and a nephelinite. Viscosity was determined as the ratio of the applied tensile stress to the observed strain rate. As reported by the authors, “before experiments, each fiber was initially annealed at a temperature for a period of 104 s. This equilibration time is greater than both the shear and volume relaxation times of a melt with 𝜂v =𝜂 s = 1012 Pa⋅s and G (shear modulus) ≈ K (bulk modulus) 10 GPa (tV = tS = 100 s).” However, the measurements had no in situ monitoring of the sample temperature and did not consider viscous heating effects. Figure 2 illustrates the apparent viscosity reduction due to viscous heating effect. As reported by Cordonnier et al. [2012], “fluids with a temperature-dependent viscosity are affected by viscous heating once this heat source overcomes the heat loss (i.e., heat transfer by conduction, convection, and/or radiation).” This effect, despite increasing shear rates, will move magmas away from the glass transition [Gonnermann and Manga 2007]. Yue and Bruckner [1996] have re-analyzed the fiber elongation observations and formulated a correction for the strain-rate dependence of viscosity that includes viscous heating effects explicitly. The question as to whether viscous heating may have had a measurable effect on previous determinations of the onset of non-Newtonian viscosity remains unsolved.

Figure 2.

Results of viscosity measurements using the Parallel plate viscometry on NIST-certified starting glass as reported in Cordonnier et al. [2012]. (A) Load applied to the sample; (B) Due to the applied load (when it is high enough), the temperature in the sample increases; (C) As temperature increases (B) viscosity drops. This viscosity drop is explained from the temperature increase and can also be predicted by the Vogel–Tamann–Fulcher equation [Hess et al. 2007].

Because magmatic suspensions are strain rate dependent, thus sensitive to the applied shear stress, it is common to define the relative viscosity 𝜂r as:

𝜂 r = 𝜂 app 𝜂 l , (3)
where 𝜂app defines the ratio of a given stress to a given strain-rate, meaning that the obtained results are only valid for specific deformation conditions and, thus, vary with strain-rate, and 𝜂l is the liquid viscosity. Particular consideration must be given to the fact that multiphase suspension deforms when subjected to stress; for magmas the most common deformation relates to shear, while other types of deformation like oscillatory, tensile, or compressive behaviors are limited or even absent [details are given in Dingwell 1995; Dingwell and Webb 1989; Webb and Dingwell 1990a, b]. It is intuitive that the viscosity of a suspension is different from the viscosity of the liquid component. Moreover, the solid phases in suspensions can vary in terms of their volume, size, and shape, and this can drastically affect the rheology of the system [Kerr and Lister 1991; Mueller et al. 2010; Pistone et al. 2016; Arzilli et al. 2019; Kolzenburg et al. 2020, and reference therein]. As an example, prolate or oblate particles in magma could be present and can act as a different opponent to the applied shear stress and, as consequence, to the flow ability (e.g., if particles’ major axes are aligned to flow direction, their influence on flow velocity is minor compared with the opposed cases as with particles’ major axes orthogonal to the flow).

Vona et al. [2011] provided a methodology to follow for the evolution of the apparent viscosity as crystallization processes proceed. Figure 3a shows an S-shape curve where relaxing, incubation, nucleation times, and crystallization processes with a final plateau (denoting the end of the crystallization process) are highlighted. This cartoon is confirmed by experimental data for different compositions (as those presented in the following paragraph). The steepness of the curve is related to the response of the melt when temperature is below the liquidus. In other words, the curve shows the time that is needed for a melt to relax and eventually to start to nucleate a solid phase. As expected, very fluid systems react much faster than more viscous ones and this is clearly evidenced in Figures 3b–d where three different melts with compositions varying from andesite (typical intermediate composition between basalt and rhyolite) to basalt and to pyroxenite (an ultramafic igneous rock consisting essentially of minerals of the pyroxene group) are compared. Obviously, considering the red curves in Figure 3 (e.g., experiments performed at shearrate = 0.1 s−1) the most depolymerized pyroxenite composition (Figure 3d) reacts faster (ca. 1800 s) when compared to basalt (ca. 2 × 104 s; Figure 3c) and basalt reacts faster when compared to andesitic composition (ca. 3 × 104 s; Figure 3b).

Figure 3.

Note: Apparent viscosity variation versus time; (a) redraw from Vona et al. [2011] while (b–d) are redrawn from Vetere and Holtz [2020] and refer to (b) andesite from Calbuco 2015 eruptions, (c) Etna basalt from 122 A.C eruption, and (d) pyroxenite from Theo’s Flow which is a proxy for Martian lava flow samples. Blue and red curves refer to experiments performed at shear rate of 1 and 0.1 s−1, respectively. Note the shear-thinning behavior (apparent viscosity decreases as shear rate increases) with nucleation and crystallization processes related to the applied deformation rates (see text for details).

The last consideration on partly crystallized system concerns the approach to the so-called maximum packing fraction (𝛷m) where the viscosity of a suspension (thus, the ability of a melt plus crystal system to flow) tends to infinity and the transition from viscous to elastic behavior might occur. 𝛷m depends on particle sizes, shapes, and distributions and its value during a flow may change due to possible crystal alignment along the flow directions [see the concept of rheological critical melt percentage, RCMP as discussed in Arzi 1978 and reviewed by Rosenberg and Handy 2005].

3. Investigated compositions

The bulk compositions reconsidered and used in this review are reported in Figure 4 (Supplementary Appendix 1), on which both kinetics-controlled crystallization events, rheological test, and/or numerical modeling have been performed.

Figure 4.

Note: Total alkalis versus silica diagram for the investigated starting compositions. VH20 refers to Vetere and Holtz [2020]; V20, 19 and 17 refer to Vetere et al. [2020, 2019, 2017, respectively]; M20 refers to Morrison et al. [2020]; K20, 19 and 18 as in Kolzenburg et al. [2020, 2019, 2018, respectively]; SW20 and 15 refer to Sehlke and Whittington [2020, 2015]; C16 refers to Campagnola et al. [2016]; S14 refers to Sehlke et al. [2014]; Vo11 refers to Vona et al. [2011]; I09 to Ishibashi [2009] and IS07 to Ishibashi and Sato [2007].

A first compositional dataset comprises rheological experimental data in which shear stress vary at isothermal conditions, ideally at a fixed undercooling degree [Morrison et al. 2020; Sehlke and Whittington 2020, 2015; Vetere and Holtz 2020; Vetere et al. 2020, 2017; Campagnola et al. 2016; Sehlke et al. 2014; Vona et al. 2011; Ishibashi 2009; Ishibashi and Sato 2007]. However, it is important to state that magmatic suspensions are complex chemical systems, where crystallizing phases have different compositions from the parental liquids [Iezzi et al. 2020]. Thereby, also at isothermal condition, the crystallization trough time determines a change in composition of both crystalline and liquid phases. The second dataset, although very limited, comprises rheological data derived upon shearing and cooling events acting on silicate liquids, and represent the closest to natural scenarios [Vetere et al. 2019; Kolzenburg et al. 2020, 2019, 2018].

These two experimental datasets were obtained with the identical methodology (i.e., concentric cylinder) and mirrors a viscosity range of 101–106 Pa⋅s, thus, partly crystallized silicate melts, where the magma dynamic regime is dictated by crystallinity (𝛷). Usually this regime varies from dilute (where 𝛷 varies between 0.05 and 0.25 and particles have low to moderate interaction mainly depending on the particle aspect ratio) to semi dilute regime (∼0.25〈𝛷〉∼0.6) where particles strongly interact with each other with a possible alignment with the flow direction [Figure 5; Hoover et al. 2001; Saar et al. 2001; Vetere and Holtz 2020]. Both regimes define non-Newtonian systems with shear thinning behavior (small or strong, depending on 𝛷) with a possible evolution of the apparent yield stress (Figure 1). Those data reveal the importance of dynamic applied to silicate melts under cooling as will be discussed below.

Figure 5.

Concentric cylinder experimental results on Calbuco andesite at shear rate of 1 s−1 and temperature of 1483 K as in Vetere and Holtz [2020]. Phases are Fe–Ti oxides, plagioclase (plg), and glass. Note the plagioclase (plg) alignment along the flow direction.

4. Numerical simulation

AlphaMELTS code was used to track the phase evolution path for compositions presented in Figure 4 as temperature cools from the liquidus to subliquidus conditions. The alphaMELTS (V. 1.9) software provide a simple text-based interface to subroutine versions of the algorithms MELTS, pMELTS, and pHMELTS [Asimow and Ghiorso 1998; Smith and Asimow 2005; Ghiorso et al. 2002; Ghiorso and Sack 1995]. Suitable for multi-component systems, it allows to calculate equilibrium assemblages along a thermodynamic path set by the user (for details on the software, see the following address: http://melts.ofm-research.org/ and https://magmasource.caltech.edu/forum/).

The evolution of residual melt compositions is reported in Figure 6a and the evolution of melt proportion versus temperature is shown in Figure 6b (as vol%). In both cases, a direct comparison with literature data is provided. The calculated liquidus temperature are also available in Supplementary Appendix 1 for all the studied compositions. In general, it shows progressively lower values as the chemistry of the systems evolves to more silicic ones (higher SiO2 content), as expected, to equilibrium conditions for the modeled melts. The evolution of SiO2 predicted in Figure 6a for samples I09, Vo11Stromboli, and Vo11Etna, as for sample M20NYI and SW15NVP–Na, shows higher values than experimental data, while the opposite is observed for sample S14, M20NYA, and, in both cases, the differences become important at lower temperatures. In general, predictions of the Al2O3 evolutions are more accurate as those for CaO and MgO oxides, while FeOtot is overestimated for M20NYA, C16 and underestimated for Vo11Stromboli, SW15NVP, and M20NYI samples, respectively. The difference between experimental data and data predicted by alphaMELTS are particularly large for M20NYI and M20NYA compositions. As reported by Morrison et al. [2020], this is related to differences in oxygen fugacities (experimental data are obtained at more reduced conditions). Thus, we will focus only on rheological data for these compositions. The observed differences can be attributed to two aspects: (i) the inherited accuracy of the thermodynamic calculations and (ii) the occurrence of disequilibrium conditions experienced by the analyzed run-products (Figure 6a).

Figure 6a.

SiO2, Al2O3, FeOt, MgO, and CaO evolution of melts with temperature as systems cool down starting from liquidus temperatures. Symbols refer to experimental data while dotted lines refer to prediction obtained by using alphaMELTS 1.9. The label on each plot refers to chemical compositions as those reported in Figure 4.

Figure 6b.

Comparison of the volume of residual melt (vol%) derived from modeling and experimental data. The label in each plot corresponds to the sample numbers provided in Figures 4 and 6a. Red curves refer to experimental data while black ones are those derived from alphaMELTS 1.9 results. The red curves fitted through the experimental data sometimes show S-shapes because the curves were drawn to fit exactly the determined melt fractions. However, such a behavior is not realistic and reflects the large uncertainty in the determination of melt/crystals ratio. Data refer to ambient pressure.

Figure 6c.

(A) Plagioclase, (B) pyroxene, and (C) olivine compositions for some of the magmatic systems studied. Note the similar compositional range of experimental data and modeling.

Figure 6d.

Crystallinity versus temperature for the compositions reported in Figure 4.

Experimental and thermodynamic calculations show that melt fraction decreases with temperature, as expected. In general, the very low effect of temperature on crystallinity predicted by to alphaMELTS 1.9 in the high temperature range is not reproduced by the experimental data (in particular, for M20NYA, V17, C16, SW15NVP–Na, SW15_EB, S14, and I09 compositions). In most cases, the volume of residual melt is much lower that the predicted volumes and the liquidus temperatures extrapolated from the experimental data are lower than predicted liquidus temperatures. As a result, the effect of temperature on the crystallinity in the experimental products is more pronounced than predicted (see also Figure 6d with the crystallinity of the individual experiments).

The composition of solid phases was compiled from the literature data and are reported in Figure 6c(A–C). The compositions of plagioclase (plg), pyroxene (px), and olivine (ol) that crystallized in the experiments are compared to predicted compositions with alphaMELTS 1.9. Predicted and experimental compositions are in general good agreement. This indicates that alphaMELTS 1.9 is successful in predicting equilibrium mineral compositions but that the prediction of phase proportions is more problematic, which may be due to calibration of the model from static experiments. Applying a shear stress to crystallization experiments is expected to promote the crystallization and reduce the delay time preceding nucleation [Rusiecka et al. 2020; Vetere et al. 2021]. Thus, the possible mechanisms to explain this difference in crystallinity might be related to the continuous stirring of the melts at high shear rate in dynamic experiments which is bringing the necessary “feed growth ingredients(elements necessary to build a crystalline structure) close to individual crystal surfaces and which facilitates crystal growth (see below).

5. Modeling experimental data by varying shear rate and crystal content

Due to the intrinsic difficulties in performing experiments, the available experimental dataset on crystallization processes under dynamic conditions are still relatively scarce, as aforementioned. Moreover, the paucity of experimental data on more evolved systems is resulting from the limitations of the experimental technique (Figures 6a6d). In particular, the range of viscosity that can be investigated with the most common apparatus to reach high stress and higher viscosity values (e.g., 𝜂 > 106 Pa⋅s) does not cover that of silicic melts at moderate temperature and even with low amounts of crystals. We collected 415 literature experimental viscosity data (see Supplementary Appendices 1 and 2) on different mafic to intermediate compositions at different temperature, shear rate, volume of crystals, and 2D aspect ratio. All data were obtained with the identical methodology by using concentric cylinder apparatus and the majority show a decrease in the apparent viscosity (at isothermal conditions) as deformation rate increases, pointing to the well-known effect of shear thinning as discussed above.

Apart from viscosity data, one general outcome for most of the investigated systems is that increasing of the strain rate promotes a fast crystal growth rate. Experiments with dynamic crystallization show that only ca. 6000 s are necessary to produce ca. 22 area % of crystals having plagioclase crystal lengths of ca. 300 mm, mirroring a growth rate value of 5 × 10−6 cm/s [Vetere and Holtz 2020], much faster than those reported in literature [i.e., 10−8–10−9 cm/s; Pupier et al. 2008].

Considering that water and other volatiles affect the rheology of silicate melts [Fiege et al. 2015, and reference therein], the crystals’ growth rate might also be dependent on strain rates in dynamic hydrous magmas. Literature provides different general models able to reproduce apparent or relative viscosity behaviors as a function of crystal content, shear rate, and aspect ratio of the solid phases. Most of them use melt viscosity in order to retrieve the relative viscosity of the suspensions (melt + crystals; see (3)), adopting general viscosity models such as proposed by Hui and Zhang [2007] or Giordano et al. [2008]. A list of the most used models is provided in Table 1.

Figure 7 provides an overview of the relative viscosity data and some of the models provided in Table 1 are plotted for comparison. All the experimental data processed at different deformation rates show a shear thinning behavior (note identical symbols with same crystal content but different 𝜂r in Figure 7) and relative viscosity differences increase with increasing the gap between low and high shear rate. For example, Sehlke and Whittington [2015] show that NVP–Na viscosity vary from 7500 to 1500 Pa⋅s as shear rate vary from 0.01 to 5 s−1 at 1473 K with crystal content 29 vol%. Vetere et al. [2017] show a decrease of 𝜂app from ca. 10200 to 1100 Pa⋅s at 1502 K and ca. 28 vol% of crystal as shear rate vary from 0.1 to 5 s−1 in synthetic silicate melt compositions analogue for Mercury (for details, please refer to Supplementary Appendix 2).

Figure 7.

Relative viscosity evolution versus crystal content in vol%. (A, B) reproduce identical dataset with different scale in order to better appreciate difference in relative viscosity induced by different shear rate. Also, note in (B) a modified KD equation [Krieger and Dougherty 1959] able to reproduce relative viscosity for early departing systems such as those deformed at relatively low shear rate.

None of the models presented here are successful for predicting the whole experimental dataset. Reasons reside on the different parameters used for modeling. For example, the fact that 𝛷m can vary during deformation due to possible particle alignment will obviously change the rheological behavior. And this, up to now, is not predictable. Also, the use of image analysis to retrieve solid phase proportions can lead to some uncertainties and a statistically representative dataset of BSEs images selected along different orientations of the investigated sample is necessary. Thus, when using few BSEs images, AS (AS = crystals long axis/short axis), and the average volume proportions could be influenced by the orientation of the 2D section in a sample with a strong 3D crystal fabric [Lanzafame et al. 2017]. Ideally, analytical data based on tomographic approaches [Polacci et al. 2018; Arzilli et al. 2019] would give the solution to this problem but, up to now, especially when using concentric cylinder apparatus, it is highly difficult and techniques need improvements. Figure 7 shows also that it is possible to modify some of the parameters used in an equation in order to get a best fit with experimental data. This is the case, for instance, for the model proposed by Krieger and Dougherty [1959] which could better reproduce experimental data by varying the parameter B and 𝛷m as done in Figure 7B and Table 1. However, a more exhaustive rheological model has to engulf rheological transitions from a low viscosity, where flow regime is mainly those of the suspending liquids to a high viscosity regime, where parameters such as particle shape, orientation, and volume dominate.

Table 1.

Most common relative viscosity models

Model Literature reference
1 𝜂r = (1 + B𝜙) Einstein [1906]; Dilute system (𝛷< 0.25)
B = 2.5
2 𝜂r = 1 + B𝛷 + B1𝛷2 Guth and Gold [1938]; Dilute system (𝛷< 0.25) and spherical particles
B = 2.5, B1 = 14.1
3 𝜂 r = 1 + 1 . 2 5 𝛷 1 𝛷 𝛷 m 2 Eilers [1941, 1943]
4 𝜂 r = e B 𝛷 1 𝛼 𝛷 Vand [1948]
5 𝜂r = 1 −𝛼𝛷B Roscoe [1952];
𝛼 = 0.135, B = 2.5
6 𝜂 r = 1 𝛷 𝛷 m 2 Maron and Pierce [1956]
7 𝜂 r = 1 𝜙 𝜙 m B 𝜙 m Krieger and Dougherty [1959]
8 𝜂 r = 1 𝜙 𝜙 m B 𝜙 m Krieger and Dougherty [1959]; B = 2.5 for spherical particles; for elongated particles B could be modified
9 𝜂 r = exp 2 . 5 + 𝜙 𝜙 m 𝜙 𝛼 𝜙 𝜙 m Gay et al. [1969];
𝛼 = 0.48
10 𝜂 r = 1 𝛼 𝜙 𝜙 m 1 𝜙 𝜙 m 2 Chong et al. [1971];
𝛼 = 0.75
11 𝜂 r = 1 1 𝛷 𝛷 m 𝛷 m Wildemuth and Williams [1984]
12 𝜂 r = 1 + 3 π 8 𝛽 𝛽 + 1 3 + 4 . 5 𝛽 + 𝛽 2 𝛽 + 1 3 1 1 𝛽 ln ( 𝛽 + 1 ) Shapiro and Probstein [1992];
𝛽 = 𝜙 𝜙 m 1 3 1 𝜙 𝜙 m 1 3
13 𝜂 r = 1 + 𝜑 𝛿 1 𝛼 erf π 2 𝛼 𝜑 1 + 𝜑 𝛾 ̇ B 𝛷 m Costa [2005];
𝜑 = 𝜙 𝜙 m
14 𝜙 m = 0 . 0 6 6 4 9 9 tanh ( 0 . 9 1 3 4 2 4 ) log ( 𝛾 ̇ ) + 3 . 8 5 0 6 2 3 + 0 . 5 9 1 8 0 6
𝛿 = 6 . 3 0 1 0 9 5 tanh [ ( 0 . 8 1 8 4 9 6 ) log ( 𝛾 ̇ ) ( 0 . 8 1 8 4 9 6 ) log ( 𝛾 ̇ ) + 2 . 8 6 ] + 7 . 4 6 2 4 0 5
𝛼 = 0 . 0 0 0 3 7 8 tanh [ ( 1 . 1 4 8 1 0 1 ) log ( 𝛾 ̇ ) + 3 . 9 2 ] + 0 . 9 9 9 5 7 2
𝛾 = 3 . 9 8 7 8 1 5 tanh [ ( 0 . 8 9 0 8 ) log ( 𝛾 ̇ ) + 3 . 2 4 ] + 5 . 0 9 9 6 4 5
Caricchi et al. [2007] on Costa [2005] with variable strain rate
15 𝜂 r = 1 𝜑 𝛿 [ 1 F ( 𝜑 𝜀 𝛾 ) ] B 𝛷 m B 𝛷 m
where: F = ( 1 𝜉 ) erf π 2 ( 1 𝜉 ) 𝜑 ( 1 + 𝜑 𝛿 ) and
𝜂 = 𝜂 ( 𝛷 m ) = 2 1 ( 1 𝜉 ) erf π ( 1 𝜉 ) B 𝛷 m
Costa et al. [2009];
𝜑 = 𝜙 𝜙 m
16 log ( 𝜂 r ) = ( 0 . 0 6 exp ( 6 𝛷 ) ) log 𝜏 𝜏 crit Cordonnier et al. [2009]; 𝜎 is the differential stress, 𝜎crit is the limiting stress (109 Pa), and (0.06 exp(6∗𝜙)) is the crystal stress dependency relation with 𝜙 the crystal fraction.
17 𝜂 r = 1 𝛷 𝛷 m 𝛼 𝛷 m 1 + 𝛽 ln ( 1 𝛷 𝛷 m ) ln 𝛾 ̇ Ishibashi [2009]
18 𝜂 r = 1 𝛷 𝛷 m 2 [ 1 𝛼 log ( 𝛾 ̇ ) ] Vona et al. [2011]
19 𝜂 r = 1 𝛷 𝛷 m 𝛼 Moitra and Gonnermann [2015];
𝛼 = 1.92
20 𝜂 r = 1 𝛷 m 𝛷 𝛷 m ( 1 𝛷 ) 2 . 5 𝛷 m 1 𝛷 m Faroughi and Huber [2015]
21 𝜂 r = 𝜂 r ( max ) 1 exp ( k 𝛷 n ) 1 exp ( k ) Liu et al. [2017];
𝜂r(max) is value for 𝛷 = 1. k and n are adjustable parameters

Please refer to the provided literature for details.

We performed diverse test on the dataset provided in Supplementary Appendix 2 for constraining a possible comprehensive apparent/relative viscosity model. The best solution was given by a polynomial approach, although it doesn’t provide excellent results. We propose an empirical approach to retrieve the apparent viscosity from (1) shear rate and volume of solid phases or (2) AS and volume of crystals in silicate melt suspensions. Results are based on the 415 experimental data presented in Supplementary Appendix 2, interpolated by a simple polynomial approach able to build a 3D surface on which log viscosity, shear rate, and volume of crystals in one case and log viscosity, AS, and volume of crystals in the other case are taken into account (Figure 8). The apparent viscosity can be estimated using the following expressions:

log 𝜂 app = a + b ln 𝛾 ̇ + c 𝛷 +d ( ln 𝛾 ̇ ) 2 + e 𝛷 2 + f 𝛷 ln 𝛾 ̇ (4)
and
log 𝜂 app = a + b lnAS + c AS 1 . 5 𝛷 +d exp ( AS ) + e 𝛷 ln 𝛷 , (5)
where 𝜂app is the apparent viscosity, 𝛾 the shear rate, 𝛷 the crystal content (vol%), and AS the crystal’s aspect ratio. The fixed parameters for (4) are a = 1.917, b = −0.034, c = 9.258, d = −0.032, e = −14.958, and f = −0.737. The parameters for (5) are a = 1.419, b = 0.313, c = −5.952, d = 10.499, and e = −6.481. These parameters were obtained by non-linear regressions giving Root Mean Square Error (RMSE) of 0.316 and 0.319 and Mean Absolute Error (MAE) of 0.257 and 0.255 for (4) and (5), respectively. The comparison of the predicted 𝜂app from (4) and (5) with experimental data is shown in Figure 9A and B, respectively.

Figure 8.

Surface modeling of apparent viscosity with respect of shear rate (𝛾 in s−1) and crystal content (𝛷; A) and aspect ratio (AS) and crystal content (𝛷; B). The model is based on 415 experimental data provided in Supplementary Appendix 2 derived from chem composition as in Figure 4 and Supplementary Appendix 1. Polynomial approaches (giving the best results in terms of R2) are provided in Table 1 together with best fit parameters and R2.

Figure 9.

log𝜂 observed versus predicted from (4) and (5). (A) refers to (4) with R2 = 0.80 while (B) refers to (5) with R2 = 0.75.

6. The effect of shear rate, crystal content, and cooling rates

The concomitant effect of shear rate and cooling on crystallinity and rheology can be retrieved from experimental data at one atmosphere pressure, mainly obtained for mafic composition. However, the rheology of systems crystallizing as a result of degassing-induced decompression as a function of decompression-rate and shear rate was never investigated experimentally. Such experiments would simulate the nearest approach to geological relevant conditions, specifically the ascent of magmas in conduits and dikes.

Vetere et al. [2020, 2019] and Kolzenburg et al. [2020, 2019, 2018, 2017] have performed experiments at one atmosphere pressure. Due to the progressive increase in crystallinity with time during cooling (contrasting to the classical isothermal viscosity experiments presented in paragraphs 4 and 5), the apparent viscosity will never reach a steady value because of the non-equilibrium conditions of the liquid plus crystal suspensions. As shown in Kolzenburg et al. [2018], crystal growth is progressively limited and potentially impeded in static systems (e.g., no shear) because the surrounding (residual intra-crystalline) melt is depleted in the compatible elements (already incorporated by the early growth events) required for further crystal growth. Also, the cooling and the first crystallizing phases enrich the residual and intra-crystalline melts in slow diffusion components (SiO2 and Al2O3), while the drop of temperature determines a progressive slowdown of chemical elements. Dendritic growth is the consequence of such situations, when the movement of chemical species is limited in a solidifying silicate liquid; if the diffusion is so sluggish compared to kinetic condition, nucleation is suppressed and a glass forms [Kirkpatrick 1981; Lofgren 1980; Shea and Hammer 2013; Vetere et al. 2013, 2015]. Shear or any other deformation acting on magmas and silicate liquids enhance the diffusion of elements by mass transfer.

As reported in Kolzenburg et al. [2017], systems cooled in the same thermal range crystallize faster (at low 𝛥T and high-T) and more abundantly with progressively increasing shear rates. This observation was also reported for isothermal experiments at sub-liquidus conditions by Kouchi et al. [1986], Vona and Romano [2013], Chevrel et al. [2015], and Vetere et al. [2017] and it is now widely accepted that increasing of shear rate increases the crystal growth rate. Also, Vetere et al. [2019, 2020] have shown that the shear rate during cooling influences the growth rate of solid phases. In particular, the growth rate of olivine, in pyroxenite melt at cooling rate of 10 K/h is ca. 5.3 × 10−6 cm/s. Clinopyroxene crystal in shoshonitic melts shows growth rates between 1.2 and 2.1 × 10−6 cm/s as cooling rates vary from 10 to 100 K/h, respectively, under a deformation rate of 1 s−1.

As shown in a simulation for magmatic conditions [Petrelli et al. 2016], mechanisms invoked as responsible for this fast crystal growth rate could be a continuous “fresh feeding ingredient” on individual crystal surfaces. The fast stirring also enhances the formation of new crystal nuclei, as suggested by Emerson [1926], Cashman et al. [1999], and Vona and Romano [2013]. Finally, oxygen fugacity also plays an important role on crystallization processes and related rheological changes [Kolzenburg et al. 2020]. For example, increasing of fO2 in magmas reaching the Earth’s surface can lead to a progressive shift of the equilibrium crystallization temperatures.

An example for the evolution of apparent viscosity with decreasing temperature and concomitant crystallization is given in Figure 10 for three different mafic magmas: pyroxenite, basalt, and shoshonite. Experimental conditions in Figure 10 are identical for pyroxenite and shoshonite (cooling rate of 10 K/h and shear rate of 1 s−1), while for basalt the applied cooling rate is slightly higher (30 K/h) and shear rate slightly lower (0.77 s−1). These differences are somewhat minimal, but in terms of reactivity to the cooling event, pyroxenite is the fastest. The kinks in the evolution of viscosity with temperature decrease are related to the onset of crystallization of mineral phases. Pyroxenite, basalt, and shoshonite compositions (composition reported in Figure 4 and Supplementary Appendix 1) have a liquidus temperature, computed by using alphaMELTS of 1543, 1444, and 1452 K, respectively. Obviously, crystallization occurs very quickly after reaching the liquidus temperature for the pyroxenite. The temperature (and time) interval between liquidus and crystallization temperature is larger for the basalt, and even larger for shoshonite. In structural terms, the chemical variability of silicate liquids might be described by the number of non-bridging oxygens (NBO) per tetrahedrally coordinated cation (T) [Mysen and Richet 2005]:

NBO T = 1 T i = 1 i n M i n + , (6)
where T is the total atomic abundance of tetrahedrally coordinated cations, M i n + is the proportion of network modifying cations, “i”, with electrical charge n+ after subtraction of the portion required for charge-balancing trivalent cations on tetravalent sites [Mysen and Richet 2005]. The resulting NBO/T values are 1.14, 0.66, and 0.36 for pyroxenite, basalt, and shoshonite melts, respectively. This shows that high NBO/T values are related with higher cation mobility favoring fast rearrangement of nuclei in order to achieve the “minimum structural unit” for a further crystal growth [Mysen and Richet 2005]. The curves in Figure 10 reflect the effect of melt structure on the nucleation and crystal growth delay during cooling, as also reported in Vetere et al. [2015] and Iezzi et al. [2017]. Nevertheless, more data are needed to possibly define a better relation between cooling rate, applied shear stress, crystal growth, and NBO/T in order to model viscosity evolution of crystallizing magmas.

Finally, it has recently become evident that the formation of Fe-oxide nanolites can affect the composition of the residual melt and induce change in crystal nucleation and growth processes and thus, on the magma rheological response [Vetere et al. 2015; Di Genova et al. 2020, and reference therein]. While these phases are too small to be analytically quantified with most of the classical techniques applied to characterize the presented dataset, the potential occurrence of these phases might influence the data in some experimental charges.

Figure 10.

Log apparent viscosity (Pa⋅s) versus temperature (K) evolution for basaltic, pyroxenite, and shoshonite melt compositions as in Figure 4 and Supplementary Appendix 1. Note the decrease of viscosity as melts cool down and crystallize.

7. Conclusion and future challenges

This review shows that silicate melts subjected to different deformation rates show shear thinning behavior both on pure melts as well as on partly crystallized systems, where crystal content does not exceed ca. 40 vol%. It is obvious that there is still a need for experiments aimed at increasing the available dataset simulating viscosity evolution in crystallizing dynamic systems. Future studies need also to address the effect of pressure, degassing, and related crystallization to expand disequilibrium rheology models to crustal magmatic conditions. A common and rationale strategy is required to study magma rheology by using a common and systematic experimental protocol and this is particularly true for the concentric cylinder approach. From a general point of view, deformation applied to silicate melts enhances nucleation and growth rates of phases when compared to static conditions. This point needs to be considered in modeling magmatic and volcanic behaviors, since rheological behavior could drastically deviate from that predicted based on static experimental approaches and equilibrium modeling.

Conflicts of interest

Authors have no conflict of interest to declare.

Acknowledgements

This study was funded by the “Fondi Ateneo of the University G. D’Annunzio,” PRIN (2009PZ47NA_003) project “Experimental determination of the glass-forming ability (GFA), nucleation and crystallization of natural silicate melts,” and PRIN (2017J277S9_003) project “Time scales of solidification in magmas: Application to Volcanic Eruptions, Silicate Melts, Glasses, Glass-Ceramics” awarded to GI. Alexander von Humbold foundation senior research grant to FV is acknowledged.


References

[Applegarth et al., 2013] L. J. Applegarth; H. Tuffen; M. R. James; H. Pinkerton Degassing-driven crystallisation in basalts, Earth Sci. Rev., Volume 116 (2013), pp. 1-16 | DOI

[Armienti, 2008] P. Armienti Decryption of igneous rock textures: crystal size distribution tools, Rev. Mineral., Volume 69 (2008), pp. 623-649 | DOI

[Arzi, 1978] A. Arzi Critical phenomena in the rheology of partially melted rocks, Tectonophysics, Volume 44 (1978), pp. 173-184 | DOI

[Arzilli et al., 2019] F. Arzilli; G. La Spina; M. R. Burton; M. Polacci; N. Le Gall; M. E. Hartley et al. Magma fragmentation in highly explosive basaltic eruptions induced by rapid crystallization, Nat. Geosci., Volume 12 (2019), pp. 1023-1028 | DOI

[Asimow and Ghiorso, 1998] P. D. Asimow; M. S. Ghiorso Algorithmic modifications extending MELTS to calculate subsolidus phase relations, Am. Mineral., Volume 83 (1998), pp. 1127-1132 | DOI

[Baker and Grove, 1985] M. B. Baker; T. L. Grove Kinetics controls on pyroxene nucleation and metastable liquid line of descent in a basaltic andesite, Am. Mineral., Volume 70 (1985), pp. 279-287

[Baker, 2008] D. R. Baker The fidelity of melt inclusions as records of melt composition, Contrib. Miner. Petrol., Volume 156 (2008), pp. 377-395 | DOI

[Campagnola et al., 2016] S. Campagnola; C. Romano; A. Vona; C. G. Giordano Crystallization kinetics and rheology of leucite-bearing tephriphonolite magmas from the Colli Albani volcano (Italy), Chem. Geol., Volume 424 (2016), pp. 12-29 | DOI

[Caricchi et al., 2007] L. Caricchi; L. Burlini; P. Ulmer; T. Gerya; M. Vassalli; P. Papale Non-Newtonian rheology of crystal-bearing magmas and implications for magma ascent dynamics, Earth Planet Sci. Lett., Volume 264 (2007), pp. 402-419 | DOI

[Cashman et al., 1999] K. Cashman; C. Thornber; J. Kauahikaua Cooling and crystallization of lava in open channels, and the transition of Pahoehoe Lava to “A”(a)over-bar, Bull. Volcanol., Volume 61 (1999), pp. 306-323 | DOI

[Cashman, 1988] K. Cashman Crystallization of Mount-St-Helens 1980–1986 dacite—a quantitative textural approach, Bull. Volcanol., Volume 50 (1988), pp. 194-209 | DOI

[Cassidy et al., 2018] M. Cassidy; M. Manga; K. Cashman; O. Bachmann Controls on explosive-effusive volcanic eruption styles, Nat. Commun., Volume 9 (2018), 2839 | DOI

[Chevrel et al., 2015] M. O. Chevrel; C. Cimarelli; L. deBiasi; J. B. Hanson; Y. Lavallée; F. Arzilli; D. B. Dingwell Viscosity measurements of crystallizing andesite from Tungurahua volcano (Ecuador), Geochem. Geophys. Geosyst., Volume 16 (2015), pp. 870-889 | DOI

[Chevrel et al., 2018] M. O. Chevrel; A. J. Harris; M. R. James; L. Calabrò; L. Gurioli; H. Pinkerton The viscosity of pāhoehoe lava: in situ syn-eruptive measurements from Kilauea, Hawaii 2, Earth Planet. Sci. Lett., Volume 493 (2018), pp. 161-171 | DOI

[Chevrel et al., 2019] M. O. Chevrel; H. Pinkerton; A. J. L. Harris Measuring the viscosity of lava in the field: A review, Earth-Sci. Rev., Volume 196 (2019), 102852 | DOI

[Chong et al., 1971] J. Chong; E. Christiansen; A. Baer Rheology of concentrated suspensions, J. Appl. Polym. Sci., Volume 15 (1971), pp. 2007-2021 | DOI

[Conte et al., 2006] A. M. Conte; C. Perinelli; R. Trigila Cooling kinetics experiments on different Stromboli lavas: Effects on crystal morphologies and phases composition, J. Volcanol. Geotherm. Res., Volume 155 (2006) no. 3–4, pp. 179-200 | DOI

[Cordonnier et al., 2009] B. Cordonnier; K. U. Hess; Y. Lavallee; D. B. Dingwell Rheological properties of dome lavas: case study of Unzen volcano, Earth Planet Sci. Lett., Volume 279 (2009), pp. 263-272 | DOI

[Cordonnier et al., 2012] B. Cordonnier; S. M. Schmalholz; K. U. Hess; D. B. Dingwell Viscous heating in silicate melts: An experimental and numerical comparison, J. Geophys. Res., Volume 117 (2012), B02203 | DOI

[Costa et al., 2009] A. Costa; L. Caricchi; N. Bagdassarov A model for the rheology of particle-bearing suspensions and partially molten rocks, Geochem. Geophys. Geosyst., Volume 10 (2009) no. 3, Q03010 | DOI

[Costa, 2005] A. Costa Viscosity of high crystal content melts: dependence on solid fraction, Geophys. Res. Lett., Volume 32 (2005), L22308 | DOI

[Del Gaudio et al., 2010] P. Del Gaudio; S. Mollo; G. Ventura; G. Iezzi; J. Taddeucci; A. Cavallo Cooling rate-induced differentiation in anhydrous and hydrous basalts at 500 MPa: implications for the storage and transport of magmas in dikes, Chem. Geol., Volume 270 (2010), pp. 164-178 | DOI

[De Paolo and Getty, 1996] D. J. De Paolo; S. R. Getty Models of isotopic exchange in reactive fluid-rock systems: implications for geochronology in metamorphic rock, Geochim. Cosmochim. Acta, Volume 60 (1996), pp. 3933-3947 | DOI

[Dingwell and Webb, 1989] D. B. Dingwell; S. L. Webb Structural relaxation in silicate melts and non-Newtonian melts rheology in geologic processes, Phys. Chem. Miner., Volume 16 (1989), pp. 508-516 | DOI

[Dingwell, 1995] D. B. Dingwell Viscosity and anelasticity of melts, Mineral Physics and Crystallography (A Handbook of Physical Constants AGU Reference Shelf), Volume 2, American Geophysical Union, Washington, 1995, pp. 209-217 | DOI

[Dingwell, 1996] D. B. Dingwell Volcanic dilemma: flow or blow?, Science, Volume 273 (1996), pp. 1054-1055 | DOI

[Di Genova et al., 2020] D. Di Genova et al. In situ observation of nanolite growth in volcanic melt: A driving force for explosive eruptions, Sci. Adv., Volume 6 (2020), eabb0413 | DOI

[Dowty, 1980] E. Dowty Crystal growth and nucleation theory and the numerical simulation of igneous crystallization, The Physics of Magmatic Processes (R. B. Hargraves, ed.), Princeton University Press, Princeton, NJ, USA, 1980, pp. 419-485 | DOI

[Eilers, 1941] V. H. Eilers Die viskosität von emulsionen hochviskoser stoffe als funktion der konzentration, Kolloid-Zeitschrift, Volume 97 (1941), pp. 313-321 | DOI

[Eilers, 1943] H. Eilers Die viskositäts-konzentrationsabhängigkeit kolloider systeme in organischen lösungsmitteln, Kolloid-Zeitschrift, Volume 102 (1943), pp. 154-169 | DOI

[Einstein, 1906] A. Einstein Eine neue Bestimmung der Molekul dimensionen, Ann. Phys., Volume 19 (1906), pp. 289-306 | DOI | Zbl

[Emerson, 1926] O. H. Emerson The formation of aa and pahoehoe, Am. J. Sci., Volume 68 (1926), pp. 109-114 | DOI

[Faroughi and Huber, 2015] S. A. Faroughi; C. Huber A generalized equation for rheology of emulsions and suspensions of deformable particles subjected to simple shear at low Reynolds number, Rheol. Acta, Volume 54 (2015), pp. 85-108 | DOI

[Fiege et al., 2015] A. Fiege; F. Vetere; G. Iezzi; A. Simon; F. Holtz The roles of decompression rate and volatiles (H 2 O + Cl ± CO 2 ± S) on crystallization in (trachy-) basaltic magma, Chem. Geol., Volume 411 (2015), pp. 310-322 | DOI

[Gay et al., 1969] E. Gay; P. Nelson; W. Armstrong Flow properties of suspensions with high solids concentration, AIChE J., Volume 15 (1969), pp. 815-822 | DOI

[Ghiorso and Sack, 1995] M. S. Ghiorso; R. O. Sack Chemical mass transfer in magmatic processes: IV. A revised and internally consistent thermodynamic model for the interpretation and extrapolation of liquid–solid equilibria in magmatic systems at elevated temperatures and pressures, Contrib. Mineral. Petrol., Volume 119 (1995), pp. 197-212 | DOI

[Ghiorso et al., 2002] M. S. Ghiorso; M. M. Hirschmann; P. W. Reiners; V. C. Kress III The pMELTS: A revision of MELTS for improved calculation of phase relations and major element partitioning related to partial melting of the mantle to 3 GPa, Geochem. Geophys. Geosyst., Volume 3 (2002) no. 5, pp. 1-36 | DOI

[Giordano et al., 2008] D. Giordano; J. K. Russell; D. B. Dingwell Viscosity of magmatic liquids: a model, Earth Planet. Sci. Lett., Volume 271 (2008), pp. 123-134 | DOI

[Giuliani et al., 2020a] L. Giuliani; G. Iezzi; H. Tyler; M. Davis; A. Elbrecht; F. Vetere; M. Nazzari; S. Mollo The onset and solidification path of a MORB by in-situ (DSC) and ex-situ investigations, Front. Earth Sci., Volume 8 (2020), 337 | DOI

[Giuliani et al., 2020b] L. Giuliani; G. Iezzi; F. Vetere; H. Behrens; S. Mollo; F. Cauti; G. Ventura; P. Scarlato Evolution of textures, crystal size distributions and growth rates of plagioclase, clinopyroxene and spinel solidified at variable cooling rates from a mid-ocean ridge basaltic liquid, Earth Sci. Rev., Volume 204 (2020), 103165 | DOI

[Giuliani et al., 2021] L. Giuliani; G. Iezzi; S. Mollo Dynamics of volcanic systems: physical and chemical models applied to equilibrium versus disequilibrium solidification of magmas, Dynamic Magma Evolution (F. Vetere, ed.) (American Geophysical Union, Geophysical Monograph), Volume 254, John Wiley & Sons, Inc., Hoboken, NJ, USA, 2021, pp. 101-132 | DOI

[Giuliani et al., 2022] L. Giuliani; F. Vetere; G. Iezzi; M. Nazzari; S. Mollo; H. Behrens; P. Scarlato; G. Ventura Crystal-chemical variations of spinel, clinopyroxene, and plagioclase in MORB basaltic melt induced by continuous cooling, Chem. Geol., Volume 594 (2022), 120765 | DOI

[Gonnermann and Manga, 2007] H. M. Gonnermann; M. Manga The fluid mechanics inside a volcano, Annu. Rev. Fluid Mech., Volume 39 (2007), pp. 321-356 | DOI | MR | Zbl

[Guth and Gold, 1938] E. Guth; O. Gold Viscosity and electroviscous effect of the AgI sol. II. Influence of the concentration of AgI and of electrolyte on the viscosity, Phys. Rev., Volume 53 (1938), p. 322

[Hammer and Rutherford, 2002] J. E. Hammer; M. J. Rutherford An experimental study of the kinetics of decompression-induced crystallization in silicic melt, J. Geophys. Res., Volume 107 (2002) no. B1, 2021 | DOI

[Hammer, 2006] J. E. Hammer Influence of fO 2 and cooling rate on the kinetics and energetics of Fe-rich basalt crystallization, Earth Planet. Sci. Lett., Volume 248 (2006), pp. 618-637 | DOI

[Hammer, 2008] J. E. Hammer Experimental studies of the kinetics and energetics of magma crystallization, Rev. Mineral. Geochem., Volume 69 (2008) no. 1, pp. 9-59 | DOI

[Herschel and Bulkley, 1926] W. Herschel; R. Bulkley Konsistenzmessungen von Gummi-Benzollösungen, Kolloid-Z., Volume 39 (1926), pp. 291-300 (Measurement of consistency as applied to rubber-benzene solutions, Proc. ASTM 82, 1925) | DOI

[Hess et al., 2007] K. U. Hess; B. Cordonnier; Y. Lavallee; D. B. Dingwell High-load, high-temperature deformation apparatus for synthetic and natural silicate melts, Rev. Sci. Instrum., Volume 78 (2007) no. 7, 075102 | DOI

[Higgins, 2006] M. D. Higgins Quantitative Textural Measurements in Igneous and Metamorphic Petrology, Cambridge University Press, Cambridge, 2006 | DOI

[Hoover et al., 2001] S. R. Hoover; K. V. Cashman; M. Manga The yield strength of subliquidus basalts—experimental results, J. Volcanol. Geotherm. Res., Volume 107 (2001), pp. 1-18 | DOI

[Hui and Zhang, 2007] H. Hui; Y. Zhang Toward a general viscosity equation for natural anhydrous and hydrous silicate melts., Geochim. Cosmochim. Acta, Volume 71 (2007) no. 2, pp. 403-416 | DOI

[Iezzi et al., 2008] G. Iezzi; S. Mollo; G. Ventura; A. Cavallo; C. Romano Experimental solidification of anhydrous latitic and trachytic melts at different cooling rates: The role of nucleation kinetics, Chem. Geol., Volume 253 (2008), pp. 91-101 | DOI

[Iezzi et al., 2009] G. Iezzi; S. Mollo; G. Ventura Solidification behaviour of natural silicate melts and volcanological implications, Volcanoes: Formation, Eruptions and Modelling (N. Lewis; A. Moretti, eds.), Nova Publishers, New York, 2009, pp. 127-151

[Iezzi et al., 2011] G. Iezzi; S. Mollo; G. Torresi; G. Ventura; A. Cavallo; P. Scarlato Experimental solidification of an andesitic melt by cooling, Chem. Geol., Volume 283 (2011), pp. 261-273 | DOI

[Iezzi et al., 2017] G. Iezzi; A. L. Elbrecht; M. Davis; F. Vetere; V. Misiti; S. Mollo; A. Cavallo Glass Stability (GS) of chemically complex (natural) sub-alkaline glasses, J. Non Cryst. Solids, Volume 477 (2017), pp. 21-30 | DOI

[Iezzi et al., 2020] G. Iezzi; J. E. Hammer; A. Whittington; D. R. Neuville Editorial: research topic crystal nucleation and growth in magmatic suspensions, Front. Earth Sci., Volume 8 (2020), 607972 | DOI

[Ishibashi and Sato, 2007] H. Ishibashi; H. Sato Viscosity measurements of subliquidus magmas: Alkali olivine basalt from the Higashi-Matsuura district, Southwest Japan, J. Volcanol. Geotherm. Res., Volume 160 (2007) no. 2007, pp. 223-238 | DOI

[Ishibashi, 2009] H. Ishibashi Non-Newtonian behavior of plagioclase-bearing basaltic magma: Subliquidus viscosity measurement of the 1707 basalt of Fuji volcano, Jpn J. Volcanol. Geotherm. Res., Volume 181 (2009) no. 1–2, pp. 78-88 | DOI

[Kerr and Lister, 1991] R. C. Kerr; J. R. Lister The effects of shape on crystal settling and on the rheology of magmas, J. Geol., Volume 99 (1991), pp. 457-467 | DOI

[Kirkpatrick, 1981] R. J. Kirkpatrick Kinetics in crystallization of igneous rocks, Rev. Mineral. Geochem., Volume 8 (1981), pp. 321-397

[Kolzenburg et al., 2017] S. Kolzenburg; D. Giordano; T. Thordarson; A. Höskuldsson; D. B. Dingwell The rheological evolution of the 2014/2015 eruption at Holuhraun, central Iceland, Bull. Volcanol., Volume 79 (2017) no. 6, 45 | DOI

[Kolzenburg et al., 2018] S. Kolzenburg; D. Giordano; K. U. Hess; D. B. Dingwell Shear rate-dependent disequilibrium rheology and dynamics of basalt solidification, Geophys. Res. Lett., Volume 45 (2018), pp. 6466-6475 | DOI

[Kolzenburg et al., 2019] S. Kolzenburg; D. Giordano; A. Di Muro; D. B. Dingwell Equilibrium viscosity and disequilibrium rheology of high magnesium basalt from Piton de la Fournaise volcano, La Reunion, Indian Ocean, France, Ann. Geophys., Volume 62 (2019) no. 2, VO218 | DOI

[Kolzenburg et al., 2020] S. Kolzenburg; K. U. Hess; K. Berlo; D. B. Dingwell Disequilibrium rheology and crystallization kinetics of basalts and implications for the phlegrean volcanic district, Front. Earth Sci., Volume 8 (2020), 187 | DOI

[Kouchi et al., 1986] A. Kouchi; A. Tsuchiyama; I. Sunagawa Effect of stirring on crystallization kinetics of basalt: Texture and element partitioning, Contrib, Mineral. Petrol., Volume 9 (1986), pp. 429-438 | DOI

[Krieger and Dougherty, 1959] I. M. Krieger; T. J. Dougherty A mechanism for non-Newtonian flow in suspensions of rigid spheres, Trans. Soc. Rheol., Volume 3 (1959), pp. 137-152 | DOI

[Lanzafame et al., 2013] G. Lanzafame; S. Mollo; G. Iezzi; C. Ferlito; G. Ventura Unravelling the solidification path of a pahoehoe “cicirara” lava from Mount Etna volcano, Bull. Volcanol., Volume 75 (2013) no. 4, pp. 1-16 | DOI

[Lanzafame et al., 2017] G. Lanzafame; G. Iezzi; L. Mancini; F. Lezzi; S. Mollo; C. Ferlito Solidification and Turbulence (Non-laminar) during magma ascent: Insights from 2D and 3D analyses of bubbles and minerals in an Etnean dyke, J. Petrol., Volume 58 (2017) no. 8, pp. 1511-1534 | DOI

[Lasaga, 1997] A. C. Lasaga Kinetic Theory in the Earth Sciences, Princeton University Press, Princeton, NJ, USA, 1997

[Lenk, 1967] R. Lenk A generalized flow theory, J. Appl. Polym. Sci., Volume 11 (1967), pp. 1033-1042 | DOI

[Liu et al., 2017] Z. Liu; L. Pandelaers; B. Blanpain; M. Guo Viscosity of heterogeneous silicate melts: assessment of the measured data and modeling, ISIJ Int., Volume 57 (2017), pp. 1895-1901 | DOI

[Lofgren, 1980] G. Lofgren Experimental studies on the dynamic crystallisation of silicate melts, The Physics of Magmatic Processes (R. B. Hargraves, ed.), Princeton University Press, Princeton, NJ, USA, 1980, pp. 487-551 | DOI

[Mader et al., 2013] H. M. Mader; E. W. Llewellin; S. P. Mueller The rheology of two-phase magmas: A review and analysis, J. Vol. Geotherm. Res., Volume 257 (2013), pp. 135-158 | DOI

[Maron and Pierce, 1956] S. H. Maron; P. E. Pierce Application of Ree–Eyring generalized flow theory to suspensions of spherical particles, J. Colloid Sci., Volume 11 (1956), pp. 80-95 | DOI

[Mezner, 2014] T. G. Mezner The Rheology Handbook, Vincentz Network, Hanover, 2014

[Moitra and Gonnermann, 2015] P. Moitra; H. M. Gonnermann Effects of crystal shape- and size-modality on magma rheology, Geochem. Geophys. Geosyst. (2015) | DOI

[Mollo and Hammer, 2017] S. Mollo; J. E. Hammer Dynamic crystallization in magmas, Mineral Reaction Kinetics: Microstructures, Textures, Chemical and Isotopic Signatures (W. Heinrich; R. Abart, eds.) (EMU Notes in Mineralogy), Volume 16, European Mineralogical Union and the Mineralogical Society of Great Britain & Ireland, Twickenham, 2017, pp. 1-46

[Mollo et al., 2010] S. Mollo; P. Del Gaudio; G. Ventura; G. Iezzi; P. Scarlato Dependence of clinopyroxene composition on cooling rate in basaltic magmas: Implications for thermobarometry, Lithos, Volume 118 (2010), pp. 302-312 | DOI

[Mollo et al., 2011a] S. Mollo; G. Lanzafame; M. Masotta; G. Iezzi; C. Ferlito; P. Scarlato Cooling history of a dike as revealed by mineral chemistry: A case study from Mt, Etna Volcano. Chem. Geol., Volume 283 (2011), pp. 261-273

[Mollo et al., 2011b] S. Mollo; K. Putirka; G. Iezzi; P. Del Gaudio; P. Scarlato Plagioclase-melt (dis)equilibrium due to cooling dynamics: implications for thermometry, barometry and hygrometry, Lithos, Volume 125 (2011), pp. 221-235 | DOI

[Morrison et al., 2020] A. Morrison; A. Whittington; B. Smets; M. Kervyn; A. Sehlke The rheology of crystallizing basaltic lavas from Nyiragongo and Nyamuragira volcanoes, D.R.C. Volcanica, Volume 3 (2020) no. 1, pp. 1-28 | DOI

[Mueller et al., 2010] S. Mueller; E. W. Llewellin; H. M. Mader The rheology of suspensions of solid particles, Proc. R. Soc. A: Math. Phys. Eng. Sci., Volume 466 (2010), pp. 1201-1228 | DOI

[Mysen and Richet, 2005] B. Mysen; P. Richet Silicate Glasses and Melts, Developments in Geochemistry, 10, Elsevier Science, Amsterdam, 2005, pp. 291-327

[Müller et al., 2004] T. Müller; L. P. Baumgartner; C. T. Foster; T. W. Vennemann Metastable prograde mineral reactions in contact aureoles, Geology, Volume 32 (2004), pp. 821-824 | DOI

[Nabelek, 2007] P. I. Nabelek Fluid evolution and kinetics of metamorphic reactions in calc-silicate contact aureoles from H 2 O to CO 2 and back, Geology, Volume 35 (2007), pp. 927-930 | DOI

[Perugini et al., 2008] D. Perugini; M. Petrelli; G. Poli A virtual voyage through 3D structures generated by chaotic mixing of magmas and numerical simulations: a new approach for understanding spatial and temporal complexity of magma dynamics, Vis. Geosci., Volume 13 (2008), pp. 1-24 | DOI

[Petrelli et al., 2016] M. Petrelli; K. El Omari; Y. Le Guer; D. Perugini Effects of chaotic advection on the timescales of cooling and crystallization of magma bodies at mid crustal levels, Geochem. Geophys. Geosyst., Volume 17 (2016) no. 2, pp. 425-441 | DOI

[Petrelli et al., 2018] M. Petrelli; K. El Omari; L. Spina; Y. Le Guer; G. La Spina; D. Perugini Timescales of water accumulation in magmas and implications for short warning times of explosive eruptions, Nat. Commun., Volume 9 (2018) no. 1, 770 | DOI

[Pistone et al., 2016] M. Pistone; B. Cordonnier; P. Ulmer; L. Caricchi Rheological flow laws for multiphase magmas: An empirical approach, J. Volcanol. Geotherm. Res., Volume 321 (2016), pp. 158-170 | DOI

[Polacci et al., 2018] M. Polacci; F. Arzilli; G. La Spina et al. Crystallisation in basaltic magmas revealed via in situ 4D synchrotron X-ray microtomography, Sci. Rep., Volume 8 (2018), 8377 | DOI

[Pupier et al., 2008] E. Pupier; S. Duchene; M. J. Toplis Experimental quantification of plagioclase crystal size distribution during cooling of a basaltic liquid, Contrib. Mineral. Petrol., Volume 155 (2008) no. 5, pp. 555-570 | DOI

[Putirka, 2008] K. D. Putirka Thermometers and barometers for volcanic systems, Reviews in Mineralogy and Geochemistry, Volume 69, Mineralogical Society of America, Virginia, USA, 2008, pp. 61-120 | DOI

[Roscoe, 1952] R. Roscoe The viscosity of suspensions of rigid spheres, J. Appl. Phys., Volume 2 (1952), pp. 267-269

[Roselle et al., 1997] G. T. Roselle; L. P. Baumgartner; J. A. Chapman Nucleation-dominated crystallization of forsterite in the Ubehebe peak contact aureole, California, Geology, Volume 25 (1997), pp. 823-826 | DOI

[Rosenberg and Handy, 2005] C. L. Rosenberg; M. R. Handy Experimental deformation of partially melted granite revisited: Implication for the continental crust, J. Metamorph. Geol., Volume 23 (2005), pp. 19-28 | DOI

[Rossi et al., 2019] S. Rossi; M. M Petrelli; D. Morgavi; F. Vetere; R. Almeev; R. Astbury; D. Perugini The role of magma mixing on pre-eruptive dynamics at the Aeolian Islands (Southern Tyrrhenian Sea, Italy), Lithos, Volume 324–325 (2019), pp. 165-179 | DOI

[Rusiecka et al., 2020] M. K. Rusiecka; M. Bilodeau; D. R. Baker Quantification of nucleation delay in magmatic systems: experimental and theoretical approach, Contrib. Mineral. Petrol., Volume 175 (2020) no. 5, 47 | DOI

[Saar et al., 2001] M. O. Saar; M. Manga; K. V. Cashman; S. Fremouw Numerical models of the onset of yield strength in crystal–melt suspensions, Earth Planet. Sci. Lett., Volume 187 (2001), pp. 367-379 | DOI

[Sehlke and Whittington, 2015] A. Sehlke; A. G. Whittington Rheology of lava flows on Mercury: An analog experimental study, J. Geophys. Res.: Planets, Volume 120 (2015), pp. 1924-1955 | DOI

[Sehlke and Whittington, 2020] A. Sehlke; A. G. Whittington Rheology of a KREEP analog magma: xperimental results applied to dike ascent through the lunar crust, Planet. Space Sci., Volume 187 (2020) no. 2020, 104941

[Sehlke et al., 2014] A. Sehlke; A. G. Whittington; B. Robert; A. Harris; L. Gurioli; E. Médard; A. Sehlke Pahoehoe to àà transition of Hawaiian lavas: an experimental study, Bull. Volcanol., Volume 76 (2014) no. 11, 876 | DOI

[Shapiro and Probstein, 1992] A. P. Shapiro; R. F. Probstein Random packings of spheres and fluidity limits of monodisperse and bidisperse suspensions, Phys. Rev. Lett., Volume 68 (1992), p. 1422 | DOI

[Shea and Hammer, 2013] T. Shea; J. E. Hammer Kinetics of cooling- and decompression-induced crystallization in hydrous mafic-intermediate magmas, J. Vocanol. Geotherm. Res., Volume 260 (2013), pp. 127-145 | DOI

[Smith and Asimow, 2005] P. M. Smith; P. D. Asimow Adiabat_1ph: A new public front-end to the MELTS, pMELTS, and pHMELTS models, Geochem. Geophys. Geosyst., Volume 6 (2005), Q02004 | DOI

[Swanson and Ottino, 1990] P. D. Swanson; J. M. Ottino A comparative computational and experimental study of chaotic mixing of viscous fluids, J. Fluid Mech., Volume 213 (1990), pp. 227-249 | DOI | MR

[Swanson, 1977] S. E. Swanson Relation of crystal-growth rate to the granitic textures, Am. Mineral., Volume 62 (1977), pp. 966-978

[Tripoli et al., 2019] B. Tripoli; M. Manga; J. Mayeux; H. Barnard The effects of deformation on the early kinetics of basaltic magmas, Front. Earth Sci., Volume 7 (2019), 250 | DOI

[Vand, 1948] V. Vand Viscosity of solutions and suspensions. I. Theory, J. Phys. Chem., Volume 52 (1948), pp. 277-299 | DOI

[Vetere and Holtz, 2020] F. Vetere; F. Holtz Rheological behavior of partly crystallized silicate melts under variable shear rate, Dynamic Magma Evolution (F. Vetere, ed.) (Geophysical Monograph Series), Wiley, Hoboken, NJ, USA, 2020 | DOI

[Vetere et al., 2013] F. Vetere; H. Sato; H. Ishibashi; R. De Rosa; P. Donato Viscosity changes during crystallization in shoshonitic magma: new insight on the lava flow emplacement, J. Mineral. Petrol. Sci., Volume 108 (2013) no. 3, pp. 144-160 | DOI

[Vetere et al., 2015] F. Vetere; G. Iezzi; H. Beherens; F. Holtz; G. Ventura; V. Misiti; A. Cavallo; S. Mollo; M. Dietrich Glass forming ability and crystallization behaviour of sub-alkaline silicate melts, Earth-Sci. Rev., Volume 150 (2015), pp. 25-44 | DOI

[Vetere et al., 2017] F. Vetere; S. Rossi; O. Namur; D. Perugini; D. Morgavi; V. Misiti; P. Mancinelli; M. Petrelli; C. Pauselli Experimental constraints on the rheology, eruption and emplacement dynamics of lavas from Mercury Northern Volcanic Plains, J. Geophys. Res.: Planets, Volume 122 (2017), pp. 1-17

[Vetere et al., 2019] F. Vetere; M. Murri; M. Alvaro; M. C. Domeneghetti; S. Rossi; A. Pisello; F. Holtz; D. Perugini Viscosity of Pyroxenite Melt and its Evolution during Cooling: an experimental approach, JGR-Planets, Volume 124 (2019), pp. 1451-1469

[Vetere et al., 2020] F. Vetere; A. Mazzeo; D. Perugini; F. Holtz Viscosity behaviour of silicate melts during cooling under variable shear rates, J. Non-Cryst. Solids, Volume 533 (2020), 119902 | DOI

[Vetere et al., 2021] F. Vetere; M. Petrelli; D. Perugini; S. Haselbach; D. Morgavi; A. Pisello; G. Iezzi; F. Holtz Rheological evolution of eruptible Basaltic-Andesite Magmas under dynamic conditions: The importance of plagioclase growth rates, J. Volcanol. Geotherm. Res., Volume 420 (2021), 107411 | DOI

[Vona and Romano, 2013] A. Vona; C. Romano The effects of undercooling and deformation rates on the crystallization kinetics of Stromboli and Etna basalts, Contrib. Mineral. Petrol., Volume 166 (2013) no. 2, pp. 491-509 | DOI

[Vona et al., 2011] A. Vona; C. Romano; D. B. Dingwell; D. Giordano The rheology of crystal-bearing basaltic magmas from Stromboli and Etna, Geochim. Cosmochim. Acta, Volume 75 (2011) no. 11, pp. 3214-3236 | DOI

[Webb and Dingwell, 1990a] S. L. Webb; D. B. Dingwell Non-newtonian rheology of igneous melts at high stresses and strain rates: experimental results for rhyolite, andesite, basalt, and nephelinite, J. Geophys. Res., Volume 95 (1990), pp. 15695-15701 | DOI

[Webb and Dingwell, 1990b] S. L. Webb; D. B. Dingwell The onset of non-Newtonian rheology of silicate melts, Phys. Chem. Minerals, Volume 17 (1990), pp. 125-132

[Wildemuth and Williams, 1984] C. Wildemuth; M. Williams Viscosity of suspensions modeled with a shear-dependent maximum packing fraction, Rheol. Acta, Volume 23 (1984), pp. 627-635 | DOI

[Yue and Bruckner, 1996] Y. Z. Yue; R. Bruckner On the different descriptions of the non-Newtonian viscosity (shear-thinning effect) of glass melts with respect to heat dissipation, Gerlands Beitr. Geophys., Volume 69 (1996) no. 6, pp. 179-185

[Zhang, 2008] Y. Zhang Geochemical Kinetics, Princeton University Press, Princeton, NJ, USA, 2008


Comments - Policy