Inference for nonlinear epidemiological models using genealogies and time series

Phylodynamics - the field aiming to quantitatively integrate the ecological and evolutionary dynamics of rapidly evolving populations like those of RNA viruses - increasingly relies upon coalescent approaches to infer past population dynamics from reconstructed genealogies. As sequence data have become more abundant, these approaches are beginning to be used on populations undergoing rapid and rather complex dynamics. In such cases, the simple demographic models that current phylodynamic methods employ can be limiting. First, these models are not ideal for yielding biological insight into the processes that drive the dynamics of the populations of interest. Second, these models differ in form from mechanistic and often stochastic population dynamic models that are currently widely used when fitting models to time series data. As such, their use does not allow for both genealogical data and time series data to be considered in tandem when conducting inference. Here, we present a flexible statistical framework for phylodynamic inference that goes beyond these current limitations. The framework we present employs a recently developed method known as particle MCMC to fit stochastic, nonlinear mechanistic models for complex population dynamics to gene genealogies and time series data in a Bayesian framework. We demonstrate our approach using a nonlinear Susceptible-Infected-Recovered (SIR) model for the transmission dynamics of an infectious disease and show through simulations that it provides accurate estimates of past disease dynamics and key epidemiological parameters from genealogies with or without accompanying time series data.     

Matching the forecast horizon with the relevant spatial and temporal processes and data sources

Most phenomenological, statistical models used to generate ecological forecasts take either a time-series approach, based on long-term data from one location, or a space-for-time approach, based on data describing spatial patterns across environmental gradients. However, the magnitude and even the sign of environment–response relationships detected using these two approaches often differs, leading to contrasting predictions about responses to future environmental change. Here we consider how the forecast horizon determines whether more accurate predictions come from the time-series approach, the space-for-time approach or a combination of the two. As proof of concept, we use simulated case studies to show that forecasts for short and long forecast horizons need to focus on different ecological processes, which are reflected in different kinds of data. First, we simulated population or community dynamics under stationary temperature using two simple, mechanistic models. Second, we fit statistical models to the simulated data using a time-series approach, a space-for-time approach or a weighted average. We then forecast the response to a temperature increase using the statistical models, and compared these forecasts to temperature effects simulated by the mechanistic models. We found that the time-series approach made accurate short-term predictions because it captured initial conditions and effects of fast processes such as birth and death. The space-for-time approach made more accurate long-term predictions because it better captured the influence of slower processes such as evolutionary and ecological selection. The weighted average made accurate predictions at all time scales, including intermediate time-scales where the other two approaches performed poorly. A weighted average of time-series and space-for-time approaches shows promise, but making this weighted model operational will require new research to predict the rate at which slow processes begin to influence dynamics.     

When mechanism matters: Bayesian forecasting using models of ecological diffusion

Ecological diffusion is a theory that can be used to understand and forecast spatio‐temporal processes such as dispersal, invasion, and the spread of disease. Hierarchical Bayesian modelling provides a framework to make statistical inference and probabilistic forecasts, using mechanistic ecological models. To illustrate, we show how hierarchical Bayesian models of ecological diffusion can be implemented for large data sets that are distributed densely across space and time. The hierarchical Bayesian approach is used to understand and forecast the growth and geographic spread in the prevalence of chronic wasting disease in white‐tailed deer (Odocoileus virginianus). We compare statistical inference and forecasts from our hierarchical Bayesian model to phenomenological regression‐based methods that are commonly used to analyse spatial occurrence data. The mechanistic statistical model based on ecological diffusion led to important ecological insights, obviated a commonly ignored type of collinearity, and was the most accurate method for forecasting.     

Breaking out from a restricted range: Alternative habitat models to assess population perspectives

Habitat conservation for restricted-range species should also consider adjacent areas, but the analytical approaches for such assessments (particularly for a future perspective) are constrained by currently observed habitat relationships. We used two conceptually different habitat modelling approaches for analysing habitat distribution for the isolated Estonian population of a species of European conservation concern, the Siberian flying squirrel (Pteromys volans (Linnaeus, 1758)). We expected that the correlative (statistical) approaches based on current location data will increasingly deviate along with the distance from the current range, compared with a mechanistic approach based on limiting factors for the species. For conservation planning, we also investigated how the current protected area network covers quality habitats around the current range. We constructed three alternative correlative models (MaxEnt; Random forest; Generalized Boosted Regression) utilizing remote-sensing (Sentinel-2; LiDAR) and forest inventory data for 1299 occurrences in the currently occupied ca. 1400 km² range. A mechanistic model was constructed as a decision tree that distinguished 11 quality classes of forest land based on the ecological prioritization of limiting factors: site type; forest cover; abundance of key tree species; stand age; patch size; and layer structure. Supporting our expectation, an overall good accordance of habitat predictions of all the correlative models and the mechanistic model (at 30 × 30 m pixel size) declined with the distance from the current range. The MaxEnt model most closely followed the full range of habitat quality classes of the mechanistic model, while the other correlative models emphasized the highest habitat-quality class. Within the current range, both MaxEnt and the mechanistic model similarly revealed habitat quality differences between occupied and unoccupied species protection areas. Delineation of habitat aggregations all over the country based on the mechanistic model revealed habitat loss both within and adjacent to the current range, which sets limits to local population recovery. For analysing wider options, we recommend complementing statistical spatial modelling of current conditions with ecologically sound mechanistic approaches. Based on our specific case, we outline how such model predictions can be assessed for management planning beyond current range.     

