Next Article in Journal
An Explorative Study of the Causal Pathogenesis of Green Liver Discoloration in Organically Reared Female Bronze Turkeys (Meleagris gallopavo) Considering the Infectious Risk Factors
Next Article in Special Issue
Regrouping in Dairy Ewes—Effects on Productive Performance and Specific Behavioral Traits
Previous Article in Journal
The Gut Microbiota of Young Asian Elephants with Different Milk-Containing Diets
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

The Effectiveness of Virtual Fencing of Bull Calves in a Holistic Grazing System

by
Søren Krabbe Staahltoft
1,*,†,
Magnus Fjord Aaser
1,†,
Jakob Nødgaard Strange Jensen
1,†,
Ismat Zadran
1,†,
Emil Birkmose Sørensen
1,†,
Anders Esbjerg Nielsen
1,†,
Aage Kristian Olsen Alstrup
2,3,
Dan Bruhn
1,4,
Anne Cathrine Linder
5,
Christian Sonne
6,
John Frikke
7 and
Cino Pertoldi
1,8
1
Department of Chemistry and Bioscience, Aalborg University, Fredrik Bajers Vej 7H, 9220 Aalborg, Denmark
2
Department of Nuclear Medicine and PET, Aarhus University Hospital, Palle Juul-Jensens Boulevard 165, 8200 Aarhus, Denmark
3
Department of Clinical Medicine, Aarhus University, Palle Juul-Jensens Boulevard 165, 8200 Aarhus, Denmark
4
Skagen Bird Observatory, Fyrvej 36, 9990 Skagen, Denmark
5
National Institute of Aquatic Resources, Technical University of Denmark, Kemitorvet 201, 2800 Kongens Lyngby, Denmark
6
Department of Ecoscience, Aarhus University, Frederiksborgvej 399, 4000 Roskilde, Denmark
7
Wadden Sea National Park, Havnebyvej 30, 6792 Rømø, Denmark
8
Aalborg Zoo, Mølleparkvej 63, 9000 Aalborg, Denmark
*
Author to whom correspondence should be addressed.
These authors contributed equally to this work.
Submission received: 13 February 2023 / Revised: 24 February 2023 / Accepted: 1 March 2023 / Published: 2 March 2023
(This article belongs to the Special Issue Ruminant Welfare Assessment)

Abstract

:

Simple Summary

Virtual fencing is a new way of enclosing livestock without the use of physical barriers. The system relies on GPS technology and works by deterring escapes using auditory warnings and electric impulses. This study examines the effectiveness of one such system in a rotational strip grazing regime with 17 Angus bull calves. This study also aimed to determine which bull calves were more likely to receive electric impulses. The system proved effective at containing the bull calves within the designated enclosure. The bull calves learned to associate the auditory warning with the electric impulse, and as such, received increasingly fewer impulses as time went on. However, in accordance with previous studies, there were clear differences in the number of warnings and impulses received between individuals. The results regarding which bull calves received the most impulses did not show any clear trends. In conclusion, the virtual fencing system proved to be effective and the animals learnt how to interact with the system while receiving very few impulses.

Abstract

Large grazers are essential for nature conservation. In order to prevent grazers from moving to unintended areas, it may be necessary to keep them inside enclosures. Physical fences present a number of problems, such as fragmenting the landscape. Virtual fencing, however, is a possible replacement for physical fencing, making it possible to enclose grazers without physical boundaries. Virtual fencing systems utilise collars with GPS technology to track animals and deliver auditory warnings and electric impulses to keep animals within predefined boundaries. This study examines how effective the virtual fencing system Nofence© is at enclosing calves in a holistically managed setting. Holistic management is a rotational grazing technique where an enclosure is grazed in small strips at a time. It is investigated whether the calves become habituated to the virtual fence and whether there is a correlation between the number of warnings received by every two calves in order to explore potential herd behaviour. Finally, this study examines which calves interact the most with the virtual fence by investigating the relationship between physical activity and number of interactions. Seventeen calves were fitted with a GPS collar from the company Nofence© and placed in a holistically managed enclosure. Data were gathered from 4 July to 30 September 2022. The study found that virtual fence was able to contain calves inside the designated enclosure, and over time the calves received notably fewer electrical impulses compared to auditory warnings. The results of Pearson’s correlation between auditory warnings received by two random calves were inconclusive, but the use of a sliding window analysis should be further explored. Lastly, the most physically active animals were the ones who received the most auditory warnings, but they did not receive more impulses. No significant correlation was found between the number of electric impulses received and the physical activity of the animals.

1. Introduction

Conservation of natural areas using grazers is required to reduce the loss of biodiversity [1,2]. This often necessitates the use of physical fencing to keep grazers inside the desired area, as well as to prevent unwanted human encounters and traffic accidents [3]. However, physical fencing presents a number of problems [3,4]. Physical fences can thus have negative impacts on wildlife by creating artificial barriers in the landscape and limiting the free movement of other animals [3,5]. These barriers can also lead to collisions, such as when low-flying bird species collide with fences with lethal consequences [3,6,7,8,9]. An example of an affected species, which are comparable in size and ecology to native European red deer and fallow deer, are the North American white-tailed deer and mule deer, as described by Harrington et al., Burkholder et al. and Bishop et al. [7,10,11]. These papers found that deer were prone to getting caught and dying when attempting to cross fences, with fence crossing success rates at around 75% and mortality rates upwards of 0.40/km/year for pasture fences [7,11]. Moreover, physical fencing can be expensive to install and maintain [5]. Even when fencing is used for nature conservation purposes, it can create new problems when trying to solve others, as physical fences fragment the landscape by preventing animals from moving around freely [3,5,12]. Additionally, physical fencing poses a logistical and practical problem in terms of relocation [6]. In some areas that are difficult to access or have challenging terrain, installing physical fences may not be feasible [13].
Virtual fencing is a possible alternative to physical fencing, without many of the aforementioned problems of physical fencing [13,14,15,16,17,18]. Virtual fencing is a GPS-based system that utilises auditory warnings and electrical impulses to keep livestock within a designated area. As the animal approaches the virtual boundary, warning signals and electrical impulses are administered by a collar fitted around the neck of each animal. This type of system is yet to be legalised in most countries in Europe, including Denmark, due to earlier concerns regarding animal welfare [6,19,20]. However, several studies suggest that welfare is not compromised when utilising virtual fencing systems and livestock are effectively kept within enclosures [13,14,18,20,21,22,23,24,25,26].
In a previous study on virtual fencing by Aaser et al. [14], it was hypothesised that a small enclosure size and regular moving of the virtual border could improve the effectiveness of virtual fencing systems. The reasoning behind the hypothesis was the cattle’s constant interaction with the boundary, through which they would be reminded of the connection between warnings and impulses [14]. According to the responsible veterinarian, the preliminary results of an experiment with holistically managed cattle on the island of Fanø raised some concern, as the animals received a high number of warnings and impulses compared to animals kept in a larger stationary enclosure. Holistic management is a rotational grazing technique where an enclosure is grazed in small strips at a time. The strip is shifted a few metres at regular intervals, usually daily. To the best of our knowledge, little to no research has been conducted to examine the effectiveness of virtual fencing in a holistically managed setting.
Cattle are, in general, gregarious and social creatures living in herds [27,28,29]. They establish hierarchies, where some of the individuals lead the herd, while others are being lead [30]. As suggested by by Aaser et al. [14], one might hypothesise that putting a collar on only some of the animals may be sufficient to keep the entire herd within a virtual enclosure, as has been shown to be the case for sheep [14,31]. This could potentially help reduce the cost even further.
The aim of this study was to investigate whether the Nofence© virtual fencing system is able to contain bull calves in a holistic management setting, whether the calves display herd behaviour and whether it is possible to predict which animals will interact most with the virtual boundary. To support this aim, the following four hypotheses were tested:
(1) A virtual fence is capable of containing bull calves in a specific area during holistic management for more than 99% of the time. (2) Bull calves will receive fewer electrical impulses compared to warnings over time. (3) There is a positive correlation between the number of warnings received by any two bull calves in any given period. (4) Bull calves with the most physical activity receive the most warnings and impulses.

2. Materials and Method

2.1. Animals and Location

