A Pilot-scale Benthic Microbial Electrochemical System (BMES) for Enhanced Organic Removal in Sediment Restoration

A benthic microbial electrochemical systems (BMES) of 195 L (120 cm long, 25 cm wide and 65 cm height) was constructed for sediment organic removal. Sediment from a natural river (Ashi River) was used as test sediments in the present research. Three-dimensional anode (Tri-DSA) with honeycomb structure composed of carbon cloth and supporting skeleton was employed in this research for the first time. The results demonstrated that BMES performed good in organic-matter degradation and energy generation from sediment and could be considered for river sediments in situ restoration as novel method. Community analysis from the soil and anode using 16S rDNA gene sequencing showed that more electrogenic functional bacteria was accumulated in anode area when circuit connected than control system.

Scientific RepoRts | 7:39802 | DOI: 10.1038/srep39802 developed with time to not only remove organics, but also for the production and recovery of value-added products from substrates 10 . Another, the merging of phototrophic organisms or plants into microbial fuel cells (MFCs) is an interesting option since they can act as efficient in situ oxygenators, thus facilitating the cathodic reaction of microbial fuel cells 9 . Experiments have also shown that through stimulating the microbial electrogenic metabolism, dibenzothiophene removal was enhanced by more than 3-fold compared to the natural attenuation 11 . On the other hand, Zhang et al. arranged anodes with two different ways for enhancing the bioremediation of contaminated soil, up to 12.5% of the total petroleum hydrocarbon (TPH) was removed in reactors with anodes horizontally arranged after 135 days, which was 95.3% higher than that in the disconnected control (6.4%) 8 . Yang et al. built a 100 L sediment microbial fuel cells (SMFC) inoculate with heavily contaminated sediments, and the total organic chemical degradation efficiency was 22.1% in the electricity generating SMFCs, which is significantly higher than that in the open-circuited SMFC (3.8%) after two years' long-term applicability without external electron donor addition 8 . It is worth noticing that in the past few years, considerable progresses have been made in MFC research, but significant challenges still remaining in river sediments restoration 10,12 .
In this study, a total volume of 195 L BMES was constructed and was expected to be employed for sediment organic removal. The aim of the present study was to supply a potential novel method to serve as an in-situ river sediment restoration.

Results and Discussion
TOC and TN variation in the sediment. The initial total nitrogen 13 and the TOC of the sediment was 3.1 and 33.7 g kg −1 in dry sediment base (Table S1). For the BMES reactor, the TOC was decreased by 5.0% during the first 15 days of operation, and then the TOC removal was 5.6% during the second 15 days, a relative lower removal of 3.6% was obtained during the third 15 days, while it only decreased by 0.3% during the last 15 d (Total decreased by 14.5%) (Fig. 1). The BMES kept a high TOC removal during the first 45 days before the removal rate sharply declined (there was no obvious change during the last 15 days) (Table S2). During the whole operation, the TOC removal efficiency of BMES (14.5%) was 1.2-and 6.9-fold in comparison with the S Control (11.6% removal) and W Control (only 2.1% removal). It indicated that the flushing method (W Control), which was usually used for sediment restoration has little effects on the organics removal. From the experimental data obtained here, it is obviously that BMES demonstrated good capability of simultaneous organic matter removal as reported 14,15 .
TN decreased by 3.2% during the first 15 days of the BMES operation, and then decreased by 1.6% during the second 15 days, a removal of 12.7% was obtained during the third 15 days, while it only decreased by 1.0% during the last 15 d (Total decreased by 18.5%). But, during the same period, there were no significant changes in S Control (only 1.9% removal) and W Control (1.0% removal) during the whole operation (Table S3). The TN removal rate of the BMES was 8.7− and 18.2− fold than that of S Control and W Control.
PAHs removal performances. In total 12 kinds of PAHs were detected in the present sediments ( Table 1).
The concentrations of both Benzo(a) pyrene (BaP) and Benzo (k) fluoranthene (BkF) (five-ringed PAHs) was higher than 12 mg/kg, accounting for 60.35% of total PAHs (TPAHs) in the initial sediments (Fig. 2). It is probably because low molecular weight PAHs, containing two-to three-ringed PAHs were less toxic to microbes in the soil and can serve as a carbon source involved in the microbial metabolisms and accumulated more in the sediment 8 . Also high molecular weight PAHs, containing four-to six-ringed PAHs were hard to be decomposed by microbes.
The total contents of BaP during the 60 days' operation were decreased by 50%, 37.4% and 30.8% in BMES, S Control and W Control respectively and BkF decreased by 50%, 28% and 21.8% in BMES, S Control and W Control. More important, BMES has higher PAHs removal efficiencies than that the S Control and W Control (1.4 fold PAH removal than S Control and 1.8 fold than W Control). Complex organic compounds, such as PAHs decomposition might need a long time leading to inhibition of TOC further removal in the sediment 16 (Fig. 1). Sorption of PAHs on natural sorbents, like sands, sandy loam soils, and silt loam soils are usually normal process and is regarded as dynamic fast step compared to the biodegradation process 17 . Former research also showed that dissolved organic matter has an positive effect of PAHs sorption and induced PAHs contents are usually proportion to the dissolved organic matter 18 . PAHs in the present system was expected to be attached on the anode surface firstly and then decomposed by the anode bacteria in BMES. But in BMES, with the microbes enrichment on the anode, PAH sorption/desorption hysteresis declined, due to EPS production of anode bacteria 19 . It is obvious that the degradation was accelerated in BMES owing to the promotion by the current generated in BMES. It can also be concluded that washing flushing has no extra effects on the PAHs removal. The attached PAHs on anode was then further degraded by the electrogensis bactetria.  Table 1. The content and ratio of different PAHs in the initial sediments from Ash River.

