조사는 전라북도 완주군 농생명로 166에 위치한 국립농업과학원에서 실시하였다. 토종꿀벌 사육장과 서양꿀벌 사육장, 그리고 서양꿀벌 사육장과 인접한 곳에 위치한 산림 세 지점에서 조사를 실시 하였다. 2009년과 비교하여, 서양꿀벌 사육장의 경우, 지형적 변화가 있었지만, 산림부나 토종꿀벌 사육장의 경우, 10년 이상 존재하던 산림이 자리잡고 있다.

조사는 2018년 5월1일부터 11월 6일까지, 2019년 4월 17일부터 11월 7일까지, 2020년 4월 6일부터 9월 21일까지 매주 1회씩 총 3년에 걸쳐 수행하였다. 태풍과 같은 천재지변의 영향을 받은 경우에는 조사 날짜를 조정하였다.

양봉꿀벌 사육장, 산림부, 토종꿀벌 사육장 각각에 말벌트랩을 6개씩 배치하여, 말벌 발생량을 모니터링 하였다. 말펄 포획을 위하여 유인액을 사용하였으며, 유인제는 조사지점당 트랩 6개 중, 3개씩 나누어 유인 물질을 달리하여 설치하였다. 유인제의 성상은 2018년에는 상업용 유인제A (장수말벌 유인용 500ml)와 꿀벌 페로몬 유인액을 사용하였으며, 2019년에는 유인제 A(장수말벌 유인용 500ml)와 포도주+가루식초(포도주 500ml+가루식초 약 20g) 를 사용하였다. 2020년에는 상업용 유인제B (등검은말벌 유인용 500ml)와 포도주+가루식초(포도주500ml+가루식초 약 20g)을 사용하였다.

환경조사를 위하여 기상청 기상자료개방포털의 종관기상관측(ASOS)서비스[11]를 이용하였다. 조사기간 동안의 평균기온(°C), 평균상대습도(%), 월강수량합(mm), 평균풍속(m/s)를 확인하였다. 관측지점은 전주관측소(N35.84092, E127.11718, Alt 60m)의 자료를 이용하였다. 인터넷 위성사진을 활용하여 각 조사지역 반경 100m를 격자(10m x 10m)로 나누어, 산림의 비중이 어느 정도를 차지하는지 분석하였다.

각 지점별 트랩에 채집된 말벌류의 군집분석을 하기 위하여 다음의 생물지수 분석방법을 이용하였다.

1) 우점도지수(Dominant index) - 제 1우점종과 제 2 우점종을 선정하여, 우점도지수를 산출(McNaughton, 1967)

2) 다양도지수(Diversity index) - Margalef(1958)의 정보 이론에 의하여 유도된 Shannon-Weaver function (H’)을 Lloyd & Gheraldi 가 변형한 공식을 사용(Pielou, 1975)

3) 종풍부도 (Richness index) - Margalef(1958)의 지수를 이용

유인물질별, 서식지별 말벌류의 포획률 평균의 차이는 Oneway ANOVA test 를 이용하였으며, 사후 검정은 DMART(Duncan’s Multiple Range test)를 이용하였다(SAS 9.4 Inc., 2008).

실험기간 내 환경으로 2018년(4월- 11월) 평균 기온 21.1°C, 평균 상대습도 78%, 평균 월 강수량 160mm, 평균풍속 1.68m/s, 2019년(4월- 11월) 평균 기온 20.7°C, 상대습도 80%, 평균 월 강수량 113mm, 평균풍속 1.54m/s, 2020년(4월- 9월) 평균 기온 20.9°C, 평균 상대습도 73%, 평균 월강수량 256.8mm, 풍속 1.65m/s로 조사되었다. 조사지점 반경 100m 내의 경관 요소 중, 산림의 비율을 비교한 결과, 산림부(Fig.1B)의 산림비율 100%, 토종꿀벌 사육장(Fig.1C)의 산림비율 59%, 서양꿀벌 사육장(Fig.1B)의 산림비율 41%로 확인되었으며, 산림부는 소나무와 밤나무가 혼재되어 있고, 토종꿀벌 사육장은 밤나무를 포함한 참나무류가 우점하고 있는 혼효림으로 확인되었다. 서양꿀벌 사육장의 경우, 서양꿀벌의 밀원수 식재에 의한 밤나무, 아까시나무, 헛개나무 등의 다양한 밀원수가 식재되어 있는 것을 확인하였다.

꿀벌 사육장 주변에 대한 3개 지점을 대상으로 총 3차년도 (2018, 2019, 2020)에 걸쳐 조사를 실시하였으며, 조사 결과, 3년동안 총 2561마리의 말벌류를 채집하였으며, 매년 평균 854마리가 채집되었다. 연도별로는 2018년에 7종 494개체, 2019년에 7종 1272개체, 2020년에 5종 795종으로 2019년에 가장 많은 개체가 포획되었다.

3년동안 채집된 개체들 중 주요종으로는 등검은말벌(Vespa velutina nigrithorax)(28%), 좀말벌(Vespa analis parallela)(23%) 순으로 우점하여 출현하는 것으로 확인되었다. 연도별로는 2018년에 꼬마장수말벌(44%), 좀말벌(25%), 2019년에 등검은말벌(42%), 장수말벌(16%), 2020년에 좀말벌(39%), 등검은말벌(20%) 순으로 우점하는 것으로 확인되었다(Table.1).

말벌류의 발생 시기로는 꼬마장수말벌의 경우, 최초 발생은 5월 말경에 확인이 되었으며, 산림 비율이 높은 토종꿀벌 사육장과 산림부에서 주로 발생하는 것을 확인하였으며, 서양꿀벌 사육장에서 현저하게 적게 발생하였다. 2018년에는 꼬마장수말벌이 8월 중순에서 9월 초까지 최성기를 나타내는 것을 알 수 있었다.

좀말벌의 경우, 4월초에 최초 발생하여 소수의 개체들이 확인되다가, 8월 중순부터 개체수가 많아지기 시작하여, 9월중순에 최성기를 이루고 10월까지도 관찰이 되었으며, 토종꿀벌 사육장에서 가장 많은 개체가 확인되었다.

말벌의 경우, 4월 말에 최초 발생하여 꾸준히 발생이 되다가 9월 초가 되면 최성기를 보이며, 주로 양봉장 인근 산림부에서 가장 많은 개체가 확인되었다.

장수말벌의 경우, 5월 초중순에 최초 발생하여 드물에 나타나다가, 8월 중순부터 개체수가 많아지다가 9월 초와 중순 사이에 최성기를 이루고 10월 초까지 활발하게 나타나며, 주로 토종꿀벌 사육장에서 가장 많은 개체가 확인되었다.

유인물질별 트랩을 비교분석한 결과, 2018년 유인물질로 사용한 유인제A와 꿀벌페로몬유인제의 경우, 꼬마장수말벌, 털보말벌, 검정말벌이 꿀벌 페로몬유인제에서 더 많은 개체가 포획되었으며, 좀말벌, 말벌, 장수말벌, 등검은말벌은 유인제 A에서 더 많이 포획되는 것을 확인하였다. 통계적으로 유의미한 결과는 얻을 수 없었다.

2019년에 유인물질로 사용한 유인제A와 포도주+가루식초의 경우, 통계적으로 유의미하지는 않지만, 유인제A에서 더 많은 개체들이 확인되었으며, 2020년에 유인물질로 사용한 유인제B와 포도주+가루식초의 경우, 통계적으로 유의미하지는 않지만, 유인제B에서 더 많은 개체들이 확인되는것을 확인하였다.

등검은말벌의 경우, 외래유입종으로 2003년 부산에서 처음으로 확인[6]된 이후, 매년 평균 9.4km씩 전국으로 확산되고 있으며[12] 꿀벌에게 심각한 피해를 주는 종이기 때문에 본 연구에서 중점적으로 분석하였다. 등검은말벌은 4월 초부터 6월 중순까지 여왕벌이 출현하면서 발생하였으며, 6월 말까지 드물게 출현하다가 9월 초에서 중순경부터 발생량이 증가하기 시작하여 10월 중순경에 최성기를 이루는 것을 확인하였다. 장소별로 확인한 결과, 여왕이 나오는 시기인 4월부터 6월까지는 서양꿀벌 사육장에서 다른 지점보다 낮은 빈도로 출현하였다. 8월이 되면서 서양꿀벌 사육장으로 출현하는 빈도가 높아지다가 9월 중순에서 10월 중순에 최성기를 보일 때는 토종꿀벌 사육장이나 산림부 보다 높은 출현 빈도를 보이는 것을 확인하였다. 또한 같은 시기에 인접부인 산림부와 유사한 피크를 보이는 것을 확인할 수 있었다.