Pushing the limits in marine species distribution modelling: lessons from the land present challenges and opportunities

Aim Species distribution models (SDMs) have been used to address a wide range of theoretical and applied questions in the terrestrial realm, but marine‐based applications remain relatively scarce. In this review, we consider how conceptual and practical issues associated with terrestrial SDMs apply to a range of marine organisms and highlight the challenges relevant to improving marine SDMs. Location We include studies from both marine and terrestrial systems that encompass many geographic locations around the globe. Methods We first performed a literature search and analysis of marine and terrestrial SDMs in ISI Web of Science to assess trends and applications. Using knowledge from terrestrial applications, we critically evaluate the application of SDMs in marine systems in the context of ecological factors (dispersal, species interactions, aggregation and ontogenetic shifts) and practical considerations (data quality, alternative modelling approaches and model validation) that facilitate or create difficulties for model application. Results The relative importance of ecological factors to be considered when applying SDMs varies among terrestrial and marine organisms. Correctly incorporating dispersal is frequently considered an important issue for terrestrial models, but because there is greater potential for dispersal in the ocean, it is often less of a concern in marine SDMs. By contrast, ontogenetic shifts and feeding have received little attention in terrestrial SDM applications, but these factors are important to many marine SDMs. Opportunities also exist for applying more advanced SDM approaches in the marine realm, including mechanistic ecophysiological models, where water balance and heat transfer equations are simpler for some marine organisms relative to their terrestrial counterparts. Main conclusions SDMs have generally been under‐utilized in the marine realm relative to terrestrial applications. Correlative SDM methods should be tested on a range of marine organisms, and we suggest further development of methods that address ontogenetic shifts and feeding interactions. We anticipate developments in, and cross‐fertilization between, coupled correlative and process‐based SDMs, mechanistic eco‐physiological SDMs, and spatial population dynamic models for climate change and species invasion applications in particular. Comparisons of the outputs of different model types will provide insight that is useful for improved spatial management of marine species.     

Wise use of statistical tools in ecological field studies

The currently dominating hypothetico-deductive research paradigm for ecology has statistical hypothesis testing as a basic element. Classic statistical hypothesis testing does, however, present the ecologist with two fundamental dilemmas when field data are to be analyzed: (1) that the statistically motivated demand for a random and representative sample and the ecologically motivated demand for representation of variation in the study area cannot be fully met at the same time; and (2) that the statistically motivated demand for independence of errors calls for sampling distances that exceed the scales of relevant pattern-generating processes, so that samples with statistically desirable properties will be ecologically irrelevant. Reasons for these dilemmas are explained by consideration of the classic statistical Neyman-Pearson test procedure, properties of ecological variables, properties of sampling designs, interactions between properties of the ecological variables and properties of sampling designs, and specific assumptions of the statistical methods. Analytic solutions to problems underlying the dilemmas are briefly reviewed. I conclude that several important research objectives cannot be approached without subjective elements in sampling designs. I argue that a research strategy entirely based on rigorous statistical testing of hypotheses is insufficient for field ecological data and that inductive and deductive approaches are complementary in the process of building ecological knowledge. I recommend that great care is taken when statistical tests are applied to ecological field data. Use of less formal modelling approaches is recommended for cases when formal testing is not strictly needed. Sets of recommendations, “Guidelines for wise use of statistical tools”, are proposed both for testing and for modelling. Important elements of wise-use guidelines are parallel use of methods that preferably belong to different methodologies, selection of methods with few and less rigorous assumptions, conservative interpretation of results, and abandonment of definitive decisions based a predefined significance level.     

Are there general mechanisms of animal home range behaviour? A review and prospects for future research

