Physical Processes Dictate Early Biogeochemical Dynamics of Soil Pyrogenic Organic Matter in a Subtropical Forest Ecosystem

Quantifying links between py OM dynamics, environmental factors and processes is central to predicting ecosystem function and response to future perturbations. In this study, changes in carbon ( TC ), nitrogen ( TN ), pH, and relative recalcitrance ( R 50 ) for pine-and cordgrass-derived py OM were measured at 3–6 weeks intervals throughout the ﬁrst year of burial in the soil. Objectives were to (1) identify key environmental factors and processes driving early-stage py OM dynamics, and (2) develop quantitative relationships between environmental factors and observed changes in py OM properties. The study was conducted in sandy soils of a forested ecosystem within the Longleaf pine range of the United States with a focus on links between changes in py OM properties, ﬁre history ( FH ), cumulative precipitation ( P cum ), average temperature ( T avg ) and soil residence time ( SRT ). P cum , SRT and T avg were the main factors controlling TC and TN accounting for 77–91% and 64–96% of their respective variability. Fire history, along with P cum , SRT and T avg , exhibited signiﬁcant controlling effects on py OM pH and R 50 —accounting for 48–91% and 88–93% of respective variability. Volatilization of volatiles and leaching of water-soluble components (in summer) and the sorption of exogenous organic matter (fall through spring) were most plausibly controlling py OM dynamics in this study. Overall, our results point to climatic and land management factors and physicochemical process as the main drivers of py OM dynamics in the pine ecosystems of the Southeastern US.


INTRODUCTION
Current evidence suggests that despite a relatively small annual global production rate (40-259 Gt/yr; Bird et al., 2015;Santin et al., 2015), carbon in the fire-derived/pyrogenic organic matter (pyOM) pool accounts for a significant quantity (2-45%) of the total organic carbon in terrestrial systems (Skjernstad et al., 1999;Skjemstad et al., 2002;Lehmann et al., 2008;Reisser et al., 2016). There is fairly widespread consensus that the significant quantity of terrestrial carbon held in the pyOM pool is largely attributable to a comparatively higher innate resistance/recalcitrance to abiotic/biotic degradation and, longer turnover times of pyrogenic carbon compared to its non-pyrogenic counterpart (Schmidt and Noack, 2000;Hammes et al., 2008). For example, turnover times for pyrogenic carbon are on the order of decades to thousands of years, compared to a few years for its nonpyrolysed counterpart (Hammes et al., 2008;Liang et al., 2008;Carvalhais et al., 2014;de la Rosa et al., 2018).
The recognition of pyOM being relatively more resistant to abiotic and biotic transformation has led to increased interest in the global pyOM pool as a potential carbon sink and carbon sequestration strategy (Lehmann et al., 2006). Current estimates for carbon sequestration potential in the pyOM pool are as high as 9.5 × 10 9 tons C/yr-comparable to estimated carbon emission from fossil fuel burning (9.1 × 10 9 tons C/yr; Lehmann et al., 2006). However, it is now well-known that this sequestration potential will vary with pyOM composition/properties as well as field/environmental conditions. The composition/properties of pyOM have been extensively studied and are known to be a heterogeneous mix of thermally-altered carbon-, hydrogen-, nitrogen-, oxygen-, and sulfur-containing organic structures with the degrees of heterogeneity or thermal alteration varying systematically across feedstock chemistry and pyrolysis conditions (e.g., Zimmerman, 2010;Harvey et al., 2012). Increasing degrees of thermal alteration favor increasingly condensed carbonand/or nitrogen-rich structures of higher innate environmental recalcitrance/stability (Knicker et al., 2008;Knicker, 2010;Zimmerman, 2010;Harvey et al., 2012).
Current understanding of the impact of variations in environmental conditions on pyOM behavior is far less clear than that for feedstock chemistry or pyrolysis conditions. Results on the biogeochemical trajectory of key pyOM properties and, most importantly, the environmental factors driving that trajectory are still unclear. For example, Mukherjee et al. (2014) reported increases of up to 124, 143, and 43% in pyOM-associated carbon, nitrogen and exchange capacity, respectively-after burying different pyOMs for 15 months in two distinct soils (a sandy Entisol used for agriculture and a forested Spodosol) in Florida. The observed differences in pyOM behavior were attributed to differences in soil composition and structure. In contrast, while Sorrenti et al. (2016) and de la Rosa et al. (2018) observed increases in pyOM-associated nitrogen (fourfold) after soil burial of pyOM for 24-48 months, they both reported decreases in pyOM-associated carbon (11-27%). Both studies were conducted in inceptisols (in the case of Sorrenti et al. (2016)-a sandy loam in Italy and in the case of de la Rosa et al.-a sandy clay loam in Spain). Except noting that hydrologic properties of pyOM were linked to surface changes induced by environmental exposure, Sorrenti et al. (2016) did not discuss specific environmental factors driving observed changes. On the other hand, de la Rosa et al. (2018) noted climatic conditions (in particular rainfall and temperature) and soil characteristics as the key drivers of the observed changes in pyOM-associated carbon and nitrogen.
Whilst it is intuitive-and of great qualitative value-that environmental factors such as ecosystem type, soil properties, temperature, nutrient availability and moisture conditions drives pyOM dynamics (Cheng et al., 2008;O'Neill et al., 2009;Khodadad et al., 2011;Bird et al., 2015); quantitative assessments of the individual and interactive contributions of these factors are necessary for predicting the behavior and ecosystem value of pyOM over time and space. Such assessments remain sparse and, to our knowledge, are currently exclusively obtained via metaanalysis type studies. Recent meta-analytical studies directly addresses important environmental factors driving; (1) the environmental persistence of pyOM, its priming effects on native soil organic matter (SOM; Wang et al., 2016) and (2) its global distribution in carbon stocks (Reisser et al., 2016). Wang et al. (2016) found that for a given pyOM, experimental duration (soil residence time, SRT) and soil clay content controlled the degradation rate of pyOM, and its priming effect-with higher pyOM degradation rates and positive pyOM priming of native SOM favored at SRT < 0.5 years and soil clay contents ≤9%. Reisser et al. (2016) found that, globally, the largest pyOMassociated carbon stocks generally occurred in areas of warmer climates (≥10 • C) and soils having high clay contents (>50%) and/or pH > 7. Higher pyOM-associated carbon in warmer climates was attributed to higher temperature-promoted biomass production and fire probability; and in high clay content, high pH soils due to protection via organo-mineral interactions leading to less pyOM degradation (Reisser et al., 2016).
As demonstrated through the work of both Wang et al. (2016) and Reisser et al. (2016), global meta-analysis-type studies can provide valuable insights into the environmental factors that drive pyOM dynamics at the global scale. It is also equally important to recognize that such studies represent the global/onaverage case and are not necessarily directly transferrable (or downscalable) to explaining pyOM dynamics at local or regional scales. That is, factors determined to be driving pyOM dynamics at the global scale may not necessarily explain variability at finer scales-especially, when spatial and/or temporal variability in potential environmental drivers differ from that captured in global assessment.
It is reasonable to expect spatial variability in the (i) fire and landuse, (ii) climatic, and (iii) pedogenic drivers of pyOM dynamics, identified by Reisser et al. (2016), at the local and regional scales. With season-to-season variability in rainfall and temperature, changes in landuse and land/fire management, as well as projected climate shifts, environmental drivers of pyOM dynamics can also be expected to change on short-and longterm temporal scales. For example, with climate change projected to induce regional shifts in temperatures and precipitation, the fire regimes, the amount of pyOM produced, the soil moisture, biota, and consequently pyOM biogeochemical cycling are expected to be altered (Flannigan et al., 2000;Santin et al., 2015). Identifying region specific drivers of pyOM dynamics and developing quantitative relationships between them and the properties and processes they drive are central to predicting impacts of environmental change on ecosystem functions such as carbon sequestration and nutrient cycling.
In this study, we tracked early dynamics of pyOM from two isotopically distinct sources at 3-6 weeks intervals throughout the first year of burial in coarse textured ultisols located in a forested ecosystem within the Longleaf pine range of the Southeastern United States. The objectives of the study were to (1) identify key environmental factors and processes driving early-stage biogeochemical alteration in pyOM in such an ecosystem, and (2) develop quantitative relationships between environmental factors and observed changes in pyOM properties to explain early-stage dynamics of pyOM. Focus was placed specifically on the relationship between fire history (FH; burnt vs. unburnt), cumulative precipitation (P cum ), average temperature (T avg ), soil residence time (SRT) and early-stage temporal variability in pyOM-associated carbon, pyOM-associated nitrogen, pyOM pH and pyOM stability/degradability. The underlying hypothesis was that in its early stages, pyOM dynamics in sandy-textured (low CEC, low organic matter) soils is analogous to initial stages of soil formation whereby climatic effects (e.g., precipitation and temperature) and the innate properties of the parent material (rock type in the case of soils; in our case, pine-vs. cordgrassderived pyOM) dictates outcome.

Study Site, Sampling, and Soil Description
The study was conducted at the Lake Thoreau Environmental Center (LTEC; 31.3466 • N, 89.4223 • W), a 121-hectare facility owned and operated by The University of Southern Mississippi in a manner that parallels historic and current forest management practices across the Southeastern US. The LTEC is located within the pine-belt of Mississippi and the wider Longleaf pine range of the Southeastern United States. It is representative of the region in terms of vegetative cover (∼52% pine forest), fire management (2 year prescribed-burning cycle vs. no burning), surface soils (coarse textured ultisols formed from coastal plain sediments; Soil Survey Staff, 2018), mean annual temperature (15-21 • C; Markewitz et al., 2002), mean annual rainfall/precipitation (1,020-1,520 mm; Markewitz et al., 2002) and Koppen-Geiger Climate Classification (humid subtropical; Kottek et al., 2006). Average seasonal temperatures at the LTEC are 32.1, 25.5, 16.2, and 25.2 • C for the summer, fall, winter, and spring seasons, respectively with average monthly precipitation values for the same seasons of 112, 105, 139, and 147 mm, respectively.
A total of six sampling locations, along two converging transects, were selected for use in this study (Figure 1). Four sampling locations were in an area (herein referred to as the Burnt Zone) with a consistent two-year prescribed burning cycle since the mid 1960's; the most recent burn occurring in January of 2014-prior to the start of the study in the Summer of 2014 (May 31, 2014). The remaining two sampling locations were in an area (herein referred to as the unburnt zone) that has no history of ever being burnt or systematically managed in the last 50 years-besides infrequent grazing from miniature cattle in a small subsection. Sampling locations in the burnt zone are designated with a "B" and those in the unburnt zone are designated with a "U" (see Figure 1 for complete sampling designations). At the time of the study, the sample sites within the unburnt zone were not being grazed and had very heavily vegetated understories and shrub layers as well as thick litter layers on the forest floor.
Elevation for the study area varied from 88 to 100 m above sea level with overall downhill sloping from NE to SW (i.e., from B4 to B1). Soils at the study site are alluvial in origin and with the exception of B2, were mapped in the US Soil Survey as being from the Freestone-McLaurin-Susquehanna association. The Freestone-McLaurin-Susquehanna association is characterized by a moderately well-drained upper profiletypical of the predominantly coarse-textured ultisols found in the Long leaf pine range. Soils at sampling location B2 were mapped as the Prentiss series (also typical of the region) and comprising also of moderately well-drained upper profile. The textures of soils in the upper 10 cm at the sample locations were determined via the hydrometer method of Gee and Bauder (1986) to be sandy clay loam (B1, B2, B3, U1, U2) or borderline sandy clay loam/clay loam (B4; Supplementary Figure 1).
At each of the six sampling location, samples of pinederived (P, Pinus taeda) and cordgrass-derived (CG, Spartina alterniflora) pyOM were placed in polypropylene litterbags and buried in a gridded design at a depth of ∼8 cm below the soil surface (Supplementary Figure 2). A total of 11 P and 11 CG pyOM-containing litterbags were buried at each sampling location. One litterbag each of P and CG pyOM were removed on each sampling excursion, which occurred at 21,42,63,84,105,126,147,217,294,322, and 350 days after burial. Samples collected on each excursion were split 2:1, with two-thirds of each sample being dried (75 • C for 24 h) and used for immediate analysis of pH, electrical conductivity (EC), relative recalcitrance (R 50 ), total carbon content and total nitrogen content, while one-third of the sample was frozen and subsampled for use in ergosterol, and δ 13 C and δ 15 N isotopic analysis.

Preparation of pyOM
Pine-derived pyOM was produced from pine needles collected from the study site, while cordgrass-derived pyOM was produced from cordgrass straw collected from USM's Gulf Coast Research Laboratory. Prior to pyrolysis, needles and grass were dried at 60 • C to constant weight and cut into ∼2 cm pieces. Pyrolysis of plant material was performed in a muffle furnace (ramp rate = 25 • C/min) under oxygen-limited conditions at 450 • C for a total of 1 h. The 450 • C was representative of average temperature conditions for a typical vegetation fire in the southeastern United States (Alexis et al., 2007). Upon cooling and removal from the furnace, the pyOM material was separated into a fine (<250 µm) and a coarse (250-500 µm) fraction. Approximately 1.5 g of the coarse pyOM material was placed into litterbags (8 × 8 cm). The litterbags used were made from needle-punched polypropylene material with apparent openings of 210 µm (small enough to contain the coarse fraction of the pyOM, yet large enough to allow environmental interaction with the pyOM).

Post-burial pyOM and Soil Analysis
Pyrogenic organic matter collected on each sampling excursion was analyzed for total carbon content (TC), total nitrogen (TN), pH and relative recalcitrance (R 50 ). Samples from selected excursions were analyzed for fungal ergosterol, base-extractable organic carbon (BEOC) and δ 13 C and δ 15 N isotopic signature.
Carbon and nitrogen content was determined in triplicate by combustion elemental analysis (Costech elemental combustion analyzer) using 2-3 mg of pyOM tightly wrapped in tin capsules combusted and compared to acetanilide as a standard (Sarkhot et al., 2012).
Base-extractable organic carbon (BEOC) of pyOM samples was determined using UV-visible spectrophotometry. Between 150 and 151 mg of pyOM was placed in 30 ml of 0.5M NaOH solution before being heated for 2 h at 80 • C. After heating, the pyOM: NaOH suspension was cooled, filtered (0.2 µm) and the absorbance measured at 365 nm then compared to a 4point standard curve of absorbance vs. dissolved organic carbon content from heated-persulfate oxidation method (Shimadzu Scientific Instruments, Columbia, Maryland; Raimbault et al., 1999;Harvey et al., 2009).
For pyOM pH and EC, 5 ml of de-ionized water was mixed with 50 mg of sample, creating a 1:100 pyOM: water suspension. This solution was allowed to equilibrate for 12 h before measuring the pH and EC of the extractant (Accumet pH/EC submersible probe).
Relative recalcitrance of the pyOM was determined by thermogravimetric analysis (TGA; SDTQ600, TA Instruments, New Castle, DE) and the use of the R 50 index as described by Harvey et al. (2012). Thermal analysis began at an oven temperature of 20 • C, increasing at a ramp rate of 10 • C per minute until 700 • C, at which point no further weight loss was observed. The R 50 value of the pyOM is a measure of the innate resistance of the pyOM material to abiotic and biotic oxidation, and was calculated as: where T 50,pyOM and T 50,graphite are the temperature values corresponding to 50% oxidation of pyOM material and graphite (<149 µm, purity 99.9995%; Alfa Aesar, Ward Hill, MA), respectively, and were obtained via use of TGA mentioned above.
The fungal biomass associated with pyOM samples were estimated using ergosterol concentrations (Kuehn et al., 2011). Ergosterol was extracted from pyOM samples (∼60 mg) by refluxing in alcoholic KOH (4% KOH in 95% methanol) for 30 min. The resultant extract was partitioned into n-pentane and evaporated to dryness under a stream of nitrogen gas. Ergosterol in dried samples was then re-dissolved by sonication in 0.5 ml of methanol and stored, tightly capped, in 1.5 ml screw cap high pressure liquid chromatography (HPLC) vials at ∼20 • C until analyzed. Separation and analysis of ergosterol was done using HPLC according to procedures from Kuehn et al. (2011).
Carbon and nitrogen (δ 13 C and δ 15 N) isotopic signature of the pyOM were determined by EA-IRMS using 2-3 mg of sample and an elemental analyzer (EA) interfaced to a Thermo-Electron Delta V Advantage isotope ratio mass spectrometer (IRMS). Values for δ 13 C were reported relative to the Vienna PeeDee Belemnite standard and values for δ 15 N to atmospheric nitrogen.
In addition to pyOM, soils were collected (to a depth of ∼3 cm) from directly beneath the litterbags on each of the following sampling excursions. The collected soil samples were air-dried to constant weight, grounded to <2 mm and used for soil pH and soil organic carbon (SOC) analyses. Soil pH and SOC was determined using the same procedures as described above for pyOM albeit samples for SOC analyses were done using 20 mg of sample. It is worth noting that in these soils TC = SOC because of the lack of inorganic carbon-verified by the lack of fizzing in response to HCl additions and supported by the acidic soil pH.

Precipitation and Temperature Data
The LTEC maintains a weather station (∼600 m west of the study site) with research grade sensors for monitoring precipitation (Tipping bucket rain gauge; model: TE525WS-Texas Electronics) and temperature (HygroClip2 HC2s3-Rotronic Instrument) as well as other climatic parameters (relative humidity, barometric pressure, solar irradiance and wind speed/direction) on 1-min temporal scales. For this study, 1-min averaged precipitation and temperature data were aggregated to obtain daily precipitation and average air temperatures for every day of the study. Daily values were further aggregated through addition or averaging to obtain cumulative precipitation (P cum ) and average maximum daily temperature (T avg ), respectively for each sampling day. It is these latter P cum and T avg that are utilized in all our analyses.

Data Analysis and Modeling
Data analysis was conducted in two phases. The first phase of analysis was executed using two-way analysis of variance and was focused on assessing the relative importance of time (as SRT) and sample site in driving observed changes in pyOM properties of both pine-and cordgrass-derived pyOM within a given fire history/zone (i.e., burnt vs. unburnt). The goal of this initial analysis was to screen if measured changes in pyOM properties at the study site were controlled primarily by temporally-or spatially-driven environmental factors (e.g., soil properties) and, therefore, data within a given fire zone could be aggregated across space or time. Two-way analysis was done using Graphpad Prism (version 7.0).
The second phase of data analysis was executed using generalized linear modeling (GLM) and was focused on identifying and then modeling statistically significant (p < 0.05) main and/or interactive effects of environmental co-factors (soil residence time, SRT (in days); cumulative precipitation, P cum (in mm); average maximum daily temperature, T avg (in • C); and fire history, FH; burnt vs. unburnt) on changes in pyOM properties (TC, TN, pH, and R 50 ). Three initial GLM models were fitted to respective datasets for cordgrass-and pine-derived pyOM: 1) A main-effects only model with a linear predictor component of the form, Y = β 0 + β 1 (SRT) + β 2 (P cum ) + β 3 (T avg ) + β 4 (FH) + ε; 2) An interaction-effects only model with a linear predictor component of the form, Y = β 0 + β 1 (SRT×P cum ) + β 2 (SRT×T avg ) + β 3 (SRT×FH) + β 4 (P cum ×T avg ) + β 5 (P cum ×FH) + β 6 (T avg ×FH) + ε; and 3) A combined-effects model with a linear predictor component of the form, where, Y is the pyOM property of interest, β 0 is the model intercept; β 1... β n are fitted values for contributions of each predictor to Y and ε is the error term. Values for FH were set to 1 for burnt zone and 0 for the unburnt zone. Each model was fitted in two iterations with a normal random component and an identity link function. The first iteration involved all possible predictors (as outlined above) while the second iteration involved only predictors that were identified from the first iteration as having significant β n values (i.e., p < 0.05). Maineffects, interaction-effects and combined-effects models from the second iteration were compared via ANOVA; with the best model being selected as that with the lowest AIC (Akiake Information Criterion) and best fit/trace (R 2 ) to the measured data. All statistical analyses were done using the GLM function from statistical package in R-Gui (R-Core Team, 2017).

