Simulating Populations VII: the Correlogram View of Complex Dispersion

In my continued quest for a more realistic statistical-mechanical theory for spatially extended population dynamics I have previously pointed out a specific property of the inter-scale spatial coefficient of variation as one of the hallmarks of scale-free dispersion (see Part III). In the present post I study another statistical property, the spatial autocorrelation, which may provide additional cue about the population’s compliance with standard or complex space use.

First, consider the standard theory, based on mean field compliant population redistribution (mixing). The following three images show a typical example, where the population is subject to 5% diffusion rate at unit (pixel) scale, net population growth of 1%, no Allée effect, and over-all population density below carrying capacity. As previously described in this series, diffusion tends to smoothen the density surface. The log(M,V) plot over a scale range typically shows y-intercept [log(a)] substantially below zero  a hallmark of fine-scale smoothness – under condition where the slope β ≈ 2. I refer to previous posts in this series for technical details.


The correlogram (above) reflects the undulating density surface.

Next, consider the following scenario under the same condition, except for 5% standard diffusion at unit scale (k=1) being replaced by 5% scale-free “zooming” with 1% pr. scale level.


The present  snapshot of population dispersion represents the population a few time increments after a general population crash (bottleneck episode), when the population is in the process of re-organizing itself. Despite the short time span since the episode (5-6 time increments) the log(M,V) plot has already adjusted itself from log(a) >>0 and b ≈ 1 (random mixing due to external forcing) to log(a)≈ 0 and b ≈ 2.

Observe that in this additional example of complex population dynamics the log(M,V) plot again satisfies log(a) ≈ 0 when b ≈ 2. Here the local effect at next time increment from overshooting local carrying capacity both locally and elsewhere is influenced by a rate of 50% remaining population at the actual location (as in previous examples) to zero individuals remaining. In other words, this hallmark of Zoomer-like dynamics is quite resilient to this modification of ecological conditions.

The interesting aspect in the present context is its spatial transect correlogram at unit scale k=1. It shows low level of autocorrelation at all spatial lags except for scale 0, which trivially illustrates that the local population correlates with itself 100% at lag zero. Despite the non-significant autocorrelation the parameter condition for log(M,V), log(a) = 0 and b = 2. 

Thus, scale-free population abundance regards both spatially autocorrelated transects (as shown in previous parts of this series) and – as shown here – non-autocorrelated transects a short period following a perturbation [during the first 3-5 increments after the event*, log(a) >> 0 and b ≈ 1].

In this respect I refer to the so-called Z-paradox, which is resolved under the Zoomer model but not under the standard framework of population modelling. Thus, my proposed model may also provide a novel approach towards the famous and controversial Taylor’s Power law.


*) Under the condition that the bottleneck condition lasts for one time increment. If the change of condition is permanent, it may typically take 20-50 time steps to restore log(a) ≈ 0 and b ≈ 2.

Simulating Populations VI: the Unrealism of Standard Models

In the previous post (Part V) I suggested that the standard theoretical framework for simulating population dynamics/kinetics is fundamentally unrealistic, since in any open environment it dooms local population abundance to approaching zero and extinction due to memory-less dispersal and further lubricated by Allée effects. The only apparent rescue conditions for a population are (a) the population lives in a closed rather than open environment, which is an unrealistic scenario in general terms; or (b) on average the net growth rate is larger than the dispersal (diffusion) rate. In this post I argue that even the latter assumption is flawed.

The traditional Coupled map lattice (CML) modelling, partial differential modelling, and other standard spin-offs from classical statistical mechanics may apparently be “rescued” with respect to model realism by defining unit spatial scale coarse enough to ensure that the net growth rate is stronger than diffusion rate at this level. Since diffusion is a scale-specific process – see below – and growth (specific birth and death rates) are not linked to a specific scale, the diffusion rate will apparently be reduced when observed from coarser spatial resolutions. For example, a ten times larger spatial resolution (unit pixel size) for a given population study on average embeds ten times more individuals. Hence, the percent of the population that disperses to the coarser-scaled neighbourhood becomes smaller from this perspective due to smaller perimeter/area, and may even be smaller than the population’s net growth rate. Paradox solved? Unfortunately not, because the reasoning is flawed.

Together we thrive. Photo: AOG.

