# Persistence in soil of Miscanthus biochar in laboratory and field conditions

Daniel P Rasse, Alice Budai, Adam O 'Toole, Xingzhu Ma, Cornélia Rumpel, Samuel Abiven

► To cite this version:

Daniel P Rasse, Alice Budai, Adam O 'Toole, Xingzhu Ma, Cornélia Rumpel, et al.. Persistence in soil of Miscanthus biochar in laboratory and field conditions. PLoS ONE, 2017, 12 (9), pp.e0184383. 10.1371/journal.pone.0184383. hal-01593642

HAL Id: hal-01593642

https://hal.sorbonne-universite.fr/hal-01593642

Submitted on 26 Sep 2017

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_2.png" /></figure>

RESEARCH ARTICLE

# Persistence in soil of *Miscanthus* biochar in laboratory and field conditions

Daniel P. Rasse¹*, Alice Budai¹,², Adam O'Toole¹, Xingzhu Ma³, Cornelia Rumpel⁴, Samuel Abiven⁵

¹ Department of Soil Quality and Climate Change, Norwegian Institute of Bioeconomy Research, Ås, Norway, ² Department of Environmental Sciences, Norwegian University of Life Sciences, Ås, Norway,

³ Institute of Soil Fertilizer and Environment Resource, Heilongjiang Academy of Agricultural Sciences, Harbin, China, ⁴ CNRS, IEES, UMR CNRS-INRA-UPMC-UPEC-IRD-ParisAgroTech, Thiverval-Grignon, France, ⁵ Department of Geography, University of Zurich, Zurich, Switzerland

* daniel.rasse@nibio.no

## Abstract

Evaluating biochars for their persistence in soil under field conditions is an important step towards their implementation for carbon sequestration. Current evaluations might be biased because the vast majority of studies are short-term laboratory incubations of biochars produced in laboratory-scale pyrolyzers. Here our objective was to investigate the stability of a biochar produced with a medium-scale pyrolyzer, first through laboratory characterization and stability tests and then through field experiment. We also aimed at relating properties of this medium-scale biochar to that of a laboratory-made biochar with the same feedstock. Biochars were made of *Miscanthus* biomass for isotopic C-tracing purposes and produced at temperatures between 600 and 700°C. The aromaticity and degree of condensation of aromatic rings of the medium-scale biochar was high, as was its resistance to chemical oxidation. In a 90-day laboratory incubation, cumulative mineralization was 0.1% for the medium-scale biochar vs. 45% for the *Miscanthus* feedstock, pointing to the absence of labile C pool in the biochar. These stability results were very close to those obtained for biochar produced at laboratory-scale, suggesting that upscaling from laboratory to medium-scale pyrolyzers had little effect on biochar stability. In the field, the medium-scale biochar applied at up to 25 t C ha⁻¹ decomposed at an estimated 0.8% per year. In conclusion, our biochar scored high on stability indices in the laboratory and displayed a mean residence time > 100 years in the field, which is the threshold for permanent removal in C sequestration projects.

## Introduction

Progress towards implementing biochar as a technology for biological carbon capture and storage is being made on several fronts. A recent analysis indicates that biochar is on average a more favorable option than other negative emission technologies in terms of required land surface, water use, soil nutrient budgets, energy requirements and costs [1]. Early fears and uncertainty about the impact of large-scale biochar deployment have been tempered by

<figure><img src="image_3.png" /></figure>

### OPEN ACCESS

**Citation:** Rasse DP, Budai A, O'Toole A, Ma X, Rumpel C, Abiven S (2017) Persistence in soil of *Miscanthus* biochar in laboratory and field conditions. PLoS ONE 12(9): e0184383. https://doi.org/10.1371/journal.pone.0184383

**Editor:** Jorge Paz-Ferreiro, RMIT University, AUSTRALIA

**Received:** June 23, 2017

**Accepted:** August 22, 2017

**Published:** September 5, 2017

**Copyright:** © 2017 Rasse et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.

**Data Availability Statement:** All datasets for laboratory and field experiments are available from the PANGAEA database (https://doi.pangaea.de/10.1594/PANGAEA.878894).

**Funding:** Funding for the research was provided by the Research Council of Norway (NFR) through the project "Advanced Techniques to Evaluate the Long-term Stability and Carbon Sequestration Potential of Different Types of Biochar" NFR197531 and "Creating a scientific basis for an integrated evaluation of soil-borne GHG emissions in Norwegian agriculture" NFR/192856; by the
Norwegian Financial Mechanism with Hungary "Green Industry Innovation" project HU09-0029-A1-2013; and by the Norwegian Ministry of Climate and Environment through the NIBIO SIS-Jordkarbon project. Alice Budai received a travel grant (EG/3958) from the MOLTER networking programme of the European Science Foundation for part of this work on BPCA at the University of Zurich. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.

**Competing interests:** The authors have declared that no competing interests exist.

**Abbreviations:** BPCA, benzene polycarboxylic acid; HTT, highest treatment temperature; MRT, mean residence time; SOM, soil organic matter; TOC, total organic carbon.

extensive work to assess possible negative impacts and tradeoffs [2–4] and the creation of industry certification protocols to ensure sustainable production of safe biochar e.g. European Biochar Certificate [5] and biochar standards of the International Biochar Initiative [6]. For C-credit accounting, biochar potentially presents the considerable advantage as compared to other soil carbon sequestration methods of relying on direct C-input accounting rather than expensive soil-based verification schemes [7]. However, C-input accounting is conditional to having an accurate estimator of the mean residence time (MRT) in soil of any given biochar source.

Persistence in soil is a fundamental quality of biochars for serving their role as C sequestration products. This persistence must exceed 100 years to match the definition of permanent removal, as defined by Noble and colleagues [8]. The bulk of plant residue biomass decomposes quickly when applied to soil, with even lignin molecules mineralizing at 90% within one year of residue application to soil [9]. The mean residence time of bulk soil organic matter (SOM) averages 50 years across studies [10]. In other words, biochar must be about 2 orders of magnitude more stable than untreated plant residues and at least twice as stable as bulk SOM to meet a 100-year MRT criteria.

Research on biochar is often carried out using laboratory-produced biochar. Due to limitations of heat transfer and the exothermic nature of pyrolysis, small-scale production offers better control and more sensitive monitoring as compared to larger scale commercial units [11]. The implementation of biochar technology is dependent on the production of biochar through larger scale commercial units. The highest treatment temperature reached during transformation is often different from the target temperature due to the endothermic and exothermic properties of the carbonization process [12], and accurate measurement of temperatures within the reactors are not always possible, especially large-scale ones. This raises the question whether biochar produced in larger reactors is of equivalent quality to that produced in the laboratory using the same feedstock and equivalent temperature.

Up to now, the vast majority of studies aiming at determining the stability of biochar in soils have been laboratory incubations. Reviews of biochar stability in soils have mostly been based on laboratory incubations and on properties of black carbon present in soils exposed to natural fires [13]. In a review of 311 papers, Gurwick and colleagues [14] found only 3 studies estimating biochar stability in the field. Similarly, less than 10% of studies presented in a recent review of biochar effects on soil respiration were based on actual field treatments [15]. Only a subset of these field treatments corresponded to CO₂ field monitoring for at least one growing season. Recently, only three isotopic field studies were available for estimating biochar decomposition and priming effects in soils, while many more came from laboratory conditions [16]. This exemplifies the need for more field evaluation of biochar, especially as its mineralization might be enhanced in field, where active roots are present [17].

One of the problems with laboratory incubations is the fact that they are usually lasting for a few weeks or months and they are addressing the timeframe of 100 years only by extrapolation of the C mineralization data [18]. Field data are needed to improve upon these extrapolations and to calibrate screening methods for biochar stability [19]. Chemical oxidation is such a screening method, which has been proposed to address long-term biochar stability [20, 21]. Another approach is based on the determination of benzene polycarboxylic acids (BPCA) as biomarkers of condensed aromatic sheets, which have been shown to isolate the most stable faction of biochar, and are therefore a promising proxy for stability [22]. Moreover, elemental composition of biochars may also be a proxy for their degradation behavior [23]. Here we considered these three types of proxies for biochar incubated under both laboratory and field conditions.
The overall objective of the present study was to investigate the stability in soils of biochar produced from *Miscanthus* feedstock. The feedstock was chosen because it is a bioenergy crop in Europe and, being a C₄-type grass, its distinct isotopic ¹³C signature can be used to trace the fate of its constitutive carbon in temperate soils. Biochars were produced with slow pyrolysis at different scales, using a medium-scale pyrolyzer (BCMED) and a laboratory unit (BCLAB), and their stabilities were analysed with different laboratory methods and compared to laboratory and field incubation results. The objectives of this study were to: 1) determine if BCMED performed as well as BCLAB in terms of carbonization, condensation, chemical stability indicators and biological stability in laboratory incubation, and 2) estimate the stability of BCMED and its feedstock in a 2-year field experiment.

