COVID-19 modeling based on real geographic and population data

BACKGROUND

: Intercity travel is one of the most important parameters for combating a pandemic. The ongoing COVID-19 pandemic has resulted in different computational studies involving intercity connections. In this study, the effects of intercity connections during an epidemic such as COVID-19 are evaluated using a new network model.

METHODS

This model considers the actual geographic neighborhood and population density data. This new model is applied to actual Turkish data by means of provincial connections and populations. A Monte Carlo algorithm with a hybrid lattice model is applied to a lattice with 8802 data points.

RESULTS

Around Monte Carlo step 70, the number of active cases in Türkiye reaches up to 8.0% of the total population, which is followed by a second wave at around Monte Carlo step 100. The number of active cases vanishes around Monte Carlo step 160. Starting with İstanbul, the epidemic quickly expands between steps 60 and 100. Simulation results fit the actual mortality data in Türkiye.

DISCUSSION

This model is quantitatively very efficient in modeling real-world COVID-19 epidemic data based on populations and geographical intercity connections, by means of estimating the number of deaths, disease spread, and epidemic termination.

A new logistic growth model applied to COVID-19 fatality data

BACKGROUND

Recent work showed that the temporal growth of the novel coronavirus disease (COVID-19) follows a sub-exponential power-law scaling whenever effective control interventions are in place. Taking this into consideration, we present a new phenomenological logistic model that is well-suited for such power-law epidemic growth.

METHODS

We empirically develop the logistic growth model using simple scaling arguments, known boundary conditions and a comparison with available data from four countries, Belgium, China, Denmark and Germany, where (arguably) effective containment measures were put in place during the first wave of the pandemic. A non-linear least-squares minimization algorithm is used to map the parameter space and make optimal predictions.

RESULTS

Unlike other logistic growth models, our presented model is shown to consistently make accurate predictions of peak heights, peak locations and cumulative saturation values for incomplete epidemic growth curves. We further show that the power-law growth model also works reasonably well when containment and lock down strategies are not as stringent as they were during the first wave of infections in 2020. On the basis of this agreement, the model was used to forecast COVID-19 fatalities for the third wave in South Africa, which was in progress during the time of this work.

CONCLUSION

We anticipate that our presented model will be useful for a similar forecasting of COVID-19 induced infections/deaths in other regions as well as other cases of infectious disease outbreaks, particularly when power-law scaling is observed.

A random walk Monte Carlo simulation study of COVID-19-like infection spread

Recent analysis of early COVID-19 data from China showed that the number of confirmed cases followed a subexponential power-law increase, with a growth exponent of around 2.2 (Maier and Brockmann, 2020). The power-law behavior was attributed to a combination of effective containment and mitigation measures employed as well as behavioral changes by the population. In this work, we report a random walk Monte Carlo simulation study of proximity-based infection spread. Control interventions such as lockdown measures and mobility restrictions are incorporated in the simulations through a single parameter, the size of each step in the random walk process. The step size l is taken to be a multiple of 〈 r 〉 , which is the average separation between individuals. Three temporal growth regimes (quadratic, intermediate power-law and exponential) are shown to emerge naturally from our simulations. For l = 〈 r 〉 , we get intermediate power-law growth exponents that are in general agreement with available data from China. On the other hand, we obtain a quadratic growth for smaller step sizes l ≲ 〈 r 〉 ∕ 2 , while for large l the growth is found to be exponential. We further performed a comparative case study of early fatality data (under varying levels of lockdown conditions) from three other countries, India, Brazil and South Africa. We show that reasonable agreement with these data can be obtained by incorporating small-world-like connections in our simulations.

Equilibrium properties of the spatial SIS model as a point pattern dynamics - How is infection distributed over space?

