“…We first construct a highly segmented model of the HFFZ (Figure 1b, Model-A), consisting of 55 partially intersecting, non-planar vertical faults, each also intersecting with the complex bathy-topography of the free surface. We integrate data from highresolution bathymetry interpretation, o shore seismic reflection campaigns in Northern Iceland (Magnúsdóttir & Brandsdóttir, 2011;Magnúsdóttir et al, 2015;Hjartardóttir et al, 2016) and relocated seismicity (Abril et al, 2018(Abril et al, , 2019 (Figure 1a). We assume vertical faults which is supported by the depth distribution of the recently relocated HFFZlocal seismicity from 1993 to 2019 considering the lateral variability of the crustal structure (Abril et al, 2018(Abril et al, , 2019(Abril et al, , 2021.…”