Influence of Scale on Biomass Growth and Nutrient ... - ACS Publications

Oct 20, 2017 - However, many of these assessment tools utilize input .... were uploaded to the MG-RAST server, and paired-end reads ..... total nitrog...
1 downloads 12 Views 2MB Size
Subscriber access provided by Gothenburg University Library

Article

Influence of scale on biomass growth and nutrient removal in an algal-bacterial leachate treatment system Kaitlyn D Sniffen, Jacob Ryan Price, Christopher M Sales, and Mira Olson Environ. Sci. Technol., Just Accepted Manuscript • DOI: 10.1021/acs.est.7b03975 • Publication Date (Web): 20 Oct 2017 Downloaded from http://pubs.acs.org on October 21, 2017

Just Accepted “Just Accepted” manuscripts have been peer-reviewed and accepted for publication. They are posted online prior to technical editing, formatting for publication and author proofing. The American Chemical Society provides “Just Accepted” as a free service to the research community to expedite the dissemination of scientific material as soon as possible after acceptance. “Just Accepted” manuscripts appear in full in PDF format accompanied by an HTML abstract. “Just Accepted” manuscripts have been fully peer reviewed, but should not be considered the official version of record. They are accessible to all readers and citable by the Digital Object Identifier (DOI®). “Just Accepted” is an optional service offered to authors. Therefore, the “Just Accepted” Web site may not include all articles that will be published in the journal. After a manuscript is technically edited and formatted, it will be removed from the “Just Accepted” Web site and published as an ASAP article. Note that technical editing may introduce minor changes to the manuscript text and/or graphics which could affect content, and all legal disclaimers and ethical guidelines that apply to the journal pertain. ACS cannot be held responsible for errors or consequences arising from the use of information contained in these “Just Accepted” manuscripts.

Environmental Science & Technology is published by the American Chemical Society. 1155 Sixteenth Street N.W., Washington, DC 20036 Published by American Chemical Society. Copyright © American Chemical Society. However, no copyright claim is made to original U.S. Government works, or works produced by employees of any Commonwealth realm Crown government in the course of their duties.

Page 1 of 31

Environmental Science & Technology

1

Influence of scale on biomass growth and nutrient removal in an algal-bacterial leachate

2

treatment system

3 4

Keywords: algae, wastewater, leachate, scale-up, biomass growth, nutrient removal

5 6

Authors: Kaitlyn D. Sniffen, Jacob R. Price, Christopher M. Sales, Mira S. Olson*

7

Affiliation of all: Drexel University, 3141 Chestnut St., Philadelphia, PA, 19104

8

Corresponding author*: [email protected]

9 10

Abstract: Data collected from experiments conducted at a flask scale are regularly used as input

11

data for life cycle assessments and techno-economic analyses for predicting the potential

12

productivities of large scale commercial facilities. This study measures and compares nitrogen

13

removal and biomass growth rates in treatment systems that utilize an algae-bacteria consortium

14

to remediate landfill leachate at three scales: small (0.25L), mid (100L) and large (1000L). The

15

mid- and large-scale vessels were run for 52 consecutive weeks as semi-batch reactors under

16

variable environmental conditions. The small-scale experiments were conducted in flasks as

17

batch experiments under controlled environmental conditions. Kolomogov-Smirnov statistical

18

tests, which compare the distributions of entire data sets, were used to determine if the ammonia

19

removal, total nitrogen removal, and biomass growth rates at each scale were statistically

20

different. Results from the Kolmogov-Smirnov comparison indicate that there is a significant

21

difference between all rates determined in the large-scale vessels compared to those in the small-

22

scale vessels. These results suggest that small-scale experiments may not be appropriate as input

23

data in predictive analyses of full scale algal processes. The accumulation of nitrite and nitrate

ACS Paragon Plus Environment

Environmental Science & Technology

24

within the reactor, observed midway through the experimental process, is attributed to high

25

relative abundances of ammonia- and nitrite-oxidizing bacteria, identified via metagenomic

26

analysis.

27 28 29

Introduction: Algae are frequently investigated for use in nutrient removal from wastewater sources, as

30

the presence of nutrients can cause costly economic and ecological damage to receiving water

31

bodies [1, 2]. While it is most commonly suggested that algae-based nutrient removal be utilized

32

in domestic wastewater treatment, the direct treatment of landfill leachate may also be favorable.

33

Landfill leachate is a high strength liquid waste, containing concentrations of ammonia regularly

34

exceeding 800mgN/L, and is most commonly trucked to and treated at domestic wastewater

35

