Highlights
- •The IPCC 2019 Tier 2 equations predict methane emission adequately.
- •The knowledge of DE improves the accuracy of CH4 prediction.
- •Hay-based diets reduce digestibility and consequently a measured DE should be applied.
- •Local values of the CH4 conversion factor (Ym) and DE are proposed for the Mediterranean area.
Abstract
The study aimed to evaluate Intergovernmental Panel on Climate Change (IPCC) Tier 2 (2006 and 2019) to predict enteric CH4 emissions from lactating cows fed Mediterranean diets. The effects of the CH4 conversion factor (Ym; CH4 energy loss as a percentage of gross energy intake) and digestible energy (DE) of the diet were evaluated as model predictors. A data set was created using individual observations derived from 3 in vivo studies on lactating dairy cows housed in respiration chambers and fed diets typical of the Mediterranean region based on silages and hays. Five models using different Ym and DE were evaluated following a Tier 2 approach: (1) average values of Ym (6.5%) and DE (70%) from IPCC (2006); (2) average value of Ym (5.7%) and DE (70.0%) from IPCC (2019; 1YM); (3) Ym = 5.7% and DE measured in vivo (1YMIV); (4) Ym = 5.7 or 6.0%, depending on dietary NDF, and DE = 70% (2YM); and (5) Ym = 5.7 or 6.0%, depending on dietary NDF, and DE measured in vivo (2YMIV). Finally, a Tier 2 model for Mediterranean diets (MED) was derived from the Italian data set (Ym = 5.58%; DE = 69.9% for silage-based diets and 64.8% for hay-based diets) and validated on an independent data set of cows fed Mediterranean diets. The most accurate models tested were 2YMIV, 2YM, and 1YMIV with predictions of 384, 377, and 377 (g of CH4/d), respectively, versus the in vivo value of 381. The most precise model was 1YM (slope bias = 1.88%; r = 0.63). Overall, 1YM showed the highest concordance correlation coefficient value (0.579), followed by 1YMIV (0.569). Cross-validation on an independent data set of cows fed Mediterranean diets (corn silage and alfalfa hay) resulted in concordance correlation coefficient of 0.492 and 0.485 for 1YM and MED, respectively. The prediction of MED (397) was more accurate than 1YM (405) in comparison with the corresponding in vivo value of 396 g of CH4/d. The results of this study showed that the average values proposed by IPCC (2019) can adequately predict CH4 emissions from cows fed typical Mediterranean diets. However, the use of specific factors for the Mediterranean area, such as DE, improved the accuracy of the models.
Graphical Abstract

