Source author record

Mark A. Lewis

Mark A. Lewis appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

16works
7topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

16 published item(s)

preprint2022arXiv

A Fisher-KPP model with a nonlocal weighted free boundary: analysis of how habitat boundaries expand, balance or shrink

In this paper, we propose a novel free boundary problem to model the movement of single species with a range boundary. The spatial movement and birth/death processes of the species found within the range boundary are assumed to be governed by the classic Fisher-KPP reaction-diffusion equation, while the movement of a free boundary describing the range limit is assumed to be influenced by the weighted total population inside the range boundary and is described by an integro-differential equation. Our free boundary equation is a generalization of the classical Stefan problem that allows for nonlocal influences on the boundary movement so that range expansion and shrinkage are both possible. In this paper we prove that the new model is well posed and possesses steady state. We show that the spreading speed of the range boundary is smaller than that for the equivalent problem with a Stefan condition. This implies that the nonlocal effect of the weighted total population on the boundary movement slows down the spreading speed of the population. While the classical Stefan condition categorizes asymptotic behavior via a spreading-vanishing dichotomy, the new model extends this dichotomy to a spreading-balancing-vanishing trichotomy. We specifically analyze how habitat boundaries expand, balance or shrink. When the model is extended to have two free boundaries, we observe the steady state scenario, asymmetric shifts, or even boundaries moving synchronously in the same direction. These are newly discovered phenomena in the free boundary problems for animal movement.

preprint2022arXiv

Calculating the timing and probability of arrival for sea lice dispersing between salmon farms

Sea lice are a threat to the health of both wild and farmed salmon and an economic burden for salmon farms. With a free living larval stage, sea lice can disperse tens of kilometers in the ocean between salmon farms, leading to connected sea lice populations that are difficult to control in isolation. In this paper we develop a simple analytical model for the dispersal of sea lice between two salmon farms. From the model we calculate the arrival time distribution of sea lice dispersing between farms, as well as the level of cross-infection of sea lice. We also use numerical flows from a hydrodynamic model, coupled with a particle tracking model, to directly calculate the arrival time of sea lice dispersing between two farms in the Broughton Archipelago, BC, in order to fit our analytical model and find realistic parameter estimates. Using the parametrized analytical model we show that there is often an intermediate inter-farm spacing that maximizes the level of cross infection between farms, and that increased temperatures will lead to increased levels of cross infection.

preprint2022arXiv

Detecting seasonal episodic-like spatiotemporal memory patterns using animal movement modelling

1. Spatial memory plays a role in the way animals perceive their environments, resulting in memory-informed movement patterns that are observable to ecologists. Developing mathematical techniques to understand how animals use memory in their environments allows for an increased understanding of animal cognition. 2. Here we describe a model that accounts for the memory of seasonal or ephemeral qualities of an animal's environment. The model captures multiple behaviors at once by allowing for resource selection in the present time as well as long-distance navigations to previously visited locations within an animal's home range. 3. We performed a set of analyses on simulated data to test our model, determining that it can provide informative results from as little as one year of discrete-time location data. We also show that the accuracy of model selection and parameter estimation increases with more location data. 4. This model has potential to identify a specific mechanism in which animals use memory to optimize their foraging, by revisiting temporally and predictably variable resources at consistent time lags.

preprint2022arXiv

Simulating how animals learn: a new modelling framework applied to the process of optimal foraging

Animal learning has interested ecologists and psychologists for over a century. Mathematical models that explain how animals store and recall information have gained attention recently. Central to this work is statistical decision theory (SDT), which relates information uptake in animals to Bayesian inference. SDT effectively explains many learning tasks in animals, but extending this theory to predict how animals will learn in changing environments still poses a challenge for ecologists. We addressed this shortcoming with a novel implementation of Bayesian Markov Chain Monte Carlo (MCMC) sampling to simulate how animals sample environmental information and learn as a result. We applied our framework to an individual-based model simulating complex foraging tasks encountered by wild animals. Simulated ``animals" learned behavioral strategies that optimized foraging returns simply by following the principles of an MCMC sampler. In these simulations, behavioral plasticity was most conducive to efficient foraging in unpredictable and uncertain environments. Our model suggests that animals prioritize highly concentrated resources even when these resources are less available overall, in line with existing knowledge on optimal foraging and ideal free distribution theory. Our innovative computational modelling framework can be applied more widely to simulate the learning of many other tasks in animals and humans.

preprint2021arXiv

A framework for studying transients in marine metapopulations

