The national inventories for greenhouse gas emissions (GHG) are considered the main instrument to formulate actions on national climate policies, such as the commitments for the Paris Agreement on climate change. However, these inventories have shown many shortcomings that need to be overcome using broader frameworks that incorporate key additional perspectives. GHG national inventories for livestock consider only direct biogenic emissions, so any mitigation measures ignore potential knock-on effects and trade-offs outside farm boundaries (e.g. pre-farm gate through fertiliser manufacturing or concentrate production) or through on-farm electricity/fuel usage. Life cycle assessment (LCA) approaches provide a more comprehensive framework to incorporate such off-farm and energy-related secondary effects. Although LCA-based approaches are very useful tools for climate and food policy, their large variation is to be considered when using average world-wide LCA-based GHG values (e.g. Gerber et al., Reference Gerber, Steinfeld, Henderson, Mottet, Opio, Dijkman, Falcucci and Tempio2013) for recommendations rather than values tailored to industrialised countries (Harwatt, Reference Harwatt2019). In order to be useful for country-level policies, LCA values must be specific enough to differentiate among production systems and countries.
There have been different studies attempting to estimate the contribution of European small ruminant systems to climate change via GHG emissions. Bellarby et al. (Reference Bellarby, Tirado, Leip, Weiss, Lesschen and Smith2013) examined overall GHG emissions in the EU27 livestock sector for the year 2007 and estimated GHG emissions from production and consumption of livestock products. They also reviewed available mitigation options and estimated their potential. Similarly, Weiss and Leip (Reference Weiss and Leip2012) investigated overall GHG emissions (as LCA from cradle to farm gate) from EU-27 livestock in the year 2004 and explored the implications for the total GHG emissions, including different estimations of land use and land use change associated CO2 emissions. Gerber et al. (Reference Gerber, Steinfeld, Henderson, Mottet, Opio, Dijkman, Falcucci and Tempio2013) investigated GHGs from livestock supply chains for one year (2010), providing disaggregated LCA-based results for different species (e.g. sheep and goats), product orientation (dairy v. meat), production systems (grassland v. mixed systems) and regions. To our knowledge, however, there is not, to date, a study that focuses on how European small ruminant dairy systems have contributed to global temperatures changes in an integrated manner, i.e. considering: (i) GHG emissions from a LCA approach, (ii) a historical time series evolution perspective and (iii) the effect on global temperatures rather than GHG emissions.
For the rest of the century, European small ruminant dairy systems, as well as other activities within the agricultural sector, are expected to pursue climate neutrality in order to contribute to the goal of the Paris Agreement on climate change regarding stabilising of global temperatures. Achieving climate neutrality in terms of metric-equivalent emissions, as is conventionally done for long-lived gases (e.g. CO2), would require that the sector progress towards net 0 cumulative CO2-e emissions by the year 2050. The Paris Agreement specifies that a balance between anthropogenic emissions by sources and removals by sinks of GHGs must be reached in the second half of this century.
Net zero is clear for CO2 and other long-lived GHGs, but not for CH4. Methane has a natural removal sink through chemical oxidation resulting in CH4 having a life of approximately 12 years in the atmosphere (IPCC, Reference Stocker, Qin, Plattner, Tignor, Allen, Boschung, Nauels, Xia, Bex and Midgley2013). Methane behaves very differently from CO2 such that if we maintain global CH4 emissions at the current level the concentrations of methane will then stabilise and not contribute further to global warming. This contrasts with CO2 where constant emissions lead to a constant warming rate and reduced emissions still lead to a reduced warming rate until they decline to zero. This problem can be avoided by re-defining climate neutrality (Ridoutt, Reference Ridoutt2020) in terms of ‘warming-equivalent’ emissions. In this sense, the new Global Warming Potential-star (GWP*) metric (Allen et al., Reference Allen, Shine, Fuglestvedt, Millar, Cain, Frame and Macey2018; Cain et al., Reference Cain, Lynch, Allen, Fuglestvedt, Frame and Macey2019; Lynch et al., Reference Lynch, Cain, Pierrehumbert and Allen2020) allows the expression of all emissions (i.e. long v. short-lived) in warming equivalents (CO2-we), which in a common cumulative framework correlates well with the resulting atmospheric temperature increase (Cain et al., Reference Cain, Lynch, Allen, Fuglestvedt, Frame and Macey2019).
In this paper we illustrate the contrasting climate impacts resulting from emissions of CH4, CO2 and N2O using GWP and GWP* methodologies. This has important implications for how we consider ‘zero emission’ or ‘climate neutral’ targets for sectors emitting different compositions of GHG.
Material and methods
Historical contribution of GHG emissions from dairy small ruminants in Europe (FAO region) and EU-27 on temperature change/warming
In order to estimate the historical warming associated with dairy sheep and goats in Europe, we need first to have their historical LCA-based GHG emissions. To our knowledge and to date, there is not such a calculation. There are, however, LCA-based GHG European calculations for single years (e.g. Weiss and Leip, Reference Weiss and Leip2012; Gerber et al., Reference Gerber, Steinfeld, Henderson, Mottet, Opio, Dijkman, Falcucci and Tempio2013) and specific LCA-based studies analysing the GHG emissions from different dairy production systems (as for Spain, e.g. in Batalla et al., Reference Batalla, Knudsen, Mogensen, Del Hierro, Pinto and Hermansen2015 or in Escribano et al., Reference Escribano, Elghannam and Mesias2020). In an attempt to derive historical European emissions, we used a combination of information: (i) LCA-based European studies for one year, (ii) Historical data on changes in productivity parameters and breeds for Spanish sheep and goats as a proxy of how European production systems have changed in the last decades and (iii) LCA-based GHG emissions for different Spanish production systems and breeds (full details are provided in the online Supplementary File: general approach in Section 1, (i) Section 2, (ii) and (iii) Section 3). We chose to develop two extrapolations in order to check if choosing different methods can affect the consistency of final results, but it is important to point out that producing an accurate historical GHG emissions calculation was beyond the scope of this study.
For assessing how this sector's emissions have contributed to global warming since 1990, we constructed a cumulative emissions framework on the basis that there is almost a linear relationship between the amount of CO2 emitted and the resulting atmospheric temperature increase (IPCC, Reference Stocker, Qin, Plattner, Tignor, Allen, Boschung, Nauels, Xia, Bex and Midgley2013). This is the basis that most countries are adopting for reductions and targets of ‘net zero’ CO2 emissions by 2050. Two approaches were used and subsequently compared: (i) Calculating emissions as CO2-e and assuming that all the 3 GHG accumulate in the atmosphere and (ii) an alternative approach, where temperatures for long-lived GHG gases (i.e. CO2 and N2O) depend on the accumulation of emissions (as approach 1), but for short-lived gases (i.e. CH4) temperature would depend on the rate of emissions. For the latter alternative approach, we used the GWP* approach developed by Allen et al. (Reference Allen, Shine, Fuglestvedt, Millar, Cain, Frame and Macey2018), which includes the differentiation rate and uses for cumulative metrics the CO2-we (CO2-warming equivalent as defined by Cain et al., Reference Cain, Lynch, Allen, Fuglestvedt, Frame and Macey2019) instead of CO2-e. A simple coefficient known as TRCE (transient climate response to cumulative carbon emissions) can then be multiplied by cumulative CO2-we to obtain an approximate estimate of temperature change due to the change in CO2-we burden experienced, as explained by Lynch et al. (Reference Lynch, Cain, Pierrehumbert and Allen2020). The TCRE metric represents the robust near-linear relationship between global warming and cumulative CO2 emissions. Essentially, the TCRE estimates for every 42 GtCO2 emitted, the global temperature will increase by between 0.009 and 0.011°C (Rogelj et al., Reference Rogelj, Shindell, Jiang, Fifita, Forster, Ginzburg, Handa, Kheshgi, Kobayashi, Kriegler, Mundaca, Séférian, Vilariño, Masson-Delmotte, Zhai, Pörtner, Roberts, Skea, Shukla, Pirani, Moufouma-Okia, Péan, Pidcock, Connors, Matthews, Chen, Zhou, Gomis, Lonnoy, Maycock, Tignor and Waterfield2018).
The GWP* methodology (full details are given in the Supplementary File: Section 4) was applied for the species, regions and extrapolations shown in Table 1.
Details on the species, regions and extrapolations basis for the different historic scenarios
Description of policy scenarios on GHG mitigation from European and EU-27 dairy small ruminant activity
We forecasted different pathways of GHG emissions from the European small ruminants' dairy system (no changes, 0.4, 0.6 and 0.8% annual GHG reductions) and assessed how these scenarios would be contributing to the global temperature targets for the 2018–2100 period using the cumulative CO2-we framework approach previously described and the TRCE coefficient. We carried out this exercise for the 6 different combinations that were used to calculate the historical contribution. It must be noted that although these scenarios do not represent concrete measures to reduce GHG emissions, they could be used to interpret known measures. For illustration purposes and as an example of these potential mitigation measures, we simulated for dairy sheep production systems in the European FAO region the gradual introduction of a feed additive (3-NOP) that can reduce about 30% CH4 from enteric sources in sheep (Martínez-Fernández et al., Reference Martínez-Fernández, Abecia, Arco, Cantalapiedra-Hijar, Martín-García, Molina-Alcaide, Kindermann, Duval and Yáñez-Ruiz2014), and do so without loss in efficiency (Hristov et al., Reference Hristov, Oh, Giallongo, Frederick, Harper, Weeks, Branco, Moate, Deighton, Williams, Kindermann and Duval2015; Lopes et al., Reference Lopes, de Matos, Harper, Giallongo, Oh, Gruen, Ono, Kinderman, Duval and Hristov2016). We assumed this 30% reduction as potentially achievable and used SHEEP-EUR-1 for 4 different scenarios varying in timeline for full introduction of the 3-NOP as a mitigation measure (20, 40, 60 and 80 years).
For all the scenarios mentioned above, we included as an alternative emissions projection, that pathway of net GHG emissions that considers CO2 offsets due to estimations of potential soil organic carbon (SOC) sequestration as approximated using data from Batalla et al. (Reference Batalla, Knudsen, Mogensen, Del Hierro, Pinto and Hermansen2015) (SOC-1), and Batalla et al. (Reference Batalla, Knudsen, Mogensen, Del Hierro, Pinto and Hermansen2015) and Escribano et al. (Reference Escribano, Elghannam and Mesias2020) (SOC-2) (full details are given in the online Supplementary File in Section 3). This approach estimates the annual SOC that can potentially be sequestered averaged over 100 years following the approach by Petersen et al. (Reference Petersen, Knudsen, Hermansen and Halberg2013). This approximation is subject to large uncertainties and requires several assumptions such as: no-tillage is carried out, SOC is not at its level of C saturation and nor does it consider climatic and soil factors. For our projections of SOC between 2020 and 2100, we are assuming that the annual C sequestration estimated for the 2018 year could be maintained until the end of the century.
Results
Contribution of historical European small ruminant activity to global temperature change (1990–2018)
Europe (FAO region): When comparing CO2-e and CO2-we for the different species and extrapolations for European dairy sheep (Fig. 1, top) and dairy goats (Fig. 1, bottom), we observe that GHG emissions differ considerably when expressed as CO2-e from when expressed as CO2-we. This is true for both annual (Fig. 1, left) and aggregated (cumulative: Fig. 1, right)values. Correspondingly, warming estimations reflect that the CO2-we metric results in much lower temperatures than using CO2-e. In addition, whereas the extrapolation method influenced the total cumulative CO2-we, it did not affect the cumulative GHG emissions expressed as CO2-e.
Corresponding annual (left panels) and cumulative (right panels) GHG emissions from dairy sheep (a:d) and goat (e:h) systems in Europe (FAO region). Solid lines: emissions calculated using the global warming potential GWP100 metric (CO2-e). Dotted lines: emissions calculated using the GWP* metric (CO2 warming equivalent, CO2-we). Values use extrapolation of GHG emissions based on Batalla et al. (Reference Batalla, Knudsen, Mogensen, Del Hierro, Pinto and Hermansen2015) (EUR-1) and Batalla et al. (Reference Batalla, Knudsen, Mogensen, Del Hierro, Pinto and Hermansen2015) and Escribano et al. (Reference Escribano, Elghannam and Mesias2020) (EUR-2).
EU-27: Similar to the previous case, a comparison of CO2-e and CO2-we for aggregated small ruminant systems in the current European Union member countries (SR-EU27-1) for the period between 1990 and 2018 also showed great differences depending on the way we express GHG emissions and associated warming (Supplementary Fig. S13).
Impacts of GHG mitigation scenarios from European small ruminant activity on global temperature change (2020–2100)
When observing the scale of warming that would result from different GHG emission reduction pathways in sheep (Fig. 2, top) and goat (Fig. 2, bottom) dairy production systems in the European FAO region, data show that keeping constant GHG emissions up until the end of the century would contribute to some net warming effect. Reducing annual GHG emissions reduces this warming effect. If annual GHG emissions are below 0.8% (dairy sheep) or 1% (dairy goats) per year, no additional warming to the atmosphere is observed in the simulation. More ambitious GHG reduction pathways would contribute to further cooling of the atmosphere in this period. When we considered emissions offsets due to potential SOC sequestration as an alternative calculation (scenarios + SOC), the pathway showed much lower associated warming than without consideration of SOC stock changes.
Warming resulting from different GHG emission reductions pathways considering the full life cycle assessment for sheep milk production (top panels) and goat milk production (bottom panels) in Europe (FAO region) in the period 2020–2100. Values show no change (left) and then 0.4, 0.6 and 0.8% annual reduction in total GHG emissions. Black lines: values use extrapolation of GHG emissions based on Batalla et al. (Reference Batalla, Knudsen, Mogensen, Del Hierro, Pinto and Hermansen2015) (EUR-1). Grey lines: values use Batalla et al. (Reference Batalla, Knudsen, Mogensen, Del Hierro, Pinto and Hermansen2015) and Escribano et al. (Reference Escribano, Elghannam and Mesias2020) (EUR-2). Solid lines are without including potential SOC sequestration, dotted lines include sequestration.
If this same observation is done at the level of the current European Union member states for aggregated small ruminant data (Supplementary Fig. S14), keeping constant GHG emissions from small ruminants' dairy production systems in the EU-27 up until the end of the century would contribute to some net warming effect. Reducing annual GHG emissions reduces this warming effect. If annual GHG emissions are reduced beyond 1.2%, no additional warming to the atmosphere was observed in the simulation. Again, more ambitious GHG reduction pathways would contribute to further cooling. There are no differences between the two approaches for extrapolation of GHG.
Mitigation example (3-NOP feed additive introduce to the European dairy sheep systems)
As expected, the faster the gradual inclusion of the 3-NOP additive in the diet was, the larger the reduction was on the warming effect of each of the dairy sheep scenarios for the European FAO region (Fig. 3). In particular, when the full implementation of the measure happened as fast as in 20 years, European dairy sheep systems would not add extra warming for the first 40 years (until 2060: Fig. 3, left). Including the potential SOC sequestration results in any of the scenarios of 3-NOP feed additive inclusion would mean there would be no additional warming for the years to come up until 2100 (Fig. 3, right).
Warming resulting from gradually (in 20, 40, 60 and 80 years) introducing the 3-NOP feed additive in dairy sheep production systems in the European FAO region for the period 2020–2100 for scenarios without (left) and including potential SOC sequestration (right). Values use extrapolation of GHG emissions based on Batalla et al. (Reference Batalla, Knudsen, Mogensen, Del Hierro, Pinto and Hermansen2015) (SHEEP-EUR-1).
Discussion
Although a number of studies indicate that milk production in Europe has led to large GHG emissions (e.g. Lesschen et al., Reference Lesschen, van den Berg, Westhoek, Witzke and Oenema2011) and there have been urgent calls to drastically reduce milk consumption (Harwatt, Reference Harwatt2019), our estimates indicate that, in terms of historical warming (1990–2018), their conclusions do not apply to small ruminants. When linking emissions to their resulting temperature we found that for both EU-27 and European FAO region dairy small ruminants have not caused additional warming in the 1990–2018 period. Conventional application of CO2 metric-equivalent emissions using GWP100 results in much larger estimated CO2 equivalent emissions than through the use of CO2 warming-equivalent emissions using GWP*.
We also observed very different trends in cumulative GHG emissions when considering the conventional GWP metric against the new GWP* metric. The differences between CO2-we and CO2-e emissions are caused by CO2-we capturing the weighting given to the reduction in CH4 emissions rate in the preceding decades which leads to a temperature decrease. In this sense, we can see that the use of CO2-e over-states the impact of CH4 emissions on this particular historic emissions trajectory.
Greenhouse gas emission trends were consistent across the different methodologies used. As expected, the use of different models (GLEAM-EUR v. CAPRI-EU27) and extrapolations (e.g. EUR1 v. EUR2) for the historical analysis resulted in different absolute GHG and associated warming values. However, their trends were similar for all their combinations, which allows us to think that our findings are valid in spite of the limitations in data availability. For the future scenarios, the differences amongst extrapolations were very small, which reinforces the robustness of this approach.
Our extrapolations indicate that there has been a reduction in total GHG emissions and a change in GHG emissions profile of the European and EU-27 C footprint of the small ruminants' milk. This is in spite of milk production from small ruminants in Europe having increased steadily, ca. 10 and 25% for sheep and goat milk production in the period 1990−2018 (data sourced from FAOSTAT). The larger expansion of goat dairy systems compared to sheep ones has resulted in differences in warming impact between the two: whereas changes in dairy goat GHG emissions have resulted in some level of increased warming, larger emission reductions in dairy sheep systems have compensated it.
In particular, there seems to be a large reduction in CH4 emissions intensity (i.e. CH4/kg milk), and in N2O emissions intensity (i.e. N2O/kg milk) and an increase in CO2 emissions intensity (i.e. CO2/kg milk, derived from fossil fuels use and land use change) (detailed in online Supplementary File in Section 5, Supplementary Fig. S7). Improvements in reproduction management, which have led to greater reproduction efficiency (smaller herds required per unit of milk produced) and a certain intensification linked with breeds that are more productive as well as higher fodder inputs of superior quality are expected to have caused both the reduction of the overall GHG footprint and the shift in GHG emission forms.
The shift between smaller CH4 (short lived) to larger CO2 (long lived) emissions implies that a fast cooling effect may be seen (reduction of CH4) at the expense of, to a certain extent, some increased level of warming for the long-run (increase of CO2). However, since reduction rates in CH4 intensity have been much larger than CO2 increases, such trade-off seems to have had limited warming implications for the future in this particular case. This should, nevertheless, be taken cautiously: the reduction in CH4 emissions by livestock no longer grazing on marginal areas may have been replaced by methane-emitting wild herbivores that would be occupying their former grazing niches in the ecosystem (Manzano and White, Reference Manzano and White2019). It is known that European wild herbivores are expanding into abandoned marginal grazing areas (Milner et al., Reference Milner, Bonenfant, Mysterud, Gaillard, Csányi and Stenseth2006), so the net amount of CH4 arising from the whole wild and domestic grazing system may have stayed around the same values.
Higher levels of intensification are expected to exacerbate this trade-off in emissions of short-lived gases to long-lived gases. This trade-off may be further imbalanced considering that N2O (a long-lived gas) has been found to follow a non-linear response to N inputs for both temperate (Cardenas et al., Reference Cardenas, Thorman, Ashlee, Butler, Chadwick, Chambers, Cuttle, Donovan, Kingston, Lane, Dhanoa and Scholefield2010) and tropical (Tully et al., Reference Tully, Abwanda, Thiong'o, Mutuo and Rosenstock2017) systems. Higher levels of intensification will lead to higher N2O emission rates, thus potentially increasing the amount of long-term warming effect. This effect has not only been recognised in the literature, but it has been indicated as potentially implementable in national inventories of GHG in the last IPCC Refinement (IPCC, 2019).
Scenarios for the small ruminant sector in Europe (for both sets of countries considered) to achieve zero warming effect for the next years to come in the twenty-first century involve modest (0.8–1.2%) annual GHG reductions considering the LCA sectorial boundaries. This approach, based on real warming rather than on CO2-e (metric-equivalent) neutrality, could induce changes in GHG reduction strategies in compliance with the Paris Agreement if compared to those arising from within a cumulative carbon framework (Lynch et al., Reference Lynch, Cain, Pierrehumbert and Allen2020). Strategies prioritising long-term climatic effects would be less concerned about biogenic CH4 targets, and design stronger efforts to reduce fossil CO2 emissions. Conversely, other strategies focused on short-term methane action, which trust effective carbon capture and storage technologies for the distant future, could be reinforced by using warming-equivalent emission metrics.
Rather than implying that these GHG reductions would require parallel reductions in milk production, mitigation measures that have been widely studied (e.g. Hristov et al., Reference Hristov, Oh, Firkins, Dijkstra, Kebreab, Waghorn, Makkar, Adesogan, Yang, Lee, Gerber, Henderson and Tricarico2013) and are already applicable could also contribute to these reductions. For example, there are some feed additives that can lower the CH4 output from the rumen. In this sense, our results for dairy sheep indicate that an introduction of the feed additive 3-NOP, assuming its full potential was achieved as a long-lasting effect, would be sufficient to reduce additional warming to close to zero in the next decades (i.e. be climatically neutral) and to maintain the same European sheep milk volume and production systems structure as those from 2018. This is especially true if we considered the potential for SOC sequestration in our calculations and/or we introduced this 3-NOP measure at a fast pace of 5% per year (in volume of milk produced). Interestingly, the speed of introduction makes a big difference for the near-term transition period but much less of a difference by the end of the century, suggesting that a slower introduction of 3-NOP would lead to almost the same, but delayed for some decades, temperature stabilisation. This would contrast with any measure that targets a reduction in a long lasting gas (e.g. CO2), which would result in much larger temperature stabilisation should its introduction be delayed. There are other specific examples of promising feed additives, such as those from red algae extracts (Kinley et al., Reference Kinley, Martinez-Fernandez, Matthews, de Nys, Magnusson and Tomkins2020) that have resulted in up to a 50% reduction in CH4 emissions.
In general, there are many other potential mitigation measures to reduce GHG emissions for small ruminants in temperate regions (Jones et al., Reference Jones, Jones and Cross2014). Those measures that improve animal fertility or longevity, which have an impact on the replacement rate of the herd, will reduce overall GHG emissions intensity (Jones et al., Reference Jones, Jones and Cross2014).There is also certain scope to reduce GHG emissions via improvement of animal health, animal productivity (through animal data interpretation for example, Belanche et al., Reference Belanche, Martín-García, Fernández-Álvarez, Pleguezuelos, Mantecón and Yáñez-Ruiz2019) and changes in the animal diet. Measures that improve the overall efficiency of the systems, such as the replacement of conventional feed by agricultural by-products (Romero-Huelva and Molina-Alcaide, Reference Romero-Huelva and Molina-Alcaide2013; Marcos et al., Reference Marcos, Carro, Fernández Yepes, Haro, Romero-Huelva and Molina-Alcaide2020), have been found to significantly reduce the total C footprint of the milk from small ruminants' systems too (Pardo et al., Reference Pardo, Martin-Garcia, Arco, Yáñez-Ruiz, Moral and Del Prado2016).
Grasslands, and consequently grassland-based small ruminant systems, have the additional advantage of acting as a carbon sink, which can be exploited for GHG mitigation (Stanley et al., Reference Stanley, Rowntree, Beede, DeLonge and Hamm2018). Different European studies (e.g. Salvador et al., Reference Salvador, Corazzin, Romanzin and Bovolenta2017; Eldesouky et al., Reference Eldesouky, Mesias, Elghannam and Escribano2018; Horrillo et al., Reference Horrillo, Gaspar and Escribano2020; Sabia et al., Reference Sabia, Kühl, Flach, Lambertz and Gauly2020) indicate that considering the potential of on-farm SOC sequestration could offset the C footprint of sheep and goat milk by up to 90%. In this study, we estimated values of SOC sequestration potential for current small ruminant systems in Europe that could offset about 17–23% of the C footprint, depending on the assumptions chosen. This SOC sequestration could partly offset CO2 emissions from fossil fuels use and land use change. Our approximation indicates that, if SOC sequestration is considered and assumed to remain at our estimated current levels (i.e. offsetting 16% or 23% of total C footprint), reductions of 0.5% or 0.6% GHG annually would be sufficient to avoid adding any additional warming in the period 2020–2100 for European sheep and goat systems, respectively. Actual SOC sequestration is, however, still subject to large uncertainties and our approximation is not meant to be considered as a robust SOC value but, rather, as a useful illustration of the importance of improving the SOC sequestration estimates. Soil organic C (SOC) tends to have large accumulation rates during the first 20 years after a change in management (assuming a particular management leads to SOC sequestration), but the rate of change slows down considerably afterwards (Powlson et al., Reference Powlson, Whitmore and Goulding2011). We have not considered this or the fact that SOC sequestration benefits can be reversed if management changes to a situation where soil is degraded or ploughed up for arable systems.
Our potential SOC sequestration approximation for the historic period of 1961–2018 indicated a large reduction of potential SOC sequestration in time (about 80% expressed as kg CO2/L milk) driven by a change in production systems structure. As the sector has intensified and improved both the quality of the feeding and the reproductive performance of the animals, enteric-CH4 has shrunk. However, this has been achieved at a certain expense of SOC sequestration, since part of the feeding improvement has been caused by increasing housing periods and decreasing reliance on permanent grasslands, especially for grazing. This type of pollution swapping issue has been previously described for cattle dairy systems by Vellinga and Hoving (Reference Vellinga and Hoving2011).
In conclusion, our study illustrates with a real example how GWP* provides an improved method to assess the impact of methane, creates actionable goals to achieve the mandated maximum 1.5°C increase in average temperature by 2050, and evaluates current policies proposed for such goals. In our case, we used a timeline for 2100. As has been highlighted by Lynch et al. (Reference Lynch, Cain, Pierrehumbert and Allen2020) and shown in our study, using GWP100 does not even indicate the correct direction of temperature change during periods of declining CH4 emissions, and new GWP* metrics offer a much more appropriate picture. This is especially important when assessing medium/long term (>100 years) periods and when comparing production systems and species with different GHG profiles. Although GWP* has raised questions of equity and fairness for its application in international climate policies (Rogelj and Schleussner, Reference Rogelj and Schleussner2019; Smith and Balmford, Reference Smith and Balmford2020), the analysis of such criticisms is beyond the scope of our paper. However, most countries with large historic agricultural CH4 emissions that are capable of reducing their additional contribution to warming (i.e. negative cumulative CO2-we), which are the centrepiece for the fairness argument, could still maintain their historically high contribution to global warming at a steady rate. This would be equivalent to a country achieving net zero CO2 emissions, as CO2 remains in the atmosphere for millennia. However, the difference is that this CH4 source is an opportunity to undo past contributions to global warming, the same way as CO2 removal and C storage would be. Moreover, these countries have abundant natural grazing ecosystems in their territory that can host many wild methane-producing herbivores. Total abandonment of grazing in natural rangelands is likely to be a very ineffective climate change policy (Manzano and White, Reference Manzano and White2019).
Supplementary material
The supplementary material for this article can be found at https://doi.org/10.1017/S0022029921000157.
Acknowledgements
This research is supported by the Spanish Government through María de Maeztu excellence accreditation 2018-2022 (Ref. MDM-2017-0714) and by the Basque Government through the BERC 2018-2021 programme. This work was also supported by the Horizon2020 SFS-01c-2015 project entitled ‘Innovation of sustainable sheep and goat production in Europe (iSAGE)’ (grant number 679302). Agustin del Prado is financed by the programme Ramon y Cajal from the Spanish Ministry of Economy, Industry and Competitiveness (RYC-2017-22143). Pablo Manzano is funded by a post-doctoral fellowship at the Helsinki Institute for Sustainability Science (HELSUS). The authors would like to acknowledge the work of two anonymous reviewers in helping us to improve the article.
Open access
