Highlights

  • Rewetting increased biomass production and supported the peatlands C sink functions

  • High biomass production compensated for decomposition losses even during a dry year

  • Root biomass was important for organic matter accumulation in graminoid fens

Introduction

Peatlands represent the most space efficient carbon (C) stores among terrestrial ecosystems. Although covering only about 3% of the terrestrial surface, their total C pool exceeds that of the world’s forests and is equivalent to about 60% of the C in the atmosphere (Joosten and Couwenberg 2008; Joosten and others 2016). In the past 200 years, human activities have reduced the area of pristine, peat accumulating peatlands by about 50% in Europe (Joosten and Clarke 2002; Tanneberger and others 2017). Drainage, for example, for agriculture and forestry, has led to aeration of the peat and its subsequent decomposition, effectively turning peatlands into C sources, jeopardizing their climate change mitigation function (Joosten 2009; Joosten and others 2016; Leifeld and others 2019). To counteract the negative effects of peatland drainage (for example, carbon loss, eutrophication, surface subsidence, loss of water retention and loss of habitat), drained sites are rewetted worldwide. However, many site characteristics have fundamentally and irreversibly changed during drainage (for example, bulk density, nutrient content, species pool) and, as a result, rewetting often does not lead to the restoration of the original ecosystem functions (Lamers and others 2015; Bonn and British Ecological Society 2016). Rewetting may result in recovery of the C sequestration function (Tuittila and others 1999; Graf and Rochefort 2015; Bérubé and Rochefort 2018; Günther and others 2015). This process is usually driven by low rates of decomposition, caused by water saturation leading to anoxia, which inhibit the activity of aerobic decomposers (that is, bacteria and fungi) (Denny 1993; Freeman and others 2001; Moore and others 2007). The balance between production and decomposition of organic matter (OM) is thus crucial for evaluating a rewetted peatland’s potential to store C again. In addition, this balance can be considerably influenced by extreme weather and climate events, as we experienced them in the hot and dry year 2018 (Schuldt and others 2020) and that are expected to occur more often in the future (IPCC 2013; Rahmstorf and Coumou 2011; Stott 2016). Increasingly warm and dry conditions introduce oxygen into these systems, stimulating microbial growth and raising the breakdown of OM and the release of C (Fenner and Freeman 2011). Evidence exists that also peatlands that often experience short-term droughts (that is, over summer) continue to preserve C depending on their vegetation (Wang and others 2015). The balance between production and decomposition depends not only on abiotic factors, such as water table depth (Wieder and others 2006) and soil temperature (Moore and others 2007), but also heavily on the vegetation present. For example, litter quality alone explains about 65% of litter decomposition rates (Thormann and Bayley 1997; Zhang and others 2008; Parton and others 2007; Djukic and others 2018), and even within a single plant species, litter decomposability can differ with roots often decomposing slower than shoots (Zhang and others 2008; Freschet and others 2013). Thus, root production is potentially more important for peat formation than aboveground production. In fact, in the extensive river valley peatlands of lowland Europe, the peat is formed primarily by roots and rhizomes and not by mosses (Carex-peat) (Succow and Joosten 2001; Michaelis and others 2020). Nevertheless, roots and their potential contribution to peat accumulation remain understudied.

Fine roots (< 2 mm in diameter, Pregitzer and others 2002) are short-lived roots whose main function is the uptake of nutrients and water from the soil (Guo and others 2008). Fine root production accounts for up to 22% of the worldwide annual net primary production (McCormack and others 2015) and controls, besides nutrient and water uptake, a large flux of C from plants into the soil in form of root exudates and fine root turnover (McCormack and others 2013; Pendall and others 2004). Despite the major role of belowground litter as a driver of ecosystem OM dynamics (Freschet and others 2013), decomposition studies mainly focus on leaf litter (Aerts 1997; Meentemeyer and Berg 1986; Zhang and others 2008) or use standard material (Domisch and others 2006; Laiho and others 2004; Minkkinen and others 1999; Lieffers 1988). However, such studies only describe the effects of environmental factors on decomposition (that is, decomposition potential), without the effect of litter quality of the actual present vegetation (Laiho 2006; Hartmann 1999). The use of local plant material, including roots of different soil depths, along with standard material potentially allows to disentangle the effect of abiotic factors and the litter quality for a given site.

Here, we studied three characteristic peatlands of the North European Plain (that is, lowlands affected by Pleistocene glaciation): an alder forest, a percolation fen and a coastal fen, each in a drained and a rewetted state. All these peatlands are fens (that is, minerotrophic peatlands that receive surface or ground water) and their peats are comprised of roots. The objective of this study was to measure above- and belowground biomass production and decomposition in both drained and rewetted sites in these three different fen types, in order to compare their peat-forming potential after rewetting. Compared to the long-term average values of the last climate reference period, the year 2018 was a very warm and dry year. So far, we lack knowledge of the effects that these drought events will have on plant production and decomposition and thus the C-storage function of rewetted fen peatlands.

We expected that rewetted sites are characterized by wetter soil conditions compared with their drained counterparts even in an exceptionally dry year. Thus, we hypothesized that rewetted sites have (i) higher biomass production due to reduced drought stress, (ii) lower decomposition rates due to longer periods of water saturation in the soil inhibiting aerobic microbial activity, and finally that (iii) rewetting increases the potential for accumulation of OM, even in an exceptionally dry year.

