Description: 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.

Description: Although offshore freshened groundwater (OFG) systems have been documented in numerous continental margins worldwide, their geometry, controls and emplacement dynamics remain poorly constrained. Here we integrate controlled-source electromagnetic, seismic reflection and borehole data with hydrological modelling to quantitatively characterise a previously unknown OFG system near Canterbury, New Zealand. The OFG system consists of one main, and two smaller, low salinity groundwater bodies. The main body extends up to 60 km from the coast and a seawater depth of 110 m. We attribute along-shelf variability in salinity to permeability heterogeneity due to permeable conduits and normal faults, and to recharge from rivers during sea level lowstands. A meteoric origin of the OFG and active groundwater migration from onshore are inferred. However, modelling results suggest that the majority of the OFG was emplaced via topographically-driven flow during sea level lowstands in the last 300 ka. Global volumetric estimates of OFG will be significantly revised if active margins, with steep coastal topographies like the Canterbury margin, are considered.

Description: 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.

Description: Slow slip events (SSEs) accommodate a significant proportion of tectonic plate motion at subduction zones, yet little is known about the faults that actually host them. The shallow depth (〈2 km) of well-documented SSEs at the Hikurangi subduction zone offshore New Zealand offers a unique opportunity to link geophysical imaging of the subduction zone with direct access to incoming material that represents the megathrust fault rocks hosting slow slip. Two recent International Ocean Discovery Program Expeditions sampled this incoming material before it is entrained immediately down-dip along the shallow plate interface. Drilling results, tied to regional seismic reflection images, reveal heterogeneous lithologies with highly variable physical properties entering the SSE source region. These observations suggest that SSEs and associated slow earthquake phenomena are promoted by lithological, mechanical, and frictional heterogeneity within the fault zone, enhanced by geometric complexity associated with subduction of rough crust.

Description: 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.

Description: 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.

Description: Sub-seabed fluid flow, gas hydrate accumulation and seafloor methane seepage are tightly interwoven processes with implications for marine biodiversity, ocean chemistry and seafloor stability. We combine long-offset seismic reflection data with high-resolution seismic data to investigate shallow structural deformation and its relationship to focused gas migration and hydrate accumulation in the southern Hikurangi subduction wedge. Anticlines, effective traps for focusing free gas, are characterized by both normal faults and vertical zones of hydraulic fracturing within the hydrate stability zone. The normal faults form as a result of sediment layer folding and gravitational collapse of ridges during uplift. We document both longitudinal (ridge-parallel) and transverse (ridge-perpendicular) extensional structures (normal faults and elongated hydraulic fracture zones) in the sub-seafloor of anticlinal ridges. Intriguingly, gas flow through ridges close to the deformation front of the wedge exploits longitudinal structures, while ridges further inboard are characterized by gas flow along transverse structures. This highlights pronounced changes in the shallow deformation of ridges in different parts of the wedge, associated with a switching of the least and intermediate principal stress directions. It is critical to understand these shallow stress fields because they control fluid flow patterns and methane seepage out of the seafloor. Key Points Gas migration through ridges occurs along both longitudinal (ridge-parallel) and transverse (ridge-perpendicular) zones of fracturing Shallow stress fields differ significantly between ridges, reflecting differences in ridge evolution and deformation Seismic reflection images of the base of gas hydrate stability and gas-water contacts are strongly affected by seismic frequency content

Description: Although marine controlled source electromagnetic (CSEM) methods are effective for investigating offshore freshened groundwater (OFG) systems, interpreting the spatial extent and salinity of OFG remains challenging. Integrating CSEM resistivity models with information on sub-surface properties, such as host-rock porosity, allows for estimates of pore-water salinity. However, deterministic inversion approaches pose challenges in quantitatively analyzing these estimates as they provide only one best-fit model with no associated estimate of model parameter uncertainty. To address this limitation, we employ a trans-dimensional Markov-Chain Monte-Carlo inversion on marine CSEM data, under the assumption of horizontal stratification, collected from the Canterbury Bight, New Zealand. We integrate the resulting posterior distributions of electrical resistivity with borehole and seismic reflection data to quantify pore-water salinity with uncertainty estimates. The results reveal a low-salinity groundwater body in the center of the survey area at varying depths, hosted by consecutive silty- and fine-sand layers approximately 20 to 60 km from the coast. These observations support the previous study’s results obtained through deterministic 2-D inversion and suggest freshening of the OFG body closer to the shore within a permeable, coarse-sand layer 40 to 150 m beneath the seafloor. This implies a potential active connection between the OFG body and the terrestrial groundwater system. We demonstrate how the Bayesian approach constrains the uncertainties in resistivity models and subsequently in pore-water salinity estimates. Our findings highlight the potential of Bayesian inversions in enhancing our understanding of OFG systems, providing crucial boundary conditions for hydrogeological modeling and sustainable water resource development.

Description: The study of offshore freshened groundwater (OFG) is gaining importance due to population growth and environmental pressure on coastal water resources. Marine controlled source electromagnetic (CSEM) methods can effectively map the spatial extent of OFG systems using electrical resistivity as a proxy. Integrating these resistivity models with sub-surface properties, such as host-rock porosity, allows for estimates of pore-water salinity. However, evaluating the uncertainty in pore-water salinity using resistivity models obtained from deterministic inversion approaches presents challenges, as they provide only one best-fit model, with no associated estimate of uncertainty. To address this limitation, we employ trans-dimensional Markov-Chain Monte-Carlo inversion on marine time-domain CSEM data, acquired in the Canterbury Bight, New Zealand. We integrate resistivity posterior probability distributions with borehole and seismic reflection data to estimate pore-water salinity with corresponding uncertainty estimates. The results highlight a low-salinity groundwater body in the center of the survey area, hosted by consecutive silty- and fine-sand layers approximately 20–60 km from the coast. The posterior probability distribution of resistivity models indicates freshening of the OFG body toward the shoreline within a permeable, coarse-sand layer 40–150 m beneath the seafloor, suggesting an active connection between the OFG body and the terrestrial groundwater system. The approach demonstrates how Bayesian inversion constrains the uncertainties in resistivity models and subsequently in pore-water salinity estimates. Our findings highlight the potential of Bayesian inversion to enhance our understanding of OFG systems and provide uncertainty constraints for hydrogeological modeling, thereby contributing to sustainable water resource development. Key Points A Bayesian workflow is employed to evaluate uncertainty in pore-water salinity estimates Offshore groundwater in Canterbury Bight stores freshened pore-water in fine-grained sediments, likely extending from the onshore aquifer Correlation between pore-water salinities and seismic-derived stratigraphy provides boundary conditions for hydrogeological modeling Plain Language Summary Geophysical methods that measure the electromagnetic properties of the Earth are effective in investigating freshwater sources beneath the seafloor. By combining the geophysical and geological information, we can better assess the quality of this groundwater. In this study, we develop a workflow that uses statistical methods to integrate electromagnetic observations with borehole and acoustic measurements along the eastern coast of the South Island of New Zealand. We aim to improve our understanding of the groundwater quality beneath the seafloor. Our research confirms the presence of freshened groundwater within the sandy seafloor up to 60 km from the coastline. Importantly, our observations indicate that the groundwater quality increases toward the coast. These findings are significant as they enhance the hydrogeological modeling of the groundwater system and suggest its potential as a source of freshwater.

Description: 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.