Tree Species of Wet Tropical Forests Differ in Their Tissue Biochemistry and Effects on Soil Carbon Dynamics

Given the hypothesized effects on soil organic matter (SOM) of polyphenols in plant tissues, differences among tree species in their biochemical composition could influence the turnover and accrual of SOM in various ways. The extent to which the biochemical composition of leaf and fine-root tissues differ among tropical tree species, and the effects on soil dynamics, are largely undocumented, however. We used cupric oxide analyses of plant tissues and soil in long-term, replicated, mono-dominant 15-year-old plantations at La Selva Biological Station, Costa Rica, to test for differences among six tree species. We related these results to companion studies in this experimental site to evaluate relationships between interspecific differences in tissue biochemistry and SOM dynamics. Newly senesced leaves and fine roots of the six species differed in their concentrations of three lignin-derived families of phenols, the cinnamyls, syringyls and vanillyls (p < 0.0001 for all tests). Cinnamyl and syringyl phenols in soil differed significantly among species (p = 0.0408, 0.0071, respectively), whereas vanillyl phenols did not (p = 0.83). The degree of decomposition of syringyl and vanillyl phenols in soil also differed (p = 0.0015, 0.0027, respectively), as evidenced by the ratio of carboxylic acid to aldehyde compounds, based on the concept that carboxylic acids are a common by-product of oxidative decomposition of lignin by microorganisms. In our study in a single site, i.e., the same soil type, climate, and growth form of vegetation, total phenols in soil ranged from 5 to 21 mg g–1 organic carbon (OC) across the 20 plots, and the endpoints were both broad-leaved evergreen species; even the means across species, 7–12 mg g–1 OC, covered half the range of values reported in another study across a broad latitudinal range of sites. This study’s tree species differed in traits that influenced at least four factors that explained their differential effects on soil organic carbon (SOC) pools: (1) Fine-root detrital inputs; (2) Fine-root syringyl concentrations; (3) Soil pH; and (4) Macroaggregate structure. This trait-based approach provides a process-based understanding of how trees species influence SOC dynamics, and the consequences for ecosystem properties, under land-use change that involves shifts in species composition.

Given the hypothesized effects on soil organic matter (SOM) of polyphenols in plant tissues, differences among tree species in their biochemical composition could influence the turnover and accrual of SOM in various ways. The extent to which the biochemical composition of leaf and fine-root tissues differ among tropical tree species, and the effects on soil dynamics, are largely undocumented, however. We used cupric oxide analyses of plant tissues and soil in long-term, replicated, mono-dominant 15-year-old plantations at La Selva Biological Station, Costa Rica, to test for differences among six tree species. We related these results to companion studies in this experimental site to evaluate relationships between interspecific differences in tissue biochemistry and SOM dynamics. Newly senesced leaves and fine roots of the six species differed in their concentrations of three lignin-derived families of phenols, the cinnamyls, syringyls and vanillyls (p < 0.0001 for all tests). Cinnamyl and syringyl phenols in soil differed significantly among species (p = 0.0408, 0.0071, respectively), whereas vanillyl phenols did not (p = 0.83). The degree of decomposition of syringyl and vanillyl phenols in soil also differed (p = 0.0015, 0.0027, respectively), as evidenced by the ratio of carboxylic acid to aldehyde compounds, based on the concept that carboxylic acids are a common by-product of oxidative decomposition of lignin by microorganisms. In our study in a single site, i.e., the same soil type, climate, and growth form of vegetation, total phenols in soil ranged from 5 to 21 mg g −1 organic carbon (OC) across the 20 plots, and the endpoints were both broad-leaved evergreen species; even the means across species, 7-12 mg g −1 OC, covered half the range of values reported in another study across a broad latitudinal range of sites. This study's tree species differed in traits that influenced at least four factors that explained their differential effects on soil organic carbon (SOC) pools: (1) Fine-root detrital inputs; (2) Fine-root syringyl concentrations; (3) Soil pH; and (4) Macroaggregate structure. This trait-based approach provides a process-based understanding of how trees species influence SOC dynamics, and the consequences for ecosystem properties, under land-use change that involves shifts in species composition.

