Empirical modeling of cutblock edge windthrow risk on Vancouver Island, Canada, using stand level information

https://doi.org/10.1016/S0378-1127(00)00620-4Get rights and content

Abstract

Winter storms routinely damage stands along clearcut boundaries on Northern Vancouver Island, disrupting forest management plans. Empirical modelling of windthrow risk is more suitable than mechanistic modelling for the western hemlock (Tsuga heterophylla (Raf.) Sarg.), redcedar (Thuja plicata (Donn)) and amabilis fir (Abies amabilis Forbes.) stands in this area which vary in composition, structure, density and age. Cutblock edge windthrow in a 900 km2 study area was mapped using 1:15000 scale aerial photographs. A geographic information system (GIS) database for 3000–50 m long by 40 m deep cutblock edge segments was assembled using readily available topographic and inventory information. Logistic regression models for predicting probability of damage for various levels of damage severity were evaluated. Key variables in these models included site quality, stocking, boundary orientation, time since harvest and topographic exposure. The models correctly predicted the outcome of 71–76% of test segments, and adequately fit using χ2 tests for test data grouped by stand characteristics. These models can be entered into GIS map calculations to produce landscape level maps of windthrow hazard, or into spreadsheets for evaluating the probability of damage for alternate cutblock design scenarios.

Introduction

Windthrow refers to the uprooting or breakage of trees by wind, and results from the interaction of climatic, topographic, stand, tree and soil factors. The risk of damage is increased by management interventions such as thinning or clearcutting which increase the exposure of residual trees to wind (e.g. Harris, 1989, Ruel, 1995). Results from a province wide census by the British Columbia (BC) Ministry of Forests indicated that a volume of timber equivalent to 4.2% of BC’s annual allowable cut was windthrown in 1 year. This damage disrupts stand and forest level management plans (Mitchell, 1995).

Documenting the relationship between windthrow occurrence and biophysical and management factors enables damage prediction and development of lower risk harvesting plans. The different approaches to windthrow risk assessment can be broadly classified as: observational, mechanistic or empirical. In observational methods, the presence of factors known to be associated with higher incidence of damage is tallied. The relative risk of windthrow is considered to increase with the number of risk indicators observed (e.g. Ruth and Yoder, 1953, Alexander, 1964, Harris, 1989, Stathers et al., 1994). Mechanistic models predict the likelihood of damage based on an evaluation of the critical windspeed for tree failure, and the probability of such a wind occurring at a given location. These models are built by winching trees to determine resistance, placing tree crowns in wind tunnels to determine drag, and using wind tunnel, or numerical modelling to determine local modification of the regional wind regime. Mechanistic models have been successfully developed for single species uniform canopied stands which become increasingly unstable as they increase in height (e.g. Smith et al., 1987, Peltola and Kellomaki, 1993). These models must be calibrated in field studies on a range of sites (Quine, 1995).

For stands with more complex structure and composition, the empirical approach is better suited. In this approach regression models are built which relate presence or magnitude of wind damage in sampling units to their environmental and management attributes. This method has been used for predicting stand edge damage in boreal black spruce strip cuts (Elling and Verry, 1978, Fleming and Crossfield, 1983) and in riparian reserves in coastal Oregon (Steinblums et al., 1984, Andrus and Froehlich, 1992). In these studies, the data were collected using field sampling, and coefficients of determination ranged from 0.57 to 0.93. Fridman and Valinger (1998) developed models for predicting snow and wind damage using tree, stand and site characteristics from permanent sample plots in Scots pine stands in Sweden. Their best model correctly predicted 81% of the undamaged and 82% of the damaged plots. Because windthrow results from the interaction of multiple factors, these equations should not be applied in locations with biophysical or management characteristics different from those used in building the model.

For the construction of empirical models, information about large numbers of sample units is necessary. Geographic information systems (GISs) provide the opportunity to assemble stand or landscape level information, generate spatial variables, analyze, and map the results. Wright and Quine (1993) used a GIS to analyze damage caused by a strong wind event in a forest in North Yorkshire. They divided the mapsheets into 50×50m2 cells, added a layer with three classes of wind damage (none, partial, total) determined through ground mapping by field foresters, and analyzed the association of windthrow with individual topographic, soil and stand variables contained in other GIS layers. While they found associations between damage and specific variables, they did not build a predictive model. Moore and Somerville (1998) in New Zealand used GIS software to run a mechanistic model with a tree/stand critical windspeed component and an airflow component. They concluded that their model provided reasonable predictions for a test area with non-complex terrain but performed poorly in more complex terrain due to inaccurate prediction of windspeeds.

The primary objective of this study is to use readily available landscape and stand level information to produce empirical models for prediction of windthrow damage along the boundaries of proposed harvest units. A secondary objective is to investigate the relationship between site, stand and management factors and damage severity.

Section snippets

Description of study area

The study area is located at the northern tip of Vancouver Island bounded by Johnston Strait to the East and Quatsino Sound to the West. Six 1:20 000 mapsheets, covering an area of 900 km2 were selected as representative of the range of biophysical conditions in this landscape. This area is in the Nahwitti Coastal Lowland physiographic region bordering on the northern extension of the Vancouver Island Ranges to the south and west (Lewis, 1985). Low pressure systems originating over the Pacific

Results

Several of the stand variables are correlated, reflecting linear relations between stand age and stocking, and height of mature stands and site quality (Table 6). Ground slope increases with elevation, and topex scores increase in areas with steeper slopes. The more mountainous condition of the west coast is reflected in the correlation of coast with elevation and TOPEX2000.