This study took place on the Danish island of Fanø. Seventeen Angus bull calves aged 9 to 12 months were placed in a holistically managed enclosure of 8.2 hectares (Figure 1 and Appendix A). Some of these calves had previously been indirectly exposed to virtual fencing as their mothers were kept in another virtually fenced enclosure, though the bull calves themselves did not wear collars. The pasture was composed of flat cultural grassland. Due to the close proximity of a somewhat busy road, an electric fence was placed around the perimeter of the pasture. An artificial shelter and a water trough were provided in the northwestern side of the enclosure. Furthermore, the neighbouring farmer to the east had a herd of cattle adjacent to the pasture.

2.2. Virtual Fencing System

Each animal was fitted with a cattle collar developed by Nofence©. At the time of the experiment, each collar cost NOK 3050, while the system, regardless of the number of collars, cost a yearly digital service fee of NOK 669 (prices provided by Nofence© on 5 December 2022). The collars were the same as described by Aaser et al. [14]. The collars were capable of registering five different types of messages with different information, which were sent to the official Nofence© server (Table 1).
The position of the animal was logged every 15 min and the activity every 30 min. Activity was measured by a step counter developed for another species (sheep) and given as a unitless number. Whenever an animal received a warning or an electric impulse, or the fence status of the collar changed, this was logged at the time of the event along with the position of the animal. The different fence statuses are outlined in Table 2. With every message, the unique serial number of the collar was included. This number was used to reference the animal to which the collar was fitted. Each animal wore the same collar throughout the experiment.

2.3. Experimental Protocol

The experiment began on 4 July 2022 and included seventeen Angus bull calves. Each calf was fitted with a collar and was let onto the western side of the field. During the initial days of the experiment, the calves were restricted to an inner enclosure of about 1 ha by an electric fence (Appendix B). This inner enclosure explains why the calves are located to the west during the beginning of the experiment (Appendix B). The virtual boundary was placed close to the electric fence, allowing the animals to receive auditory warnings when approaching the electric fence. The collars also administered electric impulses if the animals continued to move towards the electric fence after receiving a warning. On 9 July, one side of the internal fence was removed, meaning the calves would have to respect the virtual boundary to stay within the enclosure. On 10 July, the other side of the internal fence was removed, and on 11 July, all internal fences were removed. On 12 July, the enclosure size was doubled to 2 ha based on an assessment of the veterinarian associated with the project. At regular intervals, usually around every 24 h, the virtual enclosure was shifted 10–20 m. The fence was shifted eastwards in a circular motion, anchored at the water trough and resting shelters, until the physical boundary of the field was reached after 15 shifts and the procedure started over from the western boundary again (Appendix A). This method was used in order to maintain a balanced grazing intensity. Data collection began on 4 July 2022 and ended on 30 September 2022, totalling a period of 89 days.

2.4. Data and Statistical Analysis

All data were downloaded from the official Nofence© servers. Five different message types, sent by the collars, were given as point data, and the virtual enclosure history was given as a polygon for each enclosure version. The start time and end time of each version of the enclosure was recorded. Poll messages were sent every 15 min and recorded the position of the animal, the time and the battery voltage in centivolts. Seq messages were sent every 30 min and described the activity of the animal and the solar charge accumulated by the collar in the last 30 min. Warning and Zap messages were sent whenever a warning or an electric impulse had been delivered to the animal. These messages provided the time and position of the animal. Warning messages also provided the length of the warning in milliseconds. Status messages were sent whenever the fence status of the collar changed, such as when an animal had escaped the enclosure (see Table 2). A total of 241,163 data points were collected. A total of 151,576 were Poll messages, 71,990 were Seq messages, 13,906 were Warning messages, 2153 were Status messages and 1538 were Zap messages. The data were subdivided by message type and individual. QGIS [32] was used to plot and sort the data, while all statistical analyses were carried out in R version 4.0.3 [33]. Only non-parametric tests were used, as the data did not follow a normal distribution. The collar fell off once for three of the calves, all in the month of August. One was refitted on the same day, while the other two were refitted after two and six days, respectively.

2.4.1. Effectiveness of Virtual Fencing

To evaluate the effectiveness of the virtual fencing system, the number of Poll messages with the fence status Normal was compared to the number of Poll messages with fence status MaybeOutOfFence and fence status Escaped:
N o r m a l N o r m a l + M a y b e O u t O f F e n c e + E s c a p e d
Messages with the status MaybeOutOfFence were included to err on the side of caution with regard to the amount of time animals spent outside the enclosure. To visualise the movement of the animals, heatmaps were created for each version of the enclosure. Heatmaps were created for enclosures that were active for more than 12 h. This was done to exclude “intermediary” enclosures, where only one side of the enclosure had been shifted. The number of times an animal broke out of the enclosure was recorded, as well as the length of each breakout. The length of each breakout was determined as the duration of time between the Fence Status of the collar changing from Normal to Escaped and back to Normal again, rounded to the nearest half minute.

2.4.2. Habituation to Virtual Fencing

The animals’ habituation to the virtual fence was investigated by recording the number of auditory warnings and electric impulses received by each animal. The cumulative warnings and impulses were calculated and assessed, as well as the length of warnings over time. Furthermore, the number of electric impulses divided by the number of warnings was calculated over time, and the correlation between warnings and electric impulses was calculated using Pearson’s product–moment correlation ( r p ).

2.4.3. Herd Behaviour

The presence of herd behaviour was investigated by calculating the Pearson’s correlation ( r p ) between the number of warnings received by a given calf and any other calf in an 8 h period (00–08 h, 08–16 h and 16–00 h). This was calculated for every pairwise comparison of the 17 animals, resulting in 136 correlations. A sliding window analysis was performed by repeating this calculation for shorter intervals. The correlation was calculated in an 18 period window, denoting the amount of consecutive 8 h periods included in each sliding window position.

2.4.4. Sliding Window Correlation

To account for possible temporal variations in correlations between the animals, a sliding window correlation was used. A sliding window correlation is a way to assess a time series of data in smaller overlapping segments at a time, thus allowing for analysis of a long time series while still capturing short-range temporal variations in the correlation. In practice, the time series is tested one “window” at a time with a certain window size. The window is then shifted one position and the test is run again. As an example, if a window size of four is selected, the test is first run on data from day 1–4, then day 2–5, day 3–6 and so on. The results of these tests are recorded at the first point in the window, denoted the “sliding window position”. The size of the window is increased until the median p-value of all correlations is below 0.05, and this is the included result. This way, more than half of the correlations were significant.

2.4.5. Predicting the Individuals with the Most Interactions with the Fence

To determine whether the animals with the most physical activity received the most warnings, Spearman’s rank correlation ( r s ) was calculated between the summed activity of an animal for the entire period and the summed number of warnings and impulses received by the same animal. The Spearman’s correlation ( r s ) was also calculated for the summed number of warnings against the summed number of impulses the animals had received. Furthermore, the Spearman’s correlation ( r s ) was calculated between the summed number of warnings an individual had received and the same individuals mean correlation to all the other animals, with regard to the number of warnings received in 8-h periods.

3. Results

3.1. Effectiveness of Virtual Fencing

The virtual fencing system effectively kept the animals within the periodically moving enclosure during the 89-day period of data collection. In total, the animals spent 99.74% of the period inside the enclosure. The created heatmaps also indicated that the animals respected the virtual boundaries and stayed within the enclosure (Figure 2 and Appendix B).
Over the course of the 89-day period, 19 breakout events were recorded. If multiple animals escaped at the same time, this was counted as one event per escaped animal. The median escape time was 2 min and the most number of breakouts by an animal was three, with 12 of the 17 calves breaking out of the enclosure at least once (Figure 3).

3.2. Habituation to Virtual Fencing

There was a noticeable variation in the number of warnings received by the animals. The lowest total number of warnings received by a calf was 237 (an average of 2.66 warnings per day), while the highest number received was 1874 (an average of 21.06 warnings per day). The frequency of warnings received by the animals seemed to decrease from late August (Figure 4).
As with the number of warnings, the number of impulses received also varied between individuals. The lowest total number of impulses received by a single calf was 36 (an average of 0.40 impulses per day), while the highest number received was 181 (an average of 2.03 impulses per day). It is worth noting that for all calves, the frequency of impulses was highest at the start of the period and stagnated towards the end of the data collection period (Figure 5).
During the the start of the experiment, the warning duration ranged from approximately 0.01 s to 65 s, with a wide spread and many outliers. Within the first couple of weeks, the spread reduced, as well as the median warning duration, which settled around 8–10 s (Figure 6).
The median number of impulses an animal received per warning in the beginning was approximately 0.5, but immediately trended downwards before reaching a median of 0 in late August (Figure 7).
The correlation between the number of daily warnings and impulses an animal received throughout the period showed considerable variation, with the weakest correlation being r p = 0.43 (ID: 90754) and the strongest being r p = 0.85 (ID: 89495).

