Ground Motion Prediction Equations Research Articles

Physically Based and Empirical Ground Motion Prediction Equations for Multiple Intensity Measures (PGA, PGV, Ia, FIV3, CII, and Maximum Fourier Acceleration Spectra) on Sakhalin Island

In the present study, empirical attenuation relations for multiple ground motion intensity measures (PGA, PGV, Ia, FIV3, CII, and MFAS) were developed for Sakhalin Island (in the far east of Russia). A recorded strong motion dataset was used, making GMPEs applicable in active crustal regions with an earthquake magnitude range of 4–6 and a distance range of up to 150 km. The hypocentral distance was used as a basic distance metric. For the first time in the research, an analytical representation of Arias intensity (Ia) was obtained in the framework of a multi-asperity source model. Asperities are considered as sub-sources of high-frequency incoherent radiation. The physical representation of the attenuation model in our study was based on a stress drop on the asperities and the ratio of the total rupture area to the combined area of asperities. The average stress drop on asperities for the examined earthquakes was approximately 13.4 MPa, and the ratio of the total rupture area to the asperity area was 0.22, which is generally close to similar estimates for crustal earthquakes. The coefficients and statistical scattering of the attenuation models were also analyzed. Moreover, a magnitude scale based on a modified Arias intensity is proposed in the present study. The new magnitude scale has an explicit physical meaning and is characterized by its simplicity of measurement. It is associated with the acceleration source spectrum level and can be successfully used in early warning systems.

Ground-Motion Model for Hard-Rock Sites by Correction of Surface Recordings (Part 2): Correction, Mixed-Effects Regressions, and Results

ABSTRACT In the framework of site-specific seismic hazard assessment, the definition of reference motion is a crucial step. Reference motion is generally associated with hard-rock conditions, characterized by S-wave velocity exceeding 1500 m/s. However, ground motion recorded at sites with such conditions is poorly represented in existing strong-motion databases. Thus, the validity domains of most empirical ground-motion prediction equations (GMPEs) are not representative of reference rock conditions. To overcome this limitation and assess ground motion at reference conditions, the so-called “deconvolution approach” was proposed by Laurendeau et al. (2018) to correct surface recordings from theoretical 1DSH site response before GMPE developments. With the same purpose, in this article, we propose to apply the deconvolution approach using empirical site-response estimates as an alternative to theoretical ones. Using the Kiban–Kyoshin network (KiK-net) data, we estimate empirical site responses at KiK-net stations using generalized inversion techniques in addition to those from 1DSH numerical simulations, as presented in the companion article. Finally, a reference ground-motion model (RGMM) is determined based on empirically deconvolved ground motions. The advantage of using empirical rather than 1DSH site responses in the deconvolution approach is that in the former case the RGMM can be built based on records from an extensive set of sites, whereas the latter case is restricted to well-characterized sites with dominant 1D behavior. This makes the proposed approach easily exportable to different regions of the world, where precise site characterizations are not systematically available, and the knowledge of site behavior is limited.

Probabilistic seismic hazard and risk assessment of Mizoram, North East India

This paper presents seismic hazard and risk assessment for the state of Mizoram based on a classical probabilistic seismic hazard analysis and event-based probabilistic seismic risk analysis. For the seismic hazard estimation, analysis has been performed considering the areal source model, fault zone polygon, and smoothed gridded seismicity model. The earthquake activity rates for these source zones and smoothed gridded seismicity sources are estimated from the homogenized and declustered earthquake catalogue. The logic tree framework is applied in the seismic source models and Ground Motion Prediction equations (GMPEs) to account for the epistemic uncertainties. Five Next-generation attenuation (NGA) GMPEs for the active shallow region and three GMPEs for the Indo-Burma subduction zone have been used to evaluate the hazard at the reference rock condition (Vs30 ​= ​760 ​m/s). Peak Ground Acceleration (PGA) and Spectral Acceleration (SA) at 0.2 s and 1s are estimated for each eight districts headquarters of Mizoram for a 2% and 10% probability of exceedance in 50 years. The hazard curves and Uniform Hazard Spectra (UHS) are also presented. For seismic risk analysis, building exposure data are based on digitized building footprint and 2011 Housing Census data of India. All the buildings are classified into three classes, and seismic vulnerability functions are assigned to each building class. The area per building class is assigned from the digitized footprint. Building replacement costs per square meter have been chosen based on expert input and values identified from past study. Lastly, the study conducted a seismic risk analysis using the Open Quake-engine's probabilistic event-based methodology to estimate risk metrics at the district level, such as average annual losses and probability curves for loss exceedance. The study's findings provide valuable insights into the most high-risk areas, the building construction types that are most vulnerable to seismic activity, and the anticipated economic losses in the state of Mizoram. These results can serve as a guide for local government authorities in developing future city plans and implementing earthquake risk mitigation strategies.