Overall Effects of Time and Space on pyOM Dynamics
Initial properties of the pine-and cordgrass derived pyOM are summarized in Table 1. Two-way analysis of variance examining the effect of time and sampling site on changes in pyOM properties indicated that time (SRT) accounted for 80-92% of total variability in pyOM-associated organic carbon, 38-79% in pyOM-associated nitrogen and 81-91% in pyOM pH over the course of the study. By comparison, sampling site accounted for 0.005-2.53, 0.08-13, and 1-3% of variability in pyOM-associated carbon, nitrogen and pH respectively. Interactive effects between time and sampling site were also much (4-49 times) lower than time effects alone in all cases; except for TN in the unburnt zone where the effect of time alone and interactive effects were equal (Supplementary Table 1). The much higher variability in pyOM-associated carbon, nitrogen and pH accounted for by time (compared to location or interactive effects) pointed to temporally variable environmental, rather than location-based, factors being the primary drivers of pyOM dynamics in this study. This also indicated that in contrast to other studies (Kane et al., 2010;Kasin and Ohlson, 2013), variations in early-stage pyOM dynamics due to location-based parameters such as elevation, soil type and position in landscape would not be important in our study. Of such, properties associated with a given pyOM source and within a given fire zone could be considered as replicates (i.e., properties at U1 and U2 are replicates for unburnt treatment while B1, B2, B3, and B4 are replicates for burnt treatments) and can be spatially-aggregated based on pyOM feedstock and fire history without any significant loss in the amount of variability captured in the study. Temporal Trends in pyOM-Associated Carbon and Nitrogen Spatially-aggregated, temporal trends in relative total carbon and nitrogen content associated with pine-or cordgrass-derived pyOM in the unburnt and burnt zones are presented in Figures 2, 3. Over the study period, carbon content associated with pyOM varied between 0.75 and 1.25 times that of initial values in both fire zones and across both pyOM source. Interestingly, relative changes in carbon content associated with pine-derived pyOM followed that of cordgrass-derived pyOM throughout the study irrespective of fire history. Summer was characterized by a continued decrease (15-25% loss of initial) in pyOM-associated carbon across both pyOM sources and fire history. Fall through winter was characterized by an apparent re-accumulation of carbon, lost during summer, while spring was characterized by an additional 17-25% increase in pyOMassociated carbon, relative to the initial carbon content.
Like pyOM-associated carbon, temporal trends in pyOMassociated nitrogen for pine-vs. cordgrass-derived pyOM were comparable in summer-decreasing slightly to ∼0.96 times the initial (Figure 3). However, thereafter, the dynamics of pyOMassociated nitrogen for pine-and cord-grass derived pyOM followed very different trajectories within a given fire zone. For pine-derived pyOM, no further changes in associated nitrogen content were observed over fall, winter and spring in either the unburnt or burnt zones. In contrast, for cordgrass-derived pyOM, an increase of up to 1.3 times in initial pyOM-associated nitrogen was observed through Fall, peaking in winter at 1.2-1.3 times initial values and remaining stable thereafter until the end of the study in spring.
Precipitation (P cum ), soil residence time (SRT), average temperature (T avg ) and the interaction between P cum and SRT were found (through GLM analyses) to be the significant (p < 0.05) drivers of relative pyOM-associated carbon content. For observed changes in pyOM-associated nitrogen, significant drivers were P cum , SRT, the interactive effects between P cum and SRT, interactions between T avg and SRT and that between T avg and P cum . The final generalized linear model for pyOMassociated organic carbon was, TC t /TC 0 = 2.84 × 10 −4 (P cum ) − 1.64 × 10 −3 (SRT) − 5.67 × 10 −3 (T avg ) + 1.11 × 10 −6 (SRT × P cum ) + 1.123, was applicable across fire history (burnt vs. unburnt zone) and pyOM source (cordgrass vs. pine), and captured 77-91% of variability in measured values (Figure 2 and Supplementary Figure 3).

