Search

found 39 results

Research papers, University of Canterbury Library

The aim of this thesis was to examine the spatial and the temporal patterns of anxiety and chest pain resulting from the Canterbury, New Zealand earthquaeks. Three research objectives were identified: examine any spatial or termporal clusters of anxiety and chest pain; examine the associations between anxiety, chest pain and damage to neighbourhood; and determine any statistically significant difference in counts of anxiety and chest pain after each earthquake or aftershock which resulted in severe damage. Measures of the extent of liquefaction the location of CERA red-zones were used as proxy measures for earthquake damage. Cases of those who presented to Christchurch Public Hospital Emergency Department with either anxiety or chest pain between May 2010 and April 2012 were aggregated to census area unit (CAU) level for analysis. This thesis has taken a unique approach to examining the spatial and spatio-temporal variations of anxiety and chest pain after an earthquake and offers unique results. This is the first study of its kind to use a GIS approach when examining Canterbury specific earthquake damage and health variables at a CAU level after the earthquakes. Through the use of spatio-termporal scan modelling, negative and linear regression modelling and temporal linear modelling with dummy variables this research was able to conclude there are significant spatial and temporal variations in anxiety and chest pain resulting from the earthquakes. The spatio-termporal scan modelling identified a hot cluster of both anxiety and chest pain within Christchurch at the same time the earthquakes occurred. The negative binomial model found liquefaction to be a stronger predictor of anxiety than the Canterbury Earthquake Recovery Authority's (CERA) land zones. The linear regression model foun chest pain to be positively associated with all measures of earthquake damage with the exception of being in the red-zone. The temporal modelling identified a significant increase in anxiety cases one month after a major earthquake, and chest pain cases spiked two weeks after an earthquake and gradually decreased over the following five weeks. This research was limited by lack of control period data, limited measures of earthquake damage, ethical restrictions, and the need for population tracking data. The findings of this research will be useful in the planning and allocation of mental wellbeing resources should another similar event like the Canterbury Earthquakes occur in New Zealand.

Research papers, Lincoln University

Globally, the maximum elevations at which treelines are observed to occur coincide with a 6.4 °C soil isotherm. However, when observed at finer scales, treelines display a considerable degree of spatial complexity in their patterns across the landscape and are often found occurring at lower elevations than expected relative to the global-scale pattern. There is still a lack of understanding of how the abiotic environment imposes constraints on treeline patterns, the scales at which different effects are acting, and how these effects vary over large spatial extents. In this thesis, I examined abrupt Nothofagus treelines across seven degrees of latitude in New Zealand in order to investigate two broad questions: (1) What is the nature and extent of spatial variability in Nothofagus treelines across the country? (2) How is this variation associated with abiotic variation at different spatial scales? A range of GIS, statistical, and atmospheric modelling methods were applied to address these two questions. First, I characterised Nothofagus treeline patterns at a 15x15km scale across New Zealand using a set of seven, GIS-derived, quantitative metrics that describe different aspects of treeline position, shape, spatial configuration, and relationships with adjacent vegetation. Multivariate clustering of these metrics revealed distinct treeline types that showed strong spatial aggregation across the country. This suggests a strong spatial structuring of the abiotic environment which, in turn, drives treeline patterns. About half of the multivariate treeline metric variation was explained by patterns of climate, substrate, topographic and disturbance variability; on the whole, climatic and disturbance factors were most influential. Second, I developed a conceptual model that describes how treeline elevation may vary at different scales according to three categories of effects: thermal modifying effects, physiological stressors, and disturbance effects. I tested the relevance of this model for Nothofagus treelines by investigating treeline elevation variation at five nested scales (regional to local) using a hierarchical design based on nested river catchments. Hierarchical linear modelling revealed that the majority of the variation in treeline elevation resided at the broadest, regional scale, which was best explained by the thermal modifying effects of solar radiation, mountain mass, and differences in the potential for cold air ponding. Nonetheless, at finer scales, physiological and disturbance effects were important and acted to modify the regional trend at these scales. These results suggest that variation in abrupt treeline elevations are due to both broad-scale temperature-based growth limitation processes and finer-scale stress- and disturbance-related effects on seedling establishment. Third, I explored the applicability of a meso-scale atmospheric model, The Air Pollution Model (TAPM), for generating 200 m resolution, hourly topoclimatic data for temperature, incoming and outgoing radiation, relative humidity, and wind speeds. Initial assessments of TAPM outputs against data from two climate station locations over seven years showed that the model could generate predictions with a consistent level of accuracy for both sites, and which agreed with other evaluations in the literature. TAPM was then used to generate data at 28, 7x7 km Nothofagus treeline zones across New Zealand for January (summer) and July (winter) 2002. Using mixed-effects linear models, I determined that both site-level factors (mean growing season temperature, mountain mass, precipitation, earthquake intensity) and local-level landform (slope and convexity) and topoclimatic factors (solar radiation, photoinhibition index, frost index, desiccation index) were influential in explaining variation in treeline elevation within and among these sites. Treelines were generally closer to their site-level maxima in regions with higher mean growing season temperatures, larger mountains, and lower levels of precipitation. Within sites, higher treelines were associated with higher solar radiation, and lower photoinhibition and desiccation index values, in January, and lower desiccation index values in July. Higher treelines were also significantly associated with steeper, more convex landforms. Overall, this thesis shows that investigating treelines across extensive areas at multiple study scales enables the development of a more comprehensive understanding of treeline variability and underlying environmental constraints. These results can be used to formulate new hypotheses regarding the mechanisms driving treeline formation and to guide the optimal choice of field sites at which to test these hypotheses.

Research papers, University of Canterbury Library

Liquefaction-induced lateral spreading during earthquakes poses a significant hazard to the built environment, as observed in Christchurch during the 2010 to 2011 Canterbury Earthquake Sequence (CES). It is critical that geotechnical earthquake engineers are able to adequately predict both the spatial extent of lateral spreads and magnitudes of associated ground movements for design purposes. Published empirical and semi-empirical models for predicting lateral spread displacements have been shown to vary by a factor of <0.5 to >2 from those measured in parts of Christchurch during CES. Comprehensive post- CES lateral spreading studies have clearly indicated that the spatial distribution of the horizontal displacements and extent of lateral spreading along the Avon River in eastern Christchurch were strongly influenced by geologic, stratigraphic and topographic features.

Research papers, University of Canterbury Library

