Modeling Microalgae Productivity in Industrial ... - ACS Publications

Mar 29, 2018 - Environmental Science & Technology .... Panel spacings of 0.4–0.75 m at a panel height of 1 m appear most suitable .... Béchet et al...
1 downloads 0 Views 724KB Size
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