A Logic-Tree Approach for Probabilistic Seismic Hazard Assessment in the Administrative Region of Attica (Greece)

Probabilistic Seismic Hazard Assessment (PSHA) was carried out for the administrative region of Attica (Greece). Peak Ground Acceleration (PGA) and Peak Ground Velocity (PGV) values were calculated for return periods of 475 and 950 years for five sub-areas covering the entire region. PGA hazard curves and Uniform Hazard Spectra (UHS) in terms of spectral acceleration (Sa) values were generated for Athens, Methana, and the capitals of each island of Attica (Salamina, Aegina, Poros, Hydra, Spetses, Kythira, and Antikythira). Area sources were adopted from the Euro-Mediterranean Seismic Hazard Model 2013 (ESHM13) and its update, ESHM20, taking into account both crustal and slab tectonic environments. Ground Motion Prediction Equations (GMPEs) proposed for the Greek territory were ranked for PGA and PGV. Each GMPE was reconstructed as a weighted model, accounting for normal and non-normal focal mechanisms for each area source. PGA, PGV, and Sa values were computed using a logic tree, integrating the seismotectonic models as major branches and sub-logic trees, comprised of multiple ranked GMPEs for each area source, as minor branches. The results showed higher seismic hazard values in sub-areas near the Gulf of Corinth and the slab interface, which could indicate a need to revise the active building code in Attica.

Insights on the role of local site effects on damage distribution in the Izmir metropolitan area induced by the October 30, 2020 Samos earthquake

On the 30th of October 2020, a 6.6 magnitude earthquake occurred 14 km north of Samos Island, causing 119 casualties (117 in Izmir, Türkiye, and 2 in Samos, Greece) and significant damage in the 3rd biggest city of Türkiye, Izmir. Although the city is roughly 70 km far away from the epicenter, the damage was significant and concentrated in the city center settled on alluviums. This paper aims to analyze the distribution of damage in Izmir province, by crosschecking the recorded motions, the subsoil conditions and the evidence of damage as collected by an ad-hoc on-site reconnaissance. The intrinsic behavior of the Samos earthquake was investigated by employing three different ground-motion prediction equations. The results of the analyses revealed that site effects play a significant role in the amplification of ground motions, and valley effects are responsible for the concentration of damage. The damage in buildings was classified in terms of the intensity and structural typologies for the 30 districts of Izmir metropolitan area. In-depth analysis of the distribution of damages revealed that the earthquake caused damage all over the boundaries of Izmir province, and the concentration of damage in Bornova and Karşıyaka districts has a clear correlation with double resonance effects.

Earthquake Damage Assessment of Buried Pipeline Networks in the Lima Metropolitan Area

Lifelines such as drinking water and sewage systems provide the means and conveyance for daily critical services, and they are essential systems for recovery operations after a damaging earthquake. Therefore, earthquake damage to lifeline components needs to be reliably assessed for a possible future seismic scenario. This study presents an earthquake damage assessment of buried pipeline networks in the Lima Metropolitan Area (LMA). It includes seismic hazard analysis, a review of pipeline network datasets, and the selection of empirical fragility functions. Deterministic seismic hazard analysis was performed for an inter-plate earthquake scenario using ground motion prediction equations and site conditions to compute the distribution of the peak ground velocity (PGV). Recommendations are offered for an adequate selection of fragility functions developed in other regions, and a logic tree of fragility functions is proposed to be used in pipelines of LMA according to the data of pipeline damage after the 2007 Pisco Earthquake. Finally, the pipeline repair ratios and the total number of repairs are estimated for the earthquake scenario, and the results are geographically presented for each pipeline network.

Characterization of strong‐motion recording sites in the Himalayas and adjoining areas

