Special issue “The next Marmara earthquake: disaster mitigation, recovery, and early warning”

© The Author(s) 2017. This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Marmara Sea, accommodating the fault segments of a major transform fault, is well known as a seismic gap along the North Anatolian Fault (NAF), running through the northern part of Turkey and connecting the East Anatolian convergent area with the Hellenic subduction zone (e.g., Pınar 1943; Toksöz et al. 1979; Pondard et al. 2007; Şengör et al. 2014). It is obvious from historical records spanning more than 2000 years that the region is subject to frequent strong shaking that is likely associated with tsunami waves, threatening heavily populated and industrialized locations (Ambraseys 2002; Erdik et al. 2004; Hébert et al. 2005). In the twentieth century, magnitude (M) 7-class earthquakes sequentially occurred from east to west along the NAF zone, as shown in Fig. 1 (Stein et al. 1997). The last two successive events hit the eastern Marmara region in August 17 and November 12, 1999, known as Izmit (Mw = 7.5) and Düzce (Mw = 7.2) earthquakes, respectively, killing about 20,000 people and devastating the region. On the other hand, at the western edge of Marmara Sea, an M ~ 7 earthquake occurred in 1912, rupturing onshore and offshore fault segments where 4–5 m lateral displacements were measured (Armijo et al. 2005; Aksoy et al. 2010). Furthermore, in the Marmara Sea, it has been considered that the last M7 class earthquake occurred in 1766. With an average slip rate of 2 cm/year, several meters have accumulated over the past 250 years (Straub 1996; Meade et al. 2002). Based on a time-dependent model that includes the coseismic and post-seismic effects of the 1999 Kocaeli earthquake with moment magnitude Mw = 7.4, Parsons (2004) concluded that the probability of an earthquake with Mw > 7 in the Sea of Marmara near Istanbul is 35–70% in the next 30 years. According to a 2011 study, an earthquake with Mw = 7.25 on the Main Marmara Fault is expected to heavily damage or destroy 2–4% of the near 1,000,000 buildings in Istanbul, which has a population of around 13 million, with 9–15% of the buildings receiving medium damage and 20–34% of the buildings lightly damaged (Erdik 2013). The aim of this special issue is to gather information about the risk of another Marmara earthquake from the latest geophysical, geological, geotechnical, computational, and building science research results to discuss ways of mitigating disaster in advance. The collection of 12 papers constituting this special issue is based on recent research on imaging the crustal structure, the geometry of the fault segments and their microseismicity features, source characteristics of large earthquakes inferred from historical seismograms, tsunami hazard assessment and mitigation studies, site response evaluations, and development of an integrated earthquake simulation system. Polat et al. (2016) portrayed a tomographic image of the crust beneath the Sea of Marmara and its close surrounding area, derived from a dataset of more than ninety thousand Pand S-wave arrival times from local earthquakes recorded at land-based stations and ocean bottom seismographs. Their seismic velocity images illustrate positive anomalies nearby the fault segments of NAF, correlating with regions of higher seismicity. The study of Yamamoto et al. (2015) closely inspects the NAF segments beneath the Sea of Marmara, especially the fault geometry and their seismic activity carried out by the deployment of 15 OBS stations at locations close to the Main Marmara Fault (MMF). A deeper than previously known seismogenic zone, extending to the lower crust, is one of the most striking features derived from the OBS observations. The offshore seismic stations nearby MMF constrained the depth, dip, and lateral extent of the seismicity much better and, consequently, the geometry of the fault Open Access

