Abstract

Precipitation is among the more limiting meteorological factors affecting the occurrence and extent of forest fire. We examined the correlation between burned area of individual wildfires and the rainfall amounts occurring on the day of the burn and the number of consecutive dry days for a range of limiting daily rainfall amounts (0–6mm) used to define a “dry” day. Daily threshold rainfall levels that most significantly affected area burned were determined for each ecoregion in Yunnan province, a major fire-prone area, in southwestern China. Results showed that the burned area of a wildfire decreased exponentially with increasing rainfall amounts on the day of burning. Burned area was also positively correlated to the number of consecutive dry days prior to burning. The threshold rainfall value providing the highest correlation between burned area and the number of consecutive dry days prior to a burn varied between ecoregions. Consecutive dry days with rainfall less than the specified threshold predominantly affected large fires (>100 ha) rather than more frequently occurring small fires. These results will help forest managers evaluate regionalfire danger indices for forest fire prevention, particularly for catastrophic forest wildfires causing significant economic losses and threats to human life and environment.

1. Introduction

A striking increase in the number and burned area of wildfires has been reported during the past several decades in many regions of the world [14]. The number and extent of wildfires in a region are driven by many factors including weather conditions [5], human activities [6], fuel characteristics [79], fire management activities [10, 11], land use changes [12], and climatic change [1315]. Of these factors, climatic variability has been considered to be one of the major determinants of the occurrence and burned area of wildfires, particularly large fires. The effect of climate variability on wildfire occurrence and burned area remains poorly understood and difficult to predict especially under various scenarios of global climate change.

On decadal scales, climatic variability regulates the compositional and structural characteristics of the vegetation such as species populations and their drought tolerance [16] and biomass (fuel) continuity [17] and thus affecting fire regime. On interannual and shorter temporal scales, climate variability affects the flammability of live and dead forest vegetation which may change wildfire behavior [18, 19]. Previous studies have demonstrated that large-scale climate patterns related to the Pacific decadal oscillation and intense El Niño/La Niña events were more closely linked with the frequency and extent of wildfires [20, 21]. Fire severities in a given fire season were reported to be related to ocean surface temperature anomalies [22] and the Palmer drought severity index (PDSI) [5, 23]. Increases of fire occurrences in Northern Rockiesforests have been shown to be strongly associated with increased spring and summer temperatures and an earlier spring snowmelt [24]. An increase in burned area across Canada from the 1920s to 1990s was correlated with regional warming trends [25]. Both satellite dataand the Russian Federal Forest Service data showed an increase in the area burned in Russia coincided with the current warming trends and extended fire seasons, as well as with predictions of fire regime change [26]. Valendik [27] and Lelyakin et al. [28] also suggested that forest fire was expected to occur more frequently in boreal Asiaas a result of increased mean temperature.

Regional climate patterns and local weather conditions clearly play a critical role in setting conditions favorable for forest fire occurrences, and their effect on a particular fire can be inferred from antecedent climate and weather conditions at the time when a fire occurs [29]. Prolonged periods with high temperatures and no rain correlate well with periods of high forest fire danger. Prevailing dry conditions particularly affect the flammability of dead forest fuels. In the Daxinganling forest region of China, forest fires ignited by lighting mostly occur under conditions of dry thunderstorms producing little precipitation, since precipitation levels <1 mm have almost no effect on fuel moisture content [30]. The longer the period of consecutive dry days (days without rainfall or with rainfall less than some threshold value), the drier the fuel and the more readily the fuels burn [31]. A period of 3 to 5 consecutive dry days can greatly increase ignition probabilities, 7 to 10 consecutive dry days markedly increases fire spread probabilities, and the risk of catastrophic forest fire occurrence increases greatly with >15 consecutive dry days [32].

The likelihood of large fires, therefore, increases during extreme weather conditions (e.g., a long dry period). Venäläinen et al. [33] estimated that at any location in Finland approximately once every ten years a 40-day period with at most 10mm of accumulated rain will occur, and these extreme dry conditions will triple the number of fire occurrences including large fires during a fire season. Prolonged dry weather conditions (about 40 days without precipitation) can sufficiently dry live fuels and larger dead fuels to the point of carrying large, intense fires once they are ignited in the subalpine forest types. Variation in area burned per fire has been shown to be highly correlated with the moisture content of 100-hour (2.5–7.6 cm diameter) and 1000-hour (>7.6 cm diameter) dead fuels [34]. Once fuels reached critically dry moisture levels late in the season, the spatial pattern of the large, severe stand replacing fires is controlled by weather (wind direction and velocity), not by fuels, stand age, or firefighting activities [35, 36]. Variation in fuel abundance and topography (including formidable barriers such as the canyon, rivers, and roads) had little influence on the severity or direction of the fire when fuel moistures were critically low [34]. As global warming and climate change related drought events become more frequent and severe regionally and globally, the frequency and burned area of forest fires, especially large, severe, fires tend to increase correspondently [1].

Yunnan province is a major forested area with complex, rich forest resources in the southwest China. The province has experienced a large number of wildfires that caused enormous losses in terms of forest resources, human lives, and economic disruptions [37]. The death toll averages over one hundred persons per year and the average economic loss directly related to forest wildfires reached over ten million US dollars annually during the past decade [37]. Due to extensive areas dominated by flammable forest species, forest fires causing serious damage occur frequently, and fire prevention is becoming quite difficult [38]. Preliminary analyses of historical fire and weather/climate data indicate the daily precipitation and distribution throughout the fire season are an important determinant of forest fire occurrence, particularly the occurrence of large, damaging fires [39, 40]. The major objective of this study, therefore, is to quantify the relationship between precipitation and burned area of individual fires, and how this relationship may vary by ecoregion, based on historic fire data and observed precipitation data from local weather stations. We hypothesize that the size of a fire occurring on any given day will be correlated with the amount of rainfall on that day. We also predict the number of consecutive dry days with rainfall below a threshold value will be correlated with the burned area of large fires. Specifically, we will test (1) whether the distribution of area burned is different between days where precipitation occurs and dry days, (2) how area burned changes with rainfall on the days prior to fire occurrence, and (3) how the threshold precipitation values best correlated with area burned varies by ecoregion. Answers to these questions will have significant implications in fire severity forecasting and fire prevention in the region.

2. Materials and Methods

2.1. Study Area

Yunnan province is located in the southwestern area of China (21°9′–29°15′ N, 97°30′–106°E, Figure 1). The topography of the province is diverse and includes flat valley bottoms and steep mountains. General elevation ranges from approximately 4000 m in the northwest to less than 1000 m in the southeast. The region has a plateau monsoon climate characterized by moderately hot, humid summers and warm, dry winters. Mean annual air temperature averages 14–16°C, with mean temperatures of 6–8°C in January and 19–22°C in July. Winter prevails dry continental monsoon, while summer prevails wet marine monsoon, so the distribution of precipitation in the season is very uneven. Annual rainfall averages 1100 mm, and dry season extends from November to April, accounting for only 10–20% of annual precipitation. Rainy season extends from May to October, and most precipitation is 6–8 for three months, accounting for about 60% of annual precipitation (Figure 2). From north to south there is a gradient of increasing temperature, precipitation, and humidity (Table 1).

Due to the influence of tropical monsoons from the south and the decrease in altitude from north to south, Yunnan has highly diverse vegetation types ranging from alpine meadows and montane and subalpine temperate forests at higher altitudes to subtropical forests and tropical rainforests at lower altitudes. Broad-leaved species dominate forests in the southwestern portion of the province, while most other areas are dominated by coniferous forests which cover an area of approximately 4.53 million hectares and account for 48.6% of the province's forested area. The main vegetation types that are prone to forest fires are forests and woodlands dominated by Yunnan pine (Pinus yunnanensis Franch). These natural coniferous forests, along with plantations planted in recent years, are highly flammable and susceptible to varying fire ignition sources.

The study area includes five ecoregions, reflecting the spatial heterogeneity in climate, vegetation types, and topography (Table 1). Mountainous ecoregions (VI and V) have a cold and dry climate, while the lower elevation ecoregions (I, II, and III) are characterized by warmer, more humid climates. The higher-elevation ecoregions to the north have higher proportions of coniferous forests, while the lower-elevation ecoregions in the south have higher proportions of evergreen broad-leaved forest.

Ecoregion I has a tropical monsoon climate with higher temperature and abundant precipitation. The ecoregion is frost-free throughout the year; and the area has distinct wet and dry seasons, with precipitation in rainy season (May to October) accounting for 80%–90% of total annual precipitation. This ecoregion has a higher proportion of forest than the other ecoregions, and the vegetation types are tropical seasonal rainforest, evergreen monsoon forest, and mountain rainforest dominated by Pometia pinnata J. R. et G. Frost., Parashorea chinensis Wang Hsie, and Lannea coromandelica (Houtt.) Merr.

Ecoregions II and III have similar landform and climate. They have a plateau and subtropical climate, with mildseasonal temperatures and moderate precipitation. There is, however, an obvious difference in vegetation types. Ecoregion II is covered by monsoon evergreen broad-leaved forests dominated by Schima wallichii Choisy, Castanopsis delavayi Franch, and small areas of Pinus kesiya Royle ex Gordon forest, while ecoregion III is dominated by evergreen Quercus Linn. forests and large areas of Pinus yunnanensis Franchmixed with shrubs and grasslands.

Ecoregion IV has a plateau and subtropical climate, but with lower temperatures and less precipitation than ecoregions II and III. The area also has a lower percent forest cover which is dominated by evergreen Quercus semecarpifolia Smith, Cyclobalanopsis glauca (Thunb.) Oerst forest, and some Pinus yunnanensis Franch forest. Ecoregion V is highest in elevation averaging approximately 4000 meters. The area has a low temperature with a winter snow line elevation of around 5000 meters. Annual precipitation averages 600–700 mm, with dry season rainfall accounting for 10–20% of annual precipitation. Vegetation types include cold-temperate needle-leaf forest dominated by Picea asperata Mast., Abies fabri (Mast.) Craib, and Pinus densata Mast.intermixed with alpine meadow.

2.2. Climate and Fire Data

The fire season in Yunnan province extends from December to May with wildfires mainly occurring from mid-February to mid-May (Figure 2). Human-caused surface fires account for over 95% of total fires. Crown fires are rare and mostly occurred in young coniferous forests. From 1996 to 2008, in total 7,926 forest fires (135,788 ha in total) were reported in Yunnan. For each fire, the date, starting and ending times, location, burned area, and vegetation type were recorded by the Forest Fire Prevention Headquarter Office of Yunnan Province. Daily rainfall data during the same period were obtained from 35 meteorological stations of the China meteorological data sharing network (Figure 1). The daily rainfall data for the area were generated by using Daymet [41], which uses regression and digital elevation maps to interpolate data from existing weather stations over complex terrain. In this study, only 1,983 fires occurring in the fire season within the prefectures with observed rainfall data were used for data analysis.

2.3. Data Analysis

To evaluate the effect of precipitation on the distribution of fire size (area burned), forest fires were classified into two groups: with precipitation on the day of the fire () and without precipitation () on the day of the fire. Since the distribution of burned area takes a negative exponential trend, the log-transformed burned area was used in analysis of variance to compare the difference in mean fire size between these two groups. With the assumption that precipitation may significantly reduce the size of large, catastrophic fires, the ratio of fire sizes at varying percentiles of the two groups (distributions) were computed to analyze the distributional difference in fire size between the two groups.

To further quantify the change of burned area by precipitation levels on a fire day, quantile regression was used to analyze the relationship between precipitation and burned area of a forest fire. Quantile regression has several advantages over traditional (mean based) regression when the response variable (area burned) varies greatly with unequal variance and/or an abnormal distribution by the predictor variable (amount of rainfall here) and when extreme values represented by the conditional quantiles are of interest. Based on the scatter plot of area burned versus rainfall, the nonlinear quintile regression model was specified as where is the estimated conditional Q-percentile of area burned given a value of (rainfall) and and are the estimated regression coefficients. Limited by the number of forest fires by ecoregion, a pooled model was constructed for the whole study area based on the algorithm described in Koenker and Park [42].

