Tổng quan nghiên cứu

Dãy Trường Sơn là một trong những trung tâm đa dạng sinh học quan trọng bậc nhất Đông Nam Á, nơi ghi nhận sự tồn tại của nhiều loài thú đặc hữu quý hiếm. Trong số đó, Thỏ vằn Trường Sơn (Nesolagus timminsi) là loài thú cổ đại đã cách ly tiến hóa khoảng 8 triệu năm so với loài Thỏ vằn Sumatra (Nesolagus netscheri). Kể từ khi được phát hiện tại vùng biên giới Việt - Lào vào năm 1996 và mô tả khoa học chính thức vào năm 2001, loài này đã đối mặt với nguy cơ tuyệt chủng nghiêm trọng do áp lực săn bắt bằng bẫy dây phanh với quy mô lên tới hơn 10.000 bẫy bị tịch thu mỗi năm tại các khu bảo tồn trọng điểm. Hội chứng "rừng trống" cùng tình trạng suy thoái sinh cảnh rừng nguyên sinh đã đẩy quần thể loài vào tình trạng nguy cấp theo Danh lục Đỏ IUCN.

Vấn đề cốt lõi hiện nay là sự thiếu hụt trầm trọng các dữ liệu thực địa có hệ thống về tập tính sinh thái, nhịp hoạt động thời gian và mô hình dự báo phân bố không gian dưới tác động của biến đổi khí hậu. Mục tiêu cụ thể của nghiên cứu là: (1) Tổng hợp, thẩm định các điểm ghi nhận lịch sử kết hợp phỏng vấn chuyên sâu để nhận diện các vùng phân bố tiềm năng; (2) Thực hiện điều tra thực địa bằng hệ thống bẫy ảnh tự động và điều tra tuyến ban đêm để xác định đặc tính sinh thái thời gian; (3) Xây dựng mô hình phân bố loài Maxent nhằm phân tích ổ sinh thái không gian và dự báo biến động sinh cảnh đến năm 2050 theo các kịch bản phát thải khí nhà kính.

Nghiên cứu được triển khai từ tháng 1/2016 đến tháng 9/2019, bao phủ 4 khu bảo tồn trọng điểm dọc dãy Trường Sơn gồm Vườn quốc gia Bến En, Khu bảo tồn Bắc Hướng Hóa, Khu bảo tồn Sao la Huế và Khu bảo tồn Hin Nam No (Lào). Kết quả nghiên cứu cung cấp bộ dữ liệu nền tảng với độ tin cậy cao, xác lập cơ sở khoa học định lượng phục vụ công tác quy hoạch bảo tồn đa loài và quản trị rừng bền vững.

Cơ sở lý thuyết và phương pháp nghiên cứu

Khung lý thuyết áp dụng

Nghiên cứu vận dụng lý thuyết ổ sinh thái (Ecological Niche Theory) và nguyên lý Cực đại Entropy (Maximum Entropy - Maxent) do Phillips và cộng sự phát triển, cho phép mô phỏng hàm phân bố xác suất của loài dựa trên dữ liệu chỉ có điểm hiện diện (presence-only data). Khung lý thuyết này kết hợp với lý thuyết Địa lý sinh vật đảo và Cảnh quan sinh thái học để giải thích sự phân mảnh sinh cảnh và hiện tượng "rừng trống" (Empty Forest Syndrome).

Hệ thống chỉ số sinh thái định lượng cốt lõi được áp dụng bao gồm:

  1. Chỉ số phong phú tương đối (Relative Abundance Index - RAI): Tần suất ghi nhận độc lập của loài sau 100 ngày hoạt động của bẫy ảnh.
  2. Mức độ xuất hiện thô (Naïve Occupancy): Tỷ lệ phần trăm giữa số lượng trạm bẫy ảnh ghi nhận được loài trên tổng số trạm bẫy ảnh thiết lập.
  3. Độ trễ lần ghi nhận đầu (Latency to First Detection): Thời gian tích lũy tính bằng ngày từ khi kích hoạt trạm bẫy ảnh đến thời điểm ghi nhận bức ảnh đầu tiên của loài.
  4. Biến sinh khí hậu WorldClim: Tập hợp 19 biến nhiệt độ và lượng mưa chuẩn hóa toàn cầu với độ phân giải không gian 30 arc-second (khoảng 1 km²).

