Maximizing Productivity and Reducing ... - ACS Publications

Feb 29, 2016 - ABSTRACT: The ability to dynamically control algal raceway ponds to maximize biomass productivity and reduce environ- mental impacts (e...
0 downloads 0 Views 2MB Size
Subscriber access provided by UNIV OSNABRUECK

Article

Maximizing productivity and reducing environmental impacts of full-scale algal production through optimization of open pond depth and hydraulic retention time Quentin Béchet, Andy Shilton, and Benoit Guieysse Environ. Sci. Technol., Just Accepted Manuscript • DOI: 10.1021/acs.est.5b05412 • Publication Date (Web): 29 Feb 2016 Downloaded from http://pubs.acs.org on March 4, 2016

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 32

Environmental Science & Technology

1

Maximizing productivity and reducing environmental

2

impacts of full-scale algal production through optimization

3

of open pond depth and hydraulic retention time

4 5

Quentin Bécheta,b,*, Andy Shiltona, Benoit Guieyssea

6

a

7

Palmerston North 4442, New Zealand

8

b

INRIA BIOCORE, BP 93 06902 Sophia Antipolis Cedex, France

9

*

Corresponding author

School of Engineering and Advanced Technology, Massey University, Private Bag 11 222,

10

Telephone: +33 4 92 38 71 74

11

Fax: + 33 4 92 38 78 58

12

Email: [email protected]

ACS Paragon Plus Environment

1

Environmental Science & Technology

Page 2 of 32

13

Abstract

14

The ability to dynamically control algal raceway ponds to maximize biomass productivity and

15

reduce environmental impacts (e.g. land and water use) with consideration of local constraints

16

(e.g., water availability, climatic conditions) is an important consideration in algae

17

biotechnology. This paper presents a novel optimization strategy that seeks to maximize

18

growth (i.e. optimize land use), minimize respiration losses, and minimize water demand

19

through regular adjustment of pond depth and hydraulic retention time (HRT) in response to

20

seasonal changes. To evaluate the efficiency of this strategy, algae productivity and water

21

demand were simulated at 5 different climatic regions. In comparison to the standard

22

approach (constant and location-independent depth and HRT), dynamic control of depth and

23

HRT was shown to increase productivity by 0.6 - 9.9% while decreasing water demand by 10

24

- 61% depending upon the location considered (corresponding to a decrease of the water

25

footprint by 19 - 62%). Interestingly, when adding as a constraint that the water demand was

26

limited to twice the local annual rainfall, higher net productivities were predicted under

27

temperate and tropical climates (15.7 and 16.7 g-m-2-d-1, respectively) than from

28

Mediterranean and subtropical climates (13.0 and 9.7 g-m-2-d-1, respectively); while algal

29

cultivation was not economically feasible in arid climates. Using dynamic control for a full-

30

scale operation by adjusting for local climatic conditions and water constraints can notably

31

affect algal productivity. It is clear that future assessments of algae cultivation feasibility

32

should implement locally optimized dynamic process control.

ACS Paragon Plus Environment

2

Page 3 of 32

33

Environmental Science & Technology

Introduction

34 35

In the wake of massive investments to increase algae biofuel production, numerous studies

36

assessed the economic feasibility and environmental impacts of microalgae cultivation in

37

raceway ponds (Table 1). In contrast to this study, the majority of these assessments (Table 1)

38

set a constant pond depth (typically 0.2 - 0.3 m) and did not consider the hydraulic retention

39

time (HRT) despite the well-acknowledged impact of these parameters on algal

40

productivity1,2. Likewise, geographic locations have been compared on the basis of how

41

climatic conditions impact productivity using identical process parameters regardless of their

42

practical suitability for a particular climate3,4,5,6. It follows that biomass productivities in these

43

assessments may have underestimated the maximum productivity practicably achievable.

44

More importantly, underestimating productivity may negatively impact the conclusions of

45

economic and impact assessments. For example, the water footprint and land use efficiency of

46

full-scale algal cultivation may have been overestimated by previous studies. Process

47

optimization is therefore critical to any algal biotechnology applications from high-value

48

metabolite production to algae biofuel and algal wastewater treatment7.

49 50

Different optimization strategies have been investigated to maximize system productivity by

51

ensuring that growth conditions are near optimal at all times during cultivation. While

52

controlling pH and nutrient concentrations may be feasible at full-scale7, temperature control

53

is economically prohibitive for biomass production (e.g., biofuel feedstock). As a result,

54

temperature fluctuations should be considered in any attempt to optimize process operation8.

55

A second strategy is to optimize the geometry of the cultivation system by for example

56

adjusting the position of photobioreactors according to latitude and time of the year to

57

minimize mutual shading between reactors9. Ritchie and Larkum10 adjusted pond depth to

ACS Paragon Plus Environment

3

Environmental Science & Technology

Page 4 of 32

58

optimize water temperatures. While useful, these optimization strategies were based on

59

annually averaged meteorological data. Further optimization can be achieved by dynamically

60

adjusting pond depth and HRT over time as environmental conditions change. A third strategy

61

is to adjust algal concentrations to optimize light availability while minimizing biomass

62

respiration losses11. In practice, this requires constant adjustment of the harvest rate as

63

proposed by Muñoz-Tamayo et al.12 because algal growth and respiration vary with

64

meteorological conditions. Altogether, these past studies demonstrate that process

65

optimization yields quantitatively significant productivity gains but also requires models that

66

accurately predict productivity as a function of design, operational, and environmental

67

parameters. Because these previous optimization studies did not account for the impact of

68

temperature variation during full-scale cultivation and/or were not validated against full-scale

69

data, additional development is required. Furthermore, it is also critical to consider

70

environmental impacts during optimization because full-scale feasibility can be severely

71

limited by freshwater availability and other environmental concerns13.

72 73

In view on the above-mentioned uncertainties, this paper presents a novel optimization

74

strategy based on the regular adjustment of pond operation in response to seasonal

75

meteorological fluctuations to maximize algal productivity and minimize water demand. This

76

study focuses in particular on the dynamic control of pond depth and HRT that are easy to

77

regulate in practice. We examine how dynamic control yields materially different predictions

78

of productivity and water demand. For this purpose, the algal productivity and process-water

79

demand were simulated under various control strategies and 5 different climates using a

80

predictive model of algal productivity accounting for temperature and validated against full-

81

scale data14,8.

ACS Paragon Plus Environment

4

Page 5 of 32

Environmental Science & Technology

82

1

Materials and methods

83

1.1

84

This study focuses on algal cultivation in open raceway ponds mixed with paddle wheels15 as

85

these systems are the most economical for large-scale, bulk biomass cultivation16,17. Chlorella

86

vulgaris was used as a commercially relevant species18,19 given the availability of a predictive

87

productivity model validated at full-scale for this species (in particular, this model can be

88

adapted to any geometry and requires the fewest biological inputs, see S1 for further

89

discussion). The ponds were operated in a semi-continuous regime from the first day of

90

operation by withdrawing a volume of culture at 7 pm every day (to minimize respiration

91

losses at nighttime) and immediately replacing this volume with fresh medium. Fresh water

92

was added before harvesting to compensate for daily evaporation losses. The hydraulic

93

retention time (HRTi, in days) is:

94

HRTi =

95

where Vi, S, and di are the culture volume (m3), surface area (m2), and depth (m), respectively,

96

and Vs,i is the volume of culture extracted daily (m3-d-1); the subscript i indicates the day.

Cultivation system

Vi S ⋅ di (1) = Vs,i Vs,i

97 98

1.2

Water demand and water footprint

99

The annual water demand (WD, m3-yr-1) was calculated as:

100

WD = Qev + 0.25⋅ Qout (2)

101

where Qev (m3-yr-1) represents the rate of water evaporation and Qout (m3-yr-1) the rate of

102

culture withdrawal from the pond. As a portion of the liquid medium withdrawn is recycled

103

back into the pond following cell harvesting, and in the absence of consensus on the medium

104

recycling ratio that can be achieved during full-scale cultivation, we assumed 75% recycling

105

based on the literature13. The rate of water lost by leaks should be relatively minor for a pond

ACS Paragon Plus Environment

5

Environmental Science & Technology

Page 6 of 32

106

with a well-constructed liner and was therefore neglected to focus on the impact of climate

107

and process operation on water demand. The rate of evaporation Qev was calculated using the

108

model of Béchet et al.14. By assuming that open ponds are operated 365 days of the year, Qout

109

can be calculated as the sum of the volumes Vs,i (m3-d-1) extracted from the pond each day

110

(see Equation 1). The impact of precipitation was accounted for by assuming that the

111

freshwater inflow to the pond would be reduced in proportion to the rain landing on the pond.

112

The HRT was thus maintained and no freshwater resource was unnecessarily wasted. When

113

precipitations exceeded freshwater needs, the algal broth was slightly diluted and this had no

114

impact on predictions.

115 116

The water footprint (in m3-kg-1) was calculated as the ratio of the annual water demand over

117

the annual biomass productivity (section 1.3). This corresponds to the 'blue-green' water

118

footprint as defined by Gerbens-Leenes et al.20 ('grey' water footprint was excluded as it

119

requires locally dependent data).

120 121

1.3

Biomass productivity

122

Biomass productivities were predicted based on the model of Béchet et al.21, which essentially

123

computes net productivity as the difference between "photosynthetic output" (i.e., the biomass

124

produced through photosynthesis) and "respiration loss" (i.e., the biomass consumed by the

125

algae for non-growth related processes; see S1 for details). In this paper, the negative impact

126

of respiration losses on productivity is expressed as the percentage of respiration loss divided

127

by the photosynthetic output over a given period of time.

128 129

Yearly net productivities were computed at different locations by predicting (see S1 for

130

details): (1) hourly temperatures using the model developed by Béchet et al.14 and validated

ACS Paragon Plus Environment

6

Page 7 of 32

Environmental Science & Technology

131

against one year of data collected in an open pond in New Zealand; (2) hourly light intensities

132

within the culture using a model based on the Beer-Lambert law21–23; and (3) hourly

133

productivities using the model of Béchet et al.21 using the temperature and light distribution

134

profiles determined in steps 1 and 2. The pond surface area was set to 100 m2 because the

135

temperature model for ponds was validated on a system of similar size14. The estimated areal

136

productivities (in g-m-2-d-1) are considered representative of full-scale ponds. The initial algal

137

concentration was set to 0.3 g-L-1 in all simulations to ensure optimal light utilization (this

138

initial value had a negligible impact on overall productivity predictions, S2).

139 140

The productivity model of Béchet et al.21 was validated in indoor photobioreactors (accuracy

141

of +/- 15% over 163 days8) and in outdoor photobioreactors (an accuracy of +/- 8.4% over

142

148 days of cultivation21). These results indicate that the model can predict productivity for

143

different geometries, under variable weather conditions, and over the full ranges of light and

144

temperature conditions experienced in open ponds. For these reasons, the model should

145

theoretically predict productivity in open ponds with an accuracy on the order of +/- 10%.

146

Measured productivity of Chlorella sp. in an outdoor thin-layer photobioreactor located in

147

Czech Republic in summer time24,25,26 varied between 11.1-32.2 g m-2 d-1, which is in the

148

range of our predictions between 15.4 and 25.4 g m-2 d-1 in an open pond (assuming a depth of

149

0.1m and a HRT of 7d). Similarly, Moody et al.27 predicted productivities of 11.7 and 12.8

150

g.m-2-d-1 in open ponds in Arizona and Hawaii, respectively, against 12.5 and 11.2 g.m-2-d-1

151

with our model. Even if this comparison is based on systems with different geometries,

152

operational regimes, locations and most likely algal strains, the fact that very similar

153

productivities were predicted contributes to indirectly validate the productivity model in open

154

ponds (see S1 for further details).

ACS Paragon Plus Environment

7

Environmental Science & Technology

Page 8 of 32

155

The model predicts that photosynthesis of C. vulgaris over a short period of time (~1 hour) is

156

theoretically maximal at 38oC with a maximal temperature around 42oC. Experimental

157

evidence however shows temperatures over 35oC negatively affect productivity only when

158

algae are exposed to these temperatures over 1-2 days. While our biological model cannot

159

predict the impact of these ‘long duration high temperature’ events, we accounted for this

160

potential stability issues by introducing an "overheating criterion." This criterion was defined

161

as the duration of daytime exposure to temperatures higher than 35oC was therefore computed

162

along with the predicted productivities at the five climatic locations.

163 164

Simulations were performed at five locations representative of five climates: Yuma, Arizona,

165

USA (arid climate); Redding, California, USA (Mediterranean); Hilo, Hawaii, USA

166

(tropical); Sebring, Florida, USA (subtropical); and Hamilton, New Zealand (temperate; see

167

Béchet28 for more details on the climatic characteristics of each location). Hourly weather

168

data used in the productivity model were obtained from the US National Climate Data Centre

169

(NCDC) and from the New Zealand National Institute for Water and Atmospheric Research

170

(NIWA). For the US locations, the hourly relative humidity data were not available and

171

monthly averages were used instead.

172 173

1.4

Impact of depth and HRT on productivity and water demand

174

The examples shown in Boxes 1 and 2 illustrate how pond (or water-column) depth and pond

175

HRT impact the productivity and the water demand in raceway ponds. The concepts

176

introduced through these examples have already been presented in previous studies but their

177

comprehension is necessary to fully understand the novel optimization strategy presented in

178

this study.

179

ACS Paragon Plus Environment

8

Page 9 of 32

Environmental Science & Technology

180

Box 1: Impact of pond depth on productivity and water demand

181

For a given cultivation surface area, the thermal inertia is directly correlated to the mass of

182

water in the pond and, therefore, the water-column depth. Hence, the shallower the pond, the

183

greater the diurnal temperature fluctuations experienced by the algae. For a given HRT,

184

decreasing depth slightly improves overall productivity by enabling daytime temperatures to

185

climb close to the optimal temperature for photosynthesis and inducing a faster decrease of

186

temperature at nighttime, thereby reducing nighttime respiration. Productivity can therefore

187

be increased by 8% in a Mediterranean climate when the pond depth decreases from 0.5 m to

188

0.1 m (Figure 1). Decreasing pond depth, however, also increases overheating risks as shown

189

in Figure 1. The frequent occurrence of high temperatures may lead to culture collapse, which

190

would affect yearly productivity and inoculation costs. The reduction in water demand

191

evidenced in Figure 1 is largely driven by a reduction in the rate of water withdrawal from the

192

pond (see Equation 1) despite an increase in daytime evaporation losses as temperature

193

increases.

194 195

Box 2: Impact of HRT on productivity and water demand

196

Increasing HRT by reducing the rate of culture withdrawal from the pond enables algae cell

197

concentration to increase. In practice, the cell concentration at which all the available light is

198

captured is relatively low (e.g. 0.1 kg m-3) so relatively high HRT values negatively affect net

199

productivity by increasing respiration losses as algae stay in the pond longer. For example, the

200

optimal HRT was predicted to approximately 4 days in a Mediterranean climate (Figure 2), in

201

agreement with empirical observations reported in the literature11,29.

202

ACS Paragon Plus Environment

9

Environmental Science & Technology

203

Because a higher HRT entails a lower withdrawal rate from the pond, water demand

204

decreases when HRT increases as seen in Figure 2 (HRT has a negligible impact on the rate

205

of evaporation from the pond surface).

Page 10 of 32

206 207 208

1.5

Process optimization

209

Simulations were conducted for each location to (1) quantify the impact of varying depth and

210

HRT on productivity and water demand and to (2) compare "locally-optimized" operation

211

against "standard conditions" defined as a pond depth of 0.25 m and a HRT of 10 d based on

212

literature values (Table 1). For locally-optimized operation, pond depth and HRT were altered

213

on a monthly basis with the objectives of minimizing the water demand and maximizing

214

productivity by reducing overheating risks and respiration losses (the time step of one month

215

is arbitrary and may be further refined in future studies). This optimization was undertaken by

216

first optimizing depth and HRT independently as described in Table 2. The co-optimization of

217

depth and HRT was then performed by determining the optimal depth monthly profile before

218

optimizing the HRT monthly profile. A final case study investigated the co-optimization of

219

depth and HRT under the constraint of a limit on water demand (see section 2.5 for details).

220

In practice, the optimal HRT and depth monthly profiles were determined by using iterative

221

dichotomic searches.

222

223

2

Results and Discussion

224

2.1

225

The primary objective of this paper is to demonstrate a new operation optimization strategy

226

for full-scale algal cultivation in raceway ponds. Process optimization involves maximizing

Objective and criterion for optimization

ACS Paragon Plus Environment

10

Page 11 of 32

Environmental Science & Technology

227

productivity (and thus revenues) while minimizing resource use (and thus costs and

228

environmental impacts). Productivity can be severely impacted if the pond overheats or if

229

high-HRT operation entails excessive biomass respiration losses (Box 2). The following

230

simulations were therefore conducted under the criteria that "daytime pond temperature must

231

not exceed 35oC more than 20 hours per month for any month of the year" to minimize

232

overheating risks (criterion 1) and "respiration must consume less than 50% of the

233

photosynthetic output over the year" to maximize productivity (criterion 2). With regard to

234

resource use, water demand is particularly site-specific as it depends on availability and cost

235

of freshwater supply, the capacity of the receiving environment to re-assimilate process

236

wastewater, and the degree of associated local authority regulations. While these constrains

237

are variable (Guieysse et al. 2013), the impact of water use should, to a large degree, be

238

linked to the scarcity of local freshwater resources. Thus for the simulations presented in this

239

paper, we have adopted a third criterion that "yearly water demand must not exceed two times

240

the yearly rainfall at the location considered" (criterion 3).

241 242

While the first criterion is based on experimental observations21, acceptable levels of

243

respiration losses and water demand are subjective with other alternatives opened to being

244

assessed. The modeling platform used in this paper however enables examining the

245

implication of varying these criteria and this is discussed in section 2.5.

246 247

2.2

Productivity under constant depth and HRT operation

248

Productivity and water demand were simulated at five different climatic locations under

249

constant depth and HRT operation. Depending on the location, algal productivity under

250

standard conditions varied between 10.7 and 12.5 g-m-2-d-1. Interestingly, these productivities

251

are relatively low compared to values used in prior assessments found in the literature (Table

ACS Paragon Plus Environment

11

Environmental Science & Technology

Page 12 of 32

252

1). These low values are unlikely due to predictions inaccuracy as the model accuracy should

253

be around +/- 10% (see section 1.3). As discussed by Béchet et al.8, the overestimation of

254

productivity by prior studies is most likely due to the common assumption that temperature

255

has a negligible impact on productivity. Water demand was between 1.4 and 4.4 m3-m-2-yr-1

256

(Figure 3), which is within the range of values reported in the literature (Table 1).

257

Significantly lower values reported in previous studies are mostly explained by the fact that

258

water demand was computed as water evaporation and process water was not included, as

259

detailed by Guieysse et al.13.

260 261

For each depth and HRT combination, Figure 4 shows compliance against the three criteria

262

defined in section 2.1 (central part of each sub-figure). No combination of constant depth-

263

HRT operation were predicted to satisfy all three criteria for arid, Mediterranean, and

264

subtropical climates while several combinations ensuring low biomass losses, low

265

overheating risk, and low water demand were predicted for temperate and tropical climates

266

(Figure 4). These two climates therefore have potential for algal cultivation if pond depth and

267

HRT are maintained constant over the year.

268 269

2.3

Dynamic depth control

270

As explained in Box 1, decreasing pond depth efficiently minimizes water demand (related to

271

criterion 3), generates modest productivity gains (related to criterion 2), but also increases

272

overheating risks (related to criterion 1). In this simulation, pond depth was therefore

273

dynamically adjusted on a monthly basis to the lowest possible value that still ensured pond

274

temperature did not exceed 35oC more than 20 hours during the month considered (criteria 1).

275

ACS Paragon Plus Environment

12

Page 13 of 32

Environmental Science & Technology

276

As seen in Figure 4 (right column of each sub-figure), dynamic depth control was predicted to

277

remove overheating risks while still enabling operation at relatively small depths (< 0.25 m,

278

Figure 5) at all locations considered. Compared to standard operation (constant depth = 0.25

279

m, constant HRT = 10 d), dynamic depth control at a constant HRT of 10 d was predicted to

280

slightly increase productivity (4-11%) and significantly reduce water demand (15-53%) as

281

shown in Figure 3. Similarly, water footprint decreased from 19-58% depending on the

282

location. As explained in Box 1, the smaller thermal mass associated with a shallower pond

283

increased daytime temperatures (0.5-1.0oC on average over the year) and decreases nighttime

284

temperatures (1.4-2.1oC on average over the year). These changes in temperature profile are

285

predicted to cause an increase of daytime photosynthetic output and a decrease of nighttime

286

respiration loss, which explains the gain in net productivity evidenced in Figure 3. The

287

reduction in water demand was achieved by decreasing the culture volume and thus the

288

culture withdrawal rate required to maintain the HRT at 10-d.

289 290

As shown in Figure 4, dynamic depth control operation was simulated under a broad range of

291

constant HRTs. As can be seen from the "Depth control" column on the right hand side on

292

each sub-figure, dynamic depth control at fixed HRT was predicted to significantly improve

293

the feasibility of algal cultivation in hot climates (e.g., all criteria were satisfied using

294

dynamic depth control at 10-d HRT in the Mediterranean climate, Figure 4). Dynamic depth

295

control was also forecasted to expand the range of acceptable HRT values under tropical and

296

temperate climates (Figure 4).

297 298

2.4

Dynamic HRT control

299

As explained in Box 2, decreasing the HRT minimized respiration losses (criterion 2), but

300

also increased water demand (criterion 3). The optimization strategy investigated in this

ACS Paragon Plus Environment

13

Environmental Science & Technology

Page 14 of 32

301

section therefore involved dynamically adjusting the HRT on a monthly basis to the highest

302

possible value while still ensuring that respiration losses were lower than 50% of the monthly

303

photosynthetic output at the standard depth of 0.25 m (criterion 2).

304 305

At locations where the average HRT was higher than 10 d (Mediterranean and temperate

306

climates, see S5), dynamic HRT control was predicted to reduce the water demand by 6-13%

307

but slightly decreased productivity in comparison to standard 10 d HRT operation.

308

Conversely, productivity gains were obtained at the cost of a higher water demand at the other

309

climatic locations considered. Similarly, water footprint variation was either positive or

310

negative (-4.7-8.8%) depending on the average value of the HRT. These differences between

311

locations illustrate how the criterion for acceptable respiration losses (here 50%) impacts

312

water demand and productivity. For example, increasing this arbitrary threshold would lead to

313

lower productivity but also to reduce water demand. This threshold must therefore be adjusted

314

by considering the availability of freshwater resources as further discussed in section 2.5.

315 316

Dynamic HRT control was subsequently applied at different pond depths and the resulting

317

water demands and productivities were assessed. As illustrated by the upper "HRT control"

318

row of each sub-figure disclosed on Figure 4, dynamic HRT control was predicted to

319

efficiently reduce respiration losses, but the overheating risks and water demand criteria could

320

still not be met in the arid, Mediterranean, or subtropical climates.

321 322

2.5

Dynamic depth and HRT co-control

323

This section considers the possibility of dynamically adjusting both depth and HRT on a

324

regular basis. This co-control was forecasted to reduce water demand by 10-61% compared to

325

standard operation under the five climates considered (Figure 3). Productivity was also

ACS Paragon Plus Environment

14

Page 15 of 32

Environmental Science & Technology

326

increased by 0.6-9.9% depending on the location considered, which is mostly explained by

327

the monthly adjustment of the depth as described in section 2.3. The water footprint was in all

328

cases reduced by 19-62%.

329 330

Of particular interest, dynamic co-control of depth and HRT allowed low water demand under

331

Mediterranean, tropical, and temperate climates. As a result, respiration losses could be

332

further minimized at these locations while still respecting the rainfall criterion. In contrast, a

333

respiration loss higher than 50% must be accepted under arid and subtropical climates for the

334

rainfall criterion to be respected (top right corner of each subfigure of Figure 4).

335 336

In the final case study, depth and HRT co-control was used to maximize productivity while

337

ensuring that water demand met the rainfall criterion. In these simulations, the maximum

338

amount of water that can be used was fixed by the amount of available freshwater resources

339

(defined as twice the amount of annual rainfall) instead of the maximum respiration losses

340

(previously 50%). Under this constraint, the fraction of biomass lost by respiration was, as

341

expected, lower than 50% under Mediterranean, tropical, and temperate climates (Table 3). In

342

arid climates, no operation could meet the rainfall criterion because the rate of evaporation

343

alone exceeded twice the annual rainfall. Surprisingly, productivities predicted at temperate

344

and tropical locations exceed the productivities predicted at Mediterranean and subtropical

345

locations (Figure 3). This is explained by the fact that high rainfalls and low evaporation rates

346

at the tropical and temperate locations allow for low-HRT operation, which reduces

347

respiration and therefore boosts productivity. The potential impact of freshwater availability

348

was confirmed when the amount of available freshwater resources was set to 1-3 times the

349

amount of annual rainfall (see S6 for details). For example, these simulations showed that

350

allowing for a lower water use (by reducing the rainfall criterion from 2 to 1) would make

ACS Paragon Plus Environment

15

Environmental Science & Technology

351

impossible algal production in arid, Mediterranean, and subtropical climates as evaporation

352

would be higher than rainfall.

Page 16 of 32

353 354 355

2.6

Does optimization impact assessment of feasibility, environmental impacts, and economics?

356

Figure 4 illustrates how locally targeted process optimization using dynamic depth and HRT

357

co-control has the potential to reduce overheating risks and water demand. For instance, while

358

algal cultivation under standard conditions was shown to be unviable in a Mediterranean

359

climate (section 2.2), optimization through depth and HRT co-control is predicted to allow all

360

criteria to be satisfied (Figure 4) with increased productivity and reduced water demand

361

(Figure 3). In addition, accounting for environmental constraints such as water availability

362

can significantly change productivity estimates as detailed in section 2.5. While the

363

simulations presented here are dependent upon the criteria set, the results indicate that

364

attempting to assess feasibility without considering process optimization or local constraints

365

may lead to materially different outcomes. This raises the risk that prior assessments may

366

have inaccurately quantified full-scale practically achievable field outcomes.

367 368

List of supporting information

369

S1: Description of the productivity model; S2: Impact of initial concentration on productivity

370

predictions; S3: Impact of pond depth on productivity, water demand, and overheating risks at

371

five climatic locations; S4: Impact of pond HRT on productivity and water demand at five

372

climatic locations; S5: Optimal HRT profiles at constant depth at five climatic locations; S6:

373

Impact of the rainfall criterion on yearly productivity when depth and HRT are co-controlled

374

under a constraint on the water demand

375

ACS Paragon Plus Environment

16

Page 17 of 32

Environmental Science & Technology

376

Acknowledgements

377

Q.B. was supported by a Massey University Doctoral Scholarship. The world map used in the

378

Abstract and graphical (TOC) originated from the website openclipart.org. We also thank the

379

Reviewers for their contribution.

ACS Paragon Plus Environment

17

Environmental Science & Technology

380 381

Page 18 of 32

References (1) Slegers, P. M.; Lösing, M. B.; Wijffels, R. H.; van Straten, G.; van Boxtel, A. J. B.

382

Scenario evaluation of open pond microalgae production. Algal Res. 2013, 2 (4), 358–

383

368.

384

(2)

385 386

Lee, E.; Pruvost, J.; He, X.; Munipalli, R.; Pilon, L. Design tool and guidelines for outdoor photobioreactors. Chem. Eng. Sci. 2014, 106, 18–29.

(3)

Batan, L.; Quinn, J. C.; Bradley, T. H. Analysis of water footprint of a photobioreactor

387

microalgae biofuel production system from blue, green and lifecycle perspectives.

388

Algal Res. 2013, 2 (3), 196–203.

389

(4)

390 391

production scale-up. Appl. Energy 2011, 88 (10), 3377–3388. (5)

392 393

Pate, R.; Klise, G.; Wu, B. Resource demand implications for US algae biofuels

Wigmosta, M. S.; Coleman, A. M.; Skaggs, R. J.; Huesemann, M. H.; Lane, L. J. National microalgae biofuel production potential and resource demand. 2011, 47, 1–13.

(6)

Yang, J.; Xu, M.; Zhang, X.; Hu, Q.; Sommerfeld, M.; Chen, Y. Life-cycle analysis on

394

biodiesel production from microalgae: Water footprint and nutrients balance.

395

Bioresour. Technol. 2011, 102 (1), 159–165.

396

(7)

397 398

other applications: A review. Renew. Sustain. Energy Rev. 2010, 14 (1), 217–232. (8)

399 400

Mata, T. M.; Martins, A. A.; Caetano, N. S. Microalgae for biodiesel production and

Béchet, Q.; Shilton, A.; Guieysse, B. Full-scale validation of a model of algal productivity. Environ. Sci. Technol. 2014, 48 (23), 13826–13833.

(9)

Slegers, P. M. M.; van Beveren, P. J. M. J. M.; Wijffels, R. H. H.; van Straten, G.; van

401

Boxtel, a. J. B. J. B. Scenario analysis of large scale algae production in tubular

402

photobioreactors. Appl. Energy 2013, 105, 395–406.

403

(10)

404 405

Ritchie, R. J.; Larkum, A. W. D. Modelling photosynthesis in shallow algal production ponds. Photosynthetica 2012, 50 (4), 481–500.

(11)

Michels, M. H. A.; Slegers, P. M.; Vermuë, M. H.; Wijffels, R. H. Effect of biomass

ACS Paragon Plus Environment

18

Page 19 of 32

Environmental Science & Technology

406

concentration on the productivity of Tetraselmis suecica in a pilot-scale tubular

407

photobioreactor using natural sunlight. Algal Res. 2013, 4, 12–18.

408

(12)

409 410

Muñoz-Tamayo, R.; Mairet, F.; Bernard, O. Optimizing microalgal production in raceway systems. Biotechnol. Prog. 2013, 29 (2), 543–552.

(13)

Guieysse, B.; Béchet, Q.; Shilton, A. Variability and uncertainty in water demand and

411

water footprint assessments of fresh algae cultivation based on case studies from five

412

climatic regions. Bioresour. Technol. 2013, 128, 317–323.

413

(14)

Béchet, Q.; Shilton, A.; Park, J. B. K.; Craggs, R. J.; Guieysse, B. Universal

414

temperature model for shallow algal ponds provides improved accuracy. Environ. Sci.

415

Technol. 2011, 45 (8), 3702–3709.

416

(15)

Del Campo, J. A.; García-González, M.; Guerrero, M. G. Outdoor cultivation of

417

microalgae for carotenoid production: Current state and perspectives. Appl. Microbiol.

418

Biotechnol. 2007, 74 (6), 1163–1174.

419

(16)

Collet, P.; Hélias, A.; Lardon, L.; Ras, M.; Goy, R.-A.; Steyer, J.-P. Life-cycle

420

assessment of microalgae culture coupled to biogas production. Bioresour. Technol.

421

2011, 102 (1), 207–214.

422

(17)

Torres, C. M.; Ríos, S. D.; Torras, C.; Salvadó, J.; Mateo-Sanz, J. M.; Jiménez, L.

423

Microalgae-based biodiesel: A multicriteria analysis of the production process using

424

realistic scenarios. Bioresour. Technol. 2013, 147, 7–16.

425

(18)

Makareviciene, V.; Skorupskaite, V.; Andruleviciute, V. Biodiesel fuel from

426

microalgae-promising alternative fuel for the future: A review. Rev. Environ. Sci.

427

Biotechnol. 2013, 12 (2), 119–130.

428

(19)

429 430

Spolaore, P.; Joannis-Cassan, C.; Duran, E.; Isambert, A. Commercial applications of microalgae. J. Biosci. Bioeng. 2006, 101 (2), 87–96.

(20)

Gerbens-Leenes, W.; Hoekstra, A. Y.; van der Meer, T. H. The water footprint of

ACS Paragon Plus Environment

19

Environmental Science & Technology

431 432

bioenergy. Proc. Natl. Acad. Sci. U. S. A. 2009, 106 (25), 10219–10223. (21)

Béchet, Q.; Chambonnière, P.; Shilton, A.; Guizard, G.; Guieysse, B. Algal

433

productivity modeling: A step toward accurate assessments of full-scale algal

434

cultivation. Biotechnol. Bioeng. 2015, 112 (5), 987–996.

435

(22)

436 437

James, S. C.; Boriah, V. Modeling Algae Growth in an Open-Channel Raceway. J. Comput. Biol. 2010, 17 (7), 895–906.

(23)

James, S. C.; Janardhanam, V.; Hanson, D. T. Simulating pH effects in an algal-growth hydrodynamics model 1. J. Phycol. 2013, 49 (3), 608–615.

438 439

Page 20 of 32

(24)

Doucha, J.; Straka, F.; Lívanský, K. Utilization of flue gas for cultivation of microalgae

440

(Chlorella sp.) in an outdoor open thin-layer photobioreactor. J. Appl. Phycol. 2005, 17

441

(5), 403–412.

442

(25)

Doucha, J.; Lívanský, K. Productivity, CO2/O2 exchange and hydraulics in outdoor

443

open high density microalgal (Chlorella sp.) photobioreactors operated in a Middle and

444

Southern European climate. J. Appl. Phycol. 2006, 18 (6), 811–826.

445

(26)

Masojídek, J.; Kopecký, J.; Giannelli, L.; Torzillo, G. Productivity correlated to

446

photobiochemical performance of Chlorella mass cultures grown outdoors in thin-layer

447

cascades. J. Ind. Microbiol. Biotechnol. 2011, 38 (2), 307–317.

448

(27)

449 450

Moody, J. W.; McGinty, C. M.; Quinn, J. C. Global evaluation of biofuel potential from microalgae. Proc. Natl. Acad. Sci. U. S. A. 2014, 111 (23), 8691–8696.

(28)

Béchet, Q. Modelling the impact of temperature on microalgae productivity during

451

outdoor cultivation (PhD thesis), Massey University, Palmerston North, New Zealand,

452

2014.

453

(29)

Silva Benavides, A. M.; Torzillo, G.; Kopecký, J.; Masojídek, J. Productivity and

454

biochemical composition of Phaeodactylum tricornutum (Bacillariophyceae) cultures

455

grown outdoors in tubular photobioreactors and open ponds. Biomass and Bioenergy

ACS Paragon Plus Environment

20

Page 21 of 32

Environmental Science & Technology

456 457

2013, 54 (0), 115–122. (30)

Amer, L.; Adhikari, B.; Pellegrino, J. Technoeconomic analysis of five microalgae-to-

458

biofuels processes of varying complexity. Bioresour. Technol. 2011, 102 (20), 9350–

459

9359.

460

(31)

461 462

from microalgae in ponds. Bioresour. Technol. 2011, 102 (1), 50–56. (32)

463 464

Davis, R.; Aden, A.; Pienkos, P. T. Techno-economic analysis of autotrophic microalgae for fuel production. Appl. Energy 2011, 88 (10), 3524–3531.

(33)

465 466

Campbell, P. K.; Beer, T.; Batten, D. Life cycle assessment of biodiesel production

Gallagher, B. J. The economics of producing biodiesel from algae. Renew. Energy 2011, 36 (1), 158–162.

(34)

Jorquera, O.; Kiperstok, A.; Sales, E. A.; Embiruçu, M.; Ghirardi, M. L. Comparative

467

energy life-cycle analyses of microalgal biomass production in open ponds and

468

photobioreactors. Bioresour. Technol. 2010, 101 (4), 1406–1413.

469

(35)

Lardon, L.; Hélias, A.; Sialve, B.; Steyer, J.-P.; Bernard, O. Life-cycle assessment of

470

biodiesel production from microalgae. Environ. Sci. Technol. 2009, 43 (17), 6475–

471

6481.

472

(36)

Passell, H.; Dhaliwal, H.; Reno, M.; Wu, B.; Ben Amotz, A.; Ivry, E.; Gay, M.;

473

Czartoski, T.; Laurin, L.; Ayer, N. Algae biodiesel life cycle assessment using current

474

commercial data. J. Environ. Manage. 2013, 129, 103–111.

475

(37)

Resurreccion, E. P.; Colosi, L. M.; White, M. A.; Clarens, A. F. Comparison of algae

476

cultivation methods for bioenergy production using a combined life cycle assessment

477

and life cycle costing approach. Bioresour. Technol. 2012, 126, 298–306.

478

(38)

Richardson, J. W.; Johnson, M. D.; Outlaw, J. L. Economic comparison of open pond

479

raceways to photo bio-reactors for profitable production of algae for transportation

480

fuels in the Southwest. Algal Res. 2012, 1 (1), 93–100.

ACS Paragon Plus Environment

21

Environmental Science & Technology

481

(39)

Page 22 of 32

Rogers, J. N.; Rosenberg, J. N.; Guzman, B. J.; Oh, V. H.; Mimbela, L. E.; Ghassemi,

482

A.; Betenbaugh, M. J.; Oyler, G. A.; Donohue, M. D. A critical analysis of

483

paddlewheel-driven raceway ponds for algal biofuel production at commercial scales.

484

Algal Res. 2013, 4, 76–88.

485

(40)

Slade, R.; Bauen, A. Micro-algae cultivation for biofuels: Cost, energy balance,

486

environmental impacts and future prospects. Biomass and Bioenergy 2013, 53 (0), 29–

487

38.

488

(41)

Stephenson, A. L.; Kazamia, E.; Dennis, J. S.; Howe, C. J.; Scott, S. A.; Smith, A. G.

489

Life-cycle assessment of potential algal biodiesel production in the united kingdom: A

490

comparison of raceways and air-lift tubular bioreactors. Energy and Fuels 2010, 24 (7),

491

4062–4077.

492

(42)

493 494

Sun, A.; Davis, R.; Starbuck, M.; Ben-amotz, A.; Pate, R.; Pienkos, P. T. Comparative cost analysis of algal oil production for biofuels. Energy 2011, 36 (8), 5169–5179.

(43)

Wang, T.; Yabar, H.; Higano, Y. Bioresource Technology Perspective assessment of

495

algae-based biofuel production using recycled nutrient sources : The case of Japan.

496

Bioresour. Technol. 2013, 128, 688–696.

497

(44)

Zamalloa, C.; Vulsteke, E.; Albrecht, J.; Verstraete, W. The techno-economic potential

498

of renewable energy through the anaerobic digestion of microalgae. Bioresour.

499

Technol. 2011, 102 (2), 1149–1158.

500 501

ACS Paragon Plus Environment

22

Page 23 of 32

Environmental Science & Technology

502

Table 1: Pond depths and HRTs in recent assessments of algal growth for energy production

503

(NS: not specified, but accounted for in calculations; NC: not considered in the study) End product

Location/

HRT

Depth

Productivity

Water demandb

considered

climate

(d)

(m)

(g m-2 d-1)

(m3-m-2-yr-1)

30

Biofuel

NS

NC

0.3

24

NS

31

Biofuel

Australia

NC

0.3

30

7.7

16

Biogas

Europe

NC

0.3

25

28

32

Biofuel

US

NC

0.2

25

2.5

33

Biofuel

US

NC

NC

27

NS

34

Algae

NS

10

0.35

11

11.5

35

Biofuel

Europea

NC

0.3

20-30

NC

36

Biodiesel

NS

NC

NC

3-25

NS

37

Biofuel

NS

NC

NC

11.4-33.6

NC

38

Biofuel

US

NC

0.2

25

NS

39

Biofuel

US

10

0.3

15

0.13-1.1

40

Biofuel

NS

NC

NC

10-20

3-3.6

41

Biofuel

UK

10-13

0.28

36.5

0.015

42

Biofuel

NS

NS

0.2

20

NS

17

Biodiesel

Europe

NC

0.2

30

NS

43

Biofuel

Japan

NS

NS

38.9-43.9

NC

5

Biofuel

US

NC

0.3

3.8-14

0.0055-2

6

Biofuel

US

NC

NC

18

0.99

44

Biogas

NS

2

0.2

25.5-40

NS

Study

504

a

Mediterranean climate

505

b

Water demand was often not directly available from publications and was re-calculated from

506

available data provided by studies.

ACS Paragon Plus Environment

23

Environmental Science & Technology

Page 24 of 32

507

Table 2: Characteristics of the scenarios (fmax: maximum fraction of biomass lost through

508

respiration during the HRT optimization; NA: Not applicable)

509 HRT (d)

Pond depth (m)

fmax (%)

Standard conditions

10

0.25

NA

Depth control

10

Optimized

NA

HRT control

Optimized

0.25

50

Depth and HRT

Optimized

Optimized

50

Optimized

Optimized

Variablea

Scenario

control Depth and HRT control under constrained water demand 510

a

See Table 3 for details

ACS Paragon Plus Environment

24

Page 25 of 32

Environmental Science & Technology

511

Table 3: Yearly productivities and fractions of the photosynthetic output lost through

512

respiration f and fmax and in 'Standard conditions' and under water-demand control,

513

respectively. Standard conditions Climate

f (%)

Productivity

Water-demand control fmax (%)

(g-m-2-d-1)

Productivity (g-m-2-d-1)

Arid

55

12.5

NA

NA

Mediterranean

49

11.7

46

13.0

Subtropical

55

11.6

63

9.7

Tropical

51

11.2

20

16.4

Temperate

46

10.7

23

15.7

514 515

ACS Paragon Plus Environment

25

Environmental Science & Technology

Page 26 of 32

516

Figure 1: Top: Impact of pond depth on net biomass productivity; Middle: Impact of pond

517

depth on water demand; Bottom: Impact of pond depth on the number of hours when pond

518

temperature exceeds 35oC. Pond depth and HRT are held constant over the year (HRT = 10 d)

519

- Case of a Mediterranean climate (other climates are shown in S3).

520

Figure 2: Impact of HRT on net biomass productivity and water demand (WD) during algal

521

cultivation in open ponds when pond depth and HRT are held constant over the year (pond

522

depth = 0.25 m; White bar: net biomass productivity; Grey bar: respiration loss; Blue dashed-

523

line: water demand) - Case of a Mediterranean climate (other climates are shown in S4).

524

Figure 3: Net productivity and water demand of algal cultivation in open ponds operated

525

under standard conditions or different optimization strategies (Standard: HRT = 10 d; pond

526

depth = 0.25m; Depth control at 10 d HRT; HRT control at a depth of 0.25m; Depth and HRT

527

control; and depth and HRT control with constrained water demand).

528

Figure 4: Impact of pond depth and HRT on overheating risks, water demand, and respiration

529

loss during algal cultivation in open ponds (see section 2.1 for details).

530

Figure 5: Minimum pond depths ensuring pond temperature never exceeds 35oC for > 20

531

hours each month during depth control (HRT = 10 d).

532 533 534 535 536 537 538 539 540

ACS Paragon Plus Environment

26

Page 27 of 32

541

Environmental Science & Technology

Figure 1

542

ACS Paragon Plus Environment

27

Environmental Science & Technology

543 544

Page 28 of 32

Figure 2

545

ACS Paragon Plus Environment

28

Page 29 of 32

546

Environmental Science & Technology

Figure 3

547 548

ACS Paragon Plus Environment

29

Environmental Science & Technology

549

Page 30 of 32

Figure 4

550

ACS Paragon Plus Environment

30

Page 31 of 32

551

Environmental Science & Technology

Figure 5

552

ACS Paragon Plus Environment

31

Environmental Science & Technology

ACS Paragon Plus Environment

Page 32 of 32