Biochar synthesis from mineral and ash-rich waste biomass, part 2: characterization of biochar and co-pyrolysis mechanism for carbon sequestration

Abstract

The increase in mineral and ash-rich waste biomass (MWB) generation in emerging economies poses critical environmental problems and bottlenecks the solid waste and wastewater treatment systems. Transforming these MWB such as sewage sludge from wastewater treatment (SSW) to biochar can be a sustainable method for their disposal and resource recovery. However, such biochar has limited applicability due to the relatively low organic content and possibly contaminated nature of SSW. This may be offset through combined pyrolysis with other MWB, which can also support municipal solid waste management. Studies on this MWB co-pyrolysis are lacking and have not yet seen successful long-term implementation. This work is the second part of authors’ research encompassing an analytical and lab-scale investigation of biochar production from MWB through pyrolysis for the case of Chennai city, India. Here, the physicochemical properties of biochar derived from lab-scale co-pyrolysis of SSW with other MWB such as anaerobic digestate from waste to energy plants of food, kitchen or market waste fermentation, and banana peduncles (BP) collected from vegetable markets and their thermolysis mechanism are comprehensively investigated for purpose of carbon sequestration. Also, a novel preliminary investigation of the effect of sample weight (scaling effect) on the analytical pyrolysis of biomass (BP as model substrate) is undertaken to elucidate its impact on the heat of pyrolysis and carbon distribution in resultant biochar. The maximum carbon sequestration potential of the derived biochar types is 0.22 kg CO₂ kg⁻¹ biomass. The co-pyrolysis of MWB is exothermic and governed by the synergetic effects of the components in blends with emission profiles following the order CO₂ > CH₄ > CO > NH₃. Co-pyrolysis reduced the heavy metal enrichment in SSW biochar. The derived biochars can be an immediate source of N, P and S in nutrient-deficient acidic soils. The biochar has only up to 4-ring polyaromatic compounds and a residence time longer than 1 h at 500 °C is necessary to improve carbonization. The heat released during analytical pyrolysis of the model biomass and distribution of carbon in the resultant biochar are significantly influenced by scaling effects, drawing attention to the need for a more detailed scaling investigation of biomass pyrolysis.

1 Introduction

Rampant urbanization and resource utilization in emerging economies increase the generation of mineral and ash-rich waste biomass (MWB) such as sewage sludge, anaerobic digestate, and refuse lignocellulosic biomasses, which constrains the solid waste treatment (SWT) infrastructure in their urban areas. Global adoption of the waste hierarchy principle (prevention, reuse, recycling, and recovery) is essential to meet the UN sustainable development goal (SDG) 12.5 – “substantial reduction of waste generation by 2030” [1]. The integration of recycling and energy recovery from waste biomass into the national development strategies of these nations can aid in sustainable growth via a circular economy (CE). Industrialized countries have already made strides to achieve a CE. An example is the CE action plan adopted by the European Union in 2015 [2]. In such a CE development model, it is also necessary to sustainably upcycle the recycled products derived from SWT for utilization as secondary raw materials. Biomass can be valorized as biochar (carbonaceous residue) through slow pyrolysis in an inert or oxygen-deficient environment. This upcycling pathway converts waste biomass into a value-added product with multifaceted environmental applications such as carbon sequestration (UN SDG 13.2), soil amelioration (UN SDG 15.3), and adsorption of pollutants in water (UN SDG 6.3) [3]. Thus, the conversion of MWB into biochar has the potential to support SWT infrastructure in emerging economies while simultaneously tackling these four areas of the 2030 UN sustainable development agenda.

Despite the aforementioned benefits, biochar generation systems are still not widely implemented. The pyrolysis of biomass, particularly MWB, is a complex thermo-chemical process due to the interaction of various alkali and alkaline earth metals (AAEM), transition metal compounds, and other inorganics (e.g., silicates) with the lignocellulosic constituents. The inherent heterogeneity in MWB restricts the generalized modelling of their thermolysis mechanism. Sewage sludge from wastewater treatment (SSW) is an MWB that causes major bottlenecks in municipal wastewater treatment (WWT) in emerging economies. For instance, India produces ~ 3.95 Mt of dry SSW from only 32% of wastewater that is treated before disposal [4]. The commonly utilized SWT methods continue to be open-dumping and landfilling [5], which are also restrained by rising land costs, consumption and pollution trends. Recent attempts at waste-to-energy techniques for such MWB have also incurred significant failures [6]. Hence, there is an urgent need to investigate and design SSW treatment systems in this region that can also aid to mitigate climate change and support nutrient recovery in line with CE principles [7]. Pyrolytic conversion of SSW into biochar can be one such treatment method. Owing to the heterogeneous nature of SSW (even comparing two cities within the same country) and the known constraints [8] associated with SSW pyrolysis, the thermolysis and resultant biochar properties must be analyzed on a case-by-case basis to evaluate their commercial viability.

Co-pyrolyzing SSW with other organic matter can render the SWT of such MWB economically viable [9] and may also improve the quality (soil applicability) of biochar derived from SSW [10, 11]. It is generally advised to mix SSW with other additives/materials to overcome the drawbacks of its pyrolysis. Though co-pyrolysis studies of SSW with agricultural residues are reported in the literature [12, 13], reports on combined pyrolysis of SSW with other MWB are generally lacking. The pyrolytic conversion of SSW is still considered technologically immature with reported implementation failure even in industrialized nations such as Germany [14]. A detailed investigation from the analytical to industrial scale becomes crucial in implementing biochar production from SSW. In the part-1 study [15], for the case of the southern Indian city of Chennai, the authors have investigated the analytical pyrolysis of three types of MWB that are a) rich in ash – digested SSW from a wastewater treatment plant (M100), b) rich in minerals – banana peduncle with high potassium concentration (M010), and c) rich in ash and minerals – anaerobic digestate from a food waste processing facility (M001).