Phương pháp nghiên cứu

Nguồn dữ liệu của luận văn được thu thập thông qua ba trụ cột phương pháp luận đa dạng và chặt chẽ:

  • Phương pháp điều tra xã hội học và chuyên gia: Phỏng vấn bán cấu trúc 37 đối tượng gồm các nhà khoa học, cán bộ kiểm lâm và thợ săn địa phương am hiểu địa bàn, áp dụng kỹ thuật chọn mẫu quả cầu tuyết (Snowball sampling) nhằm sàng lọc thông tin tin cậy về các điểm xuất hiện chưa công bố.
  • Phương pháp điều tra thực địa theo tuyến: Thiết lập các tuyến giám sát đêm từ 19h00 đến 01h00 sáng hôm sau tại Vườn quốc gia Bến En (6 đợt, 42 ngày thực địa trong năm 2018) và Khu bảo tồn Hin Nam No (3 đợt, 30 ngày thực địa đầu năm 2019). Nhóm 3 kiểm tra viên di chuyển với vận tốc 1 km/giờ, cự ly duy trì 20 mét, sử dụng đèn đội đầu gắn kính lọc đỏ công suất cao để quét ánh lân quang mắt và xác định tọa độ qua GPS Garmin 64s.
  • Phương pháp điều tra bẫy ảnh hồng ngoại: Bố trí mạng lưới bẫy ảnh Bushnell Trophy Cam đặt ở độ cao 20 - 30 cm so với mặt đất tại Khu bảo tồn Sao la Huế (16.468 ngày bẫy hoạt động từ 01/2016 đến 01/2017) và Khu bảo tồn Bắc Hướng Hóa (1.995 ngày bẫy hoạt động từ 03/2016 đến 03/2017). Thu thập gần 300.000 bức ảnh thô, phân tích độc lập thông qua gói phần mềm camtrapR với khoảng thời gian gián đoạn 30 phút giữa các lần chụp để đảm bảo tính độc lập thống kê.
  • Phương pháp mô hình hóa không gian: Xử lý 33 điểm ghi nhận lịch sử, sử dụng gói spThin trong môi trường R để loại bỏ tự tương quan không gian ở ngưỡng cự ly 10 km, rút gọn còn 27 điểm chất lượng cao. Đánh giá đa cộng tuyến qua hệ số tương quan Pearson trong ENMTools và thuật toán Jackknife để chọn các biến môi trường tối ưu. Dự báo không gian phân bố đến năm 2050 theo 2 mô hình hoàn lưu khí hậu toàn cầu BCC-CSM1-1 và MIROC5 dưới kịch bản RCP 2.6 (phát thải thấp) và RCP 8.5 (phát thải cao).

Kết quả nghiên cứu và thảo luận

Những phát hiện chính

Thứ nhất, nghiên cứu đã xác lập chi tiết nhịp sinh học hoạt động theo thời gian của Thỏ vằn Trường Sơn. Số liệu bẫy ảnh chứng minh Nesolagus timminsi là loài hoạt động hoàn toàn về đêm. Khoảng 65% tổng số lần xuất hiện tập trung trong khung giờ từ 18h00 đến 22h00, sau đó tần suất giảm dần nhưng vẫn duy trì rải rác đến 03h00 sáng. Dữ liệu các tháng trong năm cho thấy mức độ hoạt động ghi nhận cao nhất vào mùa mưa từ tháng 9 đến tháng 12, tương ứng với chu kỳ sinh trưởng của thảm thực vật tầng thấp và cây chồi non.

Thứ hai, kết quả khảo sát thực địa đã xác nhận sự hiện diện thực tế của quần thể Thỏ vằn tại Khu bảo tồn Sao la Huế sau thời gian dài chỉ ghi nhận qua phỏng vấn gián tiếp. Ngược lại, qua 72 ngày điều tra tuyến ban đêm chuyên sâu tại Bến En và Hin Nam No, không có cá thể nào được quan sát trực tiếp ngoài tự nhiên. Điều này tái khẳng định giới hạn phân bố phía Bắc hiện tại của loài dừng lại ở khu vực Vườn quốc gia Pù Mát (Nghệ An) và Hương Sơn (Hà Tĩnh).