3.3. Herd Behaviour

The results of the correlation and sliding window correlation tests run in regards to analysing herd behaviour were inconclusive and have therefore been moved to Appendix C.

3.4. Predicting the Number of Warnings for Each Animal

A moderate positive correlation between the summed activity and summed number of warnings an animal received throughout the data collection period was observed ( r s = 0.56 , p < 0.05 ). No correlation between summed activity and summed impulses was found. Equally, no significant correlation was found between the summed number of warnings a bull calf had received and the mean correlation between the number of warnings received by the same calf and all other calves.

4. Discussion

4.1. Effectiveness of Virtual Fencing

Throughout the 89-day period, the virtual fence successfully kept the 17 calves inside the virtual enclosure zone 99.74% of the time, even with daily changes in the virtual boundary. This supports Hypothesis (1), which states that the system will contain the animals more than 99% of the time. This is comparable to another study, where 12 dairy cows were kept within a 0.83 hectare enclosure for 99% of the time in a six day period [25]. The positional data of the calves clearly showed that they respected the virtual boundary (Appendix B). The calves spent most of their time near the sheltering facilities and the water trough located in the northwest of the enclosure. Moreover, the calves grazed close to the virtual boundary to the east, where fresh grass was available (Appendix B). The fact that the calves clearly respected the virtual boundary despite the higher feed availability outside the boundary corresponds with results from other studies where virtual fencing kept cattle away from feed attractants [22,26]. A total of 19 breakouts were recorded between 12 of the animals; however, the fact that the median time before an escaped animal returned to the enclosure was just 2 min could indicate that the breakouts were not deliberate attempts to escape (Figure 3). The remaining five animals did not escape at all. According to the farmer, the calves would sometimes push each other outside the virtual boundary, which could explain some of the breakouts. The ease with which the animals were able to re-enter the enclosure after a breakout, due to the lack of a physical barrier, meant that the farmer did not have to shoo the escaped animals back into the enclosure. The gregarious nature of the bull calves could also explain the short time between an escape and the animal re-entering the enclosure, as they would seek to return to the herd [34]. The five animals that did not escape the enclosure throughout the entire experiment were also the animals that had received the fewest number number of impulses (Figure 3 and Figure 5). This could be an indication that these animals were less prone to be pushed around by their herdmates, or were simply less inquisitive regarding the enclosure border.

4.2. Habituation to Virtual Fencing

The frequency of warnings and impulses received by the animals reduced over time, with the number being notably lower in September than in July (Figure 4 and Figure 5). This indicates that the bull calves needed fewer warnings to learn where the virtual boundary was placed, despite the fact that it was shifted at regular intervals. However, in line with several other studies, the amount of warnings and impulses received varied considerably between individuals [14,22,25,26,35]. Equally, the Pearson’s correlation between the number of daily warnings and impulses received varied between individuals (Section 3.2). This variation could be an expression of different personalities among the bull calves or the hierarchical structure of the herd. Some calves continued to receive a relatively high number of warnings, albeit at a lower frequency, while others barely received any warnings in the latter half of the data collection period. This could indicate that those receiving many warnings are “testers” that determine where the new virtual boundary is placed for the rest of the herd. However, this hypothesis is not immediately supported by the observation that, throughout the data collection period, the correlation between the number of warnings received by any two calves was positive for all pair-wise comparisons (Figure A19). If negative correlations had been observed, this would clearly indicate that when these “testers” received warnings, the rest of the herd would actively avoid warnings. However, even if some animals tried to avoid warnings because they observed others interact with the virtual boundary, this would likely be masked by the gregarious nature of cattle. As such, even if they actively tried to avoid warnings by observing others, they would still receive more warnings whenever the herd moved closer to the virtual boundary, making positive correlations very likely, especially due to the narrowness of the virtual enclosure. Despite the inter-individual differences, it is clear that the bull calves became habituated to the virtual fencing system. As such, the amount of impulses compared to warnings decreased over time, which supports Hypothesis (2) (Figure 7). This hypothesis states that the animals would receive fewer impulses per warning over time. This indicates that the animals were increasingly able to react to the auditory warnings to avoid receiving an electric impulse. The median amount of impulses per warning even reached zero in late August. This is positive, as it clearly shows that the bull calves learnt to associate the warnings with impulses, and reacted favourably when receiving warnings. Another study indicated that increased stress does not occur among cows managed with virtual fencing [21]. The decreasing trend in impulses per warning is in line with other studies, which found that the number of warnings and total impulses decreased significantly over time [16,36]. The decrease in the median length of warnings supports the notion that the calves became habituated to the system and reacted to the warnings, turning away from the virtual boundary before receiving an impulse (Figure 6). The warning duration stabilised after a few weeks, which could mean the animals became so familiar with the system that they learnt how long they could safely wait before turning away from the virtual boundary without receiving an impulse.

4.3. Herd Behaviour

As a new way of analysing temporal data gathered from animals, such as from these calves, a sliding window correlation was run on the number of warnings any two calves had received in subsequent 8 h periods. A sliding window correlation makes it possible to capture any variations or trends in time, rather than relying on a single result independent of time. In this study, it was decided that due to the large window size of 8 h, the results would be too inaccurate and any discussion of these too speculative. Therefore, the results are only included in Appendix C. The authors of this paper suggest that future studies should continue to explore the use of sliding window correlations and other sliding window analyses to better capture temporal changes and trends.

4.4. Predicting the Number of Warnings for Each Animal

A moderate correlation between summed activity and summed number of warnings received was found (Section 3.4). This supports Hypothesis (4), that animals with higher activity receive more warnings. Animals that are more active tend to encounter the virtual fence more frequently than those that are more sedentary. This was expected, as an animal has to be moving to receive a warning. The fact that the correlation was only moderate may be explained by other confounders, such as risk tolerance, foraging and herd behaviour. However, no correlation was found between summed activity and summed number of impulses, indicating that overall activity is not a contributing factor in receiving more impulses (Section 3.4). This does not support our Hypothesis (4), that the most active calves receive the most impulses. Interestingly, the correlation between summed warning and activity and the lack of correlation between impulses and activity indicate that animals that move around more do not receive more impulses than more sedentary animals, despite receiving more warnings. It is possible that the more active animals are better at reacting to the warnings, as they receive them more frequently. It could also simply be that all the animals had learnt to avoid impulses, such that only the number of warnings was affected by the activity of the animal. Further research is required to ascertain which statement is most likely. Lastly, a moderate correlation between summed number of warnings and summed number of impulses was found, which was to be expected, as an animal cannot receive an impulse without first receiving a warning. The correlation is not strong because the animals can receive warnings without impulses. The correlation does support the notion that the animals are learning to react to the warnings when the lack of correlation between activity and impulses is taken into account. A correlation between warnings and impulses would indicate that more active animals receive more impulses due to receiving more warnings, but, as discussed, that is not the case.

5. Conclusions

In conclusion, the Nofence© virtual fencing system was capable of keeping the bull calves in a holistically managed enclosure for more than 99% of the time, and the calves that escaped returned of their own accord after a short time. The bull calves became habituated to the virtual fencing system, receiving fewer warnings and impulses over time. Additionally, they learnt to respond favourably to auditory warnings, receiving fewer impulses compared to warnings over time, reaching a median of zero impulses per warning for the last month of the data collection period. The median warning length decreased at the beginning of the period but then stabilised, indicating that the calves became so used to the system that they did not react to warnings immediately, but had figured out when to turn away from the boundary to avoid impulses. A sliding window correlation was used to test for short range changes or temporal trends in the median correlation between the number of warnings received by any two calves, but the results were inconclusive. Some short-range fluctuations were shown using the analysis, showcasing a need for future studies to take temporal variations into account when analysing these types of data. The calves with the highest physical activity received the most warnings, but no other significant results were found when trying to predict which animals would receive the most warnings.

Author Contributions

