Analysis of Izmit aftershocks 25 days before the November 12th 1999 Düzce earthquake, Turkey

We investigate spatial clustering of 2414 aftershocks along the Izmit M_w = 7.4 August 17, 1999 earthquake rupture zone. 25 days prior to the Düzce earthquake M_w = 7.2 (November 12, 1999), we analyze two spatial clusters, namely Sakarya (SC) and Karadere–Düzce (KDC). We determine the earthquake frequency–magnitude distribution (b-value) for both clusters. We find two high b-value zones in SC and one high b-value zone in KDC which are in agreement with large coseismic surface displacements along the Izmit rupture. The b-values are significantly lower at the eastern end of the Izmit rupture where the Düzce mainshock occurred. These low b-values at depth are correlated with low postseismic slip rate and positive Coloumb stress change along KDC. Since low b-values are hypothesized with high stress levels, we propose that at the depth of the Düzce hypocenter (12.5 km), earthquakes are triggered at higher stresses compared to shallower crustal earthquake. The decrease in b-value from the Karadere segment towards the Düzce Basin supports this low b-value high stress hypothesis at the eastern end of the Izmit rupture. Consequently, we detect three asperity regions which are correlated with high b-value zones along the Izmit rupture. According to aftershock distribution the half of the Düzce fault segment was active before the 12 November 1999 Düzce mainshock. This part is correlated with low b-values which mean high stress concentration in the Düzce Basin. This high density aftershock activity presumably helped to trigger the Düzce event (M_w = 7.2) after the Izmit M_w 7.4 mainshock.     

Heterogeneous response of hydrogeological systems to the Izmit and Düzce (Turkey) earthquakes of 1999

At four sites in Turkey and Armenia the physico-chemical properties of thermal and mineral waters were monitored continuously during the Izmit and Düzce earthquakes that occurred along the North Anatolian fault in August and November 1999. The epicentral distances between the moment magnitude (M_w) 7.6 Izmit earthquake and the monitoring locations were 313, 488, 1,161, and 1,395 km. At the most distant site, the specific electrical conductivity of mineral water from a flowing artesian well dropped co-seismically and postseismically by 7%. No changes were observed at the other sites, although the estimated earthquake strains and peak ground accelerations are much higher. A similar pattern was observed after the Düzce earthquake, which happened three months after the Izmit event. The response of a hydrogeological system seems to depend on the site characteristics rather than on the nature of the earthquake. A hydrogeological model for the sensitive observation site farthest from the Izmit earthquake explains the observations in terms of a changed mixing ratio between two fluid components. Passing seismic waves may trigger a local pore-pressure increase according to the mechanism of advective overpressure. The preconditions for this mechanism, free gas bubbles in the aquifer in combination with a trap for rising bubbles, is probably not fulfilled by the other groundwater systems. Electronic Publication     

The rupture process during the 1999 Düzce, Turkey, earthquake from joint inversion of teleseismic and strong-motion data

Teleseismic and strong-motion data are inverted to determine the rupture process during the November 1999 Düzce earthquake in NW Turkey. The fault geometry, rise time and rupture velocity are determined from the aftershock distribution and preliminary inversions of the teleseismic data. Joint inversion of the teleseismic and strong-motion data is then carried out for the slip distribution. We obtain the strike 264°, dip 64°, rake −172°, seismic moment 5.0×10¹⁹ N m (M_w 7.1), and average stress drop 7 MPa. This earthquake was characterized by bilateral fault rupture and asymmetric slip distribution. Two asperities (areas of large slip) are identified, the eastern one being 1.5 times larger than the western one. The derived slip distribution is consistent with the aftershock distribution, surface rupture and damage. The point of rupture initiation in this Düzce earthquake coincided with the eastern tip of the aftershock distribution of the August 1999 Izmit earthquake.     

Characteristics of moderate- to large-scale landslides triggered by the <Emphasis Type="Italic">M</Emphasis><Subscript><Emphasis Type="Italic">w</Emphasis></Subscript> 7.8 2015 Gorkha earthquake and its aftershocks

