# ABSTRACT:

This paper aimed to analyze the evolution of the quality of raw milk produced by producers in the Zona da Mata region, in the state of Minas Gerais, between 2012 and 2018. For this purpose, we used the linear mixed-effects model to analyze the monthly evolution of the results of milk composition indicators (fat, protein, and defatted dry extract - DDE), somatic cell count - SCC, and total bacterial count - TBC, from the official monthly registry of 94 milk producers, suppliers of a region’s dairy. Results indicate a continuous reduction in the milk composition indicators between 2014 and 2018. For the SCC and TBC indicators, we identified only one-off reductions. The supply of a larger volume of milk was associated with increased TBC. Seasonality influenced all quality indicators analyzed. According to these results, we concluded that the analyzed quality of chilled raw milk offered by producers was proven worse from 2012 to 2018, despite the actions implemented by the National Milk Quality Improvement Program - PNMQL.

Key words:
dairy; legislation; National Milk Quality Improvement Program

# RESUMO:

Palavras-chave:
laticínio; legislação; Programa Nacional de Melhoria de Qualidade do Leite

# INTRODUCTION:

Nevertheless, two decades after the implementation of PNMQL, it cannot be said that the chilled raw milk that reaches the industry is of higher quality. RIBAS et al. (2016RIBAS, N. P. et al. Total bacterial count in milk samples from bulk tank Paraná State. Archives of Veterinary Science, v.21, n.1, p.32-43, 2016. Available from: <Available from: https://revistas.ufpr.br/veterinary/article/view/41581 >. Accessed: Dec. 22, 2018. doi: 10.5380/avs.v21i1.41581.
https://revistas.ufpr.br/veterinary/arti...
), in a study with milk samples from 10 regions of Paraná State, identified a reduction in the TBC indicator between January 2012 and May 2014. However, in this and several other studies, researchers reported a high percentage of samples in non-compliance with current SCC and TBC standards (BORGES et al., 2013BORGES, L. R. et al. Milk quality according to the daily range in farm production in the Mesoregion Central Mineira and Oeste of Minas Gerais regions, Brazil. Arquivo Brasileiro de Medicina Veterinária e Zootecnia, v.65, n.4, p.1239-1246, 2013. Available from: <Available from: http://www.scielo.br/scielo.php?script=sci_abstract&pid=S0102-09352013000400038&lng=en&nrm=iso&tlng=en >. Accessed: Sep. 12, 2019. doi: 10.1590/S0102-09352013000400041.
http://www.scielo.br/scielo.php?script=s...
; MILANI et al., 2016MILANI, M. P. et al. Milk quality in different production systems, year and weather seasons. Revista Brasileira de Ciência Veterinária, v.23, n.3-4, p.206-211, 2016. Available from: <Available from: http://doi.editoracubo.com.br/10.4322/rbcv.2016.058 . Accessed: Jan. 19, 2018. doi: 10.4322/rbcv.2016.058.
http://doi.editoracubo.com.br/10.4322/rb...
; RIBAS et al., 2016RIBAS, N. P. et al. Total bacterial count in milk samples from bulk tank Paraná State. Archives of Veterinary Science, v.21, n.1, p.32-43, 2016. Available from: <Available from: https://revistas.ufpr.br/veterinary/article/view/41581 >. Accessed: Dec. 22, 2018. doi: 10.5380/avs.v21i1.41581.
https://revistas.ufpr.br/veterinary/arti...
).

Coupled with the scenario of uncertainty regarding the implementation agenda of the new standards, studies evaluating the quality of refrigerated raw milk in Brazil over the years are still scarce. However, these studies are extremely important for the sector, as they allow, among other things, the evolution of the quality of milk produced in different regions of the country to be evaluated, the factors that may be associated with the production of quality milk to be identified, and future reviews of the milk legal framework to be supported. Along these lines, this article aimed to determine whether the monthly quality of the results of composition indicators (fat, protein, and defatted dry extract - DDE), SCC, and TBC of chilled raw milk produced by producers in the Zona da Mata region, in the state of Minas Gerais, evolved between 2012 and 2018.

# MATERIALS AND METHODS:

## Data

The data come from official monthly analyses of 94 dairy farmers (individual and collective), suppliers of a dairy located in the Zona da Mata region, state of Minas Gerais. The data cover a period of 78 months from July 2012 to December 2018.

The variables considered in the study were monthly production volume (Volume), year (Year; ranging from 2012 to 2018), and seasonality (Seasonality; as a function of month), referring to analyses of fat (Fat), protein (Protein), degreased dry extract (DDE), SCC and TBC, in addition to the temperature of the milk sample analyzed (T_SCC for SCC and composition analyzes and T_TBC for TBC analysis) and the number of days elapsed between milk sample collection and analysis (Age). Values ​​of (i) SCC zero or greater than 9,999 × 103 cells/mL; (ii) TBC equal to zero or greater than 9,999 × 103 cfu/mL; (iii) Fat less than 0.5% or greater than 7.0%; (iv) Protein less than 1.0% or greater than 5.0% (OSTRENSKY, 1999OSTRENSKY, A. Efeitos de ambiente sobre a contagem de células somáticas no leite de vacas da raça holandesa no Paraná. 1999, 114p. Dissertação (Mestrado em Ciências Veterinárias) - Universidade Federal do Paraná, Curitiba, 1999. Available from: <Available from: https://acervodigital.ufpr.br/handle/1884/28158 >. Accessed: Sep. 12, 2019.
https://acervodigital.ufpr.br/handle/188...
); (v) T_SCC (or T_TBC) greater than 10 °C; and (vi) Age greater than 7 days (CASSOLI, 2005CASSOLI, L. D. Evaluation of flow cytometry as a method for total bacterial count of raw milk. 2005, 46p. Dissertação (Mestrado em Agronomia) - Escola Superior de Agricultura “Luiz de Queiroz”, Piracicaba, 2005. Available from: <Available from: https://www.teses.usp.br/teses/disponiveis/11/11139/tde-26092005-164035/publico/LaerteCassoli.pdf >. Accessed: Sep. 12, 2019.
https://www.teses.usp.br/teses/disponive...
) were removed from the analysis. We also removed the outliers referring to the monthly volume of milk supplied by producer, due to inconsistent data. In this case, we considered outliers those values ​​lower than the expression: first quartile minus 1.5 times the interquartile distance, and those higher than the expression: the third quartile plus 1.5 times the interquartile distance.

