SciELO - Scientific Electronic Library Online

 
vol.55 número1CunICA, herramienta informática para la gestión de los procesos de cría en granjas cunículas. Nota técnicaEfecto de la inclusión de materiales fibrosos en la fermentación en estado sólido de residuos poscosecha de Solanum tuberosum, inoculados con un preparado microbiano índice de autoresíndice de materiabúsqueda de artículos
Home Pagelista alfabética de revistas  

Servicios Personalizados

Revista

Articulo

Indicadores

  • No hay articulos citadosCitado por SciELO

Links relacionados

  • No hay articulos similaresSimilares en SciELO

Compartir


Cuban Journal of Agricultural Science

versión impresa ISSN 0864-0408versión On-line ISSN 2079-3480

Cuban J. Agric. Sci. vol.55 no.1 Mayabeque ene.-mar. 2021  Epub 01-Mar-2021

 

Animal Science

Effect of temperature and humidity index (THI) on the physiological responses of grazing dairy cows

J. C. Valdivia-Cruz1  * 

J. J. Reyes-González2 
http://orcid.org/0000-0002-9681-1187

G.R. Valdés-Paneque3 
http://orcid.org/0000-0002-4979-9082

1Asociación Cubana de Producción Animal (ACPA), Provincial Sancti Spíritus

2Instituto de Ciencia Animal, Apartado Postal 24, San José de las Lajas, Mayabeque, Cuba

3Universidad de Sancti Spíritus José Martí Pérez, Sancti Espíritus, Cuba

Abstract

The evaluation was performed with 25 cows from dairy 17, belonging to Dos Ríos enterprise, in order to assess the relation between the temperature and humidity index, rectal temperature and respiratory rate. The information was taken at three times of the day: morning, midday and afternoon. The results showed that in the morning the cows are in a thermo neutral zone (temperature and humidity index 69.59), while at midday they are under medium heat stress (THI 79.61), and in the afternoon under severe stress (THI 91.69). Under severe heat stress, the rectal temperature of pregnant lactating cows increased (P = 0.0136) by 0.280C with respect to the rectal temperature of non-pregnant lactating cows. The respiratory rate was not affected by the light temperature and humidity index. But in the moderate stress, pregnant cows showed higher (P = 0.0003) respiratory rates of 1.9 and 2.0 % in relation to non-pregnant cows for moderate and severe stress, respectively. The physiological response to the temperature and humidity index, according to the milking group, showed that the high production group had higher rectal temperature values (P = 0.0003) in 0.16 and 0.21 oC and the respiratory rate increased (P = 0.0024) in 9.21 and 7.89 % in relation to the average group, for moderate and severe stress, respectively. Under the conditions of this study, at midday and afternoon, the animals are under moderate and severe stress, respectively, conditions that affect more the lactating pregnant animals and the open cows or the highest production group.

Key words: rectal temperature; respiratory rate; milk production

Cuban livestock faces the great problem of the low efficiency of dairy cattle, which is due to limitations in feeding and management, as well as to climatic conditions, mainly high temperatures and relative humidity (Álvarez 2006). If the environmental temperature reaches values above 35°C for animals of tropical origin, there is a failure in the thermoregulation systems, and the rectal temperature and respiratory rate increase, join to the decrease of food intake and milk production (Corrales 2014 and Salvador 2015). In this way, the animal faces the environmental conditions of the place, using its physiological capacities to adapt to the environment (West 2003 and Polsky and von Keyserlingk 2017).

The stress term is currently used in a much broader sense for its application to animal welfare (Sejian et al. 2015). That is why codes that regulate animal welfare have been published, such as the sanitary code for terrestrial animals (OIE, 2017). In its chapter 7.11.4, this code defines the measurable criteria of animal welfare, where it includes as indicators of negative behavior the decrease in dry matter intake, the increase in respiratory rate and other abnormal behaviors, related to the effect of management the animals are submitted. The objective of this study is to evaluate the effect of temperature and humidity index (THI) on the physiological and productive indicators of Siboney de Cuba crossbred cows in dairy 17 from Dos Ríos Enterprise (E).

Materials and Methods

The study was carried out in dairy 17, from Dos Ríos Enterprise, belonging to Managuaco agricultural enterprise, located at 210 56 ’N and 790 20’ W, in Sancti Spíritus municipality, Sancti Spíritus province. The predominant soil in this region is soft brown carbonate (Hernández et al. 2015).

