Seismic behaviour of the North Anatolian Fault beneath the Sea of Marmara (NW Turkey): implications for earthquake recurrence times and future seismic hazard

Possible long-term seismic behaviour of the Northern strand of the North Anatolian Fault Zone, between western extreme of the 1999 İzmit rupture and the Aegean Sea, after 400 AD is studied by examining the historical seismicity, the submarine fault mapping and the paleoseismological studies of the recent scientific efforts. The long-term seismic behaviour is discussed through two possible seismicity models devised from M _S ≥ 7.0 historical earthquakes. The estimated return period of years of the fault segments for M1 and M2 seismic models along with their standard deviations are as follows: F4 segment 255 ± 60 and 258 ± 12; F5 segment 258 ± 60 and 258 ± 53; F6 segment 258 ± 60 and 258 ± 53; F7 segment 286 ± 103 and 286 ± 90; F8 segment 286 ± 90 and 286 ± 36. As the latest ruptures on the submarine segments have been reported to be during the 1754–1766 earthquake sequence, and the 1912 mainshock rupture has been evidenced to extend almost all over the western part of the Sea of Marmara, our results imply imminent seismic hazard and, considering the mean recurrence time, a large earthquake to strike the eastern part of the Sea of Marmara in the next two decades.     

A quantitative appraisal of earthquake hazard parameters computed from Gumbel I method for different regions in and around Turkey

Useful information concerning the earthquake hazard parameters distributed in Turkey and the adjacent areas are estimated in the present work. Based on Gumbel’s I distribution parameters we are able to estimate the hazard values of the investigated area which are the mean return periods, the most probable maximum magnitude in the time period of t-years and the probability for an earthquake occurrence of magnitude ≥M during a time span of t-years. Figures concerning the spatial distribution of probabilities and the return periods are plotted and we considered them of particular interest for mapping the earthquake hazard in Turkey and the surrounding areas. These figures effectively produce a brief earthquake hazard atlas. The quantitative appraisal of the hazard parameters is useful for engineers, planners, etc., because it provides a tool for earthquake resistant design.     

Characteristics of the 1912 co-seismic rupture along the North Anatolian Fault Zone (Turkey): implications for the expected Marmara earthquake

Evidence of right‐lateral offsets associated with the 1912 earthquake (M_w 7.4) along the North Anatolian Fault (Gaziköy–Saros segment) allow us to survey (using DGPS) the co‐seismic and cumulative slip distribution. The damage distribution and surface breaks related with the earthquake show an elongated zone of maximum intensity (X MSK) parallel to the fault rupture on land but this may extend offshore to the north‐east and south‐west. Detailed mapping of the fault using topographic maps and aerial photographs indicates the existence of pull‐apart basins and pressure ridges. At several localities, the average 1912 offset along strike is 3.5–4 m and cumulative slip is 2–6 times that of individual movement. The fault rupture geometry and slip distribution suggest the existence of three subsegments with a combined total length of 110–120 km, a fault length and maximum slip similar to those of the 1999 Izmit earthquake. The amount of slip at the north‐easternmost section and in the coastal region of the Sea of Marmara reaches an average 4 m, thereby implying the offshore extension of the 1912 rupture. The results suggest that the 1912 event generated up to 150 km of surface faulting, which would imply a M_w 7.2–7.4 earthquake and which, added with rupture lengths of the 1999 earthquakes, help to constrain the remaining seismic gap in the Sea of Marmara.     

Structural relationship between the Sea of Marmara Basin and the North Anatolian Fault Zone