First, consider the condition where the individuals are living in an open environment, which is the general condition. No species is abundant everywhere, meaning that a given population is surrounded by unoccupied space*).

The correct scale to study the diffusion rate is not an arbitrary scale (pixel size), defined as unit scale, but the actual population’s characteristic diffusion scale under the given environmental conditions. In other words, diffusion is an intrinsic property of the population, and not something that can be defined by the observer at will.

To understand this concept of characteristic diffusion rate we need to turn to statistical mechanics. The characteristic diffusion scale is given by the diffusion coefficient, which is proportional with the square of average step length by the average individual divided by the time interval for non-correlated successive directional change for the individuals’ random walk process (Brownian motion-like when we consider the statistical property of a memory-less kind of movement, which is a premise for the standard modelling framework).

The diffusion coefficient-determined rate; i.e., the square of the net length moved pr. unit time, as determined by the individual random walk properties**), is the parameter to be compared to the net growth rate of the population. Even if we consider that a diffusion in general terms is a slow process, the diffusion’s characteristic rate at its unit time scale becomes substantially larger when we re-scale the time axis to the characteristic reproduction interval for a new individual to be born (generation time divided by number of individuals in the brood, and adjusting for death rate effects). Under this biologically adjusted diffusion rate by temporal rescaling to the system’s population dynamical scale, it may easily become very much larger than the net reproduction rate. In short, such temporal rescaling ensures full “local mixing” of individuals. However, in an open environment such high rate of individual reshuffling will doom the population to extinction.

In an open environment the population is typically surrounded by suitable habitat, but for the moment it is not utilized by the species. This is in accordance to common empirical experience, not from reading Ecology text books… Thus, the population’s fringe zone “leaks” individuals to the surrounding area, due to non-zero outward dispersal in net terms. Outward-directed flow of individuals is larger than inward-directed flow, since the latter has no or negligible population source. Based on the argument above that the dispersal rate generally is substantially larger than the population’s net growth rate when diffusion is re-scaled to the characteristic time scale of population growth (see above), the population is doomed under standard model premises. Another way to look at it is to consider diffusion  rate at its correct scale given by its diffusion coefficient, and only consider the (fuzzy) zone along the population perimeter where there is a net flow of individuals. As illustrated in my previous post (Part V), unless you put up a fence the local abundance drifts towards zero!

Consider dropping a single drop of ink into a water-filled bathtub – representing our “open environment”. The black dot will over time spread out and become uniformly dispersed. Consider that the ink molecules have a limited life span, and (hypothetically!) two molecules need to meet within their respective life span to produce additional ink molecules. Since we are considering a bathtub; a constrained space, the ink “population” may still survive. However, the larger the bathtub the smaller the rate of inter-molecular encounters. The critical space size is given by the level (bathtub size) where the given diffusion rate – after temporal rescaling to match the net reproductive rate (see above) – becomes larger than the reproduction rate pr. individual of new ink molecules! This balance between diffusion rate and reproduction rate is diminished by increasing bathtub size, invoking an Allée effect in larger bathtubs!

What about a metapopulation system, where migration between sub-populations (a representative for diffusion at the scale of sub-populations) is relatively small? Such a “semi-permeable” kind of population sectioning puts constraint on the much higher diffusion rate inside sub-populations (local high-speed mixing as seen from the perspective of the net reproductive rate). In fact, in standard models for spatially extended populations it is assumed that the characteristic scale for diffusion ensures a high rate of population mixing at relatively fine spatial scales, relative to the extent of the population under study. For example, it ensures a “smooth” and thus differentiable density surface of the population’s spatial dispersion; a premise for realistic application of partial differential equations. However, even if a metapopulation system with mixing constraint on intermediate spatial scales may avoid the paradox of population extinction as outlined above, such a metapopulation system is doomed since it exists in an open environment at its fringes! Huffaker’s 1958 experiments on sub-populations of mite were run in a constrained environment, the size of the terrarium!

To conclude, in my view ecological theory for population dynamics needs a more realistic kind of modelling platform. At the individual level spatial memory and scale-free space use is now being empirically verified with a high pace. This insight needs to be reflected in theory of population level space use. As stated in my note below I propose that intraspecific cohesion, for example as implemented by the individual-level conspecific attraction property of the Zoomer model, may resolve the “doomed to extinction” paradox of standard population models for populations in an open environment. In this framework, the standard diffusion is replaced by the scale-free redistribution process termed zooming. However, this process is not compatible with the standard statistical-mechanical theory for population modelling.


