Novel Gel-push sampling was employed to obtain high quality samples of Christchurch sands from the Central Business District, at sites where liquefaction was observed in 22 February 2011, and 13 June 2011 earthquakes. The results of cyclic triaxial testing on selected undisturbed specimens of typical Christchurch sands are presented and compared to empirical procedures used by practitioners. This comparison suggests cyclic triaxial data may be conservative, and the Magnitude Scaling Factor used in empirical procedures may be unconservative for highly compressible soils during near source moderate to low magnitude events. Comparison to empirical triggering curves suggests the empirical method generally estimates the cyclic strength of Christchurch sands within a reasonable degree of accuracy as a screening evaluation tool for liquefaction hazard, however for sands with moderate to high fines content it may be significantly unconservative, highlighting the need for high quality sampling and testing on important projects where seismic performance is critical.
The last seven years have seen southern New Zealand a ected by several large and damaging earthquakes: the moment magnitude (MW) 7.8 Dusky Sound earthquake on 15 July 2009, the MW 7.1 Dar eld (Canterbury) earthquake on 4 September 2010, and most notably the MW 6.2 Christchurch earthquake on 22 February 2011 and the protracted aftershock sequence. In this thesis, we address the postseismic displacement produced by these earthquakes using methods of satellite-based geodetic measurement, known as Interferometric Synthetic Aperture Radar (InSAR) and Global Positioning System (GPS), and computational modelling. We observe several ground displacement features in the Canterbury and Fiordland regions during three periods: 1) Following the Dusky Sound earthquake; 2) Following the Dar eld earthquake and prior to the Christchurch earthquake; and 3) Following the Christchurch earthquake until February 2015. The ground displacement associated with postseismic motion following the Dusky Sound earthquake has been measured by continuous and campaign GPS data acquired in August 2009, in conjunction with Di erential Interferometric Synthetic Aperture Radar (DInSAR) observations. We use an afterslip model, estimated by temporal inversion of geodetic data, with combined viscoelastic rebound model to account for the observed spatio-temporal patterns of displacement. The two postseismic processes together induce a signi cant displacement corresponding to principal extensional and contractual strain rates of the order of 10⁻⁷ and 10⁻⁸ yr⁻¹ respectively, across most of the southern South Island. We also analyse observed postseismic displacement following the Dusky Sound earthquake using a new inversion approach in order to describe afterslip in an elasticviscoelastic medium. We develop a mathematical framework, namely the "Iterative Decoupling of Afterslip and Viscoelastic rebound (IDAV)" method, with which to invert temporally dense and spatially sparse geodetic observations. We examine the IDAV method using both numerical and analytical simulations of Green's functions. For the post-Dar eld time interval, postseismic signals are measured within approximately one month of the mainshock. The dataset used for the post-Dar eld displacement spans the region surrounding previously unrecognised faults that ruptured during the mainshock. Poroelastic rebound in a multi-layered half-space and dilatancy recovery at shallow depths provide a satisfactory t with the observations. For the post-Christchurch interval, campaign GPS data acquired in February 2012 to February 2015 in four successive epochs and 66 TerraSAR-X (TSX) SAR acquisitions in descending orbits between March 2011 and May 2014 reveal approximately three years of postseismic displacement. We detect movement away from the satellite of ~ 3 mm/yr in Christchurch and a gradient of displacement of ~ 4 mm/yr across a lineament extending from the westernmost end of the Western Christchurch Fault towards the eastern end of the Greendale East Fault. The postseismic signals following the Christchurch earthquake are mainly accounted for by afterslip models on the subsurface lineament and nearby faults.
Unreinforced masonry (URM) is a construction type that was commonly adopted in New Zealand between the 1880s and 1930s. URM construction is evidently vulnerable to high magnitude earthquakes, with the most recent New Zealand example being the 22 February 2011 Mw6.3 Christchurch earthquake. This earthquake caused significant damage to a majority of URM buildings in the Canterbury area and resulted in 185 fatalities. Many URM buildings still exist in various parts of New Zealand today, and due to their likely poor seismic performance, earthquake assessment and retrofit of the remaining URM building stock is necessary as these buildings have significant architectural heritage and occupy a significant proportion of the nation’s building stock. A collaborative research programme between the University of Auckland and Reid Construction Systems was conducted to investigate an economical yet effective solution for retrofitting New Zealand’s existing URM building stock. This solution adopts the shotcrete technique using an Engineered Cementitious Composite (ECC), which is a polyvinyl alcohol fibre reinforced mortar that exhibits strain hardening characteristics. Collaborations have been formed with a number of consulting structural engineers throughout New Zealand to develop innovative and cost effective retrofit solutions for a number of buildings. Two such case studies are presented in this paper. http://www.concrete2013.com.au/technical-program/
Live monitoring data and simple dynamic reduced-order models of the Christchurch Women’s Hospital (CWH) help explain the performance of the base isolation (BI) system of the hospital during the series of Canterbury earthquakes in 2011-2012. A Park-Wen-Ang hysteresis model is employed to simulate the performance of the BI system and results are compared to measured data recorded above the isolation layer and on the 6th story. Simplified single, two and three degree of freedom models (SDOF, 2DOF and 3DOF) show that the CWH structure did not behave as an isolated but as a fixed-base structure. Comparisons of accelerations and deflections between simulated and monitored data show a good match for isolation stiffness values of approximately two times of the value documented in the design specification and test protocol. Furthermore, an analysis of purely measured data revealed very little to no relative motion across the isolators for large events of moment magnitude scale (Mw) 5.8 and 6.0 that occurred within 3 hours of each other on December 23, 2011. One of the major findings is that the BI system during the seismic events on December 23, 2011 did not yield and that the superstructure performed as a fixed-base building, indicating a need to reevaluate the analysis, design and implementation of these structures.
A zone of active tectonism occurs in mid and north Canterbury, from the Rakaia to the Waipara Rivers, which coincides with seismicity concentrations and several Quaternary surface anomalies and is here defined as the Porters Pass Tectonic Zone. Although parallel to the Marlborough faults to the north, the lack of regional definition suggests this zone is much younger in its inception reflecting a southward movement of the plate rotation vector. The objectives of this study were to map the structures associated with this zone in the segment between the Rakaia and Waimakariri Rivers with detailed analysis concentrated in the upper Kawai Valley. Quaternary offsets on the main lineament of the Porters Pass Fault were traced through the area and evidence for the rate of movement, probable magnitudes and return periods of related seismic events was sought. The basement was found to be complicated by pre-existing deformation structures in Torlesse Group rocks which have been subsequently been re-activated or rotated by recent fault movement probably beginning in the Pleistocene. This phase is dominantly thrusting and uplift has lead to the erosion of most of the overlying sedimentary cover. Remnants of the Cret-Tertiary sediments still remain as fault-bounded packets. Evidence suggests that a change to development of a regional lateral shear associated with the Porters Pass Tectonic Zone transects the thrust system with complex interaction between the older reverse and new strike-slip faults. Offset rates along the segments of the Porters Pass Fault are not well constrained but are believed to be approximately in the range of 11-13 mm/year for at least the last 130,000 years. This rate is similar to other large faults in the Marlborough region. Two earthquake events have been identified and dated at 600 and 2000 years ago, with a magnitude of greater than 6.5. Evidence suggests characteristic earthquakes along the Porters Pass Fault are greater than Magnitude 7. This result has some major ramifications for the expected seismic hazards for nearby 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.
During 2010 and 2011, a series of major earthquakes caused widespread damage in the city of Christchurch, New Zealand. The magnitude 6.3 quake in February 2011 caused 185 fatalities. In the ensuing months, the government progressively zoned residential land in Christchurch on the basis of its suitability for future occupation (considering damage from these quakes and future earthquake risk). Over 6,000 homes were placed in the ‘red-zone’, meaning that property owners were forced to sell their land to the Crown. This study analysed patterns of residential mobility amongst thirty-one red-zone households from the suburb of Southshore, Christchurch. Drawing on interviews and surveys, the research traced their experience from the zoning announcement until they had moved to a new residence. The research distinguished between short (before the zoning announcement) and long term (post the red zone ‘deadline’) forms of household relocation. The majority of households in the study were highly resistant to short term movement. Amongst those which did relocate before the zoning decision, the desire to maintain a valued social connection with a person outside of the earthquake environment was often an important factor. Some households also moved out of perceived necessity (e.g. due to lack of power or water). In terms of long-term relocation, concepts of affordability and safety were much more highly valued by the sample when purchasing post-quake property. This resulted in a distinct patterning of post-quake housing location choices. Perceived control over the moving process, relationship with government organisations and insurance companies, and time spent in the red-zone before moving all heavily influenced participants’ disaster experience. Contrary to previous studies, households in this study recorded higher levels of subjective well-being after relocating. The study proposed a typology of movers in the Christchurch post-disaster environment. Four mobility behaviours, or types, are identified: the Committed Stayers (CSs), the Environment Re-Creators (ERCs), the Resigned Acceptors (RAs), and the Opportunistic Movers (OMs). The CSs were defined by their immobility rather than their relocation aspirations, whilst the ERCs attempted to recreate or retain aspects of Southshore through their mobility. The RAs expressed a form of apathy towards the post-quake environment, whereas, on the other hand, the OMs moved relative to pre-earthquake plans, or opportunities that arose from the earthquake itself. Possibilities for further research include examining household adaptability to new residential environments and tracking further mobility patterns in the years following relocation from the red- zone.
This study followed two similarly affected, but socio-economically disparate suburbs as residents responded to and attempted to recover from the devastating 6.3 magnitude earthquake that struck Christchurch, New Zealand, on February 22, 2011. More specifically, it focuses on the role of local churches, community-based organisations (CBOs) and non-governmental organisations (NGOs), here referred to broadly as civil society, in meeting the immediate needs of local residents and assisting with the longer-term recovery of each neighbourhood. Despite considerable socioeconomic differences between the two neighbourhoods, civil society in both suburbs has been vital in addressing the needs of locals in the short and long term following the earthquake. Institutions were able to utilise local knowledge of both residents and the extent of damage in the area to a) provide a swifter local response than government or civil defence and then help direct the relief these agencies provided locally; b) set up central points for distribution of supplies and information where locals would naturally gather; c) take action on what were perceived to be unmet needs; and d) act as a way of bridging locals to a variety of material, informational, and emotional resources. However the findings also support literature which indicates that other factors are also important in understanding neighbourhood recovery and the role of civil society, including: local leadership; a shared, place-based identity; the type and form of civil society organizations; social capital; and neighbourhood- and household-level indicators of relative vulnerability and inequality. The intertwining of these various factors seems to influence how these neighbourhoods have coped with and taken steps in recovering from this disaster. It is recommended that future research be directed at developing a better understanding of how this occurs. It is suggested that a model similar to Yasui’s (2007) Community Vulnerability and Capacity model be developed as a useful way to approach future research in this area.
During many years the analysis of some geophysical results of Charles Darwin was being carried out in Department. Darwin has connected almost 200 years ago results of catastrophic earthquakes with vertical movement of a surface of the Earth. Usually this movement less horizontal movement and its influence on destruction of cities is not considered. Earthquake hazard assessment studies were focused usually on the horizontal ground motion. Effects of the strong vertical motion were not, practically, discussed. The margins of safety against gravity-induced static vertical forces in constructed buildings usually provide adequate resistance to dynamic forces induced by the vertical acceleration during an earthquake. However, the earthquake in Christchurch is an example of the vertical seismic shock . The earthquake magnitude was rather small - nearby 6.3. However, the result was catastrophic. The same took place in 1835. It allowed to Darwin to formulate a few great ideas. Charles Darwin has explained qualitatively results of an interaction of huge seismic waves with volcanoes and the nature of volcanism and seismicity of our planet. These important data of Charles Darwin became very actual recently. It is possible to tell also the same about tsunami and extreme ocean waves described by Charles Darwin. Therefore this data were analyzed using modern mechanics, mathematics and physics in Department. In particular, the theory of catastrophic waves was developed based on Darwin's data. The theory tried to explain occurrence, evolution and distribution the catastrophic waves in various natural systems, since atoms, oceans, surfaces of the Earth and up to the very early Universe. Some results of the research were published in prestigious magazines. Later they were presented in two books devoted to Charles Darwin's anniversary (2009). Last from them was published in Russian (2011). We give here key ideas of this research which is a part of interdisciplinary researches of Department. Some ideas are discussed. Not less important purpose is very short historical review of some researches of Darwin. In particular, we underline Darwin' priority in the formulation of the bases of Dynamics Earth.
On February 22, 2011, a magnitude Mw 6.2 earthquake affected the Canterbury region, New Zealand, resulting in many fatalities. Liquefaction occurred across many areas, visible on the surface as ‘‘sand volcanoes’’, blisters and subsidence, causing significant damage to buildings, land and infrastructure. Liquefaction occurred at a number of sites across the Christchurch Boys High School sports grounds; one area in particular contained a piston ground failure and an adjacent silt volcano. Here, as part of a class project, we apply near-surface geophysics to image these two liquefaction features and determine whether they share a subsurface connection. Hand auger results enable correlation of the geophysical responses with the subsurface stratigraphy. The survey results suggest that there is a subsurface link, likely via a paleo-stream channel. The anomalous responses of the horizontal loop electromagnetic survey and electrical resistivity imaging highlight the disruption of the subsurface electrical properties beneath and between the two liquefaction features. The vertical magnetic gradient may also show a subtle anomalous response in this area, however the results are inconclusive. The ground penetrating radar survey shows disruption of the subsurface stratigraphy beneath the liquefaction features, in particular sediment mounding beneath the silt ejection (‘‘silt volcano’’) and stratigraphic disruption beneath the piston failure. The results indicate how near-surface geophysics allow the characteristics of liquefaction in the subsurface to be better understood, which could aid remediation work following liquefaction-induced land damage and guide interpretation of geophysical surveys of paleoliquefaction features.
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.
The Canterbury Region is susceptible to a variety of natural hazards, including earthquakes, landslides and climate hazards. Increasing population and tourism within the region is driving development pressures and as more and more development occurs, the risk from natural hazards increases. In order to avoid development occurring in unacceptably vulnerable locations, natural hazard assessments are required. This study is a reconnaissance natural hazard assessment of Lakes Lyndon, Coleridge and Tekapo. There is restricted potential for development at Lake Lyndon, because the land surrounding the lake is owned by the Crown and has a number of development restrictions. However, there is the potential for conservation or recreation-linked development to occur. There is more potential for development at Lake Coleridge. Most of the land surrounding the lake is privately owned and has less development restrictions. The majority of land surrounding Lake Tekapo is divided into Crown-owned pastoral leases, which are protected from development, such as subdivision. However, there are substantial areas around the lake, which are privately owned and, therefore, have potential for development. Earthquake, landslide and climate hazards are the main natural hazards threatening Lakes Lyndon, Coleridge and Tekapo. The lakes are situated in a zone of active earth deformation in which large and relatively frequent earthquakes are produced. A large number of active faults lie within 15 km of each lake, which are capable of producing M7 or larger earthquakes. Ground shaking, liquefaction, landslides, tsunami and seiches are among the consequences of earthquakes, all of which have the potential to cause severe damage to lives, lifelines and infrastructure. Landslides are also common in the landscape surrounding the lakes. The majority of slopes surrounding the lakes are at significant risk from earthquake-induced failure under moderate to strong earthquake shaking. This level of shaking is expected to occur in any 50 year period around Lakes Lyndon and Coleridge, and in any 150 year period around Lake Tekapo. Injuries, fatalities and property damage can occur directly from landslide impact or from indirect effects such as flooding from landslide-generated tsunami or from landslide dam outbreaks. Lakes Lyndon, Coleridge and Tekapo are also susceptible to climate hazards, such as high winds, drought, heavy snowfall and heavy rainfall, which can lead to landslides and flooding. Future climate change due to global warming is most likely going to affect patterns of frequency and magnitudes of extreme weather events, leading to an increase in climate hazards. Before development is permitted around the lakes, it is essential that each of these hazards is considered so that unacceptably vulnerable areas can be avoided.
The magnitude 6.2 Christchurch earthquake struck the city of Christchurch at 12:51pm on February 22, 2011. The earthquake caused 186 fatalities, a large number of injuries, and resulted in widespread damage to the built environment, including significant disruption to lifeline networks and health care facilities. Critical facilities, such as public and private hospitals, government, non-government and private emergency services, physicians’ offices, clinics and others were severely impacted by this seismic event. Despite these challenges many systems were able to adapt and cope. This thesis presents the physical and functional impact of the Christchurch earthquake on the regional public healthcare system by analysing how it adapted to respond to the emergency and continued to provide health services. Firstly, it assesses the seismic performance of the facilities, mechanical and medical equipment, building contents, internal services and back-up resources. Secondly, it investigates the reduction of functionality for clinical and non-clinical services, induced by the structural and non-structural damage. Thirdly it assesses the impact on single facilities and the redundancy of the health system as a whole following damage to the road, power, water, and wastewater networks. Finally, it assesses the healthcare network's ability to operate under reduced and surged conditions. The effectiveness of a variety of seismic vulnerability preparedness and reduction methods are critically reviewed by comparing the observed performances with the predicted outcomes of the seismic vulnerability and disaster preparedness models. Original methodology is proposed in the thesis which was generated by adapting and building on existing methods. The methodology can be used to predict the geographical distribution of functional loss, the residual capacity and the patient transfer travel time for hospital networks following earthquakes. The methodology is used to define the factors which contributed to the overall resilence of the Canterbury hospital network and the areas which decreased the resilence. The results show that the factors which contributed to the resilence, as well as the factors which caused damage and functionality loss were difficult to foresee and plan for. The non-structural damage to utilities and suspended ceilings was far more disruptive to the provision of healthcare than the minor structural damage to buildings. The physical damage to the healthcare network reduced the capacity, which has further strained a health care system already under pressure. Providing the already high rate of occupancy prior to the Christchurch earthquake the Canterbury healthcare network has still provided adequate healthcare to the community.
Following the Mw 6.2 Christchurch Earthquake on 22 February 2011, extensive ground cracking in loessial soils was reported in some areas of the Port Hills, southeast of central Christchurch. This study was undertaken to investigate the mechanisms of earthquake-induced ground damage on the eastern side of the Hillsborough Valley. A zone of extensional cracking up to 40m wide and 600m long was identified along the eastern foot-slope, accompanied by compression features and spring formation at the toe of the slope. An engineering geological and geomorphological model was developed for the eastern Hillsborough Valley that incorporates geotechnical investigation data sourced from the Canterbury Geotechnical Database (CGD), the findings of trenching and seismic refraction surveying carried out for this research, and interpretation of historical aerial photographs. The thickness and extent of a buried peat swamp at the base of the slope was mapped, and found to coincide with significant compression features. Ground cracking was found to have occurred entirely within loess-colluvium and to follow the apices of pre-1920s tunnel-gully fan debris at the southern end of the valley. The ground-cracking on the eastern side of the Hillsborough Valley is interpreted to have formed through tensile failure of the loess-colluvium. Testing was carried out to determine the tensile strength of Port Hills loess colluvium as a function of water content and density, in order to better understand the occurrence and distribution of the observed ground cracking. A comprehensive review of the soil tensile strength testing literature was undertaken, from which a test methodology was developed. Results show remoulded loess-colluvium to possess tensile strength of 7 - 28 kPa across the range of tested moisture contents (10-15%) and dry densities (1650-1900kg/m3). A positive linear relationship was observed between tensile strength and dry density, and a negative linear relationship between moisture content and tensile strength. The observed ground damage and available geotechnical information (inclinometer and piezometer records provided by the Earthquake Commission) were together used to interpret the mechanism(s) of slope movement that occurred in the eastern Hillsborough Valley. The observed ground damage is characteristic of translational movement, but without the development of lateral release scarps, or a basal sliding surface - which was not located during drilling. It is hypothesised that shear displacement has been accommodated by multiple slip surfaces of limited extent within the upper 10m of the slope. Movement has likely occurred within near-saturated colluvial units that have lost strength during earthquake shaking. The eastern Hillsborough Valley is considered to be an ‘incipient translational slide’, as both the patterns of damage and shearing are consistent with the early stages of such slide development. Sliding block analysis was utilised to understand how the eastern Hillsborough Valley may perform in a future large magnitude earthquake. Known cumulative displacements of ~0.3m for eastern Hillsborough Valley during the 2010-2011 Canterbury Earthquake Sequence were compared with modelled slope displacements to back-analyse a lower-bound yield acceleration of 0.2 - 0.25g. Synthetic broadband modelling for future Alpine and Hope Fault earthquakes indicates PGAs of approximately 0.08g for soil sites in the Christchurch area, as such, slope movement is unlikely to be reactivated by an Alpine Fault or Hope Fault earthquake. This does not take into account the possible role of strength loss due to excess pore pressure that may occur during these future events.
Since the early 1980s seismic hazard assessment in New Zealand has been based on Probabilistic Seismic Hazard Analysis (PSHA). The most recent version of the New Zealand National Seismic Hazard Model, a PSHA model, was published by Stirling et al, in 2012. This model follows standard PSHA principals and combines a nation-wide model of active faults with a gridded point-source model based on the earthquake catalogue since 1840. These models are coupled with the ground-motion prediction equation of McVerry et al (2006). Additionally, we have developed a time-dependent clustering-based PSHA model for the Canterbury region (Gerstenberger et al, 2014) in response to the Canterbury earthquake sequence. We are now in the process of revising that national model. In this process we are investigating several of the fundamental assumptions in traditional PSHA and in how we modelled hazard in the past. For this project, we have three main focuses: 1) how do we design an optimal combination of multiple sources of information to produce the best forecast of earthquake rates in the next 50 years: can we improve upon a simple hybrid of fault sources and background sources, and can we better handle the uncertainties in the data and models (e.g., fault segmentation, frequency-magnitude distributions, time-dependence & clustering, low strain-rate areas, and subduction zone modelling)? 2) developing revised and new ground-motion predictions models including better capturing of epistemic uncertainty – a key focus in this work is developing a new strong ground motion catalogue for model development; and 3) how can we best quantify if changes we have made in our modelling are truly improvements? Throughout this process we are working toward incorporating numerical modelling results from physics based synthetic seismicity and ground-motion models.
In September 2010 and February 2011 the Canterbury region of New Zealand was struck by two powerful earthquakes, registering magnitude 7.1 and 6.3 respectively on the Richter scale. The second earthquake was centred 10 kilometres south-east of the centre of Christchurch (the region’s capital and New Zealand’s third most populous urban area, with approximately 360,000 residents) at a depth of five kilometres. 185 people were killed, making it the second deadliest natural disaster in New Zealand’s history. (66 people were killed in the collapse of one building alone, the six-storey Canterbury Television building.) The earthquake occurred during the lunch hour, increasing the number of people killed on footpaths and in buses and cars by falling debris. In addition to the loss of life, the earthquake caused catastrophic damage to both land and buildings in Christchurch, particularly in the central business district. Many commercial and residential buildings collapsed in the tremors; others were damaged through soil liquefaction and surface flooding. Over 1,000 buildings in the central business district were eventually demolished because of safety concerns, and an estimated 70,000 people had to leave the city after the earthquakes because their homes were uninhabitable. The New Zealand Government declared a state of national emergency, which stayed in force for ten weeks. In 2014 the Government estimated that the rebuild process would cost NZ$40 billion (approximately US$27.3 billion, a cost equivalent to 17% of New Zealand’s annual GDP). Economists now estimate it could take the New Zealand economy between 50 and 100 years to recover. The earthquakes generated tens of thousands of insurance claims, both against private home insurance companies and against the New Zealand Earthquake Commission, a government-owned statutory body which provides primary natural disaster insurance to residential property owners in New Zealand. These ranged from claims for hundreds of millions of dollars concerning the local port and university to much smaller claims in respect of the thousands of residential homes damaged. Many of these insurance claims resulted in civil proceedings, caused by disputes about policy cover, the extent of the damage and the cost and/or methodology of repairs, as well as failures in communication and delays caused by the overwhelming number of claims. Disputes were complicated by the fact that the Earthquake Commission provides primary insurance cover up to a monetary cap, with any additional costs to be met by the property owner’s private insurer. Litigation funders and non-lawyer claims advocates who took a percentage of any insurance proceeds also soon became involved. These two factors increased the number of parties involved in any given claim and introduced further obstacles to resolution. Resolving these disputes both efficiently and fairly was (and remains) central to the rebuild process. This created an unprecedented challenge for the justice system in Christchurch (and New Zealand), exacerbated by the fact that the Christchurch High Court building was itself damaged in the earthquakes, with the Court having to relocate to temporary premises. (The High Court hears civil claims exceeding NZ$200,000 in value (approximately US$140,000) or those involving particularly complex issues. Most of the claims fell into this category.) This paper will examine the response of the Christchurch High Court to this extraordinary situation as a case study in innovative judging practices and from a jurisprudential perspective. In 2011, following the earthquakes, the High Court made a commitment that earthquake-related civil claims would be dealt with as swiftly as the Court's resources permitted. In May 2012, it commenced a special “Earthquake List” to manage these cases. The list (which is ongoing) seeks to streamline the trial process, resolve quickly claims with precedent value or involving acute personal hardship or large numbers of people, facilitate settlement and generally work proactively and innovatively with local lawyers, technical experts and other stakeholders. For example, the Court maintains a public list (in spreadsheet format, available online) with details of all active cases before the Court, listing the parties and their lawyers, summarising the facts and identifying the legal issues raised. It identifies cases in which issues of general importance have been or will be decided, with the expressed purpose being to assist earthquake litigants and those contemplating litigation and to facilitate communication among parties and lawyers. This paper will posit the Earthquake List as an attempt to implement innovative judging techniques to provide efficient yet just legal processes, and which can be examined from a variety of jurisprudential perspectives. One of these is as a case study in the well-established debate about the dialogic relationship between public decisions and private settlement in the rule of law. Drawing on the work of scholars such as Hazel Genn, Owen Fiss, David Luban, Carrie Menkel-Meadow and Judith Resnik, it will explore the tension between the need to develop the law through the doctrine of precedent and the need to resolve civil disputes fairly, affordably and expeditiously. It will also be informed by the presenter’s personal experience of the interplay between reported decisions and private settlement in post-earthquake Christchurch through her work mediating insurance disputes. From a methodological perspective, this research project itself gives rise to issues suitable for discussion at the Law and Society Annual Meeting. These include the challenges in empirical study of judges, working with data collected by the courts and statistical analysis of the legal process in reference to settlement. September 2015 marked the five-year anniversary of the first Christchurch earthquake. There remains widespread dissatisfaction amongst Christchurch residents with the ongoing delays in resolving claims, particularly insurers, and the rebuild process. There will continue to be challenges in Christchurch for years to come, both from as-yet unresolved claims but also because of the possibility of a new wave of claims arising from poor quality repairs. Thus, a final purpose of presenting this paper at the 2016 Meeting is to gain the benefit of other scholarly perspectives and experiences of innovative judging best practice, with a view to strengthening and improving the judicial processes in Christchurch. This Annual Meeting of the Law and Society Association in New Orleans is a particularly appropriate forum for this paper, given the recent ten year anniversary of Hurricane Katrina and the plenary session theme of “Natural and Unnatural Disasters – human crises and law’s response.” The presenter has a personal connection with this theme, as she was a Fulbright scholar from New Zealand at New York University in 2005/2006 and participated in the student volunteer cleanup effort in New Orleans following Katrina. http://www.lawandsociety.org/NewOrleans2016/docs/2016_Program.pdf
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.
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.
Critical infrastructure networks are highly relied on by society such that any disruption to service can have major social and economic implications. Furthermore, these networks are becoming increasingly dependent on each other for normal operation such that an outage or asset failure in one system can easily propagate and cascade across others resulting in widespread disruptions in terms of both magnitude and spatial reach. It is the vulnerability of these networks to disruptions and the corresponding complexities in recovery processes which provide direction to this research. This thesis comprises studies contributing to two areas (i) the modelling of national scale in-terdependent infrastructure systems undergoing major disruptions, and (ii) the tracking and quantification of infrastructure network recovery trajectories following major disruptions. Firstly, methods are presented for identifying nationally significant systemic vulnerabilities and incorporating expert knowledge into the quantification of infrastructure interdependency mod-elling and simulation. With application to the interdependent infrastructures networks across New Zealand, the magnitudes and spatial extents of disruption are investigated. Results high-light the importance in considering interdependencies when assessing disruptive risks and vul-nerabilities in disaster planning applications and prioritising investment decisions for enhancing resilience of national networks. Infrastructure dependencies are further studied in the context of recovery from major disruptions through the analysis of curves measuring network functionality over time. Continued studies into the properties of recovery curves across a database of global natural disasters produce statistical models for predicting the trajectory and expected recovery times. Finally, the use of connectivity based metrics for quantifying infrastructure system functionality during recovery are considered with a case study application to the Christchurch Earthquake (February 22, 2011) wastewater network response.
When the devastating 6.3 magnitude earthquake hit Christchurch, Aotearoa New Zealand, at 12.51pm on 22nd February 2011, the psychological and physical landscape was irrevocably changed. In the days and weeks following the disaster communities were isolated due to failed infrastructure, continuing aftershocks and the extensive search and rescue effort which focussed resources on the central business district. In such moments the resilience of a community is truly tested. This research discusses the role of grassroots community groups in facilitating community resilience during the Christchurch 2010/11 earthquakes and the role of place in doing so. I argue that place specific strategies for urban resilience need to be enacted from a grassroots level while being supported by broader policies and agencies. Using a case study of Project Lyttelton – a group aspiring towards a resilient sustainable future who were caught at the epicentre of the February earthquake – I demonstrate the role of a community group in creating resilience through self-organised place specific action during a disaster. The group provided emotional care, basic facilities and rebuilding assistance to the residents of Lyttelton, proving to be an invaluable asset. These actions are closely linked to the characteristics of social support and social learning that have been identified as important to socio-ecological resilience. In addition this research will seek to understand and explore the nuances of place and identity and its role in shaping resilience to such dis-placing events. Drawing on community narratives of the displacement of place identity, the potential for a progressive sense of place as instigated by local groups will be investigated as an avenue for adaptation by communities at risk of disaster and place destabilisation.
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.
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.
This dissertation addresses a diverse range of topics in the area of physics-based ground motion simulation with particular focus on the Canterbury, New Zealand region. The objectives achieved provide the means to perform hybrid broadband ground motion simulation and subsequently validates the simulation methodology employed. In particu- lar, the following topics are addressed: the development of a 3D seismic velocity model of the Canterbury region for broadband ground motion simulation; the development of a 3D geologic model of the interbedded Quaternary formations to provide insight on observed ground motions; and the investigation of systematic effects through ground motion sim- ulation of small-to-moderate magnitude earthquakes. The paragraphs below outline each contribution in more detail. As a means to perform hybrid broadband ground motion simulation, a 3D model of the geologic structure and associated seismic velocities in the Canterbury region is devel- oped utilising data from depth-converted seismic reflection lines, petroleum and water well logs, cone penetration tests, and implicitly guided by existing contour maps and geologic cross sections in data sparse subregions. The model explicitly characterises five significant and regionally recognisable geologic surfaces that mark the boundaries between geologic units with distinct lithology and age, including the Banks Peninsula volcanics, which are noted to strongly influence seismic wave propagation. The Basement surface represents the base of the Canterbury sedimentary basin, where a large impedance contrast exists re- sulting in basin-generated waves. Seismic velocities for the lithological units between the geologic surfaces are derived from well logs, seismic reflection surveys, root mean square stacking velocities, empirical correlations, and benchmarked against a regional crustal model, thus providing the necessary information for a Canterbury velocity model for use in broadband seismic wave propagation. A 3D high-resolution model of the Quaternary geologic stratigraphic sequence in the Canterbury region is also developed utilising datasets of 527 high-quality water well logs, and 377 near-surface cone penetration test records. The model, developed using geostatistical Kriging, represents the complex interbedded regional Quaternary geology by characterising the boundaries between significant interbedded geologic formations as 3D surfaces including explicit modelling of the formation unconformities resulting from the Banks Peninsula volcanics. The stratigraphic layering present can result in complex wave propagation. The most prevalent trend observed in the surfaces was the downward dip from inland to the eastern coastline as a result of the dominant fluvial depositional environment of the terrestrial gravel formations. The developed model provides a benefi- cial contribution towards developing a comprehensive understanding of recorded ground motions in the region and also providing the necessary information for future site char- acterisation and site response analyses. To highlight the practicality of the model, an example illustrating the role of the model in constraining surface wave analysis-based shear wave velocity profiling is illustrated along with the calculation of transfer functions to quantify the effect of the interbedded geology on wave propagation. Lastly, an investigation of systematic biases in the (Graves and Pitarka, 2010, 2015) ground motion simulation methodology and the specific inputs used for the Canterbury region is presented considering 144 small-to-moderate magnitude earthquakes. In the simulation of these earthquakes, the 3D Canterbury Velocity Model, developed as a part of this dissertation, is used for the low-frequency simulation, and a regional 1D velocity model for the high-frequency simulation. Representative results for individual earthquake sources are first presented to highlight the characteristics of the small-to-moderate mag- nitude earthquake simulations through waveforms, intensity measure scaling with source- to-site distance, and spectral bias of the individual events. Subsequently, a residual de- composition is performed to examine the between- and within-event residuals between observed data, and simulated and empirical predictions. By decomposing the residuals into between- and within-event residuals, the biases in source, path and site effects, and their causes, can be inferred. The residuals are comprehensively examined considering their aggregated characteristics, dependence on predictor variables, spatial distribution, and site-specific effects. The results of the simulation are also benchmarked against empir- ical ground motion models, where their similarities manifest from common components in their prediction. Ultimately, suggestions to improve the predictive capability of the simulations are presented as a result of the analysis.
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.
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.
This thesis documents the development and demonstration of an assessment method for analysing earthquake-related damage to concrete waste water gravity pipes in Christchurch, New Zealand, following the 2010-2011 Canterbury Earthquake Sequence (CES). The method is intended to be internationally adaptable to assist territorial local authorities with improving lifelines infrastructure disaster impact assessment and improvements in resilience. This is achieved through the provision of high-resolution, localised damage data, which demonstrate earthquake impacts along the pipe length. The insights gained will assist decision making and the prioritisation of resources following earthquake events to quickly and efficiently restore network function and reduce community impacts. The method involved obtaining a selection of 55 reinforced concrete gravity waste water pipes with available Closed-Circuit Television (CCTV) inspection footage filmed before and after the CES. The pipes were assessed by reviewing the recordings, and damage was mapped to the nearest metre along the pipe length using Geographic Information Systems. An established, systematic coding process was used for reporting the nature and severity of the observed damage, and to differentiate between pre-existing and new damage resulting from the CES. The damage items were overlaid with geospatial data such as Light Detection and Ranging (LiDAR)-derived ground deformation data, Liquefaction Resistance Index data and seismic ground motion data (Peak Ground acceleration and Peak Ground Velocity) to identify potential relationships between these parameters and pipe performance. Initial assessment outcomes for the pipe selection revealed that main pipe joints and lateral connections were more vulnerable than the pipe body during a seismic event. Smaller diameter pipes may also be more vulnerable than larger pipes during a seismic event. Obvious differential ground movement resulted in increased local damage observations in many cases, however this was not obvious for all pipes. Pipes with older installation ages exhibited more overall damage prior to a seismic event, which is likely attributable to increased chemical and biological deterioration. However, no evidence was found relating pipe age to performance during a seismic event. No evidence was found linking levels of pre-CES damage in a pipe with subsequent seismic performance, and seismic performance with liquefaction resistance or magnitude of seismic ground motion. The results reported are of limited application due to the small demonstration sample size, but reveal the additional level of detail and insight possible using the method presented in this thesis over existing assessment methods, especially in relation to high resolution variations along the length of the pipe such as localised ground deformations evidenced by LiDAR. The results may be improved by studying a larger and more diverse sample pool, automating data collection and input processes in order to improve efficiency and consider additional input such as pipe dip and cumulative damage over a large distance. The method is dependent on comprehensive and accurate pre-event CCTV assessments and LIDAR data so that post-event data could be compared. It is proposed that local territorial authorities should prioritise acquiring this information as a first important step towards improving the seismic resilience of a gravity waste water pipe network.
This thesis documents the development and demonstration of an assessment method for analysing earthquake-related damage to concrete waste water gravity pipes in Christchurch, New Zealand, following the 2010-2011 Canterbury Earthquake Sequence (CES). The method is intended to be internationally adaptable to assist territorial local authorities with improving lifelines infrastructure disaster impact assessment and improvements in resilience. This is achieved through the provision of high-resolution, localised damage data, which demonstrate earthquake impacts along the pipe length. The insights gained will assist decision making and the prioritisation of resources following earthquake events to quickly and efficiently restore network function and reduce community impacts. The method involved obtaining a selection of 55 reinforced concrete gravity waste water pipes with available Closed-Circuit Television (CCTV) inspection footage filmed before and after the CES. The pipes were assessed by reviewing the recordings, and damage was mapped to the nearest metre along the pipe length using Geographic Information Systems. An established, systematic coding process was used for reporting the nature and severity of the observed damage, and to differentiate between pre-existing and new damage resulting from the CES. The damage items were overlaid with geospatial data such as Light Detection and Ranging (LiDAR)-derived ground deformation data, Liquefaction Resistance Index data and seismic ground motion data (Peak Ground acceleration and Peak Ground Velocity) to identify potential relationships between these parameters and pipe performance. Initial assessment outcomes for the pipe selection revealed that main pipe joints and lateral connections were more vulnerable than the pipe body during a seismic event. Smaller diameter pipes may also be more vulnerable than larger pipes during a seismic event. Obvious differential ground movement resulted in increased local damage observations in many cases, however this was not obvious for all pipes. Pipes with older installation ages exhibited more overall damage prior to a seismic event, which is likely attributable to increased chemical and biological deterioration. However, no evidence was found relating pipe age to performance during a seismic event. No evidence was found linking levels of pre-CES damage in a pipe with subsequent seismic performance, and seismic performance with liquefaction resistance or magnitude of seismic ground motion. The results reported are of limited application due to the small demonstration sample size, but reveal the additional level of detail and insight possible using the method presented in this thesis over existing assessment methods, especially in relation to high resolution variations along the length of the pipe such as localised ground deformations evidenced by LiDAR. The results may be improved by studying a larger and more diverse sample pool, automating data collection and input processes in order to improve efficiency and consider additional input such as pipe dip and cumulative damage over a large distance. The method is dependent on comprehensive and accurate pre-event CCTV assessments and LIDAR data so that post-event data could be compared. It is proposed that local territorial authorities should prioritise acquiring this information as a first important step towards improving the seismic resilience of a gravity waste water pipe network.
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.
Beach ridge stratigraphy can provide an important record of both sustained coastal progradation and responses to events such as extreme storms, as well as evidence of earthquake induced sediment pulses. This study is a stratigraphic investigation of the late Holocene mixed sand gravel (MSG) beach ridge plain on the Canterbury coast, New Zealand. The subsurface was imaged along a 370 m shore-normal transect using 100 and 200 MHz ground penetrating radar (GPR) antennae, and cored to sample sediment textures. Results show that, seaward of a back-barrier lagoon, the Pegasus Bay beach ridge plain prograded almost uniformly, under conditions of relatively stable sea level. Nearshore sediment supply appears to have created a sustained sediment surplus, perhaps as a result of post-seismic sediment pulses, resulting in a flat, morphologically featureless beach ridge plain. Evidence of a high magnitude storm provides an exception, with an estimated event return period in excess of 100 years. Evidence from the GPR sequence combined with modern process observations from MSG beaches indicates that a paleo storm initially created a washover fan into the back-barrier lagoon, with a large amount of sediment simultaneously moved off the beach face into the nearshore. This erosion event resulted in a topographic depression still evident today. In the subsequent recovery period, sediment was reworked by swash onto the beach as a sequence of berm deposit laminations, creating an elevated beach ridge that also has a modern-day topographic signature. As sediment supply returned to normal, and under conditions of falling sea level, a beach ridge progradation sequence accumulated seaward of the storm feature out to the modern-day beach as a large flat, uniform progradation plain. This study highlights the importance of extreme storm events and earthquake pulses on MSG coastlines in triggering high volume beach ridge formation during the subsequent recovery period.
