Flight Type and Seasonal Movements Are Important Predictors for Avian Collisions in Wind Farms

Simple Summary: Wind farms are an alternative energy source mitigating environmental pollution. However, they can have adverse effects, causing an increase in mortality for wildlife through collision with wind turbines. In this work


Introduction
Many species of vertebrates are experiencing significant population declines, and determining the causes can lead to take the most appropriate management actions [1]. While the main drivers are habitat loss, invasive species and climate change, there are specific problems that cause significant direct mortality on wild birds [2]. A notable ( Figure 1). Within Spain, the administrative division Castilla y León has 243 wind farms with an amount of 5593 MW, representing 24.2% of the installed power in Spain [32]. At the beginning of the 21st century, the wind farms were mainly installed in the mountainous areas but, thanks to technological improvement, wind farms have recently been installed on the plains in the centre of the region. Due to the change, the installed turbines now affect not only birds from mountainous areas but all bird species. This study analysed birds which have been killed by wind turbines and been admitted to two rescue centres in Castilla y León. All the specimens were found dead under wind turbines after colliding with them in wind farms, and were collected by the companies that carry out the official monitoring and by environmental agents until their arrival at the rescue centers of the environmental department. The monitoring of bird mortality at the foot of each wind turbine was carried out along the line of the wind turbines, with an impact of 100 m. on each side and with a search effort of 20 min per wind turbine covered on foot. The precise places, the date and the state in which the remains of birds were found was noted, immediately notifying the environmental department of Junta de Castilla y León to proceed with the collection by its personnel. In addition, the coordinates of the place were taken by GPS. Once there, a veterinarian thoroughly examined the specimens and confirmed their death due to a collision with the wind This study analysed birds which have been killed by wind turbines and been admitted to two rescue centres in Castilla y León. All the specimens were found dead under wind turbines after colliding with them in wind farms, and were collected by the companies that carry out the official monitoring and by environmental agents until their arrival at the rescue centers of the environmental department. The monitoring of bird mortality at the foot of each wind turbine was carried out along the line of the wind turbines, with an impact of 100 m. on each side and with a search effort of 20 min per wind turbine covered on foot. The precise places, the date and the state in which the remains of birds were found was noted, immediately notifying the environmental department of Junta de Castilla y León to proceed with the collection by its personnel. In addition, the coordinates of the place were taken by GPS. Once there, a veterinarian thoroughly examined the specimens and confirmed their death due to a collision with the wind turbines. For more information about the region and the organisation of wildlife rehabilitation centres in Castilla y León, see Balmori (2019) [33].
A database was prepared including the scientific name, date, group, seasonal movements, IUCN red list category, flight type, and length, as well as the number of Spanish and European population (pairs) of the dead birds arriving at these rescue centres for 16 years (2001-2016). The sex data of some specimens were also included.
Seasons were defined considering the week number of collisions; being 49-53 and 1-8 winter, 9-22 spring, 23-35 summer, 36-48 autumn. For the seasonal movement type of each species, three categories were considered: aestival, sedentary and partial migration. Aestival is a species which only remains in Europe in spring and summer, and winters south of the Sahara Desert (long-distance migrants crossing the Sahara, e.g., the Short-toed Snake Eagle Circaetus gallicus). Wintering species in the Iberian Peninsula that remain throughout all the year are considered sedentary, e.g., the Golden Eagle (Aquila chrysaetos). Pre-Saharan migrants, that usually make shorter and more variable displacements than trans-Saharan migrants, are considered partially migrants (short-distance migrants, not crossing the Sahara, e.g., the Eurasian Sparrowhawk Accipiter nisus). The Iberian Peninsula receives many European pre-Saharan migrants (partially migrants) in winter. However, there were species with several categories and, in those cases the predominant category in the Iberian Peninsula was considered, following Peterson et al. (1973) [34].
Each species was classified according to his predominant flight type, within the following categories: flapping, active soaring, passive soaring, hovering, song-flights, and high-speed. Since some species have several flight types and use different flight speeds [35], the flight type considered more typical of the species was selected. For instance, in Common Kestrel (Falco tinnunculus) the predominant flight types are flapping and hovering, but hovering was selected. Similarly, in Skylark (Alauda arvensis), the predominant flight types are flapping and songflights, but songflights was selected.
For size, the length of each species (measurement between the tip of the bill and the tip of the tail) following Mullarney et al. (2000) [36] was used. In cases where a range of lengths existed within the same species or there was a size differences between sexes, the mean value was used. For each species, the number of breeding couples from Spanish and European populations was obtained following Heath and Hagemeijer (2000) [37].

