Giới thiệu dự án

Bối cảnh và thực trạng suy thoái tài nguyên rừng

Theo báo cáo của Tổ chức Bảo tồn Thiên nhiên Quốc tế (IUCN) và Quỹ Quốc tế Bảo vệ Thiên nhiên (WWF), độ che phủ rừng tự nhiên toàn cầu suy giảm gần 13% trong giai đoạn nửa cuối thế kỷ 20, trong đó rừng nhiệt đới Đông Nam Á chịu áp lực suy thoái lớn nhất do khai thác quá mức và chuyển đổi mục đích sử dụng đất. Tại Việt Nam, mặc dù tổng diện tích rừng tăng nhờ các chương trình trồng rừng quốc gia, diện tích rừng tự nhiên nguyên sinh chất lượng cao chỉ còn dưới 10%, phần lớn diện tích là rừng nghèo và rừng thứ sinh đang phục hồi (chiếm trên 66% tổng diện tích rừng cả nước).

Khu Bảo tồn Thiên nhiên (KBTTN) Thần Sa - Phượng Hoàng (huyện Võ Nhai, tỉnh Thái Nguyên) được thành lập theo Quyết định số 3841/QĐ-UB với diện tích 17.474,9 ha, sở hữu hệ sinh thái núi đá vôi đặc trưng cùng 1.096 loài thực vật bậc cao có mạch và nhiều loài động vật nguy cấp (như Voọc đen má trắng, Cu li lớn). Tuy nhiên, vùng lõi của Khu bảo tồn – đặc biệt là địa bàn xã Nghinh Tường – đang chịu áp lực nhân sinh khốc liệt làm suy thoái trữ lượng và đa dạng sinh thái.

+-----------------------------------------------------------------------------+
|                      KHU BẢO TỒN THẦN SA - PHƯỢNG HOÀNG                     |
|                                                                             |
|   +---------------------------------------------------------------------+   |
|   |                  VÙNG LÕI: XÃ NGHINH TƯỜNG (8.429,6 ha)             |   |
|   |                                                                     |   |
|   |  - Tổng diện tích rừng & đất LN: 8.243,3 ha (Độ che phủ 94%)       |   |
|   |  - Trạng thái IIa:   2.890,1 ha (35,1%) -> Rừng phục hồi nghèo      |   |
|   |  - Trạng thái IIb:   2.478,2 ha (30,1%) -> Rừng phục hồi trung bình |   |
|   |  - Trạng thái IIIa1: 1.498,4 ha (18,2%) -> Rừng thứ sinh giàu      |   |
|   |  - Đất LN & khác:    1.376,6 ha (17,6%)                             |   |
|   +---------------------------------------------------------------------+   |
+-----------------------------------------------------------------------------+

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

Nghinh Tường là xã vùng cao đặc biệt khó khăn với 12 xóm, trong đó có 5 xóm nằm hoàn toàn trong vùng lõi rừng đặc dụng (Nà Hẩu, Thượng Lương, Hạ Lương, Bản Rãi, Bản Nhàu). Cư dân bản địa (chủ yếu là dân tộc Tày 473 hộ, Dao 207 hộ, Nùng 26 hộ) phụ thuộc chặt chẽ vào tài nguyên rừng để duy trì sinh kế. Các điểm nghẽn kỹ thuật và quản trị gồm:

  • Thiếu hụt định lượng cấu trúc lâm phần và chỉ số ưu thế loài ($IVI$) theo từng trạng thái sinh thái.
  • Sự bất cập trong việc phân loại và định vị chính xác các nguy cơ suy thoái rừng theo chuẩn khung giám sát REDD+.
  • Áp lực kép giữa bảo tồn nghiêm ngặt vùng lõi và nhu cầu tài nguyên gỗ/củi gia dụng của cộng đồng bản địa chưa có mô hình sinh kế thay thế bền vững.

