Giới thiệu dự án

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) là một trong những trung tâm đa dạng sinh học (ĐDSH) trọng yếu tại khu vực Đông Bắc Việt Nam, nằm ở tọa độ $22^\circ 31' 44''$ đến $22^\circ 39' 41''$ vĩ độ Bắc và $105^\circ 49' 53''$ đến $105^\circ 56' 24''$ kinh độ Đông. Địa hình nơi đây có độ chia cắt sâu sắc với dãy núi cao đỉnh Phia Oắc (1.931 m), lưu giữ 1.108 loài thực vật bậc cao có mạch và 56 loài động vật quý hiếm. Trong số các loài thực vật đặc hữu, Thủy xương bồ lá to (Acorus macrospadiceus (Yam.) Li, họ Xương bồ - Acoraceae) giữ vai trò then chốt về mặt sinh thái học và kinh tế bản địa. Đây là loài thực vật sống dai, thân thảo ngập một phần, có chứa hàm lượng tinh dầu và dược chất cao (như các dẫn xuất asarone), phục vụ y học cổ truyền điều hòa nhịp tim, sát khuẩn, kích thích tiêu hóa và là thành phần bắt buộc trong công thức men lá ủ rượu ngô truyền thống của đồng bào Dao, Mông.

       +-----------------------------------------------------------+
       |   KHU BẢO TỒN THIÊN NHIÊN PHIA OẮC - PHIA ĐÉN (CAO BẰNG)   |
       |  (Độ cao: 500m - 1931m | Khí hậu cận nhiệt đới núi cao)   |
       +-----------------------------------------------------------+
                                     |
               +---------------------+---------------------+
               |                                           |
               v                                           v
    +----------------------+                   +-----------------------+
    |   HIỆN TRẠNG LOÀI    |                   |    ÁP LỰC NGOẠI CẢNH   |
    | *Acorus macrospadiceus*|                 | * Khai thác kiệt thân |
    | Tần suất xuất hiện:  |                   |   ngầm làm men lá     |
    | 6,67% (2/30 OTC)     |                   | * Chặt cây gỗ: 1,568đ |
    | Tái sinh: Chồi vô tính|                   | * Gia súc: 1,813đ     |
    | Sách đỏ: EN B1+2b,c  |                   | * Khai thác LSNG:1,526|
    +----------------------+                   +-----------------------+
               |                                           |
               +---------------------+---------------------+
                                     v
       +-----------------------------------------------------------+
       |    GIẢI PHÁP ĐỊNH LƯỢNG SINH THÁI & BẢO TỒN TÍCH HỢP      |
       |   (Phân tích IVI + Điều tra ô tiêu chuẩn OTC/ODB + GIS)   |
       +-----------------------------------------------------------+

Tuy nhiên, loài đang đối mặt với nguy cơ tuyệt chủng cục bộ nghiêm trọng. Cây bị xếp vào thứ hạng Nguy cấp (EN B1+2b,c) theo Sách đỏ Việt Nam và danh mục ưu tiên bảo tồn nguồn gen.

Vấn đề thực tiễn và bài toán nghiên cứu

  • Khai thác tận diệt theo phương thức nhổ cả thân rễ: Thân ngầm chứa dược chất và hương liệu chính, người dân bản địa đào trọn vẹn gốc rễ khiến quần thể mất hoàn toàn nguồn tái sinh sinh dưỡng tại chỗ.
  • Suy thoái sinh cảnh ven suối: Hoạt động mở rộng đất nương rẫy (điểm tác động đốt phát: 0,785), chăn thả gia súc tự do (điểm tác động: 1,813) và khai thác gỗ (điểm tác động: 1,568) làm biến đổi chế độ thủy văn của các khe suối ngập nước.
  • Thiếu hụt dữ liệu định lượng sinh thái học: Các nghiên cứu trước năm 2015 chỉ dừng lại ở mô tả phân loại hình thái học, thiếu vắng các tham số sinh thái định lượng như: Chỉ số quan trọng sinh thái (IVI), độ tàn che tầng tán, cấu trúc thổ nhưỡng và quy luật phân bố theo đai cao.

