Microalgal Growth in Paper Industry E ﬄ uent: Coupling Biomass Production with Nutrients Removal

: Paper and pulp industries produce e ﬄ uents with high phosphorus concentrations, which need to be treated before their discharge in watercourses. The use of microalgae for this purpose has attracted the attention of researchers because: (i) microalgae can assimilate phosphorus (one of the main nutrients for their growth); and (ii) growing on e ﬄ uents can signiﬁcantly reduce the costs and environmental impact of microalgal biomass production. This study evaluated the growth and ability of Chlorella vulgaris to remove the phosphorus from a secondary-treated e ﬄ uent of a Portuguese paper company. Batch experiments were performed for 11 days using di ﬀ erent dilutions of the e ﬄ uent to evaluate its inhibitory e ﬀ ect on microalgae. Results showed that the non-diluted e ﬄ uent inhibited microalgal growth, indicating that this bioremediation process is possible after a previous dilution of the e ﬄ uent. Regarding phosphorus removal, promising results were achieved, especially in the experiments conducted with the most diluted e ﬄ uent: removal e ﬃ ciencies obtained in these conditions were (54 ± 1)%. Another interesting ﬁnding of this study was microalgal growth in ﬂakes’ form (mainly due to the compounds present in the e ﬄ uent and to the pH values achieved), which can be an important economic advantage for biomass recovery after the remediation step. mol m − 2 s − 1


Introduction
Paper and pulp industries require large amounts of water during their manufacturing stages. For example, the production of 1 kg of paper requires 10 to 50 L of water [1]. At the same time, large amounts of effluents (about 2000 m 3 d −1 ) are generated, presenting as main features [1,2]: (i) high chemical oxygen demand (COD, 1000-13,000 mg O2 L −1 ); (ii) high total suspended solids contents; (iii) non-biodegradable organic materials; (iv) adsorbable organic halogens (AOX): (v) color; (vi) phenolic compounds; (vii) high total phosphorus contents; and (viii) limiting nitrogen concentrations. Due to the large volumes involved and respective compositions, discharge of these effluents without any treatment can cause several environmental problems [1,3]: (i) colored effluents can affect aesthetics, PO 4 -P; (iii) 3000.15 mg L −1 of COD; and (iv) 2944 mg L −1 of biochemical oxygen demand (BOD). The experiments, aimed at evaluating both biomass production and nutrients removal efficiencies, were performed in batch mode, for 28 days, using open ponds as cultivation system (outdoor conditions). Regarding nutrients uptake, the most promising results were obtained for the 40% dilution: (i) 65% for NO 3 -N removal; (ii) 81.3% for PO 4 -P; and (iii) 75% for COD; and (iv) 82% for BOD.
The main goal of the present study was to evaluate biomass production and phosphorus removal from a secondary-treated effluent of a Portuguese paper company using the microalga Chlorella vulgaris. Different dilutions were performed to evaluate possible inhibitory effects of the effluent on microalgal growth and phosphorus uptake ability.

Microalgae Strain and Maintenance Medium
The microalga C. vulgaris (CCAP 211/11B) was obtained from the Culture Collection of Algae and Protozoa, United Kingdom. The strain was maintained on modified Organization for Economic Co-operation and Development (OECD) culture medium [14], with the following composition (mg L −1 ): 119 KNO 3

Paper Industry Effluent and Culture Conditions
Effluent from a Portuguese paper company, collected after the secondary treatment step, was characterized (Table 1) and employed as a culture medium for microalgal growth. The methodology adopted for effluent characterization was the following: (i) COD and turbidity were determined according to the Standard Methods for the Examination of Water and Wastewater [15] (through the 5220-D and 2130-B tests, respectively); (ii) total dissolved carbon (TDC), dissolved organic carbon (DOC) and inorganic carbon (DIC) were determined using an organic carbon analyzer (TOC-V CSN , Shimadzu); and (iii) chlorides, sulfates, nitrates, nitrites and phosphates were determined through ion chromatography (ICS-2100, Dionex). Due to the low concentration of nitrogen in the effluent, when compared with typical nutritional requirements of microalgae, the effluent was supplemented with NaNO 3 to achieve N:P molar ratios ranging between 6:1 and 9:1. Ratios between 5:1 and 30:1 have been considered adequate for several microalgal species [16,17]. Table 1. Physicochemical characterization of the paper industry effluent used in this study.

