RUSLE을 활용한 금강 수변지역의 하천정비사업 전·후의 연간 표토침식량 변화 비교분석

A Comparative Analysis of Annual Surface Soil Erosion Before and After the River Improvement Project in the Geumgang Basin Using the RUSLE

Abstract

본 연구에서는 금강 유역(대청댐~금강하구둑)을 대상으로 RUSLE (Revised Universal Soil Loss Equation)를 활용하여 하천정비사업 전(2007년)과 후(2015년)의 연간 표토 침식량을 계산하였고, 이를 5단계의 class를 구분하여 해당 구간 별 표토침식면적을 산출하고 그 결과를 환경부의 토지피복 변화와 비교·분석하였다. RUSLE을 구성하는 각 인자를 생성하기 위하여 2007년과 2015년의 토지피복도, 국가기본 공간정보, 정밀 토양도, 연평균 강수량 자료 등의 다양한 공간정보자료를 활용하였다. 연구 결과, 연구지역에서의 연간 표토 침식량은 2007년과 대비하여 2015년에 class 1의 면적은 증가하였고, class 2, 3, 5의 면적은 감소하였으며, class 4의 면적은 증가하였다. 이러한 원인으로는 연구지역의 연평균 강수량의 감소, 생태공원조성과 이로 인한 인공시설물 확대, 농경지 감소 등으로 인하여 대부분의 지역에서 연평균 표토 침식량이 감소된 것으로 사료된다.

In this study, the annual surface soil erosion amount of before (2007 year) and after (2015 year) the river improvement projects were calculated using RUSLE (Revised Universal Soil Loss Equation) in the Geumgang basin (Daecheong-Dam to Geumgang Estuary-Bank). After the results were classified into five classes, the results were compared and analyzed with the results of the change in the land cover. In order to generate each factor of RUSLE, various spatial information data, such as land cover maps for 2007 and 2015 years, national basic spatial information, soil map, and average annual precipitation data were utilized. The results of the analysis are as follows: 1) annual surface soil erosion in the study area increased the area of class 1 in 2015 years compared to 2007, 2) the area of class 2, 3 and 5 decreased, 3) the area of class 4 increased. It is believed that the average annual amount of surface soil erosion decreased in most areas due to the reduction of annual average precipitation, the formation of ecological parks, the expansion of artificial facilities, and the reduction of illegal farmland.

1. 서론

표토란 지표면에서 약 20~30 cm 깊이까지의 토양층으로, 유기물 층과 용탈층 등을 포함하며 생태계 유지와 물질순환에 핵심적 구실을 한다. 표토의 침식은 물, 바람, 중력 등과 같은 자연적 현상 또는 인위적인 영향으로 인해 표토가 침식되는 현상을 말하며, 전국 표토침식량 예비조사(2012) 결과, 대한민국의 연평균 표토침식량은 국토의 30% 이상에서 1 ha당 33 톤을 기록하여 OECD평균인 1 ha당 11톤에 비하여 약 3배에 가깝다고 발표된 바 있다.

표토는 다양한 생물의 섭식활동을 위한 생활터전을 제공함과 동시에 양분 및 수분을 공급하여 생태계를 유지할 수 있는 근간이며, 오염 정화 및 물질순환 기능, 탄소의 저장으로 인한 지구온난화 방지 기능, 수자원 저수 및 홍수조절기능 등의 순기능적 측면을 가진 소중한 자원이다. 그러나 이러한 표토가 침식됨으로 인해서 발생하는 악영향으로는 첫 번째로, 다양한 생물의 서식지를 파괴하여 생물 서식환경에 악영향을 끼칠 수 있으며, 두 번째로 표토가 함유하고 있는 다량의 영양염류(N, P등)가 손실되어 토양 비옥도를 감소시켜 식물의 생장기능을 저하시킬 수 있으며, 세 번째로 손실된 영양염류가 인근 하천에 유입되어 부영양화를 초래하여 수질 오염을 일으킬 수 있다. 위와 같이 표토의 침식으로 인한 환경문제는 다각화된 방면에서 야기되고 있다.

