Low-temperature thermochronology and thermokinematic modeling of deformation, exhumation, and development of topography in the central Southern Alps, New Zealand

Apatite and zircon (U-Th)/He and fission track ages were obtained from ridge transects across the central Southern Alps, New Zealand. Interpretation of local profiles is difficult because relationships between ages and topography or local faults are complex and the data contain large uncertainties, with poor reproducibility between sample duplicates. Data do form regional patterns, however, consistent with theoretical systematics and corroborating previous observations: young Neogene ages occur immediately southeast of the Alpine Fault (the main plate boundary structure on which rocks are exhumed); partially reset ages occur in the central Southern Alps; and older Mesozoic ages occur further toward the southeast. Zircon apparent ages are older than apatite apparent ages for the equivalent method. Three-dimensional thermokinematic modeling of plate convergence incorporates advection of the upper Pacific plate along a low-angle detachment then up an Alpine Fault ramp, adopting a generally accepted tectonic scenario for the Southern Alps. The modeling incorporates heat flow, evolving topography, and the detailed kinetics of different thermochronometric systems and explains both complex local variations and regional patterns. Inclusion of the effects of radiation damage on He diffusion in detrital apatite is shown to have dramatic effects on results. Geometric and velocity parameters are tuned to fit model ages to observed data. Best fit is achieved at 9 mm a−1 plate convergence, with Pacific plate delamination on a gentle 10°SE dipping detachment and more rapid uplift on a 45–60° dipping Alpine Fault ramp from 15 km depth. Thermokinematic modeling suggests dip-slip motion on reverse faults within the Southern Alps should be highest ∼22 km from the Alpine Fault and much lower toward the southeast

Crustal structure and apparent tectonic underplating from receiver function analysis in South Island, New Zealand

We utilize seismic converted phases on more than 700 receiver functions calculated for 42 stations in the South Island, New Zealand, to infer crustal and uppermost mantle structure. We determine the crustal thickness from direct observations of conversion from the Moho interface and infer zone of the maximum thickness being located along the axis of the Southern Alps, just east from the Alpine fault. The crustal root widens from north to south in the direction perpendicular to the Alpine fault and appears to have an asymmetric structure. Stations in the alpine portion of island show evidence for prominent midcrustal conversions. Significant crustal thickening is developed in response to both the convergent component of the motion on the Alpine fault and subduction in the Fiordland region. We propose two models for a strong uppermost mantle conversion that occurs at depths between 33 and 83 km on 16 stations and forms a large continuous feature along the east coast and in the central portions of the South Island. Our preferred model attributes upper mantle conversion to tectonically underplated oceanic crust formed by late Oligocene-Miocene spreading between the Australian and Pacific plates, which was detached from the Australian plate and tectonically underplated under the South Island. An alternative model attributes the upper mantle conversions to long-lived seismic fabric created by subduction of the Gondwanaland margin

Neogene plate tectonic reconstructions and geodynamics of North Island sedimentary basins: Implications for the petroleum systems

Although the modern Australia-Pacific plate boundary through New Zealand is relatively straight, there have been significant changes in its geometry during the Neogene. Within the North Island sector there has been a fundamental transition from an Alpine Fault translation/transpression regime to a Hikurangi margin subduction regime. This transition has been accompanied by the southward encroachment of the edge of the Pacific plate oceanic slab into Australia lithosphere, shortened and thickened along its eastern margin as a consequence of the prior Alpine Fault transpression, the process now operating in South Island. The response of the Australia lithosphere at the surface to the emplacement of the subducted slab at depth, has differed in the East Coast forearc region versus the foreland in western North Island, where the depth to the slab is greater and there has been a characteristic southward migration of depocentres pinned to the leading edge of the slab. The recent publication of new rotation parameters for relative motion of the Australia, Antarctic and Pacific plates, have provided key new data from which to plot the successive emplacement history of the Pacific slab beneath North Island, thus enabling the comparisons to be made with basin stratigraphy and geohistory. These data also constrain the age of subduction initiation at various points along the present trend of the Hikurangi Trough, identifying a younging of subduction initiation to the southwest. An implication of this younging direction is that the modern accretion¬ary prism south of Cape Kidnappers can be no older than late Miocene (c. 11 Ma). The focus of this paper is on new ideas about the tectonic development of North Island and its basins, which have implications for hydrocarbon exploration

Late Holocene landscape change history related to the Alpine Fault determined from drowned forests in Lake Poerua, Westland, New Zealand

