Impacts of litter decay on organic leachate composition and reactivity

Litter decomposition produces labile and recalcitrant forms of dissolved organic matter (DOM) that significantly affect soil carbon (C) sequestration. Chemical analysis of this DOM can provide important knowledge for understanding soil DOM dynamics, but detailed molecular analyses on litter derived DOM are scarce. Here we use ultrahigh resolution mass spectrometry (FT-ICR MS) to characterize the molecular composition of DOM from fresh and progressively decomposed litter samples. We compared high reactive (HR) and low reactive (LR) litter sources with regard to changes in the chemistry and bioavailability of leachates throughout the early phase of litter decay. We show that litter reactivity is a driver of chemical changes in the leached DOM of litter species. Birch, alder and Vaccinium (i.e. HR) litter initially produced more DOM with a higher lability than that of spruce, pine and wood (i.e. LR) litter. Labile oxidized phenolic compounds were abundant in leachates produced during the initial HR litter decay stages, indicating litter lignin degradation. However, the similarity in chemistry between HR and LR leachates increased during the litter decay process as highly leachable structures in HR litter were depleted. In contrast, chemistry of leachates from LR litter changed little during the litter decay process. The oxygenated phenolic compounds from HR litter were driving the lability of HR leachates and the changes in relative abundance of molecules during DOM incubation. This appeared to result in the creation of stable aliphatic secondary microbial compounds. In LR leachates, lability was driven by labile aliphatic compounds, while more resistant phenolic compounds were associated with recalcitrance. These results show how DOM dynamics follow different paths depending on litter reactivity, which has important implications for soil biogeochemistry and C sequestration.

decomposed litter samples. We compared high reactive (HR) and low reactive (LR) litter sources with regard to changes in the chemistry and bioavailability of leachates throughout the early phase of litter decay. We show that litter reactivity is a driver of chemical changes in the leached DOM of litter species. Birch, alder and Vaccinium (i.e. HR) litter initially produced more DOM with a higher lability than that of spruce, pine and wood (i.e. LR) litter. Labile oxidized phenolic compounds were abundant in leachates produced during the initial HR litter decay stages, indicating litter lignin degradation. However, the similarity in chemistry between HR and LR leachates increased during the litter decay process as highly leachable structures in HR litter were depleted. In contrast, chemistry of leachates from LR litter changed little during the litter decay process. The oxygenated phenolic compounds from HR litter were driving the lability of HR leachates and the changes in relative abundance of molecules during DOM incubation. This appeared to result in the creation of stable aliphatic secondary microbial compounds. In LR leachates, lability was driven by labile aliphatic compounds, while more resistant phenolic compounds were associated with recalcitrance. These results show how DOM dynamics follow different paths depending on litter reactivity, which has important implications for soil biogeochemistry and C sequestration.

