Giới thiệu dự án

Cây Chò chỉ (Parashorea chinensis Wang Hsie) là loài cây gỗ lớn bản địa thuộc họ Dầu (Dipterocarpaceae), giữ vai trò xung yếu trong các hệ sinh thái rừng mưa nhiệt đới thường xanh ẩm tại Việt Nam. Theo Sách Đỏ Việt Nam, Chò chỉ được xếp vào mức độ đe dọa bậc K (loài quý hiếm, có nguy cơ tuyệt chủng) do tình trạng khai thác gỗ quá mức và suy thoái sinh cảnh tự nhiên. Mặc dù gỗ Chò chỉ có giá trị kinh tế vượt trội (tỷ trọng 0,83, xếp nhóm V, thớ mịn, chịu nước, không mối mọt), các công trình nghiên cứu định lượng về khả năng thích ứng lâm sinh và mô hình hóa sinh thái của loài này trong giai đoạn vườn ươm và rừng non còn rất hạn chế.

Vấn đề cốt lõi đặt ra là sự thiếu hụt các cơ sở dữ liệu thực nghiệm về động thái tăng trưởng cá thể ($H_{vn}$, $D_{00}$, $D_t$), quy luật phân bố cấu trúc quần thể theo thời gian, cũng như phác đồ phòng trừ sâu bệnh hại tích hợp trong điều kiện lập địa miền núi phía Bắc. Đề tài khóa luận "Thử nghiệm trồng và bước đầu đánh giá sinh trưởng phát triển cây Chò chỉ (Parashorea chinensis Wang Hsie) tại Viện Nghiên cứu và Phát triển Lâm nghiệp - Trường Đại học Nông Lâm Thái Nguyên" được triển khai nhằm thiết lập hệ số liệu thực nghiệm khoa học, phục vụ công tác bảo tồn nguồn gen và phát triển kinh doanh rừng gỗ lớn bền vững.

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

  1. Gây trồng thử nghiệm thành công 50 cây Chò chỉ mới tại khuôn viên Viện Nghiên cứu và Phát triển Lâm nghiệp (Viện NC&PTLN) và hệ thống tuyến đường nội bộ Trường Đại học Nông Lâm Thái Nguyên.
  2. Đánh giá định lượng tốc độ sinh trưởng chiều cao vút ngọn ($H_{vn}$), đường kính cổ rễ ($D_{00}$), và đường kính tán ($D_t$) qua 3 kỳ đo đếm liên tiếp trên 2 nhóm cây (trồng năm 2013 và 2015).
  3. Ứng dụng hàm phân bố thống kê toán học Weibull để mô phỏng quy luật phân bố số cây theo đường kính ($N/D_{00}$) và chiều cao ($N/H_{vn}$).
  4. Định danh thành phần sâu bệnh hại chính, tính toán mức độ hại ($R%$) và đề xuất quy trình kỹ thuật lâm sinh tác động (xử lý thực bì, bón phân, điều chỉnh không gian dinh dưỡng).

Phạm vi nghiên cứu tập trung tại xã Quyết Thắng, TP. Thái Nguyên (độ cao < 700m, lượng mưa bình quân 2.600 mm/năm, nhiệt độ $20-24^\circ\text{C}$, đất feralit đỏ vàng có hàm lượng mùn 3-7%). Giới hạn nghiên cứu là quy mô thực nghiệm cục bộ (tổng 116 cây), đánh giá biến động trong ngắn hạn từ tháng 01/2015 đến tháng 05/2015 kết hợp chuỗi dữ liệu hồi cứu từ năm 2013.


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

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

Trước đây, việc phát triển rừng trồng Chò chỉ thường gặp trở ngại do tỷ lệ nảy mầm suy giảm nhanh (>50% sức nảy mầm sau 1,5 tháng bảo quản thông thường) và cây non ưa bóng nhẹ nhưng nhạy cảm với sự chèn ép của thảm tươi tầng dưới. Các mô hình lâm sinh hiện hữu đối với cây bản địa bao gồm:

Giải pháp lâm sinh Ưu điểm Nhược điểm / Hạn chế
Trồng thuần loài đất trống Dễ cơ giới hóa, mật độ đồng đều (540 cây/ha) Cây non dễ cháy lá, xói mòn đất, nguy cơ dịch hại cao
Trồng dưới tán Keo tai tượng Có độ che bóng ban đầu Keo tai tượng hút ẩm mạnh làm đất khô cứng, Chò chỉ chết hàng loạt
Trồng rạch dưới tán rừng nghèo kiệt Bảo tồn môi trường đất, duy trì độ ẩm tự nhiên Cạnh tranh ánh sáng gay gắt, tốn công phát dọn dây leo chèn ép
Trồng thử nghiệm quy mô viện/trường (Đề tài) Kiểm soát chặt vi khí hậu, đất A giàu mùn, quản lý dịch hại vi mô Quy mô cá thể còn hẹp, đòi hỏi chăm sóc liên tục 3 năm đầu

Áp dụng phương pháp phân loại yêu cầu MoSCoW cho hệ thống kỹ thuật gây trồng Chò chỉ:

  • Must have: Bầu ươm PE $8 \times 12\text{ cm}$ ruột bầu đất tầng A (88%) + phân chuồng hoai (10%) + Supe lân (2%); xử lý ngâm hạt khử trùng bằng dung dịch $\text{KMnO}_4\ 0,1%$ trong 30 phút; giàn che sáng 50% ở giai đoạn mầm non.
  • Should have: Đào hố $40 \times 40 \times 40\text{ cm}$ bố trí so le theo hình nanh sấu; bón lót $200\text{g}$ Đạm Sunfat + $100\text{g}$ Supe lân + $100\text{g}$ Clorua kali/hố; xới váng định kỳ $10-15$ ngày/lần.
  • Could have: Bón thúc định kỳ $120, 180, 240$ ngày tuổi với dung dịch NPK tỷ lệ $3\text{N}:6\text{P}:1\text{K}$ nồng độ $0,5%$.
  • Won't have (ở giai đoạn đầu): Bón thúc đạm trước khi xuất vườn 1-2 tháng để thực hiện kỹ thuật hãm cây hóa gỗ hoàn toàn.

Thiết kế hệ thống

Hệ sinh thái kỹ thuật của đề tài tích hợp giữa quy trình canh tác lâm học tiêu chuẩn và các công cụ thống kê giải tích sinh thái học:

  • Mô hình toán học thống kê: Phân bố Weibull 2 tham số được ứng dụng để mô tả phân bố $N/D_{00}$ và $N/H_{vn}$, kết hợp kiểm định Chi-bình phương ($\chi^2$) ở mức ý nghĩa $p = 0,05$.
  • Chỉ số sinh thái đất: Độ chua $\text{pH}_{\text{KCl}} = 5 - 6$, thành phần cơ giới thịt nhẹ pha cát (sét vật lý $20 - 25%$), tầng mùn bề mặt đạt $3 - 7%$.

Methodology

Nghiên cứu áp dụng phương pháp luận Định lượng lâm học (Quantitative Silviculture) kết hợp tiếp cận "lấy không gian bù thời gian":

  • Giai đoạn 1 (T1/2015): Thiết kế lô thử nghiệm, trồng bổ sung 50 cây vụ Xuân, đo đếm ban đầu (Kỳ 1) cho toàn bộ 116 cây ($H_{vn}$ bằng sào chia vạch 0,1m; $D_{00}$ bằng thước kẹp kính; $D_t$ đo 2 hướng Đông Tây - Nam Bắc).
  • Giai đoạn 2 (T3/2015): Đo đếm định kỳ lần 2, điều tra biến động thành phần sâu hại và bệnh nấm ký sinh trên 10 lá mẫu/cây.
  • Giai đoạn 3 (T5/2015): Thu thập số liệu Kỳ 3, tính toán tăng trưởng định kỳ $Z_{nt}$, lập bảng phân cấp chất lượng, khớp hàm Weibull và kiểm định sai số $\chi^2_{\text{tn}} < \chi^2_{0,05}$.

Implementation và kết quả

Development process

