Heathcote Valley school strong motion station (HVSC) consistently recorded ground motions with higher intensities than nearby stations during the 2010-2011 Canterbury earthquakes. For example, as shown in Figure 1, for the 22 February 2011 Christchurch earthquake, peak ground acceleration at HVSC reached 1.4 g (horizontal) and 2 g (vertical), the largest ever recorded in New Zealand. Strong amplification of ground motions is expected at Heathcote Valley due to: 1) the high impedance contrast at the soil-rock interface, and 2) the interference of incident and surface waves within the valley. However, both conventional empirical ground motion prediction equations (GMPE) and the physics-based large scale ground motions simulations (with empirical site response) are ineffective in predicting such amplification due to their respective inherent limitations.
The 2010-2011 Canterbury earthquake sequence was extremely damaging to structures in Christchurch and continues to have a large economic and social impact on the city and surrounding regions. In addition to strong ground shaking (Bradley and Cubrinovski 2011 SRL; Bradley 2012 SDEE), extensive liquefaction was observed, particularly in the 4 September 2010 Darfield earthquake and the 22 February 2011 Christchurch earthquake (Cubrinovski et al. 2010 BNZSEE; 2011 SRL). Large observed vertical ground motion amplitudes were recorded in the events in this sequence, with vertical peak ground accelerations of over 2.2g being observed at the Heathcote Valley Primary School during the Christchurch earthquake, and numerous other vertical motions exceeding 1.0g (Bradley and Cubrinovski 2011 SRL; Bradley 2012 SDEE; Fry et al 2011 SRL). Vertical peak ground accelerations of over 1.2g were observed in the Darfield earthquake.
Karen Brown reports from Chch, where there seems to be no end in sight for the mental health problems caused by the earthquakes
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.
This poster discusses several possible approaches by which the nonlinear response of surficial soils can be explicitly modelled in physics-based ground motion simulations, focusing on the relative advantages and limitations of the various methodologies. These methods include fully-coupled 3D simulation models that directly allow soil nonlinearity in surficial soils, the domain reduction method for decomposing the physical domain into multiple subdomains for separate simulation, conventional site response analysis uncoupled from the simulations, and finally, the use of simple empirically based site amplification factors We provide the methodology for an ongoing study to explicitly incorporate soil nonlinearity into hybrid broadband simulations of the 2010-2011 Canterbury, New Zealand earthquakes.
This poster presents work to date on ground motion simulation validation and inversion for the Canterbury, New Zealand region. Recent developments have focused on the collection of different earthquake sources and the verification of the SPECFEM3D software package in forward and inverse simulations. SPECFEM3D is an open source software package which simulates seismic wave propagation and performs adjoint tomography based upon the spectral-element method. Figure 2: Fence diagrams of shear wave velocities highlighting the salient features of the (a) 1D Canterbury velocity model, and (b) 3D Canterbury velocity model. Figure 5: Seismic sources and strong motion stations in the South Island of New Zealand, and corresponding ray paths of observed ground motions. Figure 3: Domain used for the 19th October 2010 Mw 4.8 case study event including the location of the seismic source and strong motion stations. By understanding the predictive and inversion capabilities of SPECFEM3D, the current 3D Canterbury Velocity Model can be iteratively improved to better predict the observed ground motions. This is achieved by minimizing the misfit between observed and simulated ground motions using the built-in optimization algorithm. Figure 1 shows the Canterbury Velocity Model domain considered including the locations of small-to-moderate Mw events [3-4.5], strong motion stations, and ray paths of observed ground motions. The area covered by the ray paths essentially indicates the area of the model which will be most affected by the waveform inversion. The seismic sources used in the ground motion simulations are centroid moment tensor solutions obtained from GeoNet. All earthquake ruptures are modelled as point sources with a Gaussian source time function. The minimum Mw limit is enforced to ensure good signal-to-noise ratio and well constrained source parameters. The maximum Mw limit is enforced to ensure the point source approximation is valid and to minimize off-fault nonlinear effects.
In this paper, we perform hybrid broadband (0-10 Hz) ground motion simulations for the ten most significant events (Mw 4.7-7.1) in the 2010-2011 Canterbury earthquake sequence. Taking advantage of having repeated recordings at same stations, we validate our simulations using both recordings and an empirically-developed ground motion prediction equation (GMPE). The simulation clearly captures the sedimentary basin amplification and the rupture directivity effects. Quantitative comparisons of the simulations with both recordings and the GMPE, as well as analyses of the total residuals (indicating model bias) show that simulations perform better than the empirical GMPE, especially for long period. To scrutinize the ground motion variability, we partitioned the total residuals into different components. The total residual appears to be unbiased, and the use of a 3D velocity structure reduces the long period systematic bias particularly for stations located close to the Banks Peninsula volcanic area.
A view that was impossible six years ago. There was a group of buildings including another high-rise in the vacant area in the middle ground.
The first feeling that strikes everyone on coming to New Zealand is its intense want of animal life. Mountains, plains, rivers, – mere features without a soul; for you can hardly dignify the miserable ground lark, the wailing weka, or … Continue reading →
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.
essential systems upon which the well-being and functioning of societies depend. They deliver a service or a good to the population using a network, a combination of spatially-distributed links and nodes. As they are interconnected, network elements’ functionality is also interdependent. In case of a failure of one component, many others could be momentarily brought out-of-service. Further problems arise for buried infrastructure when it comes to buried infrastructure in earthquake and liquefaction-prone areas for the following reasons: • Technically more demanding inspections than those required for surface horizontal infrastructure • Infrastructure subject to both permanent ground displacement and transient ground deformation • Increase in network maintenance costs (i.e. deterioration due to ageing material and seismic hazard) These challenges suggest careful studies on network resilience will yield significant benefits. For these reasons, the potable water network of Christchurch city (Figure 1) has been selected for its well-characterized topology and its extensive repair dataset.
This project was initiated by ENGEO Limited and KiwiRail Holdings Limited to assess the stability of Slovens Creek Viaduct (specifically its western abutment) and a 3km section of rail corridor between Slovens Creek Viaduct and Avoca on the Midland Line (MDL). Commonly known as the scenic TranzAlpine rail journey (through Arthurs Pass National Park) the MDL connects Greymouth to Christchurch via Rolleston, where the MDL meets the Main South Line into Christchurch. The project area is approximately 40km southeast of Arthurs Pass Township, in the eastern extension of the Castle Hill Basin which is part of the Waimakariri Catchment and Canterbury Foothills. The field area is underlain by Rakaia Terrane, which is part of the Torlesse Composite Terrane forming the basement rock unit for the field area. Cretaceous-Tertiary rocks of the Castle Hill Basin overlie the basement strata and record a transgression-regression sequence, as well as mid-Oligocene submarine volcanism. The stratigraphic sequence in the Castle Hill Basin, and its eastern extension to Avoca, comprises two formations of the Eyre group, the older Broken River Formation and the younger Iron Creek Formation. Deep marine Porter Group limestones, marls, and tuffs of Oligocene age succeed the Iron Creek Formation of the Eyre Group, and probably records the maximum of the transgression. The Enys Formation lies disconformably on the Porter Group and is overlain unconformably by Late Pleistocene glacifluvial and glacial deposits. The Tertiary strata in the Slovens-Avoca rail corridor are weak, and the clay-rich tuff derived from mid-Oligocene volcanism is particularly prone to slaking. Extensive mapping carried out for this project has identified that some 90 percent of the surface along the length of the Slovens-Avoca corridor has been subject to mass movement. The landslides of the Slovens-Avoca rail corridor are clearly younger than the Last Glaciation, and Slovens Creek has been downcutting, with associated faulting and uplift, to form the present day geomorphology of the rail corridor. Deep-seated landslides in the rail corridor extend to Slovens Creek, locally deflecting the stream course, and a generic ground failure model for the rail corridor has been developed. Exploratory geotechnical investigations, including core drilling, installation of an inclinometer and a piezometer, enabled the construction of a simple ground model and cross section for the Slovens Creek Viaduct western abutment. Limit-equilibrium and pseudo-static slope stability analyses using both circular and block critical slip surface search methods were applied to the ground model for the western abutment of Slovens Creek Viaduct. Piezometric and strength data obtained during laboratory testing of core material have been used to constrain the western abutment stability assessment for one representative section line (C-C’). Prior to pseudo-static sensitivity analyses peak ground acceleration (PGA) for various Ultimate Limit State (ULS) design return periods, defined by an equation given in NZS1170.5:2004, were calculated and have been used as a calibration technique to find and compare specific PGA values for pseudo-static analyses in the Slovens Creek Viaduct area. The main purpose has been to provide an indication of how railway infrastructure could be affected by seismic events of various return periods defined by ULS design standards for the area. Limit equilibrium circular slip surface search methods, both grid search and auto refine search, indicated the slope is stable with a FoS greater than 1.0 returned from each, although one particular surface returned the lowest FoS in each. This surface is in the lower portion of the slope, adjacent to Slovens Stream and northeast of the MDL. As expected, pseudo-static analyses returned a lower FoS overall when compared to limit equilibrium analyses. The PGA analyses suggest that partial ground failure at the Slovens Creek Viaduct western abutment could occur in a 1 in 25-year return period event within materials on the slower slope beyond the immediate rail corridor. A ULS (1 in 500-year) event in the Slovens Creek Viaduct area would likely produce a PGA of ~0.9g, and the effects on the western abutment and rail infrastructure would most likely be catastrophic. Observed ground conditions for the western abutment of the Slovens Creek Viaduct suggest there is no movement within the landslide at depth within the monitoring timeframe of this project (22 May 2015 – 4 August 2015). Slope stability monitoring is recommended to be continued in two parts: (1) the inclinometer in BH1 is to be monitored on a six monthly basis for one year following completion of this thesis, and then annually unless ground movements become evident; and (2) surface movement monitoring should be installed using a fixed datum on the stable eastern abutment. Long-term stability management strategies for the Slovens Creek Viaduct western abutment are dependent upon future observed changes and ongoing monitoring. Hazard and risk assessment using the KiwiRail Qualitative Risk Assessment Framework (QRA) is recommended, and if slope stability becomes problematic for operation of the Midland Line consideration should be given to deep slope drainage. In the event of a large magnitude or high PGA earthquake all monitoring should be reviewed.
The 2015 New Zealand strong-motion database provides a wealth of new strong motion data for engineering applications. An important component of this database is the compilation of new site metadata, describing the soil conditions and site response at GeoNet strong motion stations. We have assessed and compiled four key site parameters for the ~460 GeoNet stations that recorded significant historical ground motions. Parameters include: site classification (NZS1170.5), Vs30, fundamental site period (Tsite) and depth to bedrock (Z1.0, i.e. depth to material with Vs > 1000 m/s). In addition, we have assigned a quality estimate (Quality 1 – 3) to these parameters to provide a qualitative estimate of the uncertainty. New highquality Tsite estimates have largely been obtained from newly available HVSR amplification curves and spectral ratios from inversion of regional strong motion data that has been reconciled with available geological information. Good quality Vs30 estimates, typically in urban centres, have also been incorporated following recent studies. Where site-specific measurements of Vs30 are not available, Vs30 is estimated based on surface geology following national Vs30 maps. New Z1.0 values have been provided from 3D subsurface models for Canterbury and Wellington. This database will be used in efforts to guide development and testing of new and existing ground motion prediction models in New Zealand. In particular, it will allow reexamination of the most important site parameters that control and predict site response in a New Zealand setting. Furthermore, it can be used to provide information about suitable rock reference sites for seismological research, and as a guide to site-specific references in the literature. We discuss compilation of the database, preliminary insights so far, and future directions.
Liquefaction-induced lateral spreading during earthquakes poses a significant hazard to the built environment, as observed in Christchurch during the 2010 to 2011 Canterbury Earthquake Sequence (CES). It is critical that geotechnical earthquake engineers are able to adequately predict both the spatial extent of lateral spreads and magnitudes of associated ground movements for design purposes. Published empirical and semi-empirical models for predicting lateral spread displacements have been shown to vary by a factor of <0.5 to >2 from those measured in parts of Christchurch during CES. Comprehensive post- CES lateral spreading studies have clearly indicated that the spatial distribution of the horizontal displacements and extent of lateral spreading along the Avon River in eastern Christchurch were strongly influenced by geologic, stratigraphic and topographic features.
The 2010-2011 Canterbury earthquake sequence, and the resulting extensive data sets on damaged buildings that have been collected, provide a unique opportunity to exercise and evaluate previously published seismic performance assessment procedures. This poster provides an overview of the authors’ methodology to perform evaluations with two such assessment procedures, namely the P-58 guidelines and the REDi Rating System. P-58, produced by the Federal Emergency Management Agency (FEMA) in the United States, aims to facilitate risk assessment and decision-making by quantifying earthquake ground shaking, structural demands, component damage and resulting consequences in a logical framework. The REDi framework, developed by the engineering firm ARUP, aids stakeholders in implementing resilience-based earthquake design. Preliminary results from the evaluations are presented. These have the potential to provide insights on the ability of the assessment procedures to predict impacts using “real-world” data. However, further work remains to critically analyse these results and to broaden the scope of buildings studied and of impacts predicted.
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.
With the occurrence of natural disasters on the increase, major cities around the world face the potential of complete loss of infrastructure due to design guidelines that do not consider resilience in the design. With the February 22nd, 2011 earthquake in Christchurch, being the largest insured event, lessons learnt from the rebuild will be vital for the preparation of future disasters. Therefore the objective of this research is to understand the financial implications of the changes to the waste water design guidelines used throughout the five year rebuild programme of works. The research includes a study of the SCIRT alliance model selected for the delivery that is flexible enough to handle changes in the design with minimal impact on the direct cost of the rebuild works. The study further includes the analysis and compares the impact of the three different guidelines on maintenance and replacement cost over the waste water pipe asset life. The research concludes that with the varying ground conditions in Christchurch and also the wide variety of materials in use in the waste water network up to the start of the CES, the rebuild was not a ‘one size fits all’ approach.
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.
Social media have changed disaster response and recovery in the way people inform themselves, provide community support and make sense of unfolding and past events online. During the Canterbury earthquakes of 2010 and 2011 social media platforms such as Facebook and Twitter became part of the story of the quakes in the region, as well as a basis for ongoing public engagement during the rebuild efforts in Christchurch. While a variety of research has been conducted on the use of social media in disaster situations (Bruns & Burgess, 2012; Potts, Seitzinger, Jones, & Harrison, 2011; Shklovski, Palen, & Sutton, 2008), studies about their uses in long-term disaster recovery and across different platforms are underrepresented. This research analyses networked practices of sensemaking around the Canterbury earthquakes over the course of disaster response, recovery and rebuild, focussing on Facebook and Twitter. Following a mixed methodological design data was gathered in interviews with people who started local Facebook pages, and through digital media methods of data collection and computational analysis of public Facebook pages and a historical Twitter dataset gathered around eight different earthquake-related events between 2010 and 2013. Data is further analysed through discursive and narrative tools of inquiry. This research sheds light on communication practices in the drawn-out process of disaster recovery on the ground in connecting different modes of discourse. Examining the ongoing negotiation of networked identities through technologically mediated social practices during Canterbury’s rebuild, the connection between online environments and the city of Christchurch, as a physical place, is unpacked. This research subsequently develops a new methodology to study social media platforms and provide new and detailed information on both the communication practices in issue-based online publics and the ongoing negotiation of the impact of the Canterbury earthquakes through networked digital means.
Sewerage systems convey sewage, or wastewater, from residential or commercial buildings through complex reticulation networks to treatment plants. During seismic events both transient ground motion and permanent ground deformation can induce physical damage to sewerage system components, limiting or impeding the operability of the whole system. The malfunction of municipal sewerage systems can result in the pollution of nearby waterways through discharge of untreated sewage, pose a public health threat by preventing the use of appropriate sanitation facilities, and cause serious inconvenience for rescuers and residents. Christchurch, the second largest city in New Zealand, was seriously affected by the Canterbury Earthquake Sequence (CES) in 2010-2011. The CES imposed widespread damage to the Christchurch sewerage system (CSS), causing a significant loss of functionality and serviceability to the system. The Christchurch City Council (CCC) relied heavily on temporary sewerage services for several months following the CES. The temporary services were supported by use of chemical and portable toilets to supplement the damaged wastewater system. The rebuild delivery agency -Stronger Christchurch Infrastructure Rebuild Team (SCIRT) was created to be responsible for repair of 85 % of the damaged horizontal infrastructure (i.e., water, wastewater, stormwater systems, and roads) in Christchurch. Numerous initiatives to create platforms/tools aiming to, on the one hand, support the understanding, management and mitigation of seismic risk for infrastructure prior to disasters, and on the other hand, to support the decision-making for post-disaster reconstruction and recovery, have been promoted worldwide. Despite this, the CES in New Zealand highlighted that none of the existing platforms/tools are either accessible and/or readable or usable by emergency managers and decision makers for restoring the CSS. Furthermore, the majority of existing tools have a sole focus on the engineering perspective, while the holistic process of formulating recovery decisions is based on system-wide approach, where a variety of factors in addition to technical considerations are involved. Lastly, there is a paucity of studies focused on the tools and frameworks for supporting decision-making specifically on sewerage system restoration after earthquakes. This thesis develops a decision support framework for sewerage pipe and system restoration after earthquakes, building on the experience and learning of the organisations involved in recovering the CSS following the CES in 2010-2011. The proposed decision support framework includes three modules: 1) Physical Damage Module (PDM); 2) Functional Impact Module (FIM); 3) Pipeline Restoration Module (PRM). The PDM provides seismic fragility matrices and functions for sewer gravity and pressure pipelines for predicting earthquake-induced physical damage, categorised by pipe materials and liquefaction zones. The FIM demonstrates a set of performance indicators that are categorised in five domains: structural, hydraulic, environmental, social and economic domains. These performance indicators are used to assess loss of wastewater system service and the induced functional impacts in three different phases: emergency response, short-term recovery and long-term restoration. Based on the knowledge of the physical and functional status-quo of the sewerage systems post-earthquake captured through the PDM and FIM, the PRM estimates restoration time of sewer networks by use of restoration models developed using a Random Forest technique and graphically represented in terms of restoration curves. The development of a decision support framework for sewer recovery after earthquakes enables decision makers to assess physical damage, evaluate functional impacts relating to hydraulic, environmental, structural, economic and social contexts, and to predict restoration time of sewerage systems. Furthermore, the decision support framework can be potentially employed to underpin system maintenance and upgrade by guiding system rehabilitation and to monitor system behaviours during business-as-usual time. In conjunction with expert judgement and best practices, this framework can be moreover applied to assist asset managers in targeting the inclusion of system resilience as part of asset maintenance programmes.
Liquefaction features and the geologic environment in which they formed were carefully studied at two sites near Lincoln in southwest Christchurch. We undertook geomorphic mapping, excavated trenches, and obtained hand cores in areas with surficial evidence for liquefaction and areas where no surficial evidence for liquefaction was present at two sites (Hardwick and Marchand). The liquefaction features identified include (1) sand blows (singular and aligned along linear fissures), (2) blisters or injections of subhorizontal dikes into the topsoil, (3) dikes related to the blows and blisters, and (4) a collapse structure. The spatial distribution of these surface liquefaction features correlates strongly with the ridges of scroll bars in meander settings. In addition, we discovered paleoliquefaction features, including several dikes and a sand blow, in excavations at the sites of modern liquefaction. The paleoliquefaction event at the Hardwick site is dated at A.D. 908-1336, and the one at the Marchand site is dated at A.D. 1017-1840 (95% confidence intervals of probability density functions obtained by Bayesian analysis). If both events are the same, given proximity of the sites, the time of the event is A.D. 1019-1337. If they are not, the one at the Marchand site could have been much younger. Taking into account a preliminary liquefaction-triggering threshold of equivalent peak ground acceleration for an Mw 7.5 event (PGA7:5) of 0:07g, existing magnitude-bounded relations for paleoliquefaction, and the timing of the paleoearthquakes and the potential PGA7:5 estimated for regional faults, we propose that the Porters Pass fault, Alpine fault, or the subduction zone faults are the most likely sources that could have triggered liquefaction at the study sites. There are other nearby regional faults that may have been the source, but there is no paleoseismic data with which to make the temporal link.
One of the less understood geotechnical responses to the cyclic loading from the MW6.2 Christchurch Earthquake, on the 22nd of February 2011, is the fissuring in the loessial soil-mantled, footslope positions of the north-facing valleys of the Port Hills. The fissures are characterized by mostly horizontal offset (≤500mm), with minor vertical displacement (≤300mm), and they extend along both sides of valleys for several hundred metres in an approximately contour-parallel orientation. The fissure traces correspond to extensional features mapped in other studies. Previous studies have suggested that the fissures are the headscarps of incipient landslides, but the surface and subsurface features are not typical of landslide movement. Whilst there are some features that correlate with landslide movement, there are many features that contradict the landslide movement hypothesis. Of critical importance to this investigation was the fact that there are no landslide flanks, there has been no basal shear surface found, there is little deformation in the so-called ‘landslide body’, and there have been no recorded zones of low shear strength in the soil deposit that are indicative of a basal shear surface. This thesis is a detailed geotechnical study on the fissures along part of Ramahana Road in the Hillsborough Valley, Christchurch. Shallow and deep investigation methods found that the predominant soil is loess-colluvium, to depths of ~20m, and this soil has variable geotechnical characteristics depending on the layer sampled. The factor that has the most influence on shear strength was found to be the moisture content. Direct shear-box testing of disturbed, recompacted loess-colluvium found that the soil had a cohesion of 35-65kPa and a friction angle of 38-43° when the soil moisture content was at 8-10%. However when the moisture content was at 19-20% the soil’s cohesion decreased to 3-5kPa and its friction angle decreased to 33-38°, this moisture content is at or slightly above the plastic limit. An electrical resistivity geophysical survey was conducted perpendicular to multiple fissure traces and through the compressional zone at 17 Ramahana Road. The electrical resistivity line found that there was an area of high resistivity at the toe of the slope, and an area of high conductivity downslope of this and at greater depths. This area correlated to the compressional zone recorded by previous studies. Moisture content testing of the soil in these locations showed that the soil in the resistive area was relatively dry (9%) compared to the surrounding soil (13%), whilst the soil in the conductive area was relatively wet (22%)compared to the surrounding soil (19%). Density tests of the soil in the compressional zone recorded that the resistive area had a higher dry density than the surrounding soil (~1790 kg/m3 compared to ~1650 kg/m3). New springs arose downslope of the compressional zone contemporaneously with the fissures, and it is interpreted that these have arisen from increased hydraulic head in the Banks Peninsula bedrock aquifer system, and earthquake induced-bedrock fracturing. A test pit was dug across an infilled fissure trace at 17 Ramahana Road to a depth of 3m. The fissure trace had an aperture of 450-470mm at the ground surface, but it gradually lost aperture with depth until 2.0-2.1m where it became a segmented fissure trace with 1-2mm aperture. A mixed-colluvium layer was intercepted by the fissure trace at 2.4m depth, and there was no observable vertical offset of this layer. The fissure trace was at an angle of 78° at the ground surface, but it also flattened with depth, which gave it a slightly curved appearance. The fissure trace was at an assumed angle of 40-50° near the base of the test pit. Rotational slide, translational slide and lateral spread landslide movement types were compared and contrasted as possibilities for landslide movement types, whilst an alternative hypothesis was offered that the fissures are tensile failures with a quasi-toppling motion involving a cohesive block of loessial soil moving outwards from the slope, with an accommodating compressional strain in the lower less cohesive soil. The mechanisms behind this movement are suggested to be the horizontal earthquake inertia forces from the Christchurch Earthquake, the static shear stress of the slope, and bedrock uplift of the Port Hills in relation to the subsidence of the Christchurch city flatlands. Extremely high PGA is considered to be a prerequisite to the fissure trace development, and these can only be induced in the Hillsborough Valley from a Port Hills Fault rupture, which has a recurrence interval of ~10,000 years. The current understanding of how the loess-colluvium soil would behave under cyclic loading is limited, and the mechanisms behind the suggested movement type are not completely understood. Further research is needed to confirm the proposed mechanism of the fissure traces. Laboratory tests such as the cyclic triaxial and cyclic shear test would be beneficial in future research to quantitatively test how the soil behaves under cyclic loading at various moisture contents and clay contents, and centrifuge experiments would be of great use to qualitatively test the suggested mode of movement in the loessial soil.
Following the 22nd February 2011, Mw 6.2 earthquake located along a previously unknown fault beneath the Port Hills of Christchurch, surface cracking was identified in contour parallel locations within fill material at Quarry Road on the lower slopes of Mount Pleasant. GNS Science, in the role of advisor to the Christchurch City Council, concluded that these cracks were a part of a potential rotational mass movement (named zone 11A) within the fill and airfall loess material present. However, a lack of field evidence for slope instability and an absence of laboratory geotechnical data on which slope stability analysis was based, suggested this conclusion is potentially incorrect. It was hypothesised that ground cracking was in fact due to earthquake shaking, and not mass movement within the slope, thus forming the basis of this study. Three soil units were identified during surface and subsurface investigations at Quarry Road: fill derived from quarry operations in the adjacent St. Andrews Quarry (between 1893 and 1913), a buried topsoil, and underlying in-situ airfall loess. The fill material was identified by the presence of organic-rich topsoil “clods” that were irregular in both size (∼10 – 200 mm) and shape, with variable thicknesses of 1 – 10 m. Maximum thickness, as indicated by drill holes and geophysical survey lines, was identified below 6 Quarry Road and 7 The Brae where it is thought to infill a pre-existing gully formed in the underlying airfall loess. Bearing strength of the fill consistently exceeded 300 kPa ultimate below ∼500 mm depth. The buried topsoil was 200 – 300 mm thick, and normally displayed a lower bearing strength when encountered, but not below 300 kPa ultimate (3 – 11 blows per 100mm or ≥100 kPa allowable). In-situ airfall loess stood vertically in outcrop due to its characteristic high dry strength and also showed Scala penetrometer values of 6 – 20+ blows per 100 mm (450 – ≥1000 kPa ultimate). All soils were described as being moist to dry during subsurface investigations, with no groundwater table identified during any investigation into volcanic bedrock. In-situ moisture contents were established using bulk disturbed samples from hand augers and test pitting. Average moisture contents were low at 9% within the fill, 11 % within the buried topsoil, and 8% within the airfall loess: all were below the associated average plastic limit of 17, 15, and 16, respectively, determined during Atterberg limit analysis. Particle size distributions, identified using the sieve and pipette method, were similar between the three soil units with 11 – 20 % clay, 62 – 78 % silt, and 11 – 20 % fine sand. Using these results and the NZGS soil classification, the loess derived fill and in-situ airfall loess are termed SILT with some clay and sand, and the buried topsoil is SILT with minor clay and sand. Dispersivity of the units was found using the Emerson crumb test, which established that the fill can be non- to completely dispersive (score 0 – 4). The buried topsoil was always non-dispersive (score 0), and airfall loess completely dispersive (score 4). Values for cohesion (c) and internal friction angle (φ) of the three soil units were established using the direct shear box at field moisture contents. Results showed all soil units had high shear strengths at the moisture contents tested (c = 18 – 24 kPa and φ = 42 – 50°), with samples behaving in a brittle fashion. Moisture content was artificially increased to 16% within the buried topsoil, which reduced the shear strength (c = 10 kPa, φ = 18°) and allowed it to behave plastically. Observational information indicating stability at Quarry Road included: shallow, discontinuous, cracks that do not display vertical offset; no scarp features or compressional zones typical of landsliding; no tilted or deformed structures; no movement in inclinometers; no basal shear zone identified in logged core to 20 m depth; low field moisture contents; no groundwater table; and high soil strength using Scala penetrometers. Limit equilibrium analysis of the slope was conducted using Rocscience software Slide 5.0 to verify the slope stability identified by observational methods. Friction, cohesion, and density values determined during laboratory were input into the two slope models investigated. Results gave minimum static factor of safety values for translational (along buried topsoil) and rotational (in the fill) slides of 2.4 – 4.2. Sensitivity of the slope to reduced shear strength parameters was analysed using c = 10 kPa and φ = 18° for the translational buried topsoil plane, and a cohesion of 0 kPa within the fill for the rotational plane. The only situation that gave a factor of safety <1.0 was in nonengineered fill at 0.5 m depth. Pseudostatic analysis based on previous peak ground acceleration (PGA) values for the Canterbury Earthquake Sequence, and predicted PGAs for future Alpine Fault and Hope Fault earthquakes established minimum factor of safety values between 1.2 and 3.3. Yield acceleration PGAs were computed to be between 0.8g and 1.6g. Based on all information gathered, the cracking at Quarry Road is considered to be shallow deformation in response to earthquake shaking, and not due to deep-seated landsliding. It is recommended that the currently bare site be managed by smoothing the land, installing contour drainage, and bioremediation of the surface soils to reduce surface water infiltration and runoff. Extensive earthworks, including removal of the fill, are considered unnecessary. Any future replacement of housing would be subject to site-specific investigations, and careful foundation design based on those results.
The Screw Driving Sounding (SDS) method developed in Japan is a relatively new insitu testing technique to characterise soft shallow sites, typically those required for residential house construction. An SDS machine drills a rod into the ground in several loading steps while the rod is continuously rotated. Several parameters, such as torque, load and speed of penetration, are recorded at every rotation of the rod. The SDS method has been introduced in New Zealand, and the results of its application for characterising local sites are discussed in this study. A total of 164 SDS tests were conducted in Christchurch, Wellington and Auckland to validate/adjust the methodologies originally developed based on the Japanese practice. Most of the tests were conducted at sites where cone penetration tests (CPT), standard penetration tests (SPT) and borehole logs were available; the comparison of SDS results with existing information showed that the SDS method has great potential as an in-situ testing method for classifying the soils. By compiling the SDS data from 3 different cities and comparing them with the borehole logs, a soil classification chart was generated for identifying the soil type based on SDS parameters. Also, a correlation between fines content and SDS parameters was developed and a procedure for estimating angle of internal friction of sand using SDS parameters was investigated. Furthermore, a correlation was made between the tip resistance of the CPT and the SDS data for different percentages of fines content. The relationship between the SPT N value and a SDS parameter was also proposed. This thesis also presents a methodology for identifying the liquefiable layers of soil using SDS data. SDS tests were performed in both liquefied and non-liquefied areas in Christchurch to find a representative parameter and relationship for predicting the liquefaction potential of soil. Plots were drawn of the cyclic shear stress ratios (CSR) induced by the earthquakes and the corresponding energy of penetration during SDS tests. By identifying liquefied or unliquefied layers using three different popular CPT-based methods, boundary lines corresponding to the various probabilities of liquefaction happening were developed for different ranges of fines contents using logistic regression analysis, these could then be used for estimating the liquefaction potential of soil directly from the SDS data. Finally, the drilling process involved in screw driving sounding was simulated using Abaqus software. Analysis results proved that the model successfully captured the drilling process of the SDS machine in sand. In addition, a chart to predict peak friction angles of sandy sites based on measured SDS parameters for various vertical effective stresses was formulated. As a simple, fast and economical test, the SDS method can be a reliable alternative insitu test for soil and site characterisation, especially for residential house construction.
This report provides an initial overview and gap analysis of the multi-hazards interactions that might affect fluvial and pluvial flooding (FPF) hazard in the Ōpāwaho Heathcote catchment. As per the terms of reference, this report focuses on a one-way analysis of the potential effects of multi-hazards on FPF hazard, as opposed to a more complex multi-way analysis of interactions between all hazards. We examined the relationship between FPF hazard and hazards associated with the phenomena of tsunamis; coastal erosion; coastal inundation; groundwater; earthquakes; and mass movements. Tsunamis: Modelling research indicates the worst-case tsunami scenarios potentially affecting the Ōpāwaho Heathcote catchment are far field. Under low probability, high impact tsunami scenarios waves could travel into Pegasus Bay and the Avon-Heathcote Estuary Ihutai, reaching the mouth and lower reaches of the Heathcote catchment and river, potentially inundating and eroding shorelines in sub-catchments 1 to 5, and temporarily blocking fluvial drainage more extensively. Any flooding infrastructure or management actions implemented in the area of tsunami inundation would ideally be resilient to tsunami-induced inundation and erosion. Model results currently available are a first estimate of potential tsunami inundation under contemporary sea and land level conditions. In terms of future large tsunami events, these models likely underestimate effects in riverside sub-catchments, as well as effects under future sea level, shoreline and other conditions. Also of significance when considering different FPF management structures, it is important to be mindful that certain types of flood structures can ‘trap’ inundating water coming from ocean directions, leading to longer flood durations and salinization issues. Coastal erosion: Model predictions indicate that sub-catchments 1 to 3 could potentially be affected by coastal erosion by the timescale of 2065, with sub-catchments 1-6 predicted to be potentially affected by coastal erosion by the time scale of 2115. In addition, the predicted open coast effects of this hazard should not be ignored since any significant changes in the New Brighton Spit open coast would affect erosion rates and exposure of the landward estuary margins, including the shorelines of the Ōpāwaho Heathcote catchment. Any FPF flooding infrastructure or management activities planned for the potentially affected sub-catchments needs to recognise the possibility of coastal erosion, and to have a planned response to the predicted potential shoreline translation. Coastal inundation: Model predictions indicate coastal inundation hazards could potentially affect sub-catchments 1 to 8 by 2065, with a greater area and depth of inundation possible for these same sub-catchments by 2115. Low-lying areas of the Ōpāwaho Heathcote catchment and river channel that discharge into the estuary are highly vulnerable to coastal inundation since elevated ocean and estuary water levels can block the drainage of inland systems, compounding FPF hazards. Coastal inundation can overwhelm stormwater and other drainage network components, and render river dredging options ineffective at best, flood enhancing at worst. A distinction can be made between coastal inundation and coastal erosion in terms of the potential impacts on affected land and assets, including flood infrastructure, and the implications for acceptance, adaptation, mitigation, and/or modification options. That is, responding to inundation could include structural and/or building elevation solutions, since unlike erosion, inundation does not necessarily mean the loss of land. Groundwater: Groundwater levels are of significant but variable concern when examining flooding hazards and management options in the Ōpāwaho Heathcote catchment due to variability in soils, topographies, elevations and proximities to riverine and estuarine surface waterbodies. Much of the Canterbury Plains part of the Ōpāwaho Heathcote catchment has a water table that is at a median depth of <1m from the surface (with actual depth below surface varying seasonally, inter-annually and during extreme meteorological events), though the water table depth rapidly shifts to >6m below the surface in the upper Plains part of the catchment (sub-catchments 13 to 15). Parts of Waltham/Linwood (sub-catchments 5 & 6) and Spreydon (sub-catchment 10) have extensive areas with a particularly high water table, as do sub-catchments 18, 19 and 20 south of the river. In all of the sub-catchments where groundwater depth below surface is shallow, it is necessary to be mindful of cascading effects on liquefaction hazard during earthquake events, including earthquake-induced drainage network and stormwater infrastructure damage. In turn, subsidence induced by liquefaction and other earthquake processes during the CES directly affected groundwater depth below surface across large parts of the central Ōpāwaho Heathcote catchment. The estuary margin of the catchment also faces increasing future challenges with sea level rise, which has the potential to elevate groundwater levels in these areas, compounding existing liquefaction and other earthquake associated multi-hazards. Any increases in subsurface runoff due to drainage system, development or climate changes are also of concern for the loess covered hill slopes due to the potential to enhance mass movement hazards. Earthquakes: Earthquake associated vertical ground displacement and liquefaction have historically affected, or are in future predicted to affect, all Ōpāwaho Heathcote sub-catchments. During the CES, these phenomena induced a significant cascades of changes in the city’s drainage systems, including: extensive vertical displacement and liquefaction induced damage to stormwater ‘greyware’, reducing functionality of the stormwater system; damage to the wastewater system which temporarily lowered groundwater levels and increased stormwater drainage via the wastewater network on the one hand, creating a pollution multi-hazard for FPF on the other hand; liquefaction and vertical displacement induced river channel changes affected drainage capacities; subsidence induced losses in soakage and infiltration capacities; changes occurred in topographic drainage conductivity; estuary subsidence (mainly around the Ōtākaro Avon rivermouth) increased both FPF and coastal inundation hazards; estuary bed uplift (severe around the Ōpāwaho Heathcote margins), reduced tidal prisms and increased bed friction, producing an overall reduction the waterbody’s capacity to efficiently flush catchment floodwaters to sea; and changes in estuarine and riverine ecosystems. All such possible effects need to be considered when evaluating present and future capacities of the Ōpāwaho Heathcote catchment FPF management systems. These phenomena are particularly of concern in the Ōpāwaho Heathcote catchment since stormwater networks must deal with constraints imposed by stream and river channels (past and present), estuarine shorelines and complex hill topography. Mass movements: Mass movements are primarily a risk in the Port Hills areas of the Ōpāwaho Heathcote catchment (sub-catchments 1, 2, 7, 9, 11, 16, 21), though there are one or two small but susceptible areas on the banks of the Ōpāwaho Heathcote River. Mass movements in the form of rockfalls and debris flows occurred on the Port Hills during the CES, resulting in building damage, fatalities and evacuations. Evidence has also been found of earthquake-triggered tunnel gully collapsesin all Port Hill Valleys. Follow-on effects of these mass movements are likely to occur in major future FPF and other hazard events. Of note, elevated groundwater levels, coastal inundation, earthquakes (including liquefaction and other effects), and mass movement exhibit the most extensive levels of multi-hazard interaction with FPF hazard. Further, all of the analysed multi-hazard interactions except earthquakes were found to consistently produce increases in the FPF hazard. The implications of these analyses are that multihazard interactions generally enhance the FPF hazard in the Ōpāwaho Heathcote catchment. Hence, management plans which exclude adjustments for multi-hazard interactions are likely to underestimate the FPF hazard in numerous different ways. In conclusion, although only a one-way analysis of the potential effects of selected multi-hazards on FPF hazard, this review highlights that the Ōpāwaho Heathcote catchment is an inherently multi- hazard prone environment. The implications of the interactions and process linkages revealed in this report are that several significant multi-hazard influences and process interactions must be taken into account in order to design a resilient FPF hazard management strategy.
