Abstract

We propose that changes in movement behavior may be a proximate mechanism that influences the accumulation of animals at habitat edges. We tested this idea with a combination of empirical and simulation experiments in a resource-free landscape. The movements of individual flour beetles, Tribolium confusum, were tracked across a paper arena edged with invisible tape until beetles crossed the edge. Movement behavior (step lengths and turn angles) and cumulative occupancy were analyzed according to distance from the edge. We found that beetles took smaller steps with larger turn angles near edges than in the center of the arena and that beetle distribution was highly biased towards the edge of the arena. We then tested two agent-based simulation models for each beetle: an edge-independent model and an edge-dependent model. Both models predicted less time spent at the edge than was observed. The proportion of time spent at edges depended on the propensity to cross the edge, which could not be explained by beetle body size or energetic condition. The distribution of animals with respect to habitat edges depends on many factors, but we suggest that proximate mechanisms such as movement behavior should be explicitly considered when interpreting animal distributions.

1. Introduction

Movement of individuals is a critical factor for the population ecology of most organisms, affecting energy flow, distribution, and the genetic and demographic structure of populations [14]. Knowledge of the causes of movement behavior may be key in developing a full understanding of the spatial structure and dynamics of populations [1, 510]. For example, an increasingly large body of research finds that organisms generally move quickly and directly through low-quality habitats (e.g., [4, 11, 12]). This in turn predicts there will be fewer individuals in resource-poor than in resource-rich environments [5, 13], and many studies have found such a correspondence between movement behavior and population density (e.g., [1418]).

In addition to habitat quality, another key factor affecting the movement of organisms is the degree of permeability of the interface between habitats [16, 1921]. Permeability is the tendency of the edge to inhibit or enhance organisms’ movement across it [19, 2123] and it directly affects the degree to which organisms leave a particular habitat [12, 15, 19]. Borders that are readily crossed by dispersing individuals are considered to be soft or semipermeable edges, while those that are unlikely to be crossed are known as hard or impermeable. This distinction may be a result of the environment (a physical barrier, for instance), or a behavioral response from the animal itself [12, 16, 19, 24]. Behavioral responses may vary inter- and intraspecifically according to species’ traits (e.g., specificity of diet, migratory status, and scale of movement) and according to the characteristics of individuals (e.g., age, sex, reproductive state, and condition). Individuals’ movement patterns at edges may also vary as a result of personality: the propensity of animals to explore or take on risks can be governed by whether or not they are innately bold [25, 26], aggressive, active, or social [27].

While the notions of habitat-dependent movement behavior and edge permeability are each well established, links between them remain unclear although long recognized as of interest [28]. We suggest they may share common mechanistic determinants of movement. A widely applied way of assessing changes in movement behavior is to divide a movement path into a series of step lengths and turn angles [29]. In this model, longer step lengths and a narrower distribution of turn angles characterize directed movement, whereas shorter step lengths and a wider range of turn angles typify more random movement. As noted above, directed movement can cause animals to quickly leave habitats with few resources, while more random movement behavior results in longer durations in an area. If encounters with an edge cause individuals to shorten their step lengths or increase their turn angles, then we should see an increased density of individuals at the edge. Previously, high population density at edges has been attributed to the presence of complementary or supplementary resources in two adjacent habitats, or to higher habitat quality at the edge itself (e.g., [30, 31]). We suggest that the pattern could arise from altered movement behavior at edges regardless of resource distribution.

The idea that movement behavior alone could explain that accumulation at edges has been proposed occasionally (e.g., [32, 33]), but to our knowledge there have been no empirical tests within patches where resources are known to be homogeneous. Some studies have found evidence of accumulation at edges without an associated increase in resources, but have not investigated the mechanism behind the accumulation. For example, Olofsson et al. [34] found that reindeer grazed vegetation more intensely within 50 m of a fence placed across continuous habitat than they grazed further away. Because there is no obvious change in intrinsic habitat quality with distance to the fence, this observation could be consistent with animals accumulating near the fence because of their movement behavior. However, this movement-based hypothesis has not been explicitly tested.

