Compound-V formations in shorebird flocks
Figures

Shorebird flock recording.
(a) Multi-camera videography was used to reconstruct 3D trajectories of shorebirds flying near high-tide roosts in Humboldt Bay, California. (b) Overhead and (c) profile views of an example flock. Symbol sizes reflect actual scales for birds with outstretched wings. Flock position data are available in Figure 1—source data 1.
-
Figure 1—source data 1
Three-dimensional trajectory data for godwit flock 0420–1.
- https://doi.org/10.7554/eLife.45071.004

Within-flock positioning.
(a, b) Histograms of nearest-neighbor alignment for birds flying within ±1 wingspan of elevation (godwit flock 0420–1) show a distinctive peak at a trailing distance and lateral distance of approximately one wingspan; focal birds are shown in light gray and nearest neighbors in black. Inset bird silhouettes show profile views of the birds’ relative flight elevations. (c, d) Histograms of nearest-neighbor alignment for birds flying outside ±1 wingspan of elevation for the same flock show a largely random distribution with a modal location of nearly straight above or below the focal bird. Data used for generating this figure are available in Figure 2—source data 1.
-
Figure 2—source data 1
Nearest-neighbor positioning data for flock 0420–1.
- https://doi.org/10.7554/eLife.45071.008

Modal positioning among flocks and species.
(a) Summary of modal neighbor position for nearest neighbors within ± 1 wingspan in single-species flocks of all four species, depicted in absolute metric distance and (b) the same data plotted in distances relative to the wingspan of each species. Open symbols indicate modal neighbor positions for individual flocks. Closed symbols and silhouettes show the average position for each species. Data used for generating this figure are available in Figure 3—source data 1.
-
Figure 3—source data 1
Modal neighbor positioning data for all shorebird flocks and for a flock of chimney swifts.
- https://doi.org/10.7554/eLife.45071.012

Distribution of nearest neighbors in the horizontal slice (±1 wingspan elevation) from a chimney swift roosting flock (Evangelista et al., 2017).
The silhouettes show the position of the focal bird at the origin and the lateral and horizontal position modes. The dashed line box shows the ‘aerodynamic neighbor’ region for comparison with shorebird flocks Figure 3—figure supplement 2. Note that the lateral and horizontal modal positions are calculated separately, and the result is not necessarily congruent with the most populated 2D grid cell. Data used for generating this figure are available in Figure 3—source data 1.

Distribution of nearest neighbors in the horizontal slice (±1 wingspan elevation) from all shorebird flock data described here, regardless of species.
The dashed line box shows the ‘aerodynamic neighbor’ region for comparison with shorebird flocks, 29.3% of all nearest neighbors fall within this region. Data used for generating this figure are available in Figure 3—source data 1.

Positioning in mixed-species flocks.
Data from mixed species flocks show that birds adjust their lateral spacing depending on the species (and size) of their nearest leading neighbor. (a) Godwits following conspecifics adopt a larger lateral distance than (b) godwits following the smaller dowitchers. (c) Dowitchers following conspecifics use a shorter lateral distance than (d) dowitchers following the larger godwits. These results support the hypothesis that shorebirds adopt a lateral spacing rule that is dependent on the size of their leading neighbor. Dashed lines are provided to facilitate comparison of modal lateral positions between (a) and (b) and between (c) and (d). Data used to generate this figure are available in Figure 4—source data 1.
-
Figure 4—source data 1
Neighbor position data for mixed-species flocks.
- https://doi.org/10.7554/eLife.45071.014

Godwit simple-V formation.
Incidental to our cluster flock recordings, we also recorded several instances of godwits flying in a simple-V, echelon or line formation, and the largest of these examples is shown here. (a) Overhead view of the flock; average flight direction is along the positive Y axis; blue circles show bird positions and black lines are 2D velocity vectors. All birds are within a ± 1 wingspan horizontal slice. (b) The relative location of nearest neighbors; the modal location (red circle) was at a displacement of 0.8 wingspans lateral and 0.5 wingspans trailing distance. Trailing position was more varied than lateral position. Wind speed was low (<2 m s−1) according to weather station data and the wind speed estimated from the ground speed and flight direction of the birds. The data used to generate this figure are available in Figure 3—source data 1.
-
Figure 5—source data 1
Godwit simple-V-formation position data.
- https://doi.org/10.7554/eLife.45071.016

