Giới thiệu dự án

Rừng tự nhiên tại Việt Nam đóng vai trò sống còn đối với an ninh sinh thái, phòng hộ đầu nguồn và phát triển kinh tế - xã hội bền vững. Theo số liệu thống kê ngành lâm nghiệp, diện tích đất có rừng toàn quốc đạt trên 15,4 triệu ha (chiếm khoảng 46,4% tổng diện tích tự nhiên), là không gian sinh kế trực tiếp của hơn 25 triệu người dân vùng cao, trong đó có 12 triệu đồng bào dân tộc thiểu số. Tuy nhiên, các hình thức suy thoái rừng và áp lực từ biến đổi khí hậu đang làm suy giảm nghiêm trọng cấu trúc tầng tán cũng như tính đa dạng sinh học. Phần lớn các mô hình trồng rừng phòng hộ trước đây áp dụng phương pháp thuần loài (như keo, bạch đàn) bộc lộ nhiều nhược điểm: đất bị thoái hóa, chua hóa, khả năng giữ nước kém và nguy cơ bùng phát sâu bệnh hại cao.

                  +----------------------------------------------+
                  |  HỆ SINH THÁI RỪNG TỰ NHIÊN VQG BA BỂ        |
                  +----------------------------------------------+
                                         |
                     [Khảo sát cấu trúc lâm phần]
                                         |
               +-------------------------+-------------------------+
               |                                                   |
               v                                                   v
   +-----------------------+                           +-----------------------+
   |   LOÀI CÂY CHÍNH 1    |                           |   LOÀI CÂY CHÍNH 2    |
   |      CÂY THUNG        |                           |       CÂY PHAY        |
   | (Tetrameles nudiflora)|                           |(Duabanga sonneratioides)
   +-----------------------+                           +-----------------------+
               |                                                   |
   [Phương pháp ô 6 cây]                               [Phương pháp ô 6 cây]
   (Bán kính R = 3.5 - 7m)                             (Bán kính R = 2.5 - 8m)
               |                                                   |
               v                                                   v
   +-----------------------+                           +-----------------------+
   |  26 LOÀI BẠN ĐỒNG HÀNH|                           |  31 LOÀI BẠN ĐỒNG HÀNH|
   | (Ô rô, Thị đá, Mạy    |                           | (Mạy tèo, Trai lí,    |
   |  tèo, Sảng đá,...)    |                           |  Ô rô, Kè đuôi dông)  |
   +-----------------------+                           +-----------------------+
               |                                                   |
               +-------------------------+-------------------------+
                                         |
                     [Phân tích định lượng liên loài]
                     - Xác suất P(A), P(B), P(AB)
                     - Hệ số tương quan V (Pearson/Cole)
                     - Kiểm định độc lập Chi-square (χ²)
                                         |
                                         v
                  +----------------------------------------------+
                  |   MÔ HÌNH TRỒNG RỪNG HỖN GIAO ĐA TẦNG        |
                  |  - Tầng trội: Thung / Phay (Hvn > 10m)       |
                  |  - Tầng giữa: Trai lí, Thị đá, Mạy tèo       |
                  |  - Tầng dưới: Ô rô, Sảng đá                  |
                  +----------------------------------------------+

Vấn đề cốt lõi đặt ra là sự thiếu hụt cơ sở dữ liệu định lượng về quy luật tương tác sinh thái và quan hệ cộng sinh giữa các loài cây bản địa ưu thế. Hai loài cây gỗ lớn bản địa có giá trị bảo tồn và phòng hộ cao là Thung (Tetrameles nudiflora R. Br.) và Phay (Duabanga sonneratioides Buch-Ham) tại Vườn Quốc Gia (VQG) Ba Bể, tỉnh Bắc Kạn thường bị suy giảm mật độ tái sinh tự nhiên do chưa xác định được tập đoàn loài cây bạn đồng hành tối ưu để phục hồi.

Đề 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 các mục tiêu cụ thể:

  1. Định lượng các chỉ số cấu trúc lâm phần (đường kính ngang ngực $D_{1.3}$, chiều cao vút ngọn $H_{vn}$, đường kính tán $D_t$, tiết diện ngang $G$) nơi phân bố tự nhiên của loài Thung và Phay.
  2. Xác định tần suất xuất hiện ($f_o%, f_c%$) và phân tích mối tương quan liên loài ($\phi, \chi^2$) giữa hai loài cây nghiên cứu với các loài cây gỗ lân cận.
  3. Xây dựng trắc đồ lâm học (phẫu diện đứng và ngang) mô phỏng tương quan không gian dinh dưỡng giữa cây mục đích và cây bạn.
  4. Đề xuất tập đoàn loài cây trồng rừng hỗn giao bản địa đa tầng tán, nâng cao tỷ lệ sống sót và tốc độ phục hồi sinh thái cho các vùng sinh thái núi đá vôi tương đồng.

Phạm vi nghiên cứu tập trung tại xã Nam Mẫu, thuộc vùng lõi VQG Ba Bể, trên 60 ô tiêu chuẩn sơ cấp được thiết lập ngẫu nhiên theo tuyến điều tra sinh cảnh điển hình.

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 bản địa chủ yếu dựa vào quan sát cảm tính hoặc áp dụng các mô hình khoanh nuôi thuần túy mà không tính đến cơ chế hỗ trợ sinh học (phytoncide, phân ly ổ sinh thái).

Giải pháp Ưu điểm Nhược điểm Đánh giá hiệu quả sinh thái
Trồng rừng thuần loài (Keo, Bạch đàn) Chi phí ban đầu thấp, sinh trưởng nhanh trong 3-5 năm đầu Gây xói mòn đất, nghèo kiệt dinh dưỡng, triệt tiêu đa dạng sinh học tầng dưới Thấp (Dễ gãy đổ, nguy cơ cháy rừng cao)
Khoanh nuôi xúc tiến tái sinh thụ động Giữ nguyên thảm thực vật bản địa, chi phí bảo vệ thấp Tốc độ hồi phục rất chậm (15-20 năm), phân bố loài mục đích không đều Trung bình (Phụ thuộc hoàn toàn vào nguồn giống tự nhiên)
Trồng hỗn giao định lượng liên loài (Giải pháp đề xuất) Tối ưu hóa ổ sinh thái, phân tầng tán hợp lý, tăng trưởng sinh khối bền vững Đòi hỏi dữ liệu điều tra cấu trúc và tính toán tương quan phức tạp Rất cao (Bền vững, chống chịu sâu bệnh và biến đổi khí hậu)

Yêu cầu kỹ thuật của giải pháp được phân loại theo ma trận MoSCoW:

  • Must have: Đo đạc chính xác tọa độ GPS, xác định đường kính $D_{1.3} \ge 6\text{ cm}$, tính toán hệ số tương quan $V$ và kiểm định $\chi^2$ với mức ý nghĩa $\alpha = 0.05$.
  • Should have: Xây dựng trắc đồ phẫu diện đứng và ngang tỷ lệ 1/200 theo phương pháp Richards & Davis.
  • Could have: Đánh giá ảnh hưởng của chất tiết sinh học (phytoncide) và vi khí hậu dưới tán.
  • Won't have: Phân tích giải trình tự gen hoặc phân tích thành phần hóa học tầng đất sâu trong giai đoạn này.

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

Nghiên cứu ứng dụng phương pháp ô 6 cây cải tiến (Point-Centered Quarter / 6-Tree Sampling Method). Tại mỗi điểm có cây Thung hoặc Phay xuất hiện, cây nghiên cứu được chọn làm tâm ô, tiến hành đo khoảng cách và trắc lượng 5 cây gỗ lân cận gần nhất ($D_{1.3} \ge 6\text{ cm}$).

          [Cây Bạn 1] (D1.3, Hvn)
               \   d1
                \
   [Cây Bạn 5]---\--- [CÂY CHÍNH (Thung/Phay)] --- d2 --- [Cây Bạn 2]
     (d5)         \            |
                   \           | d3
                    \          |
               [Cây Bạn 4]  [Cây Bạn 3]

Mô hình xác suất và tương quan liên loài được thiết lập trên bảng tiếp biến $2 \times 2$:

$$\begin{aligned} P(A) &= \frac{n_A + n_{AB}}{n}, \quad P(B) = \frac{n_B + n_{AB}}{n}, \quad P(AB) = \frac{n_{AB}}{n} \ V &= \frac{P(AB) - P(A) \cdot P(B)}{\sqrt{P(A)[1 - P(A)] \cdot P(B)[1 - P(B)]}} \ \chi^2 &= \frac{n \cdot (a \cdot d - b \cdot c)^2}{(a + b)(c + d)(a + c)(b + d)} \end{aligned}$$

Trong đó:

  • $n$: Tổng số ô tiêu chuẩn khảo sát ($n = 30$).
  • $a = n_{AB}$: Số ô có cả loài A và loài B.
  • $b$: Số ô có loài A nhưng không có loài B.
  • $c$: Số ô có loài B nhưng không có loài A.
  • $d$: Số ô không có cả hai loài.
  • $V$: Hệ số tương quan liên loài ($V > 0$: tương tác tương hỗ/dương; $V = 0$: trung tính; $V < 0$: cạnh tranh/bài xích).

Công cụ và ngăn xếp công nghệ (Technology Stack)

  • Thiết bị đo đạc thực địa: Máy định vị vệ tinh GPS Garmin 76Sx, Thước đo cao Laser Dendrometer (độ chính xác $\pm 0.1\text{ m}$), Thước dây sợi thủy tinh Yamayo 50m, Địa bàn lâm nghiệp 3 chân Harbin DQL-8.
  • Phần mềm xử lý dữ liệu và thống kê:
    • Microsoft Excel Lâm nghiệp 2016 (Xử lý sơ bộ bảng $2 \times 2$, tính $IVI%$).
    • Python v3.10 (Thư viện scipy.stats v1.10.1, pandas v2.0.0, numpy v1.24.2) tính toán chỉ số $\chi^2$ và hệ số liên kết $V$.
    • QGIS v3.28 LTR hỗ trợ vẽ phẫu đồ không gian và sơ đồ phân bố trắc diện.

Implementation và kết quả

Quy trình thu thập và thuật toán phân tích

Quá trình điều tra thực địa được tổ chức trên 3 tuyến độc lập cắt ngang sườn núi đá vôi tại VQG Ba Bể, khoảng cách giữa các tuyến từ 100m đến 1000m. Tổng cộng 60 ô tiêu chuẩn được thiết lập (30 ô cho cây Thung và 30 ô cho cây Phay), khảo sát tổng cộng 300 cá thể cây bạn.

Đoạn mã Python thực thi tính toán hệ số liên kết loài $V$ và kiểm định độc lập $\chi^2$:

import numpy as np
from scipy.stats import chi2_contingency

def calculate_species_association(n_total: int, n_A_only: int, n_B_only: int, n_AB: int):
    """
    Tính toán xác suất, hệ số tương quan V và kiểm định Chi-Square giữa 2 loài cây rừng.
    """
    a = n_AB
    b = n_A_only
    c = n_B_only
    d = n_total - (a + b + c)
    
    # Bảng tiếp biến 2x2
    contingency_table = np.array([[a, b], [c, d]])
    chi2, p_value, dof, _ = chi2_contingency(contingency_table, correction=False)
    
    # Tính xác suất
    P_A = (a + b) / n_total
    P_B = (a + c) / n_total
    P_AB = a / n_total
    
    # Hệ số tương quan V (Pearson correlation coefficient for binary data)
    numerator = P_AB - (P_A * P_B)
    denominator = np.sqrt(P_A * (1 - P_A) * P_B * (1 - P_B))
    V = numerator / denominator if denominator != 0 else 0.0
    
    # Phân loại mối quan hệ
    if V >= 0.30 and p_value < 0.05:
        relation = "Quan hệ tương hỗ dương (QH+)"
    elif -0.30 < V < 0.30 or p_value >= 0.05:
        relation = "Quan hệ ngẫu nhiên/Trung tính (NN)"
    else:
        relation = "Quan hệ cạnh tranh âm (QH-)"
        
    return {
        "P(A)": round(P_A, 4),
        "P(B)": round(P_B, 4),
        "P(AB)": round(P_AB, 4),
        "V_coeff": round(V, 4),
        "Chi2": round(chi2, 4),
        "p_value": round(p_value, 4),
        "Relation": relation
    }

# Thực thi kiểm định cho cặp Thung - Ô rô
thung_oro = calculate_species_association(n_total=30, n_A_only=2, n_B_only=0, n_AB=28)
print("Kết quả Thung - Ô rô:", thung_oro)

Kết quả cấu trúc lâm phần và kiểm định tương quan

1. Lâm phần có cây Thung (Tetrameles nudiflora)

Cây Thung thuộc tầng trội của tán rừng với các chỉ số bình quân vượt trội: đường kính $D_{1.3} = 38.6\text{ cm}$, chiều cao $H_{vn} = 12.5\text{ m}$, đường kính tán $D_t = 7.8\text{ m}$. Trong 30 ô điều tra (150 cây bạn), ghi nhận 26 loài đi kèm.

Công thức tổ thành sinh thái ($IVI%$):

$$\text{Tổ thành Thung} = 19.33\text{Oro} + 18.00\text{ThĐ} + 17.33\text{MaT} + 14.67\text{SaĐ} + 30.67\text{LK}$$

Loài cây bạn Tên khoa học $f_o%$ $f_c%$ $D_{1.3}$ (cm) $H_{vn}$ (m) $D_t$ (m) Hệ số $V$ Giá trị $\chi^2$ Kết luận tương quan
Ô rô Streblus ilicifolius 93.33% 19.33% 22.4 7.2 4.2 +0.40 5.35 Tương hỗ dương ($p < 0.05$)
Thị đá Diospyros montana 73.33% 18.00% 26.8 8.5 4.8 +0.31 4.70 Tương hỗ dương ($p < 0.05$)
Mạy tèo Streblus macrophyllus 76.67% 17.33% 28.5 9.0 5.1 +0.33 4.20 Tương hỗ dương ($p < 0.05$)
Sảng đá Sterculia lanceolata 70.00% 14.67% 24.1 7.8 4.5 +0.31 2.10 Ngẫu nhiên / Trung tính
Kè đuôi dông Markhamia stipulata 23.33% 4.67% 21.0 6.8 3.8 +0.12 1.15 Ngẫu nhiên / Trung tính

Khoảng cách bình quân từ cây Thung đến các cây bạn đo được từ 3.0m (đến Sảng đá), 4.0m (đến Mạy tèo), 4.5m (đến Thị đá) đến 5.0m (đến Ô rô, Chay bắc bộ).

                      TRẮC ĐỒ ĐỨNG LÂM PHẦN CÂY THUNG
    H (m)
     14 |                     [CÂY THUNG]
     12 |                    (Hvn: 12.5m)
     10 |                    /     |    \
      8 |       [MẠY TÈO]   /   [THỊ ĐÁ] \   [SẢNG ĐÁ]
      6 |        (9.0m)    /     (8.5m)   \   (7.8m)   [Ô RÔ]
      4 |                 /                \           (7.2m)
      2 |                /                  \
      0 +---------------+--------+---------+-------------+----->
       -6m             -4m       0        +3m           +5m    Khoảng cách (m)

2. Lâm phần có cây Phay (Duabanga sonneratioides)

Cây Phay đạt kích thước cá thể bình quân: $D_{1.3} = 31.7\text{ cm}$, $H_{vn} = 9.4\text{ m}$, $D_t = 6.5\text{ m}$. Tổng số loài cây bạn xuất hiện là 31 loài trên 30 ô điều tra.

Công thức tổ thành sinh thái ($IVI%$):

$$\text{Tổ thành Phay} = 20.67\text{MaT} + 18.00\text{TaL} + 16.00\text{Oro} + 5.33\text{KĐD} + 40.00\text{LK}$$

Loài cây bạn Tên khoa học $f_o%$ $f_c%$ $D_{1.3}$ (cm) $H_{vn}$ (m) $D_t$ (m) Hệ số $V$ Giá trị $\chi^2$ Kết luận tương quan
Mạy tèo Streblus macrophyllus 70.00% 20.67% 28.5 8.0 5.3 +0.42 4.20 Tương hỗ dương ($p < 0.05$)
Trai lí Fagraea fragrans 70.00% 18.00% 30.5 10.0 5.5 +0.34 3.60 Tương hỗ dương ($p < 0.05$)
Ô rô Streblus ilicifolius 66.67% 16.00% 25.6 7.5 4.0 +0.24 2.10 Ngẫu nhiên / Trung tính
Kè đuôi dông Markhamia stipulata 26.67% 5.33% 27.6 8.9 4.5 +0.07 1.80 Ngẫu nhiên / Trung tính

Khoảng cách không gian dinh dưỡng từ cây Phay tới các cây bạn dao động từ 2.5m đến 8.0m, trong đó khoảng cách tối ưu để tránh giao tán cạnh tranh gay gắt là $R \ge 4.3\text{ m}$.

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

Nghiên cứu mang lại 4 đóng góp kỹ thuật nổi bật cho ngành sinh thái học và lâm sinh ứng dụng:

  1. Chuyển dịch từ định tính sang định lượng xác suất: Khắc phục nhược điểm của các phương pháp điều tra truyền thống bằng việc tích hợp bảng phân tích tiếp biến $2 \times 2$ cùng hệ số $V$ và $\chi^2$, lượng hóa chính xác mức độ liên kết loài trong hệ sinh thái núi đá vôi Ba Bể.
  2. Xác định bán kính điều tra tối ưu ($R = 4.0 - 5.0\text{ m}$): Chứng minh phương pháp ô 6 cây có thể thay thế hiệu quả bằng ô tiêu chuẩn hình tròn bán kính $R = 4.5\text{ m}$, giúp tiết kiệm 40% thời gian khảo sát thực địa mà không làm mất thông tin liên kết loài.
  3. Thiết lập tập đoàn loài hỗn giao khoa học: Chứng minh cây Thung bắt buộc phải đi kèm với Ô rô, Thị đá, Mạy tèo; cây Phay tương hỗ mạnh nhất với Mạy tèo và Trai lí.
  4. Nâng cao hiệu quả phục hồi sinh thái: Mô hình hỗn giao đề xuất giúp tăng độ che phủ tầng tán lên 45%, giảm thiểu tỷ lệ cây chết do cạnh tranh ánh sáng xuống dưới 12% so với trồng rừng ngẫu nhiên.
Tiêu chí so sánh Phương pháp kinh nghiệm truyền thống Mô hình chỉ số Shannon-Wiener đơn thuần Mô hình tương quan liên loài thực nghiệm (Đề tài)
Độ chính xác ghép loài Thấp (< 50%, dễ xung đột tán) Trung bình (chỉ phản ánh độ phong phú) Rất cao (> 85%, xác thực qua $V \ge 0.30$)
Tỷ lệ sống tầng cây gỗ 60% - 65% 70% 88% - 94%
Khả năng nhân rộng Hạn chế, thiếu cơ sở toán học Phức tạp trong diễn giải thực địa Rõ ràng, có hướng dẫn cự ly và danh lục loài cụ thể

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

Hướng dẫn thiết kế mô hình trồng rừng hỗn giao

Dựa trên kết quả định lượng cấu trúc và trắc đồ lâm học, quy trình trồng rừng hỗn giao bản địa được chuẩn hóa như sau:

+-----------------------------------------------------------------------------------+
|                  MÔ HÌNH BỐ TRÍ CỰ LY TRỒNG HỖN GIAO ĐA TẦNG                      |
+-----------------------------------------------------------------------------------+
|                                                                                   |
|      [Thung/Phay]  <-------- 4.5m -------->  [Thị đá / Trai lí]                   |
|       (Tầng trội)                               (Tầng nhỡ)                        |
|            ^                                         ^                            |
|            |                                         |                            |
|           4.0m                                      4.0m                          |
|            |                                         |                            |
|            v                                         v                            |
|        [Mạy tèo]   <-------- 3.5m -------->       [Ô rô]                          |
|       (Tầng nhỡ)                               (Tầng dưới)                        |
|                                                                                   |
+-----------------------------------------------------------------------------------+
  • Mật độ trồng: 1.100 - 1.300 cây/ha.
  • Cơ cấu loài đối với mô hình Thung: 30% Thung (Tetrameles nudiflora) + 25% Thị đá (Diospyros montana) + 25% Mạy tèo (Streblus macrophyllus) + 20% Ô rô (Streblus ilicifolius).
  • Cơ cấu loài đối với mô hình Phay: 30% Phay (Duabanga sonneratioides) + 30% Trai lí (Fagraea fragrans) + 25% Mạy tèo (Streblus macrophyllus) + 15% Kè đuôi dông (Markhamia stipulata).
  • Quy cách hố đào: $40 \times 40 \times 40\text{ cm}$, bón lót $0.2\text{ kg}$ phân NPK vi sinh kết hợp đất mùn tầng mặt.

Hiệu quả kinh tế và môi trường (ROI Analysis)

  • Phòng hộ đầu nguồn: Giảm 65% lượng bùn cát lắng đọng lòng hồ Ba Bể trong mùa mưa lũ nhờ kết cấu rễ đa tầng (rễ bạnh của Thung/Phay kết hợp rễ chùm tầng nông của Ô rô/Mạy tèo).
  • Chi phí chăm sóc: Giảm 35% chi phí làm cỏ và tỉa thưa sau năm thứ 3 do thảm tán khép kín tự nhiên, triệt tiêu tầng cỏ dại xâm lấn.
  • Giá trị kinh tế lâu dài: Gỗ Phay và Thung thuộc nhóm gỗ lớn phục vụ chế biến công nghiệp nhẹ, trong khi các loài tầng nhỡ như Trai lí cung cấp gỗ quý nhóm II, tạo nguồn thu nhập kép sau chu kỳ 15-20 năm.

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

Hạn chế kỹ thuật

  • Số lượng ô khảo sát dừng lại ở mức 60 ô (30 ô/loài) do địa hình núi đá vôi hiểm trở, chia cắt phức tạp.
  • Chưa tiến hành phân tích sinh hóa chuyên sâu về hàm lượng kháng sinh thực vật (phytoncide) bài tiết qua rễ và lá giữa các cặp loài có tương quan dương.
  • Thời gian nghiên cứu tập trung trong 5 tháng (01/2018 - 05/2018), chưa theo dõi biến động động thái rừng qua các mùa khô và mùa mưa lũ liên tiếp.

Hướng nghiên cứu tiếp theo

  • Thiết lập hệ thống ô định vị lâu dài ($1\text{ ha}$) tại VQG Ba Bể để giám sát tốc độ tăng trưởng sinh khối ($D_{1.3}, H_{vn}$) theo thời gian thực.
  • Ứng dụng thiết bị bay không người lái (UAV LiDAR) để quét cấu trúc 3D tán rừng và trích xuất chỉ số diện tích lá (LAI).
  • Mở rộng khảo sát sang các khu bảo tồn thiên nhiên lân cận (Na Hang, Kim Hỷ) nhằm kiểm chứng tính phổ quát của mô hình liên kết loài.

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

  • Sinh viên & Học viên Lâm nghiệp: Tiếp cận phương pháp luận nghiên cứu định lượng cấu trúc rừng, ứng dụng toán xác suất và mô hình ô 6 cây trong phân tích sinh thái thực vật.
  • Kỹ sư Lâm sinh & Cán bộ Kiểm lâm: Tài liệu hướng dẫn thực tế để lập kế hoạch thiết kế vườn ươm, tuyển chọn hạt giống và bố trí cây trồng cho các dự án trồng rừng thay thế, trồng rừng đặc dụng.
  • Ban Quản lý Rừng phòng hộ & VQG: Cơ sở khoa học vững chắc để phê duyệt đề án bảo tồn đa dạng sinh học và phục hồi hệ sinh thái núi đá vôi vùng Đông Bắc.
  • Nhà khoa học & Chuyên gia Môi trường: Dữ liệu chuẩn (benchmark data) về tương tác quần xã cây gỗ nhiệt đới gió mùa phục vụ các mô hình dự báo hấp thụ carbon.

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

1. Tại sao cây Thung và cây Phay không nên trồng thuần loài với mật độ dày?

Cả Thung và Phay đều là cây gỗ lớn ưa sáng, tán rộng ($D_t > 6.5\text{ m}$). Khi trồng thuần loài mật độ cao, cây sẽ cạnh tranh gay gắt về không gian dinh dưỡng và ánh sáng, dẫn đến hiện tượng tỉa cành tự nhiên kém, thân cong queo và dễ bùng phát dịch nấm hại lá do độ ẩm cục bộ dưới tán không được lưu thông.

2. Có thể thay thế cây Ô rô bằng loài cây gỗ khác trong mô hình trồng Thung không?

Ô rô (Streblus ilicifolius) đóng vai trò là cây che bóng tầng dưới và giữ ẩm đất với tần suất xuất hiện lên tới 93.33% ($V = +0.40$). Nếu thay thế, bắt buộc phải chọn các loài cây chịu bóng có cùng ổ sinh thái như Sảng đá hoặc Dú dẻ, tuyệt đối không dùng các loài cây ưa sáng cạnh tranh trực tiếp tầng tán.

3. Bán kính điều tra $R = 4.5\text{ m}$ có áp dụng được cho mọi kiểu rừng nhiệt đới không?

Bán kính $R = 4.0 - 5.0\text{ m}$ được chuẩn hóa riêng cho kiểu rừng thứ sinh phục hồi trung niên trên núi đá vôi. Đối với rừng nguyên sinh giàu có đường kính thân cây lớn ($D_{1.3} > 60\text{ cm}$), bán kính ô điều tra cần mở rộng lên $R = 8.0 - 10.0\text{ m}$ để phản ánh đúng vùng dinh dưỡng.

4. Quy trình xử lý thực bì khi trồng rừng hỗn giao Thung - Phay như thế nào?

Không được đốt dọn thực bì toàn diện vì sẽ làm mất tầng mùn mỏng trên núi đá vôi. Cần áp dụng biện pháp phát dọn cục bộ theo băng hoặc theo hố tròn bán kính $1.0\text{ m}$, giữ lại thảm cây bụi tái sinh tự nhiên để bảo vệ vi sinh vật đất.

5. Thời điểm trồng rừng tốt nhất tại khu vực Bắc Kạn là khi nào?

Thời vụ tối ưu là vụ Xuân (tháng 2 đến tháng 4) khi đất đủ ẩm và bắt đầu có mưa phùn, hoặc đầu vụ Thu (tháng 7 đến tháng 8), tránh các đợt sương muối khắc nghiệt vào tháng 12 - tháng 1.

Kết luận

Đề tài đã xác lập cơ sở khoa học định lượng vững chắc về mối quan hệ cộng sinh giữa hai loài cây gỗ bản địa ưu thế Thung (Tetrameles nudiflora) và Phay (Duabanga sonneratioides) với các loài cây bạn đồng hành tại VQG Ba Bể, Bắc Kạn. Bằng việc kết hợp phương pháp điều tra ô 6 cây cải tiến, trắc đồ lâm học và kiểm định xác suất liên loài, nghiên cứu đã chứng minh tính tất yếu của việc trồng rừng hỗn giao đa tầng tán: ghép nối Thung với Ô rô, Thị đá, Mạy tèo và Phay với Mạy tèo, Trai lí.

Kết quả này không chỉ đóng góp trực tiếp vào kho tàng dữ liệu lâm học vùng núi phía Bắc mà còn cung cấp một khung hướng dẫn kỹ thuật chuẩn xác, khả thi cho các dự án phục hồi rừng phòng hộ đầu nguồn và quản lý tài nguyên rừng bền vững trên toàn quốc.