Graphical AbstractSummary: Intergovernmental Panel on Climate Change (IPCC) equations were evaluated to study the precision and accuracy in predicting methane (CH4) emission by lactating dairy cows. The predicted values were compared with in vivo values obtained from experimental studies conducted using cows fed diets based on hay or silage forages in the Mediterranean region. Hay-based diets differed significantly from silage-based diets having lower energy digestibility (DE). The IPCC 2019 equations predicted CH4 emission well and hence may be used for the computation of greenhouse gas inventories. However, specific equations are proposed to predict CH4 emission by dairy cows in the Mediterranean region.
Livestock and manure management contribute to 5.8% of anthropogenic-caused greenhouse gas (GHG) emissions, and considering CH4 emissions only, enteric fermentation and manure management represent 32% of anthropic CH4 emissions (
World Resources Institute, 2019
). Therefore, it is important to predict the amount of enteric CH4 emitted in a specific livestock production system (- World Resources Institute
World Greenhouse Gas Emissions: 2019.
https://www.wri.org/data/world-greenhouse-gas-emissions-2019
Date: 2019
Date accessed: November 3, 2021
Niu et al., 2018
). The Intergovernmental Panel on Climate Change (IPCC) guidelines for GHG inventories were developed to provide internationally agreed methodologies for estimating GHG emissions (- Niu M.
- Kebreab E.
- Hristov A.N.
- Oh J.
- Arndt C.
- Bannink A.
- Bayat A.R.
- Brito A.F.
- Boland T.
- Casper D.
- Crompton L.A.
- Dijkstra J.
- Eugène M.A.
- Garnsworthy P.C.
- Haque M.N.
- Hellwing A.L.F.
- Huhtanen P.
- Kreuzer M.
- Kuhla B.
- Lund P.
- Madsen J.
- Martin C.
- McClelland S.C.
- McGee M.
- Moate P.J.
- Muetzel S.
- Muñoz C.
- O'Kiely P.
- Peiren N.
- Reynolds C.K.
- Schwarm A.
- Shingfield K.J.
- Storlien T.M.
- Weisbjerg M.R.
- Yáñez-Ruiz D.R.
- Yu Z.
Prediction of enteric methane production, yield, and intensity in dairy cattle using an intercontinental database.
IPCC, 2006a
). The latest version of IPCC, 2019b
was published to refine the previous version (IPCC, 2006b
), providing supplementary methodologies and updated default values. Regarding enteric CH4 emission from ruminants, the conversion factor of gross energy intake into enteric CH4 energy (Ym) can be chosen according to the level of productivity and diet characteristics (IPCC, 2019b
). The criteria to choose the value of Ym are milk production, dietary NDF concentration, and digestible energy (DE, as % of GE intake). For example, considering the whole data set, IPCC, 2019a
suggests a value of Ym equal to 5.7% for high-yielding cows (>8,500 kg of milk/head/yr−1), associated with DE ≥70% and NDF ≤35% of diet DM. However, for this production level, the proposed values of NDF and Ym may not represent the diets used in regions with Mediterranean climate characteristics. The Ym values in - IPCC
2019 Refinement to the 2006 IPCC Guidelines for National Greenhouse Gas Inventories. Chapter 10: Emissions from livestock and manure management.
in: IPCC Guidelines for National Greenhouse Gas Inventories Volume 4: Agriculture, forestry and other land use. Intergovernmental Panel on Climate Change,
2019
IPCC, 2019a
were based on the data set of - IPCC
2019 Refinement to the 2006 IPCC Guidelines for National Greenhouse Gas Inventories. Chapter 10: Emissions from livestock and manure management.
in: IPCC Guidelines for National Greenhouse Gas Inventories Volume 4: Agriculture, forestry and other land use. Intergovernmental Panel on Climate Change,
2019
Niu et al., 2018
, which included 154 studies, 82 of them conducted in European research institutes, but none conducted in southern European countries. Diets from northern Europe (based on ryegrass and corn silages, - Niu M.
- Kebreab E.
- Hristov A.N.
- Oh J.
- Arndt C.
- Bannink A.
- Bayat A.R.
- Brito A.F.
- Boland T.
- Casper D.
- Crompton L.A.
- Dijkstra J.
- Eugène M.A.
- Garnsworthy P.C.
- Haque M.N.
- Hellwing A.L.F.
- Huhtanen P.
- Kreuzer M.
- Kuhla B.
- Lund P.
- Madsen J.
- Martin C.
- McClelland S.C.
- McGee M.
- Moate P.J.
- Muetzel S.
- Muñoz C.
- O'Kiely P.
- Peiren N.
- Reynolds C.K.
- Schwarm A.
- Shingfield K.J.
- Storlien T.M.
- Weisbjerg M.R.
- Yáñez-Ruiz D.R.
- Yu Z.
Prediction of enteric methane production, yield, and intensity in dairy cattle using an intercontinental database.
March et al., 2014
) are different from the diets of the Mediterranean region, since the latter are widely based on corn silage and alfalfa/grass hays (Gislon et al., 2020a
). For this reason, IPCC, 2019a
encourages the development of country-specific Ym factors for a more accurate estimation of enteric CH4. The present study aimed to (1) compare the CH4 emissions estimated with - IPCC
2019 Refinement to the 2006 IPCC Guidelines for National Greenhouse Gas Inventories. Chapter 10: Emissions from livestock and manure management.
in: IPCC Guidelines for National Greenhouse Gas Inventories Volume 4: Agriculture, forestry and other land use. Intergovernmental Panel on Climate Change,
2019
IPCC, 2019a
and - IPCC
2019 Refinement to the 2006 IPCC Guidelines for National Greenhouse Gas Inventories. Chapter 10: Emissions from livestock and manure management.
in: IPCC Guidelines for National Greenhouse Gas Inventories Volume 4: Agriculture, forestry and other land use. Intergovernmental Panel on Climate Change,
2019
IPCC, 2006b
equations with the values obtained in vivo by cows fed typical Mediterranean diets; (2) develop possible improvements for the IPCC parameters (Ym and DE) for a more accurate prediction of CH4 emission from lactating cows fed Mediterranean diets.A data set of 66 individual observations was created to evaluate different IPCC equations to predict CH4 emissions from Italian lactating Holstein dairy cows. The data set included individual observations from 3 in vivo studies (
Gislon et al., 2020b
; Pirondini et al., 2015
; and Colombini et al., 2015
), with 32, 16, and 18 observations, respectively. The experiments were carried out at the University of Milan Research Center, in the Po plain (Milan, Italy; 45°30′N, 9°1′E) where the climate of the region is mainly humid subtropical according to the Geiger–Köppen climate classification. The diets fed in each experiment were representative of typical diets fed in the region and were based on the following forages: alfalfa and grass silages, alfalfa and grass hays, wheat silage, corn silage, sorghum forage silage, and sorghum grain whole plant silage. Methane emissions were measured through individual open-circuit respiration chambers, and the DE was determined in vivo by total feces collection. Further details are reported in the above-cited studies, but briefly the air flow through the chambers was measured using diaphragm flow-meters at ambient temperature, pressure, and relative humidity. The air flux was then referred to standard conditions (0°C, 1 atm, relative humidity = 0) by a calculation based on temperature, pressure, and relative humidity measured in the outgoing air entering the flow-meter.Before each respiration chamber cycle, calibrations were made with a certified reference gas with methane concentrations close to the maximum values detectable by the analyzer (2,000 ppm of CH4) and pure N2 to calibrate the minimum value. The calculation of CH4 energy (kJ) was determined by multiplying CH4 production in liters at standard conditions by 39.5388. The production of CH4 in grams was obtained multiplying the CH4 production (in liters at standard conditions) by 0.71682. Five predicting equations were applied based on different Ym and DE values. The first equation (IPCC06) was based on IPCC 2006 (
IPCC, 2006b
) Tier 2 (IPCC Table 10.12) with Ym = 6.5% and DE = 70%. Two other equations were based on IPCC 2019 (IPCC, 2019a
) Tier 2 as follows: (1) the average Ym (5.7%) and DE fixed at 70% were used to calculate the ratios of net energy available in a diet for maintenance and growth to DE consumed (equations 10.14 and 10.15, respectively) and to calculate the gross energy (GE) requirements (equation 10.16; 1YM) or (2) the average Ym (5.7%) and the replacement of fixed DE with the average in vivo digestibility from the experiments (1YMIV). The final 2 equations were based on IPCC 2019 (- IPCC
2019 Refinement to the 2006 IPCC Guidelines for National Greenhouse Gas Inventories. Chapter 10: Emissions from livestock and manure management.
in: IPCC Guidelines for National Greenhouse Gas Inventories Volume 4: Agriculture, forestry and other land use. Intergovernmental Panel on Climate Change,
2019
IPCC, 2019a
; Table 10.12) as follows: for diets with NDF ≤35%, the value 5.7% of Ym was used, whereas for diets with NDF >35%, the value 6.0% of Ym was applied; for DE a fixed digestibility at 70 (2YM) or the average in vivo digestibility from the experiments was used (2YMIV). The emission factor was calculated from IPCC 2019 (equation 10.21) on the basis of the estimate of GE intake since according to - IPCC
2019 Refinement to the 2006 IPCC Guidelines for National Greenhouse Gas Inventories. Chapter 10: Emissions from livestock and manure management.
in: IPCC Guidelines for National Greenhouse Gas Inventories Volume 4: Agriculture, forestry and other land use. Intergovernmental Panel on Climate Change,
2019
IPCC, 2019a
“the Tier 2 estimate of gross energy is the preferred method” for the calculation of the emission factor.- IPCC
2019 Refinement to the 2006 IPCC Guidelines for National Greenhouse Gas Inventories. Chapter 10: Emissions from livestock and manure management.
in: IPCC Guidelines for National Greenhouse Gas Inventories Volume 4: Agriculture, forestry and other land use. Intergovernmental Panel on Climate Change,
2019
An IPCC Tier 2 model for cows fed typical Mediterranean diets (MED) was also proposed using the average value of Ym (5.58) and the average DE values of 69.9 or 64.8% for silage- or hay-based diets respectively obtained from the in vivo experiments, as described in a following section of the paper. The MED was validated on an independent data set based on individual cow observations from the study of
Enriquez-Hidalgo et al., 2020
conducted in a central region of Chile characterized by a Mediterranean climate with hot summers. The Chilean study evaluated 2 diets, one similar to diets typically fed in Italy with a forage basis including corn silage and alfalfa hay, and the other not typically Mediterranean, based on corn silage and a mixture of fresh annual ryegrass and berseem clover (MIX). The CH4 emission was measured with the SF6 technique. For this data set, the 1YM and MED equations were applied. It was not possible to apply the other equations because DE was not measured and the 2 diets had NDF <35% DM.- Enriquez-Hidalgo D.
- Teixeira D.L.
- Pinheiro Machado Filho L.C.
- Hennessy D.
- Toro-Mujica P.
- Williams S.R.O.
- Pereira F.C.
Incorporating a fresh mixed annual ryegrass and berseem clover forage into the winter diet of dairy cows resulted in reduced milk yield, but reduced nitrogen excretion and reduced methane yield.
Front. Vet. Sci. 2020; 7 (33330696)576944
The predictive equations for CH4 emission were evaluated according to
Tedeschi, 2006
. The root mean square prediction error (RMSPE) was decomposed into error due to mean bias (MB), error due to slope bias (SB), and error due to random bias (RB; Bibby and Toutenburg, 1977
).The concordance correlation coefficient (CCC) was calculated according to
where r is the Pearson correlation coefficient and Cb is a bias correction factor (a measure of accuracy), calculated as follows:
where V is a measure of scale shift (i.e., the change in standard deviation between predicted and observed values), and μ is a measure of location shift. V and μ were calculated as below:
where MO and MP are the means of observed and predicted values, respectively, and SDO and SDP are the standard deviations of observed and predicted values, respectively.
Lin, 1989
as follows:
CCC = r × Cb,
where r is the Pearson correlation coefficient and Cb is a bias correction factor (a measure of accuracy), calculated as follows:
Cb = 2/(V + 1/V + µ2),
where V is a measure of scale shift (i.e., the change in standard deviation between predicted and observed values), and μ is a measure of location shift. V and μ were calculated as below:
V = SDO/SDP,
µ = (MO − MP)/(SDO × SDP)1/2,
where MO and MP are the means of observed and predicted values, respectively, and SDO and SDP are the standard deviations of observed and predicted values, respectively.
Data were analyzed using the Proc Mixed procedure of SAS 9.4 (SAS Institute Inc.) with the model Yjkl = μ + Dj + Sk + COWl(Sk) + εjkl, where Yjkl is the dependent variable (DE and Ym); μ is the overall mean; Dj is the diet effect (j = 1, 5); Sk is the random study effect (k = 1, 3); COWl is the random animal effect nested within study (l = 1, 18), and εjkl is the residual error.
The results of IPCC predictions are in Table 1. The RMSPE (%) was similar among models (on average 10.5); all values were <13%, and according to
Kaewpila and Sommart, 2016
, a model may be considered inadequate when RMSPE is higher than 16.0%. The RMSPE reported by Niu et al., 2018
, - Niu M.
- Kebreab E.
- Hristov A.N.
- Oh J.
- Arndt C.
- Bannink A.
- Bayat A.R.
- Brito A.F.
- Boland T.
- Casper D.
- Crompton L.A.
- Dijkstra J.
- Eugène M.A.
- Garnsworthy P.C.
- Haque M.N.
- Hellwing A.L.F.
- Huhtanen P.
- Kreuzer M.
- Kuhla B.
- Lund P.
- Madsen J.
- Martin C.
- McClelland S.C.
- McGee M.
- Moate P.J.
- Muetzel S.
- Muñoz C.
- O'Kiely P.
- Peiren N.
- Reynolds C.K.
- Schwarm A.
- Shingfield K.J.
- Storlien T.M.
- Weisbjerg M.R.
- Yáñez-Ruiz D.R.
- Yu Z.
Prediction of enteric methane production, yield, and intensity in dairy cattle using an intercontinental database.
Niu et al., 2021
) were comparable to those obtained in the present study.- Niu P.
- Schwarm A.
- Bonesmo H.
- Kidane A.
- Aspeholen Åby B.
- Storlien T.M.
- Kreuzer M.
- Alvarez C.
- Sommerseth J.K.
- Prestløkken E.
A basic model to predict enteric methane emission from dairy cows and its application to update operational models for the national inventory in Norway.
Animals (Basel). 2021; 11 (34202055)1891
Table 1Results of RMSPE and CCC analysis for CH4 production for the individual lactating cow database predicted using IPCC06 (Ym = 6.5%, DE = 70%), or using
IPCC, 2019a
as follows: 1YM (Ym = 5.7%, DE = 70%), 1YMIV (Ym = 5.7%, DE measured in vivo), 2YM (Ym = 5.7 or 6.0% depending on diet NDF, DE = 70%), or 2YMIV (Ym = 5.7 or 6.0% depending on diet NDF, DE measured in vivo)- IPCC
2019 Refinement to the 2006 IPCC Guidelines for National Greenhouse Gas Inventories. Chapter 10: Emissions from livestock and manure management.
in: IPCC Guidelines for National Greenhouse Gas Inventories Volume 4: Agriculture, forestry and other land use. Intergovernmental Panel on Climate Change,
2019
1
DE = energy digestibility; 69.8% for silage-based diets, 64.8% for hay-based diets. Ym = conversion factor of gross energy intake into enteric CH4 energy. RMSPE = root mean square prediction error expressed as a percentage of the observed mean. MB = error due to bias, as a percent of total RMSPE. SB = error due to slope, as a percent of total RMSPE. RB = error due to random bias as a percent of total RMSPE. CCC = concordance correlation coefficient, where CCC = r × Cb. r = Pearson correlation coefficient. Cb = bias correction factor. V = scale shift. μ = location shift. IPCC06 = first equation based on Intergovernmental Panel on Climate Change 2006 (IPCC, 2006b).
Model | Diet NDF (%) | DE (%) | Ym (%) | CH4 (g/d) | RMSPE (%) | MB (%) | SB (%) | RB (%) | CCC | r | Cb | V | μ |
---|---|---|---|---|---|---|---|---|---|---|---|---|---|
In vivo | — | — | — | 381 | — | — | — | — | — | — | — | — | — |
IPCC06 | — | 70 | 6.5 | 428 | 12.6 | 59.7 | 2.96 | 38.0 | 0.400 | 0.630 | 0.636 | 1.18 | −1.05 |
1YM | — | 70 | 5.7 | 370 | 9.40 | 9.69 | 1.88 | 89.9 | 0.579 | 0.630 | 0.919 | 1.35 | 0.296 |
1YMIV | — | 69.8 | 5.7 | 377 | 10.1 | 0.968 | 19.2 | 81.5 | 0.569 | 0.572 | 0.995 | 1.03 | 0.091 |
— | 64.8 | 5.7 | 377 | 10.1 | 0.968 | 19.2 | 81.5 | 0.569 | 0.572 | 0.995 | 1.03 | 0.091 | |
2YM | <35 | 70 | 5.7 | 377 | 9.27 | 0.928 | 2.46 | 98.2 | 0.536 | 0.572 | 0.937 | 1.42 | 0.095 |
>35 | 70 | 6.0 | 377 | 9.27 | 0.928 | 2.46 | 98.2 | 0.536 | 0.572 | 0.937 | 1.42 | 0.095 | |
2YMIV | <35 | 69.8 | 5.7 | 384 | 11.2 | 0.614 | 28.2 | 72.8 | 0.488 | 0.490 | 0.997 | 0.968 | −0.079 |
>35 | 64.8 | 6.0 | 384 | 11.2 | 0.614 | 28.2 | 72.8 | 0.488 | 0.490 | 0.997 | 0.968 | −0.079 |
1 DE = energy digestibility; 69.8% for silage-based diets, 64.8% for hay-based diets. Ym = conversion factor of gross energy intake into enteric CH4 energy. RMSPE = root mean square prediction error expressed as a percentage of the observed mean. MB = error due to bias, as a percent of total RMSPE. SB = error due to slope, as a percent of total RMSPE. RB = error due to random bias as a percent of total RMSPE. CCC = concordance correlation coefficient, where CCC = r × Cb. r = Pearson correlation coefficient. Cb = bias correction factor. V = scale shift. μ = location shift. IPCC06 = first equation based on Intergovernmental Panel on Climate Change 2006 (
IPCC, 2006b
).The IPCC06 was the less accurate model (MB = 59.7%; Cb = 0.636) with an overestimate of CH4 emission (428 vs. 381 g/d, respectively, for IPCC06 and in vivo). Overestimation of IPCC06 was also found by
Appuhamy et al., 2016
for North American diets. Similarly, Jiménez et al., 2021
showed that the calculation of CH4 using specific emission factors for tropical climate regions is better than the IPCC 2006 default emission factors, which overestimated CH4 production. The Ym value used in the IPCC06 equation (6.5%) explains the overestimation. To the best of our knowledge, only one study (Benaouda et al., 2020
) applied the IPCC refinement (IPCC, 2019a
) to estimate cattle enteric CH4 in Latin America, and the results showed that the new factors of - IPCC
2019 Refinement to the 2006 IPCC Guidelines for National Greenhouse Gas Inventories. Chapter 10: Emissions from livestock and manure management.
in: IPCC Guidelines for National Greenhouse Gas Inventories Volume 4: Agriculture, forestry and other land use. Intergovernmental Panel on Climate Change,
2019
IPCC, 2019a
made a substantial improvement in the prediction compared with the previous IPCC Tier 2 (- IPCC
2019 Refinement to the 2006 IPCC Guidelines for National Greenhouse Gas Inventories. Chapter 10: Emissions from livestock and manure management.
in: IPCC Guidelines for National Greenhouse Gas Inventories Volume 4: Agriculture, forestry and other land use. Intergovernmental Panel on Climate Change,
2019
IPCC, 2006b
). Given the importance of Ym, for large regions such as Europe, using specific values that better represent the characteristics of local production systems is advisable. For example, a recent study showed that the predicted Ym ranged from 6.22 to 6.72% for Norway (Niu et al., 2021
), whereas in the Netherlands, a Tier 3 approach used a predicted Ym ranging from 5.88% to 6.07% (- Niu P.
- Schwarm A.
- Bonesmo H.
- Kidane A.
- Aspeholen Åby B.
- Storlien T.M.
- Kreuzer M.
- Alvarez C.
- Sommerseth J.K.
- Prestløkken E.
A basic model to predict enteric methane emission from dairy cows and its application to update operational models for the national inventory in Norway.
Animals (Basel). 2021; 11 (34202055)1891
Bannink et al., 2011
).The 1YM model was the second least accurate model (MB = 9.69; Cb = 0.919) with underestimated emissions (370 vs. 381 g/d). The underestimation was mainly due to the hay-based diets: for these diets, the use of a fixed DE (70%) rather than the lower in vivo value (64.8%) reduced the predicted GE requirement (Equation 10.16) and hence the related CH4 emission.
The models 2YMIV, 2YM, and 1YMIV all have a MB below 1% and a Cb > 0.93. Using the in vivo DE improved the accuracy of the models, as confirmed also considering the V values. The models that best predicted the in vivo variability (V close to 1) were the ones using DE measured in vivo (V = 1.03 and 0.968 for 1YMIV and 2YMIV, respectively). The importance of DE as a key factor to estimate CH4 emissions has long been known, although it is not an easily measurable parameter. Other studies suggested the use of OM digestibility to predict CH4 emissions (
Ramin and Huhtanen, 2012
; Bell et al., 2016
). However, models that include OM digestibility had lower precision and accuracy than those based on DE (Benaouda et al., 2019
).- Benaouda M.
- Martin C.
- Li X.
- Kebreab E.
- Hristov A.N.
- Yu Z.
- Yáñez-Ruiz D.R.
- Reynolds C.K.
- Crompton L.A.
- Dijkstra J.
- Bannink A.
- Schwarm A.
- Kreuzer M.
- McGee M.
- Lund P.
- Hellwing A.L.F.
- Weisbjerg M.R.
- Moate P.J.
- Bayat A.R.
- Shingfield K.J.
- Peiren N.
- Eugène M.
Evaluation of the performance of existing mathematical models predicting enteric methane emissions from ruminants: Animal categories and dietary mitigation strategies.
Anim. Feed Sci. Technol. 2019; 255114207
In contrast, using in vivo DE values decreased the model precision (SB = 19.2 and 28.2% and r = 0.572 and 0.490 for 1YMIV and 2YMIV, respectively). The most precise model was 1YM (SB = 1.88%; r = 0.63) followed by IPCC06 (SB = 2.96; r = 0.63) and 2YM (SB = 2.46; r = 0.572). Given the differences in precision and accuracy, the CCC parameter has been suggested to simultaneously account for accuracy and precision (
Tedeschi, 2006
) and the CCC evaluation showed the highest (best) value for the model 1YM (0.579) followed, with a similar value, by 1YMIV (0.569).The in vivo values of DE and Ym for cows fed Mediterranean diets are shown in Table 2. The mean DE varied across study and diet (P < 0.001), with the hay diet studied by
Gislon et al., 2020b
having significantly lower DE values than all other studies and diets except the corn silage diet within the study of Colombini et al., 2015
. Several studies (i.e., Broderick, 1995
; Gislon et al., 2020b
) reported a higher DMI for cows fed hay-based diets than silage diets, and increased intake may increase passage rate, thus decreasing the DE (Gislon et al., 2020b
). Moreover, the hay-based diet was characterized by a low NDF digestibility, further decreasing energy digestibility. Overall, there is scarce information about the DE measured in vivo of cows fed hay diets. For example, Klevenhusen et al., 2011
reported a DE of 65.3% in lactating dairy cows fed a hay-based diet, a value comparable to the Italian hay diet data set.Table 2Digestible energy (DE; %) and Ym (%), and dietary NDF and starch (% DM) determined in cows fed Mediterranean diets with a different forage basis
Item | Study | SEM | P-value | ||||||
---|---|---|---|---|---|---|---|---|---|
Gislon et al., 2020b | Pirondini et al., 2015 | Colombini et al., 2015 | |||||||
CS | AGS | WS | Hay | CS | CS | SS | |||
DE | 73.6 | 72.6 | 70.3 | 64.8 | 68.6 | 65.2 | 68.9 | 1.30 | <0.001 |
Ym | 5.67 | 5.92 | 5.82 | 5.68 | 5.41 | 5.05 | 5.52 | 0.208 | 0.007 |
NDF | 32.8 | 27.1 | 37.7 | 36.6 | 33.0 | 36.6 | 36.4 | ||
Starch | 25.9 | 26.5 | 25.9 |
a–e Least squares mean values within a row with different superscripts differ significantly at P < 0.05.
1 CS = diet based on corn silage; AGS = diet based on alfalfa and grass silage; WS = diet based on wheat silage; Hay = diet based on hay; SS = diet based on sorghum silage; Ym = conversion factor of gross energy intake into enteric CH4 energy.
The mean Ym also varied across study and diet (P = 0.007), although to a lesser extent and with no individual forage standing out as being different to the others. The IPPC refinement Tier 2 (
IPCC, 2019a
) for the European data set (milk yield >8,500 kg) assumes a Ym value of 6.0%. This value is higher than the average value (5.58%) observed in the data set of the present study. The most studied forage type for Europe is grass silage, followed by corn silage and fresh forage; diets based on grass silage are expected to have higher Ym values than diets based on corn silage due to the higher starch content of the latter, which increases the dietary starch concentration (- IPCC
2019 Refinement to the 2006 IPCC Guidelines for National Greenhouse Gas Inventories. Chapter 10: Emissions from livestock and manure management.
in: IPCC Guidelines for National Greenhouse Gas Inventories Volume 4: Agriculture, forestry and other land use. Intergovernmental Panel on Climate Change,
2019
Hassanat et al., 2013
; Benchaar et al., 2014
). Based on these results, Mediterranean coefficients of Ym and DE appropriate for the MED equation were proposed as follows: Ym = 5.58% (average of Ym values for all diets, since no effect of diet was observed) and DE = 69.9 and 64.8% for silage- and hay-based diets, respectively.The cross-validation results on the independent data set (
Enriquez-Hidalgo et al., 2020
) are in Table 3. To highlight potential differences between equations, both a typical Mediterranean-type corn silage and hay-based diet, and a corn silage and a mixture of fresh annual ryegrass and berseem clover diet (MIX) were cross-validated. The average RMSPE (28.9%) was higher than that of the Italian data set (10.5%). Part of the reason for this discrepancy could be the method used to measure CH4 emission (respiration chambers for MED model, SF6 in the data set for cross-validation). Moreover, the SF6 technique collects the animal breath CH4 emissions, but not the rectum CH4 emissions, and this can partially explain the differences observed between the in vivo and the model estimations. - Enriquez-Hidalgo D.
- Teixeira D.L.
- Pinheiro Machado Filho L.C.
- Hennessy D.
- Toro-Mujica P.
- Williams S.R.O.
- Pereira F.C.
Incorporating a fresh mixed annual ryegrass and berseem clover forage into the winter diet of dairy cows resulted in reduced milk yield, but reduced nitrogen excretion and reduced methane yield.
Front. Vet. Sci. 2020; 7 (33330696)576944
Enriquez-Hidalgo et al., 2020
modified the SF6 technique according to - Enriquez-Hidalgo D.
- Teixeira D.L.
- Pinheiro Machado Filho L.C.
- Hennessy D.
- Toro-Mujica P.
- Williams S.R.O.
- Pereira F.C.
Incorporating a fresh mixed annual ryegrass and berseem clover forage into the winter diet of dairy cows resulted in reduced milk yield, but reduced nitrogen excretion and reduced methane yield.
Front. Vet. Sci. 2020; 7 (33330696)576944
Deighton et al., 2014
, to reduce the variability of CH4 yield estimation between cows, obtaining an accuracy similar to respiration chambers, although respiration chambers remain the gold standard method for measuring CH4 emission (Garnsworthy et al., 2019
).Table 3Results of RMSPE and CCC analysis on CH4 production for an independent individual lactating cow database derived from
Enriquez-Hidalgo et al., 2020
using - Enriquez-Hidalgo D.
- Teixeira D.L.
- Pinheiro Machado Filho L.C.
- Hennessy D.
- Toro-Mujica P.
- Williams S.R.O.
- Pereira F.C.
Incorporating a fresh mixed annual ryegrass and berseem clover forage into the winter diet of dairy cows resulted in reduced milk yield, but reduced nitrogen excretion and reduced methane yield.
Front. Vet. Sci. 2020; 7 (33330696)576944
IPCC, 2019a
1YM (Ym = 5.7%, DE = 70%) or the Tier 2 coefficients derived from the present study (MED, Ym = 5.58%; DE = 69.9%)- IPCC
2019 Refinement to the 2006 IPCC Guidelines for National Greenhouse Gas Inventories. Chapter 10: Emissions from livestock and manure management.
in: IPCC Guidelines for National Greenhouse Gas Inventories Volume 4: Agriculture, forestry and other land use. Intergovernmental Panel on Climate Change,
2019
1
DE = energy digestibility; 69.8% for silage-based diets, 64.8% for hay-based diets. Ym = conversion factor of gross energy intake into enteric CH4 energy. RMSPE = root mean square prediction error expressed as a percentage of the observed mean. DE = digestible energy. MB = error due to bias, as a percent of total RMSPE. SB = error due to regression, as a percent of total RMSPE. RB = error due to random bias, as a percent of total RMSPE. CCC = concordance correlation coefficient, where CCC = r × Cb. r = Pearson correlation coefficient. Cb = bias correction factor. V = scale shift. μ = location shift.
Model | CH4 (g/d) | RMSPE (%) | MB (%) | SB (%) | RB (%) | CCC | r | Cb | V | μ |
---|---|---|---|---|---|---|---|---|---|---|
Corn silage and alfalfa hay diet | ||||||||||
In vivo | 396 | |||||||||
1YM | 405 | 30.8 | 0.547 | 52.3 | 47.1 | 0.492 | 0.858 | 0.573 | 3.16 | −0.103 |
MED | 397 | 31.5 | 0.009 | 53.2 | 46.7 | 0.485 | 0.856 | 0.565 | 3.22 | −0.013 |
Corn silage and fresh annual ryegrass/berseem herbage (MIX) | ||||||||||
In vivo | 334 | |||||||||
1YM | 366 | 26.7 | 9.70 | 1.03 | 90.3 | 0.169 | 0.277 | 0.610 | 2.64 | −0.511 |
MED | 358 | 26.7 | 5.78 | 0.904 | 93.3 | 0.172 | 0.277 | 0.621 | 2.69 | −0.392 |
1 DE = energy digestibility; 69.8% for silage-based diets, 64.8% for hay-based diets. Ym = conversion factor of gross energy intake into enteric CH4 energy. RMSPE = root mean square prediction error expressed as a percentage of the observed mean. DE = digestible energy. MB = error due to bias, as a percent of total RMSPE. SB = error due to regression, as a percent of total RMSPE. RB = error due to random bias, as a percent of total RMSPE. CCC = concordance correlation coefficient, where CCC = r × Cb. r = Pearson correlation coefficient. Cb = bias correction factor. V = scale shift. μ = location shift.
The 1YM and MED had a higher (best) accuracy (evaluated in terms of MB) and CCC for corn silage and hay diet, rather than for the MIX diet. The predicted CH4 was higher than the in vivo value but differed between diets with a greater overestimation (evaluated in terms of µ) for the MIX diet (µ = −0.452, on average of 1YM and MED) than the corn silage + hay diet (−0.058, on average of 1YM and MED). The presence in the MIX diet of berseem clover, a forage containing natural substances able to reduce methanogenesis (
Enriquez-Hidalgo et al., 2020
), partially explains the IPCC method's overestimation of CH4 emissions for this diet compared with direct in vivo measurement. In contrast, the predicted CH4 emission on the corn silage diet using our proposed Mediterranean coefficients was very similar to that measured in vivo (396 and 397 g/d, respectively, for in vivo and MED).- Enriquez-Hidalgo D.
- Teixeira D.L.
- Pinheiro Machado Filho L.C.
- Hennessy D.
- Toro-Mujica P.
- Williams S.R.O.
- Pereira F.C.
Incorporating a fresh mixed annual ryegrass and berseem clover forage into the winter diet of dairy cows resulted in reduced milk yield, but reduced nitrogen excretion and reduced methane yield.
Front. Vet. Sci. 2020; 7 (33330696)576944
The RB value was higher for the MIX diet (91.8%, on average) than for the corn silage and hay diet (46.9%, on average); however, the r values were very low for the prediction of the MIX diet (on average, 0.277) compared with the corn silage and hay diet (on average, 0.857). The high RB for the MIX diet is explained by the lower values of predicted standard deviation compared with the in vivo and the low r value.
In conclusion, the study showed that IPCC 2019 predicts CH4 emission accurately and then it can be used as a tool for the prediction of CH4 emissions for inventories; however, for Mediterranean diets specific values of Ym and DE may be preferable, especially for hay-based diets.
Notes
The study was part of the project Life+Forage4Climate–Forage systems for less GHG emissions and more soil carbon sink in continental and Mediterranean agricultural areas (LIFE15 CCM/IT/000039), funded by the EU Commission (Brussels, Belgium). The authors thank the Department of Agricultural and Environmental Science, University of Milan, for funding the paper publication, Paolo Roveda (Università degli Studi di Milano, Milano, Italy) for the technical assistance with respiration chambers, and Vincenzo D'Ardes (Università degli Studi di Milano, Milano, Italy) for the calorimetric analyses.
D. Enriquez-Hidalgo was sponsored by Rothamsted Research's Institute Strategic Programme supported by the Biotechnology and Biological Sciences Research Council (BBS/E/C/000I0320; Rothamsted Research, North Wyke, Net Zero and Resilient Farming, Okehampton, Devon, United Kingdom).
The authors have not stated any conflicts of interest.
References
- Models for predicting enteric methane emissions from dairy cows in North America, Europe, and Australia and New Zealand.Glob. Chang. Biol. 2016; 22 (27148862): 3039-3056
- A model of enteric fermentation in dairy cows to estimate methane emission for the Dutch National Inventory Report using the IPCC Tier 3 approach.Anim. Feed Sci. Technol. 2011; 166–167: 603-618
- Modelling the effect of diet composition on enteric methane emissions across sheep, beef cattle and dairy cows.Animals (Basel). 2016; 6 (27618107): 54
- Development of mathematical models to predict enteric methane emission by cattle in Latin America.Livest. Sci. 2020; 241104177
- Evaluation of the performance of existing mathematical models predicting enteric methane emissions from ruminants: Animal categories and dietary mitigation strategies.Anim. Feed Sci. Technol. 2019; 255114207
- Methane production, digestion, ruminal fermentation, nitrogen balance, and milk production of cows fed corn silage- or barley silage-based diets.J. Dairy Sci. 2014; 97 (24359826): 961-974
- Prediction and improved estimation in linear models.John Wiley & Sons, 1977
- Performance of lactating dairy cows fed either alfalfa silage or alfalfa hay as the sole forage.J. Dairy Sci. 1995; 78 (7745152): 320-329
- Substitution of corn silage with sorghum silages in lactating cow diets: In vivo methane emission and global warming potential of milk production.Agric. Syst. 2015; 136: 106-113
- A modified sulphur hexafluoride tracer technique enables accurate determination of enteric methane emissions from ruminants.Anim. Feed Sci. Technol. 2014; 197: 47-63
- Incorporating a fresh mixed annual ryegrass and berseem clover forage into the winter diet of dairy cows resulted in reduced milk yield, but reduced nitrogen excretion and reduced methane yield.Front. Vet. Sci. 2020; 7 (33330696)576944
- Comparison of methods to measure methane for use in genetic evaluation of dairy cattle.Animals (Basel). 2019; 9 (31640130): 837
- Looking for high-production and sustainable diets for lactating cows: A survey in Italy.J. Dairy Sci. 2020; 103 (32113778): 4863-4873
- Milk production, methane emissions, nitrogen, and energy balance of cows fed diets based on different forage systems.J. Dairy Sci. 2020; 103 (32622607): 8048-8061
- Replacing alfalfa silage with corn silage in dairy cow diets: Effects on enteric methane production, ruminal fermentation, digestion, N balance, and milk production.J. Dairy Sci. 2013; 96 (23684039): 4553-4567
- Chapter 1: Introduction to the 2006 guidelines.in: IPCC Guidelines for National Greenhouse Gas Inventories Volume 1: General guidance and reporting. Intergovernmental Panel on Climate Change, 2006
- Chapter 10: Emissions from livestock and manure management.in: IPCC Guidelines for National Greenhouse Gas Inventories Volume 4: Agriculture, forestry and other land use. Intergovernmental Panel on Climate Change, 2006
- 2019 Refinement to the 2006 IPCC Guidelines for National Greenhouse Gas Inventories. Chapter 10: Emissions from livestock and manure management.in: IPCC Guidelines for National Greenhouse Gas Inventories Volume 4: Agriculture, forestry and other land use. Intergovernmental Panel on Climate Change, 2019
- 2019 Refinement to the 2006 IPCC Guidelines for National Greenhouse Gas Inventories. Overview Chapter.Intergovernmental Panel on Climate Change, 2019
- Concentrate supplementation on milk yield, methane and CO2 production in crossbred dairy cows grazing in tropical climate regions.J. Anim. Behav. Biometeorol. 2021; 92118
- Development of methane conversion factor models for Zebu beef cattle fed low-quality crop residues and by-products in tropical regions.Ecol. Evol. 2016; 6 (28725409): 7422-7432
- Enteric and manure-derived methane and nitrogen emissions as well as metabolic energy losses in cows fed balanced diets based on maize, barley or grass hay.Animal. 2011; 5 (22445412): 450-461
- A concordance correlation coefficient to evaluate reproducibility.Biometrics. 1989; 45 (2720055): 255-268
- Current trends in British dairy management regimens.J. Dairy Sci. 2014; 97 (25306285): 7985-7994
- Prediction of enteric methane production, yield, and intensity in dairy cattle using an intercontinental database.Glob. Chang. Biol. 2018; 24 (29450980): 3368-3389
- A basic model to predict enteric methane emission from dairy cows and its application to update operational models for the national inventory in Norway.Animals (Basel). 2021; 11 (34202055)1891
- Effect of dietary starch concentration and fish oil supplementation on milk yield and composition, diet digestibility, and methane emissions in lactating dairy cows.J. Dairy Sci. 2015; 98 (25465540): 357-372
- Development of an in vitro method for determination of methane production kinetics using a fully automated in vitro gas system—A modelling approach.Anim. Feed Sci. Technol. 2012; 174: 190-200
- Assessment of the adequacy of mathematical models.Agric. Syst. 2006; 89: 225-247
- World Greenhouse Gas Emissions: 2019.https://www.wri.org/data/world-greenhouse-gas-emissions-2019Date: 2019Date accessed: November 3, 2021
Article info
Publication history
Published online: February 09, 2023
Accepted:
November 17,
2022
Received:
April 12,
2022
Publication stage
In Press Corrected ProofIdentification
Copyright
© 2023
User license
Creative Commons Attribution (CC BY 4.0) | How you can reuse
Elsevier's open access license policy

Creative Commons Attribution (CC BY 4.0)
Permitted
- Read, print & download
- Redistribute or republish the final article
- Text & data mine
- Translate the article
- Reuse portions or extracts from the article in other works
- Sell or re-use for commercial purposes
Elsevier's open access license policy