Assessing Risk and Economic Losses from Coconut Rhinoceros Beetle in the Americas

기술 요약: 아메리카 대륙 내 코코넛 rhinoceros beetle(CRB)의 위험성 및 경제적 손실 평가

문제 제기
코코넛 rhinoceros beetle(Oryctes rhinoceros, 이하 CRB)은 코코넛 및 오일 팜(oil palm)에 막대한 피해를 주는 파괴적인 해충으로, 역사적으로 아시아와 태평양 지역에서 상당한 생산량 손실을 초과해 왔다. 태평양 제도 전역으로 확산된 이후, 2025년 멕시코 푸에르토 바야르타(Puerto Vallarta)에서 정착된 개체군이 발견됨에 따라 아메리카 대륙에서도 처음으로 보고되었다. 멕시코 중부에서 남미 북부에 이르는 해안 및 저지대를 포함하는 신열대 아메리카 지역은 CRB의 생존에 적합한 기후 조건과 기주 식물(코코넛 및 오일 팜)을 보유하고 있다. 오일 팜 농장의 급격한 확장과 이 지역의 코코넛 생산에 대한 경제적 의존도를 고려할 때, 추가 확산 위험을 정량화하고 잠재적인 경제적 피해를 추정하여 생물 보안 전략을 수립하는 것이 시급하다.

방법론
본 연구는 아메리카 대륙 전역의 CRB 유입 및 정착에 대한 상대적 확률을 추정하기 위해 공간 명시적 위험 평가 모델을 사용한다. 모델은 5 arcminute(~9.25 km) 격자 해상도로 작동하며, 두 가지 주요 구성 요소인 **유입 압력(Introduction Pressure, $I$)**과 **정착 가능성(Establishment Potential, $E$)**을 곱하여 위험 점수($R$)를 산출한다.

1. 유입 압력 ($I$): 이 구성 요소는 자연적 분산을 제외한 인간 매개 운송 경로를 모델링한다. 이는 가중치가 부여된 세 가지 요인의 합이다.

   • 항만 근접성 (Port_Proximity, $P$): 거리와 항만 규모에 기반한 지수 감소 함수를 사용하여 해상 및 항공 경로로부터의 위험을 반영한다.
   • 인구 밀도 (Human_Density, $H$): 여행 및 우편과 관련된 위험을 나타내며, 로그 변환된 인구 밀도 데이터를 사용한다.
   • 도로 가중치 (Road_Weight, $D$): 상품 운송 위험을 나타내며, 도로 밀도에 에피센터 연결성(Epicenter_Connectivity) 요인을 곱하여 계산한다. 이 요인은 인간 이동성의 보편적 방문 법칙을 사용하며, 현재 발생 중심지(멕시코 나야리트)로부터의 거리의 역제곱에 따라 연결성이 감소한다.

2. 정착 가능성 ($E$): 이 구성 요소는 도입된 후 딱정벌레가 개체군을 형성할 가능성을 모델링한다. 이는 다음의 곱으로 이루어진다.

   • 팜 피복도 (Palm_Cover, $C$): 코코넛 및 오일 팜 밀도의 복합 지수(로그 변환 및 정규화됨).
   • 기후 적합성 (Climate_Suitability, $S$): GBIF(Global Biodiversity Information Facility)의 확인된 CRB 출현 기록에 맞춰 조정되고 WorldClim v2.1 생물 기후 변수를 대상으로 투영된 BIOCLIM 알고리즘(R 패키지 dismo 활용)을 통해 생성된 연속 표면.

팜 피복도가 0인 격자 셀은 위양성(false positives)을 방-지하기 위해 최종 위험 점수를 0으로 할당한다. 본 연구는 상세 분석을 위해 6개의 주요 팜 생산국(콜롬비아, 브라질, 멕시코, 온두라스, 과테말라, 에콰도르)에 집중한다.

주요 기여 및 결과

• 공간 위험 지도 작성: 모델은 습윤 열대 및 아열대 지역, 특히 높은 팜 피복도와 높은 인간 활동을 보이는 해안 및 저지대 지역에 집중된 고위험 구역을 식별한다. 고위험 구역은 중앙아메리카, 카리브해, 그리고 브라질의 넓은 지역을 포함한 남미 북부를 가로질러 예측된다. 북미의 온대 지역과 남미 최남단 지역은 위험도가 미미하다.
• 국가별 취약성:
  • 콜롬비아 및 에콰도르: 고위험 값은 해안 및 저지대 지역에 집중되어 있다.
    기존의 발생 중심지(나야리트, 멕시코)로부터의 거리의 역제곱에 따라 연결성이 감소한다.

