Ultrafiltration Harvesting of Microalgae Culture Cultivated in a WRRF: Long-Term Performance and Techno-Economic and Carbon Footprint Assessment

A cross-flow ultrafiltration harvesting system of a pre-concentrated microalgae culture was tested in an innovative anaerobic-based WRRF. The microalgae culture was cultivated in a membrane photobioreactor fed with the effluent from an anaerobic membrane bioreactor treating sewage. These harvested microalgae biomass was then anaerobically co-digested with primary and secondary sludge from the water line. Depending on the needs of this anaerobic co-digestion, the filtration harvesting process was evaluated intermittently over a period of 212 days for different operating conditions, mainly the total amount of microalgae biomass harvested and the desired final total solids concentration (up to 15.9 g·L1 with an average of 9.7 g·L1). Concentration ratios of 15-27 were obtained with average transmembrane fluxes ranged from 5 to 28 L·m2·h1. Regarding membrane cleaning, both backflushing and chemical cleaning resulted in transmembrane flux recoveries that were, on average, 21% higher than those achieved with backflushing alone. The carbon footprint assessment shows promising results as the GHG emissions associated with the cross-flow ultrafiltration harvesting process could be less than the emissions savings associated with the energy recovered from the biogas production from the anaerobic valorisation of the harvested microalgae.

Keywords:

Subject: Engineering - Chemical Engineering

1. Introduction

Worldwide industrial activities are facing great challenges to become more sustainable and “carbon free” by minimising emissions of greenhouse gases (GHG). Within industrial activities, large-scale wastewater treatment has scarcely considered sustainability in its processes so far, although its energy consumption is highly relevant, as it has prioritised effluents’ quality requirements. The energy consumption of urban wastewater treatment has been reported to account for above 1% of a country’s total energy consumption in Europe [1]. The current paradigm of water treatment needs to be shifted towards more sustainable and circular activities [2]. In fact, the new proposal for the updated wastewater treatment directive, which is currently being revised, aims to achieve energy neutrality by 2040 for all wastewater facilities above 10,000 p.e. [3]. For this reason, many researchers are on the lookout for alternative wastewater treatment systems that can be instrumental in the reduction of GHG emissions.

A green alternative whose interest has grown exponentially in the last decade is the cultivation of microalgae in combination with wastewater treatment. This alternative allows for the reduction of carbon emissions, the recovery of nutrients from water streams and the production of valuable microalgae biomass that contains approximately 50% of carbon, 10% of nitrogen and 1% of phosphorus [4−9]. Furthermore, wastewater valorisation through microalgae cultivation has been reported to reduce environmental impact and energy consumption compared to conventional treatment based on activated sludge systems [10].

Despite the many benefits of microalgae cultivation, current full-scale systems are scarce. One of the main issues of microalgae-based water resource recovery facilities (WRRFs) is the large surface areas that are required for successful cultivation of microalgae, as open reactors’ depths are shallow (i.e., 15-40 cm), and hydraulic retention times (HRT) are typically around 3-10 days [11−15]. To improve this, increasing number of authors are testing different options to make the technology more efficient. For instance, some authors are trying to develop instrumentation, control, and automation (ICA) systems to improve microalgae performance [16]. Other reactor configurations have been proposed to increase the photosynthetic efficiency of microalgae, for instance by using thin-layer reactors to increase light availability [17,18]. For example, Morales-Amaral et al. [19] produced higher amounts of biomass in a 2-cm-deep thin-layer reactor than in a raceway pond (43% more). Other authors have incorporated membrane separation technologies to decouple the HRT from the solids retention time (SRT) in order to increase the nutrient loading to the treatment system while maintaining the microalgae biomass in the reactor for longer periods of time, giving them more time to grow and increase productivity [20,21]. Luo et al. [22] managed to reduce the HRT up to 1 d (decoupled from SRT of 9-30 d) using membrane photobioreactor technology. This implied an increase on the nitrogen and phosphorus yields of 25-100%. In addition, Gao et al. [23] were able to increase biomass productivity when reducing the HRT from 6 to 2 d (SRT= 21 d). To achieve this decoupling, the microalgae biomass needs to be separated from the wastewater stream, which is not a simple process and requires the use of harvesting systems that are generally energetically and economically demanding [24,25]. A harvesting-step is also necessary to concentrate the microalgae biomass with the aim of obtaining added-value by-products, such as biogas from its anaerobic digestion [26], biocrude with high content of alkenes and alkanes [27], biodiesel [28], lipids [29], fatty acids methyl esters [30], nutraceutical applications as linoleic and linolenic acids [31], eicosapentaenoic acid (EPA) and docosahexanoic acid (DHA) [32] or pigments in a biorefinery approach [33]. Harvesting systems such as gravity sedimentation, coagulation-flocculation, electroflocculation, magnetic flocculation, bioflocculation, flotation, centrifugation and filtration have been widely reported (Table 1).

Filtration seems to be one of the most promising harvesting methods as it allows to obtain high concentrations (concentration ratios up to a maximum of 150, with final concentrations reaching up TSS of 170 g·L⁻¹), high quality water effluents that could be reused for irrigation; and it can avoid damage to the microalgae cells so that they could be returned to the cultivation unit [2,21,45]. Both dead-end (no retentate stream) and cross-flow (CF) ultrafiltration (UF) (retentate and permeate streams) have been reported in literature for microalgae harvesting [21, 39, 50−52] and [21, 39, 44, 45, 47, 50−53], respectively. However, the major concern with membranes is related to fouling, which significantly increases operating costs due to physical cleaning (shear stress from cross-flow velocity in cross-flow filtration, air/gas-sparging in dead-end filtration, backflushing, etc) and chemical cleaning which in turn reduces membrane lifetime [54−55]. Consequently, the use of filtration as harvesting process needs further optimisation to enhance the feasibility of microalgae cultivation technology. Zhao et al. [21] gathered energy consumption data for filtration harvesting processes from 2003 to 2022, reporting values up to 3.5 kWh·m⁻³ of treated microalgae culture, with some peak values higher, but the majority falling between 0.2 and 1 kWh·m⁻³ for permeate fluxes between 17 and 45 LMH. Particularly, energy requirements for cross-flow ultrafiltration of microalgae are reported to be between 0.17-2.23 kWh·m⁻³ for initial concentrations of the culture between 1 and 2 g·L⁻¹ and final concentration ratios ranging from 25-113, and transmembrane pressure (TMP) ranging from 0.31 to 2 bar and cross-flow velocities (CFV) from 0.17 to 2 m·s⁻¹ [21,44,47,50,53,56]. However, there is still scarce energy consumption data in membrane-based microalgae harvesting [21].

One advantage of harvesting by cross-flow ultrafiltration is that recent studies have reported that it can be used as a pre-treatment step for anaerobic valorisation of microalgae as it promotes their biodegradability (Table 1). For instance, Giménez et al. [26] concluded that the effect of the cross-flow ultrafiltration harvesting technique on the integrity of microalgae cell wall, in terms of viability and biodegradability, exhibited a noticeable effect that increased biodegradability. It was hypothesised that this effect was attributed to the induction of greater shear stress, which is not observed to the same extent in the case of dead-end ultrafiltration. The flexibility afforded by the difference in performance between the dead-end and cross-flow configurations in terms of cell wall integrity is one of the strengths of the filtration technology applied to microalgae harvesting, due to its ability to adapt to different potential uses of the harvested biomass.

The aim of this study is to evaluate the long-term performance of a cross-flow ultrafiltration process to harvest the microalgae biomass cultivated in a membrane photobioreactor for its further valorisation through anaerobic co-digestion. Specifically, this microalgae biomass is produced in an anaerobic-based WRRF platform for sulphate-rich sewage treatment [57]. In addition, this study presented a techno economic assessment and GHG emissions assessment at pilot-scale to facilitate the economic and environmental feasibility study of membrane-based harvesting. Specifically, for GHG emissions, the harvesting unit and a digestion unit for anaerobic valorisation of harvested microalgae were considered, considering different biogas valorisation scenarios.

2. Materials and Methods

2.1. Description of harvesting pilot plant

The layout of the WRRF is shown in Supplementary materials (Figure S1), which basically consisted of a combination of membrane technologies with anaerobic processes for the valorisation of organic matter from wastewater in the form of biogas with anaerobic membrane bioreactor (AnMBR) technology, and nutrient recovery using microalgae cultivation with membrane photobioreactor technology (MPBR). This MPBR unit used membranes to decouple the SRT from the HRT, which was demonstrated to promote microalgae activity significantly [45]. This initial harvesting stage consisted of two membrane tanks with a set of hollow-fibre ultrafiltration membranes (KMS Puron® Koch Membrane Systems) with a total membrane area of 3.44 m² and an average pore size of 0.03 µm. The performance of this first stage has been previously analysed in González-Camejo et al. [58] The primary microalgae species in the culture were Coelastrella and Desmodesmus, representing over 99% of the total eukaryotic cells [59].

Anaerobic co-digestion using AnMBR technology was used for the energy recovery of all waste streams from the water line (i.e., primary sludge, AnMBR sludge, and harvested microalgae). Prior to this, a second harvesting stage, using cross-flow ultrafiltration, was used to achieve the desired levels of total suspended solids (TSS) in the microalgae waste (from the MPBR unit) to be fed to the anaerobic co-digestion unit. This pilot-scale facility was located at the Carraixet WWTP (Spain) and used as influent the effluent from the pre-treatment stage of the full-scale plant. A more detailed description can be found elsewhere [57].

