The 2014-2015 lava flow field at Holuhraun: Deriving physical properties of the lava using multi remote sensing techniques and datasets

The purpose of this thesis is to employ remote sensing to study lava flow products during the 2014-2015 eruption at Holuhraun, Iceland. Multimodal remote sensing techniques and datasets were applied and developed for three study themes (1) deriving thermal properties from satellite infrared remote sensing, (2) differentiating lava surface using airborne hyperspectral remote sensing, and (3) quantifying lava surface roughness from elevation model acquired by airborne LiDAR. In the first study, we present a new approach based on infrared satellite images to derive thermal properties within the lava field during eruption and then compare the results with field measurement during the 2014-2015 eruption at Holuhraun. We develop a new spectral index for Landsat 8, named the thermal eruption index (TEI), based on the SWIR and TIR bands (bands 6 and 10). The purpose of the TEI consists mainly of two parts: (i) as a threshold for differentiating between different thermal domains; and (ii) to apply dualband technique to determine the maximum subpixel temperature (Th) of the lava. Lava surface roughness effects are accounted for by using the Hurst exponent (H), which is estimated from radar backscattering profiles. A higher H (smooth surface) generates thinner crust and high thermal flux meanwhile a lower H (rough surface) generates thicker crust and lower thermal flux. The total thermal flux peak is underestimated compared to other studies, although the trend shows good agreement with both field observation and other studies. In the second study, we focus on retrieving the lava surface types contributing to the signal recorded by airborne hyperspectral at the very top surface of the 2014-2015 lava flow field at Holuhraun. For this purpose, an airborne hyperspectral image acquired at Holuhraun with an AisaFENIX sensor onboard a NERC (Natural Environment Research Council Airborne Research Facility) campaign. For sub-pixel analysis, we used the sequential maximum angle convex cone (SMACC) algorithm to identify the spectral image endmembers and the linear spectral mixture analysis (LSMA) method was employed to retrieve the abundances. SMACC and LSMA methods offer a fast selection for volcanic product segregation. However, ground-truthing of spectra is recommended for future work. In the third study, we perform both the topographic position index (TPI) and onedimensional Hurst Exponent to derived lava flow unit roughness on the 2014-2015 lava flow field at Holuhraun using both airborne LiDAR and photogrammetry topography datasets. The roughness assessment was acquired from four lava flow features: (1) spiny pāhoehoe, (2) lava pond, (3) rubbly pāhoehoe lava, and (4) inflated channel. The TPI patterns on spiny lava and inflated channels show that the intermediate TPI values correspond to a small slope indicating a flat and smooth surface. Lava pond is characterized by low to high TPI values and forms a wave-like pattern. Meanwhile, irregular transitions patterns from low to high TPI values characterize lava with rough blocky surfaces, i.e. rubbly pāhoehoe to 'ā'a flows and lobes and their margins. These lobes and margins may give the impression of having similar roughness as the ”rough” surface on meters scale since this is an “apparent” roughness. On centimeters scale these multitudes of lobes feature coherent and smooth surfaces because they are pāhoehoe. The surface roughness of these lava features falls within the H range of 0.30 ± 0.05 to 0.76 ± 0.04. The rubbly pāhoehoe / 'ā'a has the roughest surface and the inflated lava channel along with pāhoehoe feature the smoothest surfaces among these four surface types. In general, the Hurst exponent values in the 2014-2015 lava field at Holuhraun has a strong tendency in 0.5, which is compatible with results from other study of geological surface roughness. Overall, this project provides an important insights into the application of remote sensing for monitoring and studying active lava flow fields and the techniques developed here will benefit such work in future events.Tilgangurinn með verkefninu var að rannsaka hraunrennsli og landform er urðu til í eldgosinu norðan Vatnajökuls 2014-2015 og kennt við Holuhraun. Fjölþátta fjarkönnunartækni og gögn úr gervitunglum og flugvélum voru nýtt við úrvinnslu verkefnisins. Rannsóknin sneri að þremur megin þáttum: (1) greiningu á eðli varmaútstreymis frá Holuhrauni, út frá innrauðri varmageislun sem mæld er með gervitunglagögnum (2) aðgreining á mismunandi hraunyfirborði, út frá ofur-fjölrófs mælingum úr lofti, og (3) greiningu og flokkun á yfirborðshrjúfleika Holuhrauns út frá hæðarlíkani er aflað var með LiDAR settur upp í flugvél. Fyrsti þáttur beindist að eðli varmaútstreymis á meðan á eldgosi stóð. Stuðst var við gervitunglagögn og mælingar með FLIR tækni á meðan eldgosið stóð yfir. Afraksturinn er nýr hitastuðull fyrir Landsat 8 og greiningu á eldgosum, (TEI). Hitastuðullinn TEI er unninn út frá SWIR og TIR böndum Landsat 8 (bönd 6 og 10). Með TEI næst fram tvennt: (i) að greina þröskuld milli tveggja hitasviða; og (ii) að beita tvíbanda tækni til að greina hitastig innan hverrar myndeiningar (Th) af hrauninu. Hrjúfleiki hraunsins hefur áhrif á varmaútstreymi, og er gert ráð fyrir honum með því að reikna Hurst veldisstuðulinn (H) og eru reiknuð út frá radar endurkasti hraunyfirborðs. Hátt H einkennir flatt og mjúkt yfirborð og þunna skorpu á hrauninu, á meðan að lágt H einkennir úfið yfirborð, þykka skorpu og lága varmaútgeislun. Heildar varmaútgeislun með þessari aðferð er heldur vanmetin en ofmetin í samanburði við aðrar aðferðir. Hinsvegar er góð fylgni með mælingum í mörkinni og samanburðar aðferðum. Annar hluti rannsóknarinnar sneri að túlkun ofur-fjölrófsgreininga á yfirborði Holuhrauns. Flogið var yfir Holuhraun sumarið 2015 með ofur-fjölrófsmæli (AisaFENIX) um borð í flugvél frá NERC (Natural Environment Research Council Airborne Research Facility). Við greiningu á yfirborði innan hverrar myndeiningar var, (i) notast við aðferð runubundins hámarkshorns kúptrar keilu (SMACC) til að finna útmörk ofurrófs mælinganna, (ii) blönduð línulega rófgreining (LSMA) var nýtt til að greina styrk eða gnægð innan myndeiningar. SMACC og LSMA aðferðirnar bjóða upp á mjög hraða greiningu á yfirborði og útfellingum efna á yfirborðið. Hins vegar þarf að gera fleiri rófmælingar á staðnum, til þess að auka notkunnargetu aðferðarinnar í hraungosum framtíðarinnar. Þriðji þáttur rannsóknarinnar sneri að því að greina landfræðilega stöðuvísitölu (TPI) og einvíðan Hurst veldisvísi til að meta hrjúfleika á hinu endanlega yfirborði Holuhrauns. Við þessa greiningu var notast við LiDAR mælingu af hrauninu og hæðagrunn unninn út frá ljósmyndum. Hrjúfleikinn var metinn fyrir fjögur yfirborð sem einkenna hraunið: (1) broddahraun „spiny pāhoehoe lava“, (2) hrauntjörn „lava pond“, (3) klumpahraun „rubbly pāhoehoe lava“ og (4) upptjakkaða hrauntröð „inflated lava channel“. TPI fyrir yfirborð (1) og (4) gefur meðalgildi sem einkennist af litlum halla og flötu yfirborði. Hrauntjörnin einkennist af lágum og háum TPI gildum sem endurspegla bylgjukennt mynstur. Hinsvegar einkennast hrjúfustu yfirborðin (3) og hraunjaðrar af óreglulegu mynstri lágra og hárra TPI gilda. Hrjúfleika stuðull þessara yfirborða H, er á bilinu 0.30 ± 0.05 til 0.76 ± 0.04. Mestur er hrjúfleiki kubbahrauna og minnstur er hrjúfleiki þandar hrauntraðar. Hurts veldisvísir Holuhrauns er nærri 0.5, en það er í mjög góðu samræmi við niðurstöður fyrri rannsókna á jarðfræðilegum yfirborðum. Í heild gefur verkefnið mikilvæga sýn á notagildi fjarkönnunaraðferða við rauntímaeftirlit með hraungosum, m.a. með þróun stuðla sem munu nýtast við atburði framtíðar. Þá voru tengsl hraunmyndana við ýmsa eiginleika eldgosa skýrð, sem aftur getur gefið vísbendingar um eðli fyrri atburða

