Christchurch and Canterbury suffered significant housing losses due to the earthquakes. Estimates from the Earthquake Commission (EQC) (2011) suggest that over 150,000 homes (around three quarters of Christchurch housing stock) sustained damage from the earthquakes. Some areas of Christchurch have been declared not suitable for rebuilding, affecting more than 7,500 residential properties.
This paper presents a critical evaluation of vertical ground motions observed in the Canterbury earthquake sequence. The abundance of strong near-source ground-motion recordings provides an opportunity to comprehensively review the estimation of vertical ground motions via the New Zealand Standard for earthquake loading, NZS1170.5:2004, and empirical ground motion prediction equations (GMPEs). An in-depth review of current GMPEs is carried out to determine the existing trends and characteristics present in the empirical models. Results illustrate that vertical ground motion amplitudes estimated based on NZS1170.5:2004 are significantly unconservative at short periods and near-source distances. While conventional GMPEs provide an improved prediction, in many instances they too underpredict vertical ground motion accelerations at short periods and near-source distances.
Novel Gel-push sampling was employed to obtain high quality samples of Christchurch sands from the Central Business District, at sites where liquefaction was observed in 22 February 2011, and 13 June 2011 earthquakes. The results of cyclic triaxial testing on selected undisturbed specimens of typical Christchurch sands are presented and compared to empirical procedures used by practitioners. This comparison suggests cyclic triaxial data may be conservative, and the Magnitude Scaling Factor used in empirical procedures may be unconservative for highly compressible soils during near source moderate to low magnitude events. Comparison to empirical triggering curves suggests the empirical method generally estimates the cyclic strength of Christchurch sands within a reasonable degree of accuracy as a screening evaluation tool for liquefaction hazard, however for sands with moderate to high fines content it may be significantly unconservative, highlighting the need for high quality sampling and testing on important projects where seismic performance is critical.
Liquefaction-induced lateral spreading in Christchurch and surrounding suburbs during the recent Canterbury Earthquake Sequence (2010-2011) caused significant damage to structures and lifelines located in close proximity to streams and rivers. Simplified methods used in current engineering practice for predicting lateral ground displacements exhibit a high degree of epistemic uncertainty, but provide ‘order of magnitude’ estimates to appraise the hazard. We wish to compare model predictions to field measurements in order to assess the model’s capabilities and limitations with respect to Christchurch conditions. The analysis presented focuses on the widely-used empirical model of Youd et al. (2002), developed based on multi-linear regression (MLR) of case history data from lateral spreading occurrence in Japan and the US. Two issues arising from the application of this model to Christchurch were considered: • Small data set of Standard Penetration Test (SPT) and soil gradation indices (fines content FC, and mean grain size, D50) required for input. We attempt to use widely available CPT data with site specific correlations to FC and D50. • Uncertainty associated with the model input parameters and their influence on predicted displacements. This has been investigated for a specific location through a sensitivity analysis.
The collapse of Redcliffs’ cliff in the 22 February 2011 and 13 June 2011 earthquakes were the first times ever a major failure incident occurred at Redcliffs in approximately 6000 years. This master’s thesis is a multidisciplinary engineering geological investigation sought to study these particular failure incidents, focusing on collecting the data necessary to explain the cause and effect of the cliff collapsing in the event of two major earthquakes. This study provides quantitative and qualitative data about the geotechnical attributes and engineering geological nature of the sea-cut cliff located at Redcliffs. Results from surveying the geology of Redcliffs show that the exposed lithology of the cliff face is a variably jointed rock body of welded and (relatively intact) unwelded ignimbrite, a predominantly massive unit of brecciated tuff, and a covering of wind-blown loess and soil deposit (commonly found throughout Canterbury) on top of the cliff. Moreover, detailing the external component of the slope profile shows that Redcliffs’ cliff is a 40 – 80 m cliff with two intersecting (NE and SE facing) slope aspects. The (remotely) measured geometry of the cliff face comprises of multiple outstanding gradients, averaging a slope angle of ~67 degrees (post-13 June 2011), where the steepest components are ~80 degrees, whereas the gentle sloping sections are ~44 degrees. The physical structure of Redcliffs’ cliff drastically changed after each collapse, whereby seismically induced alterations to the slope geometry resulted in material deposited on the talus at the base of the cliff. Prior to the first collapse, the variance of the gradient down the slope was minimal, with the SE Face being the most variable with up to three major gradients on one cross section. However, after each major collapse, the variability increased with more parts of the cliff face having more than one major gradient that is steeper or gentler than the remainder of the slope. The estimated volume of material lost as a result of the gradient changes was 28,267 m³ in February and 11,360 m³ in June 2011. In addition, surveys of the cliff top after the failure incidents revealed the development of fissures along the cliff edge. Monitoring 10 fissures over three months indicated that fissured by the cliff edge respond to intense seismicity (generally ≥ Mw 4) by widening. Redcliffs’ cliff collapsed on two separate occasions as a result of an accumulated amount of damage of the rock masses in the cliff (caused by weathering and erosion over time), and two Mw 6.2 trigger earthquakes which shook the Redcliffs and the surrounding area at a Peak Ground Acceleration (PGA) estimated to be around 2 g. The results of the theoretical study suggests that PGA levels felt on-site during both instances of failure are the result of three major factors: source of the quake and the site affected; topographic amplification of the ground movement; the short distance between the source and the cliff for both fault ruptures; the focus of seismic energy in the direction of thrust faulting along a path that intercepts Redcliffs (and the Port Hills). Ultimately, failure on the NE and SE Faces of Redcliffs’ cliff was concluded to be global as every part of the exposed cliff face deposited a significant volume of material on the talus at the base of the cliff, with the exception of one section on the NE Face. The cliff collapses was a concurrent process that is a single (non-monotonic) event that operated as a complex series of (primarily) toppling rock falls, some sliding of blocks, and slumping of the soil mantle on top of the cliff. The first collapse had a mixture of equivalent continua slope movement of the heavily weathered / damaged surface of the cliff face, and discontinuous slope movement of the jointed inner slope (behind the heavily weathered surface); whereas the second collapse resulted in only discontinuous slope movement on account of the freshly exposed cliff face that had damage to the rock masses, in the form of old and (relatively) new discontinuous fractures, induced by earthquakes and aftershocks leading up to the point of failure.
The purpose of this thesis is to conduct a detailed examination of the forward-directivity characteristics of near-fault ground motions produced in the 2010-11 Canterbury earthquakes, including evaluating the efficacy of several existing empirical models which form the basis of frameworks for considering directivity in seismic hazard assessment. A wavelet-based pulse classification algorithm developed by Baker (2007) is firstly used to identify and characterise ground motions which demonstrate evidence of forward-directivity effects from significant events in the Canterbury earthquake sequence. The algorithm fails to classify a large number of ground motions which clearly exhibit an early-arriving directivity pulse due to: (i) incorrect pulse extraction resulting from the presence of pulse-like features caused by other physical phenomena; and (ii) inadequacy of the pulse indicator score used to carry out binary pulse-like/non-pulse-like classification. An alternative ‘manual’ approach is proposed to ensure 'correct' pulse extraction and the classification process is also guided by examination of the horizontal velocity trajectory plots and source-to-site geometry. Based on the above analysis, 59 pulse-like ground motions are identified from the Canterbury earthquakes , which in the author's opinion, are caused by forward-directivity effects. The pulses are also characterised in terms of their period and amplitude. A revised version of the B07 algorithm developed by Shahi (2013) is also subsequently utilised but without observing any notable improvement in the pulse classification results. A series of three chapters are dedicated to assess the predictive capabilities of empirical models to predict the: (i) probability of pulse occurrence; (ii) response spectrum amplification caused by the directivity pulse; (iii) period and amplitude (peak ground velocity, PGV) of the directivity pulse using observations from four significant events in the Canterbury earthquakes. Based on the results of logistic regression analysis, it is found that the pulse probability model of Shahi (2013) provides the most improved predictions in comparison to its predecessors. Pulse probability contour maps are developed to scrutinise observations of pulses/non-pulses with predicted probabilities. A direct comparison of the observed and predicted directivity amplification of acceleration response spectra reveals the inadequacy of broadband directivity models, which form the basis of the near-fault factor in the New Zealand loadings standard, NZS1170.5:2004. In contrast, a recently developed narrowband model by Shahi & Baker (2011) provides significantly improved predictions by amplifying the response spectra within a small range of periods. The significant positive bias demonstrated by the residuals associated with all models at longer vibration periods (in the Mw7.1 Darfield and Mw6.2 Christchurch earthquakes) is likely due to the influence of basin-induced surface waves and non-linear soil response. Empirical models for the pulse period notably under-predict observations from the Darfield and Christchurch earthquakes, inferred as being a result of both the effect of nonlinear site response and influence of the Canterbury basin. In contrast, observed pulse periods from the smaller magnitude June (Mw6.0) and December (Mw5.9) 2011 earthquakes are in good agreement with predictions. Models for the pulse amplitude generally provide accurate estimates of the observations at source-to-site distances between 1 km and 10 km. At longer distances, observed PGVs are significantly under-predicted due to their slower apparent attenuation. Mixed-effects regression is employed to develop revised models for both parameters using the latest NGA-West2 pulse-like ground motion database. A pulse period relationship which accounts for the effect of faulting mechanism using rake angle as a continuous predictor variable is developed. The use of a larger database in model development, however does not result in improved predictions of pulse period for the Darfield and Christchurch earthquakes. In contrast, the revised model for PGV provides a more appropriate attenuation of the pulse amplitude with distance, and does not exhibit the bias associated with previous models. Finally, the effects of near-fault directivity are explicitly included in NZ-specific probabilistic seismic hazard analysis (PSHA) using the narrowband directivity model of Shahi & Baker (2011). Seismic hazard analyses are conducted with and without considering directivity for typical sites in Christchurch and Otira. The inadequacy of the near-fault factor in the NZS1170.5: 2004 is apparent based on a comparison with the directivity amplification obtained from PSHA.
Fine grained sediment deposition in urban environments during natural hazard events can impact critical infrastructure and properties (urban terrain) leading to reduced social and economic function and potentially adverse public health effects. Therefore, clean-up of the sediments is required to minimise impacts and restore social and economic functionality as soon as possible. The strategies employed to manage and coordinate the clean-up significantly influence the speed, cost and quality of the clean-up operation. Additionally, the physical properties of the fine grained sediment affects the clean-up, transport, storage and future usage of the sediment. The goals of the research are to assess the resources, time and cost required for fine grained sediment clean-up in an urban environment following a disaster and to determine how the geotechnical properties of sediment will affect urban clean-up strategies. The thesis focuses on the impact of fine grained sediment (<1 mm) deposition from three liquefaction events during the Canterbury earthquake sequence (2010-2011) on residential suburbs and transport networks in Christchurch. It also presents how geotechnical properties of the material may affect clean-up strategies and methods by presenting geotechnical analysis of tephra material from the North Island of New Zealand. Finally, lessons for disaster response planning and decision making for clean-up of sediment in urban environments are presented. A series of semi-structured interviews of key stakeholders supported by relevant academic literature and media reports were used to record the clean-up operation coordination and management and to make a preliminary qualification of the Christchurch liquefaction ejecta clean-up (costs breakdown, time, volume, resources, coordination, planning and priorities). Further analysis of the costs and resources involved for better accuracy was required and so the analysis of Christchurch City Council road management database (RAMM) was done. In order to make a transition from general fine sediment clean-up to specific types of fine disaster sediment clean-up, adequate information about the material properties is required as they will define how the material will be handled, transported and stored. Laboratory analysis of young volcanic tephra from the New Zealand’s North Island was performed to identify their geotechnical properties (density, granulometry, plasticity, composition and angle of repose). The major findings of this research were that emergency planning and the use of the coordinated incident management system (CIMS) system during the emergency were important to facilitate rapid clean-up tasking, management of resources and ultimately recovery from widespread and voluminous liquefaction ejecta deposition in eastern Christchurch. A total estimated cost of approximately $NZ 40 million was calculated for the Christchurch City clean-up following the 2010-2011 Canterbury earthquake sequence with a partial cost of $NZ 12 million for the Southern part of the city, where up to 33% (418 km) of the road network was impacted by liquefaction ejecta and required clearing of the material following the 22 February 2011 earthquake. Over 500,000 tonnes of ejecta has been stockpiled at Burwood landfill for all three liquefaction inducing earthquake events. The average cost per kilometre for the event clean-up was $NZ 5,500/km (4 September 2010), $NZ 11,650/km (22 February 2011) and $NZ 11,185/km (13 June 2011). The duration of clean-up time of residential properties and the road network was approximately two to three months for each of the three liquefaction ejecta events; despite events volumes and spatial distribution of ejecta. Interviews and quantitative analysis of RAMM data revealed that the experience and knowledge gained from the Darfield earthquake (4 September 2010) clean-up increased the efficiency of the following Christchurch earthquake induced liquefaction ejecta clean-up events. Density, particle size, particle shape, clay content and moisture content, are the important geotechnical properties that need to be considered when planning for a clean-up method that incorporates collection, transport and disposal or storage. The geotechnical properties for the tephra samples were analysed to increase preparedness and reaction response of potentially affected North Island cities from possible product from the active volcanoes in their region. The geotechnical results from this study show that volcanic tephra could be used in road or construction material but the properties would have to be further investigated for a New Zealand context. Using fresh volcanic material in road, building or flood control construction requires good understanding of the material properties and precaution during design and construction to extra care, but if well planned, it can be economically beneficial.