Temporal Changes in pyOM PH and Stability
Throughout summer and fall, the pH of the pine-derived pyOM remained relatively stable at 6.49 ± 0.24 in the burnt and unburnt zones, before decreasing through winter and early spring to 5.42 ± 0.18 in the burnt zone and 5.03 ± 0.10 in the unburnt zone (Figure 4). As spring progressed, pH of pinederived pyOM rebounded to 6.36 ± 0.14 by the end of the study on day 350. For cordgrass-derived pyOM, pH decreased with time from 7.89 ± 0.04 at the start of the experiment in summer to 5.94 ± 0.11 and 4.91 ± 0.07 in the burnt and unburnt zone, respectively, by early spring (day 294). As spring progressed, as with pine-derived pyOM, the pH of cordgrass-derived pyOM also rebounded-reaching 7.15 ± 0.10 in the burnt zone and 5.96 ± 0.03 in the unburnt zone by the end of the study on day 350.
Relative recalcitrance/stability, R 50 , , values of the pine-and cordgrass-derived pyOM are shown in Figure 5. These values showed an initial difference in the recalcitrance of pine-and cordgrass-derived pyOM, with R 50 values starting at 0.524 and FIGURE 2 | Temporal variation in relative pyOM-associated carbon for pine-and cordgrass-derived pyrogenic organic matter buried to a soil depth of 8 cm in either the (A) unburnt or (B) burnt fire management zone. Values are based on ash-free carbon content of pyOM at sampling, relative to that of the initial material. Lines indicate model trace from generalized linear modeling. 0.476, respectively. The relative recalcitrance, R 50 , of pinederived pyOM generally decreased with increasing soil residence time. The R 50 of pine-derived pyOM in burned and unburned areas stabilized after ∼200 days (R 50 = 0.512 and 0.518, respectively). Irrespective of fire history, the R 50 values of the cordgrass-derived pyOM increased in a logarithmic pattern with time but also converged around 0.508 and 0.518 -similar to the pine-derived pyOM. Not only were the trajectories in R 50 values of pine-and cordgrass-derived pyOM different, but the rates of change varied as well. For pine-derived pyOM the sigmoidal decrease in R 50 occurred around late summer-early fall and reaching a minimum in winter. For cordgrass-derived pyOM, the increasing segment of the R 50 -time relationship reached a maximum R 50 around late fall-early winter. Generalized linear modeling of changes in R 50 was possible for cordgrass-derived but could not be successfully applied to pinederived pyOM. In the case of cordgrass-derived pyOM, variability in R 50 in both the burnt and unburnt zone was shown to be predominantly driven by interaction-only effects (between FH and T avg , SRT and P cum , SRT and T avg , SRT and T avg ) and could be modeled as, R 50 = −3.57 × 10 −4 (FH × T avg ) − 1.06 × 10 −7 (SRT × P cum ) − 1.38 × 10 −5 (SRT × T avg ) + 7.16 × 10 −6 (P cum × T avg ) + 0.480, with the model accounting for 88-93% of variability in both cases (Supplementary Figure 6). However, such an assessment was not possible for pine-derived pyOM-because none of the tested co-factors in our main-effects only, interactions-only or mixed effects models were significant at the p = 0.05 level.