Abstract. Lake Poerua is a small, shallow lake that abuts the scarp of the Alpine Fault on the West Coast of New Zealand's South Island. Radiocarbon dates from drowned podocarp trees on the lake floor, a sediment core from a rangefront alluvial fan, and living tree ring ages have been used to deduce the late Holocene history of the lake. Remnant drowned stumps of kahikatea (Dacrycarpus dacrydioides) at 1.7–1.9 m water depth yield a preferred time-of-death age at 1766–1807 AD, while a dryland podocarp and kahikatea stumps at 2.4–2.6 m yield preferred time-of-death ages of ca. 1459–1626 AD. These age ranges are matched to, but offset from, the timings of Alpine Fault rupture events at ca. 1717 AD, and either ca. 1615 or 1430 AD. Alluvial fan detritus dated from a core into the toe of a rangefront alluvial fan, at an equivalent depth to the maximum depth of the modern lake (6.7 m), yields a calibrated age of AD 1223–1413. This age is similar to the timing of an earlier Alpine Fault rupture event at ca. 1230 AD ± 50 yr. Kahikatea trees growing on rangefront fans give ages of up to 270 yr, which is consistent with alluvial fan aggradation following the 1717 AD earthquake. The elevation levels of the lake and fan imply a causal and chronological link between lake-level rise and Alpine Fault rupture. The results of this study suggest that the growth of large, coalescing alluvial fans (Dry and Evans Creek fans) originating from landslides within the rangefront of the Alpine Fault and the rise in the level of Lake Poerua may occur within a decade or so of large Alpine Fault earthquakes that rupture adjacent to this area. These rises have in turn drowned lowland forests that fringed the lake. Radiocarbon chronologies built using OxCal show that a series of massive landscape changes beginning with fault rupture, followed by landsliding, fan sedimentation and lake expansion. However, drowned Kahikatea trees may be poor candidates for intimately dating these events, as they may be able to tolerate water for several decades after metre-scale lake level rises have occurred

Bedrock geology of DFDP-2B, central Alpine Fault, New Zealand

During the second phase of the Alpine Fault, Deep Fault Drilling Project (DFDP) in the Whataroa River, South Westland, New Zealand, bedrock was encountered in the DFDP-2B borehole from 238.5–893.2 m Measured Depth (MD). Continuous sampling and meso- to microscale characterisation of whole rock cuttings established that, in sequence, the borehole sampled amphibolite facies, Torlesse Composite Terrane-derived schists, protomylonites and mylonites, terminating 200–400 m above an Alpine Fault Principal Slip Zone (PSZ) with a maximum dip of 62°. The most diagnostic structural features of increasing PSZ proximity were the occurrence of shear bands and reduction in mean quartz grain sizes. A change in composition to greater mica:quartz + feldspar, most markedly below c. 700 m MD, is inferred to result from either heterogeneous sampling or a change in lithology related to alteration. Major oxide variations suggest the fault-proximal Alpine Fault alteration zone, as previously defined in DFDP-1 core, was not sampled

Earthquake damage in underground roadways

Earthquake damage in underground roadways and mine workings is considered, with particular application to the mines operated by Solid Energy NZ Ltd., on the West Coast of New Zealand’s South Island. The scenario considered is the effect on the mine workings of an earthquake, of moment magnitude eight, being generated by a rupture of the Alpine fault. An empirical relation from the seismology literature is used to relate earthquake magnitude, distance from the epicentre and the peak ground acceleration resulting from the seismic waves. This relation is used to estimate the likely damage at the mine site. Also, the decay scale for Rayleigh (surface) waves is calculated and the implications for the mine workings considered. The two-dimensional scattering of shear (SH) seismic waves from the mine workings is considered. Analytical solutions relevant to various mine tunnel geometries are presented with the stress and displacement amplification, due to scattering from the mine workings, calculated and discussed

Late Holocene landscape change history related to the Alpine Fault determined from drowned forests in Lake Poerua, Westland, New Zealand

Lake Poerua is a small, shallow lake that abuts the scarp of the Alpine Fault on the West Coast of New Zealand’s South Island. Radiocarbon dates from drowned podocarp trees on the lake floor, a sediment core from a rangefront alluvial fan, and living tree ring ages have been used to deduce the late Holocene history of the lake. Remnant drowned stumps of kahikatea (Dacrycarpus dacrydioides) at 1.7–1.9m water depth yield a preferred time-ofdeath age at 1766–1807 AD, while a dryland podocarp and kahikatea stumps at 2.4–2.6m yield preferred time-of-death ages of ca. 1459–1626 AD. These age ranges are matched to, but offset from, the timings of Alpine Fault rupture events at ca. 1717 AD, and either ca. 1615 or 1430 AD. Alluvial fan detritus dated from a core into the toe of a rangefront alluvial fan, at an equivalent depth to the maximum depth of the modern lake (6.7 m), yields a calibrated age of AD 1223–1413. This age is similar to the timing of an earlier Alpine Fault rupture event at ca. 1230AD±50 yr. Kahikatea trees growing on rangefront fans give ages of up to 270 yr, which is consistent with alluvial fan aggradation following the 1717AD earthquake. The elevation levels of the lake and fan imply a causal and chronological link between lake-level rise and Alpine Fault rupture. The results of this study suggest that the growth of large, coalescing alluvial fans (Dry and Evans Creek fans) originating from landslides within the rangefront of the Alpine Fault and the rise in the level of Lake Poerua may occur within a decade or so of large Alpine Fault earthquakes that rupture adjacent to this area. These rises have in turn drowned lowland forests that fringed the lake. Radiocarbon chronologies built using OxCal show that a series of massive landscape changes beginning with fault rupture, followed by landsliding, fan sedimentation and lake expansion. However, drowned Kahikatea trees may be poor candidates for intimately dating these events, as they may be able to tolerate water for several decades after metre-scale lake level rises have occurred

