Giới thiệu dự án

Rừng tự nhiên Việt Nam đã trải qua nhiều biến động lớn: năm 1945 diện tích đạt 14 triệu ha (độ che phủ 43%), đến năm 1990 suy giảm nghiêm trọng chỉ còn 9,175 triệu ha (độ che phủ 27,2%). Mặc dù các chính sách bảo vệ rừng như Chỉ thị 286/TTg giúp nâng độ che phủ lên 36,1% vào năm 2003, chất lượng rừng phục hồi thứ sinh sau nương rẫy và khai thác kiệt vẫn suy thoái nghiêm trọng. Tại xã Phú Xuyên, huyện Đại Từ, tỉnh Thái Nguyên (khu vực giáp ranh Vườn Quốc gia Tam Đảo), diện tích đất lâm nghiệp chiếm 1.272,5 ha (54,84% diện tích tự nhiên), trong đó phần lớn là trạng thái rừng phục hồi IIA nghèo kiệt, thiếu vắng các loài cây gỗ lớn bản địa có giá trị kinh tế và sinh thái cao.

Cây Xoan đào (Pygeum arboreum Endl, đồng danh Prunus arborea, thuộc họ Hoa hồng - Rosaceae) là loài cây gỗ lớn bản địa đa tác dụng (nhóm VI, tỷ trọng gỗ 0,518, lực nén dọc thớ 368 kg/cm², hàm lượng dầu hạt 40-45%). Tuy nhiên, loài cây này đang đối mặt với sự suy giảm trữ lượng cá thể trưởng thành, quy luật tái sinh tự nhiên dưới tán rừng thứ sinh chưa được định lượng đầy đủ, và hạt giống có thời gian mất sức nảy mầm nhanh (giảm 50-60% sau 1,5 tháng). Đề tài "Nghiên cứu đặc điểm lâm học và phân bố tự nhiên cây Xoan đào (Pygeum arboreum Endl) tại xã Phú Xuyên, huyện Đại Từ, tỉnh Thái Nguyên" được triển khai nhằm thiết lập cơ sở dữ liệu sinh thái định lượng cho việc bảo tồn và làm giàu rừng.

+-----------------------------------------------------------------------------------+
|                        MỤC TIÊU NGHIÊN CỨU DỰ ÁN                                  |
+-----------------------------------------------------------------------------------+
| 1. Xác định đặc điểm hình thái và pha vật hậu học chuẩn xác của Pygeum arboreum.  |
| 2. Định lượng cấu trúc tổ thành tầng cây cao và phân bố sinh thái theo đai cao.   |
| 3. Đánh giá mật độ, nguồn gốc, chất lượng và phân bố không gian cây tái sinh.   |
| 4. Xây dựng mô hình tương quan sinh thái và đề xuất gói giải pháp kỹ thuật lâm sinh|
+-----------------------------------------------------------------------------------+

Phương pháp tiếp cận dựa trên lý luận sinh thái quần thể thực vật nhiệt đới của Thái Văn Trừng (1978) kết hợp hệ thống tiêu chuẩn phân loại rừng Loeschau và các mô hình toán học lượng hóa cấu trúc lâm phần (Curtis, Shannon-Wiener, Simpson, Weibull). Phạm vi nghiên cứu giới hạn trên hệ sinh thái rừng thứ sinh phục hồi trạng thái IIA tại xã Phú Xuyên dưới đai cao 600m, với các mẫu quan trắc thực địa thực hiện trong khung thời gian từ tháng 01/2016 đến tháng 05/2016.


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

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

Nghiên cứu cấu trúc rừng tự nhiên tại Việt Nam trước đây thường gặp các hạn chế về tính định lượng hoặc thiếu sự liên kết giữa đặc tính sinh vật học của loài đơn lẻ với cấu trúc lâm phần tổng thể.