In this part-2 study, biochar is produced through the co-pyrolysis of different blends of M100 with M010 and M001 in a lab-scale reactor and the physicochemical properties are comprehensively investigated and compared concerning their carbon sequestration potential. The thermal degradation mechanism during co-pyrolysis is analyzed through simultaneous thermal analysis (TGA–DSC) coupled with infrared spectroscopy (FT-IR). Then, for the first time, a preliminary investigation is carried out to evaluate the influence of the sample weight (scaling effect) during analytical pyrolysis of biomass since the heat and mass transfer effects can affect the thermolysis mechanism [16, 17] during analytical pyrolysis, which is generally used to model industrial and lab-scale processes.

2 Materials and methods

2.1 Materials

The MWB in this study was collected from Chennai, India. These included a) anaerobic digestate from an organic waste treatment facility collected by Central Leather Research Institute, Chennai (M001). b) banana peduncles (M010) from the Koyambedu vegetable and flower market, and c) anaerobically digested sewage sludge (M100) from the Perungudi wastewater treatment plant. These constituted the unmixed baseline materials (BM) whose pyrolysis kinetics was investigated in the part-1 study [15]. Then, they were oven dried (105 °C for 48 h) and shredded to 0.2 mm particle size distribution. The elemental and proximate analysis of BM is shown in Table S1 of supplementary materials (SM). Seven blends of these three materials were prepared by keeping M100 at 50 wt% and varying the others. Table 1 shows the proportion of these MM, the labelling of different MM blends and the resulting biochar.

Table 1 List of substrates – mixed and baseline materials – used for co-pyrolysis, and the names of the biochar derived from these materials

2.2 Lab-scale slow pyrolysis

About 25 g of each B-BM and B-MM were pyrolyzed in a lab-scale furnace (Fig. S1 of SM) by linear heating (10 °C min⁻¹) from 30 to 500 °C in 0.4–0.5 L min⁻¹ nitrogen flow. At the highest treatment temperature (HTT) of 500 °C, they were kept isothermal for 60 min. During this isothermal regime, the maximum temperature attained by the materials was termed the highest pyrolysis temperature (HPT). The biochar was recovered and stored in sealed containers for transport and analysis. The name of these biochar types includes the prefix ‘B’ as shown in Table 1.

2.3 Biochar characterization

The elemental analyses of C, H, N, and O were carried out by Mikroanalytisches Laboratorium Kolbe GmbH (Fürth, Germany) after overnight drying at 105 °C. In short, biochar samples were weighed and dried overnight at 105 °C. The elements C, H, S, and N were measured in the Mikro Cube CHNS analyzer (Elementar, Frankfurt, Germany). Oxygen was measured in Vario EL Cube (Elementar). The heavy metals (HM), micro (K, Mg, Na and Ca) and macronutrients (P, S) concentration in the biochar were analyzed by inductively coupled plasma optical emission spectrometry (ICP-OES) iCAP 6000 (Thermo Fisher Scientific Corporation, Waltham, Massachusetts, USA). Infrared spectra (4000–400 cm⁻¹) of the biochar were collected using the method of Attenuated Total Reflection (ATR) using Nicolet iS50 (Thermo Fisher Scientific Corporation) having a KBr beam splitter, diamond crystal, and deuterated triglycine sulfate detector. The measurement parameters of optical velocity, signal gain, and aperture size were 0.1581 cm s⁻¹, 2.0, and 87, respectively. The FT-IR spectra were baseline-corrected. The aromatic to aliphatic ratio (A_i) was calculated according to the Eq. (1) where the region 1620–1530 cm⁻¹ represents the νC = C of the aromatic ring, and 1460–1270 cm⁻¹ includes δ_as CH₂, CH₃ and δ_s CH₃ of aliphatic structures [18].

$${\mathrm{A}}_{\mathrm{i}}=\frac{{\mathrm{Area}}_{1620-1530}}{{\mathrm{Area}}_{1460-1270}}$$

(1)

Scanning electron micrographs were captured (12 kV and 15 kV, resolution 4) using Quanta 200 (FEI, Oregon, USA) after degassing in N₂ (for 6 to 8 min) and gold coating (at 25 mA for 2 min) the samples in Edwards S150A sputter coater. The vacuum condition for the electron gun was kept below 8.7 × 10^–6 kPa. Nitrogen physisorption measurements were made at 77 K in NOVA 4000 (Quantachrome Instruments, Boynton Beach, Fl) to calculate the pore surface area (S_BET) through Brunauer–Emmett–Teller (BET) method [19]. Before analysis, each biochar sample was vacuum degassed at 40 °C until the cell pressure approached nearly zero kPa. The thermal oxidative stability of biochar was analyzed using TGA/DSC 3 + /LF (Mettler Toledo, Ohio, USA). For this, 5 mg of each sample was combusted in 70 µL alumina crucibles from 25 to 1050 °C under a synthetic air (purity of 99.999% from Linde Gases GmbH) flow of 70 mL min⁻¹ and a heating rate of 10 °C min⁻¹. Then the thermal stability was evaluated as thermal oxidative recalcitrance (TOR_i) as defined in ref [20] using MATLAB R2020b. M010 is the only naturogenic feedstock among the materials. Hence it was used as the reference for evaluating the TOR_i. The volatile matter (VM), ash content and fixed carbon (FC) (collectively referred to as proximate analysis) were also measured using TGA. At first, samples were pre-purged in an N₂ flow of 40 mL min⁻¹ for 10 min. Residual moisture was removed isothermally at 110 °C for 10 min in an N₂ atmosphere. Then, it was linearly heated (50 °C min⁻¹ with 40 mL min⁻¹ N₂ flow) and kept isothermal at 950 °C for 10 min for the removal of VM. The sample was cooled to 550 °C and kept isothermal for 10 min in synthetic air (100 mL min⁻¹) to evaluate the loss-on-ignition. Finally, the sample was cooled to 110 °C in N₂ flow (40 mL min⁻¹) and kept isothermal for ash content determination. Measurements were made in triplicates and averaged. The total organic carbon (TOC) was measured in vario TOC cube (Elementar) after ~ 10 mg of biochar was treated with 1 M HCl to remove inorganic carbonates. The zeta potential (ζ) was analyzed using ZetaPALS (BrookHaven Instruments Corporation, New York, USA) with 5 runs (10 cycles per run) for each biochar.

