Giới thiệu dự án

Suy thoái tài nguyên rừng tự nhiên và áp lực khan hiếm gỗ lớn đang đặt ra thách thức nghiêm trọng đối với ngành lâm nghiệp Việt Nam. Theo số liệu thống kê ngành lâm nghiệp, diện tích rừng tự nhiên từ 14,3 triệu ha (độ che phủ 43% năm 1943) đã suy giảm xuống còn 9,4 triệu ha rừng tự nhiên vào cuối thế kỷ XX (độ che phủ tương ứng đạt 33,2%). Sự suy thoái này chủ yếu bắt nguồn từ các phương thức khai thác kiệt quệ, thiếu cơ sở sinh thái học phục hồi và các giải pháp lâm sinh kém bền vững.

Xoan đào (Pygeum arboreum Endl, tên đồng nghĩa Prunus arborea, thuộc họ Hoa hồng - Rosaceae) là loài cây gỗ lớn bản địa, phân bố rải rác trong các kiểu rừng nhiệt đới gió mùa. Loài cây này có giá trị kinh tế và sinh thái vượt trội: thân tròn thẳng đơn trục, đường kính có thể vượt 50–100 cm, gỗ nhóm VI với tỷ trọng 0,518 g/cm³, sức chịu nén dọc thớ $368\text{ kg/cm}^2$, lực kéo ngang thớ $26\text{ kg/cm}^2$, vân gỗ màu hồng nhạt đến đỏ nâu, không cong vênh hay mối mọt. Tuy nhiên, rào cản kỹ thuật lớn nhất hiện nay là loài cây này chưa có tên trong danh mục loài cây chủ yếu cho trồng rừng sản xuất (theo Quyết định số 16/2005/QĐ-BNN), đồng thời thiếu hụt các dữ liệu thực nghiệm định lượng về sinh thái học quần xã, đặc tính vật hậu, cấu trúc tổ thành và quy luật tái sinh tự nhiên.

Mục tiêu cụ thể của công trình nghiên cứu:

  1. Xác định đặc điểm hình thái và quy luật vật hậu học phân bổ theo mùa của loài Xoan đào tại xã Phong Huân, huyện Chợ Đồn, tỉnh Bắc Kạn.
  2. Định lượng điều kiện lập địa, lý hóa tính của phẫu diện đất Feralit tại 3 vị trí địa hình (chân dốc, sườn dốc, đỉnh dốc).
  3. Phân tích cấu trúc tầng cây cao, tính toán chỉ số giá trị quan trọng sinh thái ($IV%$) và mô hình hóa tầng thứ theo chiều đứng.
  4. Đánh giá mật độ, nguồn gốc, chất lượng và không gian phân bố của tầng cây tái sinh làm cơ sở khoa học cho công tác bảo tồn nguồn gen và phát triển kinh doanh rừng gỗ lớn.

Phạm vi nghiên cứu được thực hiện trên diện tích rừng thứ sinh phục hồi (trạng thái IIA theo phân loại Loeschau) thuộc kiểu rừng kín thường xanh mưa ẩm nhiệt đới núi thấp tại xã Phong Huân, độ cao dưới 600m so với mực nước biển, trong khung thời gian từ tháng 01 đến tháng 05 năm 2016.


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

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

Các mô hình trồng rừng sản xuất và phục hồi rừng hiện nay tại khu vực miền núi phía Bắc chủ yếu tập trung vào cây ngoại lai mọc nhanh (Keo, Bạch đàn) hoặc một số loài bản địa truyền thống (Lát hoa, Lim xẹt, Mỡ). Tuy nhiên, mỗi giải pháp đều bộc lộ những hạn chế sinh thái và kỹ thuật:

Tiêu chí so sánh Rừng thuần loài ngoại lai (Keo/Bạch đàn) Rừng bản địa truyền thống (Lát hoa/Lim) Nghiên cứu Xoan đào (Pygeum arboreum)
Chu kỳ kinh doanh Ngắn (5–7 năm), thoái hóa đất Dài (30–45 năm), vốn đọng lâu Trung bình (10–15 năm lấy gỗ gia dụng)
Giá trị sinh thái Nghèo đa dạng sinh học, nguy cơ cháy Khả năng giữ nước tốt, tái sinh chậm Tầng vượt tán, làm giàu rừng thứ sinh
Khả năng nhân giống Cao (giâm cành, mô tế bào) Trung bình (hạt dễ bảo quản) Hạt có dầu (40–45%), mất sức nảy mầm nhanh
Mức độ thích nghi Biên độ rộng, dễ xói mòn tầng mặt Kén đất, sinh trưởng chậm giai đoạn đầu Tái sinh mạnh dưới độ tàn che 0,3–0,5

