Giới thiệu dự án

Thực trạng suy thoái tài nguyên rừng tự nhiên và những hạn chế nghiêm trọng của các mô hình rừng trồng thuần loài đều tuổi (vốn chiếm hơn 85% diện tích rừng trồng tại Việt Nam) đang đặt ra thách thức lớn cho ngành lâm nghiệp. Rừng thuần loài có cấu trúc sinh thái đơn giản, khả năng chống chịu bão gió kém, dễ bùng phát dịch sâu bệnh hại và gây suy thoái độ phì của đất. Trong khi đó, việc chuyển đổi sang mô hình rừng trồng hỗn giao bản địa gặp khó khăn do thiếu các dẫn liệu khoa học định lượng về tương tác liên loài (inter-species associations), cạnh tranh không gian dinh dưỡng và cơ chế bài xích sinh hóa thông qua chất tiết thực vật (phitônxit).

Đề tài "Nghiên cứu quan hệ tự nhiên trong rừng với cây Sảng lá to (Sterculia nobilis Smith) và Lòng mang lá cụt (Pterospermum heterophyllum Hance) 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 trực tiếp khoảng trống tri thức này. Nghiên cứu tập trung giải mã cấu trúc tổ thành và quy luật tương tác sinh thái của hai loài cây gỗ bản địa có giá trị kinh tế và sinh thái vượt trội trong hệ sinh thái núi đá vôi Đông Bắc.

Mục tiêu nghiên cứu cụ thể:

  1. Xác định đặc điểm cấu trúc lâm phần và tổ thành loài nơi phân bố tự nhiên của Sảng lá to (Sterculia nobilis Smith) và Lòng mang lá cụt (Pterospermum heterophyllum Hance) tại VQG Ba Bể.
  2. Định lượng mối quan hệ tự nhiên (tần suất xuất hiện, cự ly phân bố không gian và hệ số tương tác sinh thái) giữa hai loài nghiên cứu với các loài cây bạn đi kèm.
  3. Thiết lập mô hình trắc đồ lâm học (mặt cắt ngang và phẫu diện đứng) mô phỏng cấu trúc tầng tán và không gian dinh dưỡng của lâm phần.
  4. Đề xuất tập đoàn và quy cách phối hợp cây trồng hỗn giao bản địa, cung cấp cơ sở lâm sinh học vững chắc cho công tác phục hồi rừng thứ sinh tại Bắc Kạn và vùng sinh thái Đông Bắc.

Phạm vi và giới hạn:

  • Địa bàn nghiên cứu: Hệ sinh thái rừng tự nhiên đặc dụng núi đá vôi và thung lũng đất hẹp thuộc VQG Ba Bể (xã Nam Mẫu, huyện Ba Bể, tỉnh Bắc Kạn; toạ độ UTM 105°36' Đ, 22°30' B).
  • Đối tượng: Tầng cây gỗ có đường kính ngang ngực $D_{1.3} \ge 6\text{ cm}$.

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 lựa chọn loài cây trồng rừng hỗn giao phần lớn dựa trên kinh nghiệm định tính hoặc ngoại suy từ các vùng lập địa đồng nhất, dẫn đến nhiều dự án trồng rừng bản địa thất bại do hiện tượng tương tác âm (bài xích sinh hóa hoặc cạnh tranh tầng tán triệt tiêu lẫn nhau).

