Search

Highlights • Multilevel Composition is an innovative method involving color composition and co-rendering of multilevel attribute maps. • It is useful for characterizing multi-depth geological features based on their spatiotemporal distribution within three-dimensional seismic data. • The technique produces a single image map, in which inter-window/layer depth information is coded in colors for reliable representation of the actual geology. • In the eastern Nile fan, it was applied to visualize and resolve the complexities of buried clastic deep-water depositional elements. • On the Omakere Ridge, it successfully illuminated seafloor seeps and reveals their link to deeper fluid-bearing intervals. Abstract Advanced seismic data and multi-attribute visualization techniques, such as color blending of attributes, have considerably enhanced the capability of interpreters to characterize geological features in three-dimensional (3D) seismic reflection datasets. However, high resolution investigation of complex, vertically linked geological features such as channel systems and fluid conduits, remains challenging. These features may appear in the dataset as pronounced attribute anomalies, such as high-amplitude or spectrally or structurally enhanced seismic reflectivity bands, at several depth levels. Vertical linkages between these features, however, may not be readily established. We have developed an innovative method, Multilevel Composition, for an intuitive display of vertically connected features. Our method involves the composition of attribute maps from three different depth/time windows or slices onto a single map, in which inter-window/layer depth information is coded in colors. Multilevel Composition starts with the identification of suitable seismic attributes, such as high amplitudes in the examples displayed here, to map features of geological interest. At least one reference horizon is then identified and mapped in the vicinity of the target window of interest. Three sub-windows are then defined with respect to the reference horizon(s) based on the vertical and spatial distribution of the geological features. Relevant seismic attributes are computed for each of the sub-windows, and the resulting maps, one from each sub-window, are assigned basic color channels and are co-rendered to reveal multilevel linkages between these features. We demonstrate the efficacy of this method by applying it to two 3D seismic datasets, one illuminating deep-water depositional elements in the eastern Nile fan, eastern Mediterranean and the other targeting seafloor seeps and underlying gas migration systems beneath the Omakere Ridge, offshore New Zealand. The new method is simple and should be easy to implement to enhance seismic interpretation workflows.

Gas hydrates are ice-like crystalline materials that form under submarine environments of moderate pressure and low temperature. Another key factor to their formation is the abundance in gas supply from depth in addition to local biogenic gas. Detailed imaging and velocity analysis of the plumbing system of gas hydrates can provide confidence that amplitude anomalies in seismic data are related to gas hydrate accumulations. We have conducted 2D elastic full-waveform inversion (FWI) along a 14 km long segment of a 2D multichannel seismic profile to obtain a high-resolution velocity model of a hydrate system on the southern Hikurangi margin. We compare the FWI velocity model to previously published semblance- and tomography-based velocity models from the same data to explore how much more can be gained from the FWI. The FWI yielded a structurally more accurate velocity model that better delineated the low-velocity zone associated with free gas beneath the bottom simulating reflector (BSR) compared to the semblance- and tomography-based velocity models. Our results also find a lateral velocity inversion, that is, a narrow low-velocity zone surrounded by bands of higher velocities at a seaward-verging protothrust fault, which the two other methodologies failed to resolve. The FWI provides an improved lateral resolution making it an important tool when imaging the “plumbing” systems of gas hydrate reservoirs. In the southeastern limb of the anticline, our results find that the closely spaced landward-vergent protothrusts provide gas-charged fluids for hydrate formation above the BSR. Moreover, at the center of the anticline, our results find that a seaward-vergent protothrust fault appears to be acting as a conduit for gas-rich fluids into strata, although there is no accumulation of any significant hydrate above the BSR at the apex of the anticline. Our finding emphasizes the significance of densely spaced faults and fractures for providing gas for hydrate formation in the hydrate stability zone.

Glacio-eustatic cycles lead to changes in sedimentation on all types of continental margins. There is, however, a paucity of sedimentation rate data over eustatic sea-level cycles in active subduction zones. During International Ocean Discovery Program Expedition 375, coring of the upper ∼110 m of the northern Hikurangi Trough Site U1520 recovered a turbidite-dominated succession deposited during the last ∼45 kyrs (Marine Isotope Stages (MIS) 1–3). We present an age model integrating radiocarbon dates, tephrochronology, and δ18O stratigraphy, to evaluate the bed recurrence interval (RI) and sediment accumulation rate (SAR). Our analyses indicate mean bed RI varies from ∼322 yrs in MIS1, ∼49 yrs in MIS2, and ∼231 yrs in MIS3. Large (6-fold) and abrupt variations in SAR are recorded across MIS transitions, with rates of up to ∼10 m/kyr occurring during the Last Glacial Maximum (LGM), and

We constrain orientations of the horizontal stress field from borehole image data in a transect across the Hikurangi Subduction Margin. This region experiences NW‐SE convergence and is the site of recurrent slow slip events. The direction of the horizontal maximum stress is E‐W at an active splay thrust fault near the subduction margin trench. This trend changes to NNW‐SSE in a forearc trench slope basin on the offshore accretionary wedge, and to NE‐SW in the onshore forearc. Multiple, tectonic, and geological processes, either individually or in concert, may explain this variability. The observed offshore to onshore stress rotation may reflect a change from dominantly compressional tectonics at the deformation front, to a strike‐slip and/or extensional tectonic regime closer to the Taupo Volcanic Zone, further inland. In addition, the offshore stress may be affected by topography and/or stress rotation around subducting seamounts, and/or temporal stress changes during the slow slip cycle.

