*3.3. Resistivity*

Extensive resistivity surveys for geothermal exploration (dating back to the early 1970's) have been carried out in Krafla. High-temperature geothermal systems have a characteristic resistivity structure due to geothermal alteration [28,29]. Below resistive near-surface unaltered rocks, they have a shallow low-resistivity "cap" (clay cap) with conductive smectite and zeolite alteration minerals formed in the temperature range of 50–230 ◦C. The low-resistivity cap is underlain by a resistive "core" with resistive chlorite and epidote alteration minerals formed at temperatures above 230 ◦C.

Figure 9 shows a resistivity map at 200 m a.sl. (at a depth of about 300 m) based on 1D inversion of central-loop TEM soundings [30]. Areas where high resistivity (resistive core) has appeared below the conductive clay cap are gridded red. The figure also shows the 110 ka caldera, the inferred inner caldera and the transverse low-density structure.

**Figure 9.** Resistivity at 200 m a.sl. (about 300 m depth) based on 1D inversion of TEM soundings. Areas where higher resistivity is observed below low resistivity are shown as red gridded (black dots are TEM soundings (reproduced from [30])). The thick black hachured lines show the 110 ka caldera. Thinner black hachured lines show buried inner caldera and grey lines mark a buried transverse structure inferred from gravity. Coordinates are UTM WGS84 zone 28 in km.

Figure 9 shows that high-temperature geothermal alteration is confined within the inner caldera, except in the western part of the southern arm of the fissure swarm and at the rim of the outer caldera in the southern end of <sup>H</sup>ágöng, east of Mt. Krafla. Both of these anomalies have been drilled into and geothermal alteration was found in accordance with the resistivity, but the present temperature is far below that responsible for the alteration. Therefore, cooling has occurred at these places outside the inner caldera. Figure 9 shows that no mature high-temperature geothermal system has developed

outside the inner caldera except for those mentioned. It also shows that, at this depth, high-temperature alteration is absent in the transverse structure, which means a lower temperature.

To study the deeper resistivity structure, most of the TEM sounding sites (black dots in Figure 9) were later visited for MT soundings. Figure 10 shows slices through a resistivity model resulting from 3D inversion the MT soundings (static shift corrected by the TEM).

**Figure 10.** A resistivity model for Krafla from 3D inversion of MT data. Left: resistivity map at 1250 m bl. sl. The two calderas are shown as black hatched lines, purple lines mark the S-wave shadows, grey dots are earthquake epicentres and blue lines are faults and fissures. The location of the well IDDP-1 is shown as a black star. Right: 3D view from the west. The white surface shows the upper boundary of the western S-wave shadow at 3 km depth (reproduced from [31]). Coordinates are UTM WGS84 zone 28 in km.

The figure to the left, shows the resistivity at −1250 m a.sl. (about 1700–1800 m depth). It shows deep, low-resistivity bodies, mainly within the inner caldera except for a relatively low resistivity at the south-east rim of the outer caldera. The low-resistivity bodies within the inner caldera border an ESE-WNW high-resistivity structure, which roughly coincides with the S-wave shadows. This high-resistivity structure hosts most of the seismicity and the main production field is in its central and eastern part. The figure to the right shows the shallow conductive clay cap is fairly well in agreemen<sup>t</sup> with the 1D inversion of TEM and drilling, and that the two low-resistivity bodies extend down to about 6 km depth. Around that depth, the northern one seems to start to deviate to the north.

## *3.4. Seismic Studies*

ÍSOR operates a local seismic monitoring network in Krafla for Landsvirkjun (National Power Company), which operates the Krafla power plant. Even though there have been no signs of volcanic unrest since the end of Krafla Fires, Krafla turns out to be very seismically active in the central part of the inner caldera. Figure 11 shows earthquakes located in Krafla from November 2015 to November 2016 [32], a rather typical annual seimicity. There are five main clusters, one just SW of Mt. Krafla, just north of the transverse structure, two at injection wells K-26 and IDDP-1 (probably mostly due to injection) and two at and north of Leirhnjúkur. The depth sections show that maximum depth of the earthquakes, i.e., the brittle/ductile boundary, which for basalt, is estimated in the range of 600–800 ◦C [33], is at about a 2–2.3 km depth north of the transverse structure. There is also a small cluster within the transverse structure south of Leirhnjúkur. The depth sections show that seismicity there extends much deeper than in the north or to about 3.5 km. This indicates lower upper crustal temperatures within the transverse structure. This is in agreemen<sup>t</sup> with the distribution of high-temperature alteration at a 300 m depth shown in Figure 9, i.e., a lack of high-temperature alteration within the transverse structure.

**Figure 11.** Earthquakes located in Krafla from November 2015 to November 2016. The map also shows the location of reinjection wells (K-26 and IDDP-1) and the transverse gravity low (grey lines) (modified from [32]).

Recent seismic tomography study of the Krafla area [34] shows some interesting results. Figure 12 (left) shows the estimated P-wave velocity (Vp) at 0.5 km bl.sl. (at about 1 km depth) and (right) at 2 km depth below sea level (about 2.5 km depth). At 0.5 km bl.sl., the tomography indicates low Vp within the transverse structure and higher Vp north of it. At 2 km bl.sl., the situation is the opposite, i.e., high Vp in the southern part of the calderas, but low Vp and low Vp/Vs ratio within the inner caldera north of the transverse structure and above the S-wave shadows. Schuler et al. (2015) sugges<sup>t</sup> that this anomaly reflects rocks containing superheated steam.

**Figure 12.** P-wave velocity (km/s) at 0.5 km bl.sl. (left) and at 2 km bl.sl. (right) according to recent seismic tomography (reproduced from [34]). Black hatched lines show the two calderas, grey lines mark the low-density transverse structure, blue lines show faults and fissures and the purple lines mark the S-wave shadows. Coordinates are UTM WGS84 zone 28 in km.

#### **4. The Geothermal System(s)**

Through the years, a total of 43 deep exploration and production wells have been drilled in Krafla, some of which are (or have been) used for injection. Figures 13 and 14 show their location and well tracks of deviated wells. The main well field is located in the eastern part of the inner caldera. Through drilling, the stratigraphic structure at depth has been outlined as well as the physical conditions of the geothermal system [11,12,35,36].

**Figure 13.** Drilled wells in Krafla (red dot on the little inlet map to the right) and the three main production fields (right) with different thermal character. The fields are: Vtísmór and Leribotnar (blue frame), southern and western slopes of Mt. Krafla (red frame) and Hvíthólar (purple frame). To the left are shown characteristic temperature profiles in wells in each production area and a well drilled in the fissure swarm west of Hvíthólar. For well KS-01, east of Hvíthólar, see text.

**Figure 14.** Left; The locations and tracks of wells in Krafla (green lines) and the location of the cross-section to the right (black line). Right; A resistivity cross-section with alteration mineralogy and formation temperature in wells.

The drilling has revealed that different "subfields" in Krafla have quite different thermo-hydraulic characteristics (Figure 13). In part of the production area, Leirbotnar and Vítismór (blue frame in Figure 13), the geothermal system has two different thermo-hydraulic characters at different depths. Temperature profiles show that above 1000–1200 m in depth, there is a system with two phase (boiling) conditions close to the surface but then almost isothermal with depth, at about 200 ◦C. Below the isothermal part, the temperature rises almost linearly up to the boiling point curve (boiling point versus depth curve) in an underlying two-phase system. Wells in the southern and western slopes of Mt. Krafla (red frame in Figure 13) show two phase conditions from maximum drilled depth (2–2.5 km) to the surface and wells at Hvíthólar (purple frame in Figure 13) show a temperature inversion. Well KS-01 in Sandabotnaskarð, east of Hvíthólar, has a temperature profile (not shown) following the boiling point curve from the bottom (2500 m depth) and up to about 1000 m in depth, but temperatures below the boiling point curve at shallower depths (ÍSOR database). Finally, well KV-1, in the active fissure swarm west of Hvíthólar, shows extensive cooling. Resistivity surveys show high-temperature alteration there (Figure 9) and the alteration mineralogy of the well showed that temperatures there had been much higher than the current measured temperature.

Figure 14 shows a resistivity cross-section from Vítismór in the NW to the southern slopes of Mt. Krafla. Superimposed on the section is the alteration mineralogy and estimated formation temperature in wells close to the section. The figure shows an excellent correlation between resistivity and alteration, i.e., low resistivity at shallow depth in the smectite and zeolite alteration and increasing resistivity in the chlorite and chlorite/amphibole alteration [28,29]. The isotherms in the wells in the western part, west of Hveragil (a purple line in Figure 4), and also in the eastern most well (K-18), show much lower temperatures than the ones causing the alteration, so cooling has taken place.

#### **5. New Conceptual Model**

The complexity of the Krafla volcano and its geothermal system(s) has puzzled geoscientists for decades. Here, an attempt is made to shed some light on this complexity by putting it in context with the geological structures, tectonics and geological history discussed above. The first conceptual model of the geothermal system in Krafla was presented in 1977 [37]. It was based on limited data, mainly surface studies and the results for the first eleven wells in the area [35,38]. In that model, it was assumed that the shallow and the deep thermal regimes in the Vítismór-Leirbotnar area are separated by an impermeable cap-rock. It was assumed that there was a main up-flow from a deep two phase geothermal system through Hveragil (bold purple line in Figure 4) just west of Mt. Krafla and that the near-isothermal upper system in the west was due to a westward lateral flow from the up-flow under Hveragil. The difference in the contamination by volcanic gases during the Krafla Fires, discussed above, was considered to support this conceptual model [5]. This conceptual model has been considered valid until recently [39]. The author did, however, find this hypothesis unlikely and in the most recent review report, on the Krafla geothermal system [12], a quite different scenario was proposed. There is no obvious candidate for an impermeable cap-rock in the lithology in the wells in the Vítismór-Leirbotnar area, but the change in thermal character occurs at depths where intrusive rocks start to dominate. The area is within the presently active part of the fissure swarm and is likely to have very anisotropic permeability, i.e., high permeability along the fissure swarm but much lower permeability perpendicular to the faults and fissures. It is, therefore, unlikely that lateral flow, above a cap-rock from an up-flow in the east, would be to the west, perpendicular to the high permeability. An alternative explanation of the two thermo-hydraulic regimes is given below.

#### *5.1. The Role of Permeability*

Permeability plays a crucial role in the existence and development of volcanic high-temperature geothermal systems. Model calculations show that, in order for a magma intruded into the shallow crust to produce a two-phase geothermal system, permeability has to be within the rather narrow range of 0.5–5 mD (1 mD = 10−<sup>15</sup> m2) [40]. If the permeability is lower, the intrusion cools over a long time by heat conduction, but if the permeability is higher, the intrusion cools rapidly by vigorous one-phase (water) convection.

Model calculations show that if a magma intrusion is emplaced in rocks with permeability in the above range, a two-phase convective geothermal system is formed in time of the order of 1000 years and, if no further intrusions occur, the lifetime of such as system is of the order of 10 ka [40] (though depending on the size of the intrusion).

Figure 15 shows the results of two-dimensional (2D) numerical modelling using the HYDROTHERM programme [41,42]. The figure shows the thermo-hydraulic states of a geothermal system developed after the emplacement of an 800 m wide (infinitely long) "dike" intrusion with a temperature of 1100 ◦C, extending up to 2 km below the surface into rocks with a temperature gradient of 100 ◦C/km. The modelling takes into account the latent heat of the magma and the permeability created when the magma solidifies and contracts.

The figure to the left shows the developed geothermal system 5 ka after the dike injection into a host rock with a permeability of 1 mD. The figure shows a convecting two-phase geothermal system with up-flowing water and steam above the intrusion and a down flow of colder water to the sides down to the heat source. The panel in the upper-right corner shows calculated temperature profiles for hypothetical wells. Wells within the range of the dike (red and yellow) have temperature profiles corresponding to the boiling point curve (boiling point temperature vs. depth curve), while wells at the edge of or outside the dike (light blue and purple) roughly follow the boiling point curve at shallow depths but show temperature inversion at greater depths due to the in-flow of colder water towards the heat source.

Near the end of the lifetime of the system (of the order of 10 ka, if the heat source is not renewed), the heat source has lost most of its heat, and cold downwards flowing water can no longer gain enough heat to sustain two-phase conditions and the system cools down at depth. Two-phase conditions will, however, still prevail (for some hundreds to thousand years) at a shallower depth in a slowly upwards migrating "bubble", leading to temperature inversion with depth above the cooled heat source. Similar results, using di fferent modelling software have been reported [43].

The right part of Figure 15 shows the thermo-hydraulic situation 3 ka after a similar dyke is emplaced in a host rock of permeability 1 mD, underlying a 1 km thick pile of a much higher permeability (50 mD in this case). The figure shows a two-phase system in the 1 mD host-rock but in the more permeable rocks above, a second and vigorous convection cell develops. The temperature profiles in the upper-right corner show that above the centre of the dike, the temperature follows the boiling point curve near the surface and then becomes almost isothermal until just above the permeability contrast, where it rises sharply towards the boiling point curve. It is, therefore, suggested that the two di fferent thermal regimes in the Vítismór-Leirbotnar field are not due to cap-rock and lateral flow from east, but due to a much higher permeability in the fissure swarm above the intrusions at about 1 km depth. The model can also explain the di fference in contamination by volcanic gases in the upper isothermal and the deeper two-phase part during the Krafla Fires, discussed above. The vigorous convection in the high permeability isothermal part above the deeper two-phase part would dilute the gases in fluids coming from below.

It may be questioned whether the simple 2D modelling shown here is justified in this case. However, it will be argued below that the main heat sources powering the geothermal system in Krafla are ESE-WNW trending intrusions and a dike complex in the inner caldera north of the low-density transverse structure, nearly perpendicular to the fissure swarm. Since permeability above the intrusions is mainly along faults and fissures perpendicular to the heat source(s), a 2D modelling is considered as a valid approximation.

A permeability of 50 mD in the upper 1000 m may seem unrealistically high. The Svartsengi geothermal system in the Reykjanes peninsula SW Iceland, has a wide spread confined nearly isothermal reservoir of about 235–240 ◦C from about 700 to 2000 m depth [44]. In numerical reservoir modelling of the system a permeability of about 220 mD is needed to simulate this isothermal system [45].

What, from resistivity, appears as one system in the northern part of the inner caldera is actually divided into two very di fferent parts. East of Hveragil, in Vesturhlíðar and Suðurhlíðar, there are, generally speaking, two-phase conditions from the depth of drilling and to shallow depths, except in the easternmost well (K-18), which shows cooling in the upper parts although hot near the bottom. West of Hveragil (bold purple line in Figure 4), the geothermal system is divided into two very di fferent parts, a deep two-phase system, below 1 km, and an upper, almost-isothermal, one-phase condition of about 200 ◦C, except for in the near-surface, where two-phase conditions are reached. Alteration shows, however, that sometime in the past, two-phase conditions existed in the upper part.

As discussed above, until about 8 ka ago, the rifting was in the presently active eastern part of the fissure swarm with considerable volcanism, mainly within the caldera (Mt. Krafla) and to the south. Then, the rifting was shifted to the west for about 5 ka. There seems to have been amazingly little extrusive volcanic activity during that period, the main event being a phreatic eruption from the explosion crater Hvannstóð, about 5 ka ago [7] and resistivity (Figures 9 and 10) shows that no pronounced geothermal system was formed in the western part of the 110 ka caldera.

Here, it is proposed that during the 5 ka, when the rifting was in the western part, intrusions continued to be emplaced in the eastern part and a two-phase geothermal system was formed and sustained there. About 3 ka ago, the spreading shifted back to the presently active part of the fissure swarm. This greatly increased the permeability in the upper part, west of Hveragil, leading to a vigorously convecting isothermal system. At greater depth, continuing intrusions maintained lower permeability with two-phase conditions. Fissures and faults east of Hveragil were not activated when the spreading moved back to the east and two-phase conditions still prevail there.

During the time when the spreading took place in the western part of the 110 ka caldera, intrusions did probably take place in the southwestern part of the inner caldera, south of the transverse low-density structure, and a geothermal system developed. When the rifting turned back east, the permeability increased drastically and cooling took place. Intrusions have not been frequent enough to maintain the heat source(s) and an active geothermal system there. The temperature profile of KV-1 (Figure 13) shows a temperature maximum at a shallow depth and lower temperatures below. This temperature maximum is likely due to a lateral flow in shallow permeable layers. It could be that a small geothermal system still exists further to the west inside the inner caldera and south of the transverse structure.

Hvíthólar (a purple frame in Figure 13) is a small production field with three deep wells, but only one of them is productive. All three wells have temperature profiles showing inversion with depth. Hvíthólar is a small geothermal system just east of the presently active fissure swarm. It is probably isolated from the main system(s) to the north, as indicated by stable isotope ratios (see below). It is suggested that the geothermal system in Hvíthólar is near the end of its lifetime with two-phase conditions at a shallow depth but lower temperature below because the heat sources have run out of heat. Similar conditions are found in the Krýsuvík geothermal area on the Reykjanes peninsula, SW Iceland (ÍSOR database). The geological settings are similar, i.e., hyaloclastite ridges from the end of glaciation. It is suggested that at the end of the glaciation, the rapid pressure drop in the crust did initiate extensive intrusions and eruptions, which generated geothermal systems still present today, but fading out, because their heat sources have not been maintained.

Finally, well KS-1 in Sandabotnaskarð, east of Hvíthólar, encountered a deep geothermal system. The well is deviated to the east, under a fossil alteration on the surface, showing that, recently, there has been geothermal surface activity there. Shallow holes indicate lower temperatures between KS-1 and Hvíthólar (ÍSOR database). Both 1D and 3D inversion MT soundings show a deep low resistivity there (Figure 10). The system encountered in KS-01 seems to be isolated from the systems to the north and Hvíthólar because its fluids have a quite different isotope signature than the others (see below).

#### *5.2. Heat Sources, Magma Ascent and Volcanism*

Figure 16 shows a zoom-in of the residual Bouguer gravity in and around the inner caldera. The figure shows a local gravity high within the general gravity low in the inner caldera, mainly north of the transverse structure, but also across it, from Leirbotnar and towards a gravity high at Hvíthólar. These gravity highs show the same general NNE-SSW trends as the fissure swarm. The gravity highs north of the transverse structure extend from Mt. Krafla in the east and to Leirhnjúkur in the west. They are confined within the area where high-temperature alteration is found north of the transverse structure (purple lines in Figure 16, from Figure 9). It is, therefore, natural to assume that they reflect intrusions which are (and have been) heat sources for the main geothermal system within the norther part of the inner caldera.

**Figure 16.** A zoom-in on residual Bouguer gravity in and around the inner caldera. Blue lines show faults and fissures, green lines mark the boundaries of the part of the fissure swarm that was active in the Krafla Fires and an ESE-WNW fault in the southern slopes of Mt. Krafla that moved in the Krafla Fires. Yellow arcs and lines show explosion craters and eruptive fissures. Note, the abundant explosion craters west and south of Mt. Krafla. Grey lines mark the transverse structure and purple lines mark the extent of high-temperature alteration at 200 m a.sl. (from Figure 9). Coordinates are UTM WGS84 zone 28 in km.

As mentioned above, there is a local gravity high through the transverse structure from Leirbotnar to Hvíthólar. It roughly coincides with an eruptive fissure from Daleldar, about 1.1 ka old [7]. Resistivity does, however, not show shallow high-temperature alteration there. This gravity high could be due to shallow dikes that have not managed to generate a high-temperature geothermal system. The high gravity crossing the transverse structure connects to a clear gravity high south of Hvíthólar. This gravity high is a part of a more regional gravity high SE of the 110 ka caldera (Figure 7). The gravity map in Figure 16 does not show any gravity anomalies that could be associated to heat sources in Sandabotnaskarð and in the SW part of the inner caldera.

The gravity highs north of the transverse structure reflect the heat sources of the main geothermal system there. These are probably dikes. The opening component discussed above (Figure 5) makes a pathway for basaltic magma in the deeper eastern part of the calderas from the mantle. At depth, these intrusions/dikes probably have an ESE-WNW orientation, but at shallower levels, some of the magma probably migrates into the faults and fissures between Mt. Krafla and Leirhnjúkur, resulting in the NNE-SSW pattern of the gravity highs seen in Figure 16. For some reason, dike injections to the south into the transverse structure seems limited except for the gravity high extending towards Hvíthólar, discussed above. The ESE-WNW dike complex seems to be delineated in the resistivity model from a 3D inversion of MT data (Figure 10). It appears as a high-resistivity zone bordered by low-resistivity, both to the south and north. The resistive central part of the dike complex probably reflects fully crystalized rocks, while the conductive margins might be a partial melt and hence explaining the S-wave shadows (see discussion below).

The main pathway of ascending basaltic magma from the mantle is probably offset to the east of the presently active spreading as manifested by the abundant explosive craters in the western slope and south of Mt. Krafla (Figure 16). In between rifting events, basaltic intrusions may enter east of the fissure swarm without reaching surface but triggering explosions/eruptions (Figure 16) ejecting silicic tephra (see discussion below). In major rifting events with long dike injections and fissure eruptions like the Mývatn and Krafla Fires, the magma flows up and to the west towards Leirhnjúkur in the presently active fissure swarm which, on surface, looked as the centre of activity. This idea gets support from the fact that the Mývatn Fires started with an explosive eruption in Víti, just west of Mt. Krafla and about 2 km east of Leirhnjúkur and the fact that at the beginning of the Krafla Fires, fumaroles SE of Víti had large increase in volcanic gases, even larger than fumaroles in Leirhnjúkur [11]. It also gets support from the fact that an ESE-WNW trending fault in the southern part of Mt. Krafla moved in the Krafla Fires (Figures 4 and 16). As the magma ascends, it starts to degas, causing the severe gas contamination in the deeper two-phase part of the geothermal system west of Mt. Krafla. As mentioned earlier, it has been suggested [18] that there are several magma reservoirs in the volcano, based on detailed crustal deformation during an eruption in Krafla Fires. These might be the magma plumbing system conveying magma from a deep up-flow in the east and towards the fissure swarm at Leirhnjúkur.

#### *5.3. Origin of the Geothermal Fluids*

The origin of the geothermal fluids in Krafla has been studied and debated for many years. The origin of geothermal fluids is normally estimated on the bases of the deviation of ratios of stable isotopes of oxygen, δ18O, and hydrogen, δD, from that of Standard Mean Oceanic Water (SMOW) [46] and comparison with local meteoric values [47]. There can be a significant oxygen isotope exchange due to water-rock interaction but δD is assumed to be better preserved and a better indicator of the original fluid.

The δD of local rainfall in the Krafla area is about −86% [48]. The δD of well discharge in the Vítismór-Leirbotnar area is about −87% but about −92% in the southern slopes of Mt. Krafla east of Hveragil and in Hvíthólar [47]. This is considered a significant difference. In the discharge from well KS-1 in Sandabotnaskarð, δD is even lower or about −115% (ÍSOR database), similar to wells in Námafjall geothermal system south of Krafla. In [47], it is concluded that the geothermal system under Vítismór-Leirbotnar is recharged by local meteoric water and state that the geothermal systems east of Hveragil and Hvíthólar seem to have "the same source of recharge" of water from higher elevation.

In a recent study [49], it is argued that the isotope composition of the discharged fluids from the high/excess enthalpy wells east of Hveragil is distorted from that of the actual reservoir fluid towards a lower δD due to phase segregation and that the reservoir fluids in the geothermal system(s) in the inner caldera, north of the transverse structure, are of the same origin, i.e., local rainfall. This supports the idea put forward above, that what looks like two different geothermal systems east and west of Hveragil, are actually one system, the only difference being different permeability in the upper part.

The low δD in Hvíthólar and in KS-01 in Sandabotnaskarð cannot be attributed to phase segregation like in the flanks of Mt. Krafla. The enthalpy there is about 2500 kJ/kg while the enthalpy in Hvíthólar is about 1200 kJ/kg [12] and about 1500 kJ/kg in KS-01 (ÍSOR database). The geothermal fluids in KS-01 seem to originate from the general groundwater flow from higher elevation in the south with a similar δD as the fluids in Námagjall and local groundwater south of Námafjall [47]. The geothermal fluids in Hvíthólar might be a mixture of local rainwater and groundwater from south. It, therefore, seems that there is no (or very little) hydrological connection across the transverse low-density structure between the geothermal system in the northern part of the inner caldera and Hvíthólar and the system encountered by KS-01 in the southern part. Furthermore, as mentioned above, geophysical data and

the results of temperature measurements in boreholes indicate that the geothermal activity in Hvíthólar and Sandabotnaskarð are separate systems.

#### **6. The Role of Melting**

Two wells in Krafla have encountered silicic melt, K-39 and IDDP-1. K-39, which was drilled directionally towards east under the southern slopes of Krafla (for location see Figure 14), hit melt close to the bottom at about 2600 m below the surface [50]. The stuck drill string, cut loose by blasting, contained fresh rhyolite glass (69.0 to 78.8 wt.% silica). The bottom section of the hole was plugged with cement. The glass was interpreted to be quenched rhyolite magma derived by partial melting of geothermally altered basalt in a contact zone of a gabbroic intrusion [50].

IDDP-1, which was designed to be drilled vertically to 4 to 5 km in depth and into super critical fluids (for location see Figures 11 and 14), surprisingly hit silicic melt at the depth of 2104 m [51]. After drilling about 30 m with total loss of circulation, the penetration rate suddenly increased and the torque on the drilling assembly increased [52]. After pulling the drill string somewhat up and lowering it again, the drill bit go<sup>t</sup> stuck at the depth of 2095 m, so the magma had been squeezed up into the borehole and a large amount of both brown and clear glass was recovered. The drilling was then halted and a cemented casing was installed to 1960 m depth. The well was cooled by injecting cold water into it for a month. After about seven months of heating up, discharge tests started. In the beginning of the flow-test, the well discharged dark water and steam (Figure 17). After the well had discharged for some time, it discharged superheated steam at a temperature of around 450 ◦C and a pressure of 140 bars at the well head. The maximum flow was around 45–50 kg/sec with an enthalpy of about 3200 kJ/kg giving about 200 MWt and 35 MWe [53].

**Figure 17.** The beginning of a flow test of IDDP-1 (photo G. Ó. Friðleifsson).

The fluid discharged from IDDP-1 was enriched in CO2 and H2S as well as chlorine (Cl), fluorine (F) and boron (B) relative to normal geothermal fluids in Krafla [54]. The black colour of the discharge in the beginning was due to corrosion of the casing by hydrochloric acid formed when the chlorine from the superheated steam dissolved in water droplets. Several other wells in Krafla had the same character when discharged, so-called "black death", and are thought to have hit superheated steam close to magma [55].

The results of the drilling and flow-test of IDDP-1 clearly showed that right above the molten magma there was a permeable layer with superheated steam. Since "black death" has been observed in several other wells, it seems that similar conditions are found to be widespread in the reservoir inside the inner caldera, north of the transverse structure. This is in agreemen<sup>t</sup> with the observation and interpretation of the low P-wave velocity around 2.5 km depth within the ESE-WNW dike complex north of the transverse structure [34].

Detailed analysis of the glass from IDDP-1 [52,56] showed that it was quenched, near-liquidus rhyolitic magma (~76.5 wt.% SiO2) but with oxygen and hydrogen isotope values (δ<sup>18</sup> O = 3.1%; δ2H = −118%), depleted relative to mantle values. From this it was concluded that its composition is consistent with the formation by partial melting of basalt which had been hydrothermally altered by meteoric water. They estimate the in-situ temperature to be in the range of 850 to 920 ◦C and that it had degassed at a pressure of 40 MPa, which is higher than hydrostatic pressure but lower than lithostatic pressure expected at 2100 m depth. Based on crystal zoning, they further infer that the magma originated at a greater depth and migrated upward to its present position.

In a study of silicic rocks in Krafla and in the active volcanic zones in Iceland [13,57], it was concluded, based on chemical and oxygen isotopic composition, that they are derived from partial melting of older basaltic rocks that were hydrated by meteoric water. It was further concluded that the silicic magma is formed in intrusive complexes beneath the volcanoes rather than by di fferentiation in large long-lived magma chambers.

Krafla is not the only geothermal area where drilling into upper crustal magma has taken place. The first reported incidence was in 2005 in the Puna geothermal field, Big Island Hawaii [58], where a drilling hit dacite magma (67 wt.% SiO2) at the depth of 2488 m. Furthermore, in at least three wells drilled in Menengai caldera in Kenya, trachitic magma (63 wt.% SiO2) at the depth of 2040 to 2180 m was encountered [59]. A detailed petrological study has, so far, not been carried out on the quenched magma from Menengai.

All this suggests that rather shallow partial melting of old, hydrothermally altered volcanic rocks can be an important factor in volcanism and volcanic geothermal systems, at least in extensional tectonic environments. In evolved volcanos, the hydrothermal systems alter the host rocks, forming secondary hydrated minerals that lower the solidus temperature and produce more silicic early melt than the parent un-hydrated rocks. As the volcano builds up, the altered rocks subside to greater depths, particularly, in cases where a caldera structure develops. If a 1200 ◦C basaltic magma pools in this substrate, dehydration melting takes place, with partial melt segregating and buoyantly ascending. The melt is near liquidus and relatively rich in volatiles and hence with a low viscosity. As the melt migrates up, it cools, degases and the viscosity increases. Finally, it can be stagnated by minor structural obstructions, spreading out as sills. The occurrences discussed above indicate that for some reason, takes place at about 2–3 km depth.

If this is correct, re-melting of geothermally altered rocks and the ascent of the melt to shallower crustal levels can be a hitherto unrecognised but important heat transport mechanism from deep intrusions to the base of the geothermal system above. Therefore, as the volcano evolves, the geothermal system with the ascent of magma from below become a feedback system that facilitates heat flow from the mantle to the surface.

Geothermal drilling has shown that silicic intrusions are common in volcanic geothermal systems in Iceland. An exception is the Reykjanes peninsula, SW Iceland. In fact, no silicic rocks are found on surface west of the Hengill volcano SW-Iceland (this will be addressed below). In Krafla, silicic intrusions are more abundant in the eastern part of the wellfield, north of the transverse structure, where the rhyolitic explosion craters are found, than in the western part [12]. This supports the idea that basaltic intrusions at depth are more common in the eastern part of the inner caldera. Due to

the nature of their formation, the silicic intrusions are probably in isolated pockets and, therefore, not detected/resolved by the MT.

IDDP-1 showed that there is a permeable zone with superheated steam above the magma pocket. The steam is probably trapped in ductile rocks. The low P-wave velocity anomaly at depth [34] (Figure 12 right), probably reflects the extent of superheated steam above recent silicic intrusions or magma pockets and could therefore be a proxy of their distribution.

Partial re-melting of hydrothermally altered basalt can explain the bimodal behaviour of evolved volcanos such as Krafla and Askja, i.e., occasionally having large and often phreatic silicic eruptions but erupting only basalt in-between. Over time, silicic intrusions and magma pockets accumulate at shallow crustal levels. When a substantial amount of silicic magma has accumulated, major basaltic intrusion(s) can mobilise the silicic magma triggering an eruption. After they have been erupted, the cycle starts again.

As mentioned above, no silicic rocks are found on the Reykjanes peninsula west of Hengill. The only occurrences are a basalt andesite intrusion encountered at a depth of around 800 m in well SG-12 in Svartsengi [60] and some felsite veins at the depth of 4634.5 m in IDDP-2 at the tip of the peninsula [61]. This is probably because the Reykjanes peninsula is a thin and young, basically oceanic, crust and di fferent from crust further inland. It lacks the central volcanoes that have produced massive lava accumulations and caldera structures, where basalts are buried, hydrothermally altered, and subside to where they can melt.

The results of country-wide MT measurements show a deep, low-resistivity layer under most of Iceland [62], but this layer is absent under the Reykjanes peninsula. The nature of this low-resistivity layer has been a matter of discussion for a long time. It has been proposed [62] that it is partially molten rocks, but results of seismic studies, e.g., [63], show that this cannot be the case. An alternative explanation has been suggested [63] and supported by the author, namely that the deep conductive layer is due to dehydrating alteration minerals (chlorite, epidote and amphibole). Recent laboratory measurements [64,65] show an irreversible and substantial increase in electrical conductivity when rocks containing chlorite and epidote are heated above 600–700 ◦C under high pressure. This is explained by dehydration of the chlorite and epidote minerals, releasing hydrogen and water and resulting in the oxidation of Fe<sup>+</sup><sup>2</sup> to Fe<sup>+</sup><sup>3</sup> [66] and the minerals becoming electron/hole semi-conductors. If this is the correct explanation of the deep conducting layer under most of Iceland, it means that the crust is weathered/hydrothermally altered but the crust under the Reykjnes peninsula is not or at least to lesser extent. It, therefore, seems like the alteration of the crust somehow preconditions partial melting in central volcanos, but this is a subject for a separate study.

#### **7. The S-Wave Shadows**

One thing remains to be addressed. The S-wave shadows observed during the Krafla Fires have been interpreted as due to a volume containing magma. They do appear as a high resistivity in the 3D MT model, but bordered by low resistivity, both to the south and north (Figure 10). Recent measurement of resistivity of molten basalt and rhyolite from Krafla [67] give 1.2–0.8 Ωm for rhyolite melt at temperatures from 900 to 1000 ◦C and 1.7–0.5 Ωm for basalt melt at temperatures from 1170 to 1250 ◦C. These are really low values and extended volumes of magma, as indicated by the S-wave shadows, should, therefore, be observed by MT as low-resistivity bodies. The most natural explanation of the absence of such a body is that the melt is in unconnected pockets. This seems, however, in contradiction with seismological observations. Seismologists say that in order to have a substantial attenuation of S-waves, connected magma/liquid volumes have to have dimensions of the order of the wavelength of the S-waves, which, in this case, is about 0.5–1 km (e.g., [68,69]). An anonymous reviewer suggests, however, that the "S-wave shadows" could be a manifestation of the radiation pattern and network distribution or that a su fficiently dense distribution of melt pockets may be "seen" by the S-waves as a single viscous body. This needs to be studied by numerical modelling.

Here, the ESE-WNW trending high-resistivity body in the S-wave shadows is considered to be a dike complex. One possible scenario is that the low resistivity bordering it is still-molten dikes. The low resistivity south of the complex could be a dike injected during the Mývatn Fires when diking into the fissure swarm and eruptions were mainly to the south. The low resistivity on the northern side could be due to a dike injected at the northern margin of the dike complex in the Krafla Fires, where dike injections into the fissures swarm and eruptions were manly to the north. This could explain the S-wave shadows, because the S-waves from which it is inferred are mainly traveling perpendicular to such dikes. The low-resistivity bordering the dike complex could also be due to dehydration of alteration minerals as discussed above. However, this contradicting evidence from seismic and resistivity studies still remains un-resolved.

#### **8. Summary and Conclusions**

This study indicates that the geological structure of the Krafla volcano is more complex than hitherto believed. The 110 ka caldera hosts, now buried, an 80 ka inner caldera and both calderas are bisected by an ESE-WNW transverse low-density structure. Resistivity surveys show that geothermal activity has mainly been within the inner caldera, but it is cut through by the ESE-WNW structure. A di fference in the local crustal spreading directions south and north of the calderas leads to a N–S opening component, favouring the ascent of basaltic magma from the mantle and explaining why most of the volcanic activity is in the eastern part of the calderas.

The present thermo-hydraulic character of the main geothermal system can be understood by considering the geological and tectonic history. When crustal spreading moved back from the western part of the fissure swarm to the eastern part 3 ka ago, permeability above the deeper intrusions increased drastically, resulting in a vigorous almost isothermal convection in the part where faults and fractures where reactivated, but further to the east where fissures were not reactivated and two-phase conditions prevail there. Hydrogen isotope ratios in geothermal fluids in Krafla show that, north of the transverse structure, they are of a local meteoric origin, but the fluids in Sandabotnaskarð and Hvíthólar are from the general groundwater stream from higher altitude in the south. There seems to be little, if any, hydrological connection between the geothermal systems north and south of the low-density transverse structure.

The silicic magma encountered in K-39 and IDDP-1 indicates a hitherto overlooked heat transport mechanism in evolved volcanoes, i.e., an ascending silicic re-melt of altered basaltic rocks by basalt intrusions and getting stagnant magma pockets at shallow crustal levels, producing a superheated steam zone above them. Bimodal volcanic behaviour of evolved volcanos like Krafla and Askja, occasional silicic eruptions, often phreatic, but purely basaltic in-between, can be understood by considering the volcano and its geothermal system(s) as coupled systems, enhancing heat flow from the mantle. The geothermal system produces alteration for "distilling" out silicic magma. When a substantial amount of silicic intrusions/magma has accumulated, major basalt intrusion(s) may "ignite" them causing a silicic eruption and the cycle starts again.

An important conclusion from this work is that a holistic approach, considering di fferent and independent datasets and information, can shed light on complex structures of volcanoes and associated geothermal systems. Of particular importance is to study and take into account tectonics and geological history. Geophysical surveys and monitoring, sensitive to di fferent physical parameters, and geochemistry are also vital. These are inexpensive studies, but their derived models usually need to be checked by much more expensive drilling. Conceptual models, often based on limited data, should not ge<sup>t</sup> stagnant for decades like in the case of Krafla. Any new information should constantly be used to reconsider and update the conceptual models. This will make geothermal drilling and utilization more focused and cost-e ffective.

**Funding:** This research received no external funding.

**Acknowledgments:** The author wants to thank Thorbjörg Ágústdóttir and Halldór Ármannsson at ÍSOR and four anonymous reviewers for their valuable and constructive suggestions and corrections that greatly improved this paper.

**Conflicts of Interest:** The author declares no conflict of interest.