This dissertation addresses a diverse range of topics in the physics-based broadband ground motion simulation, with a focus on New Zealand applications. In particular the following topics are addressed: the methodology and computational implementation of a New Zealand Velocity Model for broadband ground motion simulation; generalised parametric functions and spatial correlations for seismic velocities in the Canterbury, New Zealand region from surface-wave-based site characterisation; and ground motion simulations of Hope Fault earthquakes. The paragraphs below outline each contribution in more detail. A necessary component in physics-based ground motion simulation is a 3D model which details the seismic velocities in the region of interest. Here a velocity model construction methodology, its computational implementation, and application in the construction of a New Zealand velocity model for use in physics-based broadband ground motion simulation are presented. The methodology utilises multiple datasets spanning different length scales, which is enabled via the use of modular sub-regions, geologic surfaces, and parametric representations of crustal velocity. A number of efficiency-related workflows to decrease the overall computational construction time are employed, while maintaining the flexibility and extensibility to incorporate additional datasets and re- fined velocity parameterizations as they become available. The model comprises explicit representations of the Canterbury, Wellington, Nelson-Tasman, Kaikoura, Marlborough, Waiau, Hanmer and Cheviot sedimentary basins embedded within a regional travel-time tomography-based velocity model for the shallow crust and provides the means to conduct ground motion simulations throughout New Zealand for the first time. Recently developed deep shear-wave velocity profiles in Canterbury enabled models that better characterise the velocity structure within geologic layers of the Canterbury sedimentary basin to be developed. Here the development of depth- and Vs30-dependent para-metric velocity and spatial correlation models to characterise shear-wave velocities within the geologic layers of the Canterbury sedimentary basin are presented. The models utilise data from 22 shear-wave velocity profiles of up to 2.5km depth (derived from surface wave analysis) juxtaposed with models which detail the three-dimensional structure of the geologic formations in the Canterbury sedimentary basin. Parametric velocity equations are presented for Fine Grained Sediments, Gravels, and Tertiary layer groupings. Spatial correlations were developed and applied to generate three-dimensional stochastic velocity perturbations. Collectively, these models enable seismic velocities to be realistically represented for applications such as 3D ground motion and site response simulations. Lastly the New Zealand velocity model is applied to simulate ground motions for a Mw7.51 rupture of the Hope Fault using a physics-based simulation methodology and a 3D crustal velocity model of New Zealand. The simulation methodology was validated for use in the region through comparison with observations for a suite of historic small magnitude earthquakes located proximal to the Hope Fault. Simulations are compared with conventionally utilised empirical ground motion models, with simulated peak ground velocities being notably higher in regions with modelled sedimentary basins. A sensitivity analysis was undertaken where the source characteristics of magnitude, stress parameter, hypocentre location and kinematic slip distribution were varied and an analysis of their effect on ground motion intensities is presented. It was found that the magnitude and stress parameter strongly influenced long and short period ground motion amplitudes, respectively. Ground motion intensities for the Hope Fault scenario are compared with the 2016 Kaikoura Mw7.8 earthquake, it was found that the Kaikoura earthquake produced stronger motions along the eastern South Island, while the Hope Fault scenario resulted in stronger motions immediately West of the near-fault region. The simulated ground motions for this scenario complement prior empirically-based estimates and are informative for mitigation and emergency planning purposes.

Research papers, University of Canterbury Library

This paper summarizes the development of a high-resolution surficial shear wave velocity model based on the combination of the large high-spatial-density database of cone penetration test (CPT) logs in and around Christchurch, New Zealand and a recently-developed Christchurch-specific empirical correlation between soil shear wave velocity and CPT. This near-surface shear wave velocity model has applications for site characterization efforts via the development of maps of time-averaged shear wave velocities over specific depths, as well as use in site response analysis and ground motion simulation.

Research papers, University of Canterbury Library

The overarching goal of this dissertation is to improve predictive capabilities of geotechnical seismic site response analyses by incorporating additional salient physical phenomena that influence site effects. Specifically, multidimensional wave-propagation effects that are neglected in conventional 1D site response analyses are incorporated by: (1) combining results of 3D regional-scale simulations with 1D nonlinear wave-propagation site response analysis, and (2) modelling soil heterogeneity in 2D site response analyses using spatially-correlated random fields to perturb soil properties. A method to combine results from 3D hybrid physics-based ground motion simulations with site-specific nonlinear site response analyses was developed. The 3D simulations capture 3D ground motion phenomena on a regional scale, while the 1D nonlinear site response, which is informed by detailed site-specific soil characterization data, can capture site effects more rigorously. Simulations of 11 moderate-to-large earthquakes from the 2010-2011 Canterbury Earthquake Sequence (CES) at 20 strong motion stations (SMS) were used to validate simulations with observed ground motions. The predictions were compared to those from an empirically-based ground motion model (GMM), and from 3D simulations with simplified VS30- based site effects modelling. By comparing all predictions to observations at seismic recording stations, it was found that the 3D physics-based simulations can predict ground motions with comparable bias and uncertainty as the GMM, albeit, with significantly lower bias at long periods. Additionally, the explicit modelling of nonlinear site-response improves predictions significantly compared to the simplified VS30-based approach for soft-soil or atypical sites that exhibit exceptionally strong site effects. A method to account for the spatial variability of soils and wave scattering in 2D site response analyses was developed and validated against a database of vertical array sites in California. The inputs required to run the 2D analyses are nominally the same as those required for 1D analyses (except for spatial correlation parameters), enabling easier adoption in practice. The first step was to create the platform and workflow, and to perform a sensitivity study involving 5,400 2D model realizations to investigate the influence of random field input parameters on wave scattering and site response. Boundary conditions were carefully assessed to understand their effect on the modelled response and select appropriate assumptions for use on a 2D model with lateral heterogeneities. Multiple ground-motion intensity measures (IMs) were analyzed to quantify the influence from random field input parameters and boundary conditions. It was found that this method is capable of scattering seismic waves and creating spatially-varying ground motions at the ground surface. The redistribution of ground-motion energy across wider frequency bands, and the scattering attenuation of high-frequency waves in 2D analyses, resemble features observed in empirical transfer functions (ETFs) computed in other studies. The developed 2D method was subsequently extended to more complicated multi-layer soil profiles and applied to a database of 21 vertical array sites in California to test its appropriate- ness for future predictions. Again, different boundary condition and input motion assumptions were explored to extend the method to the in-situ conditions of a vertical array (with a sensor embedded in the soil). ETFs were compared to theoretical transfer functions (TTFs) from conventional 1D analyses and 2D analyses with heterogeneity. Residuals of transfer-function- based IMs, and IMs of surface ground motions, were also used as validation metrics. The spatial variability of transfer-function-based IMs was estimated from 2D models and compared to the event-to-event variability from ETFs. This method was found capable of significantly improving predictions of median ETF amplification factors, especially for sites that display higher event-to-event variability. For sites that are well represented by 1D methods, the 2D approach can underpredict amplification factors at higher modes, suggesting that the level of heterogeneity may be over-represented by the 2D random field models used in this study.

Research papers, University of Canterbury Library

This paper summarizes the development of a region-wide surficial shear wave velocity model based on the combination of the large high-spatial-density database of cone penetration test (CPT) logs in and around Christchurch, New Zealand and a recently-developed Christchurch-specific empirical correlation between soil shear wave velocity and CPT. The ongoing development of this near-surface shear wave velocity model has applications for site characterization efforts via the development of maps of time-averaged shear wave velocities over specific depths, and the identification of regional similarities and differences in soil shear stiffness.

Research papers, University of Canterbury Library

In this paper we apply Full waveform tomography (FWT) based on the Adjoint-Wavefield (AW) method to iteratively invert a 3-D geophysical velocity model for the Canterbury region (Lee, 2017) from a simple initial model. The seismic wavefields was generated using numerical solution of the 3-D elastodynamic/ visco- elastodynamic equations (EMOD3D was adopted (Graves, 1996)), and through the AW method, gradients of model parameters (compression and shear wave velocity) were computed by implementing the cross-adjoint of forward and backward wavefields. The reversed-in-time displacement residual was utilized as the adjoint source. For inversion, we also account for the near source/ station effects, gradient precondition, smoothening (Gaussian filter in spatial domain) and optimal step length. Simulation-to-observation misfit measurements based on 191 sources at 78 seismic stations in the Canterbury region (Figure 1) were used into our inversion. The inversion process includes multiple frequency bands, starting from 0-0.05Hz, and advancing to higher frequency bands (0-0.1Hz and 0-0.2Hz). Each frequency band was used for up to 10 iterations or no optimal step length found. After 3 FWT inversion runs, the simulated seismograms computed using our final model show a good matching with the observed seismograms at frequencies from 0 - 0.2 Hz and the normalized least-squared misfit error has been significantly reduced. Over all, the synthetic study of FWT shows a good application to improve the crustal velocity models from the existed geological models and the seismic data of the different earthquake events happened in the Canterbury region.

