Giới thiệu dự án
Bối cảnh và thực trạng ngành lâm nghiệp
Theo số liệu từ Tổ chức Lương thực và Nông nghiệp Liên Hợp Quốc (FAO), trung bình mỗi năm thế giới mất khoảng 13 triệu ha rừng tự nhiên, làm suy giảm nghiêm trọng khả năng hấp thụ $\text{CO}_2$ và chiếm từ 12% đến 20% tổng lượng phát thải khí nhà kính toàn cầu. Tại Việt Nam, độ che phủ rừng nguyên sinh đã sụt giảm mạnh từ mức 70% ở đầu thế kỷ XX xuống còn 43% vào giữa thế kỷ và chỉ còn 24% trong giai đoạn 1979 - 1981. Dù diện tích rừng trồng đã tăng lên trong những năm gần đây, nhưng chất lượng rừng tự nhiên — đặc biệt là kiểu thảm thực vật rừng gỗ lá rộng thường xanh nhiệt đới — đang bị suy thoái nghiêm trọng, chuyển dịch phần lớn sang trạng thái rừng nghèo và rừng nghèo kiệt.
Vườn Quốc gia (VQG) Lò Gò - Xa Mát (tỉnh Tây Ninh) có tổng diện tích quản lý hơn 30.000 ha (chiếm 26% diện tích tự nhiên của tỉnh), giữ vai trò là "lá phổi xanh" và hành lang sinh thái kết nối vùng rừng ngập lũ hạ lưu sông Mê Kông với vùng bán bình nguyên Đông Nam Bộ. Tuy nhiên, dưới tác động của lịch sử khai thác trước khi thành lập khu bảo tồn (thời kỳ Lâm trường Hoà Hiệp và Tân Bình) cùng diễn thế tự nhiên, diện tích rừng gỗ lá rộng thường xanh tại Tiểu khu 21 đã bị suy giảm trữ lượng nghiêm trọng, rơi vào trạng thái nghèo kiệt, đặt ra yêu cầu cấp bách về một giải pháp nghiên cứu định lượng cấu trúc lâm học phục vụ bảo tồn và phục hồi sinh thái.
+-------------------------------------------------------------------------------+
| KHUNG ĐÁNH GIÁ CẤU TRÚC LÂM HỌC |
+-------------------------------------------------------------------------------+
| ĐIỀU TRA NGOẠI NGHIỆP XỬ LÝ DỮ LIỆU ĐỊNH LƯỢNG ỨNG DỤNG QUẢN LÝ |
| - 5 OTC (1.000 m2/ô) -> - Cấu trúc N/D, N/H -> - Phục hồi rừng |
| - 25 ÔDB (25 m2/ô) - Hệ số tổ thành IVI% - Bảo tồn loài nguy|
| - GPS Locus Map UTM - Đa dạng Shannon/Pielou cấp IUCN |
+-------------------------------------------------------------------------------+
Vấn đề cốt lõi (Problem Statement)
Công tác quản trị và phục hồi rừng tại các phân khu phục hồi sinh thái hiện nay đang đối mặt với các nút thắt kỹ thuật lớn:
- Thiếu cơ sở dữ liệu định lượng chính xác: Việc phân chia trạng thái rừng theo các phương pháp truyền thống mang tính định tính cao, thiếu các thông số cấu trúc toán học cụ thể ($N/D_{1.3}$, $N/H_{vn}$, phân bố không gian).
- Mất cân bằng cấu trúc quần thụ: Sự suy giảm các loài cây gỗ lớn họ Sao Dầu (Dipterocarpaceae) và sự xâm lấn của các loài tiên phong kém giá trị thương mại và sinh thái.
- Tắc nghẽn chu trình tái sinh tự nhiên: Thiếu dữ liệu đánh giá mối tương quan giữa tầng cây cao ($D_{1.3} \ge 8\text{ cm}$) và tầng cây tái sinh non ($H < 1\text{ m}$ đến $H > 4\text{ m}$) để đưa ra biện pháp nuôi dưỡng, xúc tiến tái sinh tự nhiên phù hợp.
Mục tiêu nghiên cứu cụ thể
- Định danh và thống kê hệ thực vật thân gỗ: Xác định toàn bộ thành phần loài, chi, họ thực vật thân gỗ ($D_{1.3} \ge 8\text{ cm}$) tại Tiểu khu 21.
- Định lượng cấu trúc lâm phần: Xác định mật độ ($N$), tiết diện ngang ($G$), trữ lượng ($M$), độ hỗn giao ($K$), chỉ số quan trọng ($IVI%$), phân bố $N/D_{1.3}$ và $N/H_{vn}$.
- Đánh giá động thái tái sinh tự nhiên: Phân tích mật độ, tổ thành loài, cấp chiều cao, chất lượng cây con và nguồn gốc tái sinh (hạt/chồi) dưới tán rừng.
- Phân tích chỉ số đa dạng sinh học và độ hiếm: Định lượng các chỉ số sinh thái Margalef ($d$), Pielou ($J'$), Shannon-Wiener ($H'$), Simpson ($\lambda'$), Caswell ($V$), Chỉ số hiếm ($RI$) và tra cứu Sách Đỏ IUCN (2022).
- Đề xuất giải pháp lâm sinh ứng dụng: Xây dựng phác đồ can thiệp lâm sinh (khoanh nuôi, xúc tiến tái sinh, làm giàu rừng) dựa trên dữ liệu thực nghiệm.
Phương pháp tiếp cận giải pháp
Đề tài tích hợp phương pháp điều tra tiêu chuẩn lâm học nhiệt đới kết hợp phần mềm phân tích sinh thái thống kê hiện đại:
- Khảo sát thực địa qua hệ thống ô tiêu chuẩn (OTC) điển hình kết hợp định vị vệ tinh viễn thám GPS/GIS.
- Ứng dụng mô hình hóa phân bố toán học (Weibull, Distance, Exponential) để lượng hóa cấu trúc quần xã thực vật.
- Chuẩn hóa quy trình phân tích đa dạng sinh học thông qua các thuật toán chỉ số sinh thái học quần xã.
Kết quả kỳ vọng với các chỉ số đo lường (Metrics)
- Xác định chính xác mật độ tầng cây đứng dao động trong ngưỡng $600 - 800\text{ cây/ha}$, trữ lượng $10 - 50\text{ m}^3\text{/ha}$ (chuẩn hóa theo Thông tư 33/2018/TT-BNNPTNT).
- Định lượng mật độ tầng cây tái sinh đạt $> 10.000\text{ cây/ha}$ với tỷ lệ tái sinh từ hạt chiếm $> 80%$.
- Khẳng định chỉ số đa dạng Shannon-Wiener ($H'$) đạt mức sinh thái cao ($> 3.0$) và chỉ số tương đồng $J' > 0.9$.
Phạm vi và giới hạn nghiên cứu
- Phạm vi không gian: Khu vực rừng gỗ lá rộng thường xanh nghèo kiệt tại Tiểu khu 21, VQG Lò Gò - Xa Mát, tỉnh Tây Ninh (Tọa độ UTM: $X: 595073 - 595417$, $Y: 1282419 - 1282617$).
- Phạm vi đối tượng: Tầng cây gỗ có đường kính ngang ngực $D_{1.3} \ge 8\text{ cm}$ và tầng cây tái sinh thân gỗ tự nhiên dưới tán rừng.
- Giới hạn kỹ thuật: Nghiên cứu thực hiện trên 5 ô tiêu chuẩn điển hình ($5.000\text{ m}^2$) và 25 ô dạng bản ($625\text{ m}^2$) trong giai đoạn từ tháng 04/2023 đến tháng 08/2023.
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 tiêu chuẩn phân loại trạng thái rừng
| Tiêu chuẩn / Phương pháp |
Cơ sở phân loại |
Ưu điểm |
Nhược điểm áp dụng thực tế |
| Loeschau (1963) |
Giai đoạn phát triển (Non, trung niên, thành thục, quá thành thục) |
Thích hợp rừng thuần loài, đều tuổi |
Khó áp dụng cho rừng tự nhiên nhiệt đới nhiều tầng tán |
| Quy phạm QP6-84 |
Nguồn gốc và trữ lượng (Nhóm I đến IV) |
Đơn giản, dễ nhận biết sơ bộ |
Mang tính định tính, biên độ trữ lượng quá rộng |
| Thông tư 34/2009/TT-BNNPTNT |
Nguồn gốc, loài cây, trữ lượng |
Chi tiết hóa phân cấp trữ lượng |
Một số chỉ tiêu chồng chéo trong điều tra diện rộng |
| Thông tư 33/2018/TT-BNNPTNT (Áp dụng) |
Trữ lượng cây đứng ($M < 10$, $10-50$, $50-100$, $101-200$, $201-300$, $>300\text{ m}^3\text{/ha}$) |
Phân lượng chính xác, định lượng hóa rõ ràng trạng thái nghèo kiệt ($10-50\text{ m}^3\text{/ha}$) |
Yêu cầu số liệu đo đếm ngoại nghiệp nghiêm ngặt ($D_{1.3}$, $H_{vn}$) |
Phân tích yêu cầu lâm học theo mô hình MoSCoW
+-------------------------------------------------------------------------------+
| MÔ HÌNH YÊU CẦU LÂM HỌC (MoSCoW) |
+-------------------------------------------------------------------------------+
| MUST HAVE (Bắt buộc) | SHOULD HAVE (Nên có) |
| - Đo đếm D1.3 >= 8cm, Hvn | - Phân tích tương quan D1.3 và Hvn |
| - Tính toán IVI% tầng cây gỗ | - Phân nhóm đường kính (D < 10cm, D >= 10cm)|
| - Đo đếm cây tái sinh 5 cấp | - Kiểm định Caswell V-test |
|--------------------------------+----------------------------------------------|
| COULD HAVE (Có thể có) | WON'T HAVE (Chưa áp dụng đợt này) |
| - Tra cứu Sách Đỏ IUCN (2022) | - Lập biểu thể tích cá thể chuyên biệt |
| - Tính chỉ số hiếm RI | - Đo sinh khối rễ dưới mặt đất |
+-------------------------------------------------------------------------------+
Thiết kế hệ thống khảo sát và xử lý dữ liệu
Cấu trúc hệ thống đo đếm và kiến trúc phần mềm
Quy trình điều tra và xử lý dữ liệu được thiết kế theo cấu trúc module hóa 3 lớp:
[ LỚP NGOẠI NGHIỆP: GPS Locus Map Free + Thước đo bách phân + Thước đo D13 ]
[ LỚP TIỀN XỬ LÝ: MS Excel 2016 (Data Analysis ToolPak) - Chuẩn hóa ma trận ]
[ LỚP PHÂN TÍCH ĐỊNH LƯỢNG: Statgraphics Centurion 15.1 + PRIMER v6.1.13 + R ]
Công nghệ và công cụ sử dụng (Tech Stack)
- Phần mềm điều hướng ngoại nghiệp: Locus Map Free v4.x (chạy trên Android OS với chip GPS đa băng tần, sai số $< 3\text{ m}$, hệ tọa độ UTM WGS84).
- Phần mềm xử lý thống kê sinh học:
- Microsoft Excel 2016 (Data Analysis ToolPak): Xử lý cơ sở dữ liệu sơ cấp, thống kê mô tả.
- Statgraphics Centurion 15.1: Kiểm định phân phối chuẩn, hồi quy phi tuyến $N/D_{1.3}$, $N/H_{vn}$.
- PRIMER v6 (Plymouth Routines in Multivariate Ecological Research v6.1.13): Phân tích đa dạng sinh thái học (Shannon, Margalef, Simpson, Caswell).
- Công cụ tra cứu phân loại học: "Sách tra cứu tên cây cỏ Việt Nam" (Võ Văn Chi, 2007) và Cơ sở dữ liệu danh lục đỏ IUCN Red List (2022-2).
Phương pháp nghiên cứu (Methodology)
[BƯỚC 1: LẬP Ô TIÊU CHUẨN]
[BƯỚC 2: THU THẬP SỐ LIỆU TẦNG CÂY GỖ LỚN]
[BƯỚC 3: THU THẬP SỐ LIỆU TẦNG CÂY TÁI SINH]
[BƯỚC 4: TÍNH TOÁN CÁC CHỈ SỐ LÂM HỌC VÀ ĐA DẠNG]
Ma trận đánh giá rủi ro và giải pháp khắc phục
| Rủi ro kỹ thuật |
Mức độ |
Biện pháp giảm thiểu (Mitigation Strategy) |
| Sai số định vị dưới tán rừng rậm |
Trung bình |
Sử dụng GPS tích hợp A-GPS và đối chiếu mốc thực địa vật lý kết hợp bản đồ vệ tinh Google Earth độ phân giải cao. |
| Nhầm lẫn phân loại loài thực vật |
Cao |
Thu thập mẫu tiêu bản lá, hoa, vỏ thân; đánh mã sp1, sp2 và tham vấn chuyên gia thực vật học đối chiếu tiêu bản chuẩn. |
| Biến động cục bộ do vi địa hình |
Thấp |
Phân bổ 5 OTC cách biệt không gian đại diện cho toàn bộ tiểu khu 21 (khoảng cách giữa các ô từ $150\text{ m} - 500\text{ m}$). |
Triển khai thực nghiệm và kết quả
Thuật toán và công thức tính toán lâm học
Dữ liệu đo đếm từ 5 OTC được xử lý bằng thuật toán tự động hóa viết trên nền tảng ngôn ngữ R/Python để chuẩn hóa các chỉ tiêu lâm học:
import numpy as np
import pandas as pd
def calculate_forestry_metrics(df_trees, total_area_ha=0.5):
"""
df_trees DataFrame chứa các cột: 'species', 'd13_cm', 'h_m', 'v_m3'
"""
total_trees = len(df_trees)
density_ha = total_trees / total_area_ha
# Tính tiết diện ngang G (m2) = pi * (D13 / 200)^2
df_trees['g_m2'] = np.pi * (df_trees['d13_cm'] / 200) ** 2
total_g_ha = df_trees['g_m2'].sum() / total_area_ha
total_v_ha = df_trees['v_m3'].sum() / total_area_ha
# Tính chỉ số tổ thành IVI% theo phương pháp Thái Văn Trừng (1999)
species_grp = df_trees.groupby('species').agg(
n_i=('d13_cm', 'count'),
g_i=('g_m2', 'sum'),
v_i=('v_m3', 'sum')
).reset_index()
species_grp['Ni_pct'] = (species_grp['n_i'] / total_trees) * 100
species_grp['Gi_pct'] = (species_grp['g_i'] / df_trees['g_m2'].sum()) * 100
species_grp['Vi_pct'] = (species_grp['v_i'] / df_trees['v_m3'].sum()) * 100
species_grp['IVI_pct'] = (species_grp['Ni_pct'] + species_grp['Gi_pct'] + species_grp['Vi_pct']) / 3
# Chỉ số đa dạng Shannon-Wiener (H') và Pielou (J')
pi = species_grp['n_i'] / total_trees
h_prime = -np.sum(pi * np.log(pi))
s_species = len(species_grp)
j_prime = h_prime / np.log(s_species)
margalef_d = (s_species - 1) / np.log(total_trees)
simpson_lambda = np.sum((species_grp['n_i'] * (species_grp['n_i'] - 1)) / (total_trees * (total_trees - 1)))
return {
'density_ha': density_ha,
'G_ha': total_g_ha,
'V_ha': total_v_ha,
'H_prime': h_prime,
'J_prime': j_prime,
'Margalef_d': margalef_d,
'Simpson_lambda': simpson_lambda,
'IVI_table': species_grp.sort_values(by='IVI_pct', ascending=False)
}
Kết quả cấu trúc tầng cây gỗ lớn ($D_{1.3} \ge 8\text{ cm}$)
Tổng hợp các chỉ số cơ bản qua 5 ô tiêu chuẩn
Kết quả điều tra trên diện tích mẫu $5.000\text{ m}^2$ ($0,5\text{ ha}$) ghi nhận tổng cộng 37 loài thuộc 23 họ thực vật thân gỗ:
| Ký hiệu OTC |
Số loài ($S$) (Loài ưu thế) |
Mật độ $N$ (cây/ha) |
Tiết diện ngang $G$ ($\text{m}^2\text{/ha}$) |
Trữ lượng $M$ ($\text{m}^3\text{/ha}$) |
Ưu hợp thực vật xác lập |
| OTC 1 |
31 (5) |
760 |
7,1 |
32,1 |
Dầu song nàng + Dầu rái + Sầm lá lớn |
| OTC 2 |
33 (7) |
800 |
7,3 |
31,7 |
Dầu rái + Bằng lăng ổi + Dẻ trắng |
| OTC 3 |
32 (6) |
690 |
6,9 |
31,9 |
Bình linh lông + Sầm lá lớn + Dầu song nàng |
| OTC 4 |
30 (4) |
800 |
7,6 |
35,2 |
Dầu song nàng + Dầu rái + Sầm lá lớn |
| OTC 5 |
29 (3) |
640 |
5,7 |
26,0 |
Cù đèn + Dầu rái + Dầu song nàng |
| Giá trị TB |
37 (3) |
738 |
6,9 |
31,4 |
Dầu song nàng + Dầu rái + Sầm lá lớn |
TỶ LỆ TỔ THÀNH TẦNG CÂY GỖ TOÀN KHU VỰC (IVI%):
Tổng cộng: 100% (3 loài ưu thế chiếm 21.9% tổng IVI toàn lâm phần)
Đặc trưng phân bố cấu trúc không gian ($N/D_{1.3}$ và $N/H_{vn}$)
| Chỉ số thống kê sinh học |
Đường kính $D_{1.3}$ (cm) |
Chiều cao vút ngọn $H_{vn}$ (m) |
| Giá trị trung bình ($\bar{X}$) |
10,7 |
10,1 |
| Sai số chuẩn ($SE$) |
0,11 |
0,10 |
| Độ lệch chuẩn ($SD$) |
2,04 |
2,01 |
| Phương sai ($S^2$) |
4,20 |
4,00 |
| Hệ số biến động ($C_v%$) |
19,0% |
19,9% |
| Hệ số độ lệch ($Sk$) |
2,00 (Lệch trái mạnh) |
0,02 (Đối xứng quanh trung bình) |
| Hệ số độ nhọn ($Ku / Ex$) |
7,70 (Đỉnh nhọn cao) |
-1,20 (Dạng bẹt phân tán đều) |
| Biên độ biến động ($R$) |
13,4 (8,0 - 21,4 cm) |
6,0 (7,0 - 13,0 m) |
- Quy luật phân bố $N/D_{1.3}$: Biểu diễn dạng đường cong giảm sâu (hàm số mũ). Cấp kính nhỏ $D \approx 10\text{ cm}$ (khoảng 8 - 12 cm) tập trung tới 578 cây/ha (chiếm 78,3% mật độ). Cấp đường kính lớn ($18 - 22\text{ cm}$) chỉ có 18 cây/ha (chiếm 2,5%), phản ánh tình trạng rừng sau khai thác chọn kiệt quệ trong quá khứ.
- Quy luật phân bố $N/H_{vn}$: Phân bố đối xứng dạng hình chuông bẹt, số cây phân bổ đồng đều ở các cấp chiều cao: Cấp 9m (29,8%), cấp 11m (29,5%), cấp 13m (29,0%). Chiều cao tán rừng bị giới hạn dưới 14m do đặc tính thổ nhưỡng đất xám bạc màu trên phù sa cổ.
PHÂN BỐ TẦNG KÍNH (N/D13) PHÂN BỐ CHIỀU CAO (N/Hvn)
Cấp D (cm) % Số cây Cấp H (m) % Số cây
10 cm [████████████] 78.3% 7 m [██] 11.7%
14 cm [███] 19.2% 9 m [█████] 29.8%
18 cm [▍] 1.4% 11 m [█████] 29.5%
22 cm [▎] 1.1% 13 m [█████] 29.0%
Kết quả điều tra động thái tầng cây tái sinh
Điều tra trên 25 ô dạng bản ($625\text{ m}^2$) ghi nhận những tín hiệu phục hồi sinh thái vượt bậc:
- Mật độ tái sinh bình quân: Đạt $12.027\text{ cây/ha}$ (thuộc 24 loài, 17 họ thực vật).
- Nguồn gốc tái sinh: Tái sinh từ hạt đạt 87,8%, tái sinh từ chồi chỉ chiếm 12,2%. Đây là minh chứng cho tiềm năng di truyền và khả năng phục hồi tự nhiên bền vững.
- Phẩm chất cây con: Cây chất lượng tốt (Cấp A) chiếm 43,9%, cây trung bình (Cấp B) chiếm 38,2%, cây xấu (Cấp C) chiếm 17,9%.
- Phân bố theo cấp chiều cao:
- Cấp 1 ($H < 1\text{ m}$): $3.456\text{ cây/ha}$ (28,7%)
- Cấp 2 ($1\text{ m} \le H < 2\text{ m}$): $3.744\text{ cây/ha}$ (31,1%)
- Cấp 3 ($2\text{ m} \le H < 3\text{ m}$): $2.528\text{ cây/ha}$ (21,0%)
- Cấp 4 ($3\text{ m} \le H < 4\text{ m}$): $1.344\text{ cây/ha}$ (11,2%)
- Cấp 5 ($H \ge 4\text{ m}$): $955\text{ cây/ha}$ (8,0%)
(Số lượng cây non tập trung chủ yếu ở tầng $< 3\text{ m}$, chiếm 80,8%, phản ánh giai đoạn bùng nổ sinh khối của lớp cây mạ và cây con).
Đánh giá đa dạng sinh học và các loài nguy cấp
- Chỉ số phong phú loài Margalef ($d$): $6,98$ (thể hiện mức độ giàu có về thành phần loài cao).
- Chỉ số đồng đều Pielou ($J'$): $0,94$ (tiệm cận mức 1,0, các loài phân bố rất đồng đều, không bị áp đảo hoàn toàn bởi một loài đơn độc).
- Chỉ số Shannon-Wiener ($H'$): $3,24$ (phản ánh hệ sinh thái có tính ổn định sinh học cao).
- Chỉ số ưu thế Simpson ($\lambda'$): $0,03$ (chỉ số ưu thế cực thấp, khẳng định tính chất đa dạng cao của rừng nhiệt đới).
- Chỉ số Caswell ($V$): Nằm trong khoảng biến thiên $[-2; 2]$, xác nhận hệ sinh thái đang trong trạng thái cân bằng bình thường, không chịu xáo trộn thảm khốc.
- Danh lục loài quý hiếm bảo tồn: Ghi nhận 9 loài nằm trong Danh lục Đỏ IUCN (2022) và Sách Đỏ Việt Nam:
- Pterocarpus pedatus (Giáng hương) - Nguy cấp (EN)
- Dipterocarpus alatus (Dầu rái) - Sắp nguy cấp (VU)
- Dipterocarpus dyeri (Dầu song nàng) - Nguy cấp (EN)
- Hopea odorata (Sao đen) - Sắp nguy cấp (VU)
- Shorea roxburghii (Sến cát) - Sắp nguy cấp (VU)
- Xylia xylocarpa (Căm xe) - Ít lo ngại (LC/Bảo tồn nguồn gen)
- Madhuca cochinchinensis (Sến nam) - Sắp bị đe dọa (NT)
- Canarium album (Trám trắng) - Cần bảo tồn
- Vitex pinnata (Bình linh lông) - Cần bảo tồn
Đổi mới và đóng góp khoa học
Những điểm cải tiến kỹ thuật nổi bật
- Lượng hóa đa chiều chỉ số ưu thế: Không chỉ dừng lại ở số lượng cá thể ($N%$), nghiên cứu tích hợp chuẩn xác 3 chiều $N%$, $G%$, $V%$ theo công thức Thái Văn Trừng để tạo lập chỉ số $IVI%$, phản ánh chân thực vị thế sinh thái của loài.
- Kiểm định môi trường bằng chỉ số Caswell ($V$): Ứng dụng mô hình thống kê lý thuyết $V = \frac{H' - E(H')}{SD(H')}$ để đánh giá mức độ ổn định của sinh cảnh rừng nghèo kiệt, chứng minh rừng nghèo kiệt tại VQG Lò Gò - Xa Mát không phải là rừng thoái hóa hoàn toàn mà đang trong giai đoạn tích lũy nội tại bền vững.
So sánh với các công trình nghiên cứu tương đương
| Chỉ tiêu lâm học |
Nghiên cứu tại Tiểu khu 21, Lò Gò - Xa Mát (2023) |
Nghiên cứu Bình Châu - Phước Bửu (Phan Minh Xuân, 2019) |
Nghiên cứu VQG Bidoup - Núi Bà (Nguyễn Trọng Bình, 2014) |
| Kiểu rừng |
Rừng thường xanh nghèo kiệt |
Rừng nghèo kiệt ven biển |
Rừng kín thường xanh á nhiệt đới |
| Mật độ tầng cây cao ($N$) |
738 cây/ha |
760 cây/ha |
916 cây/ha |
| Đường kính TB ($D_{1.3}$) |
10,7 cm |
12,4 cm |
24,8 cm |
| Chiều cao TB ($H_{vn}$) |
10,1 m |
8,5 m |
13,6 m |
| Mật độ tái sinh |
12.027 cây/ha |
8.450 cây/ha |
6.200 cây/ha |
| Tái sinh từ hạt |
87,8% |
78,5% |
65,0% |
| Chỉ số Shannon ($H'$) |
3,24 |
2,80 |
3,62 |
- Đánh giá vượt trội: Nghiên cứu ghi nhận mật độ tái sinh ($12.027\text{ cây/ha}$) và tỷ lệ tái sinh từ hạt ($87,8%$) cao hơn đáng kể so với khu vực rừng bán khô hạn Bình Châu - Phước Bửu (tăng 42,3% về mật độ tái sinh), mở ra tiềm năng phục hồi tự nhiên mà không cần can thiệp trồng mới tốn kém.
Ứng dụng thực tế và triển khai
Kịch bản ứng dụng lâm sinh thực tế
- Khoanh nuôi xúc tiến tái sinh tự nhiên không trồng bổ sung: Với mật độ cây tái sinh đạt $12.027\text{ cây/ha}$ và $87,8%$ từ hạt, biện pháp kỹ thuật tối ưu là bảo vệ nghiêm ngặt, phòng chống cháy rừng trong mùa khô (tháng 12 đến tháng 4), tỉa bớt dây leo và phát quang cây bụi cạnh tranh ánh sáng quanh các cá thể Dầu song nàng, Dầu rái con.
- Làm giàu rừng cục bộ theo băng/rạch: Áp dụng tại các ô tiêu chuẩn có mật độ thấp ($640\text{ cây/ha}$ như OTC 5), đưa các loài bản địa ưu thế họ Sao Dầu (Dipterocarpus alatus, Hopea odorata) vào làm giàu rừng.
LỘ TRÌNH QUẢN LÝ VÀ PHỤC HỒI LÂM SINH (2024 - 2028):
Năm 1 (2024): Thiết lập 5 ô định vị lâu dài (PSP) + Lập đường băng cản lửa
Năm 2 (2025): Tỉa thưa điều chỉnh mật độ tầng tái sinh (giữ lại 3.000 cây/ha cấp A)
Năm 3-4 (2026-2027): Làm giàu rừng bổ sung 200 cây/ha loài Giáng hương, Gõ đỏ
Năm 5 (2028): Đánh giá động thái trữ lượng (mục tiêu nâng lên > 50 m3/ha - Rừng nghèo)
Phân tích hiệu quả kinh tế - sinh thái (ROI sinh thái)
- Tiết kiệm chi phí đầu tư: So với phương pháp trồng rừng mới ($30 - 45\text{ triệu VNĐ/ha}$), giải pháp khoanh nuôi xúc tiến tái sinh dựa trên kết quả nghiên cứu chỉ tiêu tốn $3 - 5\text{ triệu VNĐ/ha/năm}$, tiết kiệm tới 85% ngân sách lâm sinh.
- Tiềm năng thị trường tín chỉ Carbon: Với tốc độ sinh trưởng bình quân dự kiến của họ Sao Dầu ($1,2 - 1,8\text{ cm đường kính/năm}$), trữ lượng rừng sẽ tăng từ $31,4\text{ m}^3\text{/ha}$ lên $> 70\text{ m}^3\text{/ha}$ sau 10 năm, tương đương khả năng hấp thụ bổ sung khoảng $25 - 35\text{ tấn }\text{CO}_2\text{/ha/năm}$, tạo nguồn thu bền vững từ thị trường tín chỉ carbon quốc tế.
Hạn chế và hướng phát triển
Hạn chế kỹ thuật
- Quy mô diện tích mẫu: Nghiên cứu giới hạn trên 5 OTC ($0,5\text{ ha}$) tại Tiểu khu 21, chưa bao phủ toàn bộ các kiểu lập địa trảng, bàu đất ngập nước theo mùa của VQG Lò Gò - Xa Mát.
- Tính chất thời điểm: Dữ liệu điều tra phản ánh trạng thái tĩnh tại thời điểm tháng 04 - 08/2023, chưa phản ánh đầy đủ động thái tăng trưởng chu kỳ nhiều năm (Permanent Sample Plots).
Định hướng nghiên cứu mở rộng
- Tích hợp viễn thám đa phổ vệ tinh (Sentinel-2, PlanetScope) và thiết bị bay không người lái (UAV LiDAR) để xây dựng bản đồ số hóa trữ lượng và sinh khối trên toàn bộ diện tích VQG.
- Xây dựng ô định vị quan trắc động thái sinh thái rừng vĩnh cửu (1 ha) để theo dõi tốc độ tăng trưởng đường kính ($D_{1.3}$), chiều cao ($H_{vn}$) và chu chuyển vật chất dinh dưỡng trong đất.
Đối tượng hưởng lợi
+-------------------------------------------------------------------------------+
| MA TRẬN GIÁ TRỊ CÁC BÊN LIÊN QUAN |
+-------------------------------------------------------------------------------+
| SINH VIÊN & HỌC VIÊN CAO HỌC | KỸ SƯ LÂM NGHIỆP & LẬP TRÌNH VIÊN SINH THÁI|
| - Cung cấp dữ liệu thực tế | - Thuật toán chuẩn tính IVI%, Shannon |
| - Mẫu biểu điều tra chuẩn hóa | - Script tự động hóa xử lý sinh thái R/Py |
|--------------------------------+----------------------------------------------|
| BAN QUẢN LÝ VƯỜN QUỐC GIA | CỘNG ĐỒNG KHOA HỌC & BẢO TỒN |
| - Phác đồ can thiệp lâm sinh | - CSDL đa dạng sinh học Nam Đông Dương |
| - Tối ưu 85% chi phí bảo vệ | - Cập nhật hiện trạng 9 loài IUCN |
+-------------------------------------------------------------------------------+
- Sinh viên ngành Quản lý Tài nguyên Rừng & Sinh thái: Nắm vững phương pháp luận lập ô tiêu chuẩn, kỹ năng đo đếm lâm học và quy trình xử lý số liệu trên Statgraphics, PRIMER 6.
- Ban Quản lý VQG Lò Gò - Xa Mát: Có cơ sở khoa học trực tiếp để bảo vệ và phục hồi Tiểu khu 21, làm căn cứ báo cáo định kỳ diễn biến tài nguyên rừng lên Cục Kiểm lâm và Bộ NN&PTNT.
- Các tổ chức phi chính phủ bảo tồn (IUCN, WWF, BirdLife International): Dữ liệu phân bố của 9 loài cây gỗ nguy cấp phục vụ chương trình bảo tồn loài nguy cấp tại vùng sinh cảnh chuyển tiếp Nam Đông Dương.
Câu hỏi thường gặp
1. Điều kiện kỹ thuật và trang thiết bị ngoại nghiệp để nhân rộng nghiên cứu này là gì?
Cần trang bị điện thoại thông minh tích hợp phần mềm định vị GPS (như Locus Map Free, Avenza Maps), thước dây chuyên dụng đo chu vi $C_{1.3}$ (quy đổi đường kính $D_{1.3}$), thước đo cao bách phân Haglöf (hoặc sào đo cao chuyên dụng cho tầng tái sinh), la bàn địa chất, phần mềm MS Excel (Data Analysis ToolPak) và R/PRIMER 6 để xử lý dữ liệu.
2. Vì sao rừng có trữ lượng chỉ 31,4 $\text{m}^3\text{/ha}$ nhưng chỉ số đa dạng Shannon-Wiener vẫn đạt mức cao ($H' = 3,24$)?
Trữ lượng thấp phản ánh kích thước thân cây nhỏ ($D_{1.3\text{ tb}} = 10,7\text{ cm}$) do hậu quả của khai thác trong quá khứ (trạng thái nghèo kiệt). Tuy nhiên, số lượng loài xuất hiện lại rất phong phú (37 loài thuộc 23 họ) và mức độ phân bố cá thể giữa các loài đồng đều ($J' = 0,94$, $\lambda' = 0,03$), tạo nên giá trị đa dạng loài cao.
3. Làm thế nào để tích hợp dữ liệu OTC vào hệ thống thông tin địa lý (GIS) của VQG?
Sử dụng tọa độ UTM (Zone 48N, WGS84) của 5 OTC đưa vào phần mềm QGIS hoặc ArcGIS Pro. Dữ liệu thuộc tính ($N$, $G$, $M$, $IVI%$, loài ưu thế) được liên kết qua bảng dữ liệu vector (Attribute Table) để xây dựng lớp bản đồ chuyên đề trạng thái rừng và phân vùng bảo tồn nghiêm ngặt.
4. Tần suất theo dõi và bảo trì số liệu định kỳ là bao lâu?
Đối với tầng cây gỗ lớn ($D_{1.3} \ge 8\text{ cm}$), chu kỳ đo đếm lặp lại định kỳ là 5 năm/lần. Đối với tầng cây tái sinh non ($H < 4\text{ m}$), tần suất theo dõi nên duy trì 1 năm/lần (vào cuối mùa mưa) để đánh giá tỷ lệ sống sót qua mùa khô hạn kéo dài tại Tây Ninh.
5. Chi phí khoanh nuôi phục hồi rừng nghèo kiệt so với trồng rừng mới chênh lệch thế nào?
Trồng mới rừng gỗ lớn bản địa đòi hỏi chi phí từ $30 - 45\text{ triệu VNĐ/ha}$ (cây giống, làm đất, chăm sóc 3 năm đầu). Trong khi đó, khoanh nuôi xúc tiến tái sinh tự nhiên dựa trên mật độ cây non có sẵn ($12.027\text{ cây/ha}$) chỉ tốn khoảng $3 - 5\text{ triệu VNĐ/ha/năm}$ cho công tác tuần tra bảo vệ và phát luổng cục bộ, mang lại hiệu quả chi phí - lợi ích tối ưu.
Kết luận
Tóm tắt các thành tựu cốt lõi
Khóa luận đã hoàn thành xuất sắc các mục tiêu nghiên cứu với hệ thống số liệu định lượng chi tiết:
- Xác định hệ thực vật gồm 37 loài thân gỗ thuộc 23 họ, hình thành nên ưu hợp sinh thái đặc trưng: Dầu song nàng (Dipterocarpus dyeri) + Dầu rái (Dipterocarpus alatus) + Sầm lá lớn (Memecylon harmandii).
- Định lượng chính xác các chỉ tiêu lâm học: Mật độ $N = 738\text{ cây/ha}$, Tiết diện ngang $G = 6,9\text{ m}^2\text{/ha}$, Trữ lượng $M = 31,4\text{ m}^3\text{/ha}$, chuẩn hóa đúng trạng thái rừng nghèo kiệt theo Thông tư 33/2018/TT-BNNPTNT.
- Phát hiện tiềm năng phục hồi tự nhiên vượt bậc với mật độ tái sinh đạt $12.027\text{ cây/ha}$, trong đó 87,8% có nguồn gốc từ hạt và 43,9% đạt phẩm chất tốt (Cấp A).
- Xác lập các chỉ số đa dạng sinh thái cao ($d = 6,98$; $J' = 0,94$; $H' = 3,24$; $\lambda' = 0,03$) cùng danh lục 9 loài cây gỗ nguy cấp cần ưu tiên bảo tồn theo chuẩn IUCN (2022).
Đóng góp kỹ thuật và giá trị thực tiễn
Công trình cung cấp cơ sở dữ liệu nền tảng, xóa bỏ các nhận định chủ quan về rừng nghèo kiệt, chứng minh rằng lâm phần Tiểu khu 21 tại VQG Lò Gò - Xa Mát sở hữu nội lực phục hồi mạnh mẽ. Đây là luận cứ khoa học then chốt giúp Ban Quản lý VQG Lò Gò - Xa Mát và Chi cục Kiểm lâm Tây Ninh triển khai chính sách khoanh nuôi, bảo vệ nghiêm ngặt, tiết kiệm tối đa ngân sách nhà nước so với các chương trình trồng rừng thay thế.
Kêu gọi hành động (Call to Action)
Các nhà quản lý lâm nghiệp, kỹ sư sinh thái và nghiên cứu viên quan tâm đến mô hình tính toán định lượng hoặc có nhu cầu tiếp cận bộ cơ sở dữ liệu mẫu biểu ngoại nghiệp chi tiết, vui lòng kết nối để trao đổi học thuật, chuyển giao giải pháp phần mềm phân tích sinh thái và phối hợp triển khai các hệ thống ô định vị lâu dài (PSP) bảo tồn rừng đặc dụng Việt Nam.