For months with more than one official analysis record per producer, we calculated the arithmetic mean for composition (Fat, Protein, and DDE), temperature (T_SCC and T_TBC), and Age analyzes, and geometric mean for SCC and TBC analyzes. In these cases, the use of geometric mean was due to the lognormal behavior of these data (SHOOK, 1982SHOOK, G. E. Relationship between lactation measures of somatic cell concentration and milk yield. Journal of Dairy Science , v.65, n.3, p.419-425, 1982. Available from: <Available from: https://www.sciencedirect.com/science/article/pii/S0022030282822078 >. Accessed: Jan. 05, 2019. doi: 10.3168/jds.S0022-0302(82)82207-8.
https://www.sciencedirect.com/science/ar...
). Also noteworthy it was the use of geometric mean for these variables when there were results from more than one monthly analysis, according to Normative Instruction No. 76 (BRASIL, 2018bBRASIL. Instrução Normativa n.77, de 26 de novembro de 2018b. Estabelece os critérios e procedimentos para a produção, acondicionamento, conservação, transporte, seleção e recepção do leite cru em estabelecimentos registrados no serviço de inspeção oficial, na forma desta Instrução Normativa e do seu Anexo. Diário Oficial da República Federativa do Brasil, Brasília, 30 de novembro de 2018. Seção 1, p.10.).

## Statistical analysis

Initially, we calculated the descriptive statistics of the data. Then, we used the linear mixed-effects model to identify the effect of the variables Producer, Volume, Year (2013, 2014, 2015, 2016, 2017, and 2018) and Seasonality (Station I and Station II), under the indicators of quality of chilled raw milk. In longitudinal studies, which considered repeated measures over time, the use of a linear mixed-effects model is more appropriate than Analysis of Variance - ANOVA. Using ANOVA for longitudinal studies ignores the correlation between observations and treats them as uncorrelated (FAUSTO et al., 2008FAUSTO, M. A. et al. Mixed linear regression model for longitudinal data: application to an unbalanced anthropometric data set. Cadernos de Saúde Pública, v.24, n.3, p.513-524, 2008. Available from: <Available from: http://www.scielo.br/scielo.php?script=sci_arttext&pid=S0102-311X2008000300005&lng=pt&tlng=pt >. Accessed: Sep. 12, 2019. doi: 10.1590/S0102-311X2008000300005.
http://www.scielo.br/scielo.php?script=s...
).

The linear mixed-effects model includes fixed and random effects and is commonly used in longitudinal studies in the dairy field (NIGHTINGALE et al., 2008NIGHTINGALE, C. et al. Influence of variable milk quality premiums on observed milk quality. Journal of Dairy Science , v.91, p.1236-1244, 2008. Available from: <Available from: https://www.sciencedirect.com/science/article/pii/S0022030208713821 >. Accessed: Jan, 18, 2018. doi: 10.3168/jds.2007-0609.
https://www.sciencedirect.com/science/ar...
; BOTARO et al., 2013BOTARO, B. G. et al. Quality based payment program and milk quality. Scientia Agricola, v.70, n.1, p.21-26, 2013. Available from: <Available from: http://www.scielo.br/scielo.php?script=sci_arttext&pid=S0103-90162013000100004 >. Accessed: Sep. 08, 2019. doi: 10.1590/S0103-90162013000100004.
http://www.scielo.br/scielo.php?script=s...
; FLORES-MIYAMOTO et al., 2014FLORES-MIYAMOTO, A. et al. Do farm audits improve milk quality? Journal of Dairy Science, v.97, n.1, p.1-9, 2014. Available from: <Available from: https://www.sciencedirect.com/science/article/pii/S0022030213007303 >. Accessed: Sep. 12, 2019. doi: 10.3168/jds.2012-6228.
https://www.sciencedirect.com/science/ar...
; O’CONNELL et al., 2015O’CONNELL, A. et al. Seasonal trends in milk quality in Ireland between 2007 and 2011. Journal of Dairy Science , v.98, p.3778-3790, 2015. Available from: <Available from: https://www.sciencedirect.com/science/article/pii/S002203021500212X >. Accessed: Jan. 18, 2018. doi: 10.3168/jds.2014-9001.
https://www.sciencedirect.com/science/ar...
; MILANI et al., 2016MILANI, M. P. et al. Milk quality in different production systems, year and weather seasons. Revista Brasileira de Ciência Veterinária, v.23, n.3-4, p.206-211, 2016. Available from: <Available from: http://doi.editoracubo.com.br/10.4322/rbcv.2016.058 . Accessed: Jan. 19, 2018. doi: 10.4322/rbcv.2016.058.
http://doi.editoracubo.com.br/10.4322/rb...
). In this, fixed effects are directly estimated, being analogous to the coefficients of a standard regression, while random effects are summarized according to their estimated variances and covariance (STATACORP, 2013STATACORP. Stata multilevel mixed-effects reference manual (release 13). College Station, TX: StataCorp LP. 2013. 364p. Available from: <Available from: https://www.stata.com/manuals13/me.pdf >. Accessed: Sep. 12, 2019.
https://www.stata.com/manuals13/me.pdf...
). The linear mixed-effects model can be represented by the following equation:

$y=\mathit{X.\beta }+\mathit{Z.u}+є$

where y is the known vector of observations for a particular quality indicator of chilled raw milk; β is the unknown vector of fixed effects; u is the unknown vector of random effects; ϵ is the waste vector; X and Z are known matrices, related to observations y for β and u, respectively. It is assumed that Z has a normal distribution with mean zero and variance and covariance G, regardless of ϵ , which has a normal distribution with mean zero and variance ${\sigma }^{2}$.

In this study, we analyzed the variables Year, Volume, and Seasonality as fixed effects, while the control variable Producer was analyzed as a random effect. For the Year variable, we broke down the analysis year by year by using dummy variables for each of the years. The variables Fat, Protein, DDE, SCC, and TBC were analyzed as dependent terms, each being in a distinct model. The model used in this study is represented in the equation below:

$\mathit{Indicator}=\mathit{Intercept}+\mathit{Seasonality}+\mathit{Year}+\mathit{Volume}+\mathit{Producer}\left(\mathit{random}\right)+\mathit{Error}$ Indicator = Intercept + Seasonality + Year + Volume + Producer (random) + Error>

where Seasonality is a regression term for the sine (Season I) and cosine (Season II) functions, according to NIGHTINGALE et al. (2008NIGHTINGALE, C. et al. Influence of variable milk quality premiums on observed milk quality. Journal of Dairy Science , v.91, p.1236-1244, 2008. Available from: <Available from: https://www.sciencedirect.com/science/article/pii/S0022030208713821 >. Accessed: Jan, 18, 2018. doi: 10.3168/jds.2007-0609.
https://www.sciencedirect.com/science/ar...
), which thus describe the monthly variation of the evaluated quality indicator:

$\mathit{Season I}=\mathit{sin}\left(\frac{2\mathit{.\pi .Month}}{12}\right)$

$\mathit{Season II}=\mathit{cos}\left(\frac{2\mathit{.\pi .Month}}{12}\right)$

