Applying optimal control theory to complex epidemiological models to inform real-world disease management

Infectious Disease in the Workplace: Quantifying Uncertainty in Transmission

Understanding disease transmission in the workplace is essential for protecting workers. To model disease outbreaks, the small populations in many workplaces require that stochastic effects are considered, which results in higher uncertainty. The aim of this study was to quantify and interpret the uncertainty inherent in such circumstances. We assessed how uncertainty of an outbreak in workplaces depends on i) the infection dynamics in the community, ii) the workforce size, iii) spatial structure in the workplace, iv) heterogeneity in susceptibility of workers, and v) heterogeneity in infectiousness of workers. To address these questions, we developed a multiscale model: A deterministic model to predict community transmission, and a stochastic model to predict workplace transmission. We extended this basic workplace model to allow for spatial structure, and heterogeneity in susceptibility and infectiousness in workers. We found a non-monotonic relationship between the workplace transmission rate and the coefficient of variation (CV), which we use as a measure of uncertainty. Increasing community transmission, workforce size and heterogeneity in susceptibility decreased the CV. Conversely, increasing the level of spatial structure and heterogeneity in infectiousness increased the CV. However, when the model predicts bimodal distributions, for example when community transmission is low and workplace transmission is high, the CV fails to capture this uncertainty. Overall, our work informs modellers and policy makers on how model complexity impacts outbreak uncertainty. In particular: workforce size, community and workplace transmission, spatial structure and individual heterogeneity contribute in a specific and individual manner to the predicted workplace outbreak size distribution.

Emerging Themes and Approaches in Plant Virus Epidemiology

Plant diseases caused by viruses share many common features with those caused by other pathogen taxa in terms of the host-pathogen interaction, but there are also distinctive features in epidemiology, most apparent where transmission is by vectors. Consequently, the host-virus-vector-environment interaction presents a continuing challenge in attempts to understand and predict the course of plant virus epidemics. Theoretical concepts, based on the underlying biology, can be expressed in mathematical models and tested through quantitative assessments of epidemics in the field; this remains a goal in understanding why plant virus epidemics occur and how they can be controlled. To this end, this review identifies recent emerging themes and approaches to fill in knowledge gaps in plant virus epidemiology. We review quantitative work on the impact of climatic fluctuations and change on plants, viruses, and vectors under different scenarios where impacts on the individual components of the plant-virus-vector interaction may vary disproportionately; there is a continuing, sometimes discordant, debate on host resistance and tolerance as plant defense mechanisms, including aspects of farmer behavior and attitudes toward disease management that may affect deployment in crops; disentangling host-virus-vector-environment interactions, as these contribute to temporal and spatial disease progress in field populations; computational techniques for estimating epidemiological parameters from field observations; and the use of optimal control analysis to assess disease control options. We end by proposing new challenges and questions in plant virus epidemiology.

Modeling plant diseases under climate change: evolutionary perspectives

Infectious plant diseases are a major threat to global agricultural productivity, economic development, and ecological integrity. There is widespread concern that these social and natural disasters caused by infectious plant diseases may escalate with climate change and computer modeling offers a unique opportunity to address this concern. Here, we analyze the intrinsic problems associated with current modeling strategies and highlight the need to integrate evolutionary principles into polytrophic, eco-evolutionary frameworks to improve predictions. We particularly discuss how evolutionary shifts in functional trade-offs, relative adaptability between plants and pathogens, ecosystems, and climate preferences induced by climate change may feedback to future plant disease epidemics and how technological advances can facilitate the generation and integration of this relevant knowledge for better modeling predictions.

Optimal vaccine allocation for the control of sexually transmitted infections

The burden of sexually transmitted infections (STIs) poses a challenge due to its large negative impact on sexual and reproductive health worldwide. Besides simple prevention measures and available treatment efforts, prophylactic vaccination is a powerful tool for controlling some viral STIs and their associated diseases. Here, we investigate how prophylactic vaccines are best distributed to prevent and control STIs. We consider sex-specific differences in susceptibility to infection, as well as disease severity outcomes. Different vaccination strategies are compared assuming distinct budget constraints that mimic a scarce vaccine stockpile. Vaccination strategies are obtained as solutions to an optimal control problem subject to a two-sex Kermack-McKendrick-type model, where the control variables are the daily vaccination rates for females and males. One important aspect of our approach relies on conceptualizing a limited but specific vaccine stockpile via an isoperimetric constraint. We solve the optimal control problem via Pontryagin's Maximum Principle and obtain a numerical approximation for the solution using a modified version of the forward-backward sweep method that handles the isoperimetric budget constraint in our formulation. The results suggest that for a limited vaccine supply ([Formula: see text]-[Formula: see text] vaccination coverage), one-sex vaccination, prioritizing females, appears to be more beneficial than the inclusion of both sexes into the vaccination program. Whereas, if the vaccine supply is relatively large (enough to reach at least [Formula: see text] coverage), vaccinating both sexes, with a slightly higher rate for females, is optimal and provides an effective and faster approach to reducing the prevalence of the infection.

Using 'sentinel' plants to improve early detection of invasive plant pathogens

Infectious diseases of plants present an ongoing and increasing threat to international biosecurity, with wide-ranging implications. An important challenge in plant disease management is achieving early detection of invading pathogens, which requires effective surveillance through the implementation of appropriate monitoring programmes. However, when monitoring relies on visual inspection as a means of detection, surveillance is often hindered by a long incubation period (delay from infection to symptom onset) during which plants may be infectious but not displaying visible symptoms. 'Sentinel' plants-alternative susceptible host species that display visible symptoms of infection more rapidly-could be introduced to at-risk populations and included in monitoring programmes to act as early warning beacons for infection. However, while sentinel hosts exhibit faster disease progression and so allow pathogens to be detected earlier, this often comes at a cost: faster disease progression typically promotes earlier onward transmission. Here, we construct a computational model of pathogen transmission to explore this trade-off and investigate how including sentinel plants in monitoring programmes could facilitate earlier detection of invasive plant pathogens. Using Xylella fastidiosa infection in Olea europaea (European olive) as a current high profile case study, for which Catharanthus roseus (Madagascan periwinkle) is a candidate sentinel host, we apply a Bayesian optimisation algorithm to determine the optimal number of sentinel hosts to introduce for a given sampling effort, as well as the optimal division of limited surveillance resources between crop and sentinel plants. Our results demonstrate that including sentinel plants in monitoring programmes can reduce the expected prevalence of infection upon outbreak detection substantially, increasing the feasibility of local outbreak containment.

Airline flight delays using artificial intelligence in COVID-19 with perspective analytics

This study envisages assessing the effects of the COVID-19 on the on-time performance of US-airlines industry in the disrupted situations. The deep learning techniques used are neural network regression, decision forest regression, boosted decision tree regression and multi class logistic regression. The best technique is identified. In the perspective data analytics, it is suggested what the airlines should do for the on-time performance in the disrupted situation. The performances of all the methods are satisfactory. The coefficient of determination for the neural network regression is 0.86 and for decision forest regression is 0.85, respectively. The coefficient of determination for the boosted decision tree is 0.870984. Thus boosted decision tree regression is better. Multi class logistic regression gives an overall accuracy and precision of 98.4%. Recalling/remembering performance is 99%. Thus multi class logistic regression is the best model for prediction of flight delays in the COVID-19. The confusion matrix for the multi class logistic regression shows that 87.2% flights actually not delayed are predicted not delayed. The flights actually not delayed but wrongly predicted delayed are12.7%. The strength of relation with departure delay, carrier delay, late aircraft delay, weather delay and NAS delay, are 94%, 53%, 35%, 21%, and 14%, respectively. There is a weak negative relation (almost unrelated) with the air time and arrival delay. Security delay and arrival delay are also almost unrelated with strength of 1% relationship. Based on these diagnostic analytics, it is recommended as perspective to take due care reducing departure delay, carrier delay, Late aircraft delay, weather delay and Nas delay, respectively, considerably with effect of 94%, 53%, 35%, 21%, and 14% in disrupted situations. The proposed models have MAE of 2% for Neural Network Regression, Decision Forest Regression, Boosted Decision Tree Regression, respectively, and, RMSE approximately, 11%, 12%, 11%, respectively.

Optimal Resistance Management for Mixtures of High-Risk Fungicides: Robustness to the Initial Frequency of Resistance and Pathogen Sexual Reproduction