Mục tiêu nghiên cứu

  1. Xác định hiện trạng diện tích, cấu trúc tổ thành và trữ lượng của các trạng thái rừng tự nhiên ($IIb$, $IIIa1$) tại xã Nghinh Tường.
  2. Định lượng các chỉ số sinh thái học quần xã: Chỉ số giá trị quan trọng ($IVI$), chỉ số đa dạng Shannon-Wiener ($H'$), chỉ số ưu thế Simpson ($C_d$) và chỉ số tương đồng Sorensen ($SI$).
  3. Phân tích, đánh giá các tác động nhân sinh và tự nhiên đe dọa đa dạng sinh học và trữ lượng carbon rừng.
  4. Đề xuất hệ thống giải pháp kỹ thuật lâm sinh, quản lý cộng đồng và chuyển đổi sinh kế có tính khả thi cao.

Phương cận giải pháp và kết quả kỳ vọng

Đề tài kết hợp phương pháp nghiên cứu lâm học thực địa (Hệ thống ô tiêu chuẩn lồng nhau - Nested Plots), kỹ thuật phân tích không gian Hệ thông tin Địa lý (GIS) trên phần mềm MapInfo Pro, và điều tra xã hội học nông thôn có sự tham gia của người dân (PRA - Participatory Rural Appraisal). Kết quả cung cấp bộ cơ sở dữ liệu số hóa hoàn chỉnh về 58 loài cây gỗ, xác lập công thức tổ thành sinh thái và bản đồ phân cấp nguy cơ suy thoái phục vụ công tác điều hành của Ban Quản lý KBTTN Thần Sa - Phượng Hoàng.

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

  • Phạm vi không gian: Vùng lõi và vùng đệm xã Nghinh Tường, huyện Võ Nhai, tỉnh Thái Nguyên.
  • Phạm vi học thuật: Tập trung khảo sát 6 ô tiêu chuẩn đại diện ($OTC$) tại 2 trạng thái rừng có giá trị sinh thái cao ($IIb$ và $IIIa1$); khảo sát kinh tế - xã hội trên 20 hộ gia đình tiêu biểu thuộc các xóm trọng điểm.
  • Giới hạn kỹ thuật: Do địa hình núi đá vôi hiểm trở, nghiên cứu giới hạn số lượng ô tiêu chuẩn ngoại nghiệp ở 6 vị trí chân, sườn và đỉnh núi đại diện.

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

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

Bảng so sánh các phương pháp đánh giá suy thoái rừng hiện hành

Phương pháp đánh giá Ưu điểm Nhược điểm / Hạn chế Mức độ phù hợp với KBTTN Thần Sa
Giám sát viễn thám thuần túy (Remote Sensing) Bao phủ diện tích lớn, cập nhật chu kỳ nhanh qua ảnh Landsat/Sentinel. Không phát hiện được suy thoái tầng dưới tán (degradation under canopy), sai số lớn ở núi đá vôi. Trung bình (cần kết hợp thực địa).
Kiểm kê lâm nghiệp truyền thống Số liệu đo đếm cây đứng chi tiết, độ tin cậy cơ học cao. Tốn kém nhân lực, chu kỳ 5 năm/lần, thiếu lồng ghép yếu tố áp lực xã hội nhân sinh. Thấp (thiếu tính thích ứng nhanh).
Phương pháp tích hợp Đa nguồn (PRA + Nested OTC + GIS) (Đề xuất) Định lượng chính xác $IVI$, $H'$, $SI$; tích hợp căn nguyên kinh tế - xã hội; chi phí tối ưu. Đòi hỏi cán bộ am hiểu chuyên sâu cả lâm sinh học và kỹ năng điều tra cộng đồng. Rất cao (Tối ưu cho vùng lõi KBT).

Yêu cầu hệ thống quản lý dữ liệu lâm sinh (Phân loại MoSCoW)

  • Must have: Tính toán chính xác các chỉ số cấu trúc lâm phần ($G, D_{1.3}, H_{vn}, M$), bảng giá trị quan trọng ($IVI \ge 5%$), bảng chỉ số đa dạng loài ($H', C_d, SI$).
  • Should have: Bản đồ số hóa hiện trạng 3 trạng thái rừng ($IIa, IIb, IIIa1$) trên nền MapInfo v15.0; phân tích tương quan giữa khai thác gỗ củi và suy giảm mật độ lâm phần.
  • Could have: Ma trận lượng hóa rủi ro mất rừng theo 4 cấp độ (Cao, Trung bình, Thấp, Không đáng kể).
  • Won't have (lần này): Tự động hóa giám sát bằng thiết bị bay không người lái (UAV) thời gian thực và mô hình hóa trích xuất trữ lượng carbon vi sai $CO_2$ tương đương.

Thiết kế hệ thống và quy trình xử lý dữ liệu

Technology Stack & Công cụ đo lường

  • Phần mềm GIS & Bản đồ: MapInfo Professional v15.0 (Hệ quy chiếu VN-2000, Kinh tuyến trục $106^\circ 30'$, Múi chiếu $3^\circ$).
  • Thiết bị đo đạc thực địa: GPS cầm tay Garmin GPSMAP 64s (độ chính xác $\pm 3\text{m}$), Thước đo đường kính Richter $D_{1.3}$, Thước đo cao Suunto PM-5/1520, La bàn địa hình Harbin DQL-8.
  • Môi trường xử lý thống kê sinh thái: Python v3.10 (Thư viện pandas v2.1.0, numpy v1.24.3, scipy v1.11.0) và R v4.3.1 (Package vegan v2.6-4).

Phương pháp luận (Methodology)

        SƠ ĐỒ BỐ TRÍ Ô TIÊU CHUẨN LỒNG NHAU (NESTED PLOT)
+-------------------------------------------------------------+
| OTC Tổng thể (50m x 40m = 2.000 m2)                         |
| Đo đếm cây gỗ lớn, xác định ranh giới sinh thái             |
|                                                             |
|   +-------------------------------+                         |
|   | Ô thứ cấp (25m x 20m = 500 m2)|                         |
|   | Đo đếm cây gỗ D1.3 >= 5cm     |                         |
|   |                               |                         |
|   |   +-----------------------+   |                         |
|   |   | Ô cấp 1 (10m x 20m)   |   |                         |
|   |   | Đo cây đổ, gỗ mục     |   |                         |
|   |   +-----------------------+   |                         |
|   +-------------------------------+                         |
|                                                             |
| [o] Ô dạng bản nhỏ (2m x 2m): Đo đếm tái sinh & thảm tươi   |
+-------------------------------------------------------------+
  1. Phương pháp điều tra Ô tiêu chuẩn (OTC): Áp dụng phương pháp ô tiêu chuẩn lồng nhau theo quy chuẩn lâm học Kumiatun:
    • Ô tổng thể diện tích $2.000\text{ m}^2$ ($50\text{m} \times 40\text{m}$).
    • Ô thứ cấp diện tích $500\text{ m}^2$ ($25\text{m} \times 20\text{m}$) bên trong: Thu thập toàn bộ cây thân gỗ có đường kính ngang ngực $D_{1.3} \ge 5\text{ cm}$.
    • Ô cấp 1 diện tích $200\text{ m}^2$ ($10\text{m} \times 20\text{m}$): Đo đếm cây đổ gãy, gỗ chết mục.
    • Ô dạng bản nhỏ $4\text{ m}^2$ ($2\text{m} \times 2\text{m}$): Đo đếm mật độ thảm tươi, cây bụi tái sinh.
  2. Phương pháp điều tra PRA: Phỏng vấn cấu trúc và bán cấu trúc 20 hộ dân sinh sống tại các xóm Nà Hẩu, Thượng Lương, Hạ Lương; tham vấn trực tiếp Hạt Kiểm lâm Võ Nhai và UBND xã Nghinh Tường.

Implementation và kết quả

Thuật toán và công thức toán học lâm thái học

Để định lượng cấu trúc và tương quan sinh thái học giữa các loài thực vật, các thuật toán sau được triển khai trực tiếp:

  1. Tổng tiết diện ngang ($G$) và Trữ lượng thân cây đứng ($M$): $$G = \frac{\pi}{4 \times 10.000 \times S} \sum_{i=1}^{n} D_{1.3_i}^2 \quad (\text{m}^2/\text{ha})$$ $$M = G \times \bar{H}_{vn} \times f \quad (\text{m}^3/\text{ha}) \quad (\text{với hình số } f = 0,5)$$

  2. Chỉ số giá trị quan trọng ($IVI$): $$IVI_i (%) = \frac{A_i + D_i + RF_i}{3}$$ Trong đó:

    • Độ phong phú tương đối: $A_i (%) = \frac{n_i}{N} \times 100$ ($n_i$: số cá thể loài $i$, $N$: tổng số cá thể).
    • Độ ưu thế tương đối: $D_i (%) = \frac{g_i}{G_{total}} \times 100$ ($g_i$: tổng tiết diện ngang loài $i$).
    • Tần suất xuất hiện tương đối: $RF_i (%) = \frac{F_i}{\sum F} \times 100$ ($F_i = \frac{N_i}{N_{plots}}$).

  3. Chỉ số đa dạng loài Shannon-Wiener ($H'$) & Chỉ số ưu thế Simpson ($C_d$): $$H' = - \sum_{i=1}^{s} \left( \frac{n_i}{N} \right) \ln \left( \frac{n_i}{N} \right)$$ $$C_d = \sum_{i=1}^{s} \left( \frac{n_i}{N} \right)^2$$

  4. Chỉ số tương đồng Sorensen ($SI$): $$SI = \frac{2C}{A + B}$$ Trong đó: $A, B$ là số loài trong quần xã A và B; $C$ là số loài chung xuất hiện ở cả hai quần xã.

Đoạn mã Python chuẩn hóa tính toán tự động các chỉ số sinh thái

import numpy as np
import pandas as pd

def calculate_ecological_indices(df_trees, total_plot_area_ha=0.2):
    """
    Tính toán cấu trúc lâm phần: IVI, H', Cd, G (Tổng tiết diện ngang)
    df_trees columns: ['species', 'dbh_cm', 'h_m', 'plot_id']
    """
    # 1. Tiết diện ngang từng cây (g = pi * d^2 / 4) tính theo m2
    df_trees['ba_m2'] = (np.pi * (df_trees['dbh_cm'] / 100) ** 2) / 4
    total_ba = df_trees['ba_m2'].sum()
    total_trees = len(df_trees)
    num_plots = df_trees['plot_id'].nunique()
    
    # 2. Tổng hợp theo loài
    species_grp = df_trees.groupby('species').agg(
        n_ind=('dbh_cm', 'count'),
        species_ba=('ba_m2', 'sum'),
        plot_occ=('plot_id', 'nunique')
    ).reset_index()
    
    # 3. Tính toán Ai, Di, RFi, IVI
    species_grp['Ai'] = (species_grp['n_ind'] / total_trees) * 100
    species_grp['Di'] = (species_grp['species_ba'] / total_ba) * 100
    species_grp['Fi'] = (species_grp['plot_occ'] / num_plots) * 100
    species_grp['RFi'] = (species_grp['Fi'] / species_grp['Fi'].sum()) * 100
    species_grp['IVI'] = (species_grp['Ai'] + species_grp['Di'] + species_grp['RFi']) / 3
    
    # 4. Chỉ số Đa dạng loài Shannon-Wiener (H') và Simpson (Cd)
    pi = species_grp['n_ind'] / total_trees
    h_prime = -np.sum(pi * np.log(pi))
    cd_simpson = np.sum(pi ** 2)
    
    # 5. Cấu trúc bình quân lâm phần
    g_m2_ha = total_ba / total_plot_area_ha
    d_mean = df_trees['dbh_cm'].mean()
    h_mean = df_trees['h_m'].mean()
    volume_m3_ha = g_m2_ha * h_mean * 0.5 # f = 0.5
    
    return {
        'species_table': species_grp.sort_values(by='IVI', ascending=False),
        'H_prime': round(h_prime, 4),
        'Cd_Simpson': round(cd_simpson, 4),
        'G_m2_ha': round(g_m2_ha, 2),
        'M_m3_ha': round(volume_m3_ha, 2),
        'D_mean': round(d_mean, 2),
        'H_mean': round(h_mean, 2)
    }

Kết quả cấu trúc và đa dạng sinh thái lâm phần

Bảng cấu trúc lâm học tại 6 ô tiêu chuẩn điều tra thực địa

Trạng thái Ký hiệu OTC Vị trí xóm Mật độ (cây/ha) $\bar{D}_{1.3}$ (cm) $\bar{H}_{vn}$ (m) Tổng tiết diện ngang $G$ ($\text{m}^2/\text{ha}$) Trữ lượng $M$ ($\text{m}^3/\text{ha}$) Chỉ số Simpson ($C_d$) Số lượng loài ghi nhận ($s$)
IIb OTC 01 Nà Hẩu 405 11,13 12,45 4,21 22,81 0,070 15
IIb OTC 02 Thượng Lương 315 10,82 12,04 3,23 17,37 0,074 17
IIb OTC 03 Nà Hẩu 435 11,07 13,51 4,40 25,04 0,073 18
IIIa1 OTC 04 Thượng Lương 450 16,09 15,87 8,64 53,63 0,063 19
IIIa1 OTC 05 Hạ Lương 470 17,59 13,57 9,51 57,37 0,084 18
IIIa1 OTC 06 Hạ Lương 470 16,09 12,26 8,47 49,33 0,080 18

Phân tích tổ thành loài ưu thế ($IVI \ge 5%$)

  • Trạng thái IIb (Rừng phục hồi): Tổng số ghi nhận 29 loài. Tổ thành loài ưu thế đặc trưng bởi các loài cây tiên phong ưa sáng và trung tính như: Sến ($IVI = 9,88%$), Trẩu ($9,00%$), Kẹn ($9,89%$), Kháo ($8,66%$), Lim ($8,73%$), Cọc rau ($8,11%$), Vàng trứng ($8,05%$), Phay ($8,25%$).
  • Trạng thái IIIa1 (Rừng thứ sinh giàu): Tổng số ghi nhận 46 loài (17 loài xuất hiện chung với IIb). Chiếm ưu thế bởi các loài cây gỗ lớn bản địa thuộc tầng ưu dưỡng: Đinh ($IVI = 11,26%$), Kháo gù ($10,95%$), Vàng anh ($9,38%$), Dẻ ($8,93%$), Roi rừng ($8,93%$), Phay ($8,37%$).

Ma trận tương đồng thành phần loài Sorensen ($SI$)

             [OTC 01]   [OTC 02]   [OTC 03]   [OTC 04]   [OTC 05]   [OTC 06]
[OTC 01]       1.00       0.38       0.18       0.22       0.05       0.16
[OTC 02]       0.38       1.00       0.31       0.29       0.12       0.28
[OTC 03]       0.18       0.31       1.00       0.22       0.22       0.27
[OTC 04]       0.22       0.29       0.22       1.00       0.11       0.25
[OTC 05]       0.05       0.12       0.22       0.11       1.00       0.14
[OTC 06]       0.16       0.28       0.27       0.25       0.14       1.00

Nhận xét: Hệ số tương đồng giữa các OTC trong cùng một trạng thái (như OTC 01 và OTC 02: $SI = 0,38$) cao hơn đáng kể so với giữa các trạng thái khác nhau (như OTC 01 và OTC 05: $SI = 0,05$), phản ánh sự phân hóa rõ rệt về điều kiện vi lập địa và lịch sử tác động nhân sinh.

Phân tích định lượng các nguy cơ gây suy thoái rừng

Bảng cấu trúc áp lực nhân sinh và tự nhiên đối với rừng Nghinh Tường

Nhóm tác nhân Chi tiết hoạt động Tỷ lệ tác động (%) Mức độ rủi ro Hiện trạng và Hệ quả cụ thể
Nhân sinh (90%) Khai thác gỗ củi làm chất đốt 40% Cực kỳ cao 100% hộ dân sử dụng củi đun nấu hàng ngày; cạn kiệt tầng cây bụi và gỗ tái sinh.
Phát nương làm rẫy, canh tác dốc 15% Cao Chuyển đổi đất lâm nghiệp sang trồng ngô, sắn (diện tích gieo trồng tăng từ 150 lên 232 ha).
Khai thác gỗ làm nhà & chuồng trại 10% Cao 90% hộ khai thác chọn lọc các loài gỗ quý nhóm II-IV: Táu (50%), Sến (35%), Vàng tâm (35%), Kháo (30%), Giổi (25%), Đinh (20%), Lát (15%), Nghiến (10%).
Săn bắt động vật hoang dã 10% Cao Đặt bẫy, săn bắn làm suy giảm các loài phát tán hạt giống tự nhiên.
Khai thác khoáng sản & vật liệu 10% Trung bình Khai thác quặng, đá xây dựng quy mô nhỏ gây xói mòn sườn dốc.
Định cư và xây dựng công trình 10% Trung bình Xâm lấn đất rừng phục vụ giãn dân và mở rộng cơ sở hạ tầng.
Thu hái Lâm sản ngoài gỗ (LSNG) 5% Thấp Khai thác măng, mộc nhĩ, thảo dược vào mùa sinh trưởng.
Tự nhiên (10%) Biến đổi khí hậu & Rét đậm rét hại 5% Trung bình Nhiệt độ mùa đông hạ sâu xuống $3^\circ\text{C}$, khô hạn kéo dài gây nguy cơ cháy rừng.
Thực vật ngoại lai xâm lấn 5% Thấp Sự phát triển của cỏ lào, cỏ tranh cản trở tái sinh tự nhiên sau nương rẫy.

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

Đóng góp học thuật và kỹ thuật mới

  1. Thiết lập chuẩn hóa bộ chỉ số định lượng lâm học vùng núi đá vôi Thần Sa: Cung cấp công thức tổ thành tầng cây gỗ chính xác cho 2 trạng thái $IIb$ và $IIIa1$, minh chứng rằng rừng $IIIa1$ có trữ lượng trung bình ($53,44\text{ m}^3/\text{ha}$) cao gấp $2,46$ lần rừng $IIb$ ($21,74\text{ m}^3/\text{ha}$).
  2. Tích hợp khung đánh giá REDD+ với phương pháp PRA cấp cơ sở: Lượng hóa chính xác mối liên hệ giữa áp lực sinh kế (100% dùng củi, 90% dùng gỗ xây dựng) với mức độ suy thoái thể tích gỗ đứng của các loài quý hiếm ($IVI$ của Đinh, Nghiến, Lát hoa bị suy giảm nghiêm trọng ở khu vực gần khu dân cư).

So sánh với các mô hình quản lý rừng hiện nay

Tiêu chí Mô hình Quản lý Hành chính Thuần túy Mô hình Bảo tồn Đóng kín (Strict Isolation) Mô hình Quản lý Tích hợp Dựa vào Cộng đồng (Đề xuất)
Cơ sở dữ liệu Thống kê diện tích danh nghĩa. Dữ liệu đa dạng sinh học định kỳ 5-10 năm. Số hóa GIS + Đo đếm thực địa + Ma trận áp lực PRA.
Tham gia của người dân Thụ động, nhận khoán bảo vệ đơn thuần. Xung đột lợi ích, cấm đoán tuyệt đối. Đồng quản lý (Co-management), nâng cao năng lực sinh kế.
Hiệu quả ngăn chặn suy thoái Thấp (vẫn xảy ra khai thác trộm). Trung bình (chi phí tuần tra quá cao). Cao (Giảm 60-70% áp lực khai thác củi gỗ nhờ giải pháp thay thế).

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

Giải pháp kỹ thuật và quản lý thực địa

1. Nhóm giải pháp kỹ thuật lâm sinh

  • Khoanh nuôi xúc tiến tái sinh có trồng bổ sung: Áp dụng cho $2.890,1\text{ ha}$ rừng trạng thái $IIa$ và các diện tích nương rẫy bỏ hoang. Trồng bổ sung các loài cây bản địa có giá trị sinh thái và kinh tế cao: Lát hoa (Chukrasia tabularis), Lim xanh (Erythrophleum fordii), Vàng anh (Saraca dives), Trám trắng (Canarium album).
  • Bảo tồn nguyên vị (In-situ Conservation): Khoanh bảo vệ nghiêm ngặt các dải rừng kín thường xanh núi đá vôi tại xóm Hạ Lương, Thượng Lương – nơi lưu trữ các cá thể cây họ Dầu, Re và Nghiến cổ thụ.

2. Nhóm giải pháp chuyển đổi năng lượng và sinh kế

  • Mô hình Hầm khí sinh học Biogas & Bếp đun cải tiến: Giảm ngay $60%$ nhu cầu sử dụng củi đun trong 207 hộ dân người Dao và 473 hộ người Tày.
  • Quy hoạch vùng trồng cây lấy củi và gỗ gia dụng nhanh: Trồng phân tán Keo lai, Bạch đàn, Xoan ta ở đất vườn đồi để giải tỏa áp lực đốn hạ Táu, Sến, Đinh trong vùng lõi.
  • Phát triển kinh tế dưới tán rừng: Trồng cây dược liệu bản địa (Hà thủ ô, Ba kích, Sa nhân tím) và chăn nuôi ong mật hoa rừng tự nhiên.
+-----------------------------------------------------------------------------------+
|               LỘ TRÌNH TRIỂN KHAI GIẢI PHÁP TẠI XÃ NGHINH TƯỜNG                  |
|                                                                                   |
|  [Giai đoạn 1: 0 - 6 tháng]                                                       |
|  * Thành lập 5 Tổ tuần tra rừng cộng đồng tự quản tại 5 xóm vùng lõi.            |
|  * Cắm mốc ranh giới số hóa GIS giữa vùng bảo vệ nghiêm ngặt và vùng phục hồi.    |
|                                                                                   |
|  [Giai đoạn 2: 6 - 18 tháng]                                                      |
|  * Triển khai 100 mô hình bếp tiết kiệm củi & hầm Biogas điểm.                   |
|  * Trồng bổ sung 50 ha cây bản địa (Lát hoa, Lim, Trám) tại trạng thái IIa.       |
|                                                                                   |
|  [Giai đoạn 3: 18 - 36 tháng]                                                     |
|  * Đánh giá lại chỉ số IVI, H' chu kỳ 3 năm trên 6 OTC chuẩn hóa.                 |
|  * Nhân rộng mô hình nông lâm kết hợp & sinh kế du lịch sinh thái khảo cổ học.     |
+-----------------------------------------------------------------------------------+

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

Hạn chế kỹ thuật

  • Số lượng ô tiêu chuẩn ($N=6\text{ OTC}$) còn hạn chế do ngân sách và địa hình núi karst hiểm trở, chưa bao phủ hết toàn bộ diện tích $8.243,3\text{ ha}$ đất lâm nghiệp.
  • Chưa thực hiện giải trình tự gen hoặc phân tích vi sinh vật đất để đánh giá toàn diện chu trình dinh dưỡng trong các hệ sinh thái rừng suy thoái.

Hướng nghiên cứu mở rộng

  • Ứng dụng thuật toán học máy (Random Forest / Support Vector Machine) trên ảnh viễn thám độ phân giải cao (Sentinel-2, PlanetScope) để tự động hóa phân loại suy thoái tán rừng.
  • Lượng hóa sinh khối trên mặt đất (AGB - Aboveground Biomass) và tính toán trữ lượng carbon hấp thụ để xây dựng hồ sơ tín chỉ carbon theo tiêu chuẩn Verra / Gold Standard.

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

  • Sinh viên & Học viên Lâm nghiệp / Môi trường: Cung cấp tài liệu tham khảo thực tiễn với đầy đủ quy trình thiết lập OTC lồng nhau, công thức toán học tính $IVI, H', C_d, SI$ và mã nguồn tính toán tự động.
  • Kỹ sư GIS & Quản lý Tài nguyên Rừng: Bản thiết kế kiến trúc phân tích cảnh quan không gian kết hợp số liệu điều tra kinh tế - xã hội PRA.
  • Ban Quản lý Khu Bảo tồn & Hạt Kiểm lâm: Bộ số liệu khoa học chính xác về cấu trúc lâm phần và ma trận nguyên nhân suy thoái, trực tiếp phục vụ công tác quy hoạch phân khu chức năng Thần Sa - Phượng Hoàng.
  • Chính quyền và Cộng đồng Dân cư Bản địa: Định hướng chuyển đổi sinh kế bền vững, giảm thiểu rủi ro thiên tai lũ quét, bảo vệ nguồn nước đầu nguồn.

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

1. Yêu cầu kỹ thuật cốt lõi để triển khai hệ thống ô tiêu chuẩn lồng nhau trên địa hình núi đá vôi?

Cần thiết lập ô tổng thể diện tích $2.000\text{ m}^2$ ($50\text{m} \times 40\text{m}$) căn chỉnh theo đường đồng mức bằng la bàn địa hình, các ô thứ cấp $500\text{ m}^2$ và ô cấp một $200\text{ m}^2$ được đặt lồng ở góc chuẩn. Toàn bộ cây $D_{1.3} \ge 5\text{ cm}$ phải được đánh số cố định bằng sơn chịu nhiệt và định vị tọa độ tâm ô bằng thiết bị GPS có sai số dưới $3\text{m}$.

2. Tại sao chỉ số $IVI \ge 5%$ lại được chọn làm ngưỡng xác định loài ưu thế sinh thái?

Theo quy chuẩn phân tích lâm sinh học của Thái Văn Trừng và Curtis & McIntosh, các loài có giá trị quan trọng $IVI < 5%$ thường là các loài ngẫu nhiên hoặc loài đi kèm không đóng vai trò quyết định trong việc hình thành cấu trúc tầng tán và vi khí hậu lâm phần. Ngưỡng $\ge 5%$ giúp tập trung nguồn lực phục hồi vào các loài then chốt (keystone species).

3. Phương pháp tích hợp số liệu điều tra PRA với dữ liệu không gian GIS như thế nào?

Dữ liệu PRA (vị trí nương rẫy, tuyến đường khai thác gỗ trộm, điểm chăn thả gia súc) được gắn tọa độ thực địa bằng GPS cầm tay, sau đó xuất file .kml/.shp để nạp vào MapInfo Pro. Bản đồ không gian sẽ tiến hành chồng lớp (overlay analysis) với bản đồ hiện trạng rừng để khoanh vùng các điểm nóng suy thoái (degradation hotspots).

4. Chi phí đầu tư cho giải pháp giảm thiểu khai thác gỗ củi là bao nhiêu và thời gian hoàn vốn?

Chi phí lắp đặt một hầm Biogas composite quy mô hộ gia đình khoảng 10-12 triệu VNĐ hoặc bếp đun cải tiến khoảng 300.000 - 500.000 VNĐ. Thời gian hoàn vốn thông qua việc tiết kiệm thời gian thu nhặt củi (trung bình 30-40 ngày công/năm/hộ) và tận dụng phân bón hữu cơ là khoảng 1,5 đến 2 năm.

5. Dữ liệu của đề tài có khả năng mở rộng sang các xã lân cận trong KBTTN Thần Sa không?

Hoàn toàn khả thi. Do có cùng điều kiện địa chất karst và đặc điểm kinh tế - xã hội vùng cao, quy trình khảo sát OTC lồng nhau và mô hình thuật toán Python/R tính toán cấu trúc lâm phần có thể áp dụng nguyên vẹn cho các xã lân cận như Sảng Mộc, Vũ Chấn, Phú Thượng.

Kết luận

Khóa luận tốt nghiệp của tác giả Tráng A Lầu dưới sự hướng dẫn của TS. Đàm Văn Vinh đã hoàn thành xuất sắc các mục tiêu nghiên cứu đặt ra. Đề tài không chỉ làm sáng tỏ thực trạng phân bố $8.243,3\text{ ha}$ rừng tự nhiên tại xã Nghinh Tường mà còn định lượng chính xác sự vượt trội về cấu trúc và trữ lượng của trạng thái rừng $IIIa1$ ($53,44\text{ m}^3/\text{ha}$) so với trạng thái $IIb$ ($21,74\text{ m}^3/\text{ha}$). Quan trọng hơn, nghiên cứu đã chỉ rõ áp lực nhân sinh khốc liệt từ hoạt động sử dụng củi đun (100% số hộ) và khai thác gỗ làm nhà (90% số hộ) là nguyên nhân cốt lõi gây suy thoái tài nguyên rừng vùng lõi KBTTN Thần Sa - Phượng Hoàng.

Bộ giải pháp kỹ thuật lâm sinh kết hợp chuyển đổi sinh kế cộng đồng được đề xuất là cơ sở khoa học tin cậy, đóng góp thiết thực cho chiến lược bảo tồn đa dạng sinh học và phát triển lâm nghiệp bền vững tại tỉnh Thái Nguyên. Ban quản lý khu bảo tồn và chính quyền địa phương cần sớm phê duyệt và lồng ghép các kiến nghị của đề tài vào quy chế quản lý rừng đặc dụng trong giai đoạn tới.