Giới thiệu dự án

Hệ sinh thái rừng mưa nhiệt đới trên núi đá vôi tại Việt Nam đóng vai trò thiết yếu trong bảo tồn đa dạng sinh học, điều tiết thủy văn và lưu trữ carbon. Tuy nhiên, theo số liệu từ Tổng cục Lâm nghiệp, các hệ sinh thái karst nhạy cảm đang chịu áp lực suy thoái nghiêm trọng do biến đổi khí hậu và tác động nhân sinh, làm suy giảm quần thể các loài cây gỗ bản địa quý hiếm.

Gội trắng (Aphanamixis grandiflora Blume, thuộc họ Xoan - Meliaceae) là loài cây gỗ lớn cao 35–40 m, đường kính ngang ngực ($D_{1.3}$) đạt tới 90 cm, mang giá trị đa dụng vượt trội: cung cấp gỗ thương phẩm chất lượng cao, vỏ và lá chứa hợp chất sinh học trị bệnh ngoài da, hạt giàu hàm lượng lipid phục vụ sản xuất tinh dầu và xà phòng sinh học. Dù có phân bố tự nhiên rộng nhưng loài này đang bị khai thác quá mức và chưa có quy trình lâm sinh chuyên biệt phục vụ tái thả, gây trồng hay phục hồi sinh thái.

Mục tiêu dự án

  1. Định lượng cấu trúc lâm phần: Xác định công thức tổ thành ($CTTT$), mật độ cá thể ($N/\text{ha}$), cấu trúc tầng thứ thẳng đứng và độ tàn che tán rừng tại 2 đai cao phân bố (400 m và 500 m).
  2. Đánh giá quy luật sinh trưởng: Phân tích cấu trúc phân bố số cây theo cấp đường kính ($N/D_{1.3}$) và cấp chiều cao ($N/H_{vn}$) của loài Gội trắng và toàn lâm phần.
  3. Giải mã mạng lưới tương tác sinh thái: Xác định nhóm loài cây đi kèm chủ yếu thông qua phương pháp mạng hình 6 cây lân cận.
  4. Đánh giá tiềm năng tái sinh tự nhiên: Định lượng mật độ cây tái sinh ($N_{\text{ts}}/\text{ha}$), nguồn gốc xuất hiện (hạt vs chồi), tỷ lệ cây có triển vọng và mối quan hệ với tầng cây bụi thảm tươi.

Phương pháp tiếp cận và Phạm vi

  • Phương pháp: Ứng dụng điều tra lâm học thực địa kết hợp mô hình hóa thống kê xác suất (phân bố thực nghiệm, moment bậc cao $S_k, E_x$) trên nền tảng Microsoft Excel và chuẩn hóa lâm sinh học nhiệt đới.
  • Phạm vi & Giới hạn: Nghiên cứu thực hiện tại phân khu bảo vệ nghiêm ngặt của Vườn Quốc gia Cúc Phương (Ninh Bình, 20°14'–20°24'N, 105°29'–105°44'E) trong khung thời gian từ 15/01/2018 đến 06/05/2018. Mẫu đất giới hạn ở phương pháp mục trắc do tầng đá vôi hiểm trở, rạn nứt karst sâu.

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

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

Các phương pháp đánh giá tài nguyên rừng truyền thống thường đối mặt với sự thiếu hụt dữ liệu định lượng vi mô về tương tác cá thể và khả năng tái sinh theo đai độ cao.

Tiêu chí Khảo sát lâm sinh truyền thống Mô hình hóa lâm học định lượng (Đề tài)
Độ chính xác tầng tán Mục trắc trực quan, sai số 20–30% Đo cao Blume-Leiss + Trắc đồ đứng Richards (1964)
Tương tác loài đi kèm Nhận xét cảm tính, thiếu tần số Ma trận lân cận 6 cây + Phân hạng Triệu Văn Hùng
Quy luật phân bố Ước tính bình quân số học Moment thống kê: Phương sai ($S^2$), Lệch ($S_k$), Nhọn ($E_x$)
Định hướng ứng dụng Can thiệp chung chung Thiết kế mô hình hỗn giao rừng trồng đặc dụng

