ABSTRACT
For gap-crossing agility, arboreal animals require the ability to stabilize dynamic landings on branches. Despite lacking a prehensile grip, squirrels achieve stable landings using a palmar grasp. We investigated the landing dynamics of free-ranging fox squirrels (Sciurus niger) to uncover strategies for stable, above-branch landings. Using high-speed video and force-torque measurements in the sagittal plane, we quantified landing kinetics across gap distances. Squirrels rapidly managed >80% of the landing energy with their forelimbs. With larger gaps, peak leg force and foot torque increased. Alignment between forelimbs, velocity and force also increased, likely reducing joint moment. We tested control hypotheses based on an extensible pendulum model used in a physical, hopping robot named Salto. Squirrels stabilized off-target landings by modulating leg force and foot torque. To correct for undershooting, squirrels generated pull-up torques and reduced leg force. For overshooting, squirrels generated braking torques and increased leg force. Embodying control principles in leg and foot design can enable stable landings in sparse environments for animals and robots alike, even those lacking prehensile grasps.
INTRODUCTION
Grasping capabilities can affect locomotion, feeding, social interactions and reproductive behaviors of many animals, including all tetrapod clades (Pouydebat et al., 2023; Sustaita et al., 2013). Among locomotor behaviors in tree canopies, gap crossing onto narrow and sparse branches stands out as a common dynamic activity that often requires high-impact contact and stabilization (Graham and Socha, 2020).
To stabilize dynamic motion, primates can use highly effective prehensile grasps, a type of grasp where appendages wrap around an object. They employ opposable digits such as the pollex (thumb) and hallux (big toe) to navigate arboreal environments (Toussaint et al., 2020; Nyakatura, 2019). While primates have been extensively studied, other tree-dwelling mammals such as squirrels also demonstrate remarkable adaptations for arboreal locomotion. Research on tree squirrel biomechanics has primarily focused on quantifying locomotor dynamics along branches of various inclines and sizes (Dunham et al., 2019; Essner, 2002; Wölfer et al., 2021; Orkin and Pontzer, 2011; Youlatos, 1999; Schmidt and Fischer, 2011).
Young and Chadwell (2020) directly compared the locomotor mechanics along a branch in a sciurid rodent and two platyrrhine primates. These species were selected to represent different evolutionary stages of grasping adaptations in primates, from primitive to more evolved. When presented with branches of different widths, squirrel monkeys, with their superior prehensile abilities, exhibited minimal kinematic adjustments in gait, speed, duty factor and peak impact force. In contrast, marmosets demonstrated moderate adjustments, whereas squirrels, which lack a prehensile grasp, required the greatest adjustments. Interestingly, squirrels were also characterized by the lowest values of peak rolling angular momentum over a stride. Assuming limited grasping ability to apply torques to modulate roll, the results suggest that squirrels effectively used dynamic stability as a control strategy for arboreal balance. In the present study, we hypothesize that foot friction plays a significant role in the landing and pitch balancing mechanism of squirrels, indicating a possible nuanced interplay between dynamic stability and the application of torque about a branch.
To better understand the challenges presented by arboreal environments, we can draw insight from robotic grasp taxonomies, which offer precise analytical frameworks for examining the role of foot friction. Grasp types have been defined by an object's shape and size relative to the gripper and its configuration relative to the object, which may vary in the degree of digit contact and wrapping (Cutkosky and Howe, 1990). Feix et al. (2015) define a power palmar grasp with an adducted thumb as a grasp type that relies on opposability afforded by palm and limited digit wrapping. Remarkably, squirrels employ only nonprehensile palmar grasps yet execute highly precise and stable, above-branch landings, all while maintaining readiness for dynamic parkour-like movements.
Studies on squirrel across-branch locomotion have examined jumping (Hunt et al., 2021; Boulinguez-Ambroise et al., 2023) and gliding (Byrnes et al., 2008; Paskins et al., 2007; Scheibe et al., 2007) dynamics. Hunt et al. (2021) studied free-ranging fox squirrels jumping from simulated branches to narrow perches, revealing sophisticated decision-making processes. Squirrels balanced trade-offs between branch-bending compliance and gap distance when choosing launch points and demonstrated the ability to modify impulse generation upon repeated jumps from unexpectedly compliant beams. Impressively, in over 100 trials, squirrels never missed or fell.
When far off-target, squirrels exhibited fail-safe and fault tolerant landings, skillfully swinging their center of gravity (CG) under or over the target branch (Fig. 1). As spectacular as these landings are, unstable landings can impede rapid responses, increasing susceptibility to predation. Despite squirrels using nonprehensile palmar grasps, most branch landings were direct (Fig. 1, three middle sectors) and avoided rotations of the center of gravity under (Fig. 1, darker blue sectors) or over (Fig. 1, darker red sectors) the branch, allowing squirrels to maintain their posture above the branch for subsequent maneuvers. To further understand these abilities, we aimed to quantify the dynamics of direct landings, testing hypotheses of control stabilizing mechanisms that correct for undershooting (Fig. 1, lightest blue sector) and overshooting (lightest red sector) when deviating from nominal (white sector).
Types of landing in squirrels. As defined by Hunt et al. (2021), squirrels can land in a variety of ways. Squirrels may swing under a target branch (dark blue), sometimes being able to make contact with only their front feet (reduced foot contact). Squirrels may also swing over a branch (dark red), sometimes only being able to make contact with their hind feet. Squirrels may also land directly on a branch (middle sectors) and avoid center of gravity (CG) inversion altogether, which is critical for taking another leap if necessary. Within the envelope of direct landings, squirrels may land nominally (zero landing error), undershoot (negative landing error) or overshoot (positive landing error) their target.
Types of landing in squirrels. As defined by Hunt et al. (2021), squirrels can land in a variety of ways. Squirrels may swing under a target branch (dark blue), sometimes being able to make contact with only their front feet (reduced foot contact). Squirrels may also swing over a branch (dark red), sometimes only being able to make contact with their hind feet. Squirrels may also land directly on a branch (middle sectors) and avoid center of gravity (CG) inversion altogether, which is critical for taking another leap if necessary. Within the envelope of direct landings, squirrels may land nominally (zero landing error), undershoot (negative landing error) or overshoot (positive landing error) their target.
To quantify nonprehensile, high-impact, dynamic landing in free-ranging squirrels, we conducted experiments in a eucalyptus forest using artificial branches. We measured both touchdown state and forces in the sagittal plane at various gap distances. We used high-speed video to measure the landing kinematics and we designed a force-torque sensor apparatus that could be transported to the field to enable the first landing kinetic measurements on a horizontal branch (Fig. 2A–C). This approach allowed us to test control hypotheses for direct landings postulated by Yim et al. (2025) for the physical model/robot Salto (Haldane et al., 2016; Yim and Fearing, 2018; Yim et al., 2020). We modeled the CG to front foot contact as a compressible virtual leg comparable to a simple extensible pendulum system (Fig. 2E), which provides a conceptual framework for interpreting the results.
Field apparatus and key kinematic quantities of Sciurus niger landing. (A) Mobile field apparatus in the eucalyptus forest. High-speed video recordings and force-torque measurements were made during voluntary landings using the cart supporting the aluminum frame and F/T sensor, a ramp, lights, high-speed camera and a laptop for data acquisition. (B) Diagram showing the take-off branch, the gap distance jumped, the force/torque (F/T) sensor (with lateral view in box) and high-speed camera. The take-off rod was attached to a linear rail, which allowed for variable gap distances of 50, 75 and 100 cm. (C) Center of gravity (CG) is graphically calculated by fitting a parabola between three points as in Hunt et al. (2021): nose, tailbase and the midpoint of the ventral-dorsal line (red dashed), which approximately bisects the tailbase–nose line (blue dashed). CG position and velocity were extracted to calculate CG angular momentum (Lcg) and landing error (e). The distance, h, between F/T sensor interface and the branch axis was h=40 mm. (D) Landing error (e) was calculated as the projected vertical offset of the extrapolated aerial trajectory (see Hunt et al., 2021). (E) Body angular momentum (Lb) is a function of body moment of inertia (Ib) and body pitch rate (ω). The body is modeled as a cylinder of body length (BL), diameter (body width, BW) and at the CG position. All variables are defined in Table 1.
Field apparatus and key kinematic quantities of Sciurus niger landing. (A) Mobile field apparatus in the eucalyptus forest. High-speed video recordings and force-torque measurements were made during voluntary landings using the cart supporting the aluminum frame and F/T sensor, a ramp, lights, high-speed camera and a laptop for data acquisition. (B) Diagram showing the take-off branch, the gap distance jumped, the force/torque (F/T) sensor (with lateral view in box) and high-speed camera. The take-off rod was attached to a linear rail, which allowed for variable gap distances of 50, 75 and 100 cm. (C) Center of gravity (CG) is graphically calculated by fitting a parabola between three points as in Hunt et al. (2021): nose, tailbase and the midpoint of the ventral-dorsal line (red dashed), which approximately bisects the tailbase–nose line (blue dashed). CG position and velocity were extracted to calculate CG angular momentum (Lcg) and landing error (e). The distance, h, between F/T sensor interface and the branch axis was h=40 mm. (D) Landing error (e) was calculated as the projected vertical offset of the extrapolated aerial trajectory (see Hunt et al., 2021). (E) Body angular momentum (Lb) is a function of body moment of inertia (Ib) and body pitch rate (ω). The body is modeled as a cylinder of body length (BL), diameter (body width, BW) and at the CG position. All variables are defined in Table 1.
We hypothesized that a nominal, direct landing (Fig. 1, white sector) could be attained by squirrels using only their angular momentum to swing to a stable, balanced position. Far more likely, however, we postulated that squirrels would undershoot or overshoot the target branch (Fig. 3A, Movie 1). Stabilizing off-target landings by controlling rotation about the branch could be accomplished by modulating leg force and foot torque (Fig. 3B). Specifically, control of front leg force could correct for overshooting by generating more leg-braking force, whereas undershooting could be stabilized by generating less leg-braking force. When overshooting, squirrels can also generate braking torques in the direction opposite to their forward motion, countering large landing errors and excessive rotation. Conversely, for undershooting, they can produce a pull-up torque in the direction of their forward motion, correcting for negative landing errors.
Balanced landing control hypotheses for undershooting and overshooting trajectories. (A) Snapshots of overshooting and undershooting at three different time intervals: touchdown, 20 ms and 50 ms. At touchdown, the white arrows begin at the CG and represent the velocity vector relative to the branch (see Movie 1). (B) Force-torque control hypotheses derived from Yim et al. (2025) exemplify the interplay of leg force and foot torque control for balanced, above-branch landings. When overshooting, squirrels apply greater braking forces and torque. When undershooting, squirrels apply lower braking forces and pull-up torques. The green dot represents the branch, and the dashed line, the branch level. The light blue line shows the trajectory prior to touchdown and the green dashed arrow, the projected trajectory. The gray arrow shows the touchdown velocity vector, the white arrow, the branch reaction force, and the magenta arrow, the branch reaction torque direction.
Balanced landing control hypotheses for undershooting and overshooting trajectories. (A) Snapshots of overshooting and undershooting at three different time intervals: touchdown, 20 ms and 50 ms. At touchdown, the white arrows begin at the CG and represent the velocity vector relative to the branch (see Movie 1). (B) Force-torque control hypotheses derived from Yim et al. (2025) exemplify the interplay of leg force and foot torque control for balanced, above-branch landings. When overshooting, squirrels apply greater braking forces and torque. When undershooting, squirrels apply lower braking forces and pull-up torques. The green dot represents the branch, and the dashed line, the branch level. The light blue line shows the trajectory prior to touchdown and the green dashed arrow, the projected trajectory. The gray arrow shows the touchdown velocity vector, the white arrow, the branch reaction force, and the magenta arrow, the branch reaction torque direction.
In addition to quantifying and characterizing squirrel branch landings, our study also aims to stimulate further research on the role of embodied control in foot and leg design for dynamic, high-impact interactions. Understanding the biomechanics of the elaborate morphology of squirrel foot and toe design can lead to novel use of metamaterials (e.g. Jin et al., 2023). Interdisciplinary collaboration with engineers could provide biological inspiration for the development of the next generation of legs and feet for agile legged robots, as well as nonprehensile robotic manipulators. These findings show that examining dynamic, high-impact landings and stabilization with only nonprehensile, palmar grasps can further advance the field of grasping and manipulation for both animals and engineered systems.
MATERIALS AND METHODS
To measure squirrel landing dynamics, we designed experiments to capture the kinematics, forces, and torques using an instrumented branch or rod. The experiments involved three key components: (1) the acquisition and training of free-ranging research animals; (2) the utilization of a mobile instrumented apparatus; and (3) a systematic data collection and processing procedure for subsequent analysis. The ethical treatment of animals and adherence to protocols were ensured through the approval of the University of California, Berkeley's Animal Care and Use Committee (ACUC) Protocol #AUP-2018-06-11201-1 and the California Department of Fish and Wildlife Nongame Wildlife Program.
Animal preparation
Sciurus niger Linnaeus 1758 are free-ranging arboreal squirrels known for their adept navigation in tree canopies and have adapted to urban environments. Their agility and accessibility make S. niger an ideal candidate for studying and characterizing arboreal locomotion (Hunt et al., 2021). Here, we focus on the measurement of landing forces, rod torque and high-speed video involving four female fox squirrels (729±63 g; mean±s.d.). To uniquely identify each individual, a non-toxic fur dye (Nyanzol D) was prepared and dispersed on the squirrels' fur.
Introducing free-ranging individuals to a novel apparatus required training squirrels with a shaping paradigm using peanuts as positive reinforcement for approximately 3–5 weeks. Individuals were introduced to the setup at different times over a span of 4 months based on their voluntary willingness to cooperate. After acclimation, individuals were trained to follow a feeding stick and execute jumps from a non-instrumented take-off birch perch to an instrumented target birch perch [both 1.91 cm (0.75″) diameter]. To ensure consistent behavior meeting our operational definition of direct landing, each individual underwent a training period consisting of at least five successful landings on the instrumented perch, before data collection. A successful landing was defined as a direct landing in which foot contact is restricted to the rod only without swinging under or over the landing target (Fig. 1, center sectors outlined in bold). Training was administered to each individual at each gap distance and was then followed by five more landings for which data were collected.
Field force, torque, and kinematics apparatus
We designed and implemented an experimental setup to simultaneously capture three-dimensional force/torque data and 2D high-speed video of squirrels landing on artificial and instrumented branches. The force measurement device (Fig. 2A,B) utilized a 6-axis force/torque (F/T) transducer (ATI Mini45 with SI-145-5 calibration) and data acquisition system (NI USB-6210 DAQ). The F/T sensor was connected to the DAQ and a power supply for 3D force data streaming. The DAQ was connected to a computer (Windows 11 Laptop) for data recording through NI LabVIEW. The hardware selection allowed for measurements with force resolution of 62.5 mN in all three axes, a torque resolution of 1.3 N mm in the X-axis and a sampling rate of 3000 Hz.
The mechanical structure of the apparatus was designed for high stiffness but also to be lightweight and compact for ease of transport into the field (Fig. 2A). To minimize vibrations in the F/T signal, an aluminium truss was attached to the base-side of the F/T sensor. An artificial branch was attached to the topside of the sensor, and its offset from the sensor interface was minimized to maximize the cantilever's natural frequency, which aided in minimizing signal loss when filtering the signal. Outdoors, the ground can be uneven, so two measurement levels were attached to the setup to ensure the sensor was as horizontal as possible. The take-off branch was attached to a rail, which enabled quick adjustment of gap distance, from 50 cm to 100 cm, where 100 cm is approximately four squirrel body lengths.
As shown in Fig. 2A, the instrumented setup was transported to a nearby eucalyptus grove. A ramp guided squirrels to the take-off branch substrate. A high-speed camera (Vision Research, Phantom v. 10.0) sat level and approximately 2.5 m away from the subject. Flood lights (Lowel Tota-Light Tungsten) provided a non-flickering light source. The high-speed camera was set to record at 500 frames s−1 with a resolution of 1920×1080 pixels. Force and video data were electronically synchronized through the Vision Research GUI.
Instruments of the apparatus (i.e. camera, flood lights, laptop, DAQ) were first connected to an outdoor power outlet. The F/T sensor and camera were placed and then leveled such that the landing substrate and squirrel's body were visible in the camera frame. Given the brevity of the landing event, the instruments were armed to end-trigger recordings. When available, free-ranging individuals were guided towards the setup with feeding sticks. A researcher prompted an individual to climb up the ramp and to cross the gap by leaping between rods or perches. For a given trial, another researcher end-triggered force and video data as soon as the squirrel came into contact with the F/T sensor. Data were then written to the laptop. After five trials per individual, data collection was terminated for the given experimental condition.
Feature extraction
Measured and calculated times, force/torque and kinematic variables used for characterizing each landing sequence
Variable . | Name . | Unit . | Description . |
---|---|---|---|
Timing | |||
t1 | 1st peak time | ms | First peak force |
thf | Hind feet touchdown | ms | Second contact |
T | Settling time | ms | End of landing |
Force, torque, impulse | |||
Fh | Horizontal force magnitude | N | y-axis force |
Fv | Vertical force magnitude | N | z-axis force |
τx | x-axis torque magnitude | N mm | Measured from F/T sensor |
Fnet | Net force magnitude | N | Camera-view force ![]() |
θF | Force angle w.r.t. horizontal | deg | atan2(Fv, Fh) |
τrod | Rod-axis torque magnitude | N mm | τrod=τx−Fhh, where h is defined in Fig. 2C |
J0 | Impulse from front legs | N s | ![]() |
JT | Total landing impulse | N s | ![]() |
Touchdown kinematics | |||
r | Virtual leg length from rod to CG | mm | Distance r in Fig. 2C |
v | CG speed | m s−1 | Speed v in Fig. 2C |
θcg | CG touchdown angle w.r.t. horizontal | deg | Angle θcg in Fig. 2C |
θv | Velocity angle w.r.t. horizontal | deg | Angle θv in Fig. 2C |
BL/BW | Body length and width | mm | Dimensions of squirrel body as a cylinder |
Ib | Body moment of inertia | m² | Ib=f (BW, BL, r, θcg), see Materials and Methods |
γ | Body pitch | deg | Angle w.r.t. horizontal of blue dashed line in Fig. 2C |
ω | Body pitch rate | deg s−1 | ![]() |
e | Landing error | mm | Distance e in Fig. 2D |
Lcg | CG angular momentum | m² s−1 | Lcg=v sin(θcg -θv)r in Fig. 2C |
Lb | Body angular momentum | m² s−1 | Lb=Ibω |
Ltotal | Total angular momentum | m² s−1 | Ltotal=Lcg+Lb |
Variable . | Name . | Unit . | Description . |
---|---|---|---|
Timing | |||
t1 | 1st peak time | ms | First peak force |
thf | Hind feet touchdown | ms | Second contact |
T | Settling time | ms | End of landing |
Force, torque, impulse | |||
Fh | Horizontal force magnitude | N | y-axis force |
Fv | Vertical force magnitude | N | z-axis force |
τx | x-axis torque magnitude | N mm | Measured from F/T sensor |
Fnet | Net force magnitude | N | Camera-view force ![]() |
θF | Force angle w.r.t. horizontal | deg | atan2(Fv, Fh) |
τrod | Rod-axis torque magnitude | N mm | τrod=τx−Fhh, where h is defined in Fig. 2C |
J0 | Impulse from front legs | N s | ![]() |
JT | Total landing impulse | N s | ![]() |
Touchdown kinematics | |||
r | Virtual leg length from rod to CG | mm | Distance r in Fig. 2C |
v | CG speed | m s−1 | Speed v in Fig. 2C |
θcg | CG touchdown angle w.r.t. horizontal | deg | Angle θcg in Fig. 2C |
θv | Velocity angle w.r.t. horizontal | deg | Angle θv in Fig. 2C |
BL/BW | Body length and width | mm | Dimensions of squirrel body as a cylinder |
Ib | Body moment of inertia | m² | Ib=f (BW, BL, r, θcg), see Materials and Methods |
γ | Body pitch | deg | Angle w.r.t. horizontal of blue dashed line in Fig. 2C |
ω | Body pitch rate | deg s−1 | ![]() |
e | Landing error | mm | Distance e in Fig. 2D |
Lcg | CG angular momentum | m² s−1 | Lcg=v sin(θcg -θv)r in Fig. 2C |
Lb | Body angular momentum | m² s−1 | Lb=Ibω |
Ltotal | Total angular momentum | m² s−1 | Ltotal=Lcg+Lb |
Statistical analyses
Over eight experimental sessions, squirrels performed 60 successful landings. We focused our analysis on force and high-speed video data of n=60 successful landings, five trials for four individuals at three gap conditions. Some trials were excluded for one of two reasons: (1) The camera view did not capture the squirrel's tail base at touchdown, so CG kinematics at touchdown could not be extracted (n=11), or (2) the squirrel severely undershot or overshot the target branch, resulting in an unstable landing, appendages coming into contact with other parts of the setup, and therefore invalid force-torque measurements (n=2). Then, for extraction of the landing state from high-speed video data, we analyzed n=49 landings. Means and standard deviations are reported for metrics across all trials for each gap condition. P-values and F-statistics are reported from one-way repeated-measures ANOVA comparisons across gap distances. P-values and t-statistics are also reported for the linear mixed-effects models controlling for gap and individual to show predictive power of landing error and angular momentum. Data were analyzed using MATLAB statistical software tools.
RESULTS
Landing dynamics as a function of time
Kinematics of landing
Throughout the experiments, squirrels crossed gaps following projectile motion trajectories. Forces normal to the camera view (X-axis of the load cell) were negligible compared with forces in the squirrel's sagittal plane (Y–Z plane of load cell). Squirrels exhibited a diverse range of landing touchdown states that varied in CG position and velocity with respect to the branch. In all landings, squirrels always touched down with their front feet first, quickly followed by their hind feet (Fig. 4A). Upon landing, squirrels actively controlled their body pose by rotating their joints so that their hind legs could contact the branch. Contact with their hind feet was required to reach the stable, above-branch, perched state. In a typical sequence for a given gap and individual, squirrels prepared for touchdown by extending their forelegs and front feet. The extended feet first engaged with the branch to create a reliable anchor. Then, the front limbs compress as kinetic energy decreases. A force peak always followed front feet touchdown. The first peak force induced by the front limbs was always greater than the second one, which was induced after hind feet touchdown. At the end of the landing at t=T, the hind feet provided friction for static balance while the front feet either stayed in contact or detached.
Branch landing reaction force, maximum force angle and rod torque as a function of time for a given gap and individual (n=5). (A) Mean net force and force components plotted over time. Shaded bands show one standard deviation. Squirrels showed maximum force peaks during front feet touchdown and a secondary smaller peak after hind feet touchdown. Vertical dashed lines indicate critical landing times (t1, time to peak force; thf, hind feet touchdown; T, landing completion time). Front feet (t0) or hind feet (thf) touchdown time is defined as the frame in which the feet first come into contact with the branch. Settling time (T) is defined as the time at which Fh reaches 0 N. (B) Net force angle stabilizes at approximately 50 deg and steadily increases to a vertical 90 deg after thf until t=T. (C) Maximum rod torque can be positive and negative throughout the time series, as seen in the magenta shaded area. When t>T, τrod settles to a non-zero value, which corresponds to static torque while perched.
Branch landing reaction force, maximum force angle and rod torque as a function of time for a given gap and individual (n=5). (A) Mean net force and force components plotted over time. Shaded bands show one standard deviation. Squirrels showed maximum force peaks during front feet touchdown and a secondary smaller peak after hind feet touchdown. Vertical dashed lines indicate critical landing times (t1, time to peak force; thf, hind feet touchdown; T, landing completion time). Front feet (t0) or hind feet (thf) touchdown time is defined as the frame in which the feet first come into contact with the branch. Settling time (T) is defined as the time at which Fh reaches 0 N. (B) Net force angle stabilizes at approximately 50 deg and steadily increases to a vertical 90 deg after thf until t=T. (C) Maximum rod torque can be positive and negative throughout the time series, as seen in the magenta shaded area. When t>T, τrod settles to a non-zero value, which corresponds to static torque while perched.
Branch reaction force timing during landing
Landings for a given individual were consistent. A characteristic set of leg reaction force data (n=5 trials) for one individual at the 75 cm gap distance is shown in Fig. 4A with its net force angle in Fig. 4B. Fig. 4C illustrates rod torque over time and highlights the maximum rod torque, which on average was positive (in the pull-up direction). Here, the shaded regions illustrate variability, which is defined as the standard deviation across five trials. We define t0 as t=0 ms, the time at which the front feet first come into contact with the branch. We visually examined high-speed video to extract the frame at which contact occurred. The first force peak occurred at t1, the time at which Fnet reached its global maximum. The hindfoot touchdown time, thf, is defined as the time at which the hind feet first make contact, which was visually noted on the high-speed video. Finally, we estimated the end of the landing, T, as the time at which the horizontal force, Fh, reached and settled to 0 N. Out-of-plane peak forces were ignored, as they were on average 1.34±0.69 N (or 0.19±0.09 BW). This was on average 8.7% of observed Fnet. Snapshots corresponding to t0, thf and T are displayed sequentially in Fig. 4A. Table 2 reports time variables mean and standard deviation for each event at each gap distance tested. Depending on gap distance, squirrels completed landings within T=200–350 ms. Squirrels reached peak forces within 20–65 ms. Hind feet touchdown occurred within 60–140 ms or within 29–43% of T.
Landing variable measurements as a function of gap distanced
. | Gap distance . | . | . | ||
---|---|---|---|---|---|
. | 50 cm . | 75 cm . | 100 cm . | F . | P-value . |
Timing | |||||
t1 (ms) | 45±20 | 28±7.7 | 22±9 | 30.8 | 8.04e–7 *** |
thf (ms) | 124±15 | 95±12 | 68±10 | 353.9 | 2.30e–26 *** |
T (ms) | 313±61 | 262±27 | 232±59 | 30.6 | 7.97e–7 *** |
Force and torque | |||||
Fnet/mg (BW) | 2.09±0.19 | 3.25±0.30 | 4.31±0.49 | 488 | 2.51e–29 *** |
θF (deg) | 51.7±3.4 | 47.5±4.4 | 46.1±5.2 | 38.3 | 7.47e–8 *** |
τrod/m (N mm kg−1) | 54±58 | 95±97 | 136±169 | 5.0 | 0.029 * |
|τrod|/m (N mm kg−1) | 72±31 | 114±73 | 191±97 | 31.3 | 6.86e–7 *** |
Impulse | |||||
J0 (N s) | 1.09±0.29 | 1.64±0.34 | 1.77±0.32 | 81.4 | 1.65e–12 *** |
JT (N s) | 1.67±0.30 | 2.33±0.35 | 2.77±0.26 | 324.7 | 5.66e–25 *** |
α (%) | 65.4±9.7 | 70.3±8.3 | 64.1±11.3 | 0.3 | 0.58 |
β (%) | 87.1±7.6 | 90.5±5.5 | 85.9±9.0 | 0.5 | 0.48 |
Touchdown kinematics | |||||
r (mm) | 205±9 | 208±10 | 208±7 | 5.7 | 0.021 * |
θcg (deg) | 34.9±4.0 | 34.4±4.8 | 37.7±5.6 | 2.7 | 0.10 |
θv (deg) | 23.1±4.8 | 28.4±4.9 | 34.1±4.7 | 145.4 | 5.46e–16 *** |
v (m s−1) | 1.71±0.18 | 2.29±0.23 | 2.79±0.18 | 251.7 | 1.68e–20 *** |
vr (m s−1) | 1.67±0.17 | 2.28±0.23 | 2.78±0.18 | 277.2 | 2.44e–21 *** |
vθ (m s−1) | 0.35±0.09 | 0.24±0.10 | 0.17±0.15 | 32.6 | 7.45e–7 *** |
BL (mm) | 279.2±16.3 | 270.7±13.1 | 270.0±7.8 | 4.7 | 0.035 * |
γ (deg) | 24.3±3.3 | 26.7±3.7 | 25.4±3.7 | 0.4 | 0.51 |
ω (deg s−1) | −122.4±63.2 | −59.6±87.1 | −25.0±84.2 | 14.0 | 5.06e–4 *** |
Landing error and angular momentum | |||||
e (mm) | −11.7±18.2 | −11.9±15.0 | −9.7±16.5 | 0.6 | 0.43 |
Lcg (m2 s−1) | 0.073±0.020 | 0.051±0.023 | 0.036±0.031 | 30.6 | 1.38e–6 *** |
Lb (m2 s−1) | −0.015±0.009 | −0.007±0.010 | −0.003±0.009 | 15.0 | 3.37e–4 *** |
Ltotal (m2 s−1) | 0.058±0.016 | 0.044±0.023 | 0.033±0.034 | 10.5 | 2.20e–3 ** |
. | Gap distance . | . | . | ||
---|---|---|---|---|---|
. | 50 cm . | 75 cm . | 100 cm . | F . | P-value . |
Timing | |||||
t1 (ms) | 45±20 | 28±7.7 | 22±9 | 30.8 | 8.04e–7 *** |
thf (ms) | 124±15 | 95±12 | 68±10 | 353.9 | 2.30e–26 *** |
T (ms) | 313±61 | 262±27 | 232±59 | 30.6 | 7.97e–7 *** |
Force and torque | |||||
Fnet/mg (BW) | 2.09±0.19 | 3.25±0.30 | 4.31±0.49 | 488 | 2.51e–29 *** |
θF (deg) | 51.7±3.4 | 47.5±4.4 | 46.1±5.2 | 38.3 | 7.47e–8 *** |
τrod/m (N mm kg−1) | 54±58 | 95±97 | 136±169 | 5.0 | 0.029 * |
|τrod|/m (N mm kg−1) | 72±31 | 114±73 | 191±97 | 31.3 | 6.86e–7 *** |
Impulse | |||||
J0 (N s) | 1.09±0.29 | 1.64±0.34 | 1.77±0.32 | 81.4 | 1.65e–12 *** |
JT (N s) | 1.67±0.30 | 2.33±0.35 | 2.77±0.26 | 324.7 | 5.66e–25 *** |
α (%) | 65.4±9.7 | 70.3±8.3 | 64.1±11.3 | 0.3 | 0.58 |
β (%) | 87.1±7.6 | 90.5±5.5 | 85.9±9.0 | 0.5 | 0.48 |
Touchdown kinematics | |||||
r (mm) | 205±9 | 208±10 | 208±7 | 5.7 | 0.021 * |
θcg (deg) | 34.9±4.0 | 34.4±4.8 | 37.7±5.6 | 2.7 | 0.10 |
θv (deg) | 23.1±4.8 | 28.4±4.9 | 34.1±4.7 | 145.4 | 5.46e–16 *** |
v (m s−1) | 1.71±0.18 | 2.29±0.23 | 2.79±0.18 | 251.7 | 1.68e–20 *** |
vr (m s−1) | 1.67±0.17 | 2.28±0.23 | 2.78±0.18 | 277.2 | 2.44e–21 *** |
vθ (m s−1) | 0.35±0.09 | 0.24±0.10 | 0.17±0.15 | 32.6 | 7.45e–7 *** |
BL (mm) | 279.2±16.3 | 270.7±13.1 | 270.0±7.8 | 4.7 | 0.035 * |
γ (deg) | 24.3±3.3 | 26.7±3.7 | 25.4±3.7 | 0.4 | 0.51 |
ω (deg s−1) | −122.4±63.2 | −59.6±87.1 | −25.0±84.2 | 14.0 | 5.06e–4 *** |
Landing error and angular momentum | |||||
e (mm) | −11.7±18.2 | −11.9±15.0 | −9.7±16.5 | 0.6 | 0.43 |
Lcg (m2 s−1) | 0.073±0.020 | 0.051±0.023 | 0.036±0.031 | 30.6 | 1.38e–6 *** |
Lb (m2 s−1) | −0.015±0.009 | −0.007±0.010 | −0.003±0.009 | 15.0 | 3.37e–4 *** |
Ltotal (m2 s−1) | 0.058±0.016 | 0.044±0.023 | 0.033±0.034 | 10.5 | 2.20e–3 ** |
BL, body length; α, forelimb impulse contribution; β, forelimb energy contribution; vr, radial speed; vθ, tangential speed. For all other variables, please refer to Table 1. F statistics and P-values are reported from repeated-measures ANOVA comparisons between experimental conditions. *P<0.05; **P<0.01; ***P<0.001.
Impulse reflecting energy management by forelimbs
Forelimb impulse is calculated as the area under the Fnet curve from t0 to thf (see Table 1). JT is the total impulse induced by the landing event. Impulse can also be calculated as a change in momentum (J=mΔv). We define α as the ratio between forelimb impulse and total impulse, α=Jo/JT (Eqn. 12). We found that α=67±10%, which means on average, 67% of speed was decelerated by the forelimbs alone.
We can also compute β, the ratio between kinetic energy managed by the forelimbs over the total touchdown kinetic energy. Here, we express β in terms of α, as β=2α–α2 (Eqn 15). Using this equation, we calculate that on average, squirrel forelimbs managed 88% of landing kinetic energy, regardless of gap distance. Forelimb energy absorption occurred within 60–130 ms, which on average was 36% of T. Therefore, forelimbs were responsible for managing most of the landing kinetic energy in a fraction of the landing period.
Landing dynamics as a function of gap distance
Gap distance effect on touchdown kinematics
Squirrel touchdown kinematics are defined by the CG's position and velocity at t0. As seen in Table 2, virtual leg length (r) and leg angle (θcg) did not vary significantly across gap distances (P>0.01). However, speed (v) and velocity angle (θv) both increased with gap distance. As gap distance increased, squirrels took higher (one-way repeated measures ANOVA, F1,47=145.4, P<0.001) and faster ballistic trajectories (F1,47=251.7, P<0.001). As gap distance increased, squirrel landing velocity angle, θv, also approached the 45 deg launch angle that tends to maximize horizontal travel for a particular speed.
Touchdown velocity can be broken down into radial and tangential components. Like speed, radial speed (vr) increased with gap distance (F1,47=277.2, P<0.001). Tangential speed (vθ) on the other hand, decreased with gap distance (F1,47=32.6, P<0.001). Squirrel touchdown kinematics are also defined by the body's posture and rotational speed at t0. Body length (BL) and body pitch (γ) were not significantly different across gap distances (P>0.01). However, body pitch rate (ω) decreased in magnitude as gap distance increased (F1,47=14.0, P<0.001).
For all key time points, average time decreased as gap distance increased, pointing towards a consistent force behavior over time. Hindfoot touchdown time (thf) decreased from 124 ms at 50 cm to 68 ms at 100 cm, indicating hind feet touched down faster at longer gap distances (F1,58=353.9, P<0.001; see Table 2). The landings were completed more rapidly the longer the gap distance (F1,58=30.6, P<0.001). The landing sequences themselves were approximately 80 ms shorter at the 100 cm gap compared with the 50 cm gap. Time to peak reaction force also occurred sooner, decreasing from 45 ms to 22 ms (F1,58=30.8, P<0.001).
Gap distance effects on force-torque, force angle, and impulse
Upon touchdown across all trials, squirrels exerted a peak wrench (peak force and torque) on the branch (Table 2). As expected, peak Fnet increased with gap distance more than doubling from 2.09 BW at 50 cm to 4.31 BW at the 100 cm distance (Fig. 5, F1,58=488, P<0.001). Maximum force angle (θF) decreased with increasing gap distance (F1,58=38.3, P<0.001) from 51.7 deg at 50 cm to 46.1 deg at 100 cm (see Table 2), becoming more horizontal at longer gap distances. A lower maximum force angle implies that at longer gap distances, squirrels are relying more on horizontal force rather than vertical force to decelerate their landings.
Grasping wrench space represented by peak foot torque (τrod) and peak leg force (Fnet) at each gap distance. Wrench space is the set of force and torque combinations that permit dynamic landing. Individual points and box plots illustrate data spread with median. Bars represent quartiles. Peak Fnet magnitude and variation tend to increase with gap distance (***P<0.001). Comparisons are made across gaps (***P<0.001). Mean peak τrod is statistically different across gaps (P<0.05, see Table 2). This difference is more significant when comparing magnitudes of rod torque (***P<0.001, see Table 2). Standard deviation also increases with gap distance, where a Levene's test on torque yielded P<0.05. Pink arrow, branch reaction torque direction.
Grasping wrench space represented by peak foot torque (τrod) and peak leg force (Fnet) at each gap distance. Wrench space is the set of force and torque combinations that permit dynamic landing. Individual points and box plots illustrate data spread with median. Bars represent quartiles. Peak Fnet magnitude and variation tend to increase with gap distance (***P<0.001). Comparisons are made across gaps (***P<0.001). Mean peak τrod is statistically different across gaps (P<0.05, see Table 2). This difference is more significant when comparing magnitudes of rod torque (***P<0.001, see Table 2). Standard deviation also increases with gap distance, where a Levene's test on torque yielded P<0.05. Pink arrow, branch reaction torque direction.
Peak τrod magnitude increased with gap distance (F1,58=31.3, P<0.001) more than doubling from 72 N mm kg−1 at 50 cm to 191 N mm kg−1 at 100 cm (Fig. 5). The range of peak τrod magnitude values also increased with gap distance, tripling from 31 N mm kg−1 at 50 cm to 97 N mm kg−1 at 100 cm (see Table 2; Fig. 5). A Levene's test revealed significant heterogeneity of variances across gap distances (P<0.05). We postulate that the increasing variation in peak force and torque can be explained by the variation in touchdown state.
Forelimb impulse J0 and total impulse JT increased with gap distance (F1,58=81.4, P<0.001 and F1,58=324.7, P<0.001, respectively; see Table 2). However, forelimb impulse and energy contribution, α and β, did not vary significantly across gap conditions (F1,58=0.3, P=0.58, and F1,58=0.5, P=0.48, respectively) indicating that the variation in forelimb contribution is not explained by gap distance. Forelimb impulse contribution (α) was 66.7±10.0% and forelimb energy management (β) was 87.9±7.6%.
Gap distance effect on landing error and angular momentum
The touchdown kinematics can be simplified to two variables: landing error and total angular momentum. Landing error is defined as the vertical distance between the projected ballistic trajectory of the squirrels' CG and the branch center. On average, landing error across all trials was −11.3±16.4 mm and did not increase or decrease with gap distance (P>0.05; Table 2). Normalizing using squirrel body length did not result in statistical differences in landing error.
Total angular momentum is defined as the sum of CG angular momentum, Lcg, and body angular momentum, Lb. Angular momentum of the CG about the branch was significantly different across gap distances (F1,49=30.6, P<0.001), decreasing from 0.073 m2 s−1 at 50 cm to 0.036 m2 s−1 at 100 cm. Similarly, body angular momentum decreased significantly in magnitude with increasing gap distance (F1,49=15.0, P<0.001) from −0.015 m2 s−1 at 50 cm to −0.003 m2 s−1 at 100 cm. Finally, total angular momentum also tended to decrease with increasing gap distance (F1,49=10.5, P<0.01) decreasing from 0.058 m2 s−1 at 50 cm to 0.033 m2 s−1 at 100 cm. Angular momentum was the lowest at the 100 cm gap distance.
DISCUSSION
Reaction force patterns and magnitudes during landing
Boulinguez-Ambroise and colleagues (2023) highlight in their study on tree squirrel jumping that, although jumping performance is often discussed as a pivotal aspect of early primate evolution, its quantification in arboreal mammals is lacking compared with other locomotor behaviors, such as quadrupedal walking and running on branches. The prevalent focus on locomotion along branches (Young, 2023; Young and Chadwell, 2020; Wölfer et al., 2021; Dunham et al., 2019; Hesse et al., 2015; Lammers and Gauntner, 2008; Orkin and Pontzer, 2011; Schmidt, 2011) largely stems from the research in the initial stages of primate evolution, particularly concerning the ability to grasp thin terminal branches. Our study extends beyond assessing movement along branches or saltatorial ability by examining the challenge of executing a stable landing on a narrow branch without the advantage of a prehensile grasp.
Since Bonser's (1999) review, surprisingly few studies involve landing kinetics that consider grasping or balancing. Among the studies that have explored landing dynamics, various animals such as birds, lemurs, cats and toads have been investigated. In our fox squirrels (∼750 g), balanced branch landing most often produced a bimodal branch reaction force pattern (Fig. 4A). Within just 20–65 ms after touchdown, squirrels reached peak reaction forces resulting from front foot touchdown. Hind foot touchdown occurred within 60–140 ms, corresponding to a second, smaller reaction force peak. Peak Fnet increased with gap distance, more than doubling from 2.1 BW (multiples of body weight) at 50 cm to 4.3 BW at the 100 cm distance (Fig. 5; Table 2). Peak branch torque τrod produced by the palmar grasp of the front feet occurred between their front foot touchdown and hind foot touchdown. Average peak τrod also increased with gap distance, more than doubling (72 N mm kg−1 at 50 cm to 191 N mm kg−1 at 100 cm; Fig. 5).
Birds, such as parrotlets, demonstrate an intricate use of their feet during landing on branches, employing opposable digits for effective grasping. While birds primarily utilize their wings to generate supportive aerodynamic forces, Roderick et al. (2019) observed that parrotlets (∼30 g) still exert perch reaction forces ranging from 4 to 5 BW within 5–10 ms. Unlike squirrels, which use both front and hind legs for landing, parrotlets use only their legs to apply ground reaction forces, resulting in a unimodal force pattern. Their feet follow a consistent sequence of movements when landing, including spreading, opening, pre-shaping, wrapping around the branch, and curling their claws. Their study revealed that after touchdown, the dynamics of the foot, toes and claws are crucial for a successful perch, and that these anchoring mechanisms are surface specific. These actions, particularly the toe squeeze, enhance stability upon landing, providing an advantage that squirrels lack.
Landing kinetics have been measured in leaping lemurs (∼1–5 kg) attempting to grasp a compliant vertical pole acting as a force sensor (Demes, et al., 1995, 1999). Peak landing reaction forces ranged from 5 to 11 BW, but no time course was reported. These studies show that some lemurs display a more diverse landing pattern and behavior when leaping down to a flat and horizontal force plate. They show a bimodal reaction force pattern representing the front and hind limbs (Demes et al., 2005). However, the first peak reaction force depended on whether the fore- or hindlimb struck the platform first, which differed among the two species in question. Peak vertical reaction forces ranged from 1.7 to 1.9 BW for the forelimb lander (1.9 kg) and 2.0 to 3.1 BW for the hindlimb lander (3.1 kg), both increasing with jump distance. Hindlimb landing species reached their first force peak at 100 ms and completed landing in 400 ms.
Landing kinetics using extending forelimbs, as shown by our squirrels, have also been measured in cats landing on flat surfaces instrumented with a force plate. Cats (4–5 kg) jumping down from a platform extend their forelimbs for landing impact (Zhang et al., 2014b), which is similar to our observations in squirrels. Cats generate a bimodal reaction force pattern that results from the first peak force from front leg landing followed by the rotation of the body, allowing the hind limbs to produce a second peak force (Wu et al., 2019). The bimodal pattern mirrors what we observed in squirrels, though the relative magnitudes differ. As jump height is increased further, peak reaction forces from the front limb increased from approximately 3 to 8 BW, whereas hind limb forces could increase far more, ranging from 2.5 BW to as high as 20 BW (Zhang et al., 2014a). A faster increase in hindlimb peak force implies that as platform height increases, cats use their hindlimbs more for energy management. Correspondingly, their time to peak force shifts from 10 ms to 40–50 ms. Although not tested in our study, it is likely that at comparably demanding gap distances, squirrels might increasingly rely on their hindlimbs.
Similarly, as hopping toads (250 g) prepare for landing, they also fully extend their front legs or hands (Azizi et al., 2014; Reilly et al., 2015; Cox and Gillis, 2017). Like cats and squirrels, toad landing consists of two phases defined by two peaks of vertical reaction force coinciding with the impact of two body parts. In the hand landing phase, the extended arms hit, flex and rotate as they absorb the landing energy. The vertical force peaks tend to be near 1.75 BW. The body/feet landing phase begins when the folded hindquarters (i.e. pelvis, abdomen and feet) hit the ground simultaneously. The peak vertical force for the latter landing phase was approximately equal to body weight, suggesting that, as in squirrels, the toads' forelimbs manage a significant amount of landing energy.
Force, velocity and leg alignment
Alexander (1991) emphasized the importance of aligning ground reaction force (GRF) vectors with the center of mass and joint centers for energy conservation in legged locomotion. Chen et al. (2006) and Full et al. (1991) demonstrated that animals maintaining a consistent average speed, regardless of their number of legs or posture, tend to align the force vectors along their legs. Vector alignment minimizes joint moments and reduces the work required by the limbs. The significance of alignment can become even more critical during landings at higher speeds, where the primary concern shifts from energy conservation to injury prevention – from efficiency to safety.
For landing squirrels, we used impulse estimates to conclude that the forelimbs and torso manage 80–95% of kinetic energy within 60–140 ms (Table 2). We also measured the angle differences between maximum reaction force angle, velocity angle and leg angle (Fig. 6). We found that the difference between the maximum reaction force angle and leg angle decreased with gap distance (Fig. 6C). The difference between leg angle and velocity angle also decreased with gap distance (Fig. 6D). Both of these differences decrease to zero as the gap distance increases, suggesting that alignment between ground reaction force, touchdown velocity, and virtual leg angle increases as gap distance increases. In other words, the longer the gap distance, the more vector alignment is present, which could result in a reduction of joint moments.
Force and velocity angles and angles relative to the leg at each gap distance. (A) Maximum force angle decreases as gap distance increases (*P<0.05, **P<0.01). (B) Velocity angle increases as gap distance increases. (C) Difference between maximum force and leg angle decreases as the gap distance increased (***P<0.001). (D) Difference between velocity and leg angle also decreases. As gap distance increases, maximum force and velocity align more to the forelimb such that Δθ→0.
Force and velocity angles and angles relative to the leg at each gap distance. (A) Maximum force angle decreases as gap distance increases (*P<0.05, **P<0.01). (B) Velocity angle increases as gap distance increases. (C) Difference between maximum force and leg angle decreases as the gap distance increased (***P<0.001). (D) Difference between velocity and leg angle also decreases. As gap distance increases, maximum force and velocity align more to the forelimb such that Δθ→0.
Like squirrels, cats also align ground reaction forces on jump-down landings with their forelimbs (Zhang et al., 2014b). Specifically, cats land with shallower leg angles to offload GRF to their hindlimbs, thereby reducing the peak force on the forelimbs. The reduction protects the forelimbs from damage in higher speed landings. Wu et al. (2019) have shown that at higher jump heights, cat hindlimbs play a greater role than the forelimbs in absorbing landing energy as the body rotates down and the back bends to allow hindlimb touchdown. In fact, forelimb fractures (38.5%) are less common than hindlimb fractures (61.5%) in falling cats (Zhang et al., 2014b), highlighting the tendency to use the hindlimbs at higher drop heights. Therefore, in cats, posture-dependent actuation prior to touchdown allows the animal to tune the distribution of energy absorption between forelimbs and hindlimbs after touchdown.
EMG data show that cats have a generalized motor program that is independent to drop height and is used to activate extensor muscles at the elbow joint during the pre-landing phase of self-initiated jumps (McKinley and Smith, 1983; McKinley et al., 1983). In addition to pre-touchdown muscle activity and limb coordination, cats also possess passive, post-touchdown landing mechanisms. In particular, they show a remarkable multilevel energy buffering system for shock absorption that includes paw pads, limb bones, and coordinated joints complementing each other (Wu et al., 2020). These results have inspired the design of energy dissipation pads (Lu et al., 2021) and suggestions for legged landing robot design (Xu et al., 2022). Further definition of the complete energy buffering system used in squirrel landing will likely lead to additional inspiration.
In cane toads, Azizi et al. (2014) showed that rapid modulation of hindlimb flexion during the aerial phase of a hop could shift the COM anteriorly and reduce torque by aligning the COM with the GRF vector. A similar study using the same species of cane toads (Bufo marinus) performing controlled landings found that toads use their forelimbs exclusively to decelerate and stabilize the body after impact (Cox and Gillis, 2017). By observing animals jumping from platforms of different heights, they showed that toads achieve dynamic stability across a wide range of landing conditions. Specifically, Cox and Gillis (2017) found that torques during landing could be reduced by a landing preparation motor control strategy for aligning the forelimbs with the body's instantaneous velocity vector at impact (impact angle). As in our squirrels, these two toad studies together show the importance of CG alignment with both velocity and GRF vectors.
Energy absorption can also occur during flight prior to touchdown. In flying squirrels for example, landing force is negatively correlated with glide length (Byrnes et al., 2008). Longer glides allow more time for animals to reach body orientations where they can use aerodynamics to decrease landing velocity, and thus landing forces. In fact, Paskins et al. (2007) suggested flight in flying squirrels may have been selected to control landing forces. A study on birds by Provini et al. (2014) determined that the hindlimbs of zebra finches and diamond doves produce 1.4–2.6 BW forces during landing. It was estimated that for both species, the hindlimbs reduced landing velocity by 60%, thereby contributing substantially to the absorption of kinetic energy after touchdown. The flying robot SNAG (stereotyped nature-inspired aerial grasper) incorporated an independent passive energy absorption for each leg (Roderick et al., 2021). We surmise that the rich morphology of squirrel paws almost certainly contributes to passive energy management upon landing and deserves further attention.
Landing stabilization by control of leg force and foot torque
Spring-loaded inverted pendulum (SLIP) models have long been used to test hypotheses for stable walking and running in animals (Blickhan and Full, 1993; Geyer et al., 2005; Seipel et al., 2004) and a whole host of robots (Raibert et al., 1986; McGeer, 1990). Squirrels, cats, anurans, and birds alike rely on their limbs to manage landing energy post-touchdown, regardless of their ability to control touchdown speed prior to touchdown, and therefore variations of spring-loaded inverted pendulum models could be useful to understand landing.
Zhang et al. (2014a) examined energy absorption and control by spring-mass modeling the front limb behavior of cats jumping down from 1.8 m high platforms onto a force plate. Toad (Cox and Gillis, 2017) and frog (Nauwelaerts and Aerts, 2006) landings after a hop used versions of spring-damper models to determine the alignment of forelimbs at impact and their angle for effective energy absorption. Birds have been studied using spring-mass models to design a flying robot that could perch on branches using an under-actuated, dynamic grasper (Roderick et al., 2021). Using a spring-mass model, the authors defined a ‘perching sufficiency region’. They found that the primary perching failure mode of slipping too far forward or back could be quantified by angular momentum about the branch, which is a function of mass distribution, velocity, and body angles relative to the perch (Roderick et al., 2021). To avoid toppling, balancers can generate force and torque.
Adding leg or radial force control to torque-based balance has been shown to expand disturbance rejection (Caron, 2020; Yang et al., 2021) and improve balance capture regions available to pendulum models (van Hofslot et al., 2019). Using balancing strategies with support forces and torques can also assist by adjusting linear or angular momentum (Orin et al., 2013; Lee and Goswami, 2012). We adopted the extensible pendulum model developed by Yim et al. (2025) for above-branch landing of the monopedal robot Salto. Yim et al. (submitted) proposed two adaptive control strategies that squirrels may be utilizing upon touchdown: leg force control and foot torque control.
Control hypotheses from Yim et al. (submitted) are summarized in Fig. 3B, and they are as follows. All else being equal, they postulate that when an overshoot landing error is present (Fig. 1, lightest red sector), an extensible pendulum can achieve balance by generating a greater leg force and a braking torque (Fig. 3B, overshoot). When an undershoot landing error occurs (Fig. 1, lightest blue sector), an extensible pendulum could achieve balance by producing a lesser leg force and a torque to pull up (Fig. 3B, undershoot). The extensible pendulum model was consistent with the robot's experimental data by predicting the effect of force and torque in the balanceable region. In the present study, we tested the control hypotheses proposed by Yim et al. (2025) by using the model to interpret our results on squirrels.
Testing the landing control hypothesis revealed that all else is not equal for the jumps at our three gap distances. Specifically, landing speed and velocity angle were significantly different at each gap distance, resulting in variation in both landing error and angular momentums (see Table 2). The larger variation in peak torque is likely because longer gap distances required higher take-off and landing speeds, which would yield greater variation observed in touchdown states, and consequently variation in leg force and foot torque. The space of stable above-branch landings is set by landing error and total angular momentum. A sequence of snapshots for a characteristic undershoot and overshoot trial are illustrated in Fig. 3A. The relationship between wrench and landing type indicates that fewer adjustments are necessary for a nominal landing (Fig. 1, white center sector) where the landing error is small and negative while the angular momentum is moderate. In this instance, a squirrel can passively use moderate post-impact angular momentum to compensate for the small, negative landing error, resulting in swinging up towards a balanced posture without the need for significant adjustments.
However, when the magnitude of landing error increases, greater adjustments are necessary. For example, when the landing error is large and positive such that the ballistic trajectory of the CG is above the branch, total angular momentum is also high (Fig. 7A,B in red). Under these overshoot conditions, squirrels apply a large leg force (Fig. 7A) and a braking foot torque on the branch (Fig. 7B). According to the extensible pendulum model, a large leg force is useful for preventing leg shortening or even promoting leg lengthening, which would preserve or increase inertia. Increasing inertia about the branch would decrease angular velocity, resulting in a slower rotation that squirrels can use to avoid swinging over. Likewise, the landing trajectory of the CG could be significantly below the branch (large, negative landing error) and the total angular momentum could be small (Fig. 7A,B in blue). Given these undershoot conditions, squirrels tend to generate lower leg forces (Fig. 7A) and pull-up foot torques (Fig. 7B). According to the extensible pendulum model, a lesser leg force would result in faster leg shortening, which would decrease inertia. Decreasing inertia about the branch would increase angular velocity, resulting in a faster rotation that squirrels can use to swing up on the branch. We have no evidence to determine if squirrels activate muscles to generate additional torque with their feet alone, but grasping digits must provide necessary contact friction. Furthermore, at the gap distance of 100 cm, both peak force and torque seem to be most sensitive to landing error. In other words, for the same landing error but faster landing speed, even more adjustments in force and torque may be necessary for stable above-branch landing.
Squirrel peak force data versus landing error and angular momentum as a function of gap distance. (A) Squirrels tend to apply higher peak forces the more they overshoot (linear mixed-effects model controlling for gap and individual, t1,45=−3.4, P=0.0015). (B) Squirrel peak torque data versus landing error and angular momentum plotted for each gap distance. Squirrels tend to apply higher torque the more they undershoot (t1,45=−3.7, P<0.001). (C) Separate plots from B by gap distance (50 cm, F=4.8; 75 cm, F=6.5; 100 cm, F=8.6; all P<0.05).
Squirrel peak force data versus landing error and angular momentum as a function of gap distance. (A) Squirrels tend to apply higher peak forces the more they overshoot (linear mixed-effects model controlling for gap and individual, t1,45=−3.4, P=0.0015). (B) Squirrel peak torque data versus landing error and angular momentum plotted for each gap distance. Squirrels tend to apply higher torque the more they undershoot (t1,45=−3.7, P<0.001). (C) Separate plots from B by gap distance (50 cm, F=4.8; 75 cm, F=6.5; 100 cm, F=8.6; all P<0.05).
Conclusion
In summary, our measurements on free-ranging squirrels shed light on their remarkable landing dynamics, uncovering findings that govern their agility and stability in arboreal environments. Along with general gap effects on landing kinetics, our investigation revealed three discoveries. First, squirrels exhibit rapid and precise landings, primarily utilizing their forelimbs to manage landing energy. Second, the alignment of velocity and force vectors along the limbs is more pronounced as landing speed increases, reflecting a possible adaptive strategy for managing landing energy. Third, the variability in peak force and torque is consistent with the control of overshooting and undershooting the landing target, suggesting squirrels use radial leg force and foot torque to adjust their landings actively. Specifically, squirrels employ substantial braking forces and torques when overshooting, while utilizing lesser leg force and pull-up torques to correct for undershooting.
Further exploration is warranted to deepen the understanding of squirrel landing dynamics. The extensible pendulum model and the cylindrical approximation of the squirrels' bodies used at touchdown do not include the variable curvature of the spine, the head, limbs and tail. A more complex model might reveal additional mechanisms for balancing moments and be useful for testing hypotheses about the role of the variable body and appendage inertia. XROMM data could better reveal the utility of modeling additional degrees of freedom. One example could be in revealing the stabilizing function of the shoulder and back, as suggested in cats (Zhang et al., 2014a). Kinematic analyses encompassing the entire time series could elucidate the changes in leg length over time, offering valuable insights into how leg force adjustments contribute to leg length changes and the correction of landing error. Additionally, investigating the effects of foothold parameters such as size, curvature and friction on landing control could provide a comprehensive perspective on the adaptability of these findings in different environments.
The implications of the present findings extend beyond the realm of squirrel biomechanics. Advancements in the field of aerial robotics have demonstrated robot abilities such as landing dynamically and perching on cylindrical substrates using specialized grippers (Roderick et al., 2021; Zufferey et al., 2022; Chen et al., 2022; Thomas et al., 2013). These specialized mechanisms ensure anchoring to the substrate and enhance stability during critical landing phases. Specialized gripper designs could be integrated into quadruped robots, which have demonstrated dynamic capabilities such as walking across bricks (Agrawal et al., 2022), walking along a thin walkway (Gonzalez et al., 2020) and even jumping/landing optimally on flat surfaces (Nguyen et al., 2019; Jeon et al., 2022). However, implementation of foot designs more complex than wheels or spheres remains an open challenge owing to the complexity of modeling multiple contacts for legged robot control.
One way we can complement and simplify control is by designing passive feet using compliant structures that react to substrate forces such as bistable mechanisms (Jin et al., 2023), finray designs (Manoonpong et al., 2022; Rozen-Levy et al., 2021) and multi-segment tendon-driven feet (Chatterjee et al., 2023). Our latest work examines the effect that stiffness and damping have on tendon-driven feet for dynamic branch landing (Wang et al., 2024). These passive foot designs simplify legged robot control, and they have the potential for becoming useful for dynamic grasping and detachment for agile robot locomotion on sparse terrains. Exploring the scalability of these control strategies from diverse biological systems holds promise, particularly in the domain of quadruped robot locomotion. The landing mechanisms inherent to squirrels, manifested in their body, limbs and feet, could provide inspiration for the design and control of innovative, agile, legged robots equipped with the ability to rapidly traverse sparse terrains for societal benefit.
Acknowledgements
We would like to acknowledge Lawrence Wang, Lekha Dovvoori, Dominic Pang, Frank Aliaga Auqui, Nicholas Currie, Elaine Llacuna, and Shunyu Wang for their contributions in training squirrels in the field and the data collection process. Portions of the Results and Discussion in this paper are reproduced from the PhD thesis of Sebastian David Lee (University of California, Berkeley, 2024).
Footnotes
Author contributions
Conceptualization: S.D.L., H.S.S., R.S.F., R.J.F.; Data curation: S.D.L., S.W., D.K.; Formal analysis: S.D.L., S.W., D.K., R.J.F.; Funding acquisition: H.S.S., R.S.F., R.J.F.; Investigation: S.D.L., S.W., D.K.; Methodology: S.D.L., S.W., D.K., J.K.Y., N.H.H., R.S.F., R.J.F.; Project administration: S.D.L., R.S.F., R.J.F.; Resources: S.D.L., S.W., D.K., R.J.F.; Software: S.D.L., S.W., N.H.H.; Supervision: H.S.S., R.S.F., R.J.F.; Validation: S.D.L., S.W., D.K., R.S.F., R.J.F.; Visualization: S.D.L., S.W., J.K.Y., R.J.F.; Writing – original draft: S.D.L., R.J.F.; Writing – review & editing: S.D.L., S.W., D.K., E.K.W., J.K.Y., N.H.H., H.S.S., R.S.F., R.J.F.
Funding
Work was funded by the Army Research Office through a Multidisciplinary University Research Initiatives (MURI) grant, ARO #W911NF-1810327 and ARO #W911NF-17-S-002. Open Access funding provided by University of California. Deposited in PMC for immediate release.
Data availability
All relevant data can be found within the article and its supplementary information.
References
Competing interests
The authors declare no competing or financial interests.