The North Anatolian Fault (NAF) zone is 1500 km long, extending almost up to the Greek mainland in the west. It is a seismically active right-lateral strike-slip fault that accommodates the relative motion between the Turkish block and Black Sea plate. The Sea of Marmara lies along the western part of the NAF and shows evidence of subsidence. In this area pure strike-slip motion of the fault zone changes into extensional strike-slip movement that is responsible for the creation of the Sea of Marmara and the North Aegean basins. The northern half of the Sea of Marmara is interpreted as a large pull-apart basin. This basin is subdivided into three smaller basins separated by strike-slip fault segments of uplifted blocks NE-SW. Basinal areas are covered by horizontally layered sedimentary sequences. Uplifted blocks have undergone compressional stress. All the blocks are subsiding and are undergoing vertical motions and rotations relative to one another. The uplifted blocks exhibit positive Bouguer gravity anomalies. According to gravity interpretation, there is relative crustal thinning under the Sea of Marmara. The northern side of the Sea of Marmara is marked by a distinctive deep-rooted magnetic anomaly, which is dissected and shifted southward by strike-slip faulting. The southern shelf areas of the Sea of Marmara are dominated by short-wavelength magnetic anomalies of shallow origin.     

Evaluation of earthquake potential along the Northern Anatolian Fault Zone in the Marmara Sea using comparisons of GPS strain and seismotectonic parameters

Seismotectonic parameters including the Gutenberg-Richter b-value and multifractal dimensions D₂ and D₁₅ of seismicity patterns (both spatial and temporal) were compared to GPS-derived maximum shear and dilatation strains measured in the Marmara Sea region of western Turkey along the Northern Anatolian Fault Zone (NAFZ). Comparisons of seismotectonic parameters and GPS-derived maximum shear and dilatation strain along the NAFZ in the vicinity of the 1999 M7.4 Izmit earthquake reveal a positive correlation (r = 0.5, p = 0.05) between average dilatation and the Gutenberg-Richter b-value. Significant negative correlation (r = − 0.56, p = 0.03 and r = − 0.56, p = 0.02) was also observed between the spatial fractal dimension D₂ and GPS-derived maximum geodetic and shear strain. This relationship suggests that, as maximum geodetic and shear strains increase, seismicity becomes increasingly clustered.Anomalous interrelationships are observed in the Marmara Sea region prior to the Izmit event along a bend in the NAFZ near the eastern end of the Marmara Sea known as the Northern Boundary Fault (NBF). An asperity is located near the northwest end of the NBF. Along the 50-km length of the NBF, GPS strains become slightly compressive. The correlation between b-value and GPS-derived dilatation suggests that regions in compression have increased probability of larger magnitude rupture. The NBF appears to serve as an impediment to the transfer of strain from east to west along the NAFZ. Recurrence times for large earthquakes along the NBF are larger than in surrounding areas. Temporal clustering of seismicity in the vicinity of the NBF may represent foreshocks of an impending rupture.     

Improving the level of seismic hazard parameters in Saudi Arabia using earthquake location

Saudi Arabia is characterized as largely aseismic; however, the tectonic plate boundaries that surround it are very active. To improve characterization of seismicity and ground motion hazard, the Saudi Arabian Digital Seismic Network (SANDSN) was installed in 1998 and continues to be operated by the Saudi Geological Survey (SGS) and King Abdulaziz City for Science and Technology (KACST). This article describes research performed to improve seismic hazard parameters using earthquake location and magnitude calibration of the high-quality SANDSN data. The SANDSN consists of 38 seismic stations, 27 broadband, and 11 short period. All data are telemetered in real time to a central facility at KACST in Riyadh. The SANDSN stations show low background noise levels and have good signal detection capabilities; however, some stations show cultural noise at frequencies above 1.0 Hz. We assessed the SANDSN event location capabilities by comparing KACST locations with well-determined locations derived from ground truth or global observations. While a clear location bias exists when using the global average iasp91 earth model, the locations can be improved by using regional models optimized for different tectonic source regions. The article presents detailed analysis of some events and Dead Sea explosions where we found gross errors in estimated locations. New velocity models we calculated that should improve estimated locations of regional events in three specific regions include (1) Gulf of Aqabah—Dead Sea region, (2) Arabian Shield, and (3) Arabian Platform. Recently, these models were applied to the SANDSN to improve local and teleseismic event locations and to develop an accurate magnitude scale for Saudi Arabia. The Zagros Thrust presents the most seismic hazard to eastern Saudi Arabia because of the frequent occurrence of earthquakes. Although these events are 200 km or further from the Arabian coast, wave propagation through sedimentary structure of the Gulf causes long-duration ground motions for periods between 3 and 10 s. Such ground motions could excite response in large engineered structures (e.g., tall buildings and long bridges) such as was experienced after the November 22, 2005 Qeshm Island earthquake off the southern coast of Iran.     

