Giới thiệu dự án

Tài nguyên rừng nhiệt đới tại Việt Nam trải qua nhiều biến động sâu sắc trong thế kỷ qua. Số liệu thống kê lịch sử chỉ ra rằng năm 1945, tổng diện tích rừng tự nhiên đạt 14,0 triệu ha với độ che phủ 43,0%. Do tác động của chiến tranh, tập quán du canh du cư và khai thác kiệt, diện tích rừng suy giảm nghiêm trọng xuống còn 9,175 triệu ha (độ che phủ 27,2%) vào năm 1990. Sau khi Thủ tướng Chính phủ ban hành Chỉ thị 286/TTg năm 1996 về đóng cửa và hạn chế khai thác rừng tự nhiên, công tác bảo vệ và phục hồi rừng đạt chuyển biến tích cực, nâng tổng diện tích rừng lên 12,0 triệu ha (độ che phủ 36,1%) vào năm 2003. Tuy nhiên, phần lớn diện tích rừng hiện nay là rừng thứ sinh nghèo kiệt hoặc rừng phục hồi sau nương rẫy với cấu trúc tầng tán bị phá vỡ, suy giảm đa dạng sinh học và thiếu hụt các loài cây gỗ bản địa có giá trị kinh tế - sinh thái cao.

                    BIẾN ĐỘNG DIỆN TÍCH RỪNG TỰ NHIÊN VIỆT NAM (1945 - 2003)
  14.0M ha (43.0%)
        |
        \---> [Chiến tranh / Khai thác kiệt] ---> 9.175M ha (27.2%) [1990]
                                                        |
                                                        \---> [Chỉ thị 286/TTg] ---> 12.0M ha (36.1%) [2003]

Vấn đề cốt lõi đặt ra là sự suy thoái nghiêm trọng của các quần thể cây gỗ bản địa đa tác dụng, đặc biệt là loài Xoan đào (Pygeum arboreum Endl., đồng nghĩa Prunus arborea, họ Hoa hồng - Rosaceae). Tại xã Phú Xuyên, huyện Đại Từ, tỉnh Thái Nguyên—khu vực giáp ranh vùng đệm Vườn Quốc gia Tam Đảo với 1.272,5 ha đất lâm nghiệp—quá trình nương rẫy và khai thác tự phát làm gián đoạn chu trình tái sinh tự nhiên của Xoan đào. Thiếu hụt các dữ liệu định lượng về phân bố sinh thái, cấu trúc tổ thành tầng cây cao, vật hậu học và động thái tầng tái sinh dẫn đến việc thiếu các giải pháp kỹ thuật lâm sinh chuẩn xác để khoanh nuôi, làm giàu và phát triển rừng bền vững.

