Acessibilidade / Reportar erro

Evaluation of the Factors Affecting the Mortality Rate in Poultry Transport via Panel Data Analysis

ABSTRACT

In this study, we evaluated some of the factors that affect mortality rate during transport from broiler poultry houses to slaughterhouses by conducting panel data analysis. We analyzed the data obtained from 26,599 broiler farms transported to the slaughterhouse from contracted broiler farms in 11 provinces in Turkey. Augmented Dickey-Fuller (ADF) and Phillips-Perron (PP) tests were performed to evaluate whether the series forming the dataset were stationary. To analyze individual effects, parameters were estimated using fixed and random effects models. To decide which of the two models was valid, the Hausman test and fixed effects panel data analysis were performed. The fixed effects model explained 90.93% of the changes in the mortality rate through the independent variables. In the non-periodical fixed-effect panel data analysis, the mortality rate shows a significant (p<0.01) effect between the number of animals alive, number of animals dead, mean and total live weight at poultry farm exit, mean live weight at slaughterhouse arrival, and different variables. Our results showed that the estimation equations developed using model parameters to determine the mortality rate during transportation from different provinces could contribute to effective production planning.

Keywords:
Poultry; transportation; panel data; Turkey; mortality rate

INTRODUCTION

Animals, including farm animals, are transported for various purposes, such as selling (in animal markets) and food processing (slaughterhouses). During transport, animals experience stress, which can differ between species and individuals within a species. Broiler chickens are highly sensitive to stress (Guerrero-Legarreta, 2010Guerrero-Legarreta I. Handbook of poultry science and technology primary processing. London: John Wiley & Sons; 2010.). As broiler chickens are very commonly transported worldwide, the economic losses resulting from the risk of transportation are important (EFSA, 2004).

The effect of human and environmental factors on broiler chicken transport from the poultry farm to the slaughterhouse can significantly decrease production due to stress. These losses can occur while catching individuals, loading them into vehicles, transporting them, and waiting at the slaughterhouse (Warriss et al., 1992Warriss PD, Bevis EA, Brown SN, et al. Longer journeys to processing plants are associated with higher mortality in broiler chickens. British Poultry Science 1992;33:201-6. https://doi.org/10.1080/00071669208417458
https://doi.org/10.1080/0007166920841745...
; Nijdam et al., 2004Nijdam E, Arens P, Lambooij E, et al. Factors influencing bruises and mortality of broilers during catching, transport and lairage. Poultry Science 2004;83:1610-5. https://doi.org/10.1093/ps/83.9.1610
https://doi.org/10.1093/ps/83.9.1610...
; Schwartzkopf-Genswein et al., 2012Schwartzkopf-Genswein KS, Faucitano L, Dadgar S, et al. Road transport of cattle, swine and poultry in North America and its impact on animal welfare, carcass and meat quality: a review. Meat Science 2012;92:227-43. https://doi.org/10.1016/j.meatsci.2012.04.010
https://doi.org/10.1016/j.meatsci.2012.0...
; Vieira et al., 2015Vieira FMC, Deniz M, da Silva, IJO, et al. Pre-slaughter losses of broilers: effect of time period of the day and lairage time in a subtropical climate. Semina: Ciências Agrárias 2015;36(6):3887-96. https://doi.org/10.5433/1679-0359.2015v36n6p3887
https://doi.org/10.5433/1679-0359.2015v3...
; Vecerek et al., 2016Vecerek V, Voslarova E, Conte F, et al. Negative trends in transport-related mortality rates in broiler chickens. Asian-Australasian Journal of Animal Sciences 2016;29(12):1796-804. https://doi.org/10.5713%2Fajas.15.0996
https://doi.org/10.5713%2Fajas.15.0996...
).

Broilers are rested before slaughtering to reduce stress levels caused by transport, but prolonged resting periods can increase stress due to hunger, thirst, and frequency of settlement, thus increasing the mortality rate (Nijdam et al., 2004Nijdam E, Arens P, Lambooij E, et al. Factors influencing bruises and mortality of broilers during catching, transport and lairage. Poultry Science 2004;83:1610-5. https://doi.org/10.1093/ps/83.9.1610
https://doi.org/10.1093/ps/83.9.1610...
; Warris, 2010).

