3.2 COD analysis
Hydrolysis and acidogenesis increase the solubilization of FW, which in turn increases the concentration of sCOD. Results depict that the sCOD production in mesophilic conditions increased until 24 h (20.65, 28.5, 30.87, 31.88, 32.2 g/L) after which it started decreasing till 48 h (11.3, 11.86, 12.25, 15.55, 16.38 g/L) for 2, 4, 6, 8 and 10 % GC, respectively (Fig. 1a). After 48 h, the sCOD concentrations started increasing until 96 h, in the experiments with 2, 4 and 6 % GC (15.2, 14.65, 13.78 g/L). Nevertheless, the sCOD concentrations continued to decrease for the experiments with 8% GC (13.5 g/L) and 10% GC (14 g/L) until 96 h. In the control experiment, sCOD production increased to 15 g/L after 24 h and kept increasing to a maximum of 18 g/L at 48 h, after which it decreased gradually to 7.65 g/L at 96 h. The addition of GC increased the sCOD production by 1.37 to 2.14 folds within 24 h, and the cumulative sCOD concentrations remained higher until the end of the experiment, compared to control. Macias-Corral et al. (2017) and Darimani and Pant (2020) reported that the addition of GC increased the substrate's solubilization, consequently increasing sCOD concentration. Compared to the experiments with 8 and 10 % GC, after 24 h, the sCOD concentration decreased rapidly and steadily in the digester bottles with 2, 4 and 6 % GC addition.
From Table 2, it is evident that the addition of 2, 4, 6, 8, and 10 % GC demonstrated increased COD removal efficiency of 48.86, 49.95, 50.92, 47.39 and 47.22 %, respectively, compared to the control (36.93 %). The results show that the hydrolysis and COD removal rates were minimal in the digester bottles with 8 and 10 % addition of GC because the lignin contents shield carbohydrates (cellulose and hemicelluloses) from hydrolytic enzyme activities (Ferdes et al. 2020). A tight extracellular matrix in GC and other lignocellulosic substrates prevents hydrolytic enzymes' penetration and must be pretreated thermally and chemically along with mechanical treatment to improve the rate of enzyme hydrolysis (Chakraborty and Venkata Mohan 2019; Venkata Mohan et al. 2020).
3.3 Fatty acid and alcohol profiles
After 24 h, the total VFA production was 14.09, 15.46, 16.21, 16.69, 11.50 and 8.67 g/L for 10, 8, 6, 4, 2 % of GC addition and control, respectively (Fig. 1b). However, at 48 h, the VFA concentrations decreased to 4.73, 4.80, 4.63, 4.21, 4.65 g/L in digester bottles fed with 10, 8, 6, 4, 2 % of GC, respectively. In control, the total VFA increased further to 10.14 g/L at 48 h and then reduced. Between 48 h and 96 h, the total VFA production increased minimally, except for 10 % GC. The addition of 4 and 6 % GC supported acidogenesis with good VFA production, whereas the acidogenesis process was slower with 8 and 10 % GC. Due to enhanced acidogenesis, the degree of acidification was maximum in 2, 4 and 6 % addition of GC (55.70, 58.56 and 57.83 %) than 8 % GC addition (48.49 %), 10 % GC addition (43.77 %) and in control (52.51 %).
On a similar line, it was observed that the total VFA concentration was considerably higher when sewage sludge was anaerobically co-digested with FW, wastewater and GC when compared to the sludge being used as a single substrate (Fitamo et al. 2016). Likewise, in the current investigation, the cumulative VFA concentrations were 30.18 and 32.45, 36.89, 36.43, 35.66, 35.39 g/L in control and with 2, 4, 6, 8 and 10 % of GC addition, respectively at 96 h. The addition of GC increased the solubilization of substrates, which increased the production of VFAs and other soluble metabolites such as ethanol (Eq. 1), acetic acid (Eq. 2 and Eq. 3), propionic acid (Eq. 4), butyric acid (Eq. 5). Further, the VFAs and alcohol produced during the acidogenic phase get oxidized to acetic acid, as shown in Eq. 6- Eq. 8. As could be noticed from the equations, hydrogen is generated as a byproduct from all the steps of acidogenesis.
C6H12O6 + 2H2O + 2NADH→2CH3CH2OH + 2HCO3−+2NAD++2H2 (ΔG=-234.8 kJ/mol)… Eq. 1
C6H12O6 + 2H2O → 2CH3COOH + 2CO2 + 4H2 (ΔG = -135.6 kJ/mol) … Eq. 2
C6H12O6 + 4H2O + 2NAD+→2CH3COO−+2HCO3−+2NADH + 2H2 + 6H+ (ΔG=-215.7 kJ/mol) … Eq. 3
C6H12O6→CH3COO− + CH3CH2COO−+CO2 + H2 + 2H+ (ΔG = − 287.0 kJ/mol) … Eq. 4
C6H12O6 + 2H2O→CH3CH2CH2COO−+2HCO3−+2H2 + 3H+ (ΔG = − 261.5kJ/mol) … Eq. 5
CH3CH2COO− + 3H2O → CH3COO− + H+ + HCO3− + 3H2 (ΔG = 76.1 kJ/mol) … Eq. 6
CH3CH2CH2COO−+2H2O → 2CH3COO−+H++2H2 (ΔG = 48.1 kJ/mol) … Eq. 7
CH3CH2OH + 2H2O → CH3COO− + 2H2 + H+ (ΔG = 9.6 kJ/mol) … Eq. 8
Subsequently, profiling of VFA metabolites was carried out using gas chromatography with a flame ionization detector. While ethanol production increased with time in the control, its concentration was many folds higher in the digester bottles with GC. As depicted in Fig. 2, the maximum concentration of ethanol and its respective fraction out of the total soluble products were 1.71 g/L (35.80 %) at 72 h in control, and 8.51 g/L (50.12 %), 11.56 g/L (89.88 %) and 12.31 g/L (82.05 %), respectively, in digester bottles with 6, 8 and 10 % of GC at 24 h. With 2 and 4 % GC addition, the ethanol production was lower, 1.68 g/L (36.07 %) at 48 h and 4.03 g/L (41.68 %) at 72 h, respectively. The ethanol concentration decreased after 72 h in all the digester bottles and control, mainly because of the increase in pH. In the study done by Wainaina et al. (2019), lower pH (~ 3.2–5.0) led to alcohol biosynthesis and the generation of lactic and acetic acids. Similarly, in our experiment, the pH was 3.2–4.6, which was the main reason for ethanol production.
Acetic acid was dominant at various points of time in the digester bottles with and without addition of GC: 2 % (24 h), 4 % (24 h), 6 % (96 h), 8 % (72 h), 10 % (96 h) and control (96 h) contributing 8.8, 3.10, 7.57, 0.92, 1.15 and 2.39 g/L, respectively. As depicted in Fig. 2, it can be concluded that with the increased addition of GC, acetic acid production got delayed and decreased because lignocellulosic biomass, such as GC, affect hydrolysis, the rate-limiting step of the AD process (Chakraborty and Venkata Mohan 2018; Ferdes et al. 2020; Pilarski et al. 2020). Nonetheless, the addition of 2 to 6 % of GC was most favourable for the timely production of appropriate volumes of acetic acid. The concentration of propionic acid was 4.05 g/L in the digester bottles with 2 % GC. In the rest of the digester bottles and control, the propionic acid production was < 1 g/L (Fig. 2).
It is reported in the literature that the production of VFA, including acetic and propionic acids, is maximum at pH of 6 (Dahiya et al. 2015; Frohlich-Wyder et al. 2017; Wainaina et al. 2019). Additionally, during the redox phase of AD, different groups of acidogenic microorganisms would grow well only in their respective microenvironment with specific pH. To enhance the production of a particular VFA by the acidogens, it is necessary to maintain the pH more than the pKa value of that VFA. If the pH value goes below the pKa, then the VFA accumulates in the cytosol and hinders microbial growth, which will decrease the VFA production (Sarkar et al. 2021). As acetic acid is the main product of acidogenesis, its impact on the redox environment will be more significant than the other VFAs. Correspondingly, lower pH (~ 4.6) yielded lower proportions of acetic and propionic acids in the reactors in this investigation.
The maximum butyric acid production was at 96 h in all the digester bottles with GC, but it was at 24 h in control. The butyric acid concentrations were: 1.6, 2.36, 2.25, 2.51 and 2.36 g/L for digester bottles with 2, 4, 6, 8 % GC addition and control, respectively, and absent at 10 % GC (Fig. 2). Similar to acetic acid production, which happened during the initial hours of the experiment, butyric acid production also got delayed and could be measured only at 72 and 96 h in the digester bottles with GC. Results of this study are aligned with the reports of Sarkar et al. (2016) and Karthikeyan et al. (2016a), which explain the dominance of acetate-butyrate and mixed acid pathways that augment acidogenesis and enhances the production of acetate, butyrate, ethanol, H2 and CO2. Eventually, butyrate gets converted to acetate and hydrogen; hence the concentrations of acetate and butyrate are inversely proportional to each other, in all the digester bottles with GC addition, except 10 %. The main reason for no butyrate production at 10 % GC is higher ethanol production (solventogenesis). Long-chain fatty acids (e.g., caproic acid) are generated through the β-oxidation pathway, where ethanol acts as electron donor and short-chain fatty acids (e.g., butyric acid) act as the electron acceptors (Wu et al. 2018).
Valeric (C5) and caproic (C6) acids were produced in all the digester bottles with GC, but isovaleric and isobutyric acids production was negligible. The maximum valeric acid concentrations were 1.83, 1.77, 2.02, 1.3, 2.25 and 2.34 g/L for control and the digester bottles with 2, 4, 6, 8 and 10 % addition of GC, respectively (Fig. 2). Caproic acid production increased with the addition of GC, but it was less than 1 g/L. The variation in the long-chain fatty acids was due to the acidic environment and the availability of ethanol, which promoted the β-oxidation pathway to produce long-chain fatty acids (Steinbusch et al. 2011).
3.4 Biogas production and volatile solids removal
A focus was given on the regulation of the redox environment by GC during the initial period (up to 96 h) of FW anaerobic co-digestion. Hence this investigation centred around acidogenesis, which primarily produces VFA and hydrogen (Eqs. 1 to 8) than methane. Gas chromatography analysis of biogas revealed that the maximum concentrations of hydrogen were ≈ 30 % and ≈ 25 % in experiments with GC and control, respectively, at 48 h. The specific hydrogen yield (13.56 ml/gVS, 29.68 %) was maximum in 6 % addition of GC, which was higher than 4% addition of GC (10.69 ml/gVS, 27 %), 2 % GC (10.36 ml/gVS, 27 %), 8 % GC (10.21 ml/gVS, 26.2 %), 10 % GC (10.27 ml/gVS, 26.1 %) and control (8.63 ml/gVS, 25 %) (Fig. 1c, Table 2). It has to be noted that the difference in hydrogen yields between the experiments with 2 and 6 % of GC addition is 0.82 ml/gVS. Hydrogen conversion efficiency was maximum in 6 % addition of GC (48.73 %) as compared to 2 % (47.51 %), 4 % (47.11 %), 8 % (46.26 %), 10 % addition of GC (45.51 %) and control (44.94 %).
Higher hydrogen production in the experiments with GC could be attributed to the higher production of acetate, butyrate and propionate in these digester bottles compared to control. Additionally, from the thermodynamics of anaerobic digestion, it was evident that the acetate formation pathway from Acetyl-CoA was favoured (Eqs. 6–8), which also led to hydrogen generation (Taheri et al. 2018). However, methane production was less till 96 h, mainly due to the inhibition of methanogens by the acidic environment (Han et al. 2019). During AD of FW, methane generation begins typically after 7 to 9 days from the start of the digestion due to the slow growth rate of methanogens (Chakraborty and Venkata Mohan 2018; Maus et al. 2020). For readily biodegradable substrates such as FW, accumulation of VFA and alcohols due to rapid hydrolysis and subsequent inhibition of acetogens and methanogens is a common phenomenon reported in several studies (Karthikeyan et al. 2016a; Karthikeyan et al. 2016b; Chakraborty et al. 2018b, Han et al. 2019). Additionally, due to its particulate nature, the addition of GC at higher concentrations will alter hydrolysis and VFA dynamics, which eventually hamper methanogenesis, as reported by Ferdes et al. (2020).
Compared to control, the VS removal and cumulative yields of ethanol and VFA were higher in all digester bottles with the addition of GC as co-substrate (Table 2). Digester bottles with 6 % of GC addition achieved the highest VS removal (47.37 %). The cumulative yields of ethanol (0.613 g/gVS), propionate (0.209 g/gVS) and caproate (0.064 g/gVS) were highest with 6 % of GC addition. Maximum cumulative yields of acetate (0.468 g/gVS) were with 2 % of GC, while at 4 % of GC, butyrate (0.106 g/gVS) and valerate (0.190 g/gVS) generation were maximum. Considering the VS removal, VFA and hydrogen production, the addition of 6 % of GC is recommended for anaerobic co-digestion of FW in mesophilic conditions.
3.5 VFA kinetics and statistical analyses
The acidogenic process is dynamic which gets influenced by the intermediate metabolites of anaerobic digestion, particularly the concentration and composition of the VFA (Yan 2016). Hence, VFA kinetics was evaluated as presented in Fig. 4. The maximum production of ethanol of 0.51 g/L.h (8 % GC), acetate of 0.37 g/L.h (2 % GC), propionate of 0.17 g/L.h (6 % GC), valerate of 0.08 g/L.h (10 % GC) and caproate of 0.035 g/L.h (8 % GC) were observed at 24 h, whereas for butyrate production of 0.098 g/L.h (6 % GC) was observed at 72 h. The maximum consumption of the metabolites happened at the same concentration of GC, however at different times, i.e., for ethanol − 0.215 g/L.h, acetate − 0.17 g/L.h, propionate − 0.073 g/L.h and valerate − 0.02 g/L.h at 48 h and caproate − 0.0083 g/L.h at 96 h. As seen in Fig. 4, the production and consumption rates of VFA and ethanol were higher in the digester bottles with GC addition than control. It has to be noted from Fig. 4 that the higher acetate, propionate and butyrate production in digester bottles with GC addition is also correlated with hydrogen production. A similar pattern of VFA kinetics was observed in the study by Sarkar and Venkata Mohan (2017).
To confirm the above-described relationship between ethanol and VFA in their production and consumption, a statistical correlation analysis was carried out for the digester bottles with maximum concentration (10 %) of GC and control. As explained in Sect. 3.3 for the digester bottles with higher GC (8 % and 10 %), a negative correlation was observed between the production of alcohol (increased) and acetic, butyric and propionic acids (decreased). The concentration of caproic acid was directly proportional to ethanol and inversely proportional to acetic acid. In control, acetic, butyric and caproic acids increased with reduced ethanol production, which showed a negative correlation. Propionic and butyric acids had a positive and negative correlation with acetic acid, respectively (Supplementary Table 1). It has to be noted that the dynamics of VFA production is directly dependent on the pH of the environment. In turn, the redox environment impacts the thermodynamic potential of biochemical reactions involved in acetogenesis (Eqs. 1 to 8) and the syntrophic relationship between hydrolytic, acidogenic and acetogenic microflora (Yan et al. 2016; Owusu-Agyeman et al. 2020). ANOVA single factor analysis of all experiments and control was expressed through the sum of square, mean square, F-value, and P-value, which revealed that the experiment is significant as it has a P value of ≤ 0.05 (Supplementary Table 2).