Understanding disaster risk enables Government, industry and the community to make better decisions on how to prepare for disasters and improve the resilience of communities. Geoscience Australia develops and provides fundamental data and information to understand disaster risk so that we can determine how hazards impact the things that are valuable to us.

With a population of over 250 million people, Indonesia is the fourth most populous country in the world (United Nations, 2013). Indonesia also experiences more earthquakes than any other country in the world (USGS, 2015). Its borders encompass one of the most active tectonic regions on Earth including over 18 000 km of major tectonic plate boundary, more than twice that of Japan or Papua New Guinea (Bird, 2003). The potential for this tectonic activity to impact large populations has been tragically demonstrated by the 20004 Sumatra earthquake and tsunami. In order to inform earthquake risk reduction in Indonesia, a new national earthquake hazard map was developed in 2010 (Irsyam et al., 2010). In this report historical records of damaging earthquakes from the 17th to 19th centuries are used to test our current understanding of earthquake hazard in Indonesia and identify areas where further research is needed. In this report we address the following questions: - How well does our current understanding of earthquake hazard in Indonesia reflect historical activity? - Can we associate major historical earthquakes with known active faults, and are these accounted for in current assessments of earthquake hazard? - Does the current earthquake hazard map predict a frequency and intensity of shaking commensurate with the historical record? - What would the impact of these historical earthquakes be if they were to reoccur today? To help answer questions like these, this report collates historical observations of eight large earthquakes from Java, Bali and Nusa Tenggara between 1699 and 1867. These observations are then used to: - Identify plausible sources for each event; - Develop ground shaking models using the OpenQuake Engine (GEM Foundation, 2015); - Assess the validity of the current national seismic hazard map; and - Estimate fatalities were the historical events to occur today using the InaSAFE (InaSAFE.org, 2015) software.

Many earthquakes in Indonesia have caused a large number of fatalities. Disaster risk-reduction of fatalities requires a representative fatality model derived from fatality data caused by historical earthquakes in Indonesia. We develop an empirical fatality model for Indonesia by relating macroseismic intensity to fatality rate using compiled subdistrict level fatality rate data and numerically simulated ground shaking intensity for four recent damaging events. The fatality rate data are compiled by collecting population and fatality statistics of the regions impacted by the selected events. The ground shaking intensity is numerically estimated by incorporating a finite fault model of each event and local site conditions approximated by topographically-based site amplifications. The macroseismic intensity distribution of each event is generated by using ShakeMap software with a selected pair of ground motion predictive equation (GMPE) and ground motion to intensity conversion equation (GMICE). The developed fatality model is a Bayesian generalized linear model where the fatality rate is assumed to follow a mixture of a Bernoulli and a gamma distribution. The probability of zero fatality rate and the mean non-zero fatality rate is linked to a linear function of shaking intensity by the logit and the log link functions, respectively. We estimate posterior distribution of the parameters of the model based on the Hamilton Monte Carlo algorithm. For validation of the developed model we calculate fatalities of the past events from the EXPO-CAT catalog and compare the estimates with the EXPO-CAT fatality records. While the developed fatality model can provide an estimate of the range of fatalities for future events it needs on-going refinement by incorporation of additional fatality rate data from past and future events.

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.

As part of the 2018 Tropical Cyclone Hazard Assessment (TCHA), we compiled the geospatial raster dataset that can be accessible to internal and external users via ArcGIS online and can be integrated for building additional geoprocessing applications. This web service gives more stable and easy access to data and interactive maps. With having separate geospatial layers for each recurrence interval- i.e. 5 through 10000 years, users can toggle between the layers and evaluate the changes in wind speed (km/hr) and potential areas at risk on the fly.

An evaluation of the likelihood of tropical cyclone-related extreme winds, incorporating local effects on wind speed.

The Government of Indonesia has committed to deploying a network of 500 strong-motion sensors throughout the nation. The data from these sensors have the potential to provide critical near-real-time information on the level of ground shaking and potential impact from Indonesian earthquakes near communities. We describe the implementation of real-time ‘ShakeMaps’ within Indonesia's Agency of Meteorology, Climatology and Geophysics (BMKG). These ShakeMaps are intended to underpin real-time earthquake situational awareness tools. The use of the new strong-motion network is demonstrated for two recent earthquakes in northern Sumatra: the 2 July 2013 Mw 6.1 Bener Meriah, Sumatra and the 10 October 2013 Mw 5.4 Aceh Besar earthquakes. The former earthquake resulted in 35 fatalities, with a further 2400 reported injuries. The recently integrated ShakeMap system automatically generated shaking estimates calibrated by BMKG's strong-motion network within 7 min of the Bener Meriah earthquake's origin, which assisted the emergency response efforts. Recorded ground motions are generally consistent with theoretical models. However, more analysis is required to fully characterize the attenuation of strong ground motion in Indonesia.

Hot emissions of mainly sulphur dioxide and carbon dioxide took place from a mound in Koranga open cut, near Wau, following a landslide at the end of May, 1967. Rocks of the Holocene volcano, Koranga, are exposed in the open cut. The emissions lasted about three months, and ceased on 13 August after another landslide removed the active mound. During the period of activity, recorded temperatures ranged up to 680°C; no anomalous seismic or tilt phenomena were recorded. The cause of the activity is not known, but it is thought that the high temperatures and gases may have been the result of the spontaneous combustion of reactive sulphides and carbonaceous material present in the altered rocks of Koranga volcano.

Probabilistic seismic hazard map of Papua New Guinea, in terms of Peak Ground Acceleration, is developed for return period of 475 years. The calculations were performed for bedrock site conditions (Vs30=760 m/s). Logic-tree framework is applied to include epistemic uncertainty in seismic source as well as ground-motion modelling processes. In this regard two source models, using area source zones and smoothed seismicity, are developed. Based on available geological and seismological data, defined seismic sources are classified into 4 different tectonic environments. For each of the tectonic regimes three Ground Motion Prediction Equations are selected and used to estimate the ground motions at a grid of sites with spacing of 0.1 degree in latitude and longitude. Results show high level of hazard in the coastal areas of Huon Peninsula and New Britain/ Bougainville regions and relatively low level of hazard in the southern part of the New Guinea highlands block. In Huon Peninsula, as shown by seismic hazard disaggregation results, high level of hazard is caused by modelled frequent moderate to large earthquakes occurring at Ramu-Markham Fault zone. On the other hand in New Britain/Bougainville region, the geometry and distance to the subduction zone along New Britain Trench mainly controls the calculated level of hazard. It is also shown that estimated level of PGAs is very sensitive to the selection of GMPEs and overall the results are closer to the results from studies using more recent ground-motion models.

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