We revisit the classical epidemiological SIS model as a stochastic point pattern dynamics with special focus on its spatial distribution at equilibrium. In this model, each point on a continuous space is either susceptible S or infectious I, and infection occurs with an infection kernel as a function of distance from I to S. This stochastic process has been mathematically described by the hierarchical dynamics of the probabilities that a point, a pair made by two points, and a triplet made by three points, etc., is in a specific configuration of status. Using a simple closure thereby triplet probabilities that appear in the dynamics are approximated, we show that the average singlet probabilities and the pair probabilities that describe spatial distributions of Ss and Is at equilibrium can be explicitly derived using the infection kernel; Is are spatially clustered in the same order of the infection kernel. The results highlight the advantage of point pattern approach to model spatial population dynamics in general ecology where local interactions among individuals likely depend on distance between them.

Dynamics analysis of SIR epidemic model with correlation coefficients and clustering coefficient in networks

In this paper, the correlation coefficients between nodes in states are used as dynamic variables, and we construct SIR epidemic dynamic models with correlation coefficients by using the pair approximation method in static networks and dynamic networks, respectively. Considering the clustering coefficient of the network, we analytically investigate the existence and the local asymptotic stability of each equilibrium of these models and derive threshold values for the prevalence of diseases. Additionally, we obtain two equivalent epidemic thresholds in dynamic networks, which are compared with the results of the mean field equations.

Estimating the delay between host infection and disease (incubation period) and assessing its significance to the epidemiology of plant diseases

Knowledge of the incubation period of infectious diseases (time between host infection and expression of disease symptoms) is crucial to our epidemiological understanding and the design of appropriate prevention and control policies. Plant diseases cause substantial damage to agricultural and arboricultural systems, but there is still very little information about how the incubation period varies within host populations. In this paper, we focus on the incubation period of soilborne plant pathogens, which are difficult to detect as they spread and infect the hosts underground and above-ground symptoms occur considerably later. We conducted experiments on Rhizoctonia solani in sugar beet, as an example patho-system, and used modelling approaches to estimate the incubation period distribution and demonstrate the impact of differing estimations on our epidemiological understanding of plant diseases. We present measurements of the incubation period obtained in field conditions, fit alternative probability models to the data, and show that the incubation period distribution changes with host age. By simulating spatially-explicit epidemiological models with different incubation-period distributions, we study the conditions for a significant time lag between epidemics of cryptic infection and the associated epidemics of symptomatic disease. We examine the sensitivity of this lag to differing distributional assumptions about the incubation period (i.e. exponential versus Gamma). We demonstrate that accurate information about the incubation period distribution of a pathosystem can be critical in assessing the true scale of pathogen invasion behind early disease symptoms in the field; likewise, it can be central to model-based prediction of epidemic risk and evaluation of disease management strategies. Our results highlight that reliance on observation of disease symptoms can cause significant delay in detection of soil-borne pathogen epidemics and mislead practitioners and epidemiologists about the timing, extent, and viability of disease control measures for limiting economic loss.

Endemic infections are always possible on regular networks

We study the dependence of the largest component in regular networks on the clustering coefficient, showing that its size changes smoothly without undergoing a phase transition. We explain this behavior via an analytical approach based on the network structure, and provide an exact equation describing the numerical results. Our work indicates that intrinsic structural properties always allow the spread of epidemics on regular networks.

Simplified method for including spatial correlations in mean-field approximations

Biological systems involving proliferation, migration, and death are observed across all scales. For example, they govern cellular processes such as wound healing, as well as the population dynamics of groups of organisms. In this paper, we provide a simplified method for correcting mean-field approximations of volume-excluding birth-death-movement processes on a regular lattice. An initially uniform distribution of agents on the lattice may give rise to spatial heterogeneity, depending on the relative rates of proliferation, migration, and death. Many frameworks chosen to model these systems neglect spatial correlations, which can lead to inaccurate predictions of their behavior. For example, the logistic model is frequently chosen, which is the mean-field approximation in this case. This mean-field description can be corrected by including a system of ordinary differential equations for pairwise correlations between lattice site occupancies at various lattice distances. In this work we discuss difficulties with this method and provide a simplification in the form of a partial differential equation description for the evolution of pairwise spatial correlations over time. We test our simplified model against the more complex corrected mean-field model, finding excellent agreement. We show how our model successfully predicts system behavior in regions where the mean-field approximation shows large discrepancies. Additionally, we investigate regions of parameter space where migration is reduced relative to proliferation, which has not been examined in detail before and find our method is successful at correcting the deviations observed in the mean-field model in these parameter regimes.