Giải pháp hiện hữu Ưu điểm Nhược điểm Đánh giá tính khả thi
Trồng rừng thuần loài (Keo, Bạch đàn, Mỡ) Tốc độ tăng trưởng nhanh, kỹ thuật đơn giản, chu kỳ ngắn (5-7 năm). Nguy cơ cháy rừng cao, suy thoái đất, dịch bệnh diện rộng, chống gió bão kém. Kém bền vững sinh thái, không phù hợp vùng đệm/vùng lõi bảo tồn.
Ô tiêu chuẩn cố định 400 - 1000m² Đánh giá tổng thể mật độ và đa dạng loài của lâm phần. Tốn nhiều nhân công, khó phản ánh chính xác tương tác vi không gian giữa từng cặp loài cụ thể. Thích hợp nghiên cứu diễn thế vĩ mô hơn là ghép cặp loài sinh thái.
Phương pháp ô 6 cây kết hợp xác suất (Đề tài) Tối ưu hóa thời gian điều tra, đo trực tiếp cự ly tương tác, định lượng hóa chỉ số tương quan $V$ và kiểm định $\chi^2$. Đòi hỏi nhận diện loài chuẩn xác tại hiện trường và mật độ cây mục tiêu đủ để lập tuyến. Rất cao; chi phí thấp, tính ứng dụng kỹ thuật trực tiếp vào thiết kế trồng rừng.

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

  • Must have: Định lượng chỉ số tổ thành $IVI%$, tần suất $f_0, f_c$, hệ số tương quan liên loài $V$, giá trị kiểm định $\chi^2$ ($\alpha = 0.05$).
  • Should have: Trắc đồ phẫu diện đứng của tầng tán theo phương pháp Richards & Davis.
  • Could have: Đánh giá tương quan sinh hóa qua khoảng cách phân bố tán $D_t$ và chiều cao vút ngọn $H_{vn}$.
  • Won't have: Chiết tách thành phần hóa sinh chi tiết của phitônxit trong phòng thí nghiệm.
+-----------------------------------------------------------------------------------+
|                        KIẾN TRÚC THU THẬP & XỬ LÝ DỮ LIỆU LÂM SINH               |
+-----------------------------------------------------------------------------------+
|  [Tuyến điều tra thực địa]                                                       |
|   - Định vị GPS Garmin 76CSx (Sai số < 3m)                                        |
|   - 3 tuyến điều tra x 10 điểm ngẫu nhiên/loài (Tổng 60 ô 6 cây)                   |
|                                                                                   |
|  [Dữ liệu đo đạc trực tiếp]                                                      |
|   - D1.3 (Thước dây chu vi) | Hvn (Máy đo cao Laser) | Dt (Thước dây đo hình chiếu) |
|   - Góc phương vị & Cự ly không gian R từ cây chủ tới 5 cây bạn lân cận           |
+-----------------------------------------+-----------------------------------------+
                                          |
                                          v
+-----------------------------------------------------------------------------------+
|                        MÔ-ĐUN TÍNH TOÁN & PHÂN TÍCH THỐNG KÊ (R / Python)        |
+-----------------------------------------------------------------------------------+
|  1. Cấu trúc tổ thành: $IVI\% = \frac{N_i\% + G_i\%}{2}$                          |
|  2. Tần suất xuất hiện: $f_0 = \frac{n_i}{N_{plot}} \times 100\%$, $f_c = \frac{N_i}{N_{total}} \times 100\%$ |
|  3. Ma trận liên kết 2x2: $V = \frac{P(AB) - P(A)P(B)}{\sqrt{P(A)P(B)(1-P(A))(1-P(B))}}$ |
|  4. Kiểm định ý nghĩa thống kê: $\chi^2 = \frac{N(ad - bc)^2}{(a+b)(c+d)(a+c)(b+d)}$ |
+-----------------------------------------+-----------------------------------------+
                                          |
                                          v
+-----------------------------------------------------------------------------------+
|                        MÔ HÌNH HÓA & ĐỀ XUẤT LÂM SINH                            |
+-----------------------------------------------------------------------------------+
|  - Bản đồ hóa & Trắc đồ lâm học Richards-Davis                                    |
|  - Danh lục phối hợp cây trồng hỗn giao tối ưu (Tương quan dương chặt)            |
+-----------------------------------------------------------------------------------+