To examine how movement behavior at edges might explain distribution at edges, we analyzed the movement paths of the confused flour beetle, Tribolium confusum Jacquelin du Val (Tenebrionidae) in a resource-free arena (or “patch”) bounded by a distinct but crossable edge. Tenebrionid beetles have previously been used to investigate animal responses to landscape patchiness because of their small size and propensity to move in only two dimensions (e.g., [3, 7, 35]). Furthermore, a recent study of edge responses in beetles found that nearly 60% of the 78 most common species encountered in New Zealand exhibited a significant change in abundance within two meters of the edge [36], which makes them an appropriate taxon for edge studies at small scales. We asked three main questions. First, are individuals equally likely to be in any part of the patch, or is their distribution biased with respect to the edge? Second, do movement behaviors, measured as step length and turn angles, differ between the edge and the center of the patch? And third, can beetle distribution in a bounded patch be explained simply in terms of step lengths and turn angles, or might additional mechanisms such as edge permeability or individual differences in behavior play a role?

The use of simple models in studies of movement ecology can facilitate the generalization of findings across taxa and habitats [37], but they are rarely tested with empirical data [38]. To answer the questions we posed above, we first observed the distribution and movement behaviors of individual beetles with respect to their distance from the edge of the experimental arena. We then developed an individual-based computer simulation that mimicked the conditions of the empirical experiment. The model drew step lengths and turn angles from the distributions that resulted from the experiment and simulated beetle movement using a correlated random walk. We predicted that if step length and turn angles determine beetle distribution with respect to the edge, we would be able to recreate the distributions observed in the empirical experiments in the simulation models.

2. Methods

2.1. Empirical Study

Tribolium confusum (3–5 mm in length, 0.5–1 mm in width) were obtained from a laboratory population that has been reared for fifteen years at the University of Calgary. Beetles were housed in a container of dry white wheat flour at room temperature. The experimental arena consisted of a 20 cm × 20 cm sheet of paper attached to a board on all four sides with 19 mm wide Scotch Magic tape (a semipermeable “edge” that obstructed but did not prevent movement across it). The tape differed in texture from the paper, and although it provided no vertical barrier, Morales and Ellner [35] have shown that T. confusum are unlikely to cross such tape. The arena was marked into seven concentric zones of 1.5 cm in width, which were numbered from the center outwards (Figure 1). The outer zone (zone 7a) also contained corner zones, each of which was 2.25 cm2 in area. These four corner areas were collectively known as zone 7b. Our zones allowed us to quantify time spent near and far from the edge of the arena in such a way that could be compared to a simulation model.

We conducted the experiment in a darkroom with a single incandescent light (60 W) placed 60 cm above the arena to evenly illuminate it. The light was turned on 15 min prior to the experiment to allow the temperature to stabilize; the mean temperature of trials was °C. One beetle was observed at a time, and a new experimental arena was used (Figure 2) for each trial to eliminate the potential confounding effect of residual chemical cues from previous beetles [39]. Before each trial began, the test beetle was placed in the center of the arena covered by a transparent plastic vial (2.5 cm diameter) for one minute to reduce agitation dispersal [40, 41]. The vial encompassed all of zone 1 and extended into zone 2; therefore, for the purposes of analysis, zone 1 was incorporated into zone 2. Each trial was recorded using a digital camera (DCR-TRV 230, Sony Ltd.) (60 frames sec−1) fixed directly above the arena. The trial began immediately after the one-minute acclimation period and ended after the beetle had left the experimental arena ( ) or after 600 seconds ( ), whichever occurred first. At the end of each trial, the beetle under observation was weighed to the nearest 0.1 mg and its length and width measured with a dissecting scope fitted with an ocular micrometer at 10X magnification. Beetle condition index was calculated as the residual of the regression of mass on volume (estimated as a cylinder) [42].