The Pāpaku Fault Zone, drilled at International Ocean Discovery Program (IODP) Site U1518, is an active splay fault in the frontal accretionary wedge of the Hikurangi Margin. In logging‐while‐drilling data, the 33‐m‐thick fault zone exhibits mixed modes of deformation associated with a trend of downward decreasing density, P‐wave velocity, and resistivity. Methane hydrate is observed from ~30 to 585 m below seafloor (mbsf), including within and surrounding the fault zone. Hydrate accumulations are vertically discontinuous and occur throughout the entire logged section at low to moderate saturation in silty and sandy centimeter‐thick layers. We argue that the hydrate distribution implies that the methane is not sourced from fluid flow along the fault but instead by local diffusion. This, combined with geophysical observations and geochemical measurements from Site U1518, suggests that the fault is not a focused migration pathway for deeply sourced fluids and that the near‐seafloor Pāpaku Fault Zone has little to no active fluid flow.

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.

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.

Submarine currents are a principal factor in controlling seafloor geomorphology. Herein, we investigate the role of dynamic current systems associated with the Subtropical Front in the formation an...

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.

In this study we provide evidence for methane hydrates in the Taranaki Basin, occurring a considerable distance from New Zealand's convergent margins, where they are well documented. We describe and reconstruct a unique example of gas migration and leakage at the edge of the continental shelf, linking shallow gas hydrate occurrence to a deeper petroleum system. The Taranaki Basin is a well investigated petroleum province with numerous fields producing oil and gas. Industry standard seismic reflection data show amplitude anomalies that are here interpreted as discontinuous BSRs, locally mimicking the channelized sea-floor and pinching out up-slope. Strong reverse polarity anomalies indicate the presence of gas pockets and gas-charged sediments. PetroMod™ petroleum systems modelling predicts that the gas is sourced from elevated microbial gas generation in the thick slope sediment succession with additional migration of thermogenic gas from buried Cretaceous petroleum source rocks. Cretaceous–Paleogene extensional faults underneath the present-day slope are interpreted to provide pathways for focussed gas migration and leakage, which may explain two dry petroleum wells drilled at the Taranaki shelf margin. PetroMod™ modelling predicts concentrated gas hydrate formation on the Taranaki continental slope consistent with the anomalies observed in the seismic data. We propose that a semi-continuous hydrate layer is present in the down-dip wall of incised canyons. Canyon incision is interpreted to cause the base of gas hydrate stability to bulge downward and thereby trap gas migrating up-slope in permeable beds due to the permeability decrease caused by hydrate formation in the pore space. Elsewhere, hydrate occurrence is likely patchy and may be controlled by focussed leakage of thermogenic gas. The proposed presence of hydrates in slope sediments in Taranaki Basin likely affects the stability of the Taranaki shelf margin. While hydrate presence can be a drilling hazard for oil and gas exploration, the proposed presence of gas hydrates opens up a new frontier for exploration of hydrates as an energy source.

Using automated supervised segmentation of multibeam backscatter data to delineate seafloor substrates is a relatively novel technique. Low-frequency multibeam echosounders (MBES), such as the 12-kHz EM120, present particular difficulties since the signal can penetrate several metres into the seafloor, depending on substrate type. We present a case study illustrating how a non-targeted dataset may be used to derive information from multibeam backscatter data regarding distribution of substrate types. The results allow us to assess limitations associated with low frequency MBES where sub-bottom layering is present, and test the accuracy of automated supervised segmentation performed using SonarScope® software. This is done through comparison of predicted and observed substrate from backscatter facies-derived classes and substrate data, reinforced using quantitative statistical analysis based on a confusion matrix. We use sediment samples, video transects and sub-bottom profiles acquired on the Chatham Rise, east of New Zealand. Inferences on the substrate types are made using the Generic Seafloor Acoustic Backscatter (GSAB) model, and the extents of the backscatter classes are delineated by automated supervised segmentation. Correlating substrate data to backscatter classes revealed that backscatter amplitude may correspond to lithologies up to 4 m below the seafloor. Our results emphasise several issues related to substrate characterisation using backscatter classification, primarily because the GSAB model does not only relate to grain size and roughness properties of substrate, but also accounts for other parameters that influence backscatter. Better understanding these limitations allows us to derive first-order interpretations of sediment properties from automated supervised segmentation.

The bottom simulating reflection (BSR) is widely observed along continental margins and is believed to mark the base of gas hydrate stability zone (BGHSZ). In some regions, double or multiple overlapping BSRs are observed, yet their formation mechanisms and geologic implications are not well understood. Here we present 3D seismic images from the 2018 NZ3D experiment that covers a 14 × 60 km 2 survey area on New Zealand's northern Hikurangi subduction margin. We observe double BSRs in five locations. Beneath the Tuaheni Basin in the mid-slope, a secondary BSR (BSR2) lies ∼100-360 m deeper than the primary BSR (BSR1) and its 3D geometry mimics the unconformity at the base of the basin. At three thrust ridges located 18-38 km from the deformation front, BSR2 lies ∼55-130 m below and is sub-parallel to BSR1. At another thrust ridge ∼14 km from the deformation front, BSR2 forms above the BSR1, and the two BSRs converge towards the peak of the ridge. Through 3D modeling of BGHSZ and analysis of the geometry and reflection characteristics of the double BSRs, we identify three potential mechanisms for their formation (1) rapid sedimentation, (2) tectonic uplift, (3) overpressure/heat advection caused by fluid migration. Our study demonstrates that formation of double BSRs is closely linked to subduction processes along the northern Hikurangi margin, and double BSRs may be used as indicators for areas with recent sedimentation, tectonic and/or fluid activities.

