High‐ to ultrahigh‐temperature metamorphism in the lower crust: An example resulting from Hikurangi Plateau collision and slab rollback in New Zealand

Lower crustal xenoliths erupted from an intraplate diatreme reveal that a portion of the New Zealand Gondwana margin experienced high‐temperature (HT) to ultrahigh‐temperature (UHT) granulite facies metamorphism just after flat slab subduction ceased at c. 110–105 Ma. P–T calculations for garnet–orthopyroxene‐bearing felsic granulite xenoliths indicate equilibration at ~815 to 910°C and 0.7 to 0.8 GPa, with garnet‐bearing mafic granulite xenoliths yielding at least 900°C. Supporting evidence for the attainment of HT and UHT conditions in felsic granulite comes from re‐integration of exsolution in feldspar (~900–950°C at 0.8 GPa), Ti‐in‐zircon thermometry on Y‐depleted overgrowths on detrital zircon grains (932°C ± 24°C at aTiO₂ = 0.8 ± 0.2), and correlation of observed assemblages and mineral compositions with thermodynamic modelling results (≥850°C at 0.7 to 0.8 GPa). The thin zircon overgrowths, which were mainly targeted by drilling through the cores of grains, yield a U–Pb pooled age of 91.7 ± 2.0 Ma. The cause of Late Cretaceous HT‐UHT metamorphism on the Zealandia Gondwana margin is attributed to collision and partial subduction of the buoyant oceanic Hikurangi Plateau in the Early Cretaceous. The halt of subduction caused the fore‐running shallowly dipping slab to rollback towards the trench position and permitted the upper mantle to rapidly increase the geothermal gradient through the base of the extending (former) accretionary prism. This sequence of events provides a mechanism for achieving regional HT–UHT conditions in the lower crust with little or no sign of this event at the surface.     

Subduction and back-arc activity at the Hikurangi convergent Margin,New Zealand

The Hikurangi Margin is a region of oblique subduction with northwest-dipping intermediate depth seismicity extending southwest from the Kermadec system to about 42°S. The current episode of subduction is at least 16–20 Ma old. The plate convergence rate varies along the margin from about 60 mm/a at the south end of the Kermadec Trench to about 45 mm/a at 42°S. The age of the Pacific lithosphere adjacent to the Hikurangi Trench is not known.The margin divides at about latitude 39°S into two quite dissimilar parts. The northern part has experienced andesitic volcanism for about 18 Ma, and back-arc extension in the last 4 Ma that has produced a back-arc basin onshore with high heaflow, thin crust and low upper-mantle seismic velocities. The extension appears to have arisen from a seawards migration of the Hikurangi Trench north of 39°S. Here the plate interface is thought to be currently uncoupled, as geodetic data indicate extension of the fore-arc basin, and historic earthquakes have not exceededM _s=7.South of 39°S there is no volcanism and a back-arc basin has been produced by downward flexure of the lithosphere due to strong coupling with the subducting plate. Heatflow in the basin is normal. Evidence for strong coupling comes from historic earthquakes of up to aboutM _s=8 and high rates of uplift on the southeast coast of the North Island.The reason for this division of the margin is not known but may be related to an inferred increase, from northeast to southwest, in the buoyancy of the Pacific lithosphere.     

Distribution of contourite drifts on convergent margins: Examples from the Hikurangi subduction margin of New Zealand

Contourite drift systems form a significant component of the marine clastic sedimentary record. Although contourites form in all tectonic settings, few studies have described their development along convergent margins; such characterization is needed to underpin oceanographic and palaeoenvironmental studies in active settings. This study is the first to document contourite drift development along the Hikurangi subduction margin of New Zealand. Integration of bathymetric, seismic and well data enables five classes of drift to be recognized around the subduction wedge, occurring in three principal associations: (i) an upper slope drift association of giant elongate mounded (ca 150 km long, 50 km wide and up to 1100 m thick) and plastered drifts (ca 300 km long, 8 km wide and <600 m thick), which occurs upon and inboard of a major intrabasinal thrust‐cored high, whose long axis parallels the coast; shallow bottom currents disperse sub‐parallel to this axis; (ii) a spatiotemporally discontinuous association of confined and mounded hybrid drifts (ca 500 m long, <2 km wide and up to 500 m thick) that occurs along the mid‐to‐outer slope domain of the wedge, recording the interaction of along‐slope and downslope currents within trench‐slope basins; and (iii) a trench fill assemblage that implies the passage of abyssal bottom currents across a 40 km reach of the trench‐axial Hikurangi Channel‐levée, with associated modification of the channel form and of overbank sediment waves. The fundamental presence of contourites along this margin appears to depend on the orientation and strength of oceanographic bottom currents. However, drift type and evolution vary depending on the slope gradient and the presence of irregular seafloor topography created by tectonic structures. The documented drifts are generally smaller, less continuous, and develop more intermittently than similar styles of drifts documented on passive margins; this mode of occurrence may be characteristic of contourite development on convergent margins.     

Gas hydrate formation in compressional,extensional and un-faulted structural settings – Examples from New Zealand's Hikurangi margin

