Feedstock choice, pyrolysis temperature and type influence biochar characteristics: a comprehensive meta-data analysis review

Various studies have established that feedstock choice, pyrolysis temperature, and pyrolysis type influence final biochar physicochemical characteristics. However, overarching analyses of pre-biochar creation choices and correlations to biochar characteristics are severely lacking. Thus, the objective of this work was to help researchers, biochar-stakeholders, and practitioners make more well-informed choices in terms of how these three major parameters influence the final biochar product. Utilizing approximately 5400 peer-reviewed journal articles and over 50,800 individual data points, herein we elucidate the selections that influence final biochar physical and chemical properties, total nutrient content, and perhaps more importantly tools one can use to predict biochar’s nutrient availability. Based on the large dataset collected, it appears that pyrolysis type (fast or slow) plays a minor role in biochar physico- (inorganic) chemical characteristics; few differences were evident between production styles. Pyrolysis temperature, however, affects biochar’s longevity, with pyrolysis temperatures > 500 °C generally leading to longer-term (i.e., > 1000 years) half-lives. Greater pyrolysis temperatures also led to biochars containing greater overall C and specific surface area (SSA), which could promote soil physico-chemical improvements. However, based on the collected data, it appears that feedstock selection has the largest influence on biochar properties. Specific surface area is greatest in wood-based biochars, which in combination with pyrolysis temperature could likely promote greater changes in soil physical characteristics over other feedstock-based biochars. Crop- and other grass-based biochars appear to have cation exchange capacities greater than other biochars, which in combination with pyrolysis temperature could potentially lead to longer-term changes in soil nutrient retention. The collected data also suggest that one can reasonably predict the availability of various biochar nutrients (e.g., N, P, K, Ca, Mg, Fe, and Cu) based on feedstock choice and total nutrient content. Results can be used to create designer biochars to help solve environmental issues and supply a variety of plant-available nutrients for crop growth.


Introduction
Biochars are carbon (C) rich materials typically produced via biomass pyrolysis at relatively low temperatures (300-700 °C) under limited oxygen conditions (Lehmann and Joseph 2009). Biomass, the feedstock for biochar creation, is typically derived from agricultural and forestry waste products, municipal waste, green and food waste. The creation of biochar from these products places C into a recalcitrant form which could last hundreds to over thousands of years (Spokas 2010;Kuzyakov et al. 2014;Wang et al. 2016), suggesting that biochar could aid in climate change mitigation (Tripathi et al. 2016) as one of the few negative greenhouse gas emission technologies with sustainable development co-benefits (Smith et al. 2019).
Over shorter time scales (e.g., one to several years), biochars have been proven to improve environmental quality by sorbing heavy metals and organic contaminants (e.g., Sigua et al. 2019;Cui et al. 2019;Novak et al. 2019a), positively affect soil water relations (e.g., Lentz et al. 2019;Kammann et al. 2011), reduce greenhouse gas emissions (e.g., Fuertes-Mendizábal et al. 2019;Borchard et al. 2018;Jeffery et al. 2016), and improve crop growth (e.g., Laird et al. 2017;Novak et al. 2016;Liu et al. 2013). Although creation of biochars for the above purposes may seem simply based on feedstock selection, biochar production for environmental improvements is complex. Feedstock selection, pyrolysis temperatures, and pyrolysis types can all greatly influence the final biochar product (Cao et al. 2017;Cha et al. 2016). Thus, increasing the understanding of the interaction between feedstock, pyrolysis temperature, and production technique (i.e., either fast or slow pyrolysis) would help biochar stakeholders make more well-informed choices for its use.
In terms of feedstock, understanding how initial feedstock properties influence final biochar characteristics is important. Feedstocks have been shown to play a major role in creating biochars with distinctly different chemical properties (Funke and Ziegler 2010;Novak et al. 2019b). In relative terms, wood-based biochars contain more C and lower plantavailable nutrients, manure-based biochars show opposite trends, and grass-based biochars typically fall somewhere in between woody and manure biochars (Ippolito et al. 2015). However, these properties may be altered by pyrolysis temperature and pyrolysis technique used for biochar creation.
Pyrolysis temperature and production technique play key roles in creating biochars with various chemical and structural properties. For example, nutrient availability changes drastically as pyrolysis temperature is increased (Clough et al. 2013;Nguyen et al. 2017). Specifically, with increasing pyrolysis temperature one typically observes increasing biochar C, P, K, Ca, ash content, pH, specific surface area (SSA), and decreasing N, H, and O content (e.g., Weber and Quicker 2018;Ippolito et al. 2015). These biochar characteristics are driven by forcing C into more condensed, recalcitrant forms, the creation of oxide/carbonate mineral phases (e.g., P, K, Ca) leading to greater ash content and higher pH, and the loss of N, H, and O via volatilization. Volatilization losses further concentrate other remaining elements (Kim et al. 2012;Kinney et al. 2012). When choosing pyrolysis type, in general slow pyrolysis tends to produce biochars with greater N, S, available P, Ca, Mg, surface area, and cation exchange capacity (CEC) as compared to fast pyrolysis.
The above aspects of biochar creation have led individual researchers to pay attention to biochar end-product properties. Current literature contains an untapped wealth of information regarding feedstock, pyrolysis type and temperature choices. However, there is some uncertainty in the literature with a vague description of how biochar characteristics are influenced by feedstock choice, pyrolysis type, and temperature employed. Thus, we chose to provide in this review a clearer picture by synthesizing existing literature to create a true comprehensive review of biochar properties based on feedstock choice, pyrolysis temperature, and pyrolysis type. Presenting this type of data is paramount for improving our understanding of the factors used to create biochar and characteristics within the final product. Establishing these comparisons should aid biochar practitioners and stakeholders in making well-informed decisions for biochar use as amendments in soils and for environmental mitigation purposes in mine spoils or metal contaminated soils. The objective of this work was focused on comprehensively reviewing how different feedstocks, pyrolysis temperatures, and production techniques affect biochar physicochemical properties, total and available nutrient content, other characteristics, and what these properties indicate when biochar is used for agricultural benefits.