Conceptualisation, S.K.S., M.F.A., J.N.S.J., I.Z., E.B.S., A.E.N., A.K.O.A., D.B., A.C.L., C.S., J.F. and C.P.; methodology, S.K.S., M.F.A., J.N.S.J., I.Z., E.B.S., A.E.N., A.K.O.A., D.B. and C.P.; formal analysis, S.K.S., M.F.A., J.N.S.J., I.Z., E.B.S. and A.E.N.; investigation, S.K.S., M.F.A., J.N.S.J., I.Z., E.B.S. and A.E.N.; data curation, S.K.S., M.F.A., J.N.S.J., I.Z., E.B.S. and A.E.N.; writing—original draft preparation, S.K.S., M.F.A., J.N.S.J., I.Z., E.B.S. and A.E.N.; writing—review and editing, S.K.S., M.F.A., J.N.S.J., I.Z., E.B.S., A.E.N., A.K.O.A., D.B., A.C.L., C.S., J.F. and C.P.; visualisation, S.K.S., M.F.A., J.N.S.J., I.Z., E.B.S. and A.E.N.; supervision, A.K.O.A., D.B., A.C.L., C.S., J.F. and C.P. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by 15. Juni Fonden, Hedeselskabet and Markus Jebsens Naturpulje. We are grateful for the support, which made it possible to conduct this study.

Institutional Review Board Statement

The study was conducted according to the institutional guidelines for animal research (directive 2010/63/EU), and approved by the Animal Experiments Inspectorate of Denmark (2020-15-0201-00588, approval date: 8 July 2020) prior to the start of the study.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data presented in this study are available on request from the corresponding author.

Acknowledgments

Special thanks to Dan Pode Poulsen and Michael Baun for providing data, consulting, continuously sharing information and for letting us carry out the experiment as part of the ongoing project of testing Nofence© virtual fencing on the Danish island of Fanø. This project was conducted by a partnership consisting of landowner Dan Pode Poulsen, farmer Michael Baun, The Wadden Sea National Park, The Municipality of Fanø and Hedeselskabet.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A. Maps of Virtual Enclosures

Figure A1. Virtual enclosure version 1 dating from 4 July to 10 July. The dashed black line represents the physical fence surrounding the 8.2 ha grazing area. The marked green area represents the virtual enclosure. Only virtual enclosures active for more than 12 h were included. The coordinates of the upper left corner of the map frame are 55°25′23.131″ E 8°23′53.531″ N. Map data from OpenStreetMaps http://openstreetmap.org/copyright.
Figure A1. Virtual enclosure version 1 dating from 4 July to 10 July. The dashed black line represents the physical fence surrounding the 8.2 ha grazing area. The marked green area represents the virtual enclosure. Only virtual enclosures active for more than 12 h were included. The coordinates of the upper left corner of the map frame are 55°25′23.131″ E 8°23′53.531″ N. Map data from OpenStreetMaps http://openstreetmap.org/copyright.
Animals 13 00917 g0a1
Figure A2. Virtual enclosure version 2 dating from 10 July to 18 July. The dashed black line represents the physical fence surrounding the 8.2 ha grazing area. The marked green area represents the virtual enclosure. Only virtual enclosures active for more than 12 h were included. The coordinates of the upper left corner of the map frame are 55°25′23.131″ E 8°23′53.531″ N. Map data from OpenStreetMaps http://openstreetmap.org/copyright.
Figure A2. Virtual enclosure version 2 dating from 10 July to 18 July. The dashed black line represents the physical fence surrounding the 8.2 ha grazing area. The marked green area represents the virtual enclosure. Only virtual enclosures active for more than 12 h were included. The coordinates of the upper left corner of the map frame are 55°25′23.131″ E 8°23′53.531″ N. Map data from OpenStreetMaps http://openstreetmap.org/copyright.
Animals 13 00917 g0a2
Figure A3. Virtual enclosure version 3 dating from 18 July to 26 July. The dashed black line represents the physical fence surrounding the 8.2 ha grazing area. The marked green area represents the virtual enclosure. Only virtual enclosures active for more than 12 h were included. The coordinates of the upper left corner of the map frame are 55°25′23.131″ E 8°23′53.531″ N. Map data from OpenStreetMaps http://openstreetmap.org/copyright.
Figure A3. Virtual enclosure version 3 dating from 18 July to 26 July. The dashed black line represents the physical fence surrounding the 8.2 ha grazing area. The marked green area represents the virtual enclosure. Only virtual enclosures active for more than 12 h were included. The coordinates of the upper left corner of the map frame are 55°25′23.131″ E 8°23′53.531″ N. Map data from OpenStreetMaps http://openstreetmap.org/copyright.
Animals 13 00917 g0a3
Figure A4. Virtual enclosure version 4 dating from 26 July to 3 August. The dashed black line represents the physical fence surrounding the 8.2 ha grazing area. The marked green area represents the virtual enclosure. Only virtual enclosures active for more than 12 h were included. The coordinates of the upper left corner of the map frame are 55°25′23.131″ E 8°′53.531″ N. Map data from OpenStreetMaps http://openstreetmap.org/copyright.
Figure A4. Virtual enclosure version 4 dating from 26 July to 3 August. The dashed black line represents the physical fence surrounding the 8.2 ha grazing area. The marked green area represents the virtual enclosure. Only virtual enclosures active for more than 12 h were included. The coordinates of the upper left corner of the map frame are 55°25′23.131″ E 8°′53.531″ N. Map data from OpenStreetMaps http://openstreetmap.org/copyright.
Animals 13 00917 g0a4
Figure A5. Virtual enclosure version 5 dating from 4 August to 11 August. The dashed black line represents the physical fence surrounding the 8.2 ha grazing area. The marked green area represents the virtual enclosure. Only virtual enclosures active for more than 12 h were included. The coordinates of the upper left corner of the map frame are 55°25′23.131″ E 8°23′53.531″ N. Map data from OpenStreetMaps http://openstreetmap.org/copyright.
Figure A5. Virtual enclosure version 5 dating from 4 August to 11 August. The dashed black line represents the physical fence surrounding the 8.2 ha grazing area. The marked green area represents the virtual enclosure. Only virtual enclosures active for more than 12 h were included. The coordinates of the upper left corner of the map frame are 55°25′23.131″ E 8°23′53.531″ N. Map data from OpenStreetMaps http://openstreetmap.org/copyright.
Animals 13 00917 g0a5
Figure A6. Virtual enclosure version 6 dating from 11 August to 19 August. The dashed black line represents the physical fence surrounding the 8.2 ha grazing area. The marked green area represents the virtual enclosure. Only virtual enclosures active for more than 12 h were included. The coordinates of the upper left corner of the map frame are 55°25′23.131″ E 8°23′53.531″ N. Map data from OpenStreetMaps http://openstreetmap.org/copyright.
Figure A6. Virtual enclosure version 6 dating from 11 August to 19 August. The dashed black line represents the physical fence surrounding the 8.2 ha grazing area. The marked green area represents the virtual enclosure. Only virtual enclosures active for more than 12 h were included. The coordinates of the upper left corner of the map frame are 55°25′23.131″ E 8°23′53.531″ N. Map data from OpenStreetMaps http://openstreetmap.org/copyright.
Animals 13 00917 g0a6
Figure A7. Virtual enclosure version 7 dating from 19 August to 27 August. The dashed black line represents the physical fence surrounding the 8.2 ha grazing area. The marked green area represents the virtual enclosure. Only virtual enclosures active for more than 12 h were included. The coordinates of the upper left corner of the map frame are 55°25′23.131″ E 8°23′53.531″ N. Map data from OpenStreetMaps http://openstreetmap.org/copyright.
Figure A7. Virtual enclosure version 7 dating from 19 August to 27 August. The dashed black line represents the physical fence surrounding the 8.2 ha grazing area. The marked green area represents the virtual enclosure. Only virtual enclosures active for more than 12 h were included. The coordinates of the upper left corner of the map frame are 55°25′23.131″ E 8°23′53.531″ N. Map data from OpenStreetMaps http://openstreetmap.org/copyright.
Animals 13 00917 g0a7
Figure A8. Virtual enclosure version 8 dating from 27 August to 8 September. The dashed black line represents the physical fence surrounding the 8.2 ha grazing area. The marked green area represents the virtual enclosure. Only virtual enclosures active for more than 12 h were included. The coordinates of the upper left corner of the map frame are 55°25′23.131″ E 8°23′53.531″ N. Map data from OpenStreetMaps http://openstreetmap.org/copyright.
Figure A8. Virtual enclosure version 8 dating from 27 August to 8 September. The dashed black line represents the physical fence surrounding the 8.2 ha grazing area. The marked green area represents the virtual enclosure. Only virtual enclosures active for more than 12 h were included. The coordinates of the upper left corner of the map frame are 55°25′23.131″ E 8°23′53.531″ N. Map data from OpenStreetMaps http://openstreetmap.org/copyright.
Animals 13 00917 g0a8
Figure A9. Virtual enclosure version 9 dating from 8 September to 25 September. The dashed black line represents the physical fence surrounding the 8.2 ha grazing area. The marked green area represents the virtual enclosure. Only virtual enclosures active for more than 12 h were included. The coordinates of the upper left corner of the map frame are 55°25′23.131″E 8°23′53.531″N. Map data from OpenStreetMaps http://openstreetmap.org/copyright.
Figure A9. Virtual enclosure version 9 dating from 8 September to 25 September. The dashed black line represents the physical fence surrounding the 8.2 ha grazing area. The marked green area represents the virtual enclosure. Only virtual enclosures active for more than 12 h were included. The coordinates of the upper left corner of the map frame are 55°25′23.131″E 8°23′53.531″N. Map data from OpenStreetMaps http://openstreetmap.org/copyright.
Animals 13 00917 g0a9