A total of 25 cows from Siboney de Cuba breed were used, in milking, in different lactation phases, distributed in three groups according to lactation days: a) lactation initiation (high group, 38.4 ± 11.3 d of average lactation), b) middle lactation (average group, 105.7 ± 22.5 d of average lactation) and c) lactation end (low group, 182.3 ± 31.6 d of average lactation). The animals were between the second and third lactation.

The physiological states of gestation of cows were obtained from the reproductive control cards, once the pregnancy was confirmed three months after insemination and when they did not present heat.

Environmental measurements (relative humidity and environmental temperature) and physiological responses (respiratory rate and rectal temperature) were taken at three times of the day: in the grazing, between 7:00 and 9:00 h; in the shade buildings between 11:30 and 12:30 h, and in the milking parlor between 15:00 and 16:30 h.

The measurements of environmental temperature (ET) in degrees centigrade (0C) and relative humidity (RH) in percentage value (%) were taken at 1.20 m altitude. Likewise, the meteorological data (ET and RH) were recorded at the Sancti Spíritus Meteorological Station, coinciding with the measurements date in the dairy, so that the respective comparisons could be made.

To determine the temperature and relative humidity index of the environment (THI), the formula proposed by García et al. (2007) was used:

[TeX:] THI=0.81*ET + (ET -14.4)* RH/100 + 46.4

In the animals, rectal temperature was determined by introducing a maximum clinical thermometer into the rectal vial, free of fecal matter, to eight centimeters depth. The reading was carried out after four minutes. To avoid possible negative effects on dairy production with the management proposed for the measurement, the animals were adapted for a previous 15 d before starting the procedure. Respiratory rate was obtained by direct observation of the animal and by counting the movements of the flanks for one minute.

Milk production was individually measured in both milkings (morning and afternoon) and was taken as daily yield, with biweekly frequencies.

The statistical analysis was carried out according to a general linear model with nested effect, where the gestation or not and the milking groups are nested at each moment of the day in which measurements were made. In the necessary cases, Duncan's test (1955) was applied.

The main effect of the THI interactions and the milking groups and THI and pregnancy or not is shown when it is significant.

For milk production, analysis of variance was carried out, according to a multiplicative model for pregnancy or not and milking groups, respectively. In the necessary cases, Duncan's test (1955) was applied for P <0.05.

[TeX:] Yijklm = Anb exp cn + Ti + Sj + (THI)k + GOl + eijklm

Where:

Yijklm =

milk production kg/v/d

Anb =

effect of the lactation curve

Ti =

effect of i-th treatment

Sj=

effect of the j-th sample

(THI)k=

effect of the k-th temperature-humidity index

GOl=

effect of the l-th gestation or milking group.

eijklm =

normal random error distributed with mean 0 and variance σ2

The statistical package used was the InfoStat, version 2012 (Di Rienzo et al. 2012)

Results and Discussion

Due to the high heat stress ratio in dairy cows, with THI values outside their thermoneutrality range, which can significantly reduce the ability to lose latent heat by evaporation, it is possible to use this index as an estimate of heat stress parameters. (Cerqueira et al. 2013, Callejo 2015 and Ruiz et al. 2019).

In a general sense, it is considered that, with THI values lower than 72.0 units, the animals are in thermoneutral conditions, or it is estimated that there are no stressful conditions. Faced with figures between 72.0 and 79.0, the presence of light stress is considered; from 80.0 to 89.0 is understood to be moderate, and values from 90.0 to 99.0 show severe stress. Higher figures are considered serious, and can limit the animals lives (Anzures at al. 2015 and Veissier et al. 2018).

The environmental conditions (figure 1), mainly the temperature registered in the sampling times, are typical in the province. They are characterized by being relatively low in the morning (25.8 ± 0.030C), they are highest at midday (30.1 ± 0.030C) and high in the afternoon (34.6 ± 0.040C). These temperatures, combined with high humidity (average of 91 to 93 %) throughout the study period, lead to THI values of 69.59 ± 1.59 units (light stress), 79.61 ± 1.26 units (moderate stress) and 91.69 ± 1.24 units (severe stress), for the morning, midday and afternoon, respectively.

Figure 1 Climatic variables of environmental temperature (ET), relative humidity (RH) and temperature humidity index (THI), according to sampling hours 