Some studies have shown that transportation distance and the time taken to transport broilers to the slaughterhouse at slaughter age significantly affect live weight loss. Longer transport times and distances can increase live weight loss in broilers (Ondrašovičová et al., 2008; Oba et al., 2009Oba, A, Almeida M, Pinheiro JW, et al. The effect of management of transport and lairage conditions on broiler chicken breast meat quality and DOA (Death on Arrival). Brazilian Archives of Biology and Technology 2009;52:20511. https://doi.org/10.1590/S1516-89132009000700026
https://doi.org/10.1590/S1516-8913200900...
; Aral et al., 2014Aral Y, Arikan MS, Akin AC, et al. Economic losses due to live weight shrinkage and mortality during the broiler transport. Veterinary Journal of Ankara University 2014;61:205-10. https://doi.org/10.1501/Vetfak_0000002630
https://doi.org/10.1501/Vetfak_000000263...
).

The amount of waste and the losses incurred during transport from the poultry farm to the slaughterhouse is influenced by several factors, which include the mean body weight of the chicken (Nijdam et al., 2004Nijdam E, Arens P, Lambooij E, et al. Factors influencing bruises and mortality of broilers during catching, transport and lairage. Poultry Science 2004;83:1610-5. https://doi.org/10.1093/ps/83.9.1610
https://doi.org/10.1093/ps/83.9.1610...
), the number of animals being transported (Aral et al., 2014Aral Y, Arikan MS, Akin AC, et al. Economic losses due to live weight shrinkage and mortality during the broiler transport. Veterinary Journal of Ankara University 2014;61:205-10. https://doi.org/10.1501/Vetfak_0000002630
https://doi.org/10.1501/Vetfak_000000263...
; Arıkan, 2017), the mortality rate (Lupo et al., 2009Lupo C, Le Bouquin S, Balaine L, et al. Chauvin C. Feasibility of screening broiler chicken flocks for risk markers as an aid for meat inspection. Epidemiology & Infection 2009; 137:1086-98. https://doi.org/10.1017/S095026880900209X
https://doi.org/10.1017/S095026880900209...
), lairage time at the slaughterhouse (Vieira et al., 2010Vieira FMC, Silva IJOD, Barbosa Filho JAD, et al. Productive losses on broiler preslaughter operations: effects of the distance from farms to abattoirs and of lairage time in a climatized holding area. Revista Brasileira de Zootecnia 2010;39:2471-6. https://doi.org/10.1590/S1516-35982010001100021
https://doi.org/10.1590/S1516-3598201000...
; Lupo et al., 2009; Pirompud et al., 2022Pirompud P, Sivapirunthep P, Chaosap C. Factors influencing dead on arrival and condemnation of broilers during catching, transport, lairage and slaughter. International Journal of Agricultural Technology 2022;18(6):2535-44.), and microclimate conditions (Santos et al., 2020Santos VM dos, Dallago BS, Racanicci AM, et al. Effect of transportation distances, seasons and crate microclimate on broiler chicken production losses. PLoS One 2020;15(4);e0232004. https://doi.org/10.1371/journal.pone.0232004
https://doi.org/10.1371/journal.pone.023...
).

In this study, we conducted panel data analysis to assess the factors that increase mortality rate during the transport of broiler chickens from poultry farms to the slaughterhouse. These factors included the number of live animals being transported, the number of deaths during transport, the mean and total live weight at the exit from the poultry farm, the mean live weight upon arrival at the slaughterhouse, the difference in live weight between poultry farm exit and slaughterhouse arrival, and the lairage duration at the slaughterhouse.

MATERIALS AND METHODS

Dataset

We collected the data and transport records of 83,840,909 chickens that were transported to the slaughterhouse of a company. This involved 26,599 transports from contracted broiler farms in 11 provinces affiliated with an integrated poultry company operating in Turkey in 2022.

To examine the factors affecting mortality rate (MrRt), which was considered to be a dependent variable in this study, we analyzed several independent variables, including the number of animals transported from the poultry farm to the slaughterhouse (TNa), the number of live animals arriving at the slaughterhouse (LaNa), the number of dead animals in the slaughterhouse (ExNa), the mean live weight at poultry farm exit (MLWe), the total body weight at slaughterhouse arrival (TBW), the mean live weight of the hens upon arrival at the slaughterhouse (MLWa), the live weight difference (difference) between the poultry farm and slaughterhouse, and the lairage time in the slaughterhouse (LT). The data were analyzed using the EViews 8 Enterprise Edition software (EVIEWS, 2016).

Method

Panel data, which are defined as time series of cross-sections or cross-sectional data of the time series (Greene, 2003), can also be interpreted as the expression of cross-sectional observations of specific units, such as firms, countries, or households considering the dimension of time (Baltagi, 2001Baltagi, B. Econometric analysis of panel data. 2nd ed. London: John Wiley & Sons; 2001.).

Panel regression models, in which datasets containing cross-section and time series combinations are used, include many methods such as one-way and two-way fixed effects and random effects models, dynamic panel analysis, and generalized least squares (LCC). In this study, one-way fixed effects and random effects models were used.