The M _w 7.8 2015 Gorkha earthquake and its aftershocks significantly impacted the lives and economy of Nepal. The consequences of landslides included fatalities, property losses, blockades of river flow, and damage to infrastructural systems. Co-seismic landslides triggered by this earthquake were significantly widespread and pose a major geodisaster. There were tens of thousands of landslides triggered by the earthquake, majority of which were distributed in between the epicenter of the main shock and the M _w 7.3 aftershock. Although 14,670 landslides triggered by this earthquake were identified, only approximately 23% of them were of moderate to large scale with areas greater than 100 m². Of the moderate- to large-scale landslides identified, just over 90% were triggered by the main shock and smaller aftershocks prior to the major (M _w 7.3) aftershock, while nearly 10% were triggered by the ground shaking induced by the major aftershock. Moreover, the number of landslides triggered by the 2015 Gorkha earthquake, specifically by the main shock, was slightly more than the expected number of landslides for the recorded maximum peak ground acceleration (PGA) in comparison to the co-seismic landslides triggered by 26 earthquakes. Over 90% of those moderate- to large-scale landslides were concentrated within the estimated fault rupture surface. Majority of these moderate- to large-scale landslides were disrupted failures with over 96% of which were classified as earth falls. However, the majority of small-scale landslides were rock or boulder falls. The most number of moderate- to large-scale landslides were triggered in the slate, shale, siltstone, phyllite, and schist of the Lesser Himalayan formation followed by an equally significant number in both schist, gneiss, etc. of the Higher Himalayan formation and the phyllite, metasandstone, schist, etc. of the Lesser Himalayan formation. The sizes (i.e., areas) of the landslides were lognormally distributed, with a mode area of 322.0 m². Slope inclinations of the moderate- to large-scale landslides followed a normal distribution with a mean slope inclination of 32.6° and standard deviation of 13.5°. There exists a strong correlation between the number of landslides and the peak ground acceleration within the study area, specific for different geological formations.     

Hybrid attenuation model for estimation of peak ground accelerations in the Kutch region, India

The Kutch region of Gujarat in India is the locale of one of the most devastating earthquake of magnitude (M _w) 7.7, which occurred on January 26, 2001. Though, the region is considered as seismically active region, very few strong motion records are available in this region. First part of this paper uses available data of strong motion earthquakes recorded in this region between 2006 and 2008 years to prepare attenuation relation. The developed attenuation relation is further used to prepare synthetic strong motion records of large magnitude earthquakes using semiempirical simulation technique. Semiempirical simulation technique uses attenuation relation to simulate strong ground motion records of any target earthquake. The database of peak ground acceleration obtained from simulated records is used together with database of peak ground acceleration obtained from observed record to develop following hybrid attenuation model of wide applicability in the Kutch region: $$ \begin{aligned} \ln \left( {\text{PGA}} \right) & = - 2.56 + 1.17 \, M_{\text{w}} - \, 0.015R - 0.0001\ln \left( {E + 15} \right) \\ &\quad 3.0 \le M_{\text{w}} \le 8.2;\quad 12 \le R \le 120;\quad {\text{std}} . {\text{ dev}}.(\sigma ): \pm 0.5 \\ \end{aligned} $$ ln ( PGA ) = ? 2.56 + 1.17 M w ? 0.015 R ? 0.0001 ln ( E + 15 ) 3.0 ≤ M w ≤ 8.2 ; 12 ≤ R ≤ 120 ; std . dev . ( σ ) : ± 0.5 In the above equation, PGA is maximum horizontal ground acceleration in gal, M _w is moment magnitude of earthquake, R is hypocentral distance, and E is epicentral distance in km. The standard deviation of residual of error in this relation is 0.5. This relation is compared with other available relations in this region, and it is seen that developed relation gives minimum root mean square error in comparison with observed and calculated peak ground acceleration from same data set. The applicability of developed relation is further checked by testing it with the observed peak ground acceleration from earthquakes of magnitude (M _w), 3.6, 4.0, 4.4, and 7.7, respectively, which are not included in the database used for regression analysis. The comparison demonstrates the efficacy of developed hybrid attenuation model for calculating peak ground acceleration values in the Kutch region.     

Modeling of strong motion generation areas of the 2011 Tohoku,Japan earthquake using modified semi-empirical technique

