Forecasting Pathogen Dynamics with Bayesian Model-Averaging: Application to Xylella fastidiosa.

Forecasting invasive-pathogen dynamics is paramount to anticipate eradication and containment strategies. Such predictions can be obtained using a model grounded on partial differential equations (PDE; often exploited to model invasions) and fitted to surveillance data. This framework allows the construction of phenomenological but concise models relying on mechanistic hypotheses and real observations. However, it may lead to models with overly rigid behavior and possible data-model mismatches. Hence, to avoid drawing a forecast grounded on a single PDE-based model that would be prone to errors, we propose to apply Bayesian model averaging (BMA), which allows us to account for both parameter and model uncertainties. Thus, we propose a set of different competing PDE-based models for representing the pathogen dynamics, we use an adaptive multiple importance sampling algorithm (AMIS) to estimate parameters of each competing model from surveillance data in a mechanistic-statistical framework, we evaluate the posterior probabilities of models by comparing different approaches proposed in the literature, and we apply BMA to draw posterior distributions of parameters and a posterior forecast of the pathogen dynamics. This approach is applied to predict the extent of Xylella fastidiosa in South Corsica, France, a phytopathogenic bacterium detected in situ in Europe less than 10 years ago (Italy 2013, France 2015). Separating data into training and validation sets, we show that the BMA forecast outperforms competing forecast approaches.

More pests but less pesticide applications: Ambivalent effect of landscape complexity on conservation biological control.

In agricultural landscapes, the amount and organization of crops and semi-natural habitats (SNH) have the potential to promote a bundle of ecosystem services due to their influence on ecological community at multiple spatio-temporal scales. SNH are relatively undisturbed and are often source of complementary resources and refuges, therefore supporting more diverse and abundant natural pest enemies. However, the nexus of SNH proportion and organization with pest suppression is not trivial. It is thus crucial to understand how the behavior of pest and natural enemy species, the underlying landscape structure, and their interaction, may influence conservation biological control (CBC). Here, we develop a generative stochastic landscape model to simulate realistic agricultural landscape compositions and configurations of fields and linear elements. Generated landscapes are used as spatial support over which we simulate a spatially explicit predator-prey dynamic model. We find that increased SNH presence boosts predator populations by sustaining high predator density that regulates and keeps pest density below the pesticide application threshold. However, predator presence over all the landscape helps to stabilize the pest population by keeping it under this threshold, which tends to increase pest density at the landscape scale. In addition, the joint effect of SNH presence and predator dispersal ability among hedge and field interface results in a stronger pest regulation, which also limits pest growth. Considering properties of both fields and linear elements, such as local structure and geometric features, provides deeper insights for pest regulation; for example, hedge presence at crop field boundaries clearly strengthens CBC. Our results highlight that the integration of species behaviors and traits with landscape structure at multiple scales is necessary to provide useful insights for CBC.

Dating and localizing an invasion from post-introduction data and a coupled reaction-diffusion-absorption model.

Invasion of new territories by alien organisms is of primary concern for environmental and health agencies and has been a core topic in mathematical modeling, in particular in the intents of reconstructing the past dynamics of the alien organisms and predicting their future spatial extents. Partial differential equations offer a rich and flexible modeling framework that has been applied to a large number of invasions. In this article, we are specifically interested in dating and localizing the introduction that led to an invasion using mathematical modeling, post-introduction data and an adequate statistical inference procedure. We adopt a mechanistic-statistical approach grounded on a coupled reaction-diffusion-absorption model representing the dynamics of an organism in an heterogeneous domain with respect to growth. Initial conditions (including the date and site of the introduction) and model parameters related to diffusion, reproduction and mortality are jointly estimated in the Bayesian framework by using an adaptive importance sampling algorithm. This framework is applied to the invasion of Xylella fastidiosa, a phytopathogenic bacterium detected in South Corsica in 2015, France.

Modelling Population Dynamics in Realistic Landscapes with Linear Elements: A Mechanistic-Statistical Reaction-Diffusion Approach.

We propose and develop a general approach based on reaction-diffusion equations for modelling a species dynamics in a realistic two-dimensional (2D) landscape crossed by linear one-dimensional (1D) corridors, such as roads, hedgerows or rivers. Our approach is based on a hybrid "2D/1D model", i.e, a system of 2D and 1D reaction-diffusion equations with homogeneous coefficients, in which each equation describes the population dynamics in a given 2D or 1D element of the landscape. Using the example of the range expansion of the tiger mosquito Aedes albopictus in France and its main highways as 1D corridors, we show that the model can be fitted to realistic observation data. We develop a mechanistic-statistical approach, based on the coupling between a model of population dynamics and a probabilistic model of the observation process. This allows us to bridge the gap between the data (3 levels of infestation, at the scale of a French department) and the output of the model (population densities at each point of the landscape), and to estimate the model parameter values using a maximum-likelihood approach. Using classical model comparison criteria, we obtain a better fit and a better predictive power with the 2D/1D model than with a standard homogeneous reaction-diffusion model. This shows the potential importance of taking into account the effect of the corridors (highways in the present case) on species dynamics. With regard to the particular case of A. albopictus, the conclusion that highways played an important role in species range expansion in mainland France is consistent with recent findings from the literature.