reflectance spectra measurements of mt etna a comparison with multispectral hyperspectral satellite

We present a collection of eight reflectance spectra representative of Mt. Etna volcano lava flows. The reflectance spectra were measured with a FieldSpecPro from 350 nm to 2500 nm during a fieldwork in June 2007. The reflectance has been compared with reflectance obtained by multispectral Advanced Space borne Thermal Emission and Reflection Radiometer (ASTER) and by hyper spectral EO1-Hyperion satellites. Prior the comparison, reflectance spectra have been convolved with ASTER and EO1-Hyperion spectral functions. The results show percentage errors in accordance to those present in literature in the ASTER SWIR range. However, some differences have been confirmed for the ASTER reflectance product (ASTER_07) in visible channels. Regarding EO1-Hyperion, a good agreement of reflectance against field measurement has been found resulting in 5% percentage maximum error in the VIS and up 30% in SWIR spectral range. The capacity of reproducing spectral feature by EO1-Hyperion has been checked on bright pixels (ice-snow) in the acquired image

Ten years of volcanic activity at Mt Etna: High-resolution mapping and accurate quantification of the morphological changes by Pleiades and Lidar data

Abstract The topography of Mt. Etna, Italy, is subjected to continuous modifications depending on intensity and magnitude of eruptions that frequently occur at the volcano summit and flanks. In order to make high-resolution maps of morphological changes and accurately calculate the overall volume of the erupted products (e.g., lava flows, tephra fall out, scoriae cones) in ten years, we have compared the altimetry models of Mt. Etna derived from 2005 Airborne Laser Scanning data and 2015 Pleiades stereo satellite imagery. Both models cover a common area of 400 km2 with spatial resolution of 2 m and comparable vertical accuracy (RMS

Investigations into the degassing and eruption mechanisms of Nyamuragira volcano, Democratic Republic of the Congo (Africa)

One of two active volcanoes in the western branch of the East African Rift, Nyamuragira (1.408ºS, 29.20ºE; 3058 m) is located in the D.R. Congo. Nyamuragira emits large amounts of SO2 (up to ~1 Mt/day) and erupts low-silica, alkalic lavas, which achieve flow rates of up to ~20 km/hr. The source of the large SO2 emissions and pre-eruptive magma conditions were unknown prior to this study, and 1994-2010 lava volumes were only recently mapped via satellite imagery, mainly due to the region’s political instability. In this study, new olivine-hosted melt inclusion volatile (H2O, CO2, S, Cl, F) and major element data from five historic Nyamuragira eruptions (1912, 1938, 1948, 1986, 2006) are presented. Melt compositions derived from the 1986 and 2006 tephra samples best represent pre-eruptive volatile compositions because these samples contain naturally glassy inclusions that underwent less post-entrapment modification than crystallized inclusions. The total amount of SO2 released from the 1986 (0.04 Mt) and 2006 (0.06 Mt) eruptions are derived using the petrologic method, whereby S contents in melt inclusions are scaled to erupted lava volumes. These amounts are significantly less than satellite-based SO2 emissions for the same eruptions (1986 = ~1 Mt; 2006 = ~2 Mt). Potential explanations for this observation are: 1) accumulation of a vapor phase within the magmatic system that is only released during eruptions, and/or 2) syn-eruptive gas release from unerupted magma. Post-1994 Nyamuragira lava volumes were not available at the beginning of this study. These flows (along with others since 1967) are mapped with Landsat MSS, TM, and ETM+, Hyperion, and ALI satellite data and combined with published flow thicknesses to derive volumes. Satellite remote sensing data was also used to evaluate Nyamuragira SO2 emissions. These results show that the most recent Nyamuragira eruptions injected SO2 into the atmosphere between 15 km (2006 eruption) and 5 km (2010 eruption). This suggests that past effusive basaltic eruptions (e.g., Laki 1783) are capable of similar plume heights that reached the upper troposphere or tropopause, allowing SO2 and resultant aerosols to remain longer in the atmosphere, travel farther around the globe, and affect global climates