A MATHEMATICAL MODEL FOR SPATIALLY EXPANDING INFECTED AREA OF EPIDEMICS TRANSMITTED THROUGH HETEROGENEOUSLY DISTRIBUTED SUSCEPTIBLE UNITS

Little is known about the effect of environmental heterogeneity on the spatial expansion of epidemics. In this work, to focus on the question of how the extent of epidemic damage depends on the spatial distribution of susceptible units, we develop a mathematical model with a simple stochastic process, and analyze it. We assume that the unit of infection is immobile, as town, plant, etc. and classify the units into three classes: susceptible, infective and recovered. We consider the range expanded by infected units, the infected rangeR, assuming a certain generalized relation between R and the total number of infected units k, making use of an index, a sort of fractal dimension, to characterize the spatial distribution of infected units. From the results of our modeling analysis, we show that the expected velocity of spatial expansion of infected range is significantly affected by the fractal nature of spatial distribution of immobile susceptible units, and is temporally variable. When the infection finally terminates at a moment, the infected range at the moment is closely related to the nature of spatial distribution of immobile susceptible units, which is explicitly demonstrated in our analysis.

Critical parameters for modelling the spread of foot-and-mouth disease in wildlife

A series of simulation experiments was conducted to determine how estimates of the latent and infectious periods, number of neighbours (contacts) and population size impact on the predicted magnitude and distribution of foot-and-mouth disease (FMD) outbreaks in white-tailed deer in southern Texas. Outbreaks were simulated using a previously developed and applied susceptible–latent–infected–recovered geographic automata model. There were substantial differences in the estimated predicted number of deer and locations infected, based on the model parameters used (3779–119 879 deer infected and 227–6526 locations affected). There were also substantial differences in the spatial risk of infection based on the model parameters used. The predicted spread of FMD was found to be most sensitive to the assumed latent period and the assumed number of contacts. How these parameters are estimated is likely to be critical in studies on the impact of FMD spread in situations in which wildlife reservoirs might potentially exist.

Effects of host social hierarchy on disease persistence

The effects of social hierarchy on population dynamics and epidemiology are examined through a model which contains a number of fundamental features of hierarchical systems, but is simple enough to allow analytical insight. In order to allow for differences in birth rates, contact rates and movement rates among different sets of individuals the population is first divided into subgroups representing levels in the hierarchy. Movement, representing dominance challenges, is allowed between any two levels, giving a completely connected network. The model includes hierarchical effects by introducing a set of dominance parameters which affect birth rates in each social level and movement rates between social levels, dependent upon their rank. Although natural hierarchies vary greatly in form, the skewing of contact patterns, introduced here through non-uniform dominance parameters, has marked effects on the spread of disease. A simple homogeneous mixing differential equation model of a disease with SI dynamics in a population subject to simple birth and death process is presented and it is shown that the hierarchical model tends to this as certain parameter regions are approached. Outside of these parameter regions correlations within the system give rise to deviations from the simple theory. A Gaussian moment closure scheme is developed which extends the homogeneous model in order to take account of correlations arising from the hierarchical structure, and it is shown that the results are in reasonable agreement with simulations across a range of parameters. This approach helps to elucidate the origin of hierarchical effects and shows that it may be straightforward to relate the correlations in the model to measurable quantities which could be used to determine the importance of hierarchical corrections. Overall, hierarchical effects decrease the levels of disease present in a given population compared to a homogeneous unstructured model, but show higher levels of disease than structured models with no hierarchy. The separation between these three models is greatest when the rate of dominance challenges is low, reducing mixing, and when the disease prevalence is low. This suggests that these effects will often need to be considered in models being used to examine the impact of control strategies where the low disease prevalence behaviour of a model is critical.

Representation of animal distributions in space: how geostatistical estimates impact simulation modeling of foot-and-mouth disease spread

