Nutrient Removal and Membrane Performance of an Algae Membrane Photobioreactor in Urban Wastewater Regeneration

The increase in industry and population, together with the need for wastewater reuse, makes it necessary to implement new technologies in the circular economy framework. The aim of this research was to evaluate the quality of the effluent of an algae membrane photobioreactor for the treatment of the effluent of an urban wastewater treatment plant, to characterise the ultrafiltration membranes, to study the effectiveness of a proposed cleaning protocol, and to analyse the performance of the photobioreactor. The photobioreactor operated under two days of hydraulic retention times feed with the effluent from the Los Vados wastewater treatment plant (WWTP) (Granada, Spain). The microalgae community in the photobioreactor grew according to the pseudo-second-order model. The effluent obtained could be reused for different uses of diverse quality with the removal of total nitrogen and phosphorus of 56.3% and 64.27%, respectively. The fouling of the polyvinylidene difluoride ultrafiltration membrane after 80 days of operation was slight, increasing the total membrane resistance by approximately 22%. Moreover, the higher temperature of the medium was, the lower intrinsic resistance of the membrane. A total of 100% recovery of the membrane was obtained in the two-phase cleaning protocol, with 42% and 58%, respectively.


Introduction
Microalgae are currently some of the most promising renewable raw materials to produce different subproducts [1][2][3]. From both an economic and environmental point of view, microalgae are attractive because their cultivation and processing could involve wastewater treatment and carbon dioxide capture from combustion gases [4][5][6]. Currently, microalgae-based systems are some of the most promising technologies for tertiary treatment in wastewater treatment plants (WWTPs). This is due to their associated low costs and environmental advantages [7,8].
Nowadays, microalgae cultivation has a high cost in terms of the nutrients needed. Therefore, the use of wastewater may be a viable solution because of its nutrient content, coupled with the ability of these micro-organisms to tolerate very wide conditions of temperature, pH, or salinity [9,10]. Wastewater provides the nutrients necessary for the growth of the microalgae, in turn reducing the cost of production and greenhouse gas emissions because the microalgae will absorb some of the carbon generated in the wastewater treatment that would be released as carbon dioxide [6,9,11]. This achieves the requirements of wastewater discharge standards, which specify very low nutrient concentrations to avoid eutrophication [12].
Wastewater treatment with microalgae reduces the costs of microalgae biomass production by 40% of the total cost, generating a lower energy demand for nutrient removal and, therefore, presenting lower costs compared to conventional wastewater treatment [13][14][15]. On one hand, microalgae can use inorganic phosphorus and nitrogen for their growth, To automate and control the correct operation of the photobioreactor, a vacuum and pressure gauge (Potermic, Barcelona, Spain) were installed to control the pressure during the filtering process and level probes to keep the photobioreactor at an approximately constant volume. In addition, the photobioreactors have a continuous stirring system of approximately 8 Hz (MSL 63 A, Cime Motors, Girona, Spain) and continuous air injection by means of a circular diffuser at the base of the photobioreactor. In addition to natural lighting, the photobioreactor was illuminated by fluorescent lights (L36 W/868, OSRAM DS, Munich, Germany) placed around the plant, providing an average irradiance in the photobioreactor of 13.24 μmol/m 2 ·s. A 12/12-h light/dark cycle was established by means of a timer. The fibres of the ultrafiltration membrane were taken and divided into four smaller filtration modules, adapted to the photobioreactor used in this research. To this end, the fibres that make up the new module were sealed with silicone to the PVC parts that allow the membrane to be connected to the filtration system. The length of each hollow fibre was 500 mm, and the total membrane surface area of this new membrane was 0.2 m 2 .
The permeate was extracted through the membrane by a suction-backwashing peristaltic pump (323 U, Watson-Marlow Pumps Group, Marlow, UK). The permeate extraction and backwashing periods were 9 and 1 min, respectively.
To automate and control the correct operation of the photobioreactor, a vacuum and pressure gauge (Potermic, Barcelona, Spain) were installed to control the pressure during the filtering process and level probes to keep the photobioreactor at an approximately constant volume. In addition, the photobioreactors have a continuous stirring system of approximately 8 Hz (MSL 63 A, Cime Motors, Girona, Spain) and continuous air injection by means of a circular diffuser at the base of the photobioreactor. In addition to natural lighting, the photobioreactor was illuminated by fluorescent lights (L36 W/868, OSRAM DS, Munich, Germany) placed around the plant, providing an average irradiance in the photobioreactor of 13.24 µmol/m 2 ·s. A 12/12-h light/dark cycle was established by means of a timer.
To determine the average irradiance (light intensity) in the photobioreactor, an HD 2303.2 Delta OHM light meter (Arava, Spain) was used with the DLP471 PAR probe (Arava, Spain) for photon measurement. Measurements were taken on and around the photobioreactor surface, both perpendicular and parallel to the photobioreactor.
The photobioreactor was continuously fed with treated urban wastewater from real wastewater taken from WWTP of Los Vados (Granada, Spain), and naturally occurring microalgae species were cultured. The HRT set in the photobioreactor was two days, where the transmembrane operating pressure (TMP) was approximately 0.01 bar. The experimental cycle lasted 80 days, coinciding with the SRT.

