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.
Read full abstract