The pyrolysis–gas chromatography-mass spectroscopy (Py-GC–MS) was performed in a double shot Pyrolyzer EGA/PY-3030D (Frontier Laboratories Ltd., Fukushima, Japan) connected to a Trace 1310 GC and ISQ 7000 Single quadrupole MS-System (Thermo Fisher Scientific). Unaltered ground samples (~ 25 mg) were placed in small crucibles and introduced into the furnace via autosampler (AS-1020E, F-Lab). The flash pyrolysis was performed at 800 °C and evolved gases were then directly injected into the GC/MS system for analysis. The gas chromatograph was equipped with a low-to-mid polarity Ultra Alloy® capillary column (UA⁺-5, F-Lab.) of 30 m × 250 µm × 0.25 µm film thickness. The carrier gas was helium at 1 mL min⁻¹. In the mass spectrometer, Total Ion Current chromatograms (TIC) were acquired at 70 eV ionizing energy. The acquired chromatograms were analyzed using GC–MS Postrun Analysis v4.52 (Shimadzu, Kyoto, Japan). Each chromatogram was processed with a background spectrum between 18 and 24 min and a smoothing window of 0.05 s. For peak detection and integration, a slope of 5000 min⁻¹ and a minimum height of 20 × 10⁶ units were used. The peaks were identified using the NIST database (search depth = 3, with reverse search enabled and minimum similarity index > 90).

The properties of water extractable fractions (WEF) such as pH, electrical conductivity (EC), and dissolved organic matter (DOM) were analyzed after extraction using a 0.45 μm syringe filter (CHROMAFIL® Xtra Hydrophilic PES membrane). In short, ~ 100 mg of biochar was mixed with 10 mL of ultrapure water. Then, the solution was stirred (150 to 200 rpm) for 3 h using a magnetic stirrer. After stirring, the solutions were stored in a dark container for 24 h. Then, the suspension was filtered (0.45 μm), and the pH and EC of the supernatant were measured with a multimeter WTW MultiLine® Multi 3630 IDS SET F/ SenTix® 940/ TetraCon® 925/ FDO® 925. Dissolved organic carbon (DOC) was measured as non-purgeable organic carbon using Vario TOC cube (Elementar) with a measurement uncertainty of 2%. The UV–VIS spectra of supernatant in single-use cuvettes (layer thickness, l, = 10 mm) were collected between 200 and 900 nm using DR6000 UV–VIS spectrophotometer (Hach, Colorado, USA). SUVA₂₅₄ (an indicator of DOM aromaticity) and E2/E3 (ratio of UV absorbance at 254 to 365 nm is inversely proportional to molecular weight of DOM) were evaluated. About 10 mL of this supernatant was used to measure the cations (K⁺, Na⁺, Mg²⁺, Ca²⁺, and NH₄⁺) and anions (Cl⁻, SO₄²⁻) in an Ion chromatograph (Thermo Scientific Dionex Aquion/Dionex Integrion HPIC with Dionex AS-AP autosampler).

2.4 Pyrolysis characteristics

The pyrolysis of each MM (Table 1) was investigated using a TGA/DSC 3 + LF thermal analyzer (Mettler Toledo) at six heating rates – (5, 7, 10, 12, 15 and 20 °C min⁻¹) between 25 and 1000 °C using 70 μL alumina crucibles (N₂ flow of 70 mL min⁻¹). Each measurement was preceded by a 10 min pre-purge using 100 mL min⁻¹ of N₂ and was performed in duplicates. The thermograms – weight loss/thermogravimetry (TG) curves, rate of weight loss or differential thermogravimetry (DTG) curves, reaction progress/conversions α (evaluated as in Eq. (2)) and heat flow per unit weight – were plotted against reference temperature and time.

$$\mathrm{\alpha }=\frac{\mathrm{initial \ sample \ mass}-\mathrm{instantaneous \ sample \ mass}}{\mathrm{initial \ sample \ mass}-\mathrm{ final \ sample \ mass}}$$

(2)

During the thermal analysis, the condition of linearity of heating rate (minimal difference between sample and reference temperature) is satisfied. For the calculation of heat of pyrolysis (HoP) at 15 °C min⁻¹, DSC readings were blank corrected against an empty crucible and normalized by sample weight. They were then subtracted with a spline baseline since DSC readings of such heterogenous material would also be influenced by its specific heat capacity which tends to change with temperature. The HoP at any required region was calculated by the integration of this DSC curve (W g⁻¹) against time (s). A Nicolet iS50 FT-IR (Thermo Fisher Scientific) was used (resolution, optical velocity, aperture, and signal gain of 8 cm⁻¹, 0.4747 cm s⁻¹, 100, and 1.0, respectively) for tracing the evolved gas profiles of CO₂ (2400–2250 cm⁻¹), CO (2250–2000 cm⁻¹), CH₄ (3020–2800 cm⁻¹) and NH₃ (980–920 cm⁻¹) during the pyrolysis.

2.5 Scaling effect

