Compositional stability of peat in ecosystem-scale warming mesocosms

Abstract

Peatlands historically have acted as a C sink because C-fixation rates exceeded the rate of heterotrophic decomposition. Under future warmer conditions predicted for higher latitudes, however, that balance may shift towards higher rates of heterotrophic respiration leading to the release of previously stored C as CO2 and CH4. The Spruce and Peatlands Response Under Changing Environments (SPRUCE) experiment is designed to test the response of peatlands to climate forcing using a series of warmed enclosures in combination with peat below-ground heating from 0 to +9°C above ambient conditions. This experimental design allowed a test of chemical changes occurring within peatland soils following five years of warming. We analyzed samples in the uppermost 2m of peat using Fourier Transform Infrared Spectroscopy (FT-IR) to quantify the relative abundance of carbohydrate and aromatic compounds in the peat. The peat soils were subjected to deep peat heating (DPH) beginning in June of 2014 followed by whole ecosystem warming (WEW) in August of 2015. We found that the relative amounts of labile and recalcitrant chemical compound groups across the full peat depth interval did not significantly change after five years of exposure to warming. This appears the case even though previous studies have shown that net C losses and loss of bulk peat mass to be instability over that time period. Results suggest that the current store of carbon in peatlands are largely compositionally stable leading to no changes the in the ratio of chemical moieties on the initial four-year timescale of this experiment.

About the Authors:

Mackenzie R. Baysinger

Roles Conceptualization, Data curation, Formal analysis, Investigation, Writing – original draft, Writing – review & editing

Current address: Alfred Wegner Institute, Helmholtz Centre for Polar and Marine Research, Potsdam, Germany

Affiliation: Department of Earth, Ocean and Atmospheric Sciences, Florida State University, Tallahassee, Florida, United States of America

Rachel M. Wilson

Roles Formal analysis, Investigation, Writing – review & editing

Affiliation: Department of Earth, Ocean and Atmospheric Sciences, Florida State University, Tallahassee, Florida, United States of America

Paul J. Hanson

Roles Conceptualization, Formal analysis, Funding acquisition, Investigation, Project administration, Supervision, Writing – review & editing

Affiliation: Environmental Sciences Division, Climate Change Science Institute, Oak Ridge National Laboratory, Oak Ridge, Tennessee, United States of America

Joel E. Kostka

Roles Conceptualization, Formal analysis, Funding acquisition, Investigation, Supervision, Writing – review & editing

Affiliations School of Biological Sciences, Georgia Institute of Technology, Atlanta, Georgia, United States of America, School of Earth and Atmospheric Sciences, Georgia Institute of Technology, Atlanta, Georgia, United States of America

Jeffrey P. Chanton

Roles Conceptualization, Formal analysis, Funding acquisition, Methodology, Project administration, Writing – original draft, Writing – review & editing

* E-mail: [email protected]

Affiliation: Department of Earth, Ocean and Atmospheric Sciences, Florida State University, Tallahassee, Florida, United States of America

ORCID logo https://orcid.org/0000-0002-3303-9708

Abstract

Peatlands historically have acted as a C sink because C-fixation rates exceeded the rate of heterotrophic decomposition. Under future warmer conditions predicted for higher latitudes, however, that balance may shift towards higher rates of heterotrophic respiration leading to the release of previously stored C as CO2 and CH4. The Spruce and Peatlands Response Under Changing Environments (SPRUCE) experiment is designed to test the response of peatlands to climate forcing using a series of warmed enclosures in combination with peat below-ground heating from 0 to +9°C above ambient conditions. This experimental design allowed a test of chemical changes occurring within peatland soils following five years of warming. We analyzed samples in the uppermost 2m of peat using Fourier Transform Infrared Spectroscopy (FT-IR) to quantify the relative abundance of carbohydrate and aromatic compounds in the peat. The peat soils were subjected to deep peat heating (DPH) beginning in June of 2014 followed by whole ecosystem warming (WEW) in August of 2015. We found that the relative amounts of labile and recalcitrant chemical compound groups across the full peat depth interval did not significantly change after five years of exposure to warming. This appears the case even though previous studies have shown that net C losses and loss of bulk peat mass to be instability over that time period. Results suggest that the current store of carbon in peatlands are largely compositionally stable leading to no changes the in the ratio of chemical moieties on the initial four-year timescale of this experiment.

Figures

Fig 4

Fig 5

Fig 1

Fig 2

Fig 3

Fig 4

Fig 5

Fig 1

Fig 2

Fig 3

Citation: Baysinger MR, Wilson RM, Hanson PJ, Kostka JE, Chanton JP (2022) Compositional stability of peat in ecosystem-scale warming mesocosms. PLoS ONE 17(3): e0263994. https://doi.org/10.1371/journal.pone.0263994

Editor: Dafeng Hui, Tennessee State University, UNITED STATES

Received: September 15, 2021; Accepted: February 1, 2022; Published: March 2, 2022

This is an open access article, free of all copyright, and may be freely reproduced, distributed, transmitted, modified, built upon, or otherwise used by anyone for any lawful purpose. The work is made available under the Creative Commons CC0 public domain dedication.

Data Availability: All data are publicly available at https://doi.org/10.25581/spruce.093/1820162.

Funding: JPC, RMW, JEK, & MRB were funded by DE-SC0007144, DE-SC0012088, Office of Biological and Environmental Research, Terrestrial Ecosystem Science Program, under United States DOE contracts DE-SC0007144 and DE- DE-SC0012088. https://ess.science.energy.gov/. PJH, DE-AC05-00OR22725, Oak Ridge National Laboratory is managed by UT-Battelle, LLC, for the U.S. Department of Energy under contract DE-AC05-00OR22725. https://www.ornl.gov/” Note: This manuscript has been co-authored by UT-Battelle, LLC under Contract No. DE-AC05-00OR22725 with the U.S. Department of Energy. The United States Government retains and the publisher, by accepting the article for publication, acknowledges that the United States Government retains a non-exclusive, paid-up, irrevocable, world-wide license to publish or reproduce the published form of this manuscript, or allow others to do so, for United States Government purposes. The Department of Energy will provide public access to these results of federally sponsored research in accordance with the DOE Public Access Plan (http://energy.gov/downloads/doe-public-access-plan).” The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.

Competing interests: The authors have declared that no competing interests exist.

1 Introduction

Peatlands are large terrestrial carbon banks that represent 30% of the world’s soil carbon [1]. They are of particular interest because they contain a quantity of carbon roughly equal to the quantity present in the preindustrial atmosphere [2]. High latitude peatlands have been accumulating C since deglaciation and have been important in the global carbon cycle throughout the Holocene [3, 4]. Over this period, carbon fixation and net accumulation into high latitude peatlands has outpaced decomposition due to cold temperatures, high water tables, and low O2 availability. More than 90% of peatland carbon stocks are found in high northern latitudes [2] where they are expected to experience the effects of climate change to a greater degree than peatlands closer to the equator [5]. The ability of peatlands to remain as carbon sinks in a warming climate is uncertain. With warmer climate, peatland carbon stocks could be released as CO2 and CH4 which would lead to greater atmospheric forcing, warmer temperatures and even greater subsequent soil carbon release [6–8].

Respiration and fermentation of organic material in the subsurface of peatlands is mostly anaerobic and results in the production of the end products CO2 and CH4. The decomposition rate and the ratio of the produced gases is influenced by many factors, but is largely impacted by the terminal electron accepter (TEA) which supports microbial metabolism. High energy yield TEAs increase the amount of CO2 produced, while under TEA-depleted, anaerobic conditions, CO2 itself can act as a TEA to produce CH4. However, even under completely methanogenic conditions decomposition rates and CO2/CH4 ratio can vary due to the oxidation state of the organic matter itself [9–13]. There is increasing evidence that peatlands become increasingly methanogenic under warmer conditions [14, 15]. We hypothesize that this may be due to changes in the peat composition during anaerobic metabolism that increasingly favors methane production. For example, high oxidation state organic C is thought to be highly bioavailable [12, 13] and may be used up first leaving behind an increasing peat C pool of lower oxidation state organic matter that favors CH4 over CO2 production. This effect could exacerbate climate warming since CH4 has a much stronger global warming potential than CO2 [16, 17], thus creating a positive feedback between peatlands and climate.

Because microorganism enzyme action breaks down the components of peat at different rates the relative abundance of carbohydrates and aromatic compounds is a useful metric for determining levels of humification/decomposition [7, 18–21]. Carbohydrates (O-alkyl-C) are formed in plants as cellulose and are more bioavailable than other compound classes such as lignin/aromatics [13, 22]. The Nominal Oxidation State of Carbon (NOSC, [12]) allows an estimate of the energetic potential and bioavailabity of organic matter. Degradation of high NOSC organic matter (i.e., carbohydrates) results in a greater energy yield relative to degradation of compounds with lower NOSC. Under anaerobic conditions, higher NOSC compounds appear to be used preferentially leaving behind a pool of increasingly less bioavailable, more reduced compounds (e.g., lipids and unsaturated hydrocarbons) [22, 23]. This phenomenon can be readily observed with increasing depth in the peat column; more mature peat at depth contains a lower percentage of biodegradable organic matter (carbohydrates) than the recently formed peat closer to the surface [10, 19, 20, 22, 24]. A similar trend has also been found for surface peats along a gradient of latitude; carbohydrate content declines under the warmer conditions found closer to the equator [24] because of higher decomposition rates at warmer temperatures.

Hodgkins et al., [24] postulated that the greater amount of carbohydrate in high latitude peats renders the carbon stored there vulnerable to warming. With warming such soils could undergo losses of the carbohydrate-rich fraction. To test this hypothesis, we quantified the carbohydrate fraction of peat soils under varying degrees of in situ warming, from 0 to 9°C above ambient, over five years. In this study, we sampled peat cores from an artificially heated bog (Spruce and Peatland Responses Under Changing Environments (SPRUCE) located in the Marcell Experimental Forest, MN) and measured their carbohydrate and aromatic content using Fourier Transform Infrared Spectroscopy (FTIR) analyses following the quantification method of Hodgkins et al. [24]. We hypothesized that the carbohydrate component groups would be lost from the soil organic matter at a greater rate with increasing temperature relative to other peat soil components as they are the most labile part of the sample, and would therefore be consumed preferentially by the microbial community in the soil matrix.

2 Methods

2.1 Study sites and sample collection

Ethics Statement: SPRUCE is collaboration between the USDA Forest Service and the Department of Energy through Oak Ridge National Laboratory (ORNL). The interagency agreement between the two agencies (2011-IA-123) allows sample collection for ORNL scientists and other SPRUCE collaborators. To evaluate the impacts of temperature treatment on peat carbohydrate degradation, samples were analyzed from a suite of 2 m-long peat cores collected from an ombrotrophic peatland in Minnesota (USA), where the Spruce and Peatland Responses Under Changing Environments (SPRUCE) project is being conducted (for schematic and site description see [6, 25, 26]. SPRUCE is an ecosystem-scale study that warms the vegetation and peatland soil profile to a depth of 2 m within 12.8-m diameter open-top enclosures. The peat was heated first at– 2 m for a period of 1 year starting in summer 2014 (deep peat heating phase), then heating of the overlying air began in August 2015 starting the whole-ecosystem warming phase of the experiment. Five temperature differentials have been obtained (ambient, +2.2, +4.5, +6.8, and +9°C). There were duplicate enclosures subjected to the same temperature at both ambient and ~ + 500 p.p.m.v. above ambient atmospheric CO2 concentrations [26]. Each of these enclosures were sampled for this project. One core was taken from each enclosure and divided into 11 partitions representing 10 cm depth increments up to 50cm, where they were then partitioned every 25cm to 200cm below surface. Whole-ecosystem above ground warming (WEW) was initiated August 2015, following 14 months of deep-peat heating that began June 2014. Water table was measurements using automated sensors have been taken daily beginning in 2015. Soil cores within each enclosure were collected first in August of 2014 and in August 2019, after 4 years of whole ecosystem warming and 5 years of peat-heating.

2.2 Sample analysis

After transport from the field, peat subsamples, including roots, were freeze-dried to ensure shelf stability. The peat was then ground to a homogenous powder in a Spex Sampleprep 5100 mixer-mill. To evaluate the response of peat organic carbon fractions to the ecosystem-scale heating, Fourier Transform Infrared Spectroscopy (FTIR) analyses were conducted on each depth increment (0-10cm, 10-20cm, 20-30cm, etc) of the peat cores collected from all of the 10 treatment enclosures in both 2014 and 2019, resulting in 195 individual samples. FTIR analysis was chosen because it allows us to semi-quantitatively measure the functional composition of the sample [24, 27] (Fig 1). FTIR spectra were collected with a PerkinElmer Spectrum 100 FTIR spectrometer fitted with a CsI beam splitter and a deuterated triglycine sulfate detector. Transmission-like spectra were obtained with a Universal ATR accessory with a single-reflectance system and made from a zinc selenide (ZnSe)/diamond composite. Samples were placed directly on the ATR crystal, and force was applied so that the sample came into good contact with the crystal. Spectra were acquired in % transmittance mode between 4000–650 cm −1 (wavenumber) at a resolution of 4 cm −1, and four scans were averaged for each spectrum. The spectra were ATR-corrected to account for differences in depth of beam penetration at different wavelengths, and then baseline-corrected, with the instrument software. Spectra were then converted to absorbance mode for subsequent data analysis. In this study, we followed the FTIR data processing method of Hodgkins et al [24] (available here: https://github.com/shodgkins/FTIRbaselines) that allows for more thorough quantification of compound classes. Instead of normalizing peak heights relative to other peaks via humification indices, we isolated each compound by normalizing peaks to the integrated area of the entire spectrum. Normalized peak heights are then reported and used for analysis.

[Figure omitted. See PDF.]

Fig 1. FTIR spectra of peat collected from an ambient enclosure at SPRUCE in September 2014 (depth = 25 cm).

Following Hodgkins et al., [24], we focused on the peaks at 1040 cm-1 (O-alkyl of carbohydrates), 1510 cm-1 (phenolic lignin- like structures), and 1630 cm-1 (aromatic structures because these change in systematic ways in peat soils with depth at a single location and with latitude in surface soils. By normalizing the peaks to the area below the spectral curve, we can quantitatively observe changes in the chemistry of our samples.

https://doi.org/10.1371/journal.pone.0263994.g001

Employing the calibration method from Hodgkins et al. [24], the normalized peak heights can then be used to quantify the presence of molecular groups in a given sample. Knowing the percent composition of carbohydrates, etc. in allows us to see how that composition changes over time of warming exposure. Post-treatment peat samples (2019) were compared to pre-treatment samples (2014) to measure shifts (if any) in molecular group presence.

3 Results

Carbohydrate content in the peat column generally decreased with depth down to -70 cm below which the % carbohydrate remained stable or even increased slightly (Fig 2A). The composition ranged from 23–55% at the surface and showed a negative relationship to depth from surface. The samples taken from -200 cm contained 11–31% carbohydrate content. In contrast, the aromatic content increased with depth, to about -70 cm below which the %aromatics remained relatively stable down to 2 m (Fig 2B). The aromatic content ranged from 23–26% at the surface, and steadily increased in content to 26–36% at the -200cm depth (Fig 2B). We compared the carbohydrate and aromatic content of the 2014 samples to the 2019 samples and found that the chemistry of the peat samples did not significantly change after 5 years of exposure to the ecosystem-scale warming. (carbohydrate content p-value = 0.73; aromatic content p-value = 0.84 by Wilcox-test; Fig 2A and 2B).

[Figure omitted. See PDF.]

Fig 2. Organic matter composition.

Percent carbohydrate and aromatic content in samples compared from each temperature treatment in 2019 compared to the 2014 pre-heating values, circles indicate ambient CO2 and squares indicate elevated CO2 treatments. Dashed lines indicate the locally weighted polynomial regression (LOESS) of the pre-heating (2014) values across all enclosures with the ±95% confidence interval indicated by the blue shaded area around each line (MATLAB 2020a). Results for 2019 are plotted as individual symbols for each temperature treatment as indicated in the legend. decreased with depth, with no discernable difference between pre-heating and five-years into the warming. Panel A. Carbohydrate content decreased from surface to -70cm, where carbohydrate content stabilized or slightly increased. Panel B. Aromatic content in samples increased with depth, with no discernable difference between pre-heating and five-years into the warming.

https://doi.org/10.1371/journal.pone.0263994.g002

To look at temperature-specific changes in the composition of the samples we plot the carbohydrate and aromatic content data by depth from surface for the individual enclosure treatments (Fig 3). Samples were taken in duplicate—two individual domes were subjected to identical heat treatment, as represented by the color in Fig 3. The Spruce experiment is a regression-based design (Fig 4) and includes a broad range of temperature levels to produce data to characterization response curves for application to ecosystem and Earth system models [26].

[Figure omitted. See PDF.]

Fig 3. Carbohydrate and aromatic content versus depth.

Carbohydrate content in peat cores taken from each temperature treatment in duplicate heated enclosures (two cores per temperature treatment) from the beginning of deep peat heating (2014) and after 5 years of warming (2019) (upper panel). Lower panel, aromatic content. Circles are the ambient CO2 and squares are the elevated CO2 treatments.

https://doi.org/10.1371/journal.pone.0263994.g003

[Figure omitted. See PDF.]

Fig 4. Regressions versus temperature.

Temperature regressions of percent carbohydrate (panel A) and percent aromatics (Panel B) in the surface peat (≤ 50 cm) collected in 2019. Data are plotted against the actual averaged measured temperature at each peat depth averaged for the month of August 2019. Symbol colors indicate the temperature treatment differential. Regression analyses were not significant for either % carbohydrate or % aromatics (p > 0.1).

https://doi.org/10.1371/journal.pone.0263994.g004

In the ambient and +2.25°C treatments the 2019 measurements closely matched those from the beginning of the warming treatment in 2014. In the +4.5, 6.75 and +9°C treatments the surface %carbohydrates were frequently, but not consistently, lower in 2019 relative to the starting 2014 values (Fig 4). These differences generally diminished deeper in the peat profile. Sample depth was a larger indicator of peat carbohydrate content relative to the heat treatment group for the 5 years of in situ heating.

Additionally, we plot the aromatic content by depth of sample (Fig 3). We did not see variance in the amount of carbohydrate content or aromatic content in the surface peat with differing temperature treatments (Fig 4). These findings are consistent throughout the entirety of the peat column sampled.

We compared the percent carbohydrate and percent aromatics in the elevated and ambient plots for 2019 (Fig 5). Although the median percent carbohydrate was slightly higher in the elevated CO2 treatment relative to the controls, the differences were not significant (T-test, p = 0.1). Similarly, % aromatics were not significantly different between the ambient and the elevated CO2 plots.

[Figure omitted. See PDF.]

Fig 5. Effect of CO2 treatment.

Comparison of percent carbohydrate (panel A) and percent aromatics (panel B) in the ambient and elevated CO2 treatments. Values for the boxplots were taken only from 2019, after the eCO2 treatment had begun and only from surface peat (<50 cm) where the majority of root inputs occurs. Results of T-tests comparing elevated and ambient values are given in the top of each panel.

https://doi.org/10.1371/journal.pone.0263994.g005

4 Discussion

In this study we sampled peat cores from the SPRUCE experimental treatments in northern Minnesota. We hypothesized that over that time period we would observe a decrease in the carbohydrate content of surface peats as carbohydrates preferentially fueled microbial respiration and as surface peat decomposed at a higher rate at higher temperatures [6]. In addition, we also measured changes in the aromatic content which represents the aromatic rings of quinine and amide compounds in the samples. Aromatic groups are more difficult to decompose compared with the more labile carbohydrates, thus they are expected to increase in the peat during decomposition as more labile components are preferentially removed.

The lack of chemical response of surface peat to temperature increase was surprising as previous work has clearly shown that the system is becoming increasingly methanogenic with warming [14, 28, 29]. Additionally, while early work at SPRUCE failed to find evidence of peat contribution to GHG production [14], more recent studies do show that the peat itself is being degraded as the experiment progresses [6, 29] especially at the higher temperature enclosures, and that an increasing fraction of decomposition supports methane production. We hypothesized that this increasingly methanogenic GHG production was the result of peat decomposition that preferentially removed higher NOSC compounds first leaving behind an increasingly less oxidized carbon pool for decomposition via methanogenesis. The lower NOSC residual would then result in more methane production relative to CO2, as dictated by electron balance. However, results in this paper do not yield evidence that this is the case. There does not appear to be an observable preferential loss of high NOSC (e.g., carbohydrate) compounds with temperature treatment (Fig 4) at least over a five-year period during which time methane production rates have increased relative to CO2 production rates. Therefore, the reason for this shift must be due to kinetic factors affecting rates of degradation processes differentially.

The elevated CO2 treatment stimulates primary production which increases root production and exudation especially in the rooting zone (< 30 cm) [30]. We expected that this increase in production would result in increased carbohydrate content in the shallow peat within the elevated CO2 treatments. We observed that the percent carbohydrate in the elevated CO2 treatment was slightly, though not significantly higher than in the ambient CO2 plots (Fig 5). One explanation for this result is that the large C pool in peat is shifting, but too slowly for us to observe on the timescale of this experiment. Further monitoring is necessary to determine whether this trend, of increasing percent carbohydrate in the surface peat, continues or becomes stronger with time. Increasing inputs of carbohydrates in the elevated CO2 treatments may confound our interpretation of changing carbohydrate content with temperature if new carbohydrates are being input to the system as decomposition is also stimulated to remove those carbohydrates. There were fewer carbohydrates in the +2.25°C and +4.5°C treatments, but that trend was not consistent at the higher temperature treatments (Fig 3). Recent lab-based studies show that decomposition rate at these temperatures does increase systematically across all of these temperatures (Wilson et al., 2016; Hopple et al., 2020), but a microbial preference in peat composition simply cannot be observed on the timescale of this 5-year study.

Nevertheless, the depth profiles from the peat (Figs 2 and 3) provide evidence that, in principle, as decomposition progresses, carbohydrates are preferentially removed and aromatic content increases. However, our results show that this process of microbial organisms preferentially removing carbohydrates takes more time than the scope of this experiment. Such observations are consistent with other peats where these measurements have been conducted (e.g., [20, 24]). On the longer time scales captured by our depth profiles carbohydrates are removed and aromatics increase, but the peat C pool is so large that over the timescale of our experiment we might not have the resolution to capture changes in the peat composition with temperature. Based on radiocarbon profiles at the SPRUCE site, these changes occur on a timescale of 100’s to 1000 years [31]. These rates are consistent with those shown by Hodgkins et al. [24] as well.

The stability of the peatland’s carbohydrate composition below 70 cm is not surprising, as both acrotelm and catolem peat has been shown to be stable in its waterlogged, anoxic state [14, 32]. Water table depths during the SPRUCE experiment from 2015 through 2020 have been driven by a balance between precipitation inputs and evapotranspiration losses (as expected). While those losses are higher under the warmest treatments, an abundance of ambient precipitation prior to the 2021 growing season largely sustained water tables high within the acrotelm layer of the S1-bog peatland. The net impact of this extended “wet” period is that most of the peat profile has supported anoxic conditions. Full water table data are available at [33]. That data set includes data through 2019 but will be appended with future observations for 2020, 2021 and all future years of the SPRUCE experiment. The stability of the S1-bog’s water table height throughout the entirety of the WEW measurements, has contributed to keeping the organic matter at depth in long-term storage.

Peat C decomposition is mediated by microbial enzyme activity, which in turn is strongly influenced by environmental factors such as temperature, substrate quality, and moisture [34]. Changing enzyme activities, such as the decline in β-D-glucosidase activity with depth observed by Hill et al. [35] suggest that peat C becomes increasingly less reactive with age. The enzyme β-D-glucosidase is the primary observed C-acquiring enzyme across a range of wetland types and interestingly, the abundance of this enzyme suggests that despite the large pool of C in peat, peatland microorganisms appear to be C-limited (rather than N or P limited; Hill et al., [36]) consistent with large quantities of the C being sequestered in less bioavailable forms. Nevertheless, the breakdown of compounds at depth does occur, but at a slower rate than in the upper surface layer [32]. For example, polyphenol oxidases, enzymes which break down lignin and other phenolic compounds, increase with depth in the peat profile [35], although it should be noted that the opposite trend was observed by Steinweg et al., [34]. In the peat column, depth has been shown to drive the microbial community [14, 37, 38] and resulting enzyme profiles [34]. Temperatures in the deep peat are more stable relative to the surface resulting in deep microbial enzymes that appear to be less responsive to warming temperatures suggesting that, while surface enzyme activity does appear sensitive to warmer temperatures [39], there is little change in C decomposition with temperature in the catotelm [14, 34]. Our results in Fig 2 are consistent with faster removal of carbohydrate C in the surface peat, followed by relatively stable conditions in the catotelm below 40 cm.

In conclusion, our results indicated that, while net peat C loss does occur [6], the current store of carbon in the SPRUCE peatland is compositionally stable, at least on the timescale of this experiment, as the samples did not show a significant change in absorbances at the key wavelengths from 2014 to 2019. Labile components, primarily the carbohydrate groups at the 1040 cm-1 peak, did not respond to the increased warming over the five years of applied belowground temperature treatment. Aromatic groups (1515 and 1630 cm-1) also did not show a response to the five years of applied treatment covered in this study.

Acknowledgments

The authors would like to thank Jana R. Phillips for the organization of the annual peat coring events that enabled the collection of the samples analyzed here. Samantha Bosman assisted with the preparation of the manuscript.

Citation: Baysinger MR, Wilson RM, Hanson PJ, Kostka JE, Chanton JP (2022) Compositional stability of peat in ecosystem-scale warming mesocosms. PLoS ONE 17(3): e0263994. https://doi.org/10.1371/journal.pone.0263994

References

1. Houghton RA. Balancing the Global Carbon Budget. Annu Rev Earth Planet Sci. 2007; 35(1): 313–347.

2. Yu Z, Loisel J, Brosseau DP, Beilman DW, Hunt SJ. Global peatland dynamics since the Last Glacial Maximum. Geophys Res Lett. 2010; 37(13): L13402.

3. MacDonald GM, Beilman DW, Kremenetski KV, Sheng Y, Smith LC, Velichko AA. Rapid Early Development of Circumarctic Peatlands and Atmospheric CH4 and CO2 Variations. Science. 2006; 314(5797): 285–288. pmid:17038618

4. Yu Z. Holocene carbon flux histories of the world’s peatlands:Global carbon-cycle implications. Holocene. 2011; 21(5): 761–774.

5. Collins M, Knutti R, Arblaster J, Dufresne JL, Fichefet T, Friedlingstein P, et al. Long-term climate change: projections, commitments and irreversibility. In Stocker TF, Qin D, Plattner G-K, Tignor MMB, Allen SK, Boschung J, et al, editors. Climate Change 2013-The Physical Science Basis: Contribution of Working Group I to the Fifth Assessment Report of the Intergovernmental Panel on Climate Change. Cambridge University Press; 2013. pp. 1029–1136.

6. Hanson PJ, Griffiths NA, Iversen CM, Norby RJ, Sebestyen SD, Phillips JR, et al. Rapid net carbon loss from a whole-ecosystem warmed Peatland. AGU Advances. 2020; 1(3): e2020AV000163.

7. Leifeld J, Steffens M, Galego-Sala A. Sensitivity of peatland carbon loss to organic matter quality. Geophys Res Lett. 2012; 39(14):

8. Crowther TW, Todd-Brown KE, Rowe CW, Wieder WR, Carey JC, Machmuller MB, et al. Quantifying global soil carbon losses in response to warming. Nature. 2016; 540(7631): 104–108. pmid:27905442

9. Myhre G, Samset BH, Schulz M, Balkanski Y, Bauer S, Berntsen TK, et al. Radiative forcing of the direct aerosol effect from AeroCom Phase II simulations. Atmos Chem Phys. 2013; 13(4): 1853–1877.

10. Hodgkins SB, Tfaily MM, McCalley CK, Logan TA, Crill PM, Saleska SR, et al. Changes in peat chemistry associated with permafrost thaw increase greenhouse gas production. Proc Natl Acad Sci. 2014; 111(16): 5819–5824. pmid:24711402

11. Wilson RM, Tfaily MM, Rich VI, Keller JK, Bridgham SD, Zalman CM, et al. Hydrogenation of organic matter as a terminal electron sink sustains high CO2: CH4 production ratios during anaerobic decomposition. Org Geochem. 2017; 112: 22–32.

12. Keiluweit M, Nico PS, Kleber M, Fendorf S. Are oxygen limitations under recognized regulators of organic carbon turnover in upland soils? Biogeochemistry. 2016; 127(2–3): 157–171.

13. LaRowe DE, Van Cappellen P. Degradation of natural organic matter: A thermodynamic analysis. Geochim Cosmochim Acta. 2011; 75(8): 2030–2042.

14. Wilson RM, Hopple AM, Tfaily MM, Sebestyen SD, Schadt CW, Pfeifer-Meister L, et al. Stability of peatland carbon to rising temperatures. Nat Commun. 2016; 7(1): 13723. pmid:27958276

15. Hopple AM, Wilson RM, Kolton M, Zalman CA, Chanton JP, Kostka J, et al. Massive peatland carbon banks vulnerable to rising temperatures. Nat Commun. 2020; 11(1):1–7. pmid:31911652

16. Ramaswamy V., Boucher O, Haigh J, Hauglustaine D, Haywood J; Myhre G; et al. Radiative forcing of climate change. In: Houghton JT, et al, editors. Climate Change 2001: The Scientific Basis—Contribution of Working Group I to the Third Assessment Report of the Intergovernmental Panel on Climate Change. Cambridge Univ. Press, New York; 2001. pp. 350–416.

17. Neubauer SC, Megonigal JP. Moving beyond global warming potentials to quantify the climatic role of ecosystems. Ecosystems. 2015; 18(6): 1000–1013.

18. Baldock JA, Oades JM, Nelson PN, Skene TM, Golchin A, Clarke P. Assessing the extent of decomposition of natural organic materials using solid-state 13C NMR spectroscopy. Soil Res. 1997; 35(5): 1061–1084.

19. Grover SP, Baldock JA. Carbon decomposition processes in a peat from the Australian Alps. Eur J Soil Sci. 2010; 61(2): 217–230.

20. Grover SP, Baldock JA. Carbon chemistry and mineralization of peat soils from the Australian Alps. Eur J Soil Sci. 2012; 63(2): 129–140.

21. Treat CC, Wollheim WM, Varner RK, Grandy AS, Talbot J, Frolking S. Temperature and peat type control CO 2 and CH 4 production in Alaskan permafrost peats. Glob Chang Biol. 2014; 20(8): 2674–2686. pmid:24616169

22. Tfaily MM, Cooper WT, Kostka JE, Chanton PR, Schadt CW, Hanson PJ, et al. Organic matter transformation in the peat column at Marcell Experimental Forest: humification and vertical stratification. J Geophys Res Biogeosci. 2014; 119(4): 661–675.

23. Pracht LE, Tfaily MM, Ardissono RJ, Neumann RB. Molecular characterization of organic matter mobilized from Bangladeshi aquifer sediment: tracking carbon compositional change during microbial utilization. Biogeosciences. 2018; 15(6): 1733–1747.

24. Hodgkins SB, Richardson CJ, Dommain R, Wang H, Glaser PH, Verbeke B, et al. Tropical peatland carbon storage linked to global latitudinal trends in peat recalcitrance. Nat Commun. 2018; 9(1): 1–3. pmid:29317637

25. Hanson PJ, Childs KW, Wullschleger SD, Riggs JS, Thomas WK, Todd DE, et al. A method for experimental heating of intact soil profiles for application to climate change experiments. Glob Chang Biol. 2011; 17(2): 1083–1096.

26. Hanson PJ, Riggs JS, Nettles WR, Phillips JR, Krassovski MB, Hook LA, et al. Attaining whole-ecosystem warming using air and deep-soil heating methods with an elevated CO 2 atmosphere. Biogeosciences. 2017; 14(4): 861–883.

27. Biester H, Knorr KH, Schellekens J, Basler A, Hermanns YM. Comparison of different methods to determine the degree of peat decomposition in peat bogs. Biogeosciences. 2014; 11(10): 2691–2707.

28. Wilson RM, Tfaily MM, Kolton M, Johnston ER, Petro C, Zalman CA, et al. Soil metabolome response to whole-ecosystem warming at the Spruce and Peatland Responses under Changing Environments experiment. Proc Natl Acad Sci. 2021; 118(25). pmid:34161254

29. Hopple AM, Pfeifer-Meister L, Zalman CA, Keller JK, Tfaily MM, Wilson RM, et al. Does dissolved organic matter or solid peat fuel anaerobic respiration in peatlands?. Geoderma. 2019; 349: 79–87.

30. Malhotra A, Brice DJ, Childs J, Graham JD, Hobbie EA, Vander Stel H, et al. Peatland warming strongly increases fine-root growth. Proc Natl Acad Sci. 2020; 117(30): 17627–17634. pmid:32661144

31. McFarlane KJ, Hanson PJ, Iversen CM, Phillips JR, Brice DJ. Local spatial heterogeneity of Holocene carbon accumulation throughout the peat profile of an ombrotrophic Northern Minnesota bog. Radiocarbon. 2018; 60(3): 941–962.

32. Clymo RS. Interactions of Sphagnum with water and air. In: Hutchinson TC, Meema KM, editors. Effects of atmospheric pollutants on forests, wetlands and agricultural ecosystems Springer, Berlin, Heidelberg; 1987. pp. 513–529. https://doi.org/10.1111/j.1469-8137.1987.tb00120.x pmid:33874034

33. Hanson, P.J., Phillips, J.R., Nettles, W.R., Pearson, K.J., Hook, L.A. 2020. SPRUCE Plot-Level Water Table Data Assessments for Absolute Elevations and Height with Respect to Mean Hollows Beginning in 2015. Oak Ridge National Laboratory, TES SFA, U.S. Department of Energy, Oak Ridge, Tennessee, U.S.A.

34. Steinweg JM, Kostka JE, Hanson PJ, Schadt CW. Temperature sensitivity of extracellular enzymes differs with peat depth but not with season in an ombrotrophic bog. Soil Biol Biochem. 2018; 125: 244–250.

35. Hill BH, Elonen CM, Jicha TM, Kolka RK, Lehto LL, Sebestyen SD, et al. Ecoenzymatic stoichiometry and microbial processing of organic matter in northern bogs and fens reveals a common P-limitation between peatland types. Biogeochemistry. 2014; 120(1): 203–244.

36. Hill BH, Elonen CM, Herlihy AT, Jicha TM, Serenbetz G. Microbial ecoenzyme stoichiometry, nutrient limitation, and organic matter decomposition in wetlands of the conterminous United States. Wetl Ecol Manag. 2018; 26(3):425–439. pmid:31073261

37. Lin X, Tfaily MM, Green SJ, Steinweg JM, Chanton P, Imvittaya A, et al. Microbial metabolic potential for carbon degradation and nutrient (nitrogen and phosphorus) acquisition in an ombrotrophic peatland. Appl Environ Microbiol. 2014; 80(11): 3531–3540. pmid:24682299

38. Lin X, Tfaily MM, Steinweg JM, Chanton P, Esson K, Yang ZK, et al. Microbial community stratification linked to utilization of carbohydrates and phosphorus limitation in a boreal peatland at Marcell Experimental Forest, Minnesota, USA. Appl Environ Microbiol. 2014; 80(11): 3518–3530. pmid:24682300

39. Keiser AD, Smith M, Bell S, Hofmockel KS. Peatland microbial community response to altered climate tempered by nutrient availability. Soil Biol Biochem. 2019; 137: 107561.