Scholarly article on topic 'Biochar stability in soil: meta-analysis of decomposition and priming effects'

Biochar stability in soil: meta-analysis of decomposition and priming effects Academic research paper on "Earth and related environmental sciences"

0
0
Share paper
Academic journal
GCB Bioenergy
Keywords
{""}

Academic research paper on topic "Biochar stability in soil: meta-analysis of decomposition and priming effects"

Received Date : Accepted Date : Article type :

19-Feb-2015 30-Mar-2015 Research Review

Type of Paper: Research Reviews

Biochar stability in soil: meta-analysis of decomposition and priming effects

Running head: Meta-analysis of biochar stability and priming

JINYANG WANGa b, ZHENGQIN XIONGb, YAKOV KUZYAKOVa c

a Department of Soil Science of Temperate Ecosystems, and Department of Agricultural Soil Science, University of Göttingen, 37077 Göttingen, Germany

b College of Resources and Environmental Sciences, Nanjing Agricultural University, 210095 Nanjing, China

c Institute of Environmental Sciences, Kazan Federal University, Russia Corresponding author: JINYANG WANG

This article has been accepted for publication and undergone full peer review but has not been through the copyediting, typesetting, pagination and proofreading process, which may lead to differences between this version and the Version of Record. Please cite this article as doi: 10.1111/gcbb.12266

College of Resources and Environmental Sciences, Nanjing Agricultural University, Nanjing

210095, China

Tel: +86 25 8439 5148

Fax: +86 25 8439 5210

E-mail: jyw217@gmail.com (J. Wang)

Keywords: pyrogenic organic matter, black carbon, soil respiration, priming effect, C sequestration, climate change

Abstract

The stability and decomposition of biochar is fundamental to understand its persistence in soil, its contribution to carbon (C) sequestration, and thus its role in the global C cycle. Our current knowledge about the degradability of biochar, however, is limited. Using 128 observations of biochar-derived CO2 from 24 studies with stable (13C) and radioactive (14C) carbon isotopes, we meta-analyzed the biochar decomposition in soil and estimated its mean residence time (MRT). The decomposed amount of biochar increased logarithmically with experimental duration, and the decomposition rate decreased with time. The biochar decomposition rate varied significantly with experimental duration, feedstock, pyrolysis temperature, and soil clay content. The MRTs of labile and recalcitrant biochar C pools were estimated to be about 108 days and 556 years with pool sizes of 3% and 97%, respectively. These results show that only a small part of biochar is bioavailable and that the remaining 97% contribute directly to long-term C sequestration in soil. The second database (116 observations from 21 studies) was used to

evaluate the priming effects after biochar addition. Biochar slightly retarded the mineralization of soil organic matter (SOM; overall mean: -3.8%, 95% CI=-8.1-0.8%), compared to the soil without biochar addition. Significant negative priming was common for studies with a duration shorter than half a year (-8.6%), crop-derived biochar (-20.3%), fast pyrolysis (-18.9%), the lowest pyrolysis temperature (-18.5%), and small application amounts (-11.9%). In contrast, biochar addition to sandy soils strongly stimulated SOM mineralization by 20.8%. This indicates that biochar stimulates microbial activities especially in soils with low fertility. Furthermore, abiotic and biotic processes, as well as the characteristics of biochar and soils, affecting biochar decomposition are discussed. We conclude that biochar can persist in soils on a centennial scale and that it has a positive effect on SOM dynamics and thus on C sequestration.

Introduction

Biochar, the carbonaceous residue of pyrolyzed organic materials under low oxygen conditions, has gained increasing attention in the last decade. Biochar is used for increasing soil carbon (C) sequestration (Lehmann et al., 2006; Schmidt et al., 2011), soil remediation (Ahmad et al., 2014; Zhang & Ok, 2014), greenhouse gas emissions mitigation (Woolf et al., 2010; Gurwick et al., 2013), as well as improving soil fertility and crop yields (Jeffery et al., 2011; Biederman & Harpole, 2013). Notably, its importance as a supply-side mitigation option is for the first time highlighted and summarized in the latest Intergovernmental Panel on Climate Change report (Smith et al., 2014).

A recent modeling study has demonstrated that sustainable biochar production could help reduce net greenhouse gas emissions by a maximum of 1.8 Pg CO2-Cequ annually without

penalizing ecosystem stability and food security; the inert C of biochar is the largest contributor because of its highly stable characteristic (Woolf et al., 2010). The biochemical stability of biochar is crucial mainly due to its longevity in soil after application and its long-term maintenance of fertility (Glaser et al., 2002; Lehmann et al., 2009). It remains unclear, however, to what extent biochar is degraded and what its concomitant effects on native soil organic matter (SOM) turnover and other cascading impacts are (Lehmann et al., 2011; Lorenz & Lal, 2014). Although the available studies have reviewed the fate of biochar in various environments (Schmidt & Noack, 2000; Forbes et al., 2006; Preston & Schmidt, 2006; Knicker, 2007, 2011), the current scientific understanding of biochar stability and persistence in soil remains sparse and incomplete. There are two reasons for this lack of knowledge about biochar decomposition in soil. Firstly, the contribution of biochar to total CO2 efflux from soil is too small to differentiate it from the background of high CO2 efflux from other sources, i.e. SOM, dead plant residues, DOC, etc. (Kuzyakov et al., 2009). Secondly, without information on the initial stock or input of biochar or other pyrogenic carbon in soils, any estimation of biochar decomposition based on its content in soil may be biased. Consequently, such approaches as chronosequences or 14C dating, as well as others, are strongly biased (Marschner et al., 2008; Ohlson et al., 2009; Singh et al., 2012a).

Improved approaches, such as using benzenepolycarboxylic acids (BPCAs) as molecular markers of biochar have been successfully used to quantify the amount of biochar and evaluate its sources in soils (Glaser et al., 1998; Brodowski et al., 2005; Hammes et al., 2008; Schneider et al., 2011). Although these approaches help assess potential sources of biochar in soil and the contribution of biochar compounds to SOM, they fail to quantitatively evaluate biochar

decomposition rates and dynamics in soils as BPCAs are only part of biochar C. This makes the assessment of biochar stability in soil impossible using these and the above-described approaches.

Ten years ago, a new approach based on application of 14C- (and/or 13C-labeled) biochar was suggested to trace its turnover and transformation products against the background of other CO2 sources from soil (Kuzyakov et al., 2009, 2014). Currently, changing biochar input from C3 to C4 plant material or vice versa or also applying artificially 13C-enriched or 13C-depleted biochar was introduced to estimate its turnover in various soils (e.g., Glaser & Knorr, 2008; Bruun et al., 2014; Singh et al., 2014). Compared with previously used methods to assess biochar turnover, these 13C and 14C techniques improve our understanding of biochar dynamics in soils and their cascading effects on ecosystem processes.

Nonetheless, the extent of biochar decomposition and its mean residence time (MRT) in soils-means its stability-remains nearly unknown. The overall objective of this meta-analysis is therefore to synthesize those biochar decomposition studies in which 13C or 14C isotopes were used. We first collected the published studies and established a database on the magnitude and rates of biochar decomposition. This database contains 128 observations of biochar-derived CO2 from 24 studies (Appendix Table S1). We then evaluated both the factors and processes affecting its degradation and stability. Next, we used this database to estimate the MRT of biochar in soil using the double exponential decay model. Finally, the effect of biochar addition on native SOM decomposition (e.g., priming effect) was assessed by a meta-analysis based on the studies evaluating biochar stability.

