Ramped thermal analysis for isolating biologically meaningful soil organic matter fractions with distinct residence times

In this work, we assess whether or not ramped thermal oxidation coupled with determination of the radiocarbon content of the evolved CO2 can be used to isolate distinct thermal fractions of soil organic matter (SOM) along with direct information on the turnover rate of each thermal fraction. Using a 30-year time series of soil samples from a well-characterized agronomic trial, we found that the incorporation of the bomb spike in atmospheric CO2 into thermal fractions of increasing resistance to thermal decomposition could be successfully modeled. With increasing temperature, which is proportional to activation energy, the mean residence time of the thermal fractions increased from 10 to 400 years. Importantly, the first four of five thermal fractions appeared to be a mixture of fastand increasingly slower-cycling SOM. To further understand the composition of different thermal fractions, stepped pyrolysis–gas chromatography–mass spectrometry (Py-GC/MS) experiments were performed at five temperatures ranging from 330 to 735 C. The Py-GC/MS data showed a reproducible shift in the chemistry of pyrolysis products across the temperature gradient trending from polysaccharides and lipids at low temperature to ligninand microbe-derived compounds at middle temperatures to aromatic and unknown compounds at the highest temperatures. Integrating the 14C and Py-GC/MS data suggests the organic compounds, with the exception of aromatic moieties likely derived from wildfire, with centennial residence times are not more complex but may be protected from pyrolysis, and likely also from biological mineralization, by interactions with mineral surfaces.


Introduction
Soil organic matter (SOM) consists of a spectrum of material from labile, rapidly cycling compounds to mineral-stabilized molecules that resist degradation for centuries. This spectrum of turnover rates is due to a combination of organic matter composition with varying reactivity, various degrees of interaction between organic and mineral phases and greatly varying microclimates more or less suited to microbial activity (Lehmann and Kleber, 2015;Ruamps et al., 2013;Schmidt et al., 2011). Soil scientists often deal with this complexity by 10 fractionating SOM into what are thought to be more homogenous pools in terms of reactivity, composition or microbial accessibility depending on the particulars of the study (Christensen, 2001;von Lützow et al., 2007). Similarly, soil carbon cycle models typically divide SOM into conceptual pools 15 with distinct mean turnover rates (Manzoni and Porporato, 2009). However, empiricists have been trying for decades with varying degrees of success to link physically and chemically isolated fractions of SOM to the conceptual SOM pools in carbon cycle models (Skjemstad et al., 2004;Zimmermann 20 et al., 2007).
Fractionation schemes -employing various physical, chemical, biological or thermal methods -are generally used to reduce the inherent complexity found in SOM. Both size fractionation and density separation are commonly used to 25 separate particulate from mineral-associated SOM (Elliott and Cambardella, 1991;Golchin et al., 1994;Sollins et al., 2009). Hydrolysis with strong acid has been shown to isolate SOM that is consistently hundreds to thousands of years older than the bulk SOM (Paul et al., 2001). Biological frac-5 tionation involves the modeling of mass loss or CO 2 evolution during a laboratory incubation experiment (Grandy and Robertson, 2007;Schädel et al., 2013). Baldock et al. (2013), recognizing that fire-derived pryogenic carbon (PyC) has distinct properties from plant or microbial products, combined physical size fractionation with the use of solid-state 13 C nuclear magnetic resonance spectroscopy to mathematically isolate the PyC fraction from each size fraction. There are advantages and disadvantages to each of these techniques (Poeplau et al., 2018), but even the most detailed fractionation schemes are unable to isolate homogenous SOM pools (Jastrow et al., 1996;Sanderman et al., 2013;Torn et al., 2013), which may be as much a result of methodological issues as the fact that there are multiple pathways for SOM formation (Cotrufo et al., 2013;Sokol et al., 2019). 20 Thermal analysis techniques, long used in petrochemical exploration and clay mineralogy, offer a promising alternative or complement to physically and chemically based fractionation methods and are increasingly being applied to studies of SOM stability and loss (Peltre et al., 2013;Plante et al., 25 2009; . The basic premise of this suite of tools is that by slowly heating a sample the energy needed to evolve the carbon at different temperatures, whether that comes from breaking an organic-organic bond or disrupting an organic-mineral bond or other association, can be 30 quantified and that this energy yield is somewhat related to the energy requirements for enzymatic degradation of SOM . While research suggests that more complex organic molecules have higher thermal stability (Lopez-Capel et al., 2005;Yang et al., 2006), contra-35 dictory results have also been reported (Rovira et al., 2008). Several studies have also found good correlations between thermal stability indices and biological stability (Peltre et al., 2013;Soucémarianadin et al., 2018) and model-derived stable C pools (Cécillon et al., 2018). However, other studies 40 have found that new carbon preferentially flowed into more thermally stable fractions (Helfrich et al., 2010;Schiedung et al., 2017), suggesting that the relationship between thermal stability and SOM cycling concepts may not be straightforward. 45 A recent advance in thermal analysis is the coupling of a temperature-controlled oven to a vacuum line, termed ramped pyrolysis-oxidation (RPO), enabling collection of evolved gas at distinct temperature regions for subsequent 13 C and 14 C analysis (Rosenheim et al., 2008). Unlike inter-50 pretation of other thermal indices, 14 C is a direct and powerful tracer of soil carbon cycling by providing information on the age and turnover rate of an isolated fraction of SOM (Trumbore, 2009;Trumbore et al., 1996). In the first application of the RPO system to SOM, Plante et al. (2013) found 55 that more thermally stable fractions also contained the oldest, most 14 C-depleted carbon.
In this study, we present the first use of ramped thermal analysis with a time series of soil samples to investigate the carbon cycling rate of SOM with increasing thermal stabil-60 ity. The use of a multidecadal time series of soil samples allows for modeling the uptake of the bomb spike in atmospheric 14 CO 2 due to nuclear weapons testing into SOM fractions in order to determine the turnover rate of those fractions (Baisden et al., 2013). We further couple the 14 C-based 65 turnover time estimates with parallel chemical characterization of similar thermal fractions using a stepped pyrolysisgas chromatography-mass spectrometry (Py-GC/MS) approach. With the compound-specific chemistry data we are aiming to be able to start to explain possible mechanisms 70 controlling the turnover of the thermal fractions and thereby provide new insights into the linkages between thermal stability, SOM composition and microbial cycling of SOM.

Trial and soil description 75
The soil used in this study comes from a long-term agricultural research trial evaluating alternative crop rotations located at the Waite Research Institute in South Australia (34.967 • S, 138.634 • E TS1 ), with a Mediterranean climate, where 80 % of the 626 mm of average rainfall occurs dur-80 ing the April-October growing season. The soil is classified as a Rhodoxeralf in the USDA taxonomy (Soil Survey Staff, 1999) or a Chromic Luvisol in the World Reference Base taxonomy (IUSS Working Group WRB, 2015) with a fine sandy loam texture in the upper horizons, and 85 a mean pH (H 2 O) of 5.9 and clay content of 18 % for the 0-10 cm layer (Grace et al., 1995). Full trial history including management records, agronomic performance and soil data have been reported elsewhere (Grace et al., 1995;Sanderman et al., 2017). These data, including monthly climate records, 90 can be accessed from the Commonwealth Scientific and Industrial Research Organisation (CSIRO) Data Access Portal (https://doi.org/10.4225/08/55E5165EC0D29).
For this study, we have chosen to focus on the permanent pasture treatment because previous work (Sanderman et al., 95 2017) found the greatest uptake of the bomb spike in atmospheric 14 CO 2 in this treatment. This particular trial strip (#29) was under a wheat-pea rotation from 1925 until 1950 but then converted to an improved pasture by sowing a mix of annual rye grass, subterranean clover and Phalaris in 1950 100 and was then managed consistently with simulated grazing (i.e., hand mowing) and periodic re-sowing until the end of the trial in 1996. In April of 1963April of , 1973April of , 1983April of and 1993 soil samples from the top 10 cm were collected as a composite of 20 cores taken along the center line of this 90 m long strip 105 trial. Soils were dried at 40 • C for > 48 h before being stored in glass jars prior to subsampling in 2015. Previous analy-J. Sanderman and A. S. Grandy: Thermal fractionation of soil organic matter 3 ses (Sanderman et al., 2016(Sanderman et al., , 2017 suggest that SOM and the proportion in a particulate fraction only varied slightly throughout the 1963-1993 period ( Table 1). The large change in 14 C was due to uptake and loss of the bomb spike in atmospheric 14 CO 2 .

2.2 Ramped pyrolysis-oxidation
Ramped pyrolysis-oxidation (RPO) was performed using the Dirt Burner, a custom-built evolved gas system at the National Ocean Sciences Accelerator Mass Spectrometry (NOSAMS) facility, where a sample can be linearly heated 10 under either pyrolyzing or oxidizing conditions. The evolved gases are then oxidized to CO 2 , measured on an in-line infrared gas analyzer and trapped for subsequent analysis of 13 C and 14 C composition. Initial instrument development is described by Rosenheim et al. (2008) with upgrades de-15 scribed by Plante et al. (2013). See Hemingway et al. (2017b) for a complete description of current instrument configuration and operating conditions. In this investigation we operated the Dirt Burner only in oxidizing mode. It has been demonstrated that the distribution of activation energies and 20 14 C age of thermal fractions for soil are similar under oxidizing and pyrolysis modes (Grant et al., 2019). The 1973 soil sample was initially run using a fast ramp (20 • C min −1 ) where only CO 2 concentration was recorded. Using an inversion method (Hemingway et al., 2017a), five distinct thermal 25 fractions were identified from this thermogram with temperature ranges of 100-325, 325-400, 400-445, 445-515 and > 515 • C (Fig. A1). Subsequently, 45 mg of each soil sample were combusted on the Dirt Burner using a slow temperature ramp (5 • C min −1 ), and the evolved CO 2 in each of 30 these five fractions was cryogenically collected and purified for subsequent analyses. The RPO fractions were then split for 13 C analysis on a dual-inlet isotope ratio mass spectrometer (McNichol et al., 1994b) and 14 C composition via accelerator mass spectrometry after graphitization (McNichol et al., 35 1994a). Stable isotope data are expressed in δ 13 C notation (‰) relative to the Vienna Pee Dee Belemnite (VPDB) standard. Radiocarbon data, after using 13 C data to correct for mass-dependent fractionation, are reported from NOSAMS as fraction modern (Fm, where the 1950 atmosphere is as-40 signed a value of 1.0) and subsequently converted to the geochemical 14 C notation (Stuiver and Polach, 1977) for bomb-spike modeling.

Pyrolysis-gas chromatography-mass spectrometry
We determined the relative percentages and ratios of chem- 45 ical classes using pyrolysis-mass spectrometry-gas chromatography (Py-GC/MS) using methods described previously (Grandy et al., 2009;Kallenbach et al., 2015;Wickings et al., 2011). However, in contrast to our previous studies in which we used a single pyrolysis temperature, here we used 50 a "ramp" or stepped approach by pyrolyzing the same sam-ple at five distinct, sequentially increasing temperatures: 330, 396, 444, 503 and 735 • C (Buurman et al., 2007;Hempfling and Schulten, 1990;Williams et al., 2014). Thus, the same sample was pyrolyzed five times, corresponding with each of 55 these temperatures. Temperatures were chosen to correspond with the temperature ranges from the RPO analysis. Samples were pyrolyzed at each temperature and pyrolysis products transferred to a GC, where compounds were separated on a 60 m capillary column with a starting temperature of 40 • C 60 followed by a temperature ramp of 5 • C min −1 to 270 • C followed by a final ramp (30 • C min −1 ) to 300 • C. Compounds were immediately transferred from the GC to an ion trap MS, where they were ionized at 70 eV in the electronicionization mode with the source temperature held at 200 • C, 65 detected via electron multiplier, and identified using a compound library built using the National Institute of Standards and Technology (NIST) database and published literature as reported previously (Grandy et al., , 2009. Individual compounds were classified by their source as polysaccha-70 ride, aromatic, phenolic, protein, N-bearing (non-proteins) or unknown. Proteins include pyridines, pyrroles and indoles, which have been previously identified as pyrolysis products of proteins (Schulten and Schnitzer, 1997). The compounds that are in the unknown category are identified but can be 75 derived in nature or due to pyrolysis from different sources (e.g., both protein and aromatic).

Data analysis
The turnover time of each thermal fraction was determined by modeling the incorporation of the bomb spike in atmo-80 spheric 14 CO 2 into each fraction (Baisden et al., 2013). We have applied a steady-state soil carbon turnover model to each thermal fraction where carbon inputs (C in ) are portioned into each model pool proportional to the fractional distribution of carbon in each pool (f pool ) with the 14 C/ 12 C ra-85 tio of the previous year's atmospheric CO 2 (data from Currie et al., 2011). Carbon losses follow first-order kinetics with a characteristic decay rate for each pool (k pool ), and shifts in 14 C/ 12 C ratio are also affected by radioactive decay (λ = 1.21 × 10 −4 yr −1 ). Turnover time (τ pool ) is simply 90 the inverse of the decay rate for each pool. First, we assumed that each thermal fraction is a single homogenous pool and solve the turnover model to find a single k value by minimizing the sum of square errors (SSE) between observed and predicted 14 C data for the four years of data (1963, 1973, 95 1983 and 1993) for that particular thermal fraction. Initial results suggested that most of the thermal fractions could not be represented as a single homogenous pool, so we then applied a two-pool model by assuming a fixed k for the fastcycling pool (k fast = 0.25 yr −1 ) and allowed the size (f slow ) 100 and decay rate (k slow ) of the slow-cycling pool to vary while minimizing the SSE between observed and predicted 14 C data for each thermal fraction. Given the limited number of degrees of freedom, we decided not to allow for the simulwww.soil-journal.net/6/1/2020/ SOIL, 6, 1-14, 2020 taneous optimization of all three parameters. The value for k fast was set at 0.25 yr −1 after performing a simple sensitivity analysis where k fast was varied from 0.1 to 1.0 yr −1 and f slow and k slow were optimized. This exercise suggested that the overall lowest mean root mean square error (RMSE) across 5 thermal fractions was achieved with k fast = 0.25 yr −1 . Model performance was assessed by calculating the RMSE between observed and predicted 14 C values. The stepped Py-GC/MS data are reported as percent relative abundance for each identified compound. We present 10 these data qualitatively in two ways. First, data from each thermal interval were averaged across years, and shifts in the eight major compound classes are shown. Second, all compounds that averaged > 1 % abundance across all 20 samples (4 years × 5 thermal intervals) were included in a non-metric 15 multidimensional scaling analysis (e.g., Grandy et al., 2009) after constructing a resemblance matrix using Euclidean distance from the percent abundance data.

RPO results 20
All samples produced nearly identical thermograms (Fig. 1), suggesting that the activation energies assigned to the 1973 sample ( Fig. A1) are applicable to samples from all four decades. Fractions 2 and 3 contained the greatest proportion of SOM, with the distribution across all five fractions being 25 consistent between across all years (Table B1). A greater incorporation of the bomb spike in 14 C was found in the lowertemperature fractions, with no modern carbon seen in the highest-temperature fraction (Fig. 2). There was not enough CO 2 for 14 C analysis of F5 in 1963. The 14 C value for F1 30 in 1973 is included in Fig. 2 but removed from subsequent bomb-spike modeling.

Bomb-spike turnover modeling results
The single-pool steady-state model to trace the incorporation of the bomb spike in atmospheric 14 CO 2 into the ther- 35 mal fractions indicated that turnover times increased from 37 to 386 years with increasing temperature (Table 2). How- Figure 1. Proportion of total CO 2 evolved with temperature as samples were oxidized with a ramp rate of 5 • C min −1 . Thermal fractions are shaded with colors corresponding to Figs. 2-4, and mean percent of total is given for each fraction.

Figure 2.
Variation in 14 C in the five thermal fractions across the four decades of soil sampling. There was not enough CO 2 for 14 C analysis of F5 in 1963. 14 C value from F1 from 1973 not shown in (b) as this data point was removed as an outlier from the turnover modeling. ever, with the exception of F5 a single-pool model resulted in unacceptably high error, with RMSE values ranging from 25 ‰ to 45 ‰ for the other fractions. Visually it was clear 40 that the single-pool solution could not capture the dynamics of the increase and subsequent decrease in 14 C in F1-F4 (Fig. C1). A two-pool model with a fast-cycling pool (τ = 4 years) and a variably sized slower-cycling pool was generally able to capture these dynamics with RMSE values below 13.8 ‰ (Table 2). With increasing activation energy of the thermal fraction, the proportion of fast-cycling carbon decreased (Fig. 3a) and the turnover time of the slow-cycling fraction increased (Fig. 3b).
As an independent first-order check on the reasonableness of CE1 MRT results for the thermal fractions, the inventoryweighted MRT of the bulk soil was 16.5 years, matching the 10 results of applying a three-pool turnover model to the bulk 14 C data, which varied from 11 to 20 years depending on the model structure (Sanderman et al., 2017).

Py-GC/MS results
A total of 172 individual compounds were identified in one or 15 more samples. Many of these compounds were a minor fraction of only a single sample. When these compounds were classified by source, there were sharp differences in the dominant chemical compound classes released at each pyrolysis temperature (Fig. 4a). Lipids (42.5 ± 2.74 %) and polysac-20 charides (42.4 ± 2.12 %) were the dominant classes of compounds released at the lowest pyrolysis temperature (330 • C). Relative lipid abundance was < 10 % of the pyrolysis products for all higher temperatures and < 1 % relative abundance at the highest pyrolysis temperature (735 • C). Polysaccha-25 rides (51.2 ± 1.68 %) remained an abundant component of the pyrolysis products at a pyrolysis temperature of 396 • C and 444 • C (27.4 %) but were < 10 % relative abundance at the two highest temperatures. At the fourth temperature level (503 • C), N-bearing compounds, proteins, phenols and 30 unknown compounds dominated the chemical signature and were all ∼ 15 %-20 % relative abundance. Phenols (31.4 %) were the most abundant compound class released at the highest temperature (735 • C), while proteins, aromatics and compounds of unknown origin were all > 15 % relative abun-35 dance.
We identified 23 individual compounds that averaged more than one 1 % relative abundance across all samples (Fig. 4b). These 23 compounds collectively represented from about ∼ 70 % to 90 % of the total sample relative abun-40 dance at the five different temperatures (Table D1). These compounds included three aromatic compounds, one lignin derivative, one lipid, two N-bearing compounds, two phenols, five polysaccharides, five proteins and four compounds that fell into the unknown origin. The most abundant com- 45 pound was hexadecanoic acid methyl ester (palmitic acid) with a mean relative abundance across all samples of about 10 % but composing 42 % of F1 (Table D1). The next two most abundant compounds were both phenols and included phenol and 4-methyl phenol, both increasing with increas-50 ing temperature, followed by the common polysaccharide pyrolysis product furfural, primarily found in the lowertemperature fractions.

Discussion
While this study has only examined one soil under consistent 55 management over four decades, the combined results from the modeling of the incorporation of bomb-spike 14 CO 2 into thermal fractions and the stepped Py-GC/MS analysis here suggest that evolved gas analyses can be powerful analytical tools for understanding the complexities of SOM cycling. 60 With increasing activation energy, turnover times increased from decades to centuries (Fig. 3), and there was a consistent (i.e., repeatable across time) strong shift in SOM chemistry (Fig. 4).
This investigation was partially framed in the context of 65 testing evolved gas analysis as a tool for isolating biologically meaningful fractions of SOM. The results indicate that each of the five thermal fractions contained unique information. Each of the thermal fractions represented pools with a diverse mixture of organic materials with 14 C-based 70 turnover times ranging from 10 to 400 years. A wide distribution of activation energy (Fig. A1) was needed to describe each thermal fraction. Additionally, a one-pool turnover model could not represent the dynamics of 14 C in all but the most thermally stable fraction (Table 3). This finding 75 of heterogeneity within an isolated fraction also plagues other physical and chemical techniques for isolating carbon fractions Torn et al., 2013). Here, by combining 14 C-based turnover modeling and compoundspecific chemical characterization of the thermal fractions, 80 we have found clues as to reasons for the heterogeneity of SOM within these fractions which provide insights into pathways of SOM formation and longer-term stabilization.
At the two lowest-temperature intervals, where the 14 C data suggested rapid turnover rates, the polysaccharides com-85 prised > 40 % of the identified compounds. The dominant polysaccharide products at the two lowest temperatures included furfural, 3-furaldehyde and levoglucosenone. Furfural is a pyrolysis product of hexoses and pentoses as well as uronic acids. The hexoses originate from both plant and mi-90 crobial residues, while the pentoses are primarily microbial in origin. Levoglucosenone is a pyrolysis product of neutral sugars such as glucose, galactose and mannose (Saiz-Jimenez et al., 1979;Saiz-Jimenez and De Leeuw, 1986). These polysaccharides may be part of partially decomposed 95 residues that dominate the light or particulate fractions of SOM or may be part of other soil carbon pools not strongly protected from thermal degradation by association with minerals. Across all thermal fractions, polysaccharide abundance was strongly correlated with the proportion of fast-cycling 100 carbon in that fraction.
Identifiable lignin derivatives (e.g., syringol and guaiacol) were most abundant in the second and third thermal fractions but never constituted more than 10 % of either thermal fraction. While the chemical complexity of lignin likely affords 105 some thermal stability, as suggested by the greater activation energy than the thermal fraction dominated by polysaccha-Please note the remarks at the end of the manuscript.
www.soil-journal.net/6/1/2020/ SOIL, 6, 1-14, 2020  rides, the 14 C data suggest that lignin is not a particularly persistent compound in this soil. This finding is consistent with past studies' findings that there is little recognizable lignin found in the more stable clay fractions of many soils (Baldock et al., 1997;Grandy and Neff, 2008;Schulten and 5 Leinweber, 2000), although amounts may vary across different mineral types and environmental contexts (Kramer et al., 2012). Besides polysaccharides, hexadecanoic acid methyl ester (palmitic acid) was the other dominant constituent of the 10 lowest-temperature fraction (Table D1). This compound is a major constituent of microbial cell walls, especially fungi, and has been strongly correlated to both microbial biomass and activity (Zelles et al., 1992). The relative dominance of this lipid at low pyrolysis temperatures suggests that com-15 pared to other compounds it is not protected from degradation by association with either minerals or chemical complexity. There is a growing body of evidence suggesting that a substantial portion of stable SOM in many soils is microbial in origin (Kallenbach et al., 2015;Knicker, 2011;20 Miltner et al., 2012), which could lead to the misconception that microbial products are inherently more stable than plantderived compounds. The findings here suggest that palmitic acid, the most abundant lipid that we identified, which is most likely of microbial origin, is not a very thermally sta-25 ble component in soil with a much faster than average MRT. This may arise because these lipids are not forming direct covalent bonds with mineral surfaces nor are they forming other complexes (e.g., with metal oxides) that could make them resistant to thermal and biological degradation.

30
At the third-and fourth-highest temperatures, phenols, proteins and N-bearing compounds became more relatively abundant. These phenols can originate not only from degraded lignin monomers but also other aromatics of plant or microbial origin. The proteins and other N-bearing com-35 pounds are important components of mineral-associated N. The presence of heterocyclic N (e.g., furans, pyridines and pyrroles) is a common finding of studies using Py-GC/MS (Schulten and Schnitzer, 1997;Leinweber et al., 2014). Secondary chemical reactions in soil (de Assis et al., 2012) and 40 wildfire (González-Pérez et al., 2004) as well as the pyrolysis process itself (Hatcher et al., 2001) are potential sources of heterocyclic N, but direct plant and microbial inputs also make a key contribution to this pool (Leinweber et al., 2014;Paul, 2016). The pyrazoles and pyridines, for example, may 45 largely originate from microbes, and Kallenbach et al. (2016) found in an artificial soil system with only glucose inputs that a variety of heterocyclic N compounds can be derived from the pyrolysis of microbial cells.
We have previously found that polysaccharides, proteins and N-bearing compounds are the dominant mineralassociated chemical fractions Grandy and Neff, 2008) as well as comprising an increasing proportion of these compounds with soil depth (Rothstein et al., 5 2018). Anticipating that these same compounds would also dominate the slowest-cycling highest-temperature fraction, we found that aromatic compounds and compounds from an unknown source, which included compounds such as toluene, the eleventh-most-abundant compound, which may 10 be derived from pyrolyzing aromatics or proteins, were relatively abundant. There have been reports that pyrolysis of clay fractions can produce matrix effects that inflate the abundance of aromatic compounds (Schulten and Leinweber, 1993). We have seen no detectable evidence for this in our 15 previous studies, which generally show very low abundance of aromatics in fine fractions. The relatively low clay content and dominance of illite in our soils also point to lower potential for matrix effects. More likely, the high relative abundance of stable aromatics, phenols and compounds of un-20 known origin can be attributed to the history of fire at this Mediterranean site. Previous studies at the site indicate that ∼ 30 % of the total SOM is made up of pyrogenic aromatics (Sanderman et al., 2017).
A potential limitation of this study was that we could not 25 determine if the thermal fractions isolated for 14 C analysis corresponded precisely with the stepped Py-GC/MS data. Grant et al. (2019) demonstrated that for soil samples there was minimal difference in the shape of evolved CO 2 between ramped oxidation and ramped pyrolysis on the NOSAMS 30 RPO instrumentation, but it is difficult to quantify the carbon mass evolved through the Py-GC/MS. Future research coupling the stepped pyrolysis front end with an oxidizing oven and subsequent CO 2 measurement is necessary to alleviate this uncertainty. 35 The findings from this preliminary investigation add to the growing body of literature using evolved gas analysis as a tool for understanding SOM dynamics. The two unique aspects of this work were using a time series of soil samples to calculate the MRT of different thermal fractions and relat-40 ing these MRTs to the chemical composition of each fraction. We found that mean residence time increased with increasing activation energy. Modeling the incorporation of the bomb spike in 14 CO 2 indicated that the thermal fractions were, except for the most stable fraction, heterogeneous mixtures 45 of fast-and slow-cycling SOM. Compound-specific analysis demonstrated that distinct assemblages of organic compounds were found with increasing thermal stability. These findings together suggest similar thermal activation energies may not equate to similar biological accessibility of the same 50 material, and care needs to be taken not to overinterpret results from thermal analysis alone. Important follow-up work includes analyzing just the mineral-associated OM fraction to better understand this potential disconnect between thermal and biological stability. By coupling evolved gas analy-55 sis with radiocarbon and compound-specific chemical analyses, new insights into the formation, stabilization and fate of SOM may be possible.

Appendix A
Following acquisition of an initial thermogram (i.e., the trend of CO 2 evolution versus increasing temperature) for the 1973 sample, a mathematical inversion method was employed to deconvolve the evolved gas analysis into N number of pools 5 with distinct activation energy profiles (Hemingway et al., 2017a).

Appendix C
Bomb-spike turnover modeling results for one-pool and twopool models. Figure C1. One-pool (a, c, e, g, i) and two-pool (b, d, f, h, j) solutions to a steady-state soil carbon turnover model for the five thermal fractions. Southern Hemisphere atmospheric record given in the background of each panel. For two-pool solutions the trend in 14 C is given for the fast and slow pools as well as the overall thermal fraction. The X in panels (a) and (b) refers to the 14 C measurement deemed to be an outlier for purposes of bomb-spike modeling.
the Dirt Burner and Jordon Hemingway for running the activation energy inversion model. Finally, we thank the numerous researchers at the Waite Institute for having the foresight to archive soil samples from the long-term experiment.
Review statement. This paper was edited by Steven Sleutel and 20 reviewed by two anonymous referees.