Materials and Methods

Study Area

All study sites were located in the German federal state of Mecklenburg-Western Pomerania in the North European Plain of north-eastern Germany with an average elevation of less than 6 m a.s.l. (Table 1). The study region has a maritime climate with a mean annual temperature of 8.8°C and mean annual precipitation of 565 mm. January is the coldest month with a mean temperature of 0.6°C and a mean precipitation of 40 mm, whereas July is warmest with a mean temperature of 16.7°C and a mean precipitation of 62 mm (meteorological data provided by Germany’s National Meteorological Service from 1981 to 2010). The years from 2015 to 2018 were, compared to the average climate over the last 30 years, dry and warm years (World Meteorological Organization (WMO) 2018), with annual temperatures of 10–11°C and mean annual precipitation ≤ 450 mm (Germany’s National Meteorological Service).

Table 1 General Description of the Study Sites

Experimental Set-Up

We studied three different minerotrophic peatland types (= fens) with peats mainly consisting of roots: alder forest, percolation fen, and coastal fen (Table 1). Each peatland type included a drained and a subsequently rewetted site. Within each study site, a representative study area with a size of 10 m × 35 m was set-up. The study area was fenced to exclude pasture and wild animals and a central boardwalk has been built in longitudinal direction to avoid ground disturbance (Jurasinski and others 2020). Five plots (3 × 3 m) were evenly distributed approximately 10 m from each other inside the study area with three plots located at the north side and two plots located at the south side of the boardwalk. Sites shared the same macroclimate, but differed in their microclimatic conditions (Table S1). Rewetted sites generally showed moister conditions with clearly higher water tables compared to the drained sites (Figure 1). However, during the warm and dry year 2018, water tables dropped down to − 80 cm even in the rewetted sites (Figure 1, Table S1). Still, rewetted sites experienced water saturation close to or exceeding soil surface for considerably longer periods than their drained counterparts did.

Figure 1
figure 1

Groundwater table and weather data obtained from the local weather stations for the fen types alder forest (A), percolation fen (B) and coastal fen (C). Displayed are water table depths (cm) for the drained (red line) and rewetted (blue line) sites of each fen type, as well as the air temperature (orange line), and soil temperature in two depths of 5 cm (black line), 15 cm (grey line) and 50 cm (dashed black lines). Between August 23 and October 27 of 2018 the water table in in the alder forest (A) fell below − 70 cm, exceeding the depth of the groundwater tube at this site.

Our study compared six sites with replicates within each site but not for each fen type. Appropriate, additional replicates at the fen type level are not readily available in this region, in particular for the rewetted states. Furthermore, monitoring and sampling them would have surpassed our resources. Importantly, our study sites were selected to be representative of each respective fen type in the broader region in terms of basic site characteristics (Table S2). Furthermore, geostatistical analyses in fen peatlands in the region showed that spatial auto-correlation is almost absent (Koch and others 2014, 2016). Therefore, spatial replicates in the field can be seen as independent and their variation as being representative for the respective vegetation type. Based on these aspects, we describe our statistical results as comparisons between the six study sites, but see good reasons for them being representative for the respective fen types.

Study Sites

Alder forests are deciduous fen woodlands, dominated by Alnus glutinosa (L.) Gaertn. with a poorly developed shrub layer, but a species-rich herbaceous layer with plants that prefer base-rich soil conditions, often moving groundwater and a rather constant water table (Joosten and others 2017). Alnus glutinosa can be found on peatland sites and be peat forming due to its adaptions to water saturated conditions and its tolerance to stagnant waters (McVean 1956; Succow and Joosten 2001; Barthelmes and others 2010). The two study sites are stands of the same forest, a former natural alder swamp, that was drained for wood pasture in the end of the eighteenth century and subsequently used as production forests; one site was rewetted in 2003/2004 (Jurasinski and others 2020). The tree layer of the drained alder stand is characterized by A. glutinosa with a few individuals of Fraxinus excelsior L., while the understory is dominated by Urtica dioica L. (c. 45% cover) and Rubus idaeus L. (c. 40% cover). The understory of the rewetted alder stand is dominated by Carex riparia Curtis (c. 70% cover) with occasionally occurring Glyceria fluitans (L.) R.Br. (c. 30% cover) and Solanum dulcamara L. (c. 20% cover). The two sites lie within the same basin and are about 3 km apart.

The drained and rewetted lowland percolation fens were, respectively, located in the hydrologically connected lower Recknitz and Trebel valley, about 8 km apart. In their natural state, percolation fens are characterised by high, stable water tables and a substantial water flux that percolates the whole peat body. Percolation fens are bound to landscapes where water supply is large and evenly distributed over the year and are found in river valleys (Joosten and others 2017). The drained percolation fen site is used as intensive grassland (Jurasinski and others 2020); its vegetation is mainly comprised of Holcus lanatus L. (c. 80% cover), Ranunculus repens L. (c. 60% cover), Poa trivialis L. (c. 30% cover) and Deschampsia cespitosa (L.) P.Beauv. (c. 15% cover). The rewetted percolation fen (rewetted in 1997) is protected for nature conservation and not managed. The vegetation is dominated by a vegetation with Carex acutiformis Ehrh. (c. 1 m height, 80% cover), Phalaris arundinacea L. (c. 40% cover) and scattered individuals of Epilobium hirsutum L., Equisetum fluviatile L. emend. Ehrh. and Lythrum salicaria L.