The cross-flow ultrafiltration pilot plant is shown in detail in Figure 1. The installation included a feed tank (0.7 m³), a retentate tank (0.04 m³), the cross-flow membrane module (in a vertical disposition, 1 m long), and, finally, the clean-in-place (CIP) tank (0.2 m³). This unit employed a tubular module with ultrafiltration membranes (HF 5.0-43-PM500, ROMICON® Koch Membrane Systems, USA). The pore size of the membranes had a molecular weight cut-off of 500 kDa, effectively retaining macromolecules. The total filtration surface covered 1 m², ensuring efficient separation and concentration of the microalgae biomass. The total cross-sectional area was 0.0003 m². This plant, which was operated in batch mode in cycles, allowed the removal of excess water, resulting in a highly concentrated algal biomass suitable for subsequent co-digestion processes [60].

2.2. Instrumentation and automation

The cross-flow harvesting system was automated to evaluate the filtration process properly (Figure 1). The following on-line sensors were installed: (i) one level transmitter in the feed tank LIT-1 (Waterpilot FMX167, Endress Hauser) to determine the feed tank level, LFT (m); (ii) one temperature probe in the retentate tank, TT-1 (RTD PT100 RS PRO, RS Components); (iii) one suspended solids sensor in the retentate tank, TSS-1 (SOLITAX ts-line sc LXV423.99.00100); (iv) one liquid flow meter in the retentate conduction (SITRANS FM MAG 1100 DN15, SIEMENS); (v) and three pressure gauges (PG-1 in the feed to the membrane, PG-2 in the retentate of the membrane and PG-3 in the backflushing pipe).

In terms of electromechanical equipment, there were 3 pumps: P-1 was the feed pump to the retentate tank in order to continuously replenish it with fresh culture to be concentrated and also to recirculate the fresh culture to the feed tank to guarantee proper mixing conditions; P-2 was the feed pump to the membrane module; whereas P-3 was the auxiliary pump to enable backflushing (Figure 1). Additionally, the plant was equipped with the following auxiliary equipment: an on/off solenoid valve (SOV-1) to allow the feed into the retentate tank; a level switch in the retentate tank, LS-1, acting on the solenoid valve SOV-1; and a valve (MV-1) to enable the collection of harvested microalgae after the completion of each batch's work cycle. The SOLITAX transmitter used in this study was equipped with colour correction and an automatic cleaning valve, which improved the representativeness and accuracy of the measured and recorded data. Moreover, this sensor was installed with a 30-degree deviation from the regular perpendicular angle to the photobioreactor surface to prevent probe fouling and reduce the noise of the signal. Furthermore, a comprehensive sensor cleaning protocol was implemented to improve and ensure stable data acquisition. The SOLITAX probe was calibrated once a week using laboratory data obtained from grab samples that were collected in duplicate from the retentate tank, at the beginning and end of each cycle, according to the Standard Methods [61]: method 2540 E. A linear correlation was established with the values recorded continuously by the SOLITAX probe (R2=0.9919; p-value < 0.05; n=38).

All sensors were connected to a programmable logic controller (PLC) for process automation, data acquisition and control. The PLC was communicated with a PC on which supervisory control and data acquisition software (SCADA) was installed to receive, transform, and perform calculations with the obtained data, as well as enable visualisation of all relevant process parameters.

2.3. Pilot plant operation and monitoring

The harvesting plant was operated in batch mode in cycles. In each cycle, the pre-concentrated culture from the MPBR plant stored in the feed tank was pumped (P-1) to the retentate tank until it was full as detected by the level switch (LS-1) which closed the feed valve to the retentate tank (SOV-1). The P-1 pump remained on to ensure proper mixing conditions in the feed tank. The culture in the retentate tank was then continuously fed to the cross-flow membrane producing a permeate which was stored in the CIP tank and a retentate stream which returned to the retentate tank. As the volume of the retentate tank decreased (due to the accumulation of permeate in the CIP), the feed pump (P-1) was automatically activated to replenish the retentate tank. This operation was carried out continuously until the desired level of TSS concentration in the retentate tank was reached (over 8 g TSS·L⁻¹). In some cycles the feed tank volume was exhausted, but the filtration cycle was not stop because the desired TSS concentration was not achieved. The required TSS concentration varied slightly depending on the requirements of the subsequent anaerobic co-digestion unit in the WRRF pilot plant.

To achieve this TSS concentration, the membrane was operated at a constant transmembrane pressure (TMP), whereas the transmembrane flux (J) remained variable. The TMP (bar) was calculated using Eq. 1 [62]:

T M P = \frac{P_{1} + P_{2}}{2} - P_{p e r m e a t e}

(1)

where P₁ represents the inlet gauge pressure (bar) to the membrane module, obtained by PG-1; P₂ is the outlet gauge pressure (bar), measured by PG-2 in the retentate stream; and P_permeate is the permeate gauge pressure (bar).

Cheryan [62] established that P_permeate should be zero when the permeate is open to the atmosphere, as in the present study. However, given the vertical configuration of the membrane and the outlet location at the top of the module, a constant P_permeate of 0.049 bar (0.5 m of water column) was considered.

Once the desired TSS concentration was reached, the harvested microalgae biomass was then removed, and the membrane was cleaned physically (backflushing) and sometimes chemically (every 6 operating cycles). In each cycle, the backflushing pump (P-3) was operated at its maximum capacity with a total backwash volume of 90 liters.

The cross-flow ultrafiltration plant was operated in batch cycles for 212 days (June to January), during which the WRRF pilot plant ran continuously. A total of 78 cycles were performed, two to three times per week (with some periods of non-operation), with an average duration of 25 hours per cycle. The cycle duration varied mainly according to the available pre-concentrated microalgae biomass from the MPBR pilot plant (total volume and TSS concentration), the TSS concentration required for the subsequent anaerobic co-digestion unit in the pilot plant, and membrane permeability. This initial TSS concentration (TSS_i) was measured off-line in the feed tank and varied depending on the operation of the MPBR plant, which in turn depended on several variables [58,63,64]. Table 2 summarises the average operating and outdoor conditions over the 78 cycles.

For each operating cycle, the following parameters were calculated to fully monitor the filtration process:

Cross flow velocity, CFV (m·s⁻¹), calculated as follows:

C F V = \frac{Q_{P - 2}}{A_{C F}}

(2)

where Q_P-2 is the feed flow rate to membrane module (m³·s⁻¹), A_CF is membrane cross-sectional area (m²).

2.: Permeate flow rate, Q_p (L·h⁻¹), was calculated hydraulically using L_FT data from LIT-1, i.e., feed tank level variation (which is equivalent to CIP tank level variation):

Q_{p} = \frac{|Δ L_{F T}| \cdot A_{F T}}{Δ t \cdot 1000}

(3)

where ΔL_FT was the variation of L_FT (m) in a set period within the work cycle, Δt (h); and A_FT was the feed tank area (m²).

3.: Transmembrane flux, J (LMH), was calculated as follows:

J = \frac{Q_{p}}{A_{F}}

(4)

where: Q_p was the permeate flow rate (L·h⁻¹), A_F was the filtration area of the membranes (m²).

4.: Standardised transmembrane flux at 20 °C, J₂₀ (LMH), calculated as follows:

J_{20} = J \cdot e^{- 0.0239 \cdot (T - 20)}

(5)

where T (ºC) is the temperature of the microalgae culture being concentrated in the harvesting system, measured by TT-1.

5.: Normalised transmembrane flux at 20 ºC, J₂₀:J_20,0,

J_{20} : J_{20,0} = \frac{J_{20}}{J_{20,0}}

(6)

where J_20,0 is the initial J₂₀ obtained at the start of the entire experiment (32.7 LMH).

6.: Membrane permeability standardised at 20 °C, K₂₀ (LMH·bar⁻¹) (Eq. 7):

K_{20} = \frac{J_{20}}{T M P}

(7)

7.: Backflush flow rate, Q_BF (L·min⁻¹), was calculated as follows:

Q_{B F} = \frac{V_{B F}}{t_{B F}}

(8)

where V_BF (L) was the backflush volume measured in a time t_BF (min).

8.: Transmembrane pressure during backflushing (TMP_BF) was calculated by Eq. 9:

{T M P}_{B F} = P_{B F -} (\frac{P_{1} + P_{2}}{2})

(9)

where P_BF (bar) is the discharge gauge pressure of pump P-3, measured by PG-3.

9.: Harvested microalgae culture biomass, M_TSS_HV (g), calculated as follows:

{M_T S S}_{H V} = {T S S}_{f} \cdot V_H V

(10)

where TSS_f was the TSS recorded by the TSS-1 sensor at the end of each cycle (g·m⁻³), and V_HV was the final volume in the retentate tank at the end of each cycle (m³).

10.: Harvesting rate HV_r (g TSS·m⁻²·h⁻¹), calculated as follows:

H V_r = \frac{{M_T S S}_{H V}}{{A_{F} \cdot t}_{o p e r a t i o n}}

(11)

where toperation is the cycle duration (h).

11.: Concentration ratio r (Eq. 12):

