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

Abstract

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.