Materials and methods Data compilation

The synthesis was performed on published data of biochar decomposition collected using ISI Web of Science and Google Scholar. The following criteria were applied to select appropriate studies. (1) To reliably evaluate biochar decomposition in soil, we restricted the data collection

13 13 14

to studies employing stable isotopic techniques ( C natural abundance or C labelling) or C labelling to separate biochar-derived CO2 from CO2 derived from native SOM or other C sources (e.g., plant residues). (2) If more than one publication presented results from the same experiment, we selected data only from the studies with the longest experimental duration (Stewart et al., 2013; Fang et al., 2014a; Maestrini et al., 2014a). (3) The analysis was focused solely on biochar produced by pyrolysing organic materials under dry conditions. Hence, the other types of carbonized biomass (e.g., hydrothermal char) were not considered (Malghani et al., 2013). (4) To avoid the confounding effects of temperature on biochar decomposition, we excluded incubation studies at temperatures exceeding 40 °C (Fang et al., 2014a). (5) A few studies were excluded that reported negative values of biochar decomposition because of high variation by application of 13C natural abundance or because isotopic fractionation was not considered (Cross & Sohi, 2011). After these restrictions, data were extracted from 24 studies with 128 observations that reported biochar decomposition (Appendix Table S1). The main characteristics of included studies are presented in Appendix Table S1. Experimental factors potentially affecting biochar decomposition were summarized by a number of categorical variables: type of experiment, experimental duration, feedstock, pyrolysis type and temperature of biochar production, application amount, ash content, land use, as well as soil clay content (Table 1).

In addition, a meta-analysis was performed to explore the effect of biochar-induced priming on SOM mineralization. We restricted our database of biochar-induced priming effects on those studies that simultaneously investigated biochar decomposition as described above and Jones et al. (2011). This yielded a database consisted of 116 observations from 21 studies, which was used to evaluate the priming effect after biochar addition (Appendix Table S1).

Statistical analysis

Most of the data included in our analysis were extracted from tables and figures in published studies. Few data not presented in the literature were obtained via personal communication (Farrell et al., 2013). The data presented in graphs were extracted by digitizing the figures using g3data software (http://www.frantz.fi/software/g3data.php). When means and standard errors (SE) were reported, the standard deviations (SD) were calculated as: SD = SEyfn, where n is the number of replicates. All values of biochar decomposed amounts (% of added) were recalculated and expressed as percent per day, being referred to as the biochar decomposition rate. The Kruskal-Wallis one-way analysis of variance (ANOVA) was used to investigate differences in biochar decomposition rate between various groups of each categorical variable (Table 1). A nonparametric multiple comparison test (Dunn's test) was performed to compare the difference in the sum of ranks between the groups with the expected average difference.

The MRT of biochar for each study was calculated as a reciprocal of the decomposition rate. To estimate decomposition rates and MRTs of biochar pools in soil, a double first-order exponential decay model was used to fit the experimental data of the remaining amount of

biochar against experimental duration using a least squares optimization routine (Jones et al., 2011):

y — a X expi-&1 xti -+■ b X exp*-^x°

where y is the amount of biochar remaining in the soil, t is time, a and b represent the size of labile and recalcitrant C pools of biochar, respectively, and k1 and k2 are the exponential coefficients for these two pools. The MRTs for labile and recalcitrant C pools of biochar were calculated as reciprocal of k1 and k2, respectively. These statistical analyses were performed using STATISTICA v.10 (StatSoft. Inc., Tulsa, OK, USA).

The meta-analysis followed the approach described by Hedges et al. (1999). We used the natural log-transformed response ratio: InRR = In (Xt/Xj

where Xt and Xc are the means of SOM-derived CO2 emission with and without biochar addition, respectively. The positive and negative log response ratios indicate SOM-derived CO2 emissions increased and decreased after biochar addition, respectively. Similar to previous ecological studies, a mixed effects model was used with the assumption that differences among studies within a group are due to both sampling error and random variation (Gurevitch & Hedges, 1999). For studies where neither the SD nor the SE were available in the original papers, two weighting functions were used: (i) each individual observation was weighted by the reciprocal of the pooled variance as is conventionally done in meta-analysis (Curtis & Wang, 1998); and (ii) the data were also weighted by experimental replication, using the function FN = (nt x nc)/(nt + nc), where nt and nc represent the number of replicates for the treated and control groups, respectively (Adams et al., 1997). The total heterogeneity of effect sizes among studies (QT) is

partitioned into within-group (Qw) and between-group (Qb) heterogeneity. The Q statistic approximately follows a Chi-square distribution, which allows a significance test of the null hypothesis that all response ratios are equal. A Qb larger than a critical value indicates that there has significant difference between groups. Mean effect sizes and 95% confidence intervals (CI) were generated using MetaWin 2.1 (Sinauer Associates Inc., Sunderland, MA, USA). The response ratio (RR) and CI of treatments presented were back-transformed from lnRR. For each case, biochar treatment was considered significant if the 95% CI of response ratio did not overlap with zero. Responses of groups were considered different if their 95% CIs did not overlap. Graphs were plotted in SigmaPlot v.12.0 (Systat Software Inc., San Jose, CA, USA).

Results and Discussion Decomposition amounts and rates

The decomposed amounts of biochar increased logarithmically with experimental duration (Fig. 1, left), suggesting that part of biochar is biologically available and can be slowly degraded by soil microorganisms. Nonetheless, the high uncertainty of the fitting curve does not enable predicting the decomposed amounts over periods longer than a few years. This finding is consistent with Hammes et al. (2008), who found that the proportions of less-condensed biochar C decreased when compared the quality of BPCA between archived and modern soil samples in a Russian steppe.

In contrast to the amounts, the decomposition rate of biochar showed a clear decrease pattern with experimental duration (Fig. 1, right). Overall, the biochar decomposition rate was

calculated at a mean and median of 0.013% day-1 and 0.0046% day-1, respectively. The experimental duration varied very strongly, ranging from two weeks to 8.5 years. All observations were divided for four subgroups based on the experimental duration. The mean value of biochar decomposition rate for studies lasting less than 0.5 year was 0.023% day-1, i.e. more than four times faster than that from studies lasted longer than one year (0.005% day-1; Fig. 2). Despite the major content of aromatic C, the results of short-term studies mainly represent the decomposition of the non- or less-condensed fraction of biochar. Results from several studies (Appendix Table S1) imply that this initially intensive decomposition disappeared after 2 years and was maintained at a very low level over prolonged periods (Fang et al., 2014a, b; Kuzyakov et al., 2009, 2014). Considering that the longest biochar decomposition study (8.5 years) is much shorter than its assessed MRT, all available data reflect the decomposition of only the most labile biochar compounds. The experimental duration is one of the most important factors affecting biochar decomposition rates in soils. The low determination coefficient, however, suggests that decomposition may also be affected by factors other than the duration of experiment.

Biochars derived from various sources have different chemical compositions that affect biological degradability (Lehmann et al., 2011). The crop-derived biochar decomposed faster (mean: 0.025% day-1) tha n other feedstocks investigated (Fig. 2). The grass-derived biochar is also characterized by a lower C content and is less condensed, meaning that its decomposition rate is expected to be comparable to crop-derived biochar (Hilscher et al., 2009; Knicker, 2010). In our analysis, however, the biochar derived from grasses had a much lower values (mean: 0.007% day-1) tha n those from crop. This inconsistency is probably because one third of the

group of grass-derived biochar studies was conducted with the longer period (Kuzyakov et al., 2009, 2014). In addition, wood-derived biochar had the slowest decomposition rate (mean: 0.004% day-1) as a result of its high C content (means: 66.4% vs. 59.8% for crop and 64.6% for grass; P=0.02), corresponding to a higher aryl C content relative to both crop- and grass-derived biochars (Hilscher et al., 2009; Singh & Cowie, 2014).

Increasing the pyrolysis temperature of biochar production significantly decreased its decomposition rate (Fig. 2). For example, biochar produced at 200-375 °C decomposed faster than that pyrolyzed above 375 °C. Nonetheless, based on the available dataset involving 13C and 14C studies, no significant difference between biochar groups produced at higher temperatures was obtained (except that produced at 1200 °C). The decomposition rate decreases with increasing pyrolysis temperature because of the increased portion of aromatic C and the higher degree of aromatic condensation as confirmed previously (Singh et al., 2012b; Crombie et al., 2013).

A slower decomposition rate was common in soils with the highest clay content (mean: 0.003% day-1, 40-70%), whereas no significant difference was found among soils with lower clay content (Fig. 2). Moreover, aged biochar mineralized at similar rates in a chronosequence of high-biochar-containing Anthrosols with various clay contents ranging from 9% to 36% (Liang et al., 2008).

The MRT varied greatly among studies, ranging from less than three up to 891 years, with an average of 107 years. Similar to other organic materials added into soils, biochar decomposition can be approximated well by biphasic patterns. The commonly used double first-order exponential decay model yielded rates for labile and recalcitrant biochar pools of

0.0093% day-1 and 0.0018% year-1, respectively. This corresponds to MRTs of 108 days and 556 years, respectively (Fig. 3; Table 2). The labile and recalcitrant pools represent 3% and 97% of biochar C, respectively, which are close to the results of the longest study (Kuzyakov et al., 2014). The calculated MRTs for both C pools of biochar are much shorter than previously assumed in both soils (Schmidt et al., 2002; Singh et al., 2012a) and ocean sediments (Masiello & Druffel, 1998). This is because most of the results in the database are derived from incubation studies (Appendix Table S1), the durations of which are much shorter than the MRTs of biochar. The optimal conditions and short-term duration of the experiment led 1) to faster decomposition rates than those under non-optimal field conditions, and 2) rates that were obtained mainly for labile C pool of biochar. Because the number of available studies is limited, there were no significant differences in biochar decomposition rate between lab incubations and field results (Table 1). Consequently, the available field observations do not support the standpoint that biochar can persist in soils for several thousands of years. Furthermore, current analyses of forest soils in northern Europe and Russian steppes also indicated a centennial-scale turnover of biochar (Hammes et al., 2008; Ohlson et al., 2009). Serious caution is therefore needed when extrapolating results of biochar decomposition from incubation studies to real field conditions. Despite the estimated shorter MRTs for both biochar-C pools than previous findings, nearly all studies demonstrate that all biochar types are much more refractory than their precursors and natural SOM (e.g., Bruun et al., 2008; Maestrini et al., 2014a; Yin et al., 2014).

