Seasonal activity levels of a farm-island population of striated caracaras (Phalcoboenus australis) in the Falkland Islands

Animals need adaptive strategies to cope with seasonal changes in prey availability to survive and reproduce, which can include migrating, prey-switching, or reducing metabolic needs. Human settlements can disrupt spatiotemporal patterning in resource availability, which can affect animals’ foraging success, particularly for juveniles who are behaviorally developing and learning efficient foraging skills. Here, we investigate behavioral responses of juvenile striated caracaras, a near-threatened scavenging falconid, to seasonally migratory seabird prey and a farm settlement on Saunders Island, Falklands. We deployed accelerometer–GPS dataloggers (n = 27) to measure seasonal differences in daily and hourly vectorial dynamic body acceleration (VeDBA; an activity index measured in gravitational g) and space use, and investigated seasonal variation in body mass of the tagged subset and an additional 65 caracaras. Juvenile caracaras were overall similarly active in winter and summer. However, during winter, caracaras made the most of limited daylight by increasing average daytime activity (winter males: 0.16 ± 0.03 g, summer males: 0.09 ± 0.01 g, winter females: 0.12 ± 0.02 g, summer females: 0.08 ± 0.01 g). During winter, both sexes increased the percentage of daylight spent in high activity (winter males: 35 ± 5%, summer males: 21 ± 3%, winter females: 25 ± 6%, summer females: 16 ± 3%, p < 0.001) and ranged nearly 4 times farther (95% kernel density estimate winter: 2.36 ± 0.96 km2, summer: 0.61 ± 0.20 km2; p < 0.001). Furthermore, on a daily scale, males were 21% more active than females year-round (24-h average VeDBA: males 0.07 ± 0.01 g, females 0.06 ± 0.01 g; p < 0.01). We did not observe a significant seasonal difference in mass. That caracaras’ daily activity and body mass did not vary between seasons suggests that wintering birds on Saunders are meeting resource requirements despite the absence of seasonally migratory prey. We hypothesize that human subsidies may mitigate the effect of seasonal food limitations. Further research should include studies on seasonal energetics to improve our understanding of baseline body condition, and comparative studies on other islands and including adults to understand the importance of human subsidies.

While the effects of seasonal food availability have been investigated at the population level for several species [14,[16][17][18][19][20], less is known about the impact on individual life-history strategies [14], particularly during the juvenile stage when individuals are still behaviorally developing, learning judicious time-energy management, efficient foraging skills, and competitive strategies [21,22]. Across many mammalian and bird species, juvenile mortality can be high [21,23,24], thus making it an important age class to study as impacts extend to population level effects. Here, we report on seasonal juvenile behavioral patterns of an endemic island predator in a system with fluctuating seabird prey and a relatively recent introduction of humans (i.e., since the mid-nineteenth century).
Striated caracaras (Phalcoboenus australis) are nearthreatened falconids that inhabit the southern coasts of Patagonia and the outer Falkland Islands [25][26][27]. They are long-lived [Amy Wallace pers. comm] with extended periods of adolescence [28] and exhibit reverse sexual size dimorphism [KH unpubl. data]. During austral summer, caracaras primarily feed on eggs, chicks, and dead or dying adults of breeding colonial seabirds, the afterbirth, feces, and carrion of resident pinnipeds, and invertebrates including limpets and kelp fly larvae [17,[28][29][30]. During austral winter, when most seabirds migrate offshore, caracaras' diets include native upland geese (Chloephaga picta), the feces and carrion of resident southern fur seals (Arctocephalus australis), southern sea lions (Mirounga leonina) and gentoo penguins (Pygoscelis papua), as well as beetles and subsidies provided by farms [17,28,31]. Since human occupation began in the Falklands in the mid-nineteenth century, many islands' vegetation has changed dramatically due to farming [32], resident pinniped populations suffered a population decline from the sealing industry [33], and caracaras were persecuted to the point of near extinction [25,34]. While the caracaras in the Falklands are now protected, a previous study posited that population recovery may be limited by winter food stress [31].
Animals must maintain a positive energy balance to survive [35] and partition their daily activities to do so; therefore, time-energy allocation and activity budgets provide valuable tools for assessing animals' response to their environment [36]. Advancements in animal-borne datalogger technology now allow researchers to investigate fine-scale behavioral responses of wild animals [37][38][39]. Accelerometers, for example, record high-resolution three-dimensional movement, which can be used to calculate an index of activity, such as vectorial dynamic body acceleration (VeDBA, measured in gravitational g) [40,41]. Furthermore, an animal's time-activity budget can be estimated using a hidden Markov model (HMM) to decode behavioral state from acceleration [42,43].
Our objective was to study juvenile caracaras' activity levels in the presence and absence of seasonally migratory prey on Saunders Island, Falklands, which supports seven species of breeding seabirds during summer (hereafter the "seabird colonies site") and has a farm settlement that provides human subsidies (e.g., animal feed, offal, and compost) year-round. To do this, we deployed animal-borne accelerometer and GPS dataloggers on 27 juvenile caracaras during summer and winter and investigated seasonal variation in body mass of the tagged subset and an additional 65 caracaras. If juvenile caracaras were food stressed and unable to forage optimally during winter, we predicted that caracaras would (1) have a higher daily activity index (i.e., 24-h average VeDBA); (2) engage in higher intensity daytime activity (i.e., average VeDBA during daylight); (3) range farther; and (4) have lower mass than during summer.

Daily activity budgets
Hidden Markov models (HMMs) are dependent mixture models in which several distributions are used to account for an animal's behavioral modes [42]. We found that a 4-state HMM fit the data best ( Fig. 2a) while comprising clearly interpretable states: (1) rest (e.g., sleep or a similar motionless posture); (2) rest with noise (e.g., preening or daytime rest); (3) low activity (e.g., walking or raking); and (4) high activity (e.g., flight, running, or antagonistic behavior) (Fig. 2b). When coupled with the GPS data, resting states clustered at known roost sites and the highest activity seemingly occurred in flight (e.g., coming and going from the roost site).

Discussion
Juvenile striated caracaras on Saunders Island remain as active in winter as they are in summer, though they compress their activity into winter's shorter daylight hours. These results suggest that during winter, juvenile caracaras on Saunders do not employ the energy-saving strategies that some species exhibit when faced with seasonal food scarcity [44][45][46][47]. Moreover, we found that male caracaras are more active than females in both seasons, a possible behavioral expression of the reverse sexual size dimorphism common in the Falconidae family [48,49].
In the present study, not enough males wore GPS transmitters in winter to draw a statistically significant conclusion about differences in space use by sex. However, the one male sampled during winter traveled from the farm settlement to the seabird colonies, increasing his 95% KDE twofold over the ranges of females who remained at the settlement (male 95% KDE: 3.85 km 2 , female 95% KDE: 1.98 ± 0.54 km 2 [n = 4]). On average, male caracaras in the Falklands weigh 200 g less than females, and as cost of transport scales positively with mass, we hypothesize their smaller size could facilitate their higher activity levels [50][51][52]. Since these are juvenile birds, they are not involved in breeding activity, so this cannot be due to nest provisioning. Notably, the seasonal daily ranges we measured were comparable to the multiday ranges reported in Harrington et al. (2018), which tracked 6 birds with GPS dataloggers over 5-20 days during December 2010 and 5 birds over 3-9 days during July-August 2012 and 2017.
That the caracaras' daily activity and overall body mass did not vary between the seasons suggests that wintering birds on Saunders Island are meeting their resource requirements despite the absence of seasonally migratory seabird prey. We hypothesize that resources available at the nearby farm settlement may be mitigating the effect of seasonal food limitations.
However, juvenile caracaras appear to respond atypically to human subsidies. While human subsidies are  (Fig. 2a), representative behaviors per state (Fig. 2b), and a sample 10-min time series of the state-labeled VeDBA values from a summer deployment (Fig. 2c) available year-round at the Saunders settlement, juvenile caracaras only gather and roost at the farm settlement when resources associated with the breeding seabird colonies dwindle. For other species, predictable food resources have encouraged residency (e.g., [16,53,54]). Human subsidies have also encouraged corvidsscavenging predators that fill a similar ecological niche as the caracaras-and canids to reduce and shift their home ranges closer to subsidized resources [55,56]. Caracaras on Saunders Island, on the other hand, have not become year-round residents at the settlement or reduced their home range to remain near the subsidies. Instead, they have increased winter ranges and leave the settlement as soon as their natural food supply increases (i.e., when the seabirds return to breed). This suggests that the human subsidies are less attractive to the birds than naturally available resources. Additionally, not all caracaras in the Falklands overwinter at farm settlements. For example, a 2012 austral winter survey counted 274 caracaras of all age classes on Steeple Jason Island, a nature preserve 80 km northwest of Saunders, where the caracaras were observed hunting terrestrial invertebrates, foraging in the intertidal and kelp wrack, and scavenging at gentoo colonies [57].
One limitation to our study was the use of two versions of our datalogger. While tag version was significant in modeling daily average VeDBA and daytime average VeDBA, when comparing R-squared values, we determined tag version did not explain an ecologically relevant proportion of the data. Biological factors explained 32% of the daily average VeDBA and 72% of the daytime average VeDBA, while tag version only explained an additional 13% and 6%, respectively. Additionally, when comparing between seasons measured by the same tag version, we achieved similar results. It is possible the effect was due to year and not tag version, though we were unable to include year and  (Fig. 3a) and corresponding percentage (Fig. 3b) of daytime (i.e., solar dawn to solar dusk) spent in each activity state for summer (solid) and winter (dashed). Open circles (males) and triangles (females) reveal sex differences within corresponding boxplots.
Year-round, males spent a greater proportion of daylight hours in states 3 and 4 than females. During winter, both sexes increased their proportion of time spent in states 3 and 4. Caracaras appear to spend the extra daylight of longer summer in resting states (Fig. 3a). Boxplot whiskers extend to the largest and smallest values (at most 1.5 * IQR [inter-quartile range]) tag version in the same model because we did not use both tag versions in all years.

Conclusions
We used animal-borne dataloggers to better understand behavioral response to seasonal resource fluctuation. Combining GPS data with high-resolution accelerometry is especially useful for understanding activity states in an ecological context. In addition, HMMs provide a data-driven framework for evaluating unobserved activity states. While we did not find behavioral evidence of food stress, further research should include studies on seasonal energetics to improve our understanding of baseline body condition, and comparative studies of adult caracaras and islands without human subsidies. As many island-restricted species will likely face increased variation in resource availability in response to changing conditions including rising sea levels and human population expansion, we suggest that conservation managers consider seasonal activity budgets when developing species management plans.

Study area
This study was conducted on Saunders Island (51.37°S 60.09°W), a 127-km 2 island in the northwest Falkland Islands archipelago that supports an estimated 50 to 100 caracaras year-round [31, KH unpubl. data]. During austral summer, the island supports over 25,000 breeding seabirds (see [17] for description) at the northwest end of the island. Following the breeding season, all seabirds at the colonies site-except the gentoos who roost near their breeding colonies throughout the year-migrate offshore. When this happens, caracaras vacate their summer roost and move 16 km SE to the 2.5-km 2 farm settlement, where they overwinter at a roost adjacent to the settlement. The farm is owned and operated by a family of four and 2 co-workers who maintain approximately 6500 sheep, 200 dairy and beef cattle, approximately 400 cagefree chickens and ducks, and 22 farm dogs. Five sheep are processed each week for family use, and offal is left for the scavenging birds. The family also keeps two pigs in a small, open-air pigpen, which are fed two upland geese (Chloephaga picta) daily; many caracaras attend the feeding during winter and feed on scraps as available. The pig Fig. 4 Combined 24-h kernel home range (i.e., 95% KDE, white polygon) and core area of use (i.e., 50% KDE, black polygon) during winter (n = 5; Fig. 4a) and summer (n = 7; Fig. 4b) for all individuals with GPS units. The 100% MCP (gray outer polygon) is included for comparison. The gray rectangle on Fig. 4a represents the spatial extent of Fig. 4b inset, which is the seabird colonies site. The farm settlement is indicated on Fig. 4a. Colored points represent different individuals tracked for 24-h at 2-min sampling rates feeding, weekly sheep processing, and dogs' bones appear to be the primary draw for the caracaras in winter [17]. The farmers have not observed caracaras predating the free-ranging farm animals at the settlement, although they have seen caracaras approach and pick at the nose and tails of recently birthed lambs [KH pers. comm.]. No seabirds roost or breed at the settlement.

Animal handling and instrumentation
From 2017 to 2019, we trapped 92 juvenile caracaras (Phalcoboenus australis, Gmelin 1788), 62 at the seabird colonies site during austral summer (December-February) and 30 at the farm site during austral winter (June-August). We targeted juvenile caracaras (i.e., second year birds, aged by plumage per [28]) to remove the effect of breeding on seasonal movement, activity levels, and mass. It was not possible to trap at both sites during the same season, because the birds seasonally migrate between the two sites [17]. Caracaras were trapped using a mutton-baited monofilament snare carpet that pulls tight around their metatarsus as the birds walk through the trap [58]. The trap is monitored continuously, and caracaras are removed immediately after snaring. All captured caracaras were weighed to the nearest 50 g using a hand-held scale and banded with colored alphanumeric plastic leg bands [17]. Approximately 50 μL of blood was collected from the brachial vein using a sterile 21 gauge hypodermic needle [59] and transferred immediately onto an 8-µm filter paper [60] for molecular sexing (modified protocol from [61]). Caracaras have distensible crops, thus we visually assessed crop size and categorized as not distended, partially distended, or distended, which allowed us to control for the added mass of food storage by restricting our mass analysis to birds that did not have distended crops. We observed no adverse effects of trapping and handling, with many individuals continuing to feed immediately upon release.
To 27 of those caracaras-13 at the seabird colonies site and 14 at the farm settlement-we attached custom-built Arduino-based archival dataloggers (TWLogger; [62]). We attached the loggers to the bird's two central feathers using Tesa 4651 adhesive tape (Norderstedt, Germany) [63] (Fig. 5). The dataloggers were < 2% of birds' body mass (1649.44 ± 177.08 g, mean ± SD [n = 27]); and we assumed that any unexpected behavioral effect of the dataloggers would be pervasive across our sample and thus not affect our comparative index. For all deployments, loggers recorded tri-axial acceleration (50-Hz sampling rate) data, which were used to calculate VeDBA and decode behavioral states. For 2018 and 2019 deployments, loggers also recorded GPS positions (2-min sampling rate), which were used to help validate behavioral state classifications and estimate utilization distributions and core areas of use. Prior to deployment, we spherically rotated all loggers to pass through all combinations of pitch and roll to record reference values to calibrate accelerometer data to account for inherent sensor variability (i.e., due to manufacturing) [64]. After deployment, all instrumented individuals resumed normal activity.
The present study occurred concurrently with TWLogger development, which was an iterative process. Tags deployed in 2017 (dimensions: L52 × W23 × H10 mm) used an inertial movement unit (IMU) that included a gyroscope in addition to accelerometer sensors and did not include a GPS module. Tags deployed in 2018 and 2019 were slightly larger (dimensions: L52 × W23 × H17 mm), included a GPS, and used a smaller IMU that omitted the gyroscope; this version consumed less power Fig. 5 The TWLogger inertial movement datalogger with GPS prior to encapsulating in heat shrink for deployment (Fig. 5a), and a banded caracara wearing a tail-mounted TWLogger (Fig. 5b)  To retrieve loggers, we recaptured birds using a mutton-baited single snare. Deployment durations ranged from 1.0 to 12.2 days (4.3 ± 3.6 days, mean ± SD), though logger battery durations ranged from 26 to 40 h. All but two loggers were recovered. Despite intensive efforts, we were unable to relocate one bird tagged during austral summer 2017 and another tagged during austral winter 2018. In both cases, the loggers would have been shed during their natural molt (Nov-Jan), if not before, preventing long-term consequences for the birds. Only one of 25 recovered loggers failed to record data.

Data processing
We obtained 24 complete datasets of tri-axial acceleration (winter n = 12, summer n = 12). Of these, 12 loggers also recorded GPS data (winter n = 5, summer n = 7). Each dataset represented one individual; thus, all observations were independent.
Prior to analysis, we processed acceleration data using purpose-written scripts in R (R v. 3.5.2). First, we used the spherical calibration function in the tagtools package [65] in R to calibrate the data using the reference values recorded prior to deployment. This allowed us to convert raw sensor values to acceleration in m s −2 while accounting for inherent sensor variability due to manufacturing effects. Next, we corrected any sampling rate variability, which accounted for less than 0.001% of the data. Because we were assessing hourly and daily scales and lower resolution acceleration data maintains high accuracy [66], we down sampled the data to 10-Hz to improve data processing speeds [67]. From the 2018 and 2019 deployments that included GPS units, we obtained 12 complete sets of tri-axial acceleration and GPS data. All GPS data were processed using an SDA filter (i.e., speed-distance-angle filter) [68] in R to remove aberrant positions, as environmental conditions may have affected accuracy of satellite fixes. GPS data were then rejoined to the dataset.
All dataloggers recorded data for greater than 24 h; although to ensure data were comparable on a daily scale and to minimize the effects of tag deployment, we clipped from the beginning of each dataset to make all files 24 h in length.
We classified each data point as day or night to assess diel patterning and standardize statistics such as daytime hourly VeDBA rates. We calculated solar periods (i.e., dawn, day, dusk, and night) using the maptools package in R [69] that performs sun ephemerides calculations using algorithms developed by the National Oceanic & Atmospheric Administration (NOAA). Solar dawn and dusk are defined as the time when the solar positioning is between 0 and 12 degrees below the horizon (nautical twilight). Night is defined as the time when the solar positioning is below 12 degrees on the horizon. For this analysis, "daytime" comprises dawn, day, and dusk periods. We calculated daytime length as the period between solar dawn and dusk and created 1-h solar time bins from zero to 23 beginning at solar midnight to facilitate seasonal comparisons.

VeDBA calculations
To calculate VeDBA (as described by [41]), we first calculated dynamic acceleration from each axis by subtracting static acceleration (i.e., smoothed acceleration over a 2-s window) from raw acceleration values. We then found the square root of the sum of squared dynamic acceleration of each axis. Using the 10-Hz VeDBA values, we calculated 24-h average VeDBA (i.e., daily activity index) and average daytime VeDBA per individual (i.e., VeDBA index subset to daylight hours and divided by day length). We reported VeDBA metrics in gravitational g.

State classifications
We used hidden Markov models (HMMs) to infer behavioral states from the 10-Hz acceleration data as they allow identification of relevant changes in behavioral intensity, i.e., activity versus rest, taking into account the serial dependence between observations [42,70]. HMMs require an input metric that can be used to distinguish between movements that represent likely states. We selected VeDBA as our input metric [40,42,71]. Since VeDBA values are positive and continuous, we used gamma state-dependent distributions within the HMM.
For animal behavior data, formal model selection techniques tend to select more HMM states than seem biologically reasonable (i.e., additional states compensate for the inability to capture all the structure in such complex data); therefore, we first considered a simple (2-state) model and then built up complexity until a reasonable goodness-of-fit was achieved, but such that the model states were still biologically interpretable [72].
We used ecological knowledge (i.e., 675 h of observations over 3.5 years) and GPS positions to validate the model's state classifications (e.g., known overnight roost sites where caracaras should be in a resting state). For daylight hours, we calculated absolute state durations (hours) and proportion of time spent in each state.

GPS calculations
We used the GPS data to calculate KDEs [73,74] of the 95% home range and 50% core areas of use of all recorded locations for each individual bird (winter n = 5, summer n = 7) using the kernelUD and kernel.area functions in the adehabitatHR package in R [75]. To allow comparison between seasonal home ranges of different sizes, we used a fixed value smoothing parameter h (i.e., bandwidth) for each study area. To determine our smoothing parameters, we first subset the data by season and used the href method in the kernelUD function to calculate h per individual within a season. We then used the median h per season as our smoothing parameter [76]. For comparison, we also calculated 100% MCPs [77,78] per bird using the mcp and mcp.area function from the adehabi-tatHR package.

Statistical analyses
We ran all statistical analyses in R (R v. 3.5.2). We tested for seasonal differences in 24-h average VeDBA, average daytime VeDBA, proportion of daylight hours spent in each behavioral state, and body mass using analyses of variance with season, sex, and device model as covariates. We selected the most parsimonious model to our data by comparing our initial model to subsequent models after removing non-significant fixed terms using likelihood-ratio tests. We report best model with tag version, but compare R-squared with and without tag version to understand importance. To compare KDEs and MCPs between seasons, we used two-sample t-tests and a Wilcoxon rank sum test, as our dataset contained only one year per treatment.