Switch of fungal to bacterial degradation in natural, drained and rewetted oligotrophic peatlands reflected in δ15N and fatty acid composition

For centuries European peatlands have been degrading along with drainage, land use and climate changes. Increasing pressure on peatland ecosystems calls for a more cost-efficient method to indicate the current state of peatlands and the success of restoration efforts. Metabolic pathways in peatland soils are imprinted in stable isotope compositions due to differences in microorganism communities and their metabolic pathways. Therefore, we hypothesize that depth profiles of nitrogen stable isotope values provide a promising opportunity to detect peatland decomposition or restoration. We studied five peatlands, namely Degerö Stormyr (northern Sweden), Lakkasuo (central Finland) and three mires in the Black Forest (southern Germany). At all locations, cores were taken from adjacent drained (or rewetted) and natural sites to identify δ15N trends that could indicate changes due to drainage and restoration. At all drained (and rewetted) sites we found a distinct peak (“turning point”) of the δ15N values in the center of the drained horizon. We did a fatty acids (FAs) analysis to link our results to microbial community composition. As markers, we distinguished between one fungal-derived FA (C18:2ω9c) and four bacterial-derived FAs. For bacteria, we looked for one general bacterial-derived FA (C14:0), two FAs for gram-positive bacteria (i-C15:0; a-C15:0), and one FA for gram-negative bacteria (C16:1ω9c). In accordance with other studies, our results suggest that fungi dominate the microbial metabolism in the upper aerobic peat horizon. This is reflected by depleted δ15N values. Moving downwards, the drained horizon conditions slowly switch to oxygen limitation. Consequently, fungal-derived FAs decrease whereas bacterial-derived FAs rise. The highest diversity of microbial-derived FAs is indicated by the δ15N turning point. Below the δ15N turning point, oxygen is increasingly limited and concentrations of all microbial-derived FAs are decreasing down to the onset of the permanently waterlogged anaerobic horizon. Peatland cores with restoration successes again show, above the formerly drained horizon, no depth trend of the isotopic values. Hence, we conclude that δ15N stable isotope values reflect microbial community composition, which differs between drained and natural peatlands. Published by Copernicus Publications on behalf of the European Geosciences Union. 300 M. Groß-Schmölders et al.: Switch of fungal to bacterial degradation in peatlands


