Giới thiệu dự án

Bối cảnh và thực trạng ngành lâm nghiệp

Tại Việt Nam, tài nguyên rừng đóng vai trò xung yếu trong việc cân bằng sinh thái, phòng hộ đầu nguồn, chống xói mòn rửa trôi đất, điều hòa khí hậu và bảo tồn đa dạng sinh học. Theo số liệu thống kê lâm nghiệp, mặc dù độ che phủ rừng toàn quốc đã từng bước được phục hồi đạt trên 41%, chất lượng rừng tự nhiên vẫn suy giảm nghiêm trọng. Phần lớn diện tích rừng trồng mới trong nhiều thập kỷ qua là các lâm phần thuần loài cây nhập ngoại mọc nhanh như Keo (Acacia mangium, Acacia auriculiformis), Bạch đàn (Eucalyptus urophylla), hoặc Bồ đề (Styrax tonkinensis). Mô hình thuần loài mọc nhanh bộc lộ nhiều điểm hạn chế mang tính hệ thống: khả năng giữ nước và bảo vệ đất kém, thảm mục nghèo nàn, đất nhanh bị thoái hóa và nguy cơ bùng phát dịch sâu bệnh hại quy mô lớn.

Chính sách phát triển lâm nghiệp bền vững giai đoạn hiện nay đặt trọng tâm vào việc chuyển dịch cơ cấu cây trồng sang các loài cây gỗ lớn bản địa quý hiếm, có giá trị bảo tồn nguồn gen và hiệu quả kinh tế cao. Tuy nhiên, việc gây trồng cây bản địa còn gặp nhiều rào cản do thiếu các nghiên cứu định lượng chuẩn xác về động thái sinh trưởng ban đầu, khả năng thích nghi lập địa và biện pháp kỹ thuật lâm sinh phù hợp.

+-----------------------------------------------------------------------------------+
|               MÔ HÌNH CHUYỂN DỊCH CƠ CẤU CÂY TRỒNG LÂM NGHIỆP                     |
+-----------------------------------------------------------------------------------+
|  LÂM PHẦN THUẦN LOÀI CÂY MỌC NHANH       |   LÂM PHẦN HỖN GIAO CÂY BẢN ĐỊA QUÝ    |
|  - Ưu điểm: Thu hồi vốn nhanh (5-7 năm)  |   - Ưu điểm: Bền vững sinh thái cao    |
|  - Nhược điểm: Đất thoái hóa, sâu bệnh,  |   - Đa dạng cấu trúc tầng tán, giữ đất |
|    giá trị sinh thái và gỗ xẻ thấp       |   - Giá trị gỗ xẻ, tinh dầu, dược liệu |
+-----------------------------------------------------------------------------------+

Vấn đề nghiên cứu và bài toán thực tế

Phần lớn các loài cây gỗ bản địa quý hiếm như Long não, Giổi xanh, Sao đen, Ngọc am và Gù hương đều là những loài có vùng phân bố tự nhiên ngày càng thu hẹp, nhiều loài đứng trước nguy cơ tuyệt chủng và đã được đưa vào Sách Đỏ Việt Nam cũng như Nghị định 32/2006/NĐ-CP (nhóm IIA). Khi đưa các loài này vào mô hình bảo tồn chuyển vị (ex-situ conservation) hoặc trồng rừng cảnh quan tại các vùng sinh thái mới, các nhà lâm học đối mặt với những vấn đề sau:

  • Thiếu dữ liệu sinh trắc học định lượng về tỷ lệ sống và quy luật tăng trưởng đường kính gốc ($D_{00}$), chiều cao vút ngọn ($H_{vn}$) trong giai đoạn cây con.
  • Đất đai vùng đồi núi thấp thường là đất feralit bạc màu, chua, nghèo mùn và dinh dưỡng, đòi hỏi phải xác định được loài cây có biên độ sinh thái phù hợp.
  • Thiếu quy chuẩn mã hóa, theo dõi và đánh giá thống kê sinh học theo chu kỳ thực địa định kỳ.

Mục tiêu của dự án

  1. Định lượng các chỉ tiêu sinh học: Đánh giá tỷ lệ sống ($C%$), sinh trưởng đường kính gốc ($D_{00}$), chiều cao vút ngọn ($H_{vn}$) và động thái ra lá của 5 loài cây bản địa: Long não (Cinnamomum camphora), Giổi xanh (Michelia mediocris), Sao đen (Hopea odorata), Ngọc am (Cupressus funebris), và Gù hương (Cinnamomum balansae).
  2. Đánh giá mức độ thích nghi lập địa: Xác định khả năng thích nghi của từng loài trên nền đất feralit phát triển trên sa thạch tại vườn thực vật Khoa Lâm nghiệp – Đại học Nông lâm Thái Nguyên.
  3. Đề xuất quy trình kỹ thuật lâm sinh: Xây dựng hệ thống giải pháp kỹ thuật từ khâu phân lô, mật độ trồng, làm đất, gán nhãn quản lý cho đến chế độ chăm sóc, phòng trừ sâu bệnh.

Phương pháp tiếp cận và giải pháp

Đề tài áp dụng phương pháp nghiên cứu thực nghiệm lâm sinh kết hợp thống kê sinh học (Biometrics in Forestry). Bố trí thí nghiệm theo phương thức phân lô hỗn giao theo đám với kích thước ô mạng $3\text{m} \times 3\text{m}$, kích thước hố $30\text{cm} \times 30\text{cm} \times 30\text{cm}$, áp dụng hệ thống định danh mã hóa cây hai chữ cái Latinh kết hợp thẻ nhựa chuyên dụng. Số liệu được thu thập 30 ngày/lần liên tục qua 5 đợt và được xử lý bằng các thuật toán thống kê mô tả, xác định sai số tiêu chuẩn ($S$) và hệ số biến động ($S%$).

Kết quả kỳ vọng và chỉ số đo lường

  • Tỷ lệ sống ($C%$): Đạt $\ge 90%$ trên toàn bộ lâm phần thí nghiệm.
  • Tăng trưởng đường kính ($D_{00}$): Tốc độ tăng trưởng bình quân đạt từ $0.04 - 0.06\text{ cm/tháng}$, tổng tăng trưởng sau 5 tháng đạt $\ge 0.20\text{ cm}$.
  • Hệ số biến động ($S%$): Giảm dần qua các chu kỳ đo, thể hiện tính ổn định và sự phân hóa sinh trưởng lành mạnh của quần thụ.

Phạm vi và giới hạn nghiên cứu

  • Đối tượng: 137 cây giống thuộc 5 loài cây bản địa (19 cây Long não, 25 cây Giổi xanh, 45 cây Sao đen, 18 cây Ngọc am, 30 cây Gù hương).
  • Không gian: Vườn thực vật Khoa Lâm nghiệp – Đại học Nông lâm Thái Nguyên (Xã Quyết Thắng, TP. Thái Nguyên; tọa độ địa lý vùng trung du miền núi phía Bắc, đất feralit nâu vàng trên đá sa thạch).
  • Thời gian: Theo dõi liên tục 5 đợt đo định kỳ từ ngày 10/01/2018 đến ngày 10/05/2018 (giai đoạn chuyển tiếp Đông - Xuân - Hè).

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

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

So sánh các giải pháp gây trồng lâm nghiệp hiện nay

Bảng tổng hợp dưới đây phân tích các phương thức trồng rừng hiện hành so với mô hình thực nghiệm vườn thực vật hỗn loài bản địa:

Tiêu chí phân tích Trồng rừng thuần loài mọc nhanh (Keo/Bạch đàn) Trồng bổ sung dưới tán rừng tự nhiên Mô hình vườn thực vật hỗn giao bản địa (Nghiên cứu)
Đa dạng sinh học Rất thấp (đơn loài, phá vỡ chuỗi thức ăn) Trung bình (phụ thuộc thảm rừng mẹ) Cao (tái lập cấu trúc rừng nhiều tầng tán)
Khả năng giữ đất & nước Kém, dễ gây xói mòn khi khai thác trắng Khá tốt Tối ưu nhờ phân tầng rễ và tán lá
Giá trị kinh tế lâu dài Thấp đến trung bình (gỗ dăm, bột giấy) Cao nhưng khó khai thác Rất cao (gỗ quý nhóm I-II, tinh dầu, dược liệu)
Chi phí quản lý ban đầu Thấp Cao do chi phí mở băng, rạch dọn thực bì Trung bình, tối ưu hóa nhờ mã hóa định danh
Nguy cơ dịch bệnh hại Rất cao (dễ bùng phát dịch diện rộng) Thấp Rất thấp nhờ cơ chế ức chế sinh học tự nhiên

Đánh giá điều kiện lập địa và thách thức kỹ thuật

  • Thổ nhưỡng: Đất khu vực nghiên cứu là đất feralit nâu vàng phát triển trên nền đá sa thạch, thành phần cơ giới thịt nhẹ pha cát lẫn đá vụn. Đất có phản ứng chua ($\text{pH} < 4.5$), nghèo mùn và chất hữu cơ, hàm lượng lân dễ tiêu thấp do bị cố định bởi các hợp chất sắt và nhôm.
  • Khí hậu: Mang đặc trưng nhiệt đới gió mùa vùng Đông Bắc Bắc Bộ, mùa mưa tập trung từ tháng 4 đến tháng 10, mùa khô lạnh kéo dài từ tháng 11 đến tháng 3 năm sau. Thời điểm trồng (tháng 10/2017) giúp cây con tận dụng được lượng ẩm cuối mùa mưa để bén rễ trước khi bước vào giai đoạn khô lạnh.

Phân tích yêu cầu kỹ thuật theo khung MoSCoW

+-----------------------------------------------------------------------------------+
|                    PHÂN CẤP YÊU CẦU KỸ THUẬT LÂM SINH (MoSCoW)                    |
+-----------------------------------------------------------------------------------+
|  MUST HAVE (Bắt buộc):                                                            |
|  - Cây con đủ tiêu chuẩn xuất vườn, sạch sâu bệnh, đúng xuất xứ loài.             |
|  - Quy chuẩn hố đào 30x30x30cm, mật độ trồng cự ly 3x3m đồng bộ.                  |
|  - Đánh dấu mốc chuẩn cố định ở cổ rễ bằng sơn/bút xóa chống trôi.                |
|  - Hệ thống mã hóa định danh từng cá thể cây qua thẻ nhựa dập chìm.              |
+-----------------------------------------------------------------------------------+
|  SHOULD HAVE (Nên có):                                                            |
|  - Chu kỳ đo đạc chính xác 30 ngày/lần vào cùng một ngày cố định trong tháng.     |
|  - Giám sát tình trạng phát sinh sâu bệnh theo chu kỳ 5 ngày/lần.                 |
+-----------------------------------------------------------------------------------+
|  COULD HAVE (Có thể có):                                                          |
|  - Tích hợp đo nhiệt độ và độ ẩm đất cục bộ tại thời điểm đo sinh trắc học.       |
|  - Sử dụng phần mềm trắc lượng lâm nghiệp để tự động hóa tính toán thống kê.     |
+-----------------------------------------------------------------------------------+
|  WON'T HAVE (Chưa thực hiện ở giai đoạn này):                                     |
|  - Phân tích sinh hóa hàm lượng tinh dầu và cấu trúc giải phẫu gỗ sâu.            |
+-----------------------------------------------------------------------------------+

Thiết kế hệ thống thực nghiệm

Kiến trúc không gian mô hình và bố trí thực địa

Mô hình được quy hoạch theo dạng lưới ô vuông ma trận, kết hợp hỗn giao theo đám giữa các loài nhằm tối ưu hóa quan hệ hỗ trợ sinh thái giữa các cây lân cận.

                  SƠ ĐỒ KHÔNG GIAN BỐ TRÍ THÍ NGHIỆM VÀ PHÂN LÔ
    
    [ Hàng 1 ] ---- (Ca1) ---- 3.0m ---- (Mi1) ---- 3.0m ---- (Ho1) ---- ...
                      |                    |                    |
                     3.0m                 3.0m                 3.0m
                      |                    |                    |
    [ Hàng 2 ] ---- (Cu1) ---- 3.0m ---- (Ba1) ---- 3.0m ---- (Ca2) ---- ...
                      |                    |                    |
                     3.0m                 3.0m                 3.0m
                      |                    |                    |
    [ Hàng 3 ] ---- (Mi2) ---- 3.0m ---- (Ho2) ---- 3.0m ---- (Cu2) ---- ...
    
    * Ký hiệu loài:
      Ca: Long não (Cinnamomum camphora)        - 19 cây
      Mi: Giổi xanh (Michelia mediocris)        - 25 cây
      Ho: Sao đen (Hopea odorata)               - 45 cây
      Cu: Ngọc am (Cupressus funebris)          - 18 cây
      Ba: Gù hương (Cinnamomum balansae)        - 30 cây
      Tổng dung lượng mẫu theo dõi: N = 137 cây

Công cụ đo lường và thiết bị hỗ trợ

  • Thước kẹp cơ khí Vernier Caliper (Độ chính xác $\pm 0.02\text{mm}$): Sử dụng đo đường kính sát gốc ($D_{00}$) theo hai hướng trực giao vuông góc: Đông – Tây và Nam – Bắc, lấy giá trị trung bình cộng để triệt tiêu sai số hình học của thân cây.
  • Thước dây thép bọc nhựa chuyên dụng (Độ chính xác $\pm 1\text{mm}$): Đo chiều cao vút ngọn ($H_{vn}$) từ điểm mốc cố định tại cổ rễ lên đến mầm sinh trưởng đỉnh ngọn.
  • Hệ thống nhãn định danh Composite Plastic: Bền với thời tiết nhiệt đới, gắn bằng dây rút đàn hồi tại vị trí phân cành cấp 1, không làm tổn thương tầng phát sinh vỏ cây.

Phương pháp luận nghiên cứu (Methodology)

Nghiên cứu áp dụng quy trình 4 giai đoạn chuẩn mực trong điều tra quy luật sinh trưởng rừng trồng:

+------------------+     +------------------+     +------------------+     +------------------+
|  GIAI ĐOẠN 1     |     |  GIAI ĐOẠN 2     |     |  GIAI ĐOẠN 3     |     |  GIAI ĐOẠN 4     |
|  Chuẩn bị &      | --> |  Định vị mốc &   | --> |  Thu thập dữ     | --> |  Xử lý thống kê  |
|  Phân lô thiết kế|     |  Gán mã cá thể   |     |  liệu định kỳ    |     |  & Đánh giá toán |
|  - Cự ly 3x3m    |     |  - Sơn mốc D00   |     |  - 30 ngày/lần   |     |  - Excel / R     |
|  - Đào hố 30cm3  |     |  - Mã hóa 2 chữ  |     |  - 5 đợt đo      |     |  - S, S%, Mean   |
+------------------+     +------------------+     +------------------+     +------------------+
  1. Giai đoạn chuẩn bị (Tháng 10/2017): Dọn thực bì cục bộ, cuốc xới đất, đào hố quy cách $30\text{cm} \times 30\text{cm} \times 30\text{cm}$, bón lót phân vi sinh hữu cơ, trồng cây con đảm bảo bầu không bị vỡ.
  2. Giai đoạn chuẩn hóa mốc đo (Tháng 12/2017): Định vị mốc đo $D_{00}$ tại vị trí cách mặt đất $1.0\text{cm}$ bằng bút kẻ vạch không phai. Gắn thẻ mã hóa (Ca1–Ca19, Mi1–Mi25, Ho1–Ho45, Cu1–Cu18, Ba1–Ba30).
  3. Giai đoạn thu thập số liệu thực địa (01/2018 – 05/2018): Tiến hành đo đạc chính xác vào ngày 10 hàng tháng qua 5 chu kỳ (Đợt 1: 10/01; Đợt 2: 10/02; Đợt 3: 10/03; Đợt 4: 10/04; Đợt 5: 10/05). Giám sát sâu bệnh 5 ngày/lần.
  4. Giai đoạn xử lý sinh trắc học: Làm sạch dữ liệu, loại bỏ ngoại lai (outliers), tính toán các tham số thống kê mô tả mẫu bằng công thức toán lâm sinh.

