Giới thiệu dự án

Suy thoái diện tích và chất lượng rừng nhiệt đới đang là một trong những thách thức môi trường cấp bách nhất toàn cầu. Tại Việt Nam, thống kê lịch sử lâm nghiệp cho thấy năm 1943 độ che phủ rừng đạt khoảng 43% (tương đương 14,3 triệu ha), nhưng đến năm 1999 đã suy giảm xuống còn 33,2% (10,9 triệu ha, trong đó chỉ còn 9,4 triệu ha rừng tự nhiên). Hiện tượng mất rừng và suy kiệt trữ lượng bắt nguồn từ các phương thức khai thác kiệt thiếu cơ sở sinh thái, kết hợp với tập quán du canh nương rẫy và áp lực sinh kế dân sinh tại vùng đệm các khu bảo tồn.

Khu bảo tồn loài và sinh cảnh Nam Xuân Lạc (KBTL&SCNXL), huyện Chợ Đồn, tỉnh Bắc Kạn được thành lập theo Quyết định số 342/QĐ-UB ngày 17/03/2004 với quy mô 1.788 ha, giữ vai trò hành lang sinh thái trọng yếu kết nối Vườn Quốc gia Ba Bể và Khu Bảo tồn Thiên nhiên Nà Hang. Mặc dù sở hữu nhiều loài đặc hữu nguy cấp (Voọc đen má trắng, Voọc mũi hếch, Trai lý, Nghiến, Lan hài), hệ sinh thái rừng tại đây đang chịu sức ép nghiêm trọng từ 1.709 nhân khẩu vùng đệm và 6 hộ dân sinh sống tại vùng lõi với các hoạt động khai thác gỗ lén lút, thu hái củi đốt (bình quân 300 - 350 kg/hộ/tháng), chăn thả rông gia súc và phát rẫy luân chuyển.

Vấn đề cốt lõi đặt ra là sự thiếu hụt các dữ liệu định lượng có hệ thống về quy luật tái sinh tự nhiên, cấu trúc tầng tán và mối quan hệ sinh thái giữa tầng cây cao với lớp cây tái sinh thân gỗ trên các đai cao khác nhau. Đề tài khóa luận tốt nghiệp "Nghiên cứu đặc điểm tái sinh tự nhiên của thực vật thân gỗ tại Khu bảo tồn loài và sinh cảnh Nam Xuân Lạc, huyện Chợ Đồn, tỉnh Bắc Kạn" được thực hiện nhằm giải quyết bài toán này.

+---------------------------------------------------------------------------------------------------+
|                                MỤC TIÊU VÀ PHẠM VI NGHIÊN CỨU DỰ ÁN                               |
+---------------------------------------------------------------------------------------------------+
| 1. Định lượng cấu trúc tổ thành, mật độ, chất lượng và nguồn gốc tầng cây tái sinh.               |
| 2. Phân tích quy luật phân bố cây tái sinh theo 3 cấp chiều cao (<50cm, 50-100cm, >100cm).        |
| 3. Đánh giá tương quan giữa tái sinh với tầng cây cao, thảm tươi, thổ nhưỡng và nhân sinh.       |
| 4. Đề xuất các giải pháp kỹ thuật lâm sinh và quản trị bảo tồn nguồn gen cây gỗ quý hiếm.         |
+---------------------------------------------------------------------------------------------------+

Phương pháp tiếp cận dựa trên việc khảo sát thực địa thông qua hệ thống ô tiêu chuẩn (ÔTC) và ô dạng bản (ÔDB) phân bố theo 4 trạng thái sinh thái đại diện: Rừng kín thường xanh núi đất < 600m (I.1), Rừng kín thường xanh 600 - 800m (I.2), Rừng kín thường xanh mưa mùa á nhiệt đới > 800m (I.3), và Rừng thứ sinh phục hồi sau nương rẫy (I.4). Kết quả mang lại cơ sở dữ liệu định lượng chuẩn xác phục vụ trực tiếp cho công tác khoanh nuôi, xúc tiến tái sinh tự nhiên có tác động lâm sinh.


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

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

