Supplementary Information
biodegradation Kinetics of POlycylic armotic hydrocarbons (PAHs) in soil-water systems: combined effects of DOM and biosurfactant enhanced biodegradation
Hui Yua, b, Guo-he Huanga, c. *, HuiningXiaob, Lei Wanga, Wei Chenc
aMOE Key Laboratory of Regional Energy Systems Optimization, S&C Academy of Energy and Environmental Research, North China Electric Power University, Beijing 102206, China
b Department of Chemical Engineering, University of New Brunswick, Fredericton, NB, Canada E3B 5A3
cEnvironmental Systems Engineering Program, Faculty of Engineering and Applied Science, University of Regina, Regina, SK, Canada, S4S 0A2
*Corresponding author Tel.: +1-306-585-4095. Fax: +1-306-585-4855
E-mail: (Guo-he Huang)
3tables
2 figures
DOM Characterization
Table S1shows chemical characterizations of two kinds of DOMs. The concentrations of total organic carbon in compost DOM (2128 mg C/L) was much higher than that in soil DOM (61 mg C/L); while the ash content of compost DOM (37.5%) was lower than that of soil DOM (50.1%). This indicated that organic matter in compost DOM was much higher than that in soil DOM. The N content in compost DOM (2.87%) was much higher than that in soil DOM (0.3%). As for the H/C ratio, the compost DOM exhibited a lower value than soil DOM, suggesting that compost DOM had relative higher aromatic nature. The E4/E6 ratios for soil and composting samples were 12.5 and 7.4, respectively, demonstrating that the compost DOM had a relatively higher molecular weight. When the DOM content in all samples were adjusted to 20 mg C/L, the specific UV absorptions at 254 nm were 1.050 and 0.761 Lmg−1 m−1 for DOMs derived from compost and soil, respectively. The aromatic carbon content in compost DOM was also higher than that in soil DOM.
Table S1 chemical characteristics of two DOM samples
DOM type / Total organic carbon (mgC/L) / Elemental content (%) / Atomic ratio / Ash content (%)C / H / N / H/C / C/N
Soil DOM / 61 / 35.89 / 5.678 / 0.294 / 0.16 / 112.1 / 50.1
Compost DOM / 2128 / 32.72 / 4.208 / 2.874 / 0.13 / 11.4 / 37.5
Fig. S1 shows FTIR spectra of the two kinds of DOMs. A broad band at 3200 to 3500 cm-1was assigned to H-bonds and OH groups; the peak of compost DOM showed a broader band than soil DOM. The band at 2870 to 2970 cm-1presented aliphatic carbons, of which two small peaks at 2964 and 2930 cm-1 were assigned to asymmetricalstretching of C-H in methyl and methylene groups, respectively. For soil DOM, a small peak at 2873 cm-1 was assigned to symmetrical stretching of C-H in methyl groups; a strong absorbance at 1420 cm-1 exhibited paraffinic characteristics. For the spectrum of compost DOM, there was a distinct peak at 1603 cm-1, implying polymer groups with C=C stretching vibrations in aromatic or vinyl groups and vibrations of COO- or H-bonded C=O. while there was no such kind of peak presence in the spectrum of soil DOM. This implied higher aromatic properties of compost DOM.
The 1HNMRspectraof twoDOMs are shown in figure S2. The peak shape of soil DOM showed simple and well-defined peaks, while the compost DOM presented complex and unresolved peaks. The identified peaks in the 1H NMR spectra were assigned to aliphatic H in methyl protons and main-chain methylene (0.8-1.5 ppm)
, carbonyl group in an acid or ester at β-C (1.8 ppm), ethers or hydroxyl group (3.6 ppm), esters (4.0 ppm), water (4.8 ppm) and carboxyl group (8.4 ppm). Both of these two samples showed these identified peaks with different responses. The soil DOM showed a strong peak in 3.6 ppm, which was attributed to protons to ethers or hydroxyl group; in comparison, the peak of compost DOM was unclear. Besides, for the spectrum of compost DOM, there were multiple peaks at 1.8 to 3.6 ppm. This indicated the component of high molecules of polymers; H may be bonded to aromatic C in methine group or bounded to O or N in aliphatic C. However, the detailed composition of the polymers needed to be further analyzed. In addition, there was a broad but weak peak at 6 to 8 ppm in the spectra of compost DOM, which was assigned to protons in aromatic and conjugated double-bound signals. In comparison, few signals in this band range were observed from the soil DOM. This also suggested relatively higher aromaticity of compost DOM.
(a)
(b)
Figure S1. FTIR spectra of DOM samples derived from (a)soil and (b)compost
NMR analysis
The difference of the DOMs was furtheranalyzedthrough1H NMR. The two extracted DOM samples also exhibited different peak shapes from the 1H NMR spectra (Figure S4). The peak shape of soil DOM showed simple and well-defined peaks, while the compost DOM presented complex and unresolved peaks, which indicated different compositions in the functional groups of DOM. The identified peaks in the 1H NMR spectra were assigned to aliphatic H in methyl protons and main-chain methylene (0.8-1.5 ppm), carbonyl group in an acid or ester at β-C (1.8 ppm), ethers or hydroxyl group (3.6 ppm), esters (4.0 ppm), water (4.8 ppm) and carboxyl group (8.4 ppm). Both of these two samples showed these identified peaks with different responses. The soil DOM showed a strong peak in 3.6 ppm, which was attributed to protons to ethers or hydroxyl group; in comparison, the peak of compost DOM was unclear. Besides, for the spectrum of compost DOM, there were multiple peaks at 1.8 to 3.6 ppm. This indicated the component of high molecules of polymers; H may be bonded to aromatic C in methine group or bounded to O or N in aliphatic C. However, the detailed composition of the polymers needed to be further analyzed. In addition, there was a broad but weak peak at 6 to 8 ppm in the spectra of compost DOM, which was assigned to protons in aromatic and conjugated double-bound signals. In comparison, few signals in this band range were observed from the soil DOM. This also suggested relatively higher aromaticity of compost DOM.
(a)
(b)
Figure S2. 1H NMR spectra of DOM samples derived from (a) soiland (b)compost
Table S2 Summary of biodegradation experiment design
Run # / Soil pretreatment / Solutions / Name in the paper1 / bulk / MSM( bulk contrast) / Bulk
2 / bulk / 0.1% NaN3 + MSM (sterile control) / Sterile control
3 / bulk / 100 ppm biosurfactant + MSM / 100 ppm biosurfactant
4 / bulk / 200 ppm biosurfactant + MSM / 200 ppm biosurfactant
5 / bulk / 500 ppm biosurfactant + MSM / 500 ppm biosurfactant
6 / bulk / 1000 ppm biosurfactant + MSM / 1000 ppm biosurfactant
7 / bulk / 50 ppm soil DOM + MSM / 50 ppm DOM added
8 / bulk / 50 ppm soil DOM + 100 ppm biosurfactant + MSM / 50ppm DOM+100 ppm biosurfactant
9 / bulk / 50 ppm soil DOM + 500 ppm biosurfactant + MSM / 50ppm DOM+500 ppm biosurfactant
(soil DOM+500 ppm biosurfactant)
10 / bulk / 50 ppm compost DOM + 500 ppm biosurfactant + MSM / compost DOM+100 ppm biosurfactant
11 / DOM removed / 500 ppm biosurfactant + MSM / DOM removed+500 ppm biosurfactant
12 / DOM removed / MSM
Bulk means original soil and without DOM removal
Kinetic model
For comparison of first order equation, a modified Monod equation or Michaelis-Menten kinetic model was applied to describe the PAH consumption rate [1-3].
(S1)
where C is the system concentration of PAH measured over time; t is the sampling time; Ks is the half-saturation constant for growth; KP is the equilibrium partition coefficient between soil and aqueous phase; Vmaxis the maximum degradation rate, which depends upon the total concentration of active enzyme in the bacteria, and Vmax can be expressed as, where µmax is the maximum specific growth rate and X0 is the initial cell density. The above model was the sole substrate Monod model which can be described as a no-interaction model. It represents a case where a compound in a mixture behaves as if it were the only compound without accounting for the effects from other existing substrates. When the mixture of substrates exhibit competition, the competitive inhibition kinetic model is applied[1,4].
, j ≠i(S2)
where Ci is the concentration of substrate i; Cj is the concentration of substrate j, , ,which can be defined as the inhibition constant, wherein KSi is the half-saturation constant for substrate i, and KSj is the half-saturation constant for substrate j; kj is the apparent first order rate constant; n is the number of competitive inhibitors.
EquationS2 can be extended for any number of components provided the compounds exhibit competitive inhibition kinetics. Thus, KSj represented a constant, while Ksi, qmax,i and Ci were fitting parameters. In this study, Equation S2 is specified as
(S3)
and
(S4)
Kinetic Equations S3 and S4 were integrated using the fourth-order Runge-Kutta numerical algorithm. Kinetic parameters could be derived through non-linear regression analysis using software Origin and 1stOpt, based on the Levenberg–Marquardt and Universal Global Optimization Algorithm. The biokinetic parameters of PHE and PYR fitting with two kinetic models are presented in Table S3.
Table S3 Comparison of kinetic parameters for biodegradation of PHE and PYR in soil-water systems
Biosurfactant concentration (ppm) / Soiltreatment / PAH / Modified Monod Equation / First Order Equation
Maximum degradation rate, Vmaxmg/(kg-soil·day) / Michaelis -Menten constant
Ki’, mg/kg-soil / Inhibition constant
кij / R2 / K
mg/(kg-soil·day) / R2
0 / DOM removed / PHE / 19.39 / 12.23 / 0 / 0.94 / 0.140 / 0.90
PYR / 19.39 / 150.97 / 1.37× 10-17 / 0.89 / 0.0642 / 0.96
No / PHE / 14.88 / 0 / 7.52 × 10-4 / 0.91 / 0.113 / 0.89
PYR / 8.64 / 112.3 / 1.68 × 10-15 / 0.92 / 0.0544 / 0.96
50 ppm soil DOM added / PHE / 13.84 / 0 / 8.06 × 10-17 / 0.92 / 0.100 / 0.87
PYR / 8.14 / 1 × 10-14 / 0.53 / 0.91 / 0.0397 / 0.92
100 / No / PHE / 19.34 / 11.22 / 0 / 0.95 / 0.144 / 0.89
PYR / 19.34 / 120.16 / 0.56 / 0.84 / 0.0773 / 0.96
50 ppm soil DOM added / PHE / 22.84 / 17.70 / 8.49 × 10-16 / 0.94 / 0.166 / 0.90
PYR / 22.84 / 84.43 / 0.38 / 0.84 / 0.0965 / 0.96
200 / No / PHE / 19.68 / 7.24 / 1.33× 10-15 / 0.96 / 0.161 / 0.91
PYR / 19.68 / 63.64 / 5.97× 10-16 / 0.76 / 0.0971 / 0.96
500 / DOM removed / PHE / 26.44 / 79.62 / 1.52 / 0.97 / 0.169 / 0.93
PYR / 26.44 / 478.24 / 0.76 / 0.85 / 0.105 / 0.97
No / PHE / 30.20 / 42.50 / 2.45× 10-15 / 0.90 / 0.186 / 0.94
PYR / 30.20 / 130.95 / 1.74 × 10-17 / 0.95 / 0.115 / 0.98
50 ppm soil DOM added / PHE / 30.69 / 30.03 / 2.87× 10-16 / 0.98 / 0.211 / 0.94
PYR / 30.68 / 124.88 / 1 × 10-14 / 0.89 / 0.139 / 0.98
50 ppm compost DOM added / PHE / 30.68 / 38.58 / 4.57 × 10-21 / 0.95 / 0.196 / 0.94
PYR / 30.68 / 172.39 / 5.17 × 10-16 / 0.80 / 0.121 / 0.99
1000 / No / PHE / 64.15 / 148.94 / 3.25× 10-15 / 0.98 / 0.250 / 0.99
PYR / 64.15 / 365.72 / 0.47 / 0.91 / 0.142 / 0.99
S1
References
[1]Desai AM, Autenrieth RL, Dimitriou-Christidis P, McDonald TJ. 2008. Biodegradation kinetics of select polycyclic aromatic hydrocarbon (PAH) mixtures by Sphingomonaspaucimobilis EPA505. Biodegradation 19:223-233.
[2]Lawrence AW, Mccarty PL. 1970. Unified Basis for Biological Treatment Design and Operation.J SanitEngDivAsce 96:757-778.
[3]Simkins S, Alexander M. 1984. Models for Mineralization Kinetics with the Variables of Substrate Concentration and Population-Density.Appl Environ Microbiol 47:1299-1306.
[4]Lotfabad SK, Gray MR. 2002. Kinetics of biodegradation of mixtures of polycyclic aromatic hydrocarbons.ApplMicrobiolBiot 60:361-365.
S1