This final paper for the session presents the results of the new draft earthquake hazard assessment for Australia and compares them to the previous AS1170.4 hazard values. Draft hazard maps will be presented for several spectral periods (0.0, 0.2 and 1.0 s) at multiple return periods (500, 2500 and 10,000 years). These maps will be compared with both the current earthquake hazard used in AS1170.4 and with other assessments of earthquake hazard in Australia. In general the hazard in the draft map is higher in the western cratonic parts of Australia than it is in the eastern non-cratonic parts of Australia. Where regional source zones are included, peaks in hazard values in the map are generally comparable to those in the current AS1170.4 map. When seismicity 'hotspot zones are included, as described in the previous paper, several of them produce much higher hazard peaks than any in the AS1170.4 map. However, such hotspots do not affect as large an area as many of those in the current AS1170.4 map. Finally, hazard curves for different cities will also be presented and compared to those predicted by the method outlined in AS1170.4.

Stations on the Australian continent receive a rich mixture of ambient seismic noise from the surrounding oceans and the numerous small earthquakes in the earthquake belts to the north in Indonesia, and east in Tonga-Kermadec, as well as more distant source zones. The noise field at a seismic station contains information about the structure in the vicinity of the site, and this can be exploited by applying an autocorrelation procedure to the continuous records. By creating stacked autocorrelograms of the ground motion at a single station, information on crust properties can be extracted in the form of a signal that includes the crustal reflection response convolved with the autocorrelation of the combined effect of source excitation and the instrument response. After applying suitable high pass filtering the reflection component can be extracted to reveal the most prominent reflectors in the lower crust, which often correspond to the reflection at the Moho. Because the reflection signal is stacked from arrivals from a wide range of slownesses, the reflection response is somewhat diffuse, but still sufficient to provide useful constraints on the local crust beneath a seismic station. Continuous vertical component records from 223 stations (permanent and temporary) across the continent have been processed using autocorrelograms of running windows 6 hours long with subsequent stacking. A distinctive pulse with a time offset between 8 and 30 s from zero is found in the autocorrelation results, with frequency content between 1.5 and 4 Hz suggesting P-wave multiples trapped in the crust. Synthetic modelling, with control of multiple phases, shows that a local Ppmp phase can be recovered with the autocorrelation approach. This approach can be used for crustal property extraction using just vertical component records, and effective results can be obtained with temporary deployments of just a few months.

Legacy product - no abstract available

The year was another marginally below average one with respect to the frequency and magnitude of earthquakes in Australia; there were two of magnitude 5 or more but no large earthquakes of magnitude 6 or more. It was an above average year for major earthquakes worldwide with one great earthquake of maptude M 8.2 and 16 of magnitude M 7.0 or more.

The Mw=7.8 earthquake of 15 July, 2009 occurred along a section of the subduction zone south of New Zealand, where the Puysegur Block subducts beneath the Pacific Plate. The orientation of this subduction zone suggests that tsunamis generated along it pose a significant threat to the southeast coast of Australia, but since it had not experienced megathrust rupture until the 15 July event, the question of whether it was accumulating strain energy whose release could result in a large tsunami was open. We have used seismic, tsunami, geodetic and SAR data to study this earthquake and find that it involved primarily thrust motion on a fault plane dipping east at a shallow angle, consistent with expectations for a megathrust earthquake. The ability to use multiple data types to study this earthquake lead to improved ability to resolve parameters such as rupture velocity that are often difficult to constrain with seismic data alone. Seismic array data agree with rupture modelling of broadband waveforms in their prediction of a bilateral component to the earthquake rupture. Also, a tsunami of about 10 cm peak-to-peak amplitude was recorded by two tsunameter buoys in the Tasman Sea west of the epicenter, and we find that the tsunami travel times indicated by these data suggest the earthquake was characterised by a low rupture velocity of around 1 km/s. We will also present comparisons against GPS and InSAR data that further constrain parameters of the rupture. Finally, we will discuss the potential for earthquake activity further south along the Puysegur Trench, which poses a tsunami threat particularly to the eastern coast of Tasmania.