We investigate gas hydrate formation processes in compressional, extensional and un-faulted settings on New Zealand's Hikurangi margin using seismic reflection data. The compressional setting is characterized by a prominent subduction wedge thrust fault that terminates beneath the base of gas hydrate stability, as determined from a bottom-simulating reflection (BSR). The thrust is surrounded by steeply dipping strata that cross the BSR at a high angle. Above the BSR, these strata are associated with a high velocity anomaly that is likely indicative of relatively concentrated, and broadly distributed, gas hydrates. The un-faulted setting—sedimentary infill of a slope basin on the landward side of a prominent thrust ridge—is characterized by a strong BSR, a thick underlying free gas zone, and short positive polarity reflection segments that extend upward from the BSR. We interpret the short reflection segments as the manifestation of gas hydrates within relatively coarse-grained sediments. The extensional setting is a localized, shallow response to flexural bending of strata within an anticline. Gas has accumulated beneath the BSR in the apex of folding. A high-velocity zone directly above the BSR is probably mostly lithologically-derived, and only partly related to gas hydrates. Although each setting shows evidence for focused gas migration into the gas hydrate stability zone, we interpret that the compressional tectonic setting is most likely to contain concentrated gas hydrates over a broad region. Indeed, it is the only setting associated with a deep-reaching fault, meaning it is the most likely of the three settings to have thermogenic gas contributing to hydrate formation. Our results highlight the importance of anisotropic permeability in layered sediments and the role this plays in directing sub-surface fluid flow, and ultimately in the distribution of gas hydrate. Each of the three settings we describe would warrant further investigation in any future consideration of gas hydrates as an energy resource on the Hikurangi margin.     

The lithospheric geodynamics of plate boundary transpression in New Zealand: Initiating and emplacing subduction along the Hikurangi margin, and the tectonic evolution of the Alpine Fault system

In contrast to the normal ‘Wilson cycle’ sequence of subduction leading to continental collision and associated mountain building, the evolution of the New Zealand plate boundary in the Neogene reflects the converse—initially a period of continental convergence that is followed by the emplacement of subduction. Plate reconstructions allow us to place limits on the location and timing of the continental convergence and subduction zones and the migration of the transition between the two plate boundary regimes. Relative plate motions and reconstructions since the Early to Mid-Miocene require significant continental convergence in advance of the emplacement of the southward migrating Hikurangi subduction—a sequence of tectonism seen in the present plate boundary geography of Hikurangi subduction beneath North Island and convergence in the Southern Alps along the Alpine Fault. In contrast to a transition from subduction to continental convergence where the leading edge of the upper plate is relatively thin and deformable, the transition from a continental convergent regime, with its associated crustal and lithospheric thickening, to subduction of oceanic lithosphere requires substantial thinning (removal) of upper plate continental lithosphere to make room for the slab. The simple structure of the Wadati–Benioff zone seen in the present-day geometry of the subducting Pacific plate beneath North Island indicates that this lithospheric adjustment occurs quickly. Associated with this rapid lithospheric thinning is the development of a series of ephemeral basins, younging to the south, that straddle the migrating slab edge. Based on this association between localized vertical tectonics and slab emplacement, the tectonic history of these basins records the effects of lithospheric delamination driven by the southward migrating leading edge of the subducting Pacific slab. Although the New Zealand plate boundary is often described as simply two subduction zones linked by the transpressive Alpine Fault, in actuality the present is merely a snapshot view of an ongoing and complex evolution from convergence to subduction.     

Geological controls on focused fluid flow through the gas hydrate stability zone on the southern Hikurangi Margin of New Zealand, evidenced from multi-channel seismic data

Highly concentrated gas hydrate deposits are likely to be associated with geological features that promote increased fluid flux through the gas hydrate stability zone (GHSZ). We conduct conventional seismic processing techniques and full-waveform inversion methods on a multi-channel seismic line that was acquired over a 125 km transect of the southern Hikurangi Margin off the eastern coast of New Zealand’s North Island. Initial processing, employed with an emphasis on preservation of true amplitude information, was used to identify three sites where structures and stratal fabrics likely encourage focused fluid flow into and through the GHSZ. At two of the sites, Western Porangahau Trough and Eastern Porangahau Ridge, sub-vertical blanking zones occur in regions of intensely deformed sedimentary layering. It is interpreted that increased fluid flow occurs in these regions and that fluids may dissipate upwards and away from the deformed zone along layers that trend towards the seafloor. At Eastern Porangahau Ridge we also observe a coherent bottom simulating reflection (BSR) that increases markedly in intensity with proximity to the centre of the anticlinal ridge. 1D full-waveform inversions conducted at eight points along the BSR reveal much more pronounced low-velocity zones near the centre of the ridge, indicating a local increase in the flux of gas-charged fluids into the anticline. At another anticline, Western Porangahau Ridge, a dipping high-amplitude feature extends from the BSR upwards towards the seafloor within the regional GHSZ. 1D full-waveform inversions at this site reveal that the dipping feature is characterised by a high-velocity zone overlying a low-velocity zone, which we interpret as gas hydrates overlying free gas. These results support a previous interpretation that this high-amplitude feature represents a local “up-warping” of the base of hydrate stability in response to advective heat flow from upward migrating fluids. These three sites provide examples of geological frameworks that encourage prolific localised fluid flow into the hydrate system where it is likely that gas-charged fluids are converting to highly concentrated hydrate deposits.     

Last glacial aggradation and postglacial sediment production from the non-glacial Waipaoa and Waimata catchments, Hikurangi Margin, North Island, New Zealand

The sediment flux generated by postglacial channel incision has been calculated for the 2150 km², non-glacial, Waipaoa catchment located on the tectonically active Hikurangi Margin, eastern North Island, New Zealand. Sediment production both at a sub-catchment scale and for the Waipaoa catchment as a whole was calculated by first using the tensioned spline method within ARC MAP to create an approximation of the aggradational Waipaoa-1 surface (contemporaneous with the Last Glacial Maximum), and second using grid calculator functions in the GIS to subtract the modern day surface from the Waipaoa-1 surface. The Waipaoa-1 surface was mapped using stereo aerial photography, and global positioning technology fixed the position of individual terrace remnants in the landscape. The recent discovery of Kawakawa Tephra within Waipaoa-1 aggradation gravels in this catchment demonstrates that aggradation was coincidental with or began before the deposition of this 22 600 ¹⁴C-year-old tephra and, using the stratigraphic relationship of Rerewhakaaitu Tephra, the end of aggradation is dated at ca 15 000 ¹⁴C years (ca 18 000 cal. years BP). The construction of the Waipaoa-1 terrace is considered to be synchronous and broadly correlated with aggradation elsewhere in the North Island and northern South Island, indicating that aggradation ended at the same time over a wide area. Subsequent downcutting, a manifestation of base-level lowering following a switch to postglacial incision at the end of glacial-age aggradation, points to a significant Southern Hemisphere climatic warming occurring soon after ca 15 000 ¹⁴C years (ca 18 000 cal. years BP) during the Older Dryas interval. Elevation differences between the Waipaoa-1 (c.15 ka) terrace and the level of maximum channel incision (i.e. before aggradation since the turn of the 20th century) suggest about 50% of the topographic relief within headwater reaches of the Waipaoa catchment has been formed in postglacial times. The postglacial sediment flux generated by channel incision from Waipaoa catchment is of the order of 9.5 km³, of which ~ 6.6 km³ is stored within the confines of the Poverty Bay floodplain. Thus, although the postglacial period represented a time of high terrigenous sediment generation and delivery, only ~ 30% of the sediment generated by channel incision from Waipaoa catchment probably reached the marine shelf and slope of the Hikurangi Margin during this time. The smaller adjacent Waimata catchment probably contributed an additional 2.6 km³ to the same depocentre to give a total postglacial sediment contribution to the shelf and beyond of ~ 5.5 km³. Sediment generated by postglacial channel incision represents only ~ 25% of the total sediment yield from this landscape with ~ 75% of the estimated volume of the postglacial storage offshore probably derived from hillslope erosion processes following base-level fall at times when sediment yield from these catchments exceeded storage.     

Sedimentary architecture of detached deep‐marine canyons: Examples from the East Coast Basin of New Zealand

Submarine canyons are conduits for the distribution of sediment across continental margins. Although many canyons connect directly with fluvial or marine littoral system feeders, canyons detached from direct hinterland supply are also recognized. The fill of detached canyons remains enigmatic, because their deep‐water setting restricts analysis of their evolution and stratigraphic architecture. Therefore, this study aims to investigate the sedimentary processes that infilled deep‐water canyons and the resulting architecture. Miocene outcrops of an exhumed deep‐water system from the East Coast Basin, New Zealand, are documented and compared with the morphology and seismic scale architecture of a modern detached canyon system on the same convergent margin. The outcropping system preserves the downstream margin of a sub‐basin deposited at palaeo‐water depths >700 m. A 6 km wide by 430 m deep incision is filled by heterogeneous siliciclastic sediments, 50% of which comprise graded thin‐beds with traction structures, interpreted to result from oscillatory flows. These are intercalated with concave‐up lenses, up to 15 m thick, of sigmoidally‐bedded, amalgamated sandstones, which preserve ripple casts on bed bases, interpreted as deposits at the head of a deep‐marine canyon. Palaeo‐flow was eastward, into the sub‐basin margin. On the adjacent margin of the sub‐basin down‐dip, stacked and amalgamated sandstones and conglomerates represent the fill of a submarine channel complex, at least 3 km wide. The channels are inferred to have been fed by the up‐dip canyon, which traversed the intervening structural high; similar relationships are seen in the bathymetry data. Seismic studies on this margin demonstrate that multiple phases of canyon cut and fill may occur, with downstream architectural evolution comparable to that seen at outcrop, demonstrating that detached canyons may act as sediment conduits. Breaching of developing sea‐floor structures by detached canyons can modify tortuous sediment pathways, supplying sediment to otherwise starved areas of the slope.