Quá trình xử lý số liệu toán sinh học dựa trên các hàm giải tích tăng trưởng cây rừng:

  • Tăng trưởng thường xuyên định kỳ: $Z_{nt} = T(a) - T(a-n)$
  • Tăng trưởng bình quân chung hàng năm: $\Delta_t = \frac{T(a) - T_0}{a}$
  • Suất tăng trưởng Pressler: $P_t = \frac{2(T_a - T_{a-n})}{n(T_a + T_{a-n})} \times 100%$

Thuật toán ước lượng tham số hàm phân bố Weibull và kiểm định Chi-Square được thực thi trên môi trường phân tích dữ liệu:

import numpy as np
from scipy.stats import exponweib, chisquare

def fit_weibull_silviculture(data, num_classes=6):
    """
    Mo phong phan bo cay rung N/D hoac N/H bang ham Weibull
    va kiem dinh Chi-Binh Phuong (Chi-square test)
    """
    # Uoc luong tham so phan bo (c: shape parameter alpha, scale: beta)
    shape_alpha, loc, scale_beta = exponweib.fit(data, floc=0)
    
    # Chia khoang bien thien va tinh tan so ly thuyet
    observed_freq, bin_edges = np.histogram(data, bins=num_classes)
    expected_freq = []
    total_n = len(data)
    
    for i in range(len(bin_edges) - 1):
        x1, x2 = bin_edges[i], bin_edges[i+1]
        p = np.exp(-(x1/scale_beta)**shape_alpha) - np.exp(-(x2/scale_beta)**shape_alpha)
        expected_freq.append(p * total_n)
        
    expected_freq = np.array(expected_freq)
    # Hieu chinh tan so tranh gia tri 0
    expected_freq = np.where(expected_freq == 0, 1e-5, expected_freq)
    expected_freq = expected_freq * (total_n / np.sum(expected_freq))
    
    chi2_stat, p_val = chisquare(observed_freq, f_exp=expected_freq)
    
    return {
        "alpha_shape": shape_alpha,
        "beta_scale": scale_beta,
        "chi2_calculated": chi2_stat,
        "is_accepted": p_val > 0.05
    }

Testing và validation

Toàn bộ 116 cá thể Chò chỉ được theo dõi xuyên suốt qua 3 kỳ đo đếm thực tế (T1/2015, T3/2015, T5/2015):

BẢNG DỮ LIỆU ĐỘNG THÁI TĂNG TRƯỞNG THỰC NGHIỆM QUA 3 ĐỊNH KỲ
====================================================================================
Nhóm cá thể        Chỉ tiêu đo đếm       Lần 1 (1/2015)   Lần 2 (3/2015)   Lần 3 (5/2015)
------------------------------------------------------------------------------------
Nhóm 2013          Hvn trung bình (cm)       256              268              279
(66 cây, 3 tuổi)   D00 trung bình (cm)       4,01             4,19             4,36
                   Dt trung bình (cm)        127              132              137
------------------------------------------------------------------------------------
Nhóm 2015          Hvn trung bình (cm)        98              107              116
(50 cây trồng mới) D00 trung bình (cm)       1,45             1,59             1,75
                   Dt trung bình (cm)         54               59               65
====================================================================================

Phân tích dịch hại và công thức tính mức độ hại bình quân ($R%$): $$R% = \frac{\sum (n_i \cdot v_i)}{N \cdot V} \times 100$$ Trong đó $n_i$ là số cá thể ở cấp hại $i$; $v_i$ là trị số cấp hại ($1 \to 4$); $N$ là tổng số mẫu điều tra; $V = 4$ (cấp hại cao nhất).

Ba loài dịch hại trọng điểm được định danh:

  1. Sâu cuốn lá nhỏ (Cnaphalocrocis medinalis - Pyralidae): Bướm dài $10-12\text{ mm}$, sâu non dài $19\text{ mm}$ màu xanh lá mạ chuyển hồng khi hóa nhộng. Mức độ gây hại: Cấp II (++).
  2. Câu cấu xanh (Hypomeces squamosus - Curculionidae): Côn trùng cánh cứng dài $10-14\text{ mm}$, vòi nhai ăn khuyết mép lá non, ấu trùng ăn rễ ngầm trong đất. Mức độ gây hại: Cấp II (++).
  3. Bệnh gỉ sắt lá (Uromyces appendiculatus): Nấm ký sinh tạo đốm bột màu vàng cam mặt dưới lá, lan truyền theo gió từ tháng 9 đến tháng 5 năm sau. Mức độ gây hại: Cấp I (+).

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

  1. Chỉ số sinh trưởng cá thể:
    • Cây 3 năm tuổi (nhóm 2013) đạt tốc độ tăng trưởng bình quân chung hàng năm $\Delta_t(H_{vn}) = 70\text{ cm/năm}$, $\Delta_t(D_{00}) = 1,167\text{ cm/năm}$, $\Delta_t(D_t) = 32\text{ cm/năm}$. Tăng trưởng định kỳ 2 tháng duy trì $Z_{nt}(H_{vn}) = 11 - 12\text{ cm}$, $Z_{nt}(D_{00}) = 0,17 - 0,18\text{ cm}$.
    • Nhóm cây mới trồng 2015 thích ứng tốt với giá trị tăng trưởng định kỳ $Z_{nt}(H_{vn}) = 9\text{ cm}$, $Z_{nt}(D_{00}) = 0,14 - 0,16\text{ cm}$.
  2. Mô phỏng cấu trúc quần thể bằng hàm Weibull:
    • Phân bố số cây theo đường kính $N/D_{00}$: Tham số $\alpha = 2,46$ (dạng phân bố lệch trái), giá trị $\chi^2_{\text{tn}} < \chi^2_{0,05}$. Cỡ kính $3,05 - 5,3\text{ cm}$ chiếm $75,76%$ quần thể, chứng minh lâm phần có tính chất phân hóa đồng đều ổn định.
    • Phân bố số cây theo chiều cao $N/H_{vn}$: Tham số $\alpha = 3,98$ (dạng phân bố lệch phải), giá trị $\chi^2_{\text{tn}} < \chi^2_{0,05}$. Cỡ chiều cao $2,07 - 3,15\text{ m}$ chiếm $77,27%$, phản ánh chất lượng sinh trưởng đạt mức trung bình khá.
  3. Phân cấp chất lượng lâm phần:
    • Nhóm cây 2013: Cây tốt chiếm $26%$, cây trung bình chiếm $41%$, cây xấu chiếm $33%$.
    • Nhóm cây 2015: Cây tốt chiếm $20%$, cây trung bình chiếm $50%$, cây xấu chiếm $30%$.

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

  • Định lượng hóa cấu trúc rừng non bản địa: Thay vì dừng lại ở các nhận định định tính, đề tài đã giải tích chính xác các tham số phân bố Weibull ($\alpha = 2,46$ cho đường kính và $\alpha = 3,98$ cho chiều cao), cung cấp bộ tham số thực nghiệm đầu tiên cho Chò chỉ trồng tại vùng đồi Thái Nguyên.
  • Tối ưu hóa công thức ruột bầu: Xác lập công thức phối trộn ruột bầu đất mặt giàu mùn tầng A ($88%$) + phân chuồng hoai ($10%$) + Supe lân ($2%$) giúp tỷ lệ sống đạt trên $90%$ khi đưa ra trồng thực địa.
  • So sánh hiệu quả với các nghiên cứu tiền nhiệm:
Tiêu chí so sánh Nghiên cứu Đỗ Công Huân (2015) Mô hình Cát Bà (Trần Nguyên Giảng, 1996-1998) Mô hình Lũng Đồng Đành (Trung tâm KHSX Đông Bắc Bộ)
Loài phù trợ / Không gian Khuôn viên Viện + Trồng phân tán Trồng dưới tán Keo tai tượng & Keo lá tràm Trồng dưới tán Thông mã vĩ (10 ha)
Tỷ lệ sống / Thích ứng Ổn định > 90%, sinh trưởng $70\text{ cm/năm}$ Dưới Keo tai tượng chết; dưới Keo lá tràm sống tốt Thích nghi chọn lọc trên 5 loài bản địa
Công cụ toán mô phỏng Phân bố Weibull 2 tham số + Test $\chi^2$ Đánh giá tỷ lệ sống cơ học Điều tra mô tả sinh thái

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

  1. Kịch bản ứng dụng:
    • Trồng bổ sung làm giàu rừng tự nhiên nghèo kiệt và rừng thứ sinh sau nương rẫy ở cao trình $< 700\text{ m}$.
    • Thiết lập các băng cây bản địa hỗn giao (Chò chỉ + Sấu/Re/Trám) theo đường đồng mức với mật độ 540 cây/ha (cự ly $3,5 \times 5\text{ m}$, rạch trồng $2\text{ m}$, băng chừa $3\text{ m}$).
  2. Quy trình triển khai tiêu chuẩn:
    • Xử lý đất: Hoàn thành cuốc hố $40 \times 40 \times 40\text{ cm}$ trước khi trồng 30 ngày, bón lót $200\text{g}$ Đạm Sunfat + $100\text{g}$ Supe lân + $100\text{g}$ Clorua kali.
    • Mùa vụ tối ưu: Vụ Xuân ($10/2 - 30/3$) và vụ Hè Thu ($T6 - T8$).
    • Chăm sóc 3 năm đầu: Năm 1 và Năm 2 chăm sóc 2 lần/năm (phát dọn thực bì rạch $2\text{ m}$, xới đất đường kính $50\text{ cm}$, sâu $3-4\text{ cm}$). Năm 3 loại bỏ cây chèn ép ở băng chừa, bón thúc NPK $(5:10:3)$.
  3. Phân tích kinh tế - môi trường: Chu kỳ kinh doanh cây gỗ lớn đạt 30-40 năm, cho sản lượng gỗ lõi nhóm V giá trị cao, đồng thời cải thiện khả năng giữ nước đầu nguồn và tăng cường hấp thụ carbon so với rừng thuần loài cây ngoại lai chu kỳ ngắn.

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

  • Hạn chế kỹ thuật: Dung lượng mẫu khảo sát còn ở quy mô nhỏ ($N = 116$ cây), thời gian theo dõi định kỳ liên tục kéo dài trong 5 tháng (chưa bao quát trọn vẹn chu kỳ biến thiên 12 tháng sinh khí hậu).
  • Rào cản nguồn giống: Hạt giống Chò chỉ mất sức nảy mầm nhanh, phụ thuộc mùa quả chín theo chu kỳ sai quả ($4-6$ năm/lần), đòi hỏi hoàn thiện kỹ thuật nhân giống sinh dưỡng (cắt cành giâm hom, nuôi cấy mô tế bào).
  • Định hướng tương lai:
    1. Thiết lập các ô tiêu chuẩn định vị lâu dài (OTC $1.000\text{ m}^2$) theo dõi diễn thế cấu trúc đến giai đoạn khép tán (năm thứ 5 - 8).
    2. Ứng dụng phân tích chỉ thị phân tử DNA để đánh giá mức độ đa dạng di truyền của các dòng Chò chỉ gieo ươm.

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

  • Sinh viên & Học viên Lâm nghiệp: Cung cấp tài liệu thực chứng về phương pháp đo đếm lâm học, thuật toán phân bố Weibull và phương pháp giải tích thân cây thực địa.
  • Kỹ sư Lâm sinh & Cán bộ vườn quốc gia: Nắm bắt phác đồ gieo ươm, xử lý ruột bầu và quy trình kỹ thuật trồng Chò chỉ hỗn giao theo băng/rạch.
  • Chủ rừng & Doanh nghiệp Lâm nghiệp: Cơ sở dữ liệu để xây dựng phương án chuyển hóa rừng gỗ nhỏ sang rừng gỗ lớn bản địa đa mục đích.
  • Nhà nghiên cứu Sinh thái: Bộ số liệu định lượng về mức độ thiệt hại dịch hại ($R%$) và khả năng thích ứng đất đai của chi Parashorea tại Đông Bắc Bộ.

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

  1. Yêu cầu kỹ thuật quan trọng nhất đối với cây giống Chò chỉ khi xuất vườn là gì?
    Cây con phải đạt $12-15$ tháng tuổi, chiều cao vút ngọn $60-80\text{ cm}$, đường kính cổ rễ $D_{00} \ge 0,5-0,6\text{ cm}$, thân cây đã hóa gỗ hoàn toàn, đỉnh ngọn nguyên vẹn, không bị sâu đục nõn hay rệp sáp, và tuyệt đối không đem trồng khi cây đang ra đợt lộc/lá non.

  2. Tại sao phân bố $N/D_{00}$ của lâm phần 3 tuổi lại có dạng lệch trái ($\alpha = 2,46$)?
    Tham số hình dạng $\alpha = 2,46$ (< 3,6) thể hiện đường cong phân bố lệch sang trái, nghĩa là số lượng cây có đường kính nhỏ và trung bình ($3,05-5,3\text{ cm}$) chiếm ưu thế lớn ($75,76%$), phản ánh đặc trưng của quần thể rừng trồng non đang trong giai đoạn tích lũy sinh khối ban đầu trước khi diễn ra sự phân hóa tự nhiên mạnh mẽ.

  3. Có nên trồng Chò chỉ thuần loài trên các vùng đồi trọc, đất thoái hóa không?
    Không. Trong tự nhiên, Chò chỉ là loài cây có biên độ sinh thái hẹp, ưa đất còn tính chất đất rừng feralit ẩm, giàu mùn ($3-7%$), tầng đất dày. Trồng trên đất trống đồi trọc thoái hóa sẽ làm cây bị cháy lá, sinh trưởng còi cọc và tỷ lệ chết cao. Cần trồng cây phù trợ (Cốt khí, Đậu tràm, Keo lá tràm) cải tạo đất trước khi đưa Chò chỉ vào trồng.

  4. Biện pháp xử lý và phòng trừ sâu cuốn lá và câu cấu xanh hại Chò chỉ?
    Thường xuyên kiểm tra vườn ươm và rừng non vào sáng sớm để ngắt bỏ bao lá cuốn nhộng thủ công; làm sạch cỏ dại xung quanh gốc để phá hủy nơi đẻ trứng của câu cấu non; khi mật độ sâu hại vượt ngưỡng kiểm soát, sử dụng các chế phẩm bảo vệ thực vật sinh học (BT, nấm Beauveria) hoặc thuốc phun trừ sâu gốc cúc tổng hợp liều lượng khuyến cáo.

  5. Quy cách đào hố và mật độ trồng hỗn giao khuyến nghị là bao nhiêu?
    Hố đào kích thước chuẩn $40 \times 40 \times 40\text{ cm}$, cuốc trước khi trồng 30 ngày. Mật độ trồng hỗn giao là 540 cây/ha (cự ly $3,5 \times 5\text{ m}$), trong đó gồm 270 cây Chò chỉ xen kẽ 270 cây bản địa khác (Sấu, Re) theo tỷ lệ $1:1$ theo từng hàng so le hình nanh sấu.


Kết luận

Nghiên cứu thử nghiệm gây trồng Chò chỉ (Parashorea chinensis) tại Viện NC&PTLN - Trường ĐH Nông Lâm Thái Nguyên đã chứng minh tính khả thi sinh học và giá trị lâm sinh thực tiễn của loài cây gỗ lớn bản địa này. Với tốc độ tăng trưởng chiều cao đạt $70\text{ cm/năm}$ và đường kính cổ rễ đạt $1,167\text{ cm/năm}$ ở giai đoạn 3 năm tuổi, Chò chỉ khẳng định khả năng thích ứng tốt với lập địa đồi núi thấp miền Bắc khi được áp dụng đúng quy trình kỹ thuật.

Kết quả mô hình hóa bằng hàm Weibull và giải tích dịch hại cung cấp nền tảng toán sinh học vững chắc để các nhà hoạch định lâm nghiệp nhân rộng mô hình, đóng góp thiết thực vào chiến lược bảo tồn đa dạng sinh học và phát triển kinh tế rừng bền vững tại Việt Nam.