Giới thiệu dự án

Rừng tự nhiên tại Việt Nam chiếm vai trò cốt lõi trong phòng hộ đầu nguồn, chống sạt lở đất, điều hòa thủy văn và bảo tồn đa dạng sinh học. Theo số liệu thống kê ngành lâm nghiệp, độ che phủ rừng cả nước từng chạm mức báo động 28,4% vào thập niên 1980–1990 và dù đã phục hồi lên trên 41%, diện tích rừng tự nhiên giàu chỉ còn khoảng 30%, rừng trung bình 35%, phần còn lại là rừng thứ sinh nghèo kiệt và rừng phục hồi sau nương rẫy hoặc khai thác kiệt quệ. Tại huyện Võ Nhai (tỉnh Thái Nguyên), diện tích đất lâm nghiệp chiếm tới 73,57% (61.759,07 ha), nhưng phần lớn là rừng non phục hồi với trữ lượng thấp, sinh khối và cấu trúc bị đảo lộn.

Vấn đề cấp bách đặt ra là sự thiếu hụt các nghiên cứu định lượng chi tiết về đặc tính sinh thái, mật độ, cấu trúc tổ thành và quy luật đào thải tự nhiên của lớp cây tái sinh dưới tán rừng thứ sinh trên vùng núi đá vôi hiểm trở. Nếu không có cơ sở khoa học chính xác, các dự án lâm sinh dễ rơi vào tình trạng khoanh nuôi thụ động kém hiệu quả hoặc trồng dặm sai loài, gây lãng phí nguồn lực đầu tư lâm nghiệp công.

Đề tài "Nghiên cứu một số đặc điểm tái sinh tự nhiên rừng phục hồi tại xã Nghinh Tường, huyện Võ Nhai, tỉnh Thái Nguyên" (Thực hiện bởi Trần Thanh Trường, GVHD: TS. Nguyễn Thanh Tiến, Đại học Nông Lâm Thái Nguyên) được triển khai nhằm giải quyết trực tiếp bài toán trên.

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

  1. Đánh giá thực trạng cấu trúc tầng cây gỗ (mật độ, tổ thành IVI, độ tàn che) tại trạng thái rừng phục hồi IIA/IIB xã Nghinh Tường.
  2. Định lượng các chỉ số sinh thái lớp cây tái sinh: mật độ cá thể ($N/\text{ha}$), tổ thành loài, hệ số đa dạng sinh học Shannon-Weaver ($H'$), tỷ lệ cây triển vọng (CTV), chất lượng sinh trưởng và nguồn gốc tái sinh (hạt/chồi).
  3. Xác định quy luật phân bố số cây và số loài tái sinh theo các cấp chiều cao ($<1,\text{m}$, $1\text{--}<2,\text{m}$, $\ge 2,\text{m}$).
  4. Đề 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 có tác động) phù hợp với điều kiện lập địa karst vùng cao.

Phương pháp tiếp cận dựa trên quan điểm sinh thái phát sinh quần thể của Thái Văn Trừng (1978), ứng dụng hệ thống ô tiêu chuẩn (OTC) định lượng kết hợp giải thuật xử lý thống kê sinh trắc học rừng để chuyển đổi dữ liệu thực địa thành các hướng dẫn kỹ thuật lâm sinh có tính ứng dụng cao. Phạm vi nghiên cứu tập trung vào trạng thái rừng phục hồi tại xã Nghinh Tường qua 9 OTC ($1.000,\text{m}^2/\text{ô}$) và 45 ô dạng bản ($25,\text{m}^2/\text{ô}$), không đi sâu vào phân tích lý hóa tính vi lượng của phẫu diện đất.


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 phục hồi rừng tại các địa phương vùng núi phía Bắc thường áp dụng giải pháp trồng mới thuần loài (Keo, Bạch đàn) hoặc khoanh nuôi bảo vệ nghiêm ngặt không tác động. Bảng so sánh dưới đây làm rõ ưu và nhược điểm của các hướng tiếp cận:

Phương pháp tiếp cận Ưu điểm Hạn chế Mức độ phù hợp với Nghinh Tường
Trồng rừng mới nhân tạo Tăng nhanh độ che phủ tức thời Chi phí cao ($15\text{--}25\text{ triệu VNĐ/ha}$), phá vỡ hệ sinh thái bản địa đá vôi Thấp (Địa hình chia cắt, xói mòn)
Khoanh nuôi bảo vệ thuần túy Chi phí bảo vệ thấp Tốc độ phục hồi chậm, dây leo cây bụi lấn át, cây mục đích bị đào thải Trung bình (Dễ thoái hóa thành trảng cây bụi)
Xúc tiến tái sinh tự nhiên có tác động Chi phí tối ưu ($3\text{--}5\text{ triệu VNĐ/ha}$), bảo tồn loài bản địa, hệ sinh thái bền vững Cần điều tra sinh trắc chính xác trước khi ra quyết định lâm sinh Rất cao (Giải pháp tối ưu của đề tài)

Ma trận yêu cầu giải pháp theo khung MoSCoW:

  • Must Have: Xác định chính xác mật độ tầng cây cao, mật độ cây tái sinh, chỉ số IVI, tỷ lệ cây có triển vọng ($H \ge 1,0,\text{m}$, thân thẳng, không sâu bệnh) và chỉ số Shannon-Weaver ($H'$).
  • Should Have: Phân loại nguồn gốc hạt/chồi và quy luật phân tầng theo 3 cấp chiều cao thực nghiệm.
  • Could Have: Mô hình hóa toán học hàm phân bố số cây theo cấp chiều cao (Weibull/Meyer).
  • Won't Have: Phân tích vi sinh vật đất và đo bức xạ quang hợp PAR chuyên sâu từng tháng.

Thiết kế hệ thống điều tra và xử lý dữ liệu

Hệ thống phương pháp luận định lượng sinh trắc học rừng được cấu trúc thành luồng xử lý chuẩn hóa:

Bộ công cụ và công nghệ sử dụng trong nghiên cứu:

  • Phần mềm tính toán & mô hình hóa: R Studio (v4.0.2), IBM SPSS Statistics (v22.0), Microsoft Excel ToolPak (v2016).
  • Công cụ trắc địa & sinh trắc thực địa: Thước đo cao Haga/Blume-Leiss, Thước kẹp kính Mantax Blue 500mm, GPS Garmin eTrex 64s, Thước dây sợi thủy tinh 50m.
  • Tiêu chuẩn quy chuẩn kỹ thuật: TCVN 11566-1:2016 về Điều tra kiểm kê rừng; Thông tư 33/2018/TT-BNNPTNT quy định về điều tra, kiểm kê và theo dõi diễn biến rừng.

Phương pháp nghiên cứu (Methodology)

Nghiên cứu áp dụng quy trình 4 giai đoạn chuẩn:

  1. Giai đoạn 1 - Khảo sát & Thiết kế ô mẫu: Xác định 3 tuyến điều tra dài từ $7,375\text{--}8,375,\text{km}$, rải đều 9 OTC diện tích $1.000,\text{m}^2$ ($25,\text{m} \times 40,\text{m}$) đại diện cho các sườn núi đá vôi.
  2. Giai đoạn 2 - Thu thập dữ liệu thực địa: Bố trí 5 ô dạng bản ($5,\text{m} \times 5,\text{m}$) tại 4 góc và trung tâm của mỗi OTC. Thực hiện đo đường kính ngang ngực ($D_{1.3}$), chiều cao vút ngọn ($H_{vn}$), chiều cao dưới cành ($H_{dc}$), đường kính tán ($D_t$) và phẩm chất cây tái sinh.
  3. Giai đoạn 3 - Phân tích dữ liệu sinh trắc học: Xử lý các chỉ số phức hợp sinh thái theo chuẩn thống kê toán học.
  4. Giai đoạn 4 - Đánh giá rủi ro & QA: Kiểm tra chéo mẫu tiêu bản thực vật chưa rõ loài tại phòng tiêu bản thực vật Đại học Nông Lâm Thái Nguyên; kiểm chuẩn sai số đo đếm đường kính $\le 5%$, đo cao $\le 10%$.

Implementation và kết quả

Development process & Core Algorithms

Các thuật toán và công thức sinh thái định lượng cốt lõi được lập trình và tính toán tự động qua kịch bản phân tích:

  1. Chỉ số mức độ quan trọng sinh thái (Importance Value Index - IVI): $$IVI_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$; $G_i%$ là tỷ lệ phần trăm tiết diện ngang thân cây ($g_i = \frac{\pi \times D_{1.3}^2}{4}$) của loài thứ $i$ so với tổng tiết diện ngang quần xã.

  2. Chỉ số đa dạng sinh học Shannon-Weaver ($H'$): $$H' = -\sum_{i=1}^s \left( \frac{n_i}{N} \right) \ln \left( \frac{n_i}{N} \right)$$ Trong đó: $s$ là tổng số loài cây tái sinh xuất hiện; $n_i$ là số lượng cá thể loài $i$; $N$ là tổng số cá thể quan trắc trong ô dạng bản.

  3. Tỷ lệ cây tái sinh có triển vọng ($CTV%$): $$CTV(%) = \frac{\sum n_{i,(H \ge 1,0,\text{m},,\text{tốt/TB})}}{\sum N_i} \times 100$$

Mã nguồn phân tích sinh thái bằng Python (Tương đương quy trình xử lý dữ liệu của đề tài):

import numpy as np
import pandas as pd

def calculate_ecological_metrics(species_counts, tree_diameters=None):
    total_individuals = sum(species_counts.values())
    s = len(species_counts)
    
    # 1. Tính toán Shannon-Weaver Index (H')
    shannon_h = 0.0
    composition_formula = []
    for sp, count in species_counts.items():
        pi = count / total_individuals
        if pi > 0:
            shannon_h -= pi * np.log(pi)
        
        # Hệ số tổ thành (Ki) - điều kiện tham gia công thức >= 5%
        percent = (count / total_individuals) * 100
        if percent >= 5.0:
            k_val = round((count / total_individuals) * 10, 1)
            composition_formula.append(f"{k_val} {sp}")
            
    return {
        "Total_Species": s,
        "Total_N": total_individuals,
        "Shannon_H": round(shannon_h, 3),
        "Formula": " + ".join(composition_formula) + " + LK"
    }

# Dữ liệu thực nghiệm mẫu tại OTC 07 (Nghinh Tường)
sample_data = {
    "De_gai": 12, "Nhan_rung": 8, "Khao": 7, "May_puon": 6, 
    "Lat_hoa": 5, "Phay": 5, "Soi_gai": 4, "Long_mang_cut": 3, "Tau_muoi": 3
}
results = calculate_ecological_metrics(sample_data)
print(f"Shannon H': {results['Shannon_H']} | Tổ thành: {results['Formula']}")

Kết quả đạt được và kiểm chứng thực nghiệm

1. Cấu trúc tầng cây gỗ ($D_{1.3} \ge 6,\text{cm}$)

  • Mật độ cây tầng cao dao động từ $320\text{--}410\text{ cây/ha}$ (Trung bình: $356\text{ cây/ha}$).
  • Số lượng loài cây gỗ biến động từ $16\text{--}21\text{ loài/OTC}$, số loài ưu thế ($IVI \ge 5%$) từ $5\text{--}9\text{ loài}$.
  • Các loài ưu thế tầng cao chủ yếu là loài ưa sáng mọc nhanh: Dẻ gai, Nhãn rừng, Mạy tèo, Mạy puôn, Lát hoa, Gội trắng, Phay, Sấu, Lòng mang cụt, Táu muối, Dâu da xoan.

2. Mật độ, Tổ thành và Chỉ số Đa dạng sinh học cây tái sinh

Tổng hợp dữ liệu điều tra thực nghiệm trên 9 OTC ($45\text{ ô dạng bản}$):

Chỉ tiêu sinh trắc Giá trị thực nghiệm Giá trị trung bình Đánh giá chuyên môn
Mật độ tái sinh ($N/\text{ha}$) $3.120\text{--}4.400\text{ cây/ha}$ $3.884\text{ cây/ha}$ Đạt mức mật độ tái sinh tự nhiên khá theo thang phân chia Vũ Đình Huề ($>3.000\text{ cây/ha}$)
Mật độ cây triển vọng (CTV) $1.760\text{--}2.720\text{ cây/ha}$ $2.213\text{ cây/ha}$ Chiếm tỷ lệ cao từ $47,78%\text{--}61,82%$, đủ năng lực kế thừa tầng cây cao
Chỉ số Shannon-Weaver ($H'$) $2,448\text{--}2,811$ $2,642$ Phản ánh mức độ đa dạng loài và độ đồng đều cá thể cao
Nguồn gốc tái sinh Hạt: $70,82%\text{--}81,25%$
Chồi: $18,75%\text{--}29,18%$
Hạt: $76,14%$
Chồi: $23,86%$
Tái sinh hạt chiếm ưu thế tuyệt đối, tuổi thọ và độ ổn định rễ cọc vượt trội
Chất lượng cây tái sinh Tốt: $30,23%\text{--}40,91%$
TB: $35,29%\text{--}50,00%$
Xấu: $17,02%\text{--}25,00%$
Tốt: $36,25%$
TB: $42,60%$
Xấu: $21,15%$
Cây sinh trưởng tốt và trung bình chiếm $78,85%$, gấp gần 4 lần cây phẩm chất xấu

3. Quy luật phân bố cây tái sinh theo cấp chiều cao

Dữ liệu phân bố thực nghiệm thể hiện rõ quy luật đào thải tự nhiên theo cấp sinh trưởng:

Cấp chiều cao Tỷ lệ số cá thể trung bình (%) Số loài trung bình (Loài) Tỷ lệ số loài (%) Ý nghĩa sinh thái học
$< 1,0,\text{m}$ $43,02%$ $13,44$ $79,06%$ Giai đoạn mầm mạ nảy chồi nhiều nhưng chịu áp lực cạnh tranh gay gắt
$1,0\text{--}<2,0,\text{m}$ $37,94%$ $11,44$ $67,32%$ Giai đoạn vươn cao mạnh mẽ, phân hóa tầng rõ nét
$\ge 2,0,\text{m}$ $19,04%$ $8,11$ $47,72%$ Đã vượt qua tầng cỏ rác/dây leo, bước vào lớp cây triển vọng kế cận

Đường phân bố thực nghiệm số loài theo cấp chiều cao có dạng phân bố giảm dần ($79,06% \rightarrow 67,32% \rightarrow 47,72%$). Điều này chứng minh quy luật đào thải tự nhiên diễn ra rất mạnh: chỉ những cá thể và loài có biên độ sinh thái rộng, tốc độ sinh trưởng chiều cao vượt trội mới vươn lên cấp chiều cao kế cận.


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

  1. Định lượng hóa cấu trúc tái sinh trên nền địa hình đá vôi đặc thù: Đề tài cung cấp bộ dữ liệu thực chứng chi tiết về hệ số tổ thành và phân bố tầng thứ cây tái sinh tại tiểu vùng núi cao Nghinh Tường, bổ sung khoảng trống dữ liệu lâm học vùng Đông Bắc.
  2. Chứng minh năng lực tự phục hồi bằng hạt: Với tỷ lệ tái sinh hạt đạt $76,14%$ và tỷ lệ cây triển vọng đạt $56,98%$, nghiên cứu bác bỏ giả định cần phải phá băng trồng mới nhân tạo tốn kém, chứng minh tiềm năng phục hồi tự nhiên hoàn toàn khả thi nếu được can thiệp kỹ thuật đúng cách.
  3. Cải thiện hiệu quả chi phí lâm sinh: Chuyển đổi phương thức tác động từ trồng mới sang khoanh nuôi xúc tiến tái sinh có điều chỉnh giúp tiết kiệm $65\text{--}75%$ kinh phí đầu tư công cho mỗi hecta rừng phòng hộ.

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

Dựa trên kết quả định lượng cấu trúc, đề tài đề xuất gói kỹ thuật lâm sinh tác động theo 4 nhóm giải pháp:

Lộ trình triển khai thực tế (Roadmap)

  • Năm 1: Khảo sát ranh giới, thiết lập đường băng cản lửa, phát luỗng dây leo xâm lấn cục bộ quanh các cây tái sinh triển vọng ($H \ge 1,\text{m}$).
  • Năm 2: Tỉa thưa điều chỉnh tầng tán cây gỗ che bóng quá mức ($>0,7$ độ tàn che), loại bỏ cây cụt ngọn, sâu bệnh; tra dặm bổ sung hạt/cây con loài bản địa có giá trị kinh tế cao (Nghiến, Đinh mật, Lát hoa).
  • Năm 3–5: Nghiệm thu tỷ lệ thành rừng, bàn giao khoán bảo vệ cho cộng đồng thôn bản tại xã Nghinh Tường.

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

  • Hạn chế kỹ thuật: Do địa hình đồi núi đá vôi hiểm trở, nhiều vách đá dốc đứng nên việc bố trí tuyến điều tra bị giới hạn; một số loài cây bụi tầng thấp phải định danh dựa vào hình thái định tính tại hiện trường.
  • Giới hạn thời gian: Đề tài chưa theo dõi chu kỳ biến động mùa vụ (mùa khô vs mùa mưa) và chưa đo đạc thành phần cơ giới, dung tích hấp thu của phẫu diện đất.
  • Hướng phát triển tiếp theo: Ứng dụng công nghệ ảnh viễn thám không người lái (UAV LiDAR) để quét cấu trúc 3D tán rừng và lập mô hình tự động phân loại sức khỏe thảm thực vật tái sinh bằng học máy (Machine Learning).

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

  • Sinh viên & Học viên chuyên ngành Lâm nghiệp/Lâm sinh: Tiếp cận bộ dữ liệu mẫu chuẩn mực về phương pháp điều tra OTC, ô dạng bản và công thức tính toán sinh trắc học rừng ($IVI$, $H'$, $CTV$).
  • Ban Quản lý Rừng phòng hộ & Trạm Kiểm lâm Nghinh Tường: Sở hữu cơ sở khoa học chính xác để lập hồ sơ thiết kế kỹ thuật - dự toán kinh phí khoanh nuôi phục hồi rừng hàng năm.
  • Cộng đồng dân cư vùng cao huyện Võ Nhai: Hưởng lợi từ chính sách chi trả dịch vụ môi trường rừng (DVMTR) và nguồn lâm sản ngoài gỗ khi cấu trúc rừng phục hồi ổn định.

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

1. Điều kiện kỹ thuật tối thiểu để một cây tái sinh được xếp vào nhóm "Cây triển vọng" (CTV)?

Cây triển vọng bắt buộc phải thỏa mãn đồng thời 3 tiêu chí sinh trắc: Chiều cao vút ngọn $H_{vn} \ge 1,0,\text{m}$; phẩm chất đạt loại Tốt hoặc Trung bình (thân thẳng, không cụt ngọn, không bị sâu bệnh phá hoại đỉnh sinh trưởng); và thuộc nhóm loài cây mục đích hoặc cây gỗ có giá trị sinh thái/kinh tế.

2. Tại sao tái sinh hạt lại có ý nghĩa quyết định hơn tái sinh chồi trong phục hồi rừng?

Cây mọc từ hạt hình thành hệ rễ cọc ăn sâu vào các khe nứt đá vôi, giúp cây chống chịu tốt với hạn hán mùa khô và gió bão vùng núi cao. Cây chồi thường có gốc mục từ cây mẹ, tuổi thọ ngắn và dễ bị đổ gãy khi đạt kích thước lớn.

3. Tỷ lệ cây tái sinh giảm mạnh ở cấp chiều cao $\ge 2,\text{m}$ có đáng lo ngại không?

Đây là quy luật sinh thái đào thải tự nhiên bình thường ở rừng nhiệt đới. Tuy nhiên, mật độ cây $\ge 2,\text{m}$ tại Nghinh Tường vẫn đạt trung bình $19,04%$ (tương đương $\sim 740\text{ cây/ha}$), hoàn toàn đủ mật độ cấu thành tầng cây cao ưu thế trong tương lai (khi mật độ tầng cây gỗ trưởng thành chỉ cần $350\text{--}400\text{ cây/ha}$).

4. Chi phí áp dụng giải pháp khoanh nuôi xúc tiến tái sinh có tác động là bao nhiêu?

Theo định mức kinh tế kỹ thuật lâm nghiệp hiện hành, chi phí khoanh nuôi có tác động (phát dây leo, tỉa cành, trồng bổ sung) dao động từ $3,5\text{--}5,0\text{ triệu VNĐ/ha/năm}$ trong 3 năm đầu, chỉ bằng $20\text{--}25%$ so với chi phí trồng mới ($18\text{--}22\text{ triệu VNĐ/ha}$).


Kết luận

Công trình nghiên cứu của tác giả Trần Thanh Trường đã phân tích định lượng thành công cấu trúc sinh thái và quy luật tái sinh tự nhiên tại xã Nghinh Tường, huyện Võ Nhai, tỉnh Thái Nguyên. Với các phát hiện cốt lõi: mật độ cây tái sinh đạt $3.884\text{ cây/ha}$, tỷ lệ cây triển vọng đạt trên $56%$, nguồn gốc hạt chiếm $76,14%$ và chỉ số đa dạng Shannon-Weaver đạt mức cao $2,642$, nghiên cứu khẳng định tiềm năng phục hồi tự nhiên của trạng thái rừng thứ sinh nơi đây là rất lớn.

Việc áp dụng đồng bộ các giải pháp lâm sinh can thiệp nhẹ (phát luỗng dây leo, điều chỉnh độ tàn che tầng cao và bảo tồn cây mẹ gieo giống bản địa như Nghiến, Đinh mật, Lát hoa) là chìa khóa then chốt để nâng cao chất lượng rừng, bảo tồn nguồn gen quý và đảm bảo an ninh môi trường sinh thái bền vững cho vùng Đông Bắc.