서 론

사스레피나무(Eurya japonica Thunb.)는 중국, 대만, 일 본, 인도 등에 분포하며, 우리나라 제주도 해안가와 산지, 남해 안 도서 및 남부지역에 흔히 자생하는 상록 활엽 관목 또는 소교목이다(Brickell, 2008). 반음수로 양지에서 천연하종이 잘 되며, 5~8℃ 이상에서 월동하고, 13~21℃에서 잘 자란다. 내한성은 약하나 건조하고 척박한 사질양토에 잘 자라며 내염 성, 내조성이 강하고 수분 요구도가 적어 남부 지방의 절개지 혹은 사방지 조림에도 쓰인다. 또한 내공해 수종으로 공원수, 정원수, 생울타리로도 식재된다(Chang et al., 2011). 조경수 및 약용으로 쓰이는 것 외에 사스레피나무는 환경부에서 국가 기후변화 생물지표종으로 지정되어 관리되고 있다(Ministry of Environment, 2010).

한반도의 식물대는 아고산림, 냉온대림, 난온대림, 한대림으 로 구분되는데(Yim and Kira, 1975;1976), 사스레피나무와 같은 상록활엽수가 생육하는 난온대 기후대는 연평균기온이 14℃ 이상, 한랭지수 –10℃․월 이하, 강수량 900~1,500mm로, 동쪽의 북위 35°30´과 서쪽의 북위 35°를 연결하는 선의 남쪽 이고 남해안 및 도서지방, 제주도(해발 500m), 울릉도(해발 600m 이하) 지역이 해당된다. 그러나 온실가스 배출량의 급증 으로 30년간 연속해서 지구 표면의 온난화가 심화되어 난온대 림은 계속해서 북상하고 있다. 식물의 분포 및 생육지의 확대 ·축소에 영향을 미치는 요인은 개체군 간의 경쟁, 외래종에 의 한 침입, 지형의 변화, 인간의 간섭 등 여러 가지가 있으나 대다 수 식물이 북상하는 가장 큰 요인은 기후변화이다. 이러한 기후 변화는 생물종의 분포 및 서식환경에 영향을 미치고 상록활엽 수의 종다양성 및 분포에도 변화를 불러일으킨다(Yun et al., 2011). 상록활엽수의 종다양성은 위도가 증가됨에 따라 감소 하는데(Choi, 2013), 이는 위도와 경도에 따라 난대성 상록활 엽수 분포에 영향을 주는 최한월최저기온, 겨울철 강수량, 여름 철 평균기온, 여름철 강수량 등과 같은 기후요소들이 달라지기 때문이다(Nakao et al., 2014).

우리나라 천연 상록활엽수림 분포 면적은 약 9850ha로 이는 한반도 전체 산림면적의 1.53%로 면적이 협소하며 접근이 어 려운 도서지방 등에 국지적으로 잔존하고 있다(Park et al., 2015). 상록활엽수는 진화론적 측면, 지리학적인 경계 구분, 평균기온 예측뿐만 아니라 생물유전자원 및 서식지 보존 측면 에서 중요성을 가지고 있다(Oh and Choi, 1993;Oh and Cho, 1994). 이에 따라 Uyeki(1941)는 기후변화에 따른 상록활엽 수의 서식환경 및 분포 변화를 예측하기 위해 북방한계선을 설정하였다. 한반도 난온대 상록활엽수의 북방한계선을 설정 하기 위해 난온대 상록활엽수 64종을 선정하였는데, Uyeki (1941)가 선정한 한반도 난온대림의 대표 수종으로는 가시나 무, 녹나무, 아왜나무, 붓순나무, 육박나무, 사스레피나무 등이 있다. 이 64종의 상록활엽수 분포를 조사하여 대청도-변산-영 암-죽도를 잇는 난대성 상록활엽수의 한반도 내 북방한계선을 1941년 지정하였다. 이 중 사스레피나무는 다른 상록활엽수종 에 비해 고위도에 분포하고 있다(Uyeki, 1941;Yun et al., 2011). Yun et al.(2017)은 종 분포모델 중 CT-Model 분석을 이용하여 기후변화에 따른 사스레피나무의 생육지를 예측하였 다. 현재보다 미래 기후에서 사스레피나무의 생육지가 북상할 것이며, 한반도 내륙으로 분화 확대 가능성이 다른 난온대 상록 활엽수종에 비해 높을 것이라 예측하였다. 따라서 기후변화에 따른 사스레피나무의 분포지 변화는 식물군락을 구성하는 종조 성에도 영향을 미칠 것으로 예측된다.