Impact on Soil Properties
Soils immediately below the buried pine-and cordgrass-derived pyOM did not show significant difference in pH across feedstock but were more acidic in the unburnt (pH = 5.74 ± 0.05) than burnt zone (pH = 6.00 ± 0.04; Supplementary Figure 7). By the first sampling excursion, the pH of the soil immediately beneath the litterbag in the unburnt zone increased to 5.92 ± 0.03 (in both pine-and cordgrass-derived pyOM) and in the burnt zone to 6.12 ± 0.03 (cordgrass-derived pyOM). This increase in pH was accompanied by a decrease in the electrical conductivity of the pyOM from 523 ± 2.29 to 21.7 ± 1.38 and 23.2 ± 3.40 µS cm −1 in the burnt and unburnt zone, respectivelywith 66-79% of this decrease occurring within the first 21 days (Supplementary Figure 8). The concomitant decline in soil pH and pyOM electrical conductivity continued over the fall and winter excursions. In the spring sample excursions, the pH of the soil immediately below the litterbags increased to it initial values.
As with soil pH, SOC content immediately beneath the buried pyOM litterbags were similar under pine-and cordgrass-derived pyOM; averaging 2.08 ± 0.07 wt% in the burnt zone and 1.47 ± 0.07 wt% in the unburnt zone (Supplementary Figure 9). In samples from excursions during the latter part of the fall into winter, there was evidence of a comparatively lower SOC content, compared to samples from other excursions.

