Search

Offshore aquifer research is an emerging field that is becoming increasingly important as population growth and climate change put pressure on coastal water resources. One of the largest reserves, globally, of offshore freshened groundwater (OFG) was recently identified off the South Island of New Zealand. This has highlighted the potential for OFG elsewhere in New Zealand. This study aims to: (1) screen for New Zealand coastal aquifers most likely to contain OFG and, (2) document evidence for OFG in New Zealand. An OFG-likelihood rating scheme was developed as part of the study. An application of the rating scheme used survey responses from regional councils responsible for groundwater management, in combination with national and regional-scale technical documents. The rating scheme was found to be a simple and transparent first-pass approach for highlighting areas where OFG is more or less likely at the national scale. Results are presented in a map showing the likelihood of OFG around the New Zealand coastline. Regions with aquifers where OFG likelihood is high include Greater Wellington, Canterbury, Tasman, Hawkes Bay and Marlborough. Aquifers in these regions are associated with major fluvial depositional systems, including glacial outwash gravels. Despite high dependence on groundwater in these regions and extensive groundwater extraction near the coast, there are no major reported incidences of seawater intrusion, suggesting offshore groundwater may be augmenting onshore extraction.

Shallow seabed depressions attributed to focused fluid seepage, known as pockmarks, have been documented in all continental margins. In this study we demonstrate how pockmark formation can be the result of a combination of multiple factors – fluid type, overpressures, seafloor sediment type, stratigraphy, and bottom currents. We integrate multibeam echosounder and seismic reflection data, sediment cores and pore water samples, with numerical models of groundwater and gas hydrates, from the Canterbury Margin (off New Zealand). More than 6800 surface pockmarks, reaching densities of 100 per km2, and an undefined number of buried pockmarks, are identified in the middle to outer shelf and lower continental slope. Fluid conduits across the shelf and slope include shallow to deep chimneys/pipes. Methane with a biogenic and/or thermogenic origin is the main fluid forming flow and escape features, although saline and freshened groundwaters may also be seeping across the slope. The main drivers of fluid flow and seepage are overpressure across the slope generated by sediment loading and thin sediment overburden above the overpressured interval in the outer shelf. Other processes (e.g. methane generation and flow, a reduction in hydrostatic pressure due to sea-level lowering) may also account for fluid flow and seepage features, particularly across the shelf. Pockmark occurrence coincides with muddy sediments at the seafloor, whereas their planform is elongated by bottom currents.

Turbidite palaeoseismology has produced arguably the most comprehensive multimillennial scale records of subduction-zone earthquakes but is underpinned by techniques that are vigorously debated in earthquake science. Resolving this argument requires new direct observations that test the approach’s essential assumptions. Here we present measurements from turbidites triggered by the 2016 M 7.8 Kaikōura earthquake in New Zealand, one of the most well-instrumented earthquakes in history. This natural experiment provides an ideal test for turbidite palaeoseismology because fault source, ground motions and turbidite deposition in discrete canyons are well-resolved by analysis of sediment cores combined with physics-based ground-motion modelling. We find that the Kaikōura earthquake triggered flows in ten consecutive canyon–distributary systems along a 200 km segment of the Hikurangi subduction margin where long-period (>2 s) peak ground velocities exceeded turbidity-current-triggering thresholds between 16–25 cm s. Comparison between ground motions and turbidite deposition confirm that there is a predictable relationship between earthquake source, ground motions and deposition of coseismic turbidites. We demonstrate that the patterns of triggering and resultant turbidite character may preserve evidence of fault-rupture direction along with the radiating patterns and amplification of earthquake ground motions.

This study investigates the potential for seismic seiches in Lake Tekapo, New Zealand, triggered by ground shaking from an Mw8.2 Alpine Fault earthquake. Synthetic ground motions are used as a forcing boundary to drive lake water motions by further developing a tsunami simulation model—COMCOT—and coupling it with earthquake simulation model outputs. Our modelling results reveal that lake water oscillations are mobilised immediately by the ground movement and further amplified by cross-lake seiches. Amplitudes of lake oscillations can reach up to 4.0 m in the lake’s narrow southern arm, over 1.0 m along the shore of Lake Tekapo township, and about 1.5–2.5 m along many other parts of the lake shore. Large-amplitude water oscillations quickly attenuate in the first 5–10 min after the earthquake due to their relatively short periods, while long-period oscillations continue for a long time, albeit with much smaller amplitudes. Spectral analysis clearly reveals that the ground motions trigger both fundamental and higher modes in the lake whose oscillation periods are consistent with theoretical estimates. We find that large-amplitude lake water oscillations are better correlated with low-frequency, less energetic ground motion content than with high-frequency large-amplitude ground motions. Ground motion-triggered lake oscillations are large enough to pose potential threat to tourists, residents, boats and infrastructure both in the lake water and onshore near the waterfront. In contrast, vertical co-seismic displacements in the lake area, the conventional mechanism for tsunami generation, are too small to trigger tsunami waves of concern.

The initial stages of seamount subduction and associated deformation in an overriding accretionary wedge is rarely documented. Initial subduction of Bennett Knoll seamount and faulting of the overlying strata along the Hikurangi subduction margin, New Zealand, are here studied using multibeam swath bathymetry, subbottom profiles, and regional seismic reflection lines. Our results provide new insights into the earliest stages of seamount collision at sediment-rich margins. Differential shortening along the subduction front induced by seamount subduction is initially accommodated in the accretionary wedge by conjugate strike-slip faults that straddle the buried flanks of the seamount and offset the frontal thrusts by as much as 5 km. The geometries of the strike-slip faults are controlled by the seamount’s dimensions and aspect, the obliquity of plate convergence, pore-fluid pressure, and the thickness and rheology of the incoming sedimentary section. Strike-slip faults in such settings are ephemeral and overprinted by the formation of new structures as seamount subduction advances.

Recent earthquakes involving multiple fault ruptures highlight the need to evaluate complex coastal deformation mechanisms, which are important for understanding plate boundary kinematics and seismic and tsunami hazards. We compare ages and uplift of the youngest Holocene marine terraces at Puatai Beach and Pakarae River mouth (∼10 km apart) in the northern Hikurangi subduction margin to examine whether uplift is the result of subduction earthquakes or upper-plate fault earthquakes. From stepped platform-cliff morphology, we infer uplift during 2–3 earthquakes and calculate an average uplift-per-event of 2.9±0.5 m at Puatai Beach and 2.0±0.5 m at Pakarae River mouth. Radiocarbon ages from the youngest beach deposit shells on each terrace and a tephra coverbed on one terrace constrain the timing of earthquakes to 1770–1710, 1100–910, and 420–250 cal. B.P. at Puatai Beach, and 1490–1290 and 660–530 cal. B.P. at Pakarae River mouth. The ages differ at each site indicating uplift is neither the result of subduction earthquakes nor single upper-plate fault earthquakes. A reinterpretation of new and existing bathymetry and seismic reflection data, combined with dislocation modeling, indicates that near-shore fault segmentation is more complex than previously thought and ruptures likely involve multiple upper-plate faults. Future updates of the New Zealand National Seismic Hazard Model should revise the northern Hikurangi subduction seismic sources so that rupture does not uplift Puatai Beach and Pakarae River mouth and include new near-shore upper-plate faults as multifault sources.

The Tuaheni Landslide Complex (TLC) is characterised by areas of compression upslope and extension downslope. It has been thought to consist of a stack of two genetically linked landslide units identified on seismic data. We use 3D seismic reflection, bathymetry data, and IODP core U1517C (Expedition 372), to understand the internal structures, deformation mechanisms and depositional processes of the TLC deposits. Unit II and Unit III of U1517C correspond to the two chaotic units in 3D seismic data. In the core, Unit II shows deformation whereas Unit III appears more like an in situ sequence. Variance attribute analysis shows that Unit II is split in lobes around a coherent stratified central ridge and is bounded by scarps. By contrast, we find that Unit III is continuous beneath the central ridge and has an upslope geometry that we interpret as a channellevee system. Both units show evidence of lateral spreading due to the presence of the Tuaheni Canyon removing support from the toe. Our results suggest that Unit II and Unit III are not genetically linked, that they are separated substantially in time and they had different emplacement mechanisms, but fail under similar circumstances.

Submarine landslides occur on continental margins globally and can have devastating consequences for marine habitats, offshore infrastructure and coastal communities due to potential tsunamigenic consequences. Evaluation of the magnitude and distribution of submarine landslides is central to marine and coastal hazard planning. Despite this, there are few studies that comprehensively quantify the occurrence of submarine landslides on a margin-wide scale. We present the first margin-wide submarine landslide database along the eastern margin of New Zealand comprising >2200 landslide scars and associated mass-transport deposits. Analysis of submarine landslide distribution reveals 1) locations prone to mass-failure, 2) spatial patterns of landslide scale and occurrence, and 3) the potential preconditioning factors and triggers of mass wasting across different geologic settings. Submarine landslides are widespread on the eastern margin of New Zealand, occurring in water depths from ~300 m to ~4,000 m. Landslide scars and mass transport deposits are more prevalent, and on average larger, on the active margin, compared the passive margin. We attribute higher concentrations of landslides on the active margin to the prevalence of deforming thrust ridges, related to active margin processes including oversteepening, faulting and seamount subduction. Higher sediment supply on the northernmost active margin is also likely to be a key preconditioning factor resulting in the concentration of large landslides in this region. In general, submarine landslide scars are concentrated around canyon systems and close to canyon thalwegs. This suggests that not only does mass wasting play a major role in canyon evolution, but also that slope undercutting in canyons may be a fundamental preconditioning factor for slope failure. Results of this study offer unique insights into the spatial distribution, magnitude and morphology of submarine landslides across different geologic settings, providing a better understanding of the causative factors for mass wasting in New Zealand and around the world.

17 pages, 9 figures, 1 table, Submarine canyons are deep, large-scale incisions that occur on the continental shelf and slope of all ocean margins. These landforms serve as preferential particle-transport conduits that connect the coastal zone with the deep-sea. Canyons have been studied for decades and are among the most iconic submarine geomorphic features. Advances in marine technology are increasingly providing new understanding of the physical and biological processes that take place in canyons, which are also notoriously affected by the impact of anthropogenic activities. In this article we describe the main geomorphic characteristics of submarine canyons and synthesize the past and current knowledge on their formation, evolution, global distribution and hydrodynamical and ecological significance.

Focused gas migration through the gas hydrate stability zone in vertical gas conduits is a global phenomenon. The process can lead to concentrated gas hydrate formation and seafloor gas seepage, which influences seafloor biodiversity and ocean biogeochemistry. However, much is unknown about how gas and gas hydrate co‐exist within and around gas conduits. We present seismic imaging of the gas hydrate system beneath a four‐way closure anticlinal ridge at New Zealand's southern Hikurangi subduction margin. Gas has accumulated beneath the base of gas hydrate stability to a thickness of up to ∼240 m, which has ultimately led to hydraulic fracturing and propagation of a vertical gas conduit to the seafloor. Despite the existence of an array of normal faults beneath the ridge, these structures are not exploited as long‐range gas flow conduits. Directly beneath the conduit, and extending upward from the regional base of gas hydrate stability, is a broad zone characterized by both negative‐ and positive‐polarity reflections. We interpret this zone as a volume of sediment hosting both gas hydrate and free gas, that developed due to partial gas trapping beneath a mass transport deposit. Similar highly reflective zones have been identified at the bases of other gas conduits, but they are not intrinsic to all gas conduits through gas hydrate systems. We suggest that pronounced intervening sealing units within the gas hydrate stability zone determine whether or not they form., Plain Language Summary: Gas hydrates are ice‐like substances composed of natural gas and water. They form between sediment grains underneath large regions of the Earth's seafloor. An important reason to study gas hydrates is that they partly control the way that methane gas flows through sediments and out of the seafloor. It is this flow of methane that sustains some diverse biological communities on the seafloor and affects the chemistry of the oceans. In this study, we use reflected sound waves to explore how gas flow beneath the seafloor depends on the way in which sedimentary layers are folded and fractured. Our data reveal a 240‐m thick reservoir of gas that is trapped in a large sedimentary fold, ∼500 m beneath the seafloor. The buoyancy of the gas has caused a vertical fracture zone to propagate upward to the seafloor, where gas bubbles are venting into the ocean. Further, our data suggest that a broad accumulation of gas hydrates (together with gas) has formed beneath the vertical fracture zone. This gas hydrate deposit may grow larger with time, and it will continue to influence the way that gas flows through the sediments., Key Points: 240 m thick free gas column accumulated beneath the base of hydrate stability, which led to hydraulic fracturing and gas chimney formation. Broad zone of free gas and gas hydrate formed beneath the gas chimney, extending upward from the regional base of hydrate stability. Such zones of hydrate and free gas likely form due to pronounced lithological contrasts (sealing layers) within the hydrate stability zone., New Zealand's Ministry for Business Innovation and Employment (MBIE) http://dx.doi.org/10.13039/501100004629