## Material and methods

### Biochar production and characterization

The *Miscanthus* biochar was produced in Germany in 2010 by Pyreg® Gmbh (www.pyreg.de) in a commercial prototype slow pyrolysis screw reactor operating under a continuous feeding rate of 100–150 kg dry matter per hour and a carbon efficiency of up to 60%. We define this pyrolyzer unit as being of a medium scale and refer to it hereafter as BCMED. The estimated highest treatment temperature (HTT) provided by the manufacturer was between 500–750°C. A precise temperature measurement at each phase of pyrolysis is in general difficult to obtain due to heat transfer limitations and was not possible for this machine. In order to avoid combustion risks, the biochar was moistened to approximately 35% moisture content after leaving the pyrolysis reactor. Application rates in this article were all corrected for moisture and are presented on a dry weight basis. Using the same feedstock as BCMED, we produced slow-pyrolysis biochars under controlled laboratory conditions and obtained a measured HTT of 682°C. This was performed in a muffle furnace with a heating rate of 2.5°C min⁻¹ as described by Budai and colleagues [12]. Hereafter, we will refer to this biochar as BCLAB.

The *Miscanthus* biochar was analyzed for elemental and proximate compositions. Proximate analyses for volatile matter content were conducted according to ASTM E 871 and 872 except that covered crucibles were placed at the rear of a furnace and heated for 6 minutes at 950°C, and ash content was determined according to ASTM D 1102. Specific surface area was measured by N adsorption–desorption isotherms at 77 K using a Micromeritics Tri Star 3000 instrument. Before analysis, the samples were dried at 120°C and degassed overnight in a VacPrep 061 Degasser at 0.05 mbar and 393 K. The Brunauer–Emmet–Teller equation was used to calculate the specific surface area [24]. The C and N contents were determined on a Leco CHN 1000 analyzer (Leco Corp., St. Joseph, MI, USA).

Aromaticity and condensation degree of the *Miscanthus* biochars were estimated with the method of BPCA, following Wiedemeier and colleagues [25]. BPCAs are molecular markers that originate from larger aromatic structures that compose charred biomass. The quantity and composition of the BPCA molecular markers are used to deduce information about the molecular structure of biochar. Here we used total BPCA amount in relation to organic carbon (g kg⁻¹) as an indicator of aromaticity and the ratio of B6CA per total BPCA as an indicator of condensation, as suggested by Wiedemeier and colleagues [26]. The BPCA method was carried out by digesting each ball-milled sample in 10, 15, and 20 mg aliquots for 8 hours at 170°C in quartz tubes using 2 mL of 65% nitric acid solution. The digestate was filtered through ash-free cellulose paper and a cation exchange resin, then finally freeze-dried and re-dissolved in methanol/water (1:1, v:v) before passing through a conditioned solid phase extraction column (Supelco, USA). After drying and re-dissolving in ultrapure water, the final sample was
analyzed using an Agilent 1290 Infinity HPLC system (Santa Clara, USA) according to Wiedemeier and colleagues [25].

Resistance of biochar to oxidation was tested by the acid dichromate method as described by Naisse and colleagues [20], where the total length of time for oxidation is chosen according to the time required to oxidize all of a reference material, i.e. the feedstock in this case. The method of applying fresh potassium dichromate solution and allowing for variable reaction time was applied by Rumpel and colleagues [27] and applied also at room temperature by Kuo and colleagues [28]. Here, samples of 0.3 g each were oxidized in 5 mL of 0.1 M K₂Cr₂O₇ / 2M H₂SO₄ for 1.5 or 2 hours under sonication at 70°C. Samples were recovered by centrifugation and the removal of supernatant, after which oxidation was repeated with new potassium dichromate acid solution. Oxidation was repeated until the *Miscanthus* feedstock was consumed. Total oxidation time was 15.5 hours for all samples. Remaining samples were washed three times with 5 ml distilled water and dried at 60°C for two days. Sample remnants were ground using a mortar and pestle before C and N analysis.

Incubation