Home range behaviour is a common pattern of space use, having fundamental consequences for ecological processes. However, a general mechanistic explanation is still lacking. Research is split into three separate areas of inquiry - movement models based on random walks, individual-based models based on optimal foraging theory, and a statistical modelling approach - which have developed without much productive contact. Here we review recent advances in modelling home range behaviour, focusing particularly on the problem of identifying mechanisms that lead to the emergence of stable home ranges from unbounded movement paths. We discuss the issue of spatiotemporal scale, which is rarely considered in modelling studies, as well as highlighting the need to consider more closely the dynamical nature of home ranges. Recent methodological and theoretical advances may soon lead to a unified approach, however, conceptually unifying our understanding of linkages among home range behaviour and ecological or evolutionary processes.     

Ecological hierarchies and self-organisation – Pattern analysis,modelling and process integration across scales

A continuing discussion in applied and theoretical ecology focuses on the relationship of different organisational levels and on how ecological systems interact across scales. We address principal approaches to cope with complex across-level issues in ecology by applying elements of hierarchy theory and the theory of complex adaptive systems. A top-down approach, often characterised by the use of statistical techniques, can be applied to analyse large-scale dynamics and identify constraints exerted on lower levels. Current developments are illustrated with examples from the analysis of within-community spatial patterns and large-scale vegetation patterns. A bottom-up approach allows one to elucidate how interactions of individuals shape dynamics at higher levels in a self-organisation process; e.g., population development and community composition. This may be facilitated by various modelling tools, which provide the distinction between focal levels and resulting properties. For instance, resilience in grassland communities has been analysed with a cellular automaton approach, and the driving forces in rodent population oscillations have been identified with an agent-based model. Both modelling tools illustrate the principles of analysing higher level processes by representing the interactions of basic components.The focus of most ecological investigations on either top-down or bottom-up approaches may not be appropriate, if strong cross-scale relationships predominate. Here, we propose an ‘across-scale-approach’, closely interweaving the inherent potentials of both approaches. This combination of analytical and synthesising approaches will enable ecologists to establish a more coherent access to cross-level interactions in ecological systems.     

Spatial modelling of succession-disturbance dynamics in forest ecosystems: Concepts and examples

Over the last few decades it has become increasingly obvious that disturbance, whether natural or anthropogenic in origin, is ubiquitous in ecosystems. Disturbance-related processes are now considered to be important determinants of the composition, structure and function of ecological systems. However, because disturbance and succession processes occur across a wide range of spatio-temporal scales their empirical investigation is difficult. To counter these difficulties much use has been made of spatial modelling to explore the response of ecological systems to disturbance(s) occurring at spatial scales from the individual to the landscape and above, and temporal scales from minutes to centuries. Here we consider such models by contrasting two alternative motivations for their development and use: prediction and exploration, with a focus on forested ecosystems. We consider the two approaches to be complementary rather than competing. Predictive modelling aims to combine knowledge (understanding and data) with the goal of predicting system dynamics; conversely, exploratory models focus on developing understanding in systems where uncertainty is high. Examples of exploratory modelling include model-based explorations of generic issues of criticality in ecological systems, whereas predictive models tend to be more heavily data-driven (e.g. species distribution models). By considering predictive and exploratory modelling alongside each other, we aim to illustrate the range of methods used to model succession and disturbance dynamics and the challenges involved in the model-building and evaluation processes in this arena.     

Assessing temporal scales and patterns in time series: Comparing methods based on redundancy analysis

Time-series modelling techniques are powerful tools for studying temporal scaling structures and dynamics present in ecological and other complex systems and are gaining popularity for assessing resilience quantitatively. Among other methods, canonical ordinations based on redundancy analysis are increasingly used for determining temporal scaling patterns that are inherent in ecological data. However, modelling outcomes and thus inference about ecological dynamics and resilience may vary depending on the approaches used. In this study, we compare the statistical performance, logical consistency and information content of two approaches: (i) asymmetric eigenvector maps (AEM) that account for linear trends and (ii) symmetric distance-based Moran's eigenvector maps (MEM), which requires detrending of raw data to remove linear trends prior to analysis. Our comparison is done using long-term water quality data (25 years) from three Swedish lakes. This data set therefore provides the opportunity for assessing how the modelling approach used affects performance and inference in time series modelling. We found that AEM models had consistently more explanatory power than MEM, and in two out of three lakes AEM extracted one more temporal scale than MEM. The scale-specific patterns detected by AEM and MEM were uncorrelated. Also individual water quality variables explaining these patterns differed between methods, suggesting that inferences about systems dynamics are dependent on modelling approach. These findings suggest that AEM might be more suitable for assessing dynamics in time series analysis compared to MEM when temporal trends are relevant. The AEM approach is logically consistent with temporal autocorrelation where earlier conditions can influence later conditions but not vice versa. The symmetric MEM approach, which ignores the asymmetric nature of time, might be suitable for addressing specific questions about the importance of correlations in fluctuation patterns where there are no confounding elements of linear trends or a need to assess causality.     