Subaqueous spreading is a widespread type of mass movement, which involves extensional displacement along a gliding plane and the deformation of the failing layer into a sequence of ridges and troughs. Spreading has been poorly investigated, nonetheless it poses hazard to offshore infrastructures. SubSpread is a new project that will investigate the mechanics of the spreading failure and its geological controls in the subaqueous environment. The first objective of SubSpread is to identify the topographic and sedimentary signature of subaqueous environment. We have compiled a global database of subaqueous and subaerial spreads that includes information on physiography, geomorphology, sedimentology and geotechnical properties, where available. A preliminary analysis of the database reveals that spreading morphologies occur on both passive and active margins, especially in the headwall area of translational retrogressive slides. Potential causes of spreading include seismic loading (also glacially induced), sediment loading, and increased pore pressure generated by migration of fluid or gas. The latter may induce loss of shear strength and the formation of a weak layer, particularly in gentle open slopes. Information compiled in this database will also be used to develop a numerical model that can better understand the mechanics and rheological aspects of submarine spreading, focusing on the role played by pore pressure generation. The Tuaheni slide complex in the Hikurangi Margin of New Zealand is being used as a case-study in view of the wealth of geophysical and sedimentological data that are available. The final part of the SubSpread project will test whether the morphometric and sedimentological signature of spreading can provide information on past seismicity. In this case, the test site will be Lake Tekapo in the South Island of New Zealand., peer-reviewed

Highlights • Recently acquired high-resolution seismic data and existing low-resolution industry data are presented. • Two large concentrated hydrate deposits are identified beneath Glendhu and Honeycomb ridges. • A novel method involving analysis of seismic velocity and reflectivity is used to obtain estimates of hydrate saturations. • Hydrate saturations peaks of >80% are estimated locally. • The main driving mechanism for hydrate accumulations is inferred to be along-strata gas migration. Abstract In the southern Hikurangi subduction margin, widespread gas hydrate accumulations are inferred based on the presence of bottom simulating reflections and recovered gas hydrate samples, mainly associated with thrust ridges. We present a detailed analysis of high- and medium-resolution seismic reflection data across Glendhu and Honeycomb ridges, two elongated four-way closure systems at the toe of the deformation wedge. High-amplitude reflections within the gas hydrate stability zone, coincident with high seismic velocities, suggest the presence of highly concentrated gas hydrate accumulations in the core regions of the anticlinal ridges. A novel method involving combined seismic velocity and reflectivity analysis and rock physics modelling is used to estimate hydrate saturations in localised areas. The effective medium model consistently predicts gas hydrate saturations of ~30% of the pore space at Glendhu Ridge and >60% at Honeycomb Ridge, whereas the empirical three-phases weighted equation likely underestimates the amount of gas hydrate present. We note that our estimates are dependent on the vertical resolution of the seismic data (5–14 m), and that the existence of thin layers hosting gas hydrate at higher concentrations is likely based on observations made elsewhere in similar depositional environments. A comparison between the two ridges provides insights into the evolution of thrust related anticlines at the toe of the accretionary wedge. We propose that the main driving mechanism for concentrated hydrate accumulation in the study area is along-strata gas migration. The vertical extent of these accumulations is a function of the steepness of the strata crossing the base of gas hydrate stability, and of the volume of sediments from which fluid flows into each structure. According to our interpretation, older structures situated further landward ofthe deformation front are more likely to host more extensive concentrated hydrate deposits than younger ridges situated at the deformation front and characterised by more gentle folding. The method introduced in this work is useful to retrieve quantitative estimates of gas hydrate saturations based on multi-channel seismic data.

The distribution of benthic ecosystems, dominated by filter-feeding communities, is highly influenced by the seabed geomorphology. However, the spatial variation in settlement of these species is also affected by near-bottom currents and any changes in light, nutrient concentration and food quality often associated with increases of suspended sediment concentrations within the water column. Detailed predictions of the geographic distribution of filter-feeder species and a deeper understanding of the physical processes influencing their distribution patterns is key for effective management and conservation. To date, predictive distribution modelling has been derived essentially from geomorphological parameters, mainly using spatially limited observations. In this study, seabed mapping, oceanographic modelling, hydrographic records and biological observations are integrated to provide high-resolution prediction of filter-feeder habitat distribution within Queen Charlotte Sound/Tōtaranui and Tory Channel/Kura Te Au, South Island of New Zealand. The aim is to evaluate potential suitable habitat areas for filter-feeders to inform where habitat restoration management should focus efforts to recover communities such as the horse mussel (Atrina zelandica) or the green-lipped mussel (Perna canaliculus), both of which have high economic impact in New Zealand. To accomplish this, Maximum Entropy (MaxEnt) predictive modelling was used to produce Habitat Suitability (HS) maps, using geomorphological parameters and seafloor classification information. Final HS maps also incorporated oceanographic and sediment dynamic information, showing that filter-feeder habitat distribution is highly influenced by the hydrodynamics and sedimentary processes apart from the seafloor geomorphology. Filter-feeder communities inhabit quiescent areas, limited by depth, slope and sediment type; and coincide with regions presenting low near-bottom currents and low turbidity levels. Additionally, the obtained results reveal the effects of the coastal settlements and major marine traffic routes, limiting the suitable habitats to areas with less human impact. This study demonstrates that a multidisciplinary approach is crucial to better predict the spatial distribution of benthic communities, which is key to improve benthic habitat restoration and recovery assessments.

Studies have shown the importance of submarine canyons as conduits of land-derived organic carbon beyond the coastal shelf into the deep-sea where a single obvious river source can be identified. When there is more than one river source, identifying which rivers contribute to canyon sediment organic matter is technically challenging. Here, we compare two contrasting submarine canyons: the Hokitika Canyon, a long, narrow, and gently sloping canyon on the west coast of New Zealand; and the Kaikōura Canyon, a high productivity, short, steep canyon close to shore on the east coast of New Zealand. Both canyons have multiple potential river sources, so we applied a compound specific stable isotope (CSSI) tracking technique to identify and apportion the contribution from each river at locations along the length of each canyon axis. We found that land-derived organic matter contributed between 74 and 100% of the total organic matter in the sediment of the Hokitika Canyon as far as 200 km from shore and to depths of 2000 m. However, less than 50% of the land-derived organic carbon came from the largest river closest to the canyon head. We hypothesize that longshore drift transported much of the sediment from that river past the Hokitika Canyon, while river inflows farther up-current supplied the bulk of the land-derived organic carbon. In contrast, land-derived organic matter contributed less than 50% of the total organic matter in Kaikōura Canyon sediments with land-derived organic sediment contribution decreasing steeply to less than 15% at about 24 km from shore in 1500 m water depth. Most of the land-derived organic matter (ca. 80%) came from the river with the largest suspended sediment yield, despite another (smaller) river discharging closer to the canyon head. We hypothesize that this difference in carbon source is partly due to the comparatively short and steep, and therefore dynamic, nature of Kaikōura Canyon resulting in efficient sediment through-put. The efficiency with which organic matter is captured and transferred to the deep-sea by canyons demonstrates the potential for such systems to act as natural carbon sinks driven by both geologically episodic and more regular oceanographic processes.

