Telemetry reveals strong effects of offshore wind farms on behaviour and habitat use of common guillemots (Uria aalge) during the breeding season

Seabirds have increasingly encountered offshore wind farms (OWFs) in European waters in the past 10 years, resulting in potential conflicts with offshore foraging areas. During the breeding season, seabirds are restricted in their choice of foraging habitat and are under increased pressure to find enough prey to raise their offspring. However, information on the individual reactions of seabirds towards OWFs during the breeding season is lacking. Three OWFs located 23–35 km north of the island of Helgoland have operated since October 2015. We studied their possible effects on locally breeding common guillemots (Uria aalge) using GPS tracking. GPS tags were deployed on 12 breeding guillemots from Helgoland for 8–26 days during 2016–2017. Most individuals avoided the OWFs, but one individual in each year briefly entered the OWFs on two or three occasions. Using a point process model, we revealed a 63% reduction in the resource selection of the OWF areas compared with the surroundings (lower confidence interval (CI) = 79% reduction, upper CI = 36% reduction). Furthermore, OWF avoidance was increased to 75% when the turbine blades were rotating (lower CI = 93% reduction, upper CI = 11% reduction). Guillemots mainly approached the OWFs from their eastern edge when resting or diving, and rarely approached the areas when commuting. These results provide a detailed description of guillemot reactions to OWFs during the breeding season, and the first comprehensive analysis of OWF effects on this species based on telemetry data. The strong avoidance effect for guillemots during the breeding season indicates the need to consider the presence of OWFs when interpreting future trends in the abundance and breeding success of this species.


Introduction
Seabirds increasingly encounter offshore wind farms (OWFs) in European waters (Perveen et al. 2014) as numerous OWFs were built during the last decade (Perveen et al. 2014;4COff-shoreWind 2020;BSH 2020). Further areas will be occupied in the near future as the construction of already approved OWFs and the consenting of additional areas for OWFs is still in progress (4COffshoreWind 2020; BSH 2020). A possible conflict with this development arises from the fact that seabirds, as well as other top predators such as marine mammals, thoroughly depend on offshore areas for foraging, resting and migration (Schreiber and Burger 2001;Wilson and Mittermeier 2014). The possible effect of an offshore wind farm largely depends on the local environmental conditions, the species occurring in and depending on the respective marine area (Drewitt and Langston 2006;Vanermen and Stienen 2019) as well as on the species' annual life cycle stage e.g. with regards to breeding, migrating, and moulting (Masden et al. 2009;Busch and Garthe 2016;. Depending Responsible Editor: T. A. Clay.
Reviewed by undisclosed experts.

