Education Center | Plant Disease Management Simulations
Management of Potato Late Blight: Simulation with Lateblight




Lateblight

Model Description

Simulation by
W. E. Fry, M. G. Milgroom, M. A. Doster,
J. A. Bruhn, and R. I. Bruck

Java adaptation by
Joshua M. Goldfarb and Phil A. Arneson

Lateblight source code: http://dspace.library.cornell.edu/handle/1813/2492/

© 2002 Cornell University
All rights reserved




This version of Lateblight was adapted from the simulation originally written by J. A. Bruhn et al. (1980) and modified in various ways over the years by different people working in the research program of W. E. Fry, Department of Plant Pathology, Cornell University. We have attempted to retain as much as possible the original structure of the model, but some small modifications have been necessary to enhance its pedagogical value. While the simulation is sufficiently realistic for teaching purposes, this version should not be trusted as a research tool or a management decision making aid.

In order to make "Lateblight" respond as realistically as possible to environmental and management variables, the development of the pathogen is simulated mechanistically, closely following the Phytophthora infestans life cycle. The model's state variables represent the major morphological stages in the pathogen's development:

Sporangia (landed on susceptible sites)
Zoospores (following indirect germination)
Germinated sporangia
Germinated zoospores
Latent lesions
Sporulating, expanding lesions
Sporulating, inactive lesions (not expanding)
Inactive lesions (nonsporulating, necrotic tissue)
Sporangia available for dispersal
Dispersed sporangia

The model's transfer functions describe different developmental processes:

Indirect germination
Direct germination
Infection
Lesion expansion
Sporulation
Lesion inactivation (cessation of expansion)
Necrosis (cessation of sporulation)
Dispersal of sporangia
Catch of sporangia

The model operates on a daily time step.

Temperature Effects

The effect of temperature on disease development is simulated independently for each of six different developmental processes using polynomial equations (Table 1).

Table 1. Rates of Pathogen Development as a Function of Temperature (C)
Temperature
Effect
Equation Min Max
Direct Germination TempEffDirect = 2.1 - 0.48 TempLeafWet + 0.031 TempLeafWet2 - 0.00045 TempLeafWet3 - 4.5 10-6 TempLeafWet4 14 29
Indirect Germination TempEffIndir = 0.0134 + 0.0698 TempLeafWet + 0.0185 TempLeafWet2 - 0.00223 TempLeafWet3 + 7.48 10-5 TempLeafWet4 - 7.68 10-7 TempLeafWet5 0 24
Zoospore Germination TempEffZoosp = 0.00494 + 0.364 TempLeafWet - 0.05813 TempLeafWet2 + 0.00452 TempLeafWet3 - 1.617 10-4 TempLeafWet4 + 2.068 10-6 TempLeafWet5 0 28
Infection TempEffInfec = 0.4846 + 0.3735 TempHumid90 - 0.03556 TempHumid902 + 0.001509 TempHumid903 - 2.405 10-5 TempHumid904 1.5 30
Lesion Expansion TempEffLesEx = sqrt (8.4 10-4 + 0.1261 DailyMeanTemp - 0.03079 DailyMeanTemp2 + 0.00308 DailyMeanTemp3 - 1.152 10-4 DailyMeanTemp4 + 1.414 10-6 DailyMeanTemp5) 0 33
Sporulation TempEffSporu = - 0.3882 + 0.2482 TempHumid90 - 0.03817 TempHumid902 + 0.00252 TempHumid903 - 5.32 10-5 TempHumid904 0 26
Parameters of the polynomial equations were estimated by multiple regression using the data of Crosier (1934).

Sporangial Germination

Sporangia germinate by two pathways, direct and indirect, depending on the ambient temperature during the period of leaf wetness: (See Table 1.)

GermSporangia = LandedSporangia * (DIRMAX * TempEffDirect)

Zoospores = LandedSporangia * (INDMAX * ZSPR * TempEffIndir)

where DIRMAX = 0.29 sporangia/day, INDMAX = 0.81 sporangia/day, and ZSPR = 6.567 zoopores/sporangium (Marks, 1965; Ullrich and Schober, 1972; Warren and Calhoun, 1975). The temperature effect is determined by the appropriate polynomial function in Table 1.

The zoospores germinate only while the leaves are wet and in response to the mean temperature during the period of leaf wetness:

GermZoospores = Zoospores * ZOOMAX * TempEffZoosp
where ZOOMAX = 0.96 zoospores/day.

Infection

Infection efficiency is a function of leaf wetness, mean temperature during the period of leaf wetness, the susceptibility of the potato cultivar, the type of fungicide and its residue level, and the amount of uninfected tissue remaining.

