Search

We computed a 3-D shear wave velocity model of the Marmara Sea region from ambient noise tomography. The correlations of up to 8 yr of vertical-component seismic recordings from 80 broad-band stations provided Rayleigh wave group velocity measurements in the period band 6–21 s at more than 1400 selected virtual source–receiver pairs. Rayleigh wave group velocity maps were used to derive a shear wave velocity model through simulated annealing inversion. The resulting crustal model provides coverage of the Marmara Sea along with its surrounding regional tectonic features. This allows for an investigation of the spatial extents of the Marmara Sea on a scale larger than that of basins. The low-velocity structures of the Marmara Sea and the Thrace Basins are coeval to a depth of approximately 9 km. The crustal velocities beneath the Marmara Sea basins exhibit a low vertical gradient and smooth horizontal variations. The regional tectonic structures, such as Istranca Massif, Istanbul and Sakarya Zones, display sharp velocity contrasts with the lower velocity crust beneath the Marmara Sea. The observed low crustal velocities, along with depth variations of the velocity isosurfaces (i.e. 3.4 km s−1) indicate that the Marmara region is a structural depression much deeper and larger than the three basins of the North Marmara Trough. The North Anatolian Fault Zone is unlikely to be the primary factor contributing to the origin of this significant depression, as the basin's development appears to have occurred before the fault propagated into the region. [ABSTRACT FROM AUTHOR]

The East Anatolian Fault (EAF) is a 700-km-long left-lateral transform fault located between the Anatolian and Arabian plates. The proximity of the Euler Pole to the Arabia–Anatolia Plate boundary leads to rapid changes in plate velocity along the boundary, which is manifested by the decreasing slip rates from east (10 mm yr–1) to west (∼1–4 mm yr–1). The EAF displays heterogeneous seismicity patterns with seismic gaps, localized clusters and broad diffuse zones. In this study, in order to understand the origin of these complexities and quantify the seismic hazard along the EAF, we present an improved seismicity catalogue with more than 26 000 earthquakes and 160 focal mechanisms from regional moment tensor inversion between 2007 and 2020. The focal mechanisms and seismicity show that the EAF dips towards north and forms a well-defined plate boundary in the east between Palu and Çelikhan with almost pure left-lateral motion. Further west, the boundary becomes broader with activity along subparallel faults. Focal mechanisms show heterogeneous stress orientations in consistence with geodetically determined strain rate field. The stress orientations show a transition from strike-slip to extension towards the west of Çelikhan. Amongst all segments of EAF, the Pütürge segment, which holds the near-repeating earthquakes in the vicinity of the nucleation of the 2020 M w 6.8 earthquake, is distinguished with its steady and high rate of seismicity. Further east, the neighbouring Palu segment is characterized by several distinct moderate earthquakes. We do not observe any change in the seismicity rate on these segments of the EAF following large earthquakes. In order to quantify the seismic hazard along the EAF, we calculate the recurrence time and maximum magnitude for each segment by using an extended seismicity catalogue of 150 yr including the large historical earthquakes and the geodetic strain rate. The results show ∼150 yr recurrence time with M max∼6.7–7.0 along the seismically active Palu and Pütürge segments on the east, while relatively silent western segments yield longer recurrence times; 237–772 for Pazarcık and 414–917 for Amanos segments with slightly larger magnitudes (M max ∼7–7.4). We infer that the seismicity patterns and strain-rate field along the EAF are shaped by several factors such as strong geometrical irregularities, heterogeneous coupling and complex plate motion leading to rapid change of fault slip rate. [ABSTRACT FROM AUTHOR]

In spite of growing evidence that many earthquakes are preceded by increased seismic activity, the nature of this activity is still poorly understood. Is it the result of a mostly random process related to the natural tendency of seismic events to cluster in time and space, in which case there is little hope to ever predict earthquakes? Or is it the sign that a physical process that will lead to the impending rupture has begun, in which case we should attempt to identify this process. With this aim we take a further look at the nucleation of two of the best recorded and documented strike-slip earthquakes to date, the 1999 Izmit and Düzce earthquakes which ruptured the North Anatolian Fault over ∼200 km. We show the existence of a remarkable mechanical logic linking together nucleation characteristics, stress loading, fault geometry and rupture speed. In both earthquakes the observations point to slow aseismic slip occurring near the ductile-to-brittle transition zone as the motor of their nucleation. [ABSTRACT FROM AUTHOR]

