Search

found 108 results

Research papers, University of Canterbury Library

This paper presents the ongoing development of a new 3D seismic velocity model of Canterbury, New Zealand. The model explicitly represents the Canterbury sedimentary basin, and other significant geologic horizons, which are expected to have important implications on observed ground motions. The model utilizes numerous sources of data, including 3D regional tomography with a variable-depth inferred Moho, seismic reflection survey lines, geotechnical boreholes and well logs, spectral analysis of surface waves, and CPT logs which provide velocity constraints over their respective ranges of application. The model provides P- and S-wave velocity and density (i.e. Vp, Vs and p) over a grid of input points, and is presently being utilized in broadband ground motion simulations of the 2010-2011 Canterbury earthquakes. Comparison of simulated ground motions with those observed in the 2010-2011 Canterbury earthquakes will help provide a better understanding of the salient physical processes which characterized the unique set of strong ground motions recorded in this sequence of earthquake events.

Research papers, Victoria University of Wellington

The Mѡ=7.1 Darfield (Canterbury) earthquake struck on 4 September 2010, approximately 45 km west of Christchurch, New Zealand. It revealed a previously unknown fault (the Greendale fault) and caused billions of dollars of damage due to high peak ground velocities and extensive liquefaction. It also triggered the Mw=6.3 Christchurch earthquake on 22 February 2011, which caused further damage and the loss of 185 lives. The objective of this research was to determine the relationship between stress and seismic properties in a seismically active region using manually-picked P and S wave arrival times from the aftershock sequence between 8 September 2010-13 January 2011 to estimate shear-wave splitting (SWS) parameters, VP =VS-ratios, anisotropy (delay-time tomography), focal mechanisms, and tectonic stress on the Canterbury plains. The maximum horizontal stress direction was highly consistent in the plains, with an average value of SHmax=116 18 . However, the estimates showed variation in SHmax near the fault, with one estimate rotating by as much as 30° counter-clockwise. This suggests heterogeneity of stress at the fault, though the cause remains unclear. Orientations of the principal stresses predominantly indicate a strike-slip regime, but there are possible thrust regimes to the west and north/east of the fault. The SWS fast directions (ø) on the plains show alignment with SHmax at the majority of stations, indicating stress controlled anisotropy. However, structural effects appear more dominant in the neighbouring regions of the Southern Alps and Banks Peninsula.

Research papers, University of Canterbury Library

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.

Audio, Radio New Zealand

The Temple for Christchurch is an enormous work, representing the seismic wave of the February 22 earthquake. Made mostly from salvaged wood it's taken a team of volunteers two years to build. But it will meet a fiery end on Saturday night when it is deliberately set alight. Justin Gregory meets the designer and his volunteers with two weeks left until burn day. But when he gets there, none are to be found.

Images, UC QuakeStudies

A photograph of two 'All Righties' waving in one of the Christchurch Hospital childrens' wards during All Right?'s 'Outrageous Burst of All Right: Superhero Surprise'. This event involved people dressed as superheroes surprising children at Christchurch Hospital with gifts and entertainment. All Right? posted the photograph to their Facebook Timeline on 2 July 2014 at 12:43pm with the caption, "The All Righties preparing to sneak in...".

Images, UC QuakeStudies

A scanned copy of a photograph of an optical diffraction pattern produced by the He-Ne Continuous Gas Laser used in David Lockwood's MSc research at the University of Canterbury. David explains that the photograph shows "a typical optical diffraction pattern obtained from the original red laser beam arising from a grating structure formed by the alignment of a colloid under the influence of a travelling sound wave".

Research papers, University of Canterbury Library

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

Research papers, University of Canterbury Library

