Liquefaction features and the geologic environment in which they formed were carefully studied at two sites near Lincoln in southwest Christchurch. We undertook geomorphic mapping, excavated trenches, and obtained hand cores in areas with surficial evidence for liquefaction and areas where no surficial evidence for liquefaction was present at two sites (Hardwick and Marchand). The liquefaction features identified include (1) sand blows (singular and aligned along linear fissures), (2) blisters or injections of subhorizontal dikes into the topsoil, (3) dikes related to the blows and blisters, and (4) a collapse structure. The spatial distribution of these surface liquefaction features correlates strongly with the ridges of scroll bars in meander settings. In addition, we discovered paleoliquefaction features, including several dikes and a sand blow, in excavations at the sites of modern liquefaction. The paleoliquefaction event at the Hardwick site is dated at A.D. 908-1336, and the one at the Marchand site is dated at A.D. 1017-1840 (95% confidence intervals of probability density functions obtained by Bayesian analysis). If both events are the same, given proximity of the sites, the time of the event is A.D. 1019-1337. If they are not, the one at the Marchand site could have been much younger. Taking into account a preliminary liquefaction-triggering threshold of equivalent peak ground acceleration for an Mw 7.5 event (PGA7:5) of 0:07g, existing magnitude-bounded relations for paleoliquefaction, and the timing of the paleoearthquakes and the potential PGA7:5 estimated for regional faults, we propose that the Porters Pass fault, Alpine fault, or the subduction zone faults are the most likely sources that could have triggered liquefaction at the study sites. There are other nearby regional faults that may have been the source, but there is no paleoseismic data with which to make the temporal link.
The last seven years have seen southern New Zealand a ected by several large and damaging earthquakes: the moment magnitude (MW) 7.8 Dusky Sound earthquake on 15 July 2009, the MW 7.1 Dar eld (Canterbury) earthquake on 4 September 2010, and most notably the MW 6.2 Christchurch earthquake on 22 February 2011 and the protracted aftershock sequence. In this thesis, we address the postseismic displacement produced by these earthquakes using methods of satellite-based geodetic measurement, known as Interferometric Synthetic Aperture Radar (InSAR) and Global Positioning System (GPS), and computational modelling. We observe several ground displacement features in the Canterbury and Fiordland regions during three periods: 1) Following the Dusky Sound earthquake; 2) Following the Dar eld earthquake and prior to the Christchurch earthquake; and 3) Following the Christchurch earthquake until February 2015. The ground displacement associated with postseismic motion following the Dusky Sound earthquake has been measured by continuous and campaign GPS data acquired in August 2009, in conjunction with Di erential Interferometric Synthetic Aperture Radar (DInSAR) observations. We use an afterslip model, estimated by temporal inversion of geodetic data, with combined viscoelastic rebound model to account for the observed spatio-temporal patterns of displacement. The two postseismic processes together induce a signi cant displacement corresponding to principal extensional and contractual strain rates of the order of 10⁻⁷ and 10⁻⁸ yr⁻¹ respectively, across most of the southern South Island. We also analyse observed postseismic displacement following the Dusky Sound earthquake using a new inversion approach in order to describe afterslip in an elasticviscoelastic medium. We develop a mathematical framework, namely the "Iterative Decoupling of Afterslip and Viscoelastic rebound (IDAV)" method, with which to invert temporally dense and spatially sparse geodetic observations. We examine the IDAV method using both numerical and analytical simulations of Green's functions. For the post-Dar eld time interval, postseismic signals are measured within approximately one month of the mainshock. The dataset used for the post-Dar eld displacement spans the region surrounding previously unrecognised faults that ruptured during the mainshock. Poroelastic rebound in a multi-layered half-space and dilatancy recovery at shallow depths provide a satisfactory t with the observations. For the post-Christchurch interval, campaign GPS data acquired in February 2012 to February 2015 in four successive epochs and 66 TerraSAR-X (TSX) SAR acquisitions in descending orbits between March 2011 and May 2014 reveal approximately three years of postseismic displacement. We detect movement away from the satellite of ~ 3 mm/yr in Christchurch and a gradient of displacement of ~ 4 mm/yr across a lineament extending from the westernmost end of the Western Christchurch Fault towards the eastern end of the Greendale East Fault. The postseismic signals following the Christchurch earthquake are mainly accounted for by afterslip models on the subsurface lineament and nearby faults.
The 2015 New Zealand strong-motion database provides a wealth of new strong motion data for engineering applications. An important component of this database is the compilation of new site metadata, describing the soil conditions and site response at GeoNet strong motion stations. We have assessed and compiled four key site parameters for the ~460 GeoNet stations that recorded significant historical ground motions. Parameters include: site classification (NZS1170.5), Vs30, fundamental site period (Tsite) and depth to bedrock (Z1.0, i.e. depth to material with Vs > 1000 m/s). In addition, we have assigned a quality estimate (Quality 1 – 3) to these parameters to provide a qualitative estimate of the uncertainty. New highquality Tsite estimates have largely been obtained from newly available HVSR amplification curves and spectral ratios from inversion of regional strong motion data that has been reconciled with available geological information. Good quality Vs30 estimates, typically in urban centres, have also been incorporated following recent studies. Where site-specific measurements of Vs30 are not available, Vs30 is estimated based on surface geology following national Vs30 maps. New Z1.0 values have been provided from 3D subsurface models for Canterbury and Wellington. This database will be used in efforts to guide development and testing of new and existing ground motion prediction models in New Zealand. In particular, it will allow reexamination of the most important site parameters that control and predict site response in a New Zealand setting. Furthermore, it can be used to provide information about suitable rock reference sites for seismological research, and as a guide to site-specific references in the literature. We discuss compilation of the database, preliminary insights so far, and future directions.