Introduction
In Europe 70 % of the peatlands are degraded (Joosten and Couwenberg, 2001). Leifeld and Menichetti (2018) reported that degraded peatlands account for 5 % of the anthropogenic CO 2 emissions. Despite this dramatic peat decline, we lack reliable and transferable tools that provide time-and costefficient information of the peatland hydrology status, as we define it.
We determine natural, drained and rewetted hydrological statuses of peatlands, derived by the average thickness of the aerobic horizons and the resulting degree of decomposition.
From a hydrological perspective, peatland soils can be divided in three different horizons. Most biological metabolism and nutrient cycling takes place in the acrotelm (uppermost aerobic peat horizon with living vegetation; Asada et al., 2005a;Artz, 2009;Morris et al., 2011). In the water-saturated catotelm (deeper, anaerobic horizon) organic substrates are decomposed at much smaller rates owing to anoxic conditions (Asada et al., 2005a;Artz, 2009;Lin et al., 2014). In the mesotelm, the peat horizon situated between acrotelm and catotelm, water table levels and oxygen content fluctuate, resulting in shifting aerobic and anaerobic conditions and shifting metabolism processes (Asada et al., 2005a;Artz, 2009;Lin et al., 2014). Clymo and Bryant (2008) therefore defined the mesotelm as a "transition horizon". In degraded peatlands the mesotelm is expanded and former preserved organic substrate is decomposed (Zedler and Kercher, 2005). In an expanded mesotelm conditions differ from aerobic, light and warm conditions in the upper mesotelm to semioxic, dark and cold conditions in the lower mesotelm (Artz, 2009;Lin et al., 2014). The conditions in the former mesotelm will be anaerobic, and microbial activity will be inhibited with rewetting (Andersen et al., 2006;Asada et al., 2005b;Thormann et al., 1999).
We determined the hydrological status by a vegetation analysis, the humification index (HI) after von Post (Silc and Stanek, 1977), the measurement of the water table height, and historical data of the installation of drainage channels. Natural and rewetted sites have a high water table near the surface and are mainly formed by Sphagnum mosses with low humification indices. Drained sites are characterized by low water tables, higher grades of humification, less Sphagnum and more of other moss species. However, determination of macroresiduals and their humification degree in more or less degraded peat is time consuming, needs highly specialized expert knowledge and is thus limited to a small number of samples.
Other common methods to measure peatland hydrology currently are gas emission measurements and measurement of growth heights of peatland vegetation. Gas measurements (e.g., CO 2 , N 2O and CH 4 ) provide an indirect measurement of ongoing decomposition processes (Baldocchi et al., 1988). The method is expensive and labor intensive and does not to give information on drainage history and process dynamics beyond the specific measurement time (Bubier et al., 2003). Measuring vegetation growth is connected to the following problems: (i) not only does the sole growth of mosses indicate peat growth but also the balance of growth and degradation (it is important how much vegetation material enters the catotelm and is therefore stored under anaerobic conditions); (ii) peat shrinks and swells with water supply (hence, measuring peat height at different water table heights would lead to different assumptions for peatland growth; Clymo, 1970); and (iii), peat growth, also in a natural hydrology status, is slow and an unambiguous effect on the success of restoration efforts might need decades of measurements (Clymo, 1970;Fenton, 1980).
As such, and in search of practical indicators, we measured bulk density (BD), the carbon / nitrogen ratio (C/N) and bulk stable isotope values. BD acts as an indicator for decomposition because decomposition processes lead to higher compaction of the peat soil and therefore increasing BD values (Novak et al., 2008). The C/N ratio indicates the degree of decomposition (Malmer and Holm, 1984;Kuhry and Vitt, 1996). With increasing decomposition, a preferential loss of C over N takes place and the C/N ratio decreases. Stable isotopes depth patterns of δ 13 C and δ 15 N in peat have been found to be specific for peatland hydrology (drained, rewetted or natural) in previous studies (e.g., Krüger et al., 2016;Alewell et al., 2011), but the studies were unable to find an explanation for these depth patterns. As degradation is mostly connected to drainage, we hypothesized that an increase of microbial activity is responsible for the change in isotope patterns.
Stable C and N isotopes are correlated with vegetation composition and microbial decomposition processes. As decomposition induces an enrichment of heavy isotopes ( 15 N and 13 C), vegetation is mostly more depleted in 15 N and 13 C than microbial and recycled substrate. Alewell et al. (2011) and Krüger et al. (2014) reported distinct changes in δ 13 C values for palsa peat with the onset of the decomposition of hummocks. Various authors observed the same trend with decomposition in peatlands of other climate conditions (Krüger et al., 2016;Novak et al., 1999;Hobbie et al., 2017;Biester et al., 2014). The distinct δ 13 C depth pattern is a consequence of the use of different sources by fungi and bacteria, as investigated by Kohl et al. (2015) for peat profiles. They conclude that an increasing δ 13 C signal is caused by differences in biomass synthesis and carbon sources used by fungi and bacteria, which was also reported by Lichtfouse et al. (1995) and Baumann et al. (2013). We also found distinct changes in δ 15 N with drainage. It is known that plants preferentially incorporate the lighter 14 N (Högberg, 1997), an effect that is strongly enhanced by mycorrhizal uptake of nitrogen into plants (Hobbie and Högberg, 2012). Plant rooting and the existence of mycorrhiza leads to enriched δ 15 N values in the remaining bulk material (Högberg et al., 1996) because plants and mycorrhiza preferentially process lighter 14 N (Adams and Grierson, 2001;Asada et al., 2005a;Hög-berg et al., 1996;Kohzu et al., 2003;Robinson et al., 1998). However, our study sites are open peatlands with a low occurrence of vascular plants and mycorrhiza. Hence, we assume that these mechanisms cannot be the main drivers of our observed δ 15 N depth patterns. Tfaily et al. (2014) reported that changing microbial abundance and metabolic pathways are correlated with δ 15 N values. Conversely, this would mean that δ 15 N values could reflect the hydrology status. Therefore, we assume δ 15 N values allow us to draw conclusions about whether the observed peatlands have a natural, drained or rewetted hydrology status.
Following previous studies, we use specific terms for the points of change in the stable isotope depth pattern. The points where the stable isotope signals undergo a sudden directional shift with depth are called "turning points", according to Alewell et al. (2011). Furthermore, the bottom of the mesotelm and the onset of the underlying catotelm are marked by the δ 13 C turning point.
To test the idea of changing dominant microbial communities as drivers for isotope depth patterns, we did a fatty acid (FA) analysis of four investigated sites, namely two drained and two natural sites in Degerö Stormyr (central Sweden; 70 km from Umea) and Lakkasuo (southern Finland; 14 km north of Orivesi). FAs are valid markers for indicating the abundance of specific microbial communities in the peat because they are specific and persistent compounds of cell membranes of different species (Bajerski, Wagner and Mangelsdorf, 2017;Finotti et al., 1993;Piotrowska-Seget and Mrozik 2003;Reiffarth et al., 2016). Therefore FAs enable us to make qualitative and quantitative statements about the relative abundance of different microbial communities. We will test the existence of four bacterial markers (namely C14:0 as general marker; i-C15:0 and a-C15:0 indicative for gram positive; and C16:1ω9c indicative for gram negative; Vestal and White 1989;Willers et al., 2015;Zelles, 1997) and one fungal marker (C18:2ω9c; Sundh et al., 1997;Elvert et al., 2003;Willers et al., 2015).
We hypothesize that microbial abundance and diversity are the drivers for the distinct observed δ 15 N depth pattern in natural, drained or rewetted peats. We assume the δ 15 N depth pattern can therefore be used as an inexpensive and less timeconsuming tool to obtain reliable information of peatland hydrology.