Activity level of tectonic basins,western section of the North Anatolian Fault Zone,Turkey

We employed quantitative techniques to investigate tectonic activity levels and development stages of the Bolu, Yenicaga, Dortdivan, Cerkes, Ilgaz, and Tosya structural basins along the western portions of the main trace of the North Anatolian Fault Zone (NAFZ). Our methodology incorporates six morphometric indices: basin shape (basin elongation and compactness), hypsometric integral, mountain-front sinuosity, stream length gradient index, valley floor width-to-height ratio, and asymmetry factor, obtained from the digital elevation model of the region generated from 1/25,000-scale topographic maps. These indices are integrated within the framework of an analytical hierarchy process to provide relative activity level values of the individual basins. The new analyses indicate that the basins have contrasting tectonic activity characteristics. Judging from the applied indices, the relative increasing order of the tectonic basin activity is Dortdivan, Cerkes, Yenicaga, Ilgaz, Tosya, and Bolu. Among the basins located to the north of the NAFZ, the activity decreases eastwards, whereas to the south of this profound fault zone, it decreases towards the west.     

Depth variation of seismic source scaling relations: implications for earthquake hazard in southeastern Australia

Advances in earthquake data acquisition and processing techniques have allowed for improved quantification of source parameters for local Australian earthquakes. Until recently, only hypocentral locations and local magnitudes (M_L) had been determined routinely, with little attention given to the inversion of additional source parameters. The present study uses these new source data (e.g. seismic moment, stress drop, source dimensions) to further extend our understanding of seismicity and the continental stress regime of the Australian landmass and its peripheral regions.

Earthquake activity within Australia is typically low, and the proportion of small to large events (i.e. the b value) is also low. It is observed that average stress drops for southeastern Australian earthquakes appear to increase with seismic moment to relatively high levels, up to approximately 10 MPa for M_L 5.0 earthquakes. This is thought to be indicative of high compressive crustal stress, coupled with strong rocks and fault asperities. Furthermore, the data indicates that shallow focus earthquakes (shallower than 6 km) appear to produce lower than average stress drops than deeper earthquakes (between 6 and 20 km) with similar moment.

Recurrence estimates were obtained for a discrete seismogenic zone in southeastern Australia. Decreasing b values with increasing focal depth for this zone indicate that larger earthquakes (with high stress drops) tend to occur deeper in the crust. This may offer an explanation for the apparent increase of stress drop with hypocentral depth. Consequently, earthquake hazard estimates that assume a uniform Gutenburg–Richter distribution with depth (i.e. constant b value) may be too conservative and therefore slightly overestimate seismic hazard for surface sites in southeastern Australia.     

Effect of Aftershocks on Earthquake Hazard Estimation: An Example from the North Anatolian Fault Zone

