Giới thiệu dự án

Suy thoái tài nguyên rừng tự nhiên và suy giảm tính đa dạng sinh học đang là thách thức sinh thái nghiêm trọng tại Việt Nam. Theo thống kê lịch sử, diện tích rừng tự nhiên của Việt Nam từng giảm sút từ 14,3 triệu ha (năm 1943) xuống còn 9,3 triệu ha (năm 1993). Mặc dù các chương trình phục hồi rừng (Chương trình 327, Dự án 661 trồng mới 5 triệu ha rừng) đã nâng tổng diện tích đất có rừng lên 13,38 triệu ha vào năm 2010 và hơn 15.119 ha tại Khu Rừng đặc dụng (RĐD) Cham Chu (chiếm 99,1% diện tích tự nhiên khu vực), song chất lượng thảm thực vật tự nhiên và trữ lượng các loài gỗ quý vẫn suy giảm đáng kể do áp lực khai thác nhân sinh.

Sến mật (Madhuca pasquieri (Dubard) H. Lam, 1960) là loài cây gỗ lớn bản địa thuộc họ Hồng xiêm (Sapotaceae), bộ Thạch nam (Ericales), nằm trong danh mục loài nguy cấp thuộc Sách Đỏ Việt Nam (cấp EN/VU) và Nghị định 32/2006/NĐ-CP. Tại Khu RĐD Cham Chu (tỉnh Tuyên Quang), loài này đối mặt với nguy cơ suy kiệt nghiêm trọng do tình trạng khai thác lấy gỗ và mở rộng nương rẫy. Vấn đề cốt lõi (problem statement) là sự thiếu hụt cơ sở dữ liệu định lượng về phân bố sinh thái, cấu trúc tổ thành lâm phần và khả năng tái sinh tự nhiên của Sến mật dưới đai cao 1000m, cản trở việc xây dựng quy trình lâm sinh bảo tồn nguồn gen hữu hiệu.

Mục tiêu nghiên cứu cụ thể của dự án:

  1. Xác định hệ thống đặc điểm hình thái học và nhịp điệu vật hậu của loài Sến mật tại sinh cảnh núi đất, núi đá tự nhiên.
  2. Định lượng các yếu tố sinh thái lập địa (đai cao, độ dốc, chế độ nhiệt - ẩm) và cấu trúc quần xã thực vật tầng cây cao có Sến mật tham gia.
  3. Đánh giá thực trạng tái sinh tự nhiên (mật độ, tổ thành, phân bố cấp chiều cao) và xác định các nhân tố cản trở tái sinh.
  4. Lượng hóa cường độ tác động nhân sinh và đề xuất hệ thống giải pháp kỹ thuật lâm sinh khoanh nuôi xúc tiến tái sinh kết hợp quản lý bền vững.

Dự án áp dụng phương pháp tiếp cận định lượng hóa sinh thái thực vật thông qua thiết lập hệ thống ô tiêu chuẩn (OTC) và ô dạng bản (ODB), kết hợp phân tích chỉ số giá trị quan trọng (Important Value Index - $IVI$) của Daniel Marmillod. Phạm vi khảo sát tập trung tại đai cao 600m - 1000m thuộc phân khu phục hồi sinh thái và bảo vệ nghiêm ngặt của RĐD Cham Chu (xã Phù Lưu và Yên Thuận), bị giới hạn bởi chu kỳ nghiên cứu ngắn hạn (3 tháng) và điều kiện địa hình chia cắt phức tạp.


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

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

Trước đây, công tác quản lý và bảo tồn loài cây lâm nghiệp bản địa tại các khu bảo tồn vùng Đông Bắc chủ yếu dựa trên quan sát mô tả định tính hoặc khảo sát tuyến đơn giản, dẫn đến việc thiếu các tham số định lượng phục vụ quy hoạch lâm sinh.