일반적인 표토 침식의 경우에는 집중호우, 태풍, 지진과 같은 초자연현상에 의한 표토침식과, 인간 활동의 영향으로 인한 누적된 표토침식으로 구분된다. 반면에 본 연구지역인 금강의 수변 지역의 경우, 2007년 하천정비사업의 대대적인 공사로 인해 수변 지역의 토지 피복형태가 단기간에 직접적인 형태로 변화하게 되었으며, 이러한 토지 피복변화가 표토의 침식에 미치는 영향에 대해 분석한 연구는 매우 미미한 실정이다.

표토 침식량을 정량적으로 산정하기 위해 많은 연구에서는 범용토양유실공식(RUSLE: Revised Universal Soil Loss Equation)을 사용한다. RUSLE 모델의 장점은 토지피복의 변화 및 다양한 침식 조절 요인에 대한 영향을 직접적으로 고려할 수 있고, 다른 모델에 비해 표토 침식량 산정 결과의 신뢰성이 뛰어나다고 알려진 바 있다. 또한, 수계 및 유역에 대해 적용하여 유역대책에 대한 의사결정을 지원하며, 최신 토지피복 특성을 적용하기에 적절하다는 특징이 있어, 본 연구에서는 연간 표토 침식량 산정 및 비교를 위해 RUSLE 모델을 사용하였다.

국내에서 표토 침식량을 분석한 사례로는 Kim et al. (2018c) 이 RUSLE 모델을 이용하여 금호강의 하천정비사업을 전 후로 한 연간 표토 침식량을 산정하여 비교·분석한 연구를 발표된 바 있으며, Choi et al. (2013)이 RUSLE을 이용하여 경주국립공원의 토양침식 위험지역 추정에 대한 연구결과를 발표한 바 있다. 표토는 국지적 집중호우에 의한 산사태의 발생으로 인하여 다량의 표토가 침식되는 경우가 있으며, 이러한 산사태 위험지역을 추정하기 위한 연구가 강원도 평창군을 대상으로 하여 수행된 바 있다(Kim et al., 2018b). 금강유역에 대해서 Kim et al. (2017) 의 논문에서는 한국형 RUSLE공식을 이용하여 유역 전체에 대한 토양침식 우려 지역선정에 대하여 연구한 바 있으며(Kim et al., 2018a; 2019), 논문에서는 본 연구대상지역의 대수층 지하수 생산성을 파악하기 위하여 금강 상류 옥천지역의 토양정보 및 지형정보를 이용하여 원격탐사 및 GIS 기법을 통해 해당지역의 지하수 생산성 추정 연구를 수행한 바 있다. 그러나 대부분의 국내에서 RUSLE 모델을 이용하여 표토침식에 대한 연구들의 경우 특정시점에서의 침식 우려지역을 선정하는데 그쳤으며, 하천정비 사업을 전후로 하여 토지피복 변화와 표토 침식량의 변화를 비교 분석한 선행 국내연구는 매우 희박하다.

본 연구에서는 금강 수변 지역의 하천정비 사업으로 인한 토지피복의 변화에 따른 연간 표토 침식량을 사업전(2007년)과 사업 후(2015년)로 산출하고, 그 결과를 비교·분석하고자 한다.

2. 연구지역 및 연구자료

1) 연구지역

연구대상 지역인 금강유역은 동경 126°40′25″~128° 03′53″, 북위 35°34′47″~37°03′03″사이에 위치하며 지리적으로 한반도의 중서부지역에 위치하고 있는 유역 면적 9,912.15 km², 유로연장은 397.79 km의 우리나라에서 3번째로 큰 국가 하천이다(Fig. 1).

Fig. 1. Location map of study area. (a) broad map of South Korea (Daum aerial photography), (b) detail map of study area with DEM (1:5,000).

금강 수계는 금강 본류 및 미호천, 갑천, 유등천을 포함한 7개 국가 하천 461개 지방하천으로 이루어져 있지만, 본 연구에서는 대청 댐으로부터 금강 하구둑까지의 금강 본류지역을 대상으로 하며 해당 구간의 유로 연장은 약 135.09 km이다(Lee et al., 2006). 하천정비사업으로 인한 토지피복의 변화비교를 위해 사업으로 인해 인위적으로 가장 많은 변화를 겪었던 금강 유역의 수변지역을 연구대상지역으로 하였다. 행정구역상으로는 대전광역시를 비롯하여 충청남·북도, 경기도, 전라북도, 경상남·북도에 총 1광역시, 6개도, 9개시, 18개군, 5개구를 포함하고 있다. 2006년 국토교통부에서 발행한 수자원 장기종합계획에 따른 금강유역의 연평균 기온은 11.6°C이며 연평균 강우량은 1,226.6 mm이다.