Panel Unit Root Test

To determine the stationarity of the variables, a panel unit root test was conducted using the method proposed by Im et al. (2003Im KS, Pesaran MH, Shin Y. Testing for unit roots in heterogeneous panels. Journal of Econometrics 2003;115:53-74. https://doi.org/10.1016/S0304-4076(03)00092-7
https://doi.org/10.1016/S0304-4076(03)00...
).

This test can be performed to calculate the mean test statistic of augmented Dickey-Fuller (ADF) tests for each unit in the panel (Saraçoğlu & Doğan, 2005).

For the panel unit root test, the model used was Δyit = αi + βi yi,t-1 + eit,

Here, i = 1, ..., N and t = 1, ..., T,

were defined by Im et al. (2003Im KS, Pesaran MH, Shin Y. Testing for unit roots in heterogeneous panels. Journal of Econometrics 2003;115:53-74. https://doi.org/10.1016/S0304-4076(03)00092-7
https://doi.org/10.1016/S0304-4076(03)00...
). The null hypothesis (H0) was

βi = 0 for all i,

and the alternative hypothesis (H1) was βi < 0 for i = 1, 2, ..., N1 and βi = 0 for i = N1 + 1, N1 + 2, ..., N.

Accepting the null hypothesis indicates the presence of a panel unit root, while accepting the alternative hypothesis suggests the absence of a panel unit root. Im et al. (2003Im KS, Pesaran MH, Shin Y. Testing for unit roots in heterogeneous panels. Journal of Econometrics 2003;115:53-74. https://doi.org/10.1016/S0304-4076(03)00092-7
https://doi.org/10.1016/S0304-4076(03)00...
) tested the “no unit root” hypothesis using t-bar statistics.

One-Way Fixed Effects Model

To analyze the panel data, the variables were represented using two subscripts to indicate the time and cross-sectional dimensions. This differed in terms of the time series and cross-sectional data. The model Yit = αi + X’it β + eit is a fixed-effects model, where i is the cross-section index and t is the time index. The fixed-effects model can be estimated using the within-group estimator and the least squares with a dummy variable estimator (LSDV), which can satisfy the basic assumptions (Greene, 2003).

Y i t = a i + X i t b + e i t (1)

i = 1, …………N

t = 1,…………..T

In the model, Xit, represents the vector of explanatory variables, Yit represents the dependent variable, β represents the slope coefficients, eit denotes the error term, and αi represents the constant term unit effect. We analyzed the effect of time and units, assuming that the constant term did not change with time but could vary for each unit. The constant term did not change between units, but varied over time. To determine the coefficients of the model (Yit = αi + X’it β + eit) using the within-group estimator, the mean values of the observations of each individual were subtracted from the observations of that individual. Then, using this transformed data, the least squares (LCC) method was used for estimation (Kennedy, 2006Kennedy, P. Ekonometri kilavuzu, (çev. Muzaffer Sarimeseli ve Senay Açikgöz), 5. Baski, Ankara: Gazi Kitabevi; 2006.).

The model was analyzed based on two effects: the group effect, which assumes that the constant term does not change over time but can vary across units, and the time effect, which assumes that the constant term does not change across units, but can change over time.

An LCC estimator with a dummy variable for each unit can be used as an alternative method to account for the differences in constant terms. However, this method, also known as LSDV, might result in a reduction in the degree of freedom and the problem of multicollinearity, due to the use of many dummy variables (Kennedy, 2006Kennedy, P. Ekonometri kilavuzu, (çev. Muzaffer Sarimeseli ve Senay Açikgöz), 5. Baski, Ankara: Gazi Kitabevi; 2006.). When using a dummy variable for each unit, the fixed effects model in equation (1) can be expressed as follows (Pazarlıoğlu & Gürler, 2007):

Y i t = a 1 D 1 + . + a N D N + X i t b + e i t (2)

Both models assume that the differences between units or times are due to differences in constant terms (Greene, 2003). Therefore, it is assumed that variable coefficients do not vary between units or times. To investigate the group effect, we assumed that the constant term did not change over time, but could vary between units. To investigate the time effect, the constant term was assumed to not change across units, but to potentially vary over time.

To determine the difference between the units in the fixed effects model, a group significance test was performed. The following F statistic was evaluated under the null hypothesis that the constant term did not change between units (Greene, 2003):

F ( N 1, N T N K ) = ( R L S D V 2 R Pooled  2 ) / ( N 1 ) ( 1 R L S D V 2 ) / ( N T N K ) (3)

