잣나무 인공림의 1차 및 2차 간벌에 따른 입목생장 특성 분석

본 연구대상지는 강원대학교 학술림 내 간벌강도에 따라 시업처리 후 모니터링 되고 있는 잣나무 인공림이며, 행정구역상 강원도 춘천시 동산면 봉명리와 원창리에 위치하고 있다(Figure 1). 연구대상지인 강원대학교 학술림의 전체 산림면적 3,140 ha 중 750 ha가 인공림으로 구성되어 있는데, 그 중 521 ha(69%)가 잣나무 조림지로써 인공림 면적의 대부분을 차지하고 있다(Kangwon National University, 2020). 토양은 대부분 양토 및 사질양토로 이루어져 있고, 낙엽분해산물과 광물질이 혼합되어 유기질이 풍부하여 식물생육에 적합하다(Choi, 2003). 기상청에서 발표한 춘천지역 과거 1991~2020 30년간 연별평년값 통계에 따르면, 평균기온 11.4°C, 최고 17.4°C, 최저 6.3°C, 강수량 1,341.5 mm이었다(Korea Meteorological Administration, 2021).

Figure 1. Map of the study area in research forest of Kangwon National University, Chuncheon-si, Gangwon-do, South Korea. The plot (a) presents the location of research forest, and plot (b) displays the location of thinning experimental sites for Korean white pine plantation in the research forest.

Download Original Figure

본 연구에 사용된 잣나무 간벌시험림은 총 3지역이며, 각 시험림은 강도간벌, 중도간벌, 약도간벌, 무간벌 설계에 따라 4개의 시험구로 구성되어 있다(Table 1). 제1, 2간벌시험림은 1996년에 20 m × 20 m 정방형 크기의 표준지가 임도를 사이에 두고 설치되었으며(Choi et al., 1996), 제3간벌시험림은 1999년 30 m × 30 m 정방형 크기로 표준지간 완충지역을 고려하여 4개의 시험구가 설치되었다(Choi et al., 1999). 그 중 제1, 2간벌시험림은 1996년도 19년생 때 표준지가 설치되면서 1차 간벌이, 이후 2014년 37년생 때 2차 간벌이 실시되었다(Choi et al., 2014b; 2019). 제3간벌시험림은 1999년도 21년생 때 1차 간벌이, 이후 2016년 38년생 때 2차 간벌이 실시되었다(Choi et al., 2021).

각 시험림의 1차 간벌과 2차 간벌 시업은 우량대경재로 생산될 미래목이 존치할 수 있도록 입목 간 거리, 수형급 등의 정성 및 정량을 고려한 하층간벌이 실시되었다(Choi et al., 1996; 1999; 2019; 2021). 특히 2차 간벌시업은 1차 간벌에 따라 임분밀도가 밀한 시험구에 강도로 벌채를 실시하도록 설계되었다(Table 2). 2차 간벌에서 제1, 2간벌시험림은 무간벌구가 유지되었으며, 제3간벌시험림은 무간벌구를 포함하여 간벌시업이 진행되었다. 해당 시업결과를 간벌목의 평균 흉고직경/간벌 전 전체 입목의 평균 흉고직경으로 나눈 비(Thinning quotient)로 살펴본 결과, 1차 간벌 시 0.66~0.88, 2차 간벌 시 0.70~0.87의 범위로 산출되었다(Table 2). 이는 1.00보다 낮은 수치로써 적절한 하층간벌이 시행된 것으로 분석되었다.

간벌강도에 따른 입목의 생장변화를 분석하기 위한 현장의 매목조사는 시험림에 따라 4~5회 실시하였다. 제1, 2간벌시험림은 1차 간벌이 실시된 1996년에 1차 매목조사를 시작으로 2001년, 2008년, 2014년(2차 간벌), 2019년 총 5회의 매목조사를 실시하였고, 제3간벌시험림은 1차 간벌이 실시된 1999년 1차 매목조사를 시작으로 2010년, 2016년(2차 간벌), 2021년 총 4회의 매목조사를 실시하였다(Choi et al., 1996; 1999; 2001; 2010; 2014b; 2016; 2019; 2021). 현재 각 간벌시험림은 경계부 말뚝과 입목라벨링을 이용하여 보존되고 있으며, 간벌강도에 따른 생장연구 및 차기 현장조사를 위하여 지속적으로 관리되고 있다(Figure 2).

