선택한 단어 수는 입니다.
최소 단어 이상 선택하여야 합니다.
최소 단어 이상 선택하여야 합니다.
선택한 단어 수는 30입니다.
최대 10 단어까지만 선택 가능합니다.
최대 10 단어까지만 선택 가능합니다.
문제 정의
- 최근 들어 기후변화가 산림과 임목에 미치는 영향과 그에 따른 변화를 구명하기 위하여 다양한 연구들이 진행되고 있지만 기후와 임목 간의 상관관계를 보다 정확히 분석할 수 있는 생태학적 측면에서의 연륜연대학 연구는 미흡한 실정이다. 본 연구는 한반도의 주요 산림자원이며 대표적인 한대 수종으로 기후변화에 따른 영향이 크게 나타날 것으로 예상되는 잣나무(Pinus koraiensis)를 대상으로 기후요소와 연륜생장과의 관계를 구명과 연륜기후학적 모델을 개발하고 RCP 기후변화 시나리오로 기후변화에 따른 잣나무의 연륜생장 변화를 예측하여 향후 지속될 기후변화에 따른 효과적인 대응책을 마련하기 위해 시행하였다.
제안 방법
- 잣나무 조사지역을 기후 특성을 바탕으로 상호 유사성이 높은 몇 개의 군집으로 분류하여 각 군집내의 동질성과 군집간의 이질성을 알아보기 위하여 군집분석(Cluster analysis)을 실시하였다. 권역분류는 각 조사지역의 월평균기온과 월강수량을 유클리디안 거리(Euclidean distance)로 환산하여 분류하였다. 본 연구에서는 미래의 연륜생장 변화를 예측하기 위하여 군집분석을 통해 분류한 기후군집별 연륜연대기를 이용하였다.
- 연륜기후학적 모델은 기후변수를 독립변수로 사용하여 단계적 회귀(Stepwise regression)분석으로 기후-연륜생장 관계식을 작성하였다(Fritts and Xiangding, 1986). 기후변수들 사이에 존재하는 다중공선성은 분산팽창계수(VIF)를 사용하여 제거하고 검증하였다. 생장지수 추정식은 추정식 작성에 사용되지 않은 기간의 연륜연대기를 이용하여 그 적합성을 검증하였다.
- 기후변화에 따른 잣나무 연륜생장 변화를 본 연구의 기후-연륜생장 관계식과 RCP 기후변화 시나리오를 이용하여 예측하였다. 연륜생장 변화 예측에 사용한 기후자료는 기상청 기후변화정보센터의 한반도 기후변화 시나리오 RCP 4.
- 기후와 연륜생장과의 상관 분석을 실시하여 기후변화에 따른 연륜생장 예측을 위하여 연륜기후학적 모델을 개발하였다. 연륜기후학적 모델은 기후변수를 독립변수로 사용하여 단계적 회귀(Stepwise regression)분석으로 기후-연륜생장 관계식을 작성하였다(Fritts and Xiangding, 1986).
- 본 연구는 잣나무를 대상으로 연륜생장과 기후와의 관계를 구명하기 위해 군집분석을 통해 기후 군집을 분류하였고 기후인자에 의해 연륜생장을 예측할 수 있는 생장모델을 개발하여 기후변화 시나리오 RCP 4.5와 RCP 8.5를 이용하여 연륜생장의 변화를 예측하였다. 연륜기후학적 모델은 잣나무 기후군집별 기후변수와 표준화된 잔차 연대기를 이용하여 개발하였다.
- 권역분류는 각 조사지역의 월평균기온과 월강수량을 유클리디안 거리(Euclidean distance)로 환산하여 분류하였다. 본 연구에서는 미래의 연륜생장 변화를 예측하기 위하여 군집분석을 통해 분류한 기후군집별 연륜연대기를 이용하였다.
- 기후변수들 사이에 존재하는 다중공선성은 분산팽창계수(VIF)를 사용하여 제거하고 검증하였다. 생장지수 추정식은 추정식 작성에 사용되지 않은 기간의 연륜연대기를 이용하여 그 적합성을 검증하였다. 적합성 검증을 위한 평가 통계량은 모형의 추정편의(model's estimation bias), 모형의 정도(model's precision), 그리고 측정치에 대한 오차의 평균평방화(mean square error type of measure)인 모형의 표준오차이다(Shin, 1990; Arabatzis and Burkhart, 1992; Shin et al.
- 5를 이용하여 연륜생장의 변화를 예측하였다. 연륜기후학적 모델은 잣나무 기후군집별 기후변수와 표준화된 잔차 연대기를 이용하여 개발하였다. 개발된 생장지수 추정식에는 2~4개의 기후변수가 사용되었고 결정계수는 0.
- 기후변화에 따른 잣나무 연륜생장 변화를 본 연구의 기후-연륜생장 관계식과 RCP 기후변화 시나리오를 이용하여 예측하였다. 연륜생장 변화 예측에 사용한 기후자료는 기상청 기후변화정보센터의 한반도 기후변화 시나리오 RCP 4.5와 RCP 8.5에서 2011년부터 2100년까지의 기후변화 시나리오 별 월평균기온과 월강수량을 군집분석으로 분류한 기후군집에 맞추어 산출하였고, 산출한 기온과 강수량 추정치를 연륜기후학적 모델 개발을 통해 얻어진 잣나무의 기후군집 별 생장 추정식에 적용하여 연륜생장 변화를 예측하였다.
- 49의 범위에 있는 것으로 분석되었다. 잣나무 기후군집별 생장예측은 두 개의 RCP 기후변화 시나리오를 이용해 예측되었다. 해발고도가 높은 곳에 위치한 군집 2와 군집 3에서는 연륜 생장이 증가하였지만 대조적으로 해발고도가 낮은 지역에 위치한 군집 1에서는 연륜생장이 감소할 것으로 예측되었다.
- 2m 높이에서 채취하였다(Fritts, 1976). 채취된 연륜은 건조시킨 후 연마하여 연륜측정분석기(Digital Tree System DTRS-2000)를 이용하여 1/100mm 단위까지 측정하였다. 측정된 연륜은 크로스데이팅(Crossdating)을 실시하여 각 연륜에 정확한 연도를 부여하였다(Fritts and Swetnam, 1989).
- 채취된 연륜은 건조시킨 후 연마하여 연륜측정분석기(Digital Tree System DTRS-2000)를 이용하여 1/100mm 단위까지 측정하였다. 측정된 연륜은 크로스데이팅(Crossdating)을 실시하여 각 연륜에 정확한 연도를 부여하였다(Fritts and Swetnam, 1989). 크로스데이팅은 COFECHA 프로그램(Holmes, 1983)을 이용하여 통계적으로 재확인하였다.
대상 데이터
- 연륜연대학 연구에 있어서 기후자료는 개별 측후소의 자료보다 측후소가 없는 지역까지 연구범위를 확장할 수 있고 신뢰성이 높은 지역평균값을 사용한다(Fritts, 1976). 본 연구는 당년의 원기가 형성되는 전년도 8월부터 금년도 9월까지 14개월의 월평균기온, 월강수량과 이 기간을 생장기와 비생장기로 구분한 34개 기후변수의 시군 단위 지역평균값을 사용하였다(Table 2). 시군 단위의 월평균기온과 월강수량은 기상청에서 래스터형 데이터로 제공한 전국의 일별 기온과 강수량을 GIS프로그램 ArcGIS 10.
- 연구대상지는 제4차 국가산림자원조사의 잣나무 분포지역 중 잣나무의 임목축적량이 많고 분포면적이 넓은 지역을 대상으로 생장상태를 고려하여 잣나무가 우세목과 준우세목으로 분포하는 곳에서 선정하였다. 조사지역은 강원도 춘천, 홍천, 설악산, 오대산, 경기도 가평, 광릉, 전라북도 남원, 경상남도 산청 등 8개 지역으로 조사지 개황은 Table 1과 같다.
- 연륜시료는 압축이상재에 의한 변이를 제거하기 위하여 경사와 수직방향으로 지상으로부터 1.2m 높이에서 채취하였다(Fritts, 1976). 채취된 연륜은 건조시킨 후 연마하여 연륜측정분석기(Digital Tree System DTRS-2000)를 이용하여 1/100mm 단위까지 측정하였다.
- 연구대상지는 제4차 국가산림자원조사의 잣나무 분포지역 중 잣나무의 임목축적량이 많고 분포면적이 넓은 지역을 대상으로 생장상태를 고려하여 잣나무가 우세목과 준우세목으로 분포하는 곳에서 선정하였다. 조사지역은 강원도 춘천, 홍천, 설악산, 오대산, 경기도 가평, 광릉, 전라북도 남원, 경상남도 산청 등 8개 지역으로 조사지 개황은 Table 1과 같다. 표준지는 37개소를 조사하였고 채취한 연륜시료 중 이용 가능한 283개를 분석하였다.
- 조사지역은 강원도 춘천, 홍천, 설악산, 오대산, 경기도 가평, 광릉, 전라북도 남원, 경상남도 산청 등 8개 지역으로 조사지 개황은 Table 1과 같다. 표준지는 37개소를 조사하였고 채취한 연륜시료 중 이용 가능한 283개를 분석하였다.
데이터처리
- 본 연구는 당년의 원기가 형성되는 전년도 8월부터 금년도 9월까지 14개월의 월평균기온, 월강수량과 이 기간을 생장기와 비생장기로 구분한 34개 기후변수의 시군 단위 지역평균값을 사용하였다(Table 2). 시군 단위의 월평균기온과 월강수량은 기상청에서 래스터형 데이터로 제공한 전국의 일별 기온과 강수량을 GIS프로그램 ArcGIS 10.1과 통계프로그램 R 3.0.1을 이용하여 계산하였다.
- 기후와 연륜생장과의 상관 분석을 실시하여 기후변화에 따른 연륜생장 예측을 위하여 연륜기후학적 모델을 개발하였다. 연륜기후학적 모델은 기후변수를 독립변수로 사용하여 단계적 회귀(Stepwise regression)분석으로 기후-연륜생장 관계식을 작성하였다(Fritts and Xiangding, 1986). 기후변수들 사이에 존재하는 다중공선성은 분산팽창계수(VIF)를 사용하여 제거하고 검증하였다.
- 표준화된 연대기는 자기 상관을 제거하기 위하여 ARSTAN 프로그램(Cook, 1985)을 사용하여 자기상관을 제거한 잔차연대기(residual chronology)를 작성하였다(Cook and Holmes, 1986; Speer, 2010). 연륜연대기의 기초통계자료는 평균, 연속된 두 연륜폭간 변이 크기를 나타내는 평균민감도(Mean Sensitivity), 전년도의 환경이 얼마나 수목의 생장에 영향을 주는가를 나타내는 자기상관계수(Auto-correlation), 임목내와 임목간의 차이를 나타내는 상호상관계수(Correlation coefficient), 연대기의 기간별 신호강도를 나타내는 신호대 잡음비(Signal to Noise Ratio), 모집단에 대한 지역연대기의 신호강도를 나타내는 모수신호설명(Expressed Population Signal)을 산출하였다(Fritts, 1976; Wigley et al., 1984; Cook and Kairiukstis, 1990; Seo, 1999).
- 잣나무 조사지역을 기후 특성을 바탕으로 상호 유사성이 높은 몇 개의 군집으로 분류하여 각 군집내의 동질성과 군집간의 이질성을 알아보기 위하여 군집분석(Cluster analysis)을 실시하였다. 권역분류는 각 조사지역의 월평균기온과 월강수량을 유클리디안 거리(Euclidean distance)로 환산하여 분류하였다.
- 측정된 연륜은 크로스데이팅(Crossdating)을 실시하여 각 연륜에 정확한 연도를 부여하였다(Fritts and Swetnam, 1989). 크로스데이팅은 COFECHA 프로그램(Holmes, 1983)을 이용하여 통계적으로 재확인하였다. 측정된 연륜은 수령에 따른 생물학적 생장추세, 경쟁과 교란에 의한 변동추세와 비기후인자에 의한 저주파변동을 제거하기 위하여 이중표준화법을 사용하여 표준화를 실시하였다.
이론/모형
- 크로스데이팅은 COFECHA 프로그램(Holmes, 1983)을 이용하여 통계적으로 재확인하였다. 측정된 연륜은 수령에 따른 생물학적 생장추세, 경쟁과 교란에 의한 변동추세와 비기후인자에 의한 저주파변동을 제거하기 위하여 이중표준화법을 사용하여 표준화를 실시하였다. 1차로 음지수곡선으로 수령에 따른 생물학적 생장추세를 제거하고, 2차로 스플라인곡선(50%, 반응주기: 60년)으로 임분 경쟁과 교란에 따른 변동추세를 재표준화하였다(Fritts, 1976; Cook and Kairiukstis, 1990; Seo, 1999).
- 1차로 음지수곡선으로 수령에 따른 생물학적 생장추세를 제거하고, 2차로 스플라인곡선(50%, 반응주기: 60년)으로 임분 경쟁과 교란에 따른 변동추세를 재표준화하였다(Fritts, 1976; Cook and Kairiukstis, 1990; Seo, 1999). 표준화된 연대기는 자기 상관을 제거하기 위하여 ARSTAN 프로그램(Cook, 1985)을 사용하여 자기상관을 제거한 잔차연대기(residual chronology)를 작성하였다(Cook and Holmes, 1986; Speer, 2010). 연륜연대기의 기초통계자료는 평균, 연속된 두 연륜폭간 변이 크기를 나타내는 평균민감도(Mean Sensitivity), 전년도의 환경이 얼마나 수목의 생장에 영향을 주는가를 나타내는 자기상관계수(Auto-correlation), 임목내와 임목간의 차이를 나타내는 상호상관계수(Correlation coefficient), 연대기의 기간별 신호강도를 나타내는 신호대 잡음비(Signal to Noise Ratio), 모집단에 대한 지역연대기의 신호강도를 나타내는 모수신호설명(Expressed Population Signal)을 산출하였다(Fritts, 1976; Wigley et al.
성능/효과
- 802로 가장 낮은 신호강도를 보였다. EPS는 잣나무 조사 대상지 지역 전체 평균이 0.867로 연대기 분석에 이용 가능한 것으로 나타났다(Table 4).
- Table 6은 잣나무 군집별 생장지수 추정식의 적합성 검정 결과로 각 추정식의 평균편의(MD)는 -0.0300~0.0315으로 상당히 낮은 평균편의를 보이고 모형의 정도(SDD)와 모형의 표준오차(SED)도 비교적 작은 값을 보여 본 연구에서 개발한 생장지수 추정식을 사용하는데 문제가 없을 것으로 판단된다.
- 연륜기후학적 모델은 잣나무 기후군집별 기후변수와 표준화된 잔차 연대기를 이용하여 개발하였다. 개발된 생장지수 추정식에는 2~4개의 기후변수가 사용되었고 결정계수는 0.35~0.49의 범위에 있는 것으로 분석되었다. 잣나무 기후군집별 생장예측은 두 개의 RCP 기후변화 시나리오를 이용해 예측되었다.
- 군집 2와 군집 3에서는 기후변화에 따라 생장이 2100년까지 지속적으로 증가하는 것으로 예측되었다. 군집 2에서는 RCP 4.5에서 7.9%, RCP 8.5에서 13.1% 생장률이 증가하였다. 군집 3은 RCP 4.
- 1% 생장률이 증가하였다. 군집 3은 RCP 4.5에서 14.7% 증가하였고 RCP 8.5에서 2100년에 생장률이 22.4% 증가하여 최대생장률 증가를 보였다(Table 6). 기후군집별 잣나무 연륜생장 예측결과 기후변화가 진행됨에 따라 군집 1은 연륜생장이 감소하는 경향을 보였으나 군집 2와 군집 3은 연륜생장이 증가하는 것으로 예측되었다.
- 4% 증가하여 최대생장률 증가를 보였다(Table 6). 기후군집별 잣나무 연륜생장 예측결과 기후변화가 진행됨에 따라 군집 1은 연륜생장이 감소하는 경향을 보였으나 군집 2와 군집 3은 연륜생장이 증가하는 것으로 예측되었다. 기후변화 시나리오에 따른 연륜생장 변화도 생장과 감소에서 모두 RCP 4.
- 기후군집별 잣나무 연륜생장 예측결과 기후변화가 진행됨에 따라 군집 1은 연륜생장이 감소하는 경향을 보였으나 군집 2와 군집 3은 연륜생장이 증가하는 것으로 예측되었다. 기후변화 시나리오에 따른 연륜생장 변화도 생장과 감소에서 모두 RCP 4.5에 비해 RCP 8.5에서 큰 변화를 보였다.
- 9oC로 나타났다. 기후변화에 따라 모든 군락에서 생장기 4월과 5월의 평균 기온이 상승하였지만 연륜생장이 증가하는 것으로 예측된 군락 2와 군락 3은 현재의 평균기온과 2100년까지의 기후 추정치 모두 연륜생장이 감소하는 것으로 나타난 군락 1에 비하여 낮게 나타났다. 월강수량은 기온과 다르게 경년ㆍ공간 변동성이 크고 상대적으로 증가 추세는 약하게 나타나 기온과 비교해 연륜생장에 미치는 영향이 작은 것으로 사료된다.
- 해발고도가 높은 곳에 위치한 군집 2와 군집 3에서는 연륜 생장이 증가하였지만 대조적으로 해발고도가 낮은 지역에 위치한 군집 1에서는 연륜생장이 감소할 것으로 예측되었다. 따라서 잣나무는 한대수종이지만 일정범위내의 기온 상승 조건에서는 연륜생장이 증가하는 것으로 평가되었다.
- , 2001). 이러한 결과들로 보아 수목의 생장에는 다양한 요소들이 연관되지만 해발고도가 높은 지역에서 기후변화에 의한 일정범위의 기온상승은 한대수종의 연륜생장을 촉진시킬 수 있는 것으로 판단된다. 잣나무는 현재 분포지역이 넓고 대부분의 지역에서 양호한 생장상태를 유지하고 있지만 기후변화가 진행됨에 따라 상당한 생장변화가 뒤따를 것으로 예측되었다.
- 잣나무 조사지 8개 지역의 기후 특성을 알아보기 위해 유클리디안 거리지수를 이용한 군집분석을 실시한 결과 3개의 기후군집으로 분류되었다(Fig. 1). 군집 1은 춘천, 홍천, 가평, 그리고 광릉 등 중부내륙지역, 군집 2는 오대산과 설악산의 백두대간 산악지역, 그리고 마지막으로 군집 3은 남원과 산청의 남부내륙지역으로 분류되었다.
- 잣나무의 3개 군집을 대상으로 개발한 생장지수 추정식에는 2~4개의 월평균기온 및 월강수량 독립변수가 사용되었으며, 결정계수는 0.35~0.49로 나타났다. Laroque and Smith(2005)가 캐나다 뱅쿠버섬에서 Tsugamertensiana, Chamaecyparis nootkatensis, Thuja plicata, Pseudotsuga menziesii, 그리고 Tsuga heterophylla를 대상으로 한 연륜생장 모형의 결정계수 0.
- 기후변화로 한반도의 기온이 상승하면 고산지대 및 한랭한 기후를 선호하는 한대수종인 잣나무는 생육지의 기온 상승에 따른 온도 스트레스로 부정적인 영향을 받을 것으로 예상되고 있다. 하지만 본 연구에서는 군집 1은 기후변화가 진행됨에 따라 연륜생장이 지속적으로 감소하여 부정적인 생장 경향을 보였지만 군집 2와 군집 3에서는 2100년까지 현재보다 연륜생장이 증가하는 것으로 예측되어 잣나무 분포지역에 따라 상반된 결과가 나타났다. 기후군집별 잣나무 생장지수 추정식에 사용된 기후변수의 기후값은 Table 8과 같다.
- 잣나무 기후군집별 생장예측은 두 개의 RCP 기후변화 시나리오를 이용해 예측되었다. 해발고도가 높은 곳에 위치한 군집 2와 군집 3에서는 연륜 생장이 증가하였지만 대조적으로 해발고도가 낮은 지역에 위치한 군집 1에서는 연륜생장이 감소할 것으로 예측되었다. 따라서 잣나무는 한대수종이지만 일정범위내의 기온 상승 조건에서는 연륜생장이 증가하는 것으로 평가되었다.
후속연구
- 잣나무는 현재 분포지역이 넓고 대부분의 지역에서 양호한 생장상태를 유지하고 있지만 기후변화가 진행됨에 따라 상당한 생장변화가 뒤따를 것으로 예측되었다. 따라서 추후 기후변화에 따른 보다 구체적인 대응책 마련을 위해서는 연륜기후학과 연륜생태학적 연구를 통한 수종별, 지역별 기후와 연륜 생장과의 관계 구명이 필요하다고 사료된다. 또한 기후변화의 영향이 가시화되고 있는 현실에서 임목의 생장예측은 산림경영과 시업에 있어 의사결정을 위해 필수적으로 요구되는 기본자료 중 하나로 그 활용가치가 매우 높다.
- 하지만 수많은 요소들이 상호작용한 결과인 기후나 생장과 같은 현상은 하나의 수식으로 완벽하게 표현되지 않아 미래를 항상 예측할 수는 없다는 어려움을 가지고 있어 정기적이고 장기적인 모니터링과 연구가 지속되어야 할 것으로 사료된다. 이와 같은 이유로 본 연구는 기후변화에 대응하기 위한 산림경영 목표의 수립과 관리 방법의 개선에 필요한 정보를 제공할 것으로 판단된다.
- 또한 기후변화의 영향이 가시화되고 있는 현실에서 임목의 생장예측은 산림경영과 시업에 있어 의사결정을 위해 필수적으로 요구되는 기본자료 중 하나로 그 활용가치가 매우 높다. 하지만 수많은 요소들이 상호작용한 결과인 기후나 생장과 같은 현상은 하나의 수식으로 완벽하게 표현되지 않아 미래를 항상 예측할 수는 없다는 어려움을 가지고 있어 정기적이고 장기적인 모니터링과 연구가 지속되어야 할 것으로 사료된다. 이와 같은 이유로 본 연구는 기후변화에 대응하기 위한 산림경영 목표의 수립과 관리 방법의 개선에 필요한 정보를 제공할 것으로 판단된다.
※ AI-Helper는 부적절한 답변을 할 수 있습니다.