Implementation và kết quả

Quy trình tính toán sinh trắc học lâm nghiệp

Các thuật toán và công thức toán học được ứng dụng để xử lý số liệu thô điều tra từ 137 cây mẫu:

1. Công thức tính tỷ lệ sống ($C%$)

$$C% = \frac{n}{N} \times 100$$ Trong đó: $n$ là số cây sống thực tế tại thời điểm kiểm tra; $N$ là tổng số cây trồng ban đầu.

2. Đường kính gốc trung bình ($\overline{D_{00}}$) và Chiều cao vút ngọn trung bình ($\overline{H_{vn}}$)

$$\overline{D_{00}} = \frac{\sum_{i=1}^{M} d_i}{M} \quad ; \quad \overline{H_{vn}} = \frac{\sum_{i=1}^{M} h_i}{M}$$ Trong đó: $d_i, h_i$ lần lượt là số đo đường kính gốc và chiều cao của cây thứ $i$; $M$ là tổng số cây đo được.

3. Sai tiêu chuẩn mẫu ($S$)

$$S = \sqrt{\frac{\sum_{i=1}^{n} (x_i - \bar{x})^2}{n - 1}}$$ Trong đó: $x_i$ là giá trị quan sát của cá thể thứ $i$; $\bar{x}$ là giá trị trung bình mẫu; $n$ là dung lượng mẫu.

4. Hệ số biến động sinh trắc ($S%$)

$$S% = \frac{S}{\bar{x}} \times 100$$ Hệ số biến động phản ánh mức độ đồng đều của quần thể cây trồng. $S% < 15%$ biểu thị tính đồng đều cao; $15% \le S% \le 30%$ biểu thị biến động trung bình; $S% > 30%$ biểu thị mức độ phân hóa lớn.

import numpy as np
import pandas as pd

def calculate_forestry_metrics(data_measurements: list) -> dict:
    """
    Tính toán các tham số sinh trắc học lâm nghiệp theo chuẩn điều tra rừng.
    """
    arr = np.array(data_measurements, dtype=np.float64)
    n = len(arr)
    if n <= 1:
        raise ValueError("Dung lượng mẫu n phải lớn hơn 1 để tính sai tiêu chuẩn.")
    
    mean_val = np.mean(arr)
    # ddof=1 cho mẫu ngẫu nhiên (n - 1)
    std_dev = np.std(arr, ddof=1)
    cv_percent = (std_dev / mean_val) * 100.0
    
    return {
        "Sample_Size_N": n,
        "Mean": round(float(mean_val), 4),
        "Standard_Deviation_S": round(float(std_dev), 4),
        "Coefficient_Of_Variation_S_pct": round(float(cv_percent), 3)
    }

# Dữ liệu thực nghiệm đường kính D00 (cm) của Sao Đen tại Lần 5
sao_den_d00_l5 = [0.65, 0.62, 0.58, 0.71, 0.63, 0.60, 0.66, 0.59, 0.64, 0.63] # Trích đoạn 10 cá thể mẫu
results = calculate_forestry_metrics(sao_den_d00_l5)
print("Kết quả tính toán sinh trắc học mẫu:", results)

Kết quả kiểm định thực nghiệm và số liệu đạt được

1. Đánh giá tỷ lệ sống ($C%$)

Sau 5 tháng theo dõi từ mùa đông sang đầu mùa hè, kết quả tỷ lệ sống của cả 5 loài cây bản địa đều đạt mức tuyệt đối:

STT Loài cây nghiên cứu Tên khoa học Tổng số cây ban đầu Số cây sống thực tế Số cây chết Tỷ lệ sống ($C%$)
1 Long não Cinnamomum camphora 19 19 0 100.0%
2 Giổi xanh Michelia mediocris 25 25 0 100.0%
3 Sao đen Hopea odorata 45 45 0 100.0%
4 Ngọc am Cupressus funebris 18 18 0 100.0%
5 Gù hương Cinnamomum balansae 30 30 0 100.0%
Tổng Toàn bộ lâm phần -- 137 137 0 100.0%

2. Động thái tăng trưởng đường kính sát gốc ($D_{00}$) và hệ số biến động ($S%$)

Dữ liệu sinh trưởng đường kính gốc và sự biến thiên của hệ số biến động qua 5 đợt đo được tổng hợp chi tiết tại bảng dưới đây:

STT Loài cây Chỉ tiêu sinh học Lần đo 1 (10/01) Lần đo 2 (10/02) Lần đo 3 (10/03) Lần đo 4 (10/04) Lần đo 5 (10/05) Tổng tăng trưởng ($\Delta D_{00}$)
\multirow{2}{*}{1} \multirow{2}{*}{Long não} $\overline{D_{00}}\text{ (cm)}$ 0.205 0.255 0.305 0.363 0.434 \multirow{2}{*}{+0.229 cm}
$S%\text{ (%)}$ 17.964% 14.445% 12.079% 11.093% 11.549%
\multirow{2}{*}{2} \multirow{2}{*}{Giổi xanh} $\overline{D_{00}}\text{ (cm)}$ 0.114 0.164 0.214 0.306 0.360 \multirow{2}{*}{+0.246 cm}
$S%\text{ (%)}$ 32.329% 22.473% 17.222% 9.804% 10.608%
\multirow{2}{*}{3} \multirow{2}{*}{Sao đen} $\overline{D_{00}}\text{ (cm)}$ 0.382 0.426 0.481 0.551 0.632 \multirow{2}{*}{+0.250 cm}
$S%\text{ (%)}$ 32.062% 27.739% 25.581% 22.927% 19.890%
\multirow{2}{*}{4} \multirow{2}{*}{Ngọc am} $\overline{D_{00}}\text{ (cm)}$ 0.122 0.175 0.225 0.281 0.339 \multirow{2}{*}{+0.217 cm}
$S%\text{ (%)}$ 35.001% 28.148% 21.893% 20.402% 18.632%
\multirow{2}{*}{5} \multirow{2}{*}{Gù hương} $\overline{D_{00}}\text{ (cm)}$ 0.127 0.182 0.248 0.310 0.377 \multirow{2}{*}{+0.250 cm}
$S%\text{ (%)}$ 36.991% 32.712% 32.801% 25.555% 21.414%

3. Phân tích quy luật sinh trưởng

  • Tăng tốc sinh trưởng từ tháng 3 trở đi: Trong 2 tháng đầu (tháng 1 và tháng 2), tốc độ tăng trưởng đường kính bình quân của các loài đạt khoảng $0.05\text{ cm/tháng}$ do điều kiện thời tiết khô lạnh mùa đông. Từ tháng 3 sang tháng 5, khi nhiệt độ tăng và xuất hiện mưa rào mùa xuân, tốc độ tăng trưởng tăng vọt lên $0.07 - 0.10\text{ cm/tháng}$.
  • Xu hướng thu hẹp của hệ số biến động ($S%$): Điểm đổi mới quan trọng là sự sụt giảm mạnh của hệ số biến động ở tất cả 5 loài. Cụ thể, Giổi xanh giảm từ $32.33%$ xuống $10.61%$, Ngọc am giảm từ $35.00%$ xuống $18.63%$, Long não giảm từ $17.96%$ xuống $11.55%$. Điều này chứng minh quần thể cây trồng vượt qua giai đoạn bén rễ sốc môi trường và đang phát triển đồng đều cao.
  • Tình hình sâu bệnh hại: Theo dõi định kỳ 5 ngày/lần trong suốt quá trình thí nghiệm không ghi nhận bất kỳ đợt bùng phát sâu bệnh hại nguy hiểm nào ảnh hưởng đến sinh trưởng ngọn và tán lá của 5 loài cây bản địa.

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

Các điểm cải tiến kỹ thuật nổi bật

  1. Chuẩn hóa hệ thống định danh cá thể bằng mã kép: Ứng dụng phương pháp gán nhãn khoa học dựa trên 2 chữ cái đầu của tên chi (Genus) kết hợp số thứ tự cá thể (ví dụ: Ca cho Cinnamomum, Mi cho Michelia, Ho cho Hopea, Cu cho Cupressus, Ba cho Cinnamomum balansae). Phương pháp này loại bỏ nhầm lẫn thực địa, hỗ trợ truy xuất phả hệ cây rừng lâu dài.
  2. Kỹ thuật cố định mốc đo cổ rễ ($D_{00}$): Sử dụng vạch sơn chuyên dụng chống trôi ở vị trí sát mặt đất giúp loại bỏ sai số vị trí đo kẹp thước giữa các lần đo của những người điều tra khác nhau.
  3. Mô hình hóa hỗn giao theo đám cự ly $3\text{m} \times 3\text{m}$: Tạo khoảng trống dinh dưỡng và chiếu sáng tối ưu, giúp cây con phát huy tính chịu bóng nhẹ ở giai đoạn đầu mà không bị cạnh tranh dinh dưỡng gay gắt từ tầng cây phù trợ.

So sánh với các công trình nghiên cứu trước đây