Figure 2. Photos of thinning experimental site 1 in Korean white pine plantation as an example.

Download Original Figure

각 표준지의 측정인자로써 수종, 흉고직경, 수고, 지하고, 수관폭, 수간형질 및 수관급, 입목위치, 고사여부와 같은 다양한 입목인자와 경사 및 사면방향을 포함한 지황인자를 수집하였다. 그 중 본 연구에 이용된 입목측정인자로써 흉고직경은 D-tape를 이용하여 경사 위쪽에서 입목 지상 1.2 m 높이의 직경을 측정하였고, 수고는 Haglof사 Vertax-Ⅳ를 통해 입목의 초두부까지의 높이를 측정하였다(Choi et al., 2009; 2010; 2014; 2015). 경사 및 사면방향은 Suunto사 클리노미터(Clinometer)를 이용하여 경사지 방향에 따라 측정하였다.

본 연구는 간벌강도별 생장모니터링 비교·분석을 위해 수행됨에 따라, 시험구간 완충구역이 없어 간벌효과가 중첩되는 입목을 제외시키고자 하였다. 이를 위해 제1, 2간벌시험림 표준지의 경우 당초 시험구의 크기, 입목생육상태, 존치된 잔존입목수를 고려하여 시험구마다 상·하·좌·우 2 m씩 완충구역이 설정되었다. 따라서 각 시험구 경계부에 위치한 입목들이 제외되었으며, 각 표준지의 16 m × 16 m 내 포함되는 입목을 대상으로 최종분석을 진행하 였다(Jung, 2022). 최종적으로 분석에 사용된 2차 간벌 후 입목 수는 제1시험림 11~49본, 제2시험림 18~36본이었다. 제3시험림은 설치 당시 표준지간 완충구역이 7~40 m로 존재함에 따라, 분석 시 추가적인 완충구역을 설정하지 않았다. 따라서 30 m × 30 m의 표준지 내 전체 본수에 대하여 분석이 진행되었으며, 2차 간벌 후 사용된 입목 수는 50~74본이었다.

본 연구를 위해 현장에서 수집된 측정인자를 기반으로 입목 및 임분 단위의 추가변수들이 산출되었다. 흉고단면적(BA)은 아래 식 1과 같이 산출하였으며, 개체목의 입목재적(V)은 Lee et al.(2017b)의 연구에서 보고된 이변수 재적식 중 적합성이 가장 높은 모델의 파라미터를 적용하여 식 2와 같이 계산하였다. 산출된 변수 중 단위면적당 흉고단면적 및 재적과 같은 임분제원을 산출하기 위하여 각 표준지 내 결과 값은 ha당 수치로 환산되었다.

BA: Tree basal area (m²), D: Tree diameter at breast height (cm),

π: mathematical constant approximately equal to 3.14159,

V: Tree volume (m³), H: Tree height (m).

간벌시험림별 측정시기마다 조사간격은 다소간 상이하였지만, 정기평균생장량 산출법을 이용하여 측정시기 사이의 생장량을 아래 식 3과 같이 적용·산출하였다(Burkhart et al., 2019). 이에 따른 정기평균생장량 값을 활용하여 간벌강도 및 생장기간별 효과를 비교·분석하였다.

PAI = Periodic annual increment of DBH (cm/year), tree volume (m³/year), or stand volume (m³·ha⁻¹·year⁻¹)

Y₂ = DBH, tree volume, or stand volume at end of growth period

Y₁ = DBH, tree volume, or stand volume at start of growth period

T = number of years in growth period

본 연구에서는 간벌강도에 따른 흉고직경 및 재적의 주요 생장변화를 대상으로 일원배치 분산분석(One-way ANOVA) 및 Duncan의 다중비교검정(alpha=0.05) 방법을 사용하여 시험구간 유의적인 생장차이를 분석하였다. 산점도의 생장패턴 분석 시 이해를 돕기 위해 적은 표본 수이지만 단순선형회귀분석을 사용하였으며, 본 연구의 모든 통계적 검정은 R통계소프트웨어를 이용하여 진행되었다(R Core Team, 2019).

본 시험림들에 적용된 간벌률은 시험구별 입목본수 기준 1차 간벌 18%~68%, 2차 간벌 11%~67% 범위였으며, 흉고단면적 기준 간벌률도 시업강도에 따라 다양하였다(Table 2). 국내외 간벌 연구가 입목본수 기준뿐만 아니라 흉고단면적 및 재적 기준 간벌률을 시행하는 것으로 보고되고 있어 향후 다양한 기준 간벌률에 대한 연구가 필요하다고 판단된다(Mäkinen and Isomäki, 2004a; 2004b; Lee and Choi, 2019; 2020). 본 연구의 최종 조사시기 잔존 입목본수는 잣나무림 제1시험지 657~2,035 본/ha, 제2시험지 758~1,347 본/ha, 제3시험지 618~915 본/ha으로 간벌강도에 따라 다양하였다(Table 3). 표준지 설치 후 생장변화와 임분제원현황도 간벌강도 및 조사시기에 따라 상이하였다.

1차 간벌 직후부터 2차 간벌시기까지의 평균 흉고직경 변화를 비교하면, 간벌구의 생장량은 무간벌구에 대비하여 높게 나타났다. 또한, 2차 간벌 직후부터 최종 측정시기 사이의 평균 흉고직경도 강도간벌구에서 더욱 높았다. 한편, 흉고직경과 달리 평균 수고는 ANOVA분석에서 간벌강도에 따른 유의적인 차이 및 경향이 나타나지 않았다. 해당 결과는 선행연구의 결과와 일치하였으며, 수고생장은 직경생장보다 임분밀도에 독립적이기 때문이라고 사료된다(Mäkinen and Isomäki, 2004a; 2004b). 개체목 재적은 간벌강도에 따라 증가하는 경향이었으나, 임분 흉고단면적과 재적은 무간벌구와 같이 임분밀도가 높은 표준지에서 가장 높았다. 국립산림과학원에서 관리하는 임분밀도 및 임령이 유사한 잣나무 간벌시험림과 비교하면, 본 연구대상지의 전반적인 생장현황은 유사하거나 다소 우세하였다(National Institute of Forest Science, 2021). 다만, 본 연구대상지가 잣나무 생육에 비교적 유리한 조건임을 참고하여(Choi, 2003), 후속연구에서는 다양한 지역으로부터의 간벌효과를 분석하여 토양, 기후, 임지생산성과 같은 특성이 고려된 시업지침을 마련해야 할 것으로 사료된다(Hynynen et al., 2002; Mäkinen and Isomäki, 2004b; Lee and Choi, 2014).

시험림별 경급분포의 시계열적 변화를 분석한 결과, 간벌강도가 클수록 평균경급이 커지는 경향이 관측되었다(Figure 3). 제1, 2간벌시험림의 경우, 19년생 때에는 간벌구간 경급분포의 차이가 없었으나, 이후 24년생, 31년생, 37년생까지 임령이 증가함에 따라 간벌강도별 경급분포의 차이가 존재하였다. 2차 간벌 시기에도 Plot D가 무간벌구로 지속됨에 따른 열세목 및 피압목의 경급은 명료히 낮게 분포하였다. 제3간벌시험림의 경우도, 21년생 때에는 간벌구간 경급분포의 차이가 없었으나, 이후 간벌강도에 따른 차이가 32년생과 38년생 때 관측되었다. 이후 잔존 임분밀도에 따라 2차 간벌이 실시되어, 43년생 때 간벌구간 경급분포의 차이는 다소 완화되었다.

Figure 3. DBH distribution by thinning intensity over measurement instances in each experimental site of Korean white pine plantation. DBH class was 2 cm of bandwidth. Each experimental site and plot was thinned as shown in the legend from heavy thinning treatment to light thinning treatment based on first and second thinning intensity. Specific thinning treatments were performed according to Table 2.

Download Original Figure

간벌강도가 클수록 시간이 경과함에 따라 경급이 큰 입목들이 더 많이 관측되었다(Figure 3). 따라서 간벌을 통해 임분밀도에 따른 입목 간 경쟁이 완화되면 임령이 증가하면서 잔존입목들의 평균경급은 더 큰 것으로 분석되었다. 이와 같이 흉고직경 생장촉진과 대경재 생산에 있어 간벌이 효과적인 것은 선행연구의 결과와 일치하였다(Hynynen et al., 2002; Mäkinen and Isomäki, 2004a; 2004b; Bergh et al., 2014; Rytter and Stener, 2014; Niemisto et al., 2018; Fahlvik et al., 2021). 따라서 적정 간벌강도 및 시기는 우량대경재 생산을 위한 주요 인자로써 적정 산림시업 시 더욱 경제적이고 효율적인 산림관리가 가능할 것이라 판단된다.

전체 생장기간에 대한 평균 흉고직경의 변화를 파악한 결과, 조사시기마다 간벌강도별 생장차이가 확인되었다. 간벌구별 흉고직경 범위는 제1시험림에서 19년생 때 12.2~13.9 cm에서 42년생 때 20.8~29.8 cm로, 제2시험림에서 19년생 때 12.8~13.3 cm에서 42년생 때 21.9~25.8 cm로, 제3시험림에서 21년생 때 15.8~19.2 cm에서 43년생 때 25.9~30.7 cm로 변화하였다(Figure 4). 지속가능한 산림자원관리 매뉴얼의 솎아베기 후 35년생, 400 본/ha, 목표 흉고직경 21 cm와 비교하였을 때, 본 연구의 흉고직경 생장이 다소 우세한 것으로 파악되었다(Korea Forest Research Institute, 2005).

Figure 4. Growth change of average DBH by thinning intensity in each experimental site of Korean white pine plantation. Thinning represents the thinning intensity from first thinning to second thinning. The thinning symbols are the acronyms for thinning intensity; H is heavy thinning, M is moderate thinning, L is light thinning, U is un-thinning. Specific thinning treatments were performed according to Table 2. Significant differences between thinning treatments at the final measurement are indicated by different letters based on Duncan’s multiple range test (alpha=0.05).

Download Original Figure

전반적으로 시간이 경과함에 따라 간벌강도간 평균 흉고직경 차이는 증가하였으며, 특히 무간벌구와 강도간벌구의 흉고직경 차이가 가장 크게 나타났다. 예를 들어 제1시험림의 경우 두 간벌구간 평균 흉고직경의 차이는 19년생 때(1차 간벌 직후) 1.7 cm에서 37년생 때(2차 간벌 직전) 4.8 cm로, 37년생 때(2차 간벌 직후) 7.3 cm에서 42년생 때 9.0 cm로 증가하였다. 이러한 생장차이에 대해 ANOVA분석을 실시한 결과, 제1, 2 간벌시험림의 경우 1차 간벌 후 12년이 경과된 31년생부터, 제3 간벌시험림에서는 1차 간벌 후 11년이 경과된 32년생부터 간벌구간 평균 흉고직경의 차이가 유의적으로 나타났다(p<0.05).

흉고직경 정기평균생장량 차이는 간벌강도, 간벌 후 잔존본수, 조사시기 경과에 따라 뚜렷이 구분되었으며, 간벌 직후 초기 생장기간에는 강도간벌구의 생장량이 무간벌구에 비해 높았다(Figure 5). 첫 생장기간 동안 흉고직경 정기평균생장량의 간벌구별 평균 범위는 제1간벌시험림에서 0.60~0.91 cm/year, 제2간벌시험림에서 0.48~0.74 cm/year, 제3간벌시험림에서 0.31~0.66 cm/year이었다. 이후 시간 경과에 따라 전체적으로 생장량이 감소하였고, 간벌구간 생장량의 차이도 줄어드는 경향이 나타났다(Figure 5). 특히, 2차 간벌 실시 직전, 간벌구간 생장량 차이는 가장 작았다. 해당 시기 간벌구간 생장량 범위는 제1시험림 31-37년생 때 0.22~0.30 cm/year, 제2시험림 31-37년생 때 0.22~0.27 cm/year, 제3시험림 32-38년생 때 0.28~0.43 cm/year로 시간경과에 따라 간벌효과가 많이 감소하는 경향이었다.

Figure 5. Box and whiskers plot for the periodic annual increment of DBH by thinning intensity in each experimental site of Korean white pine plantation. Thinning represents the thinning intensity from first thinning to second thinning. The thinning symbols are the acronym for thinning intensity; H is heavy thinning, M is moderate thinning, L is light thinning, U is un-thinning. Specific thinning treatments were performed according to Table 2. Significant differences between thinning treatments at each measurement are indicated by different letters based on Duncan’s multiple range test (alpha=0.05).

Download Original Figure

이후 2차 간벌이 실시되어 흉고직경 정기평균생장량은 다시 높아졌으며, 해당 시기의 생장량은 제1시험림 37-42년 때 0.29~0.67 cm/year, 제2시험림 37-42년생 때 0.30~0.62 cm/year, 제3시험림 38-43년생 때 0.35~0.44 cm/year로 각 시험림의 표준지간 생장량 차이는 통계적으로 유의하였다(p<0.05). 이와 같이 간벌강도별 흉고직경 생장량의 유의적인 차이는 선행연구의 결과와 일치하였다(Hynynen et al., 2002; Makinen, 2004a; 2004b; Bae et al., 2010; Choi et al., 2014a; Bergh et al., 2014; Rytter and Stener, 2014; Niemisto et al., 2018; Fahlvik et al., 2021).

흉고직경 정기평균생장량의 산점도를 도시한 결과, 제1, 2간벌시험림의 간벌구간 생장량 범위는 1차 간벌 직후인 19-24년생 때 0.48~0.95 cm/year, 2차 간벌 직전의 31-37년생 때 0.14~0.29 cm/year로 임령이 증가함에 따라 생장량은 감소하였다(Figure 6). 2차 간벌 이후 37-42년생 때 생장량은 0.33~0.67 cm/year로 다시 증가하였다. 제3간벌시험림의 경우, 1차 간벌 후 21-32년생 때 0.37~0.65 cm/year, 32-38년생 때 0.38~0.51 cm/year로 줄어들었다. 2차 간벌 이후 38-43년생 때 생장량은 0.34~0.48 cm/year로 이전 32-38년생 기간의 생쟝량 수준을 유지하였다.

Figure 6. Scatterplots for the periodic annual increment of DBH by stand density at all measurement instances in each experimental site of Korean white pine plantation. The simple regression lines were used to easily illustrate growth trend using the formula of Y=a×ln(X)+b notwithstanding the insufficient number of observations.

Download Original Figure

단순회귀분석의 추세선을 이용하여 간벌효과에 따른 해당 경향을 살펴보면, 각 측정시기별로 임분밀도가 낮을수록 흉고직경 정기평균생장량은 증가하였다(Figure 6). 또한, 임분밀도에 따라 1차 간벌 직후부터 2차 간벌 직전 시기까지 흉고직경 정기평균생장량은 지속적으로 감소하였으나, 2차 간벌 이후 최종 측정시기까지의 생장량은 유지되거나 다시 증가하는 경향이 관측되었다.

간벌구별 평균 개체목 재적의 범위는 제1, 2간벌시험림에서 19년생 때 0.065~0.081 m³에서 24년생 때 0.104~0.147 m³, 31년생 때 0.182~0.274 m³, 37년생 때(2차 간벌 직후) 0.230~0.372 m³, 42년생 때 0.285~0.580 m³까지 변화하였다. 제3간벌시험림에서는 21년생 때 0.127~0.190 m³, 32년생 때 0.254~0.415 m³, 38년생 때(2차 간벌 직후) 0.394~0.564 m³, 43년생 0.495~0.667 m³로 시간경과에 따른 개체목 재적 변화가 관측되었다(Figure 7).

Figure 7. Growth change of average tree volume by thinning intensity in each experimental site of Korean white pine plantation. Thinning represents the thinning intensity from first thinning to second thinning. The thinning symbols are the acronyms for thinning intensity; H is heavy thinning, M is moderate thinning, L is light thinning, U is un-thinning. Specific thinning treatments were performed according to Table 2. Significant differences between thinning treatments at the final measurement are indicated by different letters based on Duncan’s multiple range test (alpha=0.05).

Download Original Figure

시험구간 평균재적의 최대 차이를 살펴보면, 제1, 2간벌시험림의 경우 19년생 때 제1시험림은 0.017 m³, 제2시험림은 0.005 m³였으며, 21년생 때 제3시험림은 0.063 m³로 차이가 크지 않았다. 2차 간벌 직후 37년생 때 제1시험림 0.142 m³, 제2시험림 0.050 m³이었으며, 38년생 때 제3시험림 0.170 m³로 차이가 증가하였다. 최종 조사 시기에 간벌구간 개체목 평균재적의 최대 차이는 제1시험림 0.295 m³, 제2시험림 0.113 m³, 제3시험림 0.172 m³로 더욱 증가하였으며, 이에 대한 시험림별 간벌구간 ANOVA 분석에서 간벌구간 재적차이는 통계적으로 유의하였다(p<0.05).

간벌구별 개체목 재적의 정기평균생장량은 간벌강도가 커짐에 따라 증가하였으며, 모든 간벌시험림에서 해당 경향은 2차 간벌 이후에도 지속적으로 나타났다(Figure 8). 그러나 흉고직경 생장량과 달리 재적 생장량은 시간이 경과함에 따라 뚜렷하게 감소하지 않았으며, 임령이 큰 입목의 재적 생장량이 더 큰 것으로 관측되었다. 이는 흉고직경의 1차원적 생장과 달리 재적의 경우 누적 생장된 기존 입목재적에 3차원적으로 생장량이 추가되기 때문으로 사료된다(Jung et al., 2022). 제1, 2시험림의 재적 생장량은 19-24년 때의 0.008~0.013 m³/year이었으며, 최종 측정시기인 37-42년 때 0.010~0.029 m³/year 이었다. 제3시험림의 재적 생장량은 첫 생장기간인 21~32년 때 0.008~0.021 m³/year이었으며, 마지막 생장기간인 38~43년 때 0.020~0.026 m³/year이었다. 이와 같이 시간경과에 따라 간벌구의 재적생장이 증가하는 경향과 간벌강도가 높을수록 재적이 생장하는 결과는 유사 선행연구들의 결과와 일치하였다(Mäkinen and Isomäki, 2004a; 2004b; Wallentin and Nilsson, 2011). 따라서 적정간벌이 목재생산 목적에 부합하는 이용재적량을 더 증대시킬 수 있을 것으로 사료된다.

Figure 8. Box and whiskers plot for the periodic annual increment of tree volume by thinning intensity in each experimental site of Korean white pine plantation. Thinning represents the thinning intensity from first thinning to second thinning. The thinning symbols are the acronym for thinning intensity; H is heavy thinning, M is moderate thinning, L is light thinning, U is un-thinning. Specific thinning treatments were performed according to Table 2. Significant differences between thinning treatments at each measurement are indicated by different letters based on Duncan’s multiple range test (alpha=0.05).

Download Original Figure

임분재적은 간벌률이 높고 임분밀도가 소한 표준지에서 낮았으며, 시간경과에 따른 재적 증가량도 임분밀도 순으로 관측되었다(Figure 9). 제1, 2간벌시험림에서 간벌구별 임분재적의 범위는 19년생 때 67~166 m³/ha, 37년생 때(2차 간벌 직후) 204~488 m³/ha, 42년생 때 273~579 m³/ha이었다. 제3간벌시험림의 경우, 임분재적의 범위는 21년생 때 118~285 m³/ha, 38년생 때(2차 간벌 직후) 313~398 m³/ha, 43년생 때 393~523 m³/ha이었다. 전반적으로 임분재적은 임분밀도가 밀한 표준지에서 더 높았으며 본 결과는 임업선진국에서 주요 침엽수종을 대상으로 간벌에 따른 재적변화를 분석한 선행연구와 일치하였다(Mäkinen and Isomäki, 2004a; 2004b).

Figure 9. Growth change of stand volume per ha by thinning intensity in each experimental site of Korean white pine plantation. Thinning represents the thinning intensity from first thinning to second thinning. The thinning symbols are the acronyms for thinning intensity; H is heavy thinning, M is moderate thinning, L is light thinning, U is un-thinning. Specific thinning treatments were performed according to Table 2.

Download Original Figure

임분재적 정기평균생장량은 전체적으로 6~25 m³·ha⁻¹·year⁻¹ 범위로 관측되었는데, 간벌강도 및 임령에 따라 다양하였다(Figure 10). 이 결과는 잣나무의 지위지수 및 임령에 따라 3~8 m³·ha⁻¹·year⁻¹로 제시된 국내 임분수확표와 1~9 m³·ha⁻¹·year⁻¹로 제시된 현실림 임분수확표보다 높았다(Korea Forest Research Institute, 2012b; National Institute of Forest Science, 2016). 하지만, National Institute of Forest Science(2021)은 국내 잣나무 인공림의 간벌시험지를 대상으로 2004~2018년, 2004~2016년, 2010~2020년 사이에 2회 매목조사를 통해 변화된 임분현황을 제시하였는데, 간벌강도, 잔존본수, 임령에 따라 해당 연구의 임분재적 정기평균생장량은 7~24 m³·ha⁻¹·year⁻¹로 분석되었다. 본 연구결과와 선행연구를 종합하면, 임분재적은 임분밀도에 따라서 기존 수확표 상의 정기평균생장량보다 높을 수 있다고 사료된다.

