ORIGINAL RESEARCH article

Front. Earth Sci., 26 March 2021

Sec. Volcanology

Volume 9 - 2021 | https://doi.org/10.3389/feart.2021.628967

Seismicity of the Northern Volcanic Zone of Iceland

  • Institute of Earth Sciences, University of Iceland, Reykjavík, Iceland

Abstract

A half century of monitoring of the Northern Volcanic Zone of Iceland, a branch of the North America—Eurasia plate boundary, shows that the seismicity is very unevenly distributed, both in time and space. The four central volcanoes at the boundary, Þeistareykir, Krafla, Fremrinámar, and Askja, show persistent but very low-level seismicity, spatially coinciding with their high-temperature geothermal systems. On their rift structures, on the other hand, seismicity is almost absent, except during rifting episodes. Krafla went through a rifting episode in 1975–1984 with inflation, interrupted by 20 diking events with extensive rifting, eruptive activity, and intense seismicity along an 80 km long section of the rift. During inflation periods, the seismicity was contained within the caldera of the volcano, reflecting the inflation level of the magma chamber. Diking events were marked by seismicity propagating away from the volcano into the fissure swarms to the south or north of the volcano, accompanied by rapid deflation of the caldera magma chamber. These events lasted from 1 day to 3 months, and the dike length varied between 1 and 60 km. The area around the Askja volcano is the only section of the Northern Volcanic Zone that shows persistent moderate seismicity. The largest events are located between fissure swarms of adjacent volcanic systems. Detailed relative locations of hypocenters reveal a system of vertical strike-slip faults, forming a conjugate system consistent with minimum principal stress in the direction of spreading across the plate boundary. A diking event into the lower crust was identified in the adjacent fissure swarm at Upptyppingar in 2007–2008. Four nests of anomalously deep earthquakes (10–34 km) have been identified in the Askja region, apparently associated with the movements of magma well below the brittle-ductile transition. Several processes have been pointed out as possible causes of earthquakes in the deformation zone around the plate boundary. These include inflation and deflation of central volcanoes, intrusion of propagating dikes, both laterally and vertically, strike-slip faulting on conjugate fault systems between overlapping fissure swarms, migration of magma in the lower, ductile crust, and geothermal heat mining.

Introduction

Volcanism in Iceland has its roots in two different processes. The mid-Atlantic plate boundary runs through the country, separating two of the largest lithospheric plates, the North America and the Eurasia plates and spreading at a rate of about 19 mm/year (e.g., ). At the same time the area is a hotspot fed by a mantle plume centered on the central highlands of Iceland but affecting an area much larger than the present extent of the country (e.g., Wolfe et al., 1997). The two processes interact and lead to a high rate of eruptions from about 30 different and active volcanic systems. The rate fluctuates within the range of 20–50 eruptions per century, with 22 eruptions in the last 50 years, for example (). Running a modern society in a highly active area requires substantial monitoring effort. In Iceland the effort is based to a large extent on a combination of two methods to track the pathways and storage of magma in the crust, seismic monitoring and surface deformation studies (; ; ). The first step in the detection strategy is to identify the volcanic systems that are preparing for an eruption. This includes detecting deviation in behavior from the long-term background seismicity. This is not a trivial task in the tectonic environment of an active plate boundary, where earthquakes are common (). Long time series are essential. Once the active systems are identified the next step is to install instrumentation to detect the pressure changes and magma movements within the systems. Fortunately, most of the 30 volcanic systems that have been active in the Holocene are in a passive state at present and do not require intensive attention. The last step in the monitoring and warning process requires the ability to detect and correctly identify the short-term precursors to the initiation of an eruption. This last phase of the precursory activity, when the magma breaks out of its storage volume and propagates toward the surface, is marked by intensified seismic activity and rapid crustal movements. All eruptions in the last few decades in Iceland have been preceded by such detectable short-term precursors, and in 2/3 of the cases they have led to a warning to the public before the eruption broke out ().

In the present paper, we intend to give an overview of the long-term seismic behavior of one branch of the volcanically active zones of Iceland, the Northern Volcanic Zone (Figure 1), based on the experience of the last several decades, or since sensitive seismic instrumentation was installed, mostly in the early seventies. This zone includes the volcanic systems of Þeistareykir, Krafla, Fremrinámar, and Askja, all of which show different seismic characteristics. Krafla went through an episode of magmatic rifting in the period 1975–1984 with many lateral dike intrusions, nine basaltic fissure eruptions, and extensive surface rifting. The Askja area has persistent, intermediate-level seismicity that has been ascribed to several different tectonic and magmatic processes in the crust. All the systems have low-level, persistent seismicity that correlates spatially with the surface manifestation of hydrothermal systems.

FIGURE 1

), usually named by the main central volcano: Þ Þeistareykir, K Krafla, F Fremrinámar, A Askja, and Kv Kverkfjöll. Their respective fissure swarms (or rift zones) are shown in yellow. Epicenters of earthquakes for the period 1995–2020, shown with black and red dots, are from the Icelandic Meteorological Office. Black dots are epicenters of the period 1995–2010 and M ≥ 2.0, and red dots are events of 2011–2020 with M ≥ 1.5. Shown are the Northern, Eastern, and Western Volcanic Zones (NVZ, EVZ, and WVZ), the Reykjanes Peninsula Rift (RPR), the South Iceland Seismic Zone (SISZ), the Grímsey Oblique Rift (GOR), and the Hreppar Microplate (Hr). Digital elevation and bathymetric models from the National Land Survey of Iceland and the Hydrographic Survey of the Icelandic Coast Guard. Seismic stations are shown with triangles, red for the Krafla period 1974–1994 and black for the digital stations 1995–2020.

Volcanic Zones and Volcanic Systems of Iceland

A good part of the volcanism in Iceland takes place at or near the mid-Atlantic plate boundary. As it crosses the Iceland platform, the plate boundary breaks up into a series of segments, rifts, transforms and oblique segments. In some places the boundary branches out and small microplates or tectonic blocks are formed that may move independently of the large plates. The spreading across the plate boundary in South Iceland is divided between two sub-parallel branches, the Western and Eastern Volcanic Zones. The two zones are connected by a transform zone, the South Iceland Seismic Zone. Together with the Central Iceland Volcanic Zone these branches demarcate a microplate, the Hreppar Microplate (), a diamond-shaped, 100 × 100 km crustal block that appears to move independently of the major plates. The 19 mm/year total spreading rate is partitioned unevenly among the two divergent branches (; ), about one-fifth taken up by the western branch and the rest by the eastern branch. The plate divergence in the northern part of Iceland takes place along a single boundary, the Northern Volcanic Zone (e.g., ). Oblique segments also exist and contain volcanism and strike-slip tectonism in close proximity to each other. The Reykjanes Peninsula Rift (RPR) and the offshore Grímsey Oblique Rift (GOR) are examples of this type. The structure of the divergent branches of the boundary is characterized by rifting. Normal faults, fissures, and eruptive fissures are abundant and are arranged into swarms that form volcanic systems together with one or two central volcanoes, where the volcanic activity is more concentrated than elsewhere (, ).

In addition to the plate boundary volcanism there are three zones of intraplate volcanism, or flank zones. The Snaefellsnes Volcanic Zone, the South Iceland group of volcanoes, and the Öraefajökull-Snaefell zone located east of the Eastern Volcanic Zone.

Note on terminology: The term “fissure swarm” in Iceland is used for an array of fissures and normal faults. Similar structure is called “rift zone” in other areas, e.g., in Hawaii. A “rift zone” in Iceland consists of several volcanic systems, i.e., several subparallel fissure swarms.

The Northern Volcanic Zone

This branch of the plate boundary separates the two major plates, the North America and the Eurasia plates. The vector of plate separation is slightly oblique to the trend of the zone, so the volcanic systems are arranged en echelon. Most of the swarms extend along a large portion of the zone, so crossing the boundary one may have to cross a 50 km wide deformation zone with four or five semi-parallel fissure swarms. This means that the total plate velocity is partitioned between several fissure swarms and is not taken up by only one. The principal volcanic systems are Þeistareykir, Krafla, Fremrinámar, and Askja, all accompanied by impressive fissure swarms (; ). Sometimes two more are added to the list, Heiðarsporður and Hrúthálsar, both with traces of evolved rocks but limited fissure swarms ().

The seismicity of the Icelandic rift zones is highly clustered, both in time and space (). Most of the seismicity occurs at the central volcanoes, whereas the fissure swarms are generally very quiet. And the seismicity of the central volcanoes is only moderate, except during major volcano-tectonic events or episodes. Thus, the background seismicity of the Northern Volcanic Zone for the last few decades has been characterized by weak seismicity of the main central volcanoes and moderate seismicity of an area east of the Askja volcano, 15 × 20 km in dimensions. The main fissure swarms of the rift, with their impressive normal faults and fissures, are seismically virtually inactive, however, even at the microearthquake level. The weak seismicity of the central volcanoes of Þeistareykir, Fremrinámar, and part of the Askja activity, appears to be spatially linked with the high-temperature geothermal areas of these volcanoes. Geothermal heat mining and cooling contraction, or other geothermal processes, may be responsible for this seismicity (). It is also noteworthy that silicic rocks are associated with these sites and may play a part in the causal process.