Tham số so sánh Nghiên cứu của Maurand (1931) tại Trảng Bom Nghiên cứu của Phạm Xuân Hoàn (2002) tại Cát Bà Nghiên cứu của Lê Minh Cường (2007) tại Đại Lải Nghiên cứu hiện tại (Khoa Lâm nghiệp ĐHNL, 2018)
Phương thức trồng Rạch hẹp $1 - 2\text{m}$ dưới tán rừng thứ sinh Dưới tán Keo tai tượng & Keo lá tràm thuần loài Dưới tán rừng Thông mã vĩ thuần loài Vườn thực vật phân lô hỗn giao theo đám ($3\text{m} \times 3\text{m}$)
Tỷ lệ sống ($C%$) Thấp do rạch hẹp, bị tầng cây cao lấn át $79.1%$ (dưới tán Keo tai tượng); Sấu chết 100% $93.2%$ (sau khi tỉa thưa tán thông) 100.0% trên toàn bộ 137 cây mẫu
Tăng trưởng đường kính ($D_{00}$) Sinh trưởng chậm trong 2 năm đầu Tăng trưởng bình quân $0.051\text{ cm/tháng}$ Tăng trưởng bình quân $0.042\text{ cm/tháng}$ Tăng trưởng bình quân $0.043 - 0.050\text{ cm/tháng}$
Tính đồng đều ($S%$) Phân hóa mạnh, không kiểm soát được Chênh lệch lớn giữa các tiểu vùng lập địa Phân hóa theo độ tàn che tầng cây cao Hệ số biến động giảm sâu (ổn định ở mức $10.6 - 21.4%$)

Đóng góp cho khoa học và thực tiễn ngành lâm nghiệp

  • Cơ sở khoa học cho bảo tồn ngoại vi: Cung cấp bằng chứng thực nghiệm vững chắc chứng minh loài cây đặc hữu có nguy cơ tuyệt chủng cao như Gù hương (Cinnamomum balansae) và Ngọc am (Cupressus funebris) hoàn toàn có khả năng sinh trưởng khỏe mạnh trên nền đất feralit đồi núi thấp Thái Nguyên khi được chăm sóc đúng quy chuẩn.
  • Dữ liệu sinh học chuẩn hóa: Bổ sung bộ số liệu thực địa định lượng về 5 loài cây bản địa quý vào cơ sở dữ liệu nghiên cứu lâm học của các trường đại học và viện nghiên cứu chuyên ngành.

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

Kịch bản ứng dụng trong thực tiễn

                   CÁC HƯỚNG ỨNG DỤNG THỰC TẾ CỦA KẾT QUẢ NGHIÊN CỨU
    
    [ 1. CẢI TẠO RỪNG NGHÈO KIỆT ] ===> Thay thế dần keo, bạch đàn bằng Giổi xanh, Sao đen
    [ 2. QUY HOẠCH CẢNH QUAN ĐÔ THỊ ] => Trồng Long não thanh lọc không khí, hấp thụ kim loại nặng
    [ 3. PHÁT TRIỂN LÂM SẢN NGOÀI GỖ ] => Khai thác hạt Giổi (1.5 triệu/kg), tinh dầu Gù hương, Camphor
    [ 4. BẢO TỒN NGUỒN GEN TẠI CHỖ ] ==> Xây dựng hệ thống vườn thực vật tại các Vườn Quốc gia
  1. Chuyển hóa rừng gỗ nhỏ sang rừng gỗ lớn bản địa: Ứng dụng trồng làm giàu rừng dưới tán rừng Keo, Bạch đàn thoái hóa hoặc trồng xen dặm tạo lâm phần hỗn giao đa tầng tán.
  2. Quy hoạch cảnh quan đô thị và trường học: Sử dụng cây Long não và Sao đen làm cây bóng mát đường phố, khuôn viên công sở nhờ khả năng thanh lọc không khí, hấp thu ion chì và phát tán hợp chất hữu cơ dễ bay hơi (Phytoncides) kháng khuẩn.
  3. Phát triển kinh tế nông lâm kết hợp: Trồng Giổi xanh lấy hạt làm gia vị và dược liệu giá trị cao (giá trị thương phẩm đạt $\approx 1.500.000\text{ VNĐ/kg}$ hạt khô) kết hợp khai thác gỗ xẻ sau 20 năm chu kỳ kinh doanh.

Phân tích hiệu quả kinh tế và lợi ích môi trường (Cost-Benefit Analysis)

Hạch toán chi phí và doanh thu dự kiến cho 1.0 ha mô hình trồng Giổi xanh & Gù hương (Mật độ 1.100 cây/ha)

  • Chi phí đầu tư ban đầu (Năm 1 - 3):
    • Cây giống chuẩn F1 (1.100 cây $\times 25.000\text{ VNĐ/cây}$): $27.500.000\text{ VNĐ}$.
    • Nhân công dọn thực bì, đào hố $30\text{cm}^3$, bón phân lót và trồng: $35.000.000\text{ VNĐ}$.
    • Chi phí chăm sóc, làm cỏ, xới xáo (3 năm đầu): $45.000.000\text{ VNĐ}$.
    • Tổng chi phí giai đoạn thiết lập: $\approx 107.500.000\text{ VNĐ/ha}$.
  • Doanh thu dự kiến:
    • Từ năm thứ 7 đến năm 15: Thu hoạch hạt Giổi xanh bình quân $5 - 10\text{ kg hạt tươi/cây/năm}$ đối với $500$ cây Giổi $\rightarrow$ Doanh thu ước đạt $150.000.000 - 300.000.000\text{ VNĐ/ha/năm}$.
    • Sau 20 năm: Khai thác chọn gỗ Gù hương và Giổi xanh (sản lượng ước đạt $0.8 - 1.0\text{ m}^3\text{/cây}$ với giá thị trường $20.000.000 - 25.000.000\text{ VNĐ/m}^3$).
    • Hệ số hoàn vốn nội bộ (IRR) ước tính đạt $> 22.5%$, vượt trội so với chu kỳ 3 luân kỳ trồng keo thuần loài.