Dynamics of pyOM-Associated Carbon and Nitrogen
Early stage decreases in pyOM-associated carbon and nitrogen, such as that observed during the summer season of the study, has been attributed to loss of labile components of pyOM via microbially-mediated degradation, volatilization and/or precipitation-induced leaching (Rajkovich et al., 2012;Norwood et al., 2013;Mukherjee et al., 2014;Sorrenti et al., 2016). At first glance, the microbially-mediated degradation of labile components in pyOM during summer seemed plausible and was supported to some extent by the fact that most of the loss in carbon and nitrogen (between sampling days 42 and 84) coincided with a 4-5 • C increase in air temperature from 31 to 36 • C (Supplementary Figure 10). These temperatures, as well as their observed increase during periods of highest carbon and nitrogen loss, were well within the range for optimal soil microbial activity and (all else being constant) would be expected to trigger increased microbial activity (Andersson and Nilsson, 2001;Liang et al., 2003).
Volatilization and/or precipitation-induced leaching/flushing of water-soluble organic components for the observed decrease in pyOM-associated carbon and nitrogen were also plausible. With the observed increase in T avg over the summer, the volatilization of volatile organic compounds (VOCs) formed during pyOM formation can be expected; and has been previously linked to the loss of CO 2 -C from pyOM with subsequent oxidation of the pyOM surface (Bruun et al., 2008;Steiner et al., 2008;Rajkovich et al., 2012). Potential support for precipitationinduced leaching/flushing was provided by the occurrence of several rainfall events (ranging in size from 2.54 to 43 mm) during this period (Supplementary Figure 10). As noted earlier, these events were accompanied by sharp declines in the electrical conductivity of the buried pyOM, indicative of the flushing of salts from the materials (Supplementary Figure 8). It is therefore plausible that rain-induced flushing of salts from pyOM would also trigger the concomitant leaching of water-soluble carbon and nitrogen from the pyOM as well.
Direct evidence from δ 13 C isotopic signatures of the residual pyOM was more supportive of precipitation-induced leaching of water-soluble organic components (than the microbial-mediated degradation) as significant driver for observed organic carbon loss in this study. The average δ 13 C in pine-and cordgrassderived pyOM on day 21 was −28.4 ± 0.05 and −23.9 ± 0.22 respectively and did not change significantly in the unburnt or burnt zones between days 21 and 84-when pyOM-associated carbon loss was observed (Figure 6). This lack of evidence for a significant shift in the δ 13 C signature of pyOM-associated carbon was inconsistent with expectations for the microbialmediated oxidation in aerobic surface soils such as those in the study. In such soils, microbial degradation is accompanied by 12 C enrichment of released CO 2 and 13 C enrichment of the residual OM (Ågren et al., 1996;Feng, 2002;Wang et al., 2015). Conversely, the leaching of water-soluble organic components from pyOM without significant microbial transformation can be expected to have little or no effect on δ 13 C signature (Cleveland et al., 2004). The role of carbon mineralization via the volatilization of VOCs could not be ruled out as significant contributor to the observed loss of pyOM-associated carbon. In fact, given (1) the lack of isotopic support for microbial mineralization and (2) that P cum , SRT and T avg (i.e., the GLM model in Figure 2) accounted for only 6 of 15-25% of pyOMassociated carbon lost in the first 84 days (∼3 months), it was quite justifiable that of the three potential processes; volatilization was the process that would account for the greatest loss of pyOM-associated carbon observed in this study.
The scenario would be somewhat different for losses of pyOMassociated nitrogen. As the model fit in Figure 3 indicatesunlike for pyOM-associated carbon -losses in pyOM-associated nitrogen during summer could be fully accounted for by leaching (i.e., variations in P cum , SRT, and T avg ). The fact that losses in pyOM-associated nitrogen were accompanied by an approximate 0.5 enrichment in the 15 N signature of the residual pyOM (Figure 7) suggested that the water-soluble nitrogen lost from pyOM during leaching would be more depleted in 15 N than the residual pyOM fraction remaining. 15 N-enrichment in soils has also been attributed to adsorption of microbially processed organic matter and/or increased abundance of soil fungi (Natelhoffer and Fry, 1988;Etcheverría et al., 2009;Wallander et al., 2009). However, neither of these mechanisms was well-supported by our data. Firstly, the net loss of nitrogen, observed during summer, would rule out adsorption as important contributor to 15 N enrichment. Secondly, there was no evidence for significant fungal colonization of pyOM during the early stages of experiment associated with summer losses of pyOMassociated nitrogen. In fact, temporal trends in ergosterol content of pyOM pointed to summer as being the least biotically-active period-from a fungal perspective-in the study (Supplementary Figure 6). For instance, it is well-known that fungi are the primary degraders of pyOM in forest soils (Steinbeiss et al., 2009), and as such the ergosterol in fungal cells is widely used as a biomarker for fungal activity (Djajakirana et al., 1996;Mille-Lindblom et al., 2004). However, the fact that measured ergosterol content of pyOM was significantly lower in summer (2-5 µg g −1 of pyOM) than in spring (13 µg g −1 of pyOM) was consistent with the lowest fungal/biotic activity coinciding with greatest and only occurrence of losses in pyOM-associated carbon (Supplementary Figure 11).
The re-accumulation of pyOM-associated carbon and nitrogen we observed throughout fall and into spring was consistent with both the trajectory and magnitudes reported by  Mukherjee et al. (2014). In their study, Mukherjee et al. (2014) found that the burial of oak, pine and grass pyOM (produced at ≥400 • C) for a soil residence time of 15 months (May 1, 2009 to September 1, 2010) resulted in an increase of up to 124 and 143% in associated carbon and nitrogen, respectively. They attributed the observed increases to variations in the two soils and pyOM type studied. Interestingly, the experiments of Mukherjee et al. (2014) were also conducted within the Long leaf pine region about 490 miles east of our site, under comparable average temperatures and cumulative rainfall but in very different soils. That our study conditions differed from Mukherjee et al. (2014) primarily in soil types but yielded similar overall trends in pyOM dynamics brings into question the overall importance of soil type in pyOM carbon and nitrogen dynamics in the  forested ecosystems of this area and most likely the region. Particularly, since like the spodosol and entisol of Mukherjee et al. (2014) and the ultisol at our site, the surface soils of the region are predominantly coarse (sandy) textured, with low pH, low organic matter and low retention capacities (Shaw et al., 2004;Novak et al., 2009). Other studies have also reported either the accumulation of carbon (Kasin and Ohlson, 2013;Dong et al., 2017) or nitrogen (Sorrenti et al., 2016;de la Rosa et al., 2018), but not both, after the field aging of pyOM for up to 5 years. It is also interesting to note that for studies having similar trajectories in pyOM carbon and nitrogen dynamics, there were suggestions of similar climate-either similar Koppen-Geiger climate classification (continental climate; Kasin and Ohlson, 2013;Dong et al., 2017) or similar annual precipitations and temperatures (Sorrenti et al., 2016;de la Rosa et al., 2018). In all these cases, carbon or nitrogen accumulations in pyOM were attributed to the sorption of exogenous components (inorganic or organic) and/or microbial/fungal colonization.
Our data showing accumulations of pyOM-associated carbon and nitrogen, from fall through spring, was consistent with this latter explanation of exogenous carbon and nitrogen sorption to and/or the microbial colonization of the pyOM. Specifically, the climate-driven leaching of microbially "reworked" exogenous water-soluble organic matter and its subsequent sorption to the pyOM. Firstly, the fact that both carbon and nitrogen accumulations could be effectively modeled by considering only variations in P cum , T avg , and SRT points to the critical role of climatic conditions (Figures 2, 3). Secondly, the exceedance of the initial carbon and nitrogen values as shown in Figures 2, 3 for this period coupled with an accompanying 13 C depletion in carbon (from δ 13 C = −23.9 to −25.7%o; Figure 6) associated with cordgrass-derived pyOM was reflective of an exogenous source that is isotopically lighter than cordgrass. Thirdly, the presence of 15 N enrichment in cordgrass pyOM was reflective of the associated exogenous N being isotopically heavier than the residual pyOM-consistent with microbial processing of this nitrogen prior to its interaction with pyOM (Figure 7). The fact that nitrogen accumulation was only apparent in cordgrass-derived pyOM suggested the preferential association of this exogenous nitrogen with functional groups that are more prevalent in the cordgrass-than the pine-derived pyOM.
Increasing quantities of base-extracted organic carbon (BEOC) associated with both pine-and cordgrass-derived pyOM, over the course of the study, further supported an exogenous origin for the re-accumulated pyOM-associated carbon (Supplementary Figure 12). The increasing trends in BEOC with SRT also suggested that the plateau in carbon re-accumulation on pyOM around late fall or winter was more reflective of a lowered availability of "labile" carbon rather than saturation of the pyOM sorption capacity. A lowered availability of organic carbon in late fall and winter was attributable to lowered temperatures and rainfall favoring a decrease in microbial processing of organic matter, leaching and subsequently a decrease in water-soluble organic matter available for sorption. The lack of δ 13 C depletion (or enrichment) for the pine-derived pyOM (δ 13 C = −28.3‰; Figure 6) pointed to leachate from pine needles, fungal biomass and other similarly 13 C-depleted components of litter/soil organic matter as possible sources (McDowell and Likens, 1988;Qualls and Haines, 1991;Wallander et al., 2009) Increasing availability of this more 13 C-depleted carbon-due to needle shedding, senescence and biotic degradation-would account for the re-accumulation of pyOM-associated carbon in early to mid-fall (Figure 2). Lower temperatures and rainfall, as well as lower biotic activity would explain the reduce availability and hence plateauing in carbon re-accumulation during the latter parts of fall into winter. The return of higher temperatures and rainfall would explain further increases in pyOM-associated carbon content (from 1 to 1.25 times initial values) observed during spring. As temperature increases in spring, an increase in biotic activity, resulting in more water-soluble organic matter being available for association with pyOM via for sorption-as in fall. However, as opposed to fall when microbial processing of senescenced plant biomass would drive availability of exogenous carbon, in spring root exudates released to enhance nutrient uptake (as plant shift from dormancy to active growth) would be a potential source of exogenous carbon available for sorption (Mergel et al., 1998;Nardi et al., 2000;Richardson et al., 2011;Keiluweit et al., 2015).
It is worth noting that although increased pyOM-associated ergosterol was observed with time in all but our pine-derived pyOM in the unburnt zone (Supplementary Figure 11) the role of fungal colonization in observed pyOM dynamics is not unequivocal. For instance, the range of ergosterol values were well within expectations for mineral soils without fungal colonization (Sung et al., 1995;Djajakirana et al., 1996;Frostegård and Bååth, 1996). Therefore, it is quite plausible that the observed increase ergosterol with SRT was exogenous and a result of its sorption to the pyOM rather than colonization by fungi. Such persistence is not unexpected given the likely protection of fungal biomass by forest litter from sunlight (Mille-Lindblom et al., 2004).