*) In some cases this fringe zone may be easily understood from the perspective of unfriendly neighbourhood in habitat terms, but often the surroundings appear quite similar to the conditions inside the distribution range. For example, the over-all population may be spatially fragmented with respect to abundance; particularly along the core areas’ perimeters, with no apparent reason why small and large chunks of the intermediate areas should not be habitable. Pick any well-studied species, and ask an expert on its distributional range to explain population absence in some regions with apparently similar ecological conditions to the population’s present range. In an in my view unsatisfactory manner the way out of this dilemma (the “shoot from the hip” answer) is normally to point to some environmental factor still not revealed, or some kind of statistical chance effect. In my Zoomer model such apparently inexplicable “clumping” of a population is explained as an emergent property of conspecific attraction and scale free redistribution of some individuals. A given number of individuals cannot be everywhere all the time. However, in the present post I’m concerned about the basic premises of the standard framework for population dynamical modelling, not how the empirical paradox may be resolved.

**) As repeatedly underscored in my book and in this blog, individuals generally do not move in a stochastic manner. However, the animal’s path may be represented by a statistical function in over-all terms. In the standard framework this function is assumed to satisfy the parameterized random walk model of the Brownian motion type.

Simulating Populations V: Bottlenecks and Recovery

Time to simulate a stress-test of the two population-kinetic frameworks, the Coupled map lattice model and the Zoomer model! Consider a scenario where some kind of environmental event has crushed the population to about 1% of its normal carrying capacity. In addition, the remaining population has also become spatially fragmented during this catastrophe. Then consider that the condition improves to the pre-event level. What is the population’s potential to recover under the two scenaria you have become familiar with in Parts I-IV, scale-specific and scale-free kinetics?

The map to the right shows the small population’s spatial dispersion at the start of the potential recovery phase. Isopleths indicate local population density, which shows an average of 165 individuals pr. occupied cell at unit scale while the carrying capacity (CC) has been restored to a potential for 5,000 individuals at this scale. In other words, most local populations have gone extinct as a consequence of the recent crunch event.

Then the recovery phase begins to run. Starting with the standard condition of scale-specific population dynamics/kinetics (Coupled map lattice model) and setting diffusion rate at unit scale to 5% and net population growth of 2%, the following image shows the population dispersion after 20 iterations.

Since the diffusion rate is larger than the local growth rate (the general condition of spatially unconstrained animal populations) and the population is now surrounded by unoccupied area, the population is drifting towards extinction!

This faith is also facilitated by an additional model condition, Allée effects. At this low level of population abundance it is important to consider and implement three aspects: First, accelerated extinction at very low abundance levels have to be introduced. Here I set the critical level to 50 individuals pr. unit cell*). Below this level, the population is reduced by 10% pr. time increment. Second, due to the low abundance levels, one has to consider that individuals exist as discrete entities, not fractions of numbers (at high abundance the difference between discrete and continuous numbers are insignificant). Third, at very low population densities random events take its toll. I implement this as some noise level on the survival rate in the Allée zone; i.e., in cells with less than 50 individuals.

In contrast, does scale-free and memory-influenced zooming influence the population’s otherwise dire faith after the catastrophic event? Obviously it does. The Zoomer snapshots below at t=20, t= 200, t= 500 and t= 1,000 shows a population in healthy recovery, despite  being surrounded by a wide zone of unoccupied space. The 5% diffusion rate under the CML condition above is replaced by a 5% zooming rate, with 1% redistribution pr. scale level (see previous parts of this series). In other respects the conditions are similar to the CML model, including net growth being smaller than individual reshuffling at unit scale.


The log(M,V) result at t = 1,000 (as in the earlier parts of recovery; not shown) shows full compliance with intercept ≈ 0 and slope ≈ 2, as predicted by the Zoomer model.

In this manner the Zoomer model illustrates – and potentially resolves – some crucial but under-communicated issues with respect to the standard modelling framework.