Thứ ba, mô hình phân bố Maxent đạt độ chính xác xuất sắc với chỉ số AUC vượt mức 0,90. Vùng phân bố tiềm năng tối ưu của loài hiện nay trải dài dọc sườn Đông và Tây dãy Trường Sơn, bị chặn tự nhiên bởi hệ thống sông Mã ở phía Bắc và hệ thống sông Sê San - sông Ba ở phía Nam. Các yếu tố môi trường đóng góp tỷ trọng cao nhất vào mô hình gồm lượng mưa mùa khô (BIO17) và dao động nhiệt độ năm (BIO7), cho thấy loài gắn liền chặt chẽ với kiểu rừng kín thường xanh mưa ẩm nhiệt đới ở độ cao từ 300 m đến 1.200 m.

Thứ tư, phân tích biến đổi khí hậu dự báo sự thu hẹp sinh cảnh đáng kể vào năm 2050. Dưới kịch bản biến đổi khí hậu cực đoan RCP 8.5 của hai mô hình BCC-CSM1-1 và MIROC5, diện tích sinh cảnh phù hợp của Thỏ vằn Trường Sơn dự kiến suy giảm khoảng 35% đến 52% so với hiện tại. Sinh cảnh thích hợp có xu hướng dịch chuyển lên các đai cao trên 800 m, làm tăng nguy cơ phân mảnh và cô lập các tiểu quần thể cục bộ.

Thảo luận kết quả

Các phát hiện về sinh thái học không gian và thời gian có thể được trực quan hóa hiệu quả thông qua biểu đồ phân bố nhịp hoạt động 24 giờ, biểu đồ cột so sánh tần suất xuất hiện theo tháng và bản đồ nhiệt GIS mô tả xác suất hiện diện sinh cảnh. Sự trùng khớp giữa tập tính kiếm ăn đêm và cấu trúc thảm thực vật ven khe suối khẳng định giả thuyết rằng vi sinh cảnh ẩm ướt cung cấp nguồn thức ăn phong phú (hơn 44 loài thực vật thuộc 16 chi) là yếu tố quyết định nơi cư trú, thay vì bản thân nguồn nước.

So sánh với các công bố quốc tế gần đây, nghiên cứu này bổ sung bằng chứng thực nghiệm quan trọng cho thấy áp lực nhân sinh (đặc biệt là bẫy bắt không chọn lọc) có tác động phân loại mạnh mẽ hơn các biến số môi trường tự nhiên đối với mật độ quần thể. Sự suy giảm diện tích dự báo đến năm 2050 đặt ra thách thức sống còn: biến đổi khí hậu sẽ đẩy các quần thể Thỏ vằn lên các "đảo sinh thái" đỉnh núi cao, nơi diện tích rừng hẹp hơn và áp lực săn trộm dễ bao vây cô lập hơn.

Đề xuất và khuyến nghị

