Database of Ice-Rich Yedoma Permafrost (IRYP)

Vast portions of Arctic and sub-Arctic Siberia, Alaska and the Yukon Territory are covered by ice-rich silty to sandy deposits that are containing large ice wedges, resulting from syngenetic sedimentation and freezing. Accompanied by wedge-ice growth in polygonal landscapes, the sedimentation process was driven by cold continental climatic and environmental conditions in unglaciated regions during the late Pleistocene, inducing the accumulation of the unique Yedoma deposits up to >50 meters thick. Because of fast incorporation of organic material into syngenetic permafrost during its formation, Yedoma deposits include well-preserved organic matter. Ice-rich deposits like Yedoma are especially prone to degradation triggered by climate changes or human activity. When Yedoma deposits degrade, large amounts of sequestered organic carbon as well as other nutrients are released and become part of active biogeochemical cycling. This could be of global significance for future climate warming as increased permafrost thaw is likely to lead to a positive feedback through enhanced greenhouse gas fluxes. Therefore, a detailed assessment of the current Yedoma deposit coverage and its volume is of importance to estimate its potential response to future climate changes. We synthesized the map of the coverage and thickness estimation, which will provide critical data needed for further research. In particular, this preliminary Yedoma map is a great step forward to understand the spatial heterogeneity of Yedoma deposits and its regional coverage. There will be further applications in the context of reconstructing paleo-environmental dynamics and past ecosystems like the mammoth-steppe-tundra, or ground ice distribution including future thermokarst vulnerability. Moreover, the map will be a crucial improvement of the data basis needed to refine the present-day Yedoma permafrost organic carbon inventory, which is assumed to be between 83±12 (Strauss et al., 2013, doi:10.1002/2013GL058088) and 129±30 (Walter Anthony et al., 2014, doi:10.1038/nature13560) gigatonnes (Gt) of organic carbon in perennially-frozen archives. Hence, here we synthesize data on the circum-Arctic and sub-Arctic distribution and thickness of Yedoma for compiling a preliminary circum-polar Yedoma map. For compiling this map, we used (1) maps of the previous Yedoma coverage estimates, (2) included the digitized areas from Grosse et al. (2013) as well as extracted areas of potential Yedoma distribution from additional surface geological and Quaternary geological maps (1.: 1:500,000: Q-51-V,G; P-51-A,B; P-52-A,B; Q-52-V,G; P-52-V,G; Q-51-A,B; R-51-V,G; R-52-V,G; R-52-A,B; 2.: 1:1,000,000: P-50-51; P-52-53; P-58-59; Q-42-43; Q-44-45; Q-50-51; Q-52-53; Q-54-55; Q-56-57; Q-58-59; Q-60-1; R-(40)-42; R-43-(45); R-(45)-47; R-48-(50); R-51; R-53-(55); R-(55)-57; R-58-(60); S-44-46; S-47-49; S-50-52; S-53-55; 3.: 1:2,500,000: Quaternary map of the territory of Russian Federation, 4.: Alaska Permafrost Map). The digitalization was done using GIS techniques (ArcGIS) and vectorization of raster Images (Adobe Photoshop and Illustrator). Data on Yedoma thickness are obtained from boreholes and exposures reported in the scientific literature. The map and database are still preliminary and will have to undergo a technical and scientific vetting and review process. In their current form, we included a range of attributes for Yedoma area polygons based on lithological and stratigraphical information from the original source maps as well as a confidence level for our classification of an area as Yedoma (3 stages: confirmed, likely, or uncertain). In its current version, our database includes more than 365 boreholes and exposures and more than 2000 digitized Yedoma areas. We expect that the database will continue to grow. In this preliminary stage, we estimate the Northern Hemisphere Yedoma deposit area to cover approximately 625,000 km². We estimate that 53% of the total Yedoma area today is located in the tundra zone, 47% in the taiga zone. Separated from west to east, 29% of the Yedoma area is found in North America and 71 % in North Asia. The latter include 9% in West Siberia, 11% in Central Siberia, 44% in East Siberia and 7% in Far East Russia. Adding the recent maximum Yedoma region (including all Yedoma uplands, thermokarst lakes and basins, and river valleys) of 1.4 million km² (Strauss et al., 2013, doi:10.1002/2013GL058088) and postulating that Yedoma occupied up to 80% of the adjacent formerly exposed and now flooded Beringia shelves (1.9 million km², down to 125 m below modern sea level, between 105°E - 128°W and >68°N), we assume that the Last Glacial Maximum Yedoma region likely covered more than 3 million km² of Beringia. Acknowledgements: This project is part of the Action Group "The Yedoma Region: A Synthesis of Circum-Arctic Distribution and Thickness" (funded by the International Permafrost Association (IPA) to J. Strauss) and is embedded into the Permafrost Carbon Network (working group Yedoma Carbon Stocks). We acknowledge the support by the European Research Council (Starting Grant #338335), the German Federal Ministry of Education and Research (Grant 01DM12011 and "CarboPerm" (03G0836A)), the Initiative and Networking Fund of the Helmholtz Association (#ERC-0013) and the German Federal Environment Agency (UBA, project UFOPLAN FKZ 3712 41 106).

Characteristics of samples obtained during Pokhodsk 2012-2013 campaigns in the joint Russian-German POLYGON Project

Polygonal tundra, thermokarst basins and pingos are common and characteristic periglacial features of arctic lowlands underlain by permafrost in Northeast Siberia. Modern polygonal mires are in the focus of biogeochemical, biological, pedological, and cryolithological research with special attention to their carbon stocks and greenhouse-gas fluxes, their biodiversity and their dynamics and functioning under past, present and future climate scenarios. Within the frame of the joint German-Russian DFG-RFBR project Polygons in tundra wetlands: state and dynamics under climate variability in Polar Regions (POLYGON) field studies of recent and of late Quaternary environmental dynamics were carried out in the Indigirka lowland and in the Kolyma River Delta in summer 2012 and summer 2013. Using a multidisciplinary approach, several types of polygons and thermokarst lakes were studied in different landscapes units in the Kolyma Delta in 2012 around the small fishing settlement Pokhodsk. The floral and faunal associations of polygonal tundra were described during the fieldwork. Ecological, hydrological, meteorological, limnological, pedological and cryological features were studied in order to evaluate modern and past environmental conditions and their essential controlling parameters. The ecological monitoring and collection program of polygonal ponds were undertaken as in 2011 in the Indigirka lowland by a former POLYGON expedition (Schirrmeister et al. [eds.] 2012). Exposures, pits and drill cores in the Kolyma Delta were studied to understand the cryolithological structures of frozen ground and to collect samples for detailed paleoenvironmental research of the late Quaternary past. Dendrochronological and ecological studies were carried out in the tree line zone south of the Kolyma Delta. Based on previous work in the Indigirka lowland in 2011 (Schirrmeister et al. [eds.] 2012), the environmental monitoring around the Kytalyk research station was continued until the end of August 2012. In addition, a classical exposure of the late Pleistocene permafrost at the Achchaygy Allaikha River near Chokurdakh was studied. The ecological studies near Pokhodsk were continued in 2013 (chapter 13). Other fieldwork took place at the Pokhodsk-Yedoma-Island in the northwestern part of the Kolyma Delta.

Tab. 1 Mole % CO2 and 13C values of CO2 at baseline and at the final monitoring event for the eight observation wells

During CO2 storage operations in mature oilfields or saline aquifers it is desirable to trace the movement of injected CO2 for verification and safety purposes. We demonstrate the successful use of carbon isotope abundance ratios for tracing the movement of CO2 injected at the Cardium CO2 Storage Monitoring project in Alberta between 2005 and 2007. Injected CO2 had a d13C value of -4.6±1.1 per mil that was more than 10 per mil higher than the carbon isotope ratios of casing gas CO2 prior to CO2 injection with average d13C values ranging from -15.9 to -23.5 per mil. After commencement of CO2 injection, d13C values of casing gas CO2 increased in all observation wells towards those of the injected CO2 consistent with a two-source end-member mixing model. At four wells located in a NE-SW trend with respect to the injection wells, breakthrough of injected CO2 was registered chemically (>50 mol % CO2) and isotopically 1-6 months after commencement of CO2 injection resulting in cumulative CO2 fluxes exceeding 100000 m**3 during the observation period. At four other wells, casing gas CO2 contents remained below 5 mol % resulting in low cumulative CO2 fluxes (<2000 m**3) throughout the entire observation period, but carbon isotope ratios indicated contributions between <30 and 80% of injected CO2. Therefore, we conclude that monitoring the movement of CO2 in the injection reservoir with geochemical and isotopic techniques is an effective approach to determine plume expansion and to identify potential preferential flow paths provided that the isotopic composition of injected CO2 is constant and distinct from that of baseline CO2.

Gas composition, concentration and stable isotope ratios (Table 2, 3, 4)

The deployment of CCS (carbon capture and storage) at industrial scale implies the development of effective monitoring tools. Noble gases are tracers usually proposed to track CO2. This methodology, combined with the geochemistry of carbon isotopes, has been tested on available analogues. At first, gases from natural analogues were sampled in the Colorado Plateau and in the French carbogaseous provinces, in both well-confined and leaking-sites. Second, we performed a 2-years tracing experience on an underground natural gas storage, sampling gas each month during injection and withdrawal periods. In natural analogues, the geochemical fingerprints are dependent on the containment criterion and on the geological context, giving tools to detect a leakage of deep-CO2 toward surface. This study also provides information on the origin of CO2, as well as residence time of fluids within the crust and clues on the physico-chemical processes occurring during the geological story. The study on the industrial analogue demonstrates the feasibility of using noble gases as tracers of CO2. Withdrawn gases follow geochemical trends coherent with mixing processes between injected gas end-members. Physico-chemical processes revealed by the tracing occur at transient state. These two complementary studies proved the interest of geochemical monitoring to survey the CO2 behaviour, and gave information on its use.

Specific microbial biomarker signatures from Congo fan sediments of ODP site 175-1075, Congo soils, wetland and estuary sediments

Methane (CH4) is a strong greenhouse gas known to have perturbed global climate in the past, especially when released in large quantities over short time periods from continental or marine sources. It is therefore crucial to understand and, if possible, quantify the individual and combined response of these variable methane sources to natural climate variability. However, past changes in the stability of greenhouse gas reservoirs remain uncertain and poorly constrained by geological evidence. Here, we present a record from the Congo fan of a highly specific bacteriohopanepolyol (BHP) biomarker for aerobic methane oxidation (AMO), 35-aminobacteriohopane-30,31,32,33,34-pentol (aminopentol), that identifies discrete periods of increased AMO as far back as 1.2 Ma. Fluctuations in the concentration of aminopentol, and other 35-aminoBHPs, follow a pattern that correlates with late Quaternary glacial-interglacial climate cycles, with highest concentrations during warm periods. We discuss possible sources of aminopentol, and the methane consumed by the precursor methanotrophs, within the context of the Congo River setting, including supply of methane oxidation markers from terrestrial watersheds and/or marine sources (gas hydrate and/or deep subsurface gas reservoir). Compound-specific carbon isotope values of -30 per mil to -40 per mil for BHPs in ODP 1075 and strong similarities between the BHP signature of the core and surface sediments from the Congo estuary and floodplain wetlands from the interior of the Congo River Basin, support a methanotrophic and likely terrigenous origin of the 35-aminoBHPs found in the fan sediments. This new evidence supports a causal connection between marine sediment BHP records of tropical deep sea fans and wetland settings in the feeding river catchments, and thus tropical continental hydrology. Further research is needed to better constrain the different sources and pathways of methane emission. However, this study identifies the large potential of aminoBHPs, in particular aminopentol, to trace and, once better calibrated and understood, quantify past methane sources and fluxes from terrestrial and potentially also marine sources.

Methane and sulfide fluxes in permanent anoxia: in situ studies at the Dvurechenskii mud volcano (Sorokin Trough, Black Sea)

The Dvurechenskii mud volcano (DMV), located in permanently anoxic waters at 2060 m depth (Sorokin Trough, Black Sea), was visited during the M72/2 cruise with the RV Meteor to investigate the methane and sulfide release from mud volcanoes into the Black Sea hydrosphere. We studied benthic fluxes of methane and sulfide, and the factors controlling transport, consumption and production of both compounds within the sediment. The pie shaped mud volcano showed temperature anomalies as well as solute and gas fluxes indicating high fluid flow at a small elevation north of the geographical center. The anaerobic oxidation of methane (AOM) coupled to sulfate reduction (SR) was excluded from this zone due to fluid-flow induced sulfate limitation and a fresh mud flow and consequently methane escaped into the water column with a rate of 0.46 mol/m**2/d. In the outer center of the mud volcano fluid flow and total methane flux were decreased, correlating with an increase in sulfate penetration into the sediment, and with higher SR and AOM rates. Here between 50-70% of the methane flux (0.07-0.1 mol/m**2/d) was consumed within the upper 10 cm of the sediment. Also at the edge of the mud volcano fluid flow and rates of methane and sulfate turnover were substantial. The overall amount of dissolved methane released from the mud volcano into the water column was significant with a discharge of 1.4x10**7 mol/yr. The DMV maintains also high areal rates of methane-fueled sulfide production of on average 0.05 mol/m**2/d. However, we concluded that sulfide and methane emission into the hydrosphere from deep water mud volcanoes does not significantly contribute to the sulfide and methane inventory of the Black Sea.

Permanent open water surfaces in central Siberia from ENVISAT ASAR wide swath data 2003/2004, link to shapefile

Permanent water bodies not only store dissolved CO2 but are essential for the maintenance of wetlands in their proximity. From the viewpoint of greenhouse gas (GHG) accounting wetland functions comprise sequestration of carbon under anaerobic conditions and methane release. The investigated area in central Siberia covers boreal and sub-arctic environments. Small inundated basins are abundant on the sub-arctic Taymir lowlands but also in parts of severe boreal climate where permafrost ice content is high and feature important freshwater ecosystems. Satellite radar imagery (ENVISAT ScanSAR), acquired in summer 2003 and 2004, has been used to derive open water surfaces with 150 m resolution, covering an area of approximately 3 Mkm**2. The open water surface maps were derived using a simple threshold-based classification method. The results were assessed with Russian forest inventory data, which includes detailed information about water bodies. The resulting classification has been further used to estimate the extent of tundra wetlands and to determine their importance for methane emissions. Tundra wetlands cover 7% (400,000 km**2) of the study region and methane emissions from hydromorphic soils are estimated to be 45,000 t/d for the Taymir peninsula.

Physical properties of 39 sediment cores from the Bengal Fan and northeast Indian Ocean

We obtained sediment physical properties and geochemical data from 47 piston and gravity cores located in the Bay of Bengal, to study the complex history of the Late Pleistocene run-off from the Ganges and Brahmaputra rivers and its imprint on the Bengal Fan. Grain-size parameters were predicted from core logs of density and velocity to infer sediment transport energy and to distinguish different environments along the 3000-km-long transport path from the delta platform to the lower fan. On the shelf, 27 cores indicate rapidly prograding delta foresets today that contain primarily mud, whereas outer shelf sediment has 25% higher silt contents, indicative of stronger and more stable transport regime, which prevent deposition and expose a Late Pleistocene relic surface. Deposition is currently directed towards the shelf canyon 'Swatch of No Ground', where turbidites are released to the only channel-levee system that is active on the fan during the Holocene. Active growth of the channel-levee system occurred throughout sea-level rise and highstand with a distinct growth phase at the end of the Younger Dryas. Coarse-grained material bypasses the upper fan and upper parts of the middle fan, where particle flow is enhanced as a result of flow-restriction in well-defined channels. Sandier material is deposited mainly as sheet-flow deposits on turbidite-dominated plains at the lower fan. The currently most active part of the fan with 10-40 cm thick turbidites is documented for the central channel including inner levees (e.g., site 40). Site 47 from the lower fan far to the east of the active channel-levee system indicates the end of turbidite sedimentation at 300 ka for that location. That time corresponds to the sea-level lowering during late isotopic stage 9 when sediment supply to the fan increased and led to channel avulsion farther upstream, probably indicating a close relation of climate variability and fan activity. Pelagic deep-sea sites 22 and 28 contain a 630-kyear record of climate response to orbital forcing with dominant 21- and 41-kyear cycles for carbonate and magnetic susceptibility, respectively, pointing to teleconnections of low-latitude monsoonal forcing on the precession band to high-latitude obliquity forcing. Upper slope sites 115, 124, and 126 contain a record of the response to high-frequency climate change in the Dansgaard-Oeschger bands during the last glacial cycle with shared frequencies between 0.75 and 2.5 kyear. Correlation of highs in Bengal Fan physical properties to lows in the d18O record of the GISP2 ice-core suggests that times of greater sediment transport energy in the Bay of Bengal are associated with cooler air temperatures over Greenland. Teleconnections were probably established through moisture and other greenhouse-gas forcing that could have been initiated by instabilities in the methane hydrate reservoir in the oceans.

Planktonic foraminiferal oxygen isotopes of ODP Site 143-865 samples

Cool tropical sea surface temperatures (SSTs) are reported for warm Paleogene greenhouse climates based on the d18O of planktonic foraminiferal tests. These results are difficult to reconcile with models of greenhouse gas-forced climate. It has been suggested that this "cool tropics paradox" arises from postdepositional alteration of foraminiferal calcite, yielding erroneously high d18O values. Recrystallization of foraminiferal tests is cryptic and difficult to quantify, and the compilation of robust d18O records from moderately altered material remains challenging. Scanning electron microscopy of planktonic foraminiferal chamber-wall cross sections reveals that the basal area of muricae, pustular outgrowths on the chamber walls of species belonging to the genus Morozovella, contain no mural pores and may be less susceptible to postdepositional alteration. We analyzed the d18O in muricae bases of morozovellids from the central Pacific (Ocean Drilling Program Site 865) by ion microprobe using 10 ?m pits with an analytical reproducibility of ±0.34 per mil (2 standard deviations). In situ measurements of d18O in these domains yield consistently lower values than those published for conventional multispecimen analyses. Assuming that the original d18O is largely preserved in the basal areas of muricae, this new d18O record indicates Early Paleogene (~49-56 Ma) tropical SSTs in the central Pacific were 4°-8°C higher than inferred from the previously published d18O record and that SSTs reached at least ~33°C during the Paleocene-Eocene thermal maximum. This study demonstrates the utility of ion microprobe analysis for generating more reliable paleoclimate records from moderately altered foraminiferal tests preserved in deep-sea sediments.

Planktonic foraminifera, nannofossils and stable isotope record across tghe Paleocene-Eocene thermal maximum of ODP Hole 113-690B

The carbonate saturation profile of the oceans shoaled markedly during a transient global warming event known as the Paleocene-Eocene thermal maximum (PETM) (circa 55 Ma). The rapid release of large quantities of carbon into the ocean-atmosphere system is believed to have triggered this intense episode of dissolution along with a negative carbon isotope excursion (CIE). The brevity (120-220 kyr) of the PETM reflects the rapid enhancement of negative feedback mechanisms within Earth's exogenic carbon cycle that served the dual function of buffering ocean pH and reducing atmospheric greenhouse gas levels. Detailed study of the PETM stratigraphy from Ocean Drilling Program Site 690 (Weddell Sea) reveals that the CIE recovery period, which postdates the CIE onset by ~80 kyr, is represented by an expanded (~2.5 m thick) interval containing a unique planktic foraminiferal assemblage strongly diluted by coccolithophore carbonate. Collectively, the micropaleontological and sedimentological changes preserved within the CIE recovery interval reflect a transient state when ocean-atmosphere chemistry fostered prolific coccolithophore blooms that suppressed the local lysocline to relatively deeper depths. A prominent peak in the abundance of the clay mineral kaolinite is associated with the CIE recovery interval, indicating that continental weathering/runoff intensified at this time as well (Robert and Kennett, 1994). Such parallel stratigraphic changes are generally consonant with the hypothesis that enhanced continental weathering/runoff and carbonate precipitation helped sequester carbon during the PETM recovery period (e.g., Dickens et al., 1997, doi:10.1130/0091-7613(1997)025<0259:ABOGIT>2.3.CO;2 ; Zachos et al., 2005, doi:10.1126/science.1109004).

Sea surface temperatures during Marine Isotope Stage (MIS) 11

The Marine Isotope Stage (MIS) 11 (424-374 ka) was characterized by a protracted deglaciation and an unusually long climatic optimum. It remains unclear to what degree the climate development during this interglacial reflects the unusually weak orbital forcing or greenhouse gas trends. Previously, arguments about the duration and timing of the MIS11 climatic optimum and about the pace of the deglacial warming were based on a small number of key records, which appear to show regional differences. In order to obtain a global signal of climate evolution during MIS11, we compiled a database of 78 sea surface temperature (SST) records from 57 sites spanning MIS11, aligned these individually on the basis of benthic (N = 28) or planktonic (N = 31) stable oxygen isotope curves to a common time frame and subjected 48 of them to an empirical orthogonal function (EOF) analysis. The analysis revealed a high commonality among all records, with the principal SST trend explaining almost 49% of the variability. This trend indicates that on the global scale, the surface ocean underwent rapid deglacial warming during Termination V, in pace with carbon dioxide rise, followed by a broad SST optimum centered at ~410 kyr. The second EOF, which explained ~18% of the variability, revealed the existence of a different SST trend, characterized by a delayed onset of the temperature optimum during MIS11 at ~398 kyr, followed by a prolonged warm period lasting beyond 380 kyr. This trend is most consistently manifested in the mid-latitude North Atlantic and Mediterranean Sea and is here attributed to the strength of the Atlantic meridional overturning circulation. A sensitivity analysis indicates that these results are robust to record selection and to age-model uncertainties of up to 3-6 kyr, but more sensitive to SST seasonal attribution and SST uncertainties >1 °C. In order to validate the CCSM3 (Community Climate System Model, version 3) predictive potential, the annual and seasonal SST anomalies recorded in a total of 74 proxy records were compared with runs for three time slices representing orbital configuration extremes during the peak interglacial of MIS11. The modeled SST anomalies are characterized by a significantly lower variance compared to the reconstructions. Nevertheless, significant correlations between proxy and model data are found in comparisons on the seasonal basis, indicating that the model captures part of the long-term variability induced by astronomical forcing, which appears to have left a detectable signature in SST trends.

Pore water and sediment geochemistry of the Bothian Sea

Methane is a powerful greenhouse gas and its biological conversion in marine sediments, largely controlled by anaerobic oxidation of methane (AOM), is a crucial part of the global carbon cycle. However, little is known about the role of iron oxides as an oxidant for AOM. Here we provide the first field evidence for iron-dependent AOM in brackish coastal surface sediments and show that methane produced in Bothnian Sea sediments is oxidized in distinct zones of iron- and sulfate-dependent AOM. At our study site, anthropogenic eutrophication over recent decades has led to an upward migration of the sulfate/methane transition zone in the sediment. Abundant iron oxides and high dissolved ferrous iron indicate iron reduction in the methanogenic sediments below the newly established sulfate/methane transition. Laboratory incubation studies of these sediments strongly suggest that the in situ microbial community is capable of linking methane oxidation to iron oxide reduction. Eutrophication of coastal environments may therefore create geochemical conditions favorable for iron-mediated AOM and thus increase the relevance of iron-dependent methane oxidation in the future. Besides its role in mitigating methane emissions, iron-dependent AOM strongly impacts sedimentary iron cycling and related biogeochemical processes through the reduction of large quantities of iron oxides.

Simulated change in atmospheric CO2 and D14C around 14.6 kyr BP, at the onset of the northern hemispheric warming into the Bølling/Allerød

One of the most abrupt and yet unexplained past rises in atmospheric CO2 (10 p.p.m.v. in two centuries) occurred in quasi-synchrony with abrupt northern hemispheric warming into the Bølling/Allerød, 14,600 years ago. Here we use a U/Th-dated record of atmospheric D14C from Tahiti corals to provide an independent and precise age control for this CO2 rise. We also use model simulations to show that the release of old (nearly 14C-free) carbon can explain these changes in CO2 and D14C. The D14C record provides an independent constraint on the amount of carbon released (125 Pg C). We suggest, in line with observations of atmospheric CH4 and terrigenous biomarkers, that thawing permafrost in high northern latitudes could have been the source of carbon, possibly with contribution from flooding of the Siberian continental shelf during meltwater pulse 1A. Our findings highlight the potential of the permafrost carbon reservoir to modulate abrupt climate changes via greenhouse-gas feedbacks.

Stable carbon isotopes, Methane Index, and distributions of different GDGTs

Gas hydrates represent one of the largest pools of readily exchangeable carbon on Earth's surface. Releases of the greenhouse gas methane from hydrates are proposed to be responsible for climate change at numerous events in geological history. Many of these inferred events, however, were based on carbonate carbon isotopes which are susceptible to diagenetic alterations. Here we propose a molecular fossil proxy, i.e., the "Methane Index (MI)", to detect and document the destabilization and dissociation of marine gas hydrates. MI consists of the relative distribution of glycerol dibiphytanyl glycerol tetraethers (GDGTs), the core membrane lipids of archaea. The rational behind MI is that in hydrate-impacted environments, the pool of archaeal tetraether lipids is dominated by GDGT-1, -2 and -3 due to the large contribution of signals from the methanotrophic archaeal community. Our study in the Gulf of Mexico cold-seep sediments demonstrates a correlation between MI and the compound-specific carbon isotope of GDGTs, which is strong evidence supporting the MI-methane consumption relationship. Preliminary applications of MI in a number of hydrate-impacted and/or methane-rich environments show diagnostic MI values, corroborating the idea that MI may serve as a robust indicator for hydrate dissociation that is useful for studies of global carbon cycling and paleoclimate change.