In the F statistic (3), represents the coefficient of determination of the LSDV model, and represents the coefficient of determination obtained from estimating the panel data with EKK. T represents the number of observations for each unit, N represents the number of units (groups), and K represents the number of explanatory variables. When the computed F statistic is greater than the table value, the null hypothesis is rejected, indicating the presence of a group effect or a difference between units.

The same test statistic can be used to determine the difference over time. In this case, the LSDV model with a time-varying constant term was used, and the null hypothesis stated that the constant term did not vary with time.

One-Way Random Effects Model

An alternative model to consider is the random effects model. When individual effects are not related to explanatory variables and the constant terms of the units are randomly distributed across units, the model should be structured accordingly (Greene, 2003).

In random effects models, the variations associated with cross-sections and/or time are included as a component of the error term in the model. The advantage of random effects models over fixed effects models is the elimination of the loss of degrees of freedom (Baltagi, 2001Baltagi, B. Econometric analysis of panel data. 2nd ed. London: John Wiley & Sons; 2001.).

In this study, we used a one-way random effects model, which assumed that the difference between the cross-sections was a component of the error terms in the model. The model was estimated with i representing the cross-sections and t representing the time, as shown in equation (4):

Y i t = α + X i t β + ( μ i + v i t ) i = 1, N t = 1, T E ( u i ) = ( v i t ) = 0, C o v ( u t , v j t ) = σ u , v , V a r ( u i ) = σ u 2 v e E ( X i t , u t ) = 0 (4)

In this model, Xit represents the vector of explanatory variables, Yit represents the dependent variable, β represents the vector of variable coefficients, and α represents the constant term. The error terms in the model are assumed to be independently and identically distributed with zero variance. The error term µi represents unobservable random differences between the units and vit represents the remaining errors. The individual error terms µi, which express the cross-sectional effect, are not related to each other, and vit is not related to the panel error term. By combining these two error terms, a model can be obtained, assuming that the data follows a normal distribution.

Y i t = a i + X i t b + e i t (5)

e i t = µ i + v i t , (6)

The error terms in this model consist of two components, and the error term variance (6) does not exhibit constant variance and zero covariance properties. Thus, the LCC estimator cannot be applied to this model, since the error terms lack the desired properties. Instead, methods such as appropriate versions of the Generalized Least Squares (GLS) method can be used. To apply the GLS method to these models, the error term and variance components must be known. In this study, the variance components were determined using the methods described by Swamy & Arora (1972Swamy PAVB, Arora SS. The Exact Finite Sample Properties of the Esitmator of Coefficients in the Error Components Regression Model. Econometrica 1972;40(2):261-75. https://doi.org/10.2307/1909405
https://doi.org/10.2307/1909405...
) and Wallace & Hussain (1969Wallace TD, Hussain A. The use of error components models in combining cross section with time series data. Econometrica: Journal of the Econometric Society 1969;37(1);55-72. https://doi.org/10.2307/1909205
https://doi.org/10.2307/1909205...
). Swamy & Arora (1972) recommended obtaining variance components from within-group and between-group regression models (Baltagi, 2001Baltagi, B. Econometric analysis of panel data. 2nd ed. London: John Wiley & Sons; 2001.). In this study, we estimated the unit effect using the Swamy & Arora (1972) method, and the time effect using the Wallace & Hussain (1969) method.

Hausman Test

The fixed effects model is commonly used in panel data analysis with desirable statistical properties. However, if the random effects model is more effective than the fixed effects model, it should be used instead. Therefore, it is necessary to determine which of the two models is more effective, since both are consistent but may differ in effectiveness. Several studies have used this effectiveness test, known as the Hausman test, which fits the Chi-squared distribution with k degrees of freedom, to determine which model is more effective between the fixed effects model and the random effects model (Baltagi, 2001Baltagi, B. Econometric analysis of panel data. 2nd ed. London: John Wiley & Sons; 2001.).

The null hypothesis states that the coefficients obtained from the random effects model and those obtained from the fixed effects model are not different. The rejection of the null hypothesis indicates that the random effects model is more effective than the fixed effects model.

RESULTS

The details of the parameters used in the study are presented in Table 1.

Table 1
The details of the parameters used in the study.

The mean mortality rate was found to be 0.57%, and the mean live weight loss due to transport was 0.30 kg (C/A-B)-D) (Table 1). The total waste amount was calculated as the sum of the decrease in production amount due to live weight loss and deaths during transport. The values calculated for mean live weight loss, mortality loss, and total waste due to transport are presented in Table 2.

Table 2
Mean live weight loss, mortality, and total waste due to transport.