With daily precipitation data for the period prior to a fire, the number of consecutive dry days was calculated for each of the 1,983 forest fires, with “dry day” being defined as daily rainfall , or 6 mm. Spearman’s rank correlation analysis was applied to compute the correlation between the number of consecutive dry days (daily rainfall less than a specific value) and the log-transformed burned area of a subsequent forest fire for each ecoregion. The value of daily precipitation defining a dry day which maximized Spearman’s rank correlation coefficient and was significantly different from zero based on bootstrapping analysis was specified as the critical value (threshold) influencing burned areas. Furthermore, forest fires within an ecoregion were divided into six size (burned area, in ha) classes: <1.0 (small fire), 1.0–9.9 (moderately small fire), 10.0–49.9 (moderate fire), 50.0–99.9 (moderately large fire), 100.0–500.0 (large fire), and 500–1000 (extremely large fire), and Spearman’s rank correlation coefficient was computed between the number of consecutive days with daily precipitation ≤ the threshold value and the log-transformed burned areas by fire size class and ecoregion. The effect of threshold precipitation on forest fires in different size classes within an ecoregion may be evaluated based on the change of Spearman’s rank correlation coefficients.

3. Results

3.1. The Relationship between Burned Area and Rainfall Amount

Area burned followed a negative exponential distribution for forest fires occurring on both rainy days (with precipitation) and dry days (without precipitation) with mean fire size (area burned) of 43.5 ha and 76.3 ha, respectively, but the distribution for the dry day fires was more skewed, having a longer tail (larger, catastrophic fires) (Table 2). Further examination of the ratio of the percentiles of burned area on rainy days versus on dry days indicated that the burned area of 90% of the fires occurring on rainy days mostly was 0.5 to 0.8 times large as those occurring on dry days, and among the largest fires (largest 10%) the burned area on rainy days was less than one tenth large as those occurring on dry days (Figure 3(a), Table 2). ANOVA of log-transformed area burned indicated a significant difference between a dry day and a rainy day fire size (); however, when the largest 10% of fires were deleted from both samples, the difference was no longer significant at . This suggests that precipitation on the same day as a fire occurrence overall has potential to greatly reduce the infrequent, large, catastrophic fires but has a much smaller impact on the size of most frequent, small fires. Quantile regression further quantified the changes in the percentile of areas burned with precipitation on the same day as a fire occurrence (Figure 3(b), Table 3). Based on model statistics, the negative exponential model did a good job of quantifying the effect of precipitation on area burned. Precipitation alone accounted for at least 33% of variation in areas burned depending on the percentile specified (Table 3).

3.2. Correlation between Burned Area and the Number of Consecutive Dry Days prior to a Fire

Spearman’s rank correlation analysis suggests that the burned area of a fire was positively correlated with the number of consecutive days with daily precipitation less than a specified threshold value prior to the day of the fire for all ecoregions (Figure 4). Correlation coefficients generally displayed a unimodal distribution which peaked at an ecoregion-specific precipitation value. In ecoregions I and II, the highest correlation (the peak) was only found at 1 mm (, ; , ), the peak in ecoregion III was around 3 mm (, ), and in ecoregion IV and ecoregion V, the peak was located at a precipitation level of around 5 mm (, ; , ) (Figure 4). This suggested that there is a variation in the daily threshold rainfall in defining dry day among five ecoregions due to differences in some aspects such as vegetation and climate by ecoregion, with 1 mm for ecoregions I and II, 3 mm for ecoregion III, and 5 mm for ecoregions IV and V. The daily threshold rainfall in defining dry day has a practical significance for the local forest fire prevention. The number of consecutive dry days can provide fire control managers with reference for measuring fire danger severity which may be useful in planning fire control operations.

3.3. The Impact of Consecutive Dry Days on Different Sizes of Fire

Given the precipitation values (or thresholds) corresponding to the maximum correlation coefficient in each ecoregion, the rank correlation between the number of consecutive days prior to a fire with daily precipitation less than the threshold and burned area of a fire was analyzed for six fire size groups specified in the method, respectively. Within each ecoregion, the correlation coefficient increased rapidly from the small fire group (<1.0 ha) to the moderately large fire group (50.0–99.9 ha) and then maintained a gentle increase across the larger fire size groups (Figure 5). The correlation coefficients for fire sizes <50 ha were all less than 0.2 and not significantly different from zero at , while the correlation was statistically significant for all fire size groups ≥50 ha. This suggests that dry periods prior to burning have a significant influence on the size of large fires (≥ 50 ha) or may more likely cause the increase in fire size into the larger size categories, but not much of an effect on the size of small or moderate fires.