The unbroken section of the North Anatolian Fault beneath the Sea of Marmara is a major source of seismic hazard for the city of İstanbul. The northern and currently the most active branch, the Main Marmara Fault (MMF), is segmented within a shear zone and exhibits both partially creeping and locked behaviour along its 150 km length. In 2019 September, a seismic activity initiated near MMF, off-coast the town of Silivri, generating 14 earthquakes ≥  M w 3.5 in a week. The M w 5.8 Silivri earthquake, is the largest in the Marmara Sea since the 1963 M w 6.3 Çınarcık earthquake. Our analyses reveal that the activity started in a narrow zone (∼100 m) and spread to ∼7 km following an M w 4.7 foreshock within ∼2 d. The distribution of relocated aftershocks and the focal mechanisms computed from regional waveforms reveal that the M w 5.8 earthquake did not occur on the MMF, but it ruptured ∼60° north-dipping oblique strike-slip fault with significant thrust component located on the north of the MMF. Finite-fault slip model of the main shock shows 8 km long rupture with directivity toward east, where the ruptured fault merges to the MMF. The narrow depth range of the slip distribution (10–13 km) and the aftershock zone imply that the causative fault is below the deep sedimentary cover of the Marmara Basin. The distribution of aftershocks of the M w 5.8 event is consistent with Coulomb stress increase. The stress changes along MMF include zones of both stress decrease due to clamping and right-lateral slip, and stress increase due to loading. [ABSTRACT FROM AUTHOR]

The 2017 July 20, M w6.6 Bodrum–Kos earthquake occurred in the Gulf of Gökova in the SE Aegean, a region characterized by N–S extension in the backarc of the easternmost Hellenic Trench. The dip direction of the fault that ruptured during the earthquake has been a matter of controversy where both north- and south-dipping fault planes were used to model the coseismic slip in previous studies. Here, we use seismic (seismicity, main shock modelling, aftershock relocations and aftershock mechanisms using regional body and surface waves), geodetic (GPS, InSAR) and structural observations to estimate the location, and the dip direction of the fault that ruptured during the 2017 earthquake, and the relationship of this event to regional tectonics. We consider both dip directions and systematically search for the best-fitting locations for the north- and south-dipping fault planes. Comparing the best-fitting planes for both dip directions in terms of their misfit to the geodetic data, proximity to the hypocenter location and Coulomb stress changes at the aftershock locations, we conclude that the 2017 earthquake ruptured a north-dipping fault. We find that the earthquake occurred on a 20–25 km long, ∼E–W striking, 40° north-dipping, pure normal fault with slip primarily confined between 3 and 15 km depth, and the largest slip exceeding 2 m between depths of 4 and 10 km. The coseismic fault, not mapped previously, projects to the surface within the western Gulf, and partly serves both to widen the Gulf and separate Kos Island from the Bodrum Peninsula of SW Anatolia. The coseismic fault may be an extension of a mapped, north-dipping normal fault along the south side of the Gulf of Gökova. While all of the larger aftershocks are consistent with N–S extension, their spatially dispersed pattern attests to the high degree of crustal fracturing within the basin, due to rapid trenchward extension and anticlockwise rotation within the southeastern Aegean. [ABSTRACT FROM AUTHOR]

