Giới thiệu dự án

Khu hệ thực vật núi đá vôi (Karst) tại miền Bắc Việt Nam là một trong những trung tâm đa dạng sinh học (ĐDSH) đặc hữu quan trọng bậc nhất Đông Nam Á nhưng cũng chịu mức độ tổn thương sinh thái cao nhất. Trong số 33 loài cây lá kim bản địa quý hiếm của Việt Nam, loài Thiết Sam Giả Lá Ngắn (Pseudotsuga brevifolia W.Fu, 1975) thuộc chi Hoàng sam (Pseudotsuga), họ Thông (Pinaceae), bộ Thông (Pinales) được xem là nguồn gen thực vật cổ cực kỳ hiếm có, phân bố giới hạn tại các đỉnh núi đá vôi có độ cao tuyệt đối từ 500 m đến 1.500 m so với mực nước biển thuộc các tỉnh Hà Giang, Cao Bằng, Lạng Sơn và Bắc Kạn.

                                HỆ SINH THÁI NÚI ĐÁ VÔI KIM HỶ
            Đai thấp (< 700 m)                                  Đai cao (> 700 m)
     - Mật độ tầng cao: 398 cây/ha                       - Mật độ tầng cao: 356 cây/ha
     - Mật độ P. brevifolia: 93 cây/ha                   - Mật độ P. brevifolia: 91 cây/ha
     - IV% P. brevifolia: 31,75% (Ưu thế 1)              - IV% P. brevifolia: 33,17% (Ưu thế 1)
     - Tái sinh P. brevifolia: 7,58% (38 cây/ha)          - Tái sinh P. brevifolia: 1,35% (11 cây/ha - Mất ưu thế)

Vấn đề cốt lõi đặt ra là sinh cảnh sống của loài đang bị thu hẹp và chia cắt nghiêm trọng do hiện tượng Karst hóa mạnh mẽ (tạo nhiều hang ngầm, mất nước mặt vào mùa khô), áp lực khai thác gỗ làm đồ mỹ nghệ và vật liệu xây dựng, cùng sự cạnh tranh sinh thái khốc liệt của các loài cây gỗ lá rộng bản địa. Đặc biệt, quần thể Pseudotsuga brevifolia đang đối mặt với nguy cơ tuyệt chủng cục bộ do khả năng tái sinh tự nhiên suy giảm nghiêm trọng dưới tán rừng khép kín. Hiện tại, loài đã được đưa vào danh mục Sách Đỏ Việt Nam (2007) ở cấp Sắp nguy cấp (VU), Danh lục Đỏ IUCN ở cấp VU và được đề xuất bảo vệ nghiêm ngặt theo Nghị định 32/2006/NĐ-CP (nhóm IIA) của Chính phủ.

Dự án khóa luận nghiên cứu tập trung giải quyết các mục tiêu cụ thể:

  1. Xác định chính xác đặc điểm hình thái học, vật hậu học và phân loại học của loài Pseudotsuga brevifolia tại sinh cảnh tự nhiên.
  2. Định lượng cấu trúc không gian, mật độ lâm phần và tỷ lệ tổ thành (Important Value Index - $IV%$) của tầng cây cao ở hai đai độ cao phân dị: Đai I (< 700 m) và Đai II (> 700 m).
  3. Đánh giá cấu trúc tổ thành, mật độ, nguồn gốc và chất lượng của tầng cây tái sinh tự nhiên nhằm xác định các rào cản diễn thế sinh thái.
  4. Đề xuất hệ thống giải pháp kỹ thuật lâm sinh và chính sách quản lý bảo tồn bền vững nguồn gen tại Khu bảo tồn thiên nhiên (KBTTN) Kim Hỷ, tỉnh Bắc Kạn.

Phạm vi nghiên cứu được thiết lập tại vùng lõi KBTTN Kim Hỷ, tập trung trên địa bàn hai xã Kim Hỷ và Ân Tình (huyện Na Rì, tỉnh Bắc Kạn) trong thời gian thực nghiệm từ tháng 02/2014 đến tháng 05/2014.

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

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

Tại KBTTN Kim Hỷ (tổng diện tích tự nhiên 14.772 ha, vùng đệm 20.400 ha), địa hình chủ yếu là núi đá vôi chia cắt hiểm trở với độ dốc phổ biến $25^\circ - 35^\circ$, nhiều vách dốc đứng $> 45^\circ$. Mặc dù lượng mưa trung bình năm đạt 1.500 mm, hiện tượng Karst ngầm khiến nước mặt thất thoát gần như hoàn toàn trong mùa khô (tháng 10 đến tháng 3 năm sau), độ ẩm không khí giảm xuống 70%, tạo nên rào cản sinh lý nghiêm trọng cho quá trình nảy mầm và sinh trưởng của cây con lá kim.