As shown in Table 2, 95.37% of the total transport waste was attributed to the reduction in the live weight of broiler chickens (including excretion), while 4.62% was due to the loss of production caused by deaths.

We performed the Extended Augmented Dickey-Fuller (ADF) and Phillips-Perron (PP) tests to assess unit root (stationarity) in the dataset. The results of the unit root tests are presented in Table 3.

Table 3
The ADF and PP unit root test results of the variables in the Newey-West and Bartlett kernel panel data.

The panel data was stationary, as determined by the ADF and PP tests (p<0.05) (Table 3), and satisfied one of the requirements for panel data analysis.

In the panel data analysis, the parameters were estimated using fixed effects and random effects models to evaluate individual effects. The Hausman test was performed to determine the statistically valid model between the two, and the results are presented in Table 4.

Table 4
Panel data analysis of the relationship between mortality rate and the independent variables-Hausman Test.

The null hypothesis based on the Hausman test is a “random effects model” and the alternative hypothesis is a “fixed effects model”. The significance level (Probe) value was compared to the table value (α) (Table 3). As the Probe value was found to be 0.001 (i.e., p<0.050), the null hypothesis was rejected, indicating that a fixed-effects model was necessary for the panel data analysis. The estimation results obtained from the fixed-effects model are presented in Table 5.

Table 5
Fixed-effects panel data analysis of the relationship between mortality rate and several independent variables.

As shown in Table 5, the R2 value was 0.9093, indicating that the independent variables used in the model explained 90.93% of the changes in the mortality rate. Additionally, the F statistic showed that the model was significant.

The non-periodical fixed-effects panel data analysis showed a significant (p<0.01) effect between the number of living animals, number of dead animals, mean live weight at poultry farm exit, total live weight at poultry farm exit, mean live weight at slaughterhouse arrival, and a difference in the mortality rate.

The results of the number of transports and distances from poultry farms to slaughterhouses by province are presented in Table 6.

Table 6
The transportation distance from the poultry farm to the slaughterhouse by provinces.

As shown in Table 6, Bolu had the highest number of transports due to its proximity to the slaughterhouse, while Çankırı & Bilecik had the lowest number of transport as they were far from the slaughterhouse.

The independent variables included in the model explained 90.93% of the changes in the mortality rate. The estimation equations for estimating the mortality rate during transport from the poultry to the slaughterhouse in 11 provinces (Ankara, Bartın, Bilecik, Bolu, Çankırı, Düzce, Eskişehir, Karabük, Kocaeli, Sakarya, and Zonguldak), based on the model parameters, are presented in Table 7.

Table 7
Equations for estimating the mortality rate by provinces.

DISCUSSION

The transport of chickens from the broiler house to the poultry slaughterhouse incurs in financial losses, as it decreases live weight and increases the mortality rate. These losses negatively affect integrated companies aiming for high-quality, high-capacity, and sustainable production (Aral et al., 2014Aral Y, Arikan MS, Akin AC, et al. Economic losses due to live weight shrinkage and mortality during the broiler transport. Veterinary Journal of Ankara University 2014;61:205-10. https://doi.org/10.1501/Vetfak_0000002630
https://doi.org/10.1501/Vetfak_000000263...
).

Studies conducted in various countries have reported different mortality rates associated with the transport of broilers. For example, the mortality rate was reported to be 0.35% in Italy (Petracci et al., 2006Petracci M, Bianchi M, Cavani C, et al. Preslaughter mortality in broiler chickens, turkeys, and spent hens under commercial slaughtering. Poultry Science 2006;85:1660-4. https://doi.org/10.1093/ps/85.9.1660
https://doi.org/10.1093/ps/85.9.1660...
), 0.37% in the Czech Republic (Vecerek et al., 2016Vecerek V, Voslarova E, Conte F, et al. Negative trends in transport-related mortality rates in broiler chickens. Asian-Australasian Journal of Animal Sciences 2016;29(12):1796-804. https://doi.org/10.5713%2Fajas.15.0996
https://doi.org/10.5713%2Fajas.15.0996...
), and 0.41% and 0.46% in the Netherlands (Nijdam et al., 2004Nijdam E, Arens P, Lambooij E, et al. Factors influencing bruises and mortality of broilers during catching, transport and lairage. Poultry Science 2004;83:1610-5. https://doi.org/10.1093/ps/83.9.1610
https://doi.org/10.1093/ps/83.9.1610...
). The European Union (EU) has set the maximum recommended mortality rate at 0.5% (European Commission, 2005), while the USA reported a mortality rate of 0.68% (Ritz, 2005). In Turkey, studies have reported mortality rates of around 0.56%, which is between the values reported in the EU and the USA (Arikan et al., 2017Arikan MS, Akin AC, Akcay A, et al. Effects of Transportation Distance, SlaughterAge, and Seasonal Factors on Total Losses in Broiler Chickens. Brazilian Journal of Poultry Science 2017;19:421-7. https://doi.org/10.1501/Vetfak_0000002630
https://doi.org/10.1501/Vetfak_000000263...
). In this study, the mean mortality rate was found to be 0.57%, which was similar to the previously reported mortality rate in Turkey.