Modeling potential disease spread in wildlife populations is important for predicting, responding to and recovering from a foreign animal disease incursion. To make spatial epidemic predictions, the target animal species of interest must first be represented in space. We conducted a series of simulation experiments to determine how estimates of the spatial distribution of white-tailed deer impact the predicted magnitude and distribution of foot-and-mouth disease (FMD) outbreaks. Outbreaks were simulated using a susceptible-infected-recovered geographic automata model. The study region was a 9-county area (24 000 km(2)) of southern Texas. Methods used for creating deer distributions included dasymetric mapping, kriging and remotely sensed image analysis. The magnitudes and distributions of the predicted outbreaks were evaluated by comparing the median number of deer infected and median area affected (km(2)), respectively. The methods were further evaluated for similar predictive power by comparing the model predicted outputs with unweighted pair group method with arithmetic mean (UPGMA) clustering. There were significant differences in the estimated number of deer in the study region, based on the geostatistical estimation procedure used (range: 385 939-768 493). There were also substantial differences in the predicted magnitude of the FMD outbreaks (range: 1 563-8 896) and land area affected (range: 56-447 km(2)) for the different estimated animal distributions. UPGMA clustering indicated there were two main groups of distributions, and one outlier. We recommend that one distribution from each of these two groups be used to model the range of possible outbreaks. Methods included in cluster 1 (such as county-level disaggregation) could be used in conjunction with any of the methods in cluster 2, which included kriging, NDVI split by ecoregion, or disaggregation at the regional level, to represent the variability in the model predicted outbreak distributions. How animal populations are represented needs to be considered in all spatial disease spread models.

Moment Equations and Dynamics of a Household SIS Epidemiological Model

An SIS epidemiological model of individuals partitioned into households is studied, where infections take place either within or between households, the latter generally happening much less frequently. The model is explored using stochastic spatial simulations, as well as mathematical models which consist of an infinite system of ordinary differential equations for the moments of the distribution describing the proportions of individuals who are infectious among households. Various moment-closure approximations are used to truncate the system of ODEs to finite systems of equations. These approximations can sometimes lead to a system of ill-behaved ODEs which predict moments which become negative or unbounded. A reparametrization of the ODEs is then developed, which forces all moments to satisfy necessary constraints. Changing the proportion of contacts within and between households does not change the endemic equilibrium, but does affect the amount of time it takes to approach the fixed point; increasing the proportion of contacts within households slows the spread of the infection toward endemic equilibrium. The system of moment equations does describe this phenomenon, although less accurately in the limit as the proportion of between-household contacts approaches zero. The results indicate that although controlling the movement of individuals does not affect the long-term frequency of an infection with SIS dynamics, it can have a large effect on the time-scale of the dynamics, which may provide an opportunity for other controls such as immunizations to be applied.

Simulating the spatial dynamics of foot and mouth disease outbreaks in feral pigs and livestock in Queensland, Australia, using a susceptible-infected-recovered cellular automata model

We describe an approach to modelling the spatio-temporal spread of foot and mouth disease through feral animal and unfenced livestock populations. We used a susceptible-infected-recovered model, implemented in a cellular automata framework, to assess the spread of FMD across two regions of Queensland, Australia. Following a sensitivity analysis on the infectious states, scenario analyses were conducted using feral pigs only as the susceptible population, and then with the addition of livestock, and initiated in the wet season and in the dry season. The results indicate that, depending on the season the outbreak is initiated, and without the implementation of control measures, an outbreak of Foot and Mouth Disease around Winton could continue unchecked, while an outbreak around Cape York may die out naturally. The approach explicitly incorporates the spatial relationships between the populations through which the disease spreads and provides a framework by which the spread of disease outbreaks can be explored through varying the model parameters. It highlights the emergence and importance of spatio-temporal patterns, something that previous modelling of FMD in feral animal and unfenced livestock populations has lacked.

Understanding foraging behaviour in spatially heterogeneous environments