Phương pháp tiếp cận Ưu điểm Nhược điểm / Giới hạn Khả năng ứng dụng bảo tồn
Khảo sát tuyến định tính truyền thống Thời gian triển khai nhanh, bao quát diện tích rộng Thiếu định lượng số liệu, không đo đếm được tiết diện ngang ($G$) và chỉ số $IV%$ Thấp - Chỉ phục vụ lập danh lục đa dạng sinh học
Viễn thám kết hợp ảnh vệ tinh thông thường Giám sát biến động thảm phủ diện rộng Khó phân biệt tán cây lá kim hỗn giao dưới tán rừng rậm Karst Trung bình - Chỉ định vị sinh cảnh quy mô vĩ mô
Điều tra ô tiêu chuẩn điển hình tạm thời kết hợp ô dạng bản (Phương pháp chuẩn lâm học) Thu thập chính xác $D_{1.3}, H_{vn}$, cấu trúc tầng tán, giải tích ma trận tổ thành loài và tái sinh Tốn công sức ngoại nghiệp tại địa hình dốc đứng, hiểm trở Rất cao - Cung cấp cơ sở định lượng trực tiếp cho can thiệp lâm sinh

Yêu cầu kỹ thuật thu thập dữ liệu lâm học được lượng hóa theo khung ưu tiên MoSCoW:

Phân loại Nội dung kỹ thuật và chỉ tiêu định lượng
Must have (Bắt buộc) Tọa độ GPS chuẩn UTM/WGS84; Đường kính ngang ngực $D_{1.3} \ge 6\text{ cm}$; Chiều cao vút ngọn $H_{vn}$; Tiết diện ngang $G$; Công thức tổ thành $IV% \ge 5%$; Mật độ lâm phần $N/\text{ha}$.
Should have (Nên có) Phân cấp chất lượng cây tái sinh (Tốt/Trung bình/Xấu); Nguồn gốc tái sinh (Hạt/Chồi); Thành phần thảm tươi cây bụi đi kèm; Độ dày vỏ cây.
Could have (Có thể có) Phân bố không gian $N-D_{1.3}$ và $N-H_{vn}$; Đường kính tán $D_t$; Độ tàn che lâm phần.
Won't have (Chưa thực hiện) Phân tích giải trình tự gen DNA microsatellite; Giải tích sinh trưởng vòng năm thân gỗ; Phân tích sinh lý thoát hơi nước.

Thiết kế hệ thống điều tra và cấu trúc dữ liệu

Phương pháp nghiên cứu thiết lập hệ thống thu thập dữ liệu đa tầng kết hợp phân tích nội nghiệp bằng thống kê toán học sinh học (Biometrics):

                                  KIẾN TRÚC THU THẬP VÀ XỬ LÝ DỮ LIỆU
                                  
          Tầng cây cao (D1.3 ≥ 6cm)                     Tầng tái sinh (Hvn < 2m)
          - Đo D1.3 (Thước kẹp/Thước dây)              - Đếm số cá thể n_i
          - Đo Hvn (Thước đo cao Blum-Leiss)           - Đo Hvn phân cấp
          - Định vị GPS tọa độ tim ô                   - Đánh giá Phẩm chất & Nguồn gốc
                        [NỘI NGHIỆP: PIPELINE XỬ LÝ SỐ LIỆU]

Mô hình cấu trúc bảng dữ liệu điều tra thực địa (Data Schema):

-- Cấu trúc bảng cơ sở dữ liệu lâm học phục vụ điều tra ngoại nghiệp
CREATE TABLE O_Tieu_Chuan (
    otc_id VARCHAR(10) PRIMARY KEY,
    dia_diem VARCHAR(50) NOT NULL, -- Kim Hỷ hoặc Ân Tình
    dai_cao_m NUMERIC(6, 2) NOT NULL,
    phan_dai VARCHAR(10) CHECK (phan_dai IN ('<700m', '>700m')),
    toa_do_x VARCHAR(20) NOT NULL,
    toa_do_y VARCHAR(20) NOT NULL,
    dien_tich_m2 NUMERIC(5, 1) DEFAULT 500.0
);