Appendix B. Heatmaps of Cattle Movement in Virtual Enclosures

Figure A10. Heatmap series 1 visualising all Poll messages sent during the period from 4 July to 10 July. The green area shows the virtual enclosure area to which the bull calves were limited. The white to red colour indicates the number of Poll messages within a 10 m radius. The brightest colour represents the areas with the lowest density of Poll messages, while the darkest colour represents the areas with 250 or more Poll messages within the radius. The coordinates of the upper left corner of the map frame are 55°25′23.131″ E 8°23′3.531″ N. Contains data from the Agency for Data Supply and Infrastructure, Forårsbilleder Ortofoto-GeoDanmark, December 2022.
Figure A10. Heatmap series 1 visualising all Poll messages sent during the period from 4 July to 10 July. The green area shows the virtual enclosure area to which the bull calves were limited. The white to red colour indicates the number of Poll messages within a 10 m radius. The brightest colour represents the areas with the lowest density of Poll messages, while the darkest colour represents the areas with 250 or more Poll messages within the radius. The coordinates of the upper left corner of the map frame are 55°25′23.131″ E 8°23′3.531″ N. Contains data from the Agency for Data Supply and Infrastructure, Forårsbilleder Ortofoto-GeoDanmark, December 2022.
Animals 13 00917 g0a10
Figure A11. Heatmap series 2 visualising all Poll messages sent during the period from 10 July to 18 July. The green area shows the virtual enclosure area to which the bull calves were limited. The white to red colour indicates the number of Poll messages within a 10 m radius. The brightest colour represents the areas with the lowest density of Poll messages, while the darkest colour represents the areas with 250 or more Poll messages within the radius. The coordinates of the upper left corner of the map frame are 55°25′23.131″ E 8°23′53.531″ N. Contains data from the Agency for Data Supply and Infrastructure, Forårsbilleder Ortofoto-GeoDanmark, December 2022.
Figure A11. Heatmap series 2 visualising all Poll messages sent during the period from 10 July to 18 July. The green area shows the virtual enclosure area to which the bull calves were limited. The white to red colour indicates the number of Poll messages within a 10 m radius. The brightest colour represents the areas with the lowest density of Poll messages, while the darkest colour represents the areas with 250 or more Poll messages within the radius. The coordinates of the upper left corner of the map frame are 55°25′23.131″ E 8°23′53.531″ N. Contains data from the Agency for Data Supply and Infrastructure, Forårsbilleder Ortofoto-GeoDanmark, December 2022.
Animals 13 00917 g0a11
Figure A12. Heatmap series 3 visualising all Poll messages sent during the period from 18 July to 26 July. The green area shows the virtual enclosure area to which the bull calves were limited. The white to red colour indicates the number of Poll messages within a 10 m radius. The brightest colour represents the areas with the lowest density of Poll messages, while the darkest colour represents the areas with 250 or more Poll messages within the radius. The coordinates of the upper left corner of the map frame are 55°25′23.131″ E 8°23′53.531″ N. Contains data from the Agency for Data Supply and Infrastructure, Forårsbilleder Ortofoto-GeoDanmark, December 2022.
Figure A12. Heatmap series 3 visualising all Poll messages sent during the period from 18 July to 26 July. The green area shows the virtual enclosure area to which the bull calves were limited. The white to red colour indicates the number of Poll messages within a 10 m radius. The brightest colour represents the areas with the lowest density of Poll messages, while the darkest colour represents the areas with 250 or more Poll messages within the radius. The coordinates of the upper left corner of the map frame are 55°25′23.131″ E 8°23′53.531″ N. Contains data from the Agency for Data Supply and Infrastructure, Forårsbilleder Ortofoto-GeoDanmark, December 2022.
Animals 13 00917 g0a12
Figure A13. Heatmap series 4 visualising all Poll messages sent during the period from 26 July to 3 August. The green area shows the virtual enclosure area to which the bull calves were limited. The white to red colour indicates the number of Poll messages within a 10 m radius. The brightest colour represents the areas with the lowest density of Poll messages, while the darkest colour represents the areas with 250 or more Poll messages within the radius. The coordinates of the upper left corner of the map frame are 55°25′23.131″ E 8°23′53.531″ N. Contains data from the Agency for Data Supply and Infrastructure, Forårsbilleder Ortofoto-GeoDanmark, December 2022.
Figure A13. Heatmap series 4 visualising all Poll messages sent during the period from 26 July to 3 August. The green area shows the virtual enclosure area to which the bull calves were limited. The white to red colour indicates the number of Poll messages within a 10 m radius. The brightest colour represents the areas with the lowest density of Poll messages, while the darkest colour represents the areas with 250 or more Poll messages within the radius. The coordinates of the upper left corner of the map frame are 55°25′23.131″ E 8°23′53.531″ N. Contains data from the Agency for Data Supply and Infrastructure, Forårsbilleder Ortofoto-GeoDanmark, December 2022.
Animals 13 00917 g0a13
Figure A14. Heatmap series 5 visualising all Poll messages sent during the period from 3 August to 11 August. The green area shows the virtual enclosure area to which the bull calves were limited. The white to red colour indicates the number of Poll messages within a 10 m radius. The brightest colour represents the areas with the lowest density of Poll messages, while the darkest colour represents the areas with 250 or more Poll messages within the radius. The coordinates of the upper left corner of the map frame are 55°25′23.131″ E 8°23′53.531″ N. Contains data from the Agency for Data Supply and Infrastructure, Forårsbilleder Ortofoto-GeoDanmark, December 2022.
Figure A14. Heatmap series 5 visualising all Poll messages sent during the period from 3 August to 11 August. The green area shows the virtual enclosure area to which the bull calves were limited. The white to red colour indicates the number of Poll messages within a 10 m radius. The brightest colour represents the areas with the lowest density of Poll messages, while the darkest colour represents the areas with 250 or more Poll messages within the radius. The coordinates of the upper left corner of the map frame are 55°25′23.131″ E 8°23′53.531″ N. Contains data from the Agency for Data Supply and Infrastructure, Forårsbilleder Ortofoto-GeoDanmark, December 2022.
Animals 13 00917 g0a14
Figure A15. Heatmap series 6 visualising all Poll messages sent during the period from 11 August to 19 August. The green area shows the virtual enclosure area to which the bull calves were limited. The white to red colour indicates the number of Poll messages within a 10 m radius. The brightest colour represents the areas with the lowest density of Poll messages, while the darkest colour represents the areas with 250 or more Poll messages within the radius. The coordinates of the upper left corner of the map frame are 55°25′23.131″ E 8°23′53.531″ N. Contains data from the Agency for Data Supply and Infrastructure, Forårsbilleder Ortofoto-GeoDanmark, December 2022.
Figure A15. Heatmap series 6 visualising all Poll messages sent during the period from 11 August to 19 August. The green area shows the virtual enclosure area to which the bull calves were limited. The white to red colour indicates the number of Poll messages within a 10 m radius. The brightest colour represents the areas with the lowest density of Poll messages, while the darkest colour represents the areas with 250 or more Poll messages within the radius. The coordinates of the upper left corner of the map frame are 55°25′23.131″ E 8°23′53.531″ N. Contains data from the Agency for Data Supply and Infrastructure, Forårsbilleder Ortofoto-GeoDanmark, December 2022.
Animals 13 00917 g0a15
Figure A16. Heatmap series 7 visualising all Poll messages sent during the period from 19 August to 27 August. The green area shows the virtual enclosure area to which the bull calves were limited. The white to red colour indicates the number of Poll messages within a 10 m radius. The brightest colour represents the areas with the lowest density of Poll messages, while the darkest colour represents the areas with 250 or more Poll messages within the radius. The coordinates of the upper left corner of the map frame are 55°25′23.131″ E 8°23′53.531″ N. Contains data from the Agency for Data Supply and Infrastructure, Forårsbilleder Ortofoto-GeoDanmark, December 2022.
Figure A16. Heatmap series 7 visualising all Poll messages sent during the period from 19 August to 27 August. The green area shows the virtual enclosure area to which the bull calves were limited. The white to red colour indicates the number of Poll messages within a 10 m radius. The brightest colour represents the areas with the lowest density of Poll messages, while the darkest colour represents the areas with 250 or more Poll messages within the radius. The coordinates of the upper left corner of the map frame are 55°25′23.131″ E 8°23′53.531″ N. Contains data from the Agency for Data Supply and Infrastructure, Forårsbilleder Ortofoto-GeoDanmark, December 2022.
Animals 13 00917 g0a16
Figure A17. Heatmap series 8 visualising all Poll messages sent during the period from 27 August to 8 September. The green area shows the virtual enclosure area to which the bull calves were limited. The white to red colour indicates the number of Poll messages within a 10 m radius. The brightest colour represents the areas with the lowest density of Poll messages, while the darkest colour represents the areas with 250 or more Poll messages within the radius. The coordinates of the upper left corner of the map frame are 55°25′23.131″ E 8°23′53.531″ N. Contains data from the Agency for Data Supply and Infrastructure, Forårsbilleder Ortofoto-GeoDanmark, December 2022.
Figure A17. Heatmap series 8 visualising all Poll messages sent during the period from 27 August to 8 September. The green area shows the virtual enclosure area to which the bull calves were limited. The white to red colour indicates the number of Poll messages within a 10 m radius. The brightest colour represents the areas with the lowest density of Poll messages, while the darkest colour represents the areas with 250 or more Poll messages within the radius. The coordinates of the upper left corner of the map frame are 55°25′23.131″ E 8°23′53.531″ N. Contains data from the Agency for Data Supply and Infrastructure, Forårsbilleder Ortofoto-GeoDanmark, December 2022.
Animals 13 00917 g0a17
Figure A18. Heatmap series 9 visualising all Poll messages sent during the period from 8 September to 25 September. The green area shows the virtual enclosure area to which the bull calves were limited. The white to red colour indicates the number of Poll messages within a 10 m radius. The brightest colour represents the areas with the lowest density of Poll messages, while the darkest colour represents the areas with 250 or more Poll messages within the radius. The coordinates of the upper left corner of the map frame are 55°25′23.131″ E 8°23′53.531″ N. Contains data from the Agency for Data Supply and Infrastructure, Forårsbilleder Ortofoto-GeoDanmark, December 2022.
Figure A18. Heatmap series 9 visualising all Poll messages sent during the period from 8 September to 25 September. The green area shows the virtual enclosure area to which the bull calves were limited. The white to red colour indicates the number of Poll messages within a 10 m radius. The brightest colour represents the areas with the lowest density of Poll messages, while the darkest colour represents the areas with 250 or more Poll messages within the radius. The coordinates of the upper left corner of the map frame are 55°25′23.131″ E 8°23′53.531″ N. Contains data from the Agency for Data Supply and Infrastructure, Forårsbilleder Ortofoto-GeoDanmark, December 2022.
Animals 13 00917 g0a18