There is a strong consensus that selection for fungicide resistant pathogen strains can be most effectively limited by using applications of mixtures of fungicides designed to balance disease control against selection. However, how to do this in practice is not entirely characterized. Previous work indicates optimal mixtures of pairs of fungicides which are both at a high risk of resistance can be constructed using pairs of doses that select equally for both single resistant strains in the first year of application. What has not been addressed thus far is the important real-world case in which the initial levels of resistance to each fungicide differ, for example because the chemicals have been available for different lengths of time. We show how recommendations based on equal selection in the first year can be suboptimal in this case. We introduce a simple alternative approach, based on equalizing the frequencies of single resistant strains in the year that achieving acceptable levels of control is predicted to become impossible. We show that this strategy is robust to changes in parameters controlling pathogen epidemiology and fungicide efficacy. We develop our recommendation using a preexisting, parameterized model of Zymoseptoria tritici (the pathogen causing Septoria leaf blotch on wheat), which exemplifies the range of plant pathogens that predominantly spread clonally, but for which sexual reproduction forms an important component of the life cycle. We show that pathogen sexual reproduction can influence the rate at which fungicide resistance develops but does not qualitatively affect our optimal resistance management recommendation. [Formula: see text] Copyright © 2023 The Author(s). This is an open access article distributed under the CC BY 4.0 International license.

Optimal reduced-mixing for an SIS infectious-disease model

Which reduced-mixing strategy maximizes economic output during a disease outbreak? To answer this question, we formulate an optimal-control problem that maximizes the difference between revenue, due to healthy individuals, and medical costs, associated with infective individuals, for SIS disease dynamics. The control variable is the level of mixing in the population, which influences both revenue and the spread of the disease. Using Pontryagin's maximum principle, we find a closed-form solution for our problem. We explore an example of our problem with parameters for the transmission of Staphylococcus aureus in dairy cows, and we perform sensitivity analyses to determine how model parameters affect optimal strategies. We find that less mixing is preferable when the transmission rate is high, the per-capita recovery rate is low, or when the revenue parameter is much smaller than the cost parameter.

Long-term spatial and population-structured planning of non-pharmaceutical interventions to epidemic outbreaks

In this paper, we consider the problem of planning non-pharmaceutical interventions to control the spread of infectious diseases. We propose a new model derived from classical compartmental models; however, we model spatial and population-structure heterogeneity of population mixing. The resulting model is a large-scale non-linear and non-convex optimisation problem. In order to solve it, we apply a special variant of covariance matrix adaptation evolution strategy. We show that results obtained for three different objectives are better than natural heuristics and, moreover, that the introduction of an individual's mobility to the model is significant for the quality of the decisions. We apply our approach to a six-compartmental model with detailed Poland and COVID-19 disease data. The obtained results are non-trivialand sometimes unexpected; therefore, we believe that our model could be applied to support policy-makers in fighting diseases at the long-term decision-making level.

Dynamic resource allocation for controlling pathogen spread on a large metapopulation network

To control the spread of an infectious disease over a large network, the optimal allocation by a social planner of a limited resource is a fundamental and difficult problem. We address this problem for a livestock disease that propagates on an animal trade network according to an epidemiological–demographic model based on animal demographics and trade data. We assume that the resource is dynamically allocated following a certain score, up to the limit of resource availability. We adapt a greedy approach to the metapopulation framework, obtaining new scores that minimize approximations of two different objective functions, for two control measures: vaccination and treatment. Through intensive simulations, we compare the greedy scores with several heuristics. Although topology-based scores can limit the spread of the disease, information on herd health status seems crucial to eradicating the disease. In particular, greedy scores are among the most effective in reducing disease prevalence, even though they do not always perform the best. However, some scores may be preferred in real life because they are easier to calculate or because they use a smaller amount of resources. The developed approach could be adapted to other epidemiological models or to other control measures in the metapopulation setting.

Optimal strategies to protect a sub-population at risk due to an established epidemic

Epidemics can particularly threaten certain sub-populations. For example, for severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), the elderly are often preferentially protected. For diseases of plants and animals, certain sub-populations can drive mitigation because they are intrinsically more valuable for ecological, economic, socio-cultural or political reasons. Here, we use optimal control theory to identify strategies to optimally protect a ‘high-value’ sub-population when there is a limited budget and epidemiological uncertainty. We use protection of the Redwood National Park in California in the face of the large ongoing state-wide epidemic of sudden oak death (caused by Phytophthora ramorum) as a case study. We concentrate on whether control should be focused entirely within the National Park itself, or whether treatment of the growing epidemic in the surrounding ‘buffer region’ can instead be more profitable. We find that, depending on rates of infection and the size of the ongoing epidemic, focusing control on the high-value region is often optimal. However, priority should sometimes switch from the buffer region to the high-value region only as the local outbreak grows. We characterize how the timing of any switch depends on epidemiological and logistic parameters, and test robustness to systematic misspecification of these factors due to imperfect prior knowledge.

From data to complex network control of airline flight delays

Many critical complex systems and networks are continuously monitored, creating vast volumes of data describing their dynamics. To understand and optimize their performance, we need to discover and formalize their dynamics to enable their control. Here, we introduce a multidisciplinary framework using network science and control theory to accomplish these goals. We demonstrate its use on a meaningful example of a complex network of U.S. domestic passenger airlines aiming to control flight delays. Using the real data on such delays, we build a flight delay network for each airline. Analyzing these networks, we uncover and formalize their dynamics. We use this formalization to design the optimal control for the flight delay networks. The results of applying this control to the ground truth data on flight delays demonstrate the low costs of the optimal control and significant reduction of delay times, while the costs of the delays unabated by control are high. Thus, the introduced here framework benefits the passengers, the airline companies and the airports.

Data-driven methods for present and future pandemics: Monitoring, modelling and managing

This survey analyses the role of data-driven methodologies for pandemic modelling and control. We provide a roadmap from the access to epidemiological data sources to the control of epidemic phenomena. We review the available methodologies and discuss the challenges in the development of data-driven strategies to combat the spreading of infectious diseases. Our aim is to bring together several different disciplines required to provide a holistic approach to epidemic analysis, such as data science, epidemiology, and systems-and-control theory. A 3M-analysis is presented, whose three pillars are: Monitoring, Modelling and Managing. The focus is on the potential of data-driven schemes to address three different challenges raised by a pandemic: (i) monitoring the epidemic evolution and assessing the effectiveness of the adopted countermeasures; (ii) modelling and forecasting the spread of the epidemic; (iii) making timely decisions to manage, mitigate and suppress the contagion. For each step of this roadmap, we review consolidated theoretical approaches (including data-driven methodologies that have been shown to be successful in other contexts) and discuss their application to past or present epidemics, such as Covid-19, as well as their potential application to future epidemics.

State estimation-based control of COVID-19 epidemic before and after vaccine development

In this study, a nonlinear robust control policy is designed together with a state observer in order to manage the novel coronavirus disease (COVID-19) outbreak having an uncertain epidemiological model with unmeasurable variables. This nonlinear model for the COVID-19 epidemic includes eight state variables (susceptible, exposed, infected, quarantined, hospitalized, recovered, deceased, and insusceptible populations). Two plausible scenarios are put forward in this article to control this epidemic before and after its vaccine invention. In the first scenario, the social distancing and hospitalization rates are employed as two applicable control inputs to diminish the exposed and infected groups. However, in the second scenario after the vaccine development, the vaccination rate is taken into account as the third control input to reduce the susceptible populations, in addition to the two objectives of the first scenario. The proposed feedback control measures are defined in terms of the hospitalized and deceased populations due to the available statistical data, while other unmeasurable compartmental variables are estimated by an extended Kalman filter (EKF). In other words, the susceptible, exposed, infected, quarantined, recovered, and insusceptible individuals cannot be identified precisely because of the asymptomatic infection of COVID-19 in some cases, its incubation period, and the lack of an adequate community screening. Utilizing the Lyapunov theorem, the stability and bounded tracking convergence of the closed-loop epidemiological system are investigated in the presence of modeling uncertainties. Finally, a comprehensive simulation study is conducted based on Canada's reported cases for two defined timing plans (with different treatment rates). Obtained results demonstrate that the developed EKF-based control scheme can achieve desired epidemic goals (exponential decrease of infected, exposed, and susceptible people).

Optimal Control of Plant Disease Epidemics with Clean Seed Usage

The distribution and use of pathogen-free planting material ("clean seeds") is a promising method to control plant diseases in developing countries. We address the question of minimizing disease prevalence in plants through the optimal usage of clean seeds. We consider the simplest possible S-I model together with a simple economic criterion to be maximized. The static optimization problem shows a diversity of possible outcomes depending on economical and epidemiological parameters. We derive a simple condition showing to what extent subsidizing clean seeds relative to the epidemiological features of the disease may help eradicate or control the disease. Then we consider dynamic optimal control and Pontryagin's maximum principle to study the optimal usage of clean seeds to control the disease. The dynamical results are comparable to the static ones and are even simpler in some sense. In particular, the condition on the critical subsidy rate that makes clean seed usage economically viable is unchanged from the static optimization case. We discuss how these results may apply to the control of maize lethal necrosis in East-Africa.

