Giới thiệu dự án

Rừng tự nhiên lá rộng thường xanh đóng vai trò cốt lõi trong việc duy trì cân bằng sinh thái, bảo tồn nguồn gen bản địa và cung cấp các dịch vụ môi trường rừng quan trọng. Theo thống kê của ngành lâm nghiệp, mặc dù tổng diện tích rừng tại Việt Nam có xu hướng gia tăng nhờ các chính sách đóng cửa rừng tự nhiên và trồng rừng mới, chất lượng rừng tự nhiên giàu và trung bình lại đang suy giảm nghiêm trọng do tác động lịch sử của khai thác kiệt quệ và biến đổi sinh cảnh.

Tại khu vực chuyển tiếp sinh thái đặc thù giữa Tây Nguyên, Đông Nam Bộ và Đồng bằng sông Cửu Long, Vườn Quốc gia (VQG) Lò Gò - Xa Mát (tỉnh Tây Ninh) với tổng diện tích quản lý 29.997,87 ha đang đối mặt với nhiều áp lực suy thoái sinh cảnh, xâm lấn đất nông nghiệp và biến đổi khí thái. Tuy nhiên, các dữ liệu định lượng chuyên sâu về cấu trúc quần thụ và chỉ số đa dạng sinh thái của trạng thái rừng lá rộng thường xanh trung bình tại đây còn thiếu tính hệ thống.

+------------------------------------------------------------------------------------+
|                                BỐI CẢNH DỰ ÁN                                      |
|  - Địa điểm: VQG Lò Gò - Xa Mát, Tây Ninh (29.997,87 ha)                           |
|  - Trạng thái: Rừng tự nhiên lá rộng thường xanh trung bình                        |
|  - Vấn đề: Thiếu mô hình định lượng cấu trúc & chỉ số đa dạng loài phục vụ quản lý |
+------------------------------------------------------------------------------------+

Vấn đề nghiên cứu (Problem Statement)

