Massive Natural Gas Hydrate Dissociation During the Penultimate Deglaciation (~130 ka) in the South China Sea

Wei, Jiangong; Wu, Tingting; Miao, Xiaoming; Su, Pibo

doi:10.3389/fmars.2022.875374

ORIGINAL RESEARCH article

Front. Mar. Sci., 06 May 2022

Sec. Marine Biogeochemistry

Volume 9 - 2022 | https://doi.org/10.3389/fmars.2022.875374

This article is part of the Research Topic The effects of early diagenesis in various marine environments on the stable isotope records of environmental conditions and biogeochemical processes View all 7 articles

Massive Natural Gas Hydrate Dissociation During the Penultimate Deglaciation (~130 ka) in the South China Sea

Jiangong Wei^1,2,3†

Tingting Wu^1,2†

Xiaoming Miao^4,5*

Pibo Su^1,2,3,6*

¹Southern Marine Science and Engineering Guangdong Laboratory (Guangzhou), Guangzhou, China
²MLR Key Laboratory of Marine Mineral Resources, Guangzhou Marine Geological Survey, Guangzhou, China
³Academy of South China Sea Geological Science, China Geological Survey, Sanya, China
⁴College of Marine Geosciences, Ocean University of China, Qingdao, China
⁵Key Laboratory of Submarine Geosciences and Prospecting MOE China, Ocean University of China, Qingdao, China
⁶Hubei Key Laboratory of Marine Geological Resources, China University of Geosciences, Wuhan, China

There is widespread and growing scientific interest in the impact of massive gas hydrate dissociation on the global environment and climate in geological history. Based on that a global negative excursion of carbon isotopic compositions in marine (foraminifera) and terrestrial (the organic carbon and calcite) sedimentary records occurred during the penultimate deglaciation (~130 ka), we believe that methane released by hydrate dissociation may play a role in accelerating the initial increase of atmospheric methane. In order to prove that massive natural gas hydrate dissociation occurred in this period, we aim to seek for evidence of gas hydrate dissociation from seep carbonate. Here, X-ray diffraction, carbon and oxygen isotopic compositions, trace elements, and U-Th dating analyses were conducted on the deeply-buried authigenic carbonate obtained by drilling in the northern continental slope of the South China Sea. Authigenic carbonate formed at ~130 ka showed obvious characteristics of negative excursion of carbon isotope, positive excursion of oxygen isotope, and enrichment of redox sensitive elements such as Mo, U and As. These results, in particular the high oxygen isotopic compositions of carbonate, point to massive gas hydrate dissociation in the northern continental slope of the South China Sea during the penultimate deglaciation. It is further speculated that massive gas hydrate dissociation might have also occurred on a global scale, contributing to the increase of atmospheric carbon dioxide and methane concentrations during the penultimate deglaciation, and may eventually cause global carbon isotope negative excursion.

Introduction

There are huge reserves (1,000–10,000 Gt) of natural gas hydrate widely distributed in permafrost and continental margin sediments (Dickens et al., 1997). Changes in the sedimentary environment, such as temperature and sea level changes, can lead to the dissociation of natural gas hydrate and the release of methane (Crémière et al., 2016; Argentino et al., 2019; Chen et al., 2019; Deng et al., 2020). The release of these gasses can have an impact on the ocean, atmosphere, and even the global environment. Several major events that occurred in geological history may have been related to the massive gas hydrate dissociation, including the “Snowball Earth” termination event (Kennedy et al., 2008), the Permo-Triassic boundary (P/T) (Sluijs et al., 2007), the Early Toarcian oceanic anoxic event (OAE) during the Jurassic period (Hesselbo et al., 2000), the early Cretaceous (Jahren et al., 2001), the latest Palaeocene Thermal Maximum (LPTM) (Dickens et al., 1995; Dickens, 2001), and Quaternary Interstadials (Kennett et al., 2000; Dean et al., 2015).