LatentInfections = (INFMAX_SP * GermSporangia) + (INFMAX_ZSP * GermZoospores) + (TempEffInfec * LeafWetEff * Susceptibility * FungicideEffect * (1.0 - BlightProportion))
where INFMAX_SP = 0.1 infections/germinated sporangium/day and INFMAX_ZSP = 0.01 infections/germinated zoospore/day (Lapwood and McKee, 1966). The temperature effect is determined by the appropriate polynomial function in Table 1. Maximum infection efficiency occurs at 21 C with 24 hours of leaf wetness, but the effects of temperature on infection efficiency are small between 10 and 22 C. No infections occur below 1.5 C or above 30 C (Crosier, 1934). Infection occurs only if the leaves are wet. The susceptibility factors for the susceptible, moderately resistant, and resistant cultivars are 1.0, 0.9, and 0.8, respectively. The protectant and systemic fungicides have different dose responses, the systemic requiring a lower residue to achieve the same effect. (See "Fungicide effects" below.)

Lesion development

Lesions are divided into four types:

Latent lesions
Sporulating, expanding lesions
Sporulating, inactive lesions (not expanding)
Inactive lesions (nonsporulating, necrotic tissue)
All successful infections produce latent lesions. The latent periods for the susceptible, moderately resistant, and resistant cultivars are 6, 6, and 7 days, respectively. At the end of the latent period, the lesion diameter is set at 2 mm, following which the sporulating lesions expand for four days:
LesionDiam[age] = LesionDiam[age-1] + (MAXLES * TempEffLesEx * LesionExp * IncrFact * FungicideFact * (1.0 - BlightProportion))
where MAXLES = 4.35 mm/day (Guzman-N, 1964; Lapwood, 1961). the temperature effect is determined by the appropriate polynomial function in Table 1, and LesionExp is 1.5, 1.4, and 1.4 for the susceptible, moderately resistant, and resistant cultivars, respectively. IncrFact is a factor that decreases linearly from 1.0 to 0.0 over the five days following the latent period. FungicideFact is a factor between zero and one calculated from the dose response of the fungicide. (See "Fungicide effects" below.)

The total lesion area assumes circular lesions:

LesionArea[age] = LesionNumber[age] * LesionDiam[age] * LesionDiam[age] * PIBY4
where PIBY4 = 0.7854.

In the original model, the numbers of lesions of each age were passed daily into the next age class:

LesionNumber[age] = LesionNumber[age-1]
In the current version of the model, this has been modified to pass a proportion of the lesions in each class on to the next age class each day:
LesionNumber[age] = LesionNumber[age] + LesionNumber[age-1] * TRANSFER
LesionNumber[age-1] = LesionNumber[age-1] - LesionNumber[age-1] * TRANSFER
where TRANSFER = 0.38 (essentially a "tuning parameter", adjusted so that the model's output gave a good fit to disease progress observed in the field trials in 1977 and 1978 near Freeville, NY.)

Sporulation

Sporulation occurs on lesions from the end of the latent period until they are 14 days old:

AttachedSporangia = SPRMAX * SporulationFct * TempEffSporu * HumidityFct * SysFungicideFct * TotalLesionArea * (1.0 - BlightProportion)
where SPRMAX = 763.2 (Lapwood, 1961) and the SporulationFct is 0.55, 0.45, and 0.45 for the susceptible, moderately resistant, and resistant cultivars, respectively. The temperature effect is determined by the appropriate polynomial function in Table 1. Sporulation occurs only if the relative humidity is above 90% (DeWeille, 1963). The effect of the systemic fungicide on sporulation is determined by the dose response discussed below.

Lesion inactivation

On the fifteenth day following its initiation, a lesion ceases sporulation and its area is added to the inactive lesion area, which is added to the total blighted area when calculating the blight proportion. All lesions become inactive if the daily mean temperature exceeds 41 C for three consecutive days (Wallin and Hoyman, 1958).

Dispersal of sporangia

No attempt has been made in "Lateblight" to simulate the biophysical mechanisms of dispersal of sporangia. All of the spores that are produced are dispersed, but the proportion of sporangia surviving are first adjusted for temperature dependent mortality:

if (DailyMeanTemp > 21.0)
     SporeSurvival = exp (6.0 * (-0.0256 * DailyMeanTemp + 0.537)) ;
else
     SporeSurvival = 1.0
The numbers of sporangia landing on susceptible sites (and thus available for infection) depend on the hours of leaf wetness (Crosier, 1934; Glendinning, et al., 1963; Melhus, 1915), the amounts washed off by heavy rains (Waggoner and Horsfall, 1969), and by the amount of foliage present:
LandedSporangia = AttachedSporangia * SporeSurvival * LeafWetCatch * RainWashFact * (LeafAreaIndex / MaxLeafAreaIndex)
where LeafWetCatch = 0.28.0/LeafWetHrs and RainWashFact = 0.3Rainfall*10.0.

Host growth

The increase in leaf area index is simply linear up to the day determined by the characteristics of the cultivar (early-, mid-, or late-season) and reaches a maximum, also characteristic of the cultivar:

Early Mid Late
Day foliage growth stops 44 52 59
Maximum leaf area index 2.5 3.5 4.5

Tuber production is a function of temperature, cultivar, time. Tuber initiation can occur above 0 and below 30 C. A temperature factor for tuber initiation is generated by the quadratic:

TuberInitFact = - 0.00911 + 0.1038 Temp - 0.00296 Temp2
Physiological time is then accumulated at a rate that depends on the cultivar:
TuberInitTime = TuberInitTime + CvTuberFact * TuberInitFact
where CvTuberFact is 1.600, 1.404, and 1.333 for the early-, mid-, and late-season cultivars, respectively. Tuber production is initiated when 45 physiological time units have accumulated. The rate of tuber production is also a quadratic function of temperature:
TuberProdFact = - 0.3837 + 0.1261 Temp - 0.00298 Temp2
Again physiological time is accumulated for the production of tubers:
TuberProdTime = TuberProdTime + CvTuberFact * TuberProdFact
The production of tubers in tons/hectare is then simply:
Tubers = 0.56 * TuberProdTime

Fungicide effects

The details of the simulation of the attenuation of the protectant fungicide residues are described by Bruhn and Fry (1982). Basically the canopy is divided into four layers, and the fungicide residues are redistributed among the layers in response to rain. The residues of the systemic fungicide are a linear function of time after spraying only (Milgroom et al., 1988).

The only effect of the protectant fungicide is the inhibition of infection (Bruhn and Fry, 1981), whereas the systemic fungicide not only inhibits infection but also reduces the rate of lesion expansion and inhibits sporulation (Milgroom and Fry, 1988).

Weather

The following environmental variables are read from an external weather data file:

Rainfall: Daily rainfall in inches
DailyMeanTemp: Daily mean temperature (degrees Celsius)
Humid90Hours: Hours per day of relative humidity above 90%
TempHumid90: Mean temperature (C) during the above humid period

In addition, two environmental variables are approximated with the following calculations:

LeafWetHours = Humid90Hours + 3.0
TempLeafWet = DailyMeanTemp + TempHumid90) / 2.0



