Data from: Unexpected migration patterns in a high-latitude breeding songbird: Evidence from multi-sensor geolocators and isotopes
Data files
Jan 08, 2026 version files 76.51 MB
-
CG179.zip
13.34 MB
-
CG181.zip
12.89 MB
-
CG193.zip
12.65 MB
-
CG194.zip
12.84 MB
-
CG206.zip
12.84 MB
-
CG207.zip
11.92 MB
-
README.md
9.43 KB
-
Supporting_Dataset_1_Geolocator_and_Control_Bird_Sampling.csv
5.83 KB
-
Supporting_Dataset_2_Isotope_Sampling.csv
11.74 KB
Abstract
Seasonal migration allows animals to use habitat where conditions are unfavorable for part of the year but may constrain breeding ranges due to the costs of longer migrations as ranges expand poleward. In species with large ranges, high-latitude breeding populations may employ different migration strategies allowing them to persist far from other core nonbreeding areas. The myrtle warbler (Setophaga coronata coronata) has two disjunct nonbreeding ranges in North and Central America—one along the Gulf Coast and the other on the Pacific. Previous work indirectly linked birds breeding in Alaska with the Pacific nonbreeding area, suggesting that high latitude populations evolved a shorter migration route. We directly tested this hypothesis using geolocators measuring both light and atmospheric pressure to track Alaskan myrtle warbler migration in fine detail and inferred nonbreeding areas using hydrogen isotopes for a larger sample of birds breeding in Alaska, British Columbia, and Alberta. We found, contrary to expectations, that all geolocator-tracked birds—and 95% of birds with stable isotope data—migrated to the southeastern United States, a much longer migration than expected for a species commonly considered a “short-distance” migrant. We additionally demonstrate the advantages of pressure geolocation for characterizing migratory behavior at a fine scale.
https://doi.org/10.5061/dryad.cfxpnvxhp
Description of the data and file structure
Full-year multisensor geolocator data (light, atmospheric pressure, temperature) collected for six myrtle warblers breeding in Anchorage, AK in June 2022-23. Stable hydrogen isotope data generated from greater covert feathers collected from 167 myrtle warblers breeding in Anchorage, AK, northern British Columbia, and Alberta.
Files and variables
Compressed archives (CG179.zip, CG181.zip, CG193.zip, CG194.zip, CG206.zip, CG207.zip) each contain data from one geolocator
Archives are named for the geolocator ID (recorded in Geolocator_and_Control_Bird_Sampling.csv). All dates and times are GMT/UTC. Within each archive are the following directories and files:
Directory: GeoID_Raw Data – Contains raw data files produced by a geolocator
GeoID_Date_driftadj.lux
- Description: Light level data adjusted for clock drift. Data has been trimmed to exclude the period after the geolocator was activated but before it was deployed on a bird, and after it was removed from the bird but before it was deactivated.
- Variables:
- DD/MM/YYYY: Date measurement was recorded
- HH:MM:SS: Time measurement was recorded
- light(lux): Measured light level in lux
GeoID_Date.deg
- Description: Atmospheric pressure and temperature data collected by a geolocator. Data has been trimmed to exclude the period after the geolocator was activated but before it was deployed on a bird, and after it was removed from the bird but before it was deactivated.
- Variables:
- DD/MM/YYYY: Date measurement was recorded
- HH:MM:SS: Time measurement was recorded
- T('C): Temperature in degrees Celsius
- P(Pa): Atmospheric pressure in Pascals
Directory: GeoID_GeoLight – Contains output from analysis of light level data using the R packages TwGeos and GeoLight
GeoID_geolight_twilights.csv
- Description: Twilight times estimated from geolocator light-level data using the R package TwGeos and formatted for GeoLight. Outliers were removed using the twilightEdit function in TwGeos and removal criteria of window = 4, outlier.mins = 25, and stationary.mins = 20.
- Variables:
- tFirst: Date and time of first twilight event
- tSecond: Date and time of second twilight event
- type: 1 = tFirst is sunrise, 2 = tFirst is sunset
- sun_angle: Sun angle used for calibration
- zero_angle: Zero angle used for calibration
GeoID_geolight_locations_fullyear.csv
- Description: Locations estimated in GeoLight for a geolocator over the full year from twilight table (GeoID_geolight_twilights.csv)
- Variables:
- tFirst: Date and time of first twilight event
- lon: Longitude estimate in decimal degrees for a set of twilights
- lat: Latitude estimate in decimal degrees for a set of twilights
GeoID_geolight_locations_winter_spring.csv
- Description: Locations estimated in GeoLight for a geolocator over the wintering period and spring migration only. These points are plotted in Fig. S10, where fall migration coordinate estimates were omitted due to high latitudinal error around the fall equinox.
- Variables:
- tFirst: Date and time of first twilight event
- lon: Longitude estimate in decimal degrees for a set of twilights
- lat: Latitude estimate in decimal degrees for a set of twilights
Directory: GeoID_GeoPressureR – Contains output from analysis of atmospheric pressure and light level data using the R package GeoPressureR
GeoID-tag-labeled.csv
- Description: Pressure data labeled in TRAINSET with flights, elevation changes, and outliers to discard.
- Variables:
- series: Data type (here, pressure)
- timestamp: Date and time pressure was recorded
- value: Atmospheric pressure measurement (hectopascals)
- label: Blank if no label, discard=outlier to be discarded, flight=migratory flight, elev_#: label for bird changing elevation within a stationary period
GeoID-twilight-labeled.csv
- Description: Twilight data labeled in TRAINSET with stationary periods and outliers to discard.
- Variables:
- series: Twilight type (Set = sunset, Rise = sunrise)
- timestamp: Date and time pressure was recorded
- value: Twilight time converted to numerical value for plotting
- label: Numbers indicating stationary period #, "discard" for outliers to remove
Directory: GeoID_GeoPressureR_Output_PressureOnly – Contains HTML output from GeoPressureR analysis of pressure data alone
GeoID_altitude_pressonly.html
- Plot of altitude over time
GeoID_margprob_path_pressonly.html
- Marginal probability map for each stationary period
GeoID_mostlikely_path_pressonly.html
- Map of most likely migration path
GeoID_mostlikely_pressurepath_pressonly.html
- Plot of atmospheric pressure over time recorded by the geolocator (grey) vs. pressure recorded in the ERA5 weather reanalysis dataset for the most likely location of the stationary period (colored). Black dots are discarded data points and orange triangles are points considered outliers but retained in the dataset.
GeoID_pressure_likelihood_pressonly.html
- Likelihood map for each stationary period based on pressure data
Directory: GeoID_GeoPressureR_Output_PressureAndLight – Contains HTML output from GeoPressureR analysis of pressure and light data combined
Contains all files described above for the "Pressure Only" analysis, with the addition of:
GeoID_light_likelihood.html
- Likelihood map for each stationary period based on light data
File: Supporting_Dataset_1_Geolocator_and_Control_Bird_Sampling.csv
Description: Banding, location, and morphological data for each bird fitted with a geolocator or color banded as a control for assessing return rates.
Variables
- record_number: Unique ID number for each bird banded
- capture_date: Date bird was captured
- capture_time: Time bird was captured (local time)
- latitude: Latitude of banding location (decimal degrees)
- longitude: Longitude of banding location (decimal degrees)
- elevation_ft: Elevation of banding location (feet)
- band_number: Band number
- geo_or_control: Whether bird was fitted with a geolocator or only color banded (control)
- geolocator_number: Unique ID number for geolocator, N = no geolocator
- recapped_or_resighted: Y = bird was recaptured or resighted, N = bird was not resighted
- age: Age of bird estimated from molt. SY = second year, ASY = after second year
- sex: Sex of bird determined from presence of cloacal protuberance or brood patch. M = male, F = female
- cp: Cloacal protuberance. 0 = no cloacal protuberance present, 1 = small, 2 = medium, 3 = large
- bp: Brood patch. 0 = no brood patch, 2 = well developed brood patch
- tarsus_mm: Tarsus length in millimeters
- wing_mm: Wing length in millimeters
- tail_mm: Tail length in millimeters
- culmen_mm: Culmen length in millimeters
- kipps_distance_mm: Kipp's distance (distance between longest primary and first secondary feather) in millimeters
File: Supporting_Dataset_2_Isotope_Sampling.csv
Description: Data (banding, location, morphological, isotope) for each feather sample that underwent stable hydrogen isotope analysis.
Variables
- record_number: Unique ID number for each bird banded
- state_prov: State or province where bird was banded
- latitude: Latitude of banding location (decimal degrees)
- longitude: Longitude of banding location (decimal degrees)
- age: Age of bird estimated from molt. AHY = after hatch year, SY = second year, ASY = after second year
- sex: Sex of bird determined from presence of cloacal protuberance or brood patch. M = male, F = female
- wing_mm: Wing length in millimeters
- tail_mm: Tail length in millimeters
- feather: Type of feather analyzed. Basic = grown in pre-basic molt (i.e. on previous year's breeding ground). Alternate = grown in pre-alternate molt (i.e. on non-breeding ground).
- delta2H: Hydrogen isotope ratio measured in feather sample
Code/software
Light-level data analyzed using the R packages TwGeos and GeoLight and the script: GeoLight.R
Atmospheric pressure data analyzed independently and in combination with light-level data using the R package GeoPressureR, following A User Manual for GeoPressureR and using the following scripts:
- GeoPressureR_TagLabelling.qmd – Set up configuration file, create tag label file, label file using TRAINSET, and check labels
- GeoPressureR_TagTrajectory_PressureOnly.R – Generate trajectories for each geolocator using pressure data alone
- GeoPressureR_TwilightLabelling.qmd – Create twilight label file, label file using TRAINSET, and check labels
- GeoPressureR_TagTrajectory_PressureAndLight.R – Generate trajectories for each geolocator using pressure and light data combined
Stable hydrogen isotope data analyzed using the R package assignR and the script: Isotopes.R
Access information
Other publicly accessible locations of the data:
- Geolocator data is also available on Movebank (ID# 4549104814)
We deployed 30 multisensor geolocators (Migrate Tech BARP30Z11-DIP) on myrtle warblers (Setophaga coronata coronata) breeding in Far North Bicentennial Park, Anchorage, AK in June 2022 and recovered 6 of these geolocators in June 2023. The geolocators sampled light intensity every minute, recording the maximum intensity every 5 minutes, and recorded pressure and temperature every 20 minutes.
Light data (.lux) was adjusted for clock drift, and raw geolocator data files (.lux, .deg) were trimmed to remove data recorded after activation of the tag but before deployment on the bird, and after removal from the bird but before deactivation.
Light data was processed with the R packages TwGeos and GeoLight. Pressure data was processed both alone and in combination with light data using the R package GeoPressureR.
We inferred wintering areas for 167 myrtle warblers breeding in Anchorage, Alaska (n = 46), northern British Columbia (n = 54), and Alberta (n = 67) using stable hydrogen isotopes from greater covert feathers grown during the pre-alternate molt. We also assessed stable hydrogen isotopes in feathers grown during the pre-basic molt (i.e. on breeding grounds) for myrtle warblers breeding in Alaska (n = 12) and British Columbia (n = 12). Feather sample preparation and hydrogen pyrolysis were performed at the Cornell University Stable Isotope Laboratory. Stable isotope data was processed using the R package assignR.
- Szarmach, Stephanie J.; Beam, Johanna K.; Moore, Mads et al. (2025). Unexpected migration patterns in a high-latitude breeding songbird: evidence from multi-sensor geolocators and isotopes. Movement Ecology. https://doi.org/10.1186/s40462-025-00618-6