Site description
We studied five oligotrophic peatlands (Tables 1-2). All investigated sites are classified as fibric Histosols (HSf; IUSS, 2015; organic soils). Fibric Histosols are classified as soils with a cumulative organic layer and an organic matter amount of 35 % or higher in at least half of the uppermost 80-100 cm and with a high amount (two-thirds) of little decomposed plant residuals (IUSS, 2015). In addition all inves- tigated peatland soils are Sphagnum peats because of their mean annual temperatures (between +1.2 and +7 • C) and their annual precipitation between 523 mm and 1600 ppm (Eurola et al., 1984;Vitt, 2006). Degerö Stormyr (200 m above sea level -a.s.l.) is situated in northern Sweden, at the Kulbäcksliden Experimental Forest near Vindeln, between the rivers Umeälven and Vindelälven (Eurola et al., 1984). It is an acidic mire with minerotrophic conditions and consists of interconnected small mire patches divided by ridges of glacial till. Degerö Stormyr is classified as northern eccentric peatland (Eurola et al., 1984). The climate is characterized as cold, with no dry seasons and cold summers (Dfc zone, based on the Köppen-Geiger climate classification; Peel et al., 2007). In Degerö Stormyr, ditches were installed at the beginning of the 20th century, were closed in 2017 and a natural reestablishment of Sphagnum took place afterwards. The water table is at the surface in the natural part (DN; Nilsson et al., 2008) and in around 10-15 cm depths at the drained location (DD).

Soil sampling and bulk analyses
In May 2012 (Breitlohmisse), June 2012 (Rotmeer), July 2012 (Ursee) and September 2013 (Degerö Stormyr and Lakkasuo) three volumetric peat cores were drilled per site with a Russian peat corer (Eijkelkamp, the Netherlands) at a medium stage of small-scale topography. In Degerö Stormyr cores were sampled in the assumed natural center of the mire (DN) and a 1 m distance from a drainage ditch (1 m depth; DD). In Lakkasuo we took cores at the natural sites (ombrotrophic natural -LN o ; minerotrophic natural -LN m ) and the drained locations (ombrotrophic drained -LD o ; minerotrophic drained -LD m ). For Ursee two cores were taken, namely one in the natural center (UN) and one at the drained edge of the mire (UD). In Breitlohmisse and Rotmeer we took cores in a transect from natural (BN and RN) to strong drained (BD 2 and RD 2 ) sites. Each core has a composite length of 1 m. Here, we focus on the uppermost 60 cm because this part included the drained horizon, and no major changes in isotopic composition were observed at the natural sites below the mesotelm. In all investigated peatlands, the catotelm starts in the natural sites below a 10 cm depth and varied in drained sites but was always visible below a 40 cm depth. Directly after drilling, HIs were determined for each horizon with the von Post scale. The von Post scale has a range from 1 to 10. HI 1 indicates natural condition with undecomposed, completely visible vegetation residuals. HI 10 represents a strongly decomposed horizon without visible vegetation residuals (Silc and Stanek, 1977).
The cores were encased in plastic shells and covered with plastic wrap, stored in coolers, and transported to the laboratory. The cores were sliced in 2 cm sections and every second layer was analyzed, giving a 4 cm depth resolution. Samples were ovendried at 40 • C for 72 h and homogenized with a vibrating ball mill (MM400; Retsch GmbH, Haan, Germany). Stable C and N isotopic compositions were measured with an elemental analyzer combined with an isotope ratio mass spectrometer (EA-IRMS; Inegra2, Sercon Limited, Crewe, UK). Carbon isotopic composition ( 13 C/ 12 C) was expressed relative to the Vienna Pee Dee Belemnite (VPDB) standard and reported in delta notation (‰); stable nitrogen isotopes were expressed relative to the atmospheric nitrogen standard and reported in delta notation (‰). C/N was determined with the mass relationship of the measured bulk content of C and N. Bulk density was measured with volumetric samples, which were weighted before and after drying.
In Degerö Stormyr tree rings of seven individual trees were analyzed (Pinus sylvestris) to obtain information of growth conditions and, therefore, to enhance our knowledge of drainage history.