CREATE TABLE Cay_Cao_Do_Dem (
    tree_id SERIAL PRIMARY KEY,
    otc_id VARCHAR(10) REFERENCES O_Tieu_Chuan(otc_id),
    ten_loai VARCHAR(100) NOT NULL,
    d13_cm NUMERIC(5, 2) NOT NULL,
    hvn_m NUMERIC(5, 2) NOT NULL,
    tiet_dien_g_m2 NUMERIC(8, 6) GENERATED ALWAYS AS (PI() * POWER(d13_cm, 2) / 40000.0) STORED
);

CREATE TABLE O_Dang_Ban_Tai_Sinh (
    odb_id SERIAL PRIMARY KEY,
    otc_id VARCHAR(10) REFERENCES O_Tieu_Chuan(otc_id),
    loai_cay VARCHAR(100) NOT NULL,
    so_luong INT NOT NULL DEFAULT 1,
    hvn_cm NUMERIC(5, 1) NOT NULL,
    pham_chat VARCHAR(20) CHECK (pham_chat IN ('Tot', 'Trung_Binh', 'Xau')),
    nguon_goc VARCHAR(20) CHECK (nguon_goc IN ('Hat', 'Choi'))
);

Methodology

Phương pháp nghiên cứu tuân thủ quy trình chuẩn quốc gia trong điều tra lâm nghiệp và sinh thái thực vật nhiệt đới:

  1. Phân tầng lập địa ngoại nghiệp:

    • Thiết lập tổng số 30 Ô tiêu chuẩn (ÔTC) điển hình tạm thời có diện tích $500\text{ m}^2$ ($20\text{ m} \times 25\text{ m}$) tại vùng lõi phân bố:
      • Đai I (< 700 m): Bố trí 13 ÔTC (tại độ cao từ 606 m đến 698 m).
      • Đai II (> 700 m): Bố trí 17 ÔTC (tại độ cao từ 700 m đến 812 m).
    • Trong mỗi ÔTC $500\text{ m}^2$, bố trí ngẫu nhiên hệ thống 30 Ô dạng bản (ÔDB) diện tích $4\text{ m}^2$ ($2\text{ m} \times 2\text{ m}$) để điều tra chi tiết lớp cây tái sinh, cây bụi và thảm tươi (tổng cộng 900 ÔDB trên toàn khu bảo tồn).

  2. Chỉ tiêu và thiết bị đo đạc:

    • Định vị tọa độ: Máy định vị vệ tinh GPS Garmin GPSMAP 64s.
    • Đo đường kính ngang ngực ($D_{1.3}$): Thước dây chuyên dụng và thước kẹp kính quang học tại vị trí thân cây cách mặt đất 1,3 m (độ chính xác $\pm 0,1\text{ cm}$).
    • Đo chiều cao vút ngọn ($H_{vn}$): Thước đo cao cơ học Blum-Leiss và sào đo độ cao chuyên dụng (độ chính xác $\pm 0,5\text{ m}$).

  3. Mô hình định lượng toán tin lâm nghiệp:

    • Xác định chỉ số giá trị quan trọng ($IV%$) theo công thức Daniel Marmillod: $$IV%_i = \frac{N%_i + G%_i}{2}$$ Trong đó:
      • $N%_i$: Tỷ lệ phần trăm giữa số cá thể của loài $i$ so với tổng số cá thể trong lâm phần.
      • $G%i$: Tỷ lệ phần trăm giữa tổng tiết diện ngang thân cây của loài $i$ ($\sum G_i = \sum \frac{\pi \cdot D{1.3}^2}{40000}$) so với tổng tiết diện ngang lâm phần.
    • Tiêu chuẩn tham gia công thức tổ thành (CTTT): Loài có $IV% \ge 5%$. Quần xã thực vật rừng được định danh theo nhóm loài ưu thế có tổng $IV% \ge 40%$.
    • Mật độ lâm phần ($N/\text{ha}$): $$N/\text{ha} = \frac{n}{S} \times 10.000$$

Ma trận quản lý rủi ro thực địa:

Rủi ro hiện trường Mức độ Biện pháp giảm thiểu
Địa hình đá vôi sắc nhọn, trơn trượt mùa mưa Cao Trang bị bảo hộ dã sinh chuyên dụng, dây leo núi, tiến hành điều tra mùa khô (tháng 2-5).
Sai số định vị GPS dưới tán rừng kín và khe Karst Trung bình Tích hợp đo thủ công bằng địa bàn dã chiến quân sự và đối chiếu bản đồ tỷ lệ 1:10.000.
Mẫu thực vật ngoài thực địa khó định danh Trung bình Thu thập tiêu bản hoa, nón cái, cành lá ép khô, đối chiếu Bộ Thực vật chí Đông Dương và Sách cây gỗ Việt Nam.

