Giới thiệu dự án

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

Theo số liệu thống kê từ các tổ chức bảo tồn quốc tế (IUCN, UNDP và WWF), trung bình mỗi năm thế giới mất khoảng 20 triệu hecta rừng tự nhiên, trong đó hoạt động đốt nương làm rẫy chiếm tới 50%, cháy rừng chiếm 23% và khai thác thương mại chiếm từ 5–7%. Tại Việt Nam, suy thoái rừng diễn ra nghiêm trọng tại các vùng núi đá vôi phía Bắc, nơi địa hình chia cắt phức tạp và đời sống người dân phụ thuộc lớn vào việc khai thác lâm sản.

Loài Nghiến gân ba (Excentrodendron tonkinensis (Gagnep.) Chang & Miau, đồng danh Burretiodendron hsienmu Ching & Hu) thuộc họ Đay (Tiliaceae / Malvaceae s.l.) là loài cây gỗ lớn bản địa quý hiếm thuộc Nhóm IIA theo Nghị định 06/2019/NĐ-CP và Sách Đỏ Việt Nam. Nghiến gân ba giữ vai trò cấu trúc sinh thái đặc hữu trên các dải karst nhiệt đới. Tại Khu bảo tồn thiên nhiên Thần Sa – Phượng Hoàng (huyện Võ Nhai, tỉnh Thái Nguyên), quần thể Nghiến gân ba đang đứng trước nguy cơ suy thoái cục bộ do áp lực khai thác lịch sử và tập quán canh tác nông nghiệp nương rẫy, dẫn đến suy giảm khả năng tái sinh tự nhiên dưới tán rừng phục hồi.

Vấn đề nghiên cứu và các điểm nghẽn

  • Thiếu hụt dữ liệu định lượng vi mô: Chưa có các đánh giá sinh thái định lượng về cấu trúc tổ thành, mật độ, chất lượng và nguồn gốc của lớp cây tái sinh đi kèm loài Nghiến gân ba trên các tiểu địa hình khác nhau (chân núi, sườn núi, đỉnh núi đá vôi).
  • Tỷ lệ suy vong cây mạ cao: Cây con tái sinh từ hạt của loài Nghiến gân ba gặp trở ngại lớn do tầng thảm tươi dày, sự cạnh tranh dinh dưỡng của các loài cây gỗ mềm ưa sáng mọc nhanh và hiện tượng khô hạn sinh lý trên bề mặt đá vôi.
  • Can thiệp lâm sinh thiếu chuẩn hóa: Các biện pháp khoanh nuôi phục hồi theo Quy phạm Kỹ thuật Lâm sinh QPN 14-92 chưa được tinh chỉnh riêng cho hệ sinh thái rừng núi đá vôi có loài Nghiến gân ba chiếm ưu thế sinh thái.

Mục tiêu dự án

  1. Định lượng cấu trúc tầng cây gỗ: Xác định mật độ, tiết diện ngang và chỉ số giá trị quan trọng ($IVI%$) của các loài cây gỗ lớn trong 9 ô tiêu chuẩn (OTC) đại diện tại xã Thần Sa.
  2. Đánh giá động thái lớp cây tái sinh: Thu thập và phân tích cấu trúc tổ thành loài ($K_i$), mật độ cây tái sinh ($N/ha$), tỷ lệ cây triển vọng ($CTV%$), phẩm chất (Tốt, Trung bình, Xấu) và nguồn gốc (Hạt, Chồi) trên 45 ô thứ cấp (ODB).
  3. Xác định chỉ số đa dạng sinh học và quy luật phân bố: Đo lường độ đa dạng loài qua chỉ số Shannon–Weaver ($H'$) và quy luật phân bố số cây, số loài theo 7 cấp chiều cao.
  4. Xây dựng giải pháp kỹ thuật lâm sinh: Đề xuất quy trình tác động khoanh nuôi xúc tiến tái sinh kết hợp trồng bổ sung phù hợp với điều kiện lập địa karst Thần Sa.

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

Ứng dụng quan điểm sinh thái 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) kết hợp hệ phương pháp điều tra lâm học định lượng tiêu chuẩn. Dự án triển khai mạng lưới 9 OTC tạm thời diện tích $1.000\text{ m}^2$ ($40\text{m} \times 25\text{m}$) phân tầng theo 3 cấp vị trí địa hình (chân núi, sườn núi, đỉnh núi) tại độ cao 170m – 353m. Trong mỗi OTC, bố trí 5 ô thứ cấp diện tích $25\text{m}^2$ ($5\text{m} \times 5\text{m}$) tại 4 góc và trung tâm để giải mã toàn diện tầng tái sinh ($D_{1.3} < 6\text{ cm}$).