본 연구 대상지인 부산광역시는 대한민국 동남부 해안에 대 한민국 제2의 도시이자 최대의 항구도시이다. 근해의 난류 영 향을 받는 해양성 기후로 여름과 겨울의 기온 차가 내륙지역에 비해 비교적 적은 편이다. Lee et al.(2007)은 부산이 식물지리 분포상 남해안도서분포선을 잇는 최고 북위 지역으로서 상록활 엽수의 생육이 가능한 지역이라 하였다. 부산광역시의 30년 평균기온 평년값이 13.5℃(1910-1939년)에서 14.9℃(1989- 2018년)로 1.4℃ 상승하여 상록활엽수의 식생구조 및 북상에 변화가 있을 것이라 예측된다.

난온대 상록활엽수림에 관한 선행연구로는 Lee et al.(2007) 이 부산광역시를 대상으로 현장조사 및 문헌조사를 통해 상록 식물의 종조성 및 분포지를 밝혀 학술적･자원적 가치가 있는 종을 보호하기 위한 연구가 있으며, Yun et al.(2011)은 식물 분포 자료 및 기후자료를 활용하여 CT-Model 분석을 통해 기후변화에 따른 한반도 난온대 상록활엽수 분포 변화 및 분포 적지를 밝혔다. Park et al.(2018)은 한국, 일본, 중국, 대만을 대상으로 문헌조사 및 식생조사를 하여 CCA 분석을 통해 난온 대 상록활엽수림의 천이계열과 극상림을 고찰하였다. 이외에 도 난온대 상록활엽수림의 식생구조, 식생천이계열, 분포지 예 측, 복원방법 등에 관한 다양한 연구들(Oh and Choi, 1993;Oh and Kim, 1996;Park and Oh, 2002)이 진행되어 왔다. 그러나 난온대 점이지대로서 중요한 가치가 있는 부산광역시를 대상으로 하는 난온대 상록활엽수종에 관한 연구는 부족한 실 정이다. 특히 사스레피나무를 중심으로 한 연구 및 식생구조에 관한 상세연구는 전무하다. 이에 본 연구에서는 부산광역시를 대상지로 사스레피나무를 중심으로 세부 식생유형을 분류하여 생태적 특성 및 식생구조 분석을 하고자 한다.

연구방법

1. 연구대상지 선정 및 개황

연구 대상지인 부산광역시는 난온대 기후대에 위치하고 있 으며, 산지 면적은 353.86km²로 전체 면적 769.89km² 의 약 46%를 차지하고 있다. 부산광역시는 절반 이상 해안선을 끼고 있어 해양성 기후의 영향을 많이 받는다. Yun et al.(2011)에 의하면 상록활엽수의 분포를 규정하는 주요 기후요인 중 최한 월최저기온(TMC)이 –5.95℃ 이상을 난온대 상록활엽수의 분 포 적지라 하였다. 부산광역시의 2019년 연평균기온은 15.8℃, 강수량 1,778.6mm이며, 최한월평균기온 4.5℃,최한월최저기 온(TMC)은 –4.4℃로 식물지리분포상 최고 북위지 난온대 점 이대로 상록식물의 생육이 적합한 지역이다. 연구대상지 내 조 사구는 해안형과 내륙형으로 사스레피나무가 출현하는 곳을 중심으로 총 89개소의 조사구를 설치하여 식생조사를 실시하 였다.

2. 조사 및 분석 방법

1) 식생 및 환경요인 조사

식생조사는 대상지 내에 10m×10m 크기의 방형구 89개소 를 설치하고 교목층, 아교목층, 관목층으로 층위를 구분하여 수관층위별로 매목조사를 실시하였다(Park, 1985). 조사구에 출현하여 상층수관을 이루는 흉고직경(DBH) 2cm 이상의 교 목층과 수고 2m 이상, 교목층 이하 수목을 아교목층, 수고 2m 이하의 수목을 관목층으로 구분하였다. 10m×10m 방형구 내 교목층과 아교목층 수목의 흉고직경을 측정하였으며, 관목층 은 방형구의 가장자리 5m×5m 소방형구 1개소를 설치하여, 출현한 수목의 수관폭(장변×단변)을 조사하였다. 각 조사지의 환경요인은 해발고, 사면방향, 경사도, 식피율을 조사하였다.

2) 식물군락구조 조사

조사한 식생자료를 바탕으로 식물군락별 특성을 알아보기 위해 89개 조사구의 군락분류는 TWINSPAN(Two-way indicator species analysis)에 의한 군락분석(classification analysis)(Hill, 1979b) 및 DCA(Detrended correspondence analysis) ordination(Hill, 1979a)분석을 실시하여 분류하고 군락별 분포특성을 파악하였다. 분류된 군락 간 유사도를 알아 보기 위해 유사도지수를 분석하였으며, 식생조사 자료를 토대 로 각 수종의 상대적 우세를 비교하기 위하여 Curtis and McIntosh(1951)의 중요치(Importance Value; I.V.)를 통합 하여 백분율로 나타낸 상대우점치(Brower and Zar, 1977)를 수관층위별로 분석하였다. 상대우점치(Importace Percentage; I.P.)는 (상대밀도+상대피도)/2로 계산하였으며, 평균상대우점 치(Mean Importance Percentage; M.I.P.)는 층위별로 개체 간 크기(현존량)를 고려하여 가중치를 부여한 (교목층 I.P.×3+ 아교목층 I.P.×2+관목층 I.P.×1)/6의 식을 이용하였다(Park, 1985). 분류된 군락별로 Shannon의 수식(Pielou, 1975)을 이 용하여 종다양도(Species Diversity, H′), 균재도(Evenness, J′), 우점도(Dominance, D)를 계산했고, 단위면적당(100㎡) 종수 및 개체수를 분석하였다. 또한, 주요 종을 중심으로 산림 천이 양상을 간접적으로 유추할 수 있는 흉고직경급별 분석을 하였다(Harcombe and Marks, 1978).

결과 및 고찰

1. 군락분류 및 유사도지수 분석

부산광역시 일대를 중심으로 한 89개 조사구에 Classification 분석 중 TWINSPAN 분석을 통해 각 조사구에서 출현하는 수종들 중 환경요인을 간접적으로 반영하는 지표종(indicator species)을 중심으로 군락을 분류하였다(Manhas et al., 2009).

제 1단계(level 1) Division 1에서는 졸참나무(-)와 동백나무 (+)를 중심으로 분리가 되었으며, 제 2단계(level 2) Division 2에서는 곰솔(+)로 나누어졌으며, 전체적으로 3개의 군락으로 나누어졌다. 군락 Ⅰ은 24개의 조사구, 군락 Ⅱ는 49개의 조사 구, 군락 Ⅲ은 16개의 조사구로 나누어졌다(Figure 2).

TWINSPAN에 의한 군락분류를 보완하고, 군락의 분포를 알아보기 위해 조사구 간의 상이성을 바탕으로 조사구를 배치 하는(Orloci, 1978) DCA기법을 이용(Choi and Kang, 2006) 하여 전체 89개 조사구에 대해 분석을 실시하였다. DCA분석 결과 제1축과 2축의 eigenvalue가 각각 0.523, 0.317로 4개축 전체 합 1.156의 72.66%에 해당하여 total variance에 대한 집중률이 높은 편으로 나타났다.

TWINSPAN과 DCA 분석 결과와 더불어 조사구별 종 조 성과 평균상대우점치를 고려하여 89개소의 조사구는 세 개의 군락으로 나누어졌다. 군락 Ⅰ은 졸참나무-소나무-사스레피나 무 군락, 군락 Ⅱ는 곰솔-사스레피나무 군락, 군락 Ⅲ은 곰솔- 동백나무 군락으로 분류되었다(Figure 3).

Ordination 분석과 함께 분리된 3개의 군락에 대해 유사도 지수를 분석하여 백분율로 나타내었다.(Table 1). 유사도지수 가 20% 미만일 때, 이질적인 집단이고, 80% 이상일 때 생태적 으로 유사한 동질적인 집단이다(Whittaker, 1956). 군락 간 유사도지수는 군락 Ⅱ와 군락 Ⅲ이 52.20%로 가장 높게 나타 났는데 이는 두 군락의 교목층에서 곰솔의 평균상대우점치가 높게 나타났기 때문인 것으로 사료되었다. 반면에 군락 Ⅰ과 군락 Ⅲ은 사스레피나무가 출현하는 군락이지만 13.31%로 가 장 낮은 유사성을 나타내는데 이는 군락 Ⅲ의 교목층과 아교목 층의 우점종이 각각 곰솔과 동백나무로 사스레피나무의 평균 상대우점치가 다른 두 군락에 비해 상대적으로 낮기 때문인 것으로 사료되었다.

