Hydrothermal co-liquefaction of microalgae, wood, and sugar beet pulp

Hydrothermal co-liquefaction of mixed (wet and dry) biomass residue streams would greatly enhance the viability and scale up potential of the technology as platform in bioenergy and biorefinery applications. This study aims to identify possible interaction effects between three different feeds (protein-rich microalgae, lignocellulosic wood, and carbohydrate-rich sugar beet pulp) and to broaden the data set for evaluating this concept. Co-liquefaction was evaluated at 250 and 350 °C at 10 min of holding time, using 10 wt%( in water) binary mixtures (1:1 wt basis) and a (1:1:1 wt basis) ternary mixture. Results show that interaction during co-liquefaction does play a role and especially reduced the amount of biocrude produced. The biocrude yields obtained are around 15 and 40% below the estimated values for binary and ternary mixtures, on basis of linear averaging the results for the single feeds. For mixtures including algal biomass, a more than proportional nitrogen content and fraction of high molecular mass components was found in the biocrude. For the predictability of biocrude yield and composition in case of biomass mixtures, more work is needed to unravel these interactions.


Introduction
Hydrothermal processing of microalgae received in recent years a lot of attention, and not without reason; microalgae are one of the most promising bioresources for sustainable production of food and feed ingredients, high valuable chemicals, biofertilizer, and/or fuel components. The relatively high costs of cultivation can only be justified if the complete algal biomass production is valorized, either by whole-cell selling of microalgae or by producing multiple products for different markets, a so-called Algae Biorefinery. Due to the aqueous nature of the (algae) feed, hydrothermal processing is a natural choice as downstream technology for any residual algal biomass into energy carriers (biocrude and/or gaseous).
Hydrothermal liquefaction is a process in which aqueous slurries of biomass are brought to medium-high temperatures (200-375°C) at elevated pressures (up to 25 MPa) to keep the water in the liquid state. Under these conditions, the biomass converts in this hot, compressed water into gaseous compounds, dissolved organics, an organic oil phase (called biocrude), and a residual solid fraction. The biocrude is the targeted liquid phase energy carrier. The dissolved organics can be used for heterotrophic growth of aquatic biomass or gasified in hot compressed water at even higher temperatures to recover energy and nutrients. The solids produced normally contain most of the mineral matter, coke, and unconverted biomass. The gas produced consists merely of CO 2 , and recycling thereof to the algae production stage seems to be a reasonable option.
With the worldwide increase of interest in algae cultivation and processing, new attention for hydrothermal processing arose [1]. Processing whole microalgae slurries showed that high oil yields and energy recoveries (even up to 75% for lowlipid algae) are attainable and detailed analyses of the biocrude compounds produced were carried out [2,3]. The biocrude produced is not directly suitable as transportation fuel but can be considered as renewable feed for co-feeding in existing refineries.
However, hydrothermal liquefaction based on microalgae alone will not make a significant impact on the fossil crude replacement in next decade, based on the current production levels of microalgae. Therefore, it will not be easy for the hydrothermal liquefaction of microalgae to achieve an economy of scale, which is required for successful entry in a bulk market as transportation fuels.
For most other biomass resources, like forest thinning, agricultural residues, manure-derived biomass, road-side grass and straw, bagasse etc., the amount of biomass available is more abundant, but its quantity and quality are fluctuating in time and place. Therefore, for the development of hydrothermal liquefaction as a viable platform technology, it would be of great advantage if microalgal biomass could be combined with different biomass sources to create a larger scale of operation and more continuity in operation (for seasonal bound resources).
Liquefaction of different biomasses in a single step process seems an attractive and simple solution, but understanding of the fundamentals of hydrothermal liquefaction and the effects of a changing feed composition is only starting to develop. The effect of microalgal composition was studied by Biller and Ross [4] and recently by Leow et al. [5]. Biller and Ross identified that in decreasing order, lipids > proteins > carbohydrates are converted to the biocrude product. They found that proteins and lipids were converted to biocrude efficiently without additional catalyst, whereas the biocrude yield from carbohydrates increases when using sodium carbonate as catalyst [4]. Leow et al. developed an improved component additivity model to predict biocrude yield on basis of lipid, protein, and carbohydrate content. In their approach, they used a single strain of microalgae, cultivated under varying conditions to vary its biochemical composition [5].
The liquefaction of microalgae was compared with that of other biomass sources by a few researchers. Vardon et al. studied the liquefaction of Spirulina and compared this with swine manure and anaerobic digested sludge [6]. They found a significant effect of the feed composition on the product distribution and biocrude composition. From their work, it was concluded that information on feed composition and a better understanding of the chemistry during liquefaction are needed to be able to predict biocrude yield and composition. Physical mixtures of these feeds were, however, not studied.
For wood-type biomass, the crude yield from liquefaction in water was found to be strongly related with the type and amount of lignin in the biomass [7]. The crude yield decreased linearly with the lignin content, whereas char yield increased linearly. When adding alkali salts, the liquefaction crude yields increased for most woody biomass types studied, but the linearity with lignin content disappeared. Remarkably, there was little effect on char yield. These differences were attributed to the different types of lignin present, and the results show once more that predicting hydrothermal liquefaction yields is not trivial.
Co-liquefaction of wood with glycerol was demonstrated in a continuous liquefaction setup by Pedersen et al. [8]. Their product analysis did not reveal to what extent glycerol was converted during processing and contributed to the biocrude. Chen et al. studied the co-liquefaction of swine manure with mixed algal biomass from waste water streams at 300°C with 1 h of reaction time [9]. With increasing swine manure content, the char yield decreased dramatically from around 60 to 17 wt%, while the aqueous phase products showed a very distinct optimum at a 50/50 w/w feed ratio of manure/algal biomass. The biocrude product yield for a 50/50 w/w mixture could not be predicted from the experiments with algal mass and with swine manure alone. A similar result was found in the work by Xiu et al. [10] on the co-liquefaction of swine manure with vegetable oil.
In this work, new experimental results are presented for hydrothermal co-liquefaction of different biomass resources, alone and in mixtures. In this study, three distinctly different biomass types are subjected to liquefaction experiments. As biomass sources, pine wood (W), sugar extracted sugar beet pulp (S), and microalgae (A) were used. The pine wood was chosen to represent wood-type ligno-cellulosic biomass, and the (sugar extracted) sugar beet pulp is taken as carbohydrate-rich biomass. The experiments at different reaction time, different temperature, and with different feed composition are carried out to see to what extent the liquefaction of one feedstock influences the results of the other fractions in the feed mixture and to what extent the biocrude yield and composition can be predicted on basis of the results for the separate feed constituents.
This work is the start of a broader set of studies with mixed feedstocks, aiming to be able to predict the liquefaction products from mixed feeds and to optimize conversion conditions. The experimental results in this study should help to pave the way for hydrothermal liquefaction as single step process option for residue streams to produce a renewable, liquid energy carrier (biocrude).