표토의 침식현상은 해당 연구지역의 표토의 배수능과 매우 밀접한 관계를 지닌다. 국토교통부의 2011년 금강수계하천기본계획 보고서에 따르면, 금강 유역의 경우 침투율이 크고, 배수가 매우 양호한 토양군이 4,583 km²으로 전체 유역면적의 약 46%를 차지하고 있으며, 배수능이 양호한 토양군이 3,161 km²(32%), 보통의 배수능의 토양 군 1,184 km² (12%), 낮은 배수능의 토양 군 984 km²(10%)으로 구성 되어있다 (MOLIT, 2011).

2) 연구자료

표토의 침식을 정량적으로 평가하기 위해 토지피복지도, 정밀토양도, 수치표고자료(DEM), 강수량 자료를 기반으로 하여 RUSLE모델을 구성하는 각 인자들을 산정하였다. Table 1은 각 인자들을 계산하기 위해 활용된 주제도와 해당 주제도의 원 데이터 형식 및 축척을 표시하였다. 강우인자(R)는 국수자원종합정보시스템 (WAMIS)의 금강유역 14개 관측소의 월별 강수량자료를 수집하여 구축하였다. 토양침식성인자(K)를 산출하기 위해 사용되는 자료인 토양통 및 토양부호는 국립농업과학원에서 제공하는 1:5,000 정밀 토양지도를 통해 획득하였다. 지형인자(LS)는 국토지리정보원에서 제공하는 1:5,000 수치지도의 표고 및 등고선 자료를 이용하여 DEM을 제작하였고 이를 이용하여 지형인자를 산출하였다. 식생피복침식인자(C)의 경우 하천정비사업 전·후 비교를 위해 사업 전의 경우 환경부의 2007년 중분류토지피복지도, 사업 후는 2015년 중분류 토지피복지도를 이용하여 C인자를 계산하였다. 환경부 공간정보서비스(https://egis.me.go.kr)의 내용에 따르면 2007년 전국단위 중분류 토지피복지도의 경우 SPOT-5 위성에 의해 촬영된 영상을 기반으로 제작되었으며, 2015년 금강지역의 세분류 토지피복지도는 항공사진 및 KOMPSAT-2와 KOMPSAT-3, 3A 위성으로부터 촬영된 영상을 기반으로 제작되었으며, 본 연구에서는 제작된 환경부 토지피복지도를 이용하였다. 보전관리인자(P)의 경우 토지피복지도와 DEM으로부터 추출한 경사도를 중첩하여 분석하였다(Kim et al., 2018c).

Table 1. The data used in this study

3. 연구방법

1) RUSLE(범용토양유실공식)

금강유역의 표토 침식량을 산정하기 위해서, 다양한 주제도와 GIS를 기반으로 하여 RUSLE을 적용해 계산하였다. RUSLE 모델은 Renard et al. (1991)에 의하여 제안된 방법으로, 기존의 농업지역에 한하여 강우에 따른 포토의 침식량을 산정하던 방식인 USLE(Universal Soil Loss Equation)에서 계곡 지형 및 유역에서 표토의 침식량을 정량적으로 산정할 수 있도록 개선된 경험식이다. RUSLE은 표토 침식의 주요 인자를 강우(강우 침식인자), 토양의 성질(토양침식인자), 사면의 경사 및 길이(지형인자), 토지의 식생 피복 상태(식생피복침식인자), 경작지 이용상태(보전관리인자)의 5가지로 선정하여, 각 인자의 값을 산정해 정량적으로 연평균 토양 침식량을 표현한다(Fig. 2).

\(A=R \times K \times L S \times C \times P\)       (1)

Fig. 2. Flow chart of the study.

식 (1)에서 A는 연평균 토양 침식량(ton/ha/yr), R은 강우침식인자, K는 토양침식인자, LS는 사면의 경사 및 길이인자, C는 식생피복침식인자, P는 보전관리인자를 나타낸다. 본 연구에 사용된 모든 자료는 사업에 의하여 변화를 하천 범람원지역을 대상으로 하며, 하도로부터 약 200 m 구간을 버퍼로 설정하여 분석하였다.