Originally, the coastal fens of the study region were occasionally flooded with brackish sea water during storm surges. The accumulation of groundwater and the infiltration of brackish water contributed to their formation (Kirchner 1971). The vegetation of the drained (and diked) coastal site is dominated by Elymus repens (L.) Gould (c. 80% cover) and (less abundant) by Juncus effusus L. (c. 35% cover), Deschampsia cespitosa (L.) P.Beauv. (c. 20% cover), Dactylis glomerata L. (c. 20% cover) and Holcus lanatus L. (c. 10% cover). The vegetation of the rewetted site (rewetted by dike removal in 1993) is comprised of Agrostis stolonifera L. (c. 90% cover), Festuca rubra agg. L. (c. 20% cover) and occasionally Juncus gerardii Loisel. The drained and rewetted site are about 200 m apart and both sites are used for cattle grazing.

Plant Biomass Production

Aboveground

Aboveground biomass production of herbaceous plants and shrubs was measured by harvesting biomass on five square areas of 0.33 × 0.33 m within each of the five 3 × 3 m plots (in total ~ 5.5 m2 per plot) at each of the respective sites. Harvest took place in the beginning and end of July, coinciding with mowing in the coastal and drained percolation fens and thus mimicking the current use at these sites to promote the same vegetation composition as outside the fence. In addition, we harvested in the beginning of October 2018 to prevent underestimation of the biomass production because of loss of senescent leaves. Annual leaf production of the trees was measured with litter traps (n = 5 per site) that were emptied at least every 1 to 2 weeks from mid-October until the end of November 2018.

Accumulation of tree stem biomass was estimated using point dendrometer data for annual radial increment values from three trees per site. These radial increments (tree ring width) were then converted to basal area increment and ultimately to volume increment using a form factor of 0.5 (that is, ratio of the volume of the stem to the volume of the cylinder) (Pretzsch 2009). Volume was converted to biomass using a wood density of 495 kg m−3 (Meier 2015).

Belowground

