Giới thiệu dự án

Bối cảnh và thực trạng nguồn gen thực vật rừng

Theo báo cáo của Tổ chức Lương thực và Nông nghiệp Liên Hợp Quốc (FAO) và Liên minh Bảo tồn Thiên nhiên Quốc tế (IUCN), khoảng 13% đến 47% số loài thực vật trên toàn cầu đang đối mặt với nguy cơ tuyệt chủng do suy thoái sinh cảnh, biến đổi khí hậu và áp lực khai thác quá mức. Tại Việt Nam, áp lực từ việc chuyển đổi đất rừng và khai thác lâm sản đặc hữu đã làm suy giảm nghiêm trọng các quần thể cây lá kim quý hiếm thuộc nhóm Hạt trần (Gymnospermae).

Sa mộc dầu (Cunninghamia konishii Hayata) — còn gọi là Sa mộc quế phong, Sa mu dầu, hoặc Ngọc am thuộc họ Hoàng đàn (Cupressaceae) — là loài thực vật gỗ lớn đặc hữu, có giá trị khoa học và kinh tế cao. Loài được xếp vào nhóm nguy cấp trong Sách Đỏ Việt Nam (2007) ở phân hạng VU A1adC1, nằm trong Nhóm IIA của Nghị định số 32/2006/NĐ-CP, và thuộc Danh lục Đỏ IUCN cấp độ VU A1cd.

+-------------------------------------------------------------------------+
|                    QUY CÁCH PHÂN HẠNG BẢO TỒN VÀ NGUY CƠ                |
+------------------------------------+------------------------------------+
| Danh mục quản lý / Đánh giá        | Phân hạng / Tình trạng pháp lý     |
+------------------------------------+------------------------------------+
| Sách Đỏ Việt Nam (2007)            | Sẽ nguy cấp (VU A1adC1)            |
| Nghị định 32/2006/NĐ-CP             | Nhóm IIA (Hạn chế khai thác TM)    |
| Danh lục Đỏ IUCN (2007)            | Sẽ nguy cấp (VU A1cd)              |
| Công ước CITES (2008)              | Phụ lục II                         |
+------------------------------------+------------------------------------+

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

Tại Khu bảo tồn thiên nhiên (KBTTN) Tây Côn Lĩnh (tỉnh Hà Giang), với tổng diện tích 15.000 ha (trong đó phân khu bảo vệ nghiêm ngặt là 7.406 ha và phân khu phục hồi sinh thái là 7.173 ha), Sa mộc dầu đang đứng trước nguy cơ tuyệt chủng cục bộ do các nguyên nhân:

  • Áp lực khai thác triệt để: Gỗ Sa mộc dầu có vân thớ mịn, mùi tinh dầu thơm đặc trưng, chịu mối mọt tốt, bị khai thác làm quan quách, đồ mỹ nghệ cao cấp, bồn tắm và chưng cất tinh dầu thương phẩm. Người dân địa phương thậm chí đào bới cả gốc rễ để bán qua biên giới.
  • Tái sinh tự nhiên kém: Tán rừng tự nhiên rậm rạp ngăn cản ánh sáng đến tầng cỏ sạ, trong khi nón cái khi chín không tự tách hạt trên cây mà rụng cả nón xuống đất, khiến hạt khó phát tán xa bằng gió và tỷ lệ nảy mầm hình thành cây con gần như bị triệt tiêu dưới tán rừng nguyên sinh.
  • Thiếu hụt dữ liệu định lượng: Thiếu các nghiên cứu định lượng về phân bố không gian, toạ độ định vị, cấu trúc tầng thứ, hệ số tổ thành và tổ chức tái sinh tại các tiểu khu của Tây Côn Lĩnh.
                   SƠ ĐỒ NGUYÊN NHÂN SUY GIẢM QUẦN THỂ SA MỘC DẦU
                   
+---------------------------+       +---------------------------+
|  Khai thác gỗ & tinh dầu  |       |   Đặc tính nón không nở   |
|   (Thương phẩm, đào rễ)   |       |   (Hạt kẹt trong nón rơi)  |
+-------------+-------------+       +-------------+-------------+
              |                                   |
              +-----------------+-----------------+
                                |
                                v
               [ SUY GIẢM MẬT ĐỘ QUẦN THỂ TRẦM TRỌNG ]
               [ 11 CÁ THỂ / 9 Ô TIÊU CHUẨN ĐIỀU TRA ]
                                |
              +-----------------+-----------------+
              |                                   |
+-------------v-------------+       +-------------v-------------+
|  Thiếu cây con tái sinh   |       |  Suy thoái đa dạng di     |
|   (Tầng râm, thảm mục dày)|       |  truyền (Xói mòn nguồn gen)|
+---------------------------+       +---------------------------+

Mục tiêu của dự án

  1. Xác định chính xác vùng phân bố địa lý, độ cao, tọa độ GPS và diện tích sinh cảnh của loài Sa mộc dầu tại KBTTN Tây Côn Lĩnh.
  2. Định lượng các chỉ số lâm học gồm: cấu trúc mật độ ($N/ha$), đường kính ngang ngực ($D_{1.3}$), chiều cao vút ngọn ($H_{vn}$), cấu trúc tầng thứ, và xác lập công thức tổ thành loài ($N_a$).
  3. Khảo sát hình thái học, vật hậu học (chu kỳ ra nón, chín hạt) và cơ chế tái sinh tự nhiên dưới tán rừng.
  4. Đề xuất hệ thống giải pháp tích hợp: kỹ thuật lâm sinh (xúc tiến tái sinh, gieo ươm), chính sách khoanh nuôi đồng quản lý và phát triển sinh kế vùng đệm.

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

  • Thời gian thực hiện: Từ 31/12/2013 đến 20/04/2014.
  • Địa bàn: 3 xã Túng Sán (huyện Hoàng Su Phì), Tả Sử Choóng và Bản Nhùng (huyện Hoàng Su Phì, KBTTN Tây Côn Lĩnh), tập trung ở đai cao 900 m – 1.200 m thuộc các tiểu khu 194A, 194B và 204.
  • Giới hạn: Do thời gian điều tra không trùng với mùa cây trổ nón hoa (tháng 9-10), nghiên cứu tập trung vào nón quả chín, hình thái cành lá, vỏ và đo đạc sinh quyển thực địa.

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

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

