Outline
Comptes Rendus

Anaerobic co-digestion of dairy raw by-products and Ulva sp. macroalgae: effect of organic and inorganic additives
Comptes Rendus. Chimie, Sustainable Biomass Resources for Environmental, Agronomic, Biomaterials and Energy Applications 2, Volume 24 (2021) no. S1, pp. 23-37.

Abstract

An investigation has been carried out for biogas production from dairy raw materials (DRM) and Ulva sp. macroalgae as a co-substrate in order to find a usefulness for these species. Some nutrient media have been selected to optimize methane yield. The effect of the mineral additive medium (medium I), based on bicarbonate compound, and the organic additive one (medium II), based on glucose, on the anaerobic digestibility of DRM has been assessed by using a single-stage anaerobic process. Four batch experiments (batch 1, batch 2, batch 3 and batch 4) have been performed and correspond respectively to DRM without inoculum, DRM with additive medium I, DRM with additive medium II and DRM and Ulva sp. with additive medium I. Measurements of methane produced, chemical oxygen demand (COD), proteins, lipids, carbohydrates and volatile fatty acids have been carried out. It has been found that the mineral additive medium provided a higher specific methane yield compared to the organic additive one, the methane yield was about 0.208 L CH 4 /g VS . The co-digestion of DRM with Ulva sp. showed that the methane yield reached was about 0.118 L CH 4 /g VS , and revealed that the use of Ulva sp. as co-substrate reduces the concentration of volatile fatty acids, and enhances the biogas quality which contains 96% of methane and only 0.5 to 2% of hydrogen sulphide. Furthermore, a first-order kinetic model has been applied to describe proteins, lipids and carbohydrates biodegradation and showed that these three compounds follow a different kinetics according to the composition of the medium and the presence or not of additive medium.

Metadata
Online First:
Published online:
DOI: 10.5802/crchim.74
Mots-clés : Dairy raw materials, Additives, Ulva sp. macroalgae, Biogas, Volatile fatty acids

Selsabila Zara 1; Rachida Rihani 1; Walid Blel 2; Fatiha Bentahar 1