Since the values of SCC and TBC are typically lognormal (SHOOK, 1982SHOOK, G. E. Relationship between lactation measures of somatic cell concentration and milk yield. Journal of Dairy Science , v.65, n.3, p.419-425, 1982. Available from: <Available from: https://www.sciencedirect.com/science/article/pii/S0022030282822078 >. Accessed: Jan. 05, 2019. doi: 10.3168/jds.S0022-0302(82)82207-8.
https://www.sciencedirect.com/science/ar...
), they were log transformed (base 10) before statistical analyses (Log_SCC and Log_TBC, respectively). The reference unit considered for the Volume variable was 1,000 L of milk. We estimated the fixed effects and variance components of the linear mixed-effects model by the Restricted Maximum Likelihood method - REML (PATTERSON and THOMPSON, 1971PATTERSON, H. D.; THOMPSON, R. Recovery of inter-block information when block sizes are unequal. Biometrika, v.58, n.3, p.545-554, 1971. Available from: <Available from: https://www.mendeley.com/catalogue/recovery-interblock-information-block-sizes-unequal-1/ >. Accessed: Jan. 13, 2019. doi: 10.2307/2334389.
https://www.mendeley.com/catalogue/recov...
). This method provided non-biased estimates of variance, removing the addiction that exists in the estimation of the mean (RESENDE et al., 1996RESENDE, M. D. V. et al. Variance component estimation and prediction of breeding values by the restricted maximum likelihood (REML) and best linear unbiased prediction (BLUP) methods in Pinus. Boletim de Pesquisa Florestal, n.32/33, p.23-42, 1996. Available from: <Available from: https://www.alice.cnptia.embrapa.br/bitstream/doc/282154/1/mresende1.pdf >. Accessed: Sep. 12, 2019.
https://www.alice.cnptia.embrapa.br/bits...
).

# RESULTS AND DISCUSSION:

## Descriptive statistics

Table 1
Descriptive statistics of the raw milk variables of 94 producers analyzed.

The minimum values observed for Fat, Protein, and DDE, as well as the maximum values observed for SCC and TBC, indicated the existence of cases that do not comply with the legislation. The mean values of T_SCC, T_TBC, and Age are in accordance with the suggested temperature and sample age limit for the analysis of chilled raw milk (CASSOLI, 2005CASSOLI, L. D. Evaluation of flow cytometry as a method for total bacterial count of raw milk. 2005, 46p. Dissertação (Mestrado em Agronomia) - Escola Superior de Agricultura “Luiz de Queiroz”, Piracicaba, 2005. Available from: <Available from: https://www.teses.usp.br/teses/disponiveis/11/11139/tde-26092005-164035/publico/LaerteCassoli.pdf >. Accessed: Sep. 12, 2019.
https://www.teses.usp.br/teses/disponive...
). The monthly volume of milk supplied by producers was the variable with the highest standard deviation, with an average volume of 2,520 L. The evolution of the historical average of the quality indicators studied for chilled raw milk is illustrated in figure 1.

Figure 1
Monthly history of quality indicators for chilled raw milk between 2012 and 2018 Fat (a), Protein (b), DDE (c), SCC (d) and TBC (e).

From the monthly history of milk quality indicators over the period under study (Figure 1), it is possible to verify the existence of large oscillations, which represents a problem for the dairy industry. The content of milk composition indicators is closely related to the yield of various dairy products (MAGANHA, 2006MAGANHA, M. F. B. Guia técnico ambiental da indústria de produtos lácteos. São Paulo: CETESB, 2006. 95p. Available from: <Available from: http://www.cetesb.sp.gov.br >. Accessed: Jan. 1, 2015.
http://www.cetesb.sp.gov.br...
), while the SCC and TBC indicators are not only related to the yield of various dairy products (HAND et al., 2012HAND, K. J. et al. Milk production and somatic cell counts: A cow-level analysis. Journal of Dairy Science, v.95, p.1358-1362, 2012. Available from: <Available from: https://www.sciencedirect.com/science/article/pii/S0022030212000926 >. Accessed: Sep. 12, 2019. doi: 10.3168/jds.2011-4927.
https://www.sciencedirect.com/science/ar...
) but also to its validity (BARBANO et al., 2006BARBANO, D. M. et al. Influence of raw milk quality on fluid milk shelf life. Journal of Dairy Science, v.89, p.15-19, 2006. doi: 10.3168/jds.S0022-0302(06)72360-8.
https://doi.org/10.3168/jds.S0022-0302(0...
; FLORES-MIYAMOTO et al., 2014FLORES-MIYAMOTO, A. et al. Do farm audits improve milk quality? Journal of Dairy Science, v.97, n.1, p.1-9, 2014. Available from: <Available from: https://www.sciencedirect.com/science/article/pii/S0022030213007303 >. Accessed: Sep. 12, 2019. doi: 10.3168/jds.2012-6228.
https://www.sciencedirect.com/science/ar...
; CASTRO et al., 2016CASTRO, R. D. et al. Lactic acid microbiota identification in water, raw milk, endogenous starter culture, and fresh Minas artisanal cheese from the Campo das Vertentes region of Brazil during the dry and rainy seasons. Journal of Dairy Science, v.99, p.6086-6096, 2016. Available from: <Available from: https://www.sciencedirect.com/science/article/pii/S0022030216303502 >. Accessed: Sep. 12, 2019. doi: 10.3168/jds.2015-10579.
https://www.sciencedirect.com/science/ar...
) and their sensory characteristics (MA et al., 2000MA, Y. et al. Effects of somatic cell count on quality and shelf-life of pasteurized fluid milk. Journal of Dairy Science , v.83, p.264-274, 2000. Available from: <Available from: https://www.journalofdairyscience.org/article/S0022-0302(00)74873-9/abstract >. Accessed: Sep. 12, 2019. doi: 10.3168/jds.S0022-0302(00)74873-9.
https://www.journalofdairyscience.org/ar...
; BARBANO et al., 2006BRASIL. Instrução Normativa n.31, de 29 de junho de 2018a. Diário Oficial da República Federativa do Brasil, Brasília, 02 de julho de 2018. Seção 1, p.2.). Thus, the observed oscillations make it difficult to standardize production processes and the quality of the products produced.

Among the indicators of milk composition, it is possible to identify the presence of cyclic oscillations, with a slight downward trend over the years. In this case, the oscillations may be associated with the influence of the animal’s diet on the quality of milk produced (JENKINS and McGUIRE, 2006JENKINS, T. C.; McGUIRE, M. A. Major Advances in Nutrition: Impact on Milk Composition. Journal of Dairy Science , v.89, p.1302-1310, 2006. Available from: <Available from: https://www.sciencedirect.com/science/article/pii/S0022030206721981 >. Accessed: Sep. 13, 2019. doi: 10.3168/jds.S0022-0302(06)72198-1.
https://www.sciencedirect.com/science/ar...
). When it comes to grazing animals, climate change can influence the availability and quality of their food (MILANI et al., 2016MILANI, M. P. et al. Milk quality in different production systems, year and weather seasons. Revista Brasileira de Ciência Veterinária, v.23, n.3-4, p.206-211, 2016. Available from: <Available from: http://doi.editoracubo.com.br/10.4322/rbcv.2016.058 . Accessed: Jan. 19, 2018. doi: 10.4322/rbcv.2016.058.
http://doi.editoracubo.com.br/10.4322/rb...
).