+---------------------------------------------------------------------------------------------------+
|                           SO SÁNH CÁC PHƯƠNG THỨC BẢO TỒN NGUỒN GEN RỪNG                          |
+----------------------+--------------------------------+-------------------------------------------+
| Tiêu chí so sánh     | Phương thức 1: In-situ thụ động | Phương thức 2: Ex-situ thuần túy          |
|                      | (Đóng cửa rừng truyền thống)   | (Ngân hàng hạt giống / Vườn thực vật)     |
+----------------------+--------------------------------+-------------------------------------------+
| Cơ chế tác động      | Giữ nguyên hiện trạng rừng,    | Thu thập nguồn gen đưa ra ngoài           |
|                      | cấm tuyệt đối người dân vào    | sinh cảnh tự nhiên (như mô hình Đài Loan) |
+----------------------+--------------------------------+-------------------------------------------+
| Ưu điểm              | Không xáo trộn sinh thái,      | Kiểm soát tốt nảy mầm nhân tạo,           |
|                      | chi phí đầu tư ban đầu thấp    | phòng ngừa rủi ro mất trắng ngoài tự nhiên|
+----------------------+--------------------------------+-------------------------------------------+
| Nhược điểm           | Cây con không thể tự tái sinh   | Đòi hỏi công nghệ lưu giữ lạnh,           |
|                      | do thiếu sáng; lâm tặc lén lút | nguy cơ thoái hóa giống khi tái nhập nội  |
+----------------------+--------------------------------+-------------------------------------------+
| Hiệu quả bền vững    | Thấp (Quần thể già cỗi, suy tàn)| Trung bình (Cần nguồn vốn duy trì lớn)    |
+----------------------+--------------------------------+-------------------------------------------+

+---------------------------------------------------------------------------------------------------+
|                PHƯƠNG THỨC ĐỀ XUẤT: BẢO TỒN IN-SITU CHỦ ĐỘNG KẾT HỢP LÂM SINH XÃ HỘI              |
+---------------------------------------------------------------------------------------------------+
| Cơ chế tác động: Điều tra GPS xác lập ranh giới tiểu khu + Mở tán rừng cục bộ (Light Gaps)        |
|                  + Xử lý nón chín tại chỗ + Trao quyền tuần tra đồng quản lý cho người Dao/H'Mông |
| Ưu điểm:         Khắc phục triệt để điểm nghẽn tái sinh; gắn lợi ích tuần rừng với sinh kế bản địa |
| Nhược điểm:      Đòi hỏi sự phối hợp liên ngành giữa Kiểm lâm, Nhà khoa học và Cộng đồng thôn bản |
| Hiệu quả:        Tối ưu nhất cho các quần thể cực nhỏ (< 50 cá thể) có nguy cơ tuyệt chủng cao    |
+---------------------------------------------------------------------------------------------------+

Ma trận yêu cầu giải pháp theo mô hình MoSCoW

  • Must Have (Bắt buộc có):
    • Tọa độ GPS độ chính xác cao cho từng cá thể Sa mộc dầu còn sót lại.
    • Công thức tổ thành tầng cây cao ($N_a$) tại các ô tiêu chuẩn $1.000\ m^2$.
    • Biện pháp kỹ thuật thu gom nón chín (tháng 3-5) để tách hạt nhân tạo.
  • Should Have (Nên có):
    • Bản đồ phân vùng tiểu khu sinh cảnh (194A, 194B, 204) trên nền tảng QGIS.
    • Cơ chế chi trả kinh phí tuần tra rừng thường xuyên cho người dân bản địa am hiểu địa hình.
  • Could Have (Có thể có):
    • Chương trình thử nghiệm giâm hom và nhân giống mô vô tính tại vườn ươm lâm nghiệp.
    • Mô hình trồng rừng hỗn giao Sa mộc dầu với các loài cây cộng sinh bản địa (Sa mộc, Vối thuốc).
  • Won't Have (Chưa thực hiện ở giai đoạn này):
    • Khai thác thương mại tinh dầu quy mô công nghiệp từ các quần thể tự nhiên.

Thiết kế hệ thống điều tra và quản lý dữ liệu lâm học

Hệ thống điều tra thực địa được tổ chức thành cấu trúc 3 tầng phân cấp: Tuyến khảo sát (Transects) $\rightarrow$ Ô tiêu chuẩn tạm thời (Standard Plots - OTC) $\rightarrow$ Cây cá thể mục tiêu (Individual Trees).

       MÔ HÌNH CẤU TRÚC DỮ LIỆU ĐIỀU TRA LÂM HỌC VÀ ĐA DẠNG SINH HỌC
       
  +--------------------------------------------------------------------+
  |                        KBTTN TÂY CÔN LĨNH                          |
  |  Diện tích: 15.000 ha | Tọa độ: 22°39'-22°50' B, 104°39'-104°59' Đ |
  +---------------------------------+----------------------------------+
                                    |
            +-----------------------+-----------------------+
            |                       |                       |
            v                       v                       v
    [ Tuyến 1: Túng Sán ]  [ Tuyến 2: Tả Sử Choóng ]  [ Tuyến 3: Bản Nhùng ]
    (Độ cao 900-1200m)     (Độ cao 900-1200m)        (Độ cao 900-1200m)
            |                       |                       |
      +-----+-----+           +-----+-----+           +-----+-----+
      |     |     |           |     |     |           |     |     |
      v     v     v           v     v     v           v     v     v
    OTC 1 OTC 2 OTC 3       OTC 4 OTC 5 OTC 6       OTC 7 OTC 8 OTC 9
    (1000m2/OTC)            (1000m2/OTC)            (1000m2/OTC)
      |     |     |                 |                 |
      -   2 cây 2 cây               5 cây             2 cây
            |     |                 |                 |
            +-----+-----------------+-----------------+
                                    |
                                    v
                     [ TỔNG HỢP: 11 CÁ THỂ SA MỘC DẦU ]