REFERENCES

Bruhn, J. A., R. I. Bruck, W. E. Fry, P. A. Arneson, and E. V. Keokosky. 1980. Lateblight: A Disease Management Game. Computer program and manual, Department of Plant Pathology. Cornell University, Ithaca, NY.

Bruhn, J. A. and W. E. Fry. 1981. Analysis of potato late blight epidemiology by simulation modeling. Phytopathology 71:612-616.

Bruhn, J. A. and W. E. Fry. 1982. A mathematical model of the spatial and temporal dynamics of chlorothalonil residues on potato foliage. Phytopathology 72:1306-1312.

Crosier, W. 1934. Studies in the biology of Phytophthora infestans (Mont.) DeBary. Memoir 155. Cornell University Agricultural Experiment Station. March 1934. 37 pp.

DeWeille, G. A. 1963. Laboratory results regarding potato late blight and their significance in the epidemiology of blight. Eur. Potato Jour. 6:121-130.

Glendinning, D. R., J. A. MacDonald and J. Grainger. 1963. Factors affecting the germination of sporangia of Phytophthora infestans. Trans. Brit. Mycol. Soc. 45:595-603.

Guzman-N., J. 1964. Nature of partial resistance of certain clones of three Solanum species to Phytophthora infestans. Phytopathology 54:1398-1404.

Lapwood, D. H. 1961. Potato haulm resistance to Phytophthora infestans. II. Lesion production and sporulation. Ann. Appl. Biol. 49:316-330.

Lapwood, D. H. and R. K. McKee. 1966. Dose-response relationships for infection of potato leaves by zoospores of Phytophthora infestans. Trans. Brit. Mycol. Soc. 49:679-686.

Marks, G. E. 1965. The cytology of Phytophthora infestans. Chromosoma 16:681-692.

Melhus, I. E. 1915. Germination and infection with the fungus of the late blight of potato (Phytophthora infestans). Wisc. Agr. Exp. Sta. Res. Bull. 37:1-64.

Milgroom, M. G., C. E. McCulloch, and W. E. Fry. 1988 Distribution and temporal dynamics of metalaxyl in potato foliage. Phytopathology 78:555-559.

Milgroom, M. G. and W. E. Fry. 1988. A model for the effects of metalaxyl on potato late blight epidemics. Phytopathology 78:559-565.

Ullrich, J. and B. Schoger. 1972. Zoosporenzahl und Sporangiengrube bei Phytophthora infestans (Mont.) DeBary. Phytopath. Z. 74:268-271.

Waggoner, P. E. and J. G. Horsfall. 1969. EPIDEM: a simulator of plant disease written for a computer. Bull. Conn. Agr. Exp. Sta., New Haven, No. 698. 80 pp.

Wallin, J. R. and W. G. Hoyman. 1958. Influence of post-inoculum air temperaturemaxima on survival of Phytophthora infestans on potato leaves. Amer. Potato Jour. 35:769-773.

Warren, R. C. and J. Calhoun. 1975. Viability of sporangia of Phytophthora infestans in relation to drying. Trans. Brit. Mycol. Soc. 64:73-78.


Contact: Phil A. Arneson
Last updated: April 8, 2004
Copyright 2002 Cornell University