A team of earthquake geologists, seismologists and engineering seismologists from GNS Science, NIWA, University of Canterbury, and Victoria University of Wellington have collectively produced an update of the 2002 national probabilistic seismic hazard (PSH) model for New Zealand. The new model incorporates over 200 new onshore and offshore fault sources, and utilises newly developed New Zealand-based scaling relationships and methods for the parameterisation of the fault and subduction interface sources. The background seismicity model has also been updated to include new seismicity data, a new seismicity regionalisation, and improved methodology for calculation of the seismicity parameters. Background seismicity models allow for the occurrence of earthquakes away from the known fault sources, and are typically modelled as a grid of earthquake sources with rate parameters assigned from the historical seismicity catalogue. The Greendale Fault, which ruptured during the M7.1, 4 September 2010 Darfield earthquake, was unknown prior to the earthquake. However, the earthquake was to some extent accounted for in the PSH model. The maximum magnitude assumed in the background seismicity model for the area of the earthquake is 7.2 (larger than the Darfield event), but the location and geometry of the fault are not represented. Deaggregations of the PSH model for Christchurch at return periods of 500 years and above show that M7-7.5 fault and background source-derived earthquakes at distances less than 40 km are important contributors to the hazard. Therefore, earthquakes similar to the Darfield event feature prominently in the PSH model, even though the Greendale Fault was not an explicit model input.
A magnitude 6.3 earthquake struck the city of Christchurch at 12:51pm on Tuesday 22 February 2011. The earthquake caused 182 fatalities, a large number of injuries, and resulted in widespread damage to the built environment, including significant disruption to the lifelines. The event created the largest lifeline disruption in a New Zealand city in 80 years, with much of the damage resulting from extensive and severe liquefaction in the Christchurch urban area. The Christchurch earthquake occurred when the Canterbury region and its lifelines systems were at the early stage of recovering from the 4 September 2010 Darfield (Canterbury) magnitude 7.1 earthquake. This paper describes the impact of the Christchurch earthquake on lifelines by briefly summarising the physical damage to the networks, the system performance and the operational response during the emergency management and the recovery phase. Special focus is given to the performance and management of the gas, electric and road networks and to the liquefaction ejecta clean-up operations that contributed to the rapid reinstatement of the functionality of many of the lifelines. The water and wastewater system performances are also summarized. Elements of resilience that contributed to good network performance or to efficient emergency and recovery management are highlighted in the paper.
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.
On the second day of teaching for 2011, the University of Canterbury (UC) faced the most significant crisis of its 138-year history. After being shaken severely by a magnitude 7.1 earthquake on 4 September 2010, UC felt it was well along the pathway to getting back to ‘normal’. That all changed at 12:51pm on 22 February 2011, when Christchurch city was hit by an even more devastating event. A magnitude 6.3 (Modified Mercalli intensity ten – MM X) earthquake, just 13km south-east of the Christchurch city centre, caused vertical peak ground accelerations amongst the highest ever recorded in an urban environment, in some places more than twice the acceleration due to gravity. The earthquake caused immediate evacuation of the UC campus and resulted in significant damage to many buildings. Thankfully there were no serious injuries or fatalities on campus, but 185 people died in the city and many more suffered serious injuries. At the time of writing, eighteen months after the first earthquake in September, Christchurch is still experiencing regular earthquakes. Seismologists warn that the region may experience heightened seismicity for a decade or more. While writing this report we have talked with many different people from across the University. People’s experiences are different and we have not managed to talk with everyone, but we hope that by drawing together many different perspectives from across the campus that this report will serve two purposes; to retain our institutional memory of what we have learnt over the past eighteen months, and also to share our learnings with other organisations in New Zealand and around the world who, we hope, will benefit from learning about our experience.
The Canterbury earthquakes, which started with the 7.1 magnitude event on September 4, 2010, caused significant damage in the region. The September 4 earthquakes brought substantial damage to land, buildings, and infrastructure, while the 6.3 magnitude earthquake on February 22, 2011 (and its subsequent aftershocks), brought even greater property damage, but also significant loss of life in addition to the region. Thousands were injured, and 185 persons died. A national State of Emergency was declared and remained in effect until April 30, 2011. A significant number of people required immediate assistance and support to deal with loss, injury, trauma experiences, and property damages. Many had to find alternate accommodation as their houses were too damaged to stay in. Of those affected, many were already vulnerable, and others had been too traumatized by the events to effectively deal with the challenges they were faced with. A number of human service organizations in the region, from both government and non-government sectors, joined forces to be able to more effectively and efficiently help those in need. This was the start of what would become known as the Earthquake Support Coordination Service. The aim of this report is to present an evaluation of the Earthquake Support Coordination Service and its collaborative organization, based on documentation and interviews with key stakeholders of the service. The aim is also to evaluate the service based on perspectives gathered among the clients as well as the coordinators working in the service. The final aim is to offer a reflection on the service model, and on what factors enabled the service, as well as recommendations regarding aspects of the service which may require review, and aspects which may be useful in other contexts.
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:
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.
On Tuesday 22 February 2011, a 6.3 magnitude earthquake struck Christchurch, New Zealand’s second largest city. The ‘earthquake’ was in fact an aftershock to an earlier 7.1 magnitude earthquake that had occurred on Saturday 4 September 2010. There were a number of key differences between the two events that meant they had dramatically different results for Christchurch and its inhabitants. The 22 February 2011 event resulted in one of New Zealand’s worst natural disasters on record, with 185 fatalities occurring and hundreds more being injured. In addition, a large number of buildings either collapsed or were damaged to the point where they needed to be totally demolished. Since the initial earthquake in September 2010, a large amount of building-related research has been initiated in New Zealand to investigate the impact of the series of seismic events – the major focus of these research projects has been on seismic, structural and geotechnical engineering matters. One project, however, conducted jointly by the University of Canterbury, the Fire Protection Association of New Zealand and BRANZ, has focused on the performance of fire protection systems in the earthquakes and the effectiveness of the systems in the event of post-earthquake fires occurring. Fortunately, very few fires actually broke out following the series of earthquake events in Christchurch, but fire after earthquakes still has significant implications for the built environment in New Zealand, and the collaborative research has provided some invaluable insight into the potential threat posed by post-earthquake fires in buildings. As well as summarising the damage caused to fire protection systems, this paper discusses the flow-on effect for designing structures to withstand post-earthquake fires. One of the underlying issues that will be explored is the existing regulatory framework in New Zealand whereby structural earthquake design and structural design for fire are treated as discrete design scenarios.
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.
This research employs a deterministic seismic risk assessment methodology to assess the potential damage and loss at meshblock level in the Christchurch CBD and Mount Pleasant primarily due to building damage caused by earthquake ground shaking. Expected losses in terms of dollar value and casualties are calculated for two earthquake scenarios. Findings are based on: (1) data describing the earthquake ground shaking and microzonation effects; (2) an inventory of buildings by value, floor area, replacement value, occupancy and age; (3) damage ratios defining the performance of buildings as a function of earthquake intensity; (4) daytime and night-time population distribution data and (5) casualty functions defining casualty risk as a function of building damage. A GIS serves as a platform for collecting, storing and analyzing the original and the derived data. It also allows for easy display of input and output data, providing a critical functionality for communication of outcomes. The results of this study suggest that economic losses due to building damage in the Christchurch CBD and Mount Pleasant will possibly be in the order of $5.6 and $35.3 million in a magnitude 8.0 Alpine fault earthquake and a magnitude 7.0 Ashley fault earthquake respectively. Damage to non-residential buildings constitutes the vast majority of the economic loss. Casualty numbers are expected to be between 0 and 10.
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 occurred to unreinforced masonry buildings throughout the region during the mainshock and subsequent large aftershocks. Particularly 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. 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, a geotechnical reconnaissance was conducted over a period of six days (10–15 September 2010) by a team of geotechnical/earthquake engineers and geologists from New Zealand and USA (GEER team: Geo-engineering Extreme Event Reconnaissance). JGS (Japanese Geotechnical Society) members from Japan also participated in the reconnaissance team from 13 to 15 September 2010. The NZ, GEER and JGS members worked as one team and shared resources, information and logistics in order to conduct thorough and most efficient reconnaissance covering a large area over a very limited time period. This report summarises the key evidence and findings from the reconnaissance.
At 4.35am on Saturday 4 September 2010, a magnitude 7.1 earthquake struck near the township of Darfield in Canterbury leading to widespread damage in Christchurch and the wider central Canterbury region. Though it was reported no lives were lost, that was not entirely correct. Over 3,000 animals perished as a result of the earthquake and 99% of these deaths would have been avoidable if appropriate mitigation measures had been in place. Deaths were predominantly due to zoological vulnerability of birds in captive production farms. Other problems included lack of provision of animal welfare at evacuation centres, issues associated with multiple lost and found pet services, evacuation failure due to pet separation and stress impact on dairy herds and associated milk production. The Canterbury Earthquake has highlighted concerns over a lack of animal emergency welfare planning and capacity in New Zealand, an issue that is being progressed by the National Animal Welfare Emergency Management Group. As animal emergency management becomes better understood by emergency management and veterinary professionals, it is more likely that both sectors will have greater demands placed upon them by national guidelines and community expectations to ensure provisions are made to afford protection of animals in times of disaster. A subsequent and more devastating earthquake struck the region on Monday 22 February 2011; this article however is primarily focused on the events pertaining to the September 4 event.
The objective of this study is to examine the influence of near-fault motions on liquefaction triggering in Christchurch and neighboring towns during the 2010-2011 Canterbury earthquake sequence (CES). The CES began with the 4 September 2010, Mw7.1 Darfield earthquake and included up to ten events that triggered liquefaction. However, most notably, widespread liquefaction was induced by the Darfield earthquake and the Mw6.2, 22 February 2011 Christchurch earthquake. Of particular relevance to this study is the forward directivity effects that were prevalent in the motions recorded during the Darfield earthquake, and to a much lesser extent, during the Christchurch earthquake. A 2D variant of the Richart-Newmark fatigue theory was used to compute the equivalent number of cycles (neq) for the ground motions, where volumetric strain was used as the damage metric. This study is unique because it considers the contribution and phasing of both the fault-normal and fault-parallel components of motion on neq and the magnitude scaling factor (MSF). It was found that when the fault-normal and fault-parallel motions were treated individually, the former yielded a lower neq than the latter. Additionally, when the combined effects of fault-normal and fault-parallel components were considered, it was found that the MSF were higher than those commonly used. This implies that motions containing near-fault effects are less demanding on the soil than motions that do not. This may be one of several factors that resulted in less severe liquefaction occurring during the Darfield earthquake than the Christchurch earthquake.
The Amuri Earthquake of September 1, 1888 (magnitude M = 6.5 to 6.8) occurred on the Hope River Segment of the Hope Fault west of Hanmer Plains. The earthquake was felt strongly in North Canterbury and North Westland and caused considerable property damage and landsliding in the Lower Hope Valley. However, damage reports and the spatial distribution of felt intensities emphasize extreme variations in seismic effects over short distances, probably due to topographic focusing and local ground conditions. Significant variations in lateral fault displacement occurred at secondary fault segment boundaries (side-steps and bends in the fault trace) during the 1888 earthquake. This historical spatial variation in lateral slip is matched by the Late Quaternary geomorphic distribution of slip on the Hope River Segment of the Hope Fault. Trenching studies at two sites on the Hope Fault have also identified evidence for five pre-historic earthquakes of similar magnitude to the 1888 earthquake and an average recurrence interval of 134 ± 27 years between events. Magnitude estimates for the 1888 earthquake are combined with a. strong ground motion attenuation expression to provide an estimate of potential ground accelerations in Amuri District during-future earthquakes on the Hope River Segment of the Hope Fault. The predicted acceleration response on bedrock sites within 20 km of the epicentral region is between 0.23 g and 0.34 g. The close match between the historic, inferred pre-historic and geomorphic distribution of lateral slip indicates that secondary fault segmentation exerts a strong structural control on rupture propagation and the expression of fault displacement at the surface. In basement rocks at depth the spatial variations in slip are inferred to be distributed within zones of pervasive cataclastic shear, on either side of the fault segment boundaries. The large variations in surface displacement across fault segment boundaries means that one must know the geometry of the fault in order to evaluate slip-rates calculated from individual locations. The average Late Quaternary slip-rate on the Hope Fault at Glynn Wye Station is between 15.5 mm/yr and 18.25 mm/yr and the rate on the subsidiary Kakapo Fault is between 5.0 mm/yr and 7.5 mm/yr. These rates have been determined from sites which are relatively free of structural complication.
By closely examining the performance of a 22-storey steel framed building in Christchurch subject to various earthquakes over the past seven years, it is shown that a number of lessons can be learnt regarding the cost-effective consideration of non-structural elements. The first point in this work is that non-structural elements significantly affected the costs associated with repairing steel eccentrically braced frame (EBF) links. The decommissioning or rerouting of non-structural elements in the vicinity of damaged links in the case study building attributed to approximately half the total cost of their repair. Such costs could be significantly reduced if the original positioning of non-structural elements took account of the potential need to repair the EBF links. The second point highlighted is the role that pre-cast cladding apparently played on the distribution and type of damage in the building. Loss estimates obtained following the FEMA P-58 framework vary considerably when cladding is or isnt modelled, both because of changes to drift demands up the height of the building and because certain types of subsequent damage are likely to be cheaper to repair than others. Finally, costly repairs to non-structural partition walls were required not only after the moment magnitude 7.1 earthquake in 2010 but also in multiple aftershocks in the years that followed. Repair costs associated with aftershock events exceeded those from the main event, emphasizing the need to consider aftershocks within modern performance-based earthquake engineering and also the opportunity that exists to make more cost-effective repair strategies following damaging earthquakes.
The magnitude Mw7.8 ‘Kaikōura’ earthquake occurred shortly after midnight on 14 November 2016. This paper presents an overview of the geotechnical impacts on the South Island of New Zealand recorded during the postevent reconnaissance. Despite the large moment magnitude of this earthquake, relatively little liquefaction was observed across the South Island, with the only severe manifestation occurring in the young, loose alluvial deposits in the floodplains of the Wairau and Opaoa Rivers near Blenheim. The spatial extent and volume of liquefaction ejecta across South Island is significantly less than that observed in Christchurch during the 2010-2011 Canterbury Earthquake Sequence, and the impact of its occurrence to the built environment was largely negligible on account of the severe manifestations occurring away from the areas of major development. Large localised lateral displacements occurred in Kaikōura around Lyell Creek. The soft fine-grained material in the upper portions of the soil profile and the free face at the creek channel were responsible for the accumulation of displacement during the ground shaking. These movements had severely impacted the houses which were built close (within the zone of large displacement) to Lyell Creek. The wastewater treatment facility located just north of Kaikōura also suffered tears in the liners of the oxidation ponds and distortions in the aeration system due to ground movements. Ground failures on the Amuri and Emu Plains (within the Waiau Valley) were small considering the large peak accelerations (in excess of 1g) experienced in the area. Minor to moderate lateral spreading and ejecta was observed at some bridge crossings in the area. However, most of the structural damage sustained by the bridges was a result of the inertial loading, and the damage resulting from geotechnical issues were secondary.
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.
On 4 September 2010 the Magnitude 7.1 'Darfield' Earthquake marked the beginning of the Canterbury earthquake sequence. The Darfield earthquake produced strong ground shaking throughout the centralCanterbury Plains, affecting rural areas, small towns and the city of Christchurch. The event produced a 29km long surface rupture through intensive farmland, causing localised flooding and liquefaction. The central Canterbury plains were subjected to a sustained period of thousands of aftershocks in the months after the Darfield earthquake. The primary sector is a major component of the in New Zealand economy. Business units are predominantly small family-run farm organisations, though there are increasing levels of corporate farming. The agribusiness sector contributes 20 per cent of real GDP and 47 per cent of total exports for New Zealand. Of the approximately 2,000 farms that are located in the Canterbury Plains, the most common farming sectors in the region are Mixed farming (mostly comprised of sheep and/or beef farming), Dairy farming, and Arable farming (cropping). Many farms on the Canterbury Plains require some form of irrigation and are increasingly capital intensive, reliant on built infrastructure, technology and critical services. Farms are of great significance to their local rural economies, with many rural non-farming organisations dependent on the health of local farming organisations. Despite the economic significance of the sector, there have been few, if any studies analysing how modern intensive farms are affected by earthquakes. The aim of this report is to (1) summarise the impacts the Darfield earthquake had on farming organisations and outline in general terms how farms are vulnerable to the effects of an earthquake; (2) identify what factors helped mitigate earthquake-related impacts. Data for this paper was collected through two surveys of farming and rural non-farming organisations following the earthquake and contextual interviews with affected organisations. In total, 78 organisations participated in the study (Figure 1). Farming organisations represented 72% (N=56) of the sample.
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.
The previously unknown Greendale Fault was buried beneath the Canterbury Plains and ruptured in the September 4th 2010 moment magnitude (Mw) 7.1 Darfield Earthquake. The Darfield Earthquake and subsequent Mw 6 or greater events that caused damage to Christchurch highlight the importance of unmapped faults near urban areas. This thesis examines the morphology, age and origin of the Canterbury Plains together with the paleoseismology and surface-rupture displacement distributions of the Greendale Fault. It offers new insights into the surface-rupture characteristics, paleoseismology and recurrence interval of the Greendale Fault and related structures involved in the 2010 Darfield Earthquake. To help constrain the timing of the penultimate event on the Greendale Fault the origin and age of the faulted glacial outwash deposits have been examined using sedimentological analysis of gravels and optically stimulated luminescence (OSL) dating combined with analysis of GPS and LiDAR survey data. OSL ages from this and other studies, and the analysis of surface paleochannel morphology and subsurface gravel deposits indicate distinct episodes of glacial outwash activity across the Canterbury Plains, at ~20 to 24 and ~28 to 33 kyr separated by a hiatus in sedimentation possibly indicating an interstadial period. These data suggest multiple glacial periods between ~18 and 35 kyr which may have occurred throughout the Canterbury region and wider New Zealand. A new model for the Waimakariri Fan is proposed where aggradation is mainly achieved during episodic sheet flooding with the primary river channel location remaining approximately fixed. The timing, recurrence interval and displacements of the penultimate surface-rupturing earthquake on the Greendale Fault have been constrained by trenching the scarp produced in 2010 at two locations. These excavations reveal a doubling of the magnitude of surface displacement at depths of 2-4 m. Aided by OSL ages of sand lenses in the gravel deposits, this factor-of-two increase is interpreted to indicate that in the central section of the Greendale Fault the penultimate surface-rupturing event occurred between ca. 20 and 30 kyr ago. The Greendale Fault remained undetected prior to the Darfield earthquake because the penultimate fault scarp was eroded and buried during Late Pleistocene alluvial activity. The Darfield earthquake rupture terminated against the Hororata Anticline Fault (HAF) in the west and resulted in up to 400 mm of uplift on the Hororata Anticline immediately above the HAF. Folding in 2010 is compared to Quaternary and younger deformation across the anticline recorded by a seismic reflection line, GPS-measured topographic profiles along fluvial surfaces, and river channel sinuosity and morphology. It is concluded that the HAF can rupture during earthquakes dissimilar to the 2010 event that may not be triggered by slip on the Greendale Fault. Like the Greendale Fault geomorphic analyses provide no evidence for rupture of the HAF in the last 18 kyr, with the average recurrence interval for the late Quaternary inferred to be at least ~10 kyr. Surface rupture of the Greendale Fault during the Darfield Earthquake produced one of the most accessible and best documented active fault displacement and geometry datasets in the world. Surface rupture fracture patterns and displacements along the fault were measured with high precision using real time kinematic (RTK) GPS, tape and compass, airborne light detection and ranging (LiDAR), and aerial photos. This allowed for detailed analysis of the cumulative strike-slip displacement across the fault zone, displacement gradient (ground shear strain) and the type of displacement (i.e. faulting or folding). These strain profiles confirm that the rupture zone is generally wide (~30 to ~300 metres) with >50% of displacement (often 70-80%) accommodated by ground flexure rather than discrete fault slip and ground cracking. The greatest fault-zone widths and highest proportions of folding are observed at fault stepovers.
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.
Oblique-convergent plate collision between the Pacific and Australian plates across the South Island has resulted in shallow, upper crustal earthquake activity and ground surface deformation. In particular the Porters Pass - Amberley Fault Zone displays a complex hybrid zone of anastomosing dextral strike-slip and thrust/reverse faulting which includes the thrust/reverse Lees Valley Fault Zone and associated basin deformation. There is a knowledge gap with respect to the paleoseismicity of many of the faults in this region including the Lees Valley Fault Zone. This study aimed to investigate the earthquake history of the fault at a selected location and the structural and geomorphic development of the Lees Valley Fault Zone and eastern rangefront. This was investigated through extensive structural and geomorphic mapping, GPS field surveying, vertical aerial photo interpretation, analysis of Digital Elevation Models, paleoseismic trenching and optically stimulated luminescence dating. This thesis used a published model for tectonic geomorphology development of mountain rangefronts to understand the development of Lees Valley. Rangefront geomorphology is investigated through analysis of features such as rangefront sinuosity and faceted spurs and indicates the recently active and episodic nature of the uplifted rangefront. Analysis of fault discontinuity, fault splays, distribution of displacement, fault deformation zone and limited exposure of bedrock provided insight into the complex structure of the fault zone. These observations revealed preserved, earlier rangefronts, abandoned and uplifted within the eastern ranges, indicating a basinward shift in focus of faulting and an imbricate thrust wedge development propagating into the footwall of the fault zone and along the eastern ranges of Lees Valley. Fault scarp deformation analysis indicated multiple events have produced the deformation present preserved by the active fault trace in the northern valley. Vertical deformation along this scarp varied with a maximum of 11.5 m and an average of 5 m. Field mapping revealed fan surfaces of various ages have been offset and deformed, likely during the Holocene, based on expected relative surface ages. Geomorphic and structural mapping highlighted the effect of cross-cutting and inherited structures on the Lees Valley Fault, resulting in a step-over development in the centre of the eastern range-bounding trace. Paleoseismic trenching provided evidence of at least two earthquakes, which were constrained to post 21.6 ± 2.3 ka by optically stimulated luminescence dating. Single event displacements (1.48 ± 0.08 m), surface rupture earthquake magnitudes (Mw 6.7 ± 0.1, with potential to produce ≥ 7.0), and a minimum recurrence interval (3.6 ± 0.3 ka) indicated the Lees Valley Fault is an active structure capable of producing significant earthquake events. Results from this study indicate that the Lees Valley Fault Zone accommodates an important component of the Porters Pass - Amberley Fault Zone deformation and confirms the fault as a source of potentially damaging, peak ground accelerations in the Canterbury region. Remnants of previous rangefronts indicate a thrust wedge development of the Lees Valley Fault Zone and associated ranges that can potentially be used as a model of development for other thrust-fault bounded basins.
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.