Incubation was carried out using a sandy loam Inceptisol collected from an agricultural field in Rygge county, Norway (59°23'15" N; 10°46'26" E) [29]. Soil consisted of 83% sand, 11% silt, and 6% clay (Eurofins AS, Norway), had a pH of 6.8 as measured at a 1:1 soil to water ratio, TOC content of 12 g kg⁻¹ (dw), and a C/N ratio of 12. This soil does not come from our biochar field experiment, but it is a standard soil we used for laboratory incubation of our biochar series [30]. Because of this difference in soil type, our laboratory incubations of BCₘED are not directly comparable to mineralization under field conditions, but rather provide a realistic use of laboratory incubation as a proxy for field stability, where incubation would most likely not be conducted in each soil type where field application is considered. The air-dried soil was passed through a 2 mm sieve, brought to 19.8% (g g⁻¹) moisture content and pre-incubated at 20°C for 20 days. Feedstock and biochars were added to 20 g equivalent dry soil at rates of 0.025, 0.12, 0.58% for feedstock and 0.23, 1.14, 5.46% for biochars. For biochars, these rates mimicked application rates of about 6, 30 and 150 t BC ha⁻¹ within a 0.20 m soil layer of bulk density 1300 kg m⁻³ for the untreated soil. In order to adjust for the faster mineralization of the feedstock, application rate was 10% that of biochar. All mineralization rates were computed based on precise amount added to each vial. The first two rates are close to those used in our field experiments, while the high rate provides an end member for testing potential dose dependent effects on mineralization. Incubation was carried out in 120 mL incubation vials equipped with butyl rubber septa. Determination of the accumulated CO₂ concentration and ¹³CO₂ signature was conducted every 7 to 11 days according to the batch-flush method as described in Budai and colleagues [30]. In short, sampling was conducted by flushing the vial headspace with 800 mL CO₂-free air with outflow gasses collected in 1L gasbags. The concentration and ¹³C isotopic composition in the gasbag was then measured using a cavity ring-down spectrometer (G1101-i, Picarro, INC., Sunnyvale, CA, USA), which had been factory upgraded to reduce transient concentration response, water vapor interference and CH₄ interference according to Moni and colleagues [31]. In addition, a Nafion filter with desiccator was installed on-line to further reduce possible interaction with water vapor [30].

Field trial

A field trial was set up in September 2010 in Ås, Norway (59° 39' 51" N 10° 45' 40" E) in a randomized block design with 4 treatments x 4 blocks. Plots were 8 x 4 m and buffer areas between blocks were 6 m wide. The 4 treatments consisted of: (1) BC<sub>MED</sub> biochar at 8 t C ha<sup>-1</sup>
(BC8), (2) BC<sub>MED</sub> biochar at 25 t C ha<sup>-1</sup> (BC25), (3) *Miscanthus* straw (non-pyrolyzed) at 8 t C ha<sup>-1</sup> (MS8), and (4) control (neither biochar nor non-pyrolyzed *Miscanthus*). Application rates were computed per unit C so that equal quantities of C were added in BC8 and MS8 treatments. Biochar and *Miscanthus* straw were hand spread and raked out on the surface of the plots in September 2010, and immediately incorporated into the soil by inverse ploughing. Inverse ploughing to a depth of 23 cm resulted in the biochar and straw being distributed in concentrated diagonal seams in the Ap horizon in 2011. Ploughing and harrowing after harvest in 2011 and 2012 resulted in a more even distribution throughout the Ap horizon in the following years. Oats were sown in 2011 and barley in 2012. Fertilization was applied with seeding using Yaramila™ NPK 22-3-10 at 550 kg fertilizer ha<sup>-1</sup>. The fields were not treated with fungicide, herbicide or pesticide. Hand weeding was done where weeds appeared within the closed chamber collars.

Annual precipitation in 2011 was 973 mm (63% in May-Sept) and 800 mm in 2012 (47% in May-Sept). Annual average temperature was 6.7°C in 2011 (14 ± 2.6°C, May-Sept), and 5.9°C in 2012 (13.1 ± 2.3°C, May-Sept). These meteorological measurements were taken from a research weather station located on the University of Life Sciences, Ås Campus, 1.3 km from the field site. The soil of the field plots is a clay loam Epistagnic Albeluvisol (WRB classification). The clay content is 27%, silt 43% and sand 31%. pH is 6.39 (±0.18, n = 9), TOC 2.64%, total N 0.23%, and total P 0.29%.

## Field CO₂ monitoring

The CO₂ fluxes were measured during the growing seasons of 2011 and 2012, and isotopic ¹³C composition in 2012 only. Fourteen CO₂ flux measurements were undertaken from 23/05/2011 to 01/09/2011 and 11 measurements from 22/05/2012 to 04/10/2012. Measurements were conducted between 10:00 and 15:00. Thirty-two chamber collars (2 collars/plot x 16 plots) measuring 0.32m L x 0.12m W x 0.06 m H were inserted 0.05m into the soil between crop rows, leaving a water filled gutter (0.1m W x 0.1 m H) exposed at the soil surface to serve as a gas sealant for the chamber. The inter-row chambers capture the soil respiration, including root activities, but exclude the respiratory component of plant shoots, thereby increasing the signal to noise ratio of the isotopic measurements. Thirty-two rectangular aluminum closed chambers (0.30m L x 0.1m W x 0.2m H) were placed on the chamber collars immediately before measurement. There were no pressure valve tubes used on the chambers. The CO₂ flux was measured for 2 min periods with an infrared gas analyzer (IRGA) EGM-4 (PP Systems, Hitchin, UK) which cycled gases via entry and exit valves from the chamber to calculate changes in CO₂ concentration and the flux.

The δ¹³C signature of the soil CO₂ efflux was measured 6 times in 2012. Samples were taken in partially inflated 1-L gas bags [31]. Because the air inside the chamber was a mixture of atmospheric air with increasing concentration of soil-emitted CO₂, keeling plots were necessary to estimate the true δ¹³C value of the soil CO₂. The keeling plot method [32] is used to differentiate δ¹³C SOM from atmospheric δ¹³C where the linear regression plot intercept represents the δ¹³C SOM. The Keeling plot method is based on a linear relationship between δ¹³C values and the inverse of CO₂ concentrations, it is therefore not time dependent, making it a robust method even if release of soil air in the chamber might have been slightly accelerated at sampling. Preliminary tests indicated that 3-point keeling plots with sampling at 3, 8 and 1440 minutes were linear and suitable for covering a wide range of concentration necessary for proper estimates. In our analyses, any keeling plot that did not reach a significant correlation coefficient at P<0.1 (r ≥ 0.988) was excluded. Gas samples were analyzed for δ¹³C using a cavity ring down spectrometer (G1121-i, Picarro INC., Sunnyvale, CA, USA). Solid sample δ¹³C
analysis was carried out on the *Miscanthus* straw, biochar, and the C3 field soil by combusting 1–2 mg samples (3 replicates) in a combustion module connected to a cavity ring-down spectrometer (G2121-i, Picarro, INC., Sunnyvale, CA, USA). The spectrometer was controlled for drift in δ¹³C signal by including known δ¹³C standards, in this case sucrose with -11.6 ‰ and tyrosine at -23.2 ‰, within the analysis runs.

## Statistics

For the laboratory incubations, in order to determine if SOM decomposition was significantly modified by different types and quantities of biochar and feedstock amendment as compared to a non-amended control soil, we applied one-way ANOVA with the Dunnett's method for multiple comparisons vs. a control group, as implemented in SigmaPlot 12.5. Multiple comparisons effects were conducted after verifying that both normality and normal variance conditions were satisfied. Fitting of incubation data to a first-order kinetics decay model was conducted with SigmaPlot 12.5. For field data, statistical analyses of the total soil respiration and mineralization of feedstock and biochar were conducted by 2-way ANOVA, considering treatment and block effects, using the Holm-Sidak method for multiple comparisons when a main effect was detected and both normality and normal variance conditions were satisfied.

## Results

### Laboratory analysis

The BC<sub>MED</sub> produced at a reported temperature between 500–750°C appeared pyrolyzed to an equivalent extent as compared to our reference 682°C laboratory-scale slow pyrolysis biochar (BC<sub>LAB</sub>). Volatile matter content was 7.4 and 6.4% for BC<sub>MED</sub> and BC<sub>LAB</sub>, respectively (Table 1). Carbon content was 80 and 76% for BC<sub>MED</sub> and BC<sub>LAB</sub>, respectively. The H/C atomic ratio was 0.18 and 0.24 for BC<sub>MED</sub> and BC<sub>LAB</sub>, respectively.

Selective oxidations were conducted until all feedstock C was mineralized by the action of the potassium dichromate. At that time, the fraction of non-oxidized C was 75 and 74% for BC<sub>MED</sub> and BC<sub>LAB</sub>, respectively (Table 1). Total benzene poly-carboxylic acid content, i.e. the sum of B3CA, B4CA, B5CA and B6CA, was 179 and 176 g BPCA-C per kg biochar C for BC<sub>MED</sub> and BC<sub>LAB</sub>, respectively (Table 2), indicating high content of aromatic moieties in both biochars. The feedstock contained no B6CA, which is a molecular marker of condensed polyaromatic sheets. By contrast, BC<sub>MED</sub> contained 136 g B6CA-C per kg biochar C, which was 48% more than BC<sub>LAB</sub>. The ratio of B6CA/BPCA was 0.76 for BC<sub>MED</sub> and 0.53 for BC<sub>LAB</sub>, respectively. Both B6CA as percent of charcoal C [33] and the ratio B6CA/BPCA [26] have been suggested as predictors of aromatic condensation, indicating that BC<sub>MED</sub> was more condensed than BC<sub>LAB</sub>.

### Laboratory incubations

In a 90-day incubation, cumulative mineralization of feedstock, SOM and biochar approximated 45%, 1.4% and 0.12% of initial C, respectively (Table 3). This indicates that both biochar types, i.e. BC<sub>MED</sub> and BC<sub>LAB</sub>, were >300 times more stable than the *Miscanthus* feedstock and >10 times more stable than SOM within the 90-day incubation. Because we were interested in determining if the two biochar types behaved differently, we conducted a 2-way ANOVA biochar × dose for BC<sub>MED</sub> and BC<sub>LAB</sub> only (S1 Table). This analysis showed that cumulative mineralization after 90 days was significantly lower for BC<sub>MED</sub> than for BC<sub>LAB</sub> (P = 0.03, S1 Table), while there was neither significant dose effect nor significant dose × biochar interactions. This result suggests that biochar decomposed in a similar fashion whether applied at
<table><caption>Table 1. Properties of the *Miscanthus* feedstock (MS) and derived slow-pyrolysis biochars from medium-scale pyrolyzer (BC<sub>MED</sub>) and a laboratory unit (BC<sub>LAB</sub>).</caption><thead><tr><th>Property</th><th>Units</th><th>MS</th><th>BC<sub>MED</sub></th><th>BC<sub>LAB</sub></th></tr></thead><tbody><tr><td>Volatile Matter</td><td>%</td><td>78.0</td><td>7.4</td><td>6.4</td></tr><tr><td>Fixed Carbon</td><td>%</td><td>13.5</td><td>81.1</td><td>77.7</td></tr><tr><td>Ash</td><td>%</td><td>8.5</td><td>11.5</td><td>15.9</td></tr><tr><td>C</td><td>%</td><td>47.9</td><td>80.0</td><td>75.6</td></tr><tr><td>H</td><td>%</td><td>6.1</td><td>1.2</td><td>1.5</td></tr><tr><td>N</td><td>%</td><td>0.19</td><td>0.6</td><td>0.6</td></tr><tr><td>O</td><td>%</td><td>51.0</td><td>6.6</td><td>5.0</td></tr><tr><td>H/C (atomic)</td><td>-</td><td>1.51</td><td>0.18</td><td>0.24</td></tr><tr><td>O/C (atomic)</td><td>-</td><td>0.80</td><td>0.06</td><td>0.05</td></tr><tr><td>C recalcitrant to K<sub>2</sub>Cr<sub>2</sub>O<sub>7</sub></td><td>%</td><td>-</td><td>75.4</td><td>74.2</td></tr></tbody></table>

https://doi.org/10.1371/journal.pone.0184383.t001

application rates of 0.23, 1.14 or 5.46% by weight. Feedstock mineralization was also consistent across application rates of 0.03, 0.12 and 0.58% (Table 3). As there was no significant dose effect, we investigated the mineralization kinetics of BC<sub>MED</sub> vs. BC<sub>LAB</sub> averaged across doses (Fig 1). This approach shows that the shape of the mineralization curves of BC<sub>MED</sub> and BC<sub>LAB</sub> were similar, although total mineralization was significantly lower for BC<sub>MED</sub> than for BC<sub>LAB</sub> as mentioned above. Mineralization curves of our two biochars were rather quickly leveling off (Fig 1). In fact, modelling with one-pool first-order kinetics decay, predicted reactive C pool in the BC<sub>MED</sub> biochar to be 0.10%, with the remaining fraction of about 99.9% being totally inert (S1 Fig). Using a two-pool model yielded similar results. Forcing the one-pool first-order kinetics model to reach 100% mineralization yielded a decay rate of 1.25 10<sup>-5</sup> d<sup>-1</sup> or MRT of about 220 years, but the fit to the data was poor (S1 Fig). Therefore, our laboratory incubation simply indicates that BC<sub>MED</sub> is highly stable and extrapolating a precise MRT remains uncertain.

Cumulative mineralization of the indigenous SOM was significantly higher in several feedstock and biochar treatments than in the control soil, which averaged 1.38% at the end of the incubation period (Table 4). Largest difference was observed for *Miscanthus* feedstock applied at the highest gravimetric dose with additional loss of 1.08% SOM as compared to the control soil (Table 4). The BC<sub>MED</sub> applied at ten times this rate, i.e. 5.5%, resulted in a 1% priming of SOM, i.e. an increase in mineralization from 1.4% to 2.4% at the end of the 90-day incubation. All application rates at 5.5% and 1.1% for BC<sub>MED</sub> and BC<sub>LAB</sub> biochars and 0.58 and 0.12% for *Miscanthus* feedstock produced a significant increase in SOM mineralization rate as compared to the control. By contrast, no significant difference as compared to control was observed for the lowest amendment rates of 0.23% for BC<sub>MED</sub> and BC<sub>LAB</sub> biochars and 0.03% for *Miscanthus* feedstock. For the middle and high amendment doses, which produced significant

<table><caption>Table 2. Content of BPCA biomarkers in biochars from medium-scale pyrolyzer (BC<sub>MED</sub>) and laboratory unit (BC<sub>LAB</sub>).</caption><thead><tr><th></th><th>BC<sub>MED</sub></th><th>BC<sub>LAB</sub></th></tr></thead><tbody><tr><td>Total BPCA-C (g kg<sup>-1</sup> C)</td><td>179.1 (±0.7)</td><td>175.7 (±1.5)</td></tr><tr><td>B6CA (%)</td><td>76.0 (±0.7)</td><td>52.6 (±0.5)</td></tr><tr><td>B5CA (%)</td><td>14.0 (±0.1)</td><td>27.9 (±0.2)</td></tr><tr><td>B4CA (%)</td><td>10.0 (±0.5)</td><td>18.5 (±0.4)</td></tr><tr><td>B3CA (%)</td><td>0.0 (±0.3)</td><td>1.0 (±0.3)</td></tr></tbody></table>

https://doi.org/10.1371/journal.pone.0184383.t002
<table><caption>Table 3. Mineralization of C sources at the end of a 90-day incubation of feedstock (MS), BC<sub>MED</sub> and BC<sub>LAB</sub> in soils.</caption><thead><tr><th rowspan="2">C substrate</th><th>Dose</th><th colspan="2">Mineralization</th></tr><tr><th>%</th><th></th><th>%</th></tr></thead><tbody><tr><td>SOM</td><td>-</td><td>1.38</td><td>(±0.02)</td></tr><tr><td rowspan="3">MS</td><td>0.03</td><td>43.1</td><td>(±2.6)</td></tr><tr><td>0.12</td><td>44.3</td><td>(±1.6)</td></tr><tr><td>0.58</td><td>46.4</td><td>(±0.3)</td></tr><tr><td rowspan="3">BC<sub>MED</sub></td><td>0.23</td><td>0.13</td><td>(±0.04)</td></tr><tr><td>1.14</td><td>0.07</td><td>(±0.01)</td></tr><tr><td>5.46</td><td>0.09</td><td>(±0.01)</td></tr><tr><td rowspan="3">BC<sub>LAB</sub></td><td>0.23</td><td>0.14</td><td>(±0.06)</td></tr><tr><td>1.14</td><td>0.14</td><td>(±0.02)</td></tr><tr><td>5.46</td><td>0.12</td><td>(±0.00)</td></tr></tbody></table>

Data are percent losses from initial C input, standard deviations for *n* = 3 provided.

https://doi.org/10.1371/journal.pone.0184383.t003

increases in SOM decomposition, we computed the priming effect and tested for differences among both doses and treatments (Table 5). Both amendment type and dose had significant effects on the cumulative priming rates. The higher amendment dose consistently produced higher priming effects, both for feedstock and biochars (Table 5). Because of a significant amendment × dose interaction (P < 0.01, S2 Table), amendment effects were analyzed within dose. At the middle dose, MS and BC<sub>MED</sub> induced a similar priming effect, which was significantly higher than that of BC<sub>LAB</sub>. At the higher dose, priming effects were in the order MS > BC<sub>MED</sub> > BC<sub>LAB</sub>.

## Field experiment

Cumulated soil CO₂ fluxes over the course of the growing season were not significantly affected by treatment in 2011 or 2012, with non-significant highest value in the control treatment in 2011, and in BC25 treatment in 2012 (Table 6, Fig 2). Crop yields were not significantly modified by biochar treatments in either 2011 or 2012 (S3 Table), suggesting that autotrophic respiration terms were fairly similar. Across treatments, cumulated CO₂ fluxes averaged 214 g C m⁻² from May 23rd to September 1st in 2011, and 288 g C m⁻² from May 22nd to October 4th in 2012 (Fig 2).

For most sampling periods, we obtained highly linear keeling plots for estimating δ¹³C values of the soil CO₂ (e.g. S2 Fig). The average δ¹³C of soil CO₂ efflux in the plots amended with *Miscanthus* feedstock was significantly higher than in the control and biochar plots (Table 6). Neither BC8 nor BC25 displayed δ¹³C values of soil CO₂ significantly different from that of the control, although these values were consistently higher in biochar plots (Table 6). As δ¹³C of BC25 and BC8 were not significantly different from one another, we averaged them per block and compared them to the control. This analysis indicated that the increase in δ¹³C of soil CO₂ in the biochar plots as compared to the control was not significant at *P* < 0.05 but was so at *P*<0.1 (*P* = 0.06, S4 Table).

In 2012, the proportion of the soil CO₂ efflux coming from the *Miscanthus* sources ranged between 15 and 29% for straw, and between 0 and 8% for the biochar (Fig 3). This low contribution of biochar sources appear to mask potential differences between the two dose treatments, with the two curves crossing each other (Fig 3). Although proportions of *Miscanthus*-derived CO₂ varied during the growing season, no clear seasonal trend was observed (Fig 3),
<figure><img src="image_4.png" /><figcaption>**Fig 1. Cumulative mineralization in soil of biochars from medium-scale pyrolyzer (BCMED, open symbols) and a laboratory unit (BCLAB, filled symbols).** Standard errors reported for *n* = 9 (3 replicates for 3 doses). Values are in %, i.e. 100 × mineralized fraction.</figcaption></figure>

https://doi.org/10.1371/journal.pone.0184383.g001

suggesting that a season-average value for the δ¹³C was justified. Because 2012 was the only year with isotopic measurements, we estimated the proportions of *Miscanthus*-derived CO₂ for the 2012 growing season only (Table 6). Combining measured CO₂ fluxes and proportions of *Miscanthus*-derived CO₂, we estimated that our MS8 treatment applied at 800 g C m⁻² lost 67 g C m⁻² to the atmosphere during the 2012 growing season, while biochar treatments lost between 6–8 g C m⁻² during the same period (Table 6). These values correspond to a

<table><caption>**Table 4. Multiple comparison test for difference of means in SOM mineralization (from indigenous C3 source) in feedstock, and biochar amended vials vs. the non-amendment control.**</caption><thead><tr><th>Comparison</th><th>Diff of Means<br>(%)</th></tr></thead><tbody><tr><td>MS @ 0.58% vs. soil</td><td>1.08***</td></tr><tr><td>MS @ 0.12% vs. soil</td><td>0.24***</td></tr><tr><td>MS @ 0.03% vs. soil</td><td>0.06ns</td></tr><tr><td>BCMED @ 5.5% vs. Soil</td><td>1.00***</td></tr><tr><td>BCMED @ 1.1% vs. Soil</td><td>0.26***</td></tr><tr><td>BCMED @ 0.2% vs. Soil</td><td>0.06ns</td></tr><tr><td>BCLAB @ 5.5% vs. Soil</td><td>0.45***</td></tr><tr><td>BCLAB @ 1.1% vs. Soil</td><td>0.12**</td></tr><tr><td>BCLAB @ 0.2% vs. Soil</td><td>0.06ns</td></tr></tbody></table>

Treatments are *Miscanthus* feedstock (MS) and biochars from medium-scale pyrolyzer (BCMED) and laboratory unit (BCLAB) (*n* = 3 for all treatments).

** and *** indicate significant differences at *P* < 0.01 and *P* < 0.001, respectively.

ns = non significant.

https://doi.org/10.1371/journal.pone.0184383.t004
<table><caption>**Table 5.** Priming effect by *Miscanthus* feedstock (MS) and biochars from medium-scale pyrolyzer (BC<sub>MED</sub>) and laboratory unit (BC<sub>LAB</sub>).</caption><thead><tr><th rowspan="2">Amendment</th><th colspan="2">Priming effect (%)</th></tr><tr><th>Middle dose</th><th>High dose</th></tr></thead><tbody><tr><td>MS</td><td>17.1<sup>a,A</sup> (±0.9)</td><td>78.5<sup>a,B</sup> (±2.8)</td></tr><tr><td>BC<sub>MED</sub></td><td>19.0<sup>a,A</sup> (±1.6)</td><td>72.4<sup>b,B</sup> (±1.2)</td></tr><tr><td>BC<sub>LAB</sub></td><td>8.6<sup>b,A</sup> (±0.9)</td><td>32.5<sup>c,B</sup> (±0.7)</td></tr></tbody></table>

Standard errors for *n* = 3 between brackets. Within amendment, different capital letters indicate significant dose effect. Within dose, different small letters indicate significant amendment effect. Significances at *P* < 0.05 for means.

https://doi.org/10.1371/journal.pone.0184383.t005

mineralization of the MS8, BC8 and BC25 by 8.3, 0.8 and 0.3%, respectively. The average mineralization value for *Miscanthus* biochar was therefore 0.5% from May 22 to October 4 in 2012. This value translates into an annual mineralization rate of 0.8%, assuming a Q10 of 2 applied to soil temperature values measured at a depth of 2 cm at the Ås field station.

## Discussion

### Carbonization degree

Laboratory analyses pointed towards equivalent degrees of stability and aromaticity for the medium-scale and the laboratory biochars. The H/C atomic ratio of BC<sub>MED</sub> was slightly lower than that of BC<sub>LAB</sub>, i.e. 0.18 vs 0.24 (Table 1). Similar to our results, Keiluweit and colleagues [34] reported H/C atomic ratio of 0.2 for grass biochar produced at 700°C, but did not test higher HTT. However, 0.2 is not the lowest limit for biochar produced with *Miscanthus*, as Budai and colleagues [12] report H/C atomic ratio of 0.1 for biochar produced in the laboratory at 800°C. Therefore, the H/C atomic ratio suggests that BC<sub>MED</sub> reached a carbonization degree comparable to that of BC<sub>LAB</sub>, i.e. a slow-pyrolysis biochar produced in the laboratory at 682°C.

Our chemical oxidation values were close to those reported for a wheat-derived gasification char, which was resistant at 72% to chemical oxidation by potassium dichromate [20]. This latter study used a methodology similar to ours, only with a slightly shorter reaction time, i.e. 12 vs 15.5 h. In general, oxidation methods reported in the literature follow variable protocols, making it difficult to compare results among individual studies. Oxidation utilizing hydrogen peroxide and thermogravimetric analysis have also been used to estimate biochar stability [35]. Our chemical oxidation data suggest that BC<sub>LAB</sub> and BC<sub>MED</sub> were equally carbonized.

The BPCA analyses suggest that BC<sub>MED</sub> produced at a reported temperature between 500–750°C reached a higher condensation degree than our reference 682°C BC<sub>LAB</sub>. Another *Miscanthus* biochar produced by Pyreg was analyzed by Wiedner and colleagues [36] using the BPCA method. Similar to our findings, they found high levels of B6CA, i.e. 85% B6CA, 10% B5CA, 5% B4CA, 0% B3CA. The degree of condensation of this biochar was reported to be higher than all other materials tested [36]. The total BPCA content of our BC<sub>LAB</sub> and BC<sub>MED</sub> are similar to those obtained for grass biochars prepared at 700–900°C [33, 37]. Our results suggest that the medium-scale pyrolysis process affected the condensation more than the aromatization degree of BC<sub>MED</sub> vs. BC<sub>LAB</sub>.

### Stability in laboratory incubations

Laboratory incubations confirmed the high stability of BC<sub>MED</sub>, which was suggested by H/C ratio, BPCA and chemical oxidation methods. BC<sub>MED</sub> mineralized by only 0.10% after 90 days,
<table><caption>Table 6. Proportion of mineralized C₄ C source (*Miscanthus* feedstock or biochar) for the 2012 growing season based on the cumulated CO₂ flux (9 dates) and average δ¹³C values (6 dates).</caption><thead><tr><th rowspan="2">Treat</th><th>Soil CO₂ flux</th><th>δ¹³C</th><th>C₄ CO₂ in flux</th><th>C₄ CO₂ total</th><th>Mineralized C</th></tr><tr><th>g C m⁻²</th><th>‰</th><th>%</th><th>g C m⁻²</th><th>%</th></tr></thead><tbody><tr><td>C</td><td>279 (±29)</td><td>-28.2<sup>b</sup> (±0.2)</td><td></td><td></td><td></td></tr><tr><td>BC8</td><td>262 (±31)</td><td>-27.8<sup>b</sup> (±0.1)</td><td>2.5<sup>b</sup> (±0.7)</td><td>6.0<sup>b</sup> (±1.3)</td><td>0.75<sup>b</sup> (±0.16)</td></tr><tr><td>BC25</td><td>308 (±25)</td><td>-27.8<sup>b</sup> (±0.2)</td><td>2.5<sup>b</sup> (±1.0)</td><td>8.0<sup>b</sup> (±3.7)</td><td>0.32<sup>b</sup> (±0.15)</td></tr><tr><td>MS8</td><td>303 (±45)</td><td>-24.7<sup>a</sup> (±0.3)</td><td>21.9<sup>a</sup> (±1.8)</td><td>66.8<sup>a</sup> (±12.6)</td><td>8.35<sup>a</sup> (±1.58)</td></tr></tbody></table>

Treatments are control (C), 8 t biochar-C ha⁻¹ (BC8), 25 t biochar-C ha⁻¹ (BC25) and 8 t *Miscanthus*-C ha⁻¹ (MS8). Averages with different superscript letters are significantly different at *P* < 0.05 according to the Holm-Sidak method (*n* = 4 replicated blocks).

https://doi.org/10.1371/journal.pone.0184383.t006

which is consistent with results of Luo and colleagues [38] who observed a 0.16% mineralization of 700°C *Miscanthus* biochars in an 87-day incubation. Lower temperature *Miscanthus* biochars have been reported to display higher mineralization rates, from 0.73% in 87 days for a 350°C biochar [38] to 1.1% in 200 days for a 575°C biochar [39]. Here, we could not estimate a precise MRTs based on our short-time laboratory incubation, but even the most conservative first-order kinetics model suggested it to be longer than 220 years (S1 Fig). Even if a laboratory

<figure><img src="image_5.png" /><figcaption>**Fig 2.** Cumulative soil respiration measured 14 dates in 2011 (a) and 11 dates in 2012 (b). Treatments are control (C), non-pyrolyzed Miscanthus feedstock at 8 t C ha⁻¹ (MS8), biochar 8 t C ha⁻¹ (BC8), and biochar 25 t C ha⁻¹ (BC25). Data are averages of *n* = 4.</figcaption></figure>

https://doi.org/10.1371/journal.pone.0184383.g002
<figure><img src="image_6.png" /><figcaption>Fig 3. Proportion of the soil CO₂ efflux coming from the mineralization of 8 t biochar-C ha⁻¹ (BC8), 25 t biochar-C ha⁻¹ (BC25) and *Miscanthus* straw at 8 t C ha⁻¹ (MS8) (standard errors for n = 4).</figcaption></figure>

https://doi.org/10.1371/journal.pone.0184383.g003

MRT could be obtained it could not be extrapolated to field conditions, notably because incubation conditions are artificial and we used a standard soil type. Living roots can promote biochar mineralization [17] and soil type affects biochar mineralization rates [40]. What the incubations tell us is that BC<sub>MED</sub> is highly stable and therefore worthy of field investigation. Incubations are also useful to compare the decomposition kinetics of different biochars [30]. Here we show that the stability of *Miscanthus* biochar produced in a medium-scale pyrolyzer actually exceeds that of biochar produced at a laboratory scale, which suggests that the large volume of feedstock in the pyrolyzer was not a limitation for obtaining a well carbonized product.

## Mineralization of BC in a two-year field trial

Mineralization rate of BC<sub>MED</sub> in the field approximated 0.5% per growing season (Table 6), which implies that the annual rate is probably lower than 1% for the entire year under the cold-climate conditions prevailing in Norway. We acknowledge that the average 0.5% mineralization rate per growing season is only an estimate. However, we found no obvious source of bias on this estimate and therefore consider it fairly robust. Although our soil respiration fluxes were obtained with a simple manual chamber system, our results appear consistent with literature values. We measured on average a soil CO₂ efflux of about 275 g CO₂-C m⁻² over 4 months in 2012, while the annual soil respiration from all croplands averages 544 g C m⁻² yr⁻¹ [41]. Our soil respiration data appear similar or higher to those compiled for field crops in Sweden, Canada and Russia [42].

For soil respiration alone, the absence of a significant difference between our biochar treatments and the control appears consistent with recent reports. For example, Schimmelpfennig and colleagues [43] report that throughout an 18-month monitoring period, a field having received *Miscanthus* biochar had lower cumulative CO₂ emissions than biochar-free controls.
In a recent meta-analysis, Sagrilo and colleagues [15] indicate that soil CO₂ efflux from biochar treated soils are not significantly higher than from no-biochar controls when the ratio of biochar-C to SOC is lower than 2. Across application dose, these authors report no increase in soil CO₂ efflux with biochar addition when the biochar is produced with a pyrolysis retention time > 30 minutes or at a temperature above 550°C, or when it has a surface area > 50 m² g⁻¹. In addition, none of the 8 field studies included in the review of Sagrilo and colleagues [15] displayed significant higher CO₂ fluxes with biochar addition to soil. These findings suggest that biochar decomposition in the field is slow. However, actual quantification of the decomposition rate is crucial, as there is for example a large difference between a 1% and a 5% biochar decomposition rate, although both are likely to produce non-significant CO₂ responses in the field, being possibly hidden by negative priming effects and root respiration responses. Therefore, isotopic tracing of C sources is needed to estimate the actual biochar mineralization rate in the field [16], as was conducted for one growing season in the present study.

Our biochar mineralization estimates computed from δ¹³C and soil respiration measurements are in the lower range of the limited set of studies having attempted a similar assessment. A mineralization rate of 9% was reported for maize biochar after 245 days [17]. However, biochar in the latter study had an atomic H/C ratio of 0.49, which is higher than our 0.18 value. In Australia, mineralization rates of *Eucalyptus* biochar ranged from 2% to 7% per year depending on soil type and climate [40]. This high mineralization rate might be due to the high H/C ratio of the *Eucalyptus* biochar, i.e. 0.63, which is higher than the H/C threshold of 0.6 for proposed for non-stable biochars [23]. Our results are similar to those of Major and colleagues [44], who reported a biochar mineralization rate of 2.2% over 2 years, i.e. about 1% per year, in tropical conditions, using a biochar made of mango tree wood with H/C atomic ratio of 0.26. Also, Maestrini and colleagues [45] reported an *in situ* annual mineralization rate of 0.5% for pine-wood biochar in a temperate forest soil.

Estimating a MRT from the measured biochar mineralization rate in the field is the most crucial yet most uncertain step for assessing the C-storage potential of different biochar products in soil. Having measured a 2% mineralization for biochar over 12 months in an arenosol, Singh and colleagues [40] applied one-, two- and infinite-pool decomposition models and inferred that the corresponding MRT was comprised between 44 and 1079 years, which clearly exemplifies the large uncertainty associated with converting annual mineralization rates into MRT. Major and colleagues [44] observed a mineralization rate of 2.2% over two years, and extrapolated this value to a MRT of 3200 years using a two-pool model. This long MRT was a result of a three-fold decrease in biochar mineralization rate from year one to year two in their study. Our estimated mineralization rate for the 2012 season was slightly lower than that of Major and colleagues [44], i.e. 0.8 vs. 1.1% per year. However, we cannot apply a two-pool model to our results because we have no indication that such two pools actually existed in our case. Laboratory incubation (Fig 1) did not reveal any significant pool of mineralizable C for BC_MED at the beginning of the incubation. By contrast, the feedstock displayed a pronounced two-pool behavior, with 45% being mineralized in 90 days, which might explain why feedstock mineralization rates in the field in 2012 were fairly low. We used a one-pool model with constant mineralization rate of 0.8% per year, which yields a conservative MRT estimate for BC_MED of 125 years. Although this value barely exceeds the conventional 100-year threshold for permanent removal, large gains in terms of C storage in soil can still be achieved with a pyrolysis process transforming crop residues into biochar with 1% y⁻¹ mineralization rate [46].

In conclusion, our biochar produced in a medium-scale pyrolyzer: 1) scored high on stability indices in the laboratory, 2) had similar to higher stability indices than a laboratory-produced biochar, and 3) mineralized at an estimated 0.8% per year under field conditions. The corresponding MRT for field conditions exceeds 100 years, but is only an extrapolation. Based
on laboratory re-incubations, Spokas [47] argues that field-incorporated biochar might become intrinsically more susceptible to mineralization. Others have argued the opposite, that the real MRT might greatly exceed the projected MRT because biochar is not composed of one or two pools but of a continuum of increasingly recalcitrant fractions [40]. Ascertaining the long-term dynamics of this response calls for long-term monitoring of biochar field experiments having isotopic C tracing possibilities.