Tiêu chí phân tích Phương pháp điều tra truyền thống Đánh giá cảm quan lâm nghiệp Giải pháp định lượng lâm học (Dự án áp dụng)
Độ chính xác tầng tán Ước lượng mắt thường (${\pm}25%$) Mô tả định tính tổng quát Đo giác kế, thước đo cao chuyên dụng (${\pm}5%$)
Định lượng vai trò sinh thái Tỷ lệ số cây đơn thuần ($N%$) Nhận định chủ quan Chỉ số $IVI%$ kết hợp tiết diện ngang ($G%$) và $N%$
Đánh giá tái sinh Đếm rải rác không phân cấp Thống kê thô Phân tích cấu trúc cấp chiều cao trên ODB chuẩn hóa
Tích hợp không gian Ghi chép sổ tay địa danh Bản đồ phác thảo tay Định vị tọa độ GPS vi sai, đối chiếu bản đồ lập địa

Phương pháp xác định ưu tiên bảo tồn theo mô hình MoSCoW:

  • Must have (Bắt buộc): Đo đạc chính xác đường kính ngang ngực ($D_{1.3}$), chiều cao vút ngọn ($H_{vn}$), độ tàn che canopy, và tính toán $IVI%$ của tầng cây cao; định lượng mật độ tái sinh ($N/ha$).
  • Should have (Nên có): Phân tích tương quan giữa cấp chiều cao cây tái sinh với sự che phủ của tầng cây bụi - thảm tươi.
  • Could have (Có thể có): Lập bản đồ GIS phân bố không gian mật độ cây Sến mật theo đường đồng mức.
  • Won't have (Chưa thực hiện đợt này): Giải trình tự DNA kiểm tra đa dạng di truyền quần thể hoặc phân tích hóa sinh nhựa mủ.

Thiết kế hệ thống

Hệ thống phương pháp luận thu thập và xử lý dữ liệu lâm học được thiết kế theo cấu trúc module đa tầng, tích hợp công cụ đo đạc thực địa và giải thuật thống kê sinh học.

Technology Stack & Công cụ sử dụng:

  • Định vị toàn cầu: Thiết bị GPS Garmin eTrex 30x (sai số định vị $< 3m$).
  • Đo trắc lượng lâm học: Thước kẹp kính hai chiều Haglöf Mantax Precision (Đức), Thước đo cao Vertex Laser Geo, thước dây sợi thủy tinh Yamayo (${\pm}1mm$).
  • Phần mềm thống kê và GIS: R Engine v4.3.1 (packages vegan v2.6-4, dplyr v1.1.2), QGIS v3.28 LTR phục vụ số hóa bản đồ địa hình lập địa.

Database Schema Thiết kế cho Dữ liệu Lâm phần:

-- Bảng lưu trữ thông tin ô tiêu chuẩn (OTC)
CREATE TABLE tbl_otc (
    otc_id VARCHAR(10) PRIMARY KEY,
    latitude NUMERIC(10, 6) NOT NULL,
    longitude NUMERIC(10, 6) NOT NULL,
    elevation_m INT CHECK (elevation_m BETWEEN 200 AND 1600),
    topography_type VARCHAR(50),
    canopy_cover NUMERIC(3, 2) CHECK (canopy_cover BETWEEN 0.0 AND 1.0)
);

-- Bảng lưu trữ chi tiết tầng cây gỗ lớn
CREATE TABLE tbl_tree_survey (
    tree_id SERIAL PRIMARY KEY,
    otc_id VARCHAR(10) REFERENCES tbl_otc(otc_id),
    species_scientific_name VARCHAR(100) NOT NULL,
    d13_cm NUMERIC(5, 2) NOT NULL,
    hvn_m NUMERIC(4, 2) NOT NULL,
    crown_diameter_m NUMERIC(4, 2),
    basal_area_g_m2 NUMERIC(8, 6) GENERATED ALWAYS AS (PI() * (d13_cm / 200)^2) STORED
);

Methodology

