Effect of Variation in Co-Digestion Ratios of Matooke, Cassava and Sweet Potato Peels on Hydraulic Retention Time, Methane Yield and Its Kinetics

This paper presents the results of batch anaerobic co-digestion of matooke, cassava, and sweet potato peels and vines. These agricultural wastes and others form the biggest portion of household wastes in developing countries. However, they have remained an unexploited resource amidst the ever increasing needs of clean energy and waste disposal challenges. Efforts to use them individually as biogas substrates have been associated with process acidification failure resulting from their fast hydrolysis. The aim of this work was to exploit agricultural wastes is co-digestion among themselves and assess their effect on methane yield and its kinetics, pH and hydraulic retention time (HRT). Sixteen ratios of Matooke peels (MP), cassava peels (CP) and sweet potato peels (SP) were assessed in duplicate. Methane yield and its kinetics, pH and HRT demonstrated dependence on the proportion of substrates in the mixture. Depending on the ratio mixture, HRT increased to 15 days compared to less than 5 days for single substrates, hydrolysis rate constant (k) reduced to a range of 0.1 0.3 d compared to single substrates whose k-values were above 0.5 d, pH was maintained in the range of 6.38 6.43 and CH4 yield increased by 15% 200%. Ratios 2:1:0, 2:0:1, 0:1:2, 1:1:1 and 1:1:4 were consistent all through in terms of model fitting, having a positive synergetic effect on HRT, hydrolysis rate constant, lag phase and methane yield. However, more research is needed in maintaining the pH near the neutral for process stability assurance if household wastes are to be used as standalone substrates for biogas production without being co-substrates with livestock manure.


Introduction
Clean energy shortage and waste disposal challenges at household level in developing countries are a reality and solutions to them are not likely to be realised in the near future amidst the ever growing population. Biogas a clean and sustainable energy source has been developed to serve in part as a waste management mechanism. Nevertheless, biogas technology dissemination to households has been based on livestock manure, neglecting other household wastes [1], such as matooke peels, cassava peels and sweet potato peels among others, which are widely distributed in most of the developing countries. According to FAOSTAT [2], over 90% of the developing countries produce at least two of the three crops (banana, cassava and sweet potatoes) from which the peels are produced. These peels form the highest composition of household wastes [3] as well as municipal solid waste (MSW) given that over 60% of MSW comes from households [4] [5]. Because of the reliance on livestock manure, biogas impact in terms of adoption and use has not been felt since majority of the households do not have livestock. Unlike livestock manure that already contains anaerobes, household wastes and in particular of agricultural origin like matooke peels, cassava peels and sweet potato peels among others are not straightforward substrates to deal with in biogas production. They are liable to acidifying the anaerobic digestion process due to their fast degradability at the hydrolysis stage [6] [7]. In addition, washout of microorganisms is likely to occur due to their short solids retention time/hydraulic retention time associated to their quick degradability [8]. As a result of this, food waste and other household related organic wastes have always been used as co-substrates in most of research and existing working biogas projects [9] [10] [11] [12] [13]. Households however, all they generate is food waste or other agricultural organic wastes like matooke peels, cassava peels and sweet potato peels to mention but a few. Yet, little has been done to understand the performance of household wastes on biogas production process if co-digested between themselves. Co-digestion is said to influence a number of parameters which include among others hydraulic retention time, methane yield and process stability [1].
On the other hand, simple First Order Models that have been widely used in the kinetics of methane production can help in understanding the rates of digestion process [14]. This is not to mention that studies on the kinetics of methane production from the co-digestion of household agricultural wastes amongst themselves were not found in literature. Therefore, the aim of this work was to assess the effect of co-digestion ratio variations of matooke, cassava and sweet potato peels on hydraulic retention time, hydrolysis rate constant and acclimation period as well as the overall methane yield. By doing this, a substrate ratio(s) working in a stable state with controlled hydrolysis rate would be obtained to enable use of household wastes as standalone substrates for biogas production.