GNS Science's MBIE funded contracts. Grant Number: CO5X1204 and C05X0908, “Gas Hydrates Resources,” Marsden Fund project. Grant Number: GNS0902

The Hikurangi Margin, east of the North Island of New Zealand, is known to contain significant deposits of gas hydrates. This has been demonstrated by several multidisciplinary studies in the area since 2005. These studies indicate that hydrates in the region are primarily located beneath thrust ridges that enable focused fluid flow, and that the hydrates are associated with free gas. In 2009–2010, a seismic dataset consisting of 2766 km of 2D seismic data was collected in the undrilled Pegasus Basin, which has been accumulating sediments since the early Cretaceous. Bottom-simulating reflections (BSRs) are abundant in the data, and they are accompanied by other features that indicate the presence of free gas and concentrated accumulations of gas hydrate. We present results from a detailed qualitative analysis of the data that has made use of automated high-density velocity analysis to highlight features related to the hydrate system in the Pegasus Basin. Two scenarios are presented that constitute contrasting mechanisms for gas-charged fluids to breach the base of the gas hydrate stability zone. The first mechanism is the vertical migration of fluids across layers, where flow pathways do not appear to be influenced by stratigraphic layers or geological structures. The second mechanism is non-vertical fluid migration that follows specific strata that crosscut the BSR. One of the most intriguing features observed is a presumed gas chimney within the regional gas hydrate stability zone that is surrounded by a triangular (in 2D) region of low reflectivity, approximately 8 km wide, interpreted to be the result of acoustic blanking. This chimney structure is cored by a ∼200-m-wide low-velocity zone (interpreted to contain free gas) flanked by high-velocity bands that are 200–400 m wide (interpreted to contain concentrated hydrate deposits).

Small amounts of free gas in interstitial sediment pores are known to significantly lower compressional (P-) wave velocity (Vp). This effect, combined with moderately elevated Vp from the presence of gas hydrates, is usually thought to be the cause for the often observed strong negative reflection coefficients of bottom simulating reflections (BSRs) at the base of gas hydrate stability (BGHS). At several locations however, weak BSRs have been observed, which are difficult to reconcile with a presence of gas in sediment pores. We here present a rock physics model for weak BSRs on the Hikurangi Margin east of New Zealand. Thin sections of a fine-grained mudstone sample from a submarine outcrop in the vicinity of a weak BSR show macroscopic porosity in the form of fractures and intrafossil macropores. We apply the Kuster-Toksoz theory to predict seismic velocities for a rock with water-saturated interstitial micropores and gas or hydrates in macroscopic pore space simulating fractures or compliant macropores. We match field observations of a weak BSR with a reflection coefficient of −0.016 with two end-member models; (1) rocks with gas hydrate-filled voids with a concentration of

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.

Subduction erosion is a long-term, large-scale geological process that dominates the structural inventory of a large number of convergent margins. Along the Peruvian convergent margin, this process has controlled the structural deformation of the forearc as a consequence of spatio-temporal variations during the Cenozoic. However, the pace at which tectonic erosion occurred has not been linked to regional, subduction-related events. The aim of this study is to show that subduction erosion along the Central Peruvian convergent margin off Trujillo (7∘S–9∘S) was influenced by plate convergence variability since the Late Cretaceous. Drilling data, multi-channel seismic (MCS) data and subsidence analyses are used to investigate the long-term effect of subduction erosion on the tectonic evolution of the Peru margin. Our study shows that subduction erosion (1) did not occur between the Cretaceous and Early Eocene (110–56 Ma), (2) increased and fluctuated between the Eocene and Pliocene (56–5 Ma) with a maximum average trench retreat of 3 km.Myr−1 between the Middle Miocene and Late Miocene (20–5 Ma), and (3) has decreased since the Pliocene (5 Ma) (1 km.Myr−1 of average trench retreat). A remarkable finding is that shorter periods of faster subduction erosion coincide with or follow major plate reorganization events that were marked by changes in the velocity and direction of relative plate convergence as well as with well-known phases of Andean orogeny. Based on our analyses, we therefore conclude that convergence variability off Peru is a primary mechanism controlling the pace at which tectonic erosion occurred along the Trujillo Basin during the Cenozoic.

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.