This paper presents an overview of the soil profile characteristics at strong motion station (SMS) locations in the Christchurch Central Business District (CBD) based on recently completed geotechnical site investigations. Given the variability of Christchurch soils, detailed investigations were needed in close vicinity to each SMS. In this regard, CPT, SPT and borehole data, and shear wave velocity (Vs) profiles from surface wave dispersion data in close vicinity to the SMSs have been used to develop detailed representative soil profiles at each site and to determine site classes according to the New Zealand standard NZS1170.5. A disparity between the NZS1170.5 site classes based on Vs and SPT N60 investigation techniques is highlighted, and additional studies are needed to harmonize site classification based on these techniques. The short period mode of vibration of soft deposits above gravels, which are found throughout Christchurch, are compared to the long period mode of vibration of the entire soil profile to bedrock. These two distinct modes of vibration require further investigation to determine their impact on the site response. According to current American and European approaches to seismic site classification, all SMSs were classified as problematic soil sites due to the presence of liquefiable strata, soils which are not directly accounted for by the NZS1170.5 approach.

Images, UC QuakeStudies

A scanned copy of a black and white photograph of the He-Ne Continuous Gas Laser used in David Lockwood's MSc research at the University of Canterbury. David says, "The experiments in the Thesis were aimed at examining the effects of sound waves on colloids, which are very small particles that are found in clays for example (for more details see page 181 of http://www.cap.ca/PiC-PaC/static/downloads/1efdc1f3784b85c1a1b33e396b91ee8aef2072c5.pdf ). In fact, I was looking at nanoparticles of matter, which preceded the emergence of nanotechnology as a field of research by more than two decades. I needed a suitable light source to probe the alignment of the nanoparticles in the ultrasonic sound field. At that time the laser had just been invented and I realized that this was the ideal light source for my experiment. I then proceeded to construct the first home-built laser in New Zealand. This laser - a He-Ne continuous-wave gas laser - operated at 632.8 nm (in the red). This laser, over a metre long, is shown in the black-and-white photo. There were lot of problems to be overcome, but eventually, with invaluable help from Dr. Tom Seed (my MSc supervisor) and Dr. Rod Syme, I had it working".

Research papers, University of Canterbury Library

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

Research papers, University of Canterbury Library

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 papers, University of Canterbury Library

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

Research papers, University of Canterbury Library

he strong motion station at Heathcote Valley School (HVSC) recorded unusually high peak ground accelerations (2.21g vertical and 1.41g horizontal) during the February 2011 Christchurch earthquake. Ground motions recorded at HVSC in numerous other events also exhibited consistently higher intensities compared with nearby strong motion stations. We investigated the underlying causes of such high intensity ground motions at HVSC by means of 2D dynamic finite element analyses, using recorded ground motions during the 2010-2011 Canterbury earthquake sequence. The model takes advantage of a LiDAR-based digital elevation model (DEM) to account for the surface topography, while the geometry and dynamic properties of the surficial soils are characterized by seismic cone penetration tests (sCPT) and Multi-Channel Analyses of Surface Waves (MASW). Comparisons of simulated and recorded ground motions suggests that our model performs well for distant events, while for near-field events, ground motions recorded at the adopted reference station at Lyttelton Port are not reasonable input motions for the simulation. The simulations suggest that Rayleigh waves generated at the inclined interface of the surficial colluvium and underlying volcanic rock strongly affect the ground motions recorded at HVSC, in particular, being the dominant contributor to the recorded vertical motions.

Images, Alexander Turnbull Library

Inside the Christchurch City Council building someone wonders what that noise is. Outside a gigantic wave of 'public discontent' rushes towards them. Context: Relates to anger and a protest of 4000 people at what they see as poor performance by the City Council in dealing with earthquake issues. The final straw was the Council's decision on the advice of its consultants to award the Council Chief Executive Tony Marryatt a $68,000 14% pay rise mid December 2011. Quantity: 1 digital cartoon(s).

Research papers, Lincoln University

Creativity that is driven by a need for physical or economic survival, which disasters are likely to inspire, raises the question of whether such creativity fits with conventional theories and perspectives of creativity. In this paper we use the opportunity afforded by the 2010-2013 Christchurch, New Zealand earthquakes to follow and assess the creative practices and responses of a number of groups and individuals. We use in-depth interviews to tease out motivations and read these against a range of theoretical propositions about creativity. In particular, we focus on the construct of “elite panic” and the degree to which this appeared to be evident in the Christchurch earthquakes context. Bureaucratic attempts to control or limit creativity were present but they did not produce a completely blanket dampening effect. Certain individuals and groups seemed to be pre-equipped to navigate or ignore potential blocks to creativity. We argue, using Geir Kaufmann’s novelty-creativity matrix and aspects of Teresa Amabile’s and Michael G. Pratt’s revised componential theory of creativity that a special form of disaster creativity does exist.