Trước đây, công tác đánh giá tái sinh rừng thường áp dụng các phương pháp mô tả định tính hoặc chỉ thống kê trữ lượng lâm phần thuần túy, dẫn đến các biện pháp can thiệp kỹ thuật kém hiệu quả.

Tiêu chí phân tích Phương pháp truyền thống (Định tính) Giải pháp sinh thái học (Định lượng chuẩn hóa)
Độ chính xác tổ thành Ước lượng cảm quan độ phong phú Tính chỉ số quan trọng $IV_i%$ và hệ số tổ thành $K_i$
Phân tầng cây con Đếm gộp tổng số cây con Phân loại 3 cấp chiều cao ($<0,5m; 0,5-1m; >1m$)
Nguồn gốc tái sinh Không phân biệt hạt/chồi Định lượng tỷ lệ hạt và chồi theo từng trạng thái
Tương quan thảm phủ Đánh giá chung chung độ rậm rạp Định lượng độ che phủ thảm tươi ($73% - 76,7%$)
Khả năng ứng dụng Khó lập kế hoạch khoanh nuôi chi tiết Thiết kế can thiệp lâm sinh trúng đích theo loài mục đích

Phân tích yêu cầu hệ thống điều tra và quản lý lâm sinh theo mô hình MoSCoW:

  • Must have (Bắt buộc): Đo đếm chuẩn xác đường kính ngang ngực ($D_{1.3}$), chiều cao vút ngọn ($H_{vn}$), xác định công thức tổ thành theo phương pháp Daniel Marmillod, phân loại phẩm chất cây tái sinh (Tốt/Trung bình/Xấu).
  • Should have (Nên có): Định vị tọa độ GPS các ô tiêu chuẩn, phân tích tương quan giữa mật độ tầng cây cao ($N/ha$) và mật độ lớp cây mạ/cây có triển vọng.
  • Could have (Có thể có): Xây dựng biểu đồ mặt cắt cấu trúc đứng của thảm thực vật cải tiến theo phương pháp David - Richards.
  • Won't have (Không làm đợt này): Mô phỏng chuỗi biến động diễn thế đa thời gian bằng ảnh viễn thám siêu phổ vệ tinh.

Thiết kế hệ thống

Mô hình kiến trúc thu thập, xử lý và phân tích số liệu lâm sinh được chuẩn hóa theo sơ đồ sau:

[Khảo sát thực địa ÔTC (500m2) & ÔDB (25m2)] 
[Thu thập dữ liệu: D1.3, Hvn, Loài, Cấp H, Nguồn gốc, Chất lượng]
[Module xử lý số liệu & Thống kê sinh học (Biometrics / R 4.3.1 / Python 3.10)]
[Hệ thống quyết định giải pháp kỹ thuật lâm sinh & Bản đồ số hóa QGIS 3.28 LTR]
  • Công nghệ sử dụng: R phiên bản 4.3.1 (thư viện vegan, forestecology), Python 3.10 (thư viện pandas, scipy, matplotlib), phần mềm Hệ thông tin địa lý QGIS 3.28 LTR phục vụ số hóa ranh giới tiểu khu 254, 255, 257.
  • Tiêu chuẩn bảo mật và an toàn dữ liệu: Dữ liệu điều tra được lưu trữ dưới định dạng bảng chuẩn hóa CSV UTF-8, đồng bộ kiểm định chéo giữa nhật ký đi tuyến và phiếu điều tra thực địa.

Methodology

Phương pháp nghiên cứu kết hợp giữa kỹ thuật lâm học ngoại nghiệp và phân tích thống kê sinh học nội nghiệp:

  • Quy trình ngoại nghiệp:
    • Lập ÔTC điển hình tạm thời diện tích $500 m^2$ ($20m \times 25m$) cho từng trạng thái rừng.
    • Lập 5 ÔDB diện tích $25 m^2$ ($5m \times 5m$) phân bố đều trên 2 đường chéo của mỗi ÔTC để đo đếm tầng cây tái sinh, cây bụi và thảm tươi.
  • Quy trình nội nghiệp: Tính toán các chỉ số sinh thái theo thuật toán chuẩn hóa.