Parameters
Values Unit Batch experiments were performed in 1-L borosilicate glass flasks (VWR, Portugal) with a working volume of 950 mL for 11 days. The raw effluent (assay 1) and four different dilutions with freshwater (assays 2-5) were used as the culture medium for microalgal growth, with nitrogen concentrations (corresponding to the sum between nitrate-and nitrite-nitrogen) ranging between 12.7 mg N L −1 and 34.2 mg N L −1 and phosphorus concentrations (phosphate-phosphorus) ranging between 4.01 mg P L −1 and 12.3 mg P L −1 . The medium was inoculated with 250 mL of C. vulgaris inoculum to obtain an initial biomass concentration of~68 mg dw L −1 . The cultures were continuously exposed to: (i) photosynthetically active radiation between 30-40 µmol m −2 s −1 , using a 34-W white led panel; and (ii) atmospheric air filtered with 0.45-µm nylon membranes (Specanalitica, Portugal), injected at~90 L h −1 , using Trixie AP 180 air pumps (Trixie, Tarp, Germany). The experimental setup is shown in Figure 1. Two independent experiments were performed for each assay.

Parameters
Values Unit

Microalgal Growth Monitoring and Kinetic Growth Parameters
Operational parameters, such as pH and temperature, were daily monitored using a SympHony SB90M5 pH-meter (VWR, Portugal). Microalgal growth was also daily assessed through optical density measurements at 680 nm (OD 680 ) using a UV-6300 PC spectrophotometer (VWR, United States). To eliminate the interference of the effluent color on OD 680 measurements, the cells were separated from the culture medium by centrifugation (at 4000 rpm, for 10 min), the supernatant was discarded, and the cells were resuspended in an equal volume of distilled water, as described by Hodaifa, et al. [18]. This procedure was repeated twice. The relationship between OD 680 and biomass concentration (X, mg dw L −1 ) for C. vulgaris was previously established by linear regression, according to Equation (1): Biomass concentrations were used to determine the kinetic growth parameters: (i) specific growth rates (µ, d −1 ); and (ii) maximum and average biomass productivities (P max and P aver , mg dw L −1 d −1 ). The specific growth rates were determined from the first-order kinetic model, according to Equation (2): where X 1 and X 0 correspond, respectively, to biomass concentration (mg dw L −1 ) in the end (t 1 , d) and beginning (t 0 , d) of the exponential growth phase. Biomass productivities (P, mg dw L −1 d −1 ) were calculated for each pair of consecutive points, through Equation (3): where X z represents the biomass concentration (mg dw L −1 ) at time t z (d) and X z+1 corresponds to the biomass concentration (mg dw L −1 ) at time t z+1 (d). The maximum productivity was determined from the maximum value obtained from Equation (3). On the other hand, average biomass productivities were determined according to Equation (4): where X f and X i correspond, respectively, to biomass concentration (mg dw L −1 ) in the end (t f , d) and beginning (t i , d) of the cultivation period.

Nutrients Removal
Nutrients removal was evaluated in terms of nitrogen and (N) phosphorus (P) present in the culture medium/effluent. Nitrogen was assessed in the forms of nitrate and nitrite ions, whereas phosphorus was monitored through the presence of phosphate ions. From each assay, 5 mL of the microalgal suspension were periodically collected (days 0, 1, 2, 4, 7, 9, and 11). These samples were centrifuged at 4000 rpm, for 10 min, and the supernatants were filtered through 0.45-µm nylon membranes (Specanalitica, Portugal). Nitrate, nitrite and phosphate concentrations were determined in an ion chromatograph (ICS-2100, Dionex) equipped with an anion analytical column (4x 250 mm, AS11-HC) and a self-regeneration suppressor (4 mm, AERS 500). The values obtained in the first and last day of culturing were used to calculate the following removal parameters: (i) removal efficiencies (%R, %); (ii) average removal rates (RR, mg L −1 d −1 ); and (iii) mass removal (R, mg L −1 ), as shown in Equations (5), (6) and (7), respectively: where S f and S i correspond to the nitrogen (nitrate + nitrite) or phosphorus (phosphate) concentration (mg L −1 ) in the end (t f , d) and beginning (t i , d) of the cultivation period, respectively.

