Next Article in Journal
Cost Effective and Sustainable Test Methods to Investigate Vulnerabilities of EMP Attack on Existing Buildings Using Public Frequency Devices
Previous Article in Journal
A Comparison Study between Universal Design for Learning-Multiple Intelligence (UDL-MI) Oriented STEM Program and Traditional STEM Program for Inclusive Education
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Sustainability
Volume 13
Issue 2
10.3390/su13020569
sustainability-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Communication

Predicting Post-Fire Tree Mortality in a Temperate Pine Forest, Korea

by
Semyung Kwon
1,
Sanghyun Kim
1,2,3,
Jeonghwan Kim
1,4,
Wonseok Kang
1,*,
Ki-Hyung Park
1,
Chan-Beom Kim
5 and
Miguel M. Girona
2,3,6
1
Division of Forest Restoration and Resource Management, National Institute of Forest Science, Seoul 02455, Korea
2
Group of Research in Ecology-MRC Abitibi (GREMA), Forest Research Institute, University of Québec in Abitibi-Témiscamingue, Amos Campus, Amos, QC J9T 2L8, Canada
3
Center for Forest Research, University of Québec in Montréal, Montréal, QC H3C 3P8, Canada
4
Division of Forest Welfare Statistics, Korea Forest Welfare Institute, Daejeon 35236, Korea
5
Temperate Middle Part Plant Conservation Team, Sejong National Arboretum, Sejong 30106, Korea
6
Restoration Ecology Group, Department of Wildlife, Fish, and Environmental Studies, Swedish University of Agricultural Sciences (SLU), 901 83 Umeå, Sweden
*
Author to whom correspondence should be addressed.
Sustainability 2021, 13(2), 569; https://doi.org/10.3390/su13020569
Submission received: 16 September 2020 / Revised: 14 December 2020 / Accepted: 18 December 2020 / Published: 9 January 2021
(This article belongs to the Section Sustainability, Biodiversity and Conservation)
Browse Figures
Versions Notes

Abstract

:
Warmer and drier conditions in temperate regions are increasing the length of the wildfire season. Given the greater fire frequency and extent of burned areas under climate warming, greater focus has been placed on predicting post-fire tree mortality as a crucial component of sustainable forest management. This study evaluates the potential of logistic regression models to predict post-fire tree mortality in Korean red pine (Pinus densiflora) stands, and we propose novel means of evaluating bark injury. In the Samcheok region of Korea, we measured topography (elevation, slope, and aspect), tree characteristics (tree/crown height and diameter at breast height (DBH)), and bark injuries (bark scorch height/proportion/index) at three sites subjected to a surface fire. We determined tree status (dead or live) over three years after the initial fire. The bark scorch index (BSI) produced the best univariate model, and by combining this index with the DBH produced the highest predictive capacity in multiple logistic regression models. A three-variable model (BSI, DBH, and slope) enhanced this predictive capacity to 87%. Our logistic regression analysis accurately predicted tree mortality three years post fire. Our three-variable model provides a useful and convenient decision-making tool for land managers to optimize salvage harvesting of post-fire stands.

1. Introduction

Agents of natural disturbance such as wind, fire, and insects affect forest dynamics and structure by altering ecosystem functioning, creating diversified habitats, and maintaining biodiversity [1,2,3,4]. Climate change will alter the natural disturbance regimes in forest ecosystems [1,5,6]. In recent decades, warmer and drier climates have increased the length of wildfire seasons and changed the patterns of wildfires [7,8,9,10]. Fire alters regeneration patterns and weakens plant defensive systems such as resistance to herbivores [11,12]. In addition, trees may be more vulnerable to fire-caused injury after climate-mediated drought stress and these consequences affect forest management planning [13,14].
Given these critical concerns, future post-fire landscape restoration will be challenging. This is particularly true in Korea, as Korean even-aged forests are characterized by a dense growing stock (approximately 154 m3 ha−1), leading to a high fuel accumulation. Thus, these forests are quite vulnerable to the expected increase in the frequency and severity of regional fire regimes [15]. Accordingly, predicting post-fire tree mortality for various tree species and forest types will be crucial for ensuring sustainable forest management in the face of climate change.
In the Republic of Korea, the majority of wildfires occur in the dry and windy spring. Between 1980 and 1999, 287 wildfires burned an average area of 1250 ha per year. More wildfires were observed over the last 20 years (2000–2019), where 481 wildfires burned an annual average area of 2650 ha [16], thereby demonstrating a marked increase in the occurrence and area burned because of more human activities in forests and climate change [17].
These wildfires have highlighted the need for reliable tools for predicting the post-fire mortality of Korean red pine (Pinus densiflora Sieb & Zucc.), a tree accounting for 28% of all trees in Korea [18]. Most previous studies have focused mainly on high intensity wildfires in productive forests [19] and the resulting substantial tree mortality [20]. This mortality can be immediate owing to the combustion of living tissue during the fire, or it can be delayed, with the tree slowly dying over a few years because of fire damage to the crown, trunk, and roots and, sometime later, arthropod attack [20]. Research on tree mortality following less intense fires, such as surface fires, remains less common, albeit a highly relevant subject for adapting forest management to climate change.
Surface fires alter the structure and function of forest ecosystems affecting the soil environment and causing the death of trees [21]. Tree mortality is a complex phenomenon influenced by multiple factors, including tree characteristics, climate, topography, and disturbance regimes [2]. Post-fire tree mortality is related to multiple direct and indirect fire effects [22]. Direct fire effects involve heat transfer and the resulting tissue necrosis, whereas indirect fire effects include altered physiology, insect damage, and pathogenic infection [23]. Despite their major economic and ecological consequences for forest ecosystems, these effects remain poorly understood and their consequences are difficult to predict.
The fire-induced cambium/phloem necrosis and xylem damage are the main causes of post-fire tree mortality [23]. Numerous methodologies have been applied to estimate cambial mortality, from the direct sampling of the cambial tissue to indirect measurements, such as the amount or height of crown/bark scorch or bark char severity [20,24,25]. Since measuring cambial activity directly is time-consuming, indirect alternatives have been developed. Bark injuries are commonly used as a substitute for direct cambium measurements because of the simpler approach of measuring bark injuries [20]; however, the correlation between bark injury and cambium status remains unreliable [20]. Therefore, it is necessary to develop accurate and cost-effective tools for predicting post-fire tree mortality.
This study describes tree mortality three years after a surface fire and identifies the main factors affecting within-stand mortality patterns. Simple and effective models are developed for predicting Korean red pine mortality after surface fires.

2. Materials and Methods

2.1. Study Area

Our study area (37°15′53″ N, 129°1′24″ E, 700–800 m asl) is located in the mountains around Samcheok, Korea, a region having stands dominated by Korean red pine (Figure 1). This region consists of pine forests in the crown layer, and Quercus sp. in the understory layer. Pine trees in our study stands originated from a plantation and were planted in 1972. In May 2017, the study area was affected by a 765 ha wildfire, which occurred in the minimally managed mountains around Dogye. Our study sites (1.4 ha) are in areas that experienced a surface fire.
The regional climate is temperate with warm and humid summers and cold, dry winters owing to the influence from the Siberian High. The mean annual temperature is 12.6 °C, and mean total annual precipitation is 1280 mm (Mountain Meteorology Observation System, Korea Forest Service 2019). The general wind direction in the study sites during the wildfire came from the south-western direction with a mean speed of 14 m/s (Figures S1 and S2). Soils are silt loam (SiL) and silt clay loam (SiCL).