Figure 2. PAHs change in the sediment (BMES stands for the sediment in the BMES, S Control stands for the sediment in the S Control reactor, W Control stands for the sediment in the W Control reactor).
Scientific RepoRts | 7:39802 | DOI: 10.1038/srep39802 TOC and TN changes in the water layer of BMES. Total 135 L tap water (TOC 1.9 mg/L) was used for the water layer, the TOC contained in the sediment began to be flushed and released into the water layer. The initial TOC concentration in water was 60 ~70 mg L −1 , yet was reduced to < 30 mg L −1 in the first 3 days in the three research systems (Fig. 3). Resettling of pollutants to sediments 20,21 might be the main reason for the declination of TOC in the first 3 days, yet some released soluble organics from sediments remained in the water layers. In the followed 5 d to 25 d, the TOC of the BMES and S Control had similar trend. After 60 days' operation, the TOC of the BMES was approximately half than that of the S Control. The results indicated that the deployment of electrode as electron acceptor in BMES was helpful to better removal of organic pollutants than natural biodegradation in S control. The lowest TOC contained in water layer was achieved by W Control which was washed by fresh water in a continuous flow rate of 400 L h −1 .
After 135 L tap water (TN 1.9 mg/L) was used for the water layer, the initial concentration of TN in the water was 39 to 40 mg L −1 for the flushing and also the desorption. During the first 3 days, TN content in BMES and W Control dropped sharply, and in the later 30 days experiments, TN removal rates in BMES and W Control slowed down and kept at 7.47~13.3 mg L −1 and 1.1 ~3.2 mg L −1 respectively. In the first 3 days, TN in S control was at 39-45 mg L −1 and then at 17-28 mg L −1 in the following experimental period.
DO, pH and EC changes in the water layer of BMES. As the reaction proceeded, the DO in the BMES, S Control and W Control maintained a relatively stable range (4.3 ± 0.082 mg L −1 , 2.6 ± 0.21 mg L −1 and 6.5 ± 0.1 mg L −1 ). Higher DO in W Control should be related to the flushing water and keep the DO concentration similar to the DO in normal water (Fig. 4). The DO in BMES was 2 times higher than that in S control. In this study, changes in pH were observed for the entire duration. The experimental results clearly indicate that the pH of three reactors showed similar performance, and all were stable at 6.8-7.4. For the BMES, water conductivity was relatively constant at 135 μ s cm −1 , which was similar to that in the SMFC systems using fresh water as the electrolyte. This indicates a relatively low conductivity compared to the BMFC reported previously [22][23][24] .  The decomposition of organic would usually have robbed the dissolved oxygen from water which might deteriorate water quality. It is significant that the higher DO in BMES meant less DO was consumed in TOC or TN removal in BMES. When the EC of the three reactors were steady, W Control was 4-and 5-fold lower than the BMES and Control, respectively. It is significant that EC had the same trend with the TOC in the surface water. The difference on conductivity reduced with the decrease of TOC 25 .

