- About this Journal
- Abstracting and Indexing
- Aims and Scope
- Annual Issues
- Article Processing Charges
- Articles in Press
- Author Guidelines
- Bibliographic Information
- Citations to this Journal
- Contact Information
- Editorial Board
- Editorial Workflow
- Free eTOC Alerts
- Publication Ethics
- Reviewers Acknowledgment
- Submit a Manuscript
- Subscription Information
- Table of Contents
International Journal of Chemical Engineering
Volume 2013 (2013), Article ID 450542, 6 pages
Evidence for PAH Removal Coupled to the First Steps of Anaerobic Digestion in Sewage Sludge
INRA, UR050, Laboratoire de Biotechnologie de l’Environnement, Avenue des Etangs, 11100 Narbonne, France
Received 7 September 2012; Revised 11 January 2013; Accepted 28 January 2013
Academic Editor: Claire Albasi
Copyright © 2013 Glenda Cea-Barcia et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
Anaerobic degradation of polycyclic aromatic hydrocarbons has been brought to the fore, but information on removal kinetics and anaerobic degrading bacteria is still lacking. In order to explore the organic micropollutants removal kinetics under anaerobic conditions in regard to the methane production kinetics, the removal rate of 12 polycyclic aromatic hydrocarbons was measured in two anaerobic batch reactors series fed with a highly loaded secondary sludge as growth substrate. The results underscore that organic micropollutants removal is coupled to the initial stages of anaerobic digestion (acidogenesis and acetogenesis). In addition, the organic micropollutants kinetics suggest that the main removal mechanisms of these hydrophobic compounds are biodegradation and/or sequestration depending on the compounds.
Polycyclic aromatic hydrocarbons (PAHs) are nowadays considered as environmental pollutants by environmental and health agencies because of their toxic, mutagenic, and carcinogenic effects on living organisms . PAHs are mainly formed during human activities such as crude oil spillage, fossil fuel combustion, and gasoline leakage. Through the air and runoff after rainy events, the PAHs can deposit to soil, water bodies, and sewage system [2, 3]. Due to their low water solubility and high octanol-water partition coefficients, these organic micropollutants (μP) are mainly associated with hydrophobic compartments such as the organic matter in sewage sludge and river sediments or the lipids in biota, with levels between 0.001 and 10 μg/gDM for PAHs [4–6]. μP can be also sorbed irreversibly in a short time scale to the organic matter. This phenomenon named sequestration has been reported by different authors as one of the abiotic mechanisms of μP removal in digested sewage sludge and soils [7–9].
PAHs are known to be biodegraded under aerobic conditions [10, 11]. However, most contaminated environments are anaerobic. In these environments, the anaerobic digestion can occur. Anaerobic digestion is a process whereby organic matter is broken down in the absence of oxygen into methane and carbon dioxide by naturally occurring microorganisms. The anaerobic digestion consists of four stages where different microbial populations participate: (1) hydrolysis, (2) acidogenesis, (3) acetogenesis, and (4) methanogenesis. The digestion process begins with the hydrolysis of insoluble organic polymers into monomers that are available for microorganisms. Acidogenesis is the step where these monomers are converted to carbon dioxide, hydrogen, and volatile fatty acids (acetic, propionic, butyric acids, etc.) by fermentative bacteria and then, in the acetogenesis stage, these metabolites are converted to acetic acid, hydrogen, and carbon dioxide by three different groups of bacteria (hydrogenic, homoacetogenic, and sulfatoreductive). Finally, the methanogenic microorganisms (strict anaerobic archaea) convert the acetic acid, hydrogen,d and carbon dioxide into methane [12, 13].
Recent results demonstrated that some PAHs [14–17] can be degraded under anaerobic conditions. It appears that this biodegradation is mainly coupled to methanogenic terminal oxidation processes [18, 19]. A recent thermodynamic landscape study showed that the biodegradation of PAHs under methanogenic conditions is an exergonic process; it indicated also that PAHs could be degraded to acetate and H2 coupled to the conversion of these substrates to methane as the optimal pathway . In addition, observations that bromoethanesulfonic acid, an inhibitor of methanogenesis, inhibited partially the degradation of naphthalene and phenanthrene with the total elimination of methanogenic archaea and some bacteria, suggesting that the PAHs degradation could be associated not only to the methanogenesis stage but also to the first stages of anaerobic digestion .
On the other hand, microbial population and cometabolism pathway were found to influence μP removal. Indeed, an adapted inoculum is favorable to quantify μP degradation , suggesting that the abundance and diversity of μP degrading microorganisms and their associated metabolic routes are essential. PAHs are recalcitrant compounds, but their biodegradation is possible if combined with other substrate degradation and if syntrophic interactions occur. These interactions can thus be stimulated by the addition of a readily biodegradable growth substrate [16, 21]. In fact, in the case where the sewage sludge is the growth substrate, it was demonstrated that the removal of PAHs by anaerobic digestion strongly varies as a function of sludge characteristics .
In order to further explore the μP biodegradation kinetics under anaerobic conditions in regard to the methane production kinetics, both total removal and removal rate of 12 PAHs were measured in different batch reactors series fed with secondary sludge as growth substrate.
2. Materials and Methods
All solvents were purchased from J.T.Baker. Fluorene (Fluo), phenanthrene (Phe), fluoranthene (Fluor), pyrene (Pyr), benzo(a)anthracene (B(a)A), Chrysene (Chrys), benzo(b)fluoranthene (B(b)F), benzo(k)fluoranthene (B(k)F), benzo(a)pyrene (B(a)P), dibenzo(a,h)anthracene (DBA), benzo(g,h,i)perylene (B(ghi)P), and indeno(1,2,3,c,d)pyrene (Ind) were obtained from Dr. Ehrenstorfer GmbH. The spiking mix at 100 mg/L was prepared in acetonitrile except for indeno(1,2,3,c,d)pyrene (20 mg/L). The 10 mg/L HPLC standard solution of PAHs in acetonitrile was provided by Dr. Ehrenstorfer GmbH. For quantification, the standard solutions were diluted to obtain 6 calibration levels from 10 to 1000 μg/L in acetonitrile. Standards were stored at −20°C.
2.2. Sludge Samples
The experiments were carried out using a highly loaded secondary sludge sample (HLSS) which was sampled in the biological aerobic unit of an urban plant treating 250000 person equivalent (PE) with a very low sludge retention time (0.4 day); it was, therefore, not stabilized. Prior to its direct use, HLSS was stored at −20°C.
2.3. μP-Adapted Anaerobic Inoculum
The inoculum (biomass adapted to μP) was produced in a semicontinuous stirred lab-scale reactor of 5 L. It was fed with HLSS diluted with deionized water to reach gCOD/L which was spiked with a mix of 12 PAHs (Section 2.1) at 5 μg/gDM for each PAH except for indeno(1,2,3,c,d)pyrene (1 μg/gDM). This reactor was inoculated with a digested sludge coming from a mesophilic anaerobic industrial reactor adapted to PAH-polluted sludge. Temperature was regulated at 35°C using hot water circulation in the external jacket. The feed was stored at 4°C. The reactor was operated at a retention time of 20 days and an organic load of gCOD/L/day. The pH and the biogas volumetric production were monitored on line. Once a week, composite samples were taken from the feeding and outlet tanks to measure the chemical parameters and μP concentrations. The adapted-to-μP digested sludge was kept at 35°C until its use for batch experiments as inoculum.
2.4. Batch Experiments
The anaerobic digestion of HLSS in batch conditions was performed using 300 mL plasma bottles containing 250 mL of a defined mixture of sewage sludge (substrate) and inoculum (biomass). The bottles were stirred (100 rpm) and maintained at 35°C. For these experiments, the sludge to inoculum ratio was equal to 0.5 gCOD/gOrganic Matter. The inoculum was concentrated from 8 to 18 g/L of dry matter (DM) to obtain a mixture with a total DM of 20 g/L. The mixture was spiked at 15 μg/gDM for each PAH and 3 μg/gDM for Ind. Two reactors series (A and B) were carried out (Table 1). The series A was performed in order to observe the overall activities with DM removal, biogas production, and μP removal. Based on the results of series A, indicating that μP removals occur before the production of methane, the series B was later performed. This series was carried out in duplicate to accurately analyze the concomitant evolution of μP and acetate during the first steps of the anaerobic digestion. At different reaction times, one reactor (series A) or two reactors (series B) were stopped to measure the total μP concentration over time (Table 1). Indeed, to quantify the μP, it is necessary to get a large amount of dry matter to extract them in duplicate. The two series were conducted under the same initial conditions, but for the series B, sampling frequency was higher in order to observe more precisely the μP kinetics and the anaerobic digestion kinetics. The initial μP concentration was determined also in duplicate for both series. In parallel, several control reactors were realized: nonspiked HLSS reactors (one per series) and HLSS abiotic reactors (Table 1). These abiotic control reactors were prepared by autoclaving the mixture-containing plasma bottles at 121°C for 20 min.
Biogas volume was measured by liquid displacement and its composition was measured by gas chromatography. For all reactors, volumes of biogas were reported to the COD introduced in the bottle. In addition, volatile fatty acids ((VFAs), g/L) production was followed in all reactors series.
2.5. Analytical Methods
Chemical oxygen demand ((COD), g/L), dry matter ((DM), gDM/L) and organic/volatile matter ((OM), gOM/L) were determined both at the beginning and end of the experiments. DM in total mixtures was measured by weighing the sample after heating at 105°C during 24 h and the OM after heating at 550°C during 2 h. The chemical oxygen demand COD was determined thanks to Merck Spectroquant kits, according to the ISO 15705. The samples were diluted to be in the range 0–1500 mgCOD/L. The VFAs acetate, propionate, isobutyrate, butyrate, isovalerate, and valerate concentrations were determined in dissolved/colloidal phase by gas chromatography (GC800, Fisons Instruments).
At different reaction time, the total volume of each reactor (series A and B) was freeze-dried in order to extract the PAHs according to the method developed by Trably et al. . The extractions were done in duplicate. The PAHs quantification in the extracts was previously described . These measurements gave the total μP concentration in μgμP/gDM.
2.6. Data Analysis
The μP removal rates and the accumulation rates of methane and acetate were used to correlate the μP removal kinetics to the anaerobic digestion steps. μP removal rates () were calculated using discrete points by difference of concentrations between two reactors (series A) or two pairs of reactors (series B) and divided by their respective difference of reaction times. These parameters were expressed in μg/L/day. In the case of series B, the were plotted in a boxplot graph which is a descriptive representation of independent data. Methane production rate () and acetate accumulation rate () were calculated with the production of methane or accumulation of acetate in a range of time expressed in mL/gCOD/day and g/L/day, respectively. The methane production first-order kinetic constant (, day−1) was determined as following: where ( is methane volume; is maximum methane production; is first-order constant; is time).
The total μP concentrations (μg/L) were calculated using the measured concentrations (μg/gDM) and the DM concentrations of each reactor. For the abiotic reactors, the μP removals were calculated on the basis of mass balance between the initial concentration of μP in reactors (day 0) and the concentration in reactors stopped at the end of incubation.
3. Results and Discussion
3.1. Production of μP-Adapted Anaerobic Digested Sludge
The semicontinuous anaerobic digestion of spiked HLSS was carried out to produce a homogeneous μP-adapted inoculum for the batch experiments, as previously described in the Section 2.3. All the results were previously published by Barret et al. . During the anaerobic digestion, total COD and DM were reduced by % and %, respectively. μP removals for all compounds ranged from 48% (Ind) to 85% (Fluo). Since μP volatilization was undetectable in the reactor, it was assumed that the calculated removals included biodegradation and sequestration [7, 9]; therefore, a μP-adapted inoculum was obtained.
3.2. Removal of μP in HLSS Anaerobic Batch Reactors
In order to describe the μP removal kinetics under anaerobic conditions, the total μP concentrations, DM and OM concentrations, biogas production and composition were measured over time in one of the HLSS batch reactors series (series A) and its corresponding nonspiked HLSS reactor.
In both reactors, the anaerobic digestion was stable with methane content in biogas at about 80% and pH 7. In the case of spiked reactors series, the lag phase lasted approximately twice as long as the lag phase of nonspiked reactors (Figure 1). This could indicate that there was an effect of the μP concentrations on the microbial population adaptation. However, the maximum volume of biogas was similar in both spiked reactors series and nonspiked reactor, which means that there was no inhibitory effect of the μP on the overall metabolism of anaerobic digestion. In contrast, a recent study showed that the anaerobic digestion in continuous reactors of wastewaters from industrial processes that contain very high concentrations of μP (included PAHs) is inhibited with concentrations around 16 gCOD, μP/L . This result suggests that a certain μP level could have an inhibitory effect on the anaerobic digestion. To compare, in our case, the initial COD concentration was 0.001 gCOD, PAH/L, which is much lower than the μP level of Puyol et al. .
Therefore, even if no inhibition occurred on the overall metabolism, the methane first-order constant () decreased by 44% in the case of spiked reactors (Figure 1). The values of the nonspiked reactors are 40% lower than those calculated using the methane produced from a secondary sludge in batch conditions published by Mottet et al. . These lower values of the nonspiked reactors may be explained by either our initial conditions with a lower sludge to inoculum ratio and a high DM initial concentration, or a lower methanogenic activity of our μP-adapted inoculum that was maintained a long time without feeding before use. However, these results show that the presence of μP slows down the overall anaerobic digestion kinetic but leads to the same methane potential. This may be due to either the expression of different metabolic routes to degrade the substrate and the μP, or to the deviation of electron fluxes towards μP, or to the μP toxicity/inhibition on some microbial groups.
are presented in Table 2 and compared with the values obtained in the semicontinuous reactor. The are similar between the two culture conditions except for B(a)A and B(ghi)P with higher values under a semicontinuous condition. This means that the microbial potential developed to remove the μP in the semicontinuous reactor was maintained and expressed itself during the batch experiment.
Fluor and DBA removals rates (, ) together with methane production rate () over time are presented in Figure 2. In all cases, the μP removal rates increase during the first 39 days of reaction, including the lag phase and the beginning of the exponential phase, and then decrease to zero. This behavior was also found for the other 10 PAHs, and most of them have the same profile as fluor, where the decreased to zero after 45 days. However, negative values of were obtained at 45 and 53 days. This could be due to concentration variability between each reactor.
These results suggest that the reduction of μP concentration is not only linked to methanogenesis, but it may be also due to abiotic processes at the beginning of the lag phase such as volatilization, adsorption to internal glass surface of the bottles, and/or sequestration [9, 25]. Considering these results, abiotic reactors and a second HLSS reactors series (series B) were carried out to precisely describe the μP degradation kinetics during the early stages of anaerobic digestion.
3.3. μP Degradation Kinetics during HLSS Anaerobic Batch Reactors
A second HLSS reactors series (series B) was carried out according to Section 2.4 to better describe the μP removal kinetics until the beginning of the exponential phase of methane production.
The VFAs concentration and methane production were followed in each reactor. Both VFAs and methane kinetics were similar in all reactors. For most of the μP, their concentrations were approximately constant until day 8 with no significant differences between the two values of concentrations of each pair of reactors, finding the highest differences for Chrys (4%) and B(ghi)P (5%). In the case of Fluo, Phe, and Fluor, important losses were found during the first 8 days mainly at day 1 of incubation. In this period, no important acetate and methane productions were observed, indicating that these compounds could have been removed by abiotic processes like volatilization, sequestration in the sludge matrix, and adsorption on glass bottles. This behavior was confirmed by the abiotic reactors results. All abiotic controls had a very low methane and VFAs production, indicating no microbial activity. The most important abiotic losses were found for Fluo, Phe, and Fluor (, , and ) after 54 days. These results are similar with those of Chang et al. , who found abiotic losses after 60 and 84 days of incubation with spiked batch reactors of 7% Fluo and 5% Phe.
Regarding the physicochemical characteristics of these compounds and according to the model proposed by Barret et al. , in the cases of Fluo, Phe, and Fluor, their highest volatility and lowest logKow (4.18, 4.57, and 5.1) suggest that the main removal mechanism at day 1 of incubation could be the volatilization. However, the methane production at day 1 is negligible, and thus the liquid-gas transfer is low. This suggests that the drop in the concentration is not only due to volatilization but also to sludge sequestration or glass adsorption.
Different , , and were calculated according to Section 2.6. ranged from μg/L/day (Fluo) to μg/L/day (Ind) for PAHs. Maximal were gCOD/L/day and are comparable to the values found in the literature . Figure 3 shows that the kinetic of is coupled to that of and . Indeed, increases when increases. When acetate is consumed to produce methane with a sharp decrease of , the decreases also. This behavior was also found for Fluo, B(b)F, B(k)F, and B(ghi)P and in the other six compounds, although with different kinetics. In the case of Pyr, started to increase after 18 days. Furthermore, Fluor, Phe, Chrys, B(a)A, B(a)P, and Ind maintained maximum values of even at the beginning of the exponential phase (18 and 34 days). These results reveal a relation between the first stages of anaerobic digestion (acidogenesis/acetogenesis) and the removal of μP.
The anaerobic degradation of PAHs in sewage sludge has already been studied but mainly under simplified batch experiments or using complex medium in continuous reactors. Unlike the continuous reactors, our methodology allows to observe the different steps of anaerobic digestion and to link these steps to the removal of μP. The results show a strong evidence of relation between the initial stages of anaerobic digestion (acidogenesis and acetogenesis) and μP removal in HLSS reactors. This relation was demonstrated by concomitant evolution of and . On the other hand, no complete inhibition of anaerobic digestion was observed in the presence of μP with a production of methane similar to the production without μP, but lag phases occurred implying microorganisms adaptation to the μP load and also the methane production kinetic was slowed down. The evidence of microorganisms adaptation to μP concentrations supports the results that linked the first stages of anaerobic digestion with the μP removal. In addition, the results show that the μP removals are also linked to two important abiotic processes depending on the compounds: volatilization (mainly the low-molecular-weight PAHs) and sequestration. This work provided important outcomes to better understand the PAHs anaerobic removal, to know its kinetics, and to identify the most important parameters, such as cometabolic activity, for optimizing their depletion.
The authors gratefully acknowledge CONICYT (Comisión Nacional de Investigación Científica y Tecnológica), Government of Chili and CROUS International, Government of France, and the Agence de l’Environnement et de la Maîtrise de l’Energie (ADEME-Convention no. 0906C0008 from the National Program of Research on Endocrine Disruptors (PNRPE)) for their financial support that made this study feasible.
- S. K. Samanta, O. V. Singh, and R. K. Jain, “Polycyclic aromatic hydrocarbons: environmental pollution and bioremediation,” Trends in Biotechnology, vol. 20, no. 6, pp. 243–248, 2002.
- Y. Wada, T. Ozaki, and M. Muraoka, “Characteristics of Benzo (a) Pyrene in highway runoff pollution,” in Proceedings of the 10th International Conference on Urban Drainage, Copenhagen, Denmark, August 2005.
- A. Welker, “Fate and behavior of selected organic substances in the urban drainage system,” in Proceedings of the 10th International Conference on Urban Drainage, Copenhagen, Denmark, August 2005.
- E. Abad, K. Martínez, C. Planas, O. Palacios, J. Caixach, and J. Rivera, “Priority organic pollutant assessment of sludges for agricultural purposes,” Chemosphere, vol. 61, no. 9, pp. 1358–1369, 2005.
- M. Blanchard, M. J. Teil, D. Ollivon, L. Legenti, and M. Chevreuil, “Polycyclic aromatic hydrocarbons and polychlorobiphenyls in wastewaters and sewage sludges from the Paris area (France),” Environmental Research, vol. 95, no. 2, pp. 184–197, 2004.
- D. Bodzek and B. Janoszka, “Comparison of polycyclic aromatic compounds and heavy metals contents in sewage sludges from industrialized and non-industrialized region,” Water, Air, and Soil Pollution, vol. 111, no. 1–4, pp. 359–369, 1999.
- M. Barret, G. C. Barcia, A. Guillon, H. Carrère, and D. Patureau, “Influence of feed characteristics on the removal of micropollutants during the anaerobic digestion of contaminated sludge,” Journal of Hazardous Materials, vol. 181, no. 1–3, pp. 241–247, 2010.
- G. L. Northcott and K. C. Jones, “Partitioning, extractability, and formation of nonextractable PAH residues in soil. 1. Compound differences in aging and sequestration,” Environmental Science and Technology, vol. 35, no. 6, pp. 1103–1110, 2001.
- K. T. Semple, A. W. J. Morriss, and G. I. Paton, “Bioavailability of hydrophobic organic contaminants in soils: fundamental concepts and techniques for analysis,” European Journal of Soil Science, vol. 54, no. 4, pp. 809–818, 2003.
- X. Y. Lu, T. Zhang, and H. H. P. Fang, “Bacteria-mediated PAH degradation in soil and sediment,” Applied Microbiology and Biotechnology, vol. 89, no. 5, pp. 1357–1371, 2011.
- S. Y. Yuan, S. H. Wei, and B. V. Chang, “Biodegradation of polycyclic aromatic hydrocarbons by a mixed culture,” Chemosphere, vol. 41, no. 9, pp. 1463–1468, 2000.
- B. K. Ahring, “Biomethanation I,” in Advances in Biochemical Engineering/Biotechnology, T. Scheper, Ed., p. 220, Springer, Berlin, Germany, 2003.
- D. J. Batstone, J. Keller, I. Angelidaki et al., “The IWA anaerobic digestion model no. 1 (ADM1),” Water Science and Technology, vol. 45, no. 10, pp. 65–73, 2002.
- A. Bernal-Martinez, D. Patureau, J. P. Delgenès, and H. Carrère, “Removal of polycyclic aromatic hydrocarbons (PAH) during anaerobic digestion with recirculation of ozonated digested sludge,” Journal of Hazardous Materials, vol. 162, no. 2-3, pp. 1145–1150, 2009.
- B. V. Chang, S. W. Chang, and S. Y. Yuan, “Anaerobic degradation of polycyclic aromatic hydrocarbons in sludge,” Advances in Environmental Research, vol. 7, no. 3, pp. 623–628, 2003.
- D. Dionisi, L. Bertin, L. Bornoroni, S. Capodicasa, M. P. Papini, and F. Fava, “Removal of organic xenobiotics in activated sludges under aerobic conditions and anaerobic digestion of the adsorbed species,” Journal of Chemical Technology and Biotechnology, vol. 81, no. 9, pp. 1496–1505, 2006.
- E. Trably, D. Patureau, and J. P. Delgenes, “Enhancement of polycyclic aromatic hydrocarbons removal during anaerobic treatment of urban sludge,” Water Science and Technology, vol. 48, no. 4, pp. 53–60, 2003.
- W. Chang, Y. Um, and T. R. P. Holoman, “Polycyclic aromatic hydrocarbon (PAH) degradation coupled to methanogenesis,” Biotechnology Letters, vol. 28, no. 6, pp. 425–430, 2006.
- D. M. Jones, I. M. Head, N. D. Gray et al., “Crude-oil biodegradation via methanogenesis in subsurface petroleum reservoirs,” Nature, vol. 451, no. 7175, pp. 176–180, 2008.
- J. Dolfing, A. Xu, N. D. Gray, S. R. Larter, and I. M. Head, “The thermodynamic landscape of methanogenic PAH degradation,” Microbial Biotechnology, vol. 2, no. 5, pp. 566–574, 2009.
- Y. Zhong, T. Luan, X. Wang, C. Lan, and N. F. Y. Tam, “Influence of growth medium on cometabolic degradation of polycyclic aromatic hydrocarbons by Sphingomonas sp. strain PheB4,” Applied Microbiology and Biotechnology, vol. 75, no. 1, pp. 175–186, 2007.
- E. Trably, N. Delgènes, D. Patureau, and J. P. Delgènes, “Statistical tools for the optimization of a highly reproducible method for the analysis of polycyclic aromatic hydrocarbons in sludge samples,” International Journal of Environmental Analytical Chemistry, vol. 84, no. 13, pp. 995–1008, 2004.
- D. Puyol, V. M. Monsalvo, N. Garcia-Mancha, A. F. Mohedano, and J. J. Rodriguez, “Anaerobic biological treatment of wastewaters from used oil recovery,” in Proceedings of the 5th European Bioremediation Conference, p. 189, Crete, Greece, 2011.
- A. Mottet, J. P. Steyer, S. Déléris, F. Vedrenne, J. Chauzy, and H. Carrère, “Kinetics of thermophilic batch anaerobic digestion of thermal hydrolysed waste activated sludge,” Biochemical Engineering Journal, vol. 46, no. 2, pp. 169–175, 2009.
- S. M. Bamforth and I. Singleton, “Bioremediation of polycyclic aromatic hydrocarbons: current knowledge and future directions,” Journal of Chemical Technology and Biotechnology, vol. 80, no. 7, pp. 723–736, 2005.
- M. Barret, H. Carrère, E. Latrille, C. Wisniewski, and D. Patureau, “Micropollutant and sludge characterization for modeling sorption equilibria,” Environmental Science and Technology, vol. 44, no. 3, pp. 1100–1106, 2010.