Giới thiệu dự án
Hệ sinh thái rừng trên núi đá vôi (karst forest ecosystem) là một trong những sinh cảnh đặc thù và nhạy cảm nhất trong mạng lưới các khu rừng đặc dụng tại Việt Nam. Theo thống kê của ngành lâm nghiệp, diện tích núi đá vôi chiếm khoảng 1,625 triệu ha (tương đương 5,4% đến 6,1% tổng diện tích đất lâm nghiệp cả nước), nhưng chỉ có 396.200 ha (khoảng 34,4%) còn duy trì được thảm che thực vật có rừng. Nằm trong điểm nóng đa dạng sinh học toàn cầu Indo-Burma, Việt Nam sở hữu mức độ đặc hữu thực vật bậc cao lên đến 20% - 40%, tuy nhiên tốc độ suy giảm diện tích rừng nhiệt đới toàn cầu ở mức 11,3 triệu ha/năm đang tạo ra áp lực cực lớn lên nguồn gen quý hiếm.
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 tổng diện tích tự nhiên 1.788,0 ha. Khu vực này đóng vai trò hành lang sinh thái trọng yếu kết nối KBT Thiên nhiên Na Hang (Tuyên Quang) và Vườn Quốc gia Ba Bể (Bắc Kạn). Tuy nhiên, tài nguyên rừng thân gỗ tại đai cao trên 700 m đang chịu tác động nghiêm trọng bởi các hoạt động nhân sinh: khai thác chọn gỗ quý (Nghiến, Trai lý, Re hương), chăn thả gia súc tự do, và áp lực sinh kế vùng đệm.
+-----------------------------------------------------------------------------------+
| KHU BẢO TỒN LOÀI VÀ SINH CẢNH NAM XUÂN LẠC (1.788,0 ha) |
+-----------------------------------------+-----------------------------------------+
| Phân khu Bảo vệ Nghiêm ngặt | Phân khu Phục hồi Sinh thái |
| - Diện tích: 1.646,0 ha (92,06%) | - Diện tích: 142,0 ha (7,94%) |
| - Rừng nghèo IIIA1: 680,4 ha (38,1%) | - Rừng phục hồi: 77,9 ha |
| - Rừng phục hồi: 449,3 ha | - Rừng hỗn giao: 64,1 ha |
| - Rừng hỗn giao gỗ+nứa: 317,3 ha | |
| - Rừng trên núi đá: 199,0 ha (11,1%) | |
+-----------------------------------------+-----------------------------------------+
Vấn đề nghiên cứu và mục tiêu cụ thể
Thiếu hụt dữ liệu định lượng về cấu trúc tổ thành, chỉ số đa dạng sinh học và động thái tái sinh tự nhiên của tầng cây thân gỗ ở đai cao trên 700 m khiến ban quản lý gặp khó khăn trong việc xây dựng phương án lâm sinh phục hồi và phân bổ nguồn lực bảo tồn. Đề tài tập trung giải quyết các mục tiêu:
- Định lượng cấu trúc tổ thành và độ tàn che: Xác định công thức tổ thành tầng cây gỗ và chỉ số ưu thế sinh thái của các loài chủ yếu tại đai cao >700 m.
- Đo lường các chỉ số đa dạng sinh học: Tính toán chỉ số Simpson ($D_1$), Shannon-Wiener ($H'$), Margalef ($d_1$), và Menhinick ($d_2$) trên các ô tiêu chuẩn đại diện.
- Thống kê phân loại và giá trị tài nguyên: Phân loại taxon (ngành, lớp, họ, loài), xác định cơ cấu giá trị sử dụng (gỗ, dược liệu, cây cảnh, tinh dầu).
- Đánh giá nguồn gen nguy cấp: Định danh và xếp cấp các loài cây gỗ quý hiếm theo Sách đỏ Việt Nam (2007), Danh lục đỏ IUCN và Nghị định 32/2006/NĐ-CP.
- Đánh giá tiềm năng tái sinh tự nhiên: Xác định mật độ, công thức tổ thành cây tái sinh và tỷ lệ cá thể đạt phẩm chất tốt/trung bình làm căn cứ phục hồi sinh thái.
Phạm vi nghiên cứu tập trung vào tầng cây gỗ có chiều cao vút ngọn $H_{vn} \ge 7$ m, đường kính ngang ngực $D_{1.3} \ge 6$ cm và tầng cây tái sinh dưới tán tại 4 khu vực trọng điểm: Bình Trai, Lũng Trang, Lũng Lỳ và Thưa Tèo ở độ cao từ 700 m đến 1.159 m (đỉnh Tam Sao).
Phân tích và thiết kế giải pháp
Phân tích hiện trạng
Các phương pháp điều tra lâm học truyền thống và đánh giá hiện trạng tại các khu bảo tồn lân cận bộc lộ nhiều điểm hạn chế khi áp dụng cho địa hình karst chia cắt mạnh.
| Tiêu chí |
Khảo sát tuyến đơn thuần |
Phương pháp Ô tiêu chuẩn tạm thời (OTC) |
Viễn thám quang học độ phân giải TB |
| Độ chính xác định lượng |
Thấp, chỉ thống kê danh lục |
Rất cao, đo đếm chính xác $D_{1.3}, H_{vn}, G$ |
Trung bình, khó xuyên tầng tán karst |
| Đánh giá tái sinh |
Không xác định được mật độ/ha |
Chi tiết thông qua ô dạng bản (ODB) |
Không khả thi cho tầng mạ/tái sinh |
| Chi phí & độ phức tạp |
Thấp, tiến hành nhanh |
Đòi hỏi nhân lực chuyên môn cao |
Chi phí ảnh cao, xử lý dữ liệu phức tạp |
| Khả năng lặp lại |
Kém |
Tốt, chuẩn hóa theo tọa độ và diện tích |
Phụ thuộc độ phủ mây và chu kỳ vệ tinh |
So sánh tương quan đa dạng sinh học thực vật thân gỗ với các vùng lân cận cho thấy: KBT Nam Xuân Lạc sở hữu mật độ loài nguy cấp trên diện tích cao hơn Vườn Quốc gia Yok Đôn (vốn đặc trưng cho rừng khộp nghèo tầng tán) và mang tính pha trộn cao giữa luồng di cư bản địa Bắc Việt Nam - Nam Trung Hoa với các luồng Mã Lai và Himalaya.
+-------------------------------------------------------------------------------+
| MA TRẬN YÊU CẦU BẢO TỒN VÀ QUẢN LÝ DỮ LIỆU (MoSCoW) |
+-------------------------------------------------------------------------------+
| [M] MUST HAVE (Bắt buộc): |
| - Đo đếm D1.3, Hvn, Hdc, Dt tầng cây gỗ trên 7 OTC (500 m2/ô). |
| - Tính toán IV% (Important Value Index) xác định loài ưu thế sinh thái. |
| - Kiểm kê phân cấp bảo tồn theo Sách đỏ VN, IUCN, NĐ 32/2006/NĐ-CP. |
| [S] SHOULD HAVE (Nên có): |
| - Thiết lập 5 ODB (4 m2/ô) trong mỗi OTC đo đếm mật độ và phẩm chất tái sinh. |
| - Tính toán hệ chỉ số Shannon-Wiener (H') và Simpson (D1). |
| [C] COULD HAVE (Có thể có): |
| - Bản đồ hóa tọa độ phân bố tiểu quần thể Re hương (Cinnamomum parthenoxylon)|
| [W] WON'T HAVE (Chưa thực hiện trong giai đoạn này): |
| - Giải trình tự gen phân tử DNA barcoding cho toàn bộ hệ thực vật thân thảo.|
+-------------------------------------------------------------------------------+
Thiết kế hệ thống thu thập và xử lý dữ liệu
Kiến trúc xử lý dữ liệu điều tra thực địa được chuẩn hóa theo quy trình 4 giai đoạn, tích hợp các công cụ lâm học tiêu chuẩn và thuật toán sinh trắc rừng:
flowchart TD
A["Thiết lập tuyến khảo sát độ cao >700m"] --> B["Bố trí 7 Ô tiêu chuẩn (OTC: 25m x 20m = 500m2)"]
B --> C["Đo đếm Tầng Cây Gỗ (Hvn >= 7m, D1.3)"]
B --> D["Bố trí 5 Ô Dạng Bản (ODB: 2m x 2m = 4m2)"]
C --> E["Đo đếm D1.3 (Thước kẹp kính), Hvn (Blumeleiss), Dt (Thước dây)"]
D --> F["Đo đếm Chiều cao cây tái sinh, Phân cấp chất lượng (Tốt/TB/Xấu)"]
E --> G["Module Xử lý Dữ liệu Lâm học"]
F --> G
G --> H["Tính toán Tiết diện ngang G, IV%, Chỉ số Đa dạng (H', D1, d1, d2)"]
G --> I["Phân cấp Bảo tồn (Sách đỏ VN 2007, IUCN, NĐ 32)"]
H --> J["Báo cáo Cấu trúc & Đề xuất Giải pháp Lâm sinh"]
I --> J
Thiết kế Schema dữ liệu lâm học (Relational Schema)
-- Bảng dữ liệu ô tiêu chuẩn (Plot_Metadata)
CREATE TABLE Plot_Metadata (
plot_id VARCHAR(10) PRIMARY KEY,
location_name VARCHAR(100) NOT NULL,
elevation_m INT CHECK (elevation_m >= 700),
slope_degree INT CHECK (slope_degree BETWEEN 0 AND 90),
area_m2 FLOAT DEFAULT 500.0,
canopy_cover FLOAT CHECK (canopy_cover BETWEEN 0 AND 1.0)
);
-- Bảng kiểm kê tầng cây gỗ (Woody_Trees)
CREATE TABLE Woody_Trees (
tree_id INT PRIMARY KEY AUTO_INCREMENT,
plot_id VARCHAR(10),
species_scientific VARCHAR(100) NOT NULL,
species_vietnamese VARCHAR(100) NOT NULL,
dbh_cm FLOAT NOT NULL CHECK (dbh_cm > 0),
h_vn_m FLOAT NOT NULL CHECK (h_vn_m >= 7.0),
h_dc_m FLOAT,
crown_d_m FLOAT,
FOREIGN KEY (plot_id) REFERENCES Plot_Metadata(plot_id)
);
-- Bảng cây tái sinh (Regeneration_Data)
CREATE TABLE Regeneration_Data (
regen_id INT PRIMARY KEY AUTO_INCREMENT,
plot_id VARCHAR(10),
species_name VARCHAR(100) NOT NULL,
height_cm FLOAT NOT NULL,
quality_grade ENUM('Good', 'Average', 'Poor') NOT NULL,
origin ENUM('Seed', 'Sprout') DEFAULT 'Seed',
FOREIGN KEY (plot_id) REFERENCES Plot_Metadata(plot_id)
);
Phương pháp nghiên cứu và công thức sinh trắc
Áp dụng 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ương pháp tính toán giá trị quan trọng theo Daniel Marmillod:
-
Tiết diện ngang thân cây ($G$):
$$G = \sum \frac{\pi \cdot D_{1.3}^2}{4 \times 10.000} \quad (\text{m}^2/\text{ha})$$
-
Chỉ số giá trị quan trọng ($IV%$):
$$IV_i% = \frac{N_i% + G_i%}{2}$$
Trong đó: $N_i%$ là tỷ lệ phần trăm số cá thể của loài $i$; $G_i%$ là tỷ lệ phần trăm tiết diện ngang của loài $i$ so với tổng tiết diện ngang toàn ô.
-
Hệ chỉ số đa dạng sinh thái:
- Shannon-Wiener ($H'$): $H' = -\sum_{i=1}^S P_i \ln(P_i) \quad \text{với } P_i = \frac{n_i}{N}$
- Simpson ($D_1$): $D_1 = 1 - \sum_{i=1}^S P_i^2$
- Margalef ($d_1$): $d_1 = \frac{S - 1}{\ln(N)}$
- Menhinick ($d_2$): $d_2 = \frac{S}{\sqrt{N}}$
Implementation và kết quả
Quy trình triển khai thực địa và xử lý dữ liệu
Khóa luận được triển khai qua 4 giai đoạn từ ngày 15/08/2014 đến ngày 30/11/2014:
- Giai đoạn 1 (15/08 - 31/08/2014): Kế thừa tài liệu thứ cấp, bản đồ phân khu chức năng KBT Nam Xuân Lạc, khảo sát tuyến sơ thám tại đai cao >700 m.
- Giai đoạn 2 (01/09 - 15/10/2014): Lập 7 OTC ($20 \times 25$ m, diện tích $500\text{ m}^2$) trên các sườn dốc $25^\circ - 45^\circ$. Bố trí 35 ODB ($2 \times 2$ m, $4\text{ m}^2$). Đo đếm $D_{1.3}$ bằng thước kẹp kính (độ chính xác mm), $H_{vn}$ bằng thước Blumeleiss (độ chính xác dm).
- Giai đoạn 3 (16/10 - 10/11/2014): Định danh mẫu thực vật tại phòng thí nghiệm thực vật rừng Đại học Nông Lâm Thái Nguyên, đối soát Danh lục đỏ IUCN (2014) và Sách đỏ Việt Nam (2007).
- Giai đoạn 4 (11/11 - 30/11/2014): Xử lý ma trận dữ liệu sinh thái trên Microsoft Excel và tính toán thống kê.
Kịch bản Python tự động tính toán cấu trúc tổ thành và chỉ số đa dạng sinh thái
import numpy as np
import pandas as pd
def calculate_ecological_metrics(data_trees: pd.DataFrame) -> dict:
"""
Tính toán Tiết diện ngang G, IV%, Shannon-Wiener (H'), Simpson (D1),
Margalef (d1), Menhinick (d2) cho lâm phần khảo sát.
"""
# 1. Tính tiết diện ngang g_i (m2) = pi * D^2 / (4 * 10000)
data_trees['basal_area_m2'] = (np.pi * (data_trees['dbh_cm'] ** 2)) / 40000.0
total_trees = len(data_trees)
total_basal_area = data_trees['basal_area_m2'].sum()
# 2. Tổng hợp theo từng loài
species_summary = data_trees.groupby('species_scientific').agg(
ni=('dbh_cm', 'count'),
gi=('basal_area_m2', 'sum')
).reset_index()
# 3. Tính Ni%, Gi%, IV%
species_summary['Ni_percent'] = (species_summary['ni'] / total_trees) * 100.0
species_summary['Gi_percent'] = (species_summary['gi'] / total_basal_area) * 100.0
species_summary['IV_percent'] = (species_summary['Ni_percent'] + species_summary['Gi_percent']) / 2.0
# 4. Tính toán các chỉ số đa dạng
S = len(species_summary)
N = total_trees
pi = species_summary['ni'] / N
shannon_h = -np.sum(pi * np.log(pi))
simpson_d1 = 1.0 - np.sum(pi ** 2)
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_table": species_summary.sort_values(by='IV_percent', ascending=False),
"metrics": {
"S_total_species": S,
"N_total_individuals": N,
"Shannon_H": round(shannon_h, 3),
"Simpson_D1": round(simpson_d1, 3),
"Margalef_d1": round(margalef_d1, 3),
"Menhinick_d2": round(menhinick_d2, 3)
}
}
Kết quả cấu trúc tổ thành tầng cây gỗ
Kết quả phân tích 7 OTC đại diện tại đai cao trên 700 m phản ánh sự biến đổi mạnh mẽ của cấu trúc ưu thế theo độ cao và vi địa hình karst.
| OTC |
Độ cao / Địa điểm |
Công thức tổ thành tầng cây gỗ ($IV% \ge 5%$) |
Tổng tỷ lệ $IV%$ nhóm loài ưu thế |
| 01 |
~720m (Bình Trai) |
$34,48\text{Tr} + 21,86\text{Hv} + 12,42\text{N} + 9,30\text{Tbl} + 6,74\text{Thr} + 15,20\text{Lk}$ |
84,80% |
| 02 |
~750m (Lũng Trang) |
$19,97\text{K} + 13,66\text{Ddx} + 12,01\text{Cl} + 10,72\text{Cđ} + 7,23\text{S} + 7,14\text{Va} + 5,17\text{Tbl} + 25,41\text{Lk}$ |
74,59% |
| 03 |
~780m (Lũng Lỳ) |
$29,71\text{Kln} + 27,15\text{Trđ} + 8,11\text{Sb} + 6,74\text{Tbl} + 6,73\text{Cc} + 6,25\text{Tmt} + 15,33\text{Lk}$ |
84,69% |
| 04 |
~820m (Thưa Tèo) |
$46,83\text{Rh} + 18,33\text{Ôr} + 12,07\text{Nr} + 9,26\text{Ng} + 5,24\text{Tđ} + 8,27\text{Lk}$ |
91,73% |
| 05 |
~850m (Đỉnh đá vôi) |
$44,08\text{Sh} + 11,23\text{Tc} + 9,70\text{Tmt} + 6,73\text{Tbl} + 6,73\text{Ql} + 6,37\text{Gxhk} + 15,15\text{Lk}$ |
84,85% |
| 06 |
~895m (Sườn dốc) |
$29,36\text{Tl} + 9,57\text{Ng} + 9,38\text{Tm} + 7,75\text{Xn} + 6,31\text{S} + 37,63\text{Lk}$ |
62,37% |
| 07 |
921m (Đỉnh cao) |
$50,87\text{Ng} + 10,34\text{Nr} + 9,09\text{Gvn} + 8,54\text{N} + 7,46\text{Tl} + 13,72\text{Lk}$ |
86,28% |
(Ký hiệu: Tr: Trâm, Hv: Han voi, N: Nhọc, Tbl: Thôi ba lá, Thr: Thị rừng, K: Kháo, Ddx: Dâu da xoan, Cl: Cà lồ, Cđ: Chò đãi, S: Sung, Va: Vàng anh, Kln: Kháo lá nhỏ, Trđ: Trai đỏ, Sb: Sồi bang, Cc: Chân chim, Tmt: Thổ mật tù, Rh: Re hương, Ôr: Ô rô, Nr: Nhãn rừng, Ng: Nghiến, Tđ: Trai đỏ, Sh: Sếu hôi, Tc: Táo cong, Ql: Quế lợn, Gxhk: Giọt xành hồng, Tl: Trai lý, Tm: Táu muối, Xn: Xoan nhừ, Gvn: Găng Việt Nam, Lk: Loài khác).
+-------------------------------------------------------------------------------+
| BẢNG CHỈ SỐ ĐA DẠNG SINH THÁI TẠI CÁC Ô TIÊU CHUẨN (ĐAI CAO > 700 M) |
+-------+--------------------+----------------+----------------+----------------+
| ÔTC | Simpson (D1) | Margalef (d1) | Menhinick (d2) | Shannon-Wiener |
| | | | | (H') |
+-------+--------------------+----------------+----------------+----------------+
| OTC 01| 0,83 | 6,14 | 40,24 | 1,95 |
| OTC 02| 0,92 | 11,24 | 90,00 | 2,74 |
| OTC 03| 0,82 | 7,30 | 2,40 | 2,04 |
| OTC 04| 0,81 | 5,07 | 39,19 | 1,83 |
| OTC 05| 0,98 | 6,80 | 2,18 | 1,07 |
| OTC 06| 0,92 | 12,02 | 72,33 | 2,75 |
| OTC 07| 0,85 | 6,52 | 64,14 | 2,10 |
+-------+--------------------+----------------+----------------+----------------+
Đánh giá tài nguyên thực vật và động thái tái sinh
- Thống kê phân loại Taxon: Ghi nhận toàn cảnh khu bảo tồn có 1 ngành (Magnoliophyta), 3 lớp, 73 họ và 341 loài. Riêng lớp Hai lá mầm (Dicotyledones) chiếm ưu thế áp đảo với 30 họ (81,08% tổng số họ) và 39 loài thân gỗ chủ đạo chiếm 11,43% thành phần loài cây gỗ toàn khu.
- Giá trị sử dụng: Trong 39 loài thân gỗ đặc trưng: nhóm Cho gỗ chiếm 27 loài (69,23%), nhóm Làm thuốc 10 loài (25,64%), nhóm Cho quả 6 loài (15,38%), nhóm Làm cảnh 4 loài (10,25%), và nhóm Lấy dầu/nhựa 2 loài (5,12%).
- Động thái tầng cây tái sinh: Tỷ lệ cây tái sinh có chất lượng Tốt và Trung bình dao động từ 48,0% đến 77,6% (bình quân toàn lâm phần đạt 62,34%). Tầng cây mạ và cây non của các loài ưu thế (Nghiến, Trai lý, Găng Việt Nam, Kháo) có khả năng sinh trưởng cạnh tranh tốt qua tầng thảm tươi, thâm nhập vào tầng tán rừng trong tương lai.
+---------------------------------------------------------------------------------------+
| DANH MỤC VÀ PHÂN CẤP BẢO TỒN CÁC LOÀI CÂY GỖ NGUY CẤP TẠI ĐAI CAO > 700 M |
+----+---------------------+--------------------------+----------+----------+-----------+
| TT | Tên Việt Nam | Tên khoa học | Sách đỏ | IUCN | Nghị định |
| | | | VN (2007)| (2014) | 32/2006 |
+----+---------------------+--------------------------+----------+----------+-----------+
| 1 | Re hương (Vù hương) | Cinnamomum parthenoxylon | CR | DD | Nhóm IIA |
| 2 | Kháo xanh | Machilus odoratissima | EN | - | - |
| 3 | Nghiến | Excentrodendron tonkin. | EN | EN | Nhóm IIA |
| 4 | Trai lý | Garcinia fagraeoides | EN | EN | - |
| 5 | Sến mật | Madhuca pasquieri | EN | EN | - |
| 6 | Chò đãi | Annamocarya sinensis | VU | LR/NT | - |
| 7 | Chò nâu | Dipterocarpus retusus | VU | VU | - |
| 8 | Chò chỉ | Parashorea chinensis | VU | VU | - |
| 9 | Kim giao | Nageia fleuryi | - | LR/NT | Nhóm IIA |
| 10 | Giổi lông | Michelia hedyosperma | - | LR/LC | - |
| 11 | Sồi phảng | Lithocarpus fenestratus | - | LR/LC | - |
| 12 | Thung | Tetrameles nudiflora | - | - | - |
| 13 | Han voi | Laportea sp. | - | - | - |
+----+---------------------+--------------------------+----------+----------+-----------+
Đổi mới và đóng góp
-
Phát hiện và định lượng vi quần thể Re hương (Cinnamomum parthenoxylon):
Khóa luận đã phát hiện và khoanh vùng quần thể Re hương độc nhất vô nhị tại khu vực Thưa Tèo ở độ cao >820 m. Quần thể gồm 15 - 20 cá thể cổ thụ có đường kính $D_{1.3}$ từ 50 - 80 cm, chiều cao vút ngọn $H_{vn}$ đạt 25 - 30 m, chiếm lĩnh tầng vượt tán cùng với Nghiến và Lát khét ($IV% = 46,83%$). Loài này được xếp cấp Cực kỳ nguy cấp (CR) trong Sách đỏ Việt Nam 2007.
-
Chuyển đổi phương pháp luận từ định tính sang mô hình toán sinh thái:
Thay vì chỉ kiểm kê danh lục đơn thuần, nghiên cứu đã áp dụng chỉ số tổ thành kết hợp Daniel Marmillod $IV% = (N_i% + G_i%)/2$, phản ánh chính xác cả mật độ cá thể lẫn không gian dinh dưỡng thông qua tiết diện ngang thân cây ($G$).
-
Cung cấp cơ sở định lượng cho tái sinh rừng núi đá:
Khẳng định năng lực tái sinh tự nhiên đạt tỷ lệ triển vọng 62,34% (chất lượng Tốt và TB). Điều này chứng minh hệ sinh thái rừng karst >700 m tại Nam Xuân Lạc vẫn duy trì được cân bằng sinh thái nội tại và có khả năng tự phục hồi cao nếu loại bỏ được các xung đột nhân sinh.
Ứ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
- Phân vùng bảo vệ nghiêm ngặt vi điểm (Micro-Zoning): Thiết lập trạm tuần tra chuyên trách bảo vệ nghiêm ngặt quần thể Re hương Thưa Tèo và dải rừng Nghiến - Trai lý tại độ cao từ 800 m đến 921 m (OTC 04, 06, 07).
- Khoanh nuôi xúc tiến tái sinh tự nhiên (ANR - Assisted Natural Regeneration): Ứng dụng công thức tổ thành cây tái sinh để tác động các biện pháp kỹ thuật lâm sinh: phát dọn dây leo, cây bụi chèn ép quanh các cây tái sinh mục đích (Nghiến, Kháo, Trai đỏ) tại các ô có mật độ cây phẩm chất xấu cao.
+-------------------------------------------------------------------------------+
| LỘ TRÌNH TRIỂN KHAI GIẢI PHÁP LÂM SINH VÀ BẢO TỒN (2015 - 2020) |
+-------------------------------------------------------------------------------+
| Giai đoạn 1 (Năm 1): Phân định ranh giới & Định vị GPS vi quần thể Re hương. |
| Giai đoạn 2 (Năm 2-3): Triển khai giao khoán bảo vệ rừng gắn với chi trả DVMTR|
| Giai đoạn 3 (Năm 4-5): Thiết lập 10 ô tiêu chuẩn định vị vĩnh cửu (PSP) |
| theo dõi động thái diễn thế thảm thực vật đai cao. |
+-------------------------------------------------------------------------------+
Phân tích chi phí - hiệu quả (Cost-Benefit Analysis)
- Chi phí đầu tư kỹ thuật thấp: Bảo tồn tại chỗ (In-situ) dựa trên xúc tiến tái sinh tự nhiên tiết kiệm 70% ngân sách so với trồng rừng nhân tạo trên vách đá vôi (vốn có chi phí từ 35 - 50 triệu VNĐ/ha và tỷ lệ sống dưới 45%).
- Lợi ích sinh thái bền vững: Duy trì nguồn sinh thủy cho lưu vực hồ Ba Bể, hạn chế xói mòn đất rửa trôi trên địa hình dốc $25^\circ - 45^\circ$, bảo tồn nguồn gen quý hiếm quốc gia.
Hạn chế và hướng phát triển
Hạn chế kỹ thuật
- Giới hạn thời gian khảo sát: Nghiên cứu được thực hiện trong 4 tháng (mùa khô từ tháng 8 đến tháng 11/2014), chưa theo dõi được trọn vẹn chu kỳ ra hoa, kết quả và phát tán hạt giống trong cả năm.
- Dung lượng mẫu tầng cây gỗ: Quy mô 7 OTC ($3.500\text{ m}^2$) phản ánh tốt các sườn đặc trưng nhưng chưa bao quát được toàn bộ địa hình hiểm trở trên đỉnh Tam Sao (1.159 m).
Hướng phát triển đề xuất
- Xây dựng hệ thống Ô định vị vĩnh viễn (Permanent Sample Plots - PSP): Thiết lập mạng lưới quan sát sinh thái dài hạn để đo lường gia tăng thể tích và tốc độ luân chuyển cá thể (Turnover rate).
- Ứng dụng LiDAR và Viễn thám không gian: Tích hợp mô hình quét tán rừng bằng thiết bị bay không người lái (UAV-LiDAR) để trích xuất chỉ số diện tích lá (LAI) và sinh khối tầng cây gỗ trên toàn bộ 199,0 ha rừng núi đá.
Đối tượng hưởng lợi
+-------------------+-----------------------------------------------------------+
| Nhóm đối tượng | Giá trị và Lợi ích định lượng mang lại |
+-------------------+-----------------------------------------------------------+
| Sinh viên Lâm học | Tài liệu chuẩn hóa phương pháp điều tra lâm học, lập mã |
| & Quản lý TNR | biểu tính toán IV% và đo đếm chỉ số đa dạng sinh học. |
+-------------------+-----------------------------------------------------------+
| Cán bộ Ban Quản lý| Bộ dữ liệu định lượng 13 loài cây gỗ quý hiếm làm căn cứ |
| Khu Bảo tồn | phân bổ lực lượng tuần tra và xây dựng hồ sơ bảo tồn gen. |
+-------------------+-----------------------------------------------------------+
| Cộng đồng địa | Hưởng lợi từ chính sách giao khoán khoanh nuôi bảo vệ rừng|
| phương vùng đệm | và nguồn thu từ dịch vụ môi trường rừng (DVMTR). |
+-------------------+-----------------------------------------------------------+
| Nhà nghiên cứu | Dữ liệu sinh thái cơ sở về thảm thực vật karst Đông Bắc |
| Sinh thái học | phục vụ các nghiên cứu so sánh diễn thế và biến đổi khí hậu|
+-------------------+-----------------------------------------------------------+
Câu hỏi thường gặp
1. Tại sao chỉ số IV% lại là thước đo chuẩn xác hơn số lượng cây thuần túy ($N%$)?
Trong cấu trúc rừng tự nhiên nhiệt đới, một loài có thể có số lượng cây mạ hoặc cây nhỏ rất lớn nhưng sinh khối không đáng kể. Chỉ số $IV% = (N_i% + G_i%)/2$ dung hòa giữa số lượng cá thể ($N_i%$) và diện tích chiếm lĩnh không gian dinh dưỡng thực tế thông qua tiết diện ngang thân cây ($G_i%$), giúp xác định chính xác vị trí ưu thế sinh thái của loài trong lâm phần.
2. Sự khác biệt cơ bản giữa cấu trúc rừng núi đá >700 m và rừng núi đất thấp là gì?
Rừng trên núi đá vôi >700 m có tầng đất mỏng, độ tàn che gián đoạn (0,2 - 0,5), rễ cây len lỏi trong các khe nứt karst. Thành phần loài chịu hạn sinh lý tốt, giàu các loài cây lá kim và gỗ cứng như Nghiến (Excentrodendron tonkinense), Trai lý (Garcinia fagraeoides), Re hương (Cinnamomum parthenoxylon), khác biệt hoàn toàn với ưu thế họ Dầu (Dipterocarpaceae) hay họ Dẻ (Fagaceae) trên núi đất ẩm.
3. Phương pháp đánh giá chất lượng cây tái sinh được thực hiện như thế nào?
Cây tái sinh được phân thành 3 cấp: Cây Tốt (thân thẳng, ngọn vút, không sâu bệnh, sinh trưởng mạnh), Cây Xấu (cong queo, cụt ngọn, cụt tán, sâu bệnh hại), và Cây Trung bình (chất lượng trung gian). Tỷ lệ cây Tốt và Trung bình đạt 62,34% là cơ sở khẳng định tiềm năng phục hồi không cần trồng dặm tốn kém.
4. Quy trình xử lý khi phát hiện loài cực kỳ nguy cấp như Re hương tại thực địa?
Tiến hành ghi nhận tọa độ GPS chính xác, gắn mã số thẻ kim loại theo dõi cá thể, đo đếm chu vi tán lá, chiều cao dưới cành, và thiết lập các ô dạng bản xung quanh gốc cây mẹ để đánh giá khả năng nảy mầm của hạt và cây mạ dưới tán.
5. Làm thế nào để mở rộng mô hình điều tra này cho các khu bảo tồn khác?
Chuẩn hóa mẫu phiếu điều tra Form 01 (Cây gỗ) và Form 02 (Cây tái sinh), áp dụng kịch bản xử lý dữ liệu tự động bằng Python/R và tuân thủ quy chuẩn thiết lập ô tiêu chuẩn $500\text{ m}^2$ trên địa hình karst dốc $>25^\circ$.
Kết luận
Khóa luận tốt nghiệp của tác giả Giàng Quáng Tiên đã giải quyết trọn vẹn bài toán đánh giá định lượng tính đa dạng thực vật thân gỗ tại đai cao trên 700 m thuộc Khu bảo tồn Loài và Sinh cảnh Nam Xuân Lạc. Các đóng góp cốt lõi bao gồm:
- Thiết lập hệ thống 7 ô tiêu chuẩn chuẩn hóa, xác định rõ cấu trúc tổ thành và hệ chỉ số đa dạng sinh thái ($H'$ đạt tới 2,75; Simpson $D_1$ đạt tới 0,98).
- Ghi nhận 13 loài cây gỗ nguy cấp, đặc biệt định vị và đánh giá sinh thái vi quần thể Re hương (Cinnamomum parthenoxylon - Cực kỳ nguy cấp CR) và Nghiến (Excentrodendron tonkinense).
- Chứng minh khả năng tái sinh tự nhiên khả quan với 62,34% cây tái sinh đạt phẩm chất tốt và trung bình.
Công trình là nguồn tài liệu khoa học lâm nghiệp giá trị, cung cấp luận cứ vững chắc để Ban Quản lý KBT Nam Xuân Lạc và Chi cục Kiểm lâm Bắc Kạn ban hành các chính sách bảo tồn nguồn gen, phân bổ ngân sách tuần tra và triển khai các giải pháp lâm sinh phục hồi rừng bền vững.
