Unstable Rifts, a Leaky Transform Zone and a Microplate: Analogues from South Iceland

A structural analysis was undertaken in the South Iceland Seismic Zone (SISZ) transform zone, and in the Hreppar Microplate (HMP) located between the propagating Eastern Rift Zone (ERZ) and the receding Western Rift Zone (WRZ). The age of the oceanic crust in these areas is 3.4 Ma to present. About 20,000 fracture segments on aerial images reflect the dominance of NNE extensional structures in the WRZ. Around 9,000 basement faults, intrusions, secondary fractures, surface ruptures of earthquakes, and leakages were mapped in the outcrops of the HMP and the SISZ. About 23% of these fractures strike NNE, while 77% are dominantly northerly dextral and ENE sinistral, and secondarily E-W, WNW and NW sinistral strikeand oblique-slip structures, forming a Riedel shear pattern typical of a transform zone. Dyke injections into Riedel shears indicate a leaky transform zone. Fractures reactivated, accumulated slip, and re-opened for fluid flow. The ENE faults dip mostly to the southeast and could be the present boundary of the SISZ to the north. A 10 30 km wide ENE structural zone hosts a valley to the east, which could be deeper in the west. This ENE zone contains all the earthquakes, dominant ENE rivers, frequent ENE secondary fractures, and is likely the active part of the SISZ. The HMP does not show rotation since 3.4 Ma despite being between two rift segments. Future propagation/recession of the rift segments along their N55 ̊E sections would cause a migration and a clockwise rotation of the SISZ from ENE to E-W. The boundary faults of the SISZ would then be E-W, with unchanged internal Riedel shears, compensating its sinistral motion. Insights into complexities of diverging plate boundaries are critical for resource management in such tectonic contexts.


Introduction
Iceland is a surface expression of a Mid Ocean Ridge, above a hotspot, and a Tertiary member of the North Atlantic Igneous Province associated with the opening of the North Atlantic since 62 -58 Ma [1] [2]. Presently, the American and Eurasian plates separate in the direction of ~N105˚E [3] along active plate boundaries. However, divergence since Tertiary has also led to a series of extinct rift and transform segments and associated structures in Iceland (Figure 1(a)).
The study area covers portions of three structural domains, i.e., the active South Iceland Seismic Zone (SISZ), considered to be a young transform zone, the adjacent Hreppar Microplate (HMP), and the Western Rift Zone (WRZ) also called the Reykjanes-Langjökull Rift Zone (Figure 1(b)). The HMP is bounded by the SISZ, but also by the receding WRZ and the propagating Eastern Rift Zone (ERZ).
Since [4] identified that spreading is accommodated by both divergence along rift/ridges, as reflected in magnetic anomalies [5], and shear motion along transform zones, identified on the basis of earthquakes [6] [7], details have emerged regarding the complexities of these apparently simple linear boundaries. The complexities result from changes in spreading, propagation or failure of rift segments, instability of transform zones, and formation of microplates.
Initially, the processes of diverging plate boundaries in both continental and oceanic contexts were mostly reported from fast spreading centres in East Pacific, the San Andreas Fault, and the East African Rifts/Gulf of Aden. A great portion of those findings comes from offshore gravity, paleomagnetism, bathymetry, regional seismic data and analogue modelling. Therefore, observations from rare outcrops of oceanic crust in slow-spreading centres provide first hand insights into the tectonics of these plate boundaries, which ultimately control basin formation and sedimentation, block compartmentalisation, magmatism and fluid flow in fractured hydrocarbon and geothermal reservoirs, as well as related metalliferous deposits. This paper presents extensive structural observations from a 10-year geological investigation in varied regional structural domains of different characteristics and kinematics in South Iceland. The investigations were carried out for the National Power Company (Landsvirkjun), as a means to assess the seismic risk for the construction of three hydro power plants. This paper is organised as follows: 1) It first gives the key features of the deformation associated with diverging plate boundaries, as a common frame. 2) Then it presents the methods and data, followed by our observations and interpretations. 3) Finally, it discusses the significance of the interpreted deformation with respect to the regional structures, including similarities and differences with selected key tectonic models. The steps undertaken in the investigations are examples of thorough exploration methods that range from aerial images to outcrops, from cm to tens of km in scale, and from local to regional multidisciplinary interpretations. Not only does this contribution bring new insights into the geology of Iceland, but the methods and results in this case study will also benefit exploration and resource management at other diverging oceanic plate boundaries.

