Comparison of 3D MT inversions for geothermal exploration: Case studies for Krafla and Hengill geothermal systems in Iceland

Abstract

a b s t r a c tThe magnetotelluric (MT) method is important for exploration of geothermal systems. The information on the Earth's resistivity obtained with MT methods has been valuable in imaging the hydrothermal alteration of such systems. Given its ability to recover complex resistivity models for the Earth, threedimensional (3D) MT inversion has become a common practice in geothermal exploration. However, 3D inversion is a time-consuming and complicated procedure that relies on computer algorithms to search for a model that can explain the measured data to a sufficient level. Furthermore, many elements of inversion require input from the practitioner, which can easily bias the results. Consequently, final 3D MT results depend on various factors, including the inversion code, the model mesh used to represent the Earth, data quality and processing, and constraints imposed during the inversion procedure.In this paper, to explore how this variability in 3D MT modeling impacts the final model, we invert MT data sets from the Krafla and Hengill geothermal areas in Iceland, using three different inversion codes. In each case, the modelers had the freedom to select a subset of the data and implement the inversion for the respective code in an optimized way. We compare the results from all the inversion codes, as well as consider the setup and assumptions made during the inversion process, all of which helps enhance the robustness and quality of the results. The comparison is done in multiple ways, using visual comparison of the recovered resistivity models, as well as comparing the structural similarities of the models by employing a structural correlation metric based on cross-gradients and other types of metrics for structural correlation. This approach highlights structures that are common in all three models, and implies that these structures are independent of the inversion code and necessary to fit the data.All modeling results from both Krafla and Hengill are consistent to first order, recovering a conductive layer on top of a resistive core typical of high temperature geothermal systems. For Hengill, the models show strong structural agreement, with all inversions recovering a moderately layered resistivity model but adding detail to previous work done in the area. Major differences are found in areas with coarse data coverage and hence questionable model resolution. Where the recovered structures in different models coincide, our confidence that these structures are well-constrained by the data is elevated, in spite of the different setup and assumptions in the codes these structures are required; so they can be interpreted in terms of geology with more certainty. Results from Krafla are not as consistent as results for Hengill, related in part to the Krafla data being nosier than the Hengill data. The models from Krafla have coinciding larger structures, but small-scale structures there are less coherent. One of the consistent structures in all the models is a conductive zone reaching from a depth of 5 km to...

Crustal Magmatism and Anisotropy Beneath the Arabian Shield—A Cautionary Tale

“…In contrast, relatively few studies have imaged melt in extensional settings outside of mid‐ocean ridges (e.g., Baba et al, ) or continental rifts with little or no Quaternary volcanism (e.g., Feucht et al, ). Those studies of active rift systems on land are primarily in Iceland (e.g., Miensopust et al, ; Rosenkjaer et al, ), within the Taupo Volcanic Zone in New Zealand (e.g., Bertrand et al, ; Heise et al, ) or within the East African Rift system.…”

The Electrical Structure of the Central Main Ethiopian Rift as Imaged by Magnetotellurics: Implications for Magma Storage and Pathways

“…A prominent example in a continental arc setting is the massive lower crustal Altiplano-Puna magma body in the central Andes (Comeau et al, 2015), which has been interpreted as the very hydrous source of multiple andesitic and dacitic lava flows in several volcanic centers. Recent MT studies in areas with rift-related volcanism include Hengill, Iceland (Árnason et al, 2010;Rosenkjaer et al, 2015) where a good electrical conductor at 4.5 km below sea level was detected; Eyjafjallajökull, Iceland (Miensopust et al, 2014;Sigmundsson et al, 2010) with a highly conductive feature at depth and the Taupo Volcanic Zone in New Zealand (Bertrand et al, 2013(Bertrand et al, , 2015Heise et al, 2010) where both the Rohaki geothermal field and Ohaaki volcano are associated with good conductors at depth.…”

Electrical resistivity tomography and time-domain induced polarization field investigations of geothermal areas at Krafla, Iceland: comparison to borehole and laboratory frequency-domain electrical observations

“…The hydrothermal plume at KH1 may relate to a deep low-resistivity anomaly observed at 1750 m depth in 3-D inversions of MT data in the same area, named Sandabotnaskarð, around the southern caldera rim (Árnason & Magnússon 2001). However, the shallow up-flow observed with ERT data, and correlated to geological and tectonic features, is not resolved with 1-D inversion of TEM data (Rosenkjaer et al 2015). Although hydrothermal activity has faded out, KH1 is still relatively hot today at the shallow depths investigated (170 • C at 200 m depth, Fig.…”