DOI: 10.5281/zenodo.7701194
title: Automated processing of aerial imagery for geohazards monitoring: Results from Fagradalsfjall eruption, SW Iceland, August 2022

description: <strong>1-</strong> <strong>Dataset Summary</strong><br> Here we present a dataset of DEMs (Digital Elevation Models), orthomosaics, and lava area outlines for the August 2022 eruption at Fagradalsfjall, SW Iceland. The dataset consists of: (1) five aerial surveys collected over the course of the August 2022 Fagradalsfjall eruption, (2) one survey carried out on 14 August 2022 using Pléiades satellite stereo images, and (3) a larger aerial survey, covering the 2021 and 2022 eruption sites in late September 2022 after the volcanic activity concluded. <strong>2- Background</strong> The volcano at Fagradalsfjall, SW-Iceland, began erupting on 3 August 2022 at 13:20 following 10 months of quiescence. As part of the response plan, a series of photogrammetric surveys were conducted in rapid, operational mode throughout the duration of the eruption. Subsequent production of data products for natural hazards monitoring (lava maps, lava volumes, effusion rates) were calculated within hours and reported to the Icelandic Civil Defense, following a similar approach that described in Pedersen et al., 2022a and in Gouhier et al., 2022. At the start of the 2022 eruption, GCPs had not yet been placed around the new fissure, but reference data (orthomosaics and DEMs) which had been georeferenced using targets measured with differential GNSS existed of the eruption site from September 2021 from Pedersen et al. (2022b) were available to use as a reference in the new workflow instead of GCPs. Due to the urgent need from authorities for information about the new eruption, a processing method that avoids the time-consuming task of manual GCP selection using a reference image for georeferencing was preferable in this instance. Besides the acquisition of aerial photographs, the CIEST2 initiative was also re-activated to collect Pléiades stereo images in emergency mode (Gouhier et al., 2022). <strong>3 – Overview of data collection</strong> Table 1 contains the overview of the surveys collected and presented in this repository. Table 1. Summary of surveys included in this dataset, by survey date. <strong>Date & Time</strong><br> <strong>YYYYMMDD HH:MM</strong> <strong>Sensor</strong> <strong>Platform</strong> <strong>Flight alt.</strong><br> <strong>(m asl)</strong> <strong>Images</strong> <strong>Surveyed</strong><br> <strong>km<sup>2</sup></strong> 20220803 17:05 A6D TF-203* ~ 850 46 4 20220804 11:00 A6D TF-203 ~ 2100 32 35 20220813 09:00 A6D TF-203 ~ 750 123 9 20220814 13:00 Pléiades PHR1B n/a 2 14 20220815 08:15 A6D TF-203 2100 20 23 20220816 10:06 A6D TF-203 2100 19 26 20220926 12:00 A6D TF-BMW** 2100 ~20 18 * TF-203: Savannah S aircraft ** TF-BMW: Vulcanair P68 Observer 2 aircraft, operated by Garðaflug ehf. <strong>4- Methods</strong> <strong>4.1 Processing of the aerial photographs from 3-16 Aug 2022</strong><br> Throughout the eruption, aerial surveys were conducted using a Hasselblad A6D 100 MP camera with 35 mm focal lens, from a height of 750 – 2,100 m above ground over the active lava field from an ultralight aircraft with a window in the bottom to allow for vertical photos to be taken (see supplement of Pedersen et al., 2022a for details and images of the setup). The camera was manually triggered to give ~70% overlap, and approximate flight lines were prepared beforehand for use with a handheld GPS during the flight to give ~30 % side overlap.<br> <br> An automated processing pipeline was created in python, which leverages tools from the Ames Stereo Pipeline (ASP, Shean et al., 2016) and Agisoft Metashape stand-alone Python API (v. 1.8.4). The processing and georeferencing of the aerial data were done in three steps, with all steps being automated except for the digitization of lava outlines. First, using a very high-resolution reference orthomosaic and DEM created in September 2021 and georeferenced with ground control points (Pedersen et al., 2022b), interest points (IPs) in each image were matched with the reference dataset, using the ASP routine <em>ipfind. </em>This created GCPs for each image over stable terrain. Second, hillshades were created from both the reference DEM and the source dataset DEM and matches in IPs were found in both, creating a second round of ground control points to refine the georeferencing of the entire block. Finally, the alignment of the source DEM was refined using the <em>dem_align</em> (demcoreg) protocol from Shean et al. (2016) by applying a bulk linear shift in X, Y and Z which minimizes the vertical difference in stable terrain between the source and reference DEM. <strong>4.2 Processing of the Pléiades stereo images</strong><br> The Pléiades stereo images were processed using the Ames Stereo Pipeline, using the general workflow of <em>mapproject </em>and <em>parallel_stereo </em>(e.g., Deschamps-Berger et al., 2020). The <em>parallel_stereo </em>routine used default arguments, plus the following arguments: <em>--stereo-algorithm asp_mgm -t rpcmaprpc --corr-seed-mode 3 --corr-max-levels 2 --cost-mode 3 --subpixel-mode 9 --corr-kernel 7 7 --subpixel-kernel 15 15</em> We used the DEM from 4 Aug 2022 as the reference for <em>mapproject </em>and for the final DEM co-registration applied to the produced Pléiades DEM. <strong>4.3 Processing of the 26 September 2022 dataset</strong><br> The survey from 26 September 2022 was collected and processed using direct georeferencing from an on-board GPS antenna. The final alignment of the block was refined using the dem_align (demcoreg) protocol from Shean et al. (2016) by applying a bulk linear shift in X, Y and Z which minimizes the vertical difference in stable terrain between the source and reference DEM. Because this survey covered a much larger area, the reference DEM for the final coregistration was the ÍslandsDEM v.1.0 (Landmælingar Íslands, 2022). <strong>4.4. Maps of the lava outlines, lava thickness, lava volume, Time Average Effusion Rate (TADR)</strong><br> For each survey, a differential DEM (dDEM) showing elevation changes since the 2021 eruption was created by subtracting the reference DEM (ÍslandsDEM v.1.0, which includes the post-eruption DEM from Pedersen et al., 2022a) from the source DEM. Lava outlines, lava thickness lava volume, TADR and uncertainties were calculated using the methods described in Pedersen et al., 2022a. Table 2 summarizes calculations from this dataset. Table 2. Summary of survey results calculated from August 2022 Fagradalsfjall eruption DEMs and orthomosaics. <strong> Date Start</strong> <strong>Date End</strong> <strong>Time </strong> <strong>Difference</strong> <strong>Lava </strong> <strong>Area<sup>*</sup> End </strong> <strong>(km<sup>2</sup>)</strong> <strong>dh<sup>**</sup> </strong> <strong>(m)</strong> <strong>Volume<sup>+</sup></strong> <strong>End<br> (1e+6 m<sup>3</sup>)</strong> <strong>TADR<sup>++</sup><br> (m<sup>3</sup>/s)</strong> 20220803<br> 13:20 20220803<br> 17:05 0d 03h 45m 0.07 5.88 0.43 ± 0.03 32.1 ± 1.5 20220803<br> 17:05 20220804<br> 11:00 0d 21h 40m 0.14 11.13 1.57 ± 0.05 17.7 ± 0.8 20220804<br> 11:00 20220813<br> 09:00 8d 22h 00m 1.27<sup>#</sup> 6.90 10.33 ± 0.6 11.4 ± 0.7 20220813<br> 13:08 20220814<br> 13:00 0d 23h 52m 1.24 7.28 10.62 ± 0.70 2.8 ± 0.8 20220814<br> 13:00 20220815<br> 08:15 0d 19h 15m 1.26 7.46 10.99 ± 0.55 4.1 ± 0.8 20220815<br> 08:15 20220816<br> 10:16 1d 2h 01m 1.28 7.49 11.13 ± 0.53 2.0 ± 0.7 20220816<br> 10:16 20220821<br> 06:00<sup>##</sup> 4d 19h 44m 1.28 7.69 11.39 ± 0.44 0.653 ± 0.10 <sup>*</sup>Total area of the lava field since 2022-08-03 before activity started. <sup>**</sup>dh end is the mean thickness of the lava flow-field in the end of the given period. <sup>+</sup>Volume erupted since 2022-08-03 before activity started. <sup>++</sup>Time-averaged discharge rate for the given period <sup>#</sup>Extrapolated value. Survey does not cover entire active lava area. <sup>##</sup>End Time: 21 August 2022, 6:00. This time deduced from field observations from members of the Institute of Earth Sciences, University of Iceland. Values calculated from 26 September 2022 dataset. Figures and visual summaries of the processing method, resulting lava volumes, and uncertainties can be found in this poster: Fagradalsfjall August 2022.<br> <br> Orthomosaics from this dataset are viewable online at: https://atlas.lmi.is/mapview/?application=umbrotasja <strong>Data naming conventions:</strong> Data type: DEM, ortho, outline, diffDEM, lavafree, lava Acquisition date: YYYYMMDD_HHMM Platform/Sensor for data collection: Pléiades (PLE), Hasselblad A6D from aircraft (A6D) Resolution: 2x2 m (DEMs) and 30x30 cm (Orthomosaics) Folders (by survey): YYYYMMDD_HHMM_platform (during eruption) or 'posteruption'_platform (sensors/platform: A6D or PLE) <strong>Data Specifications:</strong> Cartographic projection: ISN93 / Lambert 1993 (EPSG:3057, https://epsg.io/3057) Origin of Elevation: meters above GRS80 ellipsoid (WGS84) Raster data format: GeoTIFF Raster compression system: ZSTD (http://facebook.github.io/zstd/) Vector data format: GeoPackage (https://www.geopackage.org/) Pléiades dataset includes only DEMs because the Pléiades ortho imagery is for licensed use only. Please contact the authors for further information on this.
description: Costs of overflights were paid by the Icelandic State as a part of emergency response coordinated by the National Commissioner of the Icelandic Police, Department of Civil Protection. Other funding was provided by the Institute of Earth Sciences and the National Land Survey of Iceland. The Icelandic Institute of Natural History provided camera equipment, and Físfélag Reykjavíkur provided lightweight aircrafts and pilots for the project. Pléiades images were provided by the CIEST 2 (CNES) initiative and through the Icelandic Volcanoes Supersite project supported by the Committee on Earth Observing Satellites (imagery Pléiades©CNES2021, distribution AIRBUS DS).
description: {"references": ["Deschamps-Berger, C., Gascoin, S., Berthier, E., Deems, J., Gutmann, E., Dehecq, A., Shean, D. & Dumont, M. (2020). Snow depth mapping from stereo satellite imagery in mountainous terrain: evaluation using airborne laser-scanning data. The Cryosphere, 14(9). https://doi.org/10.5194/tc-14-2925-2020", "Gouhier, M., Pinel, V., Belart, J. M. C, De Michele, M., Proy, C., Tinel, C., Berthier, E., Gu\u00e9henneux, Y., Gudmundsson, M.T., \u00d3skarsson, B.V. & Gremion, S. (2022). CNES-ESA satellite contribution to the operational monitoring of volcanic activity: The 2021 Icelandic eruption of Mt. Fagradalsfjall. Journal of Applied Volcanology, 11(1). https://doi.org/10.1051/0004-6361/202141613", "Pedersen, G. B. M., Belart, J. M. C., \u00d3skarsson, B. V., Gudmundsson, M. T., Gies, N., H\u00f6gnad\u00f3ttir, T., et al. (2022a). Volume, effusion rate, and lava transport during the 2021 Fagradalsfjall eruption: Results from near real-time photogrammetric monitoring. Geophysical Research Letters, 49, e2021GL097125. https://doi.org/10.1029/2021GL097125", "Pedersen, G. B. M., Belart, J. M. C., \u00d3skarsson, B. V., Gudmundsson, M. T., Gies, N., H\u00f6gnad\u00f3ttir, T., Hjartad\u00f3ttir, \u00c1. R., Pinel, V., Berthier, E., D\u00fcrig, T., Reynolds, H. I., Hamilton, C. W., Valsson, G., Einarsson, P., Ben-Yehoshua, D., Gunnarsson, A, & Oddsson, B. (2022b). Digital Elevation Models, orthoimages and lava outlines of the 2021 Fagradalsfjall eruption: Results from near real-time photogrammetric monitoring (v1.1) [Data set]. Zenodo. https://doi.org/10.5281/zenodo.6598466Gouhier et al., 2022", "Shean, D. E., Alexandrov, O., Moratto, Z. M., Smith, B. E., Joughin, I. R., Porter, C., & Morin, P. (2016). An automated, open-source pipeline for mass production of digital elevation models (DEMs) from very-high-resolution commercial stereo satellite imagery. ISPRS Journal of Photogrammetry and Remote Sensing, 116. http://dx.doi.org/10.1016/j.isprsjprs.2016.03.012"]}