Diverging plate boundaries consist of both rift and transform segments , regardless of whether they are above a hotspot or not. Spreading is accommodated by stretching across the rift and by shear motion along the transform zones. Classically, rifts/ridges are identified by regular magnetic anomalies, e.g. , and normal fault-controlled basins with hydrocarbon accumulation. The axial rifts/ridges are also a locus of tectonically emplaced calderas, e.g., , and geothermal resources, e.g., , where a single rift episode consists of dyke injection and related earthquakes . Transform zones, on the other hand, are primarily identified by earthquakes and shear motions, e.g., . However, the deformation becomes complex with time, and develops features such as diverse ridge geometry and failed rift segments , instability in rift and transform zones , anomalous magnetic chrons due to rift relocation, e.g., , or varied degree of rift obliquity . Due to their impacts on basin formation, petroleum system development, volcanism and hydrothermal fluid circulation, these boundaries have been studied in continental and ocean crusts of various ages, offshore as onshore.
Iceland is a window into the tectonic evolution of diverging plate boundaries that span in age from upper Tertiary to present, and in extent from land to shelf. Decades of intensive instrumental monitoring increased the understanding of active rifting, transform faulting, magma injection and crustal deformation. Studies of older structures, at a regional scale as well as offshore, have received less attention comparatively, thus leaving gaps in the literature. The present tectonic analysis contributes towards filling such a gap in an area of 16,500 km2 covering several rift, ridge, oblique rift and transform segments as well as a microplate (Figure 1). This paper presents:
・ A comprehensive and continuous regional tectonic map of the Kolbeinsey Ridge (KR), Northern Rift Zone (NRZ), the Tjörnes Fracture Zone (TFZ), and a major microplate obtained from aerial images, bathymetric and topographic maps.
・ An overview of the styles of deformation and further insights into the tectonic configuration of the plate boundaries and the microplate.
Figure 1. (a) Diverging plate boundaries, intraplate tectonic, and earthquakes in Iceland; (b) The study area in North Iceland and the earthquake swarms within the TFZ; (c) Main geological elements spanning upper Miocene-present; (d) The rift fissure swarms of the Northern Rift Zone, the three segments of the Tjörnes Fracture Zone, and focal mechanisms of selected earthquakes. Abbreviations on Figures (b) to (d) are, K-R: Kolbeinsey Ridge; Má: Mánáreyjar; Th: Theistareykir; K: Krafla; Fr: Fremrinámar; A: Askja; GOR: Grímsey Oblique Rift; HFF: Húsavík-Flatey Fault; DL: Dalvík Lineament; Stó: Stóragrunn; H: Hóllinn; Na: Nafir; Gr: Grímsey; E: Eyjafjarðaráll; Ska: Skagafjörður; Tr: Tröllaskagi; E: Eyjafjörður; Fl: Flateyjarskagi; F: Flatey; Sk: Skjálfandi; T: Tjörnes; Öx: Öxarfjörður; St: Stórihver. On Figure (c), the sub-marine hydrothermal sites of Eyjafjörður are after , the K-R and GOR after . Pockmarks offshore are after , and hydrocarbon gas seeps onshore after . Geothermal field in Öxafjörður is after , and in Theistareykir after .
・ Correlation of the observed tectonics with selected surface and depth data and their interpretation for a better understanding of the deformation.
Not only are the results of this multidisciplinary structural analysis relevant to the geology of Iceland and other rift and transform zones, but they also provide a knowledge base for exploration of geothermal or hydrocarbon resources in the context of diverging plate boundaries in general.
2. Geological Settings
2.1. Diverging Plate Boundaries at Large
As rifting, transform faulting and oblique rifting are all present in the study area, the recall of common complexities associated with these boundaries is necessary as a structural frame.
1) Rifting, transform faulting, and their instabilities are known to occur in both young and old oceanic ridges and in continental rifts. One of the complexities at these boundaries is when spreading is accommodated by a set of propagating and receding rift segments, leading to a “V” shape geometry of magnetic anomalies and pseudofaults on the flanks of the spreading centres .
2) In bookshelf tectonic between overlapping rift segments , or within transform zones, e.g., , conjugate dextral and sinistral strike-slips faults accommodate the rotation . However, regional block rotations between plate boundaries, e.g., , whether banking the spreading centres, e.g., , across transform zones, e.g., , or along rift axis, e.g., , would cause sinuosity and bend in the tectonic/magmatic structures .
3) Oblique rifting occurs - ) when the relative plate motion is not parallel to the directions of the principal strain, e.g., , and the plates spread obliquely compared to the ridge/rift trend . Rift obliquity can be moderate (15˚ - 45˚) to high (45˚ - 75˚) , but extension is accommodated by both the normal and the strike-slip faults , and characterised by en échelon segments that are surface expression of deeper structures - . Oblique rifts are considered to initiate above reactivated non-aligned major crustal weaknesses , where strain is localised , and magmatism is present at depth . Such rifts can evolve with or without transform/transfer segments . The migrating deformation due to increased rift obliquity not only leads to gradual curvature of rift/ridge axis and their internal en échelon segments , but also to stress field rotation/re-orientation/permutation in polyphase tectonics .
4) Fracturing with extensional faults and conjugate strike-slips is intense in areas under the influence of rift/transform or oblique rifting - . The en échelon segmentation and fracture sinuosity are typical fracture geometries above deeper shear faults, and even along individual strike-slip faults . The relay zones of normal faults in rifts, e.g. or the stepovers between en échelon strike-slip segments in transform zones are some of the other complexities. The stepovers could evolve into pull-apart basin (strike-slip regime) or rhomb-graben (transtension or oblique extension) , and asymmetric basins can be bounded by normal faults on one side and parallel strike-slip faults on the other side in weak transform zones subject to oblique-slips . These structures are among the most prolific for hydrocarbon and geothermal exploration. In addition, the partitioning of the damage zones of strike-slips faults control the fluid flow and permeability along individual structures.
2.2. Iceland and Its Northern Shelf
Above a hotspot , Iceland is linked to the Tertiary North Atlantic Igneous Province that is the second magmatic phase since the opening of North Atlantic about 62 - 58 Ma . The oldest igneous series in Iceland are 16 - 15 Ma , but the island has been in the making since 24 Ma or 26 Ma when the active rift jumped westwards from Aegir to the K-R and propagated onshore with a complex tectonic history.
Presently, the American and Eurasian Plates separate along the plate boundaries in Iceland and on its shelf (Figures 1(a)-(d)) at a rate of 2 cm/yr in the direction N105˚E . The active rifts are the K-R and NRZ, the Eastern Rift Zone (ERZ), and the Western Rift Zone (WRZ), also called Reykjanes-Langjökull Rift Zone. The transform zones are the TFZ, and the South Iceland Seismic Zone between the WRZ, ERZ and south of the Hreppar Microplate. Oblique rifting occurs along the GOR , the Reykjanes Peninsula and the Reykjanes Ridge with typical en échelon segmentation , V-shape geometry and pseudofaults .
Rift segments comprise fissure swarms of normal faults, open fractures and eruptive dykes, as well as calderas, acidic rocks and high temperature geothermal fields . Transform zones are loci of earthquakes , continuous seismicity , and shear motions . Intraplate volcanism and earthquakes occur outside active plate boundaries. Past rift jumps led to several extinct rifts that stretch presently from central Iceland towards the northwest , and to abandoned Tertiary transform segment . The main geological elements of study area are as follows (Figures 1(b)-(d) and Figure 2(a)):
1) The K-R is a narrow ridge, while the NRZ is 70 km wide, up to 200 km long with several fissure swarms. Only the westernmost swarm of the NRZ (Theistareykir-Mánáreyjar) falls within the study area. The TFZ is 20 km long and 70 km wide, and comprises three NW-WNW striking segments, i.e., the GOR, the Húsavík-Flatey Fault (HFF) and the Dalvík Lineament (DL). The accumulated dextral motion of the TFZ since upper Tertiary is 100 km 
Figure 2. The study area. (a) Combined Multibeam, Single beam and Spot 5 images superimposed on Digital Elevation Model (DEM) as the basis for the structural interpretation. The Multibeam bathymetric data are from , the Single beam from Icelandic Coast Guard (Sjómælingar), and the DEM (ISNET 93) from National Land Survey of Iceland (Landmælingar Íslands); (b) Distribution of the mapped structures in the study area; (c) Location of six selected areas for more detailed analysis of the deformation. Abbreviations for tectonic structures same as for Figure 1.
or 140 km , and the subsidence in the order of hundreds of metres. The offshore basins are Öxarfjörður, Skjálfandi, Eyjafjarðaráll, Eyjafjörður and Skagafjörður. The HFF cuts the Theistareykir Fissure Swarm, e.g., , stretches offshore and intersects the Skjálfandi and Eyjafjarðaráll Basins. The Holocene marine volcanic edifices align on the K-R, and in the four en échelon segments of GOR, namely Stóragrunn, Hóllinn, Nafir and Mánáreyjar, indicating dextral motion across this oblique rift, e.g., . Except a small ridge in Skjálfandi , the main highs are Grímsey and Mánáreyjar offshore, and Tjörnes, Flateyjarskagi and Tröllaskagi on land. The Grímsey High is caught between the HFF, GOR and the Eyjafjarðaráll Basin. Although the kinematics of this high is poorly known, it has been speculated that the Grímsey High is a microplate as it is caught between these plate boundary segments .
2) Onshore, the Miocene crust is 15 Ma in North Iceland and at least 10 Ma in the study area , but offshore, magnetic anomalies show an undisturbed spreading along the K-R for the past 11 - 20 Ma, e.g., . The landward continuation of this ridge is the Skagi paleo-rift (predecessor to NRZ), active at about 12 Ma . Opinions diverge on whether this paleo-rift jumped eastward to the NRZ around 6/7 Ma and the TFZ initiated then , or the NRZ became active at 8 - 8.5 Ma and the TFZ at 9 Ma . Regardless of timing, the plate boundary reorganisations caused uplift and subsidence, unconformities, flexuring, major block rotation , and intense fracturing in the study area.
3) The Miocene - present series consist of basaltic lavas, hyaloclastites, acidic rocks, and local sediments onshore. The igneous crust is covered unconformably by a 1-km thick Plio-Pleistocene sequence including the fossil rich marine Tjörnes beds and intercalated lignite, alternating lavas, thin glacial tillites and fluvio-lacustrine sediments . The Tertiary-Pliocene accumulations fill the offshore basins up to a thickness of 4 km . They are covered by the 29,000 - 14,000 yr BP moraines from the last glaciations , and by younger marine sediments. Older series are deformed at depth , while chirp and backscatter data suggest that the youngest rifting episodes are mostly between 12,000 and 3000 yr BP .
4) About 2/3 of the present transform motion across the TFZ is accommodated by GOR and 1/3 by the HFF . Micro-seismicity is continuous along the HFF, GOR, and occasionally on individual faults (Figure 1(b) and Figure 1(d)), with earthquakes reaching M6 to M7 . The focal mechanisms show dextral strike-slip on the WNW/NW and sinistral strike-slip on NNE/NE Riedel shears in the TFZ , with an occasional oblique-slip component. Earthquakes are also linked to magma injection within rift segments such as the 1976 Kópasker earthquakes ), triggered by dyke injections along the Krafla Fissure Swarm, e.g., , and the 2014 eruption in Bárðarbunga . Magmatism also occurs on strike-slip faults of the transform zone such as the 2400 yr BP eruption along a WNW blind shear fault in Theistareykir .
5) The main high temperature fields of the study area (Figure 1(c)) are in Öxarfjörður on land, e.g., , as well as in Theistareykir . Geothermal activity offshore was observed during dredging locally at K-R as well as to the west and east of Grímsey, and hydrothermal vents have also been discovered in a few places . Pockmarks are known offshore in GOR , and onshore to the north of the Krafla and Theistareykir fissure swarms where hydrocarbon gases emanate . The gases likely originated from marine sediments and lignite and maturated by the high geothermal gradient in Skjálfandi (110˚C/km) and Öxarfjörður (200˚C/km) Basins .
6) Most of the recent structural mappings, however, focused on the N-S extensional fractures and eruptive dykes of the rift/ridge . The Riedel shears were identified locally offshore, e.g., and particularly the WNW and NW sets have been mapped in places onshore . Geothermal manifestations are considered to be fracture-controlled . Their detailed tectonic settings, however, are overall poorly known, except in Theistareykir where a multidisciplinary tectonic investigation demonstrated that six sets of extensional and shear fractures belonging to the rift and transform zones control this resource .
3. Data and Method
The mapped upper Tertiary-present structures are fractures (faults, open fractures, prominent joints), and volcanic edifices. Considering the size of the area and the challenges in identifying the structures, a combination of several processed Multibeam/Single beam bathymetric maps, Digital Elevation Model (DEM), and Spot 5 satellite photos (natural, black and white, infrared channels) (Figures 2(a)-(c)) was used to detect the regional tectonic structures. As fracture identification relies on light and shadow on images, a structure may be easily detected or not at all, or its trace could appear differently on various image types. The variety of images used enhanced the chances of fracture identification. All images were reported as single layers in GIS, having true North and no distortion. The structures were handpicked and drawn directly in GIS where the range of image resolutions allowed observations down to metre-scale.
To provide an overview of the deformation in time and space, the multidisciplinary tectonic analysis is organised in this paper as follows:
・ The comprehensive regional structural pattern of the rifts, transform zone and oblique rift is presented on new regional tectonic maps with a continuous coverage (Figure 3(a)).
・ Features such as the fracture sets, trace, length, geometry, and motions are described (Figure 4), and the frequency of each set calculated in an overall statistical analysis. Some fractures are mildly curved, in which case their strikes were defined as a straight line between the two end-points of the fracture traces. Furthermore, up to 5% of the mapped fractures may not be of tectonic origin, but such a small fraction does not skew the overall frequency of the fracture sets in the statistical analysis.
・ Six areas that reflect best the variety of tectonic structures at and between rift, transform and oblique rift segments are selected on the overall tectonic map to describe in detail the styles of deformation (Figure 2(c) and Figures 5-10).
・ The extensional and shear fractures (Figure 3(a)) are interpreted in terms of the main highs, basins and volcanism (Figure 11), along with their implications for plate boundaries and a microplate (Figures 12(a)-(e)). The interpreted topographic, Multibeam/Single beam bathymetric cross-sections, as well as the seismic reflection profiles show the extent of the identified deformation from the surface down to greater crustal depths (Figures 13(a)-(g)). The correlation of the deformation and the earthquakes reflects the relation of the active deformation and the underlying major structures, with hints on fault reactivation (Figures 14(a)-(d)).
・ The main conclusions (Figures 15(a)-(k)) and some discussions summarise the contribution of this paper.
4.1. The Structural Pattern
About 2258 fracture segments were identified on the images, representing a total length of 8069 km (Figure 3(a)). The fractures cut the volcanic-sedimentary series of the upper Miocene (10 - 3.3 Ma), Plio-Pleistocene (3.3 - 0.8 Ma), upper Pleistocene-Holocene (0.8 Ma-present) and the postglacial lavas (15,000 - 2400 yr BP) onshore, but offshore they displace the postglacial formations on the seafloor. The marine Holocene volcanism are large edifices, seamounts and plateaux, while the Pleistocene - late Holocene eruptive craters onland are much smaller in size. The main features of the structures on the map are as follows:
1) As fractures are at various stages of evolution since the upper Tertiary, they have more or less strong traces in the topography. Their traces are generally stronger (sharp) in the older bedrock as faults accumulated slips, leading to a more differentiated topography. The younger Holocene faults and open fractures have subtler traces. Fractures are generally segmented along total length (Figure 3(a)) as they have not yet coalesced into single major faults.
2) Segments of an individual fault can present opposite dip directions because faults (and dykes) in Iceland are generally steeply-dipping in the range of 75˚ to 85˚ . The HFF is one of the best examples as segments of this WNW fault to the west of Flatey offshore dip towards the northwest, but dip towards southwest between Flatey and Theistareykir onshore (Figure 3(a) and Figure 4(a)).
Figure 3. Interpreted tectonic pattern in rift, transform, and oblique rift offshore and onshore North Iceland. (a) The six fracture sets, presented in different colours; (b) The statistical analysis of the observed fractures showing the frequency of each fracture set (rose diagram) and their percentage in the total fracture population (table).
Figure 4. (a) The main WNW transform segments of the HFF in the upper Tertiary outcrops to the north of Botnsvatn Lake; (b) The younger left-stepping en échelon segments of the HFF in the postglacial lavas, intersecting the N-S normal faults of Guðfinnugjá; (c) and (d) Spot 5 images showing the en échelon arrangements of open fractures, respectively, along NW dextral and ENE sinistral faults (modified from ); (e) and (f) Respectively, a saucer-shape sill (modified from ), and multiple dyke injection (Khodayar, unpublished data), in the 1-km and 1.5 km eroded crust of West Iceland (see text for explanation).
Figure 5. The main ENE shear zones of Area 1 offshore. (a) and (b) Different processing of Single beam images, presented in UTM and ISNET 93, respectively, showing that the same structure can appear differently on images; (c) The ENE sinistral shear zones cutting Skjálfandi, Mánáreyjar and Öxarfjörður, identified based on the right-stepping en échelon arrangement of short NNE segments and sinuosity of N-S and NNE segments along these zones.
Figure 6. The main NE/ENE shear zones of Area 2 offshore. (a) and (b) Different processing of Single beam images, presented in UTM and ISNET 93, respectively; (c) The NE/ENE sinistral shear zones cutting the Eyjafjarðaráll Basin, identified based on the right-stepping en échelon arrangement of the short NNE segments and sinuosity of the N-S and NNE segments. The southernmost shear zone is in the continuation of the NE horst to the north of Skagafjörður-Tröllaskagi. To the south of the Eyjafjarðaráll Basin, the NNE faults dip to the northwest within two of the shear zones.
Figure 7. The main shear zones of Area 3 offshore. (a) and (b) Multibeam/Single beam images presented in UTM and ISNET 93, respectively; (c) The N-S normal as well as the NNE/ENE sinistral and the WNW/NW dextral shear faults with sinuous and en échelon arrangements form horsts and grabens within the NW-trending GOR; (d) Interpretation of the main normal and shear zones in GOR.
Figure 8. Area 4 offshore and onshore. (a) and (b) Single beam images offshore, respectively, in UTM and ISNET 93, along with Spot 5 and DEM onshore; (c) Interpretation of the faults shown on Figure (b); (d) The same image types covering a wider structural frame for the interpretation; (e) and (f) Two alternatives for the change in the fault strike from N-S to ENE.
Figure 9. Along-axis “Ridge to Basin” configuration of the rift structures in Area 5. (a) and (b) Maps showing that the N-S ridges pass into basin along axes, and that the volcanism is mostly on the ridges and lesser within the basins; (c) Cross-sections along the axes of the N-S ridges and basins; (d) Closer look at the small N-S grabens, and the main WNW, NW and ENE faults in Theistareykir (see text for explanation).
Figure 10. Some deformation associated with the HFF in Area 6 onshore. (a) and (b) Respectively, Spot 5 and DEM showing the HFF, the surrounding faults, and a few water bodies; (c) The deformation along the left-stepping en échelon oblique-slip segments of the HFF, the rhomb-grabens occupied by lakes, and the NW horsetail splays; (d) and (e) Cross-sections showing the tectonic configuration of the rhomb-grabens and the water bodies.
Figure 11. The interpreted tectonic map of the structures at and between plate boundaries offshore and onshore. The simplified map shows the main highs, basins, shear zones, volcanism, hydrothermal sites and pockmarks offshore, as well as the main geothermal fields and hydrocarbon seeps onshore.
Figure 12. Interpretation of the regional structures in terms of plate boundaries and a microplate. (a) Distribution of volcanism in the K-R, and in the NW en échelon segments of GOR from Stóragrunn to Mánáeyjar; (b) The Öxarfjörður Basin as a likely structure of GOR; (c) The structures to the west and east of Grímsey along which the plates separate; (d) The structures at the plate boundaries, the Grímsey-Tjörnes-Dalvík major microplate, as well as the Skjálfandi and Eyjafjörður intra-microplate basins; (e) The rhomb-shaped microplate bounded by plate boundaries.
Figure 13. Interpretation of the structures at depth. (a) Location of the topographic/bathymetric cross-sections and 2D vintage seismic reflection profiles; (b) The morphostructures along cross-section A-A’; (c) and (d) The interpreted and raw seismic profile F-4-85; (e) The morphostructures along cross-section B-B’; (f) and (g) Interpretation and raw seismic profile F-6-85. Vertical scales exaggerated on all sections to show the details. On interpreted sections A-A’ and B-B’, the thick lines are the fault traces of Figure 13(a), and the thin dashed lines the detailed fracture segments of Figure 3(a).
Figure 14. Correlation of present-day earthquakes and underlying upper Tertiary-present structures. (a) Raw map of 1995-2016 earthquake swarms (SIL network); (b) Seismic lineations in the earthquake swarms and the mapped structures; (c) Selected focal mechanisms of earthquakes; (d) The stress map of North Iceland (modified from ), showing complex normal and strike-slip fault regimes; (e) The earthquake swarms of GOR may correspond either to a blind splay of GOR, or to initiation of a transform segment parallel to the HFF and DL.
Figure 15. Summary of the regional tectonic structures and the styles of deformation. (a) Simplified tectonic map of the main basins, highs and shear zones at the rift, transform and oblique rift segments offshore/onshore North Iceland; (b) The plate boundaries, the major GTD Microplate, two intra-microplate basins, and the earthquakes in GOR above a possible blind WNW structure; (c) to (e) Surface expression of deeper shear structures; (f) to (k) The range of deformation associated with the extensional and shear fractures that have fully broken the surface (see text for explanation).
3) The statistical analysis of the fractures shows the existence of six sets striking N-S, NNE, ENE, E-W, WNW and NW (Figure 3(b)). Individual NW segments are the longest, reaching up to 30 km such as the western boundary fault of the Öxarfjörður Basin (Figure 3(a)), followed by single N-S fault segments with a maximum length of 20 km. Single NNE and ENE fractures are the shortest, with a length of up to 2 km at most.
4) All fractures display dip-slip, with values ranging from 1 m to several hundreds of metres onshore, and between <10 and 100 m on the seafloor, except in GOR where the dip-slip reaches locally up to 200 m. Particularly in the upper Tertiary crust, the fractures with strongest traces are also those with the highest dip-slip (Figure 4(a) and Figure 4(b)). The segmented N-S faults (N175˚ to N20˚E) are purely extensional, without coherent en échelon arrangements. For the lack of marker horizons, the en échelon geometry and fracture sinuosity were used as criteria to identify the shear motions. The WNW (N108˚ - N135˚E) and NW (N145˚ - N170˚E) segments display systematic left-stepping en échelon arrays typical of dextral strike-slip motion, while the ENE (N51˚ - N77˚E) and NNE (N21˚ - N30˚E) segments have a right-stepping arrangement indicative of sinistral strike-slip motion (Figures 4(b)-(d)). These shear fractures are, however, often oblique-slip due to their combined strike- and dip-slips. No criterion determined satisfactorily the shear motion along the E-W set (N90˚ - 99˚E). The six described sets are clearly a combination of fractures resulting from rifting and transform faulting, which was first identified in Theistareykir in a recent structural analysis .
5) Similar to Theistareykir , the fracture population reflects here the presence of several deep-seated structural zones, expressed at the surface either as areas of higher fracture density, and/or concentration of fractures with en échelon geometry.
6) The N-S faults of the rift are the single most frequent structures with 39% of the fracture population (Figure 3(b)), but the Riedel shears of the transform zone and the oblique rift represent 61%. Among the Riedel shears, the NW set is the most frequent (21%), and the E-W set the least (3%). The E-W, NW and N-S sets are the deepest structures, as interpreted from resistivity anomalies .
4.2. Styles of Deformation
The study area has some 2000 m elevation difference between land and seafloor as the highest peak is 1289 m in Tröllaskagi onshore and the deepest part of the seafloor is about - 650 m in the Eyjafjarðaráll Basin (Figure 2(a)). Six areas are selected (Figure 2(c)) to show in detail the styles and variety of tectonic deformation, most of which are reported for the first time by this study.
4.2.1. Area 1
Area 1 is offshore and covers a great portion of the Grímsey High and the Skjálfandi Basin, as well as the totality of the Mánáreyjar High (Figure 5(a) and Figure 5(b)). A series of ENE, NNE, NW and N-S fractures fall within three parallel ENE zones that cross the Skjálfandi Basin diagonally and stretch into the south of the Grímsey High (Figure 5(c)). In these ENE zones, fractures of the four sets have more subtle traces with far less apparent dip-slips than the NW and N-S boundary faults of the Öxarfjörður and Skjálfandi Basins. The ENE and NNE/NE segments are organised in typical right-stepping en échelon arrays compatible with a sinistral strike-slip motion across the ENE zones (Figure 5(c)). The sinuosity of the N-S, NNE and NW short fractures within the two northernmost ENE zones is compatible with fault drag across the sinistral ENE zones. The traces, geometry and comparatively small dip-slips (<50 m) of these short segments indicate that the fractures are at a younger stage of development, acting very likely as surface expression of deep-seated ENE strike-slip shear zones.
4.2.2. Area 2
Area 2 in the Eyjafjarðaráll Basin displays similar deformation to Area 1 (Figures 6(a)-(c)). The Eyjafjarðaráll Basin strikes N-S over 2/3 of its length but turns to WNW at the approach of the HFF where it is bounded by faults of the same strike (Figure 2(a), Figure 3(a) and Figure 6(c)). In the N-S portion of the basin, shorter NNE and ENE segments are organised in right-stepping en échelon arrays. They are distributed in three possible parallel NE/ENE zones, each with <5 km width. At the approach of the boundaries of these ENE zones, a few of the short N-S, NE, ENE and NNE segments display gentle sinuosity, reflecting fault drags by the sinistral strike-slip motion of the ENE zones.
The southernmost NE/ENE zone in particular hosts remarkable deformation where some of the N-S faults of Eyjafjarðaráll are rather NNE (Figure 6(c)). These NNE faults have dip-slips, and have been interpreted as rift-parallel normal faults in the literature. However, they all dip towards north-northwest, which excludes their role as the boundary faults of the basin. These short segments are the surface expression of deep-seated sinistral strike-slip structure(s) crossing the basin. Furthermore, the southernmost NE zone is in exact continuation of a horst to the north of Skagafjörður-Tröllaskagi. The Eyjafjarðaráll Basin changes its strike from N-S to WNW precisely after intersecting this NE shear zone.
4.2.3. Area 3
Area 3 shows the complexity of tectonic along the GOR, part of its flanks, and the Skjálfandi Basin (Figures 7(a)-(d)). From north of Hóllinn to the north of Skjálfandi, the fracture pattern is dominated by a high number of short NW striking faults, in the presence of N-S, NNE, and a few WNW segments. All fracture sets have dip-slips, and the left- and right-stepping en échelon arrangements indicate dextral motion on NW and WNW fractures, and sinistral motion on NNE to ENE segments (Figure 7(c)). The NW, ENE, and N-S oblique-slip faults control a few local basins in Nafir and Hóllinn, along with volcanism.
The concentration of high fracture density, en échelon segments, and fracture sinuosity indicate that the deformation is partitioned in several extensional and shear zones in the GOR (Figure 7(d)). Most of the mapped fractures fall within two parallel NW deformation bands with dextral strike-slip characters, which are located towards the outer edges of the GOR, while the remaining mapped structures coincide with shorter NNE-ENE sinistral, WNW dextral and N-S extensional zones inside this oblique rift. These configurations indicate that extension at this plate boundary is accommodated by the overall dextral motion along the NW striking GOR as well as by smaller ENE sinistral, WNW dextral, and N-S extensional zones.
4.2.4. Area 4
Area 4 extends from Mánáreyjar offshore to Tjörnes and Flateyjarskagi onshore, and displays relevant features related to the rift (Figures 8(a)-(c)) and its interaction with the transform zone (Figures 8(d)-(f)). Three aspects are analysed in these relations. First is the series of N-S normal faults stretching from Mánáreyjar to the Theistareykir Fissure Swarm. No major N-S fault is identified in the western half of the Tjörnes High and the mapped N-S faults occupy the eastern halves of the Mánáreyjar and the Tjörnes Highs. Judging from the displacement of the seafloor and the topography, a great majority of these faults are eastward dipping with dip-slips in the range of <10 and to 100 m. As the faults dip to the east, they exclude a typical horst structure for these two highs. Instead, the N-S faults form the western flank of the Theistareykir/Öxarfjörður Basin (Figures 8(b)-(d)). Secondly, a portion of the offshore faults on the Mánáreyjar High and southeast of it bends from N-S to NNE at the approach of the Öxarfjörður Basin. These systematic bends are typical fault drags induced by the dextral strike-slip motion on the NW-striking western boundary fault of the Öxarfjörður Basin (Figure 8(b)-(c) and Figure 8(e)). Thirdly, a zone of tightly parallel ENE faults crosses the N-S normal faults to the north of the HFF on the Tjörnes High (Figure 8(d)). The majority of these ENE faults have pronounced trace and dips to the southeast, with apparent dip-slips up to 100 m. Towards the west of this ENE zone, a few short ENE segments display a mild right-stepping en échelon arrangement typical of sinistral motion. At Thríhyrningur to the east, however, the ENE zone isolates a small wedge composed of N-S faults, which dip mainly eastwards. Although the ENE fractures are widespread and form a separate set in the adjacent Theistareykir area to the south , the described deformation here is rather specific. Two explanatory alternatives are suggested:
1) Either the ENE zone is an independent sinistral shear zone cutting diagonally across the N-S normal faults of Tjörnes. In this case, the N-S faults of Thríhyrningur are a portion of the Mánáreyjar-Theistareykir N-S structures, detached by and dragged along the ENE sinistral zone (Figure 8(e)).
2) Or, the ENE faults were originally a part of the N-S faults of Mánáreyjar-Theistareykir, which rotated clockwise to ENE on Tjörnes north of the HFF (Figure 8(f)). Likewise, on Flateyjarskagi to the south of the HFF, a high number of Tertiary dykes are also known to have rotated clockwise . The clockwise rotation of those dykes on one side of the HFF, and of the faults on the other side are regional fault drags compatible with the dextral motion of this transform segment. In both alternatives, the N-S faults of Thríhyrningur would remain coherently a small detached wedge, caught between the ENE sinistral zone and the HFF dextral fault. As dykes feed the lavas in which they are found , and the lavas are Miocene in age here, the dyke rotation likely dates from an early stage of transform development. On the other hand, the fault rotation is, or has lasted, at least until Quaternary, i.e., the age of the youngest rocks cut by the described N-S faults.
4.2.5. Area 5
Area 5 presents interesting features related to the N-S regional structures of the K-R (Area 5a), Mánáreyjar, Tjörnes and Theistareykir (Area 5b) (Figures 9(a)-(d)). Two points should be noted prior to the analysis. Mánáreyjar High strikes overall N-S but its northern part is rather NW within the GOR. The total length of each ridge and basin can be slightly greater than shown on the maps.
The K-R and the Mánáreyjar Highs are pierced by volcanic edifices, but each of these highs evolves into a N-S basin to the south. This configuration is labelled here as “along-axis Ridge to Basin” (Figures 9(a) and Figure 9(b)). In the case of the K-R and the Eyjafjarðaráll Basin, the proportion of ridge to basin is balanced as the K-R is about 70 km long and Eyjafjarðaráll 65 km on map. An along-axis cross-section shows that the K-R is some 400 m higher in altitude than the shallowest part of the Eyjafjarðaráll Basin, but 600 m higher than the deepest part of that basin. The volcanism is concentrated on the K-R High and not in the subsided basin.
Similarly, the Mánáreyjar and Tjörnes Highs evolve along-axis into the Theistareykir low to the south of the HFF. Although the western part of Theistareykir is topographically a high (flank), the structure is overall equivalent to a basin as most of it is subsided. Additionally, it contains two narrow N-S basins with a 2.7 and 0.7 km width, cutting the postglacial and the 2400 yr BP lavas (Figures 9(b) and Figure 9(d)). Despite the general configuration “Ridge to Basin” here, there are a few differences compared to the K-R/Eyjafjarðaráll. First, the proportion of ridge to basin is uneven as the ridge/high is ~90 km long and the basin ~25 km. Secondly, the elevation is overall a mirror image of the K-R/Eyjafjarðaráll. While the top of the Tjörnes High is 400 m above the surface of the N-S basins in Theistareykir, the top of the Mánáreyjar High is 400 m lower than those basins (Figure 9(c)). Thirdly, the volcanic edifices are concentrated mostly on the Mánáreyjar High, but smaller eruptive craters exist also in the “basin” part in Theistareykir. A closer look at the structures in Theistareykir reflects interesting structural controls of some geological processes. For example, the last eruption of 2400 yr BP occurred on the Stórihver Fault (StF), which is a WNW blind dextral strike-slip fault of the transform zone and is perpendicular to this N-S rift fissure swarm . Additionally, the geothermal activity is controlled by the N-S rift-parallel faults, but more dominantly by the Riedel shears of the transform zone . These are the NW and WNW dextral oblique-slip faults (e.g., the Tjarnarás Fault and the horsetail splay of the StF), as wells as the ENE and NNE deep-seated sinistral strike-slip faults (Figure 9(d)). All of these sets are equally widespread in the present study area.
4.2.6. Area 6
Area 6 displays some detailed deformation along a 16 km length of the HFF onshore. It is emphasised that Spot 5 (Figure 10(a)) and DEM images (Figure 10(b)) provide a different view and depth to the observations, and the structural interpretation (Figure 10(c)) takes into account the average strikes and fault geometry seen on both image types. Also, only the structures relevant to the discussions are selected from the original map (Figure 3(a)) to be shown here (Figure 10(c)).
Three features exhibit the main styles of deformation. First, the HFF consists of several WNW segments, most of which dip to the southwest (Figure 10(c)). A few NW, N-S and NNE segments are also present. All fractures have dip-slip, while the left-stepping en échelon segmentation indicates dextral motion along the HFF and other WNW and NW faults. Secondly, the NW faults are a separate and widespread set. Locally, however, WNW and NW segments can act as horsetail splay of each other, as along the HFF here and in the Theistareykir area .
Thirdly, there are three water bodies (Botnsvatn, Höskuldsvatn and an unnamed pond), which are <50 m deep (Figures 10(a)-(e)). To the west, two WNW segments of the HFF dip towards each other and form a 1-km wide graben occupied by Botnsvatn. The main transform segment is to the north, dips southwest, and has ≥200 m dip-slip (Figure 4(a) and Figure 10(c)), while the southern segment dips to the northeast with ≤40 m apparent dip-slip (Figure 10(d)). A small hill, controlled by a NW fault dipping to the southwest, bounds the graben to the east. From there the graben structure decreases in width to 200 m and dies out shortly. Farther east, the elevation drops regionally about 400 m from north of the small pond towards Höskuldsvatn through a series of WNW and NW horsts and grabens (Figure 10(e)). Two of these WNW segments form a very narrow graben where the small pond is pinched in. Similar to Botnsvatn, the Höskuldsvatn is flanked to the north by a southwest dipping segment of the HFF with ≤20 m dip-slip, and to the east by a NW fault, dipping to the west. The Botnsvatn and the Höskuldsvatn are shallow lakes with likely flat bottoms, and have a configuration similar to other narrow pull-apart basins, e.g., . In fact, both these lakes have been interpreted as such by . However, Höskuldsvatn sits in a WNW shallow depression in an eroded flat land without a major fault controlling its southern flank. This configuration indicates that the lake might be the expression of a dissymmetric extensional structure. As the two lakes are at the release zones of WNW and NW strike-slip faults, they are likely rhomb-graben formed at the bends of these segmented shear faults.
5.1. The Structural Architecture
1) The detailed segmentations of individual faults (Figure 3(a)) are summarised as single fracture traces on the interpreted tectonic map (Figure 11) where the basins and highs are highlighted in colours, and the regional structures of Areas 1 to 6 are also reported. The tectonic pattern shows that the N-S normal faults parallel to the rift/ridge co-exist with the NNE/ENE sinistral and the WNW/NW dextral strike- and oblique-slip Riedel shears of the transform zone and the oblique rift. This pattern extends homogeneously from land to sea and is associated with all structures at and between the plate boundaries. In general, the long parallel ENE shear zones cut the regional offshore and onshore structures diagonally, while the other sets contribute more to basin and high formation.
2) Most of the study area consists of highs into which the basins are carved. The Eyjafjarðaráll, Skjálfandi and Öxarfjörður major basins are predominantly offshore, but as the seafloor elevates along their axes, the structures continue as topographic faulted lows onshore. The N-S faults flank most of Eyjafjarðaráll, Skjálfandi, and the grabens in Theistareykir. As to the Riedel shears, the NW dextral oblique-slip faults delimit Öxarfjörður as well as a narrow basin merging with west of Eyjarfjarðaráll. In the continuation of this latter NW graben, other NW segments extend from the Grímsey High to Flateyjarskagi onshore. After the bend of the Eyjafjarðaráll to WNW along the dextral oblique-slip of the HFF, the basin bifurcates to the south of this transform segment into narrower NE and NW graben structures. To its north, Skjálfandi bends gently to NW within GOR and merges with shorter ENE, NNE, WNW, NW and N-S basins there. But to the south, Skjálfandi splits into narrow NW basins onshore. The width of the offshore basins is rather similar, i.e., 13 km in Eyjafjarðaráll (WNW section), 12.5 km in Skjálfandi (N-S portion), and 10 - 15 km in Öxarfjörður. The NW-striking Eyjafjörður basin is the narrowest with ~3 km.
3) Due to fault block compartmentalisation, the biggest high appears to be Grímsey-Flateyjarskagi, which has an overall NW trend and extends from sea to land. The outer edges of the Grímsey High are likely bounded by N-S normal, NNE sinistral and NW dextral oblique-slip faults offshore, and those of Flateyjarskagi by NW faults onshore. The HFF and DL both cut the Grímsey-Flateyjarskagi High. Mánáreyjar-Tjörnes is the other major high bounded by N-S faults to the east. As earlier stated, the western boundary of the Tjörnes High lacks clear faulting onshore (Figure 8(c)), but offshore it is likely flanked by the N-S fault of Skjálfandi, thus defining a horst structure for the Mánáreyjar-Tjörnes High. Another high is the narrow horst to the north of Skagafjörður-Tröllaskagi, bounded by the NE sinistral oblique-slip faults (Figure 11).
4) The volcanic edifices of the K-R are organised N-S, and those of Stóragrunn, Hóllinn and Nafir in GOR are known to be within the NW left-stepping en échelon segments, reflecting a dextral motion along this oblique rift . Two volcano-tectonic features are interpreted in the present study. One is that the volcanism is in fact distributed in three NW en échelon segments, which include the totality of Mánáreyjar (Figure 12(a)). The other is the specific organisation of the volcanic edifices within each of the NW segments. The volcanic edifices in Stóragrunn and Hóllinn are likely within N-S left-stepping en échelon structures, also typical of dextral motion, while the smaller edifices in Nafir are aligned WNW and NW, or they are on shorter ENE sinistral faults. In Mánáreyjar, however, volcanism is clearly N-S to the south but WNW/NW farther north. This change of strike from N-S to WNW/NW occurs at the intersection of this volcanic segment and the ENE sinistral shear zones of Area 1.
5) Geothermal activity and pockmarks are recognised mostly onshore and found locally offshore . Although the offshore manifestations are fracture-controlled, their tectonic settings are poorly known in detail. Their position compared to the structures of the tectonic map, however, indicates that the same extensional faults and Riedel shears controlling the basins, horsts and volcanism (Figure 11) likely control also these manifestations. To reiterate, these N-S normal, NNE/ENE sinistral and WNW/NW dextral oblique-slip fault sets play the critical role in the geological processes and the geothermal activity of the prominent Theistareykir geothermal site onshore .
5.2. Configuration of Plate Boundaries and the Microplate
Chirp data (Multibeam and high resolution shallow seismic reflection) in the Skjálfandi Basin indicate that the seafloor and some 50 m of underlying Holocene sediments are displaced <100 m, vertically . The following interpretations, however, concern the deformation that has been ongoing over longer geological time scale and is more complex than the present-day or post-glacial tectonic activity.
1) The Mánáreyjar volcanic segment of the GOR is not a relay zone between the Skjálfandi and Öxarfjörður Basins (Figure 12(a)). The Öxarfjörður Basin strikes parallel to the Mánáreyjar volcanic segment, is flanked by NW dextral oblique-slip faults, and falls coherently in the continuation of the GOR. Therefore, the basin is likely one of the en échelon segments of this oblique rift (Figure 12(b)). As a consequence, the GOR would be at least 150 km long, before merging with the NRZ onshore.
2) The number of major structures on each side of Grímsey is uneven and raises questions about the configuration of the plate boundaries and the extent of the microplate. It is emphasised that due to unclear magnetic anomalies to the north of TFZ, e.g., , and the absence of geodetic measurements offshore, the following interpretation is purely speculative and relies on the geometry of the structures and kinematic considerations.
The American Plate separates from west of Grímsey across the 10 - 13 km wide K-R/Eyjafjarðaráll. To the east of Grímsey on the Eurasian Plate, spreading occurs in the presence of the Skjálfandi and Öxarfjörður Basins with a total width of ≤27.5 km (Figure 12(c)). As the GOR is at the plate boundary, contains Öxarfjörður, and extends into the active NRZ onshore (Figure 12(d)), Skjálfandi Basin falls outside of the plate boundary. The DL is the last known segment of the transform plate boundary to the south onshore. It is likely that the small NNE basin to the north of Skagafjörður and Tröllaskagi is also a part of the western plate boundary of Eyjafjarðaráll/K-R. Given these configurations, the microplate caught between these plate boundaries is much larger than the Grímsey High alone, as it would include the highs of Mánáreyjar, Tjörnes, Flateyjarskagi, and Tröllaskagi. This major microplate, labelled as Grímsey-Tjörnes-Dalvík Microplate here (GTDM), presents an overall rhomb geometry (Figure 12(e)). Skjálfandi and Eyjafjörður are two oblique-slip basins inside this microplate. It should be noted that a similar rhomb structure bounded by the K-R/Eyjafjarðaráll, GOR, NRZ, and the HFF instead of the DL, appears in the present-day kinematics monitored by GPS measurements, earthquakes and volcanism .
5.3. The Deformation at Depth
The interpreted topographic and Multibeam/Single beam bathymetric cross-sections, as well as seismic reflection profiles show the extent of the mapped structures into the upper crust (Figures 13(a)-(e)). The section A-A´ (GOR to Tröllaskagi) is 80 km long, the section B-B’ (Eyjafjarðaráll to Öxarfjörður) 120 km, but the seismic profiles cover only a portion of these sections. The stratigraphy on the seismic profiles is equivalent to the onshore Plio-Pleistocene Tjörnes sediments above the upper Tertiary basement. Both the summarised fault traces (Figure 11) and the detailed fracture segments (Figure 3(a)) are reported on the four sections to reflect the contribution of all structures at depth. The selected seismic lines are two of the six 2D vintage profiles acquired in 1985 offshore North Iceland. The six profiles were subject to an initial tectonic interpretation , and the line F-4-85 was recently re-interpreted for the NAGTEC Atlas . The interpretation of the seismic profiles below uses the same velocities for the stratigraphy as , and benefits from the comprehensive regional tectonic maps, which were not available to those previous interpretations.
1) The 1.8 seconds two-way travel time (TWT) indicative of sediment thickness on the Iceland Shelf , correlates overall with the younger sedimentary sequence showing regular reflectors above a rough basement (Figures 13(c)-(d) and Figures 13(f)-(g)). The contact between the two formations, however, undulates in detail, reflecting an erosional character. It is as shallow as 0.5 s (TWT), representing a sedimentary sequence with ~0.45 km thickness (F-4-85), and as deep as ~2.8 s (TWT) corresponding to a thicker sequence of up to ~3.5 km (F-6-85). In the upper sediments, thin formations carve locally into the regular reflectors and resemble channel infilling. But on F-6-85, three additional units, likely sedimentary, appear between the basement and the upper regular reflectors (Figure 13(c) and Figure 13(d)). The westernmost of these units (1) is without reflectors and passes laterally to the reflectors of unit 2 that is covered by the reflectors of the youngest unit 3. The units 1 and 3 are each ~1 km thick, and unit 2 ~ 0.8 km. The unit 2 in particular, is heavily broken in its top.
2) The mapped N-S normal, the NE/ENE sinistral and the WNW/NW dextral oblique-slip faults (Figure 3(a) and Figure 11) control the basins and highs down to at least 2.8 s (TWT) or 3.5 km depth (Figures 13(b)-(f)). Faults display antithetic, synthetic and hourglass geometries, and shallow basins are fully subsided in the upper 1 km of the crust. Where the Riedel shears control the basins, the grabens are within flower structures above deeper strike- to oblique-slip faults such as parts of Eyjafjörður, Skjálfandi and Öxarfjörður Basins. On all sections, the number of faults is clearly greater at depth compared to the structures mapped on the seafloor, as some of the deep faults do not reach the surface. The dip-slips of the faults are also greater at depth (hundreds of metres), but much less on the seafloor (≤100 m in general but up to 200 m on the NW boundary fault of Skjálfandi in GOR).
3) Three major faults (F1 to F3) with wide deformation zones control most of the tectonic-sedimentary processes on each seismic profile. On sections A-A’ and F-4-85, they are the HFF (F1), which dips towards the northeast offshore, as well as the NW-striking dextral oblique-slip faults F2 and F3 to the east of the Eyjafjörður Basin on Flateyjarskagi (Figure 13(c) and Figure 13(f)). On the foot wall of F1, a narrow horst is pinched between this fault and F2 and the top of the basement is uplifted to ~0.5 s (~0.45 km depth). The F2 and F3 form a graben, which has antithetic, synthetic and hourglass internal faults, and sits above a basement elevated to 1 s (0.9 km depth). On the hanging wall of F1, the basement top is at 1.8 s (~1.8 km depth), and the upper sedimentary reflectors are relatively flat on the GTDM, but severely deformed across F1. There, the WNW section of the Eyjafjarðaráll appears as a narrow basin near the surface but it widens to 5 km between 0.5 s and 1 s (0.45 and 0.9 km depth) and includes a local horst. At these depth intervals, the sedimentary formations are dragged downward by the motion of F1 and thicken within the basin, reflecting the syn-sedimentary activity of this fault.
On sections B-B’ and F-6-85, the N-S faults F1 to F3 extend from the sediments to the basement (Figure 13(e) and Figure 13(f)). The F1 and F2 dip westward and have parallel segments. The F3 is eastward dipping in its lower part, but its segments are shifted to the east above 2 s (~2.1 km) in the upper sediments. The sedimentary reflectors are flatter on the foot walls of F1 and F3 where the deformation is milder. On the platform to the east of F1, the reflectors dip gently westward where an extensional flower structure above a deep strike-slip fault zone is the most obvious structure. On the platform to the west of F3, the middle reflectors dip westward and are covered unconformably by shallower horizontal formations, but at their basis, the sedimentary unit 1 is titled eastward towards F3 and resembles a major landslide. The F1 to F3 define the Eyjafjarðaráll Basin, which is filled by sediments dipping eastward towards the F2 and F1. The F2 and F3 deepened this basin down to 2.8 s (~3.5 km) and created the accommodation space for thicker and severely deformed sediments. A syn-sedimentary rollover and onlaps of the sediments onto the basement and F2 are particularly well visible in the thickest part of the basin between 2.8 s and 1 s (3.5 km to 0.9 km depth). From 0.5 s to 2.8 s (0.45 to 3.5 km depth), the Eyjafjarðaráll appears as a half graben, i.e., a geometry likely inherited from the basement (Figure 13(f)). The half graben is ~10 km wide at depth, but above the sedimentary unit 2 (~1.7 s), the basin is narrower with <5 km width and is relatively symmetric. These configurations indicate a migration of the axis of this basin in time. The westward dipping NNE en échelon sinistral faults in the continuation of the NE horst of Skagafjörður-Tröllaskagi clearly spread outside of the Eyjafjarðaráll Basin above 1.7 s.
Although the deformation cannot be dated absolutely, the polyphase tectonic is evident from these profiles, with reactivation of the upper Miocene-lower Pliocene faults and syn-sedimentary deformation. The most severe structural phase is speculatively prior to upper Pleistocene-Holocene as it is sealed by the reflectors at 0.6 s (~0.6 km depth), above which only some of the faults remained active, bounding the basins and highs on the seafloor (Figure 13(e) and Figure 13(f)).
Additionally, the dip values of the faults and the thickness of intrusions must be commented on. For this, a comparison is made with the outcrops of the Icelandic crust in West Iceland where both the age and the level of erosion (1 and 1.5 km) are better known than in the study area in North Iceland. Such comparison provides valuable information for the structures at similar depths on the seismic profiles. First, faults dip 80˚ - 85˚ in the sediments and 75˚ - 70˚ in the basement on the seismic profiles. These dip values, which were also noted by on all the six vintage profiles, are identical to the faults in West Iceland . These dip values exclude listric faults on the profile, especially in the absence of thick clay and salt layers. Secondly, sills and dykes are important structures as dykes contribute to spreading, and both greatly to the build-up of the crust and fluid circulation. Therefore, their interpretation at depth on seismic profiles is important. However, the window of West Iceland shows that down to 1 and 1.5 km crustal depths the thickness of sills and even repeatedly injected dykes is no more than 15 m and 8.5 m, respectively (Figure 4(d) and Figure 4(e)). Such intrusions are too thin to be detected at the same depths on the vintage seismic profiles.
5.4. Selected Earthquakes
The earthquake swarms of 1994-2016 (SIL network, IMO), along with sixteen focal mechanisms from 1994-1997 , and 2013 are interpreted in light of the regional structures (Figures 14(a)-(c)). Some clarifications must be made prior to the interpretation. The earthquake swarms of 1994-2016 are up to M 4, and generally between 2 and 10 km depth. They include three categories of events. Firstly, the majority of the earthquakes in the Eyjafjarðaráll and Öxarfjörður Basins, the HFF offshore, GOR, Grímsey, Tröllaskagi, Flatey and Theistareykir reflects the transform mechanism. Secondly, some of these earthquakes are known to be associated with fluid migration such as in GOR , and in the Theistareykir geothermal field , but they can also be induced by geothermal operations, e.g., . Thirdly, to the east of the study area (Figure 14(a) and Figure 14(b)), the 1994 to 2016 earthquakes also include a N-S lineation. This lineation reflects the post-rifting deformation related to the repeated Krafla dyke injections of 1975-1984 that led to the northward migration of earthquakes , surface faulting and ground widening . Regardless of category, the correlation of the earthquakes and the identified structures brings new insights to the tectonic deformation:
1) The earthquakes are distributed along the trace of the HFF. To the west of this fault, however, the earthquakes fall within at least two WNW left-stepping en échelon segments, which onlap the fault and are conformable with the dextral motion of this transform segment. The earthquake cluster in the WNW section of the Eyjafjarðaráll Basin is within the WNW faults bounding the basin and displays both dextral and normal-slip (Figure 14(b)), while the focal mechanisms inside the basin show only dextral strike-slip motion (Figure 14(c)). These features are evidence of polyphase tectonic where the earthquakes only show the present-day activity of the underlying WNW structures. Furthermore, the earthquake cluster is denser to the east of the NE sinistral oblique-slip shear zone of north Skagafjörður-Tröllaskagi, indicating a block compartmentalisation by that zone at depth. The earthquakes are more diffuse beyond the NE shear zone in the N-S section of Eyjafjarðaráll, where they align on the mapped NNE/ENE sinistral, as well as the WNW and NW short dextral strike-slip segments (Figures 14(a) and Figure 14(b)).
2) Most focal mechanisms show either dextral slip on WNW/NW or sinistral slip on NNE fault planes along GOR. The exception is in Nafir where the sole earthquake with a N-S normal-slip focal mechanism falls on the mapped boundary fault of the Skjálfandi Basin there (Figure 14(c)).
3) The earthquakes in GOR are organised in two overall groups. One is an isolated NW cluster around Stóragrunn and relatively confined within the plate boundary. The other is the main cluster, striking WNW and extending from Nafir to Öxarfjörður. As this latter WNW cluster clearly cuts the structures of Skjálfandi, Mánáreyjar and Öxarfjörður diagonally, it indicates that earthquakes are not always ruptures of underlying older structures (Figure 14(b)). Within the Öxarfjörður Basin, the cluster delimits a WNW rectangular block in absence of clear fault traces. This may indicate the existence of WNW blind segments of the transform zone/oblique rift, which have not broken the surface.
4) The overall earthquake cluster of GOR is interrupted in three places (Figure 14(a)). The first is a local gap located to the northwest of Nafir, and coincides with a small ENE basin on short sinistral oblique-slip faults. This gap and the basin fall on the trace of the ENE sinistral shear zone that controls the NE-striking horst to the north of Skagafjörður-Tröllaskagi. Half way between Grímsey and Mánáreyjar, the cluster is again interrupted on a NNE seismic lineation that extends southwestwards to the north of the HFF at the latitude of Flatey and is within the ENE sinistral shear zone of Area 1 (Figure 5). Focal mechanism shows possible sinistral motion on this and surrounding NNE seismic lineations . The third interruption is a narrow zone farther to the southeast where the N-S faults of Mánáreyjar are dragged to NNE by the dextral motion of the NW boundary fault of the Öxarfjörður Basin. Even a local interruption of the cluster within the Öxarfjörður Basin aligns on the continuation of one of these rotated NNE faults.
The interpretation of the 1994-2016 earthquakes only exposes the recent ruptures of the pure extensional and shear fault sets. Widespread variations in the maximum horizontal stress (SHmax) are known in both the paleostress since the upper Tertiary , and in the stress map of North Iceland deduced from earthquakes, borehole measurements and older geological indicators . These variations are responsible for normal and strike-slip faulting on a range of Tertiary to present fault sets, identical to the structures here. The earthquakes reflect mostly reactivation of underlying structures but, additionally, they could indicate the activity of a new plate boundary. It is speculated that the WNW earthquake cluster cutting the NW en échelon segments of GOR in Skjálfandi, Mánáreyjar and Öxarfjörður diagonally corresponds to either a young horsetail splay of the underlying structures, or to the initiation of a new segment at GOR (Figure 14(d)). In fact, among the three segments of the transform zone, GOR is known to be at the youngest stage of development .
6. Conclusions and Discussions
The comprehensive regional tectonic map offshore and onshore North Iceland provides a continuous coverage of the structural pattern of rift, oblique rift, transform segment and a major microplate. About 2258 fractures were identified on high resolution Multibeam, Single beam and satellite images. The structures are faults, open fractures, prominent joints, eruptive craters onshore, and marine volcanic edifices offshore, but dykes are not visible on the images. Their qualitative and quantitative multidisciplinary structural analysis, and correlation with selected data, show the range and complexity of the deformation at and between these plate boundaries. The structural pattern consists of six sets of normal and shear fractures, reflecting the effect of rifting and transform faulting. These sets have been active in a polyphase tectonic since upper Tertiary during which, they have formed the basins and horsts; compartmentalised the blocks at any scale; induced fault drags, block rotation, syn-sedimentary deformation; and controlled the volcanism.
6.1. Main Conclusions
The main results regarding the structural pattern (Figure 15(a)), the configuration of plate boundaries and a microplate (Figure 15(b)), the styles of deformation above deep-seated shear zones (Figures 15(c)-(e)), or associated with structures that have fully broken the surface (Figures 15(f)-(k)) are as follows:
・ All sets have dip-slips. The N-S faults are purely extensional, but the surface expression of deeper strike-slip zones are left- and right-stepping en échelon segments, indicating dextral motion on WNW and NW, and sinistral motion on NNE to ENE faults. Shear motion could not be identified on the E-W set.
・ Faults are segmented, ranging in total length from ~1 km to a maximum of 30 km. Individual segments of a single fault can dip in opposite directions as faults are steeply-dipping (85˚ - 70˚ in the upper crust) and some have not fully coalesced, such as the HFF.
・ The N-S normal faults form 39%, and the oblique-slip Riedel shears 61% of the fracture population of which, the NW, NNE and ENE faults are the most frequent. The WNW faults, including the HFF, are though less frequent, and the E-W fractures the least of all.
・ The study area is dominantly highs onto which narrow basins are carved. The N-S, NW, and to a lesser degree the WNW and NNE/NE sets, contribute more than the NE/ENE sets to basin and horst formation such as in the Eyjafjarðaráll, Skjálfandi, Öxarfjörður, Eyjafjörður and Theistareykir Basins, as well as the Grímsey-Flateyjarskagi and Mánáreyjar-Tjörnes Highs. The NE/ENE sets are mostly long parallel shear zones cutting the regional offshore and onshore structures diagonally.
・ Within the rift, oblique rift, and the transform zone, magma injects into the N-S normal faults, the Riedel shears, or at the intersections of these fault types. Magma injection into shear faults indicates that the oblique rift and the transform zone are leaky.
・ Marine volcanism pierces the K-R and the Mánáreyjar High. A startling feature is that these ridges evolve into basins southward, a configuration labelled as “along-axis Ridge to Basin” here. The proportion of ridge (70 km) to basin (65 km) is balanced for the K-R and the Eyjafjarðaráll Basin, but uneven for Mánáreyjar (90 km) and Theistareykir (≥25 km).
・ The range of deformation associated with shear fractures includes horsetail splays (e.g., HFF and Stórihver Fault); rhomb-graben in relay zone of strike-slip faults (e.g., HFF); and the rotation of the rift structures by dextral strike-slip motions of the HFF and by the boundary fault of the Öxarfjörður Basin.
・ The mapped faults extend to 4 - 5 km depth, cutting the volcanic and sedimentary crust of upper Miocene-present (Figure 13). Listric faults are absent and the steeply dipping faults display antithetic and synthetic geometries and occasionally form extensional flower structures on deeper strike-slips.
・ The tectonic is most intense along the HFF and Eyjafjarðaráll Basin where the reactivated basement faults deepened the basin as a half graben, creating accommodation space for thicker sediments, syn-sedimentary deformation, roll-over, and sediment onlaps. Most of the structuring is prior to the upper Pleistocene-Holocene, hence an upward decrease in the number and the dip-slip of the faults. Some of the faults are still active, bounding the basins and highs on the seafloor.
・ The American plate separates along the narrower K-R/Eyjafjarðaráll, and the Eurasian Plate across a wider zone including the Skjálfandi and Öxarfjörður Basins. Öxarfjörður is likely one of the NW en échelon segments of GOR and merges with the active NRZ. The outer plate boundaries of GOR/NRZ, DL, K-R/Eyjafjarðaráll isolate a major rhomb-shaped microplate labelled as Grímsey-Tjörnes-Dalvík. The Skjálfandi and Eyjafjörður are basins on oblique-slip faults within this microplate.
・ The earthquake cluster and focal mechanisms match the overall traces and motions of the mapped structures in the HFF and GOR, but a few gaps in the cluster coincide with the ENE shear zones or with the rotated underlying faults. The long WNW earthquake cluster in GOR cuts diagonally the older structures and could correspond either to a blind splay fault, or to the initiation of a new segment of plate boundary parallel to the HFF.
This paper presents a range of deformation resulting from a polyphase tectonic since the upper Miocene, while its last stages from Holocene to the active deformation have received more attention. Some discussions on these results and other contributions are appropriate.
・ The younger mapped structural features are in good agreement with previous studies of marine volcanoes and the Holocene N-S faults offshore and onshore North Iceland . However, the six sets of extensional and shear faults identified here are not isolated, but typical of areas where diverging plate boundaries interact such as West Iceland , South Iceland , and Southwest Iceland . The key feature in all these areas is that the Riedel shears are roughly 2/3 of the fracture population and the rift-parallel structures only 1/3.
・ Despite a hotspot, the regional partitioning of rift/ridge segments, transform/transfer zones, parallel basins, fracturing and volcanism oblique to the rift seen in North Iceland are identical to diverging plate boundaries of similar age in other continental and oceanic contexts, e.g., . Similarly, specific structural features of North Iceland are also in good agreement with combined extensional and Riedel shears pattern , en échelon and sinuous fractures above deep-seated shear zones , horsetail splays , and rhomb graben in release zones described elsewhere.
・ Hydrocarbon is the most extensively exploited natural resource associated with diverging plate boundaries, where extensional and shear fractures also control the magmatic heat source, permeability and geothermal fluid flow . Although the pockmarks and the geothermal activity are occasionally detected offshore North Iceland, further exploration of these prospects requires a regional tectonic frame and detailed structural pattern, which are provided by this present contribution. Similar tectonic already proved to control the geothermal activity in Theistareykir .
A few of the obtained results, however, may diverge from previous findings:
・ The short NNE fault segments of the Eyjafjarðaráll Basin to the north of the HFF may appear as a part of the N-S boundary faults of the basin. However, these segments cannot act as boundary faults as they all dip towards north-northwest, occupy the width of the basin and expand outside of it. These NNE structures are with typical en échelon arrangement and sinuosity, and they are sinistral oblique-slip faults above deep-seated ENE shear zones (Figure 11).
・ The en échelon NW segments have been reported as surface expression of the WNW-striking HFF onshore . Such geometry was observed along the strike-slip faults of all sets in Theistareykir , and more regionally in this study.
・ The dextral offset of the HFF onshore is estimated at 100 km or 140 km , which is cumulative since the upper Tertiary. On the map (Figure 11), the HHF cuts the onshore and offshore structures of Mánáreyjar, Tjörnes, Skjálfandi, Flateyjarskagi, Eyjafjörður and the horst and graben to the north of Skagafjörður-Tröllaskagi without significant horizontal offset. The present structural configuration is, therefore, the latest stage of the activity of this transform segment and further evidence of polyphase tectonic.
・ The moraines and glacial sediments of the last glaciation cover the Iceland shelf and present a strong morphology on images. It can be argued that some of the mapped structures on the seafloor are not of tectonic origin but morphostructures of glacial origin. Considering that the glaciation ended at 14,000 yr BP , time has been sufficient for the faulting of all sediments with the seafloor during the continuous tectonic activity associated with plate boundaries. Even a few decades of earthquake monitoring shows the high frequency of fault ruptures and their widespread strikes. In fact, the offshore “morphostructures” mentioned above are faults, which extend from the seafloor into the crust and contribute to the tectonic deformation at depth (Figure 13).
The authors thank Bryndís Brandsdóttir (Institute of Earth Sciences, University of Iceland), Guðrún Helgadóttir (Marine and Freshwater Research Institute), Hilmar Helgason (Icelandic Coast Guard), and Thórarinn S. Arnarson (National Energy Authority of Iceland-Orkustofnun) for their permissions to use the Multibeam and Single beam bathymetric data. Thórarinn S. Arnarson also permitted us to use the two seismic profiles on Figure 13. At ÍSOR our special thanks go to Ögmundur Erlendsson for providing higher resolution bathymetric maps and vintage 2D seismic profiles; Albert Albertsson and Guðrún Sigríður Jónsdóttir for their support with GIS; Árni Ragnarsson for the statistical analysis and Bjarni Richter for encouragement. The authors also thank the anonymous reviewers of the International Journal of Geosciences for constructive comments, as well as the Assistant Editor, Beverly Guo, for support during publishing process.