The 4th of September 2010 Mw 7.1 Darfield (Canterbury) earthquake had generated significant ground shaking within the Christchurch Central Business District (CBD). Despite the apparently significant shaking, the observed structural damage for pre-1970s reinforced concrete (RC) buildings was indeed limited and lower than what was expected for such typology of buildings. This paper explores analytically and qualitatively the different aspects of the "apparent‟ good seismic performance of the pre-1970s RC buildings in the Christchurch CBD, following the earthquake reconnaissance survey by the authors. Damage and building parameters survey result, based on a previously established inventory of building stock of these non-ductile RC buildings, is briefly reported. From an inventory of 75 buildings, one building was selected as a numerical case-study to correlate the observed damage with the non-linear analyses. The result shows that the pre-1970s RC frame buildings performed as expected given the intensity of the ground motion shaking during the Canterbury earthquake. Given the brittle nature of this type of structure, it was demonstrated that more significant structural damage and higher probability of collapse could occur when the buildings were subjected to alternative input signals with different frequency content and duration characteristics and still compatible to the seismicity hazard for Christchurch CBD.
Natural disasters are increasingly disruptive events that affect livelihoods, organisations, and economies worldwide. Research has identified the impacts and responses of organisations to different types of natural disasters, and have outlined factors, such as industry sector, that are important to organisational vulnerability and resilience. One of the most costly types of natural disasters in recent years has been earthquakes, and yet to date, the majority of studies have focussed on the effects of earthquakes in urban areas, while rural organisational impact studies have primarily focused on the effects of meteorological and climatic driven hazards. As a result, the likely impacts of an earthquake on rural organisations in a developed context is unconstrained in the literature. In countries like New Zealand, which have major earthquakes and agricultural sectors that are significant contributors to the economy, it is important to know what impacts an earthquake event would have on the rural industries, and how these impacts compare to that of a more commonly analysed, high-frequency event. In September of 2010, rural organisations in Canterbury experienced the 4 September 2010 Mw 7.1 `Darfield' earthquake and the associated aftershocks, which came to be known as the Canterbury earth- quake sequence. The earthquake sequence caused intense ground shaking, creating widespread critical service outages, structural and non-structural damage to built infrastructure, as well as ground surface damage from ooding, liquefaction and surface rupture. Concurrently on September 18 2010, rural organisations in Southland experienced an unseasonably late snowstorm and cold weather snap that brought prolonged sub-zero temperatures, high winds and freezing rain, damaging structures in the City of Invercargill and causing widespread livestock losses and production decreases across the region. This thesis documents the effects of the Canterbury earthquake sequence and Southland snowstorm on farming and rural non-farming organisations, utilizing comparable methodologies to analyse rural organisational impacts, responses and recovery strategies to natural disasters. From the results, a short- term impact assessment methodology is developed for multiple disasters. Additionally, a regional asset repair cost estimation model is proposed for farming organisations following a major earthquake event, and the use of social capital in rural organisational recovery strategies following natural disasters is analysed.
Ground motion observations from the most significant 10 events in the 2010-2011 Canterbury earthquake sequence at near-source sites are utilized to scrutinize New Zealand (NZ)-specific pseudo-spectral acceleration (SA) empirical ground motion prediction equations (GMPE) (Bradley 2010, Bradley 2013, McVerry et al. 2006). Region-specific modification factors based on relaxing the conventional ergodic assumption in GMPE development were developed for the Bradley (2010) model. Because of the observed biases with magnitude and source-to-site distance for the McVerry et al. (2006) model it is not possible to develop region-specific modification factors in a reliable manner. The theory of non-ergodic empirical ground motion prediction is then outlined, and applied to this 10 event dataset to determine systematic effects in the between- and within-event residuals which lead to modifications in the predicted median and standard deviation of the GMPE. By examining these systematic effects over sub-regions containing a total of 20 strong motion stations within the Canterbury area, modification factors for use in region-specific ground motion prediction are proposed. These modification factors, in particular, are suggested for use with the Bradley et al. (2010) model in Canterbury-specific probabilistic seismic hazard analysis (PSHA) to develop revised design response, particularly for long vibration periods.
This poster presents preliminary results of ongoing experimental campaigns at the Universities of Auckland and Canterbury, aiming at investigating the seismic residual capacity of damaged reinforced concrete plastic hinges, as well as the effectiveness of epoxy injection techniques for restoring their stiffness, energy dissipation, and deformation capacity characteristics. This work is part of wider research project which started in 2012 at the University of Canterbury entitled “Residual Capacity and Repairing Options for Reinforced Concrete Buildings”, funded by the Natural Hazards Research Platform (NHRP). This research project aims at gaining a better understanding and providing the main end-users and stakeholders (practitioner engineers, owners, local and government authorities, insurers, and regulatory agencies) with comprehensive evidence-based information and practical guidelines to assess the residual capacity of damaged reinforced concrete buildings, as well as to evaluate the feasibility of repairing and thus support their delicate decision-making process of repair vs. demolition or replacement.
In the last two decades, New Zealand (NZ) has experienced significant earthquakes, including the 2010 M 7.2 Darfield, 2011 M 6.2 Christchurch, and 2016 M 7.8 Kaikōura events. Amongst these large events, tens of thousands of smaller earthquakes have occurred. While previous event and ground-motion databases have analyzed these events, many events below M 4 have gone undetected. The goal of this study is to expand on previous databases, particularly for small magnitude (M<4) and low-amplitude ground motions. This new database enables a greater understanding of regional variations within NZ and contributes to the validity of internationally developed ground-motion models. The database includes event locations and magnitude estimates with uncertainty considerations, and tectonic type assessed in a hierarchical manner. Ground motions are extracted from the GeoNet FDSN server and assessed for quality using a neural network classification approach. A deep neural network approach is also utilized for picking P and S phases for determination of event hypocentres. Relative hypocentres are further improved by double-difference relocation and will contribute toward developing shallow (< 50 km) seismic tomography models. Analysis of the resulting database is compared with previous studies for discussion of implications toward national hazard prediction models.
Despite the relatively low seismicity, a large earthquake in the Waikato region is expected to have a high impact, when the fourth-largest regional population and economy and the high density critical infrastructure systems in this region are considered. Furthermore, Waikato has a deep soft sedimentary basin, which increases the regional seismic hazard due to trapping and amplification of seismic waves and generation of localized surface waves within the basin. This phenomenon is known as the “Basin Effect”, and has been attributed to the increased damage in several historic earthquakes, including the 2010-2011 Canterbury earthquakes. In order to quantitatively model the basin response and improve the understanding of regional seismic hazard, geophysical methods will be used to develop shear wave velocity profiles across the Waikato basin. Active surface wave methods involve the deployment of linear arrays of geophones to record the surface waves generated by a sledge hammer. Passive surface wave methods involve the deployment of two-dimensional seismometer arrays to record ambient vibrations. At each site, the planned testing includes one active test and two to four passive arrays. The obtained data are processed to develop dispersion curves, which describe surface wave propagation velocity as a function of frequency (or wavelength). Dispersion curves are then inverted using the Geopsy software package to develop a suite of shear wave velocity profiles. Currently, more than ten sites in Waikato are under consideration for this project. This poster presents the preliminary results from the two sites that have been tested. The shear wave velocity profiles from all sites will be used to produce a 3D velocity model for the Waikato basin, a part of QuakeCoRE flagship programme 1.
Lake Taupō in New Zealand is associated with frequent unrest and small to moderate eruptions. It presents a high consequence risk scenario with immense potential for destruction to the community and the surrounding environment. Unrest associated with eruptions may also trigger earthquakes. While it is challenging to educate people about the hazards and risks associated with multiple eruptive scenarios, effective education of students can lead to better mitigation strategies and risk reduction. Digital resources with user-directed outcomes have been successfully used to teach action oriented skills relevant for communication during volcanic crisis [4]. However, the use of choose your own adventure strategies to enhance low probability risk literacy for Secondary school outreach has not been fully explored. To investigate how digital narrative storytelling can mediate caldera risk literacy, a module “The Kid who cried Supervolcano” will be introduced in two secondary school classrooms in Christchurch and Rotorua. The module highlights four learning objectives: (a) Super-volcanoes are beautiful but can be dangerous (b) earthquake (unrest) activity is normal for super-volcanoes (c) Small eruptions are possible from super-volcanoes and can be dangerous in our lifetimes (d) Super-eruptions are unlikely in our lifetimes. Students will create their digital narrative using the platform Elementari (www.elementari.io). The findings from this study will provide clear understanding of students’ understanding of risk perceptions of volcanic eruption scenarios and associated hazards and inform the design of educational resources geared towards caldera risk literacy.
This study investigates evidence for linkages and fault interactions centred on the Cust Anticline in Northwest Canterbury between Starvation Hill to the southwest and the Ashley and Loburn faults to the northeast. An integrated programme of geologic, geomorphic, paleo-seismic and geophysical analyses was undertaken owing to a lack of surface exposures and difficulty in distinguishing active tectonic features from fluvial and/or aeolian features across the low-relief Canterbury Plains. LiDAR analysis identified surface expression of several previously unrecognised active fault traces across the low-relief aggradation surfaces of the Canterbury Plains. Their presence is consistent with predictions of a fault relay exploiting the structural mesh across the region. This is characterised by interactions of northeast-striking contractional faults and a series of re-activating inherited Late Cretaceous normal faults, the latter now functioning as E–W-striking dextral transpressive faults. LiDAR also allowed for detailed analysis of the surface expression of individual faults and folds across the Cust Anticline contractional restraining bend, which is evolving as a pop-up structure within the newly established dextral shear system that is exploiting the inherited, now re-activated, basement fault zone. Paleo-seismic trenches were located on the crest of the western arm of the Cust Anticline and across a previously unrecognised E–W-striking fault trace, immediately southwest of the steeply plunging Cust Anticline termination. These studies confirmed the location and structural style of north-northeast-striking faults and an E–W-striking fault associated with the development of this structural culmination. A review of available industry seismic reflection lines emphasised the presence of a series of common structural styles having the same underlying structural drivers but with varying degrees of development and expression, both in the seismic profiles and in surface elevations across the study area. Based on LiDAR surface mapping and preliminary re-analysis of industry seismic reflection data, four fault zones are identified across the restraining bend structural culminations, which together form the proposed Oxford–Cust–Ashley Fault System. The 2010–2012 Canterbury Earthquake Sequence showed many similarities to the structural pattern established across the Oxford–Cust–Ashley Fault System, emphasising the importance of identification and characterization of presently hidden fault sources, and the understanding of fault network linkages, in order to improve constraints on earthquake source potential. Improved understanding of potentially-interactive fault sources in Northwest Canterbury, with the potential for combined initial fault rupture and spatial and temporal rupture propagation across this fault system, can be used in probabilistic seismic hazard analysis for the region, which is essential for the suitability and sustainability of future social and economic development.
Liquefaction during the 4th September 2010 Mw 7.1 Darfield earthquake and large aftershocks in 2011 (Canterbury earthquake sequence, CES) caused severe damage to land and infrastructure within Christchurch, New Zealand. Approximately one third of the total CES-induced financial losses were directly attributable to liq- uefaction and thus highlights the need for local and regional authorities to assess liquefaction hazards for present and future developments. This thesis is the first to conduct paleo-liquefaction studies in eastern Christchurch for the purpose of de- termining approximate return times of liquefaction-inducing earthquakes within the region. The research uncovered evidence for pre-CES liquefaction dated by radiocarbon and cross-cutting relationships as post-1660 to pre-1905. Additional paleo-liquefaction investigations within the eastern Christchurch suburb of Avon- dale, and the northern township of Kaiapoi, revealed further evidence for pre-CES liquefaction. Pre-CES liquefaction in Avondale is dated as post-1321 and pre-1901, while the Kaiapoi features likely formed during three distinct episodes: post-1458 and possibly during the 1901 Cheviot earthquake, post-1297 to pre-1901, and pre-1458. Evaluation of the liquefaction potential of active faults within the Can- terbury region indicates that many faults have the potential to cause widespread liquefaction within Avondale and Kaiapoi. The identification of pre-CES liquefac- tion confirms that these areas have previously liquefied, and indicates that residen- tial development in eastern Christchurch between 1860 and 2005 occurred in areas containing geologic evidence for pre-CES liquefaction. Additionally, on the basis of detailed field and GIS-based mapping and geospatial-statistical analysis, the distribution and severity of liquefaction and lateral spreading within the eastern Christchurch suburb of Avonside is shown in this study to be strongly in uenced by geomorphic and topographic variability. This variability is not currently ac- counted for in site-specific liquefaction assessments nor the simplified horizontal displacement models, and accounts for some of the variability between the pre- dicted horizontal displacements and those observed during the CES. This thesis highlights the potential applications of paleo-liquefaction investigations and ge- omorphic mapping to seismic and liquefaction hazard assessments and may aid future land-use planning decisions.
After a high-intensity seismic event, inspections of structural damages need to be carried out as soon as possible in order to optimize the emergency management, as well as improving the recovery time. In the current practice, damage inspections are performed by an experienced engineer, who physically inspect the structures. This way of doing not only requires a significant amount of time and high skilled human resources, but also raises the concern about the inspector’s safety. A promising alternative is represented using new technologies, such as drones and artificial intelligence, which can perform part of the damage classification task. In fact, drones can safely access high hazard components of the structures: for instance, bridge piers or abutments, and perform the reconnaissance by using highresolution cameras. Furthermore, images can be automatically processed by machine learning algorithms, and damages detected. In this paper, the possibility of applying such technologies for inspecting New Zealand bridges is explored. Firstly, a machine-learning model for damage detection by performing image analysis is presented. Specifically, the algorithm was trained to recognize cracks in concrete members. A sensitivity analysis was carried out to evaluate the algorithm accuracy by using database images. Depending on the confidence level desired,i.e. by allowing a manual classification where the alghortim confidence is below a specific tolerance, the accuracy was found reaching up to 84.7%. In the second part, the model is applied to detect the damage observed on the Anzac Bridge (GPS coordinates -43.500865, 172.701138) in Christchurch by performing a drone reconnaissance. Reults show that the accuracy of the damage detection was equal to 88% and 63% for cracking and spalling, respectively.
The Canterbury region of New Zealand experienced four earthquakes greater than MW 6.0 between September 2010 and December 2011. This study employs system dynamics as well as hazard, recovery and organisational literature and brings together data collected via surveys, case studies and interviews with organisations affected by the earthquakes. This is to show how systemic interactions and interdependencies within and between industry and geographic sectors affect their recovery post-disaster. The industry sectors in the study are: construction for its role in the rebuild, information and communication technology which is a regional high-growth industry, trucking for logistics, critical infrastructure, fast moving consumer goods (e.g. supermarkets) and hospitality to track recovery through non-discretionary and discretionary spend respectively. Also in the study are three urban centres including the region’s largest Central Business District, which has been inaccessible since the earthquake of 22 February 2011 to the time of writing in February 2013. This work also highlights how earthquake effects propagated between sectors and how sectors collaborated to mitigate difficulties such as product demand instability. Other interacting factors are identified that influence the recovery trajectories of the different industry sectors. These are resource availability, insurance payments, aid from central government, and timely and quality recovery information. This work demonstrates that in recovering from disaster it is crucial for organisations to identify what interacting factors could affect their operations. Also of importance are efforts to reduce the organisation’s vulnerability and increase their resilience to future crises and in day-to-day operations. Lastly, the multi-disciplinary approach to understanding the recovery and resilience of organisations and industry sectors after disaster, leads to a better understanding of effects as well as more effective recovery policy.
This paper investigates the effects of variability in source rupture parameters on site-specific physics-based simulated ground motions, ascertained through the systematic analysis of ground motion intensity measures. As a preliminary study, we consider simulations of the 22 February 2011 Christchurch earthquake using the Graves and Pitarka (2015) methodology. The effects of source variability are considered via a sensitivity study in which parameters (hypocentre location, earthquake magnitude, average rupture velocity, fault geometry and the Brune stress parameter) are individually varied by one standard deviation. The sensitivity of simulated ground motion intensity measures are subsequently compared against observational data. The preliminary results from this study indicate that uncertainty in the stress parameter and the rupture velocity have the most significant effect on the high frequency amplitudes. Conversely, magnitude uncertainty was found to be most influential on the spectral acceleration amplitudes at low frequencies. Further work is required to extend this preliminary study to exhaustively consider more events and to include parameter covariance. The ultimate results of this research will assist in the validation of the overall simulation method’s accuracy in capturing various rupture parameters, which is essential for the use of simulated ground motion models in probabilistic seismic hazard analysis.
Using case studies from the 2010-2011 Canterbury, New Zealand earthquake sequence, this study assesses the accuracies of paleoliquefaction back-analysis methods and explores the challenges, techniques, and uncertainties associated with their application. While liquefaction-based back-analyses have been widely used to estimate the magnitudes of paleoearthquakes, their uncertain efficacies continue to significantly affect the computed seismic hazard in regions where they are relied upon. Accordingly, their performance is evaluated herein using liquefaction data from modern earthquakes with known magnitudes. It is shown that when the earthquake source location and mechanism are known, back-analysis methods are capable of accurately deriving seismic parameters from liquefaction evidence. However, because the source location and mechanism are often unknown in paleoseismic studies, and because accurate interpretation is shown to be more difficult in such cases, new analysis techniques are proposed herein. An objective parameter is proposed to geospatially assess the likelihood of any provisional source location, enabling an analyst to more accurately estimate the magnitude of a liquefaction-inducing paleoearthquake. This study demonstrates the application of back-analysis methods, provides insight into their potential accuracies, and provides a framework for performing paleoliquefaction analyses worldwide.
The earthquake engineering community is currently grappling with the need to improve the post-earthquake reparability of buildings. As part of this, proposals exist to change design criteria for the serviceability limit state (SLS). This paper reviews options for change and considers how these could impact the expected repair costs for typical New Zealand buildings. The expected annual loss (EAL) is selected as a relevant measure or repair costs and performance because (i) EAL provides information on the performance of a building considering a range of intensity levels, (ii) the insurance industry refers to EAL when setting premiums, and (iii) monetary losses are likely to be correlated with loss of building functionality. The paper argues that because the expected annual loss is affected by building performance over a range of intensity levels, the definition of SLS criteria alone may be insufficient to effectively limit losses. However, it is also explained that losses could be limited effectively if the loadings standard were to set the SLS design intensity considering the potential implications on EAL. It is shown that in order to achieve similar values of EAL in Wellington and Christchurch, the return period intensity for SLS design would need to be higher in Christchurch owing to differences in local hazard conditions. The observations made herein are based on a simplified procedure for EAL estimation and hence future research should aim to verify the findings using a detailed loss assessment approach applied to a broad range of case study buildings.
In this paper we outline the process and outcomes of a multi-agency, multi-sector research collaboration, led by the Canterbury Earthquake Research Authority (CERA). The CERA Wellbeing Survey (CWS) is a serial, cross-sectional survey that is to be repeated six-monthly (in April and September) until the end of the CERA Act, in April 2016. The survey gathers self-reported wellbeing data to supplement the monitoring of the social recovery undertaken through CERA's Canterbury Wellbeing Index. Thereby informing a range of relevant agency decision-making, the CWS was also intended to provide the community and other sectors with a broad indication of how the population is tracking in the recovery. The primary objective was to ensure that decision-making was appropriately informed, with the concurrent aim of compiling a robust dataset that is of value to future researchers, and to the wider, global hazard and disaster research endeavor. The paper begins with an outline of both the Canterbury earthquake sequence, and the research context informing this collaborative project, before reporting on the methodology and significant results to date. It concludes with a discussion of both the survey results, and the collaborative process through which it was developed.
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.
This article argues that active coordination of research engagement after disasters has the potential to maximize research opportunities, improve research quality, increase end-user engagement, and manage escalating research activity to mitigate ethical risks posed to impacted populations. The focus is on the coordination of research activity after the 22nd February 2011 Mw6.2 Christchurch earthquake by the then newly-formed national research consortium, the Natural Hazards Research Platform, which included a social science research moratorium during the declared state of national emergency. Decisions defining this organisation’s functional and structural parameters are analyzed to identify lessons concerning the need for systematic approaches to the management of post disaster research, in collaboration with the response effort. Other lessons include the importance of involving an existing, broadly-based research consortium, ensuring that this consortium's coordination role is fully integrated into emergency management structures, and ensuring that all aspects of decision-making processes are transparent and easily accessed.
Geologic phenomena produced by earthquake shaking, including rockfalls and liquefaction features, provide important information on the intensity and spatiotemporal distribution of earthquake ground motions. The study of rockfall and liquefaction features produced in contemporary well- instrumented earthquakes increases our knowledge of how natural and anthropogenic environments respond to earthquakes and improves our ability to deduce seismologic information from analogous pre-contemporary (paleo-) geologic features. The study of contemporary and paleo- rockfall and liquefaction features enables improved forecasting of environmental responses to future earthquakes. In this thesis I utilize a combination of field and imagery-based mapping, trenching, stratigraphy, and numerical dating techniques to understand the nature and timing of rockfalls (and hillslope sedimentation) and liquefaction in the eastern South Island of New Zealand, and to examine the influence that anthropogenic activity has had on the geologic expressions of earthquake phenomena. At Rapaki (Banks Peninsula, NZ), field and imagery-based mapping, statistical analysis and numerical modeling was conducted on rockfall boulders triggered by the fatal 2011 Christchurch earthquakes (n=285) and compared with newly identified prehistoric (Holocene and Pleistocene) boulders (n=1049) deposited on the same hillslope. A significant population of modern boulders (n=26) travelled farther downslope (>150 m) than their most-travelled prehistoric counterparts, causing extensive damage to residential dwellings at the foot of the hillslope. Replication of prehistoric boulder distributions using 3-dimensional rigid body numerical models requires the application of a drag-coefficient, attributed to moderate to dense slope vegetation, to account for their spatial distribution. Radiocarbon dating provides evidence for 17th to early 20th century deforestation at the study site during Polynesian and European colonization and after emplacement of prehistoric rockfalls. Anthropocene deforestation enabled modern rockfalls to exceed the limits of their prehistoric predecessors, highlighting a shift in the geologic expression of rockfalls due to anthropogenic activity. Optical and radiocarbon dating of loessic hillslope sediments in New Zealand’s South Island is used to constrain the timing of prehistoric rockfalls and associated seismic events, and quantify spatial and temporal patterns of hillslope sedimentation including responses to seismic and anthropogenic forcing. Luminescence ages from loessic sediments constrain timing of boulder emplacement to between ~3.0 and ~12.5 ka, well before the arrival of Polynesians (ca AD 1280) and Europeans (ca AD 1800) in New Zealand, and suggest loess accumulation was continuing at the study site until 12-13 ka. Large (>5 m3) prehistoric rockfall boulders preserve an important record of Holocene hillslope sedimentation by creating local traps for sediment aggradation and upbuilding soil formation. Sediment accumulation rates increased considerably (>~10 factor increase) following human arrival and associated anthropogenic burning of hillslope vegetation. New numerical ages are presented to place the evolution of loess-mantled hillslopes in New Zealand’s South Island into a longer temporal framework and highlight the roles of earthquakes and humans on hillslope surface process. Extensive field mapping and characterization for 1733 individual prehistoric rockfall boulders was conducted at Rapaki and another Banks Peninsula site, Purau, to understand their origin, frequency, and spatial and volumetric distributions. Boulder characteristics and distributions were compared to 421 boulders deposited at the same sites during the 2010-2011 Canterbury earthquake sequence. Prehistoric boulders at Rapaki and Purau are comprised of two dominant lithofacies types: volcanic breccia and massive (coherent) lava basalt. Volcanic breccia boulders are found in greatest abundance (64-73% of total mapped rockfall) and volume (~90-96% of total rockfall) at both locations and exclusively comprise the largest boulders with the longest runout distances that pose the greatest hazard to life and property. This study highlights the primary influence that volcanic lithofacies architecture has on rockfall hazard. The influence of anthropogenic modifications on the surface and subsurface geologic expression of contemporary liquefaction created during the 2010-2011 Canterbury earthquake sequence (CES) in eastern Christchurch is examined. Trench observations indicate that anthropogenic fill layer boundaries and the composition/texture of discretely placed fill layers play an important role in absorbing fluidized sand/silt and controlling the subsurface architecture of preserved liquefaction features. Surface liquefaction morphologies (i.e. sand blows and linear sand blow arrays) display alignment with existing utility lines and utility excavations (and perforated pipes) provided conduits for liquefaction ejecta during the CES. No evidence of pre-CES liquefaction was identified within the anthropogenic fill layers or underlying native sediment. Radiocarbon dating of charcoal within the youngest native sediment suggests liquefaction has not occurred at the study site for at least the past 750-800 years. The importance of systematically examining the impact of buried infrastructure on channelizing and influencing surface and subsurface liquefaction morphologies is demonstrated. This thesis highlights the importance of using a multi-technique approach for understanding prehistoric and contemporary earthquake phenomena and emphasizes the critical role that humans play in shaping the geologic record and Earth’s surface processes.
Improving community resilience requires a way of thinking about the nature of a community. Two complementary aspects are proposed: the flows connecting the community with its surrounding environment and the resources the community needs for its ongoing life. The body of necessary resources is complex, with many interactions between its elements. A systems approach is required to understand the issues adequately. Community resilience is discussed in general terms together with strategies for improving it. The ideas are then illustrated and amplified by an extended case study addressing means of improving the resilience of a community on the West Coast of New Zealand to natural disasters. The case study is in two phases. The first relies on a mix of on-the-ground observations and constructed scenarios to provide recommendations for enhancing community resilience, while the second complements the first by developing a set of general lessons and issues to be addressed from observations of the Christchurch earthquakes of 2010 and 2011.
The Sendai Framework for Disaster Risk Reduction 2015-2030 finds that, despite progress in disaster risk reduction over the last decade “evidence indicates that exposure of persons and assets in all countries has increased faster than vulnerability has decreased, thus generating new risk and a steady rise in disaster losses” (p.4, UNISDR 2015). Fostering cooperation among relevant stakeholders and policy makers to “facilitate a science-policy interface for effective decisionmaking in disaster risk management” is required to achieve two priority areas for action, understanding disaster risk and enhancing disaster preparedness (p. 13, p. 23, UNISDR 2015). In other topic areas, the term science-policy interface is used interchangeably with the term boundary organisation. Both terms are usually used refer to systematic collaborative arrangements used to manage the intersection, or boundary, between science and policy domains, with the aim of facilitating the joint construction of knowledge to inform decision-making. Informed by complexity theory, and a constructivist focus on the functions and processes that minimize inevitable tensions between domains, this conceptual framework has become well established in fields where large complex issues have significant economic and political consequences, including environmental management, biodiversity, sustainable development, climate change and public health. To date, however, there has been little application of this framework in the disaster risk reduction field. In this doctoral project the boundary management framework informs an analysis of the research response to the 2010-2011 Canterbury Earthquake Sequence, focusing on the coordination role of New Zealand’s national Natural Hazards Research Platform. The project has two aims. It uses this framework to tell the nuanced story of the way this research coordination role evolved in response to both the complexity of the unfolding post-disaster environment, and to national policy and research developments. Lessons are drawn from this analysis for those planning and implementing arrangements across the science-policy boundary to manage research support for disaster risk reduction decision-making, particularly after disasters. The second aim is to use this case study to test the utility of the boundary management framework in the disaster risk reduction context. This requires that terminology and concepts are explained and translated in terms that make this analysis as accessible as possible across the disciplines, domains and sectors involved in disaster risk reduction. Key findings are that the focus on balance, both within organisations, and between organisations and domains, and the emphasis on systemic effects, patterns and trends, offer an effective and productive alternative to the more traditional focus on individual or organisational performance. Lessons are drawn concerning the application of this framework when planning and implementing boundary organisations in the hazard and disaster risk management context.
Tsunami events including the 2004 Indian Ocean Tsunami and the 2011 Tohoku Earthquake and Tsunami confirmed the need for Pacific-wide comprehensive risk mitigation and effective tsunami evacuation planning. New Zealand is highly exposed to tsunamis and continues to invest in tsunami risk awareness, readiness and response across the emergency management and science sectors. Evacuation is a vital risk reduction strategy for preventing tsunami casualties. Understanding how people respond to warnings and natural cues is an important element to improving evacuation modelling techniques. The relative rarity of tsunami events locally in Canterbury and also globally, means there is limited knowledge on tsunami evacuation behaviour, and tsunami evacuation planning has been largely informed by hurricane evacuations. This research aims to address this gap by analysing evacuation behaviour and movements of Kaikōura and Southshore/New Brighton (coastal suburb of Christchurch) residents following the 2016 Kaikōura earthquake. Stage 1 of the research is engaging with both these communities and relevant hazard management agencies, using a survey and community workshops to understand real-event evacuation behaviour during the 2016 Kaikōura earthquake and subsequent tsunami evacuations. The second stage is using the findings from stage 1 to inform an agent-based tsunami evacuation model, which is an approach that simulates of the movement of people during an evacuation response. This method improves on other evacuation modelling approaches to estimate evacuation times due to better representation of local population characteristics. The information provided by the communities will inform rules and interactions such as traffic congestion, evacuation delay times and routes taken to develop realistic tsunami evacuation models. This will allow emergency managers to more effectively prepare communities for future tsunami events, and will highlight recommended actions to increase the safety and efficiency of future tsunami evacuations.
Tsunami have the potential to cause significant disruptions to society, including damage to infrastructure, critical to the every-day operation of society. Effective risk management is required to reduce the potential tsunami impacts to them. Christchurch city, situated on the eastern coast of New Zealand’s South Island, is exposed to a number of far-field tsunami hazards. Although the tsunami hazard has been well identified for Christchurch city infrastructure, the likely impacts have not been well constrained. To support effective risk management a credible and realistic infrastructure impact model is required to inform risk management planning. The objectives of this thesis are to assess the impacts on Christchurch city infrastructure from a credible, hypothetical far-field tsunami scenario. To achieve this an impact assessment process is adopted, using tsunami hazard and exposure measures to determine asset vulnerability and subsequent impacts. However, the thesis identified a number of knowledge gaps in infrastructure vulnerability to tsunami. The thesis addresses this by using two approaches: a tsunami damage matrix; and the development of tsunami fragility functions. The tsunami damage matrix pools together tsunami impacts on infrastructure literature, and post-event field observations. It represents the most comprehensive ‘look-up’ resource for tsunami impacts to infrastructure to date. This damage matrix can inform the assessment of tsunami impacts on Christchurch city infrastructure by providing a measure of damage likelihood at various hazard intensities. A more robust approach to tsunami vulnerability of infrastructure are fragility functions, which are also developed in this thesis. These were based on post-event tsunami surveys of the 2011 ‘Tohoku’ earthquake tsunami in Japan. The fragility functions are limited to road and bridge infrastructure, but represent the highest resolution measure of vulnerability for the given assets. As well as providing a measure of damage likelihood for a given tsunami hazard intensity, these also indicate a level of asset damage. The impact assessment process, and synthesized vulnerability measures, are used to run tsunami impact models for Christchurch infrastructure to determine the probability of asset damage occurring and to determine if impact will reach or exceed a given damage state. The models suggest that infrastructure damage is likely to occur in areas exposed to tsunami inundation in this scenario, with significant damage identified for low elevation roads and bridges. The results are presented and discussed in the context of the risk management framework, with emphasis on using risk assessment to inform risk treatment, monitoring and review. In summary, this thesis A) advances tsunami vulnerability and impact assessment methodologies for infrastructure and B) provides a tsunami impact assessment framework for Christchurch city infrastructure which will inform infrastructure tsunami risk management for planners, emergency managers and lifelines groups.
This study explicitly investigates uncertainties in physics-based ground motion simulation validation for earthquakes in the Canterbury region. The simulations utilise the Graves and Pitarka (2015) hybrid methodology, with separately quantified parametric uncertainties in the comprehensive physics and simplified physics components of the model. The study is limited to the simulation of 148 small magnitude (Mw 3.5 – 5) earthquakes, with a point source approximation for the source rupture representations, which also enables a focus on a small number of relevant uncertainties. The parametric uncertainties under consideration were selected through sensitivity analysis, and specifically include: magnitude, Brune stress parameter and high frequency rupture velocity. Twenty Monte Carlo realisations were used to sample parameter uncertainties for each of the 148 events. Residuals associated with the following intensity measures: spectral acceleration, peak ground velocity, arias intensity and significant duration, were ascertained. Using these residuals, validation was performed through assessment of systematic biases in site and source terms from mixed-effects regression. Based on the results to date, initial standard deviation recommendations for parameter uncertainties, based on the Canterbury simulations have been obtained. This work ultimately provides an initial step toward explicit incorporation of modelling uncertainty in simulated ground motion predictions for future events, which will improve the use of simulation models in seismic hazard analysis. We plan to subsequently assess uncertainties for larger magnitude events with more complex ruptures, and events across a larger geographic region, as well as uncertainties due to path attenuation, site effects, and more general model epistemic uncertainties.
On 14 November 2016, the Mw 7.8 Kaikōura earthquake caused widespread damage along the east coast of the South Island, New Zealand. Kaikōura town itself was isolated from the rest of the country by landslides blocking off major roads. While impacts from the Kaikōura earthquake on large, urban population centres have been generally well documented, this thesis aims to fill gaps in academic knowledge regarding small rural towns. This thesis investigates what, where and when critical infrastructure and lifeline service disruption occurred following the 2016 Kaikōura earthquake in a selection of small towns, and how the communities in these areas adapted to disruption. Following a robust review of literature and news media, four small rural towns were selected from North Canterbury (Culverden & Waiau) and Marlborough (Seddon & Ward) in the South Island, New Zealand. Semi-structured interview sessions with a special focus on these towns were held with infrastructure managers, emergency response and recovery officials, and organisation leaders with experience or expertise in the 2016 Kaikōura earthquake. Findings were supplemented with emergency management situation reports to produce hazard maps and infrastructure exposure maps. A more detailed analysis was conducted for Waiau involving interdependence analyses and a level of service timeline for select lifeline services. The earthquake impacted roads by blocking them with landslides, debris and surface rupture. Bridges where shaken off their abutments, breaking infrastructure links such as fibre landlines as they went. Water supplies and other forms of infrastructure relied heavily on the level of service of roads, as rough rural terrain left few alternatives. Adapting to an artificial loss of road service, some Waiau locals created their own detour around a road cordon in order to get home to family and farms. Performance of dwellings was tied to socioeconomic factors as much as proximity to the epicentre. Farmers who lost water access pulled out fences to allow stock to drink from rivers. Socioeconomic differences between farmland and township residents also contributed to resilience variations between the towns assessed in this study. Understanding how small rural towns respond and adapt to disaster allows emergency management officials and policy to be well informed and flexible with planning for multiple size classes of towns.
Recently developed performance-based earthquake engineering framework, such as one provided by PEER (Deierlein et al. 2003), assist in the quantification in terms of performance such as casualty, monetary losses and downtime. This opens up the opportunity to identify cost-effective retrofit/rehabilitation strategies by comparing upfront costs associated with retrofit with the repair costs that can be expected over time. This loss assessment can be strengthened by learning from recent earthquakes, such as the 2010 Canterbury and 2016 Kaikoura earthquakes. In order to investigate which types of retrofit/rehabilitation strategies may be most cost-effective, a case study building was chosen for this research. The Pacific Tower, a 22-storey EBF apartment located within the Christchurch central business district (CBD), was damaged and repaired during the 2010 Canterbury earthquake series. As such, by taking hazard levels accordingly (i.e. to correspond to the Christchurch CBD), modelling and analysing the structure, and considering the vulnerability and repair costs of its different components, it is possible to predict the expected losses of the aforementioned building. Using this information, cost-effective retrofit/rehabilitation strategy can be determined. This research found that more often than not, it would be beneficial to improve the performance of valuable non-structural components, such as partitions. Although it is true that improving such elements will increase the initial costs, over time, the benefits gained from reduced losses should be expected to overcome the initial costs. Aftershocks do increase the predicted losses of a building even in lower intensities due to the fact that non-structural components can get damaged at such low intensities. By comparing losses computed with and without consideration of aftershocks for a range of historical earthquakes, it was found that the ratio between losses due to main shock with aftershocks to the losses due to the main shock only tended to increase with increasing main shock magnitude. This may be due to the fact that larger magnitude earthquakes tend to generate larger magnitude aftershocks and as those aftershocks happen within a region around the main shock, they are more likely to cause intense shaking and additional damage. In addition to this observation, it was observed that the most significant component of loss of the case study building was the non-structural partition walls.
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.
Recent tsunami events have highlighted the importance of effective tsunami risk management strategies (including land-use planning, structural and natural mitigation, warning systems, education and evacuation planning). However, the rarity of tsunami means that empirical data concerning reactions to tsunami warnings and evacuation behaviour is rare when compared to findings for evacuations from other hazards. More knowledge is required to document the full evacuation process, including responses to warnings, pre-evacuation actions, evacuation dynamics, and the return home. Tsunami evacuation modelling has the potential to inform evidence-based tsunami risk planning and response. However, to date, tsunami evacuation models have largely focused on the timings of evacuations, rather than behaviours of those evacuating. In this research, survey data was gathered from coastal communities in Banks Peninsula and Christchurch, New Zealand, relating to behaviours and actions during the November 14th 2016 Kaikōura earthquake tsunami. Survey questions asked about immediate actions following the earthquake shaking, reactions to tsunami warnings, pre-evacuation actions, evacuation dynamics and details on congestion. This data was analysed to characterise trends and identify factors that influenced evacuation actions and behaviour, and was further used to develop a realistic evacuation model prototype to evaluate the capacity of the roading network in Banks Peninsula during a tsunami evacuation. The evacuation model incorporated tsunami risk management strategies that have been implemented by local authorities, and exposure and vulnerability data, alongside the empirical data collected from the survey. This research enhances knowledge of tsunami evacuation behaviour and reactions to tsunami warnings, and can be used to refine evacuation planning to ensure that people can evacuate efficiently, thereby reducing their tsunami exposure and personal risk.
Disaster recovery is significantly affected by funding availability. The timeliness and quality of recovery activities are not only impacted by the extent of the funding but also the mechanisms with which funding is prioritised, allocated and delivered. This research addresses the impact of funding mechanisms on the effectiveness and efficiency of post-disaster demolition and debris management programmes. A qualitative assessment of the impacts on recovery of different funding sources and mechanisms was carried out, using the 2010 Canterbury Earthquake as well as other recent international events as case studies. The impacts assessed include: timeliness, completeness, environmental, economic and social impacts. Of the case studies investigated, the Canterbury Earthquake was the only disaster response to rely solely on a privatised approach to insurance for debris management. Due to the low level of resident displacement and low level of hazard in the waste, this was a satisfactory approach, though not ideal. This approach has led to greater organisational complexity and delays. For many other events, the potential community wide impacts caused by the prolonged presence of disaster debris means that publicly funded and centrally facilitated programmes appear to be the most common and effective method of managing disaster waste.
This dissertation addresses several fundamental and applied aspects of ground motion selection for seismic response analyses. In particular, the following topics are addressed: the theory and application of ground motion selection for scenario earthquake ruptures; the consideration of causal parameter bounds in ground motion selection; ground motion selection in the near-fault region where directivity effect is significant; and methodologies for epistemic uncertainty consideration and propagation in the context of ground motion selection and seismic performance assessment. The paragraphs below outline each contribution in more detail. A scenario-based ground motion selection method is presented which considers the joint distribution of multiple intensity measure (IM) types based on the generalised conditional intensity measure (GCIM) methodology (Bradley, 2010b, 2012c). The ground motion selection algorithm is based on generating realisations of the considered IM distributions for a specific rupture scenario and then finding the prospective ground motions which best fit the realisations using an optimal amplitude scaling factor. In addition, using different rupture scenarios and site conditions, two important aspects of the GCIM methodology are scrutinised: (i) different weight vectors for the various IMs considered; and (ii) quantifying the importance of replicate selections for ensembles with different numbers of desired ground motions. As an application of the developed scenario-based ground motion selection method, ground motion ensembles are selected to represent several major earthquake scenarios in New Zealand that pose a significant seismic hazard, namely, Alpine, Hope and Porters Pass ruptures for Christchurch city; and Wellington, Ohariu, and Wairarapa ruptures for Wellington city. A rigorous basis is developed, and sensitivity analyses performed, for the consideration of bounds on causal parameters (e.g., magnitude, source-to-site distance, and site condition) for ground motion selection. The effect of causal parameter bound selection on both the number of available prospective ground motions from an initial empirical as-recorded database, and the statistical properties of IMs of selected ground motions are examined. It is also demonstrated that using causal parameter bounds is not a reliable approach to implicitly account for ground motion duration and cumulative effects when selection is based on only spectral acceleration (SA) ordinates. Specific causal parameter bounding criteria are recommended for general use as a ‘default’ bounding criterion with possible adjustments from the analyst based on problem-specific preferences. An approach is presented to consider the forward directivity effects in seismic hazard analysis, which does not separate the hazard calculations for pulse-like and non-pulse-like ground motions. Also, the ability of ground motion selection methods to appropriately select records containing forward directivity pulse motions in the near-fault region is examined. Particular attention is given to ground motion selection which is explicitly based on ground motion IMs, including SA, duration, and cumulative measures; rather than a focus on implicit parameters (i.e., distance, and pulse or non-pulse classifications) that are conventionally used to heuristically distinguish between the near-fault and far-field records. No ad hoc criteria, in terms of the number of directivity ground motions and their pulse periods, are enforced for selecting pulse-like records. Example applications are presented with different rupture characteristics, source-to-site geometry, and site conditions. It is advocated that the selection of ground motions in the near-fault region based on IM properties alone is preferred to that in which the proportion of pulse-like motions and their pulse periods are specified a priori as strict criteria for ground motion selection. Three methods are presented to propagate the effect of seismic hazard and ground motion selection epistemic uncertainties to seismic performance metrics. These methods differ in their level of rigor considered to propagate the epistemic uncertainty in the conditional distribution of IMs utilised in ground motion selection, selected ground motion ensembles, and the number of nonlinear response history analyses performed to obtain the distribution of engineering demand parameters. These methods are compared for an example site where it is observed that, for seismic demand levels below the collapse limit, epistemic uncertainty in ground motion selection is a smaller uncertainty contributor relative to the uncertainty in the seismic hazard itself. In contrast, uncertainty in ground motion selection process increases the uncertainty in the seismic demand hazard for near-collapse demand levels.
High-quality ground motion records are required for engineering applications including response history analysis, seismic hazard development, and validation of physics-based ground motion simulations. However, the determination of whether a ground motion record is high-quality is poorly handled by automation with mathematical functions and can become prohibitive if done manually. Machine learning applications are well-suited to this problem, and a previous feed-forward neural network was developed (Bellagamba et al. 2019) to determine high-quality records from small crustal events in the Canterbury and Wellington regions for simulation validation. This prior work was however limited by the omission of moderate-to-large magnitude events and those from other tectonic environments, as well as a lack of explicit determination of the minimum usable frequency of the ground motion. To address these shortcomings, an updated neural network was developed to predict the quality of ground motion records for all magnitudes and all tectonic sources—active shallow crustal, subduction intraslab, and subduction interface—in New Zealand. The predictive performance of the previous feed-forward neural network was matched by the neural network in the domain of small crustal records, and this level of predictive performance is now extended to all source magnitudes and types in New Zealand making the neural network applicable to global ground motion databases. Furthermore, the neural network provides quality and minimum usable frequency predictions for each of the three orthogonal components of a record which may then be mapped into a binary quality decision or otherwise applied as desired. This framework provides flexibility for the end user to predict high-quality records with various acceptability thresholds allowing for this neural network to be used in a range of applications.