Electronic supplementary material
The online version of this article (https ://doi.org/10.1007/s0022 7-020-03735 -5) contains supplementary material, which is available to authorized users. on the species different behavioural reactions towards OWFs are known, ranging from complete avoidance to attraction (Garthe and Hüppop 2004;Furness et al. 2013;Dierschke et al. 2016). When avoiding OWFs species can possibly lose a former valuable (foraging) habitat or have to extend their travel route to reach their foraging areas (Fox et al. 2006) which could lead to increased energy consumption (Masden et al. 2010). When attracted to OWFs seabirds are prone to collide with the rotors or turbines which leads to an increased mortality (Drewitt and Langton 2006;Fox et al. 2006).
Common guillemots (Uria aalge) are the most abundant seabird species in the North Sea and are thus prone to experience interference with OWF at various locations. Varying responses of guillemots to OWFs are reported (overview in Dierschke et al. 2016): In Belgium, The Netherlands and Germany ship-based surveys at several OWFs revealed that guillemots largely avoided the OWF areas, with a reduced abundance of 71-75% in the OWF compared with the surrounding areas (Leopold et al. 2013;Vanermen et al. 2015;2016;Welcker and Nehls 2016). By contrast, other studies which investigated OWF effects on seabirds in the UK and were also conducted by means of ship-based surveys, either detected no changes in abundance or found that guillemots were attracted by OWFs (PMSS 2007;Vallejo et al. 2017). The aforementioned studies estimated the overall year-round avoidance reaction of guillemots towards OWFs. However, depending on the species' annual life cycle stage, e.g. the breeding compared to the non-breeding seasons, the possible effect of an OWF on seabirds can potentially be altered (Masden et al. 2009;Furness et al. 2013;Busch and Garthe 2016;. During the breeding season seabirds are restricted in their choice of foraging habitat and under increased pressure to find enough prey to raise their offspring (Orians and Pearson 1979). But few studies have explicitly investigated the reaction of seabirds towards OWFs during the breeding season (Masden et al. 2010;Thaxter et al. 2015;. The interaction of breeding lesser black-backed gulls (Larus fuscus) with OWFs was shown to be highly variable between individuals (Thaxter et al. 2015; and between years (Thaxter et al. 2015). A modelling approach revealed that the additional energetic costs due to avoidance of OWFs during the breeding season were highly species-specific and strongly increased for common guillemots (Masden et al. 2010). Due to their comparatively short foraging ranges during the breeding season (medium range between 7.8 and 37.8 km, Enstipp et al. 2006;Thaxter et al. 2012;Evans et al. 2013) and their flight characteristics (low manoeuvrability and comparatively high energetic costs of flight, Enstipp et al. 2006;Masden et al. 2010) guillemots experience large additional costs when travelling longer distances while avoiding OWFs (Masden et al. 2010). Thus, they are especially susceptible to the presence of OWFs close to their breeding colonies. However, information on the effect of OWFs on breeding guillemots, as well as on their individual movements and behaviour in and around the OWFs is largely lacking.
Studies on the effects of OWFs on seabirds have mainly applied visual or digital surveys (reviewed in Dierschke et al. 2016;Mendel et al. 2019) or radar assessments (reviewed in Dierschke et al. 2016;Fijn et al. 2015). However, behaviours in and around OWFs, and individual reactions towards the OWFs cannot be studied in detail by these methods. GPStracking methods are increasingly applied to investigate the behaviour of seabird species (Garthe et al. 2007;Ponchon et al. 2017;Thaxter et al. 2018). These methods allow researchers to follow individual birds for weeks or months. They thus provide powerful tools for studying both fineand large-scale habitat use and the behaviour of marine top predators (Garthe et al. 2007;Ponchon et al. 2017;Thaxter et al. 2018), and for monitoring alterations in these factors in relation to environmental change or anthropogenic activities (Garthe et al. 2011(Garthe et al. , 2017bParedes et al. 2014). Few studies to date have used GPS tags to investigate the effects of OWF on seabirds (Wade et al. 2014;Thaxter et al. 2015;Garthe et al. 2017a;2017b).
The only guillemot colony in the south-eastern North Sea is located on the island of Helgoland, 48 km off the German coast (Fig. 1). The breeding numbers are relatively small, comprising 3178 breeding pairs in 2017 (Dierschke et al. 2011;, and the number of breeding guillemots has remained relatively constant for the last several years (2000( , Dierschke et al. 2011.
However, three wind farms, including a total of 208 turbines and covering an area of 105 km 2 , are currently operating only 23-35 km north of Helgoland, making it necessary to assess their possible effects on breeding guillemots.
In this study, we used GPS-tracking methods to investigate if and how guillemots breeding on Helgoland interact with OWFs close to their colony; e.g. are they displaced or attracted by the OWFs? We also investigated individual differences in their interactions with the OWFs; i.e. do individuals move through the OWFs or approach them? Additionally, we examined how they behaved in the vicinity of the OWFs, e.g. if they forage close to the OWFs.

Fieldwork and data collection
Guillemots were caught on the island of Helgoland (54°11′ N, 7°55′ E) in the south-eastern North Sea during the breeding seasons in 2016 and 2017, using a noose pole. Thirteen breeding guillemots were equipped with GPS devices, of which 12 successfully transmitted data (seven in 2016, five in 2017). Each bird received a 'Uria' GPS-logger (13 g, Ecotone Telemetry, Poland) which were attached to the lower back using TESA® tape (Beiersdorf AG GmbH, Hamburg, Germany). Six birds in 2016 also received a G5 temperature-depth recorder (TDR, 2.7 g, Cefas Technology, Suffolk, UK) attached either to the lower part of the tail feathers or to a ring mounted on the bird's leg. Three of the TDR loggers were retrieved. The attached devices represented 1.8% (mean value) of the guillemot's body mass (mean: 902 g), which was below the recommended threshold of 3% (Phillips et al. 2003;Vandenabeele et al. 2012). The mean handling time was 14 min, and the birds were then released in the close vicinity of the breeding site. The bird's eggs were observed during the handling period to prevent egg predation. Individuals were then either re-caught after 2-3 weeks and the devices removed, or the devices fell off when the birds moulted. All the GPS devices recorded the date, time, and position (latitude, longitude) with a scheduled sampling interval of 10 min, and recorded if the logger was submerged or not (wet/dry sensor) every second. Depending on the battery voltage and whether the bird was submerged during a scheduled GPS recording, the GPS interval could be extended to 11-20 min (19% of the dataset) or 21-30 min (3%). The devices transmitted the data via a UHF connection to a base station. The TDR devices recorded temperature at an interval of 20 s, and depth every 1 s. The TDR devices had to be retrieved for data download.

Data analysis
To remove non-foraging activities such as resting on the sea surface, which often occurs close to the colony, trips with a duration > 12 min and ≥ 500 m distance from the breeding site were classified as foraging trips. Classification was carried out using an R code provided by Lascelles et al. (2015). Thirty-two trips were excluded from the subsequent analysis as the trip start or end times were unknown due to gaps in the GPS data. Trip statistics, i.e. trip duration (h), trip maximum distance to Helgoland, and total distance (km) were also calculated for each individual foraging trip using the R code provided by Lascelles et al. (2015).

Behavioural classification
To illustrate how the guillemots behaved in the vicinity of the OWFs, we grouped the dataset in three behavioural classes: foraging, resting and travelling. Based on the information from the wet/dry sensors, it was possible to determine when and for how long the birds were submerged. As the time of the start and end of each dive were known with an accuracy of 1 s, we interpolated the corresponding GPS data to an interval of 1 s to retrieve the most likely locations of potential dives. Here we only used GPS data which belonged to a foraging trip and only interpolated GPS intervals of maximum 10 min. Data on 6145 dives were available from the TDR data for three individuals, and detailed dive characteristics were then analysed to define the threshold duration to classify a dive recorded by the GPS device as foraging dive. Duration of dives recorded by the TDR devices started from 4 s, dives with a duration of ≥ 10 s recorded a mean water depth of ≥ 1 m during the bottom phase. Following Thaxter et al. (2010), dives with depths > 1 m can be interpreted as foraging dives, while shallower dives often occur during bathing or other activities not associated with foraging. As we aimed to exclude such events in our dataset, only dives lasting ≥ 10 s were classified as foraging dives. Furthermore, based on inspection of the speed distribution in the GPS data a bimodal pattern was revealed, with 10 m/s approximating the break between resting and travelling. Thus, birds moving at ≤ 10 m/s were classified as resting and birds moving at > 10 m/s were classified as travelling (see Fig. S3 in Online Resource 2). To visualise the areas in which the different behavioural states were shown, kernel densities of the positions assigned to each category were created in ArcGIS using the ArcMET tool (version 10.2.2v3; Wall 2014).

Point process models
In addition to the more descriptive approaches, we aimed to quantify the effect of the OWFs statistically. We thus investigated if the guillemot's resource selection of the OWF area was reduced in comparison with the areas outside the wind farm. Statistical analysis of telemetry data investigating resource selection is often challenging, and various modelling strategies have previously been developed and discussed (Hooten et al. 2017). Popular approaches include (integrated) step selection functions (Thurfjell et al. 2014;Avgar et al. 2016) and point process approaches (Johnson et al. 2013;Renner et al. 2015). Both approaches use a number of contrasting points (e.g., 'dummy points', 'pseudo-absences', or 'available steps') in addition to true tracking locations, making it possible to compare selected versus available resources. Methods using contrasting points tend to produce better results than techniques using presence points alone (Brotons et al. 2004;Elith et al. 2006;Barbet-Massin et al. 2012).
In the following analysis, we used and extended the spatio-temporal point process models (PPMs) presented by Renner et al. (2015), which naturally and automatically resolve many of the questions and pitfalls associated with alternative approaches (Warton and Shepherd 2010;Warton and Aarts 2013;Renner et al. 2015;Hooten et al. 2017). For example, the role and number of dummy points is not ad hoc but can be deduced mathematically by the efficient estimation of an integral as a part of the PPM likelihood (Warton and Shepherd 2010;Warton and Aarts 2013). Additionally, PPMs represent a generalisation of many other frequently used methods (Johnson et al. 2013;Warton and Shepherd 2010;Aarts et al. 2012). Finally, the PPM likelihood can be approximated by a mathematical method using standard generalised linear mixed modellingregression software (Johnson et al. 2013;Renner et al. 2015), ensuring flexible and individual implementation. Details of the modifications of the PPM compared with the spatio-temporal PPMs are presented by Johnson et al. (2013) (see Online Resource 1).

Preparation of covariates
ArcGIS (version 10.3; Environmental System Research Institute 2016) was used to calculate the means of the spatial covariates (owf_yn, dist_coast, dist_Helgoland, depth, slope; Table 1) for a grid with a spatial resolution of 200 × 200 m. The variables used for modelling are described in Table 1.

Model selection
We determined if habitat use by guillemots was affected by the presence of OWFs by applying a generalised additive mixed model (GAMM)-PPM to a dataset consisting of 9235 raw data points collected in 2016 and 2017. When applying the GAMM-PPM to the raw tracking data 54,715 dummy points were created.
The optimal model regarding the set of fixed-effect predictors was selected by comparing 17 different models (see Online Resource 2, Table S2) based on the Akaike Information Criterion (AIC; Akaike 1973). Models including interaction terms of the variables dist_coast and slope were not selected through AIC analysis. Inspection of the results of the best models showed that all models revealed similar patterns in the data, indicating that our main results were robust across different models. We first inspected a basic model including only autocorrelation terms, random effects, and the variables dist_Helgoland and OWFyn, which was the main focus of the analysis: where ß is the intercept and te() a tensor-product regression spline considering spatial and directional autocorrelation, where the optimal number of knots has been estimated via generalised cross-validation. Especially, log d s is the logarithm of the Euclidean step-distance, whereas angle is the turning angle. Trip id and bird id were included as random (1) Z ∼ + te log d s , angle, k = c(5, 5) + s trip id , bs = � re � + s bird id , bs = � re � + s dist Helgoland + OWFyn effects, indicated by the term s(…,bs = re). To approximate the PPM likelihood based on standard GAMM software, a weighted regression Poisson model has been fitted, using regression weights W and observations Z, where Z = 1/W has been defined for tracking points, and Z = 0 for dummy points. Especially, W are appropriate quadrature weights based on the 2D rectangle rule (for more technical details see for example Johnson et al. 2013).
We subsequently added other variables to the basic model to identify the best model for our data. We restricted the maximum number of variables added to the basic model to three, to keep the models interpretable.
The best model was then selected via the AIC (during model selection the lowest AIC value was favoured). This model enabled us to detect any reduction in resource selection inside compared with outside the OWF if the turbine blades were rotating (i.e. OWFyn = inside and OWF_ onoff = rotating, see Eq. (2); interaction term shown in bold). Notably, when interactions are included in a model, the variables included in the interaction term cannot be interpreted singularly (Field et al. 2012), but they were nevertheless included to improve the model: The second-best model was also selected because it allowed us to determine if resource selection of the OWF area was reduced in comparison with areas outside the OWFs, without analysing the combined effect of rotating turbines and OWF presence:

Model validation, numerical realisation, and software
PPM model-validation plots for the final GAMM-PPM were generated based on PPM-Pearson residuals (Baddeley and Turner 2005; . All statistical analyses were performed using the free statistical software R ( R Core Team 2017). Spatial statistics were performed using spatstat , dummy-point meshes and trapezoid rule-based quadrature weights were created using mvQuad (Weiser 2016), and GAMM and GAM fits were performed using the package mgcv (Wood 2006). All the codes were programmed to allow the main parts of the code to be run using parallel computing, using Z ∼ + te log d s , angle, k = c(5, 5) + s trip id , bs = � re � + s bird id , bs = � re � + OWFyn + s dist Helgoland + depth the parallel package and the bam() function from the mgcv package.

Overview of foraging trips
We recorded a total of 204 individual foraging trips by 12 guillemots in 2016 and 2017 (Table 2, Online Resource  2 Table S1). The mean values of the duration, and the maximum and total distances of the foraging trips did not differ between the 2 years (Table 2, Online Resource 2 Fig. S1). While the same area was generally used in both years (Fig. 1), individuals in 2017 used areas located more to the west of the colony.

Foraging trips
Most individuals completely avoided the OWFs north of Helgoland (Fig. 1). However, some individuals approached them closely, especially in the east or south. One out of seven tagged individuals in 2016 and one out of five tagged individuals in 2017 entered the OWFs on two to three occasions at night or in the evening, and stayed inside the OWFs for between 20 min and 2 h 11 min (mean = 1 h 1 min; Table 3, Fig. 1). For both these individuals, diving activity in the OWF area was recorded on one of the occasions (Table 3), but no dives were recorded or the positions of the dives were unavailable for the remaining occasions due to gaps in the GPS data. Only one individual entering the OWFs is shown in Fig. 1 because some parts of the data for the second bird did not allow trip classification. The positions of the latter individual are shown in Online Resource 2, Fig. S2.

Behaviour
Differentiating among resting, foraging, and travelling behaviours showed, that the individuals mainly used areas north and north-west of Helgoland (i.e. south-west or east of the OWFs) to commute to and from the colony (Fig. 2a), and areas to the west, south-west, and north of Helgoland for resting and diving (Fig. 2b, c). Guillemots rarely approached the OWFs when flying (Fig. 2d), and mainly approached from the east while resting, and also when diving (Fig. 2e, f).

Model results
All the models revealed significantly reduced resource selection of the OWF area compared with the surrounding areas (Tables 4,5,6). The second-best model revealed that the selection of the OWF area was reduced by 63% when the rotation of the turbine blades was not considered (Table 5; variable 'owfyn_inside ', estimate = − 0.988, p < 0.001, response = exp(estimate) = 63% reduced selection inside the OWF compared to outside, lower confidence interval (CI) = 79% reduction, upper CI = 36% reduction). The best model revealed a by 75% reduced selection of the OWF area compared with outside if the blades were rotating (

Discussion
To the best of our knowledge, this study provides the first evidence for the impact of OWFs situated in close proximity to a common guillemot breeding colony based on GPStracking data. GPS tracking combined with PPM analysis showed that guillemots almost completely avoided the OWFs north of their breeding colony, with a 63% reduction in selection of the OWF areas compared with the surrounding areas. Furthermore, OWF avoidance was increased (75% reduction) when the turbine blades were rotating. Guillemots mainly approached the OWFs at their eastern edge when resting or diving, and rarely approached them when commuting. However, one individual in each study year briefly entered the OWF on two to three occasions.
As central-place foragers, breeding seabirds have to balance their energy intake and trip distance (Chaurand and Weimerskirch 1994;Suryan et al. 2006;Bertrand et al. 2012). Displacement during the breeding season could thus increase foraging trip distance and energy expenditure, thereby affecting adult condition or survival (Masden et al. 2010) and reproductive success (Langton et al. 2014). Breeding guillemots are especially susceptible to an increase in foraging trip distance (Masden et al. 2010), thus displacement by OWFs close to their breeding colony can strongly increase their energy expenditure (Masden et al. 2010). The strong OWF avoidance detected for guillemots breeding on Helgoland, together with their comparatively small foraging distance, shows that the available foraging habitat close to the colony was greatly reduced by the OWFs. Guillemots breeding on Helgoland thus have to compensate for this lost habitat within their accessible distance. Even though no pre-construction tracking data was available for the current study, the guillemot distribution in and around the future OWFs was analysed by a recent study using aerial and ship-based surveys (Mendel et al. 2018). The OWFs north of Helgoland were not among the most intensively used areas by guillemots before their construction, but moderately high densities were still found during the breeding season (Mendel et al. 2018). The same study supports our findings that guillemots approached the OWFs closely, but generally did not enter them. Since studies on different potential prey species revealed that abundances are either not affected or positively influenced by OWFs (cf., below), the observed lack of guillemots in OWF areas is most probably not a secondary effect of a lack of prey, but can, on the contrary, be seen as a conservative estimate of the scaring effect. Other studies have revealed similar OWF avoidance in guillemots, but most focused on the entire yearly cycle and did not provide detailed information on the behaviour of breeding birds (Vanermen et al. 2015;Welcker and Nehls 2016;Vallejo et al. 2017).
The use of areas close to the OWFs for resting and diving indicates that the guillemots foraged close to the turbines on some occasions. This could be related to the so called 'reef effect' (Lindeboom et al. 2011;Leopold et al. 2013;Vanermen et al. 2015), which reflects the increased fish diversity and abundance associated with the additional benthic structures of the turbines (e.g. de Mesel et al. 2015;Stenberg et al. 2015;Vandendriessche et al. 2015). Although no information on the distribution of guillemot prey species in and around the OWFs north of Helgoland was available for the current study the influence of OWF presence on different fish species is investigated in other studies. For guillemot prey species such as sandeel (Ammodytes spec., Hyperoplus spec.), herring (Clupea harengus), sprat (Sprattus sprattus), whiting (Merlangius merlangus) or atlantic cod (Gadus morhua, Halley et al. 1995;Mitchell et al. 2004;Sonntag and Hüppop 2005) differing reactions towards OWFs are reported. Sandeel has been shown to increase on the short term but no effect was detected on the long term (Lindeboom et al. 2011;Van Deurs et al. 2012) and no increase in abundance was shown for herring or sprat (Krägefsky 2014). Whiting and Atlantic cod, however, increased in abundance in OWF areas (Reubens et al. 2013;Stenberg et al. 2015). It was furthermore shown, that also pelagic species can be attracted to foundations (here shown for Atlantic horse mackerel, Trachurus trachurus, Schröder et al. 2013). However, guillemots did not use the areas close to the OWFs very intensively but rather concentrated their foraging efforts in the west of the study area at distances of 12-18 km from the OWFs. The increased avoidance of the OWFs when the rotor blades were turning indicated a stronger reaction to moving objects, which was suggested as a possible reason for OWF avoidance by birds (Dierschke et al. 2016). However, so far only a few other studies could investigate if the movements of the rotor blades can have an effect on birds entering the OWFs (reviewed in Dierschke et al. 2016). Based on radar assessments at the OWF Egmond aan Zee, e.g., the number of bird tracks crossing the wind farm area was two to three times higher when the nearest turbine was off . The increased wind speeds while rotor blades were turning could also have influenced the guillemot's behaviour close to the OWFs. As guillemots might have less control of their flight behaviour during strong winds, their avoidance of the OWFs might increase to reduce the risk of collision. A similar behaviour was observed in wind farms on land where lower numbers of birds passed through turbine rows during higher wind speeds (Smallwood et al. 2009). At the OWF Egmond aan Zee the average number of birds migrating through the OWF also decreased with increasing wind speeds ). However, we found strong OWF avoidance across the entire dataset, irrespective of blade rotation, suggesting that the rotating blades were not the main cause of the avoidance reaction. For example, the presence of very high vertical structures which are unusual for the open marine areas and the increased vessel activity in and around the OWFs could also cause avoidance reactions (Vanermen et al. 2015;Dierschke et al. 2016;Mendel et al. 2019). Nonetheless, it could be speculated that if measures were needed to reduce the effect, deactivation of the turbines could slightly reduce the disturbance. Busch and Garthe (2016) assessed the effects of potential displacement due to the OWFs north of Helgoland on the local guillemot colony on the basis of potential biological removal rates (PBR). The PBR indicates the number of losses in addition to the natural mortality, which can be sustained each year by a certain population (Dillingham and Fletcher 2008). Transferring the findings of Busch and Garthe (2016) to our study indicated that the OWFs north of Helgoland accounted for 30% of the PBR estimated for the German guillemot population in the breeding season. This represents a substantial proportion of the PBR, bearing in mind that the effects of other human activities (e.g. shipping, fisheries, etc.) were not included in these calculations. However, it has to be kept in mind, that Busch and Garthe (2016) applied displacement values derived from changes in the species abundance in OWF areas for estimating the number of affected individuals. Thus, the changes in resource selection estimated from tracking data in our study potentially cannot directly be applied to their approach.
When interpreting the here presented findings, the duration of the study has to be considered. Studies conducted over several consecutive years will reveal if the observed pattern is consistent over time. Furthermore, information on the distribution of guillemot prey species in and around the OWFs would be useful to better understand their choice of foraging habitat.
The strong avoidance effect found in this study during the breeding season shows that, although the size of the guillemot population on Helgoland is currently quite stable (Dierschke et al. 2018), the presence of the OWFs needs to be considered when evaluating the conservation and health statuses of common guillemots at their sole breeding site in German North Sea waters, and likely also at other sites within their distribution area.
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creat iveco mmons .org/licen ses/by/4.0/.