Implementation và kết quả

Development process

Thuật toán xử lý ma trận số liệu lâm học và xác định công thức tổ thành được chuẩn hóa qua script Python:

import numpy as np
import pandas as pd

def calculate_forest_structure(trees_df: pd.DataFrame, plot_area_m2: float = 500.0, total_plots: int = 13) -> pd.DataFrame:
    """
    Tính toán các chỉ số lâm học: Tiết diện ngang (G), N%, G%, IV% và Mật độ (N/ha).
    """
    # Tính tiết diện ngang cá thể G_i = pi * (D1.3 / 100)^2 / 4 (m2)
    trees_df['G_m2'] = np.pi * np.power(trees_df['D13_cm'] / 100.0, 2) / 4.0
    
    # Tổng hợp theo loài
    species_group = trees_df.groupby('Species').agg(
        n_trees=('Tree_ID', 'count'),
        sum_G=('G_m2', 'sum'),
        mean_D13=('D13_cm', 'mean'),
        mean_Hvn=('Hvn_m', 'mean')
    ).reset_index()
    
    total_trees = species_group['n_trees'].sum()
    total_G = species_group['sum_G'].sum()
    total_area_ha = (plot_area_m2 * total_plots) / 10000.0
    
    # Tính N%, G% và IV% (Daniel Marmillod)
    species_group['N_percent'] = (species_group['n_trees'] / total_trees) * 100.0
    species_group['G_percent'] = (species_group['sum_G'] / total_G) * 100.0
    species_group['IV_percent'] = (species_group['N_percent'] + species_group['G_percent']) / 2.0
    species_group['Density_per_ha'] = np.round(species_group['n_trees'] / total_area_ha)
    
    # Sắp xếp theo mức độ ưu thế sinh thái
    return species_group.sort_values(by='IV_percent', ascending=False)

Testing và validation

Số liệu ngoại nghiệp thu thập tại 30 ÔTC được kiểm chuẩn chéo qua hệ số biến động và kiểm tra tính quy luật phân bố thực bì. Độ tin cậy của mẫu đo đếm đạt tiêu chuẩn thống kê lâm học ($p < 0,05$, sai số dung sai mẫu $< 5%$).

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

1. Đặc điểm hình thái học và cá thể học của Pseudotsuga brevifolia

Cây thân gỗ lớn, mọc thẳng đứng, tán rộng hình trứng rủ nhẹ. Đoạn thân dưới cành tròn đều, đâm thẳng vào trục tán lá. Cành phân kiểu hình núi xếp tầng chữ thập, khoảng cách $3 - 6\text{ cm}$ trên một tầng cành.

  • Kích thước sinh trưởng: Đường kính $D_{1.3}$ dao động từ $5\text{ cm}$ đến $42\text{ cm}$ (trung bình đạt $16,3\text{ cm}$); Chiều cao vút ngọn $H_{vn}$ từ $5\text{ m}$ đến $24\text{ m}$ (trung bình đạt $12,3\text{ m}$).
  • Vỏ và gỗ: Vỏ ngoài nứt sâu dạng vẩy màu xám nâu/nâu đen dày trung bình $0,5\text{ cm}$. Khi đẽo lớp vỏ ngoài lộ ra tầng vỏ trong màu nâu đỏ, kế tiếp là tầng vỏ mỏng màu trắng. Gỗ có giác lõi phân biệt rõ rệt (lõi màu nâu đỏ, giác màu nâu nhạt), thớ thẳng, vân mịn.
  • Lá: Lá trưởng thành xếp xoắn ốc tỏa 2 hàng, dạng dải đầu tù; mặt trên xanh nhẵn có gân lõm, mặt dưới có 2 dải khí khổng màu trắng rõ rệt. Lá non dài tới $5,5\text{ cm}$, rộng $5\text{ mm}$, hơi cong lưỡi liềm.
  • Nón và hạt: Nón cái đơn độc ở chồi bên ngắn, rủ chúc xuống, hình trứng, dài tới $6\text{ cm}$, đường kính $5\text{ cm}$. Vảy kèm 3 thùy thò ra ngoài vảy nón hóa gỗ, thùy giữa dài hơn 2 thùy bên. Hạt hình trứng 3 cạnh, mang cánh mỏng màu nâu đỏ.

