Figures and data

The sensorimotor model.
(A) Top view of the cave with three bats. Data are presented for the focal bat (black bat, and wide yellow trajectory). The flight trajectories are shown in thick lines, while the bat’s moment-to-moment decisions are indicated by the inner colored lines (see panel D). Green squares depict reflectors detected by the focal bat along its route. (B) A zoomed-in view of the dashed rectangular area in Panel A, where the focal bat (black) emitted one echolocation call (black filled) and received echoes from the cave walls (green hollow) and from two other bats (blue hollow). It also received conspecifics’ calls (red filled) and their reflection from the cave walls (orange hollow), as well as the reflections from other bats (red hollow). Green squares indicate points that were detected by the focal bat from this call and red x’s indicate missed points due to acoustic masking. The locations of the detected reflectors (green squares) are marked according to their localization by the bat (with simulated errors). The lines near the bats depict their flight direction. (C) The acoustic scene received by the focal bat is as depicted in B, including the emitted call and all received signals (colors as in panel B). (C1) The time-domain plot displays the envelope of signals, encompassing the emitted call and the received signals: the desired echoes from the walls and conspecifics; the calls of other bats; the echoes returning from conspecific calls and reflected off the walls and off other bats. Notably, in this example, some of the desired wall-echoes are jammed by the stronger echoes reflected off conspecifics. (C2) The spectrogram of all the received signals presented in C1; for clarity, the emitted call is not depicted. (C3) The responses of the active channels of the cochlear filter bank after de-chirping. Each channel is represented by its central frequency on the y-axis. Each black dot represents the timing of a reaction that was above the detection threshold in each channel. Note that early reactions in low-frequency channels (marked by yellow arrows) result from the stimulation of those channels caused by the higher frequencies of the downward FM chirp. However, most of these stimulations do not reach the detection threshold and are therefore not detected (see Methods). (C4) The detections of each channel are convoluted with a Gaussian kernel, summed, and compared with the detection threshold (dotted red line). The colored asterisks illustrate detections of the received signals (colors of the sources are as defined above). Panel D depicts the navigation algorithm used by the bat. The algorithm involves a correlated-random flight during the search phase (blue), collision avoidance (yellow), flying along the wall at a constant distance (green), and flying toward the center of a gap between obstacles as an indicator of a possible exit (cyan). After each echolocation call, the bat awaits an IPI (Inter Pulse Interval) period before emitting the next call. Based on the received signals, it then modifies its next call design and adjusts its direction and speed accordingly.

Key model parameters and their effects on performance metrics.
The table presents the key parameters tested, their ranges, default values, and effect sizes on various performance metrics: exit probability, jamming probability, and collision rate. The parameters comprised the number of bats, bat species (PK-Pipistrellus kuhli, RM –Rhinopoma microphyllum), integration window, nominal flight speed, call level, echo mis-identification with temporal aggregation (yes/no), masking (yes/no), and wall target strength. In each scenario, all parameters except the tested one were set to the default value. The effect sizes for each parameter on exit probability, jamming probability, and collision rate are provided (per bat). Square brackets present the minimum and maximum values of the metric across the tested range. Asterisk (*) indicates a significant impact. Each scenario was tested using Generalized Linear Models (GLMs) with number-of-bats and the tested parameters set as fixed explaining variables. Exit probability and jamming probability were treated as binomially distributed, collision rate was treated as a Poisson distributed, and all other variables were considered normally distributed. Explaining variables were set as fixed factors. The number of repetitions for each scenario was as follows: 1 bat: 240; 2 bats: 120, 5 bats: 48; 10 bats: 24; 20 bats: 12; 40 bats: 12; 100 bats: 6. ѱ A significant difference in call intensity was found only for a bat density of 100 bats/3m2, and between the group with a level of 100dB-SPL and all other groups.

