Dead wood modelling at stand-level

Transcription

1 Dead wood modelling at stand-level Jari Hynynen & Harri Mäkinen Finnish Forest Research Institute Vantaa Research Unit

2 MOTTI stand simulator Salminen et al. (2005), Hynynen et al. (2005) A stand-level analysis tool for assessing the effects of alternative forest management on development of stand and trees profitability of forest management biomass production and carbon sequestration forest biodiversity A workbench for modellers a tool for model building and testing A stand-level decision support system

3 User-defined input Stand inventory data MOTTI-system Stand simulation State update Growth prediction Management options Statistical model Process-based model Regeneration - Growth - Mortality Timber assortment criteria Economical parameters Analysis Yield Wood quality Forest economics Biomass & Carbon Biodiversity Request for reports Output of the results

4 Prediction of dead wood dynamics in Motti Prediction of mortality tree-level: survival models stand-level: self-thinning The amount of produced dead wood Decomposition of dead wood in terms of stem volume, density, mass The amount of actual dead wood Assessment of the quality of dead wood diversity index The "value" of dead wood

5 Survival models for individual trees

6 pcomp5 = Survival probability of a tree 1+ exp Mortality due to competition Haapala 1983, Hynynen et al ( ) a0 + a1 d + a2 BA + a3 BAL p comp5 d = probability of a tree to die during the next five years = breast height diameter, cm BA = stand basal area, m 2 ha -1 BAL = basal area of trees larger than subject tree, m 2 ha -1 P comp Basal area, m 2 ha d, cm 25

7 Mortality due to aging of a tree Hynynen et al Survival probability of a tree pold( age) = e 1+ e 10 age A max 10 age A Maximum tree age (A max ) max Survival probability, % Tree species Temperature sum, ddy Maximum age, years Scots pine (Pinus sylvestris) Norway spruce (Picea abies) Silver birch (Betula pendula) Pubescent birch (Betula pubescens) Aspen (Populus tremula) Grey alder (Alnus incana) Red alder (Alnus glutinosa) Mountain ash (Sorbus aucuparia) Goat-willow (Salix caprea) Other coniferous sp Other broadleaved sp FINNISH 0.5 FOREST RESEARCH 1.0 INSTITUTE 1.5 age/a max

8 Stand-level model for self-thinning

9 Model for self-thinning Hynynen 1993, Hynynen et al ( ) 0 2 ( ) 1 ( ) ln N = a + a ln SI + a ln D + b + e ij i gkij i ij N = stem number, trees/ha SI = site index, m D = mean diameter, cm Stem number, n ha Scots pine Norway spruce Silver birch Pubescent birch Mean diameter at stump height, FINNISH cm FOREST RESEARCH INSTITUTE

10 Prediction of mortality in MOTTI Survival on tree-level pcomp5 p comp5 d = 1+ exp 1 ( ) a0 + a1 d + a2 BA + a3 BAL Self-thinning on stand-level ( ) ( ) ( ) = probability of a tree to die during the next five years ln Nij = a0+ a2ln SIi + a1ln Dgkij + bi + eij = breast height diameter, cm BA = stand basal area, m 2 ha -1 BAL = basal area of trees larger than subject tree, m 2 ha -1 P comp Basal area, m 2 ha Stem number, n ha IF Scots pine Norway spruce Siver birch Pubescent birch d, cm Mean diameter at FINNISH stumpforest height, RESEARCH cm INSTITUTE

11 Models for decomposition of stems

12 Modelling data Data from permanent sample plots Scots pine (36), Norway spruce (13), birch (9) Stand treatment history is known Mortality of trees Time of death is known size of trees at the time of death is known Field measurement both living and dead trees were measured sample discs from stems Laboratory analysis wood density estimation volume, density and mass of stems Tree species Scots pine Norway spruce Birch

13 Models for dead wood dynamics Models for Scots pine, Norway spruce and birch probability that a tree remains standing as a snag after its death Models for the remaining fraction of stem volume stem density stem mass as a function of time since the death of a tree

14 Models for the probability that a tree remains standing as a snag after its death Mäkinen et al. (2006) Pr( snag ln 1 Pr( snag spti ) ) spti = g 1 + g 4 y spti + g 7 dbh spt0 + ε + ε 1s 1spt Probability Scots pine Norway spruce Birch Time since death, years

15 Models for the remaining fraction of stem volume Mäkinen et al. (2006) v v spti spt0 = exp ( exp( c + c y + c dbh + γ ) 1 4 spti 7 spt0 1spt Remaining fraction Scots pine Norway spruce Birch Time since death, years

16 Models for the remaining fraction of wood density Mäkinen et al. (2006) r r spti spt0 = exp ( exp( b + b y + b dbh + β ) 1 4 spti 7 spt0 1spt Remaining fraction Scots pine Norway spruce Birch Time since death, years

17 Models for the remaining fraction of stem mass Mäkinen et al. (2006) m m spti spt0 = exp ( exp( d + d y + d dbh + δ ) 1 4 spti 7 spt0 1spt Remaining fraction Scots pine Norway spruce Birch Time since death, years

18 Conclusions Dead wood decomposition models Decomposition rate depends on the time since death slow initial decomposition period: 0 => 5 to 10 years period of radip decomposition period of modearte decomposition rate Conifers stem mass decreased by 50 % in years stems competely decayed in years Birch stem mass decreased by 50 % in 10 years stems competely decayed in years The models provide a framework for predicting dead wood dynamics in managed and dense unthinned stands

19 Further analysis

20 Diversity Index of dead wood An index based on properties of dead wood tree species stem diameter (diameter classes) degree of decomposition (decay classes) type of dead stem (snag/log) Can be used to assess the "value" of dead wood to biodiversity

21 An example The effect of alternative management practices on the dead wood production and on the profitability of forest management

22 A simulation study for a Norway spruce stand (one rotation - 65 years) Management schedules forest management recommendations (Tapio) active production of dead wood in intermediate thinnings 10 m3/ha 20 m3/ha 40 m3/ha The effects management on the amount of dead wood diversity index of dead wood profitability of forest management

23 The volume of dead wood at the end of the rotation (Norway spruce, age 65 years) Amount of deadwood, m3/ha Log Snag Managed (Tapio) Unmanaged Management

24 Diversity index of dead wood A Norway spruce stand diversity index Management 40 Tapio thinning 10 m3/ha 20 m3/ha 20 precommercial thinning 1. thinning 40 m3/ha Unmanaged Stand age, years

25 relative NPV, % 120 Relative net present values lnterest rate 3 % Tapio 10 m3/ha 20 m3/ha 40 m3/ha Unmanaged Management

26 References Haapala, P Luonnonpoistuman ennustaminen puun kuolemistodennäköisyysmallilla. Metsäntutkimuslaitoksen puuntuotoksen tutkimussuunta. Stencil. 33 p. [In Finnish]. Hynynen, J Self-thinning models for even-aged stands of Pinus sylvestris, Picea abies and Betula pendula. Scandinavian Journal of Forest Research 8(3): Hynynen, J., Ojansuu, R., Hökkä, H., Siipilehto, J., Salminen, H. & Haapala, P Models for predicting stand development in MELA System. Finnish Forest Research Institute, Research Papers p. Hynynen, J., Ahtikoski, A., Siitonen, J., Sievänen, R. & Liski, J Applying the MOTTI simulator to analyse the effect of alternative management schedules on timber and non-timber production. Forest Ecology and Management 207: Mäkinen, H., Hynynen, J., Siitonen, J. & Sievänen, R Predicting the decomposition of Scots pine, Norway spruce, and birch stems in Finland. Ecological Applications 16(5): Salminen, H., Lehtonen, M. & Hynynen, J Reusing legacy FORTRAN in the MOTTI growth and yield simulator. Computers and Electronics in Agriculture 49(1):