AbstractIn this paper, an attempt has been made to arrive at sufficiently realistic estimates of the average shear‐wave velocity in the top 30 m of the ground () for 200 strong motion recording sites spread over the Himalayas, adjoining Indo‐Gangetic plains and northeast Indian regions. This is achieved by extending the application of the horizontal to vertical spectral ratio (HVSR) method for the strong motion recording sites with a limited number of records. The knowledge of surface geology at the recording sites and the estimates based on the topographic slope method has been utilized as supplementary information to arrive at a good first‐order approximation of for the sites with a single record also. A total of 669 three‐component acceleration records have been obtained at the 200 sites considered from 158 different earthquakes in the magnitude range of 2.3–7.9 during 1986–2015. It had not been possible to use this database fruitfully due to lack of the required site characterization. Compared to the costly and time‐taking field studies required for the estimation of exact site‐specific values, the proposed scheme can be considered to provide a convenient methodology for obtaining sufficiently reliable approximations of from the recorded strong motion data. The values obtained in this study are found to be in very good agreement with the site‐specific experimental values available for 29 of the sites, establishing the validity of the proposed scheme. An alternative site characterization is also provided in terms of the site soil and the local geological conditions, both defined in terms of three categories specified by indicator variables with values of 0, 1 and 2. Lee and Trifunac have illustrated that this characterization has stronger correlation with the ground motion amplitudes than the as well as the NEHRP site classes based on broad intervals of . The site characterization of the 200 strong motion recording sites in the present study will find useful applications in developing region‐specific empirical ground motion prediction equations with the site amplification effects modeled realistically and in selecting suitable equations from the published literature using score‐based data‐driven methods.

Ground motion hazard of the China–Pakistan Economic Corridor (CPEC) routes in Pakistan

Pakistan has seen a burst of infrastructure development recently due to the increased connection between Asia and East Europe. The China–Pakistan Economic Corridor is a project between China and Pakistan aimed to improve the regional infrastructure that would ultimately enhance the connection between Asia and Eastern Europe. However, the active tectonics of Pakistan could put this infrastructure at risk if it is not built to the highest hazard prevention standard. This study reports the ground motion hazard by using the probabilistic seismic hazard assessment approach and the areal seismic source model. The seismic hazard maps of the China–Pakistan Economic Corridor in Pakistan are derived using the Cornell–McGuire (1968–1976) approach, which takes into account all earthquakes (25AD-2020) that occurred in Pakistan and nearby regions, the newest ground motion prediction equations, and an updated seismotectonic source model of Pakistan. The final ground motion intensities are attained as peak ground acceleration and 5% damped spectral acceleration at T = 0.2 s and 1.0 s for 475- and 2475-year return periods (estimated for bedrock site conditions). The results are displayed as color-coded maps that represent the amplitude deviation of ground motion. From the spatial evaluation of the maps, a peak ground acceleration value of 0.40–0.52 g for the 475-year return period and a spectral acceleration (0.2 s) value of 1.66–2.13 g for 2475-year return period are mostly observed on the northern and western routes. The central and eastern routes are mostly characterized by a peak ground acceleration value of 0.22–0.24 g for the 475-year return period and a spectral acceleration (0.2 s) value of 0.95–1.13 g due to diffused seismicity and lower number of faults in this region. The ground motion intensity values obtained in this study can be utilized for the seismic design of all kinds of infrastructure and bridges along the CPEC routes in accordance with the Building Code of Pakistan, the International Building codes, and the load and resistance factor design codes published by American Association of the State Highway and Transportation Officials.

Updated probabilistic seismic hazard assessment of Pakistan

A probabilistic seismic hazard analysis of Pakistan is carried out to produce seismic hazard maps showing peak ground acceleration for the 2% and 10% of probability of exceedance in 50 years and spectral acceleration for the 5% and 2% probabilities of exceedance in 50 years. The main objective of this work is to define seismic hazard definition for Pakistan based on the classical approach using updated seismic hazards analysis inputs. A composite earthquake catalogue consisting of 32,700 events was compiled and used in the analysis with a magnitude range of Mw 4.0–8.3. The computations were made at a rectangular grid of 5 km in the OpenQuake engine. Ground motion parameter values have been obtained for flat rock reference seismic site conditions having a shear wave velocity of 760 m/s. Through the logic tree, epistemic uncertainties inherent in ground motion prediction equations and the maximum magnitude potential of seismic sources were estimated. Ground motion prediction equations were assigned equal weights in the logic tree while different weights were assigned to the maximum magnitude potential models. Results of the study were expressed as ground motion contour maps and median uniform hazard spectra for important cities in Pakistan. The results showed that peak ground acceleration values in Pakistan range from 0.16 to 0.54 g for the 10% of probability of exceedance, 0.23–0.72 g for the 5% probability of exceedance and 0.32–1.02 g for the 2% of probability of exceedance in 50 years. Spectral acceleration values at 0.2 s range from 0.67 to 2.19 g for the 2% chance of exceedance in 50 years, respectively, whereas spectral acceleration at 1.0 s ranges from 0.09 to 0.52 g for the 2% chance of exceedance in 50 years. The centre of Pakistan and the southern and the northern region are exposed to high hazard compared to the central-eastern part of the country. Quetta, Chitral, Peshawar, Islamabad Ziarat and Gwadar are exposed to seismic hazard in the PGA range of 0.32–0.45 g. Comparison of the results of this study with those of similar studies suggests that the seismic estimates proposed are reliable. The study provides a demonstration that the classical seismic hazard approach applied here is rational and applicable in Pakistan.