The SCC and TBC indicators, in turn, did not show such a well-defined behavior. However, it was possible to verify the presence of larger oscillations of the TBC indicator, when compared with the SCC indicator (Figure 1), which was corroborated by the coefficient of variation of the logarithm of these indicators (Table 1). This result may be associated with the fact that SCC is directly associated with the health status of the cow’s mammary gland (MACHADO et al., 2000MACHADO, P. F. et al. Somatic cells in milk of Brazilian herds. Scientia Agricola, v.57, n.2, p.359-361, 2000. Available from: <Available from: http://www.scielo.br/scielo.php?script=sci_arttext&pid=S0103-90162000000200026&lng=pt&tlng=pt >. Accessed: Sep. 12, 2019. doi: 10.1590/S0103-90162000000200026.
http://www.scielo.br/scielo.php?script=s...
), while TBC is associated, among other factors, with hygienic-sanitary milking conditions (PANTOJA et al., 2009PANTOJA, J. C. et al. Associations among milk quality indicators in raw bulk milk. Journal of Dairy Science , v.92, p.4978-4987, 2009. Available from: <Available from: https://www.sciencedirect.com/science/article/pii/S002203020970829X >. Accessed: Sep. 13, 2019. doi: 10.3168/jds.2009-2329.
https://www.sciencedirect.com/science/ar...
). Therefore, the latter tends to be more impacted by climate change, since, in rainy periods, the maintenance of hygienic milking conditions tends to be impaired (BUENO et al., 2008BUENO, V. F. F. et al. Total bacterial count: relationship to milk composition and period of the year in Goiás State, Brazil. Revista Brasileira de Ciência Veterinária, v.5, n.1, p.40-44, 2008. Available from: <http://doi.editoracubo.com.br/10.4322/rbcv.2014.194>. Accessed: Sep. 12, 2019.doi: 10.4322/rbcv.2014.194.
https://doi.org/10.4322/rbcv.2014.194....
; CASTRO et al., 2016CASTRO, R. D. et al. Lactic acid microbiota identification in water, raw milk, endogenous starter culture, and fresh Minas artisanal cheese from the Campo das Vertentes region of Brazil during the dry and rainy seasons. Journal of Dairy Science, v.99, p.6086-6096, 2016. Available from: <Available from: https://www.sciencedirect.com/science/article/pii/S0022030216303502 >. Accessed: Sep. 12, 2019. doi: 10.3168/jds.2015-10579.
https://www.sciencedirect.com/science/ar...
). In a study conducted with data from 16,491 milk samples analyzed at the Milk Quality Laboratory of the Food Research Center of the Veterinary School of the Universidade Federal de Goiás, between October 2002 and September 2003, BUENO et al. (2008)BUENO, V. F. F. et al. Total bacterial count: relationship to milk composition and period of the year in Goiás State, Brazil. Revista Brasileira de Ciência Veterinária, v.5, n.1, p.40-44, 2008. Available from: <http://doi.editoracubo.com.br/10.4322/rbcv.2014.194>. Accessed: Sep. 12, 2019.doi: 10.4322/rbcv.2014.194.
https://doi.org/10.4322/rbcv.2014.194....
reported higher values ​​for the logarithm of the TBC indicator in the rainy season when compared with the dry season.

In any case, the evolution trend of the quality indicators of chilled raw milk is not always clear, and the oscillations may be associated with a set of factors, such as seasonal variations, producer entry or exit, change in technical assistance policy by the dairy, etc. (NIGHTINGALE et al., 2008NIGHTINGALE, C. et al. Influence of variable milk quality premiums on observed milk quality. Journal of Dairy Science , v.91, p.1236-1244, 2008. Available from: <Available from: https://www.sciencedirect.com/science/article/pii/S0022030208713821 >. Accessed: Jan, 18, 2018. doi: 10.3168/jds.2007-0609.
https://www.sciencedirect.com/science/ar...
; BOTARO et al., 2013BOTARO, B. G. et al. Quality based payment program and milk quality. Scientia Agricola, v.70, n.1, p.21-26, 2013. Available from: <Available from: http://www.scielo.br/scielo.php?script=sci_arttext&pid=S0103-90162013000100004 >. Accessed: Sep. 08, 2019. doi: 10.1590/S0103-90162013000100004.
http://www.scielo.br/scielo.php?script=s...
; FLORES-MIYAMOTO et al., 2014FLORES-MIYAMOTO, A. et al. Do farm audits improve milk quality? Journal of Dairy Science, v.97, n.1, p.1-9, 2014. Available from: <Available from: https://www.sciencedirect.com/science/article/pii/S0022030213007303 >. Accessed: Sep. 12, 2019. doi: 10.3168/jds.2012-6228.
https://www.sciencedirect.com/science/ar...
). Therefore, the use of statistical models, such as linear mixed-effects models, is necessary to measure the specific contribution of each variable under the evolution observed for the indicators of milk quality.

## Evolution of chilled raw milk quality

Table 2 presents the statistics for the mixed-effects linear models generated for the quality indicators of chilled raw milk. The average number of observations recorded per producer was 53, ranging from 30 to 75, while the average number of observations recorded per model was 4,488, ranging from 4,400 to 4,514. All models were significant (P<0.0001) according to the Wald test and significantly higher than the reference linear regression model (P<0.0001) according to the likelihood ratio test. The Year and Seasonality variables were significant (P<0.05) in all models, while the Volume variable was significant only for the TBC indicator model.

Table 2
Coefficients of linear mixed effects models, referring to the variables analyzed.