Ngăn xếp công nghệ và công cụ lâm sinh

  • Thiết bị định vị & trắc địa: Máy định vị vệ tinh GPS Garmin GPSMAP 64s, Địa bàn địa chất cầm tay Harbin DQL-8, Thước dây thép 50m, Thước đo đường kính cây Richter ($D_{1.3}$).
  • Nền tảng phân tích không gian & thống kê:
    • Python 3.10+ (Thư viện pandas v2.2.0, geopandas v0.14.2, scipy v1.12.0) để xử lý hệ số tổ thành và mô hình hóa phân bố.
    • Phần mềm QGIS 3.34 LTR xử lý bản đồ không gian và phân khu bảo tồn.
  • Tiêu chuẩn thiết kế cơ sở dữ liệu lâm học:
-- Schema cấu trúc lưu trữ cơ sở dữ liệu điều tra Sa mộc dầu
CREATE TABLE plot_survey (
    plot_id VARCHAR(10) PRIMARY KEY,
    transect_name VARCHAR(50) NOT NULL,
    sub_compartment VARCHAR(20) NOT NULL,
    elevation_m INT CHECK (elevation_m BETWEEN 500 AND 2500),
    slope_percent INT CHECK (slope_percent BETWEEN 0 AND 100),
    latitude_utm VARCHAR(30) NOT NULL,
    longitude_utm VARCHAR(30) NOT NULL,
    survey_date DATE NOT NULL
);

CREATE TABLE tree_mensuration (
    tree_id SERIAL PRIMARY KEY,
    plot_id VARCHAR(10) REFERENCES plot_survey(plot_id),
    species_scientific_name VARCHAR(100) DEFAULT 'Cunninghamia konishii Hayata',
    species_local_name VARCHAR(50) DEFAULT 'Sa moc dau',
    dbh_cm NUMERIC(5, 2) NOT NULL, -- D1.3 (cm)
    height_m NUMERIC(5, 2) NOT NULL, -- Hvn (m)
    vitality_status VARCHAR(30) CHECK (vitality_status IN ('Tot', 'Trung binh', 'Xau')),
    regeneration_count INT DEFAULT 0
);

Phương pháp nghiên cứu và công thức tính toán

Hệ thống công thức toán học lâm sinh áp dụng

  1. Mật độ cây tái sinh ($N/ha$): $$N/ha = \frac{10.000 \times n}{S}$$ Trong đó: $S$ là tổng diện tích các ô điều tra tái sinh ($m^2$), $n$ là tổng số cây tái sinh ghi nhận.

  2. Tỷ lệ phân cấp chất lượng cây tái sinh ($N%$): $$N% = \left( \frac{n_i}{N_{total}} \right) \times 100$$ Trong đó: $n_i$ là số lượng cây tái sinh ở cấp phẩm chất $i$ (Tốt, Trung bình, Xấu), $N_{total}$ là tổng số cây tái sinh.

  3. Ngưỡng số cây trung bình loài ưu thế ($N_{TB}$): $$N_{TB} = \frac{N}{m}$$ Trong đó: $N$ là tổng số cây cá thể tầng cao trong OTC, $m$ là tổng số loài cây xuất hiện trong OTC. Cây chiếm vị trí ưu thế khi số cá thể $n \ge N_{TB}$.

  4. Hệ số tổ thành loài ($N_a$): $$N_a = \left( \frac{n_a}{N} \right) \times 10$$ Trong đó: $n_a$ là số lượng cá thể của loài $a$, $N$ là tổng số cá thể các loài trong OTC.

+-------------------------------------------------------------------------+
|                  QUY TẮC XÂY DỰNG CÔNG THỨC TỔ THÀNH LÂM HỌC             |
+-----------------------------+-------------------------------------------+
| Hệ số tổ thành ($N_a$)      | Cách thức ghi nhận trong biểu thức       |
+-----------------------------+-------------------------------------------+
| $N_a \ge 1.0$               | Ghi trực tiếp giá trị số + Mã loài        |
| $0.5 \le N_a < 1.0$         | Thêm dấu (+) trước mã loài                 |
| $N_a < 0.5$                 | Thêm dấu (-) trước mã loài                 |
| Các loài dưới ngưỡng TB     | Gom chung vào nhóm loài khác (LK)         |
+-----------------------------+-------------------------------------------+

           TIẾN TRÌNH KHẢO SÁT VÀ PHÂN TÍCH THEO THỜI GIAN (2013-2014)
           
 31/12/2013 - 15/01/2014        16/01/2014 - 15/03/2014       16/03/2014 - 20/04/2014
