Knowledge Exchange Fellowship: developing and sharing mathematical tools to protect urban trees and woodland from invasive pests

Lead Research Organisation: Newcastle University
Department Name: Sch of Maths, Statistics and Physics

Abstract

The loss of biodiversity due to the spread of tree diseases and harmful invasive pests within our native forests is having an enormous environmental, economic, and social impact. The threat has been escalating in recent years due to increased accidental international imports and climate change creating a more favourable environment for many pathogens and pests. In the '25 Year Environment Plan' the UK government highlights enhancing biosecurity as a key priority, through the control of existing diseases and pests, and by building forest resilience against new ones.

A key invasive pest of concern is the oak processionary moth (OPM), a destructive pest harmful to both oak trees and humans. Despite great efforts to contain the outbreak to the original infested area of south-east England, OPM continues to spread. Trees within public parks and urban areas are most affected, risking further exacerbating both the socio-economic and ethnic divides in the access to quality green space. The prediction and control of the future dynamics of the OPM population is the priority of the governmental OPM Control Programme (OPM-CP: including researchers from the Department for Environment, Food and Rural Affairs (DEFRA), the Forestry Commission (FC), and Southampton GeoData, among others).

Increasing the use of computational modelling and statistical epidemiology are priority measures identified by DEFRA to protect plant health for environmental and social wellbeing. Despite the potential impact of such interdisciplinary methods, there is a gap in translating these theoretical techniques into widespread usage by environmental organisations for developing robust forestry management and protection plans. This Fellowship aims to bridge this gap to develop mathematical models for the spread of OPM through dynamic collaboration with members of the OPM-CP.

Aligned with NERC's priority delivery areas of 'healthy environment' through protecting current woodland health and 'resilient environment' through informing strategies to ensure future woodland resilience, this Fellowship will drive the development of both a local agent-based model, and a national stochastic differential continuum model for the spread of OPM. The models will be underpinned by rigorous statistical analysis and parameter inference of real data shared by OPM-CP and used to inform policies to control OPM, particularly in deprived urban areas at high risk of infestation. The results will not only provide an insight into managing the control of OPM but also establish a framework for knowledge exchange between mathematicians and forestry partners transferable to the management of future pests and diseases.
 
