1. Introduction
The introduction of technological innovations for intensive cattle growing and the breeding of high-yielding cows is tightly connected with an increase in diseases incidence. Livestock diseases cause significant economic damage due to the high frequency of their propagation in the Russian Federation and abroad. For example, annually in Europe, up to 50% of cows are culled due to diseases of their extremities resulting in their low productivity and reduced reproductive capacity [
1,
2].
Potentially, a continuous non-invasive monitoring of cattle conditions based in various techniques can become a powerful tool for monitoring cattle’s physiological state on feedlots and in dairy production. Such monitoring can also allow the early detection of diseases; therefore, diagnoses will be made faster than with the aid of traditional methods of clinical assessment. Besides, the accuracy and the specificity of the traditional clinical assessment of cattle diseases leave much to be desired. On the other hand, the remote early detection of diseases with the aid of advanced technological systems can improve remarkably the accuracy and timeliness of disease detection. In addition, the continuous monitoring of cattle behavior based on advanced technology can provide important information useful for treatment determination or the recognition of a cow state as a temporal event in research of animals’ health and well-being [
3,
4,
5].
Changes in animals’ behavior represent one of the most important criteria for the assessment of animals’ well-being and health [
6,
7,
8]. The parameters of animals’ behavior can be used for the development of a system for the early prevention of diseases. For example, hoof injuries cause lameness, which make cows prefer laying instead of walking, thus decreasing their locomotive activity [
9,
10,
11].
A restless behavior of a sick animal can begin prior to the appearance of any of the main clinical symptoms of a disease. For example, cows react less vividly to external stimuli and their activity decreases a few hours before their temperature rises due to mastitis [
12,
13]. Besides, the circadian rhythm of cattle’s locomotive activity can decrease as early as two days before the clinical symptoms of this disease are detected [
14]. Hence, it is possible to use these behavioral deviations as early signs of diseases, detect and treat them at a very early stage, and thus reduce the cost of medicines for cattle.
Until recently, it was difficult to detect such changes in behavior on a regular basis; therefore, the employees of livestock farms had to monitor the animals continuously. Today, hi-tech sensors and tools are available for the implementation of the so-called “precision animal husbandry”. They enable continuous monitoring of the animals’ physiological state. Accelerometers can register acts of movement, positioning systems can detect animals’ positions, pH meters can measure the acidity of the animals’ rumen. and temperature sensors can provide caregivers with additional important health data [
15].
The rumen temperature measured orally (by a sensor in the animal’s mouth) has a strong correlation with the rectal temperature; thus, the rumen temperature can be considered as the body temperature [
16]. According to previous studies, the use of oral sensors for temperature measurement is a less invasive and yet reliable method to detect an increase of temperature in the cow’s body [
17,
18].
However, the question is not how to collect the data, but rather how to process them. In particular, in order to detect early signs of diseases based on cattle behavior, it is necessary to know the animals’ normal state and to predict and assess deviations from this norm. The development of algorithms for the interpretation of cattle behavior may help significantly [
19,
20,
21].
Berckmans (2013) highlighted that in the algorithm development process, different types of variables can be used and need to be defined in advance [
22].
The level of feeding is described as one of the most important indicators of animals’ health and well-being. It can be influenced by environmental factors such as food components, ambient temperature, conditions of keeping, and even group size. The level of water consumption is closely related to the feeding process [
23,
24,
25]; therefore, similar factors can also affect this indicator.
In their study, Das et al. (2016) revealed that any prolonged period of high temperature combined with high relative humidity of the air negatively affect food intake (decreased to 40%), milk production (decreased to 50%), and reproductive efficiency (decreased to 60%); at the same time, dramatically, the combination of these factors increases the probability of contracting parasitic diseases (by 43.3%) [
26].
Sejian et al. (2013) found that during heat stresses, the number of the white blood cells in cattle increased by 21–26%, and the number of the red blood cells decreased by 12–20%. This negatively affected the animals’ health and caused various diseases [
27].
This study aimed to develop algorithms for the detection of possible/incipient diseases in cattle; these algorithms are intended for digital systems for the internal monitoring of cattle physiological state.
2. Materials and Methods
The study was conducted by analyzing the available statistical data obtained in the course of practical studies conducted in real conditions in 5 large farms: “Ecocorm” Ltd. (Voronezh region, Russia), “Research Center of Traditional Technologies” Ltd. (Ivanovo region, Russia), Agricultural enterprise “Novomarkovskoye” Ltd. (Voronezh region, Russia), Agricultural enterprise “Voshchazhnikovo” Ltd. (Yaroslavl region, Russia). Thirty-seven heads of cattle were examined for 30 days with a measurement interval of 60 s. The cows’ selection for the examination was carried out in accordance with the following criteria: all cows were of the same breed (specifically, the Holstein breed), they were kept in the same conditions, and received the same unified diet. They were aged 2 or 3 years; all of them had their lactation period of 3 to 5 months, and their weight was at least 350 kg. Their average milk yield at the peak lactation was 30 (±3) kg/day. Their daily diet consisted of alfalfa hay (2.5 kg), corn silage (16.0 kg), alfalfa haylage (12.0 kg), propylene glycol (0.2 kg), glycerin (0.15 kg), and compound feeds (8.78 kg).
Prior to the start of the study, the selected cows were examined in accordance with the common clinical examination procedure [
5]. The clinical study used an algorithm involving the following sequential steps: anamnesis of the studied animals, their identification, habitus determination, examination of wool and skin, examination of visible mucous membranes, measurement of body temperature, palpation, auscultation, percussion and live weight measurement.
The average body weight of the studied cows was 517 (±2.03) kg. No clinical signs of disease were found in the first group.
For inclusion into the second group, 23 cows were selected with parameters similar to those of the cows in the first group but with the presence of some clinical signs of diseases such as encephalomyelitis, infectious enteritis, and hypothermia. In animals belonging in the second group, 3 groups of diseases were observed, namely, nervous, infectious, and parasitic diseases. The diseases of the nervous system include encephalitis, encephalomyelitis, and meningitis. Among the infectious diseases, the following ones were considered: brucellosis, infectious nodular dermatitis, and tuberculosis. As parasitic diseases, we evaluated hypodermatosis, cysticercosis and dictyocaulosis.
The input data were the measurements readings of the indicators of cattle health (pH level, rumen temperature, and locomotive activity). This preliminary data were obtained using the internal monitoring system prior to the main data collection. The measurements were carried out using sensors entered orally. Besides, microclimate parameters were measured, such as relative humidity and air temperature in the cowshed. Data monitoring was conducted via a radio transmitter on the ISM band (433 MHz). The combined database of cattle locomotive activity, pH level, and rumen temperature was transmitted via the basic station and then via a GSM module or Ethernet cable to the server.
The microclimate in the premises, where healthy cows of the first group were kept, was within the normal range (humidity 70–85%, temperature 10–20 °C). Deviations were observed within acceptable limits of the values (humidity 5–7%, temperature 3–10 °C). Meanwhile, in the parameters of the microclimate in the rooms where the sick cows were kept, deviations from the norms were periodically observed (humidity below 50%, tem-perature above 25 °C).
The measurement accuracy of the sensors was checked with the use of samples of the rumen content taken from fistulated cows. These samples were examined in our laboratory by veterinarians with the use of the standard stationary measuring devices. It was found that the deviations in the measured parameters were statistically insignificant (not exceeding 0.05).
The values of the rumen temperature, motor activity, and rumen pH for the cows in the first group were taken as the standardized boundary values of the ranged parameters; that is, when a cow presented values of the measured parameters similar to those of the group 1, the health condition of this cow was estimated as normal (as the setpoint). In contrast to the norm, The measured parameters’ values of the cows in the second group were considered as deviations from the setpoints, i.e., as deviations from the threshold values, which requested a predetermined action.
Based on the values obtained from sick cows, the main criteria were identified for 3 types of diseases (nervous, inflammatory, and parasitic).
In order to develop the algorithms, at first, we developed a program code in the software environment “Matlab R2019b” (Matrix Laboratory, the company “The MathWorks”, Natick, MA, USA). The available data were entered into this software environment, and then correlations between the measured indicators were calculated and graphs were obtained. Based on these results and taking into account the information from the analyzed scientific literature, we developed algorithms for the detection of infectious, parasitic, and nervous diseases.
3. Results
The analysis of the data obtained from the cows of the first group revealed the following facts:
- -
when cattle graze on a pasture, their motor activity is characterized by a very large amplitude of acts of movement (at least 20 units in 10 min);
- -
the duration of the increased activity is quite long (at least 2 h per day);
- -
this pattern is repeated regularly within a long period of time, more than two weeks (
Figure 1a);
- -
when walking on a pasture, cows can increase the temperature of their bodies (
Figure 1b).
The modern devices offered on the market are oriented in such a way that the y-axis represents a cow’s vertical movements (of head and legs), the x-axis is parallel to the cow body (i.e., to the direction of its walking), and the z-axis is located sideways to the cow (i.e., the z-axis reflects a cow’s movements to the left and to the right) (
Figure 2a).
When considering a healthy cows’ movements, in this case from the point of view of physical processes, it is proposed to consider the leg movements as a sinusoid of a pendulum oscillation. Thus, the physical processes of the oscillation are evaluated with the purpose to perform a diagnosis based on the cow’s motor activity (
Figure 2b).
As for walking, cows in the second group showed typical deviations corresponding to fluctuations in the graph of cows’ activity (the sinusoid was either ahead of the initial state or behind it); in contrast, cows in the first group presented a normal motor activity. A precise evaluation of this sign will allow the diagnosis of any deviation from the normal motor activity.
Figure 3 shows the graphs of the normal activity of the animals of the first group: the graphs are uniform and stable, and the activity index, on average, is not lower than 10.
In contrast,
Figure 4 shows the activity graph of one animal of the second group suffering from infectious enteritis (the activity index is below 5).
The graph in
Figure 4 regards one cow and is shown as an example; for other cows in the second group, the activity graphs are similar.
Thus, diseases contribute to decreasing animals’ motor activity level.
In addition to tracking the primary reaction of an animal when exposed to a particular factor, we also recorded the number of animals’ rest periods within a day, the average duration of one rest period, and the total duration of the resting time. The results are presented in
Table 1.
As can be seen from the data represented in the
Table 1, the shortest duration of the motor activity was observed for sick animals of the second group (801 min/day), which was on average 47 min longer than that for healthy animals. The differences between the groups are significant (
p < 0.05).
In animals belonging to the second group, three groups of diseases were observed, namely, nervous, infectious, and parasitic diseases. The diseases of the nervous system included encephalitis, encephalomyelitis, and meningitis. Among the infectious diseases, we found brucellosis, infectious nodular dermatitis, and tuberculosis. The parasitic diseases comprised hypodermatosis, cysticercosis, and distioculosis.
The input data for the algorithms were the values of the analyzed parameters of rumen temperature, rumen pH, locomotive activity, relative humidity and air temperature in the cowsheds (
Figure 1).
After data array processing by the program code, the program plotted the cross-correlation functions and also output a numeric array (
Figure 5).
In order to identify functions, dependencies, and correlations, statistical processing of the data array was carried out by the method of random variables. It was performed using the Formulas (1) and (2)
where
N is the measurements number,
Xi and
Yi are results of the
ith measurement,
and
are average values.
The normalized cross-correlation function of the random variables
X and
Y is an analogue of the linear correlation coefficient, with the only difference that one of the arrays is shifted several positions to the right or to the left:
where
is the cross-correlation function, and
j is delay (in measurements number).
From the analysis of the cross-correlation function on the Cheddock scale, we found that when measuring both the rumen temperature and the locomotive activity level at any given instant of time, a moderate negative correlation was observed, i.e., with a statistically significant probability, a decrease in one of the values meant an increase in the other one (
Figure 6a).
For measurements of rumen temperature and rumen pH, a moderate positive correlation between the values was observed, i.e., an increase in rumen temperature was more often associated with an increase in the pH than with a decrease; however, the difference in the frequency was rather insignificant (
Figure 6b).
For measurements of rumen temperature and air temperature in cowsheds, a remarkable negative correlation was observed, i.e., in an hour, any change in air temperature in the cowsheds had an impact on the rumen temperature (
Figure 6c).
For measurements of rumen pH and air temperature in the cowsheds, regardless of how they were made, i.e., simultaneously or in conditions of various delays, a moderate negative correlation was observed, i.e., the air temperature affected the rumen pH value within an hour after the measurements were taken (
Figure 6d).
The average pH of the rumen in the group under consideration was 6.08; the standard deviation (±) was 0.29. The average rumen temperature in the group under consideration was 38.92 (±0.19). The average level of the locomotive activity in the group under consideration was 8.6 (±1.4).
The highest modulus average value of the linear correlation coefficients was −0.222 (±0.117); it was reached between the locomotive activity levels and the rumen pH values.
For some animals, the linear correlation coefficient reached −0.519; this correlation was observed between pH level and locomotive activity; hence, this was a noticeable correlation and suggested a stable connection. The correlation between the parameters of rumen temperature and locomotive activity level in some cows reached −0.424, and the same was observed for the correlation between rumen temperature and pH level that reached −0.470, which is also close to what can be called a moderate correlation and can suggest a stable mutual connection.
As a result of the analysis of the mutual correlation functions, the following values were obtained: between rumen temperature and locomotive activity, a correlation value of 0.33; between temperature and rumen pH, 0.29; between rumen pH and air temperature, 0.37; between rumen pH and locomotive activity, 0.16; between rumen pH and relative humidity, 0.29; between locomotive activity level and relative humidity, 0.37; between locomotive activity level and air temperature, 0.62; between rumen temperature and relative humidity, 0.18.
It should also be noted that we found a statistically insignificant correlation between the parameters of rumen temperature and ambient air temperature (−0.05) in the absence of abnormal ambient temperature conditions.
Based on the obtained criteria and graphs of the cross-correlation functions, algorithms were developed for the detection of infectious, nervous, and parasitic diseases (
Figure 6a–d). As shown in
Figure 6, we took into account both the main health criteria, namely, rumen temperature and physical activity level, and additional ones such as rumen pH and indoor air parameters (temperature and humidity).
When issuing a message about the presence of a disease in a particular group of cows, the algorithm provides general recommendations on how to minimize the further development of the clinical signs.
To improve the detection, some additional conditions were checked, whose fulfillment increased the probability of each disease being more severe (
Figure 6d).
It should be noted that, first of all, the mathematical model (written based on the obtained algorithms) calculated the animals’ location (cattle house or pasture). This was done because the animals’ motor activity differs depending on the conditions in which they are kept (loose or stable-walking conditions (when there is some motion)); however, this difference is unrelated to any disease. In connection with this, for the two locations, the setpoint (and, as a result, the standard deviation and regression) was different. The animal location determination was done based on the analysis of motor activity amplitude. However, due to the difference in the individual characteristics of the animals, the amplitudes of their motor activity differed significantly. That is why the data on the animals’ location should be considered as approximate. In this regard, in order to obtain more accurate information, it is recommended to consider the possibility of registering the grazing time in a journal and subsequently transfer this information to the owner. In this case, it would be possible to improve and finalize the computer program.
4. Discussion
The existing algorithms for detecting deviations in indicators of the physiological state of cattle do not consider the location of the animals under study (cows are in the pasture or in the barn). That is why we conducted a comparative analysis of the motor activity index of both sick and healthy dairy cows. The cows of the second group presented deviations in frequency and amplitude of their legs’ oscillatory (pendulum) movements when walking (the sinusoid was ahead of the normal sinusoid or lagged behind it). Meanwhile, all parameters were within the normal range for the cows of the first group.
The animals’ activity index for the first group was on average not less than 10 units, while for the second group, it was below 5 units. The lowest duration of the motor activity was observed in sick animals of the second group: their rest time was 801 min/day, which was longer on average by 47 min (by 6%) than that of heathy cows. The differences between the groups are statistically significant (p < 0.05).
Davison et al. (2020) came to the conclusion that the consideration of the locations was necessary to identify deviations from the physiological state in in milk cows. It is important because, for example, heat stress significantly affects the body temperature and the locomotive activity of animals; therefore, not always deviations from the norm would be associated with diseases [
28]. Ramón-Moragues et al. (2021) revealed that under heat stress, cows noticeably decreased their locomotive activity and body temperature; besides, there was a violation of the animals’ daily routine [
29]. Rutherford et al. (2019) came to the conclusion that the rumen temperature varies greatly depending on the animals’ location. For example, when being in a pasture, due to the predominance of succulent feed in the diet, the rumen pH was below normal, while its temperature was higher than the norm [
30]. The indoor microclimate has a connection with animals’ measured values. For example, a temperature above +20 °C causes heat stress: energy consumption increases for the maintenance of body optimal temperature, as a result of which the costs of other processes decrease, and immune activity decreases [
31]. In high humidity, it is more difficult for an animal’s body to cope with infectious diseases; animals’ metabolism slows down, their blood flow becomes disordered, and the level of hemoglobin in their blood decreases [
32]. For this reason, when creating algorithms, it is necessary to take into account the animals’ location.
As for the identification of correlations between the physiological state’s parameters of cattle, we confirmed successfully the results obtained in the previous studies.
Additionally, we obtained the following data: the average value of the linear correlation coefficients for rumen temperature and locomotive activity was 0.199 (=0.088), and that between level of locomotive activity and rumen pH was −0.222 (=0.117). Between pH level and level of locomotive activity of the animals, the linear correlation coefficient reached −0.519, which suggests the presence of a stable connection between these parameters. The correlation between rumen temperature and locomotive activity level reached −0.424, and that between rumen temperature and rumen pH level was −0.470, which also suggests a stable connection. The analysis of the cross-correlation function showed a noticeable connection between level of locomotive activity and air temperature (0.622); moderate correlations were observed between rumen temperature and locomotive activity (0.33), between rumen pH and air temperature (−0.36), and between level of locomotive activity and relative humidity (−0.347).
Similar results regarding the high correlation between the studied cattle’s physiological parameters were obtained by Alzahal et al. (2008) and by Antanaitis et al. (2016). The research results by Dorokhov et al. (2021) showed a negative statistically significant correlation (r = −0.932,
p < 0.01) between rumen temperature and its pH [
33,
34,
35].
Soriani et al. revealed a significant correlation between rectal temperature and the temperature–humidity index (THI) (r = −0.41;
p < 0.01). In this regard, it is very important to take the microclimate into account, when developing algorithms for disease detection [
36].
Based on the analysis of the scientific literature, three main groups of cattle diseases were identified: infectious, invasive (parasitic), and nervous (nervous system diseases) (Bessarabov et al. 2007; Isimbekov et al. 2016). These authors described the main criteria to diagnose these diseases (
Figure 6) and how they can be detected at the initial stages of their development before the appearance of the main clinical signs; they also provided recommendations for their treatment [
37,
38].
As for the existing monitoring systems, they do not allow detecting diseases in their early stages because the algorithms used are outdated. Wang et al. (2018) highlighted, for example, that it is difficult to assess the level of feeding by using only an accelerometer; however, the combination of these data with those from other sensors would improve the assessment by as much as 20%. This principle can be applied to the rest of the measured physiological parameters of cattle [
39].