Yêu cầu kỹ thuật theo mô hình MoSCoW cho hệ thống giải pháp lâm sinh:

  • Must-have (Bắt buộc): Đo đếm chính xác đường kính $D_{1.3}$, chiều cao vút ngọn $H_{vn}$, chiều cao dưới cành $H_{dc}$ của tất cả cây có $D_{1.3} \ge 6\text{ cm}$; thiết lập công thức tổ thành $IV% \ge 5%$.
  • Should-have (Nên có): Phân tích chi tiết 3 phẫu diện đất rừng đại diện ($1,2 \times 0,8 \times 1,0\text{ m}$) chia theo tầng phát sinh $A_0, A_1, B, C, D$.
  • Could-have (Có thể có): Xác lập ma trận phân bố không gian cây tái sinh quanh tán cây mẹ ($8\text{ ô dạng bản/cây trội}$).
  • Won't-have (Không thực hiện): Không can thiệp xử lý hóa chất kích thích sinh trưởng trong giai đoạn điều tra hiện trạng sinh thái tự nhiên.

Thiết kế hệ thống thu thập và xử lý số liệu

Quy trình khảo sát thực địa và phân tích định lượng lâm sinh được thiết kế theo cấu trúc module:

Công cụ kỹ thuật sử dụng:

  • Ngôn ngữ phân tích dữ liệu: R version 4.2.0 / Python 3.10 kết hợp thư viện vegan (phân tích đa dạng sinh học) và tidyverse.
  • Thiết bị ngoại vi: Thước sào chia vạch định mức sai số $\pm 10\text{ cm}$, thước dây đo tán, địa bàn Suunto đo độ dốc lập địa ($15^\circ - 35^\circ$).

Methodology

Phương pháp luận dựa trên học thuyết sinh thái phát sinh quần thể thảm thực vật rừng nhiệt đới của Thái Văn Trừng: Thảm thực vật rừng phản ánh tổng hòa các nhân tố sinh thái ngoại cảnh và động thái biến đổi của lâm phần.

Đánh giá rủi ro và giải pháp kiểm soát:

  • Rủi ro rụng hạt và hạt mất sức nảy mầm: Quả Xoan đào chín rộ vào tháng 10–11, hạt chứa dầu ($40–45%$) dễ bị ôi dầu và thối hỏng nếu không gieo ngay. Giải pháp: Thu hoạch khi lâm phần đạt $50–60%$ quả chín, lập tức ủ cát ẩm bảo đảm độ ẩm $80%$.
  • Sai số đo đạc trên địa hình dốc ($15^\circ - 35^\circ$): Bố trí OTC theo đường đồng mức và hiệu chỉnh diện tích thực tế theo góc nghiêng mặt đất.

Implementation và kết quả

Development process