The years between 2014 and 2018 were associated with the reduction in the indicators of milk composition (Table 2). Even so, the monthly averages for Fat, Protein, and DDE levels were always higher than the minimum values established for these indicators (3.0 g/100 mL, 2.9 g/100 mL, and 8.4 g/100 mL, respectively). The only exception was made in September 2016, when the monthly average DDE was 8.3 g/100 mL, lower than the minimum value established by current legislation (BRASIL, 2011bBRASIL. Instrução Normativa n.62, de 29 de dezembro de 2011b. Aprova Regulamentos Técnicos de Produção, Identidade e Qualidade do Leite tipo A, Leite Cru Refrigerado, Leite Pasteurizado e o Regulamento Técnico da Coleta de Leite Cru Refrigerado e seu Transporte a Granel. Diário Oficial da República Federativa do Brasil, Brasília, 30 de dezembro de 2011. Seção 1, p.6-11.; BRASIL, 2018bBRASIL. Instrução Normativa n.77, de 26 de novembro de 2018b. Estabelece os critérios e procedimentos para a produção, acondicionamento, conservação, transporte, seleção e recepção do leite cru em estabelecimentos registrados no serviço de inspeção oficial, na forma desta Instrução Normativa e do seu Anexo. Diário Oficial da República Federativa do Brasil, Brasília, 30 de novembro de 2018. Seção 1, p.10.). It is still necessary to emphasize the positive correlation between the milk solids content and the productivity that the dairy will have when converting this raw material into dairy derivatives (MAGANHA, 2006MAGANHA, M. F. B. Guia técnico ambiental da indústria de produtos lácteos. São Paulo: CETESB, 2006. 95p. Available from: <Available from: http://www.cetesb.sp.gov.br >. Accessed: Jan. 1, 2015.
http://www.cetesb.sp.gov.br...
). Thus, the observed reduction in milk solids content between 2014 and 2018 is a bad result, especially in light of the PNMQL objectives.

Most of Brazil’s milk production comes from small producers (IBGE, 2019INSTITUTO BRASILEIRO DE GEOGRAFIA E ESTATÍSTICA (IBGE). Censo Agropecuário 2017 - Resultados Preliminares. 2019. Available from: <Available from: https://ibge.gov.br/ >. Accessed: Sep. 13, 2019.
https://ibge.gov.br/...
). These, in general, do not have technical assistance (OLIVAL et al., 2004OLIVAL, A. A. et al. Evaluation of limitations for improvement of milk quality in the area of Pirassununga-SP. Revista Ciência em Extensão, v.1, n.2, p.171-183, 2004. Available from: <Available from: https://ojs.unesp.br/index.php/revista_proex/article/view/173 >. Accessed: Sep. 13, 2019.
https://ojs.unesp.br/index.php/revista_p...
; OLIVEIRA et al., 2011OLIVEIRA, C. J. B. et al. Risk factors associated with selected indicators of milk quality in semiarid northeastern Brazil. Journal of Dairy Science , v.94, p.3166-3175, 2011. Available from: <Available from: https://www.sciencedirect.com/science/article/pii/S0022030211003043 >. Accessed: Sep. 13, 2019. doi: 10.3168/jds.2010-3471.
https://www.sciencedirect.com/science/ar...
), are unaware of ways to prevent milk contamination (OLIVAL et al., 2004OLIVAL, A. A. et al. Evaluation of limitations for improvement of milk quality in the area of Pirassununga-SP. Revista Ciência em Extensão, v.1, n.2, p.171-183, 2004. Available from: <Available from: https://ojs.unesp.br/index.php/revista_proex/article/view/173 >. Accessed: Sep. 13, 2019.
https://ojs.unesp.br/index.php/revista_p...
), are poorly capitalized (OLIVEIRA et al., 2011OLIVEIRA, C. J. B. et al. Risk factors associated with selected indicators of milk quality in semiarid northeastern Brazil. Journal of Dairy Science , v.94, p.3166-3175, 2011. Available from: <Available from: https://www.sciencedirect.com/science/article/pii/S0022030211003043 >. Accessed: Sep. 13, 2019. doi: 10.3168/jds.2010-3471.
https://www.sciencedirect.com/science/ar...
), have a low profit margin under the price of a liter of milk (BÁNKUTI et al., 2008BÁNKUTI, F. I. et al. Measurement and analysis of transaction costs for milk producers in formal and informal markets in the region of São Carlos, São Paulo state, Brazil. Organizações Rurais & Agroindustriais, v.10, n.3, p.343-358, 2008. Available from: <Available from: http://200.131.250.22/revistadae/index.php/ora/article/view/70/66 >. Accessed: Sep. 14, 2019.
; LOPES et al., 2012LOPES, M. A. et al. Comparative economic indicators of milk production systems with high daily volume in Minas Gerais State. Revista Ceres, v.59, n.4, p.458-465, 2012. Available from: <Available from: http://www.ceres.ufv.br/ojs/index.php/ceres/article/view/3874/1714 >. Accessed: Sep. 14, 2019.
http://www.ceres.ufv.br/ojs/index.php/ce...
; 2015LOPES, M. A. et al. Effect of technical and managable index in profitability of dairy production with different technological levels. Revista Científica de Produção Animal, v.17, n.2, p.92-102, 2015. Available from: <Available from: https://periodicos.ufpb.br/index.php/rcpa/article/view/42772 >. Accessed: Sep. 12, 2019. doi: 10.15528/2176-4158/rcpa.v17n2p92-102.
https://periodicos.ufpb.br/index.php/rcp...
), and are not financially rewarded for milk quality (VALEEVA et al., 2007VALEEVA, N. I. et al. Motivation of Dairy Farmers to Improve Mastitis Management. Journal of Dairy Science , v.90, n.9, p.4466-4477, 2007. Available from: <Available from: https://www.sciencedirect.com/science/article/pii/S0022030207719094 >. Accessed: Jan. 05, 2019. doi: 10.3168/jds.2007-0095.
https://www.sciencedirect.com/science/ar...
). These are, therefore, obstacles to improving milk quality and as such need to be targeted by PNMQL actions.

It should be noted that the herd profile also has an influence on milk solid content (MILANI et al., 2016MILANI, M. P. et al. Milk quality in different production systems, year and weather seasons. Revista Brasileira de Ciência Veterinária, v.23, n.3-4, p.206-211, 2016. Available from: <Available from: http://doi.editoracubo.com.br/10.4322/rbcv.2016.058 . Accessed: Jan. 19, 2018. doi: 10.4322/rbcv.2016.058.
http://doi.editoracubo.com.br/10.4322/rb...
). Thus, these results may be associated with a change in the herd profile over the period analyzed. However, this inference cannot be ratified due to the unavailability of monthly production data for each of the producers.

