Giới thiệu dự án

Rừng nhiệt đới đóng vai trò là "lá phổi xanh", giữ cân bằng sinh thái, điều hòa nguồn nước, bảo tồn đa dạng sinh học và giảm thiểu tác động của biến đổi khí hậu. Theo thống kê của Tổ chức Lương thực và Nông nghiệp Liên Hợp Quốc (FAO), mỗi năm trên thế giới có khoảng 11 triệu hecta rừng bị phá hủy; riêng khu vực Châu Á - Thái Bình Dương mất khoảng 1,8 triệu ha/năm (tương đương 5.000 ha rừng nhiệt đới biến mất mỗi ngày). Tại Việt Nam, độ che phủ rừng từng sụt giảm nghiêm trọng từ 43% (năm 1943) xuống còn 26% (năm 1993) do hậu quả chiến tranh, nạn du canh du cư và tập quán khai thác kiệt quệ.

Khu vực xã La Bằng (huyện Đại Từ, tỉnh Thái Nguyên) nằm ở sườn Đông dãy núi Tam Đảo (tọa độ $21^\circ 37' 34''$ B đến $105^\circ 32' 10''$ Đ), có địa hình dốc đứng từ $18^\circ$ đến $45^\circ$, độ cao 300 – 600 m. Rừng tự nhiên tại đây chịu tác động kép: khai thác chọn gỗ tự do và phát rừng làm nương rẫy canh tác chè (với diện tích lên tới 400 ha). Hệ quả là hệ sinh thái bị thoái hóa về trạng thái rừng phục hồi IIB (rừng nghèo kiệt, cấu trúc tầng tán bị xáo trộn, trữ lượng gỗ suy giảm, thảm tươi cây bụi phát triển dày đặc chèn ép cây tái sinh).

Thực trạng suy thoái -> Xáo trộn cấu trúc (N/D1.3, N/Hvn) -> Cạnh tranh sinh dưỡng gay gắt -> Chậm diễn thế tiến hóa

Để phục hồi hệ sinh thái rừng IIB theo hướng diễn thế tiến hóa, dự án xác định các mục tiêu trọng tâm:

  1. Định lượng mật độ lâm phần ($N$), tiết diện thân ($G_i$), và chỉ số giá trị quan trọng ($IVI$) của các loài cây gỗ chủ yếu tại 9 ô tiêu chuẩn (OTC) đại diện.
  2. Mô hình hóa quy luật cấu trúc không gian ngang (quan hệ số cây/cấp đường kính - $N/D_{1.3}$) và cấu trúc không gian đứng (quan hệ số cây/cấp chiều cao - $N/H_{vn}$) bằng hàm phân bố toán học Weibull.
  3. Đánh giá mức độ che phủ và tính chất cạnh tranh của tầng cây bụi thảm tươi theo thang đo Drude.
  4. Đề xuất hệ thống giải pháp lâm sinh (khoanh nuôi xúc tiến tái sinh, tra dặm loài bản địa, tỉa thưa điều chỉnh mật độ) và cơ chế quản lý giao khoán bền vững cho địa phương.
Chỉ số mục tiêu Giá trị định lượng dự kiến Phương pháp nghiệm thu
Mật độ tầng cây gỗ ($N$) 330 – 480 cây/ha Điều tra 9 OTC ($2.000,\text{m}^2/\text{OTC}$)
Tiết diện thân lâm phần ($G$) 14,0 – 33,0 $\text{m}^2/\text{ha}$ Tổng hợp tiết diện $g_i = \frac{\pi \cdot D_{1.3}^2}{4}$
Mức độ thích ứng mô hình Weibull $100%$ OTC có $\chi^2_{tn} < \chi^2_{0.05}$ Kiểm định Pearson Chi-Square ($\alpha = 0,05$)
Tầng che phủ cây bụi thảm tươi $25% - 100%$ (Cấp Cop2 đến Soc) Thang phân cấp mật độ Drude

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 đánh giá rừng phục hồi chủ yếu dựa vào phương pháp mô tả định tính hoặc áp dụng các phân bố lý thuyết cổ điển (như phân bố Chuẩn Gauss hoặc phân bố Giảm số mũ Meyer), dẫn đến sai số lớn khi lập kế hoạch tỉa thưa và dự báo trữ lượng cho rừng tự nhiên hỗn loài đa tuổi.