Inclusion of spatial memory and strategic space use – in particular the capacity for individuals to include conspecifics as part of their resource map at strategic scales – counteracts the otherwise detrimental effect of living in an open world.  At the fringe of any animal population, under the standard modelling paradigm local abundance is constantly threatened by individuals getting lost in space, by drifting away from sufficiently strong contact with conspecifics (ref: diffusion and Allée effects). The Zoomer design, by implementing spatio-temporal memory, formulates a solution to this core problem for population dynamical modelling. However, the solution requires an extended kind of statistical mechanics. Read my book – for the time being the main source (and for some parts the only source) for a theoretical overview of this approach!

In the next post I will address an expected primary objection to my quite far-fetching conclusions above, that traditional population dynamical modelling is based on shaky assumptions with respect to realism. In future posts I will also present empirical support for the Zoomer model.


*) There is nothing magic about N=50, but to avoid a more complicated formula for the Allée effect – with little or no advantage with respect to model realism in over-all terms – I have just chosen a “small abundance number” relative to the carrying capacity.

Simulating Populations IV: Environmental Heterogeneity

In the foregoing Parts I-III model complexity was increased in a stepwise manner for the sake of exploring intrinsic population behaviour one factor by the time. In this post I take one additional step by studying the overall effect of environmental heterogeneity.

For the time being I consider spatial heterogeneity only, leaving temporal fluctuations to a later post. Anyway, a new level of realism is hereby added relative to the scenaria in the previous posts: external influence is now adding to intrinsic processes with respect to variations in local population abundance.

Consider the Zoomer snapshot to the right (zero diffusion, 5% zooming over a scale range, as in previous examples), after the population has progressed 100 time steps in an environment where the local carrying capacity varied over space (CC=2942 individuals pr. cell at unit scale within the arena, on average). Due to 50% overshoot survival (see
Part I for a definition of CC), the population fluctuates between CC and CC/2.

The net growth rate is small at the defined time resolution (1%), and a given bust event at local density is passing local CC brings the local population down to CC/2 in a single time step. While it takes only one iteration to bring the local population down to CC/2, it takes many increments to bring it up towards CC again.

The standard, non-scaling Coupled map lattice condition given the same heterogeneous map for local CC variability is shown by the two images below.

The important pattern in the two sets of Figures above and below is the apparent similarity of the log(M,V) plot in the present condition of habitat heterogeneity and the condition of habitat homogeneity that was presented in Parts I-III. This similarity makes sense, since local variability as a consequence of local habitat heterogeneity needs to be analyzed at  a finer scale than the entire arena size and then compared between sections. Under all scenaria so far, the log(M,V) plots regard population abundance within the arena as a whole.

Thus, in both scenaria above the local variability is hidden; i.e., “averaged out”.

Anyway, there is a crucial difference to observe between scale-specific (CML compliant) and scale-free statistics. The intercept log(a) << 0 while the slope b ≈ 2 under the standard CML condition, and log(a) ≈ 0 when b ≈ 2 under the Zoomer condition. Thus, we can conclude that the novel indicator of complex population dynamics – self-similar population dispersion due to CV≈1 (see Part III) apparently stands the heterogeneity test!

Such resilience to environmental conditions when is comes to distinguishing standard from complex space use is of course crucial for the realism of this system property when we later on are confronting the theory with real data.

After these step-wise system introductions for the sake of revealing the respective systems’ intrinsic population kinetics, time has come to throw additional realism into the model conditions. As a starter with respect to cruising towards ecological aspects I study the populations’ response to bottleneck events (population crunches) and their ability to recover under standard and alternative statistical-mechanical premises! Look forward to Part V.

Simulating Populations III: a New Statistical Indicator of Complex Population Kinetics

In the previous parts I-II of this series I described two variants of spatially extended population dynamics, represented by a standard Coupled map lattice (CML) model and the Zoomer model. In this post I show how a specific statistical-mechanical property of scale-free space use may reveal the difference between these two space use conditions despite an apparent similar level of spatial autocorrelation below the population’s carrying capacity.