The role of stochasticity and spatial heterogeneity in foraging systems is investigated. We formulate a spatially explicit model which describes the behaviour of grazing animals in response to local information using simple stochastic rules. In particular the model reflects the biology in that decisions to move to a new location are based on visual assessment of the sward height in a surrounding neighbourhood, whilst the decision to graze the current location is based on the residual sward height and olfactory assessment of local faecal contamination. It is assumed that animals do not interact directly, but do so through modification of, and response to a common environment. Spatial heterogeneity is shown to have significant effects including reducing the equilibrium intake rate and increasing the optimal stocking density, and must therefore be taken into account by resource managers. We demonstrate the relationship between the stochastic spatial model and its non-spatial deterministic counterpart, and in the process derive a moment-closure approximation to the full process, which can be regarded as an intermediate, or pseudo-spatial model. The role of spatial heterogeneity is emphasized, and better understood by comparing the results obtained from each approach. The relative efficiency of random and directed searching behaviour in spatially heterogeneous environments is explored for both clean and contaminated pastures, and the impact of faecal avoidance behaviour assessed.

Analysis of disease progress as a basis for evaluating disease management practices

The relationship between epidemiology and disease management is long-standing but sometimes tenuous. It may seem self-evident that improved understanding of epidemic processes will lead to more effective control practices but this remains a testable proposition rather than demonstrated reality. A wide range of models differing in mathematical sophistication and computational complexity has been proposed as a means of achieving a greater understanding of epidemiology and carrying this through to improved management. The potential exists to align these modeling approaches to evaluation of control practices and prediction of the consequent epidemic outcomes, but these have yet to make a major impact on practical disease management. For the immediate future simpler pragmatic approaches for analysis of disease progress, using nonlinear growth functions and/or integrated measures such as area under disease progress curves, will play a key role in informing tactical and strategic decisions on control treatments. These approaches have proved useful in describing control effectiveness and, in some cases, optimizing or changing control practices.

Effects of dispersal mechanisms on spatio-temporal development of epidemics

The nature of pathogen transport mechanisms strongly determines the spatial pattern of disease and, through this, the dynamics and persistence of epidemics in plant populations. Up to recently, the range of possible mechanisms or interactions assumed by epidemic models has been limited: either independent of the location of individuals (mean-field models) or restricted to local contacts (between nearest neighbours or decaying exponentially with distance). Real dispersal processes are likely to lie between these two extremes, and many are well described by long-tailed contact kernels such as power laws. We investigate the effect of different spatial dispersal mechanisms on the spatio-temporal spread of disease epidemics by simulating a stochastic Susceptible-infective model motivated by previous data analyses. Both long-term stationary behaviour (in the presence of a control or recovery process) and transient behaviour (which varies widely within and between epidemics) are examined. We demonstrate the relationship between epidemic size and disease pattern (characterized by spatial autocorrelation), and its dependence on dispersal and infectivity parameters. Special attention is given to boundary effects, which can decrease disease levels significantly relative to standard, periodic geometries in cases of long-distance dispersal. We propose and test a definition of transient duration which captures the dependence of transients on dispersal mechanisms. We outline an analytical approach that represents the behaviour of the spatially-explicit model, and use it to prove that the epidemic size is predicted exactly by the mean-field model (in the limit of an infinite system) when dispersal is sufficiently long ranged (i.e. when the power-law exponent a</=2).

Combining endogenous and exogenous spatial variability in analytical population models

Analytically tractable models of dynamics in continuous space rarely incorporate both endogenous and exogenous spatial heterogeneity. We use spatial moment equations in combination with simulation models to analyze the combined effects of endogenous and exogenous variability on population viability in a simple single-population model where landscape heterogeneity and local population density both affect mortality rate. The equations partition the effects of heterogeneity into an effect of local crowding and an effect of habitat association caused by differential mortality. Exogenous heterogeneity in mortality rate increases population viability through habitat association and decreases it through increased crowding; the net effect of exogenous heterogeneity is generally to improve population viability. This result is contrary to some (but not all) conclusions in the literature, which usually focus on the effects of fragmentation rather than the benefits of refuges to short-dispersing individuals.

Pair-edge approximation for heterogeneous lattice population models