Research papers, University of Canterbury Library

The potential for a gastroenteritis outbreak in a post-earthquake environment may increase because of compromised infrastructure services, contaminated liquefaction (lateral spreading and surface ejecta), and the presence of gastroenteritis agents in the drinking water network. A population in a post-earthquake environment might be seriously affected by gastroenteritis because it has a short incubation period (about 10 hours). The potential for a gastroenteritis outbreak in a post-earthquake environment may increase because of compromised infrastructure services, contaminated liquefaction (lateral spreading and surface ejecta), and the presence of gastroenteritis agents in the drinking water network. A population in a post-earthquake environment might be seriously affected by gastroenteritis because it has a short incubation period (about 10 hours). The aim of this multidisciplinary research was to retrospectively analyse the gastroenteritis prevalence following the February 22, 2011 earthquake in Christchurch. The first focus was to assess whether earthquake-induced infrastructure damage, liquefaction, and gastroenteritis agents spatially explained the recorded gastroenteritis cases over the period of 35 days following the February 22, 2011 earthquake in Christchurch. The gastroenteritis agents considered in this study were Escherichia coli found in the drinking water supply (MPN/100mL) and Non-Compliant Free Associated Chlorine (FAC-NC) (less than <0.02mg/L). The second focus was the protocols that averted a gastroenteritis outbreak at three Emergency Centres (ECs): Burnside High School Emergency Centre (BEC); Cowles Stadium Emergency Centre (CEC); and Linwood High School Emergency Centre (LEC). Using a mixed-method approach, gastroenteritis point prevalence and the considered factors were quantitatively analysed. The qualitative analysis involved interviewing 30 EC staff members. The data was evaluated by adopting the Grounded Theory (GT) approach. Spatial analysis of considered factors showed that highly damaged CAUs were statistically clustered as demonstrated by Moran’s I statistic and hot spot analysis. Further modelling showed that gastroenteritis point prevalence clustering could not be fully explained by infrastructure damage alone, and other factors influenced the recorded gastroenteritis point prevalence. However, the results of this research suggest that there was a tenuous, indirect relationship between recorded gastroenteritis point prevalence and the considered factors: earthquake-induced infrastructure damage, liquefaction and FAC-NC. Two ECs were opened as part of the post-earthquake response in areas with severe infrastructure damage and liquefaction (BEC and CEC). The third EC (CEC) provided important lessons that were learnt from the previous September 4, 2010 earthquake, and implemented after the February 22, 2011 earthquake. Two types of interwoven themes identified: direct and indirect. The direct themes were preventive protocols and indirect themes included type of EC building (school or a sports stadium), and EC staff. The main limitations of the research were Modifiable Areal Units (MAUP), data detection, and memory loss. This research provides a practical method that can be adapted to assess gastroenteritis risk in a post-earthquake environment. Thus, this mixed method approach can be used in other disaster contexts to study gastroenteritis prevalence, and can serve as an appendage to the existing framework for assessing infectious diseases. Furthermore, the lessons learnt from qualitative analysis can inform the current infectious disease management plans, designed for a post-disaster response in New Zealand and internationally Using a mixed-method approach, gastroenteritis point prevalence and the considered factors were quantitatively analysed. A damage profile was created by amalgamating different types of damage for the considered factors for each Census Area Unit (CAU) in Christchurch. The damage profile enabled the application of a variety of statistical methods which included Moran’s I , Hot Spot (HS) analysis, Spearman’s Rho, and Besag–York–Mollié Model using a range of software. The qualitative analysis involved interviewing 30 EC staff members. The data was evaluated by adopting the Grounded Theory (GT) approach. Spatial analysis of considered factors showed that highly damaged CAUs were statistically clustered as demonstrated by Moran’s I statistic and hot spot analysis. Further modelling showed that gastroenteritis point prevalence clustering could not be fully explained by infrastructure damage alone, and other factors influenced the recorded gastroenteritis point prevalence. However, the results of this research suggest that there was a tenuous, indirect relationship between recorded gastroenteritis point prevalence and the considered factors: earthquake-induced infrastructure damage, liquefaction and FAC-NC. Two ECs were opened as part of the post-earthquake response in areas with severe infrastructure damage and liquefaction (BEC and CEC). The third EC (CEC) provided important lessons that were learnt from the previous September 4, 2010 earthquake, and implemented after the February 22, 2011 earthquake. The ECs were selected to represent the Christchurch area, and were situated where potential for gastroenteritis was high. BEC represented the western side of Christchurch; whilst, CEC and LEC represented the eastern side, where the potential for gastroenteritis was high according to the outputs of the quantitative spatial modelling. Qualitative analysis from the interviews at the ECs revealed that evacuees were arriving at the ECs with gastroenteritis-like symptoms. Participants believed that those symptoms did not originate at the ECs. Two types of interwoven themes identified: direct and indirect. The direct themes were preventive protocols that included prolific use of hand sanitisers; surveillance; and the services offered. Indirect themes included the EC layout, type of EC building (school or a sports stadium), and EC staff. Indirect themes governed the quality and sustainability of the direct themes implemented, which in turn averted gastroenteritis outbreaks at the ECs. The main limitations of the research were Modifiable Areal Units (MAUP), data detection, and memory loss. It was concluded that gastroenteritis point prevalence following the February 22, 2011 earthquake could not be solely explained by earthquake-induced infrastructure damage, liquefaction, and gastroenteritis causative agents alone. However, this research provides a practical method that can be adapted to assess gastroenteritis risk in a post-earthquake environment. Creating a damage profile for each CAU and using spatial data analysis can isolate vulnerable areas, and qualitative data analysis provides localised information. Thus, this mixed method approach can be used in other disaster contexts to study gastroenteritis prevalence, and can serve as an appendage to the existing framework for assessing infectious diseases. Furthermore, the lessons learnt from qualitative analysis can inform the current infectious disease management plans, designed for a post-disaster response in New Zealand and internationally.

Research papers, Lincoln University

Prognostic modelling provides an efficient means to analyse the coastal environment and provide effective knowledge for long term urban planning. This paper outlines how the use of SWAN and Xbeach numerical models within the ESRI ArcGIS interface can simulate geomorphological evolution through hydrodynamic forcing for the Greater Christchurch coastal environment. This research followed the data integration techniques of Silva and Taborda (2012) and utilises their beach morphological modelling tool (BeachMM tool). The statutory requirements outlined in the New Zealand Coastal Policy Statement 2010 were examined to determine whether these requirements are currently being complied with when applying the recent sea level rise predictions by the Intergovernmental Panel on Climate Change (2013), and it would appear that it does not meet those requirements. This is because coastal hazard risk has not been thoroughly quantified by the installation of the Canterbury Earthquake Recovery Authority (CERA) residential red zone. However, the Christchurch City Council’s (CCC) flood management area does provide an extent to which managed coastal retreat is a real option. This research assessed the effectiveness of the prognostic models, forecasted a coastline for 100 years from now, and simulated the physical effects of extreme events such as storm surge given these future predictions. The results of this research suggest that progradation will continue to occur along the Christchurch foreshore due to the net sediment flux retaining an onshore direction and the current hydrodynamic activity not being strong enough to move sediment offshore. However, inundation during periods of storm surge poses a risk to human habitation on low lying areas around the Avon-Heathcote Estuary and the Brooklands lagoon similar to the CCC’s flood management area. There are complex interactions at the Waimakariri River mouth with very high rates of accretion and erosion within a small spatial scale due to the river discharge. There is domination of the marine environment over the river system determined by the lack of generation of a distinct river delta, and river channel has not formed within the intertidal zone clearly. The Avon-Heathcote ebb tidal delta aggrades on the innner fan and erodes on the outer fan due to wave domination. The BeachMM tool facilitates the role of spatial and temporal analysis effectively and the efficiency of that performance is determined by the computational operating system.