In general, the climatic information suggests that in the morning the cows were in a thermoneutral zone and under acceptable conditions to express their potential for milk production. At midday and in the afternoon, they were under moderate and severe heat stress respectively, since the dairy cow begins to experience the consequences of heat stress when the THI exceeds 72 units (Anzures et al. 2015, and Polsky and von Keyserlingk 2017).

Table 1 shows the interactions between the different stress levels and the physiological responses of rectal temperature (P = 0.0136) and respiratory rate (P = 0.0003) of pregnant and non-pregnant lactating cows. With these results it is shown that, regardless of the reproductive status of the animals, as the value of the THI increases, these two physiological variables increase.

When assessing the reproductive status (pregnant or not) of lactating dairy cows, in the values of light and moderate stress, it was found that there are no differences in the physiological response of rectal temperature. However, in pregnant animals with severe stress there was a higher (P = 0.0136) rectal temperature than in non-pregnant animals in 2.73 %.

The analysis of the effect of increasing the value of stress on both reproductive states showed (P = 0.0136) that rectal temperature increased by 1.62 and 2.73, 1.84 and 1.85%, when passing from light to moderate stress, and from the latter to severe , for pregnant and non-pregnant animals, respectively.

The interaction between stress and respiratory rate of lactating cows, being pregnant or not, does not differ under light stress conditions. In this case, the animals showed respiratory rate values in the normal range established for the species (Anzures et al. 2015 and Allen et al. 2015). However, from moderate stress, pregnant cows had higher (P = 0.0003) respiratory rates of 1.96 and 2.00 % than non-pregnant cows, for moderate and severe stress, respectively. Similarly, regardless of the physiological state of cows, when changing from light to moderate stress conditions, and from the latter to severe stress, the respiratory rate increased (P = 0.003) in 42.58, 20.74 and 39.21, 20.69 % for pregnant and non-pregnant, respectively.

Table 1 Interaction between heat stress and pregnancy or not, in rectal temperature and respiratory rate of dairy cows in production 

Indicator

  • Strees/

  • Reproductive state

Light Moderate Severe Sig.
Rectal temperature (0C) Pregnant 38.96 a ±0.04 39.59 b ±0.04 40.67 d ±0.04 P=0.0136
Non-pregnant 38. 94 a ±0.04 39.66 b ±0,07 40.39 c ±0.04
Respiratory rate (No min.-1) Pregnant 53. 27 a ±0,28 75.95 c ±0,45 91.70 e ±0,27 P=0.0003
Non- pregnant 53.51 a ±0.27 74.49 b ±0.27 89.90 d ±0.28

a, b, c, d and e. Means with different letters differ at P<0.05 (Duncan 1955)

Table 2 shows the interactions of heat stress with rectal temperature (P = 0.0003) and respiratory rate (P = 0.0024) of dairy cows, according to the production group. These show that regardless of the production group of animals, as the stress value increases, these two physiological variables increase.

In the analysis of rectal temperature, face to light stress values, there were not effects on the physiological response of cows, regardless of the milking group. However, face to higher stress values, the cows from the high production group showed higher rectal temperature figures (P = 0.0003) in 0.16 and 0.21°C with respect to those of the medium group, for moderate and severe stress, respectively.

Regardless of stress, the cows in the low production group do not differ from those in the other two production groups, which may be due to the higher proportion of pregnant animals found in said group, since they end lactation. As could be observed in table 1, by itself the reproductive gestation state affects the animals equally, be under moderate or severe stress.

The rectal temperature variable of the animals according to the milking group increases (P = 0.0003) when going from light to moderate stress conditions, and from the latter to severe stress in 2.92, 1.73; 2.24, 1.61 and 2.55, 1.73 % for the high, medium and low milking groups, respectively.

Table 2 Interaction between stress and milking group on rectal temperature and respiratory rate of dairy cows in production 

Indicator

  • Stress/

  • Milking group

Light Moderate Severe Sig.
Rectal temperature (0C) High 38. 76 a ±0.05 39.89 c ±0.05 40.58 e ±0.05 P=0.0003
Medium 38. 86 a ±0.06 39.73 b ±0.06 40.37 d ±0.06
Low 38. 83 a ±0.06 39.82 bc ±0.06 40.51 de ±0.06
Respiratory rate (No min.-1) High 53.65 a ±1.25 77.78 c ±1.25 93.98 e ±1.25 P=0.0024
Medium 53.07 a ±1.28 71.22 b ±1.27 87.11 d ±1.28
Low 52.64 a ±1.32 74.92 bc ±1.32 91.08 de ±1.33