r = \frac{{T S S}_{f}}{{T S S}_{i}}

(12)

2.4. Energy and chemical reagents consumption

The energy consumption of the pumps P-1 and P-2 (W) was calculated adapting the equations proposed by Judd and Judd [65] and Ortiz-Tena et al. [66] (see Eq. 13 and Eq. 22). This energy model has already been applied to different UF membrane systems used in microalgae cultivation [58].

W_{P, j} = Q_{P, j} · {g · ρ}_{culture} · \frac{\{[{(\frac{(L + L_{eq}) · f · v^{2}}{D · 2 · g})}_{asp} + {(\frac{(L + L_{eq}) · f · v^{2}}{D · 2 · g})}_{imp}] + [Z_{1} - Z_{2}]\}}{η_{pump_j}}

(13)

where W_P,j is the power required for pump j (W), considering both the suction and discharge sections of the pump, using: the volumetric flow rate (Q_P,j in m³·s⁻¹), the acceleration of gravity (g in m·s⁻²), the density of the culture (ρ_culture in kg·m⁻³), the length of the pipeline (L in m), the pressure drops due to accidents expressed as equivalent length (L_eq in m), the velocity of the culture (v in m·s⁻¹), the friction factor (f, dimensionless), the cross-sectional diameter (D in m), the difference in elevation to overcome (Z₁-Z₂, in m), and the efficiency of pump j (η_{pump_j}, 0.60 [57]).

To calculate the ρ_culture (kg·m⁻³), the density of wet green microalgae biomass, ρ_a,w (kg·m⁻³) was calculated using Eq. 14 [67]:

ρ_{a, w} = x_{w} · ρ_{w} + (1 - x_{w}) \cdot ρ_{a, d}

(14)

where: x_w was the water mass content of wet green algal biomass, 0.82 [67]; ρ_a,d was density of dry green microalgae biomass, 1,400 kg·m⁻³ [67]; ρ_w was the water density.

The density of the culture ρ_culture can be obtained using Eq. 15:

ρ_{c u l t u r e} = \frac{ρ_{a, w} + ρ_{w}}{\frac{T S S}{ρ_{a, w}} + 1}

(15)

where: TSS (kg·m⁻³) was microalgae concentration in the culture.

Regarding the viscosity of water, μ_w (Pa·s), it was calculated using Eq. 16 [68], valid between 0 and 40 ºC, which aligned with the operating range of this study:

μ_{w} = μ_{w, 20} \cdot 10^{(\frac{20 - T}{T + 96} \cdot (1.2364 - 1.37 \cdot 10^{- 3} \cdot (20 - T) + 5.7 \cdot 10^{- 6} \cdot {(20 - T)}^{2}))}

(16)

where: μ_w,20 was 0.001002 Pa·s.

Mass fractional content in dry solids in the microalgae biomass, w_a,d, can be obtained by Eq. 17 [67]:

w_{a, d} = \frac{(1 - x_{w}) \cdot ρ_{a, d}}{ρ_{a, w}}

(17)

Microalgae volume fraction in the culture, Φ_v,a, was determined by Eq. 18 [67]:

Φ_{v, a} = \frac{T S S}{ρ_{a, w} \cdot w_{a, d}}

(18)

where TSS (kg·m⁻³) was microalgae concentration in the culture.

The viscosity of the microalgae culture can be calculated using Eq. 19, adapted from [67], which provides reliable values for the range Φ_v,a≤0.115, well above the maximum value in this study:

μ_{c u l t u r e} = μ_{w} \cdot {(1 - \frac{Φ_{v, a}}{Φ_{v, a}^{m}})}^{- 2}

(19)

where

Φ_{v, a}^{m}

was the maximum microalgae volume fraction, established in 0.637.

The friction factor (f) was calculated using the Swamee-James equation (Eq. 20):

f = \frac{0.25}{{[{l o g}_{10} (\frac{\frac{k}{D}}{3.7} + \frac{5.74}{{R e}^{0.9}})]}^{2}}

(20)

where Re is the dimensionless Reynolds number, k (m) is the internal roughness of the pipe, and D (m) is the cross-sectional diameter of the pipeline.

For the calculation of the power consumed by pump P-3 (W), used in the backflushing stages, Eq. 21 was employed:

W_{P, 3} = \frac{Q_{B F} \cdot {T M P}_{B F}}{0,6 \cdot η_{P - 3}}

(21)

Where: η_P-3 was the efficiency of the pump P-3, with a value of 0.60 Error! No se encuentra el origen de la referencia..

The total energy consumption EC (kWh) was determined (Eq. 22) by aggregating the individual energy consumptions of each pump (W_P,j in kW) multiplied by each operation time, t_j (h), for all three pumps:

E C = \sum (W_{P, j} \cdot t_{j})

(22)

Based on the total energy consumption for each cycle, three ratios of interest have been defined:

i.: the energy consumption ratio of the harvesting system (EC_m_TSS, in kWh·tTSS⁻¹) per tonne of harvested microalgae biomass (M_TSS_HV, in t) (Eq. 23);
ii.: the energy consumption ratio of the harvesting system (EC_v_HV, in kWh·m⁻³) per treated volume of pre-concentrated microalgae culture, i.e., the initial volume of the feed tank (V_i, in m⁻³) (Eq. 24);
iii.: the energy consumption ratio of the harvesting system (EC_v_WRRF, kWh·m⁻³) per treated volume of water in the WRRF pilot plant (V_WRRF treated to generate V_HV, in m⁻³) (Eq. 25):

{E C}_{m}_T S S = \frac{E C}{{M_T S S}_{H V}}

(23)

{E C}_{v}_H V = \frac{E C}{V_{i}}

(24)

{E C}_{v}_W R R F = \frac{E C}{V_W R R F}

(25)

The chemical cleanings of the membrane were effectively carried out using a dose of 200 ppm of sodium hypochlorite (NaClO) for 10 minutes, followed by rinsing with water only. As already mentioned, the chemical cleanings were performed every 150 hours of operation, i.e., 6 cycles (on average). The parameter NaOClC_v_WRRF (g Cl·m⁻³) was defined to determine the sodium hypochlorite reagent consumption ratio per m³ treated in the WRRF (Eq. 26).

{N a O C l C}_{v}_W R R F = \frac{M_C l}{{V_W R R F}_{C C}}

(26)

where: M_Cl (g Cl) is the NaClO mass used in chemical cleanings expressed as active chlorine, and V_WRRF_CC (m^–3) is the amount of water treated in the WRRF since the last chemical cleaning.

Finally, the more relevant operating costs for the cross-flow ultrafiltration pilot plant were estimated, considering energy consumption and the use of sodium hypochlorite in chemical cleanings, OPEX_EC+Cl (€·m⁻³ treated in WRRF). The cost for electricity was the average for the UE-27 in October 2023 (0,178 €·kWh^–1) [69]. The cost for hypoclorite in Europe was taken from [70] (0.013 €·gCl^–1).

2.5. GHG estimation.

To estimate the GHG emissions associated with the cross-flow ultrafiltration process, the methodology proposed by Parravacini et al. [71] was considered adapted similarly to previous works [60], and taking into account [72]. This methodology includes the estimation of both direct (GHG_direct) and indirect (GHG_indirect) emissions related to the energy consumption of the harvesting unit and the energy consumption and biogas production of the anaerobic digestion of the harvested microalgae under mesophilic conditions. The biogas production from the anaerobic digestion of harvested microalgae biomass was determined based on the experimental data from previous studies [26,73]. In addition, the increase in methane yield due to the effect of cross-flow ultrafiltration on the microalgae valorisation (i.e., higher anaerobic biodegradability), according to ¡Error! No se encuentra el origen de la referencia., was also considered. In particular, this effect was taken into account by considering the enhancement of biogas production according to the concentration ratio (r) of the harvested microalgae, and it has been presented in Supplementary Materials (Figure S2). Methane yield correlated directly with the concentration ratio (r) for r values below 13.3. In contrast, for r values above 13.3 methane yield remained constant at 0.294 Nm⁻³·kgTSS⁻¹.

As for biogas valorisation, three scenarios were established: (i) Scenario 1, involving the use of high-efficiency cogeneration technology with a combined heat and power (CHP) system; (ii) Scenario 2, considering the upgrading to biomethane using conventional membrane technology and its subsequent injection into the grid; and (iii) Scenario 3, assessing the upgrading to biomethane using microalgae culture, which increased methane production during the anaerobic digestion (AD) stage due to the additional production of algal biomass at the upgrading unit, and subsequent injection into the grid, but this technology had not the possibility of heat recovery.

To determine the net energy demand per tonne of harvested microalgae biomass, E (kWh·tTSS⁻¹) by means of (Eq. 27), was necessary to estimate: Q_demand was the total thermal energy demand (kWh·tTSS⁻¹) (Eq. 28); W_demand was the total electrical energy demand (kWh·tTSS⁻¹) (Eq. 29); and Q_BM was the biomethane production in Scenarios 2 and 3 (kWh·tTSS⁻¹), calculated by (Eq. 30) if Q_demand was positive, and by (Eq. 31) if Q_demand was zero or negative.

E (k W h \cdot {t T S S}^{- 1}) = Q_{d e m a n d} + W_{d e m a n d} + Q_{B M}

(27)