Transient dynamics can often differ drastically from the asymptotic dynamics of systems. In this paper we provide a unifying framework for analysing transient dynamics in marine metapopulations, from the choice of norms to the addition of stage structure. We use the $\ell_1$ norm, because of its biological interpretation, to extend the transient metrics of reactivity and attenuation to marine metapopulations, and use examples to compare these metrics under the more commonly used $\ell_2$ norm. We then connect the reactivity and attenuation of marine metapopulations to the source-sink distribution of habitat patches and demonstrate how to meaningfully measure reactivity when metapopulations are stage-structured.

preprint2020arXiv

A hybrid gravity and route choice model to assess vector traffic in large-scale road networks

Human traffic along roads can be a major vector for infectious diseases and invasive species. Though most road traffic is local, a small number of long-distance trips can suffice to move an invasion or disease front forward. Therefore, understanding how many agents travel over long distances and which routes they choose is key to successful management of diseases and invasions. Stochastic gravity models have been used to estimate the distribution of trips between origins and destinations of agents. However, in large-scale systems it is hard to collect the data required to fit these models, as the number of long-distance travellers is small, and origins and destinations can have multiple access points. Therefore, gravity models often provide only relative measures of the agent flow. Furthermore, gravity models yield no insights into which roads agents use. We resolve these issues by combining a stochastic gravity model with a stochastic route choice model. Our hybrid model can be fitted to survey data collected at roads that are used by many long-distance travellers. This decreases the sampling effort, allows us to obtain absolute predictions of both vector pressure and pathways, and permits rigorous model validation. After introducing our approach in general terms, we demonstrate its benefits by applying it to the potential invasion of zebra and quagga mussels (Dreissena spp.) to the Canadian province British Columbia. The model yields an R-squared value of 0.73 for variance-corrected agent counts at survey locations.

preprint2016arXiv

A unifying framework for quantifying the nature of animal interactions

Collective phenomena, whereby agent-agent interactions determine spatial patterns, are ubiquitous in the animal kingdom. On the other hand, movement and space use are also greatly influenced by the interactions between animals and their environment. Despite both types of interaction fundamentally influencing animal behaviour, there has hitherto been no unifying framework for the models proposed in both areas. Here, we construct a general method for inferring population-level spatial patterns from underlying individual movement and interaction processes, a key ingredient in building a statistical mechanics for ecological systems. We show that resource selection functions, as well as several examples of collective motion models, arise as special cases of our framework, thus bringing together resource selection analysis and collective animal behaviour into a single theory. In particular, we focus on combining the various mechanistic models of territorial interactions in the literature with step selection functions, by incorporate interactions into the step selection framework and demonstrating how to derive territorial patterns from the resulting models. We demonstrate the efficacy of our model by application to a population of insectivore birds in the Amazon rainforest.

preprint2016arXiv

Climate change and integrodifference equations in a stochastic environment

Climate change impacts population distributions, forcing some species to migrate poleward if they are to survive and keep up with the suitable habitat that is shifting with the temperature isoclines. Previous studies have analyzed whether populations have the capacity to keep up with shifting temperature isoclines, and have mathematically determined the combination of growth and dispersal that is needed to achieve this. However, the rate of isocline movement can be highly variable, with much uncertainty associated with yearly shifts. The same is true for population growth rates. Growth rates can be variable and uncertain, even within suitable habitats for growth. In this paper we reanalyze the question of population persistence in the context of the uncertainty and variability in isocline shifts and rates of growth. Specifically, we employ a stochastic integrodifference equation model on a patch of suitable habitat that shifts poleward at a random rate. We derive a metric describing the asymptotic growth rate of the linearized operator of the stochastic model and yielding a threshold criterion for population persistence. We demonstrate that the variability in the yearly shift and in the growth rate has a significant negative effect on the persistence in the sense that it decreases the threshold criterion for population persistence. Mathematically, we show how the persistence metric can be connected to the principal eigenvalue problem for a related integral operator, at least for the case where isocline shifting speed is deterministic. Analysis of dynamics for the case where the dispersal kernel is Gaussian leads to the existence of a critical shifting speed. Finally we illustrate our different results for butterfly population using numerical simulations and demonstrate how increased variances in isocline shifts and growth rates translate into decreased likelihoods of persistence.

preprint2016arXiv

State-space models' dirty little secrets: even simple linear Gaussian models can have estimation problems

