Generalized conditional intensity measure (GCIM) method is extended to ground motion selection for scenario ruptures. Using different rupture scenarios and site conditions, various aspects of the GCIM methodology are scrutinized, including: (i) implementation of different weight vectors and the composition of the IM vector; (ii) quantifying the importance of replicate selections for different number of desired ground motions; and (iii) the effect of considering bounds on the implicit causal parameters of the prospective ground motions. Using the extended methodology, representative ground motion ensembles for several major earthquake scenarios in New Zealand are developed. Cases considered include representative ground motions for the occurrence of Alpine, Hope, and Porters Pass earthquakes in Christchurch city, and the occurrence of Wellington, Wairarapa, and Ohariu fault ruptures in Wellington city. Challenges in the development of ground motion ensembles for subduction zone earthquakes are also highlighted. The selected scenario-based ground motion sets can be used to complement ground motions which are often selected in conjunction with probabilistic seismic hazard analysis, in order to understand the performance of structures for the question “what if this fault ruptures?”
Generalized conditional intensity measure (GCIM) method is extended to ground motion selection for scenario ruptures. Using different rupture scenarios and site conditions, various aspects of the GCIM methodology are scrutinized, including: (i) implementation of different weight vectors and the composition of the IM vector; (ii) quantifying the importance of replicate selections for different number of desired ground motions; and (iii) the effect of considering bounds on the implicit causal parameters of the prospective ground motions. Using the extended methodology, representative ground motion ensembles for several major earthquake scenarios in New Zealand are developed. Cases considered include representative ground motions for the occurrence of Alpine, Hope, and Porters Pass earthquakes in Christchurch city, and the occurrence of Wellington, Wairarapa, and Ohariu fault ruptures in Wellington city. Challenges in the development of ground motion ensembles for subduction zone earthquakes are also highlighted. The selected scenario-based ground motion sets can be used to complement ground motions which are often selected in conjunction with probabilistic seismic hazard analysis, in order to understand the performance of structures for the question “what if this fault ruptures?”
In 2016, the Building (Earthquake-prone Buildings) Amendment Act 2016 was introduced to address the issue of seismic vulnerability amongst existing buildings in Aotearoa New Zealand. This Act introduced a mandatory scheme to remediate buildings deemed particularly vulnerable to seismic hazard, as recommended by the 2012 Royal Commission into the Canterbury earthquake sequence of 2010–2011. This Earthquake-prone Building (EPB) framework is unusual internationally for the mandatory obligations that it introduces. This article explores and critiques the operation of the scheme in practice through an examination of its implementation provisions and the experiences of more recent seismic events (confirmed by engineering research). This analysis leads to the conclusion that the operation of the current scheme and particularly the application of the concept of EPB vulnerability excludes large numbers of (primarily urban) buildings which pose a significant risk in the event of a significant (but expected) seismic event. As a result, the EPB scheme fails to achieve its goals and instead may create a false impression that it does so
Liquefaction-induced lateral spreading in Christchurch and surrounding suburbs during the recent Canterbury Earthquake Sequence (2010-2011) caused significant damage to structures and lifelines located in close proximity to streams and rivers. Simplified methods used in current engineering practice for predicting lateral ground displacements exhibit a high degree of epistemic uncertainty, but provide ‘order of magnitude’ estimates to appraise the hazard. We wish to compare model predictions to field measurements in order to assess the model’s capabilities and limitations with respect to Christchurch conditions. The analysis presented focuses on the widely-used empirical model of Youd et al. (2002), developed based on multi-linear regression (MLR) of case history data from lateral spreading occurrence in Japan and the US. Two issues arising from the application of this model to Christchurch were considered: • Small data set of Standard Penetration Test (SPT) and soil gradation indices (fines content FC, and mean grain size, D50) required for input. We attempt to use widely available CPT data with site specific correlations to FC and D50. • Uncertainty associated with the model input parameters and their influence on predicted displacements. This has been investigated for a specific location through a sensitivity analysis.
After a high-intensity seismic event, inspections of structural damages need to be carried out as soon as possible in order to optimize the emergency management, as well as improving the recovery time. In the current practice, damage inspections are performed by an experienced engineer, who physically inspect the structures. This way of doing not only requires a significant amount of time and high skilled human resources, but also raises the concern about the inspector’s safety. A promising alternative is represented using new technologies, such as drones and artificial intelligence, which can perform part of the damage classification task. In fact, drones can safely access high hazard components of the structures: for instance, bridge piers or abutments, and perform the reconnaissance by using highresolution cameras. Furthermore, images can be automatically processed by machine learning algorithms, and damages detected. In this paper, the possibility of applying such technologies for inspecting New Zealand bridges is explored. Firstly, a machine-learning model for damage detection by performing image analysis is presented. Specifically, the algorithm was trained to recognize cracks in concrete members. A sensitivity analysis was carried out to evaluate the algorithm accuracy by using database images. Depending on the confidence level desired,i.e. by allowing a manual classification where the alghortim confidence is below a specific tolerance, the accuracy was found reaching up to 84.7%. In the second part, the model is applied to detect the damage observed on the Anzac Bridge (GPS coordinates -43.500865, 172.701138) in Christchurch by performing a drone reconnaissance. Reults show that the accuracy of the damage detection was equal to 88% and 63% for cracking and spalling, respectively.
As a result of the 4 September 2010 Darfield earthquake and the more damaging 22 February 2011 Christchurch earthquake, considerable damage occurred to a significant number of buildings in Christchurch. The damage that occurred to the Christchurch Roman Catholic Cathedral of the Blessed Sacrament (commonly known as the Christchurch Basilica) as a result of the Canterbury earthquakes is reported, and the observed failure modes are identified. A previous strengthening intervention is outlined and the estimated capacity of the building is discussed. This strengthening was completed in 2004, and addressed the worst aspects of the building's seismic vulnerability. Urgent work was undertaken post-earthquake to secure parts of the building in order to limit damage and prevent collapse of unstable parts of the building. The approach taken for this securing is outlined, and the performance of the building and the previously installed earthquake strengthening intervention is evaluated.A key consideration throughout the project was the interaction between the structural securing requirements that were driven by the requirement to limit damage and mitigate hazards, and the heritage considerations. Lessons learnt from the strengthening that was carried out, the securing work undertaken, and the approach taken in making the building "safe" are discussed. Some conclusions are drawn with respect to the effectiveness of strengthening similar building types, and the approach taken to secure the building under active seismic conditions. AM - Accepted Manuscript
In the last two decades, New Zealand (NZ) has experienced significant earthquakes, including the 2010 M 7.2 Darfield, 2011 M 6.2 Christchurch, and 2016 M 7.8 Kaikōura events. Amongst these large events, tens of thousands of smaller earthquakes have occurred. While previous event and ground-motion databases have analyzed these events, many events below M 4 have gone undetected. The goal of this study is to expand on previous databases, particularly for small magnitude (M<4) and low-amplitude ground motions. This new database enables a greater understanding of regional variations within NZ and contributes to the validity of internationally developed ground-motion models. The database includes event locations and magnitude estimates with uncertainty considerations, and tectonic type assessed in a hierarchical manner. Ground motions are extracted from the GeoNet FDSN server and assessed for quality using a neural network classification approach. A deep neural network approach is also utilized for picking P and S phases for determination of event hypocentres. Relative hypocentres are further improved by double-difference relocation and will contribute toward developing shallow (< 50 km) seismic tomography models. Analysis of the resulting database is compared with previous studies for discussion of implications toward national hazard prediction models.
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/
Ground motion observations from the most significant 10 events in the 2010-2011 Canterbury earthquake sequence at near-source sites are utilized to scrutinize New Zealand (NZ)-specific pseudo-spectral acceleration (SA) empirical ground motion prediction equations (GMPE) (Bradley 2010, Bradley 2013, McVerry et al. 2006). Region-specific modification factors based on relaxing the conventional ergodic assumption in GMPE development were developed for the Bradley (2010) model. Because of the observed biases with magnitude and source-to-site distance for the McVerry et al. (2006) model it is not possible to develop region-specific modification factors in a reliable manner. The theory of non-ergodic empirical ground motion prediction is then outlined, and applied to this 10 event dataset to determine systematic effects in the between- and within-event residuals which lead to modifications in the predicted median and standard deviation of the GMPE. By examining these systematic effects over sub-regions containing a total of 20 strong motion stations within the Canterbury area, modification factors for use in region-specific ground motion prediction are proposed. These modification factors, in particular, are suggested for use with the Bradley et al. (2010) model in Canterbury-specific probabilistic seismic hazard analysis (PSHA) to develop revised design response, particularly for long vibration periods.
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.
This study explicitly investigates uncertainties in physics-based ground motion simulation validation for earthquakes in the Canterbury region. The simulations utilise the Graves and Pitarka (2015) hybrid methodology, with separately quantified parametric uncertainties in the comprehensive physics and simplified physics components of the model. The study is limited to the simulation of 148 small magnitude (Mw 3.5 – 5) earthquakes, with a point source approximation for the source rupture representations, which also enables a focus on a small number of relevant uncertainties. The parametric uncertainties under consideration were selected through sensitivity analysis, and specifically include: magnitude, Brune stress parameter and high frequency rupture velocity. Twenty Monte Carlo realisations were used to sample parameter uncertainties for each of the 148 events. Residuals associated with the following intensity measures: spectral acceleration, peak ground velocity, arias intensity and significant duration, were ascertained. Using these residuals, validation was performed through assessment of systematic biases in site and source terms from mixed-effects regression. Based on the results to date, initial standard deviation recommendations for parameter uncertainties, based on the Canterbury simulations have been obtained. This work ultimately provides an initial step toward explicit incorporation of modelling uncertainty in simulated ground motion predictions for future events, which will improve the use of simulation models in seismic hazard analysis. We plan to subsequently assess uncertainties for larger magnitude events with more complex ruptures, and events across a larger geographic region, as well as uncertainties due to path attenuation, site effects, and more general model epistemic uncertainties.
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.
Liquefaction-induced lateral spreading in large seismic events often results in pervasive and costly damage to engineering structures and lifelines, making it a critical component of engineering design. However, the complex nature of this phenomenon leads to designing for such a hazard extremely challenging and there is a clear for an improved understanding and predicting liquefaction-induced lateral spreading. The 2010-2011 Canterbury (New Zealand) Earthquakes triggered severe liquefaction-induced lateral spreading along the streams and rivers of the Christchurch region, causing extensive damage to roads, bridges, lifelines, and structures in the vicinity. The unfortunate devastation induced from lateral spreading in these events also rendered the rare opportunity to gain an improved understanding of lateral spreading displacements specific to the Christchurch region. As part of this thesis, the method of ground surveying was employed following the 4 September 2010 Darfield (Mw 7.1) and 22 February 2011 Christchurch (Mw 6.2) earthquakes at 126 locations (19 repeated) throughout Christchurch and surrounding suburbs. The method involved measurements and then summation of crack widths along a specific alignment (transect) running approximately perpendicular to the waterway to indicate typically a maximum lateral displacement at the bank and reduction of the magnitude of displacements with distance from the river. Rigorous data processing and comparisons with alternative measurements of lateral spreading were performed to verify results from field observations and validate the method of ground surveying employed, as well as highlight the complex nature of lateral spreading displacements. The welldocumented field data was scrutinized to gain an understanding of typical magnitudes and distribution patterns (distribution of displacement with distance) of lateral spreading observed in the Christchurch area. Maximum displacements ranging from less than 10 cm to over 3.5 m were encountered at the sites surveyed and the area affected by spreading ranged from less than 20 m to over 200 m from the river. Despite the highly non-uniform displacements, four characteristic distribution patterns including large, distributed ground displacements, block-type movements, large and localized ground displacements, and areas of little to no displacements were identified. Available geotechnical, seismic, and topographic data were collated at the ground surveying sites for subsequent analysis of field measurements. Two widely-used empirical models (Zhang et al. (2004), Youd et al. (2002)) were scrutinized and applied to locations in the vicinity of field measurements for comparison with model predictions. The results indicated generally poor correlation (outside a factor of two) with empirical predictions at most locations and further validated the need for an improved, analysis- based method of predicting lateral displacements that considers the many factors involved on a site-specific basis. In addition, the development of appropriate model input parameters for the Youd et al. (2002) model led to a site-specific correlation of soil behavior type index, Ic, and fines content, FC, for sites along the Avon River in Christchurch that matched up well with existing Ic – FC relationships commonly used in current practice. Lastly, a rigorous analysis was performed for 25 selected locations of ground surveying measurements along the Avon River where ground slope conditions are mild (-1 to 2%) and channel heights range from about 2 – 4.5 m. The field data was divided into categories based on the observed distribution pattern of ground displacements including: large and distributed, moderate and distributed, small to negligible, and large and localized. A systematic approach was applied to determine potential critical layers contributing to the observed displacement patterns which led to the development of characteristic profiles for each category considered. The results of these analyses outline an alternative approach to the evaluation of lateral spreading in which a detailed geotechnical analysis is used to identify the potential for large spreading displacements and likely spatial distribution patterns of spreading. Key factors affecting the observed magnitude and distribution of spreading included the thickness of the critical layer, relative density, soil type and layer continuity. It was found that the large and distributed ground displacements were associated with a thick (1.5 – 2.5 m) deposit of loose, fine to silty sand (qc1 ~4-7 MPa, Ic 1.9-2.1, qc1n_cs ~50-70) that was continuous along the bank and with distance from the river. In contrast, small to negligible displacements were characterized by an absence of or relatively thin (< 1 m), discontinuous critical layer. Characteristic features of the moderate and distributed displacements were found to be somewhere between these two extremes. The localized and large displacements showed a characteristic critical layer similar to that observed in the large and distributed sites but that was not continuous and hence leading to the localized zone of displacement. The findings presented in this thesis illustrate the highly complex nature of lateral displacements that cannot be captured in simplified models but require a robust geotechnical analysis similar to that performed for this research.
Fine grained sediment deposition in urban environments during natural hazard events can impact critical infrastructure and properties (urban terrain) leading to reduced social and economic function and potentially adverse public health effects. Therefore, clean-up of the sediments is required to minimise impacts and restore social and economic functionality as soon as possible. The strategies employed to manage and coordinate the clean-up significantly influence the speed, cost and quality of the clean-up operation. Additionally, the physical properties of the fine grained sediment affects the clean-up, transport, storage and future usage of the sediment. The goals of the research are to assess the resources, time and cost required for fine grained sediment clean-up in an urban environment following a disaster and to determine how the geotechnical properties of sediment will affect urban clean-up strategies. The thesis focuses on the impact of fine grained sediment (<1 mm) deposition from three liquefaction events during the Canterbury earthquake sequence (2010-2011) on residential suburbs and transport networks in Christchurch. It also presents how geotechnical properties of the material may affect clean-up strategies and methods by presenting geotechnical analysis of tephra material from the North Island of New Zealand. Finally, lessons for disaster response planning and decision making for clean-up of sediment in urban environments are presented. A series of semi-structured interviews of key stakeholders supported by relevant academic literature and media reports were used to record the clean-up operation coordination and management and to make a preliminary qualification of the Christchurch liquefaction ejecta clean-up (costs breakdown, time, volume, resources, coordination, planning and priorities). Further analysis of the costs and resources involved for better accuracy was required and so the analysis of Christchurch City Council road management database (RAMM) was done. In order to make a transition from general fine sediment clean-up to specific types of fine disaster sediment clean-up, adequate information about the material properties is required as they will define how the material will be handled, transported and stored. Laboratory analysis of young volcanic tephra from the New Zealand’s North Island was performed to identify their geotechnical properties (density, granulometry, plasticity, composition and angle of repose). The major findings of this research were that emergency planning and the use of the coordinated incident management system (CIMS) system during the emergency were important to facilitate rapid clean-up tasking, management of resources and ultimately recovery from widespread and voluminous liquefaction ejecta deposition in eastern Christchurch. A total estimated cost of approximately $NZ 40 million was calculated for the Christchurch City clean-up following the 2010-2011 Canterbury earthquake sequence with a partial cost of $NZ 12 million for the Southern part of the city, where up to 33% (418 km) of the road network was impacted by liquefaction ejecta and required clearing of the material following the 22 February 2011 earthquake. Over 500,000 tonnes of ejecta has been stockpiled at Burwood landfill for all three liquefaction inducing earthquake events. The average cost per kilometre for the event clean-up was $NZ 5,500/km (4 September 2010), $NZ 11,650/km (22 February 2011) and $NZ 11,185/km (13 June 2011). The duration of clean-up time of residential properties and the road network was approximately two to three months for each of the three liquefaction ejecta events; despite events volumes and spatial distribution of ejecta. Interviews and quantitative analysis of RAMM data revealed that the experience and knowledge gained from the Darfield earthquake (4 September 2010) clean-up increased the efficiency of the following Christchurch earthquake induced liquefaction ejecta clean-up events. Density, particle size, particle shape, clay content and moisture content, are the important geotechnical properties that need to be considered when planning for a clean-up method that incorporates collection, transport and disposal or storage. The geotechnical properties for the tephra samples were analysed to increase preparedness and reaction response of potentially affected North Island cities from possible product from the active volcanoes in their region. The geotechnical results from this study show that volcanic tephra could be used in road or construction material but the properties would have to be further investigated for a New Zealand context. Using fresh volcanic material in road, building or flood control construction requires good understanding of the material properties and precaution during design and construction to extra care, but if well planned, it can be economically beneficial.
<strong>Sea level rise is one consequence of Earth’s changing climate. Century-long tide gauge records show that global-mean sea-level rise reached 11-16 cm during the twentieth century at a mean rate of 1.2 mm/y. Today, the average rate of global-mean sea-level rise is higher at 3-4 mm/y and is expected to increase in the future. This represents a hazard to low elevation coastal zones worldwide. Yet, before global sea level projections can be used to characterise future coastal flood hazard at a local scale, the effects of tectonics (and other processes) that drive vertical land motion (VLM) must be considered. VLM is defined as the vertical velocity (uplift or subsidence) of the solid surface with respect to the centre of Earth. In this study, new VLM maps are generated over coastal strips in New Zealand, using Sentinel-1 InSAR and GNSS data.</strong>In New Zealand, measuring VLM using InSAR on naturally vegetated or agricultural land is difficult due to signal decorrelation. Along the rural Bay of Plenty coastal strip, I use a persistent-scatterer approach to generate a VLM map from both east-looking ascending and west-looking descending Sentinel-1 data between 2015-2021. Using time-series data over the same time period from a dense network of 20 GNSS sensors, I tie InSAR-derived line-of-sight velocity to the 2014 ITRF reference frame. I test two different methods for measuring VLM and compare the results against GNSS vertical velocity along the Bay of Plenty coast. Best results are achieved by first removing the interpolated horizontal GNSS velocity field from each of the InSAR datasets, before averaging the two VLM estimates. Measured VLM is between -3 and 3 mm/y, with negative values (subsidence) occurring within the low-lying Rangitāiki Plain and Ōpōtiki valley, and uplift across the elevated region west of Matatā.This thesis integrates geomorphological, geological, and historical levelling VLM records with modern satellite datasets to assess VLM across timescales ranging from 10 to 100,000 years at Matatā. Uplift rate has been variable through time, with average uplift over the last 300,000 years of 1 mm/y, 4.5 mm/y since 1720 years, 2 mm/y between 1950-1978, and 10 mm/y between 2004-2011. Previous modelling has shown that the best fit to the 2004-2011 rapid uplift rates is an inflating magmatic source at ~10 km depth beneath Matatā. To reconcile all data, I present a VLM model that consists of short-lived periods (7 years) of rapid uplift (10 mm/y), separated by longer periods (30 years) of lower background uplift (3 mm/y). The episodic nature of VLM at Matatā likely reflects short-lived periods of magmatic intrusion. Episodic VLM characterised by large rates of uplift (10 mm/y) has been seen at Taupō volcano, and other volcanic centers globally. It has been 12 years since the end of the last intrusion episode; this modelling suggest one may expect to observe increased uplift rates at Matatā in the coming decades. Densely populated urban coastal strips are most at risk from the effects of relative sea-level rise. At the same time, anthropogenic activities associated with urbanization, such as groundwater withdrawal, and land reclamation can lead to local land subsidence (LLS), further exacerbating the risk to urban infrastructure. LLS refers to subsidence relative to nearby land area assumed to be stable. In this thesis, I create the first high-resolution (10 m) maps of LLS at six urban coastal strips in New Zealand, with a combined length of 285 km, using Sentinel-1 InSAR data between 2018-2021. This analysis reveals 89% of urban coastal strips are subsiding at rates of -0.5 mm/y or greater, and 11% is subsiding at higher rates of -3.0 mm/y or greater. On average, subsidence is -0.6 to -2.9 mm/y higher at the coastal strip, compared to inland areas occupied by GNSS stations. This analysis also documents highly-localised hotspots of LLS, with subsidence rates of up to -15 mm/y. In Christchurch, rapid and localised subsidence (-8 mm/y) is observed within coastal suburbs New Brighton and Southshore. In most cities, the highest subsidence rates occur on land reclaimed in the early-late twentieth century, and in areas built on Holocene sediment. Time-series analysis of LLS at sites of reclaimed land shows both linear and non-linear rates of deformation over time periods of up to 6-8 years. This thesis highlights the variable exposure to relative sea-level rise of New Zealand coastal strips, and demonstrates that in many cases current rates of VLM should be expected to continue for the next few decades.
Liquefaction during the 4th September 2010 Mw 7.1 Darfield earthquake and large aftershocks in 2011 (Canterbury earthquake sequence, CES) caused severe damage to land and infrastructure within Christchurch, New Zealand. Approximately one third of the total CES-induced financial losses were directly attributable to liq- uefaction and thus highlights the need for local and regional authorities to assess liquefaction hazards for present and future developments. This thesis is the first to conduct paleo-liquefaction studies in eastern Christchurch for the purpose of de- termining approximate return times of liquefaction-inducing earthquakes within the region. The research uncovered evidence for pre-CES liquefaction dated by radiocarbon and cross-cutting relationships as post-1660 to pre-1905. Additional paleo-liquefaction investigations within the eastern Christchurch suburb of Avon- dale, and the northern township of Kaiapoi, revealed further evidence for pre-CES liquefaction. Pre-CES liquefaction in Avondale is dated as post-1321 and pre-1901, while the Kaiapoi features likely formed during three distinct episodes: post-1458 and possibly during the 1901 Cheviot earthquake, post-1297 to pre-1901, and pre-1458. Evaluation of the liquefaction potential of active faults within the Can- terbury region indicates that many faults have the potential to cause widespread liquefaction within Avondale and Kaiapoi. The identification of pre-CES liquefac- tion confirms that these areas have previously liquefied, and indicates that residen- tial development in eastern Christchurch between 1860 and 2005 occurred in areas containing geologic evidence for pre-CES liquefaction. Additionally, on the basis of detailed field and GIS-based mapping and geospatial-statistical analysis, the distribution and severity of liquefaction and lateral spreading within the eastern Christchurch suburb of Avonside is shown in this study to be strongly in uenced by geomorphic and topographic variability. This variability is not currently ac- counted for in site-specific liquefaction assessments nor the simplified horizontal displacement models, and accounts for some of the variability between the pre- dicted horizontal displacements and those observed during the CES. This thesis highlights the potential applications of paleo-liquefaction investigations and ge- omorphic mapping to seismic and liquefaction hazard assessments and may aid future land-use planning decisions.
During many years the analysis of some geophysical results of Charles Darwin was being carried out in Department. Darwin has connected almost 200 years ago results of catastrophic earthquakes with vertical movement of a surface of the Earth. Usually this movement less horizontal movement and its influence on destruction of cities is not considered. Earthquake hazard assessment studies were focused usually on the horizontal ground motion. Effects of the strong vertical motion were not, practically, discussed. The margins of safety against gravity-induced static vertical forces in constructed buildings usually provide adequate resistance to dynamic forces induced by the vertical acceleration during an earthquake. However, the earthquake in Christchurch is an example of the vertical seismic shock . The earthquake magnitude was rather small - nearby 6.3. However, the result was catastrophic. The same took place in 1835. It allowed to Darwin to formulate a few great ideas. Charles Darwin has explained qualitatively results of an interaction of huge seismic waves with volcanoes and the nature of volcanism and seismicity of our planet. These important data of Charles Darwin became very actual recently. It is possible to tell also the same about tsunami and extreme ocean waves described by Charles Darwin. Therefore this data were analyzed using modern mechanics, mathematics and physics in Department. In particular, the theory of catastrophic waves was developed based on Darwin's data. The theory tried to explain occurrence, evolution and distribution the catastrophic waves in various natural systems, since atoms, oceans, surfaces of the Earth and up to the very early Universe. Some results of the research were published in prestigious magazines. Later they were presented in two books devoted to Charles Darwin's anniversary (2009). Last from them was published in Russian (2011). We give here key ideas of this research which is a part of interdisciplinary researches of Department. Some ideas are discussed. Not less important purpose is very short historical review of some researches of Darwin. In particular, we underline Darwin' priority in the formulation of the bases of Dynamics Earth.
Probabilistic Structural Fire Engineering (PSFE) has been introduced to overcome the limitations of current conventional approaches used for the design of fire-exposed structures. Current structural fire design investigates worst-case fire scenarios and include multiple thermal and structural analyses. PSFE permits buildings to be designed to a level of life safety or economic loss that may occur in future fire events with the help of a probabilistic approach. This thesis presents modifications to the adoption of a Performance-Based Earthquake Engineering (PBEE) framework in Probabilistic Structural Fire Engineering (PSFE). The probabilistic approach runs through a series of interrelationships between different variables, and successive convolution integrals of these interrelationships result in probabilities of different measures. The process starts with the definition of a fire severity measure (FSM), which best relates fire hazard intensity with structural response. It is identified by satisfying efficiency and sufficiency criteria as described by the PBEE framework. The relationship between a fire hazard and corresponding structural response is established by analysis methods. One method that has been used to quantify this relationship in PSFE is Incremental Fire Analysis (IFA). The existing IFA approach produces unrealistic fire scenarios, as fire profiles may be scaled to wide ranges of fire severity levels, which may not physically represent any real fires. Two new techniques are introduced in this thesis to limit extensive scaling. In order to obtain an annual rate of exceedance of fire hazard and structural response for an office building, an occurrence model and an attenuation model for office fires are generated for both Christchurch city and New Zealand. The results show that Christchurch city is 15% less likely to experience fires that have the potential to cause structural failures in comparison to all of New Zealand. In establishing better predictive relationships between fires and structural response, cumulative incident radiation (a fire hazard property) is found to be the most appropriate fire severity measure. This research brings together existing research on various sources of uncertainty in probabilistic structural fire engineering, such as elements affecting post-flashover fire development factors (fuel load, ventilation, surface lining and compartment geometry), fire models, analysis methods and structural reliability. Epistemic uncertainty and aleatory uncertainty are investigated in the thesis by examining the uncertainty associated with modelling and the factors that influence post-flashover development of fires. A survey of 12 buildings in Christchurch in combination with recent surveys in New Zealand produced new statistical data on post-flashover development factors in office buildings in New Zealand. The effects of these parameters on temperature-time profiles are evaluated. The effects of epistemic uncertainty due to fire models in the estimation of structural response is also calculated. Parametric fires are found to have large uncertainty in the prediction of post-flashover fires, while the BFD curves have large uncertainties in prediction of structural response. These uncertainties need to be incorporated into failure probability calculations. Uncertainty in structural modelling shows that the choices that are made during modelling have a large influence on realistic predictions of structural response.
On 14 November 2016, the Mw 7.8 Kaikōura earthquake caused widespread damage along the east coast of the South Island, New Zealand. Kaikōura town itself was isolated from the rest of the country by landslides blocking off major roads. While impacts from the Kaikōura earthquake on large, urban population centres have been generally well documented, this thesis aims to fill gaps in academic knowledge regarding small rural towns. This thesis investigates what, where and when critical infrastructure and lifeline service disruption occurred following the 2016 Kaikōura earthquake in a selection of small towns, and how the communities in these areas adapted to disruption. Following a robust review of literature and news media, four small rural towns were selected from North Canterbury (Culverden & Waiau) and Marlborough (Seddon & Ward) in the South Island, New Zealand. Semi-structured interview sessions with a special focus on these towns were held with infrastructure managers, emergency response and recovery officials, and organisation leaders with experience or expertise in the 2016 Kaikōura earthquake. Findings were supplemented with emergency management situation reports to produce hazard maps and infrastructure exposure maps. A more detailed analysis was conducted for Waiau involving interdependence analyses and a level of service timeline for select lifeline services. The earthquake impacted roads by blocking them with landslides, debris and surface rupture. Bridges where shaken off their abutments, breaking infrastructure links such as fibre landlines as they went. Water supplies and other forms of infrastructure relied heavily on the level of service of roads, as rough rural terrain left few alternatives. Adapting to an artificial loss of road service, some Waiau locals created their own detour around a road cordon in order to get home to family and farms. Performance of dwellings was tied to socioeconomic factors as much as proximity to the epicentre. Farmers who lost water access pulled out fences to allow stock to drink from rivers. Socioeconomic differences between farmland and township residents also contributed to resilience variations between the towns assessed in this study. Understanding how small rural towns respond and adapt to disaster allows emergency management officials and policy to be well informed and flexible with planning for multiple size classes of towns.
Depending on their nature and severity, disasters can create large volumes of debris and waste. Waste volumes from a single event can be the equivalent of many times the annual waste generation rate of the affected community. These volumes can overwhelm existing solid waste management facilities and personnel. Mismanagement of disaster waste can affect both the response and long term recovery of a disaster affected area. Previous research into disaster waste management has been either context specific or event specific, making it difficult to transfer lessons from one disaster event to another. The aim of this research is to develop a systems understanding of disaster waste management and in turn develop context- and disaster-transferrable decision-making guidance for emergency and waste managers. To research this complex and multi-disciplinary problem, a multi-hazard, multi-context, multi-case study approach was adopted. The research focussed on five major disaster events: 2011 Christchurch earthquake, 2009 Victorian Bushfires, 2009 Samoan tsunami, 2009 L’Aquila earthquake and 2005 Hurricane Katrina. The first stage of the analysis involved the development of a set of ‘disaster & disaster waste’ impact indicators. The indicators demonstrate a method by which disaster managers, planners and researchers can simplify the very large spectra of possible disaster impacts, into some key decision-drivers which will likely influence post-disaster management requirements. The second stage of the research was to develop a set of criteria to represent the desirable environmental, economic, social and recovery effects of a successful disaster waste management system. These criteria were used to assess the effectiveness of the disaster waste management approaches for the case studies. The third stage of the research was the cross-case analysis. Six main elements of disaster waste management systems were identified and analysed. These were: strategic management, funding mechanisms, operational management, environmental and human health risk management, and legislation and regulation. Within each of these system elements, key decision-making guidance (linked to the ‘disaster & disaster waste’ indicators) and management principles were developed. The ‘disaster & disaster waste’ impact indicators, the effects assessment criteria and management principles have all been developed so that they can be practically applied to disaster waste management planning and response in the future.
High-quality ground motion records are required for engineering applications including response history analysis, seismic hazard development, and validation of physics-based ground motion simulations. However, the determination of whether a ground motion record is high-quality is poorly handled by automation with mathematical functions and can become prohibitive if done manually. Machine learning applications are well-suited to this problem, and a previous feed-forward neural network was developed (Bellagamba et al. 2019) to determine high-quality records from small crustal events in the Canterbury and Wellington regions for simulation validation. This prior work was however limited by the omission of moderate-to-large magnitude events and those from other tectonic environments, as well as a lack of explicit determination of the minimum usable frequency of the ground motion. To address these shortcomings, an updated neural network was developed to predict the quality of ground motion records for all magnitudes and all tectonic sources—active shallow crustal, subduction intraslab, and subduction interface—in New Zealand. The predictive performance of the previous feed-forward neural network was matched by the neural network in the domain of small crustal records, and this level of predictive performance is now extended to all source magnitudes and types in New Zealand making the neural network applicable to global ground motion databases. Furthermore, the neural network provides quality and minimum usable frequency predictions for each of the three orthogonal components of a record which may then be mapped into a binary quality decision or otherwise applied as desired. This framework provides flexibility for the end user to predict high-quality records with various acceptability thresholds allowing for this neural network to be used in a range of applications.
The development of Digital City technologies to manage and visualise spatial information has increasingly become a focus of the research community, and application by city authorities. Traditionally, the Geographic Information Systems (GIS) and Building Information Models (BIM) underlying Digital Cities have been used independently. However, integrating GIS and BIM into a single platform provides benefits for project and asset management, and is applicable to a range of issues. One of these benefits is the means to access and analyse large datasets describing the built environment, in order to characterise urban risk from and resilience to natural hazards. The aim of this thesis is to further explore methodologies of integration in two distinct areas. The first, integration through connectivity of heterogeneous datasets where GIS spatial infrastructure data is merged with 3D BIM building data to create a digital twin. Secondly, integration through analysis whereby data from the digital twin are extracted and integrated with computational models. To achieve this, a workflow was developed to identify the required datasets of a digital twin, and develop a process of integrating those datasets through a combination of; semi-autonomous conversion, translation and extension of data; and semantic web and services-based processes. Through use of a designed schema, the data were streamed in a homogenous format in a web-based platform. To demonstrate the value of this workflow with respect to urban risk and resilience, the process was applied to the Taiora: Queen Elizabeth II recreation and sports centre in eastern Christchurch, New Zealand. After integration of as-built GIS and BIM datasets, targeted data extraction was implemented, with outputs tailored for analysis in an infrastructure serviceability loss model, which assessed potable water network performance in the 22nd February 2011 Christchurch Earthquake. Using the same earthquake conditions as the serviceability loss model, performance of infrastructure assets in service at the time of the 22nd February 2011 Christchurch Earthquake was compared to new assets rebuilt at the site, post-earthquake. Due to improved potable water infrastructure resilience resulting from installation of ductile piles, a decrease of 35.5% in the probability of service loss was estimated in the serviceability loss model. To complete the workflow, the results from the external analysis were uploaded to the web-based platform. One of the more significant outcomes from the workflow was the identification of a lack of mandated metadata standards for fittings/valves connecting a building to private laterals. Whilst visually the GIS and BIM data show the building and pipes as connected, the semantic data does not include this connectivity relationship. This has no material impact on the current serviceability loss model as it is not one of the defined parameters. However, a proposed modification to the model would utilise the metadata to further assess the physical connection robustness, and increase the number of variables for estimating probability of service loss. This thesis has made a methodological contribution to urban resilience analysis by demonstrating how readily available up-to-date building and infrastructure data can be integrated, and with tailored extraction from a Digital City platform, be used for disaster impact analysis in an external computational engine, with results in turn imported and visualised in the Digital City platform. The workflow demonstrated that translation and integration of data would be more successful if a regional/national mandate was implemented for the submission of consent documentation in a specified standard BIM format. The results of this thesis have identified that the key to ensuring the success of an integrated tool lies in the initial workflow required to safeguard that all data can be either captured or translated in an interoperable format.
The development of Digital City technologies to manage and visualise spatial information has increasingly become a focus of the research community, and application by city authorities. Traditionally, the Geographic Information Systems (GIS) and Building Information Models (BIM) underlying Digital Cities have been used independently. However, integrating GIS and BIM into a single platform provides benefits for project and asset management, and is applicable to a range of issues. One of these benefits is the means to access and analyse large datasets describing the built environment, in order to characterise urban risk from and resilience to natural hazards. The aim of this thesis is to further explore methodologies of integration in two distinct areas. The first, integration through connectivity of heterogeneous datasets where GIS spatial infrastructure data is merged with 3D BIM building data to create a digital twin. Secondly, integration through analysis whereby data from the digital twin are extracted and integrated with computational models. To achieve this, a workflow was developed to identify the required datasets of a digital twin, and develop a process of integrating those datasets through a combination of; semi-autonomous conversion, translation and extension of data; and semantic web and services-based processes. Through use of a designed schema, the data were streamed in a homogenous format in a web-based platform. To demonstrate the value of this workflow with respect to urban risk and resilience, the process was applied to the Taiora: Queen Elizabeth II recreation and sports centre in eastern Christchurch, New Zealand. After integration of as-built GIS and BIM datasets, targeted data extraction was implemented, with outputs tailored for analysis in an infrastructure serviceability loss model, which assessed potable water network performance in the 22nd February 2011 Christchurch Earthquake. Using the same earthquake conditions as the serviceability loss model, performance of infrastructure assets in service at the time of the 22nd February 2011 Christchurch Earthquake was compared to new assets rebuilt at the site, post-earthquake. Due to improved potable water infrastructure resilience resulting from installation of ductile piles, a decrease of 35.5% in the probability of service loss was estimated in the serviceability loss model. To complete the workflow, the results from the external analysis were uploaded to the web-based platform. One of the more significant outcomes from the workflow was the identification of a lack of mandated metadata standards for fittings/valves connecting a building to private laterals. Whilst visually the GIS and BIM data show the building and pipes as connected, the semantic data does not include this connectivity relationship. This has no material impact on the current serviceability loss model as it is not one of the defined parameters. However, a proposed modification to the model would utilise the metadata to further assess the physical connection robustness, and increase the number of variables for estimating probability of service loss. This thesis has made a methodological contribution to urban resilience analysis by demonstrating how readily available up-to-date building and infrastructure data can be integrated, and with tailored extraction from a Digital City platform, be used for disaster impact analysis in an external computational engine, with results in turn imported and visualised in the Digital City platform. The workflow demonstrated that translation and integration of data would be more successful if a regional/national mandate was implemented for the submission of consent documentation in a specified standard BIM format. The results of this thesis have identified that the key to ensuring the success of an integrated tool lies in the initial workflow required to safeguard that all data can be either captured or translated in an interoperable format.
Geologic phenomena produced by earthquake shaking, including rockfalls and liquefaction features, provide important information on the intensity and spatiotemporal distribution of earthquake ground motions. The study of rockfall and liquefaction features produced in contemporary well- instrumented earthquakes increases our knowledge of how natural and anthropogenic environments respond to earthquakes and improves our ability to deduce seismologic information from analogous pre-contemporary (paleo-) geologic features. The study of contemporary and paleo- rockfall and liquefaction features enables improved forecasting of environmental responses to future earthquakes. In this thesis I utilize a combination of field and imagery-based mapping, trenching, stratigraphy, and numerical dating techniques to understand the nature and timing of rockfalls (and hillslope sedimentation) and liquefaction in the eastern South Island of New Zealand, and to examine the influence that anthropogenic activity has had on the geologic expressions of earthquake phenomena. At Rapaki (Banks Peninsula, NZ), field and imagery-based mapping, statistical analysis and numerical modeling was conducted on rockfall boulders triggered by the fatal 2011 Christchurch earthquakes (n=285) and compared with newly identified prehistoric (Holocene and Pleistocene) boulders (n=1049) deposited on the same hillslope. A significant population of modern boulders (n=26) travelled farther downslope (>150 m) than their most-travelled prehistoric counterparts, causing extensive damage to residential dwellings at the foot of the hillslope. Replication of prehistoric boulder distributions using 3-dimensional rigid body numerical models requires the application of a drag-coefficient, attributed to moderate to dense slope vegetation, to account for their spatial distribution. Radiocarbon dating provides evidence for 17th to early 20th century deforestation at the study site during Polynesian and European colonization and after emplacement of prehistoric rockfalls. Anthropocene deforestation enabled modern rockfalls to exceed the limits of their prehistoric predecessors, highlighting a shift in the geologic expression of rockfalls due to anthropogenic activity. Optical and radiocarbon dating of loessic hillslope sediments in New Zealand’s South Island is used to constrain the timing of prehistoric rockfalls and associated seismic events, and quantify spatial and temporal patterns of hillslope sedimentation including responses to seismic and anthropogenic forcing. Luminescence ages from loessic sediments constrain timing of boulder emplacement to between ~3.0 and ~12.5 ka, well before the arrival of Polynesians (ca AD 1280) and Europeans (ca AD 1800) in New Zealand, and suggest loess accumulation was continuing at the study site until 12-13 ka. Large (>5 m3) prehistoric rockfall boulders preserve an important record of Holocene hillslope sedimentation by creating local traps for sediment aggradation and upbuilding soil formation. Sediment accumulation rates increased considerably (>~10 factor increase) following human arrival and associated anthropogenic burning of hillslope vegetation. New numerical ages are presented to place the evolution of loess-mantled hillslopes in New Zealand’s South Island into a longer temporal framework and highlight the roles of earthquakes and humans on hillslope surface process. Extensive field mapping and characterization for 1733 individual prehistoric rockfall boulders was conducted at Rapaki and another Banks Peninsula site, Purau, to understand their origin, frequency, and spatial and volumetric distributions. Boulder characteristics and distributions were compared to 421 boulders deposited at the same sites during the 2010-2011 Canterbury earthquake sequence. Prehistoric boulders at Rapaki and Purau are comprised of two dominant lithofacies types: volcanic breccia and massive (coherent) lava basalt. Volcanic breccia boulders are found in greatest abundance (64-73% of total mapped rockfall) and volume (~90-96% of total rockfall) at both locations and exclusively comprise the largest boulders with the longest runout distances that pose the greatest hazard to life and property. This study highlights the primary influence that volcanic lithofacies architecture has on rockfall hazard. The influence of anthropogenic modifications on the surface and subsurface geologic expression of contemporary liquefaction created during the 2010-2011 Canterbury earthquake sequence (CES) in eastern Christchurch is examined. Trench observations indicate that anthropogenic fill layer boundaries and the composition/texture of discretely placed fill layers play an important role in absorbing fluidized sand/silt and controlling the subsurface architecture of preserved liquefaction features. Surface liquefaction morphologies (i.e. sand blows and linear sand blow arrays) display alignment with existing utility lines and utility excavations (and perforated pipes) provided conduits for liquefaction ejecta during the CES. No evidence of pre-CES liquefaction was identified within the anthropogenic fill layers or underlying native sediment. Radiocarbon dating of charcoal within the youngest native sediment suggests liquefaction has not occurred at the study site for at least the past 750-800 years. The importance of systematically examining the impact of buried infrastructure on channelizing and influencing surface and subsurface liquefaction morphologies is demonstrated. This thesis highlights the importance of using a multi-technique approach for understanding prehistoric and contemporary earthquake phenomena and emphasizes the critical role that humans play in shaping the geologic record and Earth’s surface processes.
Measurement of basement seismic resonance frequencies can elucidate shallow velocity structure, an important factor in earthquake hazard estimation. Ambient noise cross correlation, which is well-suited to studying shallow earth structure, is commonly used to analyze fundamental-mode Rayleigh waves and, increasingly, Love waves. Here we show via multicomponent ambient noise cross correlation that the basement resonance frequency in the Canterbury region of New Zealand can be straightforwardly determined based on the horizontal to vertical amplitude ratio (H/V ratio) of the first higher-mode Rayleigh waves. At periods of 1-3 s, the first higher-mode is evident on the radial-radial cross-correlation functions but almost absent in the vertical-vertical cross-correlation functions, implying longitudinal motion and a high H/V ratio. A one-dimensional regional velocity model incorporating a ~ 1.5 km-thick sedimentary layer fits both the observed H/V ratio and Rayleigh wave group velocity. Similar analysis may enable resonance characteristics of other sedimentary basins to be determined. © 2013. American Geophysical Union. All Rights Reserved.
This dissertation addresses a diverse range of applied aspects in ground motion simulation validation via the response of complex structures. In particular, the following topics are addressed: (i) the investigation of similarity between recorded and simulated ground motions using code-based 3D irregular structural response analysis, (ii) the development of a framework for ground motion simulations validation to identify the cause of differences between paired observed and simulated dataset, and (iii) the illustration of the process of using simulations for seismic performance-based assessment. The application of simulated ground motions is evaluated for utilisation in engineering practice by considering responses of 3D irregular structures. Validation is performed in a code-based context when the NZS1170.5 (NZS1170.5:2004, 2004) provisions are followed for response history analysis. Two real buildings designed by engineers and physically constructed in Christchurch before the 2010-2011 Canterbury earthquake sequence are considered. The responses are compared when the buildings are subjected to 40 scaled recorded and their subsequent simulated ground motions selected from 22 February 2011 Christchurch. The similarity of recorded and simulated responses is examined using statistical methods such as bootstrapping and hypothesis testing to determine whether the differences are statistically significant. The findings demonstrate the applicability of simulated ground motion when the code-based approach is followed in response history analysis. A conceptual framework is developed to link the differences between the structural response subjected to simulated and recorded ground motions to the differences in their corresponding intensity measures. This framework allows the variability to be partitioned into the proportion that can be “explained” by the differences in ground motion intensity measures and the remaining “unexplained” variability that can be attributed to different complexities such as dynamic phasing of multi-mode response, nonlinearity, and torsion. The application of this framework is examined through a hierarchy of structures reflecting a range of complexity from single-degree-of-freedom to 3D multi-degree-of-freedom systems with different materials, dynamic properties, and structural systems. The study results suggest the areas that ground motion simulation should focus on to improve simulations by prioritising the ground motion intensity measures that most clearly account for the discrepancies in simple to complex structural responses. Three approaches are presented to consider recorded or simulated ground motions within the seismic performance-based assessment framework. Considering the applications of ground motions in hazard and response history analyses, different pathways in utilising ground motions in both areas are explored. Recorded ground motions are drawn from a global database (i.e., NGA-West2 Ancheta et al., 2014). The NZ CyberShake dataset is used to obtain simulations. Advanced ground motion selection techniques (i.e., generalized conditional intensity measure, GCIM) are used for ground motion selection at a few intensity levels. The comparison is performed by investigating the response of an example structure (i.e., 12-storey reinforced concrete special moment frame) located in South Island, NZ. Results are compared and contrasted in terms of hazard, groundmotion selection, structural responses, demand hazard, and collapse risk, then, the probable reasons for differences are discussed. The findings from this study highlight the present opportunities and shortcomings in using simulations in risk assessment. i
The Porter's Pass-Amberley Fault Zone (PPAFZ) is a complex zone of anastomosing faults and folds bounding the south-eastern edge of the transition from subducting Pacific Plate to continental collision on the Australia Plate boundary. This study combines mapping of a 2000 km2 zone from the Southern Alps northeast to the coast near Amberley, 40 km north of metropolitan Christchurch, with an analysis of seismicity and a revision of regional seismic hazard. Three structural styles: 1) a western strike-slip, and 2) a more easterly thrust and reverse domain, pass into 3) a northwest verging fold belt on the northern Canterbury Plains, reflecting the structural levels exposed and the evolving west to east propagation. Basal remnants of a Late Cretaceous-Cenozoic, largely marine sedimentary cover sequence are preserved as outliers that unconformably overlie Mesozoic basement (greywacke and argillite of the Torlesse terrain) in the mountains of the PPAFZ and are underlain by a deeply leached zone which is widely preserved. Structure contouring of the unconformity surface indicates maximum, differential uplift of c.2600 m in the southwest, decreasing to c.1200 m in the coastal fold belt to the northeast. Much lower rates (or reversal) of uplift are evident a few kilometres southeast of the PPAFZ range-front escarpment. The youngest elements of the cover sequence are basement-derived conglomerates of Plio-Pleistocene age preserved on the SE margin. The source is more distant than the intervening mountains of the PPAFZ, probably from the Southern Alps, to the west and northwest. The absence of another regional unconformity on Mesozoic basement, older than Pleistocene, indicates that this uplift is post-Pliocene. Late Pleistocene(<100 kyr) differential uplift rates of c.0.5-2.7 m/kyr from uplifted marine terraces at the east coast, and rates of 2.5-3.3 m/kyr for tectonically-induced river-down cutting further west, suggest that uplift commenced locally during the last 1 Ma, and possibly within the last 0.5 Ma, if average rates are assumed to be uniform over time. Analysis of seismicity, recorded during a 10 week regional survey of micro earthquakes in 1990, identified two seismic zones beneath North Canterbury: 1) a sub-horizontal zone of activity restricted to the upper crust (≤12 km); and 2) a seismic zone in the lower crust (below a ceiling of ≤17 km), that broadens vertically to the north and northwest to a depth of c.40 km, with a bottom edge which dips 10°N and 15°NW, respectively. No events were recorded at depths between 12 km and 17 km, which is interpreted as a relatively aseismic, mid-crustal ductile layer. Marked differences (up to 60°) in the trend of strain axes for events above and below the inferred ductile layer are observed only north of the PPAFZ. A fundamental, north-to-south increase in the Wave-length of major geological structures occurs across the PPAFZ, and is interpreted as evidence that the upper crust beneath the Canterbury Plains is coupled to the lower crust, whereas the upper crust further north is not. Most of the recorded micro earthquakes <12 km deep beneath the PPAFZ have strike-slip mechanisms. It is probable that faults splay upward into the thrusts and folds at the surface as an evolving transpression zone in response to deep shear in basement. There have been no historic surface ruptures of the PPAFZ, but the zone has been characterised historically by frequent small earthquakes. Paleoseismic data (dated landslides and surface ruptures) compiled in this study, indicate a return period of 1500-1900 years between the last two M>7-7.5 earthquakes, and 500-700 years have elapsed since the last. The magnitudes of these events are estimated at c.M7.5, which represents a probable maximum magnitude for the PPAFZ. There are insufficient data to determine whether or not the frequency of large earthquakes conforms to a recognised model of behaviour, but comparison of the paleoseismic data with the historic record of smaller earthquakes, suggests that the magnitudes of the largest earthquakes in this zone are not exponentially distributed. A seismicity model for the PPAFZ (Elder et al., 1991) is reviewed, and a b-value of 1.0 is found to be consistent with the newly acquired paleoseismic data. This b-value reduces the predicted frequency of large earthquakes (M≥7.0) in this zone by a factor of 3.5, while retaining a conservative margin that allows for temporal variations in the frequency of large events and the possibility that the geological database is incomplete, suggesting grounds for revising the hazard model for Christchurch.
Natural hazard disasters often have large area-wide impacts, which can cause adverse stress-related mental health outcomes in exposed populations. As a result, increased treatment-seeking may be observed, which puts a strain on the limited public health care resources particularly in the aftermath of a disaster. It is therefore important for public health care planners to know whom to target, but also where and when to initiate intervention programs that promote emotional wellbeing and prevent the development of mental disorders after catastrophic events. A large body of literature assesses factors that predict and mitigate disaster-related mental disorders at various time periods, but the spatial component has rarely been investigated in disaster mental health research. This thesis uses spatial and spatio-temporal analysis techniques to examine when and where higher and lower than expected mood and anxiety symptom treatments occurred in the severely affected Christchurch urban area (New Zealand) after the 2010/11 Canterbury earthquakes. High-risk groups are identified and a possible relationship between exposure to the earthquakes and their physical impacts and mood and anxiety symptom treatments is assessed. The main research aim is to test the hypothesis that more severely affected Christchurch residents were more likely to show mood and anxiety symptoms when seeking treatment than less affected ones, in essence, testing for a dose-response relationship. The data consisted of mood and anxiety symptom treatment information from the New Zealand Ministry of Health’s administrative databases and demographic information from the National Health Index (NHI) register, when combined built a unique and rich source for identifying publically funded stress-related treatments for mood and anxiety symptoms in almost the whole population of the study area. The Christchurch urban area within the Christchurch City Council (CCC) boundary was the area of interest in which spatial variations in these treatments were assessed. Spatial and spatio-temporal analyses were done by applying retrospective space-time and spatial variation in temporal trends analysis using SaTScan™ software, and Bayesian hierarchical modelling techniques for disease mapping using WinBUGS software. The thesis identified an overall earthquake-exposure effect on mood and anxiety symptom treatments among Christchurch residents in the context of the earthquakes as they experienced stronger increases in the risk of being treated especially shortly after the catastrophic 2011 Christchurch earthquake compared to the rest of New Zealand. High-risk groups included females, elderly, children and those with a pre-existing mental illness with elderly and children especially at-risk in the context of the earthquakes. Looking at the spatio-temporal distribution of mood and anxiety symptom treatments in the Christchurch urban area, a high rates cluster ranging from the severely affected central city to the southeast was found post-disaster. Analysing residential exposure to various earthquake impacts found that living in closer proximity to more affected areas was identified as a risk factor for mood and anxiety symptom treatments, which largely confirms a dose-response relationship between level of affectedness and mood and anxiety symptom treatments. However, little changes in the spatial distribution of mood and anxiety symptom treatments occurred in the Christchurch urban area over time indicating that these results may have been biased by pre-existing spatial disparities. Additionally, the post-disaster mobility activity from severely affected eastern to the generally less affected western and northern parts of the city seemed to have played an important role as the strongest increases in treatment rates occurred in less affected northern areas of the city, whereas the severely affected eastern areas tended to show the lowest increases. An investigation into the different effects of mobility confirmed that within-city movers and temporary relocatees were generally more likely to receive care or treatment for mood or anxiety symptoms, but moving within the city was identified as a protective factor over time. In contrast, moving out of the city from minor, moderately or severely damaged plain areas of the city, which are generally less affluent than Port Hills areas, was identified as a risk factor in the second year post-disaster. Moreover, residents from less damaged plain areas of the city showed a decrease in the likelihood of receiving care or treatment for mood or anxiety symptoms compared to those from undamaged plain areas over time, which also contradicts a possible dose-response relationship. Finally, the effects of the social and physical environment, as well as community resilience on mood and anxiety symptom treatments among long-term stayers from Christchurch communities indicate an exacerbation of pre-existing mood and anxiety symptom treatment disparities in the city, whereas exposure to ‘felt’ earthquake intensities did not show a statistically significant effect. The findings of this thesis highlight the complex relationship between different levels of exposure to a severe natural disaster and adverse mental health outcomes in a severely affected region. It is one of the few studies that have access to area-wide health and impact information, are able to do a pre-disaster / post-disaster comparison and track their sample population to apply spatial and spatio-temporal analysis techniques for exposure assessment. Thus, this thesis enhances knowledge about the spatio-temporal distribution of adverse mental health outcomes in the context of a severe natural disaster and informs public health care planners, not only about high-risk groups, but also where and when to target health interventions. The results indicate that such programs should broadly target residents living in more affected areas as they are likely to face daily hardship by living in a disrupted environment and may have already been the most vulnerable ones before the disaster. Special attention should be focussed on women, elderly, children and people with pre-existing mental illnesses as they are most likely to receive care or treatment for stress-related mental health symptoms. Moreover, permanent relocatees from affected areas and temporarily relocatees shortly after the disaster may need special attention as they face additional stressors due to the relocation that may lead to the development of adverse mental health outcomes needing treatment.
A zone of active tectonism occurs in mid and north Canterbury, from the Rakaia to the Waipara Rivers, which coincides with seismicity concentrations and several Quaternary surface anomalies and is here defined as the Porters Pass Tectonic Zone. Although parallel to the Marlborough faults to the north, the lack of regional definition suggests this zone is much younger in its inception reflecting a southward movement of the plate rotation vector. The objectives of this study were to map the structures associated with this zone in the segment between the Rakaia and Waimakariri Rivers with detailed analysis concentrated in the upper Kawai Valley. Quaternary offsets on the main lineament of the Porters Pass Fault were traced through the area and evidence for the rate of movement, probable magnitudes and return periods of related seismic events was sought. The basement was found to be complicated by pre-existing deformation structures in Torlesse Group rocks which have been subsequently been re-activated or rotated by recent fault movement probably beginning in the Pleistocene. This phase is dominantly thrusting and uplift has lead to the erosion of most of the overlying sedimentary cover. Remnants of the Cret-Tertiary sediments still remain as fault-bounded packets. Evidence suggests that a change to development of a regional lateral shear associated with the Porters Pass Tectonic Zone transects the thrust system with complex interaction between the older reverse and new strike-slip faults. Offset rates along the segments of the Porters Pass Fault are not well constrained but are believed to be approximately in the range of 11-13 mm/year for at least the last 130,000 years. This rate is similar to other large faults in the Marlborough region. Two earthquake events have been identified and dated at 600 and 2000 years ago, with a magnitude of greater than 6.5. Evidence suggests characteristic earthquakes along the Porters Pass Fault are greater than Magnitude 7. This result has some major ramifications for the expected seismic hazards for nearby Christchurch.
Though generally considered “natural” disasters, cyclones and earthquakes are increasingly being associated with human activities, incubated through urban settlement patterns and the long-term redistribution of natural resources. As society is becoming more urbanized, the risk of human exposure to disasters is also rising. Architecture often reflects the state of society’s health: architectural damage is the first visible sign of emergency, and reconstruction is the final response in the process of recovery. An empirical assessment of architectural projects in post-disaster situations can lead to a deeper understanding of urban societies as they try to rebuild. This thesis offers an alternative perspective on urban disasters by looking at the actions and attitudes of disaster professionals through the lens of architecture, situated in recent events: the 2010 Christchurch earthquake, the 2010 Haiti earthquake, and the 2005 Hurricane Katrina. An empirical, multi-hazard, cross-sectional case study methodology was used, employing grounded theory method to build theory, and a critical constructivist strategy to inform the analysis. By taking an interdisciplinary approach to understanding disasters, this thesis positions architecture as a conduit between two divergent approaches to disaster research: the hazards approach, which studies the disaster cycles from a scientific perspective; and the sociological approach, which studies the socially constructed vulnerabilities that result from disasters, and the elements of social change that accompany such events. Few studies to date have attempted to integrate the multi-disciplinary perspectives that can advance our understanding of societal problems in urban disasters. To bridge this gap, this thesis develops what will be referred to as the “Rittelian framework”—based on the work of UC Berkeley’s architecture professor Horst Rittel (1930-1990). The Rittelian framework uses the language of design to transcend the multiple fields of human endeavor to address the “design problems” in disaster research. The processes by which societal problems are addressed following an urban disaster involve input by professionals from multiple fields—including economics, sociology, medicine, and engineering—but the contribution from architecture has been minimal to date. The main impetus for my doctoral thesis has been the assertion that most of the decisions related to reconstruction are made in the early emergency recovery stages where architects are not involved, but architects’ early contribution is vital to the long-term reconstruction of cities. This precipitated in the critical question: “How does the Rittelian framework contribute to the critical design decisions in modern urban disasters?” Comparative research was undertaken in three case studies of recent disasters in New Orleans (2005), Haiti (2010) and Christchurch (2010), by interviewing 51 individuals who were selected on the basis of employing the Rittelian framework in their humanitarian practice. Contextualizing natural disaster research within the robust methodological framework of architecture and the analytical processes of sociology is the basis for evaluating the research proposition that architectural problem solving is of value in addressing the ‘Wicked Problems’ of disasters. This thesis has found that (1) the nuances of the way disaster agents interpret the notion of “building back better” can influence the extent to which architectural professionals contribute in urban disaster recovery, (2) architectural design can be used to facilitate but also impede critical design decisions, and (3) framing disaster research in terms of design decisions can lead to innovation where least expected. This empirical research demonstrates how the Rittelian framework can inform a wider discussion about post-disaster human settlements, and improve our resilience through disaster research.