Appendix C. Herd Behaviour

Appendix C.1. Results of Herd Behaviour

There was a positive correlation between the number of warnings received by each individual animal and all the other animals. The positive correlation indicated that one bull calf receiving an auditory warning would to some degree lead to another bull calf also receiving an auditory warning (Figure A19). There was notable variation in the correlation coefficients, with results ranging from 0.13 < = r p < = 0.6 .
Figure A19. Pearson’s correlation coefficients of the pairwise number of warnings received by the bull calves. The colour ramp indicates the degree of correlation from white (0) to blue (1) with a darker blue indicating a stronger correlation. Correlation coefficients are only shown for significant results. The numbers around the edge of the matrix are the ID numbers of the bull calves.
Figure A19. Pearson’s correlation coefficients of the pairwise number of warnings received by the bull calves. The colour ramp indicates the degree of correlation from white (0) to blue (1) with a darker blue indicating a stronger correlation. Correlation coefficients are only shown for significant results. The numbers around the edge of the matrix are the ID numbers of the bull calves.
Animals 13 00917 g0a19
The median correlation between the number of warnings received by two individual calves in 8 h periods showed some temporal fluctuations, but seemed to be more or less constant over time. A linear regression model showed a slight increase in the median correlation (Pearson’s r = 0.4 + 5.13 · 10 4 · (8 h period)). The overall regression was statistically significant (R2 = 0.11, F(1,247) = 30.85, p < 0.001). However, due to the low R-squared value, the regression did not offer sufficient prediction accuracy to be accepted. Polynomial regression did not offer any significant increase in the R-squared value (R2 = 0.12, F(2,246) = 17.3, p < 0.001).

Appendix C.2. Discussion of Herd Behaviour

A positive Pearson’s correlation between the amount of warnings received was found for every calf pair bar one (Figure A19). This supports Hypothesis (3), which states that there would be a positive correlation between the number of warnings received by two random calves in 8-h periods. The positive correlations indicate that the animals in question were more likely to receive a higher number of warnings in time periods where their herd mates received numerous warnings. This suggests that they exhibit herd behavior, which has also been documented by Hubbard et al. [37]. If a negative correlation between the number of warnings received between two calves had been observed, it would mean that the animals were less likely to receive a warning whenever another calf had received a warning. Interestingly, the pairwise correlations varied considerably (Figure A19). One possible explanation for this variation could be the hierarchical structure of the herd, with any given calf mostly moving in close proximity to specific herd mates. However, a review of the existing literature on social dominance in beef cattle by Hubbard et al. [37] found little to no evidence for herd movement being related to social dominance rank. It is not possible to tell with certainty whether a positive correlation means the calves actually follow each other, as the test was performed using long time periods of 8 h. A shorter time period might provide a clearer picture.
Another explanation could be that the less correlated calf is relatively less interested in performing the same activity as the herd. It may be more likely to rest or graze at different times compared to the other animals. As all the correlations are still positive, it would, of course, still rest and graze mostly at the same time, but it is possible that it is less affected by the actions of the herd. This could also indicate a specific position in the hierarchy, or it could be down to the personality of the animal.
As a new way of analysing temporal data gathered from animals, such as from these calves, a sliding window correlation was run on the number of warnings any two calves had received in subsequent 8-h periods. A sliding window correlation makes it possible to capture any variations or trends in time, rather than relying on a single result independent of time. In our study, the median correlation between the number of warnings received by two animals did not change significantly over time using a sliding window correlation (Figure A20). Despite that there were no significant trends in the median correlation over time, the median correlation did seem to fluctuate noticeably, indicating that a simple correlation could prevent researchers from getting the full picture when analysing temporal data. A likely explanation for these fluctuations could be the decision to sum the warnings in 8-h periods and not split the sliding window correlation into three, such that all periods from the same day were tested against each other. Additionally, it could be hypothesized that the 8-h periods themselves were too large to capture short time variations. When performing a sliding window correlation, one has to consider the trade-off between the size of the window and the significance of the results. A smaller window allows for higher fidelity, with larger windows smoothing out possible changes or fluctuations. In this study, it was decided to use the smallest window size, which gave a median significance value of less than 0.05, to ensure that more than half of the correlations run were significant. However, one might argue that the window should have been even larger to increase the predictive accuracy. Future studies should continue to explore the use of sliding window correlations and other sliding window analyses to better capture temporal changes and trends.
Figure A20. Median sliding window correlation of the number of 8 hourly warnings between individuals over time with a window of 18 periods. The blue line is the result of linear regression (Pearson’s r = 0.4 + 5.13 · 10 4 · (8 h period)). Vertical dotted grey lines show the start time of each version of the virtual enclosure.
Figure A20. Median sliding window correlation of the number of 8 hourly warnings between individuals over time with a window of 18 periods. The blue line is the result of linear regression (Pearson’s r = 0.4 + 5.13 · 10 4 · (8 h period)). Vertical dotted grey lines show the start time of each version of the virtual enclosure.
Animals 13 00917 g0a20