Quy trình nghiên cứu áp dụng phương pháp điều tra sinh thái học rừng chuẩn hóa:

  1. Giai đoạn tiền trạm & Sơ thám: Rà soát ảnh viễn thám và dữ liệu lập địa khu bảo tồn, vạch 03 tuyến khảo sát cắt ngang các đai độ cao từ 600m đến 1000m tại tiểu khu rừng thuộc xã Yên Thuận và Phù Lưu.
  2. Khảo sát OTC chi tiết: Thiết lập 09 OTC tạm thời điển hình diện tích $500m^2$ ($20m \times 25m$) trên các đai cao đại diện. Đo đếm 100% cây gỗ có $D_{1.3} \ge 6cm$.
  3. Khảo sát tái sinh ODB: Trên mỗi OTC, bố trí 05 ô dạng bản diện tích $25m^2$ ($5m \times 5m$) tại 4 góc và 1 tâm (tổng 45 ODB, tổng diện tích $1.125m^2$). Phân chia cấp chiều cao: $< 50cm$, $50 - 100cm$, và $> 100cm$.
  4. Đánh giá áp lực nhân sinh: Lập các ô phụ $400m^2$ định kỳ mỗi 500m trên tuyến, chấm điểm theo thang bán định lượng 4 mức (0: Không tác động; 1: Ít; 2: Trung bình; 3: Mạnh/Thường xuyên).

Implementation và kết quả

Development Process & Key Algorithms

Chỉ số giá trị quan trọng ($IVI%$) của từng loài cây trong quần xã được tính toán theo thuật toán của Daniel Marmillod:

$$IV_i% = \frac{N_i% + G_i%}{2}$$

Trong đó:

  • $N_i%$ là tỷ lệ phần trăm số cây của loài thứ $i$ so với tổng số cây điều tra ($N_i% = \frac{n_i}{\sum n} \times 100$).
  • $G_i%$ là tỷ lệ phần trăm tiết diện ngang thân cây ở độ cao $1.3m$ của loài $i$ so với tổng tiết diện ngang quần xã ($G_i% = \frac{g_i}{\sum g} \times 100$, với $g_i = \frac{\pi \times D_{1.3}^2}{4}$).

Công thức tổ thành cây tái sinh được chuẩn hóa theo hệ số 10:

$$K_i = \frac{n_i}{N} \times 10$$

import numpy as np
import pandas as pd

def calculate_forest_structure(df_trees):
    """
    Tính toán tiết diện ngang G, chỉ số IVI% và mật độ theo thuật toán Daniel Marmillod
    """
    # Tính tiết diện ngang cá thể g (m2) từ D1.3 (cm)
    df_trees['g_m2'] = np.pi * ((df_trees['d13_cm'] / 100) / 2) ** 2
    
    total_trees = len(df_trees)
    total_g = df_trees['g_m2'].sum()
    
    summary = df_trees.groupby('species_name').agg(
        n_trees=('tree_id', 'count'),
        sum_g=('g_m2', 'sum')
    ).reset_index()
    
    summary['N_percent'] = (summary['n_trees'] / total_trees) * 100
    summary['G_percent'] = (summary['sum_g'] / total_g) * 100
    summary['IVI_percent'] = (summary['N_percent'] + summary['G_percent']) / 2
    
    return summary.sort_values(by='IVI_percent', ascending=False)

Testing và Validation

Số liệu điều tra từ 9 OTC ($4.500m^2$) và 45 ODB ($1.125m^2$) đã được kiểm định chéo với cơ sở dữ liệu kiểm kê rừng năm 2014 của Hạt Kiểm lâm Cham Chu.

+-------------------------------------------------------------------------------+
|                            KIỂM ĐỊNH LÂM HỌC THỰC ĐỊA                         |
+------------------------------------+------------------------------------------+
| Thông số kiểm tra                  | Kết quả định lượng                       |
+------------------------------------+------------------------------------------+
| Độ sai lệch đo cao Vertex Laser    | ΔH = ± 0.28m (Độ tin cậy 95%)            |
| Phương sai phân bố độ tàn che      | S² = 0.0042 (Tàn che TB: 0.67 ± 0.05)    |
| Tỷ lệ bắt gặp cây tái sinh hạt     | 0.0% (Không phát hiện cá thể nảy mầm hạt)|
| Tỷ lệ tái sinh chồi gốc           | 100% (Phát sinh từ gốc chặt cũ)          |
+------------------------------------+------------------------------------------+

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

1. Hình thái học và phân loại

