On the swelling potential of compacted high plasticity clays

We conducted a study of the spatial distributions of seismicity and earthquake hazard parameters for Turkey and the adjacent areas, applying the maximum likelihood method. The procedure allows for the use of either historical or instrumental data, or even a combination of the two. By using this method, we can estimate the earthquake hazard parameters, which include the maximum regional magnitude Mˆ_max, the activity rate of seismic events and the well-known bˆ value, which is the slope of the frequency-magnitude Gutenberg-Richter relationship. These three parameters are determined simultaneously using an iterative scheme. The uncertainty in the determination of the magnitudes was also taken into consideration. The return periods (RP) of earthquakes with a magnitude M ≥ m are also evaluated. The whole examined area is divided into 24 seismic regions based on their seismotectonic regime. The homogeneity of the magnitudes is an essential factor in such studies. In order to achieve homogeneity of the magnitudes, formulas that convert any magnitude to an M_S-surface scale are developed. New completeness cutoffs and their corresponding time intervals are also assessed for each of the 24 seismic regions. Each of the obtained parameters is distributed into its respective seismic region, allowing for an analysis of the localized seismicity parameters and a representation of their regional variation on a map. The earthquake hazard level is also calculated as a function of the form Θ = (Mˆ_max,RP_6.0), and a relative hazard scale (defined as the index K) is defined for each seismic region. The investigated regions are then classified into five groups using these parameters. This classification is useful for theoretical and practical reasons and provides a picture of quantitative seismicity. An attempt is then made to relate these values to the local tectonics.     

Seismic hazard in mega city Kolkata, India

The damages caused by recent earthquakes in India have been a wake up call for people to take proper mitigation measures, especially the major cities that lie in the high seismic hazard zones. Kolkata City, with thick sediment deposit (∼12 km), one of the earliest cities of India, is an area of great concern as it lies over the Bengal Basin and lies at the boundary of the seismic zones III and IV of the zonation map of India. Kolkata has been affected by the 1897 Shillong earthquake, the 1906 Calcutta earthquake, and the 1964 Calcutta earthquake. An analysis on the maximum magnitude and b-value for Kolkata City region is carried out after the preparation of earthquake catalog from various sources. Based on the tectonic set-up and seismicity of the region, five seismic zones are delineated, which can pose a threat to Kolkata in the event of an earthquake. They are broadly classified as Zone 1: Arakan-Yoma Zone (AYZ), Zone 2: Himalayan Zone (HZ), Zone 3: Shillong Plateau Zone (SPZ), Zone 4: Bay of Bengal Zone (BBZ), and Zone 5: Shield Zone (SZ). The maximum magnitude (m _max) for Zones 1, 2, 3, 4, and 5 are 8.30 ± 0.51, 9.09 ± 0.58, 9.20 ± 0.51, 6.62 ± 0.43 and 6.61 ± 0.43, respectively. A probability of 10% exceedance value in 50 years is used for each zone. The probabilities of occurrences of earthquakes of different magnitudes for return periods of 50 and 100 years are computed for the five seismic zones. The Peak Ground Acceleration (PGA) obtained for Kolkata City varies from 0.34 to 0.10 g.     

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.     

Seismic slip rates,potential subsurface rupture areas and seismic potential of the Vienna Basin Transfer Fault

The Vienna Basin Transfer Fault (VBTF) is a slow active fault with moderate seismicity (I _max~8–9, M _max~5.7) passing through the most vulnerable regions of Austria and Slovakia. We use different data to constrain the seismic potential of the VBTF including slip values computed from the seismic energy release during the 20th century, geological data on fault segmentation and a depth-extrapolated 3-D model of a generalized fault surface, which is used to define potential rupture zones. The seismic slip of the VBTF as a whole is in the range of 0.22–0.31 mm/year for a seismogenic fault thickness of 8 km. Seismic slip rates for individual segments vary from 0.00 to 0.77 mm/year. Comparing these data to geologically and GPS-derived slip velocities (>1 mm/year) proofs that the fault yields a significant seismic slip deficit. Segments of the fault with high seismic slip contrast from segments with no slip representing locked segments. Fault surfaces of segments within the seismogenic zone (4–14 km depth) vary from 55 to 400 km². Empirical scaling relations show that these segments are sufficiently large to explain both, earthquakes observed in the last centuries, and the 4th century Carnuntum earthquake, for which archeo-seismological data suggest a magnitude of M ≥ 6. Based on the combination of all data (incomplete earthquake catalog, seismic slip deficits, locked segments, potential rupture areas, indications of strong pre-catalog earthquakes) we argue, that the maximum credible earthquake for the VBTF is in the range M _max = 6.0–6.8, significantly larger than the magnitude of the strongest recorded events (M = 5.7).     