Materials used
As microalgae source, the freshwater species Desmodesmus sp. was used in this work. The microalgae were obtained as dry powder from a commercial source and identical to the feed used in earlier work [2,3]. The pine wood was obtained from Rettenmaier & Söhne GmbH (Germany) and milled three times and sieved to a sieve fraction between 0.2 and 0.6 mm. The composition of the pine wood used for this research is listed in Table 1. The sugar beet pulp used was purchased from Duynie B. V (The Netherlands), where sugar beet was cut and sugar was extracted. The remaining substance was pressed, dried, and finally pressed into pellets of 6 mm. For this research, the sugar beet pulp was milled and separated by different particle size. The particle size <0.3 mm was used. The C, N, and H contents were measured by elemental analysis (EA) using a Thermo Scientific Flash 2000 apparatus for CHN-S determination. The analysis of sugar beet pulp is listed in Table 1.

Methods
HTL experiments were carried out using a 9-mL stainless steel autoclave, consisting of a cylinder with top and bottom openings. The 9-mL reactor was designed with very thin walls to allow for a fast heating rate (∼300°C/min). The autoclave can be rapidly heated up by immersing it in a hot fluidized sand bed and, after reaction, rapidly quenched again by submerging in cold water. The complete setup was located in an explosion-proof high-pressure safety room and essentially the same as used by Garcia Alba et al. [2], except for the reactor size.
The experimental procedure for liquefaction experiments starts by loading into the autoclave a pre-weighted amount of feedstock and water, typically around 0.400 g of biomass and 3.60 g water, resulting in 10 wt% biomass in water slurry. After loading the sample, the autoclave was tightly closed and connected to the pneumatic shaking piston. The gate valve from the reactor was tightly connected with a gas pipe. The autoclave was flushed three times with 1 MPa of nitrogen and pressurized to about 1 MPa. The experiment was controlled from outside the safety box via a control panel. The rotating shaking arm was turned on, and the autoclave was immersed into the top of fluidized sand bed. The moment of immersion in the sand bed marks the start of the reaction time. The heating time to the desired reaction temperature was within 1-2 min. After the desired reaction time, the reactor was lifted and quenched in a cold water bath; each time cooling was around 5 min. After the reactor was cooled, the pressure inside was noted, and the gas in the autoclave was released carefully, to not loose liquid because of dissolved gas. Then, the product was collected after disconnecting and opening the top part of the reactor. Afterward, to collect the maximum amount of product, the autoclave was filled three times of 5 mL of dichloromethane (DCM, Sigma-Aldrich 99.8% purity), closed, connected to the pneumatics shaker, and shaken in water bath (5 min each time). Then, the reactor was rinsed once more using 10 mL of solvent. Next, other parts of the setup (top of the reactor and gas pipe) were flushed with DCM to collect any remaining product.
The sample mixture was collected in a glass bottle and was vacuum filtrated over a 1-μm pore size filter (Whatman GF/B). The glass bottle used to collect the sample from the reactor was washed three times with 5 mL of DCM. Next, the filter cake was also rinsed with 5 mL of dichloromethane, ensuring that the filter consisted of just solid residue. The solids were dried in the oven at 105°C for 24 h and weighed. The filtrate was collected in a glass bottle and then was centrifuged at 9000 rpm for 10 min to optimize the separation between the dichloromethane and water phase. The aqueous top phase was separated by aspiration using a syringe.
The water-soluble organics were quantified by evaporating the water phase at 70°C for 24 h, obtaining mixture of organics and ash soluble in water. This mixture was put into the oven at 550°C for 5 h to quantify the ash content. The watersoluble organics were calculated by difference. For the organic dichloromethane fraction, the DCM was removed by evaporation by using vacuum oven with a cold trap system to collect the oil.
After the whole procedure, four products were obtained: biocrude oil (obtained from DCM-soluble fraction), gas, watersoluble organics, and solid residue. All steps are shown in Fig. 1.
Despite difficulties related to the small sample sizes used, the mass balance closure was in most cases on average 93 wt% and varied from 82 to 98 wt%. A slight tendency was observed that mass balance closure was better at higher solid content, but no trends with temperature nor feed were found.

