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.
This study investigates evidence for linkages and fault interactions centred on the Cust Anticline in Northwest Canterbury between Starvation Hill to the southwest and the Ashley and Loburn faults to the northeast. An integrated programme of geologic, geomorphic, paleo-seismic and geophysical analyses was undertaken owing to a lack of surface exposures and difficulty in distinguishing active tectonic features from fluvial and/or aeolian features across the low-relief Canterbury Plains. LiDAR analysis identified surface expression of several previously unrecognised active fault traces across the low-relief aggradation surfaces of the Canterbury Plains. Their presence is consistent with predictions of a fault relay exploiting the structural mesh across the region. This is characterised by interactions of northeast-striking contractional faults and a series of re-activating inherited Late Cretaceous normal faults, the latter now functioning as E–W-striking dextral transpressive faults. LiDAR also allowed for detailed analysis of the surface expression of individual faults and folds across the Cust Anticline contractional restraining bend, which is evolving as a pop-up structure within the newly established dextral shear system that is exploiting the inherited, now re-activated, basement fault zone. Paleo-seismic trenches were located on the crest of the western arm of the Cust Anticline and across a previously unrecognised E–W-striking fault trace, immediately southwest of the steeply plunging Cust Anticline termination. These studies confirmed the location and structural style of north-northeast-striking faults and an E–W-striking fault associated with the development of this structural culmination. A review of available industry seismic reflection lines emphasised the presence of a series of common structural styles having the same underlying structural drivers but with varying degrees of development and expression, both in the seismic profiles and in surface elevations across the study area. Based on LiDAR surface mapping and preliminary re-analysis of industry seismic reflection data, four fault zones are identified across the restraining bend structural culminations, which together form the proposed Oxford–Cust–Ashley Fault System. The 2010–2012 Canterbury Earthquake Sequence showed many similarities to the structural pattern established across the Oxford–Cust–Ashley Fault System, emphasising the importance of identification and characterization of presently hidden fault sources, and the understanding of fault network linkages, in order to improve constraints on earthquake source potential. Improved understanding of potentially-interactive fault sources in Northwest Canterbury, with the potential for combined initial fault rupture and spatial and temporal rupture propagation across this fault system, can be used in probabilistic seismic hazard analysis for the region, which is essential for the suitability and sustainability of future social and economic development.
A major hazard accompanying earthquake shaking in areas of steep topography is the detachment of rocks from bedrock outcrops that subsequently slide, roll, or bounce downslope (i.e. rockfalls). The 2010-2011 Canterbury earthquake sequence caused recurrent and severe rockfall in parts of southern Christchurch. Coseismic rockfall caused five fatalities and significant infrastructural damage during the 2011 Mw 6.2 Christchurch earthquake. Here we examine a rockfall site in southern Christchurch in detail using geomorphic mapping, lidar analysis, geochronology (cosmogenic 3He dating, radiocarbon dating, optically stimulated luminescence (OSL) from quartz, infrared stimulated luminescence from K-feldspar), numerical modeling of rockfall boulder trajectories, and ground motion prediction equations (GMPEs). Rocks fell from the source cliff only in earthquakes with interpolated peak ground velocities exceeding ~10 cm/s; hundreds of smaller earthquakes did not produce rockfall. On the basis of empirical observations, GMPEs and age chronologies we attribute paleo-rockfalls to strong shaking in prehistoric earthquakes. We conclude that earthquake shaking of comparable intensity to the strongest contemporary earthquakes in Christchurch last occurred at this site approximately 5000 to 7000 years ago, and that in some settings, rockfall deposits provide useful proxies for past strong ground motions.
Following the devastating 1931 Hawke's Bay earthquake, buildings in Napier and surrounding areas in the Hawke's Bay region were rebuilt in a comparatively homogenous structural and architectural style comprising the region's famous Art Deco stock. These interwar buildings are most often composed of reinforced concrete two-way space frames, and although they have comparatively ductile detailing for their date of construction, are often expected to be brittle, earthquake-prone buildings in preliminary seismic assessments. Furthermore, the likelihood of global collapse of an RC building during a design-level earthquake became an issue warranting particular attention following the collapse of multiple RC buildings in the February 22, 2011 Christchurch earthquake. Those who value the architectural heritage and future use of these iconic Art Deco buildings - including building owners, tenants, and city officials, among others - must consider how they can be best preserved and utilized functionally given the especially pressing implications of relevant safety, regulatory, and economic factors. This study was intended to provide information on the seismic hazard, geometric weaknesses, collapse hazards, material properties, structural detailing, empirically based vulnerability, and recommended analysis approaches particular to Art Deco buildings in Hawke's Bay as a resource for professional structural engineers tasked with seismic assessments and retrofit designs for these buildings. The observed satisfactory performance of similar low-rise, ostensibly brittle RC buildings in other earthquakes and the examination of the structural redundancy and expected column drift capacities in these buildings, led to the conclusion that the seismic capacity of these buildings is generally underrated in simple, force-based assessments.
There is very little research on total house strength that includes contributions of non-structural elements. This testing programme provides inclusive stiffness and response data for five houses of varying ages. These light timber framed houses in Christchurch, New Zealand had minor earthquake damage from the 2011 earthquakes and were lateral load tested on site to determine their strength and/or stiffness, and to identify damage thresholds. Dynamic characteristics including natural periods, which ranged from 0.14 to 0.29s were also investigated. Two houses were quasi-statically loaded up to approximately 130kN above the foundation in one direction. Another unidirectional test was undertaken on a slab-on-grade two-storey house, which was also snapback tested. Two other houses were tested using cyclic quasi-static loading, and between cycles snapback tests were undertaken to identify the natural period of each house, including foundation and damage effects. A more detailed dynamic analysis on one of the houses provided important information on seismic safety levels of post-quake houses with respect to different hazard levels in the Christchurch area. While compared to New Zealand Building Standards all tested houses had an excess of strength, damage is a significant consideration in earthquake resilience and was observed in all of the houses. http://www.aees.org.au/downloads/conference-papers/2015-2/
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 Sendai Framework for Disaster Risk Reduction 2015-2030 finds that, despite progress in disaster risk reduction over the last decade “evidence indicates that exposure of persons and assets in all countries has increased faster than vulnerability has decreased, thus generating new risk and a steady rise in disaster losses” (p.4, UNISDR 2015). Fostering cooperation among relevant stakeholders and policy makers to “facilitate a science-policy interface for effective decisionmaking in disaster risk management” is required to achieve two priority areas for action, understanding disaster risk and enhancing disaster preparedness (p. 13, p. 23, UNISDR 2015). In other topic areas, the term science-policy interface is used interchangeably with the term boundary organisation. Both terms are usually used refer to systematic collaborative arrangements used to manage the intersection, or boundary, between science and policy domains, with the aim of facilitating the joint construction of knowledge to inform decision-making. Informed by complexity theory, and a constructivist focus on the functions and processes that minimize inevitable tensions between domains, this conceptual framework has become well established in fields where large complex issues have significant economic and political consequences, including environmental management, biodiversity, sustainable development, climate change and public health. To date, however, there has been little application of this framework in the disaster risk reduction field. In this doctoral project the boundary management framework informs an analysis of the research response to the 2010-2011 Canterbury Earthquake Sequence, focusing on the coordination role of New Zealand’s national Natural Hazards Research Platform. The project has two aims. It uses this framework to tell the nuanced story of the way this research coordination role evolved in response to both the complexity of the unfolding post-disaster environment, and to national policy and research developments. Lessons are drawn from this analysis for those planning and implementing arrangements across the science-policy boundary to manage research support for disaster risk reduction decision-making, particularly after disasters. The second aim is to use this case study to test the utility of the boundary management framework in the disaster risk reduction context. This requires that terminology and concepts are explained and translated in terms that make this analysis as accessible as possible across the disciplines, domains and sectors involved in disaster risk reduction. Key findings are that the focus on balance, both within organisations, and between organisations and domains, and the emphasis on systemic effects, patterns and trends, offer an effective and productive alternative to the more traditional focus on individual or organisational performance. Lessons are drawn concerning the application of this framework when planning and implementing boundary organisations in the hazard and disaster risk management context.