treatment facilities, costing landfills hundreds of thousands to millions of dollars each year.

36

Treating the leachate directly may reduce these costs to landfills. Few studies have examined the

37

treatment of landfill leachate by algae, and even fewer have done so on a long-term and large

38

scale [3-9].

39

Data from long-term, large scale studies should be used as input for predictive analyses

40

such as life cycle assessments (LCAs) and techno-economic analyses (TEAs) which are used to

41

evaluate if a commercial investment in a technology are environmentally favorable and

42

economically feasible, respectively. However, many of these assessment tools utilize input data

43

generated from small-scale studies, carried out under conditions that are significantly different

44

than those that occur on a commercial scale [10-14]. Studies conducted by Moody et al. (2014)

45

and Quinn et al. (2015) have shown that the wildly variable results from LCAs and TEAs are due

46

to inconsistent input data and system boundaries, which define the processes included in the

ACS Paragon Plus Environment

Page 2 of 31

Page 3 of 31

Environmental Science & Technology

47

assessments. These studies emphasize that it is crucial to use input data from studies that mimic

48

the environmental conditions and type of growth vessel of the proposed commercial scale

49

process in these predictive analyses [15, 16].

50

It has repeatedly been shown in studies using algae [15, 17], E.coli [18-21], and yeast

51

[20, 22] that production efficiencies from small-scale studies using microorganisms do not scale

52

proportionally to those occurring on a commercial-scale. For example, changes in vessel

53

characteristics, specifically surface area to volume ratios, [20, 23, 24], mixing dynamics, such as

54

types of mixing used, [20, 24, 25], and environmental conditions (e.g., temperature, light

55

exposure) [26-28], can significantly affect the activity of microorganisms [20, 23, 24]. Despite

56

this discrepancy between results from small-and large-scale studies, the vast majority of LCAs

57

and TEAs for evaluating the overall environmental impact and technical feasibility of large-scale

58

algae production systems use algal growth rates and densities from small-scale studies as input

59

values for their estimation of algal growth technologies [15, 16, 29].

60

The study presented here uses a system that couples algae biomass production with the

61

treatment of landfill leachate. Specifically, this study compares algal growth and nitrogen

62

removal using untreated leachate as a nutrient source in liquid cultures of a mixed algae/bacteria

63

consortium at three different scaled systems: 0.25L, 100L, and 1000L. Algal growth and

64

nitrogen removal rates measured in small-scale (0.25L) systems are statistically compared to

65

rates to observed at mid- (100L) and large-scales (1000L) to determine if results from small-

66

scale experiments are representative of larger-scale systems.

67

The small-scale (0.25L) vessels were operated as 7-day batch experiments [30-32] in an

68

environmental chamber under controlled laboratory conditions. The flask-scale and controlled

69

growth conditions were deliberately chosen as the comparison scale for this study as this

ACS Paragon Plus Environment

Environmental Science & Technology

70

experimental set up is similar to those most frequently used to generate input data for LCAs and

71

TEAs. In contrast, the mid- (100L) and large- (1000L) scale vessels were operated on a 7-day

72

semi-batch cycle but under the uncontrolled environmental conditions of a greenhouse. These

73

operating conditions were chosen as the environmental conditions of a full-scale outdoor system

74

will not be controlled.

75 76

Materials and Methods:

77

Nutrient feed and inoculum:

78

The nutrient feed used in this study was raw, untreated landfill leachate collected from the

79

Sandtown Landfill in Felton, DE. Leachate was refrigerated at 2ºC from collection until use.

80

Over the duration of the study, the ammonia content and pH of the leachate ranged from 320 to

81

935 mgN/L and 6.48 to 7.56, respectively. A more detailed analysis of the leachate composition

82

at each collection time point, data provided by the Sandtown Landfill, is described in Table S1 in

83

the Supplemental Information. The algae-bacteria culture used in this study was first collected

84

from a fish pond on the University of Pennsylvania campus. The inoculant was grown in the

85

urban greenhouse using leachate as a nutrient feed.

86 87

Sampling, preparation, and analysis of metagenome – The microbial population’s composition

88

and structure determines a bioreactor’s kinetic behavior, which is commonly observed through

89

chemical analysis. To assist in understanding these behaviors, metagenomic sequencing was

90

applied to a sample collected from one of the growth vessels. Total DNA was extracted from

91

four biological replicates collected prior to the beginning of the experimental period of this

92

study. DNA extraction was carried out immediately after sample collection following the

ACS Paragon Plus Environment

Page 4 of 31

Page 5 of 31

Environmental Science & Technology

93