1 Laboratoire phénomènes de Transfert. Faculté de génie mécanique et génie des procédés. Université des Sciences et de la Technologie Houari Boumediene, U.S.T.H.B, BP 32 El-Alia, Bab-Ezzouar, Algérie
2 Université de Nantes, CNRS, GEPEA, UMR 6144 F-44600 Saint-Nazaire, France
License: CC-BY 4.0
Copyrights: The authors retain unrestricted copyrights and publishing rights
@article{CRCHIM_2021__24_S1_23_0,
     author = {Selsabila Zara and Rachida Rihani and Walid Blel and Fatiha Bentahar},
     title = {Anaerobic co-digestion of dairy raw by-products and {\protect\emph{Ulva}} sp. macroalgae: effect of organic and inorganic additives},
     journal = {Comptes Rendus. Chimie},
     pages = {23--37},
     publisher = {Acad\'emie des sciences, Paris},
     volume = {24},
     number = {S1},
     year = {2021},
     doi = {10.5802/crchim.74},
     language = {en},
}
TY  - JOUR
AU  - Selsabila Zara
AU  - Rachida Rihani
AU  - Walid Blel
AU  - Fatiha Bentahar
TI  - Anaerobic co-digestion of dairy raw by-products and Ulva sp. macroalgae: effect of organic and inorganic additives
JO  - Comptes Rendus. Chimie
PY  - 2021
SP  - 23
EP  - 37
VL  - 24
IS  - S1
PB  - Académie des sciences, Paris
DO  - 10.5802/crchim.74
LA  - en
ID  - CRCHIM_2021__24_S1_23_0
ER  - 
%0 Journal Article
%A Selsabila Zara
%A Rachida Rihani
%A Walid Blel
%A Fatiha Bentahar
%T Anaerobic co-digestion of dairy raw by-products and Ulva sp. macroalgae: effect of organic and inorganic additives
%J Comptes Rendus. Chimie
%D 2021
%P 23-37
%V 24
%N S1
%I Académie des sciences, Paris
%R 10.5802/crchim.74
%G en
%F CRCHIM_2021__24_S1_23_0
Selsabila Zara; Rachida Rihani; Walid Blel; Fatiha Bentahar. Anaerobic co-digestion of dairy raw by-products and Ulva sp. macroalgae: effect of organic and inorganic additives. Comptes Rendus. Chimie, Sustainable Biomass Resources for Environmental, Agronomic, Biomaterials and Energy Applications 2, Volume 24 (2021) no. S1, pp. 23-37. doi : 10.5802/crchim.74. https://comptes-rendus.academie-sciences.fr/chimie/articles/10.5802/crchim.74/

Version originale du texte intégral (Propose a translation )

1. Introduction

According to the International Energy Agency report entitled “key world energy statistics” natural gas provided about 20.3% of the total world energy consumption in 2014 [1]. The methane is the main component of natural gas and can be used as a fuel or as a feedstock in the chemical industry, depending on the quality of the processed gas. Biogas can be considered as a replacement for fossil gas fuels used in transport and power production, in order to handle environmental concerns [2, 3, 4]. It can be produced by the anaerobic digestion (AD) of organic matter contained in food processing wastes, lignocellulosic biomass, and animal manures, which contain carbohydrates, fats and proteins. Often, these feedstocks are available at a low cost, thus making the production of biogas economically more profitable [5].

The dairy industry was an important world economic driver for many years, and it remains so today, it accounts since almost $400 billion in 2015 [6]. In addition, among the food processing industries, the processing of milk consumes substantial amounts of water to ensure the quality control and hygiene standards, thus, large quantities of effluents are generated [7]. For example, the Algerian dairy industry generates nearly 47,000 t only for the whey, which can produce up to 1.41 million m3 of methane [8].

The physico-chemical characteristics of the related produced effluents are variable. They depend on the targeted dairy products and the used cleaning procedures [9]. For instance, the chemical oxygen demand (COD) contents of effluents from dairy industry vary regarding the applied processing. According to Matos et al. [5], the milk processing generate effluents with average COD value around 3000 mg⋅L−1. While, this value becomes much higher in the case of cheese production and can reach about 50,000 mg⋅L−1 [8]. Several dairy industries, particularly, those using cheese whey, discharge effluents with low alkalinity. The use of these substrates for biogas production, increases the volatile fatty acids production which could significantly inhibit the activity of methanogenic bacteria [10].

Demirel et al. [11] emphasize that the presence of certain compounds such as fats, suspended solids and nutrients (nitrogen and phosphorus) also contribute in the increase of the COD level. If dairy effluents are not properly disposed, they may result in significant environmental hazards [12]. Indeed, ineffective and inefficient treatment of these wastes may contribute to the eutrophication of waters. The latter can lead to serious water related quality problems such as loss of biodiversity and toxic effects of algal toxins on the whole ecosystems [13].

Among the available techniques for biological waste treatment, anaerobic digestion (AD) of dairy industries is suggested as a truly sustainable process which can handle the contained high organic contents. The effluents from dairy industries have a high organic matter and low alkalinity due to the by-products generated, several solutions have been advanced, such as the co-digestion by using biomasses with large buffering capacity or the addition of chemicals compounds [13, 14, 15, 16].

Macroalgae accumulation in the near-shore areas and their subsequent decomposition can constitute a particularly serious threat to the marine ecosystem [17]. Particular interest may be assigned to the macroalgae or microalgae biomasses as feedstocks to produce energy [18]. Indeed, macroalgae biomass can be used as potential materials for the biogas and bioethanol production, because of their high level of carbohydrates and low contents of lignins [3].

Ulva species, well known as sea lettuce is easily biodegradable and can provide several nutrients and minerals to sustain bacterial growth [19]. This species remains less valued than other algal biomass at industrial scale due to the seasonality of the seaweed availability and the high H2S content in the biogas produced by anaerobic digestion of this macroalgae [20, 21].

Co-digestion can improve AD performance and stability by monitoring the ratio of carbon/nutrient, overcoming deficiencies of some elements, improving the buffering effect, and reducing the effect of some inhibitors compounds [22, 23, 24, 25]. The use of Ulva sp. as co-substrate with the dairy raw effluents may enhance therefore the properties of the feedstock mixture and methane production, thus promoting the activity of methanogen bacteria compared to the sulfate-reducing ones, given that these properties were underlined during the co-digestion of this algae with other substrates [23, 26, 27]. Indeed, the algal biomass has the advantage to be biodegradable and has a high content of organic matter and other macronutrients [28].

The bicarbonate ion is often used in some applications as a buffering effect and can be added in the medium to buffer the pH variation [8, 9]. Furthermore, bacterial activity could also be improved by adding compounds that stimulate their growth [29]. The production of biogas using industrial effluents including additives media to activate methanogens bacteria remains a better way to treat these wastes, while reducing greenhouse gas emissions.

Previous studies have treated the effect of nitrogen, contained in the urea, and the oligo-elements on biogas production from sugarcane vinasse [16], and showed that a stable biogas production is obtained. This result is explained by the sufficient alkalinity of the medium even at higher OLR (9.6 gCOD⋅L−1⋅d−1) and lower HRT (2.5 days), since the addition of essential micronutrients stimulated the methanogenic activity avoiding major VFA accumulation. While the addition of phosphate in the form of KH2PO4 induced volatile fatty acids accumulation [16].

During the anaerobic digestion of raw by-products, it is necessary to add a buffering medium, since the rapid degradation of the soluble compounds causes the pH drop and could compromise the digester stability [30]. Co-digestion of a mixture containing dairy effluents is often responsible of the production of VFA. It is therefore important to test the effect of some additives media to remedy this problem. The aim of this study was to highlight the energetic valorization, using anaerobic digestion, of local wastes generated from dairy industry. The co-digestion of this effluent with the macroalgae Ulva sp. was also performed in order to explore this substrate present in significant quantities in the Algerian coast and responsible of the release of a large amount of H2S by natural fermentation. Additives media including bicarbonate or glucose for methanogen bacteria activation or pH stabilization have been used on the basis of our previous study [31] and the yields of anaerobic digestion have been compared to the co-digestion of DRM by-products and Ulva biomass. Biogas production and volatile fatty acids (VFAs) accumulation have been assessed during this study.

2. Materials and methods

2.1. Substrate and inoculum characterization

The DRM selected in this study are the effluents provided from the washing unit and the expired dairy products, characterized by their high COD. The plant is located at Boudouaou, near of Algiers (Algeria). This plant is ranked amongst the top ten in term of producing and selling of fruit juices in Algeria. The substrate to inoculum ratios (S/I) are taken 14.79 and 22.13 respectively for batch with nutrient medium I and medium II, respectively. The used Ulva sp. was collected from Bordj-El Bahri coastal marine areas, north of Algiers (Algeria). The macroalgae was dried at 70 °C during 24 h and then crushed according to Vergara-Fernández et al. [32]. The dairy raw material samples were collected fresh and stored immediately at 4 °C until subsequent use throughout the experimentation period.

2.2. Analytical methods

The organic matter degradation has been investigated by monitoring COD evolution according to the procedure reported by Wolf and Nordmann [33]. For carbohydrates characterization before and after the AD process, the phenol-sulfuric acid method has been used [34]. Proteins and lipids concentrations are determined respectively by using Lowry et al. [35] and Bligh and Dyer [36] methods. A high performance liquid chromatography (HPLC) (Agilent YL9100) has been used in order to estimate the evolution of the VFA concentrations using a C8 column (5 μm–25 × 0.46). The mobile phase is the sulphuric acid (0.005 M) at the flow rate was fixed to 0.7 mL/min. The column temperature was set at 30 °C and the detection was made spectrophotometrically at a wavelength of 210 nm. The pH, the redox potential, the conductivity and the temperature have been measured using an Inolab multi-parameter 720 device (Prolabmas, MurniSwadaya, Indonesia). Concerning the biogas composition, a coupling between gas chromatography and mass spectrometry (GC/MS) (Perkin Elmer Clarus 600, Perkin Elmer, Waltham, MA, USA) has been used. Percentages of methane and carbon dioxide (% v/v) present in the biogas have been determined by comparing the sample biogas with pure standard gases. The GC has been fitted with an Rtx®-VMS Column (60 m × 250 μm ID) (Restek Corporation, Bellefonte, PA, USA). Helium has been used as the carrier gas at a flow rate of 1 mL/min. The oven temperature has been set to 40 °C with a rate increase of 3 °C/min until the temperature of 180 °C has been reached, then it has been kept at this temperature for 3 min. The injector and detector temperatures have been set to 180 °C.

2.3. Experimental set-up

Experiments have been carried out in a cylindrical stainless steel digester at with a conical bottom (Figure 1). This digester has 15 cm of diameter and 30 cm of height with a total volume of 5 L. This geometry is well suited to perform anaerobic digestion of low-viscous effluent [37]. The mixing was promoted by an impeller of 0.14 m diameter, which was located at a distance about 1/3 above the bottom. This impeller turns at 150 rpm, inducing a radial flow inside the fermenter. The ratio of bioreactor diameter to impeller one was about 1.071. To heat and control the operating temperature during the AD process, the bioreactor was equipped wi th an integrated heating element and a thermostat. On the top of the reactor, several ports have been placed in the tight fitting lid, one serves to monitor the output of biogas, while the other control several parameters through inserted probes such as: pH, temperature, and to collect liquid samples for consequent analysis. The digester is fitted with a purge valve at the bottom of the conical arrangement for an easy draining of the effluent. A stainless steel flange has also been included in the digester in order to prevent gas leakage. The experiments have been carried out for an initial pH of 7.1 which is adjusted by adding NaOH or HCl (respectively at 2 N and 1.6 N).

Figure 1.

Experimental set-up. (1) bioreactor; (2) jacket; (3) pH control; (4) Rushton turbine; (5) electric heater; (6) drain valve; (7) recycling valve; (8) substrate; (9) gas outlet.

Four batch experiments (batch 1, batch 2, batch 3 and batch 4) were performed. The volume of dairy raw material of each digester was about 4 L. 0.5 L of inoculum contained medium I or medium II were added to the dairy raw material in the case of batches 2, 3 and 4 to activate methanogenic bacteria or to ensure stable operation, the volume choice is based on the C/N ratio and according to Zerrouki et al. [37]. The AD of the batch 1 has been performed without inoculum in order to highlight the contribution of endogenous lactic acid bacteria of the substrate in biogas production. 3 g of Ulva sp. having a COD of about 52 g/L was used as a co-substrate in batch 4. This algae mass choice is based on the corresponding I/S ratio, according to the study of Karray et al. [17], which showed that an optimal yield of biogas production was obtained at this condition.

Samples were taken daily from the sampling port to measure COD, VFA, lipids, protein, carbohydrates and pH and the tests were performed in duplicates. The digester was operated during 35 days.

2.4. Kinetic modeling

2.4.1. Biogas yield

Biogas yield was expressed as a function of the COD equivalent using the following equation:

yBiogas=VBiogasCOD(1)
where yBiogas is the biogas yield (Lbiogas⋅g−1⋅L), VBiogas the volume of produced biogas (L) and COD the degraded chemical oxygen demand (g⋅L−1).

2.4.2. Biogas production kinetics

Kinetic study of anaerobic digestion of DRM with additive mediums and Ulva sp. macroalgae as a co-substrate was carried out using a nonlinear kinetic model. Kinetic parameters are estimated by fitting the experimental cumulative biogas data with the modified Gompertz model [38]:

B(t)=PexpexpRmeP(𝜆t)+1(2)
where B(t) is the cumulative methane produced during the digestion time t (L∕gVS), P is the methane production potential (L∕gVS), Rm is the maximum methane production rate (L∕gVS⋅d), 𝜆 is the lag phase duration (d) and “e” is a mathematical constant (equal to 2.718281828).

Equation (2) has been regarded as a valuable empirical non-linear model to predict data of anaerobic digestion using different types of reactors [39, 40].

2.4.3. Kinetics of substrates biodegradation

The biogas production depends mainly of the amount of COD, presents under carbohydrates, proteins and lipids compounds. Food waste consists of this three main organic compounds, which lead into different biogas yields and bioconversion rates [41]. The removal rate of these substrates can be modelled using the following equation [42]:

C(t)=a+bekt(3)
where, C(t) is the amount of organic matter at time t, “a” is a constant, “b” is the amount of biodegradable organic matter, and “k” is the reaction rate constant.

2.4.4. Statistical analysis

Regression coefficients (R2) between experimental data and estimated during fermentation process, were estimated. Standard deviations (SD) and coefficients of variation (CV) were calculated for the kinetic parameters for each batch. The Tuckey-Kramer test (p < 0.05) was used for analysis of variance (ANOVA).

3. Results and discussion

3.1. Composition of dried Ulva sp.

Elementary analysis, provided the elementary composition of Ulva sp., is presented in Table 1. The values found revealed a high percentage of carbon and oxygen, this macroalgae contains N, S and microelements such as Cu, Ca, P, K and Fe which are considered as beneficial for the AD. In many previous works, these elements were used as additive elements to improve the AD yield [43]. Thus, the used of this algal biomass as a co-substrate is an important advantage due to its high content of nutrients and its thinner and simpler morphological structure, and a large surface area, with more efficient digestion compared to other macroalgae [44]. Moreover, this biomass offers more ecological characteristics, such as, biodegradability.

Table 1.

Elemental composition of Ulva sp.

Element wt%Elementwt%
C48.05
N5.93 Cl0.09
O31.92 K0.09
Na0.12 Ca4.84
Si0.28 Fe0.19
P0.18 Cu0.30
S4.73

Table 2.

Composition of additives media

Medium IMedium II
Component Concentration (g⋅L−1)ComponentConcentration (g⋅L−1)
NaHCO31.68C6H12O64.00
NH4Cl1.44Na2HPO40.90
MgSO4, 7H2O0.28MgSO4, 7H2O0.20
KH2PO40.60FeCl30.04
CaCl2, 2H2O0.05CaCl2, 2H2O0.10
Yeast extract0.10

Table 3.

Substrate and inoculum characterization

ExperimentBatch 1 Batch 2 Batch 3 Batch 4
Dairy raw materials (L) 4.5 4 4 4
Inoculum (L) 0.5 0.5 0.5
Ulva sp. macroalgae (g) 3
Additive mediumNutrient medium INutrient medium IINutrient medium I
C/N ratio 40.81 27.87 41.82 22.57
VS (g⋅L−1)96.6 ± 6.48999.4 ± 6.67799.5 ± 6.68455.89 ± 3.754
COD (gO 2⋅L−1)100 ± 6.717132 ± 8.867136 ± 9.13596 ± 6.448
pH 6.8 7 7 7
Conductivity (μS⋅cm−1)4.287.444.934.59
Carbohydrates (g⋅L−1)87.23 ± 5.85991.92 ± 6.17584.41 ± 5.67072.23 ± 4.852
MSS (mg⋅L−1)31 ± 2.0817 ± 0.4704 ± 0.26824 ± 1.612
Lipids (g⋅L−1)0.085 ± 0.0060.098 ± 0.0060.056 ± 0.0040.034 ± 0.002
Protein (g⋅L−1)5.347 ± 0.3595.937 ± 0.3985.285 ± 0.3553.855 ± 0.258
Ammonium Nitrogen (g⋅L−1)0.435 ± 0.0290.493 ± 0.0330.839 ± 0.0560.903 ± 0.061
Acetic acid (g⋅L−1)0.090 ± 0.0060.018 ± 0.0010.210 ± 0.0140.500 ± 0.034
Propionic acid (g⋅L−1)0.006 ± 0.00040.004 ± 0.0030.116 ± 0.0080.027 ± 0.002
Butyric acid (g⋅L−1)0.050 ± 0.0030.018 ± 0.0010.009 ± 0.00060.035 ± 0.0023

COD: Chemical Oxygen Demand; VS: Volatile solid; MSS: Dry matter in suspension.

3.2. Batch experiments

The composition and characteristics of substrates and inocula are given in Tables 2 and 3. Additives media compositions were chosen based on the synthetic salt often used in the sludge activity tests. Glucose was chosen in medium II because it was often used as a co-substrate to support microbial activity [45] and enhances methane production [37]. The iron was often supplied in the form of FeCl3 in nutrient medium to reduce the H2S content in the biogas composition and to further reduce the cost for H2S treatment. AD process required macronutrients that can be supplied in the form of Na2HPO4, MgSO4, CaCl2, they were considered as essential elements of biomass and can be used as buffering agent to keep digester stability [46, 47]. The effluents were selected according to their organic loads and C/N ratios. The initial C/N was about 40.81, 27.87, 41.82 and 22.57, respectively for batches 1, 2, 3 and 4. According to Chandra et al. [48], the optimum C/N ratio is between 20 and 40. The C/N ratio in algal biomass is low, it is about 10 which leads to ammonia accumulation. To overcome this accumulation, addition of organic materials at high COD rate is often required in order to enhance the digestion process. The initial chemical organic load was about 100, 132, 136 and 96 gO 2⋅L−1 respectively for batch 1, batch 2, batch 3 and batch 4. The analysed samples are characterized by a high COD load which corresponds to the values reported by Escalante et al. [49]. The chemical composition of effluents is characterized by high carbohydrates and proteins and low lipids contents making them more suitable for methane production.

Figure 2.

Methane yield production in batch reactors (T = 35 °C, pH = 7).

Table 4a.

Kinetic model parameters of modified Gompertz model

Batch 1Batch 2Batch 3Batch 4
Methane yield (LCH 4∕gVS)0.07920.2080.1120.118
P (LCH 4∕gVS)0.0910.2270.1120.118
Rmax (LCH 4∕gVS⋅d)0.00340.00910.00680.0086
𝜆 (d)3.0574.6250.0750.063
R20.9930.9950.9740.973
RMSE 0.0027 0.0052 0.0077 0.0082

3.3. Methane production analysis under the different substrates mixing

The experimental profiles of the cumulative biogas volume for batch reactors performed under mesophilic digestion (35 °C) are depicted in Figure 2. It can be seen that the cumulative methane tends towards a maximum value of 0.208 LCH4∕gVS in the case of batch 2, which shows the high potential of this complex substrate for methane production with longer production period. While, the cumulative biogas production exhibits an asymptotical profile for batches 1, 3 and 4. Table 4a summarizes the kinetic parameters of the modified Gompertz model. The maximum methane production rate (Rm) was more important in batch 2 than the other batches, it was about 0.0091 LCH4∕gVS⋅d, while it was only about 0.0034, 0.0068 and 0.0086 LCH4∕(gVS⋅d) in batches 1, 3 and 4, respectively. This could be explained by the high methane yield after 15 days of AD with a high slope compared to the other batches (Figure 2). This result allows to expect a high methane productivity during the co-digestion of DRM and Ulva sp. biomass. Modelled lag phase duration, 𝜆, was significantly reduced in batches 3 and 4, the values obtained were 0.075 d and 0.063 d, respectively, which indicates that the methane production was beginning immediately on the first day of digestion. Thus, the addition of nutrient medium I with Ulva sp. biomass (batch 4) or the use of nutrient medium II (batch 3) shows a low lag phase, while this parameter is higher in batches 1 and 2 (3.057 d and 4.625 d, respectively). In this case, the undissolved compounds such as proteins and fats take more time to decay into monomers [50]. In this study, the hydrolysis constant remains low compared to other studies which used the dairy by-products. It was in the range of 0.0009–0.104 d−1 for all batches (Table 4b). Indeed, the value of 0.084 d−1 was previously reported by Kafle and Chen [39] for dairy manure, while the range of 0.36–0.55 d−1 has been previously reported by Mainardis et al. [51] for some cheese whey varieties. It should be noted that a good agreement is observed between the methane yield reached (LCH 4∕gVS) and the predicted data using (2), with correlation coefficients R2 varied from 0.973 to 0.995 and the RMSE value fell within the range of 0.0027–0.0082.

Table 4b.

Kinetic model parameters of first order model

Batch 1Batch 2Batch 3Batch 4
k (d−1)0.00750.00090.08940.1040
R20.9930.9800.9790.975
RMSE 0.0024 0.0106 0.0053 0.0064

The specific methane yields gradually increasing and then levelling off at 0.208 LCH 4∕gVS, 0.112 LCH 4∕gVS, 0.118 LCH 4∕gVS and 0.079 LCH 4∕gVS, for batches 2, 3, 4 and 1, respectively after 35 days of AD. Macroalgae added in batch 4 allowed to bring essential elements and to create a synergistic effect by alleviating the imbalance of nutrients. However, lipids accumulation in batch 4 causes a slower hydrolysis rate, resulting in a lower methane yield compared to batch 2 and confirms results obtained by Rodolfi et al. [52]. The methane production is more meaningful in batch 2 than the other batches, certainly due to the presence of inorganic additive (medium I) which ensures an appropriate nutrient balance between C and N, in addition to the digester stability. This substrate has a good biodegradability of the complex organic matter without any inhibitory effect upon the biomass. These findings were consistent with BMP values of the literature for dairy wastes and Ulva sp. It was about 0.208LCH4gVS1 in the case of dairy residues, while in the case of Ulva sp., the BMP values varied between 0.094–0.177LCH4gVS1 depending upon the pretreatment techniques employed [26].

Compared to the literature, the methane production from batch 2 can be considered as a promising way forward regarding the use of combined treatment of dairy residues in AD. Lhanafi et al. [53] previously reported a maximum methane yield of about 0.093 LCH 4∕gVS. These authors added a mixture of dairy residues into inoculum which included diluted sludge of methanogenic bacteria in addition to formic, propionic and lactic acids. Some authors have underlined a methane yield of 0.197 LCH 4∕gVS from anaerobic digestion of manufacturing milk. Their yield has been improved on the basis of mixing ratios performed between food waste leachate from dairy factory and wastewater treatment sludge, the methane produced varied from 0.118 to 0.223 LCH 4∕gVS [54]. Low methane yield has been recorded from cheese whey wastewater effluents; it was only about 0.065 LCH 4∕gVS [55]. Significant methane yields have been reported by Mainardis et al. [51] which ranged from 0.353 to 0.437 NL CH4∕gVS from a variety a cheese whey. The difference may be due to the operating conditions and the composition of substrate. Indeed, these authors have mixed the sludge with dairy substrate at different ratios and the obtained mixtures are characterized by low organic content (50.10–78.73 g/L) compared to our study (96–136 g/L). Besides that, no-pretreatment of substrate has been performed in this work. Moreover, the co-digestion of the mixture of cow dung and Ulva sp. led to methane yield of 0.195 LCH 4∕gVS at Ulva sp.: cow dung ratio of 1:1 [56] which is higher to the methane yield found in this study at Ulva sp. DRM ratio of 1/10.

In addition, it was also known that the AD of this macroalgae specie gave low potential methane (0.097 L∕gVS according to Vanegas and Bartlett [57]). This low biogas yield was due to the low VS content and a late methanogenic phase [58]. However, in this study, it was shown that it is possible to valorize this type of waste, which invades the Algerian coasts, by anaerobic digestion. Even if the biogas yield was lower in the batch containing this macroalgae, it is observed that the quality of this biogas is better with higher methane content.

Figure 3.

Comparison of modelled (according to (3)) and experimental chemical oxygen demand for batch reactors (T = 35 °C, pH = 7).

The co-digestion of DRM with Ulva sp. macroalgae overcomes this problem of valorization of the algae effluent present in large quantity on the Mediterranean coast, despite the low initial COD value of the mixture (96 gO 2⋅L−1) compared to the DRM where no inoculum was added (100 gO 2⋅L−1 at batch 1). Note that in the anaerobic digestibility process of DRM, the lactic and non-lactic acid bacteria play an important role in organic compounds degradation to methanogenic substrates [59]. Then, the most biodegradable organic fraction can easily be converted into biogas. This case is observed in batch 1 where no inoculum was added. The biogas is produced only from the endogenous flora, mainly lactic acid bacteria. The dimensionless number of COD (C∕C0) of batches 1, 2, 3 and 4, were 0.28, 0.15, 0.475 and 0.291, respectively after 30, 30, 30 and 26 days of AD (Figure 3) with more regular kinetic of degradation for batches 2 and 4 (containing nutrient medium I) than the batch 3 (with nutrient medium II). The addition of organic compounds in form of glucose constitutes more carbon source for anaerobic digestion while the addition of bicarbonate as a pH buffer medium, in one side, and the C/N ratio, close to the mid-value of the recommended interval (27.87), in the other side, enhance both methane production and COD removal.

In this last case, rapid degradation kinetics were observed for the first 10 days followed by a stagnation of the C∕C0 ratio, probably due to the acidification of the medium due to the production of volatile fatty acid. The comparison between the experimental and the predicted data using (3) showed the good fitting of data with correlation coefficients R2 higher than 0.902 and the RMSE value within the range of 0.026–0.079 (Table 5). The models tested were well suited to describe the methane production behavior and organic matter content degradation in batch fermentations.

Figure 4.

Effect of nutrient media on biogas composition.

Table 5.

Kinetic parameters of COD

Batch 1Batch 2Batch 3Batch 4
a0.2970.1420.3910.230
b0.6400.8620.6090.775
kh0.0040.0030.0310.004
R20.9210.9890.9020.988
RMSE 0.066 0.026 0.079 0.032

3.4. Biogas composition

The biogas composition, given in Figure 4 shows high methane contents for all batches. It was about 92, 65, 81 and 96% for batches 1, 2, 3 and 4, respectively. The batch 4 which contains Ulva sp. macroalgae provides better methane content (about 96%) compared to the other ones.

It is noted that, under investigated mediums, the hydrogen sulphide in the biogas produced was about 0.67, 1.92, 1.82 and 1.00% in the case of batch 1, batch 2, batch 3 and Batch 4, respectively. There is a direct relationship between the inoculums, the methane production and the presence of organic and inorganic compounds. The variation in the methane content can be explained by the presence of long chain fatty acid which was resulted from the dairy fat conversion. as a consequence, the decrease in physiological activity (ATP concentration), in parallel to the decrease of methanogenic bacteria activity [60]. The addition of mineral additive, under the form of bicarbonate, ensures better methanogenic bacteria activation and the stabilization of the AD process. Moreover, the anaerobic co-digestion of dairy raw materials (DRM) with macroalgae (Ulva sp.) using the bicarbonate as a buffering medium has been selected to avoid ammonium accumulation and to enhance methane content.

Figure 5.

Evolution of organic material degradation during anaerobic digestion of dairy raw materials: (a) lipids, (b) protein and (c) carbohydrate.

3.5. Anaerobic digestion of lipid compounds

The DRM used in this study exhibits relatively high concentration of lipids (about 0.085 g⋅L−1). In Figure 5a, which depicts the dimensionless lipids concentration over time, it can be seen that the lipids content decreases from 1 until the value of 0.325 in batch 1. Addition of synthetic mediums to dairy raw materials improved the lipids uptake, which drops from 1 to 0.351 and 0.158 for batches 3 and 2, respectively. In such case, the conversion of lipids into fatty acids became easier and fast till the near total uptake. The lipid degradation of DRM was lower in the case of batch 2 compared to the batch 3, with an average degradation rate of about 65% and 84.2%, respectively. This result can be related to the effect of inorganic additive (based on bicarbonate) or organic additive (based on glucose) into the medium on the anaerobic digestibility of DRM, which in turn improved both lipids uptake and biogas production. It is noted that there was competition between bacterial degradation of lipids and proteins. It is also noted that the milk fat has an inhibiting effect on the activity of methanogenic bacteria [45]. For Batch 4, the analysis of lipids degradation distinguished two phases (Figure 5). The first phase is between 1 and 0.573 and the second one from 0.573 to 0.855 where the lipids increase after 26 d of anaerobic digestion.

This is due to the accumulation of intermediate compounds such as the long chain fatty acids (LCFAs) and the volatile fatty acids (VFAs) [61], that affect the lipid degradation process, this then represents a limiting mass transfer between substrate and microbial cells. This mass transfer depends to difficulty in transport and transform from liquid phase to solid and biological degradation of LCFAs [62]. Adsorption onto cell wall is restricted by bacteria coated LCFAs thus decrease the substrate degradation and yield of biogas [63].

3.6. Anaerobic digestion of protein compounds

The protein content is one of the most important criteria for the milk producers especially for casein and nitrogen sources in general. Figure 5b gives the protein concentration evolution for a series of batch reactors. It can be seen that at the beginning, the protein levels dropped quickly for the proposed series. This means that the lactic and non-lactic acid bacteria converted both the proteins and the lipids compounds, but the conversion rate was more important for proteins. The reaction rate constants reached values of 0.008, 0.012, 0.009 and 0.026 h−1, in the case of batches 1, 2, 3 and 4, respectively.

It can be seen also that there is a good fit between the experimental and modelled curves for the proteins uptake in the case of batches 2 and 3, where the R2 and RMSE values were respectively 0.985 and 0.032 for the first case and 0.943 and 0.079 for the second one (Table 6). However, for batches 1 and 4, the R2 values showed that the fitted data by the first order kinetic model for proteins uptake are deviated from the experimental ones (0.818 and 0.823, respectively).

Table 6.

Kinetic parameters of protein

Batch 1Batch 2Batch 3Batch 4
a0.2650.1610.1410.292
b0.7510.8450.8130.678
k0.0080.0120.0090.026
R20.8180.9850.9430.823
RMSE 0.096 0.032 0.073 0.111

Similar findings have been previously reported by [60]. In such case, the dimensionless protein concentrations decreased from 1 until the average values of 0.287, 0.183 and 0.122 after 30, 23 and 25 days of digestion in batches 1, 2 and 3 respectively.

This decrease means that the nitrogen was used as a nutrient by methanogens and lactic acid bacteria, and also converted into ammonia as previously pointed out by [64]. It has been found that a good rate of proteins assimilation by methanogens bacteria has been occurred in batches 2 and 3, particularly, when the sodium bicarbonate has been added to the medium, probably due to the effect of this compound on the pH of the fermentation medium. Indeed, according to Browne and Murphy [65], the methane production reduction can be induced by the increase of the total ammonia nitrogen rate in the digester at high VS organic loads. Thus, sodium bicarbonate addition could moderate this effect. Note that the medium II contains glucose that constitutes other carbon source, easily biodegradable, and can explain the utilization of proteins and lipids rather than carbohydrates in the AD. In the case of batch 4, the dimensionless protein concentration goes down to 0.264 after 22 days of fermentation with uptake of 74% and a maximum reaction rate constant of 0.026 h−1. This behavior is certainly due to the complex composition of DRM with Ulva sp. which contains some types of proteins, carbohydrates and lipids with different biodegradation behavior. Indeed, it is known that carbohydrates are more biodegradable in anaerobic medium than proteins and lipids, which take a longer time for hydrolysis [66]. Moreover, it has been found an increase of the ammonia concentration to the values of 127.46 mg⋅L−1 and 119.86 mg⋅L−1 in batch 4 and batch 2, respectively which remain below the value 200 mg⋅L−1 previously reported by Elbeshbishy and Nakhla [67].

3.7. Anaerobic digestion of carbohydrate compounds

Comparison of modelled and experimental carbohydrates digestion results for a series of batch reactors is given in Figure 5c. It shows that the carbohydrates uptake occurred gradually in the same way for a series of batch reactors but at different rates. The uptake of the proteins and the carbohydrates was performed exponentially and it seemed more pronounced and extremely fast in the case of proteins. The analysis of the reaction rate constant (k) confirms this remark (Table 7). Indeed, for all batches, the constants kinetic, k, of carbohydrates degradation showed lower values compared to those of proteins (0.008, 0.012, 0.009, 0.026 h−1 against to 0.002, 0.003, 0.006 and 0.011 h−1 in batches 1, 2, 3 and 4, respectively for proteins and carbohydrates). This result can be ascribed to the carbohydrates content in the batch reactors.

Figure 6.

Evolution of acids production during anaerobic digestion of dairy raw materials. (a) Acetic concentration, (b) butyric concentration.

Table 7.

Kinetic parameters of carbohydrates

Batch 1Batch 2Batch 3Batch4
a0.03400.2310.075
b0.9941.0000.6940.877
k0.0020.0030.0060.011
R20.9650.9620.9170.985
RMSE 0.057 0.063 0.072 0.036

In all cases, carbohydrates concentration decreases were recorded at the beginning of the digestion, i.e. during the first 6 days. Note that in batch 1, the carbohydrates uptake reached the efficiency of 76%. This value remained the lowest value reached in this study under the investigated conditions (without inoculum). While it was improved in batches 2 and 3 and 4 leading to uptake efficiencies of about 85, 74 and 91%, respectively. The results found indicates that the medium I was well suited to improve biogas production with respect to the carbohydrates, lipids and proteins biodegradation characteristics. Thus, the best biogas production has been reached when using the medium I which has improved the degree of organic load biodegradability, the microbial growth and their capacity of protein, carbohydrates and lipids uptake. Figure 5c also shows a close fit between the modelled and experimental profiles of carbohydrates degradation. Values of R2 and RMSE were 0.965, 0.962, 0.917, 0.985 and 0.057, 0.063, 0.072, 0.036 respectively for batches 1, 2, 3 and 4 (Table 7).

3.8. Volatile fatty acids production

Figures 6 (a, b) give the volatile fatty acids (VFAs) evolution during the digestion of dairy raw materials for a series of batch reactors. It can be seen that the main acids produced during anaerobic digestion of DRM were acetic and butyric acids. Further traces of propionic acids have been detected, although these are no direct products issued from the anaerobic digestibility of DRM. These figures show also that the concentrations of acetic acid are greater than those of butyric one over the 35 days of AD, probably due to the activity and/or the concentration of lactic acid bacteria which could not stimulate total methane production, but promote the rate of methane production at the beginning of the AD [68]. In the case of batch 4, an unexpected peak of acetic acids (1.228 g⋅L−1) was observed after 15 days of AD which results in a stabilization of the production of biogas from this same day (Figure 2). This peak may be assigned to the accumulation of hydrolysis products. The use of inorganic additive (based on bicarbonate) on the anaerobic digestibility of DRM plays an important role regarding the volatile fatty acids production and the digester stability. However, compared to the batch 2, this peak certainly has a detrimental effect on the microbial activity which would explain the yield of AD obtained for the batch 4. The same observation could be made for batch 3 with acetic acid concentration peak observed on the first day of AD. After that, concentrations of this acid varied from 0.754 to 0.041 g⋅L−1, while those of the butyric acid increased until the value of 0.221 g⋅L−1. As a consequence, the stabilization of the biogas production form 15 days of the AD. This may be due to a direct relationship between the inoculum, the methane production rate and the presence of fatty volatile acids [69].

For batch 2, the acetic acid concentration remains constant at about 2 g/L from 15 days which would explain the biogas yield obtained with this batch, even if the butyric acid concentration continues to increase (about 0.287 g⋅L−1). This can be due to the conversion of lactate from glucose by Butyri bacterium methylotrophicum as previously reported by (Shen et al. [70]). The increase in butyric acid can stimulate methanogens bacteria rather than their inhibition Stronach et al. [10]. Volatile fatty acids are the intermediate coproducts in the biogas production during the anaerobic digestion process, the levels reached are not toxic as long as batch 2 ensures buffering capacity, using additive medium I, to keep the process stability. While, the co-digestion of DRM and Ulva sp. (batch 4) resulted in the accumulation of hydrolysis products, which explains the increase in the concentration of the acetic acid.

Furthermore, for a series of batch reactors a continuous mechanization mode has been suggested to avoid accumulation of volatile acids and keep substrate concentrations constant. The continuous mode can be expected to provide better gas efficiency of the pilot-scale dairy plants.

4. Conclusion

This study has been carried out to evaluate anaerobic co-digestion performance of dairy raw materials (DRM) with macroalgae (Ulva sp.) In addition, the effect of organic additive (based on glucose) and mineral additive (based on bicarbonate) on the anaerobic digestibility of DRM was also evaluated. The dairy raw material containing the nutrient medium I (batch 2) is more suitable to improve biogas production, the highest specific methane yield reached, it was about 0.208 LCH4∕gVS after 35 days. This nutrient medium allows to improve the degree of anaerobic digestibility of DRM, the microbial growth and their capacity of protein, carbohydrates and lipids uptake, it may be a useful pattern to convert food industry streams.

The co-digestion of DRM with Ulva sp. resulted in higher biogas yield compared to DRM only-fed condition and enhances the biogas quality which contains only 0.5 to 2% of hydrogen sulphide. The methane production is more meaningful in batch 2 than the other batches, certainly due to the presence of inorganic additive (medium I) which ensures an appropriate nutrient balance between C and N, in addition to the digester stability. Lipids accumulation in batch 4 causes a slower hydrolysis rate, resulting in a lower methane yield compared to batch 2. Effect of organic and inorganic additives in a series of anaerobic laboratory tests to produce biogas from dairy residues and macroalgal biomass can give a contribution to sustainability and circular economy perspective.

The used models have well described the cumulative methane production and also the kinetic of the removal organic matters (proteins, lipid and carbohydrates). This analysis showed that these three compounds follow a different kinetics according to the composition of the medium and the addition of organic or inorganic media. Furthermore, the continuous mode and the co-digestion of DRM with Ulva sp. by optimizing the I/S ratio can be expected to provide better methane yield.

Nomenclature

VBiogas

Biogas volume produced during batch experiments (m3).

COD

Chemical oxygen demand, gO 2⋅L−1.

BP

Cumulative biogas yield during the fermentation

P

Cumulate biogas production (N⋅mL)

kh

First order kinetic constant (h−1)

t

Time (h).

Lu

Total amount of biodegradable organics matter (gCOD∕L)

MS

Mineral solids

R2

Correlation coefficient

T

Temperature (°C)

VFA

Volatile fatty acids

y(t)

Amount of organic matter at time t

a

Constant

b

An initial amount of biodegradable organic matter, and k is the reaction rate constant.


References

[1] IEA Key World Energy Statistics, International Energy Agecy, Paris, France, 2016

[2] M. Jeguirim; S. Jellali; B. Khiari C. R. Chim., 23 (2020), pp. 583-587

[3] N. Boukaous; L. Abdelouahed; M. Chikhi; Ch. Mohabeer; B. Taouk; A.H. Meniai C. R. Chim., 23 (2020), pp. 623-634

[4] C. Tanios; Y. Saadeh; M. Labaki; M. Boutros; C. Gennequin; H. L. Tidahy; A. Aboukaïs; E. Abi-Aad C. R. Chim., 23 (2020), pp. 1-19

[5] A. Matos; S. Abrahao; P. Lo Monaco; A. Sarmento; M. Matos Extraction capacity of plant grown in constructed wetland systems used to dairy wastewater treatment, 12th International Conference on Wetland Systems for Water Pollution Control (2010)

[6] MarketLine Dairy Global Industry Guide 2015-2024, 2020 (accessed 18.06.12), http://www.researchandmarkets.com/research/e7e4d9/dairy_global_indu

[7] E. Buabeng-Baidoo; N. Mafukidze; J. Pal; S. Tiwari; B. Srinivasan; T. Majozi; R. Srinivasan Chem. Eng. Res. Des., 121 (2017), pp. 81-91 | DOI

[8] A. Akbi; M. Saber; M. Aziza; N. Yassaa Renew. Sustain. Energy Rev., 72 (2017), pp. 240-245 | DOI

[9] M. Fahim; M. Hassanein; A. Abou Hadid; M. Kadah Impacts of climate change on the widespread and epidemics of some tomato diseases during the last decade in Egypt, III International Symposium on Tomato Diseases, Volume 914 (2010), pp. 317-320

[10] S. Stronach; T. Rudd; J. Lester Toxic Substances in Anaerobic Digestion, Anaerobic Digestion Processes in Industrial Wastewater Treatment (Biotechnology Monographs), Volume 2, Springer-Verlag, 1986, pp. 71-92 | DOI

[11] B. Demirel; O. Yenigun; T. T. Onay Process. Biochem., 40 (2005), pp. 2583-2595 | DOI

[12] J. Rivas; A. R. Prazeres; F. Carvalho J. Agric. Food Chem., 59 (2011), pp. 2511-2517 | DOI

[13] F. Carvalho; A. R. Prazeres; J. Rivas Sci. Total Environ., 445 (2013), pp. 385-396 | DOI

[14] G. Luo; I. Angelidaki Appl. Microbiol. Biotechnol., 97 (2013), pp. 1373-1381 | DOI

[15] S. Gao; Y. Huang; L. Yang; H. Wang; M. Zhao; Z. Xu; Z. Huang; W. Ruan Energy Convers. Manag., 93 (2015), pp. 166-174 | DOI

[16] L. Janke; S. Weinrich; A. F. Leite; H. Sträuber; C. M. Radetski; M. Nikolausz; M. Nelles; W. Stinner Energy Convers. Manag., 168 (2018), pp. 188-199 | DOI

[17] R. Karray; F. Karray; S. Loukil; N. Mhiri; S. Sayadi Waste Manag., 61 (2017), pp. 171-178 | DOI

[18] M. El Ati-Hellal; F. Hellal; M. Dachraoui; A. Hedhili C. R. Chim., 10 (2007), pp. 839-849 | DOI

[19] J. Kim; H. Kim; G. Baek; C. Lee Waste Manag., 60 (2017), pp. 322-328 | DOI

[20] G. Bonanno; M. Orlando-Bonaca Ecotoxicol. Environ. Safety, 148 (2018), pp. 44-71 | DOI

[21] P. Peu; J.-F. Sassi; R. Girault; S. Picard; P. Saint-Cast; F. Béline; P. Dabert Bioresour. Technol., 102 (2011), pp. 10794-10802 | DOI

[22] X. Dai; N. Duan; B. Dong; L. Dai Waste Manag., 33 (2013), pp. 308-316 | DOI

[23] X. Fonoll; S. Astals; J. Dosta; J. Mata-Alvarez Chem. Eng. J., 262 (2015), pp. 1268-1274 | DOI

[24] A. Khalid; M. Arshad; M. Anjum; T. Mahmood; L. Dawson Waste Manag., 31 (2011), pp. 1737-1744 | DOI

[25] J. G. Usack; L. T. Angenent Water Res., 87 (2015), pp. 446-457 | DOI

[26] X. Briand; P. Morand J. Appl. Phycol., 9 (1997), pp. 511-524

[27] A. Bruhn; J. Dahl; H. B. Nielsen; L. Nikolaisen; M. B. Rasmussen; S. Markager; B. Olesen; C. Arias; P. D. Jensen Bioresour. Technol., 102 (2011), pp. 2595-2604 | DOI

[28] M. E. Montingelli; S. Tedesco; A. G. Olabi Renew. Sustain. Energy Rev., 43 (2015), pp. 961-972 | DOI

[29] S. Arif; R. Liaquat; M. Adil Renew. Sustain. Energy Rev., 97 (2018), pp. 354-366 | DOI

[30] A. P. Flores-Mendoza; H. Hernández-García; Y. Cocotle-Ronzón; E. Hernandez-Martinez J. Chem. Technol. Biotechnol., 95 (2020), pp. 1946-1952 | DOI

[31] J. Fernández-Rodríguez; E. Robles; O. Gordobil; M. González-Alriols; J. Labidi Chem. Eng. Trans., 52 (2016), pp. 187-192

[32] A. Vergara-Fernández; G. Vargas; N. Alarcón; A. Velasco Biomass Bioenergy, 32 (2008), pp. 338-344 | DOI

[33] V. Wolf; W. Nordmann Abwasser, 24 (1977), pp. 277-279

[34] M. Dubois; K. A. Gilles; J. K. Hamilton; P. t. Rebers; F. Smith Anal. Chem., 28 (1956), pp. 350-356 | DOI

[35] O.-I. Lowry; N. Rosebrough; A. Farr; R. Randall J. Biol. Chem., 193 (1951), pp. 265-275 | DOI

[36] E. G. Bligh; W. J. Dyer Can. J. Biochem. Physiol., 37 (1959), pp. 911-917 | DOI

[37] S. Zerrouki; R. Rihani; F. Bentahar; K. Belkacemi Water Sci. Technol., 72 (2015), pp. 123-134 | DOI

[38] J.-J. Lay; Y.-Y. Li; T. Noike Water Res., 31 (1997), pp. 1518-1524 | DOI

[39] G. K. Kafle; L. Chen Waste Manag., 48 (2016), pp. 492-502 | DOI

[40] G. K. Kafle; S. H. Kim; K. I. Sung Bioresour. Technol., 127 (2013), pp. 326-336 | DOI

[41] G. Esposito; L. Frunzo; A. Giordano; F. Liotta; A. Panico; F. Pirozzi Rev. Environ. Sci. Bio./Technol., 11 (2012), pp. 325-341 | DOI

[42] K. Bernat; A. Cydzik-Kwiatkowska; M. Zielińska; I. Wojnowska-Baryła; J. Wersocka Biochem. Eng. J., 148 (2019), pp. 87-96 | DOI

[43] M. Dębowski; M. Kisielewska; J. Kazimierowicz; A. Rudnicka; M. Dudek; Z. Romanowska-Duda; M. Zieliński Energies, 13 (2020), pp. 2186-2198 | DOI

[44] J. Costa; P. Gonçalves; A. Nobre; M. Alves Bioresour. Technol., 114 (2012), pp. 320-326 | DOI

[45] M. D. Khan; N. Khan; A.-S. Nizami; M. Rehan; S. Sabir; M. Z. Khan Bioresour. Technol., 238 (2017), pp. 492-501 | DOI

[46] A. Hendriks; J. Van Lier; M. De Kreuk Biotechnol. Adv., 36 (2018), pp. 1-13 | DOI

[47] M. Romero-Güiza; J. Vila; J. Mata-Alvarez; J. Chimenos; S. Astals Renew. Sustain. Energy Rev., 58 (2016), pp. 1486-1499 | DOI

[48] R. Chandra; H. Takeuchi; T. Hasegawa Renew. Sustain. Energy Rev., 16 (2012), pp. 1462-1476 | DOI

[49] H. Escalante; L. Castro; M. Amaya; L. Jaimes; J. Jaimes-Estévez Waste Manag., 71 (2018), pp. 711-718 | DOI

[50] I. Angelidaki; S. Petersen; B. Ahring Appl. Microbiol. Biotechnol., 33 (1990), pp. 469-472 | DOI

[51] M. Mainardis; S. Flaibani; M. Trigatti; D. Goi J. Environ. Manag., 246 (2019), pp. 557-563 | DOI

[52] L. Rodolfi; G. Chini Zittelli; N. Bassi; G. Padovani; N. Biondi; G. Bonini; M. R. Tredici Biotechnol. Bioeng., 102 (2009), pp. 100-112 | DOI

[53] S. Lhanafi; Z. Anfar; B. Chebli; M. Benafqir; R. El Haouti; Y. Azougarh; M. Abbaz; N. El Alem Sustainable Environ. Res., 28 (2018), pp. 389-395 | DOI

[54] Y. Yoon; S. Lee; K. Kim; T. Jeon; S. Shin J. Mat. Cycles Waste Manag., 20 (2018), pp. 283-292 | DOI

[55] M. Carlini; S. Castellucci; M. Moneti Energy Procedia, 82 (2015), pp. 811-818 | DOI

[56] V. Akila; A. Manikandan; D. S. Sukeetha; S. Balakrishnan; P. M. Ayyasamy; S. Rajakumar Biocatalysis Agric. Biotechnol., 18 (2019), 101035

[57] C. Vanegas; J. Bartlett Environ. Technol., 34 (2013), pp. 2277-2283 | DOI

[58] M. Tomei; C. Braguglia; G. Mininni Bioresour. Technol., 99 (2008), pp. 6119-6126 | DOI

[59] J. Bohn; A. Yüksel-Dadak; S. Dröge; H. König J. Biotechnol., 244 (2017), pp. 4-15 | DOI

[60] M. Perle; S. Kimchie; G. Shelef Water Res., 29 (1995), pp. 1549-1554 | DOI

[61] S. Park; Y. Li Bioresour. Technol., 111 (2012), pp. 42-48 | DOI

[62] A. Rinzema; M. Boone; K. van Knippenberg; G. Lettinga Water Environ. Res., 66 (1994), pp. 40-49 | DOI

[63] M. Pereira; D. Sousa; M. Mota; M. Alves Biotechnol. Bioeng., 88 (2004), pp. 502-511 | DOI

[64] R. Li; S. Chen; X. Li Appl. Biochem. Biotechnol., 160 (2010), pp. 643-654 | DOI

[65] J. D. Browne; J. D. Murphy Renew. Energy, 71 (2014), pp. 69-76 | DOI

[66] Y. Miron; G. Zeeman; J. B. Van Lier; G. Lettinga Water Res., 34 (2000), pp. 1705-1713 | DOI

[67] E. Elbeshbishy; G. Nakhla Bioresour. Technol., 116 (2012), pp. 170-178 | DOI

[68] S. Menardo; P. Balsari; E. Tabacco; G. Borreani Bioenergy Res., 8 (2015), pp. 1810-1823 | DOI

[69] T. Gemechu Afr. J. Food Sci., 9 (2015), pp. 170-175

[70] G.-J. Shen; B. Annous; R. Lovitt; M. Jain; J. Zeikus Appl. Microbiol. Biotechnol., 45 (1996), pp. 355-362 | DOI


Comments - Policy