The DIOS framework for optimizing infectious disease surveillance: Numerical methods for simulation and multi-objective optimization of surveillance network architectures

Infectious disease surveillance systems provide vital data for guiding disease prevention and control policies, yet the formalization of methods to optimize surveillance networks has largely been overlooked. Decisions surrounding surveillance design parameters-such as the number and placement of surveillance sites, target populations, and case definitions-are often determined by expert opinion or deference to operational considerations, without formal analysis of the influence of design parameters on surveillance objectives. Here we propose a simulation framework to guide evidence-based surveillance network design to better achieve specific surveillance goals with limited resources. We define evidence-based surveillance design as an optimization problem, acknowledging the many operational constraints under which surveillance systems operate, the many dimensions of surveillance system design, the multiple and competing goals of surveillance, and the complex and dynamic nature of disease systems. We describe an analytical framework-the Disease Surveillance Informatics Optimization and Simulation (DIOS) framework-for the identification of optimal surveillance designs through mathematical representations of disease and surveillance processes, definition of objective functions, and numerical optimization. We then apply the framework to the problem of selecting candidate sites to expand an existing surveillance network under alternative objectives of: (1) improving spatial prediction of disease prevalence at unmonitored sites; or (2) estimating the observed effect of a risk factor on disease. Results of this demonstration illustrate how optimal designs are sensitive to both surveillance goals and the underlying spatial pattern of the target disease. The findings affirm the value of designing surveillance systems through quantitative and adaptive analysis of network characteristics and performance. The framework can be applied to the design of surveillance systems tailored to setting-specific disease transmission dynamics and surveillance needs, and can yield improved understanding of tradeoffs between network architectures.

Model predictive control to mitigate the COVID-19 outbreak in a multi-region scenario

The COVID-19 outbreak is deeply influencing the global social and economic framework, due to restrictive measures adopted worldwide by governments to counteract the pandemic contagion. In multi-region areas such as Italy, where the contagion peak has been reached, it is crucial to find targeted and coordinated optimal exit and restarting strategies on a regional basis to effectively cope with possible onset of further epidemic waves, while efficiently returning the economic activities to their standard level of intensity. Differently from the related literature, where modeling and controlling the pandemic contagion is typically addressed on a national basis, this paper proposes an optimal control approach that supports governments in defining the most effective strategies to be adopted during post-lockdown mitigation phases in a multi-region scenario. Based on the joint use of a non-linear Model Predictive Control scheme and a modified Susceptible-Infected-Recovered (SIR)-based epidemiological model, the approach is aimed at minimizing the cost of the so-called non-pharmaceutical interventions (that is, mitigation strategies), while ensuring that the capacity of the network of regional healthcare systems is not violated. In addition, the proposed approach supports policy makers in taking targeted intervention decisions on different regions by an integrated and structured model, thus both respecting the specific regional health systems characteristics and improving the system-wide performance by avoiding uncoordinated actions of the regions. The methodology is tested on the COVID-19 outbreak data related to the network of Italian regions, showing its effectiveness in properly supporting the definition of effective regional strategies for managing the COVID-19 diffusion.

Bayesian optimisation of restriction zones for bluetongue control

We investigate the restriction of animal movements as a method to control the spread of bluetongue, an infectious disease of livestock that is becoming increasingly prevalent due to the onset of climate change. We derive control policies for the UK that minimise the number of infected farms during an outbreak using Bayesian optimisation and a simulation-based model of BT. Two cases are presented: first, where the region of introduction is randomly selected from England and Wales to find a generalised strategy. This “national” model is shown to be just as effective at subduing the spread of bluetongue as the current strategy of the UK government. Our proposed controls are simpler to implement, affect fewer farms in the process and, in so doing, minimise the potential economic implications. Second, we consider policies that are tailored to the specific region in which the first infection was detected. Seven different regions in the UK were explored and improvements in efficiency from the use of specialised policies presented. As a consequence of the increasing temperatures associated with climate change, efficient control measures for vector-borne diseases such as this are expected to become increasingly important. Our work demonstrates the potential value of using Bayesian optimisation in developing cost-effective disease management strategies.

Epidemiologically and Socio-economically Optimal Policies via Bayesian Optimization