methods described in detail within Price et al 2016 [33]. The four replicate DNA extractions

94

were then assessed for quality and concentration via NanoDrop 2000 and QuBit 2.0 and then

95

pooled and submitted for 2x250 paired-end sequencing on an Illumina HiSeq sequencer. The

96

Metagenomics RAST pipeline [34] was used to analyze the taxonomic and functional

97

composition of the resulting sequences. Raw reads were uploaded to the MG-RAST server and

98

paired-end reads were joined using the join-paired-ends function. Quality control measures were

99

applied including the removal of artificially induced sequencing artifacts [35], the removal of

100

sequences derived from H. sapiens [36], and the removal of low quality sequences (minimum

101

Phred score of 15, and all sequences were trimmed such that they contained a maximum of 5 low

102

quality base calls) [37]. Metagenome identification and quality control statistics are provided in

103

Table S2 in the Supplemental Information.

104 105

Mid & large scale experimental design:

106

Setup – This set of experiments used two 100L Plexiglas aquarium tanks (ATs) and two 1000L

107

raceway ponds (RWPs) from MicroBio Engineering (San Luis Obispo, CA). The 100L ATs and

108

the 1000L RWPs are referred to herein as mid- and large-scale vessels, respectively. The

109

working volumes of the ATs and RWPs were 60L and 600L, respectively. Examples of the mid-

110

and large-scale vessels used in this study can be found in Figure 1 A-B [38]. All vessels were

111

inoculated with an algae-bacteria culture, at the beginning of this study in February 2016, and

112

housed in a greenhouse on the roof of an academic building at Drexel University in Philadelphia,

113

PA [5, 38]. These vessels were exposed to semi-ambient conditions; no additional lighting or

114

heating/cooling was provided. Daylight hours ranged from 9hr:20min to 15hr:1min throughout

115

the year, though the actual amount of light to reach the cultures in the greenhouse was

ACS Paragon Plus Environment

Environmental Science & Technology

116

substantially less due to weather and urban obstructions. Temperatures of cultures ranged from

117

7.8 to 41.7°C.

118 119

Operation – The mid and large scale vessels were operated as semi-batch reactors, on a seven –

120

day batch cycle. Initial samples were taken at the beginning, and final samples were taken at the

121

end of each weekly cycle, where ammonia-N, nitrate-N, nitrite-N, and biomass density were

122

measured and analyzed. To prepare the system for the next weekly cycle, the entire volume of

123

the mid-scale AT was pumped into the large-scale RWP. After mixing, one third of the liquid

124

and biomass volume was removed and then replaced with water and leachate. A portion of this

125

mixture was then transferred back into the AT and the next week’s cycle would begin. The

126

volumes of the mid- and large-scale vessels were mixed together in between each week in order

127

to maintain equivalent initial nutrient and biological conditions at the two scales. The volume of

128

leachate added each week was increased slowly throughout the study, ranging from 5 to 40L per

129

week. Exact weekly leachate additions can be found in Figure S1 of the Supplemental

130

Information. These vessels were operated for 52 consecutive weeks, February 2016 – February

131

2017. A more detailed description of this method can be found in Sniffen et al (2017) [38].

132 133

Environmental Monitoring – Water temperature, pH, and dissolved oxygen within the tanks and

134

raceway ponds were collected at five-minute intervals throughout the yearlong study using

135

Neptune System’s APEX controller and probes (Morgan Hill, CA). Weekly maxima, minima,

136

and averages were calculated for each of these parameters.

137 138

Oxygen production analysis – The total suspended solids contained in the weekly samples was

ACS Paragon Plus Environment

Page 6 of 31

Page 7 of 31

Environmental Science & Technology

139

made up of live and dead cellular biomass as well as particulates from the leachate additions.

140

The rate of oxygen production was measured to determine the oxygen production activity

141

relative to the total suspended solids density [39, 40]. Oxygen production rates of the biomass

142

samples were measured at the beginning and end of each week. Rates from each sample were

143

measured using 450mm Unisense O2 probes (Aarhus N, Denmark) over 10 minute intervals,

144

which were run in biological triplicates. Each sample was run at 100% and 25% dilution of the

145

initial samples. The 25% dilution samples were found to give a better representative oxygen

146

production rate than samples with high biomass densities which may have experienced

147

significant self-shading (data not shown).

148 149

Mixing dynamics – The flow rate through the raceway pond was determined using the EPA float

150

method [41]. The mixing efficiency of the ponds was determined using a salt tracer test. Five

151

liters of a 250 g/L sodium chloride salt tracer solution was poured into the raceway pond. Water

152

