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