The engines of surface deformation in the Anatolia–Aegean region are a matter of debate, including the origin of the high elevations of the Anatolian plateau. Recent publications based on geological and thermomechanical modelling emphasize the role of dynamic topography in the plateau uplift. However, quantitative estimates of the contribution of dynamic topography are affected by large uncertainties due to insufficient knowledge of the crustal structure, in particular crustal thickness and density. To reduce these uncertainties, we provide a new accurate crustal thickness map of the Anatolia–Aegean domain computed from a large volume of broadband seismic data. In addition, we display high-resolution seismic sections of the internal structure of the crust in Western and Central Anatolia. Density contrasts are derived from the same seismic data set and Bouguer gravity anomaly computed from the EGM2008 model. Our crustal thickness model is highly correlated with the topography suggesting that the Anatolian plateau is close to isostatic equilibrium. The average density difference between crust and upper mantle computed from our crustal model and Bouguer gravity anomaly is low compared to the global average, ∼0.315 ×103 kg m−3. The ratio of surface elevation to crustal thickness is lower than average, 1/9.4, which also indicates a low-density crust. Differences between isostatic topography and observed topography are overall small (<500 m). The east-to-west gradients of crustal thickness and topography changes are nearly constant in between the Taurides and Pontides at the northern and southern borders of Anatolia. The observed constant crustal thickness gradient may indicate a low viscosity lower crust supported by the thin mantle lithosphere evidenced by seismic tomography beneath the Anatolian plateau. We propose that viscous flow in the lower crust has smoothed out lateral changes in the crustal structure expected for such a heterogeneous collage of continental fragments. This flow may originate from gravitational potential energy differences between Eastern Anatolia (thick crust, high elevations) and the Aegean Sea (thin crust, low elevations), suggesting that gravity plays an integral part in the westward escape of Anatolia. [ABSTRACT FROM AUTHOR]

We report that asperities with the highest coseismic slip in the 2014 MW6.9 North Aegean earthquake persisted through the interseismic, coseismic and immediate post-seismic periods. We use GPS and seismic data to obtain the source model of the 2014 earthquake, which is located on the western extension of the North Anatolian Fault (NAF). The earthquake ruptured a bilateral, 90 km strike-slip fault with three slip patches: one asperity located west of the hypocentre and two to the east with a rupture duration of 40 s. Relocated pre-earthquake seismicity and aftershocks show that zones with significant coseismic slip were relatively quiet during both the 7 yr of interseismic and the 3-month aftershock periods, while the surrounding regions generated significant seismicity during both the interseismic and post-seismic periods. We interpret the unusually long fault length and source duration, and distribution of preand post-main-shock seismicity as evidence for a rupture of asperities that persisted through strain accumulation and coseismic strain release in a partially coupled fault zone. We further suggest that the association of seismicity with fault creep may characterize the adjacent Izmit, Marmara Sea and Saros segments of the NAF. Similar behaviour has been reported for sections of the San Andreas Fault, and some large subduction zones, suggesting that the association of seismicity with creeping fault segments and rapid relocking of asperities may characterize many large earthquake faults. [ABSTRACT FROM AUTHOR]

Western Turkey provides spectacular examples of the two end-member models of deformation of the continental lithosphere, with strain localization along the North Anatolian fault and distributed north–south extension along the Aegean coast. To provide constraints on the mechanisms of continental deformation, we present a new high-resolution image of lithospheric structure along a ∼650 km transect crossing Western Anatolia at 28°E longitude from the Black Sea to the Mediterranean. More than 2600 receiver functions are computed from records of teleseismic earthquakes at 40 broadband seismic stations with an average spacing of 15 km. Lateral variations of crustal thickness and Vp/Vs ratio are inferred from both H-k and common conversion point stacks. We observe long-wavelength variations of Moho depth from ∼31 km in the Thrace basin to ∼25 km beneath the Marmara Sea, ∼32 km beneath the Izmir–Ankara suture, ∼25 km beneath the Menderes Massif and ∼20 km on the coast of the Mediterranean. No mid- or lower-crustal interface is visible in the migrated depth section and seismic discontinuities are confined to the shallow crust. The small-amplitude and long-wavelength lateral variations of the Moho topography suggest that viscous flow in a hot lower crust has smoothed out the lateral variations of crustal thickness induced by Cenozoic continent–continent collision. The crust–mantle boundary is flat beneath the central and southern Menderes Massif. The rougher Moho topography and more heterogeneous crust of the Marmara region likely result from transtension in the North Anatolian Fault Zone superimposed onto Aegean extension. The Moho of the subducted African lithosphere is observed dipping northward between ∼40 and ∼60 km depth at the southern end of the profile. The abrupt termination of the subducted slab only 50 km to the north of the Mediterranean coast confirms the slab tear inferred from previous tomographic studies. [ABSTRACT FROM AUTHOR]

SUMMARY This study provides new constraints on the upper-mantle structure from western Greece to central Anatolia using seismic data of permanent broad-band networks recently installed in Greece and Turkey and from a two-year temporary array (SIMBAAD experiment). We used ∼200 seismic events recorded at 146 broad-band stations with a typical interstation distance of 60-100 km across the study area. The high-resolution 3-D shear wave velocity model of the mantle is obtained by inversion of fundamental-mode Rayleigh wave phase velocity maps for periods between 20 and 195 s. The tomography is based on ray tracing in heterogeneous media taking into account external propagation effects. The horizontal resolution is approximately 100 km, however small heterogeneities may suffer from some horizontal smearing and damping. The vertical resolution is approximately 100 km. The vertical smoothing is necessary to avoid unresolved spurious shear wave velocity oscillations in the upper mantle. The errors on shear wave velocities in our 3-D model (0.02-0.1 km s−1) are significantly smaller than the amplitude of Vs variations (0.3-0.5 km s−1). In spite of the vertical and horizontal smoothing, our model shows details in the upper-mantle structure never reached at regional scale in the area. The overall structure is characterized by a low-velocity zone (80-200 km depth) reflecting a slow and warm asthenosphere underlying a thin lithosphere. The southwesternmost termination of the low-velocity anomaly corresponds to the northward dipping Hellenic slab. The detailed shear velocity structure of the upper mantle beneath Anatolia appears to be far more geometrically complex than revealed in previous tomographic studies of the area. At depths larger than or equal to 160 km, velocities are overall high beneath Anatolia, partly due to the presence of dipping high-velocity anomalies which we tentatively interpret as remnant slabs. The southernmost high-velocity anomaly beneath Anatolia is separated from the eastern edge of the Hellenic slab by a major low-velocity anomaly which we interpret as the trace of asthenospheric mantle material rising inside a vertical slab tear beneath southwestern Anatolia. [ABSTRACT FROM AUTHOR]

SUMMARY This study presents an analysis of Pn traveltimes to determine lateral variations of velocity, anisotropy of uppermost mantle and crustal thickness beneath Turkey and adjacent regions. From 1999 to 2010, more than 50 000 Pn arrivals are compiled from 700 regional earthquakes and 832 stations of permanent and temporary networks operated in the study area. We used a regularized least-squares inversion to estimate the uppermost mantle parameters. The results reveal features that correlate well with the geology and the active tectonics of the region. Overall, the Pn velocities show very fast (>8.4 km s-1) and very slow (<7.6 km s-1) anomalies indicating a heterogeneous lithospheric structure. Relatively uniform Pn velocities (7.9-8.1 km s-1) are observed in western Turkey. The lowest velocities coincide with the volcanics of the easternmost Anatolia and the Central Anatolian Volcanic Zone. High Pn velocities are observed at Mediterranean Basin (>8.3 km s-1), western Black Sea Basin (>8.3 km s-1), Adriatic Sea (>8.3 km s-1) and Zagros Suture (>8.3 km s-1). Pn anisotropy has maximum amplitude of ±0.6 km s-1 in the study area. The largest and coherent anisotropic anomalies are observed in the western Anatolia, Aegean Sea and Cyprian Arc. Pn anisotropy in western Anatolia, Aegean Sea and Greece correlate well with the present deformation of N-S regional extension. No clear correlation between anisotropy and the shear deformation along the North Anatolian Fault was observed. The majority of the stations in central Anatolia show small station residuals indicating the average crustal thickness of 35 ± 2 km. Western Anatolia and the Aegean Sea have crustal thicknesses between 28 ± 2 and 33 ± 2 km. In Greece, the crustal thicknesses are increasing from 33 ± 3 km from the western coast to a maximum of 48 ± 3 km beneath Dinarides-Hellenides. The large crustal thicknesses are also observed along southern coast of Anatolia (40-48 km). [ABSTRACT FROM AUTHOR]

