In Region 1, the early warning model had a sensitivity of 62% (16 of 26 epidemics). Prediction improved dramatically for the late warning model, which had a sensitivity of 96% (25 of 26 epidemics). In Region 2, the early warning model had a sensitivity of 73% (32 of 44 epidemics). The late warning model, however, had a lower sensitivity of 66%.
These findings indicate that weather conditions at relatively coarse temporal and spatial resolutions can be used to predict RRV disease epidemics with sufficient accuracy and lag time for public health planning. Figure 7 summarizes the main results. Both regions had models with more than 70% sensitivity, and in Region 2 epidemics were predicted with 73% sensitivity as early as November (when cases typically commence).
The Influence of Rainfall Pattern
The timing, duration, and extent of rainfall events in epizootic regions are critical indicators of outbreaks. 29 Two different rainfall patterns were evident in this study: excess “winter” rainfall and excess “summer” rainfall. Breeding of the main transmission vector, C annulirostris, occurs from mid-spring to late autumn, 10 and populations typically peak in January or February. 30 Reports from previous epidemics in these regions have linked heavy rainfall or extensive flooding in the December to February period with intense breeding of Culex mosquitoes, and hence to a spillover of infection to human populations. Large epidemics in the years 1928, 1956, 1984, and 1971 demonstrate this pattern. 31,27,32,33 In the first epidemic year of the study period, 1992–1993, heavy rainfall (nearly twice the long-term mean) commenced in August, and above average rainfall was sustained throughout the spring and up until late January (see Figure 5).
Nicholls 34 has noted that severe outbreaks of RRV disease in temperate Australia appear to follow heavy summer rainfall, suggesting an association with La Niña, the cool phase of the ENSO cycle. Harley and Weinstein 35 found no association between the southern oscillation index (SOI) and RRV disease outbreak years for Australia as a whole, although some RRV disease outbreaks in southeastern Australia have been positively correlated with September values of the SOI. 36 We found August SST to be a better predictor of rainfall excess, and hence epidemics, than the SOI. This was not surprising, as calculation of the SOI incorporates air pressure values from both sides of the Pacific basin (Darwin and Tahiti), whereas SST records the rise and fall of temperatures directly in the El Niño region of the eastern Pacific (Table 2, note ∥).
Not all outbreaks in the Murray area have followed heavy summer rainfall. 4 In 1980–1981, abnormally high spring rainfalls were recorded before outbreaks. In the 1996–1997 epidemic year, total rainfall over summer was below the period average in both regions. Instead, higher than average rainfall commenced as early as June, and continued until late August in Region 1 and September in Region 2 (Figure 6). Although high rainfalls were recorded in some of these months in other years (1991–1992 and 1995–1996), they were interspersed with lower than average values in either July or August. These lower values may have been a factor in aborting a buildup in these years.
We speculate that sustained winter and spring rainfall, even in the absence of excess summer rainfall, could enhance transmission by two mechanisms. First, it would allow a longer time period for amplification of virus levels between the initiating vector (Aedes species) and host populations. Floodwater Aedes maintain the virus by transovarial transmission of infection through embryonated eggs. 37 These mosquitoes overwinter as drought-resistant eggs in mud flats and creek beds. 38 After heavy winter rainfall, the females emerge and infect the vertebrate hosts. 37 In 1996–1997, minimum temperatures during July and August were higher than average in both regions. The combination of these two factors is likely to have provided for early and prolific breeding of Aedes populations and an extended period of virus buildup, thus increasing transmission potential. Second, the prolonged heavy rainfall (which resulted in widespread flooding into October) acted to raise the water table in the regions, reducing absorption and runoff. As a consequence, ground pools were able to remain into summer, despite low summer rainfalls, providing sites for Culex and Aedes breeding. These species extended the infection to humans, and epidemics ensued.
Irrigation practices, with a special release from the Hume Reservoir, appear likely to have contributed to the maintenance of wet conditions throughout the 1996–1997 summer. The Hume modulates water levels for the Murray River, influencing flood patterns downstream. A “special release” of an excessive volume of water from the Hume was authorized in October 1996 (owing to fears of a crack in the reservoir wall), which artificially maintained the Murray River flow at flood levels for more than 30 days. Under normal circumstances, a maximum irrigation release of 25,000 ml per day is permitted, and negligible flooding occurs. Flow levels in October of 1996 were persistently above 25,000 ml per day (for 2 weeks they were more than 90,000 ml per day), and did not return to normal for 6 weeks. Under this volume of water the area inundated increased rapidly, 39 and billabongs and low depressions around the river filled up.
Extrinsic (ie, climatic) factors need to be combined with host-virus population dynamics in epidemic prediction. 40 Below-average rainfall in the spring of the preceding year was a necessary, but not sufficient, explanatory variable in all the models. We attribute this to the influence of vector activity on vertebrate host immunity. Higher than normal rainfalls during an October–December period would initiate and support mosquito virus activity and result in the infection of a proportion of the vertebrate population. As the breeding cycle of the primary vertebrate host (macropods) takes more than 1 year to complete, 41 the pool of susceptible vertebrates in the following year would be reduced, virus amplification would be minimal, and the probability of human cases would also be small. Conversely, in the year after an epidemic year, the high level of vertebrate population immunity would be sufficient to reduce the probability of successive epidemic years to very low levels, even if other climatic conditions were suitable (as occurred in 1993–1994 in the study regions). This finding has also been demonstrated by studies on RRV disease epidemic activity in arid northwestern Australia. 42
Limitations of the Study
The principal limitation of this study lies with the notification data, which were reported by place of residence rather than place of suspected infection. We made the assumption that in rural areas with relatively large SLAs, work and recreation (and hence transmission) generally occur within the same SLA. However, it is likely that a small portion of cases were misclassified, which would weaken the association. Furthermore, because of a lack of data, we were not able to control for the possible confounding effect of rice irrigation practices, which may vary by SLA. Although we understood mosquito control interventions to have been minimal throughout the period, estimates of the timing and effectiveness of sprays would have strengthened the findings.
The short span of data was a problem for some climate variables. Given the unusual number of El Niño events that occurred during the study period (1990–1994 and 1997–1998), we cannot be certain about the strength of the predictive power of sea surface temperature in the Region 2 early warning model. One La Niña (cool phase) occurred in 1998–1999, but was only of weak intensity during the critical late winter to spring months.
The cross-validation method assessed how well the variables in each of the models were able to predict epidemics, rather than the predictive performance of the estimated coefficients for the 8-year model. Epidemics occurred infrequently in this dataset (in general, in only 2 of the 8 years). In the rotating validation process, when 1 year of epidemic information was removed, only half of the data on epidemics remained to derive the parameter estimates of the predictive model for that year. Given that there are several variables in each model, it may be that some of these predicted the outcome by chance.
Climate, Climate Change, and Health
Beyond the seasonal effects studied here, our changing climate—including interannual cycles and longer-term natural and human-induced cycles—will continue to influence the extrinsic and intrinsic factors that drive vector-borne diseases. Climate projections for 2030 indicate that temperatures may rise by 0.4–2°C over most of Australia. By 2070 they may rise by 1–6°C. 43 Winter rainfall in the southeast of the country may decrease by up to 10% by 2030 (35% by 2070), and summer rainfall in these regions may increase by 10% to 20% by 2030 (35% to 60% by 2070) depending on whether the results of slab or coupled models are consulted. 43
The effect of these changes on the breeding and survival of arthropods and vertebrate host populations, even in a region as “localized” as southeastern Australia, is difficult to anticipate. Russell speculates that Aedes populations in dry areas (such as Region 1) may be adversely affected by decreased winter rainfall, possibly resulting in delayed or precluded virus activity. 6 We would expect this to result in an overall reduction in epidemics in years when the “winter pattern” might have occurred. Conversely, the predicted increase in summer rainfall may increase the availability of mosquito habitat that, combined with higher average temperatures, may lead to higher humidity, a lengthened season of abundance, and greater transmission levels. 6 Understanding how ENSO may change with global climate change is also essential in anticipating the impact of future climate on vector-borne diseases. Currently, climate models have mixed success in estimating interannual influence. 44
Weather forecasts can be used in conjunction with other surveillance techniques to identify conditions suitable for an epidemic of RRV disease. Seasonal arbovirus activity is already monitored in parts of southeastern Australia, with weekly trapping of mosquitoes throughout the Ross River season used to record population profiles and virus isolates. The multistaged approach (early and late warning models) developed in this study enables response plans to be adjusted as forecast certainty increases, allowing health authorities to make the most of limited resources. Although both weather and climate are pivotal in generating the conditions needed for epidemics of arboviral disease, human influence (through environmental modifications such as irrigation) can also contribute to events.
We acknowledge Keith Moodie and others at SILO, Queensland Department of Natural Resources and Mines, for their assistance with compiling the climate data. We acknowledge the Communicable Diseases Network of Australia and New Zealand for the disease data, and the Murray-Darling Basin Commission for providing data on releases from the Hume Reservoir. We thank staff at the Animal Health Science Unit, Department of Agriculture Fisheries and Forestry Australia, for the use of computing resources.
1. Rosen L, Gubler DG, Bennett PH. Epidemic polyarthritis (Ross River) virus infection in the Cook Islands. Am J Trop Med Hyg 1981; 30: 1294–1302.
2. Aaskov JG, Mataika JU, Lawrence GW, Rabukawaqa W, Tucker MJAR, Dalglish DA. An epidemic of Ross River virus infection in Fiji. Am J Trop Med Hyg 1981; 30: 1053.
3. Scrimgeour EM, Aaskoz JG, Matz LR. Ross River virus arthritis in Papua New Guinea. Trans R Soc Trop Med Hyg 1987; 81: 833–834.
4. Mudge PR, Aaskov JG. Epidemic polyarthritis in Australia
, 1980–1981. Med J Aust 1983: 2: 269–273.
5. Westley-Wise VJ, Beard JR, Sladden TJ, Dunn TM, Simpson J. Ross River virus infection on the north coast of New South Wales. Aust N Z J Public Health 1996; 20: 87–92.
6. Russell RC. Mosquito-borne arboviruses in Australia
: the current scene and implications of climate change
for human health. Int J Parasitol 1998; 28: 955–969.
7. Kay BH, Aaskov JG. Ross River virus (epidemic polyarthritis). In: Monath TP, ed. The Arboviruses: Epidemiology and Ecology.
Volume IV. Boca Raton, Florida: CRC Press, 1989:93–112.
8. Woodruff RE, Garner MG, Ebi K. Ross River Virus disease
: the influence of climate
on distribution and trends. Submitted for Publication.
9. Dhileepan K. Mosquito seasonality and arboviral incidence in Murray Valley, southeast Australia
. Med Vet Entomol 1996; 10: 375–384.
10. Dobrotworsky NV. The Mosquitoes of Victoria: (Diptera, Culicidae). London: Melbourne University Press, 1965.
11. Russell RC. Ross River virus: disease trends and vector ecology in Australia
. Bull Soc Vect Ecol 1994; 19: 73–81.
12. Nicholls N, Kariko A. East Australian rainfall events: interannual variation, trends, and relationships with the Southern Oscillation. J Clim 1993; 6: 1141–1152.
13. Mackenzie JS, Broom AK, Calisher CH, et al
. Diagnosis and reporting of arbovirus
infections in Australia
. Commun Dis Intell 1993; 17: 202–207.
14. Tai K, Whelan PI, Patel MS, Currie B. An outbreak of epidemic polyarthritis (Ross River virus disease
) in the Northern Territory during the 1990–1991 wet season. Med J Aust 1993: 158: 522–525.
15. MapInfo Professional
[computer program]. Version 5.0. New York: MapInfo Corporation, 1998.
16. Australian Bureau of Statistics. Integrated Regional Database (IRDB)
. Version 1998. Canberra: Commonwealth of Australia
17. Woodruff RE, Garner MG, Guest CS. Using GIS to study the epidemiology of Ross River virus, an important zoonosis in Australia
. Paper presented at GIS Veterinary Conference, Lancaster, UK, September, 2001.
18. Jeffrey SJ, Carter JO, Moodie KB, Beswick AR. Using spatial interpolation to construct a comprehensive archive of Australian climate
data. Environ Model Softw 2001; 16: 309–330.
19. Hutchinson MF. Interpolating mean rainfall using thin plate smoothing splines. Int J Geog Info Sys 1995; 9: 385–403.
20. Stata Statistical Software
. Edition 6.0. College Station: Stata Press, 1999.
21. Hosmer DW, Lemeshow S. Goodness-of-fit tests for the multiple logistic regression model. Commun Stat 1980; A9: 1043–1069.
22. Fraser JRE, Cunningham AL. Incubation time of epidemic polyarthritis. Med J Aust 1980; 1: 550–551.
23. Scott TW, McLean RG, Francy DB, Card CS. A simulation model for the vector-host transmission system of a mosquito-borne avian virus, Turlock (bunyaviridae). J Med Entomol 1983; 20: 625–640.
24. de Moor PP, Steffens FE. A computer-simulation model of an arthropod-borne virus transmission cycle, with special reference to chickungunya virus. Trans R Soc Trop Med Hyg 1970; 64: 927–934.
25. Mottram P, Kay BH, Kettle DS. The effect of temperature on eggs and immature stages of Culex annulirostris
Skuse (Diptera: Culicidae). J Aust Entom Soc 1986; 25: 131–135.
26. Weinstein P. An ecological approach to public health intervention: Ross River virus in Australia
. Environ Health Perspect 1997; 105: 364–366.
27. Hawkes RA, Boughton CR, Naim HM, Stallman ND. A major outbreak of epidemic polyarthritis in New South Wales during the summer of 1983/1984. Med J Aust 1985; 143: 330–333.
28. Tong S, Bi P, Parton K, Hobbs J, McMichael AJ. Climate
variability and transmission of epidemic polyarthritis. Lancet 1998; 351: 1100.
29. Davies FG, Linthicum KJ, James AD. Rainfall and epizootic Rift Valley Fever. Bull World Health Organ 1985; 63: 941–943.
30. Russell RC. Seasonal activity and abundance of the arbovirus
vector Culex annulirostris
skuse near Echuca, Victoria, in the Murray Valley of Southeastern Australia
1979–1985. Aust J Explor Biol Med Sci 1986; 64: 97–103.
31. Nimmo JR. An unusual epidemic. Med J Aust 1928; 1: 549–550.
32. Anderson SG, French EL. An epidemic exanthem associated with polyarthritis in the Murray Valley. Med J Aust
33. Seglenieks V, Moore BW. Epidemic polyarthritis in South Australia
: report of an outbreak in 1971. Med J Aust 1974; 2: 552–556.
34. Nicholls N. El Niño-Southern. Oscillation and vector-borne
disease. Lancet 1993; 342: 1284–1285.
35. Harley DO, Weinstein P. The Southern Oscillation Index and Ross River virus outbreaks. Med J Aust 1996; 165: 531–532.
36. Maelzer D, Hales S, Weinstein P, Zalucki M, Woodward A. El Niño and arboviral disease prediction
. Environ Health Perspect 1999; 107: 817–818.
37. Lindsay MDA, Broom AK, Wright AE, Johansen CA, Mackenzie JS. Ross River virus isolations from mosquitoes in arid regions of Western Australia
: implications of vertical transmission as a means of persistence of the virus. Am J Trop Med Hyg 1993; 49: 686–696.
38. Marshall ID, ed. The Epidemiology of Murray Valley Encephalitis in Eastern Australia
–Patterns of Arbovirus
Activity and Strategies of Arbovirus
Research in Australia
. Vol. 2. Brisbane: CSIRO Division of Animal Health and Queensland Institute of Medical Research, 1979.
39. River Murray Commission. Report on Flood Plain Management of Murray River Between Hume Dam and Yarrawonga
. Albury: Murray-Darling Basin Commission, September 1976.
40. Hay SI, Myers MF, Burke DS, et al
. Etiology of interepidemic periods of mosquito-borne disease. Proc Natl Acad Sci U S A 2000; 97: 9335–9339.
41. Caughley G, Shepherd N, Short J. Kangaroos. Their Ecology and Management in the Sheep Rangelands of Australia
. Cambridge: Cambridge University Press, 1987.
42. Mackenzie JS, Lindsay MD, Broom AK. Effect of climate
and weather on the transmission of Ross River and Murray Valley encephalitis viruses. Microbiol Aust 2000; 21: 20–24.
Impacts Group. Climate Change
Scenarios for the Australian Region. Melbourne: CSIRO Division of Atmospheric Research, 1996.
44. Hennessy KJ, Suppiah R, Page CM. Australian rainfall changes, 1910–1995. Aust Meteorol Mag 1999: 48: 1–13.
Keywords:© 2002 Lippincott Williams & Wilkins, Inc.
arbovirus; vector-borne; climate; climate change; early warning; prediction; Ross River virus disease; Australia