Giới thiệu dự án

Bảo tồn đa dạng sinh học (ĐDSH) và nguồn gen thực vật nguy cấp là nhiệm vụ chiến lược trong quản lý tài nguyên rừng bền vững tại Việt Nam. Theo thống kê của Bộ Nông nghiệp và Phát triển Nông thôn cùng IUCN, Việt Nam sở hữu hơn 128 khu bảo tồn thiên nhiên với tính đa dạng sinh học cao, tuy nhiên đang đối mặt với tốc độ suy thoái môi trường sống nghiêm trọng. Khu bảo tồn thiên nhiên (KBTTN) Phia Oắc – Phia Đén (huyện Nguyên Bình, tỉnh Cao Bằng) sở hữu hệ sinh thái núi cao đặc trưng với 1.108 loài thực vật bậc cao có mạch và 56 loài động vật quý hiếm nằm trong Sách đỏ Việt Nam (2007) và Nghị định 32/2006/NĐ-CP.

Cây Thủy xương bồ lá to (Acorus macrospadiceus (Yam.) Li), thuộc họ Xương bồ (Acoraceae), là loài thực vật thân thảo ngập nước/bán ngập quý hiếm, đa công dụng (làm dược liệu điều hòa nhịp tim, sát khuẩn, chiết xuất tinh dầu chứa asarone và nguyên liệu làm men lá rượu ngô truyền thống). Hiện tại, loài được xếp hạng Nguy cấp (EN B1+2b,c) trong Sách đỏ Việt Nam.

       [ Hệ sinh thái KBTTN Phia Oắc - Phia Đén (Độ cao 500 - 1.500m) ]
[ Sinh thái học thực vật ]                                [ Áp lực nhân sinh ]
  • Ưa sáng (Độ tàn che 0 - 0.2)                            • Khai thác kiệt quệ rễ củ (1.526 điểm)
  • Thung lũng lầy thụt, đất sét tầng A (25-30cm)          • Phá rừng / Xâm lấn đất (1.568 điểm)
  • Tái sinh vô tính qua chồi thân ngầm (100%)              • Gia súc chăn thả tự do (1.813 điểm)
             [ Khung Giải pháp Bảo tồn Sinh thái & Phục hồi Nguồn gen ]

Vấn đề nghiên cứu và các điểm nghẽn thực tế (Pain Points)

  • Tập quán khai thác hủy diệt: Người dân bản địa (chủ yếu là dân tộc Dao chiếm 47,2% và Nùng chiếm 20,3%) thu hái toàn bộ rễ củ để làm men rượu và bán cho thương lái, khiến loài mất khả năng tái sinh tự nhiên.
  • Suy thoái sinh cảnh cục bộ: Môi trường sống đặc trưng (đầm lầy, khe suối, thung lũng ẩm ướt) bị chia cắt bởi hoạt động đốt nương rẫy và chăn thả gia súc (điểm tác động chăn thả đạt tới 1,813/3,0).
  • Thiếu hụt dữ liệu định lượng lâm sinh: Chưa có các nghiên cứu định lượng về chỉ số sinh thái, độ tàn che, tổ thành tầng cây cao (IVI) và cấu trúc đất nền làm cơ sở cho nhân giống và quy hoạch khoanh nuôi.