For the SCC indicator, there was only a reduction of 0.08 log10 (cells/mL milk) in 2014 (Table 2). This result may reflect the improvement of management practices by producers, since SCC is directly associated with the health status of the cow’s mammary gland, which is the main indicator of mastitis occurrence (MACHADO et al., 2000MACHADO, P. F. et al. Somatic cells in milk of Brazilian herds. Scientia Agricola, v.57, n.2, p.359-361, 2000. Available from: <Available from: http://www.scielo.br/scielo.php?script=sci_arttext&pid=S0103-90162000000200026&lng=pt&tlng=pt >. Accessed: Sep. 12, 2019. doi: 10.1590/S0103-90162000000200026.
http://www.scielo.br/scielo.php?script=s...
). According to REZENDE et al. (2012REZENDE, E. S. J. et al. Hygienic-sanitary quality of raw milk in three regions of the state of Minas Gerais - Brazil. Revista do Instituto de Laticínios “Cândido Tostes”, v.67, n.387, p.64-69, 2012. Available from: <Available from: https://www.revistadoilct.com.br/rilct/article/view/252 >. Accessed: Jan. 12, 2019. doi: 10.5935/2238-6416.20120051.
), SCC values greater than 2.0 × 105 cells/mL indicate the occurrence of mastitis. However, it should be remembered that this year was marked by the implementation of a SCC standard for chilled raw milk, which is stricter than that in force at the time (BRASIL, 2011bBRASIL. Instrução Normativa n.62, de 29 de dezembro de 2011b. Aprova Regulamentos Técnicos de Produção, Identidade e Qualidade do Leite tipo A, Leite Cru Refrigerado, Leite Pasteurizado e o Regulamento Técnico da Coleta de Leite Cru Refrigerado e seu Transporte a Granel. Diário Oficial da República Federativa do Brasil, Brasília, 30 de dezembro de 2011. Seção 1, p.6-11.). Therefore, the result may also be associated with dairy farmers’ interest in meeting the new established standard.