Biochar-induced priming effect

The second dataset was used for the meta-analysis to evaluate the response of biochar-induced priming effect on SOM decomposition. Averaged across all studies, biochar addition induced a negative priming effect, but this effect was not statistically significant (overall mean=-3.8%, 95% CI= -8.1-0.8%; Fig. 4). This result disagrees with the overview of Maestrini et al. (2014b), who reported a 15% increase of native SOM mineralization following biochar addition. This contradiction may primarily reflect the double number of studies in our database, which included 116 observations from 21 studies. Moreover, both positive and negative priming effects can coexist, whereby the latter acts over the longer term (Maestrini et al., 2014b).

The biochar-induced priming effects on SOM mineralization depends on various factors including the experimental duration, the biochar characteristics, and soil properties. Except for the type and duration of experiment, along with land use, the between-group heterogeneity (Chi-square test) for the remaining categorical variables was significant for the biochar-induced priming effect (Table 1). Approximately half of the studies shorter than 0.5 year showed a significantly negative priming effect (-8.6%) following biochar addition, whereas no significant effect was found in studies longer than 0.5 year (Fig. 4). Preferential substrate utilization of easily available C in biochar by soil microorganisms is responsible for the significant negative priming effect (Blagodatskaya & Kuzyakov, 2008). This is also supported by the fastest decomposition rate of biochar in those studies shorter than 0.5 year (Fig. 2). Hence, we conclude that the priming after biochar addition is induced by its easily available pools. Considering the common very high application rates of biochar (Fig. 4; Jeffery et al., 2011), even

a small available pool (ca. 3%; Table 2) is sufficient to be the source for preferential substrate utilization.

The direction and magnitude of the priming effect triggered by biochar addition varied strongly with the feedstock of biochar: the crop-derived biochar addition induced the highest negative priming (-20.3%) compared to other feedstocks. As demonstrated earlier, these results are due to the higher degradability of biochar C produced from crop compared to other sources, especially wood (Fig. 2). Biochar produced either by fast pyrolysis or at the lowest pyrolysis temperature (200-375 °C), as well as a small application amount (0.1-1%), induced significant negative priming compared to other groups for each categorical variable. More bioavailable C present in biochar may result from either fast pyrolysis or lower pyrolysis temperature (Spokas et al., 2011; Bruun et al., 2012). Different ash contents in biochar resulted in contrasting priming effects, but this difference was insignificant. In contrast to most studies with positive priming, biochar applied to sandy soils (clay content < 10%) led to a significantly and strongly positive priming effect (20.8%). This indicates that biochar stimulates microbial activities especially in soils with low fertility.

We conclude that priming induced by biochar addition is slightly negative because of the preferential utilization of easily available substrates from incomplete pyrolysis. Although this portion is small (ca. 3%; Table 2), very high biochar application rates are sufficient to trigger microorganisms to switch and thereby decrease native SOM decomposition. Other potential reasons, e.g., the presence of certain toxic compounds (phenolic nature) in biochar, may be responsible for negative priming and need to be clarified in further studies. For all studies, we definitely conclude that biochar application has a positive C balance at least within 10 years

(Kuzyakov et al., 2014). This means that adding biochar to soil increases the C content even by positive priming in some soils.

Factors and processes affecting biochar stability

Biochar and soil properties

The overwhelming evidence suggests that biochar is very slow decomposed in soils. Understanding the factors affecting its degradation is thus critical to maximize the long-term C sequestration potential and for soil fertility improvements in the context of changing climate and intensified land use. The intrinsic biochemical recalcitrance of biochar is primarily determined by the feedstocks and pyrolysis conditions (Fig. 5, left). Biochar produced from residues of crops and grasses is generally more degradable than that from wood, which is probably attributed to inert properties of various feedstocks, such as the high lignin content (Bird et al., 1999). Pyrolysis temperature significantly affects the proportion of nonaromatic and aromatic C, as well as the degree of condensation of aromatic C in biochar (Preston & Schmidt, 2006). The molar ratio of oxygen to carbon (O:C) is a function of pyrolysis temperature, making the O:C ratio positively correlated with biochar degradation (Spokas, 2010; Cross & Sohi, 2013; Crombie et al., 2013). Current studies have proposed that biochar, as one of the not mineral-associated SOM components, is persistent in soils, perhaps reflecting the high aromaticity of black C (Marschner et al., 2008). Supporting this, glycolipids, phospholipids, and neutral lipids (and probably carbohydrates) were the major constitutes of biochar decomposed over 3.5 years of incubation, whereas the condensed aromatics (BPCA) remained nearly unaffected during this period (Kuzyakov et al., 2014). Hence, the pyrolysis temperature is critical for the C sequestration potential of biochar.