We present an integrated 2D model of thermal and microbial generation of methane, migration into the gas hydrate stability zone (HSZ), and formation of methane hydrates. The model reconstructs the shallow (0–20 km) thermal structure of the subduction interface between the Australian plate and the subducting Pacific plate, and the trench basin (Pegasus Basin). Modelled temperatures of less than 110 °C within Pegasus Basin constrain the generation of oil and gas. Whilst a cool thermal regime is predicted to limit thermogenic generation of gas to a burial depth of >10 km, it extends the interval where prolific microbial gas generation occurs. The modelled rate of microbial generation of methane increases beneath the HSZ and peaks at ~1600 m below seafloor. Diffusive upward migration of microbially generated methane is interpreted to lead to widespread methane hydrate formation and the presence of a semi-continuous bottom simulating reflector (BSR). Predicted average hydrate saturation within the HSZ is 0.9% for a modelled sedimentary organic matter content of 0.5% and 1.6% for 1% organic matter in fine-grained Pegasus Basin sediments. Considerably higher concentrations of methane hydrate of up to 20–70% are predicted to occur where gas migration is focussed within the frontal anticline and proto-thrust zone southeast of the modern accretionary wedge and in channel and basin floor sandstones related to the Hikurangi Channel. The Hikurangi Channel sedimentary system transported coarse clastic sediments eroded from the rising Southern Alps along the eastern margin of the Pegasus Basin since the Miocene. It provides carrier beds specifically for transport of thermogenic gas generated close to the subduction interface. A buried Mesozoic accretionary wedge originating from subduction of the Pacific Plate beneath Gondwana further focusses the migration of gas. Focussed migration of thermogenic gas leads to the highest predicted hydrate concentrations in potential channel sand reservoirs.

We present multichannel seismic data from New Zealand's Hikurangi subduction margin that show widespread evidence for gas migration into the field of gas hydrate stability. Gas migration along stratigraphic layers into the hydrate system manifests itself as highly reflective segments of dipping strata that originate at the base of hydrate stability and extend some distance toward the seafloor. The highly reflective segments exhibit the same polarity as the seafloor reflection, indicating that localized gas hydrate precipitation from gas-charged fluids within relatively permeable layers has occurred. High-density velocity analysis shows that these layer-constrained gas hydrate accumulations are underlain by thick (up to ~500 m) free gas zones, which provide the source for focused gas migration into the hydrate layer. In addition to gas being channeled along layers, we also interpret gas migration through a fault zone into the field of hydrate stability; in this case, a low-velocity layer within the hydrate stability zone extends laterally away from the fault, which might indicate that gas-charged fluids have also migrated away from the fault along strata. At this site, where both dipping strata and faulting seem to influence fluid migration, we observe anomalously high velocities at the base of hydrate stability that we interpret as concentrated gas hydrates. Our results give insight into how shallow fluid flow responds to permeability contrasts between strata, fault zones, and perhaps also the gas hydrate system itself. Ultimately, these relationships can lead to gas migration across the base of hydrate stability and the precipitation of concentrated hydrate deposits.