INTRODUCTION
Humid tropical forests account for ∼19.6 million km 2 of the Earth's area (Asner et al., 2009) and contain some of the Earth's largest soil organic carbon (SOC) stocks (Field et al., 1998;Jobbágy and Jackson, 2000). Development of these SOC stocks has been, and continues to be a function of climate, organisms, parent material, topography, time, and human factors (Jenny, 1941). Of all these factors, we focus on organisms, specifically tree species, because trees are the dominant growth form in tropical forests, which have been undergoing unprecedented rates of land-use change over the last few decades. These changes include fragmentation, degradation, regeneration of secondary forests following deforestation, and logging of remaining forests (Lewis et al., 2015). These changes generally involve shifts in tree species composition, yet it is not clear how these shifts influence SOC stocks.
Tree species can differ in a variety of traits that could influence the quantity and turnover of organic matter (OM) inputs to soil, that would in turn influence SOC stocks differentially as a result of shifts in species composition. These traits include the rate and chemistry of detrital inputs from litterfall and roots, the decomposition of these plant tissues in soil, and the fueling of microbial and fungal populations that influence SOC dynamics and stocks.
Evaluating plant traits that influence tissue decomposition, and tracing the flow of plant-based carbon (C) into soil and its stabilization there, has been a topic of much debate because so many factors can potentially contribute to the persistence of OM inputs to soil. In their review of multiple studies Schmidt et al. (2011) concluded that C residence time was only secondarily influenced by the molecular structure of detrital inputs; rather, the persistence of OM in soil was the result of complex interactions between the biota (including microbes and fungi) and abiotic factors such as climate, mineralogy, soil acidity, and water availability. In a litter-manipulation experiment, Ma et al. (2019) found that differences in temperature sensitivity of the active SOC pool were related to lignin content and the microbial community, whereas for the slow pool temperature sensitivity was related to environmental factors such as pH and soil C/N ratio.
To evaluate the relative impacts of the quantity and biochemistry of plant detrital inputs on SOC dynamics and stocks, we conducted studies in a field experiment in lowland Costa Rica that contains monocultures of six tropical tree species in four blocks, with species selected such that they differed in their root and litter attributes (Fisher, 1995). Established in 1988, this is one of the oldest replicated experiments containing broad-leaved evergreen tree species in the humid tropics. Because all tree species were planted in plots on the same soil type under similar environmental conditions and land-use history, this study focuses on the effects of the tree species on accrual of SOC after 15 years. Previous studies in this site allow for a rich database of traits to evaluate in relation to SOC stocks, including fine-root growth, litterfall, microbial studies, and various soil studies. For the study reported herein, we measured lignin-derived phenols (i.e., biomarkers that could be measured in plant tissues and soil), and evaluated SOC pools and their turnover times using a two-pool model. Given that lignin is a major constituent of forest litter, biomarkers such as cinnamyl, syringyl, and vanillyl phenols have been used to identify broad categories of vegetation in forest litter, e.g., angiosperms vs gymnosperms (Kögel, 1986), or pasture vs forest (Heim et al., 2010). The oxidative degradation of lignin with cupric oxide (CuO) has also been used to characterize the pattern of lignin decomposition (Kögel, 1986;Feng and Simpson, 2011). The acid-to-aldehyde (Ad/Al) ratio of syringyl and vanillyl units has been found to increase with the degree of lignin degradation (Kögel, 1986;Hedges et al., 1988;Opsahl and Benner, 1995;Otto and Simpson, 2006).
To compare and integrate findings about plant tissue biochemistry with SOC dynamics, long-term laboratory incubation-based methods can be used to estimate active and slow turnover pools of SOC and their residence times (Paul et al., 2001). This method has been used to compare SOC dynamics and C sequestration under different soil types and land uses and management (e.g., Six et al., 2002;McLauchlan et al., 2006), although to our knowledge laboratory incubations have not yet been coupled with biochemical analyses of tropical forest SOC.
We addressed four broad hypotheses in relation to differences in tropical tree species traits and their influence on SOC dynamics and stocks: 1. Species differ in the biochemistry of their newly senesced leaves and fine roots, and the differences persist in SOC. While differences between angiosperms and gymnosperms or grasses and trees are well-documented (Kögel, 1986;Heim et al., 2010), we are not aware of any other studies comparing broad-leaved evergreen tropical trees. 2. Species differ in the decomposability of their fine roots and soil, as indicated by the Ad/Al of their syringyl and vanillyl concentrations in soil. While this was true for studies comparing different growth forms and forest types (Guggenberger et al., 1995;Heim et al., 2010), it has not been tested among tropical tree species. 3. The size and turnover of active and passive SOC pools differ among tree species. In a previous study in this site, C mineralization rates differed among the species, as determined from 30-day laboratory incubations . However, we have not previously evaluated sizes of the active, slow and passive SOC pools and their turnover. 4. Differences among species in both the quantity and biochemistry of detrital inputs influence the dynamics and stocks of SOC. Results from previous studies in this site indicate that the tree species differ in the quantity and decay of fine-root inputs to soil Raich et al., 2009). However, we have not previously addressed the relative effects of the quantity and biochemistry of detrital inputs and their relationship with the sizes and turnover of SOC pools.

Study Site
This study was conducted in a lowland Tropical Wet Forest, sensu Holdridge life zones (Hartshorn and Hammel, 1994), within experimental plots established in 1988 at La Selva Biological Station, Costa Rica. Mean annual precipitation (MAP) is 3,960 mm, with no month receiving < 100 mm; mean annual temperature (MAT) is 25.8 • C (Sanford et al., 1994). The soil has most recently been classified as an Oxisol, a Mixed Haplic Haploperox (Kleber et al., 2007). The site is hilly, with elevations ranging from 51 to 84 m above mean sea level . Mean soil temperature (5 cm depth) was similar to air temperature (10 cm above forest floor), 24.4 and 24.1 • C, respectively, based on hourly measurements over 24-h periods on >100 days from 2004 to 2010 in these plots (Raich and Valverde-Barrantes, 2017). Of these >2,600 soil temperature observations, 50% fell within a range of −1.3 to +2.0 • C of air temperature. The field experiment was established on a 12-ha site that had been mature wet tropical forest until deforestation in 1955. The site was then converted to cattle pasture dominated by the perennial grasses Panicum maximum L. and Melinis minutiflora Pal. and grazed until abandonment in 1987. In 1988 a randomized complete block experiment was established. The design contained twelve 0.25-ha (50 m × 50 m) plots in each of four blocks; each plot was divided into four quadrants to facilitate stratified random sampling within plots. Each of the 11 plots was planted with a single tree species, with a tree spacing of 3 × 3 m; the 12th plot was left unplanted as a control (Fisher, 1995). For details about site management, and early tree growth and soil properties, see González and Fisher (1994);Fisher (1995); and Haggar et al. (1997). To provide a reference for soil property measurements, we established four 50 m × 50 m plots in the adjacent mature forest, and four 14 m × 8 m plots in the adjacent abandoned pasture.

