생물다양성협약(CBD)에서는 전 지구적으로 발생하는 환경문제의 대응방안으로 육상의 17%를 보호지역으로 지정 하도록 할 것을 권고하고 있다. 우리나라도 국제적 수준에 맞춰 보호지역을 지정하고 보호지역 지정이 가지는 의미를 설명할 필요성이 생겼으며, 이에 따라 생태계서비스에 관한 연구의 필요성이 요구되고 있다. 우리나라는 보호지역 지정을 자연공원법에 의해 위계별로 국립공원, 도립공원, 군립공원으로 지정하고 있다. 그러나 정치, 행정적 측면이 우선시되면서 생태계서비스 가치평가 연구 및 서식지 관리가 국립공원에 편중되어 도립 및 군립공원은 다소 미비하여 연구의 필요성이 있다. 본 연구에서는 국립공원에 비해 서식지 관리 및 생태계서비스 가치평가 연구가 미비한 가지산도 립공원을 대상으로 InVEST 모델 중 InVEST Habitat Quality 모델을 사용하여 서식지 질을 평가하고, 분석 결과를 16개의 산악형 국립공원과 비교하였다. 분석 결과 가지산도립공원의 서식지질 값은 0.83이었으며, 주변 지역에 비해 서식질 질이 높게 나타났다. 3개 지구별 서식지질 차이를 분석한 결과 통도사지구와 내원사지구가 0.84, 석남사지구가 0.83으로 나타났고 용도지구별로는 자연보존지구, 자연환경지구, 문화유산지구, 공원마을지구 순으로 서식지질 값이 높게 나타났다. 기존 국립공원 서식지질 분석 결과와 비교한 결과 가지산도립공원은 무등산국립공원 수준의 자연성을 나타내었다. 이러한 결과는 추후 도립공원의 생물다양성 보전 및 생태계서비스 증진을 위한 정책 수립과 관리방안 수립의 객관적인 자료로 활용될 수 있을 것이다.

The purpose of this study is to investigate a vegetation structure around Tongdosa temple forest and provincial park and to provide preliminary data. In order to look over the vegetation status, an actual vegetation map was drawn around study area. Vegetation structure survey was carried out for 6 representative communities of actual vegetation which were Quercus variavilis community, Carpinus tschonoskii community, Pinus densiflora community, P. densiflora-Broadleaf deciduous Forest community, Q. mongolica community and Broadleaf deciduous Forest community. The area of the Tongdosa district measured 29,202,262 ㎡. Actual vegetation type were divided into 35 types, and the ratio of Q. variavilis community was 32.35 % (9,447,932 ㎡). To investigate the structure of 6 representative communities, 58 plots were set up and unit area plots measured 100 ㎡. The estimated age of the forest is 50~100-years-old and the oldest tree P. densiflora is 113-years-old.

본 연구는 우리나라 3대 사찰 중 하나인 통도사를 중심으로 지정되어 있는 가지산도립공원 통도사지구의 식물군락과 환경요인과의 상관관계를 알아보고자 수행하였다. 식물군락구조를 알아보기 위해 200∼300m2 크기의 조사구 31개소를 설치하여 조사 분석하였으며, 환경요인으로는 해발고도, 경사도 등의 물리적 환경과 토양산도, 유기물함량, 전질소함량, 치환성양이온 등의 토양특성을 조사 분석하였다. 환경요인과 식물군락과의 상관관계분석을 위해 DCCA기법을 활용하였다. 식물군락구조를 분석한 결과, 소나무-갈참나무군락, 개서어나무-참나무류군락, 소나무-개서어나무군락, 소나무-굴참나무군락, 굴참나무-신갈나무군락으로 총 5개 군락으로 분리되었다. 식생천이는 소나무에서 굴참나무와 신갈나무 등 참나무류와 개서어나무로 진행되고 있었다. DCCA분석 결과, 가지산도립공원 통도사지구 내 군락분포는 해발고도에 가장 큰 영향을 받았으며, 해발고도가 비교적 낮은 곳에는 소나무-갈참나무군락, 높은 곳에는 굴참나무-신갈나무군락이 분포하였다.