The initial quality of the SOM also affects biochar degradation (Fig. 5, left). In SOC-poor soils, labile biochar compounds deliver substrate for microorganisms within the first few months to years. For SOC-rich soils, however, biochar can stimulate SOM mineralization (Cross & Sohi, 2011; Stewart et al., 2013). Biochar strongly sorbs organics in soils, especially low molecular weight organics, which then become unavailable for microorganisms (Zimmerman et al., 2011; Maestrini et al., 2014b). At the same time, the sorption of organics can block both the enzymes from accessing the biochar surface as well as oxidation and dissolution of biochar components (Kasozi et al., 2010; Bailey et al., 2011). Hence, biochar degradation depends on the soil characteristics, but the details of these interactions still require specific mechanistic investigations. Repeated addition of glucose or plant residues strongly stimulates biochar degradation (Hamer et al., 2004; Kuzyakov et al., 2009; Keith et al., 2011). The accelerated decomposition is ascribed to co-metabolism (Fig. 5, right; Kuzyakov et al., 2009) and continues over short periods of weeks to 2-3 months. Although we recognize the stimulated degradation due to other C input (priming on biochar), our current understanding of the effects of living plants and rhizodeposition on biochar decomposition is limited (Whitman et al., 2014). As biochar decomposition is stimulated by a continuous release of rhizodeposits, future studies on the mechanisms of biochar priming in soil-plant systems are required.

Abiotic factors and processes

Biochar is characterized by high hydrophobicity (corresponding to the insolubility of its compounds) and therefore affects soil hydrophobicity (Schmidt, 2004; Knicker, 2011). The hydrophobicity contributes to aggregate formation and consequently to the spatial inaccessibility of biochar particles and other compounds (Fig. 5, right). Surface runoff has been

considered as the largest biochar flux in a sandy Oxisol over 2 years, especially during the intense rain season (Major et al., 2010). Very weak biochar bounding on clay minerals and its light nature could lead to preferential surface erosion on slopes compared to other organics (Rumpel et al., 2006).

Biochar stability in soils may also be connected with spatial inaccessibility (physical protection) and organo-mineral interactions, particularly during the later decomposition phase (Fig. 5, right; Glaser et al., 2000; Lützow et al., 2006). Examining the distribution of black carbon in soils revealed that the greatest and smallest concentrations of black carbon occurred in the <53 pm and >2 mm fractions, respectively (Brodowski et al., 2006). Consequently, physical inclusion, particularly within microaggregates, at least partly contributes to the long MRTs of soil-inherent black carbon. The most biochar was recovered in particulate organic matter during the early incubation period, and this portion decreased over time (Hilscher & Knicker, 2011; Herath et al., 2014). Supporting of these findings, after 10 months in a temperate forest, 74%, 20%, and 6% of biochar C were recovered in the free light fraction, aggregate-occluded, and mineral-associated fractions, respectively (Singh et al., 2014). A fast association of biochar with the mineral fraction was found, but this portion decreased with time (Hilscher & Knicker, 2011). The formation of organo-mineral complexes is attributed to partial biochar oxidation, which provides negative surface charges of initially hydrophobic material (Glaser et al., 2000; Hilscher & Knicker, 2011). Nevertheless, more studies are necessary to prove contradictory results of organo-mineral interactions on biochar stability in soils (Glaser et al., 2000; Marschner et al., 2008).

Previous studies confirm that abiotic processes such as carbonate dissolution and chemisorption play an important role in biochar degradation (Fig. 5, right; Cheng et al., 2006; Zimmerman, 2010; Jones et al., 2011). For example, the CO2 released from incubations under abiotic conditions (sterilized sand and nutrient solution) amounted to more than half that of microbially inoculated incubations over a one-year period (Zimmerman, 2010). Nonetheless, such results remain to be proven in rigorous experiments because nearly all transformation processes in soils are driven by microorganisms (Kuzyakov & Blagodatskaya, 2015). In an artificial weathering approach, physical disintegration is suggested to pronouncedly contribute to the considerable biochar loss, though the values varied greatly with the feedstock and pyrolysis temperature (Spokas et al., 2014). Abiotic oxidation and physical disintegration can facilitate the subsequent microbial metabolization of highly aromatic biochar, such as by forming smaller particles and further increasing microbial accessibility.

Very few studies conducted at temperatures above 30 °C concluded that biochar degradation is accelerated with increasing incubation temperature (Fig. 5, left; Nguyen et al., 2010; Zimmerman et al., 2012; Fang et al., 2014a, b). Oxidation is the most important decomposition mechanism, mainly through increases in O and hydrogen (H) contents, decrease of C content, and formation of O-containing functional groups (-COOH, =C=O, -OH), providing surface negative charges. Comparing biochar mineralization in the field during different periods (e.g., 1 vs. 11 years), the relative contributions of biotic and abiotic factors to biochar-derived CO2 might be modified by temperature, but depends on the biochar age (Zimmerman et al., 2012). Furthermore, decomposition is controlled by site-specific climatic conditions, particularly temperature (Preston & Schmidt, 2006). Analyzing the historical black carbon samples along a

climosequence revealed a strongly positive relationship between black carbon oxidation and mean annual temperature (Cheng et al., 2008). Finally, the temperature sensitivity of SOC mineralization of biochar-amended soils may decrease with increasing incubation temperature (Nguyen et al., 2010; Fang et al., 2014a).

Moisture plays an important role in biochar degradation (Fig. 5, left). A modeling study of biochar movement in soil indicates that biochar mineralization was quite sensitive to variation in moisture (Foereid et al., 2011). Biochar degradation was strongly increased under unsaturated and alternating saturated-unsaturated conditions compared to saturated conditions (Nguyen & Lehmann, 2009). These authors argued that enhanced biochar mineralization was attributed to the increase in carboxylic and -OH functional groups under unsaturated and alternating conditions. Oxygen depletion contributes to the increased MRT of biochar in waterlogged ecosystems such as sediments and paddy soils (Knicker, 2011). For example, the formation of a high content of aromatic and carboxyl C in ancient Chinese paddies is attributed to the cultivation mode of paddies involving fire under drying-rewetting alternation conditions (Cao et al., 2006; Hu et al., 2009).

Biotic processes

Despite certain abiotic processes, biochar decomposition in soils is primarily microbially mediated (Fig. 5, right). Solid-state 13C NMR spectroscopy showed alteration of biochar, including the formation of O-containing functional groups and the loss of aryl C structures; this was attributed to biotic oxidation and degradation (Hilscher et al., 2009). Increasing evidence suggest microbial degradation of biochar in soil. After 3.5 years of incubation, biochar-derived C in microbial biomass ranged between 0.3% and 0.95% of the 14C input, whereby these