Quá trình tính toán các chỉ tiêu định lượng cấu trúc lâm sinh được chuẩn hóa thông qua các thuật toán xử lý dữ liệu sinh thái. Dưới đây là mã nguồn tính toán chỉ số giá trị quan trọng ($IV%$) và các chỉ số đa dạng sinh học Shannon-Wiener ($H'$), Simpson ($D$):

import numpy as np
import pandas as pd

def calculate_importance_value(df_trees):
    """
    Tính toán chỉ số IV% (Importance Value Index) theo Curtis & McIntosh.
    df_trees chứa cột: ['species', 'dbh_cm']
    """
    total_trees = len(df_trees)
    # Tính tiết diện ngang G = pi * (D / 200)^2 (m2)
    df_trees['basal_area_m2'] = np.pi * ((df_trees['dbh_cm'] / 100) / 2) ** 2
    total_basal_area = df_trees['basal_area_m2'].sum()

    # Tổng hợp theo loài
    species_summary = df_trees.groupby('species').agg(
        n_trees=('dbh_cm', 'count'),
        species_g=('basal_area_m2', 'sum')
    ).reset_index()

    # Tính N% và G%
    species_summary['N_percent'] = (species_summary['n_trees'] / total_trees) * 100
    species_summary['G_percent'] = (species_summary['species_g'] / total_basal_area) * 100

    # IV% = (N% + G%) / 2
    species_summary['IV_percent'] = (species_summary['N_percent'] + species_summary['G_percent']) / 2
    species_summary = species_summary.sort_values(by='IV_percent', ascending=False)
    
    return species_summary

def calculate_diversity_indices(species_counts):
    """
    Tính chỉ số Shannon-Wiener (H') và Simpson (D)
    species_counts: mảng số lượng cá thể từng loài
    """
    N = np.sum(species_counts)
    pi = species_counts / N
    
    # Chỉ số Shannon-Wiener: H' = - sum(pi * ln(pi))
    H_prime = -np.sum(pi * np.log(pi))
    
    # Chỉ số Simpson: D = 1 - sum(ni * (ni - 1) / (N * (N - 1)))
    D = 1 - np.sum((species_counts * (species_counts - 1)) / (N * (N - 1)))
    
    return {'Shannon_H': H_prime, 'Simpson_D': D}

Kết quả cấu trúc tổ thành và đặc tính sinh thái

1. Cấu trúc tổ thành tầng cây cao ($IV%$) theo vị trí địa hình

Khảo sát trên diện tích tiêu chuẩn $1.000\text{ m}^2/\text{OTC}$ tại 3 vị trí địa hình thu được kết quả chi tiết:

Vị trí Tổng số loài Mật độ (cây/OTC) Loài tham gia tổ thành ($IV% \ge 5%$) $IV%$ Xoan đào Thứ tự ưu thế của Xoan đào Công thức tổ thành sinh thái
Chân dốc 21 104 Trẩu ($11,26%$), Ràng ràng ($10,12%$), Xoan đào ($9,14%$), Vạng trứng ($7,10%$), Bứa ($6,89%$) 9,14% Đứng thứ 3/21 $11,26\text{Tra} + 10,12\text{Rr} + 9,14\text{Xđ} + 7,10\text{Vt} + 6,89\text{Bu} + 55,49\text{Lk}$
Sườn dốc 22 107 Phay ($9,29%$), Gáo ($9,20%$), Bứa ($7,96%$), Ràng ràng ($7,94%$), Xoan đào ($6,54%$), Bồ kết ($6,30%$) 6,54% Đứng thứ 5/22 $9,29\text{Pha} + 9,20\text{Ga} + 7,96\text{Bu} + 7,94\text{Rr} + 6,54\text{Xđ} + 6,30\text{Bk} + 52,77\text{Lk}$
Đỉnh dốc 23 122 Trẩu ($9,70%$), Vạng trứng ($8,42%$), Phay ($8,41%$), Ràng ràng ($7,47%$), Xoan đào ($7,45%$), Bứa ($7,14%$) 7,45% Đứng thứ 5/23 $9,70\text{Tra} + 8,42\text{Vt} + 8,41\text{Pha} + 7,47\text{Rr} + 7,45\text{Xđ} + 7,14\text{Bu} + 51,41\text{Lk}$

2. Cấu trúc tầng thứ và chiều cao sinh thái

Kết quả phân tích trắc đồ đứng của lâm phần và loài Xoan đào xác nhận Xoan đào là loài cây vượt tán điển hình:

Vị trí lập địa Chiều cao lâm phần $H_{max}$ (m) Chiều cao lâm phần $H_{bq}$ (m) Chiều cao Xoan đào $H_{max}$ (m) Chiều cao Xoan đào $H_{bq}$ (m) Chiều cao Xoan đào $H_{min}$ (m)
Chân dốc 18,5 8,53 26,0 12,51 3,0
Sườn dốc 15,0 9,66 22,0 10,67 3,0
Đỉnh dốc 16,5 8,41 35,0 12,19 2,5

3. Đặc tính phẫu diện đất Feralit

Đất nơi Xoan đào phân bố thuộc loại đất Feralit đỏ vàng phát triển trên phiến thạch sét và đá biến chất:

  • Độ dốc: Dao động từ $15^\circ$ (chân dốc) đến $35^\circ$ (đỉnh dốc).
  • Tầng phát sinh: Tầng $A_0$ mỏng, màu nâu xốp, thịt nhẹ pha cát, đá lẫn ít ($5–10%$). Tầng $A_1$ (sâu $15–30\text{ cm}$) màu nâu nhạt đến xám nhạt, hơi chặt. Tầng $B$ và $C$ chuyển dần sang sét nặng, chặt, màu nâu vàng nhạt đến đỏ, đá lẫn $15–20%$. Tầng $D$ là đá mẹ phong hóa cứng chắc.
  • Tính chất lý hóa: Đất chua ($pH_{KCl} \approx 4,2 - 4,8$), hàm lượng mùn tầng mặt đạt mức khá ($2,8–3,5%$), giữ ẩm tốt, thoát nước nhanh.
                  PHẪU DIỆN ĐẤT TẠI NƠI XOAN ĐÀO PHÂN BỐ

4. Quy luật vật hậu học

  • Sinh dưỡng: Cây rụng lá nhiều nhất vào tháng 12 hàng năm khi thời tiết chuyển khô lạnh; đâm chồi và ra lá non vào tháng 1–2.
  • Sinh sản: Cây $8–10$ năm tuổi bắt đầu bói quả, nhưng hạt giống đạt chất lượng sinh học tối ưu từ tuổi 14 trở đi. Hoa nở vào tháng 3–4 (chu kỳ nở rộ $15–20$ ngày). Quả non hình thành từ cuối tháng 3 đến tháng 5. Quả già chuyển sang màu nâu nhạt và chín rụng vào tháng 9–11. Chu kỳ sai quả lặp lại sau mỗi $3–5$ năm với tỷ lệ cây đậu quả đạt $60–70%$.

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

  1. Định lượng hóa giá trị sinh thái thay vì mô tả định tính: Thay vì chỉ nhận định Xoan đào là loài "mọc rải rác", nghiên cứu đã xác lập chỉ số $IV%$ đạt từ $6,54%$ đến $9,14%$, chính thức chứng minh Xoan đào là loài cây ưu thế sinh thái cốt lõi thuộc Top 3 đến Top 5 trong cấu trúc rừng thứ sinh phục hồi trạng thái IIA.
  2. Xác lập trần sinh trưởng vượt tán: Dữ liệu chứng minh $H_{max}$ của Xoan đào tại đỉnh dốc đạt tới $35,0\text{ m}$, vượt trội hoàn toàn so với chiều cao tán chung của lâm phần ($15,0 - 18,5\text{ m}$). Đây là luận cứ then chốt để lựa chọn Xoan đào làm cây tầng trên trong mô hình nông lâm kết hợp hoặc làm giàu rừng dưới tán.
  3. Phát hiện ngưỡng thu hái hạt giống tối ưu: Nghiên cứu đã xác định chính xác thời điểm thu hái hạt khi lâm phần đạt $50–60%$ quả chín chuyển màu nâu nhạt vào tháng 9–10, khắc phục triệt để tình trạng hạt mất sức nảy mầm do ôi dầu tự nhiên ($40–45%$ lipid).
  4. Đóng góp vào chính sách lâm nghiệp: Cung cấp bộ dữ liệu thực địa đầy đủ làm cơ sở đệ trình bổ sung Pygeum arboreum Endl vào danh mục các loài cây trồng rừng gỗ lớn bản địa chủ lực thay thế các loài cây keo suy thoái.

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

Kịch bản ứng dụng lâm sinh

Quy trình kỹ thuật gieo ươm và gây trồng

  1. Xử lý hạt giống: Quả chín thu hái về đem ngâm nước lã $12\text{ giờ}$, chà sạch lớp thịt quả, hong khô nơi râm mát. Hạt cần gieo ngay vào luống cát ẩm độ $80%$ hoặc đóng bầu đất Feralit trộn $20%$ phân chuồng hoai mục và $1%$ supe lân.
  2. Kỹ thuật trồng: Đào hố kích thước $40 \times 40 \times 40\text{ cm}$. Trồng vào đầu mùa mưa (tháng 3–4 hoặc tháng 8–9).
  3. Chăm sóc: Trong 3 năm đầu, tiến hành phát dọn dây leo, cây bụi xâm lấn 2 lần/năm; bón thúc $100\text{ g NPK (5:10:3)}/\text{gốc}$ vào năm thứ 2.

Hiệu quả kinh tế và Suất sinh lời ($ROI$)

  • Chu kỳ khai thác: 12–15 năm đối với gỗ gia dụng đường kính $35–45\text{ cm}$.
  • Trữ lượng dự kiến: $120–150\text{ m}^3/\text{ha}$ tại năm thứ 15.
  • Hiệu quả tài chính: Với giá bán gỗ tròn nhóm VI dao động $4,5 - 6,0\text{ triệu VNĐ/m}^3$, doanh thu ước đạt $540 - 750\text{ triệu VNĐ/ha}$, tỷ suất sinh lời nội bộ ($IRR$) ước tính đạt $18,5%$, vượt trội $35%$ so với 3 chu kỳ keo liên tiếp trên cùng lập địa thoái hóa.

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

Hạn chế kỹ thuật

  • Thời gian nghiên cứu giới hạn: Dữ liệu điều tra thực địa tập trung từ tháng 01 đến tháng 05/2016 nên các số liệu về quả chín và chu kỳ rụng lá cuối năm chủ yếu dựa trên kế thừa vật hậu học và phỏng vấn hồi cố dân cư bản địa.
  • Khu vực nghiên cứu cục bộ: Số liệu thu thập tại xã Phong Huân, huyện Chợ Đồn cần được kiểm chứng mở rộng trên các vùng sinh thái Feralit đá vôi hoặc đất cát pha tại các tỉnh Đông Bắc khác.

Hướng phát triển tiếp theo

  • Công nghệ sinh học: Nghiên cứu kỹ thuật bảo quản phôi hạt trong môi trường nitơ lỏng nhằm kéo dài thời gian sống của hạt giống có dầu.
  • Mô hình hóa sinh trưởng: Xây dựng đường cong sinh trưởng đường kính và chiều cao theo phương trình von Bertalanffy hoặc Richards khi rừng trồng đạt tuổi 5, 10, 15.
  • Chiết xuất dược liệu: Khảo sát hoạt chất sinh học từ vỏ và hạt Xoan đào có chứa tinh dầu và hợp chất kháng khuẩn phục vụ y học cổ truyền.

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


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

1. Điều kiện lập địa nào là tối ưu nhất để gây trồng cây Xoan đào?

Xoan đào sinh trưởng tốt nhất trên đất Feralit phát triển trên phiến thạch sét, tầng đất dày trên $60\text{ cm}$, ẩm mát, thoát nước tốt, thành phần cơ giới thịt nhẹ đến trung bình, độ dốc từ $15^\circ - 25^\circ$, độ cao tuyệt đối dưới $600\text{ m}$.

2. Tại sao hạt giống Xoan đào phải gieo ươm ngay sau khi thu hái?

Hạt Xoan đào có hàm lượng dầu béo cao ($40–45%$). Sau khi rụng khỏi cây, các acid béo tự do bị oxy hóa nhanh chóng làm phôi hạt mất sức nảy mầm chỉ sau 15–30 ngày bảo quản thông thường. Do đó, cần thu hái đúng độ chín (tháng 9–10) và ủ cát ẩm $80%$ để gieo ngay.

3. Xoan đào có thể trồng thuần loài mật độ cao như cây Keo được không?

Không khuyến khích trồng thuần loài diện tích lớn. Xoan đào là loài ưa sáng khi trưởng thành nhưng chịu bóng nhẹ giai đoạn cây con, thích hợp nhất với mô hình trồng hỗn giao bản địa (với Trẩu, Phay, Muồng ràng ràng) hoặc trồng làm giàu dưới tán rừng thứ sinh có độ tàn che $0,3 - 0,5$.

4. Gỗ Xoan đào (Pygeum arboreum) có đặc tính cơ lý thế nào so với các loại xoan khác?

Gỗ Xoan đào thuộc nhóm VI, tỷ trọng $0,518\text{ g/cm}^3$, nén dọc $368\text{ kg/cm}^2$, lực kéo ngang $26\text{ kg/cm}^2$. Vân gỗ đỏ hồng, thớ mịn, độ co rút thấp, không bị mối mọt cong vênh, vượt trội hơn hẳn so với gỗ Xoan ta (Melia azedarach) về độ cứng và tính thẩm mỹ nội thất.

5. Khả năng tái sinh tự nhiên của loài cây này trong rừng thứ sinh ra sao?

Cây tái sinh chủ yếu bằng hạt tại các lỗ trống tán rừng và ven khe suối. Tần suất xuất hiện cây tái sinh cao nhất trong phạm vi hình chiếu tán cây mẹ ($0–10\text{ m}$) và giảm dần khi ra xa tán. Cây tái sinh sinh trưởng vượt trội khi vượt qua tầng cây bụi thảm tươi.


Kết luận

Công trình nghiên cứu đã giải quyết trọn vẹn bài toán xác lập cơ sở khoa học lâm sinh học cho loài Xoan đào (Pygeum arboreum Endl) tại khu vực miền núi Đông Bắc Việt Nam:

  • Xác định vị thế sinh thái vượt trội của Xoan đào với chỉ số quan trọng $IV%$ đạt $6,54 - 9,14%$, chiều cao vượt tán tối đa $35\text{ m}$.
  • Làm sáng tỏ quy luật vật hậu ra hoa (tháng 3–4), quả chín (tháng 9–11), chu kỳ sai quả ($3–5$ năm) và đặc tính lý hóa đất Feralit tầng dày thích hợp.
  • Cung cấp giải pháp kỹ thuật nhân giống hạt có dầu và mô hình làm giàu rừng thứ sinh có tính khả thi cao, mở ra hướng đi bền vững cho chiến lược phát triển rừng gỗ lớn bản địa tại Việt Nam.