As part of the 2018 National Seismic Hazard Assessment (NSHA), we compiled the geographic information system (GIS) dataset to enable end-users to view and interrogate the NSHA18 outputs on a spatially enabled platform. It is intended to ensure the NSHA18 outputs are openly available, discoverable and accessible to both internal and external users. This geospatial product is derived from the dataset generated through the development of the NSHA18 and contains uniform probability hazard maps for a 10% and 2% chance of exceedance in 50 years. These maps are calculated for peak ground acceleration (PGA) and a range of response spectral periods, Sa(T), for T = 0.1, 0.2, 0.3, 0.5, 1.0, 2.0 and 4.0 s. Additionally, hazard curves for each ground-motion intensity measure as well as uniform hazard spectra at the nominated exceedance probabilities are calculated for key localities.

Seismic hazard modelling is a multi-disciplinary science that aims to forecast earthquake occurrence and its resultant ground shaking. Such seismic hazard models consist of a probabilistic framework that models the flow of uncertainty across a complex system; typically, this includes at least two model-components developed from earth science: seismic source models, and ground motion prediction models. Although there is no scientific prescription for the length of the forecasting time-window, the most common probabilistic seismic hazard analyses (PSHA hereafter) consider forecasting probabilities of ground shaking in time windows of 30 to 50 years. These types of models are the target of this review paper. Although the core methods and assumptions of such a modelling have largely remained unchanged since they were first developed more than 50 years ago, we will review the most recent initiatives which are facing the difficult task of meeting both the increasingly sophisticated demands of society and keeping pace with advances in our scientific understanding. A need for more accurate and precise hazard forecasting must be balanced with increased quantification of uncertainty and new challenges such as moving from time-independent hazard to forecasts that are time-dependent and specific to the time-period of interest. Meeting these challenges requires the development of science-driven models which integrate at best all information available, the adoption of proper mathematical frameworks to quantify the different types of uncertainties in the source and ground motion components of the hazard model, and the development of a proper testing phase of the hazard model to quantify the consistency and skill of the hazard model. We review the state-of-the-art of the national seismic hazard modeling, and how the most innovative approaches try to address future challenges.

The Philippine archipalego is tectonically complex and seismically hazardous, yet few seismic hazard assessments have provided national coverage. This paper presents an updated probabilistic seismic hazard analysis for the nation. Active shallow crustal seismicity is modeled by faults and gridded point sources accounting for spatially variable occurrence rates. Subduction interfaces are modelled with faults of complex geometry. Intraslab seismicity is modeled by ruptures filling the slab volume. Source geometries and earthquake rates are derived from seismicity catalogs, geophysical datasets, and historic-to-paleoseismic constraints on fault slip rates. The ground motion characterization includes models designed for global use, with partial constraint by residual analysis. Shallow crustal faulting near metropolitan Manila, Davao, and Cebu dominates shaking hazard. In a few places, peak ground acceleration with 10% probability of exceedance in 50 years on rock reaches 1.0 g. The results of this study may assist in calculating the design base shear in the National Structural Code of the Philippines.

The 22 September 2021 (AEST) moment magnitude MW 5.9 Woods Point earthquake was the largest in the state of Victoria’s recorded history. The ground motions were felt throughout the state of Victoria and into neighbouring states New South Wales and South Australia. Minor damage was reported in the city of Melbourne and in some regional towns close to the epicentre. This event was captured on many high-quality recorders from multiple sources, including private, university, and public stations. These recordings provide a rare opportunity to test the validity of some ground motion models thought to be applicable to the southeast region of Australia. This paper presents spectral acceleration and attenuation comparisons of the Woods Point earthquake event to some ground motion models. The results of this paper provide further evidence that the attenuation characteristics of southeastern Australia may be similar to that in central and eastern United States, particularly at shorter distances to the epicentre. This paper was presented at the Australian Earthquake Engineering Society 2021 Virtual Conference, Nov 25 – 26.

The local magnitude ML 5.4 (MW 5.1) Moe earthquake on 19 June 2012 that occurred within the Australian stable continental region was the largest seismic event for the state of Victoria for more than 30 years. Seismic networks in the southeast Australian region yielded many high-quality recordings of the moderate-magnitude earthquake mainshock and its largest aftershock (ML 4.4; MW 4.3) at a hypocentral range of 10 to 480 km. The source and attenuation characteristics of the earthquake sequence are analyzed. Almost 15,000 felt reports were received following the main shock, which tripped a number of coal-fired power generators in the region, amounting to the loss of approximately 1955 megawatts of generation capacity. The attenuation of macroseismic intensities are shown to mimic the attenuation shape of Eastern North America (ENA) models, but require an inter-event bias to reduce predicted intensities. Further instrumental ground-motion recordings are compared to ground-motion models (GMMs) considered applicable for the southeastern Australian (SEA) region. Some GMMs developed for ENA and for SEA provide reasonable estimates of the recorded ground motions of spectral acceleration within epicentral distances of approximately 100 km. The mean weighted of the Next Generation Attenuation-East GMM suite, recently developed for stable ENA, performs relatively poorly for the 2012 Moe earthquake sequence, particularly for short-period accelerations.