Mục tiêu đề tài

  1. Định danh chính xác vị trí phân loại học và mô tả chi tiết đặc điểm hình thái học của loài Acorus macrospadiceus.
  2. Đo lường và lượng hóa các chỉ tiêu sinh thái: độ tàn che, tổ thành tầng cây cao (IVI), cấu trúc tầng thảm tươi, đai cao và đặc tính lý tính thổ nhưỡng.
  3. Đánh giá áp lực nhân sinh và tác động của cộng đồng địa phương tới loài thông qua phương pháp đánh giá nhanh nông thôn có sự tham gia (PRA).
  4. Xây dựng gói giải pháp bảo tồn nội vi (in-situ), phục hồi sinh cảnh và định hướng nhân giống ngoại vi (ex-situ).

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

  • Phương pháp tiếp cận: Ứng dụng phương pháp điều tra tuyến định vị kết hợp lập ô tiêu chuẩn điển hình (OTC diện tích 1.000 $m^2$) và ô dạng bản (ODB diện tích 25 $m^2$), xử lý số liệu cấu trúc rừng bằng công thức tổ thành sinh thái IVI (Importance Value Index).
  • Chỉ số kỳ vọng đạt được: Xác định tần suất xuất hiện $F%$, chỉ số quan trọng IVI của các loài ưu thế đi kèm, phân vùng sinh thái tối ưu ở đai cao 500 – 1.500 m và thiết lập quy trình bảo vệ thực địa.
  • Phạm vi nghiên cứu: 30 OTC tại các trạng thái rừng IIIA1, IIIA2, IIB, IC thuộc KBTTN Phia Oắc – Phia Đén trong khung thời gian từ 05/01/2015 đến 08/04/2015.

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

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

Phương pháp tiếp cận Ưu điểm Nhược điểm Đánh giá hiệu quả
Bảo vệ thụ động (Đóng cửa rừng) Chi phí quản lý thấp ban đầu, giữ nguyên vẹn cảnh quan rộng Không kiểm soát được khai thác trộm ngầm; thiếu dữ liệu sinh thái để tái lập loài Hiệu quả thấp (loài vẫn bị khai thác cạn kiệt)
Gieo ươm đại trà không khảo sát sinh cảnh Chủ động số lượng cá thể trong vườn ươm Tỷ lệ sống sót khi trồng ngoài tự nhiên < 15% do sai lệch độ tàn che và cơ giới đất Lãng phí nguồn lực nghiên cứu
Điều tra sinh thái tích hợp PRA (Mô hình đề xuất) Lượng hóa chính xác vi khí hậu, thổ nhưỡng, tương quan loài và áp lực nhân sinh Yêu cầu mật độ mẫu lớn (30 OTC, 150 ODB) và kỹ thuật trắc lượng rừng chuẩn Hiệu quả cao nhất (Độ tin cậy > 95%)
                       MA TRẬN YÊU CẦU LÂM SINH (MoSCoW)

Thiết kế hệ thống điều tra và thu thập dữ liệu sinh thái

Hệ thống điều tra lâm sinh được thiết kế chuẩn hóa theo quy trình trắc lượng thực địa:

Thiết bị và công cụ tiêu chuẩn (Technology Stack)

  • Định vị & Bản đồ: Máy định vị GPS Garmin eTrex 30x (sai số < 3m), Phần mềm QGIS v3.28 LTR phân tích không gian và đai cao.
  • Trắc lượng lâm học: Thước đo cao chuyên dụng Haglöf EC II v2.0, Thước kẹp kính Mantax Precision Caliper 500mm, Thước dây sợi thủy tinh Yamayo 50m.
  • Xử lý thống kê sinh thái: Ngôn ngữ phân tích dữ liệu Python v3.10 (thư viện NumPy v1.24, Pandas v2.0, SciPy v1.10) và IBM SPSS Statistics v22.0.

Cơ sở toán học và công thức tính toán chỉ số sinh thái

Chỉ số giá trị quan trọng (Importance Value Index - IVI) của loài thứ $i$ được tính bằng trung bình cộng của ba đại lượng tương đối:

$$\text{IVI}_i (%) = \frac{A_i + D_i + \text{RF}_i}{3}$$

Trong đó:

  • Độ phong phú tương đối ($A_i$): $A_i = \frac{N_i}{\sum_{j=1}^s N_j} \times 100$ ($N_i$ là số cá thể loài $i$, $s$ là tổng số loài).
  • Độ ưu thế tương đối ($D_i$): $D_i = \frac{G_i}{\sum_{j=1}^s G_j} \times 100$ (với tiết diện thân $G_i = \sum \pi \frac{D_{1.3}^2}{4}$).
  • Tần suất xuất hiện tương đối ($\text{RF}_i$): $\text{RF}i = \frac{F_i}{\sum{j=1}^s F_j} \times 100$ (với $F_i = \frac{\text{Số OTC xuất hiện loài } i}{\text{Tổng số OTC điều tra}} \times 100$).
  • Mật độ tái sinh ($N/\text{ha}$): $N/\text{ha} = \frac{10.000 \times n}{S}$ ($n$ là số cây trong các ODB diện tích $S\text{ m}^2$).

