Risk assessment of forest fires requires an integrated estimation of fire occurrence probability and burn probability because fire spread is largely influenced by ignition locations as well as fuels, weather, topography and other environmental factors. This study aims to assess forest fire risk over a large forested landscape using both fire occurrence and burn probabilities. First, we use a spatial point processing method to generate a fire occurrence probability surface. We then perform a Monte Carlo fire spread simulation using multiple fire ignition points generated from the fire occurrence surface to compute burn probability across the landscape. Potential loss per land parcel due to forest fire is assessed as the combination of burn probability and government-appraised property values. We applied our methodology to the municipal boundary of Gyeongju in the Republic of Korea. The results show that the density of fire occurrence is positively associated with low elevation, moderate slope, coniferous land cover, distance to roads, high density of tombs and interaction among fire ignition locations. A correlation analysis among fire occurrence probability, burn probability, land property value and potential value loss indicates that fire risk in the study landscape is largely associated with the spatial pattern of burn probability.

The United States national inventory program measures a subset of tree heights in each plot in the Pacific Northwest. Unmeasured tree heights are predicted by adding the difference between modeled tree heights at two measurements to the height observed at the first measurement. This study compared different approaches for directly modeling 10-year height increment of red alder (RA) and ponderosa pine (PP) in Washington and Oregon using national inventory data from 2001–2015. In addition to the current approach, five models were implemented: nonlinear exponential, log-transformed linear, gamma, quasi-Poisson, and zero-inflated Poisson models using both tree-level (e.g., height, diameter at breast height, and compacted crown ratio) and plot-level (e.g., basal area, elevation, and slope) measurements as predictor variables. To account for negative height increment observations in the modeling process, a constant was added to shift all response values to greater than zero (log-transformed linear and gamma models), the negative increment was set to zero (quasi-Poisson and zero-inflated Poisson models), or a nonlinear model, which allows negative observations, was used. Random plot effects were included to account for the hierarchical data structure of the inventory data. Predictive model performance was examined through cross-validation. Among the implemented models, the gamma model performed best for both species, showing the smallest root mean square error (RSME) of 2.61 and 1.33 m for RA and PP, respectively (current method: RA—3.33 m, PP—1.40 m). Among the models that did not add the constant to the response, the quasi-Poisson model exhibited the smallest RMSE of 2.74 and 1.38 m for RA and PP, respectively. Our study showed that the prediction of tree height increment in Oregon and Washington can be improved by accounting for the negative and zero height increment values that are present in inventory data, and by including random plot effects in the models.

Forest wildfires consume and redistribute carbon within forest carbon pools. Because the incidence of wildfires is unpredictable, quantifying wildfire effects is challenging due to the lack of prefire data or controls from experiments over a large landscape. We explored a quasi-experimental method, propensity score matching, to estimate wildfire effects on aboveground forest woody carbon mass in Washington and Oregon, United States. Observational data, including national forest inventory plot measurements and satellite imagery metrics, were utilized to obtain a control set of unburned plots that are comparable to burned plots in terms of environmental conditions as well as spatial locations. Three matching methods were implemented: propensity score matching (PSM), spatial matching (SM), and distance-adjusted propensity score matching (DAPSM). We investigated if propensity score matching with and without spatial adjustment led to different outcomes in terms of (1) balance in covariate distributions between burned and control plots, (2) mean carbon mass obtained from the selected control plots compared to burned and all unburned plots, and (3) estimates of wildfire effects by burn severity. We found that PSM and SM, which use only the environmental covariate set or the spatial distance for estimating propensity scores, respectively, did not appear to produce a comparable set of control plots in terms of the estimated propensity scores and the outcomes of mean carbon mass. DAPSM was the preferred method both in balancing the observed covariates and in dealing with unobservable confounding variables through spatial adjustment. The average wildfire effects estimated by DAPSM showed clear evidence of redistribution of carbon among aboveground woody pools, from live to dead trees, but the consumption of total woody carbon by wildfire was not substantial. Only moderate burn severity led to significant reduction of total woody carbon mass across Washington and Oregon forests (64% of control plots remained on average). This study provides an applied example of a quasi-experimental approach to quantify the effects of a natural disturbance for which experimental settings are unavailable. The study results suggest that incorporating spatial information in addition to environmental covariates would yield a comparable set of control plots for wildfire effects quantification.

No abstract