Experimental Procedure
Samples were taken on alternate days in the influent, photobioreactor, and effluent. Once the samples were collected, the pH was measured with a Crison pH 25 meter (Hach, Barcelona, Spain), and the temperature and conductivity were measured using the multifunctional meter PCE-PHD 1 (PCE Ibérica SL, Alicante, Spain). For the analysis of Total Suspended Solids (TSS) present in the sample, 0.45 µm Millipore glass fibre filters were used following standard methods [33]. The filtrate was used for the analysis and quantification of nutrients in the samples, determined by ion chromatography using a conductivity detector (Metrohm, Metrohm AG, Herisau, Switzerland). Finally, the Heλios spectrophotometer (Thermo Spectronic, Madrid, Spain) was used to measure turbidity at a wavelength of 650 nm and the optical density of the culture at 680 nm, following the same procedure as other authors [8,34]. Table 1 shows the six kinetic models that were evaluated to adjust the kinetics of microalgae biomass growth in the photobioreactor, measured as TSS throughout the evolution of the concentration over the time (C t ): Table 1. Kinetic models.

Kinetic Model Equation
Pseudo-first-order C 0 (g·L −1 ) represents the initial biomass concentration, C t (g·L −1 ) is the biomass concentration at any time, and C e (g·L −1 ) is the biomass equilibrium concentration. The experimental time is expressed as t (days) and the model constant as k (days −1 ). The parameter h relates the biomass concentration at equilibrium and the kinetic model constant as follows: h = k·C 2 e (g 2 ·L −2 ·days −1 ). The data were fitted to the models described above, minimising the sum of squared errors between the empirical and modelled data using the Microsoft Excel Solver application.

Membrane Filtration Tests
The membrane was characterised prior to use to determine the intrinsic resistance of the clean membrane. For this purpose, a test was carried out with clean water. Subsequently, the behaviour of the clean membrane was studied with the content of the photobioreactor (microalgae culture in treated urban wastewater).
At the end of the study cycle, the membrane was again characterised with water and the content of the photobioreactor. In this way, the total resistance of the membrane after fouling was obtained.
Finally, the recovery of the membrane permeability was studied by means of a twostage cleaning protocol. In the first stage of the cleaning protocol, the membrane was inorganically cleaned with citric acid (800 ppm) for 5 h. The second cleaning step consisted of an organic cleaning with sodium hypochlorite (800 ppm) for 5 h. After each of the cleaning stages, the membrane was characterised in water to evaluate the cleaning efficiency of each stage. Each of the tests was carried out at three different temperatures (20,25, and 30 • C). To carry them out, the volumes obtained after 30 s of filtration at different values of TMP were recorded. In this way, filtration fluxes were obtained for each of the TMP sets.

Membrane Characterisation
The membrane permeability (K) was calculated from the following equation as described above [35], where J is the permeate flux (LHM, litres per square metre per hour filtered) and TMP is the measured transmembrane pressure (bar): R M represents the intrinsic resistance of the membrane, R CP is the resistance that occurs in the first moments of filtration, which is reversible, and R F corresponds to the resistance that occurs at longer times, which corresponds to fouling. R M and R F represent the total value of the intrinsic membrane resistance and the total value of the fouling resistance, respectively, by integrating the viscosity constant for each temperature.
This permeability constant encompasses the influence of temperature and the intrinsic resistance of the membrane. On the other hand, the intrinsic membrane resistance (R M ) was calculated as the inverse of the permeability, which had been obtained with the membrane cleaned in water: The fouling resistance (R F ) of the membrane was calculated as described in Equation (3), where R i represents the total membrane resistance in each of the cases studied.

