We present initial results from a set of three-dimensional (3D) deterministic earthquake ground motion simulations for the northern Canterbury plains, Christchurch and the Banks Peninsula region, which explicitly incorporate the effects of the surface topography. The simu-lations are done using Hercules, an octree-based finite-element parallel software for solving 3D seismic wave propagation problems in heterogeneous media under kinematic faulting. We describe the efforts undertaken to couple Hercules with the South Island Velocity Model (SIVM), which included changes to the SIVM code in order to allow for single repetitive que-ries and thus achieve a seamless finite-element meshing process within the end-to-end ap-proach adopted in Hercules. We present our selection of the region of interest, which corre-sponds to an area of about 120 km × 120 km, with the 3D model reaching a depth of 60 km. Initial simulation parameters are set for relatively high minimum shear wave velocity and a low maximum frequency, which we are progressively scaling up as computing resources permit. While the effects of topography are typically more important at higher frequencies and low seismic velocities, even at this initial stage of our efforts (with a maximum of 2 Hz and a mini-mum of 500 m/s), it is possible to observe the importance of the topography in the response of some key locations within our model. To highlight these effects we compare the results of the 3D topographic model with respect to those of a flat (squashed) 3D model. We draw rele-vant conclusions from the study of topographic effects during earthquakes for this region and describe our plans for future work.
Small, tight-knit communities, are complex to manage from outside during a disaster. The township of Lyttelton, New Zealand, and the communities of Corsair Bay, Cass Bay, and Rapaki to the east, are especially more so difficult due to the terrain that encloses them, which caused them to be cut-off from Christchurch, the largest city in the South Island, barely 10 km away, after the Mw 7.1 Darfield Earthquake and subsequent Canterbury Earthquake Sequence. Lyttelton has a very strong and deep-rooted community spirit that draws people to want to be a part of Lyttelton life. It is predominantly residential on the slopes, with retail space, service and light industry nestled near the harbour. It has heritage buildings stretching back to the very foundation of Canterbury yet hosts the largest, modern deep-water port for the region. This study contains two surveys: one circulated shortly before the Darfield Earthquake and one circulated in July 2011, after the Christchurch and Sumner Earthquakes. An analytical comparison of the participants’ household preparedness for disaster before the Darfield Earthquake and after the Christchurch and Sumner Earthquakes was performed. A population spatiotemporal distribution map was produced that shows the population in three-hourly increments over a week to inform exposure to vulnerability to natural hazards. The study went on to analyse the responses of the participants in the immediate period following the Chrsitchurch and Sumner Earthquakes, including their homeward and subsequent journeys, and the decision to evacuate or stay in their homes. Possible predictors to a decision to evacuate some or all members of the household were tested. The study also asked participants’ views on the events since September 2010 for analysis.
On November 14 2016 a magnitude 7.8 earthquake struck the south island of New Zealand. The earthquake lasted for just two minutes with severe seismic shaking and damage in the Hurunui and Kaikōura districts. Although these are predominantly rural areas, with scattered small towns and mountainous topography, they also contain road and rail routes that are essential parts of the national transport infrastructure. This earthquake and the subsequent recovery are of particular significance as they represent a disaster following in close proximity to another similar disaster, with the Canterbury earthquakes occurring in a neighboring district five years earlier. The research used an inductive qualitative case study to explore the nature of the Kaikōura recovery. That recovery process involved a complex interplay between the three parties; (a) the existing local government in the district, (b) central government agencies funding the recovery of the local residents and the national transport infrastructure, and (c) recovery leaders arriving with recent expertise from the earlier Canterbury disaster. It was evident that three groups: locals, government, and experts represented a multi-party governance debate in which the control of the Kaikōura earthquake recovery was shared amongst them. Each party had their own expertise, adgenda and networks that they brought to the Kaikōura recovery, but this created tensions between external expertise and local, community leadership. Recent earthquake research suggests that New Zealand is currently in the midst of an earthquake cluster, with further seismic disasters likely to occur in relatively close succession. This is likely to be compounded by the increasing frequency of other natural disasters with the effects of climate change. The present study investigates a phenomenon that may become increasingly common, with the transfer of disaster expertise from one event to another, and the interface between those experts with local and national government in directing recoveries. The findings of this study have implications for practitioners and policy makers in NZ and other countries where disasters are experienced in close spatial and temporal proximity.
This poster presents work to date on ground motion simulation validation and inversion for the Canterbury, New Zealand region. Recent developments have focused on the collection of different earthquake sources and the verification of the SPECFEM3D software package in forward and inverse simulations. SPECFEM3D is an open source software package which simulates seismic wave propagation and performs adjoint tomography based upon the spectral-element method. Figure 2: Fence diagrams of shear wave velocities highlighting the salient features of the (a) 1D Canterbury velocity model, and (b) 3D Canterbury velocity model. Figure 5: Seismic sources and strong motion stations in the South Island of New Zealand, and corresponding ray paths of observed ground motions. Figure 3: Domain used for the 19th October 2010 Mw 4.8 case study event including the location of the seismic source and strong motion stations. By understanding the predictive and inversion capabilities of SPECFEM3D, the current 3D Canterbury Velocity Model can be iteratively improved to better predict the observed ground motions. This is achieved by minimizing the misfit between observed and simulated ground motions using the built-in optimization algorithm. Figure 1 shows the Canterbury Velocity Model domain considered including the locations of small-to-moderate Mw events [3-4.5], strong motion stations, and ray paths of observed ground motions. The area covered by the ray paths essentially indicates the area of the model which will be most affected by the waveform inversion. The seismic sources used in the ground motion simulations are centroid moment tensor solutions obtained from GeoNet. All earthquake ruptures are modelled as point sources with a Gaussian source time function. The minimum Mw limit is enforced to ensure good signal-to-noise ratio and well constrained source parameters. The maximum Mw limit is enforced to ensure the point source approximation is valid and to minimize off-fault nonlinear effects.
On 22 February 2011，a magnitude Mw 6.3 earthquake occurred with an epicenter located near Lyttelton at about 10km from Christchurch in Canterbury region on the South Island of New Zealand (Figure 1). Since this earthquake occurred in the midst of the aftershock activity which had continued since the 4 September 2010 Darfield Earthquake occurrence, it was considered to be an aftershock of the initial earthquake. Because of the short distance to the city and the shallower depth of the epicenter, this earthquake caused more significant damage to pipelines, traffic facilities, residential houses/properties and multi-story buildings in the central business district than the September 2010 Darfield Earthquake in spite of its smaller earthquake magnitude. Unfortunately, this earthquake resulted in significant number of casualties due to the collapse of multi-story buildings and unreinforced masonry structures in the city center of Christchurch. As of 4 April, 172 casualties were reported and the final death toll is expected to be 181. While it is extremely regrettable that Christchurch suffered a terrible number of victims, civil and geotechnical engineers have this hard-to-find opportunity to learn the response of real ground from two gigantic earthquakes which occurred in less than six months from each other. From geotechnical engineering point of view, it is interesting to discuss the widespread liquefaction in natural sediments, repeated liquefaction within short period and further damage to earth structures which have been damaged in the previous earthquake. Following the earthquake, an intensive geotechnical reconnaissance was conducted to capture evidence and perishable data from this event. The team included the following members: Misko Cubrinovski (University of Canterbury, NZ, Team Leader), Susumu Yasuda (Tokyo Denki University, Japan, JGS Team Leader), Rolando Orense (University of Auckland, NZ), Kohji Tokimatsu (Tokyo Institute of Technology, Japan), Ryosuke Uzuoka (Tokushima University, Japan), Takashi Kiyota (University of Tokyo, Japan), Yasuyo Hosono (Toyohashi University of Technology, Japan) and Suguru Yamada (University of Tokyo, Japan).
Recent surface-rupturing earthquakes in New Zealand have highlighted significant exposure and vulnerability of the road network to fault displacement. Understanding fault displacement hazard and its impact on roads is crucial for mitigating risks and enhancing resilience. There is a need for regional-scale assessments of fault displacement to identify vulnerable areas within the road network for the purposes of planning and prioritising site-specific investigations. This thesis employs updated analysis of data from three historical surface-rupturing earthquakes (Edgecumbe 1987, Darfield 2010, and Kaikoūra 2016) to develop an empirical model that addresses the gap in regional fault displacement hazard analysis. The findings contribute to understanding of • How to use seismic hazard model inputs for regional fault displacement hazard analysis • How faulting type and sediment cover affects the magnitude and spatial distribution of fault displacement • How the distribution of displacement and regional fault displacement hazard is impacted by secondary faulting • The inherent uncertainties and limitations associated with employing an empirical approach at a regional scale • Which sections of New Zealand’s roading network are most susceptible to fault displacement hazard and warrant site-specific investigations • Which regions should prioritise updating emergency management plans to account for post-event disruptions to roading. I used displacement data from the aforementioned historical ruptures to generate displacement versus distance-to-fault curves for various displacement components, fault types, and geological characteristics. Using those relationships and established relationships for along-strike displacement, displacement contours were generated surrounding active faults within the NZ Community Fault Model. Next, I calculated a new measure of 1D strain along roads as well as relative hazard, which integrated 1D strain and normalised slip rate data. Summing these values at the regional level identified areas of heightened relative hazard across New Zealand, and permits an assessment of the susceptibility of road networks using geomorphon land classes as proxies for vulnerability. The results reveal that fault-parallel displacements tend to localise near the fault plane, while vertical and fault-perpendicular displacements sustain over extended distances. Notably, no significant disparities were observed in off-fault displacement between the hanging wall and footwall sides of the fault, or among different surface geology types, potentially attributed to dataset biases. The presence of secondary faulting in the dataset contributes to increased levels of tectonic displacement farther from the fault, highlighting its significance in hazard assessments. Furthermore, fault displacement contours delineate broader zones around dip-slip faults compared to strike-slip faults, with correlations identified between fault length and displacement width. Road ‘strain’ values are higher around dip-slip faults, with notable examples observed in the Westland and Buller Districts. As expected, relative hazard analysis revealed elevated values along faults with high slip rates, notably along the Alpine Fault. A regional-scale analysis of hazard and exposure reveals heightened relative hazard in specific regions, including Wellington, Southern Hawke’s Bay, Central Bay of Plenty, Central West Coast, inland Canterbury, and the Wairau Valley of Marlborough. Notably, the Central West Coast exhibits the highest summed relative hazard value, attributed to the fast-slipping Alpine Fault. The South Island generally experiences greater relative hazard due to larger and faster-slipping faults compared to the North Island, despite having fewer roads. Central regions of New Zealand face heightened risk compared to Southern or Northern regions. Critical road links intersecting high-slipping faults, such as State Highways 6, 73, 1, and 2, necessitate prioritisation for site-specific assessments, emergency management planning and targeted mitigation strategies. Roads intersecting with the Alpine Fault are prone to large parallel displacements, requiring post-quake repair efforts. Mitigation strategies include future road avoidance of nearby faults, modification of road fill and surface material, and acknowledgement of inherent risk, leading to prioritised repair efforts of critical roads post-quake. Implementing these strategies enhances emergency response efforts by improving accessibility to isolated regions following a major surface-rupturing event, facilitating faster supply delivery and evacuation assistance. This thesis contributes to the advancement of understanding fault displacement hazard by introducing a novel regional, empirical approach. The methods and findings highlight the importance of further developing such analyses and extending them to other critical infrastructure types exposed to fault displacement hazard in New Zealand. Enhancing our comprehension of the risks associated with fault displacement hazard offers valuable insights into various mitigation strategies for roading infrastructure and informs emergency response planning, thereby enhancing both national and global infrastructure resilience against geological hazards.
This study is a qualitative investigation into the decision-making behaviour of commercial property owners (investors and developers) who are rebuilding in a city centre after a major disaster. In 2010/2011, Christchurch, the largest city in the South Island of New Zealand, was a site of numerous earthquakes. The stronger earthquakes destroyed many buildings and public infrastructure in the commercial inner city. As a result, affected property owners lost all or most of their buildings, a significant proportion of which were old and in the last phase of their life span. They had to negotiate pay-outs with insurance companies and decide, once paid out, whether they should rebuild in Christchurch or sell up and invest elsewhere. The clear majority of those who decided to reinvest in and rebuild the city are ‘locals’, almost all of whom had no prior experience of property development. Thus, in a post-disaster environment, most of these property owners have transitioned from being just being passive investors to active property developers. Their experience was interpreted using primary data gathered from in-depth and semi-structured interviews with twenty-one “informed property people” who included commercial property owners; property agents or consultants; representatives of public-sector agencies and financial institutions. The study findings showed that the decision-making behaviour of property investors and developers rebuilding after a major disaster did not necessarily follow a strict financial or profit motive as prescribed in the mainstream or neo-classical economics property literature. Rather, their decision-making behaviour has been largely shaped by emotional connections and external factors associated with their immediate environment. The theoretical proposition emerging from this study is that after a major disaster, local urban property owners are faced with two choices “to stay” or “to go”. Those who decide to stay and rebuild are typically very committed individuals who have a feeling of ownership, belonging and attachment to the city in which they live and work. These are people who will often take the lead in commercial property development, proactively making decisions and seeking positive investment outcomes for themselves which in turn result in revitalised commercial urban precincts.
Oblique convergence of the Pacific and Australian Plates is accommodated in the northern South Island by the Marlborough Fault System. The Hope Fault is the southern of four major dextral strike-slip faults of this system. Hanmer Basin is a probable segment boundary between the Hope River and Conway segments of the Hope Fault. The Conway segment is transpressional and shows increasing structural complexity near the segment boundary at Hanmer Basin, with multiple Late Quaternary traces, and fault-parallel folding in response to across-fault shortening. Between Hossack Station and Hanmer Basin a crush zone in excess of one kilometre wide is exposed in incised streams and rivers. The crush zone has an asymmetrical geometry about the active trace of the Hope Fault, being only 100-300 metres wide south of the fault, and more than 500 metres wide north of the fault. The most intense deformation of Torlesse bedrock occurs at the south side of the fault zone, indicating that strain is accommodated against the fault footwall. North of the fault deformation is less intense, but occurs over a wider area. The wide fault zone at Hossack Station may reflect divergence of the Hanmer Fault, a major splay of the Hope Fault. At Hossack Station, the Hope Fault has accommodated at least 260 metres of dextral displacement during the Holocene. Dating of abandoned stream channels, offset by the Hope Fault, indicated a Late Holocene dextral slip-rate of 18±8 mm-¹ for the west end of the Conway segment. Using empirical formulae and inferred fault parameters, the expected magnitude of an earthquake generated by the Conway segment is M6.9 to M7.4; for an exceedence probability of 10%, the magnitude is M7.7 to M7.9. Effects associated with coseismic rupture of the Conway segment include shaking of up to MMIX along the ruptured fault and at Hanmer Basin. Uplift at the east end of Hanmer Basin, in conjunction with subsidence at the southwest margin of the basin, is resulting in the development of onlapping stratigraphy. Seismic reflection profiles support this theory. Possible along-fault migration of the basin is inferred to be a consequence of non-parallelism of the master faults.
Wellington is located on a fault line which will inevitably, one day be impacted by a big earthquake. Due to where this fault line geographically sits, the central city and southern suburbs may be cut off from the rest of the region, effectively making these areas an ‘island’. This issue has absorbed a lot of attention, in particular at a large scale by many different fields: civil engineering, architecture, infrastructure planning & design, policymaking. Due to heightened awareness, and evolved school of practice, contemporary landscape architects deal with post-disaster design – Christchurch, NZ has seen this. A number of landscape architects work with nature, following increased application of ecological urbanism, and natural systems thinking, most notably at larger scales. To create parks that are designed to flood, or implement projects to protect shorelines. A form of resilience less often considered is how design for the small scale - people’s 1:1 relationship with their immediate context in exterior space - can be influential in forming a resilient response to the catastrophe of a major earthquake. This thesis intends to provide a response to address the shift of scales, as a paradigm for preparation and recovery. After a large-scale earthquake, state and civic policies and agencies may or subsequentially not go into action. The most important thinking and acting will be what happens in the minds, and the immediate needs, of each and every person; and how they act communally. This is considered in general social terms in state and civic education programmes of civil defence, for example, but much less considered in how the physical design of the actual spaces we inhabit day-to-day can educate us to be mentally prepared to help each other survive a catastrophe. Specifically, the identification of design of typologies can provide these educative functions. Typology inherently a physical form or manipulation of a generic and substantial prototype applicable in contexts is something that exists in the mind. Working with the physical and social appearance and experience of typologies can also/will change people’s minds. Socially, and economically driven, the community-building power of community gardening is well-proven and documented, and a noticeably large part of contemporary landscape architecture. The designs of this thesis will focus on community gardening specifically to form typologies of resilience preparation and response to disaster. The foundation will remain at the small scale of the local community. The specific question this thesis poses: Can we design local typologies in landscape architecture to integrate community gardens, with public space by preparing for and acting as recovery from a disaster?
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 dissertation addresses a diverse range of applied aspects in ground motion simulation validation via the response of complex structures. In particular, the following topics are addressed: (i) the investigation of similarity between recorded and simulated ground motions using code-based 3D irregular structural response analysis, (ii) the development of a framework for ground motion simulations validation to identify the cause of differences between paired observed and simulated dataset, and (iii) the illustration of the process of using simulations for seismic performance-based assessment. The application of simulated ground motions is evaluated for utilisation in engineering practice by considering responses of 3D irregular structures. Validation is performed in a code-based context when the NZS1170.5 (NZS1170.5:2004, 2004) provisions are followed for response history analysis. Two real buildings designed by engineers and physically constructed in Christchurch before the 2010-2011 Canterbury earthquake sequence are considered. The responses are compared when the buildings are subjected to 40 scaled recorded and their subsequent simulated ground motions selected from 22 February 2011 Christchurch. The similarity of recorded and simulated responses is examined using statistical methods such as bootstrapping and hypothesis testing to determine whether the differences are statistically significant. The findings demonstrate the applicability of simulated ground motion when the code-based approach is followed in response history analysis. A conceptual framework is developed to link the differences between the structural response subjected to simulated and recorded ground motions to the differences in their corresponding intensity measures. This framework allows the variability to be partitioned into the proportion that can be “explained” by the differences in ground motion intensity measures and the remaining “unexplained” variability that can be attributed to different complexities such as dynamic phasing of multi-mode response, nonlinearity, and torsion. The application of this framework is examined through a hierarchy of structures reflecting a range of complexity from single-degree-of-freedom to 3D multi-degree-of-freedom systems with different materials, dynamic properties, and structural systems. The study results suggest the areas that ground motion simulation should focus on to improve simulations by prioritising the ground motion intensity measures that most clearly account for the discrepancies in simple to complex structural responses. Three approaches are presented to consider recorded or simulated ground motions within the seismic performance-based assessment framework. Considering the applications of ground motions in hazard and response history analyses, different pathways in utilising ground motions in both areas are explored. Recorded ground motions are drawn from a global database (i.e., NGA-West2 Ancheta et al., 2014). The NZ CyberShake dataset is used to obtain simulations. Advanced ground motion selection techniques (i.e., generalized conditional intensity measure, GCIM) are used for ground motion selection at a few intensity levels. The comparison is performed by investigating the response of an example structure (i.e., 12-storey reinforced concrete special moment frame) located in South Island, NZ. Results are compared and contrasted in terms of hazard, groundmotion selection, structural responses, demand hazard, and collapse risk, then, the probable reasons for differences are discussed. The findings from this study highlight the present opportunities and shortcomings in using simulations in risk assessment. i
This thesis focuses on the role of legal preparedness for managing large-scale urban disasters in Aotearoa New Zealand. It uses the Auckland Volcanic Field as a case study to answer the question: ‘is New Zealand’s current legal framework prepared to respond to and recover from a large-scale urban disaster?’. The Auckland Volcanic Field was chosen as the main case study because a future eruption is a low likelihood, high-impact event that New Zealand is going to have to manage in the future. Case studies are a key feature of this thesis as both New Zealand based and overseas examples are used to explore the role of legal preparedness by identifying and investigating a range of legal issues that need to be addressed in advance of a future Auckland Volcanic Field eruption. Of particular interest is the impact of legal preparedness for the recovery phase. The New Zealand case studies include; Canterbury earthquake sequence 2010-2011, the Kaikōura earthquake 2016, the Auckland flooding 2018, and the North Island Severe Weather event 2023, which encompasses both the Auckland Anniversary weekend flooding and Cyclone Gabrielle. As New Zealand has not experienced a large-scale urban volcanic eruption, overseas examples are explored to provide insights into the legal issues that are volcano specific. The overseas volcanic case studies cover eruptions in Heimaey (Iceland), the Soufrière Hills (Montserrat and the Grenadines), La Soufrière (St Vincent) and Tungurahua (Ecuador). New Zealand’s past experiences highlight a trend for introducing post-event legal frameworks to manage recovery. Consequently, the current disaster management system is not prioritising legal preparedness and instead is choosing to rely on exceptional powers. Unsurprisingly, the introduction of new post-event recovery frameworks has repercussions. In New Zealand, new post-event legal frameworks are introduced swiftly under urgency, they contain broad unstructured decision-making powers, and are often flawed. As these exceptional new frameworks sit outside the ‘normal’ legal frameworks, they in effect create a parallel “shadow system”. Based on the evidence explored in this thesis it does not appear that Auckland’s current disaster management framework is prepared to deal with a large-scale urban event caused by an Auckland Volcanic Field eruption. Following this conclusion, it is the submission of this thesis that New Zealand’s current legal framework is not prepared to respond to and recover from a large-scale urban disaster. To become legally prepared, New Zealand needs to consider the legal tools required to manage large-scale urban disasters in advance. This will prevent the creation of a legal vacuum in the aftermath of disasters and the need for new recovery frameworks. Adopting a new attitude will require a change in approach towards legal preparedness which prioritises it, rather than sidelining it. This may also require changes within New Zealand’s disaster management system including the introduction of a formal monitoring mechanism, which will support and prioritise legal preparedness. This thesis has shown that not legally preparing for future disasters is a choice which carries significant consequences. None of these consequences are inevitable when managing large-scale disasters, however they are inevitable when frameworks are not legally prepared in advance. To not legally prepare, is to prepare to fail and thus create a disaster by choice.
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.
<b>Aotearoa has undoubtedly some of the most beautiful landscapes in the world, a privilege for its inhabitants. However, as our cities have developed post-colonisation, the connection between the natural environment and its occupants has diminished. Designers play a vital role within an ever evolving world to progress the built environment in a way that reflects and restores vital values that have been deprioritised. Future practice should prioritise diversity, care for the land, enhancement of community space, and sustainable practices.</b>
This research sets out to demonstrate that new design methodologies can encourage kaitiakitanga, whilst meeting the needs of urban public space. Initially through critical analysis and literature based research, a study of Ōtautahi Christchurch, the South Island’s largest city, was undertaken. The principles of a ‘15 minute city’ were also explored and applied to the city, establishing issues within the built environment that drove the overall research direction.
Through the tools of critical reflection and a research through design methodology, a design toolkit was constructed. This toolkit sets out to provide designers with a simple streamlined method of developing urban interventions that are sustainable and beneficial for human well-being. The toolkit incorporates an abstraction of the ‘15 minute city’ ideology and introduces the concepts of evolving green transportation routes within cities. Ōtautahi Christchurch, a city with a significant history of earthquake-caused damage, was chosen as the primary site for the application of this research’s proposed toolkit. The city becomes a canvas for an urban rebuild that explores and aims to set a precedent for a progressive 21st-century city.
A key finding as the toolkit research developed was the idea of a ‘temporary’ phase or intervention, being added to traditional design methodologies prior to permanent building. The research explains how this temporary phase could more actively engage diverse user groups and create active conversations between communities and designers.
The refined toolkit sets outs proposed timeline phases, methods of site analysis and development of design drivers. Alongside this, a modular architectural system establishes a design proposal for the temporary phase of an individual site within an evolving green route. This outcome provides further opportunity for realistic testing, which would actively involve communities and aims to shift our priorities within urban development. The introduction of the ‘temporary’ phase is beneficial in mitigating psychological implications on people and limiting physical impacts on the landscape.
The final design stage of the thesis applied the toolkit process to three sites in Ōtautahi Christchurch. Through a holistic lens, the toolkit framework set out methods to collate information that provides guidance for development on the sites. While some layers are initiated simply by recognising site characteristics, others are informed through software such as GIS.
Connected by a proposed green transport route, the three initial sites are developed with temporary interventions that utilise the modular design set out previously in the research. Contextualising the interventions on real world sites tested the flexibility of the system and allowed for critical reflection on the applicability of the toolkit to Aotearoa.
The research concludes by identifying future research opportunities and speculates on possible applications of its findings within the real world. Temporary Permanence highlights the significant role that we, as designers, have in shifting urban priorities to create more holistic, sustainable, and inclusive cities for people and the planet.
This dissertation addresses a diverse range of topics in the physics-based broadband ground motion simulation, with a focus on New Zealand applications. In particular the following topics are addressed: the methodology and computational implementation of a New Zealand Velocity Model for broadband ground motion simulation; generalised parametric functions and spatial correlations for seismic velocities in the Canterbury, New Zealand region from surface-wave-based site characterisation; and ground motion simulations of Hope Fault earthquakes. The paragraphs below outline each contribution in more detail. A necessary component in physics-based ground motion simulation is a 3D model which details the seismic velocities in the region of interest. Here a velocity model construction methodology, its computational implementation, and application in the construction of a New Zealand velocity model for use in physics-based broadband ground motion simulation are presented. The methodology utilises multiple datasets spanning different length scales, which is enabled via the use of modular sub-regions, geologic surfaces, and parametric representations of crustal velocity. A number of efficiency-related workflows to decrease the overall computational construction time are employed, while maintaining the flexibility and extensibility to incorporate additional datasets and re- fined velocity parameterizations as they become available. The model comprises explicit representations of the Canterbury, Wellington, Nelson-Tasman, Kaikoura, Marlborough, Waiau, Hanmer and Cheviot sedimentary basins embedded within a regional travel-time tomography-based velocity model for the shallow crust and provides the means to conduct ground motion simulations throughout New Zealand for the first time. Recently developed deep shear-wave velocity profiles in Canterbury enabled models that better characterise the velocity structure within geologic layers of the Canterbury sedimentary basin to be developed. Here the development of depth- and Vs30-dependent para-metric velocity and spatial correlation models to characterise shear-wave velocities within the geologic layers of the Canterbury sedimentary basin are presented. The models utilise data from 22 shear-wave velocity profiles of up to 2.5km depth (derived from surface wave analysis) juxtaposed with models which detail the three-dimensional structure of the geologic formations in the Canterbury sedimentary basin. Parametric velocity equations are presented for Fine Grained Sediments, Gravels, and Tertiary layer groupings. Spatial correlations were developed and applied to generate three-dimensional stochastic velocity perturbations. Collectively, these models enable seismic velocities to be realistically represented for applications such as 3D ground motion and site response simulations. Lastly the New Zealand velocity model is applied to simulate ground motions for a Mw7.51 rupture of the Hope Fault using a physics-based simulation methodology and a 3D crustal velocity model of New Zealand. The simulation methodology was validated for use in the region through comparison with observations for a suite of historic small magnitude earthquakes located proximal to the Hope Fault. Simulations are compared with conventionally utilised empirical ground motion models, with simulated peak ground velocities being notably higher in regions with modelled sedimentary basins. A sensitivity analysis was undertaken where the source characteristics of magnitude, stress parameter, hypocentre location and kinematic slip distribution were varied and an analysis of their effect on ground motion intensities is presented. It was found that the magnitude and stress parameter strongly influenced long and short period ground motion amplitudes, respectively. Ground motion intensities for the Hope Fault scenario are compared with the 2016 Kaikoura Mw7.8 earthquake, it was found that the Kaikoura earthquake produced stronger motions along the eastern South Island, while the Hope Fault scenario resulted in stronger motions immediately West of the near-fault region. The simulated ground motions for this scenario complement prior empirically-based estimates and are informative for mitigation and emergency planning purposes.
Structures of the Lowry Peaks Range - Waikari Valley district are complex. The majority comprise three members of a predominantly WSW -ENE striking major northwards-directed, leading edge imbricate thrust system, with associated angular, asymmetric fault-propagation folds. This system forms anomalously within a large NESW trending belt of structures characterising the entire east coast of north Canterbury, both onshore and offshore and terminates westwards against N-S striking, east facing fold-fault zone. The objectives of this study address the origin, geometry and kinematics of the interaction between these diversely trending systems. Stratigraphy and small-scale structures denote three periods of deformation, namely: i) Middle Cretaceous deformation of the basement rocks, ii) weak Middle Oligocene deformation associated with the inception of the plate boundary through the South Island, and iii) major Pliocene - Recent deformation that formed the majority of the above-mentioned structures. Stress tensor analyses within competent basement and limestone cover rocks suggest two sets of sub-horizontal compression, NE-SW and NW-SE, the former likely to relate to a localised earlier period of deformation, now overprinted by the latter. NW-SE oriented sub-horizontal compression correlates well with results from other parts of north Canterbury. The result of NW-SE compression on the W-E to WSW-ENE striking structures is a large component of oblique motion, which is manifest in four ways: i) movement on two, differently oriented splays rather than a single fault strand, ii) the development of a sinuous trace for a number of the major folds, whereby the ends are oriented normal to the compression direction, the centres parallel to the strike of the faults, iii) the development of a number of cross-folds, striking NNE-SSW and iv) the apparently recent development of a strike-slip component on at least one of the major thrust faults. The origin of the W-E, or WSW-ENE striking structures may be reactivation of Late Cretaceous faults, stratigraphic evidence for the existence of a "structural high" (the Hurunui High) over the majority of the area in the Late Cretaceous to Early Eocene times suggests the formation of a W-E trending horst structure, with a corresponding asymmetric graben to the south. The junction of WSW-ENE trending structures with N-S trending structures to the west centres on an alluvial-filled depression, Waikari Flat, into which the structures of the WSW-ENE trending imbricate thrust system plunge, locally curling to the SW at their ends to link with N-S trending structures to the south. Roof thrusting on two orientations, W-E and N-S, towards to SE is currently occurring above these structures. Currently the area is not highly seismically active, although a magnitude ~6.4 Ms earthquake in historic times has been recorded. The effects of tectonics on the drainage of the area does suggest that the majority of the systems, are still potentially active, albeit moving at a comparatively slow rate. The majority of the recent motion appears to be concentrated on the roof-thrusting occurring in Waikari Flat, and uplift along the Lowry Peaks Fault System. Increasing amounts of secondary movement on back-thrusts and cross fractures is also implied for western ends of the major imbricate thrust system. In contrast, the southern-most fault system appears to be largely sustaining dextral strike-slip motion, with some local folding in central portions.
Bulk rock strength is greatly dependent on fracture density, so that reductions in rock strength associated with faulting and fracturing should be reflected by reduced shear coupling and hence S-wave velocity. This study is carried out along the Canterbury rangefront and in Otago. Both lie within the broader plate boundary deformation zone in the South Island of New Zealand. Therefore built structures are often, , located in areas where there are undetected or poorly defined faults with associated rock strength reduction. Where structures are sited near to, or across, such faults or fault-zones, they may sustain both shaking and ground deformation damage during an earthquake. Within this zone, management of seismic hazards needs to be based on accurate identification of the potential fault damage zone including the likely width of off-plane deformation. Lateral S-wave velocity variability provides one method of imaging and locating damage zones and off-plane deformation. This research demonstrates the utility of Multi-Channel Analysis of Surface Waves (MASW) to aid land-use planning in such fault-prone settings. Fundamentally, MASW uses surface wave dispersive characteristics to model a near surface profile of S-wave velocity variability as a proxy for bulk rock strength. The technique can aid fault-zone planning not only by locating and defining the extent of fault-zones, but also by defining within-zone variability that is readily correlated with measurable rock properties applicable to both foundation design and the distribution of surface deformation. The calibration sites presented here have well defined field relationships and known fault-zone exposure close to potential MASW survey sites. They were selected to represent a range of progressively softer lithologies from intact and fractured Torlesse Group basement hard rock (Dalethorpe) through softer Tertiary cover sediments (Boby’s Creek) and Quaternary gravels. This facilitated initial calibration of fracture intensity at a high-velocity-contrast site followed by exploration of the limits of shear zone resolution at lower velocity contrasts. Site models were constructed in AutoCAD in order to demonstrate spatial correlations between S-wave velocity and fault zone features. Site geology was incorporated in the models, along with geomorphology, river profiles, scanline locations and crosshole velocity measurement locations. Spatial data were recorded using a total-station survey. The interpreted MASW survey results are presented as two dimensional snapshot cross-sections of the three dimensional calibration-site models. These show strong correlations between MASW survey velocities and site geology, geomorphology, fluvial profiles and geotechnical parameters and observations. Correlations are particularly pronounced where high velocity contrasts exist, whilst weaker correlations are demonstrated in softer lithologies. Geomorphic correlations suggest that off-plane deformation can be imaged and interpreted in the presence of suitable topographic survey data. A promising new approach to in situ and laboratory soft-rock material and mass characterisation is also presented using a Ramset nail gun. Geotechnical investigations typically involve outcrop and laboratory scale determination of rock mass and material properties such as fracture density and unconfined compressive strength (UCS). This multi-scale approach is espoused by this study, with geotechnical and S-wave velocity data presented at multiple scales, from survey scale sonic velocity measurements, through outcrop scale scanline and crosshole sonic velocity measurements to laboratory scale property determination and sonic velocity measurements. S-wave velocities invariably increased with decreasing scale. These scaling relationships and strategies for dealing with them are investigated and presented. Finally, the MASW technique is applied to a concealed fault on the Taieri Ridge in Macraes Flat, Central Otago. Here, high velocity Otago Schist is faulted against low velocity sheared Tertiary and Quaternary sediments. This site highlights the structural sensitivity of the technique by apparently constraining the location of the principal fault, which had been ambiguous after standard processing of the seismic reflection data. Processing of the Taieri Ridge dataset has further led to the proposal of a novel surface wave imaging technique termed Swept Frequency Imaging (SFI). This inchoate technique apparently images the detailed structure of the fault-zone, and is in agreement with the conventionally-determined fault location and an existing partial trench. Overall, the results are promising and are expected to be supported by further trenching in the near future.
