Post-emplacement cooling and contraction of lava flows: InSAR observations and a thermal model for lava fields at Hekla volcano, Iceland

Wittmann, Werner, Sigmundsson, Freysteinn, Dumont, Stephanie and Lavallee, Yan ORCID: 0000-0003-4766-5758 (2017) Post-emplacement cooling and contraction of lava flows: InSAR observations and a thermal model for lava fields at Hekla volcano, Iceland. JOURNAL OF GEOPHYSICAL RESEARCH-SOLID EARTH, 122 (2). pp. 946-965.

Text Wittmann et al_JGR_final_withSuppl.pdf - Author Accepted Manuscript Download (7MB)

Abstract

<jats:title>Abstract</jats:title><jats:p>Lava flows contract as they cool, causing progressive subsidence of the flow surface. Here we study this process by measuring and modeling the deformation of emplaced lava flows and the surrounding substrate. The temporal trend of vertical lava movements was investigated using interferometric analysis of synthetic aperture radar (InSAR) images from the 1991 and 2000 Hekla eruptions, covering periods of 23 and 12 years, respectively. Data from six tracks from three satellites, including both ascending and descending passes, were used to create 99 interferograms, from which trends of accumulated subsidence and subsidence velocities were derived. Subsidence rates are similar for both lava flows and decay approximately exponentially from about 20 mm/yr 5 years after emplacement to about 2 mm/yr 15 years after emplacement. A one‐dimensional, semi‐analytical model was fitted to the observed subsidence rates, with subsidence due to phase change calculated analytically and subsidence due to thermal contraction calculated numerically using dilatometric constraints obtained experimentally. The initial thicknesses of the 1991 and 2000 lava fields, <jats:italic>D</jats:italic><jats:sub>1991</jats:sub> and <jats:italic>D</jats:italic><jats:sub>2000</jats:sub>, scaled thermal expansivity, <jats:italic>γ</jats:italic><jats:italic>α</jats:italic>, and thermal diffusivity, <jats:italic>κ</jats:italic>, are the crucial parameters influencing lava subsidence and subsidence rate. Inversion for these parameters reveal linear correlations between them. Best fitting results of inversions for <jats:italic>D</jats:italic><jats:sub>1991</jats:sub> range from 10 m to 27 m, for <jats:italic>D</jats:italic><jats:sub>2000</jats:sub> from 10 m to 30 m, <jats:italic>γ</jats:italic><jats:italic>α</jats:italic>= 9–24 ×10<jats:sup>−6</jats:sup> K<jats:sup>−1</jats:sup>, and <jats:italic>κ</jats:italic>= 1–7 ×10<jats:sup>−7</jats:sup> m<jats:sup>2</jats:sup> s<jats:sup>−1</jats:sup>.</jats:p>

Item Type:	Article
Depositing User:	Symplectic Admin
Date Deposited:	28 Feb 2017 10:41
Last Modified:	31 Oct 2023 02:07
DOI:	10.1002/2016JB013444
Related URLs:	Author Publisher
URI:	https://livrepository.liverpool.ac.uk/id/eprint/3006044