가지산도립공원은 한반도 남동쪽에 위치하고 있으며, 국 민휴양 및 정서 함양에 이바지하기 위한 목적으로 1979년 11월에 지정되었다. 가지산도립공원은 백두대간에서 분기 되어 남동쪽으로 뻗어있는 낙동정맥에 일부 속해 있으며, 주위의 운문산(1,190m), 천황산(1,189m), 신불산(1,209m), 간월산(1,083m), 영취산(1,059m)과 함께 영남알프스로 불 린다(Kim et al., 2010). 또한 가지산도립공원은 사찰을 중 심으로 통도사지구, 내원사지구, 석남사지구 등 3개의 지구 로 나눠 관리중이다. 가지산도립공원에 관한 연구는 석남사지구, 내원사지구 등 지구별로 수행되었지만 통도시지구에 대한 식물사회학 적 연구는 부족한 실정이다. 최근 들어, 통도사지구의 자연 환경에 대한 전반적인 조사가 이루어졌으나(Kim et al., 2014) 현존식생과 그 식생구조에 대한 현황에 그치고 있다. 현존하는 식물사회는 그들이 처한 다양한 환경 인자에 적응 하여 이룩한 결과물이라 할 수 있으므로(Ellenberg, 1992), 본 연구는 가지산도립공원 통도사지구를 중심으로 식생구 조를 파악하고 대표적 식물군락과 그 환경요인의 상관관계 를 밝혀보고자 하였다. 2013년 7월 10일부터 12일까지 3일간 현장조사를 실시 하였고, 식물군락구조 특성을 파악하기 위해 금수암, 백운 암, 반야암 등을 중심으로 Ellenberg(1956)와 Westhoff and Maarel(1973)을 참고하여 200~300㎡ 크기의 방형구 (quadrat) 31개소를 설치하였다. 식생조사는 교목층, 아교 목층, 관목층으로 나누어 수관층위별로 실시하였고, 식생조 사 자료를 토대로 상대우점치(Brower and Zar, 1977)를 분 석하였다. 이를 토대로 TWINSPAN에 의한 classification 분석(Hill, 1979a)과 DCA ordination(Hill, 1979b)분석을 실시하였으며, 종 구성의 다양한 정도를 나타내는 척도인 종다양도 분석 및 유사도를 비교 분석하였다. 가지산도립공 원 통도사지구 산림지역의 토양 이화학적 특성을 알아보기 위해 31개 조사구 내에서 임의의 지점 3곳을 골라 Ao층을 걷어내고 표층으로부터 30cm정도의 깊이의 토양을 채취․혼 합하여 실내에서 음건하였다. 분석용 토양을 토양 및 식물 체 분석법(National Institute of Agricultural Science and Technology, 2000)에 의한 pH, EC, 유기물함량, 유효인산 과 전질소함량, 치환성양이온의 함량 등을 분석하였다. 산림군락은 환경요인에 따라 그 구조가 달라지며, 군락 구조를 분석하고 환경요인과의 관계를 분석하는데 있어 ordination방법이 널리 이용되고 있다(Song and Kim, 1992). Ordination 분석은 DCA(Detrended Correspondence Analysis) 의 확장인 DCCA(Detrended Canonical Corresponded Analysis)를 사용하였으며(Hill and Gauch, 1980), DCCA 결과에서 식물군락은 점으로 분포를 나타내고 환경요인은 화살표로 나타내어 화살표의 방향과 길이로부터 식물군락 과의 유의성을 알 수 있다(Jongman et al., 1995). MVSP (Multi-Variate Statistical Package) 3.2 program을 사용하 여 분석하였으며, ordination 분석을 통하여 산림식생과 환 경 요인과의 관계를 규명하고자 하였다. 전체 31개 조사구에 대해 유형별 분류를 하기 위하여 classification 분석 중 지표종(indicator species)을 중심으로 군락을 분류할 수 있는 TWINSPAN기법을 적용한 결과, 소나무-갈참나무군락(Ⅰ), 개서어나무-참나무류군락(Ⅱ), 소나무-개서어나무군락(Ⅲ), 소나무-굴참나무군락(Ⅳ), 굴 참나무-신갈나무군락(Ⅴ) 등 총 5개 군락으로 분리되었다. Classification 분석의 상호보완적 방법으로 군락분포를 알 아보기 위해(Lee et al., 1994; Choi and Kang, 2006) ordination 분석 방법 중 DCA 기법을 적용하여 군락의 분포상 태를 살펴 본 결과, 소나무와 갈참나무가 우점하는 군락Ⅰ 과 굴참나무와 신갈나무가 우점하고 있는 군락Ⅴ를 제외하 면, 그 외 군락들은 소나무, 개서어나무, 굴참나무 등을 중심 으로 연속적으로 분포하였다.Classification 분석과 ordination 분석을 통해 분리된 5개 군락에 대해 각 군락별로 층위별 상대우점치 및 평균상대우 점치를 분석하였다. 군락Ⅰ은 교목층에서 소나무가 우점하 는 가운데 갈참나무, 졸참나무 등과 함께 경쟁하고 있었으 며, 아교목층에서는 소나무가 출현하지 않았다. 갈참나무, 노각나무, 비목나무 등이 아교목층에서 우점하고 있었고, 관목층에서는 노각나무, 비목나무 등이 우점하였다. 군락Ⅱ 는 교목층에서 개서어나무가 졸참나무, 굴참나무 등과 경쟁 하고 있었으나, 아교목층에서 개서어나무에 비해 졸참나무, 굴참나무의 세력이 크게 형성되지 못하였다. 관목층에서는 비목나무, 생강나무 등이 우점하고 있어, 군락Ⅱ는 향후에 도 개서어나무를 중심으로 한 낙엽활엽수림으로 유지될 것 으로 사료된다. 군락Ⅲ은 교목층에서 소나무와 개서어나무, 졸참나무 등이 경쟁을 하고 있었다. 아교목층에서 개서어나 무와 때죽나무가 우점하고 있어, 소나무와의 경쟁을 통해 장기적으로는 개서어나무를 중심으로 한 낙엽활엽수림으 로의 천이가 예상되었다. 군락Ⅳ는 소나무-굴참나무군락으 로 소나무와 굴참나무, 졸참나무 등의 참나무류가 경쟁하고 있어, 낙엽성 참나무류로의 천이가 예상되었다. 군락Ⅴ는 굴참나무와 신갈나무가 교목층에서 우점하고 있으며, 아교 목층에서는 쇠물푸레나무, 철쭉, 생강나무 등이 우점하였 다. 가지산도립공원 통도사지구의 31개 조사구에 대하여 토 양을 채취‧분석한 결과, 평균 수소이온농도는 pH 4.9로 경 남의 전체적인 평균 pH 5.27(Jeong et al., 2002)보다 낮게 나타났다. 수소이온농도는 여러 가지 요인에 의해 시‧공간 적으로 변화를 보이기 때문에 직접적인 비교는 어려우나 산림지역에서 pH 5.0 이하의 강산성토양 분포비율의 증가 는 환경오염물질에 의한 토양산성화물질의 기여도가 높을 가능성을 시사하고 있다(Lee and Park, 2001). 또한 침엽수 의 생육범위가 pH 4.8~5.5(Lee, 2000)인 점을 고려하고 가 지산도립공원 통도사지구에는 소나무림이 폭넓게 분포하 며 본 연구의 평균 수소이온농도 pH 4.9 이하의 조사구에서 소나무가 상당수 출현하는 점을 감안할 때, 그에 따른 영향 으로 판단된다. 유기물 함량은 약 96g/kg으로 나타났고, 전 질소함량은 0.33~0.74%의 범위로 우리나라 산림토양의 평 균값인 0.19%보다 다소 높게 나타났다. 31개 조사구의 산림식생을 5개 군락으로 분류하여 해발 고, 경사, 토양의 화학적 특성 등 11개의 환경요인과 DCCA ordination 분석한 결과, 해발고도, 경사도, 수소이온농도 (pH), 유기물함량, Ca2+, K+ 등이 식생분포에 영향을 미치는 환경인자로 나타났다. 제 1축에서 해발고 외에 수소이온농 도(pH), 유기물함량, 경사도 등이 비교적 높은 상관관계를 보였으며, 제 2축에서는 상관관계가 뚜렷하게 나타나지 않 았다. 해발고도는 식생분포에 영향을 미치는 환경인자들 중에 서 가장 중요한 인자로 알려져 있는데(Seo et al., 1995; Chung et al., 1997; Song et al., 1998; Lee et al., 2010), 본 연구에서도 해발고도가 군락의 분포에 상관이 높은 것으 로 나타났고, 굴참나무-신갈나무군락(Ⅴ)은 비교적 해발고 가 높은 지역에서 출현하였다.