In order to investigate the effect of aftershocks on earthquake hazard estimation, earthquake hazard parameters (_m, b and M_max) have been estimated by the maximum likelihood method from the main shocks catalogue and the raw earthquakes catalogue for the North Anatolian Fault Zone (NAFZ). The main shocks catalogue has been compiled from the raw earthquake catalogue by eliminating the aftershocks using the window method. The raw earthquake catalogue consisted of instrumentally detected earthquakes between 1900 and 1992, and historical earthquakes that occurred between 1000–1900. For the events of the mainshock catalogue the Poisson process is valid and for the raw earthquake catalogue it does not fit. The paper demonstrates differences in the hazard outputs if on one hand the main catalogues and on the other hand the raw catalogue is used. The maximum likelihood method which allows the use of the mixed earthquake catalogue containing incomplete (historical) and complete (instrumental) earthquake data is used to determine the earthquake hazard parameters. The maximum regional magnitude (M_max, the seismic activity rate (_m), the mean return period (R) and the b value of the magnitude-frequency relation have been estimated for the 24°–31° E, 31°–41° E, 41°–45° E sections of the North Anatolian Fault Zone from the raw earthquake catalogue and the main shocks catalogue. Our results indicate that inclusion of aftershocks changes the b value and the seismic activity rate _m depending on the proportion of aftershocks in a region while it does not significantly effect the value of the maximum regional magnitude since it is related to the maximum observed magnitude. These changes in the earthquake hazard parameters caused the return periods to be over- and underestimated for smaller and larger events, respectively.     

Seismo-Tectonic Characteristics of the North Anatolian Fault Zone Between Akyazi and Duzce (Bolu,Turkey)

The active North Anatolian fault zone (NAFZ) presents very complex seismotectonic activity. The occurrence of the Abant earthquake in 1957 (Ms = 7.1) and the Mudurnu earthquake in 1967 (Ms = 6.8) are only two examples of several seismic events associated with intense tectonic activity of the NAFZ. Statistical analyses of earthquakes in an area extending between 30° 30′ to 31° 30' E Long. and 40° 15′ to 41° 00′ N Lat. reveal that epicenters generally were shallow. However, a few deep epicenters also were located, some of which reached a depth of 30 km. The epicenters were found to concentrate in a zone lying between the Duzce and Akyazi Plain to the north of Almacik Mountain and in the Adapazari Plain. The Northern Anatolian fault displays an en echelon character in the area, except for the eastern part, where it extends as a single segment. The en echelon character of the NAFZ is interpreted as a structure distributing the potential energy and consequently reducing the intensity of earthquakes, giving rise to micro-earthquakes of magnitudes less than 4.2.     

郯庐断裂带中南段闭锁特征与地震危险性分析

利用华北地区1999～2007、2013～2017两期GPS水平运动速度场数据,采用块体负位错模型,分别反演了郯庐断裂带中南段不同段的断层闭锁程度和滑动亏损速率分布;结合地表应变结果,综合分析了郯庐断裂带前后两期的变形差异特征,并探讨了其与日本3·11地震间的可能关系。研究结果表明:日本地震后,郯庐断裂带中南段郯城以北的段落闭锁程度有所减弱,中南段东部地区主张应变率增强,处于拉张状态;日本大地震的发生对郯庐断裂带中南段的应变积累起到一定的缓解作用。2013～2017最新一期反演结果显示莒县以北断层闭锁程度仍较高,闭锁深度较深,为右旋挤压亏损,是1668年郯城地震的未破裂段;莒县以南到泗洪附近断层闭锁程度较低,无滑动亏损积累;泗洪以南到嘉山段断层闭锁程度较高,是历史地震的未破裂段,同时该地区小震不活跃,易于应力积累,地震危险性值得关注。     

1999年9月21日台湾集集地震的震源参数综合判定

本文利用半经验、半理论的方法对世界各主要地震汇集和研究机构所给出的集集主震基本震源参数进行了综合计算和评定。根据这些参数重新计算了集集主震的标量地震矩、破裂面积、应力降等。根据这些参数之间的内在联系 ,从地震分类学的角度 ,判定集集地震是发生于板块边缘的大地震。可供研究台湾岛上地震活动性以及台湾地区板块构造演化问题参考。     

1999年9月21日台湾集集地震的震源参数综合判定

本文利用半经验、半理论的方法对世界各主要地震汇集和研究机构所给出的集集主震基本震源参数进行了综合计算和评定。根据这些参数重新计算了集集主震的标量地震矩、破裂面积、应力降等。根据这些参数之间的内在联系,从地震分类学的角度,判定集集地震是发生于板块边缘的大地震。可供研究台湾岛上地震活动性以及台湾地区板块构造演化问题参考。     

Kinematic analysis and palaeoseismology of the Edremit Fault Zone: evidence for past earthquakes in the southern branch of the North Anatolian Fault Zone,Biga Peninsula,NW Turkey

The Edremit Fault Zone (EFZ) forms one of the southern segments of the North Anatolian Fault Zone (NAFZ) at the northern margin of the Edremit Gulf (Biga Peninsula, South Marmara Region, Turkey). Stratigraphic, structural and kinematic results indicate that basinward younging of the fault zone, in terms of a rolling-hinge mechanism, has resulted in at least three discrete Miocene to Holocene deformational phases: the oldest one (Phase 1) directly related to the inactive Kazda? Detachment Fault, which was formed under N–S trending pure extension; Phase 2 is characterised by a strike-slip stress condition, probably related to the progression of the NAFZ towards the Edremit area in the Plio–Quaternary; and Phase 3 is represented by the high-angle normal faulting, which is directly interrelated with the last movement of the EFZ. Our palaeoseismic studies on the EFZ revealed the occurrence of three past surface rupture events; the first one occurred before 13178 BC, a penultimate event that may correspond to either the 160 AD or 253 AD historical earthquakes, and the youngest one can be associated with the 6 October 1944 earthquake (M_w = 6.8). These palaeoseismic data indicate that there is no systematic earthquake recurrence period on the EFZ.     

Palaeolimnology of Lake Sapanca and identification of historic earthquake signals,Northern Anatolian Fault Zone (Turkey)

Lake Sapanca is located on a strand of the Northern Anatolian Fault Zone (NAFZ, Turkey), where a series of strong earthquakes (Ms >6.0) have occurred over the past hundred years. Identifying prehistoric earthquakes in and around Lake Sapanca is key to a better understanding of plate movements along the NAFZ. This study contributes to the development of palaeolimnological tools to identify past earthquakes in Lake Sapanca. To this end several promising proxies were investigated, specifically lithology, magnetic susceptibility, grain size (thin-section and laser analysis), geochemistry, pollen concentration, diatom assemblages, ¹³⁷Cs and ²¹⁰Pb. Sedimentological indicators provided evidence for reworked, turbidite-like or homogeneous facies (event layers) in several short cores (<45 cm). Other indicators of sediment input and the historical chronicles available for the area suggest that three of these event layers likely originated from the AD 1957, 1967 and 1999 earthquakes. Recent changes in sediment deposition and nutrient levels have also been identified, but are probably not related to earthquakes. This study demonstrates that a combination of indicators can be used to recognize earthquake-related event layers in cores that encompass a longer period of time.     

Assessment of the statistical earthquake hazard parameters for NW Turkey

We performed a probabilistic analysis of earthquake hazard input parameters, NW Turkey covers Gelibolu and Biga Peninsulas, and its vicinity based on four seismic sub-zones. The number of earthquakes with magnitude M ≥ 3.0 occurred in this region for the period between 1912 and 2007 is around 5130. Four seismic source sub-zones were defined with respect to seismotectonic framework, seismicity and fault geometry. The hazard perceptibility characterization was examined for each seismic source zone and for the whole region. The probabilities of earthquake recurrences were obtained by using Poisson statistical distribution models. In order to determine the source zones where strong and destructive earthquakes may occur, distribution maps for a, b and a/b values were calculated. The hazard scaling parameters (generally known as a and b values) in the computed magnitude–frequency relations vary in the intervals 4.28–6.58 and 0.59–1.13, respectively, with a RMS error percentage below 10 %. The lowest b value is computed for sub-zone three indicating the predominance of large earthquakes mostly at Gelibolu (Gallipoli) and north of Biga Peninsula (southern Marmara region), and the highest b value is computed for sub-zone two Edremit Bay (SW Marmara region). According to the analysis of each seismic sub-zone, the greatest risk of earthquake occurrence is determined for the triangle of Gelibolu–Tekirda? western part of Marmara Sea. Earthquake occurrence of the largest magnitude with 7.3 within a 100-year period was determined to be 46 % according to the Poisson distribution, and the estimated recurrence period of years for this region is 50 ± 12. The seismic hazard is pronounced high in the region extending in a NW–SE direction, north of Edremit Bay, west of Saros Bay and Yenice Gönen (southern Marmara region) in the south. High b values are generally calculated at depths of 5–20 km that can be expressed as low seismic energy release and evaluated as the seismogenic zone.     

A seismic source zone model for the seismic hazard assessment of the Italian territory

We designed a new seismic source model for Italy to be used as an input for country-wide probabilistic seismic hazard assessment (PSHA) in the frame of the compilation of a new national reference map.

We started off by reviewing existing models available for Italy and for other European countries, then discussed the main open issues in the current practice of seismogenic zoning.

The new model, termed ZS9, is largely based on data collected in the past 10 years, including historical earthquakes and instrumental seismicity, active faults and their seismogenic potential, and seismotectonic evidence from recent earthquakes. This information allowed us to propose new interpretations for poorly understood areas where the new data are in conflict with assumptions made in designing the previous and widely used model ZS4.

ZS9 is made out of 36 zones where earthquakes with M_w > = 5 are expected. It also assumes that earthquakes with M_w up to 5 may occur anywhere outside the seismogenic zones, although the associated probability is rather low. Special care was taken to ensure that each zone sampled a large enough number of earthquakes so that we could compute reliable earthquake production rates.

Although it was drawn following criteria that are standard practice in PSHA, ZS9 is also innovative in that every zone is characterised also by its mean seismogenic depth (the depth of the crustal volume that will presumably release future earthquakes) and predominant focal mechanism (their most likely rupture mechanism). These properties were determined using instrumental data, and only in a limited number of cases we resorted to geologic constraints and expert judgment to cope with lack of data or conflicting indications. These attributes allow ZS9 to be used with more accurate regionalized depth-dependent attenuation relations, and are ultimately expected to increase significantly the reliability of seismic hazard estimates.     

Bayesian estimation of seismic hazard for two sites in Switzerland

Seismic hazard analysis is based on data and models, which both are imprecise and uncertain. Especially the interpretation of historical information into earthquake parameters, e.g. earthquake size and location, yields ambiguous and imprecise data. Models based on probability distributions have been developed in order to quantify and represent these uncertainties. Nevertheless, the majority of the procedures applied in seismic hazard assessment do not take into account these uncertainties, nor do they show the variance of the results. Therefore, a procedure based on Bayesian statistics was developed to estimate return periods for different ground motion intensities (MSK scale).Bayesian techniques provide a mathematical model to estimate the distribution of random variables in presence of uncertainties. The developed method estimates the probability distribution of the number of occurrences in a Poisson process described by the parameter . The input data are the historical occurrences of intensities for a particular site, represented by a discrete probability distribution for each earthquake. The calculation of these historical occurrences requires a careful preparation of all input parameters, i.e. a modelling of their uncertainties. The obtained results show that the variance of the recurrence rate is smaller in regions with higher seismic activity than in less active regions. It can also be demonstrated that long return periods cannot be estimated with confidence, because the time period of observation is too short. This indicates that the long return periods obtained by seismic source methods only reflects the delineated seismic sources and the chosen earthquake size distribution law.     

A newly developed seismic microzonation model of Erbaa (Tokat, Turkey) located on seismically active eastern segment of the North Anatolian Fault Zone (NAFZ)

A methodology to model seismic microzonation maps is required in the hazard mitigation decision plans of the earthquake prone areas. The stage of disaster preparedness for new residential places is of great importance for detailed seismic microzonation models. The effects of local geological and geotechnical site conditions were considered in order to establish site characterization as the initial stage of the models in this study. Dynamic soil properties based on the empirical correlations between shear wave velocity (V _s) and standard penetration test blow counts were taken into account in order to define representative soil profiles extending down to the engineering bedrock. One-dimensional site response analyses were performed to analyze earthquake characteristics on the ground surface. The layers for soil classification, geology, depth to groundwater level, amplification, distance to fault, slope and aspect, and liquefaction-induced ground deformation potential of the study area were prepared in seismic microzonation models. The study area, Erbaa, is placed along the seismically active North Anatolian Fault Zone. Final seismic microzonation map of the study area was evaluated applying different GIS-based Multi-Criteria Decision Analysis (MCDA) techniques. Two of the MCDA techniques, simple additive weighting and analytical hierarchical process (AHP), are considered during the evaluation step of the final seismic microzonation map. The comparison is made in order to distinguish two different maps based on these MCDA techniques. Eventually, AHP-based seismic microzonation map is more preferable for the seismic design purposes in this study.