Many marine and terrestrial sedimentary records of carbon isotopes indicate obvious negative excursion during the penultimate deglaciation (~130 ka). These records have been collected from regions such as the East Pacific (Shackleton and Hall, 1989), the Western Pacific (Schmidt et al., 1993), the northern Atlantic (Oppo et al., 1997), the South China Sea (Li and Wang, 2006), the Indian Ocean (Farrell and Janecek, 1991), Lake Baikal (Chappellaz et al., 1990), North America (Coplen et al., 1994), and Europe (Frogley et al., 1999). In addition, an ice core taken from Vostok recorded a rapid increase in global atmospheric methane and carbon dioxide levels during that time (Chappellaz et al., 1990), indicating the release of a large quantity of methane into the atmosphere during this period. Many people believe that the methane might originate from the increased vegetation (Coplen et al., 1994; Li and Wang, 2006; Häuselmann et al., 2015) or the extended wetlands, peatlands, and inundated floodplains (Chappellaz et al., 1990; Lourantou et al., 2010; Bock et al., 2017; Schmidely et al., 2021) in deglaciation times. However, most studies in recent years have shown that the release of methane from gas hydrates has a significant impact on climate change (Kennett et al., 2000; Shakhova et al., 2010; Berndt et al., 2014), including in polar regions (Serov et al., 2017) and other deep-sea areas (Reagan and Moridis, 2007; Foschi et al., 2020). Therefore, from the perspective of methane seepage, we believe that a large amount of methane from hydrate dissociation contributes to climate change. Some obvious directions of related future research include discussing the driver of negative excursion of global carbon isotopes, and identifying the link to the massive global dissociation of natural gas hydrate.

Previous studies have shown that the dissociation of natural gas hydrate will increase methane flux and further accelerate anaerobic oxidation of methane, thereby forming a large quantity of seep authigenic carbonate (Crémière et al., 2016; Feng et al., 2016; Lu et al., 2018; Argentino et al., 2019; Argentino et al., 2020; Lu et al., 2021). Sea-floor methane seepage leaves behind carbonate that have distinct geochemical signals that can be attributed to their origin (Svensen, 2012; Hu et al., 2014; Smrzka et al., 2020; Lin et al., 2021). Therefore, the hypothesis of massive gas hydrate dissociation during ~130 ka can be tested by the study of cold-seep authigenic carbonate (Svensen, 2012). The South China Sea is considered a natural laboratory for the study of the dissociation and evolution of natural gas hydrate due to the extensive development of this resource in this region (Han et al., 2014; Feng et al., 2016; Wei et al., 2019; Ye et al., 2019; Miao et al., 2021a; Miao et al., 2021b; Miao et al., 2022). Here, the present study conducted X-ray diffraction, carbon and oxygen isotopic compositions, trace element, and U-Th dating analyses of authigenic carbonate obtained by drilling in the Qiongdongnan Basin of South China Sea (Figure 1). The aim of the present study was to identify regional evidence of massive gas hydrate dissociation during ~130 ka.

FIGURE 1

Figure 1 (A) Map of earth showing regions containing marine and terrestrial sedimentary records of carbon isotopes indicating negative excursion during the penultimate deglaciation (~130 ka). The red rectangle represents the study area of the present study in the South China Sea; (1) represents the Vostok ice core (Chappellaz et al., 1990); (2) represents the South China Sea (Li and Wang, 2006); (3) represents the Indian Ocean (Farrell and Janecek, 1991); (4) represents the Western Pacific (Shackleton and Hall, 1989); (5) represents the Eastern Pacific (Schmidt et al., 1993); (6) represents the North Atlantic (Oppo et al., 1997); (7) represents North America (Coplen et al., 1994); (8) represents Lake Baikal (Prokopenko and Williams, 2004). (B) Map showing the study area of the present study located in the South China Sea.

Materials and Methods

Materials

The GMGS5-W08 site was drilled to a depth of ~200 m below the seafloor (mbsf) at a water depth of ~1735 m in 2018 (Wei et al., 2019; Wei et al., 2020). And the bottom water temperature is ~ 3.5°C (Wei et al., 2019). A large number of gas hydrates were found in 54 mbsf, 63 mbsf, 64 mbsf and 69 mbsf and 70 mbsf (Wei et al., 2019). And continuous seep carbonate was discovered at the interval between 52–54 mbsf (Wei et al., 2020).

Methods

X-ray diffraction was used to analyze and identify the carbonate mineral composition. The samples were first dried for 7 h and then gently ground by a mortar. The ground sample was packed into X-ray sample holders for analysis. Scans were run at room temperature using a Rigaku SmartLab-9kW X-ray diffractometer with 8 kW power. The species and content of the minerals were analyzed by using the software PDXL2. The 2θ range is 3° to 75°, and speed of measurement is 7°/min. The voltage and current we used were 40 kV and 200 mA, respectively. In addition, the type of radiation we used was Cu.