Statistical Analyses
Linear models were created in order to detect factors influencing the number of collisions occurring in every species. First, a model including all species was constructed. In addition, we also fit another model using a smaller dataset in order to corroborate the most important predictors (see below and also ANOVA analyses). This second linear model included only bird of prey and passerine having more than four collisions. These groups accounted for the large majority of collisions and the purpose was determining more in-depth patterns in these cases avoiding anecdotal data. In addition, the Griffon Vulture (Gyps fulvus) was the species with the highest number of collisions (n = 1901). Due to its high number of collisions, this species was not included in the second model to avoid bias, but was considered independently. Finally, a total of 29 out of 101 species belonging to the Bird of prey and Passerine groups remained in the second model.
Quantitative variables were log10 transformed for normalizing values distributions and the test assumptions both for the linear regression (linear relationship and homoscedasticity of residuals) and the ANOVA analyses (|skeweness and kurtosis| < 1, homoscedasticity) were checked. We performed a stepwise model selection by Akaike Information Criterion (AIC) in R software [38], including the following variables, without considering interactions: taxonomic order or group (for detecting potential phylogenetic influence), seasonal movements, flight type, length, and the number of couples for the Spanish and European populations. Sex was not included in the model due to the few sexed specimens. The stepwise algorithm search was performed in both directions (backward and forward) and the model selected only relevant predictors to estimate the number of collisions. Once we knew the important predictors, the coefficients and the percent of variance explained by each relevant variable were determined by the linear model, and real values were compared with predicted values to visualize the goodness of fit. The functions were applied to the species used to construct the model. Based on the result of the model including all the species, we obtained a species collision risk index.
A two-way ANOVA was also performed with the two more important predictors selected by the model. Furthermore, repeated-measures ANOVA was carried out to study the collision seasonality and to relate it to the seasonal movement type of the different species.

Results
During the 16 years of the study period, a total of 3130 specimens belonging to 101 bird species died by collisions with wind turbines in Castilla y León. The total number of collisions and their seasonality along with the other data for each species are included in Table 1.
First, a model including all the species registered was adjusted and a total of three of the studied variables were included by AIC in the stepwise regression. Flight type was the most important variable, explaining 35% of the total variability, followed by seasonal movement type (4%) and the Spanish population (4%). The function: "Predicted number of collisions (log 10) = 0.31 + 0.28 * Flight type_Active soaring-0.67 * Flight type_Flapping + 1.01 * Flight type_Hovering + 1.03 * Flight type_Passive soaring + 0.16 * Flight type _Songflights + 0.06 * Movements_Aestival + 0.29 * Movements_Partial migration + 0.11 * Spanish population (log10 + 1)", was applied to the whole dataset of 101 species in order to estimate the predicted number of collisions. For instance, if we consider the Wood Lark (Lullula arborea) case (with Sonflights, Partial migration and a Spanish population of 900,000 individuals), the Predicted number (log10) = 0.31 + 0.16 (Flight type_Songflights) + 0.29 (Movements_Partial migration) + 0.11 * log10 (900,000 + 1) (Spanish population) will result in a log10 value of 1.41 (that is, 26 individuals; note that the slight difference with table values is due to decimal rounding performed in the model coefficients). The model (R 2 = 0.43; p < 0.001) predicted the real number of collisions for all the species quite accurately ( Figure 2a). Further, we have developed a species collision risk index based on the result of the model including all the species (predicted value, Table 1).
Birds 2023, 4, FOR PEER REVIEW we knew the important predictors, the coefficients and the percent of variance expla by each relevant variable were determined by the linear model, and real values were pared with predicted values to visualize the goodness of fit. The functions were app to the species used to construct the model. Based on the result of the model includin the species, we obtained a species collision risk index.
A two-way ANOVA was also performed with the two more important predi selected by the model. Furthermore, repeated-measures ANOVA was carried out to s the collision seasonality and to relate it to the seasonal movement type of the diffe species.

