Head and body orientation of the White Stork Ciconia ciconia during incubation: effect of wind, apex predators and power lines

Incubation behaviour is essential for understanding the reproductive success in birds. For example, the orientation of the bird is important for reducing incubation costs associated with wind or sun, but on the other hand can be modified by the perceived risk of predation. We studied the body position of incubating White Stork Ciconia ciconia in eastern Poland using a small unmanned aerial vehicle (drone). The head and body orientation of the incubating storks was non-random and modified by natural factors, mainly wind direction and speed, but also by the presence of an apex predator, the White-tailed Eagle Haliaeetus albicilla. However, head orientation during incubation in nests located on electricity poles was also modified by the presence of the power lines, probably due to disturbance in the magnetic field detected by birds. Surprisingly, although the positioning of incubating birds (mainly females) is very important for the detection of predators and for reducing energy costs, these have not previously been studied. New technologies, such as drones, make it possible to collect new, extensive information on the incubation behaviour of birds.


Introduction
Incubation behaviour is central to the reproductive success of birds (Carey 1980;Deeming and Reynolds 2015). Incubating adults must ensure that eggs are kept within a relatively narrow temperature range to promote the growth of developing embryos (Carey 1980;Deeming and Reynolds 2015). Therefore, incubation should be, and is, modified by the behaviour of the adult birds depending on weather conditions (Weller 1958;Gochfeld 1978). A particularly influential factor on the body orientation of incubating birds is wind direction and its strength, substantially modifying the behaviour of birds, and in particular turning the head and body to reduce both air resistance and the energy costs of incubation (Poussart et al. 2001;Baoqing et al. 2004). Incubation is time-consuming and takes place in a practically stationary position, by adult birds, most often the female, sitting on the nest. Such a permanent, stationary position with limited visibility (Graham 2011) makes the incubating bird relatively safe from attack by a predator (Carey 1980;Bakner et al. 2019).
In contrast, the White Stork Ciconia ciconia is a species that builds large, open nests, very often on the top of manmade structures, and less often on trees, which are clearly visible against the skyline (Creutz 1988;Tryjanowski et al. 2009;Zbyryt et al. 2020). Incubation in this species takes 32 (29-34) days and usually starts with the second egg (Bocheński and Jerzak 2006;Wuczyński 2012). Both mates incubate, but the female contributes more than the male (Bocheński and Jerzak 2006;Kosicki 2010). The White Stork is a large species (after Creutz (1988)-body mass: 2.3-4.5 kg, wing span: 155-215 cm), but sometimes is the subject of attacks by large birds of prey. However, the White Stork has to fear only a few predators such as the White-tailed Eagle Haliaetus albicilla, the golden eagle Aquila chrysaetos and Bonelli's Eagle Hieraaetus fasciatus. These raptors in the vicinity of the nest always evoke a very strong reaction. Less often the appearance of smaller raptor species, such as the Common Buzzard Buteo buteo or female of the Northern Goshawk Accipiter gentilis, sometimes elicits the anxiety stretch display (Bocheński and Jerzak 2006). In Central and Eastern Europe, confirmed cases of predating storks, both young and adults, are known only by white-tailed eagles (Zawadzka et al. 2006;Dementavičius et al. 2020).
Recently, the White Stork has increased the proportion of nests located on electricity poles (Tryjanowski et al. 2009;Vaitkuvienė and Dagys 2014), or even on special poles not connected to the power network (Kosicki 2010;Zbyryt et al. 2021). The flowing electrical current disturbs the reception of the magnetic field (Balmori 2015) and also affects the breeding parameters of birds, including the White Stork (Vaitkuvienė and Dagys 2014). It is, therefore, possible to expect, by analogy with observations of grazing, resting or defecating mammals (Burda et al. 2009;Yosef et al. 2020), changes in the alignment of the incubating bird towards the North-South axis. This may be especially observable in sitting, incubating birds on electricity poles, because their head, where the magnetoreceptors are located (Mouritsen et al. 2004), is very close to wires with a flowing electrical current.
Additionally, the White Stork is an interesting species for checking the effect of various factors on incubation behaviour from a technical point of view. As mentioned above, storks build highly visible large nest platforms and incubating birds are easy to monitor, by unmanned aerial vehicles (drones) and it is possible to obtain a large sample size in a short period of time without disturbance (Zbyryt et al. 2020).
Therefore, the main aim of the study was to describe the body orientation of incubating birds during the daytime and to establish the influence of external factors, such as wind, predator presence, and nest location (electricity poles vs. other structures) on incubation behaviour. We suggest that this information may help in better understanding breeding behaviour, ecology and, in consequence, the population trajectory of the White Stork, which are especially important in light of new nesting sites (artificial poles) provided by conservationists, both in this particular area, as well as across the geographical range of the species (Muzinic 1999;Zbyryt et al. 2021).