Stable carbon and oxygen isotopic compositions were measured using a Thermo MAT-253 isotope ratio mass spectrometer. CO₂ gas was extracted by a reaction with supersaturated phosphoric acid on a Thermo Kiel IV Carbonate Device and was introduced into the MAT-253 dual inlet system. The isotope ratios were reported relative to Vienna Peedee Belemnite (VPDB). The precisions of δ¹³C and δ¹⁸O were ±0.1‰ and ±0.25‰, respectively.

For major and trace element content analysis, bulk seep carbonates were completely dissolved by HF and HNO₃ solutions following the procedure described in Hu et al. (2014). Briefly, ~50 mg sample was weighed and transferred into a pre-cleaned Teflon beaker followed by the addition of ultra-pure 1 ml of HF and 1 ml of HNO₃ solution. The beakers were then placed in steel cans and subjected to high temperature (185°C) and high pressure. After 36 h, the solution was dried on a hotplate. The residues were fully digested using a mixture of concentrated 2 ml HNO₃ and 3 ml Milli-Q water. Thereafter, the beakers were placed into steel cans at 120°C for 5 h. After cooling, the solution was diluted to 20 ml with Milli-Q water. Major elements were analysed with ICP-OES (Optima 8300, PerkinElmer, MA, USA), and trace elements were measured via ICP-MS (X Series2, Thermo Fisher Scientific, MA, USA). The analytical precision was better than 5%.

The ²³⁰Th dating work was performed at the Isotope Laboratory, Xi’an Jaiotong University using multi-collector inductively coupled plasma mass spectrometers (MC-ICP-MS) (Thermo-Finnigan Neptune-plus). We used standard chemistry procedures to separate U and Th for dating (Edwards et al., 1987). A triple-spike (²²⁹Th–²³³U–²³⁶U) isotope dilution method was employed to correct for instrumental fractionation and determine U-Th isotopic ratios and concentrations. The instrumentation, standardization and half-lives are reported in refs. Cheng et al. (2000, 2013). All U-Th isotopes were measured on a MasCom multiplier behind the retarding potential quadrupole in the peak-jumping mode. We followed similar procedures of characterizing the multiplier as described in ref. Cheng et al. (2000). Uncertainties in U-Th isotopic data were calculated offline at 2σ level, including corrections for blanks, multiplier dark noise, abundance sensitivity, and contents of the same nuclides in spike solution. Corrected ²³⁰Th ages assume the initial ²³⁰Th/²³²Th atomic ratio of 4.4 ± 2.2 x10^-6, the values for a material at secular equilibrium with the bulk earth ²³²Th/²³⁸U value of 3.8.

Results

The carbonate content of GMGS5-W08 varied between 82.2%–100% (mean of 92.3%). Carbonates in the samples were mainly composed of aragonite (> 80 wt %) with minor amounts of calcite (Table 1).

TABLE 1

Table 1 Mineralogical compositions of seep carbonates.

The δ¹³C values in the carbonate layers of GMGS5-W08 varied from −38.1‰ to −15.2‰, with mean values of −32.9‰ (n = 10) (Figure 2 and Table 2). The δ¹⁸O values in the carbonate layers of GMGS5-W08 ranged from 3.8‰ to 5.7‰, with mean values of 4.8% (n = 10) (Figure 2 and Table 2).

FIGURE 2

Figure 2 Plots of carbon vs. oxygen stable isotope values. The data of GMGS2-08 from Chen et al., 2019.

TABLE 2

Table 2 Carbon and oxygen isotopic ratios of seep carbonate layers.

Contents of major elements and trace elements in the bulk seep carbonates are presented in Table 3. In GMGS5-W08, Al content ranged from 0.09 to 3.97 wt. % (average 2.15 wt. %, n=16). U content ranged from 4.27 to 25.50 μg/g (average 16.99 μg/g, n=16). Th content ranged from 0.13 to 3.73 μg/g (average 2.66 μg/g, n=16). Mo content ranged from 4.06 to 52.74 μg/g (average 25.82 μg/g, n=16). V content ranged from 8.96 to 54.40 μg/g (average 36.20 μg/g, n=16). Ni content ranged from 3.66 to 83.38 μg/g (average 21.96 μg/g, n=16).