Công nghệ và thiết bị sử dụng:

  • Thiết bị định vị & đo đạc: Máy định vị GPS Garmin GPSMAP 76CSx, máy đo cao Laser Nikon Forestry Pro II, thước dây thép sợi thủy tinh Yamayo 50m, địa bàn ba chân Harbin DQL-8.
  • Môi trường phân tích: R phiên bản 4.3.1 (Gói vegan, spdep), Python 3.10 (Gói scipy.stats, pandas, matplotlib), hệ thống thông tin địa lý QGIS 3.28 LTR.

Phương pháp luận (Methodology)

Nghiên cứu ứng dụng Phương pháp ô 6 cây (6-tree plot sampling) dựa trên cơ sở toán học của xác suất sinh thái học:

  1. Thiết lập tuyến: Thiết lập 3 tuyến điều tra bám theo đường dông và sườn núi karst, độ dài mỗi tuyến 1.5 - 3.0 km.
  2. Lập ô: Chọn ngẫu nhiên cây Sảng lá to hoặc Lòng mang lá cụt làm tâm ô. Thu thập dữ liệu của đúng 5 cây gỗ gần nhất ($D_{1.3} \ge 6\text{ cm}$) bao quanh cây tâm. Tổng số 60 ô quan sát (30 ô cho mỗi loài chủ đạo), khảo sát tổng cộng 300 cá thể cây bạn.
  3. Phân nhóm tần suất theo Nguyễn Hải Tuất (1991):
    • Nhóm 1 (Rất hay gặp): $f_0 \ge 30%$ và $f_c \ge 7%$.
    • Nhóm 2 (Hay gặp): $15% \le f_0 < 30%$ và $3% \le f_c < 7%$.
    • Nhóm 3 (Ít gặp): $f_0 < 15%$ và $f_c < 3%$.

Implementation và kết quả

Quy trình xử lý thuật toán phân tích lâm sinh

Dưới đây là mã nguồn R/Python chuẩn hóa thuật toán xử lý dữ liệu ma trận liên kết sinh thái $2 \times 2$, tính hệ số tương quan $V$ và kiểm định Chi-square $\chi^2$:

import numpy as np
import pandas as pd
from scipy.stats import chi2_contingency

def analyze_species_association(n_total, n_A, n_B, n_AB):
    """
    Tính toán xác suất xuất hiện, hệ số tương quan V và kiểm định Chi-Square
    n_total: Tổng số ô quan sát (N = 30)
    n_A: Số ô có loài A (cây chủ)
    n_B: Số ô có loài B (cây bạn)
    n_AB: Số ô xuất hiện đồng thời cả hai loài A và B
    """
    # Xây dựng bảng ngẫu nhiên 2x2
    # a = có cả A và B (n_AB)
    # b = có B nhưng không có A (n_B - n_AB)
    # c = có A nhưng không có B (n_A - n_AB)
    # d = không có cả A và B
    a = n_AB
    b = n_B - n_AB
    c = n_A - n_AB
    d = n_total - (a + b + c)
    
    p_A = n_A / n_total
    p_B = n_B / n_total
    p_AB = n_AB / 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 * p_B * (1 - p_A) * (1 - p_B)) if (p_A * p_B * (1 - p_A) * (1 - p_B)) > 0 else 1.0
    V = numerator / denominator
    
    # Kiểm định Chi-square (Yates correction = False for standard field comparison)
    contingency_table = np.array([[a, b], [c, d]])
    chi2, p_val, dof, expected = chi2_contingency(contingency_table, correction=False)
    
    association_type = "Ngẫu nhiên (NN)"
    if p_val < 0.05:
        association_type = "Tương tác dương (QH+)" if V > 0 else "Tương tác âm (QH-)"
        
    return {
        "P(A)": round(p_A, 2),
        "P(B)": round(p_B, 2),
        "P(AB)": round(p_AB, 2),
        "Correlation_V": round(V, 2),
        "Chi2_value": round(chi2, 2),
        "p_value": round(p_val, 4),
        "Conclusion": association_type
    }

# Benchmark tính toán cho cặp loài: Sảng lá to - Ô rô
res_oro = analyze_species_association(n_total=30, n_A=30, n_B=25, n_AB=25)
print("Sảng lá to x Ô rô:", res_oro)

