Assessment of Test-Day Milk Yield and Milk Composition in Tharparkar Cattle
Abstract
This study evaluated test-day milk yield and selected milk-quality traits in Tharparkar cattle using an observational dataset covering 91 animals and 237 animal-date recording observations collected between 07 March 2025 and 21 July 2026. The dataset contained four lactation classes (1–4) and records with and without calf suckling. Overall mean test-day milk yield was 5.48 ± 1.63 kg/day (range 2.06–11.61 kg/day). Mean milk fat, protein, SNF and lactose were 4.10%, 3.27%, 8.80% and 4.88%, respectively. One-way ANOVA did not detect a significant difference in mean test-day milk yield among lactation numbers (F=1.12, P=0.343). In contrast, records with calf suckling had higher observed milk yield than records without suckling (5.92 vs 5.19 kg/day); the unadjusted Welch test was significant (P=0.0006). A mixed-effects model accounting for repeated observations within animal also retained a significant positive association between calf-suckling status and recorded milk yield (β=0.90 kg/day, P<0.001). Milk fat showed a weak negative correlation with milk yield (r=-0.174, P=0.0074), whereas protein, SNF and lactose were not significantly correlated with yield. These findings provide a useful baseline for Tharparkar performance monitoring; however, the observational nature, unequal repeated records, and absence of detailed information on feed, days in milk, parity management, and suckling intensity limit causal interpretation.
Full Text
Tables and figures are in the full-text PDF.
Introduction
Tharparkar is an important indigenous dairy cattle breed adapted to the hot and relatively harsh environments of north-western India. Its value in dairy development is associated with adaptability, resilience and useful milk production. Published studies have reported substantial variation in lactation yield and production efficiency depending on herd, management, lactation stage and environmental conditions. For example, ICAR-NDRI reports a wet average of approximately 7.3 kg/day and a herd average of approximately 5.2 kg/day for Tharparkar cattle, with average fat and SNF values of 5.3% and 9.1%, respectively, under its farm conditions. Milk recording at the animal level is important for identifying superior animals, monitoring feeding and management responses, and supporting genetic improvement programmes. Milk composition is equally important because changes in fat, protein, lactose and SNF influence the nutritional and economic value of milk. The present dataset provides repeated field/farm observations on Tharparkar cattle and therefore offers an opportunity to quantify test-day performance and explore associations with lactation number and calf-suckling status. The specific objectives were: (i) to describe testday milk yield in Tharparkar cattle; (ii) to quantify major milk components; (iii) to compare milk yield across lactation numbers; (iv) to examine the association between calf-suckling status and recorded milk yield; and (v) to assess simple relationships between milk yield and milk composition.
Materials and Methods
Data source and study population
The analysis used the supplied Excel dataset containing 1383 component-level rows. All animals were identified as Tharparkar cattle.
Variables
The primary outcome was total milk yield (kg/day), calculated in the source data from morning, afternoon and evening entries. Milk-quality variables included fat (%), protein (%), solids-not-fat (SNF, %), lactose (%), and added water (%). Lactation number was classified as 1st, 2nd, 3rd or 4th lactation. Calf-suckling status was recorded as Yes or No.
Statistical analysis
Descriptive statistics included number of observations, number of animals, mean, standard deviation and range.
One-way analysis of variance (ANOVA) was used to compare test-day milk yield among lactation numbers. Welch's independent-samples t-test was used for the unadjusted comparison of records with and without calf suckling. Because multiple observations were available for some animals, a linear mixed-effects model was additionally fitted with AnimalID as a random intercept and lactation number and calf-suckling status as fixed effects. Pearson correlation coefficients were calculated between total milk yield and individual milk components. Statistical significance was assessed at P<0.05. The analysis is observational. Therefore, statistically significant associations should not be interpreted as evidence that calf suckling caused an increase in milk yield.
Important potential confounders such as days in milk, feed intake, season, parity-specific management, body condition and suckling frequency were not available in the supplied dataset.
Results
Table Mean test-day milk yield ranged from 5.23 kg/day in the 2nd/rd lactation group to 5.70 kg/day in the 4th lactation group. However, the overall lactation-number effect was not statistically significant (F=1.12, P=0.343). The absence of a significant overall lactation-number effect on test-day yield should not be interpreted as evidence that lactation number has no biological importance. The dataset includes only one to seven observations per animal and does not provide a complete lactation curve. A larger longitudinal dataset with days-in-milk information would be better suited to estimating lactation-specific trajectories and matureequivalent production.
The positive association between calf-suckling status and recorded milk yield differs from several controlled cow-calf studies in which suckling systems reduced saleable or machine-recorded milk and could alter milk fat. Such differences are plausible because the present field variable is simply a Yes/No status and does not describe suckling frequency, duration, timing relative to milking, or whether the recorded milk represents saleable milk, machine/milking yield, or a combined management measure. Consequently, the current finding should be treated as an empirical association that warrants targeted follow-up rather than a causal effect. The weak negative correlation between milk yield and fat is biologically plausible because dilution effects can occur as milk volume increases. Nevertheless, the correlation is small, indicating that most variation in milk yield was not explained by fat percentage alone. Protein, SNF and lactose were comparatively stable with respect to yield in the available observations.
Conclusion
The analysis of 237 animal-date observations from 91 Tharparkar cattle demonstrates an overall test-day milk yield of 5.48 ± 1.63 kg/day, with mean milk fat, protein, SNF and lactose of 4.10%, 3.27%, 8.80% and 4.88%, respectively. Lactation number was not significantly associated with test-day milk yield in the available records. Calf-suckling status was positively associated with recorded milk yield even after accounting for repeated observations within animal, but this result requires cautious interpretation because the dataset lacks information on suckling intensity and other important management and biological covariates. For future research, the recording system should capture animal-wise days in milk, parity, calving date, body condition score, feed and fodder intake, season, milking method, suckling frequency/duration, and health events. A larger longitudinal dataset analysed using mixed-effects or random-regression models would allow more robust estimation of lactation curves and management effects.
References
- Barth K. Effects of suckling on milk yield and milk composition of dairy cows in cow-calf contact systems. J Dairy Res. 2020;87(S1):133-137. doi:10.1017/S0022029920000515
- Fröberg S, Aspegren-Güldorff A, Olsson I, Marin B, Berg C, Hernández CS, et al. Effect of restricted suckling on milk yield, milk composition and udder health in cows and behaviour and weight gain in calves, in dual-purpose cattle in the tropics. Trop Anim Health Prod. 2008;[volume(issue)]:[pages not provided].
- Patel AK, Kachhawah SC, Patil NV, Chopra A. Effects of non-genetic factors on production and reproduction traits in Tharparkar cattle under arid region of Rajasthan. Indian J Anim Res. 2026;60(2):360-364. doi:10.18805/IJAR.B-4906
How this paper connects to the literature. Drag to explore, click any node to open that paper.
