Variation in Feedstock Wood Chemistry Strongly Influences Biochar Liming Potential

Chars intended for use as soil amendment (“biochars”) vary greatly in their chemical and physical properties. In the present study, 19 Canadian temperate wood feedstocks were charred across a range of pyrolysis temperatures from 300–700 ◦C. The resulting 95 biochars were tested for their physio-chemical properties and liming capacity. Data indicated increasing base cation concentrations including Ca, Mg, and K (elements that characteristically form liming compounds, i.e., carbonates) as pyrolysis temperature increased. Acidic surface functional groups were analyzed with modified Boehm titration: Carboxylic and lactonic functional group concentrations decreased and phenolic group concentration increased with pyrolysis temperature. Functional group composition also varied greatly with feedstock: In particular, conifer-derived biochars produced at pyrolysis temperatures <500 ◦C showed much higher carboxylic and lactonic functional group concentrations than did angiosperm-derived biochars. Liming capacity was assessed using soil incubation experiments and was positively related to biochar pH. Both acidic surface functional group concentration and nutrient element concentration influenced biochar pH: we developed a non-linear functional relationship that predicts biochar pH from the ratio of carboxylic to phenolic moieties, and concentrations of Ca and K. Biochar’s liming components that are inherited from feedstock and predictably modified by pyrolysis temperature provide a basis for optimizing the production of biochar with desired pH and liming characteristics.