Fatty acid analysis
Four cores (per site; one drained and one natural core) were selected to do a fatty acid analysis, namely two sites in Lakkasuo (LD o 1 and LN o 3) and two sites for Degerö Stormyr (DD 3 and DN 1). We took subsamples from all cores in the acrotelm (respectively, at the end of the mesotelm in DD) and in the catotelm. At the drained sites of DD 3 and LD o 1, we also took samples in the middle and at the end of the mesotelm. We processed 0.2-1.1 g of sample for the lipid extraction with a mixture of CH 2 Cl 2 : MeOH (9 : 1 ν/ν) in an accelerated solvent extractor (Dionex ASE 350; Thermo Fisher Scientific, Waltham, Massachusetts, USA). A total of 50 µL of an internal standard (0.4 mg mL −1 ; nonadecanoic acid) was added before processing each sample.
The total lipid extracts (TLEs) were saponified by adding 2 mL of KOH dissolved in MeOH (12 %) and putting it in the oven for 3 h at 80 • C.
Following the method of Elvert et al. (2003), TLEs were pooled afterwards with 1 mL KCl (0.1 mol), and the neutral fraction was extracted by agitating three times with hexane. Neutral fraction in the supernatant was separated, dried under a stream of N 2 and stored in the fridge for later analysis. We acidified the rest of the TLEs with fuming hydrochloric acid to a pH of 1. The acid fraction was extracted by agitating again three times with hexane. The acid fraction in the supernatant was separated and hexane was reduced to near dryness under a stream of N 2 . Then the acid fraction was methylated by adding 1 mL boron trifluoride (BF 3 ) in MeOH (12 %-14 %) and putting it in the oven for 1 h at 60 • C. Afterwards the resulting fatty acid methyl esters (FAMEs) fraction was pooled with KCl (0.1 mol) extracted by agitating again three times with hexane and transferred in 2 mL vials. The FAMEs were quantified with a gas chromatograph (Trace Ultra GC) equipped with a flame ionization detector (FID; Thermo Fisher Scientific, Waltham, Massachusetts, USA). The carrier gas (helium) had a constant flow of 1.2 mL per minute and the GC-FID was set to splitless mode. Detector temperature was 320 • C and the samples (dissolved in hexane) were injected by 300 • C. The starting temperature of the oven was 50 • C. The temperature was increased by 10 • C per minute to 140 • C. The temperature was held for 1 min before it was increased up to 300 • C. This temperature was held for 63 min.

Data evaluation and statistical analysis
As we were interested in comparing the depth trends of all single profiles with each other, we first normalized the depths of the cores. This was done using the depth of the δ 15 N turning point (see Sect. 3.1) in each drained profile as the anchor point serving as normalized depth (normD). The normalized depth of this anchor point was set to a 20 cm depth (normD = 20 cm; Fig. 1) in each single core. In the corresponding natural cores, we transferred the values from the same depth related to the drained core into the same norm depth. For example the values of the natural site (DN) in a depth of 13 cm (depth of the turning point of δ 15 N in the corresponding DD core) were set to 20 cm normD.
In a second step, because we were mainly interested in trends and not the absolute values, we normalized the isotopic values themselves because the range of δ 15 N varied considerably between the sites, whereas the trends show consistent patterns (Fig. 1). Therefore, to be able to do a meaningful comparison we set the value of δ 15 N at the turning point to zero in each profile as follows: at turning point.
Using the same procedure, all other parameters (δ 13 C, C/N and BD) were normalized using the same anchor point (e.g., δ 15 N turning point) as follows:  We also determined the slopes of each single core to acquire information on the strength of the differences of the isotopic values with depth. First, the whole peat profile of each drained core was analyzed as one trend (called "overall profile"). Second, profiles were separated into the following different horizons: (i) rewetted horizon (if present), (ii) upper mesotelm, (iii) lower mesotelm and (iv) catotelm. If the values were clearly changing with depth, then the slopes were closer to zero. In horizons with stabilized values, the slopes were distinctly higher or lower than zero.
In the following we present only the normalized data. Raw data without normalization are available in the Supplement.

