Determination of Methane and Carbon Dioxide Formation Rate Constants for Semi-continuously Fed Anaerobic Digesters

To optimize commercial-scale biogas production, it is important to evaluate the performance of each microbial step in the anaerobic process. Hydrolysis and methanogenesis are usually the rate-limiting steps during digestion of organic waste and by-products. By measuring biogas production and methane concentrations on-line in a semi-continuously fed reactor, gas kinetics can be evaluated. In this study, the rate constants of the fermentative hydrolysis step (kc) and the methanogenesis step (km) were determined and evaluated in a continuously stirred tank laboratory-scale reactor treating food and slaughterhouse waste and glycerin. A process additive containing Fe 2+ , Co 2+ and Ni 2+ was supplied until day 89, after which Ni 2+ was omitted. The omission resulted in a rapid decline in the methanogenesis rate constant (km) to 70% of the level observed when Ni 2+ was present, while kc remained unaffected. This suggests that Ni 2+ mainly affects the methanogenic rather than the hydrolytic microorganisms in the system. However, no effect was initially observed when using conventional process monitoring parameters such as biogas yield and volatile fatty acid concentration. Hence, formation rate constants can 646 reveal additional information on process performance and km can be used as a complement to conventional process monitoring tools for semi-continuously fed anaerobic digesters.


Introduction
Anaerobic digestion is the microbiological degradation of organic materials in the absence of oxygen, a process which can be exploited in commercial plants to produce biogas consisting of carbon dioxide and energy-rich methane.The anaerobic digestion process has been described in detail elsewhere [1,2], but in brief it can be separated into four different steps: (i) Hydrolysis, whereby particulate organic material is converted to mono-or oligomers by hydrolytic enzymes; (ii) acidogenesis, whereby the products from hydrolysis are converted to volatile fatty acids (VFA), alcohols, hydrogen and carbon dioxide by primary fermentative bacteria, the same microorganisms that excrete hydrolytic enzymes; (iii) acetogenesis, whereby the different products of acidogenesis are degraded by secondary (syntrophic) fermentative bacteria to acetate, carbon dioxide and hydrogen; and (iv) methanogenesis, whereby acetate is converted by acetotrophic methanogens to methane and carbon dioxide, or syntrophically oxidized to carbon dioxide and hydrogen by syntrophic acetate-oxidizing bacteria [3].The hydrogen and carbon dioxide produced during acetogenesis and syntrophic acetate oxidization are subsequently converted into methane by hydrogenotrophic methanogens.
Bio-methane is a renewable energy source that can be further used for combined heat and power or for vehicle fuel.In the past, bio-methane was generally produced and collected at wastewater treatment plants and landfills.However, there is now an increasing trend for co-digestion of different energy-rich wastes such as slaughterhouse waste, the organic fraction of municipal solid waste (OFMSW) and industrial waste, with the aim of producing methane for commercial use [4].This trend for more commercial-scale methane production has created a need for optimization of biogas processes, for example by trace element addition.Trace elements are nutrients required for microbial growth, and supplying optimal concentrations of these ensures efficient anaerobic digestion [5].For anaerobic digestion, iron (Fe), cobalt (Co) and nickel (Ni) are considered most important and the effects of adding these nutrients have been evaluated in several studies [6][7][8][9][10][11].However, optimization of the anaerobic process and evaluation of the effects of various optimization efforts, for example addition of trace elements, requires accurate methods to determine process stability (i.e., risk of process disruption) and/or process efficiency and to identify process bottlenecks.
Process stability is commonly determined and bottlenecks identified by measuring the accumulation of intermediate products in the anaerobic process, such as VFA [12] and hydrogen [13].Accumulation of VFA can result in a pH decrease, so both pH and alkalinity can also be used as stability indicators [12].Overall process efficiency is commonly determined by gas yield per unit volatile solids (VS) or chemical oxygen demand (COD), total gas production or degree of degradation (destruction of organic material, on a VS or COD basis) [14].However, most of these analyses involve time-consuming laboratory work for determination of VFA, VS or COD, leading to long response times.
An alternative approach would be to study the gas kinetics, i.e., the manner in which the gas is produced after addition of organic material to the digestion reactor.Batch experiments are generally used to determine the bio-methane potential (BMP) of substrates, but the gas production characteristics (typically decelerated exponential growth shape) can also be used to determine hydrolysis rate constants, provided that hydrolysis is the rate-limiting step [15].Hydrolysis rate constants can be determined by fitting batch experiment data to a first-order equation.In a semi-continuously fed reactor, the accumulated biogas production between feedings also shows a decelerated exponential growth shape [16].The kinetics information provided by online measurement of biogas production can be utilized to obtain additional information about the process, apart from daily gas production.Astals et al. [16,17] found that daily rate constants can vary due to the proportion of easily biodegradable organic matter in the substrate, leading to clear effects on daily gas production.
The main aim of this study was to evaluate the possibility to determine the carbon dioxide formation rate constants and methane formation rate constant between feedings at laboratory scale to observe process stability when co-digesting OFMSW, slaughterhouse waste and glycerin.To achieve this, we used the methane concentration in the digester headspace and the biogas flow to distinguish between methane and carbon dioxide formation rates.A further aim was to use the rate constants determined to identify bottlenecks in anaerobic digestion of the substrate and to determine the step in which addition of nickel had an effect.

