Livestock Research for Rural Development 28 (8) 2016 | Guide for preparation of papers | LRRD Newsletter | Citation of this paper |
The objective of this study was to evaluate the composition and sanitary quality of the milk in cooling tanks for commercial herds of Guzerat and F1 Guzolando cross in climatic seasons in the state of Ceará. The study was conducted in two private farms located in Madalena-CE. During 15 months, milk samples were collected directly from milk storing tanks, after mechanical homogenization agitation of the stored milk. The milk samples were refrigerated and sent to the Clínica do Leite, Universidade de São Paulo – ESALQ, Piracicaba, SP, Brazil, to analyse for fat, protein, lactose, somatic cell count (SCC) and total bacterial count (TBC). The data were submitted to variance analysis and "F" Fischer testing and the differences between averages were compared using the Tukey test.
The following average composition was found for Guzerat and Guzolando milk, respectively: 4.57% and 4.05% fat; 3.71% and 3.39% protein; 4.67% and 4.61% lactose. The milk fat content did not vary according to the year or seasons. A higher percentage of protein was found in the 2012 year of milk without changes according to the seasons. The lactose content was higher during 2012 and during the winter and spring. With regard to climatic season’s stations, SCC and TBC showed the lowest values during the spring, with variations in other seasons. Guzerat cows produce milk with a higher fat, protein and lactose. In the Guzolando cattle studied, the composition and sanitary quality of the milk were influenced by the years and by the seasons.
Key words: Bos indicus, Bos taurus, dairy cattle, livestock, semi-arid
The production of cow's milk has high social and economic importance in Brazil for substantially constituting the Gross Domestic Product (GDP) of agribusiness and is an activity which is conducted in almost all Brazilian municipalities. However, the low productivity of most Brazilian herds combined with the exponential increase in human population and hence the increased consumption of milk and dairy products to meet the demand for animal protein for human consumption has increasingly pushed farmers to seek more production systems adapted to Brazilian conditions and greater economic efficiency of the activity (Mello Lima 2011).
The improvement of production conditions in tropical regions depends on several factors related to the animal, the environment and production technology. In the context of the animal and based on the sustainability of production systems, an alternative has emerged by cross-breeding Zebu animals ( Bos indicus) with specialized dairy breeds (Bos taurus), mainly Holstein, exploiting the advantages of the heterosis effect. Several studies have confirmed the advantages of using Bos taurus x Bos indicus in tropical conditions (Ordaz 2009), and especially the result of such breeding is raising of animals with high productivity and high adaptability to the conditions in Brazil, especially in the Northeast (Mello Lima 2011).
Our objective was to evaluate the composition (fat, protein, lactose) and sanitary conditions (SCC, TBC) of milk in milk cooling tanks of commercial herds of Guzerat and F1 Guzolando cross (½ Guzerat : ½ Holstein) breeds in the climatic seasons in the state of Ceará, Brazil.
The study was conducted at the Fazenda Potrinha (Farm I) and São Carlos (Farm II) farms, both private farms located in the municipality of Madalena-CE, latitude 04° 51'25'' South and longitude 39° 34'37'' West, with an average altitude of 299 meters above sea level, in the Sertões Cearenses mesoregion of Ceará State belonging to the micro region of Quixeramobim backcountry. The climate is classified as AS (tropical with dry season), by the Köppen-Geiger classification. The average rainfall ranges from 600 mm to 1000 mm per year.
Due to the low rainfall during the experimental period (July 2012 to September 2013), animals were kept in confinement, receiving sorghum or corn silage as forage and a concentrate: soybean meal (18.2%), ground corn grain (75.80 %), limestone (2.80%), common salt (1.00%), mineral salt (1.70%) and livestock urea (0.50%).
The cows were milked twice a day (5am and 2pm) by milking machines at both properties. The milking parlor was equipped with a single file double-six milking machine, with a milk line channeled to the expansion tank. In 2012, the herd from Farm I was exclusively composed of 813 pure Guzerat cattle, registered by the Brazilian Association of Zebu Breeders (Associação Brasileira dos Criadores de Zebu - ABCZ), and evaluated in the Genetic Improvement Program of Zebu breeds (Programa de Melhoramento Genético das Raças Zebuínas – PMGZ, Milk control at farm and Performance Growth Control (PGC). For Farm II, the herd in 2012 and 2013 consisted of Guzolando (½ Guzerá x ½ Holstein), totaling 2,180 milking cows, which were registered at ABCZ in the category of Certified Genealogy Control (CGC).
It is important to note that during the experimental period, the average milk production was 8.36 kg/day for Guzerat cows and 15.02 kg/day for Guzolando cows.
For Farm I, only samples and analyzes carried out in 2012 were considered valid. In 2013, the Guzolando herd was introduced at the farm, leading to a mixture of milk coming from the two genetic groups in the milk tank and forcing us to ignore samples.
Data relating to 15 months of samples from July 2012 to September 2013 were analyzed. The milk sampling procedure was performed directly from the milk tank after homogenization by mechanical agitation. 214 samples were collected: 107 in 40 mL vials containing the preservative Bronopol®; and 107 in 40 mL vials with Azidiol®. Samples were identified, refrigerated between 2°C and 6°C, and sent to the Clínica do Leite, Universidade de São Paulo – ESALQ, Piracicaba, SP, Brazil, and milk components were analyzed for fat, protein and lactose, as well as for somatic cell count (SCC) and total bacterial count (TBC).
Milk composition analyzes were performed electronically by infrared absorption and somatic cell count (SCC) was determined using electronic counting by Fourier transform infrared spectroscopy (FTIR) with MILKO SCAN™ (Foss Eletric A/S, HillerØd, Denmark). TBC was analyzed by flow cytometry method, using Bactocount® (BENTLEY Instruments Inc., Chasca MN, USA). To ensure reliability of the records, null values for bacterial count and somatic cell count were excluded, as well as fat values below 2.0% and above 5.0%. For the characteristic of SCC, statistic W Shapiro-Wilk and statistic D Kolmogorov-Smirnov tests were used to check whether the observed values followed a normal distribution. As the data did not follow the normal distribution, linear score was then calculated (ECS) from the SCC data using the formula proposed by Dabdoub & Shook (1984), in which the ECS = [Log2 (SCC/100,000) + 3].
Climatic seasons and year were defined as experimental variables. Summer (Season 1), comprised of the months of January, February and March, being the start of the rainy season; Autumn (Season 2), comprising April, May and June, and the end of the rainy season; Winter (Season 3), comprised of July, August and September; and Spring (Season 4), comprising October, November and December. Data were tabulated and later submitted to variance analysis and "F" Fischer test. The differences between averages were compared using the Tukey test, at 5% probability. Statistical procedures were performed using the SAS program (Statistical Analysis System, SAS Institute, version 9.3).
The general mathematical model used was:
Where:
Yjk = milk composition characteristics and the hygienic and sanitary quality of the milk of the j-th year (j=1,2), k-th season of the year (k=1, 2, 3, 4);
j-th year (j=1,2), k-th season of the year (k=1, 2, 3, 4);
The averages for milk composition from the two experimental groups studied are arranged in Tables 1 and 2.
Table 1. Composition and sanitary quality of Guzerat cow’s milk |
||||
Variable |
Average ± SD |
CV (%) |
Minimum |
Maximum |
Composition |
||||
Fat (%) |
4.57±0.33 |
7.24 |
3.70 |
5.00 |
Protein (%) |
3.71±0.14 |
3.76 |
3.49 |
3.92 |
Lactose (%) |
4.67±0.07 |
1.44 |
4.53 |
4.82 |
Hygienic and sanitary quality of milk |
||||
SCC (10³ cls/mL) |
392±85.2 |
21.7 |
188 |
590 |
SCS (log cls/mL) |
4.34±0.241 |
5.56 |
3.63 |
4.77 |
TBC (10³ CFU/mL) |
51.6±32.4 |
62.7 |
13.0 |
43.0 |
SD = standard deviation; CV = coefficient of variation; SCS = somatic cell score. |
Data concerning fat, protein and lactose levels are similar to those described by Ribeiro et al (2009) for Guzerat, showing the influence of the breed in milk composition, especially in fat and protein and the relationship between these components.
Ribeiro Neto et al (2012) also studied crossbred cattle in the Northeast, and found lower average values for fat (3.66 ± 0.53), protein (3.16 ± 0.22) and lactose (4.41 ± 0.18) components, and higher average values for SCC (564.95 ± 653.56) and TBC (1,190.68 ± 1,384.61), compared to the evaluated herds in this study.
Minimum and maximum SCC values for Guzerat are consistent with the findings of Langoni et al (2011) who worked with raw stored milk and found values between 200,000 and 586,000 cells/ml.
Table 2 presents values that describe the quality of milk from Guzolando cows (½ Guzerat : ½ Holstein). Data are similar to those reported by Signoretti et al. (2013), when evaluating the composition of milk from crossbred cows kept in irrigated pasture of grass and supplemented with concentrate, wherein the following levels were present: fat (4.46% ± 0.06), protein (3.40% ± 0.03), lactose (4.41% ± 0.02) and SCC (360 ± 34).
Table 2. Composition and sanitary quality of Guzolando cow milk |
||||
Variable |
Average ± SD |
CV (%) |
Minimum |
Maximum |
Composition |
||||
Fat (%) |
4.05±0.26 |
6.52 |
3.44 |
5.00 |
Protein (%) |
3.39±0.15 |
4.33 |
3.14 |
3.70 |
Lactose (%) |
4.61±0.07 |
1.58 |
4.47 |
4.75 |
Hygienic and sanitary quality of milk |
||||
SCC (10³ cls/mL) |
884±412 |
46.6 |
357 |
2009 |
SCS (log cls/mL) |
5.07±0.477 |
9.41 |
4.27 |
6.00 |
TBC (10³ CFU/mL) |
139±229 |
165 |
7.00 |
1373 |
SD = standard deviation; CV = coefficient of variation; SCS = somatic cell score. |
Fat content levels obtained in this study are higher than those reported in the literature for animals specialized in milk production. Indeed, Zebu cows have a higher fat content in milk compared to Holstein cows (Fonseca and Santos, 2000). Mello Lima (2011) found 3.00% of milk fat in the Holstein cows and 3.98% in milk from Jersey cows. It is possible that the fat content found can be explained by the dilution effect of the milk components in the volume of milk produced. This effect was highlighted by Venturini et al (2007), who pointed out that milk fat content in general is inversely proportional to the milk yield. Reis et al (2013) also found fat content in milk for crossbred Bos indicus x Bos taurus to be between 3.43% and 3.63%.
Table 3 shows the tank milk composition for the Guzolando cows according to the collection years and the seasons. Milk fat content did not vary according to the year or seasons. Higher percentages of protein were found in milk in 2012 without changes according to the season. Lactose content was higher during 2012 and during winter and spring. Regarding SCC and TBC levels, it was observed in 2012 that they were well below those achieved during 2013. Concerning the seasons, SCC and TBC showed the lowest values during the spring, with variations in other seasons.
It is known that the milk fat content varies according to genetic and environmental factors (Ribeiro et al 2009), highlighting nutrition management and stressing that the quality of forage and forage to concentrate ratio interfere in milk fat percentage (Venturini et al 2007; Rosa et al 2012; Rangel et al 2014; Taffarel et al 2015).
Therefore, as there was no change in the genetic composition of the herd, and since the animals were kept in confinement during the experiment consuming silage and concentrate, it is possible to understand the lack of fat effect on farm II’s tank milk.
Table 3. Guzolando cow milk composition and health-hygiene parameters, according to the year and the seasons |
||||||||
Variable |
Year |
Season |
SEM |
p |
||||
2012 |
2013 |
Summer |
Autumn |
Winter |
Spring |
|||
Fat (%) |
4.10a |
4.03a |
4.05a |
4.07a |
3.96a |
4.17a |
0.030 |
0.21 |
Protein (%) |
3.44a |
3.30b |
3.30a |
3.37a |
3.39a |
3.36a |
0.014 |
<0.001 |
Lactose (%) |
4.65a |
4.56b |
4.58b |
4.56b |
4.59ab |
4.65a |
0.007 |
<0.001 |
SCC¹ |
683b |
1277a |
1083b |
1349a |
1010b |
657c |
29.87 |
<0.001 |
TBC² |
25.5b |
178a |
169a |
180a |
83.5ab |
28.8b |
16.60 |
0.001 |
abc
Averages followed by the same letter in the same column do not differ by Tukey test at 5% probability for type I error. |
Regarding milk protein content, the lack of effect of seasons on this variable could also be explained by the nutrition management, which did not change throughout the sample collection period. According to Rosa et al (2012), the effects of diet balance on protein levels are less apparent than changes in fat content. Thus, considering there were no changes in fat content, the variation in protein content was not expected. It is possible that this small variation in protein levels observed in accordance with the years elapsed secondary factors such as milk production and calving order (Galvão Junior et al 2010). High production values influence lower levels of fat and protein in milk (Peres et al 2001; Kgole et al 2012).
Lactose levels found in this study (Guzerat and Guzolando) are consistent with the values reported in the literature for Zebu cows (Ribeiro et al 2009; Galvão Junior et al 2010; Rangel et al 2014.) Few studies report variations in milk lactose level because there is not usually much variation in this component according to genetic and environmental factors (Signoretti et al 2013). However, in this study it is clear that variations in milk lactose level followed variations in SCC. This occurrence could be explained by the idea advocated by Silva et al (2014), who mentioned that lactose is the primary substrate for mesophilic bacteria that colonize the milk in the cooling milk tank, making it possible that the reduction and/or increased SCC has caused an increase and/or decrease in milk lactose content. The reduction could also be explained by the loss of lactose from the mammary gland to blood due to changes in the permeability of the dividing membrane, caused from the inflammatory condition of mastitis (Rangel et al 2009; Ruegg 2011; Akers and Nickerson 2011; Ballou 2012).
According to Normative Regulation 62 of the Ministry of Agriculture, Livestock and Supply (MAPA 2011), as of July 1, 2015, for the Northeast of Brazil, the maximum allowable values for SCC and CBT in refrigerated raw milk is 500,000 cells/mL and 300,000 CFU/ml, respectively. The values contained in Table 3 for such parameters are much higher than those allowed by the regulation and must be reviewed by the farmer to meet the NR 62 recommendation.
However, it is important to clarify that an infection problem of Staphylococcus aureus on farm II in dairy arrays was detected during the experimental period. Staphylococci are important causative agents of mastitis, highlighting S. aureus among them as a pathogen of mastitis classified as infectious (Sá et al 2004). This fact may justify the high values obtained for SCC, since according to Müller (2002), the extent of SCC increase is directly related to the surface of the mammary tissue affected by the inflammatory reaction. For the values of TBC, an infection in the herd also justifies the high values, because according to Beloti et al (2011), practices that prevent mastitis are largely the same that prevent contamination of the milk, so that the SCC and TBC are positively correlated.
It is worth mentioning that the results for the hygienic and sanitary parameters of milk reinforce the idea that management used in obtaining the milk can influence quality (Ruegg 2011; Silva et al 2013; Cicconi-Hogan et al 2013), particularly at farms that use machine milking. Vallin et al. (2009) found that the average TBC found at farms with machine milking was around three times higher than the average found at farms with hand milking. This situation becomes even more serious when there are highly contagious infectious agents such as S. aureus, , making it necessary to devote more attention in these cases to the process of cleaning and sanitizing the milking machines and all instruments involved in the process of collection and preservation of milk.
Akers R M and Nickerson S C 2011 Mastitis and its impact on structure and function in the ruminant mammary gland. Journal of Mammary Gland Biology and Neoplasia, Volume 16, 275-289. Retrieved January 30, 2016, from http://download.springer.com/static/pdf/211/art% 253A10.1007%252Fs10911-011-9231-3.pdf?originUrl=http%3A%2F%2Flink.sprin
Beloti V, Ribeiro Júnior J C and Tamanini R 2011 Qualidade microbiológica e físico-quimica do leite cru refrigerado produzido no município de Saporema/PR. Revista Científica Eletrônica de Medicina Veterinária, Volume 6, 2-5 Retrieved January 30, 2016, from http://faef.revista.inf.br/imagens_arquivos/arquivos_destaque/gvRfHOQjI5PmOHd_2013-6-25-16-55-49.pdf
Ballou M A 2012 Growth and development symposium: inflammation: Role in the etiology and pathophysiology of clinical mastitis in dairy cows. Journal of Animal Science, Volume 90, 1466-1478. Retrieved March 1, 2016, from https://dl.sciencesocieties.org/publications/jas/abstracts/90/5/1466?access=0&view=article
Cicconi-Hogan K M, Gamroth M, Richert R M, Ruegg P L, Stiglbauer K E and Schukken Y H 2013. Associations of risk factors with somatic cell count in bulk tank milk on organic and conventional dairy farms in the United States. Journal of Dairy Science, Volume 96, 3689-3702. Retrieved February 15, 2016, from http://jas.fass.org/
Dabdoub S M and Shook G E 1984. Phenotypic relations among milk yield, somatic cell count and clinical mastitis. Journal of Dairy Science, Volume 67 (s1).
Fonseca, L F L and Santos M V 2000. Qualidade do leite e controle da mastite. São Paulo: Lemos Editorial, 175p.
Galvão Júnior J G B, Rangel A H N, Medeiros H R, Silva J B A, Aguiar E M and Madruga R C 2010. Efeito da produção diária e da ordem de parto na composição físico-química do leite de vacas de raças zebuínas. Acta Veterinária Brasilica, Volume 4, 25-30. Retrieved February 15, 2016, from http://periodicos.ufersa.edu.br/revistas/index.php/acta
Kgole M L, Visser C and Banga C B 2012. Environmental factors influencing milk urea nitrogen in South African Holstein cattle. South African Journal of Animal Science, Volume 42, 459-463. Retrieved January 31, 2016, from http://www.sasas.co.za/journals
Langoni H, Penachio D S and Citadella J C C 2011. Aspectos microbiológicos e de qualidade do leite bovino. Pesquisa Veterinária Brasileira, Volume 31, 1059-1065. Retrieved February 20, 2016, from http://seer.sct.embrapa.br/index.php/pab
Mello Lima J A 2011. Desempenho produtivo da vacas F1 holandês/zebu submetidas ao aumento no número de ordenhas no inicio da lactação e a diferentes manejo de amamentação. 2011. 67f. Thesis (Doctor of Science - Zootecnia) Escola de Veterinária, Universidade federal de Minas Gerais.
Ministério da Agricultura, Pecuária e Abastecimento – MAPA. Instrução Normativa n.062 de 29 de dezembro 2011. Aprova os Regulamentos técnicos de produção, identidade e qualidade do leite tipo A, do leite tipo B, do leite tipo C, do leite pasteurizado e do leite cru refrigerado e o regulamento técnico da coleta do leite cru refrigerado e seu transporte a granel. Brasília, 2011. Retrieved August 23, 2014, form http://www. agricultura.gov.br/Instrução Normativa 62
Müller E E 2002. Qualidade do leite, células somáticas e prevenção da mastite. In: Simpósio sobre Sustentabilidade da Pecuária Leiteira na Região Sul do Brasil, 2., 2002. Maringá. Anais... Maringá: II Sul-Leite, UEM/CCA/DZO – NUPEL, 2002, 206-217.
Ordaz R L 2009. Producción de leche de vacas com diferente porcentaje de genes Bos taurus em El trópico mexicano. Técnica Pecuária en México, Volume 47, 435-447. Retrieved December, 2015, from http://www.redalyc.org/revista.oa?id=613
Peres J R, Gonzáles F H D, Dürr J W and Fontaneli R S 2001. Uso do leite para monitorar a nutrição e o metabolismo de vacas leiteiras. Porto Alegre: Editora da UFRGS, pp.30-45.
Rangel A H N, Medeiros H R, Silva J B, Barreto M L D J and Lima Júnior D M 2009. Correlação entre a contagem de células somáticas (CCS) e o teor de gordura, proteína, lactose e extrato secodesengordurado do leite. Revista Verde de Agroecologia e Desenvolvimento Sustentável, Volume 4, 57-60. Retrieved January 12, 2016, from http://www.gvaa.com.br/revista/index.php/RVADS/issue/archive
Rangel A H N, Araújo T P M, Medeiros H R, Lima Júnior D M, Andrade K D, Madruga R C, Bezerril R F and Novaes L P 2014. Body condition score (BCS) and calving order on milk production and composition in zebu cattle. Acta Veterinaria Brasilica, Volume 8, 247-53. Retrieved January 10, 2016, from http://periodicos.ufersa.edu.br/revistas/index.php/acta
Reis A M, Rezende Costa M, Costa R G, Suguimoto H H, Souza C H B, Aragon-Alegro L C and Santan E H W 2013. Efeito do grupo racial e do número de lactações sobre a produtividade e a composição do leite bovino. Semina: Ciências Agrárias, Volume 33 (6Supl2), 3421-3436. Retrieved January 12, 2016, from http://www.uel.br/revistas/uel/index.php/semagrarias/issue/archive
Ribeiro A B, Tinoco A F F, Ferreira G, Guilhermino M M and Rangel A H N 2009. Produção e composição do leite de vacas gir e guzerá nas diferentes ordens de parto. Revista Caatinga, Volume 22, 46-51. Retrieved March 1, 2016, from http://periodicos.ufersa.edu.br/revistas/index.php/sistema/issue/archive
Ribeiro Neto A C, Barbosa S B P, Jatobá R B, Silva A M, Silva C X, Silva M J A and Santoro K R 2012. Qualidade do leite cru refrigerado sob inspeção federal na região Nordeste. Arquivo Brasileiro de Medicina Veterinária e Zootecnia, Volume 64, 1343-1351. Retrivied January 12, 2016, from http://www.scielo.br/scielo.php?script=sci_issuetoc&pid=0102093520150006&lng=en&nrm=iso
Rosa D C, Trenentin J M, Pessoa G A, Silva C A M and Rubin M I B 2012. Qualidade do leite em amostras individuais e de tanque de vacas leiteiras. Arquivos do Instituto de Biologia, Volume 79, 485-493. Retrieved December 12, 2015, from http://www.biologico.sp.gov.br/rev_arq_todos.php
Ruegg P L 2011. Managing mastitis and producing high quality milk. In: “Dairy Cattle Production Medicine” (eds. Risco, C. and Melendez, P.), Wiley-Blackwell Publishing Ltd, Ames, Iowa, USA.
Sá M E P, Cunha M L R S, Elias A O, Vitória C and Langoni H 2004. Importância do Staphylococcus aureus nas mastites subclínicas: pesquisa de enterotoxinas e toxina do choque tóxico, e a relação com a contagem de células somáticas. Brazilian Journal of Veterinary Research and Animal Science, Volume 41, 320-326. Retrieved January 12, 2016, from http://www.revistas.usp.br/bjvras/issue/archive
SAS INSTITUTE 2011. Sas/Iml 9.3 User's Guide. Sas Institute.
Silva M A P, Brasil R B, Cabral J F, Garcia J C and Oliveira A N 2013. Qualidade do leite cru refrigerado obtido através de ordenha manual e mecânica. Revista do Instituto de Laticínios Cândido Tostes, Volume 68, 05-11. Retrieved February 14, 2016, from http://www.revistadoilct.com.br/rilct/article/view/2
Signoretti R D, Resende F D, Pessim B, Souza F H M and Souza L A 2013. Produção e composição de leite de vacas mestiças suplementadas com diferentes níveis de concentrado mantidas em pastagem de capim Tanzânia irrigada. Boletim de Indústria Animal, Volume 70, 10-19. Retrieved February 3, 2016, from http://www.iz.sp.gov.br/bia/bias.php
Silva V N, Rangel A H N, Novaes L P N, Borba L H F, Bezerril R Fand Lima Júnior D M 2014. Correlação entre a contagem de células somáticas e composição química no leite cru resfriado em propriedades do rio grande do norte. Revista do Instituto Laticínio Cândido Tostes, Volume 69, 165-172. Retrieved February 14, 2016, from http://www.revistadoilct.com.br/rilct/issue/archive
Taffarel L E, Costa P B, Tsutsumi C Y, Klosowski E S, Portugal E F and Lins A C 2015. Variação da composição e qualidade do leite em função do volume de produção, período do ano e sistemas de ordenha e de resfriamento. Semina: Ciências Agrárias, Volume 36, 2287-2300. Retrieved January 12, 2016, from http://www.uel.br/revistas/uel/index.php/semagrarias/issue/archive
Vallin V M, Beloti V, Battaglini A P P, Tamanini R, Fagnani R, Angela H L and Silva L C C 2009. Melhoria da qualidade do leite a partir da implantação de boas práticas de higiene na ordenha em 19 municípios da região central do Paraná. Semina: Ciências Agrárias, Volume 30, 181-188. Retrieved February 2 12, 2016, from http://www.uel.br/revistas/uel/index.php/semagrarias/issue/archive
Venturini K S, Sarcinelli M F and Silva L S 2007. Características do leite. Niterói: UFES, 2007. Boletim Técnico PIES: 01007.
Received 8 April 2016; Accepted 7 July 2016; Published 1 August 2016