Plant and Soil Sampling
Newly senesced, fallen leaves were collected from the ground in March 2004. The status of "newly senesced" was determined based on leaf color and integrity, as described by Raich et al. (2007). At least 20 leaves (fascicles in Pinus) were collected for each species in each plot and combined within a plot to yield n = 1 sample per plot. Fine roots (0-2 mm diameter) were sampled in June-July 2004 by carefully excavating and tracing roots to six individual, randomly selected trees within each plot. The six samples per plot were bulked to produce a single sample in each of the 24 plots. Root samples were washed before drying. All plant tissue samples were dried at 65 • C, ground in a Wiley mill, and then pulverized to a fine powder using a LABCONCO Spex 8000 Mixer/Mill.
Fine root growth over this same period had been measured in a previously reported study . To summarize, fine-root growth (0-15 cm depth) was measured using the ingrowth core method, using 20 cores per tree species. Fine-root growth was assumed to equal fine-root production because fine-root biomass did not vary significantly 2-year timeframe of that study.
Soil was sampled in the 0-15 cm layer, following the protocol of Russell et al. (2007) in June-July 2003. Within each quadrant of each plot, 20 punch-tube samples, 3.2-cm in diameter, were collected at randomly selected points and then bulked into a single sub-sample, for a total of four sub-samples per plot. For this study, the four sub-samples per plot were combined into a single sample per plot, for a total of 24 soil samples. These soil samples had been air-dried and roots were removed as the soil was passed through a 2-mm sieve. A well-mixed sub-sample was finely ground with mortar and pestle for subsequent analyses. Soil sub-samples were dried at 105 • C to determine the air-dry to dry-weight conversion factor. By 2008, it had become apparent that V. ferruginea was succumbing to a fungal disease. All of the tissue concentrations analyses we report here are for live tissues, and thus would not have been affected by the disease. However, we could not rule out the possibility that this disease could have caused premature fine-root die-off earlier than 2008, and could have thus created an artifact in the fine-root growth and soil data. Thus, we did not include V. ferruginea in any analyses involving soil or fine-root growth. This study was not designed to compare differences in land use, but for reference purposes, we used the same protocols to sample soil in the adjacent mature forest and abandoned pasture.

Lignin-derived Phenols
Cupric Oxide extractable, lignin-derived phenols were measured by adapting the procedure of Goñi and Hedges (1992). Samples of 100 mg soil or 6 mg roots were oxidized with alkaline CuO and 2 M NaOH for 3 h in a stainless steel minibomb. Bombs were loaded into a specially designed carousel in a refitted Hewlett-Packard 5890A gas chromatograph (GC). The heating rate was 4.2 • C/min for 30 min from a start temp of 27 • C, and a final temp of 150 • C was held for 150 min. Following heating, ethyl vanillin in a concentration of 240 ng µl −1 (25 µl for soils, 100 µl for roots) was added to each sample as the internal standard. The minibombs were closed again, shaken by hand, then centrifuged for 5 min at 3,500 rpm, and their supernatant transferred to test tubes. Three ml of 6 M HCl was added to each bomb to neutralize the NaOH and each test tube was vortexed. About 2.5 g NaCl was added to each test tube, then capped and mixed by hand. About 3 ml ether was added, then again capped and mixed by hand. Test tubes were shaken for 2 min to accumulate the phenols in the ether, then centrifuged for 2 min at 2,000 rpm. The ether layer containing the phenols was then carefully pipetted from the surface of each solution into a new test tube. Na 2 SO 4 (∼2.5 g) was added to each tube to remove residual H 2 O from the ether solution. Tubes were shaken, and if necessary stored under refrigeration for overnight.
For derivitization and chromatography, samples were first filtered through quartz wool then concentrated by rotary evaporation to about 1 ml volume, then transferred to 4ml vials, and dried to nearly air-dryness under a stream of Ar gas. The precipitate was then dissolved in 150 µl pyridine/50 µl methyl-3,4-dimethoxybenzoate (soils) or 300 µl pyridine/100 µl methyl-3,4-dimethoxybenzoate (roots). The methyl-3,4-dimethoxybenzoate served as a recovery standard a 50 µl aliquot of each solution was added to a new 4-ml vial and mixed with 50 µl of the derivitizing agent bis(trimethylsilyl)trifluoroacetamide and then placed in a heating block at 70 • C for 30 min, after which the samples were loaded into the GC. Gas chromatography was carried out with a Hewlett-Packard 6890 GC equipped with a derivatized, nonpacked injection liner, a HP-5 (5% crosslinked phenylmethyl siloxane) capillary column (30-m length, 0.32-mm column id., 0.25-µm film thickness) and a flame ionization detector. The following conditions were employed for phenol separation: injector temperature, 300 • C; temperature ramp, 100 • C ramped to 220 • C at 4 • C/min, 220 • C ramped to 270 • C at 25 • C/min and held for 10.1 min; detector temperature 300 • C. Results are reported in the conventional units of SOC enrichment, that is mg phenol g OC −1 .