**Supporting information**

S1 Fig. Modelling with one-pool first-order kinetics models of the mineralization rate of BCMED in laboratory incubation.

(PDF)

S2 Fig. Example of Keeling plot obtained for the determination of the δ13C value of the soil CO2 in plots having received 8 t Miscanthus-C ha-1 (MS8) or 25 t biochar-C ha-1 (BC25).

(PDF)

S1 Table. Two-way ANOVA factorial analysis of biochar type × dose effects on cumulative mineralization after 90 days.

(PDF)

S2 Table. Two-way ANOVA factorial analysis of biochar type × dose effects on cumulative priming after 90 days.

(PDF)

S3 Table. Crop yields in 2011 and 2012 for control (C), 8 t biochar-C ha⁻¹ (BC8), 25 t biochar-C ha⁻¹ (BC25) and 8 t *Miscanthus*-C ha⁻¹ (MS8).

(PDF)

S4 Table. Analysis of variance for average δ13C of soil CO2 in 2012 in the biochar plots (averaged of BC25 and BC8) vs. control plots.

(PDF)

**Acknowledgments**

The authors wish to thank Raphael Fauches for his help with field measurements.

**Author Contributions**

**Conceptualization:** Daniel P. Rasse.

**Data curation:** Alice Budai, Adam O'Toole.