Tree ring width and microscope analysis of peat
The investigation of the tree ring width of seven surrounding trees (Pinus sylvestris) in Degerö Stormyr was done with a hand-operated wood drill (Djos Increment Borers, Haglöf Sweden, Västernorrland, Sweden; 5 mm diameter). Samples were fixed on wooden carriers. The tracheids (elongated cells of the xylem of vascular plants) were cut with a sharp carbon blade and analyzed with an impinging light binocular (60×-160× amplification).
Peat samples of four study sites were analyzed using an impinging light binocular (60×-160× amplification) to obtain an overview of the vegetation assemblages and to differentiate horizons. For detailed information (and the distinction of the Sphagnum species) the samples were elutriated with water, pigmented with methyl blue and analyzed under a transmitted light microscope (100×-640× amplification).

Depth profile of vegetation assemblage and water table defining the hydrological statuses
Following our indicators (HI and vegetation assemblages), we defined the following three types of hydrological statuses: (a) natural, (b) drained up to the surface and (c) profiles with a rewetted horizon above the drained horizon (Fig. 1).
All drained sites had higher HIs even if no direct modifications in the vegetation assemblage could be documented. For type (b), there was little or no Sphagnum visible at the surface and the water table was found at lower depths (Sect. 2.1). Macroresiduals were more strongly affected by decomposition and HIs were high up to the surface. Especially the ombrotrophic-drained site (LDM) was influenced by drainage. Here, mosses of drier environments replaced Sphagnum species or mosses were completely absent (Tables 3 and S4).
For type (c), vegetation assemblages were mainly composed of Sphagnum spp. and the water table was near the surface. HIs were low in the rewetted horizon and macroresiduals were preserved well (Sect. 2.1; Tables 3 and S4). With the onset of the upper mesotelm, HIs and decomposition of macroresiduals was high. In the lower mesotelm, the HIs were decreasing and more macroresiduals were visible. In the catotelm, the quality of macroresiduals was higher than in the mesotelm and the HIs were even lower (Tables 3 and  S4).

Tree ring widths are verifying the rewetted hydrological status of Degerö Stormyr
Tree ring width is a marker for the wellbeing and/or growth rate of trees. Young trees have a small circumference, coupled with high growth rates, which leads to thicker tree rings. Tree rings get smaller with the increasing age of the tree. If there are no environmental stressors like heat, increasing wetness or drought, tree rings are bigger and the cell lumen is higher compared to trees at sites with environmental stress. With increasing environmental stress, tree ring width decreases (Stoffel et al., 2010). Before 1992, tree rings at the drained site (DD) site showed only a slightly decreasing trend, which could be due to aging (average of 1.3 mm width in the 1930s to an average width of 0.9 mm in the late 1980s). The draining ditches in Degerö Stormyr were established at the beginning of the 20th century, which supports these results, with dryer and therefore better growth conditions for trees. From 1992 onwards tree ring widths decreased, reaching 0.2 mm in 1998 and thereafter. These results suggest a restoration to a wetter, i.e., more natural hydrological status. Rewetted hydrological conditions are not favorable for tree growth and thus lead to smaller tree ring width.