The consequences of subaqueous landslides have been at the forefront of societal conscience more than ever in the last few years, with devastating and fatal events in the Indonesian Archipelago making global news. The new research presented in this volume demonstrates the breadth of ongoing investigation into subaqueous landslides, and shows that while events like the recent ones can be devastating, they are smaller in scale than those Earth has experienced in the past. Understanding the spectrum of subaqueous landslide processes, and therefore the potential societal impact, requires research across all spatial and temporal scales. This volume delivers a compilation of state-of-the-art papers covering regional landslide databases, advanced techniques for in situ measurements, numerical modelling of processes and hazards

Offshore freshened groundwater (OFG) is an important component of the global water cycle. Bodies of fresh and moderately brackish groundwater have been documented up to 100 km from modern shorelines. Their global volumetric estimates are on the order of 500,000 cubic kilometres, which is two orders of magnitude greater than the volume of groundwater extracted globally from continental aquifers since 1900. The potential use of OFG systems as a source of potable water is a main driving force for their improved understanding. OFG systems also play a fundamental role in biogeochemical fluxes to the ocean and in benthic and sub-seafloor ecology., peer-reviewed

Numerous subaqueous landslides exhibit spreading failure morphologies which are typically characterized by repetitive patterns of parallel ridges and troughs oriented perpendicular to the direction of movement. Whilst these spreading failures are commonly attributed to (i) downslope removal of material causing unloading of the temporary stable slope or (ii) significant loss of shear strength of the substratum allowing blocks of overlying sediment to detach and slide downslope, their movement rates and potential triggers remain poorly constrained. Spreading appears to be a dominant failure mechanism within the Tuaheni Landslide Complex (TLC) on the Hikurangi Subduction Margin off the coast of Gisborne, New Zealand. A combination of swath bathymetric, 2D and 3D seismic data, drilling investigations and laboratory experiments on sediments recovered from the TLC indicate that this geomorphology has been generated by translational failure. Failure could occur through episodic cycles of movement-arrest in response to either elevated pore fluid pressures or undrained loading during earthquakes. We developed numerical models that integrate this unique data set to explore the processes that lead to spreading failure and determine how large shear strains can be accommodated without accelerating to catastrophic failure. The results provide a novel approach that demonstrates how seafloor morphology can, in part, be controlled by the underlying failure processes

Highlights • Seismic depth imaging gives insight into the southern Hikurangi subduction zone. • Velocities reveal regional variations in compaction and drainage of input sediments. • Dewatering of subducted sediments might influence décollement strength. • Thrusts at the leading edge of deformation are upper-plate dewatering pathways. • Stratigraphic host of the décollement changes at the southern end of the margin. Abstract The southern end of New Zealand's Hikurangi subduction margin accommodates highly oblique convergence between the Pacific and Australian plates. We carry out two-dimensional (2D) seismic reflection tomography and pre-stack depth migrations on two seismic lines to gain insight into the nature of subducted sediments and upper plate faulting and dewatering at the toe of the wedge. We also investigate the NE to SW evolution of emergent upper plate thrust faulting using 47 seismic lines spanning an along-strike distance of ∼270 km. The upper sequence of sediments that ultimately gets subducted (the MES sequence) has an anomalously-low seismic velocity character. At the southwestern end of the margin, ∼150 km east of Kaikōura, the MES sequence has experienced greater compaction (for an equivalent effective vertical stress) than it has some 200 km further to the northeast. This difference is likely attributable to greater horizontal compression in the southwest caused by impingement of the Chatham Rise on the deformation front. Relationships between velocity and effective vertical stress suggest that the MES sequence is well-drained in the vicinity of frontal thrusts, corroborated by evidence for upper plate dewatering along those thrusts. Effective drainage of the MES sequence likely promotes interplate coupling on the southern Hikurangi margin. The décollement is generally hosted near a seismic reflector known as “Reflector 7”. East of Kaikōura, however, Reflector 7 becomes accreted, indicating that subduction slip at the southwestern end of the margin is no longer hosted at (or above) this reflector. Instead, the décollement steps down to a deeper stratigraphic level further inboard. Further to the SW, approximately in line with the lower Kaikōura Canyon, the offshore manifestation of subduction-driven compression ceases.

Hydrogeological processes influence the morphology, mechanical behavior, and evolution of subduction margins. Fluid supply, release, migration, and drainage control fluid pressure and collectively govern the stress state, which varies between accretionary and nonaccretionary systems. We compiled over a decade of published and unpublished acoustic data sets and seafloor observations to analyze the distribution of focused fluid expulsion along the Hikurangi margin, New Zealand. The spatial coverage and quality of our data are exceptional for subduction margins globally. We found that focused fluid seepage is widespread and varies south to north with changes in subduction setting, including: wedge morphology, convergence rate, seafloor roughness, and sediment thickness on the incoming Pacific plate. Overall, focused seepage manifests most commonly above the deforming backstop, is common on thrust ridges, and is largely absent from the frontal wedge despite ubiquitous hydrate occurrences. Focused seepage distribution may reflect spatial differences in shallow permeability architecture, while diffusive fluid flow and seepage at scales below detection limits are also likely. From the spatial coincidence of fluids with major thrust faults that disrupt gas hydrate stability, we surmise that focused seepage distribution may also reflect deeper drainage of the forearc, with implications for pore-pressure regime, fault mechanics, and critical wedge stability and morphology. Because a range of subduction styles is represented by 800 km of along-strike variability, our results may have implications for understanding subduction fluid flow and seepage globally., peer-reviewed

