# Height Growth of Korean Pine Seedlings Planted under Strip-Cut Larch Plantations in Northeast China

Authors: Toshiaki Owari; Shinichi Tatsumi; Liangzhi Ning; Mingfang Yin
Journal: International Journal of Forestry Research
 (2015)

---

## Abstract

To develop two-storied forest management of larch plantations in Northeast China, this study examined the height growth of Korean pine (Pinus koraiensis Sieb. et Zucc.) seedlings planted under strip-cut larch canopies. We measured the height growth of the underplanted seedlings 4 years after planting. The larch canopies were of varying stand age (12, 17, and 37 years) and strip-cut width (4.5, 6.0, and 7.5 m). We measured the seedling height growth in an open site (i.e., a site with no canopy). Underplanted seedlings had a smaller height growth (12.1–20.1 cm year−1) than the seedlings planted in the open site (23.7 cm year−1). The seedlings underplanted in the wider strip-cuts tended to have greater height growth than the seedlings underplanted in the narrowest strip-cuts. A generalized linear mixed model analysis predicted the greatest seedling height growth in the open site. A 36–47% reduction in annual height growth was predicted for the narrowest strip-cuts (4.5 m) versus the open site, while a 13–36% reduction in annual height growth was predicted for the wider strip-cuts (6.0–7.5 m) versus the open site. To maintain adequate height growth, forest managers are recommended to create wider strip-cuts (i.e., ≥6.0 m) for the purpose of underplanting Korean pine seedlings in larch plantations.


---

## Body

## 1. Introduction

Larches (Larix spp.) are important fast-growing and high-yield tree species in Northeast China, and government policy has keenly promoted larch plantations in the region [1]. According to the National Forest Resource Survey of 1988, there were 1.2 million hectares of larch plantations in the region [2]. However, short-rotation management of pure larch plantations may have some disadvantages, such as soil degradation, loss of biodiversity, high vulnerability to disease, and insect outbreaks [3, 4]. To prevent water shortages and other environmental threats, a large proportion of the planted larch forest has been designated “ecological public welfare forest” where clear-cutting is prohibited [5]. To enhance ecosystem functioning while maintaining the economic benefits of larch plantations, a pressing need exists to develop alternative forest management regimes in Northeast China [5, 6].Forest management whereby even-aged plantation monocultures are transformed to two-storied mixed-species stands has been practiced in many forests around the world [7–9]. Underplanting is combined with the partial cutting of the existing stand to promote the development of complex stand structure with multiple canopy layers and various tree sizes [9]. Two-storied forest management in which overstory larch trees are strip-cut and Korean pine (Pinus koraiensis Sieb. et Zucc.) seedlings are underplanted has been deemed a viable option in Northeast China [6, 10]. Korean pine is an ecologically key species in the region, where mixed conifer-broad-leaved forests dominated by Korean pine are naturally distributed [11, 12]. Korean pine has great economic value, because it produces high-quality timber and commercial seeds for eating and medical use [13, 14].Korean pine is shade-tolerant in the early stages of growth but becomes gradually shade-intolerant as it ages [11, 15]. In Northeast China, Korean pine has frequently been planted in broad-leaved secondary forests to restore a climax-like vegetative cover; a large area of mixed forest has been established using this approach [12, 16]. Previous studies have analyzed the growth characteristics of Korean pine planted under overstory canopies with results indicating that the light environment is an important factor influencing its growth [12, 16–19]. Consequently, additional research has examined the photosynthetic traits of Korean pine seedlings and saplings [6, 11, 16, 20–23]. Limited knowledge exists on the growth patterns of Korean pine seedlings planted under the larch canopy openings created by strip-cutting operations.The purpose of this study was to examine the height growth of Korean pine seedlings planted under the larch canopy openings created by different types of strip-cutting treatments in Northeast China. A case analysis was conducted at two-storied forest stands established in eastern Liaoning Province. We employed a statistical modeling technique to quantify the height growth of underplanted Korean pine seedlings. Overstory larch canopies and the width of strip-cut were used as surrogates of light conditions, and the possible effects on seedling size were analyzed based on the derived growth model. Our analysis aimed to find optimal levels of strip-width for securing adequate height growth. We concluded with a brief discussion of management implications for the practice of underplanting Korean pine in strip-cut larch plantations.


## 2. Materials and Methods

### 2.1. Study Site

We conducted this study at the Dabiangou Forest Farm of Qingyuan County in Liaoning Province (42°00–02′N, 125°03–04′E, 480–520 m above sea level). The forest is located in the mid to lower slopes of the mountainous area of the Changbai Mountain Longgang branch. The region has a marked continental monsoon climate featuring cold, dry winters and hot, rainy summers [24]. The mean annual temperature and mean annual precipitation are 5°C and 700–800 mm, respectively, with a maximum temperature of 36.5°C, a minimum temperature of −37.7°C, and a frost-free period of 120–125 days [6]. The major soil type is brown forest soil, with a soil depth of 30–50 cm [10]. The main plantation species in the region are conifers such asLarix kaempferi (Lamb.) Carr.,Larix olgensis Henry,Pinus koraiensis, andPinus tabulaeformis Carr. [25].We conducted field surveys in five larch plantations that had been established previously by the Dabiangou Forest Farm as experimental stands [6, 26]. The stands were composed of eitherL. olgensis orL. kaempferi. Based on the stand ages and the types of strip-cutting, we hereafter refer to these five plots as A12-W4.5, A12-W6.0, A17-W4.5, A17-W6.0, and A37-W7.5 (Table 1). The larches had originally been planted at a spacing of 1.5 m × 1.5 m in each stand, although some trees had been thinned or were dead at the time of survey.Table 1
Stand parameters and canopy openness for the plots with a larch overstory.
Plot ID
Stand age (year)
Strip-cut width (m)
Plot size (m2)
Tree density∗ (trees ha−1)
Mean DBH∗ (cm)
Mean height∗ (m)
Basal area∗ (m2 ha−1)
Mean canopy openness∗
∗ (%)
2008
2008
2011
2008
2011
2008
2011
2008
2011
2012
A12-W4.5
12
4.5
612
1683
1634
7.9 (1.8)
10.1 (2.5)
9.3 (1.3)
11.1 (1.5)
8.6
14.0
7.7a (0.8)
A12-W6.0
12
6.0
378
1429
1376
9.5 (2.1)
11.7 (2.6)
8.9 (1.1)
10.7 (1.2)
10.7
15.6
8.7a (0.7)
A17-W4.5
17
4.5
234
1282
1282
11.1 (2.3)
13.2 (2.7)
11.9 (1.6)
13.3 (1.6)
12.9
18.2
9.7a (0.7)
A17-W6.0
17
6.0
297
1246
1246
11.3 (1.9)
13.5 (2.3)
11.8 (1.2)
13.3 (1.2)
12.7
18.2
15.4b (0.4)
A37-W7.5
37
7.5
315
381
381
23.8 (3.4)
25.3 (3.5)
22.5 (1.0)
23.0 (1.0)
17.3
19.4
20.0c (2.0)
Figures in parentheses provide the standard deviations.∗Sources: measurement records by the Dabiangou Forest Farm. ∗
∗
n = 5 (A12-W4.5, A12-W6.0); n = 4 (the other plots). Canopy openness values followed by different letters differ significantly according to Tukey’s test for multiple comparisons (P < 0.001).Different types of strip-cutting were conducted in a repetitive manner in spring 2008 [6, 26]. A12-W4.5 and A17-W4.5 were strip-cut by 2-row removal with 2-row retention, A12-W6.0 and A17-W6.0 were strip-cut by 3-row removal with 3-row retention, and A37-W7.5 was strip-cut by 4-row removal with 3-row retention. A12-W4.5, A12-W6.0, A17-W4.5, and A17-W6.0 were located on west/southwest-facing slopes with slope angles of 26–37%, while A37-W7.5 was located on a steep east-facing slope with a slope angle of 49%. These stands were located at the mid part of the slopes, and the strips were orientated in the slope direction (uphill-downhill). The forest floor was mostly covered by shrubs and herbs in A17-W4.5, A17-W6.0, and A37-W7.5, while the soil surface was partly exposed with sparse vegetative cover in A12-W4.5 and A12-W6.0.In 2008, Korean pine seedlings were underplanted in the openings created by the strip-cuts [26]. The seedlings were planted by the Forest Farm according to the strip width: in a line with 2.0 m spacing at A12-W4.5, in a line with 1.2 m spacing at A17-W4.5, in two lines with 2.0 × 1.2 m spacing at A12-W6.0, and 1.2 × 1.2 m spacing at A17-W6.0. Despite the widest strip the seedlings were sparsely planted in a center line with 1.2 m spacing at A37-W7.5, because natural tree regeneration was also expected to occur.We set up one rectangular plot (234–612 m2) in each stand. The plot A12-W4.5 spanned two strip openings while the other plots had one strip, in which Korean pine seedlings had been planted. The retained rows of larches on both sides of the strip(s) were included in the plot area. In each plot, the diameter at breast height (DBH) and height for all larch trees had been measured by the Dabiangou Forest Farm every year. Table 1 summarizes the stand parameters of overstory larch trees within each plot (tree density, mean DBH, mean tree height, and stand basal area).Table1 also shows the mean canopy openness (%), which was measured in June 2012 by taking hemispherical photographs at four or five randomly selected locations above the underplanted Korean pine seedlings in each plot. The photographs were taken at breast height (130 cm) using a digital camera (GX200, Ricoh Company, Ltd., Tokyo, Japan) fitted with a fisheye lens (UWC-0195, FIT Corporation, Nagano, Japan). We used imaging software (Gap Light Analyzer ver. 2.0) to estimate the canopy openness [27]. Significant differences were found among the plots (P
<
0.001; one-way ANOVA). Generally, canopy openness was greater in the plots having older stand ages and wider strips after 4 years of the strip-cutting treatments.In addition to the five plots mentioned above, we set up one 320 m2 control plot in a young plantation hereafter referred to as OPEN to examine seedling growth in open sites. The plot OPEN was located at a distance of <5 km from the other plots and at the mid part of an east-facing slope with a slope angle of 44% where Korean pine and Korean spruce (Picea koraiensis Nakai) had been planted in 2007 alternately at a spacing of 2.0 × 2.0 m. Shrubs and herbs covered most of the forest floor.


### 2.2. Data Collection and Analysis

Field measurement was conducted in mid-June 2012; four growing seasons after the Korean pine seedlings were underplanted in the larch plantations and five growing seasons after the planting in OPEN [26]. The annual height growth of Korean pine seedlings was estimated by measuring the distance between the successive terminal bud scars (internodes or branch whorls) along the main stems downward [28, 29] to a precision of 0.1 cm. We assumed that the current-year shoot growth had completed by early June [15] and thus regarded the length of the new shoots as the height growth in 2012. The number of Korean pine seedlings in a given plot ranged from 9 to 25, with a total of 85. The number of seedlings varied among the plots due to differences in plot size, initial planting density, and mortality. The mean rate of seedling survival was 66% (85/128), with the range between 39 and 79%. The number of annual height growth measurements in a given plot ranged from 36 to 100, with a total of 339. We estimated the mean height of Korean pine seedlings for each plot and each growing season by subtracting the annual height growths backward from the current tree height. Because the homogeneity of variance assumption was questionable (P
<
0.05; Bartlett’s test), we employed a nonparametric Scheffe’s test for multiple comparisons to analyze differences between the plots.Next, a generalized linear mixed model (GLMM) analysis [30] was conducted to examine the annual height growth of Korean pine seedlings under different canopy conditions. GLMM provides a flexible modeling framework that allows the handling of data with a nonnormal distribution and the incorporation of random effects [31]. We assumed that the logarithm of the height growth followed a normal distribution, because growth is a nonnegative continuous variable [32]:(1)
log
⁡
G
i
j
k
~
Norm
log
⁡
G
-
i
j
k
,
σ
2
,where G
i
j
k is the observed annual height growth (cm year−1) of seedling i in plot j at the kth growing season, G
-
i
j
k is the expected annual height growth (cm year−1), and σ
2 is the variance.Using the logarithmic link function [31], we expressed annual seedling growth as a function of three fixed effects, seedling height at the beginning of the growing season, basal area of overstory larch trees, and strip-cut width. We assumed that the effect of strip-cut treatment on height growth would vary according to residual overstory basal area. We included random effects to take into account the pseudoreplication within individual seedlings and plots. The growing years were considered as a temporal random effect:(2)
log
⁡
G
-
i
j
k
=
α
0
+
α
1
H
i
j
k
+
α
2
+
α
3
B
j
k
W
j
+
t
r
e
e
i
+
p
l
o
t
j
+
y
e
a
r
k
,where H
i
j
k is the initial height (cm) of seedling i in plot j at the beginning of the kth growing season; B
j
k is the basal area (m2 ha−1) of overstory larch trees in plot j at the kth growing season; W
j is a categorical variable that represents the strip-cut width in plot j (i.e., 4.5, 6.0, and 7.5 m; W
j
=
0 at OPEN); t
r
e
e
i, p
l
o
t
j, and y
e
a
r
k are the random-effect parameters for individual seedling i, plot j, and the kth growing season; and α
0, α
1, α
2, and α
3 are fixed-effect parameters.Backward stepwise selection method based on Akaike’s information criterion (AIC) [33] was used to select the best model [34]. We evaluated the goodness of fit of the selected model by determining the coefficient of determination R
2 for GLMMs [35]. We further used a leave-one-out cross-validation approach, because neither independent data nor enough data to split into two groups were available [36]. At each iteration, one observation was excluded from the data set and the selected model was fitted to the remaining observations [37]. The analyses were conducted using the “lme4” [38], “lmerTest” [39], and “MuMIn” [40] packages in the “R” 3.0.3 statistical software [41].Using the selected model with the estimated parameters, the annual height growth of Korean pine seedlings was predicted. We further simulated changes in seedling height in the early stages of growth by accumulating the predicted annual growth for each year. In the simulation, the initial seedling height at the time of planting was set at 20.0 cm. Our computations were continued for 7 growing years after planting.


## 2.1. Study Site

We conducted this study at the Dabiangou Forest Farm of Qingyuan County in Liaoning Province (42°00–02′N, 125°03–04′E, 480–520 m above sea level). The forest is located in the mid to lower slopes of the mountainous area of the Changbai Mountain Longgang branch. The region has a marked continental monsoon climate featuring cold, dry winters and hot, rainy summers [24]. The mean annual temperature and mean annual precipitation are 5°C and 700–800 mm, respectively, with a maximum temperature of 36.5°C, a minimum temperature of −37.7°C, and a frost-free period of 120–125 days [6]. The major soil type is brown forest soil, with a soil depth of 30–50 cm [10]. The main plantation species in the region are conifers such asLarix kaempferi (Lamb.) Carr.,Larix olgensis Henry,Pinus koraiensis, andPinus tabulaeformis Carr. [25].We conducted field surveys in five larch plantations that had been established previously by the Dabiangou Forest Farm as experimental stands [6, 26]. The stands were composed of eitherL. olgensis orL. kaempferi. Based on the stand ages and the types of strip-cutting, we hereafter refer to these five plots as A12-W4.5, A12-W6.0, A17-W4.5, A17-W6.0, and A37-W7.5 (Table 1). The larches had originally been planted at a spacing of 1.5 m × 1.5 m in each stand, although some trees had been thinned or were dead at the time of survey.Table 1
Stand parameters and canopy openness for the plots with a larch overstory.
Plot ID
Stand age (year)
Strip-cut width (m)
Plot size (m2)
Tree density∗ (trees ha−1)
Mean DBH∗ (cm)
Mean height∗ (m)
Basal area∗ (m2 ha−1)
Mean canopy openness∗
∗ (%)
2008
2008
2011
2008
2011
2008
2011
2008
2011
2012
A12-W4.5
12
4.5
612
1683
1634
7.9 (1.8)
10.1 (2.5)
9.3 (1.3)
11.1 (1.5)
8.6
14.0
7.7a (0.8)
A12-W6.0
12
6.0
378
1429
1376
9.5 (2.1)
11.7 (2.6)
8.9 (1.1)
10.7 (1.2)
10.7
15.6
8.7a (0.7)
A17-W4.5
17
4.5
234
1282
1282
11.1 (2.3)
13.2 (2.7)
11.9 (1.6)
13.3 (1.6)
12.9
18.2
9.7a (0.7)
A17-W6.0
17
6.0
297
1246
1246
11.3 (1.9)
13.5 (2.3)
11.8 (1.2)
13.3 (1.2)
12.7
18.2
15.4b (0.4)
A37-W7.5
37
7.5
315
381
381
23.8 (3.4)
25.3 (3.5)
22.5 (1.0)
23.0 (1.0)
17.3
19.4
20.0c (2.0)
Figures in parentheses provide the standard deviations.∗Sources: measurement records by the Dabiangou Forest Farm. ∗
∗
n = 5 (A12-W4.5, A12-W6.0); n = 4 (the other plots). Canopy openness values followed by different letters differ significantly according to Tukey’s test for multiple comparisons (P < 0.001).Different types of strip-cutting were conducted in a repetitive manner in spring 2008 [6, 26]. A12-W4.5 and A17-W4.5 were strip-cut by 2-row removal with 2-row retention, A12-W6.0 and A17-W6.0 were strip-cut by 3-row removal with 3-row retention, and A37-W7.5 was strip-cut by 4-row removal with 3-row retention. A12-W4.5, A12-W6.0, A17-W4.5, and A17-W6.0 were located on west/southwest-facing slopes with slope angles of 26–37%, while A37-W7.5 was located on a steep east-facing slope with a slope angle of 49%. These stands were located at the mid part of the slopes, and the strips were orientated in the slope direction (uphill-downhill). The forest floor was mostly covered by shrubs and herbs in A17-W4.5, A17-W6.0, and A37-W7.5, while the soil surface was partly exposed with sparse vegetative cover in A12-W4.5 and A12-W6.0.In 2008, Korean pine seedlings were underplanted in the openings created by the strip-cuts [26]. The seedlings were planted by the Forest Farm according to the strip width: in a line with 2.0 m spacing at A12-W4.5, in a line with 1.2 m spacing at A17-W4.5, in two lines with 2.0 × 1.2 m spacing at A12-W6.0, and 1.2 × 1.2 m spacing at A17-W6.0. Despite the widest strip the seedlings were sparsely planted in a center line with 1.2 m spacing at A37-W7.5, because natural tree regeneration was also expected to occur.We set up one rectangular plot (234–612 m2) in each stand. The plot A12-W4.5 spanned two strip openings while the other plots had one strip, in which Korean pine seedlings had been planted. The retained rows of larches on both sides of the strip(s) were included in the plot area. In each plot, the diameter at breast height (DBH) and height for all larch trees had been measured by the Dabiangou Forest Farm every year. Table 1 summarizes the stand parameters of overstory larch trees within each plot (tree density, mean DBH, mean tree height, and stand basal area).Table1 also shows the mean canopy openness (%), which was measured in June 2012 by taking hemispherical photographs at four or five randomly selected locations above the underplanted Korean pine seedlings in each plot. The photographs were taken at breast height (130 cm) using a digital camera (GX200, Ricoh Company, Ltd., Tokyo, Japan) fitted with a fisheye lens (UWC-0195, FIT Corporation, Nagano, Japan). We used imaging software (Gap Light Analyzer ver. 2.0) to estimate the canopy openness [27]. Significant differences were found among the plots (P
<
0.001; one-way ANOVA). Generally, canopy openness was greater in the plots having older stand ages and wider strips after 4 years of the strip-cutting treatments.In addition to the five plots mentioned above, we set up one 320 m2 control plot in a young plantation hereafter referred to as OPEN to examine seedling growth in open sites. The plot OPEN was located at a distance of <5 km from the other plots and at the mid part of an east-facing slope with a slope angle of 44% where Korean pine and Korean spruce (Picea koraiensis Nakai) had been planted in 2007 alternately at a spacing of 2.0 × 2.0 m. Shrubs and herbs covered most of the forest floor.


## 2.2. Data Collection and Analysis

Field measurement was conducted in mid-June 2012; four growing seasons after the Korean pine seedlings were underplanted in the larch plantations and five growing seasons after the planting in OPEN [26]. The annual height growth of Korean pine seedlings was estimated by measuring the distance between the successive terminal bud scars (internodes or branch whorls) along the main stems downward [28, 29] to a precision of 0.1 cm. We assumed that the current-year shoot growth had completed by early June [15] and thus regarded the length of the new shoots as the height growth in 2012. The number of Korean pine seedlings in a given plot ranged from 9 to 25, with a total of 85. The number of seedlings varied among the plots due to differences in plot size, initial planting density, and mortality. The mean rate of seedling survival was 66% (85/128), with the range between 39 and 79%. The number of annual height growth measurements in a given plot ranged from 36 to 100, with a total of 339. We estimated the mean height of Korean pine seedlings for each plot and each growing season by subtracting the annual height growths backward from the current tree height. Because the homogeneity of variance assumption was questionable (P
<
0.05; Bartlett’s test), we employed a nonparametric Scheffe’s test for multiple comparisons to analyze differences between the plots.Next, a generalized linear mixed model (GLMM) analysis [30] was conducted to examine the annual height growth of Korean pine seedlings under different canopy conditions. GLMM provides a flexible modeling framework that allows the handling of data with a nonnormal distribution and the incorporation of random effects [31]. We assumed that the logarithm of the height growth followed a normal distribution, because growth is a nonnegative continuous variable [32]:(1)
log
⁡
G
i
j
k
~
Norm
log
⁡
G
-
i
j
k
,
σ
2
,where G
i
j
k is the observed annual height growth (cm year−1) of seedling i in plot j at the kth growing season, G
-
i
j
k is the expected annual height growth (cm year−1), and σ
2 is the variance.Using the logarithmic link function [31], we expressed annual seedling growth as a function of three fixed effects, seedling height at the beginning of the growing season, basal area of overstory larch trees, and strip-cut width. We assumed that the effect of strip-cut treatment on height growth would vary according to residual overstory basal area. We included random effects to take into account the pseudoreplication within individual seedlings and plots. The growing years were considered as a temporal random effect:(2)
log
⁡
G
-
i
j
k
=
α
0
+
α
1
H
i
j
k
+
α
2
+
α
3
B
j
k
W
j
+
t
r
e
e
i
+
p
l
o
t
j
+
y
e
a
r
k
,where H
i
j
k is the initial height (cm) of seedling i in plot j at the beginning of the kth growing season; B
j
k is the basal area (m2 ha−1) of overstory larch trees in plot j at the kth growing season; W
j is a categorical variable that represents the strip-cut width in plot j (i.e., 4.5, 6.0, and 7.5 m; W
j
=
0 at OPEN); t
r
e
e
i, p
l
o
t
j, and y
e
a
r
k are the random-effect parameters for individual seedling i, plot j, and the kth growing season; and α
0, α
1, α
2, and α
3 are fixed-effect parameters.Backward stepwise selection method based on Akaike’s information criterion (AIC) [33] was used to select the best model [34]. We evaluated the goodness of fit of the selected model by determining the coefficient of determination R
2 for GLMMs [35]. We further used a leave-one-out cross-validation approach, because neither independent data nor enough data to split into two groups were available [36]. At each iteration, one observation was excluded from the data set and the selected model was fitted to the remaining observations [37]. The analyses were conducted using the “lme4” [38], “lmerTest” [39], and “MuMIn” [40] packages in the “R” 3.0.3 statistical software [41].Using the selected model with the estimated parameters, the annual height growth of Korean pine seedlings was predicted. We further simulated changes in seedling height in the early stages of growth by accumulating the predicted annual growth for each year. In the simulation, the initial seedling height at the time of planting was set at 20.0 cm. Our computations were continued for 7 growing years after planting.


## 3. Results

### 3.1. Observed Height Growth of Korean Pine Seedlings

Figure1 shows the changes in the observed mean height of Korean pine seedlings over the 4 years after planting. The initial height of individual seedlings at the time of planting (2007-2008) was 18.2–24.8 cm, and there were no significant differences in mean height among the plots (P
=
0.13–1.00; Scheffe’s test). The difference gradually got larger over time and became significant at 4 years after planting. A12-W4.5 (69.8 cm), A17-W6.0 (76.9 cm), and A17-W4.5 (79.6 cm) had significantly lower mean seedling heights than OPEN (114.2 cm; P
<
0.05), while A12-W6.0 (92.1 cm; P
=
0.38) and A37-W7.5 (101.0 cm; P
=
0.85) did not have significant differences in seedling height from OPEN.Figure 1
The observed mean height of Korean pine seedlings for 4 years after planting.Significant positive correlations were found between the initial seedling height at the beginning of each growing season and the observed annual height growth at the growing season for each plot (P
<
0.01; Figure 2), with Pearson correlation coefficients of 0.40–0.77. The mean annual height growth was the largest in OPEN (23.7 cm year−1), followed by A37-W7.5 (20.1 cm year−1). All other plots (i.e., other than A37-W7.5) had significantly lower height growth (12.1–17.3 cm year−1) than OPEN (P
<
0.01; Scheffe’s test).Figure 2
Relationships between the initial height at the beginning of the growing season and the annual height growth of Korean pine seedlings (n
=
339).


### 3.2. Model Selection, Parameter Estimation, and Goodness of Fit

The GLMM inclusive of all three fixed effects (i.e., initial seedling height, stand basal area of overstory larch trees, and strip-cut width) was selected as the best model (Table2), with the lowest AIC value of 167.4. The second smallest AIC value was 172.6, when basal area of overstory larch trees was excluded as a fixed effect. According to the parameter estimates, initial seedling height had a positive effect on the annual height growth. Additionally, the height growth of underplanted seedlings in the narrowest strip-cut width (i.e., 4.5 m) was significantly smaller than that in OPEN.Table 2
Parameter estimates from the analysis of the annual height growth of Korean pine seedlings using a generalized linear mixed model.
Variable (parameter)
Level
Estimate
Standard error
t value
P value
Intercept (α
0)
2.6832
0.0807
33.24
<0.001
Initial seedling height (α
1)
0.0089
0.0012
7.52
<0.001
Strip-cut width (α
2)
4.5 m
−0.8541
0.1478
−5.78
<0.001
6.0 m
0.0650
0.2230
0.29
0.771
7.5 m
0.6021
1.1893
0.51
0.613
Basal area of larch trees × Strip-cut width (α
3)
4.5 m
0.0227
0.0104
2.18
0.031
6.0 m
−0.0284
0.0153
−1.86
0.064
7.5 m
−0.0435
0.0647
−0.67
0.502
Note: AIC = 167.4.The selected model and its associated parameter estimates provided an acceptable fit to the data, with a marginalR
G
L
M
M
(
m
)
2 (i.e., variance explained by fixed factors; [35]) of 0.49 and a conditional R
G
L
M
M
(
c
)
2 (i.e., variance explained by fixed and random factors) of 0.55. The leave-one-out cross-validation predictions fitted adequately to the observed height growth of seedlings (R
2
=
0.5
1; Figure 3), having a root mean square error of 4.7 cm.Figure 3
Relationships between observed and predicted annual height growth of Korean pine seedlings (n
=
339). Leave-one-out cross-validation was used for the prediction. The dotted line represents the point of identical observed/predicted height growth.


### 3.3. Predicted Height Growth of Korean Pine Seedlings

Figure4 shows the predicted annual height growth of Korean pine seedlings as a function of initial seedling height and stand basal area of overstory larch trees in each strip-cut width using the selected model (random-effect parameters were not used for the prediction). Seedlings planted in OPEN showed the largest predicted growth. The predicted growth of 120-cm seedlings (42.4 cm year−1) was more than twice the predicted growth of 20-cm seedlings (17.4 cm year−1). Generally, underplanted seedlings in wider strip-cuts tended to have relatively large height growth than those in narrower strip-cuts. That is, seedlings in the widest strip-cut (i.e., 7.5 m) had the largest predicted growth, whereas height growth was strongly suppressed at the narrowest strip-cut (i.e., 4.5 m). The annual height growths in the 4.5 m strip were 36−47% smaller than in OPEN. In the wider strip-cuts (6.0–7.5 m), 13–36% reductions in annual height growth were predicted versus OPEN.Figure 4
Relationships between the initial height and annual height growth of Korean pine seedlings, predicted by the generalized linear mixed model. In each strip-cut width, the basal area of overstory larch trees (BA) within a range of the plot data (Table1) was used for the prediction.
(a)
Strip width: 4.5 m
(b)
Strip width: 6.0 m
(c)
Strip width: 7.5 mFigure5 shows the predicted change in seedling height over a period of 7 years after planting in each strip-cut width using the selected model. The seedlings in OPEN were predicted to reach 301.6 cm after 7 years. The predicted seedling heights in the underplanted sites were smaller than those in OPEN. Seedling heights after 7 years in the 4.5 m strip were 109.4–112.9 cm, which corresponded to 36-37% of the predicted seedling height in OPEN. In the wider strip-cuts (6.0–7.5 m), seedlings had greater predicted heights (136.8–222.9 cm) after 7 years compared to the narrowest strip-cut.Figure 5
Predicted height of Korean pine seedlings for 7 growing years after plantation. The initial seedling height when planted was set at 20 cm for all plots, and the area of overstory larch trees (BA) was assumed to be constant over those 7 years.
(a)
Strip width: 4.5 m
(b)
Strip width: 6.0 m
(c)
Strip width: 7.5 m


## 3.1. Observed Height Growth of Korean Pine Seedlings

Figure1 shows the changes in the observed mean height of Korean pine seedlings over the 4 years after planting. The initial height of individual seedlings at the time of planting (2007-2008) was 18.2–24.8 cm, and there were no significant differences in mean height among the plots (P
=
0.13–1.00; Scheffe’s test). The difference gradually got larger over time and became significant at 4 years after planting. A12-W4.5 (69.8 cm), A17-W6.0 (76.9 cm), and A17-W4.5 (79.6 cm) had significantly lower mean seedling heights than OPEN (114.2 cm; P
<
0.05), while A12-W6.0 (92.1 cm; P
=
0.38) and A37-W7.5 (101.0 cm; P
=
0.85) did not have significant differences in seedling height from OPEN.Figure 1
The observed mean height of Korean pine seedlings for 4 years after planting.Significant positive correlations were found between the initial seedling height at the beginning of each growing season and the observed annual height growth at the growing season for each plot (P
<
0.01; Figure 2), with Pearson correlation coefficients of 0.40–0.77. The mean annual height growth was the largest in OPEN (23.7 cm year−1), followed by A37-W7.5 (20.1 cm year−1). All other plots (i.e., other than A37-W7.5) had significantly lower height growth (12.1–17.3 cm year−1) than OPEN (P
<
0.01; Scheffe’s test).Figure 2
Relationships between the initial height at the beginning of the growing season and the annual height growth of Korean pine seedlings (n
=
339).


## 3.2. Model Selection, Parameter Estimation, and Goodness of Fit

The GLMM inclusive of all three fixed effects (i.e., initial seedling height, stand basal area of overstory larch trees, and strip-cut width) was selected as the best model (Table2), with the lowest AIC value of 167.4. The second smallest AIC value was 172.6, when basal area of overstory larch trees was excluded as a fixed effect. According to the parameter estimates, initial seedling height had a positive effect on the annual height growth. Additionally, the height growth of underplanted seedlings in the narrowest strip-cut width (i.e., 4.5 m) was significantly smaller than that in OPEN.Table 2
Parameter estimates from the analysis of the annual height growth of Korean pine seedlings using a generalized linear mixed model.
Variable (parameter)
Level
Estimate
Standard error
t value
P value
Intercept (α
0)
2.6832
0.0807
33.24
<0.001
Initial seedling height (α
1)
0.0089
0.0012
7.52
<0.001
Strip-cut width (α
2)
4.5 m
−0.8541
0.1478
−5.78
<0.001
6.0 m
0.0650
0.2230
0.29
0.771
7.5 m
0.6021
1.1893
0.51
0.613
Basal area of larch trees × Strip-cut width (α
3)
4.5 m
0.0227
0.0104
2.18
0.031
6.0 m
−0.0284
0.0153
−1.86
0.064
7.5 m
−0.0435
0.0647
−0.67
0.502
Note: AIC = 167.4.The selected model and its associated parameter estimates provided an acceptable fit to the data, with a marginalR
G
L
M
M
(
m
)
2 (i.e., variance explained by fixed factors; [35]) of 0.49 and a conditional R
G
L
M
M
(
c
)
2 (i.e., variance explained by fixed and random factors) of 0.55. The leave-one-out cross-validation predictions fitted adequately to the observed height growth of seedlings (R
2
=
0.5
1; Figure 3), having a root mean square error of 4.7 cm.Figure 3
Relationships between observed and predicted annual height growth of Korean pine seedlings (n
=
339). Leave-one-out cross-validation was used for the prediction. The dotted line represents the point of identical observed/predicted height growth.


## 3.3. Predicted Height Growth of Korean Pine Seedlings

Figure4 shows the predicted annual height growth of Korean pine seedlings as a function of initial seedling height and stand basal area of overstory larch trees in each strip-cut width using the selected model (random-effect parameters were not used for the prediction). Seedlings planted in OPEN showed the largest predicted growth. The predicted growth of 120-cm seedlings (42.4 cm year−1) was more than twice the predicted growth of 20-cm seedlings (17.4 cm year−1). Generally, underplanted seedlings in wider strip-cuts tended to have relatively large height growth than those in narrower strip-cuts. That is, seedlings in the widest strip-cut (i.e., 7.5 m) had the largest predicted growth, whereas height growth was strongly suppressed at the narrowest strip-cut (i.e., 4.5 m). The annual height growths in the 4.5 m strip were 36−47% smaller than in OPEN. In the wider strip-cuts (6.0–7.5 m), 13–36% reductions in annual height growth were predicted versus OPEN.Figure 4
Relationships between the initial height and annual height growth of Korean pine seedlings, predicted by the generalized linear mixed model. In each strip-cut width, the basal area of overstory larch trees (BA) within a range of the plot data (Table1) was used for the prediction.
(a)
Strip width: 4.5 m
(b)
Strip width: 6.0 m
(c)
Strip width: 7.5 mFigure5 shows the predicted change in seedling height over a period of 7 years after planting in each strip-cut width using the selected model. The seedlings in OPEN were predicted to reach 301.6 cm after 7 years. The predicted seedling heights in the underplanted sites were smaller than those in OPEN. Seedling heights after 7 years in the 4.5 m strip were 109.4–112.9 cm, which corresponded to 36-37% of the predicted seedling height in OPEN. In the wider strip-cuts (6.0–7.5 m), seedlings had greater predicted heights (136.8–222.9 cm) after 7 years compared to the narrowest strip-cut.Figure 5
Predicted height of Korean pine seedlings for 7 growing years after plantation. The initial seedling height when planted was set at 20 cm for all plots, and the area of overstory larch trees (BA) was assumed to be constant over those 7 years.
(a)
Strip width: 4.5 m
(b)
Strip width: 6.0 m
(c)
Strip width: 7.5 m


