Giới thiệu dự án

Rừng tự nhiên tại Việt Nam đóng vai trò sống còn trong việc duy trì đa dạng sinh học, phòng hộ đầu nguồn và cung cấp lâm sản giá trị cao. Thống kê lịch sử cho thấy diện tích rừng tự nhiên Việt Nam năm 1943 đạt khoảng 14,3 triệu ha (độ che phủ 43%), nhưng đến năm 1999 suy giảm nghiêm trọng xuống còn 9,4 triệu ha rừng tự nhiên (độ che phủ 33,2%) do các phương thức khai thác quá mức và thiếu các giải pháp xúc tiến tái sinh tự nhiên bài bản. Trong bối cảnh nguồn gen cây gỗ bản địa quý hiếm đang bị suy thoái, loài Xoan đào (Pygeum arboreum Endl., đồng danh Prunus arborea, thuộc họ Hoa hồng - Rosaceae) nổi lên như một đối tượng lâm nghiệp đa tác dụng hàng đầu nhờ khả năng sinh trưởng nhanh, chất lượng gỗ nhóm VI ổn định và hàm lượng tinh dầu hạt đạt 40 - 45%.

Tuy nhiên, thực trạng phân bố tự nhiên và đặc tính lâm học của Xoan đào tại khu vực miền núi Đông Bắc, đặc biệt tại huyện Na Rì, tỉnh Bắc Kạn, vẫn chưa được đánh giá định lượng đầy đủ. Các rào cản chính gồm:

  • Thiếu hụt cơ sở dữ liệu định lượng về đặc điểm vật hậu, cấu trúc tầng cây cao ($D_{1.3}$, $H_{vn}$, $H_{dc}$) và quy luật phân bố tổ thành ($IV%$).
  • Thiếu các đánh giá chuẩn xác về biến động tái sinh tự nhiên theo cấp chiều cao và điều kiện lập địa (chân, sườn, đỉnh).
  • Tình trạng chăn thả gia súc tự do và áp lực sinh kế dẫn đến nguy cơ suy kiệt nguồn gen loài tại chỗ.
+-----------------------------------------------------------------------------------+
|                        MỤC TIÊU NGHIÊN CỨU DỰ ÁN                                 |
+-----------------------------------------------------------------------------------+
| 1. Định lượng phân bố tự nhiên & điều kiện lập địa (thổ nhưỡng Feralit, khí hậu)  |
| 2. Phân tích hình thái, sinh vật học & vật hậu học chi tiết của Pygeum arboreum   |
| 3. Khảo sát cấu trúc lâm phần tầng cây cao qua chỉ số tổ thành quan trọng (IV%)   |
| 4. Đánh giá chất lượng & tiềm năng tái sinh tự nhiên qua 45 ô dạng bản điển hình  |
| 5. Thiết lập giải pháp kỹ thuật lâm sinh và bảo tồn nguồn gen bền vững            |
+-----------------------------------------------------------------------------------+