Recovery of Membrane Characteristics
The efficiency of the membrane cleaning was calculated as shown in Equation (4), where R 0 represents the total membrane resistance with a soiled membrane, R i represents the total membrane resistance in each of the cases studied, R i−1 represents the total membrane resistance of the case before the one under study, and R M is the intrinsic membrane resistance:

Analysis of the Operational TMP
The TMP was determined by setting the flow rate necessary for the HRT in the photobioreactor of the treated urban wastewater to be two days. It was verified that the operating TMP of the membrane was in the pressure-controlled working region.
To evaluate the flux at different TMPs using the contents of the photobioreactor, tests were carried out with both clean and dirty membranes after the end of the experimental cycle in the photobioreactor. For this purpose, the equation of the series resistance model described below was used: The parameter b is a constant which relates the behaviour of the TMP to the matter transfer in the filtration process.

Photobioreactor Operation
After analyses of the samples were taken, the average operating conditions of the reactor during the experimental period were 7.40 ± 0.88 for pH, 15.54 ± 0.89 • C for temperature, and 954.30 ± 75.82 µS/cm for conductivity. transfer in the filtration process.

Photobioreactor Operation
After analyses of the samples were taken, the average operating conditions of the reactor during the experimental period were 7.40 ± 0.88 for pH, 15.54 ± 0.89 °C for temperature, and 954.30 ± 75.82 μS/cm for conductivity. Figure 2a shows the evolution of the TSS in the photobioreactor during the experimental time. The concentration of algal biomass increases as time passes until the stationary phase is reached. In this case the results obtained during the growth phase have been represented up to their maximum point, where the stationary phase, not represented, would begin. In the same way, the turbidity results obtained in the photobioreactor are shown in Figure 2b and expressed in NTU (Nephelometric turbidity unit), which are related to the TSS concentration in the photobioreactor, increasing as the concentration of microalgae biomass increases. Figure 3 shows the evolution of the optical density of the microalgae culture. These values were obtained by taking the difference between the absorbance measurement at 680 nm of the microalgae photobioreactor culture and the measurement obtained for the microalgae feed (treated urban wastewater). In the same way, the turbidity results obtained in the photobioreactor are shown in Figure 2b and expressed in NTU (Nephelometric turbidity unit), which are related to the TSS concentration in the photobioreactor, increasing as the concentration of microalgae biomass increases. Figure 3 shows the evolution of the optical density of the microalgae culture. These values were obtained by taking the difference between the absorbance measurement at 680 nm of the microalgae photobioreactor culture and the measurement obtained for the microalgae feed (treated urban wastewater).  The same growth trend was observed in the three parameters studied, with a diffe ence-in the study of the optical density, the interferences that could be caused by the TS contained in the treated urban wastewater, bacteria, and other microorganisms that influ ence the correct evaluation of the growth of the algal biomass were eliminated. Figure  The same growth trend was observed in the three parameters studied, with a differencein the study of the optical density, the interferences that could be caused by the TSS contained in the treated urban wastewater, bacteria, and other microorganisms that influence the correct evaluation of the growth of the algal biomass were eliminated. Figure 4 shows the microalgae present in the photobioreactor, which grew naturally in the treated urban wastewater. The same growth trend was observed in the three parameters studied, with a difference-in the study of the optical density, the interferences that could be caused by the TSS contained in the treated urban wastewater, bacteria, and other microorganisms that influence the correct evaluation of the growth of the algal biomass were eliminated. Figure 4 shows the microalgae present in the photobioreactor, which grew naturally in the treated urban wastewater.  The kinetic constants and correlation rates obtained for each of the kinetic models evaluated for the concentration of TSS and for the optical density obtained in the photobioreactor, which has been evaluated at a wavelength of 680 nm, are presented in Table 2:   The kinetic constants and correlation rates obtained for each of the kinetic models evaluated for the concentration of TSS and for the optical density obtained in the photobioreactor, which has been evaluated at a wavelength of 680 nm, are presented in Table 2: Table 2. Kinetic constants and correlation rates for each evaluated kinetic model.

Kinetic Model Total Suspended Solids (TSS) Optical Density (OD) Model Constants Correlation Rate R 2 Model Constants Correlation Rate R 2
Pseudo-first order Pseudo-second order C e = 1.075 g/L k = −1. Comparing the results obtained, in the case of the TSS, the kinetic model with the highest correlation index is the pseudo-second-order model. However, in the case of optical density, the best correlation rate is obtained for the zero-order model. Figures 5 and 6 show graphical representations of the adjustments that were made for the TSS and optical density (OD) data.
The kinetic constants obtained for the fit of the optical density data show, in absolute value, had higher values compared to those obtained for the TSS, except for the pseudosecond-order model. This indicates that the overall growth observed is slower than the growth of micro-organisms detected at 680 nm, including microalgae.
The best correlation rate, and therefore the best fit, was obtained for the zero-order model when evaluating the optical density values. These results are considered satisfactory because the optical density measurement eliminated interferences caused by bacteria or other microorganisms in the TSS measurement. The kinetic growth of microalgae can be characterised by a phase in which the specific growth rate exists until it reaches a maximum value, where it is maintained for a specific time. This kinetic growth model corresponds to the one described by Monod, one of the most widely used [7,36]. No substrate-dependent data are available to fit this model, as the photobioreactor is assumed to operate under perfect mixing conditions. Therefore, the nutrient concentrations (the substrate) are approximately constant throughout the experiment. For this reason, the fit of the experimental data to the Monod growth kinetic model could not be determined.
No studies evaluating microalgae growth using polynomial kinetic models have been found in the literature. However, they have been evaluated to represent the growth of microalgae attached to support materials [37][38][39]. Comparing the results obtained, in the case of the TSS, the kinetic model with the highest correlation index is the pseudo-second-order model. However, in the case of optical density, the best correlation rate is obtained for the zero-order model. Figure 5 and Figure 6 show graphical representations of the adjustments that were made for the TSS and optical density (OD) data.
(e) (f)   The kinetic constants obtained for the fit of the optical density data show, in absolute value, had higher values compared to those obtained for the TSS, except for the pseudosecond-order model. This indicates that the overall growth observed is slower than the growth of micro-organisms detected at 680 nm, including microalgae.
The best correlation rate, and therefore the best fit, was obtained for the zero-order model when evaluating the optical density values. These results are considered satisfactory because the optical density measurement eliminated interferences caused by bacteria or other microorganisms in the TSS measurement.
The kinetic growth of microalgae can be characterised by a phase in which the specific growth rate exists until it reaches a maximum value, where it is maintained for a specific time. This kinetic growth model corresponds to the one described by Monod, one of the most widely used [7,36]. No substrate-dependent data are available to fit this model, as the photobioreactor is assumed to operate under perfect mixing conditions. Therefore,

Membrane Characterisation
As described above, before starting the study in the photobioreactor, the resistance of the cleaned membrane at different temperatures was characterised.
As expected, as the temperature increases, the membrane resistance decreases, taking values between 6.21 × 10 −3 and 5.8 × 10 −3 bar·m 2 ·h/L ( Table 3). As previously studied by other authors, increasing the temperature increases the permeability of PVDF ultrafiltration membranes, or in other words, the membrane resistance decreases, as shown in our results [35,40]. Once the photobioreactor study was completed after 80 days of filtration, the membrane was characterised to assess its fouling. The data obtained are shown in Table 4, with the total membrane resistances (R i ) and the fouling resistances (R F ) for each of the temperatures tested. For this purpose, the permeate flux and transmembrane pressure results were processed to obtain the membrane permeability (k) in each case. Knowing the relationship of this parameter with the membrane resistance (Equation (2)), R i was obtained for each of the tests carried out. The R M was obtained in the same way using the data of the clean membrane. Subsequently, and taking into account Equation (3), the R F was obtained in each case. The correlation rate corresponds to the data fit by forcing the line through the origin. The R F values should be approximately the same for all temperatures, since the fouling that has occurred on the membrane is the same for all three tests. These variations are due to the small membrane washout that has occurred over time during the three characterisation tests, as well as the effect of temperature, because it affects the viscosity of the filtered fluid.
As in the case of intrinsic resistance, the total resistance of the membrane decreases as the temperature increases. The fouling resistance values obtained were between 1.67 × 10 −3 and 1.01 × 10 −3 bar·m 2 ·h/L; the membrane resistance increased due to fouling by approximately 22%. Therefore, the fouling produced by the contents of the photobioreactor (including microalgae) is small, as has been found in previous studies [27].
Low fouling ensures that the lifetime of the ultrafiltration membrane will be longer because the permeability of the membrane is not significantly affected and therefore the costs associated with pumping do not increase [41].