**Formal analysis:** Daniel P. Rasse, Alice Budai, Adam O'Toole, Xingzhu Ma.

**Funding acquisition:** Daniel P. Rasse.

**Investigation:** Daniel P. Rasse, Alice Budai, Adam O'Toole, Xingzhu Ma.

**Methodology:** Daniel P. Rasse, Alice Budai, Adam O'Toole, Cornelia Rumpel, Samuel Abiven.

**Project administration:** Daniel P. Rasse.

**Supervision:** Daniel P. Rasse, Cornelia Rumpel, Samuel Abiven.

**Writing – original draft:** Daniel P. Rasse, Alice Budai, Adam O'Toole.
**Writing – review & editing:** Daniel P. Rasse, Cornelia Rumpel, Samuel Abiven.

## References

1. Smith P. Soil carbon sequestration and biochar as negative emission technologies. Global Change Biol. 2016; 22(3):1315-24. https://doi.org/10.1111/gcb.13178 PMID: 26732128

2. Hale SE, Lehmann J, Rutherford D, Zimmerman AR, Bachmann RT, Shitumbanuma V, et al. Quantifying the Total and Bioavailable Polycyclic Aromatic Hydrocarbons and Dioxins in Biochars. Environ Sci Technol. 2012; 46(5):2830-8. https://doi.org/10.1021/es203984k PMID: 22321025