For the TBC indicator, a reduction of 0.12 log10 and 0.23 log10 (cfu/mL milk) was observed in 2013 and 2018, respectively, and an increase of 0.10 log10 (cfu/mL milk), in 2017 (Table 2). The improvement observed for the TBC indicator is indicative of improved management practices (VELTHUIS and VAN ASSELDONK, 2011VELTHUIS, A. G. J.; VAN ASSELDONK, M. A. P. M. Process audits versus product quality monitoring of bulk milk. Journal of Dairy Science , v.94, p.235-249. 2011. Available from: <Available from: https://www.sciencedirect.com/science/article/pii/S0022030210006910 >. Accessed: Jan. 05, 2019. doi: 10.3168/jds.2010-3528.
https://www.sciencedirect.com/science/ar...
), milking hygienic-sanitary conditions (PANTOJA et al., 2009PANTOJA, J. C. et al. Associations among milk quality indicators in raw bulk milk. Journal of Dairy Science , v.92, p.4978-4987, 2009. Available from: <Available from: https://www.sciencedirect.com/science/article/pii/S002203020970829X >. Accessed: Sep. 13, 2019. doi: 10.3168/jds.2009-2329.
https://www.sciencedirect.com/science/ar...

Regarding the impacts of legislation on milk quality indicators, it should be noted that the establishment of stricter standards should be carried out in conjunction with other actions in order to inform, guide, and provide conditions for producers to achieve them. For example, there is an increase in the provision of technical assistance to producers and transporters of milk by the government, industries, or producer associations and an intensification of supervision and an increase in credit supply by the government. Increased milk production in Europe in the 1970s, for example, was stimulated by a subsidy policy (BERTOLDI, 2012BERTOLDI, S. A Qualificação do setor leiteiro no RS: impactos da implementação da Instrução Normativa 51/2002. Porto Alegre, RS: UFRGS, 2012. 37p. Trabalho de Conclusão de Curso (Bacharel em Ciências Econômicas) - Universidade Federal do Rio Grande do Sul, 2012. Available from: <Available from: https://lume.ufrgs.br/handle/10183/56108 >. Accessed: Sep. 12, 2019.
https://lume.ufrgs.br/handle/10183/56108...
); while the world’s largest milk producer, the United States, has its production subsidized by a minimum price policy (CARVALHO, 2008CARVALHO, V. R. F. Restructuring of the global dairy production: an analysis of Brazilian integration. 2008. 171p. Tese (Doutorado em Ciências Econômicas) - Universidade Estadual de Campinas. Campinas 2008. Available from: <Available from: http://repositorio.unicamp.br/bitstream/REPOSIP/285788/1/Carvalho_VeraReginaFerreira_D.pdf >. Accessed: Sep. 12, 2019.
http://repositorio.unicamp.br/bitstream/...
).

Volume influenced the TBC indicator over the period analyzed. The results indicate that a 1,000 L milk size increase in a producer’s monthly output is associated with a 0.01 log10 (cfu/mL milk) increase in TBC (Table 2). This result, in turn, indicates that the larger the volume of milk supplied by the producer, i.e. its size, the higher the bacterial counts of milk supplied by the producer. In a study with quality data of chilled raw milk from a cooperative located in southern Brazil, between August 2005 and April 2008, BOTARO et al. (2013BOTARO, B. G. et al. Quality based payment program and milk quality. Scientia Agricola, v.70, n.1, p.21-26, 2013. Available from: <Available from: http://www.scielo.br/scielo.php?script=sci_arttext&pid=S0103-90162013000100004 >. Accessed: Sep. 08, 2019. doi: 10.1590/S0103-90162013000100004.
http://www.scielo.br/scielo.php?script=s...
) identified a significant effect (P<0.05) of volume on the variables SCC, TBC, and protein content. The authors also identified that higher TBC values ​​were associated with producers who provided the largest volume of milk. However, it is noteworthy that, in theory, producers who provide larger volumes of milk tend to be more automated than those who provide smaller volumes (RODRIGUES et al., 2005RODRIGUES, A. C. O. et al. Management of Wisconsin dairy herds enrolled in milk quality teams. Journal of Dairy Science , v.88, p.2660-2671, 2005. Available from: <Available from: https://www.journalofdairyscience.org/article/S0022-0302(05)72943-X/abstract >. Accessed: Dec. 21, 2018. doi: 10.3168/jds.S0022-0302(05)72943-X.
https://www.journalofdairyscience.org/ar...
; O’CONNELL et al., 2015O’CONNELL, A. et al. Seasonal trends in milk quality in Ireland between 2007 and 2011. Journal of Dairy Science , v.98, p.3778-3790, 2015. Available from: <Available from: https://www.sciencedirect.com/science/article/pii/S002203021500212X >. Accessed: Jan. 18, 2018. doi: 10.3168/jds.2014-9001.
https://www.sciencedirect.com/science/ar...
). Therefore, larger producers could be expected to have lower values ​​for TBC, as the less automated milking process tends to make milk susceptible to a greater amount of contaminants.

All quality indicators of chilled raw milk were influenced by Seasonality (Table 2). The maximum values of the indicators of milk composition were observed in the months of April, the end of the rainy season, while the minimum values of these indicators were observed in the months of November (Figures 1a, 1b, 1c), just after the dry season in the region. For the SCC and TBC indicators, the maximum values were observed in January during the summer, the period with the highest temperatures, while the minimum values were observed in July (Figures 1d, 1e), during the winter, the period with the lowest temperatures.

The rainy season, concentrated in the summer, makes it difficult to maintain hygienic conditions during milking (BUENO et al., 2008BUENO, V. F. F. et al. Total bacterial count: relationship to milk composition and period of the year in Goiás State, Brazil. Revista Brasileira de Ciência Veterinária, v.5, n.1, p.40-44, 2008. Available from: <http://doi.editoracubo.com.br/10.4322/rbcv.2014.194>. Accessed: Sep. 12, 2019.doi: 10.4322/rbcv.2014.194.
https://doi.org/10.4322/rbcv.2014.194....
; CASTRO et al., 2016CASTRO, R. D. et al. Lactic acid microbiota identification in water, raw milk, endogenous starter culture, and fresh Minas artisanal cheese from the Campo das Vertentes region of Brazil during the dry and rainy seasons. Journal of Dairy Science, v.99, p.6086-6096, 2016. Available from: <Available from: https://www.sciencedirect.com/science/article/pii/S0022030216303502 >. Accessed: Sep. 12, 2019. doi: 10.3168/jds.2015-10579.
https://www.sciencedirect.com/science/ar...
) and, consequently, contributes to the milk being susceptible to a larger number of microorganisms. In addition, according to FLORES-MIYAMOTO et al. (2014FLORES-MIYAMOTO, A. et al. Do farm audits improve milk quality? Journal of Dairy Science, v.97, n.1, p.1-9, 2014. Available from: <Available from: https://www.sciencedirect.com/science/article/pii/S0022030213007303 >. Accessed: Sep. 12, 2019. doi: 10.3168/jds.2012-6228.
https://www.sciencedirect.com/science/ar...
), the average bacterial count of milk tends to increase in summer, since temperatures are higher, hindering the process of preservation of the product. Other common sources of contamination of refrigerated raw milk are the presence of lesions on the animal’s skin and teat (CASTRO et al., 2016CASTRO, R. D. et al. Lactic acid microbiota identification in water, raw milk, endogenous starter culture, and fresh Minas artisanal cheese from the Campo das Vertentes region of Brazil during the dry and rainy seasons. Journal of Dairy Science, v.99, p.6086-6096, 2016. Available from: <Available from: https://www.sciencedirect.com/science/article/pii/S0022030216303502 >. Accessed: Sep. 12, 2019. doi: 10.3168/jds.2015-10579.
https://www.sciencedirect.com/science/ar...
), failures in teat hygiene procedures before milking (ELMOSLEMANY et al., 2010ELMOSLEMANY, A. M. et al. The association between bulk tank milk analysis for milk quality and on-farm management practices. Preventive Veterinary Medicine, v.95, p.32-40, 2010. Available from: <Available from: https://www.sciencedirect.com/science/article/pii/S0167587710000930?via%3Dihub >. Accessed: Sep. 12, 2019. doi: 10.1016/j.prevetmed.2010.03.007.
https://www.sciencedirect.com/science/ar...
; OLIVEIRA et al., 2011OLIVEIRA, C. J. B. et al. Risk factors associated with selected indicators of milk quality in semiarid northeastern Brazil. Journal of Dairy Science , v.94, p.3166-3175, 2011. Available from: <Available from: https://www.sciencedirect.com/science/article/pii/S0022030211003043 >. Accessed: Sep. 13, 2019. doi: 10.3168/jds.2010-3471.
https://www.sciencedirect.com/science/ar...
), conditions of the environment in which the animals were raised, the milking equipment, employees, tanks used for transportation (CASTRO et al., 2016CASTRO, R. D. et al. Lactic acid microbiota identification in water, raw milk, endogenous starter culture, and fresh Minas artisanal cheese from the Campo das Vertentes region of Brazil during the dry and rainy seasons. Journal of Dairy Science, v.99, p.6086-6096, 2016. Available from: <Available from: https://www.sciencedirect.com/science/article/pii/S0022030216303502 >. Accessed: Sep. 12, 2019. doi: 10.3168/jds.2015-10579.
https://www.sciencedirect.com/science/ar...
) and milk storage.

With regard to SCC, it is known that periods of heat stress, most common in summer, tend to expose animals to a higher number of pathogens. This greater exposure to periods of thermal stress, in turn, increases susceptibility to infections and, consequently, the occurrence of mastitis (SMITH et al., 1985SMITH, K. L. et al. Environmental mastitis: cause, prevalence, prevention. Journal of Dairy Science , v.68, p.1531-1553, 1985. Available from: <Available from: https://www.journalofdairyscience.org/article/S0022-0302(85)80993-0/abstract >. Accessed: Jan. 05, 2019. doi: 10.3168/jds.S0022-0302(85)80993-0.
https://www.journalofdairyscience.org/ar...
), which increases the milk SCC.

Regarding fat, protein, and DDE contents, the results (Table 2) are associated with the influence of the production system (MILANI et al., 2016MILANI, M. P. et al. Milk quality in different production systems, year and weather seasons. Revista Brasileira de Ciência Veterinária, v.23, n.3-4, p.206-211, 2016. Available from: <Available from: http://doi.editoracubo.com.br/10.4322/rbcv.2016.058 . Accessed: Jan. 19, 2018. doi: 10.4322/rbcv.2016.058.
http://doi.editoracubo.com.br/10.4322/rb...
) and diet (JENKINS and McGUIRE, 2006JENKINS, T. C.; McGUIRE, M. A. Major Advances in Nutrition: Impact on Milk Composition. Journal of Dairy Science , v.89, p.1302-1310, 2006. Available from: <Available from: https://www.sciencedirect.com/science/article/pii/S0022030206721981 >. Accessed: Sep. 13, 2019. doi: 10.3168/jds.S0022-0302(06)72198-1.
https://www.sciencedirect.com/science/ar...
) on the milk composition indicators. The dairy herd of the region under study is raised on pasture, which makes the quality of its food very dependent on the available pastures and, consequently, on the weather seasons.

In turn, the results regarding the seasonality effect on quality indicators of chilled raw milk (Table 2) corroborate those found in several other studies (NIGHTINGALE et al., 2008NIGHTINGALE, C. et al. Influence of variable milk quality premiums on observed milk quality. Journal of Dairy Science , v.91, p.1236-1244, 2008. Available from: <Available from: https://www.sciencedirect.com/science/article/pii/S0022030208713821 >. Accessed: Jan, 18, 2018. doi: 10.3168/jds.2007-0609.
https://www.sciencedirect.com/science/ar...
; BOTARO et al., 2013BOTARO, B. G. et al. Quality based payment program and milk quality. Scientia Agricola, v.70, n.1, p.21-26, 2013. Available from: <Available from: http://www.scielo.br/scielo.php?script=sci_arttext&pid=S0103-90162013000100004 >. Accessed: Sep. 08, 2019. doi: 10.1590/S0103-90162013000100004.
http://www.scielo.br/scielo.php?script=s...
; FLORES-MIYAMOTO et al., 2014FLORES-MIYAMOTO, A. et al. Do farm audits improve milk quality? Journal of Dairy Science, v.97, n.1, p.1-9, 2014. Available from: <Available from: https://www.sciencedirect.com/science/article/pii/S0022030213007303 >. Accessed: Sep. 12, 2019. doi: 10.3168/jds.2012-6228.
https://www.sciencedirect.com/science/ar...
; MILANI et al., 2016MILANI, M. P. et al. Milk quality in different production systems, year and weather seasons. Revista Brasileira de Ciência Veterinária, v.23, n.3-4, p.206-211, 2016. Available from: <Available from: http://doi.editoracubo.com.br/10.4322/rbcv.2016.058 . Accessed: Jan. 19, 2018. doi: 10.4322/rbcv.2016.058.
http://doi.editoracubo.com.br/10.4322/rb...
; RIBAS et al., 2016RIBAS, N. P. et al. Total bacterial count in milk samples from bulk tank Paraná State. Archives of Veterinary Science, v.21, n.1, p.32-43, 2016. Available from: <Available from: https://revistas.ufpr.br/veterinary/article/view/41581 >. Accessed: Dec. 22, 2018. doi: 10.5380/avs.v21i1.41581.
https://revistas.ufpr.br/veterinary/arti...
). MILANI et al. (2016)MILANI, M. P. et al. Milk quality in different production systems, year and weather seasons. Revista Brasileira de Ciência Veterinária, v.23, n.3-4, p.206-211, 2016. Available from: <Available from: http://doi.editoracubo.com.br/10.4322/rbcv.2016.058 . Accessed: Jan. 19, 2018. doi: 10.4322/rbcv.2016.058.
http://doi.editoracubo.com.br/10.4322/rb...
identified a significant effect of seasonality on fat and protein content, while NIGHTINGALE et al. (2008)NIGHTINGALE, C. et al. Influence of variable milk quality premiums on observed milk quality. Journal of Dairy Science , v.91, p.1236-1244, 2008. Available from: <Available from: https://www.sciencedirect.com/science/article/pii/S0022030208713821 >. Accessed: Jan, 18, 2018. doi: 10.3168/jds.2007-0609.
https://www.sciencedirect.com/science/ar...
and MILANI et al. (2016) MILANI, M. P. et al. Milk quality in different production systems, year and weather seasons. Revista Brasileira de Ciência Veterinária, v.23, n.3-4, p.206-211, 2016. Available from: <Available from: http://doi.editoracubo.com.br/10.4322/rbcv.2016.058 . Accessed: Jan. 19, 2018. doi: 10.4322/rbcv.2016.058.
http://doi.editoracubo.com.br/10.4322/rb...
identified a significant seasonality effect under the SCC indicator, and BOTARO et al. (2013)BOTARO, B. G. et al. Quality based payment program and milk quality. Scientia Agricola, v.70, n.1, p.21-26, 2013. Available from: <Available from: http://www.scielo.br/scielo.php?script=sci_arttext&pid=S0103-90162013000100004 >. Accessed: Sep. 08, 2019. doi: 10.1590/S0103-90162013000100004.
http://www.scielo.br/scielo.php?script=s...
, FLORES-MIYAMOTO et al. (2014)FLORES-MIYAMOTO, A. et al. Do farm audits improve milk quality? Journal of Dairy Science, v.97, n.1, p.1-9, 2014. Available from: <Available from: https://www.sciencedirect.com/science/article/pii/S0022030213007303 >. Accessed: Sep. 12, 2019. doi: 10.3168/jds.2012-6228.
https://www.sciencedirect.com/science/ar...
, MILANI et al. (2016), MILANI, M. P. et al. Milk quality in different production systems, year and weather seasons. Revista Brasileira de Ciência Veterinária, v.23, n.3-4, p.206-211, 2016. Available from: <Available from: http://doi.editoracubo.com.br/10.4322/rbcv.2016.058 . Accessed: Jan. 19, 2018. doi: 10.4322/rbcv.2016.058.
http://doi.editoracubo.com.br/10.4322/rb...
and RIBAS et al. (2016) RIBAS, N. P. et al. Total bacterial count in milk samples from bulk tank Paraná State. Archives of Veterinary Science, v.21, n.1, p.32-43, 2016. Available from: <Available from: https://revistas.ufpr.br/veterinary/article/view/41581 >. Accessed: Dec. 22, 2018. doi: 10.5380/avs.v21i1.41581.
https://revistas.ufpr.br/veterinary/arti...
identified a significant seasonality effect on the TBC indicator of chilled raw milk.

Finally, it can be inferred that the quality of chilled raw milk from the evaluated producers deteriorated between 2012 and 2018. Fat, protein and defatted dry extract levels decreased continuously between 2014 and 2018. Somatic cell count and total bacterial count indicators showed only point reductions over the analyzed period. Thus, regarding the group of producers studied, it is not possible to state that there was an improvement in milk quality. Therefore, further longitudinal studies with milk producers from other regions of the country are needed to (i) map the regions that contain the producers with the greatest difficulties in meeting the evolution of quality standards of chilled raw milk and (ii) enable the adoption of specific milk quality improvement strategies in each of these regions.

In a study of 517 dairy herds from different regions of Brazil, between January 2011 and May 2015, BUSANELLO et al. (2017BUSANELLO, M. et al. Estimation of prevalence and incidence of subclinical mastitis in a large population of Brazilian dairy herds. Journal of Dairy Science, v.100, p.6545-6553, 2017. Available from: < Available from: https://www.sciencedirect.com/science/article/pii/S0022030217305829?via%3Dihub >. Accessed: Sep. 12, 2019. doi: 10.3168/jds.2016-12042.
https://www.sciencedirect.com/science/ar...
) found a high percentage of cases of subclinical mastitis (mean prevalence of 46.4%) in all regions studied. The authors also identified an increasing tendency of these cases over the years.

# CONCLUSION:

Despite the actions implemented by PNMQL, the quality of chilled raw milk from the evaluated producers deteriorated over the period under review. Fat, protein and defatted dry extract contents decreased continuously between 2014 and 2018. The indicators of somatic cell count and total bacterial count showed only occasional reductions over the analyzed period.

Therefore, effective actions by dairy farmers and producers are needed to change this reality. Just the implementation of regulations or stricter indices for chilled raw milk quality parameters are not enough. These actions need to involve more effective technical assistance, establishment of payment for quality with bonuses, and even authomation of producers.

# ACKNOWLEDGEMENTS

To Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES) - Financing Code 001, for supporting the L.P.L. with postgraduate fellowship. To the industry that supplied the data and the Universidade Federal de Viçosa (UFV).

# REFERENCES

• CR-2019-0749.R1

# Publication Dates

• Publication in this collection
30 Mar 2020
• Date of issue
2020