Images, Alexander Turnbull Library

MP Jim Anderton who is standing for Mayor of Christchurch in the coming local body elections (9th October) stands up to his chest in a crack in the road waving a sign that reads 'Jim Anderton for mayor' and yelling "I'm still here!' Refers to the very high profile that his rival the standing mayor Bob Parker has enjoyed in the aftermath of the Christchurch earthquake of 4th September 2010. Jim Anderton has rather faded out of the limelight by comparison. Quantity: 1 digital cartoon(s).

Research papers, University of Canterbury Library

This paper provides an overview of the salient aspects of the dense array of ground motions observed in the 4 September 2010 Darfield and 22 February 2011 Christchurch earthquakes. Particular attention is given to inferred physical reasons for the observed ground motions, which include: (i) source features such as forward directivity effects; (ii) The effects of the Canterbury Plains sedimentary basin on basin-generated surface waves, and waveguide effects through the region; and (iii) the importance of local site response as evidenced by observations of large long period amplification and liquefaction. The significance of vertical ground motion intensity is also examined.

Research papers, University of Canterbury Library

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.

Research papers, University of Canterbury Library

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.

Research papers, University of Canterbury Library

Heathcote Valley school strong motion station (HVSC) consistently recorded ground motions with higher intensities than nearby stations during the 2010-2011 Canterbury earthquakes. For example, as shown in Figure 1, for the 22 February 2011 Christchurch earthquake, peak ground acceleration at HVSC reached 1.4 g (horizontal) and 2 g (vertical), the largest ever recorded in New Zealand. Strong amplification of ground motions is expected at Heathcote Valley due to: 1) the high impedance contrast at the soil-rock interface, and 2) the interference of incident and surface waves within the valley. However, both conventional empirical ground motion prediction equations (GMPE) and the physics-based large scale ground motions simulations (with empirical site response) are ineffective in predicting such amplification due to their respective inherent limitations.

Research papers, University of Canterbury Library

Beach ridge stratigraphy can provide an important record of both sustained coastal progradation and responses to events such as extreme storms, as well as evidence of earthquake induced sediment pulses. This study is a stratigraphic investigation of the late Holocene mixed sand gravel (MSG) beach ridge plain on the Canterbury coast, New Zealand. The subsurface was imaged along a 370 m shore-normal transect using 100 and 200 MHz ground penetrating radar (GPR) antennae, and cored to sample sediment textures. Results show that, seaward of a back-barrier lagoon, the Pegasus Bay beach ridge plain prograded almost uniformly, under conditions of relatively stable sea level. Nearshore sediment supply appears to have created a sustained sediment surplus, perhaps as a result of post-seismic sediment pulses, resulting in a flat, morphologically featureless beach ridge plain. Evidence of a high magnitude storm provides an exception, with an estimated event return period in excess of 100 years. Evidence from the GPR sequence combined with modern process observations from MSG beaches indicates that a paleo storm initially created a washover fan into the back-barrier lagoon, with a large amount of sediment simultaneously moved off the beach face into the nearshore. This erosion event resulted in a topographic depression still evident today. In the subsequent recovery period, sediment was reworked by swash onto the beach as a sequence of berm deposit laminations, creating an elevated beach ridge that also has a modern-day topographic signature. As sediment supply returned to normal, and under conditions of falling sea level, a beach ridge progradation sequence accumulated seaward of the storm feature out to the modern-day beach as a large flat, uniform progradation plain. This study highlights the importance of extreme storm events and earthquake pulses on MSG coastlines in triggering high volume beach ridge formation during the subsequent recovery period.

Research papers, University of Canterbury Library

