==== Front Sci Rep Sci Rep Scientific Reports 2045-2322 Nature Publishing Group UK London 78958 10.1038/s41598-020-78958-z Article Gas hydrate versus seabed morphology offshore Lebu (Chilean margin) Vargas-Cordero Iván lacruzvargas@gmail.com 1 Tinivella Umberta 2 Villar-Muñoz Lucía 34 Bento Joaquim P. 5 Cárcamo Carolina 6 López-Acevedo Diego 5 Fernandoy Francisco 6 Rivero Alessandra 6 Juan Marion San 67 1 Valparaiso, Chile 2 grid.4336.20000 0001 2237 3826Istituto Nazionale di Oceanografia e di Geofisica Sperimentale - OGS, 34010 Trieste, Italy 3 Centro de Investigación y Formación San Ignacio de Huinay, Hualaihué, Chile 4 grid.443909.30000 0004 0385 4466Departamento de Geofísica, Facultad de Ciencias Físicas y Matemáticas, Universidad de Chile, Santiago, Chile 5 grid.8170.e0000 0001 1537 5962Escuela de Ciencias del Mar, Pontificia Universidad Católica de Valparaíso, Valparaíso, Chile 6 grid.412848.30000 0001 2156 804XUniversidad Andres Bello, Viña del Mar, Chile 7 grid.443909.30000 0004 0385 4466Universidad de Chile, Santiago, Chile 14 12 2020 14 12 2020 2020 10 2192829 6 2020 1 12 2020 © The Author(s) 2020Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.Gas-hydrate occurrences along the Chilean margin have been widely documented, but the processes associated with fluid escapes caused by the dissociation of gas hydrates are still unknown. We report a seabed morphology growth related to fluid migration offshore Lebu associated with mud cones by analysing oxygen and deuterium stable water isotopes in pore water, bathymetric, biological and sedimentological data. A relief was observed at − 127 m water depth with five peaks. Enrichment values of δ18O (0.0–1.8‰) and δD (0.0–5.6‰) evidenced past hydrate melting. The orientation of the relief could be associated with faults and fractures, which constitute pathways for fluid migration. The benthic foraminifera observed can be associated with cold seep areas. We model that the mud cones correspond to mud growing processes related to past gas-hydrate dissociation. The integration of (i) the seismic data analysis performed in the surrounding area, (ii) the orientation of our studied relief, (iii) the infaunal foraminifera observed, (iv) the grain size and (v) the total organic matter and isotope values revealed that this area was formerly characterised by the presence of gas hydrates. Hence, this part of the Chilean margin represents a suitable area for investigating fluid-migration processes. Subject terms Climate sciencesEnvironmental sciencesNatural hazardsSolid Earth scienceshttp://dx.doi.org/10.13039/501100002850Fondo Nacional de Desarrollo Científico y Tecnológico11140214Vargas-Cordero Iván Italian Ministry of Education, Universities and Research under the extraordinary contribution for Italian participation in activities related to the international infrastructure PRACE—The Partnership for Advanced Computing in EuropeDecreto MIUR No. 631 dd. 8 August 2016Tinivella Umberta Oceano High-Performance Computing Cluster at PUCV, ChileFONDEQUIP Nº EQM170214issue-copyright-statement© The Author(s) 2020 ==== Body Introduction Along continental margins, morphological features associated with fluid escapes such as mud volcanoes, mud mounds, pockmarks, seeps, precipitates—carbonates and hydrates-, piping/rills and brine pools have been reported worldwide (e.g.1–7). Fluid escape morphologies related to gas escape include a) seafloor depressions as pockmarks and b) seafloor mud growing as mud volcanoes, mud cones, mud diapirs and mounds. Although numerous fluid escape features have been identified in diverse tectonic settings, there is still little information on the distinct processes and factors that form and control them. Some factors are evident, such as the chemistry of expelled fluids, and the fluid flow rates and duration. As a first approach, the accumulation of sediment in the form of mud growing, such as mud volcanoes, has been related to the rapid expulsion of fluids under overpressure. Pockmarks are associated with gas expulsion in near-surface sediments, whereas slow seepage promotes lithification of the seafloor through precipitation of a variety of mineral species (e.g.8–11). Some authors describe mud growing morphologies as mud volcanoes and mounds when they are isolated, and as mud cones and mud diapirs when they are aligned in space12,13. Others indicate that when the extruded sedimentary material is responsible for the seafloor morphology, it can be described as mud volcanoes and mud cones14, whereas when the fluidised material is not expelled, it can be defined as mud diapir1. Fluid escapes can be formed mainly by microbial and thermogenic methane gas and water15,16. This gas can give place to gas-hydrate formation in marine sediments if pressure and temperature conditions are adequate17; in this condition, the gas is trapped in a lattice of water molecules. Along continental margins, gas hydrates occur naturally within the Gas Hydrate Stability Zone (GHSZ), mostly at sub-marine sediment depths greater than 300–500 m with low temperature, high pressure and adequate amounts of sedimentary organic carbon (2–3.5%), where enough methane is present18,19. Moreover, salinity, gas composition, geological structure, fluid migration, and pore space of marine sediments influence gas-hydrate formation18–20. Gas-hydrate occurrences along the Chilean margin are distributed from 33 to 57°S as widely reported in literature21–31. On the other hand, only a few studies have documented fluid escapes related to gas-hydrate dissociation through faults and fractures (e.g.32,33). Identifying areas where gas-hydrate dissociation occurs is essential. The presence of methane, a potent greenhouse gas34, in zones of shallow fluid escapes could may contribute to (a) an increase in surface temperature in favour of global warming, (b) change in the physic-chemical conditions of seawater, (c) affect the marine microfaunal diversity pattern, and (d) affect the nucleation and the rupture propagation of earthquakes, related to slides and tsunamis35–37. The most common techniques used to identify gas-hydrate dissociation zones are biological, geochemical and geophysical analyses techniques. Biological indicators such as benthic foraminifera, bivalve and microbial communities, were related to fluid escapes38–42. For example, foraminiferal taxa reported worldwide, including Uvigerina sp., Bolivina sp., Chilostomella sp., Globobulimina sp., Quinqueloculina sp. and Nonionella sp., were associated to cold seep occurrences and the presence of methane43–45. These organisms can live in organic-rich and reducing environments, where food availability is high43–45. In different regions including the Chilean coast46–48, reported several cases of enriched stable water isotope values (as measured from pore water in marine sediments) linked to gas-hydrate dissociation. The reported values for δ18O and δD reached up to 3‰ and 10‰, respectively. Geophysical studies documented seabed morphologies in relation to fluid leakages, using acoustic data and high-resolution images4,49,50. In south-central Chile, the simultaneous interpretation of direct and indirect geophysical data is essential to determine the active structural domain located in the vicinity of the Arauco basin51,52. These geophysical data made it possible to recognise the exchange of subsidence and uplift events responsible for the accretion and erosion of the prism53,54. The geophysical data interpretation indicated that the sedimentary sequences are quite complex and distinguished by marine and continental deposit cycles. These environments are compatible with the presence of gas and carbon deposits detected in the past55,56. This study aims to understand whether past gas hydrate dissociation could be related to the morphological growth observed on the Lebu coast (Chilean margin). To address our objective, we used direct and indirect evidence of the presence of gas hydrate dissociation, including geophysical, biological and geological sampling as well as physical–chemical and geochemical analyses. Materials and methods The data, methods and techniques used in this work were based on a multi- and inter-disciplinary approach to characterize the morphology and the fluid escapes; field and laboratory data, theory and modelling approaches, were used. The approach included bathymetric data processing and sedimentological, physical–chemical, geochemical and biological analyses of seawater and marine sediment samples. Data sampling In the framework of the project entitled “Identification and quantification of gas emanations associated with gas hydrates (FONDECYT 11140214)”, sedimentological, geochemical and bathymetric studies offshore Lebu were performed (Fig. 1). In 2016 and in 2017, two oceanographic cruises onboard R/V Kay Kay II were carried out to collect bathymetric data, seawater samples and marine sediments.Figure 1 Location map of the studied area. The red star shows corer recovery and bathymetric surveys. The black dashed line shows the bathymetric profile used in Fig. 8, the white solid lines indicate faults and fractures described by several authors53,92,93 and the red dashed lines show embayments related to ancient slides reported by92. Marine sediment samples were collected about 127 m below the sea level (mbsl), using a gravity corer (GC-02; 9 cm diameter) that drilled 240 cm into the seabed. The corer was deplyed at the northern of the positive relief (at a distance of approximately 100 m; 73°44′ 25″ W/37°36′ 10″ S; Fig. 2), and the sediment-core sample was divided into four sections of 60 cm long. Each sediment-core section was frozen on board for later analysis.Figure 2 Bathymetric map indicating location corer GC-02 (red circle), water samples (white circles) and identified mud cones (MC1, MC2, MC3, MC4 and MC5). In (A) and (B) 3D images with orientation NW and SW respectively, while in (C), (D), (E) and (F) zooms for mud cones are reported. Water samples were collected using Niskin bottles at 0 m, 10 m, 20 m, and 50 m below sea surface, and above the seafloor. Temperature, conductivity, dissolved oxygen, and pH, were determined using the multiparameter Meter (IP67, model 8602). These parameters were measured at the two ends of the identified lineament (Fig. 2). The measurements were obtained near the relief to evaluate the relationship between marine sediments and the water column in the presence of ancient gas hydrates. Bathymetric data were acquired using a multibeam Reson SeaBat 7125 echosounder (400 kHz, 0.5° × 1°). Sound velocity data were acquired using a SVP90 probe, and an AML Oceanographic Model Minos X sound velocity profiler. Preliminary processing was performed on board using a PDS2000 commercial software allowing real-time bathymetric data processing using the SVP90, the AML information, and the ship motion (pitch, roll, yaw and heave). The PDS2000 software used real-time information from other sensors to calculate the Total Propagated Error (TPE), specifically the vertical base reduction, and total heave and depth. Each sounding was filtered by setting the International Hydrographic Organization (IHO) Order 2. Methods The acoustic data processing was performed using the open-source MB-System libraries57, to model the bathymetry of the region, with tide and scattering effects attenuated. Bathymetric grids were created with the nearest neighbour interpolation algorithm, using the open-source software Generic Mapping Tools (GMT;58). The algorithm constructs regularly gridded values, in which each node corresponds to the weighted mean of the samples around a search circle of 5 m. The selected grid was configured with a spatial resolution of 1 m. A squared median filter of 5 m was applied to smooth the grid. Finally, the bathymetric accuracy was off 0.2% in-depth. These values are in agreement with the accuracy values reported by other authors59. The standard deviation for depth ranged from 0.09 m to 0.27 m for 128.59–133.87 m of depth. Grain size analysis was performed by sieving method where sediments passed (by agitation) through a range of different-size meshes. Fifty grams of sediments were sieved using the following mesh sizes: 60, 80, 120 and 230. Alternatively, the pipette method was adopted to separate clay and silt fractions by selecting 15 g of mud sample. Statistical parameters were calculated in agreement with reported formulas60–62. Surface to seafloor water physical–chemical properties of temperature, pH, salinity, and dissolved oxygen to near the positive relief, were obtained using the multiparameter Meter. Temperature was measured in Celsius degrees, with an accuracy of ± 0.5 °C, and pH was directly related to the concentration ratio of hydrogen ions [H +] and hydroxyl [OH]63 with an accuracy of ± 0.1. Salinity was obtained from conductivity, with an accuracy of ± 0.1, dependent on the number of dissolved ions per unit volume and on the mobility of the ions. Dissolved oxygen was measured in % and in mg/L, with an accuracy of ± 3%. The sediment-core sample was cut in sections of 10 cm long, and physical–chemical parameters were measured, including water content (W in %), porosity (Φ in %), the content of solid material per unit volume, expressed as apparent density (ρ in g/cm3;64) and Total Organic Matter (TOM in %). Sediment samples were dried in a forced-air oven at 60 °C for 36 h and in a desiccator for 30 min. TOM content was measured by gravimetric determination of weight loss through the loss-on-ignition method65,66. Two grams of dry sediment sample were calcined in a muffle at 500 °C for 5 h and placed in a desiccator for 30 min until stable weight was obtained, in order to reduce the associated error. For foraminiferal sampling and identification, the sediment-core sample was cut into sections of 15 cm of which 50 g of material was extracted. The material was washed, dried, and sieved using 120 and 230 diameter sieves. Specimens were placed in Petri dishes and observed under binocular magnification. General morphological features were characterized using the Atlas of Benthic Foraminifera67, and the genus were identified based on the study of Chilean material (e.g.68). Pore water from the sediment-core samples was extracted using an ACME lysimeter (0.2 μm) to analyse oxygen and deuterium stable water isotopes. Pore water extraction procedure involved cutting the sediment core in sections of 5 cm long, centrifugation, and pore water extraction using Rhizon MOM with pore sizes ranging from 0.12 µm to 0.18 µm. Stable water isotope was determined by Cavity Ring-Down Spectroscopy (CRDS) method. Deuterium water isotope and oxygen content were evaluated using in-house standards LIMS69, normalised to the VSMOW-SLAP scale, and values were reported as δ-values for deuterium (δD) and oxygen (δ18O). Each sample was measured at least twice on different days. For each measurement, samples were analysed for five consecutive times. Results were accepted if the standard deviation of each single run (composed of five repetitions) was < 1‰ for δD and < 0.1‰ for δ18O. The accepted stable water isotope value of a sample was the mean of (at least) 2 different valid measurements within the range of the explained standard deviation. Results Morphology Bathymetric data revealed a positive relief with orientation N55°E at 127 mbsl. The relief showed a mean elevation of about 6 m above the seafloor, an extension of 410 m length and a width of 50 m, corresponding to an area of 14,640 m2 (Fig. 2). The relief was composed of five mud growing. They were aligned (see Introduction;12,13), thus defined as mud cones (see Fig. 2 MC1, MC2, MC3, MC4, MC5). They were characterized by peaks ranging from 3 to 9 m height, with MC5 reaching the maximum value of 9 m in height. MC1 showed 2 elongated shape cones of about 40 m length and 20 m width, with steep and smooth sides. MC2 displayed elongated and circular shape cones of about 80 m length and 50 m width, mainly exhibiting steep sides. MC3 and MC4 presented rectangular shapes cones of 60 m length and 50 m width, and their morphologies varied from smooth to steep sides. MC5 showed elongated and circular shape cones with smooth and steep sides. Physical–chemical parameters Grain size analysis showed constant values with depth. The mean grain size corresponded to sandy mud textural group, characteristic to mud cones. Silt comprised 60% of the total sediment-sample volume (Fig. 3).Figure 3 Grain size distribution in marine sediments (corer GC02). A slight variation of water content (W), ranging from 38.6 to 46.3% (mean 43.1%), porosity (ɸ), ranging from 62.7 to 69.7% (mean 66.9%) and an apparent density (ρ), ranging from 1.5 to 1.7 g/cm3 (mean 1.6 g/cm3), were detected (Table 1). TOM values showed a variable trend with a maximum value of 8.7% of total volume located at 2.2 mbsl, while the minimum value of 5.1% of total volume detected at 0.4 mbsl (Fig. 4). An opposite trend distribution was observed between porosity and apparent density.Table 1 Physical–chemical parameter distribution in marine sediments. Depth (m) W (%) ɸ (%) Ρ (g/cm3) TOM (%) 0.1 45.2 68.8 1.56 7.9 0.2 42.2 66.1 1.61 7.2 0.3 41.2 65.2 1.62 6.6 0.4 41.3 65.3 1.62 5.1 0.5 39.9 64.0 1.64 5.9 0.6 38.6 62.7 1.67 6.0 0.7 40.3 64.4 1.64 5.9 0.8 42.7 66.5 1.60 6.4 0.9 43.0 66.8 1.60 6.2 1 42.8 66.6 1.60 6.5 1.1 42.6 66.5 1.60 6.6 1.2 42.9 66.7 1.60 6.1 1.3 45.0 68.6 1.57 6.0 1.4 45.2 68.8 1.56 6.7 1.5 45.1 68.7 1.56 4.2 1.6 46.3 69.7 1.55 7.5 1.7 44.1 67.8 1.58 6.7 1.8 45.3 68.8 1.56 7.1 1.9 40.7 64.7 1.63 5.4 2 43.7 67.5 1.58 7.0 2.1 44.3 68.0 1.58 6.8 2.2 43.3 67.1 1.59 8.7 2.3 45.4 68.9 1.56 6.9 2.4 44.4 68.0 1.58 6.5 Mean 43.1 66.9 1.59 6.5 Minimum 38.6 62.7 1.55 4.2 Maximum 46.3 69.7 1.67 8.7 Figure 4 Physical–chemical parameters distribution in marine sediments (corer GC02). Isotope analysis Pore stable water isotope analysis of the sediment-core samples showed positive values ranging from 0.0 up to 1.8‰ for δ18O and 5.6‰ for δD (Fig. 5). Stable water isotope δ-values were close to 0.0 at the sediment–seawater interface, and the values increased with depth of the sediment-core samples (i.e. enrichment), with low variability (δ18O, 0.33 Std. Dev.; δD 0.95 Std. Dev.). No negative δ-value was observed in any of the sediment-core samples analysed.Figure 5 Oxygen (δ18O) and deuterium (δD) stable water isotope distribution in sediment from: (a) Depth profile of the corer; (b) Cross-isotope linear regression of pore water samples and (c) Cross-isotope relationship of pore water samples against the global meteoric water (GMWL). Biological communities Benthic foraminiferal aggregations were uncovered in the top section of the sediment-core samples (0–60 cm), comprising globose and elongated specimens. The following taxa where identified, including opportunistic species: Globobulimina sp., Bolivina plicata, Anomalinoides sp., Uvigerina peregrina, Oridorsalis tener, and Quinqueloculina vulgaris (Fig. 6).Figure 6 Benthic foraminifera. In (1a) Globobulimina sp., lateral view (10×); (1b) Globobulimina sp., lateral view (10×); (2) Bolivina plicata., lateral view (5×); (3) Anomalinoides spp., lateral view (5×); (4) Uvigerina peregrina, lateral view (5×); (5a) Oridorsalis tener, lateral view (5×); (5b) Oridorsalis tener, lateral view (5×); (6) Quinqueloculina vulgaris, lateral view (10×). Water properties Water column temperature decreases from 14 °C at the surface to 12 °C at 20 mbsl. Salinity and dissolved oxygen content also decreased with depth, from 31 to 33‰, and 60–66.2% were minimum at 0.6 mbsl. pH values ranged from 7.5 to 8.1 from surface to depth. Discussion Stable water isotope composition of pore water provides strong evidence of gas-hydrate dissociation. Heavy oxygen isotope enrichment with depth can be explained only by isotopic fractionation after hydrate dissociation (e.g.70). The fractionation process occurs during hydrate formation, which concentrates the heavy isotopes (e.g. δD and δ18O) in the hydrate layer, as a result of the dissociation pressure of the ice-lattice Thus, hydrate dissociation yields pore fluids with more positive δ18O values by comparison to seawater. Figure 5a presents the stable water isotope profile of the entire sediment-core sample, showing a clear increase in isotope value with depth, with values close to 0‰ at the seawater-sediment interface to positive values at the bottom of the corer. In addition, observational data at similar latitudes, and modelled surface water stable isotope composition for this ocean region, showed that shallow water tend to present slight negative isotope composition (∼− 0.2‰ to − 0.5‰ δ18O;71,72), which are related from the transport of Sub Antarctic Waters through the Humboldt Current System along the Chilean coast73. Negative isotopic values result mainly from the mixture of oceanic and depleted meltwater from the Antarctic Ice Sheets74. In this study, the isotopic trend was indicative of seawater mixing at the top of the sediment-core sample, and of different water source at the bottom of the core. The cross-isotope relationship of our samples shows that the stable water isotope composition of pore water has a strong positive correlation (e.g. simultaneous enrichment of δ18O and δD; Fig. 5b). Positive δ18O values are associated to clay minerals dewatering, generally related to a robust δD decrease47. Our results showed an increase in δD, as was also observed by75 in the same relief, agreeing with the hypothesis of past hydrate melting46,47,76–78. Note positive values of meteoric waters are negligible, as shown in Fig. 5c, because porous waters deviate from the global meteoric water line. The infaunal foraminifera, found in the shallower sediment samples (e.g. Bolivina sp., Globobulimina sp., Uvigerina sp.), may be associated with modern cold seeps. These taxa can metabolize seeping methane, directly or indirectly exploiting the available geochemical energy source79, as documented by several authors in cold-seep area (e.g.44,80,81), including the continental margins of California, Japan and Mexico36,43–45,82. Benthic communities are often found in high organic-content ambient and low oxygen environments, characteristic to cold seep44,45. In this study, TOM was 6.5%, suggesting that the organic-carbon stock offshore Lebu promoted the development of benthic foraminifera, as demonstrated by83 who found a positive correlation between foraminiferal distribution and at least of 1.5% TOM. Grain-size data provide crucial information about hydrodynamic conditions; in particular, mud and sand content associated to coastal and beach systems and riverine or deltaic deposits84. Recently85, indicated that grain size decreases generally with the pore size, and that sediment permeability decreases with decreasing pore size. Consequently, the formation of gas hydrate is limited in low-grain sized sediments, because of reduced flow of gas and water to the sediments. Clay- and silt-particle agglomerates may lead to an increase in both apparent grain size and pore throat size86, maintaining the relative permeability of the sediments, which promotes the formation of gas hydrate in silt-rich sediments. The depth trend of the physical–chemical parameters measured in this study (W, ɸ, ρ, TOM), in relation to grain size, permitted to describe the benthic environment offshore Lebu. Our approach was particularly relevant given the relationship between clay and silt content and high TOM values87, which are favourable to gas-hydrate formation (e.g.88,89). Seawater temperature, salinity, dissolved oxygen and pH were typical to the study region75. Water-column temperature decreased with depth and salinity and dissolved oxygen were inversely proportional63. Seawater alkalinity generally assumes pH values between 7.4 and 8.4, as obtained in our sample. Gas phases concentrations were estimated across the continental slope zone off Lebu by25, reporting 15% and 0.2% of total volume for hydrates and for free gas, respectively. As pointed out by90, the high proportion of gas hydrate amount may be generated by advective flow of dissolved hydrate or gas from depth, Other studies suggest that lateral fluid migration occurs from deep levels through faults and fractures canalising fluids that lead to the formation of mud cones and mud volcanoes (e.g.32,33). Several researchers reported faults in the proximity of our study area (Fig. 1), such as the Santa María fault of a similar orientation to the relief document in this study (N55°E;51,91,92). Gas accumulations can reach shallow areas because the base of the gas-hydrate stability zone can be very shallow on continental shelves, as indicated by theoretical modelling25. Therefore, the integration of (i) the seismic data analysis performed in the surrounding area25, (ii) the orientation of our studied relief, (iii) the infaunal foraminifera observed, (iv) the grain size and (v) the TOM and isotope values reported here, suggest that this area was characterised by the presence of gas hydrate. The gas hydrates dissociation and fluid expulsion in this region may be related to (a) regional uplift occurring since late Pleistocene-Holocene52,93,94, (b) deglaciation processes, (c) variation in heat flow caused by landslides, and (d) the vertical variability of the equatorial subsurface water mass (ESSW), transported poleward by the Peru–Chile undercurrent (PCUC) and the El Niño-Southern Oscillation (ENSO). In the Arauco basin, several authors reported a surface uplifting of about 1.5 km during the Middle Pliocene51,52. Basal accretion processes may be responsible for the characteristic uplifting of the Arauco region coast and shelf systems93,95. The constant uplifting in this area may decrease the hydrostatic pressure of the sediments and dissociate the gas-hydrate layer as reported in other areas worldwide96–99. On the other hand, deglaciation processes starting, which started 20,000 years ago in Chile100–102, lead to warmer seabed temperatures in shallow waters promoting gas-hydrate dissociation. Off the Arauco Peninsula, ancient landslides (see Fig. 1;103) may increase the heat flow as a result of fluid advection area or increased geothermal gradient29. Heat-flow rise also promotes hydrate dissociation, as reported in Mocha Island close to our study area104. Finally, in this region105, find the subsurface oceanic poleward flow, the so-called PCUC, originated in subtropical regions (6°S). This current flows from surface water to 150 mbsl, transporting the ESSW warm temperatures, high salinity, high nutrients, and low oxygen content along the continental shelf and slope off Chile up to 48°S73,106. Essential physical processes that drive the vertical variability of ESSW include mesoscale intrathermocline eddies that propagate westward105, and their transport of volume that is significantly correlated with the variability of PCUC transport of water, which in turn is forced by the ENSO equatorial signal107. Intra-seasonal coastal trapped waves and Rossby waves result mainly from the teleconnection between the equatorial Pacific and South America, particularly during El Niño events108–110, and other high-frequencies coastal internal waves forced locally111. Temperature data sampled 28 years apart revealed decadal warming of the middle to the deep-ocean layers close to the study area112, which may explain gas-hydrate dissociation in this study17,113. To understand the depth stability of gas hydrate, the theoretical base of the GHSZ was calculated assuming a geothermal gradient of 30 °C/km (as per25), and a mixture of 95% methane and 5% ethane as per measures obtained at Ocean Drilling Program Site 1235114. The theoretical base of the GHSZ was calculated as the intersection between the hydrate stability curve and the temperature/pressure curve in the sediments (e.g.115). The first curve was evaluated by using the Sloan equations17, used to model a mixture of gases in freshwater116, equations were used to shift the freshwater hydrate curve because of water salinity effect28, 3.5% in our case. The second curve was evaluated considering water density of 1040 kg/m3, as per28. It is crucial to note that in our study area, hydrate and free gas data were detected by seismic analysis, confirming the site is characterised by upward fluid flow25,93. Figure 7 shows the main geological features, including the seismic indicator of the transition between the gas hydrate and the free gas zones, i.e. bottom simulating reflector (BSR), detected in a seismic line located near the relief.Figure 7 Line drawing section corresponding to the SO161-44 seismic stacking section modified from Fig. 3 in94, in which the main geological features and the seismic reflector indicating the transition between gas hydrate and free gas (e.g., the bottom simulating reflector, BSR) are indicated. In the inset, the location map showing static model location reported in Fig. 8 (red dashed line), and mud growing zone (red circle) and the SO161-44 position (black line). Figure 8 depicts the base of GHSZ at about 400 m water depth, indicating that at shallower depth, hydrate is not stable, and only free gas is present. Our study area is located in a narrow continental shelf (15 km width), favouring fluids advection associated with gas hydrate dissociation; thus, gas accumulations migrate from the base of GHSZ to shallow areas.Figure 8 Schematic profile explaining mud growing formation (in red). The profile location is shown in Fig. 1. Dashed lines show theoretical bases of GHSZ by using the geothermal gradient of 30 °C/km for several scenarios supposing that the hydrate is formed by a mixture of 95% methane and 5% ethane. The blue dashed line indicates the actual theoretical base of the GHSZ. The dotted lines indicate the theoretical base of GHSZ supposing a decrease of the bottom temperature of 2 °C (black dotted line), 3 °C (magenta dotted line), 4 °C (green dotted line) and 5 °C (red dotted line). The solid black line indicates the seafloor. The pink arrows indicate the direction of the fluid/mud outflow. Possible faults and fractures are also reported as black lines. The rectangle in (a) shows the zoom, in which the modelling for the theoretical base of GHSZ reaches the continental shelf. At higher latitudes, a substantial reduction of the GHSZ was observed due to warming in the last 20,000 years (e.g.33,117,118). To verify such trend in our study area, we modelled the theoretical base of the GHSZ, assuming past temperature conditions as reported by paleoclimate reconstruction studies117,119, e.g. decrease of seabed temperature by 1 °C, 2 °C, 3 °C, 4 °C, and 5 °C (Fig. 8). As per117, the modelled temperatures were characteristic of about 11,000, 12,000, 16,000, 17,000 and 18,000 yr BP respectively. It was assumed that other parameters necessary to assess hydrate stability (i.e. geothermal gradient, water depth, gas composition) remained invariant. The model indicated that the origin of the mud structures observed in this study may be related to hydrate dissociation caused by past increase of seabed temperature in around 12,000 yr BP. Additional measurements are necessary to test our hypothesis. Alessandrini et al.113 recently acknowledge that this region of the Chilean margin may be critically disturbed in the long term, because of the tectonic-sedimentary configuration. Their model indicated that in the next 100 years, about 6.5% of the area where the gas hydrate is stable possibly will be affected by hydrate dissociation, induced by global warming. Hydrate dissociation may affect the Chilean margin seafloor morphology that is located close to the shoreline (less than 10 km)113. Conclusions The positive relief indicated in our multidisciplinary study may be associated to mud cones by fluid-flux supply, which may be channelized through faults and fractures, as detected by seismic data93. δ18O and δD-enrichment of pore water, resulting from gas-hydrate melting and dissociation, actively support this observation. We suggest that the dissociation of gas hydrates, by tectonic uplift (i.e. decreasing pressure) and/or climate change (i.e. increasing temperature), generates a mixture of water enriched with heavy isotopes, mud, and gas, which may be expelled through the multiple faults present in our study area offshore Lebu. The discharge seems to be intermittent in time, like pulses, supporting the particular morphology of the monticules observed in the area. The region undergoes rapid hydrate dissociation, generating substantial carbon fluxes in short period (decades), as reported in other shallow areas, such as the Gulf of Mexico and Arctic (i.e.120). Sediment-grain size analysis provided information about (i) the type of sedimentary material that migrates from deeper to shallower zones, (ii) the anomalies of heavy isotopes concentrated in pore water, and (iii) the ideal environment sustaining chemosynthetic benthic organisms in this area, such as the foraminiferal genera Bolivina, Uvigerina, and found in shallow sediments. Our research provides new information about gas-hydrates dissociation in shallow seabed features, i.e. mud cones on the continental shelf. It is critical to emphasise that fluid escapes in shelf seas may spread to the atmosphere releasing methane, thus contributing to global warming (e.g.121). Our results demonstrate that Chilean margin offshore Lebu is a suitable region for investigating such processes. Publisher's note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations. Acknowledgements The authors are grateful to Michela Giustiniani and Daniela Lazo for constructive discussions and useful comments. Special thanks to Mauricio and Daniel from the palaeontology laboratory (UNAB—Viña del Mar), who helped us with the foraminifera identification. This research was partially supported by the Italian Ministry of Education, Universities and Research (Decreto MIUR No. 631 dd. 8 August 2016) under the extraordinary contribution for Italian participation in activities related to the international infrastructure PRACE—The Partnership for Advanced Computing in Europe (www.prace-ri.eu). The authors also acknowledge computing time on the Oceano High-Performance Computing Cluster at PUCV, Chile (FONDEQUIP No. EQM170214). Special thanks to Noëlie Benoist of the National Oceanography Centre Southampton for her kind help with the grammar. This research was partially supported by the Italian Ministry of Education, Universities and Research (Decreto MIUR No. 631 dd. 8 August 2016) under the extraordinary contribution for Italian participation in activities related to the international infrastructure PRACE—The Partnership for Advanced Computing in Europe (www.prace-ri.eu). The authors also acknowledge computing time on the Oceano High-Performance Computing Cluster at PUCV, Chile (FONDEQUIP No. EQM170214). Special thanks to Noëlie Benoist of the National Oceanography Centre Southampton for her kind help with the grammar. Author contributions All authors were involved in the data processing and preparation process. All authors were involved in the discussion and revision process with section leads as follows: I.V.C, U.T., L.V.M., J.P.B., and C.C., Introduction-Discussion-Conclusions; I.V.C, U.T., L.V.M., J.P.B., C.C. , D.L.-A., F.F., A.R. and M.S.J., Materials and Methods-Results. Competing interests The authors declare no competing interests. ==== Refs References 1. Milkov AV Worldwide distribution of submarine mud volcanoes and associated gas hydrates Mar. Geol. 2000 167 29 42 10.1016/S0025-3227(00)00022-0 2. Dimitrov LI Mud volcanoes: the most important pathway for degassing deeply buried sediments Earth Sci. Rev. 2002 59 49 76 10.1016/S0012-8252(02)00069-7 3. Van Rensbergen P De Batist M Klerkx J Hus R Poort J Vanneste M Krinitsky P Sublacustrine mud volcanoes and methane seeps caused by dissociation of gas hydrates in Lake Baikal Geology 2002 30 631 634 10.1130/0091-7613(2002)030<0631:SMVAMS>2.0.CO;2 4. Loncke L Mascle J Parties FS Mud volcanoes, gas chimneys, pockmarks and mounds in the Nile deep-sea fan (Eastern Mediterranean): geophysical evidences Mar. Pet. Geol. 2004 21 669 689 10.1016/j.marpetgeo.2004.02.004 5. Hovland M Svensen H Forsberg CF Johansen H Fichler C Fosså JH Rueslåtten H Complex pockmarks with carbonate-ridges off mid-Norway: products of sediment degassing Mar. Geol. 2005 218 191 206 10.1016/j.margeo.2005.04.005 6. Lykousis V Alexandri S Woodside J De Lange G Dählmann A Perissoratis C Rousakis G Mud volcanoes and gas hydrates in the Anaximander mountains (Eastern Mediterranean Sea) Mar. Pet. Geol. 2009 26 854 872 10.1016/j.marpetgeo.2008.05.002 7. Chen SC Hsu SK Tsai CH Ku CY Yeh YC Wang Y Gas seepage, pockmarks and mud volcanoes in the near shore of SW Taiwan Mar. Geophys. Res. 2010 31 133 147 10.1007/s11001-010-9097-6 8. Hovland M Judd AG Seabed Pockmarks and Seepages: Impact on Geology, Biology and the Marine Environment 1988 London Graham and Trotman 9. Brown KM Nature and hydrogeologic significance of mud diapirs and diatremes for accretionary systems J. Geophys. Res. Solid Earth 1990 95 B6 8969 8982 10.1029/JB095iB06p08969 10. Roberts HH Paull CK Dillon WP Fluid and gas expulsion on the Northern Gulf of Mexico continental slope: Mud-prone to mineral-prone responses Natural Gas Hydrates: Occurrence, Distribution, and Detection 2001 Washington, D.C. AGU 145 161 11. Kopf AJ Significance of mud volcanism Rev. Geophys. 2002 40 2 1 10.1029/2000RG000093 12. Somoza L Diaz-del-Rio V León R Ivanov M Fernández-Puga M Gardner JM Maestro A Vasquez JT Medialdea T Fernández-Salas LM Seabed morphology and hydrocarbon seepage in the Gulf of Cadiz mud volcano area: acoustic imagery, multibeam and ultra-high resolution seismic data Mar. Geol. 2003 195 153 176 10.1016/S0025-3227(02)00686-2 13. Ceramicola, S., Dupré, S., Somoza, L. & Woodside, J. Cold seep systems. In Submarine Geomorphology 367–387 (Springer, Berlin, 2018). 14. Camerlenghi A Cita MB Della Vedova B Fusi N Mirabile L Pellis G Geophysical evidence of mud diapirism on the Mediterranean Ridge accretionary complex Mar. Geophys. Res. 1995 17 115 141 10.1007/BF01203423 15. Judd AG The global importance and context of methane escape from the seabed Geo. Mar. Lett. 2003 23 147 154 10.1007/s00367-003-0136-z 16. Judd A Hovland M Seabed Fluid Flow: The Impact on Geology, Biology and the Marine Environment 2009 Cambridge Cambridge University Press 17. Sloan ED Jr Clathrate Hydrates of Natural Gases. Revised and Expanded 1998 Boca Raton CRC Press 18. Ginsburg, G.D. & Soloviev, V.A. Submarine gas hydrates. In VNIIOkeangeologian 215 (Saint Petersburg, Russia, 1998). 19. Fehn U Snyder G Egeberg PK Dating of pore waters with 129I: relevance for the origin of marine gas hydrates Science 2000 289 2332 2335 10.1126/science.289.5488.2332 11009415 20. Dickens G On the fate of past gas: What happens to methane released from a bacterially mediated gas hydrate capacitor? Geochem. Geophys. Geosyst. 2001 2 1 10.1029/2000GC000131 21. Bangs NL Sawyer DS Golovchenko X Free gas at the base of the gas hydrate zone in the vicinity of the Chile triple junction Geology 1993 21 905 908 10.1130/0091-7613(1993)021<0905:FGATBO>2.3.CO;2 22. Froelich, P. N., Kvenvolden, K. A., Torres, M. E., Waseda, A., Didyk, B. M. & Lorenson, T. D. Geochemical evidence for gas hydrate in sediment near the Chile Triple Junction. In Proceedings of the Ocean Drilling Program, Scientific Results 141 (Texas, 1995). 23. Morales G Methane hydrates in the Chilean continental margin Electron. J. Biotechnol. 2003 6 80 84 10.4067/S0717-34582003000200002 24. Grevemeyer I Diaz-Naveas JL Ranero CR Villinger HW Heat flow over the descending Nazca plate in central Chile, 32 S to 41 S: observations from ODP Leg 202 and the occurrence of natural gas hydrates Earth Planet. Sci. Lett. 2003 213 285 298 10.1016/S0012-821X(03)00303-0 25. Vargas-Cordero I Tinivella U Accaino F Loreto MF Fanucci F Thermal state and concentration of gas hydrate and free gas of Coyhaique, Chilean Margin (44° 30′ S) Mar. Petrol. Geol. 2010 27 1148 1156 10.1016/j.marpetgeo.2010.02.011 26. Vargas-Cordero I Tinivella U Accaino F Loreto MF Fanucci F Reichert C Analyses of bottom simulating reflections offshore Arauco and Coyhaique (Chile) Geo-Mar. Lett. 2010 30 271 281 10.1007/s00367-009-0171-5 27. Vargas-Cordero I Tinivella U Villar Muñoz L Giustiniani M Gas hydrate and free gas estimation from seismic analysis offshore Chiloé island (Chile) Andean Geol. 2016 43 263 274 10.5027/andgeoV43n3-a02 28. Vargas-Cordero I Tinivella U Villar-Muñoz L Gas hydrate and free gas concentrations in two sites inside the Chilean Margin (Itata and Valdivia Offshores) Energies 2017 10 2154 2165 10.3390/en1012215 29. Villar-Muñoz L Behrmann JH Diaz-Naveas J Klaeschen D Karstens J Heat flow in the southern Chile forearc controlled by large-scale tectonic processes Geo-Mar. Lett. 2014 34 185 198 10.1007/s00367-013-0353-z 30. Villar-Muñoz L Bento JP Klaeschen D Tinivella U Vargas-Cordero I Behrmann JH A first estimation of gas hydrates offshore Patagonia (Chile) Mar. Pet. Geol. 2018 96 232 239 10.1016/j.marpetgeo.2018.06.002 31. Villar-Muñoz L Vargas-Cordero I Bento JP Tinivella U Fernandoy F Giustiniani M Behrmann JH Calderon-Diaz S Gas hydrate estimate in an area of deformation and high heat flow at the Chile triple junction Geosciences 2019 9 28 10.3390/geosciences9010028 32. Yin P Berne S Vagner P Loubrieu B Liu Z Mud volcanoes at the shelf margin of the East China Sea Mar. Geol. 2003 194 135 149 10.1016/S0025-3227(02)00678-3 33. Thatcher KE Westbrook GK Sarkar S Minshull TA Methane release from warming-induced hydrate dissociation in the West Svalbard continental margin: timing, rates, and geological controls J. Geophys. Res. Solid Earth 2013 118 22 38 10.1029/2012JB009605 34. IPCC. Contribution of Working Groups I, II and III to the Fifth Assessment Report of the Intergovernmental Panel on Climate Change. In Climate Change 2014: Synthesis Report; (eds. Pachauri, R.K., Meyer, L.A.) 151 (IPCC: Geneva, Switzerland, 2014). 35. Sibson RH Interactions between temperature and pore-fluid pressure during earthquake faulting and a mechanism for partial or total stress relief Nat. Phys. Sci. 1973 243 66 68 10.1038/physci243066a0 36. Rathburn AE Pérez ME Martin JB Day SA Mahn C Gieskes J Bahls A Relationships between the distribution and stable isotopic composition of living benthic foraminifera and cold methane seep biogeochemistry in Monterey Bay, California Geochem. Geophys. Geosyst. 2003 4 12 10.1029/2003gc000595 37. Ruppel CD Kessler JD The interaction of climate change and methane hydrates Rev. Geophys. 2017 55 126 168 10.1002/2016RG000534 38. Sen Gupta BK Aharon P Benthic foraminifera of bathyal hydrocarbon vents of the Gulf of Mexico: Initial report on communities and stable isotopes Geo-Mar. Lett. 1994 4 88 96 10.1007/BF01203719 39. Reed DW Fujita Y Delwiche ME Blackwelder DB Sheridan PP Uchida T Colwell FS Microbial communities from methane hydrate-bearing deep marine sediments in a forearc basin Appl. Environ. Microbiol. 2002 68 3759 3770 10.1128/AEM.68.8.3759-3770.2002 12147470 40. Torres ME Mix AC Kinports K Haley B Klinkhammer GP McManus J de Angelis MA Is methane venting at the seafloor recorded by δ13 C of benthic foraminifera shells? Paleoceanography 2003 18 1062 10.1029/2002PA000824 41. Chen Y Matsumoto R Paull CK Ussler W III Lorenson T Hart P Winters W Methane-derived authigenic carbonates from the northern Gulf of Mexico- MD02 Cruise J. Geochem. Explor. 2007 95 1 15 10.1016/j.gexplo.2007.05.011 42. Karstens J Haflidason H Becker LW Berndt C Rüpke L Planke S Mienert J Glacigenic sedimentation pulses triggered post-glacial gas hydrate dissociation Nat. Commun. 2018 9 635 10.1038/s41467-018-03043-z 29434195 43. Bernhard JM Buck KR Farmer MA Bowser SS The Santa Barbara Basin is a symbiosis oasis Nature 2000 403 77 80 10.1038/47476 10638755 44. Rathburn AE Levin LA Held Z Lohmann K Benthic foraminifera associated with cold methane seeps on the northern California margin: Ecology and stable isotopic composition Mar. Micropaleontol. 2000 38 247 266 10.1016/S0377-8398(00)00005-0 45. Hill TM Kennett JP Spero HJ Foraminifera as indicators of methane-rich environments: a study of modern methane seeps in Santa Barbara Channel, California Mar. Micropaleontol. 2003 49 123 138 10.1016/S0377-8398(03)00032-X 46. Tomaru H Torres Marta E Matsumoto R Borowski WS Effect of massive gas hydrate formation on the water isotopic fractionation of the gas hydrate system at Hydrate Ridge, Cascadia margin, offshore Oregon Geochem. Geophys. Geosyst. 2006 10.1029/2005GC001207 47. Hesse R Pore water anomalies of submarine gas-hydrate zones as tool to assess hydrate abundance and distribution in the subsurface: What have we learned in the past decade? Earth Sci. Rev. 2003 61 149 179 10.1016/S0012-8252(02)00117-4 48. Kvenvolden, K.A. & Kastner, M. Gas hydrates of the peruvian outer continental margin. Proce. Ocean Drill. Prog. Sci. Res.112 (1990). 49. Sager WW MacDonald IR Hou R Geophysical signatures of mud mounds at hydrocarbon seeps on the Louisiana continental slope, northern Gulf of Mexico Mar. Geol. 2003 198 97 132 10.1016/S0025-3227(03)00097-5 50. Tinivella U Accaino F Della Vedova B Gas hydrates and active mud volcanism on the South Shetland continental margin, Antarctic Peninsula Geo-Mar Lett. 2007 28 2 51. Melnick D Echtler H Inversion of forearc basins in southcentral Chile caused by rapid glacial age trench fill Geology 2006 34 709 712 10.1130/G22440.1 52. Melnick D Bookhagen B Strecker MR Echtler HP Segmentation of megathrust rupture zones from fore-arc deformation patterns over hundreds to millions of years, Arauco peninsula Chile. J. Geophys. Res. Solid Earth 2009 114 1 10.1029/2008JB005788 53. Bangs NL Cande SC Episodic development of a convergent margin inferred from structures and processes along the southern Chile margin Tectonics 1997 16 489 503 10.1029/97TC00494 54. Lohrmann, J. Identification of parameters controlling the accretive and tectonically erosive mass: transfer mode at the South-Central and North Chilean Forearc Using Scaled 2D Sandbox Experiments. PhD Thesis (2002). 55. Mordolovich, C. Geology of a part of the pacific Margin of Chile. In The Geology of Continental Margins 591–598 (1974). 56. González E Hydrocarbon Resources in the Coastal Zone of Chile, Geology of the Andes and its relation to hydrocarbon and mineral resources Earth Sci. 545 Ser. 1989 1989 11 383 404 57. Caress, D. W. & Chayes, D. N. MB-system (versión 5.5.2298) Mapping the 474 seafloor. http://www.ldeo.columbia.edu/MBSystem/html/mbsystem_home.html (2017). 58. Wessel P Smith W Scharroo R Luis JF Wobbe F Generic Mapping Tools: Improved version released EOS Trans. AGU 2013 94 409 410 10.1002/2013EO450001 59. Ondréas H Cannat M Fouquet Y Normand A Sarradin PM Sarrazin J Recent volcanic events and the distribution of hydrothermal venting at the Lucky Strike hydrothermal field, Mid-Atlantic Ridge Geochem. Geophys. Geosyst. 2009 10 2 10.1029/2008GC002171 60. Folk RL Ward WC A study in the significance of grain-size parameters J. Sediment. Petrol. 1957 27 3 26 10.1306/74D70646-2B21-11D7-8648000102C1865D 61. Carver RE Procedures in Sedimentary Petrology 1971 Hoboken Wiley 62. Scasso RA Limarino CO Petrología y diagénesis de rocas clásticas. Asociación Argentina de Sedimentología Publicación Especial 1997 1 257 63. Lozano CF Oceanografía, biología marina y pesca 1978 España Editorial Paraninfo S.A 64. Salamanca M Jara B Distribución y acumulación de plomo (Pb y 210Pb) en sedimentos de los fiordos de la XI región, Chile Revista ciencia y tecnología del mar 2003 26 61 71 65. Byers S Mills E Stewart P A comparison of methods of determining organic carbon in marine sediments, with suggestions for a standard method Hydrobiology 1978 58 43 47 10.1007/BF00018894 66. Luczak C Janquin M Kupka A Simple standard procedure for the routine determination of organic matter in marine sediment Hydrobiology 1997 345 87 94 10.1023/A:1002902626798 67. Holbourn A Andrew SH MacLeod N Atlas of benthic foraminifera 2013 Hoboken Wiley 68. Figueroa S Marchant M Giglio S Ramírez M Foraminíferos bentónicos rotalinidos del centro sur de Chile (36°S-44°S) Gayana (Concepción) 2005 69 329 363 10.4067/S0717-65382005000200013 69. Coplen T Wassenaar L LIMS for Lasers 2015 for achieving long-term accuracy and precision of δ2H, δ17O, and δ18O of waters using laser absorption spectrometry Rapid Commun. Mass Spectrom. 2015 29 2122 2130 10.1002/rcm.7372 26467224 70. Hesse R Harrison WE Gas hydrates (clathrates) causing pore-water freshening and oxygen isotope fractionation in deep-water sedimentary sections of terrigenous continental margins Earth Planet. Sci. Lett. 1981 55 453 462 10.1016/0012-821X(81)90172-2 71. Schmidt, G.A., Bigg, G.R. & Rohling, E.J. Global Seawater Oxygen-18 759 Database. https://data.giss.nasa.gov/o18data/ (1999). 72. LeGrande AN Schmidt GA Global gridded data set of the oxygen isotopic composition in seawater Geophys. Res. Lett. 2006 33 L12604 10.1029/2006GL026011 73. Silva N Rojas N Fedele A Water masses in the Humboldt Current System: Properties, distribution, and the nitrate deficit as a chemical water mass tracer for Equatorial Subsurface Water off Chile Deep Sea Res. II Top. Stud. Oceanogr. 2009 56 1004 1020 10.1016/j.dsr2.2008.12.013 74. Sharp Z Principles of Stable Isotope Geochemistry 2007 New York Prentice Hall 360 75. Rivero-Cortés AI Génesis y procesos asociados a anomalía batimétrica positiva en la zona de offshore de Lebu (37° 36’ 20” S y 73° 44’ 30” W) 2018 Chile Undergraduate Project, Universidad Andrés Bello, Chile 76. Matsumoto, R. & Borowsky, W.S. Gas hydrate estimates from newly determined oxygen isotopic fractionation (αGH–IW) and δ18O anomalies of the interstitial waters: LEG 164, Blake Ridge. In Proceedings of the Ocean Drilling Program, Scientific Results Vol. 164 (ed. Paull, C.K., Matsumoto, R., Wallace, P.J., and Dillon, W.P.) (2000). 77. Nakayama T Tomura S Ozaki M Ohmura R Mori YH Engineering investigation of hydrogen storage in the form of clathrate hydrates: conceptual design of hydrate production plants Energy Fuels 2010 24 2576 2588 10.1021/ef100039a 78. Tréhu AM Flemings PB Bangs N Chevallier J Gracia E Johnson J Liu C-S Liu X Riedel M Torres ME Feeding methane vents and gas hydrate deposits at south Hydrate Ridge Geophys. Res. Lett. 2004 31 L23310 10.1029/2004GL021286 79. Jones RW Foraminifera and their applications J. Geophys. Res. Oceans 2014 118 4811 4821 80. Panieri G Foraminiferal response to an active methane seep environment: a case study from the Adriatic Sea Mar. Micropaleontol. 2006 61 116 130 10.1016/j.marmicro.2006.05.008 81. Dessandier P-A Borrelli C Kalenitchenko D Panieri G Benthic foraminifera in arctic methane hydrate bearing sediments Front. Mar. Sci. 2019 6 765 10.3389/fmars.2019.00765 82. Bernhard JM Buck KR Barry JP Monterey Bay cold-seep biota: assemblages, abundance, and ultrastructure of living foraminifera Deep Sea Res. 2001 I 48 2233 2249 10.1016/S0967-0637(01)00017-6 83. Wollenburg JE Kuhnt W The response of benthic foraminifers to carbon flux and primary production in the Arctic Ocean Mar. Micropaleontol. 2000 40 3 189 231 10.1016/S0377-8398(00)00039-6 84. Mordojovich C Sedimentary basins of the Chilean Pacific offshore Energy. Resources of the Pacific Region Am. Assoc. Petrol. Geol 1981 2 63 82 85. Oshima M Taniguchi Y Akasaka M Abe K Ichikawa H Tabei Y Tanaka J Development of a visible marker trait based on leaf sheath-specific anthocyanin pigmentation applicable to various genotypes in rice Breed Sci. 2019 69 244 254 10.1270/jsbbs.18151 31481833 86. Dyer KR Sediment transport processes in Estuaries Dev. Sedimentol. 1995 53 423 449 10.1016/S0070-4571(05)80034-2 87. Pineda V Granulometría y geoquímica de los sedimentos marinos en el área comprendida entre el seno de Reloncaví y golfo Corcovado, Chile, Crucero CIMAR 10 fiordos Revista Ciencia y Tecnología del Mar 2009 32 1 27 47 88. Hyodo M Yoneda J Yoshimoto N Nakata Y Mechanical and dissociation properties of methane hydrate-bearing sand in deep seabed Soils Found. 2013 53 299 314 10.1016/j.sandf.2013.02.010 89. White WM Geochemistry 2013 Hoboken Wiley 90. Wallmann K Pinero E Burwicz E Haeckel M Hensen C Dale A Ruepke L The global inventory of methane hydrate in marine sediments: a theoretical approach Energies 2012 5 2449 2498 10.3390/en5072449 91. Geersen J Behrmann J Völker D Krastel S Ranero C Diaz-Naveas J Weinrebe W Active tectonics of the South Chilean marine forearc (35°S–40°S) Tectonics 2011 30 1 16 10.1029/2010TC002777 92. Becerra J Contreras-Reyes E Arriagada C Seismic structure and tectonics of the southern Arauco Basin, south-central Chile (~38S) Tectonophysics 2013 592 53 66 10.1016/j.tecto.2013.02.012 93. Vargas-Cordero I Tinivella U Accaino F Basal and frontal accretion processes versus BSR characteristics along the Chilean Margin J. Geol. Res. 2011 10.1155/2011/846101 94. Nelson AR Manley WF Holocene coseismic and aseismic uplift of Isla Mocha, south-central Chile Quatern. Int. 1992 15 61 76 10.1016/1040-6182(92)90036-2 95. Lohrmann, J., Kukowski, N., Krawczyk, C.M., Oncken, O., Sick, C., Sobiesiak, M. & Rietbro, A. Subduction channel evolution in brittle fore-arc wedges - a combined study with scaled sandbox experiments, seismological and reflectionseismic data and geological field evidence. In The Andes 237–262 (Springer, Berlin, 2006). 96. Tréhu AM Torres ME Moore GF Suess E Bohrmann G Temporal and spatial evolution of a gas hydrate–bearing accretionary ridge on the Oregon continental margin Geology 1999 27 939 942 10.1130/0091-7613(1999)027<0939:TASEOA>2.3.CO;2 97. Pecher IA Henrys SA Ellis S Chiswell SM Kukowski N Erosion of the seafloor at the top of the gas hydrate stability zone on the Hikurangi Margin, New Zealand Geophys. Res. Lett. 2005 32 24 10.1029/2005GL024687 98. Maclennan J Jones SM Regional uplift, gas hydrate dissociation and the origins of the Paleocene-Eocene Thermal Maximum Earth Planet. Sci. Lett. 2006 245 65 80 10.1016/j.epsl.2006.01.069 99. Wallmann K Riedel M Hong WL Patton H Hubbard A Pape T Hsu CW Schmidt C Johnson JE Torres ME Andreassen K Berndt C Bohrmann G Gas hydrate dissociation off Svalbard induced by isostatic rebound rather than global warming Nat. Commun. 2018 9 1 83 10.1038/s41467-017-02550-9 29311564 100. Rabassa J Clapperton CM Quaternary glaciations of the southern Andes Quatern. Sci. Rev. 1990 9 153 174 10.1016/0277-3791(90)90016-4 101. Glasser NF Jansson KN Harrison S Kleman J The glacial geomorphology and Pleistocene history of South America between 38°S and 56°S Quatern. Sci. Rev. 2008 27 365 390 10.1016/j.quascirev.2007.11.011 102. Hein AS Hulton NR Dunai TJ Sugden DE Kaplan MR Xu S The chronology of the Last Glacial Maximum and deglacial events in central Argentine Patagonia Quatern. Sci. Rev. 2010 29 1212 1227 10.1016/j.quascirev.2010.01.020 103. Geersen J Völker D Behrmann JH Reichert C Krastel S Pleistocene giant slope failures offshore Arauco peninsula, southern Chile J. Geol. Soc. 2011 168 1237 1248 10.1144/0016-76492011-027 104. Vargas-Cordero I Tinivella U Villar-Muñoz L Bento JP High Gas hydrate and free gas concentrations: an explanation for seeps offshore South Mocha Island Energies 2018 11 3062 10.3390/en11113062 105. Hormazabal S Combes V Morales CE Correa-Ramirez MA Di Lorenzo E Nuñez S Intrathermocline eddies in the coastal transition zone off central Chile (31–41°S) J. Geophys. Res. Oceans 2013 118 4811 4821 10.1002/jgrc.20337 106. Silva N Neshyba S On the southernmost extension of the Peru-Chile undercurrent Deep Sea Res. 1979 26 1387 1393 10.1016/0198-0149(79)90006-2 107. Combes V Hormazabal S Di Lorenzo E Interannual variability of the subsurface eddy field in the Southeast Pacific J. Geophys. Res. Oceans 2015 120 4907 4924 10.1002/2014JC010265 108. Shaffer G Hormazabal S Pizarro O Salinas S Seasonal and interannual variability of currents and temperature off central Chile J. Geophys. Res. Oceans 1999 104 C12 29951 29961 10.1029/1999JC900253 109. Hormazabal S Shaffer G Letelier J Ulloa O Local and remote forcing of sea surface temperature in the coastal upwelling system off Chile J. Geophys. Res. Oceans 2001 106 16657 16671 10.1029/2001JC900008 110. Hormazabal S Shaffer G Pizarro O Tropical Pacific control of intraseasonal oscillations off Chile by way of oceanic and atmospheric pathways Geophys. Res. Lett. 2002 29 6 10.1029/2001GL013481 111. Sobarzo M Shearman RK Lentz S Near-inertial motions over the continental shelf off Concepción, central Chile Prog. Oceanogr. 2007 75 348 362 10.1016/j.csr.2010.06.003 112. Schaffer G Leth O Ulloa O Bendtsen J Daneri G Dellarossa V Warming and circulation change in the eastern South Pacific Ocean Geophys. Res. Lett. 2000 27 1247 1250 10.1029/1999GL010952 113. Alessandrini G Tinivella U Giustiniani M Vargas-Cordero I Castellaro S Potential instability of gas hydrates along the Chilean margin due to ocean warming Geosciences 2019 9 234 10.3390/geosciences9050234 114. Mix, A.C., Tiedemann, R. & Blum, P. Proc. ODP, Init. Repts. 202: College 676 Station, TX (Ocean Drilling Program), Doi:10.2973/odp.proc.ir.202 (2003). 115. Tinivella U Giustiniani M Variations in BSR depth due to gas hydrate stability versus pore pressure Glob. Planet. Change 2013 100 119 128 10.1016/j.gloplacha.2012.10.012 116. Dickens G Quinby-Hunt MS Methane hydrate stability in seawater Geophys. Res. Lett. 1994 21 2115 2118 10.1029/94GL01858 117. Lamy, F. & Kaiser, J. Glacial to Holocene paleoceanographic and continental paleoclimate reconstructions based on ODP Site 1233/GeoB 3313 off southern Chile. In Past Climate Variability in South America and Surrounding Regions 129–156 (2009). 118. Westbrook GK Thatcher KE Rohling EJ Piotrowski AM Pälike H Osborne AH Hühnerbach V Escape of methane gas from the seabed along the West Spitsbergen continental margin Geophys. Res. Lett. 2009 36 15 10.1029/2009GL039191 119. Kim JH Schneider RR Hebbeln D Müller PJ Wefer G Last deglacial sea-surface temperature evolution in the Southeast Pacific compared to climate changes on the South American continent Quatern. Sci. Rev. 2002 21 2085 2097 10.1016/S0277-3791(02)00012-4 120. Reagan MT Moridis GJ Oceanic gas hydrate instability and dissociation under climate change scenarios Geophys. Res. Lett. 2007 34 L22709 10.1029/2007GL031671 121. Geersen J Scholz F Linke P Schmidt M Lange D Behrmann JH Volker D Hensen C Fault zone controlled seafloor methane seepage in the rupture area of the 2010 Maule earthquake, Central Chile Geochem. Geophys. Geosyst. 2016 17 4802 4813 10.1002/2016GC006498