Tiêu chí so sánh Phương pháp định tính Mô hình phân bố Meyer Mô hình phân bố Weibull (Giải pháp đề xuất)
Độ linh hoạt toán học Rất thấp (chỉ quan sát mắt) Thấp (chỉ mô tả đường cong giảm) Rất cao (mô tả đa dạng: lệch trái, lệch phải, chuẩn)
Độ chính xác kiểm định Không thể kiểm định Sai lệch ở cấp kính trung bình Đạt độ tin cậy $95%$ qua kiểm định $\chi^2$
Khả năng ứng dụng lâm sinh Dễ sai lệch tác động Chỉ phù hợp rừng nguyên sinh ổn định Xác định chính xác cấp kính/chiều cao cần tỉa thưa
Độ phức tạp tính toán Thấp Trung bình Yêu cầu thuật toán hồi quy phi tuyến tính

Phân tích yêu cầu lâm sinh theo mô hình MoSCoW:

  • Must Have: Xác định công thức tổ thành theo $IVI%$; mô phỏng thực nghiệm phân bố $N/D_{1.3}$ và $N/H_{vn}$ bằng hàm phân bố 3 tham số Weibull; kiểm định Pearson $\chi^2$.
  • Should Have: Phân cấp độ nhiều thảm tươi theo thang Drude trên các ô dạng bản ($25,\text{m}^2$); đo độ dốc địa hình bằng Clinometer tự chế chuẩn hóa góc nghiêng.
  • Could Have: Bản đồ hóa phân bố không gian cây rừng kết hợp tọa độ GPS và GIS.
  • Won't Have: Đo đếm sinh khối rễ ngầm và phân tích giải phẫu học từng loài gỗ trong giai đoạn này.

Thiết kế hệ thống và mô hình toán học

Cấu trúc luồng xử lý dữ liệu lâm nghiệp từ ngoại nghiệp đến mô hình hóa lý thuyết:

Hệ thống công thức toán học và giải thuật áp dụng:

  1. Quy đổi đường kính ngang ngực ($D_{1.3}$): $$D_{1.3} = \frac{P}{\pi} \quad (\text{cm})$$ (Trong đó $P$ là chu vi thân cây đo tại vị trí $1,3,\text{m}$)

  2. Chỉ số mức độ quan trọng ($IVI$): $$IVI_i (%) = \frac{A_i + D_i}{2}$$

    • Độ phong phú tương đối: $A_i (%) = \frac{N_i}{\sum N} \times 100$
    • Độ ưu thế tương đối: $D_i (%) = \frac{G_i}{\sum G} \times 100$, với $G_i = \sum \frac{\pi \cdot D_{1.3}^2}{40.000} ; (\text{m}^2)$

  3. Hàm mật độ xác suất Weibull (Weibull Probability Density Function): $$f(x) = \alpha \cdot \lambda \cdot x^{\alpha - 1} \cdot e^{-\lambda \cdot x^\alpha}$$

    • $\alpha$: Tham số hình dạng (Shape parameter)
    • $\lambda$: Tham số tỷ lệ (Scale parameter)
    • $x$: Cấp đường kính ($D_{1.3}$) hoặc cấp chiều cao ($H_{vn}$)

  4. Kiểm định tiêu chuẩn $\chi^2$ (Chi-Square Test): $$\chi^2_{tn} = \sum_{j=1}^{k} \frac{(f_{th(j)} - f_{lt(j)})^2}{f_{lt(j)}}$$ Điều kiện chấp nhận mô hình: $\chi^2_{tn} < \chi^2_{0.05}$ với bậc tự do xác định.

Methodology