Research papers, University of Canterbury Library

This report summarizes the development of a region-wide surficial soil shear wave velocity (Vs ) model based on the unique combination of a large high-spatial-density database of cone penetration test (CPT) logs in the greater Christchurch urban area (> 15, 000 logs as of 1 February 2014) and the Christchurch-specific empirical correlation between soil Vs and CPT data developed by McGann et al. [1, 2]. This model has applications for site characterization efforts via maps of time-averaged Vs over specific depths (e.g. Vs30, Vs10), and for numerical modeling efforts via the identification of typical Vs profiles for different regions and soil behaviour types within Christchurch. In addition, the Vs model can be used to constrain the near-surface velocities for the 3D seismic velocity model of the Canterbury basin [3] currently being developed for the purpose of broadband ground motion simulation. The general development of these region-wide near-surface Vs models includes the following general phases, with each discussed in separate chapters of this report. • An evaluation of the available CPT dataset for suitability, and the definition of other datasets and assumptions necessary to characterize the surficial sediments of the region to 30 m depth. • The development of time-averaged shear wave velocity (Vsz) surfaces for the Christchurch area from the adopted CPT dataset (and supplementary data/assumptions) using spatial interpolation. The Vsz surfaces are used to explore the characteristics of the near-surface soils in the regions and are shown to correspond well with known features of the local geology, the historical ecosystems of the area, and observations made following the 2010- 2011 Canterbury earthquakes. • A detailed analysis of the Vs profiles in eight subregions of Christchurch is performed to assess the variablity in the soil profiles for regions with similar Vsz values and to assess Vsz as a predictive metric for local site response. It is shown that the distrubution of soil shear wave velocity in the Christchurch regions is highly variable both spatially (horizontally) and with depth (vertically) due to the varied geological histories for different parts of the area, and the highly stratified nature of the nearsurface deposits. This variability is not considered to be greatly significant in terms of current simplified site classification systems; based on computed Vs30 values, all considered regions can be categorized as NEHRP sites class D (180 < Vs < 360 m/s) or E (Vs < 180 m/s), however, detailed analysis of the shear wave velocity profiles in different subregions of Christchurch show that the expected surficial site response can vary quite a bit across the region despite the relative similarity in Vs30

Research papers, University of Canterbury Library

This paper presents a methodology by which both site-specific and spatially distributed ground motion intensity can be obtained immediately following an earthquake event. The methodology makes use of both prediction models for ground motion intensity and its correlation over spatial distances. A key benefit of the methodology is that the ground motion intensity at a given location is not a single value but a distribution of values. The distribution is comprised of both a mean and also standard deviation, with the standard deviation being a function of the distance to nearby strong motion stations. The methodology is illustrated for two applications. Firstly, maps of conditional peak ground acceleration (PGA) have been developed for the major events in the Canterbury earthquake sequence. It is illustrated how these conditional maps can be used for post-event evaluation of liquefaction triggering criteria which have been adopted by the Department of Building and Housing (DBH). Secondly, the conditional distribution of response spectral ordinates is obtained at a specific location for the purposes of determining appropriate ground motion records for use in seismic response analyses of important structures at locations where direct recordings are absent.

Research papers, University of Canterbury Library

This paper presents site-specific and spatially-distributed ground-motion intensity estimates which have been utilized in the aftermath of the 2010-2011 Canterbury, New Zealand earthquakes. The methodology underpinning the ground motion intensity estimation makes use of both prediction models for ground motion intensity and its within-event spatial correlation. A key benefit of the methodology is that the estimated ground motion intensity at a given location is not a single value but a distribution of values. The distribution is comprised of both a mean and standard deviation, with the standard deviation being a function of the distance to nearby observations at strong motion stations. The methodology is illustrated for two applications. Firstly, maps of conditional peak ground acceleration (PGA) have been developed for the major events in the Canterbury earthquake sequence, which among other things, have been utilized for assessing liquefaction triggering susceptibility of land in residential areas. Secondly, the conditional distribution of response spectral ordinates is obtained at the location of the Canterbury Television building (CTV), which catastrophically collapsed in the 22 February 2011 earthquake. The conditional response spectra provide insight for the selection of ground motion records for use in forensic seismic response analyses of important structures at locations where direct recordings are absent.

Research papers, University of Canterbury Library

This paper presents a seismic velocity model of Canterbury, New Zealand based on 3D geologic surfaces and velocities from a range of data sources. The model provides the 3D crustal structure in the region at multiple length scales for seismic wave propagation simulations, such as broadband ground motion and shallow site response analyses related to understanding the ground motions and site responses during the 2010- 2011 Canterbury earthquakes. Pre-Quaternary geologic horizons are calculated based on the reinterpretation of a comprehensive network of seismic reflection surveys from seven different campaigns over the past 50 years, as well as point constraints across an array of petroleum industry drill holes. Particular attention is given to a detailed representation of Quaternary stratigraphy, representing shallow (z<250m) near-surface layers in the model. Seismic velocities are obtained from seismic reflection processing (for Vp) and also recently performed active and passive surface wave analyses (for Vs). Over 1,700 water wells in the region are used to constrain the complex inter-bedded Quaternary stratigraphy (gravels, sands, silts, organics etc.) near the coastline, including beneath urban Christchurch, which has resulted from fluvial deposition and marine regression and transgression. For the near-surface Springston and Christchurch Formations in the Christchurch urban area (z<50m), high-spatial resolution seismic velocities (including Vs30 ) were obtained from over 13,000 cone penetration tests combined with a recently developed CPT-Vs correlation.

Research papers, University of Canterbury Library

Liquefaction-induced lateral spreading during the 2011 Christchurch earthquake in New Zealand was severe and extensive, and data regarding the displacements associated with the lateral spreading provides an excellent opportunity to better understand the factors that influence these movements. Horizontal displacements measured from optical satellite imagery and subsurface data from the New Zealand Geotechnical Database (NZGD) were used to investigate four distinct lateral spread areas along the Avon River in Christchurch. These areas experienced displacements between 0.5 and 2 m, with the inland extent of displacement ranging from 100 m to over 600 m. Existing empirical and semi-empirical displacement models tend to under estimate displacements at some sites and over estimate at others. The integrated datasets indicate that the areas with more severe and spatially extensive displacements are associated with thicker and more laterally continuous deposits of liquefiable soil. In some areas, the inland extent of displacements is constrained by geologic boundaries and geomorphic features, as expressed by distinct topographic breaks. In other areas the extent of displacement is influenced by the continuity of liquefiable strata or by the presence of layers that may act as vertical seepage barriers. These observations demonstrate the need to integrate geologic/geomorphic analyses with geotechnical analyses when assessing the potential for lateral spreading movements.

Research papers, University of Canterbury Library

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.

Research papers, University of Canterbury Library

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.

Research papers, University of Canterbury Library