Sample Preparation and Composition Analysis
Samples of matooke peels (MP), cassava peels (CP), sweet potato peels (SP*), sweet po-tato vines (SV) were secured from the local markets in Kampala. Samples were dried and ground according to Sluiter, Hames [15]. After grinding, SP* and SV were mixed together in a ratio of 2:1 by weight to form SP which was used in combination with MP and CP. The ratio of 2:1 was based on their weights after drying and anticipated protein content. This was to enhance SP* with proteins expected from SV [16]. Augmented Simplex Lattice Design (ASLD) in Design expert 7 version was used in coming up with 16 MP:CP:SP ratio combinations in addition to SP* and SV standalone samples as presented in Table 1 and Table 2. Total Solids (TS), Volatile Solids (VS) and ash content were determined using Standard Gravimetric Method, APHA-2540D, with an oven set 103˚C -105˚C and muffle furnace at 550˚C for overnight [17]. Crude fat was extracted from a dried sample using boiling diethyl ether in a Soxhlet Apparatus (AOAC-948. 22,1990). Total and Organic Carbon and Nitrogen in % were determined using a VARIO Elementar configured for % C and N. Protein content was estimated by multiply the Nitrogen content by 6.25. Cellulose, hemicellulose and lignin contents were determined by analysis of Neutral Detergent Fibre (NDF), Acid Detergent Fibre (ADF) and Acid Detergent Lignin (ADL) following Van Soest method. That is, cellulose was computed as the difference between ADF and ADL while hemicellulose was calculated as the difference between NDF and ADF. Degradable carbohydrates were estimated by subtracting protein content, fat content and cellulose content from VS. Theoretical methane yield was determined by the composition of organic fraction method as in Equation (1), due to its ability of taking into account the easily biodegradable compounds like carbohydrates, fats and proteins as well as the poorly biodegradable compounds like cellulose [18]. × × × theo M = 415 %carbohydrates + 496 %proteins +1014 %lipids (1) where M theo is the theoretical methane yield and lipids where equated to fats.
Sample solutions for metal analysis were prepared by dry ashing standard method [19].

Inoculum and Batch Experiment Description
The Inoculum was obtained from a 0.5 m 3 household biogas digester fed with household organic waste found in Pitsane village, which is located approximately 100 km South of Gaborone City-Botswana. Inoculum pH was first measured and was found to be 7.72 before it was kept under ambient temperature for 5 days to maximize digestibility of the remaining VS. Its TS and VS were then determined and found to be 1.52% and 65.15% respectively. Inoculum pH was again measured and a value of 7.45 recorded before its use in batch digesters to ensure that it was in optimum range of 6.1 -8.5 [20].
For all the substrate ratios whose characteristics are presented in Table 1 and Table 2, approximately 0.45 g of each was weighed in to a 250 ml serum bottle. The 0.45 g of substrate was based on the VS in the inoculum and in the substrate ratios to give an inoculum to substrate ratio of 2:1. Then 80 ml of inoculum was added followed by 70 ml of clean water to make a working volume of 150 ml. The initial pH value was then measured and recorded as pH start . This was done in duplicate for each substrate ratio with limited space in the incubator being the limiting factor for more replicates. Two control serum bottles each with 120 ml of inoculum and 30 ml of water only were also prepared. All the bottles with their content were placed in an incubator set at 37˚C to run for 35 days. A data harvester of Easysense Vision model, equipped with a K-type thermocouple and a pressure sensor was set up to record any temperature fluctuations in the incubator and atmospheric pressure at an interval of 1 hour over the entire 35 days. This was to help during gas normalization to standard temperature and pressure (STP) conditions and to detect, if any, the effect of temperature changes on the batch experiment. After 35 days, pH was again recorded as pH end to determine how stable the anaerobic digestion process in each of the digesters was. Percentage reduction in VS was determined for each of the substrate ratios representing degradability using Equa- where, VS s is weight of VS fed into the digesters at the start (g), VS d is weight of VS remaining after digestion (g), VS I is weight of VS remaining from the inoculum (g).

