1. 서 론
토양유실은 강우 발생 시에 빗방울의 타격 및 지표 유 출수에 의해서 토양 표면의 토양 입자가 이동하는 현상이 다. 그러나 도시화 등과 같은 토지 이용 현황이 인위적으 로 바뀌거나 혹은 기후변화에 따른 강우 유형의 변화에 따라 토양입자의 거동은 다양하게 바뀔 수 있다. 경작지 에서의 토양유실은 토양 자원의 손실로 작용하며, 또한 유 실된 토양이 강이나 하천으로 유일될 경우에는 탁도를 증 가시킬 뿐만 아니라 토양 입자에 의해 이동하는 영양염류 로 인해 물 오염의 원인으로 작용할 수 있다.
이러한 토양유실의 양을 추정하기 위해 범용토양유실공 식(Universal Soil Loss Equation, USLE) (Wischmeier and Smith, 1965; Wischmeier and Smith, 1978)이 오 랫동안 사용되고 있다. 범용토양유실공식은 단위 면적(ha) 에서 발생할 수 있는 토양유실 발생 가능성(potential soil loss, Mg/ha/yr)을 추정하기 위한 공식으로, 강우침식능인 자(rainfall erosivity factor; R factor, MJ·mm/ha·yr·hr), 토양침식성인자(soil erodibility factor; K factor, Mg·hr/ MJ·mm), 지형인자(slope and slope steepness factor; LS factor, dimensionless), 지표피복인자(crop and cover management factor; C factor, dimensionless), 보전관리 인자(conservation practice factor; P factor, dimensionless) 의 5개 인자를 필요로 한다(Wischmeier and Smith, 1965; Wischmeier and Smith, 1978). 이 공식은 필요로 하는 5개인자의 곱만으로 토양유실 발생 가능성을 추정하기 때문에 모형의 활용성 측면에서 복잡하지 않다. 그러 나 이 공식에 의한 토양유실 발생 가능성 추정치는 적용 대상 지점에서 발생할 수 있는 실제 토양 침식량을 정의 하기 위한 것이 아니기 때문에 공식에 의해 예측된 값에 대한 유달률(sediment delivery ratio; SDR) 보정이 필요 하다(Wischmeier and Smith, 1965; Wischmeier and Smith, 1978; Park et al., 2007; Kim et al., 2017).
즉, 범용토양유실공식 및 이에 근거한 모형에 의한 토 양유실 예측값은 토양이 유실될 수 있는 가능성을 나타내 는 지표이기 때문에, 유역의 유출구까지 유실된 토양의 양 과의 비교를 위해서는 유달률을 산정하여 보정할 필요가 있다. 유달률이란 유역의 유출구에서의 유사량을 유역의 총 토양유실량으로 나눈 비율로 정의된다. 유달률은 유역 의 크기, 유역의 형상, 유역의 경사, 강우 유형, 지표유출 량, 첨두유량, 토지이용 조건, 작물의 종류 및 생장상태, 토양 속성 등과 같은 조건의 영향을 받으나 위의 모든 조건을 고려하여 유달률을 산정하는 것은 불가능하다.
유달률 산정 방법은 유역의 면적이나 경사에 계수를 곱하여 산정하는 방법(USDA, 1972; Vononi, 1975; Williams, 1977), 유역의 curve number와 유역 면적을 고려하는 방법(Williams and Berndt, 1977) 및 GIS의 격자자료별 유달률 산정 방법과 같은 매우 복잡한 형태로 개발되어(De Rosa et al., 2016) 제안되고 있다. 그러나 이러한 방법들은 토지이용현황 및 토양의 속성이 다양한 유역의 규모에 대한 적용 시에 필요한 것으로 보인다. 토 양유실을 예측하기 위한 대상 지점 하나의 토지이용현황 및 토양 속성으로 그 규모가 작을 경우에는 대상 지점의 조건을 반영할 수 있는 계수를 결정하여 적용하는 것이, 범용토양유실공식의 복잡성 유지 측면에서 볼 때, 보다 합 리적일 것으로 보인다.
우리나라에서는 토양유실에 대한 심각성 및 관리 필요성을 인지하여 2012년에 ‘표토의 침식 현황 조사에 관한 고시’(Ministry of Environment, 2012)(이하 ‘고시’)를 제 정하고 2013년부터 시행하여 국내에서의 토양유실 발생 현황을 파악하고 있다. 이 고시는 예비조사와 현장조사로 구분되며, 두 단계 조사 모두 범용토양유실공식에 근거하 도록 하고 있다. 예비조사에서는 전국을 대상으로 하여 토 양유실이 심각할 것으로 추정되는 지점을 선별하는데 주 된 목적을 두고 있으며, 현장조사는 예비조사 결과에 근 거하여 현장을 방문하여 대상 필지의 토양 속성, 지형적 인 조건(경사장 및 경사도), 토지피복현황 등을 조사하여 현장 조건에 부합하는 범용토양유실공식의 인자를 결정하 여 토양유실량을 추정하도록 하고 있다.
예비조사 단계에서는 전국을 대상으로 하기 때문에, 이 미 산정된 강우침식능인자(Jung et al., 2004), 정밀토양도 에 근거한 토양침식성인자, 수치표고모형을 이용한 지형 인자, 토지이용도를 이용한 지표피복인자, 수치표고모형 및 토지이용도를 이용한 보전관리인자를 적용하는 것이 타당할 것으로 보인다. 그러나 현장조사에서는 현장을 방 문하기 때문에 대상 지점에 대한 자료를 수집하여 각 인 자를 결정하는 것이 타당한 것으로 보이나 강우침식능인 자는 예비조사 시에 사용되었던 인자를 그대로 적용하도 록 하고 있어, 토양유실 가능성 추정에 한계가 있는 것으 로 보인다. 이와 함께 현 고시의 현장조사는 조사대상 필 지에서 발생하는 실제 토양유실량을 측정하고, 범용토양 유실공식에 의한 토양유실 발생 가능성을 추정하도록 되 어 있다. 그러나 범용토양유실공식에 의해 예측된 토양유 실 발생 가능성 추정치는 연 단위이기 때문에 1년에 각각 1개의 값만이 주어지게 된다. 따라서 범용토양유실공식으 로는 현장에서의 토양유실 실측치와 예측치를 비교, 평가 하는 과정이 이루어지기 어렵다.
본 연구는 현장 적용성이 높은 범용토양유실공식에 근거하여 토양유실 발생 가능성을 추정하고, 이를 실측 토 양유실량과 비교하고 보정할 수 있는 방안을 제시하여 환 경부 고시를 개정하고 현장에서 보다 정확한 토양유실을 예측할 수 있는 기반을 제공하는데 그 목적을 두고 있다.
2. 연구방법
2.1. 연구대상지역
환경부의 ‘표토의 침식 현황 조사에 관한 고시’에 따라 한국환경공단에서는 2014년부터 2018년까지 금강수계 2 지점, 한강 수계 2지점, 섬진강 수계 2지점, 낙동강 수계 2지점에 대한 토양유실량 모니터링 사업을 수행하였다. 이 모니터링 지점에서의 토양유실량 측정은 필지에서 발생된 토양을 전량 회수하여 유실된 토양의 양(Mg)을 산정하는 방식으로 운영되고 있다. Fig. 1(1)에서와 같이 필지에서 유실된 토양은 집수플럼부(Fig. 1(2))를 거쳐서 침사지통 (Fig. 1(3))에 수집된다. 그리고 집수플럼부에서의 유실된 토양의 양(Mg)과 침사지통에서의 유실된 토양의 양을 산 정하기 위하여 부유사 농도(mg/L)와 유출량(L)를 측정하 고, 침전된 토양의 양(Mg)을 측정하여 대상 필지에서 유 실된 토양의 양(kg)을 산정하고 있다.
본 과업에서는 각 8개 필지에서 2014년부터 측정된 토 양유실량 자료를 한국환경공단으로부터 제공받아 실측 토 양유실량 자료로 이용하였다. 전체적인 모니터링 기간은 수년에 해당하나 필지 단위 모니터링이기 때문에 강우 시 유출이 발생하지 않았거나, 지점에서 모니터링 수행 기간 이 짧아 자료의 개수가 적은 경우도 있었다. 이런 이유 때문에, 본 과업에서는 총 8개 지점 중에서 2014년부터 모니터링이 수행된 섬진강 유역 2개 지점(Fig. 2(1))에 대 해서만 토양유실양 예측이 이루어질 수 있었다. 그리고 섬 진강A와 섬진강B 지점 모두 모니터링은 2014년 7월 18 일부터 수행되었는데, 강우 시 유출과 토양유실이 발생한 경우는 섬진강A 지점이 20회이고 섬진강B지점이 39회이었다. 실측 자료의 개수가 제한적이므로 이상치 제거 없이 모든 자료를 이용하여 토양유실량 예측값과 비교하였다.
섬진강 유역 내에 위치한 모니터링 지점은 전라남도 화 순군 북면에 위치한 섬진강A(Fig. 2(2))와 전라남도 화순 군 남면에 위치한 섬진강B(Fig. 2(3))로, 섬진강A와 섬진 강B 지점들 모두 농경활동이 이루어지지 않아 초지에 가 까운 상태이며, 면적은 각각 0.0091 ha와 0.0076 ha이며, 평균 경사도는 섬진강A는 21.97%였으며 섬진강B는 20.70%였다.
2.2. 한국형 토양유실공식(KORSLE)의 인자 정의
범용토양유실공식(USLE)은 연평균 토양유실 가능성 추 정치를 예측하기 위한 공식으로 장기 평균 토양유실량 또 는 연간 토양유실량 예측에 이용될 수 있다. 이러한 경우年단위 값인 예측치를 현장 실측치와 비교하기가 어렵다. 이러한 근본적인 문제점을 해결하기 위해 월별로 토양유 실 가능성을 추정할 수 있는 한국형 토양유실공식(Korean Soil Loss Equation; KORSLE)이 최근 개발되어 적용된 바 있다(Sung et al., 2016; Kim et al., 2017; Yu et al., 2017). 이 한국형 토양유실공식은 범용토양유실공식의 한계점인 연단위의 토양유실량 발생 가능성 추정치를 월 단위로 예측할 수 있도록 하고 있으며, 또한 최근 혹은 임의 강수자료에 의해 월별 강우침식능인자를 산정할 수 있도록 하고 있다. 이와 함께 지표피복상태에 영향을 미 치는 작물의 생장을 고려할 수 있도록 4대강 수계별 및 월별 지표피복인자 데이터베이스를 제공하고 있다(Kongju National University, 2018). 범용토양유실공식에 의해 예 측되는 토양유실량 발생 가능성 추정치는 연간 1개의 값 만이 산정될 수 있기 때문에, 본 연구에서는 월별 토양유 실량 발생 가능성 추정이 가능한 한국형 토양유실공식에 의해 모니터링기간에 대한 토양유실 발생 가능성 추정치 를 산정하였다.
한국형 토양유실공식은 국내 47개 지점에 대해서 월 강 수량(mm)에 의한 월 단위 강우침식능인자를 산정할 수 있다. 본 연구의 대상지역인 섬진강A 지점은 기상청의 광 주 지점과 인접한 지역이고 섬진강B 지점은 순천 지점과 인접한 지역이므로 강수량 자료는 각각 광주와 순천 지점 의 자료를 이용하였으며, 두 기상청 지점 모두 광주 지점 에 해당하는 강우침식능인자 산정식을 적용하였다(식 1).
여기서, Ri는 월별 강우침식능인자(MJ·mm/ha·yr·hr)이며, Pi는 월 강수량(mm)이며, M은 월(1월의 경우 1, 12월의 경우 12)이다.
토양침식성인자(Mg·hr/MJ·mm)는 한국환경공단의 모니 터링 수행 결과를 적용하여 섬진강A 지점은 0.023, 섬진 강B 지점은 0.047을 적용하였다. 지형인자는 지점의 경사 도 및 경사장에 의해 산정될 수 있는데 대상 지점을 방 문하여 GPS 장비에 의해서 경위도 및 고도 자료를 수집 하여(Fig. 2(2), Fig. 2(3)) 대상 지점에 대한 1m 해상도 를 가지는 수치표고모형을 제작한 뒤에, 현재 고시에서 제 시하고 있는 지형인자 산정식(Wischmeier and Smith, 1978)(식 (2))을 이용하는 모형인 ArcGIS-based Sediment Assessment Tool for Effective Erosion Control (ArcSATEEC) (Yu et al., 2017)을 이용하여 지형인자를 산정하였다. 지표피복인자는 지표의 상태를 반영하기 위 한 인자로 섬진강A와 섬진강B 지점 모두 초지의 상태이 므로 0.04를 적용하였다. 보전관리인자는 농경지에서의 경 작 관리 형태를 반영하기 위한 것인데 대상 지점은 초지 이므로 모두 1.00을 적용하였다(Kongju National Uni-versity, 2018).
여기서, LS는 지형인자이며, λ는 경사면 투영 거리(meter), β는 경사도(radian), m은 경사에 따른 변수이다. 범용토양유실공식과 이에 근거한 모형은 적용대상 지역에서 발생할 수 있는 토양유실의 가능성을 추정하기 위한 것으로 실제 발생량과 비교하기 위해서는 보정이 필요하 다. 이러한 보정은 한국형 토양유실공식에 의해 예측된 토 양유실량에, 강우 시 필지 조건을 반영할 수 있는 계수를 결정하여 적용하고 실제로 발생하는 토양유실량과의 비교 를 통해 가능할 것으로 판단된다. 이에 본 연구에서는 한 국형 토양유실공식에 의해 예측된 토양유실 발생 가능성 추정치에 각 지점에서의 조건을 반영할 수 있는 보정계수 를 산정하여 예측 토양유실량(Mg)과 실측 토양유실량을 비교하였다.
3. 결과 및 고찰
3.1. 한국형 토양유실공식의 인자 산정 결과
본 연구에서는 실측 토양유실량과 한국형 토양유실공식 에 의한 예측 토양유실량을 비교하였는데, 한국형 토양유 실공식은 범용토양유실공식의 기본 형식을 따르면서 월별 에 대한 토양유실 가능성 예측이 가능하다. 한국형 토양 유실공식의 강우침식능인자는 기간에 대한 누적 강우량에 의해 결정되기 때문에 강우량의 크기에 따라 인자의 크기 가 변화한다. 섬진강A 지점의 경우 모니터링 시작 시점부 터 유실된 토양의 양이 측정이 된 첫 번째 날인 2014년 7월 18일의 경우, 기상청의 광주 지점으로부터 누적 강수 량은 69.6 mm였으며 식(1)에 의해서 계산된 강우침식능인 자(R)는 272.1 MJ·mm/ha·yr·hr였다. 동일한 방식으로 2017년 7월 18일부터 2017년 11월 6일까지의 강수량에 의한 강우침식능인자가 산정이 되었으며, 이로 인해 강수량에 의해 변화하는 강우침식능인자가 섬진강A와 섬진강B 지 점들에 대해서 산정되었다(Table 1, Table 2). 토양침식성 인자(K)는 고시에 제시된 바와 같이 현장조사 시에 토양 표본을 채집하여 유기물 백분율(%), 토양 유실에 대한 입 경 특성 함수, 토양 구조 지수, 토양 투수 등급, 미사 백 분율(%), 극세사 백분율(%), 점토 백분율(%)에 의하여 결정이 되는데, 본 과업에서는 한국환경공단으로부터 제공 받은 토양침식성인자를 적용하였다. 지형인자(LS)의 경우 대상 지점의 경사의 조건을 반영하기 위한 것으로 지점별 로는 다르게 적용이 되었으나, 한 지점에 대해서는 시기 와 상관없이 동일한 값이 적용되었다.
지표피복인자(C)는 작물의 생육 등에 의해 지표의 피복상태가 변화하기 때문에 한국형 토양유실공식에서는 작물 에 따라 다르게 적용이 되나, 섬진강A와 섬진강B 지점 모두 농경지가 아니기 때문에 초지에 해당하는 값이 시기 에 상관없이 동일하게 적용이 되었다. 보전관리인자(P)는 농경지가 아닌 지역이기 때문에 1.00의 값이 시기와 상관없이 동일하게 적용되었다. 즉, 본 연구에서 적용된 한국형 토양유실공식의 각 인자는 강우침식능인자는 강수량에 따라 변화하는 값이 적용이 되었으며, 토양침식성인자와 지형인자는 지점별로는 다르나 각 지점에서 시기별로는 동일하게 적용이 되었으며, 지표피복인자와 보전관리인자 는 두 지점 모두 동일하게 적용이 되었다.
3.2. 보정계수 정의
한국형 토양유실공식은 범용토양유실공식에 근거한 모 형이고, 월별 토양유실 발생 가능성을 추정하기 위한 값 을 제공하기 때문에, 실제 발생한 토양유실량과 비교하기 위해서는 예측된 값에 대한 조정이 필요하다. 따라서 대상 필지의 조건을 검토하여 한국형 토양유실공식에 의해 예측된 값을 조정할 수 있도록 하였다. 이에 각 필지에서 발생한 총 토양유실량(Net soil)에 영향을 크게 줄 수 있 는 요소를 검토하였는데, 강수량(Rainfall), 유출량(Runoff), 집수부에서의 토양유실량 측정치(Measured soil at water collecting basin; SoilWCB), 침사지에서의 토양유실량 측정 치(Measured soil at grit chamber, SoilGC), 침사지에서 의 부유사량(Suspended sediment yield, SoilSS)에 대해서 상관관계(Correlation coefficient)를 검토하였다. 섬진강A과 섬진강B 지점 모두 총 토양유실량은 강수량과 낮은 상관 관계를 보였다. 그리고 유출량은 집수플럼부 퇴적량과 낮 은 상관관계를 보였으나, 총 토양유실량을 포함하여 다른 퇴적량과는 높은 상관관계를 보였다(Table 3, Table 4). 따라서 본 연구에서는 총 토양유실량 예측에 있어 유출량 (L)에 관한 함수로 표현할 필요가 있는 것으로 판단하였다.
Table 3.
The Pearson correlation coefficients between the measured soil erosions and parameters affecting soil erosions in Seomjingang A site
Table 4.
The Pearson correlation coefficients between the measured soil erosions and parameters affecting soil erosions in Seomjingang B site
이에 Table 1과 Table 2에 제시된 한국형 토양유실공 식의 각 인자의 곱에 의해 예측된 토양유실 가능성 추정 치(AKORSLE)와 상관관계가 높은 것으로 판단된 유출량(Q) 을 보정계수(C1)에 관한 함수를 정의하였다. 그리고 이 보 정계수가 적용된 예측 토양유실량(AEst)을 실측 토양유실 량(AMsd)와 비교하면서 Nash-Stucliff Efficiency(NSE), 결 정계수(Coefficient of Determination(R2), t-test에 의한 p-value로 예측치에 대한 신뢰도를 판단하였다.
섬진강A 지점의 경우 이용 가능한 토양유실량 측정 자료는 2014년 7월 18일부터 2017년 11월 6일까지 20개 이며 측정된 토양유실량의 범위로 볼 때, 보정 및 검정 자료로 분할하기 어려운 것으로 판단하여 전체 자료를 이용하여 한국형 토양유실공식에 의한 토양유실 가능성 추 정치에 대한 조정이 이루어졌다. 섬진강B 지점의 경우, 2014년 7월 18일부터 2018년 9월 6일까지 총 39개의 자 료가 이용 가능하였으며 실측 토양유실량 값의 범위와 자 료의 개수를 고려하여 2014년 7월 18일부터 2016년 4월 26일까지 22개 측정 자료를 이용하여 한국형 토양유실공 식에 의한 토양유실 가능성 추정치에 대한 보정이 이루어 졌으며, 2016년 6월 23일부터 2018년 9월 6일까지 17개 측정 자료를 이용하여 보정과정에서 결정된 계수에 대한 검정이 이루어졌다.
섬진강A 지점에 대한 조정은 식(3)과 같이 한국형 토양 유실공식에 의해 예측된 토양유실 가능성 추정치(AKORSLE) 에 면적(Area)을 곱한 뒤에, 유량(Q)과 계수(C1)를 곱하여 토양유실량(Mg)을 예측할 수 있도록 하였는데, C1은 강우 량이 80 mm 이상일 경우에는 0.8로, 그 이하일 경우에는 3.6으로 정의되었다. 섬진강B 지점에 대한 조정은 식(4)와 같이 한국형 토양유실공식에 의해 예측된 토양유실 가능 성 추정치에 상용로그를 취한 뒤에 면적을 곱한 뒤에, 유 량과 계수(C2)를 곱하여 토양유실량(Mg)을 예측할 수 있 도록 하였다. C2는 강우량이 400 mm 이상인 경우에는 66,000, 140 ~ 400 mm인 경우에는 6,600, 140 mm 미만 인 경우에는 3,300으로 정의되었다. 이 식(3) 및 식(4)의 기본 형태와 강우량의 범위에 따른 C1과 C2의 정의는 한 국형 토양유실공식에 의한 토양유실량(Mg)과 실측 토양 유실량(Mg)을 비교하면서 시행착오법(try and error method)에 의해 결정된 것이다.
3.3. 토양유실량 보정 및 검정
한국형 토양유실공식에 근거하여 섬진강A 지점은 식(3) 에 의해 예측 토양유실량으로 결정하였으며, 섬진강B 지 점은 식(4)에 의해 예측 토양유실량으로 결정하였다. 이들 예측 유실량을 실측 토양유실량과 비교하여 KORSLE 모델의 현장 적용 가능성을 판단하였다. 모형에 의한 예측값과 실측값을 비교할 때 수문 모형에 의한 예측값의 유 의성 여부의 판단은 NSE, 결정계수, t-test에 의한 p-value 등에 의해 이루어질 수 있다. 유의성 (significance)에 대한 기준은 Ramanarayaman et al.(1997)은 NSE가 0.4 이상이고 결정계수가 0.5 이상일 때라고 제안하였으며, Donigain(2000)은 결정계수의 범위에 따라 만족 단계를 제시하였는데 0.6 이상일 때 만족 가능(Fair)으로 제안하였으며, Engel et al.(2007)은 유사량(Sediment yield)에 대해서 NSE가 0.5 이상이고 R2가 0.6 이상이며 상대오차 가 20% 이하일 때라고 제안한 바 있다. 본 과업에서는 NSE와 결정계수가 각각 0.4와 0.5 이상일 경우에 이를 도식화하고, t-test에 의한 p-value에 의해서 예측값에 대 한 유의성 여부를 결정하였다. p-value는 유의확률 95%(유의수준 α=0.05)에 대해 검토하여 p-value가 유의수준보 다 클 경우(검정통계량이 임계값의 범위 이내인 경우) 예 측값에 대해 신뢰할 수 있는 것으로 정의하였다.
섬진강A 지점은 모니터링 기간 중 토양유실이 발생한 자료의 개수는 20개로 한국형 토양유실공식에 대한 예측 정확성 향상에 중점을 두었다. 한국형 토양유실공식에 의한 토양유실 가능성 추정치는 강우량에 따라 변화하였으 며, 이 추정치에 식(3)을 적용하여 토양유실량 예측치가 산정이 되었으며 NSE가 0.8020, 결정계수가 0.8134, p-value가 0.8889로 토양유실량 실측값과 비교할 때 예측값 을 신뢰할 수 있는 것으로 나타났다(Table 5, Fig. 3). 섬진강A 지점의 경우 토양유실량이 2014년 7월 21일에 가장 높게 측정되었으며, 대상 지점에서 모니터링이 더 수 행이 되어 이 토양유실량과 비슷한 범위를 가질 경우 보 정계수에 대한 검정이 가능할 것으로 판단된다.
섬진강B 지점은 모니터링 기간 내 토양유실량이 발생 한 경우가 모형의 보정 및 검정 과정이 수행되기에 충분 하다고 판단되어 2014년 7월 18일부터 2016년 4월 26일까지 22개 자료에 대해 보정계수에 대한 판단을 통해 식(4)를 유도하였으며, 이 보정과정에서 NSE는 0.5089, 결 정계수는 0.5576, p-value는 0.9731로 한국형 토양유실공식에 대한 보정이 잘 이루어진 것으로 판단하였다(Table 6, Fig. 4). 검정과정에서는 식(4)에 대한 조정 없이 보정 과정에서 정의된 계수를 그대로 적용하여 예측값과 실측 값에 대한 비교가 이루어졌으며, NSE가 0.5157, 결정계수 가 0.5999, p-value가 0.7992로 한국형 토양유실공식 및 식(4)에 의한 토양유실량 예측에 대해 신뢰할 수 있는 것으로 판단되었다.
Table 6.
Comparison between the measured and the calibrated estimated soil erosion (Mg) in Seomjingang B site during 2014~2016
Table 7.
Comparison between the measured and the validated estimated soil erosion (Mg) in Seomjingang B site during 2016~2018
섬진강A와 섬진강B 지점 모두 모니터링 기간은 2014년부터 2018년까지이며 현 고시에서 제시하고 있는 범용토양유실공식에 의한 예측값의 신뢰 가능여부를 판단하기에는 부족하다. 그러나 연단위가 아닌 월단위 예측이 가능한 한국형 토양유실공식의 경우 같은 기간에 대해서 다수의 실측값 및 예측값이 존재하기 때문에 모형에 대한 통계치 정의가 가능하였다.
4. 요약 및 결론
본 연구는 현재 한국환경공단에서 수행하고 있는 ‘표토 의 침식 현황 조사에 관한 고시’의 현장조사 단계에서 토 양유실량 실측 및 예측에 대한 비교 및 검정 과정이 이 루어질 수 있는 방안을 제시하는 데 목적을 두고 있다. 현재 고시는 전국을 대상으로 하여 토양유실량에 대한 예 측을 수행하여 현장 방문 지점을 선별하는 예비조사와, 대상 지점(혹은 필지)을 방문하여 자료를 수집하여 토양유 실량을 보다 정확하게 측정하는 현장조사로 구분이 된다. 현 고시에서는 이 두 조사 모두 범용토양유실공식을 사용 하도록 제시하고 있는데, 전국을 대상으로 하는 예비조사 의 경우 범용토양유실공식 적용은 타당한 것으로 판단된 다. 그러나 현장을 방문하여 자료를 수집하고 대상 지점 에 대한 토양유실량 예측 정확성 향상을 위한 현장조사단계에서는 이 범용토양유실공식의 적용에 한계점이 존재한다.
첫째, 범용토양유실공식은 연평균 혹은 연단위이기 때 문에 한 개의 실측 토양유실량을 적용하기 위해서는 1년 의 기간이 필요하며 모형에 의한 예측값과의 비교를 위해서는 수십 년 동안의 자료가 수집되어야 가능하다. 두 번째 한계점으로는 현 고시의 강우침식능인자가 현재 기상 조건을 반영하기에는 어렵다는 것이다. 현 고시의 강우침 식능인자는 우리나라의 1973년부터 1996년까지의 기상자 료에 근거하여 산정되었기 때문에 현재의 기상조건을 반 영할 수 있다고 보기 어렵다(Kim et al., 2018; Kongju National University, 2018). 즉, 기상조건에 의한 토양유 실 발생 가능성이 높았던 지점이 현재에도 높다고 단정하 기에는 어렵기 때문에, 현재 혹은 최근 기상 조건을 반영 해야 토양유실량에 대한 예측이 합리적일 수 있다. 세 번 째 한계점은 범용토양유실공식에 의해 예측은 토양유실 발생 가능성에 대한 추정치이기 때문에 실측 토양유실량 과 차이가 발생하지만, 이에 대한 보정은 한계가 있다. 발생한 토양 입자가 측정 지점까지 도달할 때에는 대상 지 점의 조건에 따라 달라질 수 있으나 이를 연간 특성으로 보기에는 어렵거나 예측 과정에서 불확실성이 커질 수 있기 때문이다.
따라서 현장조사 과정에 대한 개선이 필요하며, 이는 연 단위가 아닌 월단위 예측이 가능한 한국형 토양유실공식의 적용이 필요한 것으로 판단된다. 범용토양유실공식을 개선하여 국내 여건에 부합하는 데이터베이스를 제공하고, 월 강수량에 의해 임의 기간에 대한 강우침식능인자 산정 이 가능한 한국형 토양유실공식이 고시의 현장조사 단계 에 이용이 될 경우, 수년 동안 측정된 토양유실량과의 비교가 가능하기 때문에 앞서 언급된 범용토양유실공식의 첫 번째 한계점은 해결될 수 있다. 또한 현 고시와 같이 특정 기간 동안 수집된 자료에 의해 산정된 강우침식능인 자가 아니라, 임의 기간에 대한 강우량에 대한 강우침식 능인자 산정이 가능하기 때문에 두 번째 한계점 역시 해 결이 가능하다. 그리고 실측된 토양유실량과 비교할 수 있 도록 본 연구와 같이 보정계수 산정 과정이 이루어진다면 마지막 한계점 역시 해결될 수 있을 것으로 판단된다.