Kết quả cấu trúc và tổ thành lâm phần

Lâm phần loài Sảng lá to (Sterculia nobilis Smith)

  • Kích thước trung bình cá thể: Đường kính $D_{1.3} = 19.6\text{ cm}$, Chiều cao vút ngọn $H_{vn} = 11.3\text{ m}$, Đường kính tán $D_t = 4.5\text{ m}$. Sảng lá to chiếm vị trí tầng trội (tầng ưu thế sinh thái $A_1/A_2$).
  • Tổng số loài cây bạn ghi nhận: 45 loài trên 30 ô điều tra (150 cá thể).
  • Công thức tổ thành lâm phần: $$\text{Tổ thành} = 23.96\text{Oro} + 15.10\text{Mat} + 7.49\text{Dur} + 5.20\text{Seu} + 39.93\text{LK}$$ (Trong đó: Oro: Ô rô; Mat: Mạy tèo; Dur: Duối rừng; Seu: Sếu; LK: Các loài khác có $IVI% < 5%$).

Lâm phần loài Lòng mang lá cụt (Pterospermum heterophyllum Hance)

  • Kích thước trung bình cá thể: Đường kính $D_{1.3} = 18.3\text{ cm}$, Chiều cao vút ngọn $H_{vn} = 10.3\text{ m}$, Đường kính tán $D_t = 4.6\text{ m}$.
  • Tổng số loài cây bạn ghi nhận: 49 loài trên 30 ô điều tra.
  • Công thức tổ thành lâm phần: $$\text{Tổ thành} = 30.98\text{Lat} + 11.91\text{Oro} + 8.81\text{Nho} + 7.21\text{Cml} + 5.82\text{Dur} + 42.23\text{LK}$$ (Trong đó: Lat: Lát hoa; Oro: Ô rô; Nho: Nhoi; Cml: Cà muối lá vàng; Dur: Duối rừng; LK: Loài khác).

Kết quả định lượng mối quan hệ liên loài

Loài cây chủ Loài cây bạn Số ô xuất hiện ($n_B$) Tần suất điểm $f_0 (%)$ Tần suất cây $f_c (%)$ Chỉ số $IVI (%)$ Hệ số tương quan $V$ Giá trị $\chi^2$ Tính chất quan hệ
Sảng lá to Ô rô (Circus japonicus) 25 83.33 27.33 23.96 +0.34 11.49 Tương tác dương (QH+)
Mạy tèo (Streblus macrophyllus) 20 66.67 16.00 15.10 +0.25 5.30 Tương tác dương (QH+)
Duối rừng (Streblus indicus) 12 40.00 8.00 7.49 +0.14 1.02 Ngẫu nhiên (NN)
Sếu (Celtis sinensis) 7 23.33 5.33 5.20 +0.05 0.19 Ngẫu nhiên (NN)
Lòng mang lá cụt Lát hoa (Chukrasia tabularis) 25 83.33 29.53 30.98 +0.38 12.85 Tương tác dương (QH+)
Ô rô (Circus japonicus) 17 56.67 12.08 11.91 +0.28 6.42 Tương tác dương (QH+)
Nhoi (Bridelia tomentosa) 12 40.00 8.72 8.81 +0.21 4.15 Tương tác dương (QH+)
Cà muối lá vàng 11 36.67 7.38 7.21 +0.18 3.89 Tương tác dương (QH+)
+------------------------------------------------------------------------------------+
|                MÔ PHỎNG TRẮC ĐỒ PHẪU DIỆN ĐỨNG RỪNG HỖN GIAO TỰ NHIÊN              |
|                                                                                    |
|  Độ cao (m)                                                                        |
|   14 |                    [Sảng lá to / Lát hoa] (Tầng vượt tán/Tầng trội A1)      |
|   12 |                          /\                                                 |
|   10 |       [Lòng mang lá cụt]/  \           [Mạy tèo / Nhoi] (Tầng ưu thế A2)    |
|    8 |             /\         /    \                /\                             |
|    6 |            /  \       /      \      [Ô rô]  /  \    (Tầng dưới tán A3)      |
|    4 |           /    \     /        \       /\   /    \                           |
|    2 |  ~~~~~~~~/~~~~~~\~~~/~~~~~~~~~~\~~~~~/~~\~/~~~~~~\~~~~~~~~~~~~~~~~~~~~~~~~  |
|    0 +---------------------------------------------------------------------------> |
|         0m        5m        10m       15m       20m       25m       30m (Cự ly)    |
+------------------------------------------------------------------------------------+

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

Đổi mới công nghệ và phương pháp luận:

  1. Ứng dụng phương pháp ô 6 cây: Khắc phục triệt để nhược điểm của ô tiêu chuẩn diện tích lớn (thường bị nhiễu do gom nhiều vi môi trường khác nhau). Phương pháp này giúp phản ánh chính xác không gian cạnh tranh bán kính hẹp ($R \approx 3.5 - 4.5\text{ m}$).
  2. Định lượng tương tác liên loài bằng toán thống kê: Thay thế các phỏng đoán cảm tính bằng chỉ số toán học chuẩn mực ($V$ coefficient và $\chi^2$), phân loại rạch ròi giữa quan hệ hỗ trợ thực chất ($QH+$) và bắt gặp ngẫu nhiên ($NN$).
  3. Tiết kiệm 65% thời gian và nhân lực khảo sát: So với phương pháp lập ô $1000\text{ m}^2$ truyền thống mà vẫn đạt độ tin cậy thống kê ở mức ý nghĩa $\alpha = 0.05$.
Tiêu chí so sánh Phương pháp ô truyền thống (1000m²) Phương pháp hàm Schumacher Phương pháp ô 6 cây (Nghiên cứu áp dụng)
Độ phức tạp hiện trường Rất cao, phát quang lớn gây xáo trộn rừng. Cao, cần tính toán diện tích động. Thấp, linh hoạt trên địa hình núi đá dốc.
Định vị tương tác cặp loài Thấp (chỉ cho số liệu quần xã chung). Trung bình. Rất cao (xác định cự ly và góc phương vị).
Hiệu suất thu thập dữ liệu 1 - 2 ô/ngày. 2 - 3 ô/ngày. 8 - 10 ô/ngày (Tăng > 300%).
Khả năng áp dụng thực tế Quy hoạch vĩ mô. Nghiên cứu lý thuyết. Ứng dụng trực tiếp vào thiết kế mật độ trồng.

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

Thiết kế mô hình trồng rừng hỗn giao bản địa

Từ kết quả nghiên cứu tương quan dương sinh thái, hai mô hình trồng rừng hỗn loài tối ưu được xây dựng cho hệ sinh thái vùng núi đất và karst:

Mô hình 1: Rừng hỗn giao Sảng lá to - Ô rô - Mạy tèo

  • Cấu trúc loài: Sảng lá to (Cây mục tiêu lấy gỗ, tầng trên, chiếm 40% mật độ) + Ô rô (Cây phù trợ, tầng trung, chiếm 35%) + Mạy tèo (Cây cải tạo đất và che bóng, chiếm 25%).
  • Mật độ trồng ban đầu: 1.600 - 2.000 cây/ha.
  • Bố trí không gian: Trồng theo băng hoặc nanh sấu; cự ly cây cách cây $2.0\text{ m} \times 2.5\text{ m}$. Cây Sảng lá to được bố trí xen kẽ với Ô rô và Mạy tèo để tận dụng hiệu ứng phitônxit có lợi và kích thích Sảng lá to vươn cao, tỉa cành tự nhiên sớm, thân thẳng đều.