State-space models (SSMs) are increasingly used in ecology to model time-series such as animal movement paths and population dynamics. This type of hierarchical model is often structured to account for two levels of variability: biological stochasticity and measurement error. SSMs are flexible. They can model linear and nonlinear processes using a variety of statistical distributions. Recent ecological SSMs are often complex, with a large number of parameters to estimate. Through a simulation study, we show that even simple linear Gaussian SSMs can suffer from parameter- and state-estimation problems. We demonstrate that these problems occur primarily when measurement error is larger than biological stochasticity, the condition that often drives ecologists to use SSMs. Using an animal movement example, we show how these estimation problems can affect ecological inference. Biased parameter estimates of a SSM describing the movement of polar bears (\textit{Ursus maritimus}) result in overestimating their energy expenditure. We suggest potential solutions, but show that it often remains difficult to estimate parameters. While SSMs are powerful tools, they can give misleading results and we urge ecologists to assess whether the parameters can be estimated accurately before drawing ecological conclusions from their results.

preprint2015arXiv

A mathematical approach to territorial pattern formation

Territorial behaviour is widespread in the animal kingdom, with creatures seeking to gain parts of space for their exclusive use. It arises through a complicated interplay of many different behavioural features. Extracting and quantifying the processes that give rise to territorial patterns requires both mathematical models of movement and interaction mechanisms, together with statistical techniques for rigorously extracting parameters from data. Here, we give a brisk, pedagogical overview of the techniques so far developed to tackle the problem of territory formation. We give some examples of what has already been achieved using these techniques, together with pointers as to where we believe the future lies in this area of study. This progress is a single example of a major aim for 21st century science: to construct quantitatively predictive theory for ecological systems.

preprint2015arXiv

Differentiating the Lévy walk from a composite correlated random walk

1. Understanding how to find targets with very limited information is a topic of interest in many disciplines. In ecology, such research has often focused on the development of two movement models: i) the Lévy walk and; ii) the composite correlated random walk and its associated area-restricted search behaviour. Although the processes underlying these models differ, they can produce similar movement patterns. Due to this similarity and because of their disparate formulation, current methods cannot reliably differentiate between these two models. 2. Here, we present a method that differentiates between the two models. It consists of likelihood functions, including one for a hidden Markov model, and associated statistical measures that assess the relative support for and absolute fit of each model. 3. Using a simulation study, we show that our method can differentiate between the two search models over a range of parameter values. Using the movement data of two polar bears (\textit{Ursus maritimus}), we show that the method can be applied to complex, real-world movement paths. 4. By providing the means to differentiate between the two most prominent search models in the literature, and a framework that could be extended to include other models, we facilitate further research into the strategies animals use to find resources.

preprint2015arXiv

Integrated step selection analysis: bridging the gap between resource selection and animal movement

A resource selection function is a model of the likelihood that an available spatial unit will be used by an animal, given its resource value. But how do we appropriately define availability? Step-selection analysis deals with this problem at the scale of the observed positional data, by matching each used step (connecting two consecutive observed positions of the animal) with a set of available steps randomly sampled from a distribution of observed steps or their characteristics. Here we present a simple extension to this approach, termed integrated step-selection analysis (iSSA), which relaxes the implicit assumption that observed movement attributes (i.e. velocities and their temporal autocorrelations) are independent of resource selection. Instead, iSSA relies on simultaneously estimating movement and resource-selection parameters, thus allowing simple likelihood-based inference of resource selection within a mechanistic movement model. We provide theoretical underpinning of iSSA, as well as practical guidelines to its implementation. Using computer simulations, we evaluate the inferential and predictive capacity of iSSA compared to currently used methods. Our work demonstrates the utility of iSSA as a general, flexible and user-friendly approach for both evaluating a variety of ecological hypotheses, and predicting future ecological patterns.

preprint2015arXiv

Predicting local and non-local effects of resources on animal space use using a mechanistic step selection model

