(Australia) to Earthquake-Generated Tsunamis, Storms and Sea Level Rise
Total Page:16
File Type:pdf, Size:1020Kb
OPEN The exposure of Sydney (Australia) to SUBJECT AREAS: earthquake-generated tsunamis, storms PLANETARY SCIENCE CLIMATE-CHANGE IMPACTS and sea level rise: a probabilistic Received multi-hazard approach 27 August 2014 F. Dall’Osso1, D. Dominey-Howes1, C. Moore2, S. Summerhayes3 & G. Withycombe3 Accepted 18 November 2014 1Asia – Pacific Natural Hazards Research Group, School of Geosciences, the University of Sydney, 2NOAA Pacific Marine Published Environmental Laboratory, 3Sydney Coastal Councils Group Inc. 10 December 2014 Approximately 85% of Australia’s population live along the coastal fringe, an area with high exposure to extreme inundations such as tsunamis. However, to date, no Probabilistic Tsunami Hazard Assessments Correspondence and (PTHA) that include inundation have been published for Australia. This limits the development of requests for materials appropriate risk reduction measures by decision and policy makers. We describe our PTHA undertaken for should be addressed to the Sydney metropolitan area. Using the NOAA NCTR model MOST (Method for Splitting Tsunamis), we F.D.O. simulate 36 earthquake-generated tsunamis with annual probabilities of 15100, 151,000 and 1510,000, (filippodallosso@ occurring under present and future predicted sea level conditions. For each tsunami scenario we generate a high-resolution inundation map of the maximum water level and flow velocity, and we calculate the gmail.com) or D.D.-H. exposure of buildings and critical infrastructure. Results indicate that exposure to earthquake-generated (dale.dominey- tsunamis is relatively low for present events, but increases significantly with higher sea level conditions. The [email protected]. probabilistic approach allowed us to undertake a comparison with an existing storm surge hazard au) assessment. Interestingly, the exposure to all the simulated tsunamis is significantly lower than that for the 15100 storm surge scenarios, under the same initial sea level conditions. The results have significant implications for multi-risk and emergency management in Sydney. ustralia is at risk from tsunamis but due to its geographic proximity to different seismic source regions, the tsunami hazard around Australia varies significantly1 (Figure 1a). Since the 2004 Indian Ocean tsunami disaster great effort has been made to understand the palaeo- and historic records of tsunamis affecting A 1–3 Australia, estimate future tsunami probabilities and to assess risk . For the New South Wales (NSW) coast, the most likely sources of earthquake-generated tsunamis are the New Hebrides, Puysegur, Kermadec, Tonga and Chile Trenches1. Burbidge et al.1,4 are the only team to publish a Probabilistic Tsunami Hazard Assessment (PTHA) for Australia. Burbidge et al.1 estimated tsunami offshore wave amplitudes (i.e. at the -100m bathymetric contour) corresponding to various annual probabilities of occurrence (i.e. probability of a given wave amplitude to be reached or exceeded), triggered by all possible seismic sources around Australia. The methodology is based on the approach usually adopted for Probabilistic Seismic Hazard Assessments (PSAH). Burbidge et al.1 simulated a number of earthquakes whose probability of occurrence was known for each of the selected tsunami sources (i.e. all subduction zones surrounding Australia). Other tsunami sources such as landslides and volcanic eruptions were not considered. The vertical component of the crustal deformation of each earthquake was used as the initial condition of a finite-difference numerical model to propagate each tsunami towards the coast reaching a bathy- metric contour of -100m. This produced a maximum tsunami wave amplitude (at the -100m contour) for each earthquake, associated to its probability of occurrence. The probabilistic hazard maps of Burbidge et al.1 are given at an offshore bathymetric contour of -100m. Near- shore tsunami propagation and inundation modelling were not undertaken. To the best of our knowledge, no PTHAs that incorporate onshore inundation modelling have been published in Australia. In a pilot study, Dall’Osso et al.2 explored the exposure and vulnerability of buildings to inundation that would be associated with a tsunami run-up of 7 m in the Sydney suburb of Manly. They showed that such a scenario SCIENTIFIC REPORTS | 4 : 7401 | DOI: 10.1038/srep07401 1 www.nature.com/scientificreports Figure 1 | Study area and off-shore Probabilistic Tsunami Hazard Assessment for Australia. (a) the red circle indicates the study area location: this is part of the metropolitan area of Sydney, in which most of the NSW population is clustered. The colour bars surrounding Australia represent the off-shore tsunami wave amplitude having a 1 in 1,000 chance of being exceeded per year.; (b) the case study area located south of central Sydney, NSW and includes Botany Bay, Bate Bay and Port Hacking. These areas are located within the Councils of Botany Bay, Rockdale and Sutherland Shire. Points #1, #2 and #3 located by red stars were selected to extract the tsunami wave amplitude fields shown in Figures 2 and 3. The point coordinates are: point#1 (-33.9577; 151.1930), point#2 (-33.9877; 151.1490), point #3 (-34.0503; 151.1582) (Datum: CGS WGS84). Each point has a vertical elevation of -2m AHD. Figure 1(a) was modified from Burbidge et al.1, with permission from Geoscience Australia. E Commonwealth of Australia (Geoscience Australia) 2014. Figure 1(b) was created by the first author (FD) using ESRI ArcGIS 10. SCIENTIFIC REPORTS | 4 : 7401 | DOI: 10.1038/srep07401 2 www.nature.com/scientificreports would inundate 12001 buildings. However, Dall’Osso et al.2 esti- the estuary zone of Port Hacking (Figure 1b). The PTHA mated the tsunami inundation using a deterministic, static ‘‘bathtub includes hydrodynamic modelling of tsunami generation, trans- filling’’ method that limited the scope of their work. oceanic and near-shore propagation, and inland inundation. Historically, tide-gauge records show that only small tsunamis Using the model MOST (Method of Splitting Tsunamis)12,we have affected NSW5. However, according to the ‘‘Australian simulate 36 tsunami scenarios by combining three annual prob- Megatsunami Hypothesis’’ (AMH) large palaeotsunamis may have abilities of occurrence (i.e. 1/100, 1/1,000 and 1/10,000), two occurred repeatedly throughout the Holocene although their sources tsunami source locations (i.e. New Hebrides and Puysegur remain unknown6. Although the AMH is controversial, the fact trenches) and six initial sea level states (accounting for tide fluc- remains that the exposed NSW coast (Figure 1a) is home to approxi- tuations and sea level rise estimates). For each scenario we obtain mately 80% of the State’s population. In addition to the AMH, geo- the tsunami arrival time, the wave amplitude field, the extent of logical evidence of smaller tsunamis that occurred in the mid to late the forecast area of maximum inundation, and estimates of the Holocene has been found in various estuary locations along the maximum water level and flow velocity in each point of the study southeast coast of Australia7. Although these events would have been area. smaller than those proposed by the AMH, if repeated today, their b. Assess the exposure of buildings and infrastructure to the consequences would be significant8. selected tsunami scenarios. Exposure variations across different In fact, Chen and McAneny9 estimated that in Sydney, some tsunami scenarios are then discussed and compared based on the 20,000 properties are at risk from inundation associated with a vari- tsunami probability, source location and initial sea level condi- ety of coastal processes being located ,1 km from the shoreline and tions; at no more than 3m above sea level. However, this estimate did not c. Demonstrate how the PTHA can be used in multi-hazard stud- include tsunamis. ies. This is achieved by comparing the outputs of the PTHA and The deficiency of PTHAs coupled with dynamic numerical inland the exposure assessment with a probabilistic storm surge hazard inundation modelling and building exposure assessments make it and exposure assessment undertaken by McInnes et al.11, which difficult for state and local government authorities and the emer- used the same initial sea level conditions. gency services to implement comprehensive multi-hazard risk reduction strategies, including to the threat posed by tsunamis. The paucity of PTHAs in Australia means that the risk from tsuna- Results mis cannot be fully quantified nor compared against natural hazards For each tsunami scenario, MOST calculated the time series of whose probabilities are known. Further, in the absence of this inundation depth and flow velocity reached in each cell of the information, the NSW Government10 defines tsunami risk in broad topo-bathymetric grids used as inputs. This information was used terms as ‘the entire population living within 1 km of the shoreline, to generate the following results: and below 10m above sea level’. This approach can under or over- estimate exposure, particularly in areas featuring estuaries and ria- Tsunami wave amplitude fields. Examples of the tsunami wave like coasts, such as Sydney11. amplitude field at three selected locations (Points #1, #2 and #3 on Our study produces a more accurate PTHA to assist and underpin Figure 1b) are shown in Figure 2 and 3. As expected, the observed coastal risk management and land-use planning strategies. In par- tsunami arrival time depends primarily on the distance between the ticular, this work aims to: source and the impact area (Figures 2 and 3). It is not influenced by the earthquake magnitude or by the initial sea level condition. In each a. Undertake a PTHA for a densely populated area in south Sydney, scenario, multiple tsunami waves would reach the study area. In most NSW. The study area comprises the bays of Botany and Bate, and scenarios the first wave has the largest amplitude, although there are Figure 2 | Wave amplitude fields of tsunamis triggered in New Hebrides (north of the study area), including the 1/100 (in green), the 1/1,000 (in blue) and the 1/10,000 (in red) events.