A preliminary investigation of the scaling effect was carried out for the case of M010. After a 100 mL min⁻¹ flush with N₂ for 10 min, 12 sample sizes of M010 were linearly heated (20 °C min⁻¹ with 80 mL min⁻¹ N₂ flow) in 900 μL alumina crucibles from 25 to 650 °C. Twelve sample sizes from 5 to 500 mg (labelled as scl_5, scl_500 etc.) were used. Duplicate to triplicate trials was performed for each sample (more trials for sizes less than 25 mg) and averaged to plot the thermograms. Each crucible used in each trial was blank-corrected with the same heating program to minimize errors and DSC baseline drifts. The HoP was evaluated as mentioned earlier. The biochar produced at these different scales was then investigated for aromatic carbon content and total carbon (TC). TC was measured in the vario TOC cube (Elementar) through the combustion of 6 to 15 mg of biochar sample at 950 °C. Glutamine was used as the reference standard. Duplicate to triplicate measurements were made, depending on biochar availability. The 13C solid-state NMR (ssNMR) cross-polarization (CP) and direct polarization (DP) experiments are performed on a 600 MHz SB Bruker Avance III spectrometer (Bruker, Billerica, MA, USA) equipped with commercial 3.2 mm magic angle spinning (MAS) E-free ¹H,¹³C,¹⁵N probe head. Typically, around 30 mg of biochar was packed into a thin wall ssNMR rotor. All experiments are carried out at 16 kHz MAS rate and temperature of 275 K. Typical ¹H and ¹³C 90° pulse lengths were 3.5 μs and 5.0 μs, respectively; the length of CP contact time was 2 ms. A SPINAL64 decoupling [21] of 70 kHz field strength was applied during acquisition. In total, 32,768 scans with 2 s recycle delay were acquired for each CP spectrum. In total, 2048 scans with 15 s recycle delay was acquired for each DP spectrum. The CP and DP spectra were baseline corrected and normalized to the sample mass for further semi-quantitative analysis. The chemical shifts of aryl C, alkyl C and carbonyl C were assigned 110–165 ppm, 50–110 ppm, and 165–200 ppm, respectively. A relative spin counting for the relative NMR-observable carbon (RC) at the different mass scales was calculated according to Eq. (3), where I* is the normalized total integral between 50 and 200 ppm, and TC is the total carbon content.

$$RC=\frac{{I}^{*}}{sample \ mass \times TC}$$

(3)

Though total suppression of sidebands was not performed for CP MAS acquisition, the method of mathematical spinning sideband (SSB) subtraction was avoided because it could underestimate the intensities of SSBs arising from aromatic moieties and those from the hidden ones [22].

2.6 Statistics

The descriptive statistics of the analyses were calculated in OriginPro 2021, while statistical tests were performed using IBM SPSS v28.0. For elemental analysis, the measurement uncertainty was ± 0.01% for C, H, and N, and ± 0.015 for S. The mean std deviation of S_BET and ζ potential for the biochars was ± 0.419 m² g⁻¹ and ± 3.747 mV, respectively. The weight and heat flow data during the TGA/DSC trials were reproducible with the std deviations lying within an acceptable range. For instance, the mean and standard deviation of the thermograms of the mass-scaling experimental runs of M010 is shown in Fig. S2 of SM.

3 Results & discussion

3.1 Biochar yield and pyrolysis process

The biochar yield (wt%) and HPT (°C) of BM and MM during the lab-scale pyrolysis are shown in Fig. 1. Among the BM, the yield and the HPT during the process are proportional to their ash content – M100 > M001 > M010. Inorganics are a significant portion of MWB such as M100 and M001. For the MM, the yield of M505 is the highest while that of M523 is the lowest. The yield has no significant correlation (R² = 0.171) with the proportions of M010 and M001 blended with M100. This may be because the pyrolysis of such blended biomass is influenced more by synergetic effects amongst the different blend constituents as seen in refs [23, 24].

Fig. 1

The biochar yield and highest pyrolysis temperature (HPT) reached by the feedstock materials (MM and BM) during their lab-scale pyrolysis

3.2 Biochar properties

The elemental analyses – C, H, N, and O – of biochar are shown in Table 2. The molar ratios are calculated based on the total organic carbon (C_org) to avoid interference from inorganic carbonates inherited from M001 and M100. Biochars derived at temperatures above 400 °C can have aromatic clusters. The molar ratio of H/C_org is correlated to an n*n aromatic cluster size for H/C_org ⩽ 1 as in ref [25] (Eq. (4)). B514 and B550 can have large aromatic clusters up to 3*3 while others, except B001, may have up to 2*2. Thus, these biochars do not contain high molecular weight aromatic compounds. This may be because the HTT and residence time is insufficient to cause extensive polycondensation [26]. Hence, their porosity can also be low.

Table 2 Elemental analyses (in dry wt% of biochar prepared through lab-scale pyrolysis of BM and MM at 500 °C

$$\frac{H}{{C}_{org}}=\frac{2+4n}{2{n}^{2}+4n}$$

(4)

Table 3 shows the proximate analysis results (moisture, ash, VM and FC), ash-free biochar yield and TOR_i for B-BM and B-MM. Along with water, mass loss below 120 °C during the proximate analysis can also include adsorbed alcohols and low molecular weight tar released during pyrolysis. It can be seen that a) B550 and B010 have the highest (VM) due to the insufficient carbonization of holocellulose, b) B-MM have high ash content compared to biochar from lignin-rich feedstock like wood and walnut shells, and c) VM increased with the introduction of less-polymerized organic waste (M010) during co-pyrolysis. Ash in B-MM is mainly derived from M100 (Table S1 of SM) and co-pyrolysis helps in reducing this problem. The biochars B333, B505 and B514 exhibit a ~ 23% increase in (FC) relative to their respective parent biomasses, while others register a decrease. This is because FC in biochar with ash > 20% is poorly correlated with that of their parent biomass due to the influence of inorganics during pyrolysis, and interference of ash in the estimation of FC [27]. The ash-free yield of all B-MM is ~ 20 dry wt%. B100 is the most stable (TOR_i = 0.41) owing to its highest ash while B333 is least stable with TOR_i =—0.15 (despite having similar VM relative to other B-MM) which may be from the synergetic effects of its minerals and ash during the thermal oxidation process. The stability (mass loss during thermal oxidation) of the biochars is significantly (p < 0.05) correlated to total carbon (Pearson ρ = 0.847) and ash (ρ =—0.823).

