The city of Christchurch and its surrounds experienced widespread damage due to soil liquefaction induced by seismic shaking during the Canterbury earthquake sequence that began in September 2010 with the Mw7.1 Darfield earthquake. Prior to the start of this sequence, the city had a large network of strong motion stations (SMSs) installed, which were able to record a vast database of strong ground motions. This paper uses this database of strong ground motion recordings, observations of liquefaction manifestation at the ground surface, and data from a recently completed extensive geotechnical site investigation program at each SMS to assess a range of liquefaction evaluation procedures at the four SMSs in the Christchurch Central Business District (CBD). In general, the characteristics of the accelerograms recorded at each SMS correlated well with the liquefaction evaluation procedures, with low liquefaction factors of safety predicted at sites with clear liquefaction identifiers in the ground motions. However, at sites that likely liquefied at depth (as indicated by evaluation procedures and/or inferred from the characteristics of the recorded surface accelerograms), the presence of a non-liquefiable crust layer at many of the SMS locations prevented the manifestation of any surface effects. Because of this, there was not a good correlation between surface manifestation and two surface manifestation indices, the Liquefaction Potential Index (LPI) and the Liquefaction Severity Number (LSN).
This thesis is concerned with springs that appeared in the Hillsborough, Christchurch during the 2010-2011 Canterbury Earthquake Sequence, and which have continued to discharge groundwater to the surface to the present time. Investigations have evolved, measurements of discharge at selected sites, limited chemical data on anions and isotope analysis. The springs are associated with earthquake generated fissures (extensional) and compression zones, mostly in loess-colluvium soils of the valley floor and lower slopes. Extensive peat swamps are present in the Hillsborough valley, with a groundwater table at ~1m below ground. The first appearance of the ‘new’ springs took place following the Mw 7.1 Darfield Earthquake on 4 September 2010, and discharges increased both in volume and extent of the Christchurch Mw 6.3 Earthquake of 22 February 2011. Five monitored sites show flow rates in the range of 4.2-14.4L/min, which have remained effectively constant for the duration of the study (2014-2015). Water chemistry analysis shows that the groundwater discharges are sourced primarily from volcanic bedrocks which underlies the valley at depths ≤50m below ground level. Isotope values confirm similarities with bedrock-sourced groundwater, and the short term (hours-days) influence of extreme rainfall events. Cyclone Lusi (2013-2014) affects were monitored and showed recovery of the bedrock derived water signature within 72 hours. Close to the mouth of the valley sediments interfinger with Waimakiriri River derived alluvium bearing a distinct and different isotope signature. Some mixing is evident at certain locations, but it is not clear if there is any influence from the Huntsbury reservoir which failed in the Port Hills Earthquake (22 February 2011) and stored groundwater from the Christchurch artesian aquifer system (Riccarton Gravel).
In order to provide information related to seismic vulnerability of non-ductile reinforced concrete (RC) frame buildings, and as a complementary investigation on innovative feasible retrofit solutions developed in the past six years at the University of Canterbury on pre-19170 reinforced concrete buildings, a frame building representative of older construction practice was tested on the shake table. The specimen, 1/2.5 scale, consists of two 3-storey 2-bay asymmetric frames in parallel, one interior and one exterior, jointed together by transverse beams and floor slabs. The as-built (benchmark) specimen was first tested under increasing ground motion amplitudes using records from Loma Prieta Earthquake (California, 1989) and suffered significant damage at the upper floor, most of it due to lap splices failure. As a consequence, in a second stage, the specimen was repaired and modified by removing the concrete in the lap splice region, welding the column longitudinal bars, replacing the removed concrete with structural mortar, and injecting cracks with epoxy resin. The modified as-built specimen was then tested using data recorded during Darfield (New Zealand, 2010) and Maule (Chile, 2010) Earthquakes, with whom the specimen showed remarkably different responses attributed to the main variation in frequency content and duration. In this contribution, the seismic performance of the three series of experiments are presented and compared.
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.
The structure and geomorphology of active orogens evolves on time scales ranging from a single earthquake to millions of years of tectonic deformation. Analysis of crustal deformation using new and established remote sensing techniques, and integration of these data with field mapping, geochronology and the sedimentary record, create new opportunities to understand orogenic evolution over these timescales. Timor Leste (East Timor) lies on the northern collisional boundary between continental crust from the Australian Plate and the Banda volcanic arc. GPS studies have indicated that the island of Timor is actively shortening. Field mapping and fault kinematic analysis of an emergent Pliocene marine sequence identifies gentle folding, overprinted by a predominance of NW-SE oriented dextral-normal faults and NE-SW oriented sinistral-normal faults that collectively bound large (5-20km2) bedrock massifs throughout the island. These fault systems intersect at non-Andersonian conjugate angles of approximately 120° and accommodate an estimated 20 km of orogen-parallel extension. Folding of Pliocene rocks in Timor may represent an early episode of contraction but the overall pattern of deformation is one of lateral crustal extrusion sub-parallel to the Banda Arc. Stratigraphic relationships suggest that extrusion began prior to 5.5 Ma, during and after initial uplift of the orogen. Sedimentological, geochemical and Nd isotope data indicate that the island of Timor was emergent and shedding terrigenous sediment into carbonate basins prior to 4.5 Ma. Synorogenic tectonic and sedimentary phases initiated almost synchronously across much of Timor Leste and <2 Myr before similar events in West Timor. An increase in plate coupling along this obliquely converging boundary, due to subduction of an outlying continental plateau at the Banda Trench, is proposed as a mechanism for uplift that accounts for orogen-parallel extension and early uplift of Timor Leste. Rapid bathymetric changes around Timor are likely to have played an important role in evolution of the Indonesian Seaway. The 2010 Mw 7.1 Darfield (Canterbury) earthquake in New Zealand was complex, involving multiple faults with strike-slip, reverse and normal displacements. Multi-temporal cadastral surveying and airborne light detection and ranging (LiDAR) surveys allowed surface deformation at the junction of three faults to be analyzed in this study in unprecedented detail. A nested, localized restraining stepover with contractional bulging was identified in an area with the overall fault structure of a releasing bend, highlighting the surface complexities that may develop in fault interaction zones during a single earthquake sequence. The earthquake also caused river avulsion and flooding in this area. Geomorphic investigations of these rivers prior to the earthquake identify plausible precursory patterns, including channel migration and narrowing. Comparison of the pre and post-earthquake geomorphology of the fault rupture also suggests that a subtle scarp or groove was present along much of the trace prior to the Darfield earthquake. Hydrogeology and well logs support a hypothesis of extended slip history and suggests that that the Selwyn River fan may be infilling a graben that has accumulated late Quaternary vertical slip of <30 m. Investigating fault behavior, geomorphic and sedimentary responses over a multitude of time-scales and at different study sites provides insights into fault interactions and orogenesis during single earthquakes and over millions of years of plate boundary deformation.
The Mw 7.1 Darfield earthquake generated a ~30 km long surface rupture on the Greendale Fault and significant surface deformation related to related blind faults on a previously unrecognized fault system beneath the Canterbury Plains. This earthquake provided the opportunity for research into the patterns and mechanisms of co-seismic and post-seismic crustal deformation. In this thesis I use multiple across-fault EDM surveys, logic trees, surface investigations and deformation feature mapping, seismic reflection surveying, and survey mark (cadastral) re-occupation using GPS to quantify surface displacements at a variety of temporal and spatial scales. My field mapping investigations identified shaking and crustal displacement-induced surface deformation features south and southwest of Christchurch and in the vicinity of the projected surface traces of the Hororata Blind and Charing Cross Faults. The data are consistent with the high peak ground accelerations and broad surface warping due to underlying reverse faulting on the Hororata Blind Fault and Charing Cross Fault. I measured varying amounts of post-seismic displacement at four of five locations that crossed the Greendale Fault. None of the data showed evidence for localized dextral creep on the Greendale Fault surface trace, consistent with other studies showing only minimal regional post-seismic deformation. Instead, the post-seismic deformation field suggests an apparent westward translation of northern parts of the across-fault surveys relative to the southern parts of the surveys that I attribute to post-mainshock creep on blind thrusts and/or other unidentified structures. The seismic surveys identified a deformation zone in the gravels that we attribute to the Hororata Blind Fault but the Charing Cross fault was not able to be identified on the survey. Cadastral re-surveys indicate a deformation field consistent with previously published geodetic data. We use this deformation with regional strain rates to estimate earthquake recurrence intervals of ~7000 to > 14,000 yrs on the Hororata Blind and Charing Cross Faults.
The Canterbury Earthquakes of 2010-2011, in particular the 4th September 2010 Darfield earthquake and the 22nd February 2011 Christchurch earthquake, produced severe and widespread liquefaction in Christchurch and surrounding areas. The scale of the liquefaction was unprecedented, and caused extensive damage to a variety of man-made structures, including residential houses. Around 20,000 residential houses suffered serious damage as a direct result of the effects of liquefaction, and this resulted in approximately 7000 houses in the worst-hit areas being abandoned. Despite the good performance of light timber-framed houses under the inertial loads of the earthquake, these structures could not withstand the large loads and deformations associated with liquefaction, resulting in significant damage. The key structural component of houses subjected to liquefaction effects was found to be their foundations, as these are in direct contact with the ground. The performance of house foundations directly influenced the performance of the structure as a whole. Because of this, and due to the lack of research in this area, it was decided to investigate the performance of houses and in particular their foundations when subjected to the effects of liquefaction. The data from the inspections of approximately 500 houses conducted by a University of Canterbury summer research team following the 4th September 2010 earthquake in the worst-hit areas of Christchurch were analysed to determine the general performance of residential houses when subjected to high liquefaction loads. This was followed by the detailed inspection of around 170 houses with four different foundation types common to Christchurch and New Zealand: Concrete perimeter with short piers constructed to NZS3604, concrete slab-on-grade also to NZS3604, RibRaft slabs designed by Firth Industries and driven pile foundations. With a focus on foundations, floor levels and slopes were measured, and the damage to all areas of the house and property were recorded. Seven invasive inspections were also conducted on houses being demolished, to examine in more detail the deformation modes and the causes of damage in severely affected houses. The simplified modelling of concrete perimeter sections subjected to a variety of liquefaction-related scenarios was also performed, to examine the comparative performance of foundations built in different periods, and the loads generated under various bearing loss and lateral spreading cases. It was found that the level of foundation damage is directly related to the level of liquefaction experienced, and that foundation damage and liquefaction severity in turn influence the performance of the superstructure. Concrete perimeter foundations were found to have performed most poorly, suffering high local floor slopes and being likely to require foundation repairs even when liquefaction was low enough that no surface ejecta was seen. This was due to their weak, flexible foundation structure, which cannot withstand liquefaction loads without deforming. The vulnerability of concrete perimeter foundations was confirmed through modelling. Slab-on-grade foundations performed better, and were unlikely to require repairs at low levels of liquefaction. Ribraft and piled foundations performed the best, with repairs unlikely up to moderate levels of liquefaction. However, all foundation types were susceptible to significant damage at higher levels of liquefaction, with maximum differential settlements of 474mm, 202mm, 182mm and 250mm found for concrete perimeter, slab-on-grade, ribraft and piled foundations respectively when subjected to significant lateral spreading, the most severe loading scenario caused by liquefaction. It was found through the analysis of the data that the type of exterior wall cladding, either heavy or light, and the number of storeys, did not affect the performance of foundations. This was also shown through modelling for concrete perimeter foundations, and is due to the increased foundation strengths provided for heavily cladded and two-storey houses. Heavy roof claddings were found to increase the demands on foundations, worsening their performance. Pre-1930 concrete perimeter foundations were also found to be very vulnerable to damage under liquefaction loads, due to their weak and brittle construction.
This study contains an evaluation of the seismic hazard associated with the Springbank Fault, a blind structure discovered in 1998 close to Christchurch. The assessment of the seismic hazard is approached as a deterministic process in which it is necessary to establish: 1) fault characteristics; 2) the maximum earthquake that the fault is capable of producing and 3) ground motions estimations. Due to the blind nature of the fault, conventional techniques used to establish the basic fault characteristics for seismic hazard assessments could not be applied. Alternative methods are used including global positioning system (GPS) surveys, morphometric analyses along rivers, shallow seismic reflection surveys and computer modelling. These were supplemented by using multiple empirical equations relating fault attributes to earthquake magnitude, and attenuation relationships to estimate ground motions in the near-fault zone. The analyses indicated that the Springbank Fault is a reverse structure located approximately 30 km to the northwest of Christchurch, along a strike length of approximately 16 km between the Eyre and Ashley River. The fault does not reach the surface, buy it is associated with a broad anticline whose maximum topographic expression offers close to the mid-length of the fault. Two other reverse faults, the Eyrewell and Sefton Faults, are inferred in the study area. These faults, together with the Springbank and Hororata Faults and interpreted as part of a sys of trust/reverse faults propagating from a decollement located at mid-crustal depths of approximately 14 km beneath the Canterbury Plains Within this fault system, the Springbank Fault is considered to behave in a seismically independent way, with a fault slip rate of ~0.2 mm/yr, and the capacity of producing a reverse-slip earthquake of moment magnitude ~6.4, with an earthquake recurrence of 3,000 years. An earthquake of the above characteristics represents a significant seismic hazard for various urban centres in the near-fault zone including Christchurch, Rangiora, Oxford, Amberley, Kaiapoi, Darfield, Rollestion and Cust. Estimated peak ground accelerations for these towns range between 0.14 g to 0.5 g.
Documenting earthquake-induced ground deformation is significant to assess the characteristics of past and contemporary earthquakes and provide insight into seismic hazard. This study uses airborne light detection and ranging (LiDAR) and conducts multi-disciplinary field techniques to document the surface rupture morphology and evaluate the paleoseismicity and seismic hazard parameters of the Hurunui segment of the Hope Fault in the northern South Island of New Zealand. It also documents and evaluates seismically induced features and ground motion characteristics of the 2010 Darfield and 2011 Christchurch earthquakes in the Port Hills, south of Christchurch. These two studies are linked in that they investigate the near-field coseismic features of large (Mw ~7.1) earthquakes in New Zealand and produce data for evaluating seismic hazards of future earthquakes. In the northern South Island of New Zealand, the Australian-Pacific plate boundary is characterised by strike-slip deformation across the Marlborough Fault System (MFS). The ENE-striking Hope Fault (length: ~230 km) is the youngest and southernmost fault in the MFS, and the second fastest slipping fault in New Zealand. The Hope Fault is a major source of seismic hazard in New Zealand and has ruptured (in-part) historically in the Mw 7.1 1888 Amuri earthquake. In the west, the Hurunui segment of the Hope Fault is covered by beech forest. Hence, its seismic hazard parameters and paleoearthquake chronology were poorly constrained and it was unknown whether the 1888 earthquake ruptured this segment or not and if so, to what extent. Utilising LiDAR and field data, a 29 km-long section of the Hurunui segment of the Hope Fault is mapped. LiDAR-mapping clearly reveals the principal slip zone (PSZ) of the fault and a suite of previously unrecognised structures that form the fault deformation zone (FDZ). FDZ width measurements from 415 locations reveal a spatially-variable, active FDZ up to ~500 m wide with an average width of 200 m. Kinematic analysis of the fault structures shows that the Hurunui segment strikes between 070° and 075° and is optimally oriented for dextral strike-slip within the regional stress field. This implies that the wide FDZ observed is unlikely to result from large-scale fault mis-orientation with respect to regional stresses. The analysis of FDZ width indicates that it increases with increased hanging wall topography and increased topographic relief suggesting that along-strike topographic perturbations to fault geometry and stress states increase fault zone complexity and width. FDZ width also increases where the tips of adjacent PSZ strands locally vary in strike, and where the thickness of alluvial deposits overlying bedrock increases. LiDAR- and photogrammetrically-derived topographic mapping indicates that the boundary between the Hurunui and Hope River segments is characterised by a ~850-m-wide right stepover and a 9º-14° fault bend. Paleoseismic trenching at Hope Shelter site reveals that 6 earthquakes occurred at A.D. 1888, 1740-1840, 1479-1623, 819-1092, 439-551, and 373- 419. These rupture events have a mean recurrence interval of ~298 ± 88 yr and inter-event times ranging from 98 to 595 yrs. The variation in the inter-event times is explained by (1) coalescing rupture overlap from the adjacent Hope River segment on to the Hurunui segment at the study site, (2) temporal clustering of large earthquakes on the Hurunui segment, and/or (3) ‘missing’ rupture events. It appears that the first two options are more plausible to explain the earthquake chronologies and rupture behaviour on the Hurunui segment, given the detailed nature of the geologic and chronologic investigations. This study provides first evidence for coseismic multi-segment ruptures on the Hope Fault by identifying a rupture length of 44-70 km for the 1888 earthquake, which was not confined to the Hope River segment (primary source for the 1888 earthquake). LiDAR data is also used to identify and measure dextral displacements and scarp heights from the PSZ and structures within the FDZ along the Hurunui segment. Reconstruction of large dextrally-offset geomorphic features shows that the vertical component of slip accounts for only ~1% of the horizontal displacements and confirms that the fault is predominantly strike-slip. A strong correlation exists between the dextral displacements and elevations of geomorphic features suggesting the possibility of age correlation between the geomorphic features. A mean single event displacement (SED) of 3.6 ± 0.7 m is determined from interpretation of sets of dextral displacements of ≤ 25 m. Using the available surface age data and the cumulative dextral displacements from Matagouri Flat, McKenzie Fan, Macs Knob and Hope River sites, and the mean SED, a mean slip rate of 12.2 ± 2.4 mm/yr, and a mean recurrence interval of ~320 ± 120 yr, and a potential earthquake magnitude of Mw 7.2 are determined for the Hurunui segment. This study suggests that the fault slip rate has been constant over the last ~15000 yr. Strong ground motions from the 2010 Darfield (Canterbury) earthquake displaced boulders and caused ground damage on some ridge crests in the Port Hills. However, the 2011 Christchurch earthquake neither displaced boulders nor caused ground damage at the same ridge crests. Documentation of locations (~400 m a.s.l.), lateral displacements (8-970 cm), displacement direction (250° ± 20°) of displaced boulders, in addition to their hosting socket geometries (< 1 cm to 50 cm depth), the orientation of the ridges (000°-015°) indicate that boulders have been displaced in the direction of instrumentally recorded transient peak ground horizontal displacements nearby and that the seismic waves have been amplified at the study sites. The co-existence of displaced and non-displaced boulders at proximal sites suggests small-scale ground motion variability and/or varying boulder-ground dynamic interactions relating to shallow phenomena such as variability in soil depth, bedrock fracture density and/or microtopography on the bedrock-soil interface. Shorter shaking duration of the 2011 Christchurch event, differing frequency contents and different source characteristics were all factors that may have contributed to generating circumstances less favourable to boulder displacement in this earthquake. Investigating seismically induced features, fault behaviour, site effects on the rupture behaviour, and site response to the seismic waves provides insights into fault rupture hazards.
Liquefaction-induced lateral spreading in large seismic events often results in pervasive and costly damage to engineering structures and lifelines, making it a critical component of engineering design. However, the complex nature of this phenomenon leads to designing for such a hazard extremely challenging and there is a clear for an improved understanding and predicting liquefaction-induced lateral spreading. The 2010-2011 Canterbury (New Zealand) Earthquakes triggered severe liquefaction-induced lateral spreading along the streams and rivers of the Christchurch region, causing extensive damage to roads, bridges, lifelines, and structures in the vicinity. The unfortunate devastation induced from lateral spreading in these events also rendered the rare opportunity to gain an improved understanding of lateral spreading displacements specific to the Christchurch region. As part of this thesis, the method of ground surveying was employed following the 4 September 2010 Darfield (Mw 7.1) and 22 February 2011 Christchurch (Mw 6.2) earthquakes at 126 locations (19 repeated) throughout Christchurch and surrounding suburbs. The method involved measurements and then summation of crack widths along a specific alignment (transect) running approximately perpendicular to the waterway to indicate typically a maximum lateral displacement at the bank and reduction of the magnitude of displacements with distance from the river. Rigorous data processing and comparisons with alternative measurements of lateral spreading were performed to verify results from field observations and validate the method of ground surveying employed, as well as highlight the complex nature of lateral spreading displacements. The welldocumented field data was scrutinized to gain an understanding of typical magnitudes and distribution patterns (distribution of displacement with distance) of lateral spreading observed in the Christchurch area. Maximum displacements ranging from less than 10 cm to over 3.5 m were encountered at the sites surveyed and the area affected by spreading ranged from less than 20 m to over 200 m from the river. Despite the highly non-uniform displacements, four characteristic distribution patterns including large, distributed ground displacements, block-type movements, large and localized ground displacements, and areas of little to no displacements were identified. Available geotechnical, seismic, and topographic data were collated at the ground surveying sites for subsequent analysis of field measurements. Two widely-used empirical models (Zhang et al. (2004), Youd et al. (2002)) were scrutinized and applied to locations in the vicinity of field measurements for comparison with model predictions. The results indicated generally poor correlation (outside a factor of two) with empirical predictions at most locations and further validated the need for an improved, analysis- based method of predicting lateral displacements that considers the many factors involved on a site-specific basis. In addition, the development of appropriate model input parameters for the Youd et al. (2002) model led to a site-specific correlation of soil behavior type index, Ic, and fines content, FC, for sites along the Avon River in Christchurch that matched up well with existing Ic – FC relationships commonly used in current practice. Lastly, a rigorous analysis was performed for 25 selected locations of ground surveying measurements along the Avon River where ground slope conditions are mild (-1 to 2%) and channel heights range from about 2 – 4.5 m. The field data was divided into categories based on the observed distribution pattern of ground displacements including: large and distributed, moderate and distributed, small to negligible, and large and localized. A systematic approach was applied to determine potential critical layers contributing to the observed displacement patterns which led to the development of characteristic profiles for each category considered. The results of these analyses outline an alternative approach to the evaluation of lateral spreading in which a detailed geotechnical analysis is used to identify the potential for large spreading displacements and likely spatial distribution patterns of spreading. Key factors affecting the observed magnitude and distribution of spreading included the thickness of the critical layer, relative density, soil type and layer continuity. It was found that the large and distributed ground displacements were associated with a thick (1.5 – 2.5 m) deposit of loose, fine to silty sand (qc1 ~4-7 MPa, Ic 1.9-2.1, qc1n_cs ~50-70) that was continuous along the bank and with distance from the river. In contrast, small to negligible displacements were characterized by an absence of or relatively thin (< 1 m), discontinuous critical layer. Characteristic features of the moderate and distributed displacements were found to be somewhere between these two extremes. The localized and large displacements showed a characteristic critical layer similar to that observed in the large and distributed sites but that was not continuous and hence leading to the localized zone of displacement. The findings presented in this thesis illustrate the highly complex nature of lateral displacements that cannot be captured in simplified models but require a robust geotechnical analysis similar to that performed for this research.
This study uses 44 high quality liquefaction case histories taken from 22 locations affected by the 2010-2011 Canterbury earthquake sequence to evaluate four commonly used CPT-VS correlations (i.e., Robertson, 2009; Hegazy and Mayne, 2006; Andrus et al., 2007; McGann et al., 2015b). Co-located CPT soundings and VS profiles, developed from surface wave testing, were obtained at 22 locations and case histories were developed for the Mw 7.1, 4 September 2010 Darfield and Mw 6.2, 22 February 2011 Christchurch earthquakes. The CPT soundings are used to generate VS profiles using each of four CPT-VS correlations. These correlated VS profiles are used to estimate the factor of safety against liquefaction using the Kayen et al. (2013) VS-based simplified liquefaction evaluation procedure. An error index is used to quantify the predictive capabilities of these correlations in relation to the observations of liquefaction (or the lack thereof). Additionally, the error indices from the CPT-correlated VS profiles are compared to those obtained using: (1) the Kayen et al. (2013) procedure with surface wave-derived VS profiles, and (2) the Idriss and Boulanger (2008) CPT-based liquefaction evaluation procedure. Based on the error indices, the evaluation procedures based on direct measurements of either CPT or VS provided more accurate liquefaction triggering estimates than those obtained from any of the CPT-VS correlations. However, the performance of the CPT-VS correlations varied, with the Robertson (2009) and Hegazy and Mayne (2006) correlations performing relatively poorly for the Christchurch soils and the Andrus et al. (2007) and McGann et al. (2015b) correlations performing better. The McGann et al. (2015b) correlation had the lowest error indices of the CPT-VS correlations tested, however, none of the CPT-VS correlations provided accurate enough VS predictions to be used for the evaluation of liquefaction triggering using the VS-based liquefaction evaluation procedures.
Liquefaction during the 4th September 2010 Mw 7.1 Darfield earthquake and large aftershocks in 2011 (Canterbury earthquake sequence, CES) caused severe damage to land and infrastructure within Christchurch, New Zealand. Approximately one third of the total CES-induced financial losses were directly attributable to liq- uefaction and thus highlights the need for local and regional authorities to assess liquefaction hazards for present and future developments. This thesis is the first to conduct paleo-liquefaction studies in eastern Christchurch for the purpose of de- termining approximate return times of liquefaction-inducing earthquakes within the region. The research uncovered evidence for pre-CES liquefaction dated by radiocarbon and cross-cutting relationships as post-1660 to pre-1905. Additional paleo-liquefaction investigations within the eastern Christchurch suburb of Avon- dale, and the northern township of Kaiapoi, revealed further evidence for pre-CES liquefaction. Pre-CES liquefaction in Avondale is dated as post-1321 and pre-1901, while the Kaiapoi features likely formed during three distinct episodes: post-1458 and possibly during the 1901 Cheviot earthquake, post-1297 to pre-1901, and pre-1458. Evaluation of the liquefaction potential of active faults within the Can- terbury region indicates that many faults have the potential to cause widespread liquefaction within Avondale and Kaiapoi. The identification of pre-CES liquefac- tion confirms that these areas have previously liquefied, and indicates that residen- tial development in eastern Christchurch between 1860 and 2005 occurred in areas containing geologic evidence for pre-CES liquefaction. Additionally, on the basis of detailed field and GIS-based mapping and geospatial-statistical analysis, the distribution and severity of liquefaction and lateral spreading within the eastern Christchurch suburb of Avonside is shown in this study to be strongly in uenced by geomorphic and topographic variability. This variability is not currently ac- counted for in site-specific liquefaction assessments nor the simplified horizontal displacement models, and accounts for some of the variability between the pre- dicted horizontal displacements and those observed during the CES. This thesis highlights the potential applications of paleo-liquefaction investigations and ge- omorphic mapping to seismic and liquefaction hazard assessments and may aid future land-use planning decisions.
This report presents an overview of the soil profile characteristics at a number of strong motion station (SMS) sites in Christchurch and its surrounds. An extensive database of ground motion records has been captured by the SMS network in the Canterbury region during the Canterbury earthquake sequence. However in order to comprehensively understand the ground motions recorded at these sites and to be able to relate these motions to other locations, a detailed understanding of the shallow geotechnical profile at each SMS is required. The original NZS1170.5 (SNZ 2004) site subsoil classifications for each SMS site is based on regional geological information and well logs located at varying distances from the site. Given the variability of Christchurch soils, more detailed investigations are required in close vicinity to each SMS to better understand stratigraphy and soil properties, which are important in seismic site response. In this regard, CPT, SPT and borehole data, shear wave velocity (Vs) profiles, and horizontal to vertical spectral ratio measurements (H/V) in close vicinity to the SMS were used to develop representative soil profiles at each site. NZS1170.5 (SNZ 2004) site subsoil classifications were updated using Vs and SPT N60 criteria. Site class E boundaries were treated as a sliding scale rather than as a discrete boundary to account for locations with similar site effects potential, an approach which was shown to result in a better delineation between the site classes. SPT N60 values often indicate a stiffer site class than the Vs data for softer soil sites, highlighting the disparity between the two site investigation techniques. Both SPT N60 and Vs based site classes did not always agree with the original site classifications. This emphasises the importance of having detailed site‐specific information at SMS locations in order to properly classify them. Furthermore, additional studies are required to harmonize site classification based on SPT N60 and Vs. Liquefaction triggering assessments were carried out for the Darfield and Christchurch earthquakes, and compared against observed liquefaction surface manifestations and ground motions characteristics at each SMS. In general, the characteristics of the recorded ground motions at each site correlate well with the triggering analyses. However, at sites that likely liquefied at depth (as indicated by triggering analyses and/or inferred from the characteristics of the recorded surface acceleration time series), the presence of a non‐liquefiable crust layer at many of the SMS locations prevented the manifestation of any surface effects.
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.
Recent experiences from the Darfield and Canterbury, New Zealand earthquakes have shown that the soft soil condition of saturated liquefiable sand has a profound effect on seismic response of buildings, bridges and other lifeline infrastructure. For detailed evaluation of seismic response three dimensional integrated analysis comprising structure, foundation and soil is required; such an integrated analysis is referred to as Soil Foundation Structure Interaction (SFSI) in literatures. SFSI is a three-dimensional problem because of three primary reasons: first, foundation systems are three-dimensional in form and geometry; second, ground motions are three-dimensional, producing complex multiaxial stresses in soils, foundations and structure; and third, soils in particular are sensitive to complex stress because of heterogeneity of soils leading to a highly anisotropic constitutive behaviour. In literatures the majority of seismic response analyses are limited to plane strain configuration because of lack of adequate constitutive models both for soils and structures, and computational limitation. Such two-dimensional analyses do not represent a complete view of the problem for the three reasons noted above. In this context, the present research aims to develop a three-dimensional mathematical formulation of an existing plane-strain elasto-plastic constitutive model of sand developed by Cubrinovski and Ishihara (1998b). This model has been specially formulated to simulate liquefaction behaviour of sand under ground motion induced earthquake loading, and has been well-validated and widely implemented in verifcation of shake table and centrifuge tests, as well as conventional ground response analysis and evaluation of case histories. The approach adopted herein is based entirely on the mathematical theory of plasticity and utilises some unique features of the bounding surface plasticity formalised by Dafalias (1986). The principal constitutive parameters, equations, assumptions and empiricism of the existing plane-strain model are adopted in their exact form in the three-dimensional version. Therefore, the original two-dimensional model can be considered as a true subset of the three-dimensional form; the original model can be retrieved when the tensorial quantities of the three dimensional version are reduced to that of the plane-strain configuration. Anisotropic Drucker-Prager type failure surface has been adopted for the three-dimensional version to accommodate triaxial stress path. Accordingly, a new mixed hardening rule based on Mroz’s approach of homogeneous surfaces (Mroz, 1967) has been introduced for the virgin loading surface. The three-dimensional version is validated against experimental data for cyclic torsional and triaxial stress paths.
Recent surface-rupturing earthquakes in New Zealand have highlighted significant exposure and vulnerability of the road network to fault displacement. Understanding fault displacement hazard and its impact on roads is crucial for mitigating risks and enhancing resilience. There is a need for regional-scale assessments of fault displacement to identify vulnerable areas within the road network for the purposes of planning and prioritising site-specific investigations. This thesis employs updated analysis of data from three historical surface-rupturing earthquakes (Edgecumbe 1987, Darfield 2010, and Kaikoūra 2016) to develop an empirical model that addresses the gap in regional fault displacement hazard analysis. The findings contribute to understanding of • How to use seismic hazard model inputs for regional fault displacement hazard analysis • How faulting type and sediment cover affects the magnitude and spatial distribution of fault displacement • How the distribution of displacement and regional fault displacement hazard is impacted by secondary faulting • The inherent uncertainties and limitations associated with employing an empirical approach at a regional scale • Which sections of New Zealand’s roading network are most susceptible to fault displacement hazard and warrant site-specific investigations • Which regions should prioritise updating emergency management plans to account for post-event disruptions to roading. I used displacement data from the aforementioned historical ruptures to generate displacement versus distance-to-fault curves for various displacement components, fault types, and geological characteristics. Using those relationships and established relationships for along-strike displacement, displacement contours were generated surrounding active faults within the NZ Community Fault Model. Next, I calculated a new measure of 1D strain along roads as well as relative hazard, which integrated 1D strain and normalised slip rate data. Summing these values at the regional level identified areas of heightened relative hazard across New Zealand, and permits an assessment of the susceptibility of road networks using geomorphon land classes as proxies for vulnerability. The results reveal that fault-parallel displacements tend to localise near the fault plane, while vertical and fault-perpendicular displacements sustain over extended distances. Notably, no significant disparities were observed in off-fault displacement between the hanging wall and footwall sides of the fault, or among different surface geology types, potentially attributed to dataset biases. The presence of secondary faulting in the dataset contributes to increased levels of tectonic displacement farther from the fault, highlighting its significance in hazard assessments. Furthermore, fault displacement contours delineate broader zones around dip-slip faults compared to strike-slip faults, with correlations identified between fault length and displacement width. Road ‘strain’ values are higher around dip-slip faults, with notable examples observed in the Westland and Buller Districts. As expected, relative hazard analysis revealed elevated values along faults with high slip rates, notably along the Alpine Fault. A regional-scale analysis of hazard and exposure reveals heightened relative hazard in specific regions, including Wellington, Southern Hawke’s Bay, Central Bay of Plenty, Central West Coast, inland Canterbury, and the Wairau Valley of Marlborough. Notably, the Central West Coast exhibits the highest summed relative hazard value, attributed to the fast-slipping Alpine Fault. The South Island generally experiences greater relative hazard due to larger and faster-slipping faults compared to the North Island, despite having fewer roads. Central regions of New Zealand face heightened risk compared to Southern or Northern regions. Critical road links intersecting high-slipping faults, such as State Highways 6, 73, 1, and 2, necessitate prioritisation for site-specific assessments, emergency management planning and targeted mitigation strategies. Roads intersecting with the Alpine Fault are prone to large parallel displacements, requiring post-quake repair efforts. Mitigation strategies include future road avoidance of nearby faults, modification of road fill and surface material, and acknowledgement of inherent risk, leading to prioritised repair efforts of critical roads post-quake. Implementing these strategies enhances emergency response efforts by improving accessibility to isolated regions following a major surface-rupturing event, facilitating faster supply delivery and evacuation assistance. This thesis contributes to the advancement of understanding fault displacement hazard by introducing a novel regional, empirical approach. The methods and findings highlight the importance of further developing such analyses and extending them to other critical infrastructure types exposed to fault displacement hazard in New Zealand. Enhancing our comprehension of the risks associated with fault displacement hazard offers valuable insights into various mitigation strategies for roading infrastructure and informs emergency response planning, thereby enhancing both national and global infrastructure resilience against geological hazards.