Biogeochemical parameters and hydrological status
Biogeochemical composition of peatlands strongly reflects the related hydrological status (Moore and Basiliko, 2006). As is typical for oligotrophic peatlands, our investigated natural sites have an average C/N ratio of 57 (Table S7). This is in line with results from Malmer and Holm (1984) and Kuhry and Vitt (1996) who found the C/N ratio in the acrotelm of oligotrophic peatlands to be higher than 35 (mostly between 50 and 90). The values in the mesotelm were lower compared to both acrotelm and catotelm, most likely due to higher decomposition rates and the release of CO 2 (Table S7). As is typical for peatlands, BD in our peatlands is low due to the high amount of plant residuals in the soil and low values of mineralization (Novak et al., 2008), with 0.02 kg m −3 at the surface and increased with increasing decomposition and compaction of plant material downwards to 0.04 kg m −3 in the mesotelm (Table S8). BD was also increasing in the catotelm (average of 0.05 kg m −3 ; Table S8), following the increased gravimetric pressure.
In contrast, the biogeochemical parameters of drained sites have a very different pattern. The lower C/N ratio in the acrotelm (average of 41; Table S7) and the mesotelm (average of 35; Table S7) indicates higher mineralization rates with gaseous release of carbon and nitrogen (Krüger et al., 2017). In the catotelm with natural, anaerobic conditions, the C/N ratio was in the same range as in the natural sites (average of 49; Table S7). BD of the acrotelm and mesotelm (average of 0.07 kg m −3 ; Table S8) also increased as a consequence of the enhanced decomposition processes.
These results are in line with the hydrology statuses indicated by the vegetation analysis (Sect. 3.1).

