2.1 Standardized Precipitation Index (SPI)

McKee et al. (1993, 1995)은 강수량이 감소하기 시작하면 물의 수요에 비해 상대적으로 물의 공급이 부족한 상태는 가뭄발생이 시작된다는 것을 착안하여 가뭄에 대한 크기를 분석할 수 있는 SPI를 개발하였다. SPI의 분석은 시간단위나 월단위를 기준으로 누가강수시계열을 산정하여 지속시간별 시계열을 분석한다. 지속시간별 시계열 강우량을 월별로 분석하여 변량의 누가확률을 산정하고 표준정규분포에 적용시켜 산정한다. SPI를 이용한 가뭄분석은 Gamma 확률밀도함수를 이용하여 해당 관측지점에 대한 강수량의 주기분포를 토대로 시간축별로 3개월, 6개월, 9개월, 12개월 등의 가뭄지수를 추정한다. 가뭄의 시·공간적인 변동성을 파악하기 위해 SPI를 활용하였으며 가뭄정도의 판정은 Table 1과 같다.

Table 1

Drought Classification by SPI Value and Corresponding Event Probability

SPI의 매개변수 추정은 최우도법(Maximum likelihood method)을 사용하였으며 Eq. (1)을 통해 계산된 매개변수를 이용하여 대상지점의 시간간격에 대한 강우사상의 누가확률을 분석한다.

여기서, α는 형상매개변수(Shape parameter), β는 축척매개변수(Scale parameter), x는 강우량이며 Gamma 확률밀도함수의 계수로 지정된 시간 축척별 따라 해당 관측지점별로 추정된다. 이에 대한 α, β의 추정치는 Eq. (2)~Eq. (3)에 의해서 계산할 수 있다.

그리고, n은 강수량 자료의 개수이며, 얻어진 매개변수들은 다음과 같은 누가확률분포 함수식에서 Eq. (4)와 같이 이용된다.

위의 식에서 t=x/β로 하면 Eq. (5)와 같이 정리된다.

Gamma 함수는 x = 0으로 정의되지 않으나 강수량이 0인 경우가 있으므로, 누가확률은 Eq. (6)과 같다.

q는 0인 경우의 확률로서 무강우 일수를 m, 강우관측일수를 n이라고 하면 무강우 확률 q는 q = m/n으로 나타낼 수 있다. 누가 확률 H(x)를 다시 평균 0 분산 1의 표준정규분포의 임의 변수 Z로 나타낼 수 있도록 변환하면 SPI 값이 된다(Eq. (7)~Eq. (10)).

여기서, C₀ =2.515517, C₁ =0.802853, C₂ =0.010328, d₁ =1.432788, d₂ =0.189269, d₃ =0.001308이며 x는 강수량이며 H(x)는 강수량 관측 값의 누가확률이며 C₀, C₁, C₂, d₁, d₂, d₃,는 상수 값이 된다.

2.2 농업가뭄의 피해현황

우리나라의 가뭄은 과거 약 10년을 주기로 큰 피해가 발생하였으며 2018년 ‘자연재해대책법 시행령, 2018. 06’의 개정으로 가뭄이 자연재해에 포함되어 과거 풍수해에 대한 피해가 대부분 이였다면 가뭄에 대한 체계적인 재난관리, 저감계획 등의 재해예방을 위한 제도의 개선과 피해현황에 대한 자료의 구축이 진행되고 있다. 2018년 이전의 가뭄의 피해현황은 지방정부에서 현황 파악용으로 기록된 피해자료가 전부이며 큰 가뭄피해 발생시 부처에서 ‘가뭄기록조사 보고서(1995, 2001)’, ‘이상기후보고서(2010-2018)’, ‘가뭄정보통계집(2019)’ 등의 보고서를 발간하였다. 본 연구에서는 우리나라에서 발생된 가뭄피해를 대상으로 1965년부터 2018년까지 발생한 가뭄에 대한 피해현황을 분석하고자 한다. 과거 발생한 가뭄의 피해는 대부분 농업가뭄으로서 논과 밭에 대한 피해만 제시되어 있을 뿐 생활가뭄이나 공업가뭄에 대한 피해자료는 구축되어 있지 않다. 우리나라에서 1965년부터 2018년까지 발생한 지자체별 농업가뭄 피해현황은 Fig. 1과 같다.