The 2014–2015 Lava Flow Field at Holuhraun, Iceland: Using Airborne Hyperspectral Remote Sensing for Discriminating the Lava Surface

Publisher's version (útgefin grein)The Holuhraun lava flow was the largest effusive eruption in Iceland for 230 years, with an estimated lava bulk volume of ~1.44 km3 and covering an area of ~84 km2. The six month long eruption at Holuhraun 2014–2015 generated a diverse surface environment. Therefore, the abundant data of airborne hyperspectral imagery above the lava field, calls for the use of time-efficient and accurate methods to unravel them. The hyperspectral data acquisition was acquired five months after the eruption finished, using an airborne FENIX-Hyperspectral sensor that was operated by the Natural Environment Research Council Airborne Research Facility (NERC-ARF). The data were atmospherically corrected using the Quick Atmospheric Correction (QUAC) algorithm. Here we used the Sequential Maximum Angle Convex Cone (SMACC) method to find spectral endmembers and their abundances throughout the airborne hyperspectral image. In total we estimated 15 endmembers, and we grouped these endmembers into six groups; (1) basalt; (2) hot material; (3) oxidized surface; (4) sulfate mineral; (5) water; and (6) noise. These groups were based on the similar shape of the endmembers; however, the amplitude varies due to illumination conditions, spectral variability, and topography. We, thus, obtained the respective abundances from each endmember group using fully constrained linear spectral mixture analysis (LSMA). The methods offer an optimum and a fast selection for volcanic products segregation. However, ground truth spectra are needed for further analysis.The first author was supported by the Indonesia Endowment Fund for Education (LPDP) Grant No. 20160222025516, European Network of Observatories and Research Infrastructures for Volcanology (EUROVOLC), The European Facility for Airborne Research (EUFAR) and Vinir Vatnajökuls during his Ph.D. project.Peer Reviewe