Yêu cầu hệ thống phương pháp (Phân loại MoSCoW)

  • Must have: Thu thập $D_{1.3}$, $H_{vn}$, $H_{dc}$, $D_t$ trên toàn bộ cây có $D_{1.3} \ge 6\text{ cm}$; thiết lập 6 ô tiêu chuẩn ($S = 500\text{ m}^2$) và 30 ô dạng bản ($S = 9\text{ m}^2$).
  • Should have: Phân tích tương quan độ tàn che bằng trắc đồ bằng; phân loại phẩm chất cây đứng thành 3 cấp (Tốt, Trung bình, Xấu).
  • Could have: Mở rộng khảo sát động thái vật hậu ra các vùng đệm giáp ranh Hòa Bình, Thanh Hóa.
  • Won't have: Giải trình tự gen phân tử quần thể trong phạm vi khóa luận cử nhân.

Thiết kế công cụ và quy trình thu thập dữ liệu

  • Công cụ đo đạc thực địa: Thước vanh thép dải đo 0–200 cm (sai số $\pm 0.1\text{ cm}$), thước đo cao Blume-Leiss BL6 (độ chính xác $\pm 0.5\text{ m}$), địa bàn khảo sát Suunto Tandem 360PC, GPS cầm tay Garmin GPSMAP 64s.
  • Công cụ xử lý thống kê: Bộ công cụ phân tích dữ liệu lâm nghiệp tích hợp trên Microsoft Excel theo quy chuẩn của GS. Nguyễn Hải Tuất và TS. Ngô Kim Khôi.
  • Mô hình kiến trúc điều tra:
    • Ô tiêu chuẩn chính (ÔTC): Kích thước $20\text{ m} \times 25\text{ m} = 500\text{ m}^2$.
    • Ô dạng bản (ÔDB): 5 ô/ÔTC (4 ô góc, 1 ô tâm), kích thước $3\text{ m} \times 3\text{ m} = 9\text{ m}^2$.
    • Cụm điều tra tương tác: 1 cây Gội trắng trung tâm kết nối 6 cây lân cận gần nhất theo cấu trúc lục lăng sinh thái.
               [Góc 1: ÔDB 9m²]                    [Góc 2: ÔDB 9m²]
               [Góc 3: ÔDB 9m²]                    [Góc 4: ÔDB 9m²]

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

Quy trình thực hiện gồm 4 giai đoạn logic:

  1. Khảo sát tiền trạm: Đối chiếu bản đồ thảm thực vật VQG Cúc Phương, khoanh vùng tọa độ 2 đai cao (400 m và 500 m).
  2. Ngoại nghiệp: Thu thập chỉ số sinh trưởng tầng cây cao, cây tái sinh, cây bụi thảm tươi theo mẫu biểu lâm học chuẩn hóa (Biểu 2.1 đến 2.4).
  3. Nội nghiệp: Phân loại định danh mẫu thực vật tại Bộ môn Thực vật rừng - Đại học Lâm nghiệp; tính toán hệ số tổ thành, lập bảng phân bố $N/D_{1.3}$.
  4. Kiểm chuẩn & Đánh giá: Kiểm định độ tin cậy của các moment thống kê và xác lập giải pháp kỹ thuật lâm sinh.

Implementation và kết quả

Quy chuẩn tính toán và thuật toán lâm sinh

Dữ liệu đo đếm được xử lý qua hệ thống công thức giải tích lâm học:

# Thuật toán tính toán hệ số tổ thành và đặc trưng mẫu lâm học
import math