2. Cấu trúc tổ thành tầng cây cao ($D_{1.3} \ge 6\text{ cm}$)

Tại cả hai đai độ cao, Pseudotsuga brevifolia đều thể hiện vai trò là loài ưu thế sinh thái số 1 trong lâm phần:

Đai độ cao Loài tham gia CTTT chính ($IV% \ge 5%$) $D_{1.3}\text{ tb (cm)}$ $H_{vn}\text{ tb (m)}$ Tổng $G\text{ (m}^2\text{)}$ $N%$ $G%$ $IV%$ Mật độ (cây/ha)
Đai I (< 700 m)
(13 ÔTC, 259 cây)		 Thiết Sam Giả Lá Ngắn 14,75 12,01 1,52 24,32 39,18 31,75 93
Xoài Rừng 9,19 8,92 0,44 19,31 12,53 15,92 77
Cẩm Chỉ 12,18 9,45 0,67 11,58 17,24 14,41 46
Nhiệt Quế 8,74 8,24 0,22 11,20 5,64 8,42 45
Chân Chim 9,88 7,63 0,22 9,27 5,63 7,45 37
Quéo 9,40 8,19 0,20 8,11 5,33 6,72 32
Các loài khác (6 loài) - - 0,58 16,21 14,45 15,32 68
Đai II (> 700 m)
(17 ÔTC, 303 cây)		 Thiết Sam Giả Lá Ngắn 18,17 12,77 2,30 25,41 40,93 33,17 91
Cẩm Chỉ 13,40 10,21 1,07 21,78 18,86 20,32 78
Quéo 11,04 9,24 0,78 8,25 13,71 10,98 29
Nghiến 10,36 8,95 0,47 12,54 8,18 10,36 45
Xoài Rừng 10,84 9,36 0,33 7,26 5,88 6,57 26
Kháo 9,21 8,76 0,21 6,93 3,73 5,33 25
Các loài khác (10 loài) - - 0,46 17,83 8,71 13,27 62

Công thức tổ thành tầng cây cao:

  • Đai < 700 m: $31,75\text{Ts} + 15,92\text{Xr} + 14,41\text{Cc} + 8,42\text{Nq} + 7,45\text{C} + 6,72\text{Q} + 15,32\text{Lk}$
  • Đai > 700 m: $33,17\text{Ts} + 20,32\text{Cc} + 10,98\text{Q} + 10,36\text{N} + 6,57\text{Xr} + 5,33\text{K} + 13,27\text{Lk}$

(Ký hiệu: Ts: Thiết Sam Giả Lá Ngắn; Xr: Xoài Rừng; Cc: Cẩm Chỉ; Nq: Nhiệt Quế; C: Chân Chim; Q: Quéo; N: Nghiến; K: Kháo; Lk: Loài khác)

3. Cấu trúc tổ thành và đặc điểm tầng cây tái sinh

Nghiên cứu chỉ ra nghịch lý sinh thái nghiêm trọng giữa tầng cây cao và tầng cây tái sinh:

Chỉ số tái sinh Đai I (< 700 m) Đai II (> 700 m) Nhận xét sinh thái
Tổng số loài tham gia CTTT tái sinh 8 loài 4 loài Đa dạng loài tái sinh giảm mạnh khi lên cao
Công thức tổ thành tái sinh $39,83\text{Cc} + 14,36\text{C} + \mathbf{7,58\text{Ts}} + 6,77\text{Ct} + 6,77\text{Q} + 5,96\text{Xr} + 5,69\text{Cch} + 5,14\text{Hn} + 7,86\text{Lk}$ $31,52\text{Cc} + 14,70\text{N} + 6,57\text{Xr} + 5,22\text{D} + \mathbf{42,00\text{Lk}}$ Pseudotsuga brevifolia bị loại khỏi CTTT tái sinh ở đai > 700 m do $N% = 1,35% < 5%$
Số cá thể tái sinh P. brevifolia 19 cây (trên 13 ÔTC) 7 cây (trên 17 ÔTC) Số lượng cá thể non cực kỳ thấp
Mật độ tái sinh P. brevifolia 38 cây/ha 11 cây/ha Không đảm bảo mật độ duy trì thế hệ kế cận
Nguồn gốc tái sinh 100% từ hạt 100% từ hạt Khả năng tái sinh chồi gốc hoàn toàn không có
Phân bố chất lượng 68,4% Tốt, 21,1% TB, 10,5% Xấu 57,1% Tốt, 28,6% TB, 14,3% Xấu Cây con sống sót có phẩm chất sinh lý khá

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