3. 도로 가중치 (Road_Weight, $D$): 상품 운송 위험을 나타내며, 도로 밀도에 에피센터 연결성(Epicenter_Connectivity) 요인을 곱하여 계산한다. 이 요인은 인간 이동성의 보편적 방문 법칙을 사용하며, 현재 발생 중심지(멕시코 나야리트)로부터의 거리의 역제곱에 따라 연결성이 감소한다.

4. 정착 가능성 ($E$): 이 구성 요소는 도입된 후 딱정벌레가 개체군을 형성할 가능성을 모델링한다. 이는 다음의 곱으로 이루어진다.

   • 팜 피복도 (Palm_Cover, $C$): 코코넛 및 오일 팜 밀도의 복합 지수(로그 변환 및 정규화됨).
   • 기후 적합성 (Climate_Suitability, $S$): GBIF(Global Biodiversity Information Facility)의 확인된 CRB 출현 기록에 맞춰 조정되고 WorldClim v2.1 생물 기후 변수를 대상으로 투영된 BIOCLIM 알고리즘(R 패키지 dismo 활용)을 통해 생성된 연속 표면.

팜 피복도가 0인 격자 셀은 위양성(false positives)을 방지하기 위해 최종 위험 점수를 0으로 할당한다. 본 연구는 상세 분석을 위해 6개의 주요 팜 생산국(콜롬비아, 브라질, 멕시코, 온두라스, 과테말라, 에콰도르)에 집중한다.

주요 기여 및 결과

• 공간 위험 지도 작성: 모델은 습윤 열대 및 아열대 지역, 특히 높은 팜 피복도와 높은 인간 활동을 보이는 해안 및 저지대 지역에 집중된 고위험 구역을 식별한다. 고위험 구역은 중앙아메리카, 카리브해, 그리고 남미 북부를 가로질러 예측된다. 북미의 온대 지역과 남미 최남단 지역은 위험도가 미미하다.
• 국가별 취약성:
  • 콜롬비아 및 에콰도르: 고위험 값은 해안 및 저지대 지역에 집중되어 있다.
  • 브라질: 광범위한 저-중위험을 보이며, 북부 및 멕시코만 연안에 더 높은 위험 패치들이 나타난다.
  • 멕시코: 높아진 위험은 남부 및 동남부 지역에 국한되며, 모든 최우선 순위 사이트가 이곳에 위치한다.
  • 온두라스 및 과테말라: 분포된 고순위 사이트를 포함하여 보다 넓은 영역의 중-고위험 지역을 나타낸다.
• 경제적 영향 추정: 2022/2023년 생산 데이터와 시장 가격을 사용하여, 세 가지 피해 시나리오(생산 가치의 20%, 50%, 70% 감소)에 따른 연간 경제적 손실을 추정한다:
  • 저위험 피해 (20%): 약 31억 달러.
  • 중위험 피해 (50%): 약 78억 달러.
  • 고위험 피해 (70%): 약 109억 달러.
    이 수치는 직접적인 작물 손실을 나타내며, 관광, 조경 또는 문화적 가치에 미치는 간접적인 영향은 포함하지 않는다.

의의 및 주장
본 논문은 운송 경로와 기주 및 기후 적합성을 통합한 최초의 아메리카 대륙 대상 공간 명시적 CRB 위험 평가를 제시한다. 저자들은 이 분석이 취약성에 대한 결정적이고 대륙 차원의 관점을 제공하며, 감시 및 표적 개입을 안내할 특정 고위험 격자 셀을 식별한다고 주장한다.

연구는 멕시코에서의 최근 발견을 지역적 대응을 위한 촉매제로 인용하며 선제적인 생물 보안 조치의 긴급성을 강조한다. 저자들은 일단 CRB가 지역적으로 정착되면, 영향을 받는 팜 산업의 높은 경제적 가치가 다른 침입 외래종에서 사용된 것과 유사한 박멸, 봉쇄 또는 확산 억제 전략에 대한 투자를 정당화한다고 주장한다. 저자들은 개입이 없다면 이 지역이 역사적인 태평양 지역의 발생 사례와 일치하는 경제적 손실에 직면하게 될 것이며, 이는 신흥 오일 팜 산업과 지역 생계에 잠재적으로 악영ền을 미칠 수 있다고 경고한다. 본 연구의 결과는 조기 발견, 신속 대응, 그리고 국가 및 지역 이해관계자 간의 생물 보안 프로토콜 조율을 지원하기 위한 것이다.