Lairage time, which is the period between the arrival of broiler chicks at the holding area in the slaughterhouse and their slaughter, helps animals adapt to their new environment, reduces thermal and physiological stress, and contributes to the welfare of broilers, thus improving meat quality (Hoffman & Lambrechts, 2011Hoffman LC, Lambrechts H. Bird handling, transportation, lairage and slaughter: implications for bird welfare and meat quality. In: Glatz P, Lunam C, Malecki IA, editors. The welfare of rarmed ratites. Heidelberg: Springer; 2011. p.195-235. ISBN-13: 978-3642192968). The recommended maximum waiting period is 1 h, although it might be extended to 2 h under favorable conditions (Warris et al., 1999). Exceeding this period, for example, by 4 h, significantly increases the mortality rate, as observed in previous studies (Bayliss & Hinton, 1990Bayliss PA, Hinton MH. Transportation of broilers with special reference to mortality rates. Applied Animal Behaviour Science 1990;28:93-118. https://doi.org/10.1016/0168-1591(90)90048-I
https://doi.org/10.1016/0168-1591(90)900...
). In this study, we found no significant relationship between lairage time and mortality rate; our findings differed from those of other studies (Bayliss & Hinton, 1990; Warris et al., 1999).

This discrepancy might be due to the slaughtering strategy of slaughterhouses, which prioritize herds with a low slaughter age rather than implementing a fixed waiting period in the holding area. The minimum lairage time recorded in our study was 0 min, indicating that some herds were immediately slaughtered upon arrival at the slaughterhouse.

In a study conducted in Poland in the winter, the effect of different transport distances on broiler mortality rates was examined. The mortality rates were 1.41% for a distance of 100 km, 2.65% for 200 km, and 2.36% for 300 km, compared to the control group that was not transported (Sowinska et al., 2013Sowinska J, Wojcik A, Pomianowski JF, et al. Effects of different variants of pre-slaughter transport on body weight loss and meat quality in broiler chickens. Medycyna Weterynaryjna 2013;69:420-3.). Another study in Brazil analyzed 13,937 transport operations and reported a mean mortality rate of 0.33% for distances between 24-242 km, with the highest rate (0.42%) being recorded during the summer (Vieira et al., 2015Vieira FMC, Deniz M, da Silva, IJO, et al. Pre-slaughter losses of broilers: effect of time period of the day and lairage time in a subtropical climate. Semina: Ciências Agrárias 2015;36(6):3887-96. https://doi.org/10.5433/1679-0359.2015v36n6p3887
https://doi.org/10.5433/1679-0359.2015v3...
). A study in the Czech Republic recorded a mean mortality rate of 0.247% during transport, with rates of 0.146% and 0.862% for enterprises within 50 km and over 300 km away, respectively (Vecerek et al., 2006Vecerek V, Grbalova S, Voslarova E, et al. Effects of travel distance and the season of the year on death rates in broilers transported to poultry processing plants. Poultry Science 2006;85:1881-4. https://doi.org/10.1093/ps/85.11.1881
https://doi.org/10.1093/ps/85.11.1881...
). In our study, the shortest and longest transport distances were 1.9 km and 276 km, with mortality rates of 0.45% and 1.72%, respectively.

We found a significant relationship between broiler mortality rates during transport from poultry to slaughterhouses and most variables, but not with lairage time (Table 5). The equations presented in Table 7 for each province can be used to identify critical control points that affect animal mortality and minimize losses during transport; thus, they can help reduce the mortality rate.

These equations were determined to estimate mortality rates during transport for each province, and can also be used in effective production planning for broiler integration and optimization of slaughterhouse capacity.

CONCLUSION

The mortality and live weight loss of broilers during transport from poultry farms to slaughterhouses were found to be affected by various factors, including geographical, climatic, and seasonal conditions; as well as production, slaughter, and transportation-related factors, such as transportation infrastructure, distance, time, and vehicles. Future econometric studies, paired with regular records on the factors affecting the mortality rate in poultry integrations, will reduce the mortality rate during poultry transfers.

In conclusion, studies like this should be conducted in different regions and provinces to optimize the costs of transportation, minimize mortality rates and live weight losses, optimize lairage times, select the shortest possible distances, and consider environmental and climatic conditions. Companies should consider these factors during the production planning stage.