The Australian Seismological Report 2010 provides a summary of earthquake activity for Australia for 2010. It also provides a summary of earthquakes of Magnitude 5+ in the Australian Region, as well as an summary of Magnitude 6+ earthquakes worldwide. It has dedicated state and territory earthquake information including: largest earthquakes in the year; largest earthquakes in the state; and tables detailing all earthquakes detected by Geoscience Australia during the year. There are also contributions from PIRSA describing Seismic Networks and providing Earthquake locations for South Australia.

The Australian Seismological Report 2009 provides a summary of earthquake activity for Australia for 2009. It also provides a summary of earthquakes of Magnitude 5+ in the Australian Region, as well as an summary of Magnitude 6+ earthquakes worldwide. It has dedicated state and territory earthquake information including: largest earthquakes in the year; largest earthquakes in the state; and tables detailing all earthquakes detected by Geoscience Australia during the year. There are also contributions from PIRSA describing Seismic Networks and providing Earthquake locations.

Geoscience Australia (GA) is currently undertaking the process to update the Australian National Earthquake Hazard Map using modern methods and an extended catalogue of Australian earthquakes. This map is a key component of Australia's earthquake loading code. The characterisation of strong ground-shaking using Ground-Motion Prediction Equations (GMPEs) underpins any earthquake hazard assessment. We have recently seen many advances in ground-motion modelling for active tectonic regions. However, the challenge for Australia - as it is for other stable continental regions - is that there are very few ground-motion recordings from large-magnitude earthquakes with which to develop empirically-based GMPEs. Consequently, we need to consider other numerical techniques to develop these models in the absence of these data. Recently published Australian-specific GMPEs which employ these numerical techniques are now available and these will be integrated into GA's future hazard outputs. This paper addresses several fundamental aspects related to ground-motion in Australia that are necessary to consider in the update of the National Earthquake Hazard Map, including: 1) a summary of recent advances of ground-motion modelling in Australia; 2) a comparison of Australian GMPEs against those commonly used in other stable continental regions; 3) a comparison of new GMPEs against their intensity-based counterparts used in the previous hazard map; and 4) the impact of updated attenuation factors on local magnitudes in southeastern Australia. Specific regional and temporal aspects of magnitude calculation techniques across Australia and its affects on the earthquake catalogue will also be addressed.

Regional seismic reflection data in hard rock areas contains more shallow information than might first be supposed. Here I use a subset of the 2005 Tanami Seismic Survey data to show that near surface features can be defined, including paleochannels, Palaeozoic basins and structures within the Proterozoic basement. Successful imaging depends on correct determination of refraction statics, including identification of refractor branches, and use of a floating or intermediate datum during seismic reflection processing. Recognition of steep stacking velocity gradients associated with surface referenced processing aids velocity analysis and can further delineate areas of thicker regolith in palaeochannels. The first arrival refraction analysis can also be applied in more detail to estimating thickness of regolith and depth to economic basement in areas of sedimentary cover.

This paper describes the methods used to define earthquake source zones and calculate their recurrence parameters (a, b, Mmax). These values, along with the ground motion relations, effectively define the final hazard map. Definition of source zones is a highly subjective process, relying on seismology and geology to provide some quantitative guidance. Similarly the determination of Mmax is often subjective. Whilst the calculation of a and b is quantitative, the assumptions inherent in the available methods need to be considered when choosing the most appropriate one. For the new map we have maximised quantitative input into the definition of zones and their parameters. The temporal and spatial Poisson statistical properties of Australia's seismicity, along with models of intra-plate seismicity based on results from neotectonic, geodetic and computer modelling studies of stable continental crust, suggest a multi-layer source zonation model is required to account for the seismicity. Accordingly we propose a three layer model consisting of three large background seismicity zones covering 100% of the continent, 25 regional scale source zones covering ~50% of the continent, and 44 hotspot zones covering 2% of the continent. A new algorithm was developed to calculate a and b. This algorithm was designed to minimise the problems with both the maximum likelihood method (which is sensitive to the effects of varying magnitude completeness at small magnitudes) and the least squares regression method (which is sensitive to the presence of outlier large magnitude earthquakes). This enabled fully automated calculation of a and b parameters for all sources zones. The assignment of Mmax for the zones was based on the results of a statistical analysis of neotectonic fault scarps.