Statistical Analysis
Each parameter shown in the present paper was expressed as the mean and standard deviation. The Tukey statistical test was used to investigate if the differences between the different effluent concentrations studied could be considered significant. These statistical tests were performed using Statistica 8.0 (StatSoft Inc., USA) and were carried out at a significance level (p) of 0.05.

Microalgal Growth
The C. vulgaris growth curves in raw and diluted paper industry effluent are shown in Figure 2. These results evidence the inexistence of an adaptation phase for all assays and an exponential growth phase that lasted approximately four days. In addition, no cell decay was observed during the 11-day batch culture, indicating that the experiments could be extended for a longer period. The increase of biomass concentration during the cultivation period, as well as the lack of an adaptation phase, shows that C. vulgaris was able to grow in this effluent. However, biomass concentrations achieved in non-diluted effluent (assay 1) were statistically lower (p < 0.05) than those achieved in more diluted effluents from assays 3-5.
Appl. Sci. 2020, 10, x FOR PEER REVIEW 7 of 17 cultures in different ways [19][20][21]: (i) the effluent color may act as a barrier to light penetration, thus limiting microalgal access to light and photosynthetic activity; and (ii) paper industry effluents are characterized by the presence of lignin, humic acids, furans, and dioxins and by high levels of aluminum and manganese, which exhibit toxic effects on microalgae.  cultures in different ways [19][20][21]: (i) the effluent color may act as a barrier to light penetration, thus limiting microalgal access to light and photosynthetic activity; and (ii) paper industry effluents are characterized by the presence of lignin, humic acids, furans, and dioxins and by high levels of aluminum and manganese, which exhibit toxic effects on microalgae. cultures in different ways [19][20][21]: (i) the effluent color may act as a barrier to light penetration, thus limiting microalgal access to light and photosynthetic activity; and (ii) paper industry effluents are characterized by the presence of lignin, humic acids, furans, and dioxins and by high levels of aluminum and manganese, which exhibit toxic effects on microalgae. cultures in different ways [19][20][21]: (i) the effluent color may act as a barrier to light penetration, thus limiting microalgal access to light and photosynthetic activity; and (ii) paper industry effluents are characterized by the presence of lignin, humic acids, furans, and dioxins and by high levels of aluminum and manganese, which exhibit toxic effects on microalgae. cultures in different ways [19][20][21]: (i) the effluent color may act as a barrier to light penetration, thus limiting microalgal access to light and photosynthetic activity; and (ii) paper industry effluents are characterized by the presence of lignin, humic acids, furans, and dioxins and by high levels of aluminum and manganese, which exhibit toxic effects on microalgae. cultures in different ways [19][20][21]: (i) the effluent color may act as a barrier to light penetration, thus limiting microalgal access to light and photosynthetic activity; and (ii) paper industry effluents are characterized by the presence of lignin, humic acids, furans, and dioxins and by high levels of aluminum and manganese, which exhibit toxic effects on microalgae. To complement the analysis from growth curves, microalgal growth parameters, such as specific growth rate, maximum biomass concentration, and maximum and average biomass productivities, were determined and presented in Table 2. From these data, it is possible to see a general increase in growth parameters from assay 1 to assay 5, i.e., from the non-diluted effluent to the more diluted one. Regarding specific growth rates, values ranged from (0.093 ± 0.007) d −1 to (0.16 ± 0.02) d −1 in assays 1 and 5, respectively. The highest values of maximum biomass concentrations were also obtained in more diluted effluents from assays 4 and 5: (249 ± 14) mg dw L −1 and (231 ± 31) mg dw L −1 , respectively. Similar behavior was observed for both maximum and average biomass productivities. Maximum biomass productivities/average biomass productivities obtained in assays 4 and 5 were (30 ± 3)/(16 ± 1) mg dw L −1 d −1 and (30 ± 6)/(15 ± 3) mg dw L −1 d −1 , respectively. In opposition, maximum and average biomass productivities obtained in assay 1 were (9.8 ± 0.2) mg dw L −1 d −1 and (6.2 ± 0.1) mg dw L −1 d −1 , respectively. Table 2. Specific growth rates (µ, in d −1 ), maximum biomass concentrations (X max , in mg dw L −1 ), and maximum and average biomass productivities (P max and P aver , in mg dw L −1 d −1 ) determined for C. vulgaris grown in raw and diluted secondary-treated paper industry effluent.
Values are presented as the mean ± standard deviation obtained from two independent experiments. Within the same column, mean values sharing at least one common letter (in superscript) are not statistically different (p > 0.05).
In contrast to what was observed by Gentili [13], the increment in nitrogen and phosphorus concentration did not contribute to an increase in kinetic growth parameters. Accordingly, these results may indicate inhibitory effects of the effluent on microalgae, which can influence microalgal cultures in different ways [19][20][21]: (i) the effluent color may act as a barrier to light penetration, thus limiting microalgal access to light and photosynthetic activity; and (ii) paper industry effluents are characterized by the presence of lignin, humic acids, furans, and dioxins and by high levels of aluminum and manganese, which exhibit toxic effects on microalgae.
Most studies regarding the bioremediation of paper industry effluents with microalgae focus on the removal of contaminants and only a few report biomass production yields. Polishchuk, et al. [20] reported that the maximum specific growth rate obtained for Nannochloropsis oculata grown in effluents resulting from pulp and paper industry was 0.405 d −1 . Tao, et al. [19] revealed that maximum biomass concentrations achieved by Scenedesmus acuminatus and C. vulgaris grown in paper industry effluents were 291 mg L −1 and 822 mg L −1 , respectively. Considering the values referred in the literature, microalgal growth parameters obtained in this study were significantly lower, which can be attributed to the inhibitory effects promoted by the effluent used (in assays 1-3) and to the low concentration of some essential nutrients (in more diluted effluents of assays 4 and 5). Another explanation for the low biomass concentrations and productivities achieved may be related to the phenomenon of flakes formation observed within the cultivation period (autoflocculation). Cells' agglomeration can affect the accurate measurement of OD 680 and, on the other hand, it can reduce light absorption efficiency by cells incorporated within flakes, thus resulting in lower photosynthetic activity. In this study, this phenomenon occurred due to the increase of culture pH (from 7.8 to 8.6) or due to the presence of certain compounds in the effluent, which can induce a change in the surface charge of the cells and affect suspensions' stability [22]. Despite the low microalgal growth rates, the flakes formation enables a cost-effective biomass removal after effluent remediation. The density similar to water and small size of microalgal cells difficult the harvesting process and make this step one of the most expensive within microalgal biomass production processes [22,23]. However, when cells agglomerate, an increase in their density and size is observed, contributing to higher settling rates and allowing biomass recovery using the least expensive harvesting method: sedimentation.