samples were taken at 6 locations around the pond at 10-second intervals. The conductivity of

153

these samples was measured and the concentration of salt tracer was determined using a standard

154

curve.

155 156

Small-scale experimental design:

157

Flask-scale studies were performed using the same algae culture and leachate as those used in the

158

mid- and large- scale vessels. The ranges of nutrient and biomass concentrations used in these

159

small-scale studies mimicked those used in the mid- and large-scale studies. These small-scale

160

studies used 250 mL flasks with 100mL working volume, shaken continuously in an Innova 44

161

environmental chamber (Eppendorf, Hamburg, Germany), under constant 25°C temperatures and

ACS Paragon Plus Environment

Environmental Science & Technology

162

continuous light, as depicted in Figure 1 C. Each nutrient and biomass concentration condition

163

was tested using biological triplicates. The duration of all small-scale experiments was seven

164

days. Ammonia-N, nitrate-N, nitrite-N, and biomass density were measured and analyzed at the

165

beginning and end of each of the experiments.

166 167

Analytical methods:

168

Samples collected from experiments run at all scales were analyzed using the following

169

analytical methods: Ammonia-N, nitrate-N, and nitrite-N were measured using Hach test

170

methods 10031, 10020, and 8153, respectively, with a DR2400 spectrophotometer. Removal

171

rates of each nitrogen species along with total dissolved nitrogen were calculated using

172

Equations 1-2. In Eqn. 1-2, CoN and CfN are the initial and final concentrations of the nitrogen

173

species [mgN/L], respectively, where t is time [days] [5], NT is the total inorganic nitrogen

174

concentration [mgN/L], and Ri is the removal rate of species i [mgN/L/day] for each batch

175

period. These equations were developed based on preliminary studies which monitored nitrogen

176

removal and biomass growth every 1 to 2 days during week-long studies. These preliminary

177

studies found that the removal and growth rates were linear throughout the week.

178 179

(Eqn. 1)

NT= CNH3 + C NO2- + C NO3-

(Eqn. 2)

Ri= (CoN – CfN)/t

180 181 182 183

Biomass density was measured at the beginning and end of each week by standard total

184

suspended solid protocol using a 0.45µm filter and 20mL of sample [42]. Biomass density

ACS Paragon Plus Environment

Page 8 of 31

Page 9 of 31

Environmental Science & Technology

185

samples were run using biological duplicates. Weekly biomass growth rates were calculated

186

using equation 3, where Co,Biomass and Cf, Biomass are the initial and final biomass concentrations

187

[g/L], RB is the biomass growth rate [g/L/day], and t is time [days].

188 189

(Eqn. 3)

RB= (Cf, Biomass – Co, Biomass)/t

190 191

Statistical analysis – Weekly rates (nutrient removal and biomass growth) and environmental

192

data were pooled for regression and correlation analysis to identify general trends. Weekly

193

nutrient removal and biomass growth rates were then binned by scale for rate comparisons.

194

Kolmogov-Smirnov statistical tests are used to determine if the distributions of two data sets are

195

significantly different. Data binned by scale was compared using Kolmogov-Smirnov tests to

196

determine if the weekly rates calculated at each scale were statistically different. Regressions and

197

correlation data categories included: initial and final weekly ammonia-N, nitrate-N and nitrite-N

198

concentrations and removal rates; initial and final weekly biomass concentrations and growth

199

rate; initial and final weekly oxygen production rates; as well as maximum, minimum, and

200

average pH, water temperature, dissolved oxygen, and number of daylight hours. SPSS 24[43]

201

and Microsoft Excel were used to analyze this data.

202 203

Results and Discussion:

204

Environmental conditions – Seasonal trends of the temperature, pH, and dissolved oxygen (DO)

205

in the tanks and raceway ponds, over the course of the yearlong study can be seen in Figure 2A-

206

C. Over the course of a day, the pH, temperature, and dissolved oxygen in the tanks and raceway

207

ponds rise and fall on a diurnal cycle. This daily trend is shown over a representative two-week

ACS Paragon Plus Environment

Environmental Science & Technology

208

span, depicted in Figure 3. As confirmed by observed data, the ATs and RWPs were subjected

209

to identical environmental conditions due to their colocation.

210

(Figure 2A-C)

211

(Figure 3)

212

Mid- and large-scale experimental results – The ranges in initial ammonia, total nitrogen and

213

biomass concentrations over the course of the study were 0.2-161 mgN/L, 8.12- 187.8 mgN/L,

214

and 40-1780 mgBiomass/L, respectively. A breakdown of initial conditions by scale is presented

215