Introduction
Chars are ubiquitous geosorbents that are created from incomplete pyrolysis of carbon-rich biomass in the presence of little or no oxygen. Chars applied as a soil amendment have been labeled "biochars" [1]: biochars sequester carbon in soils and their carbon is generally much more stable than carbon in the original biomass. Biochars produced at moderate to high temperatures are very resistant to microbial decomposition, and therefore can be retained in soils for mean residence times of hundreds to thousands of years [2]. In contrast, composts generally mineralize rapidly and therefore frequent application is generally needed to maintain sufficient nutrient content, cation exchange capacity, and C content [3]. In comparison to other forms of organic matter used as soil amendments, biochars also reduce waste volume, and reduces risks from pathogens, heavy metals [4], and organic pollutants [5]. Biochars have commonly been produced from wood wastes, grasses, crop residues, and animal manure feedstocks [6].
Physical and chemical properties of biochars vary greatly, potentially allowing one to produce a suitable biochar for specific applications by optimizing feedstock selection and pyrolysis conditions [6,7]. During pyrolysis, the electron cloud arrangement in the carbon skeleton varies as a function of the quantity of thermally modified carbon layers [8], resulting in significant differences in biochar properties. Biochar structure contains mostly amorphous and some crystalline parts (conjugated aromatic compounds) [9]; the latter increase with the highest treatment temperature (HTT) [10]. The biochar matrix also includes aromatic-aliphatic groups, residual volatiles, and inorganic ashes [11] with varying morphologies of voids, cracks, and pores (at micro-, meso-, and macro-scales) [12]. Pyrolysis temperature strongly affects the molecular architecture of biochar in terms of surface area, pore number, and pore size [13], as well as surface area functional groups [14]. Prior studies have suggested a linear response of biochar pH to HTT, but well-resolved temperature response curves are few (for example, see References [15,16]). In all cases, higher aromaticity is expected at higher pyrolysis temperatures [17], but feedstocks may vary in temperature response pattern. Particle size can influence biochar properties since larger particles resist the escape of primary products during pyrolysis compared to smaller biomass particles [18]. In addition, biochar yield commonly increases with increasing particle size [19].
Biochars generally have a substantial liming capacity, due in large part to mineral content carried over from feedstock [1], which makes them a possible substitute for liming agents such as dolomitic limestone or aglime. The addition of biochar (with a typical pH of [8][9] to acidic soils generally results in increased soil pH and a concomitant decrease in the mobility of cationic metal ions. This occurs due to the competition of H + ions and metal ions for cation exchange sites on the biochar surface and the liming effect on the soil matrix [20]; as with geological liming materials, the main chemical mechanism for liming involves the formation of carbonates [21]. Positive effects of biochar on crop growth and yield are most pronounced on acid soils [22], suggesting the importance of liming effects. Liming is also of importance in biochar mitigation of heavy metal bioavailability: the application of biochar to alkaline or neutral soils does not result in a significant reduction in metal mobility, suggesting that the reduction in metal mobility is not due to biochar sorption, but rather to liming effects [20,23,24]. Biochar's pH increases as acidic surface functional groups are lost during pyrolysis and replaced by neutral or basic fused aromatic moieties [25]. Alkaline elements such as Mg, Ca, and K may also increase the pH in biochars if they are abundant in the feedstock [26]. In addition to effects on metal mobility, direct nutrient release (PO 4 3− , NH 4 + , Mg 2+ , K + , and Ca 2+ ) from biochar is also strongly pH dependent [27]. The relative effects of biochar types and dosages on plant nutrient supply and responses have received considerable recent attention in the contexts of agriculture [7,28] and forestry [29,30]. A common theme in this body of research is high variability among biochars: some biochars strongly promote crop yields, some are neutral, and some are detrimental [6,7]. This variation provides a strong motivation for studies aimed at enhancing mechanistic understanding needed to produce optimized biochars for specific crop and soil systems [7]. Other critical soil properties-for example, nutrient and water retention and cation exchange capacity (CEC) [31]-are also strongly affected by biochar's porous structure and surface functional groups [32].
Wood chemistry is known to vary substantially among tree species, in terms both macromolecular constituents (cellulose, hemicellulose, and lignin) [33,34] and trace element concentrations [35]. However, surprisingly few data are available on variation in biochar properties among tree species, with few comparative studies in the literature presenting data from more than one wood species [14,29,36,37]. Hemicellulose and cellulose degrade at lower temperatures and earlier during pyrolysis than lignin [38,39]. Among common feedstocks, wood biomass contains the highest lignin content, though this varies substantially among tree species. Most conifers have higher lignin levels than do temperate hardwoods; however, data compilations indicate that hardwoods vary from 14%-48% lignin and conifers from 21%-38% lignin [33]. Lignin not only varies among tree species in terms of content, but also in terms of chemistry: angiosperm wood contains guaiacyl-syringyl linkages, while conifer wood mainly contains guaiacyl structures with a lower proportion of sinapyl units [40]. Variation in the thermal stability of oxygen-containing functional groups of these structures strongly modifies lignin's thermal decomposition range [41]. Conifer lignins are more thermally stable, while angiosperm lignins are less stable and may provide lower char yields [39]. Variation in biochar properties among wood species is thus expected to be strongly influenced by variation in lignin content and chemistry. This variation is also likely to result in large differences in responses of char properties to pyrolysis temperature.
Canadian forestry provides a wide range of wood feedstocks potentially applicable to biochar production [42]. Boreal and temperate forests together cover 417 million hectares and form the basis for Canadian forest products industries [43]. Sawmill byproducts-in particular, sawdust, lower-value wood chips, and bark residues-are important potential feedstock sources. Heschel and Klose [44] refer to these homogenous forest byproducts as "high-quality feedstocks". Lignin sludge from paper mills is another potential candidate feedstock [26]. By converting low-value or waste materials into biochar, the wood products industry could potentially both reduce environmental impacts and enhance product output and diversity. Quantification of variation in biochar physical and chemical properties is essential in assessing feedstock potential and in developing "designer" biochars for specific applications.
In the present study, we compare biochars derived from a broad range of common woods in Canada as potential biochar feedstocks, across a range of pyrolysis temperatures, focusing on physio-chemical parameters of importance to liming effects. We specifically address the following questions: (1) Is there a relationship between wood feedstock pH and the pH of the biochar produced? (2) How do concentrations of mineral nutrients and metals vary among biochars as a function of wood feedstock and pyrolysis temperature? (3) How does biochar pH increase with pyrolysis temperature, and does this relationship vary among wood feedstock species? (4) How do the concentrations of nutrient elements and acidic surface functional groups in char influence pH and liming capacity? and Pinus resinosa Sol ex. Aiton and Pinus strobus L. were supplied by Freymond Lumber, Bancroft, Ontario, Canada. All samples were sawdust with particle sizes of 63-4000 µm and were pyrolyzed with a laboratory-scale GSL-1100X programmable tube furnace (2.0" outer diameter, 1.85" inner diameter and 24" length quartz tube, and 3.5" length perforated stainless-steel sample container) after their moisture content was adjusted to 5%. All biochars were produced under a N 2 atmosphere and the final temperatures were held for 1 h. Charring temperatures were raised at rates of 10, 13, 17, 20, and 23 • C/min to 300, 400, 500, 600, and 700 • C, respectively. Following pyrolysis, samples were allowed to cool to below 100 • C before exposure to oxygen. Supplemental Table S1 summarizes the feedstock and biochar samples used. Both feedstocks and biochars were placed in airtight containers at 5 • C prior to analysis.

pH Measurement of Feedstock Biochar and Soil-Biochar Incubation
pH measurements were taken using an Oakton pH Meter with glass electrode. Prior to pH measurements, biochars were ground using a mortar and pestle. Biochar and feedstock samples were shaken in distilled water for 24 h in 0.15 g/6 mL and 0.8 g/10 mL quantities, respectively (feedstocks generally showed higher water sorption than biochar samples; the concentrations used reflected the minimum quantity of water necessary to obtain measurable liquid samples). All samples were then filtered (0.2 µm) and measured immediately.
The soil used in the incubation experiment was of the brunisolic order, with sandy texture from parent material of granite bedrock and a pre-treatment pH value of 5.3 (pH measured in 1 g/4 mL reverse osmosis water). This soil was collected from a shallow depth (<10 cm) in the Haliburton Forest and Wildlife Reserve (details on soil analysis are given in Sackett et al. [29]). Soil incubations were conducted with duplicate samples of biochar-soil mixtures (0.008 w/w ratio; 50 g soil) incubated at room temperature (~20 • C) for 18 days in 100 mL plastic cups. For soil incubation tests, the pH of the biochar-soil mixture was directly measured after stirring to homogenize the mixture. Temperatures were recorded at 1-h intervals throughout the incubation experiment (using a Logtag TRIX-8 temperature logger). Initially, 30 mL of deionized water was added to each soil-biochar mixture. Evaporated water was replaced every 3 days and before each pH measurement was taken.

Determination of Oxygen-Containing Surface Functional Groups by Boehm Titration
A modified Boehm titration was conducted following Laird et al. [45]. Briefly, prior to titration, all biochar samples were ground with a mortar and pestle. To remove acid-soluble ash components (carbonates, hydroxide, and oxides), ground biochars were treated with 1:50 biochar:solution (w/v) ratios. First, 0.05 M HCl was added to the biochars, which were then shaken for 24 h. Subsequently, the mixture was filtered and washed twice with 1 M CaCl 2 and four times with deionized water. The pre-treated samples were allowed to dry for 72 h at 50 • C. Dried samples (0.5 g) were shaken with 50 mL of 0.05 M Boehm reactants of NaOH, NaHCO 3 , and Na 2 CO 3 for 24 h. After a 24-h equilibration time, the mixtures were filtered and extracts were further treated with a sparge method [46]. Three 10-mL aliquots were taken from equilibrated extracts, and 0.05 M HCl was added at a 2:1 v/v ratio for both NaOH and NaHCO 3 and at a 3:1 v/v ratio for Na 2 CO 3 extracts. Acidified samples were sparged with N 2 for approximately 2 h to remove or minimize CO 2 prior to titration. NaOH (0.05 M) was standardized with potassium hydrogen phthalate and used to back-titrate samples where blank values were subtracted. End points were determined using phenolphthalein as an indicator.

Spectroscopy
Spectroscopic analyses were used to corroborate results of Boehm titration for a subset of chars analyzed. Solid-state nuclear magnetic resonance (SSNMR) spectroscopy analyses were performed using a Bruker 700 MHz spectrometer. Both 1D 1 H FastMAS spectra and 1D { 1 H} 13 C cross-polarization (CP) MAS spectra were acquired on an Agilent DD2 700 MHz spectrometer with an Agilent 1.6-mm FastMAS NMR probe. Only Pinus resinosa wood samples pyrolyzed at 300, 500, and 700 • C were analyzed by SSNMR. Each sample was ground into a finer powder using a mortar and pestle. The samples were packed into 1.6 mm zirconia rotors for 1D 1 H and 13 C NMR acquisition and processing.
Fourier-transform infrared (FTIR) spectroscopy was conducted using a Bruker spectrometer model Tensor 27 with a spectral resolution of 4 cm −1 within the range of 4000 to 400 cm −1 . Prior to spectra acquisition, all samples were prepared into compressed pellets from a homogenized mixture of 0.2 g Potassium bromide (KBr) and 0.25% biochar ground using a mortar and pestle.
Inductively coupled plasma-atomic emission spectroscopy (ICP-AES) analyses were performed using a PerkinElmer Optima 7300 DV with axially viewed plasma. Samples were digested with 5 mL trace metal grade (69%-70%) nitric acid following the published wet digestion method of Enders and Lehmann, 2012 [47]. Then, 0.5 g of biochar was transferred to a digestion vessel and allowed to sit for 24 h prior to digestion. Sample tubes were covered with glass stoppers both prior to and during digestion. The digestion block was then pre-heated to 120 • C and sample tubes were heated for 4 h. After cooling, 65 mL of deionized water was added to each sample to reach 5% nitric acid concentration and filtered (0.45 µm) prior to ICP analysis.

Statistical Analysis
Species and temperature effects on char properties were initially tested by two-way analysis of variance (ANOVA), where individual determinations of properties were treated as replicates. In cases where only a single measurement of char properties was made, the analyses were a non-replicated design, and thus required the assumption of no species x temperature interaction [48]. Differences in char properties between conifer and angiosperm species were also tested by one-and two-way ANOVA. Multiple linear regression was used to examine effects of feedstock properties on char properties. Non-linear least squares regression was used to fit alternative functions to relationships between pyrolysis temperature and char properties within a given species, using the nls() function in R. Where significant effects were found, the most parsimonious model was selected using an information theoretic approach [49] by minimum Akaike's information criterion (AICc, corrected for variable number). All statistical analyses were conducted using R version 3.1.2 [50].

Effects of Pyrolysis Temperature and Species on Biochar pH
Biochar pH generally increased with increasing pyrolysis temperature across all species (Figure 1a). There were statistically significant effects of both temperature and species on biochar pH (F = 354.4 1,57 and F = 6.5 1, 57 , respectively; p < 0.001 in both cases). The relationship between temperature and biochar pH was approximately linear from 300-700 • C (a second-order polynomial term in the model was not significant). In addition, there were statistically significant differences among species in feedstock pH (F = 6.9 1, 93 ; p = 0.01). Significant differences were also observed between the mean pH values for conifers and angiosperms, with higher values for both feedstocks and biochars found in angiosperms (feedstock pH: Conifers, 4.80 and angiosperms, 5.47; biochar pH: Conifers, 7.87 and angiosperms, 8.29 (Figure 1b)).

5
Species and temperature effects on char properties were initially tested by two-way analysis of variance (ANOVA), where individual determinations of properties were treated as replicates. In cases where only a single measurement of char properties was made, the analyses were a nonreplicated design, and thus required the assumption of no species x temperature interaction [48]. Differences in char properties between conifer and angiosperm species were also tested by one-and two-way ANOVA. Multiple linear regression was used to examine effects of feedstock properties on char properties. Non-linear least squares regression was used to fit alternative functions to relationships between pyrolysis temperature and char properties within a given species, using the nls() function in R. Where significant effects were found, the most parsimonious model was selected using an information theoretic approach [49] by minimum Akaike's information criterion (AICc, corrected for variable number). All statistical analyses were conducted using R version 3.1.2 [50].

Effects of Pyrolysis Temperature and Species on Biochar pH
Biochar pH generally increased with increasing pyrolysis temperature across all species ( Figure  1a). There were statistically significant effects of both temperature and species on biochar pH (F = 354.4 1, 57 and F = 6.5 1, 57, respectively; p < 0.001 in both cases). The relationship between temperature and biochar pH was approximately linear from 300-700 °C (a second-order polynomial term in the model was not significant). In addition, there were statistically significant differences among species in feedstock pH (F = 6.9 1, 93; p = 0.01). Significant differences were also observed between the mean pH values for conifers and angiosperms, with higher values for both feedstocks and biochars found in angiosperms (feedstock pH: Conifers, 4.80 and angiosperms, 5.47; biochar pH: Conifers, 7.87 and angiosperms, 8.29 (Figure 1b)).  ; biochar pH as a function of feedstock pH (c) and liming potential of biochar across five pyrolysis temperatures grouped by feedstock type; regression line drawn relating measured pH (pH of soil+biochar on incubation Day 18 = 0.15397*Biochar pH + 6.16463) (d). All analyses are from 19 temperate tree species (conifer species (n = 9) and angiosperm species (n = 10)).

Liming Potential of Biochar
Over the incubation period, soil-biochar mixtures generally showed an increase in pH values, and at day 18 and beyond, values plateaued. A strong correlation was found between biochar pH and the pH of incubated soil-biochar mixtures (at day 18: r = 0.600, p < 0.001). The relationship between biochar pH and incubated sample pH is shown in Figure 1d. There were statistically significant effects of both temperature and species on incubated soil-biochar sample pH (F = 300.0 1, 51 and F = 9.1 18, 51 , respectively; p < 0.001 in both cases). The relationship between biochar pH and incubated soil pH differed between feedstock types, i.e., conifer vs. angiosperm (F = 5.5 2, 92 ; p < 0.001).

Elemental Analysis
As pyrolysis temperature increased, biochar samples showed an increase in concentration of several alkaline elements (Ca and Mg-occurring at the highest concentrations-as well as K, Fe, Mn, Na, and Zn); in all cases, feedstocks contained the lowest concentrations of each of these elements (Figure 2a-g). ANOVA results indicated that there were significant effects of both species and pyrolysis temperature on alkaline element and metal concentrations, with the following exceptions: No relationship was found between pyrolysis temperature and either Na or Zn concentration (Table 1). Higher mean element concentrations were generally observed for biochars of angiosperms as compared to conifers; this relationship was significant for K and Mg (p < 0.001), Na (p < 0.01), and Fe (p < 0.05) ( Table 1). 6 measured pH (pH of soil+biochar on incubation Day 18 = 0.15397*Biochar pH + 6.16463) (d). All analyses are from 19 temperate tree species (conifer species (n = 9) and angiosperm species (n = 10)).

Liming Potential of Biochar
Over the incubation period, soil-biochar mixtures generally showed an increase in pH values, and at day 18 and beyond, values plateaued. A strong correlation was found between biochar pH and the pH of incubated soil-biochar mixtures (at day 18: r = 0.600, p < 0.001). The relationship between biochar pH and incubated sample pH is shown in Figure 1d. There were statistically significant effects of both temperature and species on incubated soil-biochar sample pH (F = 300.0 1, 51 and F = 9.1 18, 51, respectively; p < 0.001 in both cases). The relationship between biochar pH and incubated soil pH differed between feedstock types, i.e., conifer vs. angiosperm (F = 5.5 2, 92; p < 0.001).

Elemental Analysis
As pyrolysis temperature increased, biochar samples showed an increase in concentration of several alkaline elements (Ca and Mg-occurring at the highest concentrations-as well as K, Fe, Mn, Na, and Zn); in all cases, feedstocks contained the lowest concentrations of each of these elements (Figure 2a-g). ANOVA results indicated that there were significant effects of both species and pyrolysis temperature on alkaline element and metal concentrations, with the following exceptions: No relationship was found between pyrolysis temperature and either Na or Zn concentration (Table  1). Higher mean element concentrations were generally observed for biochars of angiosperms as compared to conifers; this relationship was significant for K and Mg (p < 0.001), Na (p < 0.01), and Fe (p < 0.05) ( Table 1).

Surface Group Analysis
Boehm titration was used to quantify the biochar samples' acidic surface functional groups (carboxylic, lactonic, and phenolic; Figure 3a-c). In all cases, there was a decrease in carboxylic and lactonic functional groups, and an increase in phenolic groups, with increasing pyrolysis temperature. ANOVA results showed significant effects of pyrolysis temperature on carboxylic, lactonic, and phenolic functional group concentrations (F = 150.0 1, 90 , F = 56.4 1,87 , and F = 48.0 1, 86 , respectively; p < 0.001), and showed a significant effect of species on lactonic group concentration only (F = 2.4 18,70 ; p < 0.01) ( Table 1). There was also a consistent difference in the concentrations of surface functional groups between conifers and angiosperms: Angiosperm-derived biochars had overall higher concentrations of carboxylic and lactonic groups, but similar concentrations of phenolic groups (Figure 3a-c).

FTIR and NMR Spectroscopy
Results obtained from FTIR and NMR spectra further demonstrated the decline of acidic groups and increase of aromatic groups with increasing pyrolysis temperature (Figure 4a). In the FTIR

FTIR and NMR Spectroscopy
Results obtained from FTIR and NMR spectra further demonstrated the decline of acidic groups and increase of aromatic groups with increasing pyrolysis temperature (Figure 4a). In the FTIR spectrum, biochars produced from 300 to 600 • C exhibited a decrease in relative intensity of the following bands as the temperature increased: 3400 cm −1 (O-H stretching of hydroxyl groups), 2900 and 2850 cm −1 (C-H asymmetric and symmetric stretching of aliphatic groups), 1620 cm −1 (C=O stretching of carboxyl mode), and 1026 cm −1 (C-O symmetric stretching of cellulose, hemicellulose, and lignin). Over this same range of pyrolysis temperatures, other bands demonstrated an increase in relative intensity: Specifically, 831 and 700 cm −1 (C-H aromatic deformation modes), 1600 cm −1 (C=C aromatic stretching and C=O stretching of conjugated ketones and quinones), and 1400 cm −1 (C-C stretching of aromatic rings). Both 1 H and 13 C NMR spectra (Figure 4b,c) detected a decrease in the aliphatic carbons and increase in the aryl carbons. In the 1 H NMR spectrum of biochar prepared at 300 °C, the methylene chemical shift at δ = 3.39 ppm decreased in intensity at 500 °C and disappeared at 700 °C. Aromatic chemical shifts at δ = ~7 ppm showed increased intensity at 500 °C and increased broadness at 700 °C. Negative chemical shifts indicate the attachment of inorganic moieties and the high-intensity peak at 700 °C suggests the presence of carbonates at high concentration. In the 13 C NMR spectrum at 300 Both 1 H and 13 C NMR spectra (Figure 4b,c) detected a decrease in the aliphatic carbons and increase in the aryl carbons. In the 1 H NMR spectrum of biochar prepared at 300 • C, the methylene chemical shift at δ = 3.39 ppm decreased in intensity at 500 • C and disappeared at 700 • C. Aromatic chemical shifts at δ =~7 ppm showed increased intensity at 500 • C and increased broadness at 700 • C. Negative chemical shifts indicate the attachment of inorganic moieties and the high-intensity peak at 700 • C suggests the presence of carbonates at high concentration. In the 13 C NMR spectrum at 300 • C, numerous carbon signals which are related to the feedstock content (cellulose, lignin) and additional chemicals inherited from the feedstock were present. However, at 500 and 700 • C, only aromatic carbons produced signals.

Modeling the Impact of Element and Surface Group Concentrations on Biochar pH
The acidic surface functional group concentrations of carboxylic, lactonic, and phenolic acid groups (measured in mmol/g), as well as the elemental nutrient concentrations of base metals Ca, Mg, and K (measured in mg/kg), were expected to have an impact on biochar pH. We examined alternative functional forms for this relationship, using a non-linear least-squares method to fit models, and comparing the fit of alternative models based on an information theoretic approach [49] (Table S3). Based on the minimum AIC, the following exponential form best predicted biochar pH as a function of the ratio of carboxylic to phenolic functional groups and concentrations of Ca and K (Equation (1)): where φ and α are 6.393 and 3.075; ctpr is the carboxylic to phenolic ratio; Ca and K are the concentrations of Ca and K in mg/kg; and the fitted coefficients β, γ, and δ are −0.5284, 1.977 × 10 −5 , and 1.496 × 10 −5 , respectively. In Figure 5, biochar pH is plotted as a function of the ratio of carboxylic to phenolic groups, with the size of points proportional to summed concentrations of Ca and K. The residual plots for the fitted model are presented in Supplemental Figure S1. additional chemicals inherited from the feedstock were present. However, at 500 and 700 °C, only aromatic carbons produced signals.

Modeling the Impact of Element and Surface Group Concentrations on Biochar pH
The acidic surface functional group concentrations of carboxylic, lactonic, and phenolic acid groups (measured in mmol/g), as well as the elemental nutrient concentrations of base metals Ca, Mg, and K (measured in mg/kg), were expected to have an impact on biochar pH. We examined alternative functional forms for this relationship, using a non-linear least-squares method to fit models, and comparing the fit of alternative models based on an information theoretic approach [49] ( Table 2). Based on the minimum AIC, the following exponential form best predicted biochar pH as a function of the ratio of carboxylic to phenolic functional groups and concentrations of Ca and K (Equation (1)): where φ and α are 6.393 and 3.075; ctpr is the carboxylic to phenolic ratio; Ca and K are the concentrations of Ca and K in mg/kg; and the fitted coefficients β, γ, and δ are −0.5284, 1.977 × 10 −5 , and 1.496 × 10 −5 , respectively. In Figure 5, biochar pH is plotted as a function of the ratio of carboxylic to phenolic groups, with the size of points proportional to summed concentrations of Ca and K. The residual plots for the fitted model are presented in Supplemental Figure S1.

Discussion
Our results indicate that pH of wood-derived biochars is strongly influenced by both pyrolysis temperature and variation in feedstock wood chemistry. Results are also consistent with the predictions that species differences in the temperature dependence of biochar properties are mainly determined by variation in wood chemical composition. The pH of biochars produced at higher pyrolysis temperatures (500-700 °C) is strongly correlated with feedstock pH, while less correlation

Discussion
Our results indicate that pH of wood-derived biochars is strongly influenced by both pyrolysis temperature and variation in feedstock wood chemistry. Results are also consistent with the predictions that species differences in the temperature dependence of biochar properties are mainly determined by variation in wood chemical composition. The pH of biochars produced at higher pyrolysis temperatures (500-700 • C) is strongly correlated with feedstock pH, while less correlation is observed between the pH of biochars produced at lower temperatures (300-400 • C) and feedstock pH. Correlations between wood feedstock pH and biochar pH indicate that biochar inherits important chemical characteristics from feedstock, and that biochar pH can in part be predicted from feedstock pH. Thermochemical lignocellulose degradation begins at approximately 130 • C, where lignin starts softening; hemicellulose degrades at 150-350 • C, cellulose at 275-350 • C, and lignin at 250-500 • C [18]. The chemical products of hemicellulose degradation are predominantly "wood vinegars" that include a mixture of acetic, propionic, and butyric acids, in addition to other low molecular weight organic acids [51]. Chars can sorb and retain these molecules [51,52]; this is likely to reduce biochar pH in a manner dependent on hemicellulose content and sorption capacity of the char and its non-pyrolyzed constituents. We speculate that the production and sorption of these strongly acidic products may account for the lack of consistent relationships between feedstock and char pH at low pyrolysis temperatures. Possibly, the respective fractions of the wood constituent's cellulose, hemicellulose, and lignin would show a relationship between feedstock and biochar pH at low pyrolysis temperatures. No prior study has examined relationships between feedstock pH and biochar pH, and, thus, there are no available data to compare to the present results. The overall lower pH of conifer woods and biochars may be influenced by the higher concentration of extractives in conifers, averaging 21% compared to 15%-16% in angiosperms [33,34].
The increase in biochar pH with pyrolysis temperature is accompanied by two parallel patterns: (1) An increasing concentration of basic elements (alkali metals and alkaline earth metals) that are inherited from the feedstock; and (2) a decrease in the concentration of acidic surface functional groups. These results are consistent with prior studies that have found an increase in elements important as plant nutrients [26] and a decrease in acidic functional groups [53] with increasing pyrolysis temperatures. These chemical changes jointly determine biochar pH, which in turn may be used to predict liming capacity. We derive an empirical non-linear expression (Equation (1)) that predicts biochar pH as a function of the ratio of carboxylic to phenolic functional groups and the concentrations of Ca and K in the char.
Elemental analyses indicated that sampled angiosperms had, on average, higher concentrations of Ca, K, Mg, Mn, Fe, Na, and Zn in both feedstocks and in biochars than did conifers (Table 1). This result is consistent with prior surveys of wood chemistry [35], but the generalization of this pattern to chars does not appear to have been previously shown. Our results indicate that biochar pH is most strongly correlated with Ca and K concentration. Studies by Novak et al. [54,55] previously examined biochar's chemical and physical characteristics, in particular pH, in relation to feedstock and pyrolysis temperature. The Novak et al. [55] study analyzed five different lignocellulosic feedstocks and a poultry litter, finding that as the pyrolysis temperature increased, the biochar's pH, ash content, and aromaticity increased. This is in agreement with our observations obtained through Boehm titration, FTIR, and NMR spectroscopy. The pine chip char produced by Novak et al. [55] had a lower pH of 6.7 at 500 • C, compared to 7.96 for both Pinus resinosa and Pinus strobus. Longer pyrolysis residence time increases the concentration of acidic functional groups [56], thus contributing to a decrease in the char's pH. Other factors such as particle size, N 2 flow, pyrolizer design, pre-and post-preparation procedures, and conditions of feedstock may also have contributed to the lower pH observed by Novak et al. [55], in addition to effects of condensed volatile matter during pyrolysis. Novak et al. [55] found no relationship between pyrolysis temperature and macronutrient concentrations in wood-derived chars. By contrast, here we observed an increase in mean nutrient concentration with pyrolysis temperature increase across all species and all elements examined, with the exceptions of Na and Zn (Figure 2f,g). This pattern seems likely to be general at pyrolysis temperatures ranging from 300-700 • C, since higher temperatures dictate a progressive loss of oxygen-containing groups with little volatilization of the nutrient elements in question [57].
Results obtained from the Boehm titration revealed a decrease in carboxylic and lactonic functional groups and an increase in phenolic groups with increasing pyrolysis temperature. This general result is consistent with earlier studies [25,56,58]. Analyses of FTIR and NMR spectra confirmed the results obtained from Boehm titration, showing that aromatic carbons associated with the feedstock biomass increased in concentration with increasing pyrolysis temperature, while the concentration of surface acidic functional groups decreased. The latter result is likewise consistent with previous studies [59,60]. Our observations are novel in that we also observed pronounced variation among wood species in the temperature-dependent changes in oxygen-containing functional groups. Higher concentrations of carboxylic and lactonic groups at a given pyrolysis temperature were generally observed among angiosperms compared to conifers. Biochar pH was negatively correlated with the concentration of carboxylic functional groups, but only weakly to the concentration of lactonic groups, and positively correlated with the concentration of phenolic groups.
Although it has commonly been assumed that more basic biochars will show higher liming capacity, this relationship has not previously been demonstrated. In the present study, all the biochars examined acted to increase soil pH from 5.3 to a range of 7.2-8.0 at a low dosage (0.8% w/w). Generally, the highest incubated soil pH was observed at the highest biochar pyrolysis temperatures. The biochar-soil mixture pH increased until it reached equilibrium on Day 18 (Table S2), during which time various chemical reactions could contribute to bringing the biochar pH to equilibrium. The flow chemistry of soluble elements found on biochars varies depending on their chemical composition, temperature, and over time [61], and the leaching mechanisms and amount depend on the biomass type and pyrolysis temperature [62]. For example, at room temperature, water-soluble divalent species of Ca 2+ and Mg 2+ have a smaller diffusion coefficient (0.792 × 10 −9 and 0.706 × 10 −9 m 2 /s, respectively) as compared to the monovalent K + and Na + (1.957 × 10 −9 and 1.334 × 10 −9 m 2 /s, respectively), because the divalent species are bonded more tightly to the biochar, resulting in slower leaching times of divalent ionic species [62,63]. Biochars produced at relatively high temperatures convert their nutrient elements into oxide and carbonate forms, both of which increase biochar alkalinity [64,65]. The water present in the soil-biochar mixture during incubation likely plays an important role in pH equilibrium by reacting with carbonates to form free ions (i.e., Ca 2+ ) [66], and this dissolution of carbonates (e.g., CaCO 3 ) results in an observed soil pH of between 7.2 and 8.5 [67]. Variable results have been obtained by prior soil incubation studies. Basso [68] observed a pH increase with soil incubation time; however, Novak et al. [69] did not. Basso suggested that the difference in their results occurred due to the ash content difference in biochars: In Basso's case, an increase in the soil pH was mainly attributed to the dissolution of carbonates, oxides, and hydroxides of the biochar ash fraction. Obia et al. [70] also found an increase in pH during incubation, specifically, a 2.3-unit increase upon the application of 1% cacao shell biochar to a sandy loam tropical acrisol soil. Obia et al. [70] suggested that carbonates are a key contributor to the alkalinisation effect.
In addition to the effects of basic elements, our results also showed consistently higher pH in both biochar and incubated soil-biochar as carboxylic and lactonic groups decline and phenolic groups increase. Phenolic groups have an aromatic character and are less acidic than carboxylic or lactonic groups. Thus, the increase in concentration of phenolic groups as pyrolysis temperature increases results in higher biochar pH and aromaticity [25]. In conventional batch-leaching systems (in the case of this experiment, soil-biochar cups), acidic species leachate contact with biomass leads to the leaching of some water-insoluble organically bound inorganic species (i.e., Cl and Si) [70]. The inorganic species leaching in a batch condition typically show two-step kinetics: an initial short-period rapid-leaching step followed by a long-period slow-leaching step [62]. Previous studies have noted that in pure water, biochars could also leach some organic matter (i.e., 2-3 and 4-5 aromatic rings) for a period of up to one month [62,[71][72][73], suggesting that organic matter leaching may also influence incubated soil-biochar pH equilibrium.
In terms of applications, the use of biochar as a soil liming agent will have a "double" effect on climate mitigation by: (1) acting as a long-term C storage in soil and decreasing greenhouse gas emissions [7,74]; and (2) reducing CO 2 emissions from geological liming agents. Application of biochar instead of other liming agents is also beneficial due to the longevity of chars in soil [2], making it possible to omit the frequent application of lime usually required to sustain nutrient availability and exchange capacity [75]. Biochar effects on soil pH in situ may also be affected by long-term weathering processes. The natural oxidation of biochar is believed to enhance phenolic and carboxylic functional groups, thus improving soil ion exchange capacity [25], as suggested by the long-lived fertility of Terra Preta soils [76].
In conclusion, our findings suggest that differences in wood chemistry are reflected in biochar's liming capacity, which is a central concern for soil amendment applications. Many prior biochar studies do not specify the actual species of wood used, and/or list only broad common names (e.g., "pine"). The differences we document among tree species have important implications for reporting in future biochar studies. Variation in pH and liming capacity of biochars is related to both feedstock and pyrolysis temperature. This work indicates that species differences in wood pH correspond to significant differences in biochar pH and provides a basis for predicting biochar pH from biochar chemistry, and from feedstock pH and pyrolysis temperature. The elevation of pyrolysis temperature corresponds to an increase in elemental nutrient concentrations and a decrease in acidic surface groups, ultimately resulting in an increase in biochar pH. Predicting biochar pH prior to production is of crucial importance because it makes possible the production of application-specific biochars for targeted purposes.