Through Australian Department of Foreign Affairs and Trade, Geoscience Australia has been working closely with the Government of Papua New Guinea technical agencies (Rabaul Volcano Observatory, Port Moresby Geophysical Observatory, and Engineering Geology Branch) since September 2010 to enhance their capabilities to monitor and assess natural hazards. The objective of this program is to support the Government of Papua New Guinea in developing fundamental information and practices for the effective response and management of natural hazard events in PNG. Earthquakes as natural hazards are one of the key focus points of this project, as they continue to cause loss of life and widespread damage to buildings and infrastructure in Papua New Guinea. The country’s vulnerability to earthquakes is evident from the significant socio-economic consequences of recent major events in Papua New Guinea, e.g., a magnitude 7.5 earthquake that occurred in the Hela Province of Papua New Guinea in 2018. Earthquake risk is likely to increase significantly in the years to come due to the growth in population and urbanization in Papua New Guinea. However, earthquake risk, unlike hazard, can be managed and minimized. One obvious example would be minimizing earthquake risk by constructing earthquake-resistant structures following building standards. The high level of earthquake hazard of Papua New Guinea has been long recognised and the suite of building standards released in 1982 contained provisions to impart adequate resilience to buildings based on the best understanding of seismic hazard available at that time. However, the building standards and incorporated seismic hazard assessment for Papua New Guinea has not been updated since the 1980s. The integration of modern national seismic hazard models into national building codes and practices provides the most effective way that we can reduce human casualties and economic losses from future earthquakes. This report aims at partially fulfilling this task by performing a probabilistic seismic hazard assessment to underpin a revision of the earthquake loading component of the building standards of Papua New Guinea. The updated assessment offers many important advances over its predecessor. It is based on a modern probabilistic hazard framework and considers an earthquake catalogue augmented with an additional four decades-worth of data. The revised assessment considers advances in ground-motion modelling through the use of multiple ground-motion models. Also, for the first time, the individual fault sources representing active major and microplate boundaries are implemented in the input hazard model. Furthermore, the intraslab sources are represented realistically by using the continuous slab volume to constrain the finite ruptures of such events. This would better constrain the expected levels of ground motion at any given site in Papua New Guinea. 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 southern part of the New Guinea Highlands Block. In comparison with the seismic zonation map in the current design standard, it can be noted that the spatial distribution 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 seismic zoning map, leading to a significant under-estimation of hazard in PNG’s second-largest city, Lae. It can also be shown that in many other regions and community localities in PNG the hazard is higher than that regulated for the design of buildings having a range of natural periods. Thus, the need for an updated hazard map for building design has been confirmed from the results of this study, and a revised map is developed for consideration in a revised building standard of Papua New Guinea.

In November, 2018 a workshop of experts sponsored by UNESCO’s Intergovernmental Oceanographic Commission was convened in Wellington, New Zealand. The meeting was organized by Working Group (WG) 1 of the Pacific Tsunami Warning System (PTWS). The meeting brought together fourteen experts from various disciplines and four different countries (New Zealand, Australia, USA and French Polynesia) and four observers from Pacific Island countries (Tonga, Fiji), with the objective of understanding the tsunami hazard posed by the Tonga-Kermadec trench, evaluating the current state of seismic and tsunami instrumentation in the region and assessing the level of readiness of at-risk populations. The meeting took place in the “Beehive” Annex to New Zealand’s Parliament building nearby the offices of the Ministry of Civil Defence and Emergency Management. The meeting was co-chaired by Mrs. Sarah-Jayne McCurrach (New Zealand) from the Ministry of Civil Defence and Emergency Management and Dr. Diego Arcas (USA) from NOAA’s Pacific Marine Environmental Laboratory. As one of the meeting objectives, the experts used their state-of-the-science knowledge of local tectonics to identify some of the potential, worst-case seismic scenarios for the Tonga-Kermadec trench. These scenarios were ranked as low, medium and high probability events by the same experts. While other non-seismic tsunamigenic scenarios were acknowledged, the level of uncertainty in the region, associated with the lack of instrumentation prevented the experts from identifying worse case scenarios for non-seismic sources. The present report synthesizes some of the findings of, and presents the seismic sources identified by the experts to pose the largest tsunami risk to nearby coastlines. In addition, workshop participants discussed existing gaps in scientific knowledge of local tectonics, including seismic and tsunami instrumentation of the trench and current level of tsunami readiness for at-risk populations, including real-time tsunami warnings. The results and conclusions of the meeting are presented in this report and some recommendations are summarized in the final section.

The Assessment of Tropical Cyclone Risks in the Pacific Region project represents a collaboration between DIICCSRTE and Geoscience Australia with PCRAFI and AIR Worldwide. Building on the expertise of each organisation, the project will deliver an assessment of the financial risks to buildings, infrastructure and agriculture arising from tropical cyclones (TCs) under current and future climate regimes. This extends previous risk assessments undertaken by incorporating the influence of climate change on the hazard (TCs) into the assessment process. The output of this study is a set of peril matrices, which detail the relative change in parameters describing TC behaviour: e.g. annual mean frequency, mean maximum intensity and mean latitude of genesis. The relative changes are evaluated as the fractional change between TC behavior in current climate GCM simulations and future climate GCM simulations.

Natural Hazards and Earth Systems Science