Study site
The research was conducted between May 17 and May 20, 2020 in north-eastern Poland. The largest population of White Storks in the country (Zbyryt et al. 2014, A. Zbyryt-unpublished data) occurs in this region.

Fieldwork
Data were collected using a small quadrocopter Phantom 4 (DJI, Shenzhen, China) with an in-built 12 Mp camera (20 mm lens), coloured white. The detailed specification of this device is described by Zbyryt et al. (2020).
Flights were made to 181 nests in which incubating storks were found: 81 were on electricity poles, 68 on free-standing poles, 23 on buildings and 9 on trees. All electricity poles with nests were connected to the low-voltage power grid (0.4 kV), and those on trees were at their top and were not sheltered by branches. Each nest was recorded only once. Take-off was performed from a distance of at least 30 m, as this has been shown to be a distance that does not cause disturbance of White Storks on nests (Zbyryt et al. 2020). An overhead photo was taken of each nest from a height of 40 m above the ground when the incubating bird was exactly in the central position of the monitor, which was helped by the arrangement of grid lines displayed on the screen of the MediaPad M3 tablet with a diagonal of 8.4″ (Huawei, Shenzhen, China). Within one flight, between 1 and 17 nests were photographed (mean = 2.4; SD = 2.7). Wind strength and direction were measured on an ongoing basis using the Meteo ICM-Numerical Weather Forecast application (University of Warsaw; www. meteo. pl), which was read for the nearest town for which weather data were given (no more than 15 km distant). The wind direction was verified in the field using a small flag and a compass built into a smartphone. During the survey period, we recorded winds only from three directions: WNW, NW, WSW and these categories were used in later analysis.

