Description: The Sea of Marmara is traversed by the North Anatolian Fault system and also presents abundant emission sites of methane gas into the water column. In order to assess the spatial relationship between gas emissions and active faults, the distribution of distances between gas emission sites and the nearest fault is calculated and compared with the distribution of distances between a uniform random distribution of points (Poisson process representing the null hypothesis of an absence of relationship between gas emissions and faults) and the nearest fault. Interestingly, the distance distribution for the Poisson process is nearly exponential, indicating that the fault map does not have a characteristic scale other than that representing the intensity of the fault network. The distance distribution for the observed gas emissions is significantly narrower than that arising from the Poisson process, with a Kolmogorov distance of 0.25 ± 0.02. The crossing point between the two distributions defines the characteristic half-width of the swath of gas emission sites around the mapped active faults. For the whole Sea of Marmara data set a characteristic half-width of 900–1000 m is found which matches the half-width of the seafloor deformation zone observed around the main active fault. When the same analysis is applied to zones covering the Western High and the Central High it is found that the swath of gas emissions is wider on the Central High (2 km half-width), and not clearly related to the seafloor deformation zone there. This difference is put in perspective with recent work showing that creep is occurring on the western segment of the Main Marmara Fault (this also causing microseismicity) while the central Istanbul-Silivri segment may have remained locked since the 1766 magnitude 7+ earthquake. This suggests that aseismic slip (and not only earthquake occurrence) effectively maintains high permeability conduits in fault zones in sediments.

Description: Understanding micro-seismicity is a critical question for earthquake hazard assessment. Since the devastating earthquakes of Izmit and Duzce in 1999, the seismicity along the submerged section of North Anatolian Fault within the Sea of Marmara (comprising the “Istanbul seismic gap”) has been extensively studied in order to infer its mechanical behaviour (creeping vs locked). So far, the seismicity has been interpreted only in terms of being tectonic-driven, although the Main Marmara Fault (MMF) is known to strike across multiple hydrocarbon gas sources. Here, we show that a large number of the aftershocks that followed the M 5.1 earthquake of July, 25th 2011 in the western Sea of Marmara, occurred within a zone of gas overpressuring in the 1.5–5 km depth range, from where pressurized gas is expected to migrate along the MMF, up to the surface sediment layers. Hence, gas-related processes should also be considered for a complete interpretation of the micro-seismicity (~M 〈 3) within the Istanbul offshore domain.

Description: High-resolution 3-D seismic data acquired in the Sea of Marmara on the Western High, along the northwestern branch of the North Anatolian Fault (also known as the Main Marmara Fault), shed new light on the evolution of the deformation over the last 500–600 ka. Sedimentary sequences in ponded basins are correlated with glacioeustatic cycles and transitions between marine and low sea/lake environments in the Sea of Marmara. In the 3 × 11 km2 of the 3-D seismic survey, deformation over the last 405–490 ka is localized along the main fault branch and north of it, where N130°–N140° trending normal faults and N40°–N50° folding accommodated strike-slip deformation associated with active argillokinesis. There is some evidence that deformation was more distributed further back in the past, at least over the depth range (〈600 m below seafloor) of our survey. A N110° basin and buried ridge system were eventually cut by the presently active fault. The southern part of the basin was then uplifted, while the northern part was folded but continued to subside along the fault. A mass transport deposits complex dated between 405–490 ka shows a lateral displacement of 7.7 ± 0.3 km, corresponding to an estimated slip rate of 15.1–19.7 mm/a. We conclude that this strand of the Main Marmara Fault on the Western High has taken up most of the strike slip motion between the Anatolian and Eurasian plates over the last 405 ka at least.

Description: The free‐air gravity in the Marmara Sea reveals that the low density of sedimentary basins is partly compensated in the lower crust. We compiled geophysical upper crust studies to determine the sediment basin geometries in and around the Marmara Sea and corrected the gravity signal from this upper crust geology with the Parker method. Then, assuming long wavelength anomalies in the residual gravity signal is caused by variations in the Moho topography, we inverted the residual to build the Moho topography. The result shows that the Moho is uplifted on an area greater than the Marmara Sea with a maximum crust thinning beneath the basins where the Moho is at about 25 km, 5 km above the reference depth. We then evaluated the Neogene extension by comparing the surface covered by our 3‐D thinned model with the surface covered by an unthinned model with same crustal volume. Comparing this surface with areal extension rate from GPS data, we found a good compatibility indicating that the extension rate averaged over the Sea of Marmara area probably remained close to its present‐day value during major changes of tectonic regime, as the incursion of the North Anatolian Fault system during the Pliocene leads to the establishment of the dominantly strike‐slip present‐day system. We also show that crustal extension is distributed over a wider domain in the lower crust than in the upper crust, and that this may be accounted for by a relatively minor component of lower crustal ductile flow.