Journal of The Royal Society Interface
You have accessResearch articles

Weber’s Law-based perception and the stability of animal groups

Andrea Perna

Andrea Perna

Life Sciences Department, University of Roehampton, London, UK

Unit of Social Ecology, Université Libre de Bruxelles, Bruxelles, Belgium

[email protected]

Google Scholar

Find this author on PubMed

,
Giulio Facchini

Giulio Facchini

Life Sciences Department, University of Roehampton, London, UK

Google Scholar

Find this author on PubMed

and
Jean-Louis Deneubourg

Jean-Louis Deneubourg

Unit of Social Ecology, Université Libre de Bruxelles, Bruxelles, Belgium

Google Scholar

Find this author on PubMed

    Abstract

    Group living animals form aggregations and flocks that remain cohesive in spite of internal movements of individuals. This is possible because individual group members repeatedly adjust their position and motion in response to the position and motion of other group members. Here, we develop a theoretical approach to address the question, what general features—if any—underlie the interaction rules that mediate group stability in animals of all species? We do so by considering how the spatial organization of a group would change in the complete absence of interactions. Without interactions, a group would disperse in a way that can be easily characterized in terms of Fick’s diffusion equations. We can hence address the inverse theoretical problem of finding the individual-level interaction responses that are required to counterbalance diffusion and to preserve group stability. We show that an individual-level response to neighbour densities in the form of Weber’s Law (a ‘universal’ law describing the functioning of the sensory systems of animals of all species) results in an ‘anti-diffusion’ term at the group level. On short timescales, this anti-diffusion restores the initial group configuration in a way that is reminiscent of methods for image deblurring in image processing. We also show that any non-homogeneous, spatial density distribution can be preserved over time if individual movement patterns have the form of a Weber’s Law response. Weber’s Law describes the fundamental functioning of perceptual systems. Our study indicates that it is also a necessary—but not sufficient—feature of collective interactions in stable animal groups.

    1. Introduction

    Gregarious and social animals are capable of forming aggregations and flocks that maintain a relatively stable configuration in spite of individual animals joining or leaving the group and continuously changing their position relative to each other. Examples of such aggregations span the whole animal kingdom from molluscs (e.g. [1]), to insects (e.g. [2]) and crustacea (e.g. [3]) to vertebrates (reviewed in [4]). All these groups can remain cohesive and maintain a stable, coherent organization over time because each member of the group continuously tracks the position or density of its surrounding neighbours and implements appropriate ‘interaction responses’ that effectively preserve group stability.

    A large number of recent studies has aimed at identifying these individual-level interaction responses empirically from tracking data of animals of various species (e.g. [510]). Together, these studies have supported the idea that group formation requires at least some form of inter-individual attraction. However, different studies have pointed to different attraction rules in different species. For instance attraction could be directed to one single individual at a time or simultaneously to multiple neighbours; in turn, these relevant neighbours could be selected based on metric or topological distance, etc.

    The problem is compounded by the fact that the interaction responses of animals of different species ultimately depend on the sensory modalities mediating the interaction (vocal calls, vision, pheromones), and on the underlying neural circuits involved, which are inevitably different from one species to another. In fact, there is large variability across the sensory systems of different animals, with only very few ‘perceptual laws’ that are shared both across taxa and across sensory modalities. One such ‘universal’ perceptual law is the Weber–Fechner Law [11], which states that the ability of a sensory system to discriminate between two physical quantities decreases in inverse proportion to the magnitude of the quantities being compared. To make an example, we can easily tell the difference between a cluster of five objects and a cluster of eight, but we cannot as easily identify the difference between–say—a cluster of 55 and one of 58 [12]. Weber’s Law has been established for humans as well as for a wide range of animals ranging from ants [13], to fish [14], to corvids [15], and to primates [16] and occurs in different sensory modalities. Because Weber’s Law controls the ability of animals to ‘count’ the number of neighbours in different directions, it is likely to play a role also in collective aggregation phenomena.

    Here, we aim at investigating the general properties of the interaction responses mediating stability in animal groups. Whereas many of the previous studies have tackled this problem starting from the analysis of tracking data, here we develop a complementary theoretical approach in which we start by considering the theoretical problem of an animal group whose member individuals do not interact. In the absence of interactions, the dynamics of the group is dominated by random motion and must resemble a diffusion process, which can be accurately characterised in terms of Fick's diffusion equations [17]. We can hence address the inverse theoretical problem of finding which individual-level interaction rules are required to precisely counter balance the effects of diffusion and preserve group stability.

    2. Evolution of densities in the absence of interactions

    In the absence of interactions, the movement in random directions of individual animals results in a change of densities over time that is well described by the diffusion equation [17]. Here, we recall it briefly for the one-dimensional case (the problem is analogous in a higher dimensional space).

    For convenience, we imagine the space to be subdivided in cells of equal width Δx and we indicate with C(x, t) the density of individuals (e.g. number of individuals over volume of the cell) at the position x and time t (figure 1). We focus on density, not on the absolute number of individuals, so that our quantities do not scale with the size of the cells. Within a certain time interval Δt a fraction 2D of individuals present in each cell move to one randomly chosen adjacent cell. This determines the following equation for the evolution of C:

    C(x,t+Δt)=C(x,t)+D[2C(x,t)+C(x+Δx,t)+C(xΔx,t)],2.1
    where the density at some position x depends on the previous density at the same position, minus the fraction of individuals that moved to the neighbour cells, plus the individuals that moved into the cell from the neighbouring ones. In the limit for Δt and Δx small, equation (2.1) becomes
    Ct=D2Cx2,2.2
    where D′ = D((Δx)2t) (see electronic supplementary material for a derivation).
    Figure 1.

    Figure 1. Diffusion and anti-diffusion in a discrete space. (a,b) Diffusion: each particle taking part in the initial distribution illustrated in (a) has a constant probability D of moving to the adjacent cell on the left, and an equal probability of moving to the cell on the right. If—as it is the case in this example—the probability is set to D = 1/6, three particles will move on average from cell x to cell x − Δx and on average one particle will move in the opposite direction from x − Δx to x. While individual particles do not show any directional preference, and they do not interact in any way, the overall process produces (in this example) a net flow of two particles from the region of higher concentration x to the region of lower concentration x − Δx. (a,c) Anti-diffusion based on Weber’s Law: here each particle has a probability of moving against the concentration gradient that is proportional to the Weber’s fraction of concentrations. If the proportionality constant is set to γ = 1/6, two particles on average will move from x − Δx to x. (d) If both diffusion and anti-diffusion take place simultaneously the net effect is that of stabilizing group configuration. (Online version in colour.)

    The solution of the diffusion equation is

    C(x,t)=C(x,0)G(0,t),2.3
    whereby the initial distribution C(x, 0) is convolved with a Gaussian G(0,t) that becomes larger and larger in time. In practice, it is as if the initial distribution was replaced by an increasingly ‘blurred’ version of itself.

    3. Individual-level responses based on Weber’s Law produce collective-level anti-diffusion

    In order to contrast diffusion, the individuals at each given location x must have an increased tendency to move to an adjacent location x + Δx if the concentration of neighbours at the target location is higher than the concentration at x. Weber’s Law states that the probability to discriminate the higher concentration C(x + Δx) from the local concentration C(x) is a monotonically increasing function F of the normalized concentration difference:

    F(|C(x+Δx)C(x)|C(x)),3.1
    In the context of the discrimination between sensory stimuli, the argument of F is often indicated as the Weber fraction of stimulus intensities. Let us consider here how this type of response would affect the evolution of local densities in the same one-dimensional discretized example of figure 1 by focusing on the following model to describe the probability of moving between x and x + Δx.
    p(xx+Δx){γ|C(x+Δx)C(x)|C(x),ifC(x+Δx)>C(x),0otherwise.3.2
    In our specific model, the probability of response is simply linear with the normalized concentration difference (F is simply replaced by the proportionality constant γ) and we do not include explicitly any noise term. We choose this particular formulation because of its nice properties in relation to the diffusion equation, which will become apparent below. In a perfectly biologically realistic scenario, probabilities of response cannot increase linearly (for instance, they cannot become larger than one) and psychophysically measured response thresholds often have a logistic shape, so our approximation is only accurate when the argument of F is small. In electronic supplementary material, section S3 we also consider briefly the case of nonlinear F.

    While we do not include an explicit noise term in equation (3.2), noise is accounted for by the fact that individuals undergo diffusion, and they move up concentration gradients (for example, the probability of actively moving in response to the concentration gradient becomes zero when C(x + Δ − x) = C(x), but individuals can still move in both directions through ‘diffusion’).

    The effect on local concentrations of individuals produced by equation (3.2) can easily be estimated by considering that there are C(x) individuals that apply the gradient–response rule. This will result in a net flow from x to x + Δx equivalent to γ(C(x + Δx) − C(x)) if x + Δx is the cell with higher concentration. If instead the cell with higher concentration is the one at x, particles will flow in the opposite direction, but their flow will turn out to be also proportional to γ(C(x + Δx) − C(x)). As a result the net change of density at x will be:

    C(x,t+Δt)=C(x,t)γ[2C(x,t)+C(x+Δx,t)+C(xΔx,t)],3.3
    whose continuous version is:
    Ct=γ2Cx2,3.4
    where γ′ = γ((Δx)2t).

    The electronic supplementary material provides another intuitive justification for the choice of focusing on Weber’s Law in the case of normally distributed densities.

    Equations (3.3) and (3.4) are identical to the diffusion equations (2.1) and (2.2), except for the sign of the diffusion coefficients D and γ. Importantly, however, their biological justification is completely different: the diffusion equation results from the total absence of interactions, while a response based on Weber’s Law implies an active decision process (see also [18] where a similar ‘anti-diffusion’ term was obtained, although in a slightly different context).

    Care should be taken however that the derivation of the anti-diffusion equation (3.3) from individual-level Weber’s Law responses (equation (3.2)) involves some simplifications. First, equation (3.3) is a mean-field approximation of equation (3.2). Second, if the number of particles or individuals in a cell is close to zero, the probability for these particles to move to an adjacent higher density cell will increase (the denominator of equation (3.2) is small), but the flow can never exceed the number of available particles. There is no control for this in equation (3.3), which means that simply applying this equation can potentially lead to some cells taking negative values.

    3.1. Restoring arbitrary distributions

    Because a response to concentration gradients based on Weber’s Law results in a change in concentrations that is analogous to an anti-diffusion, if such a response is implemented by all individuals in a group, it will have the net effect of partially restoring the distribution that was originally altered by diffusion. To visualize this, consider the example in the top row of figure 2. The example depicts the hypothetical case of animals aggregated at two high-density spots, with densities close to zero elsewhere (figure 2a). Coloured dots mark the positions of randomly selected particles. In the absence of interactions, diffusion operates on densities, which for the purpose of figure 2 we implemented by iterating 300 times a 2D version of equation (2.1) with parameter D = 0.01. The superposed particles also move to a random cell adjacent to their current position with probability D. Figure 2b shows the resulting density distribution and the trajectories of the selected particles. Applying 300 iterations of anti-diffusion response (a 2D version of equation (3.3)) results in a new spatial distribution (figure 2c) where densities are very similar to those in the original distribution.

    Figure 2.

    Figure 2. Weber’s Law based anti-diffusion restores spatial distributions after diffusion. (a) Hypothetical two-dimensional density distribution representing two ‘animal groups’ with a circular profile and homogeneous internal density. Brighter regions indicate higher density. The coloured dots indicate the position of randomly selected individuals. After a given number of iterations of random diffusion-like movements (300 in this example), individuals have moved with random trajectories producing a smoothed density distribution (b). In (c), a density distribution similar to the original (a) distribution is restored by each individual climbing the density gradient with a Weber’s Law response rule (parameter γ = D and same number of iterations: 300). (d,e,f) An identical process to that in A, B and C, but here the starting density distribution (d) is a grey-scale image. Random diffusion of the individuals that form the distribution in D has the effect of blurring the image (producing the image in e) and anti-diffusion restores a similar image to the original (in f). In both examples, the luminance map is updated based on a mean-field approach (i.e. densities are updated in proportion to the calculated probability instead of by actually drawing movement decisions from the probability distribution). (Online version in colour.)

    3.2. An analogy with the problem of deblurring in image processing

    Our problem of animals recovering a group configuration that was previously altered by diffusion is closely related to the problem of deblurring in image processing. In fact, a grey-scale image can be seen as a density distribution (for instance by assimilating bright colours to high density regions and dark colours to low-density regions). Random movements starting from that distribution correspond to blurring the image, and the task of image deblurring algorithms is that of recovering an image similar to the original one. It is well known in the human and computer vision literature that blurred images can be partially restored by subtracting their convolution with a Laplacian operator, in a similar way to what equation (3.3) does. In fact, the diffusion equation (2.2) indicates that the local change in luminance ∂C when we blur the image by a small amount ∂t is proportional to the Laplacian of the original image:

    Ct=D[2Cx2+2Cy2],3.5
    A simple approach for restoring the quality of the initial image is that of playing the diffusion process for negative timescales [19] in a way that is analogous to what is done by equation (3.3).

    Figure 2d,e illustrates the diffusion and Weber’s Law-based anti-diffusion applied to a grey-scale image. The original image (figure 2d) is initially blurred when each of the intensity levels of a pixel randomly diffuses to the neighbouring pixels (figure 2e). Implementing a Weber’s Law-based gradient-climbing response such as the one described in equation (3.2) restores much of the details that were present in the original image. The only difference between this gradient-climbing approach and classical image deblurring by subtraction of a Laplacian convolution is that our approach is based on equation (3.2), while the convolution with a Laplacian implements directly equation (3.3), which can be derived from it. Unlike gradient climbing, subtracting the Laplacian convolution can produce images with negative values; in the context of image processing this is usually resolved by rescaling the image histogram to the allowed range of pixel intensities (e.g. 0–255 for 8-bit images).

    4. Animal movements within stable density landscapes

    The previous examples focused on a somewhat artificial situation whereby particles initially diffuse without interacting and then converge again by implementing a Weber’s Law-based gradient climbing response.

    More realistic is the situation in which a group or a population of animals occupy the environment according to a non-homogeneous density distribution which remains stable over time in spite of individual animals moving in different directions.

    We model this situation by assuming that animals occupy the environment according to a non-homogeneous density landscape, which for simplicity we map to a discrete lattice. Animals are simulated in the model as particles that can move from one cell of the lattice to the adjacent cells, but we impose that the density landscape remains unchanged (see electronic supplementary material for details).

    Suppose that the number of individuals at two adjacent cells i and j is respectively Ci and Cj (figure 3). If the density landscape remains stable, Ci and Cj must remain constant over time (and this for all cells, not just for i and j). The easiest way in which this can be achieved (excluding the contribution of flows ‘running in circle’ over longer loops) is if at any given time the flow ϕ(i, j) from i to j is equal to the flow ϕ(j, i) in the opposite direction.

    Figure 3.

    Figure 3. Flow of individuals across a network of ‘sites’. (a) The simplest way in which a non-flat density landscape can remain stable over time is if the flows of individuals moving across each edge in one direction is balanced by the flow moving in the opposite direction. We built a series of density landscapes such as in (b) whereby each pixel represents a network node and it is directly connected by an edge to its four adjacent pixels. Particles move randomly across the landscape with the constraints that the density distribution is preserved at all times. (c) The probability for a particle to move from node i to node j increases linearly with the difference of density between the target Cj and the source node Ci. The measured slope of this relation, obtained from a linear fit along the dotted lines depicted in (c) is inversely proportional to the local density Ci (d). A spatial distribution of densities (e) can also be summarized in terms of its histogram (f), and a possible alternative definition of group stability is one in which the histogram is conserved over time, while the distribution itself can change. (Online version in colour.)

    Assuming that j is the cell with higher density Cj > Ci, we can say that the individuals that move from j to i do so because of diffusion alone, while the individuals that move from i to j can do it because of diffusion or because of some form of gradient climbing described by an a priori unknown function of local densities f(Ci, Cj).

    We then have

    ϕ(i,j)=CiD+Cif(Ci,Cj)=CjD=ϕ(j,i),4.1
    which gives
    f(Ci,Cj)=D(CjCi)Ci,4.2
    Weber’s Law is hence a solution.

    We tested in simulation if particles moving within an overall stable density landscape appear as if they were following Weber’s Law response to densities. First, we created a number of two-dimensional random density landscapes (obtained as low-pass filtered noise). Then, we generated random particle trajectories with the double constraint that (i) the total amount of time spent by a trajectory at a particular location is proportional to the density at that location, and that (ii) at each time step the net flow of particles between two locations is equal to the flow in the opposite direction (see electronic supplementary material for details; an example of one such density landscape and of some of the corresponding trajectories is illustrated in figure 3b). The probability for particles to move between two adjacent locations calculated over all the trajectories and binned as a function of the local densities at the source (Ci) and at the target location (Cj) is illustrated in figure 3c. While the flow of particles is identical in both directions, individual particles necessarily have a higher probability of moving from low to high-density regions. Intuitively, this is easily understood as a consequence of the fact that fewer particles are available in lower density cells, and as a consequence, they are required to have a higher probability of moving to counterbalance the flow from higher density cells. More precisely, the simulated particles responded to local densities with an apparent Weber’s Law response, whereby the slope of the probability of moving from i to j is inversely proportional to Ci (figure 3d). These simulations show that as long as the overall spatial density distribution is maintained, individual group members will appear to respond to their neighbours with a Weber’s Law type response.

    In our case scenario of a stable density landscape, a gradient response compatible with Weber’s Law results from the requirement of balancing flows for each pair of adjacent cells. This local balancing of flow means that the distribution of concentrations is stable both in space and in time. In other words, if we look at the graph of concentration versus spatial position in figure 3b, of which figure 3e represents a cross-section, this profile never changes. This assumption could be biologically realistic for some groups that remain relatively stationary in time, such as animals forming a lek or nesting in colonies. In many animal groups, such as flocks and shoals, however, the spatial distribution of a group is likely to change in time, while maintaining some conserved statistical properties. As an example, the one-dimensional distribution in figure 3e can be summarized by its histogram (figure 3f) and we can imagine that this histogram is ‘stable’ over time, while the spatial position of zones with different concentration is variable. The stability of the histogram means that the flow of individuals from areas with concentration Ci to areas with slightly larger concentration Cj is the same as the flow in the opposite direction. Given the histogram, there are Ni areas with concentration Ci and Nj areas with concentration Cj, which in a similar way as for equation (4.2) gives the following individual-level probability of moving towards a zone of higher concentration:

    p(CiCj)=D(NjCjNiCi)NiCi,4.3
    which is similar to Weber’s Law as long as NiNj, but would deviate from Weber’s Law when the histogram changes sharply. It would be impossible to generalize here our approach to more complex definitions of stable groups, including to flocks and schools, because this would require coming up with relevant definitions of ‘local density’ for a group that is translating in space, and also addressing the issue that many existing models of these systems predict that phase transitions would occur between different flocking modes (e.g. [20]). In other words, the diffusion coefficient might be itself density-dependent in these groups.

    5. Long-term dynamics of iterated diffusion and anti-diffusion

    While Weber’s Law-based responses to densities can counteract diffusive forces and are likely to be consistently observed in stable groups with an internal dynamics, can density responses based on Weber’s Law alone support group stability? Because of the linearity of the diffusion and anti-diffusion equations (equations (2.2) and (3.4)), a combination of diffusion and anti-diffusion has only one stationary state with homogeneous density. The homogeneous density state is stable—and the group disperses across the environment—when diffusion is stronger than anti-diffusion (when D > γ). When instead anti-diffusion prevails (γ > D), the homogeneous state becomes unstable and the model predicts the formation of explosively larger groups. In particular, the dynamics of the aggregation is such that the higher spatial frequencies (small spatial scales/small group sizes) are amplified or attenuated faster than lower frequencies (see electronic supplementary material for an explanation). An important consequence of this is that as soon as the regulation mechanism is noisy, the high-frequency noise is quickly amplified.

    It is important to remember that in making the considerations above we were relying on two assumptions. The first assumption is that diffusion and anti-diffusion take place simultaneously and the gradient climbing response is perfectly accurate. In reality, if the anti-diffusion response is not implemented immediately the anti-diffusion will take place on a distribution that has already been altered by some uncompensated diffusion. Similarly, if the gradient climbing response is noisy, i.e. it also diffuses a bit, this will also result in some uncompensated diffusion. The second assumption is that local densities can grow with no upper limit. In real-world aggregation phenomena, however, densities will reach a saturation at some point. Here, we incorporate these additional elements in a discrete model of diffusion and anti-diffusion. The model uncouples in time diffusion and anti-diffusion, and implements a saturating anti-diffusion response to prevent densities from increasing above a saturation point. We ran computer simulations in which an initial density distribution is altered by alternating steps of diffusion and anti-diffusion in discrete time steps. The first step is a simple diffusion step:

    C(x,y,t+1)=C(x,y,t)+D(C(x,y,t)),5.1
    where
    (C(x,y,t))=4C(x,t)+C(x+1,y,t)+C(x1,y,t)+C(x,y+1,t)+C(x,y1,t).5.2
    The anti-diffusion step is similar to all previous examples, but we also include an additional term (C(x, y, t + 1)[1 − C(x, y, t + 1)]) to keep densities within the range from zero to one:
    C(x,y,t+2)=C(x,y,t+1)γ(C(x,y,t+1))C(x,y,t+1)×[1C(x,y,t+1)],5.3
    In the examples shown in figure 4, the initial condition is a uniform random distribution and the distribution shown corresponds to 15 000 simulation steps.
    Figure 4.

    Figure 4. Pattern formation from iterated diffusion and anti-diffusion. (a,c) Spatial patterns produced on a 128×128 lattice after 15 000 iterations of the long-term dynamics model; the starting distribution is random white noise. In both figures γ = 0.9. In (a) D = 0.13, while in (b) D = 0.17. After multiple iterations of the simulation steps, a pattern emerges at a single spatial scale. This is clearly visible in (b,d), which plot the amplitude spectrum of the Fourier transform of the patterns in (a) and (c), respectively. (Online version in colour.)

    Under these conditions, we do observe the appearance of spatial patterns of a characteristic scale, at least for some parameters (figure 4). While our particular implementation of this long-term model may lack biological realism, for instance in relation to the temporal separation between diffusion and anti-diffusion periods, alternative models that also incorporate a reaction time in the response of individuals are likely to predict the formation of patterns with a characteristic scale: it is implicit in the nature of the diffusion equation that the introduction of events with a characteristic temporal scale—for instance, a fixed reaction time in the response to neighbours, also leads to the appearance of spatial patterns with a characteristic spatial scale.

    6. Discussion

    Gregarious animals can form non-homogeneous density distributions and aggregations that persist over timescales typically much longer than the timescales of movement of individual animals. In our study, we show that group stability constrains the individual-level interaction rules that mediate group cohesion, which are required to follow Weber’s Law.

    Weber’s Law describes a general functioning principle of sensory systems. Traditionally, Weber’s Law-based sensory perception has been discussed in relation to the mechanisms and to the constraints of information processing in the brain [21,22]. Here, we have shown that Weber’s Law is also an essential property of collective interactions in stable animal groups.

    Sensory perception plays a fundamental role in guiding the interactions of an animal with its environment [23], so it can reasonably be expected that the fundamental properties of sensory systems also play a role in shaping inter-individual interactions and collective behaviour. Previous studies have already established the relevance of decision-making mechanisms based on Weber’s Law for explaining collective decision-making in ants, bees and fish [13,2426]. Theoretical studies also indicate that Weber’s Law response results from optimal pooling of group-level information for collective decision-making [26].

    In the present study, we do not make any assumptions about the mechanisms or about the optimality of the decision-making rules in individual animals. Our only assumption is that animal groups remain stable over long timescales, compared to the timescale of individual animal movements. We further assumed that in the absence of interactions group dynamics would be described by diffusion. Because of the generality of these assumptions, we can conclude that Weber’s Law response to animal densities is a general property of all stable groups.

    We should, however, be careful not to conclude that all animal aggregation phenomena depend on Weber’s Law-based responses to neighbour densities. Many animal aggregations can be produced also in the complete absence of social interactions and of gregarious behaviour. This is the case for instance when animals are attracted to an environmental feature such as a water source or a particular type of vegetation. In this case, group stability would be predominantly mediated by the external attractor. Yet, our analysis shows that if we—correctly or erroneously—consider that the attraction responses are directed to other group members and not to the environmental feature, whenever the group density is stable, we would observe a response to neighbour densities in the form of Weber’s Law. In some species, the rules of interaction could potentially be extremely complicated, and involve both inter-individual attraction and repulsion. However, as long as the overall effect of these interactions is that of keeping the group stable, we would observe a Weber’s Law type of response to neighbour densities.

    Our model is not a morphogenetic model: the type of response to neighbour densities based on Weber’s Law that we describe here cannot be used to predict the size or the shape of the groups formed by animals of a particular species. We argue that this is an inevitable consequence of the ‘universality’ of Weber’s Law: animal groups present a huge variation of size and shape across species which is unlikely to be accounted for by a single shared perceptual rule. In addition, the size of animal groups is typically widely distributed also within one single species, suggesting that its regulation depends on other factors, such as merge and split phenomena, rather than being determined by sensory responses alone [2729].

    While we do actually show that an aggregation model based on Weber’s Law can lead to the formation of patterns at a fixed scale, this scale is mainly determined by the delay that we introduced between diffusion and anti-diffusion. We could try to find a resemblance between the patterns produced in our long-term iteration model and the spatial distribution of individuals in certain animal groups, such as for instance mussel beds [1]. However, we think that it is safer to argue that these pattern formation phenomena leading to a small characteristic spatial scale have little biological relevance in general, for instance, because at such small spatial scales the positioning of individuals is more affected by direct individual to individual interactions than by the generic Weber’s Law response that we consider here.

    Paradoxically, while we predict that individual-level interaction responses based on Weber’s Law should be observed in stable animal groups, Weber’s Law alone cannot explain the formation of groups of a particular form or size and is not sufficient by itself to explain the stability of pre-existing groups (because of reaction times and amplification of noise). We can conclude that Weber’s Law based interactions are a necessary but not entirely sufficient feature of stable group.

    Our study points to a general relation between Weber’s Law and collective behaviour. While Weber’s Law can explain several perceptual phenomena, there are many instances in which sensory perception deviates from Weber’s Law behaviour. Future studies should try to relate also deviations from Weber’s Law, as well as other perceptual phenomena, to collective behaviour. For example, many sensory stimuli remain undetected when their intensity falls below a perception threshold. In the context of social interactions, perceptual thresholds for responding to conspecifics can directly affect the way animals rely on private versus social information for taking decisions.

    Our approach for deriving interaction responses in animal groups is complementary to other studies based on direct observations of interacting individuals. While direct observations of behaviour inform us about how animals of a particular species interact, our approach informs us about which individual-level interactions are ‘unavoidable’ given a particular collective phenomenon. In this particular case, individual-level responses in the form of Weber’s Law are unavoidable in stable density groups. This ‘inverse problem’ approach has been previously applied to the study of other collective behaviour phenomena (e.g. [30,31]) and we are confident that in future our approach can be extended to the study of also other more complex forms of group coordination.

    Data accessibility

    This article has no data. The code used in the simulations is available at https://github.com/pernafrost/weber_law_antidiffusion.

    Authors' contributions

    Formulated the initial research idea and approach: A.P.; discussed various modelling approaches: A.P., J.-L.D.; developed the non-linear Weber’s Law models: G.F.; implemented the computer-based simulations: A.P.; wrote the article: A.P. All authors were involved in discussions on different aspects of the study.

    Competing interests

    We declare we have no competing interests.

    Funding

    A.P. was supported during part of this study by a fellowship from the Belgian National Fund for Scientific Research (FRS-FNRS). G.F. is supported by a Newton International Fellowship from the Royal Society.

    Acknowledgements

    J.-L.D. is a Research Director from the Belgian National Fund for Scientific Research (FRS-FNRS).

    Footnotes

    Electronic supplementary material is available online at http://dx.doi.org/10.6084/m9.figshare.c.4487192.

    Published by the Royal Society. All rights reserved.

    References