To increase the analytical tractability of lattice stochastic spatial population models, several approximations have been developed. The pair-edge approximation is a moment-closure method that is effective in predicting persistence criteria and invasion speeds on a homogeneous lattice. Here we evaluate the effectiveness of the pair-edge approximation on a spatially heterogeneous lattice in which some sites are unoccupiable, or "dead". This model has several possible interpretations, including a spatial SIS epidemic model, in which some sites are occupied by immobile host-species individuals while others are empty. We find that, as in the homogeneous model, the pair-edge approximation is significantly more accurate than the ordinary pair approximation in determining conditions for persistence. However, habitat heterogeneity decreases invasion speed more than is predicted by the pair-edge approximation, and the discrepancy increases with greater clustering of "dead" sites. The accuracy of the approximation validates the underlying heuristic picture of population spread and therefore provides qualitative insight into the dynamics of lattice models. Conversely, the situations where the approximation is less accurate reveals limitations of pair approximation in the presence of spatial heterogeneity.

Analytical methods for predicting the behaviour of population models with general spatial interactions

Many biologists use population models that are spatial, stochastic and individual based. Analytical methods that describe the behaviour of these models approximately are attracting increasing interest as an alternative to expensive computer simulation. The methods can be employed for both prediction and fitting models to data. Recent work has extended existing (mean field) methods with the aim of accounting for the development of spatial correlations. A common feature is the use of closure approximations for truncating the set of evolution equations for summary statistics. We investigate an analytical approach for spatial and stochastic models where individuals interact according to a generic function of their distance; this extends previous methods for lattice models with interactions between close neighbours, such as the pair approximation. Our study also complements work by Bolker and Pacala (BP) [Theor. Pop. Biol. 52 (1997) 179; Am. Naturalist 153 (1999) 575]: it treats individuals as being spatially discrete (defined on a lattice) rather than as a continuous mass distribution; it tests the accuracy of different closure approximations over parameter space, including the additive moment closure (MC) used by BP and the Kirkwood approximation. The study is done in the context of an susceptible-infected-susceptible epidemic model with primary infection and with secondary infection represented by power-law interactions. MC is numerically unstable or inaccurate in parameter regions with low primary infection (or density-independent birth rates). A modified Kirkwood approximation gives stable and generally accurate transient and long-term solutions; we argue it can be applied to lattice and to continuous-space models as a substitute for MC. We derive a generalisation of the basic reproduction ratio, R(0), for spatial models.

Spatial heterogeneity and the stability of reaction states in autocatalysis

The impact of stochasticity and spatial heterogeneity on the quadratic autocatalytic system is studied. In a nonspatial setting the reactive state of the system is found to be unstable in small volumes where internal fluctuations drive the system to the unreactive state. This phenomena is of potential importance to the stability of reactions in biological cells. A simple spatial model is constructed by linking N nonspatial models via migration of reactants controlled by a mixing rate lambda. Simulation of this stochastic process demonstrates the importance of such mixing in controlling the impact of internal fluctuations on the stability of the autocatalytic reaction. For high mixing rate the mean reactant levels in equilibrium correspond to the well-mixed deterministic system, although a significant degree of spatial heterogeneity remains. For intermediate mixing rates, mean reactant levels vary continuously with lambda, where the interaction of internal fluctuations with limited spatial mixing modifies the reactive states of the deterministic system. However, there is a threshold below which mixing is unable to control internal fluctuations which drive the system into the unreactive state. Thus a critical minimum level of communication between the cells is required to stabilize the reaction across the entire system. Approximate analytic results, obtained using moment-closure techniques, support these findings and demonstrate the relationship between the spatial stochastic and nonspatial deterministic models.

Models of Mycobacterium bovis in wildlife and cattle

Following the first model of wildlife tuberculosis (in European badgers) there has been a spate of papers modelling wildlife TB. These have looked at population parameters and disease dynamics in the badger and possum. Recent papers in particular have looked at various methods of controlling the wildlife vector to reduce the incidence of TB in cattle. The author examines the role of modelling to show what insights it has given us, which issues have not been addressed, and where the shortfalls lie. Particular attention will be paid to a comparison between models of badgers and possums, and between simple and more complex models, and possible areas of future research will be revealed.

Pair approximation for lattice models with multiple interaction scales

