# Diameter Increment

For every year of modeled growth, MGM predicts an annual diameter increment for each tree > 1.3m tall, relative to height increment, local competition, and/or other variables. Different diameter increment models are used for trees ≥ 4cm DBH and trees < 4cm DBH.

## Trees ≥ 4cm DBH

For trees ≥ 4cm DBH, annual diameter increment is predicted using 1 of 2 approaches:

*Compatible Diameter Increment Models*– Diameter increment models that include height increment as a variable.*Non-Compatible Diameter Increment Models*– Diameter increment models that do not include height increment as a variable.

#### Compatible Diameter Increment Models

The annual diameter increment models for lodgepole pine and black spruce include height increment and competition as variables (Equations D1 and D2). As a result, trees with a large height increment and low competition experience more diameter growth. Trees with a small height increment and high competition experience less diameter growth. These models apply to all lodgepole pine and black spruce ≥ 4cm DBH.

**Lodgepole Pine**

Where:

*DI _{ijk}* = Diameter increment (cm/year) of tree

*k*of species

*j*in stand

*i*

*HI _{ijk}* = Height increment (m/year) of tree

*k*of species

*j*in stand

*i*

*rDBALT _{ijk}* = DBALT

_{ijk}/ maxBA

_{i}

*rSwFBALT _{ijk}* = SwFBALT

_{ijk}/ maxBA

_{i}

*rPBALT _{ijk}* = PBALT

_{ijk}/ maxBA

_{i}

*DBALT _{ijk}* = Basal area (m²/ha) of deciduous trees larger than subject tree

*k*of species

*j*in stand

*i*

*SwFBALT _{ijk}* = Basal area (m²/ha) of white spruce-fir trees larger than subject tree

*k*of species

*j*in stand

*i*

*PBALT _{ijk}* = Basal area (m²/ha) of pine trees larger than subject tree

*k*of species

*j*in stand

*i*

*MaxBA _{i}* = maxDen

_{i}× AreaConvConst × QMD

_{i}

^{2}

*maxDen _{i}* = ((1 / QMD

_{i}+ 0.00865) / 0.001244)

^{ (1 / 0.5225)}

*AreaConvConst* = 7.8537E^{-05}

*QMD _{i}* = Quadratic mean diameter (cm) of stand

*i*for all species

*a _{0}* = -0.331,

*a*= 1.263,

_{1}*a*= 3.34,

_{2}*a*= 2.061

_{3}

**Black Spruce**

For black spruce, the annual diameter increment model (Equation D2) was developed in Oboite and Comeau (2021) using 4,139 plots, ranging from Alaska to Manitoba. Equation D2 considers the impact of height increment, diameter (i.e. tree size), and competition on annual diameter increment.

Where:

*DI _{ijk}* = Diameter increment (cm/year) of tree

*k*of species

*j*in stand

*i*

*HI _{ijk}* = Height increment (m/year) of tree

*k*of species

*j*in stand

*i*

*DBH _{ijk}* = Diameter at breast height (cm) of tree

*k*of species

*j*in stand

*i*

*DBALT _{ijk}* = Basal area (m²/ha) of deciduous trees larger than subject tree

*k*of species

*j*in stand

*i*

*SFBALT _{ijk}* = Basal area (m²/ha) of spruce-fir trees larger than subject tree

*k*of species

*j*in stand

*i*

*PBALT _{ijk}* = Basal area (m²/ha) of pine trees larger than subject tree

*k*of species

*j*in stand

*i*

*a _{0}* = 0.9285128,

*a*= -0.2663871,

_{1}*a*= 0.0218894,

_{2}*a*= 0.0076929,

_{3}*a*= 0.0910156,

_{4}*a*= 0.0027382

_{5}

#### Non-Compatible Diameter Increment Models

The annual diameter increment models for white spruce, trembling aspen, and jack pine include maximum diameter increment and competition (Equation D3). After defining maximum diameter increment, each tree is ranked by social status to determine competition-adjusted diameter increment. As a result, the largest diameter tree in each stand grows at the maximum diameter increment (Equation D3; RF_{ijk} = 1). Smaller diameter trees grow at a slower rate relative to competition (Equation D3; RF_{ijk} < 1). For mid-rotation trees, maximum diameter increment is defined using a relationship between maximum height increment, DBH, and a height-diameter curve (Equation D7). For old growth trees, maximum diameter increment is defined using a constant basal area increment model (Equation D9) and top height DBH. Competition is applied using species-specific reduction factors in Equations D10 to D13.

Where:

*DI _{ijk}* = Diameter increment (cm/year) of tree

*k*of species

*j*in stand

*i*

*MaxDI _{ijk}* = Maximum diameter increment (cm/year) of tree

*k*in species

*j*in stand

*i*(Equations D7 or D9)

*RF _{ijk}* = Reduction factor of tree

*k*of species

*j*in stand

*i*(Equations D10, D11, D12, or D13)

**Mid-Rotation Maximum Diameter Increment**

For mid-rotation trees, maximum height increment, maximum diameter increment, and DBH are linked using Equation D4 and a Weibull height-diameter curve (Equation D6). Then, Equation D4 is rearranged to solve for maximum diameter increment by multiplying the maximum potential height increment and the slope of the height-diameter curve (Equation D7). The Weibull height-diameter curve uses the parameters in Table D1.

Where:

*MaxHI _{ijk}* = Maximum potential height increment (m/year) of tree

*k*of species

*j*in stand

*i*

*MaxDI _{ijk}* = Maximum diameter increment (cm/year) of tree

*k*of species

*j*in stand

*i*

*HtDbhCurve* = Weibull height-diameter curve in Equation D6

*HD _{ijk}* = Height (m) for conifer species or diameter at breast height (cm) for deciduous species of tree

*k*of species

*j*in stand

*i*

*SiteHeight _{j}* = Height (m) defined by the provincial site index curve for species

*j*, given

*SiteIndex*and

_{ij}*BHage*

_{ijk}*SiteIndex _{ij}* = Site index (m@50 years breast height age) of species

*j*in stand

*i*

*BHage _{ijk}* = Breast height age (years) of tree

*k*of species

*j*in stand

*i*

*PH _{ijk}* = Predicted height (m) of tree

*k*of species

*j*in stand

*i*

*a _{0} – a_{2}* = Parameter estimates in Table D1.

Table D1. Parameters for the Weibull height-diameter curve (Equation D6) by species, regional variant, and subregion. The white spruce and trembling aspen parameters are from height-diameter curves fit to PSP top height trees. The jack pine parameters are from the “grand mean model” in Huang et al. (2009); Huang et al. (2009) parameter a_{2} is expressed as 1 / x to match the model form of Equation D6.

Species |
Regional Variant |
Subregion |
a_{0} |
a_{1} |
a_{2} |

White Spruce | Alberta | Montaine, Lower Foothills, Foothills Parkland | 28.54 | 1.91 | 24.34 |

Alpine, Subalpine, Upper Foothills | 25.72 | 1.68 | 24.93 | ||

Other Subregions | 25.86 | 2.23 | 21.47 | ||

British Columbia, Saskatchewan, Manitoba | All | 25.86 | 2.23 | 21.47 | |

Trembling Aspen | All | All | 34.74 | 1.14 | 29.50 |

Jack Pine | All | All | 23.54 | 1.23 | 14.36 |

**Old Growth Maximum Diameter Increment**

For old growth trees, maximum diameter increment cannot be solved using mid-rotation Equation D7. In old growth trees, maximum diameter increment approaches zero as trees approach the height-diameter asymptote (Equation D7; *MaxDI _{ijk}* → 0 as

*f ‘ (HtDbhCurve(DBH*→ 0). However, large old growth trees may persist and accrue diameter increment in late life.

_{ijk})To address this issue, a constant basal area increment (BAI) model (Equation D8) and top height DBH are used to define maximum diameter increment for old growth trees (Equation D9). Old growth trees are taller than > 80% of the height-diameter asymptote. Tree-level variables are assigned relative to a hypothetical site tree at 80% of the height-diameter asymptote. Under this approach, a constant BAI is applied to an increasing stem perimeter. This produces a maximum diameter increment that is > 0 and decreases with increasing tree diameter.

Where:

*BAI _{80%}* = Basal area increment (cm²/year) for a hypothetical tree at 80% of the height-diameter asymptote

*DBH _{80%}* = Diameter at breast height (cm) for a hypothetical tree at 80% of the height-diameter asymptote; Solved using Equation D6; If

*DBH*>

_{80%}*TopDBH*then

_{ij}*DBH*=

_{80%}*TopDBH*.

_{ij}*TopHtDBH _{ij}* = Mean top height DBH (cm) of stand

*i*of species

*j*

*MaxDI _{80%}* = Maximum diameter increment (cm/year) for a hypothetical top height tree at 80% of the height-diameter asymptote; Solved using Equation D7 and DBH

_{80%}.

**White Spruce Reduction Factors**

White spruce uses a mid-rotation reduction factor (Equation D10) and an old growth reduction factor (Equation D11). Old growth trees are taller than > 80% of the height-diameter asymptote in Equation D6.

Where:

*RF _{ijk}* = Reduction factor of tree

*k*of species

*j*in stand

*i*

*DBH* = Diameter at breast height (cm)

*TF* = Tree factor (trees/ha)

*DDLT _{ijk}* = Sum of DBH × TF values for deciduous trees larger than subject tree

*k*of species

*j*in stand

*i*

*PDLT _{ijk}* = Sum of DBH × TF values for pine trees larger than subject tree

*k*of species

*j*in stand

*i*

*SwFDLT _{ijk}* = Sum of DBH × TF values for white spruce-fir trees larger than subject tree

*k*of species

*j*in stand

*i*

*a _{0}* = -0.15,

*a*= -0.39,

_{1}*a*= -0.159

_{2}*b _{0}* = -0.07449,

*b*= -0.41993,

_{1}*b*= -0.07046

_{2}

**Trembling Aspen Reduction Factor**

Trembling aspen uses one reduction factor for mid-rotation and old growth (Equation D12). Old growth trees are taller than > 80% of the height-diameter asymptote in Equation D6.

Where:

*RF _{ijk}* = Reduction factor of tree

*k*of species

*j*in stand

*i*

*DBH* = Diameter at breast height (cm)

*TF* = Tree factor (trees/ha)

*DDLT _{ijk}* = Sum of DBH × TF values for deciduous trees larger than subject tree

*k*of species

*j*in stand

*i*

*PDLT _{ijk}* = Sum of DBH × TF values for pine trees larger than subject tree

*k*of species

*j*in stand

*i*

*SwFDLT _{ijk}* = Sum of DBH × TF values for white spruce-fir trees larger than subject tree

*k*of species

*j*in stand

*i*

*a _{0}* = 0.00002854,

*a*= 0.00010712,

_{1}*a*= 0.00001291

_{2}

**Jack Pine Reduction Factor**

Jack pine uses one reduction factor for mid-rotation and old growth (Equation D13). Old growth trees are taller than > 80% of the height-diameter asymptote in Equation D6. The jack pine reduction factor (Equation H13) was developed in Strimbu et al. (2017) using 422 plots across Alberta, Saskatchewan, and Manitoba.

Where:

*RF _{ijk}* = Reduction factor of tree

*k*of species

*j*in stand

*i*

*DBH* = Diameter at breast height (cm)

*TF* = Tree factor (trees/ha)

*DDLT _{ijk}* = Sum of DBH × TF values for deciduous trees larger than subject tree

*k*of species

*j*in stand

*i*

*PDLT _{ijk}* = Sum of DBH × TF values for pine trees larger than subject tree

*k*of species

*j*in stand

*i*

*SFDLT _{ijk}* = Sum of DBH × TF values for spruce-fir trees larger than subject tree

*k*of species

*j*in stand

*i*

*a _{0}* = 0.0005743204,

*a*= 0.0002949319,

_{1}*a*= 0.00004629732

_{2}

## Trees < 4cm DBH

For trees with a DBH < 4cm, annual diameter increment occurs under 3 cases:

*Case 1*– A tree begins the year under 1.3m tall and remains under 1.3m after the current year’s height growth.*Case 2*– A tree begins the year under 1.3m tall and exceeds 1.3m after the current year’s height growth.*Case 3*– A tree begins the year above 1.3m tall.

All trees under Case 1 receive a null DBH, identified with a “-1” in the stand worksheet. For Cases 2 and 3, annual diameter increment is modeled by MGM species group (i.e. white spruce, pine, trembling aspen, black spruce).

#### White Spruce Group

For white spruce > 1.3m tall and < 4cm DBH, annual diameter increment is modeled using Equation D14 for Case 2 and Equations D14 to D19 for Case 3. Equations D14 to D19 implement a basal area increment model. In Equation D15, the exponential adjustment for spruce-fir competition (SFBALT) is based on work by Krebs (2016), as described by Comeau and Bokalo (2020).

**Case 2**

Where:

*DBH _{ijk}* = Diameter at breast height (cm) of tree

*k*of species

*j*in stand

*i*

*H1 _{ijk}* = Height (m) of tree

*k*of species

*j*in stand

*i*after the current year’s height growth

*a _{0}* = -2.22,

*a*= 1.9877

_{1}

A minimum DBH of 0.3m is applied to any prediction < 0.3cm DBH (Equation D14).

**Case 3**

Where:

*BAI _{ijk}* = Basal area increment (cm²/year) of tree

*k*of species

*j*in stand

*i*

*DBH _{ijk}* = Diameter at breast height (cm) of tree

*k*of species

*j*in stand

*i*

*viewfactor _{ijk}* = View factor of tree

*k*of species

*j*in stand

*i*; If

*viewfactor*> 0.25 then

_{ijk}*viewfactor*= 0.25.

_{ijk}*DDenAbove _{ijk}* = Density of deciduous trees taller than (trees/ha) tree

*k*of species

*j*in stand

*i*

*DH _{ijk}* = Equation D17; If

*DH*≤ 0.1 then

_{ijk}*DH*= 0.1.

_{ijk}*MHtall10D _{i}* = Mean height (m) of the tallest 10% of deciduous trees (i.e. deciduous canopy height) in stand

*i*

*H1 _{ijk}* = Height (m) of tree

*k*of species

*j*in stand

*i*after the current year’s height growth

*SFBALT _{ijk}* = Basal area (m²/ha) of spruce-fir trees larger than subject tree

*k*of species

*j*in stand

*i*

*BAnew _{ijk}* = New basal area (cm²) of tree

*k*of species

*j*in stand

*i*for the current year

*BA _{ijk}* = Basal area (cm²) of tree

*k*of species

*j*in stand

*i*from the previous year

*DI _{ijk}* = Diameter increment (cm/year) of tree

*k*of species

*j*in stand

*i*

*a _{0}* = 1.609776,

*a*= 0.208911,

_{1}*a*= -0.0366

_{2}

#### Pine Group

For pine > 1.3m tall and < 4cm DBH, annual diameter increment is modeled using Equation D20 for Case 2 and Equations D21 to D23 for Case 3. Equations D22 and D23 apply a diameter increment reduction for high density pine stands.

**Case 2**

Where:

*DBH _{ijk}* = Diameter at breast height (cm) of tree

*k*of species

*j*in stand

*i*

*H1 _{ijk}* = Height (m) of tree

*k*of species

*j*in stand

*i*after the current year’s height growth

*a _{0}*,

*a*, and a

_{1}_{2}are subregional parameters from the Huang (1994) height-diameter model

A minimum DBH of 0.1m is applied to any prediction < 0.1cm DBH (Equation D20).

**Case 3**

Where:

*DI _{ijk}* = Diameter increment (cm/year) of tree

*k*of species

*j*in stand

*i*

*H1 _{ijk}* = Height (m) of tree

*k*of species

*j*in stand

*i*after the current year’s height growth

*H _{ijk}* = Height (m) of tree

*k*of species

*j*in stand

*i*from the previous year

*a _{0}*,

*a*, and a

_{1}_{2}are subregional parameters from the Huang (1994) height-diameter model

A minimum DBH of 0.1m is applied to any prediction < 0.1cm DBH (Equation D21).

**Case 3 – High Density Stands**

For pine under Case 3, a reduction factor is applied to high density stands above 20,000 trees/ha. This reduction factor reduces the diameter increment by a factor of 1.0 to 0.5 as density increases from 20,000 to 70,000 trees/ha (Equations D22 and D23). The reduction factor is fixed at 0.5 at densities above 70,000 trees/ha.

Where:

*DIhd _{ijk}* = Diameter increment (cm/year) of tree

*k*of species

*j*in stand

*i*under the pine high density reduction factor

*DI _{ijk}* = Diameter increment (cm/year) of tree

*k*of species

*j*in stand

*i*from Equation D21

*RF _{ijk}* = Reduction factor of tree

*k*of species

*j*in stand

*i*;

*RF*may not exceed 0.5 in Equation D22.

_{ijk}*TDen _{i}* = Total density (trees/ha) of stand

*i*

#### Trembling Aspen Group

For trembling aspen > 1.3m tall and < 4cm DBH, annual diameter increment is modeled using Equation D24 for Case 2 and Equations D25 to D28 for Case 3. Equations D25 to D28 implement a basal area increment model.

**Case 2**

Where:

*DBH _{ijk}* = Diameter at breast height (cm) of tree

*k*of species

*j*in stand

*i*

*H1 _{ijk}* = Height (m) of tree

*k*of species

*j*in stand

*i*after the current year’s height growth

*a _{0}* = -1.83,

*a*= 1.444

_{1}

A minimum DBH of 0.3m is applied to any prediction < 0.3cm DBH (Equation D24).

**Case 3**

Where:

*BAI _{ijk}* = Basal area increment (cm²/year) of tree

*k*of species

*j*in stand

*i*

*BA _{ijk}* = Basal area (cm²) of tree

*k*of species

*j*in stand

*i*from the previous year

*H1 _{ijk}* = Height (m) of tree

*k*of species

*j*in stand

*i*after the current year’s height growth

*MHtall10D _{i}* = Mean height (m) of the tallest 10% of deciduous trees (i.e. deciduous canopy height) in stand

*i*

*SFDec* = Deciduous spacing factor from Equation D26; Maximum Value = 2; If *DDen _{i}* = 0 then

*SFDec*= 2.

*DDen _{i}* = Deciduous tree density (trees/ha) of stand

*i*

*BAnew _{ijk}* = New basal area increment (cm²) of tree

*k*of species

*j*in stand

*i*

*DI _{ijk}* = Diameter increment (cm/year) of tree

*k*of species

*j*in stand

*i*

*DBH _{ijk}* = Diameter at breast height (cm) of tree

*k*of species

*j*in stand

*i*

*a _{0}* = 0.7898,

*a*= 0.6653,

_{1}*a*= 0.3452,

_{2}*a*= 0.2844

_{3}

#### Black Spruce Group

**Case 2**

For black spruce under Case 2, the initial diameter above breast height is predicted using Equation D20. This equation uses the black spruce subregional parameters in Huang (1994).

**Case 3**

For black spruce under Case 3, annual diameter increment is modeled using Equation D2 (Oboite and Comeau 2021). This relationship was developed using 4,139 plots, ranging from Alaska to Manitoba. See Equation D2 and Oboite and Comeau (2021) for more details.

^{1} This function resides in the MGM21 Workbook (VBA)

^{2} This function resides in the MGM21 Growth Engine (DLL)

Last Modified: December 13, 2021