Introduction
Boreal forests cover nearly 15% of the earth's land surface (Melillo et al. 1993;Gower et al. 2001) and contain as much as 370-1700 Pg of carbon (C), stored to * 95% below ground (Bradshaw and Warkentin 2015). Of the total net primary production (NPP) in boreal forests, 20-60% contributes to litterfall (Chen et al. 2017) and, thus, to subsequent decomposition on the forest floor, where 7-40% of litter C loss leaches as dissolved organic matter (DOM) . Through infiltrating rainfall and snowmelt water, this DOM is conveyed through soils where it plays a key role in controlling soil formation and C sequestration (Kalbitz et al. 2000Cotrufo et al. 2015). Thus, DOM is an important link between NPP and C storage in boreal regions.
Both labile and recalcitrant fractions of DOM can lead to the formation of persistent soil organic matter (SOM) (Grandy and Neff 2008;Prescott 2010;Cotrufo et al. 2015). Labile DOM is an important microbial C source contributing to the production of secondary microbial products (Kalbitz et al. 2000;Bowen et al. 2009), which have been shown to form the most stable and persistent fraction of SOM (Marschner and Kalbitz 2003;Miltner et al. 2012). The portion of DOM that is more recalcitrant, on the other hand, consists partly of aromatic molecules with high affinity for stabilization through sorption Kleber et al. 2007). Although some labile compounds such as amino acids can also be prone to sorption (Kleber et al. 2007), the decrease in decay is smaller for sorbed labile DOM than recalcitrant DOM . Thus, while the quantity of the litter to soil DOM flux is imperative to the global C cycle, the chemical composition of the DOM is key to understanding the mechanisms controlling the fate and preservation of soil DOM.
Litter chemistry has been proposed to be a key regulator of both initial rates of the mass loss and, specifically, of the loss of C through DOM leaching (Berg and Staaf 1980;Melillo et al. 1982Melillo et al. , 1984Aber et al. 1990;Cornelissen 1996;Pérez-Harguindeguy et al. 2000;Preston and Trofymow 2000;Soong et al. 2015;Campbell et al. 2016). Litter mass loss and leaching tend to go hand in hand, as early phase decay is faster in litter with a higher soluble content (Osono and Takeda 2005). Over time, soluble compounds like polyphenols, certain carbohydrates and non-lignified cellulose are exhausted, leading to decreasing rates of mass loss (Berg et al. 1982;Melillo et al. 1989;Lorenz et al. 2000Lorenz et al. , 2004. In addition, litter containing highly labile DOM compounds is subject to higher lignin degradation rates as a result of priming, with phenolic compounds released as byproducts (Klotzbücher et al. 2011b). Thus litter with rapid initial mass loss undergo greater changes in chemical composition and leach more DOM compared with litter that is initially more resistant to decay . There are good reasons to expect a corresponding difference in the composition of leachates between rapid and slow decaying litter species, but little is known about the chemical composition of DOM from different litter species and the changes throughout litter decay.
The characterization of litter-leached DOM has often not been detailed enough to provide a molecular understanding of the leaching dynamics and involved processes (Brock et al. 2020). For example, while it is well-documented that there are large variations in microbial bioavailability of DOM (BDOM) between DOM leached from different litter sources (Cleveland et al. 2004;Silveira 2011;Hensgens et al. 2020) and at different litter decay stages (Don and Kalbitz 2005;Hensgens et al. 2020), it is unknown if corresponding differences exist on a molecular composition level. This leaves the fundamental question of how BDOM relates to litter molecular DOM composition unanswered. At the same time, while soil science is largely moving away from chemical humification preservation models (Schmidt et al. 2011;Lehmann and Kleber 2015), knowledge about the chemical composition of litter derived DOM and the changes occurring during litter decay may prove instrumental in the development and testing of alternative soil C preservation models (Kaiser and Kalbitz 2012).
A variety of different techniques exist to measure the chemistry of DOM and are reviewed in (Minor et al. 2014). Fourier-transform ion cyclotron resonance mass spectrometry (FT-ICR MS) provides ultrahigh accuracy mass measurements and is widely used for DOM characterization (Nebbioso and Piccolo 2013). However, it requires ionization of the DOM molecules, and when combined with direct infusion measurements is prone to ionization artifacts like ion suppression and matrix effects (Nebbioso and Piccolo 2013). In addition, solid phase extraction (SPE), which is required for direct infusion FT-ICR MS does not recover all DOM and is known to selective enrichment of non-polar compounds from freshwater DOM (Raeke et al. 2016). Lastly, quantification of DOM compounds is not possible because there are currently no standards available for DOM. While other techniques, such as fluorescence spectroscopy, Fouriertransform infrared (FTIR) or nuclear magnetic resonance (NMR), can be used for quantification, they lack the (molecular) resolution to advance the field on the molecular scale. In this study, we used ultra-high resolution FT-ICR MS measurements on litter DOM leachates from eight different sources in a field litter decomposition experiment, with the goal to explore the dynamics in DOM chemistry as related to litter decomposition and to BDOM. The DOM was extracted from fresh, 10-, 48-and 180-day decomposed litter (Fig. 1). FT-ICR MS measurements were done before and after half-year DOM incubations, during which DOM concentrations were measured periodically and its reactivity and BDOM were calculated. We tested the idea that similar to how litter reactivity is known to govern the initial changes in litter chemistry, it also controls the corresponding initial changes in leachate chemical composition, with consequences for DOM bioavailability. In order to explore the effects of litter reactivity on DOM chemistry we conceptually divided the litter species into two separate groups based on litter mass loss and the solubility of litter C: 'high reactive' (HR) and 'low reactive' (LR) litter (Fig. 1). Early litter mass loss has been shown to be comparatively low in spruce and pine litter as compared to alder or birch litter (Berg and Ekbohm 1991). In addition, as early litter decomposition is inversely correlated with structural strength (Pérez-Harguindeguy et al. 2000), we assigned birch, alder and V. myrtillus, as HR litter and the more tough structured spruce, pine and wood litter as LR litter.
The objective of this project was to compare HR (i.e. birch, alder and Vaccinium) and LR (i.e. spruce, pine and wood) litter species with regard to changes in DOM chemistry and bioavailability of leachates throughout the early phase of litter decay. More specific, we hypothesized that as a result of highly soluble material in HR litter at the start of litter decay (1.1) the DOM reactivity is higher in HR than LR litter species and (1.2) it decreased relatively more rapid for HR litter species during litter decay, as highly soluble and reactive material in parent litter is exhausted over time. Furthermore, we predict that (2.1) as litter Fig. 1 Conceptualizing high reactive (HR) and low reactive (LR) litter species based on litter mass loss rates and C solubility. Simplified data from (Hensgens et al. 2020). Dashed lines show mass loss and DOM leaching potential (mg C g -1 ) at 10, 48 and 180 days of litter decay for HR (yellow) and LR (blue) species. In reality, LR species could show either a slight decrease or a slight increase in leachable material depending on its source litter (Don and Kalbitz 2005;Hensgens et al. 2020) decomposition progresses, leachate chemistry will show greater changes in HR as opposed to LR species, while (2.2) a higher degree of similarity is seen in leachate chemistry towards the later stages of litter decomposition. Finally, if the above hypotheses can be supported, we expect that fresh HR litter leachates contain unique chemical features that are linked to their higher overall reactivity. Based on our findings we discuss how our approach and results provide useful knowledge for a wide range of subjects from soil C respiration, SOM buildup and soil DOM leaching into the aquatic network.

Sample collection
Fresh litter samples were collected on September 25 and 28, 2015, in the Krycklan catchment in northern Sweden (Laudon et al. 2013) from the most abundant litter species: Picea abies (spruce), Pinus sylvestris (pine), Betula spp. (birch), Alnus spp. and Vaccinium myrtillus. For spruce, pine and birch species small branches were collected as wood samples. Litter samples were dried, weighed and put into the field for in situ litter decomposition. Litter samples were retrieved after 0-, 10-, 48-and 180-day field decomposition.

DOM extraction and sampling
Following collection, pure water extractions were carried out using a 48-h shaking (140 RPM) protocol at room temperature (Hensgens et al. 2020). After filtration (0.7 lm; Whatman GF/F) and dilution to 5.9 mg C l -1 , samples were taken for FT-ICR MS analyses, acidified to pH * 2 using ultrapure HCl and immediately frozen (-24°C). In order to keep factors constant besides inherent C reactivity, and in line with the finding that soil heterotrophic communities are primarily carbon limited (Soong et al. 2020), the samples for microbial incubations were diluted to a final concentration of 6 mg C l -1 , enriched with inorganic nutrients to achieve a C:N:P ratio above the Redfield ratio (25:5:1) and added the same inoculum (1% v/v of unfiltered Krycklan stream water). Incubations were done in a dark constant room at 20°C in 400 ml HDPE bottles. After the end of the incubation (180 days), samples for FT-ICR MS analyses were taken, acidified and stored. All equipment had been acid washed and samples were kept in acid washed HDPE bottles. A more detailed description of the setup can be found in a companion paper (Hensgens et al. 2020), where DOM leaching is discussed in relation to litter production of the understory versus canopy species.

DOM decay dynamics
Dark incubations were initiated within two days from DOM extraction and DOM carbon concentrations measured after 0, 28 and 182 days of incubation. For fresh litter additional measurements were made after 14, 45 and 91 days. Parallel to this, 7-day dark incubations were started in 5 ml duplicates glass vials for each DOM solution using SDR (SensorDish Readers, PreSense) nonintrusive oxygen sensors for respiration monitoring every 2 h in a high temperature precision 20°C climate chamber (Soares et al. 2018). The DOM was measured as total organic carbon (TOC) on 30-40 ml samples acidified using 50 ll of ultrapure (20%) HCl and stored at 5°C before analysis performed on an OI analytical Aurora TOC analyzer.
Changes in dissolved oxygen concentrations (O 2 ) were transformed to DOM t values according to Eq. 1, where C mass is the molar mass of C and O2 mass is the molar mass of dioxygen.
This assumes a respiratory quotient (RQ) of 1, close to RQs measured in soils (Müller et al. 2004;Dilly and Zyakun 2008). Short-term data was trimmed to 1 value per day and merged with long-term incubation measurements.
The reactivity continuum model (Eq. 2 was fitted on DOM concentrations measured throughout the incubation converted to relative DOM concentrations DOM t /DOM 0 (Koehler 2012). Fitting was done using a general linear mixed effects (glme) model, with litter decay state as fixed and species as random factor.
FT-ICR MS methods

Sample preparation
Leachate samples (30-360 ml) were extracted via solid-phase extraction using an automated sample preparation system (FreeStyle, LC Tech, Obertaufkirchen, Germany) on 50 mg styrene-divinylpolymer type sorbens (Bond Elut PPL, Agilent Technologies, Santa Clara, CA, United States) to desalt the sample for subsequent DI-ESI-MS according to Raeke et al. (2016). The SPE-DOM was eluted with 1 ml methanol (Biosolve, Valkenswaard, Netherlands), diluted to 20 ppm and mixed 1:1 (v/v) with ultrapure water immediately prior FT-ICR MS analysis. Carbon based extraction efficiency was approximately (38.9 ± 10.4)% (n = 62). SRFA measured in triplicate was used to check instrument variability and solvent and extraction blanks were prepared.

FT-ICR MS measurement
An FT-ICR mass spectrometer equipped with a dynamically harmonized analyzer cell (solariX XR, Bruker Daltonics Inc., Billerica, MA, USA) and a 12 T refrigerated actively shielded superconducting magnet (Bruker Biospin, Wissembourg, France) instrument was used in ESI negative ionization mode (capillary voltage: 5 kV). Extracts were analyzed in random order with an autosampler (infusion rate: 10 ll min -1 ). For each spectrum, 256 scans were co-added in the mass range 150-3000 m/z with 15 ms ion accumulation time and 4 MW time domain (resolu-tion@400 m/z ca. 470,000). Mass spectra were internally re-calibrated with a list of peaks (247-609.2 m/z, n [ 52) commonly present in terrestrial DOM and the mass accuracy after linear calibration was better than 0.11 ppm (n = 62). Peaks were considered if the signal-to-noise (S/N) ratio was greater than four. Raw spectra were processed with Compass DataAnalysis 4.1 (Bruker Daltonics Inc., Billerica, MA, USA).

FT-ICR MS data evaluation
Molecular formulas were assigned to peaks in the range 0-750 m/z allowing for elemental compositions C 1-60 H 0-122 N 0-3 O 0-40 S 0-2 with an error range of ± 0.5 ppm according to Lechtenfeld and others 2014. Briefly, the following rules were applied:  Kind and Fiehn (2007). Isotopologue formulas (13C, 34S) were used for quality control but removed from the final data set as they represent duplicate chemical information. Relative abundance (RA) values of assigned peaks were calculated based on the sum of all assigned peaks in each mass spectrum. The aromaticity index (AI) was calculated as Dittmar 2006, 2016). In this article we refer to molecular formulas as compounds although each formula may comprise different isomers. Compounds classes were assigned according to Table 1  Highly unsaturated and phenolic high OC B 0.5 \ 0.5 \ 1.5 Highly unsaturated and phenolic low OC B 0.5 \ 0.5 \ 1.5

Statistics
Statistics were performed using the R core program and specific statistical packages where stated (3.5.1. 2018). The Welch's two sample one-sided t-test was used to test whether the DOM reactivity leached from HR litter was higher than that of LR litter in Fig. 3. Tests were done for each litter decomposition time.
The two-sided Welch's t-test was used to test which components changed significantly in RA during litter decomposition (Fig. 4). A non-dimetric dimensional scaling (NMDS) was performed using the vegan package on the normalized RA of individual formulas (Fig. 5). Linear regression modeling was done using the standard R functions on non-transformed ( Fig. 4) and log-log transformed ( Fig. 5a:b) data. Partial least squares (PLS) analysis was done using the pls and spls packages on HR and LR leachates separately with BDOM as predictor variable. First, only formulas occurring in 50% of the dataset were selected. For both models the 3 component model was selected. Variable selection was done using the sparse PLS method, a method designed to deal with multicollinearity often seen in large genetic datasets (Chun and Keleş 2010). For the final models, the O-scores method was used together with cross-validation. The directional effect of individual molecular formulas on BDOM of the leachates was extracted from the model and simplified to negative or positive for Fig. 8 and Table 3.

Litter mass loss and DOM leaching
The cumulative mass loss of HR litter reached 4.1 ± 1.2%, 11.9 ± 4.1% and 20.6 ± 3.6% after 10, 48 and 180 days respectively. The mass loss of LR litter was relatively much lower at 0.8 ± 2.1%, 1.6 ± 2.1% and 6.2 ± 0.5% after 10, 48 and 180 days, respectively. Daily mass loss rates decreased from 0.4% day -1 during the first 10 days to 0.07% day -1 between day 48 and 180 for HR species. For LR species a decrease was seen from 0.08% day -1 during the first 10 days to 0.03% day -1 b Fig. 2 Fraction remaining DOM over incubation time as modelled by the reactivity continuum (Koehler et al. 2012) model. HR litter species are in yellow, LR in blue. Different lines denote different litter field decay times (0 to 180 days) A B Fig. 3 The decay coefficient k (a) and the total bioavailability of DOM (b) for HR and LR litter leachates related to litter decay. Significant differences between HR and LR leachates are denoted (*p \ 0.05, •p \ 0.1) between day 48 and 180. DOM leaching for HR litter decreased from 81.7 ± 42 mg g -1 to 7.4 ± 3.9 mg g -1 for fresh and 180 day decomposed litter respectively. For LR litter DOM leaching was consistently low at an average of 4.9 ± 4.1 mg g -1 for fresh and 3.4 ± 2.5 mg g -1 for 180 decomposed litter. Ranges given depict the standard deviation.

DOM decay dynamics and BDOM
The DOM decay followed the reactivity continuum model (Eq. 2), showing high initial decay rates that decrease over time. The overall DOM reactivity was high across all fresh litter sources and decreased over time as litter was decomposed, with the exception of small initial increases after 10 days in some species (Fig. 2). Overall, only two out of three wood species Ratio of hydrogen to carbon (H/C) and oxygen to carbon (O/C), double bond equivalent (DBE) and number of CHO, CHNO, CHOS and CHNOS peaks (%). Numbers are based on 3 h and 5 LR species, one standard deviation is given in brackets had lower reactive DOM than non-wood litter species (Fig. 2). However, at the start of the DOM incubation the initial decay rate (k) of fresh and 10-day decomposed HR litter was significantly higher than that of LR litter samples (Fig. 3a). This difference disappeared as litter decay progressed. A similar, though less significant pattern was found for the BDOM after 180 days of incubation (Fig. 3b).

DOM chemistry
The FT-ICR MS analysis resulted in on average 2613 to 4270 assigned molecular formulas per timepoint (Table 2). Less formulas were assigned for HR than LR litter samples. The number of assigned formulas increased after DOM incubation, although only slightly for LR (Table 2). In HR litter the assigned formulas increased during litter decay, but for LR litter this was not seen. The double bond equivalent (DBE) of HR litter decreased during litter decay and the amount of nitrogen (N) containing formulas increased. In addition, after incubation the number of N containing formulas was larger (Table 2), while the average mass of formulas was slightly lower before incubation.
As an example, spruce wood (LR species) and Vaccinium litter (HR species) are plotted on van Krevelen plots (Fig. 4). Only few compounds increase in intensity during litter decay for Vaccinium, while most compounds decrease in intensity. For spruce wood an increase and decrease can be seen across the plot, although the increasing formulas are spread along a wider OC gradient. Differences in C:H:N:O:S are shown in supplementary Figs. 2 and 3.
The relative abundance (RA) of compounds in leachate DOM samples was dominated by phenolic and aliphatic compounds (Fig. 5). Unsaturated and polyphenolic compounds were seen at lower RA followed by condensed aromatics. For HR litter, decomposition of the source material lead to a relative increase of aliphatic (p \ 0.05, R 2 = 0.33) and a decrease of high OC phenolic (p \ 0.05, R 2 = 0.36) RA (Fig. 5a). In a similar manner, no significant change was found in the composition of LR leachates. While for fresh (0 and 10 day decomposed) litter a significant difference was found between HR and LR aliphatic (p \ 0.01), phenolic high OC (p \ 0.01) and phenolic low OC (p \ 0.05), in aged (48 and 180 day decomposed) litter only phenolic high OC (p \ 0.05) compounds showed a significant difference.
These results were consistent with a nonmetric dimensional scaling model (stress = 0.042 k = 2), which showed greater distance-and thus greater change in chemistry-between sequential HR litter decomposition samples than for the LR counterparts (Fig. 6). Leachates from more decomposed HR litter were in general in closer vicinity to LR litter samples, showing greater resemblance between HR and LR leachates as litter decomposed. Changes in position were correlated with changes in DOM leaching (p \ 0.01) and litter mass loss (p \ 0.05), showing that greater changes in chemistry occurred for litter species undergoing more rapid changes in litter structure (Fig. 6b, c). Further, it shows a relatively clear separation between LR litter with low leaching, mass loss and changes in chemistry and HR litter with high leaching, mass loss and changes in chemistry. The RA of unsaturated and phenolic high OC compounds decreased for fresh HR litter while the RA of unsaturated and phenolic low OC and aliphatic compounds increased (Fig. 7a). For aged HR litter phenolic RA decreased less or only slightly increased, while aliphatic RA decreased (Fig. 7a). The RA of LR litter peaks was relatively stable for phenolic compounds and decreased for aliphatic compounds (Fig. 7b). Results for polyphenolic compounds were variable, while condensed aromatic compounds in general increased slightly in RA for both LR and HR species (Fig. 7).

The relationship between DOM chemical change and BDOM
For both HR and LR litter, a 3-component PLS model was selected using 86.93 and 94.59% of the variance in the FT-ICR MS data to explain 95.65 and 95.92% of the variance in BDOM respectively. Compounds which increased in RA with increased BDOM (Fig. 8a) were associated with either recalcitrance or microbial products formed during the decomposition process. Recalcitrance in this study is defined as molecular formulas that show a relative increased intensity over the incubation period. Lability is defined as molecular formulas that show a relative decreased intensity over the incubation period. Compounds which decreased in RA with increased BDOM (Fig. 8b) were associated with a higher relative lability. While only a marginal difference between HR and LR could be seen for recalcitrance associated molecules (Fig. 8a), a comparative clear difference was seen for molecules associated with lability, which showed relatively more oxygen compared to hydrogen and carbon (Fig. 8b). Both for HR and LR samples, molecules associated with recalcitrance were in general limited to low O and comparatively lower DBE molecules, while molecules associated with lability showed high O content with a large DBE range (Fig. 8c, d).
Recalcitrance was associated with aliphatic and unsaturated phenolic low OC compounds for HR litter (Fig. 8c). The HR PLS model was driven by compounds increasing in intensity over incubation rather than decreasing (Fig. 8c, d, Table 3). Lability was associated primarily with phenolic compounds and a lesser degree with aliphatic compounds. The LR PLS model was more equally driven by increasing and decreasing compounds than the HR model (Fig. 8g, h, Table 3). Here, low OC phenolic compounds were associated with recalcitrance and aliphatic and high OC phenolic compounds with lability.

Discussion
In line with our hypotheses, our results showed a clear difference in DOM lability between leachates from HR (i.e. birch, alder, Vaccinium) and LR (i.e. spruce, pine, wood) species. Specifically, we found a relatively higher DOM lability of HR litter that decreased as litter decay progressed, thereby eliminating the HR and LR difference (Fig. 3). Changes in DOM chemistry were positively related to changes in mass loss and DOM leaching (Fig. 6) and more profound in birch, alder and Vaccinium than spruce,   Shown are percentages of the total effect pine and wood species over the same litter decay period. Towards the later stages of litter decay the chemistry of LR and HR leachates showed a higher degree of similarity (Fig. 5, 6). Over the DOM incubation period we saw clear differences in bioavailable compounds between fresh HR and LR samples (Fig. 7, 8). Furthermore, there was a high abundance of high OC phenolic compounds associated with the high reactivity seen in fresh HR litter. These compounds seem to be derivatives of increased lignin degradation in HR source litter. These findings show a diversity in leachate chemistry and dynamics from different litter sources that so far has been understudied and stress the importance of highly reactive litter in soil biochemical processes.

Litter leached DOM chemistry and litter decay
Phenolic high OC compounds were highest in fresh HR leachates and subsequently decreased as HR litter decay progressed. Oxidized phenolic compounds are a product of lignin degradation (Klotzbücher et al. 2011a) as lignin undergoes gradual oxidative transformations during the degradation process (Kögel-Knabner 2002). A decreasing yield of lignin-derived phenols is a sign of strong lignin degradation, that has been linked to a high supply of available C (Klotzbücher et al. 2011b). Alternatively, phenolic compounds can also be found in certain tannins (Lorenz et al. 2007) and a rapid decrease of tannins has been seen in spruce litter (Lorenz et al. 2000). However, it has been pointed out that the decrease in tannins found in the latter study might be an artifact due to the formation of non-extractable complexes or structural changes that prevent their detection as litter decay progresses (Kögel-Knabner 2002). In addition, tannins can have strong antimicrobial effects (Kögel-Knabner 2002). However, the highly oxidized phenolic compounds extracted in our fresh HR litter samples decreased significantly during DOM incubation. Thus, the oxidized phenolic compounds are likely lignin derivatives produced during early litter lignin degradation in HR litter as a result of the high C leaching rates.
While the RA of aliphatic and phenolic low OC compounds were found to increase during HR litter decay, this is arguably a secondary result of the decreasing RA of phenolic formulas with high OC and consequently increase of the other RAs. Aliphatic compounds have been shown to accumulate during litter decay (Hempfling et al. 1987;Zech et al. 1987;Bonanomi et al. 2015;Cepáková and Frouz 2015), indicating a lower relative leachability than other compounds. As litter decay progresses high leachable substrates are depleted first (Lorenz et al. 2004), resulting in a proportional increase in the RA of low leachable compounds such as aliphatics. This occurs in the context of quantitatively smaller litter DOM leachates, i.e. the total amount of DOM leached is smaller. As a result, in absolute terms, aliphatic compounds are still leached at higher quantities in fresh HR than decayed litter. For LR litter no such relationships were found. This indicates that no additional lignin degradation takes place due to high C availability as LR litter shows low C leaching rates throughout the decay experiment and because no preferential leaching of certain compounds takes place. It should be noted that these results show DOM dynamics during initial litter decay only. Although lignin has a chemical structure that makes it more resistant to decay, ultimately it is broken down and oxidized (Grandy and Neff 2008), possibly resulting in the leaching of lignin derivatives in later phases of the litter decay process.

Increased chemical similarity through microbial decomposition
Microbial processing of DOM leads to an increased chemical in soils (Strid et al. 2016;Thieme et al. 2019), rivers (Mosher et al. 2015) and oceans (Lechtenfeld et al. 2015;Mentges et al. 2017). Our results show that this process starts already at the primary source-aboveground with litter microbial processing increasing the similarity of leachates between litter types. Similarly, litter chemistry converges as litter decay furthers (Preston et al. 2009b, a). This highlights the uniform direction of carbon decay processes and indicates that an increase in molecular homogeneity in soils, inland waters and oceans is not necessarily the result of increased mixing of different chemically unique sources. Instead, it is rather an effect of the accumulation of chemically similar decay products and the disappearance of unique molecular formulas as decomposition furthers.
However, while differences between species decreased, the number of molecules found in litter DOM extracts increased with HR litter decay. Thus, while microbial processing of litter leads to increased similarity in leachates between sources, at the same time it results in an interspecific increase in chemodiversity. The latter might be explained by the presence of a few high RA peaks in fresh HR litter, which could prevent the detection of less abundant compounds that as a result fall under the detection limit. Still, this suggests a less diverse composition, similar to how high abundancies of one species can prevent the detection of rare species and results in lower diversity scores (Chao and Shen 2003). It should be pointed out that while the SPE carbon recovery rate was within bounds of previous research (20-60% (Raeke et al. 2016), it varied between samples and was often higher at more decayed litter than fresh leachates (Suppl. Fig. 1). A lower SPE carbon recovery might indicate that a larger fraction or different types of molecules were not accessible to FT-ICR MS.
The chemistry of labile vs. recalcitrant litter leachates A higher BDOM was associated with a greater reduction of highly oxidized phenolic compounds in HR litter, while for LR litter it was associated with the reduction of aliphatics (Fig. 7, 8; Table 3). In soils an increase in aliphatics and a decrease in phenolics can be seen in DOM with increasing depth, although the reason for this is thought to be related to sorption mechanics. Phenolic compounds are sorbed to the soil matrix, while aliphatics remain in soluble form (Kaiser et al. 2004;Traversa et al. 2008). Our results suggest that part of the decrease in phenolic compounds seen in soils could also be due to the decomposition of the highly oxidized labile DOM leached from HR litter species. However, in situ detection of these compounds, even in the litter layer of soils, is highly situational due to the rapid decomposition of fresh HR litter leachates (Fig. 2).
A higher BDOM was associated with a greater increase of aliphatics and low oxidized phenolics for HR and LR litter leachates respectively that remained persistent over the incubation period ( Fig. 8; Table 3). The chemistry of the DOM after incubation might also reflect microbial products, as overall there were high DOM decay rates (Fig. 2) and large shifts in chemical composition. This is especially true for the fresh and 10-day decomposed HR litter leachates. As such the increase in aliphatics, as related to HR litter BDOM, might be related to the production of secondary microbial compounds, as substantial amounts of aliphatics can be found in bacteria (Kögel-Knabner 2002). In comparison, the increase of low oxidized phenolic compounds in LR litter is more likely accumulation of undecomposed litter leachates as phenolic compounds are associated more with plant than bacterial or fungal tissues (Kögel-Knabner 2002).

Importance for soil biogeochemical processing
The flux of DOM from litter to soils determines to a large extent terrestrial ecosystem C sequestration (Kalbitz et al. 2000). The bioavailability of DOM suggests that HR litter leachates are more easily processed by microbes. On a short timescale this might lead to a faster return of CO 2 through higher mineralization of the DOM. However, long term C sequestration is linked to microbial processing of C and as the DOM of HR species has a higher bioavailability this might lead to increased persistent SOM formation through the so called DOM-microbial pathway (Cotrufo et al. 2013(Cotrufo et al. , 2015. Microbially processed C has been shown to be a very stable component of SOM (Marschner and Kalbitz 2003;Miltner et al. 2012). In addition, aliphatic components have been shown to accumulate in soils and can persist for millennia (Lorenz et al. 2007). Our results indicate the production of stable aliphatic compounds as fresh HR litter leachates are degraded. In contrast, supplying subsoils with fresh plant derived C can stimulate the mineralization of aged C (Fontaine et al. 2007). However, as DOM moves through the soil it is degraded, transformed and sorbed along the way , meaning that the labile C of fresh litter leachates is less likely to arrive at subsoil depths unaltered.
Since birch, alder or Vaccinium litter initially produce far more DOM than spruce, pine or wood (i.e. LR) litter, the lability and chemistry of these leachates might have an overall larger effect on soil biogeochemistry than the abundant coniferous litter species. In our boreal research area more than 90% of potentially fresh litter leached DOM can be attributed to birch, alder and Vaccinium, while producing less than 50% of the total litterfall (Hensgens et al. 2020). This difference is bound to seasonality as for these species litter C leaching decreases sharply with increased decomposition, while for coniferous and wood litter C leaching stays more or less equal during the first 180 days. Moreover, previously we showed that approximately 70% of fresh litter leached DOM comes from the understory species V. myrtillus (Hensgens et al. 2020). The chemistry of V. myrtillus is most differentiating from the rest, showing the highest amount of oxygenized phenolic compounds released of all species. As such, the understory species V. myrtillus might have a large effect on soil biogeochemistry where it is abundant and possibly lead to higher soil C sequestration rates.

Conclusion
Our results show that following litterfall, litter decay and leaching dynamics follow different paths depending on litter species and more importantly litter reactivity. The reactive birch, alder and Vaccinium (HR) litter produces quantitatively more leachable and labile DOM. As microbial processing of litter decomposition progress, the degree of similarity in chemistry and bioavailability of the leachates increase. In addition, DOM of different litter species also increase in chemical similarity as it is decomposed. This reiterates the importance of microbial processing, not mixing of different terrestrial DOM sources, in shaping the increasingly homogenous chemistry of DOM found in soils, surface waters and oceans. The high lability of DOM across litter sources in combination with the uni-directional changes in chemistry offer an insight in microbial processing and explains why the distinct chemistry of litter sources have little effect on long-term evolving soil or terrestrially exported DOM chemistry. Environmental Research -UFZ and the funds for ProVIS establishment provided by EFRE (Europäischer Fonds für regionale Entwicklung) and the federal state of Saxony (Freistaat Sachsen). The support from Krycklan via SITES is also appreciated.
Funding This research was supported by core funding from Formas (#239-2014-698), the Royal Physiography Society Lund, Å forsk as well as the funding for the ProVIS Centre for Chemical Microscopy, which was established with the European Regional Development Funds (EFRE-Europe funds Saxony) at the Helmholtz Centre for Environmental Research and was supported by the Helmholtz Association. Additional support from Krycklan via SITES.
Data availability The data that support the findings of this study are available from the corresponding author upon reasonable request.

Declaration
Conflicts of interest The authors declare no conflict of interest.
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.