|
|
||||||||



* Puyallup Research and Extension Center, Washington State University, Puyallup 98371
Department of Animal Sciences, The Ohio State University, Columbus 43210
Department of Animal Sciences, Washington State University, Pullman 99163
Department of Dairy Science, University of Wisconsin, Madison 53706
|| Church & Dwight Co. Inc., Princeton, NJ 08543
2 Corresponding author: jhharrison{at}wsu.edu
| ABSTRACT |
|---|
|
|
|---|
Key Words: urine urinary nitrogen sodium potassium
| INTRODUCTION |
|---|
|
|
|---|
Urinary excretion is the main method of homeostatic regulation for Na and K, thus intake of these minerals will directly affect Na and K excretion (Miller, 1975). In 2001, the Natural Resources Conservation Service standard for nutrient management planning (NRCS, 2001) stated that nutrient budgets should be established for N, P, and K. Although K has not been directly linked to environmental problems, excess field application of K results in a build up of soil K and must be considered by livestock producers during development of nutrient management plans. Greater soil concentrations of K can lead to higher concentrations of K in forages because grasses are luxury consumers of K (Harrison et al., 2003). High concentrations of K in forages will lead to greater intakes of K, which is a concern for periparturient dairy cows (NRC, 2001).
Prediction equations provide a means of estimating values that are difficult to measure on dairy operations and provide information to assist with the management of manure and nutrients. Evaluation of prediction equations is necessary to determine the accuracy of equations for animals with various housing and dietary situations. Determining the usefulness of previously published equations and development of equations that more accurately estimate urinary excretion will provide opportunities to continue to refine management strategies on dairy operations.
Nitrogen excretion is an environmental concern because of its relationship to ammonia volatilization and the potential movement of N to surface or groundwater. Equations that demonstrate relationships between the amount of N excreted in urine and factors affecting excretion are beneficial for selecting management strategies to reduce urinary N excretion. Previously published prediction equations have directly linked excretion of urinary N with MUN concentrations in dairy cattle (Jonker et al., 1998; Kauffman and St-Pierre, 2001; Kohn et al., 2002). Additional equations can further refine these relationships and provide additional options for estimating urinary excretion.
Innovations in management techniques are continually being explored to decrease risks of nutrients leaving livestock operations through air or water. An example of a management strategy designed to decrease ammonia volatilization is the separation of urine and feces in animal housing facilities. Development of new management strategies that focus specifically on urine management requires information to better estimate urine excretion on dairy operations and to gain a better understanding of factors affecting urine and urinary nutrient excretion. The main objectives of this study were: 1) to develop new equations to predict excretion of urine, urinary N, urinary Na, and urinary K; 2) to evaluate newly developed urine and urinary N excretion equations; and 3) to evaluate previously published equations developed to predict excretion of urine, urinary N, and urinary minerals. During the development of prediction equations, an effort was made to include equations for use on commercial dairy operations and to gain a better understanding of dietary and metabolic factors affecting excretion.
| MATERIALS AND METHODS |
|---|
|
|
|---|
|
The independent variables included in the LACT data set were BW, DIM, DMI, milk yield (MILK), MUN, milk fat (MF), and milk protein (Table 2
). The units associated with the independent variables were kg/animal for BW, kg/d for DMI and MILK, mg/dL for MUN (which was determined by colorimetry), and MF and milk protein percentages (Table 2
). Milk CP values in the data set were converted to percent milk true protein (MTP) values using a conversion factor of 0.9345 (Mackle et al., 1999). Other dietary characteristics and chemical compositions evaluated as prediction factors for equations included N intake (NI) in g/d, and dietary concentrations of CP and NDF. The supply of RDP (g/d) was determined for each animal by entering data from each individual cow-period into the 2001 Dairy NRC model. Dependent variables included in the LACT data set included urine excretion, kg/d, and urinary N excretion, g/d (UN). The LACT data set was used to develop equations that did not require mineral data to increase the number of cows used in the development of the regression equations. The MINERAL subset was only used when mineral intake or excretion data were included in the equation. The MINERAL subset included the same variables as the LACT data set and included the following independent variables: Na intake (NaI, g/d), K intake (KI, g/d), urinary Na concentration (%), and urinary K concentration (%). The DCAD was determined using the following equation: (Na + K) (S + Cl). The MINERAL subset also included the dependent variables of urinary Na excretion (UNa, g/d), and urinary K excretion (UK, g/d).
|
Equations given for each parameter included residual standard error (SE) and interstudy SE. Equations with lower SE are expected to provide a more precise estimation of excretion and should be used when values for the input variables are available.
Evaluations of previously published prediction equations for estimating urine excretion and UN were done using both the LACT data set and the VALIDATE data set. The VALIDATE data set included 48 observations of multiparous lactating Holstein cows from a study conducted by Wattiaux and Karg (2004). Animal characteristics for the VALIDATE data set are given in Table 2
. The VALIDATE data set was also used to evaluate new equations developed using the LACT data set when the necessary variables were available in the VALIDATE data set.
Evaluation of previously published prediction equations that included mineral intakes or dietary mineral concentrations as either dependent or independent variables were evaluated with the MINERAL subset. These equations could not be evaluated with either the LACT or VALIDATE data sets because mineral data were not available for cows in these data sets. The unavailability of mineral data also precluded the evaluation of new equations developed using the MINERAL data set.
Sample Collection and Analysis
Total collection metabolism studies (16 studies) conducted at Washington State University included multiparous lactating cows (372 cow-periods). Mineral data were available for a group of these cows (115 cow-periods) included in the MINERAL subset. Feeding, sample collection, and sample analyses were conducted by methods outlined by Johnson et al. (1998) and Timmermans et al. (2000). Studies from the University of Wisconsin that were included in the VALIDATE data set included 48 observations from multiparous lactating cows. Descriptions of feeding, sampling, and collections were summarized by Wattiaux and Karg (2004).
Mineral analyses were conducted on TMR, orts, feces, and urine samples for cows in the MINERAL subset. Minerals in feeds, orts, and feces were analyzed by the University of Nebraska Soil and Plant Analytical Laboratory (Lincoln, NE) via x-ray fluorescence spectrophotometry, and urine minerals were analyzed by Dairyland Laboratories (Arcadia, WI) using an inductively coupled plasma emission spectrophotometer. Samples were analyzed for various minerals including Ca, Cl, K, Mg, Mn, Na, P, and S. Mineral analyses were not conducted on milk samples. Therefore, milk mineral contents were assumed equivalent to values outlined in the 2001 Dairy NRC. Milk Na and K were estimated at 0.63 and 1.5 g/kg of milk, respectively (NRC, 2001).
Statistical Analyses
Regression analyses were performed using PROC MIXED of SAS (SAS Institute, 1999) with the discrete effect of study included as a random variable (St-Pierre, 2001). Equations were developed by running multiple iterations in MIXED and removing the least significant effect at each iteration. For the LACT data set, variables were kept if P < 0.10 because the number of observations were greater than 200. Variables were kept if P < 0.25 for the MINERAL subset due to the fewer number of observations (less than 200) available in the subset. Significant variables were compared by determining the errors for future predictions and selecting the independent variable that resulted in the least prediction error. Adjusted observations were calculated for graphing purposes by adding the residual from each individual observation to the predicted value of the study regression (St-Pierre, 2001).
Equations were evaluated by regressing residuals (predicted values subtracted from observed values) on the predicted values (St-Pierre, 2003). Predicted values were centered by subtracting the mean of all predicted values from each prediction. This makes the slope and intercept estimates in the regression orthogonal and, thus, independent. Mean biases were assessed using the intercepts of the regression equations and the slopes of the regression equations were used to determine the presence of linear biases.
| RESULTS AND DISCUSSION |
|---|
|
|
|---|
|
Animal and dietary factors were used as independent variables to determine significant factors for predicting urine excretion in the LACT data set. Factors related to urine excretion (P < 0.01) included DMI, NI, BW, MUN, and dietary NDF concentration. Days in milk, MILK, MTP, and MF were also evaluated in the regression equations, but were not related to (P > 0.14) urine excretion. In the LACT data set, MUN was the best individual factor for predicting urine excretion [Eq.1] based on the precision of future predictions. Similarly, MUN was determined to be the best factor for predicting urine excretion by Jonker et al. (1998) and Kohn et al. (2002).
![]() | [1] |
Body weight and DMI were added to the MUN prediction equation to see if the additional variables improved the prediction of urine excretion. When BW, DMI, and MUN were evaluated, there was a slight reduction in the residual SE, from 5.8 for [1] to 5.5, but there was not any improvement in prediction across studies. An equation for estimating urine excretion with MUN as the only variable provides a simple means by which producers can estimate total urine excretion on their operations.
The VALIDATE data set was used to evaluate [1]. Evaluation of [1] resulted in an average mean bias of 5.26, (2.10) kg/d, indicating that the use of [1] for cows in the VALIDATE data set underestimated the urine excretion of these cows by an average of 5.3 kg/d. However, the error between measurements was greater than the mean bias determined during the evaluation of [1]. The linear bias was not significant (P = 0.85), demonstrating that [1] predicted urine excretion evenly across lesser and greater predicted values in the VALIDATE data set.
Effect of Minerals on Urine Excretion.
The MINERAL subset was used to predict the effect of mineral intakes and urinary mineral concentrations on urine excretion because the LACT data set did not include mineral analyses. Intakes of Na and Cl in the MINERAL subset were closely related (P < 0.01) to urine excretion. The intake of Na and the N use efficiency, as measured by MUN, were the factors that accounted for the greatest variation within studies (interstudy SE) in the MINERAL subset [Eq. 2].
![]() | [2] |
When the MINERAL subset was used to develop a simple linear equation with MUN as the only independent variable, the error of future predictions was 8.38 compared with 7.79 for [2]. Therefore, the comparison of [2] with the equation containing only MUN indicated that the precision of estimating urine excretion improved when NaI and MUN were used as variables for predicting urine excretion. Although NaI is not commonly available on many dairy operations, future improvements in measurements of intake and dietary mineral composition could increase the possibility of including NaI as a variable for estimating urine excretion on individual dairy operations.
In the MINERAL subset, there was a relationship (P < 0.01) between the DCAD and urine excretion (Figure 1
). The DCAD was the best single predictor of urine excretion between studies in the MINERAL subset [Eq.3]. Prediction of urine excretion using DCAD was 14.4% more precise than when MUN was used as the only predictor and 8.0% more accurate than [2].
|
![]() | [3] |
The effect of DCAD on urine excretion is not surprising because urinary excretion is the main homeostatic regulation for Na, K, and Cl (Miller, 1975; Maltz and Silanikove, 1996). When NI, KI, and NaI were all used as predictors of urine excretion, only NaI was significant (P < 0.01), indicating the NaI had a greater effect on urine excretion in the MINERAL subset than KI or NI.
Evaluation of KI as the only variable included for prediction of urine excretion resulted in a positive relationship of KI to urine excretion (P < 0.12). Fisher et al. (1994) also reported that KI had a significant affect on urine excretion, and Bannink et al. (1999) found that KI was significant when included with NaI. Overall, [3] had lower residual and interstudy SE than the equation that only included KI, and would be expected to provide better (more precise) estimates of urine excretion.
Concentrations of minerals in urine and urinary mineral excretion were also evaluated as predictors of urine excretion. Urinary minerals were evaluated to determine if urine excretion was more closely related to urinary concentrations of minerals or to total mineral excretion. Urine excretion in the MINERAL subset was related to both concentrations of Na, K, and S (P < 0.01) in the urine and urinary excretions of N, K, and Na (P< 0.01).
An equation for predicting urine excretion, that included urinary excretion of N and Na and the concentration of K in urine as independent variables, [Eq.4] was developed using the MINERAL subset:
![]() | [4] |
Equation 4 indicated that the excretion of UN and UNa affected the amount of urine excreted by cows in the MINERAL subset more than the concentrations of N and Na in urine. Compared with N and Na, the concentration of K in urine had a greater effect on urine excretion than did total UK. The role of Na and K in urine excretion could affect the concentration of N in the urine. Maltz and Silanikove (1996) found an increase in urea concentrations in urine when there were reductions in Na, K, and Cl ions. An increase in the urinary N concentration may affect the rate of ammonia volatilization on a dairy operation.
Equations [2], [3], and [4] could not be evaluated with the VALIDATE data set due to the unavailability of mineral data on cows in the data set. However, the MINERAL data set was used to evaluate previously published equations.
Evaluation of Published Urine Excretion Equations.
Presently, few equations for predicting urine excretion for lactating dairy cows are available in the literature. A published equation [(PE) PE1] was proposed by Fox et al. (2004) to predict urine excretion (Table 4
). When PE1 was evaluated with the VALIDATE data set, the mean bias was 8.96 (±1.10) kg/d and there was not a significant linear bias (P = 0.40). Therefore, PE1 adequately described the lesser and greater predicted values in the VALIDATE data set, but underestimated the average urine excretion.
|
In 1999, Bannink et al. developed an equation (PE2) for predicting urine excretion using NaI and KI (Table 4
). When PE2 was evaluated using the MINERAL subset, there was a linear bias (0.69; P < 0.01) and mean bias (7.3 kg/d; P < 0.01). The standard deviation of the residuals from the evaluation of PE2 was 9.08; thus, the mean bias was less than the expected variation between studies, and the use of PE2 would provide an accurate average excretion value. The linear bias resulted in an overestimation of less than 23 kg/d at the maximum predicted urine excretion (54.6 kg/d), indicating that PE2 was not accurate across the range of data in the MINERAL subset.
Daily urinary excretions of N, Na, and K have been used to predict urine excretion (Bannink et al., 1999). In 1999, Bannink et al. proposed PE3 (Table 4
) to predict urine excretion. When PE3 was evaluated using the MINERAL subset, there were mean and linear biases (P < 0.05). When the mean bias was evaluated, the standard error of the residuals was greater than the mean bias. The mean bias, even though it was significant, resulted in an average estimate of urine excretion that was less than the variation expected between studies. The value of the linear bias for the minimum predicted urine excretion of 7.2 kg/d was 8.1 kg/d, slightly greater than the standard error of 6.5 kg/d. Although PE3 underestimated urinary excretion for cows with lower predicted urinary excretion and overestimated urine excretion for cows with the greatest predicted urinary excretion, it adequately estimated the average urinary excretion for cows in the MINERAL data set.
The final equation evaluated for prediction of urine excretion was developed by Bannink et al. (1999) and combined several factors including mineral intakes, MILK, and MTP (PE4; Table 4
). Evaluation of PE4 with the MINERAL subset resulted in mean (7.06) and linear (0.693) biases (P < 0.01). Similar to the previous equations evaluated, the mean bias, although significant, was less than the variation expected between studies. The linear bias resulted in an estimated urine excretion of less than 23 kg/d greater than was observed at the maximum predicted urine excretion (55.3 kg/d) and was outside the SE of 8.3.
Urinary Nitrogen Excretion
Urinary N excretion averaged 49% of total manure N excretion for dairy cows in the LACT data set. Excess dietary N is excreted via urine in the form of urea. Urea is a soluble compound that will diffuse to various body fluids, including blood and milk, in the lactating cow (Kauffman and St-Pierre, 2001; Nousiainen et al., 2004). In addition, the concentration of MUN is directly affected by dietary CP concentrations (Nousiainen et al., 2004). Because excess dietary N is excreted mostly in the form of urea, MUN would be expected to have a direct effect on UN. As expected, MUN was positively related to UN (P < 0.01) in the LACT data set.
Equations for estimating UN were developed using the LACT data set. Independent variables related (P < 0.10) to UN included BW, NI, DMI, DIM, and MUN. The best individual predictor of UN within studies in the LACT data set was NI. An increase in UN has been associated with greater N intakes in several studies (Sannes et al., 2002; Broderick, 2003; Wattiaux and Karg, 2004). In 1998, Jonker et al. proposed an equation for predicting UN that included NI, and Nousiainen et al. (2004) found that the dietary CP concentration was a better predictor or UN than was MUN.
The RDP supply was used to evaluate the relationships between protein type and UN. The best simple linear equation [Eq. 5] for predicting UN both within and between studies in the LACT data set included RDP supply as determined by the 2001 NRC model (Figure 2
).
|
![]() | [5] |
Multiple regression equations were also evaluated for the prediction of UN [Eq. 6].
![]() | [6] |
Although the multiple regression equation [Eq. 6] decreased the variation within studies, the multiple linear regression equation did not significantly improve the prediction of UN across studies. The errors associated with each of the independent variables in the multiple linear regression equation would also decrease the preciseness of the prediction of urinary N excretion.
Equation [5], the best equation for the prediction of UN in the LACT data set, was evaluated using the VALI-DATE data set. Evaluation of [5] with the VALIDATE data set did not result in a linear bias (P = 0.90), and the mean bias was only 9.6 (±4.5) g/d and was less than the standard deviation of the residuals. The lack of biases during the evaluation of [5] indicated that this equation was an accurate estimator of urinary N excretion for cows in the VALIDATE data set.
When equation 6 was evaluated using the VALIDATE data set, there was no linear bias (P = 0.22) and the average mean bias was 43.86 (±4.20). Although multiple regression equations such as [6] may provide insights into animal metabolism, these equations are often not practical for application on dairy operations.
Evaluation of Published Urinary Nitrogen Excretion Equations.
Numerous equations have been proposed to predict UN from lactating dairy cows (Table 4
). Many of these previously published prediction equations have included MUN in the equation because equilibration of urea in the body results in a direct relationship to UN (Kauffman and St-Pierre, 2001; Nousiainen et al., 2004). An additional advantage of the use of MUN as an independent variable for prediction of UN is the ease and low cost of determining MUN. The challenge of using MUN as an independent variable is a result of the differences in the calibration and accuracy of the laboratory equipment used to determine the MUN values for many commercial dairy operations (Kohn et al., 2004).
Previously published equations for predicting UN have been evaluated by several researchers (Kauffman and St-Pierre, 2001; Kohn et al., 2002; Nousiainen et al., 2004). Besides evaluation of equations, several of these researchers proposed additional equations for prediction of UN. Evaluation of previously published equations provides information as to the usefulness of the prediction equations in various dietary and management situations.
Jonker et al. (1998) proposed the use of PE5 to predict UN. In 2001, Kauffman and St-Pierre reported that the coefficient (12.54) from PE5 should be increased to 17.6 to account for an inaccuracy in the determination of the equipment used for calculation of MUN values. The adjustment in MUN values was confirmed by Kohn et al. (2002). When PE5 was evaluated using both the LACT and VALIDATE data sets, there were mean biases (P < 0.01) of 76.7 and 55.3 g/d, respectively, which further confirmed the need for an adjustment to the coefficient of the MUN term in PE5 (Jonker et al., 1998).
Another equation for predicting UN using adjusted MUN values was published by Kauffman and St-Pierre (2001). When PE6 was evaluated with the LACT data set, there was not a mean bias (P = 0.30), but there was a linear bias (0.672; P < 0.01). The linear bias resulted in a maximum bias of 175 g/d over the entire range of the predicted values.
In 2002, Kohn et al. evaluated equations published by Jonker et al. (1998) and Kauffman and St-Pierre (2001) and proposed an additional equation for prediction of UN (PE7). Evaluation of PE7 with the LACT data set resulted in a linear bias of 0.617 (P < 0.01). Equation PE7 was comparable to PE6 for prediction of UN in the LACT and VALIDATE data sets, with a maximum bias of 135 g/d over the full range of the data.
Prediction equations for UN were developed for both Holsteins and Jerseys by Kauffman and St-Pierre (2001). The development of prediction equations across breeds led to an adjustment for BW in the UN prediction equation. Evaluation of PE8 with the LACT data set did not result in a mean bias, but there was a linear bias (P < 0.01) of 0.529 (Figure 3
). The maximum bias for PE8 was less than 122 g/d over the range of the predicted values and resulted in the least bias of the equations evaluated for prediction of UN. When the VALIDATE data set was used to evaluate PE8, the mean bias was less than the standard deviation of the residuals, and the linear bias for the minimum predicted value (129.6 g/d) was less than 65 g/d. Equation PE8, by Kauffman and St-Pierre (2001), was designed to be applicable across breeds and provides a means of determining the effects of MUN values on UN.
|
Urinary Mineral Excretion
Sodium.
Average UNa in the MINERAL subset was 84.8 g/d (Table 3
), and accounted for 71.4% of total Na excretion. Silanikove et al. (1997) reported that 65% of total Na excretion was via urine and Bannink et al. (1999) found UNa accounted for 76% of total Na excretion. In the MINERAL subset, the dietary Na concentration was above animal requirements due to the addition of sodium bicarbonate to the diets, thus resulting in greater urinary sodium concentration in the MINERAL subset than reported by Bannink et al. (1999), Shalit et al. (1991), and Tucker et al. (1988).
Urinary Na excretion was directly related to NaI in the MINERAL subset as it was in the study reported by Bannink et al. (1999). Equation 7 was developed to estimate UNa based on NaI:
![]() | [7] |
Urinary Na excretion has also been estimated by the difference in apparently digested Na and Na secreted in milk (Bannink et al., 1999). However, the use of apparently digested Na to determine Na excretion may not always be accurate due to Na recycling that occurs for ruminant animals, especially during periods of low NaI (NRC, 2001).
Potassium.
Dietary concentrations of K in the MINERAL subset averaged 1.34% of DM, lower than the dietary K concentrations in many modern dairy diets due to high K levels in forages. The lower dietary K concentration in the MINERAL subset resulted in KI that were 26% less than reported by Bannink et al. (1999). Fisher et al. (1994) reported UK of 209 g/d for diets containing 1.6% K. The lesser KI in our study resulted in approximately 50% of the UK reported by Bannink et al. (1999), and in lower concentrations of K in urine than reported by Shalit et al. (1991). In contrast, Tucker et al. (1988) reported urinary K concentrations that were similar to concentrations in our data set when the dietary concentration of (Na + K) Cl equaled 10 mEq/100 g of diet DM in their study.
Increased KI resulted in greater UK in the MINERAL subset [Eq. 8].
![]() | [8] |
In the MINERAL subset, 75% of K excretion occurred via urinary routes. Conversely, Silanikove et al. (1997) reported that ~64% of K excreted in manure was through urinary excretion, and Bannink et al. (1999) found that urinary excretion accounted for over 87% of total K excretion.
Evaluation of Published Urinary Na and K Excretion Equations.
Equation PE10 was proposed in 1999 by Bannink et al. to estimate UNa (Table 4
). When PE10 was evaluated using the MINERAL subset, the linear bias (P < 0.01) was less than 66 g/d over the range of predicted values. The mean bias was not significant (P = 0.58) and estimated the average UNa of the cows in the MINERAL subset.
Bannink et al. (1999) proposed an equation to estimate UK (PE11; Table 4
) that required an estimation of apparently digested K. There was not a mean bias when PE11 was evaluated using the MINERAL subset, but there was a linear bias (P < 0.01). The linear bias was less than 87 g/d for the maximum predicted value of 339.4 g/d. However, the amount of apparently digested K is not an easily obtainable value in practical situations and PE11 would not be practical for use on an on-farm basis.
| CONCLUSIONS |
|---|
|
|
|---|
| FOOTNOTES |
|---|
Received for publication December 22, 2004. Accepted for publication August 26, 2005.
| REFERENCES |
|---|
|
|
|---|
This article has been cited by other articles:
![]() |
J. M. Powell, G. A. Broderick, and T. H. Misselbrook Seasonal Diet Affects Ammonia Emissions from Tie-Stall Dairy Barns J Dairy Sci, February 1, 2008; 91(2): 857 - 869. [Abstract] [Full Text] [PDF] |
||||
![]() |
J. M. Powell, T. H. Misselbrook, and M. D. Casler Season and Bedding Impacts on Ammonia Emissions from Tie-stall Dairy Barns J. Environ. Qual., January 4, 2008; 37(1): 7 - 15. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |