3.1 Physicochemical characteristics of pretreated biochars
The main mechanism of organic pollutants removal for biochar is adsorption (Wang & Wang, 2019), which is highly influenced by its physiochemical properties such as surface area, functional groups, and pore size distribution. These physiochemical properties can be improved through various amendment strategies, such as chemical modifications (Wang & Wang, 2019). Chemically pre-treated biochar fibres have been found to possess increased surface areas, hydrophilicity, and oxygen-containing functional groups (Shen et al., 2008; Wang & Wang, 2019; Zhen et al., 2021). Results of the FTIR analysis shown in Fig. 1 (a) revealed that both pristine and modified fibres showed absorption in the region covering 400 to 4000 cm− 1 wavenumbers. Prominent peaks observed at 3421 cm− 1 can be assigned to –OH stretching of aromatic and aliphatic groups, while those at 2919 cm− 1 can be attributed to –CH asymmetric stretching vibration of methyl/methylene groups (Antonino et al., 2021). Shoulders related to C ≡ C and C-O-C stretching of aromatic groups can be seen at 1618 and 1085 cm− 1, respectively (Guo et al., 2021; Jiménez et al., 2020), whereas small broad peaks at 475 cm− 1 correspond to the stretching of C-H bonds of some aromatic and aliphatic functional groups (Guo et al., 2021). All fibres generally exhibited well-defined structures at similar positions. However, the acid-pretreated biochar exhibited slightly higher band intensities, indicating the relative abundance and concentration of existing functional groups.
BET results revealed a remarkable increase in surface area from 2.544 to 25.378 m2g− 1 for ARB (Table 1). In addition, the H/C ratio increased from 0.047 to 0.061. In contrast, the O/C ratio increased from 0.021 to 0.032 for the acid-rhamnolipid treated biochar (ARB), indicating that modifications with HCl improved the stability, polarity and the oxygen contents of the pretreated biochar. Finally, the SEM result of ARB at the end of the incubation period (Fig. 1b) indicated the presence of enlarged pore sizes of the biochar, likely caused by the acid modification.
Table 1
Physicochemical properties of pristine and modified biochar used in this study.
Biochar | Surface area (m2/g) | Pore size (nm) | Total pore volume (cm3/g) | Total elemental content (g/kg) | H/C ratio | O/C ratio |
C | H | N | S | O |
PB | 2.544 | 10.5 | 0.01 | 548.7 | 26.0 | 21.5 | 2.2 | 11.4 | 0.047 | 0.021 |
ARB | 25.378 | 17.8 | 0.24 | 452.1 | 27.7 | 21.3 | 1.9 | 14.6 | 0.061 | 0.032 |
3.2 TPH Degradation
Petroleum hydrocarbon removal in soils is generally characterized by rapid removal during the initial stage, followed by a slower and even plateau phase (Alexander, 2002; Qin et al., 2013). As shown in Fig. 2a, more than 51.8, 43.7, 39.3, and 25.8% of initial TPH was removed in ARB, PB + R, Rha, and PB treatment, respectively, within the first 21 days of incubation. This was mainly attributed to the metabolic consumption of the more readily available substrate by the bacterial community, facilitated by the presence of biochar or rhamnolipid in these treatments. A similar degradation pattern involving biochar and rhamnolipid has been reported by Wei et al. (2020). However, at week 4, the degradation rate slowed down (Fig. 2a) with a slight decline in removal amount from 123.3 to 96.7 mg TPH kg− 1 d− 1 in ARB, 104 to 85.3 mg TPH kg− 1 d− 1 in PB + R, 93.8 to 74.8 mg TPH kg− 1 d− 1 in Rha, and from 61.4 to 49.2 mg TPH kg− 1 d− 1 in PB treatments. After the introduction of the exogenous Scedosporium sp. ZYY at the beginning of the 5th week (day 29), degradation rate began to improve steadily in these treatments. It continued until the remaining days of the study, achieving a final removal rate of 82.4%, 73.6%, 65.7%, and 47.6% in ARB, PB + R, Rha, and PB, respectively, at 42 d. The physical dissipation of crude oil in the treatments from 1 d and at 42 d (week 6) of the experiment period is shown in Fig. S1.
Meanwhile, the control (natural attenuation) with a final TPH concentration of 4168 mg/kg did not achieve significant TPH removal due to the poor degradation capacity of the indigenous microorganisms in this treatment. This could be attributed to the slow acclimatization of the microorganisms to the presence of hydrocarbons and the overall substrate bioavailability limitation. The results suggest that the presence of biochar, especially in ARB, had a significant impact on TPH degradation within the first few weeks of the study. Biochars amended with acids result in increased surface area, improved hydrophilicity, and the introduction of oxygen-containing functional groups, which are crucial in activating the two main mechanisms (sorption and interaction between the functional groups) it employs for pollutant removal in contaminated soils. It is possible that the presence of ARB had an additional accelerative effect on hydrocarbon degradation, including providing a favorable surface for microbial colonization, which effectively shortened the contact distance between pollutants and microorganisms.
On the other hand, increased degradation rate and sediment enzyme activity (Fig. 2b) due to the addition of Scedosporium sp. ZYY highlights the effectiveness of bioaugmentation using specialized strains in bioremediation processes. Previous studies have shown fungi's advantages, such as the secretion of several low substrate specificity enzymes and the capacity of filamentous fungi to form mycelial networks to enhance access to hydrocarbon contaminants, among others, during hydrocarbons degradation in contaminated soils (Harms et al., 2011; Akbari et al., 2021). FDA activity assay conducted to monitor the degradation activities of microorganisms in different treatments revealed that microbial activities were higher in ARB, PB + R, and Rha treatments after fungi inoculation compared to the control treatment. This result corroborated with the TPH degradation results in these treatments.
Interesting findings were made in the present study, especially since no previous literature has reported on the sequential application of modified biochar and exogenous fungus in the slurry phase treatment of contaminated marine sediments. Notably, the degradation efficiencies using biochar, rhamnolipid, or bioaugmentation techniques reported in previous experiments on petroleum hydrocarbon-contaminated sediments in slurry systems are relatively lower than those obtained in the present work. For instance, Akbari et al. (2021) achieved the highest TPH removal rate of 21–40% at the end of 80 days from an initial TPH concentration of 2922 ppm, which was lower than the initial TPH amount used in this study. Few comparable studies where a higher TPH removal rate of 86% (Avona et al., 2022) and 100% (Smith et al., 2015) have been reported, involve the utilization of sophisticated laboratory bioreactor setups and longer experiment durations.
3.3 Changes in Bacterial Community Composition
High throughput analysis of the bacterial 16S rRNA gene was performed to describe the microbial dynamics in the OS and treated sediment at days 21 and 42 (Fig. 3). High-quality sequences detected from the OS and the five treatment groups revealed 95–99% bacterial composition, and 1% Archaea, while the remaining 4% were unassigned. Proteobacteria, Campilobacterota, Bacteroidota, Actinobacteriota, and Desulfobacterota represented the major phyla covering up to 96% of the total ASVs retrieved in the sample (Fig. 3a). Notably, ASVs affiliated with these phyla were more predominant in the treated samples than in the OS. Moreover, all treatments exhibited a slight increase in ASVs at 42 d for Bacteroidota, Actinobacteriota, and Desulfobacterota. However, an increased abundance of Proteobacteria at day 42 was only observed in the ARB treatment, while the abundance of Campilobacterota declined in ARB, Rha, and PB treatments. Proteobacteria are well-known efficient petroleum hydrocarbon degraders. Therefore, their high abundance in contaminated sites is an indication of their significant contribution to TPH removal (Zhen et al., 2021).
Analyzing the bacterial population at the genus level into Proteobacteria phylum also highlighted an increase in the total sequences, especially for the genera Erythrobacter sp., Sulfurovum sp., Alcanivorax sp., and Pseudomonas sp., accounting for 24.9%, 7.8%, 8.3% and 6.8% of the total ASVs, respectively (Fig. 3b). The genus Erythrobacter has been identified as a key degrader of petroleum hydrocarbons in marine sediments (Gao et al., 2015). Sulfurovum has been found to dominate microbial communities in the sediment–seawater interface, surface sediments, and deep-sea hydrothermal vents (Sun et al., 2020); while Alcanivorax and Pseudomonas genera represent two of the most studied groups of n-alkane degraders with high efficiency for crude oil bioremediation.
At the species level, a marked increase for sequences associated with Hyphomonas adhaerens (77%), Marinobacter mobilis (5%), and Erythrobacter seohaensis (4%), all belonging to Proteobacteria phylum was observed, particularly in ARB treatment at 42 d (Fig. 3c). Hyphomonas spp occupy various habitats in the marine environment and have been reported to belong to an opportunistic community of "fast" degraders in the surface environment (Tamburini et al., 2021). Members of Marinobacter spp. isolated from marine sediments produce highly stable surface-active agents for combatting marine oil spills (Raddadi et al., 2017) and have been regarded as specialists in degrading hydrocarbons, especially n-alkanes. Whereas, Erythrobacter spp. have been identified as potential degraders of PAHs (Kahla et al., 2021). Additionally, a significant (P < 0.05) increase in sequences affiliated with Methylophaga pinxianii was observed in the OS compared to the treated samples. Although previous studies have proven that Methylophaga spp. (Proteobacteria phylum) possess genes involved in the assimilation of aliphatic and aromatic compounds (Sauret et al., 2015; Song et al., 2022; Vila et al., 2010), their sparseness in the contaminated treatments is an indication that their proliferation is affected by higher degrees of oil pollution, such as the amount utilized in the present study. Overall, the abundance of hydrocarbon degraders and TPH removal in the ARB treatment suggested the presence of modified biochar and the exogenous fungus, Scedosporium sp. ZYY enhanced microbial-mediated degradation of TPH in the contaminated slurry sediment.
3.4 Diversities of Bacterial Communities during Degradation
The common and unique ASVs between the different sample groups are presented in Fig. 4. Each circle represents the various treatments at the specific sampling time. The number of circles and overlapping parts of the circle represents the number of common ASVs between the sample groups, while the numbers without overlapping parts represent the number of unique ASVs of the sample groups. Sample coverages were all > 99%, indicating the reliability of the data for diversity analysis (Fig. 5a). OS showed the highest α diversity with the largest Chao1, Shannon, and Pielou_e indexes (Fig. 5b-d). Previous studies in oil-contaminated samples have also observed decreased diversity due to crude oil (Looper et al., 2013; Wei et al., 2020). This can be attributed to the out-competition of other bacterial populations by specific oil degraders. In addition, an increase in the Chao1 index was only observed in ARB, PB, and Control treatments at 42 d compared to 21 d, while the Shannon index decreased at 42 d in all treatments except the Control and PB treatment. A high Chao1 index represents high species richness or diversity, while a high Shannon index represents high microbial community diversity (Liu et al., 2022). However, a slightly different trend in evenness was highlighted by the Pielou_e index (Fig. 5d), which increased in the Control, PB, and PB + R treatments but decreased in Rha and ARB. Overall, these results imply that Scedosporium sp. ZYY had little influence on the bacterial community evenness. Nevertheless, the presence of the fungus increased the species' diversity.
Principal coordinate analysis (PCoA) of bacterial community compositions based on the weighted UniFrac distance matrix was employed to study the β-diversity between samples. Principal coordinates accounted for 14·55% and 63·61% of the variance, respectively. As shown in Fig. 6a, treatments with high community structure similarity tend to cluster together, while treatments with significant community differences were farther apart. Samples from the OS were relatively concentrated in the area below the vertical axis' middle line, while samples from the treatment groups were concentrated in the area above, indicating the shifting bacterial community patterns between the OS and the treatment groups caused by the presence of oil. Similarly, Non-Metric Multi-Dimensional Scaling (NMDS) analysis was employed to highlight species abundance and distribution in the samples (Fig. 6b). According to NMDS analysis, the variations within the sample groups were small enough to aggregate into a cluster except in OS and some species in ARB.21 and PB.42. These findings indicated that the bacterial community composition in the original sediment was different from those in the treatments, and this corroborated with the result in Fig. 3c.
3.5 Implication of the Study
In addition to providing sufficient nutrients and conducive reproductive space for microorganisms (Zhang et al., 2021), previous studies have reported on the prospects of biochar in the remediation of soil organic pollutants (Kalsi et al., 2020). This study improved the typical adsorption properties of pristine biochar by introducing preparation conditions that influenced the adsorption and degradation of petroleum hydrocarbons. It has been reported that the diameter of most soil microorganisms is smaller than the average pore diameter of biochar (Zheng et al., 2022). Therefore, the increased rough surfaces and enlarged porous structures of biochar resulting from the modifications (Fig. 1b) enhanced soil aeration, nutrients, and water retention capacity. More so, it provided sufficient space for microbial colonization and a larger surface area for oil absorption, which eventually enhanced TPH removal (Fig. 2a).
Furthermore, previous reports have linked the significant increase in hydrocarbonoclastic bacteria and pollutant degradation to bioaugmentation with allochthonous fungi (Medaura et al., 2021). In this study, the subsequent addition of Scedosporium sp. ZYY into the treatments after 21 d of incubation improved the degradation rate by promoting the growth of TPH-degrading populations (Fig. 3c). These results suggest that biochar and exogenous fungus generally imposed synergistic functions in that biochar provided the microorganisms a conducive surface area for interfacial uptake of hydrocarbons. At the same time, the fungus promoted the proliferation of relevant species and increased their bioaccessibility to the target pollutants. This is the first study that reports the sequential application of biochar and exogenous fungus in the treatment of oil-contaminated clayey sediment. The findings of this study demonstrate the prospects of this application in the remediation of oil-contaminated clayey sediments. However, there remains a scarcity of research focused on the electron transfer mechanism of biochar in promoting microbial soil remediation and its application in removing various pollutants, which deserves further in-depth discussion.