Despite the relatively low seismicity, a large earthquake in the Waikato region is expected to have a high impact, when the fourth-largest regional population and economy and the high density critical infrastructure systems in this region are considered. Furthermore, Waikato has a deep soft sedimentary basin, which increases the regional seismic hazard due to trapping and amplification of seismic waves and generation of localized surface waves within the basin. This phenomenon is known as the “Basin Effect”, and has been attributed to the increased damage in several historic earthquakes, including the 2010-2011 Canterbury earthquakes. In order to quantitatively model the basin response and improve the understanding of regional seismic hazard, geophysical methods will be used to develop shear wave velocity profiles across the Waikato basin. Active surface wave methods involve the deployment of linear arrays of geophones to record the surface waves generated by a sledge hammer. Passive surface wave methods involve the deployment of two-dimensional seismometer arrays to record ambient vibrations. At each site, the planned testing includes one active test and two to four passive arrays. The obtained data are processed to develop dispersion curves, which describe surface wave propagation velocity as a function of frequency (or wavelength). Dispersion curves are then inverted using the Geopsy software package to develop a suite of shear wave velocity profiles. Currently, more than ten sites in Waikato are under consideration for this project. This poster presents the preliminary results from the two sites that have been tested. The shear wave velocity profiles from all sites will be used to produce a 3D velocity model for the Waikato basin, a part of QuakeCoRE flagship programme 1.
This paper summarizes the development of a high-resolution surficial shear wave velocity model based on the combination of the large high-spatial-density database of cone penetration test (CPT) logs in and around Christchurch, New Zealand and a recently-developed Christchurch-specific empirical correlation between soil shear wave velocity and CPT. This near-surface shear wave velocity model has applications for site characterization efforts via the development of maps of time-averaged shear wave velocities over specific depths, as well as use in site response analysis and ground motion simulation.
This paper summarizes the development of a region-wide surficial shear wave velocity model based on the combination of the large high-spatial-density database of cone penetration test (CPT) logs in and around Christchurch, New Zealand and a recently-developed Christchurch-specific empirical correlation between soil shear wave velocity and CPT. The ongoing development of this near-surface shear wave velocity model has applications for site characterization efforts via the development of maps of time-averaged shear wave velocities over specific depths, and the identification of regional similarities and differences in soil shear stiffness.
Bulk rock strength is greatly dependent on fracture density, so that reductions in rock strength associated with faulting and fracturing should be reflected by reduced shear coupling and hence S-wave velocity. This study is carried out along the Canterbury rangefront and in Otago. Both lie within the broader plate boundary deformation zone in the South Island of New Zealand. Therefore built structures are often, , located in areas where there are undetected or poorly defined faults with associated rock strength reduction. Where structures are sited near to, or across, such faults or fault-zones, they may sustain both shaking and ground deformation damage during an earthquake. Within this zone, management of seismic hazards needs to be based on accurate identification of the potential fault damage zone including the likely width of off-plane deformation. Lateral S-wave velocity variability provides one method of imaging and locating damage zones and off-plane deformation. This research demonstrates the utility of Multi-Channel Analysis of Surface Waves (MASW) to aid land-use planning in such fault-prone settings. Fundamentally, MASW uses surface wave dispersive characteristics to model a near surface profile of S-wave velocity variability as a proxy for bulk rock strength. The technique can aid fault-zone planning not only by locating and defining the extent of fault-zones, but also by defining within-zone variability that is readily correlated with measurable rock properties applicable to both foundation design and the distribution of surface deformation. The calibration sites presented here have well defined field relationships and known fault-zone exposure close to potential MASW survey sites. They were selected to represent a range of progressively softer lithologies from intact and fractured Torlesse Group basement hard rock (Dalethorpe) through softer Tertiary cover sediments (Boby’s Creek) and Quaternary gravels. This facilitated initial calibration of fracture intensity at a high-velocity-contrast site followed by exploration of the limits of shear zone resolution at lower velocity contrasts. Site models were constructed in AutoCAD in order to demonstrate spatial correlations between S-wave velocity and fault zone features. Site geology was incorporated in the models, along with geomorphology, river profiles, scanline locations and crosshole velocity measurement locations. Spatial data were recorded using a total-station survey. The interpreted MASW survey results are presented as two dimensional snapshot cross-sections of the three dimensional calibration-site models. These show strong correlations between MASW survey velocities and site geology, geomorphology, fluvial profiles and geotechnical parameters and observations. Correlations are particularly pronounced where high velocity contrasts exist, whilst weaker correlations are demonstrated in softer lithologies. Geomorphic correlations suggest that off-plane deformation can be imaged and interpreted in the presence of suitable topographic survey data. A promising new approach to in situ and laboratory soft-rock material and mass characterisation is also presented using a Ramset nail gun. Geotechnical investigations typically involve outcrop and laboratory scale determination of rock mass and material properties such as fracture density and unconfined compressive strength (UCS). This multi-scale approach is espoused by this study, with geotechnical and S-wave velocity data presented at multiple scales, from survey scale sonic velocity measurements, through outcrop scale scanline and crosshole sonic velocity measurements to laboratory scale property determination and sonic velocity measurements. S-wave velocities invariably increased with decreasing scale. These scaling relationships and strategies for dealing with them are investigated and presented. Finally, the MASW technique is applied to a concealed fault on the Taieri Ridge in Macraes Flat, Central Otago. Here, high velocity Otago Schist is faulted against low velocity sheared Tertiary and Quaternary sediments. This site highlights the structural sensitivity of the technique by apparently constraining the location of the principal fault, which had been ambiguous after standard processing of the seismic reflection data. Processing of the Taieri Ridge dataset has further led to the proposal of a novel surface wave imaging technique termed Swept Frequency Imaging (SFI). This inchoate technique apparently images the detailed structure of the fault-zone, and is in agreement with the conventionally-determined fault location and an existing partial trench. Overall, the results are promising and are expected to be supported by further trenching in the near future.
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).
Deep shear wave velocity (Vs) profiles (>400 m) were developed at 14 sites throughout Christchurch, New Zealand using surface wave methods. This paper focuses on the inversion of surface wave data collected at one of these sites, Hagley Park. This site is located on the deep soils of the Canterbury Plains, which consist of alluvial gravels inter-bedded with estuarine and marine sands, silts, clays and peats. Consequently, significant velocity contrasts exist at the interface between geologic formations. In order to develop realistic velocity models in this complex geologic environment, a-priori geotechnical and geologic data were used to identify the boundaries between geologic formations. This information aided in developing the layering for the inversion parameters. Moreover, empirical reference Vs profiles based on material type and confining pressure were used to develop realistic Vs ranges for each layer. Both the a-priori layering information and the reference Vs curves proved to be instrumental in generating realistic velocity models that account for the complex inter-bedded geology in the Canterbury Plains.
This presentation summarizes the development of high-resolution surficial soil velocity models in the Canterbury, New Zealand basin. Shallow (<30m) shear wave velocities were primarily computed based on a combination of a large database of over 15,000 cone penetration test (CPT) logs in and around Christchurch, and a recently-developed Christchurch-specific empirical correlation between soil shear wave velocity and CPT. Large active-source testing at 22 locations and ambient-wavefield surface wave and H/V testing at over 80 locations were utilized in combination with 1700 water well logs to constrain the inter-bedded stratigraphy and velocity of Quaternary sediments up to depths of several hundred meters. Finally, seismic reflection profiles and the ambient-wavefield surface wave data provide constraint on velocities from several hundred meters to several kilometres. At all depths, the high resolution data illustrates the complexity of the soil conditions in the region, and the developed 3D models are presently being used in broadband ground motion simulations to further interpret the observed strong ground motions in the 2010-2011 Canterbury earthquake sequence.
A 3D high-resolution model of the geologic structure and associated seismic velocities in the Canterbury, New Zealand region is developed utilising data from depthconverted 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, developed using geostatistical Kriging, explicitly represents the significant and regionally recognisable geologic surfaces that mark the boundaries between geologic units with distinct lithology and age. The model is examined in the form of both geologic surface elevation contour maps as well as vertical cross sections of shear wave velocity, with the most prominent features being the Banks Peninsula Miocene-Pliocene volcanic edifice, and the Pegasus and Rakaia late Mesozoic-Neogene sedimentary basins. The adequacy of the modelled geologic surfaces is assessed through a residual analysis of point constraints used in the Kriging and qualitative comparisons with previous geologic models of subsets of the region. Seismic velocities for the lithological units between the geologic surfaces have also been derived, thus providing the necessary information for a Canterbury velocity model (CantVM) for use in physics-based seismic wave propagation. The developed model also has application for the determination of depths to specified shear wave velocities for use in empirical ground motion modelling, which is explicitly discussed via an example.
Measurement of basement seismic resonance frequencies can elucidate shallow velocity structure, an important factor in earthquake hazard estimation. Ambient noise cross correlation, which is well-suited to studying shallow earth structure, is commonly used to analyze fundamental-mode Rayleigh waves and, increasingly, Love waves. Here we show via multicomponent ambient noise cross correlation that the basement resonance frequency in the Canterbury region of New Zealand can be straightforwardly determined based on the horizontal to vertical amplitude ratio (H/V ratio) of the first higher-mode Rayleigh waves. At periods of 1-3 s, the first higher-mode is evident on the radial-radial cross-correlation functions but almost absent in the vertical-vertical cross-correlation functions, implying longitudinal motion and a high H/V ratio. A one-dimensional regional velocity model incorporating a ~ 1.5 km-thick sedimentary layer fits both the observed H/V ratio and Rayleigh wave group velocity. Similar analysis may enable resonance characteristics of other sedimentary basins to be determined. © 2013. American Geophysical Union. All Rights Reserved.
During many years the analysis of some geophysical results of Charles Darwin was being carried out in Department. Darwin has connected almost 200 years ago results of catastrophic earthquakes with vertical movement of a surface of the Earth. Usually this movement less horizontal movement and its influence on destruction of cities is not considered. Earthquake hazard assessment studies were focused usually on the horizontal ground motion. Effects of the strong vertical motion were not, practically, discussed. The margins of safety against gravity-induced static vertical forces in constructed buildings usually provide adequate resistance to dynamic forces induced by the vertical acceleration during an earthquake. However, the earthquake in Christchurch is an example of the vertical seismic shock . The earthquake magnitude was rather small - nearby 6.3. However, the result was catastrophic. The same took place in 1835. It allowed to Darwin to formulate a few great ideas. Charles Darwin has explained qualitatively results of an interaction of huge seismic waves with volcanoes and the nature of volcanism and seismicity of our planet. These important data of Charles Darwin became very actual recently. It is possible to tell also the same about tsunami and extreme ocean waves described by Charles Darwin. Therefore this data were analyzed using modern mechanics, mathematics and physics in Department. In particular, the theory of catastrophic waves was developed based on Darwin's data. The theory tried to explain occurrence, evolution and distribution the catastrophic waves in various natural systems, since atoms, oceans, surfaces of the Earth and up to the very early Universe. Some results of the research were published in prestigious magazines. Later they were presented in two books devoted to Charles Darwin's anniversary (2009). Last from them was published in Russian (2011). We give here key ideas of this research which is a part of interdisciplinary researches of Department. Some ideas are discussed. Not less important purpose is very short historical review of some researches of Darwin. In particular, we underline Darwin' priority in the formulation of the bases of Dynamics Earth.
This report summarizes the development of a region-wide surficial soil shear wave velocity (Vs ) model based on the unique combination of a large high-spatial-density database of cone penetration test (CPT) logs in the greater Christchurch urban area (> 15, 000 logs as of 1 February 2014) and the Christchurch-specific empirical correlation between soil Vs and CPT data developed by McGann et al. [1, 2]. This model has applications for site characterization efforts via maps of time-averaged Vs over specific depths (e.g. Vs30, Vs10), and for numerical modeling efforts via the identification of typical Vs profiles for different regions and soil behaviour types within Christchurch. In addition, the Vs model can be used to constrain the near-surface velocities for the 3D seismic velocity model of the Canterbury basin [3] currently being developed for the purpose of broadband ground motion simulation. The general development of these region-wide near-surface Vs models includes the following general phases, with each discussed in separate chapters of this report. • An evaluation of the available CPT dataset for suitability, and the definition of other datasets and assumptions necessary to characterize the surficial sediments of the region to 30 m depth. • The development of time-averaged shear wave velocity (Vsz) surfaces for the Christchurch area from the adopted CPT dataset (and supplementary data/assumptions) using spatial interpolation. The Vsz surfaces are used to explore the characteristics of the near-surface soils in the regions and are shown to correspond well with known features of the local geology, the historical ecosystems of the area, and observations made following the 2010- 2011 Canterbury earthquakes. • A detailed analysis of the Vs profiles in eight subregions of Christchurch is performed to assess the variablity in the soil profiles for regions with similar Vsz values and to assess Vsz as a predictive metric for local site response. It is shown that the distrubution of soil shear wave velocity in the Christchurch regions is highly variable both spatially (horizontally) and with depth (vertically) due to the varied geological histories for different parts of the area, and the highly stratified nature of the nearsurface deposits. This variability is not considered to be greatly significant in terms of current simplified site classification systems; based on computed Vs30 values, all considered regions can be categorized as NEHRP sites class D (180 < Vs < 360 m/s) or E (Vs < 180 m/s), however, detailed analysis of the shear wave velocity profiles in different subregions of Christchurch show that the expected surficial site response can vary quite a bit across the region despite the relative similarity in Vs30
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.
The abundance of cone penetration test (CPT) data from subsurface explorations in Christchurch and the surrounding areas provides a useful source of information for a characterization of the near surface shear wave velocity ( ) profile for the region. A portion of the investigations were conducted using seismic CPT, enabling the comparison of measured shear wave velocity with CPT data, and subsequently the evaluation of existing CPT- correlations for applicability to Canterbury-specific soils. The existing correlations are shown to be biased, generally over-predicting the observed with depth, thus demonstrating the need for a Canterbury-specific CPT- correlation.
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.
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.
Recent field investigations were carried out to define the shear wave velocity (VS) profile and site periods across the Canterbury region, supplementing earlier efforts in urban Christchurch. Active source surface wave testing, ambient wave field (passive) and H/V spectral ratio methods were used to characterise the soil profile in the region. H/V spectral ratio peaks indicate site periods in the range of 5-7 seconds across much of the Canterbury Plains, broadly consistent with those based on a 1D velocity model for the region. Site periods decrease rapidly in the vicinity of the Canterbury foothills and the Banks Peninsula outcrops. In Christchurch, the Riccarton Gravels result in a significant mode of vibration that has a much shorter period than the site period of the entire soil column down to basement rock.
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.
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.
This poster provides a summary of the development of a 3D shallow (z<40m) shear wave velocity (Vs) model for the urban Christchurch, New Zealand region. The model is based on a recently developed Christchurch-specific empirical correlation between Vs and cone penetration test (CPT) data (McGann et al. 2014a,b) and the large high-density database of CPT logs in the greater Christchurch urban area (> 15,000 logs as of 01/01/2014). In particular, the 3D model provides shear wave velocities for the surficial Springston Formation, Christchurch Formation, and Riccarton gravel layers which generally comprise the upper 40m in the Christchurch urban area. Point-estimates are provided on a 200m-by- 200m grid from which interpolation to other locations can be performed. This model has applications for future site characterization and numerical modeling efforts via maps of timeaveraged Vs over specific depths (e.g. Vs30, Vs10) and via the identification of typical Vs profiles for different regions and soil behaviour types within Christchurch. In addition, the Vs model can be used to constrain the near-surface velocities for the 3D seismic velocity model of the Canterbury basin (Lee et al. 2014) currently being developed for the purpose of broadband ground motion simulation.
his poster presents the ongoing development of a 3D Canterbury seismic velocity model which will be used in physics-based hybrid broadband ground motion simulation of the 2010-2011 Canterbury earthquakes. Velocity models must sufficiently represent critical aspects of the crustal structure over multiple length scales which will influence the results of the simulations. As a result, numerous sources of data are utilized in order to provide adequate resolution where necessary. Figure 2: (a) Seismic reflection line showing P-wave velocities and significant geologic horizons (Barnes et al. 2011), and (b) Shear wave profiles at 10 locations (Stokoe et al. 2013). Figure 4: Cross sections of the current version of the Canterbury velocity model to depths of 10km as shown in Figure 1: (a) at a constant latitude value of -43.6˚, and (b) at a constant longitude value of 172.64˚. 3. Ground Surface and Geologic Horizon Models Figure 3: (a) Ground surface model derived from numerous available digital elevation models, and (b) Base of the Quaternary sediments derived from structural contours and seismic reflection line elevations. The Canterbury region has a unique and complex geology which likely has a significant impact on strong ground motions, in particular the deep and loose deposits of the Canterbury basin. The Canterbury basin has several implications on seismic wave phenomena such as long period ground motion amplification and wave guide effects. Using a realistic 3D seismic velocity model in physics-based ground motion simulation will implicitly account for such effects and the resultant simulated ground motions can be studied to gain a fundamental understanding of the salient ground motion phenomena which occurred during the Canterbury earthquakes, and the potential for repeat occurrences in the Canterbury region. Figure 1 shows the current model domain as a rectangular area between Lat=[-43.2˚,-44.0˚], and Lon=[171.5˚,173.0˚]. This essentially spans the area between the foot of the Southern Alps in the North West to Banks Peninsula in the East. Currently the model extends to a depth of 50km below sea level.
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.
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.
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.
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.
The overarching goal of this dissertation is to improve predictive capabilities of geotechnical seismic site response analyses by incorporating additional salient physical phenomena that influence site effects. Specifically, multidimensional wave-propagation effects that are neglected in conventional 1D site response analyses are incorporated by: (1) combining results of 3D regional-scale simulations with 1D nonlinear wave-propagation site response analysis, and (2) modelling soil heterogeneity in 2D site response analyses using spatially-correlated random fields to perturb soil properties. A method to combine results from 3D hybrid physics-based ground motion simulations with site-specific nonlinear site response analyses was developed. The 3D simulations capture 3D ground motion phenomena on a regional scale, while the 1D nonlinear site response, which is informed by detailed site-specific soil characterization data, can capture site effects more rigorously. Simulations of 11 moderate-to-large earthquakes from the 2010-2011 Canterbury Earthquake Sequence (CES) at 20 strong motion stations (SMS) were used to validate simulations with observed ground motions. The predictions were compared to those from an empirically-based ground motion model (GMM), and from 3D simulations with simplified VS30- based site effects modelling. By comparing all predictions to observations at seismic recording stations, it was found that the 3D physics-based simulations can predict ground motions with comparable bias and uncertainty as the GMM, albeit, with significantly lower bias at long periods. Additionally, the explicit modelling of nonlinear site-response improves predictions significantly compared to the simplified VS30-based approach for soft-soil or atypical sites that exhibit exceptionally strong site effects. A method to account for the spatial variability of soils and wave scattering in 2D site response analyses was developed and validated against a database of vertical array sites in California. The inputs required to run the 2D analyses are nominally the same as those required for 1D analyses (except for spatial correlation parameters), enabling easier adoption in practice. The first step was to create the platform and workflow, and to perform a sensitivity study involving 5,400 2D model realizations to investigate the influence of random field input parameters on wave scattering and site response. Boundary conditions were carefully assessed to understand their effect on the modelled response and select appropriate assumptions for use on a 2D model with lateral heterogeneities. Multiple ground-motion intensity measures (IMs) were analyzed to quantify the influence from random field input parameters and boundary conditions. It was found that this method is capable of scattering seismic waves and creating spatially-varying ground motions at the ground surface. The redistribution of ground-motion energy across wider frequency bands, and the scattering attenuation of high-frequency waves in 2D analyses, resemble features observed in empirical transfer functions (ETFs) computed in other studies. The developed 2D method was subsequently extended to more complicated multi-layer soil profiles and applied to a database of 21 vertical array sites in California to test its appropriate- ness for future predictions. Again, different boundary condition and input motion assumptions were explored to extend the method to the in-situ conditions of a vertical array (with a sensor embedded in the soil). ETFs were compared to theoretical transfer functions (TTFs) from conventional 1D analyses and 2D analyses with heterogeneity. Residuals of transfer-function- based IMs, and IMs of surface ground motions, were also used as validation metrics. The spatial variability of transfer-function-based IMs was estimated from 2D models and compared to the event-to-event variability from ETFs. This method was found capable of significantly improving predictions of median ETF amplification factors, especially for sites that display higher event-to-event variability. For sites that are well represented by 1D methods, the 2D approach can underpredict amplification factors at higher modes, suggesting that the level of heterogeneity may be over-represented by the 2D random field models used in this study.
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.
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.
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.
This report provides an initial overview and gap analysis of the multi-hazards interactions that might affect fluvial and pluvial flooding (FPF) hazard in the Ōpāwaho Heathcote catchment. As per the terms of reference, this report focuses on a one-way analysis of the potential effects of multi-hazards on FPF hazard, as opposed to a more complex multi-way analysis of interactions between all hazards. We examined the relationship between FPF hazard and hazards associated with the phenomena of tsunamis; coastal erosion; coastal inundation; groundwater; earthquakes; and mass movements. Tsunamis: Modelling research indicates the worst-case tsunami scenarios potentially affecting the Ōpāwaho Heathcote catchment are far field. Under low probability, high impact tsunami scenarios waves could travel into Pegasus Bay and the Avon-Heathcote Estuary Ihutai, reaching the mouth and lower reaches of the Heathcote catchment and river, potentially inundating and eroding shorelines in sub-catchments 1 to 5, and temporarily blocking fluvial drainage more extensively. Any flooding infrastructure or management actions implemented in the area of tsunami inundation would ideally be resilient to tsunami-induced inundation and erosion. Model results currently available are a first estimate of potential tsunami inundation under contemporary sea and land level conditions. In terms of future large tsunami events, these models likely underestimate effects in riverside sub-catchments, as well as effects under future sea level, shoreline and other conditions. Also of significance when considering different FPF management structures, it is important to be mindful that certain types of flood structures can ‘trap’ inundating water coming from ocean directions, leading to longer flood durations and salinization issues. Coastal erosion: Model predictions indicate that sub-catchments 1 to 3 could potentially be affected by coastal erosion by the timescale of 2065, with sub-catchments 1-6 predicted to be potentially affected by coastal erosion by the time scale of 2115. In addition, the predicted open coast effects of this hazard should not be ignored since any significant changes in the New Brighton Spit open coast would affect erosion rates and exposure of the landward estuary margins, including the shorelines of the Ōpāwaho Heathcote catchment. Any FPF flooding infrastructure or management activities planned for the potentially affected sub-catchments needs to recognise the possibility of coastal erosion, and to have a planned response to the predicted potential shoreline translation. Coastal inundation: Model predictions indicate coastal inundation hazards could potentially affect sub-catchments 1 to 8 by 2065, with a greater area and depth of inundation possible for these same sub-catchments by 2115. Low-lying areas of the Ōpāwaho Heathcote catchment and river channel that discharge into the estuary are highly vulnerable to coastal inundation since elevated ocean and estuary water levels can block the drainage of inland systems, compounding FPF hazards. Coastal inundation can overwhelm stormwater and other drainage network components, and render river dredging options ineffective at best, flood enhancing at worst. A distinction can be made between coastal inundation and coastal erosion in terms of the potential impacts on affected land and assets, including flood infrastructure, and the implications for acceptance, adaptation, mitigation, and/or modification options. That is, responding to inundation could include structural and/or building elevation solutions, since unlike erosion, inundation does not necessarily mean the loss of land. Groundwater: Groundwater levels are of significant but variable concern when examining flooding hazards and management options in the Ōpāwaho Heathcote catchment due to variability in soils, topographies, elevations and proximities to riverine and estuarine surface waterbodies. Much of the Canterbury Plains part of the Ōpāwaho Heathcote catchment has a water table that is at a median depth of <1m from the surface (with actual depth below surface varying seasonally, inter-annually and during extreme meteorological events), though the water table depth rapidly shifts to >6m below the surface in the upper Plains part of the catchment (sub-catchments 13 to 15). Parts of Waltham/Linwood (sub-catchments 5 & 6) and Spreydon (sub-catchment 10) have extensive areas with a particularly high water table, as do sub-catchments 18, 19 and 20 south of the river. In all of the sub-catchments where groundwater depth below surface is shallow, it is necessary to be mindful of cascading effects on liquefaction hazard during earthquake events, including earthquake-induced drainage network and stormwater infrastructure damage. In turn, subsidence induced by liquefaction and other earthquake processes during the CES directly affected groundwater depth below surface across large parts of the central Ōpāwaho Heathcote catchment. The estuary margin of the catchment also faces increasing future challenges with sea level rise, which has the potential to elevate groundwater levels in these areas, compounding existing liquefaction and other earthquake associated multi-hazards. Any increases in subsurface runoff due to drainage system, development or climate changes are also of concern for the loess covered hill slopes due to the potential to enhance mass movement hazards. Earthquakes: Earthquake associated vertical ground displacement and liquefaction have historically affected, or are in future predicted to affect, all Ōpāwaho Heathcote sub-catchments. During the CES, these phenomena induced a significant cascades of changes in the city’s drainage systems, including: extensive vertical displacement and liquefaction induced damage to stormwater ‘greyware’, reducing functionality of the stormwater system; damage to the wastewater system which temporarily lowered groundwater levels and increased stormwater drainage via the wastewater network on the one hand, creating a pollution multi-hazard for FPF on the other hand; liquefaction and vertical displacement induced river channel changes affected drainage capacities; subsidence induced losses in soakage and infiltration capacities; changes occurred in topographic drainage conductivity; estuary subsidence (mainly around the Ōtākaro Avon rivermouth) increased both FPF and coastal inundation hazards; estuary bed uplift (severe around the Ōpāwaho Heathcote margins), reduced tidal prisms and increased bed friction, producing an overall reduction the waterbody’s capacity to efficiently flush catchment floodwaters to sea; and changes in estuarine and riverine ecosystems. All such possible effects need to be considered when evaluating present and future capacities of the Ōpāwaho Heathcote catchment FPF management systems. These phenomena are particularly of concern in the Ōpāwaho Heathcote catchment since stormwater networks must deal with constraints imposed by stream and river channels (past and present), estuarine shorelines and complex hill topography. Mass movements: Mass movements are primarily a risk in the Port Hills areas of the Ōpāwaho Heathcote catchment (sub-catchments 1, 2, 7, 9, 11, 16, 21), though there are one or two small but susceptible areas on the banks of the Ōpāwaho Heathcote River. Mass movements in the form of rockfalls and debris flows occurred on the Port Hills during the CES, resulting in building damage, fatalities and evacuations. Evidence has also been found of earthquake-triggered tunnel gully collapsesin all Port Hill Valleys. Follow-on effects of these mass movements are likely to occur in major future FPF and other hazard events. Of note, elevated groundwater levels, coastal inundation, earthquakes (including liquefaction and other effects), and mass movement exhibit the most extensive levels of multi-hazard interaction with FPF hazard. Further, all of the analysed multi-hazard interactions except earthquakes were found to consistently produce increases in the FPF hazard. The implications of these analyses are that multihazard interactions generally enhance the FPF hazard in the Ōpāwaho Heathcote catchment. Hence, management plans which exclude adjustments for multi-hazard interactions are likely to underestimate the FPF hazard in numerous different ways. In conclusion, although only a one-way analysis of the potential effects of selected multi-hazards on FPF hazard, this review highlights that the Ōpāwaho Heathcote catchment is an inherently multi- hazard prone environment. The implications of the interactions and process linkages revealed in this report are that several significant multi-hazard influences and process interactions must be taken into account in order to design a resilient FPF hazard management strategy.
Documenting earthquake-induced ground deformation is significant to assess the characteristics of past and contemporary earthquakes and provide insight into seismic hazard. This study uses airborne light detection and ranging (LiDAR) and conducts multi-disciplinary field techniques to document the surface rupture morphology and evaluate the paleoseismicity and seismic hazard parameters of the Hurunui segment of the Hope Fault in the northern South Island of New Zealand. It also documents and evaluates seismically induced features and ground motion characteristics of the 2010 Darfield and 2011 Christchurch earthquakes in the Port Hills, south of Christchurch. These two studies are linked in that they investigate the near-field coseismic features of large (Mw ~7.1) earthquakes in New Zealand and produce data for evaluating seismic hazards of future earthquakes. In the northern South Island of New Zealand, the Australian-Pacific plate boundary is characterised by strike-slip deformation across the Marlborough Fault System (MFS). The ENE-striking Hope Fault (length: ~230 km) is the youngest and southernmost fault in the MFS, and the second fastest slipping fault in New Zealand. The Hope Fault is a major source of seismic hazard in New Zealand and has ruptured (in-part) historically in the Mw 7.1 1888 Amuri earthquake. In the west, the Hurunui segment of the Hope Fault is covered by beech forest. Hence, its seismic hazard parameters and paleoearthquake chronology were poorly constrained and it was unknown whether the 1888 earthquake ruptured this segment or not and if so, to what extent. Utilising LiDAR and field data, a 29 km-long section of the Hurunui segment of the Hope Fault is mapped. LiDAR-mapping clearly reveals the principal slip zone (PSZ) of the fault and a suite of previously unrecognised structures that form the fault deformation zone (FDZ). FDZ width measurements from 415 locations reveal a spatially-variable, active FDZ up to ~500 m wide with an average width of 200 m. Kinematic analysis of the fault structures shows that the Hurunui segment strikes between 070° and 075° and is optimally oriented for dextral strike-slip within the regional stress field. This implies that the wide FDZ observed is unlikely to result from large-scale fault mis-orientation with respect to regional stresses. The analysis of FDZ width indicates that it increases with increased hanging wall topography and increased topographic relief suggesting that along-strike topographic perturbations to fault geometry and stress states increase fault zone complexity and width. FDZ width also increases where the tips of adjacent PSZ strands locally vary in strike, and where the thickness of alluvial deposits overlying bedrock increases. LiDAR- and photogrammetrically-derived topographic mapping indicates that the boundary between the Hurunui and Hope River segments is characterised by a ~850-m-wide right stepover and a 9º-14° fault bend. Paleoseismic trenching at Hope Shelter site reveals that 6 earthquakes occurred at A.D. 1888, 1740-1840, 1479-1623, 819-1092, 439-551, and 373- 419. These rupture events have a mean recurrence interval of ~298 ± 88 yr and inter-event times ranging from 98 to 595 yrs. The variation in the inter-event times is explained by (1) coalescing rupture overlap from the adjacent Hope River segment on to the Hurunui segment at the study site, (2) temporal clustering of large earthquakes on the Hurunui segment, and/or (3) ‘missing’ rupture events. It appears that the first two options are more plausible to explain the earthquake chronologies and rupture behaviour on the Hurunui segment, given the detailed nature of the geologic and chronologic investigations. This study provides first evidence for coseismic multi-segment ruptures on the Hope Fault by identifying a rupture length of 44-70 km for the 1888 earthquake, which was not confined to the Hope River segment (primary source for the 1888 earthquake). LiDAR data is also used to identify and measure dextral displacements and scarp heights from the PSZ and structures within the FDZ along the Hurunui segment. Reconstruction of large dextrally-offset geomorphic features shows that the vertical component of slip accounts for only ~1% of the horizontal displacements and confirms that the fault is predominantly strike-slip. A strong correlation exists between the dextral displacements and elevations of geomorphic features suggesting the possibility of age correlation between the geomorphic features. A mean single event displacement (SED) of 3.6 ± 0.7 m is determined from interpretation of sets of dextral displacements of ≤ 25 m. Using the available surface age data and the cumulative dextral displacements from Matagouri Flat, McKenzie Fan, Macs Knob and Hope River sites, and the mean SED, a mean slip rate of 12.2 ± 2.4 mm/yr, and a mean recurrence interval of ~320 ± 120 yr, and a potential earthquake magnitude of Mw 7.2 are determined for the Hurunui segment. This study suggests that the fault slip rate has been constant over the last ~15000 yr. Strong ground motions from the 2010 Darfield (Canterbury) earthquake displaced boulders and caused ground damage on some ridge crests in the Port Hills. However, the 2011 Christchurch earthquake neither displaced boulders nor caused ground damage at the same ridge crests. Documentation of locations (~400 m a.s.l.), lateral displacements (8-970 cm), displacement direction (250° ± 20°) of displaced boulders, in addition to their hosting socket geometries (< 1 cm to 50 cm depth), the orientation of the ridges (000°-015°) indicate that boulders have been displaced in the direction of instrumentally recorded transient peak ground horizontal displacements nearby and that the seismic waves have been amplified at the study sites. The co-existence of displaced and non-displaced boulders at proximal sites suggests small-scale ground motion variability and/or varying boulder-ground dynamic interactions relating to shallow phenomena such as variability in soil depth, bedrock fracture density and/or microtopography on the bedrock-soil interface. Shorter shaking duration of the 2011 Christchurch event, differing frequency contents and different source characteristics were all factors that may have contributed to generating circumstances less favourable to boulder displacement in this earthquake. Investigating seismically induced features, fault behaviour, site effects on the rupture behaviour, and site response to the seismic waves provides insights into fault rupture hazards.
