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.

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.

The Earthquake Scenario Selection is an interactive tool for querying, visualising and downloading earthquake scenarios. There are over 160 sites nationally with pre-generated scenarios available. These represent plausible future scenarios that can be used for earthquake risk management and planning (see https://www.ga.gov.au/about/projects/safety/nsha for more details).

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.

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/)

A database of recordings from moderate-to-large magnitude earthquakes is compiled for earthquakes in western and central Australia. Data are mainly recorded by Australian National Seismograph Network (ANSN), complemented with data from temporary deployments, and covering the period of 1990 to 2019. The dataset currently contains 1497 earthquake recordings from 164 earthquakes with magnitudes from MW 2.5 to 6.1, and hypocentral distances up to 1500 km. The time-series data are consistently processed to correct for the instrument response and to reduce the effect of background noise. A range of ground-motion parameters in the time and frequency domains are calculated and stored in the database. Numerous near-source recordings exceed peak accelerations of 0.10 g and range up to 0.66 g, while the maximum peak velocity of the dataset exceeds 27 cm/s. In addition to its utility for engineering design, the dataset compiled herein will improve characterisation of ground-motion attenuation in the region and will provide an excellent supplement to ground-motion datasets collected in analogue seismotectonic regions worldwide. This paper was presented at the Australian Earthquake Engineering Society 2021 Virtual Conference, Nov 25 – 26.

Public concerns have been raised about the potential for induced seismicity as state and territory governments lift moratoriums on hydraulic stimulation activities for the exploration and extraction of unconventional hydrocarbons. The Scientific Inquiry into Hydraulic Fracturing in the Northern Territory articulated the need for a traffic-light system “to minimise the risk of occurrence of seismic events during hydraulic fracturing operations” within the Beetaloo Sub-basin. A temporary seismic network (Phase 1) was deployed in late 2019 to monitor baseline seismic activity in the basin. Based on the data analysed herein (November 2019 – April 2021), no seismic events were identified within the area of interest suggesting that the Beetaloo Sub-basin is largely aseismic. Observations to date indicate that there is potential to identify events smaller than ML=1.5 within the basin. The recent installation of ten semi-permanent stations for continuous real-time monitoring will contribute to ongoing baseline monitoring efforts and support the implementation of an induced seismicity traffic-light system. The outcome of this study will be used to build knowledge about potential human-induced seismic activity in the region that may be associated with unconventional hydrocarbon recovery. This paper was presented at the Australian Earthquake Engineering Society 2021 Virtual Conference, Nov 25 – 26.

Here we undertake a statistical analysis of local magnitudes (ML) calculated using the two real-time earthquake monitoring software platforms use by Geoscience Australia (GA) since 2005, Antelope and Seiscomp. We examine a database of just over 10 years duration, during a period in which both systems were in operation and over 4000 earthquakes were located and magnitudes estimated. We examine the consistency of both single-station and network ML estimates of both systems, with a view toward determining guidelines for combining them into a single catalogue, as well as for determining best practice in the for the estimation of local magnitudes for regions of sparse seismic networks. Once this guidance has been developed, it is the intention of GA to re-process magnitudes for all earthquakes using a consistent approach where digital data are available and can be integrated within the currently-used SeisComP system. This paper was presented at the Australian Earthquake Engineering Society 2021 Virtual Conference, Nov 25 – 26.

The Mwp 6.1 Petermann Ranges earthquake that occurred on 20 May, 2016 in the Central Ranges, NT, is the largest onshore earthquake to be recorded in Australia since the 1988 Tennant Creek sequence. While geodetic and geophysical analyses have characterized the extent of surface rupture and faulting mechanism respectively, a comprehensive aftershock characterization has yet to be performed. Data has been acquired from a 12-station temporary seismic network deployed jointly by the ANU and Geoscience Australia (GA), collected from five days following the mainshock to early October. Taking advantage of enhanced automatic detection techniques using the SeisComP3 real-time earthquake monitoring software within the National Earthquake Alerts Centre (NEAC) at GA, we have developed a comprehensive earthquake catalogue for this mainshock-aftershock sequence. Utilising the NonLinLoc location algorithm combined with a Tennant Creek-derived velocity model, we have preliminarily located over 5,800 aftershocks. With additional spatio-temporal analyses and event relocation, our objective will be to use these aftershocks to help delineate the geometry of the headwall rupture along the Woodroffe Thrust. These high-resolution aftershock detection techniques are intended to be implemented in real-time within the NEAC following future significant Australian intraplate earthquakes. This paper was presented at the Australian Earthquake Engineering Society 2021 Virtual Conference, Nov 25 – 26.