The research is funded by Callaghan Innovation (grant number MAIN1901/PROP-69059-FELLOW-MAIN) and the Ministry of Transport New Zealand in partnership with Mainfreight Limited. Need – The freight industry is facing challenges related to climate change, including natural hazards and carbon emissions. These challenges impact the efficiency of freight networks, increase costs, and negatively affect delivery times. To address these challenges, freight logistics modelling should consider multiple variables, such as natural hazards, sustainability, and emission reduction strategies. Freight operations are complex, involving various factors that contribute to randomness, such as the volume of freight being transported, the location of customers, and truck routes. Conventional methods have limitations in simulating a large number of variables. Hence, there is a need to develop a method that can incorporate multiple variables and support freight sustainable development. Method - A minimal viable model (MVM) method was proposed to elicit tacit information from industrial clients for building a minimally sufficient simulation model at the early modelling stages. The discrete-event simulation (DES) method was applied using Arena® software to create simulation models for the Auckland and Christchurch corridor, including regional pick-up and delivery (PUD) models, Christchurch city delivery models, and linehaul models. Stochastic variables in freight operations such as consignment attributes, customer locations, and truck routes were incorporated in the simulation. The geographic information system (GIS) software ArcGIS Pro® was used to identify and analyse industrial data. The results obtained from the GIS software were applied to create DES models. Life cycle assessment (LCA) models were developed for both diesel and battery electric (BE) trucks to compare their life cycle greenhouse gas (GHG) emissions and total cost of ownership (TCO) and support GHG emissions reduction. The line-haul model also included natural hazards in several scenarios, and the simulation was used to forecast the stock level of Auckland and Christchurch depots in response to each corresponding scenario. Results – DES is a powerful technique that can be employed to simulate and evaluate freight operations that exhibit high levels of variability, such as regional pickup and delivery (PUD) and linehaul. Through DES, it becomes possible to analyse multiple factors within freight operations, including transportation modes, routes, scheduling, and processing times, thereby offering valuable insights into the performance, efficiency, and reliability of the system. In addition, GIS is a useful tool for analysing and visualizing spatial data in freight operations. This is exemplified by their ability to simulate the travelling salesman problem (TSP) and conduct cluster analysis. Consequently, the integration of GIS into DES modelling is essential for improving the accuracy and reliability of freight operations analysis. The outcomes of the simulation were utilised to evaluate the ecological impact of freight transport by performing emission calculations and generating low-carbon scenarios to identify approaches for reducing the carbon footprint. LCA models were developed based on simulation results. Results showed that battery-electric trucks (BE) produced more greenhouse gas (GHG) emissions in the cradle phase due to battery manufacturing but substantially less GHG emissions in the use phase because of New Zealand's mostly renewable energy sources. While the transition to BE could significantly reduce emissions, the financial aspect is not compelling, as the total cost of ownership (TCO) for the BE truck was about the same for ten years, despite a higher capital investment for the BE. Moreover, external incentives are necessary to justify a shift to BE trucks. By using simulation methods, the effectiveness of response plans for natural hazards can be evaluated, and the system's vulnerabilities can be identified and mitigated to minimize the risk of disruption. Simulation models can also be utilized to simulate adaptation plans to enhance the system's resilience to natural disasters. Novel contributions – The study employed a combination of DES and GIS methods to incorporate a large number of stochastic variables and driver’s decisions into freight logistics modelling. Various realistic operational scenarios were simulated, including customer clustering and PUD truck allocation. This showed that complex pickup and delivery routes with high daily variability can be represented using a model of roads and intersections. Geographic regions of high customer density, along with high daily variability could be represented by a two-tier architecture. The method could also identify delivery runs for a whole city, which has potential usefulness in market expansion to new territories. In addition, a model was developed to address carbon emissions and total cost of ownership of battery electric trucks. This showed that the transition was not straightforward because the economics were not compelling, and that policy interventions – a variety were suggested - could be necessary to encourage the transition to decarbonised freight transport. A model was developed to represent the effect of natural disasters – such as earthquake and climate change – on road travel and detour times in the line haul freight context for New Zealand. From this it was possible to predict the effects on stock levels for a variety of disruption scenarios (ferry interruption, road detours). Results indicated that some centres rather than others may face higher pressure and longer-term disturbance after the disaster subsided. Remedies including coastal shipping were modelled and shown to have the potential to limit the adverse effects. A philosophical contribution was the development of a methodology to adapt the agile method into the modelling process. This has the potential to improve the clarification of client objectives and the validity of the resulting model.

Research papers, University of Canterbury Library

This dissertation addresses a diverse range of topics in the area of physics-based ground motion simulation with particular focus on the Canterbury, New Zealand region. The objectives achieved provide the means to perform hybrid broadband ground motion simulation and subsequently validates the simulation methodology employed. In particu- lar, the following topics are addressed: the development of a 3D seismic velocity model of the Canterbury region for broadband ground motion simulation; the development of a 3D geologic model of the interbedded Quaternary formations to provide insight on observed ground motions; and the investigation of systematic effects through ground motion sim- ulation of small-to-moderate magnitude earthquakes. The paragraphs below outline each contribution in more detail. As a means to perform hybrid broadband ground motion simulation, a 3D model of the geologic structure and associated seismic velocities in the Canterbury region is devel- oped utilising data from depth-converted seismic reflection lines, petroleum and water well logs, cone penetration tests, and implicitly guided by existing contour maps and geologic cross sections in data sparse subregions. The model explicitly characterises five significant and regionally recognisable geologic surfaces that mark the boundaries between geologic units with distinct lithology and age, including the Banks Peninsula volcanics, which are noted to strongly influence seismic wave propagation. The Basement surface represents the base of the Canterbury sedimentary basin, where a large impedance contrast exists re- sulting in basin-generated waves. Seismic velocities for the lithological units between the geologic surfaces are derived from well logs, seismic reflection surveys, root mean square stacking velocities, empirical correlations, and benchmarked against a regional crustal model, thus providing the necessary information for a Canterbury velocity model for use in broadband seismic wave propagation. A 3D high-resolution model of the Quaternary geologic stratigraphic sequence in the Canterbury region is also developed utilising datasets of 527 high-quality water well logs, and 377 near-surface cone penetration test records. The model, developed using geostatistical Kriging, represents the complex interbedded regional Quaternary geology by characterising the boundaries between significant interbedded geologic formations as 3D surfaces including explicit modelling of the formation unconformities resulting from the Banks Peninsula volcanics. The stratigraphic layering present can result in complex wave propagation. The most prevalent trend observed in the surfaces was the downward dip from inland to the eastern coastline as a result of the dominant fluvial depositional environment of the terrestrial gravel formations. The developed model provides a benefi- cial contribution towards developing a comprehensive understanding of recorded ground motions in the region and also providing the necessary information for future site char- acterisation and site response analyses. To highlight the practicality of the model, an example illustrating the role of the model in constraining surface wave analysis-based shear wave velocity profiling is illustrated along with the calculation of transfer functions to quantify the effect of the interbedded geology on wave propagation. Lastly, an investigation of systematic biases in the (Graves and Pitarka, 2010, 2015) ground motion simulation methodology and the specific inputs used for the Canterbury region is presented considering 144 small-to-moderate magnitude earthquakes. In the simulation of these earthquakes, the 3D Canterbury Velocity Model, developed as a part of this dissertation, is used for the low-frequency simulation, and a regional 1D velocity model for the high-frequency simulation. Representative results for individual earthquake sources are first presented to highlight the characteristics of the small-to-moderate mag- nitude earthquake simulations through waveforms, intensity measure scaling with source- to-site distance, and spectral bias of the individual events. Subsequently, a residual de- composition is performed to examine the between- and within-event residuals between observed data, and simulated and empirical predictions. By decomposing the residuals into between- and within-event residuals, the biases in source, path and site effects, and their causes, can be inferred. The residuals are comprehensively examined considering their aggregated characteristics, dependence on predictor variables, spatial distribution, and site-specific effects. The results of the simulation are also benchmarked against empir- ical ground motion models, where their similarities manifest from common components in their prediction. Ultimately, suggestions to improve the predictive capability of the simulations are presented as a result of the analysis.

