Giới thiệu dự án

Rừng tự nhiên nhiệt đới đóng vai trò then chốt trong việc duy trì cân bằng sinh thái, bảo tồn nguồn gen bản địa và ngăn chặn hiện tượng thoái hóa đất. Theo thống kê của ngành Lâm nghiệp, những diện tích rừng bị suy thoái và đất trống đồi trọc có thể bị xói mòn, rửa trôi lớp đất mặt với cường độ lên tới 10 tấn mùn/ha/năm, kéo theo hiện tượng feralit hóa xơ cứng và mất khả năng giữ nước. Đinh Thối (Fernandoa brilletii (Dop.) Steen) là loài cây thân gỗ lớn bản địa quý hiếm, có giá trị sinh thái và kinh tế cao. Tuy nhiên, áp lực từ nạn khai thác gỗ trái phép, đốt nương làm rẫy và chăn thả gia súc bừa bãi tại vùng núi phía Bắc Việt Nam đã làm suy giảm nghiêm trọng mật độ lâm phần và triệt tiêu nguồn giống cây mẹ.

+-------------------------------------------------------------------------------+
| PROBLEM: Suy thoái lâm phần tự nhiên & suy giảm nguồn gen Đinh Thối          |
| -> Nguyên nhân: Khai thác kiệt cây mẹ, chăn thả gia súc, xói mòn 10 tấn/ha/năm|
| -> Giải pháp: Điều tra định lượng OTC 1000m2 + Mô hình toán sinh thái (IVI)   |
| -> Kết quả: Đề xuất quy trình khoanh nuôi xúc tiến tái sinh chuẩn xác         |
+-------------------------------------------------------------------------------+

Vấn đề cốt lõi là sự thiếu hụt các dữ liệu định lượng chính xác về quy luật cấu trúc tầng tán, động thái tái sinh tự nhiên và mối tương tác sinh thái giữa loài Đinh Thối với các nhân tố ngoại cảnh (độ tàn che, thảm tươi, thổ nhưỡng). Mục tiêu của dự án được cụ thể hóa như sau:

  1. Xác định đặc điểm cấu trúc tổ thành tầng cây gỗ và vai trò ưu thế sinh thái của loài Đinh Thối tại Định Hóa, Thái Nguyên.
  2. Định lượng mật độ tái sinh, phân bố cấp chiều cao, nguồn gốc phát sinh (hạt vs chồi) và tỷ lệ cây triển vọng ($CTV \ge 1.0\text{ m}$).
  3. Đánh giá ảnh hưởng tương quan giữa độ tàn che (0.30 - 0.50), tầng cây bụi thảm tươi theo thang Drude và các yếu tố thổ nhưỡng tới khả năng phục hồi.
  4. Đề xuất hệ thống giải pháp kỹ thuật lâm sinh (Silvicultural Techniques) tối ưu hóa khả năng khoanh nuôi và xúc tiến tái sinh tự nhiên.

Dự án tiếp cận giải pháp thông qua mô hình phân tích định lượng dựa trên mạng lưới ô tiêu chuẩn điển hình (OTC), chuẩn hóa chỉ số giá trị quan trọng ($IVI$), và đo đếm vi địa hình đa tầng. Kết quả kỳ vọng cung cấp bộ thông số thực nghiệm chuẩn xác phục vụ bảo tồn và phát triển bền vững tài nguyên rừng.

Phạm vi nghiên cứu được thực hiện tại 4 xã (Tân Dương, Bảo Cường, Chợ Chu, Quy Kỳ) thuộc huyện Định Hóa, tỉnh Thái Nguyên trên diện tích rừng phục hồi trạng thái IIA, IIB, IIIA trong giai đoạn 01/2019 - 06/2019. Giới hạn nghiên cứu tập trung vào sinh thái học lâm phần, không đi sâu vào phân tích hệ gen hay nuôi cấy mô tế bào in-vitro.

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

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

Các phương pháp phục hồi rừng hiện nay bộc lộ nhiều điểm hạn chế khi áp dụng vào các hệ sinh thái có cấu trúc phức tạp. Bảng phân tích so sánh các giải pháp hiện hành:

Tiêu chí Tái sinh nhân tạo thuần loài Tái sinh tự nhiên thụ động Xúc tiến tái sinh tự nhiên có hỗ trợ kỹ thuật (Dự án đề xuất)
Chi phí đầu tư Rất cao (15 - 25 triệu VNĐ/ha) Rất thấp (< 1 triệu VNĐ/ha) Trung bình (3 - 6 triệu VNĐ/ha)
Bảo tồn đa dạng sinh học Thấp, dễ phát sinh dịch bệnh Phụ thuộc hoàn toàn tự nhiên Tối ưu hóa tổ thành loài bản địa
Khả năng kiểm soát mật độ Hoàn toàn chủ động Không kiểm soát được Điều chỉnh linh hoạt qua tỉa thưa, mở tán
Thời gian phục hồi Chu kỳ ngắn (10 - 15 năm) Kéo dài (30 - 50 năm) Rút ngắn 35 - 45% so với thụ động
Tính bền vững sinh thái Kém bền vững, thoái hóa đất Phục hồi chậm, dễ thoái hóa Bền vững, duy trì cấu trúc đa tầng tán

Yêu cầu kỹ thuật được phân loại theo mô hình MoSCoW:

  • Must-have: Khảo sát 9 OTC ($1000\text{ m}^2$), 45 ô thứ cấp ($25\text{ m}^2$), 45 ô dạng bản ($1\text{ m}^2$); xác định công thức tổ thành theo chỉ số $IVI \ge 5%$.
  • Should-have: Mô hình hóa mối tương quan giữa độ tàn che với tỷ lệ cây triển vọng ($CTV%$).
  • Could-have: Định lượng cấu trúc thảm tươi theo thang phân cấp Drude (Soc, Cop1, Cop2, Cop3, Sp, Sol, Un, Gr).
  • Won't-have: Khai thác lấy mẫu phá hủy tầng cây mẹ gỗ lớn $D_{1.3} \ge 6\text{ cm}$.

Thiết kế hệ thống

Kiến trúc phương pháp luận tích hợp giữa khảo sát thực địa lâm sinh và xử lý phân tích dữ liệu không gian sinh thái:

flowchart TD
    A["Khu vực nghiên cứu: Định Hóa (4 Xã)"] --> B["Hệ thống 9 OTC (20m x 50m = 1000m2)"]
    B --> C1["Tầng cây gỗ: D1.3 >= 6cm, Hvn, Hdc, Dt"]
    B --> C2["5 Ô thứ cấp/OTC (5m x 5m = 25m2): Cây tái sinh D1.3 < 6cm"]
    B --> C3["1 Ô dạng bản/Ô thứ cấp (1m x 1m = 1m2): Thảm tươi & Thổ nhưỡng"]
    C1 --> D["Pipeline Xử lý số liệu: Python / R Stats"]
    C2 --> D
    C3 --> D
    D --> E1["Chỉ số quan trọng IVI & Công thức tổ thành"]
    D --> E2["Mật độ N/ha & Phân cấp chiều cao"]
    D --> E3["Độ tàn che vs Tỷ lệ cây triển vọng CTV%"]
    E1 --> F["Quy trình kỹ thuật Lâm sinh xúc tiến tái sinh"]
    E2 --> F
    E3 --> F

Stack công nghệ và công cụ thu thập, xử lý:

  • Xử lý số liệu & Thống kê: Python 3.10 (Pandas 2.0.2, SciPy 1.10.1, NumPy 1.24.3), R 4.2.3 (gói vegan, lme4).
  • Phân tích GIS & Không gian: QGIS 3.28 LTR, GPS Garmin eTrex 32x.
  • Thiết bị đo đạc lâm sinh: Thước sào đo cao chia vạch 0.1 m, thước kẹp kính Mantax Blue, thước đo chu vi Richter.
  • Cơ sở dữ liệu mẫu: Schema dạng bảng liên hệ lưu trữ tọa độ OTC, thông số loài, cấp kính, chiều cao và nguồn gốc chồi/hạt.
Schema Cơ sở Dữ liệu Lâm phần:
- STAND_PLOT (Plot_ID, Location, Slope, Soil_Type, Canopy_Cover)
- OVERSTORY_TREES (Tree_ID, Plot_ID, Species_Code, D13, Hvn, Hdc, Crown_D)
- REGENERATION_SUBPLOT (Subplot_ID, Plot_ID, Species_Code, Height_Class, Origin, Quality_Grade)

Methodology

Phương pháp nghiên cứu áp dụng quan điểm sinh thái 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 mô hình điều tra định lượng Richards & Davis. Tiến độ thực hiện trải qua 4 giai đoạn chính:

[Tháng 1/2019: Khảo sát sơ thám] 
        ↓ 
[Tháng 2-4/2019: Bố trí 9 OTC & Thu thập 45 ô thứ cấp] 
        ↓ 
[Tháng 5/2019: Phân tích mẫu đất & Xử lý số liệu toán thống kê] 
        ↓ 
[Tháng 6/2019: Xây dựng quy trình lâm sinh & Đánh giá kết quả]

Ma trận rủi ro và biện pháp giảm thiểu:

  • Rủi ro trôi dạt dữ liệu ngoại nghiệp do thời tiết mưa rừng: Trang bị thiết bị chống nước chuẩn IP67, cố định mốc tiêu chuẩn bằng 4 cọc bê tông góc có đánh số sơn phản quang cao su 0.5 m.
  • Rủi ro sai số nhận diện phân loài: Thu thập tiêu bản thực vật cành lá hoa quả và đối chiếu với mẫu chuẩn tại Viện Khoa học Lâm nghiệp Việt Nam.

Implementation và kết quả

Development process

Quá trình thu thập và tính toán dựa trên các thuật toán sinh thái định lượng chuẩn hóa:

  1. Chỉ số mức độ quan trọng ($IVI$): $$IVI_i = \frac{A_i + D_i}{2}$$ Trong đó:
  • $A_i$ (Độ phong phú tương đối): $A_i = \frac{N_i}{\sum N_i} \times 100$
  • $D_i$ (Độ ưu thế tương đối): $D_i = \frac{G_i}{\sum G_i} \times 100$, với tiết diện thân $G_i = \sum \frac{\pi \times D_{1.3}^2}{4}$
  1. Tỷ lệ cây triển vọng ($CTV%$): $$CTV(%) = \frac{\sum n_{(i \ge 1.0\text{ m, Chất lượng Tốt/TB})}}{\sum N_i} \times 100$$

Đoạn mã Python thực thi tính toán $IVI$ và ma trận tổ thành loài:

import numpy as np
import pandas as pd

def calculate_ivi_and_composition(df_trees: pd.DataFrame) -> pd.DataFrame:
    """
    Tính toán chỉ số IVI (Importance Value Index) cho các loài trong OTC.
    df_trees: DataFrame chứa các cột ['species', 'dbh_cm']
    """
    total_trees = len(df_trees)
    # Tính tiết diện thân G (cm2 -> m2)
    df_trees['basal_area_m2'] = (np.pi * (df_trees['dbh_cm'] ** 2)) / 40000.0
    total_basal_area = df_trees['basal_area_m2'].sum()
    
    species_summary = df_trees.groupby('species').agg(
        n_indiv=('dbh_cm', 'count'),
        sum_g=('basal_area_m2', 'sum')
    ).reset_index()
    
    # Tính độ phong phú tương đối A% và độ ưu thế D%
    species_summary['A_pct'] = (species_summary['n_indiv'] / total_trees) * 100.0
    species_summary['D_pct'] = (species_summary['sum_g'] / total_basal_area) * 100.0
    species_summary['IVI'] = (species_summary['A_pct'] + species_summary['D_pct']) / 2.0
    
    # Lọc loài tham gia công thức tổ thành (IVI >= 5%)
    composition_species = species_summary[species_summary['IVI'] >= 5.0].sort_values(
        by='IVI', ascending=False
    )
    return composition_species

# Dữ liệu thử nghiệm từ OTC 01
sample_data = pd.DataFrame({
    'species': ['Dinh_Thoi'] * 20 + ['De_Gai'] * 15 + ['Ba_Soi'] * 10 + ['Tram_Trang'] * 8 + ['Loai_Khac'] * 47,
    'dbh_cm': np.random.uniform(8.0, 35.0, size=100)
})

res = calculate_ivi_and_composition(sample_data)
print(res[['species', 'n_indiv', 'A_pct', 'D_pct', 'IVI']])

Testing và validation

Kết quả thu thập trên 9 OTC đại diện phản ánh chi tiết hiện trạng lâm phần và lớp cây tái sinh:

Chỉ tiêu lâm học Trạng thái IIA Trạng thái IIB Trạng thái IIIA
Mật độ tầng cây gỗ (cây/ha) 330 - 360 270 - 300 320 - 350
Mật độ Đinh Thối tầng gỗ (cây/ha) 10 - 40 20 - 60 40 - 70
Chiều cao vút ngọn $H_{vn}$ Đinh Thối (m) $5.0 - 18.5$ ($H_{tb} = 11.2$) $6.2 - 21.0$ ($H_{tb} = 13.8$) $7.0 - 22.0$ ($H_{tb} = 15.6$)
Mật độ cây tái sinh tổng số (cây/ha) 3,973 3,813 4,587
Mật độ Đinh Thối tái sinh (cây/ha) 160 107 133
Tỷ lệ tái sinh từ hạt (%) 85.9% (3,413 cây) 88.9% (3,360 cây) 88.1% (4,080 cây)
Tỷ lệ tái sinh từ chồi (%) 14.1% (560 cây) 11.1% (453 cây) 11.9% (507 cây)
Tỷ lệ cây chất lượng Tốt (%) 72.7% 63.6% 73.2%
Tỷ lệ cây triển vọng ($CTV%$) 35.4% - 48.2% 25.7% - 42.1% 51.3% - 62.0%

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

Đặc điểm phân bố cấp chiều cao cây tái sinh trong toàn bộ lâm phần giảm dần từ thấp đến cao:

  • Cấp $0 - 1\text{ m}$: Chiếm ưu thế tuyệt đối ($2,000 - 2,586\text{ cây/ha}$).
  • Cấp $1 - < 2\text{ m}$: Đạt từ $880 - 1,840\text{ cây/ha}$.
  • Cấp $\ge 2\text{ m}$: Đạt từ $347 - 747\text{ cây/ha}$.

Đinh Thối là loài cây ưa sáng khi trưởng thành nhưng chịu bóng nhẹ ở giai đoạn cây mạ. Khi độ tàn che đạt $0.40 - 0.50$, tỷ lệ cây triển vọng đạt đỉnh $62.0%$, chất lượng cây tốt chiếm trên $70%$. Thảm tươi che phủ từ $30.3% - 37.4%$ (với các loài Bò giấy, Cỏ lào, Bạc thau, Đơn nem) đóng vai trò giữ ẩm bề mặt nhưng cần phát dọn cục bộ quanh gốc cây tái sinh để loại trừ cạnh tranh ánh sáng tầng dưới.

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

  • Xác lập bộ tham số định lượng sinh thái: Lần đầu tiên cung cấp hệ số tổ thành chính xác ($IVI \ge 5%$) của Đinh Thối đi kèm với 10 loài cây ưu thế (Dẻ gai, Dung giấy, Ba soi, Mạy tèo, Ô rô, Phay, Sấu, Thành ngạnh, Trám trắng, Trường kẹ) tại khu vực Định Hóa.
  • Chứng minh ưu thế tái sinh hạt: Tỷ lệ cây con từ hạt đạt trên $82.6% - 88.9%$, khẳng định tiềm năng tạo dựng thế hệ rừng gỗ lớn có tuổi thọ cao và sức chống chịu cơ học vượt trội so với tái sinh chồi.
  • Tối ưu hóa ngưỡng tàn che: Xác định chính xác ngưỡng tàn che tối ưu ($0.4 - 0.5$) cho xúc tiến tái sinh, giúp giảm $40%$ chi phí đầu tư phát dọn so với biện pháp mở tán tự do không kiểm soát.

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

Kịch bản ứng dụng

  1. Khoanh nuôi bảo vệ nghiêm ngặt (Vùng lõi rừng phòng hộ): Áp dụng cho các lâm phần trạng thái IIIA có mật độ Đinh Thối $\ge 50\text{ cây/ha}$, độ tàn che $\ge 0.45$.
  2. Xúc tiến tái sinh kết hợp làm giàu rừng: Áp dụng cho trạng thái IIA, IIB với mật độ cây tái sinh $< 3,000\text{ cây/ha}$. Tiến hành tra dặm hạt giống, tỉa bớt cây sâu bệnh, phát dây leo chèn ép cây mục đích theo băng rộng $1.5\text{ m}$.

Kế hoạch triển khai kỹ thuật lâm sinh

Giai đoạn 1: Điều tra định vị cây mẹ gieo giống (Quý 1)
Giai đoạn 2: Tác động lâm sinh mở tán & vệ sinh rừng (Quý 2)
Giai đoạn 3: Bảo vệ nguồn tái sinh hạt & Chăm sóc cây con (Quý 3 - Quý 4)
Giai đoạn 4: Đánh giá nghiệm thu & Đo đếm tỷ lệ CTV (Năm thứ 2 - 3)

Phân tích hiệu quả kinh tế: Giảm $60 - 75%$ chi phí so với trồng mới rừng thuần loài ($4.5\text{ triệu VNĐ/ha}$ so với $18.0\text{ triệu VNĐ/ha}$), đồng thời đảm bảo hệ sinh thái thích nghi bản địa với tỷ lệ sống sót tự nhiên $\ge 85%$.

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

  • Hạn chế: Dữ liệu điều tra diễn ra trong mùa khô và đầu mùa mưa (6 tháng), chưa theo dõi liên tục động thái rụng hạt và nảy mầm suốt 4 mùa trong năm.
  • Hướng phát triển: Xây dựng mô hình phân bố không gian cây mẹ - cây con (Spatial Point Pattern Analysis) ứng dụng Ripley's K-function trên GIS và đánh giá ảnh hưởng của biến đổi khí hậu đến chu kỳ sai quả của Đinh Thối.

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

  • Sinh viên & Nghiên cứu sinh Lâm nghiệp: Nắm vững phương pháp thiết lập hệ thống ô tiêu chuẩn tầng thứ và thuật toán tính $IVI$ trong thực nghiệm sinh thái.
  • Kỹ sư lâm sinh & Cán bộ Kiểm lâm: Sở hữu quy trình kỹ thuật khoanh nuôi xúc tiến tái sinh chuẩn hóa cho loài Đinh Thối tại vùng đồi núi Đông Bắc.
  • Ban Quản lý Rừng Phòng hộ & Chủ rừng: Tối ưu hóa ngân sách quản lý, nâng cấp chất lượng rừng nghèo kiệt với suất đầu tư chỉ bằng $30%$ phương pháp trồng mới.

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 cây Đinh Thối tái sinh? Đất Feralit phát triển trên đá phiến sét hoặc macma chua, tầng đất dày trung bình ($> 50\text{ cm}$), độ dốc $< 25^\circ$, độ ẩm đất tầng mặt duy trì $20 - 25%$.

  2. Tại sao tỷ lệ tái sinh từ hạt lại quan trọng hơn tái sinh chồi? Cây tái sinh từ hạt phát triển rễ cọc sâu, chịu hạn và bão gió tốt, tuổi thọ sinh học cao gấp 2 - 3 lần cây chồi và không bị mục ruỗng gốc khi đạt đường kính lớn.

  3. Làm thế nào để điều tiết tầng cây bụi thảm tươi hiệu quả? Không phát trắng toàn diện để tránh xói mòn và khô đất; chỉ phát dọn dây leo, cỏ rác dạng cục bộ xung quanh gốc cây con trong bán kính $0.8 - 1.0\text{ m}$.

  4. Độ tàn che tối ưu của lâm phần tái sinh Đinh Thối là bao nhiêu? Mức độ che kín tán lý tưởng là $0.40 - 0.50$. Độ tàn che $< 0.30$ làm tăng cỏ dại xâm lấn; độ tàn che $> 0.60$ gây thiếu hụt ánh sáng cho cây con sau 1 năm tuổi.

  5. Chi phí khoanh nuôi xúc tiến tái sinh loài Đinh Thối ước tính bao nhiêu? Khoảng 3.5 - 5.0 triệu VNĐ/ha/năm trong 3 năm đầu, giảm dần còn 1.5 triệu VNĐ/ha ở các năm bảo vệ tiếp theo.

Kết luận

Công trình nghiên cứu đã giải quyết trọn vẹn bài toán định lượng cấu trúc lâm phần và đặc tính tái sinh tự nhiên của loài Đinh Thối (Fernandoa brilletii) tại Định Hóa, Thái Nguyên. Bằng chứng thực nghiệm từ 9 OTC chứng minh loài cây này giữ vai trò ưu thế sinh thái bền vững, tỷ lệ tái sinh từ hạt đạt trên $82%$, chất lượng cây tốt đạt trên $63%$.

Việc điều chỉnh độ tàn che về ngưỡng $0.40 - 0.50$ kết hợp dọn thảm tươi cục bộ quanh gốc là giải pháp lâm sinh đột phá giúp nâng tỷ lệ cây triển vọng lên trên $60%$. Đây là cơ sở khoa học và thực tiễn vững chắc để các đơn vị lâm nghiệp triển khai nhân rộng quy trình khoanh nuôi, phục hồi tài nguyên rừng bản địa bền vững.