Organic matter quality of deep permafrost carbon – a study from Arctic Siberia

J. Strauss, L. Schirrmeister, K. Mangelsdorf, L. Eichhorn, S. Wetterich, and U. Herzschuh Alfred Wegener Institute Helmholtz Centre for Polar and Marine Research, Periglacial Research Unit Potsdam, Telegrafenberg A 43, Potsdam, Germany Helmholtz Centre Potsdam GFZ German Research Centre for Geosciences, Telegrafenberg, 14473 Potsdam, Germany International Max Planck Research School for Global Biogeochemical Cycles (Max-Planck-Institute for Biogeochemistry and Friedrich Schiller University), Burgweg 11, 07749 Jena, Germany

Abstract.The organic-carbon (OC) pool accumulated in Arctic permafrost (perennially frozen ground) equals the carbon stored in the modern atmosphere.To give an idea of how Yedoma region permafrost could respond under future climatic warming, we conducted a study to quantify the organic-matter quality (here defined as the intrinsic potential to be further transformed, decomposed, and mineralized) of late Pleistocene (Yedoma) and Holocene (thermokarst) deposits on the Buor-Khaya Peninsula, northeast Siberia.The objective of this study was to develop a stratigraphic classified organic-matter quality characterization.For this purpose the degree of organic-matter decomposition was estimated by using a multiproxy approach.We applied sedimentological (grain-size analyses, bulk density, ice content) and geochemical parameters (total OC, stable carbon isotopes (δ 13 C), total organic carbon : nitrogen (C / N) ratios) as well as lipid biomarkers (n-alkanes, n-fatty acids, hopanes, triterpenoids, and biomarker indices, i.e., average chain length, carbon preference index (CPI), and higher-plant fatty-acid index (HPFA)).Our results show that the Yedoma and thermokarst organic-matter qualities for further decomposition exhibit no obvious degradation-depth trend.Relatively, the C / N and δ 13 C values and the HPFA index show a significantly better preservation of the organic matter stored in thermokarst deposits compared to Yedoma deposits.The CPI data suggest less degradation of the organic matter from both deposits, with a higher value for Yedoma organic matter.As the interquartile ranges of the proxies mostly over-lap, we interpret this as indicating comparable quality for further decomposition for both kinds of deposits with likely better thermokarst organic-matter quality.Supported by principal component analyses, the sediment parameters and quality proxies of Yedoma and thermokarst deposits could not be unambiguously separated from each other.This revealed that the organic-matter vulnerability is heterogeneous and depends on different decomposition trajectories and the previous decomposition and preservation history.Elucidating this was one of the major new contributions of our multiproxy study.With the addition of biomarker data, it was possible to show that permafrost organic-matter degradation likely occurs via a combination of (uncompleted) degradation cycles or a cascade of degradation steps rather than as a linear function of age or sediment facies.We conclude that the amount of organic matter in the studied sediments is high for mineral soils and of good quality and therefore susceptible to future decomposition.The lack of depth trends shows that permafrost acts like a giant freezer, preserving the constant quality of ancient organic matter.When undecomposed Yedoma organic matter is mobilized via thermokarst processes, the fate of this carbon depends largely on the environmental conditions; the carbon could be preserved in an undecomposed state till refreezing occurs.If modern input has occurred, thermokarst organic matter could be of a better quality for future microbial decomposition than that found in Yedoma deposits.