Joint Occurrence of Spectral Accelerations Between Horizontal and Vertical Ground Motions

ABSTRACT This paper estimates the joint occurrence of spectral accelerations (SA) between horizontal (HGM) and vertical ground motions (VGM) based on the NGA-West2 ground motion database, which includes the correlation model and joint distribution model. In this study, the correlation model calculated via CB14 horizontal ground motion prediction equation (GMPE) and BC16 vertical GMPE is obtained. It is verified that the joint distribution of residuals follows a multivariate normal distribution via both a set of quantitative validation methods and a qualitative validation method. Conditional mean spectra of VGMs conditioned on SAs of HGMs are studied as an application of joint occurrence of SA V–SA H.

A Rapid Estimation Method for Post-earthquake Building Losses

Rapid estimation of post-earthquake building damage and loss is very important in urgent response efforts. The current approach leaves much room for improvement in estimating ground motion and correctly incorporating the uncertainty and spatial correlation of the loss. This study proposed a new approach for rapidly estimating post-earthquake building loss with reasonable accuracy. The proposed method interpolates ground motion based on the observed ground motion using the Ground Motion Prediction Equation (GMPE) as the weight. It samples the building seismic loss quantile considering the spatial loss correlation that is expressed by Gaussian copula, and kriging is applied to reduce the dimension of direct sampling for estimation speed. The proposed approach was validated using three historical earthquake events in Japan with actual loss reports, and was then applied to predict the building loss amount for the March 2022 Fukushima Mw7.3 earthquake. The proposed method has high potential in future emergency efforts such as search, rescue, and evacuation planning.

Implications of GMPE’s structure for multi-site seismic hazard

Previous research has discussed the implications of the structure of a classical ground motion prediction equation (GMPE) on single-site probabilistic seismic hazard analysis and disaggregation. Classical refers to a GMPE where local site conditions (or any other factor) are accounted for via constant (with respect to magnitude and source-to-site distance) terms added to the mean and that do not affect the distribution of the residuals. Herein, the implications of such a structure of the GMPE are briefly discussed with respect to multi-site hazard assessment that, typically, requires a large number of simulations of random fields of ground motion intensity measures. It is shown that this type of GMPEs enables to run the simulations only once, independently of the soil conditions (or any other factor modeled in a similar way) eventually assigned to each site, which can represent a significant computational advantage in the case of spatially distributed assets.

Comparison of empirically-based and physically-based analyses of coseismic landslides: A case study of the 2016 Kumamoto earthquake

This study compares the performance of empirically- and physically-based methods for predicting coseismic landslides around Aso caldera during the 2016 Mw 7.0 Kumamoto earthquake. The physically-based method couples a regional-scale wave propagation simulation with a site-scale Newmark-type sliding analysis, while the empirically-based method predicts sliding displacements based on seismic recordings and/or ground motion prediction equations. By using the well-documented landslide inventory, the predictive capacity of each method is quantitatively assessed.The case study demonstrates that the physically-based method has the advantage of simulating complicated near-fault and topographic effects. It captures around 44% of observed landslides at the caldera rim, which is significantly better than the 20% captured from the empirically-based prediction. Both methods have similar prediction performance, capturing around 56% actual landslides at the central cone region. The empirically-based method is simple to use and provides fast landslide prediction. It can be further improved by incorporating a simple prediction model for the topographic amplification effect, which increases the percentage of captured landslides by around 8% across the caldera rim.

Next generation ground-motion prediction equations for Indo-Gangetic Plains, India