proportions gradually decreased over time (Kuzyakov et al., 2014). Nonetheless, we assume that this 14C incorporation was mainly connected with microbial utilization of incompletely pyrolized products. More specifically, the utilization of biochar-C is also confirmed through its incorporation into microbial PLFA structure (Santos et al., 2012; Farrell et al., 2013). Grampositive bacteria are known to be significant consumers of aromatic and recalcitrant compounds, and they are confirmed to be responsible for the early utilization of biochar C (Santos et al., 2012; Farrell et al., 2013). This confirms that incorporation into microorganisms mainly involves the labile organics still present after pyrolysis. At the later incubation period, the incorporation of biochar C into the fungal PLFAs became dominant over bacterial PLFAs (Farrell et al., 2013). In addition, the differently promoted growth of some microbial groups may be connected with the feedstocks of biochar (Steinbeiss et al., 2009). Despite due to the lack of available data about saprophytic decomposition of biochar, there is evidence that biochar can be decomposed by saprophytic fungi (Wiedner & Glaser, 2013).

The biotic degradation mechanisms remain unclear. Considering the high molecular weight of biochar, its non-regular aromatic structure, hydrophobicity and insolubility, we expect that only exoenzymes or fungi hyphae may effect initial degradation. Oxidative exoenzymes, such as peroxidases, phenol oxidases and laccases, are commonly responsible for degradation of polyphenolic compounds such as lignins, tannins, and suberins, etc. (Baldrian, 2006; Blagodatskaya & Kuzyakov, 2008) and can thus also oxidize complex biochar molecules. Because these enzymes require abundant oxygen, biochar decomposition under anoxic conditions (e.g., flooding of paddy soils) is much slower than in well-aerated soils. Considering that fires are natural events in nearly all ecosystems, microorganisms should have at least some

adaptations for utilizing pyrogenic products. We are, however, unaware of specific enzymes released by microorganisms for decomposing pyrogenic compounds.

Indirect processes

Translocation of biochar into deeper horizons can increase its persistence in soil (Lorenz & Lal, 2014). Biochar may be moved to subsoil by bioturbation or particulate transport (Fig. 5, right; Skjemstad et al., 1999). Leifeld et al. (2007) estimated migration rates of 630-1160 mm year-1 in peat soils, where 21-69% migrated below the incorporation depth of 0.3 m, down to a maximum of 1.14 m over 95 years. Such extreme migration rates, however, may reflect strong overestimation due to unclear sources of biochar over longer periods. After 10 months in situ in a temperate forest, 3-4% of the applied biochar was recovered below the application depth, suggesting very slow downward movement (Singh et al., 2014). This redistribution mainly involves bioturbation by soil animals such as earthworms. In contrast, only 0.02-0.45% of applied biochar in a sandy Oxisol moved over 2 years below 0.3 m depth as particulate organic C and DOC (Major et al., 2010). Those inconsistencies between studies assessing downward redistribution are connected with 1) uncertainties of previous source and deposition estimations of biochar, 2) intensity of bioturbation by burrowing animals, and 3) possible cracks and pores in soil related to periodic drying/rewetting events. Last, but not the least, other indirect processes such as freezing/thawing and swelling/shrinking dynamics may significantly affect on biochar particle size (Fig. 5, right). These processes may directly contribute to the physical breakdown of biochar. Moreover, reduced particle size would facilitate further microbial decomposition by increasing the accessible surface area.

Implications of biochar stability for C sequestration in soil

The total area burned and fire intensity are projected to increase in a warming climate. This will not only increase CO2 emission, but also promote the conversion of biomass to pyrogenic organic matter from wild and anthropogenically induced fires. Furthermore, the climate-negative C budget of biochar together with its properties for fertility improvements lead to the expectation of further increasing biochar application to agricultural soils (Glaser et al., 2002; Lehmann, 2007). Considering its acknowledged importance as a C sink, the extremely slow decomposition of biochar may have direct and indirect implications for the C sequestration potential in soils. Biochar input is the direct contribution to C sequestration in soils (Fig. 6). Our current meta-analysis shows that biochar has a very low turnover rate (overall median: 0.0046% day-1). The faster turnover of the labile pool is insignificant, as this pool amounts to only 3% of the biochar C. Consequently, 97% of the added biochar can persist in soils on a centennial scale. Therefore, its high recalcitrance plays an important role in abating climatic change directly by sequestering C from biochar input (Woolf et al., 2010). Indirect effects are based on the impact of biochar on SOM decomposition, sorption of easily available organics, and higher C input by plants into the soil. Pooling all studies together, biochar addition clearly contributes to a negative priming. This is also partly supported by the finding of Maestrini et al. (2014b) that biochar amendment may result in a negative priming effect over long-term period. Consequently, besides directly adding recalcitrant C to soil, biochar indirectly affects C sequestration by lowering SOM decomposition rates, at least for several years after application.

Finally, biochar application can suppress the mineralization of labile organic matter via sorption: this gradually increases C stabilization and thereby enhances SOC stocks (Fig. 6; Liang et al., 2010; Keith et al., 2011; Maestrini et al., 2014b). Moreover, by improving soil fertility,

biochar can stimulate plant growth, boosting the belowground biomass input by roots and rhizodeposition. Last, but not least, the beneficial effects of biochar amendment on plant production may directly stimulate more C input into soils via plant residue return and rhizodeposition (Biederman & Harpole, 2013). Biochar application decreases soil bulk density and improves biotic soil properties (Lehmann et al., 2011). We therefore expect not only better root growth, but especially better root distribution and allocation of more C into subsoil. Overall, our review demonstrates that biochar remains in soils for a very long time, and that its stability could make it a promising tool to help abate climatic change by long-term sequestering C in soils.

Acknowledgements

We would like to acknowledge the work carried out by the researchers whose published data were included in this meta-analysis. We also thank two anonymous reviewers for their helpful comments to improve this manuscript. JY is grateful to the China Scholarship Council for providing fund. This work was jointly supported by the National Science Foundation of China (41171238, 41471192) and the Graduate Research and Innovation Projects for College Graduates of Jiangsu Province (CXZZ13_0298).

References

Adams DC, Gurevitch J, Rosenberg MS (1997) Resampling tests for meta-analysis of ecological data. Ecology, 78, 1277-1283.

Ahmad M, Rajapaksha AU, Lim JE et al. (2014) Biochar as a sorbent for contaminant management in soil and water: a review. Chemosphere, 99, 19-33.

Bailey VL, Fansler SJ, Smith JL, Bolton JrH (2011) Reconciling apparent variability in effects of biochar amendment on soil enzyme activities by assay optimization. Soil Biology and Biochemistry, 43, 296-301.

Baldrian P (2006) Fungal laccases-occurrence and properties. FEMS microbiology reviews, 30, 215-242.

Biederman LA, Harpole WS (2013) Biochar and its effects on plant productivity and nutrient cycling: a meta-analysis. GCB Bioenergy, 5, 202-214.

Bird MI, Moyo C, Veenendaal EM, Lloyd J, Frost P (1999) Stability of elemental carbon in a savanna soil. Global Biogeochemical Cycles, 13, 923-932.

Blagodatskaya E, Kuzyakov Y (2008) Mechanisms of real and apparent priming effects and their dependence on soil microbial biomass and community structure: critical review. Biology and Fertility of Soils, 45, 115-131.

Brodowski S, John B, Flessa H, Amelung W (2006) Aggregate-occluded black carbon in soil. European Journal of Soil Science, 57, 539-546.

Brodowski S, Rodionov A, Haumaier L, Glaser B, Amelung W (2005) Revised black carbon

assessment using benzene polycarboxylic acids. Organic Geochemistry, 36, 1299-1310.

Bruun EW, Ambus P, Egsgaard H, Hauggaard-Nielsen H (2012) Effects of slow and fast pyrolysis

biochar on soil C and N turnover dynamics. Soil Biology and Biochemistry, 46, 73-79.