Research papers, The University of Auckland Library

Critical infrastructure networks are highly relied on by society such that any disruption to service can have major social and economic implications. Furthermore, these networks are becoming increasingly dependent on each other for normal operation such that an outage or asset failure in one system can easily propagate and cascade across others resulting in widespread disruptions in terms of both magnitude and spatial reach. It is the vulnerability of these networks to disruptions and the corresponding complexities in recovery processes which provide direction to this research. This thesis comprises studies contributing to two areas (i) the modelling of national scale in-terdependent infrastructure systems undergoing major disruptions, and (ii) the tracking and quantification of infrastructure network recovery trajectories following major disruptions. Firstly, methods are presented for identifying nationally significant systemic vulnerabilities and incorporating expert knowledge into the quantification of infrastructure interdependency mod-elling and simulation. With application to the interdependent infrastructures networks across New Zealand, the magnitudes and spatial extents of disruption are investigated. Results high-light the importance in considering interdependencies when assessing disruptive risks and vul-nerabilities in disaster planning applications and prioritising investment decisions for enhancing resilience of national networks. Infrastructure dependencies are further studied in the context of recovery from major disruptions through the analysis of curves measuring network functionality over time. Continued studies into the properties of recovery curves across a database of global natural disasters produce statistical models for predicting the trajectory and expected recovery times. Finally, the use of connectivity based metrics for quantifying infrastructure system functionality during recovery are considered with a case study application to the Christchurch Earthquake (February 22, 2011) wastewater network response.

Research papers, University of Canterbury Library

In September 2010 and February 2011, the Canterbury region experienced devastating earthquakes with an estimated economic cost of over NZ$40 billion (Parker and Steenkamp, 2012; Timar et al., 2014; Potter et al., 2015). The insurance market played an important role in rebuilding the Canterbury region after the earthquakes. Homeowners, insurance and reinsurance markets and New Zealand government agencies faced a difficult task to manage the rebuild process. From an empirical and theoretic research viewpoint, the Christchurch disaster calls for an assessment of how the insurance market deals with such disasters in the future. Previous studies have investigated market responses to losses in global catastrophes by focusing on the insurance supply-side. This study investigates both demand-side and supply-side insurance market responses to the Christchurch earthquakes. Despite the fact that New Zealand is prone to seismic activities, there are scant previous studies in the area of earthquake insurance. This study does offer a unique opportunity to examine and document the New Zealand insurance market response to catastrophe risk, providing results critical for understanding market responses after major loss events in general. A review of previous studies shows higher premiums suppress demand, but how higher premiums and a higher probability of risk affect demand is still largely unknown. According to previous studies, the supply of disaster coverage is curtailed unless the market is subsidised, however, there is still unsettled discussion on why demand decreases with time from the previous disaster even when the supply of coverage is subsidised by the government. Natural disaster risks pose a set of challenges for insurance market players because of substantial ambiguity associated with the probability of such events occurring and high spatial correlation of catastrophe losses. Private insurance market inefficiencies due to high premiums and spatially concentrated risks calls for government intervention in the provision of natural disaster insurance to avert situations of noninsurance and underinsurance. Political economy considerations make it more likely for government support to be called for if many people are uninsured than if few people are uninsured. However, emergency assistance for property owners after catastrophe events can encourage most property owners to not buy insurance against natural disaster and develop adverse selection behaviour, generating larger future risks for homeowners and governments. On the demand-side, this study has developed an intertemporal model to examine how demand for insurance changes post-catastrophe, and how to model it theoretically. In this intertemporal model, insurance can be sought in two sequential periods of time, and at the second period, it is known whether or not a loss event happened in period one. The results show that period one demand for insurance increases relative to the standard single period model when the second period is taken into consideration, period two insurance demand is higher post-loss, higher than both the period one demand and the period two demand without a period one loss. To investigate policyholders experience from the demand-side perspective, a total of 1600 survey questionnaires were administered, and responses from 254 participants received representing a 16 percent response rate. Survey data was gathered from four institutions in Canterbury and is probably not representative of the entire population. The results of the survey show that the change from full replacement value policy to nominated replacement value policy is a key determinant of the direction of change in the level of insurance coverage after the earthquakes. The earthquakes also highlighted the plight of those who were underinsured, prompting policyholders to update their insurance coverage to reflect the estimated cost of re-building their property. The survey has added further evidence to the existing literature, such as those who have had a recent experience with disaster loss report increased risk perception if a similar event happens in future with females reporting a higher risk perception than males. Of the demographic variables, only gender has a relationship with changes in household cover. On the supply-side, this study has built a risk-based pricing model suitable to generate a competitive premium rate for natural disaster insurance cover. Using illustrative data from the Christchurch Red-zone suburbs, the model generates competitive premium rates for catastrophe risk. When the proposed model incorporates the new RMS high-definition New Zealand Earthquake Model, for example, insurers can find the model useful to identify losses at a granular level so as to calculate the competitive premium. This study observes that the key to the success of the New Zealand dual insurance system despite the high prevalence of catastrophe losses are; firstly the EQC’s flat-rate pricing structure keeps private insurance premiums affordable and very high nationwide homeowner take-up rates of natural disaster insurance. Secondly, private insurers and the EQC have an elaborate reinsurance arrangement in place. By efficiently transferring risk to the reinsurer, the cost of writing primary insurance is considerably reduced ultimately expanding primary insurance capacity and supply of insurance coverage.

Research papers, University of Canterbury Library

This manuscript provides a critical examination of the ground motions recorded in the near-source region resulting from the 22 February 2011 Christchurch earthquake. Particular attention is given to reconciling the observed spatial distribution of ground motions in terms of physical phenomena related to source, path and site effects. The large number of near-source observed strong ground motions show clear evidence of: forward-directivity, basin generated surface waves, liquefaction and other significant nonlinear site response. The pseudo-acceleration response spectra (SA) amplitudes and significant duration of strong motions agree well with empirical prediction models, except at long vibration periods where the influence of basin-generated surface waves and nonlinear site response are significant and not adequately accounted for in empirical SA models. Pseudo-acceleration response spectra are also compared with those observed in the 4 September 2010 Darfield earthquake and routine design response spectra used in order to emphasise the amplitude of ground shaking and elucidate the importance of local geotechnical characteristics on surface ground motions. The characteristics of the observed vertical component accelerations are shown to be strongly dependent on source-to-site distance and are comparable with those from the 4 September 2010 Darfield earthquake, implying the large amplitudes observed are simply a result of many observations at close distances rather than a peculiar source effect.

Research papers, Victoria University of Wellington

