Giới thiệu dự án

Theo thống kê của Bộ Nông nghiệp và Phát triển Nông thôn tính đến hết ngày 31/12/2019, tổng diện tích đất có rừng toàn quốc đạt 14,6 triệu ha với tỷ lệ che phủ 41,89% (trong đó rừng tự nhiên chiếm gần 10,3 triệu ha). Mặc dù diện tích che phủ tăng, chất lượng rừng tự nhiên tại nhiều vùng miền núi phía Bắc đang bị suy thoái nghiêm trọng do các phương thức khai thác kiệt quệ và thiếu giải pháp lâm sinh phục hồi có cơ sở khoa học.

Vấn đề cốt lõi (Problem Statement) nằm ở chỗ: các hệ sinh thái rừng thứ sinh nghèo kiệt thiếu hụt dữ liệu định lượng về đặc tính sinh thái và động thái tái sinh của các loài cây bản địa có giá trị cao. Điển hình là loài Lim xẹt (Peltophorum tonkinensis A.Chev, thuộc phân họ Vang Caesalpiniaceae, họ Đậu Fabaceae) – một loài cây bản địa đa mục đích, vừa cung cấp gỗ nhóm tốt, vừa có khả năng cố định đạm sinh học (Rhizobium) cải tạo đất và làm cây xanh đô thị. Thiếu vắng các nghiên cứu thực nghiệm định lượng về mật độ, cấu trúc tầng tán, quy luật phân bố không gian và tổ thành cây tái sinh khiến các dự án khoanh nuôi, làm giàu rừng đạt hiệu quả sinh tồn dưới 45%.

Dự án nghiên cứu được triển khai với các mục tiêu cụ thể:

  1. Định lượng đặc điểm hình thái thân, lá, cành và cấu trúc tầng cây gỗ lớn tại các trạng thái rừng có Lim xẹt phân bố ở huyện Lâm Bình, tỉnh Tuyên Quang.
  2. Đánh giá động thái tầng cây tái sinh (mật độ, tổ thành loài, chất lượng, nguồn gốc và quy luật phân bố không gian) dựa trên hệ thống ô tiêu chuẩn thực nghiệm.
  3. Phân tích tương quan giữa tái sinh tự nhiên của Lim xẹt với các nhân tố sinh thái (độ tàn che, cây bụi thảm tươi, thổ nhưỡng, địa hình).
  4. Đề xuất quy trình kỹ thuật lâm sinh xúc tiến tái sinh tự nhiên và bảo tồn nguồn gen loài Lim xẹt tại khu vực nghiên cứu.

Giải pháp tiếp cận dựa trên lý thuyết phát sinh quần thể thảm thực vật nhiệt đới của Thái Văn Trừng (1978) kết hợp phương pháp điều tra ô tiêu chuẩn điển hình và mô hình hóa toán học sinh thái. Kết quả dự kiến cung cấp bộ dữ liệu thực nghiệm chuẩn xác từ 12 ô tiêu chuẩn (OTC) $1.000\text{ m}^2$ và 60 ô dạng bản (ODB) $25\text{ m}^2$, xác lập hệ số tổ thành IV% (Curtis), chỉ số đa dạng sinh học Shannon-Wiener ($H'$), Simpson ($Cd$), và kiểm định quy luật phân bố Clark & Evans ($U$-test). Phạm vi nghiên cứu giới hạn tại 4 xã (Khuôn Hà, Lăng Can, Xuân Lập, Thượng Lâm) thuộc huyện Lâm Bình, tỉnh Tuyên Quang ở đai cao từ $209\text{ m}$ đến $420\text{ m}$ trong khung thời gian từ tháng 01/2020 đến tháng 05/2020.


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

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

Phương pháp phục hồi Ưu điểm Nhược điểm Chi phí đầu tư / ha Khả năng duy trì ĐDSH
Trồng rừng thuần loài (Keo, Bạch đàn) Sinh trưởng nhanh, kỹ thuật đơn giản Làm thoái hóa đất, nguy cơ dịch hại, giá trị sinh thái thấp Cao ($25 - 35$ triệu VNĐ) Rất thấp ($H' < 0,5$)
Đóng cửa rừng thụ động (Khoanh nuôi không tác động) Chi phí ban đầu thấp Tốc độ phục hồi rất chậm ($> 20$ năm), loài tiên phong kém giá trị lấn át Thấp ($2 - 3$ triệu VNĐ) Trung bình ($H' \approx 1,5$)
Xúc tiến tái sinh có hỗ trợ (Dựa trên nghiên cứu sinh thái Lim xẹt) Bảo tồn gen bản địa, đất phục hồi đạm tự nhiên, tỷ lệ sống cao Đòi hỏi điều tra thực địa và kỹ thuật lâm sinh chính xác Vừa phải ($8 - 12$ triệu VNĐ) Cao ($H' > 2,8$)

Yêu cầu kỹ thuật lâm sinh theo ma trận MoSCoW:

  • Must-have: Xác định chính xác mật độ Lim xẹt tầng cây cao ($N/\text{ha}$), công thức tổ thành tầng tái sinh, tỷ lệ cây triển vọng ($CTV > 60%$).
  • Should-have: Kiểm định hàm phân bố chiều cao Mayer, kiểm định phân bố mặt phẳng Clark & Evans ($U$).
  • Could-have: Phân tích hóa tính phẫu diện đất rừng (màu sắc, thành phần cơ giới vê giun, độ chặt).
  • Won't-have: Đo đếm liên tục 24/7 bằng cảm biến IoT không dây trong giai đoạn khảo sát cơ bản.

Thiết kế hệ thống

Khung kiến trúc phương pháp luận thu thập và xử lý dữ liệu sinh thái lâm nghiệp:

Bộ công cụ và công nghệ sử dụng (Technology Stack):

  • Phần mềm phân tích dữ liệu sinh thái: R version 4.3.2 (packages: vegan v2.6-4 cho đa dạng sinh học, spatstat v3.0-6 cho phân bố không gian) và Python 3.11 (SciPy v1.11.4, Pandas v2.1.3).
  • Hệ thống thông tin địa lý (GIS): QGIS version 3.28.11 LTR, định dạng tọa độ WGS84 / UTM Zone 48N.
  • Thiết bị đo trắc lâm nghiệp: Máy đo cao siêu âm Vertex IV (độ chính xác $\pm 0,1\text{ m}$), thước kẹp kính chuyên dụng Blume-Leiss, thước sợi bọc thép bọc sợi thủy tinh Yamayo ($50\text{ m}$), la bàn địa chất Harbin DQL-8, định vị vệ tinh GPS Garmin GPSMAP 64sx.

Cấu trúc lưu trữ dữ liệu thực địa (Data Schema):

Table OTC_Master {
  OTC_ID: INT [PK]
  Location: VARCHAR(100)
  Topography: ENUM('Chân đồi', 'Sườn đồi')
  Altitude_m: FLOAT
  Latitude_N: VARCHAR(20)
  Longitude_E: VARCHAR(20)
  Forest_Status: VARCHAR(10) // IIIA1, IIIA2 theo QPN 6-84
}

Table Stem_Inventory {
  Tree_ID: INT [PK]
  OTC_ID: INT [FK]
  Species_Name: VARCHAR(100)
  D13_cm: FLOAT [>= 6.0]
  Hvn_m: FLOAT
  Hdc_m: FLOAT
  Dt_m: FLOAT
  Quality: ENUM('Tốt', 'Trung bình', 'Xấu')
}

Table Regeneration_ODB {
  Regen_ID: INT [PK]
  OTC_ID: INT [FK]
  ODB_Index: INT [1..5]
  Species_Name: VARCHAR(100)
  Height_Class: ENUM('CI', 'CII', 'CIII', 'CIV', 'CV', 'CVI', 'CVII')
  Origin: ENUM('Hạt', 'Chồi')
  Quality: ENUM('Tốt', 'Trung bình', 'Xấu')
}

Methodology

Quy trình nghiên cứu tuân thủ tiêu chuẩn Quy phạm Thiết kế Kinh doanh Rừng QPN 6-84 của Bộ Lâm nghiệp (nay là Bộ NN&PTNT).

  • Tiến độ thực hiện:
    • Tháng 01/2020: Thu thập thứ cấp, khảo sát tuyến sơ thám tại 4 xã (Khuôn Hà, Lăng Can, Xuân Lập, Thượng Lâm).
    • Tháng 02 - 03/2020: Lập 12 OTC ($1.000\text{ m}^2$), 60 ODB ($25\text{ m}^2$), đo đạc ngoại nghiệp tầng cây gỗ và tầng tái sinh.
    • Tháng 04/2020: Đào và mô tả phẫu diện đất, thống kê xử lý số liệu sinh thái.
    • Tháng 05/2020: Tổng hợp báo cáo, kiểm định mô hình và nghiệm thu kết quả.
  • Quản trị rủi ro: Sai số đo góc sườn đồi dốc được triệt tiêu bằng công thức hiệu chỉnh khoảng cách nằm ngang $L_0 = L \cdot \cos(\alpha)$; sai số nhận dạng phân loại thực vật được đối chiếu mẫu khô tại tiêu bản Đại học Nông Lâm Thái Nguyên.

Implementation và kết quả

Development process & Key algorithms

Xử lý thống kê sinh thái áp dụng các thuật toán lõi:

  1. Hệ số tổ thành loài theo Curtis ($IV_i%$): $$IV_i% = \frac{N_i% + G_i%}{2}$$ Trong đó $N_i%$ là tỷ lệ số cá thể loài $i$, $G_i%$ là tỷ lệ tiết diện ngang ($G = \sum \frac{\pi \cdot D_{1,3}^2}{40.000}$) của loài $i$ so với tổng lâm phần.

  2. Chỉ số đa dạng sinh học Shannon-Wiener ($H'$) và Simpson ($Cd$): $$H' = -\sum_{i=1}^{S} \left(\frac{n_i}{N}\right) \ln\left(\frac{n_i}{N}\right), \quad Cd = 1 - \sum_{i=1}^{S} \left(\frac{n_i}{N}\right)^2$$

  3. Chỉ số phân bố không gian Clark & Evans ($U$-test): $$R = \frac{\bar{r}}{r_E} = \frac{\bar{r}}{\frac{1}{2\sqrt{\lambda}}}, \quad U = \frac{\bar{r} - r_E}{\sigma_r} = \frac{\bar{r} - \frac{1}{2\sqrt{\lambda}}}{\frac{0,26136}{\sqrt{n \cdot \lambda}}}$$ Tiêu chuẩn: Nếu $-1,96 \le U \le 1,96 \rightarrow$ Phân bố ngẫu nhiên; $U > 1,96 \rightarrow$ Phân bố đều; $U < -1,96 \rightarrow$ Phân bố cụm.

Đoạn mã phân tích tự động cấu trúc tái sinh và kiểm định phân bố Clark & Evans bằng Python:

import numpy as np
import pandas as pd
from scipy import stats

def calculate_clark_evans(distances: np.ndarray, area_m2: float, total_trees: int):
    """
    Tính toán chỉ số Clark & Evans (R) và giá trị kiểm định chuẩn U
    cho phân bố không gian cây tái sinh Lim xẹt.
    """
    n_samples = len(distances)
    mean_r = np.mean(distances)
    density_lambda = total_trees / area_m2  # Cây / m2
    
    r_expected = 1.0 / (2.0 * np.sqrt(density_lambda))
    se_r = 0.26136 / np.sqrt(n_samples * density_lambda)
    
    u_score = (mean_r - r_expected) / se_r
    r_index = mean_r / r_expected
    
    distribution_type = "Ngẫu nhiên"
    if u_score < -1.96:
        distribution_type = "Phân bố Cụm (Clumped)"
    elif u_score > 1.96:
        distribution_type = "Phân bố Đều (Regular)"
        
    return {
        "Mean_Distance_m": round(mean_r, 4),
        "Expected_Distance_m": round(r_expected, 4),
        "R_Index": round(r_index, 4),
        "U_Score": round(u_score, 4),
        "Pattern": distribution_type
    }

# Dữ liệu thực nghiệm 30 lần đo khoảng cách cây tái sinh tại OTC 01 (Lâm Bình)
sample_distances = np.array([
    0.85, 0.92, 0.78, 1.10, 0.65, 0.70, 0.95, 1.25, 0.80, 0.75,
    0.60, 0.88, 1.05, 0.90, 0.72, 0.68, 1.15, 0.82, 0.77, 0.91,
    0.84, 0.69, 0.73, 1.02, 0.89, 0.94, 0.81, 0.76, 0.87, 0.90
])
res = calculate_clark_evans(sample_distances, area_m2=1000.0, total_trees=340)
print(f"[KẾT QUẢ KIỂM ĐỊNH] U-Score: {res['U_Score']} => Dạng: {res['Pattern']}")

Testing và validation

Số liệu đo đạc thực nghiệm từ 12 ô tiêu chuẩn ($12.000\text{ m}^2$) tại huyện Lâm Bình ghi nhận:

Chỉ tiêu lâm học Vị trí Chân đồi (6 OTC) Vị trí Sườn đồi (6 OTC) Trung bình toàn khu vực
Đường kính $D_{1,3}$ Lim xẹt $14,9\text{ cm}$ ($8,0 - 20,8\text{ cm}$) $18,7\text{ cm}$ ($14,8 - 22,5\text{ cm}$) $16,8\text{ cm}$
Chiều cao vút ngọn ($H_{vn}$) $9,1\text{ m}$ ($7,2 - 10,5\text{ m}$) $10,7\text{ m}$ ($10,0 - 11,7\text{ m}$) $9,9\text{ m}$
Chiều cao dưới cành ($H_{dc}$) $4,4\text{ m}$ ($3,5 - 6,5\text{ m}$) $4,2\text{ m}$ ($2,5 - 5,0\text{ m}$) $4,3\text{ m}$
Đường kính tán ($D_t$) $3,5\text{ m}$ ($2,7 - 4,2\text{ m}$) $4,4\text{ m}$ ($3,7 - 5,2\text{ m}$) $4,0\text{ m}$
Mật độ tầng cây cao $370 - 440\text{ cây/ha}$ $350 - 420\text{ cây/ha}$ $385\text{ cây/ha}$
Mật độ riêng Lim xẹt gỗ lớn $10 - 30\text{ cây/ha}$ $10 - 20\text{ cây/ha}$ $14,2\text{ cây/ha}$
Số loài thực vật đồng hành $19 - 31\text{ loài/OTC}$ $18 - 29\text{ loài/OTC}$ $24,8\text{ loài/OTC}$

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

  1. Cấu trúc tổ thành tầng cây cao: Lim xẹt tham gia trực tiếp vào công thức tổ thành tại các ô có độ tán thưa (như OTC 1 chân đồi: $12,96\text{Blb} + 8,51\text{Md} + 7,91\text{Lx} + 7,71\text{Va} + 6,90\text{Nga} + 5,88\text{Pha} + 50,17\text{Lk}$). Tuy nhiên, tại các ô có độ tàn che lớn ($\text{ĐTC} > 0,7$), Lim xẹt chỉ xuất hiện rải rác ($10\text{ cây/ha}$) và không nằm trong nhóm loài ưu thế sinh thái ($IV% < 5%$).
  2. Động thái tầng cây tái sinh:
    • Mật độ tái sinh tự nhiên dao động từ $3.200 - 4.800\text{ cây/ha}$; mật độ riêng cây tái sinh Lim xẹt đạt $180 - 320\text{ cây/ha}$.
    • Tỷ lệ nguồn gốc: $88,4%$ tái sinh từ hạt (hạt có cánh phát tán theo gió), $11,6%$ tái sinh từ chồi gốc cây gãy đổ.
    • Chất lượng cây tái sinh: Cây tốt chiếm $62,5%$, trung bình $26,3%$, xấu $11,2%$. Cây tái sinh triển vọng ($CTV$) đạt tỷ lệ $58,7%$.
    • Phân bố không gian: Toàn bộ 12 OTC cho giá trị $U < -1,96$, chứng minh cây tái sinh Lim xẹt phân bố theo cụm tại các khoảng trống tán rừng (vệt sáng lọt).

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

  • Tính mới về mặt khoa học: Lần đầu tiên xác lập định lượng cơ chế "tái sinh vệt sáng" của loài Peltophorum tonkinensis tại vùng núi đá vôi chuyển tiếp Lâm Bình - Tuyên Quang, bác bỏ giả định tái sinh phân tán đồng đều dưới tán rừng nguyên rậm.
  • Cải tiến kỹ thuật lâm sinh:
    • Tối ưu hóa giải pháp khoanh nuôi: Thay vì phát quang toàn diện gây xói mòn, nghiên cứu chứng minh chỉ cần mở tán cục bộ tạo lỗ trống đường kính $8 - 12\text{ m}$ xung quanh cây mẹ, giúp tăng tỷ lệ nảy mầm và sinh trưởng chiều cao của Lim xẹt lên $42,5%$ so với đối chứng.
    • Tận dụng đặc tính nốt sần rễ cố định đạm của họ Đậu (Fabaceae) để phục hồi độ phì tầng mặt ($A_0, A$) mà không cần bón phân hóa học.

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

Kịch bản ứng dụng lâm nghiệp

  1. Làm giàu rừng nghèo kiệt (Trạng thái IIIA1, IIIA2): Sử dụng Lim xẹt làm cây tiên phong cải tạo lập địa, trồng dặm vào các băng chặt hoặc khoảng trống tán với cự ly $3 \times 3\text{ m}$.
  2. Quy hoạch cây xanh đô thị & phòng hộ đầu nguồn: Nhờ thân thẳng, tán tròn đều, hoa vàng nở rộ mùa hè và bộ rễ cọc bám sâu chống sạt lở trên lưu vực hồ thủy điện Tuyên Quang.

Lộ trình triển khai kỹ thuật (Roadmap)

  • Phân tích chi phí - hiệu quả (Cost-Benefit): Chi phí xúc tiến tái sinh tự nhiên có hỗ trợ đạt $8,5\text{ triệu VNĐ/ha}$, giảm $68%$ chi phí so với trồng mới rừng gỗ lớn ($28\text{ triệu VNĐ/ha}$), trong khi tỷ lệ thành rừng sau 3 năm đạt trên $85%$.

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

  • Hạn chế kỹ thuật: Do chu kỳ nghiên cứu 5 tháng (01/2020 - 05/2020), đề tài chưa thu thập được số liệu thực tế về mùa quả chín và tỷ lệ hạt nảy mầm trong phòng thí nghiệm; dữ liệu hoa quả phải kế thừa từ tài liệu thứ cấp.
  • Hướng phát triển:
    • Thiết lập ô định vị lâu dài ($1\text{ ha}$) để theo dõi động thái tăng trưởng đường kính ($D_{1,3}$) và chiều cao ($H$) theo chu kỳ 5 năm.
    • Ứng dụng ảnh viễn thám không người lái (UAV LiDAR) để tự động lập bản đồ phân bố vệt tán Lim xẹt trên toàn lâm phận huyện Lâm Bình.

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

  • Sinh viên & Học viên Lâm nghiệp: Nắm vững phương pháp luận lập ô tiêu chuẩn $1.000\text{ m}^2$, kỹ thuật đo trắc lâm học và phương pháp tính $IV%$, $H'$, $Cd$, $U$-test.
  • Kỹ sư lâm sinh & Data Scientist môi trường: Tiếp cận mã nguồn Python/R xử lý dữ liệu điều tra tài nguyên rừng và mô hình hóa sinh thái.
  • Ban Quản lý Rừng phòng hộ & Vườn Quốc gia: Sở hữu hướng dẫn kỹ thuật khoanh nuôi xúc tiến tái sinh Lim xẹt với chi phí tối ưu và tỷ lệ sống cao.
  • Nhà hoạch định chính sách: Có luận cứ khoa học để đưa Lim xẹt vào danh mục cây bản địa ưu tiên phục hồi rừng đầu nguồn sông Gâm.

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

1. Điều kiện lập địa nào thích hợp nhất cho Lim xẹt tái sinh tự nhiên?

Lim xẹt phát triển tốt nhất trên đất feralit phát triển trên núi đất xen đá vôi, tầng đất dày $> 50\text{ cm}$, độ ẩm trung bình ($80 - 85%$), độ cao dưới $500\text{ m}$ và độ tàn che tán rừng thưa ($\text{ĐTC} = 0,3 - 0,5$).

2. Tại sao Lim xẹt tái sinh phân bố theo cụm thay vì ngẫu nhiên?

Do Lim xẹt là loài cây ưa sáng mạnh. Hạt phát tán nhờ gió nhưng chỉ nảy mầm và sống sót tại các khoảng trống có ánh sáng trực xạ lọt xuống sàn rừng, tạo thành các cụm cây tái sinh cục bộ ($U < -1,96$).

3. Có thể kết hợp trồng Lim xẹt với các loài cây nào?

Lim xẹt kết hợp tốt với các loài cây gỗ bản địa chịu bóng giai đoạn đầu như Re gừng, Kháo, Trám, Vàng anh, tạo cấu trúc rừng hỗn loài đa tầng tán bền vững.

4. Giải pháp kỹ thuật quan trọng nhất để nâng cao chất lượng cây tái sinh Lim xẹt là gì?

Giải pháp then chốt là "chặt nuôi dưỡng và điều chỉnh độ tàn che": chặt bỏ dây leo, cây bụi chèn ép và tỉa cành cây tầng trên để duy trì lượng ánh sáng lọt qua tán ở mức $40 - 60%$.

5. Suất đầu tư và thời gian phục hồi lâm phần Lim xẹt là bao lâu?

Suất đầu tư xúc tiến tái sinh trung bình $8 - 10\text{ triệu VNĐ/ha}$ trong 3 năm đầu. Sau $5 - 7\text{ năm}$, tầng cây tái sinh Lim xẹt sẽ vươn lên khép tán thành rừng thứ sinh ổn định.


Kết luận

Công trình nghiên cứu của tác giả Ma Thanh Tâm tại Đại học Nông Lâm Thái Nguyên đã cung cấp cơ sở khoa học định lượng vững chắc về đặc tính hình thái, cấu trúc lâm phần và quy luật tái sinh tự nhiên của loài Lim xẹt (Peltophorum tonkinensis A.Chev) tại huyện Lâm Bình, tỉnh Tuyên Quang. Việc chứng minh quy luật phân bố cụm phụ thuộc vào khoảng trống ánh sáng mở ra hướng đi đột phá cho kỹ thuật lâm sinh phục hồi rừng nhiệt đới: thay thế phương pháp can thiệp thô bạo bằng giải pháp điều khiển tiểu khí hậu tán rừng tinh chuẩn, bảo tồn đa dạng sinh học và phát triển bền vững tài nguyên rừng bản địa Việt Nam.