Bruun S, Clauson-Kaas S, Bobul'ska L, Thomsen IK (2014) Carbon dioxide emissions from biochar

in soil: role of clay, microorganisms and carbonates. European Journal of Soil Science, 65, 52-59.

Bruun S, Jensen ES, Jensen LS (2008) Microbial mineralization and assimilation of black carbon: Dependency on degree of thermal alteration. Organic Geochemistry, 39, 839-845. Cao ZH, Ding JL, Hu ZY et al. (2006) Ancient paddy soils from the Neolithic age in China's Yangtze River Delta. Naturwissenschaften, 93, 232-236.

Cheng CH, Lehmann J, Thies JE, Burton SD (2008) Stability of black carbon in soils across a climatic gradient. Journal of Geophysical Research: Biogeosciences (2005-2012), 113(G2). Cheng CH, Lehmann J, Thies JE, Burton SD, Engelhard MH (2006) Oxidation of black carbon by biotic and abiotic processes. Organic Geochemistry, 37, 1477-1488.

Crombie K, Masek O, Sohi SP, Brownsort P, Cross A (2013) The effect of pyrolysis conditions on

biochar stability as determined by three methods. GCB Bioenergy, 5, 122-131.

Cross A, Sohi SP (2011) The priming potential of biochar products in relation to labile carbon

contents and soil organic matter status. Soil Biology and Biochemistry, 43, 2127-2134.

Cross A, Sohi SP (2013) A method for screening the relative long-term stability of biochar. GCB

Bioenergy, 5, 215-220.

Curtis PS, Wang X (1998) A meta-analysis of elevated CO2 effects on woody plant mass, form, and physiology. Oecologia, 113, 299-313.

Fang Y, Singh B, Singh BP, Krull E (2014a) Biochar carbon stability in four contrasting soils. European Journal of Soil Science, 65, 60-71.

Fang Y, Singh BP, Singh B (2014b) Temperature sensitivity of biochar and native carbon mineralisation in biochar-amended soils. Agriculture, Ecosystems & Environment, 191, 158-167. Farrell M, Kuhn TK, Macdonald LM et al. (2013) Microbial utilisation of biochar-derived carbon. Science of the Total Environment, 465, 288-297.

Foereid B, Lehmann J, Major J (2011) Modeling black carbon degradation and movement in soil. Plant and Soil, 345, 223-236.

Forbes MS, Raison RJ, Skjemstad JO (2006) Formation, transformation and transport of black carbon (charcoal) in terrestrial and aquatic ecosystems. Science of the Total Environment, 370, 190-206.

Glaser B, Balashov E, Haumaier L, Guggenberger G, Zech W (2000) Black carbon in density fractions of anthropogenic soils of the Brazilian Amazon region. Organic Geochemistry, 31, 669678.

Glaser B, Haumaier L, Guggenberger G, Zech W (1998) Black carbon in soils: the use of benzenecarboxylic acids as specific markers. Organic geochemistry, 29, 811-819. Glaser B, Knorr KH (2008) Isotopic evidence for condensed aromatics from non-pyrogenic sources in soils-implications for current methods for quantifying soil black carbon. Rapid Communications in Mass Spectrometry, 22, 935-942.

Glaser B, Lehmann J, Zech W (2002) Ameliorating physical and chemical properties of highly weathered soils in the tropics with charcoal-a review. Biology and Fertility of Soils, 35, 219-230. Gurevitch J, Hedges LV (1999) Statistical issues in ecological meta-analyses. Ecology, 80, 11421149.

Gurwick NP, Moore LA, Kelly C, Elias P (2013) A systematic review of biochar research, with a focus on its stability in situ and its promise as a climate mitigation strategy. PLoS one, 8, e75932.

Hamer U, Marschner B, Brodowski S, Amelung W (2004) Interactive priming of black carbon and glucose mineralisation. Organic Geochemistry, 35, 823-830.

Hammes K, Torn MS, Lapenas AG, Schmidt MW (2008) Centennial black carbon turnover observed in a Russian steppe soil. Biogeosciences, 5, 1339-1350.

Hedges LV, Gurevitch J, Curtis PS (1999) The meta-analysis of response ratios in experimental ecology. Ecology, 80, 1150-1156.

Herath HMSK, Camps-Arbestain M, Hedley MJ, Kirschbaum MUF, Wang T, Hale R (2014) Experimental evidence for sequestering C with biochar by avoidance of CO2 emissions from original feedstock and protection of native soil organic matter. GCB Bioenergy. doi: 10.1111/gcbb.12183.

Hilscher A, Heister K, Siewert C, Knicker H (2009) Mineralisation and structural changes during the initial phase of microbial degradation of pyrogenic plant residues in soil. Organic Geochemistry, 40, 332-342.

Hilscher A, Knicker H (2011) Degradation of grass-derived pyrogenic organic material, transport of the residues within a soil column and distribution in soil organic matter fractions during a 28 month microcosm experiment. Organic Geochemistry, 42, 42-54.

Hu L, Li X, Liu B, Gu M, Dai J (2009) Organic structure and possible origin of ancient charred paddies at Chuodun Site in southern China. Science in China Series D: Earth Sciences, 52, 93100.

Jeffery S, Verheijen FGA, Van Der Velde M, Bastos AC (2011) A quantitative review of the effects of biochar application to soils on crop productivity using meta-analysis. Agriculture, Ecosystems & Environment, 144, 175-187.

Jones DL, Murphy DV, Khalid M, Ahmad W, Edwards-Jones G, DeLuca TH (2011) Short-term biochar-induced increase in soil CO2 release is both biotically and abiotically mediated. Soil Biology and Biochemistry, 43, 1723-1731.

Kasozi GN, Zimmerman AR, Nkedi-Kizza P, Gao B (2010) Catechol and humic acid sorption onto a range of laboratory-produced black carbons (biochars). Environmental Science & Technology, 44, 6189-6195.

Keith A, Singh B, Singh BP (2011) Interactive priming of biochar and labile organic matter mineralization in a smectite-rich soil. Environmental Science & Technology, 45, 9611-9618. Knicker H (2007) How does fire affect the nature and stability of soil organic nitrogen and carbon? A review. Biogeochemistry, 85, 91-118.

Knicker H (2010) "Black nitrogen"-an important fraction in determining the recalcitrance of charcoal. Organic Geochemistry, 41, 947-950.

Knicker H (2011) Pyrogenic organic matter in soil: its origin and occurrence, its chemistry and survival in soil environments. Quaternary International, 243, 251-263.

Kuzyakov Y, Blagodatskaya E (2015) Microbial hotspots and hot moments in soil: Concept & review. Soil Biology and Biochemistry, 83, 184-199.

Kuzyakov Y, Bogomolova I, Glaser B (2014) Biochar stability in soil: Decomposition during eight years and transformation as assessed by compound-specific 14C analysis. Soil Biology and Biochemistry, 70, 229-236.

Kuzyakov Y, Subbotina I, Chen H, Bogomolova I, Xu X (2009) Black carbon decomposition and incorporation into soil microbial biomass estimated by 14C labeling. Soil Biology and Biochemistry, 41, 210-219.

Lehmann J (2007) A handful of carbon. Nature, 447, 143-144.

Lehmann J, Czimczik C, Laird D, Sohi S (2009) Stability of biochar in soil. In: Lehmann J, Joseph S (Eds.), Biochar for Environmental Management: Science and Technology. PP.183-205, Earthscan, London.

Lehmann J, Gaunt J, Rondon M (2006) Bio-char sequestration in terrestrial ecosystems-a review. Mitigation and Adaptation Strategies for Global Change, 11, 395-419. Lehmann J, Rillig MC, Thies J, Masiello CA, Hockaday WC, Crowley D (2011) Biochar effects on soil biota-a review. Soil Biology and Biochemistry, 43, 1812-1836.

Leifeld J, Fenner S, Müller M (2007) Mobility of black carbon in drained peatland soils. Biogeosciences, 4, 425-432.

Liang B, Lehmann J, Sohi SP et al. (2010). Black carbon affects the cycling of non-black carbon in soil. Organic Geochemistry, 41, 206-213.

Liang B, Lehmann J, Solomon D et al. (2008) Stability of biomass-derived black carbon in soils. Geochimica et Cosmochimica Acta, 72, 6069-6078.

Lorenz K, Lal R (2014) Biochar application to soil for climate change mitigation by soil organic

carbon sequestration. Journal of Plant Nutrition and Soil Science, 177, 651-670.

Lützow MV, Kögel-Knabner I, Ekschmitt K, Matzner E, Guggenberger G, Marschner B, Flessa H

(2006) Stabilization of organic matter in temperate soils: mechanisms and their relevance under

different soil conditions-a review. European Journal of Soil Science, 57, 426-445.

Maestrini B, Abiven S, Singh N, Bird J, Torn MS, Schmidt MWI (2014a) Carbon losses from

pyrolysed and original wood in a forest soil under natural and increased N deposition.

Biogeosciences, 11, 5199-5213.

Maestrini B, Nannipieri P, Abiven S (2014b) A meta-analysis on pyrogenic organic matter induced priming effect. GCB Bioenergy. doi: 10.1111/gcbb.12194.

Major J, Lehmann J, Rondon M, Goodale C (2010) Fate of soil-applied black carbon: downward migration, leaching and soil respiration. Global Change Biology, 16, 1366-1379. Malghani S, Gleixner G, Trumbore SE (2013) Chars produced by slow pyrolysis and hydrothermal carbonization vary in carbon sequestration potential and greenhouse gases emissions. Soil Biology and Biochemistry, 62, 137-146.

Marschner B, Brodowski S, Dreves A et al. (2008) How relevant is recalcitrance for the stabilization of organic matter in soils?. Journal of Plant Nutrition and Soil Science, 171, 91-110. Masiello CA, Druffel ERM (1998) Black carbon in deep-sea sediments. Science, 280, 1911-1913. Nguyen BT, Lehmann J (2009) Black carbon decomposition under varying water regimes. Organic Geochemistry, 40, 846-853.

Nguyen BT, Lehmann J, Hockaday WC, Joseph S, Masiello CA (2010) Temperature sensitivity of black carbon decomposition and oxidation. Environmental Science & Technology, 44, 33243331.

Ohlson M, Dahlberg B, 0kland T, Brown KJ, Halvorsen R (2009) The charcoal carbon pool in boreal forest soils. Nature Geoscience, 2, 692-695.

Preston CM, Schmidt MWI (2006) Black (pyrogenic) carbon: a synthesis of current knowledge and uncertainties with special consideration of boreal regions. Biogeosciences, 3, 397-420. Rumpel C, Alexis M, Chabbi A, Chaplot V, Rasse DP, Valentin C, Mariotti A (2006) Black carbon contribution to soil organic matter composition in tropical sloping land under slash and burn agriculture. Geoderma, 130, 35-46.

Santos F, Torn MS, Bird JA (2012) Biological degradation of pyrogenic organic matter in

temperate forest soils. Soil Biology and Biochemistry, 51, 115-124.

Schmidt MW (2004) Biogeochemistry: Carbon budget in the black. Nature, 427, 305-307.

Schmidt MW, Noack AG (2000) Black carbon in soils and sediments: analysis, distribution,

implications, and current challenges. Global Biogeochemical Cycles, 14, 777-793.

Schmidt MW, Skjemstad JO, Jäger C (2002) Carbon isotope geochemistry and nanomorphology

of soil black carbon: Black chernozemic soils in central Europe originate from ancient biomass

burning. Global Biogeochemical Cycles, 16, 1123. doi:10.1029/2002GB001939.

Schmidt MW, Torn MS, Abiven S et al. (2011) Persistence of soil organic matter as an ecosystem

property. Nature, 478, 49-56.

Schneider MP, Smittenberg RH, Dittmar T, Schmidt MW (2011) Comparison of gas with liquid chromatography for the determination of benzenepolycarboxylic acids as molecular tracers of black carbon. Organic Geochemistry, 42, 275-282.

Singh BP, Cowie AL (2014) Long-term influence of biochar on native organic carbon mineralisation in a low-carbon clayey soil. Scientific Reports, 4, 3687.

Singh BP, Cowie AL, Smernik RJ (2012) Biochar carbon stability in a clayey soil as a function of feedstock and pyrolysis temperature. Environmental Science & Technology, 46, 11770-11778. Singh N, Abiven S, Maestrini B, Bird JA, Torn MS, Schmidt MW (2014) Transformation and stabilization of pyrogenic organic matter in a temperate forest field experiment. Global Change Biology, 20, 1629-1642.

Singh N, Abiven S, Torn MS, Schmidt MWI (2012) Fire-derived organic carbon in soil turns over on a centennial scale. Biogeosciences, 9, 2847-2857.

Six J, Conant RT, Paul EA, Paustian K (2002) Stabilization mechanisms of soil organic matter: implications for C-saturation of soils. Plant and Soil, 241, 155-176.

Skjernstad JO, Taylor JA, Smernik RJ (1999) Estimation of charcoal (char) in soils. Communications in Soil Science & Plant Analysis, 30, 2283-2298.

Smith P, Bustamante M, Ahammad H et al. (2014) Agriculture, Forestry and Other Land Use (AFOLU). In: Climate Change 2014: Mitigation of Climate Change. Contribution of Working Group III to the Fifth Assessment Report of the Intergovernmental Panel on Climate Change (eds Edenhofer O, Pichs-Madruga R, Sokona Y, et al.), PP. 829-833. Cambridge University Press, Cambridge, United Kingdom and New York, NY, USA.

Spokas KA (2010) Review of the stability of biochar in soils: predictability of O: C molar ratios. Carbon Management, 1, 289-303.

Spokas KA, Novak JM, Masiello CA, Johnson MG, Colosky EC, Ippolito JA, Trigo C (2014) Physical disintegration of biochar: An overlooked process. Environmental Science & Technology Letters, 1, 326-332.

Spokas KA, Novak JM, Stewart CE, Cantrell KB, Uchimiya M, DuSaire MG, Ro KS (2011) Qualitative analysis of volatile organic compounds on biochar. Chemosphere, 85, 869-882. Steinbeiss S, Gleixner G, Antonietti M (2009) Effect of biochar amendment on soil carbon balance and soil microbial activity. Soil Biology and Biochemistry, 41, 1301-1310. Stewart CE, Zheng J, Botte J, Cotrufo MF (2013) Co-generated fast pyrolysis biochar mitigates greenhouse gas emissions and increases carbon sequestration in temperate soils. GCB Bioenergy, 5, 153-164.

Whitman T, Enders A, Lehmann J (2014) Pyrogenic carbon additions to soil counteract positive priming of soil carbon mineralization by plants. Soil Biology and Biochemistry, 73, 33-41. Wiednera K, Glaser B (2013) Biochar-Fungi Interactions in Soils. In: Biochar and Soil Biota (eds Ladygina N, Rineau F), PP. 69-99.CRC Press, Boca Raton, FL.

