A team of earthquake geologists, seismologists and engineering seismologists from GNS Science, NIWA, University of Canterbury, and Victoria University of Wellington have collectively produced an update of the 2002 national probabilistic seismic hazard (PSH) model for New Zealand. The new model incorporates over 200 new onshore and offshore fault sources, and utilises newly developed New Zealand-based scaling relationships and methods for the parameterisation of the fault and subduction interface sources. The background seismicity model has also been updated to include new seismicity data, a new seismicity regionalisation, and improved methodology for calculation of the seismicity parameters. Background seismicity models allow for the occurrence of earthquakes away from the known fault sources, and are typically modelled as a grid of earthquake sources with rate parameters assigned from the historical seismicity catalogue. The Greendale Fault, which ruptured during the M7.1, 4 September 2010 Darfield earthquake, was unknown prior to the earthquake. However, the earthquake was to some extent accounted for in the PSH model. The maximum magnitude assumed in the background seismicity model for the area of the earthquake is 7.2 (larger than the Darfield event), but the location and geometry of the fault are not represented. Deaggregations of the PSH model for Christchurch at return periods of 500 years and above show that M7-7.5 fault and background source-derived earthquakes at distances less than 40 km are important contributors to the hazard. Therefore, earthquakes similar to the Darfield event feature prominently in the PSH model, even though the Greendale Fault was not an explicit model input.
On 4 September 2010, a magnitude Mw 7.1 earthquake struck the Canterbury region on the South Island of New Zealand. The epicentre of the earthquake was located in the Darfield area about 40 km west of the city of Christchurch. Extensive damage occurred to unreinforced masonry buildings throughout the region during the mainshock and subsequent large aftershocks. Particularly extensive damage was inflicted to lifelines and residential houses due to widespread liquefaction and lateral spreading in areas close to major streams, rivers and wetlands throughout Christchurch and Kaiapoi. Despite the severe damage to infrastructure and residential houses, fortunately, no deaths occurred and only two injuries were reported in this earthquake. From an engineering viewpoint, one may argue that the most significant aspects of the 2010 Darfield Earthquake were geotechnical in nature, with liquefaction and lateral spreading being the principal culprits for the inflicted damage. Following the earthquake, a geotechnical reconnaissance was conducted over a period of six days (10–15 September 2010) by a team of geotechnical/earthquake engineers and geologists from New Zealand and USA (GEER team: Geo-engineering Extreme Event Reconnaissance). JGS (Japanese Geotechnical Society) members from Japan also participated in the reconnaissance team from 13 to 15 September 2010. The NZ, GEER and JGS members worked as one team and shared resources, information and logistics in order to conduct thorough and most efficient reconnaissance covering a large area over a very limited time period. This report summarises the key evidence and findings from the reconnaissance.
Research Report No.2010-03Ground motion prediction equations (GMPEs) for geometric-mean pseudo-spectral acceleration amplitudes from New Zealand (NZ) earthquakes are developed. A database of 2437 three-component ground motion records is developed by applying stringent quality criteria to the historically recorded events in NZ. Despite the large number of records, the database is deficient in empirical records from large magnitude events recorded at close distances to the fault rupture plane. As a result, the basis for the NZ-specific GMPE development is to examine the applicability of foreign GMPEs for similar tectonic regions and then modify the most applicable GMPEs based on both theoretical and statistically significant empirically-driven arguments. For active shallow crustal events, five different GMPEs are considered. It was found that the McVerry et al. (2006) model, which is the current model upon which seismic design guidelines and site-specific seismic hazard analyses in NZ are based, provided the worst fit to the NZ database, and that the Chiou et al. (2010) (C10) modification of the Chiou and Youngs (2008) model was the most applicable. Discrepancies between the C10 model and the NZ database that were empirically identified and theoretically justified were used to modify the C10 model for: (i) small magnitude scaling; (ii) scaling of short period ground motion from normal faulting events in volcanic crust; (iii) scaling of ground motions on very hard rock sites; (iv) anelastic attenuation in the NZ crust; and (v) consideration of the increased anelastic attenuation in the Taupo Volcanic Zone (TVZ). For subduction slab events, initially three models were considered. It was found that all of the models had some significant biases with respect to applicability for NZ. The Zhao et al. (2006) (Z06) model was selected because of the rigorous database upon which it was developed and modified by: (i) NZ-specific scaling at small magnitudes; (ii) path scaling at large distances; (iii) consideration of the increased TVZ attenuation; and (iv) revision of the standard deviation model. Based on these modifications the developed model showed no bias of the inter- and intra-event residuals as a function of various predictor variables. The standard deviation of the residuals using the revised standard deviation model also indicated that the model has an adequate precision. Three GMPEs were considered for subduction interface events. The Zhao et al. (2006) (Z06) model was the best performing model with only bias exhibited in the site response model, and possible over-prediction of large magnitude events. The Z06 interface model was modified to account for site response and magnitude scaling using the same functional forms as those of the developed active shallow crustal and subduction slab models. The developed model showed no bias of the inter- and intra-event residuals as a function of various predictor variables. The developed GMPEs include specific features as evident in the NZ database; consistent scaling for parameters not well constrained by the NZ database; and pseudo-spectral amplitudes for vibration periods from 0.01 to 10 seconds. Hence, these models represent a significant advance in the state-of-the art for empirical ground motion prediction in NZ.