First, a brief summary of the model conditions (for details, see Part I-II):

  • The environment is set to be homogeneous (soon to be relaxed), in order to have focus on intrinsic population kinetics.
  • The time resolution is set to be fine-grained, implying that the main driving force for change during respective time increments is individual re-shuffling rather than birth and death rates (net growth rate set to ca 1%).
  • For the CML examples (implying a scale-specific process in statistical-mechanical terms), the mean field assumption applies. Reshuffling happens by a small diffusion rate 0-1%, simplified as standard random walk at the individual level, taking place between neighbour cells (1%) at unit scale or totally hidden at finer (unobserved) scales below unit scale (0%).
  • For the Zoomer examples (implying a scale-free process), reshuffling happens by conspecific attraction over the scale range from unit scale to the extent of the simulation arena. When the environment is set to be homogeneous, conspecific attraction rate of 1% pr. scale level and 5% in total is the only remaining driving force (beyond dispersal due to local overshoot of carrying capacity) for population redistribution beyond unit scale. However, diffusion (rate 0-1%) may be present at unit scale to account for presence or absence of fine-grained stochastic noise (short term tactical response from individuals to local events, independent of the strategic zooming behaviour).

Why not including conspecific attraction under the standard CML scenario? Because one of the assumptions for mean field population kinetics, represented by a CML model, is absence of spatial memory effects on redistribution beyond the unit scale. Hence, heighbourhood re-shuffling of individuals can be represented by a standard diffusion compliant process only (assuming that advection – another scale-specific process – is also limited to finer scales than the unit grid). Like it or not! I’ll return to this topic in a later post.

While non-autocorrelated population dispersion (exponential distribution of local population density, V ∝ M)  is the typical pattern of a standard population scenario if a substantial part of the population is fluctuating close to the carrying capacity (see Part I, example 1), a similar result is difficult to achieve under scale-free space use and conspecific attraction.

Spatial autocorrelation (V ∝ M2) is swiftly restored – typically within a few time increments – following instances where a larger than normal part of the local populations are crashing simultaneously. The example above shows a snapshot close to the carrying capacity*) under condition of local overshoot crash rate of 50% (diffusion rate = 0%). A similar result appears even under a more dramatic local crash rate of 90% (not shown), which may be compared to the CML Example 1 in Part I.

Obviously, since standard and scale-free statistical mechanics for spatially extended populations are so fundamentally different i qualitative terms (for example, autocorrelation appears from different processes; growth/diffusion and intrinsically driven conspecific attraction, respectively), it should be a key goal to pinpoint statistical methods to distinguish between these two conditions for population dynamics. Here I show for the first time how spatial autocorrelation of the population’s density surface, here to be studied from the log(M,V) pattern, may provide a simple method in this respect.

Recall from Part I-II that M regards the average number of individuals pr. grid cell at a given scale, and M is changing proportionally with grid resolution (times larger linear scale gives k2 larger grid cells, which give k2 times larger M). In absence of spatial auto-correlation, the variance of M at respective scales changes proportionally with M.

To pinpoint one of the hallmarks of scale-free population dynamics, consider these log(M,V) snapshots, where I have interpolated the regression lines towards the y-intercept [giving the parameter log(a) in log(V) = log(a) + b*logM].

First, for the sake of comparison, two examples from scale-specific dynamics:


Then consider the scale-free kind of population kinetics:


In the two first of these zoomer examples, diffusion rate = 0%, while in the third example diffusion is set to 1%. All other conditions are similar to the scale-specific scenaria above.

The key difference regards the y-intercept, log(a), which – under a premise of little influence from local carrying capacity at the time of the actual snapshot – is strongly negative for a scale-specific scenario and close to zero for the zoomer conditions. The magnitude of the negative intercept from scale-specific dynamics depends on the average M; more abundant populations (from the perspective of the chosen scale range) will show a “parallel shift” towards the right-hand side in the log(M,V) plot, and thus show an even more negative log(a).

To understand this difference between scale-specific and scale-free population kinetics, observe that the power exponent b ≈ 2 implies that √V; i.e., the standard deviation (SD), is proportional with M. Further, if log(a) = 0 (≈ 1) in addition to b ≈ 2 the coefficient of variation (CV) is approximately constant over the range of M (recall that number of individuals, M,  is proportional with grid cell size, which is proportional with square root of scale, k);

CV = SD/M ≈ 1     | M ∝ √k

Only the scale-free zoomer model condition produces results under condition of spatial autocorrelation in compliance with a constant coefficient of variation over the given scale range. The population is spatially self-similar; i.e., it satisfies a statistical fractal. In contrast, the classic CML model produces a “smooth” density surface (0<CV<<1) at fine resolutions**).

