Surface ramps in normal fault zones of the Iceland plate boundary have been described in many studies, but their structure and evolution are not well understood. We show that surface ramps are manifestations of large tilted blocks (TBs) formed in dip relays of normal faults. Based on existing modeling studies, we propose three classes of TBs defined by kinematics and location of the hinge of the TB. TBs are considered a member of the family of fault relay structures that form near the surface, commonly, but not exclusively, in columnar basalts with orthotropic strength. We present high-resolution aerial vehicle–based observations of a representative set of normal faults in Iceland and compare these with geometric models we derived from modeling studies. We predict extensive tectonic cave (fluid conduit) systems under the TB, which interact with magma and groundwater flow. The general fault structure is dominated by large, subvertical open fractures reactivating cooling joints that are locally filled by basalt rubble. We propose the existence of a hybrid failure zone at larger depths before the effective vertical stress is sufficient to initiate shear fractures in intact basalt.

Close to surface, cohesive rocks fail in extension, which results in open fractures that can be several tens of meters wide, so-called massively dilatant faults. These open fractures make fault slip analysis in rifts challenging, as kinematic markers are absent. Faults in rifts often have oblique slip kinematics; however, how the amount of obliquity is expressed in the surface structure of massively dilatant faults remains enigmatic. Furthermore, the structures of oblique dilatant faults at depth is largely unconstrained. To understand the subsurface structures we need to understand how different obliquities of slip influence the surface structures and the corresponding structures at depth. We present analog models of oblique massively dilatant faults using different cohesive materials in a sandbox with adjustable basement fault slip obliquity from 0 • to 90 •. Experiments with different mean stress and material cohesion were run. Using photogrammetric 3D models, we document the final stage of the experiments and investigate selected faults by excavation. We show that fault geometry and dilatancy changes systematically with angle of obliquity. Connected open fractures occur along the entire fault to a depth of 6-8 cm, and as isolated patches down to the base of the experiments. Using the scaling relationship of our models implies that transition from mode-1 to shear fracturing occurs at depths of 250-450 m in nature. Our experiments show the failure mode transition is a complex zone and open voids may still exist at depths of at least 1 km. We apply our results to the dilatant faults in Iceland. We show that the relationship between angle of obliquity and average graben width determined on faults on Iceland matches experimental results. Similarly, fracture orientation with respect to fault obliquity as observed on Iceland and in our experiments is quantitatively comparable. Our results allow evaluation of the structure of massively dilatant faults at depth, where these are not accessible for direct study. Our finding of a complex failure mode transition zone has consequences for our understanding of fracture formation, but also influences our interpretation of fluid flow in rift systems such as magma ascent or flux of hydrothermal waters.

Abstract. Normal faults in basalts develop massive dilatancy in the upper few hundred meters below the Earth's surface with corresponding interactions with groundwater and lava flow. These massively dilatant faults (MDFs) are widespread in Iceland and the East African Rift, but the details of their geometry are not well documented, despite their importance for fluid flow in the subsurface, geohazard assessment and geothermal energy. We present a large set of digital elevation models (DEMs) of the surface geometries of MDFs with 5–15 cm resolution, acquired along the Icelandic rift zone using unmanned aerial vehicles (UAVs). Our data present a representative set of outcrops of MDFs in Iceland, formed in basaltic sequences linked to the mid-ocean ridge. UAVs provide a much higher resolution than aerial/satellite imagery and a much better overview than ground-based fieldwork, bridging the gap between outcrop-scale observations and remote sensing. We acquired photosets of overlapping images along about 20 km of MDFs and processed these using photogrammetry to create high-resolution DEMs and orthorectified images. We use this dataset to map the faults and their damage zones to measure length, opening width and vertical offset of the faults and identify surface tilt in the damage zones. Ground truthing of the data was done by field observations. Mapped vertical offsets show typical trends of normal fault growth by segment coalescence. However, opening widths in map view show variations at much higher frequency, caused by segmentation, collapsed relays and tilted blocks. These effects commonly cause a higher-than-expected ratio of vertical offset and opening width for a steep normal fault at depth. Based on field observations and the relationships of opening width and vertical offset, we define three endmember morphologies of MDFs: (i) dilatant faults with opening width and vertical offset, (ii) tilted blocks (TBs) and (iii) opening-mode (mode I) fissures. Field observation of normal faults without visible opening invariably shows that these have an opening filled with recent sediment. TB-dominated normal faults tend to have the largest ratio of opening width and vertical offset. Fissures have opening widths up to 15 m with throw below a 2 m threshold. Plotting opening width versus vertical offset shows that there is a continuous transition between the endmembers. We conclude that for these endmembers, the ratio between opening width and vertical offset R can be reliably used to predict fault structures at depth. However, fractures associated with MDFs belong to one larger continuum and, consequently, where different endmembers coexist, a clear identification of structures solely via the determination of R is impossible.

No abstract

Full-resolution UAV-based image of the joint set of the classic Lilstock benches (UK) • Layer-bound joints are fully imaged over an entire large outcrop • Up to eight sets of joints occur in a single limestone layer • Jointing is laterally heterogeneous in the same layer and different between layers • Phases of sealing accompanied the evolution of older joints at Lilstock