2.2. Experimental Design and Sampling

We established three study sites after the May 2017 wildfires to monitor post-fire tree mortality. We selected sites subjected to a lower fire intensity and only 5% to 8% of the burned area was affected by surface fire. We measured all Korean red pine trees within each study site, although not those in the 10 m buffer zones surrounding each study site. The measured trees were all standing trees having intact crowns and were at least 8 m tall with a diameter at breast height (DBH) of >15 cm. A total of 413 trees were numerically tagged and surveyed four times over three years (11 July 2018, 19 October 2018, 20 September 2019, and 8 July 2020). We categorized trees as either dead or living based on crown state (intact, browning, or absence) for delayed mortality one year after the wildfire. All the studied trees had an intact crown at the beginning of this investigation. One year post the fire (2018), we assessed site topography (elevation, slope, and aspect), tree characteristics (DBH, tree height, crown height), and bark injury (bark scorch height and proportion). We selected these eight variables to calculate rates of post-fire tree mortality using simple-to-use nondestructive methods. For each tree, we measured DBH at 1.2 m above the ground on the uphill side of the tree. We recorded tree and crown height. We also measured slope and aspect using a digital hypsometer (Haglof Vertex III and Transponder T3, Sweden) and transformed them to eight cardinal directions for aspect with raster calculator. Elevation was measured with GPS (GPS Status & Toolbox, MobiWIA-EclipSim). We collected all topographical measurements on a tree’s uphill side perpendicular to the slope at a sufficient distance to obtain an accurate value. We divided each sampled tree trunk into quadrants (centered on N, S, W, E, respectively), and for each quadrant, we measured bark scorch height (BSH) and proportion of scorched wood (BSP, Figure 2a). Then, we defined the bark scorch index (BSI) as the sum of BSH multiplied by BSP for each quadrant (Figure 2b). In Figure 2a, the left photo study tree has a BSH of 1.9 and BSP of 0.8, whereas the study tree on the right photo has a BSH of 1.3 and a BSP of 0.6. BSI was developed to determine whether the tree would be dead or living after a surface fire. We expected that combining BSH with BSP would predict post-fire tree mortality more accurately.
To consider other potential variables influencing mortality patterns, we measured the following 11 soil chemical factors, one year after the fire: pH, total nitrogen (TN, %), available phosphorus (P, P2O5, ppm), organic matter 138 (OM, %); NaCl (%); electric conductivity (EC, dS·m⁻1); exchangeable cations (K+, Na+, Ca2+, Mg2+); and 139 cation exchange capacity (CEC, cmol+/kg) [26,27]. We also analyzed inter-tree competition based on DBH and inter-tree distance [28,29]. Competition data were taken one year after the wildfire in Korean red pine stands (n = 413). The distance (Distij) and DBH of each neighboring tree (j) within a 4 m radius of the subject tree (i) were measured to calculate the Hegyi’s competition index (CIi) as:
C I i = ∑ j = 1 n ( D B H i D B H j   ×   1 D i s t i j )

2.3. Data Analysis

We applied logistic regression to develop models differentiating living trees from dead trees by applying a binary response variable (dead/living). Therefore, all sampled Korean red pines were coded as either 0 (living) or 1 (dead) based on their condition during the three years of sampling after fire. We considered p-values ≤ 0.05 to be statistically significant. To overcome the limited sample size, we performed a preselection to test all measured variables. We developed models using only one explanatory variable at a time and included in the logit model all variables having nonsignificant but relatively lower p-values (p ≤ 0.25). Our models for Korean red pine included the eight above mentioned independent variables, i.e., elevation (m), slope (°), DBH (cm), tree height (m), crown height (m), BSP (%), BSH (m), and BSI except for aspect. We selected these variables as obtaining them was relatively quick, nondestructive, and provided accurate estimates. The logistic regression equation used to model mortality was the following:
P m = 1 1   +   e − ( β 1 +   β 2 X 2   +   ⋯   +   β k X k )  
where Pm corresponds to the probability of mortality, X 2 – X k are independent fixed variables and β 1 – β k are regression coefficients estimated from the mortality data.
We added the explanatory factors one at a time, applying a statistical method to decrease the Akaike information criterion (AIC). The AIC ranges from 0 to ∞, and lower values represent a better data fit. We used the eight independent variables to obtain the best fitting model (two topographical, three tree characteristics, and three bark injury variables). After conducting several tests, we selected BSI as the variable producing the best fit. We calculated delta AICs (∆AIC) for the logistic regression models and selected the best model with the smallest AIC value (same as ∆AIC = 0). We considered models to be comparable to the best model when their ∆AIC values were <2. We considered ∆AIC values of 3–9 as less significant and values of ≥10 as nonsignificant [30]. A Hosmer–Lemeshow test determined the goodness of fit for logistic regression and how well the data fit the model. We applied a log-likelihood function to compare the fit of different coefficients and to derive the maximum likelihood estimates of the model parameters. We calculated McFadden’s pseudo-R2 to examine the goodness of fit. McFadden’s R2 is comparable to the R2 obtained from linear regression models. Moreover, we evaluated all the models’ predictive performance by determining the area under the receiver operating characteristic curve (AUC). The receiver operating characteristic (ROC) indicates model accuracy when categorizing living and dead trees. A value of 0.5 is the same as chance, and 1.0 implies a perfect fit. Values between 0.7 and 0.8 represent acceptable discrimination, 0.8 to 0.9 indicate good discrimination, and values ≥0.9 reflect excellent discrimination [31]. Multicollinearity of the predictor variables was verified using variance inflation factors (VIF) [32]. A VIF > 4 indicates a multicollinearity [33]. Assessing model accuracy was simplified by providing the threshold value 0.5 to compare models, however, this value can vary depending on management goals [34]. To determine site heterogeneity, we used Duncan’s multiple range test (DMRT) by comparing the means of 9 variables in the three study sites. All the analyses were run using R software v. 3.5.0 [35].
Our three-year post-fire mortality model should be considered to be a preliminary estimate of mortality as additional mortality in the next years is expected. We produced various models based on our logistic regression analyses to provide forest managers with options for estimating mortality using both simple and complicated combinations of variables. To provide a practical example, we presented a two-year post-fire mortality model based on BSI and DBH.

3. Results

3.1. Tree Mortality