Q_{d e m a n d} (k W h \cdot {t T S S}^{- 1}) = Q_{T O T} - Q_{r e c o v e r e d}

(28)

W_{d e m a n d} (k W h \cdot {t T S S}^{- 1}) = W_{T O T} - W_{r e c o v e r e d}

(29)

Q_{B M} (k W h \cdot {t T S S}^{- 1}) = (Q_{d e m a n d} - Q_{B G}) \cdot φ_{u p g r a d i n g}

(30)

Q_{B M} (k W h \cdot {t T S S}^{- 1}) = - Q_{B G} \cdot φ_{u p g r a d i n g}

(31)

where Q_TOT is the heat required by the anaerobic digestion process per tonne of harvested microalgae biomass (kWh·tTSS⁻¹); Q_recovered is the heat recovered in the CHP system for Scenario 1, or the heat recovered in the upgrading stage for Scenario 2 (kWh·tTSS⁻¹); W_TOT are the electrical consumptions of the equipment of the cross-flow ultrafiltration pilot plant, the anaerobic digester and the valorisation system of each scenario per tonne of harvested microalgae biomass (kWh·tTSS⁻¹); W_recovered is the electricity recovered in the CHP stage for Scenario 1 per tonne of harvested microalgae biomass (kWh·tTSS⁻¹); Q_BG is the gross production of raw biogas per tonne of harvested microalgae biomass; φ_upgrading is the efficiency of the upgrading process. Ratios, efficiencies, and emission factors are taken from [60, 71, 74−77].

Direct GHG emissions, indirect GHG emissions, and total GHG emissions per tonne of harvested microalgae biomass (M_TSS_HV, in t) were calculated using the equations 32-35.

{G H G}_{d i r e c t, s c . 1} (k g {C O}_{2} e \cdot {t T S S}^{- 1}) = M_{B G} \cdot {E F}_{{C H}_{4, s c 1}} \cdot 28

(32)

{G H G}_{d i r e c t, s c 2 & 3} (k g {C O}_{2} e \cdot {t T S S}^{- 1}) = M_{B M} \cdot {E F}_{{C H}_{4, s c, i}} \cdot 28

(33)

{G H G}_{i n d i r e c t} (k {g C O}_{2} e \cdot {t T S S}^{- 1}) = {(Q_{B M} + Q}_{d e m a n d}) \cdot {E F}_{n a t u r a l g a s} + W_{d e m a n d} \cdot {E F}_{e l e c t r i c i t y}

(34)

{G H G}_{t o t a l} = {G H G}_{d i r e c t} + {G H G}_{i n d i r e c t}

(35)

where M_BG is the calculated gross production of raw biogas per tonne of harvested microalgae biomass (kgCH₄·tTSS⁻¹); MBM is the production of biomethane per tonne of harvested microalgae biomass (kgCH₄·tTSS⁻¹);

{E F}_{{C H}_{4, i}}

is the methane losses emission factor for each scenario; Q_demand is the total thermal energy demand per tonne of harvested microalgae biomass (kWh·tTSS⁻¹) ; W_demand is the total electrical energy demand per tonne of harvested microalgae biomass (kWh·tTSS⁻¹); EF_{natural_gas} is the specific emission factor for fossil natural gas from the grid in Europe; and EF_electricity is the specific emission factor of European power companies.

3. Results

3.1. Filtration performance.

The harvesting pilot plant was in operation for a total of 1968 hours, intermittently distributed over a total of 78 cycles during a 212-day of operation of the WRRF pilot plant. Table 3 shows the average values of the main filtration performance parameters, where the resulting J and concentration ratios are in the range of other similar studies [21,50].

The filtration performance varied mainly according to the available pre-concentrated microalgae biomass from the MPBR pilot plant (TSS_i·V_i), the TSS concentration required for the subsequent anaerobic co-digestion unit (TSS_f), and the initial membrane permeability of each cycle depending on the membrane cleaning previously applied (physical or physical+chemical), in terms of normalised J₂₀:J_20,0. These operating conditions and the resulting filtration performance are shown in Figure 2a with the evolution of the TSS concentration in the retentate tank during each operating cycle, and in Figure 2b for 20 °C-standardised transmembrane flux normalised by the initial flux (J₂₀:J_20,0); concentration ratio (r) and the total mass culture to be harvested (total volume (V_i) per total suspended solids (TSS_i) in the feed tank). Both the total mass culture harvested (102-700 g) and the final TSS (8.0-15.9 g·L⁻¹) in each cycle varied significantly in some periods which strongly influenced the performance of the filtration cycle, in terms of the duration of the cycle and the final values obtained for the normalised transmembrane flux (Figure 2b) at the end of the cycle. For instance, the cycle with the highest TSSf and total mass culture harvested (cycle at 149-150 operating days) lasted 74 hours and J₂₀:J_20,0 reached the lowest values of the operating period, i.e., 0.08.

Backflushing stages were performed at the end of each cycle, and in Figure 2b the dashed lines indicate the end of the cycles when chemical cleaning was also performed (backflush followed by chemical cleaning). Within each cycle, J₂₀:J_20,0 started with a relatively high value and decreased as the cycle progresses, due to the increase in TSS concentration in the retentate tank. After each BF stage, when considering consecutive cycles with similar TSS_i, J₂₀:J_20,0 at the beginning of the cycle partially recovered compared to the initial one of the previous cycle. However, full recovery was not achieved because backflushing only removed reversible fouling [58,78]. Indeed, an accumulation of the irreversible fouling is indirectly observed with the decrease of J₂₀:J_20,0 at the beginning of the consecutive cycles between the dashed lines (between chemical cleanings). In contrast, when a backflush followed by a chemical cleaning was performed, a greater recovery of filtration capacity can be observed, expressed as J₂₀:J_20,0. The average of the initial values for J₂₀:J_20,0 obtained from cycles immediately after a combination of backflushing and chemical cleaning was 0.98±0.01, in contrast to the initial values from cycles after backflushing alone, which were 0.81±0.10.

In Figure 3, detailed figures for the period between day 54 and day 78 of operation are presented (cycles 19 to 33), in an arbitrary manner, to observe the evolution of J₂₀:J_20,0, concentration ratio, biomass to be harvested, and TSS in the retentate stream over only a few batch-operating cycles. In Figure 3b, the partial effect of the backflushing at the end of each cycle on the J₂₀:J_20,0 recovery is evident, as well as the significantly higher effect of the chemical cleaning performed after cycles 26 and 32. Figure 3a shows that the evolution of the TSS concentration in the retentate tank had an inflection point in each cycle. This inflection point corresponds to the time when the feed tank volume was exhausted, removing the dilution effect and reducing the retentate volume and the transmembrane flux as seen by the slower decline of J₂₀:J_20,0 (Figure 3b). Cycles 22 to 26 (from day 58 of operation prior to chemical cleaning) shown in Figure 3 were operated for a similar total culture mass harvested but decreasing concentration ratio, however a gradual increase in cycle time was observed as the values of J₂₀:J_20,0 decreased due to the accumulation of irreversible fouling.

One notable strength of the cross-flow ultrafiltration technology as second harvesting step lies in its robustness and reliability during operation, consistently ensuring the attainment of the 8 g·L⁻¹ target. This target was necessary for the overall functionality of the WRRF pilot plant, specifically in the anaerobic digestion process [57]. In analogous experiments involving dead-end hollow-fibre ultrafiltration membrane modules as the second harvesting step, the 8 g·L⁻¹ objective was not consistently achieved (data not shown). This discrepancy was likely attributed to variations in the TSS concentrations of the pre-concentrated feed stream or the specific of microalgae species present. However, it is well known that the energy demand of cross-flow operation is higher than that of dead-end operation. There is therefore room for improvement.

3.2. Techno-economic and carbon footprint assessment

Table 4 shows the values of the energy ratios and chemical reagent consumption rates defined in equations 24 to 27 for the entire operating period.

The energy ratios obtained in this study for an average concentration ratio of the harvested biomass of 21 are within the normal range reported in the literature for this type of biomass and technology [21]. It should be noted that the present study is a pilot plant performance of almost 2,000 h of batch operation over a 7-month operation of the WRRF pilot plant. In particular, the obtained energy demand is also in the range of other membrane-based microalgae harvesting systems, i.e., 0.17-2.23 kWh·m⁻³ [21,50]. Other harvesting technologies are even more energy intensive, such as centrifugation, which can consume up to 8 kWh·m⁻³ to achieve 90% of efficiency [21]. Important to note, is that there is room for improvement in the resulting filtration performance, for instance by means of optimising physical and chemical cleanings frequency.

Figure 4 shows, for a concentration ratio of 21, the results of the energy balance (Q_demand, W_demand and Q_BM), and the GHG assessment (GHG_direct, GHG_indirect) related to the harvesting unit and the anaerobic digestion unit for the harvested biomass for the three biogas valorisation scenarios considered.