Power generation in BMES.
Under the closed circuit (CC) mode, an average voltage of 90 mV was initially achieved in the BMES, although some voltage fluctuations were observed (Fig. 5A). The voltage decreased to 65 mV on day 11. In the present study, the internal resistance evaluated by the polarization slope method was about 40 Ω. This showed that the scale of mass transport limitations influenced the BMES performance, predominantly at the anode electrode. The potential and power density, as a function of current density, were obtained as shown in Fig. 5B. Based on the power curve, the maximum power density for the BMES was 81 mW m −2 at day 20, which is relatively higher compared with that in previous studies 26 .
For the BMES, the anode potential varied over a narrow range of − 410 to − 440 mV, while the cathode potential covered a much wider range of 240 to 180 mV, possibly indicating that its current was limited by the cathode or external resistance (Fig. 5C). A power management system (PMS) was designed that enables BMES to drive 9 red LEDs in parallel 27 . Six independent capacitor-based circuits were used to harvest electrical energy from the BMES. Each circuit consisted of capacitors (3.3 F, Panasonic Corporation, Japan) and relays controlled by programmable microcontroller (XD-J16H, Xunda Corporation, China). The capacitors in each circuit were charged in parallel by the corresponding module (1 min) and then charged capacitors connected in series to discharge (1 min) (Movie S1). Our results demonstrate that BMES can be a viable alternative renewable power source (Fig. S2).
There are also strong correlations between the current output of a simple anode-resistor-cathode device and rates of anaerobic microbial activity (TOC content) in a diversity of anoxic sediments 28 . At the initial stage of deployment, the BMES could obtain high voltage outputs, but over a long operational period, nutrient depletion, mass transfer resistance, and anode/cathode fouling, could lower the system performance 29,30 . Several factors, such as spatial and temporal variability, natural fluctuations, temperature changes, pressure variations, water flow changes, salinity and conductivity changes, and dissolved oxygen affect the BMES functions 29 . The kinetics of electron transfer from microorganisms to the anode were mainly restricted by the anode potential. Overall, the BMES revealed a high cathode potential, the lowest anode potential, and the highest cell potential (189-591 mV) 7 .
Community analysis. Analysis of the anodic bacterial community using 16S rDNA gene pyrosequencing revealed the enrichment of genera with potential exoelectrogenic capability 31 . And this has significant implications in determining the microbial community structure in SMFCs in each aquatic environment and consequently the fate and removal kinetics of organic pollutants in contaminated sediments 32 . The Shannon diversity index provided the species richness (i.e., the number of species present) and distribution of each species (i.e., the evenness of the species) among all the species in the community 33 . BMES had the highest diversity (Shannon = 7.07) that was slightly larger than that of Initial (Shannon = 6.87) and S Control (Shannon = 6.50) among the 4 communities, while S Control had the lowest diversity (Shannon = 6.09).
Earlier studies have observed that Gammaproteobacteria and Bacteroidetes occurred at high numbers within libraries from electrode biofilms, in which the Gammaproteobacteria were believed to be responsible for electron transfer and consequently power generation 23,34 . Flavobacteria and Bacteroidia have earlier been identified as the dominant bacteria for electricity generation 23,36 (Fig. 6). The microbial community composition observed in this study shows that the microbial community in BMES was dominated by exoelectrogenic bacteria 37,38 . A number of less prevalent bacteria were also detected (Fig. 6) 39 . The variety of the bacterial colonies in sediments is attributed to the process of enrichment of electricity bacteria in anodes 36 .

