1. 서 론
2. 연구 방법
2.1 연구 대상 지역
2.2 다중분광광학센서 자료
2.3 이산화탄소 플럭스 산정 방법
2.4 다중분광광학센서를 활용한 토양수분 산정 방법
2.5 농업적 가뭄 지수
3. 연구 결과
3.1 광학센서 산출물의 비교 분석
3.2 토양수분-식생지수-이산화탄소 플럭스 간 상관 분석
4. 결 론
1. 서 론
기후변화의 영향으로 인해 전 세계에 가뭄 발생이 급증하고 있다. 가뭄은 학술적인 정의에 따라 다양하게 분류 되며, 농업적 가뭄은 단기간의 강수부족으로 인해 토양수분이 감소하여 식생 활력에 영향을 준다(Wilhite et al., 2007; Sivakumar et al., 2011; Kim et al., 2017; Sur et al., 2019). 농업적 가뭄을 분석하기 위하여 원격탐사 자료 기반의 토양수분이나 식생지수를 이용하여 가뭄지수를 산정하였으며, Kogan (1997)은 식생지수를 이용하여 식생 건강 지수(Vegetation Health Index, VHI)와 식생 조건 지수(Vegetation Condition Index, VCI)를 개발하였고, Keshavarz et al. (2014)은 토양수분을 이용하여 토양 습윤 부족 지수(Soil Wetness Deficit Index, SWDI)를 개발하였다. 원격탐사 기법은 미계측 지역을 포함한 넓은 지역을 모니터링 할 수 있다는 장점 때문에 가뭄 모니터링 관련 연구에서 활용되고 있다. 원격탐사자료는 위성에 결합되어 있는 센서들의 특징에 따라 다양한 시공간 해상도로 산출물이 제공되기 때문에 연구자들은 연구목적에 따라 산출물을 선택할 수 있다.
농업적 가뭄과 관련된 인자 중 식생지수의 경우, 광학 센서에서 제공되는 산출물을 활용한다. MODerate resolution Imaging Spectroradiometer (MODIS)에서 제공되는 정규 식생 지수(Normalized Difference Vegetation Index, NDVI), 향상 식생 지수(Enhanced Vegetation Index, EVI) 등이 연구에서 활용되었다(Peters et al., 2002; Hu et al., 2019). 토양수분의 경우, passive microwave sensor에서 제공되는 산출물을 활용한다. 마이크로웨이브를 통한 토양수분 산정은 지표면의 밝기 온도를 관측하여 관계식을 이용해 추산하는 방식으로 산출된다(Njoku et al., 2003). 지표면의 물리적인 특성에 따라 토양수분을 산정하기 때문에 값의 높은 정확성의 장점이 있지만 광학센서 산출물에 비해 공간 해상도가 낮다. 따라서, 대륙 단위 공간 범위와 같이 광활한 지역에 대해 분석하기에는 적합하나 유역이나 지역 단위 범위에서는 한계가 발생한다. 이러한 한계점을 극복하기 위해 광학 센서를 통한 토양수분 산정 방법에 대한 연구들이 소개되었다(Verstraeten et al., 2006; Chang et al., 2012). 광학센서를 통한 토양수분의 산정 기법은 지표면 온도와 지표 반사도(land surface albedo)의 관계식을 통해 산정하는 방법으로 입력 자료에 따라 공간 해상도 0.5 km의 토양수분을 산정할 수 있다.
농업적 가뭄 분석과 관련된 연구들은 토양수분과 식생지수를 통해 가뭄지수들을 개발하고, 가뭄 상태를 모니터링 및 분석하였고, 식생 활력과 광합성의 관계를 분석한 연구들이 많이 수행되었다(Running and Nemani, 1988; Muraoka et al., 2013). 하지만 토양수분, 식생지수와 광합성 사이의 관계에 대한 분석과 관련된 연구는 미흡하다. 농업적 가뭄이 장기간으로 지속된다면 식생 활력이 활발하지 못하기 때문에 광합성에 부정적인 영향을 받아 대기 중으로 배출되는 이산화탄소의 양이 증가한다. 본 연구에서는 가뭄과 같이 극한 기후 조건인 경우에 한정했을 때, 인자들의 관계를 분석하고자 한다.
본 연구의 목적은 2014년부터 2015년까지 우리나라에서 발생한 가뭄 기간일 때, 다중분광 광학 센서인 MODIS 자료를 통해 생산된 토양수분과 플럭스타워에서 관측된 토양수분을 비교하여 광학센서 기반으로 산정된 토양수분의 적용성을 평가하고자 한다. 또한 MODIS 산출물인 총 일차생산량(Gross Primary Productivity, GPP)를 통해 이산화탄소 플럭스를 산정하고, 플럭스타워에서 관측된 결과와 비교 분석하여 원격탐사 기반 탄소 모니터링의 가능성을 파악하고자 한다. 농업 가뭄이 발생할 때, 토양수분이 감소하여 식생 생장에 부정적인 영향이 가해지고, 이에 따라 광합성을 활발하게 하지 못하게 되어 이산화탄소 배출이 증가하게 된다는 가설을 증명하고자 한다.
2. 연구 방법
2.1 연구 대상 지역
본 연구의 대상지역은 남한지역을 대상으로 인공위성 자료 기반의 이산화탄소 플럭스(CO2 flux)와 토양수분을 검증하기 위하여 Table 1과 같이 한국수자원조사기술원의 설마천과 청미천 플럭스타워의 CO2 flux와 토양수분 자료를 활용하였다. 설마천 플럭스타워는 경기도 파주시 적성면에 위치하며, 해발고도 293 m의 높은 고지대에 위치한 전형적인 산지지형이다. 청미천 플럭스타워는 경기도 이천시 장호원읍에 위치하며, 논에 위치하고 있다.
Table 1.
Characteristics of the selected sites in study area
2.2 다중분광광학센서 자료
토양수분, 식생지수 및 총일차생산량(GPP)을 산정하기 위해 Table 2와 같이 다중분광광학센서인 MODIS 자료를 활용하였다. MOD11A1의 지표면 온도와 MCD43A3의 지표 반사도를 이용하여 토양수분을 산정하였다. 식생지수는 MOD13A1의 NDVI를 사용하였으며, MOD17A2H의 GPP를 통해 식생 활력에 따른 CO2 발생 여부를 분석하였다. MODIS 자료는 모두 공간 해상도 0.5 km로 제공되며, 시간 해상도는 MOD11A1와 MCD43A3은 1일, MOD17A2H는 8일, MOD13A1는 16일로 제공된다. MODIS 자료는 ftp://e4ftl01.cr.usgs.gov/에서 2001년부터 자료를 획득할 수 있다.
Table 2.
Remote sensing data used in this study
2.3 이산화탄소 플럭스 산정 방법
위성영상으로부터 직접적인 CO2를 산출할 수 없기 때문에 CO2와 GPP사이의 관계를 통해 간접적으로 산출할 수 있다. CO2와 GPP의 관계를 Eq. (1)로 표현할 수 있다.
여기서, NEE (Net Ecosystem Exchange)는 생태계와 대기 간 탄소의 순 교환을 측정하는 것으로, 탄소 생태계 흡수 강도를 나타내는 주요 척도로 활용되며, 플럭스타워에서 측정된 CO2 플럭스 자료를 의미한다(Amthor et al., 2001). Re는 생태계 내에서의 식생의 호흡량을 의미하며 CO2 flux와 GPP사이의 관계는 Re 값을 통해 산출할 수 있다. Re는 Eq. (2)를 통해 산정 된다(Lloyd and Taylor, 1994).
여기서, Rref는 생태계 내에서의 기준 호흡량, Tref는 기준 온도 (10°C), E0는 활성 에너지, T0는 -46.02°C, Ta는 대기 온도이다. Rref와 E0는 경험상수이다.
2.4 다중분광광학센서를 활용한 토양수분 산정 방법
본 연구에서는 Verstraeten et al. (2006)에서 제시한 광학센서 산출물을 활용하여 토양수분을 산출하는 방법을 적용하였다.
여기서, ∆LST는 LSTD과 LSTN의 차이를 의미하며, LSTD은 낮 시간 동안의 지표면온도, LSTN은 밤 시간 대의 지표면온도를 의미한다. ATI는 겉보기 열관성(Apparent Thermal Inertia)을 의미하며, 열 관성은 온도 변화에 대한 임피던스를 설명하는 재료의 물리적 특성이다. α는 지표 반사도를 의미하며, 음수값이 적용되었기 때문에 물리적으로 흡수되는 상태를 나타낸다. ATI 값의 최대/최소값을 통해 정규화 기법을 적용하면 체적 토양수분(m3 m-3)값을 산출할 수 있다. 지표면 온도와 지표 반사도 모두 일별로 제공되기 때문에 토양수분의 경우에도 일별로 산정된다.
2.5 농업적 가뭄 지수
본 연구에서 적용한 가뭄지수는 식생조건지수(VCI)와 표준토양수분지수(Standardized Soil Moisture Index, SSMI)이며, 가뭄이 발생했을 때, 토양수분, 식생활동, 탄소발생의 관계를 분석하기 위하여 활용하였다.
농업적 가뭄상태를 분석하기 위하여 VI를 이용하였으며, MODIS 센서의 MOD13A1에서 제공되는 정규식생지수(NDVI)를 정규화하여 식생조건지수(VCI)를 산정하였다 (Kogan, 1997).
VCI를 같은 기간에 산정된 NDVI와 비교하였을 때, %값으로 표시되며, 이전 연도의 극단값(최대값과 최소값) 사이에 위치한다. VCI 값의 범위에 따라 가뭄의 등급을 결정할 수 있다.
표준 토양수분 지수(SSMI)는 체적 토양수분을 표준 정규화하여 지수로 가공하였으며, SSMI의 산정방법은 Eq. (7)과 같다.
여기서, Z는 표준 정규화를 의미하며, SMmax는 최대 토양수분 값, SMmin은 최소 토양수분 값, SMi는 i 시점의 토양수분 값을 의미한다.
3. 연구 결과
3.1 광학센서 산출물의 비교 분석
3.1.1 토양수분의 시간적 변동 및 지상관측자료와의 비교 검증
광학센서 산출물을 활용한 토양수분 산정 방법을 활용해 토양수분을 산정하고(SMMODIS) 검증 지점인 SMC와 CMC에서 관측된 토양수분(SMfluxtower)와 비교 분석하였다. Fig. 1은 2014년부터 2015년까지 SMMODIS와 SMfluxtower를 비교한 결과로 같은 기간에 관측된 강우량 값과 같이 비교하였다. 해당 기간은 남한 전역에서 극심한 가뭄상태이기 때문에 강우 발생이 평년보다 적지만, 강우가 발생했을 때 토양수분의 값도 같이 증가하는 경향을 확인하였다. 두 지역 모두 2014년 9월-10월, 2015년 7월 이후의 구간에서 상대적으로 낮은 값을 보인다.
Bias와 RMSE의 결과를 보면(Table 3), CMC에서 SMMODIS의 결과가 SMC보다 상대적으로 좋은 결과를 나타내는 반면, 상관성은 SMC에서 다소 높게 산정되었다. 두 지역의 결과 값 모두 통계적 으로 유의성을 보이고 있고. SMC에서 bias가 상대적으로 높게 관측되는 이유는 지형학적인 특징 때문이다. SMC는 산악지형으로 혼효림 지역에 위치하고 있으며, 이러한 지역은 지표 반사도값이 비균질하여 위성영상으로 측정하게 될 경우 bias가 발생한다(Sur et al., 2019).
Table 3.
Statistical comparison between estimated and in-situ measurements of soil moisture
| Site | bias | RMSE | a | Coefficients of Correlation |
| SMC | -1.36 | 3.37 | 0.81 | 0.75 |
| CMC | -0.20 | 1.71 | 0.58 | 0.66 |
3.1.2 이산화탄소 플럭스의 시간적 변동 및 지상관측자료와의 비교 검증
GPP와 CO2의 관계식에 따라 MODIS GPP로부터 CO2 (CO2_MODIS)를 산정하고, 이를 플럭스타워에서 관측된 CO2 (CO2_fluxtower)와 비교 분석하였다. CO2_MODIS와 CO2_fluxtower는 8일 단위 산출물로 구성되었다. Fig. 2는 2014년부터 2015년까지 CO2_MODIS와 CO2_fluxtower를 비교한 결과로 같은 기간에 관측된 강우량 값과 같이 비교하였다. 같은 기간의 토양수분 결과와 비교했을 때, 강우 발생 시점에서 지체시간이 발생한 후에 비슷한 양상을 나타낸다. 이는 토양수분이 식생 활동에 영향을 준 뒤 발생하는 광합성과 관련이 있기 때문에, 일정 시간 지체시간이 발생하는 것이다(Sur and Choi, 2013).
Table 4는 두 검증 지역에서의 CO2_MODIS 와 CO2_fluxtower의 통계 분석 결과이다. 통계분석 결과, CO2 flux의 결과는 토양수분의 결과와 다소 다른 점을 보인다. CO2 flux의 통계분석 결과 모두 통계적 유의성을 보이며, SMC와 CMC 두 지점 모두에서 bias, RMSE, 1:1 그래프의 기울기와 상관계수가 거의 유사하게 산정되었다. SMC의 결과가 좋게 나타난 이유는, 식생 밀집도 때문이다. 두 지점의 식생 밀집도는 잎 면적 지수(Leaf Area Index, LAI)로 판별할 수 있는데, SMC의 경우 LAI는 6-8 사이의 값을 나타내며, CMC의 경우 1-2의 LAI를 나타낸다. LAI의 차이가 다소 크지만, CMC의 경우 관개시스템이 구축되어 있는 논에 위치하고 있어 가뭄상태에 다소 덜 민감하게 반응 한다(Sur et al., 2019). 이러한 특성으로 인해 SMC와 CMC의 통계 분석 결과의 차이가 다소 줄어들었다.
Table 4.
Statistical comparison between estimated and in-situ measurements of CO2 flux
| Site | bias | RMSE | a | Coefficients of Correlation |
| SMC | 0.93 | 5.32 | 0.79 | 0.76 |
| CMC | -1.37 | 6.12 | 0.74 | 0.72 |
3.2 토양수분-식생지수-이산화탄소 플럭스 간 상관 분석
SSMI와 VCI는 60-80%의 가뭄 정확도를 나타낸다(Table 5). SSMI의 경우 토양수분으로 만든 가뭄지수로 현재 토양의 수분상태에 민감하게 반응하여 일반적인 가뭄상태일 때, 높은 정확도를 보인다. 반면 VCI의 경우, 토양수분의 부족으로 인해 식생 활력이 약해지는 것을 포착하는 지수로, 가뭄이 장기적으로 지속될 때 정확도가 높아진다. 농업적 가뭄이 얼마나 오랫동안 발생하는지에 따라 각각의 지수를 적절하게 활용할 수 있다.
Table 5.
Error analysis for the study area
가뭄 발생 기간일 때, 토양수분, 식생지수와 CO2는 밀접한 연관성이 있다. 이 인자들의 연관성 사이에는 지체 시간이 있는 것으로 판단되어 지체상관분석을 수행하였다(Table 6). 각 인자들끼리의 균질한 비교를 위하여 표준 정규화된 인자인 SSMI와 VCI를 적용하고 CO2 flux도 값의 범위를 -1에서 1까지 통일시켰다. 지체상관분석을 통해 SSMI를 기준으로 VCI와 CO2 모두 지체 시간이 발생한 후에 가장 높은 상관성을 보였다. VCI의 경우, 지체 시간 0.5-1개월일 때 가장 높은 상관성을, CO2는 1-1.5개월 지체시간에서 가장 높은 결과를 보였다. 토양수분의 변동이 가장 먼저 시작된 후, 일정시간이 지나고 식생지수와 CO2 flux의 변동이 발생하였다.
Table 6.
Lag correlation for study area
4. 결 론
본 연구에서는 토양수분, 식생 활력과 대기 중의 CO2 농도의 관계를 다중분광센서인 MODIS 산출물을 이용하여 분석하였다. 토양수분의 경우, 고해상도 자료를 이용한 분석을 위해 광학 센서를 이용한 토양수분 산정 방법을 적용하였다. 산정된 토양수분의 상관 분석 결과, 0.7 이상의 상관성을 보였다. 대기 중 CO2 농도의 경우, MODIS GPP를 이용하여 호흡량과의 관계식을 이용하여 산정하였다.
두 플럭스타워 지점에서 관측한 값과 비교 분석하여 높은 적용성을 보였다. 산정된 CO2 플럭스는 설마천 지점의 결과가 상대적으로 적용성이 높게 분석되었는데, 이는 잎면적 지수의 영향으로 설마천이 식생분포 밀집도가 높기 때문이다. 또한 국내에서 발생한 극한 가뭄사상인 2014년과 2015년도에 대하여 토양수분, 식생지수와 CO2 의 교차분석을 통해 각 인자들 사이의 관계성을 파악하였다. 분석 결과 토양수분과 식생지수 사이에는 1개월의 지체시간, 식생지수와 CO2 사이에는 0.5개월의 지체시간이 발생했을 때, 가장 높은 상관성을 보였다. 지표면의 토양수분과 식생지수를 대기 중의 CO2의 관계를 분석함으로써 대기 중의 CO2 농도가 변하게 되면 지표면 및 지표하에 저장되어 있는 탄소 농도에도 영향을 받는 것으로 사료된다.