Extended flock structure.
(a) Polar plot showing mean neighbor angle for right-aligned and left-aligned flocks over a range of distances. Shaded regions show 95% confidence intervals. (b) Overhead and profile views of an example right-aligned flock (avocet flock 1220–2). Note the many echelon formations aligned from back left to front right and the overall shape of the flock. The inset shows scale in wingspans. The data used to generate this figure are available in Figure 6—source data 1.
-
Figure 6—source data 1
Extended flock structure data.
- https://doi.org/10.7554/eLife.45071.018

Effect of positioning for aerodynamic interaction.
Here we show the effect of neighbor position on flight speed. (a) Flight speed residuals after accounting for species, flapping frequency, distance from flock edge, nearest neighbor distance in terms of wingspans and overall position along the length of the flock. (b) Flight speed residuals after accounting for just species and flapping frequency. White spaces in the heat map are bins with fewer than 20 samples, out of 2848 possible in (a) and 3306 possible in (b). Both analyses reveal a broadly similar pattern, where the positive effect of neighbor position on flight speed is strongest at a one wingspan lateral displacement and a trailing distance of 0 to 0.5 wingspans. This pattern cannot be generated by trailing birds passing leaders because the roles reverse after passing occurs, leaving no net speed difference. The data used to generate this figure are available in Figure 7—source data 1.
-
Figure 7—source data 1
Neighbor position and flight speed data.
- https://doi.org/10.7554/eLife.45071.024

Species identification in mixed flocks.
(a) Histogram of scaled pixel area of birds within a mixed-species flock. The two peaks are modeled as normal distributions. The area value where the two distributions intersect (indicated by the arrow) is used as the threshold for species identification. (b) Example section of a mixed-species flock with species identifications labeled by color.