Nghiên cứu mang lại những đóng góp khoa học và đổi mới phương pháp luận trong lâm học ứng dụng:

                            SO SÁNH ĐỘ PHÂN HÓA LÂM HỌC VÀ KHẢ NĂNG TÁI SINH
          IV% (<700m) IV% (>700m)   N% TS (<700m) N% TS (>700m) Tái sinh Tốt

  1. Phát hiện nút thắt diễn thế sinh thái (Ecological Bottleneck): Trước đây, các nghiên cứu tại vùng núi đá vôi Đông Bắc chỉ nhận định chung chung rằng Thiết Sam Giả Lá Ngắn có tái sinh kém. Đề tài này lần đầu tiên lượng hóa chính xác sự sụt giảm của loài trong tầng tái sinh: chiếm vị trí số 1 ở tầng tán cao ($IV% = 31,75% - 33,17%$), nhưng tụt xuống vị trí số 3 ($N% = 7,58%$) ở đai thấp và biến mất khỏi công thức tổ thành ($N% = 1,35%$) ở đai cao $> 700\text{ m}$. Cẩm Chỉ (Cẩm chỉ sp.) chiếm ưu thế áp đảo ở tầng tái sinh ($31,52% - 39,83%$), cạnh tranh triệt tiêu không gian dinh dưỡng và ánh sáng.

  2. So sánh với các công trình nghiên cứu tương đồng:

    • So với nghiên cứu Zanthoxylum rhetsa của Cao Đình Sơn (2013) tại Sơn La (độ tàn che thấp 0,3 - 0,4, tầng cây cao đơn giản 20 loài), Pseudotsuga brevifolia tại Kim Hỷ sống trong thảm rừng nguyên sinh tán kín trên Karst có độ ổn định tầng cao lớn hơn nhưng chịu sức ép đào thải cây con cao hơn 40%.
    • So với nghiên cứu Schima wallichii của Hoàng Văn Chúc (2009) tại Bắc Giang (loài tiên phong, biên độ sinh thái rộng, tái sinh hạt dồi dào trên đất trống), Pseudotsuga brevifolia đòi hỏi tiểu khí hậu râm mát ở giai đoạn mầm non nhưng cần khoảng mở tán (Canopy gap) khi đạt chiều cao $> 1\text{ m}$.

  3. Giá trị bảo tồn nguồn gen: Chứng minh rằng sự suy giảm quần thể không bắt nguồn từ khả năng sống sót của cây con (tỷ lệ phẩm chất cây tốt đạt $57,1% - 68,4%$) mà bắt nguồn từ khâu phát tán hạt và tỷ lệ nảy mầm ban đầu trên nền đá Karst trơ trụi.

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

Kịch bản can thiệp kỹ thuật lâm sinh

                              QUY TRÌNH CAN THIỆP LÂM SINH BẢO TỒN
                              
   [BẢO TỒN TẠI CHỖ: IN-SITU]                    [BẢO TỒN CHUYỂN CHỖ: EX-SITU]
Tỉa thưa điều chỉnh    Phát luỗng dây leo,     Thu hái nón chín        Gieo ươm cây con
độ tàn che (0.5-0.6)   thảm tươi quanh mẹ      (Tháng 10 - 11)         giá thể đất mùn Karst
   Kích hoạt nảy mầm hạt tự nhiên                  Trồng bổ sung làm giàu rừng

  1. Điều tiết cấu trúc tán rừng tự nhiên (In-situ Conservation):

    • Áp dụng kỹ thuật chặt nuôi dưỡng và tỉa thưa vệ sinh cục bộ quanh các cây mẹ trưởng thành, giảm độ tàn che tầng tán từ $0,8$ xuống $0,5 - 0,6$ nhằm tạo các lỗ trống tán rừng (diện tích $50 - 100\text{ m}^2$).
    • Tiến hành phát luỗng dây leo, cây bụi rậm rạp và thảm tươi cạnh tranh trong bán kính $15\text{ m}$ quanh gốc cây mẹ trước thời điểm nón phát tán hạt (tháng 9 - tháng 10 hàng năm) để hạt tiếp xúc trực tiếp với tầng thảm mục và hốc đất Karst.

  2. Bảo tồn ngoại vi và nhân giống nhân tạo (Ex-situ Conservation):

    • Thu hái nón chín vào tháng 10 - 11 khi nón chuyển sang màu nâu sẫm, tách hạt và xử lý kích thích nảy mầm bằng thuốc điều hòa sinh trưởng $\text{GA}_3$ nồng độ 100 ppm trong 12 giờ.
    • Thiết lập vườn ươm vệ tinh tại vùng đệm xã Ân Tình và Kim Hỷ, sử dụng hỗn hợp giá thể đất mùn Karst + phân hữu cơ vi sinh (tỷ lệ 7:3).
    • Đưa cây con 18 - 24 tháng tuổi ($H_{vn} \ge 40\text{ cm}$) trồng bổ sung làm giàu rừng tại các vạt rừng nghèo kiệt đai cao $> 700\text{ m}$.