Geoscience Australia has produced an Atlas of Australian earthquake scenarios (the Atlas) to support planning and preparedness operations for emergency management agencies. The Atlas provides earthquake scenarios represent realistic “worst-case” events that may impact population centres around Australia. Such scenarios may also support seismic risk assessments for critical infrastructure assets to inform remediation actions that could be taken to improve resilience to rare seismic events in Australia. The Atlas of seismic scenarios uses the underlying science and data of the 2018 National Seismic Hazard Assessment (NSHA18) to identify the magnitudes and epicentre locations of these hypothetical earthquakes. Locations and magnitudes of earthquake scenarios are based upon deaggregation of the NSHA18 hazard model. The USGS ShakeMap software is used to produce ground motion intensity fields with the shaking levels being modified by seismic site conditions mapped at a national scale. Fault sources are incorporated into the Atlas where the magnitude of a given scenario exceeds a threshold magnitude of 6.0 and where the rupture length is likely to be longer than 10 km. If a scenario earthquake is located near a known fault within the Australian Neotectonic Features database, a partial or full-length rupture is modelled along the mapped fault. The Atlas generated two scenarios for each of the160 localities across Australia. The scenarios are based on some of the most likely earthquake magnitude-distance combinations estimated at each site. Output products include shaking contours for a range of intensity measures, including peak acceleration and velocity, as well as response spectral acceleration for 0.3, 1.0 and 3.0 seconds. Also included are raster images and the associated metadata used for generating the scenarios.

We present earthquake ground motions based upon a paleoseismically-validated characteristic earthquake scenario for the ~ 48 km-long Avonmore scarp, which overlies the Meadow Valley Fault, east of Bendigo, Victoria. The results from the moment magnitude MW 7.1 scenario earthquake indicate that ground motions are sufficient to be of concern to nearby mining and water infrastructure. Specifically, the estimated median peak ground acceleration (PGA) exceeds 0.5 g to more than ~ 10 km from the source fault, and a 0.09 g PGA liquefaction threshold is exceeded out to approximately 50-70 kilometres. Liquefaction of susceptible materials, such as mine tailings, may occur to much greater distances. Our study underscores the importance of identifying and characterising potentially active faults in proximity to high failure-consequence dams, including mine tailings dams, particularly in light of the requirement to manage tailing dams for a prolonged period after mine closure. Paper presented at Australian National Committee on Large Dams (ANCOLD) conference 2020, online. (https://leishman.eventsair.com/ancold-2020-online/)

Prior to the development of Australian-specific magnitude formulae, the 1935 magnitude corrections by Charles Richter – originally developed for southern California – was almost exclusively used to calculate earthquake magnitudes throughout Australia prior to the 1990s. Due to the difference in ground-motion attenuation between southern California and much of Australia, many historical earthquake magnitudes are likely to be overestimated in the Australian earthquake catalogue. A method has been developed that corrects local magnitudes using the difference between the original (inappropriate) magnitude corrections and the Australian-specific corrections at a distance determined by the nearest recording station likely to have recorded the earthquake. These corrections have reduced the rates of local magnitudes of 4.5 in the historical catalogue by about 30% since 1900, while the number of magnitude 5.0 earthquakes has reduced by about 60% in the same time period. The reduction in the number of moderate-to-large-magnitude earthquakes over the instrumental period yields long-term earthquake rates that are more consistent with present-day rates, since the development of Australian-specific magnitude formulae. The adjustment of historical earthquake magnitudes is important for seismic hazard assessments, which assume a Poisson distribution of earthquakes in space and time.

We present a methodology for developing fragilities for mainshock-damaged structures, "aftershock fragility", by performing incremental dynamic analysis (IDA) with a sequence of mainshock-aftershock ground motions. The aftershock fragility herein is distinguished from a conventional fragility for an intact structure. We estimate seismic response of a mainshock-damaged building by performing nonlinear time history analysis with a sequence of mainshock and aftershock ground motions (so-called "back-to-back" dynamic analysis). We perform the back-to-back dynamic analyses for a number of levels of mainshock response/damage, and a number of sequences of mainshock and aftershock ground motions. With estimated seismic responses from the back-to-back dynamic analyses, we compute various damage state transition probabilities, the probability of exceeding a higher damage state from an aftershock given a damage state due to a mainshock. For an illustration of the methodology, we develop an aftershock fragility for a typical New Zealand 5-storey reinforced concrete moment frame building. The building is modeled using a single-degree-of-freedom (SDOF) damped nonlinear oscillator with force-deformation behavior represented by a multi-linear capacity/pushover curve with moderate pinching hysteresis and medium cyclic deterioration.

When multiple earthquakes occur within a short period of time, damage may accumulate in a building, affecting its ability to withstand future ground shaking. This study aims to quantify the post-earthquake capacity of a nonductile 4-story concrete building in New Zealand through incremental dynamic analysis of a nonlinear multipledegree-of-freedom simulation model. Analysis results are used to compute fragility curves for the intact and damaged buildings, showing that extensive damage reduces the structure’s capacity to resist seismic collapse by almost 30% percent. The damage experienced by the building in mainshock, can be compared with the ATC-20 building tagging criteria for post-earthquake inspections, the purpose of which is to ensure public safety. Extensively damaged buildings, which are likely be red tagged, pose a significant safety hazard due to decreased strength in future earthquakes. The effect of mainshock damage is also compared for multiple and simplified single-degree-of-freedom models of the same building.