We measured belowground biomass production in two soil depths (0–5 cm and 15–20 cm) with in-growth cores. For this purpose, a plastic mesh with a mesh size of 2 mm (#7 Mesh Plastic Canvas Artist Sheet, size 34.6 cm × 57.5 cm from Darice Inc., Strongsville, OH, USA) was shaped into a tube of 5 cm diameter and a length of c. 72 cm and the lower end was closed off with a plastic plug. In-growth cores were filled with autochthonous peat, free of living roots, taken in the field from each of the two respective depths (at original bulk densities). In each of five plots per site, we installed three in-growth cores at an angle of 45°, resulting a total of 90 in-growth cores. The cores were installed between 20 November and 1 December 2017 and retrieved after 1 year between 14 and 16 of November 2018. After retrieval of the in-growth cores, protruding roots were cut at the mesh and roots were subsequently washed out and dried at 70°C for at least 48 h in a drying oven before weighing. We also measured root production in 40–45 cm depth, but excluded the results, because we were not able to compare with decomposition data in this depth, as root material for litter bags was lacking at this depth at some sites.

Decomposition

Belowground Material

Autochthonous root material was collected at each site between June and October 2017 from 0 to 10 cm and 10 to 20 cm depth; roots sampled at different depths were kept separated in processing and incubation. Roots were stored at + 4°C until processing. Following the rationale of Hobbie and others (2010), we used fresh roots for our decomposition analysis. The collected roots were washed clean from peat and ‘air dried’ at 35°C in a drying oven. A subsample from each site was dried at 70°C for 48 h to determine dry weight. Root material corresponding to 0.5 g dry weight was placed into nylon mesh bags (size 5 cm × 6 cm, mesh size 0.14 mm), which were heat-sealed. To keep the litter bags in place, they were attached to a nylon mesh (tulle with a mesh size of 3 mm, 20 g m−2) that was placed between two sheets of a high-density polyethylene woven mesh (TopZeven, The Netherlands, mesh size 10 × 10 mm, 420 g m−2, final frame size 55 cm × 24 cm). The ‘litter frames’ were installed vertically into the soil by hammering metal plates into the ground, then prying them apart to create a slot large and deep enough for the frames to be slid into and be completely buried in the ground, and then removing the metal plates. The created gaps closed within minutes by lateral pressure of the peat. In each of five plots per site, three ‘litter frames’ were installed and retrieved at the same time as the in-growth cores (six sites × five plots × two depths × three replicates total). After retrieval, the root bags were air dried, cleaned and in-grown roots were removed. Finally, root remains were washed with demineralized water and dried at 70°C for 48 h before being weighed. Litter bags that were torn open or damaged, causing material to be lost, were excluded from analysis (5 bags lost, resulting in a total of 175 root litter bags used for analysis).

Aboveground Material

Freshly senesced aboveground litter was collected from all sites in October 2017. As species composition and thus phenological stages differed between sites, we collected both recently fallen leaves as well as leaves from standing biomass. In the alder forest, we differentiated between alder leaves and leaves of the understory vegetation resulting in two litter type bags for the alder forest. Litter was dried and processed similar to root litter, and cut into c. 1 cm long pieces if necessary (for example, sedge litter). Litter corresponding to c. 0.54 g of dry weight was put into nylon bags, which were then packed into frames (size 22 cm × 9 cm, otherwise: see above) and placed on the soil surface in immediate proximity to belowground litter bags (total 120 aboveground litter bags). Aboveground litter frames were installed and retrieved at the same dates as the belowground litter frames.

Tea Bags

We used commercially available tetrahedron-shaped non-woven polypropylene tea bags with sides of 5 cm containing about 2 g of green and rooibos tea (Lipton, Unilever) as standard material and followed the protocol by Keuskamp and others (2013). The two types of tea material differ in quality; green tea (Camellia sinensis; EAN no.: 8 722700 055525) with high cellulose content being easy to decompose and rooibos tea (Aspalanthus linearis; EAN no.: 8 722700 188438) with high lignin content being difficult to decompose (Keuskamp and others 2013). Tea bags were installed in the same ‘litter frames’ as below- and aboveground litter bags, with one green and one rooibos tea bag for each depth (above ground, 0–10 cm, 10–20 cm and 45–50 cm) and replicate (total of 360 green and rooibos tea bags each). After retrieval, the bags were air dried, cleaned of peat and roots, dried at 70°C for 48 h, and weighed.

The percentage mass loss (%ML) was estimated for the root, shoot and leaf material as well as green and rooibos tea:

$$ \% {\text{ML}} = 100 - \left[ {\left( {100 \times X_{1} } \right)/X_{0} } \right] $$
(1)

where X0 is the initial dry weight of the material and X1 the dry weight after 1 year.

Potential Organic Matter Accumulation

To arrive at a tentative estimate of OM accumulation potential, we combined production and decomposition estimates from this one-year study period and projected the potential remaining organic matter after one year. For this purpose, we first derived decomposition decay rates (k-values, y−1) for root, shoot and leaves from the litter bag experiments using Olson’s single-exponential equation (Wieder and Lang 1982; Olson 1963):

$$ X_{t} = X_{0} e^{ - kt} $$
(2)
$$ kt = - { \ln } \left( {X_{t} /X_{0} } \right) $$
(2a)

Then, the annually produced biomass (g m−2) as derived from in-growth cores and aboveground harvests was inserted as X0 to obtain Xt for t = 1 year. For comparison with the aboveground biomass data, we derived a tentative, but conservative value for total root production in the 0–20 cm soil layer. We took the measured production at 0–5 cm and added three times the production measured at 15–20 cm. This conservative approach will have led to an underestimation of the OM accumulation potential of roots. Results need to be interpreted with caution as they represent only the potential of OM accumulation derived from the obtained annual production and decomposition data of this study in this particular year.

We did not calculate the potential of alder wood for OM accumulation, because its decomposition rates are less meaningful over a one-year study period.

Abiotic Parameters

Soil temperature data were collected using HOBO data loggers (Onset Computer Corporation, Bourne, MA, USA) at 15-min intervals at 5 cm and 15 cm depth with six loggers per site. Air temperature was recorded by CR300 data loggers (Campbell Scientific Ltd., Bremen, Germany) and averaged over 30 min. Data were recorded from weather stations set directly at the field sites (one weather station for coastal fens and alder forest, respectively, due to close proximity of the drained and rewetted sites Table 1). Groundwater table relative to the soil surface was recorded at 15-min intervals with a CR300 logger for the two alder stands and by Dipper-PT loggers (SEBA Hydrometrie GmbH & Co. KG, Kaufbeuren, Germany) for the two percolation fen sites and the drained coastal fen and by a Baro-Dipper and a Dipper-APT logger (SEBA Hydrometrie GmbH & Co. KG, Kaufbeuren, Germany) for the rewetted coastal fen site. In the rewetted alder stand, the water table fell below − 70 cm between 23 August and 27 October 2018, which exceeded the depth of the groundwater measurement tube at this site (Figure 1).

Data Analysis

All statistical analyses and visualizations were done in R version 3.5.3 (R Core Team 2019). We tested for differences in each of the following four response variables: biomass production, mass loss, potential of OM accumulation and the ΔML(roots-rooibos tea) against the three explanatory variables hydrological status (drained and rewetted), fen type (alder forest, percolation or coastal fen), depth (0–5 cm or 15–20 cm, only for belowground data), litter type (alder leaves, shoots and roots), and their interactions with linear mixed-effects model ANOVAs [R package lmerTest, version 3.1-0, Kuznetsova and others (2017)] with plot as a random factor (because there are several measurements on the same plot, at different depths). Prior to statistical analysis we averaged the three data points from each plot to a single value. We conducted the above analysis for all aboveground biomass (shoots and leaves), belowground biomass (roots), and tree stem biomass separately. For the comparison of above- and belowground data, root biomass response variables were integrated over depth (0–20 cm, see above). Tea decomposition data were analysed above- and belowground within a single model, thus the response variable depth was ‘aboveground’, ‘0–5 cm’, ‘15–20 cm’ and ‘45–50 cm’ in this case.

A Tukey’s HSD test [R package emmeans, version 1.3.3, Lenth and others (2019)] was used to test for differences between the wet and dry sites and between the different depths. All data used in the linear models were analysed graphically to test the assumptions of normality and homogeneity of variance, that is, residuals were plotted in normal qq-plots and residual versus fitted plots of the model. If necessary, data were log (biomass production, mass loss of local plant litter and tea, potential of OM accumulation after 5 years), square root (potential of OM accumulation after 1 year) or rank transformed (resulting in a non-parametric test; ΔML(roots-rooibos tea)). Visualisations were done using the packages sciplot [version 1.1-1, Morales (2011)] and ggplot2 [version 3.1.1, Wickham (2016)]. Effects were considered significant at p < 0.05.

Results

Plant Biomass Production

The effect of rewetting on annual biomass production (g m−2 y−1) clearly differed between the three rewetted and the three drained sites (interaction fen type × hydrological status: aboveground: F = 31.7, p < 0.001; belowground: F = 11.1, p < 0.001, Figure 2A). Except for the shoot production in the coastal fens, production was either similar or higher in the rewetted sites compared to the drained sites. Root production contributed between 23 and 66% of total biomass production (excluding tree stems) across the different sites. As we did not measure root production below 20 cm depth, it is very likely that the total belowground biomass production is even higher. However, root production values from 40 to 45 cm depth showed 10–100 times lower root production at this depth (Table S3).

Figure 2
figure 2

A Annual biomass production (g m−2 y−1, excluding tree stem biomass) measured by shoot harvest, leaf litter traps, and root in-growth cores and B annual mass loss (%) measured by litter bags for the fen types alder forest, percolation fen and coastal fen. Displayed are aboveground biomass (that is, shoots = circles, alder leaves = triangles) and belowground biomass (that is, fine roots ≤ 2 mm = squares) of the two depths 0–5 cm and 15–20 cm for the hydrological status drained (white boxes, red symbols) and rewetted (grey boxes, blue symbols). Asterisks and n.s. indicate a significant (***p < 0.001, **p < 0.01, *p < 0.05) or non-significant effect (p > 0.05) of the hydrological status (drained and rewetted), respectively.

In the rewetted percolation fen, both shoot (346 g m−2 y−1) and root production (199 g m−2 y−1, 0–20 cm) were significantly higher than in the drained site (shoot: 80 g m−2 y−1, root: 43 g m−2 y−1) (Tukey HSD, p < 0.001, Table S4). In the coastal fen, shoot production was significantly higher in the drained (257 g m−2 y−1) than in the rewetted site (157 g m−2 y−1) (Tukey HSD, p = 0.045), whereas root production tended to be higher in the rewetted site at 0–5 cm depth (Tukey HSD, p = 0.058), but not at 15–20 cm depth (Figure 2A). Root production was significantly higher at the rewetted than at the drained coastal site when root biomass production was integrated over depth (Tukey HSD, p < 0.001). Root production of the rewetted coastal fen was the highest of all sites and made up 66% of the total biomass production at this site (roots: 300 g m−2 y−1, shoots: 157 g m−2 y−1, Table S4). In the alder forest, tree stem biomass production of Alnus glutinosa for the year 2018 was estimated to be more than 4 times higher in the drained site (908.29 g m−2) than in the rewetted site (215.47 g m−2) (F = 25.7, p < 0.001). Annual leaf production of Alnus measured with litter traps (drained = 223 g m−2 y−1, rewetted = 248 g m−2 y−1) was almost three times higher than understory production (drained = 84 g m−2 y−1, rewetted = 86 g m−2 y−1) (F = 34.5, p < 0.001). However, no significant differences in annual production of leaves, understory shoots or roots were found between the drained and rewetted alder stand (Figure 2A, Table S4).

Decomposition of Local Plant Material

Annual mass loss of plant litter (%ML) differed little between drained and rewetted sites. We measured higher annual %ML and, therefore, faster decomposition for shoots in the drained alder stand and percolation fen than in the respective rewetted sites (Tukey HSD, alder forest: p = 0.013; percolation fen: p = 0.014). However, there was no significant difference in %ML of the fine roots between drained and rewetted sites for any of the fen types, despite the different litter material incubated at each site and depth (Figure 2B). In the alder forest, annual %ML of shoots of the understory was higher than for Alnus leaves (Tukey HSD, p < 0.001, Figure 2B, Table S4). As expected, %ML was site-specific and differed significantly between the three fen types both above- and belowground (above: F = 21.0, p < 0.001; below: F = 47.2, p < 0.001, Figure 2B).

Roots decomposed slower than shoots: In four of the six sites litter decomposition was significantly slower for fine roots than for aboveground biomass; only the rewetted alder stand and rewetted percolation fen showed no significant difference between fine root and shoot decomposition. In the two coastal fens, %ML of aboveground litter (drained ML = 64%, rewetted ML = 56%) was about two times higher than for the fine roots (drained: ML = 29%, rewetted: ML = 31%) (Tukey HSD, p < 0.001); in the drained alder stand and percolation fen, %ML of aboveground litter was 1.5 times higher than for roots (Tukey HSD, alder forest: p = 0.002; percolation fen: p < 0.001, Table S4). There was no significant difference in root decomposition between the depths 0–5 and 15–20 cm (F = 0.18, p = 0.673; Table 2), except for the rewetted alder stand with faster decomposition in the 0–5 cm depth (Tukey HSD, p = 0.001). In the rewetted coastal fen, where production was highest (see above), %ML was lowest compared with the other sites. In the drained alder stand, the %ML of the understorey shoot litter (ML = 60%) was significantly higher than for alder leaves (ML = 42%) (Tukey HSD, p < 0.001), while in the rewetted site there was no difference between the two litter types.

Table 2 Results of the Interactions of the Linear Mixed-Effects Model ANOVAs for Biomass Production (g m2 y−1), Decomposition of Local Plant Material (Mass Loss in %), Potential Organic Matter (OM) Accumulation, and ∆ML(roots-rooibos tea) Fen Type (Alder Forest, Percolation Fen, Coastal Fen), Hydrological Status (Drained, Rewetted) and Depth (0–5 cm, 15–20 cm) for the Above (=Shoots and Alder Leaves) and the Below Litter (= Roots)

Potential Remaining Organic Matter

Modelled remaining organic matter (OM, g m−2) differed between drained and rewetted sites and among the three types of fens (Figure 3). Significant interactions between fen type and hydrological status suggest that the fen type determines in which way rewetting affects potential OM accumulation (aboveground F = 41.0, p < 0.001; belowground F = 9.3, p = 0.001). The rewetted sites showed an overall higher potential for remaining OM than the drained sites (aboveground: F = 63.7, p < 0.001; belowground: F = 13.7, p = 0.001).

Figure 3
figure 3

Potentially remaining organic matter (g m−2) after one year for the fen types alder forest, percolation fen and coastal fen and their hydrological status drained (white bars, red symbols) and rewetted (grey bars, blue symbols). An exponential decay model with the measured yearly fine root and shoot production as initial biomass and the estimated decomposition rate k based on % mass loss of the local plant material was used to project the remaining organic matter (produced over 1 year) after 1 year. Asterisks and n.s. indicate a significant (***p < 0.001, **p < 0.01, *p < 0.05) or non-significant effect (p > 0.05) with the hydrological status (drained and rewetted), respectively. Error bars show the 95% confidence intervals.

The potential for accumulating OM did not differ significantly between above- and belowground biomass (F = 1.8, p = 0.181), but this difference depended on the litter type. In both alder stands, the estimated OM remaining for alder leaves was highest compared with shoots of the understory (Tukey HSD, p < 0.001) or with roots (Tukey HSD, p < 0.001). For roots, only small or even non-significant amounts of OM were projected to remain (95% confidence intervals overlap with zero). In both percolation fens, there was no significant difference between the remains of shoot and roots. However, higher remains of OM were projected for roots than for shoots in the drained (Tukey HSD, p = 0.026) and the rewetted coastal fen (Tukey HSD, p < 0.001). The rewetted percolation fen seemed to have the overall highest potential for accumulating above- and belowground OM, while the drained percolation seemed to have the overall lowest potential.

Decomposition of Standard Material

Green tea decomposed faster than rooibos tea in all cases (Figure 4, Table S5). Against expectation, %ML was not always higher in the drained than in the rewetted sites, nor was there a clear trend with depth. Differences between the drained and rewetted sites were most pronounced in the belowground buried tea bags and less so aboveground (hydrological status × depth, F = 20.8, p < 0.001, Figure 4).

Figure 4
figure 4

Mass loss in % for the standard materials green (dots, upper panels) and rooibos tea (triangles, lower panels) for the fen types alder forest (A), percolation fen (B) and coastal fen (C) and the hydrological status drained (white boxes) and rewetted (grey boxes). Displayed are the different depths in which the tea bags were buried, including aboveground (= 0), and the soil depths (0–5 cm, 15–20 cm and 45–50 cm). Values with the same letter within a tea type are not significantly different (p ≥ 0.05, Tukey’s HSD test).

Decomposability of Different Litter Materials

The decomposability of root material versus rooibos tea depended strongly on the fen type (F = 76.5, p < 0.001; Figure 5). In the alder stands, fine roots decomposed faster and in the coastal fens slower than the rooibos tea, but we found no significant difference between the two depths (0–5 cm, 15–20 cm) or the hydrological status for these two fen types (Figure 5A, C). In the drained percolation fen roots decayed slower (0–5 cm) or faster (15–20 cm) than rooibos tea (Tukey HSD, p < 0.001; Figure 5B). In the rewetted percolation fen, we did not observe a significant difference with depth, but fine roots in the rewetted site decomposed faster than rooibos tea in the upper layer and slower in the deeper layer compared to the drained site (Tukey HSD, 0–5 cm: p < 0.001; 15–20 cm: p = 0.022; Figure 5B).

Figure 5
figure 5

Comparison of decomposition of fine roots and rooibos tea for the fen types alder forest (A), percolation fen (B) and coastal fen (C) in the two depths 0–5 cm (drained: bright red bars, rewetted: bright blue bars) and 15–20 cm (drained: dark red bars, rewetted: dark blue bars). Difference in mass loss between fine root litter and rooibos tea (ΔML(roots-rooibos tea)) showing whether decomposition was faster (> 0) or slower (< 0) for roots than for tea. Asterisks and n.s. indicate a significant effect (***p < 0.001, **p < 0.01, *p < 0.05) or non-significant effect (p > 0.05) with the hydrological status (drained and rewetted) and the depths (0–5 cm and 15–20 cm, brackets) on decomposition, respectively. Error bars show standard error.

Abiotic Parameters

Sites shared same macroclimate, but differed in their microclimatic conditions. Due to the extreme warm and drought conditions in the year 2018, water table drawdown in the rewetted sites was up to 30 m below surface in the rewetted percolation fen (Figure 1B) and even 80 cm in the rewetted alder stand and coastal fen (Figure 1A, C). However, the rewetted sites had higher water tables than the drained sites, resulting in an almost two times higher number of water saturated days (over all depths, rewetted: c 191 days; drained: c 100 days; Table S1). Mean annual air temperatures across all sites varied between 9.3 and 9.9°C, whereas soil temperatures varied between 8.1 and 10.2°C (Figure 1, Table S1).

Discussion

Plant Biomass Production

Rewetted sites had higher biomass production than their drained counterparts in the percolation fen, a higher aboveground biomass production and similar belowground biomass production in the coastal fen, and similar understory and leaf production in the alder stands (Figure 2A), supporting our hypothesis that rewetting creates conditions beneficial for plant production under dry conditions. Higher production through rewetting may seem controversial, the more so, because we compare with drained systems where vegetation was selected for high yields as part of their use as intensive grasslands (percolation fen) and production forest (alder forest). High water tables are known to inhibit plant growth by limiting the availability of oxygen needed for respiration and lowering the redox potential affecting the uptake of nutrients (Voesenek and others 2006; Lambers and others 2008). However, periodically flooded or even permanently water saturated soils with anoxic conditions induce a plant community composed of species that are adapted to these conditions (Mäkiranta and others 2018) and that can have high yields (Oehmke and Abel 2016). A common adaptation of plants to flooding is the formation of aerenchyma (that is, tissue containing enlarged gas spaces) (Moore and Garratt 2006; Voesenek and others 2006). This adaptation is evident in numerous species present in our studied peatlands (for example, Alnus glutinosa, Carex spp., Juncus sp., Phalaris arundinacea) and enables them to transport oxygen into the otherwise anoxic rhizosphere allowing for root growth deep into waterlogged soils (Evans 2004; Blom and Voesenek 1996; Schröder 1989; Voesenek and others 2006).

Such adaptations of the established plant community are most likely the reason why the rewetted percolation fen had high production, although the water table remained close to the surface for the entire measurement period. These wet conditions during the warm and dry year of 2018 rather seem to have been more favourable for biomass production than the potentially too dry conditions at the other sites. The rewetted coastal site was just as dry as the drained site during the growing season and we observed no clear difference in production here. The water table in the rewetted alder stand was considerably higher than in the drained stand, but during the summer months still (very) deep (Figure 1), not leading to any differences in biomass production between the alder stands.

The positive effect rewetting seems to have on the production of plant biomass in the percolation fen, despite of the observed water table drawdown in 2018, suggests that rewetting creates moister soil conditions beneficial for plant productivity in dry years. Rewetted peatlands may thus better cope with droughts, which are projected to increase with climate change (IPCC 2013). Mäkiranta and others (2018) found no significant change in the biomass production of sedges after moderate warming and water table drawdown in boreal fens, while forbs decreased and shrubs increased. Wet(ter) peatlands, such as our sedge-dominated percolation fen, may therefore tolerate future warmer and drier conditions better than, for example, drained forb or grass dominated fens.

Decomposition of Local Plant Biomass

We found lower shoot decomposition rates in the rewetted alder stand and percolation fen than in their drained counterparts but rates did not differ between the drained and rewetted coastal sites. Higher decomposition rates in the drained sites had been expected, as more oxygen would be available for decomposing microorganisms, facilitating the decay of OM (Denny 1993; Freeman and others 2001; Moore and Garratt 2006; Moore and others 2007). Yet, water tables were mostly below the soil surface during summer months. As decomposition of standard material and alder leaves hardly differed between the drained and rewetted sites, plant shoots in the rewetted alder stand and percolation fen are apparently less easily decomposed than shoots of the drained sites.

With respect to belowground biomass, no difference was found in root decomposition between any of the respective drained and rewetted fen types (Figure 2B), even though the rewetted sites clearly had higher water tables, resulting in an almost two times higher number of water saturated days (over all depths, rewetted: c. 191 days; drained c. 100 days, Table S1). A reason for this could be that though generally higher in rewetted sites, water tables were low during the summer months, particularly in both drained and rewetted alder and coastal sites. Nevertheless, water tables were high also during summer in the rewetted percolation fen (above − 20 cm for c. 70% of the time between April and October), while soil temperatures were lower (Figure 1B). That still no difference was found in decomposition of belowground biomass could be due to drought stress resulting in reduced decomposition in the drained site (Mäkiranta and others 2008).

Different patterns in shoot and root decomposition imply that studies focussing only on decomposition of aboveground litter miss important aspects of organic matter dynamics. For instance, we found no significant difference in decomposition of roots and leaf litter between drained and rewetted sites, but shoot decomposition was significantly higher in the drained alder stand (understory) and the drained percolation fen (Figure 2B, Table S4). Thus, the type of litter determines to some degree how fast it is decomposed (Table 2). Unlike shoots or leaves, roots of adapted plants, grow even into anoxic soil layers making them harder to decompose for aerobic microorganisms upon root death. This makes them better candidates to accumulate as OM and underlines that roots are potentially more important for peat formation than aboveground litter material.

It appears that differences in vegetation composition, and therefore litter quality, in relation to hydrologic status can have strong effects on litter decomposition. Incubation of standard material showed that environmental factors, such as temperature and water saturation, are important drivers of litter decomposition in these systems (see below), but that these are obviously offset by changes in local plant material. Decomposition rates differed between plant materials; shoots of the understory were more easily decomposed than Alnus leaves in the alder forest (Figure 2B, Table S4). Zhang and others (2008) similarly reported higher decomposition rates for grass leaves than for broadleaved litter, which they determined were caused by differences in litter quality (that is, C:N, total nutrient content, lignin content), one of the most important direct regulators of litter decomposition. In contrast to our findings, the global assessment of Freschet and others (2013) found lower decomposition rates for fine roots and stems than for leaves.

Potential Remaining Organic Matter

Based on production and decomposition rates of 2018, we tentatively modelled how much of the produced organic matter would theoretically remain of the annual produced biomass measured after one year. Note that 2018 was a dry and warm year at our sites, implying that this year does not reflect historic conditions but instead might reflect future conditions in the light of ongoing climate change (IPCC 2013). Under these conditions, the rewetted percolation fen and both coastal sites showed the highest potential for OM accumulation. Interestingly, high potentials of OM accumulation were not caused by lower decomposition rates in the respective sites, but rather by higher annual biomass production (Figure 3). The high potential in the rewetted compared to the low potential of the drained percolation fen gives an idea of how effective rewetting can be in restoring peat-forming conditions after drainage [see also Mrotzek and others (2020)]. Indeed, decomposition rates of the sites with the highest accumulation potential are only half as high as in sites with low potentials, such as the alder forest (Table S4). In fact, the alder forest’s potential in accumulating OM appears to be low compared to the other sites and equalled zero belowground after only one year.

Alder wood peats are known for their high degree of decomposition (Barthelmes and others 2006; Prager and others 2006) as a result of improved nutrient availability due to the symbiotic relationship of Alnus glutinosa with nitrogen-fixing bacteria (Frankia alni), oxygen enrichment of deeper peat layers via Alnus roots (Claessens and others 2010), and typically fluctuating water levels (Succow and Joosten 2001). This ‘fertilization’ of the soil promotes turnover and leads to high rates of decomposition (Swanston and Myrold 1997). Highest long-term rates of peat accumulation of 50 to 126 g C m−2 y−1 have been shown to occur in alder forests with water tables slightly below the surface when production is stimulated and decay not as much; those accumulation rates are even higher or comparable to the 48–96 g C m−2 y−1 of sedge peats (Prager and others 2006; Barthelmes and others 2006). However, for our dry and warm study period, higher accumulation of OM in the rewetted site could not be validated as the water table in both drained and rewetted alder stands dropped strongly.

Our projections show that the importance of litter type (that is, shoots vs. roots) for peat accumulation cannot be generalised or only projected from aboveground litter for those fen ecosystems. While almost no OM is projected to remain from roots after one year in the alder stands, roots are of special importance in the two coastal fens where their potential for OM accumulation exceeds that of shoots. This emphasizes the importance of roots for peat formation in these coastal ecosystems, at least during dry years. In fact, the potential for belowground biomass for peat formation may be even higher, because we only measured root production until 20 cm depth and interpolated conservatively between 0–5 and 15–20 cm depth.

Decomposition of Standard Material

We compared decomposition of standard material with autochthonous root material to disentangle if environmental factors or rather differences in root litter quality between soil depths cause depth–specific differences in root decomposition. We focused on roots, because we expected different roots in different depths (for example, chemical differences between deep and shallow fine roots, differing rooting depth between species).

We expected to find lower decomposition rates in rewetted sites and with increasing soil depth, due to higher probability of water saturated conditions and declining temperatures in deeper depths. Indeed, the amount of water saturated days increased and temperature decreased with increasing soil depth (Figure 1, Table S1). Yet, decomposition rates of the standard material only decreased with increasing depth in the alder forest (Figure 4A) and not in the other peatland types. While we found significant differences between drained and rewetted sites for the buried tea bags, these differences did not follow our expectations to always be higher in the drained than in the rewetted sites (Figure 4, Table S5). Yet, root litter decomposition did not differ between our drained and rewetted sites. This discrepancy indicates that the litter quality of the vegetation present may override the effects of environmental factors such as water table and soil temperature, as previously found for Carex spp. and Phragmites australis (Hartmann 1999). Indeed, how fast or slow local root material was decomposed compared to the standard material (that is, rooibos tea) depended on the fen type. Hence, local plant material differed in terms of its decomposability. Furthermore, the decomposability of the root material, but not the standard material, depended on soil depth. Such differences may be driven by a differing rooting depth between species, and thus a change in ‘species composition’ in the root material in different depths, or by changes in root chemical composition with depth (for example, greater nitrogen concentrations in shallow roots), even within the same plant species (Pregitzer and others 1997; Burton and others 2012). These aspects stress the importance of using autochthonous litter material to estimate decomposition rates, even at different soil depths, if the goal is to measure realized and not theoretical decomposition rates of the ecosystems.

Conclusions

Our findings suggest that during warm and dry conditions high rates of biomass production and not exclusively low litter decomposition are an important driver of OM accumulation in rewetted minerotrophic peatlands. Based on these results, rewetted peatlands may cope better with the extreme weather conditions that will occur more frequently in the future than their drained counterparts, emphasizing the case for rewetting those systems with respect to carbon storage even during dry periods when soils are not water saturated. The importance of belowground OM compared to aboveground OM differed between fen types in our projections, highlighting the importance of roots as main peat-forming element in graminoid-dominated fen peatlands compared to alder dominated forest peatlands and their crucial role in C cycling.