Soil C Pool Sizes and Turnover Rates
To determine sizes of active and slow pools and their respective turnover rates, we used long-term incubations of soil at 23 • C, similar to field ambient temperature, coupled with a chemical approach, described by Paul et al. (2001). We incubated 16 field replicates of soil for each species, using sieved (2 mm), air-dried subsamples of 3 g, from each of the 4 quadrants × 4 blocks, as described in Russell et al. (2007). To summarize, the subsample was mixed with acidwashed sand (1:1 w/w) and adjusted to 60% water-filled pore space, using deionized water that contained 2 mL of inoculant (5 g field-moist soil, combined from all plots, homogenized with 50 mL distilled water). We measured CO 2 -C release rates periodically over 400 day (before CO 2 concentrations reached 4%); the flask headspace was flushed through an infrared gas analyzer (LI-820, LI-COR Biosciences, Inc., Lincoln, NE, United States). The air supply was piped through a CO 2 -free water bath to keep the incubation-tube atmosphere hydrated. The mean lab temperature was very close to ambient field temperature, so we did not perform a temperature conversion of the data.
We determined the resistant fraction by refluxing 1-g soil subsamples in 6 M HCl at 115 • C for 16 h on a controlledtemperature digestion block as described by Paul et al. (2001).
The refluxed samples were then washed three times with deionized water, dried at 55 • C and finely ground in a mortar and pestle. Resistant soil C was measured by dry combustion, using a Thermo-Finnigan EA Flash (Series 1112, EA Elantech, Lakewood, NJ, United States).

Soil Pool Parameter Calculations
We used curve analyses of CO 2 evolved for each sample incubated over 400 days (n = 96) to determine the size and turnover rates of the active (C a ) and slow (C s ) SOC pools as in Paul et al. (2001). We did not have carbon dating for these soils, so we used a two-pool constrained model in the regression analysis: is the total C pool at time t, C a and k a are the size and decay rate of the active pool. C s is the size of the slow pool (SOC -C r -C a ), where C r is the size of the resistant pool as determined by acid hydrolysis and k s is the turnover of the slow pool. The mean residence times (MRT a and MRT s ) were calculated as 1/k for each of the two soil pools. We used nonlinear regression to calculate the parameters C a , k a , and k s , using PROC NLIN in the Statistical Analysis System (SAS Institute Inc [SASI], 1999). A Marquardt fit was used with initial parameter settings of: k a = 0 to 0.3 by 0.05; k s = 0.001 to 0.005 by 0.001; and C a = 100 to 2,000 by 100.

Analysis of Variance
The plot was the experimental unit; for variables with measurements at the quadrant level, the mean value for the plot was used in analyses. Soil data from V. ferruginea and the reference plots (mature forest and abandoned pasture) were not included in the analyses. To evaluate whether the tree species had effects on more than one response variable among soil phenols and SOC (pools and total), we first conducted a multivariate analysis of variance. Interactions among the response variables were significant, so we then tested response variables individually, using a generalized linear mixed effects model with both fixed and random effects in SAS (Littell et al., 1996). We tested for homogeneity of variances and normality of distributions. Variables that did not meet these assumptions included root and soil cinnamyl, and soil syringyl and vanillyl phenols, which met the assumptions when natural log transformed. Thus, analyses were conducted on the transformed data; we present the untransformed values, however. Differences were considered to be significant at α ≤ 0.05. Pairwise comparisons for significant overall F-tests were performed using p-values adjusted by Tukey's Honestly Significant Difference method.

Multiple Regression Models
We used a multiple regression model to evaluate the relative effects of the quantity, biochemistry and turnover of OM on one response variable, total SOC stocks. We used partial residual analysis to quantify the effect of a single explanatory variable on the response variable, total SOC stocks, after the effects of other possible explanatory variables in the model had been Frontiers in Forests and Global Change | www.frontiersin.org accounted for (Neter et al., 1996). With only n = 20 plots for this analysis, we were limited to only six explanatory variables in the model. There is a rich data set for this field experiment, so to select the six explanatory variables from the set of many possibilities, we used Pearson correlation analyses to evaluate the strength and significance of pair-wise comparisons of variables (Supplementary Table 1). Variables from other microbial and fungal studies in these plots were not included in the set of possible variables because previous studies had indicated a lack of difference among the tree species in microbial and fungal parameters (Kivlin and Hawkes, 2016a,b); also, microbial and fungal variables were not correlated with total SOC stocks (Russell et al., 2018). For plant tissue biochemistry, we evaluated fine roots but not leaves, given that fine roots had been more correlated with SOC stocks than were leaves in another study at La Selva (Russell et al., 2004). Based on the results of the correlational analyses, the six explanatory (=predictor) variables were included in a multiple regression model They were: (a) fineroot syringyls (mg 100 mg OC −1 ); (b) fine-root growth (Mg C ha −1 year −1 ); (c) acid/aldehyde of soil syringyls (mg 100 mg OC −1 ); (d) MRTs (year); (e) active SOC pool (mg kg −1 ); and (f) soil syringyls (mg 100 mg OC −1 ); all variables relate to the 0-15 cm soil layer.