Description Forestry Commission OPM collaboration 
Organisation Forestry Commission
Country United Kingdom 
Sector Public 
PI Contribution As part of the award we have been able to strengthen our relationship with the Forestry Commission, particularly with the OPM Control Programme. Regular meetings and discussions with the OPM Project Manager Andrew Hoppit resulted in my invitation to the OPM stakeholders meeting in January 2023 to network and discuss recent research results. We have also collaborated on a further research publication (https://doi.org/10.3390/d15040496).
Collaborator Contribution Forestry Commission have supported us through regular research meetings and data sharing, culminating in regular invitations to join the annual OPM stakeholders meeting and two research publications (https://doi.org/10.1002/ece3.8871 and https://doi.org/10.3390/d15040496).
Impact OPM Stakeholder Meeting (January 2023) - Engagement Output. Academic publication: https://doi.org/10.3390/d15040496 (March 2023)
Start Year 2022
 
Title An individual-based model with an infectious kernel describing the spatial spread of an epidemic, plus several approximations (both analytic and computational, including an inference scheme) of the reproduction number calculated for different epidemic... 
Description This code allows the user to recreate the figures from the Ecological Modelling paper: "Estimating the reproduction number, R0, from individual-based models of tree disease spread", DOI: XXXX, including model simulations and the inference scheme.Tree populations worldwide are facing an unprecedented threat from a variety of tree diseases and invasive pests. Their spread, exacerbated by increasing globalisation and climate change, has an enormous environmental, economic and social impact. Computational individual-based models are a popular tool for describing and forecasting the spread of tree diseases due to their flexibility and ability to reveal collective behaviours. In this paper we present a versatile individual-based model with a Gaussian infectivity kernel to describe the spread of a generic tree disease through a synthetic treescape. We then explore several methods of calculating the basic reproduction number $R_0$, a characteristic measurement of disease infectivity, defining the expected number of new infections resulting from one newly infected individual throughout their infectious period. It is a useful comparative summary parameter of a disease and can be used to explore the threshold dynamics of epidemics through mathematical models. We demonstrate several methods of estimating $R_0$ through the individual-based model, including contact tracing, inferring the Kermack-McKendrick SIR model parameters using the linear noise approximation, and an analytical approximation. As an illustrative example, we then use the model and each of the methods to calculate estimates of $R_0$ for the ash dieback epidemic in the UK. 
Type Of Technology Software 
Year Produced 2023 
Open Source License? Yes  
URL https://data.ncl.ac.uk/articles/software/An_individual-based_model_with_an_infectious_kernel_describ...
 
Title Quantifying invasive pest dynamics through inference of a two-node epidemic network model 
Description This code allows the user to recreate the figures from the Diversity paper: "Quantifying invasive pest dynamics through inference of a two-node epidemic network model" https://doi.org/10.3390/d15040496, including model simulations and the inference scheme. Invasive woodland pests have substantial ecological, economic, and social impacts, harming biodiversity and ecosystem services. Mathematical modelling informed by Bayesian inference can deepen our understanding of the fundamental behaviours of invasive pests and provide predictive tools for forecasting future spread. A key invasive pest of concern in the UK is the oak processionary moth (OPM). OPM was established in the UK in 2006; it is harmful to both oak trees and humans, and its infestation area is continually expanding. Here, we use a computational inference scheme to estimate the parameters for a two-node network epidemic model to describe the temporal dynamics of OPM in two geographically neighbouring parks (Bushy Park and Richmond Park, London). We show the applicability of such a network model to describing invasive pest dynamics and our results suggest that the infestation within Richmond Park has largely driven the infestation within Bushy Park. 
Type Of Technology Software 
Year Produced 2023 
Open Source License? Yes  
URL https://data.ncl.ac.uk/articles/software/Quantifying_invasive_pest_dynamics_through_inference_of_a_t...
 
Title Software supporting 'Quantifying invasive pest dynamics through inference of a two-node epidemic network model' 
Description This code allows the user to recreate the figures from the Diversity paper: "Quantifying invasive pest dynamics through inference of a two-node epidemic network model" https://doi.org/10.3390/d15040496, including model simulations and the inference scheme. Invasive woodland pests have substantial ecological, economic, and social impacts, harming biodiversity and ecosystem services. Mathematical modelling informed by Bayesian inference can deepen our understanding of the fundamental behaviours of invasive pests and provide predictive tools for forecasting future spread. A key invasive pest of concern in the UK is the oak processionary moth (OPM). OPM was established in the UK in 2006; it is harmful to both oak trees and humans, and its infestation area is continually expanding. Here, we use a computational inference scheme to estimate the parameters for a two-node network epidemic model to describe the temporal dynamics of OPM in two geographically neighbouring parks (Bushy Park and Richmond Park, London). We show the applicability of such a network model to describing invasive pest dynamics and our results suggest that the infestation within Richmond Park has largely driven the infestation within Bushy Park. 
Type Of Technology Software 
Year Produced 2023 
Open Source License? Yes  
URL https://data.ncl.ac.uk/articles/software/Quantifying_invasive_pest_dynamics_through_inference_of_a_t...
 
Title Software supporting 'Quantifying invasive pest dynamics through inference of a two-node epidemic network model' 
Description This code allows the user to recreate the figures from the Diversity paper: "Quantifying invasive pest dynamics through inference of a two-node epidemic network model" https://doi.org/10.3390/d15040496, including model simulations and the inference scheme. Invasive woodland pests have substantial ecological, economic, and social impacts, harming biodiversity and ecosystem services. Mathematical modelling informed by Bayesian inference can deepen our understanding of the fundamental behaviours of invasive pests and provide predictive tools for forecasting future spread. A key invasive pest of concern in the UK is the oak processionary moth (OPM). OPM was established in the UK in 2006; it is harmful to both oak trees and humans, and its infestation area is continually expanding. Here, we use a computational inference scheme to estimate the parameters for a two-node network epidemic model to describe the temporal dynamics of OPM in two geographically neighbouring parks (Bushy Park and Richmond Park, London). We show the applicability of such a network model to describing invasive pest dynamics and our results suggest that the infestation within Richmond Park has largely driven the infestation within Bushy Park. 
Type Of Technology Software 
Year Produced 2023 
Open Source License? Yes  
URL https://data.ncl.ac.uk/articles/software/Quantifying_invasive_pest_dynamics_through_inference_of_a_t...
 
Description Mathematical biology seminar invited speaker (Sheffield) 
Form Of Engagement Activity A talk or presentation
Part Of Official Scheme? No
Geographic Reach National
Primary Audience Other audiences
Results and Impact I was invited to give a talk at the Mathematical Biology and Ecology group at Sheffield University (around 20 in attendance, including applied mathematicians, statisticians and students). This promoted the work ongoing as part of the award.
Year(s) Of Engagement Activity 2022
 
Description OPM Stakeholders Meeting 
Form Of Engagement Activity Participation in an activity, workshop or similar
Part Of Official Scheme? No
Geographic Reach National
Primary Audience Other audiences
Results and Impact My role as PI on this award has allowed me to develop ongoing collaborations with governmental forestry groups, including the Forestry Commission. Through this ongoing relationship I was invited to join the OPM Stakeholders Meeting in London in January 2023, a one day event bringing together policymakers, foresters, land owners and researchers across all sectors to share their findings and experiences of the OPM problem within the UK. This knowledge-sharing activity provided valuable up-to-date information about the current OPM situation which will be used in my ongoing research, as well as a valuable networking opportunity to promote the research of this award and initial new collaborations (e.g., Observatree group).
Year(s) Of Engagement Activity 2023
 
Description Royal Statistical Society Annual Conference 2023 Meeting - Speaker 
Form Of Engagement Activity A talk or presentation
Part Of Official Scheme? No
Geographic Reach International
Primary Audience Professional Practitioners
Results and Impact Presented work on modelling of the oak processionary moth at the RSS Conference in 2023 to an audience of academics and industry statisticians. Post-talk discussions led to new connections with researchers at The Tree Council.
Year(s) Of Engagement Activity 2023
 
Description Southampton Statistical Sciences Research Institute seminar invited speaker 
Form Of Engagement Activity A talk or presentation
Part Of Official Scheme? No
Geographic Reach National
Primary Audience Other audiences
Results and Impact I was invited to give a talk at the Southampton Statistical Sciences Research Institute group at Southampton University (around 50 in attendance, including applied mathematicians, statisticians and students). This promoted the work ongoing as part of the award and contributed towards the generation of new research ideas.
Year(s) Of Engagement Activity 2022