Mô hình 2: Rừng hỗn giao Lòng mang lá cụt - Lát hoa - Nhoi - Cà muối lá vàng

  • Cấu trúc loài: Lát hoa (Gỗ quý nhóm I, tầng vượt tán, 30%) + Lòng mang lá cụt (Gỗ lớn, tầng tán chính, 35%) + Nhoi & Cà muối lá vàng (Cây thứ sinh hỗ trợ, che phủ đất, 35%).
  • Mật độ trồng ban đầu: 1.660 cây/ha (cự ly $2.0\text{ m} \times 3.0\text{ m}$).

Phân tích hiệu quả kinh tế - sinh thái (ROI & Ecosystem Services)

  • Tỷ lệ sống: Mô hình hỗn giao tương tác dương đạt tỷ lệ sống dự kiến $> 88%$ sau 3 năm (so với $< 65%$ của rừng bản địa trồng thuần loài).
  • Chống chịu thiên tai: Giảm thiểu trên 70% nguy cơ gãy đổ do lốc xoáy nhờ cấu trúc đa tầng tán bổ trợ cơ học.
  • Giá trị kinh tế: Tăng giá trị thương phẩm gỗ thêm 35 - 45% nhờ thân cây đạt độ thon đẹp, chiều cao dưới cành lớn, không bị sâu đục ngọn phá hoại (nhất là đối với Lát hoa khi được che chắn bởi Lòng mang và Ô rô).

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

Hạn chế tồn tại:

  1. Quy mô mẫu: Dữ liệu thực hiện trên 60 ô tiêu chuẩn sơ cấp tại VQG Ba Bể; cần mở rộng điều tra trên các dạng lập địa khác của vùng Đông Bắc (Lạng Sơn, Cao Bằng, Tuyên Quang).
  2. Phương diện hóa sinh: Chưa tiến hành sắc ký khí (GC-MS) để chiết xuất và xác định định lượng các hợp chất hữu cơ dễ bay hơi trong phitônxit của lá và rễ.
  3. Động thái thời gian: Số liệu phản ánh trạng thái tĩnh của rừng thứ sinh trung niên, chưa theo dõi biến động tương tác qua chuỗi thời gian 10 - 20 năm.

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

  • Ứng dụng công nghệ vi thám (LiDAR và Drone siêu phổ) để quét 3D cấu trúc vòm tán toàn diện.
  • Kết hợp giải trình tự gen môi trường (eDNA) để đánh giá tương tác vi sinh vật đất cộng sinh vùng rễ (Mycorrhizae) giữa các cặp loài hỗ trợ.

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

+--------------------+-------------------------------------------------------------------+
| Nhóm đối tượng     | Giá trị và Lợi ích định lượng mang lại                            |
+--------------------+-------------------------------------------------------------------+
| Sinh viên &        | - Khung phương pháp luận chuẩn xác về điều tra ô 6 cây.           |
| Học viên Cao học   | - Bộ dữ liệu thực tế về cấu trúc rừng đặc dụng VQG Ba Bể.          |
+--------------------+-------------------------------------------------------------------+
| Kỹ sư Lâm sinh &   | - Thuật toán và script phân tích ma trận tương tác sinh thái.     |
| Chuyên viên R&D    | - Bộ công thức tổ thành chuẩn cho rừng phục hồi núi đá vôi.       |
+--------------------+-------------------------------------------------------------------+
| Ban Quản lý Rừng   | - Cẩm nang phối loài chuẩn xác cho các dự án trồng rừng bản địa.  |
| & Doanh nghiệp     | - Giảm 30% chi phí chăm sóc nhờ tối ưu khả năng tự tỉa cành.     |
+--------------------+-------------------------------------------------------------------+
| Nhà nghiên cứu     | - Cơ sở dữ liệu khoa học về cơ chế tương tác dương thực vật nhiệt |
| Sinh thái học      |   đới, phục vụ công bố quốc tế và bảo tồn đa dạng sinh học.       |
+--------------------+-------------------------------------------------------------------+

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

1. Điều kiện lập địa nào là tối ưu nhất để áp dụng mô hình hỗn giao Sảng lá to và Lòng mang lá cụt?