To determine the distribution of time spent in each zone from the digital videos of the trials, we used an event recorder (Jwatcher 1.0: http://www.jwatcher.ucla.edu/). Time spent in each zone was summed for each individual. We recorded edge encounter events as the extension of the beetle’s head past the interior edge of the transparent tape. To map the paths of beetles, we used motion analysis software (Tracker 1.7.2: http://www.cabrillo.edu/~dbrown/tracker/) to determine each beetle’s coordinates at 1-second intervals from which step lengths and turn angles were calculated. Each turn angle was recorded as the angle of change in direction from the previous step. These angles were bounded between 0 and 180°. We analyzed the first and last 30 seconds of each trial but found little difference in step lengths and turn angles between the two time periods. We therefore pooled data from both periods in our analyses. We used JMP 9.0 (SAS Inc. statistical software) and Oriana 2 (http://www.kovcomp.co.uk/oriana/newver2.html) to analyze the distribution of step lengths, and turn angles respectively.

2.2. Simulation Models

We modeled each beetle’s movement using agent-based software (Netlogo 3.1.3: http://ccl.northwestern.edu/netlogo/). The conditions of the empirical experiment were replicated as closely as possible—beetles in the model began the trial in the center of an arena with the same dimensions and zones as the experimental arena. We modeled the movement of each individual beetle in our empirical trial as a correlated random walk [29, 43]. Beetles in our simulation model took a step each second, and the length of the step and turn angle relative to the previous step were chosen from the distribution of observed movement behaviors of each particular beetle. Therefore, distribution parameters were unique to each beetle. The observed distributions of step lengths and turn angles were consistent with a normal distribution when square-root and ln-transformed, respectively. Values were randomly drawn from the transformed distributions and then back-transformed; these angles were then applied relative to the beetle’s previous turn angle to move the beetle its next step.

As with the empirical trials, simulated trials ended either after a beetle had left the arena or 600 time steps (seconds) had elapsed. We initially ran two versions of the model. The edge-independent model assumed that beetles did not alter their movement behavior based on their location within the arena. For this model, distributions were determined for the step lengths and turn angles of each beetle over all zones. The edge-dependent model included different movement behaviors for the center zones (zones 2–6), the edge zone (zone 7a), and the corner zone (zone 7b) according to the observed distributions of step lengths and turn angles in these zones for each beetle. We combined center zones because we observed little difference in the movement behavior of beetles among these zones (analyses not shown). Each version of the model was replicated 50 times for each of the 54 beetles. The output of each replicate was the number of seconds the simulated beetle was in each of the 7 zones.

2.3. Data Analysis

We employed an individual-based analysis to examine the movement behavior of beetles in our empirical trials and in our models. We calculated the time spent in each zone for each beetle in the experimental study and for both models and transformed ( , where summed time) these values meet the assumptions of parametric tests, including normality, homogeneity of variance, and independent errors. The proportion of time individual beetles spent in each zone was then compared to the proportion of time expected based on the area of that zone (calculated as trial length/zone area) with a one-sample -test. We used simple linear regression to consider the effect of individual beetle and zone of the arena (center, edge, and corner) on ln-transformed step lengths and turn angles. We compared the cumulative time spent at the edge zone in the empirical study and in the models using mixed linear models with individual beetle as a random effect. Finally, to address the influence of additional factors on beetle distribution across the experimental arena, we used multiple linear regression to determine whether beetle traits (body volume, and condition index), beetle movement behavior, and the response to edge (crosses per encounter) explained the proportion of time an individual spent at the edge. Statistical significance of regression analyses was determined via -tests. For all analyses, we used a significance value of as the criterion for rejecting the null hypothesis.

3. Results

3.1. Empirical Study
3.1.1. Distribution within the Arena

Observed trial lengths varied from 8.4 to the imposed maximum of 600 sec (median = 47.0 sec, ). The time beetles spent in each zone was not proportional to the area of each zone (Figure 3(a)). Beetles spent significantly more time than expected at the edge (zone 7a: , ) and corners (zone 7b: , ). Conversely, they spent less time than expected in the middle zones (zone 4: , zone 5: , and zone 6: ; all: ). Time spent in the innermost zones was proportional to those areas of the arena (zone 2: , zone 3: , ).

3.1.2. Movement Behavior and Distribution at Edges

Beetles took smaller steps (Figure 4(a)) with larger turn angles (Figure 4(b)) in the edge zone (7a) than in the center zones (zones 2–6). The mean difference in step length (95% CI) was 1.36 mm (1.25–1.49 mm, back-transformed from analysis of ln-transformed data; , ), while turn angles decreased by 11.9° on average (95% CI: 5.6–18.2 o; , ). Individual differences were evident in that mean step lengths at the edge were positively correlated with mean step lengths in the center (Figure 4(a); , ), but there was no consistency within individuals in the relative magnitude of turn angles at the center and edge (Figure 4(b); , ).

3.2. Simulation Models

The results of our two simulation models, with edge-independent or edge-dependent behavior, differed from each other but neither replicated the distribution of time spent in each zone by beetles in our experimental study (Figures 3(b) and 3(c) versus 3(a)). The edge-dependent model resulted in a higher proportion of time being spent at the edge zone (zone 7a) than the edge-independent model ( , ). Conversely, less time was spent in the innermost zones (zones 2, 3) and in zone 4 in the edge-dependent model than in the edge-independent model. However, despite differences in magnitude, both models predicted less time spent near the edge (in zone 7a) and more time spent near the center (zones 2 through 6) than was actually observed (Figure 3).

Given the above results, we remodeled beetle movement using only the step lengths and turns angles that were most different between the center and the edge observed for any beetle in the empirical trials. We wanted to explore whether any changes in movement behavior between center and edge, within the range of those observed, would result in a distribution of beetles’ use of the arena that was similar to our empirical results. As expected from our original hypothesis, the short steps and large turn angles at the edges, relative to the center, resulted in beetles spending most of their time near the edge (Figure 3(d)). However, this modification to the edge-dependent simulation model did not succeed in replicating the distribution of beetles in our empirical experiment. This version of the model produced a distribution in which twice as much time was spent in the edge zones than beetles had spent in the empirical study (80% versus 40%) (Figure 3(d)).

To address the influence of additional factors on beetle distribution across the experimental arena, we investigated whether beetle traits, beetle movement behavior, and the response to edge (crosses per encounter) explained the proportion of time an individual spent at the edge. Beetle volume and condition did not significantly affect time at the edge (both ) and were removed from subsequent models. Some movement behaviors had small but significant effects on the proportion of time at the edge. Beetles that took larger steps with smaller turn angles at the center tended to spend more time at the edge (Table 1). However, their behavior near edges did not detectably affect the time they spent there, whether measured as the change in behavior at the edge (Table 1) or when only edge step lengths and turn angles were used in models ( , analyses not shown). The strongest predictor of the proportion of time spent at the edge was the proportion of encounters with an edge that resulted in the beetle crossing an edge (Table 1). When the probability of crossing upon encountering an edge was low, the proportion of time spent at the edge was higher (Figure 5). Furthermore, the ratio of crosses per encounter was affected neither by body volume ( , ) nor by condition index ( , ).

4. Discussion

The beetles in our experiment were surrogates for any organism capable of movement. Our experimental microlandscape was resource-free and, from the beetles’ point-of-view, was comprised of two habitat elements: the paper arena and the taped edge. For beetles born and raised in a container of flour, the arena was likely to seem foreign and exposed. It is reasonable to assume that beetles were motivated to find a better habitat (see [35]) but there was no reason to think that they should persist in trying to go in any particular direction. Thus our null hypothesis was that beetles would use the whole arena equally such that the time they spend in a given zone at some distance from the edge should be proportional to the size of that zone. The taped edge had virtually no vertical dimension and was likely not detectable from any distance. However, beetles encountering this slippery surface were unlikely to cross it to leave the experimental arena. This experimental design allowed us to consider the effect of a semipermeable boundary on the movement patterns and resulting distribution of organisms in a controlled environment on a small, manageable scale.

We found that beetles spent more time at edges and corners and less time in interior zones than expected from the area of each zone (Figure 3(a)). This occurred despite the fact that beetles were initially placed in the center of the arena. Campbell and Hagstrum [44] obtained similar results for the congeneric beetle T. castaneum, in that they found that beetles tended to be located near the edges of food patches. As in our study, they found this was due in part to beetles traveling more slowly along edges than farther away from them. The observation that animals often accumulate near edges is not new ([32, 45], and references therein), and indeed certain bird and mammal species are commonly considered to be “edge species” (but see [31, 46, 47]) for critical discussion of this term]. The underlying reasons for these distribution patterns often remain unclear because of the sheer number of variables that must be taken into consideration in natural systems. Our system suggests that an edge alone, without any evident edge-associated changes in resources, can result in greater animal activity at edges.

Our study proposed the hypothesis that disproportionate activity at edges could be explained mechanistically by changes in step lengths and turn angles. We did observe the proposed mechanism, namely, shorter step lengths and higher turn angles at edges than at the center and disproportionate time spent at edges in our putatively resource-free environment. Qualitatively, these movement patterns predict increased time spent at edges, as we observed. Quantitatively, however, we did not recreate the observed distributions using mean movement behaviors in centers and at edges in two simulation models (Figures 3(b) and 3(c)) or using “extreme” movement behaviors in a third (Figure 3(d)). Thus while it is possible that a combination of turn angles and step lengths could be found that would ultimately predict the observed distributions, it seems more likely that that there are more complex edge behaviors involved. It is important to note, however, that although the mechanism tested was not sufficient to explain the distribution of beetles with respect to edges in our study, it is a process that should be considered when integrating multiple paradigms to explain animal movement [10].

Generally, beetles encountered the edge a number of times before they crossed (Figure 5). It seems plausible that there was something fundamentally different about their behavior at the edge that we did not capture with our simple model of step lengths and turn angles. This behavioral change may be related to the beetles’ perceptual range, as they appeared not to perceive the edge until they were within millimeters of it. Many studies of animal movement have found evidence of a “boundary effect” wherein organisms perceive and respond to both semipermeable and impermeable habitat edges (e.g., [1, 4, 14, 16, 48, 49]) but few studies have clearly linked edge behavior to distribution with respect to the edge.

Crossing the boundary may be a result not of movement ability so much as movement propensity, or the probability that an individual will cross a boundary once it reaches an edge [50]. Habitat edges acting as impermeable or semi-permeable boundaries may significantly affect the trajectory of movement paths in a number of species (e.g., [1, 49]). Indeed, there is increasing evidence that animals change their movement patterns upon encountering edges (e.g., [15, 5153]). Correlated random walks may not describe movement at edges well, for example, if animals “reflect” off boundaries and take a new direction at random, if animals respond to edges by linearizing their movement paths parallel to them, or if animals otherwise shift movement modes at edges (i.e., from a correlated random walk to a biased correlated random walk, sensu Crone and Schultz [54]). Future work incorporating explicit behavioral responses at edges into individual-based models could provide more insight into how animals are distributed around edges.

Despite beetles’ general propensity to avoid crossing the edge upon their first encounter, they differed in the number of times they encountered the tape before passing over it and out of the arena (Figure 5). We interpret this to mean that some beetles considered the tape to be a hard edge, while other beetles did not perceive the tape as an edge at all. It is possible that individual beetles responded differently to edges due to factors such as sex, reproductive state, condition, and personality (i.e., boldness). We found that the number of edge encounters prior to leaving was significantly related to the time beetles spent in the edge zone (Table 1). However, we did not find any correlation between beetle condition or body size and time spent in the edge zone.

The beetles in this study exhibited a considerable amount of variation not just in their likelihood to cross the edge, but also in their movement behaviors (Figure 4). However, despite this variation, individual beetles were generally consistent within themselves in that if they took large steps in the center, they took large steps at edges (Figure 4). These results support the notion that population processes may depend more on the behaviors of a few bold individuals than on an average individual [55].

Our empirical study found that animal activity can be concentrated at edges even without a change of resources there. Thus, animal distribution at edges may not always reflect resources (cf. [30]) or as a consequence, follow an ideal free distribution. Rather, the accumulation of individuals at edges may influence resources. For example, one of the consequences of higher reindeer activity near fences is altered ecosystem processes [34]. Aggregation of individuals due to movement behavior at edges may also influence community structure [33]. Understanding the complex behavior of animals at edges independent of resources will help clarify population processes, which is critical as habitats become increasingly fragmented.

Acknowledgments

The authors thank Abbie Stewart and Josée Méthot for technical assistance and Petr Komers and Elizabeth Crone for their helpful and constructive comments on this work. This work was supported by NSERC scholarships to HCY, LEL and TL, a NSERC Discovery Grant to MLR, and the University of Calgary.