Giới thiệu dự án

Hệ sinh thái rừng trên núi đá vôi (karst ecosystem) tại khu vực Đông Nam Á được xem là một trong những trung tâm đa dạng sinh học (ĐDSH) quan trọng nhất toàn cầu nhưng cũng dễ bị tổn thương nhất trước áp lực nhân sinh và biến đổi khí hậu. Tại Việt Nam, với tổng diện tích rừng tự nhiên khoảng 12,3 triệu ha (độ che phủ 37%), các dải rừng nhiệt đới gió mùa trên nền địa chất đá vôi phân bố rải rác từ Tây Bắc đến Đông Bắc đang đối mặt với tốc độ suy thoái cục bộ nghiêm trọng. Mặc dù diện tích che phủ chung tăng lên nhờ rừng trồng và phục hồi tự nhiên, chất lượng rừng nguyên sinh chứa đựng nguồn gen quý hiếm suy giảm rõ rệt.

Khu Bảo tồn 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/3/2004 của UBND tỉnh Bắc Kạn với tổng diện tích 1.788,0 ha (vùng lõi bảo vệ nghiêm ngặt 1.646,0 ha và phục hồi sinh thái 142,0 ha), tiếp giáp hành lang sinh thái Na Hang (Tuyên Quang) và Vườn quốc gia Ba Bể. Tuy nhiên, áp lực dân sinh từ vùng đệm diện tích 7.508 ha với 1.709 nhân khẩu (tỷ lệ hộ nghèo tại xã Xuân Lạc lên tới 30,4%, tiêu thụ 300–350 kg củi khô/hộ/tháng) cùng tập quán canh tác nương rẫy, khai thác chọn lọc gỗ quý (Excentrodendron tonkinense, Garcinia fagraeoides) và săn bắt động vật hoang dã đã làm biến dạng cấu trúc tầng tán tự nhiên.

+-------------------------------------------------------------------------+
|                  KHU BẢO TỒN LOÀI VÀ SINH CẢNH NAM XUÂN LẠC             |
|                             (Tổng: 1.788 ha)                            |
+------------------------------------+------------------------------------+
|  Phân khu Bảo vệ nghiêm ngặt:     |  Phân khu Phục hồi sinh thái:      |
|  1.646,0 ha (92,06%)               |  142,0 ha (7,94%)                  |
|  - Rừng nghèo (IIIA1): 680,4 ha    |  - 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 nứa/gỗ: 317,3 ha  +------------------------------------+
|  - Rừng núi đá: 199,0 ha           |  Vùng đệm: 7.508 ha (9 thôn)       |
+------------------------------------+------------------------------------+