a, b, c, d, e. Means with different letters differ at P<0.05 (Duncan 1955)

Likewise, the lactating cows, faced to light stress effects, did not show affectations in the respiratory rate, regardless of the milking group (table 2). However, this physiological response increased (P <0.0024) by 9.21 and 7.89 % for the high production group compared to the medium group, under moderate and severe stress, respectively. There was no difference between the low production group with respect to the other two.

These affectations in the physiological variables studied under moderate or severe stress in the animals from the high production group are in agreement with what is reported in the literature. Authors such as La Manna et al. (2014) and Leucen (2014) report that cows at the beginning of lactation are more affected by high THI values, since they showed higher metabolism due to high productive levels.

The variable of the respiratory rate of the animals increased (P = 0.0024) independently of the milking group, when going from the conditions of light to moderate stress, and from the latter to severe stress in 44.98, 20.83; 34.20, 22.31 and 42.32, 21.56 % for the high, medium and low milking groups, respectively.

The respiratory rate in thermoneutral conditions has been reported to fluctuate between 40 and 55 frequencies min. for dairy cattle (Anzures et al. 2015 and Allen et al. 2015). The increase in body temperature and respiratory rate are normal mechanisms, by which cows dissipate heat to maintain their homeothermic condition in hot environments, being about 38.5 °C the body temperature considered normal in cattle (Ganaie et al. 2013).

Thus, the respiration rate is one of the most important mechanisms to consider when evaluating the heat stress level in cattle (Beauregard et al. 2018), since it is one of the main observable responses in the animal when it is exposed at temperatures that exceed your comfort threshold. This increase in the respiration rate is intended to increase respiratory heat loss, which is one of the most important for maintaining thermal balance.

Brown-Brandl et al. (2005) showed that measuring the respiration rate of animals, determining that they are in the panting process, as well as quantifying it, is the easiest and most affordable way to assess heat stress in cattle under commercial production conditions, since this does not require sophisticated equipment.

Furthermore, unlike body temperature, the response is practically immediate in the animal, and it follows almost the same pattern as environmental temperature.

These heat stress conditions lead directly to the activation of thermoregulatory, physiological and metabolic mechanisms, which cause a decrease in milk production, higher dissipation of body heat and a reduction in the metabolic heat production (Ganaie et al. 2013, Barragán et al. 2015, Beauregard et al. 2018 and Kamal et al. 2018).

When analyzing the effect of milk production according to the production group (table 3), it was found that the animals with high production produced more (P = 0.0001) milk 1.47 and 2.18 times than those with medium and low, respectively. Likewise, the average cows produced more milk (P = 0.00001) in 48.01 % than the low animals.

These milk production values according to the production group are logical, since the animals are mainly grouped by their lactation days. Cows in early lactation are in the high production group, which is why they have higher potential for milk production and therefore, a higher metabolism to achieve these productions (NRC 2001, Corrales 2014). However, as table 2 shows, this group is more susceptible to stressful environments. The adoption of measures to try to reduce this negative effect, improving management and feeding, will make possible to increase milk production (Mendoza et al. 2016).

Table 3 Effect of the milking group on the milk production of lactating dairy cows 

  • Milking group/

  • Indicator

High Medium Low Sig.
Milk prod kg cow-1 day-1 SE 2.62 c (1) (14.24) ±0.05 2.25 b (1) (9.68) ±0.03 1.80 a (1) (6.54) ±0.04 P=0.0001

Means transformed according to ln () Original means

a, b, c. Means with different letters differ at P<0.05

This research allowed assessing that under the described conditions, the midday and afternoon hours cause stress in the animals, regardless they were in the shade buildings. The effects on rectal temperature and respiratory rate under conditions of moderate or severe stress are more marked in lactating animals in gestation and in the group of high or highest production.

Acknowledgments

Thanks to the support provided by the technicians and specialists of Dos Ríos Enterprise, from “Managuaco” enterprise, as well as the technical staff of the Biomathematics group of the Institute of Animal Science.

References