The most significant seismic activity of the Northern Volcanic Zone in the last decades took place in two areas: The Krafla volcanic system during the major magmato-tectonic events of 1975–1984 and in the Askja area persistently during the whole period of observations. These events will be discussed in Sections Krafla and Askja, respectively.

The monitoring capacity for earthquakes in the Northern Volcanic Zone has improved greatly during the last half century. An overview of the history of the seismograph system is given by . The first network in North Iceland with the capability to distinguish between the activities of different volcanoes was installed in 1974, just in time to detect the beginning of the Krafla rifting episode (). The network had 6–11 stations within 80 km distance from Krafla (Figure 1), and earthquakes down to magnitude 2 could be located with uncertainties less than 1–2 km (). The network was augmented by portable seismographs during critical times if conditions allowed, which was not very often. The studies of S-wave attenuation () and of the swarms of September 1977 and July 1980 (; ) were done with a network with a few extra stations added to the permanent network. The hypocentral locations were calculated by inversion programs like HYPOELLIPSE by and subsequent inversion programs developed at the USGS in Menlo Park, CA. Crustal velocity structure was based on the extensive work of and calibrated with explosions in a crater lake (). The accuracy of depth determinations was uneven within the network. Close to the seismic stations the uncertainty was of the order of 1 km, but farther away it was not so good. Therefore, the depth determination of activity close to the central volcano and the southern fissure swarm was good, but in a large part of the uninhabited northern fissure swarm the depth is not well constrained.

The first epicentral maps for the whole zone reveal the clustered seismic activity, with most of the background seismicity occurring in the central volcanoes (). It also showed the unexplained persistent seismic activity in the volcanic zone east of the Askja volcano.

After 1990, the seismic network in the country was gradually replaced by a digital, semi-automated network () (Figure 1), and the last analog seismograph station was decommissioned in 2010. The Askja seismicity attracted attention, and in 2005 and following years, the permanent network was augmented with temporary stations around Askja, by initiative from Cambridge University. The level of detection was lowered by up to two magnitude units (), with a stream of research papers on this remarkable cluster of activity within the Northern Volcanic Zone (e.g., ; White et al., 2011; , ).

Krafla

The Krafla volcanic system consists of a central volcano with a 10 × 8 km caldera of about 110 ka age, partly buried by later lavas (; ). The caldera is transected by fissure swarms that extend about 45 km to the south and 60 km to the north of the volcano, subparallel to the plate boundary. The only known previous historic rifting episode occurred in 1724–1746, with extensive rifting, apparently mainly in the southern fissure swarm, and several fissure eruptions in the caldera and southern fissure swarm (; ). The big lessons in volcano seismicity were provided by the Krafla magmatic episode that occurred in 1974–1989 (; Tryggvason, 1984; ; ). During this time interval about twenty inflation-deflation cycles of the Krafla central volcano were detected (Figure 2 and Table 1). During inflation of the central volcano the seismicity was limited to the caldera (Figure 3) and generally increased toward the end of the inflation phase. In several cases it can be demonstrated to follow the Kaiser effect (), i.e., the seismicity was at a minimum until the inflation level exceeded the maximum level of the previous inflation. Then, it increased abruptly. An important corollary is that inflation and deflation of a caldera can occur without significant or detectable change in seismicity. This range of aseismic inflation-deflation was of the order of one meter vertical displacement at the center. A high proportion of the caldera earthquakes showed anomalously small S-waves, indicative of S-wave shadowing by a magma body (), later confirmed by seismic profiling and seismic imaging of the caldera (; ). Most of the caldera earthquakes originated in the uppermost 3 km of the crust, below which an aseismic layer extended to a depth of 7 km. Several earthquakes at 7 km depth delivered S-waves to regional seismic stations, indicating a bottom to the attenuating body ().

FIGURE 2

, and Wright et al., 2012). The graph is based on data from leveling, interpolated with the help of a water level tilt meter record. The activity ended with the Krafla volcano in an inflated state.

TABLE 1

YMDDurDirLSpeedDefl.Erupt.Prec.tVol. lava
(km)(m/s)cmMm3
1975122090 dS, N60203Y15 mSmall
197609295 dN1214
197610312 dN2457
197701202 dN1029
197704272 dS1777Y(1 h)Small
197709081 dS130.3521Y2 1/4 hSmall
197711021 dN<3<4
1978010620 dN450.2105
197807103 dN250.557
197811105 dN1864
197905135 dN2484
197912026 dS?4?3
1980021010 dS1010
198003161 dN, S200.938Y65 mSmall
198007106 dN940Y5 h30
198010185 dN730Y82 m40
198012244 dN1010
198101305 dN627Y7 h32
198111185 dN81.850Y76 m40
1984090414 dN8.50.655Y264 m80

Dikes and deflation events of the Krafla episode.

Characteristics of the 20 deflation/diking events of Krafla 1975–1984: Date of the beginning, duration of the event from the beginning of deflation until inflation sets in again, direction of propagation into the N or S fissure swarms, length of the dike measured from the center of the caldera to the distal end of the epicenter distribution (), maximum speed of propagation, vertical subsidence at the center of deflation (), eruption, and precursor time from the beginning of deflation until an eruption begins ().

FIGURE 3

). Seismic stations are shown with red triangles in (A). Digital elevation and bathymetric models from the National Land Survey of Iceland and the Hydrographic Survey of the Icelandic Coast Guard.

During the deflation part of the cycle, the fissure swarm was activated (Figure 3) with propagating earthquakes and large-scale rifting (; ).

The activity can be divided into phases, each of which was accompanied by characteristic seismicity.

The Preparatory Phase, 1974–1975

A newly installed seismograph system began picking up events originating within the caldera of Krafla in the winter of 1974–1975. This was unexpected. No significant activity had been known there before, and the increase was confirmed by local residents who had never before felt local earthquakes. This increase was limited to the caldera. The issue was politically very sensitive because a new and controversial power station was under construction within the caldera.

The Initial Rifting Event, December 1975–March 1976

A new phase of the activity began suddenly on December 20 when a small basaltic fissure eruption broke out in the center of the caldera, preceded by only 15 min of intense earthquake activity () that was detected by seismic stations throughout Iceland. The earthquake activity soon propagated out of the caldera, first to the south and then to the north. The eruption then came to an end, and only lasted about 2 h. The activity south of the caldera lasted a few hours, but the earthquakes north of the caldera kept propagating for several days until they reached the intersection of the fissure swarm with the Grímsey Oblique Rift offshore. All together about 60 km section of the fissure swarms was activated in this event. The earthquakes continued, however, and so did the deflation of the caldera. Large-scale rifting was observed along the northern end of the fissure swarm.

On January 13 a magnitude 6.4 earthquake occurred offshore at the intersection with the Grímsey Oblique Rift. The focal mechanism indicated strike-slip faulting (, ). The earthquake caused considerable damage in the surrounding area, particularly in the Kópasker village. The earthquake was clearly triggered by the widening of the Krafla fissure swarm that was still going on immediately to the south (e.g., ).

The initial rifting event ended in March, both the deflation of the caldera and the earthquake swarm and rifting at the northern end of the fissure swarm. Soon after that the caldera began inflating.

Intrusive Activity North of the Caldera, September 1976–January 1977

The re-inflation of the caldera was an indication that the activity at Krafla had not come to an end. This was confirmed when earthquake activity in the caldera began increasing in the summer of 1976 (; ). The number of detected earthquakes in the caldera increased from a few events per day to more than a hundred. This activity suddenly stopped in September, when the caldera began to slowly subside. A flurry of earthquakes was detected in the fissure swarm north of the caldera indicating an intrusion in that area. This event only lasted five days and the caldera started inflating again. Much more dramatic intrusive activity followed, in October 1976 and January 1977, after the inflation level in the caldera had reached and exceeded the previous maxima, reached before the September deflation. This activity was accompanied by intense earthquake swarms and surface rifting in the northern fissure swarm.

Intrusive Activity in the Southern Fissure Swarm, April and September 1977

The caldera seismicity was very low after the deflation of January 1977, only a few detected events per day. It remained low despite re-inflation of the caldera. When the inflation level exceeded the previous maximum, however, the seismicity increased rather abruptly. The inflation continued for about two months until finally the magma storage was ruptured, and a dike began propagating from the caldera into the southern fissure swarm. A small eruption began at the northern caldera rim but stopped as soon as the dike began propagating. Extensive surface rifting occurred inside the caldera as well as in the southern fissure swarm (Figure 4).

FIGURE 4

). Faults and fissures are from . Red lines are eruptive fissures. Digital elevation model from the National Land Survey of Iceland.

Inflation of the caldera set in as soon as the April dike became inactive. This inflation (about 1 m, Figure 2) was aseismic and continued until the level of the previous inflation was almost reached. Then, a new diking event began rather abruptly with a small eruption at the northern caldera rim, in the same place as in the previous diking event. The eruption lasted about 2 h. Also this time, the dike propagated southwards, and the eruption came to an end (). A unique event took place when the dike had propagated through the Bjarnarflag geothermal area about 10 km south of the caldera. One of the drillholes used for power production erupted, punctured a steam pipe, and erupted a few tons of basaltic pumice for a while (). The dike did not propagate as far as the one in April (Figure 4), but extensive surface rifting was observed in the fissure swarm. The diatomite factory at Bjarnarflag was seriously damaged by the surface ruptures. The hypocentral depth of the earthquakes in the swarms of April and September was also consistent with a shallow intrusion. They occurred down to 3–4 km.

