Jukuri, open repository of the Natural Resources Institute Finland (Luke) All material supplied via Jukuri is protected by copyright and other intellectual property rights. Duplication or sale, in electronic or print form, of any part of the repository collections is prohibited. Making electronic or print copies of the material is permitted only for your own personal use or for educational purposes. For other purposes, this article may be used in accordance with the publisher’s terms. There may be differences between this version and the publisher’s version. You are advised to cite the publisher’s version. This is an electronic reprint of the original article. This reprint may differ from the original in pagination and typographic detail. Author(s): A. Bougouin, A. Hristov, J. Dijkstra, M.J. Aguerre, S. Ahvenjärvi, C. Arndt, A. Bannink, A.R. Bayat, C. Benchaar, T. Boland, W.E. Brown, L.A. Crompton, F. Dehareng, I. Dufrasne, M. Eugène, E. Froidmont, S. van Gastelen, P.C. Garnsworthy, A. Halmemies-Beauchet-Filleau, S. Herremans, P. Huhtanen, M. Johansen, A. Kidane, M. Kreuzer, B. Kuhla, F. Lessire, P. Lund, E.M.K. Minnée, C. Muñoz, M. Niu, P. Nozière, D. Pacheco, E. Prestløkken, C.K. Reynolds, A. Schwarm, J.W. Spek, M. Terranova, A. Vanhatalo, M.A. Wattiaux, M.R. Weisbjerg, D.R. Yáñez-Ruiz, Z. Yu, and E. Kebreab Title: Prediction of nitrogen excretion from data on dairy cows fed a wide range of diets compiled in an intercontinental database: A meta-analysis Year: 2022 Version: Published version Copyright: The Author(s) 2022 Rights: CC BY 4.0 Rights url: http://creativecommons.org/licenses/by/4.0/ Jukuri, open repository of the Natural Resources Institute Finland (Luke) All material supplied via Jukuri is protected by copyright and other intellectual property rights. Duplication or sale, in electronic or print form, of any part of the repository collections is prohibited. Making electronic or print copies of the material is permitted only for your own personal use or for educational purposes. For other purposes, this article may be used in accordance with the publisher’s terms. There may be differences between this version and the publisher’s version. You are advised to cite the publisher’s version. Please cite the original version: A. Bougouin, A. Hristov, J. Dijkstra, M.J. Aguerre, S. Ahvenjärvi, C. Arndt, A. Bannink, A.R. Bayat, C. Benchaar, T. Boland, W.E. Brown, L.A. Crompton, F. Dehareng, I. Dufrasne, M. Eugène, E. Froidmont, S. van Gastelen, P.C. Garnsworthy, A. Halmemies-Beauchet-Filleau, S. Herremans, P. Huhtanen, M. Johansen, A. Kidane, M. Kreuzer, B. Kuhla, F. Lessire, P. Lund, E.M.K. Minnée, C. Muñoz, M. Niu, P. Nozière, D. Pacheco, E. Prestløkken, C.K. Reynolds, A. Schwarm, J.W. Spek, M. Terranova, A. Vanhatalo, M.A. Wattiaux, M.R. Weisbjerg, D.R. Yáñez-Ruiz, Z. Yu, E. Kebreab (2022). Prediction of nitrogen excretion from data on dairy cows fed a wide range of diets compiled in an intercontinental database: A meta-analysis. Journal of Dairy Science 105(9): 7462-7481. https://doi.org/10.3168/jds.2021-20885. 7462 ABSTRACT Manure nitrogen (N) from cattle contributes to ni- trous oxide and ammonia emissions and nitrate leach- ing. Measurement of manure N outputs on dairy farms is laborious, expensive, and impractical at large scales; therefore, models are needed to predict N excreted in urine and feces. Building robust prediction models requires extensive data from animals under different management systems worldwide. Thus, the study ob- jectives were (1) to collate an international database of N excretion in feces and urine based on individual lac- tating dairy cow data from different continents; (2) to determine the suitability of key variables for predicting fecal, urinary, and total manure N excretion; and (3) to develop robust and reliable N excretion prediction mod- els based on individual data from lactating dairy cows consuming various diets. A raw data set was created based on 5,483 individual cow observations, with 5,420 fecal N excretion and 3,621 urine N excretion measure- ments collected from 162 in vivo experiments conducted by 22 research institutes mostly located in Europe (n = 14) and North America (n = 5). A sequential ap- proach was taken in developing models with increasing complexity by incrementally adding variables that had a significant individual effect on fecal, urinary, or total Prediction of nitrogen excretion from data on dairy cows fed a wide range of diets compiled in an intercontinental database: A meta-analysis A. Bougouin,1* A. Hristov,2 J. Dijkstra,3 M. J. Aguerre,4 S. Ahvenjärvi,5 C. Arndt,6 A. Bannink,7 A. R. Bayat,5 C. Benchaar,8 T. Boland,9 W. E. Brown,10,11 L. A. Crompton,12 F. Dehareng,13 I. Dufrasne,14 M. Eugène,15 E. Froidmont,13 S. van Gastelen,7 P. C. Garnsworthy,16 A. Halmemies-Beauchet-Filleau,17 S. Herremans,13 P. Huhtanen,18 M. Johansen,19 A. Kidane,20 M. Kreuzer,21 B. Kuhla,22 F. Lessire,14 P. Lund,19 E. M. K. Minnée,23 C. Muñoz,24 M. Niu,1,21 P. Nozière,15 D. Pacheco,25 E. Prestløkken,20 C. K. Reynolds,12 A. Schwarm,20 J. W. Spek,7 M. Terranova,26 A. Vanhatalo,17 M. A. Wattiaux,10 M. R. Weisbjerg,19 D. R. Yáñez-Ruiz,27 Z. Yu,11 and E. Kebreab1 1Department of Animal Science, University of California, Davis 95616 2Department of Animal Science, The Pennsylvania State University, University Park 16803 3Animal Nutrition Group, Wageningen University and Research, 6700 AH Wageningen, the Netherlands 4Department of Animal and Veterinary Sciences, Clemson University, Clemson, SC 29634 5Animal Nutrition, Production Systems, Natural Resources Institute Finland (Luke), FI-31600 Jokioinen, Finland 6Mazingira Centre, International Livestock Research Institute (ILRI), 00100 Nairobi, Kenya 7Wageningen Livestock Research, Wageningen University and Research, 6700 AH Wageningen, the Netherlands 8Sherbrooke Research and Development Centre, Agriculture and Agri-Food Canada, Sherbrooke, Quebec, Canada J1M 0C8 9School of Agriculture and Food Science, University College Dublin, Belfield, Dublin 4, Ireland 10Department of Animal and Dairy Sciences, University of Wisconsin–Madison 53706-1205 11Department of Animal Sciences, The Ohio State University, Columbus 43210 12School of Agriculture, Policy and Development, University of Reading, Reading RG6 6AR, United Kingdom 13Department of Valorisation of Agricultural Products, Walloon Agricultural Research Centre, 5030 Gembloux, Belgium 14Department of Veterinary Management of Animal Resources, Faculty of Veterinary Medicine, Fundamental and Applied Research for Animal and Health (FARAH), University of Liège, 4000 Liège, Belgium 15INRAE - Université Clermont Auvergne - VetAgroSup UMR 1213 Unité Mixte de Recherche sur les Herbivores, Centre de recherche Auvergne-Rhône-Alpes, Theix, 63122 Saint-Genès-Champanelle, France 16School of Biosciences, University of Nottingham, Loughborough LE12 5RD, United Kingdom 17Faculty of Agriculture and Forestry, Department of Agricultural Sciences, University of Helsinki, 00014 Helsinki, Finland 18Department of Agricultural Science for Northern Sweden, Swedish University of Agricultural Sciences, SE-901 83, Umeå, Sweden 19Department of Animal Science, Aarhus University, AU Foulum, Dk-8830 Tjele, Denmark 20Department of Animal and Aquacultural Sciences, Norwegian University of Life Sciences, 1433 Ås, Norway 21Institute of Agricultural Sciences, ETH Zurich, 8092 Zurich, Switzerland 22Institute for Farm Animal Biology (FBN), Institute of Nutritional Physiology “Oskar Kellner,” Dummerstorf, Mecklenburg-Vorpommern, Germany 23DairyNZ Ltd., Private Bag 3221, Hamilton, New Zealand 3240 24Instituto de Investigaciones Agropecuarias, INIA Remehue, Ruta 5 S, Osorno, Chile 25Ag Research, Palmerston North 4410, New Zealand 26AgroVet-Strickhof, ETH Zurich, 8315 Lindau, Switzerland 27Estación Experimental del Zaidin, CSIC, 1, 18008 Granada, Spain J. Dairy Sci. 105:7462–7481 https://doi.org/10.3168/jds.2021-20885 © 2022, The Authors. Published by Elsevier Inc. and Fass Inc. on behalf of the American Dairy Science Association®. This is an open access article under the CC BY license (http://creativecommons.org/licenses/by/4.0/). Received June 17, 2021. Accepted April 3, 2022. *Corresponding author: aabougouin@​ucdavis​.edu https://orcid.org/0000-0001-5486-3467 https://orcid.org/0000-0002-0884-4203 https://orcid.org/0000-0003-3728-6885 https://orcid.org/0000-0001-6205-8601 https://orcid.org/0000-0001-6977-0933 https://orcid.org/0000-0002-6276-1097 https://orcid.org/0000-0001-9916-3202 https://orcid.org/0000-0002-4894-0662 https://orcid.org/0000-0002-8644-6892 https://orcid.org/0000-0002-7433-130X https://orcid.org/0000-0002-8457-7741 https://orcid.org/0000-0003-3752-4804 https://orcid.org/0000-0002-6733-4334 https://orcid.org/0000-0003-4651-3779 https://orcid.org/0000-0002-2111-0597 https://orcid.org/0000-0002-7950-4193 https://orcid.org/0000-0003-4547-8449 https://orcid.org/0000-0001-5131-3398 https://orcid.org/0000-0001-6901-1400 https://orcid.org/0000-0001-6502-7050 https://orcid.org/0000-0001-7855-7448 https://orcid.org/0000-0002-2274-8939 https://orcid.org/0000-0003-2945-7189 https://orcid.org/0000-0002-9978-1171 https://orcid.org/0000-0002-2032-5502 https://orcid.org/0000-0002-4062-6927 https://orcid.org/0000-0002-9113-4500 https://orcid.org/0000-0003-4080-9325 https://orcid.org/0000-0003-4484-4611 https://orcid.org/0000-0003-1727-8984 https://orcid.org/0000-0002-9307-9197 https://orcid.org/0000-0003-3151-6782 https://orcid.org/0000-0002-4152-1190 https://orcid.org/0000-0002-5750-2111 https://orcid.org/0000-0002-9602-6189 https://orcid.org/0000-0003-4152-8429 https://orcid.org/0000-0003-0288-8237 https://orcid.org/0000-0001-8713-1641 https://orcid.org/0000-0002-6514-9186 https://orcid.org/0000-0003-4397-3905 https://orcid.org/0000-0002-6165-8522 https://orcid.org/0000-0002-0833-1352 mailto:aabougouin@ucdavis.edu 7463 Journal of Dairy Science Vol. 105 No. 9, 2022 manure N excretion. Nitrogen excretion was predicted by fitting linear mixed models including experiment as a random effect. Simple models requiring dry matter intake (DMI) or N intake performed better for predict- ing fecal N excretion than simple models using diet nutrient composition or milk performance parameters. Simple models based on N intake performed better for urinary and total manure N excretion than those based on DMI, but simple models using milk urea N (MUN) and N intake performed even better for urinary N excretion. The full model predicting fecal N excre- tion had similar performance to simple models based on DMI but included several independent variables (DMI, diet crude protein content, diet neutral detergent fiber content, milk protein), depending on the location, and had root mean square prediction errors as a fraction of the observed mean values of 19.1% for intercontinental, 19.8% for European, and 17.7% for North American data sets. Complex total manure N excretion models based on N intake and MUN led to prediction errors of about 13.0% to 14.0%, which were comparable to models based on N intake alone. Intercepts and slopes of variables in optimal prediction equations developed on intercontinental, European, and North American bases differed from each other, and therefore region- specific models are preferred to predict N excretion. In conclusion, region-specific models that include in- formation on DMI or N intake and MUN are required for good prediction of fecal, urinary, and total manure N excretion. In absence of intake data, region-specific complex equations using easily and routinely measured variables to predict fecal, urinary, or total manure N excretion may be used, but these equations have lower performance than equations based on intake. Key words: manure nitrogen excretion, prediction model, dairy cow INTRODUCTION Ruminants are the most important source of milk for humans, providing 852 million tons of milk in 2019, of which 80% was produced by dairy cattle (OECD- FAO, 2020). Increasing demand for dairy products has led to expansion of dairy herds globally. However, the environmental impact of livestock-based food produc- tion is of concern because livestock manure generates nitrous oxide (N2O) and ammonia (NH3) emissions, and nitrate leaching into soil and ground water, con- tributing to air, water, and soil pollution (FAO, 2002). Nitrous oxide is an important greenhouse gas, with 265 times greater global warming potential than CO2 over a 100-yr period (IPCC, 2007). Emissions of N2O can occur from hydrolysis, nitrification, and denitrification of N constituents mainly in urine, in particular urea, purine derivatives, creatine, and creatinine (Dijkstra et al. 2013a). Although urea is quantitatively the most important N compound in urine, other N compounds are important as well in N2O emissions related to urine (Dijkstra et al. 2013a). Farm animals are considered the main contributors to NH3 emissions, contributing, for instance, 50% of NH3 emissions in the United States (Hristov et al., 2011). When urine and feces are mixed, urease present in feces rapidly converts urea excreted in urine into NH3 and NH4 +, resulting in NH3 vola- tilization (Bougouin et al., 2016). Therefore, urine N is the main contributor to NH3 emission from manure (Hristov et al., 2011). Overall, ruminants are relatively inefficient at utiliz- ing feed N, with 60 to 80% of dietary N being excreted mainly through urine outputs (Kebreab et al., 2001). The nondigested N, from endogenous or feed sources, is excreted in feces (Moore et al., 2014), whereas urinary N may only arise from digested and absorbed N, or mobilized body N. Although 50 to over 90% of total urine N is in the form of urea N (Bougouin et al., 2016), diet composition affects levels of various urinary N compounds and consequently susceptibility to losses after excretion (Dijkstra et al., 2013a). It has also been shown that reduced dietary N supply is accompanied by a decline in the proportion of the easily volatile urea in urine N (Dijkstra et al., 2013a). Hence, large varia- tion in urinary N excretion compared with N excretion in feces presents an opportunity to manipulate diets to reduce urinary N excretion. Predictive tools such as empirical or mechanistic models for estimating N excretion in manure are useful, especially because measurement techniques are labori- ous and costly and difficult to apply in the field (Hris- tov et al., 2019). Several models have been developed to predict total N excretion from lactating dairy cows to assess the efficiency of dairy production and to calcu- late national inventories of N2O emissions. For example, IPCC (2006), revised by IPCC (2019), methodologies for predicting N2O emissions from livestock manure using either default (tier 1 or 2) or country-specific estimates (tier 2 or 3) for total manure N excretion are recommended. However, proportions of N excreted in either urine or feces determine the fate of N as source of pollution (Tamminga, 1996). Therefore, models pre- dicting both urine and feces N excretion separately are of interest and could improve N emissions estimation, resulting in a more accurate environmental impact as- sessment. Therefore, the objectives of the present study were (1) to collate an intercontinental database with data on N intake and excretion in feces and urine for individual lactating dairy cows in different locations of the world and fed on a wide variety of diets; (2) to determine the Bougouin et al.: PREDICTION OF NITROGEN EXCRETION ACROSS DIETS Journal of Dairy Science Vol. 105 No. 9, 2022 7464 suitability of key variables for predicting N excretion in feces and urine, as well as total manure N excretion; and (3) to develop simplified but robust and reliable models that may facilitate prediction of N excretion in urine and feces. MATERIALS AND METHODS Databases and Variable Selection This study is an element of the Global Network project and the Feed and Nutrition Network, which is an activity of the Livestock Research Group of the Global Research Alliance for Agricultural Greenhouse Gases (https:​/​/​globalresearchalliance​.org). A database for analysis was created using measurements made on individual animals at research locations from different parts of the world. The following information, all on an individual animal basis, was included: (1) fecal and urinary N excretion measured by either total collection or by marker methods for fecal N, (2) DMI, (3) dietary nutrient contents [CP, NDF, ADF, ether extract (EE), and starch] and their apparent total-tract digestibility, and (4) lactation performance [milk yield (MY) and concentrations of fat, protein (MProt), lactose, and MUN in milk]. Daily intakes of diet nutrients (CP, NDF, ADF, EE, starch) were calculated based on individual DMI and dietary chemical composition within experi- ments. Intake of N was calculated as CP intake/6.25 when not reported by the collaborators. A total of 5,483 observations were compiled, with 5,420 fecal N excretion and 3,621 urine N excretion measurements, from individual lactating dairy cows of different breeds (e.g., Holstein, Ayrshire, Holstein × Jersey). The raw data set included DMI and diet chemical composition, animal characteristics, and lactation performance from the 162 in vivo experiments (see details in Supplemen- tal Table S1, https:​/​/​data​.mendeley​.com/​datasets/​ 2d3ff88t5g/​1; Kebreab, 2022) that were carried out in Europe (Belgium, Denmark, Finland, France, Ger- many, Ireland, Norway, the Netherlands, Sweden, Swit- zerland, and the United Kingdom), in North America (the United States and Canada), in Central and South America (Costa Rica and Chile), and in Oceania (New Zealand). The intercontinental data set includes all data from the locations mentioned. The data set included diets that tested effects of dietary carbohydrate source (26% of all experiments), CP content (23% of all experiments), and supplemental lipids (6% of all experiments), or were categorized as control diets (28% of all experiments). In several experi- ments, additives were tested (e.g., tannins, monensin, enzymes, yeasts, essential oils, or other plant extracts; 17% of all experiments). Measurements of N excretion in feces and in urine were conducted using total feces and urine collection or using digestibility markers [acid- insoluble ash, chromium oxide (Cr2O3), indigestible NDF, titanium dioxide (TiO2), ytterbium (Yb), and yt- terbium (III) chloride (YbCl3)] and creatinine as urine volume marker with spot sampling approaches, or using total manure collection and digestibility markers for fe- cal N excretion and calculating urinary N excretion by difference. Total manure N excretion was obtained by total collection or calculated as the sum of N excreted in feces and urine. Data Preselection for Model Development An exploratory data screening was performed to evaluate the data for completeness, consistency in no- menclature of each variable, and presence of outliers (Pyle, 1999). Data sets were screened for outliers by vi- sually checking boxplots in R (boxplot function; version 0.98.1102, R Foundation for Statistical Computing) and using the interquartile range (IQR) method (Zwillinger and Kokoska, 2000). The IQR method was applied on subsets that included intercontinental data, European data, or North America data separately. A factor of 1.5 for extremes was used in constructing markers to identify outliers, as shown in Equations [i–iii]: IQR = Third Quartile (Q3) − First Quartile (Q1); [i] Lower Fence = Q1 − IQR × 1.5; [ii] Upper Fence = Q3 + IQR × 1.5. [iii] As a result, the revised data set, used for model devel- opment, contained 5,409 observations on fecal N excre- tion (<1% removed from initial data set) and 3,584 observations on urine N excretion (1% removed) from 162 experiments. Measured variables and their sum- mary statistics are given in Table 1 for the interconti- nental data set, such information for North American and European data sets are given in Supplemental Ta- bles S2 and S3 (https:​/​/​data​.mendeley​.com/​datasets/​ 2d3ff88t5g/​1; Kebreab, 2022). To compare model per- formances on the same basis, evaluation data sets were created for fecal, urinary, and total manure N excre- tion and per region, to have complete information on N excretion, DMI, BW, dietary nutrient content (CP and NDF), nutrient intake (N and NDF), and lactation performance (MY, MProt, and MUN). As a result, the evaluation data sets containing intercontinental data for fecal, urinary, and total manure N excretion had Bougouin et al.: PREDICTION OF NITROGEN EXCRETION ACROSS DIETS https://globalresearchalliance.org https://data.mendeley.com/datasets/2d3ff88t5g/1 https://data.mendeley.com/datasets/2d3ff88t5g/1 https://data.mendeley.com/datasets/2d3ff88t5g/1 https://data.mendeley.com/datasets/2d3ff88t5g/1 7465 Journal of Dairy Science Vol. 105 No. 9, 2022 3,445, 2,489, and 2,490 individual observations, respec- tively. Summary statistics are described in Table 2. Statistical Analyses Random-Effects Model Analysis and Model Development. Linear mixed-effects models were con- structed to predict fecal, urinary, and total manure N excretion (g/d) using the final data sets with inter- continental data, European data, or North America data separately. Mixed-effect models were also devel- oped using the international database, excluding data from New Zealand, Costa Rica, and Chile, to test the models’ accuracy, as farming systems in these countries might be different from the others. No differences in accuracy were observed; thus, data from New Zealand, Costa Rica, and Chile were kept in the international database (data not shown). Random-effect meta-anal- ysis approaches (St-Pierre, 2001) were applied, and N excretion was predicted by fitting a mixed-effect model using the lmer (Bates et al., 2015) procedure of the R statistical language (R Core Team, 2017; version 3.3.0). Several models, as explained hereafter, were developed with different categories of independent variables used as fixed effects, and with experiments included as ran- dom effects. The potential predictors among the dif- ferent categories [DMI; dietary nutrient contents (CP, NDF, ADF, EE, and starch), intakes of diet nutrient contents (N, NDF), lactation performance (MY, milk fat concentration, MProt, milk lactose concentration, and MUN)] were individually tested for their effects on fecal, urinary, and total manure N excretion and se- lected when P < 0.10 (Supplemental Tables S5, S6, and S7; https:​/​/​data​.mendeley​.com/​datasets/​2d3ff88t5g/​ 1; Kebreab, 2022). To evaluate the individual effect, mixed-effect models were applied, including each in- dividual variable as a fixed effect and experiment as random effect. The first set of N excretion prediction models be- gan with simple models based on variables that had a significant individual effect on N excretion (P < 0.10). Variables where pairwise Pearson correlations for predictors had an absolute value of |r| ≥ 0.5 were not included simultaneously in model development to avoid multicollinearity, which decreased statistical power and could lead to the exclusion of significant predictor variables during model construction (Graham, 2003). For instance, DMI and milk production had a Pearson correlation coefficient of 0.73 (Supplemental Table S4; https:​/​/​data​.mendeley​.com/​datasets/​2d3ff88t5g/​1; Ke- breab, 2022) and consequently were not used together in the same model. Overall, if DMI was in the equa- Bougouin et al.: PREDICTION OF NITROGEN EXCRETION ACROSS DIETS Table 1. Summary descriptive statistics of variables used for model development1 Input variable n Mean SD Min Max Animal characteristics             DIM (d) 5,416 111 68.0 1 304   BW (kg) 4,892 641 76.6 432 850 Dietary nutrient content (% of DM)           CP 5,142 16.4 1.55 12.4 20.5   NDF 4,924 35.4 5.62 20.4 49.9   ADF 3,373 20.3 2.93 13.3 27.7   EE 5,482 4.2 1.35 1.7 8.5   Starch 5,482 20.7 5.48 7.8 33.8 Nutrient digestibilities (%)           dCP 4,285 67.3 5.93 51.0 83.5   dNDF 3,083 54.4 11.61 21.3 86.3   dADF 1,593 49.9 13.57 17.9 86.8   dEE 844 63.6 11.99 28.9 94.7   dStarch 1,378 97.6 1.51 92.8 100.0 Performance variables           DMI (kg/d) 5,452 21.4 4.19 9.9 33.0   N intake (g/d) 5,219 569.1 130.41 207.4 944.0   MY (kg/d) 5,385 32.1 9.30 7.3 56.9   MFat (%) 5,191 3.92 0.626 2.25 5.51   MProt (%) 4,813 3.26 0.310 2.42 4.07   MLact (%) 5,117 4.70 0.221 4.09 5.31   MUN (mg/dL) 4,350 11.2 4.541 1.5 24.2 Nitrogen excretion (g/d)           Fecal nitrogen 5,409 184.0 50.38 46.8 322.1   Urinary nitrogen 3,621 175.5 66.22 7.0 365.2   Total manure nitrogen 3,629 358.4 96.14 97.0 633.0 1Summary statistics of the intercontinental data set (after outlier removal) used for model development. EE = ether extract; dCP, dNDF, etc. = digestibility of respective nutrients; MY = milk yield; MFat = milk fat; MProt = milk protein; MLact = milk lactose; Min = minimum; Max = maximum. https://data.mendeley.com/datasets/2d3ff88t5g/1 https://data.mendeley.com/datasets/2d3ff88t5g/1 https://data.mendeley.com/datasets/2d3ff88t5g/1 Journal of Dairy Science Vol. 105 No. 9, 2022 7466 tion, then the other variables included did not include intake of specific nutrients (in units of mass per time) but uncorrelated nutrient concentrations (in units of mass per unit of diet DM). The simple N excretion prediction models were developed based on separate variable categories, namely DMI only (DMI_M), N intake only (N intake_M), nutrient intakes (Diet_ intake_M), dietary nutrient contents (Diet_M), and lactation performance variables (Milk_M). For the more complex models, all significant independent variables among different variable categories were used (Full_M). Another complex model was developed us- ing variables that could be easily and routinely mea- sured on field (Full_Field_M) such as dietary com- position, MY, and milk composition variables. Simple models based on MUN (MUN_M), or on MUN and N intake (MUN_Nintake_M) were only developed with urinary and total manure N excretion, because MUN reflects surplus of N to energy in the body, which is probably not related to fecal N. For each category, intercontinental and location-specific models were de- veloped: Europe (Belgium, Denmark, Finland, France, Germany, Ireland, Norway, the Netherland, Sweden, Switzerland, the United Kingdom) and North America (the United States and Canada). No models were de- veloped for the other locations (New Zealand, Costa Rica, and Chile) as there were too few observations recorded, compared with the other regions. However, observations from these 3 countries were included in the intercontinental database. The mixed-effect models development approach used in this study enabled analysis of fixed effects of inde- pendent preselected variables, as well as experiment- specific deviation of the N excretion response, which was taken into account as a random effect. The general mixed-effect model for single and multiple regressions is represented as follows: Y = β0 + β1X1 + β2X2 + . . . + βnXn + Rj + e, [iv] where β0 denotes the fixed effect of intercept; X1 to Xn denote the fixed effects of predictor variables, and β1 to βn are the corresponding slopes; Rj denotes the random effect of the jth experiment (to capture variations such as different regional weather conditions, measurement methods used, research protocols, and more); and e is the within-experiment error. All variables that had a P-value < 0.10 for their individual relationship with N excretions were further used for model development. The Bayesian information criterion (BIC) was com- puted for each model. Models with the lowest BIC were selected as the best models to predict each N excretion Bougouin et al.: PREDICTION OF NITROGEN EXCRETION ACROSS DIETS T ab le 2 . Su m m ar y de sc ri pt iv e st at is ti cs o f va ri ab le s fr om t he i nt er na ti on al d at a se t us ed f or m od el e va lu at io n1 V ar ia bl e Fe ce s N e va lu at io n da ta s et   U ri ne N e va lu at io n da ta s et   T ot al m an ur e N e v a lu at io n da ta s et n M ea n SD M in M ax n M ea n SD M in M ax n M ea n SD M in M ax D ie ta ry n ut ri en t co nc en tr at io n                               C P ( % ) 3, 44 5 16 .5 1. 48 12 .5 20 .5 2, 48 9 16 .4 1. 48 12 .5 20 .5 2, 49 0 16 .4 1. 48 12 .5 20 .5  N D F ( % ) 3, 44 5 35 .1 5. 51 20 .4 49 .8 2, 48 9 35 .2 5. 82 20 .4 49 .8 2, 49 0 34 .5 5. 84 20 .4 3 49 .8 P er fo rm an ce v ar ia bl es                               D M I (k g/ d) 3, 44 5 22 .2 3. 77 9. 9 33 .0 2, 48 9 22 .0 3. 83 9. 9 32 .8 2, 49 0 21 .9 3. 81 9. 9 32 .8  N i nt ak e (g /d ) 3, 44 5 58 4. 8 11 6. 78 22 2. 4 92 7. 2 2, 48 9 57 6. 2 11 7. 29 22 2. 4 89 5. 0 2, 49 0 57 5. 8 11 6. 92 22 2. 4 92 7. 2  M ilk y ie ld ( kg /d ) 3, 44 5 33 .5 8. 63 7. 3 56 .9 2, 48 9 33 .4 8. 25 7. 3 56 .9 2, 49 0 33 .3 8. 23 7. 3 56 .9  M P ro t (% ) 3, 44 5 3. 23 0. 29 9 2. 43 4. 06 2, 48 9 3. 21 0. 29 9 2. 45 4. 06 2, 49 0 3. 21 0. 29 8 2. 45 4. 06  M U N ( m g/ dL ) 3, 44 5 11 .4 4. 46 1. 5 24 .2 2, 48 9 10 .4 4. 11 1. 5 24 .0 2, 49 0 10 .5 4. 12 1. 5 24 .0  N it ro ge n ex cr et io n2 (g /d ) 3, 44 5 19 1. 6 48 .3 3 46 .8 32 2. 1 2, 48 9 18 4. 4 61 .6 0 7. 0 36 5. 2 2, 49 0 19 1. 6 87 .8 2 11 5. 1 63 2. 6 1 M P ro t = m ilk p ro te in ; M in = m in im um ; M ax = m ax im um . 2 N it ro ge n ex cr et io n fo r fe ca l, ur in ar y, a nd t ot al m an ur e N a cc or di ng t o th e da ta s et . 7467 Journal of Dairy Science Vol. 105 No. 9, 2022 response at each level of complexity (see below). The BIC was calculated as follows: n log (SSEp/n) + (log n) p, where p is the number of regression coefficients, n is the sample size, and SSEp is error sum of squares. A model with a smaller BIC is preferred because it reaches a balance between goodness of fit and model complexity. Variance inflation factor analysis was car- ried out to test independence of predictors in complex models, because influential cases can affect the validity and robustness of meta-analysis conclusions (Sutton et al., 2000; Viechtbauer and Cheung, 2010). Thus, residuals were visually inspected for any pattern. Stu- dentized residuals offer the possibility to identify an observation with a high leverage, which arises when an observation influences the regression model to such an extent that the estimated regression function is pulled toward that potential observation (St-Pierre, 2007). To account for this, any studentized residuals |≥3.0| were removed from the data set. Model Evaluation. Prediction accuracy of each N excretion mixed-effect model was evaluated using the revised k-fold cross-validation method (James et al., 2014), based on the evaluation region-specific data set (total of 3,445, 2,489, and 2,490 observations for fecal, urinary, and total manure N excretion in the intercon- tinental evaluation data set, for instance), with folds composed of individual experiments (k = 127, 105, and 105 experiments for fecal, urinary, and total manure N excretion, respectively). Each individual fold was treated as an evaluation set, where the prediction of N excretion of each fold was calculated using the model that was fitted from the remaining folds. In this cross- evaluation method, the predictions of all folds were used to conduct model evaluation metrics as will be described. A combination of model evaluation metrics was used to assess model performance. Root mean square of pre- diction error (RMSPE), expressed as a fraction of the observed mean was calculated, where a smaller RMSPE indicates better model predictive ability. The mean square prediction error was decomposed into mean bias (MB) and slope bias (SB) deviations to identify systematic biases. The MB and SB were calculated as shown in Equations [v] and [vi]: MB= −( )P O 2 ; [v] SB = (Sp − r × So) 2, [vi] where P and O denote the predicted and observed means, Sp denotes the standard deviation (SD) of pre- dicted values, So denotes the SD of observations, and r denotes the Pearson correlation coefficient. The concordance correlation coefficient (CCC; Lin, 1989) was calculated as follows: CCC = r × Cb, [vii], where Cb v v u = + +                       − 1 2 2 1 ; [viii] v = So/Sp; [ix] u P O o p= −( ) ×( )S S 2 , [x] where P, O, So, and Sp are as previously defined, v pro- vides a measure of scale shift, and �u provides a measure of location shift. The CCC evaluates the degree of de- viation between the best-fit line and the identity line (y = x). Therefore, a CCC of a model closer to 1 is an indication of better model performance. When using different data sets to compare the performance of mod- els, the ratio of RMSPE to SD of the data (observed values) can be used, namely RMSPE-observations SD ratio (RSR), as it takes standardized model perfor- mance relative to the variability in observations in dif- ferent data sets (Moriasi et al., 2007). The RSR was calculated as shown in Equation [xi]: RSR RMSPE S .= o [xi] Smaller RSR (<1) indicates better performance given the variability of observations. RESULTS Database Summary statistics of the final data set, which in- cluded data from 162 experiments used for model de- velopment, are presented in Table 1. Dry matter intake per cow ranged from 9.9 to 33.0 kg/d, and milk pro- duction ranged from 7.3 to 56.9 kg/d. On average, N excretion per cow in feces was 184.0 g/d (SD ± 50.38) and in urine was 175.5 g/d (SD ± 66.22), and total manure N excretion was 358.4 g/d (SD ± 96.14). Most of the experimental diets were fed as TMR or forage and concentrate separately, and included grass silage (35.9%, n = 1,969) or corn silage (32.7%, n = 1,792) as the main forage source, or a forage mix (12.4%, n = Bougouin et al.: PREDICTION OF NITROGEN EXCRETION ACROSS DIETS Journal of Dairy Science Vol. 105 No. 9, 2022 7468 680). Only 4.6% (n = 251) included grass hay as the main forage source, and only 3.9% (n = 215) included pasture grass. Some other forages were used as main source, such as cereal, clover, or alfalfa silage, as well as corn stover (5.0%; n = 275). Nutrient contents varied greatly among diets: in terms of dietary DM, CP con- tents varied between 12.4 and 20.5% with SD ± 1.55%, NDF between 20.4 and 49.9% with SD ± 5.62%, EE between 1.7 and 8.5% with SD ± 1.35%, and starch between 7.8 and 33.8% with SD ± 5.48%. Mixed-Effect Models for Fecal N Excretion Models to predict fecal N excretion are presented in Table 3. Daily fecal N excretion was positively corre- lated with individual variables such as DMI, N intake, dietary contents of CP and NDF, MY, and MProt. Negative relationships with fecal N excretion were ob- served with MUN with North American data (Equation [4]; Table 3). Positive linear relationships between N excretion in feces and DMI or N intake were observed with data from the different continents and had similar accuracy of prediction (Equations [1] and [2]; Table 3). None- theless, differences in accuracy were observed between continents. Both DMI_M and N intake_M models led to prediction errors of 17.7 and 17.4%, respectively, when models were developed from North America data. The DMI_M and N intake_M models led to greater prediction errors when based on European data, with 19.9 and 19.6% error, respectively. Similar results were observed for models based on DMI and N intake with intercontinental data (19.0 and 18.8% prediction error, respectively). The greatest RSR and the lowest CCC were observed with DMI_M models based on European data (Table 3) compared with the other locations, whereas the lowest RSR and the greatest CCC were observed with N intake_M models based on intercon- tinental data (Table 3) compared with the other loca- tions. Negligible MB and SB were observed with these models regardless of location. The DMI_M tended to slightly underpredict excretion at the lower end of N excretion and overpredict excretion at the higher end of fecal N excretion with North American and European data (Figure 1). When models were based on nutrient intake data, N intake was selected as the best predic- tor for all the different locations. The simple models Diet_M included dietary CP content in the best model with intercontinental and European data (Equation [3]; Table 3), and led to greater prediction error than mod- els using DMI or N intake (RMSPE increased for these 2 locations, with RMSPE values of 25.5 and 25.3%, respectively). With the North American data, dietary NDF content was chosen as the best predictor and led to a RMSPE of 23.3%. The RSR was greater than 1.0 compared with previously cited models, regardless of the location. In addition, MB was observed with in- tercontinental and European data (5.52 and 8.35%), and SB was observed in the North American model (11.13%). The Milk_M included MY and MProt to predict fecal N excretion for all locations (Equation [4]; Table 3), and MUN was included with North American data. Models with intercontinental and European data were of comparable accuracy, with the prediction er- ror ranging from 23.4% and 24.3%, respectively. The RMSPE was slightly lower with the model based on North American data (20.3%; Equation [4], Table 3). The RSR was <1.0 for all locations, along with negli- gible MB or SB (<3%) except for MB of 4.35% with European data. Variables used in the complex Full_M models (Equa- tion [5]; Table 3) were not the same among continents. Models based on intercontinental and European data both included DMI, dietary CP and NDF, and MProt (only with intercontinental data), whereas only DMI was used in models based on North American data. The model Full_M for North America led to similar prediction performance (RMSPE of 17.7%) compared with the simple model N intake_M for this location (17.4%). The Full_M models using intercontinental and European data led to prediction error (RMSPE of 19.1 and 19.8%, respectively) similar to the DMI_M (19.0 and 19.9% prediction error for intercontinental and Eu- ropean data) but greater than N intake models for in- tercontinental data (18.8% prediction error). Negligible bias was observed with the Full_M models (<3%), but the Full_M models tended to underpredict excretion at the lower end of fecal N excretion and overpredict excretion at the higher end of N excretion with North American and European data (Figure 2). The complex Full_Field_M model led to better pre- diction error and RSR than models based on dietary nutrient content variables, prediction error and RSR similar to models based on lactation performance, but greater prediction error than all the other models for all locations (RMSPE >20%). Variables selected in the Full_Field_M model were not the same among loca- tions. Only MY was used for all locations, whereas dietary NDF was used in models based on intercon- tinental and European data, and MProt was used for intercontinental and North American models (Equation [6]; Table 3). Mixed-Effect Models for Urinary N Excretion Models to predict urinary N excretion are presented in Table 4. In the simple models, N excretion in urine had positive relationships with DMI, intakes of N and Bougouin et al.: PREDICTION OF NITROGEN EXCRETION ACROSS DIETS 7469 Journal of Dairy Science Vol. 105 No. 9, 2022 Bougouin et al.: PREDICTION OF NITROGEN EXCRETION ACROSS DIETS T ab le 3 . P re di ct io n eq ua ti on s of f ec al N e xc re ti on ( g/ d pe r co w ) ac co rd in g to d iff er en t ca te go ri es a nd m od el p er fo rm an ce e va lu at io n M od el d ev el op m en t1 P re di ct io n eq ua ti on s2 M od el p er fo rm an ce 3 E qu at io n nu m be r   C at eg or y   L oc at io n   n R M SP E , % R SR M B SB C C C 1 D M I_ M In te rc on ti ne n t al 2. 7 (± 2. 99 ) + 8 .4 ( ± 0. 11 ) × D M I 3, 44 5 19 .0 0. 75 0. 21 0. 04 0. 60     E ur op e 2. 3 (± 3. 31 ) + 8 .4 ( ± 0. 13 ) × D M I 2, 30 6 19 .9 0. 80 0. 31 0. 16 0. 54     N or th A m er ic a 11 .1 ( ± 7. 19 ) + 8 .4 ( ± 0. 29 ) × D M I 1, 08 9 17 .7 0. 79 0. 00 0. 09 0. 56 2 N i nt ak e_ M In te rc on ti ne nt al 32 .0 ( ± 2. 90 ) + 0 .2 7 (± 0. 00 4) × N i nt ak e 3, 44 5 18 .8 0. 75 0. 20 0. 02 0. 61     E ur op e 35 .4 ( ± 3. 15 ) + 0 .2 5 (± 0. 00 4) × N i nt ak e 2, 30 6 19 .6 0. 79 0. 32 0. 28 0. 57     N or th A m er ic a 24 .5 ( ± 7. 15 ) + 0 .2 9 (± 0. 00 9) × N i nt ak e 1, 08 9 17 .4 0. 78 0. 00 0. 00 0. 56 3 D ie t_ M In te rc on ti ne nt al 13 6. 7 (± 11 .6 0) + 0 .2 5 (± 0. 07 ) × C P 3, 44 5 25 .5 1. 01 5. 52 1. 13 0. 04     E ur op e 10 5. 0 (± 11 .1 0) + 0 .3 4 (± 0. 06 ) × C P 2, 30 6 25 .3 1. 02 8. 35 1. 76 0. 06     N or th A m er ic a 95 .2 ( ± 27 .1 5) + 0 .3 7 (± 0. 08 5) × N D F 1, 08 9 23 .3 1. 04 0. 38 11 .1 3 − 0. 06 4 M ilk _M In te rc on ti ne nt al 28 .1 ( ± 8. 67 ) + 2 .3 ( ± 0. 08 ) × M Y + 2 4. 2 (± 2. 11 ) × M P ro t 3, 44 5 23 .4 0. 93 2. 98 0. 42 0. 25     E ur op e 44 .3 ( ± 10 .2 0) + 2 .2 ( ± 0. 10 ) × M Y + 1 7. 5 (± 2. 51 ) × M P ro t 2, 30 6 24 .3 0. 98 4. 35 0. 23 0. 18     N or th A m er ic a − 14 .0 ( ± 18 .2 9) + 2 .8 ( ± 0. 16 ) × M Y + 4 3. 8 (± 4. 33 ) × M P ro t − 1 .3 ( ± 0. 39 ) × M U N 1, 08 9 20 .3 0. 91 0. 16 0. 05 0. 29 5 Fu ll_ M In te rc on ti ne nt al − 18 8. 5 (± 21 .5 8) + 8 .3 ( ± 0. 21 ) × D M I + 0 .2 8 (± 0. 07 ) × C P + 0 .3 6 (± 0. 03 1) × N D F + 7 .2 ( ± 2. 36 ) × M P ro t 3, 44 5 19 .1 0. 76 0. 29 0. 09 0. 59     E ur op e − 15 4. 1 (± 22 .6 1) + 8 .0 ( ± 0. 29 ) × D M I + 0 .2 8 (± 0. 08 ) × C P + 0 .3 2 (± 0. 04 ) × N D F 2, 30 6 19 .8 0. 79 0. 15 0. 33 0. 56     N or th A m er ic a 5. 3 (± 9. 46 ) + 8 .5 ( ± 0. 32 ) × D M I 1, 08 9 17 .7 0. 79 0. 00 0. 09 0. 56 6 Fu ll_ F ie ld _M In te rc on ti ne nt al − 85 .6 ( ± 21 .3 7) + 1 .9 ( ± 0. 12 ) × M Y + 2 2. 6 (± 0. 04 ) × M P ro t + 0 .3 9 (± 0. 03 9) × N D F 3, 44 5 24 .2 0. 96 3. 25 0. 01 0. 20     E ur op e 25 .6 ( ± 19 .0 9) + 1 .4 ( ± 0. 14 ) × M Y + 0 .2 7 (± 0. 04 4) × N D F 2, 30 6 24 .7 0. 99 4. 51 0. 57 0. 15     N or th A m er ic a − 17 .7 ( ± 22 .7 7) + 2 .3 ( ± 0. 20 ) × M Y + 4 3. 7 (± 5. 49 ) × M P ro t 1, 08 9 20 .4 0. 91 0. 22 0. 01 0. 28 1 M od el c at eg or ie s ar e as f ol lo w s. S im pl e m od el s: D M I on ly ( D M I_ M ), N i nt ak e on ly ( N i nt ak e_ M ), d ie ta ry n ut ri en t co nt en ts v ar ia bl es ( D ie t_ M ), l ac ta ti on p er fo rm an ce v ar ia bl es (M ilk _M ). C om pl ex m od el s: a ll si gn ifi ca nt i nd ep en de nt v ar ia bl es a m on g va ri ab le c at eg or ie s w er e us ed f or t he f ul l m od el s el ec ti on ( Fu ll_ M ); a ll si gn ifi ca nt i nd ep en de nt v ar ia bl es ea si ly a nd r ou ti ne ly m ea su re d on f ie ld ( Fu ll_ F ie ld _M ). D ie t_ In ta ke _M ( di et ar y nu tr ie nt i nt ak e va ri ab le s) w er e id en ti ca l to N _I nt ak e_ M e qu at io ns a nd a re t he re fo re n ot s ho w n. L oc at io ns a re a s fo llo w s. I nt er co nt in en ta l: B el gi um , C an ad a, C hi le , C os ta R ic a, D en m ar k, F in la nd , Fr an ce , G er m an y, I re la nd , th e N et he rl an ds , N ew Z ea la nd , N or w ay , Sw ed en , Sw it ze rl an d, t he U K , an d th e U S. E ur op e: B el gi um , D en m ar k, F in la nd , Fr an ce , G er m an y, I re la nd , N or w ay , th e N et he rl an ds , Sw ed en , Sw it ze rl an d, a nd t he U K . N or th A m er ic a: th e U S an d C an ad a. 2 I n pa re nt he se s: ± SE . D M I (k g/ d) ; C P = d ie ta ry C P c on te nt ( g/ kg D M ); N D F = d ie ta ry N D F c on te nt ( g/ kg D M ); N in ta ke ( g/ d) ; M Y = m ilk y ie ld ( kg /d ); M P ro t = m ilk p ro te in (% ); M U N ( m g/ dL ). 3 n = n um be r of o bs er va ti on s us ed f or m od el e va lu at io n; R M SP E = r oo t m ea n sq ua re p re di ct io n er ro r, e xp re ss ed a s a pe rc en ta ge o f ob se rv ed m ea n da ily N e xc re ti on i n fe ce s; R SR = R M SP E -o bs er va ti on s SD r at io ; M B = m ea n bi as a s a pe rc en ta ge o f m ea n sq ua re p re di ct io n er ro r; S B = s lo pe b ia s as a p er ce nt ag e of m ea n sq ua re p re di ct io n er ro r; C C C = co nc or da nc e co rr el at io n co ef fic ie nt . Journal of Dairy Science Vol. 105 No. 9, 2022 7470 NDF, and dietary CP content, as well as with MY, MProt, and MUN. Only with North American data, a negative relationship was observed between urinary N excretion and dietary NDF (Equation [10]; Table 4). The simple N intake_M models for urine N excre- tion for the different locations, based on N intake, led to RMSPE of 26.1% with European data, 27.2% with intercontinental data, and 27.1% with North American data (Equation [8]; Table 4). The RSR was <1.0 for these models, and negligible MB and SB were observed (<3%). The N intake_M models tended to overpredict at the lower end of urinary N excretion and under- predict at the higher end of urinary N excretion, in particular with North American data (Figure 3). The other simple models, such as DMI_M, Diet_M, Milk_M, and MUN_M, led to lower prediction perfor- mances than N intake_M models regardless of the loca- tion, with RMSPE varying between 27.7 and 33.9%. Models based on dietary CP or both dietary CP and NDF led to MB (>3%; Equation [10]). With North American data, Equation [10] led to 17.2% SB. Models based on milk performance all included MY, MProt, and MUN (Equation [11]; Table 4), and negligible MB or SB were observed (<3%). Models based on MUN only (Equation [12]) had the worst performance of all simple models considered. The simple MUN_Nintake_M models based on MUN and N intake led to lower RMSPE, compared with previously cited simple models, which varied between 25.0% (with North American data) and 26.5% (with intercontinental data). These models led to SB with North American data (>3%, Equation [13]; Table 4). The complex Full_M models were based on N intake and MUN for intercontinental and European data, but were different from the North American data, with dietary NDF being included along with N intake and MUN as predictors. This latter model led to RMSPE of 28.9% and RSR of 0.99, along with SB (>3%; Equation [14]; Table 4), whereas intercontinental and European Full_M models led to 28.0 and 26.8% prediction error, respectively (Equation [14]; Table 4). These 3 models tended to overpredict at the lower end of urinary N excretion and underpredict at the higher end of uri- nary N excretion for North America data in particular (Figure 4). The complex Full_Field_M model (Equation [15]; Table 4) had a lower prediction error than the Full_M model, regardless of the location. The same variables were used for intercontinental and European models (MY, MProt, MUN, and dietary CP), and dietary NDF content was included in the North American model. All these models led to negligible MB and SB (<3%). Bougouin et al.: PREDICTION OF NITROGEN EXCRETION ACROSS DIETS Figure 1. Predicted vs. observed values based on fecal N excretion (g/d per cow) prediction equations with DMI_M (DMI only) for all locations. Gray lines represent the fitted regression line for the relationship between predicted and observed values, and black lines represent the identity line (y = x). 7471 Journal of Dairy Science Vol. 105 No. 9, 2022 Mixed-Effect Models for Total Manure N Excretion Models for predicting total manure N excretion are presented in Table 5. Positive relationships were observed with DMI, N intake, diet CP content, MY, MProt, and MUN. Negative relationships were observed with dietary NDF content only for the simple model. In Full_M, dietary NDF was included with a positive sign. The N intake_M models based on N intake led to the lowest prediction errors among simple models (≤14.0%; Table 5) for intercontinental, European, and North American data, along with the lowest RSR among sim- ple models, and negligible MB and SB (<3%). Overall, the N intake_M model tended to slightly overpredict at the lower end of total manure N excretion and under- predict at the higher end of total manure N excretion (Figure 5). The DMI_M had greater RMSPE (ranging from 14.5 to 17.7%) compared with N intake models for all locations. The other simple models (i.e., Diet_M, Milk_M, and MUN_M) had lower prediction ability than N intake models with larger RMSPE and RSR, along with some models having MB or SB (>3%; except for Milk_M). For models based on dietary parameters, milk per- formance, or MUN only, better results were observed with North American data compared with other loca- tions. For instance, Equation [20], based on MUN, led to an RMSPE of 19.4% with North American data, along with negligible MB and SB, whereas intercon- tinental and European models led to RMSPE of 23.3 and 24.2%, with MB for intercontinental data only. The Diet_Intake_M models all included N intake as a predicting variable. Once N intake is included along with MUN, the predic- tion performance is similar to models based on N intake only. In addition, the N intake and MUN model and the Full_M model (Equations [21] and [22]; Table 5) led to prediction error (RMSPE <14%) similar to interconti- nental, European, and North American data compared with N intake_M. In the complex Full_M models, N intake, MProt, and MUN were used as predictors with intercontinental and North American data, whereas N intake and dietary NDF content were selected with Eu- ropean data. The RSR was <0.70 for these models, and negligible MB was detected (<3%). Overall, the Full_M model tended to overpredict at the lower end of total manure N excretion and underpredict at the higher end of total manure N excretion (Figure 6). The complex Full_Field_M model had RMSPE greater (>17.0%) than the Full_M model and the N Bougouin et al.: PREDICTION OF NITROGEN EXCRETION ACROSS DIETS Figure 2. Predicted vs. observed values based on fecal N excretion (g/d per cow) prediction equations at the highest level of complexity (Full_M) for all locations. Gray lines represent the fitted regression line for the relationship between predicted and observed values, and black lines represent the identity line (y = x). Journal of Dairy Science Vol. 105 No. 9, 2022 7472Bougouin et al.: PREDICTION OF NITROGEN EXCRETION ACROSS DIETS T ab le 4 . P re di ct io n eq ua ti on s of u ri na ry N e xc re ti on ( g/ d pe r co w ) ac co rd in g to d iff er en t ca te go ri es a nd m od el p er fo rm an ce e va lu at io n M od el d ev el op m en t1 P re di ct io n eq ua ti on s2 M od el p er fo rm an ce 3 E qu at io n nu m be r   C at eg or y   L oc at io n   n R M SP E ,% R SR M B , % SB , % C C C 7 D M I_ M In te rc on ti ne nt al 57 .1 ( ± 6. 20 ) + 5 .4 ( ± 0. 22 ) × D M I 2, 48 9 30 .1 0. 90 1. 83 0. 43 0. 30     E ur op e 38 .3 ( ± 7. 46 ) + 6 .2 ( ± 0. 29 ) × D M I 1, 41 2 30 .4 0. 90 0. 85 0. 03 0. 33     N or th A m er ic a 87 .9 ( ± 11 .6 6) + 4 .3 ( ± 0. 37 ) × D M I 94 4 27 .9 0. 96 1. 98 0. 20 0. 17 8 N i nt ak e_ M In te rc on ti ne nt al 32 .9 ( ± 4. 99 ) + 0 .2 5 (± 0. 00 7) × N i nt ak e 2, 48 9 27 .2 0. 81 1. 42 0. 64 0. 48     E ur op e 20 .2 ( ± 5. 65 ) + 0 .2 8 (± 0. 00 8) × N i nt ak e 1, 22 5 26 .1 0. 78 0. 97 0. 07 0. 56     N or th A m er ic a 63 .5 ( ± 10 .7 9) + 0 .2 0 (± 0. 01 3) × N i nt ak e 94 4 27 .1 0. 93 1. 82 0. 72 0. 22 9 D ie t_ In ta ke _M In te rc on ti ne nt al 86 .5 ( ± 8. 16 ) + 1 0. 3 (± 0. 69 ) × N D F i nt ak e 2, 48 9 32 .5 0. 97 3. 61 0. 01 0. 15     E ur op e 26 .0 ( ± 12 .3 0) + 1 4. 3 (± 1. 32 ) × N D F i nt ak e 1, 22 5 31 .8 0. 94 1. 39 0. 25 0. 24     N or th A m er ic a 13 2. 2 (± 10 .2 2) + 7 .7 ( ± 0. 79 ) × N D F i nt ak e 94 4 28 .4 0. 98 1. 77 0. 48 0. 10 10 D ie t_ M In te rc on ti ne nt al − 11 0. 1 (± 13 .0 7) + 1 .7 ( ± 0. 08 ) × C P 2, 48 9 30 .4 0. 91 4. 42 0. 38 0. 31     E ur op e − 97 .7 ( ± 11 .8 6) + 1 .6 ( ± 0. 07 ) × C P 1, 22 5 28 .9 0. 86 4. 15 0. 66 0. 41     N or th A m er ic a − 14 9. 2 (± 60 .5 3) + 3 .0 ( ± 0. 30 ) × C P − 0 .5 0 (± 0. 10 1) × N D F 94 4 32 .1 1. 10 0. 90 17 .1 9 0. 06 11 M ilk _M In te rc on ti ne n t al − 17 .4 ( ± 11 .4 8) + 1 .8 ( ± 0. 12 ) × M Y + 2 0. 9 (± 2. 52 ) × M P ro t + 6 .2 ( ± 0. 26 ) × M U N 2, 48 9 31 .0 0. 93 2. 13 0. 03 0. 26     E ur op e − 31 .9 ( ± 13 .3 6) + 2 .4 ( ± 0. 15 ) × M Y + 2 0. 1 (± 2. 86 ) × M P ro t + 6 .5 ( ± 0. 30 ) × M U N 1, 22 5 31 .5 0. 94 0. 95 0. 40 0. 26     N or th A m er ic a 6. 02 ( ± 23 .3 9) + 1 .1 ( ± 0. 19 ) × M Y + 2 4. 7 (± 5. 43 ) × M P ro t + 5 .4 ( ± 0. 52 ) × M U N 94 4 27 .7 0. 95 2. 70 0. 92 0. 18 12 M U N _M In te rc on ti ne nt al 10 6. 9 (± 5. 91 ) + 6 .1 ( ± 0. 26 ) × M U N 2, 48 9 32 .6 0. 98 3. 29 0. 05 0. 15     E ur op e 99 .4 ( ± 7. 35 ) + 6 .4 ( ± 0. 30 ) × M U N 1, 22 5 33 .9 1. 01 1. 89 2. 73 0. 10     N or th A m er ic a 12 6. 6 (± 9. 67 ) + 5 .3 ( ± 0. 51 ) × M U N 94 4 28 .1 0. 96 2. 47 1. 23 0. 14 13 M U N _ N in ta ke _M In te rc on ti ne nt al − 8. 2 (± 6. 09 ) + 0 .2 5 (± 0. 00 7) × N i nt ak e + 4 .1 ( ± 0. 25 ) × M U N 2, 48 9 26 .5 0. 79 0. 64 0. 01 0. 55     E ur op e − 0. 32 ( ± 5. 82 ) + 0 .2 3 (± 0. 00 7) × N i nt ak e + 4 .2 ( ± 0. 23 ) × M U N 1, 22 5 25 .9 0. 77 0. 34 0. 18 0. 59     N or th A m er ic a 0. 59 ( ± 11 .4 6) + 0 .2 0 (± 0. 01 2) × N i nt ak e + 5 .2 ( ± 0. 46 ) × M U N 94 4 25 .0 0. 86 2. 60 3. 72 0. 37 14 Fu ll_ M In te rc on ti ne nt al − 4. 4 (± 9. 63 ) + 0 .2 2 (± 0. 01 2) × N i nt ak e + 3 .9 ( ± 0. 37 ) × M U N 2, 48 9 28 .0 0. 84 2. 58 1. 01 0. 44     E ur op e − 15 .8 ( ± 17 .5 3) + 0 .2 7 (± 0. 03 2) × N i nt ak e + 2 .0 ( ± 0. 62 ) × M U N 1, 22 5 26 .8 0. 80 2. 17 1. 50 0. 51     N or th A m er ic a 15 1. 6 (± 35 .4 4) + 0 .2 0 (± 0. 01 3) × N i nt ak e + 4 .3 ( ± 0. 50 ) × M U N − 0 .4 5 (± 0. 10 1) × N D F 94 4 28 .9 0. 99 1. 50 3. 67 0. 19 15 Fu ll_ F ie ld _M In te rc on ti ne nt al − 19 2. 8 (± 16 .2 3) + 1 .4 ( ± 0. 07 ) × C P + 1 .4 ( ± 0. 12 ) × M Y + 2 0. 1 (± 2. 79 ) × M P ro t + 3 .5 ( ± 0. 29 ) × M U N 2, 48 9 26 .5 0. 79 0. 64 0. 01 0. 55     E ur op e − 19 1. 8 (± 17 .9 7) + 1 .3 ( ± 0. 08 ) × C P + 1 .9 ( ± 0. 17 ) × M Y + 1 9. 4 (± 3. 40 ) × M P ro t + 3 .0 ( ± 0. 34 ) × M U N 1, 22 5 25 .9 0. 77 0. 34 0. 18 0. 59     N or th A m er ic a − 29 5. 6 (± 67 .2 2) + 1 .1 ( ± 0. 20 ) × M Y + 2 5. 3 (± 5. 47 ) × M P ro t + 3 .9 ( ± 0. 57 ) × M U N + 2 .5 ( ± 0. 32 ) × C P − 0 .3 3 (± 0. 12 1) × N D F 94 4 25 .7 0. 88 2. 42 0. 93 0. 35 1 M od el c at eg or ie s ar e as f ol lo w s. S im pl e m od el s: D M I on ly ( D M I_ M ), N i nt ak e on ly ( N i nt ak e_ M ), n ut ri en t in ta ke s va ri ab le s (D ie t_ in ta ke _M ), d ie ta ry n ut ri en t co nt en ts v ar ia bl es (D ie t_ M ), la ct at io n pe rf or m an ce v ar ia bl es ( M ilk _M ), M U N o nl y (M U N _M ); M U N a nd N in ta ke ( M U N _N in ta ke _M ). C om pl ex m od el s: a ll si gn ifi ca nt in de pe nd en t va ri ab le s am on g va ri ab le c at eg or ie s w er e us ed f or t he f ul l m od el s el ec ti on ( Fu ll_ M ); a ll si gn ifi ca nt i nd ep en de nt v ar ia bl es e as ily a nd r ou ti ne ly m ea su re d on f ie ld ( Fu ll_ F ie ld _M ). L oc at io ns a re a s fo llo w s. I nt er co nt in en ta l: B el gi um , C an ad a, C hi le , C os ta R ic a, D en m ar k, F in la nd , F ra nc e, G er m an y, I re la nd , t he N et he rl an ds , N ew Z ea la nd , N or w ay , S w ed en , S w it ze rl an d, t he U K , an d th e U S. E ur op e: B el gi um , D en m ar k, F in la nd , Fr an ce , G er m an y, I re la nd , N or w ay , th e N et he rl an ds , Sw ed en , Sw it ze rl an d, a nd t he U K . N or th A m er ic a: t he U S an d C an ad a. 2 I n pa re nt he se s: ± SE . D M I (k g/ d) ; C P = d ie ta ry C P c on te nt ( g/ kg D M ); N D F ( g/ kg D M ); N i nt ak e (g /d ); N D F i nt ak e (g /d ); M Y = m ilk y ie ld ( kg /d ); M P ro t = m ilk p ro te in (% ); M U N ( m g/ dL ). 3 n = n um be r of o bs er va ti on s us ed t o ev al ua te e qu at io ns ; R M SP E = r oo t m ea n sq ua re p re di ct io n er ro r, e xp re ss ed a s a pe rc en ta ge o f ob se rv ed m ea n da ily N e xc re ti on i n ur in e; R SR = R M SP E -o bs er va ti on s SD r at io ; M B = m ea n bi as a s a pe rc en ta ge o f m ea n sq ua re p re di ct io n er ro r; S B = s lo pe b ia s as a p er ce nt ag e of m ea n sq ua re p re di ct io n er ro r; C C C = c on co rd an ce c or re la ti on c oe ff ic ie nt . 7473 Journal of Dairy Science Vol. 105 No. 9, 2022 intake_M model and similar to the models based on DMI, except with North American data. Milk yield, MProt, and dietary CP were used as predictors in Full_Field_M models for all locations, whereas MUN was used in the intercontinental and North American models (Equation [23]; Table 5). Slope bias was present for European models (3.72%), whereas negligible bias was observed for the other locations (<3.0%). DISCUSSION In this study, we identified key predictor variables for fecal, urinary, and total manure N excretion in lactating dairy cows, and we evaluated the prediction accuracy of the developed models. Many different prediction models have already been published for beef and dairy cattle (Nennich et al., 2005; Spek et al., 2013a; Reed et al., 2015). However, in the current evaluation we used a larger database, encompassing many areas of the world, with individual measurements. In addition, the compilation of 162 experiments contained a large vari- ety of diets from experiments conducted primarily with Holstein cows (91.0% of all studies) and across North America, Europe, and these continents combined with Central and South America (Costa Rica, and Chile) and Oceania (New Zealand). The diets encompassed a large variety of feeds, but the major ingredients were corn silage and grass silage. Therefore, the results from this meta-analysis will be most applicable to Holstein cows fed corn silage- or grass silage-based diets. Key Predictors of N Excretion Excretion of N in feces, urine, and total manure was positively related to both DMI and N intake. The DMI is related to endogenous N and undigested microbial N, which end up in feces, leading to logical prediction ability of DMI for fecal N excretion. Increasing DMI increases N excretion simply because greater feed con- sumption means greater N consumption. However, N intake does a better job than DMI in explaining urinary N excretion and therefore total manure N excretion. It makes sense that urinary N is much more affected by N intake, as any surplus of N absorbed that cannot be used for milk N or body N accretion, will have to leave via urine. Reed et al. (2015) developed prediction equa- tions of N excretion in feces and urine in dairy cows from data collected over 30 years during USDA energy metabolism studies (United States). They also noticed a better ability of N intake to predict urinary and total Bougouin et al.: PREDICTION OF NITROGEN EXCRETION ACROSS DIETS Figure 3. Predicted vs. observed values based on urinary N excretion (g/d per cow) prediction equations with N intake_M (N intake only) for all locations. Gray lines represent the fitted regression line for the relationship between predicted and observed values, and black lines rep- resent the identity line (y = x). Journal of Dairy Science Vol. 105 No. 9, 2022 7474 manure N excretion. In addition, they showed that, for fecal N excretion, DMI was better suited than N intake. In the current study, using DMI or N intake led to very similar results in fecal N prediction. Huhtanen et al. (2008) also found DMI to be a better predictor of fe- cal N excretion than N intake, but prediction accuracy improved when both were included together. The data set used by Reed et al. (2015) was not as broad as the one used in the present study, but predic- tor coefficients were similar between studies, even with the intercontinental models. Nonetheless, prediction performance of fecal N excretion based on DMI was better in Reed et al. (2015), with prediction error of 12.6%, than in the current study (17.7% with North American data). The complex fecal N excretion model developed in Reed et al. (2015) used different variables (DMI, dietary ME, NDF and CP contents, and DIM) and led to lower RMSPE (11.4%) than in the current study (Full_M led to 17.7% with North American data). The prediction abilities were lower in the current study for total manure N excretion models than in the study by Reed et al. (2015). When N intake was used to predict urinary excretion, it led to similar prediction error, with RMSPE of 23.7% in Reed et al. (2015) and 27.1% in the current study with North American data. Endogenous N and protein synthesis from ruminal microorganisms and subsequent undigested microbial N are important sources of fecal N, as explained be- fore. In addition, undigested feed protein from forage or concentrates could also contribute to N excretion in feces. Greater DMI and N intake increase rumen microbial growth and N absorption across the diges- tive tract, but greater DMI also increases the digesta passage rate, resulting in increased fecal DM and N excretion in feces and decreased CP digestibility (Di- jkstra et al., 2013b). However, CP digestibility might have a minor contribution compared with DMI, which is directly linked to metabolic and endogenous N and subsequent N losses in feces (Huhtanen et al., 2008). Nevertheless, variation in dietary N supply particularly affects urinary N output (Huhtanen et al., 2008), and nitrogen intake has been shown to be the main driver of N losses in dairy cows (Kebreab et al., 2010), milk N efficiency (Huhtanen and Hristov, 2009), and NH3 emis- sions from manure (Hristov et al., 2011). Increasing N intake increases the amount of N excreted in urine as urea N in particular, which may lead to higher NH3 emissions from manure (Weiss et al., 2009). Colmenero and Broderick (2006) showed that any increase in N intake with diets containing more than 16.5% CP was Bougouin et al.: PREDICTION OF NITROGEN EXCRETION ACROSS DIETS Figure 4. Predicted vs. observed values based on urinary N excretion (g/d per cow) prediction equations at the highest level of complexity (Full_M) for all locations. Gray lines represent the fitted regression line for the relationship between predicted and observed values, and black lines represent the identity line (y = x). 7475 Journal of Dairy Science Vol. 105 No. 9, 2022 Bougouin et al.: PREDICTION OF NITROGEN EXCRETION ACROSS DIETS T ab le 5 . P re di ct io n eq ua ti on s of t ot al N e xc re ti on ( g/ d pe r co w ) ac co rd in g to d iff er en t ca te go ri es a nd m od el p er fo rm an ce e va lu at io n M od el d ev el op m en t1 P re di ct io n eq ua ti on s2 M od el p er fo rm an ce 3 E qu at io n nu m be r   C at eg or y   L oc at io n   n R M SP E , % R SR M B , % SB , % C C C 16 D M I_ M In te rc on ti ne n t al 54 .5 ( ± 7. 06 ) + 1 4. 4 (± 0. 26 ) × D M I 2, 49 0 16 .4 0. 71 1. 04 1. 08 0. 65     E ur op e 40 .1 ( ± 8. 56 ) + 1 4. 9 (± 0. 33 ) × D M I 1, 42 2 17 .7 0. 75 0. 62 0. 70 0. 59     N or th A m er ic a 84 .2 ( ± 12 .9 4) + 1 3. 5 (± 0. 44 ) × D M I 94 0 14 .5 0. 76 0. 52 0. 36 0. 59 17 N i nt ak e_ M In te rc on ti ne nt al 49 .0 ( ± 5. 08 ) + 0 .5 6 (± 0. 00 7) × N i nt ak e 2, 49 0 13 .4 0. 58 1. 16 1. 64 0. 79     E ur op e 43 .4 ( ± 5. 81 ) + 0 .5 7 (± 0. 00 8) × N i nt ak e 1, 42 2 13 .6 0. 57 1. 54 2. 14 0. 79     N or th A m er ic a 68 .8 ( ± 11 .2 3) + 0 .5 3 (± 0. 01 5) × N i nt ak e 94 0 13 .3 0. 69 0. 34 2. 05 0. 66 18 D ie t_ M In te rc on ti ne nt al 40 .1 ( ± 24 .6 1) + 2 .2 ( ± 0. 11 ) × C P − 0 .1 3 (± 0. 03 8) × N D F 2, 49 0 20 .8 0. 90 5. 31 1. 30 0. 33     E ur op e − 29 .0 ( ± 17 .2 4) + 2 .2 ( ± 0. 10 ) × C P 1, 42 2 20 .3 0. 86 4. 92 1. 65 0. 40     N or th A m er ic a − 80 .5 ( ± 74 .7 5) + 3 .0 ( ± 0. 44 ) × C P 94 0 19 .6 1. 02 0. 19 4. 49 0. 04 19 M ilk _M In te rc on ti ne nt al − 41 .0 ( ± 15 .9 2) + 5 .1 ( ± 0. 16 ) × M Y + 6 .0 ( ± 0. 35 ) × M U N + 5 4. 5 (± 3. 60 ) × M P ro t 2, 49 0 19 .7 0. 85 1. 90 1. 09 0. 41     E ur op e − 42 .1 ( ± 17 .9 2) + 5 .9 ( ± 0. 20 ) × M Y + 6 .7 ( ± 0. 39 ) × M U N + 4 5. 1 (± 3. 91 ) × M P ro t 1, 42 2 20 .6 0. 87 0. 70 0. 00 0. 40     N or th A m er ic a − 71 .8 ( ± 33 .6 8) + 4 .2 ( ± 0. 29 ) × M Y + 4 .1 ( ± 0. 74 ) × M U N + 8 5. 5 (± 7. 98 ) × M P ro t 94 0 17 .9 0. 93 0. 77 0. 00 0. 24 20 M U N _M In te rc on ti ne nt al 29 0. 7 (± 8. 46 ) + 6 .1 ( ± 0. 38 ) × M U N 2, 49 0 23 .3 1. 00 4. 35 0. 03 0. 07     E ur op e 26 7. 3 (± 10 .5 1) + 7 .1 ( ± 0. 47 ) × M U N 1, 42 2 24 .2 1. 02 2. 57 2. 42 0. 03     N or th A m er ic a 36 2. 1 (± 13 .4 2) + 3 .7 ( ± 0. 79 ) × M U N 94 0 19 .4 1. 01 0. 56 1. 22 0. 00 21 M U N _ N in ta ke _M In te rc on ti ne nt al 30 .0 ( ± 6. 02 ) + 0 .5 5 (± 0. 00 8) × N i nt ak e + 2 .5 ( ± 0. 25 ) × M U N 2, 49 0 13 .4 0. 58 0. 74 0. 63 0. 79     E ur op e 31 .2 ( ± 6. 84 ) + 0 .5 5 (± 0. 00 9) × N i nt ak e + 2 .1 ( ± 0. 27 ) × M U N 1, 42 2 13 .7 0. 58 1. 04 1. 16 0. 79     N or th A m er ic a 33 .1 ( ± 12 .9 4) + 0 .5 3 (± 0. 15 ) × N i nt ak e + 3 .2 ( ± 0. 55 ) × M U N 94 0 13 .0 0. 68 0. 50 2. 44 0. 67 22 Fu ll_ M In te rc on ti ne nt al − 24 .5 ( ± 20 .1 3) + 0 .5 4 (± 0. 01 7) × N i nt ak e + 2 .5 ( ± 0. 48 ) × M U N + 1 4. 8 (± 5. 00 ) × M P ro t 2, 49 0 13 .4 0. 58 0. 79 0. 69 0. 79     E ur op e − 53 .1 ( ± 41 .5 2) + 0 .4 5 (± 0. 03 4) × N i nt ak e + 0 .4 ( ± 0. 11 ) × N D F 1, 42 2 13 .5 0. 57 1. 50 2. 34 0. 79     N or th A m er ic a − 56 .1 ( ± 25 .3 8) + 0 .5 4 (± 0. 02 1) × N i nt ak e + 2 .6 ( ± 0. 68 ) × M U N + 2 6. 4 (± 6. 42 ) × M P ro t 94 0 12 .9 0. 67 0. 49 2. 53 0. 68 23 Fu ll_ F ie ld _M In te rc on ti ne nt al − 26 9. 4 (± 29 .9 8) + 4 .5 ( ± 0. 22 ) × M Y + 1 .9 ( ± 0. 51 ) × M U N + 6 3. 0 (± 5. 09 ) × M P ro t + 1 .6 ( ± 0. 13 ) × C P 2, 49 0 17 .6 0. 76 2. 47 2. 60 0. 56     E ur op e − 19 2. 9 (± 35 .1 2) + 4 .5 ( ± 0. 32 ) × M Y + 4 5. 1 (± 6. 96 ) × M P ro t + 1 .5 ( ± 0. 13 ) × C P 1, 42 2 17 .2 0. 73 2. 53 3. 72 0. 60     N or th A m er ic a − 49 7. 6 (± 79 .6 1) + 4 .1 ( ± 0. 32 ) × M Y + 3 .2 ( ± 0. 85 ) × M U N + 8 8. 0 (± 8. 73 ) × M P ro t + 2 .6 ( ± 0. 43 ) × C P 94 0 17 .9 0. 93 0. 37 0. 88 0. 28 1 M od el c at eg or ie s ar e as f ol lo w s. S im pl e m od el s: D M I on ly ( D M I_ M ), N i nt ak e on ly ( N i nt ak e_ M ), n ut ri en t in ta ke s va ri ab le s (D ie t_ in ta ke _M ), d ie ta ry n ut ri en t co nt en ts v ar ia bl es (D ie t_ M ), la ct at io n pe rf or m an ce v ar ia bl es ( M ilk _M ), M U N o nl y (M U N _M ), M U N a nd N in ta ke ( M U N _N in ta ke _M ). C om pl ex m od el s: a ll si gn ifi ca nt in de pe nd en t va ri ab le s am on g va ri ab le c at eg or ie s w er e us ed f or t he f ul l m od el s el ec ti on ( Fu ll_ M ); a ll si gn ifi ca nt i nd ep en de nt v ar ia bl es e as ily a nd r ou ti ne ly m ea su re d on f ie ld ( Fu ll_ F ie ld _M ). M od el s ba se d on in ta ke v ar ia bl es l ed t o m od el s us in g N i nt ak e, a s fo r E qu at io n [1 7] . L oc at io ns a re a s fo llo w s. I nt er co nt in en ta l: B el gi um , C an ad a, C hi le , C os ta R ic a, D en m ar k, F in la nd , Fr an ce , G er m an y, I re la nd , t he N et he rl an ds , N ew Z ea la nd , N or w ay , S w ed en , S w it ze rl an d, t he U K , a nd t he U S. E ur op e: B el gi um , D en m ar k, F in la nd , F ra nc e, G er m an y, I re la nd , N or w ay , t he N et he rl an ds , Sw ed en , Sw it ze rl an d, a nd t he U K . N or th A m er ic a: t he U S an d C an ad a. 2 I n pa re nt he se s: ± SE . D M I (k g/ d) ; C P ( g/ kg D M ); N D F ( g/ kg D M ); N i nt ak e (g /d ); N D F i nt ak e (g /d ); M Y = m ilk y ie ld ( kg /d ); M P ro t = m ilk p ro te in ( % ); M U N ( m g/ dL ). 3 n = n um be r of o bs er va ti on s us ed t o ev al ua te e qu at io ns ; R M SP E = r oo t m ea n sq ua re p re di ct io n er ro r, e xp re ss ed a s a pe rc en ta ge o f ob se rv ed m ea n to ta l m an ur e da ily N e xc re ti on ; R SR = R M SP E -o bs er va ti on s SD r at io ; M B = m ea n bi as a s a pe rc en ta ge o f r oo t m ea n sq ua re p re di ct io n er ro r; S B = s lo pe b ia s as a p er ce nt ag e of r oo t m ea n sq ua re p re di ct io n er ro r; C C C = c on co rd an ce c or re la ti on c oe ff ic ie nt . Journal of Dairy Science Vol. 105 No. 9, 2022 7476 lost mainly as urinary urea N. There are 2 main sources of urinary urea N: surplus of rumen degradable N (rela- tive to rumen available energy for microbes), leading to NH3-N absorption, and surplus absorbed amino acid-N (relative to what the body can use for milk protein and body accretion N), which is deaminated and gives rise to NH3 production (Spek et al., 2013a). Results of the present study hence agree with previous findings that N intake is the main driver of total manure N excretion. In complex models, along with N intake, Reed et al. (2015) used both dietary contents of NDF and CP as predictors of fecal, urinary, and total N excretion. In our study, dietary CP was used in simple models as well as in fecal Full_M models and urinary and total N excretion Full_field_M models, whereas dietary NDF was not included in fecal Full_M models for North American data. However, urinary Full_Field_M includ- ed dietary NDF content with North American data. Dietary NDF was included in simple fecal and urinary Diet_M models developed with North American data as well as in simple total manure N Diet_M models with intercontinental data. In the simple models, the coefficient sign for NDF was positive when predicting fecal N excretion and negative for urinary N excretion in our study but was positive in Reed et al. (2015), whereas the coefficient sign was positive for dietary CP in both studies. The effect of dietary CP content on N excretion has been well studied (Wattiaux and Karg, 2004; Colmenero and Broderick, 2006; Aguerre et al., 2010), and results generally compare well with the present study. Decreas- ing dietary CP contents from 17.5 to 12.5% reduced urine N excretion by 59% (and by 69% when additional rumen-protected methionine was supplied), total N emissions from slurry by 72%, and N2O emissions from farmyard manure with urine-rich slurry by 78% (Kröber et al., 2000; Külling et al., 2001). The present results are also in line with several studies which indicated that increasing dietary CP content linearly increased N excreted in feces and urine (Colmenero and Broderick, 2006; Weiss et al., 2009; Spek et al., 2013a). It has been shown that any increase in dietary protein or N intake would lead to substantial increases in urinary loss (Van Soest, 1994), with any surplus of N ingested in excess of requirements being excreted in urine (Castillo et al., 2000). Moreover, the true urea N part of urine N decreases more with decrease in CP content than the nonurea N part in urine N (Spek et al., 2013b). It is well known that dietary factors, such as the content of dietary CP (Marini et al., 2004; Huhtanen Bougouin et al.: PREDICTION OF NITROGEN EXCRETION ACROSS DIETS Figure 5. Predicted vs. observed values based on total manure N excretion (g/d per cow) prediction equations with N intake_M models (N intake only) developed for all locations. Gray lines represent the fitted regression line for the relationship between predicted and observed values, and black lines represent the identity line (y = x). 7477 Journal of Dairy Science Vol. 105 No. 9, 2022 and Hristov, 2009) or diet digestibility (Theurer et al., 2002), affect apparent N utilization in ruminants. Ap- parent N utilization is reported to increase and urinary N excretion to decrease with highly digestible diets containing cereal grains such as barley (Cohen et al., 2006). As discussed by Dijkstra et al. (2018), apparent N utilization and excretion are not exclusively driven by dietary protein concentration but also depend on dietary carbohydrate composition and digestible energy supply, with decreased milk N output with high-fiber diets compared with high-starch isoenergetic diets. In their complex models, Reed et al. (2015) observed a positive relationship between dietary NDF and urinary and total N excretion but a negative relationship be- tween dietary NDF and fecal N excretion. However, as mentioned before, the results of our study are different. Indeed, dietary NDF content had a positive effect on fe- cal N excretion in the Diet_M models (North American data only) and Full_M models (intercontinental and European data), and a negative effect on urinary N excretion in simple or Full_M models (North American data only) and on total manure N excretion in simple models (intercontinental data only). These differences may be due to the variability in NDF contents or na- ture within studies. Indeed, the nature of the NDF used in this study is more diverse than in that of Reed et al. (2015), due to the larger and more diverse source of data collected from all around the world, and not only in a specific region in the United States. In addition, we can speculate that N excretion is indirectly linked to dietary NDF, because both are related to DMI, and that presumably greater NDF content will coin- cide with lower DMI, leading to lower N excretion in urine and total manure. Indeed, it has been shown that greater dietary NDF content is usually linked to lower DMI (Arelovich et al., 2008), which in turn is linked to lower N excretion. However, high dietary intakes and high amounts of fiber in the diet are also associated with greater endogenous N losses from the gut (Ouellet et al., 2002), which has a major contribution to fecal N excretion (Huhtanen et al., 2008), and thus might explain the positive relationship between dietary NDF and fecal N excretion. In addition, Castillo et al. (2000) also reported that fiber-based supplements tended to increase fecal N excretion compared with starch-based concentrates. In several models, MProt was positively related to fecal, urinary, and total manure N excretion. Increasing dietary CP content has been reported to increase milk protein yield (Colmenero and Broderick, 2006), and greater level of dietary CP will most likely have greater N excretion, as mentioned before. Milk urea N was also Bougouin et al.: PREDICTION OF NITROGEN EXCRETION ACROSS DIETS Figure 6. Predicted vs. observed values based on total manure N excretion (g/d per cow) prediction equations at the highest level of com- plexity (Full_M) for all locations. Gray lines represent the fitted regression line for the relationship between predicted and observed values, and black lines represent the identity line (y = x). Journal of Dairy Science Vol. 105 No. 9, 2022 7478 positively related to both urinary and total manure N excretion, and MUN was selected in Full_M models for every region, except for the Full_M model to predict total N excretion with European data. A positive linear relationship between N excreted in urine and MUN was found also by Kauffman and St-Pierre (2001). Increas- ing N intake leads to a more pronounced increase in urinary N excretion than in fecal N excretion (Kebreab et al., 2001) and generally results in an elevated con- centration of MUN. In addition, MUN, urine volume, and urinary excretion of total N and urea N all in- creased in response to dietary CP content. Milk urea N concentration has already been suggested as a predictor of N excretion in urine by Jonker et al. (1998), who proposed a simple model of N balance from 4 in vivo studies, gathering 70 observations in dairy cows from North America. Spek et al. (2013b) compared 5 equa- tions to predict urinary N excretion based on MUN, and observed large differences between these equations in predicted urinary N excretion. Spek et al. (2013a) performed a meta-analysis using MUN as a predictor of urinary N excretion with data from Europe and North America demonstrating different relationships for both data sources. In our study, we also noticed different prediction performances between Europe and North America, which could be partly explained by location- specific diets. Simple models including only MUN led to higher pre- dictive errors for urinary N excretion compared with total manure N excretion. However, models based on MUN and predicting urinary N excretion led to 28.1 to 33.9% error. From a theoretical viewpoint, one would expect MUN to correspond better to urinary N excre- tion than to total manure N excretion, as both urine urea N excretion and MUN are a function of blood urea nitrogen concentration. However, not all urine N is urea N, and it has been shown that variation exists in the fraction of urea N in urine N (62 to 86% of urinary N excreted with urea), as discussed via a modeling ap- proach by Dijkstra et al. (2018). The accuracy of MUN as a predictor of urinary N excretion may improve when various factors that affect this relationship can be taken into account (reviewed by Spek et al., 2013b). These factors include BW, water intake, dietary CP content, and frequency of milking and feeding. Not surprisingly, MUN and N intake used simultaneously in urinary and total manure N excretion models led to better predic- tion performance than MUN alone. For urinary N ex- cretion, the MUN and N intake used simultaneously led to better prediction performance than N intake only. For total manure N excretion models, simple equations based on N intake and complex equations based on N intake and MUN led to similar prediction performance. Thus, MUN does not help in explaining more variation. This is in contrast to urinary N excretion, where MUN added to N intake explained some more variation in urinary N excretion. Application of N Excretion Models The simplest models for N excretion in feces and urine based on N intake led to about 18 to 27% predic- tion errors, depending on the location. These models led to better prediction ability for total manure N excretion by reducing the error to about 13 to 14%. Simple models based on dietary composition or milk parameters seemed more adapted for total manure N excretion than for fecal or urinary N excretion. This may be partly related to difficulties in measurement techniques to determine fecal and urinary N excretion (Hristov et al., 2019), where errors in measurement in fecal N excretion may be counterbalanced by errors in urinary N excretion and vice versa. We observed that the accuracy of prediction of fecal, urinary, or total manure N excretion was similar or best with models that included dietary intakes, either DMI or N intake, along with a few other covariates, such as dietary CP, NDF, MProt, or MUN. These other covari- ates may reflect factors that modify the partitioning of N excreted in urine and feces. Using all available vari- able information did not consistently improve predic- tion performance with the full models compared with simpler models and between locations. Models that include DMI or N intake showed better performance than models without DMI or N intake, indicating that these variables are among the main drivers of N losses in either feces or urine. One important reason that N excretion is related to N intake is simply that the varia- tion in N intake is large (200–900 g of N/d). But DMI and N intake also account for variation in energy and protein supply (i.e., DMI, N intake), or ruminal and postruminal protein supply, protein retained by the animal, and protein secreted in milk, which determines the amount of digested and undigested N as well as the fraction of digested N that is not retained by the animal and excreted in feces and urine. Models generally differed between locations, with different coefficients, prediction performance, and, in some cases, variables included in models; thus, models are location-specific. For instance, all simple models to predict fecal N excretion performed better with North American data than with intercontinental or Euro- pean data. For urinary N excretion, N intake models performed better with European data than with inter- continental or North American data. Intercontinental models were developed based on a data set contain- ing a greater proportion of European data compared with North American or other location data (73.3% Bougouin et al.: PREDICTION OF NITROGEN EXCRETION ACROSS DIETS 7479 Journal of Dairy Science Vol. 105 No. 9, 2022 vs. 23.0% and 3.6%, respectively). Thus, logically the greater proportion of European data influenced the results toward better prediction for intercontinental models compared with North American. For different models (Diet_M, Full_M, or Full_Field_M) derived from intercontinental, European, or North American data, different predicting variables were sometime used. Overall, a need exists to use location-specific models when predicting N excretion in feces, urine, or total manure, and the discrepancies are probably due to the different types of diets used among regions of the world. This is not an unexpected result, as previous analyses of enteric methane prediction models in dairy cattle or beef cattle following a similar approach to ours came to similar conclusions (Niu et al., 2018; van Lingen et al., 2019). Dufrasne et al. (2013) also pointed out the im- portance of environmental factors that would directly affect animal productivity and hence N retained by the animal, such as heat stress, disease, BW, and feeding behavior, which are indirectly linked to location. Estimation of DMI and, along with this, N intake is still challenging on farm in practice because it is not routinely monitored. Voluntary DMI prediction equa- tions require individual animal information that is in part routinely available (milk production and composi- tion, DIM), or in part not routinely available (BW; e.g., NRC, 2001; Lahart et al., 2019), or are based on dietary energy content and energy requirements for mainte- nance, milk production, body growth, and pregnancy. Nonetheless, variables not yet routinely available could be obtained with calculations based on feeding systems or from milking robot records, specifically for BW, which is recorded at each milking visit. But even then, it remains to be shown whether intake estimations are accurate enough to provide an advantage over milk parameters based on intake estimations, as milk param- eters are obtained as accurate data. In addition, current research is oriented toward the use of Fourier-transform mid-infrared spectroscopy of milk components, which shows promise for predicting N use efficiency, as shown by Grelet et al. (2020). Milk N, which is used in models developed in this study, is already routinely estimated using mid-infrared spectroscopy with a high degree of accuracy. Thus, the use of full models to predict total manure N excretion using MY, milk protein content, and diet CP content, which could be easily obtained on farm with confinement systems, is preferred only if N intake data are not available. Otherwise, simple models based on MUN and N intake or full models based on N intake, dietary NDF, MUN, and milk protein con- tent would be better suited. For urinary N excretion, in absence of intake data, full models based on diet and milk parameters would be preferred. Simple models based on milk parameters would be chosen for fecal N excretion, but, if intake data are available, then simple models based on N intake would be preferred. Thus, information on intake is required for good prediction of N excretion, in particular of total manure N excretion. In absence of intake data, milk and diet composition parameters may be used, but prediction accuracy of N excretion will decrease. CONCLUSIONS Models based on DMI or N intake can be used to predict fecal and total manure N excretion with good accuracy, and urinary N excretion with satisfactory ac- curacy. Prediction accuracy may be somewhat further improved by adding diet composition or milk parame- ters to intake parameters in complex models. In absence of intake data, models using diet composition and milk performance parameters could be used to predict fecal, urinary, and total manure N excretion, but with greater prediction error and occurrence of MB or SB, or both. Intercepts and slopes of variables in optimal prediction equations developed on intercontinental, European, and North American bases differed from each other, and region-specific models are preferred to predict N excretion. Complex, region-specific models to predict total manure N excretion should be used when inputs are available, whereas simple location-specific models based on DMI or N intake should be used for fecal and urine N excretion prediction. ACKNOWLEDGMENTS This study is part of the Global Network project, funded by the Joint Programming Initiative on Agricul- ture, Food Security and Climate Change (FACCE-JPI) and is an activity of the Feed and Nutrition Network, which is part of the Livestock Research Group of the Global Research Alliance for Agricultural Greenhouse Gases (https:​/​/​globalresearchalliance​.org; accessed May 10, 2021). Partial funding for the study was provided through USDA (Washington, DC) National Institute of Food and Agriculture (NIFA), awards 2014-67003- 21979 and 2019-67019-29400, and Federal Appropria- tions under Project PEN 04539 and Accession Number 1000803 and FONDECYT/Regular Accession Number 1191476. The authors have not stated any conflicts of interest. REFERENCES Aguerre, M. J., M. A. Wattiaux, T. Hunt, and B. R. Larget. 2010. Effect of dietary crude protein on ammonia-N emission measured by herd nitrogen mass balance in a freestall dairy barn managed under farm-like conditions. Animal 4:1390–1400. https:​/​/​doi​.org/​ 10​.1017/​S1751731110000248. Bougouin et al.: PREDICTION OF NITROGEN EXCRETION ACROSS DIETS https://globalresearchalliance.org https://doi.org/10.1017/S1751731110000248 https://doi.org/10.1017/S1751731110000248 Journal of Dairy Science Vol. 105 No. 9, 2022 7480 Arelovich, H. M., C. S. Abney, J. A. Vizcarra, and M. L. Galyean. 2008. Effects of dietary neutral detergent fiber on intakes of dry matter and net energy by dairy and beef cattle: Analysis of published data. Prof. Anim. Sci. 24:375–383. https:​/​/​doi​.org/​10​ .15232/​S1080​-7446(15)30882​-2. Bates, D., M. Mächler, B. Bolker, and S. Walker. 2015. Fitting linear mixed-effects models using lme4. J. Stat. Softw. 67:1–48. https:​/​/​ doi​.org/​10​.18637/​jss​.v067​.i01. Bougouin, A., A. Leytem, J. Dijkstra, R. S. Dungan, and E. Kebreab. 2016. Nutritional and environmental effects on ammonia emissions from dairy cattle housing: A meta-analysis. J. Environ. Qual. 45:1123–1132. https:​/​/​doi​.org/​10​.2134/​jeq2015​.07​.0389. Castillo, A. R., E. Kebreab, D. E. Beever, and J. France. 2000. A re- view of efficiency of nitrogen utilisation in lactating dairy cows and its relationship with environmental pollution. J. Anim. Feed Sci. 9:1–32. https:​/​/​doi​.org/​10​.22358/​jafs/​68025/​2000. Cohen, D. C., C. R. Stockdale, and P. T. Doyle. 2006. Feeding an energy supplement with white clover silage improves rumen fer- mentation, metabolisable protein utilisation, and milk production in dairy cows. Aust. J. Agric. Res. 57:367–375. https:​/​/​doi​.org/​10​ .1071/​AR05166. Colmenero, J. J. O., and G. A. Broderick. 2006. Effect of dietary crude protein concentration on milk production and nitrogen utilization in lactating dairy cows. J. Dairy Sci. 89:1704–1712. https:​/​/​doi​ .org/​10​.3168/​jds​.S0022​-0302(06)72238​-X. Dijkstra, J., A. Bannink, P. M. Bosma, E. A. Lantinga, and J. W. Reijs. 2018. Modelling the effect of nutritional strategies for dairy cows on the composition of excreta nitrogen. Front. Sustain. Food Syst. 2:63. https:​/​/​doi​.org/​10​.3389/​fsufs​.2018​.00063. Dijkstra, J., O. Oenema, J. W. van Groenigen, J. W. Spek, A. M. van Vuuren, and A. Bannink. 2013a. Diet effects on urine composition of cattle and N2O emissions. Animal 7:292–302. https:​/​/​doi​.org/​10​ .1017/​S1751731113000578. Dijkstra, J., C. K. Reynolds, E. Kebreab, A. Bannink, J. L. Ellis, J. France, and A. M. van Vuuren. 2013b. Challenges in ruminant nutrition: Towards minimal nitrogen losses in cattle. Pages 47–58 in Energy and Protein Metabolism and Nutrition in Sustainable Animal Production. J. W. Oltjen, E. Kebreab, and H. Lapierre, ed. Wageningen Academic Publishers. https:​/​/​doi​.org/​10​.3920/​978​-90​ -8686​-781​-3​_3. Dufrasne, I., E. Knapp, L. Istasse, D. Veselko, E. Piraux, V. Robaye, and J.-L. Hornick. 2013. Influence de facteurs environnementaux influençant la teneur en urée dans le lait de vache en Wallonie et estimation des rejets azotés (Influence of environmental fac- tors influencing the urea content in cow’s milk in Wallonia and estimation of nitrogen losses). Biotechnol. Agron. Soc. Environ. 17:251–258. FAO. 2002. World Agriculture: Towards 2015/2030. Summary Report. FAO. Accessed May 22, 2021. ftp:​/​/​ftp​.fao​.org/​docrep/​fao/​004/​ y3557e/​y3557e​.pdf. Graham, M. H. 2003. Confronting multicollinearity in ecological mul- tiple regression. Ecology 84:2809–2815. https:​/​/​doi​.org/​10​.1890/​ 02​-3114. Grelet, C., E. Froidmont, L. Foldager, M. Salavati, M. Hostens, C. P. Ferris, K. L. Ingvartsen, M. A. Crowe, M. T. Sorensen, J. A. Fernandez Pierna, A. Vanlierde, N. Gengler, and F. Dehareng. 2020. Potential of milk mid-infrared spectra to predict nitrogen use efficiency of individual dairy cows in early lactation. J. Dairy Sci. 103:4435–4445. https:​/​/​doi​.org/​10​.3168/​jds​.2019​-17910. Hristov, A. N., A. Bannink, L. A. Crompton, P. Huhtanen, M. Kreu- zer, M. McGee, P. Noziere, C. K. Reynolds, A. Bayat, D. Yáñez- Ruiz, J. Dijkstra, E. Kebreab, A. Schwarm, K. J. Shingfield, and Z. Yu. 2019. Invited review: Nitrogen in ruminant nutrition: A review of measurement techniques. J. Dairy Sci. 102:5811–5852. https:​/​/​doi​.org/​10​.3168/​jds​.2018​-15829. Hristov, A. N., M. Hanigan, A. Cole, R. Todd, T. A. McAllister, P. Ndegwa, and A. Rotz. 2011. Review: Ammonia emissions from dairy farms and beef feedlots. Can. J. Anim. Sci. 91:1–35. https:​/​ /​doi​.org/​10​.4141/​CJAS10034. Huhtanen, P., and A. N. Hristov. 2009. A meta-analysis of the effects of dietary protein concentration and degradability on milk protein yield and milk N efficiency in dairy cows. J. Dairy Sci. 92:3222– 3232. https:​/​/​doi​.org/​10​.3168/​jds​.2008​-1352. Huhtanen, P., J. I. Nousiainen, M. Rinne, K. Kytölä, and H. Khalili. 2008. Utilization and partitioning of dietary nitrogen in dairy cows fed grass silage-based diets. J. Dairy Sci. 91:3589–3599. https:​/​/​doi​ .org/​10​.3168/​jds​.2008​-1181. IPCC. 2006. IPCC Guidelines for National Greenhouse Gas Invento- ries. IGES: Intergovernmental Panel on Climate Change. IPCC. 2007. Climate Change 2007: Synthesis Report. Contribution of Working Groups I, II and III to the Fourth Assessment Report of the Intergovernmental Panel on Climate Change. IPCC. IPCC. 2019. Refinement to the 2006 IPCC Guidelines for National Greenhouse Gas Inventories. Accessed Jun. 15, 2021. https:​/​/​ www​.ipcc​-nggip​.iges​.or​.jp/​public/​2019rf/​pdf/​4​_Volume4/​19R​_V4​ _Ch11​_Soils​_N2O​_CO2​.pdf. James, G., D. Witten, T. Hastie, and R. Tibshirani. 2014. An Intro- duction to Statistical Learning: With Applications in R. Springer. Jonker, J. S., R. A. Kohn, and R. A. Erdman. 1998. Using milk urea nitrogen to predict nitrogen excretion and utilization efficiency in lactating dairy cows. J. Dairy Sci. 81:2681–2692. https:​/​/​doi​.org/​ 10​.3168/​jds​.S0022​-0302(98)75825​-4. Kauffman, A. J., and N. R. St-Pierre. 2001. The relationship of milk urea nitrogen to urine nitrogen excretion in Holstein and Jersey cows. J. Dairy Sci. 84:2284–2294. https:​/​/​doi​.org/​10​.3168/​jds​ .S0022​-0302(01)74675​-9. Kebreab, E. 2022. JDS 22 Bougouin. Mendeley Data, V1. https:​/​/​doi​ .org/​10​.17632/​2d3ff88t5g​.1. Kebreab, E., J. France, D. E. Beever, and A. R. Castillo. 2001. Ni- trogen pollution by dairy cows and its mitigation by dietary ma- nipulation. Nutr. Cycl. Agroecosyst. 60:275–285. https:​/​/​doi​.org/​ 10​.1023/​A:​1012668109662. Kebreab, E., A. B. Strathe, J. Dijkstra, J. A. N. Mills, C. K. Reynolds, L. A. Crompton, T. Yan, and J. France. 2010. Energy and protein interactions and their effect on nitrogen excretion in dairy cows. Pages 417–425 in 3rd EAAP International Symposium on Energy and Protein Metabolism and Nutrition. G. M. Crovetto, ed. Wa- geningen Academic Publishers. Kröber, T. F., D. R. Külling, H. Menzi, F. Sutter, and M. Kreuzer. 2000. Quantitative effects of feed protein reduction and methio- nine on nitrogen use by cows and nitrogen emissions from slurry. J. Dairy Sci. 83:2941–2951. https:​/​/​doi​.org/​10​.3168/​jds​.S0022​ -0302(00)75194​-0. Külling, D. R., H. Menzi, T. F. Kröber, T. F. Neftel, F. Sutter, P. Lischer, and M. Kreuzer. 2001. Emissions of ammonia, nitrous ox- ide and methane from different types of dairy manure during stor- age as affected by dietary protein content. J. Agric. Sci. 137:235– 250. https:​/​/​doi​.org/​10​.1017/​S0021859601001186. Lahart, B., S. McParland, E. Kennedy, T. Boland, T. Condon, M. Wil- liams, N. Galvin, B. McCarthy, and F. Buckley. 2019. Predicting the dry matter intake of grazing dairy cows using infrared reflec- tance spectroscopy analysis. J. Dairy Sci. 102:8907–8918. https:​/​/​ doi​.org/​10​.3168/​jds​.2019​-16363. Lin, L. I. K. 1989. A concordance correlation coefficient to evaluate reproducibility. Biometrics 45:255–268. https:​/​/​doi​.org/​10​.2307/​ 2532051. Marini, J. C., J. D. Klein, J. M. Sands, and M. E. Van Amburgh. 2004. Effect of nitrogen intake on nitrogen recycling and urea trans- porter abundance in lambs. J. Anim. Sci. 82:1157–1164. https:​/​/​ doi​.org/​10​.2527/​2004​.8241157x. Moore, K. D., E. Young, M. D. Wojcik, R. S. Martin, C. Gurell, G. E. Bingham, R. L. Pfeiffer, J. H. Prueger, and J. L. Hatfield. 2014. Ammonia measurements and emissions from a California dairy us- ing point and remote sensors. Trans. ASABE 57:181–198. https:​/​/​ doi​.org/​10​.13031/​trans​.57​.10079. Moriasi, D. N., J. G. Arnol