1. Predicting space use patterns of animals from their interactions with the environment is fundamental for understanding the effect of habitat changes on ecosystem functioning. Recent attempts to address this problem have sought to unify resource selection analysis, where animal space use is derived from available habitat quality, and mechanistic movement models, where detailed movement processes of an animal are used to predict its emergent utilization distribution. Such models bias the animal's movement towards patches that are easily available and resource-rich, and the result is a predicted probability density at a given position being a function of the habitat quality at that position. However, in reality, the probability that an animal will use a patch of the terrain tends to be a function of the resource quality in both that patch and the surrounding habitat. 2. We propose a mechanistic model where this non-local effect of resources naturally emerges from the local movement processes, by taking into account the relative utility of both the habitat where the animal currently resides and that of where it is moving. We give statistical techniques to parametrize the model from location data and demonstrate application of these techniques to GPS data of caribou in Newfoundland. 3. Steady-state animal probability distributions arising from the model have complex patterns that cannot be expressed simply as a function of the local quality of the habitat. In particular, large areas of good habitat are used more intensively than smaller patches of equal quality habitat, whereas isolated patches are used less frequently. 4. Whilst we focus on habitats in this study, our modelling framework can be readily used with any environmental covariates and therefore represents a unification of mechanistic modelling and step selection approaches to understanding animal space use.

preprint2015arXiv

Robustness of movement models: can models bridge the gap between temporal scales of data sets and behavioural processes?

Discrete-time random walks and their extensions are common tools for analyzing animal movement data. In these analyses, resolution of temporal discretization is a critical feature. Ideally, a model both mirrors the relevant temporal scale of the biological process of interest and matches the data sampling rate. Challenges arise when resolution of data is too coarse due to technological constraints, or when we wish to extrapolate results or compare results obtained from data with different resolutions. Drawing loosely on the concept of robustness in statistics, we propose a rigorous mathematical framework for studying movement models' robustness against changes in temporal resolution. In this framework, we define varying levels of robustness as formal model properties, focusing on random walk models with spatially-explicit component. With the new framework, we can investigate whether models can validly be applied to data across varying temporal resolutions and how we can account for these different resolutions in statistical inference results. We apply the new framework to movement-based resource selection models, demonstrating both analytical and numerical calculations, as well as a Monte Carlo simulation approach. While exact robustness is rare, the concept of approximate robustness provides a promising new direction for analyzing movement models.

preprint2014arXiv

A generalized residual technique for analyzing complex movement models using earth mover's distance

1. Complex systems of moving and interacting objects are ubiquitous in the natural and social sciences. Predicting their behavior often requires models that mimic these systems with sufficient accuracy, while accounting for their inherent stochasticity. Though tools exist to determine which of a set of candidate models is best relative to the others, there is currently no generic goodness-of-fit framework for testing how close the best model is to the real complex stochastic system. 2. We propose such a framework, using a novel application of the Earth mover's distance, also known as the Wasserstein metric. It is applicable to any stochastic process where the probability of the model's state at time $t$ is a function of the state at previous times. It generalizes the concept of a residual, often used to analyze 1D summary statistics, to situations where the complexity of the underlying model's probability distribution makes standard residual analysis too imprecise for practical use. 3. We give a scheme for testing the hypothesis that a model is an accurate description of a data set. We demonstrate the tractability and usefulness of our approach by application to animal movement models in complex, heterogeneous environments. We detail methods for visualizing results and extracting a variety of information on a given model's quality, such as whether there is any inherent bias in the model, or in which situations it is most accurate. We demonstrate our techniques by application to data on multi-species flocks of insectivore birds in the Amazon rainforest. 4. This work provides a usable toolkit to assess the quality of generic movement models of complex systems, in an absolute rather than a relative sense.

preprint2014arXiv

Step selection techniques uncover the environmental predictors of space use patterns in flocks of Amazonian birds

1. Anthropogenic actions cause rapid ecological changes, meaning that animals have to respond before they have time to adapt. Tools to quantify emergent spatial patterns from animal-habitat interaction mechanisms are vital for predicting the population-level effects of such changes. 2. Environmental perturbations are particularly prevalent in the Amazon rainforest, and have a profound effect on fragmentation-sensitive insectivorous bird flocks. Therefore it is important to be able to predict the effects of such changes on the flocks' space-use patterns. 3. We use a step selection function (SSF) approach to uncover environmental drivers behind movement choices. This is used to construct a mechanistic model, from which we derive predicted utilization distributions (home ranges) of flocks. 4. We show that movement decisions are significantly influenced by canopy height and topography, but not resource depletion and renewal. We quantify the magnitude of these effects and demonstrate that they are helpful for understanding various heterogeneous aspects of space use. We compare our results to recent analytic derivations of space use, demonstrating that they are only accurate when assuming that there is no persistence in the animals' movement. 5. Our model can be translated into other environments or hypothetical scenarios, such as those given by proposed future anthropogenic actions, to make predictions of spatial patterns in bird flocks. Furthermore, our approach is quite general, so could be used to predict the effects of habitat changes on spatial patterns for a wide variety of animal communities.