Screaton et al. (2019, https://doi.org/10.1029/2019GC008603) examined the role of sedimentation, sea level, and bottom water temperature (BWT) changes due to glaciation as drivers for the downward migration of the base of gas hydrate stability and gas hydrate formation. International Ocean Discovery Program (IODP) Site U1517 in the Hikurangi margin was used as a case study because data at this site document a marked increase in chloride over a broad depth range, which was attributed to recent gas hydrate formation. In a comment on Screaton et al. (2019, https://doi.org/10.1029/2019GC008603), Sultan (2020, https://doi.org/10.1029/2019gc008846) used a linear thermal profile to argue that inferences and characterization of methane hydrate at IODP Site U1517 were incorrect because some occur below his estimated base of gas hydrate stability (BGHS). Based on this apparent discrepancy, Sultan (2020, https://doi.org/10.1029/2019gc008846) further stated that low‐chloride spikes may be unreliable indicators of methane hydrate occurrence. In this reply, we emphasize that unsteady‐state, and thus nonlinear, thermal profiles are likely in areas experiencing active sedimentation and bottom‐water temperature (BWT) changes. The resulting deviation from steady‐state temperature profile shifts the BGHS downward. In addition, sedimentation has the potential to bury methane hydrate more rapidly than it dissociates, helping to explain how methane hydrate could be observed below the BGHS. We also review the supporting evidence for gas‐hydrate occurrence at Site U1517 and the criteria used for Site U1517 site selection.

Multiple (>20 \ud >20\ud ) crustal faults ruptured to the ground surface and seafloor in the 14 November 2016 M w \ud Mw\ud 7.8 Kaikōura earthquake, and many have been documented in detail, providing an opportunity to understand the factors controlling multifault ruptures, including the role of the subduction interface. We present a summary of the surface ruptures, as well as previous knowledge including paleoseismic data, and use these data and a 3D geological model to calculate cumulative geological moment magnitudes (M G w \ud MwG\ud ) and seismic moments for comparison with those from geophysical datasets. The earthquake ruptured faults with a wide range of orientations, sense of movement, slip rates, and recurrence intervals, and crossed a tectonic domain boundary, the Hope fault. The maximum net surface displacement was ∼12 m \ud ∼12 m\ud on the Kekerengu and the Papatea faults, and average displacements for the major faults were 0.7–1.5 m south of the Hope fault, and 5.5–6.4 m to the north. M G w \ud MwG\ud using two different methods are M G w \ud MwG\ud 7.7 +0.3 −0.2 \ud 7.7−0.2+0.3\ud and the seismic moment is 33%–67% of geophysical datasets. However, these are minimum values and a best estimate M G w \ud MwG\ud incorporating probable larger slip at depth, a 20 km seismogenic depth, and likely listric geometry is M G w \ud MwG\ud 7.8±0.2 \ud 7.8±0.2\ud , suggests ≤32% \ud ≤32%\ud of the moment may be attributed to slip on the subduction interface and/or a midcrustal detachment. Likely factors contributing to multifault rupture in the Kaikōura earthquake include (1) the presence of the subduction interface, (2) physical linkages between faults, (3) rupture of geologically immature faults in the south, and (4) inherited geological structure. The estimated recurrence interval for the Kaikōura earthquake is ≥5,000–10,000 yrs \ud ≥5,000–10,000 yrs\ud , and so it is a relatively rare event. Nevertheless, these findings support the need for continued advances in seismic hazard modeling to ensure that they incorporate multifault ruptures that cross tectonic domain boundaries.

In subduction zones upper-plate normal faults have long been considered a tectonic feature primarily associated with erosive margins. However, increasing data coverage has proven that similar features also occur in accretionary margins, such as Cascadia, Makran, Nankai or Central Chile, where kinematics are dominated by compression. Considering their wide distribution there is, without doubt, a significant lack of qualitative and quantitative knowledge regarding the role and importance of normal faults and zones of extension for the seismotectonic evolution of accretionary margins. We use a high-resolution 3D P-Cable seismic volume from the Hikurangi Margin acquired in 2014 to analyze the spatial distribution and mechanisms of upper-plate normal faulting. The study area is located at the upper continental slope in the area of the Tuaheni landslide complex. In detail we aim to (1) map the spatial distribution of normal faults and characterize their vertical throws, strike directions, and dip angles; (2) investigate their possible influence on fluid migration in an area, where gas hydrates are present; (3) discuss the mechanisms that may cause extension of the upper-slope in the study area. Beneath the Tuaheni Landslide Complex we mapped about 200 normal faults. All faults have low displacements ( 65°) angles. About 71% of the faults dip landward. We found two main strike directions, with the majority of faults striking 350-10°, parallel to the deformation front. A second group of faults strikes 40-60°. The faults crosscut the BSR, which indicates the base of the gas hydrate zone. In combination with seismically imaged bright-spots and pull-up structures, this indicates that the normal faults effectively transport fluids vertically across the base of the gas hydrate zone. Localized uplift, as indicated by the presence of the Tuaheni Ridge, might support normal faulting in the study area. In addition, different subduction rates across the margin may also favor extension between the segments. Future work will help to further untangle the mechanisms that cause extension of the upper continental slope.

We investigate gas hydrate system dynamics beneath a submarine canyon on New Zealand's Hikurangi subduction margin using seismic reflection data and petroleum systems modelling. High seismic velocities just above the base of gas hydrate stability (BGHS) indicate that concentrated gas hydrates exist beneath the canyon. Two-dimensional gas hydrate formation modelling shows how the process of canyon incision at this location alters the distribution and concentration of gas hydrate. The key modelling result is that free gas is trapped beneath the gas hydrate layer and then ‘captured' into a concentrated gas hydrate deposit as a result of a downward-shift in the BGHS driven by canyon incision. Our study thus provides new insight into the functioning of this process. From our data, we also conceptualise two other models to describe how canyons could significantly change gas hydrate distribution and concentration. One scenario is related to deflection of fluid flow pathways from over-pressured regions at the BGHS toward the canyon, and the other is based on relationships between simultaneous seafloor uplift and canyon incision. The relationships and processes described are of global relevance because of considerations of gas hydrate as an energy resource and the influence of both submarine canyons and gas hydrate systems on seafloor biodiversity.

The MARCAN project, launched last January, is working to fill a gap in our knowledge of how freshwater flowing underground shapes and alters the continental margins.

Coseismic coastal deformation is often used to understand slip on offshore faults in large earthquakes but in the 2016 M W 7.8 Kaikōura earthquake multiple faults ruptured across and sub-parallel to the coastline. Along ∼110 km of coastline, a rich dataset of coastal deformation comprising airborne lidar differencing, field surveying and satellite geodesy reveals highly variable vertical displacements, ranging from −2.5 to 6.5 m. These inform a refined slip model for the Kaikōura earthquake which incorporates changes to the slip on offshore faults and inclusion of an offshore reverse crustal fault that accounts for broad, low-amplitude uplift centered on Kaikōura Peninsula. The exceptional detail afforded by differential lidar and the high variability in coastal deformation combine to form the highest-resolution and most complex record of coseismic coastal deformation yet documented. This should prompt reassessment of coastal paleoseismic records that may not have considered multi-fault ruptures and high complexity deformation fields.

We provide a summary of the surface fault ruptures produced by the Mw7.8 14 November 2016 Kaikōura earthquake, including examples of damage to engineered structures, transportation networks and farming infrastructure produced by direct fault surface rupture displacement. We also provide an overview of the earthquake in the context of the earthquake source model and estimated ground motions from the current (2010) version of the National Seismic Hazard Model (NSHM) for New Zealand. A total of 21 faults ruptured along a c.180 km long zone during the earthquake, including some that were unknown prior to the event. The 2010 version of the NSHM had considered multi-fault ruptures in the Kaikōura area, but not to the degree observed in the earthquake. The number of faults involved a combination of known and unknown faults, a mix of complete and partial ruptures of the known faults, and the non-involvement of a major fault within the rupture zone (i.e. the Hope Fault) makes this rupture an unusually complex event by world standards. However, the strong ground motions of the earthquake are consistent with the high hazard of the Kaikōura area shown in maps produced from the NSHM.

The 2016 M w 7.8 Kaikōura earthquake continued a notable decade of damaging earthquake impacts in New Zealand. The effects were wide ranging across the upper South Island, and included two fatalities, tsunami, tens of thousands of landslides, the collapse of one residential building, and damage to numerous structures and infrastructure. We present a preliminary overview focused on the seismological aspects of this earthquake and the corresponding seismological response effort. The earthquake rupture was extremely complex, involving at least 13 separate faults extending over ∼150 km, from the epicenter in north Canterbury to near the Cook Strait. We use backprojection and slip inversion methods to derive preliminary insights into the rupture evolution, identifying south‐to‐north rupture, including at least three distinct southwest‐to‐northeast propagating phases. The last phase is associated with a strong second pulse of energy release in the northern half of the rupture zone ∼70 s after rupture initiation, which we associate with the Kekerengu–Needles faults where some of the largest surface displacements (dextral) were observed. The mechanism of the mainshock was oblique thrust and relocated aftershocks show a range of thrust and strike‐slip mechanisms across three dominant spatial clusters. GeoNet datasets collected during the Kaikōura earthquake will be crucial in further unraveling details of the complex earthquake rupture and its implications for seismic hazard. Ground motions during the earthquake exceeded 1 g at both ends of the rupture. Spectral accelerations exceeded 500‐year return period design level spectra in numerous towns in the upper South Island, as well as in parts of the capital city of Wellington at critical periods of 1–2 s, influenced by site/basin and directivity effects. Another important part of the response effort has been the provision of earthquake forecasts, as well as consideration of the implications of slow slip on the Hikurangi subduction interface triggered as a result of the Kaikōura earthquake.

Slow slip events (SSEs) at the northern Hikurangi subduction margin, New Zealand, are among the best-documented shallow SSEs on Earth. International Ocean Discovery Program Expeditions 372 and 375 were undertaken to investigate the processes and in situ conditions that underlie subduction zone SSEs at the northern Hikurangi Trough. We accomplished this goal by (1) coring and geophysical logging at four sites, including penetration of an active thrust fault (the Pāpaku fault) near the deformation front, the upper plate above the SSE source region, and the incoming sedimentary succession in the Hikurangi Trough and atop the Tūranganui Knoll seamount; and (2) installing borehole observatories in the Pāpaku fault and in the upper plate overlying the slow slip source region. Logging-while-drilling (LWD) data for this project were acquired as part of Expedition 372, and coring, wireline logging, and observatory installations were conducted during Expedition 375. Northern Hikurangi subduction margin SSEs recur every 1–2 y and thus provide an ideal opportunity to monitor deformation and associated changes in chemical and physical properties throughout the slow slip cycle. In situ measurements and sampling of material from the sedimentary section and oceanic basement of the subducting plate reveal the rock properties, composition, lithology, and structural character of material that is transported downdip into the SSE source region. A recent seafloor geodetic experiment raises the possibility that SSEs at northern Hikurangi may propagate to the trench, indicating that the shallow thrust fault (the Pāpaku fault) targeted during Expeditions 372 and 375 may also lie in the SSE rupture area and host a portion of the slip in these events. Hence, sampling and logging at this location provides insights into the composition, physical properties, and architecture of a shallow fault that may host slow slip. Expeditions 372 and 375 were designed to address three fundamental scientific objectives: Characterize the state and composition of the incoming plate and shallow fault near the trench, which comprise the protolith and initial conditions for fault zone rock at greater depth and which may itself host shallow slow slip; Characterize material properties, thermal regime, and stress conditions in the upper plate directly above the SSE source region; and Install observatories in the Pāpaku fault near the deformation front and in the upper plate above the SSE source to measure temporal variations in deformation, temperature, and fluid flow. The observatories will monitor volumetric strain (via pore pressure as a proxy) and the evolution of physical, hydrological, and chemical properties throughout the SSE cycle. Together, the coring, logging, and observatory data will test a suite of hypotheses about the fundamental mechanics and behavior of SSEs and their relationship to great earthquakes along the subduction interface.

Cook Strait Canyon is a submarine canyon that lies within ten kilometres of Wellington, the capital city of New Zealand. The canyon walls are covered with scars from previous landslides which could have caused local tsunamis. Palaeotsunami evidence also points to past tsunamis in the Wellington region. Furthermore, the canyon’s location in Cook Strait means that there is inhabited land in the path of both forward- and backward-propagating waves. Tsunamis induced by these submarine landslides pose hazard to coastal communities and infrastructure but major events are very uncommon and the historical record is not extensive enough to quantify this hazard. The combination of infrequent but potentially very consequential events makes realistic assessment of the hazard challenging. However, information on both magnitude and frequency is very important for land use planning and civil defence purposes. We use a multidisciplinary approach bringing together geological information with modelling to construct a Probabilistic Tsunami Hazard Assessment of submarine landslide-generated tsunami. Although there are many simplifying assumptions used in this assessment, it suggests that the Cook Strait open coast is exposed to considerable hazard due to submarine landslide-generated tsunamis. We emphasise the uncertainties involved and present opportunities for future research.

The existence of free gas and gas hydrate in the pore spaces of marine sediments causes changes in acoustic velocities that overprint the background lithological velocities of the sediments themselves. Much previous work has determined that such velocity overprinting, if sufficiently pronounced, can be resolved with conventional velocity analysis from long-offset, multichannel seismic data. We used 2D seismic data from a gas hydrate province at the southern end of New Zealand’s Hikurangi subduction margin to describe a workflow for high-resolution velocity analysis that delivered detailed velocity models of shallow marine sediments and their coincident gas hydrate systems. The results showed examples of pronounced low-velocity zones caused by free gas ponding beneath the hydrate layer, as well as high-velocity zones related to gas hydrate deposits. For the seismic interpreter of a gas hydrate system, the velocity results represent an extra “layer” for interpretation that provides important information about the distribution of free gas and gas hydrate. By combining the velocity information from the seismic transect with geologic samples of the seafloor and an understanding of sedimentary processes, we have determined that high gas hydrate concentrations preferentially form within coarse-grained sediments at the proximal end of the Hikurangi Channel. Finer grained sediments expected elsewhere along the seismic transect might preclude the deposition of similarly high gas hydrate concentrations away from the channel.

Highlights • Release of dissolved Sr2+ with low 87Sr/86Sr, as well as Ca2+ and Ba2+ suggests ongoing volcanic ash alteration. • A concurrent increase in Fe2+ and a depletion of CH4 with a decrease in C of CH4 and DIC suggest Fe-AOM. • We for the first time document the potential linkage between ash alteration and methane oxidation via Fe-AOM. • The rate of Fe-AOM is estimated to be ∼0.4 μmol cm−2 yr−1, equivalent to ∼12% of total CH4 removal. Abstract We present geochemical data collected from volcanic ash-bearing sediments on the upper slope of the northern Hikurangi margin during the RV SONNE SO247 expedition in 2016. Gravity coring and seafloor drilling with the MARUM-MeBo200 allowed for collection of sediments down to 105 meters below seafloor (mbsf). Release of dissolved Sr2+ with isotopic composition enriched in 86Sr (87Sr/86Sr minimum = 0.708461 at 83.5 mbsf) is indicative of ash alteration. This reaction releases other cations in the 30-70 mbsf depth interval as reflected by maxima in pore-water Ca2+ and Ba2+ concentrations. In addition, we posit that Fe(III) in volcanogenic glass serves as an electron acceptor for methane oxidation, a reaction that releases Fe2+ measured in the pore fluids to a maximum concentration of 184 μM. Several lines of evidence support our proposed coupling of ash alteration with Fe-mediated anaerobic oxidation of methane (Fe-AOM) beneath the sulfate-methane transition (SMT), which lies at ∼7 mbsf at this site. In the ∼30-70 mbsf interval, we observe a concurrent increase in Fe2+ and a depletion of CH4 with a well-defined decrease in C-CH4 values indicative of microbial fractionation of carbon. The negative excursions in C values of both DIC and CH4 are similar to that observed by sulfate-driven AOM at low SO concentrations, and can only be explained by the microbially-mediated carbon isotope equilibration between CH4 and DIC. Mass balance considerations reveal that the iron cycled through the coupled ash alteration and AOM reactions is consumed as authigenic Fe-bearing minerals. This iron sink term derived from the mass balance is consistent with the amount of iron present as carbonate minerals, as estimated from sequential extraction analyses. Using a numerical modeling approach we estimate the rate of Fe-AOM to be on the order of 0.4 μmol cm−2 yr−1, which accounts for ∼12% of total CH4 removal in the sediments. Although not without uncertainties, the results presented reveal that Fe-AOM in ash-bearing sediments is significantly lower than the sulfate-driven CH4 consumption, which at this site is 3.0 μmol cm−2 yr−1. We highlight that Fe(III) in ash can potentially serve as an electron acceptor for methane oxidation in sulfate-depleted settings. This is relevant to our understanding of C-Fe cycling in the methanic zone that typically underlies the SMT and could be important in supporting the deep biosphere.

Highlights • 3D seismic imaging of an entire landslide complex. • Shallow gas accumulation within and underneath Tuaheni Landslide Complex. • Imaging of a basal shear zone within a subaqueous landslide complex. Abstract The Hikurangi margin is an active continental margin east of New Zealand's North Island. It is well recognized as a seismically active zone and is known for the occurrence of free gas and gas hydrates within the shallow sediments. A variety of subaqueous landslides can be observed at the margin, including the Tuaheni Landslide Complex off Poverty Bay. This slide complex has been interpreted previously as a slowly creeping landform, as its morphology and internal deformation is comparable to terrestrial earthflows and rock glaciers. In 2014, we acquired a high-resolution 3D seismic volume covering major parts of the Tuaheni South landslide. The 3D data show a variety of fluid migration indicators, free gas accumulations and manifestations of the base of gas hydrate stability in the pre-slide sedimentary units and the lower unit of the landslide system. The data also show that the landslide system is composed of an upper and lower unit that are separated by an intra-debris negative-polarity reflection. Free gas accumulations directly beneath the landslide units suggest that the debris acts as a boundary for rising fluids and only few migration pathways to the intra-debris reflector are observed in the distal parts of the landslide. Deformation within the landslide's debris is focused in the upper landslide unit, and we interpret the intra-debris reflector as a basal shear zone or ‘glide plane’ upon which the debris has been remobilized. The origin of the intra-debris reflector is unclear, but we suggest it could be a relatively coarse-grained horizon that would be prone to fluid flow focusing and the development of excess fluid pressure. Our seismic study provides one of the most detailed examples of a subaqueous landslide system and reveals insights into the fluid flow system and potential basal shear zone development of the Tuaheni Landslide Complex.

International Ocean Discovery Program (IODP) Expedition 372 combined two research topics, slow slip events (SSEs) on subduction faults (IODP Proposal 781A-Full) and actively deforming gas hydrate-bearing landslides (IODP Proposal 841-APL). Our study area on the Hikurangi margin, east of the coast of New Zealand, provided unique locations for addressing both research topics.SSEs at subduction zones are an enigmatic form of creeping fault behavior. They typically occur on subduction zones at depths beyond the capabilities of ocean floor drilling. However, at the northern Hikurangi subduction margin they are among the best-documented and shallowest on Earth. Here, SSEs may extend close to the trench, where clastic and pelagic sediments about 1.0-1.5 km thick overlie the subducting, seamount-studded Hikurangi Plateau. Geodetic data show that these SSEs recur about every 2 years and are associated with measurable seafloor displacement. The northern Hikurangi subduction margin thus provides an excellent setting to use IODP capabilities to discern the mechanisms behind slow slip fault behaviour.

Although the global flux of sediment and carbon from land to the coastal ocean is well known, the volume of material that reaches the deep ocean—the ultimate sink—and the mechanisms by which it is transferred are poorly documented. Using a globally unique data set of repeat seafloor measurements and samples, we show that the moment magnitude (Mw) 7.8 November 2016 Kaikōura earthquake (New Zealand) triggered widespread landslides in a submarine canyon, causing a powerful “canyon flushing” event and turbidity current that traveled >680 km along one of the world’s longest deep-sea channels. These observations provide the first quantification of seafloor landscape change and large-scale sediment transport associated with an earthquake-triggered full canyon flushing event. The calculated interevent time of ~140 years indicates a canyon incision rate of 40 mm year−1, substantially higher than that of most terrestrial rivers, while synchronously transferring large volumes of sediment [850 metric megatons (Mt)] and organic carbon (7 Mt) to the deep ocean. These observations demonstrate that earthquake-triggered canyon flushing is a primary driver of submarine canyon development and material transfer from active continental margins to the deep ocean., peer-reviewed

Carbonate escarpments are submarine limestone and dolomite cliffs that have been documented in numerous sites around the world. Their geomorphic evolution is poorly understood due to difficulties in assessing escarpment outcrops and the limited resolution achieved by geophysical techniques across their steep topographies. The geomorphic evolution of carbonate escarpments in the Mediterranean Sea has been influenced by the Messinian salinity crisis (MSC). During the MSC (5.97–5.33 Ma), the Mediterranean Sea became a saline basin due to a temporary restriction of the Atlantic-Mediterranean seaway, resulting in the deposition of more than one million cubic kilometres of salt. The extent and relative chronology of the evaporative drawdown phases associated to the MSC remain poorly constrained. In this paper we combine geophysical and sedimentological data from the central Mediterranean Sea to reconstruct the geomorphic evolution of the Malta Escarpment and infer the extent and timing of evaporative drawdown in the eastern Mediterranean Sea during the MSC. We propose that, during a MSC base-level fall, fluvial erosion formed a dense network of canyons across the Malta Escarpment whilst coastal erosion developed extensive palaeoshorelines and shore platforms. The drivers of geomorphic evolution of the Malta Escarpment after the MSC include: (i) canyon erosion by submarine gravity flows, with the most recent activity taking place, This research was undertaken with funding from Marie Curie Career Integration Grant PCIG13-GA-2013-618149 (SCARP), ERC Starting Grants n°258482 (CODEMAP) and n°677898 (MARCAN), and collaborative project n°228344 (EUROFLEETS), all within the 7th European Community Framework Programme and the Horizon 2020 Programme. Financial support was also provided by the Fulbright Visiting Scholar Program, Royal Society of New Zealand (through the International Mobility Fund), New Zealand Crown Research Institute SSIF funding to NIWA, the Griffith Geoscience Awards (Department of Communications, Energy and Natural Resources under the National Geoscience Programme 2007–2013 of Ireland), and the David and Lucile Packard Foundation. MAG acknowledges INSU for cruise-related funding and the European Union FP7 project ASTARTE for post-cruise financial support. The article is based upon work from COST Action CA15103 “Uncovering the Mediterranean salt giant” (MEDSALT) supported by COST (European Cooperation in Science and Technology).

Amblas, David ... et al.-- 22 pages, 8 figures, 1 table, Submarine canyons are deep incisions observed along most of the world’s continental margins. Their topographic relief is as dramatic as that of any canyon or river valley on land but is hidden beneath the surface of the ocean. Our knowledge of canyons has therefore come primarily from remote sensing and sampling, and has involved contributions from various oceanographic disciplines. Canyons are a critical link between coastal and shelf waters and abyssal depths; water masses, sediment, nutrients, and even litter and pollutants are carried through them. Advances in technology continue to provide new insights into canyon environments by pushing the frontier of deep marine observations and measurements. In this chapter we describe the main geomorphic features of submarine canyons and what is known about their formation and the processes that shape them. We also consider submarine gullies, which are small valleys commonly found within or alongside submarine canyons on the continental slope and may represent an incipient stage of canyon development, This work was supported by the European Union’s Horizon 2020 research and innovation programme under the Marie Sklodowska-Curie grant agreement No 658358 (D. Amblas)

The Cook Strait Canyon system, located between the North and South islands of New Zealand, is a large (1800 km 2 ), multi-branching, shelf-indenting canyon on an active subduction margin. The canyon comes within 1 km of the coast, but does not intercept fluvial or littoral sediment systems and is therefore defined as a non-terrestrially coupled system. Sediment transport associated with a strong tidal stream, and seafloor disturbance related to numerous high-activity faults, is known from previous studies. Little is known, however, about the rates of sedimentary activity in the canyon and the processes driving it. A substantial dataset of EM300 multibeam bathymetry, gravity cores, 3.5 kHz seismic reflection profiles, camera and video transects and current meter data have been collected across the region between 2002 and 2011. The canyon system therefore provides an excellent study area for understanding sediment transport in a non-coupled submarine canyon system. Analysis of the data reveals a two-staged sediment transport system where: (1) oceanographic (tidal) processes mobilise sediment from the continental shelf and transport it to depocentres in the upper–central canyons, and (2) tectonic (earthquake) processes remobilise sediment that is transported through the lower canyon to the deep ocean. Tidal boundary-layer currents within the canyon reach velocities up to 0.53 m/s and are capable of mobilising fine sand in the central reach of the upper canyons. The velocity is higher at the canyon rim and capable of mobilising coarse sand. Sediment depocentres resulting from this tidally forced sediment transport have a well formed geomorphology within the mid-upper canyon arms of Cook Strait and Nicholson Canyons. Pseudo-static stability modelling, supported by sediment core analysis, indicates that sediment accumulated in the upper canyons fails during seismic events approximately every 100 years. The 100 year return period ground shaking-level (peak ground acceleration, ignoring the effect of the water column above the seabed) at this site is estimated to be 0.23 g . Fresh rock outcrops and bed-scour in the lower canyon floor indicate that remobilised material is transported to the deep ocean. The processes identified here are likely to be analogous to those occurring in many non-coupled shelf-indenting canyons on active margins globally, and provide a framework within which the biological response to geomorphic processes in submarine canyons can be assessed.