+-----------------------+      +-----------------------+      +-----------------------+
|  Giai đoạn 1: Chuẩn bị|      |  Giai đoạn 2: Ngoại   |      |  Giai đoạn 3: Phân tích|
|  - Kế thừa tài liệu   | ===> |  nghiệp thực địa      | ===> |  - Tính toán tổ thành |
|  - Thiết kế 3 tuyến   |      |  - Lập 9 OTC (1000m2) |      |  - Định vị tiểu khu   |
|  - Soạn 20 phiếu PV   |      |  - Đo D1.3, Hvn, GPS  |      |  - Hoàn thiện báo cáo |
+-----------------------+      +-----------------------+      +-----------------------+

Ma trận đánh giá và giảm thiểu rủi ro ngoại nghiệp

  • Rủi ro địa hình hiểm trở, vách đá chia cắt: Trang bị bảo hộ leo núi chuyên dụng, định tuyến song song đường đồng mức, có kiểm lâm địa bàn đi cùng hỗ trợ mở lối.
  • Rủi ro rào cản ngôn ngữ với cộng đồng bản xứ: Lựa chọn phỏng vấn có sự tham gia của già làng, trưởng bản và người dẫn đường là người H'Mông, người Dao địa phương.
  • Rủi ro sai số định vị GPS dưới tán rừng kín: Sử dụng máy định vị đa băng tần kết hợp đối chiếu mốc địa hình tự nhiên và bản đồ địa hình tỷ lệ 1:10.000.

Implementation và kết quả

Quy trình điều tra ngoại nghiệp và mã nguồn tính toán

Ngoại nghiệp được triển khai với 9 ô tiêu chuẩn tạm thời (diện tích mỗi OTC là $1.000\ m^2$, kích thước tiêu chuẩn $30\ m \times 33,33\ m$ hoặc $20\ m \times 50\ m$ tùy địa hình) tại 3 đai cao: Đai 1 (900-1000m), Đai 2 (1000-1100m), và Đai 3 (1100-1200m).

Thuật toán Python phân tích dữ liệu điều tra lâm học

Dưới đây là đoạn mã Python thực hiện tự động hóa việc tính toán hệ số tổ thành và xuất công thức lâm phần theo chuẩn lâm nghiệp:

import pandas as pd
from typing import Dict, List, Tuple

def compute_forest_composition(tree_data: Dict[str, int]) -> Tuple[str, pd.DataFrame]:
    """
    Tính toán hệ số tổ thành lâm học theo phương pháp 10 phần mười.
    
    Args:
        tree_data: Dictionary dạng {'Tên loài/Mã': Số lượng cá thể}
    Returns:
        Công thức tổ thành và DataFrame chi tiết hệ số.
    """
    df = pd.DataFrame(list(tree_data.items()), columns=['Species', 'Count'])
    total_trees = df['Count'].sum()
    num_species = len(df)
    n_mean = total_trees / num_species
    
    # Tính hệ số tổ thành theo thang 10
    df['Composition_Score'] = (df['Count'] / total_trees) * 10
    df['Is_Dominant'] = df['Count'] >= n_mean
    
    # Định dạng công thức tổ thành
    formula_parts = []
    for _, row in df.sort_values(by='Composition_Score', ascending=False).iterrows():
        score = round(row['Composition_Score'], 2)
        code = row['Species']
        if score >= 1.0:
            formula_parts.append(f"{score}{code}")
        elif 0.5 <= score < 1.0:
            formula_parts.append(f"+{score}{code}")
        else:
            formula_parts.append(f"-{score}{code}")
            
    composition_formula = " + ".join(formula_parts)
    return composition_formula, df

# Thực thi mẫu cho OTC 5 (Tuyến Tả Sử Choóng, đai 1000-1100m)
sample_otc5 = {
    'Sm': 21,   # Sa mộc (Cunninghamia lanceolata)
    'Lk': 7,    # Loài khác
    'Smd': 5,   # Sa mộc dầu (Cunninghamia konishii)
    'Mlt': 5,   # Mò lá tròn
    'Qr': 4,    # Quế rừng
    'Cn': 4     # Cọ nọt
}

formula, summary_df = compute_forest_composition(sample_otc5)
print(f"Công thức tổ thành OTC 5: {formula}")

Kiểm định và đối soát dữ liệu thực địa

Dữ liệu được kiểm chứng chéo qua 20 phiếu phỏng vấn bán cấu trúc với 4 nhóm đối tượng: Cán bộ Ban quản lý KBTTN Tây Côn Lĩnh, Kiểm lâm huyện Vị Xuyên và Hoàng Su Phì, cơ sở chưng cất tinh dầu, và các hộ dân sinh sống tại vùng đệm.

+-------------------------------------------------------------------------+
|                   KẾT QUẢ ĐỐI SOÁT ĐỘ TIN CẬY NGOẠI NGHIỆP              |
+-------------------------------+---------------------+-------------------+
| Hạng mục kiểm định            | Giá trị ghi nhận    | Ngưỡng dung sai   |
+-------------------------------+---------------------+-------------------+
| Độ lệch chuẩn vị trí GPS      | $\pm 2.8\ m$         | $< 5.0\ m$        |
| Sai số đo đường kính $D_{1.3}$| $\pm 0.5\ cm$       | $< 1.0\ cm$       |
| Sai số đo chiều cao $H_{vn}$  | $\pm 0.6\ m$        | $< 1.0\ m$        |
| Tỷ lệ hồi đáp phiếu phỏng vấn | 100% (20/20 phiếu)  | $\ge 90\%$        |
+-------------------------------+---------------------+-------------------+

Kết quả nghiên cứu đạt được

1. Dữ liệu phân bố không gian và trắc lượng cây cá thể

Trong toàn bộ 9 ô tiêu chuẩn khảo sát, chỉ phát hiện 11 cá thể Sa mộc dầu sống sót, phân bố rải rác ở 4 ô tiêu chuẩn (OTC 2, OTC 3, OTC 5, OTC 7) tại 3 xã thuộc 3 tiểu khu với tổng diện tích sinh cảnh ước tính 2.888,5 ha.