The dike intrusions of 1977 into the southern fissure swarm were accompanied by surface rifting and crustal extension within the caldera. The subsequent inflation of the caldera did therefore not lead to earthquake activity at the same inflation level as before. The correlation between the inflation and earthquake activity in the caldera became weaker ().

Intrusions in the Northern Fissure Swarm, 1978–1979

The inflation following the southern intrusions of 1977 by far exceeded the previous maximum level and was interrupted only by a small deflation event in November. In the beginning of the year 1978 the rate of inflation had decreased, and the caldera seismicity was low. A debate about the possible end of the activity was quickly silenced on January 6th when a new intrusion event began with a slow deflation at the caldera. The rate quickly increased, and earthquakes propagated into the northern fissure swarm. This event turned into the second largest dike intrusion of the Krafla rifting episode, in terms of volume, distance, and duration (Table 1). No eruption accompanied this intrusion but rifting occurred near the distal end of the dike, where it cut the main road. A re-measurement of a leveling profile there shows that the rifting has the form of a graben of 7 km width that subsided 1.2 m. The flanks were uplifted by 0.5 m (; ).

Similar events took place in July and November 1978, and May 1979 (). They were all accompanied by similar amount of deflation, 60–70 cm. The length of the dike became progressively shorter, however, and the surface rifting took place in different parts of the northern fissure swarm (Figure 4). No eruptions occurred but one drill hole within the caldera got clogged during the May 1979 intrusion. Repair and clearing of the hole revealed fresh basaltic glass.

Intrusions Into the Southern Fissure Swarm, February and March 1980

Inflation during the next seven months was interrupted by a small episode of slow deflation in the beginning of December 1979. A few small earthquakes at the southern caldera rim indicated that magma may have migrated a short distance to the south during this deflation. Inflation resumed and continued until February 10 when a new deflation event began with a swarm of small earthquakes south of the southern caldera rim (Figure 4). The activity slowly migrated southwards with varying intensity. The hypocenters of these earthquakes were unusually deep, between 5 and 7 km, deeper than during the intrusions in the same area in 1977. Very little surface rifting was observed. Apparently, this dike was intruding below the dikes formed in 1977. After ten days of deflation the caldera began inflating again.

The elevation of the caldera just about reached the previous maximum when a new deflation event began rather abruptly on March 16. This time the course of events was very rapid. The rate of deflation was unusually high and vigorous spasmodic tremor was recorded on the seismograph network with increasing amplitude, accompanied by earthquakes. An eruption broke out on a fissure within the caldera 65 min after the beginning of the deflation. The erupting fissure propagated toward the north, and in 25 min, it had attained a length of 4.5 km. The opening of the fissure was not accompanied by any noticeable increase in seismic activity. The earthquakes kept propagating to the north, however, at least for 5 more kilometers. Then, the events took an unexpected course. A new dike began propagating southwards from the caldera with larger earthquakes than before, exceeding magnitude 3 (Figure 4). This dike followed the same path as the February dike, but at shallower level. The eruption vigor decreased when the new dike began propagating and came to an end. The volume of this eruption was small, but still considerably larger than the three previous eruptions of the episode. It marked the end of the primarily intrusive part of the Krafla magmato-tectonic episode.

Eruption Phase, July and October 1980, January and November 1981, and September 1984

The remaining deflation events at Krafla were mostly eruptive, i.e., most of the mobilized magma found its way to the surface to feed eruption on long fissures within the caldera and adjacent northern fissure swarm. The only exception was a small deflation event in December 1980 when magma apparently was injected into the northern fissure swarm. The five eruptions shared several characteristics, even though the time scale and rate of change were variable. The beginning of the deflation was indistinguishable from what was observed in previous intrusive events. It was not until the dike approached the surface that a difference in seismicity was noticeable. Often low-frequency earthquakes were the last signal that an eruption was about to break out (). Usually, the earthquake activity decreased markedly when the eruption began and remained low during the eruption. The eruptive fissures of two of the eruptions (July 1980 and January 1981) were 4 and 2 km long, respectively, and located 5–7 km north of the caldera. The other three eruptions took place on much longer fissures, extending from the center of the caldera and 8–9 km to the north. The eruptions were most vigorous during the first day and then declined quickly. The activity concentrated on one or very few craters and faded out in 5–8 days. The last eruption (September 4–18, 1984) was an exception. The first few days it behaved like the previous eruptions, but after the activity had been limited to one crater the vigor of the eruption began to increase. This increase continued until the 14th day of the eruption when it suddenly ended.

Inflation, 1984–1989

The Krafla volcano began inflating again when the eruption of September 1984 was over. The inflation continued intermittently for several years, until 1989. The inflation was accompanied by a persistent background of small earthquakes, rarely exceeding magnitude 2, showing heterogeneous stress, with horizontal shear, extensional and compressional sources in close juxtaposition. A part of these events had sources with a large non-double couple component (). By 1989, the inflation had exceeded the level reached before the 1984 eruption (). Since then, a slow deflation has been the main characteristic of the crustal deformation of the volcano, most likely reflecting gradual cooling of the caldera crustal volume and harnessing of the geothermal system (). The seismicity of Krafla during this phase of slow deflation has been low, similar to that of other central volcanoes of the zone with a geothermal system, Þeistareykir, Fremrinámar and Askja, with a couple of events per year of magnitude 2 and larger.

General Pattern

The course of events was governed by magma chamber inflation that continued for at least 15 years (Figure 2), interrupted by about 20 short deflation events when magma escaped from the holding chamber in the middle of the caldera and formed a dike that propagated along the fissure swarms, either to the north, the south, or both (Table 1) (; Wright et al., 2012). These events were of variable magnitude, even though the deformation pattern in the caldera was the same. The shape of the deformation indicates a depth to the pressure source of about 3 km (; Tryggvason, 1980). Thus, since the shape of the deformation field was unchanged, the relative size of the deflation events may be quantified by the total vertical displacement in the center of deflation that coincided with the center of the caldera and scales linearly with the volume of the deflation. The central maximum subsidence ranged between 3 and 203 cm (Table 1). Similarly, the final length of the dike that was formed in the fissure swarms ranged between a few km and 60 km (Figure 2). The seismicity changed dramatically during these deflation events. The caldera seismicity stopped instantaneously when the deflation began (Figure 3). Then, a flurry of events began near the caldera, propagating outward along the fissure swarm (). The speed of propagation was high in the beginning but slowed down as the distance from the caldera increased, see Figure 5 (). The rate of deflation followed the same pattern, was high in the beginning and then slowing down. The magnitude of the fissure swarm earthquakes increased as the swarm progressed and formed a cluster in a particular part of the fissure swarm, where large scale surface rifting was also observed in the largest events. The maximum speed of propagation was quite variable between deflation events, ranging between 0.2 and 1.8 m/s (Table 1). The duration of the deflation events was also quite variable, between one day and three months.

FIGURE 5

The dikes reached the surface in nine of the cases, producing eruptions of basaltic lava from fissures of variable length, up to 9 km (). The first three eruptions were tiny, almost like geological accidents. The fourth eruption was slightly larger but of short duration. The volume of all these four eruptions was only a very small fraction of the volume of the dike or the deflation bowl. A new stage was reached with the eruption in July 1980. Then, most of the mobilized magma appears to have reached the surface. A good fraction of it flowed back into open fissures, however, and caused secondary rifting. The fissure that received the lava dilated and extended in both directions. Apparently, a dike was formed, fed from the surface, a phenomenon observed here for the first time. The flow into the fissures lasted at least 18 h, and for several of those hours, the whole lava production went this way. Subsequent dikes also led to eruptions. The largest one was the last event in this remarkable sequence. All the eruptions except the first one were anticipated or predicted based on the combination of deformation and seismicity, and the warnings were issued to the local population before them ().

During injection of dikes outside the caldera the deflation of the caldera was generally accompanied by very low caldera seismicity, sometimes with almost no detectable caldera earthquakes. Only in the two largest deflation events was there a detectable increase in the caldera seismicity. Earthquakes as large as mb 5 occurred during the first deflation event (), and a slight increase in small earthquakes was detected toward the end of the deflation event of January 1978. Deflation at the caldera center was 2 m and 1 m, respectively, in these events (Table 1). Deflation of less than 1 m apparently took place within the elastic range of the caldera floor. Dikes intruded within the caldera were clearly accompanied by seismicity during their propagation. Caldera earthquakes during inflation were significantly smaller than during deflation. They only rarely reached magnitude 4.

Attempts to estimate the volume of magma transported from depth to the shallow crust give varying results. used digital elevation models derived from aerial and spy photographs to estimate a total volume of intruded dikes 0.9–1.7 × 109 m3. They also estimate a volume of erupted lava of 0.2–0.4 × 109 m3. Values in the range 0.25–0.35 × 109 m3 were estimated from the surface flows (e.g., , ; ). Using the cumulative deflation of the caldera to estimate the total volume of dikes and flows gives 0.79–0.85 × 109 m3 ().

Askja