Reactor Performance
After an initial period of stable operation, the OLR (Organic loading rate) was increased from 3.4 to 4.7 kg• VS• m −3 •d −1 in order to stress the process and provoke process disturbances (Figure 1A).However, even though the HRT (hydraulic retention time) was as short as 18 days no process instability was observed, probably due to the continuous addition of the trace element solution including iron, cobalt and nickel.When the OLR was increased further (to 6.5 kg• VS• m −3 •d −1 ), there was still no VFA accumulation or signs of lower specific biogas production (Figure 1A).Previous studies have shown that trace element addition has positive effects, especially on process stability, i.e., reduced VFA accumulation [9,10,18].However, the effects of nickel and cobalt could not be separated from these experimental data and therefore nickel administration was terminated on day 89.This resulted in a slight increase in VFA 35 days later (on day 124), which corresponded well with the calculated time for total wash-out of nickel from the digester (Figure 1A).This indicates that nickel administration had positive effects on process stability which at higher OLR could not be achieved with cobalt administration alone.However, while a small fraction of VFA started to accumulate, no significant decrease in specific methane production was observed after nickel addition ceased.

Determination of Kinetic Rate Constants
By on-line measurement of the accumulated biogas production and the methane concentration of the gas every five minutes, the kinetic behavior of the system between daily feedings was analyzed.From the measured formed biogas and the change in headspace methane content, the specific production of methane and residual gas can separately be calculated for each five minute interval.For simplicity, the residual gas was assumed to contain 100% carbon dioxide.
When plotting the accumulated specific biogas, methane and carbon dioxide production over time, the curves displayed a decelerated exponential growth shape between feedings (Figure 2).The specific methane production rate, km (h −1 ), was calculated by fitting data to a first-order kinetic Equation (1) similarly as common practice for BMP-tests [15].m 4spec 4max CH ( ) CH ( 1) On examining the curve of carbon dioxide content between feedings (Figure 2), it was obvious that an initial peak of carbon dioxide release directly after feeding was followed by slower formation.Consequently, data was fitted to two lumped first-order (Equation ( 2)) functions to describe the total process of CO2 formation, generating a slow and fast formation rate constant, kc1 and kc2 (h −1 ), respectively:

Rate Constant of Methane Formation (km)
The specific methane production per day remained constant throughout the experiment (Figure 1A), but the rate constant of daily methane formation (km) changed considerably over time (Figure 1B).After introducing administration of additional nickel at day 0, an increase in km occurred, coinciding with increasing OLR (from day 26).The apparent relationship between OLR and rate of methane formation can only be explained speculatively.It could be expected that increasing the OLR would positively influence the total volume of gas produced, but not the rate of production.However, in order to convert a larger amount of added VS per day at a higher OLR, the rate of conversion of the products from acetogenesis would need to increase to avoid accumulation of acetate.Thus, it is possible that the difference in methane formation rate between OLR 4.7 and 6.5 kg• VS• m −3 •d −1 was a consequence of higher acetogenesis product formation, allowing for a higher methane formation rate, rather than a direct effect of nickel.This is further supported by the fact that the increase in OLR was only as additional glycerin, a trivalent carbon molecule directly degraded during acetogenesis.In this case, methanogenesis was not the rate-limiting step when nickel was present.There was an increase in km (0.060 ± 0.012 h −1 before nickel administration) reaching a plateau at about 0.12 h −1 (Figure 1B).This plateau level could be the result of the OLR being stable (from day 40), the available nickel being limited or a combination of both (Figure 1B).Nevertheless, when nickel administration was terminated at day 89 and the reactor concentration of the trace metal decreased, the plateau km level was only maintained for a short period until the availability of nickel obviously became the limiting factor (day 100).Thereafter, km rapidly decreased with the nickel concentration until a final plateau (0.037 ± 0.002 h −1 ) was reached (day 132) independently of OLR, which was stable during this period (Figure 1B).At approximately the same time as km reached its lower plateau, VFA started to accumulate, indicating a connection between these two processes.The accumulation of VFA shows that the km at the lower plateau was too low and that methanogenesis became rate-limiting for the process.These data strongly suggest that nickel is an essential metal catalyst for the formation of methane during anaerobic digestion and that without nickel, methanogenesis became the rate-limiting step in the anaerobic digestion process studied here, which was kept at a high OLR and short HRT.These results are in line with previous findings highlighting the effect of nickel, especially on methanogenesis [11,[19][20][21].They also confirm findings that nickel is an important metal in key enzymes for methanogenesis, such as methyl-coenzyme-M cofactor F430, carbon monoxide dehydrogenase and other hydrogenases [22,23].
In summary, it should be noted that the km value changed upon nickel administration/OLR increase and nickel removal long before commonly used process stability indicators such as VFA or specific methane production responded.The km value could therefore potentially be used as an early warning indicator of methanogenesis instability in semi-continuously fed anaerobic digestion processes.

Rate Constants of Carbon Dioxide Formation (kc1 and kc2)
An initial peak in carbon dioxide release directly after feeding was followed by a slower release (Figure 2).The fast, initial release of carbon dioxide is probably partly attributable to a pH decrease due to VFA in the feed [24] and partly to fast hydrolysis and fermentation of dissolved organic matter during acidogenesis and acetogenesis.The slower release of carbon dioxide could be attributable to the slower continuous processes of hydrolysis of less readily degradable organic matter and methanogenesis.
In contrast to km, after introducing nickel at day 0 there was a slight decrease in the slow carbon dioxide formation rate (kc1).kc1 decreased from 0.05 h −1 (day -25 to 0) down to 0.02 h −1 day 76 when the OLR was increased in presence of nickel (Figure 1C).After termination of nickel administration, kc1 rapidly increased from 0.04 h −1 (day 100) to 0.09 h −1 (day 109).This increase in kc1 occurred simultaneously with the start of the decrease in km.In fact, there was a significant negative correlation (p < 0.05) between kc1 and km from the start of OLR increase (day 26) until day 109 (Figure 1C).However, after day 109 the rate constants were positively correlated.Thus, it is likely that according to our hypothesis regarding carbon dioxide emanating from slow hydrolysis and methanogenesis, kc1 is affected by km.One plausible explanation would be that nickel administration stimulates the degradation of lipids, contributing to a higher methane content in the gas produced, and thus lower carbon dioxide formation rate, but higher km.Karlsson et al. [7] showed that nickel administration was clearly correlated with increased process performance and proliferation of Methanosarcinales.Lipid oxidation is highly dependent on low hydrogen pressure [13], and since Methanosarcinaceae can act as a hydrogen-consuming methanogen, lowering the hydrogen pressure [25], this would result in more favourable lipid degradation conditions.One other explanation could be that nickel stimulates hydrogenotrophic methanogens, thus resulting in a high rate of methanogenesis and consumption of carbon dioxide.Consequently, km is affected negatively and kc1 positively when nickel is removed.Assuming that kc1 is dependent on the hydrolysis of less readily degradable organic matter, it might be comparable to the lumped apparent first-order coefficient (usually denoted khyd) obtained from BMP testing.Vavilin, et al. [26] and Li, et al. [27] have determined a number of kinetic coefficients in mesophilic conditions, among which kitchen waste (0.34 and 0.18 d −1 ), food waste (0.55 d −1 ) and slaughterhouse waste (0.35 d −1 ) are relevant to compare with the substrate mix in the current paper (kc1 0.97 ± 0.49 d −1 ).Furthermore, Elbeshbishy & Nakhla [28] have reported hydrolysis constants for particulate protein (kprot 0.65 d −1 ) and carbohydrate (khyd 0.78 d −1 ).These rates are comparable to our results and it is thus possible that the khyd observed in BMP tests is represented by kc1.The reason why kc1 is slightly higher than khyd could be an effect of the (semi)-continuous operation in our study, which may lead to higher microbial activity because of the already active microorganisms and the higher inoculum: Substrate ratio.In addition, Batstone et al. [29] attempted to obtain khyd for full-scale continuous systems but the coefficient varied substantially, by several orders of magnitude, and identifiability of the constant was generally poor.The estimated rate constants, obtained from the gas flow in the full-scale system, were an order of magnitude higher (>5 d −1 ) than the constants obtained from BMP testing of the substrates.
In contrast to km and kc1, after introducing nickel there was no significant effect on the fast carbon dioxide formation constant (kc2) of increasing OLR (Figure 1C).Interestingly, kc2 also remained constant after termination of nickel administration (apart from a slight increase between days 101 and 109).There was no significant correlation (p > 0.05) between kc2 and km from the start of OLR increase (day 26) until day 109 (Figure 1C).Hence, in contrast to methanogenesis, the hydrolysis-fermentation process was not affected by the presence of nickel in the reactor.

Reactor Set-Up
One semi-continuous laboratory-scale reactor [30], with an active volume of approximately 9 L, running at 38 °C was used.The reactor was fed semi-continuously at 24-h intervals.The organic loading rate (OLR) varied according to the schedule shown in Table 1.The substrate used was OFMSW (59%-82% of OLR), slaughterhouse waste (13%-18% of OLR) and glycerin (0%-28% of OLR).The OFMSW slurry (TS, Total solids: 14.9% and VS: 13.3%) was collected from the full-scale pre-treatment plant of the city of Linköping, Sweden were the organic fraction of municipal solid waste (source-separated) was pre-treated to form substrate slurry prior to digestion.The slaughterhouse waste (TS: 11.6% and VS: 10.6%) originated from the local slaughterhouse in Linköping (Scan, Sweden) which mainly slaughter cattle and pigs.Both the OFMSW and slaughterhouse waste fraction were pasteurized in lab-scale (1 h, 70 °C); the substrates was then stored at −18 °C.The glycerine (TS: 75.0% and VS: 67.0%) was a residual product from biodiesel production; it was stored in room temperature during the experiment.The hydraulic retention time (HRT) was set to decrease from 21 to 18 days during the start-up period (day −25 to 0) (Table 1).Tap water was used to dilute the substrate mixture, enabling a constant HRT when the OLR was increased further during the experiment (by lowering the fraction of tap water).On day 54, 85 and 116, the OLR was unintentionally 4.7, 6.0 and 4.7 kg• VS• m −3 •d −1 , respectively.Table 1.Schedule of changes to organic loading rate (slaughterhouse waste, the organic fraction of municipal solid waste (OFMSW) and glycerin) and hydraulic retention time (HRT) in the semi-continuously fed anaerobic digester during the experimental period.A process additive with the goal to achieve a steady-state concentration in the reactor of 0.5 mg• L −1 Co 2+ (as CoCl2• 6H2O) and 300 mg• L −1 Fe 2+ (as FeCl2• 4H2O) dissolved in water and HCl was used for the whole experiment (day −25 to day 168) [18].The first 25-day period (day -25 to day 0) was used as a reference period with no nickel addition.From day 1 to 89, a Ni 2+ supplement (as NiCl2• 6H2O) was added to the substrate mixture.The dosage of nickel was done to reach a steady-state concentration of 0.4-0.5 mg• L −1 additional nickel in the digestate within seven days, and then keep a constant additional nickel concentration of 0.4-0.5 mg• L −1 .The substrate mixture contained a background nickel concentration of approximately 0.4 mg• kg −1 .

Experimental Monitoring and Analyses
A MGC-10 milligas counter (Ritter, Bochum, Germany) was used to measure volumetric biogas production and methane gas concentration was determined every five minutes with a gas sensor (Bluesense, Herten, Germany).VFA were analyzed using a gas chromatograph (Clarus 550, PerkinElmer, Waltham, MA, USA) equipped with a 30 m × 0.32 mm × 0.25 μm Elite-FFAP column [31].Total solids (TS) content was measured by oven-drying at 105 °C for 20 h.Content of VS was subsequently measured gravimetrically by combusting the TS sample for 3 h at 550 °C.The TS content of the substrate was compensated for the loss of VFA in order to perform calculations on OLR and specific biogas production according to Vahlberg et al. [32].All volumetric gas data presented in this paper have been converted to standard conditions at pressure 1.01325 bars and temperature 273.2 K.The daily OLR was calculated by dividing the VS weight fed into the digester by the reactor volume plus feeding volume of the previous day.

Conclusions
This study evaluated the effects of trace elements, particularly nickel, on process stability in a semi-continuously fed anaerobic digester.Use of conventional process parameters such as VFA and daily methane production failed to identify the positive effects of nickel supplementation of the process.However, an alternative method analyzing the process kinetics during one feeding cycle allowed the changes in methane and carbon dioxide formation rates to be determined over the entire experimental period.Using this method, we were able to show that nickel administration had only indirect or no effects on the hydrolysis rate constant (kc), but a great impact on the methanogenesis rate constant (km).At the high organic loading rates applied in the study, km was 98%-220% higher in the presence of nickel than when nickel was omitted, whereas general process monitoring parameters such as gas production and VFA showed only small indications of differences in process stability.This indicates that km can be used as an efficient early warning indicator of instability in semi-continuously fed anaerobic digesters.