Mass public quarantining, colloquially known as a lock-down, is a non-pharmaceutical intervention to check spread of disease. This paper presents ESOP (Epidemiologically and Socio-economically Optimal Policies), a novel application of active machine learning techniques using Bayesian optimization, that interacts with an epidemiological model to arrive at lock-down schedules that optimally balance public health benefits and socio-economic downsides of reduced economic activity during lock-down periods. The utility of ESOP is demonstrated using case studies with VIPER (Virus-Individual-Policy-EnviRonment), a stochastic agent-based simulator that this paper also proposes. However, ESOP is flexible enough to interact with arbitrary epidemiological simulators in a black-box manner, and produce schedules that involve multiple phases of lock-downs.

A theoretical framework for transitioning from patient-level to population-scale epidemiological dynamics: influenza A as a case study

Multi-scale epidemic forecasting models have been used to inform population-scale predictions with within-host models and/or infection data collected in longitudinal cohort studies. However, most multi-scale models are complex and require significant modelling expertise to run. We formulate an alternative multi-scale modelling framework using a compartmental model with multiple infected stages. In the large-compartment limit, our easy-to-use framework generates identical results compared to previous more complicated approaches. We apply our framework to the case study of influenza A in humans. By using a viral dynamics model to generate synthetic patient-level data, we explore the effects of limited and inaccurate patient data on the accuracy of population-scale forecasts. If infection data are collected daily, we find that a cohort of at least 40 patients is required for a mean population-scale forecasting error below 10%. Forecasting errors may be reduced by including more patients in future cohort studies or by increasing the frequency of observations for each patient. Our work, therefore, provides not only an accessible epidemiological modelling framework but also an insight into the data required for accurate forecasting using multi-scale models.

Preface to theme issue 'Modelling infectious disease outbreaks in humans, animals and plants: epidemic forecasting and control'

This preface forms part of the theme issue 'Modelling infectious disease outbreaks in humans, animals and plants: epidemic forecasting and control'. This theme issue is linked with the earlier issue 'Modelling infectious disease outbreaks in humans, animals and plants: approaches and important themes'.

Detection, forecasting and control of infectious disease epidemics: modelling outbreaks in humans, animals and plants

The 1918 influenza pandemic is one of the most devastating infectious disease epidemics on record, having caused approximately 50 million deaths worldwide. Control measures, including prohibiting non-essential gatherings as well as closing cinemas and music halls, were applied with varying success and limited knowledge of transmission dynamics. One hundred years later, following developments in the field of mathematical epidemiology, models are increasingly used to guide decision-making and devise appropriate interventions that mitigate the impacts of epidemics. Epidemiological models have been used as decision-making tools during outbreaks in human, animal and plant populations. However, as the subject has developed, human, animal and plant disease modelling have diverged. Approaches have been developed independently for pathogens of each host type, often despite similarities between the models used in these complementary fields. With the increased importance of a One Health approach that unifies human, animal and plant health, we argue that more inter-disciplinary collaboration would enhance each of the related disciplines. This pair of theme issues presents research articles written by human, animal and plant disease modellers. In this introductory article, we compare the questions pertinent to, and approaches used by, epidemiological modellers of human, animal and plant pathogens, and summarize the articles in these theme issues. We encourage future collaboration that transcends disciplinary boundaries and links the closely related areas of human, animal and plant disease epidemic modelling. This article is part of the theme issue 'Modelling infectious disease outbreaks in humans, animals and plants: approaches and important themes'. This issue is linked with the subsequent theme issue 'Modelling infectious disease outbreaks in humans, animals and plants: epidemic forecasting and control'.

Applying optimal control theory to a spatial simulation model of sudden oak death: ongoing surveillance protects tanoak while conserving biodiversity

Sudden oak death has devastated tree populations across California. However, management might still slow disease spread at local scales. We demonstrate how to unambiguously characterize effective, local management strategies using a detailed, spatially explicit simulation model of spread in a single forest stand. This pre-existing, parameterized simulation is approximated here by a carefully calibrated, non-spatial model, explicitly constructed to be sufficiently simple to allow optimal control theory (OCT) to be applied. By lifting management strategies from the approximate model to the detailed simulation, effective time-dependent controls can be identified. These protect tanoak-a culturally and ecologically important species-while conserving forest biodiversity within a limited budget. We also consider model predictive control, in which both the approximating model and optimal control are repeatedly updated as the epidemic progresses. This allows management which is robust to both parameter uncertainty and systematic differences between simulation and approximate models. Including the costs of disease surveillance then introduces an optimal intensity of surveillance. Our study demonstrates that successful control of sudden oak death is likely to rely on adaptive strategies updated via ongoing surveillance. More broadly, it illustrates how OCT can inform effective real-world management, even when underpinning disease spread models are highly complex.