Since the y-intercept under the Zoomer model scenaria approximates log(a) = 0, the CV≈1 condition is not sensitive to population abundance. In my book I describe other methods to study a population’s scaling properties under a premise of Multi-scaled random walk at the individual level, for example in the context of Taylor’s power law. I also provide empirical support for scale-free population dispersion.

In future blog posts I will bring the analysis further, by studying the influence of spatial heterogeneity and respective models’ responses to perturbation events.

I kindly ask you to give credit to my blog and my book when you refer to the statistical-mechanical theory for scale-free population dynamics! To my knowledge, for the moment these are the only sources for this departure from standard population dynamical modelling.


*) As mentioned in Part I (see its Note), the carrying capacity is given a special interpretation as the threshold whereby the local individuals (within the finest defined grid scale) are moving out, redistributing themselves or dying.

**) Recall from previous posts (and my book) that the standard framework to model spatially extended population dynamics – Coupled map lattice models and partial differential equations – in fact depend on such a smooth density surface at fine spatial resolutions; i.e., in the (hypothetical) theoretical limit of infinite local population density. Otherwise the population is not differentiable! If this premise fails, the standard approaches are doomed to fail with respect to the model’s realism and predictive power.

Simulating Populations II: Adding Spatial Memory and Scaling

In Part I of this series I presented a couple of spatially extended scenaria of intrinsically driven population dynamics under the standard statistical-mechanical premises (intrinsically scale-specific), using a parsimonious Coupled map lattice (CML) model for the simulations. In this Part II the framework will be extended with a scaling axis, orthogonal on space and time, to account for populations of individuals with space use satisfying the Multi-scaled random walk (MRW) properties. Using this scale-extended kind of CML design – the Zoomer model – I show how scale-free space use tend to generate spatial autocorrelation at the population level from the process of conspecific attraction.

“Intrinsically driven dynamics” implies that the model is simulated under the condition of a homogeneous environment. As soon as the behaviour (dynamics) under this condition is understood, environmental heterogeneity can be introduced for easier interpretation and improved realism. MRW-based behaviour at the individual level implies spatial memory and scaling under a premise of the parallel processing conjecture. As in Part I, the starting point for the present post is a citation from a previous post:

The Zoomer model  includes all the four standard BIDE rates (Birth, Immigration, Death and Emigration), and it is also spatially explicit. However, contrary to standard coupled map lattice models, spatial scale (the “lattice”) is implemented in a multi-scaled manner. This “scale range” approach allows for formulation of various aspects of complex population kinetics.

For example, the Zoomer model allows for explicit modelling of intraspecific cohesion (conspecific attraction), which is a complex process that depends on both temporal and spatial memory utilization at the individual level. In a simplified scenario, consider that a tendency for conspecific attraction is the main driver of the population kinetics. Further, consider that we study the process at sufficiently fine temporal and spatial scales to allow us to disregard the “slower” rate BIDE terms (it will be trivial to account for these processes as well).

I refer to the link above, and my book, for mathematical details of the model formulation. Below I take previous presentations a step further by extending the Zoomer model with fine-scale diffusion and a weak level of birth and death rates. Diffusion represents local randomization of individuals, based on memory-less, Markovian compliant moves. In the present simulations I assume that this kind of standard-model redistribution of individuals for the most part is contained at finer spatial scales than the finest grid cells in the simulation arena. Below I explore two scenaria; 0% and 1% nearest neighbour diffusion rate, respectively, at the unit scale (smallest grid cells) during a time increment. I also keep the weak growth rate of 1% (survival rate 99% and reproduction rate 2%), as in the Part I simulations under the standard CML model.

A sardinian warbler Sylvia melanocephala (Southern Spain) has tactical focus on chasing insects, watching out for raptors and relating to other inputs within its current field of perception. Simultaneously, the present model assumes that the warbler is also “background processing” past experiences in a manner which makes it able to relate to its environment over a wide range of spatial and temporal scales (spatial map utilization). One such environmental cue could be conspecific whereabouts. Photo: AOG.

Zero diffusion implies that any random walk-like movement (in a statistical-mechanical sense) of individuals takes place inside the unit grid cells only, on a 32×32 cell arena. One percent diffusion to nearest-neighbour cells makes this process substantially weaker than the zooming process from complex dynamics, where 5% of individuals at each time increment are “reshuffled” with equal weight (distributed equally) at respective scales 21, 22, 23, 24 and 25 relative to the smallest grid cells; i.e., 1% pr. scale level beyond the unit scale 20=1. This difference in model design makes sense, given the recent empirical confirmation of complex animal space use at the individual level, over a wide range of taxa and ecological conditions.

