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

Tsunami hazard maps are generated for the Mentawai Islands, West Sumatra, Indonesia, to support evacuation and disaster response planning. A random heterogeneous slip generator is used to forward model a suite of earthquake rupture scenarios on the Mentawai Segment of the Sunda Subduction Zone. A total of 1000 rupture models that fit constraints provided by coral and geodetic records of coseismic vertical deformation from great earthquakes in 1797, 1833 and 2007 are used to model inundation and define a maximum inundation zone that envelopes all of these scenarios. Results are compared with single scenario hazard assessments developed by experts and agreed through scientific consensus building processes to assess the additional value of modelling a suite of scenarios to obtain a more robust estimate of potential inundated areas by incorporating uncertainty in the earthquake source. The model presented here, like all tsunami hazard assessments, is based on assumptions about the characteristics of future events based on past events, however by sampling a range of plausible outcomes we gain a more robust estimate of which areas may be inundated during a tsunami within the bounds of our assumptions.

In June 2012 Geoscience Australia was commissioned by Commonwealth Scientific and Industrial Research Organisation (CSIRO) to undertake detailed wind hazard assessments for 14 Pacific Island countries and East Timor as part of the Pacific-Australia Climate Change Science and Adaptation Planning (PACCSAP) program. PACCSAP program follows on from work Geoscience Australia did for the Pacific Climate Change Science Program (PCCSP) looking at CMIP3 generation of climate models. The objective of this study is to improve scientific knowledge by examining past climate trends and variability to provide regional and national climate projections. This document presents results from current and future climate projections of severe wind hazard from tropical cyclones for the 15 PACCSAP partner countries describing the data and methods used for the analysis. The severe wind hazard was estimated for current (1981 to 2000) and future (2081 to 2100) climate scenarios. Tropical-cyclone like vortices from climate simulations conducted by CSIRO using six Coupled Model Intercomparison Project phase 5 (CMIP5) models (BCC-CSM1.1, NorESM1-M, CSIRO-Mk3.6, IPSL-CM5A, MRI-CGM3 and GFDL-ESM2M) as well as the International Best Track Archive for Climate Stewardship were used as input to the Geoscience Australia's Tropical Cyclone Risk Model to generate return period wind speeds for the 15 PACCSAP partner countries. The Tropical Cyclone Risk Model is a statistical-parametric model of tropical cyclone behaviour, enabling users to generate synthetic records of tropical cyclones representing many thousands of years of activity. The 500-year return period wind speed is analysed and discussed into more details in this report, since it is used as a benchmark for the design loads on residential buildings. Results indicate that there is not a consistent spatial trend for the changes in 500-year cyclonic wind speed return period when CMIP5 models are compared individually. BCC-CSM1M and IPSL-CM5A presented an increase in the annual TC frequency for East Timor, northern hemisphere and southern hemisphere. On the other hand, NorESM1M showed a decrease in the annual TC frequency for the same areas. The other three models showed a mixed of increase and decrease in their annual TC frequency. When CMIP5 models were analysed by partner county capitals for the 500-year cyclonic wind speed return period, IPSL-CM5A and GFDL-ESM2M models presented an increase in the cyclonic wind speed intensity for almost all capitals analysed with exception of Funafuti (GFDL-ESM2M), which presented a decrease of 0.7% and Honiara (IPSL-CM5A) with a decrease of 1.6%. The tropical cyclone annual frequency ensemble mean indicates an increase in the tropical cyclone frequency within all three regions considered in this study. When looking at individual capitals, a slight increase in the 500-year return period cyclonic wind speed ensemble mean varying between 0.8% (Port Vila) to 9.1% (Majuro) is noticed. A decline around 2.4% on average in the 500-year return period cyclonic wind speed ensemble mean is observed in Dili, Suva, Nukualofa and Ngerulmud. The ensemble spatial relative change did not show any particular consistency for the 500-year cyclonic wind speed. Areas where Marshall Islands and Niue are located presented an increase in the 500-year cyclonic wind speed while a decrease is observed in areas around South of Vanuatu, East of Solomon Islands, South of Fiji and some areas in Tonga. The information from the evaluation of severe wind hazard from tropical cyclones, together with other PACCSAP program outputs, will be used to build partner country capacity to effectively adapt and plan for the future and overcome challenges from climate change.

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

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.

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.

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. 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 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 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 building code of PNG. We will also discuss how the seismic hazard map of PNG is being used to underpin its building code, including what steps have been taken by GA together with the Government of PNG to promote uptake of the new hazard map by PNG’s earthquake engineering community.

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.

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.

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