3. Cowie AL, Downie AE, George BH, Singh BP, Van Zwieten L, O'Connell D. Is sustainability certification for biochar the answer to environmental risks? Pesquisa Agropecuaria Brasileira. 2012; 47(5):637-48.

4. Jeffery S, Bezemer TM, Cornelissen G, Kuyper TW, Lehmann J, Mommer L, et al. The way forward in biochar research: targeting trade-offs between the potential wins. GCB Bioenergy. 2013; 7(1):1-13. https://doi.org/10.1111/gcbb.12132

5. EBC. 'European Biochar Certificate Guidelines for a Sustainable Production of Biochar.' European Biochar Foundation (EBC), Arbaz, Switzerland. http://www.european-biochar.org/en/download. Version 6.2E of 4th February 2016, 10.13140/RG.2.1.4658.7043. 2016.

6. IBI. Standardized product definition and product testing guidelines for biochar that is used in soil. Avail. Online: http://www.biochar-international.org/characterizationstandard 2015.

7. Lehmann J. A handful of carbon. Nature. 2007; 447(7141):143-4. https://doi.org/10.1038/447143a PMID: 17495905

8. Noble IR, Apps MJ, Houghton RA, Lashof DA, Makundi W, Murdiyarso D, et al. Implications of different definitions and generic issues. In: Watson RT, Noble IR, Bolin B, Ravindranath NH, Verardo DJ, Dokken DJ, editors. Land Use, Land Use Change and Forestry. UK: Cambridge; 2000. p. 55-126.

