|
|
Bannink, A., van Lingen, H. J., Ellis, J. L., France, J., & Dijkstra, J. (2016). The contribution of mathematical modeling to understanding dynamic aspects of rumen metabolism. Frontiers in Microbiology, 7, 1820.
Abstract: All mechanistic rumen models cover the main drivers of variation in rumen function, which are feed intake, the differences between feedstuffs and feeds in their intrinsic rumen degradation characteristics, and fractional outflow rate of fluid and particulate matter. Dynamic modeling approaches are best suited to the prediction of more nuanced responses in rumen metabolism, and represent the dynamics of the interactions between substrates and micro-organisms and inter-microbial interactions. The concepts of dynamics are discussed for the case of rumen starch digestion as influenced by starch intake rate and frequency of feed intake, and for the case of fermentation of fiber in the large intestine. Adding representations of new functional classes of micro-organisms (i.e., with new characteristics from the perspective of whole rumen function) in rumen models only delivers new insights if complemented by the dynamics of their interactions with other functional classes. Rumen fermentation conditions have to be represented due to their profound impact on the dynamics of substrate degradation and microbial metabolism. Although the importance of rumen pH is generally acknowledged, more emphasis is needed on predicting its variation as well as variation in the processes that underlie rumen fluid dynamics. The rumen wall has an important role in adapting to rapid changes in the rumen environment, clearing of volatile fatty acids (VFA), and maintaining rumen pH within limits. Dynamics of rumen wall epithelia and their role in VFA absorption needs to be better represented in models that aim to predict rumen responses across nutritional or physiological states. For a detailed prediction of rumen N balance there is merit in a dynamic modeling approach compared to the static approaches adopted in current protein evaluation systems. Improvement is needed on previous attempts to predict rumen VFA profiles, and this should be pursued by introducing factors that relate more to microbial metabolism. For rumen model construction, data on rumen microbiomes are preferably coupled with knowledge consolidated in rumen models instead of relying on correlations with rather general aspects of treatment or animal. This helps to prevent the disregard of basic principles and underlying mechanisms of whole rumen function.
|
|
|
|
Hutchings, N. J., Özkan Gülzari, Ş., de Haan, M., & Sandars, D. (2018). How do farm models compare when estimating greenhouse gas emissions from dairy cattle production. Animal, 12(10), 2171–2180.
Abstract: The European Union Effort Sharing Regulation (ESR) will require a 30% reduction in greenhouse gas (GHG) emissions by 2030 compared with 2005 from the sectors not included in the European Emissions Trading Scheme, including agriculture. This will require the estimation of current and future emissions from agriculture, including dairy cattle production systems. Using a farm-scale model as part of a Tier 3 method for farm to national scales provides a more holistic and informative approach than IPCC (2006) Tier 2 but requires independent quality control. Comparing the results of using models to simulate a range of scenarios that explore an appropriate range of biophysical and management situations can support this process by providing a framework for placing model results in context. To assess the variation between models and the process of understanding differences, estimates of GHG emissions from four farm-scale models (DailyWise, FarmAC, HolosNor and SFARMMOD) were calculated for eight dairy farming scenarios within a factorial design consisting of two climates (cool/dry and warm/wet) x two soil types (sandy and clayey) x two feeding systems (grass only and grass/maize). The milk yield per cow, follower cow ratio, manure management system, nitrogen (N) fertilisation and land area were standardised for all scenarios in order to associate the differences in the results with the model structure and function. Potential yield and application of available N in fertiliser and manure were specified separately for grass and maize. Significant differences between models were found in GHG emissions at the farm-scale and for most contributory sources, although there was no difference in the ranking of source magnitudes. The farm-scale GHG emissions, averaged over the four models, was 10.6 t carbon dioxide equivalents (CO(2)e)/ha per year, with a range of 1.9 t CO(2)e/ha per year. Even though key production characteristics were specified in the scenarios, there were still significant differences between models in the annual milk production per ha and the amounts of N fertiliser and concentrate feed imported. This was because the models differed in their description of biophysical responses and feedback mechanisms, and in the extent to which management functions were internalised. We conclude that comparing the results of different farm-scale models when applied to a range of scenarios would build confidence in their use in achieving ESR targets, justifying further investment in the development of a wider range of scenarios and software tools.
|
|
|
|
Vitali, A., Felici, A., Esposito, S., Bernabucci, U., Bertocchi, L., Maresca, C., et al. (2015). The effect of heat waves on dairy cow mortality. J. Dairy Sci., 98(7), 4572–4579.
Abstract: This study investigated the mortality of dairy cows during heat waves. Mortality data (46,610 cases) referred to dairy cows older than 24 mo that died on a farm from all causes from May 1 to September 30 during a 6-yr period (2002-2007). Weather data were obtained from 12 weather stations located in different areas of Italy. Heat waves were defined for each weather station as a period of at least 3 consecutive days, from May 1 to September 30 (2002-2007), when the daily maximum temperature exceeded the 90th percentile of the reference distribution (1971-2000). Summer days were classified as days in heat wave (HW) or not in heat wave (nHW). Days in HW were numbered to evaluate the relationship between mortality and length of the wave. Finally, the first 3 nHW days after the end of a heat wave were also considered to account for potential prolonged effects. The mortality risk was evaluated using a case-crossover design. A conditional logistic regression model was used to calculate odds ratio and 95% confidence interval for mortality recorded in HW compared with that recorded in nHW days pooled and stratified by duration of exposure, age of cows, and month of occurrence. Dairy cows mortality was greater during HW compared with nHW days. Furthermore, compared with nHW days, the risk of mortality continued to be higher during the 3 d after the end of HW. Mortality increased with the length of the HW. Considering deaths stratified by age, cows up to 28 mo were not affected by HW, whereas all the other age categories of older cows (29-60, 61-96, and >96 mo) showed a greater mortality when exposed to HW. The risk of death during HW was higher in early summer months. In particular, the highest risk of mortality was observed during June HW. Present results strongly support the implementation of adaptation strategies which may limit heat stress-related impairment of animal welfare and economic losses in dairy cow farm during HW.
|
|
|
|
Carabano, M. J., Logar, B., Bormann, J., Minet, J., Vanrobays, M. L., Diaz, C., et al. (2016). Modeling heat stress under different environmental conditions. J. Dairy Sci., 99(5), 3798–3814.
Abstract: Renewed interest in heat stress effects on livestock productivity derives from climate change, which is expected to increase temperatures and the frequency of extreme weather events. This study aimed at evaluating the effect of temperature and humidity on milk production in highly selected dairy cattle populations across 3 European regions differing in climate and production systems to detect differences and similarities that can be used to optimize heat stress (HS) effect modeling. Milk, fat, and protein test day data from official milk recording for 1999 to 2010 in 4 Holstein populations located in the Walloon Region of Belgium (BEL), Luxembourg (LUX), Slovenia (SLO), and southern Spain (SPA) were merged with temperature and humidity data provided by the state meteorological agencies. After merging, the number of test day records/cows per trait ranged from 686,726/49,655 in SLO to 1,982,047/136,746 in BEL. Values for the daily average and maximum temperature-humidity index (THIavg and THImax) ranges for THIavg/THImax were largest in SLO (22-74/28-84) and shortest in SPA (39-76/46-83). Change point techniques were used to determine comfort thresholds, which differed across traits and climatic regions. Milk yield showed an inverted U-shaped pattern of response across the THI scale with a HS threshold around 73 THImax units. For fat and protein, thresholds were lower than for milk yield and were shifted around 6 THI units toward larger values in SPA compared with the other countries. Fat showed lower HS thresholds than protein traits in all countries. The traditional broken line model was compared with quadratic and cubic fits of the pattern of response in production to increasing heat loads. A cubic polynomial model allowing for individual variation in patterns of response and THIavg as heat load measure showed the best statistical features. Higher/lower producing animals showed less/more persistent production (quantity and quality) across the THI scale. The estimated correlations between comfort and THIavg values of 70 (which represents the upper end of the THIavg scale in BEL-LUX) were lower for BEL-LUX (0.70-0.80) than for SPA (0.83-0.85). Overall, animals producing in the more temperate climates and semi-extensive grazing systems of BEL and LUX showed HS at lower heat loads and more re-ranking across the THI scale than animals producing in the warmer climate and intensive indoor system of SPA.
|
|
|
|
Bernabucci, U., Biffani, S., Buggiotti, L., Vitali, A., Lacetera, N., & Nardone, A. (2014). The effects of heat stress in Italian Holstein dairy cattle. J. Dairy Sci., 97(1), 471–486.
Abstract: The data set for this study comprised 1,488,474 test-day records for milk, fat, and protein yields and fat and protein percentages from 191,012 first-, second-, and third-parity Holstein cows from 484 farms. Data were collected from 2001 through 2007 and merged with meteorological data from 35 weather stations. A linear model (M1) was used to estimate the effects of the temperature-humidity index (THI) on production traits. Least squares means from M1 were used to detect the THI thresholds for milk production in all parities by using a 2-phase linear regression procedure (M2). A multiple-trait repeatability test-model (M3) was used to estimate variance components for all traits and a dummy regression variable (t) was defined to estimate the production decline caused by heat stress. Additionally, the estimated variance components and M3 were used to estimate traditional and heat-tolerance breeding values (estimated breeding values, EBV) for milk yield and protein percentages at parity 1. An analysis of data (M2) indicated that the daily THI at which milk production started to decline for the 3 parities and traits ranged from 65 to 76. These THI values can be achieved with different temperature/humidity combinations with a range of temperatures from 21 to 36°C and relative humidity values from 5 to 95%. The highest negative effect of THI was observed 4 d before test day over the 3 parities for all traits. The negative effect of THI on production traits indicates that first-parity cows are less sensitive to heat stress than multiparous cows. Over the parities, the general additive genetic variance decreased for protein content and increased for milk yield and fat and protein yield. Additive genetic variance for heat tolerance showed an increase from the first to third parity for milk, protein, and fat yield, and for protein percentage. Genetic correlations between general and heat stress effects were all unfavorable (from -0.24 to -0.56). Three EBV per trait were calculated for each cow and bull (traditional EBV, traditional EBV estimated with the inclusion of THI covariate effect, and heat tolerance EBV) and the rankings of EBV for 283 bulls born after 1985 with at least 50 daughters were compared. When THI was included in the model, the ranking for 17 and 32 bulls changed for milk yield and protein percentage, respectively. The heat tolerance genetic component is not negligible, suggesting that heat tolerance selection should be included in the selection objectives.
|
|