## 4. Discussion

Our GLMM analysis successfully derived a statistical model that allows the assessment of annual height growth of Korean pine seedlings at the early stages of growing. The selected model and estimated parameters fit the observed height growth with an adequate predictive accuracy. The exponential growth pattern delineated in Figure5 is consistent with our field observations (Figure 1), as well as a previous study by Liu and Wang [42], which found an exponential relationship between the height and age of Korean pine seedlings for 5 growing years after planting in the same region.According to our results, the initial seedling height at the beginning of the growing year positively influenced the annual height growth of Korean pine seedlings (Table2). Although dealing with different coniferous tree species, positive correlations have also been observed at plantations of Japanese sugi (Cryptomeria japonica D. Don) and hinoki (Chamaecyparis obtusa (Sieb. et Zucc.) Endl.) when they were small in tree size [43]. Kiyono et al. [44] examined the early height growth of underplanted Japanese hinoki using a multiple regression analysis and achieved a better predictive accuracy of the model by including initial seedling height as an explanatory variable. In the United States, Moores et al. [45] found that initial seedling height had a significant effect on predicting the annual height growth of three shade-tolerant coniferous species, including balsam fir (Abies balsamea (L.) Mill), red spruce (Picea rubens Sarg.), and eastern hemlock (Tsuga canadensis (L.) Carr.). Similarly, in the case of Korean pine, the initial height of seedlings should be considered an important predictor for height growth modeling.Overstory larch canopies significantly affected the seedling height growth of underplanted Korean pine (Table2). As presented in Figure 4, seedlings planted in the open site (OPEN) showed superior height growth versus the underplanted seedlings. A relatively low light regime created by partial canopy closure may inhibit the height growth of underplanted Korean pine seedlings [12, 19]. In eastern Liaoning, Liu and Wang [42] reported a seedling height of 82.2 cm at 4 years after planting in an open site, which is relatively small compared to our control site (114.2 cm at OPEN). This is probably due to their investigations being conducted at a higher elevation site (684 m above sea level) where relatively small-sized seedlings were initially planted (13.3 cm on average).Our results also indicated that the type of strip-cut has a significant effect on the height growth of underplanted Korean pine seedlings. We observed a general increasing trend of canopy openness as the strip-cut width got wider (Table1). In the stand with the narrowest strip (i.e., 4.5 m), the annual height growth was predicted to be approximately one-half and the seedling height after 7 growing years was predicted to be one-third of those in the open site (Figures 4 and 5). Such a cramped strip-cut could not greatly improve the understory light conditions [16]. Although stand age affected canopy openness on equal strip-cut width (Table 1, probably through a change in stand density; [46]), its effect on seedling height growth was ambiguous. Despite its larger canopy openness, A17-W6.0 produced smaller height growth than A12-W6.0 (Figure 1). Available light for tree regeneration establishment (natural or artificial) can be increased through strip-cutting, but not so much as to promote competition by understory shrubs or herbaceous plants [9]. In addition to larch canopy conditions, the competitive effect of understory vegetation may be a factor for consideration in growth predictions.We assume that the seedling height growth observed in the 4.5 m strip-cut is inadequate for forest management purposes. Because larger pine trees generally produce more cone crops [47], smaller height growth may cause economic loss by delaying cone production and the harvest of nuts. To secure preferred height growth, forest managers are recommended to consider wider strip-cuts for two-storied management of larch and Korean pine plantations. Although seedlings in the wider strip-cuts (6.0–7.5 m) had smaller predicted growths than in OPEN (Figure 4), the differences were not statistically significant (Table 2). By creating 6 m or wider strip-cuts, underplanted seedlings were predicted to reach breast height within 1-2 years of seedlings planted at open sites (Figure 5). Strip-cuts that are too wide, however, may negatively affect the growth and photosynthetic capacities of Korean pine as it increases in size [16]. The planting of Korean pine seedlings in wide openings may entail risks including understory competition [9], physiological drought [48], and early stem divergence [17]. Even though our model predicted greater height growth for Korean pine seedlings in wider strip-cuts, forestry practitioners should apply this research carefully.This study was intended to quantitatively illustrate the height growth of underplanted Korean pine seedlings in strip-cut larch plantations in Northeast China. Our results, however, are limited to predicting height growth at the early growing stages; longitudinal studies within a site are needed for a more complete understanding of growth dynamics under partial larch canopies. Data of understory vegetation would be useful to examine the competitive effect on underplanted seedlings. Furthermore, site-specific conditions may have affected the results of this study. Further research should examine growth patterns under various other site conditions. Future studies should provide for a more fully developed understanding of Korean pine growth for the sustainable management of larch plantations in Northeast China.



---

*Source: 178681-2015-07-01.xml*