Seismic hazard analysis (SHA) of Indo-Gangetic Plains (IGP), India, and its proximity to the Himalayas requires reliable ground motion prediction equations (GMPEs). This study attempts to derive the GMPEs for IGP using strong-motion accelerometer data recorded from 2005 to 2015 on IGP. For regression, peak ground acceleration (PGA) and pseudo-spectral acceleration (PSA) of 5% damped linear pseudo-absolute acceleration response spectra at 27 periods ranging from 0.01–10 s were used. Two-stage nonlinear regression trains the functional form of nonlinear magnitude scaling, distance scaling, and site conditions. The model includes a regionally independent geometric attenuation finite fault distance metric, style of faulting, shallow site response, basin response, hanging wall effect, hypocentre depth, regionally dependent anelastic attenuation, site conditions, and magnitude-dependent aleatory variability. Developed GMPEs are validated for active crustal continental earthquakes for epicentral distances REPI ranging from 10–1500 km, magnitude ranging from 3.3–7.9, and focal depths 1–70 km. The GMPEs developed are compared with the Campbell and Bozorgnia 2008, 2013 and 2014, and North Indian GMPEs, which are agreed upon consistently. The model can predict the horizontal ground motion for SHA on IGP, considering possible maximum magnitude earthquakes from the seismic gaps of the Himalayas.

Study of Ground Motion Models Selection to the Concrete Frame Building Damage Probability using the 2006 Yogyakarta Earthquake Scenario

Damage probability assessment using the ground motion prediction equation (GMPE) models is paramount in seismic risk analysis. Several GMPE models have been proposed and developed worldwide. However, it has its parameters and gives different spectral acceleration results. The study aimed to investigate the impact of GMPE model selection on the estimation of building damage probability due to the scenario of Yogyakarta earthquake. In the beginning, we modeled the 3D-concrete moment frame building, followed by conducting pushover analysis to determine the building responses and calculating the building damage probability based on HAZUS method to generate fragility curve and damage probability matrix. Five GMPE models were employed in this study. The results revealed that for the modeled structure, the GMPE of Ambraseys et al. (1996) gave the highest value of the peak building response compared to the GMPE of Sadigh et al. (1997), Boore et al. (1997), Spudich et al. (1999) and Idriss (2002). The probability of moderate, extensive, and complete damage levels was also dominated by the GMPE of Ambraseys et al. (1996), which indicated that in this case, the use of this GMPE would produce the most severe damage probability value.

Spatial correlations in ground motion intensity measuring from the 2023 Türkiye earthquake

When evaluating an area's seismic risk or resilience, it is necessary to use the spatial correlation to analyze the ground motion parameters of multiple sites together in an earthquake. These two large earthquakes in Türkiye provided the possibility for spatial correlation analysis of ground motion intensity measurements in this area. Based on the strong motion records provided by The Disaster and Emergency Management Authority of Türkiye (AFAD), this study uses the local ground motion prediction equation ​in Türkiye to give spatial correlation analysis of Intensity Measurements. This study gives an exponential model based on a semivariogram and compares it with the correlation model obtained from previous studies.

New ground-motion prediction equations for significant duration of shallow crustal and upper mantle earthquakes in Japan and intra-event spatial correlation

Several models have been proposed and developed to formulate prediction equations for the significant duration (DS) of shallow crustal and subduction earthquakes. However, no equation for DS and a corresponding spatial correlation model have been proposed for upper mantle earthquakes. To address records are selected from the Kyoshin network (K-net) and Kiban Kyoshin network (KiK-net) databases in this study for shallow crustal and upper mantle earthquake with moment magnitudes of 4 ≤ Mw ≤ 7.3, rupture distances Rrup ≤ 300 km and peak ground accelerations (PGAs) ≥ 0.01 g. A new prediction equation is established for DS of shallow crustal and upper mantle earthquakes in Japan by considering the influences of Mw, Rrup, the time-averaged shear-wave velocity to 30 m VS30, the difference between the measured and predicted values of the overburden thickness for a shear wave velocity of 1000 m/s δZ1, and the depth to the top of the rupture Ztor on the DS of ground motion. The rationality and reliability of the proposed prediction model is confirmed using residual analysis and by comparison with existing prediction models. Based on intra-event residuals and single-station intra-event residuals, a semivariogram function and exponential model are used to analyze the spatial correlation of DS, which provides a reference basis for the development of GMPEs based on the nonergodic hypothesis, and regional seismic risk analysis and loss assessment.

New assessment of the probabilistic seismic hazard analysis for the greater Jakarta area, Indonesia