Dự án triển khai theo phương pháp luận sinh thái học phát sinh quần thể thảm thực vật nhiệt đới của Thái Văn Trừng (1978):

  • Thiết kế mẫu điều tra: Lập 9 ô tiêu chuẩn ngẫu nhiên hệ thống đại diện cho 3 dạng địa hình: Đỉnh (OTC 1, 4), Sườn (OTC 2, 5, 6, 8), Chân (OTC 3, 7, 9). Kích thước mỗi OTC là $40,\text{m} \times 50,\text{m} = 2.000,\text{m}^2$ (tổng diện tích điều tra $1,8,\text{ha}$).
  • Điều tra thảm tươi, tái sinh: Thiết lập 5 ô dạng bản $5,\text{m} \times 5,\text{m} = 25,\text{m}^2$ tại 4 góc và trung tâm mỗi OTC.
  • Tiến độ triển khai: 4 tháng (từ tháng 02/2015 đến tháng 05/2015), chia làm 3 giai đoạn: Chuẩn bị/Sơ thám $\rightarrow$ Ngoại nghiệp đo đếm $\rightarrow$ Nội nghiệp phân tích mô hình.

Implementation và kết quả

Development process & Thuật toán xử lý dữ liệu

Thuật toán tự động hóa tính toán chỉ số $IVI$ và ước lượng tham số hàm Weibull được viết bằng mã Python (sử dụng thư viện scipy.stats và numpy):

import numpy as np
from scipy.stats import weibull_min, chisquare
import pandas as pd

def analyze_forest_structure(diameters, observed_counts, bin_edges):
    """
    Tính toán phân bố thực nghiệm và mô hình hóa Weibull cho cấp kính D1.3
    """
    # 1. Tính toán tần số thực nghiệm f_th
    total_trees = sum(observed_counts)
    
    # 2. Ước lượng tham số Weibull liên tục
    shape, loc, scale = weibull_min.fit(diameters, floc=0)
    
    # 3. Tính tần số lý thuyết f_lt theo từng khoảng cấp kính
    theoretical_probs = []
    for i in range(len(bin_edges) - 1):
        p = weibull_min.cdf(bin_edges[i+1], shape, loc, scale) - \
            weibull_min.cdf(bin_edges[i], shape, loc, scale)
        theoretical_probs.append(p)
    
    f_lt = np.array(theoretical_probs) * total_trees
    f_th = np.array(observed_counts)
    
    # 4. Kiểm định Chi-Square
    chi2_stat = np.sum(((f_th - f_lt) ** 2) / np.where(f_lt == 0, 1e-6, f_lt))
    
    return {
        "alpha_shape": shape,
        "scale_lambda": 1 / (scale ** shape),
        "chi2_test": chi2_stat
    }

# Triển khai tính toán cho OTC 01 (N = 80 cây / 0.2 ha = 400 cây/ha)
observed_N = [0, 2, 8, 10, 10, 15, 7, 8, 5, 14, 1]
bins = np.arange(0, 60, 5)

Testing và validation

Kết quả kiểm định độ phù hợp giữa phân bố thực nghiệm và phân bố lý thuyết Weibull trên cả 9 OTC đều đạt mức ý nghĩa $\alpha = 0,05$ ($P > 0,05$), chứng minh hàm Weibull phản ánh chính xác cấu trúc lâm phần rừng IIB:

OTC Vị trí địa hình Độ dốc ($^\circ$) Mật độ $N$ (cây/ha) Tổng tiết diện $G$ ($\text{m}^2/\text{ha}$) Tham số $\alpha$ ($N/D_{1.3}$) $\chi^2_{tn}$ ($N/D_{1.3}$) $\chi^2_{0.05}$ Đánh giá kiểm định
OTC 1 Đỉnh 45 400 32,65 2,0025 11,25 15,51 Phù hợp ($\chi^2_{tn} < \chi^2_{0.05}$)
OTC 2 Sườn 32 385 26,20 2,0032 9,84 14,07 Phù hợp ($\chi^2_{tn} < \chi^2_{0.05}$)
OTC 3 Chân 20 390 14,05 2,0045 8,76 15,51 Phù hợp ($\chi^2_{tn} < \chi^2_{0.05}$)
OTC 4 Đỉnh 38 470 31,80 2,0050 7,62 14,07 Phù hợp ($\chi^2_{tn} < \chi^2_{0.05}$)
OTC 5 Sườn 29 375 24,05 2,0057 10,12 15,51 Phù hợp ($\chi^2_{tn} < \chi^2_{0.05}$)
OTC 6 Sườn 27 335 23,95 2,0038 6,45 14,07 Phù hợp ($\chi^2_{tn} < \chi^2_{0.05}$)
OTC 7 Chân 23 455 23,70 2,0027 8,91 15,51 Phù hợp ($\chi^2_{tn} < \chi^2_{0.05}$)
OTC 8 Sườn 25 395 20,20 2,0041 7,15 14,07 Phù hợp ($\chi^2_{tn} < \chi^2_{0.05}$)
OTC 9 Chân 18 415 21,21 2,0029 6,88 12,59 Phù hợp ($\chi^2_{tn} < \chi^2_{0.05}$)

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

  1. Đặc điểm cấu trúc tổ thành và ưu thế loài:

    • Công thức tổ thành điển hình:
      • OTC 1: $19,87,\text{Dex} + 17,21,\text{Bde} + 10,48,\text{Cna} + 9,30,\text{Deg} + 7,71,\text{Son} + 35,43,\text{LK}$
      • OTC 4: $19,75,\text{Vth} + 18,34,\text{Gva} + 14,56,\text{Xnh} + 10,94,\text{Gix} + 9,65,\text{Kln} + 7,61,\text{LK}$
      • OTC 5: $30,25,\text{Gix} + 24,18,\text{Soi} + 22,35,\text{Lxe} + 17,48,\text{Bde} + 5,74,\text{LK}$
    • Các loài cây ưu thế sinh thái chiếm chỉ số $IVI$ cao gồm: Bồ đề (Styrax tonkinensis), Dẻ xanh (Lithocarpus pseudosundaicus), Giổi xanh (Michelia mediocris), Chò nâu (Dipterocarpus retusus), Xoan nhừ (Choerospondias axillaris), Kháo lá nhỏ (Machilus bonii). Đây đều là những cây gỗ lớn ưa sáng, giữ vai trò khung xương ở tầng vượt tán và tầng tán chính.

  2. Cấu trúc ngang ($N/D_{1.3}$):

    • Phân bố số cây tập trung cao nhất ở các cấp đường kính từ cấp III ($10 - 15,\text{cm}$) đến cấp VI ($25 - 30,\text{cm}$).
    • Đường cong phân bố thực nghiệm có dạng một đỉnh lệch trái tại hầu hết các OTC (OTC 1, 3, 6) và tiệm cận chuẩn tại OTC 4, 8. Các cá thể có đường kính lớn ($D_{1.3} > 45,\text{cm}$) chiếm tỷ lệ dưới $5%$, phản ánh rõ dấu vết khai thác chọn gỗ trước đây.

  3. Cấu trúc đứng ($N/H_{vn}$):

    • Lâm phần phân tầng rõ rệt:
      • Tầng ưu thế sinh thái ($H_{vn} = 10 - 20,\text{m}$): Chiếm tới $60% - 75%$ tổng số cây và số loài.
      • Tầng vượt tán ($H_{vn} > 25,\text{m}$): Số lượng cá thể suy giảm mạnh, chỉ còn một vài cây tán lớn độc lập.
      • Tầng dưới tán ($H_{vn} < 5,\text{m}$): Số cây gỗ tái sinh bị thiếu hụt nghiêm trọng do bị che bóng và cạnh tranh dinh dưỡng.

  4. Đặc điểm tầng cây bụi thảm tươi:

    • Mức độ che phủ đạt từ $25%$ đến $100%$ diện tích nền rừng (cấp Cop2 đến Soc).
    • Quần xã thảm tươi chiếm ưu thế áp đảo gồm: Cây Guột (Dicranopteris linearis), Mua (Melastoma candidum), Sim (Rhodomyrtus tomentosa), Bòng bong (Lygodium microphyllum) và các loài Dương xỉ. Độ dày rậm lớn của thảm cỏ này tạo thành lớp thảm cơ giới ngăn cản hạt giống tiếp xúc với đất khoáng, kìm hãm sự tái sinh tự nhiên của cây gỗ mục đích.

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

  • Ứng dụng toán tin định lượng trong quản lý rừng: Thay thế phương pháp ước lượng cảm quan bằng việc mô hình hóa đồng thời cả 2 trục cấu trúc ($N/D_{1.3}$ và $N/H_{vn}$) bằng hàm Weibull với độ tin cậy $95%$ ($\chi^2_{tn} < \chi^2_{0.05}$ trên $100%$ mẫu điều tra).
  • Phát hiện nút thắt sinh thái học (Ecological Bottleneck): Làm sáng tỏ nguyên nhân khiến rừng IIB La Bằng chậm phục hồi không phải do thiếu hạt giống mà do sự chèn ép cơ học và ánh sáng của thảm Guột - Mua dày đặc (mức độ Soc và Cop3 chiếm $66,7%$).
  • Cơ sở khoa học cho kỹ thuật lâm sinh: Cung cấp thông số định lượng chính xác (mật độ mục tiêu 350 - 450 cây/ha; cấp kính cần nuôi dưỡng $15 - 30,\text{cm}$) để thiết kế quy trình chặt nuôi dưỡng và làm giàu rừng tại vùng đệm Vườn Quốc gia Tam Đảo.

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

[Phát luỗng dây leo, mở tán cục bộ] ---> [Tra dặm Giổi xanh, Dẻ xanh, Chò nâu] ---> [Giao khoán đất rừng theo hộ] ---> [Diễn thế phục hồi rừng IIIA]

Giải pháp kỹ thuật lâm sinh ứng dụng

  1. Khoanh nuôi xúc tiến tái sinh có tác động:
    • Tiến hành phát dọn thảm tươi cây bụi (Guột, Mua, Bòng bong) quanh gốc cây tái sinh mục đích theo rạch hoặc băng rộng $1,5 - 2,0,\text{m}$.
    • Chặt bỏ dây leo chèn ép tán, tỉa cành đối với cây cong queo, sâu bệnh để tập trung dinh dưỡng cho cây khỏe mạnh.
  2. Làm giàu rừng bằng cây bản địa:
    • Tra dặm bổ sung các loài cây bản địa có giá trị kinh tế và sinh thái cao: Giổi xanh (Michelia mediocris), Dẻ xanh (Lithocarpus pseudosundaicus), Chò nâu (Dipterocarpus retusus) với mật độ trồng dặm $300 - 500,\text{cây/ha}$.
    • Chuyển đổi các diện tích nương chè kém hiệu quả trên các đỉnh đồi có độ dốc cao ($> 30^\circ$) sang trồng rừng phòng hộ đầu nguồn.

Giải pháp quản lý và phân cấp xã hội

  • Chuyển đổi mô hình bảo vệ chung của cộng đồng sang hình thức giao khoán rừng trực tiếp đến từng hộ gia đình gắn liền với quyền lợi khai thác lâm sản ngoài gỗ (NTFPs).
  • Tổ chức tập huấn kỹ thuật lâm sinh định kỳ cho bà con dân tộc thiểu số (Tày, Nùng, Dao, Sán Dìu) trên địa bàn xã La Bằng.

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

  • Hạn chế kỹ thuật: Số lượng ô tiêu chuẩn dừng lại ở 9 OTC ($1,8,\text{ha}$) do điều kiện địa hình chia cắt phức tạp; thời gian điều tra ngắn (4 tháng) nên chưa thu thập được dữ liệu động thái sinh trưởng nhiều năm.
  • Hướng phát triển:
    • Ứng dụng công nghệ quét viễn thám LiDAR kết hợp thiết bị bay không người lái (UAV) để lập bản đồ cấu trúc tán rừng trên quy mô toàn huyện Đại Từ.
    • Xây dựng mô hình giám sát trữ lượng Carbon rừng nhằm chuẩn bị tham gia thị trường tín chỉ Carbon quốc tế (REDD+).

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

  • Sinh viên & Nghiên cứu sinh: Bộ tư liệu mẫu chuẩn mực về phương pháp lập OTC, đo đếm lâm học và code thuật toán mô phỏng hàm phân bố Weibull trên phần mềm máy tính.
  • Kỹ sư Lâm nghiệp & Cán bộ Kiểm lâm: Quy chuẩn thông số kỹ thuật ($N, G, IVI$) để xây dựng phương án điều chế rừng, thiết kế tỉa thưa nuôi dưỡng rừng phục hồi trạng thái IIB.
  • Chính quyền địa phương & Hộ dân: Cơ sở khoa học để quy hoạch chuyển đổi đất dốc, tiếp cận chính sách chi trả dịch vụ môi trường rừng (PES), nâng cao sinh kế bền vững.

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

1. Điều kiện kỹ thuật tối thiểu để triển khai mô hình điều tra cấu trúc rừng là gì?

Cần trang bị thước dây định chuẩn đường kính, thước sào đo cao $H_{vn}$, Clinometer đo độ dốc địa hình, địa bàn cầm tay định hướng và phần mềm xử lý thống kê (Microsoft Excel, R hoặc Python) để fit hàm Weibull.

2. Giới hạn quy mô áp dụng mô hình Weibull trong cấu trúc rừng là gì?

Mô hình Weibull áp dụng tối ưu nhất cho các trạng thái rừng tự nhiên phục hồi đa tuổi (IIA, IIB, IIIA) và rừng trồng hỗn giao. Đối với rừng nguyên sinh thuần thục phân bố đa tầng phức tạp, cần kết hợp thêm mô hình phân bố hỗn hợp đa đỉnh (Weibull Mixture Model).

3. Làm sao để tích hợp kết quả nghiên cứu vào hệ thống GIS địa phương?

Tọa độ 4 góc của từng OTC được xuất ra định dạng Shapefile (.shp) hoặc GeoJSON, tích hợp vào các phần mềm GIS như QGIS/ArcGIS để số hóa bản đồ hiện trạng rừng và quản lý lớp dữ liệu thuộc tính ($D_{1.3}, H_{vn}, IVI$).

4. Chi phí bảo tồn và xúc tiến tái sinh 1 ha rừng IIB ước tính bao nhiêu?

Chi phí phát dọn cục bộ, cuốc hố tra dặm cây bản địa và bảo vệ trong 3 năm đầu ước tính dao động từ $15 - 25,\text{triệu VNĐ/ha}$, thấp hơn $60%$ so với chi phí trồng mới hoàn toàn rừng gỗ lớn.

5. Khả năng nhân rộng mô hình ra các khu vực khác như thế nào?

Phương pháp luận và thuật toán mô hình hóa này có thể chuyển giao trực tiếp cho toàn bộ diện tích rừng phòng hộ vùng núi Tam Đảo (thuộc 3 tỉnh Thái Nguyên, Vĩnh Phúc, Tuyên Quang) và các khu bảo tồn thiên nhiên vùng Đông Bắc.

Kết luận

Khóa luận tốt nghiệp của tác giả Nông Hồng Hiển đã giải quyết toàn diện bài toán đánh giá cấu trúc rừng phục hồi IIB tại xã La Bằng, huyện Đại Từ, tỉnh Thái Nguyên. Bằng việc kết hợp nhuần nhuyễn giữa phương pháp điều tra thực địa kinh điển và công cụ toán học hiện đại (mô hình hóa hàm Weibull và kiểm định Pearson $\chi^2$), nghiên cứu đã cung cấp bức tranh định lượng chi tiết về mật độ ($335 - 470,\text{cây/ha}$), tiết diện thân ($14,05 - 32,65,\text{m}^2/\text{ha}$) và tổ thành loài tầng cây gỗ. Đây là luận cứ khoa học thực tiễn vững chắc giúp địa phương ban hành các biện pháp lâm sinh khoanh nuôi, tỉa thưa và làm giàu rừng, hướng tới mục tiêu phát triển lâm nghiệp bền vững và bảo vệ môi trường sinh thái lưu vực sườn Đông Tam Đảo.