Allen, J.D., Hall, L.W., Collier, R.J. & Smith, J.F. 2015. "Effect of core body temperature, time of day, and climate conditions on behavioral patterns of lactating dairy cows experiencing mild to moderate heat stress". Journal of Dairy Science, 98(1): 118-127, ISSN: 0022-0302, DOI: http://dx.doi.org/10.3168/jds.2013-7704. [ Links ]

Álvarez, A. 2006. Análisis regional de la cadena de lácteos en Querétaro. Memorias VII Congreso Latinoamericano de Sociología Rural. Quito, Ecuador. [ Links ]

Anzures, F., Macías, U., Álvarez, F.D., Correa, A., Díaz, R., Hernández, J.A. & Avendaño, L. 2015. "Effect of season (summer vs. winter) on physiological variables, milk production and antioxidant capacity of Holstein cows in an arid zone of northwestern Mexico". Archivos de Medicina Veterinaria,47(1):15-20, ISSN: 0301-732X. [ Links ]

Arias, R.A., Mader, T.L. & Escobar, P.C. 2008. "Climatic factors affecting cattle performance in dairy and beef farms". Archivos de Medicina Veterinaria, 40(1): 7-22, ISSN: 0301-732X. [ Links ]

Barragán, W.A., Mahecha, L. & Cajas, Y.S. 2015. "Variables fisiológicas metabólicas de estrés calórico en vacas bajo silvopastoreo y pradera sin árboles". Agronomía Mesoamericana, 26(2): 211-223, ISSN: 2215-3608, DOI: http://dx.doi.org/10.15517/am.v26i2.19277. [ Links ]

Beauregard, J., Prado, O., García, L., García, A., Macedo, R. & Hernández, J. 2018. "Productivity of Holstein dairy cows without shade at two seasons of the year". Abanico Veterinario, 8(3): 51-67, ISSN: 2448-6132, DOI: http://dx.doi.org/10.21929/abavet2018.83.3. [ Links ]

Brown-Brandl, T.M., Eigenberg, R.A., Nienaber, J.A. & Hahn, G.L. 2005. "Dynamic response indicators of heat stress in shaded and non-shaded feedlot cattle, Part 1: Analyses of indicators". Biosystems Engineering, 90(4): 451-462, ISSN: 1537-5110, DOI: https://doi.org/10.1016/j.biosystemseng.2004.12.006. [ Links ]

Callejo, A. 2015. El estrés calórico en vacas lecheras. Available: http://oa.upm.es/37755/1/INVE_MEM_2015_201963.pdf. [Consulted: January 27th, 2017]. [ Links ]

Cerqueira, J.O.L., Araújo, J.P., Blanco-Penedo, I., Cantalapiedra, J., Silvestre, M., Silva, S.R., 2013. Estudio de indicadores fisiológicos como predictores de estrés térmico de vacas lecheras en Norte de Portugal. XV Jornadas sobre Producción Animal, Tomo I, pp.40-42, Zaragoza, España, ISBN: 978-84-695-7684-7. [ Links ]

Corrales, I.E. 2014. El efecto del estrés calórico en bovinos de leche. Diploma Thesis. Universidad Autónoma Agraria “Antonio Narro”, Torreón, Coahuila, México. [ Links ]

Di Rienzo, J.A., Casanoves, F., Balzarini, M.G., González, L., Tablada, M. & Robledo, C.W. 2012. InfoStat, Version 2012 (Windows). Grupo InfoStat, Universidad Nacional de Córdoba, Argentina. Available: http://www.infostat.com.ar. [ Links ]

Duncan, D.B. 1955. "Multiple Range and Multiple F Tests". Biometrics, 11(1): 1-42, ISSN: 0006-341X, DOI: https://doi.org/10.2307/3001478. [ Links ]

Ganaie, A.H., Shanker, G., Bumla, N.A., Ghasura, R.S., Mir, N.A., Wani, S.A. & Dudhatra, G.B. 2013. "Biochemical and physiological changes during thermal stress in bovines: A Review". Iranian Journal of Applied Animal Science, 3(3): 423-430, ISSN: 2251-628X. [ Links ]

García-Ispierto, F., Lòpez-Gatius, G., Bech-Sabat, P., Santolaria, J.L., Yàniz, C., Nogareda, F., De Rensis, M., López-Béjar. 2007. "Climate factors affecting conception rate of high producing dairy cows in northeastern Spain". Theriogenology, 67(8): 1379-1385, ISSN: 0093-691X, DOI: https://doi.org/10.1016/j.theriogenology.2007.02.009. [ Links ]