BMES construction and operation. A 195 L of the Benthic Microbial Electrochemical Systems (BMES)
was constructed using a three-dimensional anode (Tri-TDA) and floating air cathode (FAC). The BMES consisted of a 1.20 × 0.25 m plastic channel with a depth of 0.65 m. Sediment from Ash River, Heilongjiang Province was taken and used in this research. A sediment layer of 15 cm was placed in the bottom of the BMES (Total volume of sediment was 45.0 L). Water depth was 45 cm, and the total volume of water layer was 135 L. The floating air-cathode (0.1 m 2 ) was designed using rolling pressure activated carbon as a catalyst layer and stainless steel mesh (SSM) in the middle for electron collection and transfer 40 . The three-dimensional anode (Fig. S1B) was constructed using carbon mesh heat-pressure adhered to 'honeycomb' structure supports. The TDA (with height of 10 cm) was stood in the sediment of BMES, and full covered with river sediments. The distance between the TDA and FAC was connected by titanium wire with a 10 Ω external resistor. No other inoculum was used in the system and the entire experiment only relied on naturally occurring bacteria in the sediments. Two controls were set in this research with the same size and structure of the BMES as above mentioned. S Control was designed without electrodes fabrication and operated under the same conditions as BMES to compare pollutant removal efficiencies. W Control was designed not only with no electrodes fabrication but also washed by fresh water in a continuous flow rate of 400 L h −1 to investigate the organic removal efficiencies in contrast to the BMES. Experiments and sampling. Experiments were performed in a batch mode over a period of 60 d. Water samples were taken daily for TOC, TN, pH, EC (Electrical conductivity) and DO (Dissolved oxygen) monitoring. Sediment samples were taken every 15 days for TOC, TN detection, and PAHs in sediments were quantified after 60 days' operation.
Sediment samples were collected from the BMES, S Control and W Control using gravity core sampler. Two sampling spots were set in the middle along the width direction and 40 cm and 80 cm from the influent point along the long direction (Fig. S1C). Water Sampling points for water body analysis were placed at the side of the reactor, and two points were set in the vertical direction away 10 cm and 30 cm below the water surface. Samples were taken out from the reactors for further analysis (Fig. S1A).
Analysis and calculation methods. The voltages were recorded every 30 minutes using a data acquisition board (PISO-813.ICP DAS CO, Ltd) connected to a personal computer and converted to power according to Ohm's law, P = IU, where, P = power (W), I = current (A), and U = voltage (V). Anode and cathode potentials were measured against an Ag/AgCl reference electrode (Model 13-620-45; Fisher Scientific; Hampton, NH) suspended in the overlying water. For electrochemical processes such as MES, the most relevant production rate is thatnormalized byelectrode surface area, for example gm −2 d −1 , due to the limitation of electrode surface for directing these processes 41 . Power density (mW/m 2 ) and current density (mA/m 2 ) were calculated with the function of cathodic surface area (m 2 ). In order to estimate the energy produce by MFC during the experiment, ∆ V acquired by data logger was used to calculate the current (I) from Ohms law: I = V·R −1 . Tus, energy produced by BMES during 60 days operation was calculated as kilowatt hour per cubic meter of influent wastewater (kWh m -3 ) by the following equation: Analysis of TOC, TN and nitrate nitrogen was carried out according to APHA standard methods (APHA, 1998). Electrical conductivity (EC), pH and Dissolved oxygen (DO) were monitored using hand-held meters during the research period. Inductively coupled plasma-mass spectrometry (ICP-MS, Agilent 7500, Agilent, USA) was applied for determination of the heavy metals in this work. In order to know the PAH (Polycyclic aromatic hydrocarbons) contents and removal efficiencies in the sediments, PAHs were analyzed using external standards method. A standards contained twelve 2− to 6− ring PAHs were used in this study (Table S4). Accelerated Solvent Extractor system (ASE 350 System, Thermo Scientific ™ Dionex ™ , Waltham, MA) was used to extract the PAHs from sediment samples.
Microbial Community Analysis. The total community DNA in the sediment was directly extracted from the 0.25 g homogenized samples using the Omega D5625-01 Soil DNA Kit (Omega Bio-Tek, America) according to the manufacturer's instructions. Bacterial 16S rDNA was amplified from total DNA using 63 F forward (50-CAGGCCTAACACATGCAAGTC-30) and 1387 R reverse (30-GGGCGGWGTGTACAAGGC-50) primers. In the present metagenomic analysis, a single lane of the Illumina Hiseq 2000 sequencer was used for the biofilm and soil samples. Genes with BLAST hits to the NCBI-NR database (BLAST-hit genes) were subjected to MEGAN analysis to predict their taxonomic distribution. Initial soil samples from Ash River was also used as control to compare the community changes in the system.