Among the 413 Korean red pine in our dataset, 65 (16%) died after three years. In the first year after the fire, tree mortality varied between 0% and 14% among the three sites (Figure 3). Tree mortality increased after three months post fire and continued to increase in a fairly constant manner over the three years at sites A and B. Site C did not have any mortality through to the final observation period. The increasing rate in delayed tree mortality for sites A and B was the greatest between October 2018 and September 2019. Then, the ratio of tree mortality decreased slightly from September 2019 to July 2020.
Topography, tree characteristics, and bark injury showed significant mortality patterns between dead and living trees (student’s t-test p < 0.01, Table 1). The three study sites differed significantly, in pairwise analyses, for six of the nine measured variables (Duncan’s multiple range test (DMRT) p < 0.05, Table 1). The mean slope of the terrain around dead trees (26°) was greater than that for living trees (23°), and dead trees were found at a higher mean elevation (770 m) than living trees (734 m). Dead trees also tended to be located on sites having a northern aspect. Taller trees having a smaller diameter were most susceptible to mortality from the surface wildfire. The dead trees also had a higher BSP and BSH, thus a higher BSI. Site A trees experienced the highest mortality, whereas Site C had zero mortality. Site C trees presented the least bark injury with average BSH, BSP, and BSI scores of 0.7 m, 48.7%, and 1.7, respectively. DBH was greatest in Site C, although Site C was located at the highest elevations relative to the elevations of sites A and B (p < 0.05). the slope varied between 10.8° and 32.7° across all sites, and Sites A and C shared lower slope angles than those measured at Site B (p < 0.05). The aspect varied from northeast to northwest for all trees in the study sites A and B, while 94% of trees located in the site C faced from south to west. We did not obtain a significant role for competition or the soil factors in predicting mortality.

3.2. Factors Influencing the Probability of Mortality Following a Surface Fire

Univariate logistic regression produced models of post-fire mortality of Korean red pine on the basis of topography, tree characteristics, and bark injury (Table 2). The model results showed that DBH, tree height, and bark injuries (BSH, BSP, and BSI) were significant predictors of post-fire tree mortality. Bark scorch variables were crucial for predicting mortality after a surface fire, being a factor involved in all of the strongest logistic regression models. The best univariate models for Korean post-fire mortality contained bark injury variables. The ROC curve analysis showed that the three bark injury models discriminated between living and dead trees at a reliable level (AUC = 0.81–0.86, p < 0.001). BSH and BSP produced the models having the lowest AIC and highest AUC values among single-factor models. The best univariate model was produced when both of these single factors were combined, i.e., BSI (AIC = 255). DBH and tree height each showed intermediate levels of model fitness. Models that included a measure of crown height and topography (elevation and slope) did not discriminate well, owing to the near-zero Z-values and weaker p-values. However, the results of topographical factors were likely to indicate that locations on steeper slopes, at lower elevations, and with a northern aspect, increased tree mortality. The VIF values (VIF < 4) of the eight explanatory variables illustrated that multicollinearity issues were not an issue in the univariate models.
Post-fire tree mortality increased as slope increased (Figure 4). We observed that dead trees tended to be found on north-facing slopes, with all the dead trees occurring at northern aspects (N 66%, NE 29%, and NW 5%). DBH was a significant predictor having an inverse relationship with the probability of mortality. However, trees with a DBH of 45–50 cm showed an opposite trend, likely because of the low number (n = 15) of samples. Mortality also decreased with tree height for trees taller than 20 m. For trees less than 20 m in height, the percent mortality increased as tree height increased. Crown height followed the same pattern. Tree mortality was likely to increase as BSI increased. BSI values between 6 and 8 showed the largest increase in tree mortality, increasing from 17% to 47% over three years. The trees with BSI values ≤4 were more likely to survive, having only 3% tree mortality. However, trees with a BSI value between 8 and 10 did not follow the rising trend.
The multiple logistic regression analyses used seven predictors (elevation, slope, aspect, DBH, tree height, crown height, and BSI) together to produce the best model, as well as their coefficients (Table 3 and Table 4). The results showed that the four formulas (all variables, tree characteristic variables with BSI, slope + DBH + BSI and DBH + BSI) successfully predicted post-fire tree mortality. When considering each formula of variables, the formulas that included DBH and BSI were all significant. Among the significant formulas, three variables (slope, DBH, and BSI) produced the best model for post-fire tree mortality. In particular, samples on steeper slopes (odds ratios (OR) = 1.07), with smaller tree diameters (OR = 0.91), and higher BSI values (OR = 1.45) were more likely to experience higher tree mortality (Figure 5). Although tree height was predictive of tree mortality in the univariate analysis, it was not an independent predictor in the logistic regression analysis. Conversely, slope was a significant predictor with DBH and BSI in the optimal logistic regression analysis model, although it was not predictive in the univariate analysis alone. The best multiple logistic regression model combined slope, DBH, and BSI to produce the model with the lowest AIC value and the highest AUC (Table 3).
The multiple logistic regression analyses revealed that the probability of Korean red pine mortality following surface wildfire was likely to increase with BSI, especially for trees having a low DBH (Figure 6a). According to the two-year post-fire tree mortality model, DBH and BSI were major drivers explaining delayed tree mortality after a surface wildfire. Standing trees having a 5 cm BSI and 50 cm DBH made up 4% of the mortality after fire. However, trees having a DBH = 10 cm represented 52% mortality after fire. A tree characterized by a DBH of 30 cm and a BSI of 10 experienced 56% mortality, whereas an identical DBH-valued tree with a BSI of 15 had a higher mortality (89%).
The probability table for Korean red pine mortality following a surface fire was produced, based on the two-year post-fire tree mortality. This scenario represents a practical example for forest managers on the basis of combining BSI and DBH (Figure S3). We excluded slope in the two-year post-fire tree mortality model because of the low deviance (8.2). The harvesting guideline that we established was a probability of mortality of ≥30%; a single tree had a probability of mortality of more than 30% if BSI ≥ 6 and DBH of 20 cm. However, trees having a DBH ≥ 48 cm and BSI < 13 had a probability of mortality of 1–27% mortality. This post-fire harvesting guideline, i.e., >30% is only an example and is not absolute (Figure S3); thus, it can be adjusted according to the user’s needs or the particular region.
The predictive multiple logistic regression model for the three-year post-fire tree mortality was improved by adding slope after DBH and BSI (Figure 6b). The line indicates that the trends in post-fire tree mortality, when using DBH and BSI, moved markedly to the left, as the slope increased. Post-fire tree mortality was greatly influenced by slope when holding DBH and BSI constant. For instance, a tree with a DBH of 50 cm, BSI of 10, and slope of 0° had 6% probability of tree death, whereas the same tree on 20–30° slope experienced a 23% probability of death. A BSI ≥ 25 almost ensured tree mortality regardless of DBH and slope.

4. Discussion

Forest managers require reliable predictions of post-fire tree mortality under climate change scenario (warmer and drier) to adapt harvesting plans and meet commercial and conservation goals [24]. This study provided a dataset of Korean red pine trees monitored over three years following a surface fire. We focused on the pine plantation forests, which cover a large portion of Korea’s temperate zone, and we combined topography with tree characteristics and fire-related bark injuries for our models. To the best of our knowledge, our study is the first to present a tool for fast decision making for land managers, applicable upon the first year following a surface fire in a Korean red pine forest. Our model presents a more powerful index than previous approaches and predicts post-fire tree mortality in a simpler and faster method for practitioners.

4.1. Post-Fire Korean Red Pine Mortality Modeling