TABLE 3

Table 3 Major and trace element contents of seep carbonates.

The enrichment factor (EF) was calculated as X_EF = [(X/Al)_sample/(X/Al)_PAAS], where X and Al represent the weight concentrations of elements X and Al, respectively. The samples were normalized using the Post Archean Australian Shale (PAAS) composition (Taylor and McLennan, 1985). By calculation, the Mo enrichment factors (EF) in GMGS5-W08 seep carbonate vary from 56.4 to 2710.0 (average 278.9, n= 16) (Figure 3A). The U enrichment factors in GMGS5-W08 seep carbonate vary from 16.3 to 407.9 (average 53.1, n= 16) (Figure 3A). The overall range of U/Th ratio is 5.08 to 86.87 (average 11.37, n= 16) (Figure 3B). The overall range of V/(V+Ni) ratio is 0.36 to 0.82 (average 0.63, n= 16) (Figure 3C).

FIGURE 3

Figure 3 (A) Plots of Mo_EF vs. U_EF of the seep carbonates from the sites GMGS5-W08 and GMGS2-08 (from Deng et al., 2020) taken from the South China Sea (see Figure 1). Samples were normalized by the standard Post Archean Australian Shale (PAAS) (Taylor and McLennan, 1985). The diagonal lines represent multiples of the Mo:U weight ratio of present-day sea water (SW) and were taken from Algeo and Tribovillard (2009). (B) Plots of U vs. Th contents of the seep carbonates. The diagonal line represents U:Th =1.25 and was taken from Jones and Manning (1994). (C) Plots of V vs. (V+Ni) contents of the seep carbonates. The diagonal line represent V: (V+Ni) = 0.54 and was taken from Jones and Manning (1994).

Table 4 represents the U-Th ages of the three carbonate samples. The result showed that the ages of GMGS5-W08 (52–54mbsf) are 114.6 ± 0.5–136.3 ± 3.6 ka.

TABLE 4

Table 4 U-Th isotopic data and calculated ages of seep carbonates.

Discussion

Fluid Sources and Formation Environments of Authigenic Carbonate

The carbon isotopic composition of authigenic carbonate can reveal the source of carbon during its formation, and the carbon isotopic composition is the most important indicator of methane-derived (Peckmann and Thiel, 2004; Lu et al., 2018; Deng et al., 2021; Lu et al., 2021). In general, seep carbonate shows obvious carbon isotope depletion, which is closely related to the anaerobic oxidation of methane (Peckmann and Thiel, 2004). The carbon isotopic ratios of the authigenic carbonate analyzed in the present study were relatively low, with the values of all samples, except for one, less than − 30‰ (Figure 2 and Table 2). This is very similar to the cold seep carbonates from the site GMGS2-08 formed during ~ 130ka (Chen et al., 2019; Deng et al., 2021). These results indicate that the genesis of authigenic carbonate was related to methane seepage. However, the origin of methane reflected by methane carbon isotopes may vary. Among the samples, the carbon isotopic characteristics of carbonates from the site GMGS5-W08 were more similar to those of seep carbonates in the Gulf of Cadiz (Wang et al., 2015) and the Qiongdongnan Basin (Liang et al., 2017), and mainly of thermogenic or mixed origin. Carbonates in GMGS2-08 showed lower carbon isotopic ratios, indicating an apparent microbial methane component (Chen et al., 2019; Deng et al., 2021). Oxygen isotopes are mainly used to calculate equilibrium temperatures at the time of carbonate formation (Crémière et al., 2016). The present study used the oxygen isotopic fractionation equation of the aragonite water system by Kim et al. (2007) to calculate the expected equilibrium oxygen isotopic composition of aragonite. In this paper, we assume the bottom water temperature is 3.6 °C (Wei et al., 2019) and the δ¹⁸O value of bottom water is 0‰ V-SMOW. The theoretical equilibrium values of oxygen isotopes of aragonite in GMGS5-W08 were 3.1‰ (Figure 2), higher than GMGS2-08(2.4‰, Chen et al., 2019). In addition, the values of δ¹⁸O of GMGS5-W08 exceeded the equilibrium theoretical values, showing characteristics of rich ¹⁸O (Figure 2). This is very similar to other cold seep carbonates (Han et al., 2014; Feng and Chen, 2015; Crémière et al., 2016; Argentino et al., 2019) and is thought to be caused by ¹⁸O-rich fluid of hydrate dissociation. During the period of ~130 ka, the oxygen isotope values recorded by foraminifera in the south China sea decreased significantly due to the increase of temperature (Li and Wang, 2006). Thus, the enrichment of ¹⁸O in carbonate is influenced by fluids other than seawater. Two main sources of ¹⁸O-rich fluids exist in the deep-sea sedimentary environment: (1) dehydration of clay minerals at greater depths (Hesse, 2003) and; (2) dissociation of natural gas hydrate (Chen et al., 2019). However, dehydration of clay minerals is mainly caused by transformations between smectite and illite which require higher diagenetic temperature and pressure, and no significant smectite-illite mixtures have been found in the sediments. Therefore, dehydration of clay minerals may not be the main cause of heavy oxygen enrichment. However, the conditions of the South China Sea are favorable for the formation and development of natural gas hydrate (NGH). In fact, NGH was found in both GMGS5-W08 and GMGS2-08 (Chen et al., 2019; Wei et al., 2019). Therefore, the present study proposed that the generation of ¹⁸O-rich fluid by dissociation of hydrate is the main driver of positive excursion of oxygen isotopes in the carbonate of the study area.