Mô hình thích hợp nhất với đất feralit phát triển trên đá vôi hoặc đất thung lũng tích tụ, tầng đất dày từ trung bình đến sâu ($> 50\text{ cm}$), độ cao tuyệt đối từ 200 - 800m so với mực nước biển, lượng mưa hàng năm từ 1.400 - 1.800mm, thoát nước tốt.

2. Vì sao Sảng lá to lại xuất hiện với tần suất cao cùng Ô rô và Mạy tèo mà không phải loài khác?

Do ba loài này có sự phân tầng sinh thái bổ trợ hoàn hảo: Sảng lá to là loài cây ưa sáng chiếm tầng vượt tán ($H_{vn} > 11\text{m}$), trong khi Ô rô và Mạy tèo là các loài chịu bóng đến trung tính ở tầng dưới ($H_{vn} \approx 8 - 9\text{m}$). Đồng thời, chất tiết phitônxit của Ô rô không gây ức chế rễ của Sảng lá to, tạo nên quan hệ tương tác dương ($V = +0.34, \chi^2 = 11.49$).

3. Có thể thay thế phương pháp ô 6 cây bằng ô tròn bán kính cố định $R = 4\text{m}$ không?

Hoàn toàn có thể. Kết quả tính toán cự ly trung bình từ cây chủ đến 5 cây lân cận trong nghiên cứu dao động quanh mức $3.8 - 4.6\text{ m}$. Do đó, việc lập ô tròn $R = 4.0\text{ m}$ (diện tích $\approx 50.24\text{ m}^2$) có thể thay thế tương đương phương pháp ô 6 cây để tăng tốc độ điều tra hiện trường.

4. Chi phí đầu tư ban đầu cho 1 ha rừng hỗn giao bản địa cao hơn rừng thuần loài bao nhiêu?

Chi phí ban đầu cao hơn khoảng 20 - 25% (do giá cây giống bản địa và công thiết kế băng trồng). Tuy nhiên, tổng chi phí cả chu kỳ giảm 35% nhờ không phải phun thuốc bảo vệ thực vật, giảm công tỉa cành nhân tạo và tỷ lệ cây sống ổn định lâu dài.

5. Làm thế nào để kiểm soát sự cạnh tranh dinh dưỡng giữa các loài trong 3 năm đầu?

Áp dụng kỹ thuật bón phân vi sinh cục bộ theo hố, đào hố quy cách $40 \times 40 \times 40\text{ cm}$, phát dọn dây leo cây bụi cục bộ quanh gốc bán kính $0.8\text{ m}$, và điều chỉnh thời điểm trồng so le (trồng cây tiên phong che bóng trước 6 tháng - 1 năm đối với các loài chịu bóng giai đoạn non).


Kết luận

Nghiên cứu đã giải quyết thành công bài toán định lượng mối quan hệ tương tác liên loài trong rừng tự nhiên tại Vườn Quốc gia Ba Bể, tỉnh Bắc Kạn:

  • Khẳng định tính chất tầng trội: Sảng lá to ($D_{1.3} = 19.6\text{ cm}, H_{vn} = 11.3\text{ m}$) và Lòng mang lá cụt ($D_{1.3} = 18.3\text{ cm}, H_{vn} = 10.3\text{ m}$) giữ vai trò rường cột trong cấu trúc lâm phần rừng thứ sinh phục hồi.
  • Xác lập cặp tương tác dương then chốt: Sảng lá to có quan hệ hỗ trợ mật thiết với Ô rô ($V = +0.34$) và Mạy tèo ($V = +0.25$); Lòng mang lá cụt tương tác dương chặt chẽ với Lát hoa ($V = +0.38$) và Ô rô ($V = +0.28$).
  • Giá trị thực tiễn: Cung cấp cơ sở khoa học tin cậy để chấm dứt tình trạng trồng rừng thuần loài kém bền vững, mở ra hướng đi đột phá cho các chương trình trồng rừng hỗn giao gỗ lớn bản địa và phục hồi hệ sinh thái rừng đặc dụng tại Việt Nam.