The present study reports on recent developments of the Indonesia Tsunami Early Warning System (InaTEWS), especially with respect to the tsunami modeling components used in that system. It is a dual system: firstly, InaTEWS operates a high-resolution scenario database pre-computed with the finite element model TsunAWI; running in parallel, the system also contains a supra real-time modeling component based on the GPU-parallelized linear long-wave model easyWave capable of dealing with events outside the database coverage. The evolution of the tsunami scenario database over time is covered in the first sections. Starting from the mere coverage of the Sunda Arc region, the current state contains scenarios in 15 fault zones. The study is augmented by an investigation of warning products used for early warning like the estimated wave height (EWH) and the estimated time of arrival (ETA). These quantities are determined by easyWave and TsunAWI with model specific approaches. Since the numerical setup of the models is very different, the extent of variations in warning products is investigated for a number of scenarios, where both pure database scenarios and applications to real events are considered.

Heterogeneous distribution of slip during megathrust earthquakes has been shown to significantly affect the spatial distribution of tsunami height in both numerical studies and field observations. This means that tsunami hazard maps generated using uniform slip distributions in their tsunami source models may underestimate tsunami inundation in some locations compared with real events of the same magnitude in the same location. In order to more completely define areas that may be inundated during a tsunami it is important to consider how different possible distributions of slip will impact different parts of the coastline. We generate tsunami inundation maps for the Mentawai Islands, West Sumatra, Indonesia, from a composite suite of possible source models that are consistent with current knowledge of the source region. First, a suite of earthquake source models with randomly distributed slip along the Mentawai Segment of the Sunda Subduction Zone is generated using a k-2 rupture model. From this suite we select source models that generate vertical deformation consistent with that observed in coral palaeogeodetic records of previous ruptures of the Mentawai Segment in 1797 and 1833, minus deformation observed in the 2007 Bengkulu earthquake sequence. Tsunami inundation is then modelled using high resolution elevation data for selected source models and the results compiled to generate a maximum tsunami inundation zone. This method allows us to constrain the slip distribution beneath the Mentawai Islands, where coral palaeogeodetic data is available, while allowing for greater variation in the slip distribution away from the islands, in particular near the trench where large slip events can generate very large tsunami. This method also allows us to consider high slip events on deeper portions of the megathrust between the Mentawai Islands and the Sumatran Mainland, which give greater tsunami inundation on the eastern part of the Mentawai Islands and the west coast of Sumatra compared with near-trench event. By accounting for uncertainty in slip distribution, the resulting hazard maps give a more complete picture of the areas that may be inundated compared with hazard maps derived from a single 'worst case' source model. These maps allow for more robust tsunami evacuation plans to be developed to support immediate community evacuation in response to strong or long-lasting earthquake ground shaking. From the American Geophysical Union Fall Meeting Abstracts

A component of the PNG-Australia Volcanological Services Support (VSS) Project funded by AusAID

Historical reports of earthquake effects from the period 1681 to 1877 in Java, Bali and Nusa Tenggara are used to independently test ground motion predictions in Indonesia’s 2010 national probabilistic seismic hazard assessment (PSHA). Assuming that strong ground motion occurrence follows a Poisson distribution, we cannot reject Indonesia’s current PSHA for key cities in Java at 95% confidence. However, the results do suggest that seismic hazard may be underestimated for the megacity Jakarta. Ground motion simulations for individual large damaging events are used to identify plausible source mechanisms, providing insights into the major sources of earthquake hazard in the region and possible maximum magnitudes for these sources. The results demonstrate that large intraslab earthquakes have been responsible for major earthquake disasters in Java, including a ~Mw 7.5 intraslab earthquake near Jakarta in 1699 and a ~Mw 7.8 event in 1867 in Central Java. The results also highlight the potential for large earthquakes to occur on the Flores Thrust. We require an earthquake with Mw 8.4 on the Flores Thrust to reproduce tsunami observation from Sulawesi and Sumbawa in 1820. Furthermore, large shallow earthquakes (Mw > 6) have occurred in regions where active faults have not been mapped identifying the need for further research to identify and characterize these faults for future seismic hazard assessments. <b>Citation:</b> Jonathan Griffin, Ngoc Nguyen, Phil Cummins, Athanasius Cipta; Historical Earthquakes of the Eastern Sunda Arc: Source Mechanisms and Intensity‐Based Testing of Indonesia’s National Seismic Hazard Assessment. <i>Bulletin of the Seismological Society of America </i>2018; 109 (1): 43–65. doi: https://doi.org/10.1785/0120180085

<div>The city of Lae is Papua New Guinea (PNG)’s second largest, and is the home of PNG’s largest port. Here, a convergence rate of ~50 mm/yr between the South Bismarck Plate and the Australian Plate is accommodated across the Ramu-Markham Fault Zone (RMFZ). The active structures of the RMFZ are relatively closely spaced to the west of Lae. However, the fault zone bifurcates immediately west of the Lae urban area, with one strand continuing to the east, and a second strand trending southeast through Lae City and connecting to the Markham Trench within the Huon Gulf. </div><div>The geomorphology of the Lae region relates to the interaction between riverine (and limited marine) deposition and erosion, and range-building over low-angle thrust faults of the RMFZ. Flights of river terraces imply repeated tectonic uplift events; dating of these terraces will constrain the timing of past earthquakes and associated recurrence intervals. Terrace riser heights are typically on the order of 3 m, indicating causative earthquake events of greater than magnitude 7. </div><div>Future work will expose the most recently active fault traces in trenches to assess single event displacements, and extend the study to the RMFZ north of Nadzab Airport. These results will inform a seismic hazard and risk assessment for Lae city and surrounding region.</div> Presented at the 2023 Australian Earthquake Engineering Society (AEES) Conference

Tsunami hazard modelling for Tonga shows the potential impacts of tsunami generated by a very large earthquake on the nearby Tongan Trench.