Dựa trên kết quả định lượng về sinh thái học và mô hình phân bố, bốn nhóm giải pháp chiến lược được đề xuất:

  1. Tăng cường tuần tra tháo gỡ bẫy dây phanh và thực thi pháp luật nghiêm ngặt:

    • Hành động: Thành lập các đội tuần tra chuyên trách (Forest Guard) tuần tra định kỳ tại các tiểu khu lõi thuộc Khu bảo tồn Sao la Huế và Quảng Nam.
    • Chỉ số mục tiêu: Giảm thiểu ít nhất 80% số lượng bẫy dây phanh trong vùng lõi sinh cảnh Thỏ vằn; tịch thu 100% bẫy bắt trái phép.
    • Khung thời gian: Thực hiện liên tục giai đoạn 2020 - 2023.
    • Chủ thể thực hiện: Ban quản lý các Khu bảo tồn phối hợp với Chi cục Kiểm lâm địa phương và Tổ chức WWF Việt Nam.
  2. Quy hoạch và thiết lập hành lang đa dạng sinh học liên cảnh quan:

    • Hành động: Xây dựng hành lang kết nối sinh thái xuyên biên giới giữa Khu bảo tồn Sao la Huế, Vườn quốc gia Bạch Mã và các khu rừng giáp biên giới Lào nhằm chống phân mảnh sinh cảnh.
    • Chỉ số mục tiêu: Đảm bảo tính liên tục của hơn 100.000 ha rừng kín thường xanh đai thấp và đai trung bình.
    • Khung thời gian: Giai đoạn 2021 - 2025.
    • Chủ thể thực hiện: Bộ Nông nghiệp và Phát triển Nông thôn, Ủy ban Nhân dân các tỉnh Thừa Thiên Huế, Quảng Nam, Quảng Trị.
  3. Mở rộng cơ chế Chi trả Dịch vụ Môi trường Rừng (PFES) gắn với sinh kế bền vững:

    • Hành động: Giải ngân nguồn vốn dịch vụ môi trường rừng trực tiếp cho các cộng đồng vùng đệm; hỗ trợ chuyển đổi mô hình kinh tế thay thế việc vào rừng đặt bẫy và khai thác lâm sản ngoài gỗ.
    • Chỉ số mục tiêu: 100% số thôn bản vùng giáp ranh tham gia ký cam kết bảo vệ động vật hoang dã; nâng thu nhập bình quân từ rừng bền vững lên 25%.
    • Khung thời gian: Giai đoạn 2021 - 2026.
    • Chủ thể thực hiện: Quỹ Bảo vệ và Phát triển Rừng các tỉnh phối hợp với chính quyền cấp huyện, xã.
  4. Thiết lập hệ thống trạm giám sát bẫy ảnh dài hạn và bảo tồn ngoại vi (Ex-situ):

    • Hành động: Duy trì mạng lưới 50 điểm bẫy ảnh cố định thu thập dữ liệu định kỳ 3 tháng/lần; xây dựng kế hoạch cứu hộ nhân giống bảo tồn ngoại vi phòng ngừa nguy cơ tuyệt chủng.
    • Chỉ số mục tiêu: Cập nhật cơ sở dữ liệu di truyền học và biến động quần thể hàng năm với độ tin cậy thống kê 95%.
    • Khung thời gian: Giai đoạn 2022 - 2030.
    • Chủ thể thực hiện: Viện Sinh thái và Tài nguyên sinh vật, Trường Đại học Khoa học Tự nhiên - ĐHQGHN, Viện Nghiên cứu Vườn thú và Động vật Hoang dã Leibniz (IZW).

Đối tượng nên tham khảo luận văn

  1. Cơ quan quản lý lâm nghiệp và Ban quản lý các Vườn quốc gia: Ứng dụng bản đồ dự báo phân bố tiềm năng Maxent để tối ưu hóa việc phân vùng bảo vệ nghiêm ngặt và lập kế hoạch tuần tra phá bẫy có trọng tâm.
  2. Các tổ chức phi chính phủ và quỹ bảo tồn quốc tế (WWF, FFI, IUCN, KfW): Sử dụng các chỉ số RAI và dữ liệu nhịp hoạt động làm thước đo đánh giá hiệu quả (baseline metrics) cho các dự án đầu tư bảo tồn đa dạng sinh học dãy Trường Sơn.
  3. Giảng viên, nghiên cứu sinh và học viên cao học ngành Sinh thái học, Khoa học Môi trường, Quản lý Tài nguyên: Khai thác quy trình kỹ thuật tích hợp giữa bẫy ảnh hiện đại, gói phần mềm camtrapR, xử lý không gian qua spThin và mô hình hóa sinh khí hậu WorldClim.
  4. Chuyên gia phân tích chính sách và nhà hoạch định quy hoạch không gian bảo tồn: Sử dụng kịch bản dự báo suy giảm sinh cảnh năm 2050 để tích hợp yếu tố thích ứng biến đổi khí hậu vào chiến lược bảo tồn đa dạng sinh học quốc gia.

Câu hỏi thường gặp

Loài Thỏ vằn Trường Sơn khác biệt như thế nào so với Thỏ vằn Sumatra?
Thỏ vằn Trường Sơn (Nesolagus timminsi) và Thỏ vằn Sumatra (Nesolagus netscheri) phân tách di truyền khoảng 8 triệu năm trước. Về hình thái, loài Trường Sơn có các sọc vằn đen nâu chạy dọc lưng tương đồng nhưng khác biệt rõ rệt về cấu trúc xương sọ, hình thái cung gò má và cấu tạo bộ răng thích nghi với các loại thức ăn thực vật vùng Đông Dương.

Tại sao phương pháp bẫy ảnh hồng ngoại lại vượt trội hơn điều tra quan sát trực tiếp đối với loài này?
Nesolagus timminsi là loài thú nhỏ, sống đơn độc, nhút nhát và chỉ hoạt động ban đêm từ 18h00 đến 03h00 sáng. Bẫy ảnh Bushnell Trophy Cam đặt cao 20 - 30 cm với đèn flash hồng ngoại vô hình giúp ghi nhận liên tục 24/7 mà không làm xáo trộn tập tính tự nhiên, thu về hơn 16.000 ngày bẫy dữ liệu khách quan.

Biến đổi khí hậu ảnh hưởng thế nào đến vùng phân bố của Thỏ vằn Trường Sơn vào năm 2050?
Dưới tác động của kịch bản phát thải cao RCP 8.5, nền nhiệt độ tăng và dao động lượng mưa mùa khô khiến diện tích sinh cảnh phù hợp sụt giảm từ 35% đến 52%. Vùng phân bố bị co cụm và đẩy dần lên các dãy núi cao trên 800 m, làm tăng nguy cơ suy thoái nguồn gen do cách ly địa hình.

Mối đe dọa trực tiếp nghiêm trọng nhất đối với sự tồn vong của Thỏ vằn hiện nay là gì?
Mối đe dọa nguy hiểm nhất là tình trạng bẫy bắt bằng dây phanh diễn ra tràn lan với giá thành rẻ và tính sát thương cao. Riêng tại Khu bảo tồn Sao la Huế, kiểm lâm thu giữ xấp xỉ 10.000 bẫy/năm, gây ra hội chứng "rừng trống" đe dọa trực tiếp các loài thú kiếm ăn mặt đất.

Thuật toán Maxent xử lý bài toán thiếu dữ liệu điểm ghi nhận của loài quý hiếm như thế nào?
Maxent tối ưu hóa hàm phân bố xác suất trên nền tảng chỉ cần điểm hiện diện (presence-only), không đòi hỏi điểm vắng mặt. Khi kết hợp cùng công cụ lọc không gian spThin ở bán kính 10 km, thuật toán loại bỏ hiện tượng sai lệch lấy mẫu (sampling bias), tạo ra mô hình dự báo sinh thái với chỉ số tin cậy AUC > 0,90.

Kết luận

  • Nghiên cứu đã hệ thống hóa và chuẩn hóa 27 điểm dữ liệu phân bố không gian có độ tin cậy cao của loài Thỏ vằn Trường Sơn dọc dải Trường Sơn.
  • Thực hiện chương trình giám sát thực địa quy mô lớn với 18.463 ngày bẫy ảnh kết hợp 72 ngày điều tra đêm theo tuyến, ghi nhận thành công sự hiện diện của loài tại Khu bảo tồn Sao la Huế.
  • Xác định quy luật nhịp hoạt động thời gian đặc thù của loài với 65% tần suất xuất hiện tập trung trong khung giờ từ 18h00 đến 22h00.
  • Xây dựng mô hình Maxent dự báo chính xác vùng phân bố hiện tại và cảnh báo nguy cơ mất từ 35% đến 52% sinh cảnh thích hợp vào năm 2050 do biến đổi khí hậu.
  • Đề xuất 4 nhóm giải pháp bảo tồn đồng bộ kết hợp giữa tháo gỡ bẫy bắt cơ học, bảo vệ hành lang sinh cảnh, sinh kế cộng đồng và giám sát ngoại vi dài hạn.

Các cơ quan quản lý và viện nghiên cứu cần khẩn trương đưa các vùng sinh cảnh cốt lõi của Thỏ vằn Trường Sơn vào danh mục ưu tiên bảo vệ đặc biệt trong kế hoạch hành động đa dạng sinh học giai đoạn 2021 - 2030, kết hợp ứng dụng công nghệ giám sát số để bảo tồn loài thú đặc hữu vô giá này trước nguy cơ tuyệt chủng.