Công tác quản lý bảo vệ rừng tại VQG Lò Gò - Xa Mát gặp phải các điểm nghẽn kỹ thuật cụ thể:

  1. Thiếu hụt các tham số toán học chuẩn xác mô phỏng quy luật phân bố số cây theo đường kính ($N/D_{1.3}$) và chiều cao ($N/H_{vn}$) để dự báo quy luật sinh trưởng và tái sinh tự nhiên.
  2. Chưa lượng hóa được cấu trúc tổ thành loài ($IV_i%$) và chỉ số tương quan hỗn giao tầng cây cao ($D_{1.3} \ge 6\text{ cm}$), dẫn đến các giải pháp lâm sinh nuôi dưỡng rừng còn mang tính cảm tính.
  3. Thiếu hệ số đa dạng sinh học định lượng ($d, H', J', \lambda'$) làm thước đo đánh giá độ ổn định của quần xã thực vật (QXTV).

Mục tiêu dự án

  1. Xác định kết cấu tổ thành loài và họ thực vật: Phân tích định danh thành phần loài cây gỗ lớn ($D_{1.3} \ge 6\text{ cm}$) trên hệ thống ô tiêu chuẩn (OTC).
  2. Lượng hóa các đặc trưng cấu trúc lâm phần: Tính toán mật độ bình quân ($N/\text{ha}$), tổng tiết diện ngang ($G/\text{ha}$), trữ lượng tầng cây đứng ($M/\text{ha}$) và hệ số hỗn giao ($K$).
  3. Mô hình hóa toán học cấu trúc không gian: Thử nghiệm và kiểm định mức độ tương thích của các hàm phân bố xác suất (Weibull, Phân bố Chuẩn, Meyer, Khoảng cách) đối với $N%/D_{1.3}$ và $N%/H_{vn}$.
  4. Đánh giá định lượng đa dạng sinh học: Tính toán các chỉ số Margalef ($d$), Shannon-Weiner ($H'$), Pielou ($J'$), Simpson ($\lambda'$) và phân tích tương quan quần xã.
  5. Đề xuất giải pháp lâm sinh ứng dụng: Xây dựng phương án bảo tồn và phục hồi sinh thái phù hợp với Đề án Quản lý rừng bền vững giai đoạn 2021–2030.

Giải pháp kỹ thuật và kết quả kỳ vọng

  • Phương pháp: Thiết lập 12 OTC tạm thời điển hình diện tích $1.000\text{ m}^2$ ($25\text{ m} \times 40\text{ m}$), tổng diện tích điều tra 1,2 ha đại diện cho trạng thái rừng tự nhiên lá rộng thường xanh trung bình.
  • Chỉ số kỳ vọng: Mật độ đạt từ 650–750 cây/ha; trữ lượng lâm phần $M \ge 150\text{ m}^3/\text{ha}$ (đúng tiêu chí trạng thái rừng trung bình theo Thông tư 33/2018/TT-BNNPTNT); chỉ số Shannon $H' \ge 2,5$; xác định chính xác các loài ưu thế sinh thái có $IV_i% > 5%$.
  • Phạm vi và giới hạn: Nghiên cứu tập trung vào tầng cây gỗ có đường kính ngang ngực $D_{1.3} \ge 6\text{ cm}$; không bao gồm tầng cây tái sinh dưới $6\text{ cm}$, cây bụi, thảm tươi và dây leo ký sinh.

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 PHÁP NGHIÊN CỨU HIỆN HÀNH                      |
+------------------------------------+-----------------------------+-----------------------+
| Tiêu chí so sánh                   | Phương pháp truyền thống    | Giải pháp đề tài      |
+------------------------------------+-----------------------------+-----------------------+
| Phương pháp thu thập dữ liệu       | Mục trắc, ước lượng cảm quan| OTC định vị GPS + B-L |
| Mô hình hóa toán học cấu trúc      | Biểu đồ cột tĩnh, định tính | Hàm Weibull & Normal  |
| Kiểm định độ tin cậy thống kê      | Không kiểm định hoặc sơ sài | Trắc nghiệm Chi-square|
| Đánh giá đa dạng sinh học          | Đếm số loài đơn thuần       | Bộ 4 chỉ số sinh thái |
+------------------------------------+-----------------------------+-----------------------+

Ma trận ưu - nhược điểm các phương pháp tiếp cận

Phương pháp Ưu điểm Nhược điểm Ứng dụng thực tế
Mục trắc thị giác truyền thống Chi phí thấp, tốc độ điều tra nhanh. Sai số lớn ($> 30%$), không thể mô phỏng động thái. Kiểm kê sơ bộ quy mô lớn.
Viễn thám thuần túy (Remote Sensing) Bao phủ diện tích rộng, đo lường tán rừng tốt. Sai số tầng dưới tán cao, không định danh chính xác loài. Đánh giá độ che phủ tổng quan.
Lâm trắc định lượng kết hợp mô hình hóa (Đề tài áp dụng) Độ chính xác cao, kiểm định toán học nghiêm ngặt, định lượng $IV_i%$. Tốn công lao động ngoại nghiệp thực địa. Quy hoạch lâm sinh, bảo tồn nguồn gen.

Ưu tiên yêu cầu theo mô hình MoSCoW

  • Must have (Bắt buộc): Tọa độ GPS 12 OTC; đo $D_{1.3}$, $H_{vn}$ cho 100% cá thể $D_{1.3} \ge 6\text{ cm}$; tính toán $IV_i%$, $G$, $M$, $K$; kiểm định hàm Weibull cho $N%/D_{1.3}$ qua Chi-square ($\chi^2$).
  • Should have (Nên có): Xây dựng ma trận tương quan giữa các loài; biểu đồ phân tầng $H_{vn}$ theo 4 lớp; tính toán chỉ số Caswell-$V$.
  • Could have (Có thể có): Phân tích liên hệ giữa tính chất thổ nhưỡng đất xám phù sa cổ với phân bố tầng cây gỗ lớn.
  • Won't have (Chưa thực hiện): Đánh giá sinh khối ngầm (rễ cây) và đo đếm hàm lượng carbon hấp thụ trong đất.

Thiết kế hệ thống

  [ Ngoại nghiệp: 12 OTC (1.000m2) ] --(Thu thập D1.3, Hvn, GPS)--> [ Nhập liệu & Làm sạch Excel 2010 ]
                                                                               |
                               +-----------------------------------------------+
                               |                                               |
                               v                                               v
        [ Statgraphics Centurion XV.0 ]                             [ Biodiversity Pro 2.0 ]
        - Mô phỏng Weibull (N%/D1.3)                                 - Tính Margalef (d), Shannon (H')
        - Mô phỏng Chuẩn (N%/Hvn)                                   - Tính Pielou (J'), Simpson (λ')
        - Trắc nghiệm Chi-square (χ2)                               - Đường cong K-Dominance
                               |                                               |
                               +-----------------------+-----------------------+
                                                       |
                                                       v
                                      [ MapInfo 15.0 & Global Mapper v20 ]
                                      - Số hóa bản đồ phân bố OTC
                                      - Bản đồ hiện trạng lâm phần

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

  • Microsoft Excel 2010: Xử lý cơ sở dữ liệu sơ cấp, tính toán các chỉ tiêu sinh trắc lâm học ($N, G, V, IV_i%$).
  • Statgraphics Centurion XV.0: Khớp các hàm phân bố xác suất lý thuyết (Weibull, Normal, Meyer, Poisson), tính sai số tiêu chuẩn ($S_x$), phương sai ($S^2$), hệ số biến động ($CV%$) và trắc nghiệm Chi-square ($\chi^2$).
  • Biodiversity Pro 2.0: Xử lý ma trận dữ liệu loài, tính toán các chỉ số đa dạng sinh học và lập ma trận khoảng cách tương đồng Bray-Curtis.
  • MapInfo Professional 15.0 & Global Mapper v20.0: Xây dựng bản đồ phân bố vị trí không gian các ô tiêu chuẩn và bình đồ lâm phận VQG.
  • Thiết bị hiện trường: Thước dây chuyên dụng đo chu vi/đường kính (độ chính xác 0,5 cm), thước đo cao Blume-Leiss (sai số $< 1%$), máy định vị GPS Garmin GPSMAP 64s.

Cơ sở dữ liệu và cấu trúc công thức tính toán

-- Cấu trúc bảng cơ sở dữ liệu điều tra ô tiêu chuẩn (OTC_Data)
CREATE TABLE OTC_Survey_Data (
    tree_id INT PRIMARY KEY AUTO_INCREMENT,
    otc_id INT NOT NULL,
    species_vietnamese VARCHAR(100) NOT NULL,
    species_scientific VARCHAR(100),
    family_name VARCHAR(100) NOT NULL,
    d13_cm DECIMAL(5,2) NOT NULL,
    hvn_m DECIMAL(4,2) NOT NULL,
    basal_area_m2 DECIMAL(8,6),
    stem_volume_m3 DECIMAL(8,6),
    gps_lat DECIMAL(10,7),
    gps_long DECIMAL(10,7)
);

Hệ thống công thức toán học và sinh trắc lâm học chuẩn hóa:

  1. Tiết diện ngang thân cây ($G$): $$G = \frac{\pi}{4} \times D_{1.3}^2 = \frac{\pi}{40000} \times (D_{1.3}\text{ [cm]})^2 \quad (\text{m}^2)$$
  2. Thể tích thân cây đứng ($V$): $$V = G \times H_{vn} \times f_{1.3} \quad (\text{với hình số } f_{1.3} = 0,45)$$
  3. Chỉ số giá trị quan trọng ($IV_i%$): $$IV_i% = \frac{N% + G% + F%}{3}$$ (Trong đó: $N%$ là tỷ lệ số cá thể, $G%$ là tỷ lệ tiết diện ngang, $F%$ là tỷ lệ tần suất xuất hiện qua 12 OTC).
  4. Hàm phân bố xác suất Weibull: $$f(x) = \alpha \cdot \lambda \cdot x^{\alpha-1} \cdot e^{-\lambda x^\alpha}$$ (Trong đó: $\alpha$ là tham số hình dạng biểu thị độ lệch; $\lambda$ là tham số quy mô biểu thị độ nhọn).
  5. Bộ chỉ số đa dạng sinh thái:
    • Margalef: $d = \frac{S - 1}{\ln(N)}$
    • Shannon-Weiner: $H' = -\sum_{i=1}^S \left(\frac{n_i}{N}\right) \ln\left(\frac{n_i}{N}\right)$
    • Pielou: $J' = \frac{H'}{\ln(S)}$
    • Simpson: $\lambda' = \sum_{i=1}^S \frac{n_i(n_i - 1)}{N(N - 1)}$

Phương pháp luận (Methodology)

[ Giai đoạn 1: Chuẩn bị & Thu thập ] (Tháng 9 - 10/2022) -> Bản đồ hiện trạng, kế thừa tài liệu VQG
                     |
[ Giai đoạn 2: Điều tra ngoại nghiệp ] (Tháng 11 - 12/2022) -> Lập 12 OTC (1.2 ha), đo 865 cây
                     |
[ Giai đoạn 3: Phân tích nội nghiệp ] (Tháng 01/2023) -> Chạy mô hình Weibull, Normal, tính IVi% & H'
                     |
[ Giai đoạn 4: Tổng hợp & Đề xuất ] (Tháng 02/2023) -> Báo cáo khoa học & giải pháp lâm sinh

Quản trị rủi ro ngoại nghiệp

  • Rủi ro cháy rừng mùa khô: Mùa khô tại Tây Ninh (tháng 11 đến tháng 4) có lượng bốc hơi cao ($> 65%$ cả năm), nhiệt độ đỉnh điểm $34^\circ\text{C}$. Biện pháp: Phối hợp trạm kiểm lâm tuần tra, trang bị thiết bị định vị GPS tránh lạc tuyến rừng khộp/rừng trảng.
  • Sai số đo chiều cao tán kín: Tán rừng thường xanh phân tầng phức tạp che khuất đỉnh ngọn. Biện pháp: Sử dụng thước Blume-Leiss đo mẫu 3–5 cây đại diện trên mỗi cấp kính, mục trắc kết hợp kiểm tra chéo góc ngắm.

Implementation và kết quả

Quá trình triển khai và thuật toán xử lý

Thuật toán xử lý ma trận tổ thành và kiểm định phân bố được chuẩn hóa theo quy trình sau:

import numpy as np
import pandas as pd

def calculate_forest_metrics(df_trees, total_otc_area_ha=1.2, num_otcs=12):
    # 1. Tính tiết diện ngang (G) và thể tích (V) cho từng cá thể
    f13 = 0.45
    df_trees['g_m2'] = (np.pi / 4) * (df_trees['d13_cm'] / 100) ** 2
    df_trees['v_m3'] = df_trees['g_m2'] * df_trees['hvn_m'] * f13
    
    # 2. Tổng hợp theo loài
    total_trees = len(df_trees)
    total_g = df_trees['g_m2'].sum()
    
    species_summary = df_trees.groupby('species').agg(
        n_stems=('tree_id', 'count'),
        sum_g=('g_m2', 'sum'),
        otc_presence=('otc_id', 'nunique')
    ).reset_index()
    
    # 3. Tính tỷ lệ % và IVi%
    species_summary['N_pct'] = (species_summary['n_stems'] / total_trees) * 100
    species_summary['G_pct'] = (species_summary['sum_g'] / total_g) * 100
    species_summary['F_pct'] = (species_summary['otc_presence'] / num_otcs) * 100
    species_summary['IV_pct'] = (species_summary['N_pct'] + species_summary['G_pct'] + species_summary['F_pct']) / 3
    
    # 4. Tính mật độ, G và M bình quân trên 1 ha
    stand_density_ha = total_trees / total_otc_area_ha
    stand_g_ha = total_g / total_otc_area_ha
    stand_volume_ha = df_trees['v_m3'].sum() / total_otc_area_ha
    
    return species_summary.sort_values(by='IV_pct', ascending=False), {
        'N_ha': stand_density_ha,
        'G_ha': stand_g_ha,
        'M_ha': stand_volume_ha
    }

Kiểm định và mức độ thích ứng mô hình (Testing & Validation)

Quá trình kiểm định mức độ tương thích giữa tần số thực nghiệm và tần số lý thuyết của hàm phân bố được thực hiện qua trắc nghiệm Chi-square ($\chi^2$) với bậc tự do $df = m - p - 1$:

$$\chi^2_{\text{tính}} = \sum_{i=1}^m \frac{(f_{ti} - f_{li})^2}{f_{li}}$$

  • Quy luật phân bố $N%/D_{1.3}$: Thử nghiệm qua 4 hàm (Weibull, Meyer, Khoảng cách, Chuẩn). Kết quả hàm Weibull đạt độ thích ứng cao nhất với $\chi^2_{\text{tính}} < \chi^2_{0,05; df}$, các tham số xác định: $\lambda = 0,246$; $\alpha = 0,742$. Dạng đường cong phân bố giảm dần (chữ J ngược), chứng tỏ lâm phần có mật độ cây đường kính nhỏ chiếm ưu thế tuyệt đối, khả năng phục hồi tự nhiên dồi dào.
  • Quy luật phân bố $N%/H_{vn}$: Hàm Phân bố Chuẩn (Normal Distribution) đạt độ thích ứng tối ưu với kỳ vọng toán $\mu = 12,693$ và phương sai $\sigma^2 = 18,738$. Đường cong phân bố có dạng đỉnh hơi lệch phải, phản ánh tầng tán rừng đang tập trung phát triển mạnh ở chiều cao 10–20 m.

Kết quả định lượng chi tiết

Cấu trúc tổ thành loài tầng cây cao ($IV_i%$)

Qua điều tra 12 OTC, ghi nhận tổng số 865 cá thể thuộc 55 loài và 27 họ thực vật. Họ Dầu (Dipterocarpaceae) chiếm số loài nhiều nhất (8 loài, 74 cây), kế đến là họ Thị (Ebenaceae - 4 loài, 49 cây), họ Côm (Euphorbiaceae - 3 loài, 79 cây).

STT Tên loài cây Số cây ($N$) Tiết diện ngang ($G$, $\text{m}^2$) Tần số ($F%$) $N%$ $G%$ $F%$ $IV_i%$
1 Cậy (Diospyros) 45 3,1924 75,0 5,20 12,73 3,16 7,03
2 Trâm trắng (Syzygium) 66 1,7837 100,0 7,63 7,11 4,21 6,32
3 Bứa (Garcinia) 74 1,4394 91,7 8,55 5,74 3,86 6,05
4 Dâu rừng (Baccaurea) 51 1,6818 100,0 5,90 6,71 4,21 5,60
-- Tổng 4 loài ưu thế 236 8,0973 -- 27,28 32,29 15,44 25,00
-- 51 loài còn lại 629 16,9771 -- 72,72 67,71 84,56 75,00
Tổng 55 loài 865 25,0744 -- 100,0 100,0 100,0 100,0 
TỶ LỆ TỔ THÀNH LOÀI ƯU THẾ (IVi%)
+-----------------------------------------------------------------------------+
| [Cậy: 7.03%] [Trâm trắng: 6.32%] [Bứa: 6.05%] [Dâu rừng: 5.60%]             |
| ===> 4 loài ưu thế chiếm 25.00% IVi%                                        |
| [51 loài khác đóng góp phần còn lại: 75.00% IVi%]                           |
+-----------------------------------------------------------------------------+

Đặc trưng lâm học và phân bố cấu trúc lâm phần

  • Đường kính bình quân ($D_{1.3}$): $16,3 \pm 0,68\text{ cm}$ (Độ lệch chuẩn $S = 10,26\text{ cm}$; biên độ $6,40 - 90,40\text{ cm}$).
  • Chiều cao vút ngọn bình quân ($H_{vn}$): $13,36 \pm 0,25\text{ m}$ (Độ lệch chuẩn $S = 3,81\text{ m}$; biên độ $5,0 - 26,0\text{ m}$).
  • Mật độ bình quân ($N$): $721\text{ cây/ha}$.
  • Tổng tiết diện ngang ($G$): $20,90\text{ m}^2/\text{ha}$.
  • Trữ lượng thân cây đứng ($M$): $161,06\text{ m}^3/\text{ha}$ (Đối chiếu Thông tư 33/2018/TT-BNNPTNT: Xác định chính xác trạng thái rừng tự nhiên lá rộng thường xanh trung bình).
  • Hệ số hỗn giao ($K$): $K = 0,33$ ($K < 0,5$, phản ánh độ hỗn giao ở mức thấp, số cá thể tập trung vào một số nhóm loài cục bộ).

Phân bố mật độ, tiết diện ngang và trữ lượng theo nhóm đường kính ($D_{1.3}$)

Nhóm $D_{1.3}$ (cm) Mật độ ($N$, cây/ha) Tỷ lệ $N$ (%) Tiết diện ngang ($G$, $\text{m}^2/\text{ha}$) Tỷ lệ $G$ (%) Trữ lượng ($M$, $\text{m}^3/\text{ha}$) Tỷ lệ $M$ (%)
$< 20$ 539 74,8 6,23 29,8 36,69 22,8
$20 - 40$ 155 21,5 8,78 42,0 68,13 42,3
$40 - 60$ 23 3,2 4,30 20,6 39,40 24,5
$> 60$ 3 0,5 1,58 7,6 16,84 10,5
Tổng cộng 721 100,0 20,90 100,0 161,06 100,0 
PHÂN BỐ TRỮ LƯỢNG GỖ THEO CẤP ĐƯỜNG KÍNH (m3/ha)
D1.3 < 20 cm   : [====================] 36.69 m3/ha (22.8%)
D1.3 20 - 40 cm: [==========================================] 68.13 m3/ha (42.3%)
D1.3 40 - 60 cm: [========================] 39.40 m3/ha (24.5%)
D1.3 > 60 cm   : [==========] 16.84 m3/ha (10.5%)

Phân bố theo lớp chiều cao ($H_{vn}$)

  • Lớp $H_{vn} < 10\text{ m}$: 173 cây/ha (23,9%), $G = 1,39\text{ m}^2/\text{ha}$, $M = 5,42\text{ m}^3/\text{ha}$.
  • Lớp $H_{vn} = 10 - 15\text{ m}$: 324 cây/ha (45,0%), $G = 5,58\text{ m}^2/\text{ha}$, $M = 33,59\text{ m}^3/\text{ha}$.
  • Lớp $H_{vn} = 15 - 20\text{ m}$: 198 cây/ha (27,4%), $G = 9,20\text{ m}^2/\text{ha}$, $M = 73,13\text{ m}^3/\text{ha}$ (chiếm tỷ trọng trữ lượng cao nhất: $45,4%$).
  • Lớp $H_{vn} > 20\text{ m}$: 27 cây/ha (3,7%), $G = 4,75\text{ m}^2/\text{ha}$, $M = 48,92\text{ m}^3/\text{ha}$.

Chỉ số đa dạng sinh học quần xã thực vật

  • Chỉ số phong phú loài Margalef ($d$): $5,29 \pm 0,70$ $\rightarrow$ Độ phong phú loài đạt mức trung bình khá.
  • Chỉ số đa dạng Shannon-Weiner ($H'$): $2,91 \pm 0,14$ $\rightarrow$ Thuộc phân cấp đa dạng cao ($2 < H' < 4$).
  • Chỉ số đồng đều Pielou ($J'$): $0,93 \pm 0,02$ $\rightarrow$ Mức độ đồng đều rất cao (tiệm cận 1,0), các cá thể phân bổ tương đối đồng nhất giữa các loài.
  • Chỉ số ưu thế Simpson ($\lambda'$): $0,05 \pm 0,01$ $\rightarrow$ Giá trị rất thấp, phản ánh không có hiện tượng một loài đơn lẻ độc tôn lấn át sinh cảnh.

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

Cải tiến kỹ thuật và tính mới khoa học

  1. Xác lập bộ tham số Weibull chuẩn hóa: Lần đầu tiên tham số hình dạng $\alpha = 0,742$ và quy mô $\lambda = 0,246$ được xác lập chính xác cho rừng lá rộng thường xanh trung bình tại VQG Lò Gò - Xa Mát, tạo chuẩn đối sánh cho vùng sinh thái Đông Nam Bộ.
  2. Lượng hóa chi tiết cấu trúc đa tầng: Tích hợp đồng thời hai biến số $D_{1.3}$ và $H_{vn}$ với hệ số phân bố thực nghiệm, chứng minh tỷ lệ cây tái sinh tiềm năng ($D_{1.3} < 20\text{ cm}$) đạt $74,8%$ nhưng trữ lượng lại tập trung chủ yếu ở tầng vượt tán và tán chính ($H_{vn} \ge 15\text{ m}$ chiếm $75,8%$ tổng trữ lượng).
+------------------------------------------------------------------------------------------+
|                  SO SÁNH KẾT QUẢ VỚI CÁC CÔNG TRÌNH NGHIÊN CỨU TIÊU BIỂU                 |
+--------------------------+-----------------------+-----------------------+---------------+
| Tác giả / Địa điểm       | Trữ lượng (M, m3/ha)  | Chỉ số Shannon (H')   | Độ hỗn giao K |
+--------------------------+-----------------------+-----------------------+---------------+
| Đề tài (VQG Lò Gò)       | 161,06                | 2,91                  | 0,33          |
| Trần Văn Con (Rừng khộp) | 98,50 - 120,00        | 1,85 - 2,10           | 0,21          |
| Viên Ngọc Nam (Phước Bình)| 185,40               | 3,12                  | 0,42          |
| Đặng Văn Thuyết (Bảo Lâm)| 172,30                | 2,75                  | 0,36          |
+--------------------------+-----------------------+-----------------------+---------------+

Đóng góp thực tiễn cho ngành lâm nghiệp

  • Cung cấp dữ liệu gốc chuẩn xác phục vụ lập Phương án quản lý rừng bền vững giai đoạn 2021–2030 của VQG Lò Gò - Xa Mát.
  • Chỉ ra 4 loài ưu thế sinh thái (Cậy, Trâm trắng, Bứa, Dâu rừng) làm cơ sở lựa chọn cây mẹ lấy giống phục vụ gieo ươm và làm giàu rừng bản địa.

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

Kịch bản ứng dụng lâm sinh thực địa

  • Phân vùng bảo tồn nghiêm ngặt (10.615,29 ha): Giữ nguyên vẹn cấu trúc tự nhiên, nghiêm cấm mọi tác động cơ học. Thiết lập 12 OTC thành các ô định vị cố định (Permanent Sampling Plots) để quan trắc động thái sinh trưởng định kỳ 5 năm/lần.
  • Phân khu phục hồi sinh thái (19.256,78 ha): Tận dụng mật độ cây con tự nhiên dồi dào ($539\text{ cây/ha}$ nhóm $D_{1.3} < 20\text{ cm}$), áp dụng biện pháp khoanh nuôi xúc tiến tái sinh tự nhiên có trồng bổ sung các loài họ Dầu (Dipterocarpus alatus, Hopea odorata, Shorea roxburghii) nhằm nâng cao hệ số hỗn giao từ $K = 0,33$ lên ngưỡng tối ưu $K \ge 0,50$.
+------------------------------------------------------------------------------------+
|                         LỘ TRÌNH TRIỂN KHAI ỨNG DỤNG LÂM SINH                      |
+------------------------------------------------------------------------------------+
| [2023 - 2024] -> Số hóa 12 OTC lên WebGIS & Bàn giao dữ liệu cho Ban Quản lý VQG  |
| [2025 - 2027] -> Triển khai nuôi dưỡng, tỉa thưa vệ sinh ở phân khu phục hồi       |
| [2028 - 2030] -> Tái điều tra chu kỳ II, kiểm định tăng trưởng trữ lượng & delta H'|
+------------------------------------------------------------------------------------+

Hiệu quả kinh tế - sinh thái (ROI & Ecological Benefits)

  • Tiết kiệm ngân sách phục hồi rừng: Việc xác định được trữ lượng cây tái sinh tự nhiên dồi dào giúp giảm $40 - 50%$ chi phí đầu tư trồng mới đắt đỏ (tiết kiệm khoảng 15–20 triệu VNĐ/ha so với trồng rừng mới hoàn toàn).
  • Bảo vệ nguồn nước và chống cháy rừng: Tầng tán $H_{vn} = 10 - 20\text{ m}$ khép tán tốt giúp giữ ẩm độ đất xám phù sa cổ, hạn chế tối đa nguy cơ cháy rừng lan tràn trong mùa khô.

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

Hạn chế kỹ thuật

  • Hệ thống 12 OTC là ô tiêu chuẩn tạm thời, chưa được đóng mốc bê tông cố định góc ranh giới ngoài thực địa để theo dõi sinh trưởng cá thể dài hạn.
  • Chưa điều tra chi tiết tầng cây tái sinh có chiều cao $H < 1,3\text{ m}$ và đường kính $D_{1.3} < 6\text{ cm}$ để đánh giá trọn vẹn chỉ số triển vọng phục hồi.

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

  1. Ứng dụng công nghệ LiDAR viễn thám: Kết hợp bay quét UAV LiDAR với dữ liệu 12 OTC thực địa để xây dựng bản đồ số 3D sinh khối và trữ lượng carbon toàn bộ diện tích 29.997,87 ha.
  2. Nghiên cứu động thái rễ và đất: Đánh giá mối tương quan giữa độ dày tầng đất xám loang lổ đỏ vàng với hệ số sinh trưởng của các loài cây họ Dầu.

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

  • Sinh viên và Học viên Cao học Lâm nghiệp: Tài liệu tham khảo chuẩn mực về phương pháp kết hợp sinh trắc thống kê lâm nghiệp với phần mềm chuyên dụng (Statgraphics, Biodiversity Pro).
  • Ban Quản lý Vườn Quốc gia & Cán bộ Kiểm lâm: Cung cấp cơ sở khoa học định lượng để lập hồ sơ theo dõi diễn biến rừng, thiết kế các đợt tuần tra PCCCR tại các tiểu khu trọng điểm.
  • Các Nhà Quy hoạch Lâm nghiệp & Doanh nghiệp Tín chỉ Carbon: Dữ liệu trữ lượng ($161,06\text{ m}^3/\text{ha}$) và tiết diện ngang ($20,90\text{ m}^2/\text{ha}$) là đầu vào chính xác để tính toán sinh khối (Biomass) và quy đổi tín chỉ carbon rừng theo chuẩn quốc tế.
  • Cộng đồng Nghiên cứu Sinh thái học: Bổ sung dữ liệu đa dạng sinh học của một vùng sinh thái chuyển tiếp đặc biệt tại Đông Nam Á.

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

1. Cần những trang thiết bị và kỹ năng gì để triển khai phương pháp điều tra này?

Cần trang bị máy định vị vệ tinh (GPS cầm tay sai số $< 3\text{ m}$), thước dây chuyên dụng đo chu vi/đường kính thân cây, thước đo cao Blume-Leiss, la bàn địa chất và phần mềm Statgraphics Centurion/Biodiversity Pro. Nhân lực điều tra cần kỹ năng nhận diện thực vật thực địa và xử lý thống kê sinh học.

2. Tại sao hàm Weibull lại được lựa chọn cho phân bố đường kính thay vì hàm Meyer?

Mặc dù hàm Meyer thường dùng cho rừng tự nhiên nhiệt đới, kết quả kiểm định Chi-square ($\chi^2$) thực tế trên 865 cây tại VQG Lò Gò - Xa Mát cho thấy hàm Weibull đạt mức ý nghĩa thống kê $P > 0,05$ cao hơn, mô tả chính xác hơn dạng đường cong phân bố giảm ở các cấp kính nhỏ ($D_{1.3} < 20\text{ cm}$).

3. Kết quả nghiên cứu có thể tích hợp vào hệ thống WebGIS giám sát rừng như thế nào?

Toàn bộ tọa độ 12 OTC và thuộc tính cây gỗ ($D_{1.3}, H_{vn}, V, IV_i%$) đã được số hóa trên MapInfo/Global Mapper dưới định dạng Shapefile/GeoJSON, sẵn sàng tích hợp trực tiếp vào cơ sở dữ liệu không gian của hệ thống WebGIS giám sát tài nguyên rừng VQG.

4. Hệ số hỗn giao $K = 0,33$ phản ánh thực trạng gì và cần tác động lâm sinh ra sao?

Hệ số $K = 0,33$ thuộc cấp hỗn giao thấp ($K < 0,5$), chứng tỏ số loài phân bố không đều và có sự tập trung cá thể cao ở 4 loài ưu thế. Cần áp dụng biện pháp xúc tiến tái sinh có trồng bổ sung các loài bản địa quý hiếm khác thuộc họ Dầu và họ Đậu để làm giàu thành phần loài.

5. Chi phí triển khai hệ thống ô đo đếm này trên quy mô toàn lâm phần là bao nhiêu?

Chi phí điều tra thực địa theo định mức kinh tế kỹ thuật hiện hành dao động từ 1,8 – 2,5 triệu VNĐ/OTC $1.000\text{ m}^2$. Với quy mô 30–50 OTC định vị cố định cho toàn bộ phân khu phục hồi sinh thái, tổng chi phí ước tính khoảng 80–120 triệu VNĐ với chu kỳ tái đo 5 năm/lần.

Kết luận

Nghiên cứu đã giải quyết trọn vẹn bài toán định lượng cấu trúc và đa dạng sinh thái tại VQG Lò Gò - Xa Mát với các đóng góp cốt lõi:

  1. Xác định chính xác trạng thái rừng: Trữ lượng bình quân đạt $161,06\text{ m}^3/\text{ha}$, tiết diện ngang $20,90\text{ m}^2/\text{ha}$, mật độ 721 cây/ha, khẳng định đối tượng thuộc trạng thái rừng tự nhiên lá rộng thường xanh trung bình theo Thông tư 33/2018/TT-BNNPTNT.
  2. Mô hình hóa quy luật cấu trúc: Phân bố $N%/D_{1.3}$ tuân theo hàm Weibull ($\lambda = 0,246; \alpha = 0,742$) và $N%/H_{vn}$ tuân theo hàm Phân bố Chuẩn ($\mu = 12,693; \sigma^2 = 18,738$).
  3. Đánh giá đa dạng sinh học: Quần xã đạt độ đa dạng cao ($H' = 2,91 \pm 0,14$), độ đồng đều lớn ($J' = 0,93 \pm 0,02$) với 4 loài ưu thế sinh thái gồm Cậy, Trâm trắng, Bứa, Dâu rừng ($IV_i% = 25,00%$).

Kết quả này là luận cứ khoa học thực tiễn vững chắc cho Ban Quản lý Vườn Quốc gia Lò Gò - Xa Mát trong việc triển khai các giải pháp lâm sinh thích ứng, bảo tồn đa dạng sinh học và phát triển tài nguyên rừng bền vững trong tương lai.