최종적으로 산정된 연간표토침식량결과(A)는 Roy (2019) 및 Rellini et al. (2019)가 제시한 기준을 참고하여 표토의 침식 정도에 따라 다음과 같이 5 class로 구분하였다. class 1은 표토침식 발생하지 않음(0-0.1 t/ha/yr), class 2는 표토침식 소량 발생(0.1-1 t/ha/yr), class 3은 표토침식 보통 발생(1-10 t/ha/yr), class 4는 표토침식 다량발생(10-50 t/ha/yr), class 5는 위험수준의 표토침식 발생(50 t/ha/yr 이상)하는 지역을 나타낸다.

2) 강우침식인자 R

강우침식인자(R)는 강우에 의한 토양의 침식능력을 의미하며, 강우의 운동에너지와 강우강도에 의해 결정된다. 일정 강도 이상의 강우가 지속될 경우 토양입자가 이탈되고 강우에 의해 형성된 표면 유거수에 의해 토양입자의 이송이 이루어져 토양의 침식이 발생된다. 강우침식인자를 계산하기 위한 연평균 강수량데이터는 국가수자원관리종합정보시스템(http://www.wamis.go.kr)으로부터 획득하였으며, 금강유역에 인접한 12개 시·군의 지점을 선정하여 사업 전(2007년)은 2003년~2007년, 사업 후(2015년)은 2011년~2015년까지의 5년간의 연평균 강수량 값을 이용하였다. Toxopeus (1996)의 식을 이용하여 연평균 강수량 값을 강우침식인자 값으로 변환하였다. 

\(R=3.85+0.35 \times P\)       (2)

식 (2)에서, 강우침식인자(R)은 J/m²이며, P는 연평균 강우량(mm/yr)이다. 본 연구에서 획득한 관측소 별 사업 전·후 5년간의 연평균 강우량과 강우침식인자 값(R)은 Table 2와 같다.

Table 2. Changes of Annual average precipitation and R factor between 2007 and 2015year

3)토양침식인자 K

토양침식인자(K)는 토양입자가 이루는 구조적 안정성과 투수성의 정도에 따른 토양의 침식능력을 나타내는 값이며, 극세사 또는 미사 토양과 같이 입자의 크기가 작을 수록, 유기물의 함량이 적을수록, 투수성이 낮을수록 표토침식 발생 확률이 높기 때문에 이러한 토양에서는 K값이 높게 산정된다. 본 연구에서는 Wischmeier (1976) 이 제시한 토양통별 K값 산정방법에 따라, 국립농업과학원으로부터 획득한 1:25,000 정밀토양도의 토양통값을 이용해 K값을 산정하여 토양침식인자를 계산하였다.

\(K=\frac{2.1 \times 10^{4}(12-O M) M^{1.14}+3.25\left(S_{\Gamma}-2\right)+2.5\left(P_{1}-3\right)}{100}\)       (3)

식 (3)에서, K는 토양침식인자, OM은 토양 내 유기물의 비율(%), M은 토양입자와 토사전체의 비율에 대한 함수, S₁은 토양구조의 등급(1~4), P₁은 투수에 대한 등급(1~6)이다(Lee et al., 2012).

4) 식생피복침식인자 C

식생피복침식인자(C)는 지표면의 피복상태에 따른 침식량을 설명하는 인자로서 일반적으로 식생으로 피복된 지표면의 경우 토양침식이 감소한다(Renard et al., 1997). 활엽수 또는 침엽수와 같은 교목식물의 경우 뿌리를 통해 토양입자의 분리를 막는 역할을 할 수 있으며, 특히 활엽수의 경우 넓은 잎을 통해 토양이 강우에 의해 직접적으로 침식되는 것을 상쇄시키는 역할을 하기 때문에, 산림지역과 같은 경우 C값은 매우 낮은 값을 나타낸다(Jung et al., 2004). 그러나 나대지와 같이 식생이 존재하지 않는 피복의 경우 지표면이 그대로 강우에 노출되기 때문에 1에 가까운 값을 갖는다(Choi et al., 2013). 그러나 주거지역, 교통지역, 공공시설지역의 경우 표토가 아스팔트 또는 콘크리트에 의해 포장 되어있기 때문에 식생피복침식인자 값은 0에 가깝다.