ACKNOWLEDGMENT

There are no funding sources for this article.

REFERENCES

  • Aral Y, Arikan MS, Akin AC, et al. Economic losses due to live weight shrinkage and mortality during the broiler transport. Veterinary Journal of Ankara University 2014;61:205-10. https://doi.org/10.1501/Vetfak_0000002630
    » https://doi.org/10.1501/Vetfak_0000002630
  • Arikan MS, Akin AC, Akcay A, et al. Effects of Transportation Distance, SlaughterAge, and Seasonal Factors on Total Losses in Broiler Chickens. Brazilian Journal of Poultry Science 2017;19:421-7. https://doi.org/10.1501/Vetfak_0000002630
    » https://doi.org/10.1501/Vetfak_0000002630
  • Baltagi, B. Econometric analysis of panel data. 2nd ed. London: John Wiley & Sons; 2001.
  • Bayliss PA, Hinton MH. Transportation of broilers with special reference to mortality rates. Applied Animal Behaviour Science 1990;28:93-118. https://doi.org/10.1016/0168-1591(90)90048-I
    » https://doi.org/10.1016/0168-1591(90)90048-I
  • EC - European Commission. Health and Consumer Protection Directorate General. The welfare of chickens kept for meat production (Broilers) [SANCO.B.3/AH/R15/]. Brussel; 2005 Available from: https://ec.europa.eu/food/sites/food/files/safety/docs/sci-com_scah_out39_en.pdf
    » https://ec.europa.eu/food/sites/food/files/safety/docs/sci-com_scah_out39_en.pdf
  • EFSA - European Food Safety Authority. The welfare of animals during transport [EFSA-Q-2003-094]. Parma; 2004. Available from: https://www.efsa.europa.eu/en/efsajournal/pub/44
    » https://www.efsa.europa.eu/en/efsajournal/pub/44
  • Greenee, WH. Econometric analysis. 5th ed. New York: Prentice Hall; 2003.
  • Guerrero-Legarreta I. Handbook of poultry science and technology primary processing. London: John Wiley & Sons; 2010.
  • Hoffman LC, Lambrechts H. Bird handling, transportation, lairage and slaughter: implications for bird welfare and meat quality. In: Glatz P, Lunam C, Malecki IA, editors. The welfare of rarmed ratites. Heidelberg: Springer; 2011. p.195-235. ISBN-13: 978-3642192968
  • Im KS, Pesaran MH, Shin Y. Testing for unit roots in heterogeneous panels. Journal of Econometrics 2003;115:53-74. https://doi.org/10.1016/S0304-4076(03)00092-7
    » https://doi.org/10.1016/S0304-4076(03)00092-7
  • Kennedy, P. Ekonometri kilavuzu, (çev. Muzaffer Sarimeseli ve Senay Açikgöz), 5. Baski, Ankara: Gazi Kitabevi; 2006.
  • Lupo C, Le Bouquin S, Balaine L, et al. Chauvin C. Feasibility of screening broiler chicken flocks for risk markers as an aid for meat inspection. Epidemiology & Infection 2009; 137:1086-98. https://doi.org/10.1017/S095026880900209X
    » https://doi.org/10.1017/S095026880900209X
  • Nijdam E, Arens P, Lambooij E, et al. Factors influencing bruises and mortality of broilers during catching, transport and lairage. Poultry Science 2004;83:1610-5. https://doi.org/10.1093/ps/83.9.1610
    » https://doi.org/10.1093/ps/83.9.1610
  • Oba, A, Almeida M, Pinheiro JW, et al. The effect of management of transport and lairage conditions on broiler chicken breast meat quality and DOA (Death on Arrival). Brazilian Archives of Biology and Technology 2009;52:20511. https://doi.org/10.1590/S1516-89132009000700026
    » https://doi.org/10.1590/S1516-89132009000700026
  • Ondrasovicova O, Saba S, Smirkjakova S, et al. Effects of vehicle-road transport on blood profile in broiler chickens. Medycyna Weterynaryjna 2008;64:292-3.
  • Pazarlioglu, MV, Gürler ÖK. Telekomünikasyon Yatirimlari ve Ekonomik Büyüme: Panel Veri Yaklasimi. Finans Politik & Ekonomik Yorumlar 2007;44:35-43.
  • Petracci M, Bianchi M, Cavani C, et al. Preslaughter mortality in broiler chickens, turkeys, and spent hens under commercial slaughtering. Poultry Science 2006;85:1660-4. https://doi.org/10.1093/ps/85.9.1660
    » https://doi.org/10.1093/ps/85.9.1660
  • Pirompud P, Sivapirunthep P, Chaosap C. Factors influencing dead on arrival and condemnation of broilers during catching, transport, lairage and slaughter. International Journal of Agricultural Technology 2022;18(6):2535-44.
  • Santos VM dos, Dallago BS, Racanicci AM, et al. Effect of transportation distances, seasons and crate microclimate on broiler chicken production losses. PLoS One 2020;15(4);e0232004. https://doi.org/10.1371/journal.pone.0232004
    » https://doi.org/10.1371/journal.pone.0232004
  • Saraçoglu B, Nükhet D. Avrupa birligi ülkeleri ve avrupa birligine aday ülkelerin yakinsama analizi", VII. Ulusal ekonometri ve istatistik sempozyumu. Istanbul; 2005. Available From: http://www.ekonometridernegi.org/bildiriler/o10s1.pdf/30.06.2010
    » http://www.ekonometridernegi.org/bildiriler/o10s1.pdf/30.06.2010
  • Schwartzkopf-Genswein KS, Faucitano L, Dadgar S, et al. Road transport of cattle, swine and poultry in North America and its impact on animal welfare, carcass and meat quality: a review. Meat Science 2012;92:227-43. https://doi.org/10.1016/j.meatsci.2012.04.010
    » https://doi.org/10.1016/j.meatsci.2012.04.010
  • Sowinska J, Wojcik A, Pomianowski JF, et al. Effects of different variants of pre-slaughter transport on body weight loss and meat quality in broiler chickens. Medycyna Weterynaryjna 2013;69:420-3.
  • Swamy PAVB, Arora SS. The Exact Finite Sample Properties of the Esitmator of Coefficients in the Error Components Regression Model. Econometrica 1972;40(2):261-75. https://doi.org/10.2307/1909405
    » https://doi.org/10.2307/1909405
  • Vecerek V, Grbalova S, Voslarova E, et al. Effects of travel distance and the season of the year on death rates in broilers transported to poultry processing plants. Poultry Science 2006;85:1881-4. https://doi.org/10.1093/ps/85.11.1881
    » https://doi.org/10.1093/ps/85.11.1881
  • Vecerek V, Voslarova E, Conte F, et al. Negative trends in transport-related mortality rates in broiler chickens. Asian-Australasian Journal of Animal Sciences 2016;29(12):1796-804. https://doi.org/10.5713%2Fajas.15.0996
    » https://doi.org/10.5713%2Fajas.15.0996
  • Vieira FMC, Deniz M, da Silva, IJO, et al. Pre-slaughter losses of broilers: effect of time period of the day and lairage time in a subtropical climate. Semina: Ciências Agrárias 2015;36(6):3887-96. https://doi.org/10.5433/1679-0359.2015v36n6p3887
    » https://doi.org/10.5433/1679-0359.2015v36n6p3887
  • Vieira FMC, Silva IJOD, Barbosa Filho JAD, et al. Productive losses on broiler preslaughter operations: effects of the distance from farms to abattoirs and of lairage time in a climatized holding area. Revista Brasileira de Zootecnia 2010;39:2471-6. https://doi.org/10.1590/S1516-35982010001100021
    » https://doi.org/10.1590/S1516-35982010001100021
  • Wallace TD, Hussain A. The use of error components models in combining cross section with time series data. Econometrica: Journal of the Econometric Society 1969;37(1);55-72. https://doi.org/10.2307/1909205
    » https://doi.org/10.2307/1909205
  • Warriss PD, Bevis EA, Brown SN, et al. Longer journeys to processing plants are associated with higher mortality in broiler chickens. British Poultry Science 1992;33:201-6. https://doi.org/10.1080/00071669208417458
    » https://doi.org/10.1080/00071669208417458
  • Warriss PD, Knowles TG, Brown SN, et al. Effects of lairage time on body temperature and glycogen reserves of broiler chickens held in transport modules. Veterinary Record 1999;145.8:218-22. https://doi.org/10.1136/vr.145.8.218
    » https://doi.org/10.1136/vr.145.8.218
  • Warriss PD. Meat science: an introductory text . London: CAB International; 2010.

Publication Dates

  • Publication in this collection
    13 Oct 2023
  • Date of issue
    2023

History

  • Received
    08 May 2023
  • Accepted
    13 July 2023
Fundação de Apoio à Ciência e Tecnologia Avicolas Rua Barão de Paranapanema, 146 - Sala 72, Bloco A, Bosque, Campinas, SP - 13026-010. Tel.: 19 3255-8500 - Campinas - SP - Brazil
E-mail: revista@facta.org.br