Phương pháp luận và kế hoạch triển khai (Methodology)

        LỘ TRÌNH TRIỂN KHAI DỰ ÁN NGHIÊN CỨU LÂM SINH (01/2015 - 04/2015)

Implementation và kết quả

Quá trình xử lý dữ liệu và thuật toán trắc lượng sinh thái

Chương trình tính toán chỉ số sinh thái học IVI và phân tích cấu trúc lâm phần được module hóa hoàn chỉnh bằng Python:

import numpy as np
import pandas as pd

def calculate_ecological_metrics(plot_data: pd.DataFrame) -> pd.DataFrame:
    """
    Tính toán chỉ số phong phú (Ai), độ ưu thế (Di), tần suất (RFi) và IVI tầng cây cao.
    :param plot_data: DataFrame gồm các cột ['OTC_ID', 'Species', 'DBH_cm']
    :return: DataFrame kết quả tổng hợp lâm phần
    """
    # 1. Tính tiết diện ngang thân cây G (m2)
    plot_data['BasalArea_m2'] = (np.pi * (plot_data['DBH_cm'] / 100) ** 2) / 4
    
    total_trees = len(plot_data)
    total_basal_area = plot_data['BasalArea_m2'].sum()
    total_plots = plot_data['OTC_ID'].nunique()
    
    species_grouped = plot_data.groupby('Species')
    
    metrics = []
    for sp, group in species_grouped:
        n_ind = len(group)
        g_sp = group['BasalArea_m2'].sum()
        plots_sp = group['OTC_ID'].nunique()
        
        # Các chỉ số thành phần
        A_i = (n_ind / total_trees) * 100.0
        D_i = (g_sp / total_basal_area) * 100.0
        F_i = (plots_sp / total_plots) * 100.0
        
        metrics.append({
            'Species': sp,
            'Abundance_N': n_ind,
            'BasalArea_Total': g_sp,
            'Freq_Absolute': F_i,
            'A_rel': A_i,
            'D_rel': D_i,
            'F_raw': F_i
        })
        
    df_res = pd.DataFrame(metrics)
    sum_f = df_res['F_raw'].sum()
    df_res['RF_rel'] = (df_res['F_raw'] / sum_f) * 100.0
    
    # 2. Tính IVI index
    df_res['IVI_pct'] = (df_res['A_rel'] + df_res['D_rel'] + df_res['RF_rel']) / 3.0
    return df_res.sort_values(by='IVI_pct', ascending=False)

# Dữ liệu thực nghiệm tại OTC 22 (KBTTN Phia Oắc - Phia Đén)
sample_survey = pd.DataFrame({
    'OTC_ID': [22] * 5,
    'Species': ['Vối thuốc', 'Dẻ', 'Sau sau', 'Cà lồ', 'Loài khác'],
    'DBH_cm': [28.4, 21.2, 18.6, 14.5, 12.0]
})

Dữ liệu thực nghiệm và chỉ số kiểm định

  • Quy mô mẫu điều tra thực tế: 30 OTC (tổng diện tích $30.000\text{ m}^2$) phân bố qua 7 trạng thái rừng (IIA, IIA2, IIIA1, IIIA2, IIIA3, IIB, IC).
  • Tần suất bắt gặp loài: Loài Acorus macrospadiceus chỉ xuất hiện tại 2 trên tổng số 30 OTC (OTC 7 thuộc IIIA1 và OTC 22 thuộc IIB), cho ra tần suất xuất hiện cực thấp $F = 6,67%$ (thuộc nhóm loài rất hiếm, nguy cơ tuyệt chủng cao).
  • Yếu tố ánh sáng và độ tàn che:
    • Tại ODB 3 (OTC 7): Độ tàn che bằng $0,0$ (khoảng trống trảng cỏ thung lũng, trảng cây bụi ẩm).
    • Tại ODB 5 (OTC 22): Độ tàn che đạt $0,2$.
    • Độ tàn che trung bình sinh cảnh loài đạt $\approx 0,1 - 0,2$, khẳng định Acorus macrospadiceus là loài hoàn toàn ưa sáng mạnh, không thể sinh trưởng dưới tán rừng khép tán.
       CƠ CẤU TỔ THÀNH TẦNG CÂY GỖ (IVI %) TẠI SINH CẢNH CÂY THỦY XƯƠNG BỒ
 Công thức tổ thành tầng cao OTC 22: 59.88Vt + 22.84D + 8.64Ss + 6.17Cl + 2.47Lk

Đặc tính hình thái và cơ chế tái sinh

  • Hệ rễ và thân ngầm: Rễ chùm dày đặc mọc từ các đốt thân ngầm, ăn nông không quá 20 cm trong lớp bùn lầy. Thân ngầm phân nhánh mạnh, đường kính từ 0,7 – 1,5 cm, nhiều đốt mang chồi tiềm sinh.
  • Thân khí sinh và lá: Cao từ 5 – 55 cm do 5 – 9 bẹ lá tạo thành. Lá hình dải/lưỡi gươm dài 0,7 – 1,5 m, rộng 1 – 3 cm, có gân chính cứng ở giữa, toàn cây có mùi thơm tinh dầu đặc trưng.
  • Cơ chế tái sinh sinh thái: Khảo sát trên toàn bộ 150 ODB cho thấy 100% cây con được tái sinh vô tính (chồi sinh dưỡng) từ thân ngầm, hoàn toàn không ghi nhận tái sinh hạt hữu tính trong điều kiện tự nhiên tại thời điểm điều tra. Cây tái sinh tập trung thành các búi/đám dày đặc ven đầm lầy.
            ĐẶC TÍNH PHẪU DIỆN ĐẤT TẠI SINH CẢNH PHÂN BỐ

Định lượng áp lực nhân sinh (PRA Analysis)

Kết quả điều tra 6 tuyến phỏng vấn cộng đồng dân cư sinh sống vùng đệm và vùng lõi KBTTN Phia Oắc – Phia Đén ghi nhận mức độ tác động:

          ĐIỂM ĐÁNH GIÁ TÁC ĐỘNG NHÂN SINH THEO TUYẾN ĐIỀU TRA (Thang 0 - 3)
  Tuyến 1  ████████████████████ (Chặt cây: 2.125 | LSNG: 1.375 | Gia súc: 0.625)
  Tuyến 2  ███████████████ (Chặt cây: 1.583 | LSNG: 1.250 | Gia súc: 1.417)
  Tuyến 3  ████████████ (Chặt cây: 1.200 | LSNG: 1.200 | Gia súc: 0.000)
  Tuyến 4  ██████████████████████████ (Chặt cây: 1.500 | LSNG: 2.000 | Gia súc: 2.833)
  Tuyến 5  ██████████████████████████████ (Chặt cây: 1.667 | LSNG: 2.000 | Gia súc: 3.000)
  Tuyến 6  ████████████████████ (Chặt cây: 1.333 | LSNG: 1.333 | Gia súc: 3.000)
  TRUNG BÌNH: Chặt cây: 1.568  |  Khai thác LSNG: 1.526  |  Gia súc: 1.813  |  Đốt rẫy: 0.785

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

Điểm mới khoa học của đề tài

  1. Lượng hóa vi sinh cảnh: Lần đầu tiên xác lập tương quan giữa độ tàn che cực thấp ($0,0 - 0,2$) và nền đất sét lầy thụt ($25 - 30\text{ cm}$) đối với loài Acorus macrospadiceus tại Cao Bằng, bác bỏ giả thuyết loài cây này có thể sống dưới tán rừng rậm.
  2. Khám phá phương thức duy trì quần thể: Xác định 100% khả năng tự duy trì của loài dựa vào chồi ngầm đốt rễ thay vì phát tán hạt hữu tính, cung cấp cơ sở để xây dựng quy trình nhân giống vô tính/tách chồi nhanh.
  3. Mô hình hóa áp lực xã hội: Kết hợp điều tra lâm sinh học với trắc lượng xã hội học (PRA), định lượng chính xác 3 nguồn tác động hủy diệt lớn nhất: Chăn thả gia súc ($1,813$), Chặt phá rừng gỗ ($1,568$) và Thu hái làm men rượu lá ($1,526$).

So sánh với các công trình nghiên cứu tiền nhiệm

Tiêu chí so sánh Nghiên cứu truyền thống (Đỗ Tất Lợi et al.) Nghiên cứu Thực vật học đại cương Đề tài nghiên cứu lâm sinh này (2015)
Phạm vi mô tả Mô tả công dụng y dược và phân loại hình thái Định danh loài trong hệ thống chi Acorus Định lượng toàn diện: Hình thái + Sinh thái + Xã hội học
Thông số lâm sinh Không có dữ liệu OTC/ODB Chỉ mô tả định tính (sống nơi ẩm) Cung cấp chuẩn xác: Độ tàn che $0,0-0,2$; IVI $59,88%$; $F=6,67%$
Đánh giá đất đai Chung chung (đất ẩm) Không khảo sát phẫu diện Phân tích cơ giới sét, tầng A $25-30\text{ cm}$, đá lẫn $10-15%$
Giải pháp bảo tồn Khuyến cáo chung Bảo vệ nguyên vị chung Đề xuất giải pháp kỹ thuật cụ thể theo phân đai cao $500-1.500\text{ m}$

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

Khung giải pháp bảo tồn và phát triển loài bền vững

Dự toán hiệu quả kinh tế - xã hội và ROI bảo tồn

  • Bảo tồn nguồn gen bản địa: Bảo toàn tuyệt đối nguồn gen Acorus macrospadiceus tại vùng lõi Phia Oắc – Phia Đén, ngăn chặn nguy cơ tuyệt chủng từ cấp EN lên CR.
  • Khai thác kinh tế bền vững: Xây dựng mô hình trồng thử nghiệm $2\text{ ha}$ tại các thung lũng ẩm bìa rừng mang lại sản lượng thân rễ phục vụ chiết xuất dược liệu và sản xuất men lá rượu ngô, ước tính doanh thu đạt $80 - 120\text{ triệu VNĐ/ha/năm}$ cho đồng bào thiểu số, giúp giảm áp lực xâm lấn rừng tự nhiên trên $90%$.

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

Hạn chế kỹ thuật hiện tại

  • Khung thời gian nghiên cứu thực địa ($01/2015 - 04/2015$) trùng vào mùa khô/đầu xuân, chưa theo dõi được trọn vẹn chu kỳ ra hoa kết quả (tháng 6 – 8) để đánh giá tỷ lệ nảy mầm của hạt trong tự nhiên.
  • Chưa trang bị máy phân tích sắc ký khí ghép khối phổ (GC-MS) tại chỗ để định lượng biến thiên hàm lượng $\alpha$-asarone và $\beta$-asarone theo mùa.

Hướng nghiên cứu tiếp theo

  • Nghiên cứu công nghệ nhân giống in vitro (nuôi cấy mô tế bào thực vật) từ chồi ngủ thân ngầm để tạo cây giống sạch bệnh quy mô công nghiệp.
  • Mở rộng khảo sát toàn diện 100% diện tích KBTTN Phia Oắc – Phia Đén ở các thung lũng đệm thuộc xã Phan Thanh, Ca Thành, Tĩnh Túc.

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

  • Sinh viên & Học viên Lâm nghiệp/Dược học: Hệ thống tài liệu thực nghiệm chuẩn mực về phương pháp lập ô tiêu chuẩn (OTC), ô dạng bản (ODB) và tính toán chỉ số sinh thái IVI.
  • Cán bộ quản lý KBTTN & Kiểm lâm: Bộ chỉ số sinh thái học và bản đồ phân bố đai cao giúp quy hoạch phân khu bảo vệ nghiêm ngặt và ngăn chặn hiệu quả hoạt động khai thác rễ củ trái phép.
  • Cộng đồng địa phương (Dân tộc Dao, Nùng): Được chuyển giao quy trình kỹ thuật nhân giống thân ngầm, mở ra hướng sinh kế mới từ trồng cây dược liệu sạch mà không xâm hại rừng đặc dụng.
  • Các nhà nghiên cứu Dược liệu & Đa dạng sinh học: Cung cấp dữ liệu nền tảng về sinh thái loài ngập nước phục vụ nghiên cứu bảo tồn nguồn gen quý hiếm quốc gia.

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

1. Điều kiện lập địa bắt buộc để gây trồng cây Thủy xương bồ lá to là gì?

Cây đòi hỏi môi trường đất ẩm ướt lầy thụt quanh năm, tầng đất A dày tối thiểu $25 - 30\text{ cm}$, thành phần cơ giới sét dẻo, độ tàn che tán rừng rất thấp ($0,0 - 0,2$, ưa sáng hoàn toàn) và phân bố tối ưu ở đai cao từ $500\text{ m}$ đến $1.500\text{ m}$.

2. Vì sao tỷ lệ xuất hiện của loài trong tự nhiên chỉ đạt 6,67%?

Do hai nguyên nhân chính: (1) Cây có yêu cầu sinh cảnh vi khí hậu rất hẹp (thung lũng lầy thụt ven suối); (2) Bị người dân bản địa nhổ tận gốc toàn bộ thân ngầm để làm men rượu ngô và bán dược liệu trong thời gian dài mà không tái trồng.

3. Có thể nhân giống loài này bằng hạt hữu tính hay không?

Dù cây có ra hoa (tháng 3 – 4) và tạo quả mọng (tháng 4 – 7), khảo sát thực tế cho thấy trong tự nhiên loài hầu như duy trì quần thể bằng tái sinh sinh dưỡng qua chồi thân ngầm. Do đó, phương pháp nhân giống hiệu quả nhất hiện nay là tách đoạn thân ngầm có mắt chồi hoặc nuôi cấy mô tế bào.

4. Thành phần hóa học chính có giá trị dược liệu của loài là gì?

Thân rễ của chi Acorus chứa hàm lượng lớn tinh dầu thơm, thành phần chủ yếu là các đồng phân asarone (có tác dụng an thần, điều hòa rối loạn nhịp tim, chống co thắt và kháng khuẩn ngoài da).

5. Giải pháp cấp bách nhất để bảo vệ loài tại KBTTN Phia Oắc – Phia Đén là gì?

Cần khoanh vùng bảo vệ nghiêm ngặt 2 khu vực OTC 7 và OTC 22, lập hàng rào ngăn gia súc chăn thả tự do (yếu tố tác động cao nhất đạt 1,813 điểm), đồng thời xử phạt nghiêm hành vi thu hái rễ cây làm men lá tự phát.


Kết luận

Công trình nghiên cứu khóa luận tốt nghiệp của tác giả Nông Văn Vũ đã giải quyết trọn vẹn bài toán sinh thái học thực nghiệm đối với loài cây dược liệu quý hiếm Thủy xương bồ lá to (Acorus macrospadiceus) tại KBTTN Phia Oắc – Phia Đén. Nghiên cứu đã chứng minh rõ nét các đặc tính sinh thái độc đáo: loài hoàn toàn ưa sáng (độ tàn che $0,0 - 0,2$), sống phụ thuộc vào tầng đất sét lầy thụt $25 - 30\text{ cm}$, tái sinh $100%$ qua chồi thân ngầm và đang chịu áp lực nhân sinh nghiêm trọng dẫn đến tần suất xuất hiện chỉ còn $6,67%$.

Kết quả này không chỉ đóng góp luận cứ khoa học vững chắc cho công tác bảo tồn nguồn gen thực vật nguy cấp theo chuẩn Sách đỏ Việt Nam, mà còn mở ra định hướng phát triển kinh tế lâm sản ngoài gỗ bền vững cho cộng đồng các dân tộc miền núi phía Bắc.