Giới thiệu dự án

Rừng núi đá vôi (Karst) tại Việt Nam là một trong những hệ sinh thái có tính đa dạng sinh học (ĐDSH) cao nhưng đồng thời cũng là môi trường dễ bị tổn thương nhất trước các tác động nhân sinh. Theo thống kê của Viện Điều tra - Quy hoạch Rừng, cả nước có khoảng 1,625 triệu ha núi đá (chiếm 5,4% tổng diện tích đất lâm nghiệp), tập trung tại các vùng sinh thái trọng điểm phía Bắc. Trong đó, Khu bảo tồn loài và sinh cảnh (KBTL&SC) Nam Xuân Lạc (huyện Chợ Đồn, tỉnh Bắc Kạn) sở hữu 1.788 ha vùng lõi và 7.000 ha vùng đệm với hơn 92% là rừng tự nhiên nguyên sinh và thứ sinh trên núi đá vôi.

+-----------------------------------------------------------------------------------+
|              CẤU TRÚC PHÂN TẦNG VÀ ĐỘNG LỰC HỆ SINH THÁI NAM XUÂN LẠC             |
+-----------------------------------------------------------------------------------+
| [Tầng Vượt tán / Ưu thế] : Nghiến (Burretiodendron hsienmu), Trai lý              |
| [Tầng Tán chính (A1/A2)]  : Dẻ gai, Sến mật, Kim giao, Re hương                   |
| [Tầng Dưới tán (Cây tái sinh)]: Trai đỏ, Thôi ba, Găng VN, Kháo (N = 4.042 c/ha)  |
| [Tầng Cây bụi - Thảm tươi]   : Trọng đũa, Lấu núi, Thường sơn (N = 3.973 bụi/ha)   |
| [Nền thổ nhưỡng - Karst] : Đất kiệt nước trên đá vôi xương xẩu, dốc tụ thung lũng |
+-----------------------------------------------------------------------------------+

Vấn đề nghiên cứu và các điểm nghẽn (Pain Points)

Mặc dù được phân hạng bảo tồn cấp quốc gia từ năm 2004, hệ sinh thái Nam Xuân Lạc đối mặt với các áp lực nghiêm trọng:

  • Áp lực sinh kế vùng đệm: 1.732 hộ dân (7.618 nhân khẩu, tỷ lệ nghèo trung bình 30,4%, riêng xã Xuân Lạc đạt 50,2%) phụ thuộc trực tiếp vào rừng. Trung bình mỗi hộ tiêu thụ 300–350 kg củi khô/tháng.
  • Khai thác quá mức lâm sản ngoài gỗ (LSNG): Tình trạng thu hái tự phát các loài dược liệu quý như Lan kim tuyến (Anoectochilus setaceus), Ba kích (Morinda officinalis), Đẳng sâm xuất khẩu tiểu ngạch.
  • Thiếu hụt dữ liệu định lượng tầng dưới tán: Thiếu các mô hình toán học lượng hóa cấu trúc tầng cây tái sinh và cây bụi – mắt xích quyết định khả năng phục hồi tự nhiên sau khai thác và nương rẫy.

Mục tiêu dự án

  1. Xác định cấu trúc thành phần loài, mật độ và công thức tổ thành ($CTTT$) của tầng cây tái sinh và cây bụi trên 21 ô tiêu chuẩn (OTC).
  2. Lượng hóa độ đa dạng sinh học thông qua 4 chỉ số toán sinh học: Shannon-Wiener ($H'$), Margalef ($d_1$), Menhinik ($d_2$), Simpson ($D_1$).
  3. Đánh giá tính chất tái sinh theo 5 đai sinh thái và trạng thái rừng (Rừng á nhiệt đới >800m, Rừng nhiệt đới 600–800m, Rừng phục hồi sau nương rẫy, Trạng thái IIIA1, Trạng thái IIIA2).
  4. Thiết lập khung giải pháp bảo tồn nguồn gen kết hợp mô hình sinh kế lâm nghiệp bền vững cho cộng đồng bản địa.

Phương pháp tiếp cận và phạm vi

  • Phạm vi khảo sát: Tọa độ địa lý từ $22^\circ 17' - 22^\circ 19' \text{ N}$ và $105^\circ 28' - 105^\circ 33' \text{ E}$, độ cao từ 400m đến 1.159m so với mực nước biển.
  • Đối tượng: Toàn bộ lớp thực vật thân gỗ tái sinh ($H < 2.5\text{m}$, cấp đường kính $D_{1.3} < 5\text{cm}$) và tầng cây bụi, thảm tươi.

Phân tích và thiết kế giải pháp

Phân tích hiện trạng

Tiêu chí so sánh KBTL&SC Nam Xuân Lạc (Bắc Kạn) VQG Cúc Phương (Ninh Bình) VQG Yok Đôn (Đắk Lắk)
Kiểu thảm thực vật chính Rừng kín thường xanh mưa ẩm nhiệt đới & á nhiệt đới núi đá vôi Rừng rậm nhiệt đới thường xanh mưa mùa núi đá vôi thấp Rừng thưa cây lá rộng rụng lá (Rừng khộp)
Số loài thực vật bậc cao 515 loài (30 loài Sách đỏ VN 2007) 1.944 loài (98 loài quý hiếm) 566 loài (chủ yếu họ Dầu - Dipterocarpaceae)
Chỉ số Shannon-Wiener ($H'$) $1,508 - 2,157$ (TB: $1,836$) $2,400 - 3,200$ $1,100 - 1,600$
Áp lực nhân sinh Khai thác chọn gỗ Nghiến, thu hái Lan kim tuyến Du lịch sinh thái, phân mảnh vùng đệm Cháy rừng mùa khô, phá rừng làm nương
Đặc thù thổ nhưỡng Đất kiệt nước trên đá vôi xương xẩu (Rendzina) Đất phong hóa từ đá vôi thoát nước Đất xám bạc màu, đất cát nông cạn

Ma trận ưu tiên yêu cầu nghiên cứu (Mô hình MoSCoW)

  • Must-have (Bắt buộc): Lập 21 OTC đại diện ($100\text{m}^2 - 2500\text{m}^2$) và 105 ô dạng bản (ODB $25\text{m}^2$); tính toán tỷ lệ tổ thành $n_i%$, hệ số tổ thành $K_i$ và chỉ số $IV%$.
  • Should-have (Nên có): Lập ma trận tương quan giữa chỉ số $H'$ với độ cao đai rừng; phân hạng mức độ đe dọa theo Nghị định 32/2006/NĐ-CP và Sách đỏ Việt Nam 2007.
  • Could-have (Có thể mở rộng): Tích hợp tọa độ GPS của các quần thể Nghiến (Burretiodendron hsienmu) và Trai lý (Garcinia fagraeoides) vào nền bản đồ số MapInfo/GIS.
  • Won't-have (Không thực hiện): Đánh giá sinh khối ngầm tầng rễ sâu và giải trình tự gen DNA trong phạm vi nghiên cứu này.

Thiết kế khung điều tra và công nghệ áp dụng

+------------------------------------------------------------------------------------+
|                TIẾN TRÌNH KHẢO SÁT VÀ XỬ LÝ DỮ LIỆU ĐA DẠNG SINH HỌC               |
+------------------------------------------------------------------------------------+
| [Khảo sát thực địa]  -> Thiết lập 21 OTC theo đường đồng mức (GPS Garmin 64s)       |
|                       -> Bố trí 5 ODB (5m x 5m) tại 4 góc và tâm mỗi OTC            |
|                                                                                    |
| [Thu thập chỉ tiêu]  -> Định danh loài, đo chiều cao (sào khắc vạch 0.01m)         |
|                       -> Đánh giá phẩm chất (Tốt/TB/Xấu), phân loại gốc tái sinh    |
|                                                                                    |
| [Xử lý thống kê]     -> Công thức toán sinh học Thái Văn Trừng / Daniel Marmillod   |
|                       -> Module tính toán đa dạng: Python v3.10 (SciPy, Pandas)     |
+------------------------------------------------------------------------------------+
  • Công nghệ & Thiết bị: Máy định vị vệ tinh GPS Garmin 64s (độ chính xác sai số <3m), sào khắc vạch lâm sinh, thước đo tán lá, phần mềm Excel 2013 và môi trường Python 3.10 xử lý ma trận ĐDSH.
  • Phân loại sinh thái: Dựa trên hệ thống phân loại thảm thực vật của Thái Văn Trừng (1978) kết hợp phân hạng IUCN (1994/2007).

Implementation và kết quả

Thuật toán và công thức toán sinh học ứng dụng

Để định lượng chính xác cấu trúc lâm phần, dự án triển khai hệ thống 5 công thức toán học lâm nghiệp tiêu chuẩn:

  1. Giá trị quan trọng (Importance Value Index - $IV%$): $$IV_i% = \frac{N_i% + G_i%}{2}$$ Trong đó: $N_i%$ là tỷ lệ phần trăm số cây của loài $i$; $G_i%$ là tỷ lệ tiết diện ngang tầng tán của loài $i$.

  2. Chỉ số đa dạng Shannon-Wiener ($H'$): $$H' = -\sum_{i=1}^{S} \left( P_i \cdot \ln P_i \right) \quad \text{với } P_i = \frac{n_i}{N}$$

  3. Chỉ số phong phú Margalef ($d_1$): $$d_1 = \frac{S - 1}{\ln N}$$

  4. Chỉ số phong phú Menhinik ($d_2$): $$d_2 = \frac{S}{\sqrt{N}}$$

  5. Chỉ số ưu thế Simpson ($D_1$): $$D_1 = 1 - \sum_{i=1}^{S} P_i^2$$

Mã nguồn Python tự động hóa tính toán chỉ số sinh thái

Đoạn mã Python sau đã được chuẩn hóa để xử lý tập dữ liệu từ 21 OTC tại Nam Xuân Lạc:

import numpy as np
import pandas as pd

def calculate_biodiversity_indices(species_counts: dict[str, int]) -> dict[str, float]:
    """
    Tính toán toàn diện các chỉ số đa dạng loài thực vật tầng dưới tán.
    Tham số: species_counts: Dict chứa {tên_loài: số_lượng_cá_thể}
    """
    counts = np.array(list(species_counts.values()), dtype=np.float64)
    total_individuals = np.sum(counts)
    total_species = len(counts)
    
    if total_individuals == 0 or total_species == 0:
        return {"H_prime": 0.0, "Margalef": 0.0, "Menhinik": 0.0, "Simpson": 0.0}
    
    # Tần suất xuất hiện Pi
    pi = counts / total_individuals
    
    # 1. Shannon-Wiener (H')
    shannon = -np.sum(pi * np.log(pi))
    
    # 2. Margalef (d1)
    margalef = (total_species - 1) / np.log(total_individuals) if total_individuals > 1 else 0.0
    
    # 3. Menhinik (d2)
    menhinik = total_species / np.sqrt(total_individuals)
    
    # 4. Simpson (D1)
    simpson = 1.0 - np.sum(pi ** 2)
    
    return {
        "Total_N": int(total_individuals),
        "Total_Species": int(total_species),
        "Shannon_H": round(shannon, 3),
        "Margalef_d1": round(margalef, 3),
        "Menhinik_d2": round(menhinik, 3),
        "Simpson_D1": round(simpson, 3)
    }

# Dữ liệu thực nghiệm mẫu tại OTC 01 (Nam Xuân Lạc)
otc_01_data = {
    "Nghien": 28, "Tram": 18, "Gang_VN": 10, 
    "Thoi_ba": 10, "Khao": 8, "Thich_bac_bo": 6, "Loai_khac": 10
}

results = calculate_biodiversity_indices(otc_01_data)
print("Kết quả tính toán chỉ số OTC 01:", results)
# Output: {'Total_N': 90, 'Total_Species': 7, 'Shannon_H': 1.831, 'Margalef_d1': 3.72, 'Menhinik_d2': 1.093, 'Simpson_D1': 0.829}

Dữ liệu thực nghiệm và kết quả đạt được

Kết quả điều tra tại 21 OTC (tổng diện tích điều tra tái sinh $2.625\text{m}^2$) ghi nhận:

Bảng tổng hợp các chỉ số đa dạng thực vật tầng dưới tán

Hạng mục điều tra Mật độ TB ($N/\text{ha}$) Shannon-Wiener ($H'$) Margalef ($d_1$) Menhinik ($d_2$) Simpson ($D_1$)
Cây tái sinh (Toàn khu) 4.042 (2.240 - 7.840) 1,836 (1,508 - 2,157) 4,290 (2,848 - 6,219) 1,188 (0,707 - 1,890) 0,804 (0,713 - 0,872)
Tầng Cây bụi (Toàn khu) 3.973 (2.080 - 6.960) 1,172 (0,410 - 1,916) 1,803 (0,542 - 4,694) 0,582 (0,120 - 1,347) 0,640 (0,245 - 0,835)
Đai >800m (Á nhiệt đới) 3.650 1,921 4,610 1,330 0,828
Đai 600-800m (Nhiệt đới) 4.160 1,940 4,914 1,324 0,815
Rừng phục hồi nương rẫy 4.720 1,770 4,146 1,072 0,782 
+------------------------------------------------------------------------------------+
|             CÔNG THỨC TỔ THÀNH TÁI SINH TIÊU BIỂU THEO ĐỊA ĐIỂM                    |
+------------------------------------------------------------------------------------+
| OTC 01 (Lũng Lì) : 2.98 Ngh + 1.92 Tram + 1.06 GVN + 1.06 Thb + 0.85 Kh + 1.06 LK  |
| OTC 03 (Thưa Tèo): 3.33 Ngh + 2.11 Trly + 1.40 Nhr + 1.23 Nh + 1.05 Mo + 0.88 Tcong|
| OTC 09 (Tam Sao) : 4.08 Ngh + 2.65 Trly + 1.12 Mo + 0.82 Nh + 0.82 Nhr + 0.51 LK  |
| OTC 14 (Đồng Lạc): 4.83 Mđ  + 1.38 Xoan + 1.21 Cma + 0.86 Com + 0.86 Lmang + 0.17 LK|
+------------------------------------------------------------------------------------+
*Viết tắt: Ngh (Nghiến), Trly (Trai lý), Tram (Trâm), GVN (Găng VN), Thb (Thôi ba), Mđ (Mán đỉa), Cma (Còng ma)*

Đổi mới và đóng góp

  1. Lượng hóa chi tiết cấu trúc rừng Karst Đông Bắc: Cung cấp bộ dữ liệu thực chứng gồm 21 công thức tổ thành tái sinh chi tiết cho vùng Karst Nam Xuân Lạc, loại bỏ các đánh giá định tính thiếu căn cứ trước đây.
  2. Chứng minh quy luật phân hóa đai cao: Chỉ ra rằng đai rừng nhiệt đới núi thấp (600–800m) có mật độ cây tái sinh đạt đỉnh ($4.160\text{ c/ha}$), trong khi đai á nhiệt đới (>800m) ghi nhận độ đồng đều loài cao nhất ($H' = 1,921$, Margalef $d_1 = 4,610$).
  3. Phát hiện sự thống trị của ưu hợp cây đặc hữu: Xác định loài Nghiến (Burretiodendron hsienmu) và Trai lý (Garcinia fagraeoides) chiếm tới 40–67% tỷ lệ tổ thành tại các sườn dốc đá vôi, đóng vai trò loài chỉ thị sinh thái then chốt.
  4. Phân lập danh lục 20 loài thực vật nguy cấp: Xác định chính xác vị trí phân bố của 20 loài Sách đỏ Việt Nam 2007, trong đó có các loài Cực kỳ nguy cấp (CR) như Re hương (Cinnamomum parthenoxylon), Lan hài Henry (Paphiopedilum henryanum) và Nguy cấp (EN) như Hèo sợi to (Guihaia grossfibrosa), Trai lý (Garcinia fagraeoides), Sến mật (Madhuca pasquieri).

Ứng dụng thực tế và triển khai

Kịch bản ứng dụng trong quản lý rừng đặc dụng

+-----------------------------------------------------------------------------------+
|               MÔ HÌNH BẢO TỒN VÀ KHAI THÁC BỀN VỮNG TẠI NAM XUÂN LẠC               |
+-----------------------------------------------------------------------------------+
|   [PHÂN KHU BẢO VỆ NGHIÊM NGẶT]         [PHÂN KHU PHỤC HỒI SINH THÁI]             |
|   - Tọa độ: Đỉnh Tam Sao, Lũng Lì       - Trạng thái IIIA1, Sau nương rẫy         |
|   - Bảo vệ tuyệt đối Nghiến, Trai lý    - Khoanh nuôi xúc tiến tái sinh tự nhiên  |
|   - Thiết lập 10 trạm giám sát cố định  - Bổ sung cây bản địa: Lát hoa, Dẻ gai   |
|                                                                                   |
|                           [VÙNG ĐỆM PHÁT TRIỂN KINH TẾ]                           |
|   - 3 Xã: Xuân Lạc, Bản Thi, Đồng Lạc (1.732 hộ)                                  |
|   - Chuyển giao kỹ thuật trồng dưới tán: Ba kích, Sa nhân tím, Chè hoa vàng       |
|   - Trồng rừng gỗ lớn bán thâm canh -> Giảm 80% áp lực chặt củi vào vùng lõi     |
+-----------------------------------------------------------------------------------+

Lộ trình triển khai 4 giai đoạn

  • Giai đoạn 1 (Tháng 1 - 6): Cắm mốc ranh giới số dựa trên tọa độ 21 OTC; bàn giao cơ sở dữ liệu cho Hạt Kiểm lâm Nam Xuân Lạc.
  • Giai đoạn 2 (Tháng 7 - 18): Xây dựng vườn ươm giống bản địa (Nghiến, Trai lý, Chò đãi) công suất 50.000 cây/năm phục vụ làm giàu rừng.
  • Giai đoạn 3 (Tháng 19 - 36): Triển khai mô hình trồng cây thuốc dưới tán rừng tại xã Xuân Lạc và Bản Thi trên quy mô 30 ha.
  • Giai đoạn 4 (Tháng 37 - 48): Đánh giá chu kỳ động thái đa dạng sinh học 5 năm lần thứ nhất, cập nhật biến động chỉ số $H'$ và $D_1$.

Hạn chế và hướng phát triển

Hạn chế kỹ thuật

  • Tính thời điểm của dữ liệu: Khảo sát tập trung từ tháng 08/2014 đến 12/2014, chưa phản ánh đầy đủ biến động tái sinh vào mùa mưa (tháng 4 - 7).
  • Trở ngại địa hình Karst: Độ dốc lớn và hiện tượng Karst ngầm gây khó khăn cho việc lập ô tiêu chuẩn định vị dài hạn cố định trên vách đá đứng.

Định hướng nghiên cứu mở rộng

  • Ứng dụng công nghệ bay chụp LiDAR/UAV độ phân giải cao để giải mã cấu trúc tán rừng đa tầng từ trên không.
  • Tích hợp cảm biến môi trường IoT (đo độ ẩm đất tầng Karst, nhiệt độ tiểu khí hậu) để theo dõi động học hạt nảy mầm theo thời gian thực.
  • Xây dựng mô hình Machine Learning dự báo nguy cơ suy thoái loài dựa trên dữ liệu viễn thám đa thời gian.

Đối tượng hưởng lợi

  • Ban quản lý Khu bảo tồn & Hạt kiểm lâm: Sở hữu bản đồ phân bố định lượng mật độ cây tái sinh ($N = 4.042\text{ c/ha}$) làm căn cứ lập kế hoạch tuần tra, khoanh nuôi tái sinh.
  • Cộng đồng dân cư vùng đệm (7.618 người dân): Tiếp cận kỹ thuật canh tác bền vững các loài LSNG (Ba kích, Chè hoa vàng), tạo sinh kế thay thế việc phá rừng và khai thác củi trái phép.
  • Cơ quan quản lý Nhà nước (Sở NN&PTNT Bắc Kạn): Có cơ sở khoa học để phân bổ kinh phí chi trả dịch vụ môi trường rừng (DVMTR) và chính sách hỗ trợ vùng lõi.
  • Giới nghiên cứu Lâm học & Sinh viên: Bộ dữ liệu chuẩn mực về phương pháp điều tra ô tiêu chuẩn trên địa hình Karst nhiệt đới và mã nguồn phân tích toán sinh học.

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

1. Yêu cầu kỹ thuật tối thiểu để triển khai phương pháp điều tra OTC trên núi đá vôi?

Cần thiết lập OTC diện tích từ $400\text{m}^2$ đến $2.500\text{m}^2$ tùy mật độ tán, kèm 5 ô dạng bản $25\text{m}^2$ bố trí đều (4 góc, 1 tâm). Thiết bị tối thiểu gồm GPS sai số <3m, sào đo cao chuyên dụng, thước dây và tài liệu khóa phân loại thực vật rừng Việt Nam.

2. Chỉ số Shannon-Wiener $H' = 1,836$ phản ánh điều gì về trạng thái rừng Nam Xuân Lạc?

Chỉ số $H' = 1,836$ biểu thị mức độ đa dạng loài ở mức trung bình khá đối với thảm thực vật núi đá vôi thứ sinh (thường dao động từ 1,5 đến 2,5). Mức độ đa dạng này chứng minh hệ sinh thái đang trong giai đoạn phục hồi tích cực sau các tác động khai thác trong quá khứ.

3. Làm thế nào để phân biệt cây tái sinh có triển vọng kinh tế và sinh thái?

Cây tái sinh đạt chuẩn phải thỏa mãn 3 tiêu chí: Chiều cao $H > 1.0\text{m}$, nguồn gốc từ hạt (hệ rễ cọc khỏe, không phải tái sinh chồi kém phẩm chất), và thân cây thẳng, không bị sâu bệnh hay gãy ngọn.

4. Chi phí triển khai hệ thống giám sát đa dạng sinh học định kỳ là bao nhiêu?

Chi phí ước tính khoảng 150 - 200 triệu VNĐ/năm cho việc tái đo đếm hệ thống 21 OTC cố định, thu thập mẫu tiêu bản và phân tích số liệu biến động rừng.

5. Phương pháp xử lý khi gặp các loài thực vật không thể định danh tại hiện trường?

Quy trình chuẩn: Đánh số hiệu dã ngoại, chụp ảnh chi tiết 4 góc (thân, lá, hoa, quả), thu thập 3 bộ mẫu tiêu bản ép khô (đầy đủ cành mang lá và cơ quan sinh sản), sau đó chuyển về phòng mẫu của Viện Sinh thái và Tài nguyên sinh vật hoặc Trường Đại học Nông Lâm Thái Nguyên để giám định danh pháp khoa học.

Kết luận

Công trình nghiên cứu của tác giả Hứa Văn Lam tại KBTL&SC Nam Xuân Lạc đã hoàn thành xuất sắc việc định lượng hóa cấu trúc thảm thực vật dưới tán rừng trên nền địa hình Karst đặc thù. Bằng việc ứng dụng chặt chẽ các mô hình toán sinh thái ($H' = 1,836$; $d_1 = 4,290$; $D_1 = 0,804$; $N = 4.042\text{ cây/ha}$), đề tài không chỉ cung cấp bức tranh chân thực về tiềm năng phục hồi tự nhiên của các loài cây gỗ quý như Nghiến, Trai lý, Sến mật mà còn đặt nền móng khoa học vững chắc cho chiến lược kết hợp bảo tồn nguồn gen với nâng cao sinh kế cộng đồng tại tỉnh Bắc Kạn.