Scholarly article on topic 'A comparative study of anaerobic fixed film baffled reactor and up-flow anaerobic fixed film fixed bed reactor for biological removal of diethyl phthalate from wastewater: a performance, kinetic, biogas, and metabolic pathway study'

A comparative study of anaerobic fixed film baffled reactor and up-flow anaerobic fixed film fixed bed reactor for biological removal of diethyl phthalate from wastewater: a performance, kinetic, biogas, and metabolic pathway study Academic research paper on "Environmental engineering"

Share paper
Academic journal
Biotechnol Biofuels
OECD Field of science

Academic research paper on topic "A comparative study of anaerobic fixed film baffled reactor and up-flow anaerobic fixed film fixed bed reactor for biological removal of diethyl phthalate from wastewater: a performance, kinetic, biogas, and metabolic pathway study"

Yousefzadeh et al. Biotechnol Biofuels (2017) 10:139 DOI 10.1186/s13068-017-0826-9

Biotechnology for Biofuels

RESEARCH Open Access

A comparative study of anaerobic fixed film baffled reactor and up-flow anaerobic fixed film fixed bed reactor for biological removal of diethyl phthalate from wastewater: a performance, kinetic, biogas, and metabolic pathway study

Samira Yousefzadeh1'21, Ehsan Ahmadi3,4*1, Mitra Gholami5,6*, Hamid Reza Ghaffari7, Ali Azari4, Mohsen Ansari8, Mohammad Miri9, Kiomars Sharafi10 and Soheila Rezaei11


Background: Phthalic acid esters, including diethyl phthalate (DEP), which are considered as top-priority and hazardous pollutants, have received significant attention over the last decades. It is vital for industries to select the best treatment technology, especially when the DEP concentration in wastewater is high. Meanwhile, anaerobic biofilm-based reactors are considered as a promising option. Therefore, in the present study, for the biological removal of DEP from synthetic wastewater, two different anaerobic biofilm-based reactors, including anaerobic fixed film baffled reactor (AnFFBR) and up-flow anaerobic fixed film fixed bed reactor (UAnFFFBR), were compared from kinetic and performance standpoints. As in the previous studies, only the kinetic coefficients have been calculated and the relationship between kinetic coefficients and their interpretation has not been evaluated, the other aim of the present study was to fill this research gap.

Results: In optimum conditions, 90.31 and 86.91% of COD as well as 91.11 and 88.72% of DEP removal were achieved for the AnFFBR and UAnFFFBR, respectively. According to kinetic coefficients (except biomass yield), the AnFFBR had better performance as it provided a more favorable condition for microbial growth. The Grau model was selected as the best mathematical model for designing and predicting the bioreactors' performance due to its high coefficients of determination (0.97 < R2). With regard to the insignificant variations of the calculated Grau kinetic coefficients (KG) when the organic loading rate (with constant HRT) increased, it can be concluded that both of the bioreactors can tolerate high organic loading rate and their performance is not affected by the applied DEP concentrations. Conclusions: Both the bioreactors were capable of treating low-to-high strength DEP wastewater; however, according to the experimental results and obtained kinetic coefficients, the AnFFBR indicated higher performance. Although the AnFFBR can be considered as a safer treatment option than the UAnFFFBR due to its lower DEP concentrations


+Samira Yousefzadeh and Ehsan Ahmadi contributed equally to this work

3 Students' Scientific Research Center (SSRC), Tehran University of Medical

Sciences, Tehran, Iran

5 Occupational Health Research Center (OHRC), Iran University of Medical

Sciences, Tehran, Iran

Full list of author information is available at the end of the article

O© The Author(s) 2017. This article is distributed under the terms of the Creative Commons Attribution 4.0 International License