The greater Jakarta area is located in west Java and has a high earthquake risk due to the subduction of the Australian Plate beneath the Eurasian Plate. Previous hazard studies have been carried out using various seismic probability approaches. This study presents a new probabilistic seismic hazard analysis (PSHA) by developing an updated seismic model for the crustal active fault database, using a new set of ground motion prediction equations (GMPE), and taking into consideration the effect of site amplification. The logic tree method was used to quantify the epistemic uncertainty of the source parameter components and the attenuation equation was used. We calculate the peak ground acceleration (PGA) and spectral acceleration (SA) for West Java, including metropolitan Jakarta, Bogor and Serang, which have the greatest concentrations of population and business. We found that the mean PGA reaches 1.0 g or a 2% probability of exceedance in 50 years, including site amplification factors. The seismic hazard is dominated by the shallow crustal fault near metropolitan Jakarta. The result of this study may be useful for updating the existing Indonesian hazard map.

An Evolutionary Shaking‐Forecast‐Based Earthquake Early Warning Method

AbstractHere we propose a methodology for Earthquake Early Warning (EEW) able to issue the alert based on the real‐time estimation of the epicentral area where a peak ground Intensity measure is expected to exceed a user‐set ground shaking level. The method provides in output a P‐wave‐based, time‐evolutive “early” shake map. It combines the peak ground velocity predictions available from the observed P‐wave amplitudes and from the region‐specific ground motion prediction equation, using progressively updated estimates of earthquake location and magnitude. The P‐wave displacement, velocity and acceleration amplitudes are jointly measured on a progressively expanded P‐wave time window while the earthquake location and magnitude are evaluated using the first P‐arrival time and displacement amplitudes at near source stations. A retrospective analysis of the 2016, Mw 6.5 Central Italy earthquake records shows that depending on the network density and spatial source coverage, the method naturally accounts for effects related to the earthquake rupture directivity and spatial variability of strong ground motion related to crustal wave propagation and site amplification. Within 1.5 s from the first alert (5.15 s after the origin time), the simulated performance of the system in predicting the event ground shaking is very high: in the 40 km‐radius area that suffered an Intensity MCS VIII–IX, 41 over 42 strong‐motion instrumented sites would have been successfully alerted, with only one false alarm. Even considering the calculated blind‐zone of 15 km radius, a 15–55 km wide annular area would have received the alert 2–14.5 s before the occurrence of the strong ground shaking. The proposed EEW method evolves with time in a way that it minimizes the missed alarms while increasing successful alarms and to a lesser extent false alarms, so it is necessary for the end‐user to accept these possibilities and account for them in a probabilistic decision scheme depending on the specific safety actuation measure to be undertaken in real‐time.

Probabilistic seismic hazard assessment in Namibia

Namibia is located in Southern Africa on a Stable Continental Region (SRC). Despite the fact that a number of earthquakes have been observed or recorded in Namibia since 1910, with the July 31, 2009 and the April 4, 2021 earthquakes both of moment magnitude (Mw) of 5.6, no up-to-date country-wide seismic hazard study has been conducted. As a result, the purpose of this research is to conduct a comprehensive Probabilistic Seismic Hazard Assessment (PSHA) for the entire country in order to aid in planning and seismic risk mitigation. The earthquake catalogue developed for this study spans the years 1910–2021 and includes both historical and instrumental events. Data from the International Seismological Centre (ISC), the Advanced Network of Seismological System (ANSS), and Namibia's National Seismological Network (NSN) were used to create the catalogue. The catalogue was then declustered to remove fore- and aftershocks as well as similar events; the resulting catalogue was used in conjunction with available geological evidence to identify, delineate, and characterise 11 seismic source zones as well as one fault source, the Hebron/Dreylingen fault. Four Ground Motion Prediction Equations (GMPEs) were carefully selected and used in the seismic hazard computations. The GMPEs were implemented in the calculations using a logic tree formalism with equal weights, which assisted in addressing the uncertainties associated with both the seismic sources and the ground motion models. Peak Ground Accelerations (PGA) obtained for a 10% chance of exceeding in 50 years ranged from 0.017 g to 0.149 g. The highest levels of hazard were observed in Namibia's north-western, north-central, central and southern regions. This study also produced seismic hazard maps that show the distribution of acceleration at different response periods (0s (PGA), 0.1s, 0.15s, 0.2. s, 0.3s, 0.5s and 1.0s) computed for a 10% probability of exceedance in 50 years. The findings from this Namibia's first large-scale seismic hazard study are expected to make a significant contribution for future land use planning in facility and infrastructural development in the country.