Mechanistic niche modelling: combining physiological and spatial data to predict species' ranges

Species distribution models (SDMs) use spatial environmental data to make inferences on species' range limits and habitat suitability. Conceptually, these models aim to determine and map components of a species' ecological niche through space and time, and they have become important tools in pure and applied ecology and evolutionary biology. Most approaches are correlative in that they statistically link spatial data to species distribution records. An alternative strategy is to explicitly incorporate the mechanistic links between the functional traits of organisms and their environments into SDMs. Here, we review how the principles of biophysical ecology can be used to link spatial data to the physiological responses and constraints of organisms. This provides a mechanistic view of the fundamental niche which can then be mapped to the landscape to infer range constraints. We show how physiologically based SDMs can be developed for different organisms in different environmental contexts. Mechanistic SDMs have different strengths and weaknesses to correlative approaches, and there are many exciting and unexplored prospects for integrating the two approaches. As physiological knowledge becomes better integrated into SDMs, we will make more robust predictions of range shifts in novel or non-equilibrium contexts such as invasions, translocations, climate change and evolutionary shifts.     

Integrating modelling of biodiversity composition and ecosystem function

There is increasing reliance on ecological models to improve our understanding of how ecological systems work, to project likely outcomes under alternative global change scenarios and to help develop robust management strategies. Two common types of spatiotemporally explicit ecological models are those focussed on biodiversity composition and those focussed on ecosystem function. These modelling disciplines are largely practiced separately, with separate literature, despite growing evidence that natural systems are shaped by the interaction of composition and function. Here we call for the development of new modelling approaches that integrate composition and function, accounting for the important interactions between these two dimensions, particularly under rapid global change. We examine existing modelling approaches that have begun to combine elements of composition and function, identifying their potential contribution to fully integrated modelling approaches. The development and application of integrated models of composition and function face a number of important challenges, including biological data limitations, system knowledge and computational constraints. We suggest a range of promising avenues that could help researchers overcome these challenges, including the use of virtual species, macroecological relationships and hybrid correlative‐mechanistic modelling. Explicitly accounting for the interactions between composition and function within integrated modelling approaches has the potential to improve our understanding of ecological systems, provide more accurate predictions of their future states and transform their management. Synthesis There is increasing attention from researchers and policy makers around the world on both assessing and projecting the state of the planet's biodiversity, its ecosystems and the essential services they provide to society. However, existing modelling approaches largely ignore the interactions between biodiversity composition and ecosystem function. We highlight the key challenges and potential solutions to developing integrated models of composition and function. Such models will require a new effort and focus from ecologists, yet the benefits are likely to be substantial, including better informing the management of natural systems at regional, national and international scales.     

An open source simulation model for soil and sediment bioturbation

Bioturbation is one of the most widespread forms of ecological engineering and has significant implications for the structure and functioning of ecosystems, yet our understanding of the processes involved in biotic mixing remains incomplete. One reason is that, despite their value and utility, most mathematical models currently applied to bioturbation data tend to neglect aspects of the natural complexity of bioturbation in favour of mathematical simplicity. At the same time, the abstract nature of these approaches limits the application of such models to a limited range of users. Here, we contend that a movement towards process-based modelling can improve both the representation of the mechanistic basis of bioturbation and the intuitiveness of modelling approaches. In support of this initiative, we present an open source modelling framework that explicitly simulates particle displacement and a worked example to facilitate application and further development. The framework combines the advantages of rule-based lattice models with the application of parameterisable probability density functions to generate mixing on the lattice. Model parameters can be fitted by experimental data and describe particle displacement at the spatial and temporal scales at which bioturbation data is routinely collected. By using the same model structure across species, but generating species-specific parameters, a generic understanding of species-specific bioturbation behaviour can be achieved. An application to a case study and comparison with a commonly used model attest the predictive power of the approach.     

Efficient methods for studying stochastic disease and population dynamics

Stochastic ecological and epidemiological models are now routinely used to inform management and decision making throughout conservation and public-health. A difficulty with the use of such models is the need to resort to simulation methods when the population size (and hence the size of the state space) becomes large, resulting in the need for a large amount of computation to achieve statistical confidence in results. Here we present two methods that allow evaluation of all quantities associated with one- (and higher) dimensional Markov processes with large state spaces. We illustrate these methods using SIS disease dynamics and studying species that are affected by catastrophic events. The methods allow the possibility of extending exact Markov methods to real-world problems, providing techniques for efficient parameterisation and subsequent analysis.