+-------------------------------------------------------------------------------------------------------+
|              TỔNG HỢP TỌA ĐỘ VÀ THÔNG SỐ TRẮC LƯỢNG 11 CÁ THỂ SA MỘC DẦU TẠI TÂY CÔN LĨNH             |
+-----+--------+---------------+------------+-----------+-----------+----------+-----------+------------+
| Cây | Tuyến  | Xã / Tiểu khu | Tọa độ GPS | Độ cao(m) | Độ dốc(%) | $D_{1.3}$| $H_{vn}$  | Trạng thái |
| số  | OTC    |               |            |           |           | (cm)     | (m)       | tầng tán   |
+-----+--------+---------------+------------+-----------+-----------+----------+-----------+------------+
| 1   | OTC 2  | Túng Sán      | 0319252 -  | 1.060     | 35%       | 80.00    | 33.00     | Vượt tán   |
| 2   | (TK    | (TK 194B -    | 2519095    | 1.065     | 35%       | 48.00    | 22.00     | Tán chính  |
| 3   | 194B)  | 1735 ha)      | 0322466 -  | 1.118     | 35%       | 30.25    | 18.00     | Tán chính  |
| 4   | OTC 3  |               | 2516235    | 1.161     | 35%       | 47.40    | 22.00     | Tán chính  |
+-----+--------+---------------+------------+-----------+-----------+----------+-----------+------------+
| 5   | OTC 5  | Tả Sử Choóng  | 0320038 -  | 1.062     | 30%       | 67.00    | 30.00     | Vượt tán   |
| 6   | (TK    | (TK 204 -     | 2509612    | 1.062     | 30%       | 65.00    | 28.00     | Vượt tán   |
| 7   | 204)   | 55.5 ha)      | 0320070 -  | 1.082     | 30%       | 51.00    | 25.00     | Tán chính  |
| 8   |        |               | 2509603    | 1.096     | 30%       | 25.00    | 15.00     | Tán dưới   |
| 9   |        |               | 0320276 -  | 1.096     | 30%       | 20.00    | 12.00     | Tán dưới   |
|     |        |               | 2510824    |           |           |          |           |            |
+-----+--------+---------------+------------+-----------+-----------+----------+-----------+------------+
| 10  | OTC 7  | Bản Nhùng     | 0318251 -  | 970       | 25%       | 76.00    | 27.00     | Vượt tán   |
| 11  | (TK    | (TK 194A -    | 2511536    | 976       | 25%       | 62.00    | 20.00     | Tán chính  |
|     | 194A)  | 1098 ha)      |            |           |           |          |           |            |
+-----+--------+---------------+------------+-----------+-----------+----------+-----------+------------+

               PHÂN BỐ ĐƯỜNG KÍNH (D1.3) VÀ CHIỀU CAO (Hvn) 11 CÁ THỂ
               
  Hvn (m)
   35 |        [Cây 1: D=80, H=33]
   30 |                 [Cây 5: D=67, H=30]  [Cây 6: D=65, H=28]
   25 |                          [Cây 10: D=76, H=27]   [Cây 7: D=51, H=25]
   20 |        [Cây 2: D=48, H=22]   [Cây 4: D=47.4, H=22]   [Cây 11: D=62, H=20]
   15 |                 [Cây 3: D=30.25, H=18]   [Cây 8: D=25, H=15]
   10 |                          [Cây 9: D=20, H=12]
      +------------------------------------------------------------- D1.3 (cm)
      0       10      20      30      40      50      60      70      80

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

  • Tại OTC 2 (Túng Sán, 1000-1100m, Tổng 41 cây): $$\text{Công thức tổ thành} = 6.1Sm + 0.97Vt + 0.97Nh + 0.73Mlt + 0.73Lk - 0.48Smd$$
  • Tại OTC 3 (Túng Sán, 1100-1200m, Tổng 34 cây): $$\text{Công thức tổ thành} = 3.23Vt + 2.94Sm + 1.47Lk + 1.18Mlt + 0.58Kr - 0.58Smd$$
  • Tại OTC 5 (Tả Sử Choóng, 1000-1100m, Tổng 46 cây): $$\text{Công thức tổ thành} = 4.56Sm + 1.52Lk + 1.09Smd + 1.08Mlt + 0.87Qr + 0.86Cn$$
  • Tại OTC 7 (Bản Nhùng, 900-1000m, Tổng 38 cây): $$\text{Công thức tổ thành} = 3.15Sm + 2.1Lk + 1.84Xt + 1.31Qr + 1.05Vt + 0.52Smd$$

(Ký hiệu loài: Sm: Sa mộc; Smd: Sa mộc dầu; Vt: Vối thuốc; Nh: Nhội; Mlt: Mò lá tròn; Kr: Khế rừng; Qr: Quế rừng; Cn: Cọ nọt; Xt: Xoan ta; Lk: Loài khác).