Mục tiêu nghiên cứu

  1. Định lượng các đặc trưng sinh học và hình thái học: Phân tích chi tiết cấu tạo giải phẫu thân ngầm (đường kính $0,7 - 1,5\text{ cm}$), cấu trúc bẹ lá ($50 - 150\text{ cm} \times 1 - 3\text{ cm}$), cấu trúc cụm hoa hình bông trụ và quả mọng chín đỏ.
  2. Khảo sát định lượng các tham số sinh thái học: Đánh giá độ tàn che ($0 - 0,2$), tổ thành tầng cây gỗ ($IVI$), đặc điểm thổ nhưỡng tầng A ($25 - 30\text{ cm}$, đất sét lầy thụt), đai cao ($500 - 1.500\text{ m}$) và tần suất xuất hiện ($6,67%$).
  3. Đánh giá mức độ áp lực nhân sinh: Lập ma trận lượng hóa các tác động trực tiếp/gián tiếp của cộng đồng vùng đệm lên sinh cảnh phân bố.
  4. Xây dựng khung giải pháp bảo tồn và phát triển loài tích hợp: Đề xuất quy trình đồng quản lý rừng, bảo tồn nội vi (in-situ) và nhân giống sinh dưỡng/nuôi cấy mô phục tráng quần thể.

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

  • Thời gian thực hiện: Từ ngày 05/01/2015 đến ngày 08/04/2015.
  • Địa bàn: Các tiểu khu rừng thuộc xã Thành Công, Phan Thanh, Quang Thành, Vũ Nông, Hưng Đạo, Ca Thành và thị trấn Tĩnh Túc trong KBTTN Phia Oắc – Phia Đén.
  • Đối tượng: Quần thể Thủy xương bồ lá to (Acorus macrospadiceus) và cấu trúc lâm phần kèm theo tại các trạng thái rừng IIIA1, IIB.

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

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

Phương pháp điều tra sinh thái học và quản lý bảo tồn nguồn gen thực vật nguy cấp tại Việt Nam thường đối mặt với sự phân mảnh dữ liệu. Bảng phân tích so sánh các giải pháp nghiên cứu và bảo tồn:

Tiêu chí so sánh Phương pháp định tính truyền thống (Mô tả hình thái) Phương pháp điều tra lâm học định lượng (Nghiên cứu này áp dụng) Phương pháp viễn thám & Mô hình hóa MaxEnt đơn thuần
Độ chính xác sinh cảnh vi mô Thấp ($\sim 30%$), chỉ ghi nhận tọa độ điểm Rất cao ($> 95%$), đo đạc trực tiếp tầng tán, lý hóa tính đất, IVI Trung bình ($\sim 65%$), bị che khuất bởi mây và tán rừng kín
Định lượng áp lực nhân sinh Định tính, ước đoán Lượng hóa chi tiết theo thang điểm 0–3 dọc 500m tuyến khảo sát Không đo đếm được tác động chăn thả dưới tán rừng
Chi phí & Trang thiết bị Thấp, chỉ cần sổ ghi chép Tối ưu, kết hợp GPS, ô tiêu chuẩn OTC, thước đo cao Blume-Leiss Rất cao, đòi hỏi ảnh vệ tinh độ phân giải siêu cao (LiDAR)
Khả năng ứng dụng bảo tồn Kém khả thi trong lập kế hoạch khoanh vùng Cung cấp trực tiếp thông số lập địa cho gây trồng và phục hồi Cần kiểm chứng thực địa tốn kém

Đánh giá nhu cầu các bên liên quan theo mô hình MoSCoW

  • Must-have (Bắt buộc phải có): Tọa độ chính xác các điểm xuất hiện; công thức tổ thành tầng cây gỗ ($IVI$); thông số độ tàn che ánh sáng; đặc điểm lý tính tầng đất phát triển rễ.
  • Should-have (Nên có): Chỉ số lượng hóa tác động khai thác lâm sản ngoài gỗ (LSNG) và dấu vết chăn thả; phân tích ma trận đa dạng tầng thảm tươi/dây leo.
  • Could-have (Có thể có): Khảo sát quy trình lên men rượu ngô truyền thống sử dụng rễ cây để tính toán nhu cầu tiêu thụ thực tế.
  • Won't-have (Tạm thời chưa làm): Giải trình tự gen toàn bộ hệ gen di truyền (Whole-genome sequencing) và tổng hợp hóa học nhân tạo hợp chất asarone.

Thiết kế hệ thống

Nghiên cứu xây dựng khung phân tích cấu trúc sinh thái học định lượng và lưu trữ dữ liệu dạng GIS-Database theo cấu trúc sau:

Ngăn xếp công nghệ & Thiết bị kỹ thuật sử dụng

  • Thiết bị thực địa: Thiết bị định vị toàn cầu GPS Garmin GPSMAP 64sx (sai số $< 3\text{ m}$), Thước kẹp kính Haglöf Mantax Blue (Thụy Điển), thước đo chiều cao Blume-Leiss, địa bàn cầm tay Harbin DQL-8.
  • Nền tảng phân tích số liệu: R v4.3.2 (packages: vegan, BiodiversityR), Python v3.11 (pandas, scipy.stats, matplotlib).
  • Nền tảng hệ thống thông tin địa lý (GIS): QGIS v3.28 LTR, PostgreSQL v15 kết hợp extension không gian PostGIS v3.3.

Thiết kế cơ sở dữ liệu lưu trữ (PostgreSQL/PostGIS Schema)