As Figure 4a shows, Scenario 1 would not be able to cover all the electrical consumption (W_demand), but would produce a significant surplus of heat, Q_demand. In the case of Scenario 2, heat recovery from the conventional membrane upgrading would cover the heat requirements of the mesophilic digestion unit, and the substantial W_demand requirements would be greatly exceeded by the primary energy content of the biomethane injected into the grid, Q_BM. The difference between Scenarios 2 and 3 in terms of Q_demand can be explained by the lack of thermal recovery in microalgae-based harvesting systems. In Figure 4b, GHG emissions reflect the energy assessment due to the savings related to the fossil energy sources substituted by these renewable sources. Scenario 3 would be the most favourable one both in terms of energy recovery and GHG emissions savings. This improved performance could be primarily attributed to the additional biomethane production from CO₂ absorbed in the microalgae-based upgrading stage compared to both Scenarios 1 and 2, along with the low energy consumption of the upgrading process and the high methane recovery efficiency (exceeding 97.6%), compared to a combined efficiency of approximately 68% for co-generation in Scenario 1.

Figure 5 illustrates the net GHG emissions for all batch cycles vs. the concentration ratio (r), in the case of Scenario 3, chosen as an example since the trend is proportional for all scenarios. Given that most of the employed (r) values surpass the maximum methane yield obtained from previous results, i.e. the same methane yield is considered for r values above 13.3 [26,73], the most favourable GHG emission values are observed for relatively low (r) values, around 13-17. These results are promising, as the associated emissions could remain below zero for cross-flow ultrafiltration membranes over a wide range of operating conditions, thereby resulting in net emissions savings.

4. Conclusions

Cross-flow ultrafiltration appeared to be a robust technology for harvesting pre-concentrated microalgae biomass, which operated efficiently under variable operating conditions and was able to achieve biomass concentrations of up to 16 gTSS·L⁻¹ (concentration ratio of up to 39.5). The performance of cross-flow ultrafiltration was comparable to other membrane-based harvesting systems in terms of transmembrane flux, concentration ratios and energy requirements. A backflushing step was performed after each batch cycle, and the periodic combination of backflushing with chemical cleaning was necessary for long-term operation. The recovery of transmembrane flux after physical and chemical cleaning was significant, 21% higher than that achieved with backflushing alone. Based on these results, further optimisation is required to improve membrane productivity in terms of CAPEX and OPEX by optimising physical and chemical cleaning. If the harvested microalgae biomass were valorised for biogas production, the direct and indirect GHG emissions associated with the harvesting process could be less than the emissions savings associated with the energy recovered through the production of this renewable energy source, especially if the biogas is upgraded through a microalgae cultivation process.

Supplementary Materials

The following supporting information can be downloaded at the website of this paper posted on Preprints.org, Figure S1: Flow diagram of the WRRF scheme; Figure S2: Variation of biodegradability of the microalgae biomass and methane yield production with respect to concentration ratio (r).

Author Contributions

Conceptualisation, Juan Mora-Sánchez and Maria Ruano; Data curation, Juan Mora-Sánchez and Guillermo Noriega-Hevia; Formal analysis, Juan Mora-Sánchez; Funding acquisition, Aurora Seco; Investigation, Juan Mora-Sánchez, Josué González-Camejo and Guillermo Noriega-Hevia; Methodology, Juan Mora-Sánchez and Maria Ruano; Project administration, Aurora Seco; Resources, Aurora Seco; Supervision, Maria Ruano; Validation, Juan Mora-Sánchez and Maria Ruano; Visualisation, Juan Mora-Sánchez; Writing – original draft, Juan Mora Sánchez and Josué González-Camejo; Writing – review & editing, Josué González-Camejo and Maria Ruano. All authors have read and agreed to the published version of the manuscript.

Funding

This research work was supported by the Science and Innovation Spanish Ministry (Projects CTM2014-54980-C2-1-R/C2-2-R) and the European Regional Development Fund (ERDF), which are gratefully acknowledged. Generalitat Valenciana supported this study via fellowship APOTI/2016/56 to the first author. The Science and Innovation Spanish Ministry have also supported this study via a pre-doctoral FPU fellowship to the second author (FPU14/05082). The authors would also like to acknowledge the support received from the Universitat Politècnica de València via a pre-doctoral FPI fellowship to the third author, as well as the financial aid received from the European Climate KIC association for the ‘MAB 2.0’ Project (APIN0057_2015-3.6-230_P066-05).

Data Availability Statement

The datasets generated during and/or analysed during the current study are confidential.

Acknowledgments

The support from Entitat Pública de Sanejament d’Aigues Residuals (EPSAR) of Comunitat Valenciana is gratefully acknowledged.

Conflicts of Interest

The authors declare no conflict of interest.

Glossary

A_CF	Membrane cross-sectional area
AD	Anaerobic digestion
A_F	Filtration area of the membranes
A_FT	Feed tank area
AnMBR	Anaerobic membrane bioreactor
CF	Cross-flow
CF-UF	Cross-flow ultrafiltration
CFV	Cross-flow velocity
CHP	Combined heat and power system
CIP	Clean-in-place
D	Cross-sectional diameter
E	Energy balance calculation
EC	Total energy consumption
EC_m_TSS	Energy consumption ratio of the harvesting system per harvested microalgae biomass
EC_v_HV	Energy consumption ratio of the harvesting system per treated volume of pre-concentrated microalgae culture
EC_v_WRRF	Energy consumption ratio of the harvesting system per treated volume of water in the WRRF
EF_CH4	Methane losses emission factor
EF_electricity	Specific emission factor of european power companies
EF_{natural_gas}	Specific emission factor for fossil natural gas from the grid, in Europe
f	Friction factor
FIT	Liquid flow-rate transmitter
g	Acceleration of gravity
GHG	Greenhouse gas
GHG_direct	Direct greenhouse gases emissions
GHG_indirect	Indirect greenhouse gases emissions
GHG_total	Total greenhouse gases emissions
HRT	Hydraulic retention time
HV_r	Harvesting rate
ICA	Instrumentation, control and automation
J	Transmembrane flux
J₂₀	20 ºC-standardised transmembrane flux
J_20,0	Initial 20ºC-standardised transmembrane flux at the inception of the entire experiment
J₂₀:J_20,0	Normalised transmembrane flux at 20 ºC
k	Internal roughness of the pipe
K₂₀	20 ºC-standardised permeability
l	Length of the pipeline
l_eq	Pressure drops due to accidents expressed as equivalent length
L	Level
LIT	Level indicating transmitter
LMH	Litter per square meter and hour
LS	Level switch
M_BG	Gross production of raw biogas, expressed as methane mass
M_BM	Biomethane production, expressed as methane mass
M_Cl	Sodium hypochlorite reagent mass used in chemical cleaning
M_TSS_HV	Harvested microalgae culture biomass
MPBR	Membrane photobioreactor
MV	Valve
NaOClC_v_WRRF	Sodium hypochlorite reagent consumption ratio per m3 treated in the WRRF
OPEX_EC+Cl	Operating costs for energy consumption and sodium hypochlorite for the CF-UF pilot plant per m³ treated in the WRRF
p.e.	Population equivalent
PG	Pressure gauge
P_j	Pressure at point j
P-j	Pump number j
PLC	Programmable logic controller
Q_BG	Gross production of raw biogas, expressed as primary energy
Q_BM	Biomethane production, expressed as primary energy
Q_demand	Total thermal energy demand
Q_j	Flow-rate for pump or stream j
Q_recovered	Heat recovered by the biogas valorisation system
Q_TOT	Heat required by the anaerobic co-digestion process
r	Concentration ratio
SCADA	Supervisory control and data acquisition software
SOV	On/off solenoid valve
SRT	Solids retention time
T	Temperature
TMP	Transmembrane pressure
TSS	Total suspended solids
TSS_f	Final TSS concentration for an operation cycle
TSS_i	Initial TSS concentration for an operation cycle
TT	Temperature transmitter
UF	Ultrafiltration
v	Velocity
V_HV	Final volume in the retentate tank for an operation cycle
V_WRRF	Volume of water treated in the WRRF to generate V_HV
V_WRRF_CC	Volume of water treated in the WRRF since the last chemical cleaning
V_i	Initial volume in the feed tank for an operation cycle
w_a,d	Mass fractional content in dry solids in green microalgae biomass
W_demand	Total electrical energy demand
W_P,j	Power required for pump j
W_recovered	Electricity recovered by the biogas valorisation system
WRRF	Water resource recovery facility
x_w	Water mass content of wet green microalgae biomass
Z_j	Elevation of point j
ΔL_j	Variation of level j
Δt	Variation of time
$Φ_{w, a}^{m}$	Maximum microalgae volume fraction in the culture
Φ_ω,α	Microalgae volume fraction in the culture
η_{pump_j}	Efficiency of pump j
φ_upgrading	Efficiency of the upgrading process
μ	Viscosity
ρ	Density

References