In addition to the stable carbon and oxygen isotopes compositions, the trace element composition of authigenic carbonate can be used to define the sedimentary environment and fluid geochemical characteristics at the time of their formation (Smrzka et al., 2020). A large amount of methane released from the dissociation of hydrate will accelerate the sulfate-driven anaerobic oxidation of methane (SD-AOM) (Peckmann and Thiel, 2004). Consequently, a large amount of H₂S will be released into pore water or even seawater, resulting in an environment in which carbonate deposition enters a reductive environment, thereby resulting in the enrichment of some trace elements (Hu et al., 2014; Smrzka et al., 2020). In our samples, we observed significant enrichment of Mo and U (Figure 3A). Previous studies have shown that the mechanisms responsible for U and Mo enrichment are different under anoxic conditions (Peketi et al., 2012; Sato et al., 2012). U tends to deposit in the Fe reduction zone, whereas Mo concentrates only in environments containing H₂S, and U is usually enriched earlier than Mo (Algeo and Tribovillard, 2009). However, the high methane fluxes will compress the suboxic and sulfidic zones into a narrow zone close to each other, resulting in the co-enrichment of Mo and U (Chen et al., 2016; Miao et al., 2021a; Miao et al., 2022). Therefore, the co-enrichment of Mo and U in methane seepage environment is very common and has been widely used to reconstruct carbonate formation environments (Deng et al., 2020; Smrzka et al., 2020), such as the South China Sea (Lin et al., 2021; Miao et al., 2021a) and northern Apennines (Argentino et al., 2019). As shown in Figure 4A, the carbonates of GMGS5-W08 and GMGS2-08 exhibit similar geochemical characteristics. Mo_EF of carbonates in the study area was significantly greater than U_EF, and the Mo_EF/U_EF ratio of GMGS5-W08 ranges from 2.1 to 7.6. In addition, plotted in the Mo_EF vs. U_EF diagram, they mostly fall in > 0.3×(Mo/U)_SW trend lines, indicating that the carbonates were formed in an anoxic environment or even sulfidic environment (Chen et al., 2016; Argentino et al., 2019). At the same time, the ratios of U/Th (> 1.25) and V/(V+Ni) (> 0.54) of seep carbonates also supports this view (Figures 3B, C) (Jones and Manning, 1994; Wignall and Twitchett, 1996).

FIGURE 4

Figure 4 Arsenic enrichments in the studied carbonate samples. (A), As_EF vs. Mo_EF diagram indicates a moderate correlation (R² = 0.7). (B), As_EF and Fe/Al ratios do not correlate (R² < 0.1). The data of GMGS2-08 from Deng et al., 2020.