Exit performance of P. Kuhli (PK) and R. Microphyllus (RM).
(A) Sensory interference significantly impaired the probability of exiting the cave (compare dashed black line with continuous black line). The probability of a successful exit also declined as the number of bats increased, with no significant difference observed between the species when masking interference was applied. The insert shows the spectrograms of the echolocation calls of PK (top) and RM (bottom). (B) The time-to-exit, which was calculated for successful trials only, and (C) the collision rate with the walls both increased as a function of the number of bats. (D) The probability of jamming significantly increased to about 55% and 63% with 100 bats for PK and RM, respectively. (E) The detection probability of a wall reflector at one meter or less in front of a bat decreased as a function of the number of bats. In panels (A-E), circles represent means and bars represent standard errors (see details in Table 1). Asterisks indicate significant differences between the lines in each panel.

Exit performance as a function of key sensorimotor parameters.
(A) The effect of the integration window on the probability of exiting the roost, the time-to exit, the rate of collisions with the walls, and the probability of jamming (from left to right, respectively). Each colored line shows the trend as a function of the window-size for different bat densities, with each color representing a specific density. (B) The effect of the nominal flight speed of the bats, with panels and line-colors as in panel A. An optimal speed of approximately 6 to 8 m/sec can be observed for all densities above one bat. (C) The effect of call intensity on exit performance, panels as in (A). In all panels, circles represent means and bars represent standard errors. See Table 1 for the number of simulated bats.

The impact of confusion on performance.
The figure illustrates the impact of classification confusion on roost-exit performance under various conditions. Blue lines depict trials with masking, while assuming that bats can distinguish between echoes from their own calls and those of conspecifics (referred to as “No Confusion”). Red lines depict performance where confusion between echoes is assumed. Yellow lines depict performance under the confusion condition, with the added capability of temporal aggregation process in a short-term working memory (referred to as “confusion with mitigation”, see text for further details). In all panels, circles represent means and bars indicate standard errors. (A) The probability of exiting the roost significantly decreased with masking and confusion. In conditions with confusion and no aggregation process, only 15% of bats successfully exited the roost, at a density of 40 bats/3m2. Temporal aggregation partially mitigated the confusion effect but did not eliminate it. (B) Bats with the ability to distinguish between echoes demonstrated significantly shorter exit times than those experiencing confusion. Note that time-to-exit refers only to successful attempts. (C) The collision rate with walls was highest for bats experiencing both masking and confusion but decreased significantly when without confusion. Temporal aggregation process restores performance to the “No Confusion” condition, reducing collision rates accordingly, at densities between 1 to 40 bats/3m2.


Echolocation parameters.
The table presents the echolocation parameters of the two bat species we simulated during the specified flight phases (i.e., search, approach, buzz, and final buzz). In each phase, except for the search phase, in which the parameters remain constant, the parameters for each call are determined by the distance to the closest detected object.

Decision-making in echolocation-based pathfinding
The diagram illustrates the sequential steps in a bat’s pathfinding process based on echolocation. This process starts with the emission of an echolocation call (1) and the reception of echoes and interfering signals (2), followed by sensory processing for detection, range estimation, and direction of arrival (DOA) (3). After integrating detections over a 1–10 call window (4), the bat engages in crash avoidance (5) by evaluating the proximity of conspecifics and obstacles directly ahead. If either is too close, the bat turns to the opposite direction. If no immediate threat is detected, the bat proceeds to pathfinding (6). During pathfinding, it checks for obstacles and, if an opening is detected, flies toward the gap’s center. Without the optional temporal aggregation process (green), the bat simply integrates detections and flies toward the farthest detected obstacle, interpreting it as a wall edge. If the temporal aggregation is included (9), a one-second short memory aids in clustering detections, estimating wall edges, and identifying openings, while also allowing the bat to follow walls at a constant distance. Throughout, the bat continuously adjusts echolocation parameters (8) and controls flight direction and velocity (7) based on ongoing sensory information and decision-making.