Phenols in Plant Tissues and Soils
Across the tree species in this study, mean total cinnamyl phenol concentrations (mg g OC −1 ) ranged from 1.5-3.4 in leaves, 0.8-3.7 in fine roots, and 0.9-1.7 in soil ( Table 1). Mean total syringyl phenols (mg g OC −1 ) ranged from 0.3-6.7 in leaves, 1.8-13.4 in fine-roots, and 1.5-5.6 in soil. Mean total vanillyl phenols (mg g OC −1 ) ranged from 3.8-16.0 in leaves, 8.3-17.5 in fine roots, and 3.6-5.5 in soil. Thus, for cinnamyl and separately for vanillyl phenols, concentrations in leaves and fine roots were of similar magnitudes. In contrast, the average concentration across all plots of syringyl phenols in fine roots was double that of leaves, 7.8 and 3.8 mg g OC −1 , respectively. Leaf and fine-root syringyl concentrations were significantly correlated (r = 0.62, p = 0.0098), as were vanillyl phenols (r = 0.79, p = 0.0003), but cinnamyl phenol concentrations were not significantly correlated (r = 0.18, p = 0.4936; Supplementary Table 1). The ratio of the phenolic compounds in fine roots-to-soil was calculated for each species. Across species, cinnamyl root concentrations ranged from 0.9 to 1.7 times that of soil, whereas the ratio was higher for syringyl phenols, 1.2 to 3.3 across species and for vanillyl phenols, 2.2 to 3.3.
Regarding comparisons among species, all three phenolic classes differed significantly among the tree species for both leaf and fine-root tissues (Tables 1, 2). For soil, species did not differ significantly in vanillyl phenol concentrations; for cinnamyl and syringyl phenols, the differences among species were significant, but not as strongly as for plant tissues. The sum of cinnamyl, syringyl, and vanillyl phenols (CSV) ranged from 5 to 21 mg g −1 OC across the 20 plots (published dataset), with a mean of 7-12 mg g −1 OC across species (Table 1). In comparison, the CSV sums for soil in the mature forest and abandoned pasture reference samples were 11 and 8 mg g −1 OC, respectively, with Vochysia most closely resembling mature forest values.
The Ad/Al of both syringyl and vanillyl phenols in soil differed significantly among species (Tables 1, 2), with that of Pinus higher than Hieronyma, Pentaclethra and Vochysia for syringyl phenols and Hieronyma, Pentacletha and Pinus greater than that of Virola and Vochysia for vanillyl phenols. In the mature forest, the Ad/Al for syringyl phenols in soil was lower than that of the abandoned pasture, with Pinus and Virola more closely resembling that of the abandoned pasture, and the other three species closer to the mature forest values. There was not a significant relationship of the Ad/Al of syringyl phenols in soil with either root or soil syringyl concentrations (Figure 1 and Supplementary Table 1). Acid/aldehyde ratios of roots and soil were not significantly correlated, for either syringyl or vanillyl phenols; nor was the soil Ad/Al correlated with total SOC for either biochemical family (Figures 1, 2

SOC Pools and Turnover and Fine-Root Growth
The mean size of the C a (active pool) ranged from 880 to 1,281 mg kg −1 across species; the C s (slow pool) was 25 to 32 times larger, ranging from 27,000 to 32,000 mg kg −1 ( Table 3). Both the active and slow pools were significantly larger in Vochysia and Hieronyma. The C a pool was smallest in Virola and the C s pool smallest in Pinus and Pentaclethra. The tree species differed significantly in the turnover of the slow pool, reported here as the mean residence time (MRT s ): it was longest in Virola and shortest in Hieronyma (Table 3). In contrast, the MRT a was similar across species. Fine-root growth differed significantly among the species, ranging from 1.2 to 3.9 Mg C ha −1 year −1 in Pinus and Hieronyma, respectively (Table 3; Valverde- ).

Correlation Analysis Results
Across species, fine-root growth was significantly correlated with the C a , MRT s , and total SOC (Figures 1, 2 and Supplementary  Table 1). In contrast, root syringyl phenols were not significantly correlated with these measures of SOC turnover, or with the Ad/Al of soil syringyl phenols or total SOC stocks. Within the total SOC, the relationship between the pool size and its turnover differed with the pool: C a and MRT a were significantly correlated, but C s and MRT s were not (Supplementary Table 1).

Partial Regression Results
Based on the results of the correlation analyses, six explanatory variables were included in the multiple regression model to analyze their effects on total SOC stocks: (a) Fine-root growth (FRG); (b) Fine-root syringyls (RtSyr); (c) Acid/Aldehyde ratio of Soil Syringyls (AdAlSyr); (d) mean residence time of the slow pool (MRT s ); (e) size of the active SOC pool (C a ); and (f) Soil Syringyls (SoSyr). Of those six variables, fine-root growth, fineroot syringyls and the size of the active SOC pool were significant after the other variables have been taken into account ( Table 4). Comparisons of various multiple regression models indicated that for the simplest models (only one variable), the size of the active pool and fine-root growth explained the most variability Notes. Means ± SE of n = 20 plots for leaves and roots, n = 16 for soil. Means within each comparison followed by the same letter do not differ significantly among species (Tukey's HSD, α = 0.05). Bold letters denote significant differences. Tree species are identified by four letters, the first two letters of the species and genus listed in the text. "Forest" refers to mature forest. "Pasture" refers to abandoned pasture. All root and soil data are for the 0-15 cm depth interval. ( Table 5). However, in 2-variable models, the addition of root syringyls explained substantially more of the variability, without over-fitting the model. The addition of Ad/Al soil syringyls in a 3-variable model had a similar effect. Beyond that, addition of the other variables had diminishing explanatory effects.

Effects of Tropical Tree Species on Phenolic Biochemistry
Lignin signatures based on CuO analysis have been used to distinguish between vegetation types (Hetherington and Anderson, 1998), effects of land-use change (Guggenberger et al., 1994), and invasions of a new growth form, as in shrub encroachment in grasslands in Mongolia (Zhou et al., 2018). Most studies of plant phenolic contents in relation to SOC dynamics take place in the temperate zone or higher latitudes. Moreover, studies of plant tissues usually concern aboveground plant components, despite that root-derived C dominates SOC stocks in comparison with aboveground detrital inputs, as reviewed in Schmidt et al. (2011). Comparison of a cropping rotation, grassland, mixed deciduous forest and spruce forest in a pre-Alpine region, found that the sum of lignin-derived soil phenols, CSV, ranged from 11 to 23 g/kg (Guggenberger et al., 1994). Across 18 sites that encompassed a broad range of latitudes and soil types from Saskatchewan (Canada) to Texas (United States), in which MAT ranged from 0.9 to 23.4 • C and MAP from 300 to 1,308 mm, CSV phenols in soil ranged from 9 to 26 mg g −1 OC (Amelung et al., 1999). In our study in a single site, i.e., the same soil type, climate, and growth form of vegetation, CSV phenols in soil ranged from 5 to 21 mg g −1 OC across the 20 plots; even the means across species,   Data are from long-term incubations of soil, 0-15 cm depth. SOC pools: C a = active; C s = slow; and C r = resistant. Units are: mg g −1 for all soil C pools; day −1 for k a and k s (decay rates); days for MRT a ; years for MRT s ; and Mg C ha −1 year −1 for FRG. Different letters within a row denote significant differences among species (α = 0.05) by Tukey's HSD. There were n = 20 plots in the analyses, with df for blocks = 3 and species = 4. Data for reference vegetation not included in analyses. Data for total SOC were published previously . Fine-root growth data were reported previously in Valverde- . No data = n.d. Bold letters denote significant differences.
7-12 mg g −1 OC (Table 1), covered half the range found in the broad latitudinal study by Amelung et al. (1999). The broad range found in our study was not driven by different growth forms either: values for the only gymnosperm species in the experiment fell in the middle of the range, with the endpoints of the range both represented by broad-leaved evergreen tropical tree species. In a study of 18 sites across a range of biomes, from taiga and tundra to temperate forest and tropical forest, Vancampenhout et al. (2009) concluded that tropical soils were low in lignin. However, two of three forests in that study were dry forests. Low soil lignin content could be explained by lower detrital inputs in the drier forests, given the results of another study in which net primary productivity (NPP) in tropical forests declined with MAP for sites with MAP less than ∼2,500 mm (Schuur, 2003). It is also possible that lignin degraded faster in the drier forests, given the finding that lignin was subject to accelerated degradation via photodegradation in a semi-arid environment because lignin effectively absorbed light over a wide range of wavelengths (Austin and Ballaré, 2010). Thus, the lower lignin  The p values and Type II Sums of Squares) refer to significance of the reduction in the sum of squared error when adding a parameter to the model that already has the other five parameters in it. Bold letters denote significant differences.
contents in tropical dry forests could be the result of an abiotic factor that increases lignin degradation.

SOC Pools and Their Turnover
In our study, the Ad/Al ratio of both syringyl and vanillyl phenols in soil differed significantly among species. For syringyl phenols, this ratio decreased significantly as fine-root growth increased (r = −0.51, p = 0.0210). This suggests that in a species such as Vochysia, which had a low Ad/Al ratio, the soil pool of syringyl phenols was being replenished relatively more with inputs of fineroot detritus. Hieronyma, however, also had relatively high fineroot growth, but a relatively higher Ad/Al ratio, suggesting that its fine-root detritus decomposed relatively faster. Thus, while the Ad/Al ratio provides a measure of the state of decomposition of  OM within the soil, when comparing different treatments (e.g., land use, species, and history), the interpretation of these data requires information about the quantity of detrital inputs to the soil pool in each of the treatments. The size of the active pool (0-15 cm), averaged across all species in our site, 1,078 mg kg −1 ( Table 3) was nearly double that reported for an agricultural soil (0-20 cm) at Kellogg Biological Station in Michigan, 570 mg kg −1 (Paul et al., 2001). The slow, resistant and total SOC pools in our site were > 5-, 3-, and 4fold greater, respectively, than those reported for the Michigan site. The MRT in that site for the active and slow SOC pools were 20 days and 12 years under laboratory temperatures, respectively, compared with 8 days and 10 years for MRT a and MRT s , respectively, averaged across species in this study ( Table 3). When the Michigan data were converted to field temperatures, however, the MRT a was 60 days and MRT s was 37 years. For the tropical soils in this study, the laboratory and field incubation temperatures were nearly the same; thus these tropical soils were more dynamic and had larger SOC pools.

Tree Species Effects on Accrual and Persistence of SOC
Given that lignin is generally considered to be recalcitrant, owing to its complex chemical structure that protects litter from microbial attack (Swift et al., 1979), variation in lignin content among species is expected to translate into speciesspecific variation in the decomposability of their tissues. Indeed, in multiple studies of decomposition of aboveground plant tissues, higher lignin content has been highly correlated with slower decomposition of detritus (Melillo et al., 1982;Berg et al., 2000;Cornwell et al., 2008). Lignins were also expected to degrade slowly in soil, but the evidence regarding lignin's recalcitrance within soil has been varied. A review by Vancampenhout et al. (2009) indicated a lack of persistence of lignins in soil. In a comparison of lignin biomarkers in a pasture and rainforest soil in subtropical Australia, Heim et al. (2010) found that only small portions of rainforestderived lignin persisted 100 years after conversion to pasture, indicating that lignin biomarkers didn't provide a proxy for stable SOC. In tropical soils, higher temperature is thought to increase decomposition of lignin, and thus explain their low lignin contents. As indicated above, however, the soils in our study site had lignin contents on par with temperatezone sites. Mineral-organic associations can stabilize a large fraction of SOC (Mikutta et al., 2009), however, and this could provide a mechanism for stabilizing lignin-derived phenols in our study site's Oxisol. In addition, given the high rates of NPP (Russell et al., 2010), detrital inputs were replenishing soil lignin stocks in this site. This study was not designed to distinguish between these two mechanisms of lignin retention and accrual in soil, however.
Many studies have demonstrated that leguminous trees increase SOC stocks, as reviewed by Mayer et al. (2020). Given our results that the only nodulated legume species in the study, Pentaclethra, had low SOC stocks (Table 3), we suggest that the generalization about tree species with N-fixing root associates may well hold for sites that are relatively limited in available N, but not necessarily in tropical sites where N is relatively less limiting for microbial processes (Camenzind et al., 2018). Our data also suggest that lower inputs from roots was an important factor that explained differences among tree species in the accrual of SOC stocks (Figure 1). In biomarker studies in temperate zone forests, roots have also been shown to be more important contributors to SOC than aboveground detritus (Crow et al., 2009;Spielvogel et al., 2014). Root activities in soil may afford them physical protection through mycorrhiza and interactions with metal ion (Rasse et al., 2005). By allocating C to microbes, roots may also prime microbial processes, and thus speed up decomposition of older SOC and thus, availability of nutrients for plant uptake (Fontaine et al., 2007;Kuzyakov, 2010), thereby stimulating plant growth and the detrital inputs to soil.
Tree species identity, i.e., species-specific effects, have been shown to influence SOC stocks across a range of experiments (Vesterdal and Raulund-Rasmussen, 1998;Hansson et al., 2011;Osei et al., 2021). Nevertheless, in their evaluation of a large set of factors that could provide indicators of SOC storage, Wiesmeier et al. (2019) concluded that there was no consensus among previous studies regarding the effects of tree species on SOC storage. Some studies distinguish between broad taxonomic categories of tree species, and then find differences in SOC storage between the categories, e.g., gymnosperms and broadleaved species, without differences within a category (Mayer et al., 2020). In this study, however, we found differences among the four broadleaved evergreen tropical species in their pools of SOC and turnover times. Kleber et al. (2011), in their comparisons of three different soil types, found that the turnover rate of SOC was not determined solely by the thermodynamic stability of the detrital inputs, such that chemical "recalcitrance" of SOC was not a good predictor of its turnover. Schmidt et al. (2011) emphasized that although compound chemistry is a driver of decomposition, its relative importance must be considered within the context of other environmental factors. Lehmann and Kleber (2015) found that prediction of SOC turnover required an understanding of the interactions between the substrate, microbes, and abiotic drivers; as such, OM quality was a poor predictor of OM persistence. In our study, there was not a significant correlation between total SOC stocks and the Ad/Al ratio of soil, for either syringyl or vanillyl phenols (Supplementary Table 1). These ratios can represent the state of oxidation of these phenolic compounds, and thus indicate the intensity of microbially driven oxidation or the length of time the phenols have been in soil, rather than intrinsic phenol stability. That interpretation is consistent with the results from these experimental plots that microbial composition did not differ among tree species (Kivlin and Hawkes, 2016a,b).
In contrast with phenol Ad/Al ratios in soil, fine-root syringyl concentrations did explain differences in total SOC stocks across species in our study, after other variables were taken into account in the partial residual analysis (Table 4 and Figures 1, 2). This effect of biochemistry on SOC stocks was likely promoted in part by high NPP, and accentuated by the short-term recalcitrance of syringyl phenols; these two factors together would cause syringyl input rates to exceed their degradation rates. Fineroot growth was significantly correlated with fine-root syringyl phenol content (r = 0.71, p = 0.03, and Supplementary Table 1), indicating the importance of input rates. However, of the three families of root phenols in this study, only the syringyl phenols, the "woodiest" of the three, were also identified as a significant factor; this indicates a role of their specific chemical nature or location within roots, in addition to quantity of inputs alone. In contrast, the more degradable cinnamic phenols had much lower accumulations in soil (Table 1), further suggesting that biochemistry played a role in soil organic matter dynamics.
Given the experimental design of our study, however, abiotic factors, e.g., MAT, MAP, parent material, time of soil development, relief, and previous land-use history were similar across plots. Certainly, individual species can influence microclimate and microbial communities that could influence SOC dynamics. Given that the tree species in this experiment had similar canopy coverage and height, microclimate differences were not expected to differ significantly. Other results from this study site indicated that soil bacteria richness, composition and soil fungal abundance, richness, and composition did not differ among the tree species (Kivlin and Hawkes, 2016a,b). Based on this and other studies in this field experiment, we suggest that at least four other factors explained the differences among the tree species in SOC pools and dynamics. (1) Differences among the species in fine-root productivity (Valverde-  and decay rate  resulted in differential inputs of detritus to the soil pool. (2) Differences in fine-root syringyl phenols explained some of the variation in total SOC pools across species (Tables 1, 4). (3) The species differed in their effects on soil pH, which in turn influenced the availability of cations that are limiting in this Oxisol, and hence the quantity and chemistry of detrital inputs to soil (Russell et al., 2017). As a result of Vochysia's trait of accumulating aluminum, this species de-acidified the soil. We hypothesized that increasing the soil pH above the point of zero charge in this soil with variable-charge clays dispersed soil organo-mineral colloids, thereby liberating cations available for plant uptake. In contrast, as a result of its N-fixing association, Pentaclethra reduced soil pH below the point of zero charge, which was hypothesized to reduce cation availability for uptake. (4) The species also differed in their effects on macroaggregate structure (Russell et al., 2018), which would have differentially influenced the protection of OM in soil.
Given the rates and extent of land-use change in humid tropical forests over the last several decades (Asner et al., 2009) that have generally resulted in loss of biodiversity, the need is ever more urgent to understand the effects of tree species on soil C. Our results from these studies indicate that, as Schmidt et al. (2011) suggested, the explanation for differences among the tree species in their effects on SOC stocks and dynamics is complicated and involves a mix of factors. Nevertheless, developing a framework for evaluating the effects of species advances our conceptual framework of carbon cycling (Russell et al., 2010), and offers managers an important tool, selection of tree species, for forest regeneration projects or establishment of plantation. Our study illustrates that a trait-based approach can provide a process-based understanding of how trees can influence SOC dynamics at a finer-tuned taxonomic level than the Division, e.g., Angiosperms vs. Gymnosperms, and the consequences for ecosystem properties such as SOC stocks.

CONCLUSION
We compared tropical tree species in terms of the quantity of detrital inputs and tissue biochemistry and evaluated the effects of these differences on SOC pools and turnover. The study was conducted in a 15-year-old field experiment containing six tree species, including one gymnosperm and one nodulated legume, grown in monocultures in four blocks. Even within the category of broad-leaved evergreen tree species, the species differed in their concentrations of phenols in tissues and soil. In this single site with the same soil type, MAT and MAP, the means of total phenols across the study species encompassed half the range of values found in another study that covered sites across a broad latitudinal range (Canada to Texas). We used partial residual analysis to quantify the effect of a single explanatory variable on SOC stocks; results indicated that differences among species in syringyl phenols in fine roots, and in the quantity of detrital inputs from fine roots explained differences among species in SOC stocks. Previous studies in this site indicated that the tree species also influenced SOC accumulation differentially as a result of species-specific effects on soil pH and macroaggregate structure. Together, these results highlight the complexity of tree species effects on soil, and the need for a trait-based approach in the selection of tree species when ameliorating degraded sites after land-use change.

DATA AVAILABILITY STATEMENT
The original contributions presented in the study are included in the article/Supplementary Material. The datasets for this study are available: doi: 10.5061/dryad.fbg79cntt.

AUTHOR CONTRIBUTIONS
AR and DO conceptualized the study and acquired funding. RM collected the data under supervision of AR and DO. AR analyzed the data, designed and ran the model, and prepared the figures and tables. All authors contributed to interpretation of the findings. AR and RM wrote the original draft, and all authors participated in reviewing and editing the manuscript.