In addition, organic matter, iron, and manganese oxides are potential hosts for Mo in the sediments (Algeo and Tribovillard, 2009; Scholz et al., 2011). Mo enrichment is usually evident in organic-rich sediments at modern continental margins (Scholz et al., 2011). However, in this sample, other redox sensitive elements (e.g., Ni, Cu, and Zn) that are associated with organic matter are not enriched (Table 3), which is obviously different from organic-rich sediments (Scholz et al., 2011). At the same time, the extremely low Mn content (average value was 0.059%) and anoxic environment (U/Th > 1.25 and V/(V+Ni) > 0.54) indicate that the content of Mn oxides in carbonate hardly exists. Therefore, Mo enrichment has little relation with organic matter and Fe and Mn oxides (Smrzka et al., 2020). And, considering comprehensively the characteristics of authigenic carbonate, the current study proposed that the formation process of authigenic carbonate in the study area is mainly controlled by the dissociation of hydrate.

Interestingly, arsenic (As) in the the carbonates of GMGS5-W08 and GMGS2-08 (Deng et al., 2020) has obvious enrichment characteristics (Figure 4, 0.1 to 179.8, the mean value was 20.4). This phenomenon is also common in cold-seep environments and is often used to determine the source of Mo (Hu et al., 2014; Argentino et al., 2019; Lin et al., 2021; Miao et al., 2022). The particulate shuttle process is one of the reasons for the common enrichment of Mo and As (Algeo and Tribovillard, 2009; Scott and Lyons, 2012). Iron and manganese (hydrogen) oxides can remove trace elements from water and then transfer to surface sediments (Smrzka et al., 2020). Subsequently, in the sulfidic environment, the adsorbed trace elements are released into the pore water, which is eventually scavenged again by authigenic iron sulfide minerals (Hu et al., 2014; Scholz et al., 2017). In our study, although As_EF and Mo_EF has a very high positive correlation (Figure 4, R² = 0.7), the correlation between the As_EF and the ratio of Fe/Al is very low (Figure 4, R² < 0.1). Obviously, particulate shuttle process is not the main factor controlling Mo and As enrichment. However, the common enrichment of Mo and U in carbonate indicates an overlap or at least a close proximity of the iron reduction zone and the sulfate-methane transition zones (SMTZ). Moreover, combined with bivalve shells in the samples (Wei et al., 2020; Deng et al., 2021), we believe that methane seepage is strong and SMTZ is formed at or near the seafloor. As a result, We suggest that Mo fixation in authigenic iron sulfide can be enhanced by seeps at shallow SMTZs, with seawater being the main source of Mo (Peketi et al., 2012; Lin et al., 2021).

Significance of Cold Seep Activity During the Penultimate Deglaciation (~130 ka)

U-Th dating of cold seep carbonate is of great significance for determining the age of dissociation of natural gas hydrates (Crémière et al., 2016; Chen et al., 2019). The results of U-Th dating in the current study showed that although the seep carbonates in the areas of hydrate occurrence in the South China Sea are buried in different sedimentary horizons, they were formed at roughly the same time, i.e., during the penultimate deglaciation (~130 ka) (Figure 5 and Table 4). According to the seismic reflection profiles of GMGS5-W08, the authigenic carbonate sample at 52-54 mbsf coincides with a high-amplitude reflector at 54 mbsf that extends laterally over 4500 m, which suggests that the carbonate concretions were developed not only within the gas chimney but across the paleo-seafloor (Ye et al., 2019). In addition, the Dongsha area in the South China Sea has also found cold seep carbonates formed at ~130ka (Chen et al., 2019; Deng et al., 2021). This result confirms that a massive gas hydrate dissociation event occurred in the northern continental slope of the South China Sea during this period.

FIGURE 5

Figure 5 Ages of seep carbonates versus carbon isotopic values (Shackleton and Hall, 1989; Chappellaz et al., 1990; Farrell and Janecek, 1991; Schmidt et al., 1993; Coplen et al., 1994; Prokopenko and Williams, 2004; Li and Wang, 2006) from global sites and CH4 and CO2 contents from Vostok (Chappellaz et al., 1990). The sea-level curve is modified from Rohling et al. (2009).