Hernández, A, Pérez, J.M., Bosh, D. & Castro, N. 2015. Clasificación de los suelos de Cuba. 1st Ed. Ed. Instituto Nacional de Ciencias Agrícolas, González, O. (ed.). San José de las Lajas, Mayabeque, Cuba, p. 93, ISBN: 978-959-7023-77-7. [ Links ]

Kamal, R., Dutt, T., Patel, M., Dey, A., Bharti, P.K. & Chandran, P.C. 2018. "Heat stress and effect of shade materials on hormonal and behavior response of dairy cattle: a review". Tropical Animal Health and Production, 50: 701-706, ISSN: 0049-4747, DOI: https://doi.org.10.1007/s11250-018-1542-6. [ Links ]

La Manna, A., Román, L., Bravo, R. & Aguilar, I. 2014. "Estrés térmico en vacas lecheras: con sombra y bienestar las vacas producen más". Revista INIA, 39(12): 34-39, ISSN: 1510-9011. [ Links ]

Lucena, C. 2014. "Heat stress effect parameters of lactation curves from dairy herd in the tropics". Gaceta de Ciencias Veterinarias, 19(1): 11-16, ISNN: 1690-8414. [ Links ]

Mendoza, A., Cajarville, C. & Repetto, J.L. 2016. "Short communication: Intake, milk production, and milk fatty acid profile of dairy cows fed diets combining fresh forage with a total mixed ration". Journal of Dairy Science, 99(3): 1938-1944, ISSN: 0022-0302, DOI: http://dx.doi.org/10.3168/jds.2015-10257. [ Links ]

NRC (National Research Council). 2001. Nutrient Requirements of Dairy Cattle. 7th Rev. Ed. Ed. National Academy Press, Washington D.C., U.S.A., ISBN: 978-0-309-06997-7, DOI: https://doi.org/10.17226/9825. [ Links ]

Polsky, L. & Von Keyserlingk, M.A.G. 2017. "Invited review: Effects of heat stress on dairy cattle welfare". Journal of Dairy Science, 100(11): 8645-8657, ISSN: 0022-0302, DOI: https://doi.org.10.3168/jds.2017-12651. [ Links ]

Ruiz, J.I., Vargas, B., Abarca, S. & Hidalgo, H.G. 2019. "Heat stress effect on dairy cattle production in Costa Rica". Agronomía Mesoamericana, 30(3): 733-750, ISSN: 2215-3608, DOI: https://doi.org/10.15517/am.v30i3.35984. [ Links ]

Salvador, A. 2015. Efectos del estrés calórico en vacas lecheras. Available: http://www.agritotal.com/seccion/informacion-tecnica/. [Consulted: December 28th, 2016]. [ Links ]

Sejian, V., Gaughan, J., Baumgard, L. & Prasad, C. 2015. Climate Change Impact on Livestock: Adaptation and Mitigation. Ed. Springer, New Delhi, India, p. 532, ISBN: 978-81-322-2265-1, DOI: https://doi.org/10.1007/978-81-322-2265-1. [ Links ]

Veissier, I., Van Laer, E., Palme, R., Moons, C.P.H., Ampe, B., Sonck, B., Andanson, S. & Tuyttens, F.A.M. 2018. ""Heat stress in cows at pasture and benefit of shade in a temperate climate region". International Journal of Biometeorology, 62(4): 585-595, ISSN: 1432-1254, DOI: https://doi.org/10.1007/s00484-017-1468-0. [ Links ]

West, JW, 2003. "Effects of heat-stress on production in dairy cattle". Journal of Dairy Science, 86(6): 2131-2144, ISSN: 0022-0302, DOI: https://doi.org/10.3168/jds.S0022-0302(03)73803-X. [ Links ]

Received: April 17, 2020; Accepted: October 06, 2020

*Email: jreyes@ica.co.cu

Conflict of interest: The authors declare that there are no conflicts of interests among them

Author´s contribution: J. C. Valdivia-Cruz: Design and conducting the experiment, data analysis, manuscript writing. J. J. Reyes-González: Design and conducting the experiment, data analysis, manuscript writing. G.R. Valdés-Paneque: Design and conducting the experiment, data analysis, manuscript writing.

Creative Commons License This is an open-access article distributed under the terms of the Creative Commons Attribution License