2. 군락별 대상지 개황

TWINSPAN과 DCA를 적용한 결과 전체 조사구는 3개 군 락으로 분류되었으며, 분류된 3개 군락의 대상지 개황을 나타 내었다(Table 2). 주로 해안가에 위치하였으며 해발고도는 17~370m, 주된 사면 향은 남서향, 경사도는 3~34°이다. 군락 Ⅲ의 경우 해발고도는 낮지만 경사도는 다른 두 군락에 비해 상대적으로 높은 것을 확인할 수 있는데 이는 해발고도가 낮은 동백섬에 위치하였지만 경사가 20° 정도인 급경사지에 사스레 피나무가 주로 분포했기 때문이다. 군락별 식생의 경우 교목층 평균흉고직경은 약 25~35cm, 아교목층의 평균흉고직경은 약 6~8cm, 관목층의 식피율은 20~37%로 나타났다.

3. 식생구조 분석

분류된 세 개의 군락에 대해 층위별 상대우점치(I.P.) 및 평 균상대우점치(M.I.P.)를 분석하고(Table 3), 높은 상대우점치 를 가진 우점종을 중심으로 군락구조를 분석하였다.

군락 Ⅰ은 졸참나무-소나무-사스레피나무 군락으로 교목층 에서는 졸참나무(I.P. 26.8%)와 소나무(I.P 28.4%)가 활발한 경쟁관계인 상태로 출현하고, 아교목층에서 사스레피나무(I.P. 33.5%)가 우점하고 있었다. 소나무는 아교목층과 관목층에서 출현하지 않은 반면에 졸참나무와 굴참나무 등 낙엽성 참나무 류가 활발히 출현하였다. 산림천이 양상에서 소나무가 선구수 종이라는 점을 감안하였을 때, 향후 참나무류의 세력이 점차 확장될 것이라 판단되며(Cho and Choi, 2002;Oh et al., 2014), 더 나아가 극상수중인 개서어나무(I.P. 12.4%)가 낙엽 성 참나무류 다음으로 높게 출현하였기 때문에 낙엽성 참나무 류 또는 개서어나무가 우점하는 식생으로 천이될 것으로 사료 된다(Kim, 2012).

군락 Ⅱ는 곰솔-사스레피나무 군락으로 교목층에서는 곰솔 (I.P. 85.3%)이 우점하며, 아교목층과 관목층에서는 각각 I.P. 57.9%, 11.2%로 사스레피나무가 우점하고 있었다. 관목층에 서 졸참나무, 때죽나무 등 낙엽활엽수류가 출현하고 있었으나 자연재해, 간벌, 밑깎기작업 등의 인위적인 간섭이 이루어지지 않는 이상 곰솔-사스레피나무의 세력이 유지될 것으로 예측되 었다.

군락 Ⅲ은 곰솔-동백나무 군락으로 교목층에서는 곰솔(I.P. 100.0%)만이 단일 수종으로 나타났으며, 아교목층에서는 동 백나무(I.P. 85.1%)가 우점하고 있었다. 이는 군락 Ⅲ의 조사 구가 위치한 대상지가 부산광역시 지정 기념물 제 46호인 동백 섬으로 동백나무 보식 등 인위적인 관리가 이루어지고 있었기 때문이다(Kim et al., 2014). 반면, 관리가 이루어지지 않고 경사가 급한 지역에는 사스레피나무가 우세하고 있었다. 군락 Ⅲ의 아교목층에서 높은 상대우점치를 보이는 동백나무와 사 스레피나무는 개체군간 공간적 유전 분포의 패턴이 매우 유사 하다(Oh, 1998). 따라서 동일한 생태적 지위를 가지며 경쟁관 계인 사스레피나무와 동백나무는 이전과 같은 지속적인 관리가 없을 경우, 배수가 양호하고 토심이 깊은 곳에 자라는 동백나무 는 쇠퇴하고, 관리되지 않은 건조하고 척박한 환경에 적응하여 생장하는 사스레피나무가 우점할 것으로 사료되었다. Park et al.(2018)은 한국 난온대 지역의 식생천이 중간단계 수종으로 후박나무, 생달나무 등을 꼽았으며, 양수인 곰솔보다 음수성을 띄는 상록활엽수종의 생태적 지위가 높다고 밝혔다. 따라서 군 락 Ⅲ의 관목층에서 동백나무, 사스레피나무와 더불어 후박나 무가 자라고 있어 인위적인 간섭이 없을 경우 동백나무가 쇠퇴 하고 사스레피나무가 우점하며 더 나아가 난온대상록활엽수림 으로 천이가 진행될 것으로 판단되었다.

이상 사스레피나무를 중심으로한 식생구조를 분석한 결과, 다른 두 군락에 비해 해발고도가 높고 대부분의 조사구가 내륙 에 위치한 군락 Ⅰ에서는 교목층에 낙엽성 참나무류가 소나무 경쟁관계인 상태에서 아교목층에 사스레피나무가 우점하였다. 군락 Ⅱ, 군락 Ⅲ의 조사구는 해발고도가 낮고 해안가에 인접 한 곳에 위치하여 교목층에서 곰솔이 우점하고 아교목층에 사 스레피나무가 출현하였다. 사스레피나무는 해발고도가 높은 내륙에서는 식생천이가 진행됨에 따라 낙엽성 참나무류 및 소 나무와 함께 출현하고 해발고도가 낮은 해안가에는 곰솔과 함 께 출현하였다.

4. 종수 및 개체수, 종다양도 분석

세 개 군락의 단위면적(100㎡) 당 평균 출현 종수 및 개체 수 분석을 층위별로 실시하였다(Table 4). 세 개 군락은 모두 다수의 조사구가 설정되어 포함된 조사구의 평균값을 산출하 였다. 전반적인 조사구의 층위별 평균 종수 분석결과 교목층 1.63±0.92종, 아교목층 4.24±2.37종, 관목층 8.04±3.09종으 로 전체 조사구당 11.56±3.64종이 출현하였다. 군락 Ⅰ은 졸참 나무-소나무-사스레피나무 군락으로 교목층 2.42±0.97종, 아 교목층 5.33±2.93종으로 다른 두 군락보다 출현 종수가 높게 나타났다. 이는 산림천이 중간단계로 교목층과 아교목층에서 소나무와 낙엽성 참나무류가 활발한 경쟁 상태로 다양한 수종 이 출현하기 때문으로 판단되었다. 군락 Ⅱ는 곰솔-사스레피나 무 군락으로 관목층에서 8.84±3.23종으로 세 군락 중 가장 많 은 종이 출현하였다. 이는 교목층과 아교목층에서 각각 곰솔과 사스레피나무가 출현하고 있었으나 자연재해 및 병충해 등 다 양한 교란요인으로 인해 숲 틈으로 광량이 증가하며 관목층에 서 졸참나무, 때죽나무 등 낙엽활엽수가 다수 출현하기 때문으 로 사료되었다(Yun et al., 2006). 또한 조사구 개수가 49개로 다른 두 군락에 비해 상대적으로 많아 출현 종수 또한 많은 것으로 사료되었다. 군락 Ⅲ은 곰솔-동백나무 군락으로 대부분 이 인위적인 관리가 이루어지는 지역으로 교목층은 단일수종으 로 곰솔이 출현하였고, 아교목층은 동백나무가 우점하는 가운 데 일부 비관리지역은 사스레피나무가 출현하였다. 관목층은 사스레피나무를 포함한 후박나무, 사철나무 등 상록활엽수가 다양하게 출현하였다.

전반적인 조사구의 총 평균출현 개체수는 646.75±1,893.19 개체이며, 각 층위별 평균출현 개체수는 교목층 6.02±3.17개 체, 아교목층 17.19±11.12개체, 관목층에서 623.73±1,895.24 개체가 출현하였다. 총 평균출현 개체수 및 관목층 편차가 큰 이유는 마삭줄 때문으로 마삭줄을 제외한 후 분석하였을 때 관목층 평균출현 개체수는 136.27±184.78개체, 전체 평균출현 개체수는 159.29±183.75개체로 나타났다.

3개 군락의 Shannon의 종다양도지수(H′), 균재도(J′), 우점 도(D)를 분석하였다(Table 5). 각 군락별 종다양도지수는 단위 면적(100㎡)을 기준으로 산출하였으며, 군락 Ⅰ,Ⅱ,Ⅲ 모두 다수의 조사구가 포함되어 있어 조사구별 종다양도지수의 평균 값을 산출하였다. 종다양성은 군집의 생태천이 발달과 밀접한 관계를 가지고 있으며, 군락의 안정도와 성숙도는 생태천이의 진행정도를 보여준다(Yun et al., 2006). 종다양도지수는 군락 Ⅰ이 세 군락 중 가장 낮게 나왔는데 이는 군락 내 조사구의 해발고가 다른 두 군락에 비해 높기 때문인 것으로 판단되었다 (Cho and Lee, 2011). 반면, 군락 Ⅲ이 0.8891로 세 군락 중 가장 높게 나왔다. 인위적인 관리가 이루어져 교목층은 곰솔 단일수종, 아교목층과 관목층 출현 종수가 다른 두 군락에 비해 낮은 편임에도 불구하고 교목층을 제외한 아교목층과 관목층에 서 특정 종의 우점도가 낮기 때문으로 사료되었다.

각 군락별 개체수의 분포정도를 나타내는 균재도는 1에 가 까울수록 종별 개체수가 균일함을 뜻하는데(Song et al., 2009), 군락 Ⅲ이 0.6441로 비교적 가장 균일하였다. 반면 군 락 Ⅰ은 0.1670으로 가장 낮은데 이는 마삭줄에 의한 결과로 판단되었다.

우점도는 그 값이 0.9 이상일 경우 한 종이 강하게 우세하는 것이며, 0.3~0.7 사이는 한 종이 약하게 우점하거나 두 종으로 나뉘어 우점하고, 0.1~0.3에서는 다수의 중요종에 의해 우점도 가 나누어지는 것을 의미한다(Whittaker, 1965). 각 군락별 우점도는 군락 Ⅰ이 0.8330으로 가장 높은데 교목층에서 졸참 나무와 소나무가 함께 우점하였으며, 군락 Ⅲ은 교목층에서 곰솔 단일수종이지만 아교목층과 관목층에서 다양한 상록활엽 수종이 출현하여 0.3559로 가장 낮게 나타났다.

부산광역시를 포함하는 경상권역의 전체 종다양도지수는 1.465, 균재도 0.811, 우점도 0.189인데(Park et al., 2013), 이에 비해 군락 Ⅰ, Ⅱ, Ⅲ의 종다양도지수, 균재도는 낮고 곰솔, 사스레피나무, 동백나무 등 특정 종의 우점도가 높은 것 을 알 수 있었다. 또 다른 부산광역시를 대상으로 한 본 연구와 유사한 식생구조를 가진 Kim(2012)의 연구에서는 교목층에서 곰솔과 졸참나무가 경쟁 상태이며 아교목층에서 사스레피나무 가 우점하는 군락의 종다양도는 0.8250~1.0042인데 반해 군락 Ⅰ (졸참나무-소나무-사스레피나무 군락)은 0.3141로 낮았다. 또한 교목층에 곰솔, 아교목층에 사스레피나무가 우점하는 군 락의 유사도 또한 Kim(2012)의 연구에서는 0.9657~1.3986 이었으나 군락 Ⅱ (곰솔-사스레피나무 군락)은 0.6374로 비교적 낮은 편에 속했다.

5. 흉고직경급별 분석

군락별 주요 출현종을 중심으로 수령 및 임분동태를 간접적 으로 알 수 있는 흉고직경급 분석을 실시하였다(Table 6). 본 연구에서는 국가산림자원조사 현지조사지침서의 내용을 참고 하여 흉고직경 6cm미만의 수종을 치수, 6~18cm를 소경목, 18~30cm를 중경목, 30cm 이상을 대경목으로 분류하였다.