Nutrients Removal
In this study, nitrogen (in the forms of nitrate and nitrite) and phosphorus (in the form of phosphate) concentrations were monitored within the cultivation time to evaluate the potential of C. vulgaris to uptake these nutrients from a paper industry effluent with different concentrations of both nitrogen and phosphorus. Figure 3 shows the variation of nitrogen and phosphorus concentration in each assay. Regarding nitrogen removal ( Figure 3A), this element was readily assimilated by C. vulgaris in the diluted effluents (assays 2-5). In the raw effluent (corresponding to assay 1), a two-day delay was observed in nitrogen assimilation, which may be related to the adaptation of the microalga to these conditions. Regarding the assimilation patterns observed in assays 2-5, these were approximately linear for assays 2-4, with nitrogen concentration decreasing gradually during the cultivation time. On the other hand, in assay 5, corresponding to the more diluted effluent experiments, nitrogen concentration decreased until the seventh day of culturing and then it was maintained approximately constant. This behavior may be attributed to a decrease in photosynthetic activity, as nitrogen concentration decreased, and explains the lower biomass concentrations achieved in assay 5 when compared to the one obtained in assay 4 (according to Table 2, (231 ± 31) mg dw L −1 and (249 ± 14) mg dw L −1 , respectively). Also, at the end of the cultivation time, nitrogen concentration remaining in cultures corresponding to assays 4 and 5 was approximately the same ((2.81 ± 0.05) mg N L −1 and (2.6 ± 0.2) mg N L −1 , respectively), indicating a limitation of this nutrient in the last days of assay 5. As for nitrogen concentration, phosphorus concentration also decreased within the cultivation time ( Figure 3B), but in a lesser extent, which is related with microalgal nutritional requirements, as given by its typical elemental biochemical composition: CO 0.48 H 1.83 N 0.11 P 0.01 [24]. The reduction observed in nitrogen and phosphorus concentration in the studied effluent (raw or diluted) shows that C. vulgaris can promote an efficient uptake of both nutrients. However, except for nitrogen concentration in assay 5, total depletion of these nutrients did not occur after the 11 days of culturing, reiterating what was stated in relation to cell growth, that the cultures could be extended for an increased period to further improve nutrients removal efficiencies. Another similarity with the microalgal growth parameters already described is the higher variations in nitrogen and phosphorus concentrations observed in the experiments where the effluent was previously diluted (assays 2-5), which indicate that these conditions were more favorable for C. vulgaris photosynthetic activity.
Nitrogen and phosphorus removal parameters are presented in Figures 4 and 5, respectively. As with microalgal growth parameters, a general increase in nutrients removal efficiencies was observed from assay 1 to 5, with values ranging from (24 ± 10)% to (80 ± 4)% for nitrogen ( Figure 4A) and from (13.0 ± 0.9)% to (54 ± 1)% for phosphorus ( Figure 5A). However, Figure 4A shows that there was no statistical difference (p > 0.05) in nitrogen removal efficiency between assays 4 and 5, which can be explained by the low concentration achieved in the assay 5 (the one corresponding to the most diluted effluent) that might have been limiting for microalgal growth. In fact, according to Table 2, maximum biomass concentration achieved in assay 4 was higher than that in assay 5, indicating that the highest dilution applied in this study may have contributed to nitrogen limitation to C. vulgaris, with effects on their growth and nutrients removal parameters. Regarding nitrogen removal rates ( Figure 4B) and mass removal ( Figure 4C), the highest values were determined in assays 3 and 4 and no statistical differences were observed (p > 0.05): (i) average removal rates were (1.31 ± 0.07) mg N L −1 d −1 and (1.26 ± 0.08) mg N L −1 d −1 , respectively; and (ii) mass removal values were (14.4 ± 0.8) mg N L −1 and (13.9 ± 0.9) mg N L −1 , respectively. These results are in accordance with maximum biomass concentration achieved and indicate higher photosynthetic activity of C. vulgaris in these intermediate conditions. A different behavior was observed for phosphorus. In this case, average removal rates ( Figure 5B) and mass removal values ( Figure 5C) determined for assays 1 to 4 were not statistically different (p > 0.05), but values determined for assay 5 were statistically higher (p < 0.05), reaching an average removal rate of (0.20 ± 0.01) mg P L −1 d −1 and a mass removal of (2.2 ± 0.1) mg P L −1 .
nutrients. However, except for nitrogen concentration in assay 5, total depletion of these nutrients did not occur after the 11 days of culturing, reiterating what was stated in relation to cell growth, that the cultures could be extended for an increased period to further improve nutrients removal efficiencies. Another similarity with the microalgal growth parameters already described is the higher variations in nitrogen and phosphorus concentrations observed in the experiments where the effluent was previously diluted (assays 2-5), which indicate that these conditions were more favorable for C. vulgaris photosynthetic activity.  did not occur after the 11 days of culturing, reiterating what was stated in relation to cell growth, that the cultures could be extended for an increased period to further improve nutrients removal efficiencies. Another similarity with the microalgal growth parameters already described is the higher variations in nitrogen and phosphorus concentrations observed in the experiments where the effluent was previously diluted (assays [2][3][4][5], which indicate that these conditions were more favorable for C. vulgaris photosynthetic activity. did not occur after the 11 days of culturing, reiterating what was stated in relation to cell growth, that the cultures could be extended for an increased period to further improve nutrients removal efficiencies. Another similarity with the microalgal growth parameters already described is the higher variations in nitrogen and phosphorus concentrations observed in the experiments where the effluent was previously diluted (assays [2][3][4][5], which indicate that these conditions were more favorable for C. vulgaris photosynthetic activity. did not occur after the 11 days of culturing, reiterating what was stated in relation to cell growth, that the cultures could be extended for an increased period to further improve nutrients removal efficiencies. Another similarity with the microalgal growth parameters already described is the higher variations in nitrogen and phosphorus concentrations observed in the experiments where the effluent was previously diluted (assays 2-5), which indicate that these conditions were more favorable for C. vulgaris photosynthetic activity. did not occur after the 11 days of culturing, reiterating what was stated in relation to cell growth, that the cultures could be extended for an increased period to further improve nutrients removal efficiencies. Another similarity with the microalgal growth parameters already described is the higher variations in nitrogen and phosphorus concentrations observed in the experiments where the effluent was previously diluted (assays 2-5), which indicate that these conditions were more favorable for C. vulgaris photosynthetic activity. did not occur after the 11 days of culturing, reiterating what was stated in relation to cell growth, that the cultures could be extended for an increased period to further improve nutrients removal efficiencies. Another similarity with the microalgal growth parameters already described is the higher variations in nitrogen and phosphorus concentrations observed in the experiments where the effluent was previously diluted (assays 2-5), which indicate that these conditions were more favorable for C. vulgaris photosynthetic activity.

Time (d)
). Error bars correspond to the standard deviation of the mean obtained from two independent experiments.
Appl. Sci. 2020, 10, x FOR PEER REVIEW 10 of 17 efficiencies was observed from assay 1 to 5, with values ranging from (24 ± 10)% to (80 ± 4)% for nitrogen ( Figure 4A) and from (13.0 ± 0.9)% to (54 ± 1)% for phosphorus ( Figure 5A). However, Figure   4A shows that there was no statistical difference (p>0.05) in nitrogen removal efficiency between assays 4 and 5, which can be explained by the low concentration achieved in the assay 5 (the one corresponding to the most diluted effluent) that might have been limiting for microalgal growth. In fact, according to Table 2, maximum biomass concentration achieved in assay 4 was higher than that in assay 5, indicating that the highest dilution applied in this study may have contributed to nitrogen limitation to C. vulgaris, with effects on their growth and nutrients removal parameters. Regarding nitrogen removal rates ( Figure 4B) and mass removal ( Figure 4C   Nutrients removal from paper industry effluents has already been reported in the literature. Table 3 highlights nitrogen and phosphorus removal efficiencies and removal rates obtained in these studies. According to these data, removal efficiencies reported by Tao, et al. [19] and Gentili [2] are significantly higher than those obtained in this study, whereas values reported by Usha, et al. [1] were closer to those obtained in the present study, especially in assays 3-5. The lower removal efficiencies obtained in this study when compared with those reported by Tao, et al. [19], may be associated with the higher N:P molar ratio used in the reference study, which was~66:1. On the other hand, the higher removal efficiencies reported by Gentili [13] may be associated with the use of other effluents to achieve the dual role of providing the required nutrients for microalgal growth while contributing to a reduction in the toxicity of the paper industry effluent. Another explanation for the increased efficiencies obtained in these studies is the nitrogen source used. As in the present study, Usha, et al. [1] cultivated microalgae in an effluent with nitrate-nitrogen as the main nitrogen source. On the other hand, Tao, et al. [19] tested an effluent with ammonium as the main nitrogen source (digestate obtained from the treatment of a pulp and paper industry effluent) and Gentili [13] evaluated this treatment with both nitrogen forms present. According to several studies, although nitrate-nitrogen is the most thermodynamically stable form (and the most commonly found in aquatic environments), ammonia is directly assimilated and converted into proteins by microalgae, while nitrate must be reduced to nitrite and then to ammonia before being assimilated by microalgal cells [25]. However, for an adequate comparison of nutrients removal performance, it is important to determine the average removal rate, as this parameter takes into account initial nutrients concentrations and cultivation/treatment time. Comparing average removal rates obtained in the present study and in the reference studies, values in the same order of magnitude were obtained, except in what concerns ammonium-nitrogen removal in the studies performed by Tao, et al. [19] and Gentili [2]. In these cases, the higher removal rates obtained may be associated with the higher ability of microalgae to assimilate ammonium-nitrogen than nitrate-nitrogen. Considering values of average RR, it is possible to conclude that promising results were obtained in this study. Moreover, differences found in experimental conditions used in this study and in the studies reported in the literature demonstrate that these results can be significantly enhanced. Besides increasing N:P molar ratio and providing an ammonium-nitrogen source, the increase of light supply should also be considered, as values reported in the literature correspond to cultures grown under light intensities of 130-800 µmol m −2 s −1 , whereas results reported in the present study were obtained with light intensities of 30-40 µmol m −2 s −1 .
In summary, the results obtained in this study for both nitrogen and phosphorus removal evidence that the remediation of paper industry effluents using microalgae is possible, provided that it is properly diluted to avoid inhibitory effects related to the presence of strong color or high concentrations of toxic compounds, typically associated with effluents resulting from this industrial sector [19,20]. Considering the results obtained for nitrogen removal, the dilution of the effluent to the concentrations present in assays 3 and 4 is the most adequate. In these conditions, nitrogen concentrations were significantly reduced, reaching (7.1 ± 0.7) mg N L −1 and (2.81 ± 0.05) mg N L −1 , respectively (which corresponds to the highest average removal rates: (1.31 ± 0.07) mg N L −1 d −1 and (1.26 ± 0.09) mg N L −1 d −1 , respectively). Regarding phosphorus removal, the highest removal rate was obtained for the conditions tested in assay 5: (0.20 ± 0.01) mg P L −1 d −1 . Table 3. Comparison between nutrients removal efficiencies (%R, in %) and average removal rates (RR, in mg L −1 d −1 ) obtained in this study and other studies reporting microalgal growth in effluents resulting from pulp and paper industries. Despite the promising nitrogen and phosphorus removal rates, the results obtained in this study demonstrated that the cultures were limited by nitrogen, as nitrogen and phosphorus were assimilated by C. vulgaris at a N:P molar ratio ranging from 10:1 to 24:1. Considering these results and the N:P molar ratios used in this study (between 6:1 and 9:1), nutrients uptake could be enhanced by increasing nitrogen supply. Another alternative to achieve an adequate N:P molar ratio and reduce the toxicity of this effluent would be to dilute it with other effluents, as proposed in other studies [13]. Finally, the remediation process could be further improved by modulating microalgal cultivation conditions. According to Gonçalves, et al. [25], light conditions, temperature, and pH are also important parameters that can influence microalgal growth and, hence, the efficiency of the bioremediation process. Thus, from the prospecting of this work, other studies evaluating these parameters should be carried out to further improve nitrogen and phosphorus uptake from paper industry effluents.

Conclusions
This study showed the feasibility of using C. vulgaris for the bioremediation of a paper industry effluent fortified with a nitrogen source, targeting phosphorus removal. C. vulgaris was able to grow in all studied effluent conditions (in non-diluted and diluted ones). However, it was possible to conclude that growing on non-diluted effluent resulted in lower biomass productivities, which was also reflected in nitrogen and phosphorus removal efficiencies. From microalgal growth and nitrogen removal points of view, the effluent dilutions used in assays 3 and 4 (intermediate dilutions) seem to be the most adequate, as microalgal growth was not inhibited in these conditions and nitrogen mass removal was quite satisfactory, achieving final concentrations of (7.1 ± 0.7) mg N L −1 and (2.81 ± 0.05) mg N L −1 , respectively. Regarding phosphorus removal, concentrations achieved in the last day of culturing in assays 3 and 4 were higher ((4.63 ± 0.04) mg P L −1 and (2.940 ± 0.005) mg P L −1 , respectively) than the one obtained in assay 5 ((1.85 ± 0.02) mg P L −1 ). However, the results obtained in assay 5 suggest a growth limitation, mainly related to nitrogen concentration. Accordingly, the obtained results indicate that these values can be further improved by studying different N:P molar ratios, different microalgal cultivation conditions, dilution with other effluents, among others. Improving the remediation performance can significantly contribute to the development of an effective microalgae-based remediation process of pulp and paper industry effluents.