Materials and methods
Data were compiled from literature that compared biochar macro-and micro-nutrient concentrations, pH, SSA, electrical conductivity (EC), cation and anion exchange capacity (CEC and AES, respectively), calcium carbonate equivalent (CCE), total pore volume (PV), average particle size (APS), and ash content under different biochar production temperatures and types. The literature examined was published between January 2009 through December 2016, searching electronic databases including Web of Science™ and Sco-pus™ for the keyword "biochar". Overall, 5394 publications presented reliable and valid physicochemical biochar properties which were fed into a database obtaining 50,851 individual observations (see Supporting Table 1 for references cited). We chose not to collect data from 2016 onwards because of the overwhelming points already collected from between 2009 and 2016.
Data from the selected peer-reviewed literature should provide the reader with overarching changes in biochar physicochemical characteristics, and thus where appropriate, specific inferences were made based on logical comparisons between biochar properties. The data are presented by feedstock choice, pyrolysis temperature, and pyrolysis type (i.e., either fast or slow pyrolysis). The data are further separated into nutrient availability based on feedstock choice, in order to provide biochar producers with guidance for preparing biochars to meet plant nutrient needs.

Feedstock choice
Wood feedstock data were derived from literature reporting chipped, shaved, bark, peeled, and other similar wood-based biochar products. All wood data were initially combined into one category because it has been previously shown that, in general, hardwood and softwood biochars have similar characteristics (Ippolito et al. 2015). However, this concept was revisited in the current manuscript due to the larger dataset involved as compared to Ippolito et al. (2015); wood feedstock data were further separated into either hardwood-(including bamboo) or softwood-derived biochars for specific nutrient availability analysis presented below.
Crop waste data included biochars made from a variety of agricultural crop residues, such as corn, wheat straw, rice straw and husk, potato, soybean, sugarcane and bagasse, cotton, grape, orange, peanut, and rape seed/straw. Based on the dataset collected, the major crop wastes utilized to produce biochars were corn, wheat straw, and rice straw/husk, and, therefore, this datum was further categorized for derivation of specific nutrient availability analyses presented below.
Other results from grasses were compiled from species that were typically, but not entirely, used for bioenergy, industrial production, or fodder. This dataset included grassy and leafy species such as Miscanthus (Miscanthus spp.), switchgrass (Panicatum virgatum), giant red (Arundo spp.), common reed (Phragmites australis), Typha (Typha spp.), water hyacinth (Eichhornia crassipes), and other minor grass species. Opposite of other feedstocks, this dataset was not further refined for nutrient availability analyses due to a lack of dataset robustness.
The manures/biosolids dataset included biochars made from papermill sludges, cattle and dairy manure, horse manure, swine manure, poultry manure and litter, and biosolids. Among all feedstocks, poultry manure, pig manure, cattle and dairy manure, and biosolids were the most common feedstocks employed to produce biochars. Thus, biochars made from these four feedstock choices were further separated for defining their specific nutrient availability as presented below.

Pyrolysis temperature and type
For ease of discussion, pyrolysis temperatures were grouped as follows: < 300, 300-399, 400-499, 500-599, 600-699, 700-799, and > 800 °C. A similar approach was utilized by Ippolito et al. (2015) when describing changes in biochar properties as affected by pyrolysis temperature. Pyrolysis type was simply separated into either fast or slow pyrolysis as indicated in the manuscripts evaluated.

Data interpretation
All individual data comparisons (e.g., total C, total O, total, H, etc.) between fast or slow pyrolysis were analyzed using a t test at α = 0.05. All individual data comparisons between various feedstocks or all pyrolysis temperatures were analyzed using analysis of variance (ANOVA). The Shapiro-Wilk and Levene's tests were also used to analyze for assumptions of normality and equal variance, respectively. In cases where assumptions were violated, data were transformed using either log 10 , natural log, or a square root function and re-analyzed for meeting assumptions of normality and equal variance. Tukey's Honestly Significant Difference pairwise comparisons were utilized to discern differences between various feedstocks or pyrolysis temperatures for individual (non-)transformed data at α = 0.05. The non-transformed data means are presented below.
For more specific data interpretations, we utilized the entire dataset, or specific feedstock subsets, along with SigmaPlot 13.0 for graphical interpretations. Within Sig-maPlot 13.0, we utilized regression fitting functions (linear, quadratic, exponential rise or decay, and exponential growth equations) that best fit (e.g., best R 2 values and greatest p values) the presented data. We avoided data interpretations that made no logical sense (e.g., sigmoidal equations, sinusoidal equations, etc.) yet may have had greater R 2 values than regression equations that made logical, interpretable sense. When necessary, Pearson's correlations were utilized to help drive discussion for biochar utilization from specific feedstock choices.

Pyrolysis type
In terms of the general utility in choosing pyrolysis type, slow pyrolysis uses a slow temperature heating rate (0.01-2 °C s −1 ) as compared to fast pyrolysis, and if adjusted properly can produce approximately equal ratios of solid (i.e., biochar), gas, and liquid products (Sohi et al. 2009). Fast pyrolysis uses higher heating rates (> 2 °C s −1 ) and shorter residence times (< 2 s) during the thermal conversion process, which provides greater bio-oil yields (75%) but lower gas and biochar quantities (Qambrani et al. 2017).
The influence of pyrolysis type on final biochar physicochemical properties is presented in Table 1. Fast pyrolysis led to biochars containing greater SSA and lower APS as compared to slow pyrolysis, both parameters likely a function of the relatively small initial feedstock particle size utilized during fast pyrolysis in fluidized bed reactor (Asadullah et al. 2010). Slow pyrolysis led to biochars containing greater CCE and ash content as compared to fast pyrolysis. Increasing ash content and CCE are related, as acidic functional groups are reduced (Novak et al. 2009) and mineral hydroxide and carbonate phases are increased with increasing ash content (Knicker 2007).

