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.
Following the 22nd February 2011, Mw 6.2 earthquake located along a previously unknown fault beneath the Port Hills of Christchurch, surface cracking was identified in contour parallel locations within fill material at Quarry Road on the lower slopes of Mount Pleasant. GNS Science, in the role of advisor to the Christchurch City Council, concluded that these cracks were a part of a potential rotational mass movement (named zone 11A) within the fill and airfall loess material present. However, a lack of field evidence for slope instability and an absence of laboratory geotechnical data on which slope stability analysis was based, suggested this conclusion is potentially incorrect. It was hypothesised that ground cracking was in fact due to earthquake shaking, and not mass movement within the slope, thus forming the basis of this study. Three soil units were identified during surface and subsurface investigations at Quarry Road: fill derived from quarry operations in the adjacent St. Andrews Quarry (between 1893 and 1913), a buried topsoil, and underlying in-situ airfall loess. The fill material was identified by the presence of organic-rich topsoil “clods” that were irregular in both size (∼10 – 200 mm) and shape, with variable thicknesses of 1 – 10 m. Maximum thickness, as indicated by drill holes and geophysical survey lines, was identified below 6 Quarry Road and 7 The Brae where it is thought to infill a pre-existing gully formed in the underlying airfall loess. Bearing strength of the fill consistently exceeded 300 kPa ultimate below ∼500 mm depth. The buried topsoil was 200 – 300 mm thick, and normally displayed a lower bearing strength when encountered, but not below 300 kPa ultimate (3 – 11 blows per 100mm or ≥100 kPa allowable). In-situ airfall loess stood vertically in outcrop due to its characteristic high dry strength and also showed Scala penetrometer values of 6 – 20+ blows per 100 mm (450 – ≥1000 kPa ultimate). All soils were described as being moist to dry during subsurface investigations, with no groundwater table identified during any investigation into volcanic bedrock. In-situ moisture contents were established using bulk disturbed samples from hand augers and test pitting. Average moisture contents were low at 9% within the fill, 11 % within the buried topsoil, and 8% within the airfall loess: all were below the associated average plastic limit of 17, 15, and 16, respectively, determined during Atterberg limit analysis. Particle size distributions, identified using the sieve and pipette method, were similar between the three soil units with 11 – 20 % clay, 62 – 78 % silt, and 11 – 20 % fine sand. Using these results and the NZGS soil classification, the loess derived fill and in-situ airfall loess are termed SILT with some clay and sand, and the buried topsoil is SILT with minor clay and sand. Dispersivity of the units was found using the Emerson crumb test, which established that the fill can be non- to completely dispersive (score 0 – 4). The buried topsoil was always non-dispersive (score 0), and airfall loess completely dispersive (score 4). Values for cohesion (c) and internal friction angle (φ) of the three soil units were established using the direct shear box at field moisture contents. Results showed all soil units had high shear strengths at the moisture contents tested (c = 18 – 24 kPa and φ = 42 – 50°), with samples behaving in a brittle fashion. Moisture content was artificially increased to 16% within the buried topsoil, which reduced the shear strength (c = 10 kPa, φ = 18°) and allowed it to behave plastically. Observational information indicating stability at Quarry Road included: shallow, discontinuous, cracks that do not display vertical offset; no scarp features or compressional zones typical of landsliding; no tilted or deformed structures; no movement in inclinometers; no basal shear zone identified in logged core to 20 m depth; low field moisture contents; no groundwater table; and high soil strength using Scala penetrometers. Limit equilibrium analysis of the slope was conducted using Rocscience software Slide 5.0 to verify the slope stability identified by observational methods. Friction, cohesion, and density values determined during laboratory were input into the two slope models investigated. Results gave minimum static factor of safety values for translational (along buried topsoil) and rotational (in the fill) slides of 2.4 – 4.2. Sensitivity of the slope to reduced shear strength parameters was analysed using c = 10 kPa and φ = 18° for the translational buried topsoil plane, and a cohesion of 0 kPa within the fill for the rotational plane. The only situation that gave a factor of safety <1.0 was in nonengineered fill at 0.5 m depth. Pseudostatic analysis based on previous peak ground acceleration (PGA) values for the Canterbury Earthquake Sequence, and predicted PGAs for future Alpine Fault and Hope Fault earthquakes established minimum factor of safety values between 1.2 and 3.3. Yield acceleration PGAs were computed to be between 0.8g and 1.6g. Based on all information gathered, the cracking at Quarry Road is considered to be shallow deformation in response to earthquake shaking, and not due to deep-seated landsliding. It is recommended that the currently bare site be managed by smoothing the land, installing contour drainage, and bioremediation of the surface soils to reduce surface water infiltration and runoff. Extensive earthworks, including removal of the fill, are considered unnecessary. Any future replacement of housing would be subject to site-specific investigations, and careful foundation design based on those results.
Surface rupture and slip from the Mw 7.8 2016 Kaikōura Earthquake have been mapped in the region between the Leader and Charwell rivers using field mapping and LiDAR data. The eastern Humps, north Leader and Conway-Charwell faults ruptured the ground surface in the study area. The E-NE striking ‘The Humps’ Fault runs along the base of the Mt Stewart range front, appears to dip steeply NW and intersects the NNW-NNE Leader Fault which itself terminates northwards at the NE striking Conway-Charwell Fault. The eastern Humps Fault is up to the NW and accommodates oblique slip with reverse and right lateral displacement. Net slip on ‘The Humps’ Fault is ≤4 m and produced ≤4 m uplift of the Mt Stewart range during the earthquake. The Leader Fault strikes NNW-NNE with dips ranging from ~10° west to 80° east and accommodated ≤4 m net slip comprising left-lateral and up-to-the-west vertical displacement. Like the Humps west of the study area, surface-rupture of the Leader Fault occurred on multiple strands. The complexity of rupture on the Leader Fault is in part due to the occurrence of bedding-parallel slip within the Cretaceous-Cenozoic sequence. Although the Mt Stewart range front is bounded by ‘The Humps’ Fault, in the study area neither this fault nor the Leader Fault were known to have been active before the earthquake. Fieldwork and trenching investigations are ongoing to characterise the geometry, kinematics and paleoseismic history of the mapped active faults.
In 2010 and 2011 Christchurch, New Zealand experienced a series of earthquakes that caused extensive damage across the city, but primarily to the Central Business District (CBD) and eastern suburbs. A major feature of the observed damage was extensive and severe soil liquefaction and associated ground damage, affecting buildings and infrastructure. The behaviour of soil during earthquake loading is a complex phenomena that can be most comprehensively analysed through advanced numerical simulations to aid engineers in the design of important buildings and critical facilities. These numerical simulations are highly dependent on the capabilities of the constitutive soil model to replicate the salient features of sand behaviour during cyclic loading, including liquefaction and cyclic mobility, such as the Stress-Density model. For robust analyses advanced soil models require extensive testing to derive engineering parameters under varying loading conditions for calibration. Prior to this research project little testing on Christchurch sands had been completed, and none from natural samples containing important features such as fabric and structure of the sand that may be influenced by the unique stress-history of the deposit. This research programme is focussed on the characterisation of Christchurch sands, as typically found in the CBD, to facilitate advanced soil modelling in both res earch and engineering practice - to simulate earthquake loading on proposed foundation design solutions including expensive ground improvement treatments. This has involved the use of a new Gel Push (GP) sampler to obtain undisturbed samples from below the ground-water table. Due to the variable nature of fluvial deposition, samples with a wide range of soil gradations, and accordingly soil index properties, were obtained from the sampling sites. The quality of the samples is comprehensively examined using available data from the ground investigation and laboratory testing. A meta-quality assessment was considered whereby a each method of evaluation contributed to the final quality index assigned to the specimen. The sampling sites were characterised with available geotechnical field-based test data, primarily the Cone Penetrometer Test (CPT), supported by borehole sampling and shear-wave velocity testing. This characterisation provides a geo- logical context to the sampling sites and samples obtained for element testing. It also facilitated the evaluation of sample quality. The sampling sites were evaluated for liquefaction hazard using the industry standard empirical procedures, and showed good correlation to observations made following the 22 February 2011 earthquake. However, the empirical method over-predicted liquefaction occurrence during the preceding 4 September 2010 event, and under-predicted for the subsequent 13 June 2011 event. The reasons for these discrepancies are discussed. The response of the GP samples to monotonic and cyclic loading was measured in the laboratory through triaxial testing at the University of Canterbury geomechanics laboratory. The undisturbed samples were compared to reconstituted specimens formed in the lab in an attempt to quantify the effect of fabric and structure in the Christchurch sands. Further testing of moist tamped re- constituted specimens (MT) was conducted to define important state parameters and state-dependent properties including the Critical State Line (CSL), and the stress-strain curve for varying state index. To account for the wide-ranging soil gradations, selected representative specimens were used to define four distinct CSL. The input parameters for the Stress-Density Model (S-D) were derived from a suite of tests performed on each representative soil, and with reference to available GP sample data. The results of testing were scrutinised by comparing the data against expected trends. The influence of fabric and structure of the GP samples was observed to result in similar cyclic strength curves at 5 % Double Amplitude (DA) strain criteria, however on close inspection of the test data, clear differences emerged. The natural samples exhibited higher compressibility during initial loading cycles, but thereafter typically exhibited steady growth of plastic strain and excess pore water pressure towards and beyond the strain criteria and initial liquefaction, and no flow was observed. By contrast the reconstituted specimens exhibited a stiffer response during initial loading cycles, but exponential growth in strains and associated excess pore water pressure beyond phase-transformation, and particularly after initial liquefaction where large strains were mobilised in subsequent cycles. These behavioural differences were not well characterised by the cyclic strength curve at 5 % DA strain level, which showed a similar strength for both GP samples and MT specimens. A preliminary calibration of the S-D model for a range of soil gradations is derived from the suite of laboratory test data. Issues encountered include the influence of natural structure on the peak-strength–state index relationship, resulting in much higher peak strengths than typically observed for sands in the literature. For the S-D model this resulted in excessive stiffness to be modelled during cyclic mobility, when the state index becomes large momentarily, causing strain development to halt. This behaviour prevented modelling the observed re- sponse of silty sands to large strains, synonymous with “liquefaction”. Efforts to reduce this effect within the current formulation are proposed as well as future research to address this issue.
This thesis is concerned with modelling rockfall parameters associated with cliff collapse debris and the resultant “ramp” that formed following the high peak ground acceleration (PGA) events of 22 February 2011 and 13 June 2011. The Christchurch suburb of Redcliffs, located at the base of the Port Hills on the northern side of Banks Peninsula, New Zealand, is comprised of Miocene-age volcanics with valley-floor infilling marine sediments. The area is dominated by basaltic lava flows of the Mt Pleasant Formation, which is a suite of rocks forming part of the Lyttelton Volcanic Group that were erupted 11.0-10.0Ma. Fresh exposure enabled the identification of a basaltic ignimbrite unit at the study site overlying an orange tuff unit that forms a marker horizon spanning the length of the field area. Prior to this thesis, basaltic ignimbrite on Banks Peninsula has not been recorded, so descriptions and interpretations of this unit are the first presented. Mapping of the cliff face by remote observation, and analysis of hand samples collected from the base of the debris slopes, has identified a very strong (>200MPa), columnar-jointed, welded unit, and a very weak (<5MPa), massive, so-called brecciated unit that together represent the end-member components of the basaltic ignimbrite. Geochemical analysis shows the welded unit is picrite basalt, and the brecciated unit is hawaiite, making both clearly distinguishable from the underlying trachyandesite tuff. RocFall™ 4.0 was used to model future rockfalls at Redcliffs. RocFall™ is a two-dimensional (2D), hybrid, probabilistic modelling programme for which topographical profile data is used to generate slope profiles. GNS Science collected the data used for slope profile input in March 2011. An initial sensitivity analysis proved the Terrestrial Laser Scan (TLS)-derived slope to be too detailed to show any results when the slope roughness parameter was tested. A simplified slope profile enabled slope roughness to be varied, however the resulting model did not correlate with field observations as well. By using slope profile data from March 2011, modelled rockfall behaviour has been calibrated with observed rockfall runout at Redcliffs in the 13 June 2011 event to create a more accurate rockfall model. The rockfall model was developed on a single slope profile (Section E), with the chosen model then applied to four other section lines (A-D) to test the accuracy of the model, and to assess future rockfall runout across a wider area. Results from Section Lines A, B, and E correlate very well with field observations, with <=5% runout exceeding the modelled slope, and maximum bounce height at the toe of the slope <=1m. This is considered to lie within observed limits given the expectation that talus slopes will act as a ramp on which modelled rocks travel further downslope. Section Lines C and D produced higher runout percentage values than the other three section lines (23% and 85% exceeding the base of the slope, respectively). Section D also has a much higher maximum bounce height at the toe of the slope (~8.0m above the slope compared to <=1.0m for the other four sections). Results from modelling of all sections shows the significance of the ratio between total cliff height (H) and horizontal slope distance (x), and of maximum drop height to the top of the talus (H*) and horizontal slope distance (x). H/x can be applied to the horizontal to vertical ratio (H:V) as used commonly to identify potential slope instability. Using the maximum value from modelling at Redcliffs, the future runout limit can be identified by applying a 1.4H:1V ratio to the remainder of the cliff face. Additionally, the H*/x parameter shows that when H*/x >=0.6, the percentage of rock runout passing the toe of the slope will exceed 5%. When H*/x >=0.75, the maximum bounce height at the toe of the slope can be far greater than when H*/x is below this threshold. Both of these parameters can be easily obtained, and can contribute valuable guideline data to inform future land-use planning decisions. This thesis project has demonstrated the applicability of a 2D probabilistic-based model (RocFall™ 4.0) to evaluate rockfall runout on the talus slope (or ramp) at the base of ~35-70m high cliff with a basaltic ignimbrite source. Limitations of the modelling programme have been identified, in particular difficulties with adjusting modelled roughness of the slope profile and the inability to consider fragmentation. The runout profile using RocFall™ has been successfully calibrated against actual profiles and some anomalous results have been identified.