Results
During the 16 years of the study period, a total of 3130 specimens belonging to bird species died by collisions with wind turbines in Castilla y León. The total numb collisions and their seasonality along with the other data for each species are include Table 1.
First, a model including all the species registered was adjusted and a total of thr the studied variables were included by AIC in the stepwise regression. Flight type the most important variable, explaining 35% of the total variability, followed by seas movement type (4%) and the Spanish population (4%). The function: "Predicted num of collisions (log 10) = 0.31 + 0.28 * Flight type_Active soaring-0.67 * Flight type_Flap + 1.01 * Flight type_Hovering + 1.03 * Flight type_Passive soaring + 0.16 * Flight _Songflights + 0.06 * Movements_Aestival + 0.29 * Movements_Partial migration + 0 Spanish population (log10 + 1)", was applied to the whole dataset of 101 species in o to estimate the predicted number of collisions. For instance, if we consider the Wood (Lullula arborea) case (with Sonflights, Partial migration and a Spanish populatio 900,000 individuals), the Predicted number (log10) = 0.31 + 0.16 (Flight type_Songflig + 0.29 (Movements_Partial migration) + 0.11 * log10 (900,000 + 1) (Spanish population) result in a log10 value of 1.41 (that is, 26 individuals; note that the slight difference table values is due to decimal rounding performed in the model coefficients). The m (R 2 = 0.43; p < 0.001) predicted the real number of collisions for all the species quite a rately ( Figure 2a). Further, we have developed a species collision risk index based on result of the model including all the species (predicted value, Table 1).  and real values (y-axis) for all species included in the model ((a); R 2 = 0.43) and the reduced dataset of 29 species ((b); R 2 = 0.52). Regression line and 95% confidence intervals are shown with a dark line and grey shadow. Baseline levels for each factor (High-speed for the "Flight type" predictor and Sedentary for the "Movements" predictor) have a coefficient of 0. Birds 2023, 4

93
The same three variables were included by AIC in the stepwise regression in the second model adjusted with the reduced dataset. Flight type was the most important variable, explaining 29% of the total variability, followed by seasonal movement type (17%) and the Spanish population (6%). The function: "Predicted number of collisions (log 10) = 0.23 + 0.54 * Flight type_Active soaring + 0.03 * Flight type_Flapping + 0.98 * Flight type_Hovering + 0.32 * Flight type _Songflights + 0.45 * Movements_Aestival + 0.52 * Movements_Partial migration + 0.09 * Spanish population (log10 + 1)", was applied to the 29 species used to construct the model. The model (R 2 = 0.52; p = 0.017) predicted the real number of collisions for the reduced subset of species used to construct the model quite accurately (Figure 2b).
The two-way ANOVA performed with the two most important predictors included in the regression model (flight type and seasonal movement type) for the 29 selected species showed differences between the flight types (F value = 3.177, df = 4, p = 0.039), supporting the importance of this predictor on the probability of collisions with wind turbines. Tukey post-hoc pairwise comparisons did not reveal meaningful differences between specific groups due to the fact that the p-value is adjusted for multiple comparisons. Low sample sizes in the different groups might have also affected the power of the tests. However, the model suggested that species with hovering, song-flight and active soaring flight types were most affected by collisions, while flapping and high-speed species showed a lower number of collisions ( Figure 3).
The same three variables were included by AIC in the stepwise regression in t ond model adjusted with the reduced dataset. Flight type was the most importan ble, explaining 29% of the total variability, followed by seasonal movement type and the Spanish population (6%). The function: "Predicted number of collisions (lo 0.23 + 0.54 * Flight type_Active soaring + 0.03 * Flight type_Flapping + 0.98 * type_Hovering + 0.32 * Flight type _Songflights + 0.45 * Movements_Aestival + Movements_Partial migration + 0.09 * Spanish population (log10 + 1)", was applied 29 species used to construct the model. The model (R 2 = 0.52; p = 0.017) predicted t number of collisions for the reduced subset of species used to construct the mode accurately (Figure 2b).
The two-way ANOVA performed with the two most important predictors in in the regression model (flight type and seasonal movement type) for the 29 select cies showed differences between the flight types (F value = 3.177, df = 4, p = 0.039 porting the importance of this predictor on the probability of collisions with wi bines. Tukey post-hoc pairwise comparisons did not reveal meaningful differen tween specific groups due to the fact that the p-value is adjusted for multiple compa Low sample sizes in the different groups might have also affected the power of th However, the model suggested that species with hovering, song-flight and active s flight types were most affected by collisions, while flapping and high-speed showed a lower number of collisions ( Figure 3). The repeated measures ANOVA for the 29 selected species to study the relati between collision seasonality and seasonal movement type, showed that some s have a greater amount of collision than others (F value = 8.730, df = 3, p < 0.001) an collisions in every season interact to some extent with the seasonal movement t value = 2.431, df = 6, p = 0.033) following predictions (Figure 4). Tukey contrasts re that species showing partial migration have a significant peak of collisions across and autumn. Finally, the results also indicated that sedentary species collide much summer compare to the former and suggested that aestival species have the mos sions during summer and autumn. The repeated measures ANOVA for the 29 selected species to study the relationship between collision seasonality and seasonal movement type, showed that some seasons have a greater amount of collision than others (F value = 8.730, df = 3, p < 0.001) and that collisions in every season interact to some extent with the seasonal movement type (F value = 2.431, df = 6, p = 0.033) following predictions (Figure 4). Tukey contrasts revealed that species showing partial migration have a significant peak of collisions across spring and autumn. Finally, the results also indicated that sedentary species collide much less in summer compare to the former and suggested that aestival species have the most collisions during summer and autumn.  The sex was obtained for only a very few specimens, and was not included model. However, it is interesting to note that for Skylarks, a bird with song-flight, th sexed specimens available were all males (n = 103 collided). In the rest of species wi information, no differences between sexes were detected.

Discussion
The most important species-specific factors influencing the probability of col with wind turbines appear to be the flight type, which influenced the vulnerability o species, and the seasonal movements, which influenced the period of accidents. The dance of each species also influenced the number of collisions but to a lesser exten have developed a species collision risk index based on the result of the model incl all the species and tested it against observed data showing an acceptable fit.
The large dataset of birds available in this study has allowed us to investiga main species-specific factors according to Marques et al. (2014) [29] influencing th for collisions with wind turbines. Since Castilla y León is a large region, this make ideal place to investigate this aspect, which is not possible in studies with a narrowe tial scope.
It is important to note that other species from monospecific families removed the reduced dataset to avoid statistical noise also suffer a high number of collision example Swifts (Apus apus) (n = 123) and Partridges (Alectoris rufa) (n = 43). The height of the swifts searching for aeroplankton possibly influences these results. Th portance of the flight height of the species has also been demonstrated in other work Bat mortality at wind turbines may also be linked to high-altitude feeding on mig insects that accumulate at the turbine towers [39]. In studies carried out with marin cies, it has been shown that, in addition to low maneuverability, fast flight is also a portant risk factor [14]. The reasons for the high vulnerability of Partridges shou studied more deeply but its high population size and habitat preferences may cont to this vulnerability. The sex was obtained for only a very few specimens, and was not included in the model. However, it is interesting to note that for Skylarks, a bird with song-flight, the five sexed specimens available were all males (n = 103 collided). In the rest of species with sex information, no differences between sexes were detected.

Discussion
The most important species-specific factors influencing the probability of colliding with wind turbines appear to be the flight type, which influenced the vulnerability of each species, and the seasonal movements, which influenced the period of accidents. The abundance of each species also influenced the number of collisions but to a lesser extent. We have developed a species collision risk index based on the result of the model including all the species and tested it against observed data showing an acceptable fit.
The large dataset of birds available in this study has allowed us to investigate the main species-specific factors according to Marques et al. (2014) [29] influencing the risk for collisions with wind turbines. Since Castilla y León is a large region, this makes it an ideal place to investigate this aspect, which is not possible in studies with a narrower spatial scope.
It is important to note that other species from monospecific families removed from the reduced dataset to avoid statistical noise also suffer a high number of collisions, for example Swifts (Apus apus) (n = 123) and Partridges (Alectoris rufa) (n = 43). The flight height of the swifts searching for aeroplankton possibly influences these results. The importance of the flight height of the species has also been demonstrated in other works [35]. Bat mortality at wind turbines may also be linked to high-altitude feeding on migrating insects that accumulate at the turbine towers [39]. In studies carried out with marine species, it has been shown that, in addition to low maneuverability, fast flight is also an important risk factor [14]. The reasons for the high vulnerability of Partridges should be studied more deeply but its high population size and habitat preferences may contribute to this vulnerability.

Flight Type
According to the results of this study, flight type is the most important factor in determining the probability of a bird to collide. For instance, species with hovering, songflights and active soaring flight types are more susceptible to collisions with wind farms, although differences between species depending on their seasonal movements also could be found. For partially migrating species, the three flight types mentioned showed a higher risk of collision with wind turbines. For aestival species, active soaring is the flight type with a higher risk.
Flight type, when associated with hunting and foraging strategies, seems to play an important role in collision risk with wind turbines [29]. Certain flight types, like soaring and hovering, are more suitable for looking for food. Hunting raptors, when they are on the lookout for prey, do not perceive the risk and, hence, collide [29]. On the other hand, the song-flight type is done to attract females. According to the results of this study, when the bird focuses his attention on foraging or partner attraction, collisions are more likely to occur.

Seasonality
The patterns of risky flights have a temporal component and the deaths are concentrated in specific seasons [40]. In this study, considering the evolution of impacts through the year, it is remarkable that spring and autumn are the seasons with the most collisions. We have shown that aestival species have the most collisions during summer and autumn, while partial migration species have the peak of collisions across spring and autumn. This result is expected as a consequence of the movements performed by many migratory species, and is also found in other studies [28,[41][42][43]. The location of wind farms in migration areas may cause higher mortality rates [44,45], and passerine species appear to be at the greatest risk of colliding during spring and autumn migrations [28]. The rapid expansion of wind farms in northern Spain urgently requires preventive measures to protect migratory species in particular [45].
It has been suggested that, due to their knowledge of the terrain, resident birds would be less prone to collision with structures [7]. However, comparable data on collision rates for resident birds is lacking [29]. According to the results of this study, sedentary species collide mostly in spring and winter, possibly because they make more movements in these periods, although to a lesser degree than migratory species. In some studies, resident populations have more collisions than migratory species [41], but this is probably due to the bias produced by a large number of collisions of a single species. For instance, the Griffon Vulture (although in this study is considered a partial migratory species and is removed from the reduced dataset) accounted for most of the collisions.

Species and Group
Raptors are considered the most vulnerable bird species for collisions in wind farms [28,36,41,46]. However, according to the result of this study, the bird group (passerines, raptors, etc.) does not play an important role regarding the collision risk with wind turbines.
In this work, the Griffon Vulture was the most affected species. The high susceptibility of Griffon Vultures to wind farms has already been described in other works [41,46], and their vulnerability may be a result of its flight pattern, the search for food in flight, and its large size and reduced maneuverability [44]. Therefore, this species may be causing a bias in the studies done so far by increasing the mortality of raptors by a single species. Griffon Vulture collisions occur throughout the year, but especially in spring and autumn. The increase in collisions in these periods can be explained by the increase of movements during the breeding season, a possible increased need to search for food, and immature bird dispersal in autumn.
Another species highly susceptible to collisions are the Kestrels, and this happens with both American (Falco sparverius) [28], and European species, Common Kestrel (Falco tinnunculus) [40]. Its hovering flight and concentration on prey probably influence its high mortality [29], as already discussed.
As we have pointed out, the models in this study quite accurately predicted the real number of collisions for the total number of species, and for the reduced subset of species studied. Analyses performed in this study show essential factors that should be taken into account when developing environmental impact studies associated with wind farms. Specific measures towards the identification of more vulnerable species, particularly those with higher predicted values, could help in choosing the location for wind turbines or restrict their activity the months of the year in which the largest amount of collisions occur. These actions could decrease the number of collisions for species most affected by wind farms.
Taking into account the IUCN Red List of Threatened Species, and with the results predicted by the model, the endangered species most at risk for collisions are: Cinereous Vulture (Aegypius monachus), Red Kite (Milvus milvus) (both near threatened according to IUCN), and Egyptian Vulture (Neophron percnopterus) (endangered).

Abundance and Size
Another factor traditionally considered important for the number of collisions, is the abundance of the species, although, as has been verified in this study, flight-related factors are more important than abundance [29]. Collision fatality of raptors in wind farms does not depend on their abundance [46]. In our models, the number of couples from the Spanish population explained 4% of the variability whereas the European population was not relevant. This factor can help explain the high mortality for the Griffon Vulture and the Common Kestrel, since these species represent large populations with a complete distribution on the Iberian Peninsula. Additionally, in Castilla y León, there are colonies of Griffon Vulture in several of the canyons on the main rivers, and the presence of nearby colonies is a risk factor as has been shown with Seagulls [41] and Lesser Kestrels [47].
According to the results of the model in this study, size does not seem to be an important explanatory factor on bird collisions and these results go against the knowledge that exists so far [29]. As discussed before, it could be that the griffon vulture has produced a bias due to the large number of collisions and its big size.

Sex
A sex bias on wind collision has already been found in previous works [29]. In a study from Portugal, the 22 Skylark carcasses collected revealed a higher incidence of adult males (90%), suggesting the occurrence of a sex differential mortality [48]. The carcasses were mainly found between April and May, within the breeding season, when Skylark males perform their characteristic song-flights. In this study, the Skylark was a species that collided frequently (n = 103), and had one of the highest predicted values in the model (32.1). Also, from five sexed specimens in this study, all were males. On the other hand, sex-biased collision mortality in Common Terns (Sterna hirundo) reflects differences in foraging frequency between males and females during egg-laying and incubation [49].

Modelling Approaches to Quantify the Theoretical Risk of Collision in Wind Farms and Ways to Avoid it
The results of this study can help in future modelling approaches by providing improvements in knowledge of factors influencing collision and quantifying the theoretical risk of collision in wind farms depending on the species existing in the area and their predicted values linked to flight types and seasonal movements. For many years, modelling approaches to quantify the theoretical risk for collision have been developed [31,35]. The Band model, a valuable tool for impact assessment, is widespread in the UK [50]. However, birds' behaviour is poorly understood and this can have a large effect on the model outputs [35].
There are good technical documents that could help for the adoption of measures to avoid damage to wildlife at wind farms [5,29,[51][52][53]. The location of wind farms is still the most effective measure to avoid bird fatalities [29,[54][55][56]. Unlike other non-natural causes of bird mortality, wind farm fatalities can be lowered by powering down or removing risky turbines and by placing them outside critical areas for endangered birds [6]. The spatial distribution and aggregation of some threatened species should be used as criteria for environmental planning [23]. The use of selective stopping techniques for turbines with the highest mortality rates can help to mitigate the impacts of wind farms on birds, with a minimal effect on energy production [22].

Limitations of this Study
In the human dimension of finding, catching and transporting dead birds to a rescue centre there are some assumptions necessary to let the database reflect true turbine mortality. The dead-injured ratio, the ability of an injured bird to draw attention to people, the abundance of people in the neighborhood of windfarms and the abundance of local predators might negatively affect injured numbers and biases from such data can affect the results in different ways [57]. These problems can happen especially for the small species that may have high numbers not reported or not found. On the other hand, species may be more prone to occur in some areas than others, so site-specific (orientation, location etc.) and wind farms features, are also important explanatory factors when investigating the causes of collisions with wind turbines that have not been analysed in this work [58].
The collisions with wind farms may also be influenced by behaviour associated with a specific age [29]. In this work, it has not been possible to verify the importance of age in collisions due to the absence of this data, but we predict that young birds are more vulnerable due to lack of experience and knowledge of the territory.
The dataset used for this work was a small part of the animals that die in the field, since many of them may not be found [59], or due to detection bias and loss of carcasses from scavenging [60]. On the other hand, it is possible that some specimens have arrived at the rescue centres a few weeks after the collision has occurred, although such cases should be anecdotal and should not have influenced the seasonally patterns observed. Finally, the population data to estimate the abundance of the different species in this study is the number of breeding couples [40]. However, there are species that winter but do not breed in the Iberian Peninsula, and in these cases its reproductive population is zero couples.

Conclusions
Flight type and seasonal movements explain the accidents in wind farms in a relevant way. The flight type influences the vulnerability of each species, and the seasonal movements influence the period of accidents. The abundance of each species also influences the number of collisions but to a lesser extent. Species with hovering, song-flights and active soaring flights are more susceptible to collisions with wind farms, although differences between species could also be found based on their seasonal movements. According to the results of this study, for the flight types in which the bird focuses his attention to foraging or partner attraction, collisions are more likely to occur. The results of this study can help in modelling the theoretical risk of collision with wind turbines, depending on the species existing in the area and their predicted values of vulnerability, which is linked to flight types and seasonal movements, in a similar way to that already proposed in previous studies [61].