Beach ridge stratigraphy can provide an important record of both sustained coastal progradation and responses to events such as extreme storms, as well as evidence of earthquake induced sediment pulses. This study is a stratigraphic investigation of the late Holocene mixed sand gravel (MSG) beach ridge plain on the Canterbury coast, New Zealand. The subsurface was imaged along a 370 m shore-normal transect using 100 and 200 MHz ground penetrating radar (GPR) antennae, and cored to sample sediment textures. Results show that, seaward of a back-barrier lagoon, the Pegasus Bay beach ridge plain prograded almost uniformly, under conditions of relatively stable sea level. Nearshore sediment supply appears to have created a sustained sediment surplus, perhaps as a result of post-seismic sediment pulses, resulting in a flat, morphologically featureless beach ridge plain. Evidence of a high magnitude storm provides an exception, with an estimated event return period in excess of 100 years. Evidence from the GPR sequence combined with modern process observations from MSG beaches indicates that a paleo storm initially created a washover fan into the back-barrier lagoon, with a large amount of sediment simultaneously moved off the beach face into the nearshore. This erosion event resulted in a topographic depression still evident today. In the subsequent recovery period, sediment was reworked by swash onto the beach as a sequence of berm deposit laminations, creating an elevated beach ridge that also has a modern-day topographic signature. As sediment supply returned to normal, and under conditions of falling sea level, a beach ridge progradation sequence accumulated seaward of the storm feature out to the modern-day beach as a large flat, uniform progradation plain. This study highlights the importance of extreme storm events and earthquake pulses on MSG coastlines in triggering high volume beach ridge formation during the subsequent recovery period.
Science education research shows that a traditional, stand-and-deliver lecture format is less effective than teaching strategies that are learner-centred and that promote active engagement. The Carl Wieman Science Education Initiative (CWSEI) has used this research to develop resources to improve learning in university science courses. We report on a successful adaptation and implementation of CWSEI in the New Zealand university context. This two-year project at Massey University and the University of Canterbury began by using perception and concept surveys before and after undergraduate science courses to measure students’ attitudes towards science as well as their knowledge. Using these data, and classroom observations of student engagement and corroborating focus groups, the research team worked with lecturers to create interventions to enhance student engagement and learning in those courses. Results show several positive changes related to these interventions and they suggest several recommendations for lecturers and course coordinators. The recommendations include:1. Make learning outcomes clear, both for the lecturer and the students; this helps to cull extraneous material and scaffold student learning. 2. Use interactive activities to improve engagement, develop deeper levels of thinking, and improve learning. 3. Intentionally foster “expert-like thinking” amongst students in the first few semesters of the degree programme. 4. Be flexible because one size does not fit all and contextual events are beyond anyone’s control.In addition to these recommendations, data collected at the Canterbury site during the 2010 and 2011 earthquakes reinforced the understanding that the most carefully designed teaching innovations are subject to contextual conditions beyond the control of academics.
In response to the February 2011 earthquake, Parliament enacted the Canterbury Earthquake Recovery Act. This emergency legislation provided the executive with extreme powers that extended well beyond the initial emergency response and into the recovery phase. Although New Zealand has the Civil Defence Emergency Management Act 2002, it was unable to cope with the scale and intensity of the Canterbury earthquake sequence. Considering the well-known geological risk facing the Wellington region, this paper will consider whether a standalone “Disaster Recovery Act” should be established to separate an emergency and its response from the recovery phase. Currently, Government policy is to respond reactively to a disaster rather than proactively. In a major event, this typically involves the executive being given the ability to make rules, regulations and policy without the delay or oversight of normal legislative process. In the first part of this paper, I will canvas what a “Disaster Recovery Act” could prescribe and why there is a need to separate recovery from emergency. Secondly, I will consider the shortfalls in the current civil defence recovery framework which necessitates this kind of heavy governmental response after a disaster. In the final section, I will examine how
In this paper we introduce CityViewAR, a mobile outdoor Augmented Reality (AR) application for providing AR information visualization on a city scale. The CityViewAR application was developed to provide geographical information about the city of Christchurch, which was hit by several major earthquakes in 2010 and 2011. The application provides information about destroyed buildings and historical sites that were affected by the earthquakes. The geo-located content is provided in a number of formats including 2D map views, AR visualization of 3D models of buildings on-site, immersive panorama photographs, and list views. The paper describes the iterative design and implementation details of the application, and gives one of the first examples of a study comparing user response to AR and non-AR viewing in a mobile tourism application. Results show that making such information easily accessible to the public in a number of formats could help people to have richer experience about cities. We provide guidelines that will be useful for people developing mobile AR applications for city-scale tourism or outdoor guiding, and discuss how the underlying technology could be used for applications in other areas.
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.
The Eastern Humps and Leader faults, situated in the Mount Stewart Range in North Canterbury, are two of the ≥17 faults which ruptured during the 2016 MW7.8 Kaikōura Earthquake. The earthquake produced complex, intersecting ground ruptures of these faults and the co-seismic uplift of the Mount Stewart Range. This thesis aims to determine how these two faults accommodated deformation during the 2016 earthquake and how they interact with each other and with pre-existing geological structures. In addition, it aims to establish the most likely subsurface geometry of the fault complex across the Mount Stewart Range, and to investigate the paleoseismic history of the Leader Fault. The Eastern Humps Fault strikes ~240° and dips 80° to 60° to the northwest and accommodated right- lateral – reverse-slip, with up to 4 m horizontal and 2 m vertical displacement in the 2016 earthquake. The strike of the Leader Fault varies from ~155 to ~300°, and dips ~30 to ~80° to the west/northwest, and mainly accommodated left-lateral – reverse-slip of up to 3.5 m horizontal and 3.5 m vertical slip in the 2016 earthquake. On both the Eastern Humps and Leader faults the slip is variable along strike, with areas of low total displacement and areas where horizontal and vertical displacement are negatively correlated. Fault traces with low total displacement reflect the presence of off-fault (distributed) displacement which is not being captured with field measurements. The negative correlation of horizontal and vertical displacement likely indicates a degree of slip partitioning during the 2016 earthquake on both the Eastern Humps and Leader faults. The Eastern Humps and Leader faults have a complex, interdependent relationship with the local bedrock geology. The Humps Fault appears to be a primary driver of ongoing folding and deformation of the local Mendip Syncline and folding of the Mount Stewart Range, which probably began prior to, or synchronous with, initial rupture of The Humps Fault. The Leader Fault appears to use existing lithological weaknesses in the Cretaceous-Cenozoic bedrock stratigraphy to rupture to the surface. This largely accounts for the strong variability on the strike and dip of the Leader Fault, as the geometry of the surface ruptures tend to reflect the strike and dip of the geological strata which it is rupturing through. The Leader Fault may also accommodate some degree of flexural slip in the Cenozoic cover sequence of the Mendip Syncline, contributing to the ongoing growth of the fold. The similarity between topography and uplift profiles from the 2016 earthquake suggest that growth of the Mount Stewart Range has been primarily driven by multiple (>500) discrete earthquakes that rupture The Humps and Leader faults. The spatial distribution of surface displacements across the Mount Stewart Range is more symmetrical than would be expected if uplift is driven primarily by The Humps and Leader faults alone. Elastic dislocation forward models were used to model potential sub-surface geometries and the resulting patterns of deformation compared to photogrammetry-derived surface displacements. Results show a slight preference for models with a steeply southeast-dipping blind fault, coincident with a zone of seismicity at depth, as a ‘backthrust’ to The Humps and Leader faults. This inferred Mount Stewart Fault accommodated contractional strain during the 2016 earthquake and contributes to the ongoing uplift of the Mount Stewart Range with a component of folding. Right-lateral and reverse shear stress change on the Hope Fault was also modelled using Coulomb 3.3 software to examine whether slip on The Humps and Leader faults could transfer enough stress onto the Hope Fault to trigger through-going rupture. Results indicate that during the 2016 earthquake right-lateral shear and reverse stress only increased on the Hope Fault in small areas to the west of the Leader Fault, and similar ruptures would be unlikely to trigger eastward propagating rupture unless the Hope Fault was close to failure prior to the earthquake. Paleoseismic trenches were excavated on the Leader Fault at four locations from 2018 to 2020, revealing near surface (< 4m depth) contractional deformation of Holocene stratigraphy. Three of the trench locations uncovered clear evidence for rupture of the Leader Fault prior to 2016, with fault displacement of near surface stratigraphy being greater than displacement recorded during the 2016 earthquake. Radiocarbon dating of in-situ organic material from two trenches indicate a date of the penultimate earthquake on the Leader Fault within the past 1000 years. This date is consistent with The Humps and Leader faults having ruptured simultaneously in the past, and with multi-fault ruptures involving The Humps, Leader, Hundalee and Stone Jug faults having occurred prior to the 2016 Kaikōura earthquake. Overall, the results contribute to an improved understanding of the Kaikōura earthquake and highlight the importance of detailed structural and paleoseismic investigations in determining controls on earthquake ‘complexity’.
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.
In this paper we apply Full waveform tomography (FWT) based on the Adjoint-Wavefield (AW) method to iteratively invert a 3-D geophysical velocity model for the Canterbury region (Lee, 2017) from a simple initial model. The seismic wavefields was generated using numerical solution of the 3-D elastodynamic/ visco- elastodynamic equations (EMOD3D was adopted (Graves, 1996)), and through the AW method, gradients of model parameters (compression and shear wave velocity) were computed by implementing the cross-adjoint of forward and backward wavefields. The reversed-in-time displacement residual was utilized as the adjoint source. For inversion, we also account for the near source/ station effects, gradient precondition, smoothening (Gaussian filter in spatial domain) and optimal step length. Simulation-to-observation misfit measurements based on 191 sources at 78 seismic stations in the Canterbury region (Figure 1) were used into our inversion. The inversion process includes multiple frequency bands, starting from 0-0.05Hz, and advancing to higher frequency bands (0-0.1Hz and 0-0.2Hz). Each frequency band was used for up to 10 iterations or no optimal step length found. After 3 FWT inversion runs, the simulated seismograms computed using our final model show a good matching with the observed seismograms at frequencies from 0 - 0.2 Hz and the normalized least-squared misfit error has been significantly reduced. Over all, the synthetic study of FWT shows a good application to improve the crustal velocity models from the existed geological models and the seismic data of the different earthquake events happened in the Canterbury region.
Field trips are one of the most critical pieces of learning for students in sciences like geology, biology, and geography. Virtual field trips (VFT) are being increasingly considered as sophisticated and effective forms of teaching, especially with the rise of new technologies and the growing demand for more inclusive classroom environments. This research developed a virtual field trip for Tertiary students in an introductory-level geology course (GEOL 113: Environmental Geohazards) at the University of Canterbury. This initiative was in partnership with LEARNZ – a highly esteemed virtual fieldtrip team run by CORE Education that creates successful VFTs for Primary and Secondary students in New Zealand. Key components of the Tertiary VFT include a student acting as the virtual field trip teacher interviewing experts and leading the field trip, web-based background material, online assessment, and photos. In two successive academic years, students participated in the VFT during lectures and as pre class assignments prior to a one-day earthquake hazards workshop. In 2016, the virtual field trip used the LEARNZ web platform and occurred synchronously with the class; in 2017 the virtual fieldtrip reused the video, images and word documents from the previous year with the addition of a Google Earth component and with no reliance on the LEARNZ web platform. The goals of the trip were designed to prepare students for the earthquake hazards workshop, in which students analysed earthquake impacts over varying timescales and then applied that knowledge to develop strategies for the recovery of three crucial industries (dairy, mining, or tourism) on the West Coast of New Zealand’s South Island. In both years, number of clicks data showed that students interacted with online material far more during this week of the course than any other. Following the synchronous version in 2016, the students who were surveyed reported (1) they enjoyed the trip, (2) they found background material useful for preparation for the trip and the workshop, and (3) the additional work was at the appropriate level. Despite predominantly positive responses from the students, we experienced some negative feedback from participating staff mainly associated with stress and technical difficulties in running the synchronous VFT. With the asynchronous trip in 2017, staff reported a highly positive overall experience, with a perceived enhanced interaction with class during lecture time, and an increased and enhanced engagement with course material outside of class. The student survey again showed that the majority of students surveyed enjoyed the virtual fieldtrip, and that it was useful preparation for the workshop. Additionally, they reported an improved link between earth processes and society, which was a key overarching aim for the course. We propose that the synchronous version poses more excitement and immersion in the field environment, whereas the reuse of the asynchronous version increases the utility (and hence value for money) of the trip, and minimises technical difficulties and lecturer stress. Additionally, re-using the material in the asynchronous version offered opportunities to improve and supplement the past content, such as the incorporation of following an annotated trip path in Google Earth. As recommendations for others interested in developing virtual fieldtrips, we report that the design of a virtual fieldtrip should include (1) Goal-aligned content and assessment for both practice and marks, (2) a student and instructor experience that is authentic and flexible to both the people and the place. We suggest that these aims can be achieved whatever the budget or timeframe and make our material freely available at https://serc.carleton.edu/index.html.
Farming and urban regions are impacted by earthquake disasters in different ways, and feature a range of often different recovery requirements. In New Zealand, and elsewhere, most earthquake impact and recovery research is urban focused. This creates a research deficit that can lead to the application of well-researched urban recovery strategies in rural areas to suboptimal effect. To begin to reduce this deficit, in-depth case studies of the earthquake impacts and recovery of three New Zealand farms severely impacted by the 14th November 2016, M7.8 Hurunui-Kaikōura earthquake were conducted. The initial earthquake, its aftershocks and coseismic hazards (e.g., landslides, liquefaction, surface rupture) affected much of North Canterbury, Marlborough and the Wellington area. The three case study farms were chosen to broadly represent the main types of farming and topography in the Hurunui District in North Canterbury. The farms were directly and indirectly impacted by earthquakes and related hazards. On-farm infrastructure (e.g., woolsheds, homesteads) and essential services (e.g., water, power), frequently sourced from distributed networks, were severely impacted. The earthquake occurred after two years of regional drought had already stressed farm systems and farmers to restructuring or breaking point. Cascading interlinked hazards stemming from the earthquakes and coseismic hazards continued to disrupt earthquake recovery over a year after the initial earthquake. Semi-structured interviews with the farmers were conducted nine and fourteen months after the initial earthquake to capture the timeline of on-going impacts and recovery. Analysis of both geological hazard data and interview data resulted in the identification of key factors influencing farm level earthquake impact and recovery. These include pre-existing conditions (e.g., drought); farm-specific variations in recovery timelines; and resilience strategies for farm recovery resources. The earthquake recovery process presented all three farms with opportunities to change their business plans and adapt to mitigate on-going and future risk.
© 2017 The Royal Society of New Zealand. This paper discusses simulated ground motion intensity, and its underlying modelling assumptions, for great earthquakes on the Alpine Fault. The simulations utilise the latest understanding of wave propagation physics, kinematic earthquake rupture descriptions and the three-dimensional nature of the Earth's crust in the South Island of New Zealand. The effect of hypocentre location is explicitly examined, which is found to lead to significant differences in ground motion intensities (quantified in the form of peak ground velocity, PGV) over the northern half and southwest of the South Island. Comparison with previously adopted empirical ground motion models also illustrates that the simulations, which explicitly model rupture directivity and basin-generated surface waves, lead to notably larger PGV amplitudes than the empirical predictions in the northern half of the South Island and Canterbury. The simulations performed in this paper have been adopted, as one possible ground motion prediction, in the ‘Project AF8’ Civil Defence Emergency Management exercise scenario. The similarity of the modelled ground motion features with those observed in recent worldwide earthquakes as well as similar simulations in other regions, and the notably higher simulated amplitudes than those from empirical predictions, may warrant a re-examination of regional impact assessments for major Alpine Fault earthquakes.
This paper presents an overview of the soil profile characteristics at a number of strong motion station (SMS) sites in Christchurch and its surrounds. An extensive database of ground motion records has been captured by the SMS network in the Canterbury region. However in order to comprehensively understand the ground motions recorded at these sites and to be able to relate these motions to other locations, a detailed understanding of the geotechnical profile at each SMS is required. The original NZS1170.5 (SNZ 2004) site subsoil classifications for each SMS site based on regional geological information and well logs located at varying distances from the site. Given the variability of Christchurch soils, more detailed investigations are required in close vicinity to each SMS. In this regard, CPT, SPT and borehole data, and shear wave velocity (Vs) profiles in close vicinity to the SMS are currently being used to develop representative soil profiles at each site. Site subsoil classifications based on Vs measurements performed by the authors do not always agree with the original classifications, often indicating that a softer site class is appropriate. However, SPT N values often indicate a stiffer site class than the Vs data, in some cases also disagreeing with prior assumed classifications. Hence, the recent site investigation data presented herein highlights the importance of having detailed site-specific information at SMS locations in order to properly classify them. Furthermore, additional studies are required to harmonize site classification based on SPT N and Vs-.
Tens of thousands of landslides were generated over 10, 000 km2 of North Canterbury and Marlborough as a consequence of the 14 November 2016, MW7.8 Kaikōura Earthquake. The most intense landslide damage was concentrated in 3500 km2 around the areas of fault rupture. Given the sparsely populated area affected by landslides, only a few homes were impacted and there were no recorded deaths due to landslides. Landslides caused major disruption with all road and rail links with Kaikōura being severed. The landslides affecting State Highway 1 (the main road link in the South Island of New Zealand) and the South Island main trunk railway extended from Ward in Marlborough all the way to the south of Oaro in North Canterbury. The majority of landslides occurred in two geological and geotechnically distinct materials reflective of the dominant rock types in the affected area. In the Neogene sedimentary rocks (sandstones, limestones and siltstones) of the Hurunui District, North Canterbury and around Cape Campbell in Marlborough, first-time and reactivated rock-slides and rock-block slides were the dominant landslide type. These rocks also tend to have rock material strength values in the range of 5-20 MPa. In the Torlesse 'basement' rocks (greywacke sandstones and argillite) of the Kaikōura Ranges, first-time rock and debris avalanches were the dominant landslide type. These rocks tend to have material strength values in the range of 20-50 MPa. A feature of this earthquake is the large number (more than 200) of valley blocking landslides it generated. This was partly due to the steep and confined slopes in the area and the widely distributed strong ground shaking. The largest landslide dam has an approximate volume of 12(±2) M m3 and the debris from this travelled about 2.7 km2 downslope where it formed a dam blocking the Hapuku River. The long-term stability of cracked slopes and landslide dams from future strong earthquakes and large rainstorms are an ongoing concern to central and local government agencies responsible for rebuilding homes and infrastructure. A particular concern is the potential for debris floods to affect downstream assets and infrastructure should some of the landslide dams breach catastrophically. At least twenty-one faults ruptured to the ground surface or sea floor, with these surface ruptures extending from the Emu Plain in North Canterbury to offshore of Cape Campbell in Marlborough. The mapped landslide distribution reflects the complexity of the earthquake rupture. Landslides are distributed across a broad area of intense ground shaking reflective of the elongate area affected by fault rupture, and are not clustered around the earthquake epicentre. The largest landslides triggered by the earthquake are located either on or adjacent to faults that ruptured to the ground surface. Surface faults may provide a plane of weakness or hydrological discontinuity and adversely oriented surface faults may be indicative of the location of future large landslides. Their location appears to have a strong structural geological control. Initial results from our landslide investigations suggest predictive models relying only on ground-shaking estimates underestimate the number and size of the largest landslides that occurred.
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.
Advanced seismic effective-stress analysis is used to scrutinize the liquefaction performance of 55 well-documented case-history sites from Christchurch. The performance of these sites during the 2010-2011 Canterbury earthquake sequence varied significantly, from no liquefaction manifestation at the ground surface (in any of the major events) to severe liquefaction manifestation in multiple events. For the majority of the 55 sites, the simplified liquefaction evaluation procedures, which are conventionally used in engineering practice, could not explain these dramatic differences in the manifestation. Detailed geotechnical characterization and subsequent examination of the soil profile characteristics of the 55 sites identified some similarities but also important differences between sites that manifested liquefaction in the two major events of the sequence (YY-sites) and sites that did not manifest liquefaction in either event (NN-sites). In particular, while the YY-sites and NN-sites are shown to have practically identical critical layer characteristics, they have significant differences with regard to their deposit characteristics including the thickness and vertical continuity of their critical zones and liquefiable materials. A CPT-based effective stress analysis procedure is developed and implemented for the analyses of the 55 case history sites. Key features of this procedure are that, on the one hand, it can be fully automated in a programming environment and, on the other hand, it is directly equivalent (in the definition of cyclic resistance and required input data) to the CPT-based simplified liquefaction evaluation procedures. These features facilitate significantly the application of effective-stress analysis for simple 1D free-field soil-column problems and also provide a basis for rigorous comparisons of the outcomes of effective-stress analyses and simplified procedures. Input motions for the analyses are derived using selected (reference) recordings from the two major events of the 2010-2011 Canterbury earthquake sequence. A step-by-step procedure for the selection of representative reference motions for each site and their subsequent treatment (i.e. deconvolution and scaling) is presented. The focus of the proposed procedure is to address key aspects of spatial variability of ground motion in the near-source region of an earthquake including extended-source effects, path effects, and variation in the deeper regional geology.
In this paper, we consider how religious leaders and Civil Defence authorities might collaborate to establish a two-way information conduit during the aftermath of a disaster. Using surveys and in-person interviews, clergy in different Christian denominations were asked about their roles in the earthquake, the needs of their congregations and the possibilities and obstacles to deeper collaboration with Civil Defence authorities.
This report presents an overview of the soil profile characteristics at a number of strong motion station (SMS) sites in Christchurch and its surrounds. An extensive database of ground motion records has been captured by the SMS network in the Canterbury region during the Canterbury earthquake sequence. However in order to comprehensively understand the ground motions recorded at these sites and to be able to relate these motions to other locations, a detailed understanding of the shallow geotechnical profile at each SMS is required. The original NZS1170.5 (SNZ 2004) site subsoil classifications for each SMS site is based on regional geological information and well logs located at varying distances from the site. Given the variability of Christchurch soils, more detailed investigations are required in close vicinity to each SMS to better understand stratigraphy and soil properties, which are important in seismic site response. In this regard, CPT, SPT and borehole data, shear wave velocity (Vs) profiles, and horizontal to vertical spectral ratio measurements (H/V) in close vicinity to the SMS were used to develop representative soil profiles at each site. NZS1170.5 (SNZ 2004) site subsoil classifications were updated using Vs and SPT N60 criteria. Site class E boundaries were treated as a sliding scale rather than as a discrete boundary to account for locations with similar site effects potential, an approach which was shown to result in a better delineation between the site classes. SPT N60 values often indicate a stiffer site class than the Vs data for softer soil sites, highlighting the disparity between the two site investigation techniques. Both SPT N60 and Vs based site classes did not always agree with the original site classifications. This emphasises the importance of having detailed site‐specific information at SMS locations in order to properly classify them. Furthermore, additional studies are required to harmonize site classification based on SPT N60 and Vs. Liquefaction triggering assessments were carried out for the Darfield and Christchurch earthquakes, and compared against observed liquefaction surface manifestations and ground motions characteristics at each SMS. In general, the characteristics of the recorded ground motions at each site correlate well with the triggering analyses. However, at sites that likely liquefied at depth (as indicated by triggering analyses and/or inferred from the characteristics of the recorded surface acceleration time series), the presence of a non‐liquefiable crust layer at many of the SMS locations prevented the manifestation of any surface effects.
One of the great challenges facing human systems today is how to prepare for, manage, and adapt successfully to the profound and rapid changes wreaked by disasters. Wellington, New Zealand, is a capital city at significant risk of devastating earthquake and tsunami, potentially requiring mass evacuations with little or short notice. Subsequent hardship and suffering due to widespread property damage and infrastructure failure could cause large areas of the Wellington Region to become uninhabitable for weeks to months. Previous research has shown that positive health and well-being are associated with disaster-resilient outcomes. Preventing adverse outcomes before disaster strikes, through developing strengths-based skill sets in health-protective attitudes and behaviours, is increasingly advocated in disaster research, practise, and management. This study hypothesised that well-being constructs involving an affective heuristic play vital roles in pathways to resilience as proximal determinants of health-protective behaviours. Specifically, this study examined the importance of health-related quality of life and subjective well-being in motivating evacuation preparedness, measured in a community sample (n=695) drawn from the general adult population of Wellington’s isolated eastern suburbs. Using a quantitative epidemiological approach, the study measured the prevalence of key quality of life indicators (physical and mental health, emotional well-being or “Sense of Coherence”, spiritual well-being, social well-being, and life satisfaction) using validated psychometric scales; analysed the strengths of association between these indicators and the level of evacuation preparedness at categorical and continuous levels of measurement; and tested the predictive power of the model to explain the variance in evacuation preparedness activity. This is the first study known to examine multi-dimensional positive health and global well-being as resilient processes for engaging in evacuation preparedness behaviour. A cross-sectional study design and quantitative survey were used to collect self-report data on the study variables; a postal questionnaire was fielded between November 2008 and March 2009 to a sampling frame developed through multi-stage cluster randomisation. The survey response rate was 28.5%, yielding a margin of error of +/- 3.8% with 95% confidence and 80% statistical power to detect a true correlation coefficient of 0.11 or greater. In addition to the primary study variables, data were collected on demographic and ancillary variables relating to contextual factors in the physical environment (risk perception of physical and personal vulnerability to disaster) and the social environment (through the construct of self-determination), and other measures of disaster preparedness. These data are reserved for future analyses. Results of correlational and regression analyses for the primary study variables show that Wellingtonians are highly individualistic in how their well-being influences their preparedness, and a majority are taking inadequate action to build their resilience to future disaster from earthquake- or tsunami-triggered evacuation. At a population level, the conceptual multi-dimensional model of health-related quality of life and global well-being tested in this study shows a positive association with evacuation preparedness at statistically significant levels. However, it must be emphasised that the strength of this relationship is weak, accounting for only 5-7% of the variability in evacuation preparedness. No single dimension of health-related quality of life or well-being stands out as a strong predictor of preparedness. The strongest associations for preparedness are in a positive direction for spiritual well-being, emotional well-being, and life satisfaction; all involve a sense of existential meaningfulness. Spiritual well-being is the only quality of life variable making a statistically significant unique contribution to explaining the variance observed in the regression models. Physical health status is weakly associated with preparedness in a negative direction at a continuous level of measurement. No association was found at statistically significant levels for mental health status and social well-being. These findings indicate that engaging in evacuation preparedness is a very complex, holistic, yet individualised decision-making process, and likely involves highly subjective considerations for what is personally relevant. Gender is not a factor. Those 18-24 years of age are least likely to prepare and evacuation preparedness increases with age. Multidimensional health and global well-being are important constructs to consider in disaster resilience for both pre-event and post-event timeframes. This work indicates a need for promoting self-management of risk and building resilience by incorporating a sense of personal meaning and importance into preparedness actions, and for future research into further understanding preparedness motivations.
The Leader Fault was one of at least 17 faults that ruptured the ground surface across the northeastern South Island of New Zealand during the Mw 7.8 2016 Kaikōura Earthquake. The southern ~6 km of the Leader Fault, here referred to as the South Leader Fault (SLF), ruptured the North Canterbury (tectonic) Domain and is the primary focus of this study. The main objective of the thesis is to understand the key factors that contributed to the geometry and kinematics of the 2016 SLF rupture and its intersection with The Humps Fault (HF). This thesis employs a combination of techniques to achieve the primary objective, including detailed mapping of the bedrock geology, geomorphology and 2016 rupture, measurement of 2016 ground surface displacements, kinematic analysis of slip vectors from the earthquake, and logging of a single natural exposure across a 2016 rupture that was treated as a paleoseismic trench. The resulting datasets were collected in the field, from terrestrial LiDAR and InSAR imagery, and from historical (pre-earthquake) aerial photographs for a ~11 km2 study area. Surface ruptures in the study area are a miniature version of the entire rupture from the earthquake; they are geometrically and kinematically complex, with many individual and discontinuous segments of varying orientations and slip senses which are distributed across a zone up to ~3.5 km wide. Despite this variability, three main groups of ruptures have been identified. These are: 1) NE-SW striking, shallow to moderate dipping (25-45°W) faults that are approximately parallel to Cenozoic bedding with mainly reverse dip-slip and, and for the purposes of this thesis, are considered to be part of the SLF. 2) N-S striking, steeply dipping (~85°E) oblique sinistral faults that are up to the west and part of the SLF. 3) E-NE striking, moderate to steeply dipping (45-68°N) dextral reverse faults which are part of the HF. Bedding-parallel faults are interpreted to be flexural slip structures formed during folding of the near-surface Cenozoic strata, while the steeply dipping SLF ruptured a pre-existing bedrock fault which has little topographic expression. Groups 1 and 2 faults were both locally used for gravitational failure during the earthquake. Despite this non-tectonic fault movement, the slip vectors for faults that ruptured during the earthquake are broadly consistent with NCD tectonics and the regional ~100-120° trend of the principal horizontal stress/strain axes. Previous earthquake activity on the SLF is required by its displacement of Cenozoic formations but Late Quaternary slip on the fault prior to 2016 is neither supported by pre-existing fault scarps nor by changes in topography across the fault. By contrast, at least two earthquakes (including 2016) appear to have ruptured the HF from the mid Holocene, consistent with recurrence intervals of no more than ~7 kyr, and with preliminary observations from trenches on the fault farther to the west. The disparity in paleoearthquake records of the two faults suggests that they typically do not rupture together, thus it is concluded that the HF-SLF rupture pattern observed in the Kaikōura Earthquake rarely occurs in a single earthquake.
Spatial variations in river facies exerted a strong influence on the distribution of liquefaction features observed in Christchurch during the 2010-11 Canterbury Earthquake Sequence (CES). Liquefaction and liquefaction-induced ground deformation was primarily concentrated near modern waterways and areas underlain by Holocene fluvial deposits with shallow water tables (< 1 to 2 m). In southern Christchurch, spatial variations of liquefaction and subsidence were documented in the suburbs within inner meander loops of the Heathcote River. Newly acquired geospatial data, geotechnical reports and eye-witness discussions are compiled to provide a detailed account of the surficial effects of CES liquefaction and ground deformation adjacent to the Heathcote River. LiDAR data and aerial photography are used to produce a new series of original figures which reveal the locations of recurrent liquefaction and subsidence. To investigate why variable liquefaction patterns occurred, the distribution of surface ejecta and associated ground damage is compared with near-surface sedimentologic, topographic, and geomorphic variability to seek relationships between the near-surface properties and observed ground damages. The most severe liquefaction was concentrated within a topographic low in the suburb of St Martins, an inner meander loop of the Heathcote River, with liquefaction only minor or absent in the surrounding areas. Subsurface investigations at two sites in St Martins enable documentation of fluvial stratigraphy, the expressions of liquefaction, and identification of pre-CES liquefaction features. Excavation to water table depths (~1.5 m below the surface) across sand boils reveals multiple generations of CES liquefaction dikes and sills that cross-cut Holocene fluvial and anthropogenic stratigraphy. Based on in situ geotechnical tests (CPT) indicating sediment with a factor of safety < 1, the majority of surface ejecta was sourced from well-sorted fine to medium sand at < 5 m depth, with the most damaging liquefaction corresponding with the location of a low-lying sandy paleochannel, a remnant river channel from the Holocene migration of the meander in St Martins. In the adjacent suburb of Beckenham, where migration of the Heathcote River has been laterally confined by topography associated with the volcanic lithologies of Banks Peninsula, severe liquefaction was absent with only minor sand boils occurring closest to the modern river channel. Auger sampling across the suburb revealed thick (>1 m) clay-rich overbank and back swamp sediments that produced a stratigraphy which likely confined the units susceptible to liquefaction and prevented widespread ejection of liquefied material. This analysis suggests river migration promotes the formation and preservation of fluvial deposits prone to liquefaction. Trenching revealed the strongest CES earthquakes with large vertical accelerations favoured sill formation and severe subsidence at highly susceptible locations corresponding with an abandoned channel. Less vulnerable sites containing deeper and thinner sand bodies only liquefied in the strongest and most proximal earthquakes forming minor localised liquefaction features. Liquefaction was less prominent and severe subsidence was absent where lateral confinement of a Heathcote meander has promoted the formation of fluvial stratum resistant to liquefaction. Correlating CES liquefaction with geomorphic interpretations of Christchurch’s Heathcote River highlights methods in which the performance of liquefaction susceptibility models can be improved. These include developing a reliable proxy for estimating soil conditions in meandering fluvial systems by interpreting the geology and geomorphology, derived from LiDAR data and modern river morphology, to improve the methods of accounting for the susceptibility of an area. Combining geomorphic interpretations with geotechnical data can be applied elsewhere to identify regional liquefaction susceptibilities, improve existing liquefaction susceptibility datasets, and predict future earthquake damage.
Coastal and river environments are exposed to a number of natural hazards that have the potential to negatively affect both human and natural environments. The purpose of this research is to explain that significant vulnerabilities to seismic hazards exist within coastal and river environments and that coasts and rivers, past and present, have played as significant a role as seismic, engineering or socio-economic factors in determining the impacts and recovery patterns of a city following a seismic hazard event. An interdisciplinary approach was used to investigate the vulnerability of coastal and river areas in the city of Christchurch, New Zealand, following the Canterbury Earthquake Sequence, which began on the 4th of September 2010. This information was used to identify the characteristics of coasts and rivers that make them more susceptible to earthquake induced hazards including liquefaction, lateral spreading, flooding, landslides and rock falls. The findings of this research are applicable to similar coastal and river environments elsewhere in the world where seismic hazards are also of significant concern. An interdisciplinary approach was used to document and analyse the coastal and river related effects of the Canterbury earthquake sequence on Christchurch city in order to derive transferable lessons that can be used to design less vulnerable urban communities and help to predict seismic vulnerabilities in other New Zealand and international urban coastal and river environments for the future. Methods used to document past and present features and earthquake impacts on coasts and rivers in Christchurch included using maps derived from Geographical Information Systems (GIS), photographs, analysis of interviews from coastal, river and engineering experts, and analysis of secondary data on seismicity, liquefaction potential, geology, and planning statutes. The Canterbury earthquake sequence had a significant effect on Christchurch, particularly around rivers and the coast. This was due to the susceptibility of rivers to lateral spreading and the susceptibility of the eastern Christchurch and estuarine environments to liquefaction. The collapse of river banks and the extensive cracking, tilting and subsidence that accompanied liquefaction, lateral spreading and rock falls caused damage to homes, roads, bridges and lifelines. This consequently blocked transportation routes, interrupted electricity and water lines, and damaged structures built in their path. This study found that there are a number of physical features of coastal and river environments from the past and the present that have induced vulnerabilities to earthquake hazards. The types of sediments found beneath eastern Christchurch are unconsolidated fine sands, silts, peats and gravels. Together with the high water tables located beneath the city, these deposits made the area particularly susceptible to liquefaction and liquefaction-induced lateral spreading, when an earthquake of sufficient size shook the ground. It was both past and present coastal and river processes that deposited the types of sediments that are easily liquefied during an earthquake. Eastern Christchurch was once a coastal and marine environment 6000 years ago when the shoreline reached about 6 km inland of its present day location, which deposited fine sand and silts over this area. The region was also exposed to large braided rivers and smaller spring fed rivers, both of which have laid down further fine sediments over the following thousands of years. A significant finding of this study is the recognition that the Canterbury earthquake sequence has exacerbated existing coastal and river hazards and that assessments and monitoring of these changes will be an important component of Christchurch’s future resilience to natural hazards. In addition, patterns of recovery following the Canterbury earthquakes are highlighted to show that coasts and rivers are again vulnerable to earthquakes through their ability to recovery. This city’s capacity to incorporate resilience into the recovery efforts is also highlighted in this study. Coastal and river areas have underlying physical characteristics that make them increasingly vulnerable to the effects of earthquake hazards, which have not typically been perceived as a ‘coastal’ or ‘river’ hazard. These findings enhance scientific and management understanding of the effects that earthquakes can have on coastal and river environments, an area of research that has had modest consideration to date. This understanding is important from a coastal and river hazard management perspective as concerns for increased human development around coastlines and river margins, with a high seismic risk, continue to grow.
Documenting earthquake-induced ground deformation is significant to assess the characteristics of past and contemporary earthquakes and provide insight into seismic hazard. This study uses airborne light detection and ranging (LiDAR) and conducts multi-disciplinary field techniques to document the surface rupture morphology and evaluate the paleoseismicity and seismic hazard parameters of the Hurunui segment of the Hope Fault in the northern South Island of New Zealand. It also documents and evaluates seismically induced features and ground motion characteristics of the 2010 Darfield and 2011 Christchurch earthquakes in the Port Hills, south of Christchurch. These two studies are linked in that they investigate the near-field coseismic features of large (Mw ~7.1) earthquakes in New Zealand and produce data for evaluating seismic hazards of future earthquakes. In the northern South Island of New Zealand, the Australian-Pacific plate boundary is characterised by strike-slip deformation across the Marlborough Fault System (MFS). The ENE-striking Hope Fault (length: ~230 km) is the youngest and southernmost fault in the MFS, and the second fastest slipping fault in New Zealand. The Hope Fault is a major source of seismic hazard in New Zealand and has ruptured (in-part) historically in the Mw 7.1 1888 Amuri earthquake. In the west, the Hurunui segment of the Hope Fault is covered by beech forest. Hence, its seismic hazard parameters and paleoearthquake chronology were poorly constrained and it was unknown whether the 1888 earthquake ruptured this segment or not and if so, to what extent. Utilising LiDAR and field data, a 29 km-long section of the Hurunui segment of the Hope Fault is mapped. LiDAR-mapping clearly reveals the principal slip zone (PSZ) of the fault and a suite of previously unrecognised structures that form the fault deformation zone (FDZ). FDZ width measurements from 415 locations reveal a spatially-variable, active FDZ up to ~500 m wide with an average width of 200 m. Kinematic analysis of the fault structures shows that the Hurunui segment strikes between 070° and 075° and is optimally oriented for dextral strike-slip within the regional stress field. This implies that the wide FDZ observed is unlikely to result from large-scale fault mis-orientation with respect to regional stresses. The analysis of FDZ width indicates that it increases with increased hanging wall topography and increased topographic relief suggesting that along-strike topographic perturbations to fault geometry and stress states increase fault zone complexity and width. FDZ width also increases where the tips of adjacent PSZ strands locally vary in strike, and where the thickness of alluvial deposits overlying bedrock increases. LiDAR- and photogrammetrically-derived topographic mapping indicates that the boundary between the Hurunui and Hope River segments is characterised by a ~850-m-wide right stepover and a 9º-14° fault bend. Paleoseismic trenching at Hope Shelter site reveals that 6 earthquakes occurred at A.D. 1888, 1740-1840, 1479-1623, 819-1092, 439-551, and 373- 419. These rupture events have a mean recurrence interval of ~298 ± 88 yr and inter-event times ranging from 98 to 595 yrs. The variation in the inter-event times is explained by (1) coalescing rupture overlap from the adjacent Hope River segment on to the Hurunui segment at the study site, (2) temporal clustering of large earthquakes on the Hurunui segment, and/or (3) ‘missing’ rupture events. It appears that the first two options are more plausible to explain the earthquake chronologies and rupture behaviour on the Hurunui segment, given the detailed nature of the geologic and chronologic investigations. This study provides first evidence for coseismic multi-segment ruptures on the Hope Fault by identifying a rupture length of 44-70 km for the 1888 earthquake, which was not confined to the Hope River segment (primary source for the 1888 earthquake). LiDAR data is also used to identify and measure dextral displacements and scarp heights from the PSZ and structures within the FDZ along the Hurunui segment. Reconstruction of large dextrally-offset geomorphic features shows that the vertical component of slip accounts for only ~1% of the horizontal displacements and confirms that the fault is predominantly strike-slip. A strong correlation exists between the dextral displacements and elevations of geomorphic features suggesting the possibility of age correlation between the geomorphic features. A mean single event displacement (SED) of 3.6 ± 0.7 m is determined from interpretation of sets of dextral displacements of ≤ 25 m. Using the available surface age data and the cumulative dextral displacements from Matagouri Flat, McKenzie Fan, Macs Knob and Hope River sites, and the mean SED, a mean slip rate of 12.2 ± 2.4 mm/yr, and a mean recurrence interval of ~320 ± 120 yr, and a potential earthquake magnitude of Mw 7.2 are determined for the Hurunui segment. This study suggests that the fault slip rate has been constant over the last ~15000 yr. Strong ground motions from the 2010 Darfield (Canterbury) earthquake displaced boulders and caused ground damage on some ridge crests in the Port Hills. However, the 2011 Christchurch earthquake neither displaced boulders nor caused ground damage at the same ridge crests. Documentation of locations (~400 m a.s.l.), lateral displacements (8-970 cm), displacement direction (250° ± 20°) of displaced boulders, in addition to their hosting socket geometries (< 1 cm to 50 cm depth), the orientation of the ridges (000°-015°) indicate that boulders have been displaced in the direction of instrumentally recorded transient peak ground horizontal displacements nearby and that the seismic waves have been amplified at the study sites. The co-existence of displaced and non-displaced boulders at proximal sites suggests small-scale ground motion variability and/or varying boulder-ground dynamic interactions relating to shallow phenomena such as variability in soil depth, bedrock fracture density and/or microtopography on the bedrock-soil interface. Shorter shaking duration of the 2011 Christchurch event, differing frequency contents and different source characteristics were all factors that may have contributed to generating circumstances less favourable to boulder displacement in this earthquake. Investigating seismically induced features, fault behaviour, site effects on the rupture behaviour, and site response to the seismic waves provides insights into fault rupture hazards.
