Accepted Manuscript Optimizing and Real-time Control of Biofilm Formation, Growth and Renewal in Denitrifying Biofilter Xiuhong Liu, Hongchen Wang, Feng Long, Lu Qi, Haitao Fan PII: DOI: Reference:
S0960-8524(16)30232-2 http://dx.doi.org/10.1016/j.biortech.2016.02.095 BITE 16145
To appear in:
Bioresource Technology
Received Date: Revised Date: Accepted Date:
20 November 2015 19 February 2016 22 February 2016
Please cite this article as: Liu, X., Wang, H., Long, F., Qi, L., Fan, H., Optimizing and Real-time Control of Biofilm Formation, Growth and Renewal in Denitrifying Biofilter, Bioresource Technology (2016), doi: http://dx.doi.org/ 10.1016/j.biortech.2016.02.095
This is a PDF file of an unedited manuscript that has been accepted for publication. As a service to our customers we are providing this early version of the manuscript. The manuscript will undergo copyediting, typesetting, and review of the resulting proof before it is published in its final form. Please note that during the production process errors may be discovered which could affect the content, and all legal disclaimers that apply to the journal pertain.
1
2
Optimizing and Real-time Control of Biofilm
3
Formation, Growth and Renewal in Denitrifying
4
Biofilter
5 6
Xiuhong Liu, Hongchen Wang *, Feng Long, Lu Qi, Haitao Fan
7
8 9
School of Environment & Natural Resources, Renmin University of China, Beijing, China,100872;
10
11
Corresponding author: Hongchen Wang
12
Phone: + 86-10-62510853
13
E-mail:
[email protected]
14 15
1
16
ABSTRACT
17
A pilot-scale denitrifying biofilter (DNBF) with a treatment capacity of 600 m3/d was used
18
to study real-time control of biofilm formation, removal and renewal. The results showed
19
biofilm formation, growth and removal can be well controlled using on-line monitored
20
turbidity. The status of filter layer condition can be well indicated by Turb break points on
21
turbidity profile. There was a very good linear relationship between biofilm growth degree
22
(Xbiof) and filter clogging degree (Cfilter) with R2 higher than 0.99. Filter layer clogging
23
coefficient (Yc) lower than 0.27 can be used to determine stable filter layer condition. Since
24
variations of turbidity during backwash well fitted normal distribution with R2 higher than
25
0.96, biofilm removal during backwash also can be well optimized by turbidity. Although
26
biofilm structure and nirK-coding denitrifying communities using different carbon sources
27
were much more different, DNBF was still successfully and stably optimized and real-time
28
controlled via on-line turbidity.
29
30 31
KEY WORDS
32
Biofilm; formation; turbidity; real-time control; denitrifying biofilter
33
2
34
1. Introduction
35
Wastewater reclamation is one of the most effective ways to solve both water shortage and
36
pollution problems. Many wastewater treatment plants in water deficient regions, especially
37
large and medium-sized cities, have been up-graded to produce reclaimed water, and
38
advanced nitrate removal is required. Because denitrifying biofilter (DNBF) has advantages
39
of high treatment efficiency, small footprint, and low shock loading impact, it is widely
40
used for advanced nitrate removal. However, because biofilm growth cannot be easily
41
monitored, many biofilters do not provide stable performance. Moreover, some biofilters
42
can not be normally operated due to frequent clogging (Leverenz et al., 2009; Snowball,
43
2006) (Figure S1), which limits the engineering application and development of biofilter.
44
In fixed-film wastewater treatment systems (Corona et al., 2013; Ji et al., 2014; Wang et
45
al., 2015), biofilm is used to remove pollutants in wastewater; whereas, in other fields
46
including membrane treatment processes (Miura et al., 2006; Pang et al., 2005), drinking
47
water or reclaimed water distribution systems (Mathieu et al., 2014; Yang et al., 2015), and
48
food processing industry (Simões et al., 2010), biofilm formation is harmful and should be
49
avoided. It was reported that biofilm formation comprised a series of complicated steps
50
(Giaouris et al., 2014; Simões et al., 2010) including cell deposition, cell adsorption and/or
51
desorption, cell–cell signaling and EPS production, replication and growth, secretion of
52
biofilm matrix, and biofilm detachment or sloughing. To date, there is no effective method
53
to accurately measure or determine biofilm growth condition in a biofilter, especially in a
54
DNBF. Generally, during the start-up stage of a DNBF, biofilm was naturally cultivated
55
with nitrate and organic carbon. However, since biofilm can not be practically monitored or
56
measured during treatment, it is hard to determine the biofilter run time between
3
57
backwashes (backwash cycle), the backwash strength and duration. This may result in over
58
backwash and energy cost, or biofilter clogging if backwash is not performed in time. In
59
the last decade, most studies on biofilter process were focused on optimizing backwash
60
procedures (Amburgey and Amirtharajah, 2005; Slavik et al., 2013; Snowball, 2006), while
61
limited studies have been conducted on optimizing the biofilter operation based on the
62
extent of biofilm growth and removal.
63
Many on-line biofilm monitoring methods, such as differential turbidity measurement
64
(Métadier and Bertrand-Krajewski, 2012), microwaves (Saber et al., 2013), multi-channel
65
impedimetric and amperometric sensor (Pires et al., 2013) and thermal sensors (Reyes-
66
Romero et al., 2014), have been used to monitor biofilm growth. Among these methods,
67
turbidity is suitable for industrial applications. It was used to calculate urban storm water
68
pollutant concentrations and loads (Métadier and Bertrand-Krajewski, 2012). Backwash is
69
the key operational step to control biofilm growth, removal and renewal in treatment plants
70
that employ biofilter/filter as major unit process. However, the effective method to monitor
71
biofilm growth on-line in biofilter has not been well developed, and very limited studies
72
have been conducted on real-time control and optimization of biofilter-based biofilm
73
formation and growth.
74
Therefore, this study aimed to: a) determine the relationship between biofilm formation
75
and turbidity variation during the biofilm cultivation stage of a DNBF; b) control and
76
optimize biofilm growth, removal and renewal during filtration and backwash processes in
77
the DNBF via on-line monitored turbidity; c) test the stability of real-time control of
78
biofilm growth, removal and renewal in the DNBF using different typical carbon sources.
4
79
2. Materials and methods
80
2.1 Pilot-scale and lab-scale DNBF
81
A pilot-scale DNBF located at Beijing Gaobeidian Municipal Wastewater Treatment Plant
82
(WWTP) with a maximum treatment capacity of 600 m3/d was used in this study. The
83
DNBF was packed with expanded clay particles (4-6 mm) at a bed depth of 2.5 m (Table
84
S1). The pilot-scale treatment system consisted of a DNBF reactor, a carbon dosage system,
85
and a backwash system (Figure S2).
86
The control consists of four parts: detectors, a computer, interface cards and control units.
87
Nitrate (NO3-), turbidity, pH, dissolved oxygen (DO) and pressure sensors were installed in
88
the top section and the bottom of DNBF reactor, respectively. The values of the NO3-,
89
turbidity, pH, and DO were recorded every 0.5–10 minute and then transferred to a
90
programmable logic controller (PLC) and process control system (PCS). PCS was
91
programmed according to the control logic.
92
A lab-scale DNBF with a total volume of 30 L and filter media height of 1.1m was used
93
to study biofilm formation further. Except for volume and media height, the lab-scale
94
DNBF was the same as the pilot-scale DNBF.
95
2.2 Secondary effluent composition
96
The secondary effluent of Gaobeidian municipal WWTP that employs an anoxic/oxic (A/O)
97
process was continuously fed to the DNBF. The concentrations of the soluble chemical
98
oxygen demand (SCOD), NH4+-N, NO3--N, turbidity, and SS in the secondary effluent were
99
in the range of 23.28-76.03 mg/L, 0.33-2.36 mg/L, 20.40-34.12 mg/L, 0.48-13.10 NTU,
100
and 1.33-16.15 mg/L, respectively.
101
5
102
2.3 Experimental design
103
Experiments were conducted for 210 days, which was divided into 3 phases. Experimental
104
procedures, parameters, aims, and backwash operation were shown in Table 1 and tableS2.
105
In the first phase, the natural biofilm cultivation method was used in the pilot-scale and lab-
106
scale DNBF. Since biofilm formation in the pilot-scale DNBF was only operated for one
107
time and the filter media is difficultly taken out to observation, the lab-scale DNBF was
108
operated with Filtration velocity (FV) of 2.38m/h only to further confirm the results
109
obtained in the pilot-scale DNBF. In the second phase, DNBF was firstly operated for 40
110
days to establish backwash strength. Thereafter, variations of the effluent turbidity during
111
long-term filtration and backwash operation were studied on the following 38 days and 33
112
days, respectively. In the third phase, the stability of long-term operation of DNBF with
113
real-time control using turbidity as parameter was tested using two typical carbon sources.
114
When carbon source changed from sodium acetate to methanol, nitrate removal efficiency
115
of DNBF was recovered for 30 days.
116 117
2.4 Analysis
118
COD, NH4+-N, NO3--N, NO2--N, TN, PO43--P, and total phosphorus (TP) were measured
119
according to the standard methods given in APHA (APHA, 1998). DO, turbidity, NO3--N
120
and pressure were measured on-line using oxygen, turbidity, NO3--N meters (LDO,
121
SOLITAX sc, NITRATAX plus sc, Hach Company, USA) and pressure sensor,
122
respectively.
123
Extracellular Polymeric Substances (EPS) was extracted using formaldehyde-NaOH
124
according to the method described by Liu and Fang (Liu and Fang, 2002). Pre-treatment of
6
125
filter media for SEM observation was carried out using the modified method described by
126
Wang et al. (2013) (Wang et al., 2013). Visualization of the samples was conducted using a
127
Scanning Electron Microscope (Hoskin Scientific, Tokyo, Japan).
128
2.5 DNA extraction, PCR, cloning, sequencing and phylogenetic analysis
129
Filter media samples were taken from DNBF on day 100 and 203. Biofilm was removed
130
from filter media by shaker. The collected samples were freeze-dried in Labconco Freezone
131
1 L (Labconco, USA) and stored at −20 °C.
132
0.05–0.10 g of dry sludge sample was taken to extract genomic DNA using a FastDNA
133
SPIN Kit for soil (Qiagen, CA, USA). 1ul of extracted DNA was taken to measure DNA
134
concentration by Nanodrop Spectrophotometer (ND-1000, Thermo Fisher Scientific, USA).
135
Fragments of nirK genes were amplified with the primer set FlaCu/R3Cu (473 bp) as
136
described by Hallin et al. (1999) (Hallin and Lindgren, 1999).
137
The PCR was performed in a C1000TM thermal cycler (BioRad, USA). Cycle conditions
138
for the amplification were as follows: 2 min at 94°C; 35 cycles with each cycle consisting
139
of 30 sec at 94°C, 1 min at 57°C, and 1 min at 72°C; followed by a final 10-min extension
140
at 72°C. PCR products were visualized on 1.5% (w/v) agarose gel electrophoresis to
141
confirm the product size, and then, purified by Wizard® SV Gel and PCR Clean-Up
142
System (Promega, Madison, USA). The detailed information on cloning and sequencing
143
were described by Gao et al. (2014) (Gao et al., 2014).
144
All the sequences and their reference sequences obtained from NCBI BLAST were
145
aligned using MEGA 6.0 software (Tamura et al., 2013). The sequences sharing 98%
146
similarity were grouped into the same operational taxonomic unit (OTU) using software
7
147
mother (Schloss et al., 2009). Phylogenetic trees were generated by neighbor-joining (NJ)
148
with the Jukes–Cantor correction in MEGA (Tamura et al., 2013).
149
2.6 Calculations
150
Biofilm growth degree (Xbiof), filter clogging degree (Cfilter) and filter layer clogging
151
coefficient (Yc), were calculated based on the biomass yield coefficient and total effluent
152
turbidity as described in equation(1) – (3) , respectively. NT
{
Y NO3 − N ∑ [NO3 − N ]in - [NO3 − N ]out 153
}
n ∆t
× ∆t × Q n∆t
n =1
X biofilm =
(1)
Vfilter
154
where, YNO3-N is sludge yield coefficient (0.74 gVSS/gNO3-N), [NO3-N]in and [NO3-N]out
155
are nitrate concentrations in the influent and effluent (g/m3), n is data points on the profiles
156
of on-line monitored nitrate, T is the time of on-line monitored nitrate, NT is the data points
157
at the time of T, ∆t is the on-line monitored interval (day), Q is the influent flow rate (m3/d),
158
Vfilter is the volume of filter media (m3). NT
∑ {[Turb] - [Turb ] } in
159
160
C filter =
out n ∆t
× ∆t × Q n∆t
(2)
n =1
Vfilter
where, [Turb]in and [Turb]out are turbidity concentrations in the influent and effluent(NTU). NT
161
Yc =
∑ {[Turb] - [Turb] } in
1 YNO3 − N
out n∆t
n =1 NT
∑ {[NO
3
− N ]in - [NO3 − N ]out
n =1
8
}
n∆t
(3)
162
where, Yc is the slope of Cfilter/Xbiof.
163 164
3. Results and discussion
165
3.1 Correlation between biofilm formation and effluent turbidity during start-up of
166
DNBF
167
Figure 1 showed variations of COD, NO3--N, turbidity and pressure during biofilm
168
cultivation in pilot-scale DNBF. Because almost no denitrifying bacteria was cultivated in
169
the first 2 days, very small amount of nitrate was removed and the added carbon source
170
(CODadd=CODcarbon-CODin) was almost not consumed. From the 3rd day, nitrate removal
171
efficiency increased to 55%, and CODadd was still not sufficiently utilized. Meanwhile,
172
gases emitted from the effluent, indicating that denitrification occurred. On the following 2
173
days, CODadd was sufficiently utilized with removal efficiency higher than 98%. Figure S3
174
showed variations of CODremoved, NO3--Nremoved, turbidity and pressure during biofilm
175
cultivation in the lab-scale DNBF. In the first two days, variations of the removed COD,
176
nitrate and turbidity in lab-scale DNBF were very similar to these in the pilot-scale DNBF.
177
From day 3 to 4, with COD and nitrate slightly removed, turbidity gradually increased.
178
SEM results also found cells absorbed and grew on the rough surface of filter media
179
(Figure S4). Some of cells aggregated together and formed micro colonies. These results
180
indicated that with the cells replication and growth, since some discohesive cells and
181
produced EPS were desorbed from filter media into water (Busscher et al., 2010) under the
182
shear forces of the moving water with the filter media surface, turbidity were gradually
183
increased. On the flowing four days, with the increase of COD and nitrate removal,
184
turbidity was still gradually increased. On day 8, SEM results also demonstrated the
9
185
relatively integrated and stable biofilm were formed. Meanwhile, because the loose
186
structure biofilm were detached from filter media due to the hydraulic shear force, small
187
amount of macroscopic biofilm fragments were found in the effluent. After that, since the
188
excessive growth of biofilm caused the filter layer blocked, with the matured and aged
189
biofilm slough off and new biofilm formation, the hydraulic sheer force fluctuated caused
190
turbidity fluctuation. Therefore, based on the above results in pilot-scale and lab-scale
191
DNBF, turbidity in the effluent can well indicated biofilm formation and growth. To date,
192
few studies used turbidity to monitor biofilm formation and growth in a wastewater
193
treatment system.
194
After the biofilm density gradually increased during the cultivation stage, backwash was
195
performed to scour the excessive biomass under higher hydrodynamic shear force, to renew
196
the biofilm with new denitrifying bacteria (Table S2). In Phase 2, on the first 40 days, if
197
backwash strength was controlled too low (Figure S5), during filtration, not only the
198
effluent nitrate was unstable even after adding sufficient carbon source, but also turbidity
199
was fluctuated significantly. Some large SS aggregates were also observed in the effluent
200
(Figure S6). Furthermore, because of more frequent filter clogging, backwash operation
201
became more difficult. However, high backwash intensity might also result in filter media
202
loss (Figure S7), and mix the supporting layer with filter media layer (Figure S8) in the
203
reclaimed WWTP. Therefore, appropriate backwash is the key to ensure long-term
204
operation of the biofilter.
205
Since biofilter is a closed system, the biofilm on the biofilter media can not be measured
206
directly. In addition, biofilter clogging has a time-lag property, which even increases the
10
207
difficulty to operate and maintain a biofilter. These undesirable backwash control methods
208
mainly resulted from not timely monitoring biofilm growth and filter layer clog condition.
209
In wastewater/water treatment plants, biofilters are normally operated based on operator's
210
experience, and these undesirable backwash control methods often occur, which slow the
211
application and development of biofilter. To date, there is no effective way to on-line
212
monitoring biofilm growth degree and determine filter media clog degree in fundamental
213
and applied research.
214
3.2 Correlation between effluent turbidity and biofilm renewal and removal during
215
normal operation
216
Based on the obtained results during biofilm formation, turbidity might also be used to
217
indicate biofilm growth and removal during filtration and backwash operation. To validate
218
this hypothesis, turbidity and nitrate concentrations in the effluent were monitored online
219
during filtration and backwash processes.
220
3.2.1 Variation in effluent turbidity indicating biofilm growth and filter layer condition
221
during filtration
222
Figure 2a showed typical variations of the effluent turbidity, CODadd, influent nitrate, and
223
effluent nitrate during long-term filtration. With the biodegradation of COD and the
224
reduction of nitrate, biofilm grew gradually and became thicker. Meanwhile, filter layer
225
was gradually blocked. At the beginning of filtration, the effluent nitrate varied with
226
CODadd and the effluent turbidity increase slightly. However, after continuous filtration for
227
about 49 h, since some microcolonies and/or biofilm fragments were sloughed off from the
11
228
surface of biofilm, both the effluent turbidity and nitrate increased significantly. Turb
229
break-1 and Nitrate break-1 points both appeared on the turbidity and ef-nitrate profiles,
230
respectively. Thereafter, the effluent nitrate still varied with CODadd, while the effluent
231
turbidity turned increase to decrease because of biosorption and low nitrate removal amount;
232
whereas, in the 63.25 h, because of higher hydraulic shear force caused by more serious
233
filter layer blockage, the effluent turbidity increased sharply again, and Turb break-2 point
234
appeared on the effluent turbidity. Subsequently, more and more filter layers were blocked
235
due to biofilm growth, and as a result, both the effluent turbidity and nitrate fluctuated
236
significantly. At the same time, both Turb shock and nitrate shock points appeared on the
237
turbidity and effluent nitrate profiles, respectively. Therefore, based on the above results,
238
Turb break points on turbidity profile were suitable for indicating the status of filter layer
239
condition, which can be divided into four phases, including stable, intergraded, clogged and
240
damaged phases.
241
Figure 2b showed relationships between Xbiof and Cfilter. It was found that Cfilter increased
242
with the increase of Xbiof. In each phase, there was a very good linear relationship between
243
Xbiof and Cfilter with R2 higher than 0.99. Yc kept constant in each phase, while, Yc
244
increased with the clogging status of biofilter transferred from one phase to the next.
245
Especially, Yc increased sharply from 0.54 in the clogging phase and to 2.46 in the
246
damaged phase. These results indicate that Yc can be used to determine the status of filter
247
layer clogging degree.
248
Therefore, during normal operation, if backwash cycle was controlled on-line between
249
Turb break-1 point and Turb break-2 point, not only DNBF kept stable operation and
250
avoided filter layer clogging, but also carbon dosage can be well controlled by the on-line
12
251
monitored effluent nitrate. Even when turbidity can not be monitored on-line for a long
252
time, filter layer clogging degree still can be determined by calculating Yc through
253
monitoring the effluent nitrate and turbidity for the final 2-3 h of filtration.
254
3.2.2 Effluent turbidity indicating the biofilm removal degree during backwash
255
Figure 3 showed the variations in the effluent turbidity and nitrate concentration during
256
filtration and backwash processes. During the filtration, when the Turb break-1 point lasted
257
for 30 min, the biofilter transformed from the stable phase to the intergraded phase. Then
258
the filtration was stopped and backwash was performed. During the air backwash, turbidity
259
increased significantly. Hereafter, during air + water backwash, most biofilm was detached
260
and carried out by backwash water. When the biofilm almost was completely detached,
261
turbidity reached a peak and then decreased; turbidity apex points appeared on the turbidity
262
profile. Thereafter, during water backwash, backwash water carried the residual detached
263
biofilm out of filter layer. Meanwhile, turbidity decreased sharply, followed by a gradual
264
decrease. After that, turbidity kept constant, which indicated that the detached biofilm was
265
completely removed by backwash water. Among three phases of backwash, air + water
266
backwash was the most intensive and decisive phase to keep the equilibrium between the
267
excessive biofilm growth and biofilm removal, in which, more than half of the detached
268
biofilm was removed. Furthermore, during air + water backwash, Turb apex indicated that
269
the loose biofilm was almost completely detached. If air + water backwash was
270
continuously performed, the tight formation of biofilm would be removed, which would
271
lead to the recovery of nitrate removal efficiency and biofilm renewal in the next filtration
272
operation. Therefore, Turb apex points and Turb platform states appeared on turbidity
13
273
profiles can be used to accurately determine air+water backwash time and water backwash
274
time, respectively.
275
Variations of turbidity during backwash well fitted normal distribution (R2=0.96313,
276
Figure 3). Backwash process was very similar to chromatographic separation, in which, air
277
and water, the mobile phase, carried the detached biofilm out of the filter media. Air mainly
278
enhanced filter media scour, while water flushed and carried the detached biofilm out of
279
biofilter. Hydrodynamic shear force was an effective tool to control of biofilm structures
280
(Morgenroth and Wilderer, 2000). However, in an engineering sense, as there is no
281
effective and simple method to measure the biofilm growth, the biomass removal, and the
282
strength of hydrodynamic shear force, some biofilters cannot be kept under stable condition
283
(Snowball, 2006). Furthermore, some filters had to remove filter media due to the extensive
284
clog (Figure S1). The results clearly showed that turb peak height, turbidity increase rate
285
and decrease rate indicated the strength of hydraulic shear force caused by air+water
286
velocity, air velocity and water velocity, respectively. Furthermore, turb peak area indicated
287
the removed biofilm amount. Therefore, biofilm removal during backwash process can be
288
well optimized using turbidity as control parameter.
289
3.3 Real-time control of the biofilm growth and renewal in DNBF using the effluent
290
turbidity
291
Preliminary studies found that the effluent turbidity could effectively indicate the biofilm
292
growth, renewal and removal in the DBNF, which could be used to real-time control
293
filtration and backwash for DNBF. Therefore, in Phase 3, the DNBF was operated for a
294
long period of time using the effluent turbidity to control backwash. In this test, two
14
295
widely-used carbon sources, sodium acetate and methanol, were used in DNBF,
296
respectively.
297
3.3.1 Filter performance
298
Figure 4 and Figure 5 showed the variations in the effluent nitrate and turbidity during
299
filtrations and backwashes with real-time control using sodium acetate and methanol as
300
carbon sources, respectively. During the filtrations, turbidity gradually increased and nitrate
301
gradually decreased, which indicated that biofilm gradually grew and filter layer were in
302
the stable stage. After Turb break-1 point appeared on turbidity profiles for 1.5-2.5 h
303
(extend time), the filtration was stopped. As backwash cycles were well controlled between
304
Turb break-1 point and Turb break-2 point, filter layer were well controlled at the
305
intergraded stages with Yc lower than 0.27. Turbidity and Yc during filtration in DNBF
306
using sodium acetate as carbon source were little higher than that using methanol. During
307
backwash, air+water backwash time and water backwash time was controlled at Turb apex
308
points and Turb platform points or turbidity lower than 2.5 NTU points. Each backwash
309
did not affect nitrate removal efficiency in the following filtration, which indicated
310
backwashes were well controlled without excessive backwash.
311
3.3.2 Biofilm structure and diversity of denitrifying bacteria
312
Biofilm structure plays an important role in both overall performance of filter. Figure S9
313
showed SEM images of the biofilm structure when sodium acetate (SA-biofilm) and
314
methanol (M-biofilm) were used as carbon sources. SA-biofilm and M-biofilm were both
315
composed by organisms and the biofilm matrix. The compact SA-biofilm with massive
15
316
distinct porous and channels completely covered the surface of filter media, indicating that
317
SA-biofilm was not only continuous and uneven, but also complicate and thick. Some
318
microcolonies with loose structure attached on the biofilm matrix. However, M-biofilm was
319
relatively simple and has clear structure. The distribution of M-biofilm was not only
320
discontinuous and uneven, but also thinner than SA-biofilm. Clearly, the bottom layer of
321
M-biofilm was single and regular shape bacillus, which was coated and covered by EPS.
322
SA-biofilm contained a relatively higher EPS than M-biofilm, further indicating SA-
323
biofilm was thicker than M-biofilm. The cell cluster with loose structure was also attached
324
to the biofilm matrix. The phylogenetic analysis of nirK gene sequences showed that nirK-
325
encoding denitrifiers in M-biofilm and SA-biofilm had high biodiversity (Figure 6);
326
whereas, nirK-encoding denitrifiers in M-biofilm were much more different from that in
327
SA-biofilm. Mesorhizobium and Bradyrhizobium were dominant denitrifier species in SA-
328
biofilm, while, Hyphomicrobium and Paracoccus were key denitrifying populations in M-
329
biofilm.
330
The differences of biofilm structure and diversity of nirK-encoding denitrifiers between
331
M-biofilm and SA-biofilm directly resulted in different variations of turbidity profiles using
332
sodium acetate and methanol as carbon sources. Although Turb break-1 points were both
333
appeared on turbidity profiles whether using sodium acetate or methanol as carbon sources
334
during filtration, because thick SA-biofilm has relatively more and larger microcolonies
335
with loose structure than M-biofilm, at the final stage of filtration, turbidity in the effluent
336
of DNBF when using sodium acetate as carbon source was fluctuated, and a higher value
337
than that using methanol as carbon source was observed. During the following backwash
16
338
process, turb-apex points also appeared on turbidity profiles using both sodium acetate and
339
methanol as carbon sources. Furthermore, because difference biomass growth rate caused
340
by different denitrifying communities using sodium acetate and methanol as carbon sources
341
(Gómez et al., 2000), turbidity at turb-apex point was relatively higher when sodium acetate
342
was used. Thus, because real-time control strategy via turbidity variation well indicated
343
biofilm growth and removal degree, no matter what carbon source was used, characteristic
344
points were appeared on turbidity profiles during filtration and backwash, which can be
345
used to control the filter layer and prevent excessive backwash.
346
In biofilm technologies for treating wastewater, biofilm is too thick or matured, which
347
would cause bulk biofilm slough out of carrier in the moving bed process (Huang et al.,
348
2014) or biofilter clog in the fixed bed process. However, biofilm is too thin or immature,
349
which would cause the decrease of treatment efficiency. Thus, it is difficult to control of
350
biofilm growth in biofilm system. This study suggests that the proposed real-time control
351
method can not only indicate biofilm formation and growth degree, but also well control
352
and optimize biofilm removal and renewal. More importantly, this real-time control method
353
via on-line turbidity monitoring might also be used in other fields, such as, the moving-bed
354
process, rotating biological contactor, granular system and water pipe system, which should
355
be investigated further.
356 357 358
4. Conclusions
The main conclusions obtained in the pilot-scale DNBF (600m3/d):
17
359
Biofilm formation, growth and removal can be well detected and controlled using on-
360
line monitored turbidity during biofilm cultivation, filtration and backwash processes
361
in DNBF.
362
The status of filter layer condition can be divided into stable, intergraded, clogged and
363
damaged phases based on Turb break points on turbidity profile. Yc can be used to
364
determine filter layer condition in DNBF.
365
Although biofilm structure and nirK-coding denitrifying communities using different
366
carbon sources were much more different, DNBF was still successfully and stably
367
optimized and real-time controlled via on-line turbidity.
368
Acknowledgements
369
This research was supported by National Science and Technology Major Projects on
370
Water pollution Control and Treatment (2011ZX07316-001 and 2013ZX07314-001) and
371
National Natural Science Foundation of China (51508561). The authors are grateful to
372
Jianmin Wang and Guoqiang Liu from Department of Civil, Architectural and
373
Environmental Engineering, Missouri University of Science and Technology, USA, for
374
providing language help and giving helpful advices.
375 376
Supplementary materials
377
Filter media pictures of the serious clogged DNBF in wastewater reclaimed plant (figure
378
S1), parameters, materials and schematic diagram of the pilot-scale DNBF (Table S1 and
18
379
Figure S2), correlation between biofilm formation and effluent turbidity during start-up of
380
lab-scale DNBF(Figure S3-4), backwash operation during the entire experiments (Table
381
S2), effects of low and higher backwash strength on water quality and biofilter
382
operation(Figure S5-8), and biofilm structure in DNBF using sodium acetate and methanol
383
as carbon sources(Figure S9)
384
References:
385
1. Amburgey, J.E., Amirtharajah, A., 2005. Strategic filter backwashing techniques and
386
resulting particle passage. Journal of environmental engineering-asce, 131, 535-547.
387
2. APHA, 1998. Standard Methods for Examination of Water and Wastewater. American
388
Public Health Association, Washington, DC.
389
3. Busscher, H.J., Norde, W., Sharma, P.K., van der Mei, H.C., 2010. Interfacial re-
390
arrangement in initial microbial adhesion to surfaces. Current Opinion in Colloid &
391
Interface Science, 15, 510-517.
392
4.
Corona, F., Mulas, M., Haimi, H., Sundell, L., Heinonen, M., Vahala, R., 2013.
393
Monitoring nitrate concentrations in the denitrifying post-filtration unit of a municipal
394
wastewater treatment plant. Journal of Process Control, 23, 158-170.
395
5. Gao, J., Luo, X., Wu, G., Li, T., Peng, Y., 2014. Abundance and diversity based on
396
amoA genes of ammonia-oxidizing archaea and bacteria in ten wastewater treatment
397
systems. Applied Microbiology and Biotechnology, 98, 3339-3354.
398
6. Giaouris, E., Heir, E., Hébraud, M., Chorianopoulos, N., Langsrud, S., Møretrø, T.,
399
Habimana, O., Desvaux, M., Renier, S., Nychas, G., 2014. Attachment and biofilm
400
formation by foodborne bacteria in meat processing environments: Causes, implications,
19
401
role of bacterial interactions and control by alternative novel methods. Meat Science
402
Advancing Beef Safety through Research and Innovation: Prosafebeef, 97, 298-309.
403
7. Gómez, M.A., González-López, J., Hontoria-Garcı́ a, E., 2000. Influence of carbon
404
source on nitrate removal of contaminated groundwater in a denitrifying submerged
405
filter. Journal of Hazardous Materials, 80, 69-80.
406 407
8. Hallin, S., Lindgren, P.E., 1999. PCR detection of genes encoding nitrile reductase in denitrifying bacteria. Applied and Environmental Microbiology, 65, 1652-1657.
408
9. Huang, H., Ren, H., Ding, L., Geng, J., Xu, K., Zhang, Y., 2014. Aging biofilm from a
409
full-scale moving bed biofilm reactor:Characterization and enzymatic treatment study.
410
Bioresource Technology, 154, 122-130.
411 412 413 414 415 416
10. Ji, B., Wang, H., Yang, K., 2014. Nitrate and COD removal in an upflow biofilter under an aerobic atmosphere. Bioresource Technology, 158, 156-160. 12. Leverenz, H.L., Tchobanoglous, G., Darby, J.L., 2009. Clogging in intermittently dosed sand filters used for wastewater treatment. Water Research, 43, 695-705. 13. Liu, H., Fang, H.H.P., 2002. Extraction of extracellular polymeric substances (EPS) of sludges. Journal of Biotechnology, 95, 249-256.
417
14. Mathieu, L., Bertrand, I., Abe, Y., Angel, E., Block, J.C., Skali-Lami, S., Francius, G.,
418
2014. Drinking water biofilm cohesiveness changes under chlorination or
419
hydrodynamic stress. Water Research, 55, 175-184.
420
15. Métadier, M., Bertrand-Krajewski, J.L., 2012. The use of long-term on-line turbidity
421
measurements for the calculation of urban stormwater pollutant concentrations, loads,
422
pollutographs and intra-event fluxes. Water Research, 46, 6836-6856.
20
423
16. Miura, Y., Watanabe, Y., Okabe, S., 2006. Membrane Biofouling in Pilot-Scale
424
Membrane Bioreactors (MBRs) Treating Municipal Wastewater: Impact of Biofilm
425
Formation. Environmental Science & Technology, 41, 632-638.
426 427
17. Morgenroth, E., Wilderer, P.A., 2000. Influence of detachment mechanisms on competition in biofilms. Water Research, 34, 417-426.
428
18. Pang, C.M., Hong, P., Guo, H., Liu, W., 2005. Biofilm Formation Characteristics of
429
Bacterial Isolates Retrieved from a Reverse Osmosis Membrane. Environmental
430
Science & Technology, 39, 7541-7550.
431
19. Pires, L., Sachsenheimer, K., Kleintschek, T., Waldbaur, A., Schwartz, T., Rapp, B.E.,
432
2013. Online monitoring of biofilm growth and activity using a combined multi-
433
channel impedimetric and amperometric sensor. Biosensors & Bioelectronics, 47, 157-
434
163.
435 436
20. Reyes-Romero, D.F., Behrmann, O., Dame, G., Urban, G.A., 2014. Dynamic thermal sensor for biofilm monitoring. Sensors and Actuators A: Physical, 213, 43-51.
437
21. Saber, N., Ju, Y., Hsu, H., Lee, S., 2013. A feasibility study on the application of
438
microwaves for online biofilm monitoring in the pipelines. International Journal of
439
Pressure Vessels and Piping, 111, 99-105.
440
22. Schloss, P.D., Westcott, S.L., Ryabin, T., Hall, J.R., Hartmann, M., Hollister, E.B.,
441
Lesniewski, R.A., Oakley, B.B., Parks, D.H., Robinson, C.J., Sahl, J.W., Stres, B.,
442
Thallinger, G.G., Van Horn, D.J., Weber, C.F., 2009. Introducing mothur: Open-Source,
443
Platform-Independent, Community-Supported Software for Describing and Comparing
444
Microbial Communities. Applied and Environmental Microbiology, 75, 7537-7541.
21
445 446
23. Simões, M., Simões, L.C., Vieira, M.J., 2010. A review of current and emergent biofilm control strategies. LWT - Food Science and Technology, 43, 573-583.
447
24. Slavik, I., Jehmlich, A., Uhl, W., 2013. Impact of backwashing procedures on deep bed
448
filtration productivity in drinking water treatment. Water Research, 47, 6348-6357.
449
25. Snowball, M., 2006. Water treatment: Reducing backwash with air scouring. Filtration
450
& Separation, 43, 39-40.
451
26. Tamura, K., Stecher, G., Peterson, D., Filipski, A., Kumar, S., 2013. MEGA6:
452
Molecular Evolutionary Genetics Analysis Version 6.0. Molecular Biology and
453
Evolution, 30, 2725-2729.
454
27. Wang, H., Ji, G., Bai, X., 2015. Enhanced long-term ammonium removal and its ranked
455
contribution of microbial genes associated with nitrogen cycling in a lab-scale
456
multimedia biofilter. Bioresource Technology, 196, 57-64.
457
28. Wang, Y., Guo, G., Wang, H., Stephenson, T., Guo, J., Ye, L., 2013. Long-term impact
458
of anaerobic reaction time on the performance and granular characteristics of granular
459
denitrifying biological phosphorus removal systems. Water Research, 47, 5326-5337.
460
29. Yang, G., Feng, L., Wang, S., Zhou, J., Guo, C., Xia, T., Sun, W., Jiang, Y., Sun, X.,
461
Cao, L., Xu, X., Zhu, L., 2015. Potential risk and control strategy of biofilm
462
pretreatment process treating raw water. Bioresource Technology, 198, 456-463.
463
22
Filtration
COD/(mg/L)
100
DNBFin
DNBFcarbon
DNBFout
(A)
75 50 25 0
0
1
2
3
4
5
6
(B)
20 10 3 0
(C)
Turbidity/(NTU)
Pressure
10 8
2
6
1
4 2
0 0
1
2
3
4
5
6
Pressure/(kpa)
NO3--N/(mg/L)
30
0
Time/(days)
Figure 1. Variation of COD, NO3--N, turbidity and pressure concentrations during the biofilm cultivation stage of DNBF
Figure 2. Typical variations of the effluent turbidity, CODadd, influent nitrate, and effluent nitrate during long-term filtration (a) and relationships between Xbiof and Cfilter (b) in DNBF
Backwash
Filtration Filtration Air Air+Water
Water
900 800 700
500
Y=12.58 + 897.5*exp(-2*((X-2047)/6.3)^2)
Turb Apex
600
9
500 400
6
300 200
400
3
2036 2038 2040 2042 2044 2046 2048 2050 2052 2054 2056 2058 2060 2062 2064
0
30 25 20
Turb Platform Nitrate Break
100
300 200
15
0
10
Time(min)
100
Turbidity Break
0 0
35
6
12
18 Time(h)
-
600
12
Turb Fit Curve
-
Turb(NTU)
700 Turb(NTU)
Adj.R-Square=0.96313
15
40
24
30
Figure 3. Typical variations of turbidity during filtration and backwash in DNBF
NO3 -N(mg/L)
800
-
NO3 -N
NO3 -N(mg/L)
900
5 0 36
NO3--NOUT
NO3--NIN
1000
Yc=0.26
Yc=0.18
Yc=0.25
Turbidity
35 Yc=0.21
Yc=0.17
Yc=0.20
30 800 25 20 20 15 10
Turb Break-1
Turb Break-1 Turb Break-1
Turb Break-1
Turb Break-1
Turb Break-1
NO3--N(mg/L)
Turbidity(NTU)
600
15 10
5
5
0
0 0
1
2
3
4 5 6 Time (Days)
7
8
9
Figure 4. Variations in the effluent nitrate and turbidity during filtrations and backwashes with real-time control using sodium acetate as carbon sources
-
Yc=0.21
Yc=0.13
35
Yc=0.11 Yc=0.12
Yc=0.11 Yc=0.13
800
30
600
25
400
20 15
15
-
Turbidity(NTU)
1000
Turbidity
NO3 -NOUT
NO3 -N(mg/L)
-
NO3 -NIN
Turb
Turb Break-1
10
Turb Break-1
Turb
Turb
Break-1
Break-1
Break-1
Turb Break-1
10
5
5
0
0 11
0
1
2
3
4
5
6
7
8
9
10
Time(days) Figure 5. Variations in the effluent nitrate and turbidity during filtrations and backwashes with real-time control using methanol as carbon source
OTU5 (9/27) OTU3 (1/29) OTU5 (9/27) Bradyrhizobium sp. D189(AB480442.1) Bradyrhizobium sp. D203a(AB480454.1) OTU3 (1/27) OTU16 (1/27) Bradyrhizobium japonicum USDA 110 DNA(BA000040.2) Bradyrhizobium sp. BTAi1(CP000494.1) OTU10 (1/27) Bosea sp. D257c(AB480483.1) Bosea sp.(HQ916684.1) 100 OTU11 (1/27) 100 Rhodanobacter sp. D206a(AB480456.1) Ochrobactrum sp. clone 1-40(GU136454.1) 100 Ochrobactrum sp. clone 2-80(GU136458.1) OTU15 (1/27) OTU12 (1/27) OTU1 (2/27) OTU9 (1/27) OTU6 (1/27) Mesorhizobium opportunistum WSM2075(CP002279.1) Mesorhizobium sp. GSM-373(FN600572.1) OTU9 (1/29) OTU2 (1/27) OTU14 (1/27) Chelativorans sp. BNC1(CP000390.1) OTU13 (1/27) 69 Paracoccus sp. R-26823(AM230847.1) 63 Paracoccus sp. R-28242(AM230886.1) 100 OTU5 (1/29) Paracoccus sp. R-26824(AM230857.1) OTU14 (1/29) OTU1 (1/29) OTU11 (1/29) OTU4 (10/29) 100 Rhodobacter sphaeroides forma sp.(AJ224908.1) Rhodobacter sphaeroides(U62291.1) OTU4 (3/27) Rhizobium gallicum(CP006880.1) OTU20 (1/27) Rhizobium sp. PY13(DQ096645.1) Rhizobium sp. R-24663(AM230832.1) Hyphomicrobium zavarzinii IFAM ZV-622(AJ224902.1) Hyphomicrobium nitrativorans NL23(CP006912.1) OTU14 (1/29) OTU2 (2/29) OTU7 (3/29) OTU13 (2/29) OTU8 (2/29) OTU6 (1/29) OTU10 (2/29)
100 47 46 49 67 68 49 27
90
42
13
81 13
74 25
4 15
100 46
19 29 23
24
75
100 39
92
16
100 52 60
7
100 69 67
100 62
75 67
96 55 47
Bradyrhizobium
Bosea Rhodanobacter Ochrobactrum
Mesorhizobium
Chelativorans
Paracoccus
Rhodobacter
Rhizobium
Hyphomicrobium
0.05
Figure. 6 NJ phylogenetic tree of the nirK -containing bacteria using methanol and sodium acetate as carbon sources. Sequences of nirK genes using methanol and sodium acetate as carbon sources are shown with "▲ OTU" and "♦OTU", respectively.
Table 1. Experimental procedure and aims Phases Duration Operation and main Parameters Aims (days) 1st 6/15* CS: sodium acetate Determine correlation between biofilm formation and turbidity during the FV:4.76m/h; FV:2.38m/h* start-up of the pilot-scale and lab-scale DNBF nd 2 9-120 CS: sodium acetate; FV: 4.76 m/h (40) Fixed time backwash cycle(48h) Effects of backwash on DNBF operation and establish backwash strength (38) Long-term filtration: Three times Determine the possibility of using turbidity as control parameter to indicate (Each time: filtration 5-6 days and biofilm growth and filter layer condition during long-term filtration. recovery 5days) Establish the methods of determining filter layer clogging degree. (33) Backwash cycle: Turb Break-1 point The correlation between the variations of the effluent turbidity and backwash (FV: 4.76-8.79m/h) rd 3 121-210 FV=8.79m/h Establish the real-time control approach of biofilm growth, removal and (30) CS: sodium acetate;; Real-time control renewal in DNBF via on-line monitoring turbidity. CS: methanol; Test the stability of the real-time control approach using two typical carbon (30) recovery stage sources. (30) Real-time control Analyze biofilm structure and diversity of denitrifying bacteria *in the lab-scale DNBF; CS—Carbon source; Real-time control—backwash was controlled via turbidity;
1
464 465
Optimizing and control of biofilm was studied in biofilter (600 m3/d).
466
Biofilm formation, growth and removal can be controlled using turbidity.
467
The filter layer status can be indicated by Turb break points on turbidity profile.
468
Filter layer clogging coefficient can be used to determine filter layer condition.
469 470
23