Gas seepage from marine sediments has implications for understanding feedbacks between the global carbon reservoir, seabed ecology and climate change. Although the relationship between hydrates, gas chimneys and seafloor seepage is well established, the nature of fluid sources and plumbing mechanisms controlling fluid escape into the hydrate zone and up to the seafloor remain one of the least understood components of fluid migration systems. In this study we present the analysis of new three-dimensional high-resolution seismic data acquired to investigate fluid migration systems sustaining active seafloor seepage at Omakere Ridge, on the Hikurangi subduction margin, New Zealand. The analysis reveals at high resolution, complex overprinting fault structures (i.e. protothrusts, normal faults from flexural extension, and shallow (

A TOUGH + HYDRATE flow model analogous of hydrate deposits in the Hikurangi margin has been constructed using the best available estimates of permeability, porosity and hydrate saturation. Seismic velocity and AVO analyses have been used to develop this model. Production rate forecasts were made in a single well model with homogeneous properties. In this forecast, the system was depressurised by producing gas from an underlying free gas zone included in the model. Production of 4.6 Bscf (130 × 106 m3) of gas was achieved from the well over 20 years with average flow rates of 0.65 MMscf/day (0.213 m3/s).

Gas hydrates have recently been recognised as a class of unconventional petroleum resource and the economic viability of gas production from hydrates is now being viewed as a realistic possibility within the next decade. Therefore, potential offshore hydrate accumulations in the world-class endowed gas hydrate province, the Hikurangi Margin, New Zealand, represent a significant medium- to long-term opportunity to meet the country's future energy requirements. In this paper we delineate a potential gas hydrate reservoir in the East Coast Basin, New Zealand and quantitatively estimate its gas hydrate concentrations from 2D seismic data with no well information available. The target is interesting for exploration since it shows evidence for gas-hydrate bearing sands, in particular, buried channel systems. We use a combined analysis of high-resolution velocity analysis, amplitude-versus-offset (AVO) attribute and AVO inversion to investigate whether we can identify regions that are likely to contain highly concentrated gas hydrates and whether they are likely to be sand-dominated. To estimate hydrate concentrations we apply a rock physics model. Our results indicate the presence of several – up to 200 m thick – zones that are likely to host gas hydrates, with one location predicted to consist of high-permeable channel sands and an inferred gas hydrate saturation of ∼25%. These findings suggest significant amounts of gas hydrates may be present in high-quality reservoirs on this part of the margin.

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.

Supplies of conventional natural gas and oil are declining fast worldwide, and therefore new, unconventional forms of energy resources are needed to meet the ever-increasing demand. Amongst the many different unconventional natural resources are gas hydrates, a solid, ice-like crystalline compound of methane and water formed under specific low temperature and high pressure conditions. Gas hydrates are believed to exist in large quantities worldwide in oceanic regions of continental margins, as well as associated with permafrost regions in the Arctic. Some studies to estimate the global abundance of gas hydrate suggest that the total volume of natural gas locked up in form of gas hydrates may exceed all known conventional natural gas reserves, although large uncertainties exist in these assessments. Gas hydrates have been intensively studied in the last two decades also due to connections between climate forcing (natural and/or anthropogenic) and the potential large volumes of methane trapped in gas hydrate accumulations. The presence of gas hydrate within unconsolidated sediments of the upper few hundred meters below seafloor may also pose a geo-hazard to conventional oil and gas production. Additionally, climate variability and associated changes in pressure-temperature regimes and thus shifts in the gas hydrate stability zone may cause the occurrence of submarine slope failures. Several large-scale national gas hydrate programs exist especially in countries such as Japan, Korea, Taiwan, China, India, and New Zealand, where large demands of energy cannot be met by domestic supplies from natural resources. The past five years have seen several dedicated deep drilling expeditions and other scientific studies conducted throughout Asia and Oceania to understand gas hydrates off India, China, and Korea. This thematic set of publications is dedicated to summarize the most recent findings and results of geo-scientific studies of gas hydrates in the marginal seas and continental margin of the Asia, and Oceania region.

We present indirect evidence from seismic data for gas hydrates in the deepwater Northland and northern Taranaki Basins west of New Zealand's North Island. We observe high-amplitude reflections in the regional gas-hydrate stability zone that we interpret as gas pockets and infer that the gas is likely to be surrounded by gas hydrate. We also observe a reflection that mimics the shape of the base of gas-hydrate stability and is possibly a bottom simulating reflection. These preliminary findings, which need to be corroborated by further analyses, may constitute the first evidence for gas hydrates in these basins close to New Zealand's primary petroleum province.

Dissolved methane and high resolution bathymetry surveys were conducted over the Rock Garden region of Ritchie Ridge, along the Hikurangi margin, eastern New Zealand. Multibeam bathymetry reveals two prominent, northeast trending ridges, parallel to subduction along the margin, that are steep sided and extensively slumped. Elevated concentrations of methane (up to 10 nM, 10× background) within the water column are associated with a slump structure at the southern end of Eastern Rock Garden. The anomalous methane concentrations were detected by a methane sensor (METS) attached to a conductivity‐temperature‐depth‐optical backscatter device (CTDO) and are associated with elevated light scattering and flare‐shaped backscatter signals revealed by the ship's echo sounder. Increased particulate matter in the water column, possibly related to the seepage and/or higher rates of erosion near slump structures, is considered to be the cause of the increased light scattering, rather than bubbles in the water ...

SUMMARY Gas hydrates were discovered in a mud mound in the lease block Mississippi Canyon 798, Gulf of Mexico, through piston coring. Subsequently, a seismic experiment was carried out to investigate the dynamics behind the hydrate formation. During the experiment, high-resolution multichannel seismic reflection data using a 24-channel, 240 m long streamer and wide-aperture data using six ocean bottom seismometers were collected along five lines. High-reflectivity zones (HRZs) are present in the reflection data along all lines. To better constrain the interpretation of the reflection data, the traveltimes from the multichannel and wide-aperture data sets were jointly inverted to estimate a 2-D P-wave layered velocity model for each line. A minimum-parameter/minimum-structure modelling approach yielded simple models and a comparison of the models at their intersection points shows they are consistent to within ±10 m s −1 in velocity and ±20 m in depth. In the final P-wave velocity models, the HRZs are associated with a lowering of velocity. In the reflection data, the top of the HRZs show a polarity reversal with respect to the seafloor. Presence of free gas in the HRZs best explains the velocity lowering and polarity reversal. It is speculated that the gas has deeper sources and migrates upwards through conduits formed by salt movement in the vicinity. The upward migrating gas accumulates in the axis of a channel complex and manifests itself as HRZs in the reflection data. The fluids circulating along the conduits push the base of the hydrate stability zone close to the seafloor. From the channel axis, the free gas migrates further upwards and close to the seafloor, and as it comes within the gas hydrate stability zone, it forms hydrates.

This paper presents results of shear strength and acoustic velocity (p-wave) measurements performed on: (1) samples containing natural gas hydrate from the Mallik 2L-38 well, Mackenzie Delta, Northwest Territories; (2) reconstituted Ottawa sand samples containing methane gas hydrate formed in the laboratory; and (3) ice-bearing sands. These measurements show that hydrate increases shear strength and p-wave velocity in natural and reconstituted samples. The proportion of this increase depends on (1) the amount and distribution of hydrate present, (2) differences in sediment properties, and (3) differences in test conditions. Stress-strain curves from the Mallik samples suggest that natural gas hydrate does not cement sediment grains. However, stress-strain curves from the Ottawa sand (containing laboratory-formed gas hydrate) do imply cementation is present. Acoustically, rock physics modeling shows that gas hydrate does not cement grains of natural Mackenzie Delta sediment. Natural gas hydrates are best modeled as part of the sediment frame. This finding is in contrast with direct observations and results of Ottawa sand containing laboratory-formed hydrate, which was found to cement grains (Waite et al. 2004). It therefore appears that the microscopic distribution of gas hydrates in sediment, and hence the effect of gas hydrate on sediment physical properties, differs between natural deposits and laboratory-formed samples. This difference may possibly be caused by the location of water molecules that are available to form hydrate. Models that use laboratory-derived properties to predict behavior of natural gas hydrate must account for these differences.

Recently acquired seismic reflection data across the southern edge of Ritchie Ridge, a prominent bathymetric high on the Hikurangi margin, display zones of high amplitudes and reflections that crosscut strata. We interpret the latter as bottom‐simulating reflections which are commonly associated with gas beneath gas hydrates. An analysis of reflection strength indicates the high‐amplitude zones are caused by free gas in the pore space of sediments, probably migrating upward along layers. One of the high‐amplitude regions is situated beneath the projected location of a known gas vent site. The seismic data appear to image the conduits that supply this vent site with gas. The seafloor in most of the study area is likely to be within the zone of gas hydrate stability, depending on bottom water temperatures and hydrate composition. Hence, gas appears to be venting through the gas hydrate stability zone, favouring locally high concentrations of gas hydrates.

Seismic detection of methane hydrate often relies on indirect or equivocal methods. New multichannel seismic reflection data from the Blake Ridge, located approximately 450 km east of Savannah, Georgia, show three direct seismic indicators of methane hydrate: (1) a paleo bottom‐simulating reflector (BSR) formed when methane gas froze into methane hydrate on the eroding eastern flank of the Blake Ridge, (2) a lens of reduced amplitudes and high P‐wave velocities found between the paleo‐BSR and BSR, and (3) bright spots within the hydrate stability zone that represent discrete layers of concentrated hydrate formed by upward migration of gas. Velocities within the lens (∼1910 m/s) are significantly higher than velocities in immediately adjacent strata (1820 and 1849 m/s). Conservative estimates show that the hydrate lens contains at least 13% bulk methane hydrate within a 2‐km3 volume, yielding 3.2 × 1010kg [1.5 TCF (4.2 × 1010 m3] of methane. Low seismic amplitudes coupled with high interval velocities within the lens offer evidence for possible methane hydrate “blanking.” Hydrate bright spots yield velocities as high as 2100 m/s, with bulk hydrate concentrations predicted as high as 42% in an approximately 15‐m thick layer. Our results show that, under certain circumstances, hydrate in marine sediments can be directly detected in seismic reflections but that quantification of hydrate concentrations requires accurate velocity information.

As offshore petroleum exploration and development move into deeper water, industry must contend increasingly with gas hydrate, a solid compound that binds water and a low-molecular-weight gas (usually methane). Gas hydrate has been long studied in industry from an engineering viewpoint, due to its tendency to clog gas pipelines. However, hydrate also occurs naturally wherever there are high pressures, low temperatures, and sufficient concentrations of gas and water. These conditions prevail in two natural environments, both of which are sites of active exploration: permafrost regions and marine sediments on continental slopes. In this article we discuss seismic detection of gas hydrate in marine sediments. Gas hydrate in deepwater sediments poses both new opportunities and new hazards. An enormous quantity of natural gas, likely far exceeding the global inventory of conventional fossil fuels, is locked up worldwide in hydrates. Ex-traction of this unconventional resource presents unique exploration, engineering, and economic challenges, and several countries, including the United States, Japan, Canada, India, and Korea, have initiated joint industry-academic-governmental programs to begin studying those challenges. Hydrates also constitute a potential drilling hazard. Because hydrates are only stable in a restricted range of pressure and temperature, any activity that sufficiently raises temperature or lowers pressure could destabilize them, releasing potentially large volumes of gas and decreasing the shear strength of the host sediments. Assessment of the opportunities and hazards associated with hydrates requires reliable methods of detecting hydrate and accurate maps of their distribution and concentration. Hydrate may occur only within the upper few hundred meters of deepwater sediment, at any depth between the seafloor and the base of the stability zone, which is controlled by local pressure and temperature. Hydrate is occasionally exposed at the seafloor, where it can be detected either visually or acoustically by strong seismic reflection amplitudes or high backscatter …

Morphological and seismic data from a submarine landslide complex east of New Zealand indicate flow-like deformation within gas hydrate-bearing sediment. This “creeping” deformation occurs immediately downslope of where the base of gas hydrate stability reaches the seafloor, suggesting involvement of gas hydrates. We present evidence that, contrary to conventional views, gas hydrates can directly destabilize the seafloor. Three mechanisms could explain how the shallow gas hydrate system could control these landslides. (1) Gas hydrate dissociation could result in excess pore pressure within the upper reaches of the landslide. (2) Overpressure below low-permeability gas hydrate-bearing sediments could cause hydrofracturing in the gas hydrate zone valving excess pore pressure into the landslide body. (3) Gas hydrate-bearing sediment could exhibit time-dependent plastic deformation enabling glacial-style deformation. We favor the final hypothesis that the landslides are actually creeping seafloor glaciers. The viability of rheologically controlled deformation of a hydrate sediment mix is supported by recent laboratory observations of time-dependent deformation behavior of gas hydrate-bearing sands. The controlling hydrate is likely to be strongly dependent on formation controls and intersediment hydrate morphology. Our results constitute a paradigm shift for evaluating the effect of gas hydrates on seafloor strength which, given the widespread occurrence of gas hydrates in the submarine environment, may require a reevaluation of slope stability following future climate-forced variation in bottom-water temperature.

The distribution of bottom simulating reflections (BSRs), which mark the base of the gas hydrate zone on marine seismic records across the Peru margin, is examined with respect to tectonic uplift, sedimentation, and carbon content of sediment. Where accretion of only deep ocean pelagic sediment containing less than 1% carbon is seismically imaged, no BSR is observed. No BSR is apparent on seismic records of carbon-rich (4–8%) sediment of rapidly subsiding basins. BSRs are abundant on records of carbon-rich sediment affected by tectonic uplift and/or rapid sedimentation. In areas of uplift, the bottom simulating reflectors contain small amounts of free gas in contact with sediment containing small amounts of gas hydrate. During uplift, free gas is released from gas hydrates by dissociation as the hydrate phase boundary shifts upward with respect to sea level. The strong effects of free gas in reducing seismic velocity and gas hydrate in increasing velocity can explain the relatively small total quantities of gas and hydrate indicated by observations of BSRs off Peru.

Summary Bottom simulating reflectors (BSRs) at the base of the hydrate stability zone are often observed in marine seismic reflection data. The compressional (P-)wave velocity structure of a BSR offshore Costa Rica was analysed by applying a full-waveform inversion technique. The resulting velocity profile indicates that the BSR is caused by a thin low-velocity layer, suggesting the presence of at least a small amount of free gas in the sediment pore space. In undisturbed sediment sections, BSRs offshore Costa Rica are observed over much of the continental margin. They frequently occur at shallow depths beneath the seafloor, usually ~ 100–400 m; at some locations BSRs appear to intersect the seafloor. This was interpreted as an indication that hydrates in the study area form partly from methane which was produced beneath the hydrate stability zone and migrated upwards. In the study area, the vertical movement of the hydrate stability zone relative to a given point in the sediment column, one of the potential factors leading to BSR formation or suppression, is controlled by both sedimentation and vertical tectonism. Both processes may hence play a role in controlling BSR distribution. BSRs are absent in areas which have been affected by slumping, except where the sedimentary section above the BSR remained intact during slumping. This indicates that slope failure can cause the destruction of BSRs.

Much of our knowledge of the worldwide distribution of submarine gas hydrates comes from seismic observations of Bottom Simulating Reflectors (BSRs). Full waveform inversion has proven to be a reliable technique for studying the fine structure of BSRs using the compressional wave velocity. We applied a non-linear full waveform inversion technique to a BSR at a location offshore Peru. We first determined the large-scale features of seismic velocity variations using a statistical inversion technique to maximise coherent energy along travel-time curves. These velocities were used for a starting velocity model for the full waveform inversion, which yielded a detailed velocity/depth model in the vicinity of the BSR. We found that the data are best fit by a model in which the BSR consists of a thin, low-velocity layer. The compressional wave velocity drops from 2.15 km/s down to an average of 1.70 km/s in an 18 m thick interval, with a minimum velocity of 1.62 km/s in a 6 m interval. The resulting compressional wave velocity was used to estimate gas content in the sediments. Our results suggest that the low velocity layer is a 6–18 m thick zone containing a few percent of free gas in the pore space. The presence of the BSR coincides with a region of vertical uplift. Therefore, we suggest that gas at this BSR is formed by a dissociation of hydrates at the base of the hydrate stability zone due to uplift and subsequently a decrease in pressure.

The characterization of marine gas hydrates from geophysical data has evolved significantly over the past decade. For many years, marine gas hydrate interpretation was largely limited to a determination as either “likely present” or “likely not present” based on the presence or absence of

Identification of methane sources controlling hydrate distribution and concentrations in continental margins remains a major challenge in gas hydrate research. Lack of deep fluid samples and high quality regional scale seismic reflection data may lead to underestimation of the significance of fluid escape from subducting and compacting sediments in the global inventory of methane reaching the hydrate zone, the water column and the atmosphere. The distribution of concentrated hydrate zones in relation to focused fluid flow across the southern Hikurangi subduction margin was investigated using high quality, long offset (10 km streamer), pre-stack depth migrated multichannel seismic data. Analysis of low P wave velocity zones, bright-reverse polarity reflections and dim-amplitude anomalies reveals pathways for gas escape and zones of gas accumulation. The study shows the structural and stratigraphic settings of three main areas of concentrated hydrates: (1) the Opouawe Bank, dominated by focused periodic fluid input along thrust faults sustaining dynamic hydrate concentrations and gas chimneys development; (2) the frontal anticline, with a basal set of protothrusts controlling permeability for fluids from deeply buried and subducted sediments sustaining hydrate concentrations at the crest of the anticline; and (3) the Hikurangi Channel, with buried sand dominated channels hosting significant amounts of gas beneath the base of the hydrate zone. In sand dominated channels gas injection into the hydrate zone favors highly concentrated hydrate accumulations. The evolution of fluid expulsion controlling hydrate formation offshore southern Hikurangi is described in stages during which different methane sources (in situ, buried and thermogenic) have been dominant.

[1] Multibeam swath bathymetry data from the southwest margin of the Chatham Rise, New Zealand, show gas release features over a region of at least 20,000 km2. Gas escape features, interpreted to be caused by gas hydrate dissociation, include an estimated a) 10 features, 8–11 km in diameter and b) 1,000 features, 1–5 km in diameter, both at 800–1,100 m water depth. An estimated 10,000 features, ∼150 m in diameter, are observed at 500–700 m water depth. In the latter depth range sub-bottom profiles show similar gas escape features (pockmarks) at disconformities interpreted to mark past sea-level low stands. The amount of methane potentially released from hydrates at each of the largest features is ∼7*1012 g. If the methane from a single event at one 8–11 km scale pockmark reached the atmosphere, it would be equivalent to ∼3% of the current annual global methane released from natural sources into the atmosphere.

[1] P-wave velocity increases in the presence of gas hydrates and decreases in the presence of free gas in the sediments, making it an excellent means to investigate gas hydrate systems. However, seismic velocity is typically derived from surface seismic data without consideration of seismic anisotropy. The presence of anisotropy in the hydrate bearing sediments adds an additional complexity in data analysis; however anisotropy can help reveal the distribution of hydrates. Here we report on the evidence of seismic anisotropy at Hydrate Ridge along the Cascadia convergent margin. We find that the south summit is anisotropic, while the basin side (east of south summit) is isotropic. Anisotropy is likely caused by the hydrate veins. We interpret the anisotropy parameters in terms of the distribution and fabric of gas hydrates.

[1] The dissociation of gas hydrates in sediment pores is thought to decrease seafloor strength potentially facilitating submarine slides because of the generation of overpressured gas and the “melting” of load-bearing or cementing solid hydrate. Here, we present findings that suggest gas hydrates may lead to a previously unknown mechanism of seafloor erosion. Gas-hydrate-bearing sub-sea ridges on the Hikurangi Margin are eroded close to the top of the hydrate stability field in the ocean. We hypothesize that gas hydrate instability may lead to ridge erosion by a combination of two processes; hydrate destabilization caused by depressurization during ridge uplift, and repeated pore volume expansion and contraction from hydrate dissociation and formation triggered by fluctuating water temperatures.

[1] We present results from an analysis of anisotropy in marine sediments using walkaway vertical seismic profiles from the Blake Ridge, offshore South Carolina. We encountered transverse isotropy (TI) with a vertical symmetry axis in a gas-hydrate-bearing unit of clay and claystone with Thomsen parameters e = 0.05 ± 0.02 and δ = 0.04 ± 0.06. TI increased to e = 0.16 ± 0.04 and δ = 0.19 ± 0.12 in the underlying gas zone. Rock physics modeling suggests that the observed TI is caused by a partial alignment of clay particles rather than high-velocity gas-hydrate veins. Similarly, the increase of TI in the gas zone is not caused by thin low-velocity gas layers but rather, we speculate, by the sharp contrast between seismic properties of an anisotropic sediment frame and elongated gas-bearing pore voids. Our results underscore the significance of anisotropy for integrating near-vertical and wide-angle seismic data.

[1] Subsidence of Lima Basin, part of the Peruvian forearc, is controlled by tectonic erosion by the subducting Nazca plate. Multichannel seismic reflection data coupled with age and paleowater depth constraints derived from Ocean Drilling Program (ODP) coring now allow the rates of erosion to be reconstructed through time. In trenchward locations the forearc has experienced limited recent relative uplift (700–850 m) likely due to preferential basal erosion under the center of Lima Basin. Long-term subsidence driven by basal tectonic erosion dominates and is fastest closest to the trench. Since 47 Ma (Eocene) up to 148 km of the plate margin have been lost at an average rate of up to 3.1 km myr−1. Appoximately 110 km of that total appears to be lost since 11 Ma, implying much faster average rates of trench retreat (10 km myr−1) since collision of the Nazca Ridge with the Lima Basin at 11 Ma. Although there is no clear subsidence event at ODP Site 679 during the time at which Nazca Ridge was subducting beneath this part of the forearc (4–11 Ma), the more trenchward ODP Sites 682 and 688 show significant deepening after 11 Ma indicating that subduction of the ridge accelerates tectonic erosion. Long-term rates of crustal erosion in the region of Lima Basin are greater than estimates of regional arc magmatic productivity, implying that such margins are net sinks of continental crust.

Results from Ocean Drilling Program Leg 112 indicated tectonic erosion of the Peruvian convergent margin during subduction of the Nazca Ridge followed by renewed accretion against the erosional scar. Seismic images and swath mapping morphology show a regional geology shaped as subduction of Nazca Ridge migrated 800 km along the Peru Trench. The accretionary prism is imaged at two stages of development since upper Miocene time with depth-migrated seismic records in which the “back stop” is well defined. From these data, material flux can be quantified. After the ridge crest subducted and accretion dominated over erosion, the prism grew rapidly to 10- to 15-km width. During rapid growth, the margin taper was large and about 60 percent of the sediment supply accreted. As the prism growth slowed, only 30 percent accreted despite an increased trench sediment supply. Since the convergence rate changed little, the inverse relation between sediment supply and prism growth rate suggests that other processes allow most of the sediment input to be subducted. The decreased prism growth may be partially controlled by varying the structure of the back stop. The force required to activate faults cutting across the back stop is probably greater than the force allowed by basal friction across the plate boundary. Thus the shear stress transmitted across the plate boundary is insufficient to drive “out-of-sequence” thrust faults that would thicken the back stop or raise it beyond a limited height. A constant back stop height limits accretion, and the subduction window appears to open and accommodate increased trench sediment beneath the back stop rather than inducing a transport of sediment over the top of the back stop. Nazca Ridge subduction affected the convergent margin mass flux for about 8 Ma. Mass balancing indicates an order of magnitude increase in the current sediment volume subducted in the north compared to that near the ridge crest. When growth of the accretionary prism slowed, a greater terrigenous input through the trench was accommodated by sediment subduction. The resulting range of material flux is estimated to increase the sediment interlayer between the plates about 500 m which may affect coupling.