Bandwidth determination for kernel density analysis of wildfire events at forest sub-district scale

Forest fire is regarded as one of the most significant factors leading to land degradation. While evaluating fire hazard or producing fire risk zone maps, quantitative analyses using historic fire data is often required, and during all these modeling and multi-criteria analysis processes, the fire event itself is taken as the dependent variable. However, there are two main problematic issues in analyzing historic fire data. The first difficulty arises from the fact that it is in point format, whereas a continuous surface is frequently needed for statistically analyzing the relationship of fire events with other factors, such as anthropogenic, topographic and climatic conditions. Another, and probably the most bothersome challenge is to overcome inaccuracy inherent in historic fire data in point format, since the exact coordinates of ignition points are mostly unknown. In this study, kernel density mapping, a widely used method for converting discrete point data into a continuous raster surface, was used to map the historic fire data in Mumcular Forest Sub-district in Muğla, Turkey. The historic fire data was transferred onto the digital forest stand map of the study area, where the exact locations of ignition points are unknown; however, the exact number of ignition points in each compartment of the forest stand map is known. Different random distributions of ignition points were produced, and for each random distribution, kernel density maps were produced by applying two distinct kernel functions with several smoothing parameter options. The obtained maps were compared through correlation analysis in order to illustrate the effect of randomness, choice of kernel function and smoothing parameter. The proposed method gives a range of values rather than a single bandwidth value; however, it provides a more reliable way than comparing the maps with different bandwidths subjectively by eye.

• Wildfire density is mapped by kernel density method at forest sub-district scale. • Three random distributions of ignition points are compared by correlation analysis. • Model functions are fitted onto correlation values versus bandwidth data. • Plot of model function derivative versus bandwidth is used to find final bandwidth. • Proposed method is more reliable than choosing bandwidth value subjectively by eye.