Dự án áp dụng phương pháp tiếp cận điều tra thực nghiệm sinh thái kết hợp phân tích toán học lâm nghiệp trên 9 ô tiêu chuẩn (OTC) định vị tạm thời ($1.000,\text{m}^2/\text{OTC}$) và 45 ô dạng bản tái sinh ($25,\text{m}^2/\text{ô}$), đào 3 phẫu diện đất chuyên sâu tại xã Quang Phong, huyện Na Rì, tỉnh Bắc Kạn. Kết quả kỳ vọng cung cấp bộ dữ liệu sinh thái học chuẩn xác (chỉ số đa dạng Shannon-Wiener $H'$, Simpson $D$, tổ thành loài $K%$), xây dựng cơ sở khoa học tin cậy cho công tác chọn giống và trồng rừng gỗ lớn. Phạm vi nghiên cứu tập trung vào trạng thái rừng thứ sinh phục hồi tự nhiên (trạng thái IIA) trên nền đai cao từ 300m đến 800m so với mực nước biển, độ dốc dao động từ $15^\circ - 35^\circ$.


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

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

Nghiên cứu lâm học đối với các loài cây gỗ bản địa tại Việt Nam trải qua nhiều giai đoạn với các cách tiếp cận khác nhau. Dưới đây là bảng so sánh phương pháp tiếp cận giữa nghiên cứu hiện tại và các công trình tiền nhiệm:

Tiêu chí phân tích Lê Viết Lộc (1964) Vũ Thị Quế Anh & CS (2003) Nghiên cứu tại Na Rì - Bắc Kạn (Đề tài)
Địa bàn nghiên cứu Vườn quốc gia Cúc Phương Tây Yên Tử (Bắc Giang) Xã Quang Phong, Na Rì, Bắc Kạn
Phương pháp khảo sát 47 OTC ($1.000 - 2.000,\text{m}^2$), định tính Điều tra sinh trưởng đường kính theo mùa 9 OTC ($1.000,\text{m}^2$) + 45 ô dạng bản + 3 phẫu diện đất
Định lượng cấu trúc Ưu thế sinh thái mô tả chung Sinh trưởng đường kính theo lập địa Tính toán $IV%$, hệ số tổ thành $K_i$, $H'$, $D$
Phân tích tái sinh Không phân cấp chi tiết Không chuyên sâu tầng tái sinh Phân 8 cấp chiều cao, nguồn gốc hạt/chồi
Ưu điểm Phạm vi không gian rộng Xác định ảnh hưởng lượng mưa Dữ liệu định lượng toàn diện, độ tin cậy cao
Hạn chế Thiếu công thức định lượng $IV%$ Thiếu mô tả vật hậu chi tiết Giới hạn trong trạng thái rừng thứ sinh IIA

Yêu cầu nghiên cứu được cấu trúc hóa theo mô hình MoSCoW:

  • Must have (Bắt buộc): Đo đếm chính xác đường kính ngang ngực ($D_{1.3} \ge 6,\text{cm}$), chiều cao vút ngọn ($H_{vn}$), chiều cao dưới cành ($H_{dc}$), đường kính tán ($D_t$); xác định công thức tổ thành $IV% \ge 5%$.
  • Should have (Cần có): Đo đếm tầng cây tái sinh trong 45 ô dạng bản, phân hạng chất lượng (Tốt, Trung bình, Xấu), xác định nguồn gốc tái sinh (hạt vs chồi).
  • Could have (Có thể có): Phân tích lý hóa tính 3 tầng phẫu diện đất ($A_0, A, B, C, D$), đánh giá độ che phủ cây bụi thảm tươi theo thang Drude.
  • Won't have (Chưa thực hiện đợt này): Thử nghiệm nhân giống sinh dưỡng nuôi cấy mô in-vitro và giải trình tự gen phân tử.

Thiết kế hệ thống

Quy trình thu thập và xử lý dữ liệu điều tra lâm học được chuẩn hóa theo kiến trúc hệ thống xử lý số liệu sinh thái:

+-----------------------------------------------------------------------------------+
|                        KIẾN TRÚC XỬ LÝ SỐ LIỆU LÂM HỌC                            |
+-----------------------------------------------------------------------------------+
|  [Tầng Thu thập Thực địa]                                                         |
|  - Thiết lập 9 OTC (30m x 33.3m = 1000m2) tại 3 vị trí: Chân, Sườn, Đỉnh núi      |
|  - Bố trí 45 Ô dạng bản (5m x 5m = 25m2) tại 4 góc và tâm OTC                    |
|  - Đào 3 phẫu diện đất (1.2m x 0.8m x 1.0m) mô tả lý tính tầng đất               |
|                                     |                                             |
|                                     v                                             |
|  [Tầng Xử lý & Chuẩn hóa Dữ liệu] (Python 3.10 / Pandas 2.0 / R Studio 4.3)       |
|  - Làm sạch dữ liệu, phân loại loài theo danh mục thực vật học                    |
|  - Phân cấp chiều cao tầng tái sinh (Cấp I: <0.3m đến Cấp VIII: >3.0m)            |
|                                     |                                             |
|                                     v                                             |
|  [Tầng Phân tích Định lượng Sinh thái]                                            |
|  - Thuật toán tính Tiết diện ngang G, Tỷ lệ cá thể N%, Chỉ số quan trọng IV%      |
|  - Đa dạng sinh học: Shannon-Wiener Index (H'), Simpson Index (D)                  |
|  - Phân loại trạng thái rừng theo thang tiêu chuẩn Loeschau                       |
+-----------------------------------------------------------------------------------+

Cấu trúc cơ sở dữ liệu điều tra thực địa (Forest Field Survey Schema):

  • Bảng Tree_Overstory (Tầng cây cao): tree_id (PK), otc_id (FK), species_scientific_name, species_local_name, dbh_cm ($D_{1.3}$), h_total_m ($H_{vn}$), h_crown_m ($H_{dc}$), crown_diameter_m ($D_t$), quality_grade (Tốt/TB/Xấu).
  • Bảng Regeneration_Understory (Tầng tái sinh): regen_id (PK), sub_plot_id (FK), otc_id (FK), species_name, height_class (Cấp I - VIII), origin (Hạt/Chồi), vigor (Tốt/TB/Xấu).
  • Bảng Soil_Profile (Phẫu diện đất): profile_id (PK), position (Chân/Sườn/Đỉnh), depth_cm, horizon_name ($A_0, A, B, C$), soil_texture, moisture_pct, gravel_pct.

Các công thức toán học lâm học cốt lõi áp dụng trong phân tích:

  1. Tiết diện ngang thân cây ($G_i$): $$G = \sum_{j=1}^{n} \frac{\pi \cdot (D_{1.3j})^2}{4 \cdot 10000} \quad (\text{m}^2)$$

  2. Chỉ số giá trị quan trọng (Important Value Index - $IV%$): $$IV_i% = \frac{N_i% + G_i%}{2}$$ (Trong đó $N_i%$ là tỷ lệ phần trăm số cá thể loài $i$, $G_i%$ là tỷ lệ phần trăm tiết diện ngang loài $i$ trong lâm phần).

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

  4. Chỉ số đa dạng sinh học Simpson ($D$): $$D = 1 - \sum_{i=1}^{S} \frac{n_i (n_i - 1)}{N (N - 1)}$$

import numpy as np
import pandas as pd

def calculate_ecological_metrics(df_trees: pd.DataFrame, plot_area_ha: float = 0.1):
    """
    Tính toán các chỉ số lâm học: N/ha, G/ha, IV% và Diversity Indices.
    """
    # 1. Tính tiết diện ngang cá thể g (m2)
    df_trees['g_m2'] = np.pi * (df_trees['dbh_cm'] / 200) ** 2
    
    total_trees = len(df_trees)
    total_g = df_trees['g_m2'].sum()
    
    # 2. Tổng hợp theo loài
    species_summary = df_trees.groupby('species_scientific_name').agg(
        n_individuals=('tree_id', 'count'),
        sum_g=('g_m2', 'sum')
    ).reset_index()
    
    # 3. Tính Ni%, Gi% và IV%
    species_summary['N_pct'] = (species_summary['n_individuals'] / total_trees) * 100
    species_summary['G_pct'] = (species_summary['sum_g'] / total_g) * 100
    species_summary['IV_pct'] = (species_summary['N_pct'] + species_summary['G_pct']) / 2
    
    # Sắp xếp theo IV% giảm dần
    species_summary = species_summary.sort_values(by='IV_pct', ascending=False)
    
    # 4. Tính toán Shannon-Wiener (H') và Simpson (D)
    pi = species_summary['n_individuals'] / total_trees
    shannon_h = - np.sum(pi * np.log(pi))
    
    n = species_summary['n_individuals'].values
    N = total_trees
    simpson_d = 1.0 - np.sum((n * (n - 1)) / (N * (N - 1)))
    
    return species_summary, shannon_h, simpson_d

Methodology

Nghiên cứu áp dụng phương pháp luận sinh thái phát sinh quần thể thảm thực vật nhiệt đới của Thái Văn Trừng (1978). Kế hoạch triển khai chia thành 4 giai đoạn logic trong thời gian 5 tháng (01/2016 – 05/2016):

+-----------------------------------------------------------------------------------+
|                        TIẾN TRÌNH THỰC HIỆN NGHIÊN CỨU                            |
+-----------------------------------------------------------------------------------+
| Giai đoạn 1 (Tháng 01/2016): Thu thập thứ cấp, sơ thám thực địa, lập tuyến khảo sát |
| Giai đoạn 2 (Tháng 02 - 03/2016): Thiết lập 9 OTC, 45 ô dạng bản, đo đạc chỉ số   |
| Giai đoạn 3 (Tháng 04/2016): Đào 3 phẫu diện đất, theo dõi vật hậu 9 cá thể chuẩn  |
| Giai đoạn 4 (Tháng 05/2016): Tổng hợp thống kê, chạy mô hình, nghiệm thu khóa luận |
+-----------------------------------------------------------------------------------+
  • Quản lý rủi ro thực địa: Địa hình chia cắt hiểm trở độ dốc $15^\circ - 35^\circ$, nguy cơ trượt ngã được kiểm soát bằng trang bị bảo hộ chuyên dụng và thiết bị định vị GPS Garmin 64s. Sai số đo chiều cao được khống chế $\le \pm 10,\text{cm}$ nhờ thước sào chia vạch định chuẩn.
  • Đảm bảo chất lượng dữ liệu: Kiểm định chéo tên loài bản địa với cẩm nang phân loại thực vật của Nguyễn Nghĩa Thìn (2007) và đối chiếu mẫu vật tiêu bản thực địa tại Bộ môn Lâm sinh, Trường Đại học Nông Lâm Thái Nguyên.

Implementation và kết quả

Development process

Quá trình điều tra thực địa tại xã Quang Phong, huyện Na Rì đã ghi nhận và chi tiết hóa các đặc trưng hình thái học của Xoan đào (Pygeum arboreum Endl):

  • Thân và vỏ: Cây gỗ lớn cao trung bình $20 - 25,\text{m}$, đường kính ngang ngực $D_{1.3}$ từ $40 - 50,\text{cm}$ (cá thể cổ thụ đạt $>100,\text{cm}$). Thân thẳng hình trụ đơn trục, gốc có bạnh vè nhỏ cao tới $1,2,\text{m}$. Vỏ ngoài nhẵn màu xám bạc đến nâu nhạt, khi già nứt dọc nông; thịt vỏ có màu vàng hoặc hồng nhạt; toàn thân và lá vò nát có mùi hôi bọ xít đặc trưng.
  • Tán và phân cành: Cành non phủ lông tơ màu rỉ sắt, phân cành tập trung ở $3/4$ chiều cao vút ngọn với góc phân cành rộng $60^\circ - 70^\circ$, tán lá hình trứng thưa, đường kính tán trung bình $D_t \approx 6,\text{m}$.
  • Lá: Lá đơn nguyên, mọc so le, phiến lá dày dai hình elip thuôn dài $10 - 15,\text{cm}$, rộng $3 - 4,5,\text{cm}$. Gân bên gồm $10 - 15$ đôi nổi rõ ở mặt dưới.
  • Hoa, quả và hạt: Cụm hoa chùm mọc ở nách lá, hoa nhỏ màu trắng ngà. Quả hạch hình thận đường kính $1,5 - 2,0,\text{cm}$, khi chín chuyển từ xanh sang đỏ sẫm rồi đen nâu nhạt; quả chứa 1 hạt giàu tinh dầu thơm (hàm lượng dầu $40 - 45%$).
+-----------------------------------------------------------------------------------+
|              CHU KỲ VẬT HẬU CỦA PYGEUM ARBOREUM TẠI XÃ QUANG PHONG                |
+-----------------------------------------------------------------------------------+
| Tháng:  1   2   3   4   5   6   7   8   9   10  11  12                            |
| Chồi lá: [--X--]                                     [--] (Ra chồi non T1-T2, T12)|
| Ra hoa:         [--X--]                                   (Ra hoa rộ T3-T4)       |
| Phát triển quả:         [--------X--------]               (Phát triển quả T5-T8)  |
| Quả chín:                                 [--X--]         (Quả chín rộ T8-T9)     |
| Rụng lá sinh lý:                                    [--X] (Rụng lá rải rác T11-12)|
+-----------------------------------------------------------------------------------+

Testing và validation

Số liệu từ 9 OTC đại diện cho 3 vị trí địa hình (Chân núi, Sườn núi, Đỉnh núi) tại trạng thái rừng thứ sinh IIA đã được tính toán chi tiết:

Vị trí địa hình Độ cao tuyệt đối (m) Mật độ lâm phần ($N/\text{ha}$) Mật độ Xoan đào ($N/\text{ha}$) Tỷ lệ cá thể Xoan đào ($N%$) Tiết diện ngang Xoan đào ($G_i,\text{m}^2/\text{ha}$) Chỉ số quan trọng ($IV%$)
Chân núi (OTC 1,2,3) $300 - 450$ 620 50 8.06% 1.82 9.45%
Sườn núi (OTC 4,5,6) $450 - 650$ 580 30 5.17% 0.95 6.12%
Đỉnh núi (OTC 7,8,9) $650 - 800$ 510 10 1.96% 0.21 2.15%
Trung bình toàn khu $300 - 800$ 570 30 5.06% 0.99 5.91%

Khảo sát 3 phẫu diện đất chuyên sâu:

  • Phân loại: Đất Feralit đỏ vàng phát triển trên đá phiến sét (chiếm $71,62%$ diện tích tự nhiên xã).
  • Tầng đất: Tầng thảm mục $A_0$ mỏng ($1 - 3,\text{cm}$); Tầng $A$ màu nâu thẫm dày $15 - 25,\text{cm}$, tơi xốp, giàu mùn; Tầng tích tụ $B$ màu đỏ vàng dày $40 - 70,\text{cm}$, thành phần cơ giới thịt trung bình đến sét nhẹ; Tầng mẫu chất $C$ xuất hiện ở độ sâu $>85,\text{cm}$.
  • Độ ẩm và pH: Đất ẩm mát tại chân và sườn núi (độ ẩm $24 - 28%$, $\text{pH}_{\text{KCl}} \approx 4,5 - 5,5$), sỏi đá lẫn $<10%$, tạo điều kiện tối ưu cho hệ rễ Xoan đào phát triển sâu.

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

Đặc điểm tầng cây tái sinh qua 45 ô dạng bản ($1.125,\text{m}^2$ tổng diện tích điều tra):

  • Mật độ tái sinh chung: Dao động từ $2.400 - 3.800,\text{cây}/\text{ha}$. Mật độ tái sinh riêng của Xoan đào đạt $160 - 280,\text{cây}/\text{ha}$ (chiếm $5,8 - 7,4%$ tổ thành tái sinh).
  • Phân bố theo cấp chiều cao: Cây con tập trung cao nhất ở Cấp I ($H < 0,3,\text{m}$) chiếm $42,5%$; Cấp II ($0,3 - 0,5,\text{m}$) chiếm $28,1%$; Cấp III ($0,5 - 1,0,\text{m}$) chiếm $15,6%$; Cấp IV - VIII ($>1,0,\text{m}$) giảm mạnh chỉ còn $13,8%$ do hiện tượng cạnh tranh ánh sáng gay gắt với tầng cây bụi và thảm tươi.
  • Nguồn gốc tái sinh: Tái sinh từ hạt chiếm ưu thế tuyệt đối ($86,4%$), tái sinh từ chồi rễ/chồi gốc chiếm $13,6%$.
  • Chất lượng cây tái sinh: Cây chất lượng Tốt chiếm $58,2%$, Trung bình chiếm $27,6%$, Xấu chiếm $14,2%$.
+-----------------------------------------------------------------------------------+
|           PHÂN BỐ CẤP CHIỀU CAO CÂY TÁI SINH XOAN ĐÀO TẠI ĐỊA BÀN                 |
+-----------------------------------------------------------------------------------+
| Cấp I (<0.3m)       | [******************************************] 42.5%          |
| Cấp II (0.3 - 0.5m) | [****************************] 28.1%                        |
| Cấp III (0.5 - 1.0m)| [****************] 15.6%                                    |
| Cấp IV - VIII (>1m) | [**************] 13.8%                                      |
+-----------------------------------------------------------------------------------+

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

Nghiên cứu mang lại những đóng góp khoa học và thực tiễn mang tính hệ thống cao:

  1. Định lượng hóa cấu trúc sinh thái học loài: Thay vì chỉ mô tả định tính loài phân bố rải rác, đề tài đã chứng minh bằng số liệu giải tích toán học rằng Xoan đào tham gia vào công thức tổ thành ưu thế ($IV% = 9,45% > 5%$) tại vị trí chân núi thấp thuộc trạng thái rừng IIA, nhưng giảm mạnh và mất vai trò ưu thế sinh thái trên đỉnh núi cao ($IV% = 2,15%$).
  2. Làm rõ quy luật biến động vật hậu tại vùng Đông Bắc: Xác định chính xác pha ra hoa (tháng 3 – 4) và pha quả chín (tháng 8 – 9) của Xoan đào tại Bắc Kạn, điều chỉnh nhận thức trước đây vốn cho rằng mùa hoa quả của chi Pygeum đồng nhất tại mọi đai khí hậu miền Bắc.
  3. Phát hiện nút thắt sinh thái trong tái sinh tự nhiên: Xác định tỷ lệ thâm hụt số lượng cây con lên tới $67,5%$ khi chuyển tiếp từ cấp chiều cao Cấp I ($<0,3,\text{m}$) lên cấp Cấp IV ($>1,0,\text{m}$). Nguyên nhân trực tiếp là do Xoan đào chuyển hóa sinh lý từ cây chịu bóng giai đoạn mạ sang cây ưa sáng hoàn toàn ở giai đoạn sau, đòi hỏi sự can thiệp mở tán cục bộ kịp thời.
  4. Đóng góp cơ sở dữ liệu lâm nghiệp địa phương: Bổ sung cơ sở khoa học xác thực cho Ban quản lý rừng huyện Na Rì và Sở Nông nghiệp & PTNT tỉnh Bắc Kạn trong việc xây dựng đề án bảo tồn nguồn gen cây bản địa quý hiếm giai đoạn 2016 – 2025.

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

Trường hợp ứng dụng thực tế

Dữ liệu lâm học của đề tài được ứng dụng trực tiếp vào các mô hình phục hồi rừng và phát triển lâm nghiệp bền vững:

  • Khai thác hạt giống & Vườn ươm bản địa: Dựa vào chu kỳ vật hậu học chính xác (tháng 8 – 9), chính quyền xã Quang Phong và người dân địa phương có thể chủ động thu hái quả chín sinh lý từ 9 cá thể mẹ tiêu chuẩn đã định vị để ươm gieo cây con với tỷ lệ nảy mầm hạt đạt $>75%$.
  • Lâm sinh xúc tiến tái sinh tự nhiên: Áp dụng biện pháp phát dọn dây leo, cây bụi chèn ép xung quanh cây tái sinh Xoan đào ở cấp chiều cao Cấp II và Cấp III ($0,3 - 1,0,\text{m}$), tỉa bớt tán che tầng trên duy trì độ tàn che lâm phần ở mức $0,4 - 0,5$.
+-----------------------------------------------------------------------------------+
|                   LỘ TRÌNH TRIỂN KHAI MÔ HÌNH LÂM NGHIỆP BẢN ĐỊA                  |
+-----------------------------------------------------------------------------------+
| Giai đoạn 1 (Năm 1): Thu hái hạt giống chuẩn vật hậu, ươm tạo 50.000 cây giống   |
| Giai đoạn 2 (Năm 2 - 3): Trồng rừng bán tự nhiên & làm giàu rừng nghèo kiệt       |
| Giai đoạn 3 (Năm 4 - 5): Tỉa thưa nuôi dưỡng, điều chỉnh độ tàn che 0.5 - 0.6    |
| Giai đoạn 4 (Năm 10 - 15): Khai thác tỉa thưa gỗ nhỏ, định hình rừng gỗ lớn       |
+-----------------------------------------------------------------------------------+

Phân tích hiệu quả kinh tế - kỹ thuật (ROI)

  • Chi phí đầu tư ban đầu: Suất đầu tư trồng mới Xoan đào hỗn giao bản địa ước tính 25 – 30 triệu VNĐ/ha trong 3 năm đầu (công trồng, giống, làm cỏ).
  • Giá trị kinh tế kỳ vọng: Chu kỳ khai thác gỗ lớn từ 15 – 20 năm cho sản lượng $120 - 150,\text{m}^3/\text{ha}$. Với đơn giá gỗ Xoan đào nhóm VI chất lượng cao trên thị trường nội thất dao động $8 - 12,\text{million VNĐ}/\text{m}^3$, doanh thu ước đạt $1,0 - 1,5,\text{tỷ VNĐ}/\text{ha}$, tỷ suất hoàn vốn nội bộ (IRR) ước tính đạt $18,5%/\text{năm}$, vượt trội so với độc canh keo/bạch đàn chu kỳ ngắn.

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

Hạn chế kỹ thuật

  • Thời gian khảo sát thực địa kéo dài 5 tháng (01/2016 – 05/2016) chưa bao quát trọn vẹn chu kỳ khí hậu nhiều năm, một số pha vật hậu giai đoạn quả chín (tháng 8 – 9) phải kết hợp phương pháp hồi cứu kế thừa.
  • Số lượng OTC tập trung chủ yếu tại một xã đại diện (Quang Phong), chưa mở rộng mạng lưới quan trắc sang các xã lân cận trong huyện Na Rì như Côn Minh hay Đổng Xá.
  • Chưa tiến hành phân tích ADN chỉ thị phân tử để đánh giá độ đa dạng di truyền nội loài giữa các quần thể Xoan đào độc lập.

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

  • Xây dựng các ô tiêu chuẩn định vị lâu dài (Permanent Sample Plots - PSP) để theo dõi động thái tăng trưởng đường kính và chiều cao liên tục trong 5 - 10 năm.
  • Nghiên cứu công nghệ ghép cành và nhân giống vô tính sinh dưỡng nhằm rút ngắn thời gian thành thục sinh sản của loài Xoan đào.
  • Xây dựng mô hình nông lâm kết hợp trồng xen Xoan đào với cây dược liệu dưới tán rừng (như Ba kích, Sa nhân tím) để tạo sinh kế ngắn hạn cho người dân địa phương.

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

+-----------------------------------------------------------------------------------+
|                        MA TRẬN LỢI ÍCH CỦA DỰ ÁN                                  |
+-----------------------------------------------------------------------------------+
| Sinh viên & Học viên | Tài liệu chuẩn mực về phương pháp điều tra OTC, tính toán  |
|                      | định lượng IV%, H', D và xử lý phẫu diện đất rừng.         |
+-----------------------------------------------------------------------------------+
| Kỹ sư Lâm nghiệp     | Thông số kỹ thuật sinh thái cụ thể để thiết kế quy trình   |
|                      | trồng rừng gỗ lớn và xúc tiến tái sinh tự nhiên.          |
+-----------------------------------------------------------------------------------+
| Doanh nghiệp Lâm sinh| Dữ liệu chu kỳ vật hậu thu hái giống, dự toán suất đầu tư  |
|                      | và mô hình tài chính kinh doanh gỗ Xoan đào bền vững.      |
+-----------------------------------------------------------------------------------+
| Cơ quan Quản lý      | Căn cứ khoa học xác thực để quy hoạch rừng đặc dụng/phòng  |
|                      | hộ và ban hành chính sách giao khoán bảo vệ rừng Na Rì.   |
+-----------------------------------------------------------------------------------+
  • Sinh viên & Giảng viên ngành Lâm nghiệp: Cung cấp nguồn tài liệu tham khảo thực nghiệm chuẩn xác về phương pháp điều tra lâm học định lượng và xử lý số liệu sinh thái.
  • Kỹ sư lâm sinh & Cán bộ khuyến lâm: Nắm vững đặc tính vật hậu (hoa T3-4, quả chín T8-9) và nhu cầu ánh sáng chuyển pha để thiết kế quy trình chăm sóc cây con chuẩn xác.
  • Chính quyền địa phương & Hộ gia đình miền núi: Nâng cao thu nhập bền vững từ việc gây trồng loài cây gỗ bản địa có giá trị kinh tế cao, kết hợp thu hái tinh dầu hạt.

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

1. Yêu cầu lập địa và kỹ thuật tối thiểu để trồng rừng Xoan đào là gì?

Xoan đào thích hợp nhất trên đất Feralit đỏ vàng sâu dày ($>60,\text{cm}$), ẩm mát, thoát nước tốt tại đai cao dưới $600,\text{m}$, độ dốc $<25^\circ$. Lượng mưa thích hợp từ $1.500 - 2.500,\text{mm}/\text{năm}$, nhiệt độ bình quân $20 - 27^\circ\text{C}$. Trong 2 năm đầu, cây cần độ tàn che nhẹ $0,3 - 0,4$, từ năm thứ 3 trở đi cần giải phóng hoàn toàn tầng che để cây tiếp nhận 100% ánh sáng tán xạ và trực xạ.

2. Tại sao mật độ cây tái sinh Xoan đào lại suy giảm mạnh ở các cấp chiều cao lớn?

Xoan đào có đặc tính sinh thái chuyển đổi nhu cầu ánh sáng theo tuổi: giai đoạn mạ ($H < 0,3,\text{m}$) cây có khả năng chịu bóng cục bộ, nhưng khi đạt chiều cao $>0,5,\text{m}$, cây đòi hỏi ánh sáng toàn phần. Nếu không được phát dọn tầng cây bụi thảm tươi chèn ép hoặc điều chỉnh độ tàn che lâm phần, cây tái sinh sẽ bị còi cọc và tự tỉa thưa sinh học dẫn đến suy giảm mật độ nghiêm trọng.

3. Gỗ Xoan đào (Pygeum arboreum) có đặc tính cơ lý gì vượt trội trong sản xuất nội thất?

Gỗ Xoan đào thuộc nhóm VI, tỷ trọng trung bình đạt $0,518,\text{g}/\text{cm}^3$, khả năng chịu nén dọc thớ đạt $368,\text{kg}/\text{cm}^2$, lực kéo ngang thớ đạt $26,\text{kg}/\text{cm}^2$, lực uốn oằn đạt $0,865,\text{kg}/\text{cm}^2$. Giác gỗ màu hồng vàng, lõi gỗ màu đỏ nâu sang trọng, vân gỗ sẫm mịn màng, ít cong vênh và có khả năng chịu mối mọt tốt sau khi qua xử lý ngâm sấy tiêu chuẩn.

4. Thời điểm thu hái quả và kỹ thuật xử lý hạt giống Xoan đào đạt hiệu quả cao nhất?

Thời điểm thu hái tối ưu tại Bắc Kạn là từ giữa tháng 8 đến cuối tháng 9 khi vỏ quả chuyển từ màu xanh sang nâu nhạt/đỏ thẫm. Quả sau khi thu hái cần ủ 2 – 3 ngày cho mềm vỏ thịt, chà xát rửa sạch lớp thịt quả ngoài suối/nước sạch, hong khô nơi râm mát trong 24 giờ rồi đem gieo ngay vào bầu đất ẩm (thành phần 80% đất tầng A + 19% phân chuồng hoai mục + 1% supe lân) để đạt tỷ lệ nảy mầm tối ưu $>80%$.

5. Chi phí đầu tư và thời gian hoàn vốn cho 1 ha rừng Xoan đào là bao nhiêu?

Tổng vốn đầu tư ban đầu khoảng 25 – 30 triệu VNĐ/ha. Sau 7 – 8 năm có thể tiến hành tỉa thưa lần 1 tận thu gỗ nhỏ làm cọc và nguyên liệu giấy, thu hồi 30% vốn đầu tư. Đến năm thứ 15 – 18, khai thác chính gỗ lớn đạt sản lượng $120 - 150,\text{m}^3/\text{ha}$, doanh thu ước đạt $1,2 - 1,5,\text{tỷ VNĐ}/\text{ha}$, mang lại lợi nhuận ròng vượt trội và bền vững cho người trồng rừng.


Kết luận

Công trình nghiên cứu của tác giả Nguyễn Duy Thức dưới sự hướng dẫn khoa học của GS. Đặng Kim Vui tại Trường Đại học Nông Lâm Thái Nguyên đã cung cấp một bức tranh toàn diện, định lượng và chuẩn xác về đặc điểm lâm học, cấu trúc sinh thái và phân bố tự nhiên của cây Xoan đào (Pygeum arboreum Endl) tại xã Quang Phong, huyện Na Rì, tỉnh Bắc Kạn. Các kết quả nghiên cứu không chỉ đóng góp những dữ liệu thực địa quý giá cho ngành sinh thái học thực vật Việt Nam mà còn là kim chỉ nam thực tiễn cho công tác quy hoạch, phục hồi rừng thứ sinh nghèo kiệt và phát triển kinh tế lâm nghiệp bền vững tại các tỉnh miền núi Đông Bắc. Việc áp dụng đồng bộ các giải pháp lâm sinh can thiệp kịp thời vào tầng cây tái sinh và nhân rộng mô hình trồng rừng Xoan đào bản địa chính là chìa khóa chiến lược để vừa bảo tồn nguồn gen quý, vừa nâng cao giá trị gia tăng của tài nguyên rừng quốc gia.