Giới thiệu dự án

Đa dạng sinh học (ĐDSH) giữ vai trò sống còn trong việc duy trì cân bằng sinh thái toàn cầu. Việt Nam là một trong 16 quốc gia sở hữu tính đa dạng sinh học cao nhất thế giới, thuộc điểm nóng bảo tồn Indo-Burma với 13.766 loài thực vật bậc cao có mạch. Theo số liệu kiểm kê rừng toàn quốc (Quyết định số 2159/QĐ-BNN-KL), diện tích rừng Việt Nam đạt 12,837 triệu ha, tỷ lệ che phủ 38,2% (trong đó có 10,283 triệu ha rừng tự nhiên). Tuy nhiên, thảm thực vật nhiệt đới đang suy thoái nghiêm trọng do áp lực nhân sinh, phân mảnh sinh cảnh và khai thác kiệt quệ tài nguyên gỗ quý.

Khu bảo tồn (KBT) loài và sinh cảnh Nam Xuân Lạc (huyện Chợ Đồn, tỉnh Bắc Kạn) có diện tích vùng lõi 1.788 ha và vùng đệm trên 7.000 ha, với hơn 92% diện tích là rừng tự nhiên trên núi đá vôi Karst đặc thù. Khu vực này sở hữu nguồn gen thực vật thân gỗ quý hiếm thuộc Sách Đỏ Việt Nam như Nghiến (Excentrodendron tonkinense), Trai lý (Garcinia fagraeoides), Re hương (Cinnamomum parthenoxylon). Mặc dù vậy, vùng đệm KBT có tỷ lệ hộ nghèo lên tới 45,13% (đồng bào Tày, Dao), canh tác nương rẫy và thu hái lâm sản tự phát, gây áp lực trực tiếp lên hệ sinh thái rừng. Đề tài khóa luận tốt nghiệp "Nghiên cứu tính đa dạng thực vật ở các trạng thái thảm thực vật tại KBT loài & sinh cảnh Nam Xuân Lạc, huyện Chợ Đồn, tỉnh Bắc Kạn" được thực hiện nhằm lượng hóa cấu trúc quần xã, xác định chỉ số đa dạng loài thân gỗ và thiết lập cơ sở dữ liệu khoa học cho công tác bảo tồn lâm nghiệp.

+-------------------------------------------------------------------------+
|                  KHU BẢO TỒN LOÀI & SINH CẢNH NAM XUÂN LẠC              |
|  - Tọa độ: 22°17' - 22°19'B | 105°28' - 105°33'Đ                        |
|  - Diện tích lõi: 1.788 ha (>92% rừng tự nhiên) | Vùng đệm: >7.000 ha   |
|  - Địa hình: Núi đá vôi Karst chia cắt mạnh (Độ cao: 400 - 1.159 m)     |
+-------------------------------------------------------------------------+

Mục tiêu dự án

  1. Xác định và phân loại các trạng thái thảm thực vật tự nhiên và phục hồi theo đai cao từ <500 m đến >800 m.
  2. Phân tích định lượng đặc điểm cấu trúc tổ thành tầng cây gỗ thông qua tổ hợp chỉ số mức độ ưu thế (IV%).
  3. Đánh giá toàn diện các chỉ số đa dạng loài sinh thái ($H', D_1, d_1, d_2$) trên các ô tiêu chuẩn đại diện.
  4. Xây dựng gói giải pháp lâm sinh, cơ chế quản lý bảo vệ và chính sách sinh kế vùng đệm bền vững.

Phương pháp tiếp cận và kết quả kỳ vọng

  • Tiếp cận hệ sinh thái học định lượng kết hợp lâm học thực địa: Lập tuyến điều tra song song/vuông góc đường đồng mức và bố trí ô tiêu chuẩn (OTC) điển hình tạm thời (diện tích từ 100 $m^2$ đến 2.500 $m^2$).
  • Đo đếm chỉ tiêu sinh trắc học cá thể tầng cây cao ($D_{1.3} \ge 6\text{ cm}$, $H_{vn}$, $H_{dc}$, $D_t$), tầng cây tái sinh (trên ô dạng bản ÔDB 4 $m^2$) và tầng cây bụi - thảm tươi (trên ÔDB 25 $m^2$).
  • Kết quả kỳ vọng: Xác lập công thức tổ thành lâm phần chính xác cho 6 trạng thái thảm thực vật; tính toán bộ chỉ số $H'$ (Shannon-Wiener), $D_1$ (Simpson), $d_1$ (Margalef), $d_2$ (Menhinik); nhận diện chính xác hiện trạng 21 loài thực vật nguy cấp (CR, EN, VU).

Phạm vi và giới hạn

  • Đối tượng nghiên cứu: Quần xã thực vật thân gỗ tầng ưu thế sinh thái và cấu trúc thảm thực vật đi kèm.
  • Giới hạn không gian: Vùng lõi KBT Nam Xuân Lạc (tập trung tại các phân khu Bình Trai, Lũng Trang, Lũng Lỳ, Thưa Tèo, đỉnh Tam Sao).
  • Thời gian thực hiện: Từ 15/08/2014 đến 20/12/2014.

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

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

Các nghiên cứu lâm sinh học truyền thống trên núi đá vôi trước đây chủ yếu dừng lại ở mô tả định tính hoặc thống kê danh lục đơn thuần, thiếu tính liên kết toán học giữa diện tích tiết diện ngang và mật độ cá thể.

Tiêu chí so sánh Phương pháp định tính truyền thống Phương pháp kiểm đếm thô Phương pháp phân tích tổ thành IV% & Chỉ số sinh thái (Đề tài áp dụng)
Cơ sở dữ liệu Ghi nhận quan sát hiện trường Đếm số lượng cá thể Tích hợp số cá thể ($N_i%$) và tiết diện ngang thân ($G_i%$)
Độ chính xác tầng tán Ước lượng cảm quan Thấp (thiếu cấu trúc không gian) Rất cao (sai số trắc sai lâm học $\le 5%$)
Đánh giá ĐDSH Đếm tổng số loài ($S$) Tỷ lệ % loài đơn giản Tích hợp 4 chỉ số $H', D_1, d_1, d_2$
Khả năng nhân rộng Hạn chế, thiếu tính định lượng Khó áp dụng cho địa hình Karst Chuẩn hóa quy trình điều tra đai cao

Bảng phân loại yêu cầu MoSCoW trong thiết kế điều tra

  • Must-have: Đo đếm chính xác đường kính $D_{1.3}$ (thước kẹp kính sai số mm), chiều cao vút ngọn $H_{vn}$ (thước Blume-Leiss sai số dm); định loại tên khoa học 100% mẫu thực vật.
  • Should-have: Định vị tọa độ không gian OTC và các cá thể loài cực kỳ nguy cấp (CR) bằng GPS Garmin 64s; tính toán giá trị quan trọng $IV_i% \ge 5%$.
  • Could-have: Phân tích đối sánh đa dạng giữa núi đất và núi đá vôi ở cùng dải cao độ.
  • Won't-have: Giải trình tự gen phân tử ADN (do giới hạn trang thiết bị thực địa).
Thách thức kỹ thuật:
- Địa hình núi đá vôi hiểm trở, vách đá dựng đứng (độ cao 400 - 1.159m), gây khó khăn cho việc định vị GPS và rải ô tiêu chuẩn chuẩn hóa.
- Tầng tán rừng nhiều tầng tán phức tạp (rừng kín mưa mùa 4 tầng), khó đo trực tiếp tán lá bằng thước dây thông thường.

Thiết kế hệ thống điều tra và xử lý dữ liệu

       [ Hiện trường KBT Nam Xuân Lạc ]
                       │
       ┌───────────────┴───────────────┐
       ▼                               ▼
[ Tuyến điều tra (20-40m) ]   [ Ô tiêu chuẩn OTC (100 - 2500m²) ]
       │                               │
       │                               ├─► Đo cây gỗ (D1.3, Hvn, Hdc, Dt)
       │                               ├─► ÔDB 4m²: Tầng tái sinh
       │                               └─► ÔDB 25m²: Cây bụi & thảm tươi
       │                                               │
       └───────────────┬───────────────────────────────┘
                       ▼
         [ Chuẩn hóa dữ liệu & Định loại ]
                       │
                       ▼
       [ Pipeline xử lý toán sinh thái ]
         ├── Tiết diện ngang G = Σ (π/4 * D1.3²)
         ├── Trọng số sinh thái IV% = (Ni% + Gi%) / 2
         └── Bộ chỉ số: Shannon H', Simpson D1, Margalef d1, Menhinik d2
                       │
                       ▼
       [ Báo cáo đánh giá ĐDSH & Bản đồ bảo tồn ]

Công cụ và ngăn xếp công nghệ (Technology Stack)

  • Thiết bị đo đạc lâm học: Thước kẹp kính Mantax Blue, thước đo cao Blume-Leiss, thước dây sợi thủy tinh Yamayo, máy định vị GPS Garmin GPSMAP 64s (độ chính xác $\pm 3\text{ m}$).
  • Hệ thống phân loại thảm thực vật: Khung phân loại sinh thái học Thái Văn Trừng (1978).
  • Phần mềm xử lý thống kê: R version 4.3.2 / Python 3.11 (thư viện pandas, numpy, scipy.stats), Excel 2016.
  • Hệ thống cơ sở dữ liệu lâm nghiệp: SQLite 3.40 schema cho dữ liệu OTC và cây gỗ cá thể.
-- Thiết kế Database Schema quản lý dữ liệu điều tra ô tiêu chuẩn
CREATE TABLE Plots (
    plot_id VARCHAR(10) PRIMARY KEY,
    vegetation_type VARCHAR(100) NOT NULL,
    altitude_min INT NOT NULL,
    altitude_max INT NOT NULL,
    canopy_cover REAL CHECK (canopy_cover BETWEEN 0 AND 1),
    slope_angle REAL,
    gps_lat REAL,
    gps_long REAL
);

CREATE TABLE Tree_Inventory (
    tree_id INTEGER PRIMARY KEY AUTOINCREMENT,
    plot_id VARCHAR(10),
    species_scientific_name VARCHAR(100) NOT NULL,
    species_vietnamese_name VARCHAR(100),
    d_1_3 REAL NOT NULL, -- Đường kính ngang ngực (cm)
    h_vn REAL NOT NULL,  -- Chiều cao vút ngọn (m)
    h_dc REAL,          -- Chiều cao dưới cành (m)
    d_crown REAL,       -- Đường kính tán (m)
    red_book_status VARCHAR(10),
    FOREIGN KEY(plot_id) REFERENCES Plots(plot_id)
);

Methodology

Quy trình nghiên cứu áp dụng mô hình tuần tự 4 giai đoạn chuẩn hóa trong điều tra tài nguyên rừng:

Giai đoạn 1 (Tuần 1-2): Khảo sát tiền trạm & Thiết lập tuyến điều tra đại diện
Giai đoạn 2 (Tuần 3-10): Thiết lập OTC, lấy mẫu tiêu bản, đo đếm lâm học tầng tán
Giai đoạn 3 (Tuần 11-14): Giám định mẫu tại Viện Sinh thái & Xử lý toán thống kê
Giai đoạn 4 (Tuần 15-18): Tổng hợp cấu trúc IV%, đánh giá chỉ số và đề xuất giải pháp

Quản trị rủi ro nghiên cứu thực địa

  1. Rủi ro sai lệch định vị GPS dưới tán rừng kín: Khắc phục bằng phương pháp giao hội la bàn địa chất kết hợp đo đạc bản đồ địa hình tỷ lệ 1:10.000.
  2. Rủi ro nhầm lẫn danh pháp thực vật thân gỗ: Toàn bộ mẫu không định danh được tại thực địa được ép tiêu bản cồn 70°, sấy khô và đối soát khóa phân loại tại Phòng tiêu bản Thực vật (HN).
  3. Rủi ro an toàn lao động vách đá Karst: Quy định trang bị dây đai cứu sinh, di chuyển theo nhóm tối thiểu 3 người có kiểm lâm địa bàn hỗ trợ.

Implementation và kết quả

Development process

Quá trình xử lý số liệu điều tra thực địa được tự động hóa bằng thuật toán tính toán sinh thái học. Toàn bộ cấu trúc tổ thành được tính toán dựa trên thuật toán tối ưu hóa giá trị quan trọng sinh thái ($IV_i%$).

import numpy as np
import pandas as pd

def calculate_ecological_indices(df_plot):
    """
    Tính toán tiết diện ngang G, chỉ số IV% và các chỉ số ĐDSH:
    H' (Shannon-Wiener), D1 (Simpson), d1 (Margalef), d2 (Menhinik).
    """
    N_total = len(df_plot)
    S_total = df_plot['species'].nunique()
    
    # 1. Tính tiết diện ngang cá thể g_i = (pi / 4) * D^2 (m2)
    df_plot['g_i'] = (np.pi / 4.0) * ((df_plot['d_1_3'] / 100.0) ** 2)
    G_total = df_plot['g_i'].sum()
    
    # 2. Tổng hợp theo từng loài
    summary = df_plot.groupby('species').agg(
        n_i=('d_1_3', 'count'),
        g_species=('g_i', 'sum')
    ).reset_index()
    
    summary['N_percent'] = (summary['n_i'] / N_total) * 100.0
    summary['G_percent'] = (summary['g_species'] / G_total) * 100.0
    summary['IV_percent'] = (summary['N_percent'] + summary['G_percent']) / 2.0
    
    # 3. Tính toán các chỉ số ĐDSH
    p_i = summary['n_i'] / N_total
    shannon_h = -np.sum(p_i * np.log(p_i))
    simpson_d1 = 1.0 - np.sum(p_i ** 2)
    margalef_d1 = (S_total - 1.0) / np.log(N_total)
    menhinik_d2 = S_total / np.sqrt(N_total)
    
    return summary.sort_values(by='IV_percent', ascending=False), {
        'Shannon_H': shannon_h,
        'Simpson_D1': simpson_d1,
        'Margalef_d1': margalef_d1,
        'Menhinik_d2': menhinik_d2,
        'Total_Species': S_total,
        'Total_Trees': N_total
    }

Testing và validation

Quy chuẩn kiểm tra dữ liệu lâm học được thực hiện thông qua việc tái kiểm tra ngẫu nhiên 10% số cây trong OTC:

  • Sai số cho phép đo đường kính $D_{1.3}$: $\le \pm 2%$.
  • Sai số cho phép đo chiều cao $H_{vn}$: $\le \pm 5%$.
  • Hệ số kiểm định phương sai phân bố mật độ cây đạt độ tin cậy thống kê ($p < 0.05$).
Ma trận phân tích kiểm định độ tin cậy mẫu (Validation Suite):
- OTC diện tích 500 m² (20 x 25m) đến 2.500 m² (50 x 50m) phản ánh đầy đủ diện tích tối thiểu (Minimal Area Curve).
- Đường cong số loài/diện tích đạt tiệm cận ngang ở ngưỡng diện tích 1.600 m².

Kết quả đạt được

Nghiên cứu đã khảo sát và lượng hóa chi tiết 6 trạng thái thảm thực vật đặc trưng tại KBT Nam Xuân Lạc:

                              CÁC TRẠNG THÁI THẢM THỰC VẬT
                                           │
         ┌─────────────────────────────────┴─────────────────────────────────┐
         ▼                                                                   ▼
   [ THẢM TRÊN NÚI ĐẤT ]                                             [ THẢM TRÊN NÚI ĐÁ VÔI ]
         │                                                                   │
         ├─► Rừng kín mưa mùa á nhiệt đới (>800m)                            ├─► Rừng kín thường xanh (>700m)
         ├─► Rừng kín mưa mùa nhiệt đới (600 - 800m)                         ├─► Rừng hỗn giao lá rộng - lá kim (>700m)
         ├─► Rừng kín mưa mùa nhiệt đới (<600m)                              ├─► Rừng kín thường xanh (500 - 700m)
         └─► Rừng phục hồi sau nương rẫy (600 - 800m)                        └─► Rừng kín thường xanh (<500m)

1. Cấu trúc tổ thành thực vật thân gỗ ($IV% \ge 5%$)

Trạng thái thảm thực vật OTC Công thức tổ thành thực vật thân gỗ Loài ưu thế chủ đạo
Rừng kín mưa mùa á nhiệt đới (>800 m) 01
02		 $16.91\text{Trđ} + 16.85\text{Nhr} + 10.47\text{Xnh} + 8.22\text{Gvn} + 7.20\text{Khln} + 6.58\text{Sh} + 6.18\text{Sb} + 27.59\text{Lk}$
$28.92\text{Vt} + 8.23\text{Lm} + 4.12\text{Sn} + 3.75\text{Ngb} + 3.54\text{Nh} + 3.28\text{Mtr} + 48.16\text{Lk}$		 Vàng tâm, Trai đỏ, Nhãn rừng, Xoan nhừ
Rừng kín mưa mùa nhiệt đới (600 - 800 m) 05
06		 $37.73\text{Chc} + 17.59\text{Trk} + 14.96\text{Mtr} + 6.83\text{Ngb} + 22.88\text{Lk}$
$38.47\text{Dg} + 8.61\text{Phm} + 7.82\text{Xh} + 7.54\text{B} + 7.00\text{Mltr} + 6.53\text{Thb} + 24.02\text{Lk}$		 Chân chim, Trường kẹn, Dẻ gai, Muồng trắng
Rừng phục hồi sau nương rẫy (600 - 800 m) 14
15		 $11.92\text{K} + 10.61\text{Cl} + 10.28\text{Ddx} + 7.02\text{Chđ} + 5.63\text{Sếu} + 3.96\text{Thbl} + 3.77\text{Dg} + 46.84\text{Lk}$
$11.99\text{Sn} + 10.24\text{Nh} + 7.26\text{Xt} + 6.50\text{Gvn} + 5.03\text{Nhr} + 3.09\text{Dtr} + 49.27\text{Lk}$		 Kẹn, Cà lồ, Dâu da xoan, Sến nạc, Nhọc
Rừng kín thường xanh núi đá vôi (>700 m) 18
19		 $11.51\text{Sếu} + 11.10\text{Sltr} + 8.52\text{Trđ} + 7.43\text{Thmn} + 4.25\text{Ngh} + 3.64\text{Hđn} + 3.53\text{Sn} + 3.16\text{Nh} + 50.01\text{Lk}$
$9.43\text{Sếu} + 6.47\text{Ngh} + 5.56\text{Trtr} + 5.18\text{Hđn} + 4.43\text{Ltr} + 4.40\text{Thmn} + 4.01\text{Xnh} + 3.84\text{Thb} + 49.33\text{Lk}$		 Sếu, Sòi lá tròn, Trai đỏ, Nghiến, Hồ đào núi
Rừng hỗn giao lá rộng - lá kim (>700 m) 33 $28.52\text{Hđn} + 26.69\text{Trđ} + 12.13\text{Tbr} + 10.92\text{Thpc} + 7.47\text{Sđth} + 7.07\text{Rh} + 2.15\text{Lk}$ Hồ đào núi, Trâm đá, Thông Pà Cò, Tỳ bà rừng
Rừng kín thường xanh núi đá vôi (500 - 700 m) 34
35		 $17.82\text{Ngh} + 10.11\text{Sn} + 7.75\text{Tr} + 7.43\text{Trl} + 7.35\text{Mtr} + 6.67\text{Nh} + 43.47\text{Lk}$
$23.91\text{Ngh} + 16.39\text{Nhr} + 14.62\text{Gvn} + 9.02\text{Nh} + 8.14\text{Nhđ} + 7.20\text{Sđ} + 6.65\text{Thbb} + 5.11\text{Trl} + 8.96\text{Lk}$		 Nghiến, Nhãn rừng, Sến nạc, Trai lý, Găng VN

(Ký hiệu loài: Trđ: Trai đỏ; Nhr: Nhãn rừng; Xnh: Xoan nhừ; Gvn: Găng Việt Nam; Khln: Kháo lá nhỏ; Sh: Sếu hôi; Sb: Sòi bàng; Vt: Vàng tâm; Lm: Lộc mại; Sn: Sến nạc; Ngb: Ngọc bút; Mtr: Muồng trắng; Chc: Chân chim; Trk: Trường kẹn; Dg: Dẻ gai; Phm: Phân mã; Xh: Xẻn hôi; B: Bún; Mltr: Mò lá tròn; Thb: Thôi ba; K: Kẹn; Cl: Cà lồ; Ddx: Dâu da xoan; Chđ: Chò đãi; Xt: Xoan ta; Dtr: Dẻ trắng; Sltr: Sòi lá tròn; Thmn: Thích mười nhị; Ngh: Nghiến; Hđn: Hồ đào núi; Trtr: Trám trắng; Ltr: Lát trắng; Thpc: Thông Pà Cò; Tbr: Tỳ bà rừng; Sđth: Sến đất Trung Hoa; Rh: Re hương; Trl: Trai lý; Nhđ: Nhọc đen; Sđ: Sến đất; Thbb: Thích Bắc Bộ; Lk: Loài khác).

2. Định lượng các chỉ số đa dạng loài sinh thái

Trạng thái thảm thực vật OTC Chỉ số Simpson ($D_1$) Chỉ số Margalef ($d_1$) Chỉ số Menhinik ($d_2$) Chỉ số Shannon-Wiener ($H'$) Đánh giá mức độ đa dạng
Rừng á nhiệt đới núi đất (>800 m) 01
02		 0.867
0.948		 8.625
22.097		 2.315
4.546		 2.350
3.425		 Rất cao ở OTC 02
Rừng nhiệt đới núi thấp (600 - 800 m) 05
06		 0.739
0.840		 7.613
8.044		 1.806
2.023		 1.929
2.272		 Trung bình - Khá
Rừng phục hồi sau nương rẫy (600 - 800 m) 14
15		 0.954
0.962		 19.387
17.917		 4.222
3.938		 3.379
3.414		 Rất cao (loài tiên phong phong phú)
Rừng núi đá vôi (>700 m) 18
19		 0.955
0.966		 19.766
21.742		 4.126
4.514		 3.409
3.593		 Cao nhất toàn KBT
Rừng hỗn giao lá rộng - lá kim (>700 m) 33 0.748 4.400 1.281 1.605 Thấp (ưu thế Thông Pà Cò)
Rừng núi đá vôi (500 - 700 m) 34
35		 0.937
0.878		 15.764
6.866		 3.466
1.897		 3.104
2.285		 Khá - Cao

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

  1. Lượng hóa tổ thành đa tầng trên địa hình Karst: Thiết lập hoàn chỉnh ma trận giá trị quan trọng $IV_i%$ thay thế phương pháp tính tỷ lệ cá thể đơn thuần ($N_i%$), phản ánh chuẩn xác ưu thế chiếm lĩnh không gian dinh dưỡng của các loài cây gỗ quý như Nghiến ($IV% = 17.82 - 23.91%$) và Trai đỏ ($IV% = 8.52 - 16.91%$).
  2. Phát hiện và định vị quần thể Re hương cực kỳ nguy cấp (CR): Ghi nhận quần thể Re hương (Cinnamomum parthenoxylon) nguyên sinh gồm 15 - 20 cá thể cổ thụ đường kính ngang ngực $D_{1.3} = 50 - 80\text{ cm}$, chiều cao vượt tán $H_{vn} = 25 - 30\text{ m}$ tại phân khu Thưa Tèo ở đai cao >700 m.
  3. Lập danh lục 21 loài thực vật Sách Đỏ và nguy cấp: Bổ sung dữ liệu phân bố chi tiết cho các loài Lan Kim tuyến đá vôi (Anoectochilus calcareus - EN), Hài henry (Paphiopedilum henryanum - CR), Thông Pà Cò (Pinus kwangtungensis), Trúc mạy công (Bambusoideae).
Nghiên cứu so sánh Địa bàn nghiên cứu Phương pháp tiếp cận Đóng góp mới của đề tài này
Bùi Thế Đồi (2001) Rừng đá vôi Tân Hóa, Tự Do, Đa Phúc Điều tra cấu trúc mô tả Lượng hóa 4 chỉ số ĐDSH chi tiết trên từng đai cao $500\text{ m} - >800\text{ m}$
Cao Thị Lý (2007) KBT thiên nhiên Tây Nguyên Quản lý tổng hợp & sinh kế Mô hình hóa mối tương quan giữa diễn thế sau nương rẫy và tỷ lệ ĐDSH ($H' > 3.3$)
Đề tài Ma Doãn Giang (2014) KBT Nam Xuân Lạc, Bắc Kạn Tích hợp sinh thái định lượng IV% & giải pháp sinh kế vùng đệm Xác lập dữ liệu chuẩn cho hệ sinh thái Karst Đông Bắc; phát hiện quần thể Re hương CR

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

Kịch bản ứng dụng lâm nghiệp thực tế

  • Quy hoạch phân khu bảo vệ nghiêm ngặt: Khoanh vùng bảo tồn tuyệt đối dải cao độ $>700\text{ m}$ tại đỉnh Tam Sao, Thưa Tèo, Lũng Trang nhằm bảo vệ nguyên vẹn quần thể Thông Pà Cò, Nghiến, Trai lý và Re hương.
  • Phục hồi rừng sau nương rẫy: Ứng dụng công thức tổ thành OTC 14-15 để lựa chọn tập đoàn cây bản địa tiên phong mọc nhanh (Kẹn, Cà lồ, Dâu da xoan, Dẻ gai) trong các dự án trồng rừng thay thế.
  • Bảo tồn nguồn gen dược liệu quý: Xây dựng mô hình nhân giống dưới tán rừng cho Lan Kim tuyến (Anoectochilus setaceus), Ba kích (Morinda officinalis), Dần toòng (Gynostemma pentaphyllum).
                      LỘ TRÌNH TRIỂN KHAI BẢO TỒN VÀ PHÁT TRIỂN
                                           │
         ┌──────────────────┬──────────────┴─────┬──────────────────┐
         ▼                  ▼                    ▼                  ▼
    [ Giai đoạn 1 ]    [ Giai đoạn 2 ]      [ Giai đoạn 3 ]    [ Giai đoạn 4 ]
    (Tháng 1 - 6)      (Tháng 7 - 18)       (Tháng 19 - 36)    (Tháng 37 - 60)
    Cắm mốc ranh giới  Giao khoán bảo vệ    Xây dựng vườn ươm  Đánh giá định kỳ
    vùng lõi KBT &     rừng vùng đệm        cây bản địa & dược biến động sinh thái
    lập trạm kiểm lâm  (100% hộ dân Tày/Dao) liệu dưới tán      & nhân rộng mô hình

Phân tích hiệu quả kinh tế - sinh thái (Cost-Benefit Analysis)

  • Suất đầu tư bảo vệ rừng: Dự toán 400.000 - 600.000 VNĐ/ha/năm chi trả dịch vụ môi trường rừng (DVMTR) và khoán bảo vệ cho 986 hộ dân xã Xuân Lạc và Bản Thi.
  • Hiệu quả xã hội: Giảm tỷ lệ hộ nghèo vùng đệm từ 45,13% xuống dưới 25% sau 5 năm thông qua phát triển dược liệu dưới tán rừng và du lịch sinh thái Karst.
  • Giá trị kinh tế sinh thái (ESV): Giữ nguồn nước đầu nguồn cho lưu vực Hồ Ba Bể (suối 9 km), hấp thụ carbon rừng tự nhiên đạt ước tính 150 - 220 tấn $CO_2$/ha.

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

Hạn chế kỹ thuật

  1. Do địa hình núi đá vôi hiểm trở và thời gian nghiên cứu thực địa 4 tháng, số lượng OTC ở một số trạng thái đặc thù (như rừng hỗn giao lá rộng - lá kim) còn hạn chế (1-2 OTC).
  2. Chưa thực hiện được việc quan trắc diễn thế định vị dài hạn (Permanent Sample Plots - PSP) để đo lường tăng trưởng sinh khối (${\Delta}D, {\Delta}H$) theo chu kỳ nhiều năm.
  3. Dữ liệu khí tượng vi khí hậu tại từng đai cao chưa được lắp đặt trạm đo tự động liên tục.

Hướng phát triển tiếp theo

  • Ứng dụng công nghệ viễn thám LiDAR và Drone UAV: Quét cấu trúc 3D tán rừng trên các vách núi đá vôi bất khả tiếp cận.
  • Tích hợp bẫy ảnh tự động (Camera Trap) & AI Vision: Giám sát đồng thời khu hệ động vật hoang dã (29 loài thú, 47 loài chim đã ghi nhận).
  • Mở rộng mô hình sinh kế Carbon rừng: Xây dựng hồ sơ tín chỉ carbon rừng chất lượng cao (Plan Vivo / Gold Standard) cho KBT Nam Xuân Lạc.

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

+-------------------------------------------------------------------------------+
| Sinh viên Lâm nghiệp / Môi trường:                                           |
| -> Tài liệu chuẩn về phương pháp lập OTC, đo đếm sinh trắc và tính IV%, H'.   |
+-------------------------------------------------------------------------------+
| Cán bộ Kiểm lâm & Ban quản lý KBT:                                            |
| -> Bản đồ phân bố 6 trạng thái thảm và danh lục 21 loài nguy cấp cần bảo vệ.  |
+-------------------------------------------------------------------------------+
| Cộng đồng dân cư vùng đệm (986 hộ dân):                                       |
| -> Cơ chế giao khoán rừng và hướng dẫn kỹ thuật canh tác dược liệu dưới tán.  |
+-------------------------------------------------------------------------------+
| Nhà nghiên cứu Sinh thái học Karst:                                           |
| -> Dữ liệu đối sánh cấu trúc tổ thành và đa dạng loài hệ sinh thái Đông Bắc.  |
+-------------------------------------------------------------------------------+

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

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

Cần thiết lập OTC diện tích chuẩn từ 400 $m^2$ đến 2.500 $m^2$ tùy theo trạng thái rừng, sử dụng thước kẹp kính đo $D_{1.3} \ge 6\text{ cm}$ theo 2 hướng vuông góc, đo chiều cao bằng thước Blume-Leiss và định vị tọa độ bằng GPS đa băng tần có độ nhạy cao.

2. Giới hạn sinh thái của rừng hỗn giao lá rộng - lá kim tại Nam Xuân Lạc là gì?

Kiểu rừng này chỉ phân bố cục bộ ở đai cao $>700\text{ m}$ trên các đỉnh núi đá vôi lởm chởm, nơi có loài Thông Pà Cò (Pinus kwangtungensis) chiếm ưu thế với chỉ số đa dạng $H' = 1.605$, độ tàn che thấp (0.2 - 0.3) và tầng đất mỏng.

3. Phương pháp tích hợp dữ liệu điều tra thực địa vào hệ thống GIS quản lý rừng?

Số liệu tọa độ OTC, thuộc tính tầng tán ($IV%$, $D_{1.3}$, $H_{vn}$) và danh lục loài được cấu trúc hóa dưới dạng bảng dữ liệu không gian (PostGIS/Shapefile) để chồng lớp lên bản đồ số độ cao (DEM) và hiện trạng rừng trên QGIS/ArcGIS.

4. Công tác duy trì và giám sát diễn thế rừng phục hồi sau nương rẫy cần thực hiện như thế nào?

Cần duy trì các dải băng khoanh nuôi xúc tiến tái sinh tự nhiên, loại bỏ dây leo chèn ép cây tái sinh mục đích (Dẻ gai, Sến nạc, Chò đãi), và kiểm tra định kỳ mật độ cây tái sinh trên các ô dạng bản 4 $m^2$.

5. Khung phân bổ ngân sách và thời gian thu hồi vốn sinh thái cho KBT?

Ngân sách bảo tồn được huy động từ Quỹ Bảo vệ và Phát triển rừng (chi trả DVMTR), vốn ngân sách Nhà nước và các chương trình mục tiêu quốc gia. Hiệu quả sinh thái - nguồn nước được xác lập ngay từ năm thứ 2, hiệu quả sinh kế vùng đệm phát huy sau 3 - 5 năm triển khai.

Kết luận

Khóa luận tốt nghiệp của tác giả Ma Doãn Giang đã hoàn thành xuất sắc các mục tiêu nghiên cứu, lượng hóa toàn diện tính đa dạng thực vật thân gỗ tại KBT loài và sinh cảnh Nam Xuân Lạc. Nghiên cứu đã xác lập công thức tổ thành chính xác cho 6 trạng thái thảm thực vật, chứng minh chỉ số đa dạng loài đạt mức rất cao trên núi đá vôi đai cao $>700\text{ m}$ ($H' = 3.593, D_1 = 0.966$), đồng thời phát hiện quần thể Re hương quý hiếm cấp CR. Các giải pháp lâm sinh kết hợp hỗ trợ sinh kế vùng đệm là cơ sở khoa học và thực tiễn vững chắc phục vụ công tác quy hoạch bảo tồn tài nguyên rừng bền vững tại tỉnh Bắc Kạn và vùng núi Karst Đông Bắc Việt Nam.