4. Discussion

Variability in intensity and timing of precipitation (e.g., daily rainfall amount and the number of consecutive days without rain or with rainfall amounts less than certain thresholds) strongly influence the extent (area burned) of wildfires [39, 43]. Precipitation generally reduces the likelihood of large, catastrophic fires (Table 2). In Yunnan province, for instance, the area burned of the largest fire that occurred on a fire day without precipitation is more than ten times larger than that of a fire on a day with precipitation prior to the fire occurrence. D. X. Viegas and M. T. Viegas [43] reported an exponential relationship between areas burned and total precipitation during fire season. One explanation for the relationship is, presumably, that precipitation increases fuel moisture content thus decrease the potential for fire spread [44, 45].

Song et al. [46] used a fire model to simulate the effect of different rainfall probabilities on the frequency-size distribution of forest fires and found that rainfall had a larger impact on large forest fires than on small ones. This is consistent with our study using real fire and precipitation data. As we have shown, area burned is, on average, significantly reduced with an increase in precipitation (Figure 3(b)). As the percentile of area burned increased from 5% to 95%, the percentile regression lines became steeper, meaning precipitation had a substantial impact on the size (or extent) of larger fires. Based on the fitted equations, one can calculate that the area burned of fires in the 99th percentile would be reduced from around 852 ha (the coefficient in Table 3) at a very small amount of precipitation (e.g., 0.01 mm) to 298 ha at 6 mm of precipitation, a reduction of nearly two-thirds. Similarly, the area burned for the 95th percentile of large fires would be reduced from an average of 57 ha to 51 ha, a decrease of approximately 10% for the same range of variability in precipitation.

The coefficient reflects the th percentile of the area burned when the precipitation approaches zero (no rain), and the coefficient reflects the rate of reduction in the area burned as the precipitation changes. The larger the is, the lower the rate is. From the regression models in Table 3, precipitation on a fire day is especially effective to reduce the likelihood of large fires (100–499 ha) and extremely large fires (≥500 ha) but has a limited effect on reducing the average size of moderate and small fires (<100 ha). Knapp [47], Littell et al. [48], and others have indicated the potential impact of other climatic factors such as temperature and vegetation cover on area burned. Because we lacked the data for constructing more comprehensive quantile regression models, we built a set of models for the entire province. However we recognize that, because of the great variability in temperatures and vegetation composition among ecoregions, distinct ecoregion models should be developed to better quantify the precipitation effect. This would provide better information for decision making on presuppression activities.

Fire occurrence is closely related to meteorological preconditions [49, 50]. In Yunnan province, the driest months (November and December through February, with mean daily rainfall amounts <1 mm) are followed by the peak fire months (March through May), suggesting a lag effect for periods of little precipitation on fire incidence and extent (Figure 2). The mechanism for this lag effect may be related to the interception storage capacities of forest canopies, ground evaporation, and gradual changes in soil and fuel moisture conditions. Rutter [51] reported interception storage capacities of 1.0 to 2.1 mm for coniferous forests and 1.0 mm for deciduous forest in summer. A similar value of 1.3 mm was found by Simpson et al. [52] for a Douglas-fir forest in western Washington. Spittlehouse and Black [53] used 3.0 mm for a Douglas-fir forest in northwest of Courtenay on Vancouver Island, but this value includes some loss due to evaporation. Flannigan and Harrington [39] defined limiting precipitation of 1.5 mm as a dry day, assuming this value closely to fit the estimates of maximum water holding capacity of various forest canopies. Hall [54] estimated a precipitation threshold of 2 mm for both shrubs and trees and a value of 1 mm in grass vegetation, based on an analysis of the difference in precipitation amounts between wildfires associated with lightning and lightning without ignitions. In national fire danger rating system (NERDS) used in the United States, a drought index is calculated based on the number of consecutive days with precipitation ≤5.0 mm [44]. In the Canadian fire danger rating system (CFDRS), the definition of the threshold precipitation for drought was 2.8 mm [55].