In general, the pressure change caused by sea level fall and temperature change caused by bottom water temperature rise are the main factors that trigger hydrate dissociation in continental margin (Kennett et al., 2000; Shakhova et al., 2010; Dean et al., 2015; Deng et al., 2021). In Figure 5, we find that the carbonates were mainly formed in MIS6/5e. During the MIS6/5e transition period, the sea level was higher than the modern sea level (Rohling et al., 2009), which increased the stability of seabed hydrate and inhibited the dissociation of hydrate. Therefore, sea level change was not the main factor of hydrate dissociation during this period. However, Chen et al. (2019) found that the bottom water temperature of the South China Sea increased by 1.8-4.5°C during this period, which was enough to trigger the dissociation of hydrate. Because every 1°C increase in the temperature of bottom water is enough to trigger the dissociation of local hydrates (Reagan and Moridis, 2007). Therefore, we believe that the dissociation of hydrate in this period was caused by the rise of bottom water temperature. During the transition from glacial to interglacial, the temperature of bottom water increased in almost all sea areas (Rohling et al., 2014). For example, during deglaciation, the temperature of the bottom water in the Atlantic increased by 3- 4.5°C (Dwyer et al., 1995). Therefore, we believe that hydrate dissociation may also occur in other hydrate regions around the world. In Bock et al. (2017), when the concentration of CH₄ starts to increase at the initial stage of the penultimate deglaciation, the δD and δ¹³C values of CH₄ are also increased slightly. Such synchronous variations may indicate methane source from gas hydrate dissociation. This result combined with the observed negative excursion of global carbon isotopes and the rapid increase in CH₄ content in the atmosphere during this period (Figure 5) indicates that the hydrate dissociation event during the penultimate deglaciation period (~130 ka) may have occurred at a global scale. According to “ The Clathrate Gun Hypothesis ” hypothesis (Kennett et al., 2003), we believe that this event may be one of the driver of the global negative excursion of carbon isotopes and the increase in atmospheric methane and carbon dioxide content. Throughout the late Pleistocene, we find that this phenomenon is quite common. Hydrate dissociation events occurred in MIS10/9 (Tong et al., 2013), MIS4/3 (Han et al., 2014; Yang et al., 2018) and MIS2/1 (Wei et al., 2020; Deng et al., 2021) periods.

In addition, age of carbonate (114.6–136.3ka) indicate that methane seepage occurred not only during MIS6/5e, but throughout MIS5e (Figure 5). Methane is a very important greenhouse gas, and its increase in the atmosphere is bound to cause global warming (Hesselbo et al., 2000; Kennett et al., 2000). At the same time, we believe that the climate warming caused by the continuous release of methane may be one of the reasons for maintaining the temperature of interglacial and delaying the arrival of glacial.

Conclusion and Outlook

The geochemical characteristics and U-Th dating of deep authigenic carbonate in different areas of the South China Sea provide strong evidence for a massive gas hydrate dissociation event during the penultimate deglaciation (~130 ka). Although the present study focused on carbonate in the South China Sea, this phenomenon should also exist in the seeping hydrates of other marine areas worldwide, which should be confirmed by future studies. At the same time, the current study proposes that the massive gas hydrate dissociation during this period was related to changes relating to the glacial-interglacial period. The rise of sea water temperature that occurred at the end of the glacial period and beginning of the interglacial period was the main driver of the dissociation of hydrate (Chen et al., 2019; Deng et al., 2021). Moreover, large amounts of methane released by hydrate decomposition may have entered the ocean and atmosphere, which would have contributed positively to the sudden warming of the climate during the penultimate deglaciation.

Data Availability Statement

The original contributions presented in the study are included in the article/supplementary material. Further inquiries can be directed to the corresponding authors.

Author Contributions

JW collected samples. XM formal analysis. JW, TW, XM, and PS writing–review and editing. PS and JW funding acquisition. JW and TW wrote the paper with contributions from all the co-authors. All authors contributed to the article and approved the submitted version.

Funding

This research was funded by Key Special Project for Introduced Talents Team of Southern Marine Science and Engineering Guangdong Laboratory (Guangzhou) (GML2019ZD0201), Project of Hubei Key Laboratory Marine Geological Resources (MGR202002), China Geological Survey Project (No. DD20160227).

Conflict of Interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Publisher’s Note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

References

Algeo T., Tribovillard N. (2009). Environmental Analysis of Paleoceanographic Systems Based on Molybdenum Uranium Covariation. Chem. Geol. 268, 211–225. doi: 10.1016/j.chemgeo.2009.09.001