Biogas and Methane Measurements
Two separate measurements were performed interchangeably at separate intervals. On the first measurement, the volume of biogas was measured by downward displacement of NaCl solution to minimize CO 2 diffusing/dissolving. The second measurement after an interval of six hours, biogas was passed through NaOH solution allowing all CO 2 to be absorbed. The remaining gas that displaced NaOH solution was assumed to be CH 4 [21]. This interchanging measurement of biogas volume and CH 4  Since both biogas and methane were collected over a liquid, they were normalized to STP conditions as in Equation (3) [23]. To correct for the gas that was generated from the inoculum, Equation (4) was used.
where, V o is normalized volume of dry gas (Nml), V w is wet volume of gas as measured over the liquid (ml), P a is ambient pressure in the room at the time of gas measurement (Pa), P w is vapour pressure of water as a function of ambient temperature in the incubator T o is standard temperature (273 K) T a is ambient temperature in the incubator (K), P o is standard pressure (Pa), V oc is the corrected normalized cumulative methane yield (Nml), is the normalized cumulative methane yield inclusive of methane from the inoculum (Nml), I s is the quantity of inoculum added in the substrate (ml), I c is the quantity of inoculum in the control (ml), t is the digestion time (days), i is the number of days the experiment takes.

Kinetic Models
Three models were selected based on their parameters to predict hydrolysis rate constant and lag phase plus their accuracy in predicting parameters in literature [22] [24] [25]. The Modified Gompertz Model (MGM), Exponential model and the Cone model in Equations (5), (6) and (7) respectively were selected for the kinetics of methane production.
( ) where, M is cumulative methane production (m 3 /kg VS), λ is lag-phase time (days), P is methane production, potential (m 3 /kg VS), R is methane production rate (m 3 /kg VS day), t is digestion time (days), k is hydrolysis rate constant (day −1 ) and n is a dimensionless shape factor.

Data Analysis
Non-linear regression analysis was used in analysis of the models. Solver tool for optimisation in Excel was used to fit the kinetic models already presented in Equations (5), and (7) onto the experimental data to predict the lag-phase time, hydrolysis rate constant, rate of methane production and methane potential. This was achieved by minimizing the sum of square deviations between the measured cumulative methane and the simulated methane values represented in Equation (8). To determine how well the models fitted the data, coefficients of determination were determined through regression analysis using data analysis tools in Excel. Anova single factor analysis was performed to determine whether there was significant difference between the different substrate ratios in terms of methane content in biogas and methane yield. Where a significant difference was found at 5% significant level, a post-hoc t-test of two sample assuming equal variances was used to differentiate the substrate ratios that were significantly different from one another.
= min (8) where, f is the objective function, ε represents other parameters the models depend upon, V m is the simulated cumulative methane values.

Substrate Composition
Characteristics of the substrates that were used in the batch experiment are presented in Table 1 and Table 2. The composition of K in MP (1:0:0) was extremely high which agrees with data in literature [26]. However, K composition in CP (0:1:0) and SP (0:0:1) was less than 2500 mg/kg which is in the acceptable range for anaerobic digestion [27]. Conversely, Co and Mn were not detected in MP while in CP and SP were above the recommended range. On a similar note, Ca, Na, Mg and Cu were not detected in CP while in MP and SP they were either low, high or within the acceptable range. Such imbalance of nutrients in standalone substrates is said to be rectified when they are mixed together. For instance Mg was not detected in CP (0:1:0) but when mixed with SP (0:0:1) whose Mg composition was 52.89 mg/kg, the composition of the mixtures ranged from 9.88 to 13.22 mg/kg. On the other hand, Na was not detected in CP and SP and when they were mixed together in different proportions, still it was not detected suggesting that selection of substrates to include in the mixture is very important.
These results show that it would be good to use sweet potato peels (SP*) alone without mixing them with SV due to their low %cellulose of 5.94 and low lignin (%ADL) of 2.69% compared to SV with 27.88% of cellulose and 5.3% of lignin. However, SP* shows to be low in protein content (3.13%) due to its low %N (0.5) leading to high C:N ratio of 81.15 compared to SV. Therefore, mixing SV and SP* to form SP, represented by the ratio of 0:0:1, improved the C:N ratio to 45.67 which is near the 20 -30 range said to be good for anaerobic digestion [28]. Though, [6] says a ratio of C:N ranging from 20 to 50 is acceptable depending on the substrate. Proximate composition results for MP (1:0:0) obtained in this study are comparable to those of Anhwange, Ugye [29] and [30], although they were looking at proximate composition of banana peels from other banana species and varieties. Similarly, proximate composition results for CP (0:1:0) were comparable to those of Oboh [31], in particular the crude protein and carbohydrates, while crude fibre and fat content compare to results determined by Ofoefule and Uzodinma [32]. It is important to note that even when parameters said to influence proximate composition like differences in cultivars, time of harvest or maturity and location among others, were not considered in the study, the findings are still comparable to other varieties from different areas in different countries. Generally mixing MP, CP and SP showed to balance nutrient composition in the mixtures favourable for anaerobic digestion.

