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.

Abstract: Modelling N transformations within cropping systems is crucial for N management optimization in order to increase N use efficiency and reduce N losses. Such modelling remains challenging because of the complexity of N cycling in soil–plant systems. In the current study, the uncertainties of six widely used process-based models (PBMs), including APSIM, CROPSYST, DAISY, FASSET, HERMES and STICS, were tested in simulating different N managements (catch crops (CC) and different N fertilizer rates) in 12-year rotations in Western Europe. Winter wheat, sugar beet and pea were the main crops, and radish was the main CC in the tested systems. Our results showed that PBMs simulated yield, aboveground biomass, N export and N uptake well with low RMSE values, except for sugar beet, which was generally less well parameterized. Moreover, PBMs provided more accurate crop simulations (i.e. N export and N uptake) compared to simulations of soil (N mineralization and soil mineral N (SMN)) and environmental variables (N leaching). The use of multi-model ensemble mean or median of four PBMs significantly reduced the mean absolute percentage error (MAPE) between simulations and observations to less than 15% for yield, aboveground biomass, N export and N uptake. Multi-model ensemble also significantly reduced the MAPE for net N mineralization and annual N leaching to around 15%, while it was larger than 20% for SMN. Generally, PBMs well simulated the CC effects on N fluxes, i.e. increasing N mineralization and reducing N leaching in both short-term and long-term, and all PBMs correctly predicted the effects of the reduced N rate on all measured variables in the study. The uncertainties of multi-model ensemble for N mineralization, SMN and N leaching were larger, mainly because these variables are influenced by plant-soil interactions and subject to cumulative long-term effects in crop rotations, which makes them more difficult to simulate. Large differences existed between individual PBMs due to the differences in formalisms for describing N processes in soil–plant systems, the skills of modelers and the model calibration level. In addition, the model performance also depended on the simulated variables, for instance, HERMES and FASSET performed better for yield and crop biomass, APSIM, DAISY and STICS performed better for N export and N uptake, STICS provided best simulation for SMN and N leaching among the six individual PBMs in the study, but all PBMs met difficulties to well predict either average or variance of soil N mineralization. Our results showed that better calibration for soil N variables is needed to improve model predictions of N cycling in order to optimize N management in crop rotations.