On 22 February 2011，a magnitude Mw 6.3 earthquake occurred with an epicenter located near Lyttelton at about 10km from Christchurch in Canterbury region on the South Island of New Zealand (Figure 1). Since this earthquake occurred in the midst of the aftershock activity which had continued since the 4 September 2010 Darfield Earthquake occurrence, it was considered to be an aftershock of the initial earthquake. Because of the short distance to the city and the shallower depth of the epicenter, this earthquake caused more significant damage to pipelines, traffic facilities, residential houses/properties and multi-story buildings in the central business district than the September 2010 Darfield Earthquake in spite of its smaller earthquake magnitude. Unfortunately, this earthquake resulted in significant number of casualties due to the collapse of multi-story buildings and unreinforced masonry structures in the city center of Christchurch. As of 4 April, 172 casualties were reported and the final death toll is expected to be 181. While it is extremely regrettable that Christchurch suffered a terrible number of victims, civil and geotechnical engineers have this hard-to-find opportunity to learn the response of real ground from two gigantic earthquakes which occurred in less than six months from each other. From geotechnical engineering point of view, it is interesting to discuss the widespread liquefaction in natural sediments, repeated liquefaction within short period and further damage to earth structures which have been damaged in the previous earthquake. Following the earthquake, an intensive geotechnical reconnaissance was conducted to capture evidence and perishable data from this event. The team included the following members: Misko Cubrinovski (University of Canterbury, NZ, Team Leader), Susumu Yasuda (Tokyo Denki University, Japan, JGS Team Leader), Rolando Orense (University of Auckland, NZ), Kohji Tokimatsu (Tokyo Institute of Technology, Japan), Ryosuke Uzuoka (Tokushima University, Japan), Takashi Kiyota (University of Tokyo, Japan), Yasuyo Hosono (Toyohashi University of Technology, Japan) and Suguru Yamada (University of Tokyo, Japan).
The region in and around Christchurch, encompassing Christchurch city and the Selwyn and Waimakariri districts, contains more than 800 road, rail, and pedestrian bridges. Most of these bridges are reinforced concrete, symmetric, and have small to moderate spans (15–25 m). The 22 February 2011 moment magnitude (Mw) 6.2 Christchurch earthquake induced high levels of localized ground shaking (Bradley and Cubrinovski 2011, page 853 of this issue; Guidotti et al. 2011, page 767 of this issue; Smyrou et al. 2011, page 882 of this issue), with damage to bridges mainly confined to the central and eastern parts of Christchurch. Liquefaction was evident over much of this part of the city, with lateral spreading affecting bridges spanning both the Avon and Heathcote rivers.
The paper discusses modelling of cyclic stress-strain behaviour of soil, in particular a simple model that can produce a desired stiffness and hysteretic damping for a given strain level as observed in laboratory testing is formulated. The unloading-reloading relationship is developed for total stress seismic site response analysis with appropriate damping at large strain. The constitutive model employs a hyperbolic equation as the backbone curve, and uses a modification of the extended Masing unloading-reloading relationship leading to correct measured modulus reduction and damping curves simultaneously. A quasi-static cyclic loading of increasing amplitude is used to demonstrate the model’s performance and its capability to allow improved modelling of the magnitude of energy dissipation based on an experimental program on native sandy soils from Christchurch, New Zealand.
Geomorphic, structural and chronological data are used to establish the late Quaternary paleoseismicity of the active dextral-oblique Northern Esk Fault in North Canterbury, New Zealand. Detailed field mapping of the preserved c. 35 km of surface traces between the Hurunui River and Ashley Head reveals variations in strike ranging from 005° to 057°. Along with kinematic data collected from fault plane striae and offset geomorphic markers along the length of the fault these variations are used to distinguish six structural subsections of the main trace, four dextral-reverse and two dextral-normal. Displacements of geomorphic markers such as minor streams and ridges are measured using differential GPS and rangefinder equipment to reveal lateral offsets ranging from 3.4 to 23.7 m and vertical offsets ranging from < 1 to 13.5 m. Characteristic single event displacements of c. 5 m and c. 2 m have been calculated for strike-slip and reverse sections respectively. The use of fault scaling relationships reveals an anomalously high displacement to surface rupture length ratio when compared to global data sets. Fault scaling relationships based on width limited ruptures and magnitude probabilities from point measurements of displacement imply earthquake magnitudes of Mw 7.0 to 7.5. Optically Stimulated Luminescence (OSL) ages from displaced Holocene alluvial terraces at the northern extent of the active trace along with OSL and radiocarbon samples of the central sections constrain the timing of the last two surface rupturing events (11.15 ±1.65 and 3.5 ± 2.8 ka) and suggest a recurrence interval of c. 5612 ± 445 years and late Quaternary reverse and dextral slip rates of c. 0.31 mm/yr and 0.82 mm/yr respectively. The results of this study show that the Northern Esk Fault accommodates an important component of the c. 0.7 – 2 mm/yr of unresolved strain across the plate boundary within the North Canterbury region and affirm the Esk Fault as a source of potentially damaging ground shaking in the Canterbury region.
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.
Following the magnitude 6.3 aftershock in Christchurch, New Zealand, on 22 February 2011, a number of researchers were sent to Christchurch as part of the New Zealand Natural Hazard Research Platform funded “Project Masonry” Recovery Project. Their goal was to document and interpret the damage to the masonry buildings and churches in the region. Approximately 650 unreinforced and retrofitted clay brick masonry buildings in the Christchurch area were surveyed for commonly occurring failure patterns and collapse mechanisms. The entire building stock of Christchurch, and in particular the unreinforced masonry building stock, is similar to that in the rest of New Zealand, Australia, and abroad, so the observations made here are relevant for the entire world.
The University of Canterbury Dept. of Chemistry has weathered the Canterbury Earthquake of September 4, 2010 very well due to a combination of good luck, good planning and dedicated effort. We owe a great deal to university Emergency Response Team and Facilities Management Personnel. The overall emergency preparedness of the university was tested to a degree far beyond anything else in its history and shown to be well up to scratch. A strong cooperative relationship between the pan-campus controlling body and the departmental response teams greatly facilitated our efforts. Information and assistance was provided promptly, as and when we needed it without unnecessary bureaucratic overheads. At the departmental level we are indebted to the technical staff who implemented the invaluable pre-quake mitigation measures and carried the majority of the post-quake clean-up workload. These people put aside their personal concerns and anxieties at a time when magnitude-5 aftershocks were still a regular occurrence.
The devastating magnitude M6.3 earthquake, that struck the city of Christchurch at 12:51pm on Tuesday 22 February 2011, caused widespread damage to the lifeline systems. Following the event, the Natural Hazard Research Platform (NHRP) of New Zealand funded a short-term project “Recovery of Lifelines” aiming to: 1) coordinate the provision of information to meet lifeline short-term needs; and to 2) facilitate the accessibility to lifelines of best practice engineering details, along with hazards and vulnerability information already available from the local and international scientific community. This paper aims to briefly summarise the management of the recovery process for the most affected lifelines systems, including the electric system, the road, gas, and the water and wastewater networks. Further than this, the paper intends to discuss successes and issues encountered by the “Recovery of Lifelines” NHRP project in supporting lifelines utilities.
This dissertation addresses a diverse range of topics in the physics-based broadband ground motion simulation, with a focus on New Zealand applications. In particular the following topics are addressed: the methodology and computational implementation of a New Zealand Velocity Model for broadband ground motion simulation; generalised parametric functions and spatial correlations for seismic velocities in the Canterbury, New Zealand region from surface-wave-based site characterisation; and ground motion simulations of Hope Fault earthquakes. The paragraphs below outline each contribution in more detail. A necessary component in physics-based ground motion simulation is a 3D model which details the seismic velocities in the region of interest. Here a velocity model construction methodology, its computational implementation, and application in the construction of a New Zealand velocity model for use in physics-based broadband ground motion simulation are presented. The methodology utilises multiple datasets spanning different length scales, which is enabled via the use of modular sub-regions, geologic surfaces, and parametric representations of crustal velocity. A number of efficiency-related workflows to decrease the overall computational construction time are employed, while maintaining the flexibility and extensibility to incorporate additional datasets and re- fined velocity parameterizations as they become available. The model comprises explicit representations of the Canterbury, Wellington, Nelson-Tasman, Kaikoura, Marlborough, Waiau, Hanmer and Cheviot sedimentary basins embedded within a regional travel-time tomography-based velocity model for the shallow crust and provides the means to conduct ground motion simulations throughout New Zealand for the first time. Recently developed deep shear-wave velocity profiles in Canterbury enabled models that better characterise the velocity structure within geologic layers of the Canterbury sedimentary basin to be developed. Here the development of depth- and Vs30-dependent para-metric velocity and spatial correlation models to characterise shear-wave velocities within the geologic layers of the Canterbury sedimentary basin are presented. The models utilise data from 22 shear-wave velocity profiles of up to 2.5km depth (derived from surface wave analysis) juxtaposed with models which detail the three-dimensional structure of the geologic formations in the Canterbury sedimentary basin. Parametric velocity equations are presented for Fine Grained Sediments, Gravels, and Tertiary layer groupings. Spatial correlations were developed and applied to generate three-dimensional stochastic velocity perturbations. Collectively, these models enable seismic velocities to be realistically represented for applications such as 3D ground motion and site response simulations. Lastly the New Zealand velocity model is applied to simulate ground motions for a Mw7.51 rupture of the Hope Fault using a physics-based simulation methodology and a 3D crustal velocity model of New Zealand. The simulation methodology was validated for use in the region through comparison with observations for a suite of historic small magnitude earthquakes located proximal to the Hope Fault. Simulations are compared with conventionally utilised empirical ground motion models, with simulated peak ground velocities being notably higher in regions with modelled sedimentary basins. A sensitivity analysis was undertaken where the source characteristics of magnitude, stress parameter, hypocentre location and kinematic slip distribution were varied and an analysis of their effect on ground motion intensities is presented. It was found that the magnitude and stress parameter strongly influenced long and short period ground motion amplitudes, respectively. Ground motion intensities for the Hope Fault scenario are compared with the 2016 Kaikoura Mw7.8 earthquake, it was found that the Kaikoura earthquake produced stronger motions along the eastern South Island, while the Hope Fault scenario resulted in stronger motions immediately West of the near-fault region. The simulated ground motions for this scenario complement prior empirically-based estimates and are informative for mitigation and emergency planning purposes.
Oblique convergence of the Pacific and Australian Plates is accommodated in the northern South Island by the Marlborough Fault System. The Hope Fault is the southern of four major dextral strike-slip faults of this system. Hanmer Basin is a probable segment boundary between the Hope River and Conway segments of the Hope Fault. The Conway segment is transpressional and shows increasing structural complexity near the segment boundary at Hanmer Basin, with multiple Late Quaternary traces, and fault-parallel folding in response to across-fault shortening. Between Hossack Station and Hanmer Basin a crush zone in excess of one kilometre wide is exposed in incised streams and rivers. The crush zone has an asymmetrical geometry about the active trace of the Hope Fault, being only 100-300 metres wide south of the fault, and more than 500 metres wide north of the fault. The most intense deformation of Torlesse bedrock occurs at the south side of the fault zone, indicating that strain is accommodated against the fault footwall. North of the fault deformation is less intense, but occurs over a wider area. The wide fault zone at Hossack Station may reflect divergence of the Hanmer Fault, a major splay of the Hope Fault. At Hossack Station, the Hope Fault has accommodated at least 260 metres of dextral displacement during the Holocene. Dating of abandoned stream channels, offset by the Hope Fault, indicated a Late Holocene dextral slip-rate of 18±8 mm-¹ for the west end of the Conway segment. Using empirical formulae and inferred fault parameters, the expected magnitude of an earthquake generated by the Conway segment is M6.9 to M7.4; for an exceedence probability of 10%, the magnitude is M7.7 to M7.9. Effects associated with coseismic rupture of the Conway segment include shaking of up to MMIX along the ruptured fault and at Hanmer Basin. Uplift at the east end of Hanmer Basin, in conjunction with subsidence at the southwest margin of the basin, is resulting in the development of onlapping stratigraphy. Seismic reflection profiles support this theory. Possible along-fault migration of the basin is inferred to be a consequence of non-parallelism of the master faults.
Unreinforced masonry (URM) buildings have repeatedly been shown to perform poorly in large magnitude earthquakes, with both New Zealand and Australia having a history of past earthquakes that have resulted in fatalities due to collapsed URM buildings. A comparison is presented here of the URM building stock and the seismic vulnerability of Christchurch and Adelaide in order to demonstrate the relevance to Australian cities of observations in Christchurch resulting from the 2010/2011 Canterbury earthquake swarm. It is shown that the materials, architecture and hence earthquake strength of URM buildings in both countries is comparable and that Adelaide and other cities of Australia have seismic vulnerability sufficient to cause major damage to their URM buildings should a design level earthquake occur. Such an earthquake is expected to cause major building damage, and fatalities should be expected.
This paper examines the consistency of seismicity and ground motion models, used for seismic hazard analysis in New Zealand, with the observations in the Canterbury earthquakes. An overview is first given of seismicity and ground motion modelling as inputs of probabilistic seismic hazard analysis, whose results form the basis for elastic response spectra in NZS1170.5:2004. The magnitude of earthquakes in the Canterbury earthquake sequence are adequately allowed for in the current NZ seismicity model, however the consideration of ‘background’ earthquakes as point sources at a minimum depth of 10km results in up to a 60% underestimation of the ground motions that such events produce. The ground motion model used in conventional NZ seismic hazard analysis is shown to provide biased predictions of response spectra (over-prediction near T=0.2s , and under-predictions at moderate-to-large vibration periods). Improved ground motion prediction can be achieved using more recent NZ-specific models.
On 4 September 2010, a magnitude Mw 7.1 earthquake struck the Canterbury region on the South Island of New Zealand. The epicentre of the earthquake was located in the Darfield area about 40 km west of the city of Christchurch. Extensive damage was inflicted to lifelines and residential houses due to widespread liquefaction and lateral spreading in areas close to major streams, rivers and wetlands throughout Christchurch and Kaiapoi. Unreinforced masonry buildings also suffered extensive damage throughout the region. Despite the severe damage to infrastructure and residential houses, fortunately, no deaths occurred and only two injuries were reported in this earthquake. From an engineering viewpoint, one may argue that the most significant aspects of the 2010 Darfield Earthquake were geotechnical in nature, with liquefaction and lateral spreading being the principal culprits for the inflicted damage. Following the earthquake, an intensive geotechnical reconnaissance was conducted to capture evidence and perishable data from this event. This paper summarizes the observations and preliminary findings from this early reconnaissance work.
The UC CEISMIC Canterbury Earthquakes Digital Archive was built following the devastating earthquakes that hit the Canterbury region in the South Island of New Zealand from 2010 – 2012. 185 people were killed in the 6.3 magnitude earthquake of February 22nd 2011, thousands of homes and businesses were destroyed, and the local community endured over 10,000 aftershocks. The program aims to document and protect the social, cultural, and intellectual legacy of the Canterbury community for the purposes of memorialization and enabling research. The nationally federated archive currently stores 75,000 items, ranging from audio and video interviews to images and official reports. Tens of thousands more items await ingestion. Significant lessons have been learned about data integration in post-disaster contexts, including but not limited to technical architecture, governance, ingestion process, and human ethics. The archive represents a model for future resilience-oriented data integration and preservation products.
On February 22, 2011, Christchurch-based journalists were jolted out of their normal work routine by a large 6.3 magnitude earthquake that killed 185 people, wrecked the city and forced reporters to reappraise their journalism. This study considers how the earthquake affected journalists’ relationship to the community, their use of sources and news selection. A theory of collective trauma is used to explain the changes that journalists made to their reporting practice. Specifically, Christchurch journalists had a greater identification and attachment to their audience post-earthquake. Journalists viewed themselves as part of the earthquake story, which prompted them to view sources differently, use those sources differently and see advocacy as a keystone of their news work after the disaster. This study adds to a growing scholarship about journalists and trauma, but focuses on what the event meant for local reporters’ choice of sources and news selection rather than measuring rates of psychological distress.
This report provides an understanding of the nature of Canterbury subcontracting businesses operating in the space of earthquake reconstruction in Christchurch. It offers an in-depth look at the factors that influence the development of their capacity and capability to withstand the impact of volatile economic cycles, including the 2008 global financial crisis and the subsequent 2010/11 Canterbury earthquakes. There have been significant changes to the business models of the 13 subcontracting businesses studied since the earthquakes. These changes can be seen in the ways the case study subcontractors have adapted to cope with the changing demands that the rebuild posed. Apart from the magnitude of reconstruction works and new developments that directly affect the capacity of subcontracting businesses in Canterbury, case studies found that subcontractors’ capacity and capability to meet the demand varies and is influenced by the: subcontractors’ own unique characteristics, which are often shaped by changing circumstances in a dynamic and uncertain recovery process; and internal factors in relation to the company’s goal and employees’ needs
Two days after the 22 February 2011 M6.3 earthquake in Christchurch, New Zealand, three of the authors conducted a transect of the central city, with the goal of deriving an estimate of building damage levels. Although smaller in magnitude than the M7.1 4 September 2010 Darfield earthquake, the ground accelerations, ground deformation and damage levels in Christchurch central city were more severe in February 2011, and the central city was closed down to the general public. Written and photographic notes of 295 buildings were taken, including construction type, damage level, and whether the building would likely need to be demolished. The results of the transect compared favourably to Civil Defence rapid assessments made over the following month. Now, more than one year and two major aftershocks after the February 2011 earthquake these initial estimates are compared to the current demolition status to provide an updated understanding of the state of central Christchurch.
The 2011, 6.3 magnitude Christchurch earthquake in New Zealand caused considerable structural damage. It is believed that this event has now resulted in demolition of about 65-70% of the building stock in the Central Business District (CBD), significantly crippling economic activities in the city of Christchurch. A major concern raised from this event was adequacy of the current seismic design practice adopted for reinforced concrete walls due to their poor performance in modern buildings. The relatively short-duration earthquake motion implied that the observed wall damage occurred in a brittle manner despite adopting a ductile design philosophy. This paper presents the lessons learned from the observed wall damage in the context of current state of knowledge in the following areas: concentrating longitudinal reinforcement in wall end regions; determining wall thickness to prevent out-of-plane wall buckling; avoiding lap splices in plastic hinge zones; and quantifying minimum vertical reinforcement. http://www.2eceesistanbul.org/
The paper proposes a simple method for quick post-earthquake assessment of damage and condition of a stock of bridges in a transportation network using seismic data recorded by a strong motion array. The first part of the paper is concerned with using existing free field strong motion recorders to predict peak ground acceleration (PGA) at an arbitrary bridge site. Two methods are developed using artificial neural networks (a single network and a committee of neural networks) considering influential parameters, such as seismic magnitude, hypocentral depth and epicentral distance. The efficiency of the proposed method is explored using actual strong motion records from the devastating 2010 Darfield and 2011 Christchurch earthquakes in New Zealand. In the second part, two simple ideas are outlined how to infer the likely damage to a bridge using either the predicted PGA and seismic design spectrum, or a broader set of seismic metrics, structural parameters and damage indices.
The Canterbury earthquakes are unique in that the there have been a series of major earthquakes, each with their own subsequent aftershock pattern. These have extended from the first large earthquake in September 2010 to currently, at the time of writing, two years later. The last significant earthquake of over magnitude 5.0 on the Richter scale was in May on 2012, and the total number of aftershocks has exceeded 12,000. The consequences, in addition to the loss of life, significant injury and widespread damage, have been far reaching and long term, with detrimental effects and still uncertain effects for many. This provides unique challenges for individuals, communities, organisations and institutions within Canterbury. This document reviews research-based understandings of the concept of resilience. A conceptual model is developed which identifies a number of the factors that influence individual and household resilience. Guided by the model, a series of recommendations are developed for practices that will support individual and household resilience in Canterbury in the aftermath of the 2010-2011 earthquakes.
The city of Christchurch has experienced over 10,000 aftershocks since the 4th of September 2010 earthquake of which approximately 50 have been greater than magnitude 5. The damage caused to URM buildings in Christchurch over this sequence of earthquakes has been well documented. Due to the similarity in age and construction of URM buildings in Adelaide, South Australia and Christchurch (they are sister cities, of similar age and heritage), an investigation was conducted to learn lessons for Adelaide based on the Christchurch experience. To this end, the number of URM buildings in the central business districts of both cities, the extent of seismic strengthening that exists in both cities, and the relative earthquake hazards for both cities were considered. This paper will report on these findings and recommend strategies that the city of Adelaide could consider to significantly reduce the seismic risk posed by URM buildings in future earthquake.
The Canterbury earthquake series of 2010/2011 has turned the city of Christchurch into a full scale natural laboratory testing the structural and non-structural response of buildings under moderate to very severe earthquake shaking. The lessons learned from this, which have come at great cost socially and economically, are extremely valuable in increasing our understanding of whole building performance in severe earthquakes. Given current initiatives underway on both sides of the Tasman towards developing joint Australasian steel and composite steel/concrete design and construction standards that would span a very wide range of geological conditions and seismic zones, these lessons are relevant to both countries. This paper focusses on the performance of steel framed buildings in Christchurch city, with greatest emphasis on multi-storey buildings, but also covering single storey steel framed buildings and light steel framed housing. It addresses such issues as the magnitude and structural impact of the earthquake series, importance of good detailing, lack of observed column base hinging, the excellent performance of composite floors and it will briefly cover research underway to quantify some of these effects for use in design.
On 14 November 2016 a magnitude Mw 7.8 earthquake struck the upper South Island of New Zealand with effects also being observed in the capital city, Wellington. The affected area has low population density but is the largest wine production region in New Zealand and also hosts the main national highway and railway routes connecting the country’s three largest cities of Auckland, Wellington and Christchurch, with Marlborough Port in Picton providing connection between the South and North Islands. These transport facilities sustained substantial earthquake related damage, causing major disruptions. Thousands of landslides and multiple new faults were counted in the area. The winery facilities and a large number of commercial buildings and building components (including brick masonry veneers, historic masonry construction, and chimneys), sustained damage due to the strong vertical and horizontal acceleration. Presented herein are field observations undertaken the day immediately after the earthquake, with the aim to document earthquake damage and assess access to the affected area.
The Christchurch earthquake sequence has been on-going since September 4th 2010. The largest two earthquakes, magnitude (M) 7.1 on September 4th and the M 6.3 on February 22nd 2011 caused immediate and significant damage to the city of Christchurch. As a consequence of the earthquakes, the tourism sector in the Canterbury region has been heavily impacted, with broader impacts being felt throughout the South Island. Resilient Organisations and the University of Canterbury began a series of quantitative investigations into the recovery and response of key business sectors to the earthquakes. The purpose of this study was to build on this work by exploring the outcomes of the earthquakes on the tourism sector, a critical economic driver in the region. Two postal surveys were sent to 719 tourism business managers; the first to businesses in the ‘Impact Zone’ defined as areas that experienced Modified Mercalli intensities greater than 6. The second survey was sent to the remaining businesses throughout the Canterbury region (‘Rest of Canterbury’). Response rates were 46% response for the Impact Zone, and 29% for the Rest of Canterbury. Key findings:
Ground motion observations from the most significant 10 events in the 2010-2011 Canterbury earthquake sequence at near-source sites are utilized to scrutinize New Zealand (NZ)-specific pseudo-spectral acceleration (SA) empirical ground motion prediction equations (GMPE) (Bradley 2010, Bradley 2013, McVerry et al. 2006). Region-specific modification factors based on relaxing the conventional ergodic assumption in GMPE development were developed for the Bradley (2010) model. Because of the observed biases with magnitude and source-to-site distance for the McVerry et al. (2006) model it is not possible to develop region-specific modification factors in a reliable manner. The theory of non-ergodic empirical ground motion prediction is then outlined, and applied to this 10 event dataset to determine systematic effects in the between- and within-event residuals which lead to modifications in the predicted median and standard deviation of the GMPE. By examining these systematic effects over sub-regions containing a total of 20 strong motion stations within the Canterbury area, modification factors for use in region-specific ground motion prediction are proposed. These modification factors, in particular, are suggested for use with the Bradley et al. (2010) model in Canterbury-specific probabilistic seismic hazard analysis (PSHA) to develop revised design response, particularly for long vibration periods.
The majority of Christchurch’s stormwater has historically been discharged untreated directly into urban surface waterways. These receiving waterways have become adversely affected by the contaminants carried in the stormwater, particularly sediment and heavy metals. An event-based contaminant load model was developed to identify the distribution and magnitude of contaminant loads entering the waterway, as well as to assess the reduction in TSS and heavy metal loads that can be achieved by various stormwater management options. The GIS-Excel based model estimates contaminant loads from an individual storm event based on different contributing impervious surfaces and key rainfall characteristics (rainfall intensity, duration, pH and antecedent dry days). It then calculates contaminant reduction loads that could be achieved through source reduction (e.g. green roofs, repainting) as well as from treatment (e.g. raingardens, wet ponds) applied to different surfaces within the catchment. This model differs from other annual load models as it is event-based and accounts for storm characteristics in its calculation of contaminant loads. Christchurch is a valuable case setting due the unique opportunity for retrofitting improved stormwater management in the post-earthquake rebuild. It is anticipated that this modelling approach could later be adapted for use in other urban settings outside of Christchurch.
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.
We present ground motion simulations of the Porters Pass (PP) fault in the Canterbury region of New Zealand; a major active source near Christchurch city. The active segment of the PP fault has an inferred length of 82 km and a mostly strike-slip sense of movement. The PP fault slip makes up approximately 10% of the total 37 mm/yr margin-parallel plate motion and also comprises a significant proportion of the total strain budget in regional tectonics. Given that the closest segment of the fault is less than 45 km from Christchurch city, the PP fault is crucial for accurate earthquake hazard assessment for this major population centre. We have employed the hybrid simulation methodology of Graves and Pitarka (2010, 2015), which combines low (f<1 Hz) and high (f>1 Hz) frequencies into a broadband spectrum. We have used validations from three moderate magnitude events (𝑀𝑤4.6 Sept 04, 2010; 𝑀𝑤4.6 Nov 06, 2010; 𝑀𝑤4.9 Apr 29, 2011) to build confidence for the 𝑀𝑤 > 7 PP simulations. Thus far, our simulations include multiple rupture scenarios which test the impacts of hypocentre location and the finite-fault stochastic rupture representation of the source itself. In particular, we have identified the need to use location-specific 1D 𝑉𝑠/𝑉𝑝 models for the high frequency part of the simulations to better match observations.
The Porter's Pass-Amberley Fault Zone (PPAFZ) is a complex zone of anastomosing faults and folds bounding the south-eastern edge of the transition from subducting Pacific Plate to continental collision on the Australia Plate boundary. This study combines mapping of a 2000 km2 zone from the Southern Alps northeast to the coast near Amberley, 40 km north of metropolitan Christchurch, with an analysis of seismicity and a revision of regional seismic hazard. Three structural styles: 1) a western strike-slip, and 2) a more easterly thrust and reverse domain, pass into 3) a northwest verging fold belt on the northern Canterbury Plains, reflecting the structural levels exposed and the evolving west to east propagation. Basal remnants of a Late Cretaceous-Cenozoic, largely marine sedimentary cover sequence are preserved as outliers that unconformably overlie Mesozoic basement (greywacke and argillite of the Torlesse terrain) in the mountains of the PPAFZ and are underlain by a deeply leached zone which is widely preserved. Structure contouring of the unconformity surface indicates maximum, differential uplift of c.2600 m in the southwest, decreasing to c.1200 m in the coastal fold belt to the northeast. Much lower rates (or reversal) of uplift are evident a few kilometres southeast of the PPAFZ range-front escarpment. The youngest elements of the cover sequence are basement-derived conglomerates of Plio-Pleistocene age preserved on the SE margin. The source is more distant than the intervening mountains of the PPAFZ, probably from the Southern Alps, to the west and northwest. The absence of another regional unconformity on Mesozoic basement, older than Pleistocene, indicates that this uplift is post-Pliocene. Late Pleistocene(<100 kyr) differential uplift rates of c.0.5-2.7 m/kyr from uplifted marine terraces at the east coast, and rates of 2.5-3.3 m/kyr for tectonically-induced river-down cutting further west, suggest that uplift commenced locally during the last 1 Ma, and possibly within the last 0.5 Ma, if average rates are assumed to be uniform over time. Analysis of seismicity, recorded during a 10 week regional survey of micro earthquakes in 1990, identified two seismic zones beneath North Canterbury: 1) a sub-horizontal zone of activity restricted to the upper crust (≤12 km); and 2) a seismic zone in the lower crust (below a ceiling of ≤17 km), that broadens vertically to the north and northwest to a depth of c.40 km, with a bottom edge which dips 10°N and 15°NW, respectively. No events were recorded at depths between 12 km and 17 km, which is interpreted as a relatively aseismic, mid-crustal ductile layer. Marked differences (up to 60°) in the trend of strain axes for events above and below the inferred ductile layer are observed only north of the PPAFZ. A fundamental, north-to-south increase in the Wave-length of major geological structures occurs across the PPAFZ, and is interpreted as evidence that the upper crust beneath the Canterbury Plains is coupled to the lower crust, whereas the upper crust further north is not. Most of the recorded micro earthquakes <12 km deep beneath the PPAFZ have strike-slip mechanisms. It is probable that faults splay upward into the thrusts and folds at the surface as an evolving transpression zone in response to deep shear in basement. There have been no historic surface ruptures of the PPAFZ, but the zone has been characterised historically by frequent small earthquakes. Paleoseismic data (dated landslides and surface ruptures) compiled in this study, indicate a return period of 1500-1900 years between the last two M>7-7.5 earthquakes, and 500-700 years have elapsed since the last. The magnitudes of these events are estimated at c.M7.5, which represents a probable maximum magnitude for the PPAFZ. There are insufficient data to determine whether or not the frequency of large earthquakes conforms to a recognised model of behaviour, but comparison of the paleoseismic data with the historic record of smaller earthquakes, suggests that the magnitudes of the largest earthquakes in this zone are not exponentially distributed. A seismicity model for the PPAFZ (Elder et al., 1991) is reviewed, and a b-value of 1.0 is found to be consistent with the newly acquired paleoseismic data. This b-value reduces the predicted frequency of large earthquakes (M≥7.0) in this zone by a factor of 3.5, while retaining a conservative margin that allows for temporal variations in the frequency of large events and the possibility that the geological database is incomplete, suggesting grounds for revising the hazard model for Christchurch.
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.