Effect of MP:CP:SP Ratio Variation on Methane Yield
Results of batch digestion tests for the varying proportions of MP, CP and SP are presented in Table 3. Matooke peels (MP) represented by ratio 1:0:0 produced 0.28  [16]. On the other hand, SP produced 0.14 Nm 3 CH 4 /kg VS. Given that SP was a mixture of SP* and SV, data in literature to compare it with was not found. On a similar note, the rest of the ratios were not compared with any literature as it could not be found and are therefore discussed with respect to their parent substrates.
Varying MP and CP, showed a significant difference (P < 0.00) in terms of methane yield between the substrate ratios. It is noted that 1:1:0 and 1:2:0 ratios were no better than their parent substrates (0:1:0 and 1:0:0). This suggests that mixing two or more substrates in co-digestion is not a guarantee to achieve higher methane yields than the original/parent materials in the mix. There has to be an optimal ratio. Therefore ratio 2:1:0 that produced 0.32 Nm 3 /kg VS of methane, though not significant from 0.28 Nm 3 /kg VS of methane produced from 1:0:0 ratio at 5% significant level, is said to have improved methane production of MP by 15% and that of CP by 79% (Figure 1(a)).
Unlike MP and CP, variation of MP and SP showed that SP had a synergetic effect on MP at 1:0:1 and 2:0:1 ratios. Ratio 1:0:1 produced 0.38 Nm 3 /kg VS of methane signifying an increase of 36% for MP and 177% for SP while, ratio 2:0:1 produced 0.41 Nm 3 /kg VS implying an increase in methane yield of 47% for MP and 199% for SP (Figure 1(b)).  Again it can be seen that ratio 1:0:2 produced 0.21 Nm 3 /kg VS of methane which is 23% less compared to 0.28 Nm 3 /kg VS from MP (1:0:0). Since SP was a combination of SV and SP*, Figure 1(b) shows that methane yield increased by over 200% for SV and 29% for SP* at 2:0:1 ratio and 40% for SV and 19% for SP* at 1:0:1 ratio. On the other hand, varying CP and SP had a synergetic effect on each other. That is to say, all the ratios 0:1:1, 0:1:2 and 0:2:1 produced methane that was significant (P < 0.00) compared to their parent substrates. This can generally be seen clearly in Figure 1(c) which shows great increase of over 100% in methane yield for the three ratios referenced to their parent substrates. Similarly, Figure 1( It is important to note that no buffer for pH control like the common Na 2 CO 3 was used in this batch co-digestion experiment. This was intended to assess whether the three substrates in question would have a buffering effect towards one another for pH control. Although it is generally seen in Table 3 that pH had decreased from pH start to pH end by the end of the experiment, it is worthy to note that pH end values for a number of ratios were still within the pH range of 6.1 -8.3 said to favour anaerobic digestion [20]. This is not to mention that with the exception of ratio 0:1:2, pH end values for the other ratios where CP and SP were varied as well as ratio 1:4:1 were below the acceptable pH range. In addition, ratios where CP or SP was more than performed poorly in terms of methane yield except ratio 0:2:1, though its pH end value was 5.9. This is out of the acceptable 6.1 -8.3 pH range. This poor performance of the ratios having more CP than MP or SP can be partly explained by the kinetics of methane production presented in section 3.3. However, Panichnumsin, Nopharatana [38], also noted a decrease in methane production when CP exceeded 60% in co-digestion of CP and pig manure. Unlike methane yield that showed significant difference at 5% significant level be-  Figure 2, Figure 3 and Figure 4 show MGM, exponential and Cone models respectively fitted onto experimental methane production data for the respective MP:CP:SP ratios. Table 4 presents the kinetic parameters estimated by the respective models. To guide the discussion, any model whose coefficient of determination (R 2 ) was 0.99 and above (explaining 99% and above of the variations in experimental data) was considered to have adequately fitted the data. Otherwise it was taken not to have adequately fitted the data given that hydrolysis rate of any substrate can be represented by its methane production only when there is no accumulation of intermediary compounds,   according to Veeken and Hamelers [24]. Therefore, any MP:CP:SP ratio whose variations in methane production data were not adequately explained by the models (R 2 < 0.99) was assumed to have had inhibition from intermediary compounds.

Estimation of Kinetics of Methane Production Parameters
Considering the parent substrate MP, CP, SP, SP* and SV) alone, it can be generally said that MGM and Cone model adequately fitted methane production from MP and SP with 99% of the variations explained by the models. However, considering their hydrolysis rate constants (k), MP, CP and SP* had relatively high k-values of 0.547, 0.96 and 0.912 respectively compared to SV and SP. For SV having a low k-value of 0.254 could be attributed to its high content of cellulose compared to the other parent substrates (Table 2 Section 3.1). On the other hand, the relatively low k-value for SP of 0.394 could be due to the synergetic effect resulting from mixing SV and SP*. The high k-values for MP, CP and SP* affirm to very fast hydrolysis processes during digestion which may lead to accumulation of intermediary products beyond which the methanogenesis process can handle, thus resulting into process instability. This is in agreement with other researchers [6] [7]. It is therefore hoped that using them in mixed mode would reduce on their respective hydrolysis rates thus stabilizing the digestion process.
When MP and CP were varied, only methane production from ratio 2:1:0 was adequately fitted by the MGM and Cone model. The hydrolysis rate constant k for 2:1:0 ratio estimated by the Cone model was reduced to 0.277 d −1 which is in the range of 0.005 NB: 0:0:1 ratio is a mixture of sweet potato peels (SP*) and vines (SV), M is cumulative methane production (m 3 /kg VS), λ is lag-phase time (days), P is methane production potential (m 3 /kg VS), R is methane production rate (m 3 /kg VS day), t is digestion time (days), k is hydrolysis rate constant (day −1 ) and n is a dimensionless shape factor. effect on hydrolysis rate when two parts of MP are mixed with one part of CP for anaerobic digestion. Similarly, there was a synergetic effect on lag phase (λ) for ratio 2:1:0 as it was reduced to 0.175 days compared to 0.731 and 0.188 days for MP and CP respectively. Although ratios 1:2:0 and 1:1:0 had also k and λ values synergised, none of the models fitted them, an indication of inhibition evidenced by their pH end in Table 3.
However, when MP and SP were varied, all the three ratios (   For this study, HRT is taken as the time when over 90% of the cumulative specific methane yield is attained for any given MP:CP:SP ratio [42]. It is taken as the time when cumulative specific methane production curve for any substrate levels off or reaches its plateau according to Zhang, El-Mashad [43]. Therefore, from Figures 5(a)-(c) it can be Figure 5. Illustration of temperature effect on digestion process and substrate ratio effect on HRT. NB: Lines 1 and 2 represent HRT for the respective substrate ratios when cumulative specific methane production curves level off. The dotted line T = 26 deg shows temperature shock effect on methane production process.

Effect of MP:CP:SP Ratio Variation on Hydraulic Retention Time
seen that the parent substrates (MP, CP, SP and SP*) with the exception of SV had reached over 90% of their methane yield by day five as indicated by vertical line 1. This short HRT of less than 5 days is attributed to their fast digestion evidenced by their high hydrolysis rate constants presented in Table 4. Using such substrates with such a short HRT are susceptible to washout of microorganisms which may result into volatile fatty acids build up leading to acidification [8] [26]. Similarly, this short HRT of less than 5 days was exhibited by ratios 1:0:1 and 1:0:2 in Figure 5(b) and ratio 4:1:1 in Figure 5(d). This again is due to the quick digestion of these ratios as supported by their high hydrolysis rate constants of 0.684, 0.779 and 0.798 d −1 respectively.
Nevertheless, co-digestion has shown that HRT for ratio 2:1:0 went up to more than 15 days when MP and CP were varied. When MP and SP were varied, HRT for ratio 2:0:1 also increased to more than 15 days as indicated by vertical line 2 compared to its parent substrates. The same is true for ratios 0:1:1, 0:1:2 and 0:2:1 when CP and SP were varied and ratios 1:1:1 and 1:1:4 when MP, CP and SP varied. Results reported by Heo, Park [12] also show that when food waste with HRT of less than 5 days was co-digested with waste activated sludge at a ratio of 1:1, HRT reached 13 days for optimum digester performance. Figure 6 shows temperature and pressure profiles recorded during the batch digestion experiment of MP:CP:SP ratio variations. Looking at both Figure 5 and Figure 6, Figure 5 shows that between the first day and third day, the cumulative methane curves were generally running smooth when temperature fluctuations in Figure 6 were going below 30˚C but not reaching to 26˚C between the same period. Smooth cumulative methane curves are an indication that the rate of biogas production was constant in the stated period even when temperature reached to 27˚C. However, when temperature fluctuations in the incubator reached to 26˚C after the third day, shown by the "star marker" in Figure 6, there was a setback in the digestion process. This setback is evidenced by the kinks in the cumulative methane curves in Figure 5 at points shown by the dotted lines labeled T = 26 deg which is as a result of a reduction in biogas production rate. Nevertheless, when the temperature rose again, the curves smoothened for most of the substrate ratios with the exception of ratios 1:1:0 and 1:2:0 in Figure 5(a),

Temperature Effect on Anaerobic Co-Digestion of MP, CP and SP Process Stability
implying that microorganisms were able to adjust after temperature shocks. This adjusting of microorganisms after temperature shock was also noted by Chae, Jang [44].
Datta [45], also found that methane production was inhibited by temperature ≤ 25˚C.
In addition, the kinks in curves were more pronounced in substrate ratios consisting of CP ( Figure 5(a), Figure 5(c), Figure 5(d)) than those where MP and SP were varying ( Figure 5(b)). This suggests that temperature could be substrate specific and in this case to CP, given that substrate ratios consisting of cassava peels were more susceptible to instability when temperature dropped to 26˚C. Therefore, there is need to maintain Figure 6. Incubator temperature and ambient pressure changes. temperature above 26˚C for effective anaerobic digestion of substrates consisting MP, CP and SP, but more particularly CP, though more investigations are needed. Unlike methane yield, kinetic parameters and HRT that showed to be dependent on the proportions of the substrates in the mixtures, temperature effect was not, but specific to a particular substrate in the mixture. Generally, ratios 2:1:0, 2:0:1, 1:0:1, 0:1:1, 0:1:2, 0:2:1 1:1:1 and 1:1:4 were able to with stand temperature shocks, which are almost the same ratios that had a synergetic effect on HRT, kinetics of methane production, pH and methane yield.

Conclusions
1) Results from this study have shown that MP, CP, SP, and SP* are digested very fast due to their high rate of hydrolysis compared to SV whose hydrolysis rate is low due to its high content of cellulose. As a result of this fast digestion, they have a short HRT of less than 5 days which may result into washout of microorganisms leading to accumulation of intermediary products that inhibit methanogenesis process.
2) Co-digestion of MP, CP and SP has shown to have synergetic effect on HRT, methane yield and its kinetics (hydrolysis rate constant and lag phase). In addition, HRT, methane yield and its kinetics have shown to be dependent on the proportion of any given parent substrate in a mixture. lag phase and methane yield, withstanding temperature shock and maintaining a relatively favourable pH for methanogens. Therefore, anyone having at least two of MP, CP and SP has an optimum ratio that favours methane yield and the factors that may affect its production.
6) However, more investigations are needed when using bigger samples as there could have been sampling errors while using small samples for the tested substrates. In addition, ways of how to maintain pH in the optimum range of 6.6 -7.3 without using the usual Na 2 CO 3 as a buffering agent is another area for research to enable use of household waste as a potential standalone feedstock for biogas production.
grant Training the next generation of scientists provided by Carnegie Cooperation of New York through the Regional Universities Forum for Capacity Building in Agriculture (RUFORUM). Peter is a METEGA sponsored graduate student at University of Botswana.