Predictive modelling provides an efficient means to analyse the coastal environment and generate knowledge for long term urban planning. In this study, the numerical models SWAN and XBeach were incorporated into the ESRI ArcGIS interface by means of the BeachMMtool. This was applied to the Greater Christchurch coastal environment to simulate geomorphological evolution through hydrodynamic forcing. Simulations were performed using the recent sea level rise predictions by the Intergovernmental Panel on Climate Change (2013) to determine whether the statutory requirements outlined in the New Zealand Coastal Policy Statement 2010 are consistent with central, regional and district designations. Our results indicate that current land use zoning in Greater Christchurch is not consistent with these predictions. This is because coastal hazard risk has not been thoroughly quantified during the process of installing the Canterbury Earthquake Recovery Authority residential red zone. However, the Christchurch City Council’s flood management area does provide an extent to which managed coastal retreat is a real option. The results of this research suggest that progradation will continue to occur along the Christchurch foreshore due to the net sediment flux retaining an onshore direction and the current hydrodynamic activity not being strong enough to move sediment offshore. However, inundation during periods of storm surge poses a risk to human habitation on low lying areas around the Avon-Heathcote Estuary and the Brooklands lagoon.
Seismically vulnerable buildings constitute a major problem for the safety of human beings. In many parts of the world, reinforced concrete (RC) frame buildings designed and constructed with substandard detailing, no consideration of capacity design principles, and improper or no inclusion of the seismic actions, have been identified. Amongst those vulnerable building, one particular typology representative of the construction practice of the years previous to the 1970’s, that most likely represents the worst case scenario, has been widely investigated in the past. The deficiencies of that building typology are related to non-ductile detailing in beam column joints such as the use of plain round bars, the lack of stirrups inside the joint around the longitudinal reinforcement of the column, the use of 180° end hooks in the beams, the use of lap splices in potential ‘plastic hinge’ regions, and substandard quality of the materials. That type of detailing and the lack of a capacity design philosophy create a very fragile fuse in the structure where brittle inelastic behaviour is expected to occur, which is the panel zone region of exterior beam column joints. The non-ductile typology described above was extensively investigated at the University of Canterbury in the context of the project ‘Retrofit Solutions for New Zealand Multi-Storey Buildings’ (2004-2011), founded by the ‘Foundation for Research, Science and Technology’ Tūāpapa Rangahau Pūtaiao. The experimental campaign prior to the research carried out by the author consisted of quasi-static tests of beam column joint subassemblies subjected to lateral loading regime, with constant and varying axial load in the column. Most of those specimens were representative of a plane 2D frame (knee joint), while others represented a portion of a space 3D frame (corner joints), and only few of them had a floor slab, transverse beams, and lap splices. Using those experiments, several feasible, cost-effective, and non-invasive retrofit techniques were developed, improved, and refined. Nevertheless, the slow motion nature of those experiments did not take into account the dynamical component inherent to earthquake related problems. Amongst the set of techniques investigated, the use of FRP layers for strengthening beam column joints is of particular interest due to its versatility and the momentum that its use has gained in the current state of the practice. That particular retrofit technique was previously used to develop a strengthening scheme suitable for plane 2D and space 3D corner beam column joints, but lacking of floor slabs. In addition, a similar scheme was not developed for exterior joints of internal frames, referred here as ‘cruciform’. In this research a 2/5 scale RC frame model building comprising of two frames in parallel (external and internal) joined together by means of floor slabs and transverse beams, with non-ductile characteristics identical to those of the specimens investigated previously by others, and also including lap splices, was developed. In order to investigate the dynamic response of that building, a series of shake table tests with different ground motions were performed. After the first series of tests, the specimen was modified by connecting the spliced reinforcement in the columns in order to capture a different failure mode. Ground motions recorded during seismic events that occurred during the initial period of the experimental campaign (2010) were used in the subsequent experiments. The hierarchy of strengths and sequence of events in the panel zone region were evaluated in an extended version of the bending moment-axial load (M-N) performance domain developed by others. That extension was required due to the asymmetry in the beam cross section introduced by the floor slab. In addition, the effect of the torsion resistance provided by the spandrel (transverse beam) was included. In order to upgrade the brittle and unstable performance of the as-built/repaired specimen, a practical and suitable ad-hoc FRP retrofit intervention was developed, following a partial retrofit strategy that aimed to strengthen exterior beam column joints only (corner and cruciform). The ability of the new FRP scheme to revert the sequence of events in the panel zone region was evaluated using the extended version of the M-N performance domain as well as the guidelines for strengthening plane joints developed by others. Weakening of the floor slab in a novel configuration was also incorporated with the purpose of reducing the flexural capacity of the beam under negative bending moment (slab in tension), enabling the damage relocation from the joint into the beam. The efficacy of the developed retrofit intervention in upgrading the seismic performance of the as-built specimen was investigated using shake table tests with the input motions used in the experiments of the as-built/repaired specimen. Numerical work aimed to predict the response of the model building during the most relevant shake table tests was carried out. By using a simple numerical model with concentrated plasticity elements constructed in Ruaumoko2D, the results of blind and post-experimental predictions of the response of the specimen were addressed. Differences in the predicted response of the building using the nominal and the actual recorded motions of the shake table were investigated. The dependence of the accuracy of the numerical predictions on the assumed values of the parameters that control the hysteresis rules of key structural members was reviewed. During the execution of the experimental campaign part of this thesis, two major earthquakes affected the central part of Chile (27 of February 2010 Maule earthquake) and the Canterbury region in New Zealand (22 February 2011 Canterbury earthquake), respectively. As the author had the opportunity to experience those events and investigate their consequences in structures, the observations related to non-ductile detailing and drawbacks in the state of the practice related to reinforced concrete walls was also addressed in this research, resulting in preliminary recommendations for the refinement of current seismic code provisions and assessment guidelines. The investigations of the ground motions recorded during those and other earthquakes were used to review the procedures related to the input motions used for nonlinear dynamic analysis of buildings as required by most of the current code provisions. Inelastic displacement spectra were constructed using ground motions recorded during the earthquakes mentioned above, in order to investigate the adequacy of modification factors used to obtain reduced design spectra from elastic counterparts. Finally a simplified assessment procedure for RC walls that incorporates capacity compatible spectral demands is proposed.
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.
The Mw 7.1 Darfield earthquake generated a ~30 km long surface rupture on the Greendale Fault and significant surface deformation related to related blind faults on a previously unrecognized fault system beneath the Canterbury Plains. This earthquake provided the opportunity for research into the patterns and mechanisms of co-seismic and post-seismic crustal deformation. In this thesis I use multiple across-fault EDM surveys, logic trees, surface investigations and deformation feature mapping, seismic reflection surveying, and survey mark (cadastral) re-occupation using GPS to quantify surface displacements at a variety of temporal and spatial scales. My field mapping investigations identified shaking and crustal displacement-induced surface deformation features south and southwest of Christchurch and in the vicinity of the projected surface traces of the Hororata Blind and Charing Cross Faults. The data are consistent with the high peak ground accelerations and broad surface warping due to underlying reverse faulting on the Hororata Blind Fault and Charing Cross Fault. I measured varying amounts of post-seismic displacement at four of five locations that crossed the Greendale Fault. None of the data showed evidence for localized dextral creep on the Greendale Fault surface trace, consistent with other studies showing only minimal regional post-seismic deformation. Instead, the post-seismic deformation field suggests an apparent westward translation of northern parts of the across-fault surveys relative to the southern parts of the surveys that I attribute to post-mainshock creep on blind thrusts and/or other unidentified structures. The seismic surveys identified a deformation zone in the gravels that we attribute to the Hororata Blind Fault but the Charing Cross fault was not able to be identified on the survey. Cadastral re-surveys indicate a deformation field consistent with previously published geodetic data. We use this deformation with regional strain rates to estimate earthquake recurrence intervals of ~7000 to > 14,000 yrs on the Hororata Blind and Charing Cross Faults.
This thesis considers the presence and potential readings of graffiti and street art as part of the wider creative public landscape of Christchurch in the wake of the series of earthquakes that significantly disrupted the city physically and socially. While documenting a specific and unprecedented period of time in the city’s history, the prominence of graffiti and street art throughout the constantly changing landscape has also highlighted their popularity as increasingly entrenched additions to urban and suburban settings across the globe. In post-quake Christchurch, graffiti and street art have often displayed established tactics, techniques and styles while exploring and exposing the unique issues confronting this disrupted environment, illustrating both a transposable nature and the entwined relationship with the surrounding landscape evident in the conception of these art forms. The post-quake city has afforded graffiti and street art the opportunity to engage with a range of concepts: from the re-activation and re-population of the empty and abandoned spaces of the city, to commentaries on specific social and political issues, both angry and humorous, and notably the reconsideration of entrenched and evolving traditions, including the distinction between guerrilla and sanctioned work. The examples of graffiti and street art within this work range from the more immediate post-quake appearance of art in a group of affected suburbs, including the increasingly empty residential red-zone, to the use of the undefined spaces sweeping the central city, and even inside the Canterbury Museum, which housed the significant street art exhibition Rise in 2013-2014. These settings expose a number of themes, both distinctive and shared, that relate to both the post-disaster landscape and the concerns of graffiti and street art as art movements unavoidably entangled with public space.
The Amuri Earthquake of September 1, 1888 (magnitude M = 6.5 to 6.8) occurred on the Hope River Segment of the Hope Fault west of Hanmer Plains. The earthquake was felt strongly in North Canterbury and North Westland and caused considerable property damage and landsliding in the Lower Hope Valley. However, damage reports and the spatial distribution of felt intensities emphasize extreme variations in seismic effects over short distances, probably due to topographic focusing and local ground conditions. Significant variations in lateral fault displacement occurred at secondary fault segment boundaries (side-steps and bends in the fault trace) during the 1888 earthquake. This historical spatial variation in lateral slip is matched by the Late Quaternary geomorphic distribution of slip on the Hope River Segment of the Hope Fault. Trenching studies at two sites on the Hope Fault have also identified evidence for five pre-historic earthquakes of similar magnitude to the 1888 earthquake and an average recurrence interval of 134 ± 27 years between events. Magnitude estimates for the 1888 earthquake are combined with a. strong ground motion attenuation expression to provide an estimate of potential ground accelerations in Amuri District during-future earthquakes on the Hope River Segment of the Hope Fault. The predicted acceleration response on bedrock sites within 20 km of the epicentral region is between 0.23 g and 0.34 g. The close match between the historic, inferred pre-historic and geomorphic distribution of lateral slip indicates that secondary fault segmentation exerts a strong structural control on rupture propagation and the expression of fault displacement at the surface. In basement rocks at depth the spatial variations in slip are inferred to be distributed within zones of pervasive cataclastic shear, on either side of the fault segment boundaries. The large variations in surface displacement across fault segment boundaries means that one must know the geometry of the fault in order to evaluate slip-rates calculated from individual locations. The average Late Quaternary slip-rate on the Hope Fault at Glynn Wye Station is between 15.5 mm/yr and 18.25 mm/yr and the rate on the subsidiary Kakapo Fault is between 5.0 mm/yr and 7.5 mm/yr. These rates have been determined from sites which are relatively free of structural complication.
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.
Post-tensioned timber technology was originally developed and researched at the University of Canterbury (UC) in New Zealand in 2005. It can provide a low-damage seismic design solution for multi-storey mass timber buildings. Since mass timber products, such as cross-laminated timber (CLT), have high in-plane stiffness, a post-tensioned timber shear wall will deform mainly in a rocking mechanism. The moment capacity of the wall at the base is commonly determined using the elastic form of the Modified Monolithic Beam Analogy (MMBA). In the calculation of the moment capacity at the wall base, it is critical to accurately predict the location of the neutral axis and the timber compressive stress distribution. Three 2/3 scale 8.6m tall post-tensioned CLT walls were experimentally tested under quasi-static cyclic loading – both uni-directional and bi-directional- in this study. These specimens included a single wall, a coupled wall, and a C-shaped core-wall. The main objective was to develop post-tensioned C-shaped timber core-walls for tall timber buildings with enhanced lateral strength and stiffness. To better understand the timber compressive stress distributions at the wall base, particle tracking technology (PTT) technology was applied for the first time to investigate the behaviour of the compression toe. Previous post-tensioned timber testing primarily used the displacement measurements to determine the timber compressive behavior at the wall base or rocking interfaces. However, by using PTT technology, the timber strain measurements in the compression zone can be much more accurate as PTT is able to track the movement of many particles on the timber surface. This paper presents experimental testing results of post-tensioned CLT walls with a focus on capturing timber compressive behavior using PTT. The PTT measurements were able to better capture small base rotations which occurred at the onset of gap opening and capture unexpected phenomena in core-wall tests. The single wall test result herein presented indicates that while the MMBA could predict the moment rotation behavior with reasonable accuracy, the peak strain response was under predicted in the compression toe. Further detailed study is required to better understand the complex strain fields generated reflective of the inherent cross-thickness inhomogeneity and material variability of CLT.
Liquefaction features and the geologic environment in which they formed were carefully studied at two sites near Lincoln in southwest Christchurch. We undertook geomorphic mapping, excavated trenches, and obtained hand cores in areas with surficial evidence for liquefaction and areas where no surficial evidence for liquefaction was present at two sites (Hardwick and Marchand). The liquefaction features identified include (1) sand blows (singular and aligned along linear fissures), (2) blisters or injections of subhorizontal dikes into the topsoil, (3) dikes related to the blows and blisters, and (4) a collapse structure. The spatial distribution of these surface liquefaction features correlates strongly with the ridges of scroll bars in meander settings. In addition, we discovered paleoliquefaction features, including several dikes and a sand blow, in excavations at the sites of modern liquefaction. The paleoliquefaction event at the Hardwick site is dated at A.D. 908-1336, and the one at the Marchand site is dated at A.D. 1017-1840 (95% confidence intervals of probability density functions obtained by Bayesian analysis). If both events are the same, given proximity of the sites, the time of the event is A.D. 1019-1337. If they are not, the one at the Marchand site could have been much younger. Taking into account a preliminary liquefaction-triggering threshold of equivalent peak ground acceleration for an Mw 7.5 event (PGA7:5) of 0:07g, existing magnitude-bounded relations for paleoliquefaction, and the timing of the paleoearthquakes and the potential PGA7:5 estimated for regional faults, we propose that the Porters Pass fault, Alpine fault, or the subduction zone faults are the most likely sources that could have triggered liquefaction at the study sites. There are other nearby regional faults that may have been the source, but there is no paleoseismic data with which to make the temporal link.
The National Science Challenge ‘Building Better Homes, Towns and Cities’ is currently undertaking work that, in part, identifies and analyses the Waimakariri District Council’s (WMK/Council) organisational practices and process tools. The focus is on determining the processes that made the Residential Red Zone Recovery Plan, 2016 (RRZRP) collaboration process so effective and compares it to the processes used to inform the current Kaiapoi Town Centre Plan - 2028 and Beyond (KTC Plan). This research aims to explore ‘what travelled’ in terms of values, principles, methods, processes and personnel from the RRZRP to the KTC planning process. My research will add depth to this research by examining more closely the KTC Plan’s hearings process, reviewing submissions made, analysing background documents and by conducting five semi-structured interviews with a selection of people who made submissions on the KTC Plan. The link between community involvement and best recovery outcomes has been acknowledged in literature as well as by humanitarian agencies (Lawther, 2009; Sullivan, 2003). My research has documented WMK’s post-quake community engagement strategy by focusing on their initial response to the earthquake of 2010 and the two-formal plan (RRZRP and KTC Plan) making procedures that succeeded this response. My research has led me to conclude that WMK was committed to collaborating with their constituents right through the extended post-quake sequence. Iterative face to face or ‘think communications’ combined with the accessibility of all levels of Council staff – including senior management and elected members - gave interested community members the opportunity to discuss and deliberate the proposed plans with the people tasked with preparing them. WMK’s commitment to collaborate is illustrated by the methods they employed to inform their post-quake efforts and plans and by the logic behind the selected methods. Combined the Council’s logic and methods best describe the ‘Waimakariri Way’. My research suggests that collaborative planning is iterative in nature. It is therefore difficult to establish a specific starting point where collaboration begins as the relationships needed for the collaborative process constantly (re)emerge out of pre-existing relationships. Collaboration seems to be based on an attitude, which means there is no starting ‘point’ as such, rather an amplification for a time of a basic attitude towards the public.
A city’s planted trees, the great majority of which are in private gardens, play a fundamental role in shaping a city’s wild ecology, ecosystem functioning, and ecosystem services. However, studying tree diversity across a city’s many thousands of separate private gardens is logistically challenging. After the disastrous 2010–2011 earthquakes in Christchurch, New Zealand, over 7,000 homes were abandoned and a botanical survey of these gardens was contracted by the Government’s Canterbury Earthquake Recovery Authority (CERA) prior to buildings being demolished. This unprecedented access to private gardens across the 443.9 hectares ‘Residential Red Zone’ area of eastern Christchurch is a unique opportunity to explore the composition of trees in private gardens across a large area of a New Zealand city. We analysed these survey data to describe the effects of housing age, socio-economics, human population density, and general soil quality, on tree abundance, species richness, and the proportion of indigenous and exotic species. We found that while most of the tree species were exotic, about half of the individual trees were local native species. There is an increasing realisation of the native tree species values among Christchurch citizens and gardens in more recent areas of housing had a higher proportion of smaller/younger native trees. However, the same sites had proportionately more exotic trees, by species and individuals, amongst their larger planted trees than older areas of housing. The majority of the species, and individuals, of the larger (≥10 cm DBH) trees planted in gardens still tend to be exotic species. In newer suburbs, gardens in wealthy areas had more native trees than gardens from poorer areas, while in older suburbs, poorer areas had more native big trees than wealthy areas. In combination, these describe, in detail unparalleled for at least in New Zealand, how the tree infrastructure of the city varies in space and time. This lays the groundwork for better understanding of how wildlife distribution and abundance, wild plant regeneration, and ecosystem services, are affected by the city’s trees.
The Acheron rock avalanche is located in the Red Hill valley almost 80 km west of Christchurch and is one of 42 greywacke-derived rock avalanches identified in the central Southern Alps. It overlies the Holocene active Porters Pass Fault; a component of the Porters Pass-Amberley Fault Zone which extends from the Rakaia River to beyond the Waimakariri River. The Porters Pass Fault is a dextral strike-slip fault system viewed as a series of discontinuous fault scarps. The location of the fault trace beneath the deposit suggests it may represent a possible source of seismic shaking resulting in the formation of the Acheron rock avalanche. The rock mass composition of the rock avalanche source scar is Torlesse Supergroup greywacke consisting of massive sandstone and thinly bedded mudstone sequences dipping steeply north into the centre of the source basin. A stability analysis identified potential instability along shallow north dipping planar defects, and steep south dipping toppling failure planes. The interaction of the defects with bedding is considered to have formed conditions for potential instability most likely triggered by a seismic event. The dTositional area of the rock avalanche covers 7.2 x 105 m2 with an estimated volume of 9 x 10 m3 The mobilised rock mass volume was calculated at 7.5 x 106 m3• Run out of the debris from the top of the source scar to the distal limit reached 3500m, descending over a vertical fall of almost 700m with an estimated Fahrboschung of 0.2. The run out of the rock avalanche displayed moderate to high mobility, travelling at an estimated maximum velocity of 140-160 km/hour. The rapid emplacement of the deposit is confirmed by highly fragmented internal composition and burial of forest vegetation New radiocarbon ages from buried wood retrieved from the base of Acheron rock avalanche deposit represents an emplacement age closely post-dating (Wk 12094) 1152 ± 51 years B.P. This differs significantly from a previous radiocarbon age of (NZ547) 500 ± 69 years B.P. and modal lichenometry and weathering-rind thickness ages of approximately 460 ± 10 yrs and 490 ± 50 years B.P. The new age shows no resemblance to an earthquake event around 700- 500 years B.P. on the Porters Pass-Amberley Fault Zone. The DAN run out simulation using a friction model rheology successfully replicated the long run out and velocity of the Acheron rock avalanche using a frictron angle of 27° and high earth pressure coefficients of 5.5, 5.2, and 5.9. The elevated earth pressure coefficients represent dispersive pressures derived from dynamic fragmentation of the debris within the mobile rock avalanche, supporting the hypothesis of Davies and McSaveney (2002). The DAN model has potential applications for areas prone to large-scale instability in the elevated slopes and steep waterways of the Southern Alps. A paleoseismic investigation of a newly identified scarp of the Porters Pass Fault partially buried by the rock avalanche was conducted to identify any evidence of a coseismic relationship to the Acheron rock avalanche. This identified three-four fault traces striking at 078°, and a sag pond displaying a sequence of overbank deposits containing two buried soils representing an earthquake event horizon. A 40cm vertical offset of the ponded sediment and lower buried soil horizqn was recorded, which was dated to (Wk 13112 charcoal in palosol) 653 ± 54 years B.P. and (Wk 13034 palosol) 661 ± 34 years B.P. The evidence indicates a fault rupture occurred along the Porters Pass Fault, west of Porters Pass most likely extending to the Red Lakes terraces, post-dating 700 years B.P., resulting in 40cm of vertical displacement and an unknown component of dextral strike slip movement. This event post dates the event one (1000 ± 100 years B.P) at Porters Pass previously considered to represent the most recent rupture along the fault line. This points to a probable source for resetting of the modal weathering-rind thicknesses and lichen size populations in the Red Hill valley and possibly the Red Lakes terraces. These results suggest careful consideration must be given to the geomorphic and paleoseismic history of a specific site when applying surface dating techniques and furthermore the origin of dates used in literature and their useful range should be verified. An event at 700-500 years B.P did not trigger the Acheron rock avalanche as previously assumed supporting Howard's conclusions. The lack of similar aged rupture evidence in either of the Porters Pass and Coleridge trenches supports Howard's hypothesis of segmentation of the Porters Pass Fault; where rupture occurs along one fault segment but not along another. The new rock avalanche age closely post-dating 1200-1100 years B.P. resembles the poorly constrained event one rupture age of 1700-800 years B.P for the Porters Pass Fault and the tighter constrained Round Top event of 1010 ± 50 years B.P. on the Alpine Fault. Eight other rock avalanche deposits spread across the central Southern Alps also resemble the new ages however are unable to be assigned specific earthquake events due to the large associated error bars of± 270 years. This clustering of ages does represent compelling lines of evidence for large magnitude earthquake events occurring over the central Southern Alps. The presence of a rock avalanche deposit does not signify an earthquake based on the historical evidence in the Southern Alps however clustering of ages does suggest that large Mw >7 earthquakes occurred across the Southern Alps between 1200-900 years BP.
The Canterbury Region is susceptible to a variety of natural hazards, including earthquakes, landslides and climate hazards. Increasing population and tourism within the region is driving development pressures and as more and more development occurs, the risk from natural hazards increases. In order to avoid development occurring in unacceptably vulnerable locations, natural hazard assessments are required. This study is a reconnaissance natural hazard assessment of Lakes Lyndon, Coleridge and Tekapo. There is restricted potential for development at Lake Lyndon, because the land surrounding the lake is owned by the Crown and has a number of development restrictions. However, there is the potential for conservation or recreation-linked development to occur. There is more potential for development at Lake Coleridge. Most of the land surrounding the lake is privately owned and has less development restrictions. The majority of land surrounding Lake Tekapo is divided into Crown-owned pastoral leases, which are protected from development, such as subdivision. However, there are substantial areas around the lake, which are privately owned and, therefore, have potential for development. Earthquake, landslide and climate hazards are the main natural hazards threatening Lakes Lyndon, Coleridge and Tekapo. The lakes are situated in a zone of active earth deformation in which large and relatively frequent earthquakes are produced. A large number of active faults lie within 15 km of each lake, which are capable of producing M7 or larger earthquakes. Ground shaking, liquefaction, landslides, tsunami and seiches are among the consequences of earthquakes, all of which have the potential to cause severe damage to lives, lifelines and infrastructure. Landslides are also common in the landscape surrounding the lakes. The majority of slopes surrounding the lakes are at significant risk from earthquake-induced failure under moderate to strong earthquake shaking. This level of shaking is expected to occur in any 50 year period around Lakes Lyndon and Coleridge, and in any 150 year period around Lake Tekapo. Injuries, fatalities and property damage can occur directly from landslide impact or from indirect effects such as flooding from landslide-generated tsunami or from landslide dam outbreaks. Lakes Lyndon, Coleridge and Tekapo are also susceptible to climate hazards, such as high winds, drought, heavy snowfall and heavy rainfall, which can lead to landslides and flooding. Future climate change due to global warming is most likely going to affect patterns of frequency and magnitudes of extreme weather events, leading to an increase in climate hazards. Before development is permitted around the lakes, it is essential that each of these hazards is considered so that unacceptably vulnerable areas can be avoided.
Mitigating the cascade of environmental damage caused by the movement of excess reactive nitrogen (N) from land to sea is currently limited by difficulties in precisely and accurately measuring N fluxes due to variable rates of attenuation (denitrification) during transport. This thesis develops the use of the natural abundance isotopic composition of nitrate (δ15N and δ18O of NO₃-) to integrate the spatialtemporal variability inherent to denitrification, creating an empirical framework for evaluating attenuation during land to water NO₃- transfers. This technique is based on the knowledge that denitrifiers kinetically discriminate against 'heavy' forms of both N and oxygen (O), creating a parallel enrichment in isotopes of both species as the reaction progresses. This discrimination can be quantitatively related to NO₃- attenuation by isotopic enrichment factors (εdenit). However, while these principles are understood, use of NO₃- isotopes to quantify denitrification fluxes in non-marine environments has been limited by, 1) poor understanding of εdenit variability, and, 2) difficulty in distinguishing the extent of mixing of isotopically distinct sources from the imprint of denitrification. Through a combination of critical literature analysis, mathematical modelling, mesocosm to field scale experiments, and empirical studies on two river systems over distance and time, these short comings are parametrised and a template for future NO₃- isotope based attenuation measurements outlined. Published εdenit values (n = 169) are collated in the literature analysis presented in Chapter 2. By evaluating these values in the context of known controllers on the denitrification process, it is found that the magnitude of εdenit, for both δ15N and δ18O, is controlled by, 1) biology, 2) mode of transport through the denitrifying zone (diffusion v. advection), and, 3) nitrification (spatial-temporal distance between nitrification and denitrification). Based on the outcomes of this synthesis, the impact of the three factors identified as controlling εdenit are quantified in the context of freshwater systems by combining simple mathematical modelling and lab incubation studies (comparison of natural variation in biological versus physical expression). Biologically-defined εdenit, measured in sediments collected from four sites along a temperate stream and from three tropical submerged paddy fields, varied from -3‰ to -28‰ depending on the site’s antecedent carbon content. Following diffusive transport to aerobic surface water, εdenit was found to become more homogeneous, but also lower, with the strength of the effect controlled primarily by diffusive distance and the rate of denitrification in the sediments. I conclude that, given the variability in fractionation dynamics at all levels, applying a range of εdenit from -2‰ to -10‰ provides more accurate measurements of attenuation than attempting to establish a site-specific value. Applying this understanding of denitrification's fractionation dynamics, four field studies were conducted to measure denitrification/ NO₃- attenuation across diverse terrestrial → freshwater systems. The development of NO₃- isotopic signatures (i.e., the impact of nitrification, biological N fixation, and ammonia volatilisation on the isotopic 'imprint' of denitrification) were evaluated within two key agricultural regions: New Zealand grazed pastures (Chapter 4) and Philippine lowland submerged rice production (Chapter 5). By measuring the isotopic composition of soil ammonium, NO₃- and volatilised ammonia following the bovine urine deposition, it was determined that the isotopic composition of NO₃ - leached from grazed pastures is defined by the balance between nitrification and denitrification, not ammonia volatilisation. Consequently, NO₃- created within pasture systems was predicted to range from +10‰ (δ15N)and -0.9‰ (δ18O) for non-fertilised fields (N limited) to -3‰ (δ15N) and +2‰ (δ18O) for grazed fertilised fields (N saturated). Denitrification was also the dominant determinant of NO₃- signatures in the Philippine rice paddy. Using a site-specific εdenit for the paddy, N inputs versus attenuation were able to be calculated, revealing that >50% of available N in the top 10 cm of soil was denitrified during land preparation, and >80% of available N by two weeks post-transplanting. Intriguingly, this denitrification was driven by rapid NO₃- production via nitrification of newly mineralised N during land preparation activities. Building on the relevant range of εdenit established in Chapters 2 and 3, as well as the soil-zone confirmation that denitrification was the primary determinant of NO₃- isotopic composition, two long-term longitudinal river studies were conducted to assess attenuation during transport. In Chapter 6, impact and recovery dynamics in an urban stream were assessed over six months along a longitudinal impact gradient using measurements of NO₃- dual isotopes, biological populations, and stream chemistry. Within 10 days of the catastrophic Christchurch earthquake, dissolved oxygen in the lowest reaches was <1 mg l⁻¹, in-stream denitrification accelerated (attenuating 40-80% of sewage N), microbial biofilm communities changed, and several benthic invertebrate taxa disappeared. To test the strength of this method for tackling the diffuse, chronic N loading of streams in agricultural regions, two years of longitudinal measurements of NO₃- isotopes were collected. Attenuation was negatively correlated with NO₃- concentration, and was highly dependent on rainfall: 93% of calculated attenuation (20 kg NO₃--N ha⁻¹ y⁻¹) occurred within 48 h of rainfall. The results of these studies demonstrate the power of intense measurements of NO₃- stable isotope for distinguishing temporal and spatial trends in NO₃ - loss pathways, and potentially allow for improved catchment-scale management of agricultural intensification. Overall this work now provides a more cohesive understanding for expanding the use of NO₃- isotopes measurements to generate accurate understandings of the controls on N losses. This information is becoming increasingly important to predict ecosystem response to future changes, such the increasing agricultural intensity needed to meet global food demand, which is occurring synergistically with unpredictable global climate change.
As the future of the world’s oil reserves becomes progressively more uncertain, it is becoming increasingly important that steps are taken to ensure that there are viable, attractive alternatives to travel by private motor vehicle. As with many of New Zealand’s major urban centres, Christchurch is still exceptionally reliant on private motor vehicles; although a significant proportion of the population indicate that they would like to cycle more, cycling is still an underutilised mode of transport. Following a series of fatal earthquakes that struck the city in 2010 and 2011, there has been the need to significantly redevelop much of the city’s horizontal infrastructure – subsequently providing the perfect platform for significant changes to be made to the road network. Many of the key planning frameworks governing the rebuild process have identified the need to improve Christchurch’s cycling facilities in order to boost cycling numbers and cyclist safety. The importance of considering future growth and travel patterns when planning for transport infrastructure has been highlighted extensively throughout literature. Accordingly, this study sought to identify areas where future cycle infrastructure development would be advantageous based on a number of population and employment projections, and likely future travel patterns throughout the city. Through the use of extensive GIS analysis, future population growth, employment and travel patterns for Christchurch city were examined in order to attain an understanding of where the current proposed major cycleways network could be improved, or extended. A range of data and network analysis were used to derive likely travel patterns throughout Christchurch in 2041. Trips were derived twice, once with a focus on simply finding the shortest route between each origin and destination, and then again with a focus on cyclist safety and areas where cyclists were unlikely to travel. It was found that although the proposed major cycleways network represents a significant step towards improving the cycling environment in Christchurch, there are areas of the city that will not be well serviced by the current proposed network in 2041. These include a number of key residential growth areas such as Halswell, Belfast and Prestons, along with a number of noteworthy key travel zones, particularly in areas close to the central city and key employment areas. Using network analysis, areas where improvements or extensions to the proposed network would be most beneficial were identified, and a number of potential extensions in a variety of areas throughout the city were added to the network of cycle ways. Although it has been found that filling small gaps in the network can have considerable positive outcomes, results from the prioritisation analysis suggested that initially in Christchurch demand is likely to be for more substantial extensions to the proposed major cycleways network.
Millions of urban residents around the world in the coming century will experience severe landscape change – including increased frequencies of flooding due to intensifying storm events and impacts from sea level rise. For cities, collisions of environmental change with mismatched cultural systems present a major threat to infrastructure systems that support urban living. Landscape architects who address these issues express a need to realign infrastructure with underlying natural systems, criticizing the lack of social and environmental considerations in engineering works. Our ability to manage both society and the landscapes we live in to better adapt to unpredictable events and landscape changes is essential if we are to sustain the health and safety of our families, neighbourhoods, and wider community networks. When extreme events like earthquakes or flooding occur in developed areas, the feasibility of returning the land to pre-disturbance use can be questioned. In Christchurch for example, a large expanse of land (630 hectares) within the city was severely damaged by the earthquakes and judged too impractical to repair in the short term. The central government now owns the land and is currently in the process of demolishing the mostly residential houses that formed the predominant land use. Furthermore, cascading impacts from the earthquakes have resulted in a general land subsidence of .5m over much of eastern Christchurch, causing disruptive and damaging flooding. Yet, although disasters can cause severe social and environmental distress, they also hold great potential as a catalyst to increasing adaption. But how might landscape architecture be better positioned to respond to the potential for transformation after disaster? This research asks two core questions: what roles can the discipline of landscape architecture play in improving the resilience of communities so they become more able to adapt to change? And what imaginative concepts could be designed for alternative forms of residential development that better empower residents to understand and adapt the infrastructure that supports them? Through design-directed inquiry, the research found landscape architecture theory to be well positioned to contribute to goals of social-ecological systems resilience. The discipline of landscape architecture could become influential in resilience-oriented multi disciplinary collaborations, with our particular strengths lying in six key areas: the integration of ecological and social processes, improving social capital, engaging with temporality, design-led innovation potential, increasing diversity and our ability to work across multiple scales. Furthermore, several innovative ideas were developed, through a site-based design exploration located within the residential red zone, that attempt to challenge conventional modes of urban living – concepts such as time-based land use, understanding roads as urban waterways, and landscape design and management strategies that increase community participation and awareness of the temporality in landscapes.
A buckling-restrained braced frame (BRBF) is a structural bracing system that provides lateral strength and stiffness to buildings and bridges. They were first developed in Japan in the 1970s (Watanabe et al. 1973, Kimura et al. 1976) and gained rapid acceptance in the United States after the Northridge earthquake in 1994 (Bruneau et al. 2011). However, it was not until the Canterbury earthquakes of 2010/2011, that the New Zealand construction market saw a significant uptake in the use of buckling-restrained braces (BRBs) in commercial buildings (MacRae et al. 2015). In New Zealand there is not yet any documented guidance or specific instructions in regulatory standards for the design of BRBFs. This makes it difficult for engineers to anticipate all the possible stability and strength issues within a BRBF system and actively mitigate them in each design. To help ensure BRBF designs perform as intended, a peer review with physical testing are needed to gain building compliance in New Zealand. Physical testing should check the manufacturing and design of each BRB (prequalification testing), and the global strength and stability of each BRB its frame (subassemblage testing). However, the financial pressures inherent in commercial projects has led to prequalification testing (BRB only testing) being favoured without adequate design specific subassemblage testing. This means peer reviewers have to rely on BRB suppliers for assurances. This low regulation environment allows for a variety of BRBF designs to be constructed without being tested or well understood. The concern is that there may be designs that pose risk and that issues are being overlooked in design and review. To improve the safety and design of BRBFs in New Zealand, this dissertation studies the behaviour of BRBs and how they interact with other frame components. Presented is the experimental test process and results of five commercially available BRB designs (Chapter 2). It discusses the manufacturing process, testing conditions and limitations of observable information. It also emphasises that even though subassemblage testing is impractical, uniaxial testing of the BRB only is not enough, as this does not check global strength or stability. As an alternative to physical testing, this research uses computer simulation to model BRB behaviour. To overcome the traditional challenges of detailed BRB modelling, a strategy to simulate the performance of generic BRB designs was developed (Chapter 3). The development of nonlinear material and contact models are important aspects of this strategy. The Chaboche method is employed using a minimum of six backstress curves to characterize the combined isotropic and kinematic hardening exhibited by the steel core. A simplified approach, adequate for modelling the contact interaction between the restrainer and the core was found. Models also capture important frictional dissipation as well as lateral motion and bending associated with high order constrained buckling of the core. The experimental data from Chapter 2 was used to validate this strategy. As BRBs resist high compressive loading, global stability of the BRB and gusseted connection zone need to be considered. A separate study was conducted that investigated the yielding and buckling strength of gusset plates (Chapter 4). The stress distribution through a gusset plate is complex and difficult to predict because the cross-sectional area of gusset plate is not uniform, and each gusset plate design is unique in shape and size. This has motivated design methods that approximate yielding of gusset plates. Finite element modelling was used to study the development of yielding, buckling and plastic collapse behaviour of a brace end bolted to a series of corner gusset plates. In total 184 variations of gusset plate geometries were modelled in Abaqus®. The FEA modelling applied monotonic uniaxial load with an imperfection. Upon comparing results to current gusset plate design methods, it was found that the Whitmore width method for calculating the yield load of a gusset is generally un-conservative. To improve accuracy and safety in the design of gusset plates, modifications to current design methods for calculating the yield area and compressive strength for gusset plates is proposed. Bolted connections are a popular and common connection type used in BRBF design. Global out-of-plane stability tends to govern the design for this connection type with numerous studies highlighting the risk of instability initiated by inelasticity in the gussets, neck of the BRB end and/or restrainer ends. Subassemblage testing is the traditional method for evaluating global stability. However, physical testing of every BRBF variation is cost prohibitive. As such, Japan has developed an analytical approach to evaluate out-of-plane stability of BRBFs and incorporated this in their design codes. This analytical approach evaluates the different BRB components under possible collapse mechanisms by focusing on moment transfer between the restrainer and end of the BRB. The approach have led to strict criteria for BRBF design in Japan. Structural building design codes in New Zealand, Europe and the United States do not yet provide analytical methods to assess BRB and connection stability, with prototype/subassemblage testing still required as the primary means of accreditation. Therefore it is of interest to investigate the capability of this method to evaluate stability of BRBs designs and gusset plate designs used in New Zealand (including unstiffened gusset connection zones). Chapter 5 demonstrates the capability of FEA to study to the performance of a subassemblage test under cyclic loading – resembling that of a diagonal ground storey BRBF with bolted connections. A series of detailed models were developed using the strategy presented in Chapter 3. The geometric features of BRB 6.5a (Chapter 2) were used as a basis for the BRBs modelled. To capture the different failure mechanisms identified in Takeuchi et al. (2017), models varied the length that the cruciform (non-yielding) section inserts into the restrainer. Results indicate that gusset plates designed according to New Zealand’s Steel Structures Standard (NZS 3404) limit BRBF performance. Increasing the thickness of the gusset plates according to modifications discussed in Chapter 4, improved the overall performance for all variants (except when Lin/ Bcruc = 0.5). The effect of bi-directional loading was not found to notably affect out-of-plane stability. Results were compared against predictions made by the analytical method used in Japan (Takeuchi method). This method was found to be generally conservative is predicting out-of-plane stability of each BRBF model. Recommendations to improve the accuracy of Takeuchi’s method are also provided. The outcomes from this thesis should be helpful for BRB manufacturers, researchers, and in the development of further design guidance of BRBFs.
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.
Recent global tsunami events have highlighted the importance of effective tsunami risk management strategies (including land-use planning, structural and natural defences, warning systems, education and evacuation measures). However, the rarity of tsunami means that empirical data concerning reactions to tsunami warnings and tsunami evacuation behaviour is rare when compared to findings about evacuations to avoid other sources of hazard. To date empirical research into tsunami evacuations has focused on evacuation rates, rather than other aspects of the evacuation process. More knowledge is required about responses to warnings, pre-evacuation actions, evacuation dynamics and the return home after evacuations. Tsunami evacuation modelling has the potential to inform evidence-based tsunami risk planning and response. However to date tsunami evacuation models have largely focused on timings of evacuations, rather than evacuation behaviours. This Masters research uses a New Zealand case study to reduce both of these knowledge gaps. Qualitative survey data was gathered from populations across coastal communities in Banks Peninsula and Christchurch, New Zealand, required to evacuate due to the tsunami generated by the November 14th 2016 Kaikōura Earthquake. Survey questions asked about reactions to tsunami warnings, actions taken prior to evacuating and movements during the 2016 tsunami evacuation. This data was analysed to characterise trends and identify factors that influenced evacuation actions and behaviour. Finally, it was used to develop an evacuation model for Banks Peninsula. Where appropriate, the modelling inputs were informed by the survey data. Three key findings were identified from the results of the evacuation behaviour survey. Although 38% of the total survey respondents identified the earthquake shaking as a natural cue for the tsunami, most relied on receiving official warnings, including sirens, to prompt evacuations. Respondents sought further official information to inform their evacuation decisions, with 39% of respondents delaying their evacuation in order to do so. Finally, 96% of total respondents evacuated by car. This led to congestion, particularly in more densely populated Christchurch city suburbs. Prior to this research, evacuation modelling had not been completed for Banks Peninsula. The results of the modelling showed that if evacuees know how to respond to tsunami warnings and where and how to evacuate, there are no issues. However, if there are poor conditions, including if people do not evacuate immediately, if there are issues with the roading network, or if people do not know where or how to evacuate, evacuation times increase with there being more bottlenecks leading out of the evacuation zones. The results of this thesis highlight the importance of effective tsunami education and evacuation planning. Reducing exposure to tsunami risk through prompt evacuation relies on knowledge of how to interpret tsunami warnings, and when, where and how to evacuate. Recommendations from this research outline the need for public education and engagement, and the incorporation of evacuation signage, information boards and evacuation drills. Overall these findings provide more comprehensive picture of tsunami evacuation behaviour and decision making based on empirical data from a recent evacuation, which can be used to improve tsunami risk management strategies. This empirical data can also be used to inform evacuation modelling to improve the accuracy and realism of the evacuation models.
Current research in geotechnical engineering at the University of Canterbury includes a number of laboratory testing programmes focussed on understanding the behaviour of natural soil deposits in Christchurch during the 2010-2011 Canterbury Earthquake Sequence. Many soils found in Christchurch are sands or silty sands with little to no plasticity, making them very difficult to sample using established methods. The gel-push sampling methodology, developed by Kiso-Jiban Consultants in Japan, was developed to address some of the deficiencies of existing sampling techniques and has been deployed on two projects in Christchurch. Gel push sampling is carried out with a range of samplers which are modified versions of existing technology, and the University of Canterbury has acquired three versions of the tools (GP-S, GP-Tr, GP-D). Soil samples are extracted from the bottom of a freshly drilled borehole and are captured within a liner barrel, close to 1m in length. A lubricating polymer gel coats the outside of the soil sample as it enters the liner barrel. The frictional rubbing which normally occurs on the sides of the soil samples using existing techniques is eliminated by the presence of the polymer gel. The operation of the gel-push samplers is significantly more complicated than conventional push-tube samplers, and in the initial trials a number of operational difficulties were encountered, requiring changes to the sampling procedures. Despite these issues, a number of high quality soil samples were obtained on both projects using the GP-S sampler to capture silty soil. Attempts were made to obtain clean sands using a different gel-push sampler (GP-TR) in the Red Zone. The laboratory testing of these sands indicated that they were being significantly disturbed during the sampling and/or transportation procedures. While it remains too early to draw definitive conclusions regarding the performance of the gel-push samplers, the methodology has provided some promising results. Further trialling of the tools are required to refine operating procedures understand the full range of soil conditions which can be successfully sampled using the tools. In parallel with the gel-push trials, a Dames and Moore fixed-piston sampler has been used by our research partners from Berkeley to obtain soil samples at a number of sites within Christchurch. This sampler features relatively short (50cm), thin-walled liner barrels which is advanced into the ground under the action of hydraulic pressure. By reducing the overall length of the soil being captured, the disturbance to the soil as it enters the liner barrel is significantly reduced. The Dames and Moore sampler is significantly easier to operate than the gel-push sampler, and past experience has shown it to be successful in soft, plastic materials (i.e. clays and silty clays). The cyclic resistance of one silty clay obtained using both the gel-push and Dames & Moore samplers has been found to be very similar, and ongoing research aims to establish whether similar results are obtained for different soil types, including silty materials and clean sands.
4th September 2010 a 7.1 magnitude earthquake strikes near Christchurch, New Zealand’s second largest city of approximately 370,000 people. This is followed by a 6.3 magnitude quake on 22nd February 2011 and a 6.4 on 13th June. In February 181 people died and a state of national emergency was declared from 23 February to 30th April. Urban Search and Rescue teams with 150 personnel from New Zealand and 429 from overseas worked tirelessly in addition to Army, Police and Fire services. Within the central business district 1,000 buildings (of 4,000) are expected to be demolished. An estimated 10,000 houses require demolition and over 100,000 were damaged. Meanwhile the over 7,000 aftershocks have become part of the “new normal” for us all. During this time how have libraries supported their staff? What changes have been made to services? What are the resourcing opportunities? This presentation will provide a personal view from Lincoln University, Te Whare Wanaka o Aoraki, Library Teaching and Learning. Lincoln is New Zealand's third oldest university having been founded in 1878. Publicly owned and operated it is New Zealand's specialist land-based university. Lincoln is based on the Canterbury Plains, 22 kilometres south of Christchurch. On campus there was mostly minor damage to buildings while in the Library 200,000 volumes were thrown from the shelves. I will focus on the experiences of the Disaster Team and on our experiences with hosting temporarily displaced staff and students from the Christchurch Polytechnic Institute of Technology, Library, Learning & Information Services. Experiences from two other institutions will be highlighted: Christchurch City Libraries, Ngā Kete Wānanga-o-Ōtautahi. Focusing on the Māori Services Team and the Ngā Pounamu Māori and Ngāi Tahu collections. The Central library located within the red zone cordon has been closed since February, the Central library held the Ngā Pounamu Māori and Ngai Tahu collections, the largest Māori collections in the Christchurch public library network. The lack of access to these collections changed the way the Māori Services Team, part of the larger Programmes, Events and Learning Team at Christchurch City Libraries were able to provide services to their community resulting in new innovative outreach programmes and a focus on promotion of online resources. On 19th December the “temporary” new and smaller Central library Peterborough opened. The retrieved Ngā Pounamu Māori and Ngai Tahu collections "Ngā rakau teitei e iwa”, have since been re-housed and are once again available for use by the public. Te Rūnanga o Ngāi Tahu. This organisation, established by the Te Rūnanga o Ngāi Tahu Act 1996, services the statutory rights for the people of Ngāi Tahu descent and ensures that the benefits of their Treaty Claim Settlement are enjoyed by Ngāi Tahu now and in the future. Ngāi Tahu are the indigenous Māori people of the southern islands of New Zealand - Te Waipounamu. The iwi (people) hold the rangatiratanga or tribal authority to over 80 per cent of the South Island. With their headquarters based in the central business they have also had to be relocated to temporary facilities. This included their library/archive collection of print resources, art works and taonga (cultural treasures).
Eccentrically Braced Frames (EBFs) are a widely used seismic resisting structural steel system. Since their inception in the late 1970s, they have been a viable option with an available stiffness that is between simple braced systems and moment resisting systems. A similar concept, the linked column frame (LCF), uses shear links between two closely spaced columns. In both cases, the key component is the active link or the shear link, and this component is the objective of this study. The performance of high rise EBF buildings in the 2010 and 2011 Christchurch earthquakes was beyond that which was expected, especially considering the very high accelerations recorded. As the concrete high-rises were torn down, two EBF buildings remained standing and only required some structural repair. These events prompted a renewed interest in bolted shear links, as well as their performance. While some research into replaceable shear links had already been done (Mansour, 2011), the objectives of this study were to improve on the shear link itself, with the consideration that links built in the future are likely to be bolted. The main components of this study were to: 1. Reduce or eliminate the requirements for intermediate web stiffeners, as they were suspected of being detrimental to performance. Furthermore, any reduction in stiffening requirements is a direct fabrication cost saving. Links with low web aspect ratios were found to achieve exceptional ductilities when no stiffeners were included, prompting new design equations. 2. Ensure that the stresses in the ends of links are adequately transferred into the endplates without causing fractures. Although most of the experimental links had web doubler plates included, four had varied lengths of such doubler plates from 0.0 in. to 8.0 in. The link without any doubler plates performed to a similar level to its peers, and thus it is likely that links with quality end details may not need web doubler plates at all. 3. Evaluate the performance of a link with double sided stiffeners without the use of web welds, as opposed to conventional single sided, welded stiffeners. This link performed well, and web-weld-less double sided stiffeners may be an economical alternative to conventional stiffeners for deeper sections of links. 4. Evaluate the performance of a link with thin endplates that are made efficient with the use of gusset plates. This link performed to an acceptable level and provides evidence for a cost effective alternative to thick endplates, especially considering the high overstrength end moments in links, typically requiring 16-bolt connections. 5. Examine the potential use of an alternative EBF arrangement where the collector beam is over sized, and the link section is formed by cutting out parts of the beam's web. After running a series of finite element models each with a unique variation, a number of approximate design rules were derived such that future research could develop this idea further experimentally. 6. Ensure that during testing, the secondary elements (members that are not the shear link), do not yield and are not close to yielding. None of the instrumented elements experienced any unexpected yielding, however the concerns for high stresses in the collector beam panel zone during design were warranted. The use of an existing New Zealand design equation is recommended as an extra check for design codes worldwide. The above objectives were mainly conducted experimentally, except: the data set for item 1 was greatly expanded through the use of a calibrated numerical model which was then used in an extensive parametric study; item 5 was purely finite element based; and, a small parametric study was included for item 3 in an attempt to expand on the trends found there.
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.
The structure and geomorphology of active orogens evolves on time scales ranging from a single earthquake to millions of years of tectonic deformation. Analysis of crustal deformation using new and established remote sensing techniques, and integration of these data with field mapping, geochronology and the sedimentary record, create new opportunities to understand orogenic evolution over these timescales. Timor Leste (East Timor) lies on the northern collisional boundary between continental crust from the Australian Plate and the Banda volcanic arc. GPS studies have indicated that the island of Timor is actively shortening. Field mapping and fault kinematic analysis of an emergent Pliocene marine sequence identifies gentle folding, overprinted by a predominance of NW-SE oriented dextral-normal faults and NE-SW oriented sinistral-normal faults that collectively bound large (5-20km2) bedrock massifs throughout the island. These fault systems intersect at non-Andersonian conjugate angles of approximately 120° and accommodate an estimated 20 km of orogen-parallel extension. Folding of Pliocene rocks in Timor may represent an early episode of contraction but the overall pattern of deformation is one of lateral crustal extrusion sub-parallel to the Banda Arc. Stratigraphic relationships suggest that extrusion began prior to 5.5 Ma, during and after initial uplift of the orogen. Sedimentological, geochemical and Nd isotope data indicate that the island of Timor was emergent and shedding terrigenous sediment into carbonate basins prior to 4.5 Ma. Synorogenic tectonic and sedimentary phases initiated almost synchronously across much of Timor Leste and <2 Myr before similar events in West Timor. An increase in plate coupling along this obliquely converging boundary, due to subduction of an outlying continental plateau at the Banda Trench, is proposed as a mechanism for uplift that accounts for orogen-parallel extension and early uplift of Timor Leste. Rapid bathymetric changes around Timor are likely to have played an important role in evolution of the Indonesian Seaway. The 2010 Mw 7.1 Darfield (Canterbury) earthquake in New Zealand was complex, involving multiple faults with strike-slip, reverse and normal displacements. Multi-temporal cadastral surveying and airborne light detection and ranging (LiDAR) surveys allowed surface deformation at the junction of three faults to be analyzed in this study in unprecedented detail. A nested, localized restraining stepover with contractional bulging was identified in an area with the overall fault structure of a releasing bend, highlighting the surface complexities that may develop in fault interaction zones during a single earthquake sequence. The earthquake also caused river avulsion and flooding in this area. Geomorphic investigations of these rivers prior to the earthquake identify plausible precursory patterns, including channel migration and narrowing. Comparison of the pre and post-earthquake geomorphology of the fault rupture also suggests that a subtle scarp or groove was present along much of the trace prior to the Darfield earthquake. Hydrogeology and well logs support a hypothesis of extended slip history and suggests that that the Selwyn River fan may be infilling a graben that has accumulated late Quaternary vertical slip of <30 m. Investigating fault behavior, geomorphic and sedimentary responses over a multitude of time-scales and at different study sites provides insights into fault interactions and orogenesis during single earthquakes and over millions of years of plate boundary deformation.
Following the Mw 6.2 Christchurch Earthquake on 22 February 2011, extensive ground cracking in loessial soils was reported in some areas of the Port Hills, southeast of central Christchurch. This study was undertaken to investigate the mechanisms of earthquake-induced ground damage on the eastern side of the Hillsborough Valley. A zone of extensional cracking up to 40m wide and 600m long was identified along the eastern foot-slope, accompanied by compression features and spring formation at the toe of the slope. An engineering geological and geomorphological model was developed for the eastern Hillsborough Valley that incorporates geotechnical investigation data sourced from the Canterbury Geotechnical Database (CGD), the findings of trenching and seismic refraction surveying carried out for this research, and interpretation of historical aerial photographs. The thickness and extent of a buried peat swamp at the base of the slope was mapped, and found to coincide with significant compression features. Ground cracking was found to have occurred entirely within loess-colluvium and to follow the apices of pre-1920s tunnel-gully fan debris at the southern end of the valley. The ground-cracking on the eastern side of the Hillsborough Valley is interpreted to have formed through tensile failure of the loess-colluvium. Testing was carried out to determine the tensile strength of Port Hills loess colluvium as a function of water content and density, in order to better understand the occurrence and distribution of the observed ground cracking. A comprehensive review of the soil tensile strength testing literature was undertaken, from which a test methodology was developed. Results show remoulded loess-colluvium to possess tensile strength of 7 - 28 kPa across the range of tested moisture contents (10-15%) and dry densities (1650-1900kg/m3). A positive linear relationship was observed between tensile strength and dry density, and a negative linear relationship between moisture content and tensile strength. The observed ground damage and available geotechnical information (inclinometer and piezometer records provided by the Earthquake Commission) were together used to interpret the mechanism(s) of slope movement that occurred in the eastern Hillsborough Valley. The observed ground damage is characteristic of translational movement, but without the development of lateral release scarps, or a basal sliding surface - which was not located during drilling. It is hypothesised that shear displacement has been accommodated by multiple slip surfaces of limited extent within the upper 10m of the slope. Movement has likely occurred within near-saturated colluvial units that have lost strength during earthquake shaking. The eastern Hillsborough Valley is considered to be an ‘incipient translational slide’, as both the patterns of damage and shearing are consistent with the early stages of such slide development. Sliding block analysis was utilised to understand how the eastern Hillsborough Valley may perform in a future large magnitude earthquake. Known cumulative displacements of ~0.3m for eastern Hillsborough Valley during the 2010-2011 Canterbury Earthquake Sequence were compared with modelled slope displacements to back-analyse a lower-bound yield acceleration of 0.2 - 0.25g. Synthetic broadband modelling for future Alpine and Hope Fault earthquakes indicates PGAs of approximately 0.08g for soil sites in the Christchurch area, as such, slope movement is unlikely to be reactivated by an Alpine Fault or Hope Fault earthquake. This does not take into account the possible role of strength loss due to excess pore pressure that may occur during these future events.
Following the 22 February 2011, MW 6.2 earthquake located on a fault beneath the Port Hills of Christchurch, fissuring of up to several hundred metres in length was observed in the loess and loess-colluvium of foot-slope positions in north-facing valleys of the Port Hills. The fissuring was observed in all major valleys, occurred at similar low altitudes, showing a contour-parallel orientation and often accompanied by both lateral compression/extension features and spring formation in the valley floor below. Fissuring locations studied in depth included Bowenvale Valley, Hillsborough Valley, Huntlywood Terrace–Lucas Lane, Bridle Path Road, and Maffeys Road–La Costa Lane. Investigations into loess soil, its properties and mannerisms, as well as international examples of its failure were undertaken, including study of the Loess Plateau of China, the Teton Dam, and palaeo-fissuring on Banks Peninsula. These investigations lead to the conclusion that loess has the propensity to fail, often due to the infiltration of water, the presence of which can lead to its instantaneous disaggregation. Literature study and laboratory analysis of Port Hills loess concluded that is has the ability to be stable in steep, sub-vertical escarpments, and often has a sub-vertically jointed internal structure and has a peak shear strength when dry. Values for cohesion, c (kPa) and the internal friction angle, ϕ (degrees) of Port Hills loess were established. The c values for the 40 Rapaki Road, 3 Glenview Terrace loess samples were 13.4 kPa and 19.7 kPa, respectively. The corresponding ϕ values were thought unusually high, at 42.0° and 43.4°.The analysed loess behaved very plastically, with little or no peak strength visible in the plots as the test went almost directly to residual strength. A geophysics resistivity survey showed an area of low resistivity which likely corresponds to a zone of saturated clayey loess/loess colluvium, indicating a high water table in the area. This is consistent with the appearances of local springs which are located towards the northern end of each distinct section of fissure trace and chemical analysis shows that they are sourced from the Port Hills volcanics. Port Hills fissuring may be sub-divided into three categories, Category A, Category B, and Category C, each characterised by distinctive features of the fissures. Category A includes fissures which display evidence of, spring formation, tunnel-gullying, and lateral spreading-like behaviour or quasi-toppling. These fissures are several metres down-slope of the loess-bedrock interface, and are in valleys containing a loess-colluvium fill. Category B fissures are in wider valleys than those in Category A, and the valleys contain estuarine silty sediments which liquefied during the earthquake. Category C fissures occurred at higher elevations than the fissures in the preceding categories, being almost coincident with bedrock outcropping. It is believed that the mechanism responsible for causing the fissuring is a complex combination of three mechanisms: the trampoline effect, bedrock fracturing, and lateral spreading. These three mechanisms can be applied in varying degrees to each of the fissuring sites in categories A, B, and C, in order to provide explanation for the observations made at each. Toppling failure can describe the soil movement as a consequence of the a three causative mechanisms, and provides insight into the movement of the loess. Intra-loess water coursing and tunnel gullying is thought to have encouraged and exacerbated the fissuring, while not being the driving force per se. Incipient landsliding is considered to be the least likely of the possible fissuring interpretations.
Landslides are significant hazards, especially in seismically-active mountainous regions, where shaking amplified by steep topography can result in widespread landsliding. These landslides present not only an acute hazard, but a chronic hazard that can last years-to-decades after the initial earthquake, causing recurring impacts. The Mw 7.8 Kaikōura earthquake caused more than 20,000 landslides throughout North Canterbury and resulted in significant damage to nationally significant infrastructure in the coastal transport corridor (CTC), isolating Kaikōura from the rest of New Zealand. In the years following, ongoing landsliding triggered by intense rainfall exacerbated the impacts and slowed the recovery process. However, while there is significant research on co-seismic landslides and their initial impacts in New Zealand, little research has explored the evolution of co-seismic landslides and how this hazard changes over time. This research maps landslides annually between 2013 and 2021 to evaluate the changes in pre-earthquake, co-seismic and post-earthquake rates of landsliding to determine how landslide hazard has changed over this time. In particular, the research explores how the number, area, and spatial distribution of landslides has changed since the earthquake, and whether post-earthquake mitigation works have in any way affected the long-term landslide hazard. Mapping of landslides was undertaken using open-source, medium resolution Landsat-8 and Sentinel-2 satellite imagery, with landslides identified visually and mapped as single polygons that capture both the source zone and deposit. Three study areas with differing levels of post-earthquake mitigation are compared: (i) the northern CTC, where the majority of mitigation was in the form of active debris removal; (ii) the southern CTC, where mitigation was primarily via passive protection measures; and (iii) Mount Fyffe, which has had no mitigation works since the earthquake. The results show that despite similar initial impacts during the earthquake, the rate of recovery in terms of landslide rates varies substantially across the three study areas. In Mount Fyffe, the number and area of landslides could take 45 and 22 years from 2021 respectively to return to pre-earthquake levels at the current rate. Comparatively, in the CTC, it could take just 5 years and 3-4 years from 2021 respectively. Notably, the fastest recovery in terms of landslide rates in the CTC was primarily located directly along the transport network, whereas what little recovery did occur in Mount Fyffe appeared to follow no particular pattern. Importantly, recovery rates in the northern CTC were notably higher than in the southern CTC, despite greater co-seismic impacts in the former. Combined, these results suggest the active, debris removal mitigation undertaken in the northern CTC may have had the effect of dramatically reducing the time for landslide rates to return to pre-earthquake levels. The role of slope angle and slope aspect were explored to evaluate if these observations could be driven by local differences in topography. The Mount Fyffe study area has higher slope angles than the CTC as a whole and landslides predominantly occurred on slightly steeper slopes than in the CTC. This may have contributed to the longer recovery times for landsliding in Mount Fyffe due to greater gravitational instability, however the observed variations are minor compared to the differences in recovery rates. In terms of slope aspect, landslides in Mount Fyffe preferentially occurred on north- and south-facing slopes whereas landslides in the CTC preferred the east- and south-facing slopes. The potential role of these differences in landslide recovery remains unclear but may be related to the propagation direction of the earthquake and the tracking direction of post-earthquake ex-tropical cyclones. Finally, landslides in the CTC are observed to be moving further away from the transport network and the number of landslides impacting the CTC decreased significantly since the earthquake. Nevertheless, the potential for further landslide reactivation remains. Therefore, despite the recovery in the CTC, it is clear that there is still risk of the transport network being impacted by further landsliding, at least for the next 3-5 yrs.
Christchurch City Council (Council) is undertaking the Land Drainage Recovery Programme in order to assess the effects of the earthquakes on flood risk to Christchurch. In the course of these investigations it has become better understood that floodplain management should be considered in a multi natural hazards context. Council have therefore engaged the Jacobs, Beca, University of Canterbury, and HR Wallingford project team to investigate the multihazards in eastern areas of Christchurch and develop flood management options which also consider other natural hazards in that context (i.e. how other hazards contribute to flooding both through temporal and spatial coincidence). The study has three stages: Stage 1 Gap Analysis – assessment of information known, identification of gaps and studies required to fill the gaps. Stage 2 Hazard Studies – a gap filling stage with the studies identified in Stage 1. Stage 3 Collating, Optioneering and Reporting – development of options to manage flood risk. This present report is to document findings of Stage 1 and recommends the studies that should be completed for Stage 2. It has also been important to consider how Stage 3 would be delivered and the gaps are prioritised to provide for this. The level of information available and hazards to consider is extensive; requiring this report to be made up of five parts each identifying individual gaps. A process of identifying information for individual hazards in Christchurch has been undertaken and documented (Part 1) followed by assessing the spatial co-location (Part 2) and probabilistic presence of multi hazards using available information. Part 3 considers multi hazard presence both as a temporal coincidence (e.g. an earthquake and flood occurring at one time) and as a cascade sequence (e.g. earthquake followed by a flood at some point in the future). Council have already undertaken a number of options studies for managing flood risk and these are documented in Part 4. Finally Part 5 provides the Gap Analysis Summary and Recommendations to Council. The key findings of Stage 1 gap analysis are: - The spatial analysis showed eastern Christchurch has a large number of hazards present with only 20% of the study area not being affected by any of the hazards mapped. Over 20% of the study area is exposed to four or more hazards at the frequencies and data available. - The majority of the Residential Red Zone is strongly exposed to multiple hazards, with 86% of the area being exposed to 4 or more hazards, and 24% being exposed to 6 or more hazards. - A wide number of gaps are present; however, prioritisation needs to consider the level of benefit and risks associated with not undertaking the studies. In light of this 10 studies ranging in scale are recommended to be done for the project team to complete the present scope of Stage 3. - Stage 3 will need to consider a number of engineering options to address hazards and compare with policy options; however, Council have not established a consistent policy on managed retreat that can be applied for equal comparison; without which substantial assumptions are required. We recommend Council undertake a study to define a managed retreat framework as an option for the city. - In undertaking Stage 1 with floodplain management as the focal point in a multi hazards context we have identified that Stage 3 requires consideration of options in the context of economics, implementation and residual risk. Presently the scope of work will provide a level of definition for floodplain options; however, this will not be at equal levels of detail for other hazard management options. Therefore, we recommend Council considers undertaking other studies with those key hazards (e.g. Coastal Hazards) as a focal point and identifies the engineering options to address such hazards. Doing so will provide equal levels of information for Council to make an informed and defendable decision on which options are progressed following Stage 3.
The Canterbury Earthquake Sequence (CES) of 2010-2011 produced large seismic moments up to Mw 7.1. These large, near-to-surface (<15 km) ruptures triggered >6,000 rockfall boulders on the Port Hills of Christchurch, many of which impacted houses and affected the livelihoods of people within the impacted area. From these disastrous and unpredicted natural events a need arose to be able to assess the areas affected by rockfall events in the future, where it is known that a rockfall is possible from a specific source outcrop but the potential boulder runout and dynamics are not understood. The distribution of rockfall deposits is largely constrained by the physical properties and processes of the boulder and its motion such as block density, shape and size, block velocity, bounce height, impact and rebound angle, as well as the properties of the substrate. Numerical rockfall models go some way to accounting for all the complex factors in an algorithm, commonly parameterised in a user interface where site-specific effects can be calibrated. Calibration of these algorithms requires thorough field checks and often experimental practises. The purpose of this project, which began immediately following the most destructive rupture of the CES (February 22, 2011), is to collate data to characterise boulder falls, and to use this information, supplemented by a set of anthropogenic boulder fall data, to perform an in-depth calibration of the three-dimensional numerical rockfall model RAMMS::Rockfall. The thesis covers the following topics: • Use of field data to calibrate RAMMS. Boulder impact trails in the loess-colluvium soils at Rapaki Bay have been used to estimate ranges of boulder velocities and bounce heights. RAMMS results replicate field data closely; it is concluded that the model is appropriate for analysing the earthquake-triggered boulder trails at Rapaki Bay, and that it can be usefully applied to rockfall trajectory and hazard assessment at this and similar sites elsewhere. • Detailed analysis of dynamic rockfall processes, interpreted from recorded boulder rolling experiments, and compared to RAMMS simulated results at the same site. Recorded rotational and translational velocities of a particular boulder show that the boulder behaves logically and dynamically on impact with different substrate types. Simulations show that seasonal changes in soil moisture alter rockfall dynamics and runout predictions within RAMMS, and adjustments are made to the calibration to reflect this; suggesting that in hazard analysis a rockfall model should be calibrated to dry rather than wet soil conditions to anticipate the most serious outcome. • Verifying the model calibration for a separate site on the Port Hills. The results of the RAMMS simulations show the effectiveness of calibration against a real data set, as well as the effectiveness of vegetation as a rockfall barrier/retardant. The results of simulations are compared using hazard maps, where the maximum runouts match well the mapped CES fallen boulder maximum runouts. The results of the simulations in terms of frequency distribution of deposit locations on the slope are also compared with those of the CES data, using the shadow angle tool to apportion slope zones. These results also replicate real field data well. Results show that a maximum runout envelope can be mapped, as well as frequency distribution of deposited boulders for hazard (and thus risk) analysis purposes. The accuracy of the rockfall runout envelope and frequency distribution can be improved by comprehensive vegetation and substrate mapping. The topics above define the scope of the project, limiting the focus to rockfall processes on the Port Hills, and implications for model calibration for the wider scientific community. The results provide a useful rockfall analysis methodology with a defensible and replicable calibration process, that has the potential to be applied to other lithologies and substrates. Its applications include a method of analysis for the selection and positioning of rockfall countermeasure design; site safety assessment for scaling and demolition works; and risk analysis and land planning for future construction in Christchurch.
Sewerage systems convey sewage, or wastewater, from residential or commercial buildings through complex reticulation networks to treatment plants. During seismic events both transient ground motion and permanent ground deformation can induce physical damage to sewerage system components, limiting or impeding the operability of the whole system. The malfunction of municipal sewerage systems can result in the pollution of nearby waterways through discharge of untreated sewage, pose a public health threat by preventing the use of appropriate sanitation facilities, and cause serious inconvenience for rescuers and residents. Christchurch, the second largest city in New Zealand, was seriously affected by the Canterbury Earthquake Sequence (CES) in 2010-2011. The CES imposed widespread damage to the Christchurch sewerage system (CSS), causing a significant loss of functionality and serviceability to the system. The Christchurch City Council (CCC) relied heavily on temporary sewerage services for several months following the CES. The temporary services were supported by use of chemical and portable toilets to supplement the damaged wastewater system. The rebuild delivery agency -Stronger Christchurch Infrastructure Rebuild Team (SCIRT) was created to be responsible for repair of 85 % of the damaged horizontal infrastructure (i.e., water, wastewater, stormwater systems, and roads) in Christchurch. Numerous initiatives to create platforms/tools aiming to, on the one hand, support the understanding, management and mitigation of seismic risk for infrastructure prior to disasters, and on the other hand, to support the decision-making for post-disaster reconstruction and recovery, have been promoted worldwide. Despite this, the CES in New Zealand highlighted that none of the existing platforms/tools are either accessible and/or readable or usable by emergency managers and decision makers for restoring the CSS. Furthermore, the majority of existing tools have a sole focus on the engineering perspective, while the holistic process of formulating recovery decisions is based on system-wide approach, where a variety of factors in addition to technical considerations are involved. Lastly, there is a paucity of studies focused on the tools and frameworks for supporting decision-making specifically on sewerage system restoration after earthquakes. This thesis develops a decision support framework for sewerage pipe and system restoration after earthquakes, building on the experience and learning of the organisations involved in recovering the CSS following the CES in 2010-2011. The proposed decision support framework includes three modules: 1) Physical Damage Module (PDM); 2) Functional Impact Module (FIM); 3) Pipeline Restoration Module (PRM). The PDM provides seismic fragility matrices and functions for sewer gravity and pressure pipelines for predicting earthquake-induced physical damage, categorised by pipe materials and liquefaction zones. The FIM demonstrates a set of performance indicators that are categorised in five domains: structural, hydraulic, environmental, social and economic domains. These performance indicators are used to assess loss of wastewater system service and the induced functional impacts in three different phases: emergency response, short-term recovery and long-term restoration. Based on the knowledge of the physical and functional status-quo of the sewerage systems post-earthquake captured through the PDM and FIM, the PRM estimates restoration time of sewer networks by use of restoration models developed using a Random Forest technique and graphically represented in terms of restoration curves. The development of a decision support framework for sewer recovery after earthquakes enables decision makers to assess physical damage, evaluate functional impacts relating to hydraulic, environmental, structural, economic and social contexts, and to predict restoration time of sewerage systems. Furthermore, the decision support framework can be potentially employed to underpin system maintenance and upgrade by guiding system rehabilitation and to monitor system behaviours during business-as-usual time. In conjunction with expert judgement and best practices, this framework can be moreover applied to assist asset managers in targeting the inclusion of system resilience as part of asset maintenance programmes.