References

  1. Ceausu, S.; Hofmann, M.; Navarro, L.M.; Carver, S.; Verburg, P.H.; Pereira, H.M. Mapping opportunities and challenges for rewilding in Europe. Conserv. Biol. 2015, 29, 1017–1027. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  2. Fraser, M.; Vallin, H.; Roberts, B. Animal board invited review: Grassland-based livestock farming and biodiversity. Animal 2022, 16, 100671. [Google Scholar] [CrossRef]
  3. Boone, R.B.; Hobbs, N.T. Lines around fragments: Effects of fencing on large herbivores. Afr. J. Range Forage Sci. 2004, 21, 147–158. [Google Scholar] [CrossRef]
  4. Tzanidakis, C.; Tzamaloukas, O.; Simitzis, P.; Panagakis, P. Precision Livestock Farming Applications (PLF) for Grazing Animals. Agriculture 2023, 2, 288. [Google Scholar] [CrossRef]
  5. Hayward, M.W.; Kerley, G.I. Fencing for conservation: Restriction of evolutionary potential or a riposte to threatening processes? Biol. Conserv. 2009, 142, 1–13. [Google Scholar] [CrossRef]
  6. Umstatter, C. The evolution of virtual fences: A review. Comput. Electron. Agric. 2011, 75, 10–22. [Google Scholar] [CrossRef]
  7. Harrington, J.L.; Conover, M.R. Characteristics of Ungulate Behavior and Mortality Associated with Wire Fences. Wildl. Soc. Bull. 2006, 34, 1295–1305. [Google Scholar] [CrossRef]
  8. Baines, D.; Summers, R.W. Assessment of Bird Collisions with Deer Fences in Scottish Forests. J. Appl. Ecol. 1997, 34, 941–948. [Google Scholar] [CrossRef] [Green Version]
  9. Mbaiwa, J.; Mbaiwa, O. The Effects of Veterinary Fences on Wildlife Populations in Okavango Delta, Botswana. Int. J. Wilderness 2006, 12, 17–23. [Google Scholar]
  10. Burkholder, E.N.; Jakes, A.F.; Jones, P.F.; Hebblewhite, M.; Bishop, C.J. To Jump or Not to Jump: Mule Deer and White-Tailed Deer Fence Crossing Decisions. Wildl. Soc. Bull. 2018, 42, 420–429. [Google Scholar] [CrossRef]
  11. MacDonald, A.M.; Jones, P.F.; Hanlon, J.A.; Martin, B.H.; Jakes, A.F. How did the deer cross the fence: An evaluation of wildlife-friendlier fence modifications to facilitate deer movement. Front. Conserv. Sci. 2022, 3, 991765. [Google Scholar] [CrossRef]
  12. Woodroffe, R.; Hedges, S.; Durant, S.M. To Fence or Not to Fence. Science 2014, 344, 46–48. [Google Scholar] [CrossRef] [PubMed]
  13. Brunberg, E.; Bøe, K.; Sørheim, K. Testing a new virtual fencing system on sheep. Acta Agric. Scand. Sect. A—Animal Sci. 2015, 65, 168–175. [Google Scholar] [CrossRef]
  14. Aaser, M.F.; Staahltoft, S.K.; Korsgaard, A.H.; Trige-Esbensen, A.; Alstrup, A.K.O.; Sonne, C.; Pertoldi, C.; Bruhn, D.; Frikke, J.; Linder, A.C. Is Virtual Fencing an Effective Way of Enclosing Cattle? Personality, Herd Behaviour and Welfare. Animals 2022, 12, 842. [Google Scholar] [CrossRef] [PubMed]
  15. Log, T.; Gjedrem, A.M.; Metallinou, M. Virtually Fenced Goats for Grazing Fire Prone Juniper in Coastal Norway Wildland-Urban Interface. Fire 2022, 5, 188. [Google Scholar] [CrossRef]
  16. Confessore, A.; Aquilani, C.; Nannucci, L.; Fabbri, M.C.; Accorsi, P.A.; Dibari, C.; Argenti, G.; Pugliese, C. Application of Virtual Fencing for the management of Limousin cows at pasture. Livest. Sci. 2022, 263, 105037. [Google Scholar] [CrossRef]
  17. Eftang, S.; Vas, J.; Holand, O.; Bøe, K. Goats are able to adapt to virtual fencing; A field study in commercial goat herds on Norwegian farms. Appl. Anim. Behav. Sci. 2022, 256, 105755. [Google Scholar] [CrossRef]
  18. Hamidi, D.; Grinnell, N.; Komainda, M.; Riesch, F.; Horn, J.; Ammer, S.; Traulsen, I.; Palme, R.; Hamidi, M.; Isselstein, J. Heifers do not care: No evidence of negative impact on animal welfare of growing heifers when using virtual fences compared to physical fences for grazing. Animal 2022, 16, 100614. [Google Scholar] [CrossRef]
  19. Stampa, E.; Zander, K.; Hamm, U. Insights into German consumers’ perceptions of virtual fencing in grassland-based beef and dairy systems: Recommendations for communication. Animals 2020, 10, 2267. [Google Scholar] [CrossRef]
  20. Campbell, D.L.M.; Lea, J.M.; Keshavarzi, H.; Lee, C. Virtual Fencing Is Comparable to Electric Tape Fencing for Cattle Behavior and Welfare. Front. Vet. Sci. 2019, 6, 445. [Google Scholar] [CrossRef] [Green Version]
  21. Sonne, C.; Alstrup, A.K.O.; Pertoldi, C.; Frikke, J.; Linder, A.C.; Styrishave, B. Cortisol in Manure from Cattle Enclosed with Nofence Virtual Fencing. Animals 2022, 12, 3017. [Google Scholar] [CrossRef]
  22. Lee, C.; Henshall, J.M.; Wark, T.J.; Crossman, C.C.; Reed, M.T.; Brewer, H.G.; O’Grady, J.; Fisher, A.D. Associative learning by cattle to enable effective and ethical virtual fences. Appl. Anim. Behav. Sci. 2009, 119, 15–22. [Google Scholar] [CrossRef]
  23. Lee, C.; Prayaga, K.; Reed, M.; Henshall, J. Methods of training cattle to avoid a location using electrical cues. Appl. Anim. Behav. Sci. 2007, 108, 229–238. [Google Scholar] [CrossRef]
  24. McSweeney, D.; O’Brien, B.; Coughlan, N.E.; Férard, A.; Ivanov, S.; Halton, P.; Umstatter, C. Virtual fencing without visual cues: Design, difficulties of implementation, and associated dairy cow behaviour. Comput. Electron. Agric. 2020, 176, 105613. [Google Scholar] [CrossRef]
  25. Lomax, S.; Colusso, P.; Clark, C.E.F. Does virtual fencing work for grazing dairy cattle? Animals 2019, 9, 429. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  26. Campbell, D.L.M.; Lea, J.M.; Haynes, S.J.; Farrer, W.J.; Leigh-Lancaster, C.J.; Lee, C. Virtual fencing of cattle using an automated collar in a feed attractant trial. Appl. Anim. Behav. Sci. 2018, 200, 71–77. [Google Scholar] [CrossRef]
  27. Costa, J.; von Keyserlingk, M.; Weary, D. Invited review: Effects of group housing of dairy calves on behavior, cognition, performance, and health. J. Dairy Sci. 2016, 99, 2453–2467. [Google Scholar] [CrossRef] [Green Version]
  28. Grignard, L.; Boissy, A.; Boivin, X.; Garel, J.; Le Neindre, P. The social environment influences the behavioural responses of beef cattle to handling. Appl. Anim. Behav. Sci. 2000, 68, 1–11. [Google Scholar] [CrossRef]
  29. Bouissou, M. Social relationships in domestic cattle under modern management techniques. Ital. J. Zool. 1980, 47, 343–353. [Google Scholar] [CrossRef]
  30. Doyle, R.; Moran, J. Cow Talk: Understanding Dairy Cow Behaviour to Improve Their Welfare on Asian Farms. CSIRO: Canberra, Australia, 2015. [Google Scholar]
  31. Marini, D.; Kearton, T.; Ouzman, J.; Llewellyn, R.; Belson, S.; Lee, C. Social influence on the effectiveness of virtual fencing in sheep. PeerJ 2020, 8, e10066. [Google Scholar] [CrossRef]
  32. QGIS. QGIS Desktop 3.26. 2022. Available online: http://www.qgis.org/ (accessed on 21 September 2022).
  33. R Core Team. R: A Language and Environment for Statistical Computing; R Foundation for Statistical Computing: Vienna, Austria, 2022. [Google Scholar]
  34. Reinhardt, V.; Mutiso, F.; Reinhardt, A. Social behaviour and social relationships between female and male prepubertal bovine calves (Bos indicus). Appl. Anim. Ethol. 1978, 4, 43–54. [Google Scholar] [CrossRef]
  35. Verdon, M.; Lee, C.; Marini, D.; Rawnsley, R. Pre-Exposure to an Electrical Stimulus Primes Associative Pairing of Audio and Electrical Stimuli for Dairy Heifers in a Virtual Fencing Feed Attractant Trial. Animals 2020, 10, 217. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  36. Markus, S.B.; Bailey, D.W.; Jensen, D. Comparison of electric fence and a simulated fenceless control system on cattle movements. Livest. Sci. 2014, 170, 203–209. [Google Scholar] [CrossRef]
  37. Hubbard, A.J.; Foster, M.J.; Daigle, C.L. Social dominance in beef cattle—A scoping review. Appl. Anim. Behav. Sci. 2021, 241, 105390. [Google Scholar] [CrossRef]
