3.1. General cultivation condition
The temperature range (22°C to 28°C) is conducive to the growth of Chlorella vulgaris (Dahmani et al. 2016), as indicated in Table 3, although an optimum temperature was suggested in the range 15°C to 20°C (Zhang et al. 2020). Exceeding a temperature of 30°C might render the conditions less favorable for microalgae growth (González-Camejo et al. 2019). Media temperature has an impact the on solubility of CO2, which serves as the primary carbon source (Al-Anezi et al. 2008).
Table 3. Environmental parameters measured periodically during the cultivation
Treatment
|
Measured parameters
|
Temperature (°C)
|
Salinity (ppt)
|
pH
|
Control (BBM)
|
22-28
|
20-25
|
7.3 - 8.1
|
A1
|
22-28
|
20-23
|
5.5 - 5.8
|
A2
|
22-28
|
20-24
|
5.5 - 6.0
|
A3
|
22-28
|
20-26
|
5.6 - 6.2
|
B1
|
22-28
|
20-23
|
5.5 - 5.8
|
B2
|
22-28
|
20-24
|
5.5 - 6.4
|
B3
|
22-28
|
20-26
|
5.5 - 6.5
|
C1
|
22-28
|
20-23
|
5.2 - 6.0
|
C2
|
22-28
|
20-24
|
5.6 - 6.3
|
C3
|
22-28
|
20-26
|
5.5 - 6.2
|
The pH levels of the growth media ranged from 5.2 to 6.3 for HF and from 7.3 to 8.1 for the control. All tested media were in acidic pH, with higher HF concentrations tended to give higher pH value. No pH buffer was applied to the media with an intention of observing the authentic response of microalgae, thereby reducing the additional cultivation expenses, particularly in economically disadvantaged regions. Numerous studies have indicated a preference for Chlorella vulgaris for alkaline pH levels (Choi and Lee 2011; Yu et al. 2018), with few reports of acidic pH (Sarker and Salam 2020). The pH value observed in BBM demonstrated a more ideal for promoting the growth of Chlorella vulgaris.
The cultivation commenced at 20 ppt and gradually rose to 26 ppt. Chlorella vulgaris is known to thrive in saline environments, with the optimal salinity range from 25 to 30 ppt (Adenan et al. 2013; Pandit et al. 2017). Evaporation due to constant aeration during the cultivation process might be the reason for this salinity increase.
3.2. Growth of the microalgae
Figure 2 illustrates the appearance of the microalgae during different time intervals for treatment in A3, B3, and C3. On day 30, microalgae exhibited a notably high cell density. Among treatments in Concentration-3, A3 entered the death phase earlier than the other two.
Figure 3 illustrates the growth of Chlorella vulgaris under varying HF ratios and concentrations with their non-linear fitting included. Microalgae cultivation in BBM lasted for 19 days, with the highest density observed on day 16, which was slightly different from those in treatments A2, B2, and C2. The shortest time required to reach the highest density was observed in treatments A1, B1 and C1 (12 days), while longer ones were observed in all treatments of Concentration-3.
The non-linear fitting using the Gompertz model conducted on the data shows the correlation coefficient (R2) ranged from 0.927 to 0.991, with an average R2 value of each treatment was 0.961. All parameter values in Equation 1 were obtained by the software.
In general, the duration of microbial growth is notably impacted by the concentration of the growth media, regardless of the specific HF ratio employed. Higher concentrations of the added HF extended the cultivation period prior to entering the death phase, as evidenced by the extended duration required to reach peak density in each HF ratio. Cultivation in BBM media continued for 19 days, with the peak biomass concentration observed on day 16. This pattern only slightly differed from the results in treatments A2, B2, and C2, where the average time ranged from 15 to 17 days to reach the peak biomass concentration. Conversely, the shortest time required to reach peak cell density was observed in treatments A1 and C1, occurring on day 12, with a one-day delay for treatment B1.
Comparing the growth rate value resulting from Equation 1, it shows that the growth rate of microalgae cultivated in control media (0.166.day-1) exceeded that in most of the tested media, except for treatments B1, B2, and C1 (Figure 4). While examining individual treatments, it was noted that the highest growth rate was observed in one of the B1 samples (0.313 day-1). However, when having a look at the average growth rate in each treatment, B2 shows slightly higher than that of B1, with both showing no significant difference. Specifically, the growth rate for B1 was measured at 0.2625, while B2 had a growth rate of 0.2644. A long treatment duration tends to have a negative impact on the overall growth rate, as evidenced by much lower growth rates in all treatments at Concentration-3.
By applying the ANOVA test (Table 4), ratio of the HF significantly influences the growth rate of the microalgae, which is shown by P Value of 2.09E-04. Although the overall ANOVA test results shows that the concentration and the interaction between the two tested parameters did not significantly influence the growth rate, the Tukey’s post-hoc test comparing a more detail comparison within the same factor shows that the growth rate of all Concentration-3 treatment shows significantly lower, as can be conclude visually by the data. The detail of this post-hot test can be found in the Supplement 2.
Table 4. Two-way ANOVA test for the growth rate
|
DF
|
Sum of Squares
|
Mean Square
|
F Value
|
P Value
|
Ratio
|
3
|
0.0534
|
0.0178
|
11.3318
|
2.09E-04
|
Concentration
|
2
|
0.0067
|
0.0033
|
2.1309
|
0.1477
|
Interaction
|
6
|
0.0179
|
0.003
|
1.8928
|
0.1374
|
Model
|
11
|
0.0994
|
0.009
|
5.7511
|
5.77E-04
|
Error
|
18
|
0.0283
|
0.0016
|
--
|
--
|
Corrected Total
|
29
|
0.1277
|
--
|
--
|
--
|
Figure 5 illustrates that the earliest inflection point was observed in treatment B2, occurring at 3.4 days, followed closely by treatment B1 at 3.6 days. Conversely, all Concentration-3 treatments exhibited the latest times to reach the inflection point in each HF ratio, with values slightly exceeding 10 days. Given that the inflection point signifies the phase of fastest microalgal growth, making it a crucial indicator in microalgal culture for determining the optimal time to harvest cell biomass, particularly when biomass is the primary product (Richmond 2003).
3.3. Biomass yield
The calibration graph presenting DBW shows the linear regression equation y = 3.1388x + 1.9755 with an R2 value of 0.974 (Supplement 1). Figure 6 illustrates the biomass yield at its highest density. Treatment B3 exhibits the highest biomass yield, approximately 9.75g.l-1. Unsurprisingly, all treatments in Concentration-3 (i.e. A3, B3, and C3) demonstrate distinctively higher biomass yield, which can be attributed to their extended cultivation period as a result of more available nutrients. In a comparable time to reach the peak density of the BBM (15-16 days), only treatment A2 slightly surpassed the biomass yield of the control. The two-way ANOVA followed by Tukey's post hoc test concluded that while the ratio did not exert a significant impact, the concentration significantly affected the biomass yield (Supplement 3).
3.4. Proximate value (protein, lipid, carbohydrate, and ash)
The proximate analysis of microalgal cells identifies major components in dried biomass, essential for assessing nutritional value and suitability for consumption. This analysis typically determines key nutrients: carbohydrate, protein, and lipid. In this study, the proximate analysis was conducted only for BBM and Concentration-3 of each HF ratio.
Chlorella vulgaris cultivated in Ratio-B of HF demonstrates the highest protein content, accounting for approximately 33.93% of the dry weight of biomass (Figure 7). This figure does not significantly differ from those of other treatments, which ranged between 30% and 34%. Microalgae cultivated in BBM exhibit the highest lipid (35.30%), slightly surpassing that of microalgae cultivated in HF. As for carbohydrates, microalgae cultivated in BBM contained the lowest percentage (28%), contrasting with other treatments that ranged from 35% to 42%. Similarly, the ash content in the various treatments exhibits relatively comparable percentages.
3.5. Cultivation cost
The cultivation cost discussed in this study refers to the expenses involved in producing one kilogram of dried microalgae biomass. This cost estimation is tailored for SCP production within households, small and medium-scale enterprises. This estimation was originally conducted in the local currency of Indonesia (IDR), considering the specific cultivation location, with a specific emphasis on the East Java Province. To enhance the clarity, the cost figures presented in this paper have been converted into Euro currency (EUR). The 'dried biomass' in this context refers to biomass with no moisture content, unlike dried microalgae powder available in the market that has moisture content ranges from 2% to 9% (Amin et al. 2021; Zhang et al. 2022). In other words, the production cost per unit mass could potentially be lower if the moisture content of market-available microalgae powder were applied as the benchmark.
The cultivation cost estimation primarily focuses on two pivotal components: fertilizer and energy cost. A 10-liter effective working volume of a vertical tubular photobioreactor serves as the basis for the estimation. The energy cost is associated with 24-hour daily aeration and LED lighting. Notably, expenditure categories such as investment, labor, and downstream processing cost have been intentionally excluded from this analysis, considering for less developed regions, where home production or micro- or middle enterprises may handle the production process, or if it is intended for direct consumption.
Figure 8 shows the cultivation cost estimation of all treatments. More detail of each cost component, number of photobioreactors and energy required for the cultivation can be seen in Supplement 4. In comparison to the control media cost (EUR 8.24), all tested media exhibited lower cultivation costs, primarily attributed to the considerably lower price of fertilizer. While there was no significant difference in fertilizer costs among the tested samples, the lowest fertilizer cost was presented by all treatments in Ratio-A, with the standard HF ratio (A1) presented the lowest one (EUR 1.42).
The cost for energy of all treatments was significantly higher than that of the fertilizer cost, except for cultivation in BBM that was only slightly lower, which was EUR 4.39 and EUR 3.86 for fertilizer and energy cost, respectively. Energy cost was influenced by the number of photobioreactors utilized for the cultivation, as each requires its own aerator and a set of LED light strips. Another influencing factor is the duration of treatment, which was proportional to the cost incurred for lighting and aeration. While a higher biomass yield may reduce the number of bioreactors needed to produce a kilogram of biomass, an extended cultivation duration negatively impacted the total energy cost (Table 5). Comparing the energy costs of all treatments, the control media shows a moderate cost. For a more detail estimation of cultivation cost, please refer to Supplement 4. The energy cost of the control was comparable to that of most of the tested media, except for the concentration-3 of HF that shows relatively higher (> EUR 4.5).
Table 5. Factors influence the cultivation cost
Parameters
|
BBM
|
Ratio-A
|
Ratio-B
|
Ratio-C
|
A1
|
A2
|
A3
|
B1
|
B2
|
B3
|
C1
|
C2
|
C3
|
Biomass productivity (g DW.L-1.day-1)
|
0.42
|
0.47
|
0.42
|
0.31
|
0.44
|
0.49
|
0.36
|
0.48
|
0.42
|
0.33
|
Biomass density at peak (g DW. L-1)
|
6.71
|
5.65
|
7.07
|
9.00
|
5.78
|
6.83
|
9.75
|
5.77
|
6.34
|
8.96
|
Mass of biomass in 10l PBR
|
67.11
|
56.54
|
70.68
|
89.95
|
57.85
|
68.33
|
97.45
|
57.73
|
63.36
|
89.65
|
Number of PBR required to produce 1 kg DW of Biomass
|
14.90
|
17.69
|
14.15
|
11.12
|
17.29
|
14.63
|
10.26
|
17.32
|
15.78
|
11.16
|
Time to reach the peak density (day)
|
16
|
11
|
17
|
29
|
13
|
14
|
27
|
11
|
15
|
27
|