Module 1: Flight Biomechanics
Eagles are among the most efficient large fliers ever evolved. Their wing-loading, aspect ratio and slotted primary feathers permit a 10:1 glide ratio, thermal-soaring costs below 3 W/kg, and stoop speeds exceeding 300 km/h. This module derives Pennycuick’s flight-power equation, extracts the minimum-power speed \( V_{mp} \) and maximum-range speed \( V_{mr} \), explains slotted wingtip aerodynamics, and models the drag-minimising stoop dive used by Aquila chrysaetos.
1. Wing Morphometrics
Three dimensionless groups capture most of the aerodynamic character of a bird’s wing: wing loading \( W/S \), aspect ratio \( AR = b^2/S \), and Reynolds number \( Re = \rho v c / \mu \). For the golden eagle:
\[ \frac{W}{S} \approx 70\,\mathrm{N/m^2} \;(\;7\text{--}10\,\mathrm{kg/m^2}),\qquad AR = \frac{b^2}{S} \approx 7\text{--}8 \]
The wing loading of large eagles sits midway between that of shorebirds (\( W/S \lesssim 20 \) N/m2) and the albatross (\( W/S \gtrsim 150 \) N/m2). Intermediate values permit thermaling in narrow thermals while still being able to use flapping flight for attack and migration. An aspect ratio near 7 is the evolutionary compromise between the high-AR gliders (albatross AR ≈ 15) and the low-AR forest hawks (Accipiter AR ≈ 5).
Slotted primary feathers — Graham 1934
At the wingtip, the outer 5–7 primary feathers separate into distinct narrow wings during slow flight and soaring — this is known as wingtip fingering or slotting. The individual feathers function as winglets that reduce induced drag by vertically spreading the trailing vortex into smaller vortices, following the principle described by Graham (1934) and later formalised by Tucker (1993, 1995).
\[ D_\mathrm{ind} = \frac{L^2}{\tfrac{1}{2}\rho v^2 \pi b_\mathrm{eff}^2 e}, \qquad b_\mathrm{eff} = b\sqrt{1 + \Delta_\mathrm{slot}} \]
\( \Delta_\mathrm{slot} \approx 0.12\text{--}0.18 \) for a golden eagle; equivalent to ~12% reduction in induced drag at thermal-soaring speeds.
Wing-shape classification
Norberg (1990) introduced a two-axis classification of avian wings using aspect ratio vs. wing loading. In this scheme eagles occupy the upper-central region (high AR, medium WL) — ideal for dynamic soaring and thermal soaring. Short-winged forest hawks (Accipiter) sit in the lower-left (low AR, low WL), which maximises manoeuvrability in cluttered habitat.
Norberg wing-morphology space
2. Pennycuick’s Flight-Power Curve
The mechanical power needed for horizontal flapping flight has three components: induced, parasite, and profile. Pennycuick (1975, 1989) combined lifting-line theory with Rayleigh body drag and empirical profile-power scaling to obtain:
\[ P(v) = \underbrace{\frac{2k\,(mg)^2}{\rho\,v\,\pi\,b^2}}_{P_\mathrm{ind}} + \underbrace{\tfrac{1}{2}\rho\,v^3\,S_b\,C_{D,b}}_{P_\mathrm{par}} + \underbrace{1.15\,P_\mathrm{ind}^{\min}}_{P_\mathrm{pro}} \]
The minimum-power speed \( V_{mp} \) is obtained by \( dP/dv = 0 \):
\[ V_{mp} = \left(\frac{2k\,(mg)^2}{3\,\rho^{2}\,\pi\,b^2\,\tfrac{1}{2}S_b\,C_{D,b}}\right)^{1/4} \approx \left(\frac{W}{S}\right)^{1/2}\cdot C_1 \]
The maximum-range speed \( V_{mr} \) is determined by the tangent from the origin to \( P(v) \), equivalent to minimising cost of transport \( C_t = P/(mg v) \):
\[ V_{mr} = \left.\frac{P(v)}{v}\right|_{\min} \;\Rightarrow\; V_{mr} \approx 1.32\, V_{mp} \]
For an Aquila with m = 4.5 kg, b = 2.1 m, S = 0.65 m2: Vmp ≈ 11 m/s, Vmr ≈ 15 m/s (40–55 km/h).
Glide polar and L/D
In a steady glide, the drag polar \( C_D = C_{D,0} + C_L^2/(\pi\,AR\,e) \) yields the sink rate:
\[ V_s = V_\infty \cdot \frac{C_D}{C_L}, \qquad V_\infty = \sqrt{\frac{2\,mg}{\rho\,S\,C_L}} \]
The best glide ratio for a golden eagle is \( (L/D)_\mathrm{max} \approx 10 \) (Tucker 1987), compared to 14 for a turkey vulture (lighter WL, higher AR) and 23 for the albatross (Diomedea exulans). Wing morphing in eagles permits a remarkably wide CL range via dihedral, twist, and alula deployment, so the bird can maintain high L/D at speeds varying by a factor of three.
Thermal vs. ridge soaring
Static soaring uses deflected wind over ridges; the bird maintains altitude if the vertical wind \( w_0 \ge V_s \). Dynamic soaring extracts energy from wind-shear profiles, principally \( dU/dz \ne 0 \) over the sea — the technique that albatrosses perfected but that eagles employ only transiently. Most eagle long-distance travel uses thermal soaring: the bird spirals inside convective thermal cores at \( V_s < w(r) \) and glides between them. MacCready theory gives the optimal inter-thermal speed:
\[ V_{opt} : \left.\frac{dV_s}{dV}\right|_{V_{opt}} = \frac{V_s(V_{opt}) + w_c}{V_{opt}} \]
where \( w_c \) is the next thermal’s expected climb rate. Eagles implement this strategy with remarkable accuracy (Sherub et al. 2016 GPS tracks).
3. The High-Speed Stoop
The stoop dive is the signature attack of large raptors. A golden eagle (Aquila chrysaetos) folds its wings into a teardrop shape and accelerates under gravity toward a prey. Tucker (1998) instrumented a trained golden eagle named “Savannah” and measured a peak speed of 320 km/h (89 m/s) during a ~900 m dive. Higher unverified estimates (up to 240 mph in peregrine falcons) come from radar and GPS logger studies.
The terminal speed in a steep dive is set by the balance of gravity and drag:
\[ v_\mathrm{term} = \sqrt{\frac{2\,m\,g}{\rho\,(C_D A)}} \]
Extended-wing CDA ≈ 0.02 m2 gives 95 km/h; teardrop CDA ≈ 0.0025 m2 gives 270 km/h; a fully optimised dive with favourable air density yields 320+ km/h.
Three-phase wing kinematics
- Acquisition phase: bird identifies prey and enters shallow dive with wings partially extended for targeting control.
- Acceleration phase: progressive wing-fold reduces CDA and permits rapid velocity gain. The bird adjusts pitch control via tail fanning and minimal wing actuation.
- Terminal/strike phase: teardrop configuration with primaries tight against body; feet deployed ~0.5 s before contact; wings flared to brake only if overshoot risk.
Equation of motion (planar)
\[ m\,\ddot{\mathbf{r}} = m\,\mathbf{g} - \tfrac{1}{2}\rho\,|\mathbf{v}|\,\mathbf{v}\,C_D(t)\,A(t) \]
Drag area CDA(t) is piecewise constant corresponding to the three wing-fold phases. Numerical integration (see Sim 2) reproduces observed field velocity profiles within 8%.
Rufous-chested falcon comparison
The rufous-chested falcon (Falco deiroleucus, neotropical) achieves similar stoop speeds despite being 5× lighter. Scaling analysis shows that terminal speed scales as \( v_t \propto \sqrt{m/A} \), and because the effective cross-section A scales as \( m^{2/3} \), we obtain \( v_t \propto m^{1/6} \) — a weak dependence, which is why both small falcons and large eagles reach similar stoop peaks.
\[ v_t \propto m^{1/6} \]
3b. Flapping Kinematics and Unsteady Aerodynamics
Pennycuick’s equation is a time-averaged statement; the actual wing performs unsteady motion producing lift via leading-edge vortices (LEVs), wake capture, and added-mass effects. The Strouhal number \( St = f\,A_\mathrm{wing}/U \) organises the flapping-flight design space; cruising birds converge to \( St \approx 0.2\text{--}0.4 \)(Taylor et al. 2003).
\[ St = \frac{f\,A}{U},\qquad 0.2 \le St \le 0.4 \]
For a golden eagle at cruise: f ≈ 2.5 Hz, A ≈ 1.4 m (peak-to-peak vertical wingtip displacement), U ≈ 15 m/s giving St ≈ 0.23.
The instantaneous lift coefficient during a flap stroke may exceed the steady-state stall value by 30–50%, due to the dynamic delay of leading-edge separation (Dickinson 1999). Eagles exploit this during take-off and slow-flight capture but rely almost entirely on quasi-steady lift during cruise.
Reynolds number and flow regime
At cruise, the chord-based Reynolds number is:
\[ Re = \frac{\rho\,U\,\bar{c}}{\mu} \approx \frac{1.1 \cdot 15 \cdot 0.31}{1.8 \times 10^{-5}} \approx 2.8 \times 10^5 \]
This sits in the transitional regime where laminar separation bubbles affect the pressure distribution — wing camber, covert feather roughness, and the alula all act to promote favourable early transition.
4. Muscle Power and Energetics
The flight muscles (pectoralis + supracoracoideus) comprise roughly 15–17% of body mass in large eagles (vs. 25% in pigeons). Specific muscle power \( P_\mathrm{musc} \approx 100\text{--}200 \) W/kg gives a maximum sustainable power envelope of ~160 W for a 4.5 kg golden eagle. Because \( P(V_{mp}) \approx 80 \) W, the eagle operates at \( \sim 50\% \) of peak aerobic muscle output during cruise — comfortably sustainable for hours.
Oxygen supply and mitochondrial density
Eagle flight muscle shows exceptional aerobic capacity: myoglobin content \( \sim 20 \) mg/g, mitochondrial volume density \( \sim 30\% \) (compared to 2–5% in human vastus). This enables sustained flapping at high altitude — golden eagles have been tracked at 6,500 m over the Himalayas with an oxygen partial pressure \( \sim 50\% \) of sea level.
Metabolic cost scaling
\[ P_\mathrm{met} = \alpha\, m^{0.75}\,(1 + f\,v/V_{mp}) \]
With allometric exponent 0.75 (Kleiber) on basal, and a speed-dependent correction factor f ≈ 1.5 for flapping flight.
Using these parameters, the range of a golden eagle with 500 g of stored fat is roughly 1100 km without feeding — well matched to the longest observed non-stop migrations (e.g. Steller’s sea eagle over the Bering Sea).
Heart output and cardiovascular integration
Raptor heart mass scales as \( m_\mathrm{heart} \approx 0.01\,m \) (Else & Hulbert 1985), double the mammalian benchmark. Cardiac output at peak flight \( Q = f_h \, V_s \) reaches 30–40 L/min for a 4.5 kg golden eagle (heart rate ~350 bpm at peak, stroke volume ~100 mL). This supports the muscular oxygen demand of sustained flapping flight at altitude, discussed in more detail in Module 5 (migration physiology).
5. Attack Aerodynamics and Pullout
At the end of a stoop, the eagle must decelerate and match prey velocity. The pullout manoeuvre imposes sustained load factor \( n = L/(mg) \); measured values reach n ≈ 8–12 in large raptors over a 0.4 s arc, corresponding to centripetal accelerations of 80–120 m/s2.
\[ R_\mathrm{turn} = \frac{v^2}{g\,(n^2-1)^{1/2}},\qquad n_\mathrm{max} = \frac{C_{L,\max}\,\rho\,v^2\,S}{2\,m\,g} \]
At v = 80 m/s with \( C_{L,\max} \approx 1.4 \), the maximum load factor exceeds 25 — but the bird is structurally limited (bones & thorax) to roughly 12 g before injury (Warrick & Dial 1998).
Tucker (2000) documented that peregrine falcons stoop along a curved path that keeps the prey image on the deep fovea, rather than a straight line — trading ~7% extra path length for a factor of 2 reduction in visual slewing. The golden eagle uses a modified variant of this tactic, deflecting during the final phase so that the retinal image sits on the central deep fovea.
Curved-path stoop geometry (Tucker 2000)
The curved dive increases path length by ~7% but keeps the prey image fixed on the deep fovea, eliminating the visual slewing that would arise from a straight-line dive with the eye cocked sideways. See Module 2 for the retinal geometry.
Grip-strike transfer of momentum
At contact, the eagle must convert its kinetic energy into either prey-kill impulse or recoverable braking force. For a 4.5 kg eagle striking a 2 kg rabbit at 30 m/s relative velocity:
\[ J_\mathrm{strike} = m_\mathrm{eag} \Delta v \approx 4.5\cdot 30 = 135\,\mathrm{N\,s} \]
spread over ~30 ms contact time, peak normal force on talons ~4.5 kN — bearing resemblance to the pressure loads discussed in Module 3 (talon biomechanics).
Simulation 1: Pennycuick Power Curve & Glide Polar
Full Pennycuick 1975 decomposition of flight power into induced, parasite and profile components for a 4.5 kg golden eagle, with explicit extraction of \( V_{mp} \), \( V_{mr} \) and the glide polar yielding best L/D.
Click Run to execute the Python code
Code will be executed with Python 3 on the server
Simulation 2: Stooping Dive with Variable Wing-Fold
Numerical integration of Newton’s equations with a piecewise CDA(t) corresponding to extended, partial, and teardrop wing-fold geometries. Output reproduces Tucker’s 1998 320 km/h measurement and quantifies the energy partitioning between kinetic, potential and dissipated contributions.
Click Run to execute the Python code
Code will be executed with Python 3 on the server
Key References
• Pennycuick, C. J. (1975). “Mechanics of flight.” In Avian Biology, volume 5, Academic Press, pp. 1–75.
• Pennycuick, C. J. (1989). Bird Flight Performance: A Practical Calculation Manual. Oxford University Press.
• Pennycuick, C. J. (2008). Modelling the Flying Bird. Elsevier.
• Tucker, V. A. (1987). “Gliding birds: the effect of variable wing span.” Journal of Experimental Biology, 133, 33–58.
• Tucker, V. A. (1993). “Gliding birds: reduction of induced drag by wing-tip slots between the primary feathers.” Journal of Experimental Biology, 180, 285–310.
• Tucker, V. A. (1998). “Gliding flight: speed and acceleration of ideal falcons during diving and pull out.” Journal of Experimental Biology, 201, 403–414.
• Graham, R. R. (1934). “The silent flight of owls.” Journal of the Royal Aeronautical Society, 38, 837–843.
• Norberg, U. M. (1990). Vertebrate Flight: Mechanics, Physiology, Morphology, Ecology and Evolution. Springer.
• Sherub, S., Fiedler, W., Duriez, O. & Wikelski, M. (2016). “Bio-logging, new technologies to study conservation physiology on the move: a case study on annual survival of Himalayan vultures.” Journal of Comparative Physiology A, 203, 531–542.
• Watanabe, Y. Y. (2016). “Flight mode affects allometry of migration range in birds.” Ecology Letters, 19, 907–914.