Technological and Energetic Aspects of Multi-Component Co-Digestion of the Beverage Industry Wastes and Municipal Sewage Sludge

: In the present study, the co-digestion effectiveness of the selected beverage wastes and municipal sewage sludge in two- and three-component mixtures was evaluated. Orange peels and orange pulp, as well as brewery spent grain were applied as co-substrates to sewage sludge at the following doses: 1.5 and 3.0 g of orange peels, 2.5 and 5 g of orange pulp, and 1.5 g brewery spent grain. Mono-digestion of sewage sludge was used as a control. The experiments were performed under mesophilic conditions in batch reactors. As compared to the control, only in the presence of the highest dose of pulp, brewery spent grain and sewage sludge was the increased methane production of 395 mL CH 4 g − 1 VS accompanying an additional energy proﬁt of 82% observed. Moreover, in this case, the enhanced volatile solids removal and lower accumulation of p-cymene were found. These results were despite the increased limonene and phenol content in the feedstock, conﬁrming a synergistic effect at the highest dose of pulp, brewery spent grain and sewage sludge.


Introduction
Currently, industrial food production generates considerable amounts of by-products, the effective management of which still remains a technological and economical challenge. A significant share of the global market is represented by beverage manufacturing (both alcoholic and non-alcoholic). Despite major technological improvements, this industry is recognised as the most energy-and water-consuming. Moreover, its production is related to substantial CO 2 emissions. Presently, the global beverage market is dominated by the orange juice and beer industry [1][2][3]. When processing orange juice, about 50-60% of the fresh fruit becomes waste, which includes peel, seeds, pulp and membrane residues [4][5][6]. In contrast, about 85% of all by-products generated by the beer industry are represented by brewery spent grain [7]. This solid fraction is composed mainly by the husk and pericarpseed coat that covers the barley grain [8]. Disposal of both solid residues indicates many difficulties for companies. Due to significant moisture content, their properties deteriorate rapidly, making long-term storage impossible; long-distance transport should also be avoided because of the significant costs. Additionally, both types of waste contain a high organic load; for this reason, they cannot be directly discharged into the environment. Therefore, their reuse, recovery and recycling are limited [4,9].
In turn, the afore-mentioned wastes comprise numerous valuable components. Orange waste contains various nutrients, including vitamins A, B and C; minerals (phosphorus, calcium, and potassium); triterpenes; amino acids; phenolic acids and carotenoids [10]. In contrast, brewery spent grain exhibits a significant content of lipids, proteins, minerals to high energy consumption, while the mono-digestion of sewage sludge results in low methane production, as it corresponds to the substrate characteristics-i.e., deficiency of organic matter, imbalanced C/N ratio (carbon to nitrogen ratio) and presence of inhibitors (heavy metals, pharmaceuticals and pathogens). At many WWTPs, the digesters are usually oversized (even up to 30%), which enables introducing an additional substrate without any investments [18]. The implementation of both beverage wastes as co-substrates for sewage sludge anaerobic digestion might improve the feedstock composition by ensuring a significant amount of organic matter. It should be pointed out that such a study has not yet been conducted.
The basic novelty of the work is such a multi-aspect investigation, estimating the effect of another substrate (brewery spent grain) addition on the efficiency of orange waste and sewage sludge co-digestion. It is worth mentioning that many studies conducted in this field have concerned only two-component mixtures. The works involving additional components are rare-although the use of an additional substrate may allow for the efficient anaerobic bioconversion of orange waste in anaerobic co-digestion, with sewage sludge providing an additional energy gain for the WWTPs. Orange waste could increase the C/N ratio, while brewery spent grain could provide the required buffering capacity. Importantly, each of these indicates a significant biogas potential and contains the microelements stimulating the activity of methanogens, and thus contributing to the increase in biogas production. However, anaerobic co-digestion performance and its operational conditions should be constantly monitored due to the presence and effect of inhibitors typical of the co-substrates mentioned.
The aim of this study is to evaluate the effectiveness of the co-digestion of selected beverage wastes and municipal sewage sludge in tertiary mixtures, comparing the results with those for two-substrate configurations. Orange peels and pulp as well as brewery spent grain have been applied as co-substrates. The influence of orange waste and brewery spent grain addition on sewage sludge anaerobic digestion is examined on the basis of biogas/methane yields, kinetics (using four different models) and volatile solids removal (VSR). Additionally, the effects of limonene, phenols and p-cymene are evaluated by Fourier-transform infrared photoacoustic spectroscopy (FT-IR/PAS) analysis. To evaluate the possible energy profits resulting from the implementation of these substrates to the existing digester, the energy analysis is also provided. It should be noted that the proposed technology might be easily implemented on a technical scale without significant financial outlays.

