Abstract: Croplands and grasslands are agricultural systems that contribute to land–atmosphere exchanges of carbon (C). We evaluated and compared gross primary production (GPP), ecosystem respiration (RECO), net ecosystem exchange (NEE) of CO2, and two derived outputs – C use efficiency (CUE=-NEE/GPP) and C emission intensity (IntC= -NEE/Offtake [grazed or harvested biomass]). The outputs came from 23 models (11 crop-specific, eight grassland-specific, and four models covering both systems) at three cropping sites over several rotations with spring and winter cereals, soybean and rapeseed in Canada, France and India, and two temperate permanent grasslands in France and the United Kingdom. The models were run independently over multi-year simulation periods in five stages (S), either blind with no calibration and initialization data (S1), using historical management and climate for initialization (S2), calibrated against plant data (S3), plant and soil data together (S4), or with the addition of C and N fluxes (S5). Here, we provide a framework to address methodological uncertainties and contextualize results. Most of the models overestimated or underestimated the C fluxes observed during the growing seasons (or the whole years for grasslands), with substantial differences between models. For each simulated variable, changes in the multi-model median (MMM) from S1 to S5 was used as a descriptor of the ensemble performance. Overall, the greatest improvements (MMM approaching the mean of observations) were achieved at S3 or higher calibration stages. For instance, grassland GPP MMM was equal to 1632 g C m−2 yr-1 (S5) while the observed mean was equal to 1763 m-2 yr-1 (average for two sites). Nash-Sutcliffe modelling efficiency coefficients indicated that MMM outperformed individual models in 92.3 % of cases. Our study suggests a cautious use of large-scale, multi-model ensembles to estimate C fluxes in agricultural sites if some site-specific plant and soil observations are available for model calibration. The further development of crop/grassland ensemble modelling will hinge upon the interpretation of results in light of the way models represent the processes underlying C fluxes in complex agricultural systems (grassland and crop rotations including fallow periods).

Abstract: The European Union Effort Sharing Regulation (ESR) will require a 30% reduction in greenhouse gas (GHG) emissions by 2030 compared with 2005 from the sectors not included in the European Emissions Trading Scheme, including agriculture. This will require the estimation of current and future emissions from agriculture, including dairy cattle production systems. Using a farm-scale model as part of a Tier 3 method for farm to national scales provides a more holistic and informative approach than IPCC (2006) Tier 2 but requires independent quality control. Comparing the results of using models to simulate a range of scenarios that explore an appropriate range of biophysical and management situations can support this process by providing a framework for placing model results in context. To assess the variation between models and the process of understanding differences, estimates of GHG emissions from four farm-scale models (DailyWise, FarmAC, HolosNor and SFARMMOD) were calculated for eight dairy farming scenarios within a factorial design consisting of two climates (cool/dry and warm/wet) x two soil types (sandy and clayey) x two feeding systems (grass only and grass/maize). The milk yield per cow, follower cow ratio, manure management system, nitrogen (N) fertilisation and land area were standardised for all scenarios in order to associate the differences in the results with the model structure and function. Potential yield and application of available N in fertiliser and manure were specified separately for grass and maize. Significant differences between models were found in GHG emissions at the farm-scale and for most contributory sources, although there was no difference in the ranking of source magnitudes. The farm-scale GHG emissions, averaged over the four models, was 10.6 t carbon dioxide equivalents (CO(2)e)/ha per year, with a range of 1.9 t CO(2)e/ha per year. Even though key production characteristics were specified in the scenarios, there were still significant differences between models in the annual milk production per ha and the amounts of N fertiliser and concentrate feed imported. This was because the models differed in their description of biophysical responses and feedback mechanisms, and in the extent to which management functions were internalised. We conclude that comparing the results of different farm-scale models when applied to a range of scenarios would build confidence in their use in achieving ESR targets, justifying further investment in the development of a wider range of scenarios and software tools.

Abstract: Impaired animal health causes both productivity and profitability losses on dairy farms, resulting in inefficient use of inputs and increase in greenhouse gas (GHG) emissions produced per unit of product (i.e. emissions intensity). Here, we used subclinical mastitis as an exemplar to benchmark alternative scenarios against an economic optimum and adjusted herd structure to estimate the GHG emissions intensity associated with varying levels of disease. Five levels of somatic cell count (SCC) classes were considered namely 50,000 (i.e. SCC50), 200,000, 400,000, 600,000 and 800,000 cells/mL (milliliter) of milk. The effects of varying levels of SCC on milk yield reduction and consequential milk price penalties were used in a dynamic programming (DP) model that maximizes the profit per cow, represented as expected net present value, by choosing optimal animal replacement rates. The GHG emissions intensities associated with different levels of SCC were then computed using a farm-scale model (HolosNor). The total culling rates of both primiparous (PP) and multiparous (MP) cows for the five levels of SCC scenarios estimated by the model varied from a minimum of 30.9% to a maximum of 43.7%. The expected profit was the highest for cows with SCC200 due to declining margin over feed, which influenced the DP model to cull and replace more animals and generate higher profit under this scenario compared to SCC50. The GHG emission intensities for the PP and MP cows with SCC50 were 1.01 kg (kilogram) and 0.95 kg carbon dioxide equivalents (CO2e) per kg fat and protein corrected milk (FPCM), respectively, with the lowest emissions being achieved in SCC50. Our results show that there is a potential to reduce the farm GHG emissions intensity by 3.7% if the milk production was improved through reducing the level of SCC to 50,000 cells/mL in relation to SCC level 800,000 cells/mL. It was concluded that preventing and/or controlling subclinical mastitis consequently reduces the GHG emissions per unit of product on farm that results in improved profits for the farmers through reductions in milk losses, optimum culling rate and reduced feed and other variable costs. We suggest that further studies exploring the impact of a combination of diseases on emissions intensity are warranted.