Geodynamics of Diverging Plate Boundaries
As oblique rifting is absent in the study area, an overview of the key features associated with rift/ridges, transform zones, and microplates is presented below.  Continental rifts and oceanic ridges trend orthogonal to spreading. From initial continental break-up, extension is accommodated by rift boundary normal faults, but intra-rift complexities appear during rift evolution with half-grabens, subsided central grabens/axial valleys, uplifted rift shoulders, rotated fault blocks and detachments [8] [9], as well as ridge segmentation [10] and ridge instability [11] [12]. During their evolution as slow (~20 mm/yr) or fast (80 -145 mm/yr) spreading centres [3], the crust thins out under stretching, leading to volcanism [1] and the build-up of the oceanic crust by magma injections from the mantle into the axes of rift/ridge segments. Acidic and basaltic magmatism appears either as eruptions in central volcanoes and from fissures, or as injections of deeper sills and dykes in and adjacent to rifts/ridges. Seismic swarms accompany rifting episodes and magma injections [13] [14] [15]. High temperature geothermal activity is present in rifts onland, and as hydrothermal vents and black smokers along ridges offshore. The geometry, the rate and direction of spreading are reflected in magnetic chrons, which are the basis for modelling of the rifts and their instabilities [11] [16].  Transform zones offset two ridge segments. Their sense of horizontal motion differs from Anderson's fault theory as transform zones are parallel to the relative plate displacement and accommodate the spreading [4] [17]. At oceanic stages, the transform zones separate the oceanic lithospheres [18], but their origin is debatable. They are seen initiating as crustal weaknesses during continental break-up [4] [19] and are continuously active past the oceanisation stage [20], appearing only after oceanic accretion has started [21] [22], predating oceanic accretion [23], or developing after ridge propagation stops [24]. The crust is neither built nor destroyed, except in leaky transform zones where volcanism produces locally new crust under oblique-extension [25] [26] or magma is injected as deep dykes and sills [27]. The depth of earthquakes in transform zones is generally 10 -15 km in the brittle crust, but limited by the 600˚C isotherm [28]. While some consider the crust to be cold in transform zones [29] [30], others predict an elevated heat flow [31] [32]. Sea water contributes to serpentinisation, and along with hydrothermal fluid flow, they change the total stress distribution [33]. High seismicity due to rupture of strike-and oblique-slip faults, and internal bookshelf tectonics with counter-clockwise block rotation on parallel strike-slips [34] [35] [36] are two key features within transform zones. However, the spreading and slip rates, the length of the transform zones, as well as the presence of a viscous layer determine the geometry and structural pattern of these shear zones [37]. Oceanic transform zones are knife cut valleys of ~10 to 30 km width in slow spreading ridges, and >100 km in fast spreading centres [9] [24], with inner and outer bends at the intersections with ridges, and offsets ranging from 50 km to >350 km [29] [38]. Narrow transform zones are controlled by parallel boundary and internal faults, while the overall shear motion of wider transform zones is accommodated mostly by secondary faults oblique to the zones [9] [37] [39]. At a young stage, the internal oblique faults present en échelon arrangements above deeper strike-slip faults [40] [41] [42]. Continental transform zones are sharp strike-slip faults with internal stepovers and splays [43]  [51], as well as small en échelon basins [52]. Some of these internal structures are common in both oceanic transform zones [53] [54] [55], and in transfer zones that connect continental rift segments [56] [57]. Transfer zones equally undergo seismicity and local magmatism [58] [59] [60]. In all of these contexts, partitioning of the damage zones of strike-slips faults [8] [61] [62] [63] [64] controls the fluid flow and permeability along individual structures.  Transform zones are generally orthogonal to rift/ridge segments but their angular relation is an important feature. When initiating from crustal weaknesses, particularly in leaky transform zones, a complex transtensional Riedel shear zone develops at an angle of ~73˚ to the rift segment, but later rotates to the stable orthogonal configuration at an oceanic stage [37] [65] [66].  Plate boundaries become unstable due to changes in the rate and direction of plate motion, and readjust through overlapping rifts, rift-jumps, migrating and rotating transform zones, as well as the formation and rotation of microplates [24] [35] [36] [67] [68]. When spreading is accommodated by overlapping rifts/ridges, especially near hotspots, the propagating segment lengthens along strike, while the active length of the receding segment shortens and it eventually fails [11] [69] [70]. Spreading and propagation cause the V-shaped geometry of the magnetic chrons and the pseudofaults on each side of the propagator, a transfer of narrow lithospheric block from one plate to another by rift-jump, and a migration of the transform zone [11] [12] [71]. Hotspots produce an elevated topography, sometimes above sea level [72], thus facilitating the propagation of rifts/ridges away from the hotspot [68] [73]. However, examples to the contrary also exist where ridges have a tendency to propagate towards hotspots and halt there [74].  Overlapping rifts generally display inward bends at the tips of propagating segments, with complexities in the linkage area [75]. Typically, in slow spreading centres where a major transform fault is absent in the linkage area, extension is accommodated by bookshelf faulting on, and a horizontal rotation of small blocks along strike-slip faults parallel to rift/ridge segments [76] [77] [78] [79]. In fast oceanic spreading, some overlap areas of propagating/receding ridges are broad "non-transform" shear zones. There, seismicity occurs on pre-existing faults of the failed rifts, which were tectonically rotated into new oblique trends [68]. In a fast spreading context, when the overlap area becomes too wide, a microplate develops by a transfer of lithosphere from one plate to another, and rotates around a vertical axis to accommodate independently the spreading [67] [80]. The amount of rotation is a function of the rates of spreading and propagation [81]. Microplate rotation produces an oblique seafloor, expressed as a fan-shaped pattern in magnetic chrons or rotated bathymetric fabrics [68]. Additional deformation appears in the rotated microplate if adjacent to a transform zone. Magnetic chrons would show seafloor rotation in the microplate if the deformation is accomplished by bookshelf faulting and internal rotation [35]. But no rotation occurs if the deformation is accomplished by transform-parallel shears [82]. The study area contains similar structural domains. Therefore, the outcrop observations and interpretations in this paper give a means to check which of the processes described above have taken place in South Iceland. As diverging plate boundaries host both the geothermal fields and the petroleum systems, the deformation at rifts, transform zones and microplates have direct implications for both resources. However, the above processes occur over geological time-scales, of the order of million years, and are not fully detectable by short-term present-day crustal deformation measurements.

