We examined the stratigraphy of alluvial fans formed at the steep range front of the Southern Alps at Te Taho, on the north bank of the Whataroa River in central West Coast, South Island, New Zealand. The range front coincides with the Alpine Fault, an Australian-Pacific plate boundary fault, which produces regular earthquakes. Our study of range front fans revealed aggradation at 100- to 300-year intervals. Radiocarbon ages and soil residence times (SRTs) estimated by a quantitative profile development index allowed us to elucidate the characteristics of four episodes of aggradation since 1000 CE. We postulate a repeating mode of fan behaviour (fan response cycle [FRC]) linked to earthquake cycles via earthquake-triggered landslides. FRCs are characterised by short response time (aggradation followed by incision) and a long phase when channels are entrenched and fan surfaces are stable (persistence time). Currently, the Te Taho and Whataroa River fans are in the latter phase. The four episodes of fan building we determined from an OxCal sequence model correlate to Alpine Fault earthquakes (or other subsidiary events) and support prior landscape evolution studies indicating ≥M7.5 earthquakes as the main driver of episodic sedimentation. Our findings are consistent with other historic non-earthquake events on the West Coast but indicate faster responses than other earthquake sites in New Zealand and elsewhere where rainfall and stream gradients (the basis for stream power) are lower. Judging from the thickness of fan deposits and the short response times, we conclude that pastoral farming (current land-use) on the fans and probably across much of the Whataroa River fan would be impossible for several decades after a major earthquake. The sustainability of regional tourism and agriculture is at risk, more so because of the vulnerability of the single through road in the region (State Highway 6).
This study investigates evidence for linkages and fault interactions centred on the Cust Anticline in Northwest Canterbury between Starvation Hill to the southwest and the Ashley and Loburn faults to the northeast. An integrated programme of geologic, geomorphic, paleo-seismic and geophysical analyses was undertaken owing to a lack of surface exposures and difficulty in distinguishing active tectonic features from fluvial and/or aeolian features across the low-relief Canterbury Plains. LiDAR analysis identified surface expression of several previously unrecognised active fault traces across the low-relief aggradation surfaces of the Canterbury Plains. Their presence is consistent with predictions of a fault relay exploiting the structural mesh across the region. This is characterised by interactions of northeast-striking contractional faults and a series of re-activating inherited Late Cretaceous normal faults, the latter now functioning as E–W-striking dextral transpressive faults. LiDAR also allowed for detailed analysis of the surface expression of individual faults and folds across the Cust Anticline contractional restraining bend, which is evolving as a pop-up structure within the newly established dextral shear system that is exploiting the inherited, now re-activated, basement fault zone. Paleo-seismic trenches were located on the crest of the western arm of the Cust Anticline and across a previously unrecognised E–W-striking fault trace, immediately southwest of the steeply plunging Cust Anticline termination. These studies confirmed the location and structural style of north-northeast-striking faults and an E–W-striking fault associated with the development of this structural culmination. A review of available industry seismic reflection lines emphasised the presence of a series of common structural styles having the same underlying structural drivers but with varying degrees of development and expression, both in the seismic profiles and in surface elevations across the study area. Based on LiDAR surface mapping and preliminary re-analysis of industry seismic reflection data, four fault zones are identified across the restraining bend structural culminations, which together form the proposed Oxford–Cust–Ashley Fault System. The 2010–2012 Canterbury Earthquake Sequence showed many similarities to the structural pattern established across the Oxford–Cust–Ashley Fault System, emphasising the importance of identification and characterization of presently hidden fault sources, and the understanding of fault network linkages, in order to improve constraints on earthquake source potential. Improved understanding of potentially-interactive fault sources in Northwest Canterbury, with the potential for combined initial fault rupture and spatial and temporal rupture propagation across this fault system, can be used in probabilistic seismic hazard analysis for the region, which is essential for the suitability and sustainability of future social and economic development.
Saltwater Forest is a Dacrydium cupressinum-dominated lowland forest covering 9000 ha in south Westland, South Island, New Zealand. Four thousand hectares is managed for sustainable production of indigenous timber. The aim of this study was to provide an integrated analysis of soils, soil-landform relationships, and soil-vegetation relationships at broad and detailed scales. The broad scale understandings provide a framework in which existing or future studies can be placed and the detailed studies elucidate sources of soil and forest variability.
Glacial landforms dominate. They include late Pleistocene lateral, terminal and ablation moraines, and outwash aggradation and degradation terraces. Deposits and landforms from six glacial advances have been recognised ranging from latest Last (Otira) Glaciation to Penultimate (Waimea) Glaciation. The absolute ages of landforms were established by analysis of the thickness and soil stratigraphy of loess coverbeds, augmented with radiocarbon dating and phytolith and pollen analysis.
In the prevailing high rainfall of Westland soil formation is rapid. The rate of loess accretion in Saltwater
Forest (ca. 30 mm ka⁻¹) has been low enough that soil formation and loess accretion took place contemporaneously. Soils formed in this manner are known as upbuilding soils. The significant difference between upbuilding pedogenesis and pedogenesis in a topdown sense into an existing sediment body is that each subsoil increment of an upbuilding soil has experienced processes of all horizons above. In Saltwater Forest subsoils of upbuilding soils are strongly altered because they have experienced the extremely acid environment of the soil surface at some earlier time. Some soil chronosequence studies in Westland have included upbuilding soils formed in loess as the older members of the sequence. Rates and types of processes inferred from these soils should be reviewed because upbuilding is a different pedogenic pathway to topdown pedogenesis.
Landform age and morphology were used as a primary stratification for a study of the soil pattern and nature of soil variability in the 4000 ha production area of Saltwater Forest. The age of landforms (> 14 ka) and rapid soil formation mean that soils are uniformly strongly weathered and leached. Soils include Humic Organic Soils, Perch-gley Podzols, Acid Gley Soils, Allophanic Brown Soils, and Orthic or Pan Podzols. The major influence on the nature of soils is site hydrology which is determined by macroscale features of landforms (slope, relief, drainage density), mesoscale effects related to position on landforms, and microscale influences determined by microtopography and individual tree effects. Much of the soil variability arises at microscales so that it is not possible to map areas of uniform soils at practical map scales. The distribution of soil variability across spatial scales, in relation to the intensity of forest management, dictates that it is most appropriate to map soil complexes with boundaries coinciding with landforms.
Disturbance of canopy trees is an important agent in forest dynamics. The frequency of forest disturbance in the production area of Saltwater Forest varies in a systematic way among landforms in accord with changes in abundance of different soils. The frequency of forest turnover is highest on landforms with the greatest abundance of extremely poorly-drained Organic Soils. As the abundance of better-drained soils increases the frequency of forest turnover declines. Changes in turnover frequency are reflected in the mean size and density of canopy trees (Dacrydium cupressinum) among landforms. Terrace and ablation moraine landforms with the greatest abundance of extremely poorly-drained soils have on average the smallest trees growing most densely. The steep lateral moraines, characterised by well drained soils, have fewer, larger trees. The changes manifested at the landform scale are an integration of processes operating over much shorter range as a result of short-range soil variability. The systematic changes in forest structure and turnover frequency among landforms and soils have important implications for sustainable forest management.
The previously unknown Greendale Fault was buried beneath the Canterbury Plains and ruptured in the September 4th 2010 moment magnitude (Mw) 7.1 Darfield Earthquake. The Darfield Earthquake and subsequent Mw 6 or greater events that caused damage to Christchurch highlight the importance of unmapped faults near urban areas. This thesis examines the morphology, age and origin of the Canterbury Plains together with the paleoseismology and surface-rupture displacement distributions of the Greendale Fault. It offers new insights into the surface-rupture characteristics, paleoseismology and recurrence interval of the Greendale Fault and related structures involved in the 2010 Darfield Earthquake. To help constrain the timing of the penultimate event on the Greendale Fault the origin and age of the faulted glacial outwash deposits have been examined using sedimentological analysis of gravels and optically stimulated luminescence (OSL) dating combined with analysis of GPS and LiDAR survey data. OSL ages from this and other studies, and the analysis of surface paleochannel morphology and subsurface gravel deposits indicate distinct episodes of glacial outwash activity across the Canterbury Plains, at ~20 to 24 and ~28 to 33 kyr separated by a hiatus in sedimentation possibly indicating an interstadial period. These data suggest multiple glacial periods between ~18 and 35 kyr which may have occurred throughout the Canterbury region and wider New Zealand. A new model for the Waimakariri Fan is proposed where aggradation is mainly achieved during episodic sheet flooding with the primary river channel location remaining approximately fixed. The timing, recurrence interval and displacements of the penultimate surface-rupturing earthquake on the Greendale Fault have been constrained by trenching the scarp produced in 2010 at two locations. These excavations reveal a doubling of the magnitude of surface displacement at depths of 2-4 m. Aided by OSL ages of sand lenses in the gravel deposits, this factor-of-two increase is interpreted to indicate that in the central section of the Greendale Fault the penultimate surface-rupturing event occurred between ca. 20 and 30 kyr ago. The Greendale Fault remained undetected prior to the Darfield earthquake because the penultimate fault scarp was eroded and buried during Late Pleistocene alluvial activity. The Darfield earthquake rupture terminated against the Hororata Anticline Fault (HAF) in the west and resulted in up to 400 mm of uplift on the Hororata Anticline immediately above the HAF. Folding in 2010 is compared to Quaternary and younger deformation across the anticline recorded by a seismic reflection line, GPS-measured topographic profiles along fluvial surfaces, and river channel sinuosity and morphology. It is concluded that the HAF can rupture during earthquakes dissimilar to the 2010 event that may not be triggered by slip on the Greendale Fault. Like the Greendale Fault geomorphic analyses provide no evidence for rupture of the HAF in the last 18 kyr, with the average recurrence interval for the late Quaternary inferred to be at least ~10 kyr. Surface rupture of the Greendale Fault during the Darfield Earthquake produced one of the most accessible and best documented active fault displacement and geometry datasets in the world. Surface rupture fracture patterns and displacements along the fault were measured with high precision using real time kinematic (RTK) GPS, tape and compass, airborne light detection and ranging (LiDAR), and aerial photos. This allowed for detailed analysis of the cumulative strike-slip displacement across the fault zone, displacement gradient (ground shear strain) and the type of displacement (i.e. faulting or folding). These strain profiles confirm that the rupture zone is generally wide (~30 to ~300 metres) with >50% of displacement (often 70-80%) accommodated by ground flexure rather than discrete fault slip and ground cracking. The greatest fault-zone widths and highest proportions of folding are observed at fault stepovers.
Geologic phenomena produced by earthquake shaking, including rockfalls and liquefaction features, provide important information on the intensity and spatiotemporal distribution of earthquake ground motions. The study of rockfall and liquefaction features produced in contemporary well- instrumented earthquakes increases our knowledge of how natural and anthropogenic environments respond to earthquakes and improves our ability to deduce seismologic information from analogous pre-contemporary (paleo-) geologic features. The study of contemporary and paleo- rockfall and liquefaction features enables improved forecasting of environmental responses to future earthquakes. In this thesis I utilize a combination of field and imagery-based mapping, trenching, stratigraphy, and numerical dating techniques to understand the nature and timing of rockfalls (and hillslope sedimentation) and liquefaction in the eastern South Island of New Zealand, and to examine the influence that anthropogenic activity has had on the geologic expressions of earthquake phenomena. At Rapaki (Banks Peninsula, NZ), field and imagery-based mapping, statistical analysis and numerical modeling was conducted on rockfall boulders triggered by the fatal 2011 Christchurch earthquakes (n=285) and compared with newly identified prehistoric (Holocene and Pleistocene) boulders (n=1049) deposited on the same hillslope. A significant population of modern boulders (n=26) travelled farther downslope (>150 m) than their most-travelled prehistoric counterparts, causing extensive damage to residential dwellings at the foot of the hillslope. Replication of prehistoric boulder distributions using 3-dimensional rigid body numerical models requires the application of a drag-coefficient, attributed to moderate to dense slope vegetation, to account for their spatial distribution. Radiocarbon dating provides evidence for 17th to early 20th century deforestation at the study site during Polynesian and European colonization and after emplacement of prehistoric rockfalls. Anthropocene deforestation enabled modern rockfalls to exceed the limits of their prehistoric predecessors, highlighting a shift in the geologic expression of rockfalls due to anthropogenic activity. Optical and radiocarbon dating of loessic hillslope sediments in New Zealand’s South Island is used to constrain the timing of prehistoric rockfalls and associated seismic events, and quantify spatial and temporal patterns of hillslope sedimentation including responses to seismic and anthropogenic forcing. Luminescence ages from loessic sediments constrain timing of boulder emplacement to between ~3.0 and ~12.5 ka, well before the arrival of Polynesians (ca AD 1280) and Europeans (ca AD 1800) in New Zealand, and suggest loess accumulation was continuing at the study site until 12-13 ka. Large (>5 m3) prehistoric rockfall boulders preserve an important record of Holocene hillslope sedimentation by creating local traps for sediment aggradation and upbuilding soil formation. Sediment accumulation rates increased considerably (>~10 factor increase) following human arrival and associated anthropogenic burning of hillslope vegetation. New numerical ages are presented to place the evolution of loess-mantled hillslopes in New Zealand’s South Island into a longer temporal framework and highlight the roles of earthquakes and humans on hillslope surface process. Extensive field mapping and characterization for 1733 individual prehistoric rockfall boulders was conducted at Rapaki and another Banks Peninsula site, Purau, to understand their origin, frequency, and spatial and volumetric distributions. Boulder characteristics and distributions were compared to 421 boulders deposited at the same sites during the 2010-2011 Canterbury earthquake sequence. Prehistoric boulders at Rapaki and Purau are comprised of two dominant lithofacies types: volcanic breccia and massive (coherent) lava basalt. Volcanic breccia boulders are found in greatest abundance (64-73% of total mapped rockfall) and volume (~90-96% of total rockfall) at both locations and exclusively comprise the largest boulders with the longest runout distances that pose the greatest hazard to life and property. This study highlights the primary influence that volcanic lithofacies architecture has on rockfall hazard. The influence of anthropogenic modifications on the surface and subsurface geologic expression of contemporary liquefaction created during the 2010-2011 Canterbury earthquake sequence (CES) in eastern Christchurch is examined. Trench observations indicate that anthropogenic fill layer boundaries and the composition/texture of discretely placed fill layers play an important role in absorbing fluidized sand/silt and controlling the subsurface architecture of preserved liquefaction features. Surface liquefaction morphologies (i.e. sand blows and linear sand blow arrays) display alignment with existing utility lines and utility excavations (and perforated pipes) provided conduits for liquefaction ejecta during the CES. No evidence of pre-CES liquefaction was identified within the anthropogenic fill layers or underlying native sediment. Radiocarbon dating of charcoal within the youngest native sediment suggests liquefaction has not occurred at the study site for at least the past 750-800 years. The importance of systematically examining the impact of buried infrastructure on channelizing and influencing surface and subsurface liquefaction morphologies is demonstrated. This thesis highlights the importance of using a multi-technique approach for understanding prehistoric and contemporary earthquake phenomena and emphasizes the critical role that humans play in shaping the geologic record and Earth’s surface processes.