Determining flock edge and maximum radius.
An overhead view of an example flock of avocets (flock 1220–2). Because flocks were always spread out in the horizontal direction, a compact hull is fitted to the XY-coordinates to create a boundary. The minimum horizontal distance of each bird to the hull is the bird’s edge distance. The median edge distance is taken as the flock’s maximum radius for computing alignment metrics (Figure 6). Here, birds within the maximum edge distance (6.5 wingspans or 4.55 m) are labeled edge, and birds beyond the maximum edge distance are labeled core.
Tables
Flock parameters.
https://doi.org/10.7554/eLife.45071.005Flock | Species | N birds | N frames | †, ‡Nearest neighbor distance (m). | §Nearest neighbor power | †Ground speed (m·s−1) | †Airspeed (m·s−1) | Wind speed (m·s−1) | ¶Wind direction (deg.) | †Z-speed (m·s−1) | †Turnrate (° s−1) |
---|---|---|---|---|---|---|---|---|---|---|---|
0417–2 | Godwit | 286 | 147 | 1.30 (0.79–2.15) 1.67 | 0.361 | 5.23 (3.59–8.23) | 9.36 (8.04–10.54) | 5.73 | 4.3 | 1.11 (0.16–1.91) | 40.5 (18.4–80.9) |
0417–2 | Dowitcher | 309 | 147 | 1.16 (0.72–1.89) 2.23 | 0.385 | 5.15 (3.39–8.42) | 9.25 (7.94–10.47) | 5.73 | 4.63 | 1.26 (0.30–2.00) | 43.0 (19.4–81.9) |
0417–3 | Godwit | 474 | 278 | 1.81 (1.01–2.94) 2.32 | 0.428 | 3.32 (1.71–6.55) | 7.87 (6.42–9.39) | 5.73 | 5.3 | −0.18 (−1.27 – 0.76) | 26.8 (8.0–80.6) |
0417–4 | Godwit | 803 | 397 | 1.71 (1.02–2.58) 2.19 | 0.391 | 6.30 (3.04–9.44) | 8.82 (7.68–9.93) | 5.73 | 47.2 | 0.22 (−0.73 – 1.03) | 14.3 (4.66–40.8) |
0417–4 | Dowitcher | 74 | 397 | 1.57 (0.94–2.37) 3.01 | 0.382 | 6.51 (3.00–9.56) | 8.65 (7.49–9.65) | 5.73 | 48.1 | 0.39 (−0.53 – 1.21) | 19.8 (6.6–49.2) |
0420–1 | Godwit | 639 | 177 | 1.19 (0.59–1.94) 1.52 | 0.408 | 7.05 (5.12–9.91) | 10.54 (7.91–12.98) | 4.06 | 9.2 | −0.26 (−1.83 – 1.03) | 22.6 (7.39–75.9) |
0420–2 | Godwit | 309 | 147 | 1.54 (0.83–2.38) 1.97 | 0.43 | 8.95 (7.09–11.59) | 10.50 (8.76–12.11) | 4.06 | 59.7 | −0.06 (−0.95 – 0.68) | 11.9 (3.77–37.7) |
0427–2 | Dowitcher | 354 | 397 | 1.07 (0.61–2.00) 2.06 | 0.424 | 6.22 (3.06–9.76) | 9.03 (5.53–12.77) | 3.50 | 23.5 | 0.56 (−0.89 – 1.79) | 20.2 (6.68–60.1) |
0427–3 | Dowitcher | 391 | 217 | 1.15 (0.61–2.06) 2.21 | 0.463 | 10.9 (8.64–13.3) | 10.98 (7.99–13.98) | 3.50 | 73.3 | 0.32 (−2.10 – 2.15 | 24.8 (8.30–70.9) |
0427–5 | Dowitcher | 511 | 170 | 1.23 (0.70–1.98) 2.37 | 0.421 | 5.31 (4.23–6.49) | 7.14 (4.85–9.18) | 4.60 | 26.9 | 0.73 (−0.24 – 1.52) | 20.8 (7.62–57.4) |
1230–1 | Dunlin | 351 | 198 | 1.08 (0.58–1.78) 3.18 | 0.465 | 6.98 (6.03–7.85) | 7.65 (6.67–8.34) | 1.20 | 44.1 | −0.23 (−0.52 – 0.09) | 27.3 (11.9–46.6) |
1230–2 | Dunlin | 592 | 75 | 0.80 (0.47–1.23) 2.35 | 0.387 | 6.61 (6.00–7.30) | 7.60 (5.84–8.48) | 1.10 | 22.6 | −0.12 (−0.56 – 0.32) | 11.4 (3.7–26.8) |
1230–3 | Dunlin | 477 | 125 | 0.86 (0.49 1.37 2.53 | 0.392 | 6.71 (5.85–7.73) | 7.54 (5.52–8.36) | 1.08 | 35.5 | 0.11 (−0.49 – 0.84) | 19.6 (5.27–42.9) |
1230–4 | Dunlin | 189 | 73 | 0.89 (0.50–1.50) 2.62 | 0.502 | 8.28 (7.44–9.70) | 7.47 (5.27–8.53) | 1.08 | 36.5 | −0.03 (−0.57 – 0.45 | 24.6 (7.5–52.9) |
0101–1 | Dunlin | 1039 | 228 | 1.03 (0.59–1.64) 3.03 | 0.41 | 8.39 (7.41–9.98) | 7.46 (5.74–8.64) | 1.63 | 118.4 | −0.23 (−0.73 – 0.27) | 17.7 (4.9–42.1) |
0101–3 | Dunlin | 961 | 188 | 0.92 (0.52–1.50) 2.71 | 0.416 | 8.61 (7.70–9.42) | 7.74 (6.25–8.76) | 1.63 | 117.6 | −0.02 (–0.33 – 0.36) | 13.0 (3.6–32.2) |
0101–4 | Dunlin | 269 | 340 | 1.00 (0.35–2.12) 2.94 | 0.45 | 6.56 (4.96–8.12) | 7.63 (5.12–8.72) | 1.63 | 39.3 | −0.14 (−0.70 – 0.30) | 18.3 (4.4–46.9) |
1220–1 | Avocet | 323 | 90 | 1.09 (0.70–1.69) 1.51 | 0.429 | 6.02 (4.55–7.57) | 8.18 (6.21–9.31) | 2.39 | 1.4 | 0.33 (−0.98 – 0.90) | 25.2 (10.4–48.8) |
1220–2 | Avocet | 321 | 245 | 1.19 (0.72–1.90) 1.65 | 0.432 | 6.96 (5.26–9.22) | 8.00 (6.88–8.96) | 2.39 | 8.2 | 0.10 (−1.49 – 0.74) | 30.2 (12.7–54.5) |
1220–3 | Avocet | 281 | 280 | 1.30 (0.78–2.10) 1.81 | 0.472 | 7.50 (5.32–8.93) | 7.93 (6.28–8.88) | 2.39 | 22.3 | 0.33 (−0.62 – 0.83 | 23.9 (6.18–49.1) |
-
†Values are medians and (in brackets) 10th-90th percentiles of values extracted at one-wingbeat intervals from all individuals of each flock.
‡Values in italics are in wingspan units instead of metric units.
-
§Exponent of power law fit to distance of 10 nearest neighbors.
¶Wind direction is relative to the overall flight direction where 0° is a pure headwind and 180° a pure tailwind.
-
Note that data are presented separately in consecutive rows for each species in mixed-species flocks (0417–2 and 0417–4). Data used for generating this table are available in Table 1—source data 1.
-
Table 1—source data 1
Flock parameter data.
- https://doi.org/10.7554/eLife.45071.006
Flock orientation.
https://doi.org/10.7554/eLife.45071.019Variable | Test | N | R2 | t/F | P |
---|---|---|---|---|---|
Wind direction | Circular correlation | 18 | 0.29 | 2.29 | 0.02 |
Turn direction | Linear regression | 18 | 0.00 | 0.04 | 0.83 |
Camera direction | Circular correlation | 18 | 0.02 | 0.61 | 0.54 |
-
Tests of the relationship between flock left-right orientation (Figure 6) and environmental factors. The data usedto generate this table are available in Table 2—source data 1.
-
Table 2—source data 1
Flock orientation data.
- https://doi.org/10.7554/eLife.45071.020
Flock biomechanics.
https://doi.org/10.7554/eLife.45071.021Wingbeat frequency predictors | Estimate | S.E. | T | d.f. | P |
---|---|---|---|---|---|
Intercept (dowitchers) | 8.82 | 0.132 | 66.5 | 2817 | <0.00001 |
Godwit | −2.19 | 0.035 | −63.5 | 2817 | <0.00001 |
Avocet | −1.91 | 0.040 | −47.5 | 2817 | <0.00001 |
Airspeed (m s−1) | −0.05 | 0.013 | −4.3 | 2817 | 0.00002 |
Flock position | 0.13 | 0.032 | 4.0 | 2817 | 0.00006 |
Nearest neighbor distance (wingspans) | −0.03 | 0.010 | −3.2 | 2817 | 0.00153 |
Nearest neighbor species | −0.45 | 0.038 | 11.9 | 2817 | <0.00001 |
Z-speed (m s−1) | 0.29 | 0.019 | 15.1 | 2817 | <0.00001 |
Airspeed predictors | Estimate | S.E. | T | d.f. | P |
Intercept (dowitchers) | 10.69 | 0.225 | 47.5 | 2832 | <0.00001 |
Godwit | −0.25 | 0.067 | −3.7 | 2832 | 0.00022 |
Avocet | −1.51 | 0.067 | −22.6 | 2832 | <0.00001 |
n.n. distance (wingspans) | −0.08 | 0.015 | −5.2 | 2832 | <0.00001 |
Edge distance (wingspans) | −0.03 | 0.003 | −10.3 | 2832 | <0.00001 |
Wingbeat frequency (H z) | −0.12 | 0.025 | −5.0 | 2832 | <0.00001 |
Flock position | −0.24 | 0.045 | −5.4 | 2832 | <0.00001 |
Aerodynamic neighbor | 0.23 | 0.040 | 5.7 | 2832 | <0.00001 |
-
Results from linear mixed-effects models relating wingbeat frequency and airspeed to other measured variables. Nearest neighbor only defined when this bird is leading the focal bird. Godwit and avocet are dummy variables coding species differences relative to dowitchers. Nearest neighbor species is coded −1 for a smaller neighbor, 0 for same the species, and 1 for a larger neighbor. Flock position is continuously scaled from 0.0 (front) to 1.0 (back). Aerodynamic neighbor was coded 1 for birds flying within 0.7–1.5 wingspans lateral distance and within two wingspans distance from their nearest leading neighbor, 0 otherwise. Models were selected using Bayesian information criteria. The data used to generate this table are available in Table 3—source data 1. d.f., degrees of freedom.
-
Table 3—source data 1
Flock biomechanical data.
- https://doi.org/10.7554/eLife.45071.022
Additional files
-
Transparent reporting form
- https://doi.org/10.7554/eLife.45071.027