본 연구에서는 금강지역의 하천정비사업 전·후 비교를 위해 환경부의 2007년, 2015년 토지피복지도를 중분류로 수준을 통일시켜 이용하였으며, Kim et al. (2018c)이 사용한 미국 농무성(USDA)의 토지피복 별 식생피복침식인자값을 이용하여 값을 산정하였다(Table 3).

Table 3. Value of C factor according to land cover types

5) 보전관리인자 P

보전관리인자(P)는 농경지가 경작형태에 따라 토양침식에 미치는 영향을 설명하는 인자로서 경작형태의 경우 보통 등고선 경작(contour system), 계단식 경작(terrace systems)으로 구분된다. 등고선 경작은 밭, 과수원 등의 경작형태를 띄며, 계단식 경작은 논과 같이 경사가 없는 경작 형태를 띈다. 등고선 경작형태는 계단식 경작과 비교했을 때, 일정한 방향과 경사를 갖기 때문에 강우가 발생할 경우 토양침식에 취약하며, 경사가 높아질수록 P값은 1에 가까워진다(Hwang et al., 2010). 본 연구에서는 국토지리정보원으로부터 획득한 표고점과 1:5,000 등고선으로부터 ArcGIS software를 이용해 DEM을 제작하였으며, 경작지에 관한 정보는 환경부의 토지피복지도 중분류 정보를 활용하였다. 경사도에 따른 경작형태별 보전침식인자 값은 Table 4에 표기되어 있다. 표토의 침식에 미치는 영향을 나타낸다. 사면의 경사길이는 강우에 의한 흐름이 시작되는 지점으로부터 퇴적이 일어나는 곳까지의 길이를 말하며, 경사길이가 길어지고 경사도가 높아질 표면의 강우 흐름이 증가함에 따라 더 많은 에너지가 전달되기 때문에 침식량이 증가하게 된다.

Table 4. Value of P factor according to cultivation type and slope grade

\(L_{(i, j)}=\frac{\left(A_{i j}+D^{2}\right)^{m+1}-A_{(i, j)}^{m+1}}{D^{m+2} \times x_{(i, j)}^{m} \times 22.13^{m}}\)       (4)

\(m=\frac{\beta}{\beta+1}, \beta=\frac{(\sin \theta / 0.0896)}{\left(2.96 \times \sin ^{0.79} \theta+0.56\right)}\)       (5)

식 (4)와 (5)에서, L_(i,j)은 (i, j)번째 셀의 경사길이인자, A_(i,j)는 (i, j) 번째 셀이 나누어 받게 되는 상류기여면적, D는 셀의 크기(m²), m은 경사길이지수, x_(i,j)는 흐름방향에 직교하는 등고선의 길이이다(Ozsahin et al., 2018). 복잡하고 불규칙한 사면의 한 지점에서 흐름의 차이는 그 지점에 기여하는 상부의 유역면적에 따라 달라진다는 특징에 따라 ArcGIS 10.7과 Saga GIS프로그램의 L-S factor모델을 이용하여 지형인자를 추출하였다(Ozsahin et al., 2018).

4. 연구결과

1) 하천정비사업 전후 토지피복변화 비교

금강의 수변지역은 2008년부터 시행된 하천정비정비사업으로 인하여 불법으로 경작되던 농경지를 제거하고 근린·친수 생태공원으로 새롭게 조성되었다. 따라서 논, 밭 또는 하우스의 형태로 경작되던 농경지가 일부는 공공시설지(자전거도로, 광장, 화장실 등)로 일부는 일년생 또는 다년생 초본식생지로, 일부는 식생천이가 진행되지 않은 나대지로 변화하였다.

Fig. 3은 금강의 세종시 지역으로 사업 전인 2007년과 사업 후인 2014년의 KOMPSAT-2 & 3 위성의 RGB영상과, 이를 이용하여 제작된 환경부의 2007년과 2015년의 토지피복지도를 나타내었다. (a)와 (c)의 2007년의 경우 하천에 인접한 대부분 지형의 토지피복이 경작지(agriculture) 로 이용되거나 또는 나지(barren)의 상태였는데, (b)와 (d)의 2014~2015년의 경우 초지(grass)로 변화하였음을 확인할 수 있다. 국립생태원의 보고서(2018)에 따르면 금강의 하천정비사업 1~7공구(대청댐~금강하구 둑) 구간에서 2007년과 2017년을 비교한 결과, 농경지의 약 92%가 감소하였다고 분석한 바 있다(NIE, 2019).