군락 Ⅰ에서 졸참나무는 소경목에서 대경목에 해당하는 DBH 2cm 미만~52cm 구간에 고르게 나타나며 그 중 DBH 12~17cm에서 14개체, 관목층에 64개체가 출현하였으며, 굴 참나무는 DBH 2~47cm 사이 및 관목층에 분포하였다. 개서어 나무의 경우 DBH 2~52cm 이상의 대경목까지 고르게 출현하 였다. 소나무는 DBH 7~47cm에서 나타났으며 곰솔은 DBH 12~47cm에 나타났다. 사스레피나무는 DBH 2~32cm에 고르 게 분포하며, DBH 2~7cm 사이에 42개체가 나타났다. 졸참나 무와 굴참나무 등 참나무류가 관목층부터 대경목까지 분포한 반면 소나무와 곰솔은 어린 치수목이 없어 산림천이 양상에 따라 곰솔·소나무의 수세는 약해지며 참나무 세력이 확장될 것으로 판단되었다. 더 나아가 치수에서 대경목까지 극상수종 중 하나인 개서어나무가 고르게 출현하고 있어 낙엽활엽수림에 서 극상단계의 숲으로 발전할 것으로 예측되었다(Lee et al., 1996;Kim et al., 2011;Lee et al., 2013).

군락 Ⅱ에서는 곰솔이 DBH 7~52cm로 소경목부터 대경목 까지 고르게 분포하였다. 졸참나무는 DBH 2~42cm에 나타났 으며, DBH 2~7cm 구간에 40개체, 관목층에 316개체로 치수 와 관목층에 가장 많이 분포하였다. 사스레피나무의 경우 치수 부터 소경목에 해당하는 2cm 미만부터 17cm 사이에 분포하였 으며, DBH 2~7cm, 7~12cm 구간에서 각각 473, 133개체가 출현하였다. 병․충해, 자연재해 등의 교란요인이 없는 한 곰솔- 사스레피나무의 세력이 지속될 것으로 판단되었다.

군락 Ⅲ에서 곰솔 DBH 17~52cm 이상의 구간에 분포하였 으며, 동백나무와 사스레피나무는 소경목에서 중경목에 해당 하는 DBH 2~22cm 구간에 출현하였으며, 두 수종 모두 관목 층에 집중적으로 분포하였다. 소경목부터 중경목에 해당하는 구간에 동백나무는 210개체, 사스레피나무는 23개체가 출현하 였다. 앞서 언급한 군락 Ⅲ의 대상지는 동백섬으로 동백나무에 대한 집중적인 관리가 이루어지고 있어 사스레피나무보다 동백 나무가 훨씬 더 많이 출현하였다.

6. 종간상관관계

사스레피나무가 출현한 전체 89개의 조사구에서 126종이 출 현하였다. 이 중 출현빈도 각각 10%, 20% 이상인 수종과 조사 구 내 출현하는 126수종 전체를 대상으로 나누어서 SPSS 25 를 통해 상관분석을 실시하였다(Table 7).

전체 126수종 중 출현빈도 20% 이상인 수종을 대상으로 사 스레피나무와 상관관계 분석을 한 결과, 유의한 값을 가지는 수종은 때죽나무, 광나무, 마삭줄, 팔손이, 동백나무가 있었다. 5개의 수종 모두 유의확률 0.05 이하로 통계적으로 유의한 상 관관계이며, 사스레피나무와 모두 부(-)의 상관관계를 이루고 있었다.

전체 126수종 중 출현빈도 10% 이상을 기준으로 사스레피나 무와 상관관계 분석을 하였을 때, 출현빈도 20% 이상에서 사 스레피나무와 부(-)의 상관관계를 이루는 5개의 수종을 포함하 여 정(+)의 상관관계를 이루는 돈나무, 자금우가 추가 되었다.

126수종 전체에서 사스레피나무와 상관관계분석 결과, 유의 한 관계를 가지는 수종은 때죽나무, 광나무, 마삭줄, 팔손이, 동백나무, 돈나무, 자금우, 조릿대, 떡갈나무로 총 9 수종이 있었다. 이 중 사스레피나무와 정(+)의 상관관계를 갖는 수종 은 돈나무, 자금우, 조릿대, 떡갈나무가 있으며, 부(-)의 상관관 계를 갖는 수종은 때죽나무, 광나무, 마삭줄, 팔손이, 동백나무 가 있다. 사스레피나무와 부(-)의 상관관계를 갖는 수종은 전부 출현빈도가 20% 이상인 수종이었으며, 아교목층에 주로 분포 하여 사스레피나무와 동일한 생태적 지위를 가져 경쟁 상태에 있었다. 이 중 동백나무와의 상관계수가 -0.384**로 부(-)의 상관관계가 다른 4개의 수종에 비해 높은 것을 알 수 있었다. 이로 인해 군락 Ⅲ의 경우 인위적인 관리가 없을 경우 아교목 층에서 사스레피나무와 동백나무의 경쟁이 있을 것으로 예상 되었다.

