Giới thiệu dự án

Rừng nhiệt đới đóng vai trò là "lá phổi xanh" điều hòa khí hậu toàn cầu, lưu trữ trên 45% lượng carbon trên cạn và duy trì hơn 50% tính đa dạng sinh học của thế giới. Tại Việt Nam, theo số liệu từ Viện Điều tra Quy hoạch Rừng, độ che phủ rừng từng giảm mạnh từ 43% (năm 1945) xuống mức đáy 27,2% vào năm 1990 trước khi có sự can thiệp của Chỉ thị 286/TTg. Quá trình khai thác tài nguyên kiệt quệ và tập quán canh tác nương rẫy kéo dài đã đẩy phần lớn diện tích rừng tự nhiên vào các trạng thái suy thoái thứ sinh, phổ biến nhất là trạng thái rừng phục hồi IIA (theo thang phân loại Loeschau, 1966).

+-------------------------------------------------------------------------+
|                  HỆ THỐNG DIỄN THẾ PHỤC HỒI RỪNG THỨ SINH               |
|                                                                         |
|  [Rừng Suy Thoái/Đất Nương Rẫy]                                         |
|  [Trạng thái Rừng Phục hồi IIA] (Mật độ cây gỗ > 1000 cây/ha, D1.3 > 10cm)|
| [Tái sinh tự nhiên]  [Cạnh tranh sinh thái]                             |
| (Hạt: 56.49%)        (Tầng thảm tươi & Cây bụi chèn ép)                 |
+-------------------------------------------------------------------------+

Vấn đề nghiên cứu và Hiện trạng (Problem Statement)

Khu vực xã Yên Lãng, huyện Đại Từ, tỉnh Thái Nguyên có 2.647,45 ha đất lâm nghiệp (chiếm 68,20% tổng diện tích tự nhiên), tiếp giáp Vườn Quốc gia Tam Đảo. Tuy nhiên, hệ sinh thái rừng phục hồi trạng thái IIA tại đây đối mặt với 3 thách thức lớn:

  1. Mất cân bằng cấu trúc sinh thái: Tầng cây cao chủ yếu gồm cây tiên phong ưa sáng giá trị kinh tế kém (Thành ngạnh, Thẩu tấu, Hu đay), thiếu vắng các loài cây gỗ lớn bản địa đạt chuẩn.
  2. Áp lực cạnh tranh sinh dưỡng khốc liệt: Tầng cỏ lào, cây bụi và dây leo sinh trưởng dày đặc (độ nhiều Drude mức Cop1 - Soc), ngăn chặn hạt nảy mầm và gây suy giảm tỷ lệ cây mạ bước sang cấp chiều cao lớn hơn.
  3. Thiếu cơ sở dữ liệu định lượng phục vụ khoanh nuôi: Các can thiệp lâm sinh trước đây chủ yếu mang tính định tính, thiếu các chỉ số đo lường sinh thái chuẩn xác (như chỉ số mức độ quan trọng IVI, đa dạng Shannon-Weaver, tỷ lệ cây triển vọng CTV).

