Papua New Guinea (PNG) lies in a belt of intense tectonic activity that experiences high levels of seismicity. Although this seismicity poses significant risks to society, the Building Code of PNG and its underpinning seismic loading requirements have not been revised since 1982. This study aims to partially address this gap by updating the seismic zoning map on which the earthquake loading component of the building code is based. We performed a new probabilistic seismic hazard assessment for PNG using the OpenQuake software developed by the Global Earthquake Model Foundation (Pagani et al. 2014). Among other enhancements, for the first time together with background sources, individual fault sources are implemented to represent active major and microplate boundaries in the region to better constrain the earthquake-rate and seismic-source models. The seismic-source model also models intraslab, Wadati–Benioff zone seismicity in a more realistic way using a continuous slab volume to constrain the finite ruptures of such events. The results suggest a high level of hazard in the coastal areas of the Huon Peninsula and the New Britain – Bougainville region, and a relatively low level of hazard in the southwestern part of mainland PNG. In comparison with the seismic zonation map in the current design standard, it can be noted that the spatial distribution of seismic hazard used for building design does not match the bedrock hazard distribution of this study. In particular, the high seismic hazard of the Huon Peninsula in the revised assessment is not captured in the current building code of PNG. <b>Citation:</b> Ghasemi, H., Cummins, P., Weatherill, G. <i>et al.</i> Seismotectonic model and probabilistic seismic hazard assessment for Papua New Guinea. <i>Bull Earthquake Eng, </i><b>18</b>, 6571–6605 (2020). https://doi.org/10.1007/s10518-020-00966-1

Seismic hazard models, commonly produced through probabilistic seismic hazard analysis, are used to establish earthquake loading requirements for the built environment. However, there is considerable uncertainty in developing seismic hazard models, which require assumptions on seismicity rates and ground-motion models (GMMs) based on the best evidence available to hazard analysts. This paper explores several area-based tests of long-term seismic hazard forecasts for the Australian continent. ShakeMaps are calculated for all earthquakes of MW 4.25 and greater within approximately 200 km of the Australian coastline using the observed seismicity in the past 50 years (1970-2019). A “composite ShakeMap” is generated that extracts the maximum peak ground acceleration “observed” in this 50-year period for any site within the continent. The fractional exceedance area of this composite map is compared with four generations of Australian seismic hazard maps for a 10% probability of exceedance in 50 years (~1/500 annual exceedance probability) developed since 1990. In general, all these seismic hazard models appear to be conservative relative to the observed ground motions that are estimated to have occurred in the last 50 years. To explore aspects of possible prejudice in this study, the variability in ground-motion exceedance was explored using the Next Generation Attenuation-East GMMs developed for the central and eastern United States. The sensitivity of these results is also tested with the interjection of a rare scenario earthquake with an expected regional recurrence of approximately 5,000 - 10,000 years. While these analyses do not provide a robust assessment of the performance of the candidate seismic hazard for any given location, they do provide—to the first order—a guide to the performance of the respective maps at a continental scale. This paper was presented at the Australian Earthquake Engineering Society 2021 Virtual Conference, Nov 25 – 26.