Simulation-based PSHA for the Canterbury region (Cybershake v17.8)
Research papers, University of Canterbury Library
None
None
To identify key ground characteristics that led to different liquefaction manifestations during the Canterbury earthquakes
None
None
Objectives • To develop a system dynamics model of Christchurch post-quake reconstruction process that captures all the critical dynamics influencing its pathway • To investigate the implications of current rebuild pathway • To build a reconstruction module to be integrated in MERIT (Measuring the Economics of Resilient Infrastructure Tool)
Abstract This study provides a simplified methodology for pre-event data collection to support a faster and more accurate seismic loss estimation. Existing pre-event data collection frameworks are reviewed. Data gathered after the Canterbury earthquake sequences are analysed to evaluate the relative importance of different sources of building damage. Conclusions drawns are used to explore new approaches to conduct pre-event building assessment.
Background Liquefaction induced land damage has been identified in more than 13 notable New Zealand earthquakes within the past 150 years, as presented on the timeline below. Following the 2010-2011 Canterbury Earthquake Sequence (CES), the consequences of liquefaction were witnessed first-hand in the city of Christchurch and as a result the demand for understanding this phenomenon was heightened. Government, local councils, insurers and many other stakeholders are now looking to research and understand their exposure to this natural hazard.
The term resilience‘’is increasingly being used in a multitude of contexts. Seemingly the latest buzz‘’word, it can mean many things to many people, in many different situations. In a natural hazard context, the terms sustainable planning‘’, and resilience‘planning are now’being used, often interchangeably. This poster provides an overview of resilience and sustainability within a land use planning and natural hazard context, and discusses how they are interrelated in the situation of the earthquake impacted city of Christchurch, New Zealand.
he 2016 Building (Earthquake Prone Building) Amendment Act aims to improve the system for managing earthquake-prone buildings. The proposed changes to the Act were precipitated by the Canterbury earthquakes, and the need to improve the seismic safety of New Zealand’s building stock. However, the Act has significant ramifications for territorial authorities, organisations and individuals in small New Zealand towns, since assessing and repairing heritage buildings poses a major cost to districts with low populations and poor rental returns on commercial buildings.
Many buildings with relatively low damage from the 2010-2011 Canterbury were deemed uneconomic to repair and were replaced [1,2]. Factors that affected commercial building owners’ decisions to replace rather than repair, included capital availability, uncertainty with regards to regional recovery, local market conditions and ability to generate cash flow, and repair delays due to limited property access (cordon). This poster provides a framework for modeling decision-making in a case where repair is feasible but replacement might offer greater economic value – a situation not currently modeled in engineering risk analysis.
Unreinforced masonry churches in New Zealand, similarly to everywhere else in the word have proven to be highly vulnerable to earthquakes, because of their particular construction features. The Canterbury (New Zealand) earthquake sequence, 2010-2011 caused an invaluable loss of local architectural heritage and of churches, as regrettably, some of them were demolished instead of being repaired. It is critical for New Zealand to advance the data collection, research and understanding pertaining to the seismic performance and protection of church buildings, with the aim to:
Nel presente articolo si illustra una procedura per il processamento automatizzato di prove CPT, il calcolo di vari indici di liquefazione e la rappresentazione dei dati su mappa. La procedura è applicata al caso studio del terremoto di Christchurch, Nuova Zelanda, del 22 febbraio 2011 (magnitudo momento, Mw = 6.2). Dall’analisi spaziale dei risultati emerge una buona correlazione tra le mappe ottenute per l’indicatore degli effetti al suolo e i danni osservati (su terreni e strutture). Tuttavia, per confermare la validità di tale procedura, sarà necessario esaminare ulteriori casi studio nel mondo.
This study examines the performance of nonlinear total-stress wave-propagation site response analysis for modelling site effects in physics-based ground motion simulations of the 2010-2011 Canterbury, New Zealand earthquake sequence. This approach allows for explicit modeling of 3-dimensional ground motion phenomena at the regional scale, as well as detailed site effects and soil nonlinearity at the local scale. The approach is compared to a more commonly used empirical VS30 (30 m time-averaged shear wave velocity)-based method for computing site amplification as proposed by Graves and Pitarka (2010, 2015).
Motivation This poster aims to present fragility functions for pipelines buried in liquefaction-prone soils. Existing fragility models used to quantify losses can be based on old data or use complex metrics. Addressing these issues, the proposed functions are based on the Christchurch network and soil and utilizes the Canterbury earthquake sequence (CES) data, partially represented in Figure 1. Figure 1 (a) presents the pipe failure dataset, which describes the date, location and pipe on which failures occurred. Figure 1 (b) shows the simulated ground motion intensity median of the 22nd February 2011 earthquake. To develop the model, the network and soil characteristics have also been utilized.
This poster aims to present fragility functions for pipelines buried in liquefaction-prone soils. Existing fragility models used to quantify losses can be based on old data or use complex metrics. Addressing these issues, the proposed functions are based on the Christchurch network and soil and utilizes the Canterbury earthquake sequence (CES) data, partially represented in Figure 1. Figure 1 (a) presents the pipe failure dataset, which describes the date, location and pipe on which failures occurred. Figure 1 (b) shows the simulated ground motion intensity median of the 22nd February 2011 earthquake. To develop the model, the network and soil characteristics have also been utilized
Existing unreinforced masonry (URM) buildings are often composed of traditional construction techniques, with poor connections between walls and diaphragms that results in poor performance when subjected to seismic actions. In these cases the application of the common equivalent static procedure is not applicable because it is not possible to assure “box like” behaviour of the structure. In such conditions the ultimate strength of the structure relies on the behaviour of the macro-elements that compose the deformation mechanisms of the whole structure. These macroelements are a single or combination of structural elements of the structure which are bonded one to each other. The Canterbury earthquake sequence was taken as a reference to estimate the most commonly occurring collapse mechanisms found in New Zealand URM buildings in order to define the most appropriate macroelements.
The latest two great earthquake sequences; 2010- 2011 Canterbury Earthquake and 2016 Kaikoura Earthquake, necessitate a better understanding of the New Zealand seismic hazard condition for new building design and detailed assessment of existing buildings. It is important to note, however, that the New Zealand seismic hazard map in NZS 1170.5.2004 is generalised in effort to cover all of New Zealand and limited to a earthquake database prior to 2001. This is “common” that site-specific studies typically provide spectral accelerations different to those shown on the national map (Z values in NZS 1170.5:2004); and sometimes even lower. Moreover, Section 5.2 of Module 1 of the Earthquake Geotechnical Engineering Practice series provide the guidelines to perform site- specific studies.
SeisFinder is an open-source web service developed by QuakeCoRE and the University of Canterbury, focused on enabling the extraction of output data from computationally intensive earthquake resilience calculations. Currently, SeisFinder allows users to select historical or future events and retrieve ground motion simulation outputs for requested geographical locations. This data can be used as input for other resilience calculations, such as dynamic response history analysis. SeisFinder was developed using Django, a high-level python web framework, and uses a postgreSQL database. Because our large-scale computationally-intensive numerical ground motion simulations produce big data, the actual data is stored in file systems, while the metadata is stored in the database.
This paper presents on-going challenges in the present paradigm shift of earthquakeinduced ground motion prediction from empirical to physics-based simulation methods. The 2010-2011 Canterbury and 2016 Kaikoura earthquakes are used to illustrate the predictive potential of the different methods. On-going efforts on simulation validation and theoretical developments are then presented, as well as the demands associated with the need for explicit consideration of modelling uncertainties. Finally, discussion is also given to the tools and databases needed for the efficient utilization of simulated ground motions both in specific engineering projects as well as for near-real-time impact assessment.
Overview of SeisFinder SeisFinder is an open-source web service developed by QuakeCoRE and the University of Canterbury, focused on enabling the extraction of output data from computationally intensive earthquake resilience calculations. Currently, SeisFinder allows users to select historical or future events and retrieve ground motion simulation outputs for requested geographical locations. This data can be used as input for other resilience calculations, such as dynamic response history analysis. SeisFinder was developed using Django, a high-level python web framework, and uses a postgreSQL database. Because our large-scale computationally-intensive numerical ground motion simulations produce big data, the actual data is stored in file systems, while the metadata is stored in the database. The basic SeisFinder architecture is shown in Figure 1.
This work investigates the possibility of developing a non-contact, non-line of sight sensor to measure interstorey drift through simulation and experimental validation. • The method uses frequency-modulated continuous wave (FMCW) radar to measure displacement. This method is commonly in use in a number of modern applications, including aircraft altimeters and automotive parking sensors. • The technique avoids numerous problems found in contemporary structural health monitoring methods, namely integral drift errors and structural modification requirements. • The smallest achievable detection error in displacement was found to be as low as 0.26%, through simulated against the displacement response of a single degree of freedom structure subject to ground motion excitation. • This was verified during experimentation, when a corner-style reflector was placed on a shake table running ground motion data taken from the 4th September 2010 earthquake in Christchurch. These results confirmed the conclusions drawn from simulation.
Context of the project: On 4 September 2010, 22 February 2011, 13 June 2011 and 23 December 2011 Christchurch suffered major earthquakes and aftershocks (well over 10,000) that have left the central city in ruins and many of the eastern suburbs barely habitable even now. The earthquakes on 22 February caused catastrophic loss of life with 185 people killed. The toll this has taken on the residents of Christchurch has been considerable, not least of all for the significant psychological impact and disruption it has had on the children. As the process of rebuilding the city commenced, it became clear that the arts would play a key role in maintaining our quality of life during difficult times. For me, this started with the children and the most expressive of all the art forms – music.
The Canterbury Earthquake Sequence 2010-2011 (CES) induced widespread liquefaction in many parts of Christchurch city. Liquefaction was more commonly observed in the eastern suburbs and along the Avon River where the soils were characterised by thick sandy deposits with a shallow water table. On the other hand, suburbs to the north, west and south of the CBD (e.g. Riccarton, Papanui) exhibited less severe to no liquefaction. These soils were more commonly characterised by inter-layered liquefiable and non-liquefiable deposits. As part of a related large-scale study of the performance of Christchurch soils during the CES, detailed borehole data including CPT, Vs and Vp have been collected for 55 sites in Christchurch. For this subset of Christchurch sites, predictions of liquefaction triggering using the simplified method (Boulanger & Idriss, 2014) indicated that liquefaction was over-predicted for 94% of sites that did not manifest liquefaction during the CES, and under-predicted for 50% of sites that did manifest liquefaction. The focus of this study was to investigate these discrepancies between prediction and observation. To assess if these discrepancies were due to soil-layer interaction and to determine the effect that soil stratification has on the develop-ment of liquefaction and the system response of soil deposits.
1. Background and Objectives This poster presents results from ground motion simulations of small-to-moderate magnitude (3.5≤Mw≤5.0) earthquake events in the Canterbury, New Zealand region using the Graves and Pitarka (2010,2015) methodology. Subsequent investigation of systematic ground motion effects highlights the prediction bias in the simulations which are also benchmarked against empirical ground motion models (e.g. Bradley (2013)). In this study, 144 earthquake ruptures, modelled as point sources, are considered with 1924 quality-assured ground motions recorded across 45 strong motion stations throughout the Canterbury region, as shown in Figure 1. The majority of sources are Mw≥4.0 and have centroid depth (CD) 10km or shallower. Earthquake source descriptions were obtained from the GeoNet New Zealand earthquake catalogue. The ground motion simulations were performed within a computational domain of 140km x 120km x 46km with a finite difference grid spacing of 0.1km. The low-frequency (LF) simulations utilize the 3D Canterbury Velocity Model while the high-frequency (HF) simulations utilize a generic regional 1D velocity model. In the LF simulations, a minimum shear wave velocity of 500m/s is enforced, yielding a maximum frequency of 1.0Hz.
This paper presents preliminary results of an experimental campaign on three beam-column joint subassemblies extracted from a 22-storey reinforced concrete frame building constructed in late 1980s at the Christchurch’s Central Business District (CBD) area, damaged and demolished after the 2010-2011 Canterbury earthquakes sequence (CES). The building was designed following capacity design principles. Column sway (i.e., soft storey) mechanisms were avoided, and the beams were provided with plastic hinge relocation details at both beam-ends, aiming at developing plastic hinges away from the column faces. The specimens were tested under quasi-static cyclic displacement controlled lateral loading. One of the specimens, showing no visible residual cracks was cyclically tested in its as-is condition. The other two specimens which showed residual cracks varying between hairline and 1.0mm in width, were subjected to cyclic loading to simulate cracking patterns consistent with what can be considered moderate damage. The cracked specimens were then repaired with an epoxy injection technique and subsequently retested until reaching failure. The epoxy injection techniques demonstrated to be quite efficient in partly, although not fully, restoring the energy dissipation capacities of the damaged specimens at all beam rotation levels. The stiffness was partly restored within the elastic range and almost fully restored after the onset of nonlinear behaviour.
Geological research in the immediate aftermath of the 2016 Kaikōura Earthquake, New Zealand, was necessary due to the importance and perishability of field data. It also reflects a real desire on the part of researchers to contribute not only to immediate scientific understanding but also to the societal recovery effort by enhancing knowledge of the event for the benefit of affected communities, civil defence organizations and regional and national decision makers. This commitment to outreach and engagement is consistent with the recent IAPG statement of Geoethics. More immediately, it was informed by experience of the 2010-2011 Canterbury Earthquake sequence. After that earlier disaster, intense interactions between researchers and various response agencies as well as local communities informed the development and dissemination of a set of ethical guidelines for researchers immediately following the Mw7.8 14 November 2016 Kaikōura Earthquake. In this presentation, I argue that ethical engagement of this kind is the key to gathering high quality research data immediately after the event. Creating trusting and mutually respectful, mutually beneficial relationships is also vital to ongoing engagement to facilitate further “in depth” research in collaboration with communities.
The focus of the study presented herein is an assessment of the relative efficacy of recent Cone Penetration Test (CPT) and small strain shear wave velocity (Vs) based variants of the simplified procedure. Towards this end Receiver Operating Characteristic (ROC) analyses were performed on the CPT- and Vs-based procedures using the field case history databases from which the respective procedures were developed. The ROC analyses show that Factors of Safety (FS) against liquefaction computed using the most recent Vs-based simplified procedure is better able to separate the “liquefaction” from the “no liquefaction” case histories in the Vs liquefaction database than the CPT-based procedure is able to separate the “liquefaction” from the “no liquefaction” case histories in the CPT liquefaction database. However, this finding somewhat contradicts the assessed predictive capabilities of the CPT- and Vs-based procedures as quantified using select, high quality liquefaction case histories from the 20102011 Canterbury, New Zealand, Earthquake Sequence (CES), wherein the CPT-based procedure was found to yield more accurate predictions. The dichotomy of these findings may result from the fact that different liquefaction field case history databases were used in the respective ROC analyses for Vs and CPT, while the same case histories were used to evaluate both the CPT- and Vs-based procedures.
The Mw 7.8 Kaikōura earthquake ruptured ~200 km at the ground surface across the New Zealand plate boundary zone in the northern South Island. This study was conducted in an area of ~600 km2 in the epicentral region where the faults comprise two main non-coplanar sets that strike E-NE and NNE-NW with mainly steep dips (60о-80°). Analysis of the surface rupture using field and LiDAR data provides new information on the dimensions, geometries and kinematics of these faults which was not previously available from pre-earthquake active faults or bedrock structure. The more northerly striking fault set are sub-parallel to basement bedding and accommodated predominantly left-lateral reverse slip with net slips of ~1 and ~5 m for the Stone Jug and Leader faults, respectively. The E-NE striking Conway-Charwell and The Humps faults accrued right-lateral to oblique reverse with net slips of ~2 and ~3 m, respectively. The faults form a hard-linked system dominated by kinematics consistent with the ~260° trend of the relative plate motion vector and the transpressional structures recorded across the plate boundary in the NE South Island. Interaction and intersection of the main fault sets facilitated propagation of the earthquake and transfer of slip northwards across the plate boundary zone.
The south Leader Fault (SLF) is a newly documented active structure that ruptured the surface during the Mw 7.8 Kaikoura earthquake. The Leader Fault is a NNE trending oblique left lateral thrust that links the predominantly right lateral ‘The Humps’ and Conway-Charwell faults. The present research uses LiDAR at 0.5 m resolution and field mapping to determine the factors controlling the surface geometries and kinematics of the south Leader Fault ruptures at the ground surface. The SLF zone is up to 2km wide and comprises a series of echelon NE-striking thrusts linked by near-vertical N-S striking faults. The thrusts are upthrown to the west by up to 1 m and dip 35-45°. Thrust slip surfaces are parallel with Cretaceous-Cenozoic bedding and may reflect flexural slip folding. By contrast, the northerly striking faults dip steeply (65° west- 85° east), and accommodate up to 3m of oblique left lateral displacement at the ground surface and displace Cenozoic bedding. Some of the SLF has been mapped in bedrock, although none were known to be active prior to the earthquake or have a strong topographic expression. The complexity of fault rupture and the width of the fault zone appears to reflect the occurrence of faulting and folding at the ground surface during the earthquake.
Surface rupture and slip from the Mw 7.8 2016 Kaikōura Earthquake have been mapped in the region between the Leader and Charwell rivers using field mapping and LiDAR data. The eastern Humps, north Leader and Conway-Charwell faults ruptured the ground surface in the study area. The E-NE striking ‘The Humps’ Fault runs along the base of the Mt Stewart range front, appears to dip steeply NW and intersects the NNW-NNE Leader Fault which itself terminates northwards at the NE striking Conway-Charwell Fault. The eastern Humps Fault is up to the NW and accommodates oblique slip with reverse and right lateral displacement. Net slip on ‘The Humps’ Fault is ≤4 m and produced ≤4 m uplift of the Mt Stewart range during the earthquake. The Leader Fault strikes NNW-NNE with dips ranging from ~10° west to 80° east and accommodated ≤4 m net slip comprising left-lateral and up-to-the-west vertical displacement. Like the Humps west of the study area, surface-rupture of the Leader Fault occurred on multiple strands. The complexity of rupture on the Leader Fault is in part due to the occurrence of bedding-parallel slip within the Cretaceous-Cenozoic sequence. Although the Mt Stewart range front is bounded by ‘The Humps’ Fault, in the study area neither this fault nor the Leader Fault were known to have been active before the earthquake. Fieldwork and trenching investigations are ongoing to characterise the geometry, kinematics and paleoseismic history of the mapped active faults.