Natural Hazards and Earth Systems Science

Reports of bitumen stranding on the ocean beaches of southern Australia date back to the early days of European settlement. Previous investigations have shown that this ‘coastal bitumen’ comprises three categories of stranded petroleum: waxy bitumen, asphaltite and oil slicks. All three varieties are physically and chemically distinct from each other, and bear no geochemical resemblance to any indigenous Australian crude oil. This study focuses on the most common variety, waxy bitumen, which accounted for 90% of the strandings on six South Australian beaches repeatedly surveyed during 1991–1992. Geochemical analysis of 96 individual specimens collected from these survey sites and other beaches in South Australia and western Victoria has shown them to be variously weathered high-wax crude oils of paraffinic to aromatic-intermediate bulk composition. Elemental, isotopic and biomarker differences allow their assignment to at least five oil families with inferred source facies that range from deep freshwater lacustrine through paludal and deltaic to euxinic marine, possibly deposited within different sedimentary basins. Family 1, 2 and 3 waxy bitumens all contain biomarkers derived from the freshwater alga Botryococcus sp. and tropical angiosperms (notably dipterocarps). Similar biomarker assemblages are unknown in Australian sedimentary basins but are common in Cenozoic crude oils and source rocks throughout western Indonesia. Family 4 waxy bitumens lack these biomarkers, but do contain dinosterane and 24-n-propylcholestane, indicative of a marine source affinity, while the carbon isotopic signatures and high pristane/phytane (Pr/Ph) ratios of Family 5 waxy bitumens are consistent with their origin from coal-rich source rocks deposited in fluvial to deltaic sedimentary successions. The majority of these waxy bitumens represent an oceanic influx of non-indigenous, Southeast Asian crude oils carried into the waters of southern Australia by the Leeuwin Current. Although they are likely to originate from natural seepage within the Indonesian Archipelago, it is unknown whether the parent oils emanate from submarine seeps or from inland seepages which are then carried to the sea by rivers. The common practice of tanker cleaning operations in the Java and Banda seas may augment the supply of natural bitumen to the beaches of Australia.

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

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.

Indonesia is one of the most disaster prone countries in the world due to its hazard profile and high population exposure. Despite its risk profile, disaster management has not traditionally been informed by best available information. Since 2008, the Australian and Indonesian governments have partnered to increase the use of science and technology in Indonesia to support decision making in disaster management. Our partnership has concentrated on strengthening the evidence base for informed disaster management by improving: - hazard information for earthquake, tsunami, volcano and flood - spatial data for exposure (population, building, roads and infrastructure) - decision support tools such as InaSAFE that assist disaster managers to combine hazard and exposure data to inform disaster response and management. We have worked alongside technical and disaster management agencies, universities, non-government organisations and the private sector to develop Indonesian capacity to manage disasters and institutionalise best-practices. Partnerships have facilitated science-to-policy and science-to-programming outcomes in disaster management that help people prepare for, respond to and recover from crises. Ten years is a good time for a partnership to form, blossom and deliver effective and sustainable changes. The achievements over this time are too many to list in entirety. Science and Technology for Disaster Management merely scratches the surface to highlight the most significant achievements. There is a tendency in doing so to focus on achievements over the last three years. In most cases the achievements of this later period of the program have only been possible because they have built on and extended the achievements of the former seven years. Science and Technology for Disaster Management demonstrates how our collaboration has increased the use of science and technology in Indonesian disaster management by developing new knowledge and enhancing capacity, both within the scientists and within policy and decision makers. In doing so, tangible policies have been developed and implemented and practices have changed. The program has helped to strengthen relationships between agencies within Indonesia and also between Indonesia and Australia in the areas of disaster management and science more broadly. New relationships manifest as distinct government-to-government collaborations and strong peer-to-peer links.

On the 30th September 2009 a magnitude 7.6 earthquake struck West Sumatra in the Padang and Pariaman regions. It caused widespread damage to buildings and resulted and an estimated 1,117 fatalities. Thankfully the event was not accompanied by a tsunami that could have had additional devastating impacts and a greatly increased mortality. Under its mandate the AIFDR responded to the earthquake event with the objective of deriving an understanding of the factors that had contributed to outcome. It supported a team of Indonesian and international engineers and scientists who collected and analysed damage information that could subsequently be used for future disaster risk reduction in West Sumatra and Indonesia more broadly. The activity was jointly led by the Centre for Disaster Mitigation at the Institut Teknologi Bandung (ITB) and Geoscience Australia. This report provides a background to the region, describes the nature of the earthquake and its impacts, details the survey activity and outlines the significant outcomes that has come from it. Importantly, it makes several recommendations to assist in the regional reconstruction after the event and to guide future development in the Padang region and Indonesia more generally.

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.

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.

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.