Fig. 3. Satellite images of KOMPSAT-2 at 2007 year (a) and KOMPSAT-3 at 2014 year (b), and the Land cover map at 2007 year(c), at 2015 year (d) from the Ministry of Environment, respectively.

Table 5는 연구지역에 대해 환경부로부터 획득한 2007년 토지피복지도와 2015년의 토지피복지도의 각 토지피복에 따른 면적을 비교한 결과이다. 초지와 나지가 약 350%와 약 118% 증가로 가장 큰 증가율을 나타났고, 밭과 논이 약 56%와 약 35% 감소로 가장 큰 감소율을 나타냈다. 이밖에 수역, 인공시설지, 활엽수림, 과수원 면적이 2007년 대비 2015년에 증가하였고침엽수림과 혼효림의 피복면적은 감소하였다. 이러한 토지피복의 면적변화는 결과적으로 연간표토침식량 결과를 정성적, 정량적으로 해석하는데 있어 매우 중요하다.

Table 5. Changes in area by type of land cover between 2007 and 2015year

2) 금강유역의 연간표토 침식량 결과

Fig. 4는 금강유역의 연간표토침식량(A)을 산정하여 하천정비공사 전 2007년 (a)과 후 2015년 (b)의 결과를 나타냈으며, (a)와 (b)의 우측은 공주시의 2007년과 2015년의 표토침식량의 변화를 육안확인이 가능하도록 확대 표현하였다

Fig. 4. The results of RUSLE between 2007year (a) and 2015year (b), respectively.

Table 6의 결과에 따르면, 하천정비사업 전인 2007년과 비교하여 2015년에 class 1은 495.6 ha(약 2.50%) 증가, class 2는 440.4 ha(약 12.46%) 감소, class 3은 81.6 ha(약 13.25%) 감소, class 4는 26.6 ha(약 53.95%) 증가, class 5는 0.1 ha (약 32%) 감소로 나타났다. Class1에 해당하는 지역의 토지피복은 주로 수역, 습지, 도심지이며, class 2에 해당하는 지역의 토지피복은 농경지로 나타났으며, class 3에 해당하는 지역의 토지피복은 초지이며, class 4, 5에 해당하는 지역의 토지피복은 주로 나대지로 나타났다.

Table 6. The results of annual average of surface soil erosion between 2007 and 2015year

5. 토의

본 연구의 연간 표토 침식량 결과는 다음과 같이 정리할 수 있다. 1) class 1의 증가, 2) class 2, 3, 5의 감소, 3) class 4의 증가.

첫 번째로, class 1의 증가는 Table 5의 River 및 Commercial area의 증가와 관련이 깊다고 사료된다. 2012년 금강의 하천정비 사업 이후 세종보, 공주보, 백제보가 건설되었고, 일정수량 이상의 관리수위를 유지함에 따라서 2007년 당시 강에 근접한 부분의 모래톱 및 나지가 수몰되었기 때문에 2015년 해당지역의 표토 침식량은 0의 값을 나타낸다. 또한, 사업 전의 초지와 나지가 공원조성사업에 의해 도로 및 인공시설물에 포장되었기 때문에 표토 침식량 값이 0을 나타낸다. 위와 같은 원인으로 인하여 class 1의 표토 침식면적이 증가하였다고 사료된다.

두 번째로, class 2, 3, 5의 면적감소는 연평균 강수량의 감소, 불법 경작지의 감소, 초지의 증가와 관련이 깊다고 사료된다. Table 2에서 나타난 바와 같이 지역별 하천정비사업 전후 5년 기간의 평균 강수량 비교 시, 보은군을 제외한 연구지역 전 지역에서 약 10% 이상의 연평균 강수량감소를 확인할 수 있다. 따라서, 연평균 강수량의 감소는 곧 표토침식에 큰 영향을 미치는 R 인자값의 감소와 비례하기 때문에, 전체적인 구간에서 연평균 표토 침식량이 감소하였음을 추정할 수 있다. 또한, Table 5에서 나타난 바와 같이 불법 경작지(논, 밭, 과수원)의 합산 면적이 약 49%이상 감소되고 초지 면적이 약 350% 이상 증가하였는데, 초지의 C 인자 값은 0.05이고 논, 밭, 과수원의 C 인자 값은 각각 0.3, 0.24, 0.2로, 표토침식에 더 큰 영향을 미치는 경작지가 감소하고 상대적으로 영향이 작은 초지가 증가하였기 때문에 이러한 원인으로 인해 class 2, 3, 5의 전체적인 표토침식면적은 감소하였다고 사료된다.