Walker, N.L. , Williams, A.P., Styles, D. Pitfalls in international benchmarking of energy intensity across wastewater treatment utilities. J. Environ. Manag. 2021, 300, 113613. [Google Scholar] [CrossRef] [PubMed]
Foglia, A. , González-Camejo, J., Radini, S., Sgroi, M., Li, K., Eusebi, A.L., Fatone, F., 2023. Transforming wastewater treatment plants into reclaimed water facilities in water-unbalanced regions. An overview of possibilities and recommendations focusing on the Italian case. Journal of Cleaner Production, Volume 410, 137264. [CrossRef]
European Comission. Proposal for a DIRECTIVE OF THE EUROPEAN PARLIAMENT AND OF THE COUNCIL concerning urban wastewater treatment, 2022. https://environment.ec.europa.eu/publications/proposal-revised-urban-wastewater-treatment-directive_en.
Yuan, S. , Lei, W., Cen, Y., Liu, Q., Liu, J., Fu, J., Han, Y. Economic analysis of global microalgae biomass energy potential. Science of The Total Environment, 2023, Volume 899, 165596. [CrossRef]
Khan, S.; Thaher, M.; Abdulquadir, M.; Faisal, M.; Mehariya, S.; Al-Najjar, M.A.A.; Al-Jabri, H.; Das, P. Utilization of Microalgae for Urban Wastewater Treatment and Valorization of Treated Wastewater and Biomass for Biofertilizer Applications. Sustainability 2023, 15, 16019. [Google Scholar] [CrossRef]
Olabi, A.G. , Shehata, N., Sayed, E.T., Rodriguez, C., Anyanwu, R.C., Russell, C., Abdelkareem, M.A. Role of microalgae in achieving sustainable development goals and circular economy. Science of The Total Environment, 2023, Volume 854, 158689. [CrossRef]
Shilton, A. Pond Treatment Technology. IWA Publishing, 2006. [CrossRef]
Reynolds, C.S. The ecology of phytoplankton (ecology, biodiversity and conservation). Cambridge: Cambridge University Press; 2006. [Google Scholar] [CrossRef]
Solovchenko, A. , Ismagulova, T. , Lukyanov, A., Vasilieva, S., Konyukhov, I., Pogosyan, S., Lobakova, E., Gorelova, O.A. Luxury phosphorus uptake in microalgae. J Appl Phycol, 2019, 31, 2755–2770. [Google Scholar] [CrossRef]
Acién Fernández, F.G. , Gómez-Serrano, C., Fernández-Sevilla, J.M. Recovery of Nutrients From Wastewaters Using Microalgae. Frontiers in Sustainable Food Systems, 2018, 2:59. [CrossRef]
Nagarajan, D. , Lee, D.-J., Chen, C.-Y., Chang, J.-S. Resource recovery from wastewaters using microalgae-based approaches: A circular bioeconomy perspective. Bioresource Technology, 2020, Volume 302,122817. [CrossRef]
Arbib, Z. , de Godos, I. , Ruiz, J., Perales, J.A. Optimization of pilot high rate algal ponds for simultaneous nutrient removal and lipids production. Sci. Total Environ., 2017, 589, 66–72. [Google Scholar] [CrossRef] [PubMed]
García-Galán, M.J. , Gutiérrez, R. , Uggetti, E., Matamoros, V., García, J., Ferrer, I. Use of full-scale hybrid horizontal tubular photobioreactors to process agricultural runoff. Biosyst. Eng., 2018, 166, 138–149. [Google Scholar] [CrossRef]
Rossi, S. , Mantovani, M., Marazzi, F., Bellucci, M., Casagli, F., Mezzanotte, V., Ficara, E. Microalgae cultivation on digestate: Process efficiency and economics. Chemical Engineering Journal, 2023, Volume 460, 141753. [CrossRef]
Xu, X. , Gu, X. , Wang, Z., Shatner, W., Wang, Z. Progress, challenges and solutions of research on photosynthetic carbon sequestration efficiency of microalgae. Renew. Sust. Energy Rev., 2019, 110, 65–82. [Google Scholar] [CrossRef]
Oruganti, R.K. , Biji, A.P., Lanuyanger, T., Show, P.L., Sriariyanun, M., Upadhyayula, V.K.K., Gadhamshetty, V., Bhattacharyya, D. Artificial intelligence and machine learning tools for high-performance microalgae wastewater treatment and algal biorefinery: A critical review. Science of The Total Environment, 2023, Volume 876, 162797. [CrossRef]
Clagnan, E. , Dell'Orto, M., Štěrbová, K., Grivalský, T., Manoel, J.A.C., Masojídek, J., D'Imporzano, G., Acién-Fernández, F.G., Adani, F. Impact of photobioreactor design on microalgae-bacteria communities grown on wastewater: Differences between thin-layer cascade and thin-layer raceway ponds. Bioresource Technology, 2023, Volume 374,128781. [CrossRef]
De Vree, J.H. , Bosma, R., Janssen, M. et al. Comparison of four outdoor pilot-scale photobioreactors. Biotechnol Biofuels 8, 215 (2015). [CrossRef]
Morales-Amaral, M.M. , Gómez-Serrano, C., Acién, F.G., Fernández-Sevilla, J.M., Molina-Grima, E. Outdoor production of Scenedesmus sp. in thin-layer and raceway reactors using centrate from anaerobic digestion as the sole nutrient source. Algal Research, 2015, Volume 12, Pages 99-108. [CrossRef]
Luo, Y. , Le-Clech, P., Henderson, R.K. Assessment of membrane photobioreactor (MPBR) performance parameters and operating conditions. Water Research, 2018, Volume 138, Pages 169-180. [CrossRef]
Zhao, Z. , Muylaert, K., Vankelecom, I.F.J. Applying membrane technology in microalgae industry: A comprehensive review. Renewable and Sustainable Energy Reviews, 2023, Volume 172, 113041. [CrossRef]
Luo, Y. , Le-Clech, P., Henderson, R.K. Assessment of membrane photobioreactor (MPBR) performance parameters and operating conditions. Water Res., 2018 138, 169-180. [CrossRef]
Gao, F. , Peng, Y.Y., Li, C., Cui, W., Yang, Z.-H., Zeng, G.-M., 2018. Coupled nutrient removal from secondary effluent and algal biomass production in membrane photobioreactor (MPBR): Effect of HRT and long-term operation. Chemical Engineering Journal, Volume 335, Pages 169-175. [CrossRef]
Alkarawi, M.A.S. , Caldwell, G.S., Lee, J.G.M. Continuous harvesting of microalgae biomass using foam flotation. Algal Research, 2018, Volume 36, Pages 125-138. [CrossRef]
Liu, Z. , Hao, N., Hou, Y., Wang, Q., Liu, Q., Yan, S., Chen, F., Zhao, L. Technologies for harvesting the microalgae for industrial applications: Current trends and perspectives. Bioresource Technology, Volume 387, 2023, 129631. [CrossRef]
Giménez, J.B. , Bouzas, A., Carrere, H., Steyer, J.-P., Ferrer, J., Seco, A., 2018. Assessment of cross-flow filtration as microalgae harvesting technique prior to anaerobic digestion: Evaluation of biomass integrity and energy demand. Bioresource Technology, Volume 269, Pages 188-194. [CrossRef]
Das, P. , Khan, S., Thaher, M., AbdulQuadir, M., Hoekman, S.K., Al-Jabri, H. Effect of harvesting methods on the energy requirement of Tetraselmis sp. biomass production and biocrude yield and quality. Bioresource Technology, 2019, Volume 284, Pages 9-15. [CrossRef]
Liu, Q. , Zhang, M., Lv, T., Chen, H.J., Chika, A.O., Xiang, C.L., Guo, M.X., Wu, M.H., Li, J. J., Jia, L.S. Energy-producing electro-flocculation for harvest of Dunaliella salina. Bioresour. Technol. 2017, 241, 1022–1026. [Google Scholar] [CrossRef] [PubMed]
Kim, K. , Shin, H., Moon, M., Ryu, B.G., Han, J.I., Yang, J.W., Chang, Y.K. Evaluation of various harvesting methods for high-density microalgae, Aurantiochytrium sp. KRS101. Bioresour. Technol. 2015, 198, 828–835. [Google Scholar] [CrossRef] [PubMed]
Mahesh, R. , Naira, V.R., Maiti, S.K. Concomitant production of fatty acid methyl ester (biodiesel) and exopolysaccharides using efficient harvesting technology in flat panel photobioreactor with special sparging system via Scenedesmus abundans. Bioresour. Technol. 2019, 278, 231–241. [Google Scholar] [CrossRef] [PubMed]
Lv, J.P. , Zhao, F., Feng, J., Liu, Q., Nan, F.R., Xie, S.L. Extraction of extracellular polymeric substances (EPS) from a newly isolated self-flocculating microalga Neocystis mucosa SX with different methods. Algal Res. 2019, 40, 101479. [Google Scholar] [CrossRef]
Chen, Z.H. , Shao, S. S., He, Y.J., Luo, Q.Q., Zheng, M.M., Zheng, M.Q., Chen, B.L., Wang, M.Z. Nutrients removal from piggery wastewater coupled to lipid production by a newly isolated self-flocculating microalga Desmodesmus sp. PW1. Bioresour. Technol., 2020, 302, 122806. [Google Scholar] [CrossRef] [PubMed]
Rajesh-Banu, J. , Preethi, Kavitha, S., Gunasekaran, M., Kumar, G. Microalgae based biorefinery promoting circular bioeconomy-techno economic and life-cycle analysis. Bioresource Technology, 2020, Volume 302, 122822. [CrossRef]
Baroni, E.G. , Yap, K.Y., Webley, P.A., Scales, P.J., Martin, G.J.O.The effect of nitrogen depletion on the cell size, shape, density and gravitational settling of Nannochloropsis salina, Chlorella Sp. (marine) and Haematococcus pluvialis. Algal Research, 2019; 39, 101454. [Google Scholar] [CrossRef]
Japar, A.S. , Azis, N.M., Takriff, M.s., Haiza, M.Y.N., Yasin, M.H.M. Application of Different Techniques to Harvest Microalgae. Trans Sci Technol, 2017; 4, 98–108. [Google Scholar]
Razzak, S.A. , Ali, S.A.M., Hossain, M.M., deLasa, H. Biological CO2 fixation with production of microalgae in wastewater – A review. Renewable and Sustainable Energy Reviews, 2017, Volume 76, Pages 379-390. [CrossRef]
Sun, J. , Jiang, S., Yang, L., Chu, H., Peng, B.-Y., Xiao, S., Wang, Y., Zhou, X., Zhang, Y., 2023. Microalgae wastewater recycling: Suitability of harvesting methods and influence on growth mechanisms. Science of The Total Environment, Volume 859, Part 2, 160237. [CrossRef]
Dai, D., Qv, M., Liu, D., Tang, C., Wang, W., Wu, Q., Yin, Z., Zhu, L. Structural insights into mechanisms of rapid harvesting of microalgae with pH regulation by magnetic chitosan composites: A study based on E-DLVO model and component fluorescence analysis. Chemical Engineering Journal 2023, Volume 456, 141071. [CrossRef]
De Morais, E.G. , Sampaio, I.C.F., Gonzalez-Flo, E., Ferrer, I., Uggetti, E., García, J. Microalgae harvesting for wastewater treatment and resources recovery: A review. New Biotechnology, 2023, Volume 78, Pages 84-94. [CrossRef]
Rossi, S. , Visigalli, S., Cascino, F.C., Mantovani, M., Mezzanotte, V., Parati, K., Canziani, R., Turolla, A., Ficara, E. Metal-based flocculation to harvest microalgae: A look beyond separation efficiency. Science of The Total Environment, 2021, Volume 799, 149395. [CrossRef]
Figueiredo, D. , Ferreira, A., Quelhas, P, Schulze, P.S.C., Gouveia, L. Nannochloropsis oceanica harvested using electrocoagulation with alternative electrodes – An innovative approach on potential biomass applications. Bioresource Technology, 2022, Volume 344, Part B, 126222. [CrossRef]
Coward, T. , Lee, J.G.M., Caldwell, G.S. Harvesting microalgae by CTAB-aided foam flotation increases lipid recovery and improves fatty acid methyl ester characteristics. Biomass and Bioenergy, 2014, Volume 67, Pages 354-362. [CrossRef]
Fasaei, F. , Bitter, J.H., Slegers, P.M., van Boxtel, A.J.B. Techno-economic evaluation of microalgae harvesting and dewatering systems. Algal Res. 31, 2018, 347-362. [CrossRef]
Gerardo, M.L. , Oatley-Radcliffe, D.L, Lovitt, R.W. Minimizing the Energy Requirement of Dewatering Scenedesmus sp. by Microfiltration: Performance, Costs, and Feasibility. Environmental Science & Technology, 2014; 48, 845–853. [Google Scholar] [CrossRef]
Zhang, M. , Yao, L., Maleki, E., Liao, B.Q., Lin, H. Membrane technologies for microalgae cultivation and dewatering: Recent progress and challenges. Algal Res., 2022 44, 101686. [CrossRef]
Wang, L. , Pan, B., Gao, Y., Li, C., Ye, J., Yang, L., Chen, Y., Hu, Q., Zhang, X. Efficient membrane microalgae harvesting: Pilot-scale performance and techno-economic analysis. Journal of Cleaner Production, 2019, Volume 218, Pages 83-95. [CrossRef]
Hwang, T. , Park, S-J., Oh, Y-K., Rashid, N., Han, J-I. Harvesting of Chlorella sp. KR-1 using a cross-flow membrane filtration system equipped with an anti-fouling membrane. Bioresource Technology, 2013, Volume 139, Pages 379-382. [CrossRef]
Elcik, H. , Cakmakci, M., Ozkaya, B. The fouling effects of microalgae cells on crossflow membrane filtration. Journal of Membrane Science, 2016, Volume 499, Pages 116-125. [CrossRef]
Monte, J. , Sá, M., Galinha, C.F., Costa, L., Hoekstra, H., Brazinha, C., Crespo, J.G. Harvesting of Dunaliella salina by membrane filtration at pilot scale. Separation and Purification Technology, 2018, Volume 190, 2018. Pages 252-260. [CrossRef]
Mo, W. , Soh, L., Werber, J.R., Elimelech, M., Zimmerman, J.B. Application of membrane dewatering for algal biofuel. Algal Research, 2015, Volume 11, Pages 1-12. [CrossRef]
Bilad, M.R. , Discart, V., Vandamme, D., Foubert, I., Muylaert, K., Vankelecom, I.F.J. Harvesting microalgae biomass using a magnetically induced membrane vibration (MMV) system: Filtration performance and energy consumption. Bioresource Technology, 2013, Volume 138, Pages 329-338. [CrossRef]
De Baerdemaeker, T. , Lemmens, B., Dotremont, C., Fret, J., Roef, L., Goiris, K., Diels, L. Benchmark study on algae harvesting with backwashable submerged flat panel membranes. Bioresource Technology, 2013, Volume 129, Pages 582-591. [CrossRef]
Bhave, R. , Kuritz, T., Powell, L., Adcock, D. Membrane-Based Energy Efficient Dewatering of Microalgae in Biofuels Production and Recovery of Value Added Co-Products. Environ. Sci. Technol. 2012, 46, 10–5599. [Google Scholar] [CrossRef] [PubMed]
González-Camejo, J. , Andreola, C., Maceratesi, V., Toscano, G., Eusebi, A.L., Fatone, F., 2023. Biorefineries to improve water and resource recovery in the seafood-processing industry. Advanced Technologies in Wastewater Treatment: Food Processing Industry, pp. 127 - 154. [CrossRef]
Segredo-Morales, E. , González-Martín, C., Vera, L., González, E., 2023. Performance of a novel rotating membrane photobioreactor based on indigenous microalgae-bacteria consortia for wastewater reclamation. Journal of Industrial and Engineering Chemistry, Volume 119, Pages 586-597. [CrossRef]
Zhang, X. , Hu, Q., Sommerfeld, M., Puruhito, E., Chen, Y. Harvesting algal biomass for biofuels using ultrafiltration membranes. Bioresource Technology, 2010, Volume 101, Issue 14, Pages 5297-5304. [CrossRef]
Seco, A.; Aparicio, S.; González-Camejo, J.; Jiménez-Benítez, A.; Mateo, O.; Mora-Sánchez, J.F.; Noriega-Hevia, G.; Sanchis-Perucho, P.; Serna-García, R.; Zamorano-López, N.; et al. Resource recovery from sulphate-rich sewage through an innovative anaerobic-based water resource recovery facility (WRRF). Water Sci. Technol. 2018, 78, 1925–1936. [Google Scholar] [CrossRef] [PubMed]
González-Camejo, J. , Aparicio, S., Jiménez-Benítez, A., M. Pachés, M.V. Ruano, L. Borrás, R. Barat, A. Seco. Improving membrane photobioreactor performance by reducing light path: Operating conditions and key performance indicators. Water Research, 2020, Volume 172, 115518. [CrossRef]
Aparicio, S. , González-Camejo, J., Seco, S., Borrás, L., Robles, A., Ferrer, J. Integrated microalgae-bacteria modelling: Application to an outdoor membrane photobioreactor (MPBR). Science of The Total Environment, 2023, Volume 884, 163669. [CrossRef]
Mora-Sánchez, J.F.; Serna-García, R.; Bouzas, A.; Seco, A.; Ruano, M.V. Anaerobic Membrane Bioreactor for Microalgae and Primary Sludge Co-Digestion at Pilot Scale: Instrumentation, Control and Automation Implementation, and Performance Assessment. Water 2023, 15, 3225. [Google Scholar] [CrossRef]
APHA. American Public Health Association, American Water Works Association, Water Environment Federation. Standard methods for the examination of water and wastewater, 2012, 22nd., Washington, USA.
M Cheryan. Ultrafiltration and Microfiltration: Second edition. CRC Press LLC, 1998.
González-Camejo, J. , Barat, R., Aguado, D., Ferrer, J., 2020. Continuous 3-year outdoor operation of a flat-panel membrane photobioreactor to treat effluent from an anaerobic membrane bioreactor. Water Research, Volume 169, 115238. [CrossRef]
Rossi, S. , Sforza, E., Pastore, M., Bellucci, M., Casagli, F., Marazzi, F., Ficara, E. Photo-respirometry to shed light on microalgae-bacteria consortia—A review. Rev Environ Sci Biotechnol, 2020; 19, 43–72. [Google Scholar] [CrossRef]
Judd, S. , Judd, C., 2006. The MBR Book. Chapter 3 - Design. Elsevier Science. Pages 123-162. ISBN 9781856174817. [CrossRef]
Ortiz-Tena, F. , Ranglová, K., Kubač, D., Steinweg, C., Thomson, C., Masojidek, J., Posten, C., 2021. Characterization of an aerated submerged hollow fiber ultrafiltration device for efficient microalgae harvesting. Eng Life Sci.; 21: 607–622. [CrossRef]
Soulies, A. , Pruvost, J., Legrand, J.M., Castelain, C., Burghelea, T. Rheological properties of suspensions of the green microalga Chlorella vulgaris at various volume fractions. Rheologica Acta, 2013, 52 (6), pp.589-605. [CrossRef]
Kestin, J. , Sokolov, M., Wakeham, W.A. Viscosity of liquid water in the range −8 °C to 150 °C. J. Phys. Chem. Ref. Data, 1978; 7, 941–948. [Google Scholar] [CrossRef]
Ministry of Industry and Tourism; Net price of electricity for domestic use and industrial use: Madrid, Spain, 2023. Available online (in Spanish): https://www.mincotur.gob.es/es-es/IndicadoresyEstadisticas/BoletinEstadistico/Energ%C3%ADa%20y%20emisiones/4_12.pdf. (accessed on 23 October 2023).
Bunzl; Floquet Aqua: Sant Boi Llobregat, Spain, 2023. Available online (in Spanish): https://www.bunzlspain.com/hipoclorito-sodico-floquet-aqua-150g.html. (accessed on 23 October 2023).
Parravicini, V.; Nielsen, P.H.; Thornberg, D.; Pistocchi, A. , 2022. Evaluation of greenhouse gas emissions from the European urban wastewater sector, and options for their reduction. Sci. Total Environ., 838, 156322. [CrossRef]
Eriksson, O. Environmental Technology Assessment of Natural Gas Compared to Biogas. In Natural Gas; Potocnik, P., Ed.; InTechOpen: London, UK, 2010; pp. 127–146. [Google Scholar] [CrossRef]
Greses, S. , Degradación anaerobia de microalgas procedentes del tratamiento del efluente de un reactor anaerobio de membranas sumergidas (Anaerobic degradation of microalgae from the effluent treatment of a submerged anaerobic membrane reactor), PhD Thesis, Universitat de València, Valencia (Spain), 2017. (In Spanish). [Google Scholar]
Bremond, U. , Bertrandias, A., de Buyer, R., Latrille, E., Jimenez, J, Escudié, R., Steyer, J.-P., Bernet, N., Carrère, H., 2021. Recirculation of solid digestate to enhance energy efficiency of biogas plants: Strategies, conditions and impacts. Energy Conversion and Management, 231, pp.113759. [CrossRef]
Muller, O.; Bricker, P.; Laurin, B.; Mouton, L. Climate Change and Electricity, European Carbon Factor, Benchmarking of CO2 Emissions by Europe’s Largest Electricity Utilities, 21st ed.; INIS-FR--22-0804; International Atomic Energy Agency: Vienna, Austria, 2022. [Google Scholar]
Bright Biomethane, 2023. PUREPAC GRAND. Enschede, the Netherlands. Retrieved June 28th, 2023, from https://www.brightbiomethane.com/.
Mendez, L. , Garcia D. , Perez E., Blanco S., Munoz R. Photosynthetic upgrading of biogas from anaerobic digestion of mixed sludge in an outdoors algal-bacterial photobioreactor at pilot scale. Journal of Water Process Engineering, 2022, 48, art no 102891. [Google Scholar] [CrossRef]
Du, X.; Shi, Y.; Jegatheesan, V.; Haq, I.U. A Review on the Mechanism, Impacts and Control Methods of Membrane Fouling in MBR System. Membranes 2020, 10, 24. [Google Scholar] [CrossRef] [PubMed]

Figure 1. Flow diagram of the cross-flow ultrafiltration-based harvesting system.

Figure 2. Evolution of filtration performance regarding: (a) TSS concentration in the retentate and biomass to be harvested (TSS_i·V_i); (b) 20 °C- standardised transmembrane flux (J₂₀) normalised by the initial flux (J_20,0); concentration ratio (r) and biomass to be harvested (TSS_i·V_i)/10. Dashed lines indicate the performance of chemical cleanings.

Figure 3. Details of the filtration performance for an aleatory period of 24 days of operation (cycles 19 to 33): (a) TSS concentration in the retentate and biomass to be harvested (TSSi·Vi); (b) 20 °C- standardised transmembrane flux (J₂₀) normalised by the initial flux (J_20,0); concentration ratio (r) and biomass to be harvested (TSS_i·V_i)/10. Dashed lines indicate the performance of chemical cleanings.

Figure 4. Energy balance and GHG assessment for a concentration ratio of 21 for three different biogas valorisation scenarios: Scenario 1: co-generation; Scenario 2: membrane upgrading; Scenario 3: microalgae upgrading: (a) energy balance in terms of Q_demand, W_demand, and Q_BM; (b) direct and indirect GHG emissions.

Figure 5. Total GHG emissions for Scenario 3 (microalgae upgrading to biomethane).