경상남도 밀양 소재 가지산도립공원에서 채집한 317개체의 거미류는 20과 59속 66종으로 동정되었다. 종지수(species richness)는 깡충거미과 12종(20.0%), 꼬마거미과 8종(12.1%), 접시거미과 5종(7.6%), 왕거미과, 코리나거미과, 수리거미과 및 게거미과가 각 4종(6.1%), 갈거미과, 비탈거미과 및 새우게거미과가 각 3종(4.5%)의 순이었다. 동물지리구에 따른 종은 범열대종 1종(1.5%), 구북구종 4종(6.1%) 및

낙동정맥의 가지산～능동산 구간 능선부 주변의 식생구조를 파악하기 위하여 10개 조사지에 50개 조사구를 설치하여 조사를 실시하였다. Ordination기법 중 DCA 분석을 실시한 결과 군락의 종조성이 유사하여 연속적으로 나타났다. 각 조사지별로 식생구조를 분석한 결과 능선부의 주요 식생은 교목층에서 떡갈나무, 소나무 신갈나무이고 아교목층에서는 철쭉꽃, 관목층에서는 미역줄나무였다. 단위면적당(500㎥) 종다양도는 1.5315～2.4005였으며, 연륜분석결과 약 30～40년생 숲으로 밝혀졌다.