Modification in the semi-empirical technique for the simulation of strong ground motion has been introduced to incorporate the strong motion generation areas (SMGA) in the modeled rupture plane. Strong motion generation areas identified within the rupture plane of the Tohoku earthquake of March 11, 2011 (M _w = 9.0), have been modeled using this modified technique. Two different source models having four and five SMGAs, respectively, are considered for modeling purpose. Strong motion records using modified semi-empirical technique have been simulated at two near-field stations located at epicentral distance of 137 and 140 km, respectively, using two different source models. Comparison of the observed and simulated acceleration waveforms is made in terms of root mean square error (RMSE) at both stations. Minimum root mean square error of the waveform comparison has been obtained at both the stations for source model having five SMGAs. Simulations from same rupture model have been made at other four stations lying at epicentral distance between 154 and 249 km. Comparison of observed and simulated records has been made in terms of RMSE in acceleration records, velocity records and response spectra at each six station. Simulations have been made at six other stations to obtain distribution of peak ground acceleration and peak ground velocity with hypocentral distance. Peak ground acceleration and velocity from simulated and observed records are compared at twelve stations surrounding the source of Tohoku earthquake. Comparison of waveforms and parameters extracted from observed and simulated strong motion records confirms the efficacy of the developed modified technique to model earthquake characterized by SMGAs.     

Seismic hazard across Bulgaria and neighbouring areas: regional and site-specific maximum credible magnitudes and earthquake perceptibility

A probabilistic seismic hazard assessment is developed here using maximum credible earthquake magnitude statistics and earthquake perceptibility hazard. Earthquake perceptibility hazard is defined as the probability a site perceives ground shaking equal to or greater than a selected ground motion level X, resulting from an earthquake of magnitude M, and develops estimates for the most perceptible earthquake magnitude, M _P(max). Realistic and usable maximum magnitude statistics are obtained from both whole process and part process statistical recurrence models. These approaches are extended to develop relationships between perceptible earthquake magnitude hazard and maximum magnitude recurrence models that are governed by asymptotic and finite return period properties, respectively. Integrated perceptibility curves illustrating the probability of a specific level of perceptible ground motion due to all earthquakes over the magnitude range extending from ?∞ to a magnitude M _i are then developed from reviewing site-specific magnitude perceptibility. These lead on to achieving site-specific annual probability of exceedance hazard curves for the example cities of Sofia and Thessaloniki for both horizontal ground acceleration and ground velocity. Both the maximum credible earthquake magnitude M ₃ and the most perceptible earthquake magnitude M _P(max) are of importance to the earthquake engineer when approaching anti-seismic building design. Both forms of hazard are illustrated using contoured hazard maps for the region bounded by 39°–45°N, 19°–29°E. Patterns are observed for these magnitude hazard estimates—especially M _P(max) specific to horizontal ground acceleration and horizontal ground velocity—and compared to inferred patterns of crustal deformation across the region. The full geographic region considered is estimated to be subject to a maximum credible earthquake magnitude M ₃—estimated using cumulative seismic moment release statistics—of 7.53 M _w, calculated from the full content of the adopted earthquake catalogue, while Bulgaria’s capital, Sofia, is estimated a comparable value of 7.36 M _w. Sofia is also forecast most perceptible earthquake magnitudes for the lowest levels considered for horizontal ground acceleration of M _PA(50) = 7.20 M _w and horizontal ground velocity of M _PV(5) = 7.23 M _w for a specimen focal depth of 15 km.     

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.     

Deterministic seismic microzonation of Kolkata city

This paper presents the deterministic seismic microzonation of densely populated Kolkata city situated on the world’s largest delta island with very soft and thick soil deposit in the surficial layers. A fourth-order accurate staggered-grid finite-difference algorithm for SH-wave propagation simulation in visco-elastic medium is used for the linear computation of ground motion amplifications in sedimentary deposit. Different maps such as for fundamental frequency (F ₀), peak ground acceleration (PGA), peak ground velocity, and peak ground displacement are developed for variety of end-users communities, including structural and geotechnical engineers for performance-based designs, building officials, emergency managers, land-use planners, private businesses, and the general public. The scenario of simulated amplification factors in the different frequency bands revealed that the Kolkata city is very much prone to severe damage even during a moderate earthquake and very selective damage may occur at some of the localities during local and distant earthquakes. The deterministically predicted PGA at bedrock level is 0.0844 g and the maximum PGA predicted at the free surface is 0.6 g in Kolkata city due to maximum credible earthquake (M _w = 5.4) associated with Eocene Hinge Zone at a depth of 36 km. The seismic microzonation of Kolkata city reveals that the Nager Bazar and Nimtala areas are the safest regions with earthquake point of view.     

Andean earthquakes triggered by the 2010 Maule,Chile (Mw 8.8) earthquake: Comparisons of geodetic,seismic and geologic constraints