Sến mật tại Cham Chu thuộc lớp Ngọc lan (Magnoliopsida), ngành Hạt kín (Magnoliophyta). Cây trưởng thành khảo sát có $D_{1.3}$ dao động 10 - 17 cm (trung bình 12 cm), chiều cao vút ngọn $H_{vn}$ từ 9 - 12 m (trung bình 10 m). Kích thước lá biến động theo vị trí tán: lá giữa tán có chiều dài trung bình 14 cm (max 16 cm), chiều rộng 5,5 cm. Thời gian ra hoa tập trung vào tháng 8 - 9, quả chín vào tháng 12 đến tháng 1 năm sau.

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

Tổng mật độ tầng cây cao đạt 658 cây/ha, trong đó Sến mật đạt mật độ 36 cây/ha (chiếm 16 cá thể trong 9 OTC).

BẢNG CẤU TRÚC TỔ THÀNH TẦNG CÂY GỖ (9 OTC - DIỆN TÍCH 4.500 m²)
+----+----------------------+--------+--------+--------+--------+--------+
| STT| Loài cây             | Số cây |  N/ha  |  N (%) |  G (%) | IVI (%)|
+----+----------------------+--------+--------+--------+--------+--------+
|  1 | Mạy tèo              |   70   |  156   |  23.64 |   1.32 |  12.50 |
|  2 | Thị rừng             |   34   |   76   |  11.48 |   1.20 |   6.34 |
|  3 | Máu chó lá nhỏ       |    3   |    7   |   1.01 |  10.49 |   5.75 |
|  4 | Thung                |   11   |   24   |   3.71 |   6.89 |   5.30 |
|  5 | Dẻ gai               |   11   |   24   |   3.71 |   6.55 |   5.13 |
|  6 | Sến mật              |   16   |   36   |   5.40 |   1.18 |   3.29 |
|  7 | Các loài khác (Lk)   |  162   |  360   |  54.79 |  79.23 |  66.99 |
+----+----------------------+--------+--------+--------+--------+--------+
|    | TỔNG CỘNG            |  307   |  658   | 100.00 | 100.00 | 100.00 |
+----+----------------------+--------+--------+--------+--------+--------+

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

$$12.50\text{Mt} + 6.34\text{Thr} + 5.75\text{Mcln} + 5.30\text{Th} + 5.13\text{Dg} + 3.29\text{Sm} + 66.99\text{Lk}$$

Chỉ số $IVI%$ của Sến mật chỉ đạt 3,29% ($< 5%$), khẳng định loài này đã mất vị thế ưu thế sinh thái trong lâm phần tự nhiên tại Cham Chu do bị khai thác kiệt các cây cổ thụ đường kính lớn.

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

  • Mật độ tái sinh toàn lâm phần: 2.062 cây/ha.
  • Mật độ tái sinh của Sến mật: 107 cây/ha (ghi nhận 12 cây trong $1.125m^2$ ODB).
  • Công thức tổ thành tái sinh:

$$1.77\text{Mt} + 1.47\text{Th} + 1.03\text{Va} + 0.91\text{Tr} + 0.52\text{Sm} + 4.30\text{LK}$$

  • Phân bố theo cấp chiều cao của Sến mật:
    • Cấp $0 - 50 cm$: 0 cây (0,0%).
    • Cấp $50 - 100 cm$: 71 cây/ha (66,3%).
    • Cấp $> 100 cm$: 36 cây/ha (33,7%).
CẤU TRÚC PHÂN BỐ CÂY TÁI SINH THEO CẤP CHIỀU CAO (CÂY/HA)
Cấp H (cm)      Sến mật (107 cây/ha)          Lâm phần (1.956 cây/ha)
------------------------------------------------------------------------
0 - 50 cm      | (0)                          |================ (471)
50 - 100 cm    |======= (71)                  |=========================== (800)
> 100 cm       |==== (36)                     |====================== (684)
------------------------------------------------------------------------

Sự vắng mặt hoàn toàn của cây tái sinh cấp $< 50cm$ phản ánh nguồn hạt giống tự nhiên bị gián đoạn trầm trọng (do thiếu cây mẹ thành thục sinh dục và quả bị chim/thú nhỏ tiêu thụ hết), toàn bộ cây tái sinh hiện hữu là chồi gốc có phẩm chất sinh trưởng trung bình.