in Table S3 in the Supplemental Information. Leachate input volumes were approximately 20

216

L/week during the first half of this study, while the input volumes during the second half were

217

approximately 30-40L/week. During times of high nitrite concentrations in the systems, no

218

leachate was added. Weekly leachate input volumes during this study can be found in Figure S1

219

of the Supplemental Material.

220 221

The dissolved nitrogen content of the leachate was 99% in the form of ammonia. Ammonia and

222

total nitrogen removal rates along with biomass growth rates varied over the course of the study;

223

the average ± standard deviation of these rates for the large-, mid- and small-scales are presented

224

in Table 1. The ammonia and total nitrogen removal rates for each of the vessels are presented

225

in Figure 4 A-B.

226 227

Table 1: Average ± standard deviation of ammonia removal, total nitrogen removal, and biomass growth rates, at each scale. Large-scale Mid-scale Small-scale Total N Removal Rate 1.63 ± 2.95 1.82 ± 3.97 3.16 ± 2.73 (mgN/L/day) Ammonia Removal Rate 2.33 ± 3.05 2.68 ± 3.90 2.96 ± 2.71 (mgN/L/day) Net Biomass Growth Rate 12 ± 24 4 ± 38 18 ± 21 (mg/L/day)

228

ACS Paragon Plus Environment

Page 10 of 31

Page 11 of 31

229

Environmental Science & Technology

(Figure 4 A-B)

230 231

The maximum biomass density measured over the course of this study was 2100 mg/L. Biomass

232

growth rates ranged from -106 to 199 mg/L/day. A negative growth rate means that there was

233

biomass loss over the course of the week. Biomass loss was seen in 70 out 205 weeks’ worth of

234

data from all vessels. Negative biomass growth rates did not statistically favor any scale, season,

235

or environmental condition. Detailed discussion of this point can be found in Sniffen et al. (2017,

236

submitted). When accounting for only weeks with positive growth rates the average and standard

237

deviation was 28 ± 32 mg/L/day. This range of growth rate is similar to rates reported for algae

238

growth fed with domestic wastewater [44-47]. Weekly biomass growth rates for all vessels can

239

be found in Figure S2 of the Supplemental Information.

240 241

Metagenomics results – Alpha diversity within the metagenomic sample was estimated to be 376

242

taxa. Annotations from the NCBI Reference Sequence Database [48-50] and the M5nr protein

243

sequence database [51], acquired via the MG-RAST interface, were used to analyze and interpret

244

the taxonomic composition of the microbial community. The vast majority of reads were

245

attributed to the domains of Bacteria and Eukaryota (Table S4 in the Supplemental Information).

246

Proteobacteria and Bacteroidetes were the first and second most abundant phyla respectively

247

(Table S5 in Supplemental Information). Under the annotations for both reference databases,

248

Chlorophyta and Cyanobacteria comprised about 2% of the annotated reads (Table S5).

249

Chlamydomonas reinhardtii, Volvox carteri, and Scenedesmus obliquus accounted for the bulk

250

of the reads attributed to Chlorophyta (Table S6 in the Supplemental Information). The

251

distribution of reads amongst species falling within the phyla Cyanobacteria was much more

ACS Paragon Plus Environment

Environmental Science & Technology

252

uniform at the species level (Table S7 in the Supplemental Information); agglomerating taxa by

253

Genus indicates that the genus Synechococcus accounts for roughly 30% of the reads within

254

Cyanobacteria, followed by Cyanothece spp. at 19%, and Nostoc spp. at 12% within this phyla.

255

AOB and NOB were observed to comprise roughly 12% and 1% of the total microbial

256

population (Table S8 in Supplemental Information). The anammox bacteria Candidatus Kueneia

257

was also discovered in the M5nr annotations. This indicates that anammox may be a potential

258

nitrogen transformation pathway within this system (Table S8).

259

Low abundances for eukaryotes, such as Chlorophyta, within this study, are common in

260

metagenomic analyses, arising from their being underrepresented within reference databases,

261

including those used by MG-RAST [52]. This induces bias towards the identification and

262

annotation of bacterial sequences, while sequences from eukaryotic sources may fail to be

263

annotated. Because of this, the abundances of Chlorophyta reported here should be viewed as a

264

floor or minimum and should not be directly compared to the relative abundance of the kingdom

265

Bacteria, or the taxa therein.

266 267

Statistical analysis – A Pearson’s correlation matrix was calculated for all of the categories listed

268

in the statistical methods section. Strong, statistically significant correlations were found

269

between: initial ammonia concentration and ammonia removal rate (r=0.854, p