Pair approximation has frequently proved effective for deriving qualitative information about lattice-based stochastic spatial models for population, epidemic and evolutionary dynamics. Pair approximation is a moment closure method in which the mean-field description is supplemented by approximate equations for the frequencies of neighbor-site pairs of each possible type. A limitation of pair approximation relative to moment closure for continuous space models is that all modes of interaction between individuals (e.g., dispersal of offspring, competition, or disease transmission) are assumed to operate over a single spatial scale determined by the size of the interaction neighborhood. In this paper I present a multiscale pair approximation which allows different sized neighborhoods for each type of interaction. To illustrate and test the approximation I consider a spatial single-species logistic model in which offspring are dispersed across a birth neighborhood and established individuals have a death rate depending on the population density in a competition neighborhood, with one of these neighborhoods nested inside the other. Analysis of the steady-state equations yields several qualitative predictions that are confirmed by simulations of the model, and numerical solutions of the dynamic equations provide a close approximation to the transient behavior of the stochastic model on a large lattice. The multiscale pair approximation thus provides a useful intermediate between the standard pair approximation for a single interaction neighborhood, and a complete set of moment equations for more spatially detailed models.

Coupling Disease-Progress-Curve and Time-of-Infection Functions for Predicting Yield Loss of Crops

ABSTRACT A general approach was developed to predict the yield loss of crops in relation to infection by systemic diseases. The approach was based on two premises: (i) disease incidence in a population of plants over time can be described by a nonlinear disease progress model, such as the logistic or monomolecular; and (ii) yield of a plant is a function of time of infection (t) that can be represented by the (negative) exponential or similar model (zeta(t)). Yield loss of a population of plants on a proportional scale (L) can be written as the product of the proportion of the plant population newly infected during a very short time interval (X'(t)dt) and zeta(t), integrated over the time duration of the epidemic. L in the model can be expressed in relation to directly interpretable parameters: maximum per-plant yield loss (alpha, typically occurring at t = 0); the decline in per-plant loss as time of infection is delayed (gamma; units of time(-1)); and the parameters that characterize disease progress over time, namely, initial disease incidence (X(0)), rate of disease increase (r; units of time(-1)), and maximum (or asymptotic) value of disease incidence (K). Based on the model formulation, L ranges from alphaX(0) to alphaK and increases with increasing X(0), r, K, alpha, and gamma(-1). The exact effects of these parameters on L were determined with numerical solutions of the model. The model was expanded to predict L when there was spatial heterogeneity in disease incidence among sites within a field and when maximum per-plant yield loss occurred at a time other than the beginning of the epidemic (t > 0). However, the latter two situations had a major impact on L only at high values of r. The modeling approach was demonstrated by analyzing data on soybean yield loss in relation to infection by Soybean mosaic virus, a member of the genus Potyvirus. Based on model solutions, strategies to reduce or minimize yield losses from a given disease can be evaluated.

An effective sample size for predicting plant disease incidence in a spatial hierarchy

ABSTRACT For aggregated or heterogeneous disease incidence, one can predict the proportion of sampling units diseased at a higher scale (e.g., plants) based on the proportion of diseased individuals and heterogeneity of diseased individuals at a lower scale (e.g., leaves) using a function derived from the beta-binomial distribution. Here, a simple approximation for the beta-binomial-based function is derived. This approximation has a functional form based on the binomial distribution, but with the number of individuals per sampling unit (n) replaced by a parameter (v) that has similar interpretation as, but is not the same as, the effective sample size (n(deff) ) often used in survey sampling. The value of v is inversely related to the degree of heterogeneity of disease and generally is intermediate between n(deff) and n in magnitude. The choice of v was determined iteratively by finding a parameter value that allowed the zero term (probability that a sampling unit is disease free) of the binomial distribution to equal the zero term of the beta-binomial. The approximation function was successfully tested on observations of Eutypa dieback of grapes collected over several years and with simulated data. Unlike the beta-binomial-based function, the approximation can be rearranged to predict incidence at the lower scale from observed incidence data at the higher scale, making group sampling for heterogeneous data a more practical proposition.