Difference between revisions of "Climate Model"

Latest revision as of 09:07, 7 December 2022

[DRAFT]. Last updated 2016-03-23.

Introduction

In Heureka, the effect of global warming on forest growth can be accounted for based on results from a process-based model called BIOMASS. BIOMASS was developed by McMurtrie et. al (1990)^[1] and has been modified and validated for Swedish conditions in a number of studies ^[2] ^[3] ^[4] ^[5]. BIOMASS is computationally demanding to run and not suitable to include directly in Heureka’s growth simulator. Therefore BIOMASS has been run for different parts of the country, forest conditions and climate scenarios, and the results have been used to construct an approximation model of BIOMASS.

How the model works

A special approach has also been developed to adapt the response of the growth models used in Heureka. The principal idea is to augment or shift the growth function when climate changes. This is done by changing tree ages; the size to age relation of trees is a critical variable in all growth functions used in Heureka. In the system, we differ between actual tree age and biological tree age. The biological age is the one subject to adjustment.

The climate response model affects the following variables directly:

The biological age or basal area depending on user setting
Site index (SIS)
Vegetation index
Temperature sum

Main calculation steps

The selected climate scenario has information on relative growth difference (α, see below) between an unchanged climate and a changed climate. This is called the growht response.
At a given year t, Heureka calculates the growth from year t to t+5, for an unchanged climate (ΔG_empir).
Calculate G_clim = α · ΔG_empir = The expected growth when there is a climate change.
A new growth calculation is done, with climate-adjusted values for site index and vegetation index (see below)
Depending on whether tree ages should be adjusted or not (see control table Climate Model, one of the following is applied:
Age adjustment should be made

By using binary search, all tree ages are adjusted with a certain factor c, until the obtained growth equals G_clim (within a certain tolerance).

If no age adjustment should be made

The growth projection length is extended (or shortened if the climate model predicts decreased growth due to drought) until the expected growth is obtained. Non-linear interpolation is used to handle that this will unlikely coincide with an exact five-year interval. All tree attributes in at time t+5 (diameter, height etc) except tree ages are set to the values obtained in the extended (or shortened period).

Calculation of growth response (changed growth due to climate change)

Step 1: Uncorrected growth response ( $\alpha_{BIOMASS})$

A "climate scenario" file imported to Heureka is not an actually scenario but a table of coeffiicients used to calculate growth modification factors that are then used in the Heureka "climate model".

For a given prediction unit that has soil moisture class j, the growth correction factor is calculated as function of the parameters a, b and c, which depend on scenario, geographic locations, tree species and soil moisture code.

$\alpha_{BIOMASS}(s) = a(s,j) LAI(s)^2 + b(s,j) LAI(s) + c(s,j)$

where

α_BIOMASS(i,j) = Growth correction factor for species j according to model BIOMASS, and
LAI(s) = Leaf area index for species s (m² leaf area / m² forest floor area)

Since leaf area is not calculated in Heureka, foliage biomass is used instead, and multiplied with a conversion factor (SLA) to get LAI.

$\displaystyle LAI(s) = \frac{fbm(s) SLA(s)}{f(s)}$

where
fbm(s) = Foliage biomass (kg/m², dry matter / forest floor area), and
SLA(s) = Conversion factor for biomass foliage to leaf area for species s (m² projected one-side leaf area / kg dry matter),
f(s) = Species proportion of species s in the stand, used as indicator for how much of the forest floor area is occupied by species s. The division with the species distribution is done because the divisor (forest floor area) should only include the forest floor occupied by the subject trees.

Step 2: Corrected response (β)

$\alpha_{BIOMASS}$ reflects "optimal" stand conditions and should therefore be modified. The response modifier is a linear function of vegetation index.

$\displaystyle \beta(s) = c_0 +(VIX-VIX_{min}(s))\frac{1-c_0(s)}{VIX_{max}(s)-VIX_{min}(s)}$

where

β(s) = Response modifier for species s, and
VIX = Vegetation index, and

c₀(s) = Intercept for species s (see control table), and
VIX_min = Minimum vegetation index for species s, and
VIX_max = Maximum vegetation index for species s.

The modified response for each species s is:

$\alpha_{c1}(s)=1 + \beta(s)(\alpha_{BIOMASS}(s)-1)$

Step 3: Carbon dioxide response

Step 1 and 2 can either included the combined effect or water, temperature and carbon dioxide, or add the carbon dioxide component in a third step. In that case, other coefficients are used in step 1. The carbon dioxide response modifier is calculated with the same type of formula as α(s), but with other coefficients.

$\displaystyle \gamma(s) = d(s,j) LAI(s)^2 + e(s,j) LAI(s) + f(s,j)$

This value is then multiplied with the calculated growth response.

$\alpha_{c2}(s)= \gamma(s) \alpha_{c1}(s)$

Temperature sum change

The change in temperature over time caused by a changing climate is expressed as an average temperature sum change factor (λ) for five years, and is part of the climate scenario input data.

Site index change

The site index (sis) change from one five-year period to the next is calculated as

$\Delta sis = \lambda (21.9947 \Delta TS - 2.1384 \Delta TS2)$

where

$\lambda$ = Temperature change factor for a five-year period
$TS_{t} = \lambda \cdot TS_{t-1}$ (temperature sum in °C)
$\Delta TS$ = Change in temperature sum (1000 °C) between two five-year periods: $0.001\ (TS_{t} - TS_{t-1})$ ) where t is period index, and
$\Delta TS2$ = Change between squared temperature sums (in 1000 °C) from period t-1 to period t: $(0.001 TS_{t})^2 - (0.001 TS_{t-1})^2$

Vegetation index change

The vegetation index is a mapping of the categorial variable vegetation type code, to a nominal variable. The vegetation index is used as one of the explanatory in the Whole-stand growth model (see page 39 in the "Growth modelling in Heureka" document).

The climate model modifies the vegetation index. The vegetation index increases when the climate model predicts increased growth. The vegetation index is updated in each projection period as a linear function of ΔTS2:

$\Delta VIX = 1.24 \cdot \Delta TS2 \cdot \frac{periodLength}{5}$

Growth model

Case 1: Whole-stand growth model

In Heureka a stand-level growth model is used (as default) to calculate basal area growth per ha on each plot, in combination with a single-tree growth model that distributes the growth to individual tree objects. Mean age is used as one of the explanatory variables. Different species are not differentiated, and it is not possible to make an adjustment for each species separately. Therefore, the response for all species together is weighted proportionally to the species proportions:

$\displaystyle \alpha_{tot} = \sum_{s=1}^{m} f(s) \alpha_{c2}(s)$

where

m = number of species, and
f(p) = volume proportion of species s

Case 2: Single-tree growth model

Single-tree or other growth model that used species-wise ages as function parameter. In this case, the growth adjustment is differentiated with respect to species (binary search can be applied to each species group separately when adjusting tree ages), using indexed \alpha_{c2}(s) instead of alpha_{tot}.

Input data

With the installation of Heureka, there are three climate model scenarios available:

MPI 4.5: Based on Max Planck Institute MPI-ESM model using radiation scenario RCP 4.5, which assumes that radiative forcing stabilises at 4.5 W/m² before the year 2100. See also RCP4.5 at SMHI's homepage.

MPI 8.5: Based on Max Planck Institute MPI-ESM model using radiation scenario RCP 4.5, which assumes that radiative forcing stabilises at 8.5 W/m² before the year 2100.. See also RCP8.5 at SMHI's homepage

ECHAMS_A1B: Based on Max Planck Institute climate model ECHAM using emission scenario SRES A1B.

About the climate scenarios

About SMHI climate scenarios

Climate scenarios used in Heureka in the SKA 15 project

Import climate scenario in Heureka Helpdoc.

Model settings that a user can modify

See Climate Model control table.

References

↑ McMurtrie R. E., Rook D. A. & Kelliher F. M. 1990. Modelling the yield of Pinus radiate on a site limited by water and nitrogen. Forest Ecology and Management 30: 381–413.
↑ Bergh J., McMurtrie R. E. & Linder S. 1998. Climatic factors controlling the productivity of Norway spruce: a modelbased analysis. Forest Ecology and Management 110: 127–139.
↑ Freeman M & Linder S. 2001. Boreal forests. In: Long-term effects of climate change on carbon budgets of forests in Europe (eds. Kramer, K. & Mohren, G.M.J.) pp. 197–203. Alterra-report 194. Alterra, Green World Research, Wageningen, 2001
↑ Bergh J., Freeman M., Sigurdsson B. D., Kellomäki S., Laitinen K., Niinistö S., Peltola, H. & Linder S. 2003. Modelling the shortterm effects of climate change on the productivity of selected tree species in Nordic countries. Forest Ecology and Management 183:327–340
↑ Freeman M, Morén A-S, Strömgren M & Linder S. 2005. Climate Change Impacts on Forests in Europe: Biological Impact Mechanisms. In: Management of European Forests under Changing Climatic Conditions (eds. Kellomäki, S. and Leinonen, S.). ResearchNotes 163, University of Joensuu, Forest Faculty, pp. 45-115. ISBN 952-458-652-5

[BIOMASS_ref-1] McMurtrie R. E., Rook D. A. & Kelliher F. M. 1990. Modelling the yield of Pinus radiate on a site limited by water and nitrogen. Forest Ecology and Management 30: 381–413.

[2] Bergh J., McMurtrie R. E. & Linder S. 1998. Climatic factors controlling the productivity of Norway spruce: a modelbased analysis. Forest Ecology and Management 110: 127–139.

[3] Freeman M & Linder S. 2001. Boreal forests. In: Long-term effects of climate change on carbon budgets of forests in Europe (eds. Kramer, K. & Mohren, G.M.J.) pp. 197–203. Alterra-report 194. Alterra, Green World Research, Wageningen, 2001

[4] Bergh J., Freeman M., Sigurdsson B. D., Kellomäki S., Laitinen K., Niinistö S., Peltola, H. & Linder S. 2003. Modelling the shortterm effects of climate change on the productivity of selected tree species in Nordic countries. Forest Ecology and Management 183:327–340

[5] Freeman M, Morén A-S, Strömgren M & Linder S. 2005. Climate Change Impacts on Forests in Europe: Biological Impact Mechanisms. In: Management of European Forests under Changing Climatic Conditions (eds. Kellomäki, S. and Leinonen, S.). ResearchNotes 163, University of Joensuu, Forest Faculty, pp. 45-115. ISBN 952-458-652-5

[1]

[2]

[3]

[4]

[5]

@@ Line 1: / Line 1: @@
 [[Category:Model]]
-TODO.
+[DRAFT]. Last updated 2016-03-23.
 ==Introduction==
@@ Line 6: / Line 6: @@
 <ref>Bergh J., McMurtrie R. E. & Linder S. 1998. Climatic factors controlling the
   productivity of Norway spruce: a modelbased analysis. Forest Ecology and
-  Management 110: 127–139.</ref>,
+  Management 110: 127–139.</ref>
 <ref>Freeman M & Linder S. 2001. Boreal forests. In: Long-term effects of climate change
   on carbon budgets of forests in Europe (eds. Kramer, K. & Mohren, G.M.J.) pp.
-–203. Alterra-report 194. Alterra, Green World Research, Wageningen, 2001</ref>,
+–203. Alterra-report 194. Alterra, Green World Research, Wageningen, 2001</ref>
 <ref>Bergh J., Freeman M., Sigurdsson B. D., Kellomäki S., Laitinen K., Niinistö S., Peltola, H. & Linder S. 2003. Modelling the shortterm effects of climate change on the
   productivity of selected tree species in Nordic countries. Forest Ecology and
-  Management 183:327–340</ref>,
+  Management 183:327–340</ref>
 <ref>Freeman M, Morén A-S, Strömgren M & Linder S. 2005. Climate Change Impacts on
   Forests in Europe: Biological Impact Mechanisms. In: Management of European
@@ Line 18: / Line 18: @@
   ResearchNotes 163, University of Joensuu, Forest Faculty, pp. 45-115.
   <nowiki>ISBN 952-458-652-5</nowiki></ref>.
-BIOMASS is computationally demanding to run and not suitable to include directly in Heureka’s growth simulator. Therefore BIOMASS was run for different parts of the country and different forest conditions, and the results have been used to construct an approximation model of BIOMASS.
+BIOMASS is computationally demanding to run and not suitable to include directly in Heureka’s growth simulator. Therefore BIOMASS has been run for different parts of the country, forest conditions and climate scenarios, and the results have been used to construct an approximation model of BIOMASS.
 ==How the model works==
@@ Line 24: / Line 24: @@
 The climate response model affects the following variables directly:
-*The biological age or basal area (gt) depending on user setting
+*The biological age or basal area depending on [[ControlTable_Climate_Model#Apply Age Adjustment| user setting]]
 *Site index (SIS)
 *Vegetation index
 *Temperature sum
+====Main calculation steps====
+#The selected climate scenario has information on relative growth difference (&alpha;, see below) between an unchanged climate and a changed climate. This is called the growht response.
+#At a given year ''t'', Heureka calculates the growth from year ''t'' to ''t''+5, for an unchanged climate (&Delta;G<sub>empir</sub>).
+#Calculate G<sub>clim</sub> = &alpha; &middot; &Delta;G<sub>empir</sub> = The expected growth when there is a climate change.
+#A new growth calculation is done, with climate-adjusted values for site index and vegetation index (see below)
+#Depending on [[ControlTable_Climate_Model#Apply Age Adjustment| whether tree ages]] should be adjusted or not (see control table [[ControlTable_Climate_Model#Apply Age Adjustment|Climate Model]], one of the following is applied:
+#;[[ControlTable_Climate_Model#Apply Age Adjustment|Age adjustment]] should be made: By using binary search, all tree ages are adjusted with a certain factor c, until the obtained growth equals G<sub>clim</sub> (within a certain tolerance).
+#;If no age adjustment should be made: The growth projection length is extended (or shortened if the climate model predicts decreased growth due to drought) until the expected growth is obtained. Non-linear interpolation is used to handle that this will unlikely coincide with an exact five-year interval. All tree attributes in at time ''t''+5 (diameter, height etc) except tree ages are set to the values obtained in the extended (or shortened period).
+==Calculation of growth response (changed growth due to climate change)==
 ====Step 1: Uncorrected growth response (<math>\alpha_{BIOMASS})</math>====
 A "climate scenario" file imported to Heureka is not an actually scenario but a table of coeffiicients used to calculate growth modification factors that are then used in the Heureka "climate model".
@@ Line 76: / Line 86: @@
 <math>\alpha_{c2}(s)= \gamma(s) \alpha_{c1}(s)</math>
-====Temperature sum change====
+==Temperature sum change==
 The change in temperature over time caused by a changing climate is expressed as an average temperature sum change factor (&lambda;) for five years, and is part of the climate scenario input data.
-====Site index change====
+==Site index change==
 The site index (sis) change from one five-year period to the next is calculated as
@@ Line 87: / Line 97: @@
 <math>\lambda</math> = Temperature change factor for a five-year period
-<br><math>TS_{t} = \lambda \cdot TS_{t-1}</math>
+<br><math>TS_{t} = \lambda \cdot TS_{t-1}</math> (temperature sum in &deg;C)
-<br><math>\Delta TS</math> = Change in temperature sum between two five-year periods (<math>TS_{t} - TS_{t-1}</math>) and t is period index, and
+<br><math>\Delta TS</math> = Change in temperature sum (1000 &deg;C) between two five-year periods<tt>:</tt> <math>0.001\ (TS_{t} - TS_{t-1})</math>) where ''t'' is period index, and
-<br><math>\Delta TS2</math> = Change between squared temperature sums from period t-1 to period t (<math>TS_{t}^2 - TS_{t-1}^2</math>)
+<br><math>\Delta TS2</math> = Change between squared temperature sums (in 1000 &deg;C) from period ''t''-1 to period ''t''<tt>:</tt> <math>(0.001 TS_{t})^2 - (0.001 TS_{t-1})^2</math>
-====Vegetation index change====
+==Vegetation index change==
 The vegetation index is a mapping of the categorial variable [[Definition:VegetationTypeCode | vegetation type code]], to a nominal variable. The vegetation index is used as one of the explanatory in the Whole-stand growth model (see page 39 in the [http://heurekaslu.org/mw/images/9/93/Heureka_prognossystem_(Elfving_rapportutkast).pdf "Growth modelling in Heureka" document]).
-The climate model modifies the vegetation index. The vegetation index increases when the climate model predicts increased growth. The vegetation index is updated in each projection period.
+The climate model modifies the vegetation index. The vegetation index increases when the climate model predicts increased growth. The vegetation index is updated in each projection period as a linear function of &Delta;TS2<tt>:</tt>
+<math>\Delta VIX = 1.24 \cdot \Delta TS2 \cdot \frac{periodLength}{5}</math>
+==Growth model==
 ====Case 1: Whole-stand growth model====
 In Heureka a stand-level growth model is used (as default) to calculate basal area growth per ha on each plot, in combination with a single-tree growth model that distributes the growth to individual tree objects. Mean age is used as one of the explanatory variables. Different species are not differentiated, and it is not possible to make an adjustment for each species separately. Therefore, the response for all species together is weighted proportionally to the species proportions: <br>
-<math> \displaystyle \alpha_{tot} = \sum_{s=1}^{m} f(s) \alpha_c(s)</math>
+<math> \displaystyle \alpha_{tot} = \sum_{s=1}^{m} f(s) \alpha_{c2}(s)</math>
 where
@@ Line 107: / Line 120: @@
 ====Case 2: Single-tree growth model====
-In case the whole-stand growth model is not used
+Single-tree or other growth model that used species-wise ages as function parameter. In this case, the growth adjustment is differentiated with respect to species (binary search can be applied to each species group separately when adjusting tree ages), using indexed \alpha_{c2}(s) instead of alpha_{tot}.
 ==Input data==
 With the installation of Heureka, there are three climate model scenarios available:
-:'''MPI 4.5''': Based on Max Planck Institute [http://www.mpimet.mpg.de/en/science/models/mpi-esm/ MPI-ESM model] using radiation scenario RCP 4.5, which assumes that radiative forcing stabilises at 4.5 W/m² before the year 2100. See also [http://www.smhi.se/en/climate/climate-scenarios#area=eur&dnr=0&sc=rcp45&seas=ar&var=t RCP4.5 at SMHI's homepage].
+:'''MPI 4.5''': Based on Max Planck Institute [http://www.mpimet.mpg.de/en/science/models/mpi-esm/ MPI-ESM model] using radiation scenario RCP 4.5, which assumes that radiative forcing stabilises at 4.5 W/m² before the year 2100. See also [https://www.smhi.se/en/climate/future-climate/basic-climate-change-scenario-service/sverige/medeltemperatur/rcp45/2071-2100 RCP4.5 at SMHI's homepage].
-:'''MPI 8.5''': Based on Max Planck Institute [http://www.mpimet.mpg.de/en/science/models/mpi-esm/ MPI-ESM model] using radiation scenario RCP 4.5, which assumes that radiative forcing stabilises at 8.5 W/m² before the year 2100.. See also [http://www.smhi.se/en/climate/climate-scenarios#area=eur&dnr=0&sc=rcp85&seas=ar&var=t RCP8.5 at SMHI's homepage]
+:'''MPI 8.5''': Based on Max Planck Institute [http://www.mpimet.mpg.de/en/science/models/mpi-esm/ MPI-ESM model] using radiation scenario RCP 4.5, which assumes that radiative forcing stabilises at 8.5 W/m² before the year 2100.. See also [https://www.smhi.se/en/climate/future-climate/basic-climate-change-scenario-service/sverige/medeltemperatur/rcp85/2071-2100 RCP8.5 at SMHI's homepage]
 :'''ECHAMS_A1B''': Based on Max Planck Institute climate model [http://www.mpimet.mpg.de/en/science/models/echam.html ECHAM] using emission scenario SRES A1B.
-====See also====
+====About the climate scenarios====
 :[http://www.smhi.se/en/climate/climate-scenarios/haag_en.html About SMHI climate scenarios]
 :[http://shop.skogsstyrelsen.se/shop/9098/art35/31485335-a7bfe8-Klimat_webb.pdf Climate scenarios used in Heureka in the SKA 15 project]
 :[http://heurekaslu.org/help/en/index.html?importera_klimatscenario.htm Import climate scenario in Heureka Helpdoc].
 ==Model settings that a user can modify==
-====Selected climate scenario====
+See [[ControlTable_Climate_Model|Climate Model control table]].
-====Model parameters====
-Several model parameters can be modified via the [[ControlTable_Climate_Model|Climate Model control table]].
-==Age adjustments==
 ==References==
 <references />
--->

Difference between revisions of "Climate Model"

Latest revision as of 09:07, 7 December 2022

Contents

Introduction

How the model works

Main calculation steps

Calculation of growth response (changed growth due to climate change)

Step 1: Uncorrected growth response ( $\alpha_{BIOMASS})$

Step 2: Corrected response (β)

Step 3: Carbon dioxide response

Temperature sum change

Site index change

Vegetation index change

Growth model

Case 1: Whole-stand growth model

Case 2: Single-tree growth model

Input data

About the climate scenarios

Model settings that a user can modify

References

Navigation menu

Search