Woolf D, Amonette JE, Street-Perrott FA, Lehmann J, Joseph S (2010) Sustainable biochar to mitigate global climate change. Nature communications, 1, 56. doi: 10.1038 / ncomms1053. Yin YF, He XH, Gao R, Ma HL, Yang YS (2014) Effects of rice straw and its biochar addition on soil labile carbon and soil organic carbon. Journal of Integrative Agriculture, 13, 491-498. Zhang M, Ok YS (2014) Biochar soil amendment for sustainable agriculture with carbon and contaminant sequestration. Carbon Management, 5, 255-257.

Zimmerman AR (2010) Abiotic and microbial oxidation of laboratory-produced black carbon (biochar). Environmental Science & Technology, 44, 1295-1301.

Zimmerman AR, Gao B, Ahn MY (2011) Positive and negative carbon mineralization priming effects among a variety of biochar-amended soils. Soil Biology and Biochemistry, 43, 1169-1179. Zimmermann M, Bird MI, Wurster C et al. (2012) Rapid degradation of pyrogenic carbon. Global Change Biology, 18, 3306-3316.

Table 1 Statistical results of comparisons for biochar decomposition rate (% day-1) by Kruskal-Wallis ANOVA and for the weighted response ratios of the priming effect response to biochar addition by Chi-Square test among groups for selected categorical variables.

Biochar decomposition rate

Priming effect

Categorical variable Group n. obs Kruskal-Wallis test n. obs Chi-Square test

Type of experiment Incubation 121 0.59 113 0.37

Field 7 3

Experimental duration <0.5 55 ** * 37.85 59 7.07

(year) 0.5-1 8 5

1-2 46 37

>3 19 15

Feedstock Crop 51 ** * 32.56 48 * ** 31.99

Wood 44 43

Grass 28 20

Manure 4 4

Pyrolysis type Slow 108 0.26 96 * ** 14.21

Fast 20 20

Pyrolysis temperature 200-375 22 ** * 34.53 17 ** 14.10

(°C) 400-475 41 33

500-550 48 56

600-700 11 6

1200 6 4

Ash content <20 21 0.06 22 * 6.37*

(%) >20 22 20

Application amount 0.1-1 50 3.62 38 * * 10.20

(%) 1.1-3 50 54

5-20 28 24

Land use Agriculture 103 1.16 96 0.71

Forest 14 11

Other 11 9

Soil clay content <10 13 ** * 20.42 12 * ** 27.14

(%) 10-20 26 16

20-40 49 48

40-70 27 27

Values in bold indicate significant levels at P<0.05, P<0.01 and P<0.001.

Table 2

Kinetic parameters of the double first-order exponential decay model describing biochar decomposition in soils. Values represent means ± standard errors.

Mean residence

Size Decomposition rate time

Labile C pool 3±0.6% 0.0093% day-1 108±196 days

Recalcitrant C pool 97±0.6% 0.0018% year-1 556±483 years

Figure captions

Fig. 1 Relationships between the decomposed amount (% of added, left) and rates of biochar (% day-1, right) and experimental duration (days). Data are grouped based on biochar feedstocks. The dotted line indicates the 95% confidence band.

Fig. 2 Comparisons of biochar decomposition rates among different groups of the significant explanatory variables based on Table 1. Within each panel, significant difference among groups are indicated by different lower case letters (Kruskal-Wallis ANOVA, P<0.05). The black solid and red dashed lines, lower and upper edges, bars and black circles represent the median and mean values, 25th and 75th, 10th and 90th percentiles, and outliers of all data, respectively. Number of studies included from each group is inserted next to x-axis.

Fig. 3 The percent of biochar remaining in soil plotted against experimental duration. Line represents fit of double first-order exponential decay model to the experimental data. The dotted line indicates the 95% confidence band.

Fig. 4 Comparisons of percent change of priming effects response to biochar addition among different groups of the explanatory variables of experimental duration, feedstock, pyrolysis type and temperature, application amount, and ash content of biochar, and soil clay content. Means and 95% CI are given. Number of studies included from each group is given in parentheses on the left y-axis.

Fig. 5 Overview of factors and processes affecting biochar decomposition in soil. Internal (biochar) and external (environment/soil) factors are included (left). Abiotic, biotic and indirect processes are presented (right). The plus (+) and minus (-) signs indicate the positive and negative effects of different processes on biochar decomposition. See text for further explanation.

Fig. 6 Schematics for potential contributions of biochar addition to soil C sequestration. The solid and dashed green lines indicate the equilibrium and steady state of C in soils after organic matter input, respectively. The black line indicates the possible curve of C equilibrium following biochar addition. Modification from Six et al. (2002) and Lehmann et al. (2006).

Supporting Information captions

Table S1 Characteristics of studies included in this meta-analysis.

O C<1»

O ÙIH

V U*ru*

□ KVscd

¿-oie

»5 7S0 10&5 1460 1825 2190 2555 2920

3 DC 5p

¡5 0 00

¿=022

0 »5 730 1095' 460 1825 2190 2555 2«20

Experimental duration (days)

Experimental duration (year)

0 04 -

0.00 -

n-12 3 7 3 _I_I_I_L_

Feedstock

n-11 8

Pyroiysi» temperature ("C)

• •

n=e 13 10 •» 2

_i_i_i_I_i_

Soil clay content (%)

«•5 9 11 4

__I_I—

o>-g .£ "o c T3 S 18 ¡o

w — ■

P,E Q g

0 ° ...... 0

o 8 □ ..........8

SB ? r2 • 0.16

o° □ A

O Crop

- O»»

o V Manure

o A MM

_ _1_1_ _i_i □ Wood ■ '

0 365 730 1095 1460 1825 2190 2555 2920 Experimental duration (days)

Overall |2t I

«OS 112 >0&-<1 (2 >1—c2 (6 »3 (2

Crop (9 Wood (8 Grast(S Manure (1

Slaw I IB F«»| (3

200-375 (5 400-476 111 500-360 (9 800-700 (3 1200(2

<20(5 >20(3

01-1(12 1 1-3(6 5-20(4

«10(5 »10-<20 (7 >20-<40 (8 >40-<70 (4

r v >'. I ■ . «

Iy«H) --

/Vayui ifjm

vyay*» --ritiivwaiw rc:

Ml eMail! ,'HI

«Hinni

Son dar caiitM (%)

-60 -40 -20 0 20 40 Priming effects <% change)

Factors

Internal Biochar

• Feedstock

• Pyrolysis temperature

• Slow or fast pyrolysis

External Environment/soil

• Abiotic

Temperature Moisture • Texture o Clay Ä minerals SOC

• Biotic . Fi

Biochar

ungi/bactena Saprophytic fungi Soil animals

Afriobc Abiotic oxidation Carbonate dissolution Chem ¡sorption Org/min interactions Physical disintegration Biotic Microbial degradation Pnming effect of biochar Toxic effects on microbes

If. 0 v .00 %

Surface runoi Bioturbation Particulate transport Freerin^thawing Swelling/shrinking

^ / • Decomposition <« input

¿Jff' / • Negative pomiog effect on SOC

/' • Suppression of decompos*ion of other organies

/ • Improvement of soil fertility and crop growth

< • increased belowground C input ' Steady state "

'decomposition ■ input—--— "

C input