Advertisement
Research Article|Articles in Press

Taking the diet of cows into consideration in designing payments to reduce enteric methane emissions on dairy farms.

Open AccessPublished:May 23, 2023DOI:https://doi.org/10.3168/jds.2022-22766

      ABSTRACT

      Enteric fermentation from dairy cows is a major source of methane. Significantly and rapidly reducing those emissions would be a powerful lever to mitigate climate change. For a given productivity level, introducing fodder with high sources of omega 3 content such as grass or linseed in the feed ration of dairy cows both improves the milk nutritional profile and reduces enteric methane emissions per liter. Changing cows' diet may represent additional costs for dairy farmers and calls for the implementation of payments for environmental services to support the transition. This paper analyzes 2 design elements affecting the effectiveness of a payment conditioned to the reduction of enteric methane emissions: (i) the choice of emission indicator capturing the effect of farmers' practices (ii) the payment amount relative to the extra milk production costs incurred. Using representative farm-level economic data from the French farm accountancy data network, we compare enteric methane emissions per liter of milk calculated with an Intergovernmental Panel on Climate Change Tier 2 method, to baseline emissions from a Tier 3 method accounting for diet effects. We also quantify the extra milk production costs of integrating more grass in the fodder systems by estimating variable cost functions for different dairy systems in France. Our results show the relevance of using an emission indicator sensitive to diet effects, and that the significance and direction of the extra-costs for producing milk with more grass differ according to the production basin and the current share of grasslands in the fodder crop rotation. We stress the importance of developing payments for environmental services with well-defined environmental indicators accounting for the technical problem addressed, and the need to better characterize heterogeneous funding requirements for supporting a large-scale adoption of more environment-friendly practices by farmers.

      Key words

      INTERPRETIVE SUMMARY: This study provides insights on the design of effective payments for environmental services conditioned to enteric methane emissions reduction. Two indicators of enteric emissions were compared, and the extra production costs for adding more grass in fodder systems as an emission reduction lever were computed for representative French medium and large dairy farms. Results show that an emission indicator accounting for diet effects is crucial to accurately monitor farmers' effort. Moreover, farms face heterogeneous extra-costs for changing cows' feed ration, and it should be taken into account in the definition of the payment amount to ensure a large-scale adoption.

      INTRODUCTION

      The agricultural sector is a major source of greenhouse gases (GHG), accounting for 10% of EU-KPs (European Union, United Kingdom and Iceland) in 2018 (
      • Citepa
      Inventaire des émissions de polluants atmosphériques et de gaz à effet de serre en France - Format Secten. National inventory of GHG and air pollutant emissions in France (in French).
      ;
      • EEA
      Annual European Union greenhouse gas inventory 1990 – 2018 and inventory report 2020. Submission to the UNFCCC Secretariat.
      ). 81% of agricultural methane emissions result from enteric fermentation, and 39% of those 81% are produced by dairy cows (
      • EEA
      Annual European Union greenhouse gas inventory 1990 – 2018 and inventory report 2020. Submission to the UNFCCC Secretariat.
      ). Methane is the second contributor to radiative forcing. Currently, the global warming potential of methane is set at 28 times higher than the global warming potential of carbon dioxide over 100 years and 84 over 20 years (
      • Myhre G.
      • Shindell D.
      • Bréon F.-M.
      • Collins W.
      • Fuglestvedt J.
      • Huang J.
      • Koch D.
      • Lamarque J.-F.
      • Lee D.
      • Mendoza B.
      • Nakajima T.
      • Robock A.
      • Stephens G.
      • Takemura T.
      • Zhang H.
      Anthropogenic and Natural Radiative Forcing.
      ). As methane is a short-lived climate pollutant continuously destroyed in the atmosphere, its effect on climate change depends mostly on short-term emission rates. In theory, decreasing the methane emissions rate below its natural destruction rate would have a cooling effect (
      • Cain M.
      • Lynch J.
      • Allen M.R.
      • Fuglestvedt J.S.
      • Frame D.J.
      • Macey A.H.
      Improved calculation of warming-equivalent emissions for short-lived climate pollutants.
      ). Therefore, a significant reduction in methane emissions, in particular from agricultural activities, would rapidly mitigate climate change and is a powerful lever to meet the European Union's 2050 climate targets (
      • Dupraz P.
      Policies for the ecological transition of agriculture: the livestock issue.
      ).
      Enteric fermentation is identified as the first source of GHG emissions from dairy farms in both developed and developing countries (
      • Jayasundara S.
      • Worden D.
      • Weersink A.
      • Wright T.
      • VanderZaag A.
      • Gordon R.
      • Wagner-Riddle C.
      Improving farm profitability also reduces the carbon footprint of milk production in intensive dairy production systems.
      ;
      • Wilkes A.
      • Wassie S.
      • Odhong' C.
      • Fraval S.
      • van Dijk S.
      Variation in the carbon footprint of milk production on smallholder dairy farms in central Kenya.
      ). The quantity of methane produced during the digestive process of ruminants depends highly on the characteristics of the animal itself, such as the breed, body weight and age (
      • Gavrilova O.
      • Leip A.
      • Dong H.
      • Macdonald J.D.
      • Gomez Bravo C.A.
      • Amon B.
      • Barahona Rosales R.
      • Agustin del Prado A.
      • Aparecida de Lima M.
      • Oyhantcabal W.
      • van der Weerden T.J.
      • Widiawati Y.
      Emissions from Livestock and Manure Management.
      ). However, enteric methane emissions are also directly related to farming practices, in particular the amount of feed intake and composition, and the proportion of carbohydrates that feed ration contains (
      • Martin C.
      • Morgavi D.P.
      • Doreau M.
      Methane mitigation in ruminants: from microbe to the farm scale.
      ). In particular, for a given productivity level, enteric methane emissions decline as dairy cow feed is enriched with α-linolenic acid ALA (polyunsaturated fatty acid of the omega-3 family), for which the main natural sources are linseed and grass fodders (
      • Dong Y.
      • Bae H.D.
      • McAllister T.A.
      • Mathison G.W.
      • Cheng K.J.
      Lipid induced depression of methane production and digestibility in the artificial rumen system (RUSITEC).
      ;
      • Martin C.
      • Morgavi D.
      • Doreau M.
      • Jouany J.P.
      Comment réduire la production de méthane chez les ruminants? How can the production of methane by ruminants be reduced?.
      ,
      • Martin C.
      • Rouel J.
      • Jouany J.P.
      • Doreau M.
      • Chilliard Y.
      Methane output and diet digestibility in response to feeding dairy cows crude linseed, extruded linseed, or linseed oil.
      ,
      • Martin C.
      • Morgavi D.P.
      • Doreau M.
      Methane mitigation in ruminants: from microbe to the farm scale.
      ,
      • Martin C.
      • Promiès D.
      • Ferlay A.
      • Rochette Y.
      • Martin B.
      • Chilliard Y.
      • Morgavi D.
      • Doreau M.
      Methane output and rumen microbiota in dairy cows in response to long-term supplementation with linseed or rapeseed of grass silage or pasture based diets.
      ;
      • Chilliard Y.
      • Martin C.
      • Rouel J.
      • Doreau M.
      Milk fatty acids in dairy cows fed whole crude linseed, extruded linseed, or linseed oil, and their relationship with methane output.
      ;
      • Grainger C.
      • Beauchemin K.A.
      Can enteric methane emissions from ruminants be lowered without lowering their production?.
      ). Moreover, as productivity per cow increases, methane emissions per kilogram of milk decrease (
      • Martin C.
      • Morgavi D.
      • Doreau M.
      • Jouany J.P.
      Comment réduire la production de méthane chez les ruminants? How can the production of methane by ruminants be reduced?.
      ). Animal productivity can also be improved through nutrition, as well as through herd's management and genetics (
      • Boadi D.
      • Benchaar C.
      • Chiquette J.
      • Massé D.
      Mitigation strategies to reduce enteric methane emissions from dairy cows: Update review.
      ). To accurately monitor the evolution of enteric methane emissions on dairy farms, one must consider both dimensions (productivity and feeding). Authors show enteric methane emissions can differ significantly from one indicator to another and recommend considering both production intensity and feed usage for more accurate estimates (
      • Hagemann M.
      • Hemme T.
      • Ndambi A.
      • Alqaisi O.
      • Sultana M.N.
      Benchmarking of greenhouse gas emissions of bovine milk production systems for 38 countries.
      ;
      • Sauvant D.
      • Giger-Reverdin S.
      • Serment A.
      • Broudiscou L.
      Influences des régimes et de leur fermentation dans le rumen sur la production de méthane par les ruminants.
      ).
      Changing management practices, and in particular cows diet, to decrease enteric methane emissions may be costly for many farmers. Economic incentives are currently being developed in the agricultural sector to support the reduction of GHG emissions. While studies show taxes and emission permits are the most efficient mitigation policy options, they are unpopular tools among producers and consumers (
      • Key N.
      • Tallard G.
      Mitigating methane emissions from livestock: A global analysis of sectoral policies.
      ;
      • Douenne T.
      • Fabre A.
      French attitudes on climate change, carbon taxation and other climate policies.
      ;
      • Funke F.
      • Mattauch L.
      • van den Bijgaart I.
      • Godfray H.C.J.
      • Hepburn C.
      • Klenert D.
      • Springmann M.
      • Treich N.
      Toward Optimal Meat Pricing: Is It Time to Tax Meat Consumption?.
      ). The agricultural sector remains exempted from carbon pricing schemes in most developed countries (
      • Wirsenius S.
      • Hedenus F.
      • Mohlin K.
      Greenhouse gas taxes on animal food products: Rationale, tax scheme and climate mitigation effects.
      ;
      • Henderson B.
      • Verma M.
      Global assessment of the carbon leakage implications of carbon taxes on agricultural emissions.
      ). In the EU, the agri-environmental policy has historically be designed as positive support instruments (agri-environmental schemes) (
      • Baylis K.
      • Peplow S.
      • Rausser G.
      • Simon L.
      Agric.-environmental policies in the EU and United States: A comparison.
      ). Alternatively, there is growing interest in supporting financially changes of management practices through payments for environmental services (PES), as more acceptable tools on the short-term. Targeting enteric methane emissions is of main interest for potential contributors, due to the high and rapid mitigation potential. In particular, supporting diet change levers such as increasing grass use, is also highlighted by life cycle assessment approaches has having the potential to reduce the total global warming potential and biodiversity damage of milk production (
      • O'Brien D.
      • Shalloo L.
      • Patton J.
      • Buckley F.
      • Grainger C.
      • Wallace M.
      Evaluation of the effect of accounting method, IPCC v. LCA, on grass-based and confinement dairy systems' greenhouse gas emissions.
      ;
      • Guerci M.
      • Knudsen M.T.
      • Bava L.
      • Zucali M.
      • Schönbach P.
      • Kristensen T.
      Parameters affecting the environmental impact of a range of dairy farming systems in Denmark, Germany and Italy.
      ). Developing PES schemes for the reduction of enteric methane emissions raises the question of the choice of a practical emission indicator sensitive to both diet and productivity effects, easily applicable and measurable on farm. While numerous other motivations may encourage farmers to join a PES program, such as improving milk quality, environmental quality, zootechnical performance and the image of agriculture, economic interests are likely to be crucial factors. An efficient payment level of a PES scheme targeting GHG emissions is equal to the socially optimal carbon price. In the EU, carbon tax levels differ largely among countries, illustrating the difficulty in pricing GHG emissions and therefore the social value of mitigation efforts. In 2021, it was set at 120€/teqCO2 in Sweden and 45€/teqCO2 in France, while the EU ETS (European Trading System) market price was 44€/teqCO2 (
      • The World Bank
      Carbon Pricing Dashboard | Up-to-Date Overview of Carbon Pricing Initiatives.
      ).
      In this paper, we aim to provide insights into the design of PES schemes targeting the reduction of enteric methane emissions in dairy farms by examining 2 aspects of a payment mechanism for which failing to consider the feeding dimension could undermine its effectiveness in cutting enteric methane emissions: the choice of emissions indicator defining the environmental service, and the level of payment. By comparing an indicator constructed using a methodology applied in a real-life PES case study from France to the Tier 2 indicator currently used in the French annual GHG emissions inventory which considers productivity only, we examine how the diet of dairy cows influences the enteric methane emissions attributed to farms. Changing the diet of cows to improve the milk fatty acid profile can generate additional production costs that are not yet evaluated. The second contribution is to quantify the additional cost of a change in the diet of cows at the farm level to evaluate the economic incentives needed for improving dairy systems toward more environment-friendly practices. We estimate the variable cost function of dairy farms at the scale of France and for different fodder systems.

      BACKGROUND ON ENTERIC METHANE EMISSION INDICATORS

      Many indicators of enteric methane emissions have been developed and adapted to specific constraints, often related to the scale of their application and the data available for estimating emissions (
      • Kebreab E.
      • Clark K.
      • Wagner-Riddle C.
      • France J.
      Methane and nitrous oxide emissions from Canadian animal agriculture: A review.
      ;
      • Ellis J.L.
      • Kebreab E.
      • Odongo N.E.
      • McBride B.W.
      • Okine E.K.
      • France J.
      Prediction of methane production from dairy and beef cattle.
      ;
      • Negussie E.
      • de Haas Y.
      • Dehareng F.
      • Dewhurst R.J.
      • Dijkstra J.
      • Gengler N.
      • Morgavi D.P.
      • Soyeurt H.
      • van Gastelen S.
      • Yan T.
      • Biscarini F.
      Invited review: Large-scale indirect measurements for enteric methane emissions in dairy cattle: A review of proxies and their potential for use in management and breeding decisions.
      ). For result-based PES, emissions should regularly and easily be monitored on farm. The emission indicator chosen must therefore be easily measurable by the farmers and/or the paying agent using a simple methodology, representative of the environmental target and reliable (based on strong and reliable scientific evidence) (
      • Allen B.
      • Hart K.
      • Radley G.
      • Tucker G.
      • Keenleyside C.
      • Oppermann R.
      • Underwood E.
      • Menadue H.
      • Poux X.
      • Beaufoy G.
      • Herzon I.
      • Povellato A.
      • Vanni F.
      • Prazan J.
      • Hudson T.
      • Yellachich N.
      Biodiversity protection through results based remuneration of ecological achievement. Report Prepared for the European Commission, DG Environment, Contract No ENV.B.2/ETU/2013/0046.
      ). In addition, it must be sensitive to the different dimensions on which farmers can act to decrease emissions, in particular cows' diet and productivity.
      The Intergovernmental Panel on Climate Change (IPCC) defines 3 families of methods, to be applied for national inventories according to data availability. Tier 1 methods attribute default yearly enteric methane emissions factor per dairy cow. Tier 1 methods provide aggregate estimates, and are not adequate for monitoring changes over time and taking into account for the variability of dairy farming practices. Tier 2 methods improve the accuracy of emission factors by including feed intake estimates of a representative diet and dairy cow (
      • Gavrilova O.
      • Leip A.
      • Dong H.
      • Macdonald J.D.
      • Gomez Bravo C.A.
      • Amon B.
      • Barahona Rosales R.
      • Agustin del Prado A.
      • Aparecida de Lima M.
      • Oyhantcabal W.
      • van der Weerden T.J.
      • Widiawati Y.
      Emissions from Livestock and Manure Management.
      ). Finally, Tier 3 methods require a precise characterization of cows' diet to account for digestibility. Both Tier 2 and Tier 3 approaches for France are presented in
      • Eugène M.
      • Sauvant D.
      • Nozière P.
      • Viallard D.
      • Oueslati K.
      • Lherm M.
      • Mathias E.
      • Doreau M.
      A new Tier 3 method to calculate methane emission inventory for ruminants.
      .
      Several recent studies applied those indicators.
      • Stetter C.
      • Sauer J.
      Greenhouse Gas Emissions and Eco-Performance at Farm Level: A Parametric Approach.
      applied a Tier 1 enteric methane emission indicator at the micro scale, but as part of an overall assessment of the relative GHG emissions and economic performance of farms. Life cycle assessment analyzes calculating the carbon footprint of dairy products tend to apply Tier 2 methods (
      • Jayasundara S.
      • Worden D.
      • Weersink A.
      • Wright T.
      • VanderZaag A.
      • Gordon R.
      • Wagner-Riddle C.
      Improving farm profitability also reduces the carbon footprint of milk production in intensive dairy production systems.
      ;
      • Wilkes A.
      • Wassie S.
      • Odhong' C.
      • Fraval S.
      • van Dijk S.
      Variation in the carbon footprint of milk production on smallholder dairy farms in central Kenya.
      ), sometimes also including diet composition data (
      • Hagemann M.
      • Hemme T.
      • Ndambi A.
      • Alqaisi O.
      • Sultana M.N.
      Benchmarking of greenhouse gas emissions of bovine milk production systems for 38 countries.
      ;
      • Gollnow S.
      • Lundie S.
      • Moore A.D.
      • McLaren J.
      • van Buuren N.
      • Stahle P.
      • Christie K.
      • Thylmann D.
      • Rehl T.
      Carbon footprint of milk production from dairy cows in Australia.
      ). Bioeconomic models estimating emissions abatement costs in the agricultural sector use Tier 3 indicators of enteric methane emissions precise enough to capture both productivity and diet effects (
      • Lengers B.
      • Britz W.
      • Holm-Müller K.
      Comparison of GHG-Emission indicators for dairy farms with respect to induced abatement costs, accuracy, and feasibility.
      ;
      • Mosnier C.
      • Britz W.
      • Julliere T.
      • De Cara S.
      • Jayet P.A.
      • Havlík P.
      • Frank S.
      • Mosnier A.
      Greenhouse gas abatement strategies and costs in French dairy production.
      ). Most Tier 3 and “individualized” Tier 2 approaches mentioned above require a large amount of detailed individual farm-level data on cows' feed composition, making it too costly to be applied as monitoring indicators on a large amount of farms. However, some of the Tier 3 methods can use information from milk analyzes and thus, be more easily integrated in farm routines.

      A TIER 3 INDICATOR APPLIED IN THE FRENCH ECO-METHANE RESULT-BASED SCHEME

      Numerous studies have been carried out to understand the connection between dairy cows' enteric methane emissions, fat intake and milk composition (
      • Dong Y.
      • Bae H.D.
      • McAllister T.A.
      • Mathison G.W.
      • Cheng K.J.
      Lipid induced depression of methane production and digestibility in the artificial rumen system (RUSITEC).
      ;
      • Martin C.
      • Morgavi D.
      • Doreau M.
      • Jouany J.P.
      Comment réduire la production de méthane chez les ruminants? How can the production of methane by ruminants be reduced?.
      ,
      • Martin C.
      • Rouel J.
      • Jouany J.P.
      • Doreau M.
      • Chilliard Y.
      Methane output and diet digestibility in response to feeding dairy cows crude linseed, extruded linseed, or linseed oil.
      ,
      • Martin C.
      • Morgavi D.P.
      • Doreau M.
      Methane mitigation in ruminants: from microbe to the farm scale.
      ,
      • Martin C.
      • Promiès D.
      • Ferlay A.
      • Rochette Y.
      • Martin B.
      • Chilliard Y.
      • Morgavi D.
      • Doreau M.
      Methane output and rumen microbiota in dairy cows in response to long-term supplementation with linseed or rapeseed of grass silage or pasture based diets.
      ;
      • Chilliard Y.
      • Martin C.
      • Rouel J.
      • Doreau M.
      Milk fatty acids in dairy cows fed whole crude linseed, extruded linseed, or linseed oil, and their relationship with methane output.
      ;
      • Grainger C.
      • Beauchemin K.A.
      Can enteric methane emissions from ruminants be lowered without lowering their production?.
      ). In particular, experimental research shows that enteric methane emissions in gCH4/L (Methane) can be calculated from milk productivity in kg/cow/year (Productivity) and the ratio of the sum of fatty acids with 16 carbon atoms or less (FA ≤ C16) to the total amount of fatty acids (TotalFA) in the milk they produce (1). This ratio has a strong biological causal relationship with methanogenesis in the rumen and is significantly reduced by more sources of omega-3 in cows' diet.
      Methane=11.368×productivity0.4274×FA16totalFA
      [1]


      This formula was coinvented by teams from the animal feed manufacturing company Valorex (P. Weill and G. Chesneau) and the French National Institute for Agricultural Research (INRA) (Y. Chilliard, M. Doreau and C. Martin), and received a patent under the title “Method for evaluating the quantity of methane produced by a dairy ruminant and method for decreasing and controlling such quantity” (WO2009156453A1) (

      Weill, P., G. Chesneau, Y. Chilliard, M. Doreau, and C. Martin. 2009. Method to evaluate the quantity of methane produced by a dairy ruminant and method for decreasing and controlling this quantity. Valorex, assignee. Pat. No. WO 2009/156453 A1.

      ). The necessary data to implement it are relatively easy to collect on farm. One needs milk productivity and an analysis of the fatty acid composition of milk, obtained with infrared spectroscopy. Milk infrared spectroscopy is relatively simple to integrate into the milk analyzes routines of dairy farms and involves low costs. Since dairy farms already undertake milk analyzes on a regular basis, it makes it a relatively cheap indicator to monitor emissions. In addition, the equation calculates enteric methane emissions per unit of product by taking into account both milk productivity and feed quality, hence captures both dimensions on which farmers can dedicate effort to reduce emissions.
      The indicator described above is used since 2011 to evaluate the reduction of enteric methane emissions of dairy farms participating in the PES scheme Eco-Methane in France, implemented and coordinated by the Bleu-Blanc-Coeur association (BBC). Eco-Methane meets the reference PES definition by
      • Wunder S.
      Revisiting the concept of payments for environmental services.
      . In the scheme, private stakeholders (service users) give financial support to volunteer dairy farmers (service providers) for actions that contribute to climate change mitigation (environmental services). The payment is conditional and proportional to the reduction of CO2eq, making the scheme a result-based PES. The reduction is calculated relative to baseline monthly emission levels attributed to the farm that are representative of its type of fodder system. Through the definition of different baseline scenarios, the PES design partially considers the variability in the potentiality of environmental services provision according to the production basin and the fodder system. Hence, rather than rewarding farms that produce the least emissions per unit of product (which would tend to favor the most productive farms), Eco-Methane supports all emission reduction efforts.
      The Eco-Methane scheme uses 11 scenarios of baseline emission levels representative of 11 fodder systems of French specialized dairy farms. These fodder systems were characterized by the French Dairy Interbranch Organization (CNIEL) in collaboration with the French Livestock Institute (IDELE) in 2009, based on large production basins (plain and mountainous areas), and the proportion of corn in fodder crop rotation systems (
      • CNIEL
      Observatoire de l'alimentation des vaches laitières: description des 8 principaux systèmes d'élevage. Edition 2015–2018. Dairy cows' feed observatory: description of the 8 main farming systems. 2015–2018 edition (in French).
      ). The baseline emissions for each scenario are determined by BBC using (1).
      Eco-Methane brings together more than 600 farmers whose emissions reduction was estimated at 11% on average in 2017 (
      • Bleu-Blanc-Coeur
      Démarche Environnementale : La Démarche Éco-Méthane de Bleu-Blanc-Cœur. Environmental Engagement: The Eco-Methane Programme of Bleu-Blanc-Coeur (in French).
      ). BBC pays farmers according to their reduction of methane emissions in CO2eq with a financial envelope made of donations from private actors (15€/tCO2eq on average in 2017). The main strengths of the scheme lie in the strong scientific foundations of the method for quantifying emissions and the easy participation procedure for dairy farmers. Each contract signatory commits to provide a monthly milk analysis to the association and to include feed with a high content of sources of omega-3 in dairy cow rations (alfalfa, extruded linseed, grass). This data is used to estimate enteric methane emissions using (1). Eco-Methane is recognized by the United Nations as a GHG emission reduction project eligible for issuing carbon credits (
      • UNFCCC
      Joint Implementation Project FR1000365: “Réduction des émissions de méthane d'origine digestive par l'apport dans l'alimentation des vaches laitières de sources naturelles en Acide Alpha Linolénique (ALA)” (in French).
      ).

      MATERIALS AND METHODS

      Data

      Observations of a balanced panel of 735 French farms specialized in dairy milk production (OTEXE 45) from France's farm accountancy data network (FADN) for the years 2016 to 2018 was selected for the analysis. This European database is freely accessible online and representative of socioeconomic and accountancy information of medium and large dairy farms contributing to more than to more than 90% of the gross production and utilized agricultural area. Due to this characteristic, it is a particularly relevant data set to investigate emission abatement costs at the national level. As the compositions of the feed ration and milk are not surveyed, information on the diet of dairy cows is limited and in particular the fatty acid profile is unknown. Instead, data on the fodder crop rotation systems are used to assess a change in crop rotation and approximate a change in feed composition. Descriptive statistics of the sample are presented (Table 1).
      Table 1Description of the sample (2205 observations)
      Variable1st quartileMedianMean3rd quartile
      Utilized Agricultural Area (ha)
      Information on surfaces and the number of dairy cows available in the database are ranges of values. We constructed the variables used in the analysis by taking the lower value of the range for each observation. Source: The authors, based on French FADN data.
      50.080.087.4110.0
      Fodder area (ha)
      Information on surfaces and the number of dairy cows available in the database are ranges of values. We constructed the variables used in the analysis by taking the lower value of the range for each observation. Source: The authors, based on French FADN data.
      40.060.067.480.0
      Corn silage area (ha)
      Information on surfaces and the number of dairy cows available in the database are ranges of values. We constructed the variables used in the analysis by taking the lower value of the range for each observation. Source: The authors, based on French FADN data.
      1.010.014.120.0
      Pasture area (permanent and temporary) (ha)
      Information on surfaces and the number of dairy cows available in the database are ranges of values. We constructed the variables used in the analysis by taking the lower value of the range for each observation. Source: The authors, based on French FADN data.
      26.040.050.361.0
      Productivity (L/cow)5,593.46,676.46,707.97,851.1
      Number of dairy cows
      Information on surfaces and the number of dairy cows available in the database are ranges of values. We constructed the variables used in the analysis by taking the lower value of the range for each observation. Source: The authors, based on French FADN data.
      35555870
      Agricultural Work Unit1.02.01.82.1
      Purchase of cattle feed concentrates (€)14,326.024,996.532,853.243,645.0
      1 Information on surfaces and the number of dairy cows available in the database are ranges of values. We constructed the variables used in the analysis by taking the lower value of the range for each observation.Source: The authors, based on French FADN data.

      Attribution of enteric methane emissions

      We use the Tier 2 method used in the French annual inventory of GHG emissions to define an indicator of enteric methane emissions sensitive to productivity effects. Dairy cows emission factors are calculated from equation (2) (
      • Citepa
      Organisation et méthodes des inventaires nationaux des émissions atmosphériques en France OMINEA - 17 ème édition. Inventory methodology report of GHG and air pollutant emissions in France - 17th edition (in French).
      ). The parameters of the formula are such that they are representative of breeding, feeding and genetic conditions in France. Based on this calibration, the emission factor per cow varies according to milk productivity.
      EF=0.0105×MilkproductionNcows+48.971
      (2)


      The emission factor EF (kgCH4/cow/year) can be easily calculated for each farm of the FADN from the milk production (kg/year) of the herd (Milk production) and the number of dairy cows (Ncows). We then derive an emission indicator (TIER2) per liter of milk (gCH4/L), capturing variability according to milk productivity (L/cow/year) (3).
      TIER2=EF×1,000productivity
      (3)


      Due to the absence of data on dairy cows' diet, a Tier 3 method cannot be applied to evaluate individual emissions of FADN farms and capture both productivity and diet effects. In particular, data are too limited to estimate individual enteric methane emissions of French farms using (1). They are, however, sufficient to identify their Eco-Methane scenario and therefore the baseline emissions corresponding to their fodder system. Baseline emissions from the 11 scenarios are available per month and were obtained from BBC. We calculate the annual average to define the Eco-Methane baseline emissions indicator (Eco-Methane baseline) (Table 2).
      Table 2Characteristics of the eleven baseline scenarios used in the Eco-Methane scheme
      ScenarioCorn in the fodder areaProduction basinEco-Methane baseline (gCH4/L)
      1More than 30%Plains outside the western region15.75
      2Plains of the western region15.92
      3Between 10 and 30%Plains outside the western region15.83
      4Plains of the western region16.43
      5Less than 10%Plains outside the western region16.56
      6Plains of the western region17.38
      7More than 10%Mountains15.96
      8Less than 10%Mountains of the Massif Central17.13
      9Mountains of the Northern Alps17.83
      10Mountains of Franche-Comté16.22
      11Other mountains17.20
      Source: The authors, based on BBC data.
      Source: The authors, based on French FADN and BBC data.
      An individual baseline scenario was assigned to each farm of the sample based on 2 criteria: the location and the share of corn silage in the fodder area of the farm. In the FADN database, the farm location variable corresponds to the 21 old French administrative regions (the administrative divisions were changed in 2015), while the Eco-Methane scenarios are defined according to large production basins built from a lower administrative level (departments). It was therefore necessary to allocate a production basin to each administrative region. For the regions with departments belonging to different production basins, we allocated the basin of the departments producing the highest volumes of milk to the entire region. This attribution was made using the 2018 annual dairy survey (

      Agreste. 2019. Enquête annuelle laitière 2018. French 2018 annual dairy survey (in French). Agreste Chiffres et Données 13.

      ) and is detailed in the Appendix (Table A1) .

      Estimation of the extra-cost of milk production with more grass in the fodder system

      In France, the closest financial tool to a carbon tax is the Climate and Energy Contribution proportional to the carbon dioxide content of energy products (fossil fuels) (
      • Rogissart L.
      • Postic S.
      • Grimault J.
      La Contribution Climat Energie en France : fonctionnement, revenus et exonérations. The Climate-Energy Contribution in France: operation, revenues and exemptions.
      ). The contribution level was 30€/tCO2eq in 2017 and increased to 40€/tCO2eq in 2018 and 2019. Farmers participating in Eco-Methane received an average of 15€/tCO2eq in 2017, suggesting that the payment of the scheme is suboptimal and provides little incentive to participate (
      • Bleu-Blanc-Coeur
      Démarche Environnementale : La Démarche Éco-Méthane de Bleu-Blanc-Cœur. Environmental Engagement: The Eco-Methane Programme of Bleu-Blanc-Coeur (in French).
      ). Evaluating farmers' willingness to accept is necessary to define a more efficient price and trigger a large scale adoption of practices decreasing enteric emissions.
      Since the composition of the feed ration and milk of cows are not available in the FADN, the effect of an improvement of the fatty acid profile on milk production costs cannot be analyzed directly. Instead, an evolution of the fodder crop rotation is assumed. As grass is a high source of omega-3 fatty acids strongly encouraged in Eco-Methane, we assume that a commitment to the program would lead to an increase in grassland surfaces in farms. This hypothesis is quite strong and implies that the estimation of extra costs does not consider either the strategy of supplementing the ration with other feeds with high omega-3 content such as extruded linseed or the optimization of grazing increasing grass yield and quality without necessarily increasing grassland surfaces.
      Based on dual production theory (
      • McFadden D.
      Cost, Revenue and Profit Functions.
      ), we estimate a variable cost function describing expenditures in variable production factors x with exogenous input prices w that minimize variable costs given the production level y targeted by the farmer and available quasifixed inputs z such as land, labor and equipment that are assumed to be predetermined in the short term. The cost minimization approach is motivated by the fact that French dairy farmers are constrained in the quantity of milk they produce in the terms of their contract with dairy plants (
      • Lambaré P.
      • Dervillé M.
      • You G.
      What will be the conditions of market access for dairy farmers after the end of dairy quotas?.
      ). It is confirmed in our data (see Table A2 in the Appendix), as we observe a low variability of milk volumes across years for a same farm. The variable cost function (4) meets theoretical properties that we empirically check.
      VC(w,y,z)=minxwxsubjecttoyf(x,z)
      (4)


      It must be concave, nondecreasing and homogeneous of degree 1 in variable input prices, decreasing with (binding) fixed factors of production, and monotonic according to output levels.
      We estimate a system of equations comprising a homogeneous translog cost function (5) in which variable costs VC correspond to intermediate consumption and the variable input cost share functions (6) and (7), derived from Shephard's lemma. The translog functional form is commonly used in the literature on the cost structure and efficiency of dairy farms because of its flexibility and the possibility of imposing homogeneity of degree 1 (
      • Moschini G.
      The Cost Structure of Ontario Dairy Farms: A Microeconometric Analysis.
      ;
      • Alvarez A.
      • Arias C.
      Diseconomies of size with fixed managerial ability.
      ;
      • Mosheim R.
      • Lovell C.A.K.
      Scale economies and inefficiency of U.S. dairy farms.
      ;
      • Nehring R.
      • Gillespie J.
      • Sandretto C.
      • Hallahan C.
      Small U.S. dairy farms: can they compete?.
      ;
      • Sobczyński T.
      • Klepacka A.M.
      • Revoredo-Giha C.
      • Florkowski W.J.
      Dairy farm cost efficiency in leading milk-producing regions in Poland.
      ;
      • Tsionas E.G.
      • Kumbhakar S.C.
      • Malikov E.
      Estimation of Input Distance Functions : A System Approach.
      ;
      • Singbo A.
      • Larue B.
      Scale economies, technical efficiency, and the sources of total factor productivity growth of Quebec dairy farms.
      ;
      • Wimmer S.
      • Sauer J.
      Diversification economies in dairy farming - Empirical evidence from Germany.
      ). i and t are indices for individuals and years, respectively. We assume that dairy farms produce one output, the quantity Y1 of the milk of cows produced per year (milk production). Considering the joint production of milk and a bad output (enteric methane emissions) would be of main interest. This is the approach used in (
      • Njuki E.
      • Bravo-Ureta B.E.
      The economic costs of environmental regulation in U.S. Dairy farming: a directional distance function approach.
      ;
      • Le S.
      • Jeffrey S.
      • An H.
      Greenhouse Gas Emissions and Technical Efficiency in Alberta Dairy Production : What Are the Trade-Offs?.
      ) for instance. Implementing it in our study would require having enough information to compute an individualised enteric methane emission indicator capturing both cows' productivity and diet. The FADN data set does not provide enough information on cows' diet to compute the Eco-Methane indicator for each individual farm and we cannot precisely compute the “bad” output.
      We consider 2 variable inputs, fuel X1 and cattle feeding stuffs X2, for which the expenses represent a high share of intermediate consumption. The choice of including fuel is motivated by the possibility of calculating farm-level fuel prices and therefore capturing more heterogeneity. The price of fuel W1 (fuel price) is calculated from the nonroad gas oil expenses and volumes. As individual cattle feeding stuff prices are not available in the data, W2 (feed price) is measured by the index of purchase prices of the means of agricultural production (IPAMPA) for adult cattle feeding stuffs of year t-1, available for each French current administrative region. While in practice, a bargaining effect may create endogeneity between milk production and feed prices, we assume it is sufficiently low not to overestimate too much milk marginal production costs (
      • Duvaleix-Tréguer S.
      • Gaigné C.
      On the nature and magnitude of cost economies in hog production.
      ). Marginal costs tend to be overestimated when the adjustment of the feed unit prices to a change of farm size is not controlled for. In our case, capturing this endogeneity issue would require individual farm feed price data, information we do not have.
      • Duvaleix-Tréguer S.
      • Gaigné C.
      On the nature and magnitude of cost economies in hog production.
      found evidence of a price bargaining effect for hog production, but lower than the technology effect. For hog production, feed represents a larger share of the costs (about 60% of total cost) than for dairy production (25% of variable cost, see Table A2 of the appendix). Therefore, the price bargaining effect is expected to be lower for dairy farms.
      We also include 3 quasifixed inputs. Grassland surface Z1 (grassland) includes permanent and temporary pastures, alfalfa for dehydration and other artificial fodders. We add 2 other quasifixed factors of production: machinery and constructions fixed assets Z2 (capital) and annual work units Z3 (labor). Since sample farms are all specialized in dairy milk production, we do not consider them as multi-output farms. We verify the production of crops and other livestock products represents a much lower share of total gross products (see Table A2 in the Appendix). Yet, we add the aggregated volume Y2 of the other products of the farm (other productions) as a control variable to capture heterogeneity linked to diversification. Y2 is calculated as the total gross product of the year (crop products, livestock products and other products) net of animal purchases and milk production of the cows, deflated by the French agricultural producer price index (API) of year t. We further control for the farm total utilized agricultural area Q1 (utilized agricultural area), 20 regional dummies Dri approximating the pedoclimatic conditions, and time fixed effects with 2 yearly dummies Tt. The coefficients α, β, δ, ν, ρ, ζ, θ, γ, and µ are the parameters to be estimated and describe the effects of the covariates on the independent variables. eit, u1it and u2it are the error terms assumed to follow a normal distribution.
      lnVCitW1it=α0+β1lnY1it+α2lnW2itW1it+12α22(lnW2itW1it)2+h=13δhlnZhit+12β11lnY1it2+12h=13k=13δhklnZhitlnZkit+h=13ν2hlnW2itW1itlnZhith=13ρ1hlnY1itlnZhit+ζ12lnY1itlnW2itW1it+β2Y2it+θ1Q1i+r=120γrDri+t=1617μtTt+et,
      (5)


      X1itW1itVCit=1α2α22lnW2itW1ith=13ν2hlnZhitζ12lnY1it+u1it,
      (6)


      X2itW2itVCit=α2+α22lnW2itW1it+h=13ν2hlnZhit+ζ12lnY1it+u2it.
      (7)


      To account for the correlations between the error terms of the different equations, we use a 3-stage least squares regression analysis. The system of equations (5) + (6) + (7) is estimated for all the farms of the sample and then for the 3 major production basins and groups of Eco-Methane scenarios defined in Table 2 to identify potential differences in extra costs according to the type of dairy system. The different steps to derive the specification of the system are described in the Appendix. Descriptive statistics of the model variables of the sample and subsamples are presented in Table 3.
      Table 3Mean values of the model variables (standard deviation in parentheses)
      VariableFrance n = 2,205Western plains n = 645Plains outside the western region n = 965Mountains n = 585Plains, > 30% of corn in the fodder area n = 767Plains, 10-30% of corn in the fodder area n = 574Plains, <10% of corn in the fodder area n = 279
      Scenarios1–112, 4, 61, 3, 57–111–23–45–6
      Variable Costs (€/year)128,073.5139,405.3136,399.496,533.0170,654.1124,656.090,248.5
      (105,612.7)(107,167.1)(112,800.6)(80,772.0)(124,937.8)(71,877.3)(57,036.7)
      Milk production (L/year)398,594.1446,128.0402,355.2306,703.0526,707.0396,281.9156,204.9
      (297,513.2)(328,781.5)(296,923.0)(226,416.1)(333,580.8)(213,270.9)(266,216.2)
      Other productions (€/base 100/year)498.6581.3554.6274.4750.2488.6317.0
      (671.5)(642.3)(770.2)(413.0)(871.1)(432.8)(349.0)
      Fuel price (€/L)0.600.590.590.610.590.590.60
      (0.10)(0.10)(0.11)(0.10)(0.10)(0.10)(0.10)
      Feed price (base 100)96.696.596.696.896.696.596.6
      (2.4)(2.4)(2.4)(2.3)(2.4)(2.4)(2.4)
      Grassland (ha)51.142.152.765.534.250.666.1
      (41.0)(31.5)(38.6)(47.8)(25.8)(32.5)(36.1)
      Capital (1000€)171.0160.5179.4179.5189.1144.8163.4
      (155.9)(146.8)(163.7)(160.8)(161.8)(141.1)(153.5)
      Labor (AWU)1.81.91.91.72.01.71.7
      (1.0)(1.0)(1.0)(1.0)(1.1)(0.8)(0.9)
      UAA (ha)87.484.392.686.393.385.579.3
      (58.1)(55.3)(60.0)(57.8)(60.6)(48.2)(42.1)
      Source: The authors, based on French FADN data.

      RESULTS AND DISCUSSION

      First, we present the allocation of the FADN dairy farms in the Eco methane scenarios, then we compare the 2 indicators Tier 2 and the Eco-Methane one. Finally, we discussed the extra-costs associated with an increase in farm grassland area based on estimated marginal costs.

      Allocation of Eco-Methane scenarios and distinction of 3 milk production basins

      Following the allocation of the Eco-Methane baseline scenarios to the farms of the sample, their proportion within each former administrative region can be observed (Figure 1). Farms in the regions of the western plains basin (Brittany, Pays de la Loire and Basse-Normandie) represent 44% of the sample and are characterized by a strong dominance of corn silage and few grasslands in forage crop rotation systems. For example, 68% of farms in Brittany would be assigned scenario 2 with more than 30% of corn in the forage crop rotation and 28% scenario 4 with 10 to 30%. The administrative regions of the plain production basin outside the western region (31% of the sample) and of the mountainous areas (24% of the sample) are quite different. Some plain regions, such as Rhône-Alpes, contain a high proportion of grazing systems dominated by grasslands (less than 10% of corn silage in the fodder crop rotation), while others, such as the Centre, have more intensive systems dominated by corn silage. All the observations from Languedoc-Roussillon correspond to grazing systems, while the observations from Midi-Pyrénées have a relatively small proportion (32%). Due to the missing information on departments in the FADN data set and our scenario allocation procedure, some farms have been allocated to a plain system scenario, while, in reality, they are located in a mountainous department and vice versa. This allocation might partly explain the large share of farms with grazing systems in Rhône-Alpes. Nevertheless, those farms produce relatively low volumes of milk in comparison with plain farms from the same region.
      Figure thumbnail gr1
      Figure 1Distribution of Eco-Methane baseline scenarios among French old administrative regions.

      Enteric emissions: relation to productivity and fodder system

      The mean of the Tier 2 indicator for the sample is 18.5 gCH4/L, while it is 16.3 gCH4/L for the Eco-Methane baseline indicator, suggesting that considering the fodder cropping system in the calculation revises enteric emissions downwards. Both indicators show a decrease in emissions per liter of milk as productivity increases (Table 4). Dairy systems with lower enteric emissions per liter of milk are productive systems with more than 30% of corn silage in the fodder crop rotation, while grass-based systems have the highest emissions. The comparison of the 2 indicators allows to nuance the performance of the different systems in terms of diet and productivity effects.
      Table 4Average enteric emissions of the sample according to the two indicators
      ScenariosPlains outside the western regionWestern plainsMountains
      1, 3, 51352, 4, 62467–1178 to 11
      % corn in fodder area>3030–10<10>3030–10<10≥10<10
      Sample share (%)31.410.09.911.544.323.716.83.824.36.517.8
      Productivity (1,000L/cow)6.727.656.945.726.987.336.795.596.206.915.94
      TIER2 (gCH4/L)18.55d17.46l18.22fhj19.78a18.23ehi17.80k18.38defg20.22a19.07c18.19gij19.40b
      Eco-Methane baseline (gCH4/L)16.07q15.75u15.83t16.56n16.24p15.92s16.43°17.38l16.55n15.96r16.76m
      Mean values with different superscripts differ (p-value <0.05).
      Source: The authors, based on French FADN and BBC data.
      According to indicator Tier 2, farms in mountains emit significantly more methane per liter of milk on average than farms in the plains which can be explained by their lower productivity. The same observation is made with the Eco-Methane indicator, but the difference between the groups is less (Table 4 and Figure 2). This finding shows that a diet effect partly compensates high emissions from low productivity in mountainous areas. The average share of grass in mountain fodder crop rotation systems is 90%, providing a diet rich in sources of omega-3 to dairy cows and contributing to the mitigation of enteric emissions. In fact, we observe that the difference between the 2 indicators is large for low productivity and grass-based dairy systems, and low for productive systems with a large share of corn silage in the fodder area (see Figure A1 in the Appendix).
      Figure thumbnail gr2
      Figure 2Average enteric emissions in French old administrative regions according to the 2 indicators Source: The authors, based on French FADN and BBC data.
      In plain areas, the average methane emissions per liter of milk of dairy farms with a high share of corn silage in the fodder area are higher in western regions because of lower productivity (Table 4). However, the difference between the 2 plain production basins is not significantly different for a fodder system with an intermediate share (scenario 3 vs 4), or a low share (scenario 5 vs 6) of corn silage. Therefore, computing enteric emissions with the Tier 2 indicator results in no significant difference between the 2 plain production basins for a same fodder system. However, the Eco-Methane baseline attributed to plain farms outside the western region is lower. In western regions, farms have a significantly lower share of grass in the fodder rotation system (P-value < 0.01). It suggests that cows receive less omega-3 rich fodder in their ration.
      Our observations are in line with the literature (
      • Martin C.
      • Morgavi D.
      • Doreau M.
      • Jouany J.P.
      Comment réduire la production de méthane chez les ruminants? How can the production of methane by ruminants be reduced?.
      ;
      • Grainger C.
      • Beauchemin K.A.
      Can enteric methane emissions from ruminants be lowered without lowering their production?.
      ). A recent meta-analysis of life cycle assessments also highlighted the negative relationships between milk yield and enteric methane emissions on the one hand and pasture intake and enteric methane emissions on the other hand (
      • Lorenz H.
      • Reinsch T.
      • Hess S.
      • Taube F.
      Is low-input dairy farming more climate friendly? A meta-analysis of the carbon footprints of different production systems.
      ). Other authors also show that enteric methane emissions particularly differ from one indicator to another in grazing systems (
      • Hagemann M.
      • Hemme T.
      • Ndambi A.
      • Alqaisi O.
      • Sultana M.N.
      Benchmarking of greenhouse gas emissions of bovine milk production systems for 38 countries.
      ). Choosing an adequate indicator of enteric methane emissions is the topic of ongoing debate. This analysis illustrates the importance of the choice of environmental indicator when designing a PES scheme targeting the reduction of enteric methane emissions in dairy farms, as it is likely to affect its environmental performance. We show the relevance of using an emission indicator sensitive to diet when defining economic incentives for the reduction of GHG emissions.
      Detailed indicators are the most cost-effective, but due to heavy data collection needs (precise feed digestibility and composition), this advantage decreases when applied at a large scale (
      • Lengers B.
      • Britz W.
      • Holm-Müller K.
      Comparison of GHG-Emission indicators for dairy farms with respect to induced abatement costs, accuracy, and feasibility.
      ). An indicator such as the one used in Eco-Methane presents several advantages to be implemented on a large scale, and be a better proxy compared with the most widely used Tier 2 and data demanding Tier 3 indicators. It is precise enough to capture the efforts of farmers on both cow productivity and feed ration composition, and it considers the potential of an omega-3-rich diet as a climate change mitigation practice. This feeding strategy is already implemented in dairy systems integrating a large share of grasslands in their fodder crop rotation systems, with the side provision of other environmental benefits (biodiversity maintenance). Additionally, the fatty acid composition of milk provides information on the complementary health benefits for consumers of an increase in sources of omega-3 fatty acids in the diet of dairy cows (
      • Weill P.
      • Schmitt B.
      • Chesneau G.
      • Daniel N.
      • Safraou F.
      • Legrand P.
      Effects of introducing linseed in livestock diet on blood fatty acid composition of consumers of animal products.
      ). The accuracy of indicators based on milk analyzes could be further improved by controlling for factors likely to affect the correlation between milk fatty acid composition and enteric emissions such as the lactation stage (
      • Negussie E.
      • de Haas Y.
      • Dehareng F.
      • Dewhurst R.J.
      • Dijkstra J.
      • Gengler N.
      • Morgavi D.P.
      • Soyeurt H.
      • van Gastelen S.
      • Yan T.
      • Biscarini F.
      Invited review: Large-scale indirect measurements for enteric methane emissions in dairy cattle: A review of proxies and their potential for use in management and breeding decisions.
      ).
      Furthermore, the environmental performance of PES schemes specifically targeting enteric methane emissions depends on the absence of negative spillovers on other factors of GHG emissions in dairy farms (fertilization management, machinery…). On dairy farms, other important sources of GHG emissions are energy use, manure production and management, fertilizer production and application, and concentrate feed production (
      • Chen W.
      • Holden N.M.
      Tiered life cycle sustainability assessment applied to a grazing dairy farm.
      ;
      • Jayasundara S.
      • Worden D.
      • Weersink A.
      • Wright T.
      • VanderZaag A.
      • Gordon R.
      • Wagner-Riddle C.
      Improving farm profitability also reduces the carbon footprint of milk production in intensive dairy production systems.
      ). Using more grass is a lever that can also contribute to lower emissions from concentrate feed and to increase carbon storage in agricultural soils (
      • Mosnier C.
      • Britz W.
      • Julliere T.
      • De Cara S.
      • Jayet P.A.
      • Havlík P.
      • Frank S.
      • Mosnier A.
      Greenhouse gas abatement strategies and costs in French dairy production.
      ;
      • Wilkes A.
      • Wassie S.
      • Odhong' C.
      • Fraval S.
      • van Dijk S.
      Variation in the carbon footprint of milk production on smallholder dairy farms in central Kenya.
      ). Moreover, dairy farms with more grasslands also tend to emit less GHG emissions in total at the farm level (
      • Senga Kiessé T.
      • Corson M.S.
      • Wilfart A.
      Analysis of milk production and greenhouse gas emissions as a function of extreme variations in forage production among French dairy farms.
      ). However, nitrous oxide emissions from soil can increase from more manure deposition on pasture (
      • Del Prado A.
      • Mas K.
      • Pardo G.
      • Gallejones P.
      Modelling the interactions between C and N farm balances and GHG emissions from confinement dairy farms in northern Spain.
      ;
      • Wilkes A.
      • Wassie S.
      • Odhong' C.
      • Fraval S.
      • van Dijk S.
      Variation in the carbon footprint of milk production on smallholder dairy farms in central Kenya.
      ). Decreasing enteric methane emissions per liter of milk is a lever particularly interesting to rapidly reduce the contribution of dairy farming to climate change while also considering the food security dimension as it refers to the quantity of food produced. But it should be complementary to further farm-level and area-based assessments of GHG emissions to support effective mitigation strategies.

      Impact of increasing the grassland area on the marginal cost of milk production

      Several model specifications were tested, and the results are robust to a change (single variable costs equation, system of equations with and without imposing constraints on the parameters for ensuring homogeneity of degree 1). The estimation of variable input shares provides additional information and improves the quality of the variable cost estimation (measured by R2). Consistent with the hypothesis of cost minimization, imposing restrictions on the parameters across equations also improved the variable cost estimation quality. Therefore, we present the results of the constrained system estimation.
      The variable cost functions estimated are homogeneous of degree 1 in variable input prices by specification, and we verify they are concave and nondecreasing with input prices (positive estimated variable input cost shares), and monotonic in the milk production level. However, we observe that variable costs increase with the level of fixed factors of production for some observations (see Table A3 in the Appendix). In particular, variable costs are decreasing with more capital for a significant share of farms in 4 models out of 7. Previous studies also found evidence of the violation of theoretical properties on quasi-fixed inputs, in particular capital (
      • Mosheim R.
      • Lovell C.A.K.
      Scale economies and inefficiency of U.S. dairy farms.
      ;
      • Singbo A.
      • Larue B.
      Scale economies, technical efficiency, and the sources of total factor productivity growth of Quebec dairy farms.
      ;
      • Wimmer S.
      • Sauer J.
      Diversification economies in dairy farming - Empirical evidence from Germany.
      ).
      The first-order derivative of the variable cost function (5) gives the marginal cost function (8), in which parameter ρ11 corresponds to the effect of grassland surfaces on the marginal cost of milk.
      VCitY1it=VCitY1it(β1+β11lnY1it+h=13ρ1hlnZhit+ζ12lnW2itW1it)
      (8)


      The results presented in the following paragraphs are calculated from the regression results detailed in the Appendix (Tables A4-A10).
      When applied to all dairy farms of the sample and on the subsamples of the western plains production basin and mountainous areas, the model suggests that producing milk with more grass does not significantly affect variable costs (Table 5). There are low but significant extra costs per additional hectare of grassland in plain dairy farms outside the western regions (+0.6€/1000L/ha). This finding illustrates the natural competitive advantage that dairy farms from the western regions or mountainous areas have in producing more milk with grass, the former because of favorable temperate climate (
      • Hennessy D.
      • Delaby L.
      • van den Pol-van Dasselaar A.
      • Shalloo L.
      Increasing grazing in dairy cow milk production systems in Europe.
      ), the latter because of climatic and technical limitations for another land use (
      • Coppa M.
      • Chassaing C.
      • Sibra C.
      • Cornu A.
      • Verbič J.
      • Golecký J.
      • Engel E.
      • Ratel J.
      • Boudon A.
      • Ferlay A.
      • Martin B.
      Forage system is the key driver of mountain milk specificity.
      ). A first implication for the design of a PES is that dairy farms from plain areas with a natural disadvantage for producing milk with grass require higher levels of economic incentives in comparison with other production basins. Other authors found evidence of differences in emissions abatement costs among dairy farms according to their geographical location (
      • Njuki E.
      • Bravo-Ureta B.E.
      The economic costs of environmental regulation in U.S. Dairy farming: a directional distance function approach.
      ).
      Table 5Extra-costs of milk production with an increase in grassland area
      FrancePlains of the western regionPlains outside the western regionMountainsPlains with more than 30% of corn silage in the fodder areaPlains with 10 to 30% of corn silage in the fodder areaPlains with less than 10% of corn silage in the fodder area
      Marginal cost (€/1000L)241.5230.5251.0261.1225.7225.1262.0
      Extra-cost (€/1000L/ha)−0.02−0.480.58+−0.41−1.38+0.28−1.82***
      Variable cost regression R20.910.910.910.910.900.910.89
      +P < 0.10, *P < 0.05, ** P < 0.01, *** P < 0.001.
      Source: The authors, based on French FADN and BBC data.
      Behind the heterogeneous extra costs found in plain production basins also lies a disparity depending on the type of fodder system (Table 5). Considering plain dairy farms, we compare those with a share of corn silage in the fodder area greater than 30%, between 10 and 30% and less than 10%. Extra-costs are significantly negative when the share of corn is higher than 30% (−1.4€/1000L/ha) or lower than 10% (−1.8€/1000L/ha), while they are not significant for an intermediate share between 10 and 30%. This second result reveals the economic opportunity of producing milk with more grass in terms of higher feed self-sufficiency (lower cost feed) (
      • Hanrahan L.
      • McHugh N.
      • Hennessy T.
      • Moran B.
      • Kearney R.
      • Wallace M.
      • Shalloo L.
      Factors associated with profitability in pasture-based systems of milk production.
      ). For plain grass-based systems with low milk productivity, achieving high self-sufficiency compensates otherwise relatively high marginal costs of production. For “corn silage” plain systems, reducing dependency on high protein content complements represents an important source of cost savings. Current dominant intensive fodder systems involve high expenditures on specific corn inputs (seeds, herbicides, etc.) and high protein content complements (soya, rapeseed) to balance dairy cow feed rations. Synergies between the reduction of GHG emissions and the economic performance of intensive dairy farms have already been pointed out in the literature (
      • Borreani G.
      • Coppa M.
      • Revello-Chion A.
      • Comino L.
      • Giaccone D.
      • Ferlay A.
      • Tabacco E.
      Effect of different feeding strategies in intensive dairy farming systems on milk fatty acid profiles, and implications on feeding costs in Italy.
      ;
      • Jayasundara S.
      • Worden D.
      • Weersink A.
      • Wright T.
      • VanderZaag A.
      • Gordon R.
      • Wagner-Riddle C.
      Improving farm profitability also reduces the carbon footprint of milk production in intensive dairy production systems.
      ;
      • Le S.
      • Jeffrey S.
      • An H.
      Greenhouse Gas Emissions and Technical Efficiency in Alberta Dairy Production : What Are the Trade-Offs?.
      ). In contrast, no evidence of this economic opportunity is found for dairy systems with an intermediate share of corn silage in the fodder area. For those fodder systems, higher variable costs from an increase in energy consumption (machinery) and other expenses (seeds, fertilizers, etc.) related to additional pastures and alfalfa management may limit feed cost savings. A second implication of our findings is that currently, a PES for producing milk with more grass with a low level of payment favors adoption by plain fodder systems with a large share of corn silage on the one hand, and fodder systems with a large share of grasslands on the other hand, over intermediate systems combining both corn silage and grass.
      Our results suggest that the financial needs for dairy farms to incorporate more grass in their fodder crop rotation are different from one system to another. In particular, we identify the existence and magnitude of extra production costs linked to the modification of fodder crop rotations (increase in grassland area). In the prospect of designing payments for environmental services such as the Eco-Methane program, and aiming to integrate the maximum number of existing farms into the scheme at the lowest cost, we provide evidence of the variability of dairy systems regarding their willingness to accept. With the perspective of accompanying dairy farms to trigger the adoption of mitigation practices at large scale, it is relevant to take into account this heterogeneity, for instance with a differentiated payment or a higher payment for all. Not considering the feeding strategy of dairy farmers, and in particular the type of fodder system, would lower the attractiveness of a payment scheme for some systems.
      Given the low level of payment in current PES schemes (for instance in the Eco-Methane program), it seems reasonable to assume that participating farms already had good economic profitability and/or feeding practices compatible with emission reductions when entering the program. In the prospect of engaging more farmers and more liters of milk in the feed ration transition, we identify 3 types of dairy farms that could integrate a PES for the reduction of enteric methane emissions. Farms for which reducing enteric emissions is already profitable (negative extra costs) (type 1), farms for which it requires a financial incentive (no extra costs) (type 2), and farms for which it requires high financial support (positive extra costs) (type 3). Our study suggests that French plain-intensive farms broadly correspond to type 1. Although their individual willingness to accept is likely to be low, the program would still need important financial means to offer a payment given the large number of unit of emissions reductions to compensate (high milk productivity). Grass-based plain dairy systems also belong to type 1 farms, but their lower milk productivity involves lower units of emission reduction to compensate. In contrast, dairy farms located in mountainous areas are more likely to correspond to type 2, and require a payment level high enough to make the adoption of more grass-based diet the profitable option over keeping their current practices. For farms with already a high share of grasslands, increasing milk productivity may be an additional profitable lever to reduce enteric emissions per liter. We also classify plain dairy farms with 10 to 30% of corn silage in the fodder area as type 2 farms involving a large number of emission reduction units to compensate given the large number of potential participants and productivity levels. Finally, we identify plain dairy farms from outside the western region with a natural disadvantage in producing milk with more grass as belonging to type 3 on average. While they represent the most expensive dairy farms to support, they are associated with lower enteric methane emissions per liter of milk according to the Eco-Methane baseline indicator. Higher payment levels would recognize the mitigation efforts they already provide.
      However, in the long term, adapting payments to this heterogeneity might not lead to the most optimal allocation to minimize emissions per liter of milk and may introduce a distortion of competition in favor of the most polluting. An optimal scheme should equalise the marginal cost of methane abatement across all farms, for example through a tax on methane emissions. Knowledge of the heterogeneity of abatement costs can then be used to temporarily support the most polluting farms or with a natural disadvantage in improving their performance or in reorienting their activity if they are unable to compete with such a tax.
      The positive effect of grass on methane emissions is likely to be partly offset by a drop in the productivity of the cows. Other sources of omega-3 fatty acids such as linseed could also be integrated in the feed ration to maintain productivity (
      • Fuentes M.C.
      • Calsamiglia S.
      • Sánchez C.
      • González A.
      • Newbold J.R.
      • Santos J.E.P.
      • Rodríguez-Alcalá L.M.
      • Fontecha J.
      Effect of extruded linseed on productive and reproductive performance of lactating dairy cows.
      ). Feed complementation is however likely to be a more expensive lever to implement for most farmers, as they are often not produced in farms. Due to the absence of data, the extra-costs of increasing the use of omega-3 sources rich complements could not be estimated. Furthermore, for the Eco-Methane payment to efficiently subsidize the reduction of enteric methane emissions through an increase in grassland areas, the payment will have to cover both the additional costs per liter of milk and the other extra costs per hectare of grass. Beyond impacting production costs per unit of milk, a new hectare of grassland can have a direct effect on farm costs. In this study, we consider only variable costs (intermediate consumption). There may also be fixed costs (specific machinery for grass cultivation, buildings for storage) or other constraints (access to land) increasing the overall extra costs of participation. A study considering all farm costs (variable and fixed costs) found higher GHG emissions abatement costs per liter of milk in large farms (with a high number of dairy cows) compared with smaller farms (
      • Njuki E.
      • Bravo-Ureta B.E.
      • Mukherjee D.
      The good and the bad: environmental efficiency in northeastern U.S. dairy farming.
      ).

      CONCLUSIONS

      Two aspects of PES design for reducing enteric emissions per liter of milk are particularly crucial to favor its environmental performance. First, the choice of emission indicator should measure farmers' efforts on both cows' productivity and diet, while being easily and regularly implemented on farm at low cost. Second, evaluating farmers' willingness to accept and its variability according to farm type is necessary to define the optimal payment level and scheme's budget ensuring sufficient participation. For a given productivity level, producing milk with more grass have a different impact on milk variable production costs depending on the production basin and fodder system. This research provide more insights into the impact of methane emissions reduction on the production costs of livestock farms, and how to improve support for pressing abatement measures and contribute effectively to achieving climate targets.

      ACKNOWLEDGMENTS

      This research is funded by the Horizon 2020 programme of the European Union (EU) under Grant Agreement No. 817949 (CONSOLE project, https://console-project.eu/). The authors gratefully acknowledge the Bleu-Blanc-Cœur association for providing reference data on the Eco-Methane programme. We also thank the 2 anonymous referees for their constructive comments that greatly helped improving the revised manuscript. On behalf of all authors, the corresponding author states that there is no conflict of interest.

      APPENDIX

      Cost function with the translog specification:
      lnVC(y,w,Z)=α0+r=1,2βrlnYr+i=1,2αilnwi+h=1,2,3δhlnZh+β11lnY1lnY1+i=1,2j=1,2αijlnwilnwj+h=1,2,3n=1,2,3δhnlnZhlnZn+i=1,2ζi1lnwilnY1+k=1,2,3ρk1lnZklnY1+i=1,2k=1,2,3νiklnwilnZk


      Steps to impose homogeneity of degree 1 in input prices:
      lnVC(y,λw,X)=lnλ+lnVC(y,w,X)


      lnVC(y,λw,Z)=α0+r=1,2βrlnYr+h=1,2,3δhlnZh+i=1,2αi(lnλ+lnwi)+β11lnY1lnY1+i=1,2j=1,2αij(lnλ+lnwi)(lnλ+lnwj)+h=1,2,3n=1,2,3δhnlnZhlnZn+i=1,2ζi1(lnλ+lnwi)lnY1+k=1,2,3ρk1lnZklnY1+i=1,2k=1,2,3νik(lnλ+lnwi)lnZk


      lnVC(y,λw,Z)=i=1,2αilnλ+i=1,2j=1,2αij(lnλ)2+i=1,2j=1,2αijlnλlnwj+i=1,2j=1,2αijlnλlnwi+i=1,2ζi1lnλlnY1+i=1,2k=1,2,3νiklnλlnZk+α0+r=1,2βrlnYr+h=1,2,3δhlnZh+i=1,2αilnwi+β11lnY1lnY1+i=1,2j=1,2αijlnwilnwj+h=1,2,3n=1,2,3δhnlnZhlnZn+i=1,2ζi1lnwilnY1+k=1,2,3ρk1lnZklnY1+i=1,2k=1,2,3νiklnwilnZk


      lnVC(y,λw,Z)=lnVC(y,w,Z)+i=1,2αilnλ+i=1,2j=1,2αij(lnλ)2+i=1,2j=1,2αijlnλlnwj+i=1,2j=1,2αijlnλlnwi+i=1,2ζi1lnλlnY1+i=1,2k=1,2,3νiklnλlnZk


      Parameter constraints to impose to our system:
      i=1,2αi=1;i(1,2),j=1,2αij=0;(i,j)(1,2),αij=αji;i=1,2ζi1=0;k(1,2,3),i=1,2νik=0


      Homogeneous function to estimate:
      lnVC(y,w,Z)w1=α0+r=1,2βrlnYr+α2lnw2w1+h=1,2,3δhlnZh+β11lnY1lnY1+α22(lnw2w1)2+h=1,2,3n=1,2,3δhnlnZhlnZn+ζ12lnY1(lnw2w1)+k=1,2,3ρk1lnZklnY1+k=1,2,3ν2klnw2w1lnZk


      Input demand functions: applying Shephard' lemma:
      VC(y,w,Z)wi=Xi


      lnVC(y,w,Z)lnwi=VC(y,w,Z)wiwiVC(y,w,Z)=XiwiVC(y,w,Z)=Si


      S1=1α22α22lnw2w1ζ12lnY1k=1,2,3ν2klnZk


      S2=α2+2α22lnw2w1+ζ12lnY1+k=1,2,3ν2klnZk


      Figure thumbnail fx1
      Figure A1Ratio of the Eco-Methane baseline and Tier 2 enteric methane emissions indicators according to milk productivity for different Eco-Methane scenarios. Source: The authors, based on French FADN and BBC data.
      Table A1.Eco-Methane scenarios attribution to the old French administrative Regions using the 2018 annual dairy survey
      Administrative RegionDepartmentScenariosProduction share (volume)Attributed scenarios
      Ile de France1, 3 or 51, 3 or 5
      Champagne Ardennes1, 3 or 51, 3 or 5
      Picardie1, 3 or 51, 3 or 5
      Haute Normandie1, 3 or 51, 3 or 5
      Centre1, 3 or 51, 3 or 5
      Basse Normandie2, 4, 62, 4 or 6
      Bourgogne1, 3 or 51, 3 or 5
      Nord Pas De Calais1, 3 or 51, 3 or 5
      Lorraine541, 3 or 559%1, 3 or 5
      551, 3 or 5
      571, 3 or 5
      887 or 1141%
      Alsace1, 3 or 51, 3 or 5
      Franche-Comté257 or 1076%7 or 10
      397 or 10
      701, 3 or 524%
      Pays de la Loire2, 4 or 62, 4 or 6
      Bretagne1, 3 or 51, 3 or 5
      Poitou-Charentes1, 3 or 51, 3 or 5
      Aquitaine1, 3 or 51, 3 or 5
      Midi-Pyrénées467 or 860%7 or 8
      127 or 8
      91, 3 or 540%
      311, 3 or 5
      321, 3 or 5
      651, 3 or 5
      811, 3 or 5
      821, 3 or 5
      Limousin1, 3 or 51, 3 or 5
      Rhône Alpes11, 3 or 575%1, 3 or 5
      71, 3 or 5
      261, 3 or 5
      381, 3 or 5
      421, 3 or 5
      691, 3 or 5
      737 or 825%
      74
      Auvergne7 or 87 or 8
      Languedoc-Roussillon111, 3 or 510%7 or 8
      301, 3 or 5
      341, 3 or 5
      661, 3 or 5
      487 or 890%
      Provence Alpes Côte d'Azur57 or 1187%7 or 11
      41, 3 or 513%
      61, 3 or 5
      131, 3 or 5
      831, 3 or 5
      841, 3 or 5
      Source: The authors, based on the French 2018 annual dairy survey and BBC data.
      Table A2.Distribution of the variation of milk production volumes and of the share of farm outputs in the total gross product
      Q1Q2MeanQ3
      % variation of milk production volume relative to the farm mean over 2016–20181.643.494.856.63
      Share of milk production in the total gross product69.1977.7476.5785.15
      Share of crop production in the total gross product0.203.996.3410.98
      Share of other livestock products in the total gross product0.000.000.870.00
      Share of fuel expenses in intermediate consumption3.604.604.886.00
      Share of cattle feed expenses in intermediate consumption19.1025.2024.6331.00
      Share of fertilizer expenses in intermediate consumption3.295.295.787.76
      Source: The authors, based on French FADN data.
      Table A3.Violation of the theoretical properties of the variable cost function (%)
      FrancePlains of the western regionPlains outside the western regionMountainsPlains with more than 30% of corn silage in the fodder areaPlains with 10 to 30% of corn silage in the fodder areaPlains with less than 10% of corn silage in the fodder area
      VC/VCY1<0Y1<00.00.00.00.00.00.00.0
      VC/VCZ1<0Z1<00.22.34.226.321.551.024.4
      VC/VCZ2<0Z2<086.019.786.568.745.059.964.5
      VC/VCZ3<0Z3<038.341.429.216.939.221.371.8
      VC/VCW1>0W1>00.00.00.00.00.00.00.0
      VC/VCW2>0W2>00.00.00.00.00.00.00.0
      2VC/2VCW12<0W12<00.00.00.00.00.00.00.0
      2VC/2VCW22<0W22<00.00.00.00.00.00.00.0
      Table A4.Result of the system estimation for French dairy farms (n = 2,205)
      CoefficientEstimateVariable cost equationFuel cost share functionFeed cost share function
      α03.398*** (0.658)Constant
      β11.156*** (0.202)Ln (milk production)
      α20.457*** (0.029)Ln (feed price/fuel price)1- constantconstant
      δ10.007 (0.108)Ln (grassland)
      δ20.086 (0.065)Ln (capital)
      δ30.739*** (0.214)Ln (labor)
      β11−0.054 (0.042)0.5 *Ln (milk production)2
      ρ11−0.002 (0.020)Ln (milk production) *Ln (grassland)
      ρ12−0.026+ (0.015)Ln (milk production) *Ln (capital)
      ρ13−0.151*** (0.040)Ln (milk production) *Ln (labor)
      ζ120.024*** (0.003)Ln (milk production) *Ln (feed price/fuel price)-Ln (milk production)Ln (milk production)
      α220.002 (0.002)0.5 *Ln (feed price/fuel price)2-Ln (feed price/fuel price)Ln (feed price/fuel price)
      ν210.004* (0.002)Ln (feed price/fuel price) *Ln (grassland)-Ln (grassland)Ln (grassland)
      ν22−0.003** (0.001)Ln (feed price/fuel price) *Ln (capital)-Ln (capital)Ln (capital)
      ν23−0.009** (0.003)Ln (feed price/fuel price) *Ln (labor)-Ln (labor)Ln (labor)
      δ11−0.014 (0.013)0.5 *Ln (grassland)2
      δ120.001 (0.009)Ln (grassland) *Ln (capital)
      δ13−0.014 (0.023)Ln (grassland) *Ln (labor)
      δ220.014* (0.007)0.5 *Ln (capital)2
      δ230.035* (0.016)Ln (capital) *Ln (labor)
      δ330.091 (0.061)0.5 *Ln (labor)2
      β20.000*** (0.000)other productions
      θ10.002*** (0.000)utilized agricultural area
      γ1−0.199 (0.326)Champagne-Ardennes
      γ2−0.146 (0.325)Picardie
      γ3−0.121 (0.325)Haute Normandie
      γ4−0.086 (0.326)Centre
      γ5−0.119 (0.325)Basse Normandie
      γ6−0.107 (0.329)Bourgogne
      γ7−0.107 (0.325)Nord-Pas-De-Calais
      γ8−0.119 (0.325)Lorraine
      γ9−0.094 (0.327)Alsace
      γ10−0.123 (0.325)Franche-Comté
      γ11−0.156 (0.325)Pays de la Loire
      γ12−0.207 (0.324)Bretagne
      γ13−0.111 (0.326)Poitou-Charentes
      γ14−0.080 (0.326)Aquitaine
      γ15−0.257 (0.325)Midi-Pyrénées
      γ16−0.200 (0.327)Limousin
      γ17−0.104 (0.325)Rhône-Alpes
      γ18−0.247 (0.325)Auvergne
      γ19−0.233 (0.327)Languedoc-Roussillon
      γ20−0.281 (0.332)Provence-Alpes-Côte d'Azur
      Baseline: Ile-de-France
      μ160.024* (0.011)Year 2016
      μ170.003 (0.010)Year 2017
      Baseline: 2018
      R20.912−306.535−14.425
      Standard errors in parentheses. + P < 0.10, * P < 0.05, ** P < 0.01, *** P < 0.001.
      Source: The authors, based on French FADN data.
      Table A5.Result of the system estimation for plain dairy farms in the western region (n = 645)
      CoefficientEstimateVariable cost equationFuel cost share functionFeed cost share function
      α05.422*** (1.075)Constant
      β10.413 (0.419)Ln (milk production)
      α20.423*** (0.048)Ln (feed price/fuel price)1- constantconstant
      δ10.423* (0.210)Ln (grassland)
      δ2−0.078 (0.146)Ln (capital)
      δ30.744+ (0.400)Ln (labor)
      β110.015 (0.100)0.5 *Ln (milk production)2
      ρ11−0.046 (0.043)Ln (milk production) *Ln (grassland)
      ρ120.044 (0.038)Ln (milk production) *Ln (capital)
      ρ13−0.100 (0.085)Ln (milk production) *Ln (labor)
      ζ120.041*** (0.005)Ln (milk production) *Ln (feed price/fuel price)- Ln (milk production)Ln (milk production)
      α220.001 (0.004)0.5 *Ln (feed price/fuel price)2-Ln (feed price/fuel price)Ln (feed price/fuel price)
      ν21−0.018*** (0.003)Ln (feed price/fuel price) *Ln (grassland)-Ln (grassland)Ln (grassland)
      ν22−0.005* (0.002)Ln (feed price/fuel price) *Ln (capital)-Ln (capital)Ln (capital)
      ν230.003 (0.005)Ln (feed price/fuel price) *Ln (labor)-Ln (labor)Ln (labor)
      δ11−0.047 (0.042)0.5 *Ln (grassland)2
      δ120.004 (0.021)Ln(grassland)*Ln(capital)
      δ130.030 (0.051)Ln(grassland)*Ln(labor)
      δ22−0.031+ (0.018)0.5*Ln(capital)2
      δ23−0.063* (0.032)Ln(capital)*Ln(labor)
      δ330.143 (0.106)0.5*Ln(labor)2
      β20.000*** (0.000)other productions
      θ10.002*** (0.000)utilized agricultural area
      γ11−0.043* (0.019)Pays de la Loire
      γ12−0.099*** (0.019)Bretagne
      Baseline: Basse Normandie
      μ160.003 (0.019)Year 2016
      μ17−0.018 (0.017)Year 2017
      Baseline: 2018
      R20.915−362.617−18.850
      Standard errors in parentheses. + P < 0.10, * P < 0.05, ** P < 0.01, *** P < 0.001.
      Source: The authors, based on French FADN data.
      Table A6.Result of the system estimation for plain dairy farms outside the western region (n = 975)
      CoefficientEstimateVariable cost equationFuel cost share functionFeed cost share function
      α05.254*** (0.929)Constant
      β10.668* (0.293)Ln (milk production)
      α20.404*** (0.043)Ln (feed price/fuel price)1- constantconstant
      δ1−0.310* (0.156)Ln (grassland)
      δ20.101 (0.087)Ln (capital)
      δ31.558*** (0.321)Ln (labor)
      β110.012 (0.057)0.5*Ln (milk production)2
      ρ110.047+ (0.027)Ln (milk production) *Ln (grassland)
      ρ12−0.038* (0.019)Ln (milk production) *Ln (capital)
      ρ13−0.297*** (0.056)Ln (milk production) *Ln (labor)
      ζ120.033*** (0.004)Ln (milk production) *Ln (feed price/fuel price)- Ln (milk production)Ln (milk production)
      α220.001 (0.004)0.5 *Ln (feed price/fuel price)2-Ln (feed price/fuel price)Ln (feed price/fuel price)
      ν210.009*** (0.002)Ln (feed price/fuel price) *Ln (grassland)-Ln (grassland)Ln (grassland)
      ν22−0.001 (0.002)Ln (feed price/fuel price) *Ln (capital)-Ln (capital)Ln (capital)
      ν23−0.032*** (0.005)Ln (feed price/fuel price) *Ln (labor)-Ln (labor)Ln (labor)
      δ11−0.006 (0.016)0.5 *Ln (grassland)2
      δ12−0.001 (0.012)Ln (grassland) *Ln (capital)
      δ13−0.040 (0.031)Ln (grassland) *Ln (labor)
      δ220.024* (0.009)0.5 *Ln (capital)2
      δ230.087*** (0.022)Ln (capital) *Ln (labor)
      δ330.022 (0.089)0.5 *Ln (labor)2
      β20.000*** (0.000)other productions
      θ10.001*** (0.000)utilized agricultural area
      γ1−0.175 (0.283)Champagne-Ardennes
      γ2−0.100 (0.282)Picardie
      γ3−0.078 (0.282)Haute Normandie
      γ4−0.010 (0.283)Centre
      γ6−0.067 (0.285)Bourgogne
      γ7−0.079 (0.282)Nord-Pas-De-Calais
      γ8.−0.078 (0.282)Lorraine
      γ9−0.059 (0.283)Alsace
      γ13−0.052 (0.283)Poitou-Charentes
      γ14−0.063 (0.283)Aquitaine
      γ16−0.163 (0.285)Limousin
      γ17−0.093 (0.282)Rhône-Alpes
      Baseline: Ile-de-France
      μ160.057** (0.017)Year 2016
      μ170.017 (0.016)Year 2017
      Baseline: 2018
      R20.911−276.665−12.377
      Standard errors in parentheses. + P < 0.10, * P < 0.05, ** P < 0.01, *** P < 0.001.
      Source: The authors, based on French FADN data.
      Table A7.Result of the system estimation for mountain dairy farms (n = 585)
      CoefficientEstimateVariable cost equationFuel cost share functionFeed cost share function
      α01.397 (1.152)Constant
      β11.285** (0.427)Ln (milk production)
      α20.408*** (0.058)Ln (feed price/fuel price)1- constantconstant
      δ1514* (0.257)Ln (grassland)
      δ20.270* (0.136)Ln (capital)
      δ3−0.193 (0.437)Ln (labor)
      β110.001 (0.094)0.5 *Ln (milk production)2
      ρ11−0.060 (0.047)Ln (milk production) *Ln (grassland)
      ρ12−0.055+ (0.032)Ln (milk production) *Ln (capital)
      ρ13−0.093 (0.083)Ln (milk production) *Ln (labor)
      ζ120.025*** (0.006)Ln (milk production) *Ln (feed price/fuel price)- Ln (milk production)Ln (milk production)
      α220.006 (0.005)0.5 *Ln (feed price/fuel price)2-Ln (feed price/fuel price)Ln (feed price/fuel price)
      ν210.019*** (0.004)Ln (feed price/fuel price) *Ln (grassland)-Ln (grassland)Ln (grassland)
      ν22−0.013*** (0.003)Ln (feed price/fuel price) *Ln (capital)-Ln (capital)Ln (capital)
      ν23−0.016* (0.006)Ln (feed price/fuel price) *Ln (labor)-Ln (labor)Ln (labor)
      δ11−0.062 (0.049)0.5 *Ln (grassland)2
      δ12−0.014 (0.020)Ln (grassland) *Ln (capital)
      δ130.081 (0.055)Ln (grassland) *Ln (labor)
      δ220.025* (0.012)0.5* Ln (capital)2
      δ230.078* (0.039)Ln (capital)* Ln (labor)
      δ33−0.074 (0.145)0.5 *Ln (labor)2
      β20.000*** (0.000)other productions
      θ10.002*** (0.000)utilized agricultural area
      γ100.109*** (0.030)Franche-Comté
      γ180.006 (0.025)Auvergne
      γ190.032 (0.050)Languedoc-Roussillon
      γ200.017 (0.079)Provence-Alpes-Côte d'Azur
      Baseline: Midi-Pyrénées
      μ160.032 (0.023)Year 2016
      μ170.030 (0.020)Year 2017
      Baseline: 2018
      R20.910−289.950−12.903
      Standard errors in parentheses. + P < 0.10, * P < 0.05, ** P < 0.01, *** P < 0.001.
      Source: The authors, based on French FADN data.
      Table A8.Result of the system estimation for plain dairy farms with more than 30% of corn in the fodder system (n = 767)
      CoefficientEstimateVariable cost equationFuel cost share functionFeed cost share function
      α06.811*** (1.111)Constant
      β1−0.265 (0.449)Ln (milk production)
      α20.388*** (0.046)Ln (feed price/fuel price)1- constantconstant
      δ10.283 (0.233)Ln (grassland)
      δ20.320* (0.140)Ln (capital)
      δ31.445** (0.437)Ln (labor)
      β110.262* (0.105)0.5 *Ln (milk production)2
      ρ11−0.092+ (0.049)Ln (milk production) *Ln (grassland)
      ρ12−0.078* (0.037)Ln (milk production) *Ln (capital)
      ρ13−0.235* (0.096)Ln (milk production) *Ln (labor)
      ζ120.046*** (0.006)Ln (milk production) *Ln (feed price/fuel price)- Ln (milk production)Ln (milk production)
      α220.002 (0.004)0.5*Ln(feed price/fuel price)2-Ln (feed price/fuel price)Ln (feed price/fuel price)
      ν21−0.013*** (0.003)Ln (feed price/fuel price) *Ln (grassland)-Ln (grassland)Ln (grassland)
      ν22−0.007** (0.002)Ln (feed price/fuel price) *Ln (capital)-Ln (capital)Ln (capital)
      ν23−0.008 (0.005)Ln (feed price/fuel price) *Ln (labor)-Ln (labor)Ln (labor)
      δ11−0.001 (0.019)0.5* Ln (grassland)2
      δ120.063** (0.019)Ln (grassland) *Ln (capital)
      δ13−0.061 (0.040)Ln(grassland)*Ln(labor)
      δ22−0.005 (0.016)0.5*Ln(capital)2
      δ230.018 (0.031)Ln(capital)*Ln(labor)
      δ330.237* (0.102)0.5*Ln(labor)2
      β20.000*** (0.000)other productions
      θ10.002*** (0.000)utilized agricultural area
      γ1−0.194 (0.300)Champagne-Ardenne
      γ2−0.241 (0.285)Picardie
      γ3−0.236 (0.285)Haute Normandie
      γ4−0.178 (0.286)Centre
      γ5−0.214 (0.284)Basse Normandie
      γ6−0.163 (0.293)Bourgogne
      γ7−0.185 (0.284)Nord-Pas-De-Calais
      γ8−0.199 (0.287)Lorraine
      γ9−0.217 (0.289)Alsace
      γ11−0.228 (0.284)Pays de la Loire
      γ12−0.273 (0.284)Bretagne
      γ13−0.240 (0.287)Poitou-Charentes
      γ14−0.178 (0.286)Aquitaine
      γ16−0.236 (0.308)Limousin
      γ17−0.120 (0.290)Rhône-Alpes
      Baseline: Ile-de-France
      μ160.034+ (0.017)Year 2016
      μ170.004 (0.016)Year 2017
      Baseline: 2018
      R20.902−320.701−17.652
      Standard errors in parentheses. + P < 0.10, * P < 0.05, ** P < 0.01, *** P < 0.001.
      Source: The authors, based on French FADN data.
      Table A9.Result of the system estimation for plain dairy farms with 10 to 30% of corn in the fodder system (n = 574)
      CoefficientEstimateVariable cost equationFuel cost share functionFeed cost share function
      α06.926*** (1.244)Constant
      β1−0.008 (0.460)Ln (milk production)
      α20.426*** (0.051)Ln (feed price/fuel price)1- constantconstant
      δ1−0.071 (0.289)Ln (grassland)
      δ20.100 (0.147)Ln (capital)
      δ30.992* (0.445)Ln (labor)
      β110.083 (0.105)0.5 *Ln (milk production)2
      ρ110.031 (0.063)Ln (milk production) *Ln (grassland)
      ρ120.003 (0.035)Ln (milk production) *Ln (capital)
      ρ13−0.190* (0.088)Ln (milk production) *Ln (labor)
      ζ120.035*** (0.006)Ln (milk production) *Ln (feed price/fuel price)- Ln (milk production)Ln (milk production)
      α22−0.001 (0.004)0.5 *Ln (feed price/fuel price)2-Ln (feed price/fuel price)Ln (feed price/fuel price)
      ν21−0.004 (0.004)Ln (feed price/fuel price) *Ln (grassland)-Ln (grassland)Ln (grassland)
      ν22−0.003 (0.002)Ln (feed price/fuel price) *Ln (capital)-Ln (capital)Ln (capital)
      ν230.009 (0.006)Ln (feed price/fuel price) *Ln (labor)-Ln (labor)Ln (labor)
      δ110.009 (0.060)0.5 *Ln (grassland)2
      δ12−0.023 (0.022)Ln (grassland) *Ln (capital)
      δ13−0.008 (0.063)Ln (grassland) *Ln (labor)
      δ22−0.002 (0.017)0.5 *Ln (capital)2
      δ230.008 (0.035)Ln (capital) *Ln (labor)
      δ330.076 (0.119)0.5 *Ln (labor)2
      β20.000*** (0.000)other productions
      θ10.001* (0.000)utilized agricultural area
      γ1−0.021 (0.053)Champagne-Ardenne
      γ20.024 (0.071)Picardie
      γ3−0.039 (0.047)Haute Normandie
      γ40.130 (0.119)Centre
      γ5−0.055+ (0.029)Basse Normandie
      γ6−0.036 (0.066)Bourgogne
      γ7−0.011 (0.045)Nord-Pas-De-Calais
      γ8−0.020 (0.039)Lorraine
      γ90.004 (0.059)Alsace
      γ11−0.063* (0.029)Pays de la Loire
      γ12−0.163*** (0.029)Bretagne
      γ130.029 (0.053)Poitou-Charentes
      γ14−0.038 (0.053)Aquitaine
      γ16−0.037 (0.062)Limousin
      Baseline: Ile-de-France
      μ160.006 (0.019)Year 2016
      μ17−0.027 (0.018)Year 2017
      Baseline: 2018
      R20.908−341.414−18.898
      Standard errors in parentheses. + P < 0.10, * P < 0.05, ** P < 0.01, *** P < 0.001.
      Source: The authors, based on French FADN data.
      Table A10.Result of the system estimation for plain dairy farms with less than 10% of corn in the fodder system (n = 279)
      CoefficientEstimateVariable cost equationFuel cost share functionFeed cost share function
      α0−5.524* (2.228)Constant
      β12.652*** (0.786)Ln (milk production)
      α20.441*** (0.096)Ln (feed price/fuel price)1- constantconstant
      δ12.645*** (0.552)Ln (grassland)
      δ20.086 (0.210)Ln (capital)
      δ3−1.532* (0.731)Ln (labor)
      β11−0.190 (0.172)0.5*Ln (milk production)2
      ρ11−0.269*** (0.080)Ln (milk production) *Ln (grassland)
      ρ12−0.008 (0.054)Ln (milk production) *Ln (capital)
      ρ130.163 (0.157)Ln (milk production) *Ln (labor)
      ζ120.039*** (0.010)Ln (milk production) *Ln (feed price/fuel price)- Ln (milk production)Ln (milk production)
      α22−0.001 (0.009)0.5 *Ln (feed price/fuel price)2-Ln (feed price/fuel price)Ln (feed price/fuel price)
      ν21−0.002 (0.008)Ln (feed price/fuel price) *Ln (grassland)-Ln (grassland)Ln (grassland)
      ν22−0.002 (0.003)Ln (feed price/fuel price) *Ln (capital)-Ln (capital)Ln (capital)
      ν23−0.039*** (0.011)Ln (feed price/fuel price) *Ln (labor)-Ln (labor)Ln (labor)
      δ11−0.260* (0.118)0.5 *Ln (grassland)2
      δ12−0.055 (0.047)Ln (grassland) *Ln (capital)
      δ130.280** (0.105)Ln (grassland) *Ln (labor)
      δ220.037+ (0.019)0.5 *Ln (capital)2
      δ230.004 (0.057)Ln (capital) *Ln (labor)
      δ33−0.492*** (0.182)0.5 *Ln (labor)2
      β20.000*** (0.000)other productions
      θ10.002+ (0.001)utilized agricultural area
      γ1−0.184** (0.060)Champagne-Ardennes
      γ2−0.131+ (0.076)Picardie
      γ3−0.084 (0.086)Haute Normandie
      γ4−0.108 (0.156)Centre
      γ50.028 (0.052)Basse Normandie
      γ7−0.249** (0.086)Nord-Pas-De-Calais
      γ8−0.088 (0.054)Lorraine
      γ90.368+ (0.208)Alsace
      γ11−0.164*** (0.047)Pays de la Loire
      γ12−0.060 (0.069)Bretagne
      γ14−0.012 (0.191)Aquitaine
      γ16−0.436** (0.158)Limousin
      Baseline: Ile-de-France
      μ160.006 (0.034)Year 2016
      μ170.004 (0.030)Year 2017
      Baseline: 2018
      R20.894−282.326−9.027
      Standard errors in parentheses. + P < 0.10, * P < 0.05, ** P < 0.01, *** P < 0.001.
      Source: The authors, based on French FADN data.

      REFERENCES

      1. Agreste. 2019. Enquête annuelle laitière 2018. French 2018 annual dairy survey (in French). Agreste Chiffres et Données 13.

        • Allen B.
        • Hart K.
        • Radley G.
        • Tucker G.
        • Keenleyside C.
        • Oppermann R.
        • Underwood E.
        • Menadue H.
        • Poux X.
        • Beaufoy G.
        • Herzon I.
        • Povellato A.
        • Vanni F.
        • Prazan J.
        • Hudson T.
        • Yellachich N.
        Biodiversity protection through results based remuneration of ecological achievement. Report Prepared for the European Commission, DG Environment, Contract No ENV.B.2/ETU/2013/0046.
        Institute for European Environmental Policy, London2014
        • Alvarez A.
        • Arias C.
        Diseconomies of size with fixed managerial ability.
        Am. J. Agric. Econ. 2003; 85: 134-142
        • Baylis K.
        • Peplow S.
        • Rausser G.
        • Simon L.
        Agric.-environmental policies in the EU and United States: A comparison.
        Ecol. Econ. 2008; 65: 753-764
        • Bleu-Blanc-Coeur
        Démarche Environnementale : La Démarche Éco-Méthane de Bleu-Blanc-Cœur. Environmental Engagement: The Eco-Methane Programme of Bleu-Blanc-Coeur (in French).
        • Boadi D.
        • Benchaar C.
        • Chiquette J.
        • Massé D.
        Mitigation strategies to reduce enteric methane emissions from dairy cows: Update review.
        Can. J. Anim. Sci. 2004; 84: 319-335
        • Borreani G.
        • Coppa M.
        • Revello-Chion A.
        • Comino L.
        • Giaccone D.
        • Ferlay A.
        • Tabacco E.
        Effect of different feeding strategies in intensive dairy farming systems on milk fatty acid profiles, and implications on feeding costs in Italy.
        J. Dairy Sci. 2013; 96 (24011944): 6840-6855
        • Cain M.
        • Lynch J.
        • Allen M.R.
        • Fuglestvedt J.S.
        • Frame D.J.
        • Macey A.H.
        Improved calculation of warming-equivalent emissions for short-lived climate pollutants.
        npj Climate and Atmospheric Science. 2019; 29: 1-7
        • Chen W.
        • Holden N.M.
        Tiered life cycle sustainability assessment applied to a grazing dairy farm.
        J. Clean. Prod. 2018; 172: 1169-1179
        • Chilliard Y.
        • Martin C.
        • Rouel J.
        • Doreau M.
        Milk fatty acids in dairy cows fed whole crude linseed, extruded linseed, or linseed oil, and their relationship with methane output.
        J. Dairy Sci. 2009; 92 (19762838): 5199-5211
        • Citepa
        Inventaire des émissions de polluants atmosphériques et de gaz à effet de serre en France - Format Secten. National inventory of GHG and air pollutant emissions in France (in French).
        Citepa, Paris2020
        • Citepa
        Organisation et méthodes des inventaires nationaux des émissions atmosphériques en France OMINEA - 17 ème édition. Inventory methodology report of GHG and air pollutant emissions in France - 17th edition (in French).
        Citepa, Paris2020
        • CNIEL
        Observatoire de l'alimentation des vaches laitières: description des 8 principaux systèmes d'élevage. Edition 2015–2018. Dairy cows' feed observatory: description of the 8 main farming systems. 2015–2018 edition (in French).
        CNIEL, Paris2015
        • Coppa M.
        • Chassaing C.
        • Sibra C.
        • Cornu A.
        • Verbič J.
        • Golecký J.
        • Engel E.
        • Ratel J.
        • Boudon A.
        • Ferlay A.
        • Martin B.
        Forage system is the key driver of mountain milk specificity.
        J. Dairy Sci. 2019; 102 (31495613): 10483-10499
        • Del Prado A.
        • Mas K.
        • Pardo G.
        • Gallejones P.
        Modelling the interactions between C and N farm balances and GHG emissions from confinement dairy farms in northern Spain.
        Sci. Total Environ. 2013; 465 (23601287): 156-165
        • Dong Y.
        • Bae H.D.
        • McAllister T.A.
        • Mathison G.W.
        • Cheng K.J.
        Lipid induced depression of methane production and digestibility in the artificial rumen system (RUSITEC).
        Can. J. Anim. Sci. 1997; 77: 269-278
        • Douenne T.
        • Fabre A.
        French attitudes on climate change, carbon taxation and other climate policies.
        Ecol. Econ. 2020; 169106496
        • Dupraz P.
        Policies for the ecological transition of agriculture: the livestock issue.
        Rev. Agric. Food Environ. Stud. 2021; 101: 529-538
        • Duvaleix-Tréguer S.
        • Gaigné C.
        On the nature and magnitude of cost economies in hog production.
        Agric. Econ. 2016; 47: 465-476
        • EEA
        Annual European Union greenhouse gas inventory 1990 – 2018 and inventory report 2020. Submission to the UNFCCC Secretariat.
        EEA, Copenhagen2020
        • Ellis J.L.
        • Kebreab E.
        • Odongo N.E.
        • McBride B.W.
        • Okine E.K.
        • France J.
        Prediction of methane production from dairy and beef cattle.
        J. Dairy Sci. 2007; 90 (17582129): 3456-3466
        • Eugène M.
        • Sauvant D.
        • Nozière P.
        • Viallard D.
        • Oueslati K.
        • Lherm M.
        • Mathias E.
        • Doreau M.
        A new Tier 3 method to calculate methane emission inventory for ruminants.
        J. Environ. Manage. 2019; 231 (30602259): 982-988
        • Fuentes M.C.
        • Calsamiglia S.
        • Sánchez C.
        • González A.
        • Newbold J.R.
        • Santos J.E.P.
        • Rodríguez-Alcalá L.M.
        • Fontecha J.
        Effect of extruded linseed on productive and reproductive performance of lactating dairy cows.
        Livest. Sci. 2008; 113: 144-154
        • Funke F.
        • Mattauch L.
        • van den Bijgaart I.
        • Godfray H.C.J.
        • Hepburn C.
        • Klenert D.
        • Springmann M.
        • Treich N.
        Toward Optimal Meat Pricing: Is It Time to Tax Meat Consumption?.
        Rev. Environ. Econ. Policy. 2022; 16: 219-240
        • Gavrilova O.
        • Leip A.
        • Dong H.
        • Macdonald J.D.
        • Gomez Bravo C.A.
        • Amon B.
        • Barahona Rosales R.
        • Agustin del Prado A.
        • Aparecida de Lima M.
        • Oyhantcabal W.
        • van der Weerden T.J.
        • Widiawati Y.
        Emissions from Livestock and Manure Management.
        in: 2019 Refinement to the 2006 IPCC Guidelines for National Greenhouse Gas Inventories. IPCC, Geneva2019: 1-209
        • Gollnow S.
        • Lundie S.
        • Moore A.D.
        • McLaren J.
        • van Buuren N.
        • Stahle P.
        • Christie K.
        • Thylmann D.
        • Rehl T.
        Carbon footprint of milk production from dairy cows in Australia.
        Int. Dairy J. 2014; 37: 31-38
        • Grainger C.
        • Beauchemin K.A.
        Can enteric methane emissions from ruminants be lowered without lowering their production?.
        Anim. Feed Sci. Technol. 2011; 166–167: 308-320
        • Guerci M.
        • Knudsen M.T.
        • Bava L.
        • Zucali M.
        • Schönbach P.
        • Kristensen T.
        Parameters affecting the environmental impact of a range of dairy farming systems in Denmark, Germany and Italy.
        J. Clean. Prod. 2013; 54: 133-141
        • Hagemann M.
        • Hemme T.
        • Ndambi A.
        • Alqaisi O.
        • Sultana M.N.
        Benchmarking of greenhouse gas emissions of bovine milk production systems for 38 countries.
        Anim. Feed Sci. Technol. 2011; 166–167: 46-58
        • Hanrahan L.
        • McHugh N.
        • Hennessy T.
        • Moran B.
        • Kearney R.
        • Wallace M.
        • Shalloo L.
        Factors associated with profitability in pasture-based systems of milk production.
        J. Dairy Sci. 2018; 101 (29525299): 5474-5485
        • Henderson B.
        • Verma M.
        Global assessment of the carbon leakage implications of carbon taxes on agricultural emissions.
        OECD FOOD, AGRICULTURE AND FISHERIES PAPER, 2021
        • Hennessy D.
        • Delaby L.
        • van den Pol-van Dasselaar A.
        • Shalloo L.
        Increasing grazing in dairy cow milk production systems in Europe.
        Sustainability (Basel). 2020; 12: 1-15
        • Jayasundara S.
        • Worden D.
        • Weersink A.
        • Wright T.
        • VanderZaag A.
        • Gordon R.
        • Wagner-Riddle C.
        Improving farm profitability also reduces the carbon footprint of milk production in intensive dairy production systems.
        J. Clean. Prod. 2019; 229: 1018-1028
        • Kebreab E.
        • Clark K.
        • Wagner-Riddle C.
        • France J.
        Methane and nitrous oxide emissions from Canadian animal agriculture: A review.
        Can. J. Anim. Sci. 2006; 86: 135-158
        • Key N.
        • Tallard G.
        Mitigating methane emissions from livestock: A global analysis of sectoral policies.
        Clim. Change. 2012; 112: 387-414
        • Lambaré P.
        • Dervillé M.
        • You G.
        What will be the conditions of market access for dairy farmers after the end of dairy quotas?.
        Écon. Rurale (Paris). 2018; 2 (in French): 55-71
        • Le S.
        • Jeffrey S.
        • An H.
        Greenhouse Gas Emissions and Technical Efficiency in Alberta Dairy Production : What Are the Trade-Offs?.
        J. Agric. Appl. Econ. 2020; 52: 177-193
        • Lengers B.
        • Britz W.
        • Holm-Müller K.
        Comparison of GHG-Emission indicators for dairy farms with respect to induced abatement costs, accuracy, and feasibility.
        Appl. Econ. Perspect. Policy. 2013; 35: 451-475
        • Lorenz H.
        • Reinsch T.
        • Hess S.
        • Taube F.
        Is low-input dairy farming more climate friendly? A meta-analysis of the carbon footprints of different production systems.
        J. Clean. Prod. 2019; 211: 161-170
        • Martin C.
        • Morgavi D.
        • Doreau M.
        • Jouany J.P.
        Comment réduire la production de méthane chez les ruminants? How can the production of methane by ruminants be reduced?.
        Fourrages (Versailles). 2006; 187 (in French): 283-300
        • Martin C.
        • Morgavi D.P.
        • Doreau M.
        Methane mitigation in ruminants: from microbe to the farm scale.
        Animal. 2010; 4 (22443940): 351-365
        • Martin C.
        • Promiès D.
        • Ferlay A.
        • Rochette Y.
        • Martin B.
        • Chilliard Y.
        • Morgavi D.
        • Doreau M.
        Methane output and rumen microbiota in dairy cows in response to long-term supplementation with linseed or rapeseed of grass silage or pasture based diets.
        in: Proceedings of the New Zealand Society of Animal Production. volume 71. New Zealand Society of Animal Production, Invercargill2011: 243-247
        • Martin C.
        • Rouel J.
        • Jouany J.P.
        • Doreau M.
        • Chilliard Y.
        Methane output and diet digestibility in response to feeding dairy cows crude linseed, extruded linseed, or linseed oil.
        J. Anim. Sci. 2008; 86 (18469051): 2642-2650
        • McFadden D.
        Cost, Revenue and Profit Functions.
        in: Fuss M. McFadden D. Production Economics: A Dual Approach to Theory and Applications, Volume I: The Theory of Production. North Holland, Amsterdam1978: 2-109
        • Moschini G.
        The Cost Structure of Ontario Dairy Farms: A Microeconometric Analysis.
        Can. J. Agric. Econ. 1988; 36: 187-206
        • Mosheim R.
        • Lovell C.A.K.
        Scale economies and inefficiency of U.S. dairy farms.
        Am. J. Agric. Econ. 2009; 91: 777-794
        • Mosnier C.
        • Britz W.
        • Julliere T.
        • De Cara S.
        • Jayet P.A.
        • Havlík P.
        • Frank S.
        • Mosnier A.
        Greenhouse gas abatement strategies and costs in French dairy production.
        J. Clean. Prod. 2019; 236: 1-11
        • Myhre G.
        • Shindell D.
        • Bréon F.-M.
        • Collins W.
        • Fuglestvedt J.
        • Huang J.
        • Koch D.
        • Lamarque J.-F.
        • Lee D.
        • Mendoza B.
        • Nakajima T.
        • Robock A.
        • Stephens G.
        • Takemura T.
        • Zhang H.
        Anthropogenic and Natural Radiative Forcing.
        in: Stocker T.F. Qin D. Plattner G.-K. Tignor M. Allen S.K. Boschung J. Nauels A. Xia Y. Bex V. Midgley P.M. Climate Change 2013 the Physical Science Basis: Working Group I Contribution to the Fifth Assessment Report of the Intergovernmental Panel on Climate Change. Cambridge University Press, Cambridge, New York2013: 659-740
        • Negussie E.
        • de Haas Y.
        • Dehareng F.
        • Dewhurst R.J.
        • Dijkstra J.
        • Gengler N.
        • Morgavi D.P.
        • Soyeurt H.
        • van Gastelen S.
        • Yan T.
        • Biscarini F.
        Invited review: Large-scale indirect measurements for enteric methane emissions in dairy cattle: A review of proxies and their potential for use in management and breeding decisions.
        J. Dairy Sci. 2017; 100 (28161178): 2433-2453
        • Nehring R.
        • Gillespie J.
        • Sandretto C.
        • Hallahan C.
        Small U.S. dairy farms: can they compete?.
        Agric. Econ. 2009; 40: 817-825
        • Njuki E.
        • Bravo-Ureta B.E.
        The economic costs of environmental regulation in U.S. Dairy farming: a directional distance function approach.
        Am. J. Agric. Econ. 2015; 97: 1087-1106
        • Njuki E.
        • Bravo-Ureta B.E.
        • Mukherjee D.
        The good and the bad: environmental efficiency in northeastern U.S. dairy farming.
        Agric. Resour. Econ. Rev. 2016; 45: 22-43
        • O'Brien D.
        • Shalloo L.
        • Patton J.
        • Buckley F.
        • Grainger C.
        • Wallace M.
        Evaluation of the effect of accounting method, IPCC v. LCA, on grass-based and confinement dairy systems' greenhouse gas emissions.
        Animal. 2012; 6 (23031525): 1512-1527
        • Rogissart L.
        • Postic S.
        • Grimault J.
        La Contribution Climat Energie en France : fonctionnement, revenus et exonérations. The Climate-Energy Contribution in France: operation, revenues and exemptions.
        Point Climat. 2018; 56 (in French): 1-7
        • Sauvant D.
        • Giger-Reverdin S.
        • Serment A.
        • Broudiscou L.
        Influences des régimes et de leur fermentation dans le rumen sur la production de méthane par les ruminants.
        INRA Prod. Anim. 2011; 24 (in French): 433-446
        • Senga Kiessé T.
        • Corson M.S.
        • Wilfart A.
        Analysis of milk production and greenhouse gas emissions as a function of extreme variations in forage production among French dairy farms.
        J. Environ. Manage. 2022; 307 (35078066)114537
        • Singbo A.
        • Larue B.
        Scale economies, technical efficiency, and the sources of total factor productivity growth of Quebec dairy farms.
        Can. J. Agric. Econ. 2016; 64: 339-363
        • Sobczyński T.
        • Klepacka A.M.
        • Revoredo-Giha C.
        • Florkowski W.J.
        Dairy farm cost efficiency in leading milk-producing regions in Poland.
        J. Dairy Sci. 2015; 98 (26476947): 8294-8307
        • Stetter C.
        • Sauer J.
        Greenhouse Gas Emissions and Eco-Performance at Farm Level: A Parametric Approach.
        Environ. Resour. Econ. 2022; 81: 617-647
        • The World Bank
        Carbon Pricing Dashboard | Up-to-Date Overview of Carbon Pricing Initiatives.
        The World Bank, 2021
        • Tsionas E.G.
        • Kumbhakar S.C.
        • Malikov E.
        Estimation of Input Distance Functions : A System Approach.
        Am. J. Agric. Econ. 2015; 97: 1478-1493
        • UNFCCC
        Joint Implementation Project FR1000365: “Réduction des émissions de méthane d'origine digestive par l'apport dans l'alimentation des vaches laitières de sources naturelles en Acide Alpha Linolénique (ALA)” (in French).
        UNFCC, Rio de Janeiro, New York2016
      2. Weill, P., G. Chesneau, Y. Chilliard, M. Doreau, and C. Martin. 2009. Method to evaluate the quantity of methane produced by a dairy ruminant and method for decreasing and controlling this quantity. Valorex, assignee. Pat. No. WO 2009/156453 A1.

        • Weill P.
        • Schmitt B.
        • Chesneau G.
        • Daniel N.
        • Safraou F.
        • Legrand P.
        Effects of introducing linseed in livestock diet on blood fatty acid composition of consumers of animal products.
        Ann. Nutr. Metab. 2002; 46 (12378041): 182-191
        • Wilkes A.
        • Wassie S.
        • Odhong' C.
        • Fraval S.
        • van Dijk S.
        Variation in the carbon footprint of milk production on smallholder dairy farms in central Kenya.
        J. Clean. Prod. 2020; 265: 1-15
        • Wimmer S.
        • Sauer J.
        Diversification economies in dairy farming - Empirical evidence from Germany.
        Eur. Rev. Agric. Econ. 2020; 47: 1338-1365
        • Wirsenius S.
        • Hedenus F.
        • Mohlin K.
        Greenhouse gas taxes on animal food products: Rationale, tax scheme and climate mitigation effects.
        Clim. Change. 2011; 108: 159-184
        • Wunder S.
        Revisiting the concept of payments for environmental services.
        Ecol. Econ. 2015; 117: 234-243