Introduction
During the late Quaternary, the rate of organic-matter decomposition in the Arctic was slower than plant growth, sedimentation, and freezing rates.Thus, a large pool of organic carbon (OC) accumulated in the Arctic and was sequestered deep in the permafrost.Hugelius et al. (2014) estimates an OC storage of 1300 Gt for the circum-Arctic permafrost region with ∼ 850 Gt OC sequestered in permafrost.This is approximately the carbon stored in the modern atmosphere (Dlugokencky and Tans, 2014).During warming and permafrost thawing, this formerly cryo-sequestered OC gradually entered the modern biogeochemical cycle by microbial turnover.By thawing and microbial activity, the permafrost deposits can turn from a carbon sink to a source (Schuur et al., 2009), releasing greenhouse gases such as carbon dioxide and methane to the atmosphere.In addition to the near-surface carbon pool representing the uppermost 3 m below surface, deep OC pools, especially those held in ice-rich permafrost deposits in the Yedoma region, are of great significance for current concerns about the effects of global warming because of rapid permafrost thaw, such as thermokarst and thermoerosion.According to Strauss et al. (2013) and Hugelius et al. (2014), the Yedoma region is defined as the area of potential distribution of late Pleistocene ice-rich and organic-rich silty deposits (Yedoma), covering large areas in Siberia and Alaska.Estimates of OC stored in the Yedoma region amount to 83 ± 12 Gt for late Pleistocene Yedoma deposits (ages shown in Table 1).Due to Holocene warming, subsequent ground ice melt and surface subsidence, thermokarst basins formed and were partly occupied by lakes.Holocene thermokarst deposits (ages shown in Table 1) contain 130 ± 29 Gt organic carbon.In total, the Yedoma region extends to an area of about 1 387 000 km 2 , of which about 70 % is already affected by permafrost degradation (thermokarst; Strauss et al., 2013).Kuhry et al. (2009) and Schirrmeister et al. (2011a) showed that Yedoma deposits accumulated at fast rates, implying a short time for the organic matter to decay before it became locked into a perennially frozen state.Therefore, the organic-matter availability for microorganisms is expected to be excellent, resulting in great vulnerability to warming ground conditions (Mu et al., 2014).To elucidate how the Yedoma region permafrost could respond under conditions of future climatic warming, we studied the organic-matter degradation state of Yedoma and its Holocene degradation features (called thermokarst deposits) on the Buor-Khaya Peninsula, eastern Laptev Sea.As mentioned above, Strauss et al. (2013) found that thermokarst deposits contain the quantitatively more important carbon pool, but the unsolved question is this: is the thermokarst organic-matter pool as degradable as the frozen late Pleistocene Yedoma, or has the most labile carbon already been emitted due to thermokarst degradation processes?In both kinds of deposits the OC was incorporated deep (deeper than 3 m) into permafrost (Schirrmeister et al., 2013;Strauss et al., 2013).As shown by models and extrapolation from recent observations, the more southern portions of Yedoma deposits thawed during the last deglaciation, resulting in large emissions of greenhouse gases to the atmosphere (Walter et al., 2007a;Ciais et al., 2012;Walter Anthony et al., 2014).Recent ground warming has been observed in the permafrost zone (Romanovsky et al., 2010), and incubation experiments reveal that permafrost warming is accompanied by a substantial outgassing of greenhouse gases (Lee et al., 2012;Knoblauch et al., 2013;Schädel et al., 2014).As an illustration of the important influence of ground temperature on organic-matter degradation, a higher respiration rate at greater depth close to the permafrost table (Mangelsdorf et al., 2009;Waldrop et al., 2010) was found inside the seasonally thawed active layer and interpreted as a greater lability of the organic matter close to the perennially frozen ground.Focusing on permafrost deposits in the Laptev Sea region, which includes our Buor-Khaya study site, Schirrmeister et al. (2011a) characterize the Yedoma region permafrost organic matter as weakly decomposed.
Biomarkers are used for the paleoenvironmental reconstruction of terrestrial permafrost (Andersson et al., 2011) or the characterization of permafrost organic-matter degradation (Andersson and Meyers, 2012;Vonk et al., 2013;Routh et al., 2014).In our study we estimate molecular markers (n-alkanes, n-fatty acids, hopanes, and triterpenoids) and use biomarker proxies and indices (absolute lipid concentration, average chain length (ACL), carbon preference index (CPI), hop-17(21)-ene, higher-plant fatty-acid (HPFA) index, and an oleanene ratio) to test whether they are useful mirrors of organic-matter decomposition, i.e., of the organicmatter state of degradation in permafrost deposits.Fairly established methods, both cryolithological (grain-size analyses, bulk density, ice content) and biogeochemical (total organic carbon (TOC wt % ), stable carbon isotope ratios (δ 13 C in TOC), total nitrogen (TN), and TOC wt % / TN (C / N) ratios), are applied to our sample set.Finally, principal components analysis (PCA) highlights the relationships between different organic-matter degradation proxies.
Because the future feedback from the Yedoma region permafrost OC to climate forcing is driven by both (1) the pool size, estimated to be ∼ 211 Gt (Strauss et al., 2013), and (2) the state of degradation of OC stored in the studied deposits, the objective of this study is the development of a stratigraphically differentiated organic-matter quality characterization using sample material representative of widespread Yedoma and thermokarst permafrost.We hypothesize increased organic-matter degradation during thermokarst processes but also increased organic-matter input during climatically favorable Holocene times.
Table 1.Radiocarbon dating on plant macroremains.Calibrations were done by using the Calib 6.0 software and the IntCal09 calibration curve (Stuiver et al., 2010).Depth is given in meters below surface level (m b.s.l.) and height in meters above sea level (m a.s.l. 2 Materials and methods

Study area
The Buor-Khaya Peninsula study site (71 • 34 N, 132 • 12 E) is located in the northeastern part of Siberia (Fig. 1).The Buor-Khaya Peninsula is framed by the Laptev Sea, a shallow epicontinental part of the Arctic Ocean, and geologically by two rift structures (Drachev et al., 1998).Buor-Khaya is underlain by continuous permafrost with ground temperatures of less than −11 • C (Drozdov et al., 2005).The permafrost thickness is estimated to be between 450 and 650 m (Romanovskii et al., 2004).Stratigraphically, outcrops from two sediment units are distinguished and studied: (1) ice-rich permafrost, called Yedoma deposits, and (2) deposits in permafrost rapid-thaw features, generalized as thermokarst deposits.Three profiles of thermokarst deposits (in a thermokarst basin: Buo-01 and Buo-05; initial thermokarst on top of a Yedoma hill: Buo-03) and two profiles of Yedoma deposits (Buo-02, Buo-04) were studied and sampled.Figure 1 shows an overview of the sampled profiles and their position relative to each other.

Field work
Field studies were undertaken in summer 2010 at outcrops situated on the western coast of the Buor-Khaya Peninsula.
The sediment of the profiles and sub-profiles, exposed at the cliff wall or partly in thermokarst mounds in thaw slumps, was dug with spades and cleaned with hacks.The cry-olithology, sediment characteristics, and visible organic matter in the sediments of the sequences chosen were surveyed and described.Moreover, the profiles were photographed and sketched.Sub-profiles were stacked together to create composite profiles.Sampling positions in neighboring subprofiles were correlated by height estimation using measuring tape.The upper edge of each profile was calibrated with tacheometer measurements (Günther et al., 2012).In the field laboratory, all sample volumes were measured with a balance following Archimedes principle, and the absolute ice content was determined by drying the sample.In total, 91 samples were taken and cooled for transport to laboratories for further analysis.The detailed sampling positions for each profile are shown in Strauss and Schirrmeister (2011).

Indicators of organic-matter quality for further decomposition
To validate and to extend the sedimentological approach used and to estimate the organic-matter quality for further decomposition, lipid biomarkers were measured in order to estimate the degree of organic-matter degradation.For biomarker studies we used a "fingerprint" approach by focusing on identifiable markers related to the state of organicmatter degradation.Below, the geochemical indicators and biomarkers used are described.

Grain-size analyses
Grain sizes were analyzed using a laser particle sizer (LS 200, Beckmann-Coulter) with a measuring range of 0.375 to 1000 µm (Figs. 2, S1 in the Supplement).Grain-size calculations were done following Folk and Ward (1957) using Gradistat v8 (Blott and Pye, 2001).A detailed description of these analytical techniques is given in the Supplement (Sect.S1.1 in the Supplement).

Elemental composition
To determine the total elemental carbon and total nitrogen (TN) content, the samples were measured by a carbonnitrogen-sulfur analyzer (Vario EL III, Elementar).TOC wt % was measured with a TOC analyzer (Vario Max C, Elementar).The volumetric TOC content (TOC kg m −3 ) was calculated according to Strauss et al. (2013).A detailed description of this technique is given in the Supplement (Sect.S1.2).
The TOC / TN (C / N) ratio has been used as a general indicator of the degree of organic-matter decomposition (Stevenson, 1994).Based on the assumption that organicmatter components are degraded selectively, degradation modifies elemental compositions and hence C / N in deposits.Because a decrease in the C / N ratio has been ob-served in aerated deposits with microbial immobilization of TN (nitrogen stays in the system) accompanied by the remineralization of TOC (Sollins et al., 1984) and CO 2 emission, this ratio is used in the following way: the higher the C / N ratio, the lower the degree of decomposition.

Bulk density and volumetric carbon content
BD (bulk density) was calculated using Eq.(1).
Estimating the BD is required in order to convert the measured-weight-based TOC wt % content per sample to a volume-based value.Thus, the TOC kg m −3 was calculated according to Eq. ( 2): TOC wt% 100 . (2)

Carbon isotope studies
Stable TOC carbon isotopes were determined with a Finnigan MAT Delta-S mass spectrometer combined with a FLASH elemental analyzer and a ConFlo III gas mixing system.A detailed methodology is given in the Supplement (Sect.S1.4).The stable carbon isotopes of OC reflect (1) an initial contribution from different plant species and plant components and (2) subsequent degradation processes (Gundelwein et al., 2007).Assuming constant photosynthetic isotope fractionation in source plants in the region (C 3 plants are ubiquitous in the Arctic; Tieszen, 1973), we use δ 13 C ratios as a degradation proxy.According to Heyer et al. (1976), decomposition discriminates against the lighter isotope ( 12 C), resulting in more negative δ 13 C ratios.Thus, this proxy is used in the following way: lower (more negative) δ 13 C values are connected to less degraded material, while higher (less negative) δ 13 C values reflect greater decomposition.Ages were determined by radiocarbon dating of selected macroscopic plant remains performed at the Poznań Radiocarbon Laboratory, Poland (Goslar et al., 2004).The presented radiocarbon ages are uncalibrated ages; Table 1 includes calibrated ages as well.Radiocarbon ages are given in year before present (yr BP).

Lipid biomarkers
To look more closely at the molecular composition, we used specific lipid biomarkers.Molecular fossils or biomarkers were studied by chromatography methods coupled with mass spectrometers.Characteristic fractions like n-alkanes, n-fatty acids, sterols, and hopanes were isolated.Because the TOC wt % in the profiles is not equally distributed, we calculated and visualized the biomarker concentration as µg gTOC −1 wt % and µg gSediment −1 (µg gSed −1 ).For the results, we focus on µg gTOC −1 wt % .Every radiocarbon-dated sample and additional samples were used for biomarker analysis.In total 25 biomarker samples were analyzed.Independently of TOC wt % , the sample selection for biomarkers was based on their stratigraphic position with the aim of covering the maximum time period.

Extraction and fraction separation
For lipid biomarker analyses, 2-12 g of ground sediment was weighed in an extraction cell with an accelerated solvent extractor (ASE 200,Dionex).Samples were extracted with dichloromethane / methanol (99 : 1).Each sample was held in a static phase for 20 min at 75 • C (after 5 min heating, no preheating) at a pressure of 5 MPa.Afterwards, the dissolved compounds were concentrated with a TurboVap (Zymark) closed-cell concentrator and further dried by evaporating the solvent in a stream of nitrogen gas.After that, internal standards (5α-androstane for the aliphatic fraction, ethylpyrene for the aromatic fraction, 5α-androstan-17-on for nitrogen-, sulfur-, and oxygen-(NSO-) containing compounds and erucic acid for the NSO fatty-acid fraction) were added.The amount of internal standards depended on the TOC wt % content (< 10 wt%: 8 µg; > 10 to ≤ 25 wt%: 20 µg; > 25 wt%: 50 µg).After the removal of the n-hexane-insoluble fraction (by the addition of a large excess of n-hexane, called "asphaltene" precipitation), the hexane-soluble portion of the extract was separated by medium-pressure liquid chromatography (MPLC; Radke et al., 1980) into fractions of different polarity (aliphatic and aromatic hydrocarbons as well as polar hetero-(NSO) components).Afterwards, the NSO fraction was split into a fatty acids and an alcohol fraction using a KOH-impregnated silica gel column (Schulte et al., 2000).For this study, the focus was placed on the aliphatic fraction (containing n-alkanes and triterpenoid compounds) and the NSO fraction (containing n-fatty acids).The fractions were measured by gas chromatography-mass spectrometry (GC-MS).All compounds of interest were identified using the Xcalibur software (Thermo Fisher Scientific).

GC-MS measurement and compound quantification
The n-alkanes, n-alcohols, hopenes (hop-17(21)-ene), and other triterpenoids (β-amyrin (olean-12-en-3β-ol), olean-12-ene, and olean-13(18)-ene) were measured with a GC-MS system (GC -Trace GC Ultra; MS -DSQ; both Thermo Fisher Scientific).Prior to the measurements, the n-fatty acids were methylated with diazomethane and the alcohols were silylated with N-methyl-N -trimethylsilyltrifluoroacetamide (MSTFA).The GC was equipped with a programmable temperature vaporization (PTV) injector system (starting temperature of 50 • C, heating rate of 10 • C s −1 to 300 • C, isothermal holding time of 10 min, operated in splitless mode) and a fused silica capillary column (SGE BPX5, 50 m length, 0.22 mm inner diameter, 0.25 µm film thickness).For the measurements the GC oven was programmed with a starting temperature of 50 • C, a heating rate of 3 • C min −1 to 310 • C, and an isothermal holding time of 30 min.Helium with a constant flow rate of 1 mL min −1 was used as a carrier gas.For the n-fatty-acid fraction a different temperature program (starting temperature of 50 • C, 1 min isotherm, heating rate of 3 • C min −1 to 350 • C, isothermal holding time 25 min) was used.For compound identification, the gas chromatograph was linked to a mass spectrometer, which was operated in electron impact ionization mode at 70 eV.The temperature of the ion source was set to 230 • C. Full scan mass spectra were recorded from m/z 50 to 600 Da at a scan rate of 2.5 scans s −1 .For the nfatty-acid fraction, the scan rate was m/z 50 to 650 Da.
The quantification of n-alkanes, n-fatty acids, and βamyrin was done in the GC-MS total ion current chromatogram by relating the peak area of the target compound to the peak area of an internal standard of known concentration.Other triterpenoids like olean-12-ene, olean-13(18)-ene, and hopene were quantified using the m/z 191 mass trace relative to the peak area of the β, β-diploptene (in the m/z 191 mass trace), the concentration of which was calculated in the to-tal ion current chromatogram relative to the internal standard (5α-androstane).

Absolute lipid concentration
The absolute lipid concentration is used as a rough estimator of organic-matter quality for degradation in the following way: the higher the concentration, the better the conservation of the lipid and the better the quality of the organic matter.

Carbon preference index
The CPI was introduced by Bray and Evans (1961) as the ratio of odd-to neighboring even-numbered alkanes; this is a measure of the alteration of organic matter.Here we use the improved formula from Marzi et al. (1993).In addition, we also applied the CPI for fatty acids in which evennumbered fatty acids predominate over adjacent odd n-fatty acids (Glombitza et al., 2009).
where n is the starting dominating chain length/2, m is the ending dominating chain length/2, i is the index (carbon number), and C is concentration.
The CPI is used as a degradation and alteration proxy by quantifying the odd over even (n-alkanes, Fig. S2) or even over odd (n-fatty acids, Fig. S3) predominance of the carbon chains (Bray and Evans, 1961;Glombitza et al., 2009).A low CPI means mature or degraded organic matter (e.g., CPI of crude oil ∼ 1).

Average chain length
As introduced by Poynter (1989), the n-alkane ACL value is the concentration-weighted mean of different carbon chain lengths in a geological sample.For n-alkanes we use the C 23 -C 33 interval, for n-fatty acids the C 20 -C 34 : where i is the index (carbon number) and C is concentration.
The ACL is a rough OC source parameter.A schematic showing different chain lengths in different organisms is given in Fig. S4.The higher C 3 land plants are expected to have an ACL of ∼ 28-29.

Higher-plant fatty-acid index
The ratio of the major even wax alcohols over the sum of major odd wax alkanes plus even alcohols was introduced by Poynter (1989) as the higher-plant alcohol (HPA) index.It is applied as an indicator for chemical degradation of the wax components.Based on this index, but using fatty acids instead of alcohols, we developed the HPFA index.The general assumption for this index is that it reflects the preservation degree of organic matter due to the higher lability of n-fatty acids in relation to n-alkanes.The HPFA ratio cannot be considered an absolute index of degradation but is an indicator of the relative amounts of the more labile fatty acids that remain in a sample.Since n-alkanes are preserved preferentially compared to n-fatty acids, a decrease in this index indicates increased decomposition (the more degraded, the lower the HPFA index).

Acetate
Pore water was obtained from each sample by centrifugation in specific pore water tubes.Water extracts were analyzed twice using ion chromatography with conductivity detection (ICS 3000, Dionex).An analytical column (AS 11 HC, 2 × 250 mm, Dionex) was used at a constant 35 • C. The sample was eluted with a KOH solution of varying concentration over time.The initial concentration was 1.4 mM.Between 0 and 6 min, the KOH solution was increased at a constant rate to 1.6 mM.Between 6 and 12 min the solution was increased to 10.0 mM KOH, and a concentration of 15.0 mM KOH was reached at 22 min.After 32 min, 60.0 mM KOH concentration was achieved, maintained for 1 min, and followed by a rapid decrease to 1.4 mM after 33 min, when samples were fixed for 45 min to equilibrate the system.For the quantification of acetate, standards containing the investigated compound were measured.The standard deviation of the sample and of standard quantification was < 5 %.Because acetate can act as excellent feedstock for microbes (Smith and Mah, 1980;Vieth et al., 2008) and it has been shown that acetate was rapidly consumed in the presence of oxygen and nitrate (Kuesel and Drake, 1995), we use the acetate pore water concentrations in the different deposits as a parameter to assess the quality of the organic matter and to compare the potential of the different deposits for future microbial degradation.

Significance testing
For testing the samples concerning their statistical distribution, the Shapiro-Wilk normality test was applied.
Because of non-normal distribution, we used the Mann-Whitney-Wilcoxon test for significance testing of Yedoma and thermokarst samples.In order to compare the five different profiles, we used the Kruskal-Wallis rank sum test.

Principal component analysis
Multivariate statistical techniques, such as the PCA used here, allow the analysis of multiple variables in order to investigate connections between the different degradation proxies.Prior to the PCA, concentration data were transformed using a log (x+1) transformation.As the square root transformation is commonly applied to count data, especially if the values are mostly rather small, we decided to use this weaker (compared to logarithm) transformation for the TOC ( wt % and kg m −3 ) data.Both transformations were applied to reduce right skewness and to put the parameters on the same scale.We performed three PCA runs.First, a PCA of the sediment parameters was implemented to infer differences between Yedoma and thermokarst deposits.Second, a PCA of biomarker proxies was performed.For this purpose, other characteristics were added as supplementary variables (TOC wt % , TOC kg m −3 , C / N, δ 13 C, grain size, BD, ice content, and depth) without being included in the PCA calculation.These supplementary variables have no influence on the PCA and were plotted afterwards in the PCA biplot.Third, a PCA was conducted on samples of the major odd n-alkanes to infer possible changes in the source organisms with the same supplementary variables as described above.Computations were performed using the "vegan" package of R software (Oksanen, 2013).

Results
Stratigraphically, there are two types of deposition units at the study site.The first unit is composed of Yedoma deposits.The second unit represents thermokarst deposits re-sulting from thermal degradation of Yedoma.Grain-size distributions (Figs. 2, S1) and PCA of sediments illustrate that thermokarst deposits are made up of degraded Yedoma sediments.According to Gubin and Veremeeva (2010) and Zanina et al. (2011), the Yedoma deposit soil types are mainly less developed cryopedoliths containing more developed paleocryosol parts (Figs. 3 and 4, labeled and grey shaded areas).
3.1 Organic-matter quality of Yedoma deposits

Sedimentological and biogeochemical proxies
The radiocarbon ages (Table 1, Fig. S5) of the Yedoma deposits range from infinite ages (> 55 000 yr BP) at the very bottom to 30 100 yr BP at the uppermost sampled Yedoma unit.This is comparable to other Yedoma sequences in the region (Schirrmeister et al., 2011b).The mean grain sizes show a decreasing trend in the Buo-04 lower Yedoma profile, from 28 µm at the bottom to 11 µm in the upper part of Buo-04-A.The Buo-02 Yedoma profile shows no trend but exhibits a more heterogeneous mean grain size, including three maxima at 22.5 m a.s.l.(32 µm), 23.7 m a.s.l.(34 µm), and 25.5 m a.s.l.(33 µm).Nevertheless, all Yedoma deposit samples are classified as poorly sorted medium-to-coarse silts with a stable low clay fraction (< 15 %).
The TOC wt % contents vary from 0.2 wt% at 5 m a.s.l to 24.0 wt% in a peaty paleocryosol horizon at 24 m a.s.l.(Fig. 3).The mean TOC wt % content is 2.4 wt% (median 0.97 wt%).Calculating the TOC kg m −3 according to Strauss et al. (2013), by utilizing the BD (between 0.1 and 1.5 10 3 kg m −3 ; 10 3 kg m −3 = g cm −3 ) and ice content (without ice wedges; 21 to 90 vol %), the Yedoma sediments contain from 3 to 46 kg C m −3 with a mean of 14 kg C m −3 (median 9 kg C m −3 ).The maxima correspond to the peaty horizons with large TOC wt % contents and a low BD.Within the paleocryosol horizons, located at 6.8, 24.0 to 24.5, 24.8, and 27.8 to 28.9 m a.s.l., maxima in the C / N ratio are observable.The C / N range in these horizons is 8 to 31.In the cryopedolith profile parts, the C / N maximum is reached in the lowermost Buo-04-C sub-profile (17.7 and 16.7).The C / N of the rest of the Yedoma profile falls between 4.1 (at Buo-02-C, 23.7 m a.s.l.) and 14.3 (below the paleocryosol at 23.5 m a.s.l.) The δ 13 C of the Yedoma deposits ranges from −29.0 to −24.7 ‰.The minima fit the maxima of the C / N ratio well in the paleocryosol horizons at 6.8, 24.0 to 24.5, 24.8, and 27.8 to 28.9 m a.s.l.The minimum C / N of the Buo-02-C sub-profile corresponds approximately to the δ 13 C maximum (−25.0 to −24.7 ‰).

Biomarker proxies and indices
A series of long-chain n-alkanes that exhibit a strong oddcarbon preference ranging from n-C 21 to n-C 33 are iden-  tified in all Yedoma samples (Fig. S2).Moreover, the nalkanes show a unimodal distribution maximizing at the C 27 , C 29 , or C 31 n−alkane (Fig. S2).The n-fatty acids show strong even over odd carbon number predominance and a bimodal distribution ranging from C 14 to C 30 (Fig. S3).The maxima are generally located at n-C 16 in the lower carbon number range and at n-C 24 in the higher carbon number range.Total n-alkane and n-fatty-acid concentrations related to TOC wt % and sediment weight show a homogeneous pattern similar to that of the TOC wt % and C / N values.
The n-alkane concentration ranges from 3 to 75 µg gSed −1 (mean 20 µg gSed −1 ) and from 387 to 1715 µg TOC −1 wt % (mean 1132 µg TOC −1 wt % ).The n-fatty acids range from 4 to 306 µg gSed −1 (mean 51 µg gSed −1 ) and from 475 to 4669 µg TOC −1 (mean 2196 µg TOC −1 wt % ).This Yedoma series shows a distinct preference between even and odd carbon.The mean CPI values of the n-alkanes (12.2, ranging from 8.3 to 15.9) are higher than the CPI values of the n-fatty acids (4.9, ranging from 3.8 to 7.6).Because n-fatty acids are functional compounds (including a functional group, e.g., a carboxyl group), their degradation rates are much higher compared to those of n-alkanes (Poynter and Eglinton, 1990).This statement is also based on the assumption of similar sources.and n-fatty acids is very stable at around 28.4 (range: 27.6 to 29.2) and 25.0 (range: 23.8 to 25.6), respectively.Relatively higher hop-17(21)-ene concentrations are used as an indicator for a lower organic-matter degradation state.In the lower Yedoma profile, the hop-17(21)-ene ranges from 0.0 µg gTOC −1 in the lowermost and uppermost samples (4.3 and 18.5 m a.s.l.) to the overall maximum of 19.4 µg gTOC −1 at the Buor-Khaya site (9.1 m a.s.l.).At Buo-02, the hop-17(21)-ene concentration is lower compared to the other Yedoma profile, with a mean of 1.9 µg gTOC −1 wt % and a maximum of 7.7 µg gTOC −1 wt % in the potentially Holocenecontaminated uppermost sample.The HPFA ratio for the Yedoma samples is very stable around the mean value of 0.50 (median 0.54), with a minimum in 18.5 m a.s.l.(0.15) and a maximum at the uppermost sample (0.69) at 29.7 m a.s.l.For Yedoma, the oleanene ratio is 0.0 (except for a ratio of 10.0 in the uppermost sample).The acetate content of the Yedoma sample is between 0.6 and 57.5 mg L −1 , with a mean of 6.7 mg L −1 (median 1.2 mg L −1 ).

Sedimentological and biogeochemical proxies
The radiocarbon dating shows Holocene ages between 8140 ± 50 and 3665 ± 35 yr BP (Fig. S6, Table 1).The lowermost Buo-05-C profile shows an age inversion for the two samples (0.3 and 2.2 m a.s.l.).The mean grain size at Buo-05 from the bottom to 6.7 m a.s.l. is 13 µm.Above, the mean grain size increases to 19 µm.The Buo-05 clay fraction is stable at a low level (< 15 %).The Buo-01 profile shows a very scattered grain size ranging from 4 to 44 µm mean grain size.For the whole data set, there is a maximum in the clay fraction (35 %) in the peat horizon at 8.7 m a.s.l.Buo-03 shows a slight decrease from 18 to 11 µm.All thermokarst deposits are classified as (very) poorly sorted silts.Similarly to the Yedoma deposits, the BD of the thermokarst deposits is between 0.1 and 1.5 × 10 3 kg m −3 and the ice content (without the ice wedges) is 23 to 87 wt% (Fig. 4).
The mean TOC wt % contents of the thermokarst deposits, 4.7 wt% (median 1.7 wt%), are higher compared to Yedoma deposits, varying between 0.2 and 43.0 wt%.Minimum and maximum TOC wt % both occur at Buo-01 and exhibit the same scatter as in the grain sizes.TOC kg m −3 ranges between 2.8 and 93.5 kg C m −3 (mean 24, median 19 kg C m −3 ).
At Buo-05 the C / N ratio is stable at around 9 to 10 (Fig. 4), except for a paleocryosol horizon at 9.2 m a.s.l. that shows a value of 22.At Buo-01, the C / N ratio below the paleocryosol horizon is remarkably low, between 2 and 9, followed by the overall maximum in the peaty horizon with a ratio of 34.The Buo-03 cryopedolith samples show C / N ratios around 10, while the paleocryosol samples exhibit C / N ratios from 16 to 19.The δ 13 C values range between −29.5 and −25.0 ‰, with minima corresponding to the C / N max- ima at the paleocryosol horizons (anticorrelated to the C / N, Fig. 5a).
A series of long-chain n-alkanes were recognized in all thermokarst samples with a strong odd carbon number preference ranging from n-C 21 to n-C 33 .Nearly all samples show a unimodal distribution of n-alkanes maximized at C 27 , C 29 , or C 31 (Fig. S2).Sample Buo-03-A-03 alone does not fit into this scheme because it maximizes at n-C 25 .Compared to Yedoma samples, the short-chain fraction < n-C 27 is more pronounced (Fig. S2).The n-fatty acids show a strong evencarbon-number preference and a bimodal distribution between n-C 14 and n-C 30 (Fig. S3), but the n-C 16 is less pronounced than in the Yedoma deposits.An exception to this is found in sample Buo-01-A-02, where the C 16 monomer reaches the overall maximum of the distribution.Apart from that, the maxima are generally located at the C 24 n-fatty acid.
The n-alkane CPI of thermokarst averages 9.6 (median 9.3) and is lower compared to the Yedoma deposits, although the CPI values are in the same range (between 7.0 and 15.3).The CPI of the fatty acids ranges from 4.0 to 9.0 (mean 5.3, median 4.9).The ACL of n-alkanes and fatty acids reveal a homogeneous signal between 27.2 and 29.2 (mean 28.3) for n-alkanes and 23.6 to 25.6 (mean 24.8) for n-fatty acids.
The acetate content of the thermokarst samples is between 0.4 and 109.4 mg L −1 , with a mean of 23.5 mg L −1 (median 2.8 mg L −1 ).Large acetate contents are found especially in the middle part of Buo-05, from 3.4 m a.s.l.(74.1 mg L −1 ) to 6.1 m a.s.l.(109.4 mg L −1 ), and in the uppermost Buo-03 sample (35.3 mg L −1 ).

Significance testing
Except for the Yedoma CPI, the Yedoma HPFA and the thermokarst hop-17( 21)-ene Shapiro-Wilk normality tests reveal a non-normal distribution.Based on this, we chose non-parametric significance testing.This reveals significant differences for TOC, C / N, δ 13 C, and HPFA on the stratigraphical level (Yedoma vs. thermokarst, Mann-Whitney-Wilcoxon test, Table S1 in the Supplement).On the profile level, we found significant differences for C / N, δ 13 C, and HPFA by applying the Kruskal-Wallis test (Table S1).

Principal component analyses
The first PCA diagram (Fig. 6a) shows that thermokarst sediments, especially at Buo-05, could not be separated from Yedoma deposits.This diagram, including the first two principal components, explains 79 % (pc1 57 %, pc2 22 %) of the total data set variance.The second PCA diagram (Fig. 6b) illustrates that biomarker quality estimators in Yedoma samples have slightly lower variability because they cluster in an area at pc1 and pc2 > 0, while the thermokarst samples do not cluster.In this diagram 53 % of the data set variance is explained.Moreover, this PCA shows that there is good consistency between the CPI alkane quality estimator and the C / N ratio (Fig. 6b).The PCA of the n-alkane chain length (Fig. 6c) shows that the best separating variables for thermokarst are the shorter-chain n-alkanes (C 17 , C 19 , and C 21 ), contrary to C 29 for the majority of the Yedoma samples.The pc1 explains 39 and pc2 explains 29 % (total 68 %) of the data set variance.

Discussion
The Buor-Khaya Peninsula is a typical Yedoma hillthermokarst basin landscape of the Yedoma region (Strauss et al., 2013).The Yedoma deposits cover ∼ 15 % of the peninsula (Günther et al., 2013), which is less than the Yedoma region mean of 30 %, but within the overall range of Yedoma deposit coverage (Grosse et al., 2013;Strauss et al., 2013).Thus, the current study of Yedoma and thermokarst deposits is representative of an area covered by similar permafrost deposits of late Pleistocene and Holocene age.

Sediment facies
The grain-size distribution curves (Figs. 2, S1) indicate a constant deposition environment for the Yedoma sequences.According to Strauss et al. (2012), there have been stable deposition conditions during Yedoma accumulation; this hypothesis is supported by the data presented here.The three thermokarst profiles include three different kinds of thermokarst deposits.Buo-05 is dominated by a lake facies containing valves of two freshwater ostracod taxa: Cytherissa lacustris and Cypria sp.Moreover, shells have been found in Buo-05 (Strauss and Schirrmeister, 2011).An ice wedge is located next to Buo-01, which points to subaerial conditions, such as a polygon mire.Buo-03 is interpreted as initial thermokarst on top of a Yedoma hill.Thus, the grain-size distributions of Buo-05 and Buo-01 reveal that the thermokarst is granulometrically composed of the same material as Yedoma.The grain-size distributions in Buo-03 paint a different picture.This distribution is likely caused by the early state of thermokarst development dominated by peat aggradation.This peat can act like a selective sediment trap influencing the grain-size distributions, e.g., by producing a less distinct coarse-silt-fine-sand peak.

Organic-matter degradation
The organic-matter proxies of Yedoma deposits are less variable than those of thermokarst deposits (Buo-01 and 03).Except for the paleocryosols, the cryopedolith parts of the Yedoma and the Buo-05 thermokarst profile reveal a rather homogeneous picture (Figs. 3,4,S5,S6).Constant grainsize distributions, less TOC wt % , and smaller absolute lipid concentration scattering reveal that the OC stored in the Yedoma deposits has likely been kept perennially frozen since incorporation.The organic-matter signatures (Figs. 4,S2,and S3) as well as the grain-size distributions (Figs. 2, S1) of thermokarst deposits, especially in Buo-01 and Buo-03, show broader variations.This is caused by a more complex degradation and redeposition history due to reworking.The degradation markers of organic matter found in the paleocryosol parts of all profiles reveal a less degraded state, indicating that the organic matter in these portions is the best preserved.
The mean TOC wt % content for Yedoma deposits is comparable to other sites (Fig. S7) in the Yedoma region (Schirrmeister et al., 2011b(Schirrmeister et al., , 2013)).Intense accumulation and frozen preservation of plant remains (14 kgC m −3 for Yedoma and 24 kgC m −3 for thermokarst deposits) is caused by syngenetic permafrost formation in polygonal tundra landscapes over long periods in the Quaternary (Schirrmeister et al., 2013).However, comparing the studied deposits to the overall Yedoma region mean (19 kgC m −3 for Yedoma deposits and 33 kgC m −3 (disregarding wedge-ice content) for thermokarst deposits; Strauss et al., 2013) on the Buor-Khaya Peninsula reveals that both deposit types contain less OC.Nevertheless, these numbers show that such deposits comprise a large pool of dormant carbon, which could be reactivated due to permafrost thawing.Moreover, thermokarst deposits seem to be the quantitatively more important OC pool (Yedoma : thermokarst carbon ratio ∼ 2 : 3).The higher carbon inventory in thermokarst deposits is partially related to a concentration effect for reworked Yedoma OC due to thaw subsidence progression, including ground ice loss plus the input of Holocene OC.Together with ecosystem recovery, thermokarst basins can act as a local sink for portions of the carbon released from thawing permafrost deposits (van Huissteden and Dolman, 2012).Nevertheless, at the same time, thermokarst lakes also promote intense organic-matter degradation, including methane production in the anaerobic environments of organically rich lake sediments and unfrozen deposits (Walter et al., 2007b;Shirokova et al., 2013).To answer the question raised in the introduction of whether the thermokarst organic-matter pool is as degradable as the frozen late Pleistocene Yedoma, we visualized the stratigraphically differentiated main proxies in Fig. 7.
In our study the C / N data show some overlap (Fig. 7b).The average values are relatively close together for all profiles, but the differences are statistically significant (Table S1).Thus, the C / N medians and means suggest a lower degradation state or better organic-matter quality of thermokarst deposits (especially Buo-05 and Buo-03).Moreover, in both Yedoma and thermokarst deposits, the same pattern is visible: a positive linear relationship exists between TOC wt % and C / N ratios (Fig. 5b).In soil science literature, it is agreed that the elemental composition of organic matter is affected by the degree of humification and microbial activity that metabolize the organic matter (Kumada, 1987).Ongoing organic-matter decomposition will release stored C to the atmosphere and N to the soil (Weintraub and Schimel, 2005), resulting in a lower C / N ratio for more degraded deposits (Gundelwein et al., 2007).This was found in (sub-) Arctic peat deposits and soils, where the C / N ratio decreases with depth (Kuhry and Vitt, 1996;McKane et al., 1997;Ping et al., 1998).Because a high TN content can promote the stabilization of organic matter at late stages of decomposition (Berg, 2000), this further supports the interpretation that a low C / N ratio indicates recalcitrant or matured organic matter (Rumpel and Kögel-Knabner, 2011).Schädel et al. (2014) found, in incubation studies, that the C / N ratio is a good estimator for organic-matter decomposability and vulnerability.Although the C / N ratios are lower than in Arctic peat deposits (Hugelius et al., 2012;Routh et al., 2014), the ratios are still in the range of or higher than those found in many other deep mineral soils of the temperate zone (Jenkinson et al., 2008;Rumpel and Kögel-Knabner, 2011).Thus, both Yedoma and thermokarst deposits show relatively good organic-matter quality for microbial decay after becoming available by thaw.The C / N ratios, especially for the paleocryosols, suggest that good quality organic matter was preserved (by temperature near or below 0 • C during thermokarst processes) for future decomposition.This is shown by the δ 13 C ratio as well.Neglecting the influence of different sources of organic matter on the δ 13 C ratio, which is justified by constant ACL values of > 28 (higher land plants, Fig. S4) for Yedoma and thermokarst deposits, the δ 13 C ratio is an appropriate proxy for estimating the intrinsic state of degradation.Therefore, the δ 13 C values indicate a significantly lower organic-matter degradation for the thermokarst samples, implying a better quality than that found in Yedoma samples.The high CPI values of the thermokarst and the Yedoma organic matter (around 9 and higher) indicate fresh and less degraded terrigenous organic matter (Brassell et al., 1978) for both deposits (Fig. 7d).The CPI values of the Yedoma deposit organic matter, which are (significantly) higher than in the thermokarst deposits, indicate a better quality for further decomposition (Fig. 7d).Routh et al. (2014) states that other, more labile compounds like n-alcohols and n-fatty acids are degraded to n- alkanes.Thus, an increase in n-alkanes (Figs. 3 and 4, absolute lipid concentration column) is an indicator for cumulative decay.We do not see a decreasing trend, which points to a constantly low decomposition state.In addition, no increasing n-fatty-acid CPI with depth (as was shown in an Arctic peat by Andersson and Meyers, 2012) was obvious.Andersson and Meyers (2012) interpreted this to reflect fatty-acid production during humification, but we do not see this humification effect in our data, neither in Yedoma nor in the thermokarst deposits.Moreover, as indicated by the dominance of long-chain n-alkane and n-fatty-acid compounds vs. compounds of shorter chain length (Höfle et al., 2013), we confirm the interpretation that there is good organicmatter preservation in both Yedoma and thermokarst deposits.At first view, the hop-17(21)-ene (Fig. 7e) concentration does not show a significant preservation difference between both kinds of deposits because the Buo-04 Yedoma profile contains hopene concentrations in the same range as those found in thermokarst deposits.However, if we focus on the median values, the Yedoma deposits again appear to be slightly more degraded than the thermokarst deposits.With the exception of Buo-01, the HPFA index (Fig. 7f) also suggests lower degradation and better organic-matter quality in the thermokarst deposit profiles (Buo-05 and Buo-03).Our HPFA index, introduced based on the HPA index of Poynter (1989), which was tested in the Arctic environment by Routh et al. (2014), is an appropriate indicator of the relative amount of the labile fatty acids that remain in a sample.The uppermost samples just below the surface at Buo-04, Buo-05, and Buo-01 and with lower HPFA values are clearly an exception and suggest the entrainment of higher proportions of material influenced by Holocene degradation.This is likely caused by the fairly recent influence of an active layer or transient layer and warmer permafrost temperatures.The oleanene ratio shows a separation of Yedoma and thermokarst deposits, but this ratio is dominated by numerous 0.0 measurements in the Yedoma deposits.These results might be caused not only by the transformation of β-amyrin to olean-12-ene (by losing the hydroxyl group) or to olean-13(18)-ene (by losing the double bond) but also by processes so far unknown in the Yedoma deposits.Thus, because of sparse data and to avoid overinterpretation, this proxy is interpreted only as a weak indication of better Yedoma organicmatter quality for further decomposition.Summing up Fig. 7, thermokarst organic matter is partly less degraded compared to the organic matter sequestered in Yedoma deposits (see Table S1, significance for C / N, δ 13 C, and the HPFA index).The CPI points in the other direction (Fig. 7 and Table S1).For hop-17(21)-ene, we do not see significant differences.Nevertheless, the interquartile ranges show an overlap for most proxies.We interpret this in the following way: compared to unaltered Yedoma deposits, degradation during thermokarst processes, but also heightened amounts of OC input during climatically more favorable Holocene times, balance each other out with regard to the organic-matter quality for future degradation.Nevertheless, as there is more carbon stored in the thermokarst basins (Strauss et al., 2013), thermokarst deposits imply a higher intrinsic potential to contribute greenhouse gases in a warmer future.This is supported by the acetate data indicating a higher mean content for the thermokarst deposits.Acetate is an excellent substrate for microbial turnover, e.g., acetoclastic methanogenesis (Kotsyurbenko et al., 2004).The PCA confirms this picture that shows little difference between the organic-matter preservation of the Yedoma and the thermokarst samples.Especially Fig. 6a, supported by Fig. 2, reveals that Yedoma and thermokarst deposits are composed of similar sediments.The Buo-05 thermokarst profile is very similar to both Yedoma profiles.The PCA of the degradation proxies (Fig. 6b) also shows no clusters but exhibits slightly better separation between both kinds of deposits.Figure 6b reveals that the C / N ratio, the δ 13 C ratio, and the CPI are correlated.This is also separately illustrated in Fig. 5a and  c.Thus, these proxies seem to confirm each other.The PCA of the n-alkane chain length points to a potential dominance of longer-chain alkanes in Yedoma and shorter-chain alkanes in thermokarst, indicating better quality for further decomposition of Yedoma samples (Höfle et al., 2013).Exceptions are the Buo-05-A-01 and Buo-03-A-03 thermokarst samples, which point in the same direction as the n-C 35 concentration.
The abovementioned overlap of the interquartile range (Fig. 7) and especially the PCA of the biomarkers (Fig. 6b  and c) show that the organic-matter degradation and decomposition vulnerability is heterogeneous and depends on different decomposition trajectories and differing former decomposition and incorporation histories.This is likely shown in both Yedoma and thermokarst deposits, covering the whole range of degradation proxy values (Fig. 7b, c, e).To elucidate this was one of the benefits of the applied multiproxy approach.With the addition of biomarker data, it is possible to show that the permafrost organic-matter degradation is not a linear function of age or sediment facies but likely a combination of (interrupted) degradation cycles and a cascade of degradation steps.In particular, the reasonably good organic-matter preservation of thermokarst deposits reveals that the sediment degradation processes do not necessarily degrade the organic matter.Potentially, the loss of labile OC during thermokarst processes was compensated for by high rates of Holocene OC accumulation in, e.g., lake sediments.Nutrient release from thawing permafrost could have stimulated lake productivity, whereas decomposition was slow because of low lake temperatures, resulting in cold anoxic lake environments (Boike et al., 2013;Walter Anthony et al., 2014).When the lake drained, permafrost formation rapidly recovered the sediments (Jones et al., 2011), including any possibly newly accumulated OC.

Fate of organic matter
The permafrost OC resilience or vulnerability is a topic of recent research (Schuur and Abbott, 2011;Knoblauch et al., 2013;Hodgkins et al., 2014;Li et al., 2014;Mu et al., 2014).Any warming permafrost is potentially vulnerable to thawing.The remaining important question is this: what is the fate of the organic matter exposed to degradation after permafrost has thawed?The lipid biomarker data discussed (CPI etc.) indicate that the organic matter in the sediments was, after initial degradation processes, relatively quickly protected against microbial alteration by freezing.This is confirmed by an absent degradation-depth trend which reveals good organic-matter quality independent of age.Thus, the very old frozen organic matter is also vulnerable to degradation after thawing.This interpretation fits results from studies of permafrost-affected Arctic peats (Hugelius et al., 2012;Routh et al., 2014).Walter Anthony et al. (2014) found a net accumulation in thermokarst basins since the last deglaciation but predict that these will change to a large carbon source when permafrost thaws and the OC will be available for oxidation.Due to ongoing climate warming in the Arctic, Grosse et al. (2011b) supposed an increasing occurrence in and magnitude of disturbance processes, especially fire and thermokarst, which will accelerate permafrost degradation.Because our sedimentological and biomarker proxies show a low degradation state, especially for the paleocryosol sequences, we expect a significant vulnerability to microbial degradation after thawing.As evidence that the OC is vulnerable when thawed, Gaglioti et al. (2014) found that ∼ 10 times more ancient OC found in permafrost was made available for degradation during warm times of the Holocene (Holocene Thermal Maximum (11.7-9.0 ka) and Bølling-Allerød periods) than is available today.Through increased disturbances, such as deep surface subsidence caused by thawing and the draining of excess water from melting ice in a warmer climate, the Yedoma and, to a lesser degree because of lower excess ice content, the thermokarst organic matter could become bioavailable deep in the sediment.The wedge-ice volume is estimated to be up to ∼ 60 vol % for Yedoma and up to ∼ 10 vol % for thermokarst deposits (Ulrich et al., 2014).When added to segregated ice, ∼ 80 and ∼ 65 vol % mean sedimentary ice volume exists in Yedoma and thermokarst, respectively (Strauss et al., 2013).When it becomes available and is exported as dissolved OC to, e.g., river systems, Vonk et al. (2013) and Mann et al. (2014) found that dissolved OC (< 0.45 µm) in ancient Yedoma is exceptionally biolabile.However, if it is not dissolved, the suspended (> 0.45 µm) eroded ancient organic matter could be protected from extensive degradation by organo-mineral bonds, which stabilize the organic matter (Höfle et al., 2013) and, in an aquatic environment, promote rapid settling because they weigh down the organic matter (Vonk et al., 2010).
From the modeling perspective, global-scale models have, so far, been limited because they implement one-dimensional vertical thaw only (Koven et al., 2011;Schneider von Deimling et al., 2012;Schaphoff et al., 2013).Thus, the potentially labile Yedoma and thermokarst deep OC pool described in this study is not realistically implemented in these models because the models disregard rapid phenomena, such as thermokarst processes.Thermokarst processes, despite being local in nature, are widespread on the regional scale (Grosse et al., 2011a) and may constitute the crucial process making the deep OC studied here microbiologically available.

Conclusions
Being freeze-locked, a great amount of organic matter in the studied sediments is highly decomposable.Generally, in all applied proxies there is no degradation-depth trend obvious, revealing that permafrost acts like a freezer, preserving the organic matter after freezing.Based on interpreting the mean values of the C / N ratio, isotope ratio (δ 13 C), and the HPFA index, the thermokarst organic matter is less degraded and of better quality for degradation after thawing compared to the organic matter sequestered in Yedoma deposits.The CPI data suggest less degradation of the organic matter from both deposits, with a higher value for Yedoma organic matter.For the hop-17(21)-ene concentration, no significant difference was found.We do not see any conflict between these two determinations because the interquartile ranges overlap for most proxies.We interpret this as indicating a comparable magnitude of organic-matter quality in both kinds of deposits, but with a likely better thermokarst organic-matter quality for further degradation.For a modeling approach, this conclusion could be extrapolated to the Laptev Sea region as the studied deposits are akin to other Yedoma and thermokarst deposits of the northeast Siberian Arctic (Schirrmeister et al. 2011a).
The fate of mobilized Yedoma deposit OC depends largely on the environmental conditions that exist during the thermokarst processes and in the resulting thermokarst basin.In conclusion, when the conditions during the thermokarst processes are good for organic-matter preservation, for example cold (slightly above 0 • C) or anoxic (lake) conditions, and when reworked fossil organic matter can rapidly refreeze to permafrost, good-quality organic matter for further decomposition can be maintained and inputs likely compensate for losses due to thermokarst degradation.
The Supplement related to this article is available online at doi:10.5194/bg-12-2227-2015-supplement.
analysis and interpretation.J. Strauss carried out the laboratory analyses, except for one profile, which was analyzed by L. Eichhorn.U. Herzschuh designed the statistical analyses.J. Strauss planned and wrote the publication with input from all co-authors.

Figure 1 .
Figure 1.Location of the Buor-Khaya Peninsula and the study area.The square black box in the globe inset indicates the area shown in the map below.The profile diagram and the photographs below it show the profiles and their positions relative to each other.Adapted from Strauss and Schirrmeister (2011); pictures taken by J. Strauss.

Figure 2 .
Figure 2. Three-dimensional grain-size distributions of (a) Yedoma and (b) thermokarst profiles.To avoid an overlap of Buo-05 and Buo-01 in (b), the altitude axis was adapted and does not ascend consistently.A two-dimensional grain-size plot is shown in Fig. S1.

Figure 3 .
Figure 3. Summary of sedimentological, biogeochemical, and biomarker parameters for the Buo-04 and Buo-02 Yedoma profiles.All diagrams are drawn in such a way as to show more degraded samples on the left and less degraded samples on the right side.Thus, the axis of δ 13 C shows descending values.In the text, the paleocryosol parts are reported with altitude measurements from the lowest to the highest sample of each paleocryosol.The grey shaded areas are for visualization, not for exact height estimations of the paleocryosols.

Figure 4 .
Figure 4. Summary of sedimentological, biogeochemical, and biomarker parameters for the Buo-05, Buo-01, and Buo-03 thermokarst profiles.The n-alkane and n-fatty-acid symbols are explained in Fig. 3.All diagrams are drawn in such a way as to show more degraded samples on the left and less degraded samples on the right side (descending axis of δ 13 C values).In the text, the paleocryosol parts are reported with altitude measurements from the lowest to the highest sample of each paleocryosol.The grey shaded areas are for visualization, not for exact height estimations of the paleocryosols.

Figure 5 .
Figure 5. Scatterplots of selected degradation markers.In all three plots, the x axis shows the C / N ratio.Yedoma deposits are shown as black triangles, thermokarst deposits as grey diamonds.Regression equations, the r 2 , the sample number (n), and the p value are given to the right of the graph.

Figure 6 .
Figure 6.Ordination plots of the principal component analyses (PCA).In diagram (a) the sedimentological parameters are plotted.In (b) a PCA of biomarker proxies is shown.Supplementary variables (in blue: TOC wt % , C / N, δ 13 C, grain size, BD, ice content) were added without including them in the PCA calculation.In diagram (c) the PCA of the major odd n-alkanes is visualized using the same supplementary variables as in (b).

Figure 7 .
Figure 7. Conceptual scheme of the organic-matter degradation state, estimated using the different applied proxies with boxplots.The merged profiles of Yedoma deposit boxplots (white boxes) are shown below the thermokarst deposits (grey boxes).The whiskers illustrate the data range, and the box ends indicate the 25th and the 75th quartile (interquartile range).The vertical lines inside each box show the median (= 50th quartile), including the 95 % confidence intervals, represented by notches.Grey numbers show the sample size.All diagrams are drawn in such a way as to show more degraded samples on the left and less degraded samples on the right side.Thus, the axis of δ 13 C shows descending values.