Measures of bark injuries were the most important predictors of Korean red pine mortality following a surface fire. We measured the bark injuries as a surrogate for cambium status; it remains uncertain, however, how well the bark injuries reflect cambium status. Previous studies used more complex codes or classification systems of bark or bole injuries to predict post-fire tree mortality [20,36]. Our novel index (BSI) combined BSH and BSP to better represent cambium status and predict post-fire tree mortality. We showed that BSI was the best univariate model among all measured variables.
Tree size is used extensively as a critical variable in evaluating tree resistance to fire damage [37]. However, DBH was the only significant variable in our models among the measured tree characteristic variables, because tree and crown height did not directly relate to trees affected by the surface fire. Our models predict mortality to decrease with increasing DBH, thereby confirming a pattern observed in previous studies [38,39]. The DBH is closely related to bark thickness and is used in many models instead of bark thickness to indicate resistance to basal heating [19,37,40]. However, some post-fire tree mortality models have reported that DBH was not a significant predictor of post-fire tree mortality [24]. DBH’s failure as a predictor likely occurred because DBH could not represent bark thickness very well in the sampled conifer species in those studies. Therefore, tree diameter in Korean red pine may be a significant predictor of post-fire mortality because of its close correlation with bark thickness.
Topographical factors also influenced Korean red pine mortality following a surface fire. Dead trees are likely to be at lower elevations, on steeper slopes, and with a northern aspect. For instance, Site C, having only living trees, had opposite trends for topographical traits relative to the other sites (A and B). Lydersen and North (2012) [41] found that southwestern aspects (Northern Hemisphere) had significantly greater proportions of bark scorch, whereas differences between northeast- and southwest-facing slopes were generally small. Likewise, we expected a similar pattern that most dead trees had northwest to northeast aspects, due to the wind blowing from the south-western direction. However, this question was not confirmed in this study because the aspects between dead and living trees were not significantly different, probably due to high variation among sites and the low number of replications in this study (only three sites). Moreover, stands burned by low and moderate fires are found on less steep plots (slope <18%), probably owing to lower stand densities on these sites [42]. Topography can directly influence fire intensity [43], fire frequency [44], and forest vegetation by altering environmental factors, including soil moisture [45], soil thickness [46], and microclimate [47]. Previous studies [41] have also demonstrated the effect of topography on tree density [48] and species composition [44].
Our results from the logit models confirmed that both tree diameter and bark injury are important predictors of conifer mortality following a fire. These findings confirm previous studies [20,25,36]. Our models showed that BSI and DBH were both robust individual predictors of post-fire tree mortality and performed nearly as well as the best model. The best set of predictive variables for our logit models included slope, BSI, and DBH; however, excluding slope from the model would depend on the practitioner’s preferences when it is a gentle slope. Although we applied our accuracy assessment with a 0.5 cutoff to simply model comparisons, practitioners can customize the post-fire tree mortality guidelines by selecting different or more applicable cutoff levels and by examining ecological elements.

4.2. Potential Factors Influencing Post-Fire Tree Mortality

The delayed mortality of bark-scorched trees could also be potentially influenced by root injury, inter- and intraspecific competition, insects, and pathogens. As fires can markedly decrease fine-root biomass, which would reduce water uptake [49], indirect tree mortality can also occur through hydraulic limitation caused by root injuries [23]. Competition with neighboring trees can also influence post-fire tree mortality by limiting essential soil resources, because slower-growing trees are more likely to die than faster-growing trees for trees having similar levels of fire-induced injury [50]. Insects and pathogens can also alter delayed tree mortality by reducing growth and impairing conductivity [23,51,52]. In this study, we also analyzed the inter-tree competition and soil properties. However, we did not obtain a significant result in predicting mortality, likely because of the relatively low fire intensity.
Prediction of post-fire tree mortality can also be influenced by tree species. The previous models predicting Douglas-fir (Pseudotsuga menziesii (Mirb.) Franco) and ponderosa pine (Pinus ponderosa) mortality following wildfire were optimal when they included variables of crown volume, cambium kill rating, and the presence of beetles [53]. The best model for predicting sugar pine (Pinus lambertiana) mortality included percent crown length, cambium kill rating, and beetle attacks, whereas the best white fir model also included DBH [24]. Predicting post-fire tree mortality has also relied on fire intensity, bole scorch height, bark thickness, and tree height [54,55]. Thus, it is important to identify those variables that are optimal for a given tree species.
Estimating a cambium kill rating is laborious, and therefore multiple post-fire tree mortality models lack this measurement [56]. Models without a cambium kill rating are indeed less accurate and have lower discriminatory power, thereby demonstrating the trade-off between cost/labor and reliability [53]. Nonetheless, field workers may prefer using simpler measurements for cost-efficiency purposes to have broader applicability of these models. For instance, the Ryan and Reinhardt model has been extensively applied in fire management settings despite lacking a cambium kill rating component [57]. Although our model did not produce the highest accuracy relative to similar models, our use of DBH and BSI to represent cambium status provided an acceptable average discriminatory ability of 87%. Further studies with our model should include a cambium kill rating and many replicates to improve our post-fire tree mortality predictions.

4.3. Issues Related to Our Experimental Design

Most post-fire tree mortality occurs within the first two years after a fire, and, in North American forests, patterns of fire-related tree death stabilize by the third year [20,24,55]. We also found a similar tree mortality pattern in Korean red pine dominated forests, i.e., an intensive pulse of mortality in the first two years after the surface fire and a reduced mortality rate over the third year. We expect tree mortality to continue rising, although at a lower rate over the following years. However, a longer-term study of tree mortality is necessary, as we did not observe a flattened mortality curve (rate at zero) by the end of our study, and patterns of tree survival can be predicted at least ten years after a fire [58].
The major determinants of fire behavior and burn severity were not estimated in the burned stands because we confined our study to surface fires, and Korean land managers consider crown-damaged trees as trees to be harvested. Thus, we focused on trees with intact crowns after a surface wildfire to provide practical knowledge to Korean forest managers. Moreover, constraining models to certain types of fire can heighten the accuracy of post-fire tree mortality predictions.
The regional characteristics may also indirectly influence site effects on post-fire tree mortality. Our data is specific to the Samcheok region, located near the East Sea, and is an area where the summer climate is hot and humid [59] with very cold winters, relative to comparable temperate regions, because of the strong influence of the Siberian High [60]. Measurements of microclimate variables, such as wind, temperature, relative humidity, and local fire fuel, i.e., accumulated deadwood, can enhance model predictability and can clarify the complicated interactions among topography, microclimate, fuels, and fire impacts [42]. Furthermore, similar studies of other tree species and sites are needed to determine the general importance of slope, BSI, and DBH as predictors of mortality after a surface fire.

5. Conclusions

The increasing number of large wildfires in Korean red pine forests has focused attention on predicting post-fire tree mortality. The development of tree mortality models using physiological factors such as cambial status is recommended [61]; however, such measurements are costly and have not yet produced reliable models. This study provides simple and effective models using easy-to-measure bark injuries as a proxy of direct cambium sampling. This new tool can help forest managers to establish criteria for selecting those trees available for salvage post-fire. Our univariate models for Korean red pine relied on the use of the novel BSI index to produce the best univariate models of post-fire tree mortality. A logit model combining slope, DBH, and BSI offered the best predictions. The simple nature of the three selected variables demonstrated that practitioners could easily and reliably estimate the three-year tree mortality after a surface fire. Further studies should test whether slope, DBH, and BSI are the three best predictors of Korean red pine mortality in other regions, regardless of differences in fire intensity and fire timing. Our models will help forest managers to reduce the number of trees salvaged after a fire, and therefore help forest ecosystems remain resilient following more frequent and severe fire disturbances expected under a warmer and drier climate change scenario.