9. Rasse DP, Dignac MF, Bahri H, Rumpel C, Mariotti A, Chenu C. Lignin turnover in an agricultural field: from plant residues to soil-protected fractions. Eur J Soil Sci. 2006; 57(4):530-8. https://doi.org/10.1111/j.1365-2389.2006.00806.x

10. Schmidt MWI, Torn MS, Abiven S, Dittmar T, Guggenberger G, Janssens IA, et al. Persistence of soil organic matter as an ecosystem property. Nature. 2011; 478(7367):49-56. https://doi.org/10.1038/nature10386 PMID: 21979045

11. Ripberger GD, Jones JR, Paterson A, Holt R. Is it possible to produce biochar at different highest treatment temperatures in the pyrolysis range?—The exothermic nature of pyrolysis [online]. Asia Pacific Confederation of Chemical Engineering Congress 2015: APCCHE 2015, incorporating CHEMECA 2015. 2015:1950-7.

12. Budai A, Wang L, Gronli M, Strand LT, Antal MJ, Abiven S, et al. Surface Properties and Chemical Composition of Corncob and Miscanthus Biochars: Effects of Production Temperature and Method. J Agric Food Chem. 2014; 62(17):3791-9. https://doi.org/10.1021/jf501139f PMID: 24720814

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

14. Gurwick NP, Moore LA, Kelly C, Elias P. A Systematic Review of Biochar Research, with a Focus on Its Stability in situ and Its Promise as a Climate Mitigation Strategy. Plos One. 2013; 8(9). https://doi.org/10.1371/journal.pone.0075932 PMID: 24098746

15. Sagrilo E, Jeffery S, Hoffland E, Kuyper TW. Emission of CO2 from biochar-amended soils and implications for soil organic carbon. GCB Bioenergy. 2014; 7(6):1294-304. https://doi.org/10.1111/gcbb.12234

16. Wang J, Xiong Z, Kuzyakov Y. Biochar stability in soil: meta-analysis of decomposition and priming effects. GCB Bioenergy. 2016; 8(3):512-23. https://doi.org/10.1111/gcbb.12266

17. Ventura M, Alberti G, Viger M, Jenkins JR, Girardin C, Baronti S, et al. Biochar mineralization and priming effect on SOM decomposition in two European short rotation coppices. GCB Bioenergy. 2014; 7(5):1150-60. https://doi.org/10.1111/gcbb.12219

18. Naisse C, Girardin C, Lefevre R, Pozzi A, Maas R, Stark A, et al. Effect of physical weathering on the carbon sequestration potential of biochars and hydrochars in soil. Global Change Biology Bioenergy. 2015; 7(3):488-96. https://doi.org/10.1111/gcbb.12158

19. Brewer CE, Unger R, Schmidt-Rohr K, Brown RC. Criteria to Select Biochars for Field Studies based on Biochar Chemical Properties. Bioenerg Res. 2011; 4(4):312-23. https://doi.org/10.1007/s12155-011-9133-7