3. Cấu trúc tầng thứ thẳng đứng

  • Tầng vượt tán ($H > 20\ m$): Các cá thể Sa mộc dầu cổ thụ lớn ($D_{1.3}$ lên đến 80 cm, $H_{vn}$ 33 m), tán hình tháp thưa, phân cành cao.
  • Tầng tán chính ($H = 10 - 20\ m$): Chiếm ưu thế bởi Vối thuốc (Schima wallichii), Sa mộc (Cunninghamia lanceolata), Quế rừng (Cinnamomum iners).
  • Tầng dưới tán ($H = 5 - 10\ m$): Các loài cây chịu bóng trung bình như Cọ nọt, Lá nến, Sổ.
  • Tầng thảm tươi cây bụi ($H < 5\ m$): Cỏ nhật, guột, dương xỉ, chè rừng và cây mạ chịu bóng.
                            MÔ HÌNH LÁT CẮT TẦNG THỨ RỪNG
                            
  Cao (m)
   35 |           / \  [Sa mộc dầu - Tầng vượt tán: Cây số 1 (H=33m, D=80cm)]
   30 |          /   \
   25 |         /     \                / \  [Sa mộc dầu - Cây số 5 (H=30m)]
   20 |        |       |    / \       /   \         / \  [Vối thuốc]
   15 |        |       |   /   \     /     \       /   \
   10 |        |       |  [Sa mộc]  [Quế rừng]    |     |  <-- TẦNG TÁN CHÍNH (10-20m)
    5 |        |       |  [Cọ nọt / Lá nến]       |     |  <-- TẦNG DƯỚI TÁN (5-10m)
    0 +========+=======+==========================+=====+===== <-- THẢM TƯƠI / CÂY BỤI

4. Thực trạng tái sinh tự nhiên

Trong toàn bộ diện tích $9.000\ m^2$ điều tra tại 9 ô tiêu chuẩn, không ghi nhận bất kỳ cây con tái sinh tự nhiên nào của Sa mộc dầu ở cấp chiều cao triển vọng ($H > 1.0\ m$).

+-------------------------------------------------------------------------+
|                  ĐÁNH GIÁ CHỈ SỐ TÁI SINH TỰ NHIÊN SA MỘC DẦU           |
+------------------------------------+------------------------------------+
| Chỉ số lâm học                     | Giá trị điều tra thực tế           |
+------------------------------------+------------------------------------+
| Mật độ cây mạ ($H < 30\ cm$)       | Xuất hiện rất ít rải rác           |
| Mật độ cây con triển vọng ($H > 1m$)| 0 cây/ha                           |
| Tỷ lệ cây con tốt                  | 0.0%                               |
| Tỷ lệ nảy mầm ngoài thảm mục dày   | < 1.0%                             |
+------------------------------------+------------------------------------+

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

Điểm mới về sinh học và lâm học

  1. Làm rõ đặc điểm sinh thái nón và hạt dị biệt: Phát hiện hạt Sa mộc dầu có cánh mỏng ($5\ mm \times 4\ mm$), nhưng hạt không bung tách khỏi nón trên cây khi chín. Nón rơi nguyên xuống nền rừng; nếu gặp thảm mục dày ẩm ướt, hạt sẽ nảy mầm thành cụm ngay trên quả nón cũ. Điều này giải thích hiện tượng loài chỉ mọc theo cụm và không thể mở rộng vùng phân bố tự nhiên dưới tán rừng kín.
  2. Khác biệt hình thái lá theo giai đoạn phát triển và giới tính cành:
    • Lá non: Mềm, mọc xoắn ốc, mép có răng cưa nhỏ, dài 5-7 cm, rộng 0.1-0.2 cm, màu xanh nhạt.
    • Lá mang nón cái: Mọc xoắn ốc dày đặc, xu hướng mọc thẳng đứng hướng thiên, dài 2.5-3 cm, rộng 0.2-0.3 cm, 2 dải khí khổng trắng rất rõ ở mặt dưới.
    • Lá mang nón đực: Xếp 2 dãy phẳng, dài 3-4 cm, rộng 0.2-0.25 cm, đầu lá tù mềm.
+---------------------------------------------------------------------------------------------------+
|                        SO SÁNH CÁC CÔNG TRÌNH NGHIÊN CỨU VỀ SA MỘC DẦU                            |
+----------------------+--------------------------+------------------------+------------------------+
| Thông số khảo sát    | Nghiên cứu này           | VQG Pù Mát             | Đài Loan               |
|                      | (Tây Côn Lĩnh, Hà Giang) | (Nguyễn Văn Sinh, 2009)| (Huang et al., 2008)   |
+----------------------+--------------------------+------------------------+------------------------+
| Đai cao phân bố      | 900 - 1.200 m            | 1.100 - 1.800 m        | 1.300 - 2.800 m        |
| Mật độ quần thể      | Cực kỳ thưa (11 cá thể)  | Rải rác theo cụm nhỏ   | Phân tán do khai thác  |
| Đặc tính tái sinh    | 0% cây con triển vọng    | Tái sinh mạ yếu        | Tái sinh nhân tạo 100% |
| Phương thức hạt nảy  | Nảy mầm trên nón rụng    | Nảy mầm trên nón rụng  | Xử lý tách hạt gieo ươm|
| Nhiệt độ thích ứng   | 20 - 23 °C               | 18 - 22 °C             | 17 - 22 °C             |
+----------------------+--------------------------+------------------------+------------------------+

Đóng góp cho khoa học và thực tiễn lâm nghiệp

  • Cung cấp cơ sở dữ liệu định vị địa lý (GPS) chính xác để Ban quản lý KBTTN Tây Côn Lĩnh đưa 11 cá thể cây mẹ vào chế độ bảo vệ nghiêm ngặt cấp 1.
  • Định hình hướng xử lý cơ học lâm sinh: Phá vỡ thảm mục, tạo khoảng trống ánh sáng (canopy gaps) và bắt buộc phải can thiệp thu hái nón để tách hạt nhân tạo trước khi gieo ươm.

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

Kịch bản và giải pháp lâm sinh can thiệp

       QUY TRÌNH KỸ THUẬT BẢO TỒN VÀ PHÁT TRIỂN NGUỒN GEN SA MỘC DẦU
       