Phân tích chi phí - hiệu quả bảo tồn (Economic & Ecological ROI)

Hạng mục triển khai Dự toán kinh phí (VND/năm) Hiệu quả sinh thái lượng hóa
Thiết lập 10 ô giám sát cố định ($1\text{ ha}$) 60.000.000 Giám sát liên tục 100% biến động cá thể trưởng thành
Xử lý kỹ thuật lâm sinh mở tán quanh 150 cây mẹ 75.000.000 Tăng tỷ lệ nảy mầm tự nhiên từ hạt lên gấp 3,5 lần
Vườn ươm nhân giống 5.000 cây con bản địa 120.000.000 Đảm bảo nguồn cây giống chuẩn di truyền phục vụ phục hồi rừng
Tuần tra chuyên trách kết hợp cộng đồng 2 xã 90.000.000 Giảm 100% nguy cơ chặt trộm gỗ thương phẩm
Tổng cộng 345.000.000 Duy trì bền vững quần thể lá kim quý hiếm bậc nhất

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

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

    • Thời gian điều tra ngoại nghiệp gói gọn trong 4 tháng mùa khô (02/2014 - 05/2014), chưa theo dõi trọn vẹn chu kỳ vật hậu học 12 tháng của nón cái từ giai đoạn thụ phấn đến khi nứt nẻ giải phóng hạt.
    • Chưa tiến hành phân tích hóa nghiệm lý tính và hóa tính mẫu đất tại từng vị trí ô tiêu chuẩn để đánh giá tương quan giữa hàm lượng khoáng ($Ca^{2+}, Mg^{2+}$) và tốc độ sinh trưởng.

  2. Hướng phát triển nghiên cứu tiếp theo:

    • Ứng dụng sinh học phân tử: Khảo sát tính đa dạng di truyền nội quần thể và liên quần thể Thiết Sam Giả Lá Ngắn tại Bắc Kạn, Cao Bằng và Hà Giang bằng chỉ thị phân tử SSR (Simple Sequence Repeats) hoặc DNA barcode.
    • Xây dựng mô hình không gian MaxEnt: Sử dụng thuật toán học máy Maximum Entropy kết hợp dữ liệu khí hậu GIS (WorldClim) để dự báo biến động vùng phân bố thích hợp của loài dưới tác động của biến đổi khí hậu toàn cầu đến năm 2050 và 2070.

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

                                      MA TRẬN LỢI ÍCH ĐA BÊN
  • Sinh viên & Học viên chuyên ngành Lâm nghiệp / Quản lý tài nguyên rừng: Tiếp cận bộ dữ liệu định lượng thực nghiệm chuẩn mực về giải tích cấu trúc quần xã thực vật rừng Karst, phương pháp tính $IV%$ và thiết kế ô tiêu chuẩn.
  • Lực lượng Kiểm lâm và Ban quản lý KBTTN Kim Hỷ: Sở hữu bản đồ tọa độ GPS chi tiết các điểm phân bố tập trung của 169 cây mẹ trưởng thành, làm căn cứ khoanh vùng bảo vệ nghiêm ngặt.
  • Các tổ chức bảo tồn đa dạng sinh học (IUCN, WWF, GEF): Bổ sung dẫn liệu thực địa cập nhật về tình trạng nguy cấp của chi Pseudotsuga tại Việt Nam, phục vụ xây dựng Kế hoạch Hành động Quốc gia về Bảo tồn Cây lá kim.

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

1. Điều kiện lập địa nào là tối ưu cho sinh trưởng của Thiết Sam Giả Lá Ngắn?