Mục tiêu dự án

  1. Xác định đặc điểm cấu trúc tầng tán và thành phần tổ thành loài thực vật thân gỗ thuộc kiểu rừng kín thường xanh trên núi đá vôi ở hai đai cao (< 500m và 500–700m).
  2. Định lượng các chỉ số đa dạng sinh học ($H'$, $D_1$, $d_1$, $d_2$, $IV%$) của tầng cây cao và tầng cây tái sinh.
  3. Đánh giá hiện trạng phân bố các loài cây gỗ quý hiếm có giá trị bảo tồn cao trong Sách Đỏ Việt Nam và Danh lục Đỏ IUCN.
  4. Xây dựng gói giải pháp lâm sinh và quản lý tài nguyên rừng bền vững dựa vào cộng đồng địa phương.

Phạm vi và giới hạn nghiên cứu

  • Không gian: Đai cao dưới 700m thuộc KBT Nam Xuân Lạc (tọa độ chân núi Tam Sao, Lũng Lỳ và Kéo Nàng).
  • Đối tượng: Quần xã thực vật thân gỗ tầng cây cao ($D_{1.3} \ge 6\text{ cm}$, $H_{vn} > 7\text{ m}$) và tầng cây tái sinh ($H_{ts} \ge 0,5\text{ m}$).
  • Thời gian điều tra thực địa: Từ 18/08/2014 đến 31/11/2014.

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

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

Trước đây, các nghiên cứu kiểm kê tài nguyên rừng chủ yếu dựa trên phương pháp thống kê mô tả định tính hoặc viễn thám độ phân giải thấp, không phản ánh được mức độ suy thoái ngầm dưới tán rừng karst chia cắt mạnh (độ dốc 25°–45°).

Phương pháp tiếp cận Ưu điểm Nhược điểm Đánh giá tính khả thi tại Nam Xuân Lạc
Quan trắc định tính truyền thống Triển khai nhanh, chi phí thấp Không định lượng được chỉ số đa dạng, sai số lớn do chủ quan Không đáp ứng yêu cầu phân cấp bảo tồn
Viễn thám quang học đơn thuần Bao quát diện rộng, định vị nhanh Sai số tầng tán cao do tán đứt quãng (độ tàn che 0,3–0,7) Chỉ phân loại được trạng thái thảm phủ
Điều tra ô tiêu chuẩn kết hợp định lượng sinh thái Độ chính xác cao, định lượng $IV%$, cấu trúc tái sinh Đòi hỏi nhân lực chuyên môn, tốn thời gian thực địa Lựa chọn tối ưu (Được triển khai trong đề tài)

Thiết kế hệ thống thu thập và phân tích dữ liệu

[Khảo sát 04 Tuyến thực địa] ---> [Thiết lập 09 OTC (500m2/ô)] ---> [Thu mẫu & Định danh]
                                              |                         (Brummitt 1992)
                                              v                               |
                                  [45 Ô Dạng bản ODB (25m2)]                  v
                                              |                    [Số hóa MS Excel / R Engine]
                                              v                               |
                                  [Đo đếm D1.3, Hvn, Hdc, Dt]                 v
                                              |                    [Tính toán IV%, H', D1, d1, Ki]
                                              +------------------------------>|
                                                                              v
                                                                 [Báo cáo Đa dạng & Bảo tồn]

Yêu cầu kỹ thuật và công nghệ áp dụng

  • Thiết bị định vị: GPS Garmin GPSMAP 64s (độ chính xác sai số < 3m trong điều kiện tán rừng karst).
  • Thiết bị đo đắc lâm sinh: Thước Blumeleiss v1.0 (đo $H_{vn}$, $H_{dc}$ độ chính xác dm), thước kẹp kính 0–100 cm (đo $D_{1.3}$ theo 2 hướng Đông - Tây, Nam - Bắc), thước dây mềm (đo đường kính tán $D_t$).
  • Chuẩn phân loại sinh học: Hệ thống phân loại thực vật Brummitt (1992) và Luật danh pháp quốc tế về thực vật (ICBN Tokyo Code 1994).
  • Môi trường xử lý dữ liệu: Microsoft Excel 2016 / Python v3.10 (thư viện scipy, pandas) phục vụ mô hình hóa chỉ số đa dạng.

Implementation và kết quả

Quy trình tính toán chỉ số sinh thái

Chỉ số giá trị quan trọng ($IV%$) của loài thứ $i$ được xác định theo mô hình Daniel Marmillod: $$IV_i% = \frac{N_i% + G_i%}{2}$$ Trong đó:

  • $N_i%$: Tỷ lệ phần trăm số cá thể của loài $i$ so với tổng số cá thể trong lâm phần ($N_i% = \frac{N_i}{N} \times 100$).
  • $G_i%$: Tỷ lệ phần trăm tiết diện ngang thân cây loài $i$ so với tổng tiết diện ngang $\sum G$ ($G_i = \sum \frac{\pi D_{1.3}^2}{4 \times 10000}$).

Hệ thống chỉ số đa dạng sinh thái áp dụng:

  • 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}}$
  • Hệ số tổ thành cây tái sinh ($K_i$): $K_i = \frac{N_i}{N} \times 10$
import numpy as np
import pandas as pd

def calculate_biodiversity_metrics(species_counts: dict, basal_areas: dict):
    total_n = sum(species_counts.values())
    total_g = sum(basal_areas.values())
    s = len(species_counts)
    
    metrics = []
    pi_sq_sum = 0.0
    h_prime = 0.0
    
    for sp, n_i in species_counts.items():
        g_i = basal_areas[sp]
        ni_pct = (n_i / total_n) * 100
        gi_pct = (g_i / total_g) * 100
        iv_pct = (ni_pct + gi_pct) / 2.0
        
        pi = n_i / total_n
        pi_sq_sum += pi ** 2
        h_prime -= pi * np.log(pi) if pi > 0 else 0
        
        metrics.append({"Species": sp, "Ni": n_i, "IV_pct": round(iv_pct, 2)})
        
    simpson_d1 = 1.0 - pi_sq_sum
    margalef_d1 = (s - 1) / np.log(total_n)
    menhinick_d2 = s / np.sqrt(total_n)
    
    return {
        "Shannon_H": round(h_prime, 3),
        "Simpson_D1": round(simpson_d1, 3),
        "Margalef_d1": round(margalef_d1, 3),
        "Menhinick_d2": round(menhinick_d2, 3),
        "Species_Table": pd.DataFrame(metrics).sort_values(by="IV_pct", ascending=False)
    }

Kết quả định lượng cấu trúc và chỉ số đa dạng

1. Đai cao 500m – 700m (Chân núi Tam Sao - Lũng Lỳ)

  • Cấu trúc tầng tán: Phân hóa thành 4 tầng rõ rệt: Tầng 1 (tán vượt ưu thế: $H = 25\text{–}35\text{ m}$, $D_{1.3} = 45\text{–}65\text{ cm}$, mật độ 40–100 cây/ha, độ tàn che 0,6–0,7), Tầng 2 (tầng dưới tán: $H = 12\text{–}20\text{ m}$), Tầng 3 (cây bụi: $H = 4\text{–}6\text{ m}$), Tầng 4 (thảm tươi và dây leo).
  • Công thức tổ thành tầng cây cao: $$29,36\text{ Ngh} + 9,90\text{ Nh} + 9,19\text{ Gvn} + 6,88\text{ Trl} + 44,68\text{ LK}$$ (Ngh: Nghiến, Nh: Nhọc, Gvn: Găng Việt Nam, Trl: Trai lý, LK: Loài khác)
[Chỉ số đa dạng tầng cây cao đai 500 - 700m]
OTC 01: D1=0.94 | d1=15.764 | d2=3.466 | H'=3.09
OTC 02: D1=0.88 | d1=6.866  | d2=1.897 | H'=2.29
OTC 03: D1=0.86 | d1=7.710  | d2=2.166 | H'=2.24
OTC 04: D1=0.86 | d1=7.542  | d2=2.081 | H'=2.19
OTC 05: D1=0.94 | d1=14.216 | d2=4.016 | H'=2.95

2. Đai cao dưới 500m (Khu vực Kéo Nàng)

  • Cấu trúc tầng tán: Cấu trúc đơn giản hóa còn 3 tầng do chịu tác động khai thác trong quá khứ. Tầng cây gỗ ($H = 10\text{–}25\text{ m}$, $D_{1.3} = 25\text{–}45\text{ cm}$, độ tàn che 0,6–0,8), tầng cây bụi và thảm tươi/dây leo.
  • Công thức tổ thành tầng cây cao: $$24,10\text{ Ddx} + 15,38\text{ Mtr} + 12,06\text{ Th} + 8,44\text{ Sb} + 7,15\text{ C} + 6,73\text{ Thc} + 26,14\text{ LK}$$ (Ddx: Dâu da xoan, Mtr: Muồng trắng, Th: Thung, Sb: Sòi bàng, C: Côm, Thc: Thôi chanh)

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

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

  1. Phân tầng sinh thái theo đai cao nghiêm ngặt: Xác định ngưỡng chuyển vị cấu trúc tại mốc 500m. Tại đai 500–700m, ưu thế tuyệt đối thuộc về các loài gỗ cứng nguyên sinh bản địa (Excentrodendron tonkinense chiếm $IV = 29,36%$). Ngược lại, đai dưới 500m bị chi phối bởi các loài cây tiên phong ưa sáng thứ sinh (Allospondias lakonensis chiếm $IV = 24,10%$).
  2. Định lượng đa dạng kết hợp 4 chỉ số chuẩn hóa: Khắc phục sai lệch của phương pháp đơn chỉ số khi phân tích quần xã rừng karst có độ phân tán không gian cao.
Tiêu chí so sánh Nghiên cứu Nam Xuân Lạc (2015) Nghiên cứu Hữu Liên (Nguyễn Nghĩa Thìn) Nghiên cứu Tân Hóa - Quảng Bình (Bùi Thế Đồi)
Độ cao khảo sát Phân tầng rõ rệt (< 500m & 500–700m) Khảo sát bình diện tổng quát 48 km² Khảo sát điểm không phân đai cao chi tiết
Chỉ số đa dạng Tích hợp $H'$ (2.19–3.09), $D_1$ (0.86–0.94), $d_1$, $d_2$ Thống kê số lượng loài và chi thuần túy Tính toán chỉ số đa dạng tầng cây gỗ chung
Phân tích tái sinh Lập 45 ODB (25 m²), tính $K_i$ và mật độ $N/\text{ha}$ Mô tả chung thảm tươi cây bụi Phân tích tái sinh chưa gắn liền với $IV%$ tầng mẹ

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

Kịch bản quản trị và phân vùng bảo tồn

  • Phân khu bảo tồn nghiêm ngặt đai cao 500–700m: Thiết lập ranh giới bảo vệ tuyệt đối các lâm phần có Nghiến và Trai lý chiếm $IV% > 20%$. Cấm hoàn toàn việc chặt hạ cây tái sinh và thu hái lâm sản phụ.
  • Phân khu phục hồi sinh thái đai dưới 500m: Thực hiện các biện pháp khoanh nuôi xúc tiến tái sinh tự nhiên có làm giàu rừng bằng các loài cây bản địa (Markhamia stipulata, Garcinia fagraeoides), tỉa thưa dây leo cạnh tranh dinh dưỡng.
+-----------------------------------------------------------------------------------+
|                  LỘ TRÌNH TRIỂN KHAI BẢO TỒN VÀ QUẢN LÝ TÀI NGUYÊN                |
+--------------------------+----------------------------+---------------------------+
| Giai đoạn 1 (Tháng 1-3)  | Giai đoạn 2 (Tháng 4-8)    | Giai đoạn 3 (Tháng 9-12)  |
| - Cắm mốc 09 OTC định vị | - Tập huấn kỹ thuật khoanh | - Đánh giá chỉ số sinh học|
| - Cài đặt GPS giám sát   |   nuôi & chuyển đổi bếp đun| - Nghiệm thu mô hình sinh |
| - Họp dân 9 thôn vùng đệm| - Giao khoán 1.646 ha rừng |   kế vùng đệm (Cây hồng)  |
+--------------------------+----------------------------+---------------------------+

Đánh giá hiệu quả kinh tế - sinh thái (Cost-Benefit & ROI)

  • Giảm áp lực củi đun: Triển khai bếp đun cải tiến tiết kiệm năng lượng cho 410 hộ vùng đệm giúp cắt giảm 40% lượng củi tiêu thụ (tiết kiệm ~1.600 tấn củi/năm), bảo vệ trực tiếp hàng chục hecta rừng tái sinh.
  • Chuyển đổi sinh kế: Hỗ trợ mở rộng diện tích Hồng không hạt (hiện có 18 ha) và thâm canh lúa ruộng nước giúp ổn định thu nhập, giảm 85% tần suất xâm nhập rừng trái phép sau 2 năm.

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

Hạn chế kỹ thuật

  1. Quy mô mẫu điều tra: Do địa hình karst tai mèo hiểm trở, số lượng OTC mới dừng lại ở 09 ô (tổng diện tích $4.500\text{ m}^2$), chưa bao phủ toàn bộ đỉnh dốc trên 700m (đỉnh Tam Sao cao 1.159m).
  2. Yếu tố chu kỳ thời gian: Dữ liệu điều tra tập trung vào mùa thu - đông (tháng 8–11/2014), chưa theo dõi biến động động thái tái sinh qua mùa mưa nóng ẩm (tháng 4–7).

Hướng nghiên cứu mở rộng

  • Ứng dụng công nghệ LiDAR quét tán rừng kết hợp máy bay không người lái (UAV) để lập bản đồ sinh khối và trữ lượng carbon ($C$) trên toàn bộ 1.788 ha.
  • Thiết lập ô định vị vĩnh viễn (Permanent Sample Plot - PSP 1 ha) để quan trắc diễn thế sinh thái dài hạn 5–10 năm.

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

+-------------------------------------------------------------------------------+
|                             ĐỐI TƯỢNG HƯỞNG LỢI                               |
+-----------------------+-------------------------------------------------------+
| Sinh viên &           | - Bộ dữ liệu thực nghiệm chuẩn mực về lâm phần karst. |
| Giảng viên Lâm nghiệp | - Giáo trình sống về phương pháp điều tra OTC & ODB.  |
+-----------------------+-------------------------------------------------------+
| Ban Quản lý KBT       | - Cơ sở dữ liệu khoa học để xây dựng phương án QLRBV. |
| & Cán bộ Kiểm lâm     | - Bản đồ ưu tiên tuần tra tại các điểm nóng sinh thái.|
+-----------------------+-------------------------------------------------------+
| Cộng đồng dân cư      | - Hưởng lợi từ chính sách chi trả DVMTR & giao khoán. |
| 9 thôn vùng đệm       | - Giảm thiểu nguy cơ sạt lở đất, lũ quét đầu nguồn suối|
+-----------------------+-------------------------------------------------------+
| Giới nghiên cứu       | - Dữ liệu so sánh đa dạng sinh học hệ núi đá Đông Bắc.|
| ĐDSH quốc tế          | - Đóng góp bằng chứng bảo tồn loài đặc hữu Việt Nam.  |
+-----------------------+-------------------------------------------------------+

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

1. Yêu cầu kỹ thuật tối thiểu để thiết lập ô tiêu chuẩn trên địa hình núi đá vôi là gì?

Cần trang bị máy định vị vệ tinh GPS có độ nhạy cao bắt sóng dưới tán lá dày, thước dây thép chuyên dụng không co giãn, thước đo chiều cao Blumeleiss và tuân thủ nguyên tắc lập ô chữ nhật $25\text{ m} \times 20\text{ m}$ bám sát đường đồng mức để giảm thiểu sai số đo đạc tiết diện ngang trên sườn dốc $> 30^\circ$.

2. Sự khác biệt cốt lõi về chỉ số sinh thái giữa đai rừng 500–700m và đai dưới 500m là gì?

Đai 500–700m có tính đa dạng loài và độ ổn định sinh thái vượt trội ($H'$ đạt tối đa 3.09, $D_1$ đạt 0.94), cấu trúc 4 tầng với sự thống trị của nhóm loài cây gỗ bản địa quý hiếm (Nghiến $IV=29,36%$). Đai dưới 500m chỉ có 3 tầng tán, chỉ số $H'$ thấp hơn và bị lấn át bởi các loài cây thứ sinh (Dâu da xoan $IV=24,10%$).

3. Làm thế nào để giải quyết mâu thuẫn giữa sinh kế người dân vùng đệm và bảo tồn vùng lõi KBT?

Áp dụng mô hình đồng quản lý: Kết hợp 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 (DVMTR), hỗ trợ kỹ thuật thâm canh nông nghiệp vùng thung lũng, cung cấp bếp đun tiết kiệm củi và quy hoạch vùng chăn thả gia súc cố định để chấm dứt thả rông trong phân khu bảo vệ nghiêm ngặt.

4. Quy trình định loại mẫu thực vật không rõ tên ngoài hiện trường được thực hiện như thế nào?

Cây chưa rõ tên được đánh mã hiệu hiện trường, thu thập cành có lá, hoa/quả (tiêu bản khô ép cồn 70°) kết hợp chụp ảnh macro cấu trúc vỏ, lá, sau đó tra cứu tài liệu chuyên khảo theo hệ thống Brummitt (1992), Cây cỏ Việt Nam (Phạm Hoàng Hộ) và đối chiếu mẫu chuẩn tại các viện nghiên cứu chuyên ngành.

5. Chi phí triển khai gói giải pháp lâm sinh và thời gian hoàn vốn sinh thái?

Tổng kinh phí ước tính cho việc lập mốc ranh giới, giám sát sinh thái và hỗ trợ bếp đun cải tiến dao động từ 150–200 triệu VNĐ/năm. Thời gian phục hồi cấu trúc tán rừng thứ sinh đai dưới 500m dự kiến từ 5–7 năm, trong khi giá trị phòng hộ đầu nguồn và bảo tồn nguồn gen quý có giá trị bền vững lâu dài.


Kết luận

Đề tài "Đánh giá đa dạng thực vật thân gỗ kiểu rừng kín thường xanh trên núi đá độ cao dưới 700m thuộc Khu bảo tồn Loài và Sinh cảnh Nam Xuân Lạc, huyện Chợ Đồn, tỉnh Bắc Kạn" đã cung cấp hệ thống số liệu định lượng chuẩn xác về cấu trúc tầng tán, tổ thành loài và tiềm năng tái sinh của hệ sinh thái rừng núi đá vôi Đông Bắc.

Kết quả nghiên cứu khẳng định tính cấp thiết phải thiết lập chế độ bảo vệ nghiêm ngặt đối với đai rừng cao từ 500–700m – nơi lưu giữ nguồn gen vô giá của các loài cây gỗ quý như Nghiến (Excentrodendron tonkinense), Trai lý (Garcinia fagraeoides). Đồng thời, gói giải pháp đồng quản lý kết hợp cải thiện sinh kế vùng đệm chính là chìa khóa then chốt để duy trì tính toàn vẹn của hệ sinh thái cảnh quan Nam Xuân Lạc trong tương lai.