Radiant and mass fluxes in multi-platform, multi-payload satellite-based volcano monitoring

The enormous amount of remote sensing (RS) data available today at a range of temporal and spatial resolutions aid emergency management in volcanic crises. RS provides a technological solution for bridging critical gaps in volcanic hazard assessment and risk mitigation. Detection and measurement of high-temperature thermal anomalies enable eruption monitoring and new lava flow propagation forecasts, for example. The accuracy of such thermal estimates relies on the knowledge of input parameters, such as emissivity - the efficiency with which surfaces radiate thermal energy at various wavelengths and temperatures. Emissivity is directly linked to the measurement of radiant flux and therefore affects the mass flux estimate as well as any model-based prediction of lava flow behaviour. Emissivity is not commonly measured across the range of volcanic lava compositions and temperatures, and it is generally assumed to have a constant value between 1.0 and 0.80 for basaltic lava. There is a lack of field and laboratory-based emissivity data for robust, more realistic modelling. To address this deficit, experiments on ‘aa’ lava samples were performed using data from Mount Etna (Italy), representing the range of its eruptive behaviour. In three sequential stages, emissivity was measured over the widest range of temperatures (294 – 1373 K) and wavelengths (2.17 - 15.0 μm) executable in the laboratory environment.The results show that emissivity is temperature, composition and wavelength dependent. Measured emissivity increases non-linearly with temperature decrease (cooling), exhibiting significant variations above 900 K with values considerably lower than the typically assumed 0.80. The measured and modelled emissivity values were applied to various remote sensing applications as input parameters for physical modelling of lava flows. This new evidence has significant impact on the computation of radiant heat flux from spaceborne data, as well as on modelling of lava flow ‘distance-to-run’ simulations. Furnished with improved input parameters (multicomponent emissivity), the novel approach developed here can be used to test an improved version of an unsupervised multi-platform, multi-payload volcano monitoring system

An Analysis of Proximal Volcanic Ash Emissions

Volcanic ash is a product of explosive volcanic eruptions and refers to those particles < 2 mm in size that are ejected from a volcano. Once erupted into the atmosphere, ash can be transported vast distances. Satellite remote sensing has provided us with the tools to map and monitor these plumes. These methods are restricted in viewing only that portion that is optically transparent, due to the requirement to observe the interaction of ground upwelling radiation with the plume. For this study, methods have been developed to map not only the opaque portion of the plume, but in higher spatial resolution than has been attempted previously. This has been done by using the unique emissivity spectra that are produced by materials in the thermal infrared portion of the electromagnetic spectrum. By using an end-member linear unmixing model with a spectral emissivity library of different volcanic ash types, opaque plume bearing pixels of Advanced Spaceborne Thermal Emission and Reflection Radiometer (ASTER) data have been processed. These end-member results have shown the applicability of this technique, with sampled volcanoes mapped with very low root mean square (RMS) errors. Volcanoes with unknown composition were also mapped with the library with varying degrees of success. Further expansion of the laboratory spectral library will allow more accurate assessment of these volcanoes. The collection of per pixel emissivity spectra of opaque plumes was also attempted using an experimental multispectral FLIR thermal camera apparatus, in order to collect data on the rising volcanic ash column. These data showed that the camera can produce emissivity spectra, however further calibration and correction is required in order to make a volcanological assessment of the eruptive products. Finally, a method of better tracking disconnected ash plumes is assessed. In cases where a plume is detected disconnected from or distally from the source volcano, the application of geostatistical methods to backward trajectory model data can yield the source location. This could then be used to helped point the ASTER sensor off-axis towards erupting targets, thus providing a greater quantity of data to be analyzed