Fig. 1

Number of Agricultural Drought Damage Occurrences by Municipality (1965-2018)

대한민국에서 1965년부터 2018년까지 논과 밭에서 작물의 고사와 시듦 등의 피해가 발생된 농업가뭄은 총 32회가 발생하였다. 농업가뭄의 피해는 최소 1개에서 많게는 12개의 지자체에서 피해가 발생하였으며 울산광역시만 피해가 발생하지 않았다. 농업가뭄의 피해는 최소 1개에서 많게는 12개의 지방정부에서 피해가 발생하였으며 울산광역시는 농업가뭄의 피해가 발생하지 않았다. 지자체에서 발생한 농업가뭄의 피해발생 횟수는 울산광역시를 제외하고 최소 1회에서 많게는 20회의 피해가 발생하였다. 농업가뭄의 피해는 특별시도와 광역시에서는 최소 0회에서 최대 2회가 발생하였으며 8개의 시도에서 최소 7회에서 최대 20회가 발생하였다.

가뭄피해가 10회 이상 발생한 지자체는 강원도가 15회, 경기도가 12회, 경상남도가 20회, 경상북도가 18회, 전라남도가 18회, 전라북도가 12회, 충청남도가 16회, 충청북도가 14회로 특광역시를 제외한 지자체에서 대부분 농업가뭄의 피해가 발생하였다. 농업가뭄의 피해현황은 1965년부터 1982년까지는 매년 지속적인 피해가 발생하였다면 그 이후에는 1994년과 1995년, 2000년, 2001년, 2006년, 2013년 ~ 2018년에 가뭄피해가 발생하였다. 농업가뭄의 피해지역은 총 32회의 농업가뭄 피해기간 동안 145개의 지방정부에서 발생하였으며 평균적으로 4개 ~ 5개의 지역이 농업가뭄의 피해영향에 포함되었다.

3.1 대상지역 및 강우분석

본 연구에서는 과거 우리나라에서 1965년부터 2018년까지 발생한 농업가뭄피해를 고려한 가뭄지수를 분석하고자 한다. 대상지역은 17개 지자체 중에 총 농업가뭄 피해 32회 중 50%인 16회 이상이고 지속적인 가뭄피해가 5년 이상 지속된 경상남도를 대상지역으로 Fig. 2(a)와 같이 선정하였다. 우리나라는 지자체뿐만 아니라 다양한 부처에서 강우관측소를 운영하고 있으나. 설계기준이나 재해저감계획 등에서는 기상청에서 운영하고 30년 이상의 관측자료를 확보하고 있는 강우관측소를 사용하고 있다.

Fig. 2

Study Area and Raingauge

경상남도를 대상으로 강우량을 분석하기 위해서 경상남도 지역에 포함된 강우관측소와 인근지역에 위치한 강우관측소 15개를 선정하여 티센다각형을 작도하였다. 경상남도의 강우량을 산정하기 위해서는 다양한 설계기준에서 활용하고 있는 티센다각형 방법을 사용하였으며 경상남도와 강우관측소에 대한 티센다각형은 Fig. 2(b)와 같이 분석되었다. 경상남도의 강우량 산정은 티센다각형의 면적비를 가중치로 부여하여 평균 강우량을 산정하였다. 기상청에서 운영하고 있는 강우관측소의 관측기간은 1965년 ~ 2018년까지 일단위 강우량을 수집하였으며 15개 강우관측소의 지점번호, 관측소명, 티센가중치는 Table 2와 같다.

Table 2

Thiessen Weight by Rainfall Station

3.2 가뭄지수 분석