Table 3 Proximate analyses results (water content, ash, VM, and FC), ash-free biochar yield, thermal oxidative recalcitrance (TOR_i), pore surface area (S_BET) and zeta (ζ) potential of B-BM and B-MM

For such mineral and ash-rich biochar (MAB), a conventional van Krevelen diagram is error-prone (since the determination of oxygen requires ash correction which is not reliable [28]) and a plot of H/C vs Ash/C is recommended [29]. Figure 2 shows the plot of H/C vs Ash/C and the van Krevelen diagram for B-MM and B-BM. Here, B333, B532, B550 and B514 fall under the definition of MAB. Literature suggests that pyrolysis temperatures between 350 to 500 °C are sufficient for the carbonization of biomass through slow pyrolysis [26]. Here, even the cellulose-rich M010 (H/C_org = 1.09) has not been sufficiently transformed, implying that the residence time is the limiting factor. Furthermore, except for B541 and B505, co-pyrolysis has generally increased the carbonization (H/C_org) relative to B100 as seen in the case of sludge and wetland biomass waste [30].

Fig. 2

A van Krevelen diagram (b) of the biochar types under study; and plot of H/C vs Ash/C (a) showing mineral and ash-rich biochar (MAB) defined by H/C ⩽ 0.7 and Ash/C ⩽ 4.1

Due to high ash, the O/C_org does not have a statistically significant correlation with VM (R² = 0.349 and p = 0.072). Similarly, there is no correlation between H/C_org and VM. The consensus that measured VM represents the labile C, –O and –N functional groups in biochar [31] is not always valid when the pyrolysis HTT is low compared to the temperature at which VM is measured because a) carbonates and some polyaromatics can also decompose [32] during the measurement of VM and b) variability in measurement techniques [33]. Thus, VM is not completely representative of the labile fraction in MAB. As per International Biochar Initiative (IBI) standards [34], based on H/C_org, the most carbonized biochar types here are B550 and B514, and they belong to class 2 (10% ⩽ C_org ⩽ 30%). At least 50% of their carbon will be stable in soil for 100 years (BC₊₁₀₀). The carbon sequestration potential of B550 and B514 (calculated as per ref [35]) are 0.21 and 0.22 kg_CO2 kg_biomass⁻¹ respectively. They can mitigate the N₂O emissions from the soil by ~ 40 ± 16% [36]. Despite having a similar VM, B550 is more carbonized compared to B010 at 500 °C due to the synergetic influences during co-pyrolysis. This further confirms that proximate analysis cannot be used to predict the stability of MAB.

The pores in these MAB are not well developed with S_BET < 10 m² g⁻¹ (Table 3) due to large ash concentrations, tar clogging the pores (also seen in the literature [37, 38]) and inadequate carbonization as seen above. B505 has the highest S_BET of only 8.1 m² g⁻¹, which tends to decrease with the ratio of M010 to M001 in the MM. This can be either from the coking of volatiles introduced by M010 that subsequently block pores [39] or from the volatiles that has not sufficiently carbonized and then occupies the porous network of minerals present in M100 and M001. Except for B505 and B100, others are colloidally stable (ζ < -30 mV). The negative surface charge of these biochar types stems from the acidic oxygenated functional groups from the uncarbonized organics that deprotonate in water. The molar ratios O/C and (N + O)/C should not be correlated with hydrophobicity and polarity index, respectively, as seen in ref [40] since most O is bonded in the form of silicates in ash and inorganic carbonates. Figure S3 (SM) shows the IR spectra of B-BM and B-MM and their calculated A_i. The peak regions 1470–1430 cm⁻¹, and 3600–3000 cm⁻¹ are from functional groups of O – CH₃, and –OH, respectively. The low A_i of B100 is due to its low C_org, while that of B010 is due to a lack of carbonization. The A_i does not follow the trend of H/C_org ratios of biochars. This is because the surface measurement nature of the ATR technique may not reveal the complete A_i in MWB biochars. However, it can be confirmed that, unlike in ref [41], co-pyrolysis of sludge has increased the aliphatic functional groups present in biochar. This is because, as in the aforementioned ref, up to 2 h of residence time may be required for sludge co-pyrolysis in this reactor.

Table S2 in SM shows the heavy metal, micronutrients (K, Mg, Na and Ca) and macronutrient (P, S) content (mg g⁻¹) in the biochar measured by ICP-OES. The synergy during co-pyrolysis influences the concentration of these elements (Fig. S4 of SM). M001 and M010 were able to reduce the Pb and Cu enrichment in M100, which is a major constraint in the application of sludge biochar [42]. This effect is usually attributed to the dilution effect that occurs during co-pyrolysis [39]. However, the Cu in B-MM is still higher than the land application limit of 0.10 mg g⁻¹. The enrichment of Cr, Ni, and Zn in B-MM are also above their limiting values of 0.09, 0.03, and 0.40 mg g⁻¹, respectively, as specified by the European Biochar Certificate [43]. Hence, B-MM is not suitable for direct land application, unless mixed with other soil ameliorators such as compost or organic fertilizers. However, despite the net HM enrichment, such biochar can have low exchangeable (acid soluble and carbonate bound) and a reducible fraction (bound to Fe and Mn oxides) – the most bioavailable forms of HMs. The exchangeable fraction of HMs can further transform into less bioavailable forms with increasing residence time. On the other hand, B-MM has higher concentrations of P and S which are essential plant nutrients. It is not recommended to increase the HTT to improve the carbonization of BM since it would concurrently raise the HM enrichment and pH. But, for biochars produced at HTTs around 500 °C, increasing the residence time may be another option to improve carbonization [44] and provide sufficient energy for the transformation of the acid exchangeable fraction of HMs to less bio-available forms.

