Agronomía Mesoamericana
Artículo científico
Volumen 32(3):939-948. Septiembre-diciembre, 2021
e-ISSN 2215-3608, doi:10.15517/am.v32i3.45471
https://revistas.ucr.ac.cr/index.php/agromeso/index
Pablo Montero-Prado2, Giancarlo Antonio Ruiz-Morales2, Anibal Fossatti-Carrillo2
1 Recepción: 21 de enero, 2021. Aceptación: 12 de abril, 2021. This work was part of a Project sponsored by the Secretaria Nacional de Ciencia, Tecnología de Innovación de Panamá (SENACYT); Code IAP17-P-047.
2 Universidad Tecnológica de Panamá, Centro de Innovación y Transferencia Tecnológica. Grupo de Investigación en Innovación, Desarrollo y Transferencia de Agrotecnología (IDTAT). 0819-07289, El Dorado, Panamá, República de Panamá. pablo.montero@utp.ac.pa (autor para correspondencia; https://orcid.org/0000-0002-1471-2741); giancarlo.ruiz@utp.ac.pa (https://orcid.org/0000-0002-0862-9968); anibal.fossatti@utp.ac.pa (https://orcid.org/0000-0003-4284-3500).
Introduction. Currently, there are no bibliographic or statistical records that allow establishing a correlation between the physical-chemical quality characteristics of raw cow’s milk produced in the Republic of Panama and the classification it receives. Objective. To characterize and quantify the physical-chemical parameters of the quality of raw cow’s milk and correlate these characteristics with the classification that milk receives in Panama. Materials and methods. A stratified sampling was carried out according to the classification of raw cow´s milk (grade A, B, and C), 48 of grade A, 71 of grade B, and 393 of grade C. Random samples were taken in Cocle, Veraguas, Herrera and Los Santos during the rainy season from June to September 2018. In the Chiriquí province, the sampling was conducted during the dry season, from January to March 2019. In addition, it was necessary to have at least 512 producers or farms, so that the analysis reached a confidence level of 95 %, a standard deviation of 5 %, and a limit of error of 5 %. Results. The values found allow us to consider that the milk produced in Panama has the physical-chemical characteristics according to international standards. The correlation analysis indicated that the protein had a high relationship with solids non-fat (SNF), density (D), and lactose (L) (0.93, 0.91, and 0.96, respectively); also L is strongly related to SNF and D (0.94 and 0, 91, respectively). Conclusions. The analysis of variance did not allow to establish significant differences between the three existing classifications in the country, which allow a distribution of milk according to the grades of these parameters. However, it was possible to establish the existing correlation between some of these parameters.
Keywords: unprocessed milk, dairy herd, physic-chemical parameters, milk production, grade.
Introducción. Actualmente, no existen registros bibliográficos o estadísticos que permitan establecer una correlación entre las características de calidad físico-químicas de la leche cruda de vaca producida en la República de Panamá y la clasificación que recibe. Objetivo. Caracterizar y cuantificar los parámetros físico-químicos de la calidad de la leche cruda de vaca y correlacionar estas características con la clasificación que recibe la leche en Panamá. Materiales y métodos. Se realizó un muestreo estratificado según la clasificación de la leche cruda de vaca (tipos A, B y C), 48 del grado A, 71 del grado B y 393 del grado C. Se tomaron muestras aleatorias en las provincias de Coclé, Veraguas, Herrera y Los Santos durante la época lluviosa, de junio a septiembre de 2018. En la provincia de Chiriquí, el muestreo se realizó durante la época seca, de enero a marzo de 2019. Además, fue necesario contar con 512 productores o fincas, para que el análisis alcanzara un nivel de confianza del 95 %, una desviación estándar del 5 % y un límite de error del 5 %. Resultados. Los valores encontrados permiten considerar que la leche producida en Panamá tiene las características físico-químicas acordes a normas internacionales. El análisis de correlación indicó que la proteína tuvo una alta relación con los sólidos no grasos (SNF), la densidad (D) y la lactosa (L) (0,93, 0,91 y 0,96, respectivamente); también L está fuertemente relacionada con SNF y D (0,94 y 0, 91, respectivamente). Conclusiones. El análisis de varianza no permitió establecer diferencias significativas entre las tres clasificaciones existentes en el país, que permiten una distribución de la leche de acuerdo a los grados de estos parámetros. Sin embargo, fue posible establecer la correlación existente entre algunos de estos parámetros.
Palabras claves: leche no procesada, hato lechero, parámetros físico-químicos, producción láctea.
Considering that milk is the only material produced by nature to function exclusively as a food source, since it constitutes a nutritional source complete, not surpassed by any other known to humans (Jay, 2009).
This food consumption may represent differences in the quality of diets of people without access to adequate nutrition, especially among children. In this sense, in recent years, world milk production exceeds 600 million tons per year according to the Food and Agriculture Organization Statistics Division (Food and Agriculture Organization, 2020). The main producers are the countries of the European Union, the United States and the Russian Federation, contributing around 45 %.
In the last 10 years, the production of fresh milk has suffered negative fluctuations, in 2015, the Republic of Panama reached a production of 125 million liters, reflecting a 12 % increase in production (Instituto Nacional de Estadística y Censo, 2020). With this increase, the capacity to satisfy local consumption is far from being achieved as approximately 50 % of the demand is generated.
In Panama, milk production occurs mainly on small farms, mostly carried as a form of family income. According to current regulations, the milk production receives a classification based on the quality of this, milk grade A, grade B and grade C or Industrial. This classification (A, B and C) is mainly related to the bacterial load that the milk has: a) grade A milk <200,000 bacteria ml-1, b) grade B milk with <500,000 bacteria ml-1, and c) grade C milk, milk that does not meet the parameters required by grade A and B milk (Ministerio de Gobierno y Justicia, 1970).
Actually, there are approximately 5,000 milk producers. Around the 8 % of them, generate grade A milk, while 12 and 80 % provide grade B and C, respectively. Of the total milk production, the Centrals provinces (Cocle, Los Santos, Herrera and Veraguas) and Chiriquí contribute about 85 % of national production (Cadena Agroalimentaria de Leche, 2017).
Currently, there is no type of publication or reference at the national grades that indicates a physical-chemical characterization of the quality of the raw cow´s milk (RCM), which allows it to be correlated with the classification it receives according to current regulations. Actually, the concept of milk quality involves, among others, physicochemical variables such as density (D), pH, fat (F), solids not-fat (SNF), total solids (TS), lactose (L) and protein (P).
Milk has a wide variety of components such as minerals, vitamins, etc. These components can vary even in different milkings of the same cow, mainly due to natural causes such as: genetic factors, lactation stage, health status of the cow, diet, etc. In addition to these factors, it is possible that the exposure of the milk to the environment, an inadequate management of the cold chain, hygiene of the facilities and storage time, can influence the characteristics of raw milk (Bernabucci et al., 2015; Hogenboom et al., 2019; Linn, 1988; Walstra et al., 2006).
There is the possibility of correlating the classification that RCM currently receives with the physicochemical profile of some of the characteristics considered important in milk. Given that in Panama there are no bibliographic or statistical records that allow the traceability and characterization of these quality characteristics of the milk produced in the country to be established, which incidentally, are carried out only in some cases and only by the receiving companies of milk for processing. The objective of this work was characterize and quantify the physical-chemical parameters of the quality of raw cow’s milk and correlate these characteristics with the classification that milk receives in Panama.
In order to maintain a true representation of the diversity of existing producers, the minimum representative sample necessary was calculated through stratified sampling, considering the classification of milk grade A, grade B and grade C or Industrial, according to the presence or not of bacteria: a) grade A milk with <200,000 bacteria ml-1, b) grade B milk with <500,000 bacteria ml-1, and c) grade C milk, milk that does not meet the parameters required by grade A and B milks (Ministerio de Gobierno y Justicia, 1970) as a stratum. With a confidence level of 95 %, a standard deviation of 5 % and a limit of error of 5 %. These considerations identified the need to have a minimum of 359 samples, 30 from grade A, 44 from grade B and 285 from grade C. However, a total of 512 milk samples were collecgted, distributed as follows: 48 grade A, 71 grade B, and 393 grade C. These samples were randomly selected in the provinces of Cocle (A = 1, B = 3, C = 14), Veraguas (A = 3, B = 4, C = 57), Herrera (A = 3, B = 7, C = 110) and Los Santos (A = 3, B = 12, C = 97), during the rainy season corresponding to the months of June to September 2018. In the Chiriqui province (A = 38, B = 45, C = 115), the sampling was carried out during the dry season, corresponding to the months from January to March 2019. This selection was made through a call to existing producers and processing companies in the different Provinces. The type of milk produced (A, B and C) was previously classified by the collection centers of the raw milk processing companies. If samples were obtained through tours of various farms or at collection points for processing plants. It is necessary to point out that the majority of dairy producers that have a great diversity of breeds in their herds, so this aspect was not considered. They were collected and kept preserved in coolers prepared for this purpose, at temperatures between 4 and 6 degrees Celsius (Martínez et al., 2007) in 30 ml sterile plastic bottles.
Once the interested producers have been identified, visits will be made to the producing farms in order to carry out the specific analyzes, these measurements contemplate evaluating the following variables: density (D), pH, fat (F), solids not-fat (SNF), total solids (TS), lactose (L), and protein (P).
The physico-chemical parameters under study were measured using the Speedy Lab analyzer, based on ultrasonic technique, sensitivity, specificity (Astori Tecnica snc, Brescia, Italy). It only requires conditioning the sample temperature (35 to 40 C°) for its measurement. The pH measurement was performed with the Oakton Instruments pH 11 meter, Vernon Hills, USA.
The data obtained was subjected to various statistical processes to obtain the best interpretation, among them, analysis of variance (ANOVA), correlation test, individually and cluster evaluation using the milk grade and geographic origin as differential parameters. For this, the statistical software such as:Statistical Package for the Social Sciences 15.0 (SPSS) and with the software XLSTAT 2019.2.1 were used.
The samples were obtained in two ways, collected directly from the production farms and from collection centers of milk processing companies. A total of 512 milk samples were collected, distributed as follows: 48 samples grade A, 71 samples grade B and 393 samples grade C.
There is a wide variety of works and publications that deal with the quality of RCM, but, they do not show standardized values that indicate specific ranges that producers must meet to consider that the milk produced can be considered or not as quality milk. It was possible to generate a table of ranges of suggested values for each of the physical-chemical characteristics evaluated (Table 1).
Table 1. Mean values of the parameters studied in raw milk according to the class of raw cow’s milk (A, B and C) and the global average (mean). The samples were collected in the provinces of Cocle, Herrera, Los Santos, Veraguas, from June to September 2018 and in the province of Chiriqui, from January to March 2019. Panama Republic. Range of values of the physical-chemical parameters to be evaluated according to the bibliography.
Cuadro 1. Valores medios de los parámetros estudiados en la leche cruda según la clase de leche cruda de vaca (A, B y C) y el promedio global (mean). Las muestras se recolectaron en la provincia de Coclé, Herrera, Los Santos, Veraguas, de junio a septiembre 2018 y en la provincia de Chiriquí, de enero a marzo de 2019. República de Panamá. Rango de valores de los parámetros físico-químicos a evaluar según la bibliografía.
The descriptive analysis of the data obtained, allowed to establish ranges of indicative values for each one of the studied physico-chemical parameters. As can be seen in Table 1, the values of solids non-fat, total solids and proteins were very close to the minimum reference values established. As expected, these grades influence the density results, which changes with a tendency towards minimum values. On the other hand, fats and lactose showed values with a tendency towards the upper limits considered in the bibliography. Considering that Panama does not have an official standard that establishes adequate ranges for these milk parameters, ranges of values identified from a previous bibliographic review were pre-established and based on standards from other countries and regions.
For all the physical-chemical parameters, it was possible to observe a tendency for these values to decrease as the grades of raw milk decreases. In the case of FAT, the values vary from 4.61 %, 4.12 %, and 4.08 % for grade A, B, and C, respectively. For the SNF, 9.19 %, 8.32 %, and 7.92 %, respectively for A, B and C. In respective order of A, B, and C, the TS was 13.80 %, 12.44 %, and 12.0 %. The D showed values in the same order, respectively for A, B, and C, 1032, 1029 and 1028 kg m3. The P was 3.36 %, 3.04 %, and 2.90 % for milk A, B, and C. Finally, the L presented the same behavior, 5.05 %, 4.57 %, and 4.36 %. The pH maintains constant values, of 6.56, 6.58, and 6.55. This indicates that the grades of the physico-chemical conditions influence the quality of the RCM.
The measurements of the parameters studied reflect very close values between the classes of raw milk produced, for which an evaluation of the means (ANOVA) was carried out with the combined effect of the class of milk produced. As can be seen in Table 2, in all cases the F statistic and its level of significance (p<0.05) that accepts the validity of the hypothesis of equality of the means, that is, under the conditions of this study, no significant differences were found between the milk grades considered. This condition may be due to the high percentage of diversity of breeds and crosses existing in a large part of the dairy herds involved.
Table 2. Analysis of variance with the combined effect of the class of raw cow’s milk produced. Milk grade A, grade B, and grade C or Industrial, according to the presence or not of bacteria a) grade A milk with <200,000 bacteria ml-1, b) grade B milk with <500,000 bacteria ml-1, c) grade C milk, milk that does not meet the parameters required by grade A and B milks (Ministerio de Gobierno y Justicia, 1970). The samples were collected in the provinces of Cocle, Herrera, Los Santos, Veraguas, from June to September 2018 and in the province of Chiriqui, from January to March 2019. Panama Republic.
Cuadro 2. Análisis de varianza con el efecto combinado de la clase de la leche cruda de vaca producida. Leche grado A, grado B y grado C o industrial, de acuerdo con la presencia o no de bacterias a) grado A <200 000 bacterias ml-1, b) grado B con <500 000 bacterias ml-1, c) grado C, la leche que no reúne los parámetros requeridos por los grados A y B (Ministerio de Gobierno y Justicia, 1970). Las muestras se recolectaron en la provincia de Coclé, Herrera, Los Santos, Veraguas, de junio a septiembre 2018 y en la provincia de Chiriquí, de enero a marzo de 2019. República de Panamá.
Given the result of the absence of significant between variables regarding grade classification, the data was subjected to Pearson correlation analysis. In the Table 3, shows the grades resulting similarities, the P has a high level of correlation with the SNF, with the D and L (0.93, 0.91, and 0.96, respectively), and a positive relationship intermediate with TS (0.65), conditions to be expected in this type of matrix.
Table 3. Pearson Correlation Matrix of the parameters obtained from the samples of raw cow’s milk analyzed. Provinces of Cocle, Herrera, Los Santos, Veraguas, from June to September 2018 and in the Province of Chiriquí, from January to March 2019. Panama Republic.
Cuadro 3. Matriz de correlación de Pearson de los parámetros obtenidos de las muestras analizadas de leche cruda de vaca analizadas. Provincias de Coclé, Herrera, Los Santos, Veraguas, de junio a septiembre 2018 y en la Provincia de Chiriquí, de enero a marzo de 2019. República de Panamá.
L was very strongly related to SNF and D (0.94, and 0.91, respectively. A positive intermediate relationship can also be observed between the TS with SNF, D, P, and L (0.70; 0.43; 0.65, and 0.67, respectively). In all cases, it can be seen that the pH values and the grades to which belong the samples of milk reflected almost no interaction.
A negative intermediate relationship showed that relates the grades with SNF (-0.50), D (-0.48), P (-0.50), and L (-0.50). Under the conditions carried out in this work, it is not possible to establish significant correlations between the physical-chemical characteristics evaluated that can differentiate the classes of milk, or the place where the samples come from.
It is necessary to consider that due to the notable complexity of the sample and the large number of variables that may intervene in the characteristics of milk production (species, diet, climate, geographic location, animal health, etc.) (Allore et al., 1997; Dominguez-Salas et al., 2019; Food and Agriculture Organization, 2013; Manterola, 2011). Much of the objective of the study was achieved, to characterize and quantify the physical-chemical parameters of the quality of raw cow’s milk in the different provinces studied.
The results obtained were consistent with what was observed by Hogenboom et al. (2019), where he points out that compared milk samples obtained with automatic milking system and traditional milking systems in herds of different sizes, in different lactation stages and in different periods of year, found that the milking system did not significantly affected the F, P, L, SNF, and casein. This indicates that what was found in the study is possible, considering that the samples used represented a homogenization of the milk of the herds, possibly involving animals of different breeds, lactation states, health, diet, factors already defined as those responsible for the conditions. physical-chemical and microbiological quality of raw cow’s milk.
The behavior of the results obteined from fat, SNF, TS, density, P, L and pH not present significant differences between the different grades of raw milk, it did represented an indication of the possible impact that the different variables studied represent on the quality grade of raw cow’s milk. Condition mentioned in different studies carried out (Bernabucci et al., 2015; Hogenboom et al., 2019; Linn, 1988; Walstra et al., 2006).
Additionally, according to the values found in the present study, it was possible to affirm that RCM in the regions under study have the profile that adjusts to the wide variety of studies carried out and regulations used by different countries (Comunidad Andina, 2009; East African Community, 2006; European Economic Community, 1992). Since, it was not possible to consider that a cross-sectional analysis collects enough information to discriminate the quality of the RCM within a specific classification.
The current regulations require deepening the microbiological aspect and consider novel aspects such as the presence of somatic cells as an interpretation of the health status of the cow’s udder, without neglecting the physical-chemical condition of the same, such as tools to establish the quality grades of the milk produced.
All the values of the different variables analyzed were within the ranges of values pre-established as satisfactory within the bibliographic review and existing regulations in other countries.
In the conditions carried out in this study, the values found, showed a slight tendency that agrees with the quality grade of the RCM, did not present significant differences that allow them to be considered as descriptors of the quality. The correlation analysis showed that there are strong positive relationship between the parameters studied.
Having carried out a sample greater than the minimum suggested according to the methodology, the results do not yield conclusive information to significantly validate the distribution of each sample in the respective current classification. It is advisable to carry out a broaders studies to establish standardized indicators that allow identifying and locating national milk as a quality product, considering sampling in consecutive periods to establish statistical significance for different climatic seasons (rainy and dry), genetic characteristics, health states herds, geographical areas, etc.
The authors thank the National Research System (SNI) for financing the initial phase of the Project. To the National Secretariat of Science, Technology and Innovation (SENACYT) for financing the Project (IAP17-P-047) through the Public Call for Innovation Projects to the Agricultural Sector of Panama, 2017. Special thanks to the producers and companies that facilitated the infrastructure and collaboration for the collection of the samples, among them, Ramón Batista, Arturo Fuentes, Prolacsa, Quesos Doña Mery, Estrella Azul, and many others.
Allore, H. G., Oltenacu, P. A., & Erb, H. N. (1997). Effects of Season, herd size, and geographic region on the composition and quality of milk in the Northeast1. Journal of Dairy Science, 80(11), 3040-3049. http://www.sciencedirect.com/science/article/pii/S0022030297762714
Bernabucci, U., Basiricò, L., Morera, P., Dipasquale, D., Vitali, A., Piccioli Cappelli, F., & Calamari, L. (2015). Effect of summer season on milk protein fractions in Holstein cows. Journal of Dairy Science, 98(3), 1815–1827. https://doi.org/10.3168/jds.2014-8788
Bondan, C., Folchini, J. A., Noro, M., Machado, K. M., Muhls, E., González, F. H. D., Bondan, C., Folchini, J. A., Noro, M., Machado, K. M., Muhls, E., & González, F. H. D. (2019). Variation of cow’s milk composition across different daily milking sessions and feasibility of using a composite sampling. Ciência Rural, 49(6), Article 2018004. https://doi.org/10.1590/0103-8478cr20181004
Cadena Agroalimentaria de Leche. (2017). Desafíos para la consolidación de la cadena de leche. Ministerio de Desarrollo Agropecuario. https://www.mida.gob.pa/upload/documentos/leche2012.pdf
Closa, S. J., de-Landeta, M. C., Andérica, D., Pighín, A., & Cufré, J. A. (2003). Contenido de nutrientes minerales en leches de vaca y derivados de Argentina. Archivos Latinoamericanos de Nutrición, 53(3), 320-324. http://ve.scielo.org/scielo.php?script=sci_abstract&pid=S0004-06222003000300016&lng=es&nrm=iso&tlng=es
Comunidad Andina. (2009). Norma Andina 0063:2009. Leche Cruda. Requisitos. Comunidad Andina.
Dominguez-Salas, P., Galiè, A., Omore, A., Omosa, E., & Ouma, E. (2019). Contributions of milk production to food and nutrition security. In P. Ferranti, E. M. Berry, & J. R. Anderson (Eds.), Encyclopedia of Food Security and Sustainability (pp. 278–291). Elsevier. https://doi.org/10.1016/B978-0-08-100596-5.21526-6
East African Community. (2006). East African standard. Raw cow milk-Specification. https://law.resource.org/pub/eac/ibr/eas.67.2006.html
European Economic Community. (1992, September 14). Council directive 92/46/EEC laying down the health rules for the production and placing on the market of raw milk, heat-treated milk and milk-based products. Official Journal of the European Communities.
Food and Agriculture Organization of the United Nations. (2013). Dairy production and products: milk composition. http://www.fao.org/dairy-production-products/products/milk-composition/en/
Food and Agriculture Organization of the United Nation. (2020). FAOSTAT statistical database. http://www.fao.org/faostat/en/#data/QC
Franzoi, M., Niero, G., Visentin, G., Penasa, M., Cassandro, M., & De-Marchi, M. (2019, abril). Variation of detailed protein composition of cow milk predicted from a large database of mid-infrared spectra. Animals, 9(4), Article 176. https://www.mdpi.com/2076-2615/9/4/176
Gallardo, M. R., Onetti, S. G., Castillo, A. R., & Nari, J. O. (1996). Proteína en leche y su relación con el manejo nutricional (Informe Técnico No 56). Instituto Nacional de Tecnología Agropecuaria.
Hogenboom, J. A., Pellegrino, L., Sandrucci, A., Rosi, V., & D’Incecco, P. (2019). Invited review: Hygienic quality, composition, and technological performance of raw milk obtained by robotic milking of cows. Journal of Dairy Science, 102(9), 7640–7654. https://doi.org/10.3168/jds.2018-16013
Instituto Nacional de Estadística y Censo. (2020). Estimación de la producción de leche de vaca en la república. https://www.inec.gob.pa/archivos/P7421312-13.pdf
Jay, J. M. (2009). Microbiología moderna de los alimentos. (5a Ed.). Editorial ACRIBIA S.A.
Linn, J. G. (1988). Factors affecting the composition of milk from dairy cows. National Academies Press. https://www.ncbi.nlm.nih.gov/books/NBK218193/
Manterola, H. (2011). Manejo nutricional y composición de la leche el desafío de incrementar los sólidos totales en la leche: Una necesidad de corto plazo. VirtualPro.co. https://www.virtualpro.co/biblioteca/manejo-nutricional-y-composicion-de-la-leche-el-desafio-de-incrementar-los-solidos-totales-en-la-leche-una-necesidad-de-corto-plazo
Martínez, L. R. B., Ortega, A. E., Flores, J. E., & Ortega, O. A. C. (2007). Determinación de la calidad fisicoquímica de la leche cruda producida en sistemas campesinos en dos regiones del Estado de México. Veterinaria México, 38(4), 395–407. https://www.redalyc.org/articulo.oa?id=42338402
Ministerio de Gobierno y Justicia (1970, febrero 4). Decreto de Gabinete 229 del 16 de julio de 1969, por el cual se dictan algunas disposiciones relacionadas con el aspecto sanitario y calidad de la leche y de los productos lácteos y se subrogan algunos aertículos del Decreto N° 256 del 13 de junio de 1962. Gaceta Oficial del Estado https://www.gacetaoficial.gob.pa/gacetas/16537_1970.pdf
Ocampo, R., Gomez, C., Retrepo, D., & Cardona, H. (2016). Estudio comparativo de parámetros composicionales y nutricionales en leche de vaca, cabra y búfala, Antioquia, Colombia. Revista Colombiana de Ciencia Animal, 8(2), 177–186. https://doi.org/10.24188/recia.v8.n2.2016.185
Oliszewski, R., Rabasa, A. E., Fernández, J. L., Poli, M. A., & Núñez-de-Kairúz, M. S. (2002). Composición química y rendimiento quesero de la leche de cabra Criolla Serrana del noroeste argentino. Zootecnia Tropical, 20(2), 179-189. http://ve.scielo.org/scielo.php?script=sci_abstract&pid=S0798-72692002000200003&lng=es&nrm=iso&tlng=es
Ponce, P. (2009). Composición láctea y sus interrelaciones: Expresión genética, nutricional, fisiológica y metabólica de la lactación en las condiciones del Trópico. Revista de Salud Animal, 31(2), 69-76. http://scielo.sld.cu/scielo.php?script=sci_abstract&pid=S0253-570X2009000200001&lng=es&nrm=iso&tlng=es
Walstra, P., Wouters, J., & Geurts, T. (2006). Dairy science and technology (2nd Ed.). Taylor & Francis Group. https://www.routledge.com/Dairy-Science-and-Technology/Walstra-Walstra-Wouters-Geurts/p/book/9780824727635