+-------------------------+     +-------------------------+     +-------------------------+
| Bước 1: Thu hái nón chín|     | Bước 2: Tách hạt & Xử lý|     | Bước 3: Gieo ươm mô hình|
| - Thời gian: Tháng 3-5  | ==> | - Phơi nhẹ trong bóng râm| ==> | - Giá thể: Đất mùn rừng |
| - Hái nón chuyển nâu xám|     | - Tách hạt thủ công     |     | - Che bóng 50-70% giai  |
| - Không để nón chạm đất |     | - Xử lý nấm Trichoderma |     |   đoạn cây mạ non       |
+-------------------------+     +-------------------------+     +-------------------------+
                                                                             |
                                                                             v
+-------------------------+     +-------------------------+     +-------------------------+
| Bước 6: Trồng hỗn giao  |     | Bước 5: Mở tán rừng     |     | Bước 4: Xúc tiến in-situ|
| - Trồng kèm Sa mộc,     | <== | - Phát dọn thực bì 50m2 | <== | - Đưa cây con 1-2 tuổi  |
|   Vối thuốc theo tỷ lệ  |     |   quanh cây mẹ          |     |   trở lại sinh cảnh     |
|   tổ thành 3:4:3        |     | - Tạo lỗ thủng ánh sáng |     | - Đóng cọc rào bảo vệ   |
+-------------------------+     +-------------------------+     +-------------------------+

Phân tích chi phí - lợi ích và giá trị kinh tế

+---------------------------------------------------------------------------------------------------+
|                        BẢNG ĐỊNH LƯỢNG HIỆU QUẢ KINH TẾ VÀ XÃ HỘI                                 |
+-------------------------------+-----------------------------------+-------------------------------+
| Hạng mục                      | Khai thác tự phát (Hiện trạng)    | Quản lý bảo tồn & Trồng rừng  |
+-------------------------------+-----------------------------------+-------------------------------+
| Giá trị gỗ thương phẩm        | 50 - 80 triệu VNĐ/m3 (Gỗ lậu)     | Bảo tồn nguồn gen vĩnh viễn   |
| Tinh dầu chưng cất            | Khai thác tận diệt rễ củ          | Khai thác tỉa cành lá định kỳ |
| Nguồn thu cho cộng đồng       | Bất hợp pháp, rủi ro pháp lý cao  | Chi trả Dịch vụ Môi trường Rừng|
| Nguy cơ suy thoái sinh thái   | 100% cạn kiệt trong 5-10 năm      | Tăng 300% mật độ cây con      |
| Suất đầu tư gieo ươm bảo tồn  | 0 VNĐ                             | ~ 150 triệu VNĐ/năm/trạm      |
+-------------------------------+-----------------------------------+-------------------------------+

Lộ trình triển khai khuyến nghị

                           LỘ TRÌNH 3 GIAI ĐOẠN TRIỂN KHAI
                           
  2014 - 2016                   2016 - 2018                   2018 - 2022
+-----------------------+     +-----------------------+     +-----------------------+
|  GIAI ĐOẠN 1: BẢO VỆ  |     |  GIAI ĐOẠN 2: NHÂN TẠO|     |  GIAI ĐOẠN 3: NHÂN RỘNG|
|  - Cắm mốc 11 cây mẹ  | ==> |  - Vườn ươm 5.000 cây | ==> |  - Trồng 50 ha rừng   |
|  - Thành lập tổ tuần  |     |  - Thử nghiệm giâm hom|     |    hỗn giao đệm       |
|    rừng thôn bản      |     |  - Mở tán quanh cây mẹ|     |  - Phát triển sinh kế |
+-----------------------+     +-----------------------+     +-----------------------+

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

Hạn chế kỹ thuật và điều kiện nghiên cứu

  • Thiếu hụt dữ liệu hình thái hoa: Do khung thời gian nghiên cứu kết thúc vào tháng 4, chưa bao quát được thời điểm thụ phấn và ra nón đực/nón cái (tháng 9 - tháng 10).
  • Mẫu khảo sát hạn chế do tính chất quý hiếm: Tổng số cá thể tìm thấy chỉ dừng ở 11 cây mẹ, chưa đủ lớn để xây dựng mô hình hồi quy sinh trưởng phi tuyến đa biến.
  • Rào cản tiếp cận lập địa: Các dông núi hiểm trở ở độ cao $> 1.200\ m$ có vách đá dựng đứng chưa thể bao phủ toàn diện 100% diện tích các tiểu khu.

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

  • Ứng dụng chỉ thị phân tử sinh học (ADN Barcoding, chỉ thị microsatellite/SSR) để đánh giá độ đa dạng di truyền và mức độ cận huyết của 11 cá thể còn sót lại.
  • Thử nghiệm công nghệ nuôi cấy mô tế bào (In-vitro tissue culture) nhằm nhân bản nhanh các dòng cây mẹ có hoạt chất tinh dầu cao phục vụ bảo tồn ngoại vi.
  • Xây dựng trạm quan trắc tiểu khí hậu tự động tại xã Tả Sử Choóng để đo đạc chính xác biên độ nhiệt ẩm liên tục theo mùa.

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

+---------------------------------------------------------------------------------------------------+
|                           PHÂN TÍCH ĐỐI TƯỢNG HƯỞNG LỢI VÀ GIÁ TRỊ TÁC ĐỘNG                       |
+------------------------+--------------------------------------------------------------------------+
| Nhóm đối tượng         | Giá trị thụ hưởng định lượng & Định tính                                 |
+------------------------+--------------------------------------------------------------------------+
| Sinh viên & Giảng viên | Cung cấp tài liệu thực địa chuẩn xác, phương pháp lập 9 OTC tiêu chuẩn   |
| ngành Lâm nghiệp       | và bộ công thức toán học tính toán tổ thành $N_a$ phục vụ đào tạo.        |
+------------------------+--------------------------------------------------------------------------+
| Cán bộ Kiểm lâm &      | Bộ dữ liệu tọa độ GPS 11 cây mẹ, phân vùng 2.888,5 ha sinh cảnh trọng    |
| Ban quản lý KBTTN      | điểm thuộc 3 tiểu khu (194A, 194B, 204) phục vụ kiểm soát lâm tặc.       |
+------------------------+--------------------------------------------------------------------------+
| Cộng đồng người Dao,   | Cơ chế tham gia bảo vệ rừng có hưởng thù lao; hướng tới mô hình trồng     |
| H'Mông địa phương      | cây bản địa xen canh, tạo nguồn sinh kế thay thế việc đào rễ bán lậu.   |
+------------------------+--------------------------------------------------------------------------+
| Giới nghiên cứu sinh học| Giải mã tập tính nón quả không tự nở - tiền đề cho các công trình công   |
| & Công nghệ sinh học   | nghệ nhân giống loài Hạt trần nguy cấp tại Việt Nam.                     |
+------------------------+--------------------------------------------------------------------------+

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

1. Tại sao Sa mộc dầu không thể tự tái sinh tốt dưới tán rừng nguyên sinh Tây Côn Lĩnh?

Do hai nguyên nhân sinh thái chính: Thứ nhất, Sa mộc dầu là loài cây ưa sáng ở giai đoạn sau cây mạ, nhưng tán rừng kín tại đai cao 1000-1100m có độ tàn che lớn làm thiếu hụt bức xạ mặt trời. Thứ hai, nón cái khi chín không giải phóng hạt trên cành mà rơi nguyên nón xuống thảm mục dày, hạt bị nấm mốc hoặc thối rữa trước khi chạm được vào tầng đất khoáng.

2. Sự khác biệt giữa Sa mộc (Cunninghamia lanceolata) và Sa mộc dầu (Cunninghamia konishii) là gì?

Sa mộc (C. lanceolata) là loài cây trồng phổ biến, sinh trưởng nhanh, lá lớn và cứng hơn, nón dễ nở giải phóng hạt. Sa mộc dầu (C. konishii) là loài bản địa quý hiếm, gỗ có hàm lượng tinh dầu cao, mùi thơm ngát, lá có 2 dải khí khổng trắng rõ rệt ở mặt dưới và chu kỳ nón chín đặc thù kéo dài qua mùa đông.

3. Làm thế nào để phân biệt lá mang nón đực và lá mang nón cái của Sa mộc dầu?

Lá mang nón cái mọc xoắn ốc dày đặc, ngắn hơn ($2.5 - 3\ cm \times 0.2 - 0.3\ cm$), có xu hướng dựng thẳng đứng hướng lên ngọn. Lá mang nón đực xếp phẳng gần như 2 dãy, dài hơn ($3 - 4\ cm \times 0.2 - 0.25\ cm$) và có đầu tù mềm mại.

4. Đai cao và điều kiện địa hình tối ưu của Sa mộc dầu tại Hà Giang là bao nhiêu?

Loài phân bố tập trung tối ưu ở đai cao 1.000 m – 1.100 m, trên các sườn dốc từ 25% đến 35%, tại các lập địa rừng kín thường xanh mưa ẩm á nhiệt đới núi trung bình có nhiệt độ bình quân năm từ 20 – 23 °C.

5. Cần áp dụng biện pháp kỹ thuật cấp bách nào để nhân giống loài cây này?

Cần thu hái quả nón vào tháng 3 đến tháng 5 khi vỏ nón chuyển sang màu nâu xám, tách hạt thủ công, xử lý nấm bệnh bằng dung dịch sát khuẩn nhẹ, sau đó gieo trên luống đất cát pha mùn giàu dinh dưỡng có mái che bớt 50% ánh sáng trong 6 tháng đầu.

Kết luận

Đề tài “Nghiên cứu phân bố và đặc điểm lâm học của loài cây Sa mộc dầu (Cunninghamia konishii Hayata) tại Khu bảo tồn thiên nhiên Tây Côn Lĩnh - tỉnh Hà Giang” đã giải quyết trọn vẹn các yêu cầu khoa học và thực tiễn:

  • Khảo sát và lập hồ sơ quản lý cho 11 cá thể cây mẹ cổ thụ với đường kính thân $D_{1.3}$ đạt tới 80 cm và chiều cao $H_{vn}$ lên đến 33 m, định vị trên diện tích sinh cảnh 2.888,5 ha thuộc các tiểu khu 194A, 194B và 204.
  • Xây dựng thành công hệ thống công thức tổ thành tầng cây cao tại 9 ô tiêu chuẩn $1.000\ m^2$, xác định Sa mộc (Sm) và Vối thuốc (Vt) là các loài cây ưu thế sinh thái đồng hành quan trọng nhất.
  • Phát hiện cơ chế nghẽn sinh học trong quá trình phát tán nón quả và tình trạng suy kiệt cây con tái sinh ($0\ cây/ha$ có triển vọng).
  • Thiết lập khung giải pháp 3 trục: Lâm sinh kỹ thuật (mở tán, thu nón tách hạt) - Pháp lý kiểm soát (giao khoán tiểu khu) - Sinh kế cộng đồng vùng đệm.

[!IMPORTANT] Công tác bảo tồn Sa mộc dầu tại Tây Côn Lĩnh đòi hỏi hành động can thiệp kỹ thuật ngay lập tức trong mùa thu hái nón (tháng 3 - 5) kết hợp nguồn ngân sách ổn định để chi trả thù lao tuần rừng cho người dân bản địa, ngăn chặn nguy cơ tuyệt chủng của loài cây quý hiếm này.