Modelling of faults in LoopStructural 1.0

Without properly accounting for both fault kinematics and observations of a faulted surface, it is challenging to create 3D geological models of faulted geological units. Geometries where multiple faults interact, where the faulted surface geometry significantly deviate from a flat plane and where the geological interfaces are poorly characterised by sparse datasets are particular challenges. There are two existing approaches for incorporating faults into geological surface modelling. One approach incorporates the fault displacement into the surface description but does not incorporate fault kinematics and in most cases will produce geologically unexpected results such as shrinking intrusions, fold hinges without offset and layer thickness growth in flat oblique faults. The second approach builds a continuous surface without faulting and then applies a kinematic fault operator to the continuous surface to create the displacement. Both approaches have their strengths; however, neither approach can capture the interaction of faults within complicated fault networks, e.g. fault duplexes, flower structures and listric faults because they either (1) impose an incorrect (not defined by data) fault slip direction or (2) require an over-sampled dataset that describes the faulted surface location. In this study, we integrate the fault kinematics into the implicit surface, by using the fault kinematics to restore observations, and the model domain prior to interpolating the faulted surface. This new approach can build models that are consistent with observations of the faulted surface and fault kinematics. Integrating fault kinematics directly into the implicit surface description allows for complexly faulted stratigraphy and fault–fault interactions to be modelled. Our approach shows significant improvement in capturing faulted surface geometries, especially where the intersection angle between the faulted surface and the fault surface varies (e.g. intrusions, fold series) and when modelling interacting faults (fault duplex).

Slow-slip, earthquake-swarms and fault-interactions at the western-end of the Hellenic Subduction System precede the Mw 6.9 Zakynthos Earthquake, Greece

<p>The month-to-year-long deformation of the Earth’s crust where active subduction zones terminate is poorly explored. Here we report on a multidisciplinary dataset that captures the synergy of slow-slip events, earthquake swarms and fault-interactions during the ~5 years leading up to the 2018 M<sub>w</sub> 6.9 Zakynthos Earthquake at the western termination of the Hellenic Subduction System (HSS). It appears that this long-lasting preparatory phase initiated due to a slow-slip event that lasted ~4 months and released strain equivalent to a ~M<sub>w</sub> 6.3 earthquake. We propose that the slow-slip event, which is the first to be reported in the HSS, tectonically destabilised the upper 20-40 km of the crust, producing alternating phases of seismic and aseismic deformation, including intense microseismicity (M<4) on neighbouring faults, earthquake swarms in the epicentral area of the M<sub>w</sub> 6.9 earthquake ~1.5 years before the main event, another episode of slow-slip immediately preceding the mainshock and, eventually, the large (M<sub>w </sub>6.9) Zakynthos Earthquake. Tectonic instability in the area is evidenced by a prolonged (~4 years) period of overall suppressed b-values (<1) and strong earthquake interactions on discrete strike-slip, thrust and normal faults. We propose that composite faulting patterns accompanied by alternating (seismic/aseismic) deformation styles may characterise multi-fault subduction-termination zones and may operate over a range of timescales (from individual earthquakes to millions of years).</p>

Fault interactions in a complex fault system: insight from the 1936–1997 NE Lut earthquake sequence

SUMMARY Calculations of Coulomb stress changes have shown that moderate to large earthquakes may increase stress at the location of future earthquakes. Coulomb stress transfers have thus been widely accepted to explain earthquake sequences, especially for sequences occurring within parallel or collinear fault systems. Relating, under this framework, successive earthquakes occurring within more complex fault systems (i.e. conjugate fault system) is more challenging. In this study, we assess which ingredients of the Coulomb stress change theory are decisive for explaining the succession of three large (Mw 7+) earthquakes that occurred on a conjugate fault system in the NE Lut, East Iran, during a 30-yr period. These earthquakes belong to a larger seismic sequence made up of 11 earthquakes (Mw 5.9+) from 1936 to 1997. To reach our goal, we calculate, at each earthquake date, the stress changes generated by the static deformation of the preceding earthquakes, the following post-seismic deformation due to the viscoelastic relaxation of the lithosphere, and the interseismic deformation since 1936. We first show that accurately modelling the source and receiver fault geometry is crucial to precisely estimating Coulomb stress changes. Then we show that 7 out of 10 earthquakes of the NE Lut sequence, considering the uncertainties, are favoured by the previous earthquakes. Furthermore, the last two M7+ earthquakes of the sequence (1979 and 1997) have mainly been favoured by the moderate Mw ∼ 6 earthquakes. Finally, we investigate the link between the Coulomb stress changes due to previous earthquakes and the rupture extension of the next earthquake and show that a correlation does exist for some earthquakes but is not systematic.