Giới thiệu dự án

Suy thoái đa dạng sinh học (ĐDSH) và mất rừng tự nhiên đang diễn ra với tốc độ báo động trên quy mô toàn cầu. Theo báo cáo của Tổ chức Lương thực và Nông nghiệp Liên Hợp Quốc (FAO) và Quỹ Quốc tế Bảo vệ Thiên nhiên (WWF), thế giới mất trung bình khoảng 11,3 triệu ha rừng nhiệt đới mỗi năm, kéo theo nguy cơ tuyệt chủng của 20% đến 50% các loài sinh vật liên đới. Tại Việt Nam, mặc dù độ che phủ rừng đã tăng từ 37,3% (năm 2004) lên 38,2% với tổng diện tích 12,837 triệu ha, nhưng thực trạng chất lượng rừng lại ở mức đáng lo ngại: hơn 66% diện tích là rừng nghèo hoặc rừng phục hồi non trẻ, trong khi diện tích rừng giàu nguyên sinh chỉ chiếm từ 3,4% đến 4,6%. Hầu hết các diện tích rừng giàu còn lại đều bị chia cắt, cô lập thành những "ốc đảo" sinh thái phân tán.

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 thành lập theo Quyết định số 342/QĐ-UB ngày 17/03/2004 của UBND tỉnh Bắc Kạn với diện tích 1.788 ha. Đây là hành lang sinh thái trọng yếu nối liền Vườn quốc gia Ba Bể với Khu bảo tồn Thiên nhiên Nà Hang (Tuyên Quang). Mặc dù có tính đa dạng sinh thái cao, song khu vực rừng thường xanh trên núi đất tại đai cao 600 - 800 m đang chịu sức ép nặng nề từ các hoạt động nhân sinh của cộng đồng dân tộc thiểu số vùng đệm (chiếm 89,5% dân số) như khai thác lâm sản trái phép, du canh du cư và chuyển đổi đất rừng.

+-------------------------------------------------------------------------+
|                  KHU BẢO TỒN THIÊN NHIÊN NÀ HANG                        |
+------------------------------------+------------------------------------+
                                     | (Hành lang sinh thái rừng núi đất)
+------------------------------------v------------------------------------+
|       KBT LOÀI VÀ SINH CẢNH NAM XUÂN LẠC (1.788 ha - Đai 600-800m)     |
+------------------------------------+------------------------------------+
                                     | (Kết nối bảo tồn nguồn gen)
+------------------------------------v------------------------------------+
|                      VƯỜN QUỐC GIA BA BỂ                                |
+-------------------------------------------------------------------------+

Vấn đề nghiên cứu (Problem Statement)

Thiếu hụt dữ liệu định lượng có hệ thống về cấu trúc tầng tán, tổ thành loài, hệ số đa dạng toán học và khả năng tái sinh tự nhiên của quần xã thực vật thân gỗ trên lập địa núi đất ở đai cao 600 - 800 m. Việc quản lý hiện nay chủ yếu dựa trên quan sát định tính, dẫn đến việc thiếu các giải pháp khoanh nuôi, bảo tồn có cơ sở dữ liệu xác thực cho các loài cây gỗ quý hiếm có nguy cơ tuyệt chủng.

Mục tiêu dự án

  1. Định lượng cấu trúc tổ thành và hình thái tầng tán: Xác định các chỉ số sinh thái học của lâm phần rừng kín thường xanh mưa ẩm nhiệt đới núi đất tại đai cao 600 - 800 m.
  2. Đánh giá định lượng tính đa dạng thực vật: Thiết lập hệ thống chỉ số đa dạng loài bao gồm Simpson ($D_1$), Shannon-Wiener ($H'$), Margalef ($d_1$) và Menhinick ($d_2$).
  3. Thống kê phân loại học và giá trị sử dụng: Thống kê taxon theo hệ thống Brummitt (1992) và phân loại giá trị kinh tế/sinh thái (lấy gỗ, làm thuốc, cây cảnh, lấy quả).
  4. Phân tích động thái tái sinh tự nhiên: Đo đếm mật độ, nguồn gốc và chất lượng cây con tái sinh, xác định các loài có giá trị bảo tồn ưu tiên (Sách Đỏ Việt Nam, Danh lục Đỏ IUCN, Nghị định 32/2006/NĐ-CP).
  5. Đề xuất chiến lược bảo tồn ứng dụng: Xây dựng mô hình can thiệp lâm sinh và quản lý dựa vào cộng đồng cho KBT Nam Xuân Lạc.

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

  • Phương pháp tiếp cận: Ứng dụng điều tra sinh thái rừng chuẩn hóa qua hệ thống ô tiêu chuẩn (OTC) điển hình kết hợp thuật toán thống kê sinh học (Biostatistics) để xử lý dữ liệu lâm học.
  • Phạm vi nghiên cứu: Tầng cây cao ($H > 7\text{ m}$, $D_{1.3} \ge 6\text{ cm}$) và tầng cây tái sinh thân gỗ tại đai cao 600 - 800 m trên núi đất thuộc KBT Nam Xuân Lạc. Các dạng sống thân thảo, dây leo, tre nứa không thuộc phạm vi định lượng chính.

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

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

Tiêu chí Phương pháp điều tra truyền thống (Định tính) Giải pháp nghiên cứu sinh thái định lượng (Đề tài)
Cơ chế thu thập Đi tuyến quan sát, ghi chép mô tả trực quan Thiết lập 07 OTC ($500\text{ m}^2/\text{OTC}$) kết hợp 35 ô dạng bản (ODB)
Chỉ số cấu trúc Ước lượng độ tàn che, ghi nhận loài bắt gặp Tính chỉ số quan trọng ($IV%$), tiết diện ngang ($G$), mật độ ($N/\text{ha}$)
Đánh giá đa dạng Đếm số lượng loài tổng thể đơn thuần Phân tích tổ hợp đa dạng ($D_1, H', d_1, d_2$) phản ánh cả độ phong phú và độ đồng đều
Tái sinh rừng Đánh giá cảm quan tốt/xấu Định lượng hệ số tổ thành tái sinh ($K_i$), phân cấp chất lượng, nguồn gốc chồi/hạt
Khả năng nhân rộng Thấp, phụ thuộc kinh nghiệm chuyên gia Rất cao, chuẩn hóa dạng số liệu, dễ tích hợp vào GIS và CSDL lâm nghiệp

Ma trận ưu tiên yêu cầu nghiên cứu (MoSCoW)

  • Must-have: Xác định chính xác công thức tổ thành tầng cây cao ($IV% \ge 5%$), đo đếm $D_{1.3}$, $H_{vn}$, $H_{dc}$, xác lập danh lục taxon hạt kín.
  • Should-have: Tính toán 4 chỉ số đa dạng toán sinh học, đánh giá mật độ và chất lượng tầng cây tái sinh trên 35 ODB.
  • Could-have: Định vị tọa độ GPS vi mô cho từng cá thể loài quý hiếm (Anoectochilus calcareus, Cycas balansae).
  • Won't-have: Giải trình tự gen DNA phân tử, đo đếm sinh khối tươi/khô tuyệt đối ngoài hiện trường (do hạn chế thời gian và thiết bị phá hủy mẫu).

Thiết kế hệ thống thu thập và xử lý dữ liệu

+----------------------------------------------------------------------------+
|                         THU THẬP HIỆN TRƯỜNG                               |
|  - Định vị GPS vệ tinh 7 OTC (500m2) ở đai cao 600-800m                     |
|  - Đo đếm cây gỗ (D1.3, Hvn, Hdc, Dt) + Lập 5 ODB (25m2) điều tra tái sinh|
+-------------------------------------+--------------------------------------+
                                      |
                                      v
+----------------------------------------------------------------------------+
|                          XỬ LÝ DỮ LIỆU LÂM HỌC                             |
|  - Chuẩn hóa danh pháp khoa học (Brummitt 1992, Tokyo Code 1994)            |
|  - Tính toán Tiết diện ngang G (m2/ha), Độ ưu thế sinh thái IV%             |
|  - Tính toán chỉ số đa dạng (Simpson D1, Shannon H', Margalef d1, Menhinick)|
+-------------------------------------+--------------------------------------+
                                      |
                                      v
+----------------------------------------------------------------------------+
|                        PHÂN TÍCH VÀ RA QUYẾT ĐỊNH                          |
|  - Xác lập nhóm loài ưu thế lâm phần                                       |
|  - Dự báo động thái suy thoái / phục hồi qua tầng cây con                  |
|  - Phân vùng ưu tiên bảo tồn nguồn gen (Sách Đỏ, IUCN, NĐ 32/2006)        |
+----------------------------------------------------------------------------+

Ngăn xếp công cụ và phiên bản kỹ thuật (Tech Stack)

  • Hệ thống phân loại thực vật: Hệ thống Brummitt (1992), Luật danh pháp quốc tế Tokyo (1994), Hệ thống phân loại thảm thực vật Thái Văn Trừng (1978).
  • Công cụ đo đạc thực địa: Thước đo cao Blume-Leiss (độ chính xác $\pm 0.1\text{ dm}$), Thước kẹp kính cơ khí (độ chính xác $\pm 1\text{ mm}$), Máy định vị vệ tinh GPS Garmin GPSMAP 64s, Thước dây thép bọc nhựa.
  • Nền tảng phân tích số liệu: Microsoft Excel 2013 / Office Analysis Toolpak, Scripting thống kê sinh học Python 3.10 (thư viện numpy, pandas, scipy).

Cấu trúc cơ sở dữ liệu mẫu điều tra (Data Schema)

[OTC_Plots]
  - plot_id (PK): INT (1 to 7)
  - elevation: INT (627m - 786m)
  - slope: FLOAT (25 - 45 deg)
  - coordinates: VARCHAR (Lat/Long)
       |
       |-- 1:N --< [Tree_Inventory]
       |             - tree_id (PK): INT
       |             - plot_id (FK): INT
       |             - species_name: VARCHAR
       |             - d_1_3: FLOAT (cm)
       |             - h_vn: FLOAT (m)
       |             - h_dc: FLOAT (m)
       |             - d_crown: FLOAT (m)
       |
       |-- 1:N --< [Regeneration_ODB]
                     - odb_id (PK): INT
                     - plot_id (FK): INT
                     - species_name: VARCHAR
                     - height: FLOAT (cm)
                     - quality_grade: ENUM ('A - Tot', 'B - Trung binh', 'C - Xau')
                     - origin: ENUM ('Hat', 'Choi')

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

Quy trình thực hiện được chia làm 3 giai đoạn:

  1. Khảo sát tuyến tiền trạm: Lập các tuyến điều tra cắt ngang và song song với đường đồng mức, bao quát dải độ cao từ 600 m đến 800 m. Băng quan sát có chiều rộng $20 - 40\text{ m}$ (mỗi bên $10 - 20\text{ m}$).
  2. Thiết lập ô tiêu chuẩn điển hình tạm thời: Bố trí 07 OTC có diện tích $500\text{ m}^2$ ($20\text{ m} \times 25\text{ m}$) tại các tọa độ và độ cao đại diện.
  3. Điều tra cấu trúc chi tiết:
    • Tầng cây gỗ: Đo toàn bộ cây có chiều cao $H > 7\text{ m}$, đường kính $D_{1.3} \ge 6\text{ cm}$.
    • Tầng tái sinh: Thiết lập 05 ô dạng bản (ODB) diện tích $5\text{ m}^2$ tại 4 góc và tâm của mỗi OTC (tổng diện tích điều tra tái sinh: $25\text{ m}^2/\text{OTC}$).
    • Tầng cây bụi và thảm tươi: Điều tra trên 05 ODB diện tích $25\text{ m}^2$ ($5\text{ m} \times 5\text{ m}$) phân bố tại các góc và tâm OTC.

Implementation và kết quả

Quy trình phân tích số liệu và thuật toán

Toàn bộ số liệu thô từ phiếu điều tra được số hóa và tính toán qua các công thức toán sinh thái học:

  1. Độ ưu thế sinh thái (Importance Value Index - $IV%$): $$IV% = \frac{N_i% + G_i%}{2}$$ Trong đó:

    • $N_i% = \frac{n_i}{N} \times 100$: Tỷ lệ phần trăm số cá thể của loài $i$ so với tổng số cá thể trong OTC.
    • $G_i% = \frac{g_i}{\sum G} \times 100$: Tỷ lệ phần trăm tiết diện ngang của loài $i$. Tiết diện ngang cá thể: $g = \frac{\pi \times D_{1.3}^2}{40.000}$ ($\text{m}^2$).
    • Nhóm loài có $IV% \ge 5%$ được đưa vào công thức tổ thành; dừng gom nhóm khi tổng $IV% \ge 50%$.

  2. Các chỉ số đa dạng sinh thái loài:

    • Chỉ số Simpson ($D_1$): $D_1 = 1 - \sum_{i=1}^{S} P_i^2$
    • Chỉ số Shannon-Wiener ($H'$): $H' = -\sum_{i=1}^{S} P_i \ln(P_i)$
    • Chỉ số Margalef ($d_1$): $d_1 = \frac{S - 1}{\ln(N)}$
    • Chỉ số Menhinick ($d_2$): $d_2 = \frac{S}{\sqrt{N}}$ (Với $S$ là tổng số loài, $N$ là tổng số cá thể, $P_i = n_i / N$ là tỷ lệ cá thể loài $i$)
import numpy as np
import pandas as pd

def calculate_ecological_metrics(tree_data: pd.DataFrame):
    """
    Tính toán tiết diện ngang G, IV%, Simpson (D1), Shannon (H'), 
    Margalef (d1) và Menhinick (d2) cho lâm phần.
    """
    # 1. Tính tiết diện ngang g (m2) cho từng cây
    tree_data['g_m2'] = (np.pi * (tree_data['d_1_3'] ** 2)) / 40000.0
    
    total_trees = len(tree_data)
    total_g = tree_data['g_m2'].sum()
    
    # 2. Tổng hợp theo loài
    species_group = tree_data.groupby('species_name').agg(
        ni=('tree_id', 'count'),
        gi=('g_m2', 'sum')
    ).reset_index()
    
    # 3. Tính Ni%, Gi% và IV%
    species_group['Ni_pct'] = (species_group['ni'] / total_trees) * 100.0
    species_group['Gi_pct'] = (species_group['gi'] / total_g) * 100.0
    species_group['IV_pct'] = (species_group['Ni_pct'] + species_group['Gi_pct']) / 2.0
    
    # Sắp xếp theo IV% giảm dần
    species_group = species_group.sort_values(by='IV_pct', ascending=False)
    
    # 4. Tính toán các chỉ số đa dạng
    S = len(species_group)
    N = total_trees
    pi = species_group['ni'] / N
    
    simpson_d1 = 1.0 - np.sum(pi ** 2)
    shannon_h = -np.sum(pi * np.log(pi))
    margalef_d1 = (S - 1) / np.log(N) if N > 1 else 0
    menhinick_d2 = S / np.sqrt(N) if N > 0 else 0
    
    return {
        "species_summary": species_group,
        "metrics": {
            "S": S,
            "N": N,
            "Total_G_m2": total_g,
            "Simpson_D1": simpson_d1,
            "Shannon_H": shannon_h,
            "Margalef_d1": margalef_d1,
            "Menhinick_d2": menhinick_d2
        }
    }

Kết quả nghiên cứu thực nghiệm

1. Cấu trúc tầng tán lâm phần

Rừng tự nhiên trên núi đất đai cao 600 - 800 m tại KBT Nam Xuân Lạc duy trì cấu trúc 4 tầng điển hình của rừng mưa nhiệt đới nguyên sinh/ít bị tác động:

  • Tầng vượt tán và tán chính (A1/A2): Chiều cao vút ngọn $H_{vn} = 20 - 30\text{ m}$, đường kính ngang ngực $D_{1.3} = 40 - 50\text{ cm}$, mật độ $60 - 100\text{ cây/ha}$. Đại diện ưu thế: Muồng trắng (Zenia insignis), Kháo (Machilus platycarpa, Machilus thunbergii), Mọ (Deutzianthus tonkinensis), Xoan nhừ (Choerospondias axillaris), Dẻ (Castanopsis spp.), Sồi (Quercus spp., Lithocarpus spp.).
  • Tầng dưới tán (A3): Chiều cao $H_{vn} = 8 - 15\text{ m}$. Loài ưu thế: Phân mã (Archidendron chevalieri), Nhãn rừng (Dimocarpus fumatus), Trai đỏ (Garcinia bracteata), Mạy tèo (Streblus macrophyllus).
  • Tầng cây bụi (B): Chiều cao $3 - 5\text{ m}$. Gồm Lấu núi (Psychotria montana), Cam gai (Zanthoxylum avicenniae), Găng Việt Nam (Rothmannia vietnamensis).
  • Tầng thảm tươi và dây leo (C/D): Gồm các loài họ Gừng (Zingiberaceae), họ Ráy (Araceae), Dương xỉ (Adiantaceae, Dryopteridaceae). Đặc biệt xuất hiện các loài quý hiếm: Lan kim tuyến đá vôi (Anoectochilus calcareus), Hoàng tinh cách (Disporopsis longifolia), Thiên tuế (Cycas balansae), Trân châu chen (Lysimachia chenii).

2. Cấu trúc tổ thành thực vật tầng cây gỗ

+-------------------------------------------------------------------------+
|                  CÔNG THỨC TỔ THÀNH CHUNG LÂM PHẦN                     |
|                                                                         |
|   16.10 Cl + 15.27 Mt + 8.45 Chc + 6.82 Dg + 53.63 Lk                   |
|                                                                         |
|   - Cl:  Cà lồ (16.10%)          - Dg: Dẻ gai (6.82%)                   |
|   - Mt:  Muồng trắng (15.27%)    - Lk: Loài khác (53.63%)               |
|   - Chc: Chân chim (8.45%)                                              |
+-------------------------------------------------------------------------+

Chi tiết công thức tổ thành tại 07 OTC khảo sát:

  • OTC 1 (695 m): $37.73\text{Chc} + 17.59\text{Trk} + 14.96\text{Mt} + 6.83\text{Ngb} + 22.89\text{Lk}$
  • OTC 2 (699 m): $54.26\text{Mt} + 9.11\text{Chc} + 8.83\text{Sn} + 6.38\text{Dv} + 5.05\text{Nc} + 16.37\text{Lk}$
  • OTC 3 (700 m): $15.90\text{Bn} + 9.49\text{Ldđ} + 9.38\text{Ngl} + 7.90\text{Sv} + 6.62\text{Tb} + 6.38\text{Du} + 5.99\text{Cn} + 5.29\text{Dg} + 33.05\text{Lk}$
  • OTC 4 (756 m): $38.47\text{Dg} + 8.61\text{Phm} + 7.82\text{Xnh} + 7.54\text{Bn} + 7.00\text{Mld} + 6.53\text{Thb} + 24.03\text{Lk}$
  • OTC 5 (786 m): $12.32\text{Thb} + 10.70\text{Xn} + 9.41\text{Sog} + 8.71\text{Gt} + 7.91\text{Trat} + 6.38\text{Dh} + 44.57\text{Lk}$
  • OTC 6 (627 m): $41.39\text{Cl} + 14.08\text{Sg} + 13.62\text{Mt} + 6.93\text{Dg} + 6.31\text{Trt} + 17.67\text{Lk}$
  • OTC 7 (643 m): $41.05\text{Cl} + 10.88\text{Sg} + 7.01\text{Trt} + 5.55\text{Kn} + 5.21\text{Mt} + 30.30\text{Lk}$

3. Bộ chỉ số đa dạng sinh thái học qua 07 OTC

OTC Độ cao (m) Simpson ($D_1$) Margalef ($d_1$) Menhinick ($d_2$) Shannon-Wiener ($H'$) Đánh giá mức độ đa dạng
OTC 1 695 0.74 7.6135 1.8058 1.93 Trung bình
OTC 2 699 0.84 7.3436 2.2937 2.12 Khá
OTC 3 700 0.91 12.8827 3.0729 2.81 Rất cao
OTC 4 756 0.84 8.0443 2.0226 2.27 Khá
OTC 5 786 0.94 11.7125 2.9406 2.98 Cao nhất
OTC 6 627 0.85 6.2851 1.7614 2.12 Khá
OTC 7 643 0.88 9.1414 2.5725 2.40 Khá

Nhận xét: Chỉ số $H'$ đạt giá trị cao nhất tại OTC 5 ($H' = 2.98$) và OTC 3 ($H' = 2.81$), chỉ số Simpson dao động từ 0.74 đến 0.94 chứng minh độ phong phú loài lớn, không bị ưu thế tuyệt đối bởi một loài đơn lẻ, phản ánh trạng thái rừng nguyên sinh ổn định cao.

4. Thống kê Taxon và giá trị sử dụng

Tổng điều tra đã ghi nhận 72 loài thực vật thân gỗ, thuộc 30 họ, 100% thuộc Lớp Hai lá mầm (Magnoliopsida), Ngành Ngọc lan (Magnoliophyta).

TỔNG HỢP GIÁ TRỊ SỬ DỤNG CỦA 72 LOÀI THỰC VẬT THÂN GỖ:
+---------------------------------------------------------------+
| Cây lấy gỗ:       61 loài (84.72%) [████████████████████░░░░] |
| Cây làm thuốc:     7 loài ( 9.72%) [██░░░░░░░░░░░░░░░░░░░░] |
| Cây làm cảnh:      2 loài ( 2.78%) [█░░░░░░░░░░░░░░░░░░░░░] |
| Cây ăn quả:        2 loài ( 2.78%) [█░░░░░░░░░░░░░░░░░░░░░] |
+---------------------------------------------------------------+

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

Đổi mới phương pháp luận

Khóa luận đã chuyển hóa phương thức đánh giá rừng từ quan sát cảm tính truyền thống sang mô hình toán sinh học định lượng chính xác:

  • Tích hợp đa chỉ số: Kết hợp đồng thời 4 chỉ số $D_1, H', d_1, d_2$ giúp loại bỏ sai số lấy mẫu do kích thước quần thể khác biệt.
  • Phân tầng sinh thái chi tiết: Lần đầu tiên phân lập rõ ranh giới sinh thái đai cao 600 - 800 m trên núi đất trong quần thể bảo tồn Nam Xuân Lạc vốn chủ yếu được biết đến với địa hình núi đá vôi.

So sánh với các công trình nghiên cứu trong nước

Tiêu chí KBT Nam Xuân Lạc (Đề tài) VQG Ba Vì (Đai cao) VQG Yok Đôn (Rừng khộp)
Kiểu phụ thảm thực vật Rừng kín thường xanh mưa ẩm nhiệt đới núi đất Rừng kín thường xanh và rừng hỗn hợp lá rộng - lá kim Rừng thưa cây lá rộng rụng lá khô hạn
Số lượng loài cây gỗ 72 loài / 30 họ (trên 7 OTC) 417 loài / 136 họ (toàn bộ đai cao) 566 loài / 108 họ
Loài ưu thế chính Cà lồ ($16.1%$), Muồng trắng ($15.3%$) Re vàng, Vàng tâm, Bời lời Ba Vì Dầu trà beng, Dầu đồng, Cẩm liên
Chỉ số Shannon ($H'$) $1.93 - 2.98$ (Trung bình $2.40$) $2.10 - 2.65$ $1.20 - 1.85$ (Đơn giản tầng thứ)
Độ đa dạng taxon 100% hạt kín hai lá mầm Hạt kín ($90.4%$) + Hạt trần ($9.6%$) Hạt kín ($98.8%$) + Hạt trần ($0.1%$)

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

Kịch bản ứng dụng trong quản trị tài nguyên rừng

  1. Phân vùng bảo vệ nghiêm ngặt (Vùng lõi): Thiết lập tọa độ bảo vệ ưu tiên cho các phân khu tại OTC 3 (700 m) và OTC 5 (786 m) - nơi có chỉ số $H' > 2.80$ và xuất hiện mật độ cao của các loài cần bảo tồn gen như Muồng trắng (Zenia insignis), Trai đỏ (Garcinia bracteata).
  2. Kế hoạch xúc tiến tái sinh tự nhiên: Áp dụng các giải pháp lâm sinh can thiệp nhẹ (phát dây leo, mở tán cục bộ) quanh cây mẹ gieo giống của các loài Cà lồ, Muồng trắng và Dẻ gai để nâng cao tỷ lệ sống của cây mạ.
  3. Mô hình đồng quản trị rừng với cộng đồng vùng đệm: Chuyển giao kỹ thuật nhân giống cây lâm sản ngoài gỗ và cây dược liệu (Anoectochilus calcareus) dưới tán rừng cho 410 hộ dân tại 9 thôn vùng đệm, thay thế tập quán phá rừng làm nương rẫy.
LỘ TRÌNH TRIỂN KHAI ỨNG DỤNG BẢO TỒN (2015 - 2020)

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

Hạn chế kỹ thuật

  • Quy mô lấy mẫu: Nghiên cứu mới dừng lại ở 07 ô tiêu chuẩn ($3.500\text{ m}^2$) đại diện cho đai cao 600 - 800 m, chưa bao phủ toàn bộ diện tích thung lũng xen kẽ núi đá.
  • Chu kỳ thời gian: Dữ liệu điều tra phản ánh trạng thái tĩnh trong giai đoạn 2014 - 2015, chưa theo dõi được biến động mùa vụ sinh học (phenology) qua nhiều năm.
  • Giới hạn taxon: Chưa định lượng sâu các chỉ số sinh thái của lớp thực vật phụ sinh, phong lan và nấm rừng.

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

  • Tích hợp công nghệ viễn thám và LiDAR: Sử dụng ảnh vệ tinh đa phổ độ phân giải cao kết hợp thiết bị bay không người lái (UAV/LiDAR) để lập bản đồ sinh khối và trữ lượng carbon rừng.
  • Giám sát động thái đa dạng sinh học dài hạn: Chuyển đổi 07 OTC thành hệ thống ô định vị sinh thái lâu dài (Permanent Sample Plots - PSP) theo chuẩn mạng lưới CTFS-ForestGEO.
  • Nghiên cứu sinh học phân tử: Ứng dụng mã vạch DNA (DNA barcoding) phục vụ giám định các loài gỗ quý và dược liệu đặc hữu tại Nam Xuân Lạc.

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

+-------------------------------------------------------------------------------+
|                      MA TRẬN LỢI ÍCH ĐỊNH LƯỢNG                                |
+-----------------------+-------------------------------------------------------+
| Sinh viên & Học viên  | Kế thừa bộ dữ liệu mẫu 7 OTC, công thức tính toán     |
| Lâm nghiệp            | và bộ tiêu bản 72 loài làm tài liệu học tập thực tế.  |
+-----------------------+-------------------------------------------------------+
| Cán bộ Kỹ thuật &     | Ứng dụng thuật toán Python/Excel để tự động hóa       |
| Data Analyst          | quy trình tính toán IV% và chỉ số Shannon/Simpson.    |
+-----------------------+-------------------------------------------------------+
| Ban Quản lý KBT       | Sở hữu cơ sở khoa học định lượng phục vụ quy hoạch    |
| Nam Xuân Lạc          | vùng lõi và lập đề án xin tài trợ bảo tồn quốc tế.    |
+-----------------------+-------------------------------------------------------+
| Nhà nghiên cứu        | Cung cấp dữ liệu so sánh chuẩn về hệ sinh thái rừng   |
| Đa dạng sinh học      | thường xanh núi đất trong cảnh quan Karst miền Bắc.   |
+-----------------------+-------------------------------------------------------+

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

1. Yêu cầu kỹ thuật và công cụ cần thiết để nhân rộng mô hình điều tra này là gì?

Cần trang bị thiết bị định vị GPS (độ chính xác $< 3\text{ m}$), thước đo cao Blume-Leiss hoặc Vertex IV, thước kẹp kính cây gỗ chuyên dụng, thước dây $50\text{ m}$. Về nhân lực, một tổ điều tra tiêu chuẩn gồm 3 - 4 người (1 chuyên gia phân loại thực vật, 1 kỹ thuật viên đo đếm, 1 nhân viên ghi chép và 1 người dẫn đường địa phương).

2. Tại sao chỉ số đa dạng Shannon-Wiener ($H'$) ở OTC 3 và OTC 5 lại cao vượt trội?

OTC 3 ($H' = 2.81$) và OTC 5 ($H' = 2.98$) nằm ở đai cao trên $700\text{ m}$ đến $786\text{ m}$, khu vực sườn dốc hiểm trở, xa các khu dân cư thôn Nà Dạ và Bản Khang nên ít chịu tác động chặt hạ của con người. Cấu trúc rừng ở đây giữ được độ khép tán nguyên sinh, cho phép nhiều tầng cây gỗ cùng chung sống mà không có loài nào chiếm ưu thế tuyệt đối áp đảo.

3. Phương pháp tính toán $IV%$ có thể áp dụng cho các kiểu rừng khác không?

Hoàn toàn có thể. Công thức $IV% = (N_i% + G_i%) / 2$ là phương pháp chuẩn mực của sinh thái học rừng quốc tế, áp dụng hiệu quả cho mọi kiểu rừng tự nhiên nhiệt đới hỗn loài, rừng ngập mặn hay rừng rụng lá.

4. Làm thế nào để giải quyết mâu thuẫn giữa sinh kế người dân và bảo tồn rừng tại Nam Xuân Lạc?

Cần kết hợp 3 giải pháp đồng bộ: (1) Giao khoán bảo vệ rừng gắn với chính sách chi trả dịch vụ môi trường rừng; (2) Hỗ trợ chuyển đổi giống cây trồng năng suất cao ở diện tích đất nông nghiệp hiện có (thay vì mở rộng nương rẫy); (3) Phát triển mô hình trồng cây dược liệu dưới tán rừng có kiểm soát.

5. Dữ liệu từ đồ án có thể tích hợp vào các hệ thống GIS như thế nào?

Bảng thuộc tính cây gỗ và chỉ số đa dạng có thể xuất ra định dạng .csv hoặc Geodatabase, liên kết với lớp bản đồ ranh giới tiểu khu của KBT Nam Xuân Lạc thông qua trường khóa plot_id trên phần mềm ArcGIS hoặc QGIS để nội suy không gian (Interpolation IDW/Kriging) phân bố đa dạng sinh học.

Kết luận

Công trình nghiên cứu của tác giả Trương Phúc Cường đã giải quyết trọn vẹn mục tiêu định lượng hóa tính đa dạng và cấu trúc lâm phần kiểu rừng thường xanh trên núi đất đai cao 600 - 800 m tại KBT Loài và Sinh cảnh Nam Xuân Lạc. Kết quả đã xác lập danh lục 72 loài cây gỗ thuộc 30 họ, xây dựng công thức tổ thành chuẩn $16.10\text{Cl} + 15.27\text{Mt} + 8.45\text{Chc} + 6.82\text{Dg} + 53.63\text{Lk}$, đồng thời minh chứng trạng thái rừng có tính đa dạng sinh học cao ($H'$ đạt tới 2.98). Đây là luận cứ khoa học thực tiễn vững chắc phục vụ công tác quy hoạch bảo tồn nguồn gen quý hiếm và quản lý bền vững tài nguyên rừng tại tỉnh Bắc Kạn.

Hãy tham khảo các quy chuẩn điều tra sinh thái định lượng trong đồ án này để xây dựng các dự án giám sát ĐDSH và ứng dụng công nghệ trong quản lý lâm nghiệp hiện đại.