Table 1. Pros and cons of microalgae biomass separation technology.

Technology	Advantages	Disadvantages	References
Sedimentation	Simple. Low CAPEX and OPEX.	Poor settling capacity. Low effluent quality. Low biomass concentration. Biomass losses. Time consuming.	[24,34,37]
Coagulation- Flocculation	Proper settling rate. Proper effluent quality. Effective as pre-concentrating step.	High doses of chemical reagents. Possible photosynthesis inhibition from metallic flocculants. Metal presence in harvested biomass. Poor effluent quality for reuse.	[28,34,37−40]
Electroflocculation	Flocculants not required.	Metal presence in harvested biomass. Emerging technology (low TRL).	[25,28,41]
Magnetic flocculation	Fast, scalable, and efficient. Natural polymers can be added to coat magnetic particles.	Low efficiency for small particles. Low magnetic capacity lifetime. Emerging technology (low TRL).	[25,38,39]
Bioflocculation	Toxic chemicals are not required. Prevention of microalgae contamination. Bioflocculants are produced from biomass. Genetic engineering can enhance bioflocculants production.	Factors affecting bioflocculant release remain unclear. Emerging technology (low TRL).	[25,39]
Flotation	Low CAPEX. Low HRT. High cost-effective method. Small footprint. Possible disruption of microalgae (pro for anaerobic valorisation).	Use of reagents can imply extra cost. Possible disruption of microalgae (con if integrity of biomass is needed).	[24,39,42]
Centrifugation	Quick and simple. Most algal cell types can be harvested. Can be used as second-step harvesting process.	High energy demand. High CAPEX. Shear stress to microalgae. Low EPS removal.	[25,35,36,39,43]
Filtration	High effluent quality. Small footprint. Can be combined with pre-harvesting steps. Moderately to high biomass concentration. Depending on the technology and configuration employed, can ensure the integrity of the microalgae or promote the hydrolysis of microalgae.	Membrane fouling and clogging. Membrane productivity in terms of CAPEX and OPEX. High energy requirements.	[21,25,26,44−49]

Table 2. Average operating and outdoor conditions over the 78 cycles of cross-flow ultrafiltration.

Parameter	Unit	Mean ± SD
Temperature	ºC	22 ± 6
Inlet gauge pressure (P1)	bar	2.11 ± 0.04
Outlet gauge pressure (P2)	bar	0.48 ± 0.02
Transmembrane pressure (TMP)	bar	1.25 ± 0.02
Initial volume in the feed tank (V_i)	L	450 ± 65
Initial TSS in the feed tank (TSS_i)	g·L⁻¹	0.50 ± 0.12

Table 3. Average filtration performance over the entire operating period (78 cycles).

Parameter	Unit	Mean ± SD
TSS_f	g·L⁻¹	9.7 ±1.7
CFV	m·s⁻¹	1.2 ± 0.3
J₂₀	LMH	16 ± 8
Normalised flux (J_20:J_20,0)	-	0.49 ± 0.23
K₂₀	LMH·bar⁻¹	13 ± 6
TMP_BF	bar	1.21 ± 0.08
Q_BF	L·min⁻¹	31 ± 8
HV_r	g·m⁻²·h⁻¹	9 ± 3
r	-	21 ± 6
Cycle duration	h	25 ± 9

CFV: cross flow velocity; HV_r: biomass harvesting rate; J₂₀: transmembrane flux standardised at 20 °C; J_20,0: initial J₂₀; K₂₀: membrane permeability standardised at 20 ºC; Q_BF: backflush flow-rate; r: concentration ratio; TMP_BF: backflush stage transmembrane pressure; TSS_f: total suspended solids of the harvested biomass.

Table 4. Energy and chemical reagents consumption of the cross-flow ultrafiltration unit.

Parameter	Acronym	Unit	Mean ± SD
Energy consumption ratio per harvested algal biomass	EC_m_TSS	kWh·kg⁻¹	1.51 ±0.64
Energy consumption ratio per treated volume of pre-concentrated microalgae culture	EC_v_HV	kWh·m⁻³	0.76 ± 0.32
Energy consumption ratio per total treated volume in the WRRF	EC_v_WRRF	kWh·m⁻³	0.39 ± 0.16
Sodium hypochlorite consumption per total treated volume in the WRRF	NaOClC_V_WRRF	g Cl·m⁻³	0.97 ± 0.39
Operating expenses of energy and sodium hypochlorite per total treated volume in the WRRF	OPEX_EC+Cl	€·m⁻³	0.082 ± 0.034

Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

© 2023 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (http://creativecommons.org/licenses/by/4.0/).

Copyright: This open access article is published under a Creative Commons CC BY 4.0 license, which permit the free download, distribution, and reuse, provided that the author and preprint are cited in any reuse.

MDPI Initiatives

Important Links

Choose an area of interest and we will send you notifications of new preprints at your preferred frequency.

Disclaimer