Analysis
The carbon, hydrogen, and nitrogen contents of biocrude products were determined with an elemental analyzer. The oxygen The energy recovery, from biomass to biocrude, is calculated using the following: Gel permeation chromatography (GPC) is used to determine the molecular weight distribution of the bio-oil obtained after HTL. The bio-oil was dissolved in THF, filtered through 0.2 μm syringe filter, and injected into the system. GPC analysis is performed using an Agilent 1200 series HPLC system with 3 GPC PLgel 3 μm MIXED-E columns connected in series. The column temperature was 40°C, with a flow of 1 mL/min, with tetrahydrofuran (THF) being the solvent.
The biocrude was divided into the fractions referred to as Bheavy^and Blight^compounds. The light fraction is defined as all of the components in the biocrude with a MW GPC lower than 200 Da and a heavy fraction, comprising compounds with a MW GPC higher than 1000 Da. The RID detector signal is plotted versus the elution time, and the area corresponding to the heavy or light fractions is divided by the total area.
Functional group compositional analysis was carried out in mid-infrared wavenumber range using FTIR spectrometry (Bruker), equipped with an attenuated total reflection system (ATR).

Product yields
For each of the biomass types, separate liquefaction experiments were performed, prior to mixing the biomass types. All experiments were measured in duplo or triplo, and the results presented in Fig. 2 are the mean values with standard deviation. For all experiments, the biomass content of the slurry is 10 wt%.
The results from Fig. 2 show already some clear differences between the different feeds. Wood produced the largest amount of solids and a relatively low amount of biocrude. Especially, the solids collected for the 250°C run visually resembled the (unconverted) wood feed material. Microalgae clearly produced the largest amount of biocrude. It is also remarkable to see that some of the trends are different with respect to an increase in temperature. For all feeds, the biocrude and gas yield increase when going from 250 to 350°C. However, the aqueous phase organics (AP) decrease for algae and beet pulp, but increase for wood. Whereas the solid residue increased for sugar beet pulp, it strongly decreased for wood and algae as feed.
It seems that with increasing temperature for microalgae, the solids and aqueous phase organics are more converted into additional biocrude, whereas for beet pulp, more solids are formed probably because of secondary polymerization of the carbohydrate-derived decomposition products. For wood, the AP yield increases with temperature, which is probably related to an incomplete primary conversion of the biomass feed at 250°C. The results show that although sugar beet seems more easily converted (at lower temperature) than microalgae, the selectivity towards biocrude of the latter is much more favorable. Apparently, wood is most recalcitrant towards liquefaction, most likely related to its lignin content. The attractiveness of microalgae for biocrude production via HTL is clearly subscribed by the product distribution in Fig. 2, where the specific biocrude yield for microalgae at 350°C is more than double the values obtained for pine wood and sugar beet pulp.
To evaluate possible interactions during hydrothermal liquefaction, 50/50 wt%/wt% mixtures of different combinations of the three feeds were made and used for liquefaction experiments. The results are presented in Fig. 3. The HTL product yields for these mixtures are compared in Fig. 3 with the weight-averaged yields based on the results of the separate feeds, see Eq. (3), as presented in Fig. 2.
In Eq. (3), Y p,mix is the estimated product yield (weight basis) for the mixture, with P = gas, aqueous phase organics, biocrude, or solids, and w feed,i is the weight fraction of that feed type in the total feed mixture. Y p,i is the corresponding product yield (for p = gas/biocrude/AP/solids) as obtained for feedstock i (here: i = microalgae, sugar beet pulp, or pine wood). To correct for the (small) differences in the mass balance closure (MB) for the different experiments, the This Blinear averaging^model is found to be reasonable successful in predicting the product yields at 250°C, but less successful in predicting the product yields for the 350°C series. For the biocrude yield, which is the targeted product, the model significantly overpredicts the yield with a relative error of around 15%, calculated via (Y p,mix − Y p,exp. ) / Y p,exp. , with Y p,exp. the experimental value.
In Table 2, the average deviations (for each product class, for all binary mixtures) between model predictions and actual yields are presented. From Table 2, it is clear that the prediction is less accurate at the higher temperature. The underpredicted AP yield and the overpredicted solid yield for the algae-sugar beet pulp combination (Fig. 3a) suggest an enhanced biomass decomposition in the presence of microalgae liquefaction products. For the pine wood-algae mixture, less biocrude was found than predicted at both temperatures, while the AP organics and gaseous product yields were higher than estimated. This suggests that the algae liquefaction products assist in breaking down the fragments produced during liquefaction of other biomasses.
Remarkably, the observed decrease in solid residue in Fig. 5a is not accompanied by an (expected) increase in biocrude yield. On the contrary, instead of a higher than expected (by linear averaging) biocrude yield, the biocrude yield is also overpredicted in these cases. The product class which benefits most from the reduced solids and biocrude yield is the aqueous phase organics.
For the mixture of the three feeds together (microalgae, pine wood, and sugar beet pulp; each one third on mass basis), the experimental results as well as the predicted yields by the linear model are presented in Fig. 4. Next to linear averaging of the results for the separate feeds, also a linear averaging of the results for the three binary mixtures was tested, but this did not significantly improve the predictions. As the experimental procedures for the separation of biocrude from the solids might have played a role, the estimation for these grouped classes is reported in Fig. 4b. It is found that for the 350°C data, the higher than predicted yield of the fraction AP + gas is visible, beyond possible effects of incomplete solids/biocrude separation.
Most striking result in Fig. 4a is the strong overprediction (39% at 250°C and 47% at 350°C) of the crude oil yield. Since the mass balance closure of these experiments with the ternary mixtures was at the lower end of all experiments (90 wt% versus overall 94 wt% ± 4% standard deviation), one is tempted to suggest that the incomplete mass balance closure is merely due to the crude fraction. There is, however, no reason to assume that in these mixed feed experiments, the error in the crude oil yield determination is larger than in the other experiments. Also, the reproducibility of these experiments was very good. A possible explanation for the lower mass balance closure could be that the mixture leads to the production of more light and volatile products in the aqueous phase (AP) and/or in the biocrude that are lost during solvent evaporation, but this would not affect the conclusion with respect to the biocrude yield.
The energy recovery, in the form of the targeted biocrude yield and its heating value, is for all mixtures less than expected by averaging the results of the single feeds, see Fig. 5. This effect is most pronounced for the ternary biomass mixed feed. It can be concluded that there is no synergetic effect observed for the energy recovery for mixing these biomass streams. The observed anti-synergetic effect underlines the importance of testing combinations of biomass feeds and improving our

Product analysis
The biocrudes produced by the hydrothermal liquefaction at different temperatures for the different feeds were subjected to GPC analysis, elemental analysis, and FT-IR spectroscopy. In Fig. 6, next to the crude yields, the nitrogen content and the higher heating values are reported in a comprehensive   overview. It is clear that a high crude yield and high HHV are reached using algae as (part of the) feed. At these conditions, also the nitrogen level in the crude is relatively high. It furthermore shows that the application of microalgae in the feed leads to a relatively high content of Bheavies,^defined as GPC detectables with an apparent molar weight exceeding 1000 Da.
In Fig. 7, some characteristics of the biocrudes produced are compared with those for the biocrudes produced for the separate biomass feeds used. In Fig. 7a, the nitrogen content of the biocrude is presented. It is clear to see that the N content of the actual biocrudes exceeds the estimated value based on linear averaging of the contribution of the separate feeds for mixtures  in which algal biomass is present. This suggests that the algalbased contribution to the biocrude is more than actual fraction of the algal biomass in the feed mixture. Especially when considering that the biocrude yields are (on average) some 15% less than estimated by the same linear averaging method.
In Fig. 7b, for the fraction heavies in the biocrude, this effect is even stronger. Biocrude from algal biomass containing mixtures has around 40% more heavies than estimated via linear averaging of the biocrudes for the separate feeds. This confirms that the algal-based contribution to the biocrude is more than the actual  fraction of the algal biomass in the feed mixture. It shows that (part of) the biocrude constituents in biocrude from algal biomass are larger than obtained from the other feeds tested, and considering the similarity between Fig. 7a, b, this heavies fraction has a significant N-content (protein-derived compounds). For the wood-sugar beet pulp mixture, the heavies fraction and Ncontent are very similar to the estimated values at both temperatures from the separate crudes, suggesting little interaction. Figure 8a, b shows that the GPC profile of the crude obtained for the ternary mixture (all feeds) resembles most the one obtained for the crude obtained from the micro algae feed alone. This holds for the both the 250°C HTL crudes as well as for the 350°C HTL crudes. Similar findings were observed in the GPC analysis of the binary mixtures. In mixtures in which algae were used as part of the feed, the GPC profile of the biocrude from the mixture is very similar to that of crude obtained by HTL of microalgae under the same conditions. For the sugar beet-pine wood mixture, hence without the microalgae, the profile was more an Baverage^of both. Comparing the FT-IR signals showed a more Bmixed response^of the crudes produced from mixed feedstocks, in comparison with the FT-IR signals for the crudes based on the separate feeds.

Discussion
The results for the product yield show that for the different feeds, clearly different results are obtained with respect to product distribution (crude, solids, AP, and gas) and with respect to crude composition (as indicated by N content and heavies fraction). With increasing temperature under the HTL conditions studied (10-min reaction time and 10 wt% feed) in most cases, the solid yield for the feed mixtures is decreasing, suggesting that the solid feeds are more readily converted to smaller compounds, especially aqueous phase organics and gaseous compounds. From Fig. 7, it appears that the strongest interaction is for the feed mixtures containing algal biomass. The reason for this could be that microalgae (as it contains no recalcitrant lignin) decompose at lower temperatures than sugar beet pulp and pine wood. Hence, liquefaction products from the microalgae liquefaction may enhance liquefaction of the other biomasses by, e.g., catalyzing hydrolysis. Acetic acid is one of the compounds that may cause this effect. In an earlier study by Ramos-Tercero, the production of acetic acid from microalgae liquefaction is shown [14]. It was also shown that for a single biomass feed of microalgae, there was no correlation between acetic concentration and biocrude production. However, the hydrolysis of, e.g., compounds in the sugar beet pulp (Fig. 3a) to aqueous phase organics or the production of gaseous compounds from pine wood, may be enhanced by the higher acid concentration [15]. Considering that algal biomass-derived biocrude is more than proportionally represented in the biocrude of the mixture of algal biomass with pine wood (see Fig. 7a, b), it is even more surprising that the biocrude yield itself is less than based on linear averaging and illustrates the complexity of the interaction.
The results from Fig. 7 on the N-content and on heavies fraction, in combination with the reduced biocrude yield, clearly show that the co-liquefaction of algal biomass leads to enhanced decomposition of the other biomass components (pine wood, sugar beet pulp) tested into fragments that do not report to the DCM soluble biocrude fraction. More experimental work and more in-depth analysis of the chemical composition of biocrude and of the aqueous phase organics are required to identify more clearly the nature of the interaction.

Conclusions
Hydrothermal liquefaction of three different biomass types (pine wood, microalgae, and extracted sugar beet pulp) at a relatively low (250°C) and high (350°C) temperature showed clearly that interactions occur during co-liquefaction of different biomass feed types. Mixing feeds resulted in all cases in lower than expected biocrude yields (on average 15% reduction for the binary mixtures and 40% for the ternary mixture), when comparing with linear averaging of the results for the separate feeds. From the feeds studied, the use of microalgae gives the highest crude yields (almost double the amount for pine wood), with the largest fraction of heavies but also the highest nitrogen content in the oil. When used in mixtures, a more than proportional part of the biocrude formed appears algal biomass based. The significant deviations of the predictions in biocrude yield by the linear averaging model underline that more indepth studies are needed to unravel the interactions occurring during co-feeding different biomasses. Despite the lack of an accurate model, mixing different feeds seems a simple method to overcome the seasonality and the variations in availability and quality of biomass resources, thereby enabling hydrothermal liquefaction to become a biorefinery platform technology, able to benefit from Beconomy of scale^effects.
Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http:// creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made.