The major stand types are also associated with specific site features. The χ2 tests indicated that the distribution of

Discussion

Boundary orientation relative to damaging winds is commonly associated with higher windthrow risk (e.g. Ruel, 1995). There is a clear directionality to the damage in this area of Northern Vancouver Island, with boundary orientation and position of segment relative to opening centroid both contributing to increased damage. The vulnerability of boundaries is consistent with the storm wind directions recorded at the local AES stations with a more easterly component on the east side of the island

Conclusion

Windthrow results from complex interactions of climatic, site and management factors. GISs combined with aerial photograph detection of wind damage provide an opportunity to develop windthrow risk models and maps for large landscape units. The large datasets produced in this approach enable detection of the association of damage with site and management features which is not easy to determine from field observations of a limited number of sites. By using endemic damage associated with cutblock

Acknowledgements

Funding for this project was provided by Forest Renewal British Columbia. Western Forest Products Limited provided logistical support, GIS and aerial photograph information. Arnold Moy and Eric Aubert assisted with programming. Jerry Maedel of the UBC GIS/FIRMS lab provided advice and assistance throughout the project.

References (32)

  • R.R Alexander

    Minimizing windfall around clear cuttings in spruce–fir forests

    For. Sci.

    (1964)
  • C.J Anderson et al.

    Root extraction force measurements for Sitka spruce

    Forestry

    (1989)
  • Andrus, C., Froehlich, H., 1992. Wind damage in streamside buffers and its effect on accelerated sedimentation in...
  • J Becquey et al.

    L’existence de zones de stabilite des peuplements. Consequences sur la gestion

    Revue Forestiere Francaise

    (1987)
  • K.W Cremer et al.

    Effects of stocking and thinning on wind damage in plantations

    NZ J. For. Sci.

    (1982)
  • A.E Elling et al.

    Predicting wind-caused mortality in strip cut stands of peatland black spruce

    For. Chron.

    (1978)
  • Environment Canada, 1994. Canadian Monthly Climate Data and 1961–1990 Climate Normals. CDROM. Environment Canada...
  • ESRI, 1996. Using ArcView GIS. Environmental Systems Research International Inc., Redlands, CA, 340...
  • Fleming, R.L., Crossfield, R.M., 1983. Strip cutting in shallow-soil upland black spruce near Nipigon, Ontario. III....
  • J Fridman et al.

    Modelling probability of snow and wind damage using tree, stand, and site characteristics from Pinus sylvestris sample plots

    Scand. J. For. Res.

    (1998)
  • B.A Gardiner et al.

    Field and wind tunnel assessments of the implications of respacing and thinning for tree stability

    Forestry

    (1997)
  • Hannah, P., Palutikof, J.P., Quine, C.P., 1995. Predicting windspeeds for forest areas in complex terrain. In: Coutts,...
  • Harris, A.S., 1989. Wind in the forests of southeast Alaska and guides for reducing damage. USDA Forest Service...
  • Hosmer, D.W., Lemeshow, S., 1989. Applied Logistic Regression. Wiley, Toronto, Ont., 307...
  • O Laiho

    Susceptibility of forest stands to windthrow in Southern Finland

    Folia Forestalia

    (1987)
  • Lewis, T., 1985. Ecosystems of the Quatsino Tree Farm Licence (T.F.L. 6). Unpublished Report prepared for Western...
  • Cited by (97)

    • Predicting the risk of tree fall onto railway lines

      2024, Forest Ecology and Management
    • Long-term cumulative impacts of windthrow and subsequent management on tree species composition and aboveground biomass: A simulation study considering regeneration on downed logs

      2021, Forest Ecology and Management
      Citation Excerpt :

      Windthrow risk is influenced not only by stand age but also by the wind direction, wind intensity, forest structures, topography, and so on. Previous studies have reported that the risk of windthrow is greater in individual trees with larger diameters at breast height or tree heights (Rich et al., 2007) and in stands that have homogeneous structures (Jalkanen and Mattila, 2000; Mitchell et al., 2001; Morimoto et al., 2019b). On the other hand, predicting windthrow risk is very challenging because it varies greatly among regions and among windthrow events (Dobor et al., 2019).

    • An evidence-based review of the effectiveness of riparian buffers to maintain stream temperature and stream-associated amphibian populations in the Pacific Northwest of Canada and the United States

      2021, Forest Ecology and Management
      Citation Excerpt :

      Wind damage is strongly associated with buffer orientation relative to the predominant storm direction (i.e., southeast, south, southwest in the Pacific Northwest) and local conditions (including wind fetch length resulting from the size of clearcut units; Kramer et al., 2001; Mitchell et al., 2001; Beese et al., 2019). Reductions in windthrow mortality are feasible when site and landscape factors are considered in harvest unit plans (Kramer et al., 2001; Mitchell et al., 2001; Beese et al., 2019). The temperature regimes of headwater streams were spatially and temporally variable.

    • Risk assessment of forest disturbance by typhoons with heavy precipitation in northern Japan

      2021, Forest Ecology and Management
      Citation Excerpt :

      The estimated value of the probability by modelling did not indicate the actual probability of forest disturbance at this event but rather its relative probability. As for explanatory variables, we first selected the crucial factors identified by previous studies focused on windthrow risk assessments in mountainous regions dominated by cool-temperate and boreal forests (Kramer et al., 2001; Mitchell et al., 2001; Nakajima et al., 2009). In addition, we selected a series of precipitation indicators identified by a previous study focused on shallow landslides in forests caused by heavy rainfall (Dou et al., 2019).

    View all citing articles on Scopus
    View full text