Dynamics of pyOM pH and Stability
The differences observed in initial pH in the study was attributable to differences in source material (i.e., pine-vs. cordgrass-derived pyOM; de la Rosa et al., 2018). The generally lower pH values and lower pH-sensitivity of pine-derived pyOM, compared to its cordgrass-derived counterpart, was reflective of the acidic conditions typical of the soils in pine forest ecosystems. Cordgrass, on the other hand, is native to coastal regions where pH is usually more basic and more reflective of seawater pH (around pH 8-8.2). The pH of a basic material in an acidic environment is expected to decline as it adjusts to equilibrium conditions, thus the observed decline in pH in cordgrass-derived pyOM. Declines in pyOM pH with aging/SRT has been widely observed and is attributable to the weathering-induced oxidation of basic to acidic functionalities (e.g., alcohols to carboxyls), the dissolution of carbonates or (oxy)hydroxides and/or the leaching and exchange of basic cations with H + or acidic cations (Joseph et al., 2010;Yao et al., 2010;Spokas et al., 2012;Sorrenti et al., 2016;de la Rosa et al., 2018). All three of these mechanisms were quite plausible in our study given the observed flushing of salts and observations by others indicating increased ion exchange with aging of pyOM (Cheng et al., 2008;Mukherjee et al., 2014). The fact that the largest declines in pyOM pH (1-2 pH units; fall through early spring) coincided with period of accumulation for exogenous carbon and nitrogen further points to the acidic nature of the exogenous organic matter being leached from litter and partitioned to the pyOM.
The fact that the decline in pH at day 217 was seen in both types of pyOM and is associated with the transition from winter to spring, suggested that the change was likely linked to a seasonally-driven shift in environmental conditions. The most likely trigger is the occurrence of several large precipitation events throughout winter into early spring (between days 210 and 240; Supplementary Figure 10). This period represented the most intense period of precipitation, resulting in an influx of more acidic exogenous materials and the observed decline in pH of both pine-and cordgrass-derived pyOM between sampling day 217 and 294. This would also be consistent with results from GLM indicative of P cum , T avg , and SRT being significant parameters in explaining observed variations in the pH of pyOM. The rebound of pyOM pH back to more circumneutral conditions (pH 6.0-7.2) by the end of the study was congruent with eventual loss of acidic exogenous materials via their neutralization through the ingress of basic cations mobilized by root exudates into the pyOM. The difference in the absolute magnitudes of pH declines and subsequent rebound observed across fire history can be reasonably explained as a fire-induced soil buffering effect; where previous burning imparts a higher soil buffering capacity to soils in the burnt zone (Mukherjee et al., 2014;Wang et al., 2014;Shi et al., 2017Shi et al., , 2018. The source of this buffering would be increased cation exchange capacity and basic cation content arising from higher soil carbon and ash-contained basic cations contents accumulating over progressive burns. Therefore, these soils are better able to resist decreases in pH arising from ingress of exogenous organic matter and/or H + . This fire-induced buffering effect would also explain the significance of FH (specifically its interaction with T avg ), as identified via GLM, as an important parameter in explaining observed variability in pH. We posit that as T avg decreases and plants shift into dormancy, the interaction of precipitation with senesced/senescing plant matter results in the leaching of exogenous organic matter (and H + ) that flushes basic cations and gets sorbed to pyOM reducing its pH. A rebounding of T avg in spring, the return of active plant growth and the associated release of root exudates triggered the release of basic cations from minerals with the leaching and ingress of these basic cations into pyOM triggering its increase in pH observed in the latter part of the study.
Measured R 50 values provided a proxy for the bulk recalcitrance (relative to graphite) of the pyOM remaining in the soil (Harvey et al., 2012). Since the exogenous carbon and nitrogen interacting with the pyOM was soil dissolved organic matter leached from the same source; the observed changes in R 50 were reflective of changes in the innate recalcitrance of the original pyOM material and/or the induced effects of the exogenous OM. Of such, observed R 50 provided insights into the overall recalcitrance or so called "quality" of the remaining pyOM.
The different temporal trajectories in R 50 observed in pinevs. cordgrass-derived pyOM can be reasonably explained by considering differences in feedstock structure and effect of associated water-soluble and exogenous OM. The interactions of these factors and subsequent impacts on recalcitrance (R 50 ) are summarized in Figure 8. Initial differences observed in R 50 were attributable to structural differences between the pine-and cordgrass-derived pyOM (mainly the higher fraction of lignin in wood vs. grass feedstocks) that has been wellrecognized in a number of previous studies (e.g., McBeath and Smernik, 2009;Harvey et al., 2012). Over the summer, the lack of a significant change in R 50 values for pinederived pyOM indicated that the associated flushing of watersoluble OM (and salts) from this pyOM had no significant effect on its relative recalcitrance. This was consistent with the leached water-soluble components of pine-derived pyOM having comparable recalcitrance as the remaining residual pyOM component. By comparison, the increase in R 50 values for the cordgrass-derived pyOM-over the same period-was consistent with the flushed water-soluble OM components of this pyOM being of lower recalcitrance than the remaining residual pyOM.
The decline in R 50 of pine-derived pyOM with the onset and subsequent increase in the accumulation of exogenous OM (in the pyOM), over fall and winter, was reflective of this exogenous OM being of lower recalcitrance than the residual pyOM component. Therefore as the amount of exogenous OM associated with the pine-derived pyOM increased, its negative effect on R 50 was enhanced until an apparent plateauing point was reached in the latter part of winter. In contrast, the corresponding increase in R 50 for cordgrass-derived pyOM observed over fall and winter was congruent with exogenous OM having a positive effect (i.e., R 50 of exogenous OM is greater than that for residual pyOM) on the overall recalcitrance of the cordgrass-derived pyOM; with this positive effect on recalcitrance also reaching a plateauing point in late winter.
Similarities in the timing and R 50 values at which plateauing occurred in cordgrass-and pine-derived pyOM, within a given fire history, pointed to a common "equilibrium" R 50 (R 50, eq ). In addition to being reflective of the overall effect of residual pyOM and exogenous OM on pyOM recalcitrance, this R 50, eq provide some interesting insights into overall pyOM recalcitrance as it pertains to, 1) fire history at the study site and 2) potential priming effects of exogenous OM on pyOM recalcitrance. Lower values for R 50, eq in the burnt vs. unburnt zone suggested that pyOM was likely to have a comparatively lower persistence in soils with longer fire histories (Figure 8). Lower persistence of pyOM in soils with a history of fires have been attributed to the presence of microbial populations that are well-adapted to utilizing fire-derived carbon (Khodadad et al., 2011;Masiello et al., 2013). However, our results suggested that leaching and the modifying effects of exogenous OM are also important factors in pyOM persistence. For instance, in the case of cordgrass-derived pyOM the difference in R 50, eq between the burnt and unburnt zones was established during the pyOM flushing stage-indicative of leaching differences. By comparison, in pine-derived pyOM, the precedence for the observed difference in R 50, eq between the burnt and unburnt zones was established during the stage of accumulation of exogenous OM-indicative of differences in the interaction between pyOM and exogenous OM in the fire zones being an important factor. In terms of potential priming effects, the higher residual R 50 (R 50, res ) at the end of summer for pine-derived materials compared to the corresponding R 50, eq (i.e., R 50, res > R 50, eq ) is consistent with exogenous OM reducing the overall recalcitrance of pyOM and hence favor positive priming (Figure 8). In contrast, the reverse situation (R 50, res < R 50, eq ) existed for the cordgrass-derived pyOM, indicative of exogenous OM in the study area increasing overall recalcitrance of pyOM and hence is most likely to result in negative priming of pyOM. The higher absolute magnitude of R 50, res minus R 50, eq in burnt vs. unburnt fire zones suggest that the priming effects of exogenous OM on pyOM degradability would be more pronounced in soils with an history of fire/pyOM input.

Temporal Changes in Soil Within the Vicinity of the pyOM
The overall acidic pH conditions of the soils were consistent with expectations, for soils within a pine forest ecosystem (Falkengren-Grerup, 1987) while the higher pH conditions observed in burnt vs. unburnt soils were consistent with the presence of basic salts (e.g., CaO and CaCO 3 )-produced as components in the ash byproduct of plant biomass pyrolysis (Bodí et al., 2014;Wang et al., 2014 and references therein). The concomitant increase in the pH of the soil immediately beneath the litter bags, with loss in pyOM electrical conductivity, after the first excursion was consistent with in-situ flushing of basic salts from ash (13-14 wt%) within the pyOM into the immediate soil. The role of ash components in early soil pH modification was well-supported by evidence from temporal changes in the electrical conductivity of CG-derived pyOM consistent with enhanced in-situ flushing of the pyOM within the first 21 days of the study (Supplementary Figure 8). For example, over the course of the study the electrical conductivity of the CG-derived pyOM decreased from 523 ± 2.29 µS cm −1 to 21.7 ± 1.38 and 23.2 ± 3.40 µS cm −1 in the burnt and unburnt zone, respectively. However, 66-79% of this decrease in electrical conductivity occurred within the first 21 days-congruent with the vast majority of basic salts being leached from the pyOM during this period into the soil below.
The continued decline in soil pH and pyOM electrical conductivity observed for the fall and winter excursions reflected continued leaching and the loss of basic salts from the soils immediately below the litterbags. The subsequent increase in soil pH for the spring sample excursions suggested either a latter influx of basic salts into, or loss of acidic components from, the soil immediately below the pyOM. However, the lack of an accompanying decrease in pyOM electrical conductivity evidence would exclude ash from the pyOM as the source of the observed increase in soil pH in the spring. The actual source of this increase in soil pH is unclear at this time but is likely linked to temperature-and precipitation-enhanced process(es)in particular the root exudate enhanced mobilization of base cations from minerals as plants return to active growth from dormancy.
Higher SOC content in the burnt zone was consistent with the migration and storage of pyOM from previous burning events below the soil surface (Supplementary Figure 9). The lower SOC content in samples from late fall and winter, compared to the other seasons, was congruent with decreasing soil pH trends and suggest an association of leaching with lower SOC immediately below the buried litterbag (Supplementary Figure 4).

SUMMARY AND CONCLUDING REMARKS
In an effort to better understand the factors driving earlystage pyOM dynamics we studied the chemical evolution of freshly produced cordgrass-and pine-derived pyOM over the first year post-burial in pine-forest soils with similar origins (parent material and texture) but very different fire management histories. Results from the study indicated that spatially-driven factors had little to no effect on pyOM dynamics in the shortterm (up to 1 year). By comparison, temporally-variable factors such SRT, P cum , T avg and to a lesser degree the land management factor, FH were found to be the primary contributors to observed variation in the biogeochemical dynamics-specifically pyOM-associated carbon content, nitrogen content, pH and bulk recalcitrance-of pyOM from different sources. Interestingly, we found no evidence to support microbially-driven mineralization of pyOM as an important driver of pyOM dynamics in this study. Instead, the volatilization of VOCs and/or the leaching of water-soluble components from the pyOM were most plausibly controlling pyOM dynamics in the first 3 months (i.e., summer) of the study; while in the next 9 months (i.e., fall through spring), the dynamics of pyOM was likely controlled by availability of exogenous organic matter and its interaction with pyOM.
The volatilization of VOCs and the flushing of water-soluble components from pyOM (including water-soluble organic matter) during the first 3 months resulted in (1) a decrease in pyOM-associated carbon, (2) minimal effect on pyOM-associated nitrogen or the overall stability of pine-derived pyOM but (3) increased the overall stability of cordgrass-derived pyOM. The increasing availability of labile exogenous organic matter (from needle fall and senescence in fall and root exudates in spring) resulted in a re-accumulation/accumulation of pyOMassociated carbon and nitrogen (for cordgrass-derived pyOM) in fall through to spring. Interaction of this exogenous organic matter with pyOM caused an initial decrease and increase in the recalcitrance of pine-derived and cordgrass-derived pyOM, respectively -with an eventual convergence in recalcitrance for both pyOMs by mid-winter. A higher recalcitrance convergence point in the unburnt zone was congruent with a longer potential persistence of pyOM in these soils compared to that in the burnt zone. The decreasing effect of exogenous organic matter on pyOM recalcitrance suggests that native non-pyrogenic organic matter in the study area is likely to have a positive priming effect on any eventual microbial degradation.
The overall lack of evidence to support a significant microbiological influence-contrasted against that for volatilization, leaching of water-soluble pyOM, and sorption of exogenous organic matter-pointed overwhelmingly to physicochemical (rather than microbiological) processes as the primary drivers of early pyOM dynamics in our study system. Such results suggest that models of pyOM behavior (especially in the early stages) in humid subtropical systems may need to adopt a multifactorial approach in determining the factors dictating pyOM dynamics.