We present the analysis of a Love-wave dispersion study performed across Turkey and the surroundings. Group velocity dispersion curves were obtained from the local and regional earthquakes recorded at both permanent and temporary stations operated in the region. Love-wave group velocity maps in the period range of 10–50 s were computed using a tomographic inversion method. The study shows the existence of significantly different crustal types in the area. The group velocity anomalies are associated with the major geological structures in the region. Low group velocities at shorter periods (10–20 s) are observed in the local sedimentary basins, the eastern Mediterranean and the Black Sea. The eastern Anatolia region is also characterized by low group velocities while Pontides and Bitlis–Pötürge massif display higher group velocities. The central Anatolia exhibits uniform velocity distribution indicating more homogenous crust. The Isparta angle is marked by a wedge-shaped low group velocity anomaly. The low velocities observed in the Isparta angle are related to crustal thickening and subduction-related complex. High velocities observed on the maps are associated with metamorphics, magmatic arcs along the orogenic belts of Pontides, Pötürge massif and crustal thinning in the Aegean region. At larger periods (40–50 s) the Anatolian block shows low and uniform group velocity distribution while its surroundings display higher group velocities with the exception of the eastern Mediterranean region. [ABSTRACT FROM AUTHOR]

Three main shocks M-1, M-2 and M-3 (17 October 2005 at 05:45 UTC, Mw 5.4; 17 October at 09:46 UTC, Mw 5.8 and 20 October at 21:40 UTC, Mw 5.9) and their associated aftershocks within the Gulf of S iğac ik, 50 km southwest of Izmir, Turkey were studied in detail. A temporary seismic network deployed during the activity allowed the hypocentre of M-3 and subsequent aftershocks to be determined with high accuracy. A relative relocation technique was used to improve the epicentres of M-1 and M-2. All three main shocks have strike-slip mechanisms which agree with the linear trends of the aftershock locations. Two distinct zones were illuminated by the aftershock locations. The zones contain clear echelon patterns with slightly different orientations from the trend of the aftershock distribution. M-2 and M-3 ruptured along of the eastern rupture zone which aligns N45°E. However the strike direction of M-1 is not clearly identified. The alignment of the two rupture zones intersect at their southern terminus at an angle of 90°. The fault zones form conjugate pair system and static triggering is considered as a probable mechanism for the sequential west to east occurrence of M-1, M-2 and M-3. This earthquake sequence provides seismological evidence for conjugate strike-slip faulting co-existing within a region dominated by north–south extension and well-developed east–west trending normal faults. [ABSTRACT FROM AUTHOR]

We present observations of ground accelerations recorded at a small array close to the fault during the Düzce earthquake and its early aftershocks. The records show the strong spatial variability of ground acceleration over distances of only a few hundred metres. During the main shock, the peak horizontal acceleration values ranged from 0.3 to about 1.0 g at stations distant of 1.5 km only. We attribute this spatial variability to a fault zone site effect as peak ground acceleration steadily increases as the distance to the fault trace decreases. The spectral ratio between the ground motion recorded near the fault and the one outside the fault zone shows a shift of the spectral peak to lower frequencies with increasing peak accelerations. Such an observation suggests a non-linear behaviour of the fault zone due to the strong ground shaking. As much as a 45 per cent reduction in the shear wave velocity is necessary for the observed shifts. The opening of pre-existing cracks throughout the fault zone is the proposed mechanism to account for the observed shear wave reductions. The observation that elastic fault zone properties are soon recovered following episodes of large strains shows that cracks and fissures close rapidly after the strong shaking is over. [ABSTRACT FROM AUTHOR]

Focuses on the use of an indirect boundary element method in computing the seismic wavefields in layered media. Advantages of the technique; Suitability of the technique for models with high aspect ratios and a limited number of layers; Features of the technique.