세 번째로, class 4의 증가는 나대지의 증가와 관련이 깊다고 사료된다. 있다. 위에서 언급한 바와 같이 Table 5에서 나타난 바와 같이 불법 경작지면적은 약 49% 이상 감소되었고 이중 일부 지역은 생태공원으로 조성되어 초지가 면적이 증가하였으나, 일부 지역은 불법 경작지를 제거한 상태 그대로 방치되어 나지면적이 증가하였다. 이러한 정량적 결과는 Table 5에서 나지면적이 하천정비 사업전과 대비하여 사업 후 약 117% 증가하였음을 확인할 수 있다. 불법 경작지(논, 밭, 과수원)의 C 인자 값은 각각 0.3, 0.24, 0.2이고, 나지의 C 인자 값은 1.0으로, 표토침식에 영향이 상대적으로 작은 경작지 면적이 감소하고, 가장 큰 영향을 미치는 나지면적이 증가하였기 때문에, 결과적으로 class 4의 표토침식면적은 증가하였다고 사료된다.

본 연구의 결과를 도출하고 분석하는 과정에서 한가지 한계점이 존재하였는데, 이는 연구에 사용된 환경부 토지피복지도의 축척이 다르기 때문에 정량적인 비교에서 오차가 발생할 가능성이 존재한다는 점이다. 해당 지역에서 2007년의 1:5,000의 세분류 토지피복지도는 아직 제작되기 전이기 때문에 2015년의 세분류 지도를 1:25,000 중분류 수준으로 down-sampling하여 사용하였으나 그 과정에서 오차가 존재할 수 있다. 지역이 방대하여 동일한 시기의 영상수집이 어려워 환경부의 토지피복지도를 사용하였으나, 향후에는 KOMPSAT영상을 활용하여 직접 2007년과 2015년의 토지피복지도를 제작하여 연구를 수행한다면 보다 정확한 표토 침식량결과를 산출 가능할 것이라 생각된다.

6. 결론

본 연구에서는 금강 유역(대청댐~금강하구 둑)을 대상으로 RUSLE를 활용하여 하천정비사업 전(2007년)과 후(2015년)의 연간 표토 침식량을 계산하였고, 이를 5단계의 class로 구분하여 해당 구간 별 표토침식면적을 산출하고 그 결과를 토지피복의 변화와 비교·분석하였다. 결과적으로 본 연구지역에서의 연간 표토 침식량은 2007년과 대비하여 2015년에 class 1의 면적은 증가하였고, class 2, 3, 5의 면적은 감소하였으며, class 4의 면적은 증가하였다. 하천정비 사업을 통해 하천근처의 불법 경작지를 제거한 뒤 생태공원을 조성하여 초지와 산림이 증가한 지역은 연간 표토 침식량이 감소하였고, 경작지 제거 후 방치되어 나지로 변한 지역은 연간 표토 침식량이 증가하였으며, 도로 및 인공시설물이 조성되었거나 다목적보의 건설로 인해 수량이 증가하여 물에 잠긴 지역은 표토 침식량이 0에 가깝게 감소하였음을 본 연구결과를 통해 알 수 있었다. 표토의 침식은 생물의 서식공간문제, 영양염류 유입에 의한 수질문제, 하천의 미세플라스틱 유입 문제 등 다양한 환경문제를 야기할 수 있다. 국가의 하천정비사업으로 인해 수변 공간의 토지피복이 인위적으로 변화되었고, 이에 따른 영향으로 변화된 표토의 침식량을 정량적으로 비교 및 분석한 사례는 국제적으로도 찾아보기 힘들며, 매우 큰 의미가 있는 연구결과 자료가 될 수 있다. 본 연구의 결과는 향후 국가 또는 지방하천 정비사업에 있어 다양한 환경문제의 원인이 되는 표토의 침식 및 관리에 대한 의사결정을 지원하는 도구로 활용될 수 있을 것이다.