RELATIONSHIPS OF THE SOUTHERN OSCILLATION OF THE PACIFIC OCEAN AND DROUGHTS AND WILDFIRES IN THE OHIO RIVER VALLEY Daniel A. Yaussy and Elaine Kennedy Sutherland Research Forester and Research Ecologist, respectively, USDA Forest Service, Northeastern Forest Experiment Station, 359 Main Road, Delaware, OH 43015 ABSTRACT: Changes in global climate might alter the frequency, duration, and intensity of catastrophic weather events such as droughts, floods, and violent storms. As the frequency and severity of drought increase, so does the probability of wildfires in forest ecosystems and the area burned. This study investigates the effect of the El Nino/Southern Oscillation (ENSO) climatic phenomenon on the occurrence of drought and wildfire in the Ohio River Valley, USA. The strength of the ENSO phenomenon is measured by the Southern Oscillation Index (SOI) and the Palmer Drought Severity Index (PDSI). We used the area burned by wildfires to associate PDSI and SOI with fire because it is less dependent on human population density than the numbers of fires reported. Although droughts could not be related to the Southern Oscillation, 25 percent of the variation in area burned by spring and fall fires could be explained by SOI and variables computed from PDSI measured more than 5 months before the fire season. INTRODUCTION Changes in global climate might alter the frequency, duration, and intensity of catastrophic weather events such as droughts, floods, violent storms, etc. As the frequency and severity of drought increase, so does the probability of wildfires in forest ecosystems and the area burned. The ability to foresee these catastrophic events would better equip forest managers to make contingency plans and develop fire management policies. Teleconnections between large-scale atmospheric circulation patterns and regional weather patterns might provide forest managers with an important planning tool. For example, in a broad regionally-based study, Simard et al. (1985) found a high correlation with fires in the southern states and the El Nino/Southern Oscillation (ENSO), but found little or no correlation in the northeastern or southwestern United States. Simard and Main (1987) used the ENSO to develop a regression equation to predict the severity of the fire season in the South. Although these preliminary studies showed no relation between the ENSO and wildfires in the southwest, Swetnam and Betancourt (1990) found significant teleconnections between the ENSO and the occurrence of wildfires in the National Forests of Arizona and New Mexico. In this study we investigated the interactions of the ENSO climatic phenomenon, the occurrence of drought, and the area burned by wildfire in the Ohio River Valley, USA. We have loosely defined the Ohio River Valley to include the states of Illinois, Indiana, Kentucky, Missouri, Ohio, and West Virginia. METHODS Our analysis was based on graphical interpretation and the use of the Pearson Correlation test. Correlations were calculated among the Palmer Drought Severity Index (PDSI), Southern Oscillation Index (SOI), and the area burned by wildfire. To provide predictive capabilities, we calculated the correlations between the variable of interest and values of the other variables in the preceding 12 months. DATA We chose the PDSI to measure drought intensity (Palmer 1965). Based on climate factors (temperature, season, region), water supply (precipitation, stored soil moisture), and water demand (evapotranspiration, runoff, percolation), the PDSI is slow to respond to sudden changes in the weather, so the value from one month is highly correlated to the value of the previous month. Monthly PDSIs were obtained from data supplied by the National Climatic Data Center for 1895 through 1988 for the Hydrologic regions of the states in the Ohio River Valley ( Figure 1). The SOI is the calibrated difference between atmospheric sea-level pressure at Tahiti, French Polynesia and Darwin, Australia. When the SOI is low, the atmospheric sea-level pressure is low in Tahiti and high in Darwin. This is called El Nino and corresponds to wet weather in the southwestern United States (Swetnam and Betancourt, 1990). When the SOI is high, it is termed a La Nina and is associated with droughts in the Southwest. Monthly SOIs in this study are those developed and standardized by Wright (1989) for 1895 through 1988. We used a 3-month moving average to smooth the SOI time series to eliminate noise without masking trends, such that: SOIAV(t) = {SOI(t) + SOI(t-1) + SOI(t-2)}/3. This corresponds to 3-month seasonal averages used by others, such as Swetnam and Betancourt (1992). In recent history, fires in the eastern United States have been considered a cultural artifact ( Figure 2); that is, the number of fires is heavily influenced by human population dynamics. For example, the higher the population density, the larger the number of fires set and reported. In an area of constant population, the number of fires should be distributed uniformly, and the area burned by these fires should respond to climatic influences and be less uniform in distribution. Therefore, we chose to study the area burned by wildfires, which is more indicative of the severity of a fire season. Daily records of fire events were available for the Daniel Boone National Forest (DBNF) in Kentucky from January 1970 through November 1992. These data were summed to monthly totals to assess the correlations with the SOIAV and PDSI data sets. RESULTS PDSIs for each hydrologic region for each month were compared to the SOIAV time series lagged from 1 to 12 months. The resulting correlations were low (r < 0.1) but highly significant (p < 0.0001) due to the large number of data points (n = 49,631). On inspecting the data during extreme drought years (1931 and 1953-54), we concluded that the hydrologic regions bordering the Ohio River may have been more affected by drought than those farther from the river ( Figure 3). We recalculated the correlations, including only those regions adjacent to the river. This raised the correlations slightly (maximum r = .18), but there was no meaningful relationship between PDSI and combinations of lagged SOIAV, latitude and longitude by region, or month. The highest correlation between PDSI and SOIAV was revealed when considering extreme drought events (PDSI more than two standard deviations from the mean). There was a 0.56 correlation (n = 39) between an extreme PDSI and SOIAV from the previous 2 months. Considering Figure 4 and the small number of observations, we doubt the usefulness of this statistic. Values for acreage burned obtained from the daily record of fire events for the DBNF were summed by month to correspond to the monthly records of SOIAV and PDSI ( Figure 5). We correlated acreage burned (AB) with the SOIAV at lags up to 12 months. The most significant correlation was at 5 months' lag, r = -0.20, p < 0.0018 with n = 228 (only the years 1970 through 1988 were present in all data sets). The 5-month lag corresponds to a similar lag found by Swetnam and Betancourt (1992) for the Southwest. They found a relationship between the June fire season in the Southwest and SOIAV from February. The correlation between lagged PDSIs and AB decreased steadily as lag increased; that is, the relationship between PDSI and AB was most significant in the month preceding the month of the fires (r = -0.20, p < 0.0026, n = 228). This relationship is of little long-range predictive use. Figure 6 shows that the DBNF has two fire season each year -- March-April and October-November. Since the AB for the fire seasons is of greatest interest in fire management planning, these values for March and April and those for October and November were summed and analyzed as representing the spring and fall fire season, respectively. We found that the fall AB was most highly correlated to the SOIAV of the previous June (r = -0.61, p < 0.0075, n = 18), the 5-month lag reported earlier. However, the spring AB was most significantly correlated to the SOIAV from 10 months previous -- again, June (r = -.039, p < 0.0879, n = 20). While re-examining Figure 5, we noticed that the large values of AB preceded the lowest values of PDSI. This indicated that only several months of declining PDSIs were needed to create conditions favorable for fire, which is quite logical. The larger values of AB followed periods of high PDSI, which would stimulate lush periods of foliage growth (fine fuels). The severe fire seasons of the fall of 1981 and the spring of 1988 preceded the most severe periods of drought in 1982 and 1988 ( Figure 5). To capture this relationship, a variable (DUR) was created that measured the length of time (in months) that the PDSI had been increasing or decreasing. If the trend changed from increasing or decreasing or vice versa, DUR was set to 0. If PDSI was increasing (decreasing), 1 (-1) was added to the previous month's DUR. DUR should represent the length of the trend in PDSI, but does not account for the high initial PDSI that would stimulate growth. DURMAX was set equal to DUR times the maximum PDSI of the preceding 12 months to quantify the potential fine-fuel growth. Correlating these variables, lagged from 1 to 12 months, with spring AB and fall AB provided the link between PDSI and AB. Fall AB was most highly correlated with the previous June's DUR (r = -0.34, p < 0.1212, n = 22) and DUR from 12 months earlier, November (r = -0.34, p < 0.1170, n = 22). Spring AB was most significantly associated with the DUR from the previous June (r = -0.50, p < 0.0160, n = 23), DURMAX of the previous October (r = -0.50, p < 0.0213, n = 21), and DURMAX of the previous June (r = -0.61, p < 0.0025, n = 22). Figure 7 shows that these correlations are due to one or two highly influential observations. We usually deleted these obvious outliers, but these extreme fire seasons are the events we are trying to predict. Table 1 lists the results of multiple linear regressions using the variables mentioned. Two independent variables explain 33 percent of the variation in the amount of AB during the fall fire season, and three variables explain 51 percent of the variation of AB in the spring fire season. Table 1. Regression coefficients and measures of fit for area burned (AB), in hectares, for spring and fall fire seasons. Variable Estimate Std. error of estimate Significance -------- -------- ---------------------- ------------ Fall AB: Adj. R-sq = 0.370, F = 6.002, p < 0.012, n = 18 Intercept 464.214 555.388 0.416 SOIAV Jun -1571.743 551.436 0.012 DUR Nov 495.446 346.032 0.173 Spring AB: Adj. R-sq = 0.512, F = 7.647, p < 0.002, n = 20 Intercept 459.947 91.134 0.000 SOIAV Jun -204.159 92.762 0.043 DURMAX Oct -18.655 10.649 0.099 DURMAX Jun -69.369 24.513 0.012 CONCLUSIONS Although fires in the Ohio River Valley almost always are started by humans, there is a relationship among severe fire seasons, climate, and the SOI. These connections offer some potential for predicting the severity of a fire season. A larger data set is needed to further test these connections, and we soon will be gaining access to the large USDA Forest Service fire data base and other data bases at the state level. Weather patterns and systems that develop and generally move from west to east are influenced by the Southern Oscillation and the earth's surface (mountains, bodies of water, soil moisture, vegetation, and so on). Weather and climate in the Ohio River Valley is strongly influenced by this, the Atlantic Ocean, and Gulf moisture. We hope to acquire a data base of the North Atlantic Oscillation to enhance our ability to predict droughts and fire seasons. ACKNOWLEDGEMENTS We would like to thank Kathleen Kennedy, Ginger Brudevold, Warren Heilman, and Thomas Swetnam for providing data and insight for this project. LITERATURE CITED Palmer, W. C. 1965. Meteorological drought. Res. Pap. No. 65, U.S. Weather Bureau, Office of Climatology, Washington, DC. 58 p. Simard, A. J. and W. A. Main. 1987. Global climate change, the potential for changes in wildland fire activity in the southeast. Proceedings of the symposium on climate change in the southern United States: future impacts and present policy issues, May 28-29, 1987, New Orleans, LA. U.S. Environmental Protection Agency, Washington, DC. pp. 280-308. Simard, A. J., D. A. Haines and W. A. Main. 1985. Relations between El Nino/ Southern Oscillation anomalies and wildland fire activity in the United States. Agricultural and Forest Meteorology 36:93-104. Swetnam, T. W. and J. L. Betancourt. 1990. El Nino-Southern Oscillation (ENSO) phenomena and forest fires in the southwestern United States. In: Proceedings of the sixth annual Pacific climate (PACLIM) workshop, March 5-8, 1989, Asilomar, CA. Edited by J. L. Betancourt and A. M. MacKay. California Department of Water Resources, Interagency Ecological Studies Program Technical Report 23. pp. 129-134. Swetnam, T. W. and J. L. Betancourt. 1992. Temporal patterns of El Nino/ Southern Oscillation -- wildfire teleconnections in the southwestern United States. In: El Nino -- historical and paleoclimatic aspects of the Southern Oscillation. Edited by H. F. Diaz and V. Markgraf. Cambridge Univ. Press. pp. 259-270. Wright, P. B. 1989. Homogenized long-period southern-oscillation indices. International Journal of Climatology 9:33-54.