본 연구에서는 경상남도를 대상으로 15개 강우관측소에서 1965년부터 2018년까지 관측된 강우자료를 이용하여 지속기간 3개월, 6개월, 9개월, 12개월의 SPI를 분석하였다. SPI는 각 월별 지속시간에 따라 SPI 3은 가뭄지수가 최소 –3.23에서 최대 2.90, SPI 6은 최소 –2.83에서 최대 2.42, SPI 9는 최소 –2.36에서 최대 2.41, SPI 12는 최소 –2.60에서 최대 2.82가 Fig. 3과 같이 분석되었다.

Fig. 3

Analysis of SPI by Duration

SPI는 가뭄지수가 –1.0 이하일 때 가뭄으로 진입한다고 판단하며 –1.0 ~ -1.5는 약한가뭄, -1.5 ~ -2.0은 보통가뭄, -2.0 이하는 심한가뭄이 발생되는 것으로 판단한다. 경상남도의 SPI 3에서 –1.0 이하의 가뭄은 총 40회가 분석되었으며 –1.5 이상의 약한가뭄은 18회, -2.0 이상의 보통가뭄은 9회, -2.0 이하의 심한가뭄은 13회가 발생하였다. SPI 6에서 –1.0 이하의 가뭄은 총 29회가 분석되었으며 –1.5 이상의 약한가뭄은 12회, -2.0 이상의 보통가뭄은 9회, -2.0 이하의 심한가뭄은 8회가 발생하였다. SPI 9에서 –1.0 이하의 가뭄은 총 18회가 분석되었으며 –1.5 이상의 약한가뭄은 6회, -2.0 이상의 보통가뭄은 5회, -2.0 이하의 심한가뭄은 7회가 발생하였다. SPI 12에서 –1.0 이하의 가뭄은 총 11회가 분석되었으며 –1.5 이상의 약한가뭄은 3회, -2.0 이상의 보통가뭄은 4회, -2.0 이하의 심한가뭄은 4회가 발생하였다.

4.1 농업가뭄 피해특성을 고려한 가뭄지수 분석

경상남도에서 1965년부터 2018년까지 발생한 농업가뭄의 피해현황을 조사하였으며 각 가뭄피해에 대한 재해발생기간을 분석하였다(Fig. 4). 경상남도의 농업가뭄은 총 20회가 발생하였으며 대부분의 재해발생기간은 5월부터 9월 사이에 대부분의 피해가 발생하였으며 장마기간을 전후로 해갈되었다. 그 외 2000년은 2월부터 5월까지 봄 가뭄이 2008년은 9월부터 12월까지 가을 가뭄이 발생하였다. 가뭄의 지속기간은 1980년 이전에는 6월부터 9월까지 4개월이 1980년 이후에는 5월부터 8월 사이에 짧게는 1개월 길게는 2개월 동안 가뭄이 지속되었다.

Fig. 4

Agricultural Drought Damage Period in Gyeongsangnam-do (1965-2018)

우리나라의 농업가뭄피해에 대한 지속기간은 4개월 이하이며 대부분 여름철 장마전선의 영향으로 가뭄이 종료된다. 이러한 영향으로 가뭄의 피해가 대부분 단기가뭄의 형태로 발생되기 때문에 지속기간이 긴 SPI를 분석하면 SPI 3에서는 가뭄이 발생하였는데 SPI 6에서는 가뭄이 해갈되는 결과가 분석된다. 따라서, 본 연구에서는 SPI에 대한 가뭄피해기간을 적용할 때 가뭄지수가 0.0 이상의 경우 분석자료에서 제외하였다. 농업가뭄 피해발생기간에 대하여 지속시간별 SPI의 분석결과를 Fig. 5와 같이 산정하였다.

Fig. 5

Establishment of Drought Disaster Period considering Damage of Agricultural Drought