Possible Jurassic age for part of Rakaia Terrane: implications for tectonic development of the Torlesse accretionary prism

Greywacke sandstone and argillite beds comprising Rakaia Terrane (Torlesse Complex) in mid Canterbury, South Island, New Zealand, are widely regarded as Late Triassic (Norian) in age based on the occurrence of Torlessia trace fossils, Monotis, and other taxa. This paleontological age assignment is tested using published 40Ar/39Ar mica and U-Pb zircon ages for these rocks and published and new zircon fission track (FT) ages. The youngest U-Pb zircon ages in the Rakaia Terrane rocks in mid Canterbury are Norian, whereas 10-20% of the 40Ar/39Ar muscovite ages are younger than Norian. Numerical modelling of these mica ages shows that they cannot have originated from partial thermal overprinting in the Torlesse prism if the thermal maximum was short-lived and early in the prism history (210-190 Ma), as commonly inferred for these rocks. The young component of mica ages could, however, be explained by extended residence (200-100 Ma) at 265-290deg.C in the prism. Early Jurassic (c. 189 Ma) zircon FT ages for sandstone beds from Arthur's Pass, the Rakaia valley, and the Hermitage (Mt Cook) are interpreted not to have experienced maximum temperatures above 210deg.C, and therefore cannot have been reduced as a result of partial annealing in the Torlesse prism. This is based on identification of a fossil Cretaceous, zircon FT, partial annealing zone in low-grade schists to the west, and the characteristics of the age data. The Early Jurassic zircon FT ages and the young component of 40Ar/39Ar mica ages are regarded therefore as detrital ages reflecting cooling in the source area, and constrain the maximum depositional age of parts of the Rakaia Terrane in mid Canterbury. The zircon FT data also show the initiation (c. 100 Ma) of marked and widespread Late Cretaceous cooling of Rakaia Terrane throughout Canterbury, which is attributed to uplift and erosion of inboard parts of the Torlesse prism due to continuing subduction accretion at its toe. The critical wedge concept is proposed as a new framework for investigating the development of the Torlesse Complex. The Rakaia Terrane may have formed the core of an accretionary wedge imbricated against the New Zealand margin during the Middle or Late Jurassic. Late Jurassic nonmarine sediments (e.g., Clent Hills Formation) accumulated upon the inner parts of the prism as it enlarged, emerged, and continued to be imbricated. Exhumation of Otago Schist from c. 135 Ma may mark the development of a balance (steady state) between sediments entering the prism at the toe and material exiting at the inboard margin. The enlargement of the area of exhumation to all of Canterbury from c. 100 Ma may reflect a dynamic response to widening of the prism through the accretion of Cretaceous sediments. The model of a dynamic critical wedge may help to explain the various expressions of the Rangitata Orogeny

Metamorphism, argon depletion, heat flow and stress on the Alpine fault

The Alpine fault of New Zealand is a major continental transform fault which was uplifted on its southeast side 4 to 11 km within the last 5 m.y. This uplift has exposed the Haast schists, which were metamorphosed from the adjacent Torlesse graywackes. The Haast schists increase in metamorphic grade from prehnite-pumpellyite facies 9-12 km from the fault through the chlorite and biotite zones of the greenschist facies to the garnet-oligoclase zone amphibolite facies within 4 km of the fault. These metamorphic zone boundaries are subparallel to the fault for 350 km along the strike. The K-Ar and Rb-Sr ages of the schists increase with distance from the fault: from 4 m.y. within 3 km of the fault to approximately 110 m.y. 20 km from the fault. Field relations show that the source of heat that produced the argon depletion aureole was the fault itself

Earthquake source parameters of the 2009 Mw 7.8 Fiordland (New Zealand) earthquake from L-band InSAR observations

The 2009 MW7.8 Fiordland (New Zealand) earthquake is the largest to have occurred in New Zealand since the 1931 Mw 7.8 Hawke’s Bay earthquake, 1 000 km to the northwest. In this paper two tracks of ALOS PALSAR interferograms (one ascending and one descending) are used to determine fault geometry and slip distribution of this large earthquake. Modeling the event as dislocation in an elastic half-space suggests that the earthquake resulted from slip on a SSW-NNE orientated thrust fault that is associated with the subduction between the Pacific and Australian Plates, with oblique displacement of up to 6.3 m. This finding is consistent with the preliminary studies undertaken by the USGS using seismic data