This study uses 44 high quality liquefaction case histories taken from 22 locations affected by the 2010-2011 Canterbury earthquake sequence to evaluate four commonly used CPT-VS correlations (i.e., Robertson, 2009; Hegazy and Mayne, 2006; Andrus et al., 2007; McGann et al., 2015b). Co-located CPT soundings and VS profiles, developed from surface wave testing, were obtained at 22 locations and case histories were developed for the Mw 7.1, 4 September 2010 Darfield and Mw 6.2, 22 February 2011 Christchurch earthquakes. The CPT soundings are used to generate VS profiles using each of four CPT-VS correlations. These correlated VS profiles are used to estimate the factor of safety against liquefaction using the Kayen et al. (2013) VS-based simplified liquefaction evaluation procedure. An error index is used to quantify the predictive capabilities of these correlations in relation to the observations of liquefaction (or the lack thereof). Additionally, the error indices from the CPT-correlated VS profiles are compared to those obtained using: (1) the Kayen et al. (2013) procedure with surface wave-derived VS profiles, and (2) the Idriss and Boulanger (2008) CPT-based liquefaction evaluation procedure. Based on the error indices, the evaluation procedures based on direct measurements of either CPT or VS provided more accurate liquefaction triggering estimates than those obtained from any of the CPT-VS correlations. However, the performance of the CPT-VS correlations varied, with the Robertson (2009) and Hegazy and Mayne (2006) correlations performing relatively poorly for the Christchurch soils and the Andrus et al. (2007) and McGann et al. (2015b) correlations performing better. The McGann et al. (2015b) correlation had the lowest error indices of the CPT-VS correlations tested, however, none of the CPT-VS correlations provided accurate enough VS predictions to be used for the evaluation of liquefaction triggering using the VS-based liquefaction evaluation procedures.

Research papers, University of Canterbury Library

This poster presents work to date on ground motion simulation validation and inversion for the Canterbury, New Zealand region. Recent developments have focused on the collection of different earthquake sources and the verification of the SPECFEM3D software package in forward and inverse simulations. SPECFEM3D is an open source software package which simulates seismic wave propagation and performs adjoint tomography based upon the spectral-element method. Figure 2: Fence diagrams of shear wave velocities highlighting the salient features of the (a) 1D Canterbury velocity model, and (b) 3D Canterbury velocity model. Figure 5: Seismic sources and strong motion stations in the South Island of New Zealand, and corresponding ray paths of observed ground motions. Figure 3: Domain used for the 19th October 2010 Mw 4.8 case study event including the location of the seismic source and strong motion stations. By understanding the predictive and inversion capabilities of SPECFEM3D, the current 3D Canterbury Velocity Model can be iteratively improved to better predict the observed ground motions. This is achieved by minimizing the misfit between observed and simulated ground motions using the built-in optimization algorithm. Figure 1 shows the Canterbury Velocity Model domain considered including the locations of small-to-moderate Mw events [3-4.5], strong motion stations, and ray paths of observed ground motions. The area covered by the ray paths essentially indicates the area of the model which will be most affected by the waveform inversion. The seismic sources used in the ground motion simulations are centroid moment tensor solutions obtained from GeoNet. All earthquake ruptures are modelled as point sources with a Gaussian source time function. The minimum Mw limit is enforced to ensure good signal-to-noise ratio and well constrained source parameters. The maximum Mw limit is enforced to ensure the point source approximation is valid and to minimize off-fault nonlinear effects.

Research papers, University of Canterbury Library

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.

Research papers, University of Canterbury Library

This work investigates the possibility of developing a non-contact, non-line of sight sensor to measure interstorey drift through simulation and experimental validation. • The method uses frequency-modulated continuous wave (FMCW) radar to measure displacement. This method is commonly in use in a number of modern applications, including aircraft altimeters and automotive parking sensors. • The technique avoids numerous problems found in contemporary structural health monitoring methods, namely integral drift errors and structural modification requirements. • The smallest achievable detection error in displacement was found to be as low as 0.26%, through simulated against the displacement response of a single degree of freedom structure subject to ground motion excitation. • This was verified during experimentation, when a corner-style reflector was placed on a shake table running ground motion data taken from the 4th September 2010 earthquake in Christchurch. These results confirmed the conclusions drawn from simulation.