def calculate_stand_composition(species_counts, total_trees):
    """
    Tính hệ số tổ thành ki cho từng loài trong ô tiêu chuẩn
    ki = (ni / N) * 10
    """
    composition = {}
    mean_trees_per_species = total_trees / len(species_counts)
    
    for sp, count in species_counts.items():
        if count >= mean_trees_per_species:
            ki = round((count / total_trees) * 10, 2)
            composition[sp] = ki
    return composition

def statistical_moments(data_points):
    """
    Tính toán trung bình mẫu, phương sai, hệ số biến động, độ lệch Sk, độ nhọn Ex
    """
    n = len(data_points)
    mean = sum(data_points) / n
    variance = sum((x - mean) ** 2 for x in data_points) / (n - 1)
    std_dev = math.sqrt(variance)
    cv_percent = (std_dev / mean) * 100
    
    # Moment bậc 3 và 4
    m3 = sum((x - mean) ** 3 for x in data_points) / n
    m4 = sum((x - mean) ** 4 for x in data_points) / n
    skewness = m3 / (std_dev ** 3)
    kurtosis = (m4 / (std_dev ** 4)) - 3
    
    return {
        "Mean": mean, "Variance": variance, "StdDev": std_dev,
        "CV%": cv_percent, "Skewness": skewness, "Kurtosis": kurtosis
    }

Các công thức toán học chuyên ngành:

  1. Hệ số tổ thành loài: $$k_i = \frac{n_i}{N} \times 10$$ Trong đó: $n_i$ là số cá thể loài $i$, $N$ là tổng số cá thể quan sát.

  2. Hệ số biến động đường kính ($S%$): $$S = \sqrt{\frac{\sum (X_i - \bar{X})^2 n_i}{N - 1}}, \quad S% = \frac{S}{\bar{X}} \times 100%$$

  3. Chỉ số tương tác loài đi kèm (Phân hạng Triệu Văn Hùng): $$P_0 = \frac{m}{M} \times 100%, \quad P_c = \frac{n_i}{N_{\text{total}}} \times 100%$$

    • Nhóm I (Rất hay gặp): $P_0 > 30%$ và $P_c > 7%$
    • Nhóm II (Hay gặp): $15% \le P_0 \le 30%$ và $3% \le P_c \le 7%$
    • Nhóm III (Ít gặp): $P_0 < 15%$ và $P_c < 3%$

Dữ liệu thực nghiệm và Kiểm định thống kê

1. Cấu trúc tổ thành tầng cây cao

Tổng số loài tham gia vào công thức tổ thành tại đai 400 m là 9 loài, đai 500 m là 7 loài. Vàng anh (Saraca dives) và Sấu (Dracontomelon duperreanum) đóng vai trò ưu thế sinh thái.

Đai cao ÔTC Mật độ ÔTC ($N$) Số loài ($A$) $N_{\text{Gội trắng}}$ Công thức tổ thành tầng cây cao ($CTTT$)
400 m ÔTC 1 19 10 3 $2.10\text{VA} + 1.57\text{GT} + 1.05\text{NR} + 1.05\text{PM} + 1.05\text{N} + 2.08\text{LK}$
ÔTC 2 21 7 3 $3.81\text{VA} + 1.42\text{GT} + 1.42\text{S} + 1.42\text{SA} + 0.95\text{M} + 0.94\text{LK}$
ÔTC 3 20 10 3 $2.00\text{VA} + 1.50\text{GT} + 1.50\text{NR} + 1.50\text{S} + 1.00\text{Đ} + 2.50\text{LK}$
500 m ÔTC 4 24 7 3 $3.30\text{PM} + 1.60\text{VA} + 1.60\text{NV} + 1.25\text{GT} + 1.25\text{THV} + 0.82\text{LK}$
ÔTC 5 20 11 2 $2.50\text{VA} + 2.00\text{NV} + 1.00\text{GT} + 1.00\text{PM} + 3.50\text{LK}$
ÔTC 6 17 9 2 $2.35\text{VA} + 1.76\text{Đ} + 1.17\text{S} + 1.17\text{NV} + 1.17\text{GT} + 2.32\text{LK}$