The Maule, Chile, (M_w 8.8) earthquake on 27 February 2010 triggered deformation events over a broad area, allowing investigation of stress redistribution within the upper crust following a mega-thrust subduction event. We explore the role that the Maule earthquake may have played in triggering shallow earthquakes in northwestern Argentina and Chile. We investigate observed ground deformation associated with the M_w 6.2 (GCMT) Salta (1450 km from the Maule hypocenter, 9 h after the Maule earthquake), M_w 5.8 Catamarca (1400 km; nine days), M_w 5.1 Mendoza (350 km; between one to five days) earthquakes, as well as eight additional earthquakes without an observed geodetic signal. We use seismic and Interferometric Synthetic Aperture Radar (InSAR) observations to characterize earthquake location, magnitude and focal mechanism, and characterize how the non-stationary, spatially correlated noise present in the geodetic imagery affects the accuracy of our parameter estimates. The focal mechanisms for the far-field Salta and Catamarca earthquakes are broadly consistent with regional late Cenozoic fault kinematics. We infer that dynamic stresses due to the passage of seismic waves associated with the Maule earthquake likely brought the Salta and Catamarca regions closer to failure but that the involved faults may have already been at a relatively advanced stage of their seismic cycle. The near-field Mendoza earthquake geometry is consistent with triggering related to positive static Coulomb stress changes due to the Maule earthquake but is also aligned with the South America-Nazca shortening direction. None of the earthquakes considered in this study require that the Maule earthquake reactivated faults in a sense that is inconsistent with their long-term behavior.     

City of Bingöl in May 2003: Assessment of strong ground motion records

The main shock of Bingöl earthquake (M _W = 6.4) recorded by six accelerometers in the area occurred at 03:27 local time on May 1, 2003. The largest acceleration value of north–south component was recorded as 545.5 cm/s² at the nearest station which it is 12 km away from the epicenter of earthquake. Especially, 0.15 s short period was observed when high spectral acceleration value occurred. An acceleration value greater than 50 gal was recorded at the BNG (Bingöl) station and structural damage occurred within 6.5 s was very important for the near source and strong ground motion seismology. The recorded peak acceleration values were greater than the estimated empirical acceleration values. However, the structural damage was not as high and widespread as expected. This occurrence was explained by considering the factors of earthquake source, frequency content, effective duration, effective acceleration value, local soil conditions, rupture direction and attenuation.     

Effect of frequency-dependent radiation pattern in the strong motion simulation of the 2011 Tohoku earthquake,Japan, using modified semi-empirical method

We perform a strong ground motion simulation using a modified semi-empirical technique (Midorikawa in Tectonophysics 218:287–295, 1993), with frequency-dependent radiation pattern model. Joshi et al. (Nat Hazards 71:587–609, 2014) have modified the semi-empirical technique to incorporate the modeling of strong motion generation areas (SMGAs). A frequency-dependent radiation pattern model is applied to simulate high-frequency ground motion more precisely. Identified SMGAs (Kurahashi and Irikura in Earth Planets Space 63:571–576, 2011) of the 2011 off the Pacific coast of Tohoku earthquake (M _w  = 9.0) were modeled using this modified technique. We analyzed the effect of changing seismic moment values of SMGAs on the simulated acceleration time series. Final selection of the moment values of SMGAs is based on the root-mean-square error (RMSE) of waveform comparison. Records are simulated for both frequency-dependent and constant radiation pattern function. Simulated records for both cases are compared with observed records in terms of peak ground acceleration, peak ground velocity and pseudo-acceleration response spectra at different stations. Comparison of simulated and observed records in terms of RMSE suggests that the method is capable of simulating record, which matches in a wide frequency range for this earthquake and bears realistic appearance in terms of shape and strong motion parameters. The results confirm the efficacy and suitability of rupture model defined by five SMGAs for the developed modified technique.     

Near-field and far-field simulation of accelerograms of Sikkim earthquake of September 18, 2011 using modified semi-empirical approach

The semi-empirical approach for modeling of strong ground motion given by Midorikawa (Tectonophysics 218:287?C295, 1993) has been modified in the present paper for component wise simulation of strong ground motion. The modified approach uses seismic moment in place of attenuation relation for scaling of acceleration envelope. Various strong motion properties like directivity effect and dependence of peak ground acceleration with respect to surface projection of source model have been studied in detail in the present work. Recently, Sikkim earthquake of magnitude 6.9 (M _w ) that occurred on September 18, 2011 has been recorded at various near-field and far-field strong motion stations. The modified semi-empirical technique has been used to confirm the location and parameters of rupture responsible for this earthquake. Strong motion record obtained from the iterative modeling of the rupture plane has been compared with available strong motion records from near as well as far-field stations in terms of root mean square error between observed and simulated records. Several possibilities of nucleation point, rupture velocity, and dip of rupture plane have been considered in the present work and records have been simulated at near-field stations. Final selection of model parameters is based on root mean square error of waveform comparison. Final model confirms southward propagating rupture. Simulations at three near-field and twelve far-field stations have been made using final model. Comparison of simulated and observed record has been made in terms of peak ground acceleration and response spectra at 5?% damping. Comparison of simulated and observed record suggests that the method is capable of simulating record which bears realistic appearance in terms of shape and strong motion parameters. Present work shows that this technique gives records which matches in a wide frequency range for Sikkim earthquake and that too from simple and easily accessible parameters of the rupture plane.     

Spatial characteristics of landslides triggered by the 2015 M<Subscript>w</Subscript> 7.8 (Gorkha) and M<Subscript>w</Subscript> 7.3 (Dolakha) earthquakes in Nepal

High magnitude earthquakes trigger numerous landslides and their occurrences are mainly controlled by terrain parameters. We created an inventory of 15,551 landslides with a total area of 90.2 km² triggered by the 2015 M_w 7.8 (Gorkha) and M_w 7.3 (Dolakha) earthquakes in Nepal, through interpretation of very high resolution satellite images (e.g. WorldView, Pleiades, Cartosat-1 and 2, Resourcesat-2). Our spatial analysis of landslide occurrences with ground acceleration, slope, lithology and surface defomation indicated ubiquitous control of steep slope on landslides with ground acceleration as the trigger. Spatial distribution of landslides shows increasing frequency away from the Gorkha earthquake epicentre up to 130 km towards east, dropping sharply thereafter, which is an abnormal phenomenon of coseismic landslides. Landslides are laterally concentrated in three zones which matches well with the seismic rupture evolution of Gorkha earthquake, as reported through teleseismic measurements.     

Triggered/migrated seismicity due to the 2001 M w7.7 Bhuj earthquake, Western India

Paper describes triggered seismicity to 200?km distance and for a decade due to the 2001 M _w7.7 Bhuj earthquake. The Kachchh region is seismically one of the most active intraplate regions of the World due to the occurrence of two large earthquakes 1819 (M _w7.8) and 2001 (M _w7.7). Though, it has high hazard but was known to have low seismicity in view of the occurrence of fewer smaller shocks. However, the status seems to have changed after 2001. Besides the strong aftershock activity for over a decade, seismicity has spread to nearby faults in Kachchh peninsula and at several places southward for 200?km distance in Saurashtra peninsula. Beyond the rupture zone of the 2001 Bhuj earthquake, more than 40 mainshocks of M _w?~?3?C5 have occurred at 20 different locations, which is unusual. The increased seismicity is inferred to be caused by stress perturbation due to the 2001 Bhuj earthquake by viscoelastic process. In Saurashtra, over and above the viscoelastic stress increase, the transient stress increase by water table rise in monsoons seems to be affecting the timing of mainshocks and associated sequences of earthquakes.     

Aftershock sequences of two great Sumatran earthquakes of 2004 and 2005 and simulation of the minor tsunami generated on September 12, 2007 in the Indian Ocean and its effect

We present the seismic energy, strain energy, frequency–magnitude relation (b-value) and decay rate of aftershocks (p-value) for the aftershock sequences of the Andaman–Sumatra earthquakes of December 26, 2004 (M _w 9.3) and March 28, 2005 (M _w 8.7). The energy released in aftershocks of 2004 and 2005 earthquake was 0.135 and 0.365% of the energy of the respective mainshocks, while the strain release in aftershocks was 39 and 71% for the two earthquakes, respectively. The b-value and p-value indicate normal value of about 1. All these parameters are in normal range and indicate normal stress patterns and mechanical properties of the medium. Only the strain release in aftershocks was considerable. The fourth largest earthquake in this region since 2004 occurred in September 2007 off the southern coast of Island of Sumatra, generating a relatively minor tsunami as indicated by sea level gauges. The maximum wave amplitude as registered by the Padang, tide gauge, north of the earthquake epicenter was about 60 cm. TUNAMI-N2 model was used to investigate ability of the model to capture the minor tsunami and its effect on the eastern Indian Coast. A close comparison of the observed and simulated tsunami generation, propagation and wave height at tide gauge locations showed that the model was able to capture the minor tsunami phases. The directivity map shows that the maximum tsunami energy was in the southwest direction from the strike of the fault. Since the path of the tsunami for Indian coastlines is oblique, there were no impacts along the Indian coastlines except near the coast of epicentral region.     

The Gornozavodsk earthquake of August 17(18), 2006, in the south of Sakhalin Island

The August 17 (18), 2006, Gornozavodsk earthquake (M_w = 5.6) in the southwestern part of Sakhalin was preceded by a number of anomalous seismological and geophysical phenomena. The extensive data recorded by a network of digital seismic stations make it possible to track the aftershock dynamics of the process within 24 hours after the main event. The paper describes various manifestations of the earthquake.     

The Seismogenic Structure and Deformation Mechanism of the Lushan (Mw 6.6) Earthquake,Sichuan, China

On April 20 th, 2013, an earthquake of magnitude MW 6.6 occurred at Lushan of Sichuan on the southern segment of the Longmenshan fault zone, with no typical coseismic surface rupture. This work plotted an isoseismal map of the earthquake after repositioning over 400 post–earthquake macro–damage survey points from peak ground acceleration(PGA) data recorded by the Sichuan Digital Strong Earthquake Network. This map indicates that the Lushan earthquake has a damage intensity of IX on the Liedu scale, and that the meizoseismal area displays an oblate ellipsoid shape, with its longitudinal axis in the NE direction. No obvious directivity was detected. Furthermore, the repositioning results of 3323 early aftershocks, seismic reflection profiles and focal mechanism solutions suggests that the major seismogenic structure of the earthquake was the Dayi Fault, which partly defines the eastern Mengshan Mountain. This earthquake resulted from the thrusting of the Dayi Fault, and caused shortening of the southern segment of the Longmenshan in the NW–SE direction. Coseismal rupture was also produced in the deep of the Xinkaidian Fault. Based on the above seismogenic model and the presentation of coseismic surface deformation, it is speculated that there is a risk of more major earthquakes occurring in this region.     

Fault rupture and stress changes associated with the November 27, 2005 Qeshm Island (Iran) earthquake (<Emphasis Type="Italic">M</Emphasis><Subscript>W</Subscript> = 6.0)

The November 27, 2005 Qeshm Island earthquake (M_w 6.0) occurred along the Zagros Thrust and Fold Belt which accommodates about half of the deformation caused by the Arabian and Eurasian Plates convergence. As typical for the belt, the earthquake was associated with buried reverse faulting and produced no surface rupture. Here, teleseismic broadband P velocity waveforms of the earthquake are inverted to obtain coseismic finite-fault slip distribution of the earthquake. It is obtained that rupture was controlled by failure of a single asperity with largest displacement of approximately 0.6 m, which occurred at a depth of 9 km. The slip model indicated radial rupture propagation from the hypocentre and confirmed blind reverse faulting within deeper part (below the depth of 6 km) of the sedimentary cover above the Hormuz Salt, lying between the cover and the basement, releasing a seismic moment of about 1.3?×?10¹⁸ Nm (M_W?=?6.0). The results also confirm that the Hormuz Salt behaves as a barrier for rupture propagation to the basement below and occurrence of the aftershock activity downdip from the rupture within the Hormuz Salt. Calculated Coulomb stress variations caused by the coseismic rupture indicates stress coupling between the 2005 Qeshm Island earthquake and both the largest aftershock several hours later and the 2008 Qeshm Island earthquake (M_W?=?5.9). The stress calculations further indicated stress load at the depth range (15–20 km) of the well-located aftershocks, corresponding to depths of the Hormuz Salt and top of the basement and providing plausible explanation for occurrence of the aftershocks within those layers.     

Damage risk assessment of breakwaters under tsunami attack

Turkey was struck by two major events on 17 August and 12 November 1999, named Izmit (M _w = 7.4) and Düzce (M _w = 7.2) earthquakes, respectively. Rubble mound breakwaters in Izmit Bay experienced little damage, as forecasted by the new risk assessment model in which tsunami occurrence risk was included in the damage estimations. In order to determine the occurrence probability of structural damage under design conditions, including the environmental loading parameters of tsunami and storm waves, tidal range and storm surge, the Conditional Expections Monte Carlo simulation was applied in the risk assessment model developed in this study for the Esenköy Fishery Harbour, Turkey. A tsunami was not the key design parameter when compared to storm waves for the main breakwater of the harbour, however, in places with great seismic activity, the tsunami risk should be important depending on the occurrence probability and magnitude of the tsunami.