During the analytical pyrolysis, the evolved products from the three B-BM are listed in Table 4. Since the biomass was pyrolyzed at 500 °C HTT, the evolved compounds during analytical pyrolysis at 800 °C also include the new pyrolysis products formed from any remaining uncarbonized biomass fraction. Peak tailing and double peaks were seen for some evolved compounds due to the large sample mass and the release of trapped volatiles from pores and the surface of biochar. The initial broad peaks are attributed to CO₂. Methoxylated mono-phenols and cresols are common products of lignin pyrolysis. Key markers of cellulose are levoglucosan (C₆H₁₀O₅) and 1,4:3,6-dianhydro-α-d-glucopyranose, while that of hemicellulose are ketones such as acetone and 1-hydroxy-2-butanone. The absence of these markers indicates negligible holocellulose remnants in B010 [45]. The absence of furans (mainly formed during holocellulose pyrolysis) also confirms this inference [46]. Long-chained aliphatic hydrocarbons (LAH) such as eicosane, dodecane, pentadecane and 1-heptadecene and recalcitrant organic aromatic compounds such as indole are formed from chain length shortening of aliphatic C in biomass. They represent an intermediate stage before aromatic condensation (through hydrogen abstraction and carbon addition) and are commonly present in biochar prepared at HTT below 500 °C [47]. BTEX compounds such as benzene, toluene, ethylbenzene, and xylene are markers of pyrogenic carbon from carbohydrates [48] and are also formed during analytical pyrolysis itself through the demethoxylation and/or dehydrogenation of methoxy or phenolic compounds. Chromatogram peaks of pyridine and benzonitrile that are more prominent in B100 and B001 signify that nitrogen functionalities have not been eliminated during the devolatilization of BM. The relative prominence of polyaromatic hydrocarbons (PAH, e.g., fluorene, anthracene, fluoranthene, benz[a]anthracene) and alkyl-PAH over LAH are seen in B100 and B001. More N-compounds (such as quinoline, isoquinoline and carbazole) are also seen in B100 and B001. PAH in B010 and B100 is limited to 4-rings, while B001 also shows 5-ring compounds – perylene and benzo[j]fluoranthene. For B-MM, the evolved compounds constitute more LAH with an increasing contribution of M010 in the blends (Table S3 in SM) and aromatics are limited to a maximum of 4-ring benz[a]anthracene. This further confirms that with more residence time these LAH can be transformed further, thereby increasing the carbonization of B-MM. Another discernible difference is that nitrous oxide seems to evolve from biochar derived from those blends where all BM were present – B514, B523, B532, B541 and B333.

Table 4 List of identified compounds, their relative %-area, and similarity index during the analytical pyrolysis (using Py-GC–MS) of B-BM

The properties of WEF from the biochar are shown in Table 5. The EC and pH of WEF have a strong negative correlation (ρ = -0.888, p < 0.01 and ρ = -0.791, p < 0.01, respectively) with the ash content, and a strong positive correlation (ρ = 0.933, p < 0.01 and ρ = 0.771, p < 0.01, respectively) with VM in biochar, respectively, similar to that seen in refs [49, 50]. Usually, when prepared at high HTTs, EC of biochar tends to increase with higher ash content due to the dissolved ions from the ash fraction. However, when biochars have large VM content, the dissolved ions from the VM dictate the EC. Here, for instance, EC amongst the B-BM is influenced mainly by K⁺ and Cl⁻ present in VM incorporated by M010. Except for B100, biochars are alkaline due to the carbonate and chloride forms of AAEM [50]. B100 is acidic due to the hydrolysis of Fe compounds present that releases H⁺ ions [49]. B100 is also saline (EC ⪆ 4 dS m⁻¹) due to the high concentrations of Cl⁻ and Na⁺ ions that may induce salinization problems, if applied to agriculture in arid regions. Water extractable N, P, K, Mg, Ca, and S are nutrients directly available to plants in the short term. B-MM are rich in SO₄²⁻ and K⁺, while water-leachable NO₃⁻–N is low and below the quantification limit of the instrument. After devolatilization, the remaining N must be in aromatic C structures that are non-water leachable as seen during the analytical pyrolysis. The low water-available form of P is common among biochars [31] and arises from the immobilization of P as minerals by multivalent metal cations of Al, Mg, Ca, and Fe under alkaline pH. Though this P is not readily mobilizable it can become a long-term source, which does not get lost through surface runoffs. Ions of Ca and Mg are much lower in the WEF due to the relatively low solubility of their carbonates and oxides in water. Hence, the B-MM and B-BM are not balanced commercial fertilizers [51]. The DOC from B010 is the highest since M010 has the most VM. However, SUVA₂₅₄ shows that this DOC mainly constitutes hydrophilic aliphatic compounds (since SUVA₂₅₄ < 3 L mg⁻¹ m⁻¹, except in B523, an outlier). N, P and S are important soil macronutrients and the maintenance of optimal carbon-nutrient ratios of C:N, C:P and C:S of biochar are essential for healthy soil microbial growth. These ratios of B-BM and B-MM compared with that of average soil [52] are shown in Fig. S5 (SM). Under the assumption of homeostasis and carbon-controlled microbial communities, B-MM tends to increase the carbon use efficiency of microbes as the nutrient ratios are at or below the soil limits. Thus, microbes can be energy limited and mobilize the excess N, P, and S, improving plant availability of these macronutrients in the short term due to low microbial demand. Hence, B-MM is suitable for addition in nutrient-deficient soil systems.

Table 5 Properties of water extractable fraction (WEF) from the B-BM and B-MM derived from lab-scale pyrolysis

3.3 Co-pyrolysis characteristics