농업가뭄 피해기간의 SPI 3은 평균 –1.52로 약한가뭄이며 최대 –0.36에서 최소 –2.20, SPI 6은 평균 –1.65로 보통가뭄이며 최대 –0.65에서 최소 –2.83, SPI 9는 평균 –1.31로 약한가뭄이며 최대 –0.18에서 최소 –2.36, SPI 12는 평균 –1.06로 약한가뭄이며 최대 –0.07에서 최소 –2.45가 분석되었다(Table 3). 지속시간별 SPI의 분석결과 농업가뭄의 피해범위는 최소 –2.0 이하, 최대 –1.0 이상에서 발생하는 것은 비슷하였다. 그러나 지속시간별 SPI의 평균은 SPI 3에서 SPI 6은 증가하는 반면 SPI 9에서 SPI 12는 감소하는 가뭄지수를 나타내는 것으로 보아 우리나라 장마기간의 호우영향이 반영된 것으로 판단된다.

Table 3

SPI by Duration of Agricultural Drought Damage Period

4.2 재해예측을 위한 가뭄지수 기준 설정

지속시간별 SPI에 대하여 농업가뭄의 피해기간동안 발생된 연도별 최소 가뭄지수와 경상남도의 –1.0 이하로 약한가뭄이 발생되는 연도별 최소 가뭄지수에 대하여 상관계수(Coefficient of correlation)를 Fig. 6과 같이 분석하였다. SPI 3의 –1.0 이하로 약한가뭄이 발생되는 경상남도의 연도별 최소 가뭄지수 40개와 농업가뭄 피해가 발생된 연도의 최소 가뭄지수 20개의 상관계수는 0.749로 양의 상관성이 높은 것을 분석되었다. SPI 6의 –1.0 이하로 약한가뭄이 발생되는 경상남도의 연도별 최소 가뭄지수 29개와 농업가뭄 피해가 발생된 연도의 최소 가뭄지수 16개의 상관계수는 0.697로 양의 상관성이 높은 것을 분석되었다. SPI 9의 –1.0 이하로 약한가뭄이 발생되는 경상남도의 연도별 최소 가뭄지수 18개와 농업가뭄 피해가 발생된 연도의 최소 가뭄지수 15개의 상관계수는 0.645로 양의 상관성이 높은 것을 분석되었다. SPI 12의 –1.0 이하로 약한가뭄이 발생되는 경상남도의 연도별 최소 가뭄지수 12개와 농업가뭄 피해가 발생된 연도의 최소 가뭄지수 16개의 상관계수는 0.001로 무상관으로 분석되었다. 따라서, 우리나라에서 발생된 농업가뭄의 피해특성을 고려한 재해예측은 SPI 3, SPI 6, SPI 9는 영향이 높으며 SPI 12는 농업가뭄의 피해와는 영향이 미미한 것으로 분석되었다.

Fig. 6

Correlation Coefficient of SPI by Duration

경상남도에서 발생한 농업가뭄의 피해기간과 지속기간별 SPI의 상관관계는 높은 것으로 분석되었으나 재해예측을 위한 가뭄지수의 기준을 제안하기에는 어려움이 있다. 어떤 지속기간별 SPI가 농업가뭄의 피해발생을 예측하였는지 판단하기에는 지속기간별 SPI의 차이가 크게 발생하지 않는다. 따라서, 본 연구에서는 농업가뭄의 재해예측을 위한 기준을 과거 농업가뭄의 피해가 발생된 기간의 지속시간별(3개월, 6개월, 9개월, 12개월) SPI에 최소 가뭄지수를 연도별 대표 가뭄지수로 산정하여 통합적 기준을 Fig. 7과 같이 제안하고자 한다. 연도별 과거 농업가뭄의 지속시간별 최소 SPI는 피해횟수 총 20회 중에 SPI 3이 12회, SPI 6이 4회, SPI 9가 2회, SPI 12가 2회로 분석되었다. 농업가뭄의 재해예측을 위한 가뭄지수의 기준은 과거 피해가 발생한 연도의 지속기간별 SPI의 평균인 –1.64로 설정하였다.

Fig. 7

Establishment of Drought Index Criteria for Disaster Prediction