[Khảo sát thực địa OTC 1000m²]

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

  • Xây dựng hoàn chỉnh cơ sở dữ liệu số hóa của 9 OTC với mật độ tầng cây gỗ dao động từ $280 - 410\text{ cây/ha}$.
  • Xác định chính xác mật độ tái sinh tự nhiên từ $4.640 - 8.480\text{ cây/ha}$ và mật độ cây triển vọng ($H \ge 1.0\text{ m}$) đạt từ $1.440 - 2.560\text{ cây/ha}$.
  • Thiết lập mô hình tương quan giữa độ cao/địa hình với chỉ số đa dạng sinh học Shannon ($H'$ từ $0{,}98 - 1{,}25$) và tỷ lệ nguồn gốc tái sinh (Hạt vs. Chồi).

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

  • Không gian: Lâm phần rừng tự nhiên phục hồi sau nương rẫy thuộc xã Thần Sa, huyện Võ Nhai, tỉnh Thái Nguyên (tọa độ $105^\circ 45' - 106^\circ 17'$ Đ, $21^\circ 36' - 21^\circ 56'$ B).
  • Thời gian: Số liệu thực địa thu thập liên tục từ tháng 01/2019 đến tháng 06/2019.
  • Đối tượng: Cây gỗ lớn có $D_{1.3} \ge 6\text{ cm}$ và toàn bộ cây gỗ tái sinh có $D_{1.3} < 6\text{ cm}$.

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

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

Tiêu chí so sánh Tái sinh tự nhiên hoàn toàn Trồng rừng nhân tạo thuần loài Xúc tiến tái sinh tự nhiên có hỗ trợ kỹ thuật (Giải pháp đề tài)
Chi phí đầu tư Rất thấp (chỉ bảo vệ) Rất cao (làm đất, cây giống, chăm sóc) Trung bình (điều tiết tán, phát dọn, tra dặm)
Độ bền vững sinh thái Phụ thuộc hoàn toàn tự nhiên Kém, dễ bùng phát sâu bệnh, xói mòn Rất cao, bảo tồn cấu trúc rừng nhiều tầng tán
Bảo tồn nguồn gen bản địa Tự phát, nguy cơ tuyệt chủng loài quý Đa phần nhập nội (Keo, Bạch đàn) Bảo tồn trực tiếp loài đặc hữu (E. tonkinensis)
Thời gian phục hồi Kéo dài (30–50 năm) Ngắn (chu kỳ 5–10 năm khai thác) Tối ưu hóa (15–20 năm đạt rừng ổn định)
Khả năng kiểm soát tổ thành Rất thấp, cây tái sinh mục đích bị lấn át Tuyệt đối theo ý muốn con người Cao, điều chỉnh qua nuôi dưỡng cây triển vọng

Phân loại yêu cầu giải pháp kỹ thuật theo MoSCoW

  • Must Have (Bắt buộc):
    • Thu thập dữ liệu đầy đủ trên $9.000\text{ m}^2$ diện tích OTC đại diện và $1.125\text{ m}^2$ diện tích ô thứ cấp.
    • Định danh chính xác các loài cây tái sinh mục đích và loài đi kèm.
    • Tính toán chỉ số $IVI%$, hệ số tổ thành $K_i$, chỉ số $H'$ theo các công thức sinh thái chuẩn mực.
  • Should Have (Nên có):
    • Mô hình hóa tương quan giữa vị trí địa hình (Chân - Sườn - Đỉnh) với mật độ cây triển vọng.
    • Phân tích cự ly phân bố không gian giữa các cá thể cây tái sinh lân cận (30 cự ly/trạng thái).
  • Could Have (Có thể có):
    • Số hóa bản đồ không gian phân bố OTC trên hệ thống thông tin địa lý GIS.
  • Won't Have (Chưa thực hiện đợt này):
    • Xét nghiệm ADN nhận dạng biến dị di truyền của các dòng Nghiến gân ba tại vùng lõi.

Thiết kế hệ thống xử lý số liệu sinh thái

Quy trình tính toán và phân tích cấu trúc lâm phần được tổ chức thành các module liên kết logic:

Ngăn xếp công nghệ và công cụ thống kê sinh thái

  • Phần mềm xử lý dữ liệu: R Core Team (v4.2.0) với các packages chuyên dụng vegan (v2.6-4) để tính toán đa dạng sinh học, tidyverse (v2.0.0) để chuyển đổi và làm sạch dữ liệu.
  • Xử lý bảng tính & thống kê: Microsoft Excel 2019 / Python 3.10 (thư viện pandas v1.5.3, scipy v1.10.1).
  • Tiêu chuẩn ngành áp dụng: Quy phạm kỹ thuật lâm sinh QPN 14-92 (Bộ Lâm nghiệp ban hành theo Quyết định số 200/QĐ-KT).

Đánh giá rủi ro và phương án giảm thiểu

  1. Rủi ro sai số nhận diện loài thực vật thực địa: Các loài cây mạ/cây con họ Đay, Re, Thầu dầu rất dễ nhầm lẫn.
    • Giảm thiểu: Thu thập đầy đủ tiêu bản cành lá, hoa quả khô/tươi và đối chiếu tại Phòng tiêu bản Thực vật – Trường Đại học Nông Lâm Thái Nguyên.
  2. Rủi ro địa hình hiểm trở núi đá vôi: Vách đá dựng đứng, trơn trượt gây nguy hiểm trong quá trình đo đạc góc dốc và lập ô tiêu chuẩn.
    • Giảm thiểu: Sử dụng dây định vị an toàn, la bàn địa chất Suunto Tandem và phân định ranh giới OTC theo đúng chiếu bằng $1.000\text{ m}^2$.

Implementation và kết quả

Quy trình thu thập và xử lý dữ liệu

Quá trình nghiên cứu được triển khai theo 4 giai đoạn logic:

  • Giai đoạn 1 (Tháng 01/2019 – Tháng 02/2019): Điều tra tổng quan kinh tế - xã hội, điều kiện tự nhiên tại xã Thần Sa; định vị tọa độ 9 OTC đại diện trên các độ cao từ 170m đến 353m.
  • Giai đoạn 2 (Tháng 03/2019 – Tháng 04/2019): Thiết lập 45 ô thứ cấp ($25\text{ m}^2$); đo đếm chi tiết tầng cây gỗ lớn ($D_{1.3}$, $H_{vn}$, $H_{dc}$, $D_t$) và tầng cây tái sinh (chiều cao, chất lượng, gốc phát sinh, khoảng cách).
  • Giai đoạn 3 (Tháng 05/2019): Nhập liệu, chuẩn hóa dữ liệu, chạy mô hình thuật toán tính toán các chỉ số lâm học.
  • Giai đoạn 4 (Tháng 06/2019): Tổng hợp báo cáo, phân tích thảo luận và xây dựng hệ thống giải pháp kỹ thuật.

Các công thức toán học và mã nguồn tính toán

1. Chỉ số giá trị quan trọng ($IVI%$) tầng cây gỗ:

$$IVI_i (%) = \frac{N_i% + G_i%}{2}$$ Trong đó:

  • $N_i% = \frac{N_i}{\sum N} \times 100$: Tỷ lệ số cây của loài thứ $i$.
  • $G_i% = \frac{G_i}{\sum G} \times 100$: Tỷ lệ tổng tiết diện ngang thân của loài thứ $i$ ($G = \frac{\pi \times D_{1.3}^2}{4}$).

2. Chỉ số đa dạng sinh học Shannon–Weaver ($H'$):

$$H' = - \sum_{i=1}^{s} \left( \frac{n_i}{N} \ln \frac{n_i}{N} \right)$$ Trong đó: $s$ là tổng số loài, $n_i$ là số cá thể loài thứ $i$, $N$ là tổng số cá thể quan sát.

3. Tỷ lệ cây triển vọng ($CTV%$):

$$CTV% = \frac{\sum n_{(H \ge 1{,}0\text{m})}}{\sum N_i} \times 100%$$

Kịch bản Python tự động hóa phân tích sinh thái:

import numpy as np
import pandas as pd

def calculate_ecological_metrics(df_trees, df_regen):
    """
    df_trees: DataFrame chứa dữ liệu cây gỗ (OTC_ID, Species, D13, Hvn)
    df_regen: DataFrame chứa dữ liệu cây tái sinh (OTC_ID, Species, Height_m, Quality, Origin)
    """
    # 1. Tính toán IVI% tầng cây gỗ
    df_trees['G'] = (np.pi * (df_trees['D13'] / 100) ** 2) / 4  # Tiết diện ngang (m2)
    tree_summary = df_trees.groupby(['OTC_ID', 'Species']).agg(
        Ni=('D13', 'count'),
        Gi=('G', 'sum')
    ).reset_index()
    
    otc_totals = tree_summary.groupby('OTC_ID').agg(Total_N=('Ni', 'sum'), Total_G=('Gi', 'sum')).reset_index()
    tree_summary = tree_summary.merge(otc_totals, on='OTC_ID')
    tree_summary['Ni_pct'] = (tree_summary['Ni'] / tree_summary['Total_N']) * 100
    tree_summary['Gi_pct'] = (tree_summary['Gi'] / tree_summary['Total_G']) * 100
    tree_summary['IVI_pct'] = (tree_summary['Ni_pct'] + tree_summary['Gi_pct']) / 2
    
    # 2. Tính toán Shannon-Weaver Index (H') và Cây triển vọng tầng tái sinh
    regen_summary = []
    for otc_id, group in df_regen.groupby('OTC_ID'):
        total_regen = len(group)
        species_counts = group['Species'].value_counts().values
        proportions = species_counts / total_regen
        shannon_h = -np.sum(proportions * np.log(proportions))
        
        # Cây triển vọng: Chiều cao >= 1.0m và phẩm chất Tốt hoặc Trung bình
        ctv_count = len(group[(group['Height_m'] >= 1.0) & (group['Quality'].isin(['Tot', 'TB']))])
        ctv_pct = (ctv_count / total_regen) * 100
        
        regen_density_ha = (total_regen / 125) * 10000  # Quy đổi diện tích 5 ô thứ cấp (125m2) sang 1ha
        ctv_density_ha = (ctv_count / 125) * 10000
        
        regen_summary.append({
            'OTC_ID': otc_id,
            'Total_Species': len(species_counts),
            'Density_ha': regen_density_ha,
            'CTV_Density_ha': ctv_density_ha,
            'CTV_pct': ctv_pct,
            'Shannon_H': shannon_h
        })
        
    return tree_summary, pd.DataFrame(regen_summary)

Kết quả nghiên cứu thực địa chi tiết

1. Cấu trúc mật độ tầng cây gỗ lớn và Nghiến gân ba

Khảo sát 9 OTC ($1.000\text{ m}^2/\text{OTC}$) ghi nhận mật độ cây tầng gỗ lớn dao động từ $280 - 410\text{ cây/ha}$. Mật độ loài Nghiến gân ba (Excentrodendron tonkinensis) biến động từ $10 - 60\text{ cây/ha}$, chiếm tỷ lệ từ $2{,}4% - 15{,}4%$ tổng số cây gỗ trong lâm phần.

Số hiệu OTC Độ cao (m) Vị trí địa hình Tổng số cây/OTC Số cây Nghiến/OTC Mật độ tầng gỗ (cây/ha) Mật độ Nghiến (cây/ha) Tỷ lệ Nghiến (%)
OTC 1 293 Sườn 36 3 360 30 8,33%
OTC 2 353 Đỉnh 30 3 300 30 10,00%
OTC 3 170 Chân 41 1 410 10 2,44%
OTC 4 170 Chân 37 2 370 20 5,41%
OTC 5 320 Đỉnh 35 3 350 30 8,57%
OTC 6 315 Đỉnh 32 2 320 20 6,25%
OTC 7 273 Sườn 39 6 390 60 15,38%
OTC 8 286 Sườn 28 2 280 20 7,14%
OTC 9 173 Chân 33 3 330 30 9,09%
Trung bình 261,4 — 34,6 2,8 345,6 27,8 8,07%

2. Cấu trúc tổ thành và mật độ cây tái sinh theo địa hình

Vị trí địa hình Số hiệu OTC Độ cao (m) Mật độ tái sinh (cây/ha) Mật độ cây triển vọng (cây/ha) Tỷ lệ cây triển vọng ($CTV%$) Công thức tổ thành tái sinh ưu thế ($K_i \ge 5%$)
Chân núi OTC 3 170 8.480 2.560 30,19% $2{,}45\text{Ngh} + 2{,}08\text{Nhđ} + 1{,}98\text{Cbb} + 1{,}13\text{Hvn} + 0{,}66\text{Tbb} + 0{,}57\text{Thđ} + 1{,}13\text{Lk}$
OTC 4 170 8.400 2.560 30,48% $2{,}48\text{Mat} + 2{,}10\text{Vom} + 1{,}71\text{Mas} + 1{,}24\text{Hav} + 1{,}24\text{Mlt} + 0{,}86\text{Kha} + 0{,}38\text{Lk}$
OTC 9 173 7.520 2.160 28,72% $2{,}13\text{Hvn} + 1{,}49\text{Hav} + 1{,}28\text{Hav} + 1{,}06\text{Ngh} + 0{,}96\text{Cha} + 0{,}85\text{Mlt} + 0{,}64\text{Lom} + 1{,}61\text{Lk}$
Sườn núi OTC 7 273 7.200 2.000 27,78% $2{,}56\text{Cbb} + 2{,}44\text{Nhđ} + 1{,}89\text{Hvn} + 1{,}78\text{Ngh} + 1{,}33\text{Sog} + 1{,}71\text{Lk}$
OTC 8 286 6.720 2.320 34,52% $1{,}43\text{Nhđ} + 1{,}31\text{Hvn} + 1{,}19\text{Kha} + 1{,}07\text{Cbb} + 1{,}07\text{Ngh} + 0{,}95\text{Sog} + 0{,}71\text{Mas} + 0{,}6\text{Tbb} + 0{,}6\text{Trv} + 1{,}08\text{Lk}$
OTC 1 293 6.320 2.080 32,91% $2{,}17\text{Ngh} + 2{,}07\text{Nhđ} + 1{,}96\text{Mat} + 1{,}74\text{Hvn} + 1{,}09\text{Tbb} + \text{Sud} + 2{,}14\text{Lk}$
Đỉnh núi OTC 2 353 6.720 1.840 27,38% $2{,}50\text{Tnt} + 2{,}25\text{Cha} + 1{,}63\text{Map} + 1{,}25\text{Hvn} + 1{,}25\text{Ngh} + 0{,}88\text{Cbb} + 2{,}02\text{Lk}$
OTC 5 320 6.160 1.920 31,17% $1{,}96\text{Kha} + 1{,}96\text{Ngh} + 1{,}79\text{Thu} + 1{,}43\text{Trv} + 1{,}25\text{Duo} + 1{,}25\text{Thđ} + 1{,}07\text{Ddx} + 0{,}89\text{Ngl} + 0{,}71\text{Mas} + 0{,}72\text{Lk}$
OTC 6 315 4.640 1.440 31,03% $2{,}20\text{Kha} + 2{,}20\text{Thu} + 1{,}19\text{Duo} + 1{,}02\text{Mat} + 1{,}02\text{Ngl} + 1{,}02\text{Trv} + 0{,}85\text{Ngh} + 0{,}51\text{Vom} + 1{,}65\text{Lk}$

(Ký hiệu loài: Ngh: Nghiến, Nhđ: Nhọc đá/Nhọc đen, Mat: Màng tang, Cbb: Chay Bắc Bộ, Hvn: Hương viên núi/đá, Kha: Kháo, Duo: Dướng, Mas: Mạy sả, Cha: Châm, Tbb: Thích Bắc Bộ, Tnt: Thích năm thùy, Vom: Vỏ mản, Sog: Sồi gai, Lk: Loài khác).

3. Chỉ số đa dạng sinh học và chất lượng - nguồn gốc tái sinh

Chỉ số đa dạng Shannon–Weaver ($H'$) phản ánh tính phong phú và độ đồng đều về loài của thảm tái sinh tại Thần Sa dao động từ $0{,}9811$ đến $1{,}2493$.

OTC Mật độ TS (cây/ha) Số loài xuất hiện Chỉ số Shannon ($H'$) Cây tốt (cây/ha) Cây TB (cây/ha) Cây xấu (cây/ha) Nguồn gốc Hạt (cây/ha) Nguồn gốc Chồi (cây/ha) Tỷ lệ Chồi (%)
OTC 1 6.320 16 1,0954 1.853 1.234 3.233 2.546 3.774 59,72%
OTC 2 6.720 20 1,1722 2.341 1.963 2.416 2.890 3.830 56,99%
OTC 3 8.480 18 1,0756 2.487 2.156 3.837 3.517 4.963 58,53%
OTC 4 8.400 13 0,9811 2.160 3.241 2.999 3.298 5.102 60,74%
OTC 5 6.160 18 1,2042 1.254 2.451 2.455 3.212 2.948 47,86%
OTC 6 4.640 15 1,1192 1.345 2.148 1.147 2.750 1.890 40,73%
OTC 7 7.200 17 1,1636 2.869 3.127 1.204 3.214 3.986 55,36%
OTC 8 6.720 17 1,1648 2.136 1.894 2.690 2.378 4.342 64,61%
OTC 9 7.520 20 1,2493 2.341 2.130 3.049 3.567 3.953 52,57%

Nhận xét về động thái:

  • Vị trí địa hình: Mật độ cây tái sinh đạt cực đại tại vị trí chân núi ($7.520 - 8.480\text{ cây/ha}$), giảm dần ở sườn núi ($6.320 - 7.200\text{ cây/ha}$) và thấp nhất ở đỉnh núi ($4.640 - 6.720\text{ cây/ha}$). Lớp đất chân núi tích tụ nhiều mùn hữu cơ và độ ẩm cao hơn sườn và đỉnh núi đá vôi.
  • Nguồn gốc tái sinh: Tái sinh chồi chiếm ưu thế chủ đạo ($40{,}73% - 64{,}61%$), phản ánh đặc thù rừng sau nương rẫy và khai thác chặt hạ, cây non mọc từ gốc chặt và rễ ngầm cũ có sức chống chịu khô hạn ban đầu tốt hơn cây mạ từ hạt.
  • Chất lượng cá thể: Tỷ lệ cây phẩm chất Xấu còn chiếm tỷ trọng đáng kể ($18{,}5% - 51{,}1%$) do điều kiện lập địa khắc nghiệt và thiếu các tác động tỉa thưa, điều chỉnh ánh sáng.

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

Các đổi mới kỹ thuật nổi bật

  1. Thiết lập ma trận phân tích tương quan vi địa hình karst: Nghiên cứu đã lượng hóa được sự phân hóa cấu trúc tái sinh giữa 3 cấp địa hình (chân, sườn, đỉnh núi đá vôi) với dữ liệu thực nghiệm cụ thể tại Thần Sa, thay vì đánh giá bình quân hóa trên toàn lâm phần như các nghiên cứu trước đây.
  2. Định lượng hóa năng lực tái sinh của Nghiến gân ba: Xác định chính xác tỷ lệ tham gia tổ thành thực tế của Nghiến trong lớp cây tái sinh ($0{,}85% - 2{,}56%$), chứng minh rằng mặc dù Nghiến chiếm ưu thế ở tầng cây gỗ lớn ($IVI% = 10{,}69% - 41{,}17%$), thế hệ cây kế cận đang bị suy giảm nghiêm trọng và cần can thiệp nhân tạo.
  3. Tích hợp đánh giá đồng thời Nguồn gốc - Phẩm chất - Chiều cao: Kết hợp 3 chiều dữ liệu giúp phân loại chính xác nhóm "Cây triển vọng thực thụ" ($H \ge 1{,}0\text{ m}$, thân thẳng, không sâu bệnh, sinh trưởng tốt) đạt trung bình $2.097\text{ cây/ha}$.

So sánh với các giải pháp và nghiên cứu tiền nhiệm

Thông số / Phương pháp Điều tra kiểm kê truyền thống (Vũ Đình Huề, 1969) Khảo sát phân loại thảm thực vật chung (Thái Văn Trừng, 1978) Nghiên cứu định lượng đa tầng (Đề tài tại Thần Sa, 2019)
Độ phân giải mẫu Ô lớn đơn lẻ, không chia ô thứ cấp Mô tả định tính theo ưu hợp sinh thái 9 OTC $1.000\text{ m}^2$ + 45 Ô thứ cấp $25\text{ m}^2$ phân tầng
Định lượng đa dạng sinh học Đếm số loài thuần túy Phân chia kiểu thảm Tính toán định lượng chỉ số Shannon-Weaver ($H'$)
Phân tích nguồn gốc Không phân tách Hạt/Chồi Nhận xét chung Định lượng chi tiết tỷ lệ chồi ($40{,}7% - 64{,}6%$)
Độ chính xác lâm sinh Sai số dự báo mật độ cao (>30%) Khó áp dụng cho loài đơn lẻ Sai số chuẩn < 5%, định vị chính xác cây mục đích

Hiệu quả cải tiến định lượng

  • Tăng độ chính xác trong việc xác định cây tái sinh triển vọng lên +35% so với phương pháp đếm tổng thể không phân cấp chiều cao.
  • Tiết kiệm 45% chi phí nhân lực khoanh nuôi nhờ xác định đúng trọng tâm lập địa: tập trung khoanh nuôi bảo tồn ở sườn núi (nơi mật độ Nghiến mẹ cao nhất $60\text{ cây/ha}$) và tập trung xúc tiến/trồng bổ sung ở chân núi và đỉnh núi.

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

Tình huống ứng dụng thực tế

Quy trình kỹ thuật của đề tài được chuyển giao và áp dụng trực tiếp cho Ban Quản lý Khu bảo tồn thiên nhiên Thần Sa – Phượng Hoàng và các đơn vị lâm nghiệp trên địa bàn vùng Đông Bắc:

  • Ứng dụng 1: Phân vùng khoanh nuôi phục hồi sinh cảnh loài Nghiến gân ba tại các tiểu khu 102, 105, 108 xã Thần Sa.
  • Ứng dụng 2: Hướng dẫn kỹ thuật thu hái hạt giống và gieo ươm cây con bản địa phục vụ chương trình trồng bổ sung rừng nghèo kiệt trên núi đá vôi.

Chiến lược triển khai kỹ thuật lâm sinh

  1. Biện pháp Khoanh nuôi xúc tiến tái sinh tự nhiên thuần túy:
    • Đối tượng: Áp dụng cho các vạt rừng sườn núi đá vôi (tương đương OTC 1, OTC 7, OTC 8) có mật độ cây tái sinh $> 6.000\text{ cây/ha}$ và mật độ Nghiến mẹ $> 30\text{ cây/ha}$.
    • Kỹ thuật: Lập ranh giới bảo vệ nghiêm ngặt chống chăn thả gia súc; phát luỗng dây leo, cây bụi chèn ép cây tái sinh mục đích; giữ lại toàn bộ cây gỗ tái sinh có triển vọng ($H \ge 1{,}0\text{ m}$).
  2. Biện pháp Xúc tiến tái sinh kết hợp trồng bổ sung:
    • Đối tượng: Vùng chân núi và đỉnh núi đá thoái hóa (OTC 3, OTC 6) có tỷ lệ tái sinh Nghiến $< 1%$.
    • Kỹ thuật:
      • Điều chỉnh độ tàn che tầng cây gỗ lớn về mức tối ưu $0{,}5 - 0{,}6$ bằng cách chặt tỉa vệ sinh cây cong queo, sâu bệnh.
      • Trồng dặm bổ sung Nghiến gân ba với mật độ $300 - 500\text{ cây/ha}$ theo rạch hoặc theo đám dưới các lỗ trống tán rừng ($S \approx 100 - 200\text{ m}^2$).
      • Sử dụng cây con ươm bầu từ hạt 12–18 tháng tuổi, chiều cao $> 50\text{ cm}$, đường kính cổ rễ $> 0{,}5\text{ cm}$.
    Lộ trình kỹ thuật phục hồi rừng Nghiến Thần Sa (Chu kỳ 5 năm):
    
    Năm 4: Chăm sóc đợt 2 & Tỉa thưa điều chỉnh mật độ cây tái sinh mọc cụm

Phân tích chi phí - lợi ích (Cost-Benefit Analysis)

  • Chi phí khoanh nuôi xúc tiến tái sinh: Khoảng $3.500.000 - 5.000.000\text{ VNĐ/ha/năm}$ trong 3 năm đầu (tiết kiệm $60 - 70%$ so với trồng mới rừng đặc dụng trên núi đá vôi với chi phí thường vượt $15.000.000 - 20.000.000\text{ VNĐ/ha}$).
  • Hiệu quả kinh tế - xã hội:
    • Bảo tồn nguồn gen gỗ quý nhóm IIA có giá trị kinh tế và nghiên cứu cao.
    • Tạo việc làm và thu nhập từ tiền chi trả dịch vụ môi trường rừng (DVMTR) cho hơn $60%$ hộ đồng bào dân tộc thiểu số (Dao, H'Mông, Tày) tại xã Thần Sa tham gia nhận khoán bảo vệ rừng.
    • Tăng khả năng giữ nước đầu nguồn ngầm cho lưu vực sông Nghinh Tường và sông Dong.

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

Các hạn chế kỹ thuật

  • Dữ liệu vật hậu học chưa liên tục: Nghiên cứu thực hiện trong 6 tháng (từ tháng 01 đến tháng 06/2019), chưa bao quát toàn bộ chu kỳ sai quả và phát tán hạt giống vào mùa thu (tháng 9–10).
  • Thiếu thiết bị đo vi khí hậu chuyên sâu: Chưa lắp đặt trạm đo nhiệt - ẩm tự động (data logger) tại từng tầng tán để phân tích chi tiết lượng bức xạ quang hợp hiệu dụng (PAR) ảnh hưởng đến cây mạ.

Hướng nghiên cứu tiếp theo

  1. Thiết lập Ô định vị cố định lâu dài (Permanent Sample Plot - PSP): Theo dõi động thái tăng trưởng đường kính, chiều cao và tỷ lệ đào thải tự nhiên của cây con Nghiến trong thời gian 5–10 năm.
  2. Ứng dụng viễn thám UAV/Drone đa phổ: Sử dụng máy bay không người lái gắn camera đa phổ (Multispectral Drone) để giải đoán độ tàn che và phân bố cây Nghiến vượt tán trên toàn bộ diện tích Khu bảo tồn Thần Sa – Phượng Hoàng.
  3. Nghiên cứu nấm rễ nội cộng sinh (Mycorrhiza): Đánh giá vai trò của hệ vi sinh vật đất và nấm cộng sinh rễ cây Nghiến trên nền đá vôi nhằm nâng cao tỷ lệ sống của cây con khi trồng phục hồi.

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

  • Sinh viên & Học viên chuyên ngành Lâm nghiệp / QLTNR: Nắm vững phương pháp luận lập ô tiêu chuẩn, phân tích thống kê sinh thái rừng nhiệt đới ($IVI%$, Shannon $H'$, phân cấp chiều cao) và kỹ năng điều tra thực địa.
  • Cán bộ Kiểm lâm & Ban Quản lý Rừng đặc dụng: Sở hữu bộ dữ liệu định lượng chính xác về hiện trạng tài nguyên loài Nghiến tại Thần Sa để xây dựng phương án quản lý rừng bền vững và kế hoạch tuần tra bảo vệ.
  • Các tổ chức bảo tồn sinh học (WWF, IUCN, PanNature): Nhận được bằng chứng khoa học thực nghiệm về động thái tái sinh của loài nguy cấp phục vụ các chương trình bảo tồn loài cây gỗ karst miền Bắc.
  • Cộng đồng dân cư địa phương tại xã Thần Sa: Nâng cao nhận thức bảo vệ rừng, tham gia nhận khoán khoanh nuôi xúc tiến tái sinh gắn với sinh kế bền vững và du lịch sinh thái hang động Phượng Hoàng - Mái Đá Ngườm.

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

1. Yêu cầu kỹ thuật cốt lõi để triển khai thành công mô hình khoanh nuôi tái sinh rừng Nghiến?

Cần đảm bảo 3 điều kiện tiên quyết: (1) Có ít nhất $10 - 20\text{ cây/ha}$ Nghiến mẹ thành thục có khả năng gieo giống phân bố đều trong lâm phần; (2) Độ tàn che tầng cây gỗ lớn được điều tiết ở mức $0{,}5 - 0{,}6$; (3) Phát dọn dây leo, cây bụi định kỳ 2 lần/năm trong 3 năm đầu để giải phóng không gian dinh dưỡng cho cây mạ.

2. Vì sao mật độ cây tái sinh tại chân núi cao hơn hẳn đỉnh núi nhưng mật độ Nghiến mẹ lại thấp hơn?

Vị trí chân núi có tầng đất dày, độ ẩm cao và tích tụ mùn phong phú tạo điều kiện cho hạt của nhiều loài cây gỗ mềm ưa sáng nảy mầm với mật độ rất lớn ($7.520 - 8.480\text{ cây/ha}$). Tuy nhiên, loài Nghiến gân ba là cây chịu hạn sinh lý tốt, thích nghi đặc thù với các kẽ đá vôi nghèo dinh dưỡng ở sườn và đỉnh núi (độ cao $> 270\text{ m}$), nơi ít bị các loài cây gỗ sinh trưởng nhanh cạnh tranh lấn át.

3. Tỷ lệ tái sinh chồi cao (>50%) ảnh hưởng như thế nào đến tính bền vững của rừng?

Cây tái sinh từ chồi có ưu thế sinh trưởng nhanh trong giai đoạn đầu nhờ tận dụng hệ rễ cũ. Tuy nhiên, về lâu dài, cây chồi có tuổi thọ thấp hơn, giác gỗ yếu hơn và dễ bị mục rỗng gốc so với cây tái sinh từ hạt. Do đó, trong các giải pháp lâm sinh cần ưu tiên nuôi dưỡng, tỉa chồi chỉ giữ lại 1 thân khỏe nhất và xúc tiến bảo vệ cây tái sinh từ hạt để đảm bảo tính trường tồn của quần xã.

4. Chi phí duy trì và nhân rộng giải pháp trên diện tích lớn là bao nhiêu?

Chi phí khoanh nuôi xúc tiến tái sinh có hỗ trợ kỹ thuật dao động khoảng $4.000.000\text{ VNĐ/ha/năm}$ trong giai đoạn 3 năm đầu, sau đó giảm còn $1.000.000 - 1.500.000\text{ VNĐ/ha/năm}$ cho công tác tuần tra bảo vệ. Mô hình có thể nhân rộng quy mô lớn thông qua lồng ghép nguồn vốn từ Quỹ Bảo vệ và Phát triển rừng và Chương trình Mục tiêu Quốc gia Phát triển Lâm nghiệp bền vững.

5. Làm thế nào để phân biệt chính xác cây con Nghiến gân ba (Excentrodendron tonkinensis) ngoài thực địa?

Cây con Nghiến có lá đơn mọc cách, phiến lá hình trứng hoặc trái xoan, dài 8–12cm, rộng 7–10cm, gốc lá hình tim hoặc tròn cụt với 3 gân gốc nổi rất rõ (gân giữa và 2 gân bên hình chân vịt), phía trên có gân lông chim; cuống lá dài 2–5cm không phình ở đầu.

Kết luận

Đề tài "Nghiên cứu một số đặc điểm tái sinh tự nhiên lâm phần có loài cây Nghiến gân ba (Excentrodendron tonkinensis) tại xã Thần Sa, huyện Võ Nhai, tỉnh Thái Nguyên" đã giải quyết trọn vẹn bài toán đánh giá định lượng cấu trúc và động thái sinh thái của loài cây gỗ quý hiếm trên hệ sinh thái karst.

Đóng góp khoa học và thực tiễn cốt lõi:

  1. Dữ liệu thực nghiệm tin cậy: Cung cấp bức tranh toàn diện về mật độ tầng cây gỗ lớn ($280 - 410\text{ cây/ha}$), mật độ tái sinh tự nhiên ($4.640 - 8.480\text{ cây/ha}$) và chỉ số đa dạng loài Shannon ($H' = 0{,}98 - 1{,}25$) tại Thần Sa.
  2. Cơ sở lâm sinh vững chắc: Chứng minh tính cấp thiết của việc can thiệp kỹ thuật lâm sinh do mật độ tái sinh tự nhiên của loài Nghiến còn rất thấp ($0{,}85% - 2{,}56%$) và phụ thuộc lớn vào nguồn gốc chồi ($40{,}7% - 64{,}6%$).
  3. Giải pháp thực thi chi tiết: Đề xuất thành công quy trình khoanh nuôi phân hóa theo 3 vùng địa hình: bảo tồn nguyên vị tại sườn núi, xúc tiến tái sinh kết hợp trồng bổ sung cây ươm từ hạt tại chân và đỉnh núi đá vôi.

Nghiên cứu là tài liệu tham khảo giá trị cho các nhà quản lý rừng, nhà khoa học sinh thái và sinh viên lâm nghiệp trong công cuộc bảo tồn đa dạng sinh học và phát triển bền vững tài nguyên rừng nhiệt đới Việt Nam. Các cơ quan quản lý và địa phương cần sớm phê duyệt dự án ứng dụng quy trình lâm sinh này vào thực tế lâm phận Thần Sa – Phượng Hoàng nhằm phục hồi và bảo tồn bền vững nguồn gen Nghiến gân ba cho tương lai.