Geological Settings of Iceland and the Study Area
With a spreading rate of 20 mm/yr [83], Iceland is a slow spreading centre, which displays active and extinct plate boundaries since the Tertiary ( Figure  1(a)), as well as Quaternary intraplate earthquakes [84] and volcanism [85] [86]. The active rift segments are the Northern Rift Zone (NRZ) and the Eastern Rift Zone (ERZ), as well as the Western Rift Zone (WRZ) (Figure 1(a)). The active transform zones, initially identified from earthquakes [6] [87] [88], are the Tjörnes Fracture Zone (TFZ) and the South Iceland Seismic Zone (SISZ). Oblique rifting occurs in the TFZ [89], and on the Reykjanes Peninsula [90] before joining the Reykjanes Ridge (RR) offshore.
Relocations of rifts (rift jumps) since 24 Ma [91] have left a series of parallel extinct rift segments that stretch from the centre of the hotspot towards the northwest (Figure 1(a)). Opinions diverge on many aspects of rifting and transform faulting. The Northwest Rift Zone, active from 24 to 15 Ma [92], shifted to the Snaefellsnes-Skagi Rift at 15 Ma [93] [94], and from the Snaefellsnes Rift Zone (SRZ) to the WRZ at 5/7 Ma [95]. These two latter were considered to be overlapping rifts, with the SRZ receding and the WRZ propagating southwestward until the latter reached its present geometry [93] [96]. Some consider an eastward jump of the Skagi Rift Zone to the NRZ at about 12 Ma [94] [97], or 8 -8.5 Ma [98]. However, others propose an active predecessor to the NRZ located to its east, which jumped westward at 12 Ma and at 6/7 Ma, stabilised to its present location and propagated southward to the ERZ [99] [100]. An E-W transform zone around 65˚N has been speculated during Tertiary rifting, connecting the SRZ, SKR and NRZ [93] [94] [100] [101], with a dextral [102] or a sinistral motion and an internal deformation similar to the SISZ [103]. During plate boundary reorganisations, crustal blocks have been transferred from one plate to another, resulting in the formation of the Grímsey-Tjörnes-Dalvík and the Hreppar Microplates. Many areas were subject to uplift and subsidence, unconformities, flexuring, major block rotation and intense fracturing [53] [93] [102] [104] [105]. Despite past and present overlapping rifts, the V-shaped geometry is demonstrated only offshore along the RR [71], although that geometry does not appear in some bathymetric/paleomagnetic studies along the RR [106] [107].
Rift segments consist of fissure swarms of parallel normal faults, eruptive fissures and dykes emitting basaltic and intermediate lavas, as well as calderas hosting acidic rocks [97] [108] [109]. These swarms strike N-S in the NRZ and NNE/NE in the WRZ and ERZ (Figure 1(a)). Reworked, glacio-lacustrine or local marine sediments are intercalated in the series spanning 15 Ma, while subglacial hyaloclastite ridges are within formations of <0.7 Ma age. The oceanic crust is eroded down to 1.5 km in Tertiary rift flanks, and down to 0.7 km in Hreppar [110] [111] [112]). Both of the active transform segments undergo continuous earthquakes, as large as M L 6 -7 [113] due to ruptures of strike-and oblique-slip faults [83] [114] [115] [116] [117]. High temperature geothermal fields (200˚C -300˚C) are at the intersection of rift and transform segments [118] [119] and medium to low-temperature fields (<150˚C/km) are common off rifts. Both types of resources are fracture-controlled.
In the study area, the main structural elements of the WRZ, ERZ, HMP, and SISZ are:  The late Pleistocene-Holocene lavas are tholeiitic basalt within the overlapping ERZ and WRZ, but transitional alkalic and alkalic at the southwestern tip of the ERZ (Figure 1(b)). Furthermore, the results of continuous GPS measurements indicate a higher horizontal velocity east of the ERZ, holding the American Plate fixed [120]. This is in agreement with the estimation that only up to 35% of the total spreading rate takes place along the WRZ [121]) and 65% along the ERZ (Figure 1(b)). The petrology and the spreading rate indicate a southwestward propagation of the ERZ.  The HMP, which formed at the plate boundary and shifted away, is caught between the WRZ, ERZ and SISZ, but its northern boundary is poorly defined. The trace of the suspected Tertiary transform zone at around 65˚N coincides with the northern boundary of this microplate. However, that boundary in central Iceland presently displays an active tholeiitic volcanism and fissure swarms similar to rifts (Figure 1(a)). The fracture pattern consists mostly of NNE rift-parallel faults and dykes, although ENE and WNW sets have been also identified [122]. Another structure of the HMP is an antiform, labelled as an "anticline" in the literature. Its proposed axis is mostly NE, parallel to the rifts, but bends to N-S in the south adjacent to the SISZ [122]. Such a structure does not result from shortening due to compression. It forms due to loading and extension associated with rifting [123] [124], as the lava pile tilts first towards the active rift zone (WRZ), then towards the new rift (ERZ).  The study area is covered by the 3.4 -0.7 Ma Hreppar formation, which extends from the HMP into the SISZ [96]. Those series consist of basaltic lavas and hyaloclastites, with acidic rocks limited to Stóra-Laxá (SLCV) and Þjórsárdalur central volcanoes (ÞCV) [125], and thin sills [126]. The 8,600-yrs old postglacial Þjórsá lava covers the Hreppar Formation in places and is the longest Holocene tholeiitic lava in Iceland, with some 22 m thickness and ~25 km 3 volume [127]. Superficial formations overlie both the Hreppar formations and the Þjórsá lava. Although some models predict elevated temperature and thin lithosphere at depth under transform zones [31] [32], regional tomography in South Iceland shows that the crust is thin in the WRZ, and thicker (20 -24 km) under the SISZ until its intersection with the ERZ near the centre of the hotspot [128] [129]. Temperature gradients also indicate a cold crust (50˚C -150˚C/km) in South Iceland, with the highest values being near the rift segments [130].  The SISZ, connecting the WRZ and ERZ, is proposed as an overall E-W seismic zone with a sinistral motion and a constant slip rate of 19 ± 1 mm/yr [83] [116]. Its width is ~25 km, corresponding to the extent of the present-day earthquakes, and its length is considered 80 km between Hengill and Hekla. However, the results of detailed GPS measurements indicate that the SISZ could be longer, extending to Torfajökull farther east [131]. Opinions also diverge on the overall trend of this transform zone as some suggest it to strike WNW [88] [100] [132]. Regardless of its trend, the SISZ is thought to be a nascent transform zone. In the absence of boundary faults, the sinistral motion of this zone is achieved by periodic earthquakes up to M L 7, generally shallow, but as deep as 10 km such as in the 2008 earthquake doublet [116]. Earthquakes occur primarily on parallel internal N-S dextral faults, with an anticlockwise block rotation or bookshelf tectonics [83] [133], causing wide damage zones that extend into the HMP [134]. Outcrop mapping [126], and relative locations [135] [136] [137], however, show ruptures on additional albeit less frequent strike-slip fault sets striking ENE, WNW/NW and E-W, indicating that the SISZ is a Riedel shear zone [138]. The ENE sinistral and the N-S dextral faults form the main conjugate source faults of earthquakes. The N-S faults can reach up to 20 -25 km in length and their surface expression is left-stepping en échelon sink-holes, open fractures and push-ups. They are nearly vertical at depth [139], and have a horizontal offset up to 2 m [116] and dip-slips < 1.5 m in a single earthquake [117]. The ENE sinistral faults appear as shorter segments at the surface and have a right-stepping en échelon arrangement [140] [141].

Methods and Data
The data and methods presented in this study address various scales of tectonic observations, reflecting rifting and transform faulting in formations of different ages. In the following chapters:  We first present a lineament map made from aerial photographs that covers the WRZ, HMP, and SISZ for an overview of the regional fracture pattern ( Figure 1 Figure   14).  As the results have regional significance, the deformation is further interpreted in Chapter 6 with respect to selected common processes of diverging plate boundaries. For this, a morphostructural analysis is carried out using topographic cross-sections through the rift segments, the HMP and the SISZ, to address the geometry and the boundary of the SISZ ( Figure 15). Finally, the regional lineament map, as well as the results of the latest available continuous GPS measurements, are interpreted to better understand the longand short-term deformation of the HMP ( Figure 16 and Figure 17). The lineaments cover an area of ~3500 km 2 , and cut the upper Tertiary bedrocks to late Holocene lavas in the WRZ, HMP and SISZ. As the map was prepared to provide a regional overview of the fracture pattern, the structures are drawn as undifferentiated faults, dykes, and prominent joints, irrespective of the magnitude of fault displacement, thickness of dykes, and whether a segment is major or minor.