Photo processing
The orientation of the incubating bird (an axis along the bird's back and beak relative to north) was determined by means of a circular overlay on the photo in the form of a circle with an angular scale and auxiliary lines, made in vector graphics. The orientation of power lines was determined in the same way. Each electricity pole on which the nest was located had from 1 to 3 electrical lines (each of 3-5 parallel cables, Fig. 1) "connecting" under the nest. To determine the impact of breeding white-tailed eagles on the position of an incubating stork, data on the location of eagle nests in the area covered by the study was obtained from the database of the Eagle Conservation Committee. White-tailed eagles usually hunt within 4 km (range 2.5-15 km) of their nest in the breeding season (Mizera 2015;Dementavičius et al. 2020). We assumed that these raptors could cause persistent reactions in storks nesting within this radius around the eagle nest, covering the area that eagles most frequently use (mentioned 4 km). Assuming that the White-tailed Eagle would approach in a straight line from their nest, the storks nesting in raptor territory were classified into two categories: (1) with the eagle direction in the field of view, or (2) in the blind sector. The blind sector in the White Stork at the back of its head amounts to 72° (Graham 2011) and we used this range in our study. The location and time of day of each nest image were read from the Exchangeable Image File Format.

Statistical analysis
We started our analysis by building the model with body orientation as the dependent variable. As independent variables we tested: Wind direction (three levels: WNW, NW, WSW), wind gusts (in meters per second, range 11-17), wind speed (metres per second, range 4-10), eagle presence (two levels: yes, no), nest type (four levels: electricity pole, free-standing pole, tree, building) and time of day. In total we analyzed 181 records. The orientations of the incubating birds are typical circular data whose main characteristic is the periodic nature (0° is identical to 360°) that makes them fundamentally different from interval data which would treat the difference between 0° and 360° as 360° (Cremers and Klugkist 2018). In our first approach, we used a projected normal circular general linear model (PN GLM) (Cremers and Klugkist 2018) which is a Bayesian method. To estimate model parameters we had to define parameters for the Markov Chain Monte Carlo (MCMC) algorithm (number of iteration, its = 10,000; burn-in period, burn = 100; lag, n.lag = 3). All models were performed in the package bpnreg (Cremers 2018) in R software (R Development Core Team 2018). Circular plots were made in the package circular (Agostinelli and Lund 2017). We compared competing models using the deviance information criterion (DIC). A better fitting model is represented by lower values of DIC; two competing models can be compared where a difference ≥ 2 DIC indicates a significantly better model.
For the subset of nests on electricity poles (n = 81) we wanted to analyze whether birds oriented themselves along power lines. The electrical lines consisted of 3 to 5 parallel cables, and further we did not distinguish a number of wires on each of the line, but we counted only the number of lines "connecting" under the nest (see Fig. 1B for an example of two lines of two cables each). Only 11 nests had one line, 40 had two lines, and 30 nests had three lines. In the latter two cases, we had to select one of the lines for line orientation. We chose this according to its closeness to the head of the bird, because the magnetic sensory organ is located in the head in migratory birds (Mouritsen et al. 2004;Heyers et al. 2007). However, the selection of the closest line could have biased our results, e.g. where there were 3 lines we selected that which was closest to the bird's head. To check this approach, we tested whether randomly generated bird orientations were correlated with the orientation of the selected line. If the randomly generated orientation was correlated with that of the selected lines it means that the observed correlation of bird orientation could be the result of chance not bird preference. To test it we used circular correlation (circular version of the Pearson's product-moment correlation in the package circular). The correlation between random bird orientations and those of power lines was not significant (Pearson's r = 0.004, t-test = 0.04, P < 0.970). Finally, we used the same statistical procedure as for the full data set. We compared the full model with lines as a new explanatory variable (M9 in Table 9) with the model without this variable (M10).

Results
The head orientation of the incubating storks was positioned in all compass directions ( Fig. 2A) but were modified by external factors. The best model included wind direction, wind speed and eagle presence. We found that the type of nest did not have a significant effect on the bird's orientation (M1 vs M2). Wind direction significantly increased the fit of the model (M1 vs M5). If the wind was blowing from the NWW the mean orientation of the birds was 267° and the mean resultant length (0 indicates that the spread of data points was very large, 1 indicates that the data were aggregated as a single value) was 0.56 (Fig. 2B). If the wind was blowing from the NW the mean orientation of the birds was 331° and the mean resultant length was 0.45 (Fig. 2C). If the wind was blowing from the SWW the mean orientation of the birds was 166° and the mean resultant length was 0.15 (Fig. 2D). Model parameters indicate that NWW was significantly different from NW and SWW. Further, the model fit decreased significantly (M5 vs M3) by adding the variable time of day. We found that the model with the variable gusts was a significantly better fit (M5 vs M6), however, by adding the variable wind speed (which was highly correlated with gusts; Spearman Rank Correlation r s = 0.64, P < 0.001) improved the model fit (M7 vs M5) so we included this variable in the final model. However, we did not find any circular coefficient of wind speed to be significant (highest posterior density interval of coefficient included zero); only the linear coefficient (β = 0.167 ± 0.076 SD) was significant and it can be interpreted that increased wind speed increased the radius of the bird orientation (birds aligned more north-south). The fit of the model significantly increased (M8 vs M7) by adding eagle presence. However, there was no difference between the estimated circular means. This is due to a different spread of data points; if an eagle was not present (Fig. 3A) the mean orientation was 287° and the mean resultant length was 0.74. If an eagle was present (Fig. 3B) the mean orientation was 298° and the mean resultant length was 0.24 which indicates a high spread of data points.
In the second data set where the lines variable was added we found that this model was significantly better than without (M9 vs M10; Table 1). However, the proposed projected normal circular general linear model (PN GLM) did not account for the circular-circular relation. To show the relationship we performed a simple circular correlation. The correlation between bird orientations and those of power lines was statistically significant (Pearson's r = 0.606, t-test = 5.54, P < 0.001) and can be interpreted as birds aligning with the power lines (Fig. 4).

Discussion
Given the fact that the White Stork is a well-studied species, it is surprising that the phenomenon we have examined had not previously been investigated. The results of our studies, although it may seem in some cases to be quite obvious, have interesting behavioural and conservation implications, as discussed below.
The head and body orientation of the incubating storks was non-random. As expected, several factors modified it. The strongest effect was caused by the prevailing wind, and depending on its strength and direction, the birds aligned their head and body in such a way as to provide the least possible air resistance, taking aerodynamic positions, as has also been shown in other species, e.g. the Common Tern Sterna hirundo (Gochfeld 1978), the Greater Snow Goose Chen caerulescens atlantica (Poussart et al. 2001), all species nesting on exposed sites, where winds can be particularly severe. This phenomenon must be very common among White Storks, considering that they nest on elevated structures in an open landscape. Such behaviour must be adaptive to such conditions, the more so since the stork's nest is quite flat and does not provide protection against the wind. This is confirmed by opportunistic observations made by the authors and other stork researchers, but this phenomenon has been described here for the first time. However, that aerodynamic position is probably also more common for birds using cup-nests located in exposed places, such as the top of trees, buildings or electricity poles (e.g. ospreys, vultures, white-tailed eagles, ravens).
In our research, we also showed that the orientation of the incubating female (predominantly) can be modified by the appearance of an apex predator, in this case, the White-tailed Eagle (Dementavičius et al. 2020). We see some limitations in the results of our research. An incubating bird does not need to see the incoming predator quickly, but its effect is noticed by other storks, including the partner and is communicated by a specific nervous behaviour (Jakubiec and Peterson 2005;Bocheński and Jerzak 2006). Individuals living in high densities do not have to constantly react to the threat from the predator when other factors come into play constantly (e.g. wind). They can be based on the manyeyes hypothesis, which assumes that the cost of scanning the environment is spread over many individuals (Lima 1995). The predator can approach from different directions, so it is difficult to expect a directional reaction, and the incubating bird has only a limited field of view (Graham 2011), so the body layout of the incubating female was also modified by other factors such as the wind direction (which interacts constantly and permanently). However, at least a few cases are known from the study area of the attack and direction of attack of white-tailed eagles on storks in the nests. The raptor came in a straight line from their nest site, grabbed the adult or chick and returned in a straight line to their nest.
The situation was repeated until all the chicks were predated (A. Zbyryt-unpublished data). This is rather unsurprising, taking into account the fact that the nests of White Storks are clearly visible and often function for decades in one place. Raptors can easily learn their location, distribution, occupation or breeding stage. This would also confirm our methodological assumption of how eagles depredated storks in their nests. In our study, the occupied nests that were within the range of the White-tailed Eagle were highly dispersed. We do not know if it is an effect of raptor pressure or a feature of the local population of White Storks. However, in this case, we can assume that storks have to rely on their own risk assessment or possibly their partner when he is near the nest. Storks usually forage within one km around the nest (Dziewiaty 1992;Alonso et al. 1994;Ożgo and Bogucki 1999), but in search of food they can fly up to 5 km, which depends on agricultural activity (Johst et al. 2001). Therefore, in conditions of increased predation pressure, the reliance on the partner as an additional source of threat detection is very limited. Then the incubating female may turn her head more often to widen the field of view. The problem, however, is how often this should be done to avoid being caught off guard by the eagle. Under such conditions, such behaviour can lead to increased stress (Boonstra 2013). To reduce it, the incubating female can arrange her body in such a way as to be able to constantly monitor the area from which a predator may likely come. We believe that it is the much higher cost associated with the stress caused by the predator (a permanent factor) that contributes to the fact that an incubating stork responds more often with an appropriate body position in response to this factor than the negative impact of the wind (an occasional factor). Hence, it seems to us that the result obtained by us shows the real effect of shaping a new behavioural response to predatory pressure.  This also opens up an interesting area for discussion of how two endangered bird species, each of which is targeted with practical conservation assistance (Treinys et al. 2011), affect one another. The growing population of White-tailed Eagle in central and eastern Europe is already provoking discussions about the impact on other birds of prey and on the black stork Ciconia nigra (Zawadzka et al. 2006;Treinys et al. 2011). A similar effect is now also shown in the White Stork (Dementavičius et al. 2020), a species nesting in human settlements and traditionally considered a species under low predatory pressure (Creutz 1988;Kosicki 2010). Perhaps, however, this low level of pressure has only been a feature of the last few decades, when birds of prey were persecuted and therefore their impact on stork populations was not apparent. However, from a historical and evolutionary point of view, the trophic relations (predator (eagle)prey (stork)) are probably very old and are still documented in Africa (Parsons 1977). In the case of the population we studied and nearby, we have evidence that white-tailed eagles hunt White Storks (Dementavičius et al. 2020), from occasional observations by ornithologists to photos of the raptors' prey from their nests, where photo traps have been installed. The study of the food composition of the Whitetailed Eagle in the area of our study showed that the White Stork (mostly adults) is an important component of the raptor's diet-it was found in half of the nests. White stork was second in food composition in terms of biomass after Common Carp (Cyprinus carpio). White-tailed eagles hunt for storks mainly in spring, i.e. in the period that overlaps with the period of incubation of eggs by White Storks (Komar and Mirski 2019).
This also shows the formation of a "landscape of fear" among the local stork population-birds, despite the lack of direct exposure to attack, begin to modify their behaviour to avoid predator attack (Beauchamp 2017). Future research is needed to understand how common this phenomenon is, how it is learned to avoid predators, and how this information is spread between individuals (Berger 2001). The fact that only birds in the immediate vicinity of the White-tailed Eagle react may indicate that this is a very local phenomenon, limited to specific pairs and other individuals nesting within their visual range. Storks are strongly attached to their nests, which may limit the transfer of such knowledge if they do not live in high density.
Many animals tend to align their body axes in the geomagnetic North-South direction (Burda et al. 2009;Yosef et al. 2020). In the case of a long process of incubating eggs by storks, magnetic alignment could be expected, due to the fact that birds, especially long-distance migrants, have the ability to detect the Earth's magnetic field (Wiltschko and Wiltschko 1996;Heyers et al. 2007). However, we did not notice this phenomenon, but we have found other interesting issues that the low-voltage power lines influence the positioning of the bird's body on the nest. Electrical power lines generate an alternating magnetic field. It has been proven so far that productivity is the lowest in White Stork nests located directly on low-voltage power lines. The effect was stronger the closer the nest was to the wires (Vaitkuvienė and Dagys 2014). The explanation for why incubating birds align their heads to power lines may be related to the fact that migratory birds use their visual system to perceive the reference compass direction of the geomagnetic field. In other words, migratory birds, which include the White Stork, "see" the reference compass direction provided by the geomagnetic field (Heyers et al. 2007). It is true that the birds did not migrate at the time of or study, but this ability must be permanent and somehow cause a strong behavioural response of storks to this phenomenon. In our case, the alignment also had a stronger effect, as in terms of breeding parameters (Vaitkuvienė and Dagys 2014). This also shows how complex and multifaceted the impact of eclectic smog on animals is (Balmori 2015). Finally, it is also worth noting that the orientation during incubation in conditions of greater population density, that is, in the case of colony nesting of the White Stork (Turjeman et al., 2016) could be further modified, at least at the beginning of incubation when females are still in their fertile period, by the desire to guard the partner and thus avoid extra-pair copulations. However, this would require research in colonies of this species, for example on the Iberian Peninsula or in northern Poland (e.g. Cabodevilla and Aguirre 2019; Zbyryt et al. 2014). In our case, it did not matter, because to avoid this effect, we carried out our research in a narrow time window in the final phase of egg incubation, just before hatching of the first chicks.