Site-specific earthquake response study for hazard assessment in Kolkata city,India

The assessment of local site effects on seismic ground motions is of great importance in earthquake engineering practice. Several destructive earthquakes in the past have demonstrated that the amplification of ground motion and associated damage to structures due to local site conditions is a significant consideration in earthquake hazard analysis. A recent paper published in this journal highlights the hazard posed by earthquakes in the megacity of Kolkata in India due to its seismic and geological settings. The seismic hazard assessment study speculates that the deep alluvial deposit in the city may increase the seismic hazard probably due to the amplification of the seismic energies. This paper focuses on the seismic response studies of the various soil strata (i.e. for local subsurface conditions) obtained from various construction sites in the city for predicted earthquake. It is very well recognized that site response studies (a part of seismic microhazard zonation for urban areas) are the first step towards performance-based foundation design or seismic risk analysis and mitigation strategy. One of the problems for carrying out site-specific study in Kolkata is the lack of recorded strong motion data in the city. Hence, this paper outlines a methodology to carry out site-specific study, where no strong motion data or seismic data are available. The methodology uses wavelet-based spectrum compatibility approach to generate synthetic earthquake motions and equivalent linear method for seismic site response analysis. The Mega City of Kolkata has been considered to explain the methodology. Seismic hazard zonation map by the Bureau of Indian Standards classifies the City of Kolkata as moderate seismic zone (Zone III) with a zone factor 0.16. On the other hand, GSHAP(Global Seismic Hazard Assessment Program) map which is based on 10% probability of exceedance in 50 years specifies a maximum peak ground acceleration (PGA) of 1.6 m/s² (0.163 g) for this region. In the present study, the seismic response has been carried out based on GSHAP. The results of the analysis indicate the amplification of ground motion in the range of 4.46–4.82 with the fundamental period ranging from 0.81 to 1.17 s. Furthermore, the maximum spectral accelerations vary in the range of 0.78–0.95 g.     

Probabilistic seismic hazard of Pakistan, Azad-Jammu and Kashmir

The seismic hazard study for Pakistan and Azad Jammu and Kashmir has been conducted by using probabilistic approach in terms of peak ground acceleration (PGA) in m/s² and also seismic hazard response spectra for different cities. A new version of Ambraseys et al. (Bull Earthq Eng 3:1–53, 2005) ground acceleration model is used, and parameterization is based on most recent updated earthquake catalogs that consisted of 14,000 events. The threshold magnitude was fixed at M _w 4.8, but seismic zones like northern Pakistan–Tajikistan, Hindukush and northern Afghanistan–Tajikistan border had M _w 5.2. The average normalized ‘a’ and ‘b’ values for all zones are 6.15 and 0.95, respectively. Seismicity of study area was modeled, and ground motion was computed for eight frequencies (0.025, 0.1, 0.2, 0.5, 1.0, 1.5, 2.0, 2.5 s) for different annual exceedance rates of 0.02, 0.01, 0.005, 0.002 and 0.001 (return periods 50, 100, 200, 500 and 1,000 years) for stiff rocks at the gridding of 0.1° × 0.1°. Seismic hazard maps based on computed PGA for 0.02, 0.01 and 0.002 annual exceedance are prepared. These maps indicate the earthquake hazard of Pakistan and surrounding areas in the form of acceleration contour lines, which are in agreement with geological and seismotectonic characteristics of the study area. The maximum seismic hazard values are found at Muzaffarabad, Gilgit and Quetta areas.     

Precursors and fluid flows in the case of the 1996, ML = 5.2 Saint-Paul-de-Fenouillet earthquake (Pyrenees,France): A complete pre-, co- and post-seismic scenario

Earthquake precursors are now regularly described but often detected only after a major or moderate seismic event. Presence and influence of fluids in the seismogenic processes are often observed at the time of earthquake studies. Even today, the understanding of the physical processes involved in the source region is a real challenge for seismic hazard assessment. Here, the aftershock sequence of the M_L = 5.2, 1996 Saint-Paul-de-Fenouillet (Eastern Pyrenees, France) earthquake is first re-examined with P-wave cross-correlations, resulting in extracting three multiplets and in determining new locations. Multiplets and spatio-temporal distribution analysis of the aftershocks allow for quantifying the hydraulic diffusivity D at a maximum value of 5 m²/s and the permeability K at 10^? 15 m² in the upper Pyrenean crust. Second, a model is established in order to explain the hydrogeochemical transient anomalies, which occurred during the 15 day-period preceding the 1996 earthquake. These anomalies consist on a temporal and spatial sequence of gas emissions in the epicentral area and on chloride and lead concentration variations in a bottled mineral water 25 km north to the main shock epicenter. The proposed model processed in a standard elastic half-space, consists of creep on a low-angle crustal normal-fault, generating volumetric strain field changes over a distance of 25 km from the epicentral area. This model is able to constrain not only the mechanisms and the locations of the geochemical anomalies, but also their timing and probable casual links to the triggering of the impending major event. Also, the active extension proposed here is compatible with seismological observations in the Pyrenees. Thus, the possibility of such creep, which can be considered as a slow-slip event, is discussed in the Pyrenean tectonic and geological context. The model is discussed and compared to previous proposed models on precursor processes of earthquakes, especially concerning the preparation zone concept. Finally, a complete seismic scenario over the period beginning 15 days before the quake and ending 5 days after is proposed and discussed.     

Analysis of seismic hazard in landslide-prone regions: criteria and example for an area of Daunia (southern Italy)

In relation to the assessment of earthquake-induced landslide hazard, this paper discusses general principles and describes implementation criteria for seismic hazard estimates in landslide-prone regions. These criteria were worked out during the preparation of a hazard map belonging to the official Italian geological cartography and they are proposed as guidelines for future compilation of similar maps. In the presented case study, we used a procedure for the assessment of seismic hazard impact on slope stability adopting Arias intensity Ia as seismic shaking parameter and critical acceleration a _c as parameter representing slope strength to failures induced by seismic shaking. According to this procedure, after a preliminary comparison of estimated historical maximum values of Ia with values proposed in literature as landslide-triggering thresholds, a probabilistic approach, based on the Newmark’s model, is adopted: it allows to estimate the minimum critical acceleration a _c required for a slope to keep under a prefixed value, the probability of failures induced by seismic shakings expected in a given time interval. In this way, one can prepare seismic hazard maps where seismic shaking is expressed in an indirect way through a parameter (the critical acceleration) representing the “strength” that seismic shakings mobilise in slope materials (strength demand) with a prefixed exceedance probability. This approach was applied to an area of Daunia (Apulia—southern Italy) affected by frequent landslide phenomena. The obtained results indicate that shakings with a significant slope destabilisation potential can be expected particularly in the north-western part of the area, which is exposed to the seismic activity of Apennine tectonic structures.     

Seismic Hazard and Risk Assessment for Tulbagh,South Africa: Part I – Assessment of Seismic Hazard

In this part of our study the probabilistic seismic hazard analysis (PSHA) for Tulbagh was performed. The applied procedure is parametric and consists essentially of two steps. The first step is applicable to the area in the vicinity of Tulbagh and requires an estimation of the area-specific parameters, which, in this case, is the mean seismic activity rate, , the Gutenberg-Richter parameter, b, and the maximum regional magnitude, m_max. The second step is applicable to the Tulbagh site, and consists of parameters of distribution of amplitude of the selected ground motion parameter. The current application of the procedure provides an assessment of the PSHA in terms of peak ground acceleration (PGA) and spectral acceleration (SA). The procedure permits the combination of both historical and instrumental data. The historical part of the catalogue only contains the strongest events, whereas the complete part can be divided into several subcatalogues, each assumed complete above a specified threshold of magnitude. In the analysis, the uncertainty in the determination of the earthquake was taken into account by incorporation of the concept of `apparent magnitude'. The PSHA technique has been developed specifically for the estimation of seismic hazard at individual sites without the subjective judgement involved in the definition of seismic source zones, when the specific active faults have not been mapped or identified, and where the causes of seismicity are not well understood. The results of the hazard assessment are expressed as probabilities that specified values of PGA will be exceeded during the chosen time intervals, and similarly for the spectral accelerations. A worst case scenario sketches the possibility of a maximum PGA of 0.30g. The results of the hazard assessment can be used as input to a seismic risk assessment.     

汶川震区成兰铁路关键段多尺度滑坡危险性评估

以穿越汶川震区的成兰铁路龙门山关键段为例, 探索提出了强震扰动背景下重大工程场区多尺度滑坡危险性评估方法。利用信息量模型反演评估了汶川地震诱发的同震滑坡空间分布特征, 以此为前提开展了区域和局地两种空间尺度的滑坡危险性预测评估。在区域廊带尺度上, 分别利用可能最大降雨量预测方法和信息量模型, 进行了日超越概率10%的最大降雨量时空分布预测及其诱发滑坡的危险性评估; 同时, 结合地震危险性区划成果, 开展了50年超越概率10%的基本地震动诱发滑坡的危险性评估。在局地场站尺度上, 利用基于崩塌运动过程模拟的Rockfall Analyst软件, 开展了柿子园大桥周边崩塌运动学特征(Runout)模拟和危险性评估。滑坡和崩塌危险性评估的结果分别为铁路规划选线和场站防护设计提供了不同尺度的地质安全依据。      

Probabilistic seismic hazard assessment in the Constantine region,Northeast of Algeria

This paper presents a seismic hazard evaluation and develops an earthquake catalogue for the Constantine region over the period from 1357 to 2014. The study contributes to the improvement of seismic risk management by evaluating the seismic hazards in Northeast Algeria. A regional seismicity analysis was conducted based on reliable earthquake data obtained from various agencies (CRAAG, IGN, USGS and ISC). All magnitudes (M _l, m _b) and intensities (I ₀, I _MM, I _MSK and I _EMS) were converted to M _s magnitudes using the appropriate relationships. Earthquake hazard maps were created for the Constantine region. These maps were estimated in terms of spectral acceleration (SA) at periods of 0.1, 0.2, 0.5, 0.7, 0.9, 1.0, 1.5 and 2.0 s. Five seismogenic zones are proposed. This new method differs from the conventional method because it incorporates earthquake magnitude uncertainty and mixed datasets containing large historical events and recent data. The method can be used to estimate the b value of the Gutenberg-Richter relationship, annual activity rate λ(M) of an event and maximum possible magnitude M _max using incomplete and heterogeneous data files. In addition, an earthquake is considered a Poisson with an annual activity rate λ and with a doubly truncated exponential earthquake magnitude distribution. Map of seismic hazard and an earthquake catalogue, graphs and maps were created using geographic information systems (GIS), the Z-map code version 6 and Crisis software 2012.     

Seismic hazard assessment in Jordan and its vicinity

The earthquake hazard in Jordan and its vicinity is assessed on the basis of probabilistic methods. For this purpose, an updated earthquake catalog is compiled which covers the period between AD 1–1989. The earthquakes lie between latitudes 27.0°-35.5° N and longitudes 32.0°-39.0° E. Thirteen seismic zones are defined on a regional seismic and tectonic map presented for the area. Point-source and line-source models are used. The seismic hazard parameters, namely, theb-parameter (of the Gutenberg-Richter relation),m ₁ (the upper bound magnitude), and ₄ (the annual rate of occurrence of earthquakes with local magnitudeM _L 4.0) are calculated for each zone. The results of the seismic hazard assessment are displayed as iso-acceleration contours expected to be exceeded during typical economic life times of structures, i.e. 50 and 100 years. For each model, two seismic hazard maps are derived. In order to determine the importance of the South-eastern Mediterranean zone and the north part of the Red Sea zone from a seismic hazard point of view for Jordan, one seismic hazard map which corresponds to 50 years' economic life for every model, excluding the seismicity of these zones, is derived.     

The July 2007 rock and ice avalanches at Mount Steele, St. Elias Mountains, Yukon, Canada

A large rock and ice avalanche occurred on the north face of Mount Steele, southwest Yukon Territory, Canada, on July 24, 2007. In the days and weeks preceding the landslide, several smaller avalanches initiated from the same slope. The ice and rock debris traveled a maximum horizontal distance 5.76 km with a maximum vertical descent of 2,160 m, leaving a deposit 3.66 km² in area on Steele Glacier. The seismic magnitude estimated from long-period surface waves (M _s) is 5.2. Modeling of the waveforms suggests an estimated duration of approximately 100 s and an average velocity of between 35 and 65 m/s. This landslide is one of 18 large rock avalanches known to have occurred since 1899 on slopes adjacent to glaciers in western Canada. We describe the setting, reconstruct the event chronology and present a preliminary characterization of the Mount Steele ice and rock avalanches based on field reconnaissance, analysis of seismic records and an airborne LiDAR survey. We also present the results of a successful dynamic simulation for the July 24 event.     

Estimation of earthquake hazard parameters for incomplete and uncertain data files

The maximum likelihood estimation of earthquake hazard parameters (maximum regional magnitudem _max, activity rate , and theb parameter in the Gutenberg-Richter distribution) is extended to the cases of incomplete and uncertain data. The method accepts mixed data containing only large (extreme) events and a variable quality of complete data with different threshold magnitude values. Uncertainty of earthquake magnitude is specified by two values, the lower and upper magnitude limits. It is assumed that such an interval contains the real unknown magnitude. The proposed approach allows the combination of different quality catalog parts, e.g. those where the assignment of magnitude is questionable and those with magnitudes precisely determined.As an illustration of the method, the seismic hazard analysis for western Norway and adjacent sea area (4–8°E, 58–64°N) is presented on the basis of the strongest earthquakes felt during the period 1831–1889 and three complete catalog parts, covering the period 1890–1987.Paper presented at the 21st General Assembly of the European Seismological Commission held in Sofia, 1988.On leave from Institute of Geophysics, Polish Academy of Sciences, Warsaw, Poland.     

Evidences of landslide earthquake triggering due to self-excitation process

The basin-like setting of stiff bedrock combined with pre-existing landslide masses can contribute to seismic amplifications in a wide frequency range (0–10 Hz) and induce a self-excitation process responsible for earthquake-triggered landsliding. Here, the self-excitation process is proposed to justify the far-field seismic trigger of the Cerda landslide (Sicily, Italy) which was reactivated by the 6th September 2002 Palermo earthquake (M _s = 5.4), about 50 km far from the epicentre. The landslide caused damage to farm houses, roads and aqueducts, close to the village of Cerda, and involved about 40 × 10⁶ m³ of clay shales; the first ground cracks due to the landslide movement formed about 30 min after the main shock. A stress–strain dynamic numerical modelling, performed by FDM code FLAC 5.0, supports the notion that the combination of local geological setting and earthquake frequency content played a fundamental role in the landslide reactivation. Since accelerometric records of the triggering event are not available, dynamic equivalent inputs have been used for the numerical modelling. These inputs can be regarded as representative for the local ground shaking, having a PGA value up to 0.2 m/s², which is the maximum expected in 475 years, according to the Italian seismic hazard maps. A 2D numerical modelling of the seismic wave propagation in the Cerda landslide area was also performed; it pointed out amplification effects due to both the structural setting of the stiff bedrock (at about 1 Hz) and the pre-existing landslide mass (in the range 3–6 Hz). The frequency peaks of the resulting amplification functions (A(f)) fit well the H/V spectral ratios from ambient noise and the H/H spectral ratios to a reference station from earthquake records, obtained by in situ velocimetric measurements. Moreover, the Fourier spectra of earthquake accelerometric records, whose source and magnitude are consistent with the triggering event, show a main peak at about 1 Hz. This frequency value well fits the one amplified by the geological setting of the bedrock in correspondence with the landslide area, which is constituted of marly limestones and characterised by a basin-like geometry.     

Neotectonics at the Arabian plate margins

Opening of the Red Sea is accompanied by convergence between the Arabian plate and Eurasia. Regional topography and structure favour gravity glide as the main driving force of plate translation. At the leading edge of the plate, the Zagros Mountains undergo coseismic serial folding which is equivalent to Holocene shortening by ∼20 mm/year and which has led to major episodes of coastal uplift of which the last was ∼1700 years BP. At the Jordan Rift transform, which bounds the Arabian plate on the west, a recurrence interval of ∼1600 years is reported for events of M_L≥5.5. The palaeomagnetic record for the last 3.2 Ma indicates an average spreading rate for the Red Sea of ∼20 mm/year; there is some evidence that hydrothermal activity in the Red Sea is pulsatory, with a period of ∼2000 year, and that it reflects discontinuous spreading. The Holocene neotectonic records of the Zagros, the Jordan Rift and the Red Sea are the product of complex plate interactions and of the accumulation and release of strain in the crust along the plate margins. But they also reflect elastic strain energy storage and release within the Arabian plate, whence parallels in the period of major deformation episodes in the three deforming zones and the apparent discrepancy between the seismic moment predicted by plate kinematics and that recorded in the Zagros. Any associated intraplate deformation, if detected geodetically, would thus help the assessment of seismic hazard.     

First Order Seismic Microzonation of Delhi,India Using Geographic Information System (GIS)

A first order seismic microzonation map of Delhi is prepared using five thematic layers viz., Peak Ground Acceleration (PGA) contour, different soil types at 6 m depth, geology, groundwater fluctuation and bedrock depth, integrated on GIS platform. The integration is performed following a pair-wise comparison of Analytical Hierarchy Process (AHP), wherein each thematic map is assigned weight in the 5-1 scale: depending on its contribution towards the seismic hazard. Following the AHP, the weightage assigned to each theme are: PGA (0.333), soil (0.266), geology (0.20), groundwater (0.133) and bedrock depth (0.066). The thematic vector layers are overlaid and integrated using GIS. On the microzonation theme, the Delhi region has been classified into four broad zones of vulnerability to the seismic hazard. They are very high (> 52%), high (38–52%), moderate (23–38%) and less ( < 23%) zones of seismic hazard. The “very high” seismic hazard zone is observed where the maximum PGA varies from 140 to 210 gal for a finite source model of M_w 8.5 in the central seismic gap. A site amplification study from local and regional earthquakes for Delhi region using Delhi Telemetry Network data shows a steeper site response gradient in the eastern side of the Yamuna fluvial deposits at 1.5 Hz. The ‘high’ seismic hazard zone occupies most of the study area where the PGA value ranges from 90 to 140 gal. The ‘moderate’ seismic hazard zone occurs on either side of the Delhi ridge with PGA value varying from 60 to 90 gal. The ‘less’ seismic hazard zone occurs in small patches distributed along the study area with the PGA value less than 60 gal. Site response studies, PGA distribution and destruction pattern of the Chamoli earthquake greatly corroborate the seismic hazard zones estimated through microzonation on GIS platform and also establishes the methodology incorporated in this study.     

Recent evaluation of the assessment seismic hazards for Nuweiba, Gulf of Aqaba

The Gulf of Aqaba is considered seismically as one of the most active zones of the Dead Sea Transform region. The main shock of the 1995 Gulf of Aqaba earthquake sequence is considered as the largest shock in the (surface wave magnitude Ms?=?7.2) since the sixteenth century. The present study is a trial to detect the probabilistic seismic hazard analysis (PSHA) for Nuweiba site. Data used for this study was a combination of both historical and recent instrumental data. Results of the hazard assessment, expressed as in the worst case scenario, reveal that Nuweiba is exposed to the occurrence of a maximum credible earthquake of magnitude $ m_{{\max }} ~ = ~7.4 \pm 0.31 $ , at hypocentral distance of 15.6?±?10 km. For structure with the return period of 100 years, with a 90% probability of exceedance, the maximum expected earthquake magnitude (M_L) is 5.9 in this lifetime. The possibility of the maximum peak ground acceleration at the Nuweiba site is 0.41 g. Results of the hazard assessment can be used as an input data to assess the seismic risk for site of interest.