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 Lake Coleridge Rock Avalanche Deposits (LCRADs) are located on Ryton Station in the middle Rakaia Valley, approximately 80 km west of Christchurch. Torlesse Supergroup greywacke is the basement material and has been significantly influenced by both active tectonics and glaciation. Both glacial and post-glacial processes have produced large volumes of material which blanket the bedrock on slopes and in the valley floors. The LCRADs were part of a regional study of rock avalanches by WHITEHOUSE (1981, 1983) and WHITEHOUSE and GRIFFITHS (1983), and a single rock avalanche event was recognised with a weathering rind age of 120 years B.P. that was later modified to 150 ± 40 years B.P. The present study has refined details of both the age and the sequence of events at the site, by identifying three separate rock avalanche deposits (termed the LCRA1, LCRA2 and LCRA3 deposits), which are all sourced from near the summit of Carriage Drive. The LCRA1 deposit is lobate in shape and had an estimated original deposit volume of 12.5 x 10⁶ m³, although erosion by the Ryton River has reduced the present day debris volume to 5.1 x 10⁶ m³. An optically stimulated luminescence date taken from sandy loess immediately beneath the LCRA1 deposit provided a maximum age for the rock avalanche event of 9,720 ± 750 years B.P., which is believed to be realistic given that this is shortly after the retreat of Acheron 3 ice from this part of the valley. Emplacement of rock avalanche material into an ancestral Ryton riverbed created a natural dam with a ~17 M m³ lake upstream. The river is thought to have created a natural spillway over the dam structure at ~557 m (a.s.l), and to have existed for a number of years before any significant downcutting occurred. Although a triggering mechanism for the LCRA1 deposit was poorly constrained, it is thought that stress rebound after glacial ice removal may have initiated failure. Due to the event occurring c.10,000 years ago, there was a lack of definition for a possible earthquake trigger, though the possibility is obvious. The LCRA₂ event had an original deposit volume of 0.66 x 10⁶ m³, and was constrained to the low-lying area adjacent to the Ryton River that had been created by river erosion of the LCRA1 deposit. Further erosion by the Ryton River has reduced the deposit volume to 0.4 x 10⁶ m³. A radiocarbon date from a piece of mānuka found within the LCRA2 deposit provided an age of 668 ± 36 years B.P., and this is thought to reliably date the event. The LCRA2 event also dammed the Ryton River, and the preservation of dam-break outwash terraces downstream from the deposit provides clear evidence of rapid dam erosion and flooding after overtopping, and breaching by the Ryton River. Based on the mean annual flow of the Ryton River, the LCRA2 lake would have taken approximately two weeks to fill assuming that there were no preferred breach paths and the material was relatively impermeable. The LCRA2 event is thought to have been coseismic with a fault rupture along the western segment of the PPAFZ, which has been dated at 600 ± 100 years B.P. by SMITH (2003). The small LCRA3 event was not able to be dated, but it is believed to have failed shortly after the LCRA2 event and it may in fact be a lag deposit of the second rock avalanche event possibly triggered by an aftershock. The deposit is only visible at one locality within the cliffs that line the Ryton River, and its lack of geomorphic expression is attributed to it occurring closely after the LCRA2 event, while the Ryton River was still dammed from the second rock avalanche event. A wedge-block of some 35,000 m³ of source material for a future rock avalanche was identified at the summit of Carriage Drive. The dilation of the rock mass, combined with unfavourably oriented sub-vertical bedding in the Torlesse Supergroup bedrock, has allowed toppling-style failure on both of the main ridge lines around the source area for the LCRADs. In the event of a future rock avalanche occurring within the Ryton riverbed an emergency response plan has been developed to provide a staged response, especially in relation to the camping ground located at the mouth of the Ryton River. A long-term management plan has also been developed for mitigation measures for the Ryton riverbed and adjacent floodplain areas downstream of a future rock avalanche at the LCRAD site.
<strong>Sea level rise is one consequence of Earth’s changing climate. Century-long tide gauge records show that global-mean sea-level rise reached 11-16 cm during the twentieth century at a mean rate of 1.2 mm/y. Today, the average rate of global-mean sea-level rise is higher at 3-4 mm/y and is expected to increase in the future. This represents a hazard to low elevation coastal zones worldwide. Yet, before global sea level projections can be used to characterise future coastal flood hazard at a local scale, the effects of tectonics (and other processes) that drive vertical land motion (VLM) must be considered. VLM is defined as the vertical velocity (uplift or subsidence) of the solid surface with respect to the centre of Earth. In this study, new VLM maps are generated over coastal strips in New Zealand, using Sentinel-1 InSAR and GNSS data.</strong>In New Zealand, measuring VLM using InSAR on naturally vegetated or agricultural land is difficult due to signal decorrelation. Along the rural Bay of Plenty coastal strip, I use a persistent-scatterer approach to generate a VLM map from both east-looking ascending and west-looking descending Sentinel-1 data between 2015-2021. Using time-series data over the same time period from a dense network of 20 GNSS sensors, I tie InSAR-derived line-of-sight velocity to the 2014 ITRF reference frame. I test two different methods for measuring VLM and compare the results against GNSS vertical velocity along the Bay of Plenty coast. Best results are achieved by first removing the interpolated horizontal GNSS velocity field from each of the InSAR datasets, before averaging the two VLM estimates. Measured VLM is between -3 and 3 mm/y, with negative values (subsidence) occurring within the low-lying Rangitāiki Plain and Ōpōtiki valley, and uplift across the elevated region west of Matatā.This thesis integrates geomorphological, geological, and historical levelling VLM records with modern satellite datasets to assess VLM across timescales ranging from 10 to 100,000 years at Matatā. Uplift rate has been variable through time, with average uplift over the last 300,000 years of 1 mm/y, 4.5 mm/y since 1720 years, 2 mm/y between 1950-1978, and 10 mm/y between 2004-2011. Previous modelling has shown that the best fit to the 2004-2011 rapid uplift rates is an inflating magmatic source at ~10 km depth beneath Matatā. To reconcile all data, I present a VLM model that consists of short-lived periods (7 years) of rapid uplift (10 mm/y), separated by longer periods (30 years) of lower background uplift (3 mm/y). The episodic nature of VLM at Matatā likely reflects short-lived periods of magmatic intrusion. Episodic VLM characterised by large rates of uplift (10 mm/y) has been seen at Taupō volcano, and other volcanic centers globally. It has been 12 years since the end of the last intrusion episode; this modelling suggest one may expect to observe increased uplift rates at Matatā in the coming decades. Densely populated urban coastal strips are most at risk from the effects of relative sea-level rise. At the same time, anthropogenic activities associated with urbanization, such as groundwater withdrawal, and land reclamation can lead to local land subsidence (LLS), further exacerbating the risk to urban infrastructure. LLS refers to subsidence relative to nearby land area assumed to be stable. In this thesis, I create the first high-resolution (10 m) maps of LLS at six urban coastal strips in New Zealand, with a combined length of 285 km, using Sentinel-1 InSAR data between 2018-2021. This analysis reveals 89% of urban coastal strips are subsiding at rates of -0.5 mm/y or greater, and 11% is subsiding at higher rates of -3.0 mm/y or greater. On average, subsidence is -0.6 to -2.9 mm/y higher at the coastal strip, compared to inland areas occupied by GNSS stations. This analysis also documents highly-localised hotspots of LLS, with subsidence rates of up to -15 mm/y. In Christchurch, rapid and localised subsidence (-8 mm/y) is observed within coastal suburbs New Brighton and Southshore. In most cities, the highest subsidence rates occur on land reclaimed in the early-late twentieth century, and in areas built on Holocene sediment. Time-series analysis of LLS at sites of reclaimed land shows both linear and non-linear rates of deformation over time periods of up to 6-8 years. This thesis highlights the variable exposure to relative sea-level rise of New Zealand coastal strips, and demonstrates that in many cases current rates of VLM should be expected to continue for the next few decades.
Deformational properties of soil, in terms of modulus and damping, exert a great influence on seismic response of soil sites. However, these properties for sands containing some portion of fines particles have not been systematically addressed. In addition, simultaneous modelling of the modulus and damping behaviour of soils during cyclic loading is desirable. This study presents an experimental and computational investigation into the deformational properties of sands containing fines content in the context of site response analysis. The experimental investigation is carried on sandy soils sourced from Christchurch, New Zealand using a dynamic triaxial apparatus while the computational aspect is based on the framework of total-stress one-dimensional (1D) cyclic behaviour of soil. The experimental investigation focused on a systematic study on the deformational behaviour of sand with different amounts of fines content (particle diameter ≤ 75µm) under drained conditions. The silty sands were prepared by mixing clean sand with three different percentages of fines content. A series of bender element tests at small-strain range and stress-controlled dynamic triaxial tests at medium to high-strain ranges were conducted on samples of clean sand and silty sand. This allowed measurements of linear and nonlinear deformational properties of the same specimen for a wide strain range. The testing program was designed to quantify the effects of void ratio and fines content on the low-strain stiffness of the silty sand as well as on the nonlinear stress-strain relationship and corresponding shear modulus and damping properties as a function of cyclic shear strains. Shear wave velocity, Vs, and maximum shear modulus, Gmax, of silty sand was shown to be significantly smaller than the respective values for clean sands measured at the same void ratio, e, or same relative density, Dr. However, the test results showed that the difference in the level of nonlinearity between clean sand and silty sands was small. For loose samples prepared at an identical relative density, the behaviour of clean sand was slightly less nonlinear as compared to sandy soils with higher fines content. This difference in the nonlinear behaviour of clean sand and sandy soils was negligible for dense soils. Furthermore, no systematic influence of fines content on the material damping curve was observed for sands with fines content FC = 0 to 30%. In order to normalize the effects of fines on moduli of sands, equivalent granular void ratio, e*, was employed. This was done through quantifying the participation of fines content in the force transfer chain of the sand matrix. As such, a unified framework for modelling of the variability of shear wave velocity, Vs, (or shear modulus, Gmax) with void ratio was achieved for clean sands and sands with fines, irrespective of their fines content. Furthermore, modelling of the cyclic stress-strain behaviour based on this experimental program was investigated. The modelling effort focused on developing a simple constitutive model which simultaneously models the soil modulus and damping relationships with shear strains observed in laboratory tests. The backbone curve of the cyclic model was adopted based on a modified version of Kondner and Zelasko (MKZ) hyperbolic function, with a curvature coefficient, a. In order to simulate the hysteretic cycles, the conventional Masing rules (Pyke 1979) were revised. The parameter n, in the Masing’s criteria was assumed to be a function of material damping, h, measured in the laboratory. As such the modulus and damping produced by the numerical model could match the stress-strain behaviour observed in the laboratory over the course of this study. It was shown that the Masing parameter n, is strain-dependent and generally takes values of n ≤ 2. The model was then verified through element test simulations under different cyclic loadings. It was shown that the model could accurately simulate the modulus and the damping simultaneously. The model was then incorporated within the OpenSees computational platform and was used to scrutinize the effects of damping on one-dimensional seismic site response analysis. For this purpose, several strong motion stations which recorded the Canterbury earthquake sequence were selected. The soil profiles were modelled as semi-infinite horizontally layered deposits overlying a uniform half-space subjected to vertically propagating shear waves. The advantages and limitations of the nonlinear model in terms of simulating soil nonlinearity and associated material damping were further scrutinized. It was shown that generally, the conventional Masing criteria unconservatively may underestimate some response parameters such as spectral accelerations. This was shown to be due to larger hysteretic damping modelled by using conventional Masing criteria. In addition, maximum shear strains within the soil profiles were also computed smaller in comparison to the values calculated by the proposed model. Further analyses were performed to study the simulation of backbone curve beyond the strain ranges addressed in the experimental phase of this study. A key issue that was identified was that relying only on the modulus reduction curves to simulate the stress-strain behaviour of soil may not capture the actual soil strength at larger strains. Hence, strength properties of the soil layer should also be incorporated to accurately simulate the backbone curve.