Evaluation of the Effectiveness of the Cleaning Protocol and Recovery of Membrane Characteristics
To recover the filtering characteristics of the membrane, the cleaning protocol described above was carried out. After each cleaning step, the membrane was characterised in water to calculate the efficiency of each cleaning step (Table 5) according to Equation (4).
In all cases, the initial resistance of the membrane was recovered, and the overall effectiveness of the cleaning protocol was 100%. As can be seen from the results, in the first cleaning step, a recovery of approximately 42% was obtained. On the other hand, in the second cleaning step, the remaining 58% was recovered, reaching the initial R M values.
Both steps were effective and therefore necessary for complete cleaning of the fouled membrane in the photobioreactor characteristics. In the tests carried out at different temperatures, similar results were obtained, with variation in the values of the membrane resistance because of temperature. Therefore, the optimal temperature for the cleaning protocol is considered to be 30 • C, where a lower R M value is obtained. This protocol would be cost-effective at an industrial level because it is a discontinuous cleaning process which can be sized according to the membranes in large installations to minimise costs, limiting and spacing the cleaning to four times per year. This also highlights the usefulness of this cleaning protocol, which can be used on a small scale with very good results. Initially, only organic cleaning with sodium hypochlorite could be carried out, minimising costs, since from the trials it is expected to be at least 50% effective. Previous authors applied another cleaning protocol after the use of the membranes in a microalgae culture. Their protocol was less effective than the one applied in our study because the inorganic fouling was not successfully removed, so they suggested the use of citric acid to complete the cleaning, and this has been done in our study [27].

Analysis of the Operational TMP
To determine the optimal working TMP for our study system, fluxes (J) at different TMPs were analysed to fit the data to the Equation (5) model. In no case was it possible to fit the proposed equation, because the studied TMPs allowed by the filtration device were in the pressure-controlled zone, without reaching the values of the zone controlled by matter transfer. Therefore, it was decided to fit these data to the linear part of the equation. Figure 7 shows the data obtained for the three temperatures studied with the membrane clean before starting the study cycle and after finishing it (with the membrane soiled). From these representations, the data shown in Table 6 were obtained: The fouling that occurs is so small that the R M values obtained with water and with the photobioreactor contents are similar both at the beginning and at the end of the cycle. The low fouling obtained in the membrane may be due to the consortium of bacteria and algae that exist in the MPBR, decreasing the organic load present in the medium and therefore decreasing fouling [42]. Again, the behaviour is similar in the three studies at different temperatures, obtaining a variation only due to the influence of temperature. Membrane resistance increases between 28 and 13% with temperature.
Again, the resistance values obtained for the fouled membrane should be approximately the same for all temperatures, as the fouling that has occurred on the membrane is the same for all three tests. These variations are due to the small amount of membrane washout that has occurred over time during the three characterisation tests. As before, these differences are observed because of the effect of temperature on the viscosity of the filtered fluid.  The fouling that occurs is so small that the RM' values obtained with water and with the photobioreactor contents are similar both at the beginning and at the end of the cycle. The low fouling obtained in the membrane may be due to the consortium of bacteria and algae that exist in the MPBR, decreasing the organic load present in the medium and therefore decreasing fouling [42]. Again, the behaviour is similar in the three studies at different temperatures, obtaining a variation only due to the influence of temperature. Membrane resistance increases between 28 and 13% with temperature.
Again, the resistance values obtained for the fouled membrane should be approximately the same for all temperatures, as the fouling that has occurred on the membrane is the same for all three tests. These variations are due to the small amount of membrane As can be seen, the selected working TMP, which is 0.01 bar, is in the linear zone and therefore within the recommended range for the operation of this system, being in the pressure-controlled zone. This TMP is similar to that used in other studies in which photobioreactors with microalgae are used [20]. Working under these TMP conditions favours less fouling of the membrane. This effect is explained by the fact that at higher TMP values, greater membrane resistance is caused by the compression of the layer of cells deposited on it [43].

Effluent Quality
The quality of the effluent from the photobioreactor was assessed by analysing the TSS, turbidity, and quantification of nitrogen and phosphorus. Table 7 shows the average values of temperature, pH, conductivity, TSS, and the turbidities of influent and effluent. According to Regulation (EU) 2020/741 of 25 May 2020, which establishes minimum requirements for water reuse in agricultural irrigation, the effluent from the photobioreactor can be reused as reclaimed water quality class A, B, C, and D [44]. In other words, it complies with all the minimum requirements set out in that regulation.
According to Royal Decree 1620/2007, of 7 December 2007, which establishes the legal regime for the reuse of treated water and in relation to the results obtained for the concentration of TSS and turbidity, the effluent from our photobioreactor could be reused for urban uses of quality 1.1 and 1.2, agricultural uses of quality 2.1, 2.2, and 2.3, industrial uses of quality 3.1, recreational uses of quality 4.1 and 4.2, and environmental uses of quality 5.1, 5.2, and 5.3 [45]. Therefore, the effluent obtained satisfies the requirements established for the reuse of treated wastewater in all of the uses regulated in this Royal Decree, except for quality 3.2 (industrial uses in cooling towers and evaporative condensers). Table 8 shows the results obtained in the nutrient analysis: the average phosphorus and nitrogen in the influent and effluent, the removal yields, and the consumption of each nutrient. These last two parameters were calculated considering that the HRT in the photobioreactor is 2 days. The total nitrogen removal performance was slightly lower than that obtained by other authors [7]. This effect is probably due to the presence of autotrophic microorganisms in the wastewater. These consume the ammonium (reducing it in the photobioreactor) but generate nitrates, increasing their concentration [20]. However, the total nitrogen removal rate in the photobioreactor is higher than that obtained by other authors [46].

Nitrogen and Phosphorus Removal
Other authors report that working times at high SRTs influence the feeding preference of microalgae, observing that for long SRTs ammonium is preferred to nitrate, raising nitrate levels in the effluent [28]. The preference for nitrate or ammonium will depend on the species of microalgae grown in the photobioreactor and its operational conditions [7].
One of the possible solutions to improve the nutrient yield in the photobioreactor would be to increase the HRT. In this way, the dilution of the culture is not so great, and therefore, the renewal of the substrate is not so fast, thus improving the elimination yields, as the microalgae have more time to assimilate the nutrients [12].
On the other hand, it has been observed that phosphorus concentrations have been quite low in the influent of the photobioreactor. This event has probably conditioned the growth of the microalgae, acting as a limiting nutrient for the culture. To achieve higher nutrient removal yields under the same operating conditions, a higher amount of biomass would be necessary, and this would require a higher concentration of phosphorus in the medium. Despite this, the average removal efficiency obtained for phosphorus is higher than that obtained in other studies [25]. Low phosphate levels cannot be ignored as this effect can reduce nitrate uptake, limiting microalgae growth and nutrient removal performance in the photobioreactor [7,47].

Conclusions
The following conclusions were reached after 80 days under continuous operation (two days of HRT) of a microalgae membrane photobioreactor for the treatment of real urban wastewater effluent by culturing the microalgae present in it.
On one hand, by analyzing the results obtained, it was observed that optical density was the parameter that best described the growth kinetics of microalgae present in the photobioreactor. The zero-order model had the highest correlation rate of 0.9776, and a value of 8.921 × 10 −3 days −1 was obtained for the kinetic constant k.
On the other hand, with respect to the membrane, there was an increase in the total membrane resistance of approximately 22%. Therefore, under the operative conditions tested, the fouling was low. A two-phase cleaning protocol (inorganic and organic cleaning) was proposed to recover membrane permeability. The recoveries were approximately 42 and 58% in the inorganic and organic phases, respectively, for the temperatures tested (20,25, and 30 • C), but R M values were lower as the temperature increased, as would be expected. At all temperatures, the overall effectiveness of the cleaning protocol was 100%.
However, the optimum temperature tested for the cleaning protocol is considered to be 30 • C, where a lower R M value is obtained.
Regarding the effluent obtained and the removal of the nutrients present in the treated wastewater, total nitrogen and phosphorus removals in the algae membrane photobioreactor were 56.3 and 64.27%, respectively. Therefore, taking into account other water quality parameters, the effluent of the microalgae membrane photobioreactor could be reused for agricultural irrigation, according to Regulation (EU) 2020/741 of 25 May 2020.
Considering the above, a microalgae membrane photobioreactor could be an appropriate technology for urban wastewater regeneration to be included in a conventional activated sludge treatment plant. Funding: This research was funded by the Junta de Andalucía and the European Regional Development Fund (ERDF), grant number P18-TP-4732.

Data Availability Statement:
The data presented in this study are available on request from the corresponding author.