Regional Tectonic Lineaments
The main features of the fracture distribution and frequency at the regional scale are as follows: Open Journal of Geology  The statistical analysis of the entire fracture population was carried out using a 10˚ interval in the strikes of the lineaments. Results show that the most frequent fractures strike N20˚ -N30˚E, followed by another peak at N90˚ -N100˚E, although fractures with other strikes are also widely present ( Figure  2(b)).  The visual inspection of the lineament map in the three structural domains indicates that fracturing is denser in the older Hreppar formation of the HMP to the east of the antiform, and in a few similar outcrops in the SISZ. But it is less dense in the Pleistocene-Holocene formations of the WRZ, and least dense in the SISZ where Þjórsá lava covers the older Hreppar formation.  There is also a continuity between the widely striking fractures in the WRZ and the HMP, while fractures striking northerly and ENE are more prominent in the SISZ. Although the boundary of the SISZ is more complex in detail, the Þjórsá River acts as a general limit between the HMP and the SISZ. Another feature on the lineament map is that some of the fractures continue from the SISZ into the HMP.  Whether in the HMP or in the SISZ, the Hreppar formation has been generated partly in the WRZ and partly in the ERZ, with some of the outcrops studied here still at or near the rift segments before shifting away. Therefore, the NNE rift-parallel extensional fractures are expected to dominate the tectonic pattern in these older bedrocks. Although fracture types per set will be dis-

Outcrop Geology
Four selected outcrops (Areas 1 to 4 on Figure 2(a)) were investigated for detailed information because of their relevance for the planned hydro power plants. Area 1 is Flúðir/Stóra-Laxá, while areas 2 to 4 covering Skarðsfjall/Hagafjall, Árnes/Akbraut, and Urriðafoss ( Figure 5 and Figure 6) are directly at the sites of Hvammsvirkjun, Holtavirkjun, and Urriðafossvirkjun power plants. Results of these field investigations are presented below.