The Marmara Sea, accommodating the fault segments of a major transform fault, is well known as a seismic gap along the North Anatolian Fault (NAF), running through the northern part of Turkey and connecting the East Anatolian convergent area with the Hellenic subduction zone (e.g., Pınar 1943;Toksöz et al. 1979;Pondard et al. 2007;Şengör et al. 2014). It is obvious from historical records spanning more than 2000 years that the region is subject to frequent strong shaking that is likely associated with tsunami waves, threatening heavily populated and industrialized locations (Ambraseys 2002;Erdik et al. 2004;Hébert et al. 2005). In the twentieth century, magnitude (M) 7-class earthquakes sequentially occurred from east to west along the NAF zone, as shown in Fig. 1 (Stein et al. 1997). The last two successive events hit the eastern Marmara region in August 17 and November 12, 1999, known as Izmit (Mw = 7.5) and Düzce (Mw = 7.2) earthquakes, respectively, killing about 20,000 people and devastating the region. On the other hand, at the western edge of Marmara Sea, an M ~ 7 earthquake occurred in 1912, rupturing onshore and offshore fault segments where 4-5 m lateral displacements were measured (Armijo et al. 2005;Aksoy et al. 2010). Furthermore, in the Marmara Sea, it has been considered that the last M7 class earthquake occurred in 1766. With an average slip rate of 2 cm/year, several meters have accumulated over the past 250 years (Straub 1996;Meade et al. 2002). Based on a time-dependent model that includes the coseismic and post-seismic effects of the 1999 Kocaeli earthquake with moment magnitude Mw = 7.4, Parsons (2004) concluded that the probability of an earthquake with Mw > 7 in the Sea of Marmara near Istanbul is 35-70% in the next 30 years. According to a 2011 study, an earthquake with Mw = 7.25 on the Main Marmara Fault is expected to heavily damage or destroy 2-4% of the near 1,000,000 buildings in Istanbul, which has a population of around 13 million, with 9-15% of the buildings receiving medium damage and 20-34% of the buildings lightly damaged (Erdik 2013).
The aim of this special issue is to gather information about the risk of another Marmara earthquake from the latest geophysical, geological, geotechnical, computational, and building science research results to discuss ways of mitigating disaster in advance. The collection of 12 papers constituting this special issue is based on recent research on imaging the crustal structure, the geometry of the fault segments and their microseismicity features, source characteristics of large earthquakes inferred from historical seismograms, tsunami hazard assessment and mitigation studies, site response evaluations, and development of an integrated earthquake simulation system. Polat et al. (2016) portrayed a tomographic image of the crust beneath the Sea of Marmara and its close surrounding area, derived from a dataset of more than ninety thousand P-and S-wave arrival times from local earthquakes recorded at land-based stations and ocean bottom seismographs. Their seismic velocity images illustrate positive anomalies nearby the fault segments of NAF, correlating with regions of higher seismicity. The study of Yamamoto et al. (2015) closely inspects the NAF segments beneath the Sea of Marmara, especially the fault geometry and their seismic activity carried out by the deployment of 15 OBS stations at locations close to the Main Marmara Fault (MMF). A deeper than previously known seismogenic zone, extending to the lower crust, is one of the most striking features derived from the OBS observations. The offshore seismic stations nearby MMF constrained the depth, dip, and lateral extent of the seismicity much better and, consequently, the geometry of the fault Open Access *Correspondence: horit@jamstec.go.jp segments. The geometry between the fault segments and the orientation of the maximum horizontal compressive stress axis is correlated through the frictional strength of faults in the study presented by Pinar et al. (2016). They use the frequency distribution of P-and T-axes of the focal mechanisms of events taking place around three fault segments in the eastern Marmara. A low frictional coefficient for the Princes' Islands segment and high frictional coefficients for the Yalova-Çınarcık and Yalova-Hersek segments are derived by the authors.
Further information on the fault segmentation of NAF in the Marmara Sea and north Aegean Sea regions is determined by modeling the waveforms at broadband land stations to retrieve the source properties of moderate to large events, estimating the CMT parameters and source time functions (Nakano et al. 2016). The source time function of the 2014 North Aegean earthquake (Mw = 6.9) provides evidence on ruptures of different fault segments associated with the mainshock. Moreover, the CMT solutions derived for the events taking place along the MMF fault indicate fault parallel T-axis orientations, suggesting segmentation and development of local extension features. Baştürk et al. (2016) relocated two ancient moderate-sized earthquakes, using historical seismograms of events occurring in 1935 in the proximity of Marmara Island. Applying modern techniques to the old seismograms, CMT parameters for the two events are retrieved that indicate predominantly normal faulting mechanisms taking place at the upper crust. Similar source parameters are obtained for the 1963 Çınarcık earthquake, whose location was long a subject of debate. The location of the fault segment ruptured by the 1963 earthquake and its source mechanism is of great importance in seismic hazard assessment studies. Mert et al. (2016) tried to predict the strong ground motions to be generated by a rupture on the Princes' Islands segment, which lies beneath the Çınarcık basin, using physically based probabilistic seismic hazard analysis (PSHA) methodology where broadband strong ground motion simulations are conducted. To generate the high-frequency (0.5-20 Hz) part of the broadband earthquake simulation, real, small-magnitude earthquakes recorded by a local seismic array were used as empirical Green's functions. As for frequencies below 0.5 Hz, the synthetic Green's functions are calculated by an explicit 2D/3D elastic finite difference wave propagation routine. Aytore et al. (2016) used NAMI-DANCE code to carry out high-resolution tsunami simulations in the Marmara Sea, focusing on Haydarpaşa Port in the megacity of Istanbul, and computed the tsunami parameters in and around the port. They observed that the stability of the breakwaters is one of the major factors that influence whether agitation and inundation can be diminished in the event of a tsunami in Haydarpaşa Port, as harbor protection structures have not been designed to withstand tsunamis. The flow depth, momentum fluxes, and current pattern were identified as the other factors that cause unexpected circulation and uncontrolled movements of objects on land and vessels in the sea. Cankaya et al. (2016) applied a new methodology in Yenikapı region as a case study for tsunami vulnerability assessment, based on high-resolution coastal inundation modeling. Using NAMI-DANCE code, they constructed vulnerability at location and evacuation resilience maps, using the analytical hierarchical process (AHP) method of multicriteria decision analysis (MCDA) to define the tsunami human vulnerability of the region. The vulnerability at location map is composed of metropolitan use, geology, elevation, and distance from shoreline layers, whereas the evacuation resilience map considers slope, distance within flat areas, distance to buildings, and distance to road networks layers. Following this, they computed the tsunami risk map from the proposed new relationship, which uses flow depth maps, vulnerability at location Fig. 1 Surface rupture extent associated with major earthquakes along NAF maps, and evacuation resilience maps. Necmioglu (2016) indicated that in the absence of adequate post-earthquake assembly areas, especially in heavily urbanized Istanbul, citizens would be rushing to landfill assembly and recreational areas in the coastal parts of the city after a major earthquake. To address this, he proposed a model for a tsunami warning system specific for the Marmara region that is strongly coupled with the earthquake early warning system and stakeholders of tsunami mitigation activities, such as the local and regional components of disaster and emergency management and civil protection units. This would ensure that the citizens would stay away from the coastline in case of a large earthquake. Necmioglu also discussed associated challenges, such as decoupled earthquake and tsunami mitigation activities in the Marmara region.
Istanbul is the largest metropolitan city expected to be hit by the impending Marmara Sea earthquake. Rapid response systems are effective in mitigating the loss of life and property. Zülfikar et al. (2017) describe how realtime ground motion shaking maps are constructed from the strong motion stations distributed throughout the densely populated areas of the city. Building damage estimation is then computed by using grid-based building inventory, and the related loss is estimated. Zülfikar et al. (2017) further inform how rapidly estimated data enable public and private emergency management authorities to take action to allocate and prioritize resources, minimizing casualties in urban areas during immediate postearthquake periods. Site response plays an important role in generating shake maps. Karagoz et al. (2015) explore the S-wave velocity structure of shallow soils using microtremors to estimate site responses in Tekirdag and the surrounding area. They collected data at 44 sites in Tekirdag, Marmara Ereglisi, Corlu, and Muratlı, and estimated phase velocities of Rayleigh waves from the microtremor data using a spatial autocorrelation method. A hybrid genetic simulated annealing algorithm was applied by the authors to obtain a 1D S-wave velocity structure at each site. All the studies constituting this special issue contribute to some extent to an integrated earthquake simulation (IES) system that was developed for Istanbul by Sahin et al. (2016). The IES for Istanbul is built in MATLAB and includes site response analysis as well as structural seismic response analysis of existing buildings; building models are made by using GIS databases. An initial application is performed in the Zeytinburnu District of Istanbul, and the results are expressed in the form of spatial distribution of ground motion and building responses. The IES analysis illustrates a nonuniform distribution of seismic responses, indicating the possibility that there are more severely damaged areas in the district compared with others.
The Marmara region, long ago identified as a seismic gap, is waiting for the impending major earthquake to rupture the MMF segments. Meanwhile, several studies have been carried out to determine the seismic hazard level, the fault geometry and segmentation, mitigation of seismic risk, rapid response, and early warning. All these subjects are tackled in this special issue, making a valuable contribution on the existing information that is constrained by high-quality data acquired in the frame of the MarDIM project, a joint effort of Turkish and Japanese scientists financially supported by JICA, JST and the Ministry of Development in Turkey.
Author details 1 Japan Agency for Marine-Earth Science and Technology, Yokohama, Japan. 2 Kandilli Observatory and Earthquake Research Institute, Bogazici University, Istanbul, Turkey. 3 Earthquake Research Institute, The University of Tokyo, Tokyo, Japan. 4 Hydrographic and Oceanographic Department, Japan Coast Guard, Tokyo, Japan.

Publisher's Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.