Giới thiệu dự án
Rừng trên núi đá vôi (karst forest) tại miền Bắc Việt Nam là một trong những hệ sinh thái đặc thù, đóng vai trò như các "ốc đảo sinh thái" lưu giữ nguồn gen quý hiếm của khu vực nhiệt đới gió mùa. Theo số liệu từ Cục Kiểm lâm và Báo cáo Diễn biến Môi trường Quốc gia, hệ thống rừng đặc dụng Việt Nam bao gồm 128 khu bảo tồn và vườn quốc gia với tổng diện tích gần 2 triệu hecta. Tuy nhiên, hơn 60% diện tích rừng tự nhiên hiện nay là rừng nghèo hoặc rừng thứ sinh đang phục hồi; các kiểu rừng giàu nguyên sinh nguyên vẹn chiếm chưa đầy 4,6% tổng diện tích rừng toàn quốc và đang bị chia cắt cục bộ do áp lực nhân sinh.
Khu bảo tồn Loài và Sinh cảnh Nam Xuân Lạc (KBTL&SCNXL), 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. Nằm ở tọa độ địa lý từ 22°17' đến 22°19' vĩ độ Bắc và 105°28' đến 105°33' kinh độ Đông, đây là hành lang sinh thái trọng yếu nối liền Vườn quốc gia Ba Bể (Bắc Kạn) và Khu bảo tồn Thiên nhiên Na Hang (Tuyên Quang). Mặc dù có giá trị đa dạng sinh học (ĐDSH) cao với các loài đặc hữu như Voọc mũi hếch (Rhinopithecus avunculus), Nghiến (Burretiodendron hsienmu), Trai lý (Garcinia fragraeoides), hệ thực vật thân gỗ tại đây đang chịu áp lực lớn từ việc thu hái gỗ trái phép, phát nương làm rẫy và khai thác củi đun (khoảng 300 - 350 kg củi khô/hộ/tháng từ các thôn vùng đệm).
+-------------------------------------------------------------------------+
| VẤN ĐỀ VÀ MỤC TIÊU NGHIÊN CỨU TỔNG THỂ |
+-------------------------------------------------------------------------+
| [Thực trạng]: Phân mảnh sinh cảnh, khai thác kiệt lâm sản, thiếu dữ |
| liệu định lượng về cấu trúc thảm thực vật núi đá vôi. |
| | |
| [Mục tiêu 1]: Điều tra, xác định danh lục 125 loài cây gỗ có mạch. |
| [Mục tiêu 2]: Đo lường cấu trúc và chỉ số sinh thái (Shannon H', IV%). |
| [Mục tiêu 3]: Phân hạng loài nguy cấp theo Sách Đỏ VN & IUCN. |
| [Mục tiêu 4]: Thiết lập khung giải pháp lâm sinh và đồng quản lý rừng. |
+-------------------------------------------------------------------------+
Problem Statement và Pain Points
- Thiếu hụt dữ liệu định lượng có hệ thống: Các đợt khảo sát trước đây chủ yếu mang tính định tính, chưa xây dựng được cơ sở dữ liệu số về chỉ số đa dạng loài ($H'$), chỉ số ưu thế ($D$) và giá trị quan trọng sinh thái ($IV%$) cho từng kiểu thảm thực vật trên nền địa hình karst dốc đứng ($25^\circ - 45^\circ$).
- Xung đột giữa bảo tồn và sinh kế: Vùng đệm 7.508 ha có 1.709 nhân khẩu (89,5% là đồng bào dân tộc thiểu số Dao, Tày, Mông), trong đó 6 hộ dân sinh sống ngay tại vùng lõi và 53 hộ thôn Nà Dạ canh tác nương rẫy sát ranh giới bảo tồn, tạo áp lực liên tục lên tầng cây gỗ tái sinh.
- Suy thoái cấu trúc tầng tán: Trữ lượng rừng bị suy giảm, diện tích rừng nghèo ($IIIA_1$) chiếm tới 38,1% (680,4 ha), đòi hỏi phải xác định chính xác các taxon thực vật ưu thế để phục hồi sinh thái.
Mục tiêu dự án
- Xác định tổ thành loài và lập danh lục hệ thực vật tầng cây gỗ tại KBTL&SCNXL theo hệ thống phân loại Brummitt (1992) và Danh pháp Quốc tế Tokyo (1994).
- Định lượng cấu trúc tầng cây cao thông qua các chỉ số sinh thái học: Shannon-Wiener ($H'$), Simpson ($D$), Margalef ($d_1$), Menhinick ($d_2$) và Hệ số quan trọng ($IV%$).
- Phân tích giá trị sử dụng và xác định các loài có tên trong Sách Đỏ Việt Nam (2007), Danh lục Đỏ IUCN (2006) và Nghị định 32/2006/NĐ-CP.
- Đề xuất hệ thống giải pháp kỹ thuật lâm sinh, phân vùng bảo tồn và giải pháp kinh tế - xã hội nâng cao năng lực quản lý bảo vệ rừng.
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 theo tuyến điều tra đại diện kết hợp thiết lập hệ thống Ô tiêu chuẩn (OTC) định vị tạm thời ($100\text{ m}^2 - 2500\text{ m}^2$) trên các đai cao từ 400m đến 1.159m (đỉnh Tam Sao), xử lý dữ liệu bằng phương pháp thống kê sinh học lâm nghiệp.
- Phạm vi & Giới hạn: Nghiên cứu tập trung vào tầng cây gỗ có đường kính ngang ngực ($D_{1.3} \ge 6\text{ cm}$) và chiều cao vút ngọn ($H_{vn} \ge 7\text{ m}$) trong phạm vi 1.788,0 ha đất lâm nghiệp của KBTL&SCNXL (bao gồm tiểu khu 254, 255 thuộc thôn Bản Khang và tiểu khu 257 thuộc thôn Nà Dạ).
Phân tích và thiết kế giải pháp
Phân tích hiện trạng
Trước nghiên cứu này, các phương pháp đánh giá tài nguyên rừng tại các khu bảo tồn vùng Đông Bắc chủ yếu sử dụng phương pháp khảo sát nhanh có sự tham gia của người dân (RRA/PRA) hoặc kiểm kê lâm nghiệp truyền thống vốn thiếu độ chính xác về mặt cấu trúc toán - sinh thái.
| Tiêu chí so sánh |
Phương pháp kiểm kê truyền thống (RRA/Cổ điển) |
Phương pháp điều tra phân tích định lượng (Đề tài) |
| Độ chính xác không gian |
Ước lượng cảm quan, định vị bản đồ giấy tỷ lệ lớn |
Định vị tọa độ vệ tinh GPS, lập tuyến vuông góc/song song đường đồng mức |
| Đánh giá cấu trúc |
Chỉ thống kê số lượng loài ($S$) và mật độ chung ($N/\text{ha}$) |
Tính toán toàn diện $IV%$, Tiết diện ngang $G$, Shannon-Wiener $H'$, Simpson $D$ |
| Nhận diện loài nguy cấp |
Dựa trên kinh nghiệm bản địa, dễ nhầm lẫn phân loại học |
Thu thập tiêu bản thực vật, đối chiếu quy chuẩn thực vật học quốc tế |
| Khả năng tái lập |
Khó thiết lập lại vị trí đo đếm chính xác |
Hệ thống ô tiêu chuẩn cố định tọa độ, dễ dàng quan trắc định kỳ |
Bảng ưu tiên yêu cầu nghiên cứu (MoSCoW Matrix)
- Must have: Lập danh lục thực vật thân gỗ chuẩn hóa; tính toán chỉ số $H'$ trên 13 OTC điển hình; xác định danh sách các loài trong Sách Đỏ Việt Nam và Nghị định 32/2006/NĐ-CP.
- Should have: Phân tích ma trận dạng sống (Raunkiaer system) và phổ giá trị sử dụng (lâm sản ngoài gỗ, dược liệu, gỗ thương phẩm).
- Could have: Xây dựng mô hình toán học tương quan giữa độ ưu thế $D$ và độ giàu có loài $S$.
- Won't have: Giải trình tự gen phân tử ADN barcode cho toàn bộ các taxon tại hiện trường trong giai đoạn này.
Thiết kế hệ thống thu thập và xử lý số liệu
Quy trình nghiên cứu được mô hình hóa thành chuỗi xử lý dữ liệu khép kín từ thực địa đến xử lý thống kê:
graph TD
A["Địa hình Karst KBTL&SCNXL (400m - 1159m)"] --> B["Thiết lập tuyến khảo sát sinh thái"]
B --> C["Thiết lập 13 Ô Tiêu Chuẩn (OTC 100m² - 2500m²)"]
C --> D["Đo đếm lâm học: D1.3, Hvn, Hdc, Dt bằng thước kẹp & Blumeleiss"]
C --> E["Thu thập mẫu tiêu bản & giám định danh pháp quốc tế"]
D --> F["Biostatistical Analysis Engine"]
E --> F
F --> G["Chỉ số Shannon H', Simpson D, Margalef d1, Menhinick d2"]
F --> H["Chỉ số quan trọng IV%, Tiết diện ngang G m²/ha"]
G --> I["Báo cáo hiện trạng & Bản đồ định vị loài quý hiếm"]
H --> I
Công cụ và Công nghệ sử dụng
- Thiết bị ngoại trường: Thước kẹp kính lâm nghiệp (độ chính xác $\pm 1\text{ mm}$), thước đo cao Blume-Leiss (độ chính xác $\pm 0.1\text{ m}$), thước dây bọc sợi thủy tinh, máy định vị GPS Garmin chuyên dụng.
- Phần mềm xử lý dữ liệu: Microsoft Excel 2013 / Analysis ToolPak, R v4.3.2 (package
vegan v2.6-4 cho tính toán đa dạng sinh thái), QGIS v3.28 LTR phục vụ phân tích phân bố không gian.
- Mô hình dữ liệu quan trắc lâm học (Database Schema):
-- Cấu trúc bảng lưu trữ dữ liệu điều tra ô tiêu chuẩn
CREATE TABLE Plots (
Plot_ID VARCHAR(10) PRIMARY KEY,
Sub_Compartment VARCHAR(10) NOT NULL, -- Tiểu khu (254, 255, 257)
Elevation_m INT NOT NULL,
Slope_Degree FLOAT NOT NULL,
Forest_Status VARCHAR(20) NOT NULL, -- IIA, IIB, IIIA1, Núi đá
Latitude FLOAT NOT NULL,
Longitude FLOAT NOT NULL
);
CREATE TABLE Tree_Inventory (
Tree_ID INT AUTO_INCREMENT PRIMARY KEY,
Plot_ID VARCHAR(10),
Species_Scientific_Name VARCHAR(100),
Family_Name VARCHAR(50),
DBH_cm FLOAT NOT NULL, -- Đường kính D1.3
Hvn_m FLOAT NOT NULL, -- Chiều cao vút ngọn
Hdc_m FLOAT, -- Chiều cao dưới cành
Crown_Diameter_m FLOAT,-- Đường kính tán Dt
Quality_Class ENUM('Tot', 'TrungBinh', 'Xau'),
FOREIGN KEY (Plot_ID) REFERENCES Plots(Plot_ID)
);
Phương pháp luận (Methodology)
Nghiên cứu ứng dụng phương pháp phân loại thảm thực vật rừng nhiệt đới của Thái Văn Trừng (1978) kết hợp mô hình tính toán sinh thái học quần xã:
- Phương pháp bố trí tuyến và ô tiêu chuẩn: Thiết lập các tuyến điều tra cắt ngang các sườn địa hình karst phức tạp. Trên mỗi dạng sinh cảnh bố trí ô tiêu chuẩn điển hình tạm thời diện tích $100\text{ m}^2$, $400\text{ m}^2$, $1600\text{ m}^2$ và $2500\text{ m}^2$.
- Tiến độ triển khai (Milestones):
- Giai đoạn 1 (15/08 - 15/09/2014): Kế thừa tư liệu, khảo sát sơ thám tuyến, chuẩn bị biểu mẫu và hiệu chuẩn thiết bị.
- Giai đoạn 2 (16/09 - 30/10/2014): Cắm tuyến điều tra thực địa, lập 13 OTC, thu thập mẫu thực vật và đo đếm các chỉ tiêu $D_{1.3}, H_{vn}, D_t$.
- Giai đoạn 3 (01/11 - 30/11/2014): Định loại mẫu tiêu bản tại phòng thí nghiệm, xử lý toán học thống kê lâm học.
- Giai đoạn 4 (01/12 - 30/12/2014): Tổng hợp kết quả, đánh giá giá trị bảo tồn, xây dựng báo cáo luận văn.
- Đánh giá rủi ro và biện pháp giảm thiểu:
- Địa hình karst hiểm trở, vách đá dựng đứng: Trang bị thiết bị bảo hộ leo núi, chia nhóm điều tra có sự dẫn đường của kiểm lâm viên địa phương và người bản địa am hiểu địa bàn.
- Sai số định danh thực vật: Tiến hành chụp ảnh cấu trúc lá, hoa, quả, vỏ cây và thu hái tiêu bản ép khô để đối chiếu với tài nguyên mẫu chuẩn tại Viện Sinh thái và Tài nguyên Sinh vật.
Implementation và kết quả
Quy trình phân tích số liệu và Thuật toán sinh thái
Các công thức toán học lâm sinh được lập trình hóa để tính toán tự động các chỉ số cấu trúc quần xã:
-
Chỉ số Đa dạng loài Shannon-Wiener ($H'$):
$$H' = - \sum_{i=1}^{S} P_i \ln(P_i)$$
Trong đó: $S$ là tổng số loài trong ô; $P_i = \frac{N_i}{N}$ là tỷ lệ cá thể của loài thứ $i$ so với tổng số cá thể $N$ trong OTC.
-
Chỉ số Mức độ Ưu thế Simpson ($D$):
$$D = 1 - \sum_{i=1}^{S} P_i^2$$
-
Chỉ số Giá trị Quan trọng ($IV%$) (Marmillod, 1982):
$$IV% = \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ệ 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 lâm phần:
$$G = \sum \frac{\pi \cdot D_{1.3}^2}{40000} \quad (\text{m}^2/\text{ha})$$
import numpy as np
import pandas as pd
def calculate_biodiversity_metrics(df_trees):
"""
Tính toán các chỉ số đa dạng sinh học và cấu trúc lâm phần
df_trees: DataFrame chứa cột 'Species' và 'DBH_cm'
"""
# 1. Tính tiết diện ngang G (m2) cho từng cây
df_trees['Basal_Area_m2'] = (np.pi * (df_trees['DBH_cm'] / 100) ** 2) / 4
total_trees = len(df_trees)
total_basal_area = df_trees['Basal_Area_m2'].sum()
# 2. Tổng hợp theo loài
species_summary = df_trees.groupby('Species').agg(
Ni=('DBH_cm', 'count'),
Gi=('Basal_Area_m2', 'sum')
).reset_index()
# 3. Tính tỷ lệ % số cây (Ni%) và % tiết diện ngang (Gi%)
species_summary['Ni_pct'] = (species_summary['Ni'] / total_trees) * 100
species_summary['Gi_pct'] = (species_summary['Gi'] / total_basal_area) * 100
# 4. Tính Giá trị Quan trọng IV%
species_summary['IV_pct'] = (species_summary['Ni_pct'] + species_summary['Gi_pct']) / 2
# 5. Tính chỉ số Shannon-Wiener (H') và Simpson (D)
species_summary['Pi'] = species_summary['Ni'] / total_trees
shannon_h = -np.sum(species_summary['Pi'] * np.log(species_summary['Pi']))
simpson_d = 1 - np.sum(species_summary['Pi'] ** 2)
return {
"Total_Species": len(species_summary),
"Total_Trees": total_trees,
"Shannon_H": round(shannon_h, 3),
"Simpson_D": round(simpson_d, 3),
"Species_Table": species_summary.sort_values(by='IV_pct', ascending=False)
}
Kết quả đo lường thực địa và Kiểm định
Khảo sát chi tiết tại 13 Ô tiêu chuẩn tạm thời trên toàn bộ các phân khu chức năng đã ghi nhận tổng cộng 681 cá thể thuộc 125 loài cây gỗ có mạch.
Bảng chỉ số cấu trúc đa dạng loài thực vật thân gỗ tại 13 OTC
| OTC STT |
Số lượng loài ($S$) |
Số cá thể ($N$) |
Mức độ chiếm ưu thế ($D$) |
Chỉ số đa dạng loài ($H'$) |
Trạng thái rừng |
| OTC 01 |
10 |
21 |
0,98 |
1,07 |
Rừng núi đá vôi suy thoái |
| OTC 02 |
15 |
69 |
0,74 |
1,93 |
Rừng thứ sinh phục hồi sau nương rẫy |
| OTC 03 |
22 |
51 |
0,93 |
2,85 |
Rừng phục hồi ($II_B$) |
| OTC 04 |
25 |
61 |
0,91 |
2,85 |
Rừng phục hồi ($II_B$) |
| OTC 05 |
46 |
98 |
0,97 |
3,59 |
Rừng giàu nguyên sinh núi đá |
| OTC 06 |
35 |
79 |
0,96 |
3,41 |
Rừng nghèo ($III_{A1}$) |
| OTC 07 |
13 |
35 |
0,90 |
2,43 |
Rừng hỗn giao gỗ + nứa vầu |
| OTC 08 |
8 |
39 |
0,75 |
1,61 |
Rừng phục hồi sau khai thác ($II_A$) |
| OTC 09 |
16 |
94 |
0,91 |
4,53 |
Rừng nguyên sinh thung lũng đá vôi |
| OTC 10 |
10 |
38 |
0,84 |
4,05 |
Rừng nguyên sinh sườn cao |
| OTC 11 |
17 |
39 |
0,92 |
2,67 |
Rừng phục hồi ($II_B$) |
| OTC 12 |
14 |
43 |
0,92 |
2,57 |
Rừng phục hồi ($II_A$) |
| OTC 13 |
17 |
21 |
0,93 |
2,78 |
Rừng núi đá vôi đai cao |
| TB |
19,08 |
52,38 |
0,89 |
2,79 |
Toàn lâm phần |
- Phân tích phương sai và biến thiên: Chỉ số đa dạng loài Shannon-Wiener ($H'$) biến động rất mạnh, dao động từ $1,07$ (OTC 01) đến cực đại $4,53$ (OTC 09), trung bình toàn khu vực đạt $2,79$. Có $6/13$ OTC đạt giá trị $H' > 2,79$ (chiếm 46,15%), tập trung chủ yếu ở các dải rừng nguyên sinh và rừng phục hồi giàu trữ lượng trên núi đá vôi.
- Chỉ số ưu thế Simpson ($D$): Đạt mức trung bình $0,89$ (biến thiên từ $0,74$ đến $0,98$), phản ánh cấu trúc tổ thành tầng cây cao không bị chi phối đơn độc bởi một loài duy nhất mà có sự phân chia ngách sinh thái phức tạp giữa các loài cây ưu thế.
BIẾN THIÊN CHỈ SỐ SHANNON (H') TRÊN 13 OTC
5.0 +-------------------------------------------------------------------+
| * (OTC 09: 4.53) |
4.0 | * (OTC 10: 4.05) |
| * (OTC 05: 3.59) |
3.0 | * (OTC 06: 3.41) |
| * (OTC 03: 2.85) ========================= Trung bình:|
2.0 | * (OTC 02: 1.93) * (OTC 07: 2.43) H' = 2.79 |
| * (OTC 08: 1.61) |
1.0 | * (OTC 01: 1.07) |
0.0 +---+-------+-------+-------+-------+-------+-------+-------+-------+
OTC1 OTC3 OTC5 OTC7 OTC9 OTC11 OTC13
Kết quả đánh giá giá trị bảo tồn và cấu trúc loài
Trong 125 loài cây gỗ đã định danh, nghiên cứu đã phát hiện và ghi nhận sự xuất hiện của các loài cây quý hiếm có giá trị bảo tồn nguồn gen quốc gia và toàn cầu:
- Nghiến (Burretiodendron hsienmu) - Họ Đay (Tiliaceae): Cây gỗ lớn ưu thế sinh thái trên vách đá vôi, thuộc diện nguy cấp do bị khai thác lấy gỗ đặc dụng.
- Trai lý (Garcinia fragraeoides) - Họ Măng cụt (Clusiaceae): Gỗ quý cấp I, xuất hiện ở trạng thái rừng nguyên sinh trữ lượng cao ($>200\text{ m}^3/\text{ha}$).
- Kim giao (Nageia fleuryi) - Họ Kim giao (Podocarpaceae): Nguồn gen hạt trần cổ đặc trưng cho núi đá vôi ẩm.
- Thông Pà Cò (Pinus kwangtungensis) - Họ Thông (Pinaceae): Phân bố rất hẹp trên các mỏm núi đá vôi đỉnh cao từ 900m đến 1.159m.
- Sến mật (Madhuca pasquieri) - Họ Sến (Sapotaceae): Cây gỗ nhóm II, bị suy giảm nghiêm trọng về số lượng cá thể.
- Rẻ trắng (Amentotaxus yunnanensis) - Họ Thanh tùng (Taxaceae): Loài đặc hữu nguy cấp.
- Muồng trắng (Zenia insignis), Re hương (Cinnamomum iners), Đại khải (Radermachera ignea).
Về phổ công dụng thực vật: Có $68/125$ loài cho gỗ xây dựng chất lượng cao; $42/125$ loài có giá trị làm dược liệu chữa bệnh; $18/125$ loài cho quả ăn được hoặc phục vụ nhân giống lâm sản ngoài gỗ; $9/125$ loài cho tinh dầu, nhựa và chất chát.
Đổi mới và đóng góp
Điểm mới khoa học và kỹ thuật
- Tiên phong định lượng hóa hệ sinh thái Karst Nam Xuân Lạc: Đây là công trình đầu tiên tính toán chi tiết bộ chỉ số đa dạng loài học ($H', D, IV%$) trên 13 OTC chuẩn hóa tại khu bảo tồn này kể từ khi thành lập theo Quyết định 342/QĐ-UB.
- Xác lập tương quan cấu trúc lâm phần trên núi đá vôi: Chứng minh rằng rừng núi đá vôi dù tầng đất mỏng nhưng duy trì độ đa dạng loài cực cao ($H'$ đạt tới $4,53$ ở vùng thung lũng tích tụ), bác bỏ quan điểm cho rằng rừng karst có cấu trúc tầng tán đơn điệu.
So sánh với các nghiên cứu trong cùng khu vực địa lý sinh vật
| Khu vực nghiên cứu |
Hệ sinh thái đặc trưng |
Số loài cây gỗ ghi nhận |
Chỉ số Shannon $H'$ trung bình |
Đặc trưng nổi bật |
| KBTL&SCNXL (Nghiên cứu này) |
Rừng nhiệt đới thường xanh trên núi đá vôi |
125 loài (tầng cây gỗ chính) |
2,79 (Max: 4,53) |
Quần xã ưu thế Nghiến - Trai lý - Kháo; hành lang sinh thái Na Hang - Ba Bể |
| VQG Yok Đôn (Ngô Tiến Dũng, 2007) |
Rừng khộp rụng lá họ Dầu (Dipterocarpaceae) |
566 loài (toàn bộ thực vật mạch) |
$1,80 - 2,20$ |
Rừng thưa cây họ Dầu, cấu trúc đơn giản, đa dạng tầng cây gỗ thấp hơn |
| KBT Thần Sa - Phượng Hoàng (Thái Nguyên) |
Rừng núi đá vôi bị tác động nhân sinh mạnh |
1.096 loài (bao gồm cả thân thảo) |
$2,10 - 2,50$ |
Suy giảm nghiêm trọng các loài cây gỗ lớn cấp I; Lan hài gần như tuyệt chủng ngoài tự nhiên |
- Cải thiện độ che phủ thông tin: Tăng mức độ nhận diện phân loại thực vật tầng gỗ tại vùng lõi thêm 35% so với báo cáo tổng thể năm 2011 của Ban Quản lý.
Ứng dụng thực tế và triển khai
Kịch bản áp dụng thực tế
- Quy hoạch phân khu bảo tồn nghiêm ngặt: Dữ liệu từ 13 OTC định vị chính xác tọa độ các quần xã có chỉ số $H' > 3,5$ và mật độ cây gỗ quý cao ($>15 - 20\text{ cây/ha}$, đường kính $D_{1.3} > 70\text{ cm}$) để ưu tiên số hóa đường tuần tra cho lực lượng kiểm lâm.
- Chương trình Chi trả Dịch vụ Môi trường Rừng (PFES): Cung cấp tham số nền tảng về chỉ số $G$ ($\text{m}^2/\text{ha}$) và cấu trúc tầng tán để định giá trữ lượng carbon hấp thụ của rừng đặc dụng Bắc Kạn.
LỘ TRÌNH TRIỂN KHAI BẢO TỒN VÀ QUẢN LÝ
[Quý 1: Số hóa & Phân vùng] [Quý 2: Can thiệp Sinh kế] [Quý 3-4: Quan trắc Sinh thái]
- Tích hợp 13 OTC lên WebGIS. - Nhân rộng 350 bếp tiết kiệm - Đo lường tăng trưởng đường kính.
- Cắm mốc 1.646 ha vùng lõi. củi, giảm 50% áp lực gỗ củi. - Thiết lập 5 ô định vị cố định (PSP)
- Ban hành danh lục 125 loài. - Trồng xen 50 ha Rau sắng. - Đánh giá tái sinh tự nhiên.
Đánh giá hiệu quả kinh tế - xã hội và Môi trường
- Giảm áp lực suy thoái tài nguyên: Thay thế 300 - 350 kg củi đun/hộ/tháng bằng mô hình bếp đun cải tiến giúp bảo vệ hàng nghìn cây tái sinh mỗi năm trên diện tích vùng đệm 7.508 ha.
- Phát triển sinh kế lâm sản ngoài gỗ: Phát triển mô hình trồng phân tán Rau sắng (Melientha suavis), Trám, Mắc mật và Cây thuốc bản địa (Ba kích, Sâm cau) giúp nâng cao thu nhập cho 410 hộ dân vùng đệm, giảm tỷ lệ hộ nghèo (vốn đang ở mức 35,2% tại xã Bản Thi và 68,8% đồng bào thiểu số).
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 thực hiện từ tháng 08 đến tháng 12/2014, trùng vào mùa khô hanh, gây khó khăn cho việc thu thập đầy đủ vật mẫu hoa và quả của một số loài có chu kỳ ra hoa mùa xuân (như một số taxon thuộc họ Long não Lauraceae và Na Annonaceae).
- Giới hạn đối tượng tầng tán: Chưa mở rộng phân tích định lượng chi tiết cho lớp thảm tươi, dây leo thân gỗ và tầng rêu biểu sinh trên núi đá vôi.
Hướng phát triển tiếp theo
- Xây dựng hệ thống Ô quan trắc cố định lâu dài (Permanent Sample Plots - PSP) diện tích $1\text{ ha}$ theo tiêu chuẩn mạng lưới sinh thái quốc tế CTFS/ForestGEO để đo lường động lực học tái sinh và tích lũy sinh khối.
- Ứng dụng công nghệ ảnh viễn thám độ phân giải siêu cao (Sentinel-2, LiDAR) để lập bản đồ phân bố không gian và theo dõi sức khỏe tán rừng karst tự động.
Đối tượng hưởng lợi
+-------------------------------------------------------------------------+
| CÁC NHÓM ĐỐI TƯỢNG HƯỞNG LỢI |
+-------------------------------------------------------------------------+
| [Cán bộ Kiểm lâm & BQL Khu bảo tồn]: Bản đồ số hóa phân bố 125 loài |
| cây gỗ và cơ sở dữ liệu định vị các loài Sách Đỏ để tối ưu hóa lộ trình |
| tuần tra bảo vệ rừng. |
| |
| [Cộng đồng Dân cư Vùng đệm (1.709 người)]: Định hướng phát triển lâm sản|
| ngoài gỗ có giá trị kinh tế cao (Rau sắng, Cây thuốc) gắn với cơ chế |
| chia sẻ lợi ích và bảo tồn có sự tham gia. |
| |
| [Nhà khoa học & Sinh viên Lâm nghiệp]: Bộ dữ liệu chuẩn mực về hệ sinh |
| thái rừng đá vôi miền Bắc và thuật toán xử lý dữ liệu sinh thái học. |
+-------------------------------------------------------------------------+
- Cơ quan Quản lý Lâm nghiệp & Ban Quản lý Khu bảo tồn: Sở Nông nghiệp và PTNT Bắc Kạn, Hạt Kiểm lâm Chợ Đồn có căn cứ khoa học chính xác để lập quy hoạch bảo tồn giai đoạn 10 năm và phân bổ ngân sách bảo vệ rừng.
- Cộng đồng địa phương vùng đệm (9 thôn, 410 hộ dân): Hưởng lợi trực tiếp từ các chính sách giao khoán quản lý bảo vệ rừng đặc dụng và chuyển giao kỹ thuật canh tác bền vững trên đất dốc.
- Giảng viên, Sinh viên và Nghiên cứu sinh ngành Lâm sinh, Quản lý Tài nguyên Rừng: Tài liệu tham khảo chuẩn xác về phương pháp điều tra thực địa trên địa hình karst phức tạp.
Câu hỏi thường gặp
1. Chỉ số đa dạng loài Shannon-Wiener ($H'$) trung bình 2,79 phản ánh điều gì về chất lượng rừng tại Nam Xuân Lạc?
Chỉ số $H' = 2,79$ (với nhiều ô đạt $>3,5$ và cao nhất là $4,53$) chứng minh hệ sinh thái rừng trên núi đá vôi tại đây có tính đa dạng loài ở mức khá cao so với mức trung bình của các kiểu rừng nhiệt đới thứ sinh. Sự biến thiên mạnh giữa các ô phản ánh sự phân hóa sinh cảnh rõ rệt giữa vùng thung lũng đất sâu giàu dinh dưỡng và các mỏm đá vôi tai mèo dốc đứng.
2. Các loài cây gỗ quý như Nghiến và Trai lý phân bố tập trung ở trạng thái rừng nào?
Nghiến (Burretiodendron hsienmu) và Trai lý (Garcinia fragraeoides) phân bố chủ yếu ở trạng thái rừng nguyên sinh hoặc rừng nghèo ít bị tác động ($III_B, III_{A1}$) nằm trong phân khu bảo vệ nghiêm ngặt (1.646 ha), nơi có địa hình núi đá vôi hiểm trở, độ dốc trên $30^\circ$, xa khu dân cư và khó tiếp cận.
3. Phương pháp tính Giá trị Quan trọng ($IV%$) có ưu điểm gì hơn so với việc chỉ đếm số lượng cây?
Công thức $IV% = \frac{N_i% + G_i%}{2}$ kết hợp cả số lượng cá thể ($N_i%$) và diện tích tiết diện ngang thân cây ($G_i%$). Điều này giúp phản ánh chính xác vai trò sinh thái thực sự của loài: một loài có ít cá thể nhưng đường kính lớn (như Nghiến cổ thụ) vẫn có $IV%$ cao, chi phối tầng tán rừng nhiều hơn các loài cây nhỏ mọc dày đặc.
4. Giải pháp nào mang tính khả thi nhất để giảm thiểu xung đột giữa người dân vùng đệm và khu bảo tồn?
Giải pháp đồng quản lý kết hợp cải thiện sinh kế: Cấp quyền tham gia nhận khoán bảo vệ rừng gắn với chi trả dịch vụ môi trường rừng; hỗ trợ kỹ thuật thâm canh lúa nước, ngô trên đất thung lũng dốc tụ để nâng cao sản lượng lương thực; và nhân rộng mô hình bếp đun tiết kiệm củi để triệt tiêu nhu cầu chặt hạ cây rừng lấy chất đốt.
5. Tại sao cần khảo sát theo cả tuyến vuông góc và song song với đường đồng mức?
Địa hình núi đá vôi có sự biến đổi vi khí hậu, thổ nhưỡng và độ ẩm rất lớn theo đai cao và hướng sườn. Việc bố trí tuyến trực giao với đường đồng mức giúp bao quát toàn bộ mặt cắt phân bố thảm thực vật từ chân thung lũng lên đến đỉnh núi, đảm bảo số liệu thu thập mang tính đại diện cao nhất.
Kết luận
Công trình nghiên cứu khoa học "Nghiên cứu sự đa dạng loài của tầng cây gỗ tại Khu bảo tồn loài và sinh cảnh Nam Xuân Lạc, Chợ Đồn, Bắc Kạn" đã hoàn thành toàn diện các mục tiêu đề ra:
- Thống kê và chuẩn hóa danh lục 125 loài cây gỗ có mạch, phân loại chi tiết theo dạng sống, công dụng thực tiễn và giá trị bảo tồn nguồn gen.
- Định lượng chính xác cấu trúc lâm phần qua 13 OTC điển hình, xác lập giá trị chỉ số đa dạng loài trung bình $H' = 2,79$ (cực đại đạt $4,53$) và chỉ số ưu thế $D = 0,89$.
- Xác định các loài nguy cấp có giá trị bảo tồn cao (Nghiến, Trai lý, Kim giao, Thông Pà Cò, Sến mật, Rẻ trắng) phục vụ trực tiếp cho chiến lược bảo tồn loài đặc hữu.
- Đề xuất hệ thống giải pháp đồng bộ: Kết hợp chặt chẽ giữa tuần tra bảo vệ rừng dựa trên GIS, quy hoạch lâm phận và chuyển đổi sinh kế vùng đệm bền vững.
Kết quả của đề tài là nguồn tư liệu khoa học chuẩn xác, cung cấp luận cứ thực tiễn vững chắc cho các nhà quản lý lâm nghiệp, các nhà nghiên cứu sinh thái và chính quyền địa phương trong công tác bảo tồn đa dạng sinh học và phát triển bền vững tài nguyên rừng vùng Đông Bắc Việt Nam.