As the heating rate increases from 5 to 20 °C min⁻¹, the rate of decomposition increases, and max peak decomposition shifts towards higher temperature in all reaction stages. This is because biomass is a poor conductor of heat and, even at small sample weights (1 to 10 mg), higher heating rates cause temperature gradients within the substrate owing to the ineffective heat transfer amongst the particles. The Figs. S6–S12 in SM show the thermograms – TG, DTG, conversions, and DSC curves – of all MM at six heating rates (5, 7, 10, 12, 15 and 20 °C min⁻¹) between 25 and 1000 °C. Furthermore, multiple reactions can concurrently initiate at higher heating rates (and cause higher degradation rates) thereby impeding the contact time of evolved volatiles with the substrate. This lowers the secondary char formation reactions and reduces the char yield as seen in the TG curves of all MM. Such changes in char yield are more pronounced during high-temperature gasification in presence of various AAEM and ash constituents of MM, which also renders this region relatively endothermic. The biochar yields are statistically independent (Sect. 4.2 in SM) of heating rates except for M532 and M541 due to the sample heterogeneity induced by M001 and M100.

The thermograms of all MM at the heating rate of 15 °C min⁻¹ are shown in Fig. 3. By 500 °C ~ 69–78% of the reaction is complete in all MM with a mean biochar yield of 52%. The identified pyrolysis reaction stages are (i) moisture removal ≤ 145 °C, (ii) organic matter decomposition – a combination of active pyrolysis and secondary heterogeneous reactions, (iii) decomposition of inorganic carbonates, and (iv) decomposition of other minerals and gasification above ~ 830 °C. In all but M532, stage 1 has an additional peak at ~ 130 °C possibly from the removal of chemically bound water. Detailed information on these stages is in Table S4 and S5 (SM). The increasing proportion of M010 in the co-pyrolysis blends lowers the peak onset temperature of stage (ii) from 207 to 187 °C due to the higher hemicellulose and catalytic effect of K as seen in another biomass thermolysis [53]. In this stage, the max peak decomposition temperature is delayed to 319 °C for M505 due to low synergetic influence from M010 and the max peak decomposition rate increases with the ratio of M010 to M001 due to the increase in organics with M532 as an outlier.

Fig. 3

The thermograms – TG, DTG, conversions, and DSC curves – during the pyrolysis of MM at 15 °C min⁻¹. In DSC, curves below 0 W g⁻¹ denote endothermic heat flow

From part-1 study [15] and Table S5, mixing M010 with M100 has shifted the active pyrolysis of sludge to lower temperatures, reduced the biochar yield and increased the exothermic heat release. An increase in temperature to 600 °C makes the process more energy intensive. This is in contrast to co-pyrolysis of SSW with wheat straw [23] and SSW with wood powder [54] where the exothermicity occurs mainly after 700 °C due to the late onset of secondary pyrolysis of char, which was also seen in the individual pyrolysis of M100 and M101 in the previous study. This difference depends on the temperature range at which the exothermic char condensation reactions takes place. When these reactions occur predominantly in the active pyrolysis regime, exothermic heat of pyrolysis is seen at lower temperatures as in co-pyrolysis of SSW with sugar cane bagasse [55]and with rice husk [56]. Here, this is another reason (in addition to HM enrichment) for limiting the pyrolysis of MM to 500 °C. Stages (iii) and (iv) are endothermic (DSC curve < 0) for all MM. However, between 700–900 °C, there is a small exothermic peak as in ref [23] due to interactions between the ash components and the biochar. This peak was absent in the pyrolysis of M100 [15]. The evaluation of deviations in weight loss (ΔW) between theoretical and experimental TG curves of the MM (Section S4.3 and Fig. S13 in SM) confirms the synergy between the components in the MM as seen in ref [13]. After stage 1, M505, M532, and M550 show a statistically significant (min std deviations) catalytic synergy during pyrolysis.

The relative emission intensities during active pyrolysis follow the order CO₂ > CH₄ > CO > NH₃ (Fig. S14 of SM). Here, CO₂ mainly evolves from the fracture of carboxyl and carbonyl groups of holocellulose. Therefore, M505 and M514, with their low holocellulose content, have the smallest CO₂ peaks. For M550, M541, and M532, the high K concentration has a catalytic effect on hemicellulose decomposition. Thus, the CO₂ shoulder peak from hemicellulose (at ~ 14 min) occurs at lower temperatures compared to other blends. The CH₄ release profile arises from the cracking of methoxyl groups in lignin and ceases after active pyrolysis with M532 as an outlier. The NH₃ (a NO_x precursor) profile diminishes with increasing amounts of M010 in the blends, while CO emitted during cellulose devolatilization is proportional to the M010 concentrations. The second peak of CO₂ and CO visible at ~ 44 min corresponds to the carbonate decomposition region. The evolution of CO beyond 700 °C is from the gasification of char in presence of CO₂.

3.4 Scaling effect during biomass pyrolysis

To the extent of the authors’ knowledge, the literature has rarely investigated whether the process characteristics and biochar properties usually studied via analytical pyrolysis (sample mass < 10 mg) at a particular HTT of lignocellulosic biomass are comparable to pyrolysis investigations using larger sample masses that induce heat and mass transfer effects. Among the BM, M010 is used for this preliminary scaling investigation since it has relatively lower ash and paramagnetic impurities that can influence the thermal gradient within the substrate during pyrolysis and the ¹³C ssNMR analysis, respectively at increasing mass scales. Figure 4 shows the rate of mass loss and heat flow during the pyrolysis of M010 at different mass scales between 5 and 500 mg in the TGA. In the normalized DTG curve, with increasing scales, the max peak devolatilization rate increases and the contribution of total mass loss shifts towards higher temperatures. This stems from the increased contact of evolved volatiles with the substrate leading to exothermic char-forming reactions. This is also seen in the total HoP (Fig. S15 in SM) which tends to shift from endothermic to exothermic after scl_75. There is also a statistically significant difference in heat flow between different scales (Sect. 4.4 in SM).