Images, UC QuakeStudies

A tribute taped to a window of a house on Tasman Place. The tribute reads, "Our red zoned house. When we bought you years ago, you looked a bit tired and sad, but overall you weren't too bad. We spruced you up with paint and love and asked for a blessing from above. The years went by, family and friends celebrations under your roof, your 'veggie' garden gave us kai. We felt safe within your wall, then one dreadful September night, the shaking earth made you fall. You tried with all your groaning might to keep us from harm. Because you were strongly built we held onto the door, while a wave of terror buckled the floor and outside the garden flooded with silt. Now you are near the end, sunken walls and windows bend. We say goodbye today and let you go, Our spirit and heart feels low. You are more than just mortar and brick. For us you were a gift, a safe haven where we once lived".

Research papers, University of Canterbury Library

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.

Research papers, University of Canterbury Library

Well-validated liquefaction constitutive models are increasingly important as non-linear time history analyses become relatively more common in industry for key projects. Previous validation efforts of PM4Sand, a plasticity model specifically for liquefaction, have generally focused on centrifuge tests; however, pore pressure transducers installed at several free-field sites during the Canterbury Earthquake Sequence (CES) in Christchurch, New Zealand provide a relatively unique dataset to validate against. This study presents effective stress site response analyses performed in the finite difference software FLAC to examine the capability of PM4Sand to capture the generation of excess pore pressures during earthquakes. The characterization of the subsurface is primarily based on extensive cone penetration tests (CPT) carried out in Christchurch. Correlations based on penetration resistances are used to estimate soil parameters, such as relative density and shear wave velocity, which affect liquefaction behaviour. The resulting free-field FLAC model is used to estimate time histories of excess pore pressure, which are compared with records during several earthquakes in the CES to assess the suitability of PM4Sand.

Images, UC QuakeStudies

A motion-blurred photograph of houses, with the Port Hills in the background. The photographer comments, "This I hope gives you a feel of what it feels like in an earthquake. When you spend your whole life thinking that you and your home are built on solid ground, it can be quite a shock when you find it is not. You can feel the house shaking like a dog with a toy, rising up violently underneath you or the most gentle form which is when the ground moves gently like a wave moving under a rowing boat. It is not just the movement, you often get a rumbling sound which can precede a violent shake or can result in no movement at all. This means that some vehicles can sound like the rumbling initially and in the early days would get your heart racing. Another form of stress is when big excavators as heavy as a tank move as you can feel the ground shake from streets away, but you do not always hear the engine. For most of us the problem when the shaking starts, is wondering if this is the start of an extremely violent earthquake or will it peter out".

Research papers, University of Canterbury Library

1. Background and Objectives This poster presents results from ground motion simulations of small-to-moderate magnitude (3.5≤Mw≤5.0) earthquake events in the Canterbury, New Zealand region using the Graves and Pitarka (2010,2015) methodology. Subsequent investigation of systematic ground motion effects highlights the prediction bias in the simulations which are also benchmarked against empirical ground motion models (e.g. Bradley (2013)). In this study, 144 earthquake ruptures, modelled as point sources, are considered with 1924 quality-assured ground motions recorded across 45 strong motion stations throughout the Canterbury region, as shown in Figure 1. The majority of sources are Mw≥4.0 and have centroid depth (CD) 10km or shallower. Earthquake source descriptions were obtained from the GeoNet New Zealand earthquake catalogue. The ground motion simulations were performed within a computational domain of 140km x 120km x 46km with a finite difference grid spacing of 0.1km. The low-frequency (LF) simulations utilize the 3D Canterbury Velocity Model while the high-frequency (HF) simulations utilize a generic regional 1D velocity model. In the LF simulations, a minimum shear wave velocity of 500m/s is enforced, yielding a maximum frequency of 1.0Hz.