Figure 10. Periodic annual increment of stand volume by plot in each experimental site of Korean white pine plantation. Specific thinning treatments by plot were performed according to Table 2.

Download Original Figure

제2시험림 24-37년생 기간과 제3시험림 21-38년생 기간에 무간벌구(Plot D)의 임분재적은 당시 약도간벌구(Plot C)보다 낮았는데, 이는 많은 고사목이 발생함과 임분밀도별 생장차이에 따른 영향으로 파악되었다(Figure 9). 또한, 간벌강도 및 생장시기별 임분재적 정기평균생장량은 전반적으로 강도간벌구에서 낮고 약도 및 중도간벌구의 임분재적 생장량이 높게 관측되었다(Figure 10). Mäkinen and Isomäki(2004b)에 따르면, 전체 총 임분재적은 약도간벌과 무간벌에서 큰 차이가 없을 수 있지만, 무간벌구의 자연고사를 고려하면 최소 약도간벌 수준 이상이 필요하다고 보고하였는데, 본 연구의 제2, 3시험림 임분재적 및 정기평균생장량 분석결과에서도 임분 재적의 최대 증가량을 위해 적정 간벌이 필요한 것으로 파악되었다(Figures 9, 10).

임목생장학 측면에서 간벌에 따른 개체목과 임분재적 변화를 종합하면, 간벌강도가 클수록 임분재적의 총량은 다소 낮을 수 있지만 개체목 재적생장은 촉진되는 것으로 분석되었다(Figures 7, 9). 따라서 일정한 산림경영계획 기간 내 형질이 우수한 대경재 생산을 촉진하기 위해서는 1차 및 2차 간벌의 적정 산림시업이 필요한 것으로 판단된다(Hynynen et al., 2002; Huuskonen and Hynynen, 2006; Huuskonen et al., 2014; Niemistö et al., 2018). 또한, 자연고사로 인한 임목손실을 고려하면 최소 약도간벌 수준 이상의 산림시업이 임분재적 증가량 측면에서도 더욱 용이할 것으로 사료된다(Mäkinen and Isomäki, 2004b). 추후 잔존 입목본수, 임분 흉고단면적, 상대밀도(Relative density), 임분밀도지수(Stand density index), 임목축적도(Stocking percent)와 같은 다양한 밀도진단 및 평가를 통해 임분관리 목표별 최적 산림시업이 제시될 수 있도록 지속적인 연구가 필요하다(Lee and Choi, 2019; 2020; Jung, 2022).

제1, 2시험림의 고사목 발생은 임분밀도가 높았던 무간벌구에서 많이 관측되었다. 특히 전체 누적 고사본수가 간벌구 Plot A, B, C에서 350 본/ha 이하였던 것에 반하여, 무간벌구 Plot D에서는 900 본/ha 이상으로 매우 많이 발생하였다(Figure 11). 이와 같은 경향은 제3간벌시험림에서도 나타났는데, 간벌구 Plot A에서는 고사목 누적본수가 전혀 발생하지 않았던 것에 대비하여, 무간벌구 Plot D에서는 600 본/ha이 발생하였다. 따라서 무간벌구에서 간벌구보다 고사목이 많이 발생한다고 보고한 선행연구들과 본 연구결과가 일치하였다(Mäkinen and Isomäki, 2004a; 2004b; Zhang et al., 2013).

Figure 11. Mortality by plot over time in each experimental site of Korean white pine plantation. Specific thinning treatments by plot were performed according to Table 2.

Download Original Figure

특히, 제1,2시험림은 2차 간벌시기에도 Plot D가 무간벌로 유지되었는데, 이에 따라 37-42년생 때 고사목의 발생은 더욱 많았다. 반면에 제3시험림에서는 2차 간벌 시 강도간벌구를 제외한 모든 간벌구에서 간벌이 시행됨에 따라 이후 최종 생장기간인 38-43년생 때 고사목은 거의 발생하지 않았다. 본 연구결과, 고사목 발생은 간벌시업에 따른 임분밀도에 크게 영향을 받는 것으로 사료된다.