Figure 1. The island of Fanø located in the Wadden Sea of Denmark. The dashed black line represents the physical fence surrounding the 8.2 ha pasture. The coordinates of the upper left corner of the map frame are 55°25′23.131″ E 8°23′53.531″ N. This dataset includes Intellectual Property from European National Mapping and Cadastral Authorities and is licensed on behalf of these by EuroGeographics. Original dataset is available for free at https://www.mapsforeurope.org. Terms of the licence available at https://www.mapsforeurope.org/licenc All attribution statements can be found here.
Figure 1. The island of Fanø located in the Wadden Sea of Denmark. The dashed black line represents the physical fence surrounding the 8.2 ha pasture. The coordinates of the upper left corner of the map frame are 55°25′23.131″ E 8°23′53.531″ N. This dataset includes Intellectual Property from European National Mapping and Cadastral Authorities and is licensed on behalf of these by EuroGeographics. Original dataset is available for free at https://www.mapsforeurope.org. Terms of the licence available at https://www.mapsforeurope.org/licenc All attribution statements can be found here.
Animals 13 00917 g001
Figure 2. Heatmap visualising all Poll messages sent during the period from 25 September to 3 October. The green area shows the virtual enclosure area to which the bull calves were limited. The white to red colour scale indicates the number of Poll messages within a 10 m radius on a discrete scale, ranging from 1 point to above 250 points. The coordinates of the upper left corner of the map frame are 55°25′23.131″ E 8°23′53.531″ N. Contains data from the Agency for Data Supply and Infrastructure, Forårsbilleder Ortofoto-GeoDanmark, December 2022.
Figure 2. Heatmap visualising all Poll messages sent during the period from 25 September to 3 October. The green area shows the virtual enclosure area to which the bull calves were limited. The white to red colour scale indicates the number of Poll messages within a 10 m radius on a discrete scale, ranging from 1 point to above 250 points. The coordinates of the upper left corner of the map frame are 55°25′23.131″ E 8°23′53.531″ N. Contains data from the Agency for Data Supply and Infrastructure, Forårsbilleder Ortofoto-GeoDanmark, December 2022.
Animals 13 00917 g002
Figure 3. Boxplot of the duration of time spent outside the virtual enclosure for each breakout event during the 89-day period. The points represent each breakout and the colours indicate the escaped individual’s ID number. The boxplot shows min, Q1, median, Q3 and max. Outliers are calculated as 1.5 times IQR. The points are mainly clustered between 0 and 2 min with a few breakouts ranging from 13 to 33.5 min and a single breakout being an outlier at 66 min.
Figure 3. Boxplot of the duration of time spent outside the virtual enclosure for each breakout event during the 89-day period. The points represent each breakout and the colours indicate the escaped individual’s ID number. The boxplot shows min, Q1, median, Q3 and max. Outliers are calculated as 1.5 times IQR. The points are mainly clustered between 0 and 2 min with a few breakouts ranging from 13 to 33.5 min and a single breakout being an outlier at 66 min.
Animals 13 00917 g003
Figure 4. Cumulative number of warnings received by each bull calf between 4 July 2022 and 30 September 2022. The ID Number refers to each of the 17 calves.
Figure 4. Cumulative number of warnings received by each bull calf between 4 July 2022 and 30 September 2022. The ID Number refers to each of the 17 calves.
Animals 13 00917 g004
Figure 5. Cumulative number of impulses received by each bull calf between 4 July 2022 and 30 September 2022. The ID Number refers to each of the 17 calves.
Figure 5. Cumulative number of impulses received by each bull calf between 4 July 2022 and 30 September 2022. The ID Number refers to each of the 17 calves.
Animals 13 00917 g005
Figure 6. Warning duration in seconds over time for all animals represented as a series of boxplots showing min, Q1, median, Q3 and max. Outliers calculated as 1.5 times IQR.
Figure 6. Warning duration in seconds over time for all animals represented as a series of boxplots showing min, Q1, median, Q3 and max. Outliers calculated as 1.5 times IQR.
Animals 13 00917 g006
Figure 7. Number of impulses per warning over time. Ratio shown as boxplots showing min, Q1, median, Q2 and max. Outliers calculated as 1.5 times IQR.
Figure 7. Number of impulses per warning over time. Ratio shown as boxplots showing min, Q1, median, Q2 and max. Outliers calculated as 1.5 times IQR.
Animals 13 00917 g007
Table 1. Overview of the different message types and description of frequency, data type and number of observations.
Table 1. Overview of the different message types and description of frequency, data type and number of observations.
TypeFrequencyData TypeNo. of Obs.
PollEvery 15 minPositional151,576
SeqEvery 30 minActivity71,990
WarningAfter receiving warningPositional13,906
StatusUpon status change (see Table 2)Status2153
ZapUpon receiving impulsePositional1538
Table 2. Overview and description of the different fence statuses.
Table 2. Overview and description of the different fence statuses.
Fence StatusDescription
NormalThe virtual enclosure is active
Not StartedThe virtual enclosure is not active
Maybe Out Of FenceThe collar position is outside the enclosure, but the GPS accuracy is not adequate to a trigger warning
EscapedThe animal has escaped the enclosure
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Staahltoft, S.K.; Aaser, M.F.; Jensen, J.N.S.; Zadran, I.; Sørensen, E.B.; Nielsen, A.E.; Alstrup, A.K.O.; Bruhn, D.; Linder, A.C.; Sonne, C.; et al. The Effectiveness of Virtual Fencing of Bull Calves in a Holistic Grazing System. Animals 2023, 13, 917. https://0-doi-org.brum.beds.ac.uk/10.3390/ani13050917

AMA Style

Staahltoft SK, Aaser MF, Jensen JNS, Zadran I, Sørensen EB, Nielsen AE, Alstrup AKO, Bruhn D, Linder AC, Sonne C, et al. The Effectiveness of Virtual Fencing of Bull Calves in a Holistic Grazing System. Animals. 2023; 13(5):917. https://0-doi-org.brum.beds.ac.uk/10.3390/ani13050917

Chicago/Turabian Style

Staahltoft, Søren Krabbe, Magnus Fjord Aaser, Jakob Nødgaard Strange Jensen, Ismat Zadran, Emil Birkmose Sørensen, Anders Esbjerg Nielsen, Aage Kristian Olsen Alstrup, Dan Bruhn, Anne Cathrine Linder, Christian Sonne, and et al. 2023. "The Effectiveness of Virtual Fencing of Bull Calves in a Holistic Grazing System" Animals 13, no. 5: 917. https://0-doi-org.brum.beds.ac.uk/10.3390/ani13050917

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop