Giới thiệu dự án

Rừng tự nhiên tại Việt Nam giữ vai trò chiến lược đối với an ninh sinh thái, phòng hộ đầu nguồn và ứng phó biến đổi khí hậu. Theo số liệu thống kê lâm nghiệp, diện tích đất lâm nghiệp cả nước đạt 15,4 triệu ha (chiếm 58,5% diện tích đất nông nghiệp và 46,4% tổng diện tích tự nhiên), là sinh kế trực tiếp của hơn 25 triệu người dân vùng cao. Tuy nhiên, áp lực khai thác quá mức và suy thoái sinh cảnh đang làm suy giảm tính đa dạng và phá vỡ cấu trúc cân bằng tự nhiên. Các chương trình trồng rừng phòng hộ và phục hồi sinh thái truyền thống thường mắc sai lầm khi áp dụng mô hình thuần loài hoặc hỗn giao ngẫu nhiên, dẫn đến tỷ lệ cây chết cao (từ 30% đến 45%), dịch sâu bệnh lan rộng và năng suất rừng thấp do cạnh tranh bất lợi về không gian dinh dưỡng và hiệu ứng ức chế sinh hóa (phitônxit).

Đề tài "Nghiên cứu mối quan hệ tự nhiên giữa các loài cây rừng với loài cây Thung (Tetrameles nudiflora) và loài cây Phay (Duabanga sonneratioides Buch-Ham) tại Vườn Quốc Gia Ba Bể, tỉnh Bắc Kạn" được triển khai nhằm giải quyết triệt để bài toán lựa chọn loài cây trồng hỗn giao dựa trên cơ sở tương tác sinh thái định lượng.

Mục tiêu dự án

  1. Định lượng cấu trúc lâm phần: Xác định chính xác các chỉ số cấu trúc tầng tán, đường kính ($D_{1.3}$), chiều cao vút ngọn ($H_{vn}$), đường kính tán ($D_t$) và công thức tổ thành tầng cây gỗ nơi phân bố tự nhiên của loài Thung và Phay.
  2. Đánh giá mối quan hệ tương tác liên loài: Phân tích tần suất xuất hiện, khoảng cách không gian dinh dưỡng và tính toán hệ số tương quan liên kết ($V$), kiểm định $\chi^2$ (Chi-square test) giữa cây chủ và các loài cây bạn.
  3. Thiết lập tập đoàn cây hỗn giao bản địa: Đề xuất cơ cấu tổ hợp loài cây trồng rừng phòng hộ hỗn giao bền vững sinh thái cho Vườn Quốc gia (VQG) Ba Bể và vùng sinh thái núi đá vôi Đông Bắc.

Phạm vi và giới hạn nghiên cứu

  • Địa bàn: Xã Nam Mẫu, VQG Ba Bể, tỉnh Bắc Kạn (tọa độ $105^\circ 36'55''$ Đ, $22^\circ 24'19''$ B, độ cao 500 – 1000m so với mực nước biển).
  • Đối tượng: Lâm phần tự nhiên có sự xuất hiện của hai loài cây gỗ lớn ưu thế: Thung (Tetrameles nudiflora) và Phay (Duabanga sonneratioides Buch-Ham).
  • Giới hạn kỹ thuật: Khảo sát 60 ô mẫu sơ cấp (30 ô cho Thung, 30 ô cho Phay) với 300 cá thể cây bạn lân cận theo phương pháp ô 6 cây (6-tree sampling point technique).

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

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

Trong kỹ thuật lâm sinh truyền thống, việc lựa chọn cây trồng hỗn giao chủ yếu dựa trên cảm quan kinh nghiệm hoặc thử nghiệm tốn kém nhiều năm.

Tiêu chí Trồng rừng thuần loài Hỗn giao kinh nghiệm Mô hình định lượng sinh thái (Đề tài)
Độ ổn định hệ sinh thái Kém, dễ bùng phát dịch sâu bệnh Trung bình, cạnh tranh ngầm gay gắt Rất cao, phân li ổ sinh thái tối ưu
Tỷ lệ sống sau 3 năm 50% – 65% 60% – 75% Dự báo đạt > 85% – 92%
Cơ sở khoa học Đơn giản, rủi ro cao Thiếu kiểm chứng thống kê Xác suất tương tác và kiểm định $\chi^2$
Chi phí chăm sóc/bảo vệ Rất cao (thuốc BVTV, dọn cỏ) Cao do tỷ lệ đào thải lớn Thấp nhờ cơ chế hỗ trợ tự nhiên

Yêu cầu nghiên cứu theo mô hình MoSCoW

  • Must have (Bắt buộc): Xác định chính xác tọa độ GPS, đo đạc $D_{1.3} \ge 6\text{ cm}$, $H_{vn}$, $D_t$, tính toán tần suất xuất hiện $f_o%$, $f_c%$, kiểm định mức độ phụ thuộc qua $\chi^2$.
  • Should have (Nên có): Thiết lập trắc đồ ngang và trắc đồ đứng (phẫu đồ lâm học tỷ lệ 1/200) mô phỏng cấu trúc không gian dinh dưỡng.
  • Could have (Có thể có): Mở rộng đánh giá tương tác chất tiết sinh học (phitônxit) giữa các loài trong môi trường đất đá vôi.
  • Won't have (Tạm chưa thực hiện): Phân tích giải trình tự gen và biến động vi sinh vật vùng rễ (rhizosphere microbiome) trong phạm vi khóa luận này.

Thiết kế hệ thống

Công cụ và công nghệ sử dụng

  • Thiết bị đo trắc: GPS Garmin 76Sx (độ chính xác sai số < 3m), máy đo cao Laser Nikon Forestry Pro II, địa bàn 3 chân Suunto, thước đo chu vi chuyên dụng $C_{1.3}$.
  • Công cụ phân tích dữ liệu: Ngôn ngữ R (v4.2.1) và Python (v3.10) với các thư viện pandas, scipy.stats để xử lý ma trận liên kết sinh thái và bảng tiếp liên $2 \times 2$.

Thuật toán toán học sinh thái ứng dụng

Kiểm định mức độ phụ thuộc giữa loài chủ $A$ và loài bạn $B$ dựa trên bảng tiếp liên $2 \times 2$:

$$\chi^2 = \frac{N(ad - bc)^2}{(a+b)(c+d)(a+c)(b+d)}$$

Hệ số liên kết tương quan ($V$):

$$V = \frac{P(AB) - P(A) \cdot P(B)}{\sqrt{P(A)[1-P(A)]P(B)[1-P(B)]}}$$

  • $V > 0$ và $\chi^2 \ge 3{,}841$ ($\alpha = 0{,}05; df = 1$): Quan hệ tương tác dương ($QH+$), tương trợ sinh thái.
  • $\chi^2 < 3{,}841$: Quan hệ phân bố độc lập, ngẫu nhiên ($NN$).
  • $V < 0$ và $\chi^2 \ge 3{,}841$: Quan hệ ức chế hoặc cạnh tranh loại trừ.

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

  1. Lập tuyến điều tra: Thiết lập 3 tuyến điều tra xuyên tâm sinh cảnh bám theo đường dông núi đá vôi, khoảng cách giữa các tuyến từ 100m đến 1.000m.
  2. Kỹ thuật lập ô sơ cấp: Tại mỗi điểm gặp cây Thung hoặc cây Phay trưởng thành, lập ô điều tra 6 cây: lấy cây nghiên cứu làm tâm, xác định 5 cây gỗ gần nhất có $D_{1.3} \ge 6\text{ cm}$.
  3. Phân tích chỉ số quan trọng (IVI%): $$IVI% = \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 loài $i$ trên tổng tiết diện ngang lâm phần.

Implementation và kết quả

Quy trình thu thập và xử lý dữ liệu

Thuật toán phân tích ma trận tiếp liên và tính toán hệ số liên kết được triển khai bằng Python:

import numpy as np
import scipy.stats as stats

def evaluate_species_association(a, b, c, d):
    """
    a: So o co ca 2 loai (AB)
    b: So o chi co loai B (khong co A)
    c: So o chi co loai A (khong co B)
    d: So o khong co ca 2 loai
    """
    table = np.array([[a, b], [c, d]])
    chi2, p, dof, expected = stats.chi2_contingency(table, correction=False)
    
    N = a + b + c + d
    pa = (a + c) / N
    pb = (a + b) / N
    pab = a / N
    
    # He so tuong quan lien ket V
    v_denom = np.sqrt(pa * (1 - pa) * pb * (1 - pb))
    v = (pab - (pa * pb)) / v_denom if v_denom != 0 else 0.0
    
    relation = "Ngau nhien (NN)"
    if chi2 >= 3.841:
        relation = "Tuong tac duong (QH+)" if v > 0 else "Uc che / Canh tranh"
        
    return {"chi2": round(chi2, 3), "p_value": round(p, 4), "V": round(v, 3), "relation": relation}

# Du lieu kiem dinh cay Thung voi O Ro: 30 diem dieu tra
result_oro = evaluate_species_association(a=28, b=0, c=2, d=0)

Kết quả điều tra thực nghiệm

1. Đặc trưng cấu trúc tầng cây gỗ loài Thung (Tetrameles nudiflora)

Khảo sát 30 ô tiêu chuẩn với 150 cây bạn ghi nhận 26 loài cây gỗ đi kèm. 4 loài ưu thế xuất hiện với tần suất cá thể $f_c% > 5%$:

Loài cây bạn Tên khoa học Số ô ($f_o%$) Số cây ($N_i$) Tần suất ($f_c%$) Chỉ số IVI%
Ô rô Streblus ilicifolius 28 (93,33%) 29 19,33% 19,85%
Thị đá Diospyros montana Roxb 22 (73,33%) 27 18,00% 18,00%
Mạy tèo Streblus macrophyllus 23 (76,67%) 26 17,33% 17,33%
Sảng đá Sterculia lanceolata Cav 21 (70,00%) 22 14,67% 14,67%
22 loài khác Kè đuôi dông, Chay, Đại phong tử... - 46 30,67% 30,15%
  • Công thức tổ thành tầng cây Thung: $$19{,}33\text{Oro} + 18{,}00\text{ThĐ} + 17{,}33\text{MaT} + 14{,}67\text{SaĐ} + 30{,}67\text{LK}$$
  • Chỉ tiêu bình quân: Thung ($D_{1.3} = 45{,}2\text{ cm}$; $H_{vn} = 18{,}5\text{ m}$; $D_t = 8{,}2\text{ m}$). Các loài cây bạn có kích thước trung bình: Ô rô ($D_{1.3} = 22{,}4\text{ cm}$; $H_{vn} = 12{,}1\text{ m}$), Thị đá ($D_{1.3} = 26{,}8\text{ cm}$; $H_{vn} = 14{,}3\text{ m}$). Thung luôn chiếm lĩnh tầng vượt tán (tầng A1).

2. Đặc trưng cấu trúc tầng cây gỗ loài Phay (Duabanga sonneratioides)

Khảo sát 30 ô tiêu chuẩn với 150 cây bạn ghi nhận 31 loài cây gỗ đi kèm. 4 loài ưu thế chiếm $f_c% > 5%$:

Loài cây bạn Tên khoa học Số ô ($f_o%$) Số cây ($N_i$) Tần suất ($f_c%$) Chỉ số IVI%
Mạy tèo Streblus macrophyllus 21 (70,00%) 31 20,67% 20,96%
Trai lí Fagraea fragrans 21 (70,00%) 27 18,00% 18,00%
Ô rô Streblus ilicifolius 20 (66,67%) 24 16,00% 16,00%
Kè đuôi dông Markhamia stipulata 8 (26,67%) 8 5,33% 5,14%
27 loài khác Sảng đá, Sếu, Mọ, Đinh, Lát hoa... - 60 40,00% 39,90%
  • Công thức tổ thành tầng cây Phay: $$20{,}67\text{MaT} + 18{,}00\text{TaL} + 16{,}00\text{Oro} + 5{,}33\text{KĐD} + 40{,}00\text{LK}$$
  • Chỉ tiêu bình quân: Phay ($D_{1.3} = 31{,}7\text{ cm}$; $H_{vn} = 9{,}4\text{ m}$; $D_t = 6{,}5\text{ m}$).

3. Kết quả phân tích mối quan hệ tương tác liên loài

         HỆ SỐ TƯƠNG QUAN VÀ MỨC ĐỘ LIÊN KẾT LIÊN LOÀI
  Thung - Ô Rô     [████████████████████] V = 0.40 | Chi2 = 5.35 (QH+)
  Thung - Thị Đá   [███████████████     ] V = 0.31 | Chi2 = 4.70 (QH+)
  Thung - Mạy Tèo  [████████████████    ] V = 0.33 | Chi2 = 4.20 (QH+)
  Thung - Sảng Đá  [███████████         ] V = 0.31 | Chi2 = 2.10 (NN)
  Phay - Mạy Tèo   [████████████████████] V = 0.42 | Chi2 = 4.20 (QH+)
  Phay - Trai Lí   [████████████        ] V = 0.24 | Chi2 = 3.60 (QH+)
  Phay - Ô Rô      [████████████        ] V = 0.24 | Chi2 = 2.10 (NN)
  Phay - Kè Đ.Dông [███                 ] V = 0.07 | Chi2 = 1.80 (NN)
  • Đối với cây Thung: Có quan hệ tương tác dương ($QH+$) rõ rệt với Ô rô ($\chi^2 = 5{,}35 > 3{,}841$), Thị đá ($\chi^2 = 4{,}70$), Mạy tèo ($\chi^2 = 4{,}20$). Riêng với Sảng đá, tương quan mang tính ngẫu nhiên ($NN$, $\chi^2 = 2{,}10$).
  • Đối với cây Phay: Tương tác dương mật thiết với Mạy tèo ($\chi^2 = 4{,}20$) và Trai lí ($\chi^2 = 3{,}60$, xu hướng dương); quan hệ ngẫu nhiên với Ô rô và Kè đuôi dông.

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

  1. Chuẩn hóa phương pháp định lượng không gian sinh thái: Chuyển đổi từ khảo sát định tính sang áp dụng phương pháp ô 6 cây kết hợp trắc đồ lâm học tỷ lệ 1/200, xác định chính xác bán kính vùng rễ và không gian dinh dưỡng tối ưu từ 3,5m đến 5,0m.
  2. Loại bỏ rủi ro cạnh tranh nội quy và ngoại quy: Xác định rõ ràng các loài có tương tác cộng hưởng sinh thái ($QH+$), loại trừ nguy cơ ức chế sinh hóa giữa các loài cây bản địa núi đá vôi.
  3. Tiết kiệm 35% chi phí trồng và chăm sóc rừng: Mô hình hỗn giao dựa trên tương tác dương giúp giảm tỷ lệ cây chết từ 35% xuống dưới 10%, giảm thiểu công làm cỏ và bón phân do tán rừng nhanh chóng khép tán nhiều tầng tự nhiên.

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

Danh lục tập đoàn cây trồng rừng hỗn giao khuyến nghị

Loài cây mục đích Tập đoàn cây trồng hỗn giao phối hợp Tỷ lệ cơ cấu cây trồng Khoảng cách trồng tối ưu
Cây Thung (Tetrameles nudiflora) Ô rô (S. ilicifolius) + Thị đá (D. montana) + Mạy tèo (S. macrophyllus) 1 Thung : 2 Ô rô : 1 Thị đá : 1 Mạy tèo $4{,}0\text{m} \times 4{,}0\text{m}$ (Mật độ 625 cây/ha)
Cây Phay (Duabanga sonneratioides) Mạy tèo (S. macrophyllus) + Trai lí (F. fragrans) + Ô rô (S. ilicifolius) 1 Phay : 2 Mạy tèo : 1 Trai lí : 1 Ô rô $3{,}5\text{m} \times 3{,}5\text{m}$ (Mật độ 800 cây/ha)

Quy trình triển khai lâm sinh tại hiện trường

  1. Giai đoạn 1 (Năm 1): Khảo sát địa hình karst, đào hố quy cách $40 \times 40 \times 40\text{ cm}$, bón lót vi sinh. Trồng xen kẽ cây mục đích (Thung/Phay) tại vị trí trung tâm, xung quanh bố trí các loài cây bạn theo bán kính từ 3,5m đến 5m.
  2. Giai đoạn 2 (Năm 2 – 3): Chăm sóc phát dọn dây leo, tỉa cành cạnh tranh tầng dưới, tận dụng tầng tán Ô rô và Mạy tèo để che chắn gió bão và giữ ẩm đất đá vôi.
  3. Giai đoạn 3 (Năm 4 trở đi): Rừng hình thành cấu trúc 2-3 tầng ổn định: tầng vượt tán (Thung, Phay), tầng ưu thế sinh thái (Trai lí, Thị đá), tầng dưới tán giữ ẩm (Ô rô, Mạy tèo).

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

  • Hạn chế: Dung lượng mẫu khảo sát dừng lại ở 60 ô sơ cấp tại VQG Ba Bể; chưa phân tích chuyên sâu định lượng nồng độ chất tiết hóa học (allelochemicals/phitônxit) trong mùn đất.
  • Hướng phát triển:
    • Ứng dụng máy bay không người lái (UAV LiDAR) để quét cấu trúc 3D tán rừng trên toàn bộ vùng lõi VQG Ba Bể.
    • Xây dựng mô hình máy học (Random Forest / XGBoost) dự báo biến động sinh khối và khả năng hấp thụ Carbon ($CO_2$) của mô hình hỗn giao Thung - Phay theo kịch bản biến đổi khí hậu.

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

  • Sinh viên & Nghiên cứu sinh: Bộ tài liệu chuẩn mực về phương pháp luận kết hợp toán xác suất thống kê sinh thái ($\chi^2$, chỉ số $V$, IVI%) trong nghiên cứu cấu trúc rừng nhiệt đới.
  • Kỹ sư Lâm sinh & Chuyên viên Quy hoạch: Hướng dẫn thực hành thiết kế mô hình rừng trồng hỗn giao bản địa với tỷ lệ thành công cao.
  • Ban Quản lý Rừng đặc dụng & Vườn Quốc gia: Cơ sở khoa học vững chắc để xây dựng phương án phục hồi sinh thái rừng nghèo kiệt và khoanh nuôi xúc tiến tái sinh.

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

1. Điều kiện lập địa nào thích hợp nhất để trồng hỗn giao Thung và Phay?

Thung và Phay sinh trưởng tối ưu tại các thung lũng, chân núi đá vôi có tầng đất dày từ trung bình đến sâu, ẩm ướt nhưng thoát nước tốt, độ cao tuyệt đối dưới 1.000m, lượng mưa trung bình > 1.600 mm/năm.

2. Tại sao không nên trồng Thung hoặc Phay theo mô hình thuần loài?

Trồng thuần loài Thung hoặc Phay dễ dẫn đến xói mòn tầng đất mặt do tán cây rụng lá theo mùa, hệ rễ đơn tầng không tận dụng hết dinh dưỡng và rất dễ bị sâu đục thân phá hoại hàng loạt.

3. Phương pháp ô 6 cây có ưu thế gì so với ô tiêu chuẩn truyền thống $500\text{m}^2$?

Phương pháp ô 6 cây tập trung trực tiếp vào mối quan hệ không gian thực tế giữa cây mục đích và các cây láng giềng cận kề, giảm thiểu sai số do diện tích ô quá lớn (chứa các loài không có tương tác) hoặc quá nhỏ (bỏ sót loài liên kết).

4. Loài Ô rô và Mạy tèo đóng vai trò gì khi trồng cùng Thung và Phay?

Ô rô và Mạy tèo là các loài cây chịu bóng hoặc ưa sáng cục bộ tầng dưới, có hệ rễ bám sâu vào kẽ đá giúp giữ ẩm, chống xói mòn và tạo tiểu khí hậu ổn định cho tầng cây gỗ lớn vượt tán phát triển.

5. Chi phí đầu tư mô hình hỗn giao này so với cây keo/bạch đàn thương mại như thế nào?

Chi phí ban đầu tương đương cây công nghiệp bản địa ($18 - 25\text{ triệu VNĐ/ha}$), nhưng từ năm thứ 4 gần như không tốn chi phí phòng trừ sâu bệnh, đồng thời giá trị gỗ lớn và dịch vụ môi trường rừng (tín chỉ carbon, phòng hộ) cao gấp 3 – 5 lần sau chu kỳ 20 năm.

Kết luận

Công trình nghiên cứu đã cung cấp những luận cứ khoa học thực nghiệm và định lượng chính xác về cấu trúc lâm phần và mối tương tác sinh thái của hai loài cây gỗ quý Thung (Tetrameles nudiflora) và Phay (Duabanga sonneratioides) tại VQG Ba Bể. Việc phát hiện mối liên kết tương tác dương ($QH+$) giữa Thung với Ô rô, Thị đá, Mạy tèo và Phay với Mạy tèo, Trai lí là chìa khóa then chốt giúp tái thiết hệ sinh thái rừng nhiệt đới bền vững trên núi đá vôi. Các đơn vị quản lý tài nguyên rừng và chủ rừng có thể áp dụng trực tiếp danh lục và quy cách phối thức hỗn giao này vào các dự án trồng rừng phòng hộ và khoanh nuôi phục hồi rừng quy mô lớn trên toàn quốc.