(Ký hiệu: VA: Vàng anh, GT: Gội trắng, NR: Nhãn rừng, PM: Phân mã, N: Ngát, S: Sấu, SA: Sâng, M: Mang, Đ: Đăng, NV: Nhò vàng, THV: Trà hoa vàng, LK: Loài khác)

2. Đặc trưng sinh trưởng đường kính ($D_{1.3}$) của lâm phần

Đai cao ÔTC $\bar{D}_{1.3}\text{ (cm)}$ $S^2$ $S$ $S%\text{ (%)}$ Độ lệch $S_k$ Độ nhọn $E_x$ Phạm vi $R\text{ (cm)}$
400 m ÔTC 1 25.00 66.58 8.16 36.28 0.36 -1.13 31.50
ÔTC 2 42.00 301.37 17.36 65.39 1.23 0.10 67.70
ÔTC 3 45.33 281.23 16.77 76.06 2.35 5.00 81.70
500 m ÔTC 4 58.50 631.01 25.12 101.50 5.77 21.04 105.00
ÔTC 5 21.00 74.30 8.62 51.50 -2.41 0.52 30.50
ÔTC 6 35.00 282.24 16.80 57.07 -3.66 2.68 55.90
Tần số N
     6-12  13-18 19-24 25-30 >30

3. Cấu trúc tầng thứ và Độ tàn che

  • Tầng vượt tán ($H > 15\text{ m}$): Gồm các cá thể ưu thế như Nhãn rừng, Vàng anh, Chò nhai và Gội trắng thành thục.
  • Tầng tán chính ($H = 8 - 12\text{ m}$): Vàng anh, Gội trắng trung gian, Phân mã.
  • Tầng dưới tán ($H = 3 - 8\text{ m}$): Trà hoa vàng, Sâng, Đại khải.
  • Độ tàn che: Trung bình đạt $0.56$ (biến thiên từ $0.46$ ở ÔTC 4 đến $0.63$ ở ÔTC 3). Lỗ thủng tán tại ÔTC 4 ($0.46$) kích thích mật độ cây bụi thảm tươi phát triển mạnh, tạo áp lực cạnh tranh gay gắt với cây tái sinh non.

4. Phân tích loài đi kèm (Neighboring Species Network)

Từ khảo sát ô 6 cây xung quanh Gội trắng, xác định được 16 loài cây cộng sinh:

  • Nhóm I (Cộng sinh mật thiết): Vàng anh ($P_0 = 100%$, $P_c = 30.55%$), Nhãn rừng ($P_0 = 50%$, $P_c = 8.34%$), Đăng ($P_0 = 50%$, $P_c = 8.34%$).
  • Nhóm II (Cộng sinh trung bình): Ngát, Phân mã, Nhò vàng ($P_0 = 33.34%$, $P_c = 5.56%$).
  • Nhóm III (Cộng sinh ngẫu nhiên): 10 loài gồm Sấu, Chò nhai, Sâng, Thừng mực, Ba gạc, Trà hoa vàng ($P_c = 2.78%$).

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

  1. Ứng dụng cấu trúc mạng 6 cây lục lăng: Thay vì chỉ khảo sát mật độ bình quân diện tích, nghiên cứu xác lập ma trận tương tác 6 cây lân cận cho loài Aphanamixis grandiflora, chứng minh tính tương thích sinh thái cao giữa Gội trắng và Vàng anh (Saraca dives) với tần số xuất hiện $P_0 = 100%$.
  2. So sánh với các nghiên cứu tiền nhiệm:
Chỉ số lâm sinh Nguyễn Đắc Triển et al. (2012) - VQG Xuân Sơn Lê Minh Thuấn (2009) - Tây Nguyên Đề tài nghiên cứu (Thu Hà, 2018) - VQG Cúc Phương
Đối tượng Tái sinh rừng thường xanh Vối thuốc răng cưa Gội trắng (A. grandiflora)
Mật độ cá thể loài Rải rác, $10 - 25\text{ cây/ha}$ $30 - 45\text{ cây/ha}$ 40 - 60 cây/ha (ưu thế ở đai 400m)
Đường kính TB ($\bar{D}_{1.3}$) $18.5\text{ cm}$ $24.2\text{ cm}$ 21.30 - 30.00 cm
Độ tàn che tối ưu $0.6 - 0.7$ $0.7 - 0.8$ 0.50 - 0.63
Động thái phân bố Không tính toán $S_k, E_x$ Phân bố giảm lũy thừa $S_k > 0$ ở 400m; $S_k < 0$ ở 500m
  1. Cơ sở khoa học cho phục hồi rừng karst: Cung cấp bộ thông số chuẩn về đặc trưng sinh thái và lập địa, khẳng định đai cao 400 m là vùng phân bố tối ưu của Gội trắng tại khối núi đá vôi Bắc Bộ.

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

Kịch bản ứng dụng

  • Phục hồi rừng đầu nguồn vùng Karst: Bố trí Gội trắng làm cây tầng cao ($H = 15 - 20\text{ m}$) xen kẽ Vàng anh và Nhãn rừng theo tỷ lệ mật độ $2:3:1$.
  • Lâm nghiệp đô thị & Cảnh quan sinh thái: Tận dụng đặc tính xanh quanh năm, tán lá dày, hoa chùm thơm và khả năng hấp thụ bụi của Gội trắng để phát triển hệ thống công viên cây xanh.
  • Vùng nguyên liệu dược liệu - tinh dầu: Khai thác hạt có hàm lượng dầu béo cao phục vụ công nghiệp hóa mỹ phẩm và y dược học cổ truyền.

Lộ trình kỹ thuật lâm sinh (5 giai đoạn)

Năm 2: Trồng dặm mô hình hỗn giao (Gội trắng + Vàng anh + Nhãn rừng) mật độ 500 cây/ha
Năm 3-4: Phát luỗng dây leo, điều tiết độ tàn che tầng cây che bóng đạt 0.5 - 0.6
Năm 5+: Tỉa thưa định hình tầng tán, nuôi dưỡng cây tái sinh triển vọng

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

  • Hạn chế kỹ thuật: Điều kiện địa hình núi đá vôi karst tai mèo hiểm trở gây khó khăn trong việc đào phẫu diện đất chuyên sâu; chưa phân tích định lượng vi khí hậu dưới tán theo từng tháng trong năm.
  • Hướng phát triển tiếp theo:
    1. Xây dựng mô hình tương quan sinh khối allometric ($AGB = f(D_{1.3}, H_{vn}, \rho)$) riêng cho chi Aphanamixis.
    2. Ứng dụng công nghệ GIS viễn thám độ phân giải cao (Sentinel-2, LiDAR) để lập bản đồ phân vùng sinh thái tiềm năng cho loài Gội trắng trên toàn bộ dãy núi đá vôi Hòa Bình - Ninh Bình - Thanh Hóa.
    3. Phân tích thành phần hóa thực vật (phytochemical screening) từ dịch chiết lá và vỏ cây phục vụ sàng lọc dược chất kháng khuẩn.

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

  • Sinh viên & Học viên Lâm nghiệp: Cung cấp bộ dữ liệu thực nghiệm chuẩn mực về phương pháp điều tra ô tiêu chuẩn, phân tích moment thống kê lâm học ($S_k, E_x$) và xử lý cấu trúc tổ thành.
  • Kỹ sư Lâm sinh & Quản lý Dự án: Nắm bắt chính xác công thức hỗn giao tối ưu ($2.1\text{VA} + 1.57\text{GT}$) để thiết kế hồ sơ kỹ thuật trồng rừng phòng hộ và làm giàu rừng tự nhiên.
  • Ban Quản lý Vườn Quốc gia & Khu Bảo tồn: Căn cứ khoa học phục vụ khoanh nuôi xúc tiến tái sinh, thiết lập khu bảo tồn nguồn gen bản địa tại VQG Cúc Phương.
  • Nhà khoa học Sinh thái: Cơ sở dữ liệu nền về biến thiên cấu trúc quần xã thực vật nhiệt đới trên nền địa hình núi đá vôi có tính chất karst xâm thực đặc thù.

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