import pandas as pd
import numpy as np

def calculate_forest_indices(df_trees, df_regen, area_otc_m2=500, area_odb_m2=125):
    """
    Thuật toán tính toán chỉ số IVI% và mật độ cây tái sinh theo Daniel Marmillod.
    """
    # 1. Tính toán tầng cây cao
    df_trees['G'] = (np.pi / 4) * (df_trees['D13_cm'] / 100) ** 2  # Tiết diện ngang (m2)
    total_G = df_trees['G'].sum()
    total_N = len(df_trees)
    
    species_summary = df_trees.groupby('Species').agg(
        Ni=('Species', 'count'),
        Gi=('G', 'sum')
    ).reset_index()
    
    species_summary['Ni_pct'] = (species_summary['Ni'] / total_N) * 100
    species_summary['Gi_pct'] = (species_summary['Gi'] / total_G) * 100
    species_summary['IVI_pct'] = (species_summary['Ni_pct'] + species_summary['Gi_pct']) / 2
    
    # Lọc loài tham gia công thức tổ thành (IVI% >= 5%)
    composition = species_summary[species_summary['IVI_pct'] >= 5.0].sort_values(
        by='IVI_pct', ascending=False
    )
    
    # 2. Tính toán mật độ cây tái sinh (N/ha)
    regen_count = len(df_regen)
    density_ha = (regen_count / area_odb_m2) * 10000
    
    return composition, density_ha

Implementation và kết quả

Development process

Quá trình nghiên cứu được chia làm 3 giai đoạn chặt chẽ:

  1. Giai đoạn 1 (Chuẩn bị và tiền trạm): Kế thừa bản đồ địa hình, bản đồ quy hoạch ranh giới 1.788 ha, phân tích dữ liệu khí tượng thủy văn trạm Chợ Đồn (lượng mưa bình quân $>1500 mm$, nhiệt độ trung bình $20 - 22^\circ\text{C}$).
  2. Giai đoạn 2 (Thi công thực địa): Khảo sát 12 ÔTC tại 4 trạng thái rừng trên núi đất từ độ cao chân núi đến đỉnh Tam Sao ($1.159 m$). Giám định tiêu bản các loài khó định danh như Machilus platycarpa, Zenia insignis, Excentrodendron tonkinense.
  3. Giai đoạn 3 (Xử lý dữ liệu): Tổng hợp phiếu điều tra, phân loại chất lượng cây con theo tiêu chuẩn sinh lực (thân thẳng, không cụt ngọn, không sâu bệnh).

Testing và validation

Dữ liệu điều tra tại các trạng thái rừng được kiểm tra tính đại diện và kiểm chuẩn sai số mẫu ngẫu nhiên:

Trạng thái rừng điều tra Ký hiệu Độ cao tuyệt đối (m) Số ÔTC lập Tổng diện tích ÔDB ($m^2$) Mật độ tầng cao ($N/ha$)
Kín thường xanh núi thấp I.1 $< 600$ 1 125 1.260
Kín thường xanh núi trung bình I.2 $600 - 800$ 7 875 1.540
Kín thường xanh á nhiệt đới I.3 $> 800$ 1 125 1.980
Thứ sinh phục hồi sau nương rẫy I.4 $400 - 700$ 3 375 1.786

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

Nghiên cứu đã thu thập và lượng hóa chi tiết các chỉ số sinh thái học của thảm thực vật tái sinh:

+---------------------------------------------------------------------------------------------------+
|                       TỔNG HỢP MẬT ĐỘ VÀ CHẤT LƯỢNG TÁI SINH THEO TRẠNG THÁI                      |
+---------------------------------------------------------------------------------------------------+
| Trạng thái  | Mật độ TSR  | Cấp <50cm   | Cấp 50-100cm | Cấp >100cm  | Tốt (%) | Hạt (%) | Chồi (%) |
|-------------+-------------+-------------+--------------+-------------+---------+---------+----------|
| I.1 (<600m) | 4.560 c/ha  | 1.200 c/ha  | 2.880 c/ha   | 480 c/ha    | 35,5%   | 38,2%   | 61,8%    |
| I.2 (600-800| 3.592 c/ha  | 1.714 c/ha  | 1.348 c/ha   | 822 c/ha    | 9,4%    | 36,4%   | 63,6%    |
| I.3 (>800m) | 2.560 c/ha  | 480 c/ha    | 400 c/ha     | 1.680 c/ha  | 63,0%   | 42,1%   | 57,9%    |
| I.4 (Rẫy)   | 2.880 c/ha  | 5.546 c/ha* | 1.973 c/ha   | 1.573 c/ha  | 30,2%   | 39,5%   | 60,5%    |
+---------------------------------------------------------------------------------------------------+
(*Ghi chú: Cấp <50cm ở trạng thái I.4 đạt mật độ hạt nảy mầm/mạ cục bộ cao nhưng tỷ lệ đào thải sinh học lớn).
  • Tổ thành loài tái sinh:
    • Trạng thái I.1: Tham gia công thức tổ thành gồm 6-7 loài, chiếm ưu thế là Hà nu, Muồng trắng, Dẻ gai; xuất hiện cây con loài quý hiếm Trai lý (Garcinia fragraeoides).
    • Trạng thái I.2: Độ đa dạng cao nhất (có ô xuất hiện 8 loài tham gia tổ thành), ưu thế gồm Cà lồ, Hà nu, Nhọc, Găng Việt Nam.
    • Trạng thái I.3: Độ phong phú loài giảm, ưu thế bởi các loài chịu lạnh vùng núi cao: Hồ đào núi, Kháo (Machilus), Trai đỏ, có sự xuất hiện của Nghiến (Excentrodendron tonkinense).
    • Trạng thái I.4: Ưu thế các loài tiên phong ưa sáng: Bùm bụm, Kẹn, Mán đỉa, Trâm, Xoan nhừ.
  • Đặc điểm thảm tươi, cây bụi: Độ che phủ trung bình dao động từ $73,0%$ (I.3) đến $76,7%$ (I.2). Chiều cao thảm tươi phổ biến $< 100 cm$, đóng vai trò giữ ẩm đất, chống xói mòn bề mặt trên sườn dốc $25 - 45^\circ$, không gây cản trở gay gắt đến tầng cây tái sinh đã vượt qua giai đoạn mạ.

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

  • Định lượng hóa quy luật phân bố sinh thái theo đai cao: Chứng minh rõ quy luật mật độ cây tái sinh giảm dần từ chân núi lên đỉnh núi (từ $4.560$ cây/ha ở đai $<600m$ xuống $2.560$ cây/ha ở đai $>800m$), nhưng tỷ lệ cây có chất lượng tốt và nguồn gốc từ hạt lại tăng vọt ở đai cao (đạt $63,0%$ chất lượng tốt và $42,1%$ tái sinh hạt ở đỉnh Tam Sao).
  • So sánh với các nghiên cứu tiền nhiệm:
    • So với Thái Văn Trừng (1978): Nghiên cứu đã kiểm chứng mô hình cấu trúc 4-5 tầng tại vùng núi đất Bắc Kạn, đồng thời lượng hóa được chỉ số $IV_i%$ chi tiết cho từng phân tầng thay vì chỉ phân loại định tính.
    • So với Lê Đồng Tấn (1999) & Phạm Ngọc Thường (2003): Khẳng định tính tương đồng về sự suy giảm mật độ theo sườn dốc nương rẫy, nhưng bổ sung phát hiện mới về sự phục hồi mạnh mẽ của cây tái sinh mục đích (đạt $1.573$ cây triển vọng $>1m$/ha tại trạng thái I.4 sau 18-19 năm bỏ hóa).
  • Đóng góp thực tiễn cho bảo tồn nguồn gen: Định vị chính xác các tiểu quần thể tái sinh tự nhiên của loài cực kỳ nguy cấp như Nghiến, Trai lý, Lan kim tuyến đá vôi (Anoectochilus calcareus), Hoàng tinh cách, Giổi lụa (Tsoongiodendron odorum).

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

Tình huống triển khai thực địa

Kết quả nghiên cứu cung cấp căn cứ khoa học trực tiếp cho Ban Quản lý Khu bảo tồn Loài và Sinh cảnh Nam Xuân Lạc triển khai các biện pháp lâm sinh phân vùng:

[Phân vùng can thiệp lâm sinh]

Phân tích hiệu quả kinh tế - xã hội và ROI

  • Giảm chi phí trồng mới: Chi phí trồng rừng mới trên núi đá/núi đất dốc dao động từ 25 - 30 triệu VNĐ/ha. Áp dụng giải pháp khoanh nuôi xúc tiến tái sinh tự nhiên có làm giàu rừng chỉ tốn khoảng 6 - 8 triệu VNĐ/ha (tiết kiệm $\approx 70%$ ngân sách lâm nghiệp).
  • Giảm áp lực sinh kế: Đề xuất mô hình chuyển giao kỹ thuật trồng cây dược liệu tán rừng (Lan kim tuyến, Ba kích, Gừng đá) và nhân rộng 100% bếp đun cải tiến tiết kiệm củi giúp giảm lượng khai thác gỗ củi từ 350 kg/hộ/tháng xuống dưới 120 kg/hộ/tháng cho 410 hộ vùng đệm.

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

  • Hạn chế kỹ thuật: Thời gian theo dõi điều tra thực địa tập trung trong giai đoạn ngắn (tháng 09 đến tháng 12 năm 2014), chưa phản ánh hết biến động động thái tái sinh theo chu kỳ sai quả - mất mùa hạt giống qua các năm khác nhau.
  • Ràng buộc tài nguyên: Chưa thực hiện được việc gắn thẻ cố định (tagging) để theo dõi tốc độ tăng trưởng đường kính và chiều cao cá thể cây con dài hạn theo ô định vị cố định (Permanent Sample Plot - PSP).
  • Hướng nghiên cứu tiếp theo:
    • Thiết lập mạng lưới ô định vị cố định PSP 1 ha để quan trắc biến động quần xã thực vật liên tục 5 - 10 năm.
    • Ứng dụng mô hình hóa không gian GIS và cảm biến viễn thám LiDAR để dự báo sinh khối tích lũy carbon của tầng cây tái sinh.

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

+---------------------------------------------------------------------------------------------------+
|                                 GIÁ TRỊ CHUYỂN GIAO CHO CÁC NHÓM ĐỐI TƯỢNG                         |
+---------------------------------------------------------------------------------------------------+
| 1. Sinh viên ngành Lâm nghiệp / Sinh thái: Bộ dữ liệu mẫu chuẩn về phương pháp điều tra ÔTC/ÔDB.   |
| 2. Kỹ sư lâm sinh & Nhà nghiên cứu: Thuật toán và công thức tính toán chỉ số IVI%, Ki chuẩn hóa.   |
| 3. Ban Quản lý Rừng đặc dụng: Bản đồ hiện trạng và giải pháp kỹ thuật khoanh nuôi trúng đích.    |
| 4. Cộng đồng dân cư vùng đệm: Cơ chế hưởng lợi từ chính sách chi trả DVMTR và giao khoán bảo vệ.  |
+---------------------------------------------------------------------------------------------------+
  • Sinh viên & Học viên cao học: Cung cấp tài liệu tham khảo chuẩn mực về cấu trúc rừng và quy trình xử lý số liệu điều tra tái sinh tự nhiên nhiệt đới.
  • Ban Quản lý Khu bảo tồn & Chi cục Kiểm lâm: Có cơ sở khoa học phê duyệt kế hoạch tác nghiệp lâm sinh giai đoạn 2025 - 2030, tối ưu hóa nguồn vốn chương trình phát triển lâm nghiệp bền vững.
  • Cộng đồng địa phương: Tạo điều kiện phân bổ nguồn kinh phí khoanh nuôi bảo vệ rừng từ quỹ dịch vụ môi trường rừng (DVMTR), hỗ trợ ổn định sinh kế cho 89,5% đồng bào dân tộc thiểu số (Dao, Tày, Mông).

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

1. Điều kiện kỹ thuật và mật độ tối thiểu để khoanh nuôi tái sinh tự nhiên thành rừng là gì?

Mật độ cây tái sinh mục đích có chiều cao $>1m$ (cây có triển vọng) phải đạt tối thiểu từ $500 - 1.000$ cây/ha, phân bố tương đối đồng đều trên lâm phần, và độ tàn che của tầng cây cao đạt từ $0,3 - 0,5$.

2. Tại sao mật độ cây tái sinh ở đai cao > 800m thấp nhất nhưng chất lượng tốt lại cao nhất?

Đai cao $>800m$ có điều kiện khí hậu khắc nghiệt (nhiệt độ mùa đông xuống tới $4^\circ\text{C}$, gió mạnh, sương muối), áp lực chọn lọc tự nhiên cao đã đào thải sớm các cây mạ yếu ớt. Đồng thời, khu vực này hầu như không chịu tác động chặt phá hay chăn thả của con người, giúp cây con phát triển thân thẳng, không bị tổn thương cơ học, đạt tỷ lệ cây tốt $63,0%$.

3. Giải pháp nào ngăn chặn hiện tượng khai thác củi đun ảnh hưởng đến tầng cây con?

Triển khai đồng bộ 3 giải pháp: (1) Hỗ trợ xây dựng bếp đun cải tiến tiết kiệm $65%$ lượng củi; (2) Quy hoạch vùng thu nhặt củi khô gãy mục tại rừng sản xuất; (3) Nghiêm cấm hành vi chặt hạ cây gỗ tái sinh để tạo củi khô mục tiêu thụ sau.

4. Chi phí đầu tư hệ sinh thái số hóa phục vụ quản lý tái sinh rừng là bao nhiêu?

Sử dụng hoàn toàn nền tảng phần mềm mã nguồn mở (QGIS, R, Python) với chi phí bản quyền 0 VNĐ. Thiết bị phần cứng chỉ yêu cầu GPS cầm tay tiêu chuẩn và máy tính bảng dã sinh thu thập số liệu ngoại nghiệp.

5. Khả năng nhân rộng quy trình điều tra này cho các khu bảo tồn khác như thế nào?

Quy trình điều tra ÔTC $500 m^2$ kết hợp 5 ÔDB $25 m^2$ theo 2 đường chéo hoàn toàn tương thích với Tiêu chuẩn Quốc gia (TCVN) về điều tra rừng nhiệt đới, có thể áp dụng ngay cho toàn bộ hệ thống rừng đặc dụng trên núi đất và núi đá vôi tại vùng Đông Bắc và Tây Bắc Việt Nam.


Kết luận

Công trình nghiên cứu đã hoàn thành xuất sắc các mục tiêu đề ra, làm sáng tỏ các đặc trưng cấu trúc lâm phần và quy luật tái sinh tự nhiên của thực vật thân gỗ tại Khu bảo tồn loài và sinh cảnh Nam Xuân Lạc. Các phát hiện định lượng then chốt gồm: mật độ tái sinh tự nhiên dao động từ $2.560 - 4.560$ cây/ha; đai cao $>800m$ sở hữu phẩm chất cây con tốt nhất ($63,0%$); và trạng thái rừng phục hồi sau nương rẫy có khả năng tự tái thiết mạnh mẽ với $1.573$ cây triển vọng/ha.

Nghiên cứu khẳng định tiềm năng phục hồi to lớn của hệ sinh thái rừng Nam Xuân Lạc nếu được áp dụng đúng đắn các biện pháp khoanh nuôi xúc tiến tái sinh kết hợp hỗ trợ sinh kế cộng đồng vùng đệm. Đây là tài liệu khoa học - kỹ thuật có giá trị ứng dụng thực tiễn cao cho công tác lâm sinh học, quản lý rừng bền vững và bảo tồn đa dạng sinh học tại Việt Nam.