Stable nitrogen isotope depth trends as indicators for the hydrology status
While mineral soils have been shown to have continuous increasing values of δ 15 N (Nadelhoffer et al., 1996;Högberg et al., 1997), we found increasing δ 15 N and δ 13 C values with depth down to particular isotopic-specific turning points in drained peatland soils (Fig. 1). The trends of the single eight out of nine studied drained peatlands and the average trend confirm the existence of a δ 15 N turning point. We determined a significant difference with p < 0.05 between δ 15 N in the center of the mesotelm compared to the δ 15 N values in undrained horizons (Table S1), with a nonsignificant difference for one drained site, namely Breitlohmisse (BD 1 ). The latter was most likely related to generally higher δ 15 N values of the natural site in Breitlohmisse (BN) compared to a smaller increase in δ 15 N at the related drained site (BD 1 ). The depth of δ 15 N turning point (center of the mesotelm) differs from δ 13 C turning point (end of the mesotelm) for all investigated sites (Fig. 2).
Changed slope values of the separated horizons indicate significant trend changes (Table S3). In anaerobic conditions (natural and catotelm) with stabilized isotopic values with depth, slopes were distinctly different to 0 (cm ‰ −1 ). δ 15 N values seem to change rapidly within the mesotelm and slope values were closer to zero. Most interesting was a switch to negative trend values at the δ 15 N turning point in all investigated drained sites, which marks the beginning of the lower mesotelm (Table S3). In a natural hydrological status -type (a) -all investigated parameters had a low variability and indicated a natural, wet mire hydrology status (Fig. 1). There were two exceptions, namely Breitlohmisse natural (BN; 40-60 cm normD) and Rotmeer natural (RN; 30-50 cm normD), with trend instabilities of δ 15 N. This might indicate some minor drainage or disturbance in the wetland sites we classified as "natural" (Fig. 1).
In contrast, the values of the drained sites showed significant trends. We found two different trend types in the drained sites, namely type (b) and (c; Fig. 1). For type (b) we distinguished six sites, namely Lakkasuo ombrotrophic drained (LD o ); Breitlohmisse natural dry (BNd d ); Breitlohmisse drained 1 (BD1); Breitlohmisse 4 (BD 2 ); Rotmeer drained 1 (RD 1 ); and Rotmeer drained 2 (RD 2 ) with clear signs of decomposition up to the surface. Type (c) was visible in three sites, namely drained site Degerö Stormyr (DD), minerotrophic drained site Lakkasuo (LD m ), and Ursee 1 (UD). At type (c) sites the isotopic values, C/N and BD, were stabilized again above the mesotelm. Therefore, they are assumed to be in a "new" natural status (Figs. 1-2).
Below 8 cm (normD; average profile) all drained profiles showed the typical signs of the upper mesotelm with increasing values of δ 15 N, δ 13 C and BD down to the δ 15 N turning point and decreasing C/N. Below the δ 15 N turning point, in the lower mesotelm, δ 15 N values were decreasing. In this horizon, δ 13 C values, C/N and BD were increasing. The end of the lower mesotelm was mostly linked to a clear shift in δ 13 C trend to either stable values or a slow decreasing trend; hence, we called this point δ 13 C turning point (28 cm normD; average profile; e.g., Krüger et al., 2014). Constant C/N, BD and δ 15 N values below the δ 13 C turning point also served as indicators for reduced compaction and decomposition. Most likely the δ 13 C turning point marked the onset of permanent waterlogged anaerobic conditions (e.g., Krüger et al.,

Changing microbial FAs and nitrogen stable isotope depth pattern
Fungal-derived FAs (80 % of all microbial-derived FAs) were the dominant fraction near the surface. In the catotelm the microbial-derived FA values decreased down to 30 %, compared to the acrotelm and the mesotelm, with a clear dominance of bacterial-derived FAs (98 %), as a consequence of the anaerobic conditions (Fig. 3). The latter is congruent with the results of Thormann et al. (1999); fungi will be outcompeted by bacteria with increasing depth and changing hydrological conditions (darker, less oxygen). In the acrotelm of the natural sites, 70 % less microbial-derived FAs, compared to the acrotelm of the drained sites, confirmed the clear link between microbial abundance and the hydrological status. In contrast, we found similar values of microbial FAs in the catotelm for drained and natural sites. This suggests that drainage did not affect the catotelm. In the drained sites the enhanced microbial-derived FAs, abundance could be caused by the improved conditions for metabolic processes by drainage, i.e., enhanced oxygen abundance and relatively high nutrient availability of the prior conserved plant material (Peltoniemi et al., 2009). In the acrotelm and the upper mesotelm fungal-derived FAs were dominating (77 %). At the δ 15 N turning point, lower values of fungal markers (23 %) and increased bacterial- derived FAs (67 %) could be found. In the lower mesotelm, the abundance of microbial-derived FAs was generally decreased and 69 % of the detected FAs were bacterial derived (Fig. 3).

Microbial metabolism mirrored by stable isotope patterns
Our findings suggest that nitrogen stable isotope values are linked to microbial abundance and diversity. We found a clear correlation for stable isotope depth patterns and microbialderived FAs in all sites (r 2 = 0.4; Fig. 3), with high values of nitrogen stable isotopes being linked to high amounts of microbial-derived FAs.
Generally, plants are depleted in 15 N compared to atmospheric nitrogen (which is, by definition, 0 ‰ because air is used as the nitrogen isotopic standard) due to the general preference of plants for the lighter isotope 14 N. As such, the average signal of the relatively undecomposed peat (e.g., the acrotelm of the natural/rewetted sites; the catotelm) is −10 ‰ to −4 ‰. These plant signals are imprinted in the acrotelm (average of −6.09 ‰; Table S5). Furthermore, δ 15 N values of plants (here mostly Sphagnum mosses) are lower than the values of microbes and bulk material (Aldous, 2002;Lichtfouse et al., 1995). Microbes prefer to mineralize the lighter 14 N and plants incorporate (and therefore extract) the microbial mineralized lighter nitrogen (Dijkstra et al., 2006;Novák et al., 1999). Contrary to plants, microbial biomass is enriched in 15 N, probably as the result of processing and releasing the lighter 14 N during mineralization and hence sequestering the remaining heavier 15 N (Roswell, 1976). In addition, as caused by the preferential mineralization of lighter nitrogen, the heavier 15 N might be also enriched in the remaining humic substances (Novák et al., 1999). The effect of the latter on δ 15 N bulk values is probably also enhanced due to the loss of 15 N-depleted material during leaching (Damman, 1988;Niemen, 1998), denitrification and the release of gaseous nitrogen (Kohzu et al., 2003;Niemen, 1998). Our values confirm these reported patterns, with highest δ 15 N values in the mesotelm (average of −3.63 ‰; Table S5) and the correspondence of high microbial activity (reflected by the highest values of microbial-derived FAs) to the δ 15 N turning point . In acid bogs under aerobic conditions, fungi will dominate the general metabolism in upper peat soils (Thormann et al., 2003). This is pictured by the highest amount of fungal-derived FAs in the acrotelm and the upper mesotelm (Fig. 4). Fungi are the preferred decomposers of primary plant material (Wallander et al., 2009;Thormann et al., 2004); hence, the depleted plant isotopic signal is relatively preserved in the upper most aerobic horizons. Furthermore, fungi have a relatively low nitrogen demand compared to bacteria (Myers et al., 2012). With increasing depth and increasing oxygen limitation, fun- gal metabolism decreases (Thormann, 2011). In parallel, the amount of bacterial-derived FAs increases (Fig. 4) as Lin et al. (2014), Hu et al. (2011) andBauersachs et al. (2009) also reported. They found evidence for bacterial-dominated decomposition in hypoxic conditions. This is in line with the findings of Kohl et al. (2015) and Schmidt and Bölter (2002), who also reported a switch from fungal to bacterial dominance in the mesotelm. Andersen et al. (2013), Wallander et al. (2009), Winsborough and Basiliko (2010 and Myers et al. (2012) also stated that fungal biomass is decreasing in peatland soils with depth. In addition, bacterial metabolism is generally faster than fungal metabolism and needs higher amounts of nitrogen (Brunner et al., 2013). We assume that bacteria and fungi compete most over decomposable substrates (not only nitrogen) at the δ 15 N turning point, resulting in the highest turnover rates with an enrichment of δ 15 N in the remaining peat, similar to reports from mineral soils with aerobic decomposition (Alewell et al., 2011;Nadelhofer et al., 1996). As such, we assume that besides the highest microbial activity, the diversity of microbial metabolism also peaks at the δ 15 N turning point (Fig. 5). This would also be related to the highest δ 15 N values because (1) different microbial communities prefer different sources (Dijkstra et al., 2006;Drollinger et al., 2019) and (2) with increasing bacterial abundance, fungi have to also use recalcitrant (isotopi-cally lighter) sources because bacterial metabolism will outcompete fungi for the easily degradable substances (Rousk and Bååth, 2007;Winsborough and Basiliko, 2010). Hence, with increasing microbial diversity, the diversity of the mineralized organic fractions also increases (Thormann, 2006).
To summarize, with an increased diversity of utilized nitrogen sources, more release of lighter 14 N is possible and the δ 15 N values in the remaining substrate should increase (Dijkstra et al., 2008;Figs. 3-4). However, because of the faster and more complete decomposition with increasing microbial activity (Damman, 1988), metabolism of 15 N increases as well and fractionation will be less (Lerch et al., 2011). These contrasting patterns must lead to only small increases in the δ 15 N values of the bulk material, because if all nitrogen is used, fractionation will be lower at the δ 15 N turning point.
In the lower mesotelm, oxygen limitation increases, leading to a general decrease in microbial metabolism and related concentrations of microbial-derived FAs (Fig. 4). The decreasing microbial metabolism leads to simultaneously decreasing δ 15 N values because an increasing amount of intact vegetation (with low δ 15 N values) will be conserved .
Finally, with the establishment of permanently waterlogged anaerobic conditions in the catotelm (also indicated by the δ 13 C turning point), FA concentration decreases sharply to near-zero values. Here, decomposition processes are largely inhibited, which leads to stable δ 15 N and δ 13 C values, close to the original vegetation signals (Alewell et al., 2011;Krüger et al., 2015;Figs. 1 and 4).

Conclusions
Our results confirmed that the nitrogen isotopic depth trends of peatlands are suitable indicators of the natural, drained or rewetted hydrological status. We validated our isotopic hypothesis with microscope analysis of the vegetation remains in the cores and the investigation of tree rings as indicators for changed hydrological status in the past. An analysis of gram-positive and gram-negative bacterial-derived FAs versus fungal-derived FAs underpinned our hypothesis with the expected changes in microbial abundance with depth. The aerobic acrotelm was characterized by a high fungal abundance with low nitrogen demand and turnover. The upper mesotelm was the transition to a mixture of decreasing fungal and increasing bacterial abundance, competing for organic substrates, and resulting in an enrichment of δ 15 N values. In the lower mesotelm microbial decomposition generally decreased but was dominated by bacterial abundance, and finally, microbial metabolism was strongly impeded and δ 15 N values were stabilized in the anaerobic catotelm.
Carbon isotope compositions were also changed with drainage, but they are neither a suitable indicator for a switch in microbial abundance within the drained horizon nor for the trend induced by the rewetting of the peatland. In summary, δ 15 N depth profiles in peat might give more insights into a switch of microbial metabolism because they reflect more precisely different microbial abundance than carbon isotope compositions. Therefore, we conclude that δ 15 N depth profiles could act as a reliable and efficient tool to obtain fast and easy information about the hydrological status, restoration success and drainage history. Data availability. All used data are available in the Supplement.
Author contributions. MGS conducted the sampling, measurements, evaluation and analysis of the data and cowrote the paper with CA. JPK did the sampling and measuring, and AB assisted with the measurements and helped with the analytics. KK added to the discussion. PvS provided the peat microscopy and vegetation analysis. JL and CA had the project idea, supervised and added to the discussion.