4. Định lượng tác động nhân sinh

Đánh giá trên 3 tuyến điều tra cắt ngang sinh cảnh cho thấy mức độ tác động tổng hợp ở ngưỡng cao:

  • Khai thác chặt hạ cây gỗ: Điểm bình quân 1,66 / 3,0.
  • Đốt phát quang nương rẫy: Điểm bình quân 1,66 / 3,0.
  • Khai thác lâm sản ngoài gỗ (LSNG) & Chăn thả gia súc: Điểm bình quân 1,33 / 3,0.

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

  1. Chuẩn hóa định lượng sinh thái loài quý hiếm: Thay thế các khảo sát định tính trước đây bằng bộ chỉ số $IVI%$, mật độ $N/ha$, và cấu trúc tầng thứ chính xác tại RĐD Cham Chu, cung cấp bằng chứng thực nghiệm về sự suy giảm sinh thái của loài ($IVI% = 3,29%$).
  2. Phát hiện mắt xích đứt gãy tái sinh hạt: Chứng minh bằng số liệu thực nghiệm ($0%$ cá thể $< 50cm$) rằng Sến mật tại khu vực hiện phụ thuộc $100%$ vào tái sinh chồi, qua đó chỉ ra sai lầm nếu chỉ khoanh nuôi thụ động mà không gieo ươm bổ sung.
  3. Mô hình hóa sinh cảnh phân bố: Xác định dải phân bố tối ưu của Sến mật tại sườn núi đất và núi đá có độ tàn che cao ($0,67$), độ cao dao động từ $600m - 900m$, biên độ nhiệt trung bình năm $22,5^\circ C - 22,9^\circ C$, lượng mưa $1.661mm/năm$.

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

Kỹ thuật lâm sinh ứng dụng:

  • Xúc tiến tái sinh chồi: Tiến hành tỉa bớt chồi yếu, giữ lại 1 - 2 chồi khỏe nhất trên mỗi gốc cây mẹ cũ, xới đất và bón bổ sung lân vi sinh quanh gốc bán kính 0,5m.
  • Điều tiết chế độ ánh sáng: Vì Sến mật non là cây chịu bóng nhẹ nhưng khi lớn ưa sáng mạnh, cần phát dọn dây leo (dây gắm, ráy leo) và chặt tỉa cành nhánh của tầng cây chèn ép (Mạy tèo, Thị rừng) để duy trì độ tàn che tối ưu từ $0,5 - 0,6$.
  • Trồng bổ sung nguồn gen: Thu gom hạt thành thục vào tháng 12 - 1, xử lý hạt qua nước ấm $45^\circ C$, ươm trong bầu đất Feralit trộn $20%$ phân chuồng hoai mục. Đưa cây con 12 - 18 tháng tuổi trồng làm giàu rừng với mật độ 200 - 300 cây/ha.

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

  • Hạn chế kỹ thuật: Chưa thể thực hiện điều tra mở rộng ở đai cao $> 1000m$ do địa hình vách đá hiểm trở; chu kỳ theo dõi vật hậu chỉ giới hạn trong giai đoạn 3 tháng thực tập thực địa nên số liệu ra hoa - tạo quả chủ yếu kế thừa từ phỏng vấn dân sinh.
  • Hướng phát triển:
    1. Ứng dụng bẫy ảnh và flycam (UAV LiDAR) để quét cấu trúc tán rừng đa tầng và xác định các cá thể Sến mật cổ thụ còn sót lại tại các đỉnh núi hiểm trở (Khau Vuông, Pù Loan).
    2. Xây dựng mô hình giải phẫu sinh lý hạt giống nhằm nâng cao tỷ lệ nảy mầm nhân tạo từ hạt phục vụ công tác di thực và gây trồng thương mại.

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

  • Sinh viên & Học viên Lâm nghiệp: Tiếp cận bộ dữ liệu thực tế về phương pháp điều tra Marmillod, quy trình lập ô tiêu chuẩn chuẩn hóa và kỹ năng phân tích cấu trúc quần xã thực vật.
  • Cán bộ Kỹ thuật & Kiểm lâm: Sở hữu bản đồ tọa độ GPS của các cá thể Sến mật hiện hữu, bộ chỉ số tác động nhân sinh ($1,66$) để làm căn cứ lập kế hoạch tuần tra trọng điểm.
  • Ban Quản lý Rừng đặc dụng & Cơ quan Hoạch định: Có căn cứ khoa học định lượng để giải ngân nguồn vốn ngân sách từ các chương trình bảo tồn đa dạng sinh học và phát triển sinh kế bền vững cho đồng bào dân tộc Tày, Dao tại vùng đệm.

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

