This study examines the performance of nonlinear total-stress wave-propagation site response analysis for modelling site effects in physics-based ground motion simulations of the 2010-2011 Canterbury, New Zealand earthquake sequence. This approach allows for explicit modeling of 3-dimensional ground motion phenomena at the regional scale, as well as detailed site effects and soil nonlinearity at the local scale. The approach is compared to a more commonly used empirical VS30 (30 m time-averaged shear wave velocity)-based method for computing site amplification as proposed by Graves and Pitarka (2010, 2015).
SeisFinder is an open-source web service developed by QuakeCoRE and the University of Canterbury, focused on enabling the extraction of output data from computationally intensive earthquake resilience calculations. Currently, SeisFinder allows users to select historical or future events and retrieve ground motion simulation outputs for requested geographical locations. This data can be used as input for other resilience calculations, such as dynamic response history analysis. SeisFinder was developed using Django, a high-level python web framework, and uses a postgreSQL database. Because our large-scale computationally-intensive numerical ground motion simulations produce big data, the actual data is stored in file systems, while the metadata is stored in the database.
Overview of SeisFinder SeisFinder is an open-source web service developed by QuakeCoRE and the University of Canterbury, focused on enabling the extraction of output data from computationally intensive earthquake resilience calculations. Currently, SeisFinder allows users to select historical or future events and retrieve ground motion simulation outputs for requested geographical locations. This data can be used as input for other resilience calculations, such as dynamic response history analysis. SeisFinder was developed using Django, a high-level python web framework, and uses a postgreSQL database. Because our large-scale computationally-intensive numerical ground motion simulations produce big data, the actual data is stored in file systems, while the metadata is stored in the database. The basic SeisFinder architecture is shown in Figure 1.
The 2010–2011 Canterbury earthquakes and their aftermath have been described by the Human Rights Commission as one of New Zealand's greatest contemporary human rights challenges. This article documents the shortcomings in the realisation of the right to housing in post-quake Canterbury for homeowners, tenants and the homeless. The article then considers what these shortcomings tell us about New Zealand's overall human rights framework, suggesting that the ongoing and seemingly intractable nature of these issues and the apparent inability to resolve them indicate an underlying fragility implicit in New Zealand's framework for dealing with the consequences of a large-scale natural disaster. The article concludes that there is a need for a comprehensive human rights-based approach to disaster preparedness, response and recovery in New Zealand.
Following the Canterbury earthquake sequence of 2010-11, a large and contiguous tract of vacated ‘red zoned’ land lies alongside the lower Ōtākaro / Avon River and is known as the Avon-Ōtākaro Red Zone (AORZ). This is the second report in the Ecological Regeneration Options (ERO) project that addresses future land uses in the AORZ. The purpose of this report is to present results from an assessment of restoration opportunities conducted in April 2017. The objectives of the assessment were to identify potential benefits of ecological restoration activities across both land and water systems in the AORZ and characterise the key options for their implementation. The focus of this report is not to provide specific advice on the methods for achieving specific restoration endpoints per se. This will vary at different sites and scales with a large number of combinations possible. Rather, the emphasis is on providing an overview of the many restoration and regeneration options in their totality across the AORZ. An additional objective is to support their adequate assessment in the identification of optimum land uses and adaptive management practices for the AORZ. Participatory processes may play a useful role in assessment and stakeholder engagement by providing opportunities for social learning and the co-creation of new knowledge. We used a facilitated local knowledge based approach that generated a large quantity of reliable and site specific data in a short period of time. By inviting participation from a wide knowledge-holder network inclusivity is improved in comparison to small-group expert panel approaches. Similar approaches could be applied to other information gathering and assessment needs in the regeneration planning process. Findings from this study represent the most comprehensive set of concepts available to date to address the potential benefits of ecological regeneration in the AORZ. This is a core topic for planning to avoid missed opportunities and opportunity costs. The results identify a wide range of activities that may be applied to generate benefits for Christchurch and beyond, all involving aspects of a potential new ecology in the AORZ. These may be combined at a range of scales to create scenarios, quantify benefits, and explore the potential for synergies between different land use options. A particular challenge is acquiring the information needed within relatively short time frames. Early attention to gathering baseline data, addressing technical knowledge gaps, and developing conceptual frameworks to account for the many spatio-temporal aspects are all key activities that will assist in delivering the best outcomes. Methodologies by which these many facets can be pulled together in quantitative and comparative assessments are the focus of the final report in the ERO series.
The Canterbury Earthquake Sequence 2010-2011 (CES) induced widespread liquefaction in many parts of Christchurch city. Liquefaction was more commonly observed in the eastern suburbs and along the Avon River where the soils were characterised by thick sandy deposits with a shallow water table. On the other hand, suburbs to the north, west and south of the CBD (e.g. Riccarton, Papanui) exhibited less severe to no liquefaction. These soils were more commonly characterised by inter-layered liquefiable and non-liquefiable deposits. As part of a related large-scale study of the performance of Christchurch soils during the CES, detailed borehole data including CPT, Vs and Vp have been collected for 55 sites in Christchurch. For this subset of Christchurch sites, predictions of liquefaction triggering using the simplified method (Boulanger & Idriss, 2014) indicated that liquefaction was over-predicted for 94% of sites that did not manifest liquefaction during the CES, and under-predicted for 50% of sites that did manifest liquefaction. The focus of this study was to investigate these discrepancies between prediction and observation. To assess if these discrepancies were due to soil-layer interaction and to determine the effect that soil stratification has on the develop-ment of liquefaction and the system response of soil deposits.
The Kaikoura earthquake in November 2016 highlighted the vulnerability of New Zealand’s rural communities to locally-specific hazard events, which generate regional and national scale impacts. Kaikoura was isolated with significant damage to both the east coast road (SH1) and rail corridor, and the Inland Road (Route 70). Sea bed uplift along the coast was significant – affecting marine resources and ocean access for marine operators engaged in tourism and harvesting, and recreational users. While communities closest to the earthquake epicentre (e.g., Kaikoura, Waiau, Rotherham and Cheviot) suffered the most immediate earthquake damage, the damage to the transport network, and the establishment of an alternative transport route between Christchurch and Picton, has significantly impacted on more distant communities (e.g., Murchison, St Arnaud and Blenheim). There was also considerable damage to vineyard infrastructure across the Marlborough region and damage to buildings and infrastructure in rural settlements in Southern Marlborough (e.g., Ward and Seddon).
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.
This article presents a quantitative case study on the site amplification effect observed at Heathcote Valley, New Zealand, during the 2010-2011 Canterbury earthquake sequence for 10 events that produced notable ground acceleration amplitudes up to 1.4g and 2.2g in the horizontal and vertical directions, respectively. We performed finite element analyses of the dynamic response of the valley, accounting for the realistic basin geometry and the soil non-linear response. The site-specific simulations performed significantly better than both empirical ground motion models and physics based regional-scale ground motion simulations (which empirically accounts for the site effects), reducing the spectral acceleration prediction bias by a factor of two in short vibration periods. However, our validation exercise demonstrated that it was necessary to quantify the level of uncertainty in the estimated bedrock motion using multiple recorded events, to understand how much the simplistic model can over- or under-estimate the ground motion intensities. Inferences from the analyses suggest that the Rayleigh waves generated near the basin edge contributed significantly to the observed high frequency (f>3Hz) amplification, in addition to the amplification caused by the strong soil-rock impedance contrast at the site fundamental frequency. Models with and without considering soil non-linear response illustrate, as expected, that the linear elastic assumption severely overestimates ground motions in high frequencies for strong earthquakes, especially when the contribution of basin edge-generated Rayleigh waves becomes significant. Our analyses also demonstrate that the effect of pressure-dependent soil velocities on the high frequency ground motions is as significant as the amplification caused by the basin edge-generated Rayleigh waves.
The Canterbury region of New Zealand experienced a sequence of strong earthquakes during 2010-2011. Responses included government acquisition of many thousands of residential properties in the city of Christchurch in areas with severe earthquake effects. A large and contiguous tract of this ‘red zoned’ land lies in close proximity to the Ōtākaro / Avon River and is known as the Avon-Ōtākaro Red Zone (AORZ). The focus of this study was to provide an overview of the floodplain characteristics of the AORZ and review of international experience in ecological restoration of similar river margin and floodplain ecosystems to extract restoration principles and associated learnings. Compared to pre-earthquake ground levels, the dominant trend in the AORZ is subsidence, together with lateral movement especially in the vicinity of waterway. An important consequence of land subsidence in the lower Ōtākaro / Avon River is greater exposure to flooding and the effects of sea level rise. Scenario modelling for sea level rise indicates that much of the AORZ is exposed to inundation within a 100 year planning horizon based on a 1 m sea level rise. As with decisions on built infrastructure, investments in nature-based ‘green infrastructure’ also require a sound business case including attention to risks posed by climate change. Future-proofing of the expected benefits of ecological restoration must therefore be secured by design. Understanding and managing the hydrology and floodplain dynamics are vital to the future of the AORZ. However, these characteristics are shared by other floodplain and river restoration projects worldwide. Identifying successful approaches provides a useful a source of useful information for floodplain planning in the AORZ. This report presents results from a comparative case study of three international examples to identify relevant principles for large-scale floodplain management at coastal lowland sites.
Post-traumatic stress symptoms are a common reaction to experiencing a traumatic event such as a natural disaster. Young children may be at an increased risk for such mental health problems as these catastrophic events may coincide with developmentally sensitive periods of development. Treatments currently recommended for children with post-traumatic stress symptoms insufficiently acknowledge the role of neurobiological stress related systems responsible for these symptoms. As such, alternative approaches to the treatment of posttraumatic symptoms have been explored, with nature-based interventions offering a potential alternative based on two different theories that uphold the stress reducing benefits of natural environments. To date, there are a limited number of experimental studies that have explored the use of nature-based interventions with children, and no known research that has used a simulated nature experience with child participants. The purpose of this study was to investigate the effects of a simulated nature experience on the physiological and behavioural responses of children with post-traumatic stress symptoms that experienced the Christchurch earthquakes. A single-case research design with repeated measures of heart rate and teacherreported behaviour was gathered across a 20-day period. Heart rate data was collected before and after participants watched a 10-minute nature video, while data from a teacher rating scale provided information about the participants’ behaviours in the 30-minute period after they watched the nature video. Comparisons made to data collected during two different baseline phases indicated that the nature video intervention had no recognisable effects on the participants’ physiological and behavioural stress responses. Limitations to the current study are discussed as possible reasons for the incompatibility between the current study’s results and the findings from previous research. Suggestions are made for any future replications of the study.
We present initial results from a set of three-dimensional (3D) deterministic earthquake ground motion simulations for the northern Canterbury plains, Christchurch and the Banks Peninsula region, which explicitly incorporate the effects of the surface topography. The simu-lations are done using Hercules, an octree-based finite-element parallel software for solving 3D seismic wave propagation problems in heterogeneous media under kinematic faulting. We describe the efforts undertaken to couple Hercules with the South Island Velocity Model (SIVM), which included changes to the SIVM code in order to allow for single repetitive que-ries and thus achieve a seamless finite-element meshing process within the end-to-end ap-proach adopted in Hercules. We present our selection of the region of interest, which corre-sponds to an area of about 120 km × 120 km, with the 3D model reaching a depth of 60 km. Initial simulation parameters are set for relatively high minimum shear wave velocity and a low maximum frequency, which we are progressively scaling up as computing resources permit. While the effects of topography are typically more important at higher frequencies and low seismic velocities, even at this initial stage of our efforts (with a maximum of 2 Hz and a mini-mum of 500 m/s), it is possible to observe the importance of the topography in the response of some key locations within our model. To highlight these effects we compare the results of the 3D topographic model with respect to those of a flat (squashed) 3D model. We draw rele-vant conclusions from the study of topographic effects during earthquakes for this region and describe our plans for future work.
This thesis seeks to examine how the integration of play, small toys specifically, and the use of solution-focused brief therapy techniques can affect the outcomes for primary school aged children undergoing counselling. The setting is a counselling agency in Christchurch, New Zealand. A qualitative research approach is used and the data analysed using a narrative inquiry approach. The context of this study is the counselling service of an agency where young children, adolescents and their families are helped and supported through a variety of life issues. The counselling the participants are offered uses a combination of a solution-focused and play therapy where the purpose is to encourage clients to find exceptions to their presenting problems and identify their preferred future. The aim of this study is to help the children navigate their problem through a better understanding of and the gaining of personal skills and strengths. Participants were invited to be part of this study through the agency waiting list. The four included presented with a variety of reasons for coming to counselling yet these proved similar to that which the agency has been routinely presented with in the aftermath of the Canterbury earthquakes from 2011 to present day. Each participant had the consent of their parents or caregivers to engage in this project. The participants themselves separately agreed to engage in a solution- focused counselling process where the counsellor also integrated the use of small toys as part of the course. Counselling sessions were audiotaped, aspects photographed and analysed with a specific focus on client engagement. Four key themes emerged as the participants explored their personal narrative. Firstly, the “I’m OK” theme depicted in their first scaling activity, secondly a recognition that things could indeed be better and they needed help. Thirdly, a realisation of their own strengths and skills and finally that the future was an optimistic place to look forward to. These themes are described and explained through descriptions of the participant’s stories as well as self-reflection by the researcher. Transcriptions of sessions are included as are excerpts from the research journal and photographs of the use of the small toys by the children.
The purpose of this paper is to empirically investigate the effects of a major disaster on the management of human resources in the construction sector. It sets out to identify the construction skills challenges and the factors that affected skills availability following the 2010/2011 earthquakes in Christchurch. It is hoped that this study will provide insights for on-going reconstruction and future disaster response with respect to the problem of skills shortages. Design/methodology/approach A triangulation method was adopted. The quantitative method, namely, a questionnaire survey, was employed to provide a baseline description. Field observations and interviews were used as a follow-up to ascertain issues and potential shortages over time. Three focus groups in the form of research workshops were convened to gain further insight into the feedback and to investigate the validity and applicability of the research findings. Findings The earthquakes in Christchurch had compounded the pre-existing skills shortages in the country due to heightened demand from reconstruction. Skills shortages primarily existed in seismic assessment and design for land and structures, certain trades, project management and site supervision. The limited technical capability available nationally, shortage of temporary accommodation to house additional workers, time needed for trainees to become skilled workers, lack of information about reconstruction workloads and lack of operational capacity within construction organisations, were critical constraints to the resourcing of disaster recovery projects. Research limitations/implications The research findings contribute to the debate on skills issues in construction. The study provides evidence that contributes to an improved understanding of the industry’s skills vulnerability and emerging issues that would likely exist after a major disaster in a resource-limited country such as New Zealand. Practical implications From this research, decision makers and construction organisations can gain a clear direction for improving the construction capacity and capability for on-going reconstruction. Factors that affected the post-earthquake skills availability can be considered by decision makers and construction organisations in their workforce planning for future disaster events. The recommendations will assist them in addressing skills shortages for on-going reconstruction. Originality/value Although the study is country-specific, the findings show the nature and scale of skills challenges the construction industry is likely to face following a major disaster, and the potential issues that may compound skills shortages. It provides lessons for other disaster-prone countries where the resource pool is small and a large number of additional workers are needed to undertake reconstruction.
Christchurch City Council (Council) is undertaking the Land Drainage Recovery Programme in order to assess the effects of the earthquakes on flood risk to Christchurch. In the course of these investigations it has become better understood that floodplain management should be considered in a multi natural hazards context. Council have therefore engaged the Jacobs, Beca, University of Canterbury, and HR Wallingford project team to investigate the multihazards in eastern areas of Christchurch and develop flood management options which also consider other natural hazards in that context (i.e. how other hazards contribute to flooding both through temporal and spatial coincidence). The study has three stages: Stage 1 Gap Analysis – assessment of information known, identification of gaps and studies required to fill the gaps. Stage 2 Hazard Studies – a gap filling stage with the studies identified in Stage 1. Stage 3 Collating, Optioneering and Reporting – development of options to manage flood risk. This present report is to document findings of Stage 1 and recommends the studies that should be completed for Stage 2. It has also been important to consider how Stage 3 would be delivered and the gaps are prioritised to provide for this. The level of information available and hazards to consider is extensive; requiring this report to be made up of five parts each identifying individual gaps. A process of identifying information for individual hazards in Christchurch has been undertaken and documented (Part 1) followed by assessing the spatial co-location (Part 2) and probabilistic presence of multi hazards using available information. Part 3 considers multi hazard presence both as a temporal coincidence (e.g. an earthquake and flood occurring at one time) and as a cascade sequence (e.g. earthquake followed by a flood at some point in the future). Council have already undertaken a number of options studies for managing flood risk and these are documented in Part 4. Finally Part 5 provides the Gap Analysis Summary and Recommendations to Council. The key findings of Stage 1 gap analysis are: - The spatial analysis showed eastern Christchurch has a large number of hazards present with only 20% of the study area not being affected by any of the hazards mapped. Over 20% of the study area is exposed to four or more hazards at the frequencies and data available. - The majority of the Residential Red Zone is strongly exposed to multiple hazards, with 86% of the area being exposed to 4 or more hazards, and 24% being exposed to 6 or more hazards. - A wide number of gaps are present; however, prioritisation needs to consider the level of benefit and risks associated with not undertaking the studies. In light of this 10 studies ranging in scale are recommended to be done for the project team to complete the present scope of Stage 3. - Stage 3 will need to consider a number of engineering options to address hazards and compare with policy options; however, Council have not established a consistent policy on managed retreat that can be applied for equal comparison; without which substantial assumptions are required. We recommend Council undertake a study to define a managed retreat framework as an option for the city. - In undertaking Stage 1 with floodplain management as the focal point in a multi hazards context we have identified that Stage 3 requires consideration of options in the context of economics, implementation and residual risk. Presently the scope of work will provide a level of definition for floodplain options; however, this will not be at equal levels of detail for other hazard management options. Therefore, we recommend Council considers undertaking other studies with those key hazards (e.g. Coastal Hazards) as a focal point and identifies the engineering options to address such hazards. Doing so will provide equal levels of information for Council to make an informed and defendable decision on which options are progressed following Stage 3.
This dissertation addresses several fundamental and applied aspects of ground motion selection for seismic response analyses. In particular, the following topics are addressed: the theory and application of ground motion selection for scenario earthquake ruptures; the consideration of causal parameter bounds in ground motion selection; ground motion selection in the near-fault region where directivity effect is significant; and methodologies for epistemic uncertainty consideration and propagation in the context of ground motion selection and seismic performance assessment. The paragraphs below outline each contribution in more detail. A scenario-based ground motion selection method is presented which considers the joint distribution of multiple intensity measure (IM) types based on the generalised conditional intensity measure (GCIM) methodology (Bradley, 2010b, 2012c). The ground motion selection algorithm is based on generating realisations of the considered IM distributions for a specific rupture scenario and then finding the prospective ground motions which best fit the realisations using an optimal amplitude scaling factor. In addition, using different rupture scenarios and site conditions, two important aspects of the GCIM methodology are scrutinised: (i) different weight vectors for the various IMs considered; and (ii) quantifying the importance of replicate selections for ensembles with different numbers of desired ground motions. As an application of the developed scenario-based ground motion selection method, ground motion ensembles are selected to represent several major earthquake scenarios in New Zealand that pose a significant seismic hazard, namely, Alpine, Hope and Porters Pass ruptures for Christchurch city; and Wellington, Ohariu, and Wairarapa ruptures for Wellington city. A rigorous basis is developed, and sensitivity analyses performed, for the consideration of bounds on causal parameters (e.g., magnitude, source-to-site distance, and site condition) for ground motion selection. The effect of causal parameter bound selection on both the number of available prospective ground motions from an initial empirical as-recorded database, and the statistical properties of IMs of selected ground motions are examined. It is also demonstrated that using causal parameter bounds is not a reliable approach to implicitly account for ground motion duration and cumulative effects when selection is based on only spectral acceleration (SA) ordinates. Specific causal parameter bounding criteria are recommended for general use as a ‘default’ bounding criterion with possible adjustments from the analyst based on problem-specific preferences. An approach is presented to consider the forward directivity effects in seismic hazard analysis, which does not separate the hazard calculations for pulse-like and non-pulse-like ground motions. Also, the ability of ground motion selection methods to appropriately select records containing forward directivity pulse motions in the near-fault region is examined. Particular attention is given to ground motion selection which is explicitly based on ground motion IMs, including SA, duration, and cumulative measures; rather than a focus on implicit parameters (i.e., distance, and pulse or non-pulse classifications) that are conventionally used to heuristically distinguish between the near-fault and far-field records. No ad hoc criteria, in terms of the number of directivity ground motions and their pulse periods, are enforced for selecting pulse-like records. Example applications are presented with different rupture characteristics, source-to-site geometry, and site conditions. It is advocated that the selection of ground motions in the near-fault region based on IM properties alone is preferred to that in which the proportion of pulse-like motions and their pulse periods are specified a priori as strict criteria for ground motion selection. Three methods are presented to propagate the effect of seismic hazard and ground motion selection epistemic uncertainties to seismic performance metrics. These methods differ in their level of rigor considered to propagate the epistemic uncertainty in the conditional distribution of IMs utilised in ground motion selection, selected ground motion ensembles, and the number of nonlinear response history analyses performed to obtain the distribution of engineering demand parameters. These methods are compared for an example site where it is observed that, for seismic demand levels below the collapse limit, epistemic uncertainty in ground motion selection is a smaller uncertainty contributor relative to the uncertainty in the seismic hazard itself. In contrast, uncertainty in ground motion selection process increases the uncertainty in the seismic demand hazard for near-collapse demand levels.
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.