General feedstock choice
Various feedstock sources also influence final biochar properties (Table 1). Wood-based biochars typically contained a greater SSA and PV as compared to other feedstock choices. In terms of SSA and PV, this is a function of reducing relatively large wood-based cell structures to smaller pores and thus increasing overall SSA (Weber and Quicker 2018;Downie et al. 2009) and, consequently, PV. Increasing SSA and PV may also be associated with gas or water volatilization processes during pyrolysis and the loss of micro-molecule organic compounds, both of which can create voids within the biochar matrix during pyrolysis (Guo and Lua 1998). Contrarily, biochars produced from manures/ biosolids typically exhibited relatively low SSA, likely due to the development of deformation, structural cracking or micropore blockage (Ahmad et al. 2014;Lian et al. 2011), along with less distinct porous structures in the feedstock as compared to wood-based biochars.
Crop waste, other grasses, and manures/biosolids biochars had a greater CEC and pH compared to wood-based biochars. CEC can be generated during pyrolysis, as oxygenated surface and inorganic functional groups are formed (Briggs et al. 2012). The increased CEC may also be attributed to elevated pH leading to pH-dependent charge, or insoluble precipitates present in the ash that act as reaction sites . Others have related increased CEC to increases in SSA (Kloss et al. 2012;Qambrani et al. 2017;Liang et al. 2006), yet the data presented here contradict these findings. Concomitant with increasing ash content from wood, to crop, to other leafy-grassy, to manures/biosolids biochars, the CCE is also greatest in manure-biosolidsderived biochars. This is likely due to oxide and hydroxide mineral phase precipitation during pyrolysis (e.g., Ippolito et al. 2012).
Specific surface area increased with increasing pyrolysis temperature, as shown by others (Domingues et al. 2017;Luo et al. 2015). This is a function of shrinking the solid matrix, causing relatively large pores to become smaller and thus increasing overall SSA (Weber and Quicker 2018;de Mendonça et al. 2017). Surface area has been previously correlated with sorption/retention of nutrients and contaminants, while pore volume (which, in general increases with 1 3 pyrolysis temperature) is assumed to affect water availability and soil aeration (Ajayi and Horn 2016;Qambrani et al. 2017). Increasing pyrolysis temperature also clearly increased biochar ash content and pH, as solid phase hydroxide and carbonate phases increase within the ash, causing pH values to concomitantly increase. This concept is described in more detail in the biochar physicochemical correlations, Sect. 3.1.4 below.

Biochar physicochemical correlations
Correlations between biochar pH and ash content, pyrolysis temperature and pH, and pyrolysis temperature and SSA, based on all pertinent data collected, are shown in Fig. 1. Increasing biochar pH generally correlates with an increase in ash content ( Fig. 1a) as also shown by Lehmann (2007). This correlation could be useful for practitioners designing pyrolysis systems for creating biochars needed for soil liming purposes (e.g., in acidic or acid generating mines soils; Ippolito et al. 2017;Sigua et al. 2019;Godlewska et al. 2017). Increasing pyrolysis temperature caused biochar pH to increase (Fig. 1b), likely due to loss of acidic functional groups (Novak et al. 2009) and the formation of Ca-Mg-, Na-, and K-bearing oxide, hydroxide, and carbonate mineral phases (Cao and Harris 2010) that can raise the pH to ranges extending from 9.9 to 13 (aqion 2019). Ngatia et al. (2017) also illustrated biochar pH correlations to alkali earth elements such as those mentioned above. Increasing pyrolysis temperature also increased biochar specific surface area ( Fig. 1c), as observed by others (e.g., Lu et al. 2012;Hass et al. 2012). Increasing specific surface area appears to be a function of decreases in cell pore diameter (Ahmad et al. 2012), along with tar, oils, H and O removal with increasing pyrolysis temperatures (Kloss et al. 2012;Chen et al. 2008).

Pyrolysis type, general feedstock choice, and pyrolysis temperature
Changes in biochar total macro-elements as a function of pyrolysis type, general feedstock choice, and pyrolysis temperature, are shown in Table 2. Fast pyrolysis favored total S, K, Ca, and Mg content as compared to slow pyrolysis. However, in general, pyrolysis type had little influence on total macro-elements as compared to feedstock choice and temperature; total macro-elements have been previously correlated with feedstock choice and pyrolysis temperature (Zhang et al. 2017b). In terms of feedstock choice, woodbased feedstocks led to greater C content over other feedstock choices. If C storage was a goal, wood-based biochars, especially with high aromatic-C character  and low O/C and H/C ratios (Spokas 2010) should be applied to soils. Grass-based feedstocks contained relatively elevated K and Ca content as compared to other feedstocks, while manures/biosolids-based feedstocks contained the greatest N, S, P, Ca, and Mg concentrations over biochars from other feedstocks, similar to findings by Amoah-Antwi et al. (2020). At first glance, it appears that if N-P-K fertilizer use was a goal, grass-or manure/biosolids-based biochars could be utilized (Lin et al. 2017;Gondek and Mierzwa-Hersztek 2016;Wang et al. 2013;Brewer et al. 2012); this concept is revised within the context of specific feedstock choices in Sect. 6.5 below. In terms of pyrolysis temperature, increasing temperature increased C, P, Ca, and Mg concentrations, also as suggested by Chen et al. (2016). Specifics, with respect to biochar total elemental analysis, are also discussed below.