Lộ trình triển khai nhân rộng (Roadmap)

+-------------------+     +-------------------+     +-------------------+
|  NĂM 1 - NĂM 2    |     |  NĂM 3 - NĂM 5    |     |  NĂM 6 - NĂM 20   |
|  - Nhân giống &   | --> |  - Tỉa thưa mở    | --> |  - Thu hoạch hạt  |
|    hoàn thiện mô  |     |    tán, điều      |     |    lâm sản phụ    |
|    hình vườn ươm  |     |    chỉnh tầng che |     |  - Khai thác chọn |
|  - Trồng 10 ha    |     |  - Mở rộng 50 ha  |     |    gỗ lớn bền vững|
+-------------------+     +-------------------+     +-------------------+

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

Các hạn chế kỹ thuật của đề tài

  1. Khung thời gian theo dõi ngắn: Chuỗi dữ liệu 5 tháng mới chỉ phản ánh được động thái sinh trưởng trong giai đoạn cây con mới trồng, chưa bao quát được toàn bộ chu kỳ mùa vụ qua các năm biến động thời tiết cực đoan.
  2. Hạn chế trong chỉ tiêu đếm số lá đối với loài lá kim: Đối với loài Ngọc am (Cupressus funebris), do đặc điểm hình thái lá vảy áp sát cành nên không thể áp dụng phương pháp đếm số lượng lá trực tiếp như các loài lá rộng (Long não, Giổi, Sao đen, Gù hương).
  3. Quy mô diện tích thực nghiệm: Thử nghiệm mới dừng lại ở quy mô vườn thực vật thuộc trường đại học, chưa triển khai đa dạng trên các dạng lập địa đồi dốc cao hoặc đất sau nương rẫy thoái hóa nặng.

Định hướng nghiên cứu mở rộng

  • Nghiên cứu sinh sản vô tính: Khảo nghiệm phương pháp nhân giống sinh dưỡng bằng giâm hom kết hợp xử lý chất kích thích sinh trưởng (IBA, NAA nồng độ $500 - 1500\text{ ppm}$) cho loài Gù hương nhằm giải quyết triệt để vấn đề tái sinh tự nhiên kém từ hạt.
  • Theo dõi chu kỳ dài hạn (Long-term ecological monitoring): Tiếp tục quan trắc các chỉ tiêu sinh học ở các mốc tuổi 3, 5, 8 và 10 năm để xây dựng đường cong tăng trưởng chuẩn ($Growth\ Curves$) cho từng loài bản địa.
  • Đánh giá tương tác hóa sinh học (Allelopathy): Nghiên cứu ảnh hưởng ức chế hoặc kích thích tương hỗ giữa rễ cây Long não, Gù hương đối với vi sinh vật đất và các loài cây bụi liên kết.

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

+-----------------------------------------------------------------------------------+
|                    BẢNG ĐỐI TƯỢNG HƯỞNG LỢI VÀ GIÁ TRỊ MANG LẠI                   |
+-----------------------------------------------------------------------------------+
|  1. SINH VIÊN & HỌC VIÊN CAO HỌC NGHÀNH LÂM NGHIỆP:                               |
|     - Tiếp cận mô hình thực địa sống ngay trong khuôn viên trường đại học.        |
|     - Nắm vững phương pháp luận đo đạc, mã hóa và xử lý sinh trắc học chuẩn xác.  |
+-----------------------------------------------------------------------------------+
|  2. CÁC KỸ SƯ LÂM SINH & TRẠM KHUYẾN NÔNG ĐỊA PHƯƠNG:                             |
|     - Sở hữu tài liệu hướng dẫn kỹ thuật gán nhãn và phân lô hỗn giao 3x3m.        |
|     - Căn cứ khoa học để tư vấn cơ cấu cây trồng rừng phòng hộ và đặc dụng.      |
+-----------------------------------------------------------------------------------+
|  3. CÁC DOANH NGHIỆP & HỘ GIA ĐÌNH TRỒNG RỪNG KINH TẾ:                            |
|     - Giảm thiểu rủi ro đầu tư nhờ số liệu tỷ lệ sống 100% của 5 loài cây bản địa.|
|     - Nắm bắt bài toán kinh tế giá trị cao từ Giổi xanh, Gù hương và Long não.    |
+-----------------------------------------------------------------------------------+
|  4. CƠ QUAN QUẢN LÝ BẢO TỒN THIÊN NHIÊN (Vườn Quốc gia / Khu bảo tồn):           |
|     - Phương pháp bảo tồn ngoại vi (Ex-situ) khả thi cho các loài nguy cấp.      |
|     - Dữ liệu phục vụ việc tái thả và làm giàu rừng tự nhiên bị suy thoái.        |
+-----------------------------------------------------------------------------------+

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