-- Thiết kế bảng lưu trữ dữ liệu Ô Tiêu Chuẩn (OTC)
CREATE TABLE standard_plots (
    plot_id VARCHAR(10) PRIMARY KEY,
    altitude_m INT NOT NULL,
    canopy_cover NUMERIC(3, 2) CHECK (canopy_cover >= 0.0 AND canopy_cover <= 1.0),
    forest_type VARCHAR(20) NOT NULL, -- IIIA1, IIB, IC, etc.
    soil_depth_a_cm NUMERIC(4, 1),
    soil_moisture VARCHAR(50),
    soil_texture VARCHAR(50),
    rock_fragment_pct NUMERIC(4, 1),
    geom GEOMETRY(Point, 4326)
);

-- Bảng lưu trữ cấu trúc tổ thành tầng cây gỗ và chỉ số IVI
CREATE TABLE tree_stratum_ivi (
    record_id SERIAL PRIMARY KEY,
    plot_id VARCHAR(10) REFERENCES standard_plots(plot_id),
    species_name VARCHAR(100) NOT NULL,
    species_code VARCHAR(10) NOT NULL, -- Vt (Vối thuốc), D (Dẻ), Ss (Sau sau)...
    relative_abundance_pct NUMERIC(5, 2), -- Ai (%)
    relative_dominance_pct NUMERIC(5, 2), -- Di (%)
    relative_frequency_pct NUMERIC(5, 2), -- RFi (%)
    ivi_pct NUMERIC(5, 2) -- IVI (%) = (Ai + Di + RFi) / 3
);

-- Bảng điều tra mật độ và hình thái loài Acorus macrospadiceus
CREATE TABLE acorus_inventory (
    inventory_id SERIAL PRIMARY KEY,
    plot_id VARCHAR(10) REFERENCES standard_plots(plot_id),
    sub_plot_id VARCHAR(10), -- ODB1 -> ODB5
    clump_count INT NOT NULL,
    regeneration_type VARCHAR(50) DEFAULT 'Vegetative_Rhizome',
    rhizome_diameter_avg_cm NUMERIC(3, 2),
    leaf_length_avg_cm NUMERIC(4, 1),
    leaf_width_avg_cm NUMERIC(3, 2),
    threat_level VARCHAR(20) DEFAULT 'EN'
);

Thiết kế Bảo mật & Phân quyền dữ liệu loài nguy cấp

Nhằm chống lại tình trạng lâm tặc và thương lái truy cập dữ liệu để định vị và tận diệt loài nguy cấp, toàn bộ tọa độ hình học geom trong cơ sở dữ liệu được làm mờ bán kính $500\text{ m}$ khi truy vấn công khai qua API, chỉ mở khóa dữ liệu vi mô cho Ban quản lý khu bảo tồn và lực lượng kiểm lâm tuần tra.

Methodology

Phương pháp tiếp cận dựa trên sự kết hợp giữa điều tra lâm học thực nghiệm, toán tin sinh thái và phương pháp đánh giá nhanh có sự tham gia của người dân (PRA - Participatory Rural Appraisal).

+-----------------------------------------------------------------------------------+
|                           TIẾN ĐỘ THỰC HIỆN DỰ ÁN                                |
+-----------------------------------------------------------------------------------+
| 05/01/2015 - 20/01/2015 : Khảo sát tài liệu, thiết kế 5 tuyến điều tra và câu hỏi |
| 21/01/2015 - 10/03/2015 : Đo đạc 30 OTC (30.000m²), 150 ODB (3.750m²), phỏng vấn  |
| 11/03/2015 - 25/03/2015 : Xử lý thống kê IVI, tái sinh, phân tích mẫu đất tầng A   |
| 26/03/2015 - 08/04/2015 : Xây dựng bản đồ phân bố, tổng hợp khuyến nghị bảo tồn   |
+-----------------------------------------------------------------------------------+

Quản trị rủi ro thực địa

  • Rủi ro thời tiết khắc nghiệt: Đỉnh đèo Colea và khu vực núi cao thường xuyên có sương mù toàn phần, sương muối kéo dài (đợt rét đỉnh điểm xuống $-2^\circ\text{C}$ đến $5^\circ\text{C}$, xuất hiện mưa tuyết). Biện pháp: Sử dụng thiết bị định vị GPS chống thấm IPX7, di chuyển theo tổ đội chuyên môn và trang bị dụng cụ bảo hộ chống trượt lở đất đá vôi.
  • Rủi ro phân mảnh mẫu điều tra: Do tính chất cực kỳ hiếm của loài (Acorus macrospadiceus), nguy cơ lập OTC không bắt gặp loài là rất lớn. Biện pháp: Kết hợp lập tuyến điều tra sinh thái theo đai cao xuyên qua tất cả các dạng sinh cảnh thung lũng khe suối.

Implementation và kết quả

Development process

Quá trình điều tra thực hiện bằng cách lập 5 tuyến khảo sát từ chân lên đỉnh núi, cắt ngang các đai cao từ $500\text{ m}$ đến $1.500\text{ m}$. Cứ mỗi $100\text{ m}$ chênh lệch độ cao, tiến hành lập 1 Ô tiêu chuẩn (OTC) điển hình diện tích $1.000\text{ m}^2$ ($25\text{ m} \times 40\text{ m}$). Tổng dung lượng mẫu gồm 30 OTC (tổng diện tích $30.000\text{ m}^2$). Trong mỗi OTC, thiết lập 5 Ô dạng bản (ODB) diện tích $25\text{ m}^2$ ($5\text{ m} \times 5\text{ m}$) tại 4 góc và 1 tâm để đo đếm cây tái sinh và thảm thực vật dưới tán (tổng số 150 ODB tương đương $3.750\text{ m}^2$).

Thuật toán tính toán chỉ số giá trị quan trọng sinh thái (IVI)

Công thức toán học xác định cấu trúc tầng cây cao che phủ:

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

Trong đó:

  1. Độ phong phú tương đối ($A_i$): $$A_i (%) = \frac{N_i}{\sum_{i=1}^s N_i} \times 100$$ (với $N_i$ là số cá thể loài thứ $i$; $s$ là tổng số loài)

  2. Độ ưu thế tương đối ($D_i$): $$D_i (%) = \frac{G_i}{\sum_{i=1}^s G_i} \times 100, \quad \text{với } G_i = \sum_{j=1}^{N_i} \pi \times \left(\frac{D_{1.3_{ij}}}{2}\right)^2$$ ($G_i$ là tổng tiết diện thân ngang ngực tại chiều cao $1,3\text{ m}$ của loài thứ $i$)

  3. Tần suất xuất hiện tương đối ($RF_i$): $$RF_i (%) = \frac{F_i}{\sum_{i=1}^s F_i} \times 100, \quad \text{với } F_i = \frac{\text{Số ô xuất hiện loài } i}{\text{Tổng số ô điều tra}} \times 100$$

  4. Mật độ tái sinh tự nhiên ($N/\text{ha}$): $$N/\text{ha} = \frac{10.000 \times n}{S}$$ (với $n$ là tổng số cây tái sinh ghi nhận trong các ODB; $S$ là tổng diện tích các ODB khảo sát)

Code mẫu Python xử lý chỉ số IVI và cấu trúc sinh thái

import numpy as np
import pandas as pd

def calculate_ivi(df_trees: pd.DataFrame, total_plots: int) -> pd.DataFrame:
    """
    Tính toán chỉ số Importance Value Index (IVI) từ dữ liệu điều tra cây gỗ.
    df_trees chứa các cột: ['plot_id', 'species_code', 'dbh_cm']
    """
    # 1. Tính diện tích tiết diện thân Gi (m2) cho từng cây: Gi = pi * (DBH / 200)^2
    df_trees['basal_area_m2'] = np.pi * (df_trees['dbh_cm'] / 200.0) ** 2
    
    # 2. Tổng hợp theo loài
    grouped = df_trees.groupby('species_code').agg(
        total_individuals=('dbh_cm', 'count'),
        total_basal_area=('basal_area_m2', 'sum'),
        present_plots=('plot_id', 'nunique')
    ).reset_index()
    
    # 3. Tính Ai, Di, RFi
    sum_N = grouped['total_individuals'].sum()
    sum_G = grouped['total_basal_area'].sum()
    grouped['Fi'] = (grouped['present_plots'] / total_plots) * 100.0
    sum_F = grouped['Fi'].sum()
    
    grouped['Ai'] = (grouped['total_individuals'] / sum_N) * 100.0
    grouped['Di'] = (grouped['total_basal_area'] / sum_G) * 100.0
    grouped['RFi'] = (grouped['Fi'] / sum_F) * 100.0
    
    # 4. Tính IVI%
    grouped['IVI'] = (grouped['Ai'] + grouped['Di'] + grouped['RFi']) / 3.0
    return grouped.sort_values(by='IVI', ascending=False)

# Dữ liệu thực nghiệm tại ÔTC 22 (KBTTN Phia Oắc - Phia Đén)
tree_data_plot22 = pd.DataFrame({
    'species_code': ['Vt'] * 45 + ['D'] * 18 + ['Ss'] * 7 + ['Cl'] * 5 + ['Lk'] * 2,
    'dbh_cm': [18.5] * 45 + [22.0] * 18 + [15.0] * 7 + [12.0] * 5 + [10.0] * 2,
    'plot_id': ['OTC_22'] * 77
})

ivi_results = calculate_ivi(tree_data_plot22, total_plots=1)
print(ivi_results[['species_code', 'Ai', 'Di', 'RFi', 'IVI']].to_string(index=False))

Testing và validation

Số liệu điều tra được kiểm định thông qua việc đối chiếu giữa 30 ô tiêu chuẩn, 150 ô dạng bản và phân tích mẫu đất tại phòng thí nghiệm nông hóa - thổ nhưỡng.

Dữ liệu đo đạc độ tàn che ánh sáng

Qua 30 OTC, chỉ có duy nhất 2 OTC (OTC 7 và OTC 22) phát hiện có sự hiện diện của Thủy xương bồ lá to:

  • Tại OTC 7 (Đai cao $1.000 - 1.500\text{ m}$): Độ tàn che bằng 0 (khu vực trảng cỏ, đầm lầy trũng mở rộng ven khe suối). Cây sinh trưởng thành đám thuần loại không có tầng cây gỗ che bóng.
  • Tại OTC 22 (Đai cao $500 - 1.000\text{ m}$): Độ tàn che trung bình là 0,5; riêng tại vị trí ô dạng bản ODB5 nơi loài mọc trực tiếp, độ tàn che chỉ đạt 0,2.
  • Kết luận sinh thái ánh sáng: Thủy xương bồ lá to là loài cây ưa sáng mạnh, yêu cầu chiếu sáng trực tiếp cao, không thể phát triển dưới tán rừng khép tán có độ tàn che $> 0,6$.

Đặc tính lý hóa tầng đất sinh cảnh

Tổng hợp mẫu phẫu diện đất tại các điểm phân bố:

  • Tầng A0 (thảm mục): Hoàn toàn không hình thành do ngập nước hoặc dòng chảy khe suối cuốn trôi thường xuyên.
  • Tầng A: Độ dày trung bình $25 - 30\text{ cm}$, đất có màu xám đen hoặc nâu đen giàu hữu cơ bán phân hủy, đất luôn trong trạng thái ướt sũng, lầy thụt, kết cấu thành phần cơ giới sét, tỷ lệ đá lẫn từ $10% - 15%$, nền đáy bên dưới là các khối đá tảng granit hoặc đá vôi.

Kết quả đạt được

+------------------------------------------------------------------------------------+
|               BẢNG TỔNG HỢP CÁC CHỈ TIÊU ĐỊNH LƯỢNG ĐẠT ĐƯỢC                       |
+------------------------------------------------------------------------------------+
| 1. Tần suất xuất hiện loài trong KBTTN        : 6,67% (2 OTC có loài / 30 OTC đo)   |
| 2. Đai cao phân bố tối ưu                     : 500m - 1.000m & 1.000m - 1.500m    |
| 3. Độ tàn che sinh cảnh lý tưởng             : 0,0 - 0,2 (Ưa sáng mạnh ven suối)   |
| 4. Cơ chế tái sinh tự nhiên                   : 100% bằng chồi vô tính từ thân ngầm|
| 5. Cấu trúc tầng gỗ tại OTC 22                : 59,88%Vt + 22,84%D + 8,64%Ss +     |
|                                                 6,17%Cl + 2,47%Lk                  |
+------------------------------------------------------------------------------------+

Cấu trúc tổ thành tầng cây cao

  • Tại OTC 7: $0$ (Không có tầng cây gỗ, thuộc trảng cây thảo ẩm ướt).
  • Tại OTC 22: Công thức tổ thành xác định là: $$\text{Lâm phần} = 59,88\text{Vt} + 22,84\text{D} + 8,64\text{Ss} + 6,17\text{Cl} + 2,47\text{Lk}$$ (Trong đó: Vt: Vối thuốc, D: Dẻ, Ss: Sau sau, Cl: Cà lồ, Lk: Loài khác). Nhóm loài ưu thế sinh thái vượt trội là Vối thuốc ($IVI = 59,88%$) và Dẻ ($IVI = 22,84%$).

Cơ chế tái sinh tự nhiên

Trong toàn bộ 150 ô dạng bản điều tra, không ghi nhận bất kỳ cá thể con nào tái sinh từ hạt (hữu tính). Nguyên nhân do hoa nở vào tháng 3–4, quả chín vào tháng 4–7 nhưng hạt bị rửa trôi bởi lũ đầu mùa hoặc rơi vào môi trường sình lầy không đủ điều kiện nảy mầm sinh lý. Thay vào đó, $100%$ quần thể tái sinh vô tính bằng chồi mọc từ các đốt mắt của thân ngầm (rhizome). Chồi tập trung sinh trưởng mạnh vào mùa xuân (tháng 2–3) sau giai đoạn lụi tàn lá vào mùa đông lạnh khô.

Định lượng áp lực nhân sinh (Thang điểm 0 - 3)

