정주지 단풍나무의 탄소저장량 추정 상대생장식

본 연구는 경상국립대학교 종합농장(35°09'N, 128°04'E) 내 조경용으로 식재되고 정상적으로 생육하고 있는 단풍나무 20본을 대상으로 하였다. 선정된 표본목은 접목이 실시된 홍단풍나무와 청단풍나무가 고르게 섞여 있었으며, 가지치기나 초두부 상처 부위의 유합 상태를 고려할 때[Figure 1(a)] 최소 5년 이전에 수형조절을 실시한 것으로 사료되었다.

Figure 1. A pruned and grafted sample tree (a) with branch and leaf biomass (b).

Download Original Figure

탄소저장량 추정 상대생장식 개발을 위한 표본목은 ‘산림 바이오매스 및 토양 탄소 조사·분석 표준’(Korea Forest Research Institute, 2007)에 의거 흉고직경이 고르게 분포하도록 하고 수형 등을 고려하여 총 20본을 선정하였다. 선정된 표본목은 잎의 생장이 완료된 2021년 5월 하순 지상부 0.2 m 높이를 기계톱을 이용하여 벌채하고 줄자를 이용하여 수고, 수관폭, 지하고를 측정하였다. 표본목의 연령은 13년에서 26년생까지 분포하였으며, 평균 연령은 17년이었다. 평균 흉고직경과 수고는 각 11.7 cm와 4.9 m였고, 평균 수관폭과 지하고는 각 3.3 m와 1.2 m 정도로 나타났다(Table 1).

벌도된 표본목의 줄기로부터 가지를 분리하고 줄기의 0.2 m에서 1.2 m까지는 1 m, 1.2 m 이후부터는 2 m 간격(3.2 m, 5.2 m, 7.2 m …)으로 절단하였다. 줄기의 생중량은 10 g 단위까지 측정이 가능한 전자저울을 이용하였으며, 건중량과 탄소 농도 측정을 위해 각 절단부의 원구로부터 5 cm 두께의 원판 시료를 채취하였다. 가지 부위는 가지와 잎으로 분리하고[Figure 2(b)] 5 g 단위까지 측정이 가능한 전자저울로 생중량을 측정한 후 건중량 환산과 탄소 농도 분석을 위하여 일부 시료를 채취하였다. 뿌리의 굴취는 굴삭기를 이용하였고 뿌리에 있는 흙을 최대한 제거한 다음 전자저울을 이용하여 10 g 단위까지 측정하였다. 뿌리 분석 시료는 굵은 뿌리와 잔뿌리의 비율을 고려하여 고루 섞어 채취하였다.

Figure 2. Scatter plots of tree component biomass in Acer palmatum (**P<0.001, ***P<00.0001).

Download Original Figure

실험실로 운반된 원판 시료는 85°C 온풍 건조기에서 항량에 도달할 때까지 건조하였다. 건조한 원판 시료는 수피 칼을 이용하여 목질부와 수피로 분리하고 건중량을 환산하였다. 부위별 건중량 시료는 분쇄기로 곱게 갈은 후 40메쉬(0.420 mm) 체를 이용하여 탄소 농도 분석용 시료를 조제하였다. 부위별 탄소 농도 측정은 원소분석기(Vario Macro Cube, Langensebold, Germany)를 이용하였고, 탄소 농도 값에 해당 부위의 건중량을 곱하여 각 부위의 탄소저장량으로 환산하였다.

단풍나무 식재지 주위의 토양특성 조사를 위해 400 cm³ 원통형 채토기를 이용하여 0~10 cm 깊이의 토양 10점을 채취하였다. 채취한 토양 시료는 65°C의 온풍건조기에서 항량에 도달할 때까지 건조 후 2 mm 체를 이용하여 토양 물리·화학적 특성 분석용 시료를 조제하였다. 토양 입경 분포는 비중계법, 토양 탄소 농도는 원소분석기(Vario Macro Cube, Langensebold, Germany)를 이용하였다. 유효인산, 칼륨, 칼슘, 마그네슘은 토양 2.5 g을 암모늄아세테이트 용액 55 ml를 첨가하고 16시간 동안 진공 추출 후 ICP (Perkin Elmer Optima 5300, USA)로 분석하였다(Kalra and Maynard, 1991). 토양 pH와 EC는 토양시료 5 g과 증류수 25 ml의 혼합용액을 이온전극(Istec 735-P, Korea)과 EC 메터(Orion 3 star conductivity, Singapore)로 측정하였다.

본 연구에서는 DBH와 DBH²×H 를 독립변수(X)로 하고, 줄기 목질부(stem wood), 줄기 수피(stem bark), 가지(branch), 잎(leaf), 뿌리(root), 지상부+지하부(aboveground+ root) 탄소저장량을 종속변수(Y)로 하는 상대생장식(log₁₀Y=a+b×log₁₀X)을 개발하였다. 한편, 본 연구에서 조사된 표본목은 연령과 식재지가 다양하고, 표본목 수가 20본 정도로 적어 추정 변수에 대한 정규분포의 가정과 상대생장식의 잔차 이분산성의 문제점이 발생할 수 있다(Bi et al., 2004; Wang et al., 2018). 이러한 문제점을 해결하기 위해 SAS의 Proc Model Procedure(SAS Institute Inc., 2003)에 도구변수(instrumental variable)로 표본목의 연령과 수관폭을 포함하여 GMM 방법에 의한 상대생장식을 개발하였다. 또한, GMM 상대생장식과 비교를 위해 탄소저장량 추정에 일반적으로 사용하는 Log 상대생장식을 SAS의 Proc Reg Procedure(SAS Institute Inc., 2003)로 계산하였고, Log 변환 상대생장식에서 나타날 수 있는 편향(bias)에 대한 보정 인자(correction factors, CF)도 제시하였다(Sprugel, 1983). 표본목의 부위별 탄소농도는 연령을 공변이로 하여 SAS의 Proc GLM Procedure로 유의성 검정 후(P<0.05), Least square means에 의해 평균 탄소 농도의 차이를 비교하였다.

조사지의 용적밀도는 0.81 g cm⁻³(Table 2)로 본 연구 조사지와 인접한 조경용 느티나무 식재지의 용적밀도 0.87 g cm⁻³(Ha et al., 2022)과 유사하였다. 토성은 사양토, 석력 함량, 토양 pH, 전기전도도는 각각 0.24 g g⁻¹, 6.10, 311 ㎲ cm⁻¹로 나타났다(Table 2). 유기탄소 농도는 1.23%로 진주시 4개의 도시공원 토양 0~10 cm 깊이에 유기탄소 농도 1.12%(Ahn et al., 2022)와 큰 차이는 없었다. 교환성 칼륨, 칼슘, 마그네슘은 각각 0.24, 5.81, 0.75 cmolc kg⁻¹로 타 연구에서 보고된 단풍나무 식재지 토양(Jo and Ahn, 2012)의 0.23, 1.19, 0.60 cmolc kg⁻¹에 비해 칼슘 함량은 높으나 칼륨과 마그네슘 함량은 큰 차이가 없었다.

표본목의 목질부, 뿌리, 지상부+지하부 총건중량은 흉고직경의 증가와 함께 강한 양의 상관관계를 보이고 있으나, 줄기 수피, 가지, 소지, 잎 등은 표본목 사이에 변이 폭이 크게 나타나 상관관계가 비교적 약하게 나타났다(Figure 2). 이는 가지치기나 초두부 절단 같은 수관부의 훼손에 의해 이들 부위의 건중량이 흉고직경 생장과 밀접한 관계가 나타나지 않은 것으로 사료된다. 한편, 본 연구에서 조사된 단풍나무 표본목 당 바이오매스 총량은 12 kg~158 kg로, 뿌리 바이오매스가 높게 나타난 일부 표본목을 제외하고 Jo and Ahn(2012)이 보고한 흉고직경 4.9~19.6 cm 범위의 조경용 단풍나무 단목 당 바이오매스 5~133 kg에 속하고 있다.

표본목의 부위별 바이오매스 분포 비율은 뿌리(41.2%), 줄기 목질부(35.4%), 가지(17.1%), 잎(4.0%), 줄기 수피(2.3%) 순이었다(Figure 2). 타 연구에서 보고된 조경용 단풍나무, 느티나무, 왕벚나무의 평균 바이오매스 분포 비율(Jo and Ahn, 2012)은 줄기(41~44%), 뿌리(26~30%), 가지(21~26%), 잎(4~7%)으로 본 연구 결과와 비교할 때 뿌리의 비율은 높고, 줄기 목질부와 가지 바이오매스는 낮게 나타났다. 본 연구에서 줄기 비율이 낮은 것은 수관을 이루는 줄기, 가지, 잎 등의 일부가 제거되어 뿌리 바이오매스 분포비가 상대적으로 높게 나타났거나, 접목에 따른 지상부와 지하부 바이오매스 축적량의 차이로 사료된다. 한편, 잎 바이오매스는 총 바이오매스의 4%로 정주지에 식재된 느티나무 4.2%, 은행나무 4.4%(Ha et al., 2022)와 큰 차이는 없었다. 도시에 식재된 수목의 바이오매스 분포비율은 생육지의 입지환경, 수형 조절을 위한 가지치기 강도(pruning intensity) 및 초두부 제거 작업 같은 관리방법에 따라 차이가 나타날 수 있는 것으로 알려져 있다(Park and Kang, 2010; Ha et al., 2022). 표본목의 바이오매스 분포에 대한 변이계수는 가지>잎>줄기 수피>뿌리>줄기 목질부 순으로(Figure 3) 수관부를 구성하는 가지와 잎의 변이 폭이 크게 나타났다.

Figure 3. Proportion of biomass distribution in tree components in Acer palmatum. Vertical bars indicate standard errors. CV: coefficient variation.

Download Original Figure

단풍나무 목질부의 탄소농도는 49.8%로 잎 47.2%, 가지 48.8%, 줄기 수피 46.5%에 비해 높은 값을 보였으나, 뿌리의 탄소 농도 49.3%와는 유의적 차이가 없었다(Figure 4). 수목의 부위별 탄소 농도는 리그닌 함량 등과 함께 생리적·화학적 요인에 따라 차이가 있는 것으로 알려져 있다(Herrero de Aza et al., 2011; Thomas and Martin, 2012; Morhart et al., 2016). 본 연구에서 나타난 목질부의 높은 탄소 농도는 이 부위의 리그닌 함량이 타 수목 부위에 비해 높기 때문일 수 있다(Novaes et al., 2010). 그러나 잎이나 수피의 탄소 농도는 타 수목 부위에 비해 낮은 값을 보이고 있으며, 이는 잎이나 형성층을 포함하는 수피는 주로 대사작용에 관여하는 기관으로 조직 내 양분 함량이 타 수목 부위에 비해 높게 나타나기 때문으로 사료된다. 예를 들면, 수목의 조직 내 양분 함량과 탄소 농도 사이는 부의 상관이 보고된 바 있다(Kim, 2019). 한편, 본 연구에서 조사된 단풍나무의 평균 탄소 농도는 48.4%로 IPCC에서 보고된 온대지역 낙엽활엽수의 평균 탄소농도 48% (IPCC, 2006)와는 차이가 없었다.

Figure 4. Carbon concentration (%) of tree components (SW: stem wood; SB: stem bark; BR: branch; LF: leaf; RT: root) in Acer palmatum. The box represents the median and the 25^th and 75^th percentiles. × represents the arithmetic mean. The solid lines extend to 1.5 of the interquartile range and the values outside this range are indicated by circle. Vertical bars indicate standard error. Different letters on the bar of each tree component represent a significant difference at P<0.05.

Download Original Figure

수목 부위별 탄소저장량(줄기 목질부, 줄기 수피, 가지, 잎, 뿌리, 지상부, 지상부+지하부 총량)은 DBH를 독립변수로 하는 Log 상대생장식의 결정계수가 0.3858~0.8450, GMM 상대생장식은 0.3027~0.8087로 두 방법 사이에 큰 차이는 없었으나 Log 상대생장식의 결정계수 값이 약간 높았다(Table 3). DBH²×H를 독립변수로 하는 Log 상대생장식의 결정계수는 0.5619~0.9069, GMM 상대생장식의 결정계수는 0.5055~0.8837로 DBH²×H가 DBH를 독립변수로하는 상대생장식에 비해 결정계수 값이 약간 높게 나타났다. 유사한 결과로 Wang et al.(2018)은 중국 갈참나무(Quercus aliena Bl. var. acuteserrata) 바이오매스 추정을 위한 Log 상대생장식의 결정계수는 0.9516, GMM 상대생장식은 0.9414로 두 방법 사이에 큰 차이 없음을 보고한 바 있다. 본 연구 결과는 탄소저장량 추정 상대생장식에 있어서 표본목 수가 적거나, 탄소저장량 추정 관련 변수들이 정규분포하지 않은 경우, Log 상대생장식으로부터 발생할 수 있는 잔차의 이분산성의 해결을 위해 도구변수로 표본목의 연령과 수관폭을 사용한 GMM 상대생장식의 이용이 가능함을 시사한다.

한편, 본 연구의 상대생장식 결정계수는 중부지역 조경수 재배지의 단풍나무(Jo and Ahn, 2012)나, 본 연구 조사지와 인접지에 조경수로 식재된 느티나무, 은행나무, 왕벚나무(Ha et al., 2022)의 DBH 또는 DBH²×H 독립변수의 탄소저장량 추정 상대생장식의 결정계수 0.97~0.99에 비해 낮은 값을 보였다. 이는 타 연구의 경우 표본목의 수가 5~11본으로 표본목 수가 적어 회귀식의 설명력이 높게 나타났거나(Chen et al., 2021), 표본목의 흉고직경 사이 범위가 넓어 직경 생장에 따른 임목 부위별 탄소저장량의 상대생장 관계가 분명하였기 때문으로 사료된다. 그러나 DBH 또는 DBH²×H를 독립변수로 한 Log나 GMM 상대생장식 모두 가지나 잎의 결정계수는 0.3027~0.6380으로 타 임목 부위에 비해 매우 낮은 값을 보였다. 이는 도시에서 생육하는 수목은 수형조절을 위한 가지치기와 초두부 제거 작업과 같은 인위적인 관리로 인하여 수관을 구성하는 잎이나 가지의 변이 폭이 크게 나타나 해당 부위의 상대생장식 결정계수 값이 낮게 나타나는 결과로 해석된 바 있다(Ha et al., 2022). 본 연구에서도 가지와 잎의 변이계수는 34% 이상으로 타 수목부위 21% 이하에 비해 높은 값을 보였다(Figure 3). 한편, Log 변환 상대생장식에 대한 보정 인자(CF)의 경우 줄기 목질부, 줄기 수피, 뿌리, 총 탄소저장량 등은 1.05보다 작은 값을 보여 탄소저장량 추정 시 발생할 수 있는 변이(variation)가 작지만, 가지나 잎 등은 1.10 이상으로 크게 나타나 이 부위의 탄소저장량 추정 상대생장식에 CF의 사용이 필요한 것으로 나타났다.