1. Yêu cầu kỹ thuật tối thiểu về cây giống khi xuất vườn để đạt tỷ lệ sống 100% như nghiên cứu là gì?

Cây con giống phải đạt từ $6 - 12$ tháng tuổi, chiều cao tối thiểu $30 - 45\text{ cm}$, đường kính cổ rễ $\ge 0.3\text{ cm}$, thân cây hóa gỗ hoàn toàn, đỉnh sinh trưởng không bị cụt hay sâu bệnh, bộ rễ phát triển cân đối và bầu đất nguyên vẹn (kích thước túi bầu tối thiểu $12\text{cm} \times 18\text{cm}$).

2. Mô hình phân lô hỗn giao $3\text{m} \times 3\text{m}$ có bị cạnh tranh dinh dưỡng khi cây bước vào giai đoạn khép tán không?

Cự ly $3\text{m} \times 3\text{m}$ (mật độ $1.111\text{ cây/ha}$) là khoảng cách tối ưu cho $3 - 5$ năm đầu. Khi cây đạt $6 - 8$ năm tuổi và bước vào giai đoạn giao tán mạnh, bắt buộc phải áp dụng biện pháp kỹ thuật tỉa thưa lâm sinh định kỳ, chặt bỏ các cây sinh trưởng kém hoặc tỉa cành nhánh thấp để duy trì mật độ cuối cùng khoảng $400 - 500\text{ cây/ha}$ lấy gỗ lớn.

3. Có thể trồng xen các loài cây này trực tiếp dưới tán rừng Keo 5 - 7 năm tuổi được không?

Hoàn toàn khả thi. Khi trồng dưới tán keo, cần mở băng rạch có chiều rộng từ $2 - 3\text{m}$, tỉa bớt cành nhánh cây keo để đảm bảo độ tàn che tầng trên duy trì ở mức $0.4 - 0.5$. Sau $2 - 3$ năm khi cây bản địa đã bén rễ vững chắc và bắt đầu cần nhiều ánh sáng, tiến hành khai thác tỉa dần cây keo tầng trên.

4. Chi phí duy trì và bảo vệ vườn thực vật trong 3 năm đầu gồm những hạng mục nào?

Chi phí tập trung vào 3 hoạt động chính:

  1. Phát dọn dây leo, cỏ dại xung quanh gốc cây theo đường kính $1.0\text{m}$ ($2 - 3\text{ lần/năm}$).
  2. Xới đất quanh gốc sâu $10\text{cm}$, kết hợp bón thúc $0.1 - 0.2\text{kg}$ phân NPK vi sinh vào đầu mùa mưa.
  3. Kiểm tra định kỳ phòng trừ mối cắn rễ và sâu ăn lá trong mùa xuân hè.

5. Khả năng nhân rộng mô hình ra các tỉnh miền núi phía Bắc khác có gặp rào cản khí hậu không?

Cả 5 loài nghiên cứu (Long não, Giổi, Sao đen, Ngọc am, Gù hương) đều có biên độ sinh thái rộng, thích hợp với khí hậu nhiệt đới gió mùa có mùa đông lạnh vùng đồi núi trung du phía Bắc (như Thái Nguyên, Bắc Kạn, Tuyên Quang, Phú Thọ, Lạng Sơn). Yếu tố quyết định là việc chọn lập địa đất sâu dày, thoát nước tốt và tránh ngập úng cục bộ.


Kết luận

Nghiên cứu sinh trưởng của 5 loài cây bản địa (Long não, Giổi xanh, Sao đen, Ngọc am, Gù hương) tại mô hình vườn thực vật Khoa Lâm nghiệp – Đại học Nông lâm Thái Nguyên đã đạt được các kết quả mang tính khoa học và ứng dụng thực tiễn cao:

  • Tỷ lệ sống tuyệt đối: Đạt $100.0%$ trên toàn bộ 137 cây mẫu sau 5 tháng theo dõi, chứng minh tính thích ứng tốt của các loài cây bản địa quý hiếm trên nền đất feralit đồi núi thấp.
  • Tốc độ tăng trưởng ổn định: Tăng trưởng đường kính gốc ($D_{00}$) bình quân đạt từ $0.217\text{ cm}$ đến $0.250\text{ cm}$, trong đó Sao đen và Gù hương đạt mức tăng trưởng cao nhất ($+0.250\text{ cm}$).
  • Tính đồng đều quần thể vượt trội: Hệ số biến động ($S%$) giảm mạnh qua các chu kỳ đo (Giổi xanh giảm từ $32.33%$ xuống $10.61%$, Long não ổn định ở mức $11.55%$), khẳng định sự thích nghi sinh thái và phân hóa sinh trưởng tích cực.

Kết quả của đề tài là minh chứng khoa học vững chắc cho chiến lược chuyển dịch cơ cấu cây trồng lâm nghiệp từ thuần loài mọc nhanh sang rừng hỗn loài gỗ lớn bản địa, đóng góp thiết thực vào công tác bảo tồn đa dạng sinh học và phát triển kinh tế lâm nghiệp bền vững tại Việt Nam.