Substrates and Inoculum
The sewage sludge and inoculum were collected from a mechanical-biological municipal WWTP in Lublin (Hajdów), Poland. The inoculum originated as effluent sludge from mesophilic digesters. Immediately after collection it was transported to the laboratory and then it was sieved to remove fibres and particles above 5 mm. Subsequently, it was introduced into batch reactors, where it was stored under anaerobic conditions at 37 • C to diminish the non-specific biogas production. The procedure of inoculum preparation was adopted from a study conducted by Angelidaki et al. [38]. This sample was characterised by total solids (TS) and volatile solids (VS) of 18.85 ± 0.23 and 9.95 ± 0.16 g kg −1 , respectively, as well as a pH of 7.9 ± 0.1.
The sewage sludge was a main substrate; it was prepared as a mixture of primary and waste sludge, both thickened, and taken separately from Hajdów WWTP. Immediately after collection, the samples were sieved and mixed in a volumetric ratio of 60:40 v/v under laboratory conditions. Such a proportion is recommended as beneficial for an effective biogas production from sewage sludge.
In the present study, the following beverage solid residues were used as additives for co-digestion with sewage sludge: orange peels, pulp and brewery spent grain. The first two substrates were prepared under laboratory conditions as by-products from a bench-scale orange juice processor. These were mechanically blended to obtain a particle size of 2-5 mm. Then, these samples were frozen and stored at −25 • C in a laboratory freezer. Before each experiment, they were thawed. Brewery spent grain originated from a craft brewery located in Lublin, Poland. Before supplying the reactor, it was dried at a temperature of 55 • C for at least for two hours in a laboratory dryer. Such preparation of the substrates was adopted to maintained unchanged co-substrate properties during the experiment. This procedure allowed for the elimination of the substrate composition variability's influence on the obtained results. The composition of the sewage sludge and co-substrates was controlled every time while introducing them into the reactors (Table 2).

Laboratory Installation and Operational Set-Up
The batch tests were conducted with the specially designed BioReactor Simulator (Bioprocess Control AB, Lund, Sweden). It compromised 6 reactors, with an active volume of 2.0 L. Each digester was equipped with a mechanical stirrer operating on a cycle of 5 min mixing and 25 min resting. In order to maintain a proper temperature, reactors were kept in a water bath. The biogas production was monitored and recorded automatically by using a wet gas flow-measuring unit that operated on the principle of liquid displacement and buoyancy and was integrated with the system for continuous data acquisition. Before the experiments, all reactors were flushed with inert gas (N 2 ) to achieve anaerobic conditions.
Each bottle was supplied with 1.4 L of inoculum and 0.4 L of sewage sludge. The substrate-to-inoculum ratio was adopted based on a procedure given by Angelidaki et al. [38]. In order to evaluate the co-digestion efficiency, a separate control series (R1) was provided, wherein no co-substrates were added. In R2-R5, the anaerobic co-digestion of sewage sludge and orange waste in a two-component system was studied. In turn, in R6-R8, three-component co-digestion was evaluated. The adopted operational set-up is shown in Table 1. The study was conducted under mesophilic conditions (37 • C) and lasted 21 days. The experiments were repeated three times under unchanged operational conditions; for each one, new samples of inoculum and sewage sludge were taken and prepared in the same manner.

Analytical Methods
Feedstock and digestate composition were monitored by analysing the following parameters: total chemical oxygen demand (COD), TS, VS, total nitrogen (TN) and total phosphorus (TP). The soluble chemical oxygen demand (SCOD), VFA, alkalinity (ALK), ammonia nitrogen (NH 4 + −N) and orthophosphate phosphorus (PO 4 3− −P) were measured in supernatant. The soluble fraction was obtained by centrifugation of the samples at 4000 r min −1 for 30 min. The composition of substrates was examined after their collection. For each one, the analogous parameters, as listed above, were included. TS and VS were determined using the Standard Methods for the Examination of Water and Wastewater [39]. Other analyses were performed using Hach Lange UV-VIS DR 3900 (Hach, Loveland, CO, USA) and s cuvette tests. The pH values were monitored using a HQ 40D Hach Lange multimeter (Hach, Loveland, CO, USA).
In order to examine the biogas composition, a ThermoTrace GC-Ultra gas chromatograph provided by Thermo Fisher Scientific (Milan, Italy) was applied. It was equipped with a conductivity detector fitted with divinylbenzene (DVB) packed columns (RTQ-Bond). The parameters for the analyses were: 50 • C for the injector and 100 • C for the detector; the employed carrier gas was helium, with a flow rate of 1.5 cm 3 min −1 .
Trace GC Ultra PolarisQ (Thermo Electron, Italy) was used to analyse the limonene and p-cymene contents. The parameters of the capillary column were: Supelco (30 m × 0.25 mm ID × 0.25 µm). The spectra were taken under the following conditions: transporter gas helium, with a flow rate of 1.5 cm 3 min −1 ; injection volume of 1 mL. The inlet temperature was initially 40 • C. It was then raised by 5 • C/min to reach the level of 250 • C, where it was kept for 3 min. In order to estimate the peak areas, the CHROM-CARD software was involved.

Kinetic Modelling
The kinetic evaluation was performed using four different models: Hyperbolic, modified Gompertz, cone and logistic growth. Their equations are presented below: • Hyperbolic • Modified Gompertz • Cone • Logistic growth model where: • P(t) is the cumulative methane production (mL CH 4 g VS −1 ); • P m is the methane production (mL CH 4 g −1 VS); • R m is the is the maximum methane production rate (mL CH 4 g −1 VS d −1 ); • e is a constant (2.71828); • λ is the lag phase (d); • n is the dimensionless shape factor; • k is rate constant (d −1 ). To indicate the suitability of the model, the root-mean-square error (RMSE) was applied. It was estimated using the following equation [40]: where: • P exp is the methane production achieved in the experiment (mL CH 4 g −1 VS); • P mod is the methane production achieved in kinetic evaluation (mL CH 4 g −1 VS); • N is a number of measurements.

FT-IR/PAS Analysis
The samples of feed and digestate were collected from the reactors for FT-IR/PAS (Fourier-transform infrared photoacoustic spectroscopy). First, they were dried via sublimation (i.e., lyophilised). The 50 mL samples were frozen at a temperature of −25 • C in the laboratory freezer and subsequently dried at −56 • C by means of an Alpha 1-4 freeze dryer. The FT-IR/PAS spectra of the lyophilised samples were obtained by using a Nicolet 6700 spectrometer and a Gasera PA301 photoacoustic cell in the 3800-700 cm −1 range, resolution 4 cm −1 , at room temperature. Dry helium was utilized to purge the photoacoustic cell prior to the collection of data. By computing the ratio of a sample spectrum to the spectrum of a carbon black standard, the spectra were normalised. A stainless-steel cup with a diameter of 10 mm was filled with the sample that had a thickness of less than 3 mm. Interferograms of 1024 scans were averaged for the spectrum; thus, a good signal-to-noise (S/N) ratio was provided. The smoothing functions were not used. All spectral measurements were performed with at least three repetitions.