1. Điều kiện lập địa nào thích hợp nhất để gây trồng Gội trắng?

Cây Gội trắng sinh trưởng tối ưu trên đất feralit phát triển trên đá vôi hoặc đất renzin màu nâu vàng, độ xốp đất cao (60–65%), giàu mùn hữu cơ (4–5%), ở đai độ cao 300–400 m với lượng mưa bình quân trên 2.000 mm/năm và độ ẩm không khí $>88%$.

2. Tại sao độ lệch $S_k$ về đường kính của Gội trắng ở đai 400 m lại lớn hơn 0?

Chỉ số $S_k > 0$ biểu thị đồ thị phân bố lệch trái so với giá trị trung bình. Điều này phản ánh cấu trúc rừng tự nhiên lành mạnh, trong đó số lượng cây ở cấp kính nhỏ và trung bình chiếm đa số, đảm bảo tính liên tục của quá trình tái sinh và thay thế tầng tán.

3. Loài cây nào là bạn đồng hành tốt nhất của Gội trắng trong mô hình rừng trồng?

Kết quả điều tra lân cận 6 cây khẳng định Vàng anh (Saraca dives) là loài đi kèm tối ưu nhất với tần số xuất hiện $P_0 = 100%$ và tỷ lệ cá thể $P_c = 30.55%$, kế tiếp là Nhãn rừng và Đăng.

4. Thời điểm thu hái hạt giống Gội trắng tốt nhất vào tháng mấy?

Gội trắng ra hoa vào tháng 3–4 và quả chín sinh lý vào tháng 5–6. Cần thu hái quả mọng khi vỏ chuyển sang màu đỏ tươi, hạt có áo hạt màu đỏ sẫm để đạt tỷ lệ nảy mầm cao nhất trong vườn ươm.

5. Biện pháp kỹ thuật nào giúp nâng cao tỷ lệ sống của cây tái sinh Gội trắng?

Cần điều tiết độ tàn che tán rừng duy trì ở mức $0.5 - 0.6$. Tại các lỗ hở có độ tàn che $<0.5$, cần tiến hành phát dọn dây leo, cây bụi thảm tươi xâm lấn nhằm giảm áp lực cạnh tranh ánh sáng và dinh dưỡng khoáng cho cây non.


Kết luận

Công trình nghiên cứu của tác giả Phạm Thị Thu Hà (2018) đã bổ sung cơ sở dữ liệu lâm học định lượng quan trọng về loài Gội trắng (Aphanamixis grandiflora Blume) tại Vườn Quốc gia Cúc Phương. Thông qua hệ thống 6 ô tiêu chuẩn và 30 ô dạng bản, đề tài đã chứng minh Gội trắng phân bố hiệu quả nhất ở đai cao 400 m với mật độ 60 cây/ha, đường kính trung bình đạt 30 cm, độ tàn che ổn định $0.55 - 0.63$ và có mối quan hệ cộng sinh mật thiết với Vàng anh, Nhãn rừng.

Kết quả là nền tảng vững chắc phục vụ chuyển giao kỹ thuật gieo ươm, quy hoạch trồng rừng hỗn giao bản địa và bảo tồn nguồn gen thực vật nguy cấp tại các hệ sinh thái rừng núi đá vôi nhiệt đới Việt Nam. Các đơn vị nghiên cứu và quản lý rừng có thể ứng dụng ngay mô hình tổ thành và giải pháp kỹ thuật nuôi dưỡng để hiện thực hóa các mục tiêu lâm nghiệp bền vững.