20. Naisse C, Alexis M, Plante AF, Wiedner K, Glaser B, Pozzi A, et al. Can biochar and hydrochar stability be assessed with chemical methods? Org Geochem. 2013; 60:40-4. https://doi.org/10.1016/j.orggeochem.2013.04.011
21. Ngo P-T, Rumpel C, Ngo Q-A, Alexis M, Vargas GV, Mora Gil MdLL, et al. Biological and chemical reactivity and phosphorus forms of buffalo manure compost, vermicompost and their mixture with biochar. *Bioresour Technol.* 2013; 148:401–7. http://dx.doi.org/10.1016/j.biortech.2013.08.098. https://doi.org/10.1016/j.biortech.2013.08.098 PMID: 24071441

22. Kuzyakov Y, Bogomolova I, Glaser B. Biochar stability in soil: Decomposition during eight years and transformation as assessed by compound-specific C-14 analysis. Soil Biol Biochem. 2014; 70:229–36. https://doi.org/10.1016/j.soilbio.2013.12.021

23. Schimmelpfennig S, Glaser B. One Step Forward toward Characterization: Some Important Material Properties to Distinguish Biochars. J Environ Qual. 2012; 41(4):1001–13. https://doi.org/10.2134/jeq2011.0146 PMID: 22751042

24. Brunauer S, Emmett PH, Teller E. Adsorption of Gases in Multimolecular Layers. J Am Chem Soc. 1938; 60(2):309–19. https://doi.org/10.1021/ja01269a023

25. Wiedemeier DB, Hilf MD, Smittenberg RH, Haberle SG, Schmidt MWI. Improved assessment of pyrogenic carbon quantity and quality in environmental samples by high-performance liquid chromatography. J Chromatogr. 2013; 1304:246–50. https://doi.org/10.1016/j.chroma.2013.06.012 PMID: 23880465

26. Wiedemeier DB, Abiven S, Hockaday WC, Keiluweit M, Kleber M, Masiello CA, et al. Aromaticity and degree of aromatic condensation of char. Org Geochem. 2015; 78:135–43. https://doi.org/10.1016/j.orggeochem.2014.10.002

27. Rumpel C, Alexis M, Chabbi A, Chaplot V, Rasse DP, Valentin C, et al. Black carbon contribution to soil organic matter composition in tropical sloping land under slash and burn agriculture. Geoderma. 2006; 130(1–2):35–46. https://doi.org/10.1016/j.geoderma.2005.01.007

28. Kuo LJ, Herbert BE, Louchouarn P. Can levoglucosan be used to characterize and quantify char/charcoal black carbon in environmental media? Org Geochem. 2008; 39(10):1466–78. https://doi.org/10.1016/j.orggeochem.2008.04.026

29. O'Toole A, de Zarruk KK, Steffens M, Rasse DP. Characterization, Stability, and Plant Effects of Kiln-Produced Wheat Straw Biochar. J Environ Qual. 2013; 42(2):429–36. https://doi.org/10.2134/jeq2012.0163 PMID: 23673835

30. Budai A, Rasse DP, Lagomarsino A, Lerch TZ, Paruch L. Biochar persistence, priming and microbial responses to pyrolysis temperature series. Biol Fertility Soils. 2016:1–13. https://doi.org/10.1007/s00374-016-1116-6

31. Moni C, Rasse DP. Detection of simulated leaks from geologically stored CO2 with C-13 monitoring. International Journal of Greenhouse Gas Control. 2014; 26:61–8. https://doi.org/10.1016/j.ijggc.2014.04.010

32. Keeling CD. The concentration and isotopic abundances of atmospheric carbon dioxide in rural areas. Geochim Cosmochim Acta. 1958; 13(4):322–34. http://dx.doi.org/10.1016/0016-7037(58)90033-4.

33. Schneider MPW, Smittenberg RH, Dittmar T, Schmidt MWI. Comparison of gas with liquid chromatography for the determination of benzenepolycarboxylic acids as molecular tracers of black carbon. Org Geochem. 2011; 42(3):275–82. https://doi.org/10.1016/j.orggeochem.2011.01.003

34. Keiluweit M, Nico PS, Johnson MG, Kleber M. Dynamic Molecular Structure of Plant Biomass-Derived Black Carbon (Biochar). Environ Sci Technol. 2010; 44(4):1247–53. https://doi.org/10.1021/es9031419 PMID: 20099810

35. Yang F, Zhao L, Gao B, Xu X, Cao X. The Interfacial Behavior between Biochar and Soil Minerals and Its Effect on Biochar Stability. Environ Sci Technol. 2016; 50(5):2264–71. https://doi.org/10.1021/acs.est.5b03656 PMID: 26828311

36. Wiedner K, Naisse C, Rumpel C, Pozzi A, Wieczorek P, Glaser B. Chemical modification of biomass residues during hydrothermal carbonization—What makes the difference, temperature or feedstock? Org Geochem. 2013; 54:91–100. https://doi.org/10.1016/j.orggeochem.2012.10.006

37. Wurster CM, Saiz G, Schneider MPW, Schmidt MWI, Bird MI. Quantifying pyrogenic carbon from thermosequences of wood and grass using hydrogen pyrolysis. Org Geochem. 2013; 62:28–32. https://doi.org/10.1016/j.orggeochem.2013.06.009

38. Luo Y, Durenkamp M, De Nobili M, Lin Q, Brookes PC. Short term soil priming effects and the mineralisation of biochar following its incorporation to soils of different pH. Soil Biol Biochem. 2011; 43(11):2304–14. https://doi.org/10.1016/j.soilbio.2011.07.020

39. Bai M, Wilske B, Buegger F, Esperschutz J, Kammann CI, Eckhardt C, et al. Degradation kinetics of biochar from pyrolysis and hydrothermal carbonization in temperate soils. Plant Soil. 2013; 372(1–2):375–87. https://doi.org/10.1007/s11104-013-1745-6

40. Singh BP, Fang Y, Boersma M, Collins D, Van Zwieten L, Macdonald LM. In Situ Persistence and Migration of Biochar Carbon and Its Impact on Native Carbon Emission in Contrasting Soils under
Managed Temperate Pastures. PLoS ONE. 2015; 10(10):e0141560. https://doi.org/10.1371/journal.pone.0141560 PMID: 26509506

41. Raich JW, Schlesinger WH. The global carbon dioxide flux in soil respiration and its relationship to vegetation and climate. Tellus B. 1992; 44(2):81-99. https://doi.org/10.1034/j.1600-0889.1992.t01-1-00001.x

42. Raich JW, Tufekciogul A. Vegetation and soil respiration: Correlations and controls. Biogeochemistry. 2000; 48(1):71-90. https://doi.org/10.1023/a:1006112000616

43. Schimmelpfennig S, Muller C, Grunhage L, Koch C, Kammann C. Biochar, hydrochar and uncarbonized feedstock application to permanent grassland-Effects on greenhouse gas emissions and plant growth. Agriculture Ecosystems & Environment. 2014; 191:39-52. https://doi.org/10.1016/j.agee.2014.03.027

44. Major J, Lehmann J, Rondon M, Goodale C. Fate of soil-applied black carbon: downward migration, leaching and soil respiration. Global Change Biol. 2010; 16(4):1366-79. https://doi.org/10.1111/j.1365-2486.2009.02044.x

45. Maestrini B, Abiven S, Singh N, Bird J, Torn MS, Schmidt MWI. Carbon losses from pyrolysed and original wood in a forest soil under natural and increased N deposition. Biogeosciences. 2014; 11(18):5199-213. https://doi.org/10.5194/bg-11-5199-2014

46. O'Toole A, Rasse D. Biochar: Soil Carbon and Fertility. Encyclopedia of Soil Science, Third Edition: CRC Press; 2016. p. 193-7.

47. Spokas KA. Impact of biochar field aging on laboratory greenhouse gas production potentials. GCB Bioenergy. 2013; 5(2):165-76. https://doi.org/10.1111/gcbb.12005