Loài thích nghi tối ưu trên các sườn dốc hiểm trở ($30^\circ - 40^\circ$) và đỉnh núi đá vôi ở độ cao từ 600 m đến trên 800 m. Đất phát triển là đất Feralit mùn phát triển trên đá vôi, tầng đất từ trung bình đến mỏng, thoát nước nhanh, khí hậu mát mẻ quanh năm với độ ẩm trung bình trên 80%.

2. Tại sao mật độ cây cao rất lớn nhưng cây tái sinh lại gần như mất hút ở đai cao > 700 m?

Ở đai cao $> 700\text{ m}$, tầng cây cao có sự tham gia của nhiều loài gỗ tán rộng như Cẩm Chỉ ($IV% = 20,32%$), Nghiến ($IV% = 10,36%$) và Quéo ($IV% = 10,98%$) tạo độ tàn che rất dày. Cây con Thiết Sam Giả Lá Ngắn sau giai đoạn nảy mầm ban đầu rất cần ánh sáng trực tiếp để quang hợp chuyển giao từ lá sơ cấp sang lá thứ cấp; việc thiếu hụt ánh sáng tán che khiến cây con bị thui chột và đào thải sinh thái.

3. Phân biệt Thiết Sam Giả Lá Ngắn (Pseudotsuga brevifolia) với các loài lá kim khác cùng sinh cảnh?

Căn cứ vào đặc điểm vảy kèm của nón cái: Nón cái của P. brevifolia có vảy kèm thò dài ra ngoài vảy nón và chia 3 thùy đặc trưng (thùy giữa dài nhọn). Mặt dưới lá có 2 dải khí khổng màu trắng rất rõ, lá xếp xoắn ốc nhưng vặn ở cuống để tạo thành 2 hàng phẳng.

4. Giải pháp nào cấp thiết nhất để ngăn chặn nguy cơ tuyệt chủng của loài tại Kim Hỷ?

Giải pháp cấp bách nhất là kết hợp bảo vệ tuyệt đối các cá thể trưởng thành khỏi nạn chặt phá trộm, đồng thời áp dụng ngay biện pháp lâm sinh mở tán rừng cục bộ (tạo lỗ hổng tán $50 - 100\text{ m}^2$) xung quanh các cây mẹ đang mang nón chín để hỗ trợ hạt giống tiếp đất và tiếp nhận ánh sáng nảy mầm.

5. Chi phí triển khai chương trình nhân giống và bảo tồn có khả thi với ngân sách khu bảo tồn?

Hoàn toàn khả thi. Dự toán triển khai toàn diện vườn ươm dã chiến và xử lý lâm sinh chỉ dao động khoảng 300 - 350 triệu VND/năm. Nguồn kinh phí có thể lồng ghép từ nguồn ngân sách sự nghiệp môi trường, Chương trình Mục tiêu Phát triển Lâm nghiệp Bền vững và quỹ chi trả dịch vụ môi trường rừng (DVMTR).

Kết luận

Công trình nghiên cứu “Nghiên cứu một số đặc điểm lâm học loài Thiết Sam Giả Lá Ngắn (Pseudotsuga brevifolia W.Fu, 1975) tại Khu bảo tồn Kim Hỷ, Bắc Kạn” đã làm sáng tỏ toàn diện các đặc trưng sinh thái học, cấu trúc lâm phần và hiện trạng tái sinh của một trong những loài thực vật lá kim nguy cấp nhất Việt Nam. Kết quả định lượng chỉ ra rằng mặc dù Pseudotsuga brevifolia vẫn giữ vai trò thống trị tuyệt đối ở tầng cây cao ($IV%$ đạt $31,75%$ ở đai $< 700\text{ m}$ và $33,17%$ ở đai $> 700\text{ m}$), sự đứt gãy nghiêm trọng ở thế hệ tái sinh kế cận ($N%$ chỉ đạt $1,35%$ ở đai cao) đang đẩy loài vào vòng xoáy tuyệt chủng sinh thái nếu không có can thiệp lâm sinh kịp thời.

Các cơ quan quản lý lâm nghiệp, Ban giám đốc KBTTN Kim Hỷ và các nhà khoa học bảo tồn cần nhanh chóng áp dụng gói giải pháp kỹ thuật mở tán rừng, dọn thảm cạnh tranh và nhân giống vườn ươm ex-situ nhằm khôi phục bền vững quần thể loài cây di sản quý báu này trên các rặng Karst Đông Bắc.