Carbon, hydrogen, and oxygen
For biochars, total C content is often elevated as most feedstocks contain appreciable C concentrations, yet feedstock choice significantly affects biochar C content (Table 2). Wood-based biochars contained greater C as compared to biochars made from other feedstock choices, simply due to a lack of other elements (e.g., N, S, P, K, Ca, and P) leading to a smaller C-dilution effect in wood-based biochars. Biochar C compounds can be grouped into relatively condensed (i.e., stable, non-mineralizable) aromatic C compounds, compared to easier-degradable, micro-molecular or water-soluble C (Lehmann et al. 2011). During pyrolysis, dehydration, cleavage, and polymerase reactions cause easily degradable C compounds to be restructured, while other elements may be lost to volatilization, and thus overall biochar total C content increases with increasing pyrolysis temperature (Weber and Quicker 2018;Antal and Grønli 2003; Table 2). Biochar C bioavailability is temperature dependent, with higher pyrolysis temperatures related to larger non-labile C fractions (Nelissen et al. 2012). Greater pyrolysis temperatures tend to create relatively stable, aromatic C compounds and silicate-carbon complexes that are usually regarded as recalcitrant to microbial oxidation (Guo and Chen 2014). When placed in soils, recalcitrant C could last hundreds to thousands of years, and as such, biochar land application may play a role in climate mitigation (Schmidt et al. 2019;Werner et al. 2018;Woolf et al. 2018;Bolan et al. 2012).
Increasing pyrolysis temperature significantly increased biochar C content via volatilization losses of other elements, especially H and O ( Table 2). As pyrolysis temperature increases, water, organic surface functional groups, and tars are lost, all of which contain H and O atoms (Ahmad et al. 2014;Cantrell and Martin 2012). During pyrolysis, O is released at a greater rate than H, with the final biochar product characterized by a decrease in H/C-ratio and containing low oxygen-content. Thus, understanding the ratios between H, C, and O, as depicted in the Van-Krevelen diagram, may be useful for discerning biochar longevity in soil systems (Schimmelpfennig and Glaser 2012;Spokas 2010).
Specifically, the atomic ratios of O/C to H/C have been used to describe the pyrolysis carbonization process (Weber and Quicker 2018), and more specifically to describe several key biochar environmental longevity factors as a function of these ratios. Figure 2 illustrates how the O/C and H/C ratios are altered as a function of pyrolysis temperature. Spokas (2010) described biochars with O/C ratios: (a) < 0.2 as highly stable (half-life > 1000 years); (b) between 0.2 and 0.6 as moderately stable (half-life between 100 and 1000 years); and (c) > 0.6 as relatively unstable (half-life less than 100 years). Meanwhile, biochars with H/C ratios < 0.7 have greater fused aromatic ring structures and have been thermochemically altered as compared to biochars with H/C ratios > 0.7 (IBI 2015). Biochar H/C ratio has been also found to be a key factor defining its N 2 O mitigation potential, with biochars with H/C ratios < 0.3 being the most effective (Cayuela et al. 2015). Based on the above information and the data presented in Fig. 2, biochars created at pyrolysis temperatures > 600 °C would typically be the most recalcitrant when placed in soils and, at the same time, among the most effective in mitigating N 2 O emissions (e.g., see: Borchard et al. 2018;Cayuela et al. 2015). Biochars created between 500 °C and 599 °C would typically have half-lives of 100-1000 years when placed in soils (Pariyar et al. 2020). Biochars created at temperatures < 500 °C should have relatively lower half-lives when placed in soils due to these materials only being partially thermochemically altered (e.g., see Table 2 in Spokas 2010). Similar observations have been found by Wolf et al. (2013) when combusting various feedstocks over increasing temperatures. Comparing 128 observations from 24 studies, Wang et al. (2016) found that the mean residence time of recalcitrant biochar C pools was 556 years. It is obvious that understanding pyrolysis temperature is important for the C sequestration potential of biochar ).

Total macro-nutrient concentrations
Total N and P contents are comparable between slow and fast pyrolysis, yet fast pyrolysis favored greater S, K, Ca, and Mg concentrations (Table 2). Feedstock choice dictated total N content, in the order of wood biochars < crop waste biochars < other grasses < manures/biosolids biochars. Greater biochar N content is likely a function of greater amino acids and proteins present in these materials (Tsai et al. 2012). As suggested by , this provides an advantage to greater N-containing biochars, such as manures/biosolids biochars, to act as a fertilizer (e.g., N) source. Contrarily, substantial N amounts are released as gaseous forms (e.g., NO, N 2 O, NO 2 , NH 3 , N 2 ; > 60%) with increasing pyrolysis temperature as illustrated in Table 2. The same response, at the same temperatures, were reported by Ippolito et al. (2015) when comparing pyrolysis temperature effects on N within a smaller biochar dataset (n < 100). Feedstock choice has been noted to typically be the primary determinant on biochar total S content (Cheah et al. 2014); our meta-analysis results suggest that only manures/ biosolids-based biochars contain a greater S content as compared to biochars made from other feedstocks (Table 2). Biochars containing elevated S content (e.g., manures; Sager 2012) have been shown to increase soil S bioavailability (Blum et al. 2013). Relatively lower pyrolysis temperatures (e.g., < 500 °C) have been shown to keep total biochar S content intact, yet greater temperatures can lead to gaseous S losses (Wang et al. 2013). However, data presented in Table 2 do not support this contention, showing no significant difference in S content when varying pyrolysis temperature. Pyrolysis temperature has also been shown to be a main determinant of the S forms found in biochar. Biochar S is mainly found in the organic fraction, such as dibenzothiophene (59% of total sulfur) and dibenzyldisulfide (14% of total sulfur), both of which can be bound to biochar-borne C (Cheah et al. 2014). As pyrolysis temperature increases, these organic-S bearing phases are converted to gases and lost (Carpenter et al. 2010). However, some S still remains, likely in a recalcitrant form bound to Ca, K, Mg, and Si at higher pyrolysis temperatures (> 800 °C) (Knudsen et al. 2004). Perhaps S data presented in Table 2 reflect these recalcitrant S forms present with respect to increasing pyrolysis temperature.
Biochar total P, K, Ca, and Mg concentrations are also affected by feedstock choice (Table 2). Manures/biosolidsbased biochars typically contained greater P, Ca, and Mg concentrations, other grass-based biochars contained appreciable K and Ca, while crop wastes-and wood-based biochars contained the least quantities of these four elements. In general, increasing pyrolysis temperature increased P and Ca content, had mixed effects with respect to Mg content, and had no significant impact on K concentration. The P, K, Ca, and Mg feedstock and pyrolysis temperature findings are likely a result of the initial elemental concentrations present in various feedstock, as well as increasing ash content with increasing pyrolysis temperature (Table 1). As ash content increases, Ca oxides, hydroxides, and carbonate mineral phases precipitate, which leads to increased total Ca concentration. In terms of P, as pyrolysis temperature increases biochar P content typically increases as well (e.g., Ippolito et al. 2015), likely due to P forming associations with carbonate phases. Although not strongly illustrated in Table 2, it is important to note that P-containing compounds can volatilize above 760 °C (Knicker 2007) causing total P concentrations to decrease.

Total micro-nutrient concentrations
Biochar total micro-nutrient concentrations, based on pyrolysis type, feedstock source, and pyrolysis temperature, are shown in Table 3. Only total Cu and Zn concentrations were affected by pyrolysis type, with fast pyrolysis favoring greater concentrations over slow pyrolysis. As with total macro-nutrient concentrations, total micro-nutrient concentrations were almost always significantly greater in manures/biosolids-based biochars as compared to other biochars. This is a function of manures/ biosolids feedstocks being fortified with micro-nutrients that were not assimilated by livestock or as a waste product from municipal/industrial sources (Sistani and Novak 2006). Thus, manure-based feedstocks contain inherently greater micro-nutrient concentrations than other feedstocks. This could also be a function of ash content present based on feedstock choice (Table 1) creating micronutrient oxide, hydroxide, and carbonate phases (e.g., CuO, Cu(OH) 2 , or CuCO 3 ; Ippolito et al. 2012;Zhang et al. 2017a). Likewise, increasing pyrolysis temperature was shown to affect ash content (Table 1), with relatively similar connections to increasing micro-nutrient concentrations as a function of temperature (Table 3).
Increasing micro-nutrient concentrations within biochars suggest that they may be used as a micro-nutrient fertilizer source. Chang et al. (2015) showed that a biochar containing Fe, Cu, Zn, Mn, and Co could potentially be used to supply micro-nutrients to plants. More specifically, Sigua et al. (2016) utilized poultry litter biochar addition to soil, observing an increase in extractable Fe, Mn, Cu and Zn by 19, 46, 68 and 32%, respectively, as compared to a control. Additionally, Zhao et al. (2018) reported that although a pig manure biochar contained excessive total Mn (1230 mg kg −1 ), Cu (780 mg kg −1 ), and Zn (1010 mg kg −1 ), their bioavailability was up to several orders of magnitude lower.

Biochar available nutrient analysis as a function of pyrolysis type, feedstock, and pyrolysis temperature
Biochar total elemental analysis only describes maximum nutrient concentrations present, which is not truly indicative of the nutrient availability when applied to soil. Available nutrients are those elements that may be readily absorbed by plants, as determined using various extractants (e.g., water, 1 M KCl, 0.01 M CaCl 2 , DTPA, Morgan, Mehlich, or Olsen extracting solutions). Although these extractants were originally developed for use in soils, they have been widely used with other soil amendments, including biochars and in biochar-amended soils (e.g., Lehmann et al. 2003;Lentz and Ippolito 2012). Changes in biochar available macro-and micro-elements, as a function of pyrolysis type, feedstock choice, and pyrolysis temperature, are shown in Table 4.

Pyrolysis type
Out of the ten available nutrients studied, only two were affected by pyrolysis type (Table 4). Biochars created via slow pyrolysis favored greater Fe and NO 3 availability over  fast pyrolysis. Depending on the extractant used, nutrient availability may be either over-or under-estimated (e.g., the use of Morgan or Mehlich extractants, which contain weak acids, for use with alkaline-containing materials could overestimate availability). In the current study, we grouped all data together and did not model effect of extractants used for nutrient availability determination; this would be worth exploring in the future.

Feedstock choice
In general, manure/biosolids-based biochars typically contained the greatest available nutrient concentrations as compared to other biochars (Table 4). This is not surprising given that these feedstocks contain, overall, greater nutrient contents as compared to other feedstocks (e.g., Williams et al. 2017). All other feedstock-based biochars generally were grouped together in terms of nutrient availability, potentially a consequence of lower ash content leading to lower nutrient retention (e.g., Williams et al. 2016).

Pyrolysis temperature
Increasing pyrolysis temperature increased only the macroelements P, K, Mg, and Ca (Table 4). As mentioned in the biochar total elemental analysis discussion (Sect. 4, above), as pyrolysis temperature increases, water, volatile bio-oil compounds, acids, organic surface functional groups and tars are lost, all of which contain H and O resulting in their losses as well (Ahmad et al. 2014;Cantrell and Martin 2012;Antal and Grønli 2003). This concentrates other elements in the final product, with macro-nutrient availability appearing to be a function of pyrolysis temperature. This concept is explored in more detail below. Ippolito et al. (2015) suggested that total elements present in biochars cannot accurately predict nutrient availability. However, their dataset utilized was based only on approximately 80 published biochar articles. It is worth to revisit this concept, as well as other correlations, based on nearly 5400 published articles reviewed in the current study. Thus, Pearson correlations between biochar available and total nutrient contents, pyrolysis temperature, pyrolysis type, or feedstock choice were first determined (Table 5). Most coefficients were not significant. However, some exceptions were present for predicting biochar nutrient availability based on total nutrient content, such as for K and Cu. Predicting biochar nutrient content solely on pyrolysis temperature or production technique appears poor; predicting nutrient availability based on feedstock choice appears somewhat more promising. However, in order to parse the data into something meaningful for the end-user, presented below are comparisons between total and available nutrients with respect to all biochar data, fast or slow pyrolysis, pyrolysis temperature, and feedstock choice. Only data where both total and available nutrients were presented (N, P, K, Ca, Mg, Fe, and Cu) in previously published works are shown.

All data: effects on nutrient availability
The total versus available N, P, K, Ca, Mg, Fe, and Cu data, for the complete dataset for each element, are presented in Figures S1A through S1G. The data suggest that predicting biochar plant-available concentrations from total elemental content is relatively poor for most elements (R 2 = 0.19,0.01,0.35,0.11,no fit,and 0.09 for N,P,K,Ca,Mg,and Fe,respectively). However, predicting plant-available biochar Cu from total Cu content had an R 2 = 0.97 utilizing an exponential growth model. Unfortunately, most biochar producers and end-product users are likely not as concerned about applying plant-available Cu to soils via biochar application as compared to other macro-element availability.
In terms of macro-element availability, most end-product users would likely be concerned about N availability and whether they need to add supplemental N fertilizers with biochars to offset potential negative crop N responses (e.g., Borchard et al. 2014;Lentz and Ippolito 2012). For example, in the western US, irrigated winter wheat or corn could require 85 and 235 kg N ha −1 , respectively, to maximize yields (Davis and Westfall 2014a, b). The NO 3 and NH 4 data in Table 4, along with data shown in Figure S1A, suggest that most biochars (except manures/biosolids biochars) contain available N contents below ~ 200 mg kg −1 ; at 200 mg of available N kg −1 biochar, an end-user might need to apply between 425 and 1175 Mg biochar ha −1 to supply the N 1 3 demands of winter wheat or corn, respectively. Unfortunately, microbial immobilization should also be accounted for, which would likely limit N availability, as described by Borchard et al. (2014). Furthermore, these biochar application rates to soils would equal ~ 19% and 52% by weight, well over the general application guidelines many biochar researchers are considering (~ 0.5% to 1% by weight, or 112 to 224 Mg ha −1 ). If biochar costs are considered (e.g., $500-$1000 USD Mg −1 ), crop-land applications would be economically unrealistic if it was just for an N fertilization effect (that is not given in woody biochars anyway due to the N being bound into the biochar matrix). Biochars would need to contain 2000 mg of available N kg −1 , or greater, to bring application rates near the 0.5-1% application level. According to data presented in Table 4 and Figure S1A, very few untreated biochars meet this requirement. Apart from N, taking a closer look into varying pyrolysis type, temperature, or feedstock choice may provide additional insight into overall biochar nutrient availability prediction. There may be ways to increase the ability of biochar to provide plant-available N, even though NO 3 and NH 4 content of freshly produced biochar is near-zero (e.g. Kammann et al. 2015;Haider et al. 2016). Biochar has the ability to capture NO 3 when co-composted (Hagemann et al. 2017a, (b); wood biochars have been shown to release up to 5000 mg NO 3 kg −1 following co-composting when repeatedly extracted (Kammann et al. 2015;Haider et al. 2016;Hagemann et al. 2017b). Biochar can also capture available N when biochar resides in soils for longer time periods (Haider et al. 2016(Haider et al. , 2017. In organic environments with alternating moisture regimes, an organo-mineral coating forms on biochar surfaces (Hagemann et al. 2017a). The coating is enriched in functional groups that improve the ability of biochar to capture and release nitrate (and to a lesser extent ammonium) (Kammann et al. 2015;Hagemann et al. 2017a). The mechanism, however, is not yet completely understood (Joseph et al. 2017) and deserves further attention.

Fast pyrolysis effects on nutrient availability
The effect of fast pyrolysis on the correlation between total and available biochar nutrients are presented in Figures S2. As compared to all data ( Figure S1), separating data by fast pyrolysis did not increase the predictive fit for most elements, and in some cases decreased the degree of fit. This suggests that predicting biochar nutrient availability based on fast pyrolysis is relatively poor. However, predicting available Cu based on total Cu in fast pyrolysis biochars still has a relatively good fit (R 2 = 0.74; exponential growth model). Unfortunately, as previously mentioned, most biochar end-users likely would not be concerned with Cu availability.

Slow pyrolysis effects on nutrient availability
The effect of slow pyrolysis on the correlation between total and available biochar nutrients are presented in Figures S3. As compared to all data ( Figure S1), separating data by slow pyrolysis increased the predictive fit for most elements. Furthermore, slow pyrolysis tended to fit data to a greater degree than fast pyrolysis. This may simply be due to the larger number of biochars created from slow pyrolysis as compared to fast pyrolysis, with gaps in the data reduced simply due to more literature information. Fitting biochar available nutrient to total nutrient content from fast to slow pyrolysis improved the predictive R 2 function for P (0.05-0.13), K (0.32-0.47), Ca (0.25-0.46), Mg (0.03-0.29), and Cu (0.74-0.98). Although these fits might not be perfect, the data suggest that predicting slow pyrolysis biochar's plant-available nutrient content, based on total nutrient content, is not out of the question (e.g., K, Ca, and Cu).

Pyrolysis temperature effects on nutrient availability
The effects of pyrolysis temperature (from < 300 °C to > 800 °C) on the correlation between total and available biochar nutrients are presented in Figures S4. In general, available K, Mg, and Fe content appear to increase with increasing pyrolysis temperature, then decrease with further rising pyrolysis temperatures > 800 °C. It also appears that, in general, temperatures < 300 °C or > 800 °C produce biochars with low available nutrient content.

Specific feedstock choice effects on nutrient availability
Based on Pearson correlations (Table 5), predicting biochar nutrient availability on feedstock choice appears somewhat promising. The effect of feedstock choice on the correlation between total and available biochar primary macro-nutrients (N, P, and K) are presented in Table 6, while all nutrient data (where both total and available nutrients were reported) are presented in Figures S5. The discussion below focuses solely on the predicted availability of N, P, and K in biochars created from specific feedstocks. However, the reader is encouraged to visit the supplemental figures for details regarding predicting N, P, K, Ca, Mg, Fe, and Cu availability.

Wood-based feedstocks
Previous research has focused a great deal of attention on utilizing hardwood and softwood feedstocks for biochar creation as compared to other materials. Based on a smaller dataset analyzed, Ippolito et al. (2015) showed that hardwood biochars would supply less than 0.002% of the total N, 2.2% of the total P, and 17% of the total K present. The authors also showed that softwoods could supply 27% and 6% of the total P and K present, respectively. Others have shown increases in available N, P, and K with softwoodbased biochars (Dieguez-Alonso et al. 2018). Based on the current dataset, if producing biochars from hardwoods, K availability can be predicted relatively well based on total content (R 2 = 0.87) as compared to N and P. If, however, one produces biochars from softwoods, the availability of both N and K may be predicted relatively well based on total content (R 2 = 0.68 and 0.54, respectively). In Sect. 6.1 above, it was stated that based on targeted biochar application rates of 0.5% to 1% (by wt.), one might need to apply at least 2000 mg of available N kg −1 to meet some of the major crop N demands (e.g., corn). Based on the predicted softwood N availability function (Table 6), it would require softwood total N contents to be ~ 7%. None of the softwood biochar data collected exceeded 4.5% ( Figure  S5A), and thus supplemental N fertilizers would be needed when utilizing softwood biochars applied to high N-demanding crops such as corn. Softwood biochars could likely be co-composted in order to capture available N, as described in Sect. 6.1 above. In addition, less N-demanding crops could potentially have N supplied by softwood biochars applied at reasonable application rates (i.e., 0.5-1% by wt.).
Applying K fertilizers to meet crop demands depends on soil K availability. Potassium applications could be as low as -56 kg K ha −1 (corn; Davis and Westfall 2014a) to as high as 250 kg ha −1 (corn or alfalfa; University of Delaware 2019; Lissbrant et al. 2009). We can predict soil K availability well with both hardwood (R 2 = 0.87) and softwood (R 2 = 0.54) biochars, with wood-based biochars averaging 1660 mg K kg −1 (Table 4). Using 1660 mg K kg −1 as a starting point, a biochar application rate of between 0.8 and 7% (by wt.) would be required to meet the aforementioned low-and high-end crop K demands. However, available K in hardwoods and softwoods can range from 0 up to either 13,000 or 4000 mg kg −1 , respectively ( Figure S5G). At almost 4-10 times greater than the average K concentration in wood-based biochars, application rates could easily be 4-10 times lower than the estimated 0.8-7% by wt. biochar application rates needed to meet crop K demands. Utilizing total K content in wood-based biochars could be used as a decision-making tool for K application and could specifically make hardwood and softwood biochar land applications attractive for supplying crop K demands globally.

Crop waste feedstocks
The use of agricultural crop waste products for biochar creation has occurred extensively throughout the world. The dataset created for the current study was dominated by a variety of crops, with the three major crops globally grown (i.e., corn, wheat, rice) utilized for predicting N, P, and K availability. One could potentially predict P availability marginally-to-well when utilizing any one of these three biochars (R 2 ranged from 0.48 to 0.82). Utilizing wheat straw biochars as a K source could be predicted fairly well based on total K present (R 2 = 0.66), while estimating N availability from these three crop-waste biochars appears to be weak at best (R 2 = 0.02-0.41).
As with any fertilizer, applying P fertilizers to meet crop demands is dependent on soil P availability. Phosphorus applications could be as low as 9 kg P ha −1 (band application to winter wheat; Davis and Westfall 2014b) to as high as 130 kg P ha −1 (broadcast application to corn; University of Delaware 2019). We can predict P availability fairly well with corn, wheat, or rice-based biochars, with crop waste averaging ~ 520 mg P kg −1 (Table 4). Using 520 mg P kg −1 as a starting point, a biochar application rate of between 0.8% and 11% (by wt.) would be required to meet the aforementioned low-and high-end crop P demands. Even though the available P in these three crop wastes ranges from 0 to 1800 mg kg −1 (Figure S5E), biochar applications to supply crop P needs would likely only be warranted for those crops with low P requirements. We can, however, predict wheat straw K availability fairly well. Wheat straw available K concentrations ranged from ~ 6000 to 30,000 mg kg −1 (Figure S5H), which were greater than those for hardwoods/softwoods. These results should make applying wheat straw-based biochars attractive for supplying most crop K demands globally and may solve a problem related to waste use in certain areas of the globe.

Manures/biosolids feedstocks
The use of manures or biosolids feedstocks for biochar creation has also been extensively studied globally. This may be driven by the fact that these feedstocks typically contain greater nutrient contents as compared to other feedstocks (e.g., Williams et al. 2017), and that pyrolyzing these materials represents a sound form of hygienization. Biochars made from these feedstocks have been shown to release available N (Liu et al. 2014;Sigua et al. 2016) in contrast to woody biochars. In the current study, N availability based on total biochar N content was predicted very well when utilizing either pig or cattle manure (R 2 = 0.99 and 0.99, respectively), but not with poultry or biosolids biochars. This suggests to biochar producers that, if they wanted to create a biochar that would act as an N fertilizer, that they either have to blend it with organic N-rich fertilizer materials (e.g. liquid manures, biogas digestate, vinasse, etc.) and use it as a carrier in low doses in the root zone  where 0.8-2 Mg ha −1 may suffice, or utilize either pig or cattle manure as a feedstock which might be also promising. As outlined in Sect. 6.1 above, based on a somewhat realistic biochar application rate (0.5-1.0% by wt.), one might need to apply at least 2000 mg of available N kg −1 to meet some of the major crop N demands (e.g., corn). Based on the predictive N functions in Table 6, pig or cattle manure would need to contain 3.8% or 11.7% total N, respectively. This total N content may be potentially feasible in pig manure ( Figure S5C) if the material were sufficiently dried to concentrate total N (swine manures typically have low dry matter contents of 15-20%; e.g., McFarland et al. 2012). However, reaching 11.7% total N in cattle manure is likely nearly impossible given total N content ranges in the material and the relatively high dry matter content present (50-80%, e.g., McFarland et al. 2012). This leaves pig manure as the only biochar feedstock in the current study to have the potential to sufficiently supply N for crop growth without mixing the biochar with organic fertilizer materials, or use it in their treatment (e.g. composting; see: Godlewska et al. 2017;meta-study: Zhao et al. 2020).
It may be somewhat surprising that the prediction of available N from poultry manure and biosolids biochars was relatively low given that these feedstocks may contain greater initial N contents as compared to other biochar feedstocks ( Figure S5C). However, the lack of prediction may simply be due to pyrolysis itself. Clough et al. (2013) noted that manures and biosolids pyrolysis can result in increasing aromatic and heterocyclic N structures within biochars, which are more difficult to degrade (i.e., recalcitrant).
It is also surprising to note that predicting either P or K availability was, at best, marginal with manures or biosolidsbased biochars. Again, this may simply be due to pyrolysis itself. Biochars derived from manures and biosolids contain greater ash content (Table 1). As previously mentioned, this ash contains oxides, hydroxides, carbonates, as well as silicate phases (Novak et al. 2019a), that can form recalcitrant associations with K. Manure and biosolids derived biochars also contain greater Ca concentrations as compared to other biochars (Table 1), which, in combination with elevated alkalinity, pH, and aromatic C has been shown to reduce P solubility (Ngatia et al. 2017). This may also be the reason for lack of P fit for other biochars above (Wang et al. 2013).
Finally, the information presented above may lend itself for creating biochars from a combined variety of feedstocks. For example, mixing pig manure with agricultural crop wastes and hardwood or softwood, then pyrolyzing the combined materials, might effectively supply N, P, and K, respectively. Other researchers have followed a similar approach, albeit post biochar creation. Novak et al. (2014) blended manure biochar containing excessive P, with a nutrient poor biochar to achieve an end-product with a more balanced nutrient content. Sigua et al. (2016) utilized a 50:50 mix of softwood biochar with poultry manure biochar, observing a 670% and 830% increase in soil P and K, respectively. The promise of mixing feedstocks for balancing nutrient availability in biochars could be potentially realized based on the data presented in the current manuscript.

Conclusions
Understanding the influences that pyrolysis type, pyrolysis temperature, and initial feedstocks have on final biochar properties can help researchers and practitioners create biochars to meet agricultural environmental demands. Based on ~ 5400 published articles and over 50,000 individual observations, this project makes inferences to further our understanding of biochar physicochemical properties from the broad to specific and minute perspective. As compared to fast pyrolysis, slow pyrolysis leads to biochars containing greater SSA, CCE, ash content, available Fe and NO 3 concentrations.
Pyrolysis temperature influences biochar stability, with temperatures > 500 °C generally leading to longer-term half-lives (> 1000 years). This, in combination with greater pyrolysis temperatures promoting more stable C structures, greater SSA, and potential improvements in soil aeration, percolation, infiltration, and overall structure, potentially suggests that greater pyrolysis temperatures may lead to long-term soil improvements and C storage.
Perhaps the most important influence on final biochar properties is feedstock choice. Wood-based feedstocks typically led to biochars containing the greatest SSA as compared to other feedstocks; this, in combination with pyrolysis temperature could greatly influence soil improvements. Crop-, grass-, and manures/biosolids-based feedstocks led to biochars containing elevated CECs as compared to wood-based biochars, which could affect nutrient sorption following land application. Based on the complete dataset collected, it appears possible to predict some plant-available biochar nutrients simply from total nutrient analysis. The collected data showed that we could reasonably predict (1) available N from softwood, corn, pig manure, and cattle manure biochars; (2) available P from corn, wheat, and rice straw/husk biochars; and (3) available K from hardwood, softwood, and wheat-derived biochars. This latter information could be useful when creating designer biochars for specific nutrient applications, simply by blending several feedstocks together. Based on this information, future research should test whether the available nutrient predictive functions, in combination with created mixed feedstock biochars, would hold true when placed within nutrient-poor soils. were funded by FACCE-CSA no 276610/MIT04-DESIGN-UPVASC and IT-932-16, MLC thanks the Spanish Ministry of Science, Innovation and Universities, project #RTI2018-099417-B-I00, cofinanced with EU FEDER funds and US colleagues (JN, JI and KS) were funded by The USDA-National Institute of Food and Agriculture (Project # 2014-35615-21971), USDA-ARS CHARnet and GRACENet programs -D4F greatly stimulated discussions. Any opinions, findings, or recommendation expressed in this publication are those of the authors and do not necessarily reflect the view of the USDA. This work was also partially supported by the National Natural Science Foundation of China under a Grant number of 41501339, 21677119, 21277115, 41301551, 21407123, Jiangsu Province Science Foundation for Youths under a grant number of BK20140468, sponsored by Qing Lan Project.
Funding Open access funding provided by Natural Resources Institute Finland (LUKE).
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://creat iveco mmons .org/licen ses/by/4.0/.