In this study, predryness conditions were estimated by the number of the consecutive dry days with daily rainfall amounts across a range of values (<0, 1, 2, 3, 4, 5, and 6 mm). These measures of “dryness” were correlated, by ecoregion, with area burned of forest fires. The results suggested that there is variation in the threshold rainfall (1 mm for ecoregions I and II, 3 mm for ecoregion III, and 5 mm for ecoregions IV and V) that was maximally correlated with the area burned. There is a decreasing trend for the threshold rainfall value with an increase in humidity across ecoregions, most likely due to differences in soil moisture content determined by wet conditions across ecoregions (not reported here).

Drought influences forest fire by means of changing fuel moisture content [56]. In arid regions, with perennial drought, lower soil water content typically means that surface fuels are relatively dry. There tends to be a greater difference between fuel moisture content and the moisture of extinction, and subsequently the threshold rainfall value is larger. Conversely, in humid regions, due to more precipitation, soils have consistently higher moisture contents and surface fuel moistures are higher. Consequently, the amount of precipitation needed to raise fuel moisture contents to the moisture of extinction is less and the threshold rainfall value is correspondingly smaller.

This study showed that the number of consecutive days has a greater association with the size of large fires (>100 ha) but has less influence on the size of small fires. This may be explained by the fact that most fires in the study area are caused by humans. Small fires have a greater degree of randomness, largely influenced by human factors, so the correlation between small fires and preconditions of dryness is weaker. Concerning the relationship between environmental conditions and major forest fires, Zhang et al. [57] found a close relationship between the size of major forest fires and consecutive dry days prior to fire occurrence. Brotak and Reifsnyder [58] and Crimmins [59] demonstrated that extended dryness is associated with large fire occurrence, which is often dominated by synoptic weather conditions of dry. Many other studies had shown previous drought to have a significant impact on large fires occurrence, and forecast the probability of large fire events using drought index [6062]. Besides, the dryness of fuel is not only associated with rainfall but also with other meteorological factors and fuel characteristics [5, 7, 9]. In our study, differences in fuel types among ecoregions will inevitably affect the threshold rainfall of burned areas; therefore, whether the conclusion of this study can be applied to other areas requires further research.

5. Conclusions

Based on our analysis of the 1996–2008 wildfire records and precipitation data in Yunnan province, southwestern China, the following conclusions were drawn.(1)The size of a fire is closely associated with the rainfall amounts on the day of the burn. The relationship between burned area and rainfall amounts displayed a significant exponential decay. Large amounts of rainfall have the potential to reduce fire size with increases in fuel moisture contents and slow rates of spread.(2)Burned area is related to the number of consecutive dry days prior to fire occurrence. The strength of this correlation varies depending on the daily rainfall amount used to define “dry day,” reaching a maximum value at a threshold daily rainfall amount which changes between ecoregions. This suggests that different thresholds of daily rainfall amount should be used to define “dry day” for the ecoregions in Yunnan province.(3)The increasing trend in correlation between burned area (fire size) and the number of consecutive dry days defined by the identified threshold of daily rainfall amount from the small fire group to large fire group manifests that the threshold of daily rainfall identified for each ecoregion will be useful for presuppression activities of large, severe fires (e.g., > 100 ha). One application will be, among others, evaluating the current fire danger rating system in Yunnan province.

Conflict of Interests

The authors declare no conflict of interests.

Authors’ Contribution

Feng Chen, Zhaofei Fan, and Shukui Niu contributed equally to this work.

Acknowledgments

The authors would like to acknowledge Forest Fire Prevention Headquarter Office of Yunnan Province for providing fire data. Funding for this research was provided by the Forestry Public Welfare Project (200804001), which was initiated by SFA (State Forestry Administration). Dr. Scott D. Roberts from the Department of Forestry, Mississippi State University, and Dr. Michael Crosby from Shorter University reviewed this paper. The authors thank both of them.

Supplementary Materials

The table presents the code, name, and geographical location of 36 meteorological stations used for this study in Yunnan province, with one station in one prefecture. The observation time of these stations is from 1960 to 2008.

  1. Supplementary Table