(, which permits unrestricted use, distribution, and reproduction in any medium,

provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license,

and indicate if changes were made. The Creative Commons Public Domain Dedication waiver ( publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.

in sludge, the UAnFFFBR had lower VSS/TSS ratio and sludge yield, which could make it more practical for digestion. Finally, both the bioreactors showed considerable methane yield; however, compared to the UAnFFFBR, the AnFFBR had more potential for bioenergy production. Although both the selected bioreactors achieved nearly 90% of DEP removal, they can only be considered as pre-treatment methods according to the standard regulations and should be coupled with further technology.

Keywords: Biofilm, Diethyl phthalate, Wastewater, Anaerobic treatment, Organic loading rate, Hydraulic loading rate, Biodegradation, Phthalic acid esters


Environmental and water pollution has become an issue of serious international concern in recent years [1-5]. In this regard, legislation requirements for discharging wastewater have recently become much stricter [6].

Emerging synthetic chemicals have entered into the environment through industrial activities [6, 7]. One of these emerging synthetic compounds is phthalic acid esters (PAEs), commonly named phthalates [8], which are widely used as plasticizers [9, 10] and globally applied in large quantities to make products such as plastics, pesticides, adhesives, paints, and cosmetics [11-13]. PAEs are ubiquitous pollutants as they have extensive applications in industrial processes, are not chemically bonded to products, and can migrate to environments [6, 14].

In recent years, PAEs have been considered as hazardous environmental pollutants [8]. They can have carcinogenic, mutagenic, and teratogenic potential impacts due to their ability to disrupt endocrine function [15-17]. Moreover, they can cause adverse health problems in humans' reproductive and developmental systems [1821] and have been associated with type 2 diabetes [22].

Therefore, some regulatory organizations including China National Environmental Monitoring Center (CNEMC), the United State Environment Protection Agency (USEPA), and the European Union have classified PAEs as top-priority pollutants [23].

Diethyl phthalate (DEP) is one of the important phtha-lates, the concentration of which was reported to be about 500 mg L-1 in discharged industrial effluent [8].

Although phthalates may not be significantly eliminated by physico-chemical processes like hydrolysis and photolysis, and other advanced treatment processes such as membrane and oxidation techniques are not advisable due to their high investment and operation costs, it has been shown that microbial biodegradation can play an important role in removing these pollutants from various environments [8, 24-26].

Meanwhile, a safe and cost-effective biological waste-water treatment is an important issue to consider [27, 28], especially for the industries which have to find financially feasible treatment methods to meet the permissible discharge levels [29]. Anaerobic biodegradation of

organic matters in wastewater is an economic way for wastewater remediation [30].

The logical reasons to use anaerobic treatment processes can be explained by considering the advantages of these processes which are related to their low energy demand, less sludge production, bioenergy (methane and hydrogen) production, and low nutrient requirement [31-36]. Furthermore, industrial wastewaters with high chemical oxygen demand (COD) containing high concentration of pollutants cannot be easily biotreated by aerobic bioreactors; this is while anaerobic bioreac-tors have been preferably and successfully applied to such wastewater [30, 31] and organic compounds of the mentioned wastewaters are ideal for green bioenergy production [37].

Attached growth (biofilm) anaerobic bioremediation has additional advantages over suspended growth bio-reactors, including more stability in operation, higher tolerance for pH, temperature, and toxic shocks, higher utilization rate of substrate, and ability to retain high biomass concentration under shorter hydraulic retention time (HRT) and overcome wash-out issues [30, 31, 37-39].

As highlighted above, microbial biodegradation can play an important role in removing PAEs and it is vital for industries to select the best treatment technology for all the wastewater with high concentration of phthalates [6]. Therefore, the main aim of the present study was to evaluate the performance of two different anaerobic biofilm bioreactors including anaerobic fixed film baffled reactor (AnFFBR) and up-flow anaerobic fixed film fixed bed reactor (UAnFFFBR), which have promising performance and have been recently used for industrial wastewater [40-42]. Moreover, as kinetic coefficients and mathematical modeling are two crucial tools that can be used to predict the performance of bioreactors and optimize their design [43, 44], they were used for comparing the selected bioreactors, since they can represent the microbial activity in different environmental conditions provided by the selected bioreactors. As kinetic and mathematical modeling directly relies on removal efficiency and microbial mass, they preferably, practically, and easily can be used by the designers and operators of

treatment plants, while the recognition of microbial type and their number cannot be used for predicting the performance of bioreactors [30, 45]. In addition, there are many analytical and conceptual challenges for profiling the diversity of microbial communities and diagnosis of their species in complex microbial communities such as wastewater [46, 47].

While many studies have been conducted to find a practicable way for evaluating the behavior of bioreac-tors applied to wastewater bioremediation, they could not completely evaluate the bioreactors' performance in terms of kinetic coefficients and mathematical modeling. Moreover, other studies have been only focused on calculating critical kinetic coefficients and computing removal efficiency, and have not investigated and evaluated their variations [48-50]. In addition, selecting the best biore-actors is traditionally done with their performance like COD removal, while their kinetic coefficients and bioen-ergy production are not usually considered very well. It should be noted that we previously applied mathematical modeling and kinetic evaluation for comparing the bio-degradability of different PAEs as substrates [51]; to our best knowledge, it is the first study for comparing the performance of the selected bioreactors by the mentioned methods. It should also be noted that there are other important parameters like the produced sludge quality and bioenergy production which need to be considered. Therefore, to fill the mentioned research gaps, this study was mainly aimed to compare two biofilm-based reactors including anaerobic fixed film baffled reactor (AnFFBR) and anaerobic fixed film fixed bed reactor (UAnFFFBR) for DEP removal from synthetic wastewater.


Implementation and operation of bioreactors

Two laboratory-scaled rectangular-shape bioreactors with the identical 6 L operating volume were used in this study. The reactors were made of 4 mm-thick Plexiglas. The AnFFBR was divided into three equal compartments by vertical baffles (each part contained 2 L operating volume and split by a further baffle) which were connected internally.

The dimensions of UAnFFFBR were 10 (length) x 10 (width) x 70 cm (height) (Additional file 1: Figure S1).

Both bioreactors were seeded with an anaerobic sludge obtained from a full-scale municipal wastewa-ter treatment plant (Ekbatan, Iran). Synthetic waste-water was continuously pumped from the feed tank into the AnFFBR using a dosing pump (Etatron, Italy). Another identical dosing pump was also used to feed the UAnFFFBR. For inoculating both the bioreactors (startup phase), further dosing pumps were used to recycle the washed-out sludge from the settling tank.

The bioreactors were operated at 25 ± 2 °C and fitted with the heaters to maintain the reactors' temperature stable. The bioreactors were filled by high-density polyethylene (HDPE) carriers, which acted as a fixed bed for biofilm (microbial) growth. The carriers had approximately 0.95-0.98 g cm-3 and 535 m2 m-3 of density and active surface area, respectively. The bioreactors were filled by 1.6 m2 of the carriers' available surface area (50% filling ratio).

In the acclimation step, glucose corresponding to 600 mg L-1 of chemical oxygen demand (COD) was used as the primary substrate and sole carbon source. Then, after reaching the steady-state condition, DEP was added stepwise and glucose concentration was reduced in parallel until DEP formed the entire carbon source. The experiments were conducted in two different study steps. The effects of HRT and organic loading rate (different concentrations of DEP) on the performance of the selected bioreactors were evaluated in study steps (A) and (B), respectively.

To have the chemical oxygen demand/nitrogen/phosphorous ratio of 350/5/1, ammonium chloride (NH4Cl) and ammonium bicarbonate (NH4HCO3) (for nitrogen source) along with mono-potassium phosphate (KH2PO4) (for phosphorous source) were used as nutrients in all the study steps. The composition of trace elements in synthetic wastewater was selected as follows: CaCl2-2H2O (14 mg), MgSO4-7H2O (90 mg), and 0.3 mL of trace solution per liter of synthetic wastewater. The following compounds were dissolved per liter to prepare the trace solution: KI (0.18 g), MnCl2-H2O (0.12 g), FeCl3-6H2O (1.5 g), CuSO4-5H2O (0.03 g), H3BO3 (0.15 g), CoCl2-6H2O (0.15 g), ZnSO4-7H2O (0.12 g), Na2MoO4-2H2O (0.06 g), and EDTA (10 g) [52]. Furthermore, NaHCO3 was applied to adjust pH at 7.5 ± 0.2.

Kinetics and mathematical modeling

Critical kinetic parameters and mathematical models have important impacts for predicting and designing biological wastewater treatment plants.

In the present study, three common substrate removal models, namely Stover-Kincannon, first order, and Grau (second order), were used to design the bioreactors and predict substrate removal rate. Under steady-state conditions (ddT = 0), if the first-order model prevails, the substrate consumption rate can be predicted by Eq. (1). Furthermore, Eq. (2) is the simplified form of Eq. (1):

So - S HRT

Q • So Q • S

= KiS,

- Kl S

where K1 is the kinetic constant for the first-order model (day-1), S0 and S are the influent and effluent substrate concentrations (mg L-1), respectively, Q is the inflow rate (L day-1), V is the reactor volume (liter), and HRT is the hydraulic retention time (day).

In the modified Stover-Kincannon model, the substrate utilization rate for the biofilm-based bioreactors can be determined by Eq. (3) which is based on organic loading rate. Equation (4) is also the linearized form of Eq. (3):

Kb + (q x s0)

Q x (S0 - S)

Q x So U

where is the maximum removal rate of the sub-

-1 day-1 day-1)

and KB is the saturation value

strate (mg COD L

constant (mg L-1 day-1). After the linear regression of Eq. (4), these coefficients can be calculated from the intercept and slope of the linear graph, respectively.

Equations (5) and (6) represent the second-order substrate removal model proposed by Grau. Equation (6) is also the integrated form of Eq. (5):

=ST = kg * x *( S0

= (n x HRT) + m

where n (dimensionless) and m (day-1) are the constants for a second-order model which can be determined from the intercept and slope of the plotted line of HRT versus Moreover, (X) is the suspended biomass con-(E) is the fractional substrate



centration (mg VSS L removal efficiency (dimensionless), and KG is the Grau second-order substrate removal rate constant (day-1). KG can be obtained from the following equation [53]:

Kg x X

Other kinetic coefficients including (KS) and (K) in the biofilm-based reactors can be obtained by combining mass balance equation [Eq. (8)] with Monod equation [Eq. (9)]. rs is the substrate utilization rate (g m-2 day-1).

V17 = QSo - QS - A(rs) at

r _ d5 _ K • S • Xatt

ks + s

In steady-state conditions, the substrate concentration changes in Eq. (8) can be ignored (-jf = 0) and Eqs. (8) and (9) given above can be combined as the following equation:

k f A ■ X(

Ks V Q(So - S)

Finally, KS as the half saturation constant (mg L-1) and (K) as the overall reaction rate (day-1) can be calculated from the linear regression of plotting | versus q^-S) line. Furthermore, Xatt which is the attached mass of biofilm (g VS) can be calculated by multiplying A (total available area, m-2) by XA:

Xatt = Xa x A.

Moreover, (Y) and (Kd) as the biomass yield coefficient (g VS produced/g substrate consumed) and biomass decay rate (day-1), respectively, can be calculated by mass balance equation (Eq. 12) and Monod growth kinetic (Eq. 13) [8]:

v dX = QXo - QXe + A(rg)

rg = Y (fsu) - Kd • A • X(a)

where r and Xe are the specific growth rate (g VSS m-2 day-1) and sloughed VS from the reactor (g VS day-1). If steady-state conditions are achieved (ddX = 0), Eq. (14) can be obtained by combining Eqs. (12) and (13) as follows:

(So - S) = Kd / Xatt

Xe " Y VQ • Xe

Finally, by the linear regression of Eq. (14), (Y) and (Kd) can be determined.

Subsequently, (^m) coefficient which is the maximum specific growth rate (maximum specific growth rate, day-1) can be calculated by multiplying (K) and (Y) coefficients expressed in Eq. (15). Eventually, influent substrate utilization rate can be calculated from Eq. (16), obtained by combining Eq. (9) with Eq. (14) [30, 49]:

^m = K • Y

Rsu —

[¿m • X • S0

Y (Ks + S0) '

Analytical methods

COD and biofilm characteristics including total solids (TS) and volatile solids (VS) were analyzed by the analytical techniques of standard methods [54].

Total organic carbon (TOC) samples were filtered by a 0.45 |jm filter and, subsequently, analyzed by TOC-Vcsh (Shimadzu, Japan) to measure the extend of mineralization.

Biogas measurements were performed as described in the study of Lay et al. [55].

Moreover, the analytical measurement of diethyl phthalate was performed using a gas chromatograph (GC) which was coupled with a flame ionization detector (FID) and equipped with a HP-5 capillary column.

To measure the DEP concentration, 10 mL of the wastewater sample was filtered through a glass fiber filter with 0.7 |im pore size and, subsequently, extracted by 2 mL of w-hexane solution. Finally, 2 |L of the extracted sample was injected into GC-FID. Concentration of the samples was measured by comparing them to the standard curve prepared at five points. Naphthalene was added as the internal standard.

GC temperature program was as follows: the oven's initial temperature was set at 70 °C for 1 min and, then, followed by a 10 °C min-1 ramp to 250 °C as its final temperature and was maintained at this temperature for 2 min. The injector and detector temperatures were set at 250 and 260 °C, respectively. To quantify the DEP concentration in effluent sludge (sloughed biofilm), 1 g of the freeze-dried sludge sample was ground and, subsequently, extracted by 10 mL of w-hexane.

To detect the metabolites of biodegradation, a gas chromatograph system equipped with a mass spectrometer (MS) detector and liquid chromatography tandemmass spectrometry (LC-MS/MS) were used.

Results and discussion

Evaluating performance of bioreactors

Experimental results for the AnFFBR and UAnFFFBR performance in terms of removing diethyl phthalate under different operation conditions including different HRTs and influent DEP concentrations are shown in Tables 1 and 2, respectively. It should be noted that data presented in Table 2 [study step (A)] which obtained under different HRTs have previously been used for comparing biodegradability of different PAEs [51]. The concentrations of COD, TOC, and DEP in the effluent of the bioreactors were determined as the bioreactor responses. The initial analysis for the COD and TOC experiments showed that each 1 g DEP L-1 produced 1.85 g COD L-1 and 0.62 g TOC L-1, respectively.

Analyzing impacts of hydraulic retention time

In study step (A), the influent concentration of DEP was maintained constant to 300 mg L-1 and HRTs were changed from 12 to 36 h. At this stage, the obtained results showed that the DEP removal rate for the AnFFBR can increase from 73.46 to 90.26% with the retention time increase from 12 to 36 h. Pirsaheb et al. [56] found similar results for the suspended growth anaerobic baffled reactor (ABR) and realized that higher HRTs resulted in more COD removal. The related COD and TOC removal efficiencies which showed similar behaviors and

Table 1 Performance of anaerobic fixed film baffled reactor (AnFFBR) in removing DEP

Study step Study phase A B

1 2 3 4 5 1 2 3 4

Influent DEP concentration (mg L-1) 300 300 300 300 300 400 500 600 700

HRT (h) 12 18 24 30 36 36 36 36 36

¿org (g COD m-2 day-1)a 4.162 2.775 2.081 1.665 1.387 1.85 2.312 2.775 3.237

Attached mass [TS (mg) of biofilm] 6180 4860 4180 3880 3710 4560 5400 6260 7170

Attached mass [VS (mg) of biofilm] 4500 3300 2820 2520 2400 3120 3900 4560 5340

VS/TS ratio 0.728 0.679 0.674 0.649 0.647 0.684 0.722 0.728 0.744

Effluent total suspended solids 334.3 205 137.7 111.4 90 141.6 168.8 239.5 276.2

(mg day-1)

Effluent volatile suspended solids 224.4 130.2 84.5 64.2 48.4 81 98.6 144.3 168.9

(mg day-1)

DEP concentration in TSSe (mg g-1)b 8.5 8.2 6.6 6.2 4.1 7.3 7.5 7.5 7.9

SRT (day) 18.48 23.7 30.35 34.83 41.22 32.2 31.99 26.13 25.96

Methane production (L/g CODrem) 0.21 (0.19)c 0.28 (0.25) 0.32 (0.29) 0.37 (0.34) 0.44 (0.4) 0.38 (0.35) 0.35 (0.32) 0.34 (0.31) 0.31 (0.28)

Methane percentage (%) 42.2 44.2 52.5 56.4 64.8 63.3 62.1 61.7 61.1

DEP removal (%) 73.46 78.43 81.4 86.33 90.26 90.62 90.64 90.83 91.11

COD removal (%) 66.66 70.1 77.33 82.29 87.01 88.04 89.48 89.6 90.31

TOC removal (%) 51.82 57.95 66.82 75.7 83.17 83.06 81.90 82.5 82.57

a Organic loading rate b Effluent total suspended solids

c Numbers in brackets are methane yield at STP (1 bar and 273.15 °K)

Table 2 Performance of up-flow anaerobic fixed film fixed bed reactor (UAnFFFBR) in removing DEP

Study step Study phase A B

1 2 3 4 5 1 2 3 4

Influent DEP concentration (mg L-1) 300 300 300 300 300 400 500 600 700

HRT (h) 12 18 24 30 36 36 36 36 36

¿0rg (g COD m-2 day-1) 4.162 2.775 2.081 1.665 1.387 1.85 2.312 2.775 3.237

Attached mass [TS (mg) of biofilm] 5840 4370 4110 3880 3720 4340 5440 6080 6840

Attached mass [VS (mg) of biofilm] 4080 2940 2700 2460 2280 2880 3660 4320 4980

VS/TS ratio 0.698 0.672 0.657 0.634 0.613 0.663 0.672 0.71 0.728

Effluent total suspended solids (mg day-1) 326.8 227.3 169.3 129.6 112.2 144.5 181.4 220 270

Effluent volatile suspended solids 164.4 102.4 74.4 55.7 48 76.8 96 117.6 151.8

(mg day-1)

DEP concentration in TSSe (mg g-1) 9.7 9.3 8.6 8.1 5.2 7.4 7.5 8.4 9.1

SRT (day) 17.87 19.22 24.27 29.94 33.15 30.03 29.98 27.63 25.33

Methane production (L/g CODrem) 0.22 (0.2) 0.24 (0.22) 0.28 (0.26) 0.34 (0.31) 0.42 (0.38) 0.34 (0.31) 0.32 (0.29) 0.30 (0.27) 0.27 (0.24)

Methane percentage (%) 41.1 42.3 52.8 56.1 61.5 61.0 59.5 59.6 57.3

DEP removal (%) 67.9 72.5 76.46 80.96 87.86 88.15 88.3 88.36 88.72

COD removal (%) 57.91 60.86 71.85 78.79 83.55 84.25 86.02 86.27 86.91

TOC removal (%) 45.32 53.97 62.58 70.05 78.33 78.34 78.25 78.79 79.65

Lorg organic loading rate, TSSe total solids, US volatile solids, TSSe effluent concentration of total suspended solids

maximum removal of COD and TOC was measured as 87.01 and 83.17%, respectively, in 36 h of HRT. The DEP removal in the AnFFBR was considerably higher than its mineralization (COD and TOC removal), which can be related to the presence of benzene ring (aromatic structure) known to be more refractory to biodegradation [8, 57]. The experimental results of study phase (A-1) which compared DEP and TOC removal confirmed this theory. Nearly 21.64% higher removal rate of DEP than TOC showed that the side chains of DEP were biodegraded faster than its benzene ring. This difference reduced to 7.1% by increasing HRT to 36 h.

Similar observations were recorded for the UAnFFFBR. The improved effluent quality of both bioreactors at higher HRTs can be related to the higher contact time of microbial mass to DEP. Another critical factor is solid retention time (SRT), which shows the presence time of microbial mass to fulfill the organic substrate biodegradation. In study step (A), the SRTs of both reactors notably increased by reducing the organic loading rate and increasing HRT. The highest SRTs for study step (A) were observed in 36 h HRT, which were equal to 41.22 and 33.15 (days) for the AnFFBR and UAnFFFBR, respectively. These obtained SRTs were considerably higher than other suspended growth processes (e.g., conventional aerobic activated sludge) ranging from 3 to 30 days. It was also known that SRTs between 5 and 50 days are ideal for xenobiotic compounds to be biodegraded which can be considered as an advantage for both of the studied reactors [30].

The effluent quality of the UAnFFFBR, especially in study phases (A-1) and (A-2), was significantly less than that in the AnFFBR. This result revealed the AnFFBR could tolerate higher organic and hydraulic loading rates than the UAnFFFBR. Greater performance and tolerance of the AnFFBR can be the result of phase separation known as the significant advantage of anaerobic baffled reactors (ABR) [58]. This phenomenon may cause acido-genesis and methanogenesis bacteria to be divided into two different phases in the AnFFBR, which subsequently make it possible for these bacteria to grow and synthesize under their favorable conditions. In this regard, the growth rates of the attached biofilm mass (as mg VS and TS) were higher for the AnFFBR than UAnFFFBR. Furthermore, as the volatile solids (VS)-to-total solids (TS) ratio of the biofilm were higher for the AnFFBR than UAnFFFBR, it can be concluded that the microorganisms had greater activity in the AnFFBR.

Moreover, higher SRTs in the AnFFBR are due to greater persistence of microbial mass to the attained conditions, which consequently leads to more effective bio-reaction in the AnFFBR. However, both the bioreactors performed well at 300 mg influent DEP concentration per L and 36 h HRT, resulting in 90.26 and 87.86% of DEP removal to be achieved for the AnFFBR and UAnFFFBR, respectively.

Analyzing effect of organic loading rate

The bioreactors' performance under different diethyl phthalate concentrations and organic loading rates

(stable hydraulic loading rate of 4 L day-1) was evaluated in study step (B). According to the bioreactors' best performance observed in study phase (B-4), 90.31 and 86.91% of COD and 91.11 and 88.72% of DEP were removed by the AnFFBR and UAnFFFBR, respectively. In this step, the variation of removal rate slightly increased as the substrate and organic loading rate increased from 400 to 700 mg DEP L-1 and 1.85 to 3.237 (g COD/m2arrier day), respectively. It can be concluded that HRT had a greater impact than organic loading rate on these biofilm reactors. Although the bioreactors' removal efficiencies increased in study step (B), the COD and DEP concentrations slightly increased in the effluent. Moreover, VS/TS ratio and active biofilm mass are two parameters which can affect the bioreactors' performance in study step (B). It can be concluded that, by increasing organic loading rate, both the attached microbial masses and their volatile portion increased and, consequently, led to higher biodegradation capacity of the bioreactors (Tables 1, 2). It should be noted that better effluent concentration and performance of the AnFFBR than those of the UAnFFFBR in study step (B) can be related to its higher biofilm mass and VS/TS ratio. Another factor is phase separation which was previously mentioned. Moreover, it has been known that anaerobic reactors' performance in higher organic loading rates and for strong wastewa-ters (more than 1000 mg L-1 of influent COD) is better compared to when they are used for low strength waste-water treatment [31]. The experimental results of this study (which indicated stable performance of the anaerobic bioreactors) along with Farzadkia et al.'s investigation can demonstrate the advantage of these anaerobic biofilm reactors over aerobic bioreactors [59]. In the study by Farzadkia et al., high organic loading rate of the substrate had an adverse effect on the biodegradation and metabolic activity of aerobic fixed bed activated sludge hybrid reactor, and the removal efficiency of the reactor dropped from about 96-79% when organic loading rate was increased from about 1 to 4.5 kg COD m- ~3 day-1 [ ].

It should be also noted that the statistical t test analyses for DEP, COD, and TOC removal efficiencies showed the mentioned parameters to be statistically different in most of the study phases when compared between the bioreactors (P value <0.05).

A summary of some studies conducted on the evaluation of biological removal of PAEs is presented in Table 3. It should be noted that, although all these biological methods indicate that bioremediation can play an important role in removing PAEs, the type and condition of bioreactors including (anaerobic, aerobic, and anoxic conditions), their operational conditions (particularly HRT and SRT), and finally type of the selected phthalate can impact the performance of bioreactors [8, 60-65].

It can be observed that these bioreactors achieve more removal rates as they increase the HRT and SRT, which are in agreement with the observed results of the present study [8, 60, 64]. In addition, there is a diverse relation between the length of alkyl-side chains of phthalates and their biodegradability [8, 61].

Biogas production

Methane, carbon dioxide, hydrogen, and hydrogen sulfide were observed as the main gases in biogas. Both methane production rate and its percentage for bioreactors were affected by HRT and organic loading rate. According to the results, the methane production increased by increasing HRT (study step A) and, then, decreased by increasing organic loading rates (study step B). It is commonly known that operational parameters can critically affect the biogas production rates [66]. The maximum methane production rates were observed in study phase (A-5) with 36 h of HRT as 0.44 and 0.42 L CH4 g COD,-^ at 1 bar and 298.15 °K (or 0.4 and 0.36 L CH4 g COD-! at 1 bar and 273.15 °K) for the AnFFBR and UAnFFFBR, respectively.

For most of the study phases, the methane produced by the AnFFBR was considerably higher than that produced by the UAnFFFBR, which may be the result of phase separation provided by the baffles of AnFFBR that consequently help the acetogenic bacteria to utilize the volatile fatty acids (VFAs) produced by acidogenic bacteria before these metabolites reach methanogenic bacteria and inhibit their activities. This issue can be confirmed by the results obtained in the study by Wang et al., indicating that some VFAs like propionic acid can significantly inhibit the methanogenesis process [67]. This hypothesis can be supported when the methane yield of the UAnFFFBR in study step B more significantly decreased than that of the AnFFBR by increasing organic loading rate. This phenomenon can be attributed to greater accumulation of some inhibitory compounds (e.g., propionic acid) in the UAnFFBR. Moreover, the study results of step (A) demonstrated that methanogenic bacteria needed more time to convert the VFAs and other metabolites into methane. It can be also observed that increasing HRTs from 12 to 36 h for each of the biore-actors can lead to the increase of the methane production for about twice. These results could be attributed to slow growth and synthesis rate of methanogen bacteria, which results in more SRTs requirement for optimum growth and the substrate concentration increase needed by these bacteria (e.g., acetic acid) in higher HRTs. However, some of the obtained methane yields in the present study, especially in study phase (A-5), were higher for both of the bioreactors compared to some other studies which have commonly reported the methane yields of

Table 3 Summary of some studies conducted for biological removal of PAEs

Treatment method

Substrate type (PAEs)

Experimental conditions and removal efficiency

Important observation

Anaerobic/anoxic/oxic (AAO) treatment system

Trickling filter

Cyclic activated sludge technology (CAST), anoxic/oxic (AO), and anaerobic/anoxic/oxic (AAO) processes

Moving bed biofilm reactor (MBBR)




Up-flow anaerobic sludge blanket reactor (IJASB) DMP

Anaerobic/anoxic/oxic treatment system (AAO) DnBP

Submerged membrane bioreactor (aerobic condition)

The optimal HRT and SRT for DMP and nutrients removal were 18 h and 15 days, respectively, and the degradation rates of anaerobic, anoxic and aerobic zones for DMP were 13.4,13.0, and 67.7%, respectively

Trickling filter achieved 94-99% of DEP and 44% of DEHP removal

The overall removal efficiency of all the selected PAEs was more than 72% in CAST while AO and AAO only achieved about 30% of PAEs removal

In optimum conditions with HRT of 9 h, about 95 and 94% of DEP and DAP were removed, respectively. In addition, more than 92% of COD removal was achieved for both phthalates

More than 99% of DMP and 93% of COD were removed from the wastewater containing 600 mg L-1 of DMP (corresponding to 3 g COD L-1 day-1) of organic loading rate of) at 8 h of HRT

The optimal HRT and SRT for DnBP removal were 18 h and 15 days, respectively. In the mentioned conditions, about 72.66% of DnBP was degraded by the selected process, 24.44% was accumulated in the system, 2.44% was released in the effluent, and 0.5% remained in the waste sludge

The removal efficiency of DEHP under HRTs of 4 and 6 h, SRT of 140 days, and sludge concentration of 11.5 and 15.8 g VS L-1 ranged between 91 and 98%

The biodégradation process of DMP by the [60]

selected method was in accordance with the first-order kinetics model. Under the optimal conditions, about 73.8,5.8,19.3, and 1.1% of DMP was biodegraded, released in the effluent, accumulated in the system, and remained in the waste sludge, respectively

DEHP was the most recalcitrant among the [61]

selected phthalates and DEP with less molecular weight was biodegraded with higher rate

The better performance of the CAST process was [62] attributed to its better indoor-conditions for bacterial community

MBBR tolerated both of the selected phthalates [8] with the influent concentrations of 100-300 mg L-1. Sludge yield constants (V) varied from 0.3 to 0.54.The maximum specific growth rates were 0.37 and 0.34 day-1 for DEP and DAP, respectively

The sludge yield was estimated as 0.08 g [63]

VSS g COD-1. DMP was first de-esterified to mono-methyl phthalate (MMP) and, then, to phthalate (or phthalic acid) and, subseguently de-aromatized and converted into methane and C02.The maximum specific degradation rates of DMP, MMP, and phthalate were 415,88, and 36 mg (g VSS day)-1, respectively

Increasing SRT from lOto 15-25 days resulted in [64] increasing removal efficiency from 90 to 95%. Higher SRT improved DnBP biodégradation efficiency. The removal efficiencies of anaerobic, anoxic and oxic reactors were 17.14,15.02, and 63.46% of the total DnBP removal, respectively

The removal of DEHP was closely dependent on [65] the membrane pore size and about 74% of inlet DEHP was biodegraded

DMP dimethyl phthalate, DEHP bis(2-ethylhexyl) phthalate, DnBP di-n-butyl phthalate, BBP butyl benzyl phthalate, OOP di-n-octyl phthalate, D/\Pdiallyl phthalate

about 0.27 to 0.36 L CH4/g COD [68, 69]. Moreover, the methane production rate is usually expected to be placed around 0.35 L CH4/g COD in STP condition (or 0.382 L CH4/g COD at 1 bar and 298.15 °K), which is the theoretical methane yield [34, 69, 70]. The obtained higher methane yields in this study might be the result of further conversion of biosolids produced from dead biofilm cells and their lysis into methane and other biogases. It can be confirmed that the effluent volatile suspended solids (VSSe) and effluent total suspended solids ( TSSe)/VSSe ratios were considerably lower for both of the bioreactors in study phase (A-5), as compared to other study phases. As an example, for the AnFFBR, the VSSe and VSSe/ TSSe ratios in study phase (A-5) were 48.4 mg day-1 and 0.537, respectively, which were considerably lower than the VSSe and VSSe/TSSe ratios of 224.4 mg day-1 and 0.671, respectively, obtained in study phase (A-1). This point illustrates that the remained tissue of dead bacteria was converted into other inert solids and final products including methane.

However, it should be noted that there are many studies which have reported the methane production rates between 0.4 and 0.47 L CH4/g COD [71-73].

Sludge quality

Another essential factor which is particularly important for land application is the concentration of DEP adsorbed by biomass, especially as it is known that the anaerobic digestion of biomass can produce residues that are rich in nutrients and have the potential to be used as the fertilizer [74]. According to the results, DEP can be considerably adsorbed by biofilm and, consequently, observed in effluent TSS. Moreover, the adsorbed concentration of DEP in sloughed biofilm depends on influent DEP concentration and organic loading rate. Minimum DEP concentrations in TSS for both the bioreactors were observed in study phase (A-5) which had the minimum organic loading rate and influent DEP concentration and were as 4.1 and 5.2 mg DEP g TSS-1 for the AnFFBR and UAnFFFBR, respectively. It should be noted that, for most of the study phases, the DEP concentrations in TSS were lower in the AnFFBR than UAnFFFBR, demonstrating that the sludge of the AnFFBR can be used more safely. This fact becomes more important when the produced sludge is used for land application and especially for agriculture purposes, because the previous research has shown that some phthalates can reduce the plant growth and cause chloroplast disintegrating [9]; others have stated that they can be introduced into the food chain and, consequently, produce human exposure. However, both the bioreactors' sludge should further be treated like other sludge. The performed research for evaluating digestion of phthalates in sludge has indicated that such phthalates,

particularly short-chains phthalates [e.g., DEP and dimethyl phthalate (DMP)], can be removed and mineralized in significant amounts. However, some other phthalates with long chains like di(2-ethylhexyl) phthalate (DEHP) may show resistance to digestion [6]. Although the produced sludge in the AnFFBR had less adsorbed DEP, the UAnFFFBR had sludge with less VSS/TSS ratio (effluent volatile suspended solids/effluent total suspended solids), which can make it more feasible for digestion and further uses. The minimum and maximum VSS/TSS ratios were observed as 0.537 and 0.671 for the AnFFBR, and 0.427 and 0.562 for the UAnFFFBR [both the minimum ratios were observed in study phase (A-5)]. The minimum VSS/ TSS ratio for the UAnFFFBR (0.427) was considerably less than the one reported for the aerobic moving bed biofilm reactor which ranged from about 0.52 to 0.68 [8] and can be considered as an important further sludge management advantage for this anaerobic reactor.

Metabolic pathways of DEP biodegradation

According to the observed by-products, DEP in both bio-reactors was primarily and mainly decomposed to mono-ethyl phthalate (MEP) by the de-alkylation of its first side chain and, then, the bioreaction continued to produce phthalic acid (PA) by removing the other side chain of MEP. This enzymatic reaction is known as de-esterifi-cation and commonly reported as the main phthalates' biodegradation route [75]. In addition, the trace amount of other center metabolites including dimethyl phthalate (DMP) and mono-methyl phthalate (MMP) can be the result of de-methylation which is a less common and predominant biodegradation pathway [75, 76]. Amir et al. [76] reported similar observation and stated that during sludge composting of di(2-ethylhexyl) phthalate (DEHP), DMP can be produced as the by-products of diethyl phthalate biodegradation due to de-methylation pathway. Similar observed metabolites for both bioreactors have indicated that the biodegradation pathway is not different and microbial activity in both bioreactors' condition can lead to the same by-products.

The other important degradation metabolites produced before ring cleavage included protocatechuic acid, 4-hydroxyphthalic acid, benzoic acid, 4,5-dihy-droxyphthalic acid, catechol, and 3,4-dihydroxy benzoic. Then, the biodegradation of the remained benzene ring of the mentioned by-products can result in the production of observed volatile fatty acids (VFAs), which have commonly been reported in the anaerobic biodegradation of organic compounds [77]. The final products of biodegradation were methane, carbon dioxide, dihy-drogen oxide, and hydrogen as mentioned before in the biogas. It should be noted that the VFAs and some of the other by-products can be produced in each of the other

biodegradation steps. Simplified biodegradation pathway is presented in Fig. 1.

Modeling and analysis of kinetic coefficients

As shown in Fig. 2, the first-order model for both of the bioreactors had low coefficients of determination, which were 0.377 and 0.305 for the AnFFBR and UAnFFFBR, respectively. Consequently, this model cannot be used for predicting performance of both of the bioreactors.

For computing saturation value constant (KB) and maximum substrate removal rate (Umax), Eq. (4) was plotted for both bioreactors in Fig. 3. Moreover, semi-strong values of coefficient of determination were obtained as 0.929 and 0.858 for the AnFFBR and UAnFFFBR, respectively, indicating that this model can be applied for designing and predicting the selected bioreactors. Higher Umax and KB of 4.04 (mg COD L-1 day-1) and 4.296 (mg L-1 day-1) for the AnFFBR compared to 2.404 (mg COD L-1 day-1) and 2.507 (mg L-1 day-1) for the UAnFFFBR, respectively, demonstrated that the microbial community had better conditions for biodegrading DEP and stabilizing COD in the AnFFBR than UAnFFFBR.

Equation (6) was plotted in Fig. 4 for computing KG coefficient. After calculating (m) and (n) coefficients, KG was determined using Eq. (7). Table 4 presents the Grau coefficient (KG) values. This model clearly had a high degree of precision with the coefficients of determination of 0.987 and 0.976 for the AnFFBR and UAnFFFBR,

respectively. In this regard, the Grau model can be used for predicting the performance of both of the bioreac-tors. Evaluating the KG values demonstrated this model to be in compliance with the obtained results, including COD, DEP, and TOC removal, because the KG values for the AnFFBR were higher than those for the UAnFFFBR. In study step (A), as HRT increased, the KG values in the reactors were significantly raised, which may be due to the bioavailability of microorganisms to the substrate. In study step (B), the KG values did not considerably change. With regard to the increasing DEP concentration and organic loading rate in study step (B) and stability of the calculated KG, it can be concluded that both of the bio-reactors can tolerate high organic loading rate and their performance was not affected by DEP concentration. This point can be the advantage of these anaerobic systems over aerobic or other treatment processes for industrial wastewater treatment [30].

According to Table 5 [correlation of methane yields (L CH4/g CODrem) and Grau coefficients], although in study step (B) for both bioreactors, strong correlation (R2 > 0.7) between methane production rates and KG values was not observed and, therefore, KG values could not be used for predicting methane yield with a high degree of precision, strong coefficients of determination (R2) of study step (A) as 0.791 and 0.937 for the UAnFFFBR and AnFFBR, respectively, can be used for the mentioned prediction. Positive slopes of the obtained equations in study

Fig. 1 Biological degradation pathway of diethyl phthalate

900 800 700 600 500 400 300 200 100 0

♦ AnFFBR ^ y = 2.7508x + 205.59 (Rz = 0.3773)

■ UAnFFFBR ^ y = 1.765x + 227.07 (Rz = 0.3056)

---- AnFFBR


♦ г ■ Bioreactor Ki(d-1)

AnFFBR 2.75

UAnFFFBR 1.765

Fig. 2 First-order model for the AnFFBR and UAnFFFBR

35 ♦ AnFFBR ^ y = 1.0634x + 0.2475 (Rz = 0.9298) UAnFFFBR ^ y = 1.0432x + 0.416 (Rz = 0.8588)

3 "-----AnFFBR

25 - -UAnFFFBR

o 2 -Ü2.

£ 1.5 ->

i -0.5 -0--

0 0.5 1 1.5 2 2.5


Fig. 3 Stover-Kincannon model for the AnFFBR and UAnFFFBR

Bioreactor Umax (mgCOD/L.d) Kb (mg/L.d)

AnFFBR 4.04 4.296

UAnFFFBR 2.404 2.507

♦ AnFFBR ^ У = 0.8913x + 0.363 (Rz = 0.9879) ■ UAnFFFBR ^ y = 0.8211x + 0.5362 (Rz = 0.9764)

---- AnFFBR


2 1.8 1.6 1.4 1.2 1 0.8 0.6 0.4 0.2 0

0 0.2 0.4 0.6 0.8 1 1.2 1.4 1.6

HRT (day)

. 4 Second-order (Grau) model for the AnFFBR and UAnFFFBR

Bioreactor m (d-1) n (dimentionless)

AnFFBR 0.3636 0.8913

UAnFFFBR 0.5362 0.8211

Table 5 Correlation of methane yields (L CH4/g CODrem) and Grau coefficients


Study phase A Y = 0.117* - 0.041 Y = 0.156* - 0.05

(R2 = 0.937) (R2 = 0.791)

Study phase B Y = -0.39* + 1.892 Y = -0.560* + 1.917

(R2 = 0.571) (R2 = 0.390)

step (A) (0.156 and 0.117 for the UAnFFFBR and AnFFBR, respectively) can demonstrate that HRT had a positive impact on both KG and methane yield values and these dependent variables changed with a similar direction. Higher calculated slope for the UAnFFFBR than AnFFBR can demonstrate that HRT had a more impact on methane production of the UAnFFFBR and, considering lower obtained methane yields for the UAnFFFBR, the AnFFBR had more advantages due to its less dependence on HRT. Furthermore, with regard to the insignificant variations of the calculated KG and negative slopes of study step (B) as 0.56 and 0.3G9 for the UAnFFFBR and AnFFBR, respectively, it can be concluded that methanogens activities were affected more than other non-methanogens by higher DEP loading rates. In addition, higher negative slope for the UAnFFFBR can demonstrate that this biore-actor was more sensitive to higher organic and DEP loading rates and its methane yield was more affected.

The kinetic coefficients including (KS) and (K) were computed using Eq. (10), which is plotted in Fig. 5. Half saturation constant and overall reaction rate were computed as 31.34 mg L-1 and 1.13 day-1 for the AnFFBR, and 24.87 mg L-1 and 1.03 day-1 for the UAnFFBR, respectively. It should be noted that the anaerobic waste-water treatment has slower substrate utilization rate as one of their disadvantages [31], which can be confirmed by the obtained overall reaction rates. As an example, for an aerated submerged fixed film reactor with glucose as the substrate, the overall reaction rate has been computed about 2.7 day-1, while for other conventional aerobic treatments, it is reported up to 12 day-1 [30, 49]. Half saturation constant presents the substrate concentration at half of the maximum substrate utilization rate. The previous studies have declared that a substrate with higher biodegradability has more KS [8]. Under low

Table 4 Grau second-order kinetic coefficients (KG) for AnFFBR and UAnFFFBR

Study phase A-1 A-2 A-3 A-4 A-5 B-1 B-2 B-3 B-4

KG for AnFFBR 2.038 2.779 3.253 3.64 3.822 3.92 3.92 4.023 4.01

KG for UAnFFFBR 1.522 2.112 2.300 2.524 2.723 2.875 2.828 2.875 2.909

0.016 -| 0.014 -0.012 -0.01 0.008 -0.006 0.004 0.002 0

AnFFBR ^ y = 0.0361X - 0.0319 {R2 = 0.9203) UAnFFFBR ^ y = 0.0417X - 0.0402 (R2 = 0.9686) AnFFBR UAnFFFBR

Bioreactor Ks (mg/L) K (d-1)

AnFFBR 31.34 1.131

UAnFFFBR 24.87 1.035

0.6 0.8 1 Xatt/(Q.(S0-S)

Fig. 5 Mathematical calculation of (KS) and (K) for the AnFFBR and UAnFFFBR S

substrate availabilities and at the same maximum specific growth rate in the two microbial groups or even two different bioreactors, smaller KS is obtained for microbial mass or a bioreactor with higher affinity to the substrate [78]. As the calculated maximum growth rate (^m) for the selected bioreactors is different, KS coefficients should not be directly applied to compare the bioreactors' performances. By contrast, if the difference between two values is ignored, the UAnFFFBR can be stated to have better performance under low substrate concentrations and, therefore, the microbial growth will be less influenced by substrate concentration. It is a very important point of view, especially for industries having different work shifts, as it does not continuously generate waste-water or the substrate concentration is not stable in influent wastewater.

The computed biomass yield coefficients (Y) plotted in Fig. 6 were obtained as 0.156 and 0.146 (g VS produced/g COD utilized). The (Y) coefficient is an important parameter for sludge management and its subsequent disposal [77, 79]. From this point of view and without other

removal performances, the UAnFFFBR can be preferably used when strong wastewater with high hydraulic loading rate must be treated. Under this circumstance, the insignificant difference observed for (Y) between the selected bioreactors results in sludge production to differ significantly. Meanwhile, the biomass yields for both bioreactors are very desirable when compared to those for aerobic based reactors. Experimental results for the biodegradation of similar compounds in aerobic biore-actors have shown that they have biomass yields more than 0.5 or even up to 0.78, which is more than five times greater than the calculated (Y) in this study [8, 48]. The Kd coefficients, which present specific decay rate (and are expressed as g VS of microbial mass loses/g VS of presents mass day or day-1) were computed as 0.107 and 0.100 day-1 for the AnFFBR and UAnFFFBR, respectively. The difference of the values, however, was observed to be insignificant between the bioreactors. These values placed in the typical range of 0.06-0.15 day-1 for the conventional treatment [30].

Consequently, maximum specific growth rate (day-1) values were calculated from Eq. (15) as 0.176 and 0.151 day-1 for the AnFFBR and UAnFFFBR, respectively. This coefficient had good compliance with other experimental results (including COD removal) and with the Stover-Kincannon and second-order (Grau) models, indicating that the AnFFBR had better performance. Considering the same inoculation used for both of the bioreactors, demonstrated that the microbial mass had more favorable conditions in the AnFFBR. Finally, (rsu) coefficients were computed for the study phases with maximum COD removal (phase B-4 for both of the bioreactors). Higher rsu of 938.11 mg L-1 day-1 for the AnFFBR compared to 842.26 mg L-1 day-1 for the UAnFFFBR confirmed other experimental results and proved the AnFFBR to have better performance due to its faster utilization rate.

AnFFBR ^ y = 0.685X + 6.4058 (R2 = 0.9884)

UAnFFFBR ^ y = 0.6877X + 6.8403 (R2 = 0.9912)



Bioreactor Y Kd (d-1)

AnFFBR 0.156 0.107

UAnFFFBR 0.146 0.100


Fig. 6 Mathematical calculation of (Y and (Kd) for the AnFFBR and UAnFFFBR


Both the bioreactors were capable of treating low-to-high strength DEP wastewater; however, the AnFFBR was preferred, since it could achieve 90.31% COD removal at 36 h of HRT. This study suggested the Grau and Stover-Kincannon models for predicting bioreactors due to their suitable coefficients of determination and good conformity of their kinetic parameters to the obtained results. Moreover, the experimental results and obtained kinetic coefficients indicated that the AnFFBR had better performance than the UAnFFFBR. Although both these bioreactors can achieve nearly 90% of DEP removal, they are promising only as pre-treatment methods and, due to standard regulations, should be coupled with further technology. Although from the sludge management

standpoint, the AnFFBR can be considered as a safer treatment option than the UAnFFFBR due to its lower DEP concentrations in sludge, the UAnFFFBR had less VSS/TSS ratio, which makes it more practical for digestion. Moreover, the AnFFBR had more sludge yield which should also be taken into consideration. Finally, both the bioreactors showed considerable methane yield; however, compared to the UAnFFFBR, the AnFFBR had more potential for bioenergy production, which could result in saving more energy and costs.

Additional file


DEP: diethyl phthalate; AnFFBR: anaerobic fixed film baffled reactor; UAnFFFBR: and up-flow anaerobic fixed film fixed bed reactor; CNEMC: China National Environmental Monitoring Center; USEPA: United State Environment Protection Agency; COD: chemical oxygen demand; MWTP: municipal wastewater treatment plant; HDPE: high-density polyethylene; HRT: hydraulic retention time; GC: gas chromatograph; FID: flame ionization detector; MS: mass spectrometer; TOC: total organic carbon; VS: volatile solids; TS: total solids; VSS: volatile suspended solids; TSS: total suspended solids; MEP: mono-ethyl phthalate; PA: phthalic acid; MMP: mono-methyl phthalate; DEHP: di(2-ethylhexyl) phthalate; DMP: dimethyl phthalate; VFAs: volatile fatty acids; STP: standard temperature and pressure; Lorg: organic loading rate; SRT: solid retention time.

Authors' contributions

EA, SY, MG, MA, and AA designed and performed the experiments. EA, SY, and HRG did the mathematical modeling and kinetic calculations. EA, SY, HRG, and KS drafted the manuscript. EA, MM, SR, HRG, and MG revised the manuscript and contributed in the improvement of paper. MM, SR, and EA prepared the additional file. SY and EA analyzed the data, prepared all the figures, and supervised the research. EA, AA, and MA operated the bioreactors and performed the sampling. EA, MG, SY, AA, and HRG wrote the manuscript. All authors read and approved the final manuscript.

Author details

1 Department of Environmental Health Engineering, Aradan School of Public Health and Paramedical, Semnan University of Medical Sciences, Semnan, Iran.

2 Department of Environmental Health Engineering, School of Public Health, Tehran University of Medical Sciences, Tehran, Iran. 3 Students' Scientific Research Center (SSRC), Tehran University of Medical Sciences, Tehran, Iran.

4 Department of Environmental Health Engineering, School of Public Health, Kashan University of Medical Sciences, Kashan, Iran. 5 Occupational Health Research Center (OHRC), Iran University of Medical Sciences, Tehran, Iran. 6 Department of Environmental Health Engineering, School of Public Health, Iran University of Medical Sciences, Tehran, Iran. 7 Department of Environmental Health Engineering, Faculty of Health, Hormozgan University of Medical Sciences, Bandar Abbas, Iran. 8 Environmental Science and Technology Research Center, Department of Environmental Health, School of Public Health, Shahid Sadoughi University of Medical Sciences, Yazd, Iran. 9 Department of Environmental Health, School of Public Health, Sabzevar University of Medical Sciences, Sabzevar, Iran. 10 Research Center for Environmental Determinants of Health (RCEDH), Kermanshah University of Medical Sciences, Kermanshah, Iran. 11 Social Determinants of Health Research Center, Yasuj University of Medical Sciences, Yasuj, Iran.


The authors are grateful to the financial support by Occupational Health Research Center (OHRC), Iran University of Medical Sciences. The authors would like to thank Prof. Alireza Mesdaghinia and Prof. Kamyar Yaghmaeian for fruitful discussions. The authors would like to express their

gratitude to Prof. Parvin Pasalar for her constant inspiration, encouragement, and support throughout the whole period of this project.

Competing interests

The authors declare that they have no competing interests. Availability of supporting data

All data generated or analyzed during this study are included in this published article with Additional file 1: Figure S1.


This work was financially supported by grants from Occupational Health Research Center (OHRC), Iran University of Medical Sciences (Grant Number 92-03-132-23829).

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Received: 8 January 2017 Accepted: 22 May 2017 Published online: 31 May 2017


1. Chen M, Wang P, Chen Q, Wu J, Chen X. A clustering algorithm for sample data based on environmental pollution characteristics. Atmos Environ. 2015;107:194-203.

2. Pirsaheb M, Khosravi T, Sharafi K, Babajani L, Rezaei M. Measurement of heavy metals concentration in drinking water from source to consumption site in Kermanshah—Iran. World Appl Sci J. 2013;21(3):416-23.

3. Atharizade M, Miranzadeh MB. Evaluation of efficacy of advanced oxidation processes fenton, fenton-like and photo-fenton for removal of phenol from aqueous solutions. J Chem Soc Pak. 2015;37(02):266.

4. Azari A, Kakavandi B, Kalantary R, Ahmadi E, Gholami M, Torkshavand Z, Azizi M. Rapid and efficient magnetically removal of heavy metals by magnetite-activated carbon composite: a statistical design approach. J Porous Mater. 2015;22(4):1083-96.

5. Farzadkia M, Gholami M, Kermani M, Yaghmaeian K. Biosorption of hexavalent chromium from aqueous solutions by chemically modified brown algae of Sargassum sp. and dried activated sludge. Asian J Chem. 2012;24(11):5257.

6. Abdel daiem MM, Rivera-Utrilla J, Ocampo-Pérez R, Méndez-Díaz JD, Sánchez-Polo M. Environmental impact of phthalic acid esters and their removal from water and sediments by different technologies—a review. J Environ Manag. 2012;109:164-78.

7. Thomaidi VS, Stasinakis AS, Borova VL, Thomaidis NS. Is there a risk for the aquatic environment due to the existence of emerging organic contaminants in treated domestic wastewater? Greece as a case-study. J Hazard Mater. 2015;283:740-7.

8. Ahmadi E, Gholami M, Farzadkia M, Nabizadeh R, Azari A. Study of moving bed biofilm reactor in diethyl phthalate and diallyl phthalate removal from synthetic wastewater. Bioresour Technol. 2015;183:129-35.

9. Wu Y, Si Y, Zhou D, Gao J. Adsorption of diethyl phthalate ester to clay minerals. Chemosphere. 2015;119:690-6.

10. Zhou X, Zhang Y, Huang Z, Lu D, Zhu A, Shi G. Ionic liquids modified gra-phene oxide composites: a high efficient adsorbent for phthalates from aqueous solution. Sci Rep. 2016;6:38417.

11. Prasad B, Suresh S. Biodegradation of dimethyl phthalate ester using free cells, entrapped cells of Variovorax sp. BS1 and cell free enzyme extracts: a comparative study. Int Biodeterior Biodegrad. 2015;97:179-87.

12. Benjamin S, Pradeep S, Josh MS, Kumar S, Masai E. A monograph on the remediation of hazardous phthalates. J Hazard Mater. 2015;298:58-72.

13. Li B, Xu X, Zhu Y, Cao J, Zhang Y, Huo X. Neonatal phthalate ester exposure induced placental MTs, FATP1 and HFABP mRNA expression in two districts of southeast China. Sci Rep. 2016;6:21004.

14. Ma P, Liu X, Wu J, Yan B, Zhang Y, Lu Y, Wu Y, Liu C, Guo J, Nanberg E, et al. Cognitive deficits and anxiety induced by diisononyl phthalate in mice and the neuroprotective effects of melatonin. Sci Rep. 2015;5:14676.

Additional file 1: Figure S1. Schematic diagram of the (A) UAnFFFBR; and (B) AnFFBR.

Whangsuk W, Sungkeeree P, Nakasiri M, Thiengmag S, Mongkolsuk S, Loprasert S. Two endocrine disrupting dibutyl phthalate degrading esterases and their compensatory gene expression in Sphingobium sp. SM42. Int Biodeterior Biodegrad. 2015;99:45-54.

Tsai H-J, Wu C-F, Tsai Y-C, Huang P-C, Chen M-L, Wang S-L, Chen B-H, Chen C-C, Wu W-C, Hsu P-S, et al. Intake of phthalate-tainted foods and serum thyroid hormones in Taiwanese children and adolescents. Sci Rep. 2016;6:30589.

Tan W, Zhang Y, He X, Xi B, Gao R, Mao X, Huang C, Zhang H, Li D, Liang Q, et al. Distribution patterns of phthalic acid esters in soil particle-size fractions determine biouptake in soil-cereal crop systems. Sci Rep. 2016;6:31987.

Liu X, Shi J, Bo T, Li H, Crittenden JC. Occurrence and risk assessment of selected phthalates in drinking water from waterworks in China. Environ Sci Pollut Res. 2015;22(14):10690-8.

Wang C, Yang L, Wang S, Zhang Z, Yu Y, Wang M, Cromie M, Gao W, Wang S-L. The classic EDCs, phthalate esters and organochlorines, in relation to abnormal sperm quality: a systematic review with meta-analysis. Sci Rep. 2016;6:19982.

Lea RG, Byers AS, Sumner RN, Rhind SM, Zhang Z, Freeman SL, Moxon R, Richardson HM, Green M, Craigon J, et al. Environmental chemicals impact dog semen quality in vitro and may be associated with a temporal decline in sperm motility and increased cryptorchidism. Sci Rep. 2016;6:31281.

Nerin C, Ubeda J, Alfaro P, Dahmani Y, Aznar M, Canellas E, Ausejo R Compounds from multilayer plastic bags cause reproductive failures in artificial insemination. Sci Rep. 2014;4:4913.

Shapiro GD, Dodds L, Arbuckle TE, Ashley-Martin J, Fraser W, Fisher M, Taback S, Keely E, Bouchard MF, Monnier P, et al. Exposure to phthalates, bisphenol A and metals in pregnancy and the association with impaired glucose tolerance and gestational diabetes mellitus: the MIREC study. Environ Int. 2015;83:63-71.

Chen X, Zhang X, Yang Y, Yue D, Xiao L, Yang L. Biodegradation of an endocrine-disrupting chemical di-n-butyl phthalate by newly isolated Camelimonas sp. and enzymatic properties of its hydrolase. Biodegradation. 2015;26(2):171-82.

Dehghani S, Jafari AJ, Farzadkia M, Gholami M. Sulfonamide antibiotic reduction in aquatic environment by application of fenton oxidation process. Iran J Environ Health Sci Eng. 2013;10(1):1-5. Mohammadi H, Sabzali A, Gholami M, Dehghanifard E, Mirzaei R Comparative study of SMBR and extended aeration activated sludge processes in the treatment of high-strength wastewaters. Desalination. 2012;287:109-15.

Shahamat YD, Farzadkia M, Nasseri S, Mahvi AH, Gholami M, Esrafili A. Magnetic heterogeneous catalytic ozonation: a new removal method for phenol in industrial wastewater. J Environ Health Sci Eng. 2014;12(1):1-12. Kotlarska E, tuczkiewicz A, Pisowacka M, Burzynski A. Antibiotic resistance and prevalence of class 1 and 2 integrons in Escherichia coli isolated from two wastewater treatment plants, and their receiving waters (Gulf of Gdansk, Baltic Sea, Poland). Environ Sci Pollut Res. 2015;22(3):2018-30. Boboescu IZ, Ilie M, Gherman VD, Mirel I, Pap B, Negrea A, Kondorosi E, Biro T, Maroti G. Revealing the factors influencing a fermentative biohy-drogen production process using industrial wastewater as fermentation substrate. Biotechnol Biofuels. 2014;7(1):1.

Guven D, Hanhan O, Aksoy EC, Insel G, Cokgor E. Impact of paint shop decanter effluents on biological treatability of automotive industry wastewater. J Hazard Mater. 2017;330:61-7.

Tchobanoglous G, Stensel H, Tsuchihashi R, Burton F, Abu-Orf M, Bowden G, Pfrang W, Metcalf and Eddy I AECOM. Wastewater engineering: treatment and resource recovery. In: New York: McGraw-Hill Book Company; 2014.

Chan YJ, Chong MF, Law CL, Hassell DG. A review on anaerobic-aerobic treatment of industrial and municipal wastewater. Chem Eng J. 2009;155(1-2):1-18.

Marques IP, Gil L, La Cara F. Energetic and biochemical valorization of cork boiling wastewater by anaerobic digestion. Biotechnol Biofuels. 2014;7(1):1.

Liu W, He Z, Yang C, Zhou A, Guo Z, Liang B, Varrone C, Wang A-J. Microbial network for waste activated sludge cascade utilization in an integrated system of microbial electrolysis and anaerobic fermentation. Biotechnol Biofuels. 2016;9(1):1.

España-Gamboa EI, Mijangos-Cortés JO, Hernández-Zárate G, Maldonado JAD, Alzate-Gaviria LM. Methane production by treating vinasses from hydrous ethanol using a modified UASB reactor. Biotechnol Biofuels. 2012;5(1):1.

Oosterkamp MJ, Méndez-García C, Kim C-H, Bauer S, Ibáñez AB, Zimmerman S, Hong P-Y, Cann IK, Mackie RI. Lignocellulose-derived thin stillage composition and efficient biological treatment with a high-rate hybrid anaerobic bioreactor system. Biotechnol Biofuels. 2016;9(1):1. Bakonyi P, Nemestóthy N, Simon V, Bélafi-Bakó K. Review on the start-up experiences of continuous fermentative hydrogen producing bioreactors. Renew Sustain Energy Rev. 2014;40:806-13. Sivagurunathan P, Sen B, Lin C-Y. High-rate fermentative hydrogen production from beverage wastewater. Appl Energy. 2015;147:1-9. Kumar G, Sivagurunathan P, Park J-H, Park J-H, Park H-D, Yoon J-J, Kim S-H. HRT dependent performance and bacterial community population of granular hydrogen-producing mixed cultures fed with galactose. Biore-sour Technol. 2016;206:188-94.

Sivagurunathan P, Anburajan P, Kumar G, Kim S-H. Effect of hydraulic retention time (HRT) on biohydrogen production from galactose in an up-flow anaerobic sludge blanket reactor. Int J Hydrogen Energy. 2016;41(46):21670-7.

Limkhuansuwan V, Wangnai C, Chaiprasert P Evaluation of performance characteristics in an anaerobic baffled-fixed film reactor (ABFFR) treating palm oil mill effluent. J Biotechnol. 2010;150(Supplement):267. Kocadagistan B, Kocadagistan E, Topcu N, Demircioglu N. Wastewater treatment with combined upflow anaerobic fixed-bed and suspended aerobic reactor equipped with a membrane unit. Process Biochem. 2005;40(1):177-82.

Mohanty A, Yadav AK, Chaudhury GR Removal of nitrate nitrogen and chemical oxygen demand in upflow fixed bed reactor using heterotrophic microorganisms. J Clean Prod. 2016;127:573-8. Yetilmezsoy K. Integration of kinetic modeling and desirability function approach for multi-objective optimization of UASB reactor treating poultry manure wastewater. Bioresour Technol. 2012;118:89-101. Leyva-Díaz JC, González-Martínez A, González-López J, Muñío MM, Poyatos JM. Kinetic modeling and microbiological study of two-step nitrification in a membrane bioreactor and hybrid moving bed biofilm reactor-membrane bioreactor for wastewater treatment. Chem Eng J. 2015;259:692-702.

Leyva-Díaz JC, Calderón K, Rodríguez FA, González-López J, Hontoria E, Poyatos JM. Comparative kinetic study between moving bed biofilm reactor-membrane bioreactor and membrane bioreactor systems and their influence on organic matter and nutrients removal. Biochem Eng J. 2013;77:28-40.

Tebbe CC, Dohrmann AB, Hemkemeyer M, Nather A. Microbial Community Profiling:SSCP and T-RFLP Techniques. In: McGenity TJ, Timmis KN, Nogales B, editors. Hydrocarbon and Lipid Microbiology Protocols: Microbial Quantitation, Community Profiling and Array Approaches. Berlin, Heidelberg: Springer; 2017. p. 101-26. Astudillo-Garcia C, Bell JJ, Webster NS, Glasl B, Jompa J, Montoya JM, Taylor MW. Evaluating the core microbiota in complex communities: a systematic investigation. Environ Microbiol. 2017;19(4):1450-62. Pirsaheb M, Mesdaghinia A-R, Shahtaheri SJ, Zinatizadeh AA. Kinetic evaluation and process performance of a fixed film bioreactor removing phthalic acid and dimethyl phthalate. J Hazard Mater. 2009;167(1-3):500-6.

Nabizadeh R, Mesdaghinia A. Simulation of microbial mass and its variation in biofilm systems using STELLA. J Chem Technol Biotechnol. 2006;81(7):1209-17.

Srinivasan G, Subramaniam R, Kumar N. Kinetic evaluation of fixed film fixed bed anaerobic reactor by using dairy wastewater. Am Eurasian J Sci Res. 2009;4(3):213-8.

Ahmadi E, Yousefzadeh S, Ansari M, Ghaffari HR, Azari A, Miri M, Mesdaghinia A, Nabizadeh R, Kakavandi B, Ahmadi P, et al. Performance, kinetic, and biodegradation pathway evaluation of anaerobic fixed film fixed bed reactor in removing phthalic acid esters from wastewater. Sci Rep. 2017;7:41020.

Kishida N, Kim J, Tsuneda S, Sudo R. Anaerobic/oxic/anoxic granular sludge process as an effective nutrient removal process utilizing denitrifying polyphosphate-accumulating organisms. Water Res. 2006;40(12):2303-10.

53. Borghei SM, Sharbatmaleki M, Pourrezaie P, Borghei G. Kinetics of organic removal in fixed-bed aerobic biological reactor. Bioresour Technol. 2008;99(5):1118-24.

54. Rice EW, Bridgewater L, Association APH. Standard methods for the examination of water and wastewater. Washington, DC: American Public Health Association; 2012.

55. Lay J-J, Lee Y-J, Noike T. Feasibility of biological hydrogen production from organic fraction of municipal solid waste. Water Res. 1999;33(11):2579-86.

56. Pirsaheb M, Rostamifar M, Mansouri AM, Zinatizadeh AAL, Sharafi K. Performance of an anaerobic baffled reactor (ABR) treating high strength baker's yeast manufacturing wastewater. J Taiwan Inst Chem Eng. 2015;47:137-48.

57. Juretic D, Kusic H, Dionysiou DD, Rasulev B, Peternel I, Bozic AL. Prediction of key structural features responsible for aromaticity of single-benzene ring pollutants and their photooxidative intermediates. Chem Eng J. 2015;276:261-73.

58. Zhu G-F, Li J-Z, Wu P, Jin H-Z, Wang Z. The performance and phase separated characteristics of an anaerobic baffled reactor treating soybean protein processing wastewater. Bioresour Technol. 2008;99(17):8027-33.

59. Farzadkia M, Rezaei Kalantary R, Mousavi G, Jorfi S, Gholami M. The effect of organic loading on propylene glycol removal using fixed bed activated sludge hybrid reactor. Chem Biochem Eng Q. 2010;24(2):227-34.

60. Zhang T, Huang Z, Chen X, Huang M, Ruan J. Degradation behavior of dimethyl phthalate in an anaerobic/anoxic/oxic system. J Environ Manag. 2016;184(2):281-8.

61. Oliver R, May E, Williams J. The occurrence and removal of phthalates in a trickle filter STW. Water Res. 2005;39(18):4436-44.

62. Gao D, Li Z, Wen Z, Ren N. Occurrence and fate of phthalate esters in full-scale domestic wastewater treatment plants and their impact on receiving waters along the Songhua River in China. Chemosphere. 2014;95:24-32.

63. Liang D-W, Zhang T, Fang HHP Anaerobic degradation of dimethyl phthalate in wastewater in a UASB reactor. Water Res. 2007;41(13):2879-84.

64. Huang M-Z, Ma Y-W, Wang Y, Wan J-Q, Zhang H-P The fate of di-n-butyl phthalate in a laboratory-scale anaerobic/anoxic/oxic wastewater treatment process. Bioresour Technol. 2010;101(20):7767-72.

65. Zolfaghari M, Drogui P, Seyhi B, Brar SK, Buelna G, Dube R. Occurrence, fate and effects of di(2-ethylhexyl) phthalate in wastewater treatment plants: a review. Environ Pollut. 2014;194:281-93.

66. De Vrieze J, Regueiro L, Props R, Vilchez-Vargas R, Jauregui R, Pieper DH, Lema JM, Carballa M. Presence does not imply activity: DNA and RNA patterns differ in response to salt perturbation in anaerobic digestion. Biotechnol Biofuels. 2016;9(1):244.

67. Wang Y, Zhang Y, Wang J, Meng L. Effects of volatile fatty acid concentrations on methane yield and methanogenic bacteria. Biomass Bioenergy. 2009;33(5):848-53.

68. Maya-Altamira L, Baun A, Angelidaki I, Schmidt JE. Influence of wastewater characteristics on methane potential in food-processing industry wastewaters. Water Res. 2008;42(8-9):2195-203.

69. Fuchs W, Binder H, Mavrias G, Braun R. Anaerobic treatment of wastewater with high organic content using a stirred tank reactor coupled with a membrane filtration unit. Water Res. 2003;37(4):902-8.

70. Ghasimi DS, Tao Y, Kreuk M, Zandvoort MH, Lier JB. Microbial population dynamics during long-term sludge adaptation of thermophilic and mesophilic sequencing batch digesters treating sewage fine sieved fraction at varying organic loading rates. Biotechnol Biofuels. 2015;8(1):1.

71. Cavaleiro AJ, Pereira MA, Alves M. Enhancement of methane production from long chain fatty acid based effluents. Bioresour Technol. 2008;99(10):4086-95.

72. Kujgu ÖS, Sponza DT. Performance of anaerobic baffled reactor (ABR) treating synthetic wastewater containing p-nitrophenol. Enzyme Microb Technol. 2005;36(7):888-95.

73. Faisal M, Unno H. Kinetic analysis of palm oil mill wastewater treatment by a modified anaerobic baffled reactor. Biochem Eng J. 2001;9(1):25-31.

74. Sun L, Liu T, Müller B, Schnürer A. The microbial community structure in industrial biogas plants influences the degradation rate of straw and cellulose in batch tests. Biotechnol Biofuels. 2016;9(1):128.

75. Liang D-W, Zhang T, Fang HH, He J. Phthalates biodegradation in the environment. Appl Microbiol Biotechnol. 2008;80(2):183-98.

76. Amir S, Hafidi M, Merlina G, Hamdi H, Jouraiphy A, El Gharous M, Revel JC. Fate of phthalic acid esters during composting of both lagooning and activated sludges. Process Biochem. 2005;40(6):2183-90.

77. Tchobanoglous G, Stensel H, Tsuchihashi R, Burton F, Abu-Orf M, Bowden G, Pfrang W. Wastewater engineering: treatment and resource recovery. 5th ed. New York: McGraw-Hill Book Company; 2014.

78. Arnaldos M, Amerlinck Y, Rehman U, Maere T, Van Hoey S, Naessens W, Nopens I. From the affinity constant to the half-saturation index: understanding conventional modeling concepts in novel wastewater treatment processes. Water Res. 2015;70:458-70.

79. Zhou A, Zhang J, Wen K, Liu Z, Wang G, Liu W, Wang A, Yue X. What could the entire cornstover contribute to the enhancement of waste activated sludge acidification? Performance assessment and microbial community analysis. Biotechnol Biofuels. 2016;9(1):241.

Submit your next manuscript to BioMed Central and we will help you at every step:

• We accept pre-submission inquiries

• Our selector tool helps you to find the most relevant journal

• We provide round the clock customer support

• Convenient online submission

• Thorough peer review

• Inclusion in PubMed and all major indexing services

• Maximum visibility for your research

Submit your manuscript at * Central