7. 종합고찰

부산광역시의 사스레피나무가 출현하는 89개 조사구를 대 상으로 식생조사를 실시하고 분석하였다. TWINSPAN과 DCA 분석결과 3개의 군락으로 나누어졌다. 사스레피나무는 기후변화, 토양조건 등에 영향을 받으며 교목층 및 아교목층에 우점하고 있는 수종에 따라서도 차이를 보였다. 대부분의 조사 구는 해안가를 따라 위치하였으나 군락 Ⅰ은 졸참나무-소나무- 사스레피나무 군락으로 다른 군락에 비해 고도가 높은 내륙 곳곳에 많이 분포하였다. 이러한 특징으로 군락 Ⅰ내에는 양호 한 생태적 적지에서는 다른 수종과의 경쟁에서 밀리고 건조하 고 척박한 능선과 같은 곳에 생육하는 소나무가 많이 분포하였 다(Cho and Lee, 2011). 군락 Ⅰ의 평균상대우점치와 흉고직 경급별 분포를 보면 교목층에 졸참나무와 소나무가 경쟁하는 가운데 극상수종에 해당하는 개서어나무가 고르게 분포하고 있어 당분간은 현 상태를 유지할 것으로 예측되나 장기적으로 봤을 때 졸참나무, 개서어나무 등을 포함한 낙엽활엽수림으로 발전할 것으로 예측되었다. 사스레피나무는 졸참나무 및 개서 어나무와 유의한 관계가 없었지만 떡갈나무와는 0.222로 정 (+)의 상관관계를 가지고 있어 향후 천이가 진행되어도 현 세력 을 유지할 것으로 예상되었다. 군락 Ⅱ는 곰솔-사스레피나무 군락으로 간혹 내륙에 분포하였지만 해안가를 따라 분포하는 비율이 더 높았다. 교목층에 곰솔이 압도적으로 우점하고 있었 으며, 아교목층에는 사스레피나무가 우점하는 전형적인 곰솔- 사스레피나무 군락의 형태로 특별한 교란요인이 없는 한 현 상태가 유지 될 것으로 예측되었다. 군락 Ⅲ은 곰솔-동백나무 군락으로 조사구는 동백섬에 위치해 있어 동백나무를 중심으로 인위적인 관리가 이뤄지고 있었다. 관리가 이루어지는 곳은 동 백나무가 우점하고 있었으나 비관리지역은 사스레피나무가 우 점하고 있었다. 사스레피나무는 동일한 생태적 지위를 가진 동 백나무와 -0.384 수치로 부(-)의 상관관계를 가지는 것으로 나 타났다. 동백나무를 중심으로 한 지속적인 관리가 없을 경우 사스레피나무가 건조하고 척박한 토양 환경에 유리하여 우점할 것으로 예상되었다. 관목층에는 후박나무, 돈나무, 호랑가시나 무, 아왜나무 등 다수의 상록활엽수가 출현하고 있는데 그 중 후박나무는 상록활엽수림 식생천이 과정의 수종이며, 계속해 서 기온 상승이 이루어진다면 곰솔이 쇠퇴하고 상록활엽수림이 우점하는 식생으로 천이될 것이라 판단되었다(Lee et al., 2016). 또한 군락 Ⅲ은 관리로 인해 교목층은 단일수종, 아교 목층과 관목층도 두 군락에 비하면 종수가 낮음에도 불구하고 다른 두 군락에 비해 종다양도 및 균재도가 높았다. 이는 곰솔 이 단일수종으로 나타난 교목층을 제외하고 아교목층과 관목층 에 출현한 종수가 조사구 면적대비 균일하게 분포하며 관리가 미비한 아교목 및 관목층에서 팽나무, 뽕나무, 후박나무 등 다 양한 수종들이 출현하였기 때문으로 사료되었다.

기후변화 및 서식환경의 변화에 따라 사스레피나무는 해안 림뿐만 아니라 내륙으로 점차 북상할 것으로 예상되며(Yun et al., 2011;Yun et al., 2017), 교목층에 우점하는 종에 따라 서도 수세에 차이를 보이므로 이에 따라 국가 기후변화 생물지 표종인 사스레피나무를 중심으로 한 부산광역시 상록활엽수림 에 대한 지속적인 모니터링이 필요할 것으로 판단된다.