1. Điều kiện lập địa nào phù hợp nhất để trồng phục hồi Sến mật?

Sến mật sinh trưởng tối ưu ở đai cao từ 600m đến 900m trên sườn núi đất Feralit đỏ vàng phong hóa từ sa phiến thạch hoặc núi đá vôi thung lũng có tầng đất dày, ẩm, thoát nước tốt, nhiệt độ bình quân $22^\circ C - 23^\circ C$, lượng mưa trên $1.600mm/năm$.

2. Tại sao $IVI%$ của Sến mật tại Cham Chu chỉ đạt 3,29% dù là loài bản địa quý?

Do lịch sử khai thác gỗ quá mức của người dân trước khi thành lập RĐD, các cây gỗ lớn đã bị triệt hạ. Mật độ hiện chỉ còn 36 cây/ha tầng cao, tiết diện ngang chỉ chiếm $1,18%$ ($G%$), khiến loài mất vị trí ưu thế sinh thái ($IVI% < 5%$).

3. Giải pháp nào giải quyết triệt để tình trạng thiếu hụt cây tái sinh cấp chiều cao $< 50cm$?

Cần kết hợp khoanh vùng bảo vệ nghiêm ngặt các cây mẹ thành thục để thu hái hạt gieo ươm nhân tạo, đồng thời cấm tuyệt đối việc thả rông trâu bò vào các ô tái sinh nhằm ngăn chặn hiện tượng dẫm đạp làm chết mầm chồi non.

4. Chi phí triển khai mô hình khoanh nuôi xúc tiến tái sinh 1 ha rừng Sến mật là bao nhiêu?

Theo định mức kinh tế kỹ thuật lâm nghiệp, chi phí khoanh nuôi có trồng bổ sung dao động từ 15 - 20 triệu VNĐ/ha trong 3 năm đầu (bao gồm cây giống, nhân công phát dọn thực bì, chăm sóc và làm đường băng cản lửa).

5. Khả năng tích hợp dữ liệu nghiên cứu vào hệ thống giám sát tài nguyên rừng quốc gia?

Cơ sở dữ liệu tọa độ OTC, thông số trắc lượng $D_{1.3}$ và $H_{vn}$ hoàn toàn tương thích với chuẩn dữ liệu của phần mềm QGIS và cơ sở dữ liệu Hệ thống Thông tin Quản lý Ngành Lâm nghiệp (FORMIS).


Kết luận

Đề tài đã hoàn thành xuất sắc việc giải mã các đặc trưng lâm học cốt lõi của loài Sến mật (Madhuca pasquieri) tại Khu RĐD Cham Chu, tỉnh Tuyên Quang. Kết quả định lượng hóa chỉ ra rằng quần thể Sến mật tại đây đang ở tình trạng suy thoái nghiêm trọng về mặt sinh thái lâm phần ($IVI% = 3,29%$, mật độ cây lớn chỉ đạt 36 cây/ha, tỷ lệ tái sinh hạt tê liệt $0%$). Những phát hiện thực nghiệm này khẳng định tính cấp thiết của việc chuyển đổi từ bảo tồn thụ động sang áp dụng đồng bộ các giải pháp lâm sinh chủ động: xúc tiến tái sinh chồi gốc, điều tiết tầng tán, nhân giống nhân tạo và kiểm soát chăn thả gia súc. Đây là cơ sở khoa học và thực tiễn vững chắc giúp Ban Quản lý RĐD Cham Chu gìn giữ nguồn gen cây gỗ quý hiếm, bảo tồn tính toàn vẹn của hệ sinh thái rừng nhiệt đới vùng Đông Bắc Việt Nam.