Slip-rate variability and distributed deformation in the Marmara Sea fault system

Abstract

The slip rate along a fault controls the accumulation of strain that is eventually released during an earthquake. Estimates from a three-dimensional geomechanical model of the slip rate on the main Marmara fault near Istanbul, Turkey reconcile geodetic and geological observations and indicate smaller values and higher variability than previously thought. The slip rate along a fault controls the accumulation of strain that is eventually released during an earthquake1. Along a 150-km-long stretch of the North Anatolian fault near Istanbul, Turkey, strain has been building up2 since the last large earthquake in 1766. Estimates of the geodetic slip rates along the main Marmara fault vary widely, ranging between 17 and 27.9 mm yr−1 (refs 2, 3, 4, 5). This slip rate is difficult to quantify because of the lack of satellite observations offshore and the complexity of the submarine fault system that includes the main Marmara fault2,6,7. Here we estimate the right-lateral slip rate on the main Marmara fault using a three-dimensional geomechanical model that incorporates these structural complexities. From our simulations we infer slip rates between 12.8 and 17.8 mm yr−1; our estimates are smaller and more variable than previous results, primarily because of slip partitioning and internal deformation. Our model results reconcile geodetic observations and geological fault slip rates8,9,10, which had been considered conflicting previously. We suggest that the inferred variability in slip rate on the main Marmara fault favours segmented release of seismic moment during consecutive events over the failure of the whole seismic gap in one large earthquake.