This study provides an initial examination of source parameter uncertainty in a New Zealand ground motion simulation model, by simulating multiple event realisations with perturbed source parameters. Small magnitude events in Canterbury have been selected for this study due to the small number of source input parameters, the wealth of recorded data, and the lack of appreciable off-fault non-linear effects. Which provides greater opportunity to identify systematic source, path and site effects, required to robustly investigate the causes of uncertainty.
This poster presents the on-going development of a new 3D seismic velocity model of Canterbury, New Zealand. The intention of the model is to provide the 3D crustal structure in the region at multiple length scales for seismic wave propagation simulations, both broadband ground motion and more localized shallow site response analyses.
A preliminary case study assessing the seismic sustainability of two reinforced concrete structures, a frame structure and a wall structure, was conducted to determine which structural system is more seismically sustainable. The two structures were designed to the same standards and were assumed to be located in Christchurch, New Zealand. A component-based probabilistic seismic loss assessment, considering direct losses only, was conducted for two ground motion records, regarded to approximately represent a 1 in 500 year earthquake event and a 1 in 2500 year earthquake event, respectively. It is shown that the wall structure results in lower direct losses than the frame structure in the less severe ground motion scenario. However, in the more severe ground motion scenario, the frame structure results in lower direct losses. Hence, this study demonstrates that which structural system has the lower direct losses depends on the ground motion intensity level.
During the 2011 M7.8 Kaikōura earthquake, ground motions recorded near the epicentre showed a significant spatial variation. The Te Mara farm (WTMC) station, the nearest to the epicentre, recorded 1g and 2.7g of horizontal and vertical peak ground accelerations (PGA), respectively. The nearby Waiu Gorge (WIGC) station recorded a horizontal PGA of 0.8g. Interestingly, however, the Culverden Airlie Farm (CULC) station that was very closely located to WIGC recorded a horizontal PGA of only 0.25g. This poster demonstrates how the local geological condition could have contributed to the spatially variable ground motions observed in the North Canterbury, based on the results of recently conducted geophysical investigations. The surficial geology of this area is dominated by alluvial gravel deposits with traces of silt. A borehole log showed that the thickness of the sediments at WTMC is over 76 metres. Interestingly, the shear wave velocity (Vs) profiles obtained from the three strong motion sites suggest unusually high shear wave velocity of the gravelly sediments. The velocity of sediments and the lack of clear peaks in the horizontal-to-vertical (H/V) spectral ratio at WTMC suggest that the large ground motion observed at this station was likely caused by the proximity of the station to the causative fault itself; the site effect was likely insignificant. Comparisons of H/V spectral ratios and Vs profiles suggest that the sediment thickness is much smaller at WIGC compared with CULC; the high PGA at WIGC was likely influenced by the high-frequency amplification caused by the response of shallow sediments.
8-pages Special devices can be used to minimize structural damage by energy dissipation or seismic isolation. This research considers High Force-to-Volume (HF2V), Symmetric Friction Connection (SFC), Asymmetric Friction Connection (AFC) and Linear-Elastic Isolators (LEI). Device architectures connecting column-to-deck and ground-to-deck are also compared. Bridge columns are assumed to remain elastic. Performance of bridge columns (peak and residual displacement) under 20 probabilistically scaled ground motions is assessed in spectral analysis (0.1-5.0sec) using reduction factors compared to a fixed, no-device case. Energy dissipating devices have minimum column displacement reduction factors when placed between the column and the deck for rigid connection system periods up to ~2.5s. Above that fundamental period, dissipating devices connecting ground-to-deck provide the optimum configuration. Residual displacements obtained when the energy dissipators are placed between the column and the deck are larger than those of the ground to deck case for periods below ~3.7s. Above this 3.7s, frictional dissipators in the column to deck case are more efficient, but HF2V devices connecting ground to deck remain as the best alternative with no residual displacements. The performance curves obtained in this research provide design guidelines for the best device and configuration applicable to a broad range of bridge structures.
This study investigates the uncertainty of simulated earthquake ground motions for smallmagnitude events (Mw 3.5 – 5) in Canterbury, New Zealand. 148 events were simulated with specified uncertainties in: event magnitude, hypocentre location, focal mechanism, high frequency rupture velocity, Brune stress parameter, the site 30-m time-averaged shear wave velocity (Vs30), anelastic attenuation (Q) and high frequency path duration. In order to capture these uncertainties, 25 realisations for each event were generated using the Graves and Pitarka (2015) hybrid broadband simulation approach. Monte-Carlo realisations were drawn from distributions for each uncertainty, to generate a suite of simulation realisations for each event and site. The fit of the multiple simulation realisations to observations were assessed using linear mixed effects regression to generate the systematic source, path and site effects components across all ground motion intensity measure residuals. Findings show that additional uncertainties are required in each of the three source, path, and site components, however the level of output uncertainty is promising considering the input uncertainties included.
Our poster will present on-going QuakeCoRE-founded work on strong motion seismology for Dunedin-Mosgiel area, focusing on ground motion simulations for Dunedin Central Business District (CBD). Source modelling and ground motion simulations are being carried out using the SCEC (Southern California Earthquakes Center) Broad Band simulation Platform (BBP). The platform computes broadband (0-10 Hz) seismograms for earthquakes and was first implemented at the University of Otago in 2016. As large earthquakes has not been experienced in Dunedin in the time of period of instrumental recording, user-specified scenario simulations are of great value. The Akatore Fault, the most active fault in Otago and closest major fault to Dunedin, is the source focused on in the present study. Simulations for various Akatore Fault source scenarios are run and presented. Path and site effects are key components considered in the simulation process. A 1D shear wave velocity profile is required by SCEC BBP, and this is being generated to represent the Akatore-to-CBD path and site within the BBP. A 3D shear velocity model, with high resolution within Dunedin CBD, is being developed in parallel with this study (see Sangster et al. poster). This model will be the basis for developing a 3D shear wave velocity model for greater Dunedin-Mosgiel area for future ground motion simulations, using Canterbury software (currently under development).
Research Report No.2010-03 Ground motion prediction equations (GMPEs) for geometric-mean pseudo-spectral acceleration amplitudes from New Zealand (NZ) earthquakes are developed. A database of 2437 three-component ground motion records is developed by applying stringent quality criteria to the historically recorded events in NZ. Despite the large number of records, the database is deficient in empirical records from large magnitude events recorded at close distances to the fault rupture plane. As a result, the basis for the NZ-specific GMPE development is to examine the applicability of foreign GMPEs for similar tectonic regions and then modify the most applicable GMPEs based on both theoretical and statistically significant empirically-driven arguments. For active shallow crustal events, five different GMPEs are considered. It was found that the McVerry et al. (2006) model, which is the current model upon which seismic design guidelines and site-specific seismic hazard analyses in NZ are based, provided the worst fit to the NZ database, and that the Chiou et al. (2010) (C10) modification of the Chiou and Youngs (2008) model was the most applicable. Discrepancies between the C10 model and the NZ database that were empirically identified and theoretically justified were used to modify the C10 model for: (i) small magnitude scaling; (ii) scaling of short period ground motion from normal faulting events in volcanic crust; (iii) scaling of ground motions on very hard rock sites; (iv) anelastic attenuation in the NZ crust; and (v) consideration of the increased anelastic attenuation in the Taupo Volcanic Zone (TVZ). For subduction slab events, initially three models were considered. It was found that all of the models had some significant biases with respect to applicability for NZ. The Zhao et al. (2006) (Z06) model was selected because of the rigorous database upon which it was developed and modified by: (i) NZ-specific scaling at small magnitudes; (ii) path scaling at large distances; (iii) consideration of the increased TVZ attenuation; and (iv) revision of the standard deviation model. Based on these modifications the developed model showed no bias of the inter- and intra-event residuals as a function of various predictor variables. The standard deviation of the residuals using the revised standard deviation model also indicated that the model has an adequate precision. Three GMPEs were considered for subduction interface events. The Zhao et al. (2006) (Z06) model was the best performing model with only bias exhibited in the site response model, and possible over-prediction of large magnitude events. The Z06 interface model was modified to account for site response and magnitude scaling using the same functional forms as those of the developed active shallow crustal and subduction slab models. The developed model showed no bias of the inter- and intra-event residuals as a function of various predictor variables. The developed GMPEs include specific features as evident in the NZ database; consistent scaling for parameters not well constrained by the NZ database; and pseudo-spectral amplitudes for vibration periods from 0.01 to 10 seconds. Hence, these models represent a significant advance in the state-of-the art for empirical ground motion prediction in NZ.
Deformational properties of soil, in terms of modulus and damping, exert a great influence on seismic response of soil sites. However, these properties for sands containing some portion of fines particles have not been systematically addressed. In addition, simultaneous modelling of the modulus and damping behaviour of soils during cyclic loading is desirable. This study presents an experimental and computational investigation into the deformational properties of sands containing fines content in the context of site response analysis. The experimental investigation is carried on sandy soils sourced from Christchurch, New Zealand using a dynamic triaxial apparatus while the computational aspect is based on the framework of total-stress one-dimensional (1D) cyclic behaviour of soil. The experimental investigation focused on a systematic study on the deformational behaviour of sand with different amounts of fines content (particle diameter ≤ 75µm) under drained conditions. The silty sands were prepared by mixing clean sand with three different percentages of fines content. A series of bender element tests at small-strain range and stress-controlled dynamic triaxial tests at medium to high-strain ranges were conducted on samples of clean sand and silty sand. This allowed measurements of linear and nonlinear deformational properties of the same specimen for a wide strain range. The testing program was designed to quantify the effects of void ratio and fines content on the low-strain stiffness of the silty sand as well as on the nonlinear stress-strain relationship and corresponding shear modulus and damping properties as a function of cyclic shear strains. Shear wave velocity, Vs, and maximum shear modulus, Gmax, of silty sand was shown to be significantly smaller than the respective values for clean sands measured at the same void ratio, e, or same relative density, Dr. However, the test results showed that the difference in the level of nonlinearity between clean sand and silty sands was small. For loose samples prepared at an identical relative density, the behaviour of clean sand was slightly less nonlinear as compared to sandy soils with higher fines content. This difference in the nonlinear behaviour of clean sand and sandy soils was negligible for dense soils. Furthermore, no systematic influence of fines content on the material damping curve was observed for sands with fines content FC = 0 to 30%. In order to normalize the effects of fines on moduli of sands, equivalent granular void ratio, e*, was employed. This was done through quantifying the participation of fines content in the force transfer chain of the sand matrix. As such, a unified framework for modelling of the variability of shear wave velocity, Vs, (or shear modulus, Gmax) with void ratio was achieved for clean sands and sands with fines, irrespective of their fines content. Furthermore, modelling of the cyclic stress-strain behaviour based on this experimental program was investigated. The modelling effort focused on developing a simple constitutive model which simultaneously models the soil modulus and damping relationships with shear strains observed in laboratory tests. The backbone curve of the cyclic model was adopted based on a modified version of Kondner and Zelasko (MKZ) hyperbolic function, with a curvature coefficient, a. In order to simulate the hysteretic cycles, the conventional Masing rules (Pyke 1979) were revised. The parameter n, in the Masing’s criteria was assumed to be a function of material damping, h, measured in the laboratory. As such the modulus and damping produced by the numerical model could match the stress-strain behaviour observed in the laboratory over the course of this study. It was shown that the Masing parameter n, is strain-dependent and generally takes values of n ≤ 2. The model was then verified through element test simulations under different cyclic loadings. It was shown that the model could accurately simulate the modulus and the damping simultaneously. The model was then incorporated within the OpenSees computational platform and was used to scrutinize the effects of damping on one-dimensional seismic site response analysis. For this purpose, several strong motion stations which recorded the Canterbury earthquake sequence were selected. The soil profiles were modelled as semi-infinite horizontally layered deposits overlying a uniform half-space subjected to vertically propagating shear waves. The advantages and limitations of the nonlinear model in terms of simulating soil nonlinearity and associated material damping were further scrutinized. It was shown that generally, the conventional Masing criteria unconservatively may underestimate some response parameters such as spectral accelerations. This was shown to be due to larger hysteretic damping modelled by using conventional Masing criteria. In addition, maximum shear strains within the soil profiles were also computed smaller in comparison to the values calculated by the proposed model. Further analyses were performed to study the simulation of backbone curve beyond the strain ranges addressed in the experimental phase of this study. A key issue that was identified was that relying only on the modulus reduction curves to simulate the stress-strain behaviour of soil may not capture the actual soil strength at larger strains. Hence, strength properties of the soil layer should also be incorporated to accurately simulate the backbone curve.
Following the 22 February 2011, MW 6.2 earthquake located on a fault beneath the Port Hills of Christchurch, fissuring of up to several hundred metres in length was observed in the loess and loess-colluvium of foot-slope positions in north-facing valleys of the Port Hills. The fissuring was observed in all major valleys, occurred at similar low altitudes, showing a contour-parallel orientation and often accompanied by both lateral compression/extension features and spring formation in the valley floor below. Fissuring locations studied in depth included Bowenvale Valley, Hillsborough Valley, Huntlywood Terrace–Lucas Lane, Bridle Path Road, and Maffeys Road–La Costa Lane. Investigations into loess soil, its properties and mannerisms, as well as international examples of its failure were undertaken, including study of the Loess Plateau of China, the Teton Dam, and palaeo-fissuring on Banks Peninsula. These investigations lead to the conclusion that loess has the propensity to fail, often due to the infiltration of water, the presence of which can lead to its instantaneous disaggregation. Literature study and laboratory analysis of Port Hills loess concluded that is has the ability to be stable in steep, sub-vertical escarpments, and often has a sub-vertically jointed internal structure and has a peak shear strength when dry. Values for cohesion, c (kPa) and the internal friction angle, ϕ (degrees) of Port Hills loess were established. The c values for the 40 Rapaki Road, 3 Glenview Terrace loess samples were 13.4 kPa and 19.7 kPa, respectively. The corresponding ϕ values were thought unusually high, at 42.0° and 43.4°.The analysed loess behaved very plastically, with little or no peak strength visible in the plots as the test went almost directly to residual strength. A geophysics resistivity survey showed an area of low resistivity which likely corresponds to a zone of saturated clayey loess/loess colluvium, indicating a high water table in the area. This is consistent with the appearances of local springs which are located towards the northern end of each distinct section of fissure trace and chemical analysis shows that they are sourced from the Port Hills volcanics. Port Hills fissuring may be sub-divided into three categories, Category A, Category B, and Category C, each characterised by distinctive features of the fissures. Category A includes fissures which display evidence of, spring formation, tunnel-gullying, and lateral spreading-like behaviour or quasi-toppling. These fissures are several metres down-slope of the loess-bedrock interface, and are in valleys containing a loess-colluvium fill. Category B fissures are in wider valleys than those in Category A, and the valleys contain estuarine silty sediments which liquefied during the earthquake. Category C fissures occurred at higher elevations than the fissures in the preceding categories, being almost coincident with bedrock outcropping. It is believed that the mechanism responsible for causing the fissuring is a complex combination of three mechanisms: the trampoline effect, bedrock fracturing, and lateral spreading. These three mechanisms can be applied in varying degrees to each of the fissuring sites in categories A, B, and C, in order to provide explanation for the observations made at each. Toppling failure can describe the soil movement as a consequence of the a three causative mechanisms, and provides insight into the movement of the loess. Intra-loess water coursing and tunnel gullying is thought to have encouraged and exacerbated the fissuring, while not being the driving force per se. Incipient landsliding is considered to be the least likely of the possible fissuring interpretations.
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.
This study explicitly investigates uncertainties in physics-based ground motion simulation validation for earthquakes in the Canterbury region. The simulations utilise the Graves and Pitarka (2015) hybrid methodology, with separately quantified parametric uncertainties in the comprehensive physics and simplified physics components of the model. The study is limited to the simulation of 148 small magnitude (Mw 3.5 – 5) earthquakes, with a point source approximation for the source rupture representations, which also enables a focus on a small number of relevant uncertainties. The parametric uncertainties under consideration were selected through sensitivity analysis, and specifically include: magnitude, Brune stress parameter and high frequency rupture velocity. Twenty Monte Carlo realisations were used to sample parameter uncertainties for each of the 148 events. Residuals associated with the following intensity measures: spectral acceleration, peak ground velocity, arias intensity and significant duration, were ascertained. Using these residuals, validation was performed through assessment of systematic biases in site and source terms from mixed-effects regression. Based on the results to date, initial standard deviation recommendations for parameter uncertainties, based on the Canterbury simulations have been obtained. This work ultimately provides an initial step toward explicit incorporation of modelling uncertainty in simulated ground motion predictions for future events, which will improve the use of simulation models in seismic hazard analysis. We plan to subsequently assess uncertainties for larger magnitude events with more complex ruptures, and events across a larger geographic region, as well as uncertainties due to path attenuation, site effects, and more general model epistemic uncertainties.
In this article we utilize grounded theory to explore women’s experiences in the unique construction industry context that followed the 2010 Canterbury (New Zealand) earthquakes. Data were obtained from 36 semi-structured interviews conducted with women working in a variety of occupations in the construction industry. We identify three inter-related categories: capitalizing on opportunity, demonstrating capability and surface tolerance, which together represent a response process that we label ‘deferential tailoring’. The deferential tailoring process explains how women intentionally shape their response to industry conditions through self-regulating behaviors that enables them to successfully seize opportunities and manage gender-related challenges in the working environment. Our findings challenge existing research which suggests that women adopt submissive coping strategies to conform to androcentric norms in the construction industry. Instead, we argue that the process of deferential tailoring can empower women to build positive workplace relationships, enhance career development, and help shift perceptions of the value of their work in the industry.
Following the Mw 6.2 Christchurch Earthquake on 22 February 2011, extensive ground cracking in loessial soils was reported in some areas of the Port Hills, southeast of central Christchurch. This study was undertaken to investigate the mechanisms of earthquake-induced ground damage on the eastern side of the Hillsborough Valley. A zone of extensional cracking up to 40m wide and 600m long was identified along the eastern foot-slope, accompanied by compression features and spring formation at the toe of the slope. An engineering geological and geomorphological model was developed for the eastern Hillsborough Valley that incorporates geotechnical investigation data sourced from the Canterbury Geotechnical Database (CGD), the findings of trenching and seismic refraction surveying carried out for this research, and interpretation of historical aerial photographs. The thickness and extent of a buried peat swamp at the base of the slope was mapped, and found to coincide with significant compression features. Ground cracking was found to have occurred entirely within loess-colluvium and to follow the apices of pre-1920s tunnel-gully fan debris at the southern end of the valley. The ground-cracking on the eastern side of the Hillsborough Valley is interpreted to have formed through tensile failure of the loess-colluvium. Testing was carried out to determine the tensile strength of Port Hills loess colluvium as a function of water content and density, in order to better understand the occurrence and distribution of the observed ground cracking. A comprehensive review of the soil tensile strength testing literature was undertaken, from which a test methodology was developed. Results show remoulded loess-colluvium to possess tensile strength of 7 - 28 kPa across the range of tested moisture contents (10-15%) and dry densities (1650-1900kg/m3). A positive linear relationship was observed between tensile strength and dry density, and a negative linear relationship between moisture content and tensile strength. The observed ground damage and available geotechnical information (inclinometer and piezometer records provided by the Earthquake Commission) were together used to interpret the mechanism(s) of slope movement that occurred in the eastern Hillsborough Valley. The observed ground damage is characteristic of translational movement, but without the development of lateral release scarps, or a basal sliding surface - which was not located during drilling. It is hypothesised that shear displacement has been accommodated by multiple slip surfaces of limited extent within the upper 10m of the slope. Movement has likely occurred within near-saturated colluvial units that have lost strength during earthquake shaking. The eastern Hillsborough Valley is considered to be an ‘incipient translational slide’, as both the patterns of damage and shearing are consistent with the early stages of such slide development. Sliding block analysis was utilised to understand how the eastern Hillsborough Valley may perform in a future large magnitude earthquake. Known cumulative displacements of ~0.3m for eastern Hillsborough Valley during the 2010-2011 Canterbury Earthquake Sequence were compared with modelled slope displacements to back-analyse a lower-bound yield acceleration of 0.2 - 0.25g. Synthetic broadband modelling for future Alpine and Hope Fault earthquakes indicates PGAs of approximately 0.08g for soil sites in the Christchurch area, as such, slope movement is unlikely to be reactivated by an Alpine Fault or Hope Fault earthquake. This does not take into account the possible role of strength loss due to excess pore pressure that may occur during these future events.
This thesis documents the development and demonstration of an assessment method for analysing earthquake-related damage to concrete waste water gravity pipes in Christchurch, New Zealand, following the 2010-2011 Canterbury Earthquake Sequence (CES). The method is intended to be internationally adaptable to assist territorial local authorities with improving lifelines infrastructure disaster impact assessment and improvements in resilience. This is achieved through the provision of high-resolution, localised damage data, which demonstrate earthquake impacts along the pipe length. The insights gained will assist decision making and the prioritisation of resources following earthquake events to quickly and efficiently restore network function and reduce community impacts. The method involved obtaining a selection of 55 reinforced concrete gravity waste water pipes with available Closed-Circuit Television (CCTV) inspection footage filmed before and after the CES. The pipes were assessed by reviewing the recordings, and damage was mapped to the nearest metre along the pipe length using Geographic Information Systems. An established, systematic coding process was used for reporting the nature and severity of the observed damage, and to differentiate between pre-existing and new damage resulting from the CES. The damage items were overlaid with geospatial data such as Light Detection and Ranging (LiDAR)-derived ground deformation data, Liquefaction Resistance Index data and seismic ground motion data (Peak Ground acceleration and Peak Ground Velocity) to identify potential relationships between these parameters and pipe performance. Initial assessment outcomes for the pipe selection revealed that main pipe joints and lateral connections were more vulnerable than the pipe body during a seismic event. Smaller diameter pipes may also be more vulnerable than larger pipes during a seismic event. Obvious differential ground movement resulted in increased local damage observations in many cases, however this was not obvious for all pipes. Pipes with older installation ages exhibited more overall damage prior to a seismic event, which is likely attributable to increased chemical and biological deterioration. However, no evidence was found relating pipe age to performance during a seismic event. No evidence was found linking levels of pre-CES damage in a pipe with subsequent seismic performance, and seismic performance with liquefaction resistance or magnitude of seismic ground motion. The results reported are of limited application due to the small demonstration sample size, but reveal the additional level of detail and insight possible using the method presented in this thesis over existing assessment methods, especially in relation to high resolution variations along the length of the pipe such as localised ground deformations evidenced by LiDAR. The results may be improved by studying a larger and more diverse sample pool, automating data collection and input processes in order to improve efficiency and consider additional input such as pipe dip and cumulative damage over a large distance. The method is dependent on comprehensive and accurate pre-event CCTV assessments and LIDAR data so that post-event data could be compared. It is proposed that local territorial authorities should prioritise acquiring this information as a first important step towards improving the seismic resilience of a gravity waste water pipe network.
This thesis documents the development and demonstration of an assessment method for analysing earthquake-related damage to concrete waste water gravity pipes in Christchurch, New Zealand, following the 2010-2011 Canterbury Earthquake Sequence (CES). The method is intended to be internationally adaptable to assist territorial local authorities with improving lifelines infrastructure disaster impact assessment and improvements in resilience. This is achieved through the provision of high-resolution, localised damage data, which demonstrate earthquake impacts along the pipe length. The insights gained will assist decision making and the prioritisation of resources following earthquake events to quickly and efficiently restore network function and reduce community impacts. The method involved obtaining a selection of 55 reinforced concrete gravity waste water pipes with available Closed-Circuit Television (CCTV) inspection footage filmed before and after the CES. The pipes were assessed by reviewing the recordings, and damage was mapped to the nearest metre along the pipe length using Geographic Information Systems. An established, systematic coding process was used for reporting the nature and severity of the observed damage, and to differentiate between pre-existing and new damage resulting from the CES. The damage items were overlaid with geospatial data such as Light Detection and Ranging (LiDAR)-derived ground deformation data, Liquefaction Resistance Index data and seismic ground motion data (Peak Ground acceleration and Peak Ground Velocity) to identify potential relationships between these parameters and pipe performance. Initial assessment outcomes for the pipe selection revealed that main pipe joints and lateral connections were more vulnerable than the pipe body during a seismic event. Smaller diameter pipes may also be more vulnerable than larger pipes during a seismic event. Obvious differential ground movement resulted in increased local damage observations in many cases, however this was not obvious for all pipes. Pipes with older installation ages exhibited more overall damage prior to a seismic event, which is likely attributable to increased chemical and biological deterioration. However, no evidence was found relating pipe age to performance during a seismic event. No evidence was found linking levels of pre-CES damage in a pipe with subsequent seismic performance, and seismic performance with liquefaction resistance or magnitude of seismic ground motion. The results reported are of limited application due to the small demonstration sample size, but reveal the additional level of detail and insight possible using the method presented in this thesis over existing assessment methods, especially in relation to high resolution variations along the length of the pipe such as localised ground deformations evidenced by LiDAR. The results may be improved by studying a larger and more diverse sample pool, automating data collection and input processes in order to improve efficiency and consider additional input such as pipe dip and cumulative damage over a large distance. The method is dependent on comprehensive and accurate pre-event CCTV assessments and LIDAR data so that post-event data could be compared. It is proposed that local territorial authorities should prioritise acquiring this information as a first important step towards improving the seismic resilience of a gravity waste water pipe network.
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 paper provides a summary of initial research results investigating systematic site effects from the prediction residuals of empirical- and physics-based ground-motion models (GMMs) for small magnitude (i.e., 3.5 ≤ MW ≤ 5) active shallow crustal earthquakes in New Zealand (NZ). Advancing ground-motion predictability through physics-based GMMs is an iterative process and requires addressing fundamental questions like: Is there salient physics which has been overlooked? Which geographic regions have predictions that significantly deviate from observations and why? Which sites exhibit systematic prediction residuals and how can the attributes influencing them be identified? This preliminary study examines these questions by classifying 171 sites from the Canterbury and Wellington regions into four geomorphic categories: basin, basin-edge, hill, and valley, following the categorisation by Nweke et al. (2022). Trends in the site-to-site residuals for each geomorphic category indicate apparent differences between the four categories, with residuals for valley sites illustrating a clear dependence with the inferred fundamental site period. Computed residuals from both empirical- and physics-based GMMs also provided insight into the role of site-specific attributes vs. the different prediction methods, assisting to understand the salient causes of these residuals.
The city of Christchurch and its surrounds experienced widespread damage due to soil liquefaction induced by seismic shaking during the Canterbury earthquake sequence that began in September 2010 with the Mw7.1 Darfield earthquake. Prior to the start of this sequence, the city had a large network of strong motion stations (SMSs) installed, which were able to record a vast database of strong ground motions. This paper uses this database of strong ground motion recordings, observations of liquefaction manifestation at the ground surface, and data from a recently completed extensive geotechnical site investigation program at each SMS to assess a range of liquefaction evaluation procedures at the four SMSs in the Christchurch Central Business District (CBD). In general, the characteristics of the accelerograms recorded at each SMS correlated well with the liquefaction evaluation procedures, with low liquefaction factors of safety predicted at sites with clear liquefaction identifiers in the ground motions. However, at sites that likely liquefied at depth (as indicated by evaluation procedures and/or inferred from the characteristics of the recorded surface accelerograms), the presence of a non-liquefiable crust layer at many of the SMS locations prevented the manifestation of any surface effects. Because of this, there was not a good correlation between surface manifestation and two surface manifestation indices, the Liquefaction Potential Index (LPI) and the Liquefaction Severity Number (LSN).
© 2019, Springer-Verlag GmbH Germany, part of Springer Nature. Prediction of building collapse due to significant seismic motion is a principle objective of earthquake engineers, particularly after a major seismic event when the structure is damaged and decisions may need to be made rapidly concerning the safe occupation of a building or surrounding areas. Traditional model-based pushover analyses are effective, but only if the structural properties are well understood, which is not the case after an event when that information is most useful. This paper combines hysteresis loop analysis (HLA) structural health monitoring (SHM) and incremental dynamic analysis (IDA) methods to identify and then analyse collapse capacity and the probability of collapse for a specific structure, at any time, a range of earthquake excitations to ensure robustness. This nonlinear dynamic analysis enables constant updating of building performance predictions following a given and subsequent earthquake events, which can result in difficult to identify deterioration of structural components and their resulting capacity, all of which is far more difficult using static pushover analysis. The combined methods and analysis provide near real-time updating of the collapse fragility curves as events progress, thus quantifying the change of collapse probability or seismic induced losses very soon after an earthquake for decision-making. Thus, this combination of methods enables a novel, higher-resolution analysis of risk that was not previously available. The methods are not computationally expensive and there is no requirement for a validated numerical model, thus providing a relatively simpler means of assessing collapse probability immediately post-event when such speed can provide better information for critical decision-making. Finally, the results also show a clear need to extend the area of SHM toward creating improved predictive models for analysis of subsequent events, where the Christchurch series of 2010–2011 had significant post-event aftershocks.
As a consequence of the 2010 – 2011 Canterbury earthquake sequence, Christchurch experienced widespread liquefaction, vertical settlement and lateral spreading. These geological processes caused extensive damage to both housing and infrastructure, and increased the need for geotechnical investigation substantially. Cone Penetration Testing (CPT) has become the most common method for liquefaction assessment in Christchurch, and issues have been identified with the soil behaviour type, liquefaction potential and vertical settlement estimates, particularly in the north-western suburbs of Christchurch where soils consist mostly of silts, clayey silts and silty clays. The CPT soil behaviour type often appears to over-estimate the fines content within a soil, while the liquefaction potential and vertical settlement are often calculated higher than those measured after the Canterbury earthquake sequence. To investigate these issues, laboratory work was carried out on three adjacent CPT/borehole pairs from the Groynes Park subdivision in northern Christchurch. Boreholes were logged according to NZGS standards, separated into stratigraphic layers, and laboratory tests were conducted on representative samples. Comparison of these results with the CPT soil behaviour types provided valuable information, where 62% of soils on average were specified by the CPT at the Groynes Park subdivision as finer than what was actually present, 20% of soils on average were specified as coarser than what was actually present, and only 18% of soils on average were correctly classified by the CPT. Hence the CPT soil behaviour type is not accurately describing the stratigraphic profile at the Groynes Park subdivision, and it is understood that this is also the case in much of northwest Christchurch where similar soils are found. The computer software CLiq, by GeoLogismiki, uses assessment parameter constants which are able to be adjusted with each CPT file, in an attempt to make each more accurate. These parameter changes can in some cases substantially alter the results for liquefaction analysis. The sensitivity of the overall assessment method, raising and lowering the water table, lowering the soil behaviour type index, Ic, liquefaction cutoff value, the layer detection option, and the weighting factor option, were analysed by comparison with a set of ‘base settings’. The investigation confirmed that liquefaction analysis results can be very sensitive to the parameters selected, and demonstrated the dependency of the soil behaviour type on the soil behaviour type index, as the tested assessment parameters made very little to no changes to the soil behaviour type plots. The soil behaviour type index, Ic, developed by Robertson and Wride (1998) has been used to define a soil’s behaviour type, which is defined according to a set of numerical boundaries. In addition to this, the liquefaction cutoff point is defined as Ic > 2.6, whereby it is assumed that any soils with an Ic value above this will not liquefy due to clay-like tendencies (Robertson and Wride, 1998). The method has been identified in this thesis as being potentially unsuitable for some areas of Christchurch as it was developed for mostly sandy soils. An alternative methodology involving adjustment of the Robertson and Wride (1998) soil behaviour type boundaries is proposed as follows: Ic < 1.31 – Gravelly sand to dense sand 1.31 < Ic < 1.90 – Sands: clean sand to silty sand 1.90 < Ic < 2.50 – Sand mixtures: silty sand to sandy silt 2.50 < Ic < 3.20 – Silt mixtures: clayey silt to silty clay 3.20 < Ic < 3.60 – Clays: silty clay to clay Ic > 3.60 – Organics soils: peats. When the soil behaviour type boundary changes were applied to 15 test sites throughout Christchurch, 67% showed an improved change of soil behaviour type, while the remaining 33% remained unchanged, because they consisted almost entirely of sand. Within these boundary changes, the liquefaction cutoff point was moved from Ic > 2.6 to Ic > 2.5 and altered the liquefaction potential and vertical settlement to more realistic ii values. This confirmed that the overall soil behaviour type boundary changes appear to solve both the soil behaviour type issues and reduce the overestimation of liquefaction potential and vertical settlement. This thesis acts as a starting point towards researching the issues discussed. In particular, future work which would be useful includes investigation of the CLiq assessment parameter adjustments, and those which would be most suitable for use in clay-rich soils such as those in Christchurch. In particular consideration of how the water table can be better assessed when perched layers of water exist, with the limitation that only one elevation can be entered into CLiq. Additionally, a useful investigation would be a comparison of the known liquefaction and settlements from the Canterbury earthquake sequence with the liquefaction and settlement potentials calculated in CLiq for equivalent shaking conditions. This would enable the difference between the two to be accurately defined, and a suitable adjustment applied. Finally, inconsistencies between the Laser-Sizer and Hydrometer should be investigated, as the Laser-Sizer under-estimated the fines content by up to one third of the Hydrometer values.
Recurrent liquefaction in Christchurch during the 2010-2011 Canterbury earthquake sequence created a wealth of shallow subsurface intrusions with geometries and orientations governed by (1) strong ground motion severity and duration, and (2) intrinsic site characteristics including liquefaction susceptibility, lateral spreading severity, geomorphic setting, host sediment heterogeneity, and anthropogenic soil modifications. We present a suite of case studies that demonstrate how each of these characteristics influenced the geologic expressions of contemporary liquefaction in the shallow subsurface. We compare contemporary features with paleo-features to show how geologic investigations of recurrent liquefaction can provide novel insights into the shaking characteristics of modern and paleo-earthquakes, the influence of geomorphology on liquefaction vulnerability, and the possible controls of anthropogenic activity on the geologic record. We conclude that (a) sites of paleo-liquefaction in the last 1000-2000 years corresponded with most severe liquefaction during the Canterbury earthquake sequence, (b) less vulnerable sites that only liquefied in the strongest and most proximal contemporary earthquakes are unlikely to have liquefied in the last 1000-2000 years or more, (c) proximal strong earthquakes with large vertical accelerations favoured sill formation at some locations, (d) contemporary liquefaction was more severe than paleoliquefaction at all study sites, and (e) stratigraphic records of successive dike formation were more complete at sites with severe lateral spreading, (f) anthropogenic fill suppressed surface liquefaction features and altered subsurface liquefaction architecture.
Liquefaction-induced lateral spreading in large seismic events often results in pervasive and costly damage to engineering structures and lifelines, making it a critical component of engineering design. However, the complex nature of this phenomenon leads to designing for such a hazard extremely challenging and there is a clear for an improved understanding and predicting liquefaction-induced lateral spreading. The 2010-2011 Canterbury (New Zealand) Earthquakes triggered severe liquefaction-induced lateral spreading along the streams and rivers of the Christchurch region, causing extensive damage to roads, bridges, lifelines, and structures in the vicinity. The unfortunate devastation induced from lateral spreading in these events also rendered the rare opportunity to gain an improved understanding of lateral spreading displacements specific to the Christchurch region. As part of this thesis, the method of ground surveying was employed following the 4 September 2010 Darfield (Mw 7.1) and 22 February 2011 Christchurch (Mw 6.2) earthquakes at 126 locations (19 repeated) throughout Christchurch and surrounding suburbs. The method involved measurements and then summation of crack widths along a specific alignment (transect) running approximately perpendicular to the waterway to indicate typically a maximum lateral displacement at the bank and reduction of the magnitude of displacements with distance from the river. Rigorous data processing and comparisons with alternative measurements of lateral spreading were performed to verify results from field observations and validate the method of ground surveying employed, as well as highlight the complex nature of lateral spreading displacements. The welldocumented field data was scrutinized to gain an understanding of typical magnitudes and distribution patterns (distribution of displacement with distance) of lateral spreading observed in the Christchurch area. Maximum displacements ranging from less than 10 cm to over 3.5 m were encountered at the sites surveyed and the area affected by spreading ranged from less than 20 m to over 200 m from the river. Despite the highly non-uniform displacements, four characteristic distribution patterns including large, distributed ground displacements, block-type movements, large and localized ground displacements, and areas of little to no displacements were identified. Available geotechnical, seismic, and topographic data were collated at the ground surveying sites for subsequent analysis of field measurements. Two widely-used empirical models (Zhang et al. (2004), Youd et al. (2002)) were scrutinized and applied to locations in the vicinity of field measurements for comparison with model predictions. The results indicated generally poor correlation (outside a factor of two) with empirical predictions at most locations and further validated the need for an improved, analysis- based method of predicting lateral displacements that considers the many factors involved on a site-specific basis. In addition, the development of appropriate model input parameters for the Youd et al. (2002) model led to a site-specific correlation of soil behavior type index, Ic, and fines content, FC, for sites along the Avon River in Christchurch that matched up well with existing Ic – FC relationships commonly used in current practice. Lastly, a rigorous analysis was performed for 25 selected locations of ground surveying measurements along the Avon River where ground slope conditions are mild (-1 to 2%) and channel heights range from about 2 – 4.5 m. The field data was divided into categories based on the observed distribution pattern of ground displacements including: large and distributed, moderate and distributed, small to negligible, and large and localized. A systematic approach was applied to determine potential critical layers contributing to the observed displacement patterns which led to the development of characteristic profiles for each category considered. The results of these analyses outline an alternative approach to the evaluation of lateral spreading in which a detailed geotechnical analysis is used to identify the potential for large spreading displacements and likely spatial distribution patterns of spreading. Key factors affecting the observed magnitude and distribution of spreading included the thickness of the critical layer, relative density, soil type and layer continuity. It was found that the large and distributed ground displacements were associated with a thick (1.5 – 2.5 m) deposit of loose, fine to silty sand (qc1 ~4-7 MPa, Ic 1.9-2.1, qc1n_cs ~50-70) that was continuous along the bank and with distance from the river. In contrast, small to negligible displacements were characterized by an absence of or relatively thin (< 1 m), discontinuous critical layer. Characteristic features of the moderate and distributed displacements were found to be somewhere between these two extremes. The localized and large displacements showed a characteristic critical layer similar to that observed in the large and distributed sites but that was not continuous and hence leading to the localized zone of displacement. The findings presented in this thesis illustrate the highly complex nature of lateral displacements that cannot be captured in simplified models but require a robust geotechnical analysis similar to that performed for this research.
One of the current challenges in physics-based ground-motion simulations is to refine the modeling of local site effects. These effects require a finer spatial resolution in the material modeling than that generally considered in regional-scale simulations. Because of this, empirical amplification factors are typically applied to capture these unmodeled phenomena. The ergodic nature of this approach suggests that there is room for improvement. In this study, the predictive capability of simulations is evaluated using alternative methods for capturing local site effects. In addition to the conventional empirical approach, two methods are examined that allow for more site-specific information to be incorporated: the square-root impedance method and the 1D time domain site-response analysis. The three approaches are tested using 1000+ observed ground motions from 150+ small-magnitude events (3.5 ≤ Mw ≤ 5.0), recorded at 20 strong-motion stationsin the Canterbury, New Zealand, region. These 20 well-characterized sites represent a wide range of soil conditions, including stiff gravels with Vs30 values greater than 500 m/s, and sand and silt deposits with Vs30 valuesless than 200 m/s. Multiple intensity measures are computed and prediction residuals are partitioned using mixed-effects regression to rigorously assess the relative performance of the different approaches considered. The results indicate that the benefit of using more sophisticated methods is highly dependent on the characteristics of the site. Key site parameters and trends are identified and discussed in light of the assumptions and limitations of each approach.
When the devastating 6.3 magnitude earthquake hit Christchurch, Aotearoa New Zealand, at 12.51pm on 22nd February 2011, the psychological and physical landscape was irrevocably changed. In the days and weeks following the disaster communities were isolated due to failed infrastructure, continuing aftershocks and the extensive search and rescue effort which focussed resources on the central business district. In such moments the resilience of a community is truly tested. This research discusses the role of grassroots community groups in facilitating community resilience during the Christchurch 2010/11 earthquakes and the role of place in doing so. I argue that place specific strategies for urban resilience need to be enacted from a grassroots level while being supported by broader policies and agencies. Using a case study of Project Lyttelton – a group aspiring towards a resilient sustainable future who were caught at the epicentre of the February earthquake – I demonstrate the role of a community group in creating resilience through self-organised place specific action during a disaster. The group provided emotional care, basic facilities and rebuilding assistance to the residents of Lyttelton, proving to be an invaluable asset. These actions are closely linked to the characteristics of social support and social learning that have been identified as important to socio-ecological resilience. In addition this research will seek to understand and explore the nuances of place and identity and its role in shaping resilience to such dis-placing events. Drawing on community narratives of the displacement of place identity, the potential for a progressive sense of place as instigated by local groups will be investigated as an avenue for adaptation by communities at risk of disaster and place destabilisation.
Predicting building collapse due to seismic motion is critical in design and more so after a major event. Damaged structures can appear sound, but collapse under following major events. There can thus be significant risk in decision making after a major seismic event concerning the safe occupation of a building or surrounding areas, versus the unknown impact of unknown major aftershocks. Model-based pushover analyses are effective if the structural properties are well understood, which is not valid post-event when this risk information is most useful. This research combines Hysteresis Loop Analysis (HLA) structural health monitoring (SHM) and Incremental Dynamic Analysis (IDA) methods to determine collapse capacity and probability of collapse for a specific structure, at any time, a range of earthquake excitations to ensure robustness. The nonlinear dynamic analysis method presented enables constant updating of building performance predictions using post-event SHM results. The resulting combined methods provide near real-time updating of collapse fragility curves as events progress, quantifying the change of collapse probability or seismic induced losses for decision-making - a novel, higher resolution risk analysis than previously available. The methods are not computationally expensive and there is no requirement for a validated numerical model. Results show significant potential benefits and a clear evolution of risk. They also show clear need for extending SHM toward creating improved predictive models for analysis of subsequent events, where the Christchurch series of 2010-2011 had significant post-event aftershocks after each main event. Finally, the overall method is generalisable to any typical engineering demand parameter.
The seismic performance of soil profiles with potentially liquefiable deposits is a complex phenomenon that requires a thorough understanding of the soil properties and ground motion characteristics. The limitations of simplified liquefaction assessment methods have prompted an increase in the use of non-linear dynamic analysis methods. Focusing on onedimensional site response of a soil column, this thesis validated a soil constitutive model using in-situ pore pressure measurements and then assessed the influence of input ground motion characteristics on soil column response using traditional and newly developed metrics. Pore pressure recordings during the Canterbury Earthquake Sequence (CES) in New Zealand were used to validate the PM4Sand constitutive model. Soil profile characterization was key to accurate prediction of excess pore pressure response and accounting for any densification during the CES. Response during multiple earthquakes was captured effectively and cross-layer interaction demonstrated the model capability to capture soil response at the system-level. Synthetic and observed ground motions from the Christchurch earthquake were applied to the validated soil column to quantify the performance of synthetic motions. New metrics were developed to facilitate a robust comparison to assess performance. The synthetic input motions demonstrated a slightly larger acceleration and excess pore pressure response compared to the observed input motions. The results suggest that the synthetic motions may accumulate higher excess pore pressure at a faster rate and with fewer number of cycles in the shear response. This research compares validated soil profile subject to spectrally-matched pulse and non-pulse motions, emphasizing the inclusion of pulse motions with distinctive characteristics in ground motion suites for non-linear dynamic analysis. However, spectral matching may lead to undesired alterations in pulse characteristics. Cumulative absolute velocity and significant duration significantly differed between these two groups compared to the other key characteristics and contributed considerably to the liquefaction response. Unlike the non-pulse motions, not all of the pulse motions triggered liquefaction, likely due to their shorter significant duration. Non-pulse motions developed a greater spatial extent of liquefaction triggering in the soil profile and extended to a greater depth
Damage distribution maps from strong earthquakes and recorded data from field experiments have repeatedly shown that the ground surface topography and subsurface stratigraphy play a decisive role in shaping the ground motion characteristics at a site. Published theoretical studies qualitatively agree with observations from past seismic events and experiments; quantitatively, however, they systematically underestimate the absolute level of topographic amplification up to an order of magnitude or more in some cases. We have hypothesized in previous work that this discrepancy stems from idealizations of the geometry, material properties, and incident motion characteristics that most theoretical studies make. In this study, we perform numerical simulations of seismic wave propagation in heterogeneous media with arbitrary ground surface geometry, and compare results with high quality field recordings from a site with strong surface topography. Our goal is to explore whether high-fidelity simulations and realistic numerical models can – contrary to theoretical models – capture quantitatively the frequency and amplitude characteristics of topographic effects. For validation, we use field data from a linear array of nine portable seismometers that we deployed on Mount Pleasant and Heathcote Valley, Christchurch, New Zealand, and we compute empirical standard spectral ratios (SSR) and single-station horizontal-to-vertical spectral ratios (HVSR). The instruments recorded ambient vibrations and remote earthquakes for a period of two months (March-April 2017). We next perform two-dimensional wave propagation simulations using the explicit finite difference code FLAC. We construct our numerical model using a high-resolution (8m) Digital Elevation Map (DEM) available for the site, an estimated subsurface stratigraphy consistent with the geomorphology of the site, and soil properties estimated from in-situ and non-destructive tests. We subject the model to in-plane and out-of-plane incident motions that span a broadband frequency range (0.1-20Hz). Numerical and empirical spectral ratios from our blind prediction are found in very good quantitative agreement for stations on the slope of Mount Pleasant and on the surface of Heathcote Valley, across a wide range of frequencies that reveal the role of topography, soil amplification and basin edge focusing on the distribution of ground surface motion.
The purpose of this thesis is to evaluate the seismic response of the UC Physics Building based on recorded ground motions during the Canterbury earthquakes, and to use the recorded response to evaluate the efficacy of various conventional structural analysis modelling assumptions. The recorded instrument data is examined and analysed to determine how the UC Physics Building performed during the earthquake-induced ground motions. Ten of the largest earthquake events from the 2010-11 Canterbury earthquake sequence are selected in order to understand the seismic response under various levels of demand. Peak response amplitude values are found which characterise the demand from each event. Spectral analysis techniques are utilised to find the natural periods of the structure in each orthogonal direction. Significant torsional and rocking responses are also identified from the recorded ground motions. In addition, the observed building response is used to scrutinise the adequacy of NZ design code prescriptions for fundamental period, response spectra, floor acceleration and effective member stiffness. The efficacy of conventional numerical modelling assumptions for representing the UC Physics Building are examined using the observed building response. The numerical models comprise of the following: a one dimensional multi degree of freedom model, a two dimensional model along each axis of the building and a three dimensional model. Both moderate and strong ground motion records are used to examine the response and subsequently clarify the importance of linear and non-linear responses and the inclusion of base flexibility. The effects of soil-structure interaction are found to be significant in the transverse direction but not the longitudinal direction. Non-linear models predict minor in-elastic behaviour in both directions during the 4 September 2010 Mw 7.1 Darfield earthquake. The observed torsional response is found to be accurately captured by the three dimensional model by considering the interaction between the UC Physics Building and the adjacent structure. With the inclusion of adequate numerical modelling assumptions, the structural response is able to be predicted to within 10% for the majority of the earthquake events considered.
