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.
Impact between structures of bridge sections can play a major, unexpected role in seismic structural damage. Linear and non-linear models are developed to analyze structural impact and response of two single-degree-of-freedom structures, representing adjacent buildings or bridge sections. The analyses presented assess probability of impact, displacement change due to impact, and the probability of increased displacement due to impact. These are assessed over a matrix of structural periods for each degree-of-freedom, different impact coefficients of restitution, and a probabilistically scaled suite of earthquake events. Linear versus non-linear effects are assessed using a Ramberg-Osgood non-linear model for column inelasticity. The normalized distance, or gap-ratio (GR), defined as a percentage of the summed spectral displacements, is used to create probabilistic design requirements. Increasing GR and structural periods that are similar (T2/T1~0.8-1.25) significantly decrease the likelihood of impact, and vice-versa. Including column inelasticity and decreasing coefficient of restitution decrease displacement increases due to impact and thus reduce potential damage. A minimum GR~0.5-0.9 ensures that any displacement increases will be less than 10% for 90% of ground motions over all structural period combinations (0.2-5.0sec). These results enable probabilistic design guidelines to manage undesirable effects of impact– an important factor during the recent Canterbury, New Zealand Earthquakes.
A linear and non-linear model are developed to analyze the structural impact and response of two single degree of freedom structures, representing adjacent buildings or bridge sections. Different impact coefficients of restitution, normalized distances between structures and a range of different structural periods are considered. The probability of impact and the displacement changes that can result from these collisions are computed. The likelihood of an increase in displacement is quantified in a probabilistic sense. A full matrix of response simulations are performed to individually investigate and delineate the effects of inter-structure gap-ratio, period ratios, structural non-linearity and impact elasticity. Column inelasticity is incorporated through the use of a Ramberg-Osgood type hysteresis rule. The minimum normalized distance, or gap-ratio, required between two structures to ensure that the likelihood of increased displacement of more than 10% for either structure for 90% of the given earthquake ground motions is assessed as one of many possible design risk bounds. Increased gap ratio, defined as a percentage of spectral displacement, is shown to reduce the likelihood of impact, as well as close structural periods. Larger differences in the relative periods of the two structures were seen to significantly increase the likelihood of impact. Inclusion of column inelasticity and higher plasticity of impact reduce displacement increases from impact and thus possible further damage to the structures. Such information can be used as a guideline to manage undesirable effects of impact in design - a factor that has been observed to be very important during the recent Canterbury, New Zealand Earthquakes.
Following the 22 February 2011, MW 6.2 earthquake located on a fault beneath the Port Hills of Christchurch, fissuring of up to several hundred metres in length was observed in the loess and loess-colluvium of foot-slope positions in north-facing valleys of the Port Hills. The fissuring was observed in all major valleys, occurred at similar low altitudes, showing a contour-parallel orientation and often accompanied by both lateral compression/extension features and spring formation in the valley floor below. Fissuring locations studied in depth included Bowenvale Valley, Hillsborough Valley, Huntlywood Terrace–Lucas Lane, Bridle Path Road, and Maffeys Road–La Costa Lane. Investigations into loess soil, its properties and mannerisms, as well as international examples of its failure were undertaken, including study of the Loess Plateau of China, the Teton Dam, and palaeo-fissuring on Banks Peninsula. These investigations lead to the conclusion that loess has the propensity to fail, often due to the infiltration of water, the presence of which can lead to its instantaneous disaggregation. Literature study and laboratory analysis of Port Hills loess concluded that is has the ability to be stable in steep, sub-vertical escarpments, and often has a sub-vertically jointed internal structure and has a peak shear strength when dry. Values for cohesion, c (kPa) and the internal friction angle, ϕ (degrees) of Port Hills loess were established. The c values for the 40 Rapaki Road, 3 Glenview Terrace loess samples were 13.4 kPa and 19.7 kPa, respectively. The corresponding ϕ values were thought unusually high, at 42.0° and 43.4°.The analysed loess behaved very plastically, with little or no peak strength visible in the plots as the test went almost directly to residual strength. A geophysics resistivity survey showed an area of low resistivity which likely corresponds to a zone of saturated clayey loess/loess colluvium, indicating a high water table in the area. This is consistent with the appearances of local springs which are located towards the northern end of each distinct section of fissure trace and chemical analysis shows that they are sourced from the Port Hills volcanics. Port Hills fissuring may be sub-divided into three categories, Category A, Category B, and Category C, each characterised by distinctive features of the fissures. Category A includes fissures which display evidence of, spring formation, tunnel-gullying, and lateral spreading-like behaviour or quasi-toppling. These fissures are several metres down-slope of the loess-bedrock interface, and are in valleys containing a loess-colluvium fill. Category B fissures are in wider valleys than those in Category A, and the valleys contain estuarine silty sediments which liquefied during the earthquake. Category C fissures occurred at higher elevations than the fissures in the preceding categories, being almost coincident with bedrock outcropping. It is believed that the mechanism responsible for causing the fissuring is a complex combination of three mechanisms: the trampoline effect, bedrock fracturing, and lateral spreading. These three mechanisms can be applied in varying degrees to each of the fissuring sites in categories A, B, and C, in order to provide explanation for the observations made at each. Toppling failure can describe the soil movement as a consequence of the a three causative mechanisms, and provides insight into the movement of the loess. Intra-loess water coursing and tunnel gullying is thought to have encouraged and exacerbated the fissuring, while not being the driving force per se. Incipient landsliding is considered to be the least likely of the possible fissuring interpretations.