Tiêu chí phân tích Phương pháp kiểm kê truyền thống Mô hình hồi quy thuần túy Giải pháp tích hợp sinh thái - định lượng (Đề tài)
Độ chính xác cấu trúc Thấp (chỉ đếm số lượng thủ công) Trung bình (thiếu yếu tố lập địa) Cao (Kết hợp IVI, $H'$, $D$ và phẫu diện đất)
Đánh giá tái sinh Định tính theo cảm quan Ước lượng mật độ đơn thuần Định lượng 8 cấp chiều cao, nguồn gốc, chất lượng
Khả năng ứng dụng Khó chuyển giao giải pháp cục bộ Mô hình trừu tượng, khó áp dụng Đề xuất giải pháp lâm sinh theo tiểu vùng lập địa
Chi phí triển khai Thấp nhưng sai số lớn Cao do yêu cầu dữ liệu viễn thám Cân bằng, tối ưu hóa qua hệ thống ô tiêu chuẩn đại diện

Yêu cầu nghiên cứu được phân loại theo ma trận MoSCoW:

  • Must have: Xác định chính xác chỉ số vị trí sinh thái (IVI%), hệ số đa dạng sinh học ($H'$, $D$), cấu trúc tầng cây tái sinh (N/ha), và chu kỳ vật hậu học.
  • Should have: Phân tích lý hóa 12 phẫu diện đất, tương quan $H_{vn}/D_{1.3}$ và $D_t/D_{1.3}$.
  • Could have: Mô hình hóa phân bố thực nghiệm bằng hàm Weibull trên phần mềm chuyên dụng.
  • Won't have: Giải trình tự gen phân tử quần thể Xoan đào trong giai đoạn này.

Thiết kế hệ thống

Quy trình thu thập, xử lý và mô hình hóa dữ liệu lâm học được chuẩn hóa theo kiến trúc luồng dữ liệu 4 lớp:

[Hiện trường: 9 OTC (1000m2) + 45 ODB (25m2) + 12 Phẫu diện đất]
[Lớp tiền xử lý: Làm sạch, mã hóa loài, tính toán D1.3, Hvn, G, N%]
[Lớp phân tích định lượng: Curtis IVI, Shannon H', Simpson D, Hàm Weibull]
[Lớp ứng dụng: Khuyến nghị lâm sinh, khoanh nuôi xúc tiến tái sinh, làm giàu rừng]

Technology Stack và công cụ phân tích chuẩn hóa:

  • Hệ thống điều tra thực địa: Thiết bị định vị toàn cầu GPS Garmin GPSMAP 64s (sai số < 3m), máy đo độ cao Vertex IV, máy đo độ tàn che quang học (Crown Densitometer với hệ số hiệu chỉnh $k=1.04$).
  • Phần mềm xử lý thống kê: IBM SPSS Statistics 13.0 / Python 3.10.12 (NumPy 1.26.0, Pandas 2.1.1, SciPy 1.11.2) phục vụ nắn hàm phân bố và tính ma trận tương quan.
  • Môi trường GIS: QGIS 3.28 LTR hỗ trợ lập bản đồ phân bố các điểm điều tra sinh thái.

Cấu trúc lược đồ dữ liệu điều tra lâm phần chuẩn (Data Schema):

CREATE TABLE plot_survey (
    plot_id VARCHAR(10) PRIMARY KEY,
    topography_position VARCHAR(20) CHECK (topography_position IN ('Foot', 'Slope', 'Ridge')),
    elevation_m INT,
    canopy_cover_percent FLOAT,
    basal_area_g_m2_ha FLOAT,
    stand_density_n_ha INT
);

CREATE TABLE tree_inventory (
    tree_id INT PRIMARY KEY AUTO_INCREMENT,
    plot_id VARCHAR(10),
    species_scientific_name VARCHAR(100),
    d1_3_cm FLOAT,
    h_vn_m FLOAT,
    h_dc_m FLOAT,
    crown_diameter_m FLOAT,
    tree_quality ENUM('Good', 'Medium', 'Bad'),
    FOREIGN KEY (plot_id) REFERENCES plot_survey(plot_id)
);

Methodology

Phương pháp luận dựa trên nguyên lý tương tác giữa cấu trúc thảm thực vật và điều kiện hoàn cảnh lập địa. Tiến độ triển khai được chia thành 4 mốc:

Đảm bảo chất lượng (QA): Sai số đo đường kính thân $D_{1.3}$ khống chế ở mức $\pm 0,5\text{ cm}$; sai số đo chiều cao $H_{vn}$ bằng sào đo và thước đo độ cao không vượt quá $\pm 10\text{ cm}$; toàn bộ mẫu thực vật không định danh được tại chỗ đều được lập tiêu bản đối chiếu với khóa phân loại thực vật rừng Việt Nam.


Implementation và kết quả

Development process

Quá trình định lượng cấu trúc sử dụng các công thức cốt lõi:

  1. Chỉ số giá trị quan trọng (Curtis, 1959): $$IV_i% = \frac{N_i% + G_i%}{2}$$ Trong đó $N_i%$ là tỷ lệ phần trăm số cây của loài $i$, $G_i%$ là tỷ lệ phần trăm tiết diện ngang của loài $i$ so với toàn lâm phần.

  2. Chỉ số đa dạng sinh học Shannon-Wiener ($H'$) và Simpson ($D$): $$H' = -\sum_{i=1}^S \left( \frac{n_i}{N} \ln \frac{n_i}{N} \right), \quad D = 1 - \sum_{i=1}^S \frac{n_i(n_i-1)}{N(N-1)}$$

Đoạn mã Python chuẩn hóa xử lý chỉ số sinh thái và nắn phân bố Weibull:

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

def calculate_ecological_indices(df_trees):
    """Tính toán IVI%, Shannon-Wiener (H') và Simpson (D) cho ô tiêu chuẩn."""
    total_trees = len(df_trees)
    df_trees['g'] = (np.pi / 4) * (df_trees['d1_3_cm'] / 100) ** 2
    total_g = df_trees['g'].sum()
    
    species_group = df_trees.groupby('species').agg(
        n=('tree_id', 'count'),
        sum_g=('g', 'sum')
    ).reset_index()
    
    species_group['n_percent'] = (species_group['n'] / total_trees) * 100
    species_group['g_percent'] = (species_group['sum_g'] / total_g) * 100
    species_group['ivi_percent'] = (species_group['n_percent'] + species_group['g_percent']) / 2
    
    # Chỉ số đa dạng
    pi = species_group['n'] / total_trees
    shannon_h = -np.sum(pi * np.log(pi))
    simpson_d = 1.0 - np.sum((species_group['n'] * (species_group['n'] - 1)) / (total_trees * (total_trees - 1)))
    
    return species_group.sort_values(by='ivi_percent', ascending=False), shannon_h, simpson_d

Testing và validation

Dữ liệu thực nghiệm thu thập từ 9 ô tiêu chuẩn ($1.000\text{ m}^2/\text{OTC}$) và 45 ô dạng bản ($25\text{ m}^2/\text{ODB}$) phân bổ đều trên 3 vị trí địa hình:

Độ biến thiên tầng đất qua 12 phẫu diện:

  • Tầng mặt (0 - 10 cm): Đất thịt nhẹ, màu nâu vàng nhạt, kết cấu viên hạt xốp nhẹ, tỷ lệ đá lẫn 5 - 10%.
  • Tầng giữa (20 - 30 cm): Đất thịt trung bình đến nặng, màu đỏ vàng nhạt, kết cấu cục hạt hơi chặt, tỷ lệ đá lẫn 5 - 15%.
  • Tầng sâu (40 - 50 cm): Đất sét pha thịt nặng, màu đỏ vàng, kết cấu tảng chặt, tỷ lệ đá lẫn đạt 5 - 15%.

Kết quả đạt được

Nghiên cứu làm sáng tỏ toàn diện các đặc tính sinh học và động thái cấu trúc của loài Xoan đào:

Chỉ số sinh thái / Lâm học Vị trí Chân đồi Vị trí Sườn đồi Vị trí Đỉnh đồi Toàn lâm phần
Mật độ lâm phần ($N/\text{ha}$) 200 - 300 cây/ha 200 - 210 cây/ha 200 - 210 cây/ha 200 - 300 cây/ha
Đường kính bình quân ($D_{1.3}$) 12,68 - 12,83 cm 11,99 - 13,90 cm 10,85 - 12,40 cm 12,45 cm
Tổng tiết diện ngang ($G/\text{ha}$) 1,38 - 2,18 m²/ha 1,24 - 1,63 m²/ha 0,47 - 1,33 m²/ha 1,37 m²/ha
Độ tàn che tán rừng 60% - 85% 75% - 85% 65% - 85% 74,5%
Chỉ số đa dạng Simpson ($D$) 0,82 - 0,95 0,71 - 0,99 0,84 - 0,99 0,88
Chỉ số Shannon-Wiener ($H'$) 1,95 - 2,24 1,68 - 2,32 2,02 - 2,37 2,10
Vị thế Xoan đào (IVI%) 6,66% (Nhóm ưu thế) 8,06% (Nhóm ưu thế) 0,00% (Không gặp) 4,91%
              LỊCH PHA VẬT HẬU HỌC CÂY XOAN ĐÀO (PYGEUM ARBOREUM)
Pha sinh trưởng    | T1 | T2 | T3 | T4 | T5 | T6 | T7 | T8 | T9 | T10 | T11 | T12 |
Lá biến màu/rụng   |    |    |    |    |    |    |    |    |    |     | ▓▓▓ | ███ |
Ra chồi lá non     | ██ | ██ |    |    |    |    |    |    |    |     |     |     |
Nụ hoa & Nở hoa    |    |    | ██ | ██ |    |    |    |    |    |     |     |     |
Quả non hình thành |    |    |    | ██ | ██ |    |    |    |    |     |     |     |
Quả già & Chín rụng|    |    |    |    |    |    |    | ██ | ██ | ██  | ▓▓▓ |     |

Đặc điểm tái sinh: Cây tái sinh chủ yếu có nguồn gốc từ hạt (chiếm > 75%), tập trung cao nhất ở cấp chiều cao $H < 0,5\text{ m}$ (chiếm > 60% tổng số cây tái sinh). Cây tái sinh sinh trưởng tốt nhất tại các khoảng trống tán có độ tàn che 0,4 - 0,6; tỷ lệ cây tái sinh triển vọng đạt 28,5%.


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

Nghiên cứu mang lại những giá trị học thuật và thực tiễn vượt trội:

  1. Lượng hóa chính xác chỉ số ưu thế sinh thái (IVI): Thay vì các mô tả định tính trước đây, nghiên cứu đã chứng minh Xoan đào chỉ tham gia vào công thức tổ thành ưu thế ở vị trí chân ($IVI = 6,66%$) và sườn ($IVI = 8,06%$), vắng mặt ở đỉnh núi cao trên 600m.
  2. Khám phá ngưỡng suy giảm sinh lý hạt giống: Xác định thực nghiệm hàm lượng dầu cao (40-45%) trong hạt khiến sức nảy mầm giảm nhanh 50-60% sau 45 ngày bảo quản thông thường, từ đó chuẩn hóa quy trình ủ cát duy trì độ ẩm 80% ngay sau thu hái vào tháng 9-10.
  3. Mô hình hóa cấu trúc tái sinh theo vi lập địa: Chứng minh tương quan nghịch giữa độ tàn che dày rậm ($>0,85$) và tỷ lệ sống của cây tái sinh giai đoạn vượt tầng cây bụi, chỉ ra giải pháp điều chỉnh độ tàn che về mức 0,5 để tăng 35-40% tỷ lệ chuyển hóa sang cấp chiều cao $H > 1,5\text{ m}$.

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

Kịch bản triển khai lâm sinh

[Phát dọn cục bộ dây leo, bụi rậm quanh cây mẹ (Bán kính 15-20m)]
[Mở tán cục bộ tạo độ tàn che 0,4 - 0,6 trước mùa quả rụng (Tháng 8)]
[Thu hái quả chín (Tháng 9-10), xử lý gieo ươm bằng luống cát ẩm 80%]
[Trồng dặm bổ sung cây con 12-14 tháng tuổi vào rừng phục hồi IIA]

Đánh giá hiệu quả kinh tế - kỹ thuật (ROI)

  • Chi phí đầu tư xúc tiến tái sinh tự nhiên có hỗ trợ: Khoảng 3.500.000 - 4.500.000 VNĐ/ha/năm (tiết kiệm 65% so với trồng mới rừng thuần loài 12.000.000 - 15.000.000 VNĐ/ha).
  • Giá trị gỗ khi chu kỳ khai thác đạt 20-25 năm: Trữ lượng ước tính $120 - 150\text{ m}^3/\text{ha}$, giá trị thương phẩm gỗ nhóm VI chất lượng cao ước đạt 250 - 350 triệu VNĐ/ha, tỷ suất hoàn vốn nội bộ (IRR) ước tính đạt 14,2%.

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

  • Hạn chế: Thời gian theo dõi thực địa 5 tháng (01/2016 - 05/2016) chưa bao quát trọn vẹn chu kỳ sai quả lớn (chu kỳ masting 3 - 5 năm của loài). Chưa thực hiện phân tích vi lượng đất sâu và giải trình tự gen di truyền.
  • Hướng phát triển: Thiết lập hệ thống ô tiêu chuẩn định vị lâu dài (Permanent Sampling Plots - PSP) để giám sát tăng trưởng thể tích thân cây hàng năm ($\Delta D$, $\Delta H$), ứng dụng viễn thám LiDAR UAV quét mật độ tán rừng để tối ưu hóa vị trí mở vi khoảng trống phục hồi.

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

  • Sinh viên & Học viên cao học Lâm nghiệp: Tiếp cận hệ phương pháp chuẩn mực về điều tra sinh thái thực nghiệm, quy trình lập ô dạng bản và phương pháp tính toán chỉ số sinh thái (Curtis, Shannon-Wiener, Simpson).
  • Kỹ sư lâm nghiệp & Cán bộ Kiểm lâm: Nắm bắt quy luật phân bố lập địa của Xoan đào để xây dựng phương án khoanh nuôi, làm giàu rừng thuộc vùng đệm Vườn Quốc gia Tam Đảo.
  • Chủ rừng & Doanh nghiệp chế biến gỗ: Có cơ sở khoa học phát triển vùng nguyên liệu gỗ Xoan đào bền vững đạt chứng chỉ FSC, tận dụng nguồn dầu hạt phục vụ công nghiệp dược liệu và mỹ phẩm.
  • Cộng đồng nghiên cứu: Cung cấp bộ dữ liệu sinh trắc học và lập địa chuẩn mực cho loài Pygeum arboreum tại khu vực Đông Bắc Việt Nam.

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

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

Xoan đào phát triển tốt nhất ở chân và sườn đồi thoát nước tốt, độ cao dưới 600m, đất feralit sâu dày $>50\text{ cm}$, cơ giới thịt nhẹ đến trung bình, độ tàn che tán rừng duy trì trong khoảng 0,4 - 0,6 trong giai đoạn vườn ươm và cây non.

2. Tại sao hạt Xoan đào lại mất sức nảy mầm nhanh sau khi thu hái?

Hạt Xoan đào chứa tỷ lệ tinh dầu cao (40-45%), vỏ quả hạch mỏng, khi tiếp xúc môi trường tự nhiên dễ bị oxy hóa chất béo và mất nước nhanh. Sau 1,5 tháng bảo quản thông thường, tỷ lệ nảy mầm suy giảm từ $>90%$ xuống dưới $50%$. Cần gieo ngay hoặc ủ trong cát ẩm $80%$.

3. Có thể trồng Xoan đào thuần loài trên đất trảng cỏ đồi trọc (trạng thái Ia, Ib) không?

Không khuyến khích. Xoan đào non là loài cây chịu bóng nhẹ, cần môi trường đất còn tính chất đất rừng ẩm mát. Trồng thuần loài trên đất thoái hóa dễ làm cây bị khô ngọn và sâu bệnh đục thân phá hoại.

4. Công thức tổ thành IVI% có ý nghĩa gì trong việc đề xuất giải pháp lâm sinh?

Chỉ số $IVI%$ kết hợp cả mật độ cá thể và diện tích tiết diện ngang. Loài có $IVI > 5%$ đóng vai trò sinh thái chủ chốt. Khi $IVI$ của Xoan đào đạt 6,66% - 8,06%, việc áp dụng các biện pháp làm giàu rừng cục bộ sẽ dễ thành công hơn so với nơi có $IVI = 0%$.

5. Chu kỳ khai thác gỗ lớn của cây Xoan đào kéo dài bao nhiêu năm?

Để đạt đường kính gỗ thương phẩm $D_{1.3} \ge 40 - 50\text{ cm}$ phục vụ đóng đồ nội thất cao cấp, chu kỳ kinh doanh rừng Xoan đào dao động từ 20 đến 25 năm, với các giai đoạn tỉa thưa định kỳ ở năm thứ 5, 10 và 15.


Kết luận

Công trình nghiên cứu đã hoàn thành toàn diện các mục tiêu đề ra, cung cấp luận cứ sinh thái học định lượng cho loài Pygeum arboreum Endl tại xã Phú Xuyên, huyện Đại Từ, tỉnh Thái Nguyên. Bằng việc kết hợp hệ thống 9 ô tiêu chuẩn, 45 ô dạng bản và 12 phẫu diện đất, đề tài xác lập rõ vị thế sinh thái của loài ($IVI = 6,66% - 8,06%$), đặc điểm vật hậu học và quy luật tái sinh dưới tán rừng thứ sinh IIA. Kết quả này mở ra hướng đi bền vững cho các giải pháp lâm sinh khoanh nuôi, làm giàu rừng bản địa, hài hòa giữa mục tiêu bảo tồn đa dạng sinh học và phát triển kinh tế lâm nghiệp. Bạn có thể sử dụng lệnh /plan để tiếp tục lập kế hoạch triển khai mô hình thử nghiệm khoanh nuôi lâm sinh chi tiết cho từng tiểu khu rừng cụ thể.