A major hazard accompanying earthquake shaking in areas of steep topography is the detachment of rocks from bedrock outcrops that subsequently slide, roll, or bounce downslope (i.e. rockfalls). The 2010-2011 Canterbury earthquake sequence caused recurrent and severe rockfall in parts of southern Christchurch. Coseismic rockfall caused five fatalities and significant infrastructural damage during the 2011 Mw 6.2 Christchurch earthquake. Here we examine a rockfall site in southern Christchurch in detail using geomorphic mapping, lidar analysis, geochronology (cosmogenic 3He dating, radiocarbon dating, optically stimulated luminescence (OSL) from quartz, infrared stimulated luminescence from K-feldspar), numerical modeling of rockfall boulder trajectories, and ground motion prediction equations (GMPEs). Rocks fell from the source cliff only in earthquakes with interpolated peak ground velocities exceeding ~10 cm/s; hundreds of smaller earthquakes did not produce rockfall. On the basis of empirical observations, GMPEs and age chronologies we attribute paleo-rockfalls to strong shaking in prehistoric earthquakes. We conclude that earthquake shaking of comparable intensity to the strongest contemporary earthquakes in Christchurch last occurred at this site approximately 5000 to 7000 years ago, and that in some settings, rockfall deposits provide useful proxies for past strong ground motions.
Unrestrained unreinforced clay brick masonry (URM) parapets are found atop a large number of vintage URM buildings. Parapets are typically non-structural cantilevered wall elements that form a fire barrier and in most cases form decorative and ornamental features of vintage URM buildings. Parapets are considered to be one of the most vulnerable elements that are prone to out-of-plane collapse when subjected to earthquake induced shaking. An in-depth analysis of the damage database collected following the 2010/2011 Canterbury earthquakes was performed to obtain information about the distribution, characteristics and observed performance of both the as-built and retrofitted parapets in the Christchurch region. Results, statistical interpretation and implications are presented herein. http://www.aees.org.au/downloads/conference-papers/2015-2/
Recent earthquakes have shown that liquefaction and associated ground deformations are major geotechnical hazards to civil engineering infrastructures, such as pipelines. In particular, sewer pipes have been damaged in many areas in Christchurch as a result of liquefaction-induced lateral spreading near waterways and ground oscillation induced by seismic shaking. In this paper, the addition of a flexible AM liner as a potential countermeasure to increase sewer pipe capacity was investigated. Physical testing through 4-point loading test was undertaken to characterise material properties and the response of both unlined pipe and its lined counterpart. Next, numerical models were created using SAP2000 and ABAQUS to analyse buried pipeline response to transverse permanent ground displacement and to quantify, over a range of pipe segment lengths and soil parameters, the effectiveness of the AM liner in increasing displacement capacity. The numerical results suggest that the addition of the AM liner increases the deformation capacity of the unlined sewer pipe by as much as 50 times. The results confirmed that AM liner is an effective countermeasure for sewer pipes in liquefied ground not only in terms of increased deformation capacity but also the fact that AM-Liner can prevent influx of sand and water through broken pipes, making sewer pipes with liner remaining serviceable even under severe liquefaction condition.
Recurrent liquefaction in Christchurch during the 2010-2011 Canterbury earthquake sequence created a wealth of shallow subsurface intrusions with geometries and orientations governed by (1) strong ground motion severity and duration, and (2) intrinsic site characteristics including liquefaction susceptibility, lateral spreading severity, geomorphic setting, host sediment heterogeneity, and anthropogenic soil modifications. We present a suite of case studies that demonstrate how each of these characteristics influenced the geologic expressions of contemporary liquefaction in the shallow subsurface. We compare contemporary features with paleo-features to show how geologic investigations of recurrent liquefaction can provide novel insights into the shaking characteristics of modern and paleo-earthquakes, the influence of geomorphology on liquefaction vulnerability, and the possible controls of anthropogenic activity on the geologic record. We conclude that (a) sites of paleo-liquefaction in the last 1000-2000 years corresponded with most severe liquefaction during the Canterbury earthquake sequence, (b) less vulnerable sites that only liquefied in the strongest and most proximal contemporary earthquakes are unlikely to have liquefied in the last 1000-2000 years or more, (c) proximal strong earthquakes with large vertical accelerations favoured sill formation at some locations, (d) contemporary liquefaction was more severe than paleoliquefaction at all study sites, and (e) stratigraphic records of successive dike formation were more complete at sites with severe lateral spreading, (f) anthropogenic fill suppressed surface liquefaction features and altered subsurface liquefaction architecture.
Severe liquefaction was repeatedly observed during the 2010 - 2011 C hristchurch earthquake s , particularly affecting deposits of fine sands and silty sands of recent fluvial or estuarine origin. The effects of liquefaction included major sliding of soil tow ard water bodies ( i.e. lateral spreading ) rang ing from centimetres to several metres. In this paper, a series of undrained cyclic torsional shear tests were conducted to evaluate the liquefaction and extremely large deformation properties of Christchurch b oiled sand . In these tests, the simple shear conditions were reproduced in order to apply realistic stress conditions that soil s experience in the field during horizontal seismic shaking. Several hollow cylindrical medium dense specimens ( D r = 50%) were pr epared by pluviation method, isotropically consolidated at an effective stress of 100 kPa and then cyclically sheared under undrained conditions up to 10 0% double amplitude shear strain (γ DA ) . The cyclic strength at different levels of γ DA of 7.5%, 15%, 3 0 % and 6 0%, development of extremely large post - liquefaction deformation and shear strain locali s ation properties were assessed from the analysis of the effective stress paths and stress - strain responses . To reveal possible distinctiveness, the cyclic undra ined behaviour of CHCH boiled sand was compared with that of Toyoura sand previously examined under similar testing conditions
The Global Earthquake Model’s (GEM) Earthquake Consequences Database (GEMECD) aims to develop, for the first time, a standardised framework for collecting and collating geocoded consequence data induced by primary and secondary seismic hazards to different types of buildings, critical facilities, infrastructure and population, and relate this data to estimated ground motion intensity via the USGS ShakeMap Atlas. New Zealand is a partner of the GEMECD consortium and to-date has contributed with 7 events to the database, of which 4 are localised in the South Pacific area (Newcastle 1989; Luzon 1990; South of Java 2006 and Samoa Islands 2009) and 3 are NZ-specific events (Edgecumbe 1987; Darfield 2010 and Christchurch 2011). This contribution to GEMECD represented a unique opportunity for collating, comparing and reviewing existing damage datasets and harmonising them into a common, openly accessible and standardised database, from where the seismic performance of New Zealand buildings can be comparatively assessed. This paper firstly provides an overview of the GEMECD database structure, including taxonomies and guidelines to collect and report on earthquake-induced consequence data. Secondly, the paper presents a summary of the studies implemented for the 7 events, with particular focus on the Darfield (2010) and Christchurch (2011) earthquakes. Finally, examples of specific outcomes and potentials for NZ from using and processing GEMECD are presented, including: 1) the rationale for adopting the GEM taxonomy in NZ and any need for introducing NZ-specific attributes; 2) a complete overview of the building typological distribution in the Christchurch CBD prior to the Canterbury earthquakes and 3) some initial correlations between the level and extent of earthquake-induced physical damage to buildings, building safety/accessibility issues and the induced human casualties.
As part of a seismic retrofit scheme, surface bonded glass fiber-reinforced polymer (GFRP) fabric was applied to two unreinforced masonry (URM) buildings located in Christchurch, New Zealand. The unreinforced stone masonry of Christchurch Girls’ High School (GHS) and the unreinforced clay brick masonry Shirley Community Centre were retrofitted using surface bonded GFRP in 2007 and 2009, respectively. Much of the knowledge on the seismic performance of GFRP retrofitted URM was previously assimilated from laboratory-based experimental studies with controlled environments and loading schemes. The 2010/2011 Canterbury earthquake sequence provided a rare opportunity to evaluate the GFRP retrofit applied to two vintage URM buildings and to document its performance when subjected to actual design-level earthquake-induced shaking. Both GFRP retrofits were found to be successful in preserving architectural features within the buildings as well as maintaining the structural integrity of the URM walls. Successful seismic performance was based on comparisons made between the GFRP retrofitted GHS building and the adjacent nonretrofitted Boys’ High School building, as well as on a comparison between the GFRP retrofitted and nonretrofitted walls of the Shirley Community Centre building. Based on detailed postearthquake observations and investigations, the GFRP retrofitted URM walls in the subject buildings exhibited negligible to minor levels of damage without delamination, whereas significant damage was observed in comparable nonretrofitted URM walls. AM - Accepted Manuscript
The Canterbury region experienced widespread damage due to liquefaction induced by seismic shaking during the 4 September 2010 earthquake and the large aftershocks that followed, notably those that occurred on 22 February, 13 June and 23 December 2011. Following the 2010 earthquake, the Earthquake Commission directed a thorough investigation of the ground profile in Christchurch, and to date, more than 7500 cone penetration tests (CPT) have been performed in the region. This paper presents the results of analyses which use a subset of the geotechnical database to evaluate the liquefaction process as well as the re-liquefaction that occurred following some of the major events in Christchurch. First, the applicability of existing CPT-based methods for evaluating liquefaction potential of Christchurch soils was investigated using three methods currently available. Next, the results of liquefaction potential evaluation were compared with the severity of observed damage, categorised in terms of the land damage grade developed from Tonkin & Taylor property inspections as well as from observed severity of liquefaction from aerial photography. For this purpose, the Liquefaction Potential Index (LPI) was used to represent the damage potential at each site. In addition, a comparison of the CPT-based strength profiles obtained before each of the major aftershocks was performed. The results suggest that the analysis of spatial and temporal variations of strength profiles gives a clear indication of the resulting liquefaction and re-liquefaction observed in Christchurch. The comparison of a limited number of CPT strength profiles before and after the earthquakes seems to indicate that no noticeable strengthening has occurred in Christchurch, making the area vulnerable to liquefaction induced land damage in future large-scale earthquakes.
Following the Mw 6.2 Christchurch Earthquake on 22 February 2011, extensive ground cracking in loessial soils was reported in some areas of the Port Hills, southeast of central Christchurch. This study was undertaken to investigate the mechanisms of earthquake-induced ground damage on the eastern side of the Hillsborough Valley. A zone of extensional cracking up to 40m wide and 600m long was identified along the eastern foot-slope, accompanied by compression features and spring formation at the toe of the slope. An engineering geological and geomorphological model was developed for the eastern Hillsborough Valley that incorporates geotechnical investigation data sourced from the Canterbury Geotechnical Database (CGD), the findings of trenching and seismic refraction surveying carried out for this research, and interpretation of historical aerial photographs. The thickness and extent of a buried peat swamp at the base of the slope was mapped, and found to coincide with significant compression features. Ground cracking was found to have occurred entirely within loess-colluvium and to follow the apices of pre-1920s tunnel-gully fan debris at the southern end of the valley. The ground-cracking on the eastern side of the Hillsborough Valley is interpreted to have formed through tensile failure of the loess-colluvium. Testing was carried out to determine the tensile strength of Port Hills loess colluvium as a function of water content and density, in order to better understand the occurrence and distribution of the observed ground cracking. A comprehensive review of the soil tensile strength testing literature was undertaken, from which a test methodology was developed. Results show remoulded loess-colluvium to possess tensile strength of 7 - 28 kPa across the range of tested moisture contents (10-15%) and dry densities (1650-1900kg/m3). A positive linear relationship was observed between tensile strength and dry density, and a negative linear relationship between moisture content and tensile strength. The observed ground damage and available geotechnical information (inclinometer and piezometer records provided by the Earthquake Commission) were together used to interpret the mechanism(s) of slope movement that occurred in the eastern Hillsborough Valley. The observed ground damage is characteristic of translational movement, but without the development of lateral release scarps, or a basal sliding surface - which was not located during drilling. It is hypothesised that shear displacement has been accommodated by multiple slip surfaces of limited extent within the upper 10m of the slope. Movement has likely occurred within near-saturated colluvial units that have lost strength during earthquake shaking. The eastern Hillsborough Valley is considered to be an ‘incipient translational slide’, as both the patterns of damage and shearing are consistent with the early stages of such slide development. Sliding block analysis was utilised to understand how the eastern Hillsborough Valley may perform in a future large magnitude earthquake. Known cumulative displacements of ~0.3m for eastern Hillsborough Valley during the 2010-2011 Canterbury Earthquake Sequence were compared with modelled slope displacements to back-analyse a lower-bound yield acceleration of 0.2 - 0.25g. Synthetic broadband modelling for future Alpine and Hope Fault earthquakes indicates PGAs of approximately 0.08g for soil sites in the Christchurch area, as such, slope movement is unlikely to be reactivated by an Alpine Fault or Hope Fault earthquake. This does not take into account the possible role of strength loss due to excess pore pressure that may occur during these future events.
As a consequence of the 2010 – 2011 Canterbury earthquake sequence, Christchurch experienced widespread liquefaction, vertical settlement and lateral spreading. These geological processes caused extensive damage to both housing and infrastructure, and increased the need for geotechnical investigation substantially. Cone Penetration Testing (CPT) has become the most common method for liquefaction assessment in Christchurch, and issues have been identified with the soil behaviour type, liquefaction potential and vertical settlement estimates, particularly in the north-western suburbs of Christchurch where soils consist mostly of silts, clayey silts and silty clays. The CPT soil behaviour type often appears to over-estimate the fines content within a soil, while the liquefaction potential and vertical settlement are often calculated higher than those measured after the Canterbury earthquake sequence. To investigate these issues, laboratory work was carried out on three adjacent CPT/borehole pairs from the Groynes Park subdivision in northern Christchurch. Boreholes were logged according to NZGS standards, separated into stratigraphic layers, and laboratory tests were conducted on representative samples. Comparison of these results with the CPT soil behaviour types provided valuable information, where 62% of soils on average were specified by the CPT at the Groynes Park subdivision as finer than what was actually present, 20% of soils on average were specified as coarser than what was actually present, and only 18% of soils on average were correctly classified by the CPT. Hence the CPT soil behaviour type is not accurately describing the stratigraphic profile at the Groynes Park subdivision, and it is understood that this is also the case in much of northwest Christchurch where similar soils are found. The computer software CLiq, by GeoLogismiki, uses assessment parameter constants which are able to be adjusted with each CPT file, in an attempt to make each more accurate. These parameter changes can in some cases substantially alter the results for liquefaction analysis. The sensitivity of the overall assessment method, raising and lowering the water table, lowering the soil behaviour type index, Ic, liquefaction cutoff value, the layer detection option, and the weighting factor option, were analysed by comparison with a set of ‘base settings’. The investigation confirmed that liquefaction analysis results can be very sensitive to the parameters selected, and demonstrated the dependency of the soil behaviour type on the soil behaviour type index, as the tested assessment parameters made very little to no changes to the soil behaviour type plots. The soil behaviour type index, Ic, developed by Robertson and Wride (1998) has been used to define a soil’s behaviour type, which is defined according to a set of numerical boundaries. In addition to this, the liquefaction cutoff point is defined as Ic > 2.6, whereby it is assumed that any soils with an Ic value above this will not liquefy due to clay-like tendencies (Robertson and Wride, 1998). The method has been identified in this thesis as being potentially unsuitable for some areas of Christchurch as it was developed for mostly sandy soils. An alternative methodology involving adjustment of the Robertson and Wride (1998) soil behaviour type boundaries is proposed as follows: Ic < 1.31 – Gravelly sand to dense sand 1.31 < Ic < 1.90 – Sands: clean sand to silty sand 1.90 < Ic < 2.50 – Sand mixtures: silty sand to sandy silt 2.50 < Ic < 3.20 – Silt mixtures: clayey silt to silty clay 3.20 < Ic < 3.60 – Clays: silty clay to clay Ic > 3.60 – Organics soils: peats. When the soil behaviour type boundary changes were applied to 15 test sites throughout Christchurch, 67% showed an improved change of soil behaviour type, while the remaining 33% remained unchanged, because they consisted almost entirely of sand. Within these boundary changes, the liquefaction cutoff point was moved from Ic > 2.6 to Ic > 2.5 and altered the liquefaction potential and vertical settlement to more realistic ii values. This confirmed that the overall soil behaviour type boundary changes appear to solve both the soil behaviour type issues and reduce the overestimation of liquefaction potential and vertical settlement. This thesis acts as a starting point towards researching the issues discussed. In particular, future work which would be useful includes investigation of the CLiq assessment parameter adjustments, and those which would be most suitable for use in clay-rich soils such as those in Christchurch. In particular consideration of how the water table can be better assessed when perched layers of water exist, with the limitation that only one elevation can be entered into CLiq. Additionally, a useful investigation would be a comparison of the known liquefaction and settlements from the Canterbury earthquake sequence with the liquefaction and settlement potentials calculated in CLiq for equivalent shaking conditions. This would enable the difference between the two to be accurately defined, and a suitable adjustment applied. Finally, inconsistencies between the Laser-Sizer and Hydrometer should be investigated, as the Laser-Sizer under-estimated the fines content by up to one third of the Hydrometer values.
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.