The last seven years have seen southern New Zealand a ected by several large and damaging earthquakes: the moment magnitude (MW) 7.8 Dusky Sound earthquake on 15 July 2009, the MW 7.1 Dar eld (Canterbury) earthquake on 4 September 2010, and most notably the MW 6.2 Christchurch earthquake on 22 February 2011 and the protracted aftershock sequence. In this thesis, we address the postseismic displacement produced by these earthquakes using methods of satellite-based geodetic measurement, known as Interferometric Synthetic Aperture Radar (InSAR) and Global Positioning System (GPS), and computational modelling.  We observe several ground displacement features in the Canterbury and Fiordland regions during three periods: 1) Following the Dusky Sound earthquake; 2) Following the Dar eld earthquake and prior to the Christchurch earthquake; and 3) Following the Christchurch earthquake until February 2015.  The ground displacement associated with postseismic motion following the Dusky Sound earthquake has been measured by continuous and campaign GPS data acquired in August 2009, in conjunction with Di erential Interferometric Synthetic Aperture Radar (DInSAR) observations. We use an afterslip model, estimated by temporal inversion of geodetic data, with combined viscoelastic rebound model to account for the observed spatio-temporal patterns of displacement. The two postseismic processes together induce a signi cant displacement corresponding to principal extensional and contractual strain rates of the order of 10⁻⁷ and 10⁻⁸ yr⁻¹ respectively, across most of the southern South Island.  We also analyse observed postseismic displacement following the Dusky Sound earthquake using a new inversion approach in order to describe afterslip in an elasticviscoelastic medium. We develop a mathematical framework, namely the "Iterative Decoupling of Afterslip and Viscoelastic rebound (IDAV)" method, with which to invert temporally dense and spatially sparse geodetic observations. We examine the IDAV method using both numerical and analytical simulations of Green's functions.  For the post-Dar eld time interval, postseismic signals are measured within approximately one month of the mainshock. The dataset used for the post-Dar eld displacement spans the region surrounding previously unrecognised faults that ruptured during the mainshock. Poroelastic rebound in a multi-layered half-space and dilatancy recovery at shallow depths provide a satisfactory t with the observations.  For the post-Christchurch interval, campaign GPS data acquired in February 2012 to February 2015 in four successive epochs and 66 TerraSAR-X (TSX) SAR acquisitions in descending orbits between March 2011 and May 2014 reveal approximately three years of postseismic displacement. We detect movement away from the satellite of ~ 3 mm/yr in Christchurch and a gradient of displacement of ~ 4 mm/yr across a lineament extending from the westernmost end of the Western Christchurch Fault towards the eastern end of the Greendale East Fault. The postseismic signals following the Christchurch earthquake are mainly accounted for by afterslip models on the subsurface lineament and nearby faults.

Research papers, University of Canterbury Library

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.

Research papers, University of Canterbury Library

Collective identity construction in organisations engaged in an inter-organisational collaboration (IOC), especially temporary IOCs set up in disaster situations, has received scant attention in the organisational studies literature yet collective identity is considered to be important in fostering effective IOC operations. This doctoral study was designed to add to our understanding about how collective identity is constituted throughout the entire lifespan of a particular temporary coopetitive (i.e., simultaneously collaborative and competitive) IOC formed in a post-disaster environment. To achieve this purpose, a qualitative case study of the Stronger Christchurch Infrastructure Rebuild Team (SCIRT), a time-bound coopetition formed to repair the horizontal infrastructure in Christchurch, New Zealand after the devastating 2011 Canterbury earthquakes, was undertaken. Using data from semi-structured interviews, field observations, and organisational documents and other artefacts, an inductive analytic method was employed to explore how internal stakeholders engaged with and co- constructed a collective SCIRT identity and reconciled this with their home organization identity. The analysis revealed that the SCIRT collective identity was an ongoing process, involving the interweaving of social, temporal, material and geospatial dimensions constructed through intersecting cycles of senior managers’ sensegiving and employees’ sensemaking across SCIRT’s five and a half years of existence. Senior management deliberately undertook identity work campaigns that used organisational rituals, artefacts, and spatial design to disseminate and encourage a sense of “we are all SCIRT”. However, there was no common sense of “we-ness”. Identification with SCIRT was experienced differently among different groups of employees and across time. Employees’ differing senses of collective identity were accounted for by their past, present, and anticipated future relationships with their home organisation, and also (re)shaped by the geosocial environments in which they worked. The study supports previous research claiming that collective identity is a process of recursive sensegiving and sensemaking between senior managers and employees. However, it extends the literature by revealing the imbricated nature of collective identity, how members’ sense of “who we are” can change across the entire lifetime of a temporary IOC, and how sociomateriality, temporality, and geosocial effects strongly intervene in employees’ emerging senses of collective identity. Moreover, the study demonstrates how the ongoing identity work can be embedded in a time-space frame that further accentuates the influence of temporality, especially the anticipated future, organisational rituals, artefacts, and the geosocial environment. The study’s primary contribution to theory is a processual model of collective identity that applies specifically to a temporary IOC involving coopetition. In doing so, it represents a more finely nuanced and situational model than existing models. At a practical level, this model suggests that managers need to appreciate that organisational artefacts, rituals, and the prevailing organisational geosocial environment are inextricably linked in processes that can be manipulated to enhance the construction of collective identity.

Research papers, University of Canterbury Library

Bulk rock strength is greatly dependent on fracture density, so that reductions in rock strength associated with faulting and fracturing should be reflected by reduced shear coupling and hence S-wave velocity. This study is carried out along the Canterbury rangefront and in Otago. Both lie within the broader plate boundary deformation zone in the South Island of New Zealand. Therefore built structures are often, , located in areas where there are undetected or poorly defined faults with associated rock strength reduction. Where structures are sited near to, or across, such faults or fault-zones, they may sustain both shaking and ground deformation damage during an earthquake. Within this zone, management of seismic hazards needs to be based on accurate identification of the potential fault damage zone including the likely width of off-plane deformation. Lateral S-wave velocity variability provides one method of imaging and locating damage zones and off-plane deformation. This research demonstrates the utility of Multi-Channel Analysis of Surface Waves (MASW) to aid land-use planning in such fault-prone settings. Fundamentally, MASW uses surface wave dispersive characteristics to model a near surface profile of S-wave velocity variability as a proxy for bulk rock strength. The technique can aid fault-zone planning not only by locating and defining the extent of fault-zones, but also by defining within-zone variability that is readily correlated with measurable rock properties applicable to both foundation design and the distribution of surface deformation. The calibration sites presented here have well defined field relationships and known fault-zone exposure close to potential MASW survey sites. They were selected to represent a range of progressively softer lithologies from intact and fractured Torlesse Group basement hard rock (Dalethorpe) through softer Tertiary cover sediments (Boby’s Creek) and Quaternary gravels. This facilitated initial calibration of fracture intensity at a high-velocity-contrast site followed by exploration of the limits of shear zone resolution at lower velocity contrasts. Site models were constructed in AutoCAD in order to demonstrate spatial correlations between S-wave velocity and fault zone features. Site geology was incorporated in the models, along with geomorphology, river profiles, scanline locations and crosshole velocity measurement locations. Spatial data were recorded using a total-station survey. The interpreted MASW survey results are presented as two dimensional snapshot cross-sections of the three dimensional calibration-site models. These show strong correlations between MASW survey velocities and site geology, geomorphology, fluvial profiles and geotechnical parameters and observations. Correlations are particularly pronounced where high velocity contrasts exist, whilst weaker correlations are demonstrated in softer lithologies. Geomorphic correlations suggest that off-plane deformation can be imaged and interpreted in the presence of suitable topographic survey data. A promising new approach to in situ and laboratory soft-rock material and mass characterisation is also presented using a Ramset nail gun. Geotechnical investigations typically involve outcrop and laboratory scale determination of rock mass and material properties such as fracture density and unconfined compressive strength (UCS). This multi-scale approach is espoused by this study, with geotechnical and S-wave velocity data presented at multiple scales, from survey scale sonic velocity measurements, through outcrop scale scanline and crosshole sonic velocity measurements to laboratory scale property determination and sonic velocity measurements. S-wave velocities invariably increased with decreasing scale. These scaling relationships and strategies for dealing with them are investigated and presented. Finally, the MASW technique is applied to a concealed fault on the Taieri Ridge in Macraes Flat, Central Otago. Here, high velocity Otago Schist is faulted against low velocity sheared Tertiary and Quaternary sediments. This site highlights the structural sensitivity of the technique by apparently constraining the location of the principal fault, which had been ambiguous after standard processing of the seismic reflection data. Processing of the Taieri Ridge dataset has further led to the proposal of a novel surface wave imaging technique termed Swept Frequency Imaging (SFI). This inchoate technique apparently images the detailed structure of the fault-zone, and is in agreement with the conventionally-determined fault location and an existing partial trench. Overall, the results are promising and are expected to be supported by further trenching in the near future.