While the “zoomers” during a given time increment are distributed equally over the actual scale range, they are “zooming in” from scale level k+1 to to their neighbourhood at respective scales k-1 where the local population density (from that scale’s perspective) is highest at that moment. Hence, it is assumed that the individuals have the cognitive capacity to interpret respective scales’ general level of conspecific density with the help of historic experience from the environment (as supported by empirical tests and anecdotal support of the MRW model). In other words, the memory map is providing the raw data for these individual considerations what neighbourhood to zoom in to.

An important aspect of the zooming process is that individuals moving (“zooming out” from a locally high density region will have a tendency (statistically) to “zoom in” towards the same region, while individuals starting from a relatively low density region will tend to zoom into the neighbourhood with higher density.

In summary, standard diffusion-compliant reshuffling of individuals (tactical, memory-independent behaviour) is constrained to very fine spatial resolutions, while multi-scaled reshuffling is distributed over a scale range covered by the model’s grain to extent range.

The intriguing property of this construction is model coherence between scale-free distribution of displacement lengths at the individual level (a power law, with tail exponent β ≈ -2; i.e., Lévy walk-like if we ignore the spatial memory aspect) and scale-free inter-scale redistribution of the density surface at the population level. In my book and in several blog posts I have provided many examples supporting the present model assumption that the tail of individual movement lengths can span a decently large population range (for example, recall this post for the lesser kestrel example, and this post for the Florida snail kite).

Since the zooming rate between scale k and scale k+1 is 1% for all zoomer scales 1, 2, …, k, …, 5 included in the present model, this implies that individuals during an increment on average has a 1% chance of moving a distance of magnitude (k+1) relative to distance (k); which contains 22 = 4 times as large cells. However, this lower probability for long-distance moves is compensated for by (21)2 = 4 times as many individuals embedded in a k+1 grid cell than a cell at scale k (Gautestad and Mysterud, 2005). Hence, a constant cross-scale zoomer rate at the population level reflects power exponent β ≈ -2 in the step length distribution at the individual level, under the assumption that individuals on average are putting “equal weight” into space use at respective scales.

Further, in the present illustrative simulations with weight on intrinsic population dynamics the main driver for spatial redistribution from zoooming is a tendency for conspecific attraction. Avoiding getting lost in space relative to the local population is obviously an important driver for animal space use. Zoomers to a specific level k+1 during Δt are assumed to prefer to move to the highest local population density in the individual’s neighbourhood at scale k. Again, I have to refer to my papers, my book and other blog posts for details.

Finally, let’s turn to the Zoomer results. Activating zooming to all levels beyond the unit scale, and setting diffusion rate to 0%, leads to the following snapshots after a couple hundred iterations:


Setting diffusion rate to 1% and confined to scale k=1, leads to a somewhat smoother density surface:


As was the case for a non-zooming condition (Part I), also the Zoomer model shows a steep log(M,V) plot, which reflects a high degree of spatial autocorrelation in the density surface of the population. Small-scale diffusion tend to smoothen the density surface and reduce the y-intercept [log(a) becomes somewhat smaller]. However, as long as the population is by the most part undisturbed by perturbations from the local carrying capacity ceiling the power exponent remains close to b=2, as for the non-complex CML scenaria in Part I.

However, there is a specific difference between auto-correlation from a scale-specific (CML model) scenario and a scale-free (Zoomer model) scenario! In the next post, Part III, this intrinsic difference in population kinetics will be explored.

The main take-home message from this Part II is the inclusion of conspecific attraction at “strategic” scales; i.e., at spatial resolutions coarser than the individuals’ field of perception. Such a behaviour is well documented, but it has been problematic to implement in traditional model designs due to (a) the classic statistical mechanical framework is void of spatial memory effects, which strategic conspecific attraction depends on; and (b) conspecific attraction should be implemented in a scale-free manner to make the model coherent with complex space use at the individual level.


Gautestad, A. O., and I. Mysterud. 2005. Intrinsic scaling complexity in animal dispersion and abundance. The American Naturalist 165:44-55.