The Askja volcanic system consists of a central volcano and fissure swarms extending to the north and south. The northern swarm extends all the way to the north coast. The extent of the southern swarm is more uncertain. It interfingers with the swarms of the Bárðarbunga volcano and disappears beneath the Vatnajökull glacier. The structure of the Askja volcano is dominated by at least four overlapping caldera structures (; ). Seismic tomographic studies reveal bodies of anomalous VP/VS ratios underlying the calderas at discrete levels as deep as 20 km (). The most recent caldera was formed in 1874 and following years, after a major rifting episode with basaltic eruptions in the fissure swarm and a large silicic explosive eruption of the central volcano (e.g., ; ). The activity since then has been limited to the central volcano, most recently in 1961. Since about 1973, the main caldera has been slowly deflating. The cumulative subsidence of the caldera center is approaching 0.4 m (Tryggvason, 1989; ; ; ). The cause of the subsidence is generally thought to be a combination of a cooling and contracting magma chamber at a divergent plate boundary.

The Askja area in general is moderately seismically active (Figures 1, 6), and only a part of the activity appears to be related to the volcano itself (; ). Several processes have been suggested as sources of the different classes of earthquakes, including geothermal heat mining, transfer of spreading between overlapping fissure swarms, intrusion of a dike in the lower crust, and deep magma movements.

FIGURE 6

, seismogenic, conjugate faults east of Askja volcano are from , the blue ones are left-lateral, and the yellow ones are right-lateral. Locations of deep earthquake nests from are shown with orange patches. The Upptyppingar dipping dike is shown with the red-orange patch, from . The high earthquake activity in the SW-corner of the map accompanied the dike propagation and caldera collapse of the Bárðarbunga volcano in 2014–1015 (; ). Digital elevation model from the National Land Survey of Iceland.

Geothermal Heat Mining

A persistent cluster of shallow level hypocenters is spatially correlated to the high-temperature geothermal field in the SE corner of the main caldera, suggesting a causal relationship between geothermal heat mining and brittle failure. This is also an area of silicic rock outcrops, and in 2014, this part of the caldera fault collapsed in a major rockslide, causing a large tsunami in the caldera lake (). Small swarms of shallow earthquakes occasionally occur in other parts of the caldera and to the N and NW of it. Most likely these events are related to the gentle deflation of the caldera.

Transfer of Spreading Between the Askja and Kverkfjöll Fissure Swarms

A large majority of the seismicity in the Askja area occurs east of Askja and appears to have no direct connection to the volcano. This activity was first described by and . One NE-striking seismogenic fault was defined by the epicentral distribution, highly oblique to the fissure swarms nearby. Fault plane solutions showed strike-slip faulting (). More detailed work using a dense temporary network by and revealed more active faults in the area. Relative locations of numerous events show a system of conjugate faults bridging the gap between the fissure swarms of Askja and Kverkfjöll (), NE-striking faults with left-lateral strike slip on them, and N-striking faults of right-lateral faulting (Figure 6). None of these faults are visible at the surface. The conjugate faults suggest a causative stress field with the least compressive principal stress in a WNW-ESE direction, consistent with the direction of plate spreading. The depth of the hypocenters is less than 8 km, well within the brittle part of the crust. This activity therefore appears to be caused by transfer of spreading between the two overlapping fissure swarms.

Upptyppingar Lower Crustal Intrusion

An unusual sequence of small earthquakes was detected in the Northern Volcanic Zone, near the hyaloclastite mountain Upptyppingar east of the Askja volcano, beginning in February 2007 (Figure 6). The events were ordinary in appearance, but originated at unusually large depth, 10–20 km. The swarm was very persistent and spread out along a planar structure (; White et al., 2011; ) beginning at depth and then spreading upwards and eastwards. After a year the activity faded away and stopped. The swarm was accompanied by uplift of a large area. Modeling of the surface deformation revealed a planar intrusion with a strike of 80°E and dip of 45° to 49°coinciding with the planar structure outline by the hypocentral distribution (). The volume was (42–47) × 106 m3 ().

The intrusion is located within the Kverkfjöll fissure swarm (Figures 1, 6), about 40 km north of the Kverkfjöll central volcano. No deflation was detected of the Kverkfjöll central volcano in association with this intrusion, nor was the Askja volcano affected, the nearest neighboring central volcano. Apparently, the intrusion was fed directly from below. It is noteworthy that the orientation of the dike plane does not conform to any known structure at the surface. Fracture orientations are generally consistent with the stress field at a divergent plate boundary with the least compressive principal stress in the direction of relative plate motion (). suggested that the stress field implied by the dike intrusion was modified by the presently ongoing retreat of the ice caps in Iceland, leading to enhanced capture of magma within the crust.

Abnormally “Deep” Swarms

A special class of abnormally deep earthquakes was first identified and described by in the area within and around the Askja volcano. These events originate in the depth range 10–34 km, i.e., well below the brittle-ductile transition. They are typically small, M < 2, and occur in dense and short-lived swarms of a few minutes’ duration. The rapid succession of the events within a swarm may sometimes give the swarm the appearance of a tremor burst. The sources are quite persistent, on a time scale of years (). Four nests of deep events have been located in the Askja area, two in the area of the central volcano, one about 10 km east of the volcano and one about 15 km NNE of the volcano (Figure 6). suggested that these events represent magmatic injections along the fissure swarm of Askja, and suggested that they represent magma transport through the lower crust.

Discussion

Lateral Diking

The observations at Krafla of propagating earthquakes away from the central volcano and into the fissure swarms on either side of the volcano, accompanied by deflation of the caldera and extensive rifting in the fissure swarms, were interpreted as a consequence of dikes propagating laterally from a magma chamber within the caldera and into the rift zone (e.g., ; ; Tryggvason, 1980). The interpretation was inspired by similar events at Kilauea volcano in Hawaii (e.g., ), despite differences in tectonic setting. The horizontal extensional stress, providing the conditions for lateral propagation, was generated by plate divergence in Iceland but by gravitational sliding of the flank of Kilauea in Hawaii. There was no general consensus about this interpretation, however. Alternative interpretations were proposed (e.g., ). Cases of lateral diking accompanied by deflation of a magma storage or even collapse of a caldera have been documented since the Krafla activity began. These include the Piton de la Fournaise in 2007 (), Miyakejima in 2000 (), Asal rifting of 1978, Dabbahu in Afar 2005–2010 (Wright et al., 2012), Bárðarbunga in Iceland in 2014–2015 (; ), and Kilauea, Hawaii 2018 (). Lateral draining of magma reservoirs is an important mechanism to form collapse calderas, involving processes like external tectonic stress, magma density and buoyancy in the crust, magma viscosity, and viscoelasticity in the crust (; ).

Magmatic and Non-magmatic Rifting

One of the remarkable features of the Icelandic plate boundary rift zones is their highly variable seismicity, even on a time scale of decades. The rifting segment of the Eastern Volcanic Zone is, for example, almost totally aseismic, and so are segments of the Northern Volcanic Zone. Yet the relative plate divergence rate across these zones is 15–20 mm/year. Other parts of the divergent plate boundary, such as the segments of the Western Volcanic Zone, are moderately seismically active despite much lower spreading rate. , suggested that this discrepancy was the result of cyclic activity. The deviatoric stress in a rift zone immediately following a rifting episode is low and the deformation due to plate spreading occurs within the elastic limit of the crust. The seismicity will be very low. With increasing stress, the deformation will gradually lead to failure and the background seismicity will increase. If magma becomes available, the failure limit of the crust is lowered and rifting may ensue even if the background seismicity is still low (; ). This was the case before the Krafla rifting episode in the early seventies. If, on the other hand, magma is not available, the stress may increase to the dry failure limit of the crust and failure takes place by normal faulting without the help of magma, as seems to be the case in a segment of the WVZ. This may be the explanation for the different level of seismicity in different segments of the Northern Volcanic Zone, evident in Figure 1. The recent large rifting episode of the Askja system in 1874–1875 relieved stress in the central part of the zone and the stress has not yet built up to seismogenic levels.

Faulting Between Overlapping Rifting Structures

Rifting in the rift zones of Iceland occurs within fissure swarms, apparently mostly by diking during rifting episodes. Modern examples of this type of activity include the Krafla rifting episode of 1975–1984 and the Bárðarbunga dike of 2014 (; ). Each rift zone includes several fissure swarms that are arranged within the zone parallel to each other or en-echelon. Each dike or rifting event may only affect a limited length of one fissure swarm. This leads to accumulated horizontal shear stress between adjacent fissure swarms that is released in some kind of transfer motion on sets of conjugate strike-slip faulting. Two examples of such tectonism exist in the Icelandic rifts, the persistent seismicity east of the Askja volcano in the Northern Volcanic Zone, and the seismic activity in the Reykjanes Peninsula Oblique Rift (; ). In both areas, the orientation of the inferred least principal stress is horizontal and in the general direction of plate spreading.

“Deep” Crustal and Upper Mantle Earthquakes

Majority of earthquakes in Iceland originate in the upper, brittle part of the crust, in the depth range 0–10 km, and are apparently linked with relative plate movements across the plate boundaries. There are, however, occasional deeper earthquakes, originating in material that would, at normal strain rates, deform in a ductile way under the pressure and temperature conditions corresponding to the depth. Two classes of deep earthquakes may be identified: (1) earthquake sequences occurring immediately preceding or during eruptions and (2) small events occurring in short-lived swarms in persistently active “nests” beneath some parts of the volcanic zones.

Examples of the first class of deep events are provided by the Heimaey eruption of 1973 and the Eyjafjallajökull eruptions of 2010. Most earthquakes during the Heimaey eruption occurred at 15–25 km depth and most likely also the swarm of events immediately preceding the eruption (; ). The Eyjafjallajökull eruptions of 2010 were preceded by a series of shallow-level sill intrusions in 1994, 1999, 2009, and 2010 (, ; ), and a period of anomalously deep earthquakes in 1996 (; ). Some of the earthquakes during the eruptions also occurred at lower-crustal depths and were clearly related to the feeding channels of the eruptions (,b). The prolonged earthquake sequence during the Upptyppingar dipping dike intrusion of 2007–2008 may also be included in this class, although no eruption ensued.

The second class of deep events, first identified on the temporary Askja network by in the area around the Askja volcano, occur in short-lived swarms, of a few minutes’ duration. In addition to the four nests of deep events located in the Askja area, at least three additional nests of deep swarms have been identified and are regularly recorded by the permanent network. The most active one of them is located 15–20 km ESE of Bárðarbunga volcano, slightly south of the dike that propagated laterally from the volcano in August 2014 and fed the Holuhraun eruption (). The relationship of this nest to the dike and the plumbing of the Bárðarbunga volcanic system is still unclear (). Several swarms per month are detected in this nest at the present time and the activity is not decreasing. Another nest with a spatial relationship to the Bárðarbunga system is located beneath the Trölladyngja lava shield, within the NE-fissure swarm of Bárðarbunga. Only one deep nest is yet known in South Iceland. It is located beneath the eastern rim and flank of the Katla volcano (Vogfjörð et al., 2009).

The detection of the deep nests depends critically on the seismic network. A seismic station in the nearfield of the earthquakes is required for sufficient depth control on the location of the events. It is rather likely that more of these nests will be identified in the future. The known nests are all, except one, located in the central area of the Iceland Hotspot. Some parts of the plate boundary areas seem to be devoid of deep nests. None have been found in the Reykjanes Peninsula Oblique Rift or in the Krafla volcanic system despite relatively dense networks in operation there. The deep swarms in Iceland show great similarity with recently discovered events in the Eifel volcanic field in Germany (; ).

Stages Within Eruptive Cycle

A majority of the eruptions of Icelandic volcanoes within the last several decades occurred after a time of inflation of the respective volcano (; ). The eruption is then accompanied by deflation and followed by repeated inflation. This includes the nine Krafla eruptions of 1975–1984, the three eruptions of Grímsvötn 1998–2011, and the three eruptions of Hekla 1980–2000. In the cases of Krafla and Grímsvötn, the seismicity follows the cyclic behavior; it increases during the inflation part of the cycle and is greatly reduced when the deflation sets in. If a dike propagates away from the inflated area, its propagation can often be traced by the earthquakes that are generated around the dike tip. The re-inflation is accompanied by little seismicity during its first stage. It may increase rather suddenly if the inflation level exceeds the previous maximum, a phenomenon called Kaiser effect (). The deflation itself takes place almost aseismically as long as the strain remains within the elastic limits. In the case of Krafla, the elastic limit was exceeded only twice. The first deflation was about 2 m and was accompanied by considerable seismicity, with maximum magnitude 5.0. The second case was the January 1978 deflation (Table 1) when the subsidence in the caldera reached 1 m. There was a slight increase of the daily number of small caldera earthquakes during the last days of deflation. Based on these observations, it may be concluded that the elastic range of the Krafla caldera is about 1 m. Earthquakes associated with the inflation part of the cycle only reached magnitude 3.5.

Similar behavior is observed in Grímsvötn. The seismic activity of the volcano is generally very low but increases slowly in the latter half of the inflation period. The maximum magnitude does not exceed 3. The deflation has so far always been within the elastic limit.

The large Bárðarbunga 65 m deflation of 2014–2015 was accompanied by numerous earthquakes, more than 70 of which exceeded magnitude 5 (). The caldera seismicity decreased greatly when the deflation stopped in February 2015. The volcano began inflating again after the eruption and deflation ended but the seismicity remained low until September 2015 (). Then, both the earthquake frequency and the maximum magnitude increased and remained high for several years (). The maximum magnitude exceeded 4 but has remained below 5, however. All three volcanoes appear to have a similar seismic behavior with respect to inflation and deflation. Inflation-related earthquakes are generally smaller than the earthquakes accompanying deflation.

Geothermal Activity and Seismicity

It has been noted in many publications that there is a general correlation between persistent microearthquake activity within the volcanic zones of Iceland and the high-temperature geothermal areas (e.g., Ward and Björnsson, 1971). It may be argued that temperature changes associated with geothermal systems and the transfer of heat from the heat source to the colder bodies of rock leads to stress changes and thus to micro-cracking in the brittle part of the crust. The non-double-couple part of the earthquake source mechanisms identified by ; , and in the Hengill and Krafla areas may be taken as evidence for this. This mechanism has also been suggested as one of the possible sources of the seismicity of the Katla and Askja calderas, Torfajökull volcano, and the subglacial Loki Ridge (; , ; ; ). It is rather unlikely, however, that all the background seismicity of these locations is due to the geothermal mining of heat. The case of the Kerlingarfjöll area within the Hofsjökull volcanic system serves as a counterexample. This central volcano has a very powerful geothermal system, and yet, the seismicity is very low. This case suggests that the geothermal effect is likely to act as a trigger rather than the main driving force of the seismicity. The Kerlingarfjöll volcano is located outside the main deformation zone of the plate boundary. The driving force may therefore be missing even though the trigger is there.

Elevated pressure in the pore fluid may be another mechanism to release rock stress in a geothermal system close to a volcano by lowering the effective stress on pre-stressed faults. The validity of this mechanism is demonstrated in cases where the pore pressure change is induced by man, e.g., at the Hellisheiði Geothermal Power Station in SW-Iceland, where seismicity and crustal deformation is induced by injection of wastewater into the crust ().

We point out that many of the high-temperature geothermal systems are spatially correlated with outcrops of silicic rocks. This raises the question whether shallow, unstable, low-density bodies of silicic rocks (cryptodomes) may be the source of both the geothermal activity and the earthquakes.

Conclusion

Seismicity serves as the most important element in the long-term monitoring of the 30 + active volcanic systems of Iceland. About half of them have shown evidence of magma movements in the last 46 years of coverage by sensitive seismographs.

Processes identified as causes of earthquakes within the volcanic zones include volcano inflation and deflation, caldera collapse, dike and sill intrusion, amagmatic rifting by normal faulting, geothermal heat mining, strike-slip faulting on systems of conjugate faults, and unspecified movements of magma in the deep crust. Dome and cryptodome activity may also play a role.

Seismicity within the volcanically active areas of Iceland is very unevenly distributed, both in time and space. Persistent seismic sources are mostly associated with central volcanoes. The fissure swarms of the volcanic systems are usually inactive, even at the microearthquake level, except during rifting episodes when dikes propagate within them. The major rifting episode in the Krafla area 1975–1984 showed lateral injection of dikes, a model that was verified by the Bárðarbunga dike injection and rifting of 2014.

The long-term behavior of the different volcanic systems of the Northern Volcanic Zone is quite varied, ranging from very low-level background activity like Fremrinámar and Þeistareykir, to almost continuous, intermediate-level seismicity, like that of the area around the Askja volcano. For successful monitoring it is therefore essential to have a long time-history of the seismicity. A part of the seismic activity east of Askja is generated by strike-slip faulting on conjugate fault sets, taking up differential strain between overlapping rift zones of adjacent volcanic systems.

The rifting events at Krafla in 1975–1984 occurred during a phase of cyclic behavior, with periods of magma accumulation at depth and inflation punctuated by deflation, dike intrusions, rifting and eruptions. The cycle of inflation and deflation was reflected in seismicity that provides an additional parameter to apply in the monitoring effort. The examples of this process so far suggest that the largest earthquakes during deflation are significantly larger than the inflation-related earthquakes of the same volcano.

Persistent activity of small earthquakes in all four volcanic systems appears to be spatially correlated with high-temperature hydrothermal systems. The causal relationship is not always obvious. Cooling of hot rock bodies and consequent contraction and associated stress change may provide part of the explanation, also pore pressure fluctuations within the geothermal system. It is also pointed out that both phenomena, seismicity and geothermal activity, in some cases could have a common cause, such as a rising cryptodome.

Statements

Data availability statement

The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.

Author contributions

Both authors contributed equally to the subject of the study. PE with emphasis on the text and BB with emphasis on data analysis and construction of the figures. Both authors have worked for several decades on seismic monitoring and tectonics of Iceland and were actively involved in the seismic monitoring of the Krafla rifting episode.

Acknowledgments

This study is based on data and active cooperation of many people and institutions, including the Lamont-Doherty Earth Observatory and Cambridge University. The monitoring effort during the Krafla rifting episode relied to a large extent on the work of devoted local attendants of the analog seismograph system. The seismic data since 1990 are obtained from the Icelandic Meteorological Office. Seismic monitoring in the Northern Volcanic Zone has been partly supported by Landsvirkjun, the National Power Company of Iceland. The Generic Mapping Tools software (Wessel et al., 2019) was used to generate most of the Figures. Comments from the Editor and two reviewers improved the study considerably.

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

References

  • 1

    AcocellaV.TrippaneraD. (2016). How diking affects the tectonomagmatic evolution of slow spreading plate boundaries: overview and model.Geosphere12867–883. 10.1130/GES01271.1

  • 2

    ÁrnadóttirT. H.LundB.JiangW.GeirssonH.BjörnssonH.EinarssonP.et al (2009). Glacial rebound and plate spreading: results from the first countrywide GPS observations in Iceland.Geophys. J. Int.177691–716. 10.1111/j.1365-246X.2008.04059.x

  • 3

    BjörnssonA.EysteinssonH. (1998). Breytingar á landhæð við Kröflu 1974-1995. Samantekt á landhæðarmælingum (Changes in Elevation at Krafla 1974-1995. Compilation of Elevation Measurements, in ICELANDIC).Reykjavík: Nordic Volcanological Institute. National Energy Authority, Report OS-98002, Report 9801.

  • 4

    BjörnssonA.SæmundssonK.EinarssonP.TryggvasonE.GrönvoldK. (1977). Current rifting episode in north Iceland.Nature266318–323. 10.1038/266318a0

  • 5

    BjörnssonS.EinarssonP. (1974). “Seismicity of Iceland,” in Geodynamics of Iceland and the North Atlantic Area, ed.KristjánssonL. (Dordrecht: Reidel Publ. Co), 225–239. 10.1007/978-94-010-2271-2_16

  • 6

    BrandsdóttirB. (1992). Historical accounts of earthquakes associated with eruptive activity in the Askja volcanic system.Jökull421–12.

  • 7

    BrandsdóttirB.EinarssonP. (1979). Seismic activity associated with the September 1977 deflation of the Krafla central volcano in NE-Iceland.J. Volc. Geotherm. Res.6197–212. 10.1016/0377-0273(79)90001-5

  • 8

    BrandsdóttirB.MenkeW. H. (1992). Thin low-velocity zone within the Krafla caldera, NE-Iceland attributed to a small magma chamber.Geophys. Res. Lett.192381–2384. 10.1029/92gl02541

  • 9

    BuckW. R.EinarssonP.BrandsdóttirB. (2006). Tectonic stress and magma chamber size as controls on dike propagation: constraints from the 1975-1984 Krafla rifting episode.J. Geophys. Res.111:B124041.

  • 10

    DahmT.BrandsdóttirB. (1997). Moment tensors of microearthquakes from the Eyjafjallajökull volcano in South Iceland.Geophys. J. Int.130183–192. 10.1111/j.1365-246x.1997.tb00997.x

  • 11

    DahmT.StillerM.MechieJ.HeimannS.HenschM.WoithH.et al (2020). Seismological and geophysical signatures of the deep crustal magma systems of the Cenozoic volcanic fields beneath the Eifel, Germany.G321:e2020GC009062. 10.1029/2020GC009062

  • 12

    de Zeeuw-van DalfsenE.RymerH.SturkellE.PedersenR.HooperA.SigmundssonF.et al (2013). Geodetic data shed light on ongoing caldera subsidence at Askja, Iceland.Bull. Volcanol.75:709. 10.1007/s00445-013-0709-2

  • 13

    DrouinV.SigmundssonF.ÓfeigssonB. G.HreinsdóttirS.SturkellE.EinarssonP. (2017a). Deformation in the Northern Volcanic Zone of Iceland 2008-2014: an interplay of tectonic, magmatic, and glacial isostatic deformation.J. Geophys. Res. Solid Earth1223158–3178. 10.1002/2016JB013206

  • 14

    DrouinV.SigmundssonF.VerhagenS.ÓfeigssonB. G.SpaansK.HreinsdóttirS. (2017b). Deformation at Krafla and Bjarnarflag geothermal areas, Northern Volcanic Zone of Iceland, 1993-2015.J. Volcan. Geotherm. Res.34492–105. 10.1016/j.jvolgeores.2017.06.013

  • 15

    EinarssonP. (1978). S-wave shadows in the Krafla caldera in NE-Iceland, evidence for a magma chamber in the crust.Bull. Volc.411–9.

  • 16

    EinarssonP. (1986). “Seismicity along the eastern margin of the North American Plate,” in The Western North Atlantic Region: Geological Society of America, edsVogtP. R.TucholkeB. E. (Boulder, CO: The Geology of North America), 99–116. 10.1130/dnag-gna-m.99

  • 17

    EinarssonP. (1987). “Compilation of earthquake fault plane solutions in the North Atlantic and Arctic Oceans,” in Recent Plate Movements and Deformation, Geodynamics Series, Vol. 20ed.KasaharaK. (Washington, DC: America Geophysical Union), 47–62. 10.1029/gd020p0047

  • 18

    EinarssonP. (1991a). Earthquakes and present-day tectonism in Iceland.Tectonophysics189261–279. 10.1016/0040-1951(91)90501-i

  • 19

    EinarssonP. (1991b). “The Krafla rifting episode 1975-1989,” in Náttúra Mývatns, edsGarðarssonA.EinarssonA. (Reykjavík: Icelandic Nature Science Society), 97–139. 10.1016/0377-0273(84)90036-2

  • 20

    EinarssonP. (2008). Plate boundaries, rifts and transforms in Iceland.Jökull5835–58.

  • 21

    EinarssonP. (2018). Short-term seismic precursors to Icelandic eruptions 1973-2014.Front. Earth Sci. Volcanol. Res. Topic6:45. 10.3389/feart.2018.00045

  • 22

    EinarssonP.BrandsdóttirB. (1980). Seismological evidence for lateral magma intrusion during the July 1978 deflation of the Krafla volcano in NE-Iceland.J. Geophys.47160–165.

  • 23

    EinarssonP.JakobssonS. (2020). Scanning the analog seismogram archives of Iceland.Jökull70:16. 10.33799/jokull2020

  • 24

    EinarssonP.SaemundssonK. (1987). “Earthquake epicentres 1982-1985 and volcanic systems in Iceland: A map,” in Í hlutarins eðli, Festschrift for Þorbjörn Sigurgeirsson, ed.SigfússonÞ. (Reykjavík: Menningarsjóður).

  • 25

    FiskeR. S.JacksonE. D. (1972). Orientation and growth of Hawaiian volcanic rifts: the effect of regional structure and gravitational stresses.Proc. R. Soc. Lond.329299–326. 10.1098/rspa.1972.0115

  • 26

    FoulgerG. R. (1988). Hengill triple junction, SW Iceland; 2. Anomalous earthquake focal mechanisms and implications for process within the geothermal reservoir and at accretionary plate boundaries.J. Geophys. Res.93507–523.

  • 27

    FoulgerG. R.LongR. E.EinarssonP.BjörnssonA. (1989). Implosive earthquakes at the active accretionary plate boundary, Northern Iceland.Nature337640–642. 10.1038/337640a0

  • 28

    GeshiN.ShimanoT.ChibaS.NakadaS. (2002). Caldera collapse during the 2000 eruption of Miyakejima Volcano.Jpn. Bull. Volcanol.6455–68. 10.1007/s00445-001-0184-z

  • 29

    GreenR. G.GreenfieldT.WhiteR. S. (2015). Triggered earthquakes suppressed by an evolving stress shadow from a propagating dyke.Nat. Geosci.8629–632. 10.1038/ngeo2491

  • 30

    GreenR. G.WhiteR. S.GreenfieldT. (2014). Motion in the north Iceland volcanic rift zone accommodated by bookshelf faulting.Nat. Geosci.729–33. 10.1038/ngeo2012

  • 31

    GreenfieldT.WhiteR. S. (2015). Building icelandic ignieous crust by repeated melt injections.J. Geophys. Res.1201–14. 10.1016/j.tecto.2007.11.024

  • 32

    GreenfieldT.WhiteR. S.RoeckerS. (2016). The magmatic plumbing system of the Askja central volcano, Iceland, as imaged by seismic tomography.J. Geophys. Res. Solid Earth1217211–7229. 10.1002/2016JB013163

  • 33

    GreenfieldT.WhiteR. S.WinderT.ÁgústsdóttirT. H. (2018). Seismicity of the Askja and Bárðarbunga volcanic systems of Iceland, 2009-2015.J. Volcanol. Geotherm. Res.391106432. 10.1016/j.jvolgeores.2018.08.010

  • 34

    GudmundssonA. (1995). Infrastructure and mechanics of volcanic systems in Iceland.J. Volcanol. Geotherm. Res.641–22. 10.1016/0377-0273(95)92782-q

  • 35

    GudmundssonM. T.JónsdóttirK.HooperA.HolohanE. P.HalldórssonS. A.ÓfeigssonB. G.et al (2016). Gradual caldera collapse at Bárdarbunga volcano, Iceland, regulated by lateral magma outflow.Science353:aaf8988. 10.1126/science.aaf8988

  • 36

    GylfadóttirS. S.KimJ.HelgasonJ. K.BrynjólfssonS.HöskuldssonÁJóhannessonT.et al (2017). The 2014 Lake Askja rockslide-induced tsunami: Optimization of numerical tsunami model using observed data.J. Geophys. Res. Oceans122:496. 10.1002/2016JC012496

  • 37

    HarrisA.MurrayJ.AriesS.DaviesM.FlynnL.WoosterM.et al (2000). Effusion rate trends at Etna and Krafla and their implications for eruptive mechanisms.J. Volcanol. Geotherm. Res.102237–270. 10.1016/s0377-0273(00)00190-6

  • 38

    HeimissonE. R.EinarssonP.SigmundssonF.BrandsdóttirB. (2015). Kilometer- scale Kaiser effect identified in Krafla volcano, Iceland.Geophys. Res Lett.427958–7965. 10.1002/2015GL065680

  • 39

    HenschM.DahmT.RitterJ.HeimannS.SchmidtB.StangeS.et al (2019). Deep low-frequency earthquakes reveal ongoing magmatic recharge beneath Laacher See Volcano (Eifel, Germany).Geophys. J. Int.2162025–2036. 10.1093/gji/ggy532

  • 40

    HjaltadóttirS. (2009). Use of Relatively Located Microearthquakes to Map Fault Patterns and Estimate the Thickness of the Brittle Crust in Southwest Iceland. M. Sc. thesis, University of Iceland, Reykjavík. Available online at: https://skemman.is/bitstream/1946/3990/1/Sigurlaug-MS-fixed.pdf

  • 41

    HjaltadóttirS.VogfjörðK.SlungaR. (2009). Seismic signs of magma pathways Through the Crust in the Eyjafjallajökull Volcano, South Iceland. Icelandic Meteorological Office, Report VÍ 2009-013, 3. Available online at: https://www.vedur.is/media/vedurstofan/utgafa/skyrslur/2009/VI_2009_013.pdf

  • 42

    HjartardóttirÁR.EinarssonP. (2012). The Kverkfjöll fissure swarm and the eastern boundary of the Northern Volcanic Rift Zone, Iceland.Bull. Volcanol.74143–162. 10.1007/s00445-011-0496-6

  • 43

    HjartardóttirÁR.EinarssonP.BramhamE.WrightT. J. (2012). The Krafla fissure swarm, Iceland, and its formation by rifting events.Bull. Volcanol.742139–2153. 10.1007/s00445-012-0659-0

  • 44

    HjartardóttirÁR.EinarssonP.MagnúsdóttirS.BjörnsdóttirÞBrandsdóttirB. (2016a). “Fracture systems of the Northern Volcanic Rift Zone, Iceland - An onshore part of the Mid-Atlantic plate boundary,” in Magmatic Rifting and Active Volcanism, edsWrightT. J.AyeleA.FergusonD. J.KidaneT.Vye-BrownC. (London: Geological Society), 420. 10.1144/SP420.1

  • 45

    HjartardóttirÁ. R.EinarssonP.BjörnsdóttirS. G. (2016b). Fissure swarms and fracture systems within the Western Volcanic Zone, Iceland – effects of spreading rates.J. Struct. Geol.9139–53, 10.1016/j.jsg.2016.08.007

  • 46

    HjartardóttirÁR.EinarssonP.SigurðssonH. (2009). The fissure swarm of the Askja volcanic system along the divergent plate boundary of N Iceland.Bull. Volc.71961–975. 10.1007/s00445-009-0282-x

  • 47

    HollingsworthJ.LeprinceS.AyoubF.AvouacJ. P. (2013). New constraints on dike injection and fault slip during the 1975-1984 Krafla rift crisis, NE Iceland.J. Geophys. Res.: Solid Earth1183707–3727. 10.1002/jgrb.50223

  • 48

    HooperA.OfeigssonB. G.SigmundssonF.LundB.GeirssonH.EinarssonP.et al (2011). Increased capture of magma in the crust promoted by ice-cap retreat in Iceland.Nat. Geosci.4783–786. 10.1038/NGEO1269

  • 49

    HorálekJ.FischerT.EinarssonP.JakobsdóttirS. S. (2015). “Earthquake swarms,” in Encyclopedia of Earthquake Engineering, edsBeerM.PatelliE.KougioumtzoglouI. A.AuS.-K. (Berlin: Springer), 10.1007/978-3-642-36197-5_294-1

  • 50

    HudsonT. S.WhiteR. S.GreenfieldT.ÁgústsdóttirT. H.BrisbourneA.GreenR. G. (2017). Deep crustal melt plumbing of Bárðarbunga volcano, Iceland.Geophys. Res. Lett.448785–8794. 10.1002/2017GL074749

  • 51

    JakobsdóttirS. (2008). Seismicity in Iceland: 1994-2007.Jökull5875–100.

  • 52

    JohnsenG.BjörnssonA.SigurðssonS. (1980). Gravity and elevation changes caused by magma movement beneath the Krafla caldera, Northeast Iceland.J. Geophys.47132–140.

  • 53

    JuncuD.ÁrnadóttirT.GeirssonH.GuðmundssonG. B.LundB.GunnarssonG.et al (2020). Injection-induced surface deformation and seismicity at the Hellisheidi geothermal field, Iceland.J. Volc. Geotherm. Res.391:106337. 10.1016/j.jvolgeores.2018.03.019

  • 54

    KeyJ.WhiteR. S.SoosaluH.JakobsdóttirS. S. (2011). Multiple melt injection along a spreading segment at Askja, Iceland.Geophys. Res. Lett.38:L05301. 10.1029/2010GL046264

  • 55

    La FeminaP. C.DixonT. H.MalservisiR.ÁrnadóttirT. H.SturkellE.SigmundssonF.et al (2005). Geodetic GPS measurements in south Iceland: strain accumulation and partitioning in a propagating ridge system.J. Geophys. Res.110:B11405. 10.1029/2005JB003675

  • 56

    LahrJ. J. C.WardP. L. (1976). Hypoellipse: A Computer Program for Determining Local Hypocentral Parameters, Magnitude and First Motion Pattern.Reston, VA: U. S. Geological Survey. U. S. Geological Survey Open File Report.

  • 57

    LarsenG.GrönvoldK.ThorarinssonS. (1979). Volcanic eruption through a geothermal borehole at Námafjall, Iceland.Nature278707–710. 10.1038/278707a0

  • 58

    MartensH. R.WhiteR. S.KeyJ.DrewJ.SoosaluH.JakobsdóttirS. S. (2010). Dense seismic network provides new insight into the 2007 Upptyppingar dyke intrusion.Jökull6047–66.

  • 59

    MillerA. D.JulianB. R.FoulgerG. R. (1998). Three-dimensional seismic structure and moment tensors of non-double-couple earthquakes at the Hengill-Grensdalur volcanic complex, Iceland.Geophys. J. Int.133309–325. 10.1046/j.1365-246x.1998.00492.x

  • 60

    NealC. A.BrantleyS. R.AntolikL.BabbJ. L.BurgessM.CallesK.et al (2019). The 2018 rift eruption and summit collapse of Kilauea Volcano.Science363367–374. 10.1126/science.aav7046

  • 61

    PagliC. F.SigmundssonL.ÁrnadóttirT.EinarssonP.SturkellE. (2006). Deflation of the Askja volcanic system: constraints on the deformation source from combined inversion of satellite radar interferograms and GPS measurements.J. Volc. Geotherm. Res.15297–108. 10.1016/j.jvolgeores.2005.09.014

  • 62

    PálmasonG. (1963). Seismic refraction investigation of the basalt lavas in northern and eastern Iceland.Jökull1340–60.

  • 63

    PálssonM.EinarssonP.HrafnkelssonB. (2019). Variation in b-value of caldera earthquakes during recent activity of the Bárðarbunga volcano in Iceland.Jökull6971–82. 10.33799/jokull2019.69.071

  • 64

    ParameswaranR. M.ThorbjarnardóttirB. S.StefánssonR.BjarnasonI. T. H. (2020). Seismicity on conjugate faults in Ölfus, South Iceland: case study of the 1998 Hjalli-Ölfus earthquake.J. Geophys. Res. Solid Earth125:e2019JB019203. 10.1029/2019JB019203

  • 65

    PassarelliL.MaccaferriF.RivaltaE.DahmT.BokuE. A. (2013). A probabilistic approach for the classification of earthquakes as ‘triggered’or ‘not triggered’.J. Seismol.17165–187. 10.1007/s10950-012-9289-4

  • 66

    PedersenR.SigmundssonF. (2004). InSAR based sill model links spatially offset areas of deformation and seismicity for the 1994 unrest episode at Eyjafjallajökull volcano, Iceland.Geophys. Res. Lett.31:L14610. 10.1029/2004GL020368

  • 67

    PedersenR.SigmundssonF. (2006). Temporal development of the 1999 intrusive episode in the Eyjafjallajökull volcano, Iceland, derived from InSAR images.Bull. Volcanol.68377–393. 10.1007/s00445-005-0020-y

  • 68

    PedersenR.SigmundssonF.EinarssonP. (2007). Controlling factors on earthquake swarms associated with magmatic intrusions; constraints from Iceland.J. Volcanol. Geotherm. Res.16273–80. 10.1016/j.jvolgeores.2006.12.010

  • 69

    RooyakkersS. M.StixJ.BerloK.BarkerS. J. (2020). Emplacement of unusual rhyolitic to basaltic ignimbrites during collapse of a basalt-dominated caldera: The Halarauður eruption, Krafla (Iceland).GSA Bull.1321881–1902. 10.1130/B35450.1

  • 70

    SchulerJ.GreenfieldT.WhiteR. S.RoeckerS. W.BrandsdóttirB.StockJ. M.et al (2015). Seismic imaging of the shallow crust beneath the Krafla central volcano, NE Iceland.J. Geophys. Res. Solid Earth1207156–7173. 10.1002/2015JB012

  • 71

    SchulerJ.PughD. J.HaukssonE.WhiteR. S.StockJ. M.BrandsdóttirB. (2016). Focal mechanisms and size distribution of earthquakes beneath the Krafla central volcano, NE Iceland.J. Geophys. Res. Solid Earth1215152–5168. 10.1002/2016jb013213

  • 72

    SgattoniG.LucchiF.EinarssonP.De AstisG.TranneC. A. (2019). The 2011 unrest at Katla volcano: seismicity and geological context.Jökull6953–70. 10.33799/jokull2019.69.053

  • 73

    SigmundssonF. (2006). Iceland Geodynamics: Crustal Deformation and Divergent Plate Tectonics.Chichester: Springer-Praxis, 209.

  • 74

    SigmundssonF.EinarssonP.BilhamR.SturkellE. (1995). Rift-transform kinematics in south Iceland: deformation from Global Positioning system measurements, 1986 to 1992.J. Geophys. Res.1006235–6248. 10.1029/95jb00155

  • 75

    SigmundssonF.HooperA.HreinsdóttirS.VogfjördK. S.ÓfeigssonB. G.HeimissonE. R.et al (2015). Segmented lateral dyke growth in a rifting event at Bárðarbunga volcanic system, Iceland.Nature517191–195. 10.1038/nature14111

  • 76

    SigmundssonF.HreinsdóttirS.HooperA.ÁrnadóttirT. H.PedersenR.RobertsM. J.et al (2010). Intrusion triggering of the 2010 Eyjafjallajökull explosive eruption.Nature468426–430. 10.1038/nature09558

  • 77

    SigmundssonF.ParksM.PedersenR.JónsdóttirK.ÓfeigssonB. G.GrapenthinR.et al (2018). “Magma movements in volcano plumbing systems and their associated ground deformation and seismic pattern. Chapter 11,” in Volcanic and Igneous Plumbing Systems, ed.BurchardtS. (Amsterdam: Elsevier), 285–322. 10.1016/B978-0-12-809749-6.00011-X

  • 78

    SigmundssonF.PinelV.GrapenthinR.HooperA.HalldórssonS.EinarssonP.et al (2020). Unexpected large eruptions from buoyant magma bodies within viscoelastic crust. Nature Commun.11:2403. 10.1038/s41467-020-16054-6

  • 79

    SigurdssonH.SparksS. R. J. (1978). Rifting episode in North Iceland in 1874-1875 and the eruptions of Askja and Sveinagja.Bull. Volcanol.41149–167. 10.1007/bf02597219

  • 80

    SigurdssonO. (1980). Surface deformation of the Krafla fissure swarm in two rifting events.J. Geophys.47154–159.

  • 81

    SoosaluH.EinarssonP. (1996). “Torfajökull seismicity - Hydrothermal cooling of a magma body?,” in Seismology in Europe, ed.ThorkelssonB. (Reykjavík: European Seismological Commission), 467–472.

  • 82

    SoosaluH.EinarssonP. (1997). Seismicity around the Hekla and Torfajökull volcanoes, Iceland, during a volcanically quiet period, 1991-1995.Bull. Volcanol.5936–48. 10.1007/s004450050173

  • 83

    SoosaluH.KeyA. J.WhiteR. S.KnoxC.EinarssonP.JakobsdóttirS. S. (2009). Lower-crustal earthquakes caused by magma movement beneath Askja volcano on the north Iceland rift.Bull. Volcanol.7255–62. 10.1007/s00445-009-0297-3

  • 84

    StaudacherT.FerrazziniV.PeltierA.KowalskiP.BoissierP.CatherineP.et al (2009). The April 2007 eruption and the Dolomieu crater collapse, two major events at Piton de la Fournaise (La Réunion Island, Indian Ocean).J. Volcanol. Geotherm. Res.184126–137. 10.1016/j.jvolgeores.2008.11.005

  • 85

    StefánssonR.BöðvarssonR.SlungaR.EinarssonP.JakobsdóttirS.BungumH.et al (1993). Earthquake prediction research in the South Iceland seismic zone and the SIL project.Bull. Seismol. Soc. Am.83696–716.

  • 86

    SturkellE.EinarssonP.SigmundssonF.GeirssonH.PedersenR.Van DalfsenE.et al (2006). Volcano geodesy and magma dynamics in Iceland.J. Volc. Geotherm. Res.15014–34. 10.1016/j.jvolgeores.2005.07.010

  • 87

    SturkellE.SigmundssonF.GeirssonH.ÓlafssonH.TheodórssonT. (2008). Multiple volcano deformation sources in a post-rifting period: 1989-2005 behaviour of Krafla, Iceland constrained by levelling, tilt and GPS observations.J. Volcanol. Geotherm. Res.177405–417. 10.1016/j.jvolgeores.2008.06.013

  • 88

    SaemundssonK. (1974). Evolution of the axial rifting zone in Northern Iceland and the Tjörnes fracture zone.Geol. Soc. Am. Bull.85495–504. 10.1130/0016-7606(1974)85<495:eotarz>2.0.co;2

  • 89

    SaemundssonK. (1978). “Fissure swarms and central volcanoes of the neovolcanic zones in Iceland,” in Crustal Evolution in Northwestern Britain and Adjacent Regions, Vol. 10edsBowesD. R.LeakeB. E. (Washington, DC: Geological Society), 415–432.

  • 90

    SaemundssonK. (1982). “Calderas in the Neovolcanic zones of Iceland,” in Eldur er í norðri (Fires in the North), Festschrift for Sigurður Þórarinsson, eds H. Þórarinsdóttir, Ó. H. Óskarsson, S. SteinÞórsson, and Þ. Einarsson (Reykjavík: Sögufélag), 221–239.

  • 91

    SaemundssonK. (1991). “Geology of the Krafla area,” in Náttúra Mývatns, (The Nature of lake Mývatn), edsGardarssonA.EinarssonA. (Reykjavík: Icelandic Nature Science Society), 25–95.

  • 92

    TarasewiczJ.BrandsdóttirB.WhiteR. S.HenschM.ThorbjarnardóttirB. (2012a). Using microearthquakes to track repeated magma intrusions beneath the Eyjafjallajökull stratovolcano, Iceland.J. Geophys. Res.117:B00C06. 10.1029/2011JB008751

  • 93

    TarasewiczJ.WhiteR. S.WoodsA. W.BrandsdóttirB.GudmundssonM. T. (2012b). Magma mobilization by downward-propagating decompression of the Eyjafjallajökull volcanic plumbing system.Geophys. Res. Lett.39:L19309. 10.1029/2012GL053518

  • 94

    TryggvasonE. (1980). Subsidence events in the Krafla caldera, North Iceland, 1975-1979.J. Geophys.47141–153.

  • 95

    TryggvasonE. (1984). Widening of the Krafla fissure swarm during the 1975-1981 volcano-tectonic episode.Bull. Volcanol.4747–69. 10.1007/bf01960540

  • 96

    TryggvasonE. (1989). Ground deformation in Askja, Iceland: Its source and possible relation to flow of the mantle plume.J. Volcanol. Geotherm. Res.3961–71. 10.1016/0377-0273(89)90021-8

  • 97

    VogfjörðK. S.HjaltadóttirS.KjartanssonE.RobertsM. J.SlungaR. (2009). “Interpreting seismic signals from Icelandic volcanoes,” in Proceedings of the Conference on Earthquakes and Pre-Earthquake Processes, Reykjavík.

  • 98

    WardP. L.BjörnssonS. (1971). Microearthquakes, swarms and the geothermal areas of Iceland.J. Geophys. Res.763953–3982. 10.1029/jb076i017p03953

  • 99

    WesselP.LuisJ. F.UiedaL.ScharrooR.WobbeF.SmithW. H. F.et al (2019). The generic mapping tools version 6.Geochem. Geophys. Geosyst.205556–5564. 10.1029/2019GC008515

  • 100

    WhiteR. S.DrewJ.MartensH. R.KeyJ.SoosaluH.JakobsdóttirS. S. (2011). Dynamics of dyke intrusion in the mid-crust of Iceland.Earth Planet. Sci. Lett.304300–312. 10.1016/j.epsl.2011.02.038

  • 101

    WolfeC. J.BjarnasonI. T. H.VanDecarJ. C.SolomonS. C. (1997). Seismic structure of the Iceland mantle plume.Nature385245–247. 10.1038/385245a0

  • 102

    WrightT. J.SigmundssonF.AyeleA.BelachewM.BrandsdottirB.CalaisE.et al (2012). Geophysical constraints on the dynamics of spreading centres from rifting episodes on land.Nat. Geosci.5242–250. 10.1038/NGEO1428

Summary

Keywords

Iceland, Northern Volcanic Zone, earthquakes, Krafla volcano, Askja volcano, dike intrusion, rifting, inflation–deflation cycle

Citation

Einarsson P and Brandsdóttir B (2021) Seismicity of the Northern Volcanic Zone of Iceland. Front. Earth Sci. 9:628967. doi: 10.3389/feart.2021.628967

Received

13 November 2020

Accepted

24 February 2021

Published

26 March 2021

Volume

9 - 2021

Edited by

Derek Keir, University of Southampton, United Kingdom

Reviewed by

Finnigan Illsley-Kemp, Victoria University of Wellington, New Zealand; Luigi Passarelli, Université de Genève, Switzerland

Updates

Copyright

*Correspondence: Páll Einarsson,

This article was submitted to Volcanology, a section of the journal Frontiers in Earth Science

Disclaimer

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.

Outline

Figures

Cite article

Copy to clipboard


Export citation file


Share article

Article metrics