Fig. 4

Rate of mass loss (a) and heat flow curves (b) of M010 during pyrolysis at 20 °C min⁻¹ at different mass scales inside the TG reactor plotted against pyrolysis time. The legends scl-5, scl_10 etc., denote different mass scales inside TGA such as 5 mg, and 10 mg, respectively

The distribution of carbon moieties – aryl, alkyl, and carbonyl – based on CP and DP spectra, in the biochar produced in the TGA at different scales (at constant HTT 650 °C) are presented in Fig. 5 while Fig. 6 shows the relative observable carbon and total carbon in the biochar. The SSB occur in the CP and DP (Figs. S16 and S17 in SM) due to insufficient MAS compared to the anisotropy of the chemical shift signals in biochar. Furthermore, the non-protonated fraction of C has not been isolated from the total C signal and measured separately using the method of DP/MAS with dipolar dephasing of ¹H from carbon nuclei. Hence, the determination of C moieties based on CP and DP here is only semi-quantitative. The CP and DP spectra show aryl resonance at ~ 130 ppm. The assigned aryl C region can also include phenols and aromatic O–C such as those in recalcitrant lignin. The aryl C estimated by CP/MAS is higher than DP/MAS at all mass scales due to the stronger spillover of alkyl bands (sp³ O – C – O) to the aromatic region in biochar carbon matrix with higher protonation (H/C molar ratio). However, it can be observed that the aryl C content does vary with scales in the same pattern in both CP and DP, with the lowest values at scl_100 and scl_200. The relative observable carbon tends to decrease with scale in both CP and DP spectra despite DP directly measuring the nuclear spin polarization of the C nuclei. This could hint at the enrichment of paramagnetic impurities such as Fe, K and Mg that changes with scales. The total carbon seems to increase by more than ~ 16% from scl_400. It can be hypothesized that for a given pyrolysis reactor type (such as TGA) and HTT, there is a certain mass scale beyond which M010 show a predominant influence of char forming secondary reactions that maximize its exothermicity with carbon content and its distribution in the derived biochar also dependent on mass scales. The determination of the correlation between these secondary reactions and carbon content deems further structural analyses such as X-ray diffraction, N₂ physisorption, and Raman spectroscopy, which are beyond the scientific and financial scope of this study but explored in a follow-up investigation.

Fig. 5

The biochar produced from different mass scales of pyrolysis of M010 in TGA showing the distribution of carbon content (%) in biochar as a percentage of.¹³C NMR (a) cross polarization (CP) signal intensity and (b) direct polarization (DP)

Fig. 6

Biochar derived from M010 at different mass scales (mg) showing a) the relative observable carbon (RC) in arbitrary units (a.u) and b) total carbon content (%)

4 Conclusions

Mineral and/or ash-rich biomass such as sewage sludge from WWT constitutes a large proportion of the influx to solid waste treatment facilities in urban areas, especially in emerging economies. Their combined conversion into biochar for carbon sequestration can support local waste management and global decarbonization efforts. To this end, the biochar derived from the co-pyrolysis of digested sewage sludge with refuse lignocellulosic waste (banana peduncle) and anaerobic digestate from organic waste processing facilities are investigated in a lab-scale reactor to understand their capacities for carbon sequestration in soils. From the same materials and their blends, biochar was also produced and studied in a TGA/DSC-FTIR analyser to investigate the process characteristics and possible scaling effects. The major conclusions are summarized as follows:

• i. Yield, HPT and thermal oxidative recalcitrance of biochar during lab-scale co-pyrolysis are independent of blending ratios of M100, M010 and M001. Compared to parent feedstock, blended MWB tends to show their characteristics during pyrolysis due to synergy amongst the biomass components.

• ii. Co-pyrolysis of these MWB at 500 °C for 60 min was sufficient to depolymerize the lignocellulose. However, complete carbonization at this HTT is not achieved, which can be mitigated by raising the residence time in the lab-scale reactor. However, HTT above 500 °C can render the process energy intensive and increase HM enrichment in biochar.

• iii. Compared to SSW-derived biochar, better carbonization and reduced ashing were seen in the biochar when SS was co-pyrolyzed with other MWB. Amongst the biochars produced in lab-scale at 500 °C, B550 and B514 are the most carbonized (class 2 according to IBI). Except for B505 and B100, the biochar can be deposited in acidic soils due to the presence of NPK nutrients, and hydrophilic DOC.

• iv. The general pyrolysis emission profiles follow the order CO₂ > CH₄ > CO > NH₃. When the proportion of M010 in MM is increased, peak temperatures and NH₃ emissions during co-pyrolysis under 500 °C are lowered. The addition of M010 improved the exothermicity of active pyrolysis of M100.

• v. To minimize energy consumption, CO and NH₃ emissions, and HM aggregation the optimal temperature for co-pyrolysis is 500 °C. For the first time, it is experimentally proven that the scaling effect during analytical pyrolysis influences the char formation reactions, the carbon distribution (aryl, alkyl, and carbonyl C) and the total carbon content in biochar. For future modelling of pyrolysis systems, especially to forecast the heat of pyrolysis, this scaling effect must be considered.

The investigation confirmed a maximum carbon sequestration potential (up to 0.22 kg CO₂ kg⁻¹ biomass) of biochar derived from co-pyrolysis of SSW with other MWB, thereby aiding the solid waste management efforts in emerging urban areas such as Chennai. Further in-depth investigations into the plant-available heavy metals in this biochar (sequential speciation), long-term field tests in soil, technology costing, and the emission profile during pyrolysis scale-up are necessary. In the follow-up studies, authors explore the influence of scaling effects on analytical pyrolysis and biochar production from these MWB using an industrial-scale pyrolysis reactor in Chennai.