Mục tiêu dự án

  1. Xác định đặc trưng cấu trúc tầng cây gỗ và thảm thực vật tái sinh tự nhiên tại 3 vị trí địa hình (chân, sườn, đỉnh) của trạng thái rừng IIA tại xã Yên Lãng.
  2. Đo lường định lượng các chỉ số sinh thái: Mật độ ($N/ha$), Hệ số tổ thành ($K_i$), Chỉ số quan trọng sinh thái ($IVI$), Đa dạng loài ($H'$), Nguồn gốc tái sinh (Hạt vs Chồi), và Phân bố theo 7 cấp chiều cao ($H_{vn}$).
  3. Thiết lập hệ thống giải pháp kỹ thuật lâm sinh ứng dụng phương thức khoanh nuôi xúc tiến tái sinh tự nhiên có tác động chọn lọc, định hướng phục hồi rừng bền vững.

Phạm vi và Giới hạn

  • Địa điểm: Xã Yên Lãng, huyện Đại Từ, tỉnh Thái Nguyên (tọa độ tiếp giáp dãy núi Tam Đảo và vùng mỏ than Núi Hồng).
  • Phạm vi khảo sát: 3 ô tiêu chuẩn lớn ($OTC = 2.500\text{ m}^2$) và 12 ô thứ cấp ($25\text{ m}^2$), 12 ô dạng bản ($1\text{ m}^2$).
  • Thời gian thực hiện: 01/2015 – 05/2015.

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

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

Phương pháp tiếp cận Ưu điểm Nhược điểm Đánh giá hiệu quả kinh tế/sinh thái
Trồng rừng thuần loài (Keo, Bạch đàn) Tăng độ che phủ nhanh, dễ quản lý chu kỳ ngắn (5-7 năm) Đơn giản hóa cấu trúc hệ sinh thái, xói mòn đất, nghèo kiệt đa dạng sinh học Chi phí đầu tư cao (~15-20 triệu VNĐ/ha), rủi ro sâu bệnh lớn
Khoanh nuôi bảo vệ thụ động (Đóng cửa rừng thuần túy) Chi phí đầu tư thấp Diễn thế kéo dài 30-50 năm, tầng thảm tươi chèn ép khiến cây gỗ mục đích bị đào thải Hiệu quả phục hồi chậm, mật độ cây chất lượng tốt không đồng đều
Khoanh nuôi xúc tiến tái sinh định lượng (Đề xuất) Tận dụng ngân hàng hạt tự nhiên, bảo tồn đa dạng di truyền bản địa, rút ngắn thời gian thành rừng 40% Đòi hỏi điều tra sinh thái và kỹ thuật điều tiết tán chính xác Tiết kiệm 65% ngân sách so với trồng mới, gia tăng giá trị gỗ và phòng hộ lâu dài

Đánh giá yêu cầu theo mô hình MoSCoW

  • Must Have: Đánh giá chính xác mật độ tầng cây gỗ ($D_{1.3} \ge 6\text{ cm}$), cấu trúc mật độ tầng cây tái sinh ($D_{1.3} < 6\text{ cm}$), xác định tỷ lệ cây triển vọng ($CTV \ge 1.0\text{ m}$).
  • Should Have: Tính toán chỉ số mức độ quan trọng $IVI$ của từng loài gỗ cao, phân tích tương quan chỉ số Shannon-Weaver ($H'$) theo 3 vị trí lập địa (Chân - Sườn - Đỉnh).
  • Could Have: Đánh giá chi tiết trữ lượng hạt giống tầng thảm mục ($g/\text{m}^2$) và khảo sát phẫu diện đất $120 \times 80 \times 120\text{ cm}$.
  • Won't Have: Xây dựng mô hình diễn thế vi khí hậu tự động hóa bằng cảm biến IoT (để dành cho pha mở rộng).

Thiết kế hệ thống khảo nghiệm và công nghệ phân tích

+-----------------------------------------------------------------------------------+
|               SƠ ĐỒ BỐ TRÍ KHÔNG GIAN ĐIỀU TRA LÂM HỌC ĐỊNH LƯỢNG                |
|                                                                                   |
|  +-----------------------------------------------------------------------------+  |
|  | Ô Tiêu Chuẩn Sơ Cấp (OTC): 50m x 50m (S = 2.500 m²)                         |  |
|  | Đánh giá tầng cây gỗ D1.3 ≥ 6cm, Hvn, Dtán, Phẫu diện đất                    |  |
|  |                                                                             |  |
|  |   [OTC Thứ cấp 1] (5m x 5m)                     [OTC Thứ cấp 2] (5m x 5m)  |  |
|  |                                                                             |  |
|  |                        [OTC Thứ cấp 5] (5m x 5m)                            |  |
|  |                                                                             |  |
|  |   [OTC Thứ cấp 4] (5m x 5m)                     [OTC Thứ cấp 3] (5m x 5m)  |  |
|  +-----------------------------------------------------------------------------+  |
+-----------------------------------------------------------------------------------+

Technology Stack và Công cụ xử lý

  • Thước đo chuyên dụng: Thước kẹp kính Haglöf Mantax Precision (Đức), Sào đo cao lồng rút 15m, Địa bàn địa chất Harbin DQL-8 xác định phương vị và góc dốc.
  • Phân loại sinh thái: Hệ thống phân loại thảm thực vật Thái Văn Trừng (1978), Loeschau (1966), thang phân cấp độ nhiều Drude.
  • Nền tảng xử lý dữ liệu: Python 3.10 runtime (sử dụng thư viện pandas v2.1.0, scipy v1.11.2, matplotlib v3.8.0) phục vụ mô hình hóa cấu trúc phân bố thống kê.

Implementation và kết quả

Quy trình thu thập và Thuật toán tính toán chỉ số sinh thái

Quy trình chuẩn hóa dữ liệu lâm học bao gồm 4 công đoạn cốt lõi:

Dưới đây là mã nguồn Python triển khai thuật toán tính toán các chỉ số tổ thành sinh thái ($IVI$), chỉ số đa dạng Shannon-Weaver ($H'$), và phân bố cây triển vọng ($CTV$):

import numpy as np
import pandas as pd

def calculate_ecological_metrics(tree_df, total_plot_area_ha=0.25):
    """
    Tính toán chỉ số IVI, Shannon-Weaver H' và tỷ lệ cây triển vọng CTV.
    tree_df columns: ['species', 'd13_cm', 'hvn_m', 'quality', 'origin']
    """
    # 1. Tính tiết diện thân G (m2)
    tree_df['basal_area_m2'] = np.pi * ((tree_df['d13_cm'] / 100) / 2) ** 2
    total_basal_area = tree_df['basal_area_m2'].sum()
    total_individuals = len(tree_df)
    
    # 2. Tổng hợp theo loài
    species_summary = tree_df.groupby('species').agg(
        n_i=('species', 'count'),
        g_i=('basal_area_m2', 'sum')
    ).reset_index()
    
    # 3. Tính độ phong phú tương đối A% và độ ưu thế tương đối D%
    species_summary['A_pct'] = (species_summary['n_i'] / total_individuals) * 100
    species_summary['D_pct'] = (species_summary['g_i'] / total_basal_area) * 100
    
    # 4. Tính Importance Value Index (IVI)
    species_summary['IVI_pct'] = (species_summary['A_pct'] + species_summary['D_pct']) / 2
    
    # 5. Tính chỉ số đa dạng Shannon-Weaver H'
    p_i = species_summary['n_i'] / total_individuals
    shannon_h = -np.sum(p_i * np.log(p_i))
    
    # 6. Tính tỷ lệ cây triển vọng CTV (Cao >= 1.0m, chất lượng Tốt/TB)
    prospective_trees = tree_df[(tree_df['hvn_m'] >= 1.0) & (tree_df['quality'].isin(['Tot', 'TrungBinh']))]
    ctv_rate = (len(prospective_trees) / total_individuals) * 100
    
    density_per_ha = total_individuals / total_plot_area_ha
    
    return {
        "species_composition": species_summary.sort_values(by='IVI_pct', ascending=False),
        "shannon_index": round(shannon_h, 3),
        "density_per_ha": round(density_per_ha, 1),
        "ctv_rate_pct": round(ctv_rate, 2)
    }

Kết quả nghiên cứu thực nghiệm

1. Cấu trúc tầng cây gỗ ($D_{1.3} \ge 6\text{ cm}$)

  • Mật độ tầng cây gỗ: Biến động từ $308 - 336\text{ cây/ha}$ (trung bình vị trí chân: $324\text{ cây/ha}$, sườn: $328\text{ cây/ha}$, đỉnh: $320\text{ cây/ha}$).
  • Tổ thành sinh thái: Số loài gỗ xuất hiện trong các OTC dao động từ $16 - 24\text{ loài}$, trong đó có $6 - 10\text{ loài}$ có $IVI \ge 5%$. Các loài chiếm ưu thế sinh thái gồm: Thành ngạnh (Cratoxylum prunifolium), Chẹo tía (Engelhardtia chrysolepis), Kháo đá (Machilus bonii), Lim xẹt (Peltophorum dasyrrhachis), Dẻ gai (Castanopsis indica), Xồi (Lithocarpus dealbatus).

2. Cấu trúc thảm thực vật tái sinh tự nhiên ($D_{1.3} < 6\text{ cm}$)

  • Mật độ tái sinh: Trung bình toàn lâm phần đạt $4.133\text{ cây/ha}$.
    • Vị trí Sườn đạt mật độ cao nhất: $4.640 - 5.040\text{ cây/ha}$ (trung bình $4.860\text{ cây/ha}$).
    • Vị trí Chân: $3.680 - 4.160\text{ cây/ha}$ (trung bình $3.860\text{ cây/ha}$).
    • Vị trí Đỉnh: $3.520 - 3.840\text{ cây/ha}$ (trung bình $3.680\text{ cây/ha}$).
  • Chỉ số đa dạng sinh học Shannon-Weaver ($H'$): Biến động từ $2.64 - 2.81$, khẳng định tính đa dạng loài ở mức khá cao và tương đối đồng đều giữa các điều kiện vi địa hình.
Mật độ cây tái sinh theo địa hình (cây/ha):
Chân : ███████████████████▌ 3.860
Sườn : ████████████████████████ 4.860
Đỉnh : ██████████████████ 3.680

3. Chất lượng, nguồn gốc và tỷ lệ cây triển vọng (CTV)

  • Nguồn gốc tái sinh: Tái sinh từ hạt chiếm tỷ lệ vượt trội $56,49% - 91,67%$ (mật độ $2.720 - 3.440\text{ cây/ha}$), tái sinh chồi chiếm $8,33% - 43,51%$.
  • Phân cấp chất lượng: Cây tốt ($29,77% - 45,65%$) và cây trung bình ($29,39% - 42,86%$) chiếm tổng cộng trên $75%$, gấp 2 lần tỷ lệ cây xấu ($13,04% - 28,26%$).
  • Cây triển vọng (CTV): Chiếm từ $32,61% - 41,67%$ tổng số cây tái sinh, đảm bảo nguồn bổ sung hoàn hảo cho tầng tán chính.

4. Quy luật phân bố số cây và số loài theo 7 cấp chiều cao

Cấp chiều cao ($H_{vn}$) Số cây TB (cây/ha) Tỷ lệ số cây (%) Số loài TB (loài/OTC) Tỷ lệ số loài (%) Đặc tính cạnh tranh sinh thái
< 0,5 m 873 21,12% 11 54,08% Giai đoạn mạ, cạnh tranh ánh sáng với thảm cỏ rậm
0,5 – 1,0 m 933 22,57% 12 58,16% Sinh trưởng sung mãn, tích lũy sinh khối
1,0 – 1,5 m 620 15,00% 8 38,78% Vượt qua tầng thảm tươi, bắt đầu phân hóa cá thể
1,5 – 2,0 m 567 13,72% 7 35,20% Cạnh tranh không gian dinh dưỡng tầng nhỡ
2,0 – 2,5 m 480 11,61% 6 30,10% Đào thải tự nhiên mạnh mẽ đối với cây chịu bóng yếu
2,5 – 3,0 m 407 9,84% 5 25,51% Bắt đầu tiếp cận tầng tán dưới của cây gỗ
> 3,0 m 273 6,61% 3 17,35% Cây đã ổn định, chuẩn bị gia nhập tầng tán chính
Tổng cộng 4.133 100,00% 20 - Quy luật phân bố giảm liên tục theo chiều cao

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

Đột phá khoa học và Đóng góp kỹ thuật

  1. Ứng dụng tổ hợp chỉ số định lượng ($IVI$ + $H'$ + Cấp $H_{vn}$): Thay thế hoàn toàn cách tiếp cận mô tả cảm tính trước đây bằng hệ thống ma trận số liệu sinh thái chuẩn mực, cho phép xác định chính xác thời điểm rừng bước vào giai đoạn cạnh tranh khốc liệt (ngưỡng chuyển tiếp $1,0 - 1,5\text{ m}$).
  2. Làm rõ vai trò cấu trúc tái sinh hạt bản địa: Chứng minh tỷ lệ tái sinh hạt tự nhiên đạt trên $56,49%$, với sự xuất hiện thực tế của các loài gỗ có giá trị như Dẻ gai, Chò chỉ, Mỡ, Re, Bứa... Đây là cơ sở thực chứng bác bỏ quan điểm cho rằng rừng IIA suy thoái không còn khả năng tự phục hồi mà bắt buộc phải phá bỏ để trồng cây ngoại lai.
  3. Mô hình hóa tác động của vi địa hình: Xác định sườn dốc tích tụ mật độ cây tái sinh cao nhất ($4.860\text{ cây/ha}$), trong khi đỉnh dốc chịu hạn chế về tầng đất mặt và độ ẩm dẫn đến số loài ở cấp $>3\text{ m}$ tiệm cận 0 loài nếu không có tác động tỉa thưa cục bộ.

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

Quy trình kỹ thuật lâm sinh can thiệp chọn lọc

+-----------------------------------------------------------------------------------+
|               QUY TRÌNH KỸ THUẬT LÂM SINH KHÔI PHỤC RỪNG TRẠNG THÁI IIA           |
|                                                                                   |
|  [Bước 1: Điều tra Định vị Ô Tiêu Chuẩn & Cây Triển Vọng (CTV)]                   |
|  [Bước 2: Xử lý Tầng Thảm tươi & Dây leo Chèn ép (Drude > Cop2)]                  |
|  [Bước 3: Tỉa thưa Vệ sinh Tầng Tán Cây Tiên phong Suy thoái]                    |
|  [Bước 4: Trồng Dặm Bổ sung Cây Gỗ Lớn Bản địa (Mật độ 300-500 cây/ha)]           |
|  [Bước 5: Bảo vệ Nghiêm ngặt & Giám sát Tăng trưởng Định kỳ 12 tháng/lần]        |
+-----------------------------------------------------------------------------------+

Phân tích Hiệu quả Đầu tư và Tài chính (ROI)

  • Suất đầu tư khoanh nuôi xúc tiến tái sinh có trồng dặm: Ước tính khoảng $4.500.000 - 6.000.000\text{ VNĐ/ha/3 năm}$.
  • So sánh với trồng mới rừng phòng hộ/sản xuất: Giảm thiểu $60 - 70%$ chi phí đầu tư ban đầu (trồng mới bình quân $18.000.000 - 22.000.000\text{ VNĐ/ha}$).
  • Giá trị kinh tế lâu dài: Chuyển hóa rừng phục hồi nghèo kiệt thành rừng gỗ lớn sau 15-20 năm, giá trị trữ lượng ước đạt $150 - 200\text{ m}^3\text{/ha}$, tăng nguồn thu bền vững từ lâm sản ngoài gỗ và dịch vụ chi trả môi trường rừng (PES).

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

Hạn chế kỹ thuật

  1. Phạm vi chu kỳ thời gian: Dữ liệu điều tra phản ánh đặc trưng sinh thái trong mùa khô sang đầu mùa mưa (5 tháng), chưa bao quát trọn vẹn biến động nảy mầm theo mùa rụng hạt của tất cả các loài gỗ đỉnh cực.
  2. Số lượng ô phẫu diện đất: Mỗi OTC mới đào 1 phẫu diện đại diện, chưa phân tích đầy đủ các chỉ số hóa - lý chuyên sâu (như hàm lượng mùn tổng số, N-P-K dễ tiêu và vi sinh vật đất).

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

  • Ứng dụng GIS & Viễn thám (Remote Sensing): Xây dựng bản đồ phân vùng mật độ cây tái sinh toàn huyện Đại Từ bằng ảnh vệ tinh Sentinel-2 và máy bay không người lái (UAV).
  • Thiết lập ô định vị sinh thái lâu dài (Permanent Sample Plots - PSP): Theo dõi động thái tăng trưởng đường kính ($D_{1.3}$) và chiều cao ($H_{vn}$) của cây tái sinh triển vọng theo chu kỳ 3 năm và 5 năm.

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

  • Sinh viên & Học viên Lâm nghiệp: Nắm vững phương pháp luận lập ô tiêu chuẩn, phân cấp 7 cấp chiều cao và thuật toán tính toán chỉ số $IVI$, Shannon-Weaver.
  • Kỹ sư Lâm sinh & Cán bộ Kiểm lâm: Sở hữu sổ tay hướng dẫn kỹ thuật điều tiết tán che ($0,6 - 0,7$) và phát dọn dây leo cục bộ quanh cây CTV.
  • Ban Quản lý Rừng & Doanh nghiệp Lâm nghiệp: Tiết kiệm ngân sách phục hồi rừng hàng tỷ đồng nhờ áp dụng phương án khoanh nuôi định lượng thay cho trồng mới cơ học.
  • Chính quyền địa phương (UBND xã Yên Lãng): Có luận cứ khoa học để quy hoạch 2.647 ha đất lâm nghiệp gắn liền với bảo tồn vùng đệm Vườn Quốc gia Tam Đảo.

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

1. Rừng phục hồi trạng thái IIA có những đặc điểm nhận biết lâm học nào ngoài thực địa?

Theo phân loại của Loeschau (1966), rừng IIA là giai đoạn phục hồi sau nương rẫy hoặc sau khai thác kiệt, tầng cây gỗ có mật độ $>1.000\text{ cây/ha}$ với $D_{1.3} > 10\text{ cm}$, trữ lượng gỗ thấp dưới $30\text{ m}^3\text{/ha}$, thành phần chiếm ưu thế là các loài tiên phong ưa sáng (Thành ngạnh, Thẩu tấu, Hu đay) xen lẫn cây con chịu bóng tái sinh dưới tán.

2. Tại sao mật độ cây tái sinh ở cấp chiều cao > 3m lại suy giảm mạnh (chỉ còn 6,61%)?

Đây là hệ quả tất yếu của quá trình đào thải tự nhiên và cạnh tranh không gian sống khốc liệt. Ở cấp $<1\text{ m}$, cây con sinh trưởng ồ ạt; khi vươn tới cấp $>1,5\text{ m}$, nhu cầu ánh sáng và dinh dưỡng tăng vọt, trong khi tầng tán cây tiên phong che khuất khiến những cây yếu ớt bị chết dần.

3. Tỷ lệ cây tái sinh hạt đạt 56,49% có ý nghĩa gì đối với công tác lâm sinh?

Cây mọc từ hạt có bộ rễ cọc phát triển sâu, khả năng chống chịu bão gió, hạn hán và sâu bệnh vượt trội hơn cây mọc từ chồi. Tỷ lệ hạt $>50%$ khẳng định lâm phần có tính ổn định di truyền cao và đủ năng lực hình thành tầng rừng gỗ lớn bền vững.

4. Khi nào cần tiến hành trồng dặm bổ sung trong rừng phục hồi IIA?

Trồng dặm chỉ áp dụng khi mật độ cây tái sinh mục đích triển vọng ($H_{vn} \ge 1\text{ m}$) đạt dưới $500\text{ cây/ha}$ hoặc tại các khoảng trống trảng cỏ/dây leo lớn không có nguồn giống tự nhiên phân tán tới.

5. Chi phí khoanh nuôi xúc tiến tái sinh có tiết kiệm hơn trồng rừng mới không?

Có. Khoanh nuôi xúc tiến tái sinh tự nhiên tận dụng tối đa cây con có sẵn và đất rừng tự nhiên, giúp tiết kiệm $60 - 70%$ chi phí đầu tư ban đầu so với việc xử lý thực bì, cuốc hố và mua cây giống trồng mới hoàn toàn.


Kết luận

Nghiên cứu về đặc điểm tái sinh tự nhiên trạng thái rừng phục hồi IIA tại xã Yên Lãng, huyện Đại Từ, tỉnh Thái Nguyên đã cung cấp hệ thống luận cứ khoa học và số liệu thực chứng chuẩn xác về hệ sinh thái rừng thứ sinh nhiệt đới. Với mật độ tái sinh tự nhiên trung bình đạt $4.133\text{ cây/ha}$, tỷ lệ cây tái sinh nguồn gốc hạt chiếm đa số ($56,49% - 91,67%$) cùng tỷ lệ cây triển vọng đạt trên $32%$, trạng thái rừng IIA tại khu vực có tiềm năng tự phục hồi rất to lớn. Việc áp dụng đồng bộ các giải pháp lâm sinh can thiệp chọn lọc (chặt tỉa cây sâu bệnh chèn ép, phát dọn thảm tươi quanh gốc cây mục đích và bảo vệ nghiêm ngặt) sẽ tối ưu hóa tốc độ phục hồi, nâng cao chất lượng rừng và bảo tồn đa dạng sinh học bền vững cho vùng trung du miền núi phía Bắc.