Supplementary Materials

The following are available online at https://www.mdpi.com/2071-1050/13/2/569/s1, Figure S1: The wind rose chart during the Samcheok wildfire in 2017. Figure S2: Meteorological data observed during the Samcheok wildfire in 2017. The figure represents wind speed and relative humidity. Figure S3: A Korean practical example of harvesting criterion: 30% based on the multiple logistic regression model with a binary response variable for predicting the two-year post-fire tree mortality in relation to the diameter at breast height (DBH) and the bark scorch index (BSI). The green area represents trees to retain while brown is for harvesting.

Author Contributions

Conceptualization and methodology, W.K.; validation, S.K. (Semyung Kwon), S.K. (Sanghyun Kim), and M.M.G.; formal analysis, S.K. (Semyung Kwon); investigation, S.K. (Semyung Kwon), W.K., M.M.G., S.K. (Sanghyun Kim), J.K., C.-B.K., and K.-H.P.; data curation, S.K. (Semyung Kwon) and S.K. (Sanghyun Kim); writing—original draft preparation, S.K. (Sanghyun Kim); writing—review and editing, S.K. (Sanghyun Kim) and M.M.G.; visualization, S.K. (Semyung Kwon); supervision, M.M.G. and W.K.; project administration, W.K.; resources: W.K. and M.M.G.; funding acquisition, W.K. and M.M.G. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the “Long-term Monitoring of Ecosystem Changes on Natural Regenerated Forest and Plantation in Forest Fire Damaged Area (FE-0500-2016-03)” of the National Institute of Forest Science, Republic of Korea and the silvicultural research grant of MRC-Abitibi (https://mrcabitibi.qc.ca/) obtained by Miguel Montoro Girona.

Data Availability Statement

Data available on request due to patent issues. The data presented in this study are available on request from the corresponding author. The data are not publicly available due to the future patent.

Acknowledgments

We gratefully appreciate all the field assistants. We extend special thanks to M. Hay for the English editing and the Group in Ecology at the MRC Abitibi (GREMA) for the support.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Montoro Girona, M.; Navarro, L.; Morin, H. A Secret Hidden in the Sediments: Lepidoptera Scales. Front. Ecol. Evol. 2018, 6, 1–5. [Google Scholar] [CrossRef] [Green Version]
  2. Montoro Girona, M.; Morin, H.; Lussier, J.-M.; Ruel, J.-C. Post-cutting Mortality Following Experimental Silvicultural Treatments in Unmanaged Boreal Forest Stands. Front. For. Glob. Chang. 2019, 2, 4. [Google Scholar] [CrossRef] [Green Version]
  3. Bergeron, Y.; Gauthier, S.; Kafka, V.; Lefort, P.; Lesieur, D. Natural fire frequency for the eastern Canadian boreal forest: Consequences for sustainable forestry. Can. J. For. Res. 2001, 31, 384–391. [Google Scholar] [CrossRef]
  4. Martin, M.; Girona, M.M.; Morin, H. Driving factors of conifer regeneration dynamics in eastern Canadian boreal old-growth forests. PLoS ONE 2020, 15, e0230221. [Google Scholar] [CrossRef] [PubMed]
  5. Lavoie, J.; Girona, M.M.; Morin, H. Vulnerability of conifer regeneration to spruce budworm outbreaks in the Eastern Canadian boreal forest. Forests 2019, 10, 850. [Google Scholar] [CrossRef] [Green Version]
  6. Labrecque-Foy, J.P.; Morin, H.; Girona, M.M. Dynamics of territorial occupation by North American beavers in Canadian boreal forests: A novel dendroecological approach. Forests 2020, 11, 221. [Google Scholar] [CrossRef] [Green Version]
  7. Navarro, L.; Morin, H.; Bergeron, Y.; Girona, M.M. Changes in Spatiotemporal Patterns of 20th Century Spruce Budworm Outbreaks in Eastern Canadian Boreal Forests. Front. Plant Sci. 2018, 9, 1–15. [Google Scholar] [CrossRef] [Green Version]
  8. Seidl, R.; Thom, D.; Kautz, M.; Martin-Benito, D.; Peltoniemi, M.; Vacchiano, G.; Wild, J.; Ascoli, D.; Petr, M.; Honkaniemi, J.; et al. Forest disturbances under climate change. Nat. Clim. Chang. 2017, 7, 395–402. [Google Scholar] [CrossRef] [Green Version]
  9. Lucash, M.S.; Scheller, R.M.; Sturtevant, B.R.; Gustafson, E.J.; Kretchun, A.M.; Foster, J.R. More than the sum of its parts: How disturbance interactions shape forest dynamics under climate change. Ecosphere 2018, 9, e02293. [Google Scholar] [CrossRef]
  10. Westerling, A.L.; Hidalgo, H.G.; Cayan, D.R.; Swetnam, T.W. Warming and Earlier Spring Increase Western, U.S. Forest Wildfire Activity. Science 2006, 313, 940–943. [Google Scholar] [CrossRef] [Green Version]
  11. Bentz, B.J.; Rgnire, J.; Fettig, C.J.; Hansen, E.M.; Hayes, J.L.; Hicke, J.A.; Kelsey, R.G.; Negron, J.F.; Seybold, S.J. Climate change and bark beetles of the western United States and Canada: Direct and indirect effects. Bioscience 2010, 60, 602–613. [Google Scholar] [CrossRef]
  12. Wan, H.Y.; Rhodes, A.C.; St. Clair, S.B. Fire severity alters plant regeneration patterns and defense against herbivores in mixed aspen forests. Oikos 2014, 123, 1479–1488. [Google Scholar] [CrossRef]
  13. Fettig, C.J.; Reid, M.L.; Bentz, B.J.; Sevanto, S.; Spittlehouse, D.L.; Wang, T. Changing climates, changing forests: A western North American perspective. J. For. 2013, 111, 214–228. [Google Scholar] [CrossRef]
  14. Marques, S.; Marto, M.; Bushenkov, V.; McDill, M.; Borges, J.G. Addressing wildfire risk in forest management planning with multiple criteria decision making methods. Sustainability 2017, 9, 298. [Google Scholar] [CrossRef] [Green Version]
  15. Korea Forest Service. Statistical Yearbook of Forestry; Korea Forest Service: Daejeon, Korea, 2018. [Google Scholar]
  16. Lee, S.-Y.; An, S.-H. Comparative Analysis of Mortality in Species of Trees after Surface Forest Fire. J. Korean Soc. Hazard Mitig. 2009, 9, 39–43. [Google Scholar]
  17. Korea Forest Service. Forest Fire Statistical Yearbook 2019; Korea Forest Service: Daejeon, Korea, 2019. [Google Scholar]
  18. Korea Forest Service. Forestry Management Survey; Korea Forest Service: Daejeon, Korea, 2019. [Google Scholar]
  19. Stevens-Rumann, C.S.; Sieg, C.H.; Hunter, M.E. Ten years after wildfires: How does varying tree mortality impact fire hazard and forest resiliency? For. Ecol. Manag. 2012, 267, 199–208. [Google Scholar] [CrossRef]
  20. Hood, S.M.; Smith, S.L.; Cluck, D.R. Delayed Conifer Tree Mortality Following Fire in California; USDA Forest Service: Albany, CA, USA, 2007; Volume PSW-GTR-20. [Google Scholar]
  21. Haugaasen, T.; Barlow, J.; Peres, C.A. Surface wildfires in central Amazonia: Short-term impact on forest structure and carbon loss. For. Ecol. Manag. 2003, 179, 321–331. [Google Scholar] [CrossRef]
  22. Michaletz, S.T.; Johnson, E.A. A biophysical process model of tree mortality in surface fires. Can. J. For. Res. 2008, 38, 2013–2029. [Google Scholar] [CrossRef]
  23. Bär, A.; Michaletz, S.T.; Mayr, S. Fire effects on tree physiology. New Phytol. 2019, 223, 1728–1741. [Google Scholar] [CrossRef] [Green Version]
  24. Hood, S.M.; Smith, S.L.; Cluck, D.R. Predicting mortality for five California conifers following wildfire. For. Ecol. Manag. 2010, 260, 750–762. [Google Scholar] [CrossRef]
  25. Vega, J.; Jimenez, E.; Vega, D.; Ortiz, L.; Pérez, J.R. Pinus pinaster Ait. tree mortality following wildfire in Spain. For. Ecol. Manag. 2011, 261, 2232–2242. [Google Scholar] [CrossRef]
  26. Kim, S.; Axelsson, E.P.; Girona, M.M.; Senior, J.K. Continuous-cover forestry maintains soil fungal communities in Norway spruce dominated boreal forests. For. Ecol. Manag. 2021, 480, 118659. [Google Scholar] [CrossRef]
  27. Han, S.H.; Kim, J.H.; Kang, W.S.; Hwang, J.H.; Park, K.H.; Kim, C.B. Monitoring Soil Characteristics and Growth of Pinus densiflora Five Years after Restoration in the Baekdudaegan Ridge. Korean J. Environ. Ecol. 2019, 33, 453–461. [Google Scholar] [CrossRef]
  28. Biging, G.S.; Dobbertin, M. Evaluation of Competition Indices in Individual Tree Growth Models. For. Sci. 1995, 41, 360–377. [Google Scholar]
  29. Girona, M.M.; Morin, H.; Lussier, J.M.; Walsh, D. Radial growth response of black spruce stands ten years after experimental shelterwoods and seed-tree cuttings in boreal forest. Forests 2016, 7, 240. [Google Scholar] [CrossRef] [Green Version]
  30. Burnham, K.; Anderson, D. Model Selection and Multimodel Inference; Springer: New York, NY, USA, 2004. [Google Scholar]
  31. Hosmer, D.W.; Lemeshow, S. Applied Logistic Regression; Wiley: New York, NY, USA, 2000; ISBN 9780471722144. [Google Scholar]
  32. Graham, M.H. Confronting multicollinearity in ecological multiple regression. Ecology 2003, 84, 2809–2815. [Google Scholar] [CrossRef] [Green Version]
  33. Kalantar, B.; Pradhan, B.; Naghibi, S.A.; Motevalli, A.; Mansor, S. Assessment of the effects of training data selection on the landslide susceptibility mapping: A comparison between support vector machine (SVM), logistic regression (LR) and artificial neural networks (ANN). Geomat. Nat. Hazards Risk 2018, 9, 49–69. [Google Scholar] [CrossRef]
  34. Freeman, E.; Gretchen, M. PresenceAbsence: An R Package for Presence Absence Analysis. J. Stat. Softw. 2007, 23, 31. [Google Scholar]
  35. Rs. Team. Rstudio: Integrated Development for R; PBC: Boston, MA, USA, 2020; Available online: http://www.rstudio.com/ (accessed on 28 July 2020).
  36. Hood, S.M.; Cluck, D.R.; Smith, S.L.; Ryan, K.C. Using Bark Char Codes to Predict Post-fire Cambium Mortality. Fire Ecol. 2008, 4, 57–73. [Google Scholar] [CrossRef]
  37. Ryan, K.C.; Reinhardt, E.D. Predicting postfire mortality of seven western conifers. Can. J. For. Res. 1988, 18, 1291–1297. [Google Scholar] [CrossRef]
  38. Stephens, S.L.; Finney, M.A. Prescribed fire mortality of Sierra Nevada mixed conifer tree species: Effects of crown damage and forest floor combustion. For. Ecol. Manag. 2002, 162, 261–271. [Google Scholar] [CrossRef]
  39. Mutch, L.S.; Parsons, D.J. Mixed Conifer Forest Mortality and Establishment Before and After Prescribed Fire in Sequoia National Park, California. For. Sci. 1998, 44, 341–355. [Google Scholar]
  40. Kelley, D.I.; Harrison, S.P.; Prentice, I.C. Improved simulation of fire–vegetation interactions in the Land surface Processes and eXchanges dynamic global vegetation model (LPX-Mv1). Geosci. Model Dev. 2014, 7, 2411–2433. [Google Scholar] [CrossRef] [Green Version]
  41. Lydersen, J.; North, M. Topographic Variation in Structure of Mixed-Conifer Forests Under an Active-Fire Regime. Ecosystems 2012, 15, 1134–1146. [Google Scholar] [CrossRef]
  42. Lentile, L.B.; Smith, F.W.; Shepperd, W.D. Influence of topography and forest structure on patterns of mixed severity fire in ponderosa pine forests of the South Dakota Black Hills, USA. Int. J. Wildland Fire 2006, 15, 557–566. [Google Scholar] [CrossRef]
  43. Beaty, R.M.; Taylor, A.H. Spatial and Temporal Variation of Fire Regimes in a Mixed Conifer Forest Landscape, Southern Cascades, California, USA. J. Biogeogr. 2001, 28, 955–966. [Google Scholar] [CrossRef]
  44. Taylor, A.H.; Skinner, C.N. Spatial patterns and controls on historical fire regimes and forest structure in the Klamath Mountains. Ecol. Appl. 2003, 13, 704–719. [Google Scholar] [CrossRef]
  45. Scholl, A.E.; Taylor, A.H. Fire regimes, forest change, and self-organization in an old-growth mixed-conifer forest, Yosemite National Park, USA. Ecol. Appl. 2010, 20, 362–380. [Google Scholar] [CrossRef]
  46. Meyer, M.D.; North, M.P.; Gray, A.N.; Zald, H.S.J. Influence of soil thickness on stand characteristics in a Sierra Nevada mixed-conifer forest. Plant Soil 2007, 294, 113–123. [Google Scholar] [CrossRef]
  47. Abella, S.R.; Denton, C.W. Spatial variation in reference conditions: Historical tree density and pattern on a pinus ponderosa landscape. Can. J. For. Res. 2009, 39, 2391–2403. [Google Scholar] [CrossRef]
  48. Underwood, E.C.; Viers, J.H.; Quinn, J.F.; North, M. Using topography to meet wildlife and fuels treatment objectives in fire-suppressed landscapes. Environ. Manag. 2010, 46, 809–819. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  49. Smirnova, E.; Bergeron, Y.; Brais, S.; Granström, A. Postfire root distribution of Scots pine in relation to fire behaviour. Can. J. For. Res. 2008, 38, 353–362. [Google Scholar] [CrossRef]
  50. Nesmith, J.C.B.; Das, A.J.; O’Hara, K.L.; van Mantgem, P.J. The influence of prefire tree growth and crown condition on postfire mortality of sugar pine following prescribed fire in Sequoia National Park. Can. J. For. Res. 2015, 45, 910–919. [Google Scholar] [CrossRef]
  51. Hood, S.M.; Varner, J.M.; Van Mantgem, P.; Cansler, C.A. Fire and tree death: Understanding and improving modeling of fire-induced tree mortality. Environ. Res. Lett. 2018, 13, 113004. [Google Scholar] [CrossRef]
  52. Kane, J.M.; Varner, J.M.; Metz, M.R.; van Mantgem, P.J. Characterizing interactions between fire and other disturbances and their impacts on tree mortality in western U.S. Forests. For. Ecol. Manag. 2017, 405, 188–199. [Google Scholar] [CrossRef]
  53. Ganio, L.M.; Progar, R.A. Mortality predictions of fire-injured large Douglas-fir and ponderosa pine in Oregon and Washington, USA. For. Ecol. Manag. 2017, 390, 47–67. [Google Scholar] [CrossRef] [Green Version]
  54. Rigolot, E. Predicting postfire mortality of Pinus halepensis Mill. and Pinus pinea L. Plant Ecol. 2004, 171, 139–151. [Google Scholar] [CrossRef]
  55. Sieg, C.H.; McMillin, J.D.; Fowler, J.F.; Allen, K.K.; Negron, J.F.; Wadleigh, L.L.; Anhold, J.A.; Gibson, K.E. Best predictors for postfire mortality of ponderosa pine trees in the Intermountain West. For. Sci. 2006, 52, 718–728. [Google Scholar]
  56. Thies, W.G.; Westlina, D.J.; Loewen, M.; Brenner, G. Prediction of delayed mortality of fire-damaged ponderosa pine following prescribed fires in eastern Oregon, USA. Int. J. Wildland Fire 2006, 15, 19–29. [Google Scholar] [CrossRef]
  57. Lutes, D.C. FOFEM User Guide; USDA Forest Service: Albany, CA, USA, 2020. [Google Scholar]
  58. Angers, V.A.; Gauthier, S.; Drapeau, P.; Jayen, K.; Bergeron, Y. Tree mortality and snag dynamics in North American boreal tree species after a wildfire: A long-term study. Int. J. Wildland Fire 2011, 20, 751–763. [Google Scholar] [CrossRef]
  59. Gangwon Province. Nature Environment (in Korean). 2018. Available online: http://www.provin.gangwon.kr/gw/portal/sub03_01 (accessed on 16 July 2020).
  60. Korea Meteorological Administration. Korean Climate: Gangwon-do. 2020. Available online: https://www.weather.go.kr/weather/climate/average_regional01.jsp (accessed on 16 July 2020).
  61. Hood, S.M.; Morgan, J. Post-fire Tree Mortality. In Encyclopedia of Wildfires Wildland-Urban Interface Fires; Springer: Berlin/Heidelberg, Germany, 2019; p. 10. [Google Scholar]
Sustainability 13 00569 g001 550
Figure 1. (a) Location (37°15′53″ N, 129°1′24″ E, 700–800 m asl) and topographical map of the study area near Samcheok, the Republic of Korea. Each point represents a tagged tree that was scorched by a surface wildfire in May 2017; (b) A pair of satellite photos showing the study area before and after the May 2017 Samcheok wildfire.
Figure 1. (a) Location (37°15′53″ N, 129°1′24″ E, 700–800 m asl) and topographical map of the study area near Samcheok, the Republic of Korea. Each point represents a tagged tree that was scorched by a surface wildfire in May 2017; (b) A pair of satellite photos showing the study area before and after the May 2017 Samcheok wildfire.
Sustainability 13 00569 g001
Sustainability 13 00569 g002 550
Figure 2. Calculations of the bark scorch index (BSI). (a) Photos showing measurements of two fire-scarred trunks. The left photo study tree has a bark scorch height (BSH) of 1.9 and a bark scorch proportion (BSP) of 0.8. These values for the trunk on the right photo are 1.3 and 0.6, respectively; (b) Bark scorch index (BSI) for a trunk is calculated by multiplying BSP and BSH for each quadrant, and then the four obtained values are summed. For this example, the BSI is (2.0 × 0.8) + (1.5 × 0.45) + (1.2 × 0.65) + (0.7 × 0.5) = 3.405.
Figure 2. Calculations of the bark scorch index (BSI). (a) Photos showing measurements of two fire-scarred trunks. The left photo study tree has a bark scorch height (BSH) of 1.9 and a bark scorch proportion (BSP) of 0.8. These values for the trunk on the right photo are 1.3 and 0.6, respectively; (b) Bark scorch index (BSI) for a trunk is calculated by multiplying BSP and BSH for each quadrant, and then the four obtained values are summed. For this example, the BSI is (2.0 × 0.8) + (1.5 × 0.45) + (1.2 × 0.65) + (0.7 × 0.5) = 3.405.
Sustainability 13 00569 g002
Sustainability 13 00569 g003 550
Figure 3. Three years of cumulative tree mortality after the May 2017 surface fire. Measurements took place on 11 July 2018, 19 October 2018, 20 September 2019, and 8 July 2020.
Figure 3. Three years of cumulative tree mortality after the May 2017 surface fire. Measurements took place on 11 July 2018, 19 October 2018, 20 September 2019, and 8 July 2020.
Sustainability 13 00569 g003
Sustainability 13 00569 g004 550
Figure 4. Bar charts showing the percentage of cumulative tree mortality by (a) slope; (b) diameter at breast height (DBH); (c) tree height; and (d) bark scorch index (BSI) value.
Figure 4. Bar charts showing the percentage of cumulative tree mortality by (a) slope; (b) diameter at breast height (DBH); (c) tree height; and (d) bark scorch index (BSI) value.
Sustainability 13 00569 g004
Sustainability 13 00569 g005 550
Figure 5. Odds ratio plot of the final model of Korean red pine mortality following a surface wildfire where DBH and BSI correspond to the diameter at breast height and bark scorch index, respectively. DBH, diameter at breast height. * Represents significant values according to *** p < 0.001 and ** p < 0.01.
Figure 5. Odds ratio plot of the final model of Korean red pine mortality following a surface wildfire where DBH and BSI correspond to the diameter at breast height and bark scorch index, respectively. DBH, diameter at breast height. * Represents significant values according to *** p < 0.001 and ** p < 0.01.
Sustainability 13 00569 g005
Sustainability 13 00569 g006 550
Figure 6. Multiple logistic regression models with a binary response variable, (a) For predicting the two-year post-fire tree mortality in relation to diameter at breast height (DBH) and the bark scorch index; (b) For predicting the three-year post-fire tree mortality in relation to the bark scorch index (BSI), DBH, and various slopes between 0° and 30°.
Figure 6. Multiple logistic regression models with a binary response variable, (a) For predicting the two-year post-fire tree mortality in relation to diameter at breast height (DBH) and the bark scorch index; (b) For predicting the three-year post-fire tree mortality in relation to the bark scorch index (BSI), DBH, and various slopes between 0° and 30°.
Sustainability 13 00569 g006
Table
Table 1. Summary of variables for tree mortality following the May 2017 surface wildfire. Values are the mean and standard deviation, measured in July 2018.
Table 1. Summary of variables for tree mortality following the May 2017 surface wildfire. Values are the mean and standard deviation, measured in July 2018.
CategoriesVariablesDead TreesLive TreesSite ASite BSite C
TopographyElevation (m)734.4 (±15.7) **769.6 (±86.0) **726.3 (±6.1) b737.6 (±19) b849.5 (±5.6) a
Slope (°)26.3 (±7.8) *23.1 (±8.0) *21.2 (±8.9) b26.4 (±6.3) a19.0 (±8.2) b
AspectNNENNESW
Growth featureDBH (cm)30.6 (±6.7) **34.3 (±7.5) **33.5 (±7.7) b32.1 (±6.8) b36.9 (±7.6) a
TH (m)18.9 (±2.4) **17.3 (±3.2) **22.0 (±1.6) a17.9 (±2.8) b15.3 (±1.9) c
CH (m)11.5 (±2.4) **9.5 (±5.3) **13.7 (±2.3) a10.72 (±5.6) b6.7 (±1.5) c
Bark injuryBSH (m)2.6 (±1.1) **1.3 (±1.0) **2.7 (±1.7) a1.7 (±0.9) b0.7 (±0.4) c
BSP (%)78.9 (±14.8) **57.9 (±20.4) **74.19 (±18.3) a65.6 (±20.1) b48.7 (±18.1) c
BSI (ns)8.6 (±3.9) **3.6 (±3.2) **8.8 (±6.0) a5.0 (±3.0) b1.7 (±1.3) c
Student’s t-test, ** p < 0.01, * p < 0.05; Parameter values followed by the same letter are not significantly different (Duncan, p < 0.05). Duncan’s multiple range test (DMRT) of the variables by sites. DBH, diameter at breast height; TH, tree height; CH, crown height; BSH, bark scorch height; BSP, bark scorch proportion; BSI, bark scorch index (BSH × BSP). Aspect is categorical for eight cardinal directions and was not tested for significant differences among dead and live trees, nor sites.
Table
Table 2. The univariate logistic regression models in the selection of topographical, tree characteristic, and bark injury variables for Korean red pine mortality following surface wildfire.
Table 2. The univariate logistic regression models in the selection of topographical, tree characteristic, and bark injury variables for Korean red pine mortality following surface wildfire.
CategoryVariableEstimateStd. ErrorZpAIC∆AICAUCVIF
TopographyElevation−0.0050.002−2.40.018329.373.90.6821.034
Slope0.0510.0192.80.005330.174.70.5971.249
Growth featureDBH−0.070.02−3.5<0.001325.269.80.6411.426
TH0.1690.0473.6<0.001324.368.90.6522.449
CH0.0850.0412.10.039331.275.80.7122.16
Bark injuryBSH1.0860.1557<0.001268.312.90.8381.214
BSP21.2632.9977.1<0.001270.515.10.8081.145
BSI4.0751.2853.2<0.001255.400.862°
DBH, diameter at breast height; TH, tree height; CH, crown height; BSH, bark scorch height; BSP, bark scorch proportion; BSI, bark scorch index (BSH × BSP). The best models were determined according to Akaike model selection approach (∆AIC < 2). VIF, variance inflation factors.
Table
Table 3. Multiple logistic regression analyses of independent predictors for Korean red pine mortality over three years following a surface wildfire.
Table 3. Multiple logistic regression analyses of independent predictors for Korean red pine mortality over three years following a surface wildfire.
ModelX2Log-LikelihoodMcFadden’s Pseudo R2AIC∆AICAUC
M15.32−115.720.31245.455.490.87
M211.52−159.820.04325.6485.70.673
M312.7 **−122.280.27252.5612.60.864
M431.62 ***−150.940.09309.8869.90.733
M57.02−119.780.28249.569.60.867
M67.27−115.980.31239.9600.87
M77.24−119.840.28245.685.70.867
M816.16 **−122.640.27251.2811.30.862
X2, The Hosmer–Lemeshow test is a statistical test of goodness of fit for logistic regression models. Bold lettering refers to the best model based on the ∆AIC values. Model factors are in Table 4. *** p < 0.001 and ** p < 0.01.
Table
Table 4. Coefficients in the multiple logistic regression models.
Table 4. Coefficients in the multiple logistic regression models.
VariablesModelFormula
AllM1−1.6278 + (−0.0014 × Elevation) + (0.0645 × Slope **) + (−0.0909 × DBH **) + (0.0201 × TH)
+ (−0.0010 × CH) + (0.3632 × BSI ***)
TopographyM20.7224 + (−0.0043 × Elevation *) + (0.0645 × Slope *)
M3−3.5282 * + (−0.0017 × Elevation) + (0.0536 × Slope *) + (0.3561 × BSI ***)
Growth featuresM4−2.3565 * + (−0.0882 × DBH ***) + (0.2063 × TH ***) + (0.0001 × CH)
M5−1.0877 + (−0.0861 × DBH **) + (0.0251 × TH) + (−0.0111 × CH) + (0.3546 × BSI ***)
Optimal
models		M6−2.4525 * + (0.0649 × Slope **) + (−0.0892 × DBH ***) + (0.3733 × BSI ***)
M7−0.9018 + (−0.0824 × DBH ***) + (0.3607 × BSI ***)
M8−4.9233 *** + (0.0547 × Slope **) + (0.3621 × BSI ***)
DBH, diameter at breast height; BSI, bark scorch index; CH, crown height; TH, tree height. *** p < 0.001, ** p < 0.01, and * p < 0.05.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Kwon, S.; Kim, S.; Kim, J.; Kang, W.; Park, K.-H.; Kim, C.-B.; Girona, M.M. Predicting Post-Fire Tree Mortality in a Temperate Pine Forest, Korea. Sustainability 2021, 13, 569. https://doi.org/10.3390/su13020569

AMA Style

Kwon S, Kim S, Kim J, Kang W, Park K-H, Kim C-B, Girona MM. Predicting Post-Fire Tree Mortality in a Temperate Pine Forest, Korea. Sustainability. 2021; 13(2):569. https://doi.org/10.3390/su13020569

Chicago/Turabian Style

Kwon, Semyung, Sanghyun Kim, Jeonghwan Kim, Wonseok Kang, Ki-Hyung Park, Chan-Beom Kim, and Miguel M. Girona. 2021. "Predicting Post-Fire Tree Mortality in a Temperate Pine Forest, Korea" Sustainability 13, no. 2: 569. https://doi.org/10.3390/su13020569

APA Style

Kwon, S., Kim, S., Kim, J., Kang, W., Park, K. -H., Kim, C. -B., & Girona, M. M. (2021). Predicting Post-Fire Tree Mortality in a Temperate Pine Forest, Korea. Sustainability, 13(2), 569. https://doi.org/10.3390/su13020569

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop