ORIGINAL RESEARCH article
Rifting Kinematics Produced by Magmatic and Tectonic Stresses in the North Volcanic Zone of Iceland
- 1Department of Earth and Environmental Sciences, University of Milan-Bicocca, Milan, Italy
- 2CRUST- Interuniversity Center for 3D Seismotectonics with Territorial Applications, Chieti Scalo, Italy
- 3Department of Human and Innovation Sciences, Insubria University, Varese, Italy
- 4Institute of Earth Sciences, University of Iceland, Reykjavik, Iceland
In the North Volcanic Zone of Iceland, we studied with the greatest possible detail the complete structural architecture and kinematics of the whole Theistareykir Fissure Swarm (ThFS), an N-S-trending, 70 km long active rift. We made about 7500 measurements along 6124 post-Late Glacial Maximum (LGM) extension fractures and faults, and 685 pre-LGM structures. We have collected the data over the last 6 years, through extensive field surveys and with the aid of drone mapping with centimetric resolution. In the southern sector of the study area, extension fractures and faults strike mainly N10°-20°, the opening direction is about N110°, and the dilation amount is in the range 0.1–10 m. In the central sector, faults and extension fractures strike mainly N00-10°, the opening direction is N90-100°, and the dilation amount is 0.1–9 m. In the northern sector, extension fractures and faults strike N30-40°, the opening direction is about N125°, and the dilation amount is 0.1–8 m. The variations in strike are attributable to two processes: the interaction with the WNW-ESE-striking Husavik-Flatey transform fault and Grímsey Oblique Rift (Grímsey lineament), and the structural inheritance of older NNE- to NE-striking normal faults. Most extension fractures show a minor strike-slip component: a systematic right-lateral component can be accounted for by the interaction with the WNW-ESE-striking fault zones and the regional, oblique opening of the rift. We regard dyke propagation as a possible cause for the more complex strike-slip components measured at several other fractures. Cumulated dilation and fracture frequency decrease along the rift with distance away from the Theistareykir volcano, situated in the central sector of the ThFS. This is interpreted as a decrease in the number of dykes that are capable of reaching great distances after being injected from the magma chamber.
Rifting is a widespread process in both submarine and terrestrial environments, the study of which enables enhancing knowledge about how plate separation works and about related geohazards such as volcanism and seismicity. Although classical rifting is assumed to take place as plate separation along a direction perpendicular to plate boundary, frequently this process is complicated by an oblique component that results in shear along the rift zone (Withjack and Jamison, 1986; Tron and Brun, 1991). Oblique opening brings about a complex network of normal as well as strike-slip faults and transtensional structures, whose kinematics and orientation are guided by the angle between the direction of relative motion of the two diverging plates and the rift trend (Corti, 2012). This complexity is made even more so by the influence of magma-induced stress that produces deformation at the surface, especially in relation to shallow dyke emplacement (Mastin and Pollard, 1988; Rubin and Pollard, 1988; Rubin, 1992; Gudmundsson et al., 2008; Tibaldi, 2015). Recently, the development of strike-slip motions has been described as taking place along fractures induced by shallow dyking at the 2014–2015 Bárðarbunga eruption in southern Iceland (Ágústsdóttir et al., 2016; Ruch et al., 2016). Consequently, it seems that both tectonic and magma-induced stresses are capable of inducing strike-slip components along rift structures. This setting is made even more complicated by the fact that rift zones with rift-parallel shear components can result in different fault patterns: on one side, oblique rifting typically produces en-échelon fault systems with segments not perpendicular to the direction of plate separation (Fournier and Petit, 2007; Figure 1A). On the other side, Karson (2017) proposed the partitioning of motions with strike-slip faults at the sides of a propagating rift, within which normal faults do occur (Figure 1B). In the latter case, the presence of strike-slip faults does not necessarily result from oblique plate separation, but it can be the effect of the outward rift propagation from a hotspot center.
Figure 1. Sketches of rift zones with rift-parallel shear component that takes place through different fault patterns: (A) oblique rifting with en-échelon faults, the strike of which is not perpendicular to plate separation direction (modified after Fournier and Petit, 2007). (B) Partitioning of motions with strike-slip faults at the sides of a rift that is propagating outward from a hotspot center (modified after Karson, 2017).
In order to shed light on the possible causes of shear components at rifts, it is necessary to carry out detailed surveys of fracture location, strike, age, and kinematics. While at continental rifts the surface rupture processes allow for a direct observation of fracture systems, in oceanic rifting settings the processes mostly occur at great depth below sea level, thus hindering detailed data collection. The North Volcanic Zone (NVZ) of Iceland represents the surface expression of the Mid-Atlantic Ridge, where plate spreading occurs at about 2 cm/yr (Geirsson et al., 2006; DeMets et al., 2010). Along the Icelandic rift system it is possible to study ongoing processes of seismic deformation and volcanism; this geodynamic setting provides a formidable natural laboratory to investigate the complex interactions between magma- and tectonic-induced stresses (e.g., Einarsson and Brandsdóttir, 1980; Björnsson, 1985; Heki et al., 1993; Jonsson et al., 1997; Cattin et al., 2005; Wright et al., 2012; Tibaldi et al., 2013; Pasquaré Mariotto et al., 2020). In the NVZ, previous authors suggested the presence of shear components along the Theistareykir Fissure Swarm (ThFS) (Karson, 2017; Tibaldi et al., 2019; Bonali et al., 2019a, b) and interaction with transform faults (Tibaldi et al., 2016a, b, 2018); this makes the ThFS an outstanding site for gaining major insights into deformation patterns at complex diverging plate boundaries.
With the purpose of contributing to a better understanding of the above, we hereby provide a complete picture of the ThFS, composed of thousands of extension fractures and faults extending for about 70 km (Figure 2). The ThFS represents the westernmost rift of the NVZ. This has been active since 8–9 Ma (Saemundsson, 1974; Young et al., 1985; Bergerat and Angelier, 2008) and is suggested to have been propagating northwards during the past million years (Gudmundsson, 2007). The goal of our work has been to define, with an unprecedented detail, the kinematics of a rift, by performing systematic field and drone surveys (Figure 3) aimed at assessing the geometry, direction and amount of opening of faults as well as extension fractures. Measuring several sites along each structure has allowed us to collect a great deal of data (about 7500 measures), instrumental for defining strike-slip components and distinguishing the different roles played by tectonic and magma-induced forces in affecting post-Late Glacial Maximum (LGM) deformation within the ThFS. These data are meant to integrate the first investigations made along the same rift by Bonali et al. (2019b), who focused only on its northern part. We are now able to present new data also from the central and southern parts of the rift; in detail, we focused on: (i) the geometric characterization of all the structures, with special reference to the strike and length of faults and fractures; (ii) the quantification of the opening directions of extension fractures and the possible presence of shear components; (iii) the estimation of the strain amount along the rift; (iv) the assessment of the origin of the structures; (v) the evaluation of the overall working of the rift.
Figure 2. Tectonic setting of northeastern Iceland. Within the North Volcanic Zone, the Mid-Atlantic Ridge is offset by the Húsavık-Flatey Fault (HFF) and by the Grimsey Lineament (GRL), which, together with the Dalvík Lineament (DAL), compose the Tjörnes Fracture Zone. Yellow stars locate historical earthquakes with M > 6 (after Stefansson et al., 2008). Orange stripes represent volcano-tectonic rift zones, also termed “fissure swarms,” where the main normal faults, extension fractures and volcanic fissures are shown (after Magnúsdóttir and Brandsdóttir, 2011; Hjartardóttir et al., 2012, 2015); Quaternary central volcanoes are highlighted as white triangles. Box locates (Figure 4A). The inset shows the location of the area and the main volcano-tectonic rift zones of Iceland (after Einarsson and Saemundsson, 1987; Hjartardóttir et al., 2016a); the black rectangle locates the Tjörnes Fracture Zone and the NVZ. The cartographic data is from the National Land Survey of Iceland AFS, Askja Fissure Swarm; DAL, Dalvík Lineament; FFS, Fremrinámar Fissure Swarm; KFS, Krafla Fissure Swarm; KR, Krafla central volcano (triangle); KOR, Kolbeinsey Ridge; HU, town of Húsavik; MA, Mánáreyjar eruptive site of 1867; SB, Skjálfandi Bay; TH, Theistareykir central volcano; ThFS, Theistareykir Fissure Swarm. Spatial reference: datum WGS 84/Geographic coordinates.
Figure 3. Examples of structures studied in the field. (A) Graben structures, with the two opposite-dipping normal faults. (B) Normal fault, person for scale. (C) UAV-captured picture of an extension fracture field, person for scale. (D) Example of piercing points along pure extension fracture, with indication of opening direction and amount of opening. (E) Example of piercing points along an oblique extension fracture, with a right lateral component. (F) Example of extension fracture with indication of real and wrong dilation, the latter caused by a funnel-shaped profile.
The NVZ is a 200 km long volcano-tectonic structure extending from the central portion of the Icelandic hotspot as far the Tjörnes Fracture Zone to the north (Figure 2). The latter is a zone of shear divided into three main, N120°-striking segments that are, from north to south: the Grímsey Oblique Rift (also known as Grímsey Lineament), the Húsavík-Flatey Fault and the Dalvík Lineament. All these three zones are seismically active and connect with the Kolbeinsey Ridge (KOR) through a complex oceanic transtensional rifting zone (Einarsson, 1991; Magnúsdóttir and Brandsdóttir, 2011).
The NVZ is composed of five, about N-S-trending, volcanic systems: the ThFS, Krafla, Fremrinámar, Askja and Kverkfjöll (Hjartardóttir et al., 2016b). Each consists of 5–20 km-wide and 60–100 km-long fracture swarms and a central volcano (Saemundsson, 1974). Other features are normal faults, eruptive fissures and extension fractures that strike parallel to the rift (Figure 2). The ThFS is an active, N-S trending, 70 km long volcano-tectonic rift that forms the westernmost edge of the NVZ, marked by the presence of normal faults and extension fractures, parallel to sub-parallel to each other, striking mostly N00-30° (Angelier et al., 1997; Acocella et al., 2000; Dauteuil et al., 2001; Hjartardóttir et al., 2012, 2016a; Pasquaré Mariotto et al., 2015). The ThFS is also dotted with a few volcanic vents, but their number is too limited to enable any statistical analysis. The rift is also home to the major Theistareykir central volcano. Along this rift, some structural changes can be noticed, such as changes in strike, most of which occur in the sector where the ThFS intersects the Husavik-Flatey Fault, the latter characterized by dominant right-lateral strike-slip movements (Gudmundsson et al., 1993; Fjäder et al., 1994; Magnúsdóttir and Brandsdóttir, 2011; Pasquaré Mariotto et al., 2015; Tibaldi et al., 2016a, b).
In the southernmost portion of the ThFS, normal faults dip mostly to the west, whereas approaching the Theistareykir volcano, faults mostly dip east and then vanish in correspondence of the youngest products of this volcano (Figure 4A; Hjartardóttir and Einarsson, 2015; Khodayar, 2018; Tibaldi et al., 2019). Between the volcano and the ThFS – Husavik-Flatey Fault intersection, normal faults dominate with a predominant westward dip direction in the eastern part of the rift and an eastward dip in the western part (Tibaldi et al., 2016a, b). The ThFS – Husavik-Flatey Fault junction was studied in detail by Gudmundsson et al. (1993), Fjäder et al. (1994), Magnúsdóttir and Brandsdóttir (2011), Pasquaré Mariotto et al. (2015), and Tibaldi et al. (2016a, b, 2018), while the Husavik-Flatey Fault was first studied by Einarsson (1958) and mapped by Saemundsson (1974), Gudmundsson et al. (1993), Garcia et al. (2002), Garcia and Dhont (2005), Gudmundsson (2007), Magnúsdóttir and Brandsdóttir (2011), Pasquaré Mariotto et al. (2015), and Tibaldi et al. (2016a, b). Hjartardóttir et al. (2012) suggested that the Husavik-Flatey Fault extends further east than the ThFS, even as far as the Krafla rift. This hypothesis is supported by Tibaldi et al. (2016a, b, 2018), who observed an interaction between the Husavik-Flatey Fault and the structures of the ThFS. North of this junction, the faults terminate and only extension fractures are observed (Tibaldi et al., 2016b; Bonali et al., 2019b). Further north, the ThFS continues offshore within the Öxarfjördur Basin (Einarsson, 2008; Hjartardóttir, 2013) and it connects to the Grimsey Lineament (Figure 2).
Figure 4. (A) Simplified geological map with indication of normal faults (red), strike-slip faults (black) and extension fractures (blue) identified in the study area. Ticks on downthrown side of normal faults. Geological units after Saemundsson et al. (2012). Spatial reference: WGS 84/UTM zone 38N. Rose diagrams show: (B) the strike of all the strike-slip segments belonging to the Husavik-Flatey Fault (HFF), (C) the strike of all structures affecting pre-LGM units, (D) the strike of all structures affecting post-LGM units, (E) the strike of faults only affecting post-LGM units, and (F) the strike of extension fractures only affecting post-LGM units. (G) Geological-structural sections following parts of transects 1, 17, and 30 of Figure 8, the exact traces of sections are represented in Figure 4A, geological units as in Figure 4A.
Along the ThFS, several depositional units crop out, represented by sedimentary sequences from Miocene to Pleistocene times, as well as lava flows of Pliocene, Pleistocene and Holocene age (Saemundsson, 1974; Garcia et al., 2002). The latest eruption here took place about 2.4 ka BP, with the emission of the “Theistareykjahraun” lava flows between the central volcano and the Husavik-Flatey Fault (Saemundsson et al., 2012; Figure 4A). Active vertical deformation at Theistareykir volcano, measured through GPS-based methods (Metzger et al., 2011), was interpreted as induced by inflation of the underlying magma chamber. During historical times, the following rifting and volcanic unrest episodes are documented: (i) between 1867 and 1868, an M 6 earthquake hit in the northern part of the ThFS, accompanied by volcanic eruptions in the northern offshore sector of the rift, within the Mánáreyjar volcanic system (Figure 2; Thoroddsen, 1925; Halldórsson, 2005; Magnúsdóttir and Brandsdóttir, 2011); (ii) between 1884 and 1885, an earthquake swarm culminated in an M 6.3 earthquake along the rift, near the coast to the north (Halldórsson, 2005); (iii) between 2007 and 2009 an 8 cm uplift in correspondence of an area centered on the Theistareykir volcano was recorded, and it was interpreted as the result of an inflation of the magma chamber below (Metzger et al., 2011; Metzger and Jónsson, 2014).
Materials And Methods
All the structures documented in the present work were first mapped on aerial photos, and then they were studied and classified during fieldwork. In order to collect data along the identified structures, both field surveys and Unmanned Aerial Vehicles (UAV) surveys were employed, the latter integrated with the aid of the Structure-from-Motion technique, enabling to produce Orthomosaics and DSMs with centimetric resolution. The structures investigated in the field were classified as: (i) normal faults, when a systematic and continuous vertical displacement > 50 cm is observed (Figures 3A,B) and (ii) extension fractures, which are dilatant fissures with a vertical component <50 cm (Figure 3C). In light of the above, the term “fracture” is used in the present work to represent the overall population composed of faults (with a vertical component > 50 cm) and extension fractures, the latter including structures with pure extensional opening (Figure 3D) or with a minor strike-slip component (Figure 3E).
In the field, wherever possible, we verified that most faults have pure dip-slip motions or a very minor strike-slip component. We also observed the presence of fissures in correspondence of the footwall block at normal faults. Wherever the fissures were located close to the associated faults, they were not included in the calculation of the opening amount; in fact, they can be regarded as due to gravitation instability induced by unbuttressing along the fault scarp. These fissures are usually quite short and can be classified as belonging to the fault structure.
UAV Use and Structure From Motion
Some areas along the rift were surveyed by means of UAVs; in particular, we used three different types of quadcopter, which is the most used and versatile type: the DJI Phantom 4, the DJI Phantom 4 PRO and the DJI Spark. The DJI Phantom 4 and the DJI Phantom 4 PRO are both equipped with a high-quality camera together with GPS coordinates, provided by the integrated Satellite Positioning Systems GPS/GLONASS (referred to the WGS84 datum). Photos are taken at 12.4 and 20 Megapixels, respectively. The UAVs can fly for a long time (28–30 min approximately for each battery). Thanks to a built-in, high precision 3-axis gimbal, flights are smoother and the collection of the images is more stable.
We also made use of the smaller DJI Spark for surveys over more limited areas. This quadcopter is equipped with a camera taking photos at 12 Megapixels, and is capable of flying for as long as 16 min.
The UAV-captured photos have been taken in such a way as to attain a 90–85% overlap along the path and a 80–75% one in a lateral direction: this is crucial to obtain a good alignment of the images and to reduce the distortions on the resulting Orthomosaics. With the DJI Phantom 4 (and the PRO version), we collected pictures between 70 and 100 m of height from the ground. With the DJI Spark, pictures were collected between 20 and 30 m of height. The collected images were then processed through Agisoft Metashape, a commercial Structure-from-Motion (SfM) software.
In order to allow for the co-registration of datasets and the calibration of the models (orthomosaics and DSMs) resulting from SfM processing, World Geodetic System (WGS84) coordinates of four artificial Ground Control Points (G) were fixed near each corner within each surveyed area; one was fixed in the central part (e.g., James and Robson, 2012; Turner et al., 2012; Westoby et al., 2012), with the purpose of reducing the “doming” effect resulting from SfM processing (e.g., James et al., 2017). The GCPs were surveyed by way of low-cost single frequency receivers in RTK configuration (with centimeter-level accuracy). The SfM technique enabled us to identify matching features in different images, collected along a defined flight path, and to combine them to create a sparse and dense cloud, a mesh, an Orthomosaic, a DTM/DSM (further details in Stal et al., 2012; Westoby et al., 2012). Particularly, our Orthomosaics and DSMs are characterized by a resolution in the following ranges: 1.0–4.5 and 2.0–18.0 cm/pixel, respectively.
For each structure, both normal faults and extension fractures, we assessed the overall strike and length in a GIS environment, where we plotted all the measured structures, mapped and classified as explained above. Other measurements, which will be described below, were made only for structures longer than 20 m.
As far as faults are concerned, in the field we determined the dip, dip-slip and strike-slip offsets by systematically measuring displacements of piercing points such as stream beds, crests, lava flows and deposits. As the downdip slip is the dominant one, we measured the vertical offset in the field at intervals of 50 m, or continuously through drones, and we were able to come up with a fault slip profile, following Manighetti’s approach (Manighetti et al., 1997, 2001). In the field, vertical offsets were measured by tape and/or laser rangefinder, where fault scarps were in the order of meters, and by GPS measurement where the scarps were in the order of tens of meters; furthermore, we assessed offsets on DSMs by measuring the difference in altitude along a topographic profile traced orthogonally to the fault scarp.
Regarding extension fractures, both in the field and on DSMs we measured strike, opening direction and net opening amount. In the field, we measured the opening direction with a compass, where corresponding piercing points on each side of the extension fracture could be located (Figures 3D,E); the net opening amount was measured along the opening direction, by means of a tape or a laser rangefinder. In case fractures were filled by sediments, we measured the amount of opening across their lower portions and only where fracture rock walls were detectable, so as to avoid any possible overestimation of the dilation amount due to gravity, frost weathering, or a combination of both, which can create a funnel-shaped profile (Figure 3F). On the orthomosaics obtained through UAV surveys, the opening direction and net amount of opening were quantified by assessing the direction and measuring the length of a line connecting two piercing points on each side of a fracture (in a GIS environment). All these measurements were made only along extension fractures that are at least 50 m away from a scarp, so as to avoid local gravity effects. Moreover, opening amounts were also collected along N109°-striking transects, each spaced 1 km in a N-S direction; this was done with the purpose of evaluating along-rift variations in the extensional strain parallel to the average direction of plate spreading affecting this part of Iceland (DeMets et al., 1994, 2010; Metzger et al., 2011, 2013; Hjartardóttir et al., 2012; Drouin et al., 2017).
Pattern of Faults and Extension Fractures
In the study area, the structures are found in pre-LGM lithostratigraphic units and in post-LGM units (Figure 4A). The 685 structures belonging to the first group, located in the western part of the area, strike mainly from N-S to NNE-SSW, and subordinately NE-SW and NW-SE (Figures 4B,C). They are represented by 679 normal faults, mostly dipping E and W, and subordinately SE to S, and 6 extension fractures, striking NNE-SSW and NW-SE. The majority of post-LGM structures (total = 6124) strike NNE-SSW, whereas minor sets strike N-S and NE-SW (Figures 4A,D). The same orientations can be recognized both analyzing extension fractures only (total = 4181) and faults only (total = 1943) (Figures 4E,F). The structures belonging to the Husavik Flatey Fault strike mostly NW-SE, and subordinately WNW-ESE and NNW-SSE (Figure 4B). In the northern part of the study area, most normal faults have downthrown to the east, as exemplified in Transect 1 of Figure 4G. In the central part, most normal faults have downthrown in the opposite direction (Transect 17 in Figure 4G). South of the Theistareykir volcano, most normal faults dip to the east again (e.g., Transect 30 in Figure 4G), while further south both fault dips are developed.
The length of pre-LGM structures is shown in Figure 5A, resulting to be generally longer than the structures affecting post-LGM units (Figure 5B); the longest pre- and post-LGM structures measure 7652 and 3062 m, respectively, with an average length of 488 and 39 m. The minimum length of both pre-LGM and post-LGM structures can be 20 m or shorter. Among post-LGM structures, normal faults present a maximum length of 3062 m and an average length of 45 m (Figure 5C), while extension fractures are 634 m long, with an average length of 36 m (Figure 5D).
Figure 5. (A) Length distribution of pre-LGM structures. (B) Length distribution of post-LGM structures. (C) Length distribution of post-LGM faults. (D) Length distribution of post-LGM extension fractures.
As can be observed in Figure 6A, the sites of measurement for a total of 1633 are distributed homogeneously in order to have a complete areal coverage of the whole rift zone. The resulting range of measures is from N 27° to N 174°, with an average of N 111° and a standard deviation of 17. By analyzsing the distribution of opening directions with respect to Easting (Figure 6B), it is possible to observe the following pattern: in the western part of the rift, directions are mostly in the N80-140° range, whereas they trend mostly between N100° and N160° in the eastern part. As regards the distribution of opening directions with respect to Northing (Figure 6C), the southern and central sections of the rift show similar ranges, mostly in the order of N80-140°, whereas the northern part shows a range mostly of N100-160°. Thus, combining the two graphs, it is possible to highlight a clockwise rotation of opening directions from the southwestern to the northeastern part of the ThFS.
Figure 6. (A) Map of opening directions (lines) obtained along extension fractures, (B) graph of opening directions (Y axis) vs. Easting (X axis), (C) graph of opening directions (X axis) vs. Northing (Y axis), the intersection with the HFF is shown, (D) graph of opening directions (Y axis) vs. extension fracture strikes (X axis). Spatial reference: WGS 84/UTM zone 38N.
The graph of opening directions with respect to extension fracture strikes shows that, as the latter rotate clockwise from NW-SE to NE-SW, so do opening directions, which rotate also clockwise from around N60° to N160° (Figure 6D). Although this indicates that the two are linked together, as will be discussed in the related section, it is worth underscoring the wide data scatter, which may indicate a more complex relation.
Figure 7A shows the strike-slip component of opening, which was observed at 1129 out of 1633 structural stations along extension fractures. This results in the majority of extension fractures being open obliquely to their strike. The structures striking from N60°W to N10°W (−60° to −10° in Figure 7A) show only right-lateral components, whereas left-lateral components are present along fractures striking N60-80°. Those striking N10°W to N60°can have both components, although it is possible to notice a gradual shift from right-lateral to left-lateral components, along with the clockwise rotation of fracture strike. It is noteworthy that only 31% of the entire extension fracture population shows pure extension. We regard as pure extension fractures those with a lateral component < 5°, whereas the others show a predominant right-lateral or left-lateral component: in view of the above, out of our 1633 total data, 504 are pure extension fractures (31%), 854 fractures show a right-lateral component (52%) and 275 show a left-lateral component (17%). The average value of the strike-slip component is + 6.4°, and the standard deviation is 14°, suggesting a predominance of the right-lateral component. The red stripes of Figures 7A,B indicate how the lateral component is expected to change in case of variations in fracture strike, according to the plate separation vector. The latter has been obtained as the average value (N109°) out of the available vectors; in the study area, based on geomagnetic data (as far back in time as 2.6 Ma) the long-term spreading vector is N104-106° (DeMets et al., 1994, 2010; Hjartardóttir et al., 2012). Based on GPS data from different time windows (1997–2011, 2006–2010, 2008–2014), the short-term spreading vectors are N109-115-112° (Metzger et al., 2011, 2013; Drouin et al., 2017).
Figure 7. (A) Graph of the lateral component of opening directions (Y axis) with respect to fracture strike (X axis), (B) opening direction (X axis) vs. lateral component (Y axis). In (A,B), red stripes indicate how the lateral component is expected to change in case of variations of the values along the X axis, according to plate separation vector range (N104-115°). (C) Graph of lateral component (Y axis) vs. Easting (X axis), (D) graph of lateral component (X axis) vs. Northing (Y axis); the intersection with the HFF is shown. Black dashed lines enclose pure extension fractures, with lateral component values <5°. Spatial reference: WGS 84/UTM zone 38N.
By plotting the opening direction azimuth with respect to the strike-slip component (Figure 7B), it is possible to note a shift from left-lateral to right-lateral component, along with the clockwise rotation of opening directions. The comparison of the strike-slip component with respect to Easting shows a slightly higher frequency of right-lateral components along the westernmost part of the rift (Figure 7C). The comparison of the strike-slip component with respect to Northing indicates the more frequent occurrence of right-lateral components northwards (Figure 7D).
We measured the extensional strain by cumulating the contribution from normal faults (Figure 8A) and extension fractures (Figure 8B). Opening amounts were measured in the field as well as on UAV-derived high-resolution images across the study area, totaling 2745 measurements (1698 along extension fractures and 1047 along normal faults). In order to better assess extensional deformation, we cumulated extension values also along N109°-striking paths, each spaced 1 km in a N-S direction, parallel to the average long-term spreading direction, obtained as explained above. At extension fractures, as already mentioned, we quantified the net opening. At normal faults, we calculated the horizontal extension based on the downdip slip and the dip of the fault plane, averaged at a value of 75° (Forslund and Gudmundsson, 1991). Extension fracture opening is in the 0.1–10 m range, with an average of 1.2 m (Figure 9A). The higher values are in the south-central part of the rift. Horizontal dilation values along normal faults are in the 0.1–55 m range, with an average of 2.4 m (Figure 9B).
Figure 8. (A) Distribution of measurements of vertical offset values at normal faults and (B) net opening amount values at extension fractures. Black lines show transects along which extension was cumulated. Spatial reference: datum WGS 84/UTM zone 28N.
Figure 9. (A) Graph showing the frequency distribution of dilation at extension fractures; (B) graph showing the frequency distribution of horizontal dilation at normal faults; (C) cumulated dilation and (D) number of fractures, (E) average fracture spacing and (F) deformation zone width; (G) stretch ratio, calculated both considering only Post-LGM units and considering also Pre-LGM units along WNW-ESE transects (shown as thin black lines in Figure 8A).
In regard to the extensional strain measured along the N109° transects, values have been calculated both considering only fractures affecting Post-LGM units and considering also fractures affecting Pre-LGM units. In the first case, the cumulated dilation is in the range of 3–6 m north of transect 11 (Figure 9C). A minimum is observed at transect 12, and then a rather steady increase can be noticed southward as far as transect 21, where the cumulated dilation is = 35 m. Further south, the cumulated dilation decreases. If we consider also fractures affecting Pre-LGM units, the cumulated dilation is in the 16–38 m range north of transect 11, a minimum is reached at transect 15, and then an increase is observed up to transect 28, where the cumulated dilation is 111.6 m. At transects 18–22, the values coincide with the ones obtained considering only Post-LGM units. Moving south, dilation decreases at transect 29, whereas another peak can be observed at transect 30 (dilation = 106.5 m). Dilation values remain greater than 70 m between transect 29 and 31, showing a gradual decrease further south.
Regarding the number of fractures calculated along the various transects, it ranges mostly around 1–20 in Post-LGM units (Figure 9D). This range is also valid in the rift section where there is the greatest value of cumulated dilation (i.e., between transect 20 and 22). This means that in this area the same amount of fractures takes up a larger deformation than in the other sections of the rift. Considering also Pre-LGM units, the number of fractures increases showing a range between 4 and 28.
The average fracture spacing is 0.40 km, with a maximum value of 1.35 km if we consider only fractures affecting Post-LGM units, whereas it presents an average value of 0.67 km and a maximum value of 1.52 km considering also Pre-LGM units (Figure 9E).
In Post-LGM units, the area affected by dilation is wider in correspondence of the central sector of the ThFS (transects 17–19) (Figure 9F). From here, it slightly increases northward from 0.9 km to 3 km, whereas it decreases southward from 4.5 to 2 km. In the southernmost sector of the rift, each transect (36–40) intersects only one fracture, setting to 0 km the average spacing between fractures and the area of deformation. Taking into account also the Pre-LGM units, the width of the area affected by dilation stays the same between transects 17–19, whereas it increases going northward up to 11.5 km (transect 1); going southward it remains greater than 6 km up to transect 37, decreasing further south.
The stretch ratio is a measure of the extensional strain of a linear element, which can be defined as the ratio between the final length and the initial length. In Post-LGM units (Figure 9G), it has been calculated on the basis of the maximum width of the area affected by fractures (as wide as 9.8 km); it ranges up to a maximum of 1.0035, with an average of 1.00084. Being generally constant in the northern part, it reaches a peak in the central part (at transect 21), and it gradually decreases in the southern part. Considering also Pre-LGM units, the maximum width of the area affected by fractures is 11.5 km, which results in a maximum stretch ratio of 1.107, with an average of 1.033. Two peaks are observed at transects 28 and 30, with a decrease in the southern part, whereas in the northernmost part an increase can be noticed.
Interaction Between Transform Faults and the ThFS Rift
The most relevant changes of fault and extension fracture strikes are observed in correspondence of the possible prolongation of the Husavik-Flatey Fault and in the northernmost part of the rift, near the sea coast in front of the Grimsey Lineament (Figure 4A). Approaching the Husavik-Flatey Fault, the rift structures bend from a dominant NNE-SSW orientation to a NNW-SSE orientation. This can be better understood in Figure 10A, in which the faults and extension fractures located along the prolongation of the Husavik-Flatey Fault are plotted, with the exception of those striking NNE-SSW, which are more frequent north and south of this area. In light of the above, our observations corroborate the preliminary results provided by Tibaldi et al. (2018) who suggested that this anticlockwise re-orientation of the rift structures is due to the possible role played by the right-lateral strike-slip Husavik-Flatey Fault underneath lava flows of Holocene age. As attested by analog modeling as well (Figure 10B), shallow right-lateral normal faults, striking NNW-SSE to N-S, can form at the surface with an en-échelon geometry, on top of a WNW-ESE transtensional right-lateral wrench fault affecting the basement.
Figure 10. (A) Plot of extension fractures and faults located exclusively along the prolongation of the Husavik-Flatey Fault (HFF) and having a different strike from the remaining structures of the ThFS, (B) analog model of shallow, right-lateral normal faults above a transtensional right-lateral wrench fault. Note the similarity with (A). Redrawn after Tibaldi et al. (2016b).
The interaction between the structures of the ThFS and the Grimsey Lineament, in the northernmost part of the rift, can be interpreted in a different manner. Here, the rift structures tend to rotate clockwise, so much as to attain a preferential NE-SW orientation. The right-lateral, transtensional Grimsey Lineament is not located below the studied ThFS structures (as opposed to the above outlined setting, marked by the intersection of the structures with the Husavik-Flatey fault), but it is found further north. Hence, we can suggest a process, according to which the structures of the ThFS are dragged along the southern part of the shear zone, attaining a heterogeneous simple shear with a smoothly increasing strain gradient northwards (Figure 11A). This setting, in fact, may result in a clockwise rotation of the structures as a consequence of the shearing deformation. In theory, an alternative interpretation may call for a rigid clockwise block rotation, that may also take place in correspondence of a right-lateral shear zone (Figure 11B; e.g., Nelson and Jones, 1987). However, in the latter case, the book-shelf rotation produces secondary zones of left-lateral shear in correspondence of the sides of the rotating blocks. On the contrary, our hundreds of field data point to a dominant right-lateral component of slip, whereas no fractures or faults with a dominant left-lateral strike-slip have been observed here. Although left-lateral strike-slip faulting has been suggested by Stefansson et al. (2008), to account for the 1885 earthquake, it is worth noting that there is no indication of fault direction in the micro-earthquake pattern of this seismic event, and only its location and magnitude have been obtained from intensity observations by Halldórsson (2005). Stefansson et al. (2008) suggested left-lateral slip along a NNE-SSW plane only on the basis of the distribution of small events associated with the large M6.2 1976 earthquake (Rögnvaldsson et al., 1998). Actually, this event took place 25 km further NE than the 1885 earthquake, in an offshore environment. The focal mechanism solution of the 1976 event might suggest left-lateral strike-slip along a NE-SW plane or right-lateral strike-slip along a NW-SE plane parallel to the Grimsey Lineament (Einarsson, 1987). Through mapping, mainly with the aid of aerial photos and satellite images, lineaments and fractures in the ThFS, Khodayar et al. (2018) suggested that the local, right-stepping en-échelon arrangement of some of the NNE-SSW fractures might indicate left-lateral strike-slip displacements. However, our detailed survey along all these fractures does not confirm this kinematics. We thus can conclude that there are no data documenting book-shelf faulting in the onshore, northern portion of the ThFS.
Figure 11. (A) Sketch of dragging of the about N-S structures located south of the Grimsey Lineament right-lateral shear zone by heterogeneous simple shear with a smoothly increasing strain gradient. (B) Sketch of rigid clockwise block rotation (book-shelf type) that also may occur in correspondence of a right-lateral shear zone (after Nelson and Jones, 1987). Model B is not reconciled with field data collected in the northern part of the study area.
Finally, also at other locations along the ThFS, a few segments of faults and fractures can be found, striking preferentially NNE to NE. In agreement with Bonali et al. (2019b), we also believe that these changes in strike orientation may derive from the structural inheritance of an older fault pattern located below the Holocene lavas. This hypothesis is based on the fact that within the conterminous, older rocks of Upper Miocene-Pleistocene age that crop out immediately west of post-LGM units (Figure 4A), two main sets of faults can be noticed, striking around N-S and NE-SW. The scarps of the NE-SW striking faults, in particular, stop abruptly at the stratigraphic contact with the post-LGM lava flows, indicating that the faults must prolong below the younger lavas. Such substrate faults may have locally guided the propagation of fractures upward with a similar NE-SW strike.
Our extremely detailed survey of the opening directions along the entire set of structures has made it possible to unveil a quite peculiar pattern. From Figure 6D it can be seen that, as the fracture strike rotates, their opening direction changes as well. This may be claimed to reflect a process that is linked to fracture orientation, and only subordinately to the remote tectonic stresses. If the tectonic least principal stress (σ3) would dominate, the opening direction should remain quite stable at fractures with different orientations. We thus suggest that the genesis of most extension fractures in the ThFS is linked to a strong local dilatant stress that forces the fracture walls to open along a direction from perpendicular to sub-perpendicular to the fracture strike. This type of stress can be produced by magma migrating in the crust along shallow-depth dykes. In fact, according to a commonly accepted model of dyke emplacement at shallow depth, the upward advancement of a dyke tip line should induce, at the surface, the development of a narrow graben or extension fractures above it (Pollard and Holzhausen, 1979; Mastin and Pollard, 1988; Rubin, 1992; Bonafede and Olivieri, 1995; Chadwick and Embley, 1998). In Iceland, such process has been observed, for example, during dyke emplacement events at Krafla (Opheim and Gudmundsson, 1989) and Holuhraun (Hjartardóttir et al., 2016b; Ruch et al., 2016). In particular, Gudmundsson (1995), thanks to field observations in Iceland, suggested that dyke emplacement does not necessarily lead to graben formation but, instead, results more often in the development of extension fractures above it. We believe that, when magma propagates horizontally (i.e., laterally) along a shallow, vertical dyke, it can produce also fractures with strike-slip components as well. This is testified to by the onset of seismic events, contemporaneous to dyke intrusion, with double-couple focal mechanism solutions showing strike-slip components of deformation, as observed in the past at Long Valley caldera (Savage and Cockerham, 1984), during the 2014 event (Shelly et al., 2016) at Yellowstone caldera (Russo et al., 2017, 2020), and during the 2014 Bárðarbunga dyking episode in Central Iceland (Ágústsdóttir et al., 2016).
Figures 6A,B show that the southern and central sections of the rift display similar opening directions, in the N80-140° range, whereas the northern part shows a systematic clockwise rotation of opening directions, that attain a N100-160° range. This is reflected in both Northing and Easting variation because the rift is oriented NNE-SSW and, therefore, the measurement sites positioned further east correspond also to those located at the northern termination of the rift. We regard this clockwise rotation as induced by the same process described in the previous chapter, i.e., the effect of shear on the northernmost onshore part of the rift, produced by the transtensional right-lateral Grimsey Lineament.
The graphs of Figures 7A,C,D demonstrate the predominance of opening directions rotated in a clockwise direction relative to fracture strike, which are expressed as the dominant right-lateral component. From a quantitative point of view, 52% of extension fractures show a right-lateral component, 31% pure extension, and 17% a left-lateral component. Figure 7A is particularly insightful, as it shows that the presence of the right-lateral component cannot be accounted for in terms of a rotation of the strike of fractures. The dominant right-lateral component, instead, is systematically observed whatever the orientation of the fractures. The red stripes indicating how the lateral component is expected to change in case of variations in the fracture strike with respect to the plate spreading direction, separates a large amount of values where the right-lateral component is present. The reason for this opening pattern, marked by a widespread right-lateral component, can be found in the model of Figure 1A, where oblique rifting is normally found at rifts with structures whose strike is not perpendicular to the direction of plate separation (Fournier and Petit, 2007). A dextral shear along a fossil dyke swarm was found at Álftafjörður, in eastern Iceland, by Eriksson et al. (2015), and was interpreted as due to a paleo-spreading direction of N120–125°, which is rotated slightly more northerly than present-day. In regard to the ThFS, another possibility is that its kinematics is also guided by the presence of the nearby Krafla rift and the interposed tectonic block. The latter has a wedge shape with the apex toward the south, where the two rifts come closer. A further complication may arise from different opening directions and/or amount along the Krafla rift. In order to look into these possibilities, we set out to map, with a very high detail, the Krafla structures as well.
Along-Rift Variation of Dilation
Along the ThFS, the cumulated dilation tends to decrease northward and southward from the rift’s central area, as can be observed both considering deformation only in Post-LGM units and considering also Pre-LGM units (Figure 9C). The variation in the number of fractures is more complex, but a general decrease outward from the rift’s center is recognizable, though there are other peaks at transects 10–31 (Figure 9D). This combination of the outward-directed decrease in fracture number and overall opening indicates that the fracturing process is more intense at the rift center near the Theistareykir volcano. In map view, the exact location of the volcano (situated between transects 25 and 27) corresponds to a lower amount of fractures, but this is due to the local presence of historic lavas that hide most structures. In Post-LGM units, the width of the deformation zone is also maximum in correspondence of the central part of the rift, while a minimum is reached in correspondence of the intersection with the HFF; going southward, values decrease (Figure 9F). Considering also Pre-LGM units, the width of the area affected by dilation is more stable, with a decrease in the southernmost sector, but an important decrease can still be noticed in correspondence of the intersection with the HFF. The outward decrease of all these features can be explained in terms of a process of faulting and fracturing whose intensity varies in space, under the influence of tectonic and magmatic stresses. In regard to tectonics, plate separation here can contribute to triggering faulting processes, especially in the case of the studied normal faults that display vertical offsets up to 200 m. For major normal faults within rift zones, Acocella and Trippanera (2016) and Trippanera et al. (2015a, b) suggested that several dyke events can induce a cumulated amount of vertical displacement at the surface, in the order of several tens of meters. However, according to the present level of knowledge of the ThFS rift, we cannot rule out the possibility that such major fault displacements could be compatible with the long-term effects of multiple rifting events induced in part by amagmatic tectonic plate spreading (fault slip not related to dyke intrusion) and in part by shallow dyke intrusion events.
Deformation may be concentrated in correspondence of the central zone of the rift because, here, the heat flux from the magma chamber of the Theistareykir volcano tends to focus the tectonic stresses, thus promoting fault slip near the volcano, as well as producing the above cited, outward-tapering pattern. A reorientation of rifting in correspondence of shallow magma chambers was documented at oceanic ridges, such as, for example at the East Pacific Rise where Rea (1978) and Martinez et al. (1997) documented a shift in the spreading center axis location. As observed also at the Afar depression (Lahitte et al., 2003), active magma chambers can lead to the concentration of the maximum tensile stress, due to high temperatures that increase ductility and reduce the yield strength and the thickness of the brittle lithosphere.
Regarding extension fractures, several may have been produced by multiple shallow dyking events. As dykes can propagate horizontally from the magma chamber, magma pressure decreases with distance, as magma is drawn from the chamber. This, in turn, brings about a decrease in the number of dykes that can reach far distances. As a consequence of this, the intensity of dyke-triggered deformation decreases outward from the central part of the rift (Buck et al., 2006).
The width variation of the deformation zone is different in a northward direction; in Post-LGM units, first there is a drop from 9 to 1 km width, and then a slight increase from 1 to 3 km toward the sea coast (Figure 9F). The same pattern can be observed if also Pre-LGM units are considered, with a drop from 9 to 3 km width, and then an increase from 3 to 11,5 km going northward. This geometry can be explained in two ways: on one side, it can depend on the fact that with distance from the magma chamber, dykes tend to become deeper and the width of the related surface deformation zone widens, as observed at the Krafla rift in Iceland by Hollingsworth et al. (2013). According to a second hypothesis, the northward-directed widening of the deformation zone may be the consequence of other dyke intrusion episodes here; in fact, the presence of a submarine volcanic system (Mánáreyjar, Figure 2) north of the ThFS (Hjartardóttir et al., 2009; Magnúsdóttir et al., 2015; Khodayar and Björnsson, 2018) may account for the southward dyke propagation from Mánáreyjar along the ThFS during the Holocene.
The calculated stretch ratio in the study area is as much as 1.0035 considering only Post-LGM units, and 1.107 if we consider also the deformation in Pre-LGM units. The first value is very similar to the stretch ratio of 1.009 estimated along the Krafla Fissure Swarm, resulting from the deformation accumulated during the last 10 ka (Dauteuil et al., 2001). Over this period, both tectonic stresses and stresses produced by dyking events have contributed to the total cumulated surface deformation. The importance of multiple dyking events in producing stretching at the surface was documented also by Paquet et al. (2007), who calculated a stretch ratio of 1.036–1.046 across the Tertiary Alftafjordur dyke swarm, produced solely by intrusion events during about 1 Ma.
Our work documents, for the first time and with unprecedented detail, the overall deformation pattern along one of the most active rifts of Iceland, the Theistareykir Fissure Swarm (ThFS). We mapped the whole population of 2622 faults and 4187 extension fractures by field surveys at 1:5000 scale and Unmanned Aerial Vehicle surveys with cm resolution. We identified 685, N-S to NE-SW striking structures affecting pre-LGM units across the western part of the study area. We focused most of our attention on the 6124, post-LGM structures, so as to compare structures of the same age.
Most post-LGM faults and extension fractures strike N-S to NNE-SSW, but in the central part of the ThFS and at its northernmost part, they rotate in anticlockwise (striking NNW-SSE) and clockwise (striking NE-SW) fashion, respectively. The first case is explained in terms of the presence, underneath Holocene lava units, of the Husavik-Flatey right-lateral strike-slip fault, which brings about fault rotation at the shallowest crustal level. The NE-SW strike of the structures in the northern ThFS is interpreted as due to the dragging effect along the block located south of the transtensional, right-lateral Grimsey Lineament, and to the structural inheritance of older faults in the substrate with the same NE-SW orientation that guided successive fracture propagation up to the surface, through post-LGM lava flows.
Along the ThFS, the cumulated dilation, expressed as the sum of the opening measured at all fractures, decreases southward and northward with respect to the central part of the rift, where the Theistareykir central volcano is located. Outward from this volcano, there is also an along-rift decrease in the number of structures, and the deformation zone becomes narrower. The above evidence can be explained by taking into consideration the outward propagation of dykes from the Theistareykir magmatic chamber; with distance, magma pressure decreases as magma is drawn away from the chamber, and this results in a decrease in the number of dykes that can reach far distances. As a consequence of this, the intensity of dyke-triggered deformation decreases outward from the central part of the rift.
The increase in the number of faults and in their dip-slip offsets at the central part of the rift, with special reference to the large faults located west of the Theistareykir volcano, is related to the heat flux connected with its magma chamber; the higher temperature increases ductility and reduces the yield strength and the thickness of the brittle lithosphere, focusing tectonic stresses and favoring fault slip. The proximity to the magma chamber is also associated with greater magma pressure, which may contribute to increasing dyke-induced fault slip especially along the rift axis, which is to say immediately north and south of the Theistareykir volcano.
Most extension fractures show components of right-lateral and left-lateral shear, with deviations as much as 40° from pure extension. This kinematics may be explained with the earlier proposed models related to different processes: (i) repeated horizontal propagation of vertical dykes that produces strike-slip components of motions, accompanied by magma inflation and separation of the fracture walls; (ii) dragging along the Husavik-Flatey Fault and Grimsey Lineament; (iii) tectonic events related to regional stresses and plate motions. Cumulated episodes may lead to a resulting very complex kinematics in a volcano-tectonic rift zone, where local processes interact with far-field tectonic stress.
Data Availability Statement
The datasets generated for this study are available on request to the corresponding author.
AT, FB, and FP wrote the first draft of the manuscript. NC, ER, PE, and ÁH wrote sections of the manuscript. All authors contributed to manuscript revision, read, and approved the submitted version.
This work was a contribution to the International Lithosphere Program – Task Force II. Funding for fieldwork were from ILP and from the MIUR project ACPR15T4_00098 – Argo3D (http://argo3d.unimib.it/). Furthermore, this study has been carried out in the framework of the Project MIUR – Dipartimenti di Eccellenza 2018-2022, and European Space Agency project no. 38829. Finally, this effort was also an outcome of GeoVires, the Virtual Reality Lab for Earth Sciences hosted at the Department of Earth and Environmental Sciences, University of Milan Bicocca (Italy). Agisoft Metashape was also acknowledged for photogrammetric data processing.
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.
We wish to acknowledge the reviewers DT and SS for their constructive insights and suggestions on an earlier version of the manuscript.
Acocella, V., Gudmundsson, A., and Funiciello, R. (2000). Interaction and linkage of extension fractures and normal faults: examples from the rift zone of Iceland. J. Struct. Geol 22, 1233–1246. doi: 10.1016/s0191-8141(00)00031-6
Ágústsdóttir, T., Woods, J., Greenfield, T., Green, R. G., White, R. S., Winder, T., et al. (2016). Strike-slip faulting during the 2014 Bárðarbunga-Holuhraun dike intrusion, central Iceland. Geophys. Res. Lett. 43, 1495–1503. doi: 10.1002/2015gl067423
Angelier, J., Bergerat, F., Dauteuil, O., and Villemin, T. (1997). Effective tension-shear relationships in extensional fissure swarms, axial rift zone of northeastern Iceland. J. Struct. Geol. 19, 673–685. doi: 10.1016/s0191-8141(96)00106-x
Bergerat, F., and Angelier, J. (2008). Immature and mature transform zones near a hot spot: the South Iceland Seismic Zone and the Tjörnes Fracture Zone (Iceland). Tectonophysics 447, 142–154. doi: 10.1016/j.tecto.2006.05.046
Bonafede, M., and Olivieri, M. (1995). Displacement and gravity anomaly produced by a shallow vertical dyke in a cohesionless medium. Geophys. J. Int. 123, 639–652. doi: 10.1111/j.1365-246x.1995.tb06880.x
Bonali, F. L., Tibaldi, A., Marchese, F., Fallati, L., Russo, E., Corselli, C., et al. (2019b). UAV-based surveying in volcano-tectonics: an example from the Iceland rift. J. Struct. Geol. 121, 46–64. doi: 10.1016/j.jsg.2019.02.004
Bonali, F. L., Tibaldi, A., Pasquaré Mariotto, F., Saviano, D., Meloni, A., and Sajovitz, P. (2019a). Geometry, oblique kinematics and extensional strain variation along a diverging plate boundary: the example of the northern Theistareykir Fissure Swarm, NE Iceland. Tectonophysics 756, 57–72. doi: 10.1016/j.tecto.2019.02.018
Buck, W. R., Einarsson, P., and Brandsdóttir, B. (2006). Tectonic stress and magma chamber size as controls on dike propagation: constraints from the 1975–1984 Krafla rifting episode. J. Geophys. Res. B12404. doi: 10.1029/2005JB003879
Cattin, R., Doubre, C., De Chabalier, J. B., King, G., Vigny, C., Avouac, J. P., et al. (2005). Numerical modelling of quaternary deformation and post-rifting displacement in the Asal–Ghoubbet rift (Djibouti, Africa). Earth Planet. Sci. Lett. 239, 352–367. doi: 10.1016/j.epsl.2005.07.028
Chadwick, W. W. Jr., and Embley, R. W. (1998). Graben formation associated with recent dike intrusions and volcanic eruptions on the mid-ocean ridge. J. Geophys. Res. Solid Earth 103, 9807–9825. doi: 10.1029/97jb02485
Corti, G. (2012). Evolution and characteristics of continental rifting: analog modelling inspired view and comparison with examples from the East African Rift System. Tectonophysics 522–523, 1–33. doi: 10.1016/j.tecto.2011.06.010
Dauteuil, O., Angelier, J., Bergerat, F., Verrier, S., and Villemin, T. (2001). Deformation partitioning inside a fissure swarm of the northern Icelandic rift. J. Struct. Geol. 23, 1359–1372. doi: 10.1016/s0191-8141(01)00002-5
DeMets, C., Gordon, R. G., Argus, D. F., and Stein, S. (1994). Effect of recent revisions to the geomagnetic reversal time scale on estimates of current plate motions. Geophys. Res. Lett.21, 2191–2194. doi: 10.1029/94gl02118
Drouin, V., Sigmundsson, F., Ófeigsson, B. G., Hreinsdóttir, S., Sturkell, E., and Einarsson, P. (2017). Deformation in the Northern Volcanic Zone of Iceland 2008–2014: an interplay of tectonic, magmatic, and glacial isostatic deformation. J. Geophys. Res.: Solid Earth 122, 3158–3178. doi: 10.1002/2016jb013206
Einarsson, P. (1987). “Compilation of earthquake fault plane solutions in the North Atlantic and Arctic oceans,” in Recent Plate Movements and Deformation, Vol. 20, ed. K. Kasahara (Washington, DC: AGU), 47–62.
Einarsson, P., and Saemundsson, K. (1987). “Earthquake epicentres 1982-1985 and volcanic systems in Iceland (map),” in Í hlutarins eðli, Festschrift for Thorbjörn Sigurgeirsson, ed. Th. I. Sigfússon (Reykjavík: Menningarsjóður).
Eriksson, P. I., Riishuus, M. S., and Elming, S. Å. (2015). Magma flow and palaeo-stress deduced from magnetic fabric analysis of the Alftafjörður dyke swarm: implications for shallow crustal magma transport in Icelandic volcanic systems. Geol. Soc. Lond. Spec. Publ. 396, 107–132. doi: 10.1144/sp396.6
Fjäder, K., Gudmundsson, A., and Forslund, T. (1994). Dikes, minor faults and mineral veins associated with a transform fault in North Iceland. J. Struct. Geol. 16, 109–119. doi: 10.1016/0191-8141(94)90022-1
Fournier, M., and Petit, C. (2007). Oblique rifting at oceanic ridges: relationship between spreading and stretching directions from earthquake focal mechanisms. J. Struct. Geol. 29, 201–208. doi: 10.1016/j.jsg.2006.07.017
Garcia, S., and Dhont, D. (2005). Structural analysis of the Húsavík-Flatey transform fault and its relationships with the rift system in Northern Iceland. Geodinam. Acta 18, 31–41. doi: 10.3166/ga.18.31-41
Garcia, S., Angelier, J., Bergerat, F., and Homberg, C. (2002). Tectonic analysis of an oceanic transform fault zone revealed by faulte slip data and earthquake focal mechanisms: the Husavik-Flatey Fault, Iceland. Tectonophysics 344, 157–174. doi: 10.1016/s0040-1951(01)00282-7
Geirsson, H., Árnadóttir, T., Völksen, C., Jiang, W., Sturkell, E., Villemin, T., et al. (2006). Current plate movements across the Mid-Atlantic Ridge determined from 5 years of continuous GPS measurements in Iceland. J. Geophys. Res. Solid Earth 111:B09407. doi: 10.1029/2005JB003717
Gudmundsson, A., Brynjólfsson, S., and Jónsson, M. T. (1993). Structural analysis of a transform fault-rift zone junction in North Iceland. Tectonophysics 220, 205–221. doi: 10.1016/0040-1951(93)90232-9
Heki, K., Foulger, G., Julian, B., and Jahn, C. H. (1993). Plate dynamics near divergent boundaries: geophysical implications of postrifting crustal deformation in NE Iceland. J. Geophys. Res. 98, 14279–14297. doi: 10.1029/93jb00781
Hjartardóttir, Á. R., and Einarsson, P. (2015). The interaction of fissure swarms and monogenetic lava shields in the rift zones of Iceland. J. Volcanol. Geotherm. Res. 299, 91–102. doi: 10.1016/j.jvolgeores.2015.04.001
Hjartardóttir, Á. R., Einarsson, P., and Sigurdsson, H. (2009). The fissure swarm of the Askja volcanic system along the divergent plate boundary of N Iceland. Bull. Volc. 71, 961–975. doi: 10.1007/s00445-009-0282-x
Hjartardóttir, Á. R., Einarsson, P., Bramham, E., and Wright, T. J. (2012). The Krafla fissure swarm, Iceland, and its formation by rifting events. Bull. Volcanol. 74, 2139–2153. doi: 10.1007/s00445-012-0659-0
Hjartardóttir, Á. R., Einarsson, P., Gudmundsson, M. T., and Högnadóttir, T. (2016b). Fracture movements and graben subsidence during the 2014 Bárðarbunga dike intrusion in Iceland. J. Volcanol. Geotherm. Res. 310, 242–252. doi: 10.1016/j.jvolgeores.2015.12.002
Hjartardóttir, Á. R., Einarsson, P., Magnusdóttir, S., Bjornsdóttir, ϸ, and Brandsdóttir, B. (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. Geological Society, Vol. 420, eds T. J. Wright, A. Ayele, D. J. Ferguson, T. Kidane, and C. Vye-Brown (London: Special Publications), 297–314. doi: 10.1144/sp420.1
Hjartardóttir, Á. R., Hjaltadóttir, S., Einarsson, P., Vogfjörd, K., and Belart, J. M. C. (2015). Fracturing and earthquake activity within the Prestahnúkur fissure swarm in the Western Volcanic Rift Zone of Iceland. J. Geophys. Res. Solid Earth 120, 8743–8757. doi: 10.1002/2015jb012341
Hollingsworth, J., Leprince, S., Ayoub, F., and Avouac, J. P. (2013). New constraints on dike injection and fault slip during the 1975–1984 Krafla rift crisis, NE Iceland. J. Geophys. Res. Solid Earth 118, 3707–3727. doi: 10.1002/jgrb.50223
James, M. R., and Robson, S. (2012). Straightforward reconstruction of 3D surfaces and topography with a camera: accuracy and geoscience application. J. Geophys. Res. Earth Surface 117:F03017. doi: 10.1029/2011JF002289
James, M. R., Robson, S., d’Oleire-Oltmanns, S., and Niethammer, U. (2017). Optimising UAV topographic surveys processed with structure-from-motion: ground control quality, quantity and bundle adjustment. Geomorphology 280, 51–66. doi: 10.1016/j.geomorph.2016.11.021
Jonsson, S., Einarsson, P., and Sigmundsson, F. (1997). Extension across a divergent plate boundary, the Eastern Volcanic Rift Zone, south Iceland, 1967–1994, observed with GPS and electronic distance measurements. J. Geophys. Res. 102, 11913–11929. doi: 10.1029/96jb03893
Karson, J. A. (2017). The Iceland plate boundary zone: propagating rifts, migrating transforms, and rift-parallel strike-slip faults. Geochem. Geophys. Geosyst. 18, 4043–4054. doi: 10.1002/2017gc007045
Khodayar, M., and Björnsson, S. (2018). Structures and styles of deformation in rift, ridge, transform zone, oblique rift and a microplate offshore/onshore North Iceland. Int. J. Earth Sci. 9, 461–511. doi: 10.4236/ijg.2018.98029
Khodayar, M., Björnsson, S., Kristinsson, S. G., Karlsdóttir, R., Ólafsson, M., and Víkingsson, S. (2018). Tectonic control of the Theistareykir geothermal field by rift and transform zones in North Iceland: a multidisciplinary approach. Open J. Geol. 8, 543–584. doi: 10.4236/ojg.2018.86033
Lahitte, P., Gillot, P. Y., and Courtillot, V. (2003). Silicic central volcanoes as precursors to rift propagation: the Afar case. Earth Planet. Sci. Lett. 207, 103–116. doi: 10.1016/s0012-821x(02)01130-5
Magnúsdóttir, S., Brandsdóttir, B., Driscoll, H., and Detrick, R. (2015). Postglacial tectonic activity within the Skjálfandadjúp Basin, Tjörnes Fracture Zone, offshore Northern Iceland, based on high resolution seismic stratigraphy. Mar. Geol. 367, 159–170.
Manighetti, I., King, G. C. P., Gaudemer, Y., Scholz, C. H., and Doubre, C. (2001). Slip accumulation and lateral propagation of active normal faults in Afar. J. Geophys. Res. Solid Earth 106, 13667–13696. doi: 10.1029/2000jb900471
Manighetti, I., Tapponnier, P., Courtillot, V., Gruszow, S., and Gillot, P.-Y. (1997). Propagation of rifting along the Arabia–Somalia plate boundary: the gulfs of Aden and Tadjoura. J. Geophys. Res. 102, 2681–2710. doi: 10.1029/96jb01185
Martinez, F., Hey, R. N., and Johnson, P. D. (1997). The east ridge system 28.5∘32”S East Pacific Rise: implication for overlapping spreading center development. Earth Planet. Sci. Lett. 151, 13–31. doi: 10.1016/s0012-821x(97)00095-2
Mastin, L. G., and Pollard, D. D. (1988). Surface deformation and shallow dike intrusion processes at Inyo Craters, Long Valley, California. J. Geophys. Res. 13, 13221–13235. 221–13,235. doi: 10.1029/jb093ib11p13221
Metzger, S., and Jónsson, S. (2014). Plate boundary deformation in North Iceland during 1992-2009 revealed by InSAR time-series analysis and GPS. Tectonophysics 634, 127–138. doi: 10.1016/j.tecto.2014.07.027
Metzger, S., Jónsson, S., and Geirsson, H. (2011). Locking depth and slip-rate of the Húsavík Flatey fault, North Iceland, derived from continuous GPS data 2006-2010. Geophys. J. Int. 187, 564–576. doi: 10.1111/j.1365-246x.2011.05176.x
Metzger, S., Jónsson, S., Danielsen, G., Hreinsdóttir, S., Jouanne, F., Giardini, D., et al. (2013). Present kinematics of the Tjörnes Fracture Zone, North Iceland, from campaign and continuous GPS measurements. Geophys. J. Int. 192, 441–455. doi: 10.1093/gji/ggs032
Opheim, J. A., and Gudmundsson, A. (1989). Formation and geometry of fractures, and related volcanism, of the Krafla fissure swarm, northeast Iceland. Geol. Soc. Am. Bull. 101, 1608–1622. doi: 10.1130/0016-7606(1989)101<1608:fagofa>2.3.co;2
Pasquaré Mariotto, F. A., Bonali, F. L., Tibaldi, A., Rust, D., Oppizzi, P., and Cavallo, A. (2015). Holocene displacement field at an emerged oceanic transform-ridge junction: the Husavik-Flatey Fault-Gudfinnugja Fault system, North Iceland. J. Struct. Geol. 75, 118–134. doi: 10.1016/j.jsg.2015.03.011
Rögnvaldsson, S. T., Gudmundsson, A., and Slunga, R. (1998). Seismotectonic analysis of the Tjörnes Fracture Zone, an active transform fault in north Iceland. J. Geophys. Res. 103, 30117–30129. doi: 10.1029/98jb02789
Russo, E., Tibaldi, A., Waite, G. P., Bonali, F. L., Massin, F., and Farrell, J. (2020). Unrevealing the complex deformation pattern at Yellowstone plateau through seismicity and fracture analysis. Tectonophysics 778:228352. doi: 10.1016/j.tecto.2020.228352
Russo, E., Waite, G. P., and Tibaldi, A. (2017). Evaluation of the evolving stress field of the Yellowstone volcanic plateau, 1988 to 2010, from earthquake first-motion inversions. Tectonophysics 700, 80–91. doi: 10.1016/j.tecto.2017.02.009
Saemundsson, K., Hjartarson, A., Kaldal, I., Sigurgeirsson, M. A., Kristinsson, S. G., and Vikingsson, S. (2012). Geological map of the Northern Volcanic Zone, Iceland. Northern Part 1: 100.000. Reykjavik: Iceland GeoSurvey and Landsvirkjun.
Savage, J. C., and Cockerham, R. S. (1984). Earthquake swarm in Long Valley caldera, California, January 1983: evidence for dike inflation. J. Geophys. Res. Solid Earth 89, 8315–8324. doi: 10.1002/9781118782095.ch9
Shelly, D. R., Ellsworth, W. L., and Hill, D. P. (2016). Fluid-faulting evolution in high definition: connecting fault structure and frequency-magnitude variations during the 2014 Long Valley Caldera, California, earthquake swarm. J. Geophys. Res. Solid Earth 121, 1776–1795. doi: 10.1002/2015jb012719
Stal, C., Bourgeois, J., De Maeyer, P., De Mulder, G., De Wulf, A., Goossens, R., et al. (2012). “Test case on the quality analysis of structure from motion in airborne applications,” in Proceedings of the 32nd EARSeL Symposium: Advances in Geosciences, Mykonos Island.
Stefansson, R., Gudmundsson, G. B., and Halldorsson, P. (2008). Tjörnes fracture zone. New and old seismic evidences for the link between the North Iceland rift zone and the Mid-Atlantic ridge. Tectonophysics 447, 117–126. doi: 10.1016/j.tecto.2006.09.019
Tibaldi, A., Bonali, F. L., and Pasquaré Mariotto, F. A. (2016a). Interaction between transform faults and rift systems: a combined field and experimental approach. Front. Earth Sci. 4:33. doi: 10.3389/feart.2016.00033
Tibaldi, A., Bonali, F. L., Einarsson, P., Hjartardóttir, Á. R., and Pasquarè Mariotto, F. A. (2016b). Partitioning of Holocene kinematics and interaction between the Theistareykir Fissure Swarm and the Husavik-Flatey Fault, North Iceland. J. Struct. Geol. 83, 134–155. doi: 10.1016/j.jsg.2016.01.003
Tibaldi, A., Bonali, F. L., Pasquaré Mariotto, F. A., and Ranieri Tenti, L. M. (2019). The development of divergent margins: insights from the North Volcanic Zone, Iceland. Earth Planet. Sci. Lett. 509, 1–8. doi: 10.1016/j.epsl.2018.12.019
Tibaldi, A., Bonali, F. L., Pasquaré Mariotto, F. A., and Russo, E. (2018). Interplay between inherited rift faults and transcurrent structures: insights from analogue models and field data from Iceland. Glob. Planet. Change 171, 88–109. doi: 10.1016/j.gloplacha.2018.03.009
Tibaldi, A., Bonali, F. L., Pasquaré, F., Rust, D., Cavallo, A., and D’Urso, A. (2013). Structure of regional dykes and local cone sheets in the Midhyrna-Lysuskard area, Snaefellsnes Peninsula (NW Iceland). Bull. Volcanol. 75, 764–780.
Trippanera, D., Acocella, V., Ruch, J., and Abebe, B. (2015a). Fault and graben growth along active magmatic divergent plate boundaries in Iceland and Ethiopia. Tectonics 34, 2318–2348. doi: 10.1002/2015tc003991
Trippanera, D., Ruch, J., Acocella, V., and Rivalta, E. (2015b). Experiments of dike-induced deformation: insights on the long-term evolution of divergent plate boundaries. J. Geophys. Res. Solid Earth 120, 6913–6942. doi: 10.1002/2014jb011850
Turner, D., Lucieer, A., and Watson, C. (2012). An automated technique for generating georectified mosaics from ultra-high resolution unmanned aerial vehicle (UAV) imagery, based on structure from motion (SfM) Point clouds. Remote Sens. 4, 1392–1410. doi: 10.3390/rs4051392
Westoby, M. J., Brasington, J., Glasser, N. F., Hambrey, M. J., and Reynolds, J. M. (2012). ‘Structure from-Motion’ photogrammetry: a low-cost, effective tool for geoscience applications. Geomorphology 179, 300–314. doi: 10.1016/j.geomorph.2012.08.021
Wright, T. J., Sigmundsson, F., Ayele, A., Belachew, M., Brandsdottir, B., Calais, E., et al. (2012). Geophysical constraints on the dynamics of spreading centres from rifting episodes on land. Nat. Geosci. 5, 242–250. doi: 10.1038/NGEO1428
Keywords: rift, Iceland, fault, fracture, stress, dyke
Citation: Tibaldi A, Bonali FL, Pasquaré Mariotto F, Corti N, Russo E, Einarsson P and Hjartardóttir ÁR (2020) Rifting Kinematics Produced by Magmatic and Tectonic Stresses in the North Volcanic Zone of Iceland. Front. Earth Sci. 8:174. doi: 10.3389/feart.2020.00174
Received: 27 February 2020; Accepted: 06 May 2020;
Published: 02 June 2020.
Edited by:Andrea Billi, National Research Council (CNR), Italy
Reviewed by:Daniele Trippanera, King Abdullah University of Science and Technology, Saudi Arabia
Susanta Kumar Samanta, Jadavpur University, India
Copyright © 2020 Tibaldi, Bonali, Pasquaré Mariotto, Corti, Russo, Einarsson and Hjartardóttir. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Federico Pasquaré Mariotto, firstname.lastname@example.org