Statistical Analysis
All performed analyses were performed in triplicates; the average values and standard deviation are given in tables and figures. For statistical analysis, ANOVA (Shapiro-Wilk's, Levene's and Tukey's tests were included) was performed with Statsoft Statistica software (v 13). The significance level was assumed at p < 0.05. All kinetic constants were established using a non-linear regression method. The strength of the relationships between the results was determined using Pearson's correlation coefficient (R) and determination coefficient (R 2 ).

Characteristics of Substrates and Organic Removal
The characteristics of the applied substrates are presented in Table 2. The used orange wastes indicated similar composition. Both orange peels and pulp were characterised by a low pH, relatively high organic content and comparable nutrient contents. The major differences were observed for TS, VS and phenol and limonene content. It should be pointed out that the orange peels presented an enhanced concentration of inhibitors: phenol and limonene, which might be problematic in anaerobic bioconversion. In turn, brewery spent grain was characterised by the highest organic content as well as nitrogen and phosphorus concentrations. Moreover, to other applied substrates, it indicated an enhanced alkalinity.
The composition of feedstock and digestate is shown in Figure 1. Introducing cosubstrates resulted in an increase in the VS, TS and COD contents as compared to the control (sewage sludge only), and the highest enhancements were observed for VS and TS in three-component anaerobic co-digestion. This observation was related to the enrichment of feedstock with brewery spent grain ( Table 2) by 18, 14 and 16% for VS, as well as 15, 12 and 12.5% for TS in R6, R7 and R8, respectively. Considering COD, in two-component anaerobic co-digestion, the major increase of 7% was found in the presence of both orange pulp doses, while in threecomponent series it was 10% when involving 5 g of orange pulp and 1.5 g of brewery spent grain (R8). Similarly, the highest growths for SCOD of approx. 50% occurred in R7 and R8 at both doses of orange pulp mixed with brewery spent grain. In relation to sewage sludge mono-digestion, the VSR and TS removal were improved every time (Figure 1a,b), although to a lesser extent in the series with brewery spent grain. The highest increases were found in the lowest co-substrates doses (with comparable VS loads to control). An increased organic load in feedstock might result in This observation was related to the enrichment of feedstock with brewery spent grain (Table 2) by 18, 14 and 16% for VS, as well as 15, 12 and 12.5% for TS in R6, R7 and R8, respectively. Considering COD, in two-component anaerobic co-digestion, the major increase of 7% was found in the presence of both orange pulp doses, while in threecomponent series it was 10% when involving 5 g of orange pulp and 1.5 g of brewery spent grain (R8). Similarly, the highest growths for SCOD of approx. 50% occurred in R7 and R8 at both doses of orange pulp mixed with brewery spent grain.
In relation to sewage sludge mono-digestion, the VSR and TS removal were improved every time (Figure 1a,b), although to a lesser extent in the series with brewery spent grain. The highest increases were found in the lowest co-substrates doses (with comparable VS loads to control). An increased organic load in feedstock might result in poor process performance and promote the process instability caused by VFA accumulation. In particular, this effect might occur in anaerobic co-digestion with easily biodegradable substrates such as fruit and other food wastes [41][42][43]. Nevertheless, such an effect was not observed when analysing the VFA in digestate (Table 3). The results in terms of VSR were consistent with other studies. Serrano et al. [44] performed the mesophilic anaerobic co-digestion of sewage sludge and orange peels mixed at a proportion of 70:30 (wet weight), achieving the VSR of 76%; in sewage sludge, the monodigestion was only 53%. Martínez et al. [22] evaluated the anaerobic co-digestion of sewage sludge and orange peels, as well as sewage sludge, orange peels and biochar addition. At HRT of 20 days, VS was removed by 63 and 49% in the two-and three-component systems, respectively.
Considering the COD and SCOD, their removal was also improved as compared to the control, wherein a significant drop referred to the latter and reached 36 and 64% in R7 and R8, respectively. The increases might point to the effective consumption of soluble organics via anaerobic co-digestion. The exceptions were the reactors supplied by the orange peels and brewery spent grain (R6), as well as the highest orange peel dose (R5), wherein a declining trend occurred. It seems to indicate that the process inhibition was caused by the high concentration of p-cymene (Table 3) as a product of limonene conversion. Importantly, in R5 and R6, the highest concentrations of this compound in the feedstock were found.

Inhibitors Presence
To evaluate the influence of the co-substrate application on the process stability, the following parameters were analysed: the pH value, ALK, VFA concentration and the VFA/ALK ratio (Table 3).
It should be noted that all of the analysed parameters reached the values recommended for adequate methanogenic activity [45]. Stable process performance was also confirmed by the VFA/ALK ratio, which reached values much below 0.3, ranging from 0.052 to 0.068 in all anaerobic co-digestion series, with the highest values found in the two-component system.
Due to the possible inhibitory effect, the contents of the phenolic compounds in the feedstock and digestate were also analysed [46]. Their negative influence on the AD process is connected to the damage of microbial cells through changing the membrane permeability, which might lead to the leakage of intra-cellular components and the inactivation of essential enzymatic systems [47]. However, its toxicity is dependent upon several factors, such as: concentration, micro-organisms' consortium and acclimatization period [48].
The application of orange pulp and orange peels increased the content of phenols in the feedstock; however, the difference was of no statistical significance. The significant changes were found in the three-component anaerobic co-digestion. In R6, the concentration was enhanced twice, while in R7 and R8, more than three-fold increases were noticed. It was not surprising due to the composition of substrates used and their doses (Tables 1 and 2).
Compared to the feedstock compositions in R1, R3 and R4, a slight increase in the phenol content in the digestate was found. In the remaining series, the concentration of phenols decreased, indicating their metabolization through anaerobic co-digestion. Importantly, the major removal efficiency of approx. 57% was noted in the anaerobic co-digestion of orange pulp and brewery spent grain. Some studies have indicated that the presence of easily biodegradable substrates might improve the degradation of recalcitrant compounds, such as phenol; this effect is related to the occurring co-metabolism and presence of enzymes involved in these mechanisms [49][50][51]. This observation also confirmed the synergistic effect occurring between substrates.
In the present study, the limonene and p-cymene contents were also analysed. This was due to the possible inhibitory effect of limonene, which is related to the occurrence of cymene which is generated during biotransformation of the former in AD [21]. Inhibition is expected because essential oils might accumulate in the cell membrane and change its fluidity, resulting in increased permeability and finally lead to leakage of the cell contents [52].
The application of both orange wastes resulted in a significant increase in the limonene concentration, as compared to sewage sludge. There were 3.8, 4.7, 10.7, 21.1, 16.1, 9.9 and 13.3-fold increases in feedstock for R2, R3, R4, R5, R6, R7 and R8, respectively-the greatest increase being in the presence of orange peels. The major increment was obviously related to the application of the highest orange peel dose. In turn, the p-cymene content was at a low level. The introduction of co-substrates slightly increased its concentration in the feedstock, with the exception of R2. As shown in Table 3, in all series, limonene was efficiently degraded in the AD process. The related removal efficiencies were 67, 85, 85, 94, 95, 97, 93 and 96% in R1, R2, R3, R4, R5, R6, R7 and R8, respectively. Simultaneously, a significant increase in the p-cymene content of the digestate was noted. The exception was a control series fed with sewage sludge, wherein its removal was observed. Furthermore, the major increases indicating the accumulation of p-cymene through anaerobic co-digestion were observed in the presence of orange peels, and they corresponded with both the high content of limonene in the feedstock (R4, R5 and R6) and the orange peel dose.
In three-component anaerobic co-digestion, using orange pulp and brewery spent grain as co-substrates (R7 and R8), and despite a significant content of limonene in the feedstock, a lower accumulation of p-cymene has been recorded. This seems to indicate the greater flexibility of such a configuration and its usefulness in metabolizing the by-products generated through the biotransformation of limonene. The synergistic effect in anaerobic co-digestion of orange pulp, brewery spent grain and sewage sludge was confirmed. The obtained results are consistent with other studies that suggested the co-digestion strategy for reducing the toxicity of individual substrates [53]. The previous studies have indicated that anaerobic micro-organisms exhibit some adaptability to the inhibitory effect of essential oils, while its extent depends on various factors such as the type of substrate and kind of inoculum, as well as the temperature [54].
The beneficial effect observed in the batch system needs to be confirmed each time in the semi-flow configuration, where operational factors (i.e., organic loading rate and hydraulic retention time) significantly affect the anaerobic co-digestion performance. In such systems, the process inhibition might occur more rapidly due to the continuous supply of the substrate, the accumulation of toxicants and insufficient biomass adaptation. Accordingly, studies in future would continue to involve a semi-flow mode.

Biogas and Methane Productions
The results of biogas and methane production are presented in Figure 2 and Table 4. As is shown in Figure 2, in all reactors, the peak related to decomposition of carbohydrates occurred within the first three days of the experiment. The highest initial growth was found at both orange peel doses. As compared to the control, the decrease in daily biogas and methane yields were observed only in the presence of the lowest dose of orange pulp (R2), with the subsequent enhancement between 4 and 10 days of the experiment. Moreover, between 17 and 20 days, a minor additional peak was noticed in the cases of R3, R4 and R5.
Such fluctuations might result from the decomposition of the more complex compounds such as proteins and lignocellulose presented in the co-substrates [55]. occurred within the first three days of the experiment. The highest initial growth was found at both orange peel doses. As compared to the control, the decrease in daily biogas and methane yields were observed only in the presence of the lowest dose of orange pulp (R2), with the subsequent enhancement between 4 and 10 days of the experiment. Moreover, between 17 and 20 days, a minor additional peak was noticed in the cases of R3, R4 and R5. Such fluctuations might result from the decomposition of the more complex compounds such as proteins and lignocellulose presented in the co-substrates [55]. The cumulative biogas production was improved in almost all of the co-digestion series. The statistically significant increments of 23, 20 and 18% were achieved in R4, R5 and R8, respectively, as compared to sewage sludge mono-digestion. Moreover, it should be pointed out that in R7, supplied by a lower dose of orange pulp and brewery spent grain, the enhancement of 12% occurred as compared to the control reactor (R1); however  The cumulative biogas production was improved in almost all of the co-digestion series. The statistically significant increments of 23, 20 and 18% were achieved in R4, R5 and R8, respectively, as compared to sewage sludge mono-digestion. Moreover, it should be pointed out that in R7, supplied by a lower dose of orange pulp and brewery spent grain, the enhancement of 12% occurred as compared to the control reactor (R1); however the observed change was of no statistical significance. The achieved increases were attributed to the introduction of co-substrates rich in organic matter, improving C/N and providing deficient micro-and macro-elements. This was favourable for the anaerobic co-digestion performance, which was indicated by the enhanced values of the biogas production rate (GPR) ( Table 4). Ensuring optimal conditions for anaerobic micro-organisms contributed to enhancing the microbial activity, resulting in increased biogas production.
A different trend was observed in the tertiary mixtures in the presence of sewage sludge, orange peels and brewery spent grain (R6); therein, a comparable result to sewage sludge mono-digestion was found. This tendency occurred despite improved feedstock composition, as well as high VSR. It might indicate process inhibition related to the lowest pH in the feedstock and the presence of both limonene and phenol (Table 3). Moreover, in this case, micro-organisms were additionally involved in the transformation of complex compounds present in brewery spent grain, also contributing to the lower biogas production. This seems to be confirmed by the GPR value, which was the lowest among all the experiments ( Table 4). The explanation might be the fact that cellulose, hemicellulose and lignin, as well as the interactions of these components, create a highly resistant structure for AD microbes [25]. Importantly, in the presence of both orange pulp and brewery spent grain, there was no adverse effect observed, which testifies to their proper selection for anaerobic co-digestion with sewage sludge.
It is not a surprise that the application of carbohydrate-rich substrates resulted in a reduction of the methane content in the biogas (Table 4). This obviously influenced the results in terms of methane production ( Figure 2).
As compared to sewage sludge mono-digestion, an improvement of 12% was observed only in the presence of the highest dose of orange pulp and brewery spent grain (R8), accompanied by the highest MPR value. In turn, a major, statistically significant drop of 26% was found in three-component anaerobic co-digestion with orange peels and brewery spent grain, also confirming the process inhibition visible in the lowest MPR and cumulative methane production (Figure 2).
The results in terms of biogas/methane production, as well as the process disturbances, correspond to the previous studies referring to the application of citrus wastes in the anaerobic digestion process. The methane potential of various citrus peels varied between 353.9 and 398.2 mL CH 4 g −1 VS [21]. However, in this case, process inhibition was observed for a limonene content higher than 200 mg kg −1 . In turn, the AD of pre-treated orange peels by steam distillation resulted in the production of 230 mL CH 4 g −1 VS. The higher values of 400-600 mL CH 4 g −1 VS were observed under thermophilic conditions [6]. In a semi-flow system, Serrano et al. [44] evaluated the mesophilic anaerobic co-digestion of sewage sludge and orange peels mixed at a proportion of 70:30 (wet weight). The achieved methane yield was 165 mL CH 4 g −1 VS, and the VSR reached 76%. Nevertheless, in this case, at the highest organic loads, the methane production decreased due to ammonia inhibition. Martinez et al. [22] supplemented the two-component mixture with biochar to overcome the limitation of the anaerobic degradation of sewage sludge and citrus peels.
The enhanced values of 500 and 704 mL CH 4 g −1 VS were observed in the presence of 10 and 30 g of biochar, respectively. The methane potential of sewage sludge and orange peels mixed at a VS ratio of 1:1 was only 298 mL CH 4 g −1 VS. The increases in biochar presence were explained by the adsorption of inhibitors, as well as provided the high surface area that promoted the adhesion and growth of micro-organisms.

Kinetic Evaluation
According to Kainthola et al. [56], to analyse the results of a kinetic study, the percentage error between the experimental and predicted values should not exceed 25%. As was shown in Table 4, this condition has been fulfilled for all the applied models. Moreover, the results indicated that the best fitted model for selected substrates were the modified Gompertz, cone and logistic growth models; this fact was also confirmed by high R2 values as well as a low value of RMSE. Previous studies indicated that the modified Gompertz showed the best fit for anaerobic co-digestion systems [57,58]. In turn, the cone and logistic growth models have been frequently applied for various food wastes [58][59][60][61].
In both the hyperbolic and cone models, the k constant was improved in almost all of the anaerobic co-digestion series as compared to control, indicating enhanced rates of degradation in the presence of co-substrates [58,62]. The highest k values were found in the presence of both doses of orange peels and sewage sludge (R4 and R5), corresponding to increased biogas production. However, in R8, supplied by 5 g of orange pulp and 1.5 g of brewery spent grain, a comparable result to the control was found. This fact might be related to both the introduction of hardly biodegradable brewery spent grain with a significant lignin content and the increased dose of orange pulp [59,63]. A similar trend, indicating a drop in k during the anaerobic digestion of fruit residues with significant lignin content, was also observed by Zhao et al. [59].
In R2, fed by sewage sludge and the lowest dose of orange pulp, a decreased k was observed. This seems to confirm to the process inhibition caused by toxic compounds. It should be pointed out that the obtained k constants were higher than for those from other studies [64,65]. This is probably due to the fact that the value of k is directly related to the feedstock characteristic, as well as to the adopted process conditions (mono-or co-digestion, temperature, substrate doses and type of inoculum). It typically varies between a wide range of 0.02-0.49 d −1 [66]. Moreover, the higher k values were found in the cone model, though the same observation was found in the study performed by Zhao et al. [59].
While considering R8, there was no concordance between the tendency observed for the experimental data (MPR) and the k constant from the models listed above. It was shown that in the experiment an increased MPR value was found as compared to control. This seems to indicate some inadequacy in using such models for a complex, three-component feedstock of sewage sludge, orange pulp and brewery spent grain. It is commonly known that models are usually well verified for mono-substrates with strictly defined composition, although they might be inadequate for describing the complexity of the co-digestion of sewage sludge with other waste.
Regarding the maximum methane production rate, the increased values were found in the presence of both doses of orange peels (R4 and R5), as well as for orange pulp and brewery spent grain (R7 and R8) in the two-and three-component systems, respectively. Importantly, as compared to control, the decrease in this parameter occurred for both doses of pulp (R2 and R3) and for orange peels and brewery spent grain (R6), which is consistent with the trend observed for MPR (Table 4).
Lag phase describes the microbial adaptation to the substrates. In both modified Gompertz and logistic growth models, this parameter reached values below zero, indicating a sufficient acclimatization of micro-organisms to degrade the substrates and a rapid start to the reaction. The results also indicated the absence of the lag phase. Similar findings were reported by Zhao et al. [59] and Liew et al. [66]. The lowest values of λ corresponded to the highest dose of easily biodegradable orange pulp (R3). In contrast, the increased λ was obtained in the presence of co-substrates with a complex structure, such as brewery spent grain.
The values of the shape factor in the cone model were comparable to the results reported in literature (1.3-3.6) [58,60,64,67]. The increased values occurred in the presence of both the pulp and brewery spent grain co-substrates.

FT-IR/PAS Analysis
The FT-IR/PAS analysis was used to evaluate both the feed and digestate quality and provided insight into the transformations of organic matter resulting from anaerobic co-digestion.
The FT-IR/PAS spectra of both the raw and digested sewage sludge, as well as the spectra of the sewage sludge with the addition of pulp, orange pulp and brewery spent grain (both raw and digested) are shown in Figure 3a. The intensities of the IR bands are lower in the spectra of the digested samples than in the raw ones. The smallest differences were observed in the raw and digested sewage sludge spectra without an added co-substrate (Figure 3a). The bands indicating the presence of nitrile -C=N and isonitrile -N=C bands (2135 and 2062 cm -1 , respectively) almost completely disappeared for the digested samples. The intensity changes of bands within 2956-2852 cm -1 (C-H stretching) also indicate the decomposition processes of organic components and can be used to evaluate the anaerobic digestion efficiency [68].
The addition of orange pulp and orange peels causes the appearance of new bands in the FT-IR/PAS spectra (Figure 3, R2-F-R8-F) and increases the intensity of other bands. Broad bands with maxima at 3400 and~3340 cm −1 are typical of the alcohol hydroxyl. The band at 1733 cm -1 (C=O stretching) is more intense and sharper than in the spectrum of sewage sludge (Figure 3, R1-F). A similar situation occurs in the case of vibration of the C-H groups (2956-2852 cm -1 ). The band at~1572 cm −1 is due to the presence of a C=O stretching vibration in flavonoids [69], aromatic C=C [70], COOand/or amide II. This band may also indicate the presence of C=C vibration in the terminal methylene group of limonene, which may be confirmed by the presence of the band at 898 cm −1 (out-of-plane bending of the terminal methylene group) [71]. The band at 1460 cm −1 , as well as the band at 1412 cm −1 , may have contributions from more than one group vibration: aliphatic C-H, aromatic C=C and C-O in alcohols. The first of the mentioned bands is characteristic for limonene and p-cymene [72,73]. The band at 1373 cm −1 is probably due to the combination of the O-H deformation vibration and C-O stretching vibration in the phenolic compounds. The band at 1203 cm −1 is due to C-O stretching vibration in alcohols, as well as the bands at 1105 and 1052 cm −1 [69]. The band at 1052 cm −1 may also correspond to C-O-H or C-O-R in alcohols or esters [74].  The addition of orange pulp and orange peels causes the appearance of new bands in the FT-IR/PAS spectra (Figure 3, R2-F-R8-F) and increases the intensity of other bands. Broad bands with maxima at 3400 and ~3340 cm −1 are typical of the alcohol hydroxyl. The band at 1733 cm -1 (C=O stretching) is more intense and sharper than in the spectrum of sewage sludge (Figure 3, R1-F). A similar situation occurs in the case of vibration of the C-H groups (2956-2852 cm -1 ). The band at ~1572 cm −1 is due to the presence of a C=O stretching vibration in flavonoids [69], aromatic C=C [70], COOand/or amide II. This band may also indicate the presence of C=C vibration in the terminal methylene group of limonene, which may be confirmed by the presence of the band at 898 cm −1 (out-of-plane bending of the terminal methylene group) [71]. The band at 1460 cm −1 , as well as the band at 1412 cm −1 , may have contributions from more than one group vibration: aliphatic C-H, aromatic C=C and C-O in alcohols. The first of the mentioned bands is characteristic for limonene and p-cymene [72,73]. The band at 1373 cm −1 is probably due to the combination of the O-H deformation vibration and C-O stretching vibration in the phenolic compounds. The band at 1203 cm −1 is due to C-O stretching vibration in alcohols, as well The decrease in intensity of the bands in the FT-IR/PAS spectra may be evidence of the reduction of organic compounds in the digestate. In the case of the R1-F and R1-D samples (sewage sludge without co-addition before and after anaerobic digestion, respectively), these changes are not as significant as in the cases of the subsequent samples, although the decrease in peak intensities is observable for O-H stretching of hydroxyl and =C-H in aromatic structures groups (3600-3080 cm −1 ), C-H and -OH of COOH in hydroxyl or amine compounds (1412 cm −1 ), C=C in aromatic rings, C-H and C-O-C in ethers (1254-1226 cm −1 ) and C-O-C (1074-1034 cm −1 ). As mentioned before, the addition of co-substrates changes the spectral characteristics. The additional bands visible in the spectrum of the R2-F-R8-F samples disappear completely for digested ones (R2-D-R8-D). Additionally, the intensities of the other bands also decrease. In the R5-F and R6-F spectra, the flavonoids, limonene and p-cymene bands (1572, 1538, 1460, 898 cm −1 ) are the most intense, which suggests that their highest contents are in these samples. These bands, but of lower intensity, are also visible in the spectra of other samples with the pulp and orange peels and with brewery spent grain co-additions. However, the addition of the latter does not significantly affect the shape of the IR spectrum. The absence of the mentioned bands in the spectra of the samples after digestion (R3-D-R8-D) suggests the decomposition of flavonoids, p-cymene and limonene during the AD process. The band at 1572 cm −1 and 898 cm −1 disappear completely, same as the bands of the hydrogen-bonded -OH groups in alcohols (3400, 3340, 1203, 1105 and 1052 cm −1 ), while the bands of O-H deformation vibration and C-O stretching vibration in phenolic compounds (1373 cm −1 ) have much lower intensity. In the spectra of all digested samples, except R1 (Figure 3a), a decrease in the intensity of the aliphatic C-H groups bands can be observed within the 2956-2852 cm −1 range.
The decrease in the intensity of all bands, and thus the highest degree of decomposition, can be observed for the three-component systems (R6-R8, Figure 3f-h). These observations were consistent with the results presented and discussed in the previous sections.

Energy Balance
To estimate the potential energy profits in full-scale facilities, the energy balance for the most favourable two-and three-component series (R5 and R8) was performed. This evaluation was performed on the basis of the procedure described by Szaja and Montusiewicz [32]. The detailed calculations are presented in Table 5. The evaluation was made for an existing digester with a volume of 2500 m 3 , operating at an HRT of 20 days under mesophilic conditions (35 • C) (located at the WWTP in Lublin, Poland).
In co-digestion series R5 and R8, the obtained thermal energy would completely cover the digester energetic demand necessary for heating the feedstock as well as the heat losses. Importantly, in both cases, an additional energy surplus was found. As compared to the control, the energy profit was improved by 15.4 and 81.2% in R5 and R8, respectively. Moreover, the related thermal/electric power production was enhanced by 5.4 and 23.3%. The beneficial results for R8 were attributed to the synergistic effect of the substrates used; thus, significantly higher methane production accompanied the increased VS load in the feedstock. The results in terms of energy balance, regarding the implementation of this technology in an existing digester, are particularly important. Such a solution might provide energy neutrality to WWTPs. However, consideration should be taken regarding the seasonality of these wastes or, eventually, the necessity of their storage, as well as their transport cost.

Conclusions
The conducted research has investigated the effectiveness of the co-digestion of selected beverage wastes and municipal sewage sludge in two and tertiary mixtures. As compared to the control, the increased methane production was achieved only in threecomponent co-digestion in the presence of the highest dose of orange pulp and brewery spent grain (R8). This was accompanied by an increased VSR, comparable kinetics and a lower accumulation of p-cymene. Importantly, therein was the essential improvement of energy balance by 81.2% as compared to control also observed. These were achieved despite the increased limonene and phenol content in the feedstock, confirming that a synergistic effect occurred in this case. Moreover, the FT-IR/PAS analysis also seems to indicate such a three-component anaerobic co-digestion as a beneficial choice.