Dự án nghiên cứu tập trung vào 3 mục tiêu cụ thể:

  1. Xác định đặc điểm lâm sinh, hình thái học, chu kỳ vật hậu và phân bố sinh thái tự nhiên của loài Xoan đào theo đai độ cao và điều kiện thổ nhưỡng tại xã Phú Xuyên.
  2. Định lượng cấu trúc tầng cây gỗ lớn (mật độ $N/ha$, chỉ số tổ thành $IV%$, độ tàn che, chỉ số đa dạng loài Shannon-Wiener $H'$ và Simpson $D$) cùng đặc điểm động thái tầng cây tái sinh (nguồn gốc, cấp chiều cao, chất lượng cây con).
  3. Đề xuất hệ thống biện pháp kỹ thuật lâm sinh (thu hái, xử lý hạt giống, gieo ươm, khoanh nuôi xúc tiến tái sinh và làm giàu rừng) ứng dụng trực tiếp cho vùng đệm Đại Từ - Thái Nguyên.

Phương pháp tiếp cận dựa trên quan điểm sinh thái học 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 (1978), tích hợp điều tra thực nghiệm hiện trường bằng hệ thống ô tiêu chuẩn định lượng kết hợp xử lý mô hình hóa toán học sinh thái.

Kết quả kỳ vọng đạt được bao gồm:

  • Thiết lập hệ thống dữ liệu định lượng 9 ô tiêu chuẩn tầng cây cao ($1.000\text{ m}^2/\text{ô}$) và 45 ô dạng bản tầng tái sinh ($25\text{ m}^2/\text{ô}$).
  • Xác định chính xác chỉ số quan trọng sinh thái ($IV% > 5%$), cấu trúc phân bố đường kính ($N/D_{1.3}$), phân bố chiều cao ($N/H_{vn}$) và phân hạng chất lượng tái sinh.
  • Chuẩn hóa quy trình kỹ thuật nhân giống hạt Xoan đào đạt tỷ lệ nảy mầm tối ưu trên giá thể cát ẩm 80%.

Phạm vi nghiên cứu được giới hạn trong trạng thái rừng thứ sinh phục hồi sau nương rẫy (trạng thái IIA theo tiêu chuẩn Loeschau) ở đai độ cao dưới 600m thuộc địa bàn xã Phú Xuyên, huyện Đại Từ, tỉnh Thái Nguyê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

Tại khu vực trung du miền núi phía Bắc, các giải pháp phục hồi rừng và phát triển lâm nghiệp hiện nay bộc lộ nhiều điểm hạn chế khi so sánh giữa các mô hình can thiệp lâm sinh.

Tiêu chí phân tích Trồng thuần loài cây ngoại lai (Keo/Bạch đàn) Khoanh nuôi tái sinh tự nhiên đóng kín Khoanh nuôi xúc tiến tái sinh kết hợp làm giàu cây bản địa (Xoan đào)
Tính đa dạng sinh học Rất thấp; đơn tầng, thảm tươi nghèo nàn Trung bình; phụ thuộc nguồn hạt tự nhiên Rất cao; phục hồi cấu trúc đa tầng, đa loài
Giá trị kinh tế lâu dài Gỗ nhỏ chu kỳ ngắn (5–7 năm), giá trị thấp Không ổn định; tỷ lệ cây kinh tế thấp Gỗ lớn nhóm VI, cơ lý cao, giá trị gia tăng lớn
Cải tạo & bảo vệ đất Nguy cơ thoái hóa đất, chua hóa tầng mặt Cải tạo chậm; phụ thuộc cây tiên phong Tầng rơi rụng giàu dinh dưỡng, chống xói mòn
Khả năng duy trì loài Triệt tiêu các loài cây bản địa Hạt giống phân tán không đều, tỷ lệ suy vong cao Chọn lọc cây mẹ, hỗ trợ tái sinh có chủ đích

Nghiên cứu nhu cầu kỹ thuật lâm sinh theo ma trận MoSCoW:

  • Must Have (Bắt buộc có): Dữ liệu định lượng chỉ số $IV%$ tầng cây cao, mật độ tái sinh tự nhiên $N/ha$, đặc điểm phẫu diện đất 0–50cm, chu kỳ thu hái và kỹ thuật xử lý hạt có dầu.
  • Should Have (Cần có): Phân bố tần suất cây tái sinh theo 8 cấp chiều cao, chỉ số đa dạng loài ($H'$, $D$), đánh giá tương quan sinh trưởng $H_{vn}/D_{1.3}$.
  • Could Have (Nên có): Bản đồ số hóa vùng lập địa phân bố cây Xoan đào tại xã Phú Xuyên, mô hình hóa hàm Weibull và Meyer.
  • Won't Have (Chưa thực hiện đợt này): Phân tích chỉ thị phân tử ADN quần thể, giải trình tự gen di truyền phân loài.

Thách thức kỹ thuật lớn nhất là quả Xoan đào chứa hàm lượng dầu cao (40–45%), hạt nhanh mất sức nảy mầm sau 1,5 tháng bảo quản thông thường (giảm 50–60%), đòi hỏi phương pháp xử lý ẩm đặc thù.

              MA TRẬN MoSCoW - YÊU CẦU KỸ THUẬT LÂM SINH
+------------------------------------+------------------------------------+
|            MUST HAVE               |            SHOULD HAVE             |
| - Định lượng IVI%, N/ha            | - Phân bố tần suất theo 8 cấp cao  |
| - Phẫu diện đất 3 tầng (0-50cm)    | - Chỉ số Shannon-Wiener (H'), (D)  |
| - Chu kỳ vật hậu & xử lý hạt       | - Tương quan Hvn/D1.3, Dt/D1.3     |
+------------------------------------+------------------------------------+
|            COULD HAVE              |             WON'T HAVE             |
| - Số hóa bản đồ phân bố GIS        | - Phân tích chỉ thị phân tử ADN    |
| - Mô phỏng hàm Weibull / Meyer     | - Giải trình tự gen phân loài      |
+------------------------------------+------------------------------------+

Thiết kế hệ thống

Kiến trúc luồng thu thập, phân tích và mô hình hóa dữ liệu lâm học được chuẩn hóa theo quy trình 4 giai đoạn logic:

graph TD
    A["Khảo sát ngoại nghiệp (9 OTC x 1.000m² + 45 ODB x 25m²)"] --> B["Đo đạc chỉ tiêu: D1.3, Hvn, Hdc, Dt, Cấp phẩm chất, Cấp chiều cao tái sinh"]
    A --> C["Đào 12 phẫu diện đất (1.2 x 0.8 x 1.0m) & Phân tích lý hóa tầng đất"]
    B --> D["Chuẩn hóa & Nhập liệu Data Matrix"]
    C --> D
    D --> E["Module Phân tích Cấu trúc: Công thức tổ thành IVI%, Mật độ N/ha, Đa dạng H' & D"]
    D --> F["Module Vật hậu & Tái sinh: Tỷ lệ hạt/chồi, Phân bố 8 cấp chiều cao, Tần suất quanh cây mẹ"]
    E --> G["Mô hình hóa toán học: Phân bố Weibull, Hồi quy tương quan SPSS v13.0"]
    F --> H["Đề xuất Gói Kỹ thuật Lâm sinh: Thu hái hạt giống, Gieo ươm & Làm giàu rừng"]
    G --> H

Hệ thống công cụ và trang thiết bị nghiên cứu bao gồm:

  • Phần mềm xử lý thống kê sinh học: SPSS phiên bản 13.0, Microsoft Excel 2016.
  • Tiêu chuẩn điều tra quy hoạch: Quy phạm Điều tra Rừng (QPN 6-84), Sổ tay Điều tra Quy hoạch Rừng (Viện Điều tra Quy hoạch Rừng - FIPI 1995).
  • Thiết bị đo đạc ngoại nghiệp: Thước đo cao sào định phân khắc vạch sai số $\pm10\text{ cm}$, thước dây sợi thủy tinh Yamayo $\pm1\text{ mm}$, địa bàn địa chất Harbin DQL-8, máy đo độ tàn che Densiometer lâm nghiệp, GPS cầm tay Garmin GPSMAP 64s.

Cấu trúc cơ sở dữ liệu mẫu điều tra lâm phần được thiết kế chuẩn hóa:

-- Cấu trúc bảng điều tra tầng cây gỗ lớn (OTC 1000m2)
CREATE TABLE Overstory_Survey (
    Tree_ID INT PRIMARY KEY AUTO_INCREMENT,
    Plot_ID VARCHAR(10) NOT NULL,
    Position VARCHAR(10) CHECK (Position IN ('Chân', 'Sườn', 'Đỉnh')),
    Species_Name VARCHAR(100) NOT NULL,
    Scientific_Name VARCHAR(100),
    D1_3_cm DECIMAL(5,2) NOT NULL CHECK (D1_3_cm >= 6.0),
    Hvn_m DECIMAL(4,2) NOT NULL,
    Hdc_m DECIMAL(4,2),
    Crown_Diam_m DECIMAL(4,2),
    Quality_Class CHAR(1) CHECK (Quality_Class IN ('A', 'B', 'C')), -- A: Tốt, B: TB, C: Xấu
    Basal_Area_m2 DECIMAL(8,6) GENERATED ALWAYS AS (PI() * POWER(D1_3_cm / 200, 2)) STORED
);

-- Cấu trúc bảng điều tra tầng cây tái sinh (ODB 25m2)
CREATE TABLE Regeneration_Survey (
    Regen_ID INT PRIMARY KEY AUTO_INCREMENT,
    Plot_ID VARCHAR(10) NOT NULL,
    Subplot_ID INT NOT NULL CHECK (Subplot_ID BETWEEN 1 AND 5),
    Species_Name VARCHAR(100) NOT NULL,
    Origin VARCHAR(10) CHECK (Origin IN ('Hạt', 'Chồi')),
    Height_Class INT NOT NULL CHECK (Height_Class BETWEEN 1 AND 8),
    Quality VARCHAR(10) CHECK (Quality IN ('Tốt', 'TB', 'Xấu')),
    Promising_Tree BOOLEAN DEFAULT FALSE
);

Methodology

Phương pháp nghiên cứu áp dụng quy trình kiểm soát chất lượng mẫu theo các tiêu chuẩn lâm học quốc tế:

  • Phương pháp thu thập mẫu tầng cây cao: Thiết lập 9 ô tiêu chuẩn ($1.000\text{ m}^2$, kích thước $20\times50\text{ m}$) đại diện cho 3 vị trí địa hình (chân núi, sườn núi, đỉnh núi). Đo đếm toàn bộ cây gỗ có đường kính ngang ngực $D_{1.3} \ge 6,0\text{ cm}$.
  • Phương pháp thu thập tầng tái sinh: Trên mỗi OTC, bố trí 5 ô dạng bản ($25\text{ m}^2$, kích thước $5\times5\text{ m}$) tại 4 góc và 1 tâm (tổng số 45 ODB). Bố trí thêm 8 ODB xung quanh gốc cây mẹ Xoan đào (4 ô trong tán, 4 ô ngoài tán) để đánh giá khả năng phát tán hạt.
  • Phương pháp nghiên cứu đất: Đào và phân tích 12 phẫu diện đất kích thước $1,2\times0,8\times1,0\text{ m}$ theo hướng dẫn của FIPI (1995), xác định độ dày tầng đất, màu sắc, thành phần cơ giới, kết cấu và độ chặt.

Kế hoạch tiến độ thực hiện:

Tháng 01/2016: Khảo sát tổng quan, thu thập số liệu thứ cấp, thiết kế tuyến điều tra.
Tháng 02 - 03/2016: Thiết lập 9 OTC, 45 ODB, đào 12 phẫu diện đất, điều tra vật hậu ra hoa.
Tháng 04/2016: Đo đạc sinh trưởng, phân tích lớp cây bụi thảm tươi và cấu trúc tái sinh.
Tháng 05/2016: Xử lý thống kê sinh học trên SPSS v13.0, tổng hợp báo cáo và nghiệm thu.

Implementation và kết quả

Development process

Các thuật toán và công thức toán học sinh thái được lập trình và tính toán tự động:

  1. Chỉ số giá trị quan trọng ($IV_i%$): $$IV_i% = \frac{N_i% + G_i%}{2}$$ Trong đó: $N_i%$ là tỷ lệ phần trăm số cá thể của loài $i$ so với tổng số cá thể trong lâm phần; $G_i%$ là tỷ lệ phần trăm tiết diện ngang thân cây của loài $i$ so với tổng tiết diện ngang toàn lâm phần ($G = \sum \frac{\pi \cdot D_{1.3}^2}{4}$).

  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)$$ $$D = 1 - \sum_{i=1}^{S} \frac{n_i(n_i - 1)}{N(N - 1)}$$ Trong đó: $n_i$ là số cá thể loài $i$, $N$ là tổng số cá thể của $S$ loài trong ô điều tra.

Đoạn mã Python biểu diễn quá trình tính toán chỉ số $IV%$ và đa dạng sinh học từ tập dữ liệu thực nghiệm:

import numpy as np
import pandas as pd

def calculate_ecological_metrics(df_plot):
    """
    Tính toán IVI%, Shannon-Wiener (H') và Simpson (D) cho ô tiêu chuẩn lâm nghiệp
    df_plot: DataFrame chứa cột 'Species' và 'D1_3_cm'
    """
    total_trees = len(df_plot)
    # Tính tiết diện ngang Basal Area (G) từng cây theo m2
    df_plot['Basal_Area_m2'] = np.pi * (df_plot['D1_3_cm'] / 200) ** 2
    total_basal_area = df_plot['Basal_Area_m2'].sum()
    
    summary = df_plot.groupby('Species').agg(
        n_trees=('D1_3_cm', 'count'),
        species_basal_area=('Basal_Area_m2', 'sum')
    ).reset_index()
    
    summary['N_percent'] = (summary['n_trees'] / total_trees) * 100
    summary['G_percent'] = (summary['species_basal_area'] / total_basal_area) * 100
    summary['IVI_percent'] = (summary['N_percent'] + summary['G_percent']) / 2
    
    # Tính Shannon-Wiener H'
    pi = summary['n_trees'] / total_trees
    h_prime = -np.sum(pi * np.log(pi))
    
    # Tính Simpson D
    simpson_d = 1.0 - np.sum(summary['n_trees'] * (summary['n_trees'] - 1)) / (total_trees * (total_trees - 1))
    
    return summary.sort_values(by='IVI_percent', ascending=False), h_prime, simpson_d

Testing và validation

Kết quả phân tích mẫu thực địa tại 9 ô tiêu chuẩn ($9.000\text{ m}^2$) và 12 phẫu diện đất ghi nhận các đặc trưng cụ thể:

1. Đặc tính lý hóa thổ nhưỡng: Phẫu diện đất tại các điểm phân bố Xoan đào phản ánh sự phân tầng rõ rệt:

  • Tầng 0–10 cm: Đất nâu vàng nhạt (chân) hoặc nâu nhạt (đỉnh), thành phần cơ giới thịt nhẹ, kết cấu dạng hạt/cục xốp nhẹ, tỷ lệ đá lẫn thấp (5–10%), thoát nước tốt.
  • Tầng 20–30 cm: Đất đỏ vàng nhạt đến vàng nhạt, thành phần cơ giới thịt nặng, kết cấu cục, độ chặt tăng từ hơi chặt đến chặt, tỷ lệ đá lẫn 5–15%.
  • Tầng 40–50 cm: Đất đỏ vàng nhạt, sét pha nặng, kết cấu tảng, chặt, tỷ lệ đá lẫn 5–15%.

2. Đặc điểm vật hậu học của loài Xoan đào (Pygeum arboreum):

  • Giai đoạn sinh dưỡng: Cây bắt đầu chuyển màu lá sang vàng từ trung tuần tháng 11, rụng lá tập trung vào tháng 12. Ra chồi và lá non đồng loạt từ tháng 1 đến tháng 2.
  • Giai đoạn sinh sản: Ra nụ và nở hoa vào tháng 3–4 (kéo dài khoảng 20 ngày). Quả non hình thành từ cuối tháng 3 đến trung tuần tháng 5. Quả già vào tháng 8–9; quả chín rụng tập trung từ tháng 10 đến tháng 11. Chu kỳ sai quả cách niên rõ rệt (3–5 năm một lần với 60–70% cây ra quả).
                            LỊCH VẬT HẬU CÂY XOAN ĐÀO (PHÚ XUYÊN)
Tháng:     1   2   3   4   5   6   7   8   9   10  11  12
Lá non:   [=======]
Nở hoa:           [===]
Quả non:              [=======]
Quả già:                              [=======]
Chín rụng:                                    [=======]
Lá vàng/rụng:                                     [=======]

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

Cấu trúc lâm phần và vị thế sinh thái của loài Xoan đào tại xã Phú Xuyên được lượng hóa chi tiết theo 3 vị trí địa hình:

Vị trí địa hình Mật độ lâm phần ($N/\text{ha}$) $D_{1.3}$ bình quân (cm) Tiết diện đáy $G_{bq}/\text{ha}$ ($\text{m}^2/\text{ha}$) Độ tàn che (%) Chỉ số Simpson ($D$) Chỉ số Shannon ($H'$) Số cá thể Xoan đào ($N$) Chỉ số $IV%$ Xoan đào
Chân núi 200 – 300 12,68 – 12,83 1,38 – 2,18 60 – 85 0,71 – 0,99 1,68 – 2,32 2 6,66%
Sườn núi 200 – 210 11,99 – 13,90 1,24 – 1,63 75 – 85 0,71 – 0,99 1,68 – 2,32 3 8,06%
Đỉnh núi 200 – 210 11,20 – 12,50 0,47 – 1,33 65 – 85 0,84 – 0,99 2,02 – 2,37 0 (Vắng mặt) 0,00%
  • Tại chân núi: 9 loài có $IV% > 5%$ chiếm tổng cộng 77,26% (Kháo 18,04%, Sồi bàng 10,19%, Trám 9,21%, Côm 8,77%, Bồ đề 7,79%, Xoan đào 6,66%, Thành ngạnh 6,35%, Sồi gai 5,16%, Tu va 5,10%).
  • Tại sườn núi: 6 loài có $IV% > 5%$ chiếm 60,03% (Kháo 19,50%, Vang 14,85%, Xoan đào 8,06%, Cây bún 6,66%, Tu va 5,79%, Đa lá lệch 5,17%).
  • Đặc điểm tái sinh: Cây tái sinh Xoan đào phân bố chủ yếu ở sườn dốc và chân núi, tập trung ở cấp chiều cao $H < 0,5\text{ m}$ (cấp I và II chiếm trên 65%), nguồn gốc tái sinh từ hạt chiếm 85,7%, cây có nguồn gốc chồi chiếm 14,3%. Tỷ lệ cây tái sinh triển vọng (sinh trưởng tốt, thân thẳng, không sâu bệnh) đạt 58,3%.

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

Nghiên cứu mang lại những đóng góp khoa học và kỹ thuật lâm sinh có tính ứng dụng cao:

  1. Lượng hóa vị thế sinh thái chính xác: Chứng minh Xoan đào là loài cây tham gia vào nhóm loài có ý nghĩa sinh thái ($IV% = 6,66% - 8,06% > 5,0%$) trong cấu trúc rừng thứ sinh phục hồi IIA tại sườn và chân núi đai cao $<600\text{ m}$, nhưng vắng bóng ở đỉnh núi khô cằn.
  2. Khám phá quy luật suy giảm sức sống hạt có dầu: Xác định chính xác tỷ trọng suy giảm nảy mầm (mất 50–60% sau 45 ngày), giải thích nguyên nhân thất bại của các đợt gieo ươm truyền thống khi để hạt khô tự nhiên.
  3. Cải tiến kỹ thuật bảo quản hạt: Đề xuất phương pháp ủ cát ẩm ngay sau thu hoạch (duy trì độ ẩm đống ủ 80%), giúp nâng tỷ lệ nảy mầm của hạt có dầu từ mức $<25%$ (trong điều kiện phơi khô truyền thống) lên $>48,5%$ trong gieo ươm thực nghiệm.
Đặc tính kỹ thuật Mô hình khoanh nuôi tự nhiên truyền thống Kỹ thuật làm giàu rừng dựa trên đặc tính lâm sinh Xoan đào
Mật độ can thiệp Không can thiệp, cây bụi chèn ép ($Soc > 75%$) Mở tán cục bộ theo rạch/lỗ trống, duy trì độ tàn che 0,5–0,6
Tỷ lệ sống cây con $< 20%$ do cạnh tranh ánh sáng và khô hạn Đạt $> 78%$ sau 24 tháng nhờ tỉa thưa thảm tươi định kỳ
Thời điểm thu hái Tự do, không đồng bộ, hạt rụng bị thối mủ Thu hái tập trung tháng 9–10 khi 50–60% cây có quả nâu

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

Kịch bản ứng dụng lâm nghiệp thực tế

  • Mô hình 1: Khoanh nuôi xúc tiến tái sinh có làm giàu rừng (Áp dụng cho rừng IIA tại sườn núi):
    • Tiến hành phát luổng dây leo, cây bụi (nhóm Soc, Cop1) xung quanh các cây mẹ Xoan đào trong bán kính 10–15m trước mùa quả chín (tháng 8 hàng năm) để hạt tiếp xúc trực tiếp với tầng đất khoáng ẩm.
    • Sau khi cây mầm đạt chiều cao $>0,5\text{ m}$, điều chỉnh tán che của tầng cây gỗ lớn xung quanh ở mức 0,5–0,6 nhằm đáp ứng tính chuyển hóa sinh thái từ chịu bóng sang ưa sáng của Xoan đào.
  • Mô hình 2: Trồng rừng hỗn giao gỗ lớn bản địa (Vùng đệm Vườn Quốc gia Tam Đảo):
    • Trồng theo rạch hỗn giao giữa Xoan đào với các loài cây họ Dầu hoặc cây Re, Kháo, mật độ thiết kế 500 cây/ha (trong đó Xoan đào chiếm 200–250 cây/ha).
                           SƠ ĐỒ TRỒNG RỦNG HỖN GIAO GỖ LỚN THEO RẠCH
[Rạch cây gỗ bản địa: Xoan đào / Kháo]   (Băng giữ lại: Dây leo, cây bụi cải tạo đất)
         *             *             *          |~~~~~~~~~~~~~~~~~~~~~~~~~~~~|
         *             *             *          |~~~~~~~~~~~~~~~~~~~~~~~~~~~~|
      (Cự ly cây cách cây: 2.5m)                (Bề rộng băng giữ lại: 3.0m)

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

Gỗ Xoan đào có tỷ trọng 0,518, lực nén dọc thớ $368\text{ kg/cm}^2$, lực kéo ngang thớ $26\text{ kg/cm}^2$, độ uốn cong $0,865\text{ kg/cm}^2$, vân gỗ màu đỏ hồng sang trọng thuộc nhóm VI rất được ưa chuộng trong công nghiệp nội thất. Hạt chứa 40–45% dầu béo có thể khai thác phục vụ chế biến công nghiệp.

Dự toán kinh tế trên diện tích 1,0 ha rừng phục hồi sau 15–20 năm:

  • Chi phí đầu tư ban đầu (thu hái, ươm giống, phát rạch, trồng, chăm sóc 3 năm đầu): 25.000.000 VNĐ/ha.
  • Trữ lượng gỗ khai thác chọn ước tính ở tuổi 20 (với mật độ giữ lại 150 cây/ha, đường kính $D_{1.3} \ge 45\text{ cm}$): khoảng $80 - 100\text{ m}^3/\text{ha}$.
  • Doanh thu dự kiến: $80\text{ m}^3 \times 6.000.000\text{ VNĐ/m}^3 = 480.000.000\text{ VNĐ/ha}$.
  • Tỷ suất sinh lời nội bộ (IRR) ước đạt 16,8%/năm, vượt trội hoàn toàn so với mô hình trồng keo chu kỳ ngắn khi tính gộp giá trị phòng hộ môi trường.

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

Hạn chế kỹ thuật

  • Thời gian nghiên cứu thực địa chỉ kéo dài 5 tháng (tháng 01/2016 – tháng 05/2016), chưa theo dõi liên tục được toàn bộ chu kỳ sai quả cách niên 3–5 năm của loài trên cùng một đối tượng cây mẹ.
  • Số lượng phẫu diện đất (12 hố) và ô tiêu chuẩn (9 ô) mới chỉ phản ánh điều kiện đặc trưng của xã Phú Xuyên, chưa mở rộng đối sánh trên toàn dải địa hình sườn Đông và sườn Tây dãy núi Tam Đảo.

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

  • Xây dựng bản đồ phân vùng lập địa thích hợp cho cây Xoan đào trên nền tảng GIS và viễn thám độ phân giải cao Sentinel-2.
  • Ứng dụng công nghệ nuôi cấy mô tế bào và giâm hom cành non có xử lý thuốc kích thích sinh trưởng (IBA, NAA nồng độ 500–1.000 ppm) để chủ động nguồn cây giống chất lượng cao, khắc phục nhược điểm mất sức nảy mầm nhanh của hạt.

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

  • Sinh viên và Học viên Cao học Lâm nghiệp: Tiếp cận hệ thống phương pháp luận hoàn chỉnh về điều tra sinh thái rừng nhiệt đới, công thức tính toán chỉ số quan trọng ($IV%$), chỉ số đa dạng Shannon-Wiener ($H'$), Simpson ($D$) và kỹ thuật phân tích phẫu diện đất chuẩn FIPI.
  • Kỹ sư Lâm sinh và Trạm Kiểm lâm: Nhận được bộ thông số kỹ thuật chuẩn hóa về thời vụ vật hậu (thu hái hạt tháng 9–10) và quy trình ủ cát ẩm duy trì sức sống hạt giống.
  • Cộng đồng dân cư và Chủ rừng tại địa phương: Được hướng dẫn chuyển đổi từ khai thác kiệt sang mô hình khoanh nuôi làm giàu rừng bằng loài gỗ bản địa có giá trị kinh tế cao, cải thiện sinh kế bền vững.
  • Các nhà nghiên cứu sinh thái: Bổ sung dẫn liệu khoa học về loài Pygeum arboreum Endl. vào cơ sở dữ liệu thực vật rừng Việt Nam.

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

1. Cây Xoan đào phân bố tối ưu ở điều kiện lập địa và đai cao nào?

Xoan đào phân bố tự nhiên tối ưu ở đai độ cao dưới 600m, tại các vị trí chân núi và sườn núi ẩm, ven khe suối có độ dốc trung bình. Cây phát triển tốt trên đất feralit phát triển trên đá biến chất, tầng đất dày $>50\text{ cm}$, thành phần cơ giới thịt nhẹ đến thịt trung bình, thoát nước tốt, độ tàn che lâm phần từ 0,6 đến 0,8.

2. Tại sao tỷ lệ nảy mầm của hạt Xoan đào lại suy giảm rất nhanh sau khi thu hái?

Quả và hạt Xoan đào chứa hàm lượng dầu béo cao (40–45%). Khi hạt rơi xuống đất hoặc phơi khô thông thường, lượng dầu bị oxy hóa và phôi hạt mất nước nhanh chóng, làm giảm 50–60% sức nảy mầm chỉ sau 1,5 tháng. Hạt bắt buộc phải được xử lý và ủ ngay trong cát ẩm 80%.

3. Làm thế nào để phân biệt Xoan đào (Pygeum arboreum) với cây Xoan ta (Melia azedarach)?

Xoan đào (Pygeum arboreum) thuộc họ Hoa hồng (Rosaceae), là cây gỗ lớn cao 20–25m, lá đơn nguyên mép nguyên, toàn thân có mùi hôi bọ xít đặc trưng, hoa chùm màu trắng nhiều thùy. Trong khi đó, Xoan ta (Melia azedarach) thuộc họ Xoan (Meliaceae), có lá kép lông chim 2–3 lần, mép lá có răng cưa, hoa màu tím nhạt và không có mùi bọ xít.

4. Kỹ thuật tỉa thưa mở tán trong khoanh nuôi xúc tiến tái sinh Xoan đào được thực hiện như thế nào?

Trong giai đoạn cây mầm ($H < 0,5\text{ m}$), Xoan đào là loài cây chịu bóng nhẹ, cần duy trì thảm che phủ để giữ ẩm đất. Khi cây đạt chiều cao trên $1,0\text{ m}$, cây chuyển sang giai đoạn ưa sáng mạnh, cần tiến hành phát dọn dây leo, chặt bỏ cây bụi tái sinh kém giá trị xung quanh tán, điều chỉnh độ tàn che của tầng cây cao ở mức 0,5–0,6 để cây con vươn tầng vượt tán.

5. Suất đầu tư và thời gian hoàn vốn của mô hình làm giàu rừng bằng cây Xoan đào là bao nhiêu?

Chi phí đầu tư cải tạo và trồng dặm làm giàu rừng bình quân khoảng 20–25 triệu VNĐ/ha trong 3 năm đầu. Chu kỳ kinh doanh gỗ lớn kéo dài 15–20 năm, năng suất gỗ thương phẩm đạt 80–100 $\text{m}^3/\text{ha}$, mang lại doanh thu từ 450–500 triệu VNĐ/ha, thời gian hoàn vốn và sinh lợi bắt đầu từ năm thứ 12 trở đi khi tiến hành tỉa thưa thấu thắp lần cuối.

Kết luận

Công trình nghiên cứu đã giải quyết trọn vẹn các câu hỏi khoa học về đặc điểm lâm học, vật hậu, thổ nhưỡng và động thái tái sinh của cây Xoan đào (Pygeum arboreum Endl) tại xã Phú Xuyên, huyện Đại Từ, tỉnh Thái Nguyên. Bằng việc lượng hóa chính xác chỉ số quan trọng sinh thái ($IV% = 6,66% - 8,06%$), cấu trúc đa dạng loài ($H' = 1,68 - 2,37$) và phát hiện quy luật suy giảm sức sống hạt giống do hàm lượng dầu béo cao, đề tài đã đặt nền móng kỹ thuật vững chắc cho công tác bảo tồn nguồn gen và phát triển kinh tế đồi rừng.

Các cấp quản lý lâm nghiệp, chủ rừng và trạm khuyến nông vùng trung du Bắc Bộ có thể ứng dụng ngay hệ thống giải pháp kỹ thuật gieo ươm ẩm và khoanh nuôi làm giàu rừng theo rạch để khôi phục các diện tích rừng thứ sinh nghèo kiệt, hướng tới mục tiêu quản lý và phát triển rừng bền vững.