Regional Strikes and Dips
The series mapped in the four areas consist of the Hreppar formation (3.4 -0.7 Ma), Open Journal of Geology      the interglacial lavas (≤0.7 Ma), and formations < 0.7 Ma. Before presenting the detailed geology, a regional overview of the series with respect to some major structural elements is necessary. Area 1 has been primarily investigated for additional regional dips of the series and some cone-sheets (dip 20˚ to 50˚) around the SLCV. All measured dips in the four areas are compiled on (Figure 3 The third series strike northerly and dip to the west. The dips are as high as 35˚ at Miðfell in Flúðir (Figure 4(a)), and decrease to the north (10˚ -18˚) outside of the SLCV and to the east (4˚ -18˚) within the volcano (Figure 3(a)). An exact relation between the northerly, ENE and WNW striking series cannot be established as the mapped area lacks critical exposure in places, and the ages of the rocks are not always well known. Nevertheless, two major observations can be made (Figure 3(b)):  Based on the successions of the series and the topography, the ENE striking series appears older than the series striking WNW. Furthermore, the series striking WNW cover the ENE striking formations unconformably in a gully to the east of Minni-Núpur. The WNW striking series are higher in altitude than the northerly striking formations filling the SLCV. They are thus likely younger. However, the relation between the northerly and the ENE striking series cannot be deduced from these data, nor the relation of the WNW striking series in Skarðsfjall and the ENE series on each side of the Þjórsá River.  Previous interpretations of the axis of the Hreppar antiform and the extent of the SLCV [145] are also reported on Figure 3. If the axis of the antiform is NE, with a bend to N-S near the SLCV, the series on each side of the antiform should strike NE and N-S with dips away from the axis. However, the mapped dips do not reflect such an axis as the WNW striking series dip to the northeast on both sides of the supposed axis, and the NE striking series towards northwest in the eastern flank of this antiform. The only anomaly in the area near the suggested NE axis of the antiform by [145] is the contact between the northerly and the WNW striking series. That contact is located at some 3 to 5 km to the west of the proposed axis, and is a structural unconformity, coinciding partly with the NNE course of the Stóra-Laxá River in the same locality. Furthermore, the elongation of the Langholtsfjall hill above the flat land, as well as the dips of the series, point to a wider and more circular contour of the SLCV. The structural unconformity comes to a halt to the north of this volcano (Figure 3(b)).

The Stratigraphy of Areas 2 to 4
Areas 2 to 4 were subject of detailed stratigraphic and structural mapping. Area 2 is the largest and on a higher topographic level. It incorporates all stratigraphic units from Hreppar formation to the series < 0.7 Ma ( Figure 5). Areas 3 ( Figure   6(a)) and 4 ( Figure 6(b)) are in the low-land within the SISZ and are covered for the most by late glacial and Holocene Þjórsá lava although some units of the Hreppar formation crop out locally. The mapped series in the three areas are as follows:  Hreppar formation has limited outcrops of ignimbrites and fine-grained lavas in the eastern part of the SLCV ( Figure 5). Otherwise, the series consist of accumulated lavas and intercalated reworked sediments, primary hyaloclastites, tillites and fluvio-glacial sediments (Figures 4(d)-(i)), to the south of HMP and in Skarðsfjall within the SISZ. Tholeiite or olivine tholeiite lavas range from fine to coarse-grained, with or without vesicles, and are occasionally plagioclase-porphyritic. Individual flows can be columnar-jointed ( Figure  4(d)), massif up to 10 m (Figure 4(e)), or as thin as 0.5 m. To the northeast of Ásólfsstaðir, the sheeted olivine tholeiite lavas of a shield volcano seem to be filling a N-S paleo-valley (Figure 4(f)), while in Urriðafoss, a scoria cone indicates a possible nearby eruptive site. The tillites and reworked sediments are fine to coarse-grained (Figure 4(i)), reaching up to tens of metres. Primary hyaloclastites are thick or thin layers, or mountains with fine-grained tuffs in the lower part and pillow lavas at their tops (Figure 4(g) and Figure  4(h)). The most geothermally altered area is the SLCV where intense intrusions and gabbro are present. However, zeolites are present throughout the series and reveal the level of erosion. While the regional alteration may be down to the stilbite zone, it is mostly within the chabazite/thomsonite zone, indicating an erosional level down to 500 -700 m in Areas 1 and 2.  (Figure 4(k)). The age of these valley infilling sediments is unknown, but considering their positions in the outcrop, they could either be the upper part of the Hreppar formation, or contemporaneous to the interglacial lavas.  The Holocene formations consist of sediments and the Þjórsá lava. The sedimentary formations are eolian, fluvio-glacial, or silt and end-moraines (Figure 4(q)), 0.5 to 15 m thick. In places, the Holocene soil contains distinct Hekla tephra layers. The 8,600 yrs Þjórsá lava covers vast surfaces in the lowland, particularly within the SISZ (Figure 4(r)). The lava is of a porphyritic, aa-type with a rough surface (Figure 4(s)), and columnar-jointed at the basis (Figure 4(t)). It reaches locally up to 20 m in outcrops. Numerous basement faults, dykes, sills, and secondary fractures, as well as surface ruptures of earthquakes and seepages of water were also mapped along with the stratigraphy in the bedrock and Holocene formations of Areas 2 to 4. They are detailed below.

Major and Secondary Fractures in the Bedrock
Note that on all rose diagrams, the letter (N) indicates the total number of mea-

Faults, Dykes and Striated Planes
In Areas 2 to 4, basement faults are those cutting the Hreppar formation and the interglacial lavas. Dykes and a few cone-sheets are limited to the Hreppar formation under the interglacial lavas, while striae were found along few faults and dykes in the basement rocks of the Hreppar formation (Figures 7(a) 5-7), and the thickness of those shown is exaggerated. Dykes dip 50˚ to almost 90˚ (Figures 8(d)-(i)), sills 20˚ -0˚ (Figures 8(i)-(l)), and cone-sheets 50˚ -20˚ (Figure 8(m)). Sills could be fed from minor cone-sheets or dykes, or be conspicuous and fed by major dykes. All intrusions are columnar jointed, with black or brown chilled margins of 1 to a maximum of 10 cm thickness (Figure 8(k)). Despite most intrusions are in the proximity of the SLCV and few even inside the volcano, hand specimens showed no acidic intrusion. They are all basaltic, similar to the tholeiite and olivine tholeiite lavas. Dykes can be dense and fine-grained, or coarse-grained and vesicular.
Their geometries range from straight to undulating, with en échelon segmentations, bifurcating, and ending upwards or downwards. They shape equally the landscape and the waterfalls (Figure 8(d)). The dominant strike of the dykes is northerly, then NNE, ENE, and the least E-W. Similar to faults, NNE rift-parallel dykes are not dominant, and many N-S, ENE and WNW/NW dykes are injected into faults of the same directions. Dyke density is unusually high in Skarðsfjall where most of the WNW and NW dykes are concentrated (Figure 7

Mineral Veins, Joints and Breccia
Measured secondary fractures are discussed separately below.  About 1528 mineral veins were measured mostly in the lower and middle parts of the Hreppar formation, with a few directly in faults zones and along dykes (Figure 9(a) and Figure 9(b)). The veins are longer (several metres) and thicker (up to ~10 cm) near the SLCV (Figure 9(c)) and the ÞCV (Figure 9(d)), and are filled with zeolites. They cut the lavas, reworked sediments and tuffs, and are generally short (up to a metre) and thin (<2 cm).
Hand specimens indicate that their infilling is zeolites such as thomsonite and chabazite (Figure 9(c)), clay, iron oxide, and even calcite and opal (Figures 9(e)-(g)). The strike of the veins varies from outcrop to outcrop, reflecting the adjacent and often eroded major faults or dykes to which they are associated. However, their cumulative strikes (Figure 9(a)) show the same sets as in other fracture types but with a different frequency. Unlike faults, dykes and striated planes, the mineral veins are dominantly ENE and NNE, then northerly and WNW. No clear relation appears between the strikes of the veins and their infilling, but the thickest veins are in the ignimbrites of Hamraholt and in the Hreppar formation between Minni-Núpur and Hagafjall , and they strike respectively, WNW and ENE.  Some 6314 tectonic joints and secondary fractures were measured in Areas 2 to 4 (Figure 10(a) and Figure 10(b)). Of these, 5119 joints are in lavas, sediments, and hyaloclastites of the Hreppar formation, as well as in fault zones cutting those series, while 1195 of the joints are in Holocene formations. The 5119 measured joints at individual sites in the Hreppar formation reflect the strikes of faults and dykes to which they are associated. Their combined strikes in the Hreppar formation are dominantly northerly and NNE/NE in the fault zones cutting those formations, but northerly in Holocene formations. However, the cumulative strikes of joints in all formations show that contrary to faults, dykes, and striated planes, tectonic joints strike dominantly ENE, then northerly, and finally NNE and WNW. These dominant regional strikes are, therefore, more comparable to mineral veins (Figure 9(a) and Figure 10). The joints have a variety of geometry and spacing. They can be regularly spaced along intersecting sets (Figure 11(a) and Figure 11(b)). Along single sets, they can be tightly parallel and separated by unbroken rocks (Figure 11(c)), or densely populated adjacent major fractures ( Figure  11(d)). Regardless of strikes, joints display extension and opening, generally a few millimetres to centimetres in the basement rocks (Figures 11(a)-(d)), but exceptionally up to 0.7 m, particularly in cliffs above the slopes when a part of the opening could be gravitational (Figure 11(e)).  Fault breccia was found in a few fault zones, along the fault trace, or on fault planes themselves (Figures 11(f)-(h)). Joints in fault breccia are generally parallel in case of pure normal faults, but complex in case of strike-and oblique-slip faults, or when fractures of various sets intersect. Riedel shears are the main secondary fractures resulting from strike-slip motions. Their inte-ractions with the main faults can cut the blocks into wedges (Figure 11(i)), otherwise they appear as various sets of secondary fractures associated with the main strike-slip structure (Figure 11(j) and Figure 11(k)). Regardless of sets, the width of the fault breccia ranges between 0.5 and 5 m, with the highest thickness associated with strike-and oblique-slip faults of the basement, particularly if the faults have been reactivated. The matrix of breccia along major faults could be clay, oxidised rock, with or without the presence of zeolites, silica and calcite. But in case of minor structures, the breccia is rather homogeneous, lacking matrix.

Holocene Structures and Leakages
In addition to the 1195 joints in the Holocene formations (Figure 10(a)), about 317 segments of Holocene surface ruptures, and 181 points of leakages (cold and warm water) were mapped in Areas 2 to 4 (Figures 12(a)-(c)). The surface ruptures and leakages are primarily in the low land of the SISZ although some were also found at higher altitude such as in Skarðsfjall or in Minni-Núpur.  (Figures 13(g)-(j)). They cut the soil, moss, scree, aeolian sediments, semi consolidated sand, the Þjórsá lava, and the Hreppar formation. Sink-holes are circular, oval or elongated, and their collapsed geometry is conical at depth. Despite earthquakes with M L 6 -7 in the SISZ, push-ups are rare in the study area. Generally, it is difficult to identify the surface ruptures in the Þjórsá lava due to subtle similarities with lava structures. Combined analysis of aerial photographs and outcrop mapping indicates that the mapped surface ruptures, particularly in Area 4 (Urriðafoss), are of tectonic origin and not lava tubes, degassing craters, river channels, or related to the Þjórsá lava front. Other Holocene structures are minor (horseshoe) landslides aligned above/near deeper source faults of earthquakes ( Figure  13(k)), or major landslides occurring after earthquakes such as in 2000 ( Figure 13(l)). Surface ruptures display evidence of extension and strike-slip, and are surface expressions of deeper shear faults. They are organised in individual segments of a few to hundreds of metres long, separated by gaps. Segments can extend over 20 km, indicating the maximum length of the underlying source faults. Extension appears as local openings of a few centimetres to 0.5 m and rarely up to 4 m, while normal-slips reach a maximum of 0.5 -1 m. In the lack of marker horizons, the shear motions are deduced from the en échelon geometry of the segments, where a left-stepping array indicates dextral motion, and a right-stepping array sinistral motion ( Figure  12). The main conjugate source faults are the N-S dextral and ENE sinistral strike-and oblique-slip faults, although surface ruptures also align on WNW, NW, and to a lesser degree on short E-W segments, all with sinistral motion in general. Surface ruptures also appear along NNE rift-parallel segments, without being frequent and having coherent en échelon arrangements. The   . Interpretation of the observed structures with respect to regional tectonic elements. (a) and (b) respectively, models of a propagating/receding rifts (modified from [157]), and the rotation of a microplate between two overlapping rifts (modified from [67]); (c) Selected individual fracture segments on aerial images are highlighted to check for their continuity. They do not show a rotation in their strikes within the microplate; (d) Horizonal velocities due to plate spreading obtained by continuous GPS measurements (modified from [120]) also show no rotation in the short term in the HMP. five sets with strike-slip motions form a Riedel shear pattern. All these fracture sets also appear in the present-day seismicity of the SISZ [135] [137]. The mapped structures in the study area are identical to those of the 2008 earthquake doublets in Hengill [117], and to the 2000 and historic earthquakes mapped in this part of the SISZ by [140] [141]. For comparison and complement, the faults mapped by [141] are also reported on Figure 12(a) and Figure 12(b), as they are in continuation of those mapped in this study.  A number of cold springs, small ponds, wet areas, and local low temperature geothermal areas and warm springs were mapped (Figure 12(a) and Figure  12(b)). Wet areas covered by green moss contrasting the grass have the lowest discharge (Figure 13(m)). The springs burst through Holocene sediments or soil, have low to moderate discharges and are accompanied by swamp and algae. Ponds are often seepages through underlying young fractures ( Figure  13(n)), or are sometimes sink-holes filled with water. In basement rocks, seepages are either through single fracture on or near fault zones, with very low discharges, white deposition and green moss. Or they are at fault intersections, especially when one fault is a strike-slip ( Figure 13(o)). Low-temperature geothermal activity was found only in the three localities of Laugar, Lón, and the island of Árnes, with very low discharges, associated or not with cold leakages (Figure 12(a) and Figure 12(b)). The springs of Laugar are 16.4˚C to 47.2˚C, in an area where drilled wells yielded up to 63.5˚C hot water. In Lón, the warm springs are 10.1˚C and 24.5˚C, and they appear between the aligned rows of cold-water leakages. Due to a low water level at the time of mapping, active warm springs were not found in Árnes. However, a large zone of geothermal alteration with present-day mineralisation was mapped ( Figure 13(p)), near a borehole yielding 50˚C hot water. The three local low temperature manifestations are aligned ENE in Laugar, WNW and NW in Árnes, and northerly in Lón.  The main structural feature of the Holocene joints, surface ruptures and the leakages is their identical trends. As their respective rose diagrams show, they strike dominantly northerly then ENE, and secondarily NNE, ENE, E-W and WNW ( Figure 10 and Figure 12(c)). The similarity in the strikes of the three structural features reflects the latest status of tectonic deformation. It indicates that the joints likely result from the activity of the source faults of earthquakes during Holocene, and that all the dominant and secondary shear fractures are permeable host fractures channelling cold and warm water (Figure 12(c)).

Interpretation
A total of 8942 faults, dykes, mineral veins, tectonic joints, surface ruptures, leakages of cold and warm water, and 110 striated planes were studied in the outcrops of Areas 2 to 4 (Figure 14(a)). They cut the 500 -700 m eroded Hreppar formation, the interglacial lavas, the Holocene and superficial formations of the HMP and the SISZ. Due to their high number, these observations represent the tectonics of the upper oceanic crust, with insights from the scale of individual fractures to that of diverging plate boundaries. For interpretation, all major structures are compiled onto a single map and their traces highlighted regionally (Figure 14(a)). To provide a deeper understanding, the historic surface ruptures mapped by [140] [141], some active faults deduced from micro-earthquakes [113] [147], and a few local magnetic structures [148] are also reported on Figure 14(a) and Figure 14(b). To reiterate, "older" refers to structures with higher dip-slips and sharper traces, indicative of a longer life-span, and "younger" to fractures cutting the Holocene formations such as surface ruptures, and structures with finer traces and less or no dip-slips.

The Structural Pattern
While the NNE fractures appear dominant on the tectonic lineament map where almost half of those regional structures is within the WRZ reflect rifting ( Figure   2), the rose diagrams of fracture types from outcrops show the presence of six sets in the SISZ and the HMP. Basement faults, dykes, joints in Holocene formations, surface ruptures and leakages strike dominantly northerly, then ENE, and to a lesser degree NNE, E-W, WNW and NW (Figure 14(b)). However, mineral veins and tectonic joints, which are generally associated with basement faults and dykes, strike dominantly ENE.
The dextral and sinistral motions of the mapped older and younger N-S and ENE structures, identified from their en échelon geometry and striae, are identical to those of the source faults of historic [141] and the 2000 earthquakes [149], as well as to the focal mechanisms in the SISZ (e.g. [116] [137]). However, the present study identifies shorter E-W, WNW, and NW shear fractures as additional sets, all of which have sinistral motion. Although micro-earthquakes align primarily on the active N-S faults in the study area, some do also align on WNW/NW structures, supporting the existence of these shear fractures ( Figure   14). The five sets with strike-slip motions form a Riedel shear pattern where the main N-S dextral and ENE sinistral structures act as (R') and (R), respectively, the short NNE subtle structures are tension fractures or (T), and the E-W to NW structures would be PS, P and X (Figure 14(c)). Such a pattern is typical of a shear zone, which in this case is the SISZ transform zone. However, the overall mapped structural pattern is a combination of the Riedel shears and the NNE rift-parallel extensional fractures where the shear fractures dominate in the SISZ and this part of the HMP (Figure 14(b) and Figure 14(c)).
The observed strikes and motions fit with a maximum principal stress N45˚E to N60˚E [138], in agreement with a S Hmax ranging between N29˚E and N62˚E in South Iceland as deduced from borehole, earthquake and geological data [150].
The attempt to interpret all striated planes of the present study, including the opposite motions in the basement rocks, the complex cross-cuttings, and the rare reverse-slip slickensides (Figure 7 Open Journal of Geology same as in the study of [151], who interpreted the complex fault-slips with several stress fields and compression parallel and perpendicular to rift axis. The present study proposes kinematic alternatives for the complex motions, which are mostly located on planes between Minni-Núpur and the SLCV (Figure 7(a)).
Some of the opposite strike-slip motions along the same sets may be due to block adjustments related to the collapse of the volcanic caldera. The normal-and strike-slip striae cutting each other alternatively could indicate simultaneous motions. Finally, the few reverse-slip striae mainly on dykes' edges could result from friction between fracture planes during repeated magma injections. Reykjanes [143], and even during the 2014 Bárðarbunga fissure eruption north of Vatnajökull [153]. While the damage zones of historic earthquakes in South Iceland show reactivation of the same source faults over a period of 100 years [134] [154], the present study demonstrates that fracture reactivation is an ongoing process over longer geological time scales.

A Major ENE Structural Zone
In addition to local tectonic and magmatic behaviour of fractures in the oceanic crust, the present study also provides insights into regional deformation.
One prominent aspect is the relevant features of the ENE fractures. As mentioned above, ENE surface ruptures group into short individual segments, but the segments correlate well together over distances, indicating much longer ENE structures within the SISZ. Young ENE sinistral segments of surface ruptures that fall in the continuation of older ENE faults (and dykes) of the HMP are tightly parallel on both sides of the Þjórsá River, forming an ENE structural weak zone with at least 5 km width. A few ENE segments are outside this zone, farther to the south (Figure 14(b)). This major ENE structural zone is at least 45 km long, i.e., the extent of the land studied here, but it could be longer. Within this zone, the great majority of the older ENE strike-and oblique-slip sinistral faults dip to the southeast, similar to shorter E-W, WNW and NW faults ( Figure  14(b)). The Þjórsá River flows within this structural zone, and in many places its course is controlled by the individual young and older ENE mapped fractures. The fact that the prominent strike of mineral veins, tectonic joints and striated planes associated with basement fractures is ENE does reflect the strong presence of this set. Local surveys have also confirmed the presence of ENE magnetic anomalies coinciding with mapped ENE faults and dykes of Area 4 at Urriðafoss [148].
The ENE surface ruptures have been considered to be a secondary set, conjugate to the main N-S dextral source faults of earthquakes. However, the existence of this major ENE structural zone with most of the older ENE faults dipping to the southeast is symptomatic of a more severe deformation within the SISZ. The ENE fractures can indeed be deeper, rupturing locally during earthquakes, but they play an equal role to the N-S faults in accommodating the overall sinistral motion of the SISZ. This major ENE structural zone is further discussed in the following chapter as it has more regional implications.

Discussion
As the observed structures reflect the deformation related to rifting, transform faulting and a microplate, their significance needs to be addressed in relation to these types of regional structures. Of the many tectonic aspects, we only interpret three features, briefly, within the space of this paper. These are the geometry and boundary of the SISZ transform zone, the rotation of the HMP, and the evolution of rift and transform segments in time. Known models are used as a basis for discussions despite being mostly from fast spreading centres.

The Boundary and Trend of the SISZ
Although a WNW trend has been suggested for the SISZ [88] [100] [132], its E-W trend is mostly proposed based on the extent of the active and historic N-S faults [83], and the results of continuous GPS measurements [131]. As oceanic transform zones of slow spreading centres are characterised by narrow knife-cut valleys of ~10 to 30 km width [9] [24], and controlled by parallel boundary faults and internal faults oblique to the shear zone [9] [37], three cross-sections are made into the study area to check for a similar valley (Figure 15(a) and Figure  15(b)). The cross-section A-A' shows that the HMP rises in altitude gradually from the WRZ towards the ERZ, and that a subtle valley exists in the eastward continuation of the SISZ where the Þjórsá River flows (Figure 15(b)). The cross-sections B-B' and C-C' cut the SISZ obliquely. Although a sharp narrow valley on section C-C´ coincides with the easternmost earthquakes of the SISZ, no valley appears along section B-B' to the west where the earthquakes of the transform zone are denser (Figure 15(b)).
An alternative is proposed here for the geometry and strike of the SISZ based on the ENE sinistral structural weak zone (Figure 14(b)), the extent of the N-S faults and the earthquakes, as well as the morphostructural elements ( Figure   15(c)). The narrow ENE valley between east of Búrfell and Árnes on sections A-A' and C-C' is partly carved by a glacier and filled with Holocene Þjórsá lava ( Figure 5 and Figure 6). However, the origin of this valley is likely tectonic and its northern boundary controlled by the ENE sinistral oblique-slip faults dipping mostly to the southeast that are identified in this study (Figure 15(c)). The dense earthquakes between Árnes and west of Hestfjall appear on the map in the westward continuation of the same ENE valley. However, a relatively flat topography to the west may indicate a deeper ENE transform zone there, where a valley could be present at depth and fully covered by superficial formations (Figure 1(b)). The ENE boundary of the SISZ is likely irregular and could consist of major right-stepping ENE sinistral oblique-slip segments, cut locally by the active N-S dextral faults. Judging from the extent of the earthquakes, the minimum width of the ENE transform zone is 10 to 15 km. But if the southernmost isolated earthquakes around Hrútsvatn belong to the SISZ, then the width of the ENE transform zone is up to 30 km (Figure 15(c)). The Þjórsá and the Ytri-Rangá are the main rivers within this ENE shear zone, both of which have a

Structural Fabric and Rotation of the HMP
The location of the HMP between the WRZ and the ERZ calls for an assessment of the deformation in the presence of two rift segments and their linkage area. In case of propagating/receding rifts in fast spreading centres [11] [12] [68], the main structural features are the V-shaped seafloor fabric, the pseudofaults, a transferred lithosphere, and an orthogonal transform zone (Figure 16(a)).
However, the V-shaped geometry is also reported from oceanic slow spreading such as the RR [71]. Although these features appear also in presence of two overlapping rifts that propagate in opposite directions, generally, the deformation in the linkage area is either through bookshelf tectonics [78], or anticlockwise rotations about vertical axes both inside the rigid microplate and at its edge near the overlapping rifts [67] (Figure 16(b)). Such a rotation differs from the bookshelf tectonics within transform zones that occur along the main internal faults of the shear zone.
The above structural features should also be present in the HMP. Although these features develop over millions of years, we assess both the older ( Figure   16(c)) and the present-day deformation of the study area ( Figure 16(d)). The regional tectonic lineament is a key as it cuts the 3.4 Ma to present series covering the axis and both sides of the Hreppar antiform as well as the SISZ [96].  earthquakes [120]. Although the measurements cover a very short time-span and focus mostly on the SISZ and the ERZ, the velocities do not show any sensible rotation, neither within nor by the edge of the HMP and the SISZ.  [36]. In slow spreading centres, transform zones may develop oblique to rifts [156], and later rotate to a stable orthogonal configuration [65]. However, particularly in case of leaky or weak transform zones, these ductile shear zones may be orthogonal at depth but oblique to the rifts near the surface [66]. Regardless of the process, their obliquity and rotation lead to complex fracturing [37].

Migration and Rotation of Unstable Rifts and the Leaky SISZ
Similarly, an assessment of the angular relation of the SISZ and the two rift segments, and the instability of these boundaries in time is necessary. The following interpretations are speculative, assuming that the spreading rate remains constant in the future as in the past 3 Ma, with 35% across the WRZ and 65% across the ERZ, and that both rift segments continue to propagate and recede towards the southwest. However, an important parameter would be the strike of the rift segments, which is a less discussed criterion as a possible contributor to plate boundary instabilities.
Based on the distribution of the tholeiitic lavas in the fissure swarms, the presently southwestward propagating ERZ and receding WRZ trend N43˚E and N38˚E, respectively (Figure 17(a)). Our results suggest that the active part of the SISZ connecting these rift segments strikes ~N75˚E, and the transform zone has a maximum width of 30 km. The high number of dykes striking similar to the shear fractures, particularly in the Hreppar formation of Skarðsfjall in the SISZ ( Figure 5 and Figure 6 and Figure 14) indicates that the transform zone is likely leaky and unstable. Even if at an unstable stage the SISZ may be E-W at depth, its present near-surface expression is a strong ENE oblique-slip sinistral shear zone at an angle to the rifts. The boundaries of the transform zone are controlled by the ENE parallel faults with its internal fabric consisting of oblique Riedel shears (Figure 17(b)). Figure 17(c) suggests the future evolution of the plate boundaries in the southern half of Iceland, which would occur over hundreds of thousands or millions years, thus undetectable by present-day measurements. That evolution would be dominated by those sections of the two rift segments that strike N55˚E. In fact, the N55˚E sections are already in place between Hengill and the Reykjanes Ridge (RR) in the WRZ, as well as from the area of the transitional crust in the ERZ to the Westmann Islands. As a consequence of the activity of the N55˚E sections of the rifts in the future, the direction of propagation and recession would turn slightly westward . To connect the active sections of the rifts, the transform zone would migrate, as it has done in the past [138], but also rotate clockwise. The 15˚ -17˚ clockwise rotation of the active rift segments would be the main drive for the 15˚ clockwise rotation of the SISZ from N75˚E to N90˚E (Figure 17(c)). Such a rotation applies to the entire transform zone and differs from its internal bookshelf tectonics. After this rotation, the transform zone sta-bilises more, with parallel E-W boundary faults but unchanged internal Riedel shears, all compensating its sinistral motion (Figure 17(d)). As only one rift segment can propagate away from a hotspot, ultimately, the ERZ would connect to the RR via this E-W transform zone, and the WRZ would cease entirely.

Conclusions
A large number of fractures were mapped from aerial images in the SISZ, the HMP and the WRZ, along with the stratigraphy and dips of the series as well as various fracture types in the outcrops of the HMP and the SISZ. The mapped fractures from aerial images are 20,127 undifferentiated segments in the three structural domains. The outcrop data from the SISZ and the HMP are 262 basement faults, 340 dykes and cone-sheets, 110 striated planes, 1528 mineral veins, 6314 tectonic joints, 317 surface ruptures of earthquakes, and 181 present-day leakages from fractures. The structures cut series of 3.4 Ma to present, are of various relative ages, and represent the tectonics of the eroded upper oceanic crust in this slow spreading site. The main findings are:  The regional dips show a complex antiform in the HMP, whose axis is some 5 km farther to the west than previously thought. Almost half of the fractures on the regional lineament map are in the WRZ, reflecting the dominance of the NNE rift-parallel extensional fractures. But this set represents only 23% of the total fracture population in the outcrops of the SISZ and the HMP where fractures of all types group into five sets with strike-and oblique-slip motions. Both within the active transform zone and in the microplate, these sets are dominantly N-S dextral and ENE sinistral, and secondarily WNW, NW and E-W sinistral, forming a Riedel shear pattern typical of a major shear zone. Striated planes generally support the fault motions. However, the few incoherent fault slip-data can be explained as a result of block adjustments near the SLCV, or as friction between fracture planes related to dyke injection.  Fractures are reactivated over time, as evidenced by: 1) The extent of the young surface ruptures from the SISZ into the HMP where they align with older faults and dykes; 2) Present-day leakages of water from the Riedel shears; 3) Cumulative dip-slips, which are higher in the older HMP (35 m or more) and lower in the SISZ (0.5 -1 m). Furthermore, dykes of the SISZ and the HMP align dominantly along the Riedel shears, indicating that magma injections into shear fractures are not a rare occurrence, and that the transform zone is leaky.  The ENE fractures are tightly parallel on both sides of the Þjórsá River and most of them dip to the southeast. They control the boundary of a narrow ENE valley that is typical of transform zones in slow spreading centres. This valley is sharp in the northeast of SISZ where it hosts the last earthquakes of the transform zone. But it does not appear to the west where earthquakes are denser, likely because the ENE transform zone is deeper there and a valley is buried under younger formations. The structural pattern, the extent of the earthquakes, and the dominant ENE rivers indicate a 10 to 30 km width for the SISZ, with an overall ENE strike highly oblique to the rifts. Its ENE sinistral oblique-slip boundaries are irregular and cut by the active N-S dextral faults.  Although the HMP is between this ENE transform zone and two rift segments, the structural fabric of the microplate does not show evidence of rotation since 3.4 Ma. However, the future propagation of the ERZ and recession of the WRZ would be along their N55˚E southern sections, causing a migration of the SISZ and its clockwise rotation from ENE to E-W. At that stage, E-W boundary faults would control the SISZ, with unchanged internal Riedel shears, all compensating its sinistral motion. As only one rift can propagate away from a hotspot, the ERZ will ultimately connect to the RR via this E-W transform zone, and the WRZ would cease entirely. Results of this research contribute to a better understanding of the geology of Iceland but also to the complexities of unstable rifts, transform zones, and microplates in slow-spreading centres in general. Such complexities control the geological processes of these varied structural domains at diverging plate boundaries where many geothermal fields and petroleum systems develop, particularly in case of fractured reservoirs. Therefore, case studies that provide observations from cm to tens of km, a multidisciplinary approach, and interpretations in time and space such as in the present paper, are critical knowledge for exploration and the choice of the targets, as well as for decisions concerning the management of these resources.