서 론

El Niño는 적도의 태평양 온도 변화를 나타내며 이에 따라 다양한 지역의 온도와 강수에 영향을 미친다(Trenberth, 1997; Naylor et al., 2001; Barrios-perez et al., 2021). 이로 인해 El Niño가 벼의 수량에 영향을 줄 것으로 예상되지만(Zubair, 2002; Deng et al., 2010; Xie & Fang, 2020), Shim et al. (2016)의 연구에 따르면 El Niño 발생 시작 연도, 발생 종료 연도, 미 발생 연도로 분류하여 전국 또는 지역별 쌀 수량을 비교한 결과 통계적으로 유의미한 관계가 없다고 보고되었다. 그러나 2015-2016, 1997-1998, 1982-1983과 같이 역사적으로 강한 El Niño가 발생한 시기에는 쌀 수량이 상당히 높았다(https://kosis.kr/). 예를 들어, 2015년은 통계청 공식 기록이 시작된 1965년 이후 가장 높은 수량(542 kg/10a)을 기록했으며, 1997년은 1990년대 최고 수량(518 kg/10a)을 기록했다. 1982년의 경우 수량(438 kg/10a)은 1980년대 최고 수량은 아니었지만 전년보다는 높았다(https://kosis.kr/). 또한 2023/24년 역시 강한 El Niño가 될 것으로 예측되고 있는데 2023년 수량(523 kg/10a) 역시 전년(518 kg/10a)보다 높아 수량이 낮은 해라고 볼 수는 없어 El Niño가 쌀수량에 영향을 미친 사례로 추가되었다. 즉 역사적으로 강한 El Niño가 발생한 년도에는 쌀수량이 증가하는 특성을 보였다. 이는 벼 생육기간의 El Niño와 쌀 생산량을 변화에 대해서 다시 한번 분석할 필요성을 보여주고 있다. 따라서, 발생 시작 연도의 El Niño, 특히 벼 생육기간 동안의 El Niño와 쌀 생산량 간의 관계를 다시 분석함으로써 El Niño 발생을 활용하여 쌀 수량을 예측하는 데 도움이 될 수 있음을 확인하고자 본 연구를 수행하였다.

재료 및 방법

국내기상 및 El Niño

국내 기상자료는 기상청 기상자료개방포털(data.kma.go.kr)에서 1965년부터 2023년까지의 기간 동안 관측된 13개 지역의 종관기상관측 자료를 활용하였다. 이 자료에는 평균온도, 일조시간, 강수량 등이 포함되어 있다. 이 중에서 벼의 수량분석에 필수적인 일조시간이 지속적으로 기록된 지점 중, 1965년부터 벼의 생산 자료가 공개된 13개 지역이 선택되었다(Table 1).

El Niño 자료는 미 해양대기청(NOAA, National Oceanic and Atmospheric Administration)의 기후관측자료(https://www.ncei.noaa.gov/access/monitoring/enso/sst) 중 ONI (Oceanic Niño Index)에 따른 정의를 활용하였다. ONI에서는 Niño 3,4 지역 태평양 해안의 표면 온도의 평년차 3개월 이동평균이 0.5 이상인 경우를 El Niño로 정의한다. ONI 자료는 기상자료의 기간과 동일한 1965년 자료부터 사용하였으며, 벼의 작기 중 발생한 El Niño 현상이 이번 연구의 범위이므로 ONI 지수는 5․6․7월 이동평균값인 MJJ부터 8․9․10월 이동평균값인 ASO까지 모두 0.5 이상인 년도를 선정하였다. 해당 년도는 1965, 1972, 1982, 1987, 1991, 1997, 2002, 2015, 2023년이다. 이러한 분류는 Shim et al. (2016)의 El Niño 발생 년도에 따른 쌀 생산량을 분석한 결과를 참고하였다.

쌀수량 자료 및 추세 제거

쌀수량 자료는 통계청 국가통계포털(https://kosis.kr/)에서 획득하였으며, 통계청 공식 통계자료가 존재하는 1965년부터 2023년까지의 전국평균수량을 사용하였다. 시계열상의 쌀수량은 기상에 의한 변동 뿐만 아니라 품종 및 재배기술의 발전에 따른 수량증가 추세가 반영되어 있기 때문에 이를 기상과 수량과의 관계를 명확히 하기 위해서는 추세를 제거해야 한다(Kim et al., 2019). 추세 제거를 위해 5년 이동평균을 이용하였으며, 추세가 제거된 수량은 5년 이동평균값 대비 증감율(%)로 나타내었다(Fig. 1). 이동평균 중 5년을 이용하는 이유는 현재 농촌진흥청 작황보고(Rural Development Administration, 2022)와 한국농촌경제연구원 농업관측 보고(Rural Development Administration Agricultural Observatory Center, 2023)에서 5개년 이동 평균을 활용하기 때문에 동일한 기준을 적용하였다. 다만 농촌진흥청 작황보고는 5개년 평균을 한국농촌경제원 농업관측에서는 최대 최소값 제외한 평균을 활용하는데, 본 연구에서는 농촌진흥청의 5개년 평균값을 적용하였다. 1965년부터 1969년까지의 자료는 추세제거 과정에서 제외되었기 때문에 해당연도의 기상자료와 El Niño 자료 역시 차후 분석에서 제외하였다.

i : 해당연도

MA_i-1 : i-5부터 i-1까지의 수량평균

Yield_i : i 년도의 수량

수량 자료 이외에도 쌀 생산의 감소를 설명하기 위해 재해연보 자료(Ministry of Construction, 1982; Ministry of Construction, 1987; Ministry of the Interior, 1991; Ministry of Public Administration and Security Central Disaster Response Headquarters, 1997; Ministry of Public Administration and Security Central Disaster Response Headquarters, 2002; Ministry of Public Safety and Security, 2015)를 사용하였으며 재해의 지표로는 농작물 피해 면적 중 답(畓)의 면적만을 이용하였다. 다만 1972년은 재해연보 발간이전으로 한국의 홍수 자료를 대신 사용하였다(Ministry of Construction, 1972). 이 자료는 피해면적이 아닌 피해량으로 통계가 기록되어 있어, 1972년 평균 단수를 이용하여 면적을 역산한 값을 사용하였다. 또한 1991년 이전 자료는 농작물 피해면적 단위가 ha가 아닌 정보였으나, 1 ha와 1정보는 같은 단위로 가정하여 처리하였다. 또한 2023년은 공식 자료가 나오지 않아 우선 제외하였다.

Fig. 1.

National mean yield of milled rice (a) and detrended milled rice yield (b) in South Korea.

통계처리

ONI 지수와 추세 제거 수량 간의 관계를 규명하고, 이를 작황 예측의 활용을 위해 선형 회귀 분석을 수행하였다. 이 과정에서는 R.4.2.1 (Venables et al., 2022)의 l m 함수를 이용하였다.

결과 및 고찰

쌀수량의 추세제거

쌀수량의 전체적 수량 변화 추이는 Fig. 1a와 같이 통계가 작성된 1965년 이래 지속적으로 증가하고 있다. 이러한 추세는 Kim et al. (2019)에 따르면 결국 품종과 재배기술의 발전에 기인한 것이기 때문에 5년 이동 평균을 이용하여 추세를 제거한 결과는 Fig. 1b와 같다.

우리나라의 쌀수량의 진폭은 1980년 이후로는 비교적 일정한 것으로 보이며 수량이 음의 극값을 보이는 연도는 1980, 1993, 2003, 2012, 2020년이었으며 반대로 양의 극값을 보이는 연도는 1984, 1996, 2008, 2015년이었다. 음의 극값은 약 13년의 주기에서 8년의 주기로 주기가 짧아지고 있으며, 양의 극값은 대략 12년 주기에서 7년 주기로 역시 짧아지고 있다. 음의 극값의 주기는 1980년과 1993년 그리고 2003년이 있는 것으로 볼 때 냉해와 관련 있는 것으로 보이며, 양의 극값을 나타내는 시점의 기상학적 원인은 일조시간의 증가에 따른 것으로 보인다(Kim et al., 2017). 2015년을 제외한 위의 극값이 발생한 연도들은, 벼 작기 중 ONI가 지속적으로 0.5이상이었던 1965, 1972, 1982, 1987, 1991, 1997, 2002, 2015, 2023년과는 일치하지 않아 본 연구에서는 활용하지 않았지만, 수량 변동의 주기와 기상학적 요소와 관련해서는 이 부분에 대한 추가적인 연구가 필요할 것으로 보인다.

추세가 제거된 쌀 수량과 ONI와의 관계

추세가 제거된 쌀 수량은 결국 기상에 의한 벼의 생육 반응이며, 그에 따른 쌀수량의 변동이라고 할 수 있다. 이 변동에 대한 El Niño와의 관계를 각 시기별 ONI와 추세가 제거된 수량 간의 관계로 나타내었다(Fig. 2). ONI가 커질수록 추세가 제거된 수량이 증가하는 경향은 있었다. 각 시기별 결정 계수는 MJJ (May, Jun and July) 에서 0.3639로 가장 낮고 ASO (August, September and October)에서 가장 높아 0.6838이어 ONI와 수량과의 일정한 관계가 있을 수 있음을 시사하였으나 작황을 예측하기에는 설명력이 부족하였다.

Fig. 2.

Relationship between Oceanic Niño Index (ONI) and detrended rice yield. a: ONI of 3-month moving average for May, June and July. b: ONI of 3-month moving average for Jun, July and August, c: ONI of 3-month moving average for July, August and September, d: ONI of 3-month moving average for August, September and October.

또한, 8개의 발생 연도 중 3번은 수량이 감소하였고, 5번은 증가하여 명확한 일관성이 발견되지 않았다. 이를 분석하기 위해 수량이 감소한 연도와 증가한 연도의 기상을 분리하여 1991-2020년 평년 기상차를 비교하였다(Fig. 3).

벼 재배기간동안의 온도를 먼저 살펴보면, 결과적으로 추세 제거된 수량이 증가한 연도에서는 전반적으로 온도가 높았고, 특히 7월 하순과 8월 상순의 온도가 평년대비 1°C가량 높았다. 반면 추세 제거된 수량이 감소한 연도에서는 6월 중순 이후 지속적으로 온도가 낮아져 8월 상순에는 2.5°C까지 낮아진 것으로 나타났으며(Fig. 3a), 이는 지속적인 강우로 흐린 날씨가 지속되어 발생한 것으로 보인다. 그에 대한 근거로 추세 제거된 수량이 감소한 연도들의 평균 강수량은 7월 중순부터 8월까지 대체로 증가하는 반면, 증가하였던 연도는 강수량이 감소하였다(Fig. 3b). 일조시간의 경우에는 수량이 감소한 연도에서는 7월과 8월 일조시간이 감소하였고, 8월 초순에는 2.5시간 이상 감소하였다(Fig. 3C). 그러나 수량이 증가한 연도에서는 일조시간이 7월부터 9월까지 평년보다 높게 유지되었다(Fig. 3C). Kim et al. (2014)은 우리나라에서 벼가 풍년이 되기 위해서는 7월의 일조와 8월 하순부터 9월 중순까지의 일조시간이 동시에 높아야 한다고 하였다. 본 연구에서도 수량이 증가한 연도에서는 7월부터 9월까지의 일조가 높아 추세 제거 수량이 양의 값을 보였던 반면, 수량이 낮은 연도에서는 반대로 7월과 8월 일조가 낮은 경향을 보였다(Fig. 3C). 한편, 태풍 및 호우 발생에 의한 차이 역시 수량의 증감을 결정하는 것으로 보인다. 본 연구에서 선정된 8개 El Niño 발생 연도 중 추세 제거 수량이 음의 값을 보였던 3개년과 양의 값을 보였던 5개년에 대해서 각각 태풍 피해에 따른 농작물(답)의 피해 면적을 조사한 결과 음수였던 해의 피해가 더 큰 것을 알 수 있었다(Fig. 2; Fig. 4).

Fig. 3.

The averaged meteorological factor for years with negative and positive detrended yields. Meteorological factors: mean temperature (a), precipitation (b) and sunshine hours (c).

Fig. 4.

Box plot of damaged rice area in years with negative and positive yields.

연도별로 보면 2002년이 1987년과 1991년에 비해 높았는데, 2002년은 태풍 RUSA, 1987년은 태풍 THELMA로 2014년 국민안전처(Ministry of Public Safety and Security, 2014)에서 정리한 우리나라 기상재해 피해액 순위에서 1위와 9위를 차지할 만큼 큰 기록적인 태풍이 발생한 해였다. 1991년 역시 태풍 GLADYS에 의한 피해가 발생하여 16위에 위치했다. 반면 추세 제거 수량이 양의 값을 보였던 해에는 태풍이 발생했으나 큰 피해가 없었다(Table 2).

이는 El Niño의 발생에 따라 태풍의 크기, 태풍이 오는 시기, 태풍의 이동 경로, 풍속 등이 변화되어 벼의 수량 감소 정도가 달라진 것으로 볼 수 있었다. 이에 대해서는 차후 추가적인 연구가 필요할 것으로 보인다. 종합하면 추세 제거 수량이 양의 값을 보이는 해에는 7월과 8월 여름 및 초가을 9월의 강우량이 적고, 가을인 10월 강우가 많아 벼의 출수 및 등숙 초기에는 일조가 많은 특성을 보이며, 또한 태풍이 오더라도 벼에는 적은 피해를 입혔다고 할 수 있다(Fig. 2; Fig. 3; Fig. 4; Table 2).

반면에 추세 제거 수량이 음의 값을 보이는 해는 7월과 8월 여름의 강수량이 높아, 벼의 출수 및 등숙 초기 일조가 적고, 반대로 9월과 10월 가을에는 강수량이 낮고 일조시간이 높았으며, 태풍의 영향이 크다고 할 수 있었다(Fig. 2; Fig. 3; Fig. 4; Table 2). 이와 관련하여 El Niño의 2가지 종류에 대해서 검토해 볼 필요가 있을 것으로 생각되었다.

El Niño 종류에 따른 쌀수량의 증감

Kug et al. (2009)는 El Niño를 Warm pool El Niño와 Cold tongue El Niño로 분류하였다. 1972/73, 1976/77, 1982/83과 1997/98을 Cold tongue El Niño으로, 1977/78, 1990/91, 1994/95, 2002/03과 2004/05를 Warm pool El Niño로 분류하였으며, 1986/87, 1987/88,과 1991/92은 이 두 종류의 중간 정도로 분류하였다(Kug et al., 2010). Kug et al. (2009)의 정의대로라면 2015년은 Cold tongue El Niño, 2023년 역시 Cold tongue El Niño로 분류할 수 있다.

Kug et al. (2009)의 정의는 간단히 설명하면, El Niño가 발생한 연도의 9월부터 그 다음 해 2월까지의 SST (Sea surface temperature)에 대한 과거 35년 평년값에 대해 평균을 사용하며, 그 결과 SST의 이동평균이 평년값 대비 약 0.5 이상의 고온을 보이는 현상의 중심이 Nino 3에 있으면 Cold tongue El Niño으로 분류한다고 한다. 그런데 이들 Cold tongue El Niño인 해에서는 대부분 추세 제거 수량이 양의 값을 보였다는 점이 확인되었다(Fig. 2). 1972, 1982, 1997, 2015, 2023년이 본 연구의 양의 값을 보인 연도이자 Cold tongue인 연도였다(Fig. 2). 한편 SST의 이상고온 현상의 중심이 Nino 4에 있으면 Warm pool El Niño로 분류한다. Warm pool El Niño에서는 Cold tongue El Niño와 반대로 2002년이 음의 수량을 보이는 연도와 동일하였고, Kug et al. (2010)의 분류에서 중간적인 특성을 보인 연도인 1987년, 1991년이 음의 수량을 보이는 연도와 동일하였다. 즉, Cold tongue El Niño인 연도와 그렇지 않은 연도로 나눌 경우, 양의 수량을 보이는 연도와 음의 수량을 보이는 연도가 구분될 수 있다는 것이다.

두 종류의 El Niño 중 Warm pool El Niño인 경우에는 강우량과 관련해서 명확한 관련성이 없다고 하였다(Kug et al., 2010). 일반적으로 El Niño는 발생 시기 이전 연도 6월부터를 발달기로 정의하고 이듬해 쇠태기를 7월로 정의하는데, 이를 적용하면 벼의 생육기는 발달기 6월부터로 볼 수 있다. Kug et al. (2010)의 결과에 따르면 Cold tongue El Niño시에는 강우량에 유의한 증가가 있다고 하였다. 11월의 강수량이 높고 9월의 강수량이 낮았는데 보통 강수량과 일조시간이 반대인 것을 생각하면 11월 일조가 낮고 9월 일조가 높을 가능성이 높다(Kug et al., 2010). 이는 벼의 중후기 등숙에 유리한 조건을 제공할 수 있다. 앞서 정리한 추세 제거한 수량이 양의 값을 보였던 연도의 특성과 일정부분 일치한다(Fig. 2). 만일 El Niño 종류에 따른 분류를 Fig. 2에 적용하여, 음의 수량값을 보이는 연도를 Cold tongue El Niño가 아닌 것으로 분류하여 제거하면 Fig. 5와 같이 좀 더 높은 설명력을 보이는 관계식을 얻을 수는 있다. El Niño가 자주 발생하는 현상이 아니며 벼의 작기중에 지속적으로 ONI 지수가 0.5를 넘는 경우는 많이 발생하지 않기 때문에 단정지을 수는 없으나, 벼의 작기중 발생한 El Niño 즉 6월 발달기의 El Niño가 Cold tongue El Niño로 판정된다면 작황 예측에도 활용이 가능할 것으로 보인다. 만일 작황예측에 활용된다면 El Niño의 종류에 대해 판별한 후 5, 6, 7월 이동평균 값이 나오는 8월 또는 6, 7, 8월 이동평균 값이 산출되는 9월에 예측이 가능할 수 있을 것이다.

그러나, 이를 위해서는 향후 작기 중 발생한 El Niño 종류에 대한 정의와 판별이 우선적으로 필요하다. 본 연구에서는 Kug et al. (2009)과 Kug et al. (2010)을 인용하였으나, 9월이 후의 자료를 활용하여 벼의 작기 중에 판별이 용이하지 않는 문제가 있다. 따라서 이를 판별하기 위한 다른 방법이 필요하다. 그런데 El Niño의 종류에 대한 정의 방법이 연구자 마다 다르기 때문에(Kim et al., 2014), 각각의 정의 방법별로 추세 제거 수량이 양이 나오는 연도와 음이 나오는 연도에 대해서 검토하고, 이때의 기상자료를 벼의 생육 단계에 따라 세부적으로 분석할 필요가 있다. 또한 향후 발생할 El Niño자료를 지속적으로 수집해서 지속적으로 검증을 해야 할 것이다.

Fig. 5.

The relationship between Oceanic Niño Index (ONI) and *** (DR), excluding the three years in which yields decreased due to typhoons. a: ONI of 3-month moving average for May, Jun and July. b: ONI of 3-month moving average for Jun, July and August, c: ONI of 3-month moving average for July, August and September, d: ONI of 3-month moving average for August, September and October.

적 요

1.2015년 이후 가장 강한 El Niño가 2023년 발생하여 El Niño가 농업부분에 미치는 영향에 대해서 관심이 높아지고 있으며, 특히 우리나라의 쌀 수량에 미치는 영향을 다시 검토할 필요가 있다. 벼의 작기중에 발생한El Niño가 쌀 수량의 증감에 영향을 주는지 규명하기 위해, 추세가 제거된 쌀 수량 자료와 벼 작기중ONI 값이 0.5이상으로 유지된 연도의 자료를 이용하여 이들의 관계를 조사하였다. 그 결과 ONI 강도가 높을수록 추세가 제거된 쌀 수량이 높게 나오는 경향을 확인하였다.

2.또한 추세 제거 쌀 수량 지수가 양의 값으로 나오는 연도와 음의 값으로 나오는 연도에는 여름 및 가을의 일조시간 및 강수량이 다를 뿐만 아니라 태풍의 영향이 다르다는 것을 확인하였다. 이는 El Niño를 Cold tongue El Niño인 경우와 그렇지 않은 경우로 분류할 때와 유사하였다. Cold tongue El Niño인 경우에는 추세 제거 수량이 양의 값을 보였으며, Cold tongue El Niño이 아닌 경우에는 모두 음의 값을 보였다. 이를 활용하여 Cold tongue인 해에 ONI강도에 따라 수량의 변화를 일정수준 예측이 가능할 것으로 보인다.