Research papers, University of Canterbury Library

As a global phenomenon, many cities are undergoing urban renewal to accommodate rapid growth in urban population. However, urban renewal can struggle to balance social, economic, and environmental outcomes, whereby economic outcomes are often primarily considered by developers. This has important implications for urban forests, which have previously been shown to be negatively affected by development activities. Urban forests serve the purpose of providing ecosystem services and thus are beneficial to human wellbeing. Better understanding the effect of urban renewal on city trees may help improve urban forest outcomes via effective management and policy strategies, thereby maximising ecosystem service provision and human wellbeing. Though the relationship between certain aspects of development and urban forests has received consideration in previous literature, little research has focused on how the complete property redevelopment cycle affects urban forest dynamics over time. This research provides an opportunity to gain a comprehensive understanding of the effect of residential property redevelopment on urban forest dynamics, at a range of spatial scales, in Christchurch, New Zealand following a series of major earthquakes which occurred in 2010 – 2011. One consequence of the earthquakes is the redevelopment of thousands of properties over a relatively short time-frame. The research quantifies changes in canopy cover city-wide, as well as, tree removal, retention, and planting on individual residential properties. Moreover, the research identifies the underlying reasons for these dynamics, by exploring the roles of socio-economic and demographic factors, the spatial relationships between trees and other infrastructure, and finally, the attitudes of residential property owners. To quantify the effect of property redevelopment on canopy cover change in Christchurch, this research delineated tree canopy cover city-wide in 2011 and again in 2015. An object-based image analysis (OBIA) technique was applied to aerial imagery and LiDAR data acquired at both time steps, in order to estimate city-wide canopy cover for 2011 and 2015. Changes in tree canopy cover between 2011 and 2015 were then spatially quantified. Tree canopy cover change was also calculated for all meshblocks (a relatively fine-scale geographic boundary) in Christchurch. The results show a relatively small magnitude of tree canopy cover loss, city-wide, from 10.8% to 10.3% between 2011 and 2015, but a statistically significant change in mean tree canopy cover across all the meshblocks. Tree canopy cover losses were more likely to occur in meshblocks containing properties that underwent a complete redevelopment cycle, but the loss was insensitive to the density of redevelopment within meshblocks. To explore property-scale individual tree dynamics, a mixed-methods approach was used, combining questionnaire data and remote sensing analysis. A mail-based questionnaire was delivered to residential properties to collect resident and household data; 450 residential properties (321 redeveloped, 129 non- redeveloped) returned valid questionnaires and were identified as analysis subjects. Subsequently, 2,422 tree removals and 4,544 tree retentions were identified within the 450 properties; this was done by manually delineating individual tree crowns, based on aerial imagery and LiDAR data, and visually comparing the presence or absence of these trees between 2011 and 2015. The tree removal rate on redeveloped properties (44.0%) was over three times greater than on non-redeveloped properties (13.5%) and the average canopy cover loss on redeveloped properties (52.2%) was significantly greater than on non-redeveloped properties (18.8%). A classification tree (CT) analysis was used to model individual tree dynamics (i.e. tree removal, tree retention) and candidate explanatory variables (i.e. resident and household, economic, land cover, and spatial variables). The results indicate that the model including land cover, spatial, and economic variables had the best predicting ability for individual tree dynamics (accuracy = 73.4%). Relatively small trees were more likely to be removed, while trees with large crowns were more likely to be retained. Trees were most likely to be removed from redeveloped properties with capital values lower than NZ$1,060,000 if they were within 1.4 m of the boundary of a redeveloped building. Conversely, trees were most likely to be retained if they were on a property that was not redeveloped. The analysis suggested that the resident and household factors included as potential explanatory variables did not influence tree removal or retention. To conduct a further exploration of the relationship between resident attitudes and actions towards trees on redeveloped versus non-redeveloped properties, this research also asked the landowners from the 450 properties that returned mail questionnaires to indicate their attitudes towards tree management (i.e. tree removal, tree retention, and tree planting) on their properties. The results show that residents from redeveloped properties were more likely to remove and/or plant trees, while residents from non- redeveloped properties were more likely to retain existing trees. A principal component analysis (PCA) was used to explore resident attitudes towards tree management. The results of the PCA show that residents identified ecosystem disservices (e.g. leaf litter, root damage to infrastructure) as common reasons for tree removal; however, they also noted ecosystem services as important reasons for both tree planting and tree retention on their properties. Moreover, the reasons for tree removal and tree planting varied based on whether residents’ property had been redeveloped. Most tree removal occurred on redeveloped properties because trees were in conflict with redevelopment, but occurred on non- redeveloped properties because of perceived poor tree health. Residents from redeveloped properties were more likely to plant trees due to being aesthetically pleasing or to replace trees removed during redevelopment. Overall, this research adds to, and complements, the existing literature on the effects of residential property redevelopment on urban forest dynamics. The findings of this research provide empirical support for developing specific legislation or policies about urban forest management during residential property redevelopment. The results also imply that urban foresters should enhance public education on the ecosystem services provided by urban forests and thus minimise the potential for tree removal when undertaking property redevelopment.

Research papers, University of Canterbury Library

Asset management in power systems is exercised to improve network reliability to provide confidence and security for customers and asset owners. While there are well-established reliability metrics that are used to measure and manage business-as-usual disruptions, an increasing appreciation of the consequences of low-probability high-impact events means that resilience is increasingly being factored into asset management in order to provide robustness and redundancy to components and wider networks. This is particularly important for electricity systems, given that a range of other infrastructure lifelines depend upon their operation. The 2010-2011 Canterbury Earthquake Sequence provides valuable insights into electricity system criticality and resilience in the face of severe earthquake impacts. While above-ground assets are relatively easy to monitor and repair, underground assets such as cables emplaced across wide areas in the distribution network are difficult to monitor, identify faults on, and repair. This study has characterised in detail the impacts to buried electricity cables in Christchurch resulting from seismically-induced ground deformation caused primarily by liquefaction and lateral spread. Primary modes of failure include cable bending, stretching, insulation damage, joint braking and, being pulled off other equipment such as substation connections. Performance and repair data have been compiled into a detailed geospatial database, which in combination with spatial models of peak ground acceleration, peak ground velocity and ground deformation, will be used to establish rigorous relationships between seismicity and performance. These metrics will be used to inform asset owners of network performance in future earthquakes, further assess component criticality, and provide resilience metrics.