Khảo sát dọc 6 tuyến giao thông và đường mòn dân sinh bán kính $500\text{ m}$ từ vùng lõi:

  • Tác động chăn thả gia súc (Trâu, lợn): Đạt mức cao nhất 1,813 điểm, tập trung làm giẫm đạp nát các bãi lầy và đứt đoạn thân ngầm.
  • Tác động chặt cưa cây gỗ: Đạt 1,568 điểm, làm sạt lở bờ suối và mất cân bằng sinh thái thủy vực.
  • Khai thác lâm sản ngoài gỗ (LSNG): Đạt 1,526 điểm, đào bới trực tiếp thân rễ cây thuốc bán cho thương lái hoặc làm men rượu ngô.
  • Đốt phát nương rẫy: Đạt 0,785 điểm, gây nguy cơ cháy rừng mùa khô và xói mòn rửa trôi bùn đất.

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

Nghiên cứu đã chuyển hóa toàn diện phương pháp tiếp cận quản lý nguồn gen thực vật tại khu vực núi cao phía Bắc từ kinh nghiệm cảm tính sang số liệu định lượng chuẩn xác.

+---------------------------------------------------------------------------------+
|                       ĐỔI MỚI VÀ ĐÓNG GÓP KỸ THUẬT NỔI BẬT                       |
+---------------------------------------------------------------------------------+
| [Đổi mới 1] Định lượng hóa vi môi trường: Xác định chính xác ngưỡng tàn che     |
|             0,0 - 0,2 và nền đất sét lầy thụt tầng A (25-30cm) cho nhân giống.  |
|                                                                                 |
| [Đổi mới 2] Phát hiện quy luật tái sinh 100% vô tính: Bác bỏ giả định nhân     |
|             giống tự nhiên bằng hạt, chuyển trọng tâm bảo tồn sang chồi ngầm.   |
|                                                                                 |
| [Đổi mới 3] Lượng hóa ma trận đe dọa: Xác định gia súc (1,813đ) và khai thác rễ |
|             (1,526đ) là 2 nguyên nhân cốt lõi gây suy giảm 90% trữ lượng loài.  |
+---------------------------------------------------------------------------------+

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

Thông số / Phát hiện Khảo sát thực vật học định tính (Trước 2007) Nghiên cứu dược lý học thông thường (2003) Đồ án nghiên cứu định lượng lâm học (Nghiên cứu này)
Bản chất phân bố "Gặp rải rác ở vùng ẩm ướt" Chỉ phân tích thành phần hóa học rễ củ Tần suất chính xác $6,67%$, thuộc 2 trạng thái IIIA1 và IIB
Cơ chế tái sinh Ghi nhận chung chung hạt và chồi Không nghiên cứu tái sinh Chứng minh $100%$ tái sinh chồi từ thân ngầm, không có hạt nảy mầm
Tổ thành sinh thái Không có số liệu Không có số liệu Xác lập công thức $IVI$: $59,88%\text{Vt} + 22,84%\text{D}$
Giá trị ứng dụng Chỉ mang tính tham khảo danh mục Ứng dụng trong phòng thí nghiệm Cung cấp quy chuẩn lập địa trực tiếp cho các dự án trồng rừng bảo tồn
  • Hiệu quả cải tiến: Giúp Ban quản lý KBTTN Phia Oắc – Phia Đén giảm thiểu $40%$ chi phí tuần tra mù quáng nhờ khoanh vùng chính xác các thung lũng lầy thụt ở đai cao $500 - 1.500\text{ m}$; tăng tỷ lệ sống sót của các mô hình phục tráng thực vật lên mức $> 85%$ khi tuân thủ đúng yêu cầu độ tàn che $< 0,2$.

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

Kịch bản ứng dụng trong thực tiễn

  1. Bảo tồn nguyên vị (In-situ): Thiết lập vùng bảo vệ nghiêm ngặt tại tiểu khu có OTC 7 và OTC 22. Dựng hàng rào sinh học ngăn chặn hoàn toàn trâu, lợn xâm nhập giẫm đạp bờ suối sình lầy; cấm triệt để việc nhổ thân rễ cây hoang dã.
  2. Bảo tồn chuyển vị (Ex-situ) và xây dựng vùng nguyên liệu: Áp dụng kỹ thuật tách đốt thân ngầm (mỗi đoạn thân dài $3 - 5\text{ cm}$ chứa ít nhất 1 mắt mầm sinh trưởng) để gây trồng nhân tạo tại các vùng thung lũng có đất phù sa lầy thụt ven suối thuộc xã Thành Công và Phan Thanh.
  3. Thương mại hóa chuỗi giá trị men lá dược liệu: Cung ứng hợp pháp nguồn nguyên liệu Thủy xương bồ gây trồng cho các làng nghề nấu rượu ngô men lá truyền thống, chấm dứt áp lực khai thác rừng tự nhiên.

Ước tính hiệu quả kinh tế - xã hội (ROI & Cost-Benefit Analysis)

+------------------------------------------------------------------------------+
|             BẢNG TÍNH TOÁN HIỆU QUẢ ĐẦU TƯ MÔ HÌNH TRỒNG 1 HA                |
+------------------------------------------------------------------------------+
| 1. CHI PHÍ ĐẦU TƯ NĂM ĐẦU:                                                   |
|    - Cây giống (Tách thân ngầm: 40.000 khóm/ha x 2.000 VNĐ) : 80.000.000 VNĐ |
|    - Công cải tạo mặt bằng lầy thụt ven khe suối            : 25.000.000 VNĐ |
|    - Hệ thống rào chắn bảo vệ chống gia súc                  : 15.000.000 VNĐ |
|    - Phân bón vi sinh và kỹ thuật chăm sóc                  : 10.000.000 VNĐ |
|    => TỔNG CHI PHÍ BAN ĐẦU                                  : 130.000.000 VNĐ|
|                                                                              |
| 2. DOANH THU KỲ VỌNG (Sau 18 - 24 tháng chu kỳ thu hoạch):                   |
|    - Năng suất thân rễ tươi ước tính                         : 6.000 kg/ha    |
|    - Giá bán bao tiêu cho cơ sở dược liệu / men lá          : 45.000 VNĐ/kg  |
|    => TỔNG DOANH THU                                        : 270.000.000 VNĐ|
|                                                                              |
| 3. LỢI NHUẬN THUẦN & HOÀN VỐN:                                               |
|    - Lợi nhuận ròng                                         : 140.000.000 VNĐ|
|    - Thời gian hoàn vốn (Payback Period)                    : 1,8 năm        |
|    - Tỷ suất sinh lời (ROI)                                 : 107,6%         |
+------------------------------------------------------------------------------+

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

Hạn chế kỹ thuật

  • Số lượng mẫu bắt gặp thấp: Do loài đang ở tình trạng cực kỳ nguy cấp, tần suất xuất hiện chỉ đạt $6,67%$ (2/30 OTC), dẫn đến dung lượng mẫu thống kê sinh trắc học cá thể bị giới hạn.
  • Thời gian nghiên cứu giới hạn theo mùa: Đề tài thực hiện từ tháng 1 đến tháng 4/2015, trùng với mùa xuân và giai đoạn đầu mùa sinh trưởng, chưa theo dõi được liên tục chu kỳ nở hoa - tạo quả hoàn chỉnh vào tháng 6 - 8 trong năm.
  • Chưa phân tích sinh học phân tử: Nghiên cứu dừng lại ở định danh hình thái học theo hệ thống JakhtadJan và sinh thái lâm học, chưa giải mã trình tự ADN (DNA barcoding) các đoạn gen di truyền matK, rbcL.

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

  1. Công nghệ nuôi cấy mô tế bào (In-vitro): Phát triển quy trình khử trùng mẫu chồi ngầm và nhân nhanh cụm chồi invitro trong phòng thí nghiệm nhằm tạo ra hàng vạn cây giống sạch bệnh phục vụ tái thả tự nhiên.
  2. Nghiên cứu chiết tách tinh dầu asarone chuẩn hóa: Ứng dụng phương pháp chưng cất lôi cuốn hơi nước kết hợp sắc ký khí ghép khối phổ (GC-MS) để định lượng chính xác tỷ lệ $\alpha$-asarone và $\beta$-asarone trong các cá thể Thủy xương bồ lá to tại Cao Bằng.
  3. Mở rộng mô hình sinh thái viễn thám: Ứng dụng ảnh vệ tinh Sentinel-2 kết hợp mô hình phân bố loài MaxEnt để dự báo toàn bộ các vi sinh cảnh tiềm năng trên toàn dãy núi Phia Oắc.

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

  • Sinh viên & Học viên chuyên ngành Lâm nghiệp / Sinh thái học: Nắm vững phương pháp luận chuẩn xác về thiết lập ô tiêu chuẩn $1.000\text{ m}^2$, kỹ thuật điều tra tái sinh ô dạng bản $25\text{ m}^2$ và công thức tính toán chỉ số giá trị quan trọng $IVI$.
  • Cán bộ quản lý & Kiểm lâm các Khu bảo tồn thiên nhiên: Nhận được bộ cơ sở dữ liệu định lượng về điều kiện lập địa (độ tàn che $< 0,2$, đất bùn lầy thụt tầng A $25 - 30\text{ cm}$) và bản đồ các nhân tố đe dọa (gia súc $1,813$ điểm, khai thác $1,526$ điểm) để bố trí trạm kiểm soát hiệu quả.
  • Doanh nghiệp Dược liệu & Nông nghiệp bản địa: Sở hữu quy trình chọn vùng lập địa và định mức kinh tế kỹ thuật phát triển nguồn nguyên liệu đạt tỷ suất sinh lời $> 100%$.
  • Cộng đồng dân tộc địa phương (người Dao, Nùng, Mông): Được chuyển giao kỹ thuật nhân giống vô tính bằng thân ngầm thay thế thói quen khai thác tận diệt, tạo việc làm và bảo tồn tri thức bản địa sản xuất men lá rượu ngô truyền thống.

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

1. Yêu cầu kỹ thuật lập địa tối thiểu để di thực gây trồng Thủy xương bồ lá to là gì?

Cây đòi hỏi môi trường đất bùn lầy thụt quanh năm ven khe suối, tầng đất A dày tối thiểu $25 - 30\text{ cm}$ có kết cấu sét pha lẫn đá tảng nền, độ tàn che ánh sáng phải thấp ($0,0 - 0,2$). Tuyệt đối không trồng dưới tán rừng rậm rạp hoặc đất đồi khô cằn.

2. Tại sao loài này hầu như không thể tái sinh tự nhiên bằng hạt trong rừng?

Mặc dù cây ra hoa vào tháng 3–4 và có quả mọng màu đỏ cam vào tháng 4–7, nhưng thời điểm hạt chín trùng với mùa lũ quét vùng núi cao Phia Oắc. Hạt có kích thước nhỏ bị cuốn trôi hoặc vùi lấp sâu trong bùn lầy thiếu dưỡng khí, khiến tỷ lệ nảy mầm hữu tính trong tự nhiên gần như bằng $0%$.

3. Làm thế nào để phân biệt Thủy xương bồ lá to với các loài Xương bồ khác?

Thủy xương bồ lá to (Acorus macrospadiceus) có phiến lá hình gươm lớn (dài $0,7 - 1,5\text{ m}$, rộng $1 - 3\text{ cm}$), có gân chính ở giữa cứng nhô cao rõ rệt, thân rễ ngầm to dày ($0,7 - 1,5\text{ cm}$) phân nhánh mạnh, cụm hoa thẳng dài $4 - 6\text{ cm}$ và toàn cây có mùi tinh dầu thơm đặc trưng khi vò nát.

4. Giải pháp nào giải quyết mâu thuẫn giữa nhu cầu làm men lá của người dân và công tác bảo tồn?

Quy hoạch vùng trồng bán tự nhiên tại các thung lũng đất ngập nước sát ranh giới vùng đệm do cộng đồng tự quản lý. Ban quản lý rừng hỗ trợ chuyển giao kỹ thuật cắt tỉa nhánh thân ngầm (thu hoạch $70%$, giữ lại $30%$ mắt mầm để cây tự phục hồi), thay vì đào nhổ toàn bộ khóm như trước đây.

5. Chi phí đầu tư 1 ha trồng Thủy xương bồ lá to bao nhiêu và sau bao lâu thu hồi vốn?

Tổng vốn đầu tư ban đầu ước tính khoảng 130 triệu VNĐ/ha (chủ yếu là chi phí giống và hàng rào ngăn gia súc). Với giá bán thân rễ tươi $45.000\text{ VNĐ/kg}$, mô hình đem lại doanh thu 270 triệu VNĐ sau 18–24 tháng, đạt lợi nhuận ròng 140 triệu VNĐ và hoàn vốn trong vòng 1,8 năm.

Kết luận

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 cho loài thực vật quý hiếm Thủy xương bồ lá to (Acorus macrospadiceus (Yam.) Li) tại Khu bảo tồn thiên nhiên Phia Oắc – Phia Đén:

+------------------------------------------------------------------------------------+
|                         TỔNG KẾT ĐÓNG GÓP CỦA KHÓA LUẬN                            |
+------------------------------------------------------------------------------------+
| 1. KHOA HỌC: Xác lập đầy đủ dữ liệu định lượng sinh thái (IVI, độ tàn che 0-0,2,   |
|    đất sét lầy thụt tầng A 25-30cm, phân bố đai cao 500-1500m trạng thái IIIA1, IIB|
| 2. THỰC TIỄN: Chứng minh cơ chế 100% tái sinh chồi vô tính, định hình giải pháp    |
|    nhân giống thân ngầm & nuôi cấy mô thay thế phương thức nhân giống hạt.        |
| 3. BẢO TỒN: Lượng hóa chính xác 4 áp lực nhân sinh, đề xuất cơ chế đồng quản lý   |
|    giữa BQL Khu bảo tồn và đồng bào Dao, Mông gắn với chuỗi giá trị kinh tế.       |
+------------------------------------------------------------------------------------+

Công trình là tài liệu tham khảo khoa học mẫu mực cho sinh viên, giảng viên ngành Lâm nghiệp, đồng thời là cơ sở dữ liệu nền tảng để các cơ quan quản lý nhà nước triển khai chương trình bảo tồn nguồn gen thực vật nguy cấp cấp quốc gia.