Subscriber access provided by Kaohsiung Medical University
Energy and the Environment
Modeling Microalgae Productivity in IndustrialScale Vertical Flat Panel Photobioreactors Christian Hermann Endres, Arne Roth, and Thomas Bartholomaeus Brück Environ. Sci. Technol., Just Accepted Manuscript • DOI: 10.1021/acs.est.7b05545 • Publication Date (Web): 29 Mar 2018 Downloaded from http://pubs.acs.org on March 30, 2018
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 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 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.
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
84x47mm (300 x 300 DPI)
ACS Paragon Plus Environment
Environmental Science & Technology
1
TITLE
2
Modeling Microalgae Productivity in Industrial-Scale Vertical Flat Panel Photobioreactors
Page 2 of 32
3 4
AUTHORSHIP
5
Christian H. Endresa,b,*, Arne Rotha, Thomas B. Brückb
6 7 8
AFFILIATIONS
9
a
10
Bauhaus Luftfahrt, Willy-Messerschmitt-Str. 1, 82024 Taufkirchen, Germany, Tel. +49-
11
89307484948, E-mail:
[email protected] 12 13
b
14
Werner Siemens-Chair of Synthetic Biotechnology, Department of Chemistry, Technical
15
University of Munich, Lichtenbergstraße 4, 85748 Garching, Germany
16 17 18
ABSTRACT
19
Potentially achievable biomass yields are a decisive performance indicator for the economic
20
viability of mass cultivation of microalgae. In this study, a computer model has been
21
developed and applied to estimate the productivity of microalgae for large-scale outdoor
22
cultivation in vertical flat panel photobioreactors. Algae growth is determined based on
23
simulations of the reactor temperature and light distribution. Site-specific weather and
24
irradiation data are used for annual yield estimations in six climate zones. Shading and
25
reflections between opposing panels and between panels and the ground are dynamically
26
computed based on the reactor geometry and the position of the sun. The results indicate that 1 ACS Paragon Plus Environment
Page 3 of 32
Environmental Science & Technology
27
thin panels (≤ 0.05 m) are best suited for the assumed cell density of 2 g L-1 and that reactor
28
panels should face in north-south direction. Panel spacings of 0.4 – 0.75 m appear most
29
suitable for commercial applications. Under these preconditions, yields of around 10 kg m-2 a-
30
1
31
have to be expected, as extreme reactor temperatures limit overall productivity.
are possible for most locations in the U.S. Only in hot climates significantly lower yields
32 33
TOC ABSTRACT
34 35 36
INTRODUCTION
37
Microalgae have attracted substantial interest as feedstock for a broad range of value added
38
products, encompassing pharmaceuticals, pigments, proteins and biofuels
39
show high growth rates, by far exceeding those of terrestrial plants
40
cultivation does not depend on arable land and can be operated using seawater
41
waste water 11 as growth medium.
1–5
. Microalgae
1,6–9
. In addition, 10
or even
42
Various reactor designs for mass outdoor cultivation of microalgae are currently
43
pursued. In principle, microalgae can be cultivated in open systems or closed photobioreactors
44
(PBRs)
45
engineering, closed PBRs offer high controllability of growth conditions, albeit at the cost of
46
increased system complexity
12
. While open systems, particularly open ponds, tend to require less elaborate
9,13,14
. Vertical flat panel PBRs, an example of closed reactor
2 ACS Paragon Plus Environment
Environmental Science & Technology
Page 4 of 32
15–19
47
systems, represent a popular reactor concept for academic research
and commercial
48
activities 20–22 and has therefore been selected as cultivation system for our study.
49
Despite of the widely acknowledged potentials of microalgae as feedstock for the
50
production of fuels and other bioproducts, large-scale cultivation, harvesting and processing
51
of algae biomass are mostly immature technologies. With regard to the economic viability of
52
algae production, the achievable biomass yield represents a crucial performance indicator.
53
However, experiences on microalgal biomass yields are usually based on scientific
54
experiments at laboratory scale
55
for industrial-scale cultivation systems, due to the small number of operational production
56
plants and the lack of publicly-available information on these facilities.
23,24
, while only little knowledge exists about potential yields
57
In the absence of an empirical knowledge base on microalgae biomass yields,
58
computational simulations are indispensable to assess potential biomass productivities, to
59
screen potential cultivation sites and to evaluate optimum reactor geometries and reactor field
60
configurations
61
and on the temperature of the culture medium. Consequently, it is crucial for any
62
computational assessment of algae productivity to thoroughly model light and temperature as
63
a function of local external conditions 27–29.
25,26
. The growth of algae fundamentally depends on the availability of light
64
Only a limited number of simulations of outdoor cultivation of microalgae in large-
65
scale facilities has been reported to date. Most of this work has been devoted to open pond
66
systems, which are easier to describe mathematically in a computer model
67
few studies exist that cover yield estimations for closed PBRs under industrially relevant
68
cultivation conditions. In two such studies, algae growth was simulated in arrayed tubular
69
and flat panel PBRs 25. The light distribution was dynamically modeled based on the position
70
of the sun and the reactor geometry. Furthermore, shading caused by the reactors was taken
71
into account. A major drawback of these studies is the fact that the influence of the reactor
72
temperature on algae growth has been completely neglected. In contrast, in another study
30–34
, while only
35
3 ACS Paragon Plus Environment
Page 5 of 32
Environmental Science & Technology
26
73
temperature as well as light were dynamically simulated for a bubble column reactor
.
74
However, as only a single reactor was considered, effects like mutual shading and radiation
75
transfer were not taken into account. Consequently, the results of this study are only of
76
limited relevance for industrial applications.
77
In summary, to the best of our knowledge, no simulation of the productivity of
78
microalgae cultivated in arrayed closed PBRs has been reported to date that includes a
79
detailed temperature model. The work presented here closes this gap: For the first time, a
80
sophisticated temperature model is combined with an elaborate simulation of light availability
81
to assess potential biomass yields for microalgae cultivation in an arrayed set-up of vertical
82
flat panel PBRs. Shading effects and first-order reflections are taken into account; the light
83
propagation within the reactors is dynamically simulated based on the actual angles of
84
incidence and on the geometry and relative arrangement of reactor panels. Site-specific data
85
on weather and solar irradiation are used to simulate microalgae growth and to estimate
86
potential biomass yields for a complete year of cultivation at six U.S. locations in different
87
climate zones. Important design parameters for industrial-scale cultivation of microalgae in
88
vertical flat panels PBRs are discussed.
89
The presented productivity model provides a valuable tool to assess potential
90
microalgae yields at any given location under local growing conditions and to select the most
91
suitable set-up of the reactor panel field, based on the site-specific economic boundary
92
conditions.
93 94
METHODOLOGY
95
Examined Cultivation System. The examined cultivation system consists of vertical
96
flat plate PBRs that are arranged in parallel panel rows (see abstract graphic). Panels at the
97
edges of an array are neglected since they are exposed to higher levels of irradiation and are
98
thus not representative for the majority of the reactors. The loss of accuracy associated with 4 ACS Paragon Plus Environment
Environmental Science & Technology
Page 6 of 32
99
this simplification is only small, considering that even for a small field of 100 m × 100 m and
100
a panel width of 2 m, a panel distance of 0.5 m and a panel thickness of 0.05 m, only about
101
5% of all panels would be located at the edges of the field. Furthermore, only the back and
102
front surfaces of the reactor panel are considered for light capture and heat transfer. The areas
103
at the sides, the top and the bottom of the reactor are neglected as they are small in
104
comparison to the total surface area of the reactor. In addition, the respective areas are often at
105
least partially blocked by the ground (bottom), adjacent reactors (sides), equipment or an
106
attachment to a frame (top), further reducing light input through these areas. The impact of
107
neglecting these small areas on the reported outcomes is therefore considered small.
108
Reactor panels of different thicknesses are examined, while the width and height of the
109
reactors are kept constant at 2 and 1 m, respectively. However, yields generated for this
110
standard reactor can be transferred with good approximation to reactors of different
111
dimensions, given that the new reactor shows the same ratio of panel height to panel distance
112
and is identical with respect to panel thickness (Supporting Information, Section 2.2). For this
113
case, yields per panel are directly proportional to the reactor’s surface area (front and back).
114
Yields per panel of a reactor system that is 2 m high with a panel distance of 1 m is, for
115
example, twice as high as yields corresponding to reactors that are 1 m high with a panel
116
spacing of 0.5 m. Areal yield can be calculated from yields per panel by simply dividing the
117
latter by the occupied ground area ((panel distance + panel thickness) × panel width). Hereby,
118
the panel distance is defined as the distance between two opposing panel surfaces of adjacent
119
panel rows (see also length d in Figure S7, Supporting Information).
120
The culture medium is continuously mixed by pneumatic aeration. Temperature
121
gradients within the reactor are therefore assumed to be negligible. A constant cell
122
concentration of 2 g L-1 is assumed throughout this. Quasi-constant cell concentrations can,
123
for example, be reached by dynamically harvesting the algae according to the current biomass
124
productivity. 5 ACS Paragon Plus Environment
Page 7 of 32
Environmental Science & Technology
125 126
General Approach. The calculation of microalgae productivity for outdoor
127
cultivation plants comprises three steps (Figure 1). First, cultivation temperature is simulated
128
as a function of local climate, reactor geometry and solar irradiation. Second, based on the
129
reactor geometry and solar irradiation, the light distribution within the reactors is determined.
130
Third, microalgae productivities are obtained from the cultivation temperature and local light
131
intensities by applying a mechanistic growth model.
132
133 134
Figure 1. Basic approach for calculating algae biomass yields in closed PBRs.
135 136
Details Regarding Temperature and Light Simulation. The implemented 36
137
temperature model is described in detail in a previous publication
138
calculation of the heat balance in the reactors as a function of the geometry, distance and
139
orientation of reactor panels, of local climatic conditions and irradiation (including all first-
140
order reflections) and thermal radiation from the ground. Importantly, mutual shading is
141
calculated as a dynamic function of the sun’s position. The original temperature model has
142
been slightly adapted to integrate it into the comprehensive model presented here.
. It is based on the
6 ACS Paragon Plus Environment
Environmental Science & Technology
Page 8 of 32
143
Specifically, the fixed biomass fixation rate simplistically used in the original model has now
144
been replaced by a mechanistic growth model, as described in the Supporting Information
145
(Section S.1.1).
146
For an infinite array of flat panel PBRs, the light distribution does not change along
147
the panel rows. The light distribution therefore is only determined for two dimensions, i.e. the
148
cross section of a single panel (panel thickness × height). This cross section is divided by a
149
fine grid into cells of 10 mm in height and 0.5 mm in width. The light intensity is determined
150
at the center of each cell. Six components are taken into account to determine local light
151
intensities: direct and diffuse sunlight as well as the respective reflections of sunlight at the
152
panels and the ground. Basic assumptions and principles regarding the calculation of local
153
light intensities are described in the following. For more detailed information and underlying
154
equations, please refer to the Supporting Information (Section S1).
155
Light attenuation in the reactor is calculated using a modified version of the Lambert37
156
Beer law
157
usually independent of the concentration of the considered substance and represents only light
158
absorption. The modified law, however, also accounts for the scattering effect of algae. As a
159
result, the influence of the cell concentration has to be considered. For the simulation, an
160
extinction coefficient of 100 m2 kg-1 is used for a cell concentration of 2 g L-1. This value is in
161
accordance with experimentally derived values for Chlorella vulgaris at high cell densities 37.
162
. The expression “modified” here refers to the extinction coefficient, which is
Reflection losses at the panel surfaces are calculated by applying Snell’s law assuming
163
refractive indices of 1.0 for air 38, 1.5 for the panel wall (glass, plastic)
39–41
164
culture medium 38. The reflectivity of the ground was assumed to be 0.3 corresponding to dry
165
sandy soil without vegetation 42,43. Most natural soils show albedos between 0.1 and 0.45. The
166
selected value thus represents an intermediate with respect to this range. The impact of the
167
ground’s reflectivity on biomass yields is illustrated and discussed in the Supporting
and 1.3 for the
7 ACS Paragon Plus Environment
Page 9 of 32
Environmental Science & Technology
168
Information (Section S3.4). Shading of the panels is dynamically computed based on the
169
sun’s position, panel geometry and orientation of the panels.
170 171
Translation of Temperature and Light into Cell Growth. A wide variety of
172
mechanistic growth models exists. Research in this field has been comprehensively reviewed
173
28,44
174
can be used to determine algae productivity. For the simulation of industrial-scale outdoor
175
cultivation, however, several criteria should be met. Cultivation temperature has a strong
176
effect on algae growth. A respective productivity model should therefore take the cultivation
177
temperature into account. In contrast, the influence of nutrient concentrations can be
178
neglected for most commercial applications, assuming that nutrients are added in an optimum
179
amount. Incoming sunlight can be considered by using a spatially resolved light distribution
180
or by using average or incident light intensities. Growth models using the latter options,
181
however, do not allow the assessment of various reactor dimensions and reactor field
182
configurations as their utilization is generally limited to a specific reactor geometry. For the
183
present work, thus a model is selected that allows productivity simulations based on the
184
distribution of light that is individually modeled for the respective reactor geometry. Lastly, a
185
model for simulating outdoor cultivation should comprise a term for endogenous cell
186
respiration, to account for biomass losses in poorly illuminated areas of the reactor and during
187
the night.
188
; other growth models have been introduced recently 45,46. In principle, any growth model
The growth model introduced by Béchet et al.
37
includes all of the above mentioned
189
aspects and was therefore adopted for the present study. Specific growth parameters for the
190
microalgae Chlorella vulgaris were provided in the same publication. The respective wild
191
type strain (GenBank rbcL sequence: EF589154), which was isolated in New Zealand, is very
192
robust and grows over a wide range of temperatures. The strain is therefore an ideal candidate
8 ACS Paragon Plus Environment
Environmental Science & Technology
Page 10 of 32
193
for outdoor applications. Validations of the growth model were conducted in laboratory 37 and
194
in outdoor experiments 26.
195
The total productivity, P, of algae within a reactor panel is determined by integrating
196
local productivities, µloc, over the reactor volume, VR, and by multiplying them by the cell
197
concentration, X (Equation 1).
198 (1) 199 200
Iloc is the local light intensity and TR the reactor temperature. The calculation of µloc is
201
described in the Supporting Information (Section S1.3).
202
Béchet et al. reported positive cell growth for the examined Chlorella strain in the
203
temperature range between 0 and 42.1 °C. However, it was not reported to which extent the
204
cell culture is harmed at temperatures exceeding 42.1 °C. Lee et al. observed Chlorella-
205
cultures resisting peak temperatures of 49 °C in outdoor PBRs 47. Other eukaryotic microalgae
206
were reported to survive even higher temperatures 48. Based on these references, the following
207
assumptions are made for the applied productivity model: Algae cells survive temperatures
208
between 42.1 to 50 °C; however, biomass production is stopped for the time the algae are
209
exposed to these temperatures. If exceeding 50 °C, massive cell death is expected, resulting in
210
a collapse of the cell culture. As a consequence, production is put to a halt for seven days after
211
the occurrence of temperatures above 50 °C, representing the time algae cells need to recover
212
or a new culture needs to reach the original cell concentration. If temperatures exceed 50 °C
213
again while still in recovery, the phase of zero biomass production is extended by further
214
seven days.
215
Subzero cultivation temperatures are not covered by the model, as this case does not
216
represent a practical application. In the model, the cultivation temperature is therefore 9 ACS Paragon Plus Environment
Page 11 of 32
Environmental Science & Technology
217
artificially kept at 0 °C, even though the thermal balance could result in ice formation. When
218
conditions improve, the reactor temperature starts to rise again from 0 °C.
219 220
Selected Cultivation Sites and Input Data. Six U.S. locations have been selected as
221
cultivation sites for productivity simulation (Table 1). The selected locations represent distinct
222
climate zones that can also be found at other places worldwide. Meteorological and solar
223
irradiation data used for the simulation originates from the National Solar Radiation Data
224
Base
225
locations in different climate zones within the U.S. and provides the required data in a
226
consistent quality for all selected locations.
49
. The database is specifically intended for solar conversion systems covering various
227 228
Table 1. U.S. locations selected as cultivation sites for productivity simulations with
229
respective climate classification and average air temperature (Tavg, annual average; Tcold,
230
average coldest month; Thot, average hottest month). Location
Forks, WA
Climatea
Temperate; without dry
Tavg
Tcold
Thot
[°C]
[°C]
[°C]
9.6
5.0
15.2
10.6
-3.0
23.4
15.5
7.5
23.5
season; warm summer
Boston, MA
Cold; without dry season; warm/hot summer
Sacramento, CA
Temperate; dry and hot summer
Phoenix, AZ
Arid; desert; hot
23.8
11.7
35.6
New Orleans, LA
Temperate; without dry
20.4
9.9
30.0
23.1
21.7
24.6
season; hot summer Hilo, HI
231
a
Tropical; rainforest
According to Köppen-Geiger climate classification 50.
10 ACS Paragon Plus Environment
Environmental Science & Technology
Page 12 of 32
232 233
Computation Details. All models were developed, implemented and executed in
234
MATLAB® (The MathWorks®, Inc., Natick, MA). The calculations of algae productivity (as
235
well as of temperature and light distribution) were performed for every minute of a complete
236
year. The time needed for simulating a complete year of cultivation amounts to approximately
237
30 min (Intel Core i5 2.70 GHz, 4 GB RAM).
238 239
RESULTS AND DISCUSSION
240
Characterization of Light Distribution and Local Biomass Production. A cross
241
sectional cut of the reactor is generated to provide a visual impression of the light distribution
242
and corresponding local biomass production (Figure 2). For this purpose, an exemplary
243
reactor situated in Sacramento, CA is examined. Reactor panels face in north and south
244
direction. The distance between the panels is set to 0.5 m and the thickness of the panels to
245
0.05 m. The images show a day in late spring (day 100) at noon (12:00).
246
11 ACS Paragon Plus Environment
Page 13 of 32
Environmental Science & Technology
247 248
Figure 2. (A) Distribution of light and (B) biomass production in the PBRs at noon (12:00).
249
Different scalings of the vertical and horizontal axes are used to adequately display light and
250
productivity distributions (for both subfigures: location, Sacramento, CA; panel distance,
251
0.5 m; panel thickness, 0.05 m; orientation, north-south; considered day of the year, 100).
252 253
The light distribution is dominated by direct irradiation which enters the reactor on the
254
left side of the panel (Figure 2A). The significant change of light intensity 0.2 m above
255
ground is caused by shading from the opposing panel row. Diffuse irradiation is most intense
256
at the upper edges of the reactor; however, it is only visible in the right corner of the image, as
257
it is superimposed by direct irradiation in the left corner. Direct light reflected at the opposing
258
panel wall can be seen in the lower third on the right side of the panel. Other types of
259
irradiation cannot be recognized in the illustration as the respective light intensities are too
260
small. The respective irradiation types are nevertheless still considered for productivity
261
simulations. The major fraction of incoming sunlight is absorbed close to the panel wall and 12 ACS Paragon Plus Environment
Environmental Science & Technology
Page 14 of 32
262
its intensity quickly decreases when moving towards the center of the reactor. As a
263
consequence, large areas of the reactor are poorly illuminated. It therefore could be concluded
264
that algae growth is limited to the edges of the reactor and, furthermore, that the selected
265
panel thickness of 5 cm is far too large for efficient biomass production.
266
This assumption, however, has to be corrected when looking at the image showing the
267
local biomass production (Figure 2B). As can be seen, algae growth is not limited to the
268
edges, but stretches far deeper into the reactors. Even areas appearing not to be illuminated at
269
all can significantly contribute to the overall productivity. The reason for this lies in the
270
relation between growth and light intensity. Algae are very efficient at low irradiation values
271
and convert a relatively large fraction of incoming photons suitable for photosynthesis into
272
biomass 51. At high light intensities, algae usually enter a state of light saturation. At this state,
273
light conversion is very inefficient. Even higher light intensities may cause photoinhibition;
274
however, this aspect is not covered by the model.
275
Based on the presented results, it has to be noted that even light sources weakly
276
contributing to the overall solar energy received by the reactors can significantly enhance
277
overall productivity (see diffuse irradiation and reflected direct irradiation). Biomass
278
production is further not limited to highly irradiated areas close to the panel wall but stretches
279
far deeper into the reactor. Only in the center of the reactor light intensity is too weak and cell
280
respiration results in a net loss of biomass productivity (Figure 2B, area of negative net
281
productivity). For the presented example at the given point in time, the total loss of biomass
282
productivity in this specific area of the reactor amounts to 0.09 g min-1. 2.61 g min-1 are
283
produced in the parts of positive biomass production, resulting in an overall productivity of
284
2.52 g min-1 for the whole panel.
285
In terms of productivity, dark zones are not the only criterion for optimizing reactor
286
design and plant layout. Thick panels are less prone to temperature fluctuations and in
287
particularly reduce the risk of algae being exposed to extremely high reactor temperatures 13 ACS Paragon Plus Environment
Page 15 of 32
Environmental Science & Technology
288
during summer. Losses that result from operating reactors in an unfavorable temperature
289
regime may thus exceed those associated with to dark zones in the reactor. For these
290
considerations, however, biomass yields must be analyzed on a yearly basis with respect to
291
panel distance and thickness, as discussed in detail further below.
292 293
Impact of the Various Types of Irradiation on Microalgae Productivity. The
294
results described in the previous section indicate that even though direct irradiation may
295
dominate an irradiation profile, weaker types of irradiation may still have a significant
296
influence on the overall productivity. The impact of various irradiation types on the overall
297
productivity is therefore further examined (Table 2). One year of cultivation is selected as
298
time frame for the analysis. The reactor configuration is identical with the one described in
299
the previous section (location Sacramento, CA; panel distance, 0.05 m; panel thickness, 0.5
300
m; orientation, north-south).
301
In a first step, the shares of various irradiation types in the total amount of absorbed
302
solar energy are examined (Table 2, second column). For this purpose, a certain type of light
303
is neglected and the remaining received solar energy is quantified. The comparison clearly
304
shows that direct irradiation dominates the yearly energy input. Neglecting direct irradiation
305
consequently results in a reduction of 60% of captured sunlight. This is followed by diffuse
306
irradiation which covers around 30% of incoming sunlight. Reflections of direct and diffuse
307
light at opposing panels and the ground cover the residual 10%. Of the total available sunlight
308
(7.1 GJ panel-1 a-1) only about two thirds (4.4 GJ panel-1 a-1) are actually captured by the
309
reactors of the given geometry. Factors responsible for this loss, are light absorbed by the
310
ground as well as reflections at the ground/panels that are not captured by another panel. In
311
addition, the top surfaces of the reactors are not accounted for in the model as described in the
312
methodology section. A comprehensive analysis of captured light as a function of panel
313
spacing and thickness is presented in our previous publication 36. In general, it has to be noted 14 ACS Paragon Plus Environment
Environmental Science & Technology
Page 16 of 32
314
that the limited efficiency of light capture represents a drawback of vertical systems that must
315
be outweighed by their advantages compared to other reactor systems.
316 317
Table 2. Impact of neglecting specific types of irradiation on total captured solar energy and
318
annual biomass yield. Values in brackets indicate the relative decrease of captured solar
319
energy and productivity with respect to the case of regular irradiation (location, Sacramento,
320
CA; panel distance, 0.5 m; panel thickness, 0.05 m; orientation, north-south). Neglected irradiation type
Captured sunlight
Yieldsa
[-]
[GJ panel-1 a-1]
[kg panel-1 a-1]
None (regular irradiation)
4.4
9.0
Direct
1.8 (-60.0%)
4.0 (-55.6%)
Diffuse
3.1 (-29.7%)
3.0 (-66.8%)
Direct refl. by panels
4.4 (-2.0%)
8.6 (-5.1%)
Diffuse refl. by panels
4.4 (-1.2%)
8.8 (-2.8%)
Direct refl. by ground
4.2 (-5.0%)
8.1 (-9.9%)
Diffuse refl. by ground
4.4 (-2.1%)
8.5 (-5.2%)
All reflections
4.0 (-10.4%)
6.8 (-24.1%)
321
a
322
do not add up to 100%. For the same reason, production losses corresponding to the case that all reflections are
323
neglected cannot be determined from the sum of losses corresponding individual neglected reflections.
Please note that due to the non-linear relationship between light intensity and algae growth, production losses
324 325
In a second step, a sensitivity analysis is performed to quantify the significance of the
326
various types of irradiation on annual biomass yields (Table 2, third column). For the case of
327
regular irradiation, 4.4 GJ of captured sunlight (full solar spectrum) translates into 9.0 kg of
328
algae biomass. Assuming a heating value of 22 MJ kg-1, this corresponds to an overall
329
photosynthetic efficiency of 4.5%, which is in good correspondence to values cited for
330
microalgae
52,53
and plant leaves 54. The most striking result of the sensitivity analysis relates 15 ACS Paragon Plus Environment
Page 17 of 32
Environmental Science & Technology
331
to the different impacts of direct and diffuse irradiation: While far more light energy is
332
received by direct irradiation than by diffuse irradiation, diffuse irradiation has a stronger
333
influence on annual biomass yields. Neglecting diffuse irradiation leads to a drop of
334
productivity by approximately 67%, while neglecting direct sunlight only results in a drop of
335
56%. This can be explained by the fact that direct light typically illuminates only a small
336
fraction of the reactor’s surface, albeit with high local intensities of the incident radiation.
337
High light intensities, however, cannot be processed efficiently by microalgae due to light
338
saturation
339
expected. For comparison, diffuse irradiation represents a weaker source of energy, but the
340
incident light fully covers both sides of the panel. The broader spatial distribution of diffuse
341
light explains the strong influence on biomass yields.
51
. This is the reason why biomass yields related to direct sunlight are lower than
342
Reflected light illuminates areas in the reactor that are only poorly reached by direct
343
and diffuse irradiation, i.e. close to the ground. Reflections thus substantially contribute to a
344
more homogenous distribution of light in the reactor. Consequently, reflections have a
345
stronger impact on annual biomass yields (yield reductions of 24%) than might be anticipated
346
from their limited contribution to the total energy input (10% of captured sunlight).
347
Yield simulations for other panel distances and locations lead to similar results
348
(Supporting Information, Section S3.1). Therefore, it can be concluded that diffuse and
349
reflected light generally play a more important role for algae growth than could be expected
350
from their share in the total solar energy captured per reactor panel. This aspect must be kept
351
in consideration for future simulations and the process of integrated reactor and plant design.
352 353
Trade-Off Between Area-Specific Yield and Yields per Reactor Panel. The annual
354
yield with respect to the ground area (Figure 3, black lines) is a measure of how much
355
biomass can be produced on a certain area of land. The higher this area-specific yield, the less
356
land has to be acquired to produce a given quantity of biomass. Furthermore, piping systems 16 ACS Paragon Plus Environment
Environmental Science & Technology
Page 18 of 32
357
and infrastructure can be designed in a compact way for plants with high areal biomass
358
outputs, thus reducing investment costs. Alternatively, microalgae yield can be related to the
359
number of reactor panels in which it is achieved (Figure 3, green lines). A higher panel-
360
specific yield implies that fewer reactors are required for producing a certain amount of
361
biomass. By reducing the total number of required reactors, significant reductions of capital
362
costs can be expected.
363
364 365
Figure 3. Annual areal biomass yields (left y-axis, black lines, filled markers) and panel-
366
specific biomass yields (right y-axis, green lines, empty markers) as function of panel distance
367
and thickness. Square, circle, diamond and triangle markers indicate panel thicknesses of
368
0.025, 0.05, 0.1 and 0.15 m (location, Sacramento, CA; orientation, north-south).
369 370
High areal yields are achievable at small distances below 1 m. At larger panel
371
distances more and more light falls on the ground between the panel rows and is lost for
372
biomass production if not reflected back to the panels. Areal productivity therefore decreases
373
with increasing panel distance. For very small panel distances (< 0.3 m) a sharp decline in
374
areal productivity can be observed. Two major reasons are responsible for this decline. First, 17 ACS Paragon Plus Environment
Page 19 of 32
Environmental Science & Technology
375
even though more light with respect to the ground area is received by the panels, the available
376
light is distributed among a larger number of reactors. The individual reactor may therefore
377
receive only very little light turning large areas in the reactor virtually black. Algae cells in
378
these areas are not capable of producing biomass via photosynthesis but rather consume
379
biomass via cell respiration. As a result, areal yields are declining even though a large fraction
380
of incoming sunlight is captured by the panels. The second reason for the abrupt decline of
381
areal productivity is related to a simplification made in the model: Light input through the top
382
surface area is neglected for the simulation. The reason for this simplification is a reduced
383
model complexity and the fact that the top surface may be blocked for incoming light by
384
equipment or by the attachment to a frame. Small panel distances in the range of the panel
385
thickness, however, represent a case usually not relevant for industrial applications, as very
386
small panel distances correspond with a large number of reactors, implying high investment
387
costs.
388
With increasing panel distance the exposure to sunlight increases for the individual
389
reactor, typically resulting in higher productivities per panel. The impact of the increasing
390
panel distance on biomass yield per panel is strongest at small panel distances and loses
391
influence at larger panel distances. From a panel distance of approximately 2 m onwards,
392
productivity gains are often negligible. The major reason for this behavior is the fact that at
393
large panel distances only little additional light capture can be expected with increasing panel
394
distance. In addition, reactors can reach very high temperatures at large panel distances,
395
caused by a lack of shading
396
temperatures negatively impact biomass productivity of algae. As a consequence, maximum
397
productivities per panel are not necessarily found at the highest panel distance examined (5
398
m), where irradiation is strongest for an individual panel.
36
, in particular if located in a warm climate. Very high
399
For commercial applications, both, high areal yields and high yields per reactor unit,
400
are desirable. Therefore, a trade-off is necessary with respect to panel distance. Based on our 18 ACS Paragon Plus Environment
Environmental Science & Technology
Page 20 of 32
401
results, panel distances in the range of 0.4 to 0.75 m appear most promising. This range offers
402
also the highest potential in case of the other examined locations, as shown in the Supporting
403
Information (Figure S22). However, for a specific determination of the optimum
404
configuration of the production plant, a detailed techno-economic assessment is required,
405
based on, i.a., information regarding costs for land and reactor acquisition. Such an economic
406
analysis is beyond the scope of the present work.
407 408
Optimal Thickness and Orientation of the Reactor Panels for High Biomass
409
Outputs. When concentrating on the optimal panel thickness, thin panels (0.025 m, 0.05 m)
410
appear most favorable for Sacramento (square and circle markers, Figure 3) as well as for the
411
other locations examined in this work (Supporting Information, Figure S22). For the present
412
study, an algae concentration of 2 g L-1 with a corresponding extinction coefficient of 100 m2
413
kg-1 was assumed. Under these conditions, light cannot penetrate very far into the reactor
414
(light intensity is reduced by 95% after 0.015 m). Nevertheless, thicker panels can be of
415
interest, if focus lies on maximizing yields per panel. For this case, increased panel distances
416
are required to enable a sufficient light exposure of the panels and to minimize dark zones
417
within the cultivation medium, however at the expense of reduced areal productivity. At the
418
same time, thicker panels are less prone to overheating. When going on to greater panel
419
thicknesses (above 0.1 m), however, biomass production rates are strongly reduced due to
420
poor light availability in large parts of the reactor volume.
421
In this context it has to be noted that the correlation of biomass productivity with
422
thickness and distance of the reactor panels depends on the optical density of the algae
423
suspension, as given by the cell density and the specific extinction coefficient. While we
424
generally used a fixed optical density for our simulations, the impact of algae cell
425
concentration on annual biomass yields was examined in a sensitivity analysis for the case of
426
cultivation in Sacramento, CA (Supporting Information, Section S3.4). It shows that a higher 19 ACS Paragon Plus Environment
Page 21 of 32
Environmental Science & Technology
427
cell concentration of 3 g L-1 generally results in lower or unchanged yields, compared to the
428
standard concentration of 2 g L-1. In contrast, a reduced cell concentration of 1 g L-1 enables
429
higher yields in thicker panels (0.05 m and more), as a consequence of reducing dark zones
430
and the associated biomass degradation through cell respiration. For example, a rise by 120%
431
of the annual yield per panel is observed for a panel thickness of 0.1 m and a panel distance of
432
1 m. For future work in this field it would be interesting to analyze biomass productivities for
433
optimal combinations of cell density and panel thickness. Such an optimal combination would
434
depend on the location, the orientation of the panels, the season of the year and other
435
conditions. However, it has to be noted that lower cell concentrations can negatively affect
436
downstream processing, e.g. cell harvest. In any case, a careful techno-economic assessment
437
is required to identify the most suitable combination of panel thickness, panel distance and
438
optical density of the algae suspension for the given site-specific physical and economic
439
boundary conditions.
440
The influence of the orientation of reactor panels on biomass yields are examined for
441
all selected U.S. locations (Supporting Information, Figure S22). Reactors are assumed to face
442
east and west or north and south. The analysis shows that for the considered locations and
443
ranges of panel thickness and distances, a north-south orientation of reactors tends to result in
444
higher annual biomass yields. The reasons for this finding are complex and relate to the
445
duration of light exposure, local light intensities, temperature regimes and seasonal variations.
446
Even though a general superiority of a north-south orientation was observed for the studied
447
locations, a thorough assessment of the most suitable orientation is required for any specific
448
plant set-up and local conditions of a potential cultivation site.
449 450
Influence of Geographic Location. Knowledge of the influence of the geographic
451
location on biomass yield is particularly important in the planning phase of a commercial
452
cultivation plant. Thus, biomass yields were analyzed for six U.S. locations in different 20 ACS Paragon Plus Environment
Environmental Science & Technology
Page 22 of 32
453
climate zones (Table 1). Results of this analysis are presented in Table 3 for panel distances
454
between 0.4 and 0.75 m, a panel thickness of 0.025 m and north-south orientation of the
455
panels. These design parameters are considered most suitable for commercial algae
456
cultivation as described above. Biomass yields simulated for a wider spectrum of design
457
parameters can be found in the Supporting Information (Figure S21).
458 459
Table 3. Simulated algae biomass yields for different geographic locations and various panel
460
distances, d (panel thickness, 0.025 m; orientation, north-south). Location
Areal yield
Yield per panel
[kg m-2 a-1]
[kg panel-1 a-1]
d = 0.4 m
d = 0.5 m d = 0.75 m
d = 0.4 m
d = 0.5 m d = 0.75 m
Forks, WA
9.9
9.6
8.4
8.4
10.1
13.0
Boston, MA
10.3
9.4
7.7
8.8
9.9
12.0
Sacramento, CA
11.5
10.0
6.6
9.8
10.6
10.3
6.7
5.5
3.8
5.7
5.6
6.0
New Orleans, LA
10.1
9.3
6.9
8.5
9.8
10.7
Hilo, HI
11.3
10.7
8.1
9.6
11.3
12.6
Phoenix, AZ
461 462
With the exception of Phoenix, areal yields of about 10 kg m-2 a-1 in combination with
463
yields per panel of around 10 kg panel-1 a-1 can be achieved at all of the examined locations.
464
Biomass yields in Phoenix are substantially lower compared to other locations, as the hot and
465
arid climate leads to overheating of the reactors. With extremely high reactor temperatures,
466
production is brought to halt for elongated periods during summer, indicating that microalgae
467
cultivation in closed PBRs is not a suitable option in such climates (exemplary temperature
468
and productivity profiles are displayed in the Supporting Information, Section S3.3).
21 ACS Paragon Plus Environment
Page 23 of 32
Environmental Science & Technology
469
In contrast, it is astonishing that northern locations, such as Forks and Boston, offer
470
the potential for comparably high biomass yields. This observation is also linked to the local
471
prevailing climates. The Köppen-Geiger climate classification
472
Forks as temperate with warm summers, while the climate in the Boston area is classified as
473
cold with warm/hot summers. Under these climatic conditions, winter is not long enough to
474
severely reduce annual biomass output. During summer, however, northern locations such as
475
Boston and Forks experience long periods of daylight with reactors being exposed to
476
moderately warm ambient temperatures, posing a low risk of overheating of the cultivation
477
medium. In this context, it is important to point out that the used mechanistic growth model is
478
based on specific growth parameters of an algae strain that maintains relatively high growth
479
rates even at cultivation temperatures below 15 °C
480
growing under cold conditions, yield expectations have to be reduced for northern locations.
37
50
describes the climate in
. For strains that are less adapted to
481
The probably most suitable locations for algae cultivation correspond with tropical
482
(Hilo, HI) and Mediterranean climates (Sacramento, CA) with potential biomass yields
483
exceeding 11 kg m-2 a-1.
484 485
Implications for Commercial Algae Cultivation in PBRs. In the present work,
486
potential biomass outputs for microalgae grown in vertical flat panel PBRs were analyzed for
487
a wide variety of plant designs and climate zones. Realistic large-scale outdoor cultivation
488
conditions were simulated in high detail.
489
The results show that areal microalgae yields of around 10 kg m-2 a-1 (100 t ha-1 a-1)
490
and more can be potentially achieved. This by far exceeds biomass yields of land-based
491
energy crops, for example of fast growing trees cultivated in short-rotation coppice
492
plantations in Europe (5 to 18 t ha-1 a-1) 55, North America (10 to 12 t ha-1 a-1) 56 or Brazil (20
493
to above 30 t ha-1 a-1) 57. To make use of these potentials, producers should limit the thickness
494
of the reactor panels to 0.05 m and use a panel spacing of 0.4 to 0.75 m (other panel distances 22 ACS Paragon Plus Environment
Environmental Science & Technology
Page 24 of 32
495
might be required if the reactor height deviates from the standard height of 1 m used
496
throughout this study; please refer to the Methodology section). Our results also indicate that
497
a north-south orientation of reactor panels is preferable to an east-west orientation for most
498
locations.
499
The reactor temperature is critical for outdoor biomass generation. Even in moderately
500
warm areas, reactor temperature can repeatedly exceed 40 °C during one year of cultivation.
501
Thus, it is mandatory to focus on temperature-robust strains or even thermophilic species. Hot
502
and arid areas are often considered promising for algae cultivation, as they offer cheap land
503
that is not used for food production, coupled with high levels of solar irradiation. Our results,
504
however, indicate that overheating of the reactors may substantially reduce biomass yields in
505
these areas. Active cooling of reactors, e.g. spray cooling, represents an option to avoid
506
overheating and mitigate yield losses. However, the installation of heat exchangers, pipes and
507
pumps would significantly impact the economic viability. Pumping of the required volumes
508
of cooling medium would affect the energy balance of cultivation, and the increased
509
consumption of fresh water would pose severe environmental and social risks, particularly in
510
arid regions where water is scarce and the need for cooling is most urgent. We analyzed the
511
effect of active cooling, limiting reactor temperatures to 40 °C, for the locations Sacramento,
512
CA, and Phoenix, AZ (Supporting Information, Section 3.6). The analysis shows that yields in
513
Phoenix, representing a location in a hot and arid region, can be strongly increased through
514
cooling especially when maximizing yields per panel. The effect of cooling is less
515
pronounced for Sacramento, as the risk of overheating is substantially lower at that location.
516
In the context of the potential benefit of cooling, we note that conclusions regarding the
517
environmental, economic and social sustainability of cooling can only be drawn on the basis
518
of a detailed life cycle analysis and cost assessment, which is beyond the scope of the present
519
work.
23 ACS Paragon Plus Environment
Page 25 of 32
Environmental Science & Technology
520
As the developed simulation tool has been shown to be a valuable asset in the process
521
of integrated plant design, future work should be directed to the refinement of the model. In
522
this context, the model should be extended to cover further promising reactor types and algae
523
species. Most importantly, the simulation tool should be validated against experimental data.
524
We did not have access to a suitable reactor array in the course of the present study. An
525
experimental validation of the model, however, would add to the value of the generated
526
results and may lead to further refinements. The experimental validation should consequently
527
be the next step in future research efforts.
528 529
ASSOCIATED CONTENT
530
MATLAB model can be obtained from the corresponding author upon request; Supporting
531
Information: Description of the productivity model; validation; additional results
532 533
ACKNOWLEDGEMENTS
534
We gratefully acknowledge the financial support by the German Federal Ministry of
535
Education and Research (Project: Advanced Biomass Value, 03SF0446C) and the support
536
granted by the TUM Graduate School. We further thank Christoph Falter, Valentin Batteiger
537
and Andreas Sizmann for their advice and many fruitful discussions.
538 539
REFERENCES
540
(1)
Chisti, Y. Biodiesel from microalgae. Biotechnol. Adv. 2007, 25 (3), 294–306.
541
(2)
Guedes, A. C.; Amaro, H. M.; Malcata, F. X. Microalgae as sources of high added-
542
value compounds – a brief review of recent work. Biotechnol. Prog. 2011, 27 (3), 597–
543
613.
544
(3)
Harun, R.; Singh, M.; Forde, G. M.; Danquah, M. K. Bioprocess engineering of 24 ACS Paragon Plus Environment
Environmental Science & Technology
545
microalgae to produce a variety of consumer products. Renew. Sustain. Energy Rev.
546
2010, 14 (3), 1037–1047.
547
(4)
548 549
Norsker, N.; Barbosa, M. J.; Vermuë, M. H.; Wijffels, R. H. Microalgal production – a close look at the economics. Biotechnol. Adv. 2011, 29 (1), 24–27.
(5)
550 551
Clark, J. H.; Luque, R.; Matharu, A. S. Green chemistry, biofuels, and biorefinery. Annu. Rev. Chem. Biomol. Eng. 2012, 3 (1), 183–207.
(6)
Mutanda, T.; Ramesh, D.; Karthikeyan, S.; Kumari, S.; Anandraj, A.; Bux, F.
552
Bioprospecting for hyper-lipid producing microalgal strains for sustainable biofuel
553
production. Bioresour. Technol. 2011, 102 (1), 57–70.
554
Page 26 of 32
(7)
Scott, S. A.; Davey, M. P.; Dennis, J. S.; Horst, I.; Howe, C. J.; Lea-Smith, D. J.;
555
Smith, A. G. Biodiesel from algae: challenges and prospects. Curr. Opin. Biotechnol.
556
2010, 21 (3), 277–286.
557
(8)
558 559
Wijffels, R. H.; Barbosa, M. J. An outlook on microalgal biofuels. Science 2010, 329 (5993), 796–799.
(9)
Tredici, M. R. Mass production of microalgae: photobioreactors. In Handbook of
560
microalgal culture: biotechnology and applied phycology; Richmond, A., Ed.;
561
Blackwell Science Ltd: Oxford, Ames, Carlton, 2004; pp 178–214.
562
(10)
Parmar, A.; Singh, N. K.; Pandey, A.; Gnansounou, E.; Madamwar, D. Cyanobacteria
563
and microalgae: a positive prospect for biofuels. Bioresour. Technol. 2011, 102 (22),
564
10163–10172.
565
(11)
566 567 568
Pittman, J. K.; Dean, A. P.; Osundeko, O. The potential of sustainable algal biofuel production using wastewater resources. Bioresour. Technol. 2011, 102 (1), 17–25.
(12)
Apel, A. C.; Weuster-Botz, D. Engineering solutions for open microalgae mass cultivation and realistic indoor simulation of outdoor environments. Bioprocess 25 ACS Paragon Plus Environment
Page 27 of 32
Environmental Science & Technology
569 570
Biosyst. Eng. 2015, 38 (6), 995–1008. (13)
571 572
Algae for biofuels and energy; Borowitzka, M. A., Moheimani, N. R., Eds.; Springer Netherlands: Dordrecht, 2013.
(14)
Olivieri, G.; Salatino, P.; Marzocchella, A. Advances in photobioreactors for intensive
573
microalgal production: configurations, operating strategies and applications. J. Chem.
574
Technol. Biotechnol. 2014, 89 (2), 178–195.
575
(15)
Zittelli, G. C.; Rodolfi, L.; Bassi, N.; Biondi, N.; Tredici, M. R. Photobioreactors for
576
microalgal biofuel production. In Algae for biofuels and energy; Borowitzka, M. A.,
577
Moheimani, N. R., Eds.; Springer Netherlands: Dordrecht, 2013; pp 115–131.
578
(16)
Zhang, K.; Kurano, N.; Miyachi, S. Outdoor culture of a cyanobacterium with a
579
vertical flat-plate photobioreactor: effects on productivity of the reactor orientation,
580
distance setting between the plates, and culture temperature. Appl. Microbiol.
581
Biotechnol. 1999, 52 (6), 781–786.
582
(17)
Sierra, E.; Acién, F. G.; Fernández, J. M.; García, J. L.; González, C.; Molina, E.
583
Characterization of a flat plate photobioreactor for the production of microalgae. Chem.
584
Eng. J. 2008, 138 (1–3), 136–147.
585
(18)
Rodolfi, L.; Chini Zittelli, G.; Bassi, N.; Padovani, G.; Biondi, N.; Bonini, G.; Tredici,
586
M. R. Microalgae for oil: strain selection, induction of lipid synthesis and outdoor mass
587
cultivation in a low-cost photobioreactor. Biotechnol. Bioeng. 2009, 102 (1), 100–112.
588
(19)
589
Richmond, A.; Cheng-Wu, Z. Optimization of a flat plate glass reactor for mass production of Nannochloropsis sp. outdoors. J. Biotechnol. 2001, 85 (3), 259–269.
590
(20)
Algenol LLC http://www.algenol.com/ (accessed Oct 16, 2017).
591
(21)
Subitec GmbH http://subitec.com/de (accessed Oct 16, 2017).
592
(22)
Tredici, M. R.; Rodolfi, L.; Biondi, N.; Bassi, N.; Sampietro, G. Techno-economic 26 ACS Paragon Plus Environment
Environmental Science & Technology
Page 28 of 32
593
analysis of microalgal biomass production in a 1-ha Green Wall Panel (GWP®) plant.
594
Algal Res. 2016, 19, 253–263.
595
(23)
Quinn, J. C.; Davis, R. The potentials and challenges of algae based biofuels: a review
596
of the techno-economic, life cycle, and resource assessment modeling. Bioresour.
597
Technol. 2015, 184, 444–452.
598
(24)
Sniffen, K. D.; Price, J. R.; Sales, C. M.; Olson, M. S. Influence of scale on biomass
599
growth and nutrient removal in an algal-bacterial leachate treatment system. Environ.
600
Sci. Technol. 2017, 51 (22), 13344–13352.
601
(25)
602 603
for flat panel photobioreactors. Appl. Energy 2011, 88, 3342–3353. (26)
604 605
Slegers, P. M.; Wijffels, R. H.; van Straten, G.; van Boxtel, A. J. B. Design scenarios
Béchet, Q.; Shilton, A.; Guieysse, B. Full-scale validation of a model of algal productivity. Environ. Sci. Technol. 2014, 48 (23), 13826–13833.
(27)
Carvalho, A. P.; Malcata, F. X. Kinetic modeling of the autotrophic growth of Pavlova
606
lutheri: study of the combined influence of light and temperature. Biotechnol. Prog.
607
2003, 19 (4), 1128–1135.
608
(28)
Béchet, Q.; Shilton, A.; Guieysse, B. Modeling the effects of light and temperature on
609
algae growth: state of the art and critical assessment for productivity prediction during
610
outdoor cultivation. Biotechnol. Adv. 2013, 31 (8), 1648–1663.
611
(29)
612 613
Bernard, O.; Rémond, B. Validation of a simple model accounting for light and temperature effect on microalgal growth. Bioresour. Technol. 2012, 123, 520–527.
(30)
Wigmosta, M. S.; Coleman, A. M.; Skaggs, R. J.; Huesemann, M. H.; Lane, L. J.
614
National microalgae biofuel production potential and resource demand. Water Resour.
615
Res. 2011, 47 (3), 1–13.
616
(31)
Béchet, Q.; Shilton, A.; Guieysse, B. Maximizing productivity and reducing 27 ACS Paragon Plus Environment
Page 29 of 32
Environmental Science & Technology
617
environmental impacts of full-scale algal production through optimization of open
618
pond depth and hydraulic retention time. Environ. Sci. Technol. 2016, 50 (7), 4102–
619
4110.
620
(32)
Orfield, N. D.; Keoleian, G. a.; Love, N. G. A GIS based national assessment of algal
621
bio-oil production potential through flue gas and wastewater co-utilization. Biomass
622
and Bioenergy 2014, 63, 76–85.
623
(33)
Sukenik, A.; Levy, R. S.; Levy, Y.; Falkowski, P. G.; Dubinsky, Z. Optimizing algal
624
biomass production in an outdoor pond: a simulation model. J. Appl. Phycol. 1991, 3
625
(3), 191–201.
626
(34)
627 628
Ritchie, R. J.; Larkum, A. W. D. Modelling photosynthesis in shallow algal production ponds. Photosynthetica 2012, 50 (4), 481–500.
(35)
Slegers, P. M.; van Beveren, P. J. M.; Wijffels, R. H.; van Straten, G.; van Boxtel, A. J.
629
B. Scenario analysis of large scale algae production in tubular photobioreactors. Appl.
630
Energy 2013, 105, 395–406.
631
(36)
Endres, C. H.; Roth, A.; Brück, T. B. Thermal reactor model for large-scale algae
632
cultivation in vertical flat panel photobioreactors. Environ. Sci. Technol. 2016, 50 (7),
633
3920–3927.
634
(37)
Béchet, Q.; Chambonnière, P.; Shilton, A.; Guizard, G.; Guieysse, B. Algal
635
productivity modeling: a step toward accurate assessments of full-scale algal
636
cultivation. Biotechnol. Bioeng. 2015, 112 (5), 987–996.
637
(38)
638 639 640
Hecht, E. Die Ausbreitung des Lichts. In Optik; Oldentburg Verlag: München, 2009; pp 151–251.
(39)
Meschede, D. Lichtstrahlen. In Optik, Licht und Laser; Vieweg + Teubner: Wiesbaden, 2008; pp 1–33. 28 ACS Paragon Plus Environment
Environmental Science & Technology
641
(40)
642 643
Filmetrics. Brechungsindex von PE, Polyethene http://www.filmetrics.de/refractiveindex-database/Polyethylene/PE-Polyethene (accessed Oct 16, 2017).
(41)
GoodFellow Corporation. Polyethylene - Low Density (LDPE)
644
http://www.goodfellow.com/E/Polyethylene-Low-Density.html (accessed Oct 16,
645
2017).
646
(42)
Soil physics companion, 1st ed.; Warrick, A. W., Ed.; CRC Press: Boca Raton, 2001.
647
(43)
Muneer, T. Solar radiation and daylight models, 2nd ed.; Elsevier Butterworth-
648 649
Heinemann: Burlington, 2004. (44)
650 651
Lee, E.; Jalalizadeh, M.; Zhang, Q. Growth kinetic models for microalgae cultivation: a review. Algal Res. 2015, 12, 497–512.
(45)
Huesemann, M.; Crowe, B.; Waller, P.; Chavis, A.; Hobbs, S.; Edmundson, S.;
652
Wigmosta, M. A validated model to predict microalgae growth in outdoor pond
653
cultures subjected to fluctuating light intensities and water temperatures. Algal Res.
654
2016, 13, 195–206.
655
(46)
Huesemann, M.; Dale, T.; Chavis, A.; Crowe, B.; Twary, S.; Barry, A.; Valentine, D.;
656
Yoshida, R.; Wigmosta, M.; Cullinan, V. Simulation of outdoor pond cultures using
657
indoor LED-lighted and temperature-controlled raceway ponds and Phenometrics
658
photobioreactors. Algal Res. 2017, 21, 178–190.
659
Page 30 of 32
(47)
Lee, T. M.; Tseng, Y. F.; Cheng, C. L.; Chen, Y. C.; Lin, C. S.; Su, H. Y.; Chow, T. J.;
660
Chen, C. Y.; Chang, J. S. Characterization of a heat-tolerant Chlorella sp. GD mutant
661
with enhanced photosynthetic CO2 fixation efficiency and its implication as lactic acid
662
fermentation feedstock. Biotechnol. Biofuels 2017, 10 (1), 1–12.
663 664
(48)
Kobayashi, Y.; Harada, N.; Nishimura, Y.; Saito, T.; Nakamura, M.; Fujiwara, T.; Kuroiwa, T.; Misumi, O. Algae sense exact temperatures: small heat shock proteins are 29 ACS Paragon Plus Environment
Page 31 of 32
Environmental Science & Technology
665
expressed at the survival threshold temperature in cyanidioschyzon merolae and
666
chlamydomonas reinhardtii. Genome Biol. Evol. 2014, 6 (10), 2731–2740.
667
(49)
National Renewable Energy Laboratory (NREL). National solar radiation data base
668
(1991 – 2005 update: typical meteorological year 3)
669
http://rredc.nrel.gov/solar/old_data/nsrdb/ (accessed Oct 16, 2017).
670
(50)
671 672
Peel, M. C.; Finlayson, B. L.; McMahon, T. A. Updated world map of the KöppenGeiger climate classification. Hydrol. Earth Syst. Sci. 2007, 11, 1633–1644.
(51)
Masojídek, J.; Koblížek, M.; Torzillo, G. Photosynthesis in microalgae. In Handbook
673
of microalgal culture: biotechnology and applied phycology; Richmond, A., Ed.;
674
Blackwell Science Ltd.: Oxford, Ames, Carlton, 2004; pp 20–39.
675
(52)
676 677
Norsker, N. H.; Barbosa, M. J.; Vermuë, M. H.; Wijffels, R. H. Microalgal production a close look at the economics. Biotechnol. Adv. 2011, 29 (1), 24–27.
(53)
Janssen, M.; Tramper, J.; Mur, L. R.; Wijffels, R. H. Enclosed outdoor
678
photobioreactors: light regime, photosynthetic efficiency, scale-up, and future
679
prospects. Biotechnol. Bioeng. 2003, 81 (2), 193–210.
680
(54)
681 682
Cambridge, 1999. (55)
683 684
Hall, D. O.; Rao, K. K. Photosynthesis, 6th ed.; Cambridge University Press:
Dimitriou, I.; Rutz, D. Sustainable short rotation coppice; Rutz, D., Ed.; WIP Renewable Energies: München, 2015.
(56)
Guidi, W.; Pitre, F. E.; Labrecque, M. Short-rotation coppice of willows for the
685
production of biomass in Eastern Canada. In Biomass now – sustainable growth and
686
use; Matovic, M. D., Ed.; InTech, 2013; pp 421–448.
687 688
(57)
Stape, J. L.; Binkley, D.; Ryan, M. G.; Fonseca, S.; Loos, R. A.; Takahashi, E. N.; Silva, C. R.; Silva, S. R.; Hakamada, R. E.; de A. Ferreira, J. M.; et al. The Brazil 30 ACS Paragon Plus Environment
Environmental Science & Technology
689
Eucalyptus potential productivity project: influence of water, nutrients and stand
690
uniformity on wood production. For. Ecol. Manage. 2010, 259 (9), 1684–1694.
Page 32 of 32
691
31 ACS Paragon Plus Environment