Giới thiệu dự án

Rừng tự nhiên đóng vai trò thiết yếu trong việc bảo tồn đa dạng sinh học (ĐDSH), điều hòa thủy văn và hấp thụ khí nhà kính toàn cầu. Tại Việt Nam, diện tích rừng tự nhiên từng trải qua giai đoạn suy giảm nghiêm trọng từ 43% năm 1943 xuống còn 28,4% năm 1990 do áp lực chuyển đổi đất nông nghiệp, du canh du cư và khai thác gỗ quá mức. Mặc dù các chính sách giao đất giao rừng và đóng cửa rừng tự nhiên theo Chỉ thị 286/TTg đã nâng độ che phủ rừng toàn quốc lên 38,7% vào năm 2008 (tương đương 13,1 triệu ha), chất lượng rừng tự nhiên vẫn tiếp tục suy giảm. Rừng giàu và rừng kín chỉ chiếm 4,6% tổng diện tích, trong khi 2/3 diện tích còn lại là rừng nghèo, rừng non tái sinh hoặc rừng thứ sinh phân mảnh.

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 ngày 01/12/1999 với tổng diện tích 11.220 ha. Đây là khu vực đặc trưng cho hệ sinh thái rừng trên núi đá vôi độc đáo của vùng Đông Bắc, lưu trữ nguồn gen phong phú với 1.096 loài thực vật bậc cao có mạch cùng 212 loài động vật có xương sống, trong đó có nhiều loài nằm trong Sách Đỏ Việt Nam và Danh lục Đỏ IUCN. Tuy nhiên, vùng lõi của khu bảo tồn—đặc biệt là xã Thần Sa (chiếm trên 80% diện tích đất lâm nghiệp)—đang đối mặt với nguy cơ suy thoái cấu trúc lâm phần và suy giảm đa dạng sinh học do áp lực nhân sinh phức tạp.

+-----------------------------------------------------------------------------------+
|                  BỐI CẢNH SUY THOÁI RỪNG TẠI XÃ THẦN SA                           |
+-----------------------------------------------------------------------------------+
|  1. Hệ sinh thái núi đá vôi: Đất mỏng, tái sinh tự nhiên cực chậm                |
|  2. Kinh tế tự cấp: 100% hộ dân dùng củi, phụ thuộc khai thác lâm sản            |
|  3. Lâm sản quý bị tận diệt: Nghiến, Lát hoa, Trai lý suy giảm tầng cây ưu thế   |
|  4. Thiếu hụt công cụ định lượng: Chưa chuẩn hóa dữ liệu cấu trúc cho REDD+       |
+-----------------------------------------------------------------------------------+

Vấn đề nghiên cứu và Điểm nghẽn thực tế

Hệ sinh thái rừng trên núi đá vôi tại xã Thần Sa có tính dễ tổn thương sinh thái cao: tầng đất mặt mỏng, địa hình chia cắt mạnh mẽ (độ dốc bình quân > 35°), khả năng tái sinh tự nhiên sau tác động chặt phá diễn ra rất chậm. Các điểm nghẽn chính bao gồm:

  • Khai thác chọn lọc tận diệt: Các loài cây gỗ có giá trị thương phẩm cao như Nghiến (Burretiodendron hsienmu), Lát hoa (Chukrasia tabularis), Trai lý (Garcinia fagraeoides) bị đốn hạ triệt để bằng cưa máy cầm tay, làm phá vỡ tầng tán và mất nguồn hạt giống mẹ.
  • Áp lực sinh kế bản địa: Cộng đồng 6 dân tộc (Tày, Dao, H'Mông, Sán Dìu, Nùng, Kinh) với tập quán sử dụng gỗ làm nhà sàn, khai thác củi đun nấu 100% và chăn thả gia súc tự do trong rừng đặc dụng.
  • Thiếu hụt dữ liệu định lượng phục vụ cơ chế REDD+: Công tác quản lý bảo vệ rừng tại địa phương thiếu các chỉ số đo lường sinh thái chuẩn xác (Chỉ số giá trị quan trọng IVI, trữ lượng $M$, tiết diện ngang $G$, hệ số tương đồng $SI$) để xây dựng phương án chi trả dịch vụ môi trường rừng và giảm phát thải khí nhà kính từ suy thoái rừng (REDD+).

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

  1. Xác định hiện trạng phân bố diện tích và cấu trúc tầng cây gỗ: Thống kê định lượng 6 trạng thái thảm thực vật tự nhiên (IA, IB, IC, IIA, IIB, IIIA1) trên diện tích 9.537,30 ha đất lâm nghiệp tại xã Thần Sa.
  2. Định lượng đa dạng thực vật thân gỗ và chất lượng trữ lượng: Thiết lập hệ thống ô tiêu chuẩn (OTC) lồng nhau và ô vệ tinh để tính toán mật độ ($N$), đường kính ngang ngực ($D_{1.3}$), chiều cao vút ngọn ($H_{vn}$), tổng tiết diện ngang ($G$), trữ lượng gỗ ($M$), chỉ số giá trị quan trọng ($IVI$), chỉ số đa dạng Shannon-Wiener ($H'$) và chỉ số ưu thế Simpson ($Cd$).
  3. Phân tích nguy cơ suy thoái và lượng hóa áp lực nhân sinh: Sử dụng công cụ Đánh giá nhanh nông thôn có sự tham gia (PRA) để phân loại 7 nhóm tác động chính và xác định chuỗi nguyên nhân gốc rễ.
  4. Đề xuất hệ thống giải pháp quản lý đa tầng: Xây dựng mô hình phối hợp quản lý giữa Ban quản trị KBT, chính quyền xã và cộng đồng thôn bản gắn liền với cải thiện sinh kế.

Phương pháp tiếp cận và Phạm vi

  • Phương pháp tiếp cận: Kết hợp giữa điều tra trắc lượng lâm nghiệp dã ngoại bằng hệ thống ô tiêu chuẩn cải tiến (Zoehrer, 1980), phân tích không gian cảnh quan GIS (MapInfo Pro v15.0) và phân tích xã hội học môi trường (PRA).
  • Phạm vi không gian: Toàn bộ địa giới hành chính xã Thần Sa với diện tích tự nhiên 10.144 ha ($101,44\text{ km}^2$), tọa độ địa lý từ $105^\circ 51'12''$ đến $106^\circ 08'38''$ Độ Đông và $21^\circ 45'12''$ đến $21^\circ 56'30''$ Vĩ Độ Bắc.
  • Phạm vi đối tượng: Tầng cây gỗ có đường kính ngang ngực $D_{1.3} \ge 5\text{ cm}$ thuộc các trạng thái rừng phục hồi (IIB) và rừng trung bình/giàu nghèo kiệt (IIIA1) trên núi đá vôi, cùng mẫu điều tra xã hội học 21 hộ gia đình tại các thôn trọng điểm (Xuyên Sơn, Ngọc Sơn 1, Ngọc Sơn 2).

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

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

Phương pháp đánh giá suy thoái rừng truyền thống thường dựa vào phân loại giải đoán ảnh viễn thám thô hoặc báo cáo kiểm kê diện tích mang tính định tính, dẫn đến sai số lớn trong việc nhận diện suy thoái tầng tán rừng kín trên núi đá vôi.

Tiêu chí so sánh Phương pháp Viễn thám thuần túy Phương pháp Thống kê lâm bạ cũ Phương pháp Tích hợp OTC Lồng nhau + GIS + PRA (Đề tài)
Độ phân giải cấu trúc Thấp (Chỉ thấy độ che phủ tán) Trung bình (Số liệu tổng hợp vĩ mô) Rất cao (Đo đếm chi tiết từng cá thể $D_{1.3} \ge 5\text{ cm}$)
Đánh giá suy thoái chất lượng Kém (Không phân biệt được loài ưu thế) Yếu (Chỉ ước lượng trữ lượng bình quân) Chính xác qua chỉ số $IVI$, cấu trúc tầng tán và $Cd$
Phân tích áp lực nhân sinh Không xác định được nguyên nhân Định tính, chung chung Định lượng chi tiết qua 21 phiếu PRA và ma trận tác động
Khả năng ứng dụng cho REDD+ Cần kiểm chứng mặt đất phức tạp Thiếu độ tin cậy khoa học Đáp ứng chuẩn đo lường MRV (Đo đạc, Báo cáo, Thẩm định)

Ma trận Yêu cầu Kỹ thuật theo Khung MoSCoW

  • Must have (Bắt buộc): Xác định chính xác tọa độ GPS của mạng lưới OTC; hiệu chỉnh chiều dài cạnh ô theo độ dốc địa hình bằng thước đo độ dốc Suunto Clinometer; tính toán đầy đủ các chỉ số $N, G, M, IVI, H', Cd, SI$.
  • Should have (Nên có): Phân bố ô tiêu chuẩn vệ tinh 4 hướng đối với sườn núi đá vôi dốc đứng; phân tích tương quan giữa thành phần loài bị khai thác và thu nhập kinh tế hộ dân.
  • Could have (Có thể mở rộng): Xây dựng bản đồ phân vùng nhạy cảm suy thoái tỷ lệ 1:25.000 trên hệ thống MapInfo.
  • Won't have (Không thực hiện): Khảo sát chi tiết tầng cây tái sinh dưới $D_{1.3} < 5\text{ cm}$ và thảm tươi cây bụi trong các ô tiêu chuẩn thứ cấp $1\text{ m} \times 1\text{ m}$.

Thiết kế hệ thống điều tra và thu thập dữ liệu

Quy trình nghiên cứu được thiết kế theo cấu trúc module dữ liệu chuẩn hóa, tích hợp các công cụ phần cứng lâm sinh và phần mềm phân tích thống kê chuyên dụng.

                    SƠ ĐỒ KHUNG PHÂN TÍCH VÀ ĐIỀU TRA TÀI NGUYÊN
+-----------------------------------------------------------------------------------+
|                           1. NGOẠI NGHIỆP DÃ NGOẠI                                |
|  +------------------------------+     +----------------------------------------+  |
|  |     Định vị GPS & MapInfo    | --> | Thiết lập OTC Lồng nhau (50m x 40m)    |  |
|  | (Garmin 64s / VN-2000 / WGS) |     | Hiệu chỉnh góc dốc Suunto Clinometer   |  |
|  +------------------------------+     +----------------------------------------+  |
|                                                           |                       |
|                                                           v                       |
|  +------------------------------+     +----------------------------------------+  |
|  |     Phiếu PRA Xã hội học     |     | Đo đếm Trắc lượng Tầng Cây Gỗ          |  |
|  | (21 hộ dân / 7 nhóm nguy cơ) |     | (D1.3 >= 5cm, Hvn, Tên loài thực vật)  |  |
|  +------------------------------+     +----------------------------------------+  |
+-----------------------------------------------------------------------------------+
                                            |
                                            v
+-----------------------------------------------------------------------------------+
|                           2. NỘI NGHIỆP TÍNH TOÁN                                 |
|  +-----------------------------------------------------------------------------+  |
|  | Module R / MS Excel: Tính Tiết diện ngang G, Trữ lượng M, Chỉ số IVI, H', Cd |  |
|  | Module MapInfo Pro: Phân lớp không gian 6 trạng thái rừng (IA -> IIIA1)     |  |
|  | Module PRA Matrix: Xếp hạng nguyên nhân suy thoái & Đề xuất giải pháp       |  |
|  +-----------------------------------------------------------------------------+  |
+-----------------------------------------------------------------------------------+

Ngăn xếp Công nghệ và Thiết bị (Technology & Tool Stack)

  • Thiết bị đo trắc lượng dã ngoại:
    • Máy định vị vệ tinh cầm tay: Garmin GPSMAP 64s (Độ chính xác định vị sai số $< 3.5\text{ m}$).
    • Thước đo độ dốc quang học: Suunto Clinometer PM-5/360 PC (Độ chính xác $\pm 0.25^\circ$).
    • Dụng cụ trắc lượng thân cây: Thước dây thép 50m, thước kẹp kính lúp rèn sợi kim loại (đo đường kính $D_{1.3}$ sai số $0.1\text{ cm}$), thước đo cao bọt nước Blum-Leiss.
  • Phần mềm xử lý dữ liệu và Không gian:
    • Hệ thông tin địa lý (GIS): MapInfo Professional v15.0 (64-bit) (Số hóa ranh giới thảm phủ, biên tập bản đồ chuyên đề thảm thực vật tỷ lệ 1:25.000).
    • Phân tích thống kê sinh thái: R v3.3.0 Engine kết hợp Microsoft Excel 2016 (Xử lý ma trận tổ thành, tính toán $IVI$, tương đồng Sorensen).

Cấu trúc Dữ liệu Đo đếm Tầng Cây Gỗ (Data Schema)

Mỗi bản ghi cây gỗ trong ô tiêu chuẩn được chuẩn hóa theo định dạng bảng sau:

Table: Forest_Inventory_Records
-----------------------------------------------------------------
Field Name      | Data Type    | Unit      | Description
-----------------------------------------------------------------
Plot_ID         | String(10)   | None      | Mã ô tiêu chuẩn (vd: IIB_OTC1)
Subplot_Index   | Integer      | None      | Chỉ số ô thứ cấp (1 đến 4)
Tree_ID         | Integer      | None      | Số thứ tự cá thể cây trong OTC
Species_VN      | String(50)   | None      | Tên tiếng Việt (vd: Vàng anh)
Species_Sci     | String(100)  | None      | Tên khoa học (vd: Saraca dives)
DBH_1_3         | Float        | cm        | Đường kính ngang ngực tại 1.3m
Height_Top      | Float        | m         | Chiều cao vút ngọn (Hvn)
Basal_Area_g    | Float        | m^2       | Tiết diện thân cá thể (pi * D^2 / 40000)
Volume_v        | Float        | m^3       | Thể tích thân cây (g * Hvn * f)
Status_Health   | String(20)   | None      | Tình trạng sinh trưởng (Tốt/Cụt ngọn)
-----------------------------------------------------------------

Methodology

Quy trình phương pháp luận được thiết kế thành 3 giai đoạn chặt chẽ:

+------------------------------------------------------------------------------------+
|                           QUY TRÌNH PHƯƠNG PHÁP NGHIÊN CỨU                         |
+------------------------------------------------------------------------------------+
|  Giai đoạn 1: Chuẩn bị & Thu thập thứ cấp (Tháng 02/2016)                          |
|  - Kế thừa bản đồ địa chính, số liệu khí tượng thủy văn, báo cáo KBT Thần Sa       |
|  - Rút mẫu điều tra: Tỷ lệ rút mẫu 2% diện tích từng khối trạng thái rừng          |
+------------------------------------------------------------------------------------+
|  Giai đoạn 2: Ngoại nghiệp thực địa (Tháng 03 - Tháng 04/2016)                     |
|  - Thiết lập OTC tập trung 2.000 m2 (50m x 40m) & OTC vệ tinh (Zoehrer, 1980)      |
|  - Hiệu chỉnh góc dốc sườn núi đá: L' = L / cos(alpha)                             |
|  - Phỏng vấn PRA 21 hộ dân tại 3 xóm trọng điểm (Xuyên Sơn, Ngọc Sơn 1 & 2)        |
+------------------------------------------------------------------------------------+
|  Giai đoạn 3: Nội nghiệp & Xử lý thống kê (Tháng 05/2016)                          |
|  - Chỉnh lý số liệu, tính toán IVI, Shannon-Wiener (H'), Simpson (Cd), Sorensen    |
|  - Phân tích nguyên nhân suy thoái, xây dựng ma trận đánh giá áp lực nhân sinh     |
|  - Hoàn thiện báo cáo khoa học và khuyến nghị chính sách đồng quản lý               |
+------------------------------------------------------------------------------------+

Ma trận Quản lý Rủi ro Thực địa (Risk Assessment)

  • Rủi ro trượt ngã địa hình núi đá tai mèo: Độ dốc trung bình trên 35°, vách đá hiểm trở. Biện pháp giảm thiểu: Thay thế OTC tập trung bằng hệ thống 4 ô tiêu chuẩn vệ tinh $10\text{ m} \times 50\text{ m}$ theo các hướng chính hoặc hướng phụ (Đông Nam, Tây Nam) khi gặp địa hình chắn.
  • Sai số đo đạc chiều dài cạnh ô do độ dốc: Biện pháp giảm thiểu: Đo góc dốc $\alpha$ bằng máy nghiêng bọt nước Clinometer và tra bảng hiệu chỉnh cự ly nằm ngang theo công thức $L' = L / \cos(\alpha)$.
  • Sai số thiên vị trong phỏng vấn PRA: Người dân có tâm lý giấu hành vi khai thác gỗ lậu. Biện pháp giảm thiểu: Kết hợp phỏng vấn bán cấu trúc với trưởng thôn, kiểm lâm địa bàn và quan sát trực tiếp dấu vết chặt hạ tại hiện trường dã ngoại.

Implementation và kết quả

Development process & Thuật toán sinh thái định lượng

Toàn bộ dữ liệu trắc lượng thu thập từ 6 ô tiêu chuẩn đại diện (diện tích mỗi OTC = $2.000\text{ m}^2 = 0,2\text{ ha}$, tổng diện tích đo đếm trực tiếp $12.000\text{ m}^2$) thuộc 2 trạng thái IIB (Rừng phục hồi) và IIIA1 (Rừng nghèo kiệt sau khai thác) được xử lý qua hệ thống công thức toán học sinh thái:

1. Công thức Hiệu chỉnh Chiều dài Cạnh Ô theo Góc dốc

$$L' = \frac{L}{\cos(\alpha)}$$ Trong đó: $L'$ là chiều dài thực tế cần kéo thước trên sườn dốc ($\text{m}$); $L$ là chiều dài hình chiếu nằm ngang chuẩn ($50\text{ m}$ hoặc $40\text{ m}$); $\alpha$ là góc dốc địa hình đo bằng Clinometer ($^\circ$).

2. Công thức Chỉ số Giá trị Quan trọng (Importance Value Index - IVI)

$$IVI_i (%) = \frac{A_i + D_i + RF_i}{3} \quad \text{hoặc} \quad IVI_i = A_i + D_i + RF_i$$

  • Độ phong phú tương đối ($A_i$): $$A_i (%) = \left(\frac{n_i}{\sum_{j=1}^s n_j}\right) \times 100$$
  • Độ ưu thế tương đối ($D_i$): $$D_i (%) = \left(\frac{g_i}{\sum_{j=1}^s g_j}\right) \times 100 \quad \text{với} \quad g_i = \sum \frac{\pi \cdot D_{1.3}^2}{40.000}$$
  • Tần suất xuất hiện tương đối ($RF_i$): $$RF_i (%) = \left(\frac{F_i}{\sum_{j=1}^s F_j}\right) \times 100 \quad \text{với} \quad F_i = \left(\frac{N_i}{N_{\text{total_sub}}}\right) \times 100$$

3. Các Chỉ số Đa dạng và Tương đồng Sinh thái

  • Mật độ cá thể trên một hecta ($N$): $$N = \frac{n \times 10.000}{S} \quad (\text{cây/ha})$$
  • Chỉ số Đa dạng loài Shannon-Wiener ($H'$): $$H' = -\sum_{i=1}^s \left(P_i \ln P_i\right) \quad \text{với} \quad P_i = \frac{n_i}{N_{\text{total}}}$$
  • Chỉ số Tập trung Ưu thế Simpson ($Cd$): $$Cd = \sum_{i=1}^s \left(\frac{n_i}{N_{\text{total}}}\right)^2$$
  • Chỉ số Tương đồng Thành phần Loài Sorensen ($SI$): $$SI = \frac{2a}{2a + b + c}$$ Trong đó: $a$ là số loài xuất hiện ở cả hai trạng thái; $b$ là số loài duy nhất chỉ có ở trạng thái 1; $c$ là số loài duy nhất chỉ có ở trạng thái 2.
# Script R tính toán cấu trúc tổ thành và đa dạng sinh thái lâm phần
calculate_ecology_indices <- function(data_tree) {
  # data_tree: data.frame(Species, DBH, Subplot_ID)
  total_trees <- nrow(data_tree)
  total_g <- sum(pi * (data_tree$DBH^2) / 40000)
  
  species_summary <- aggregate(DBH ~ Species, data = data_tree, FUN = function(x) {
    c(n = length(x), g = sum(pi * (x^2) / 40000))
  })
  
  species_df <- data.frame(
    Species = species_summary$Species,
    n = species_summary$DBH[, "n"],
    g = species_summary$DBH[, "g"]
  )
  
  # Tính toán Ai và Di
  species_df$Ai <- (species_df$n / total_trees) * 100
  species_df$Di <- (species_df$g / total_g) * 100
  
  # Tính toán RFi dựa trên 4 ô thứ cấp
  n_subplots <- length(unique(data_tree$Subplot_ID))
  freq <- sapply(species_df$Species, function(sp) {
    length(unique(data_tree$Subplot_ID[data_tree$Species == sp])) / n_subplots
  })
  species_df$Fi <- freq * 100
  species_df$RFi <- (species_df$Fi / sum(species_df$Fi)) * 100
  
  # Chỉ số IVI
  species_df$IVI <- (species_df$Ai + species_df$Di + species_df$RFi) / 3
  
  # Chỉ số Shannon-Wiener (H') và Simpson (Cd)
  pi_vector <- species_df$n / total_trees
  H_prime <- -sum(pi_vector * log(pi_vector))
  Cd <- sum(pi_vector^2)
  
  return(list(Composition = species_df, H = H_prime, Cd = Cd))
}

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

1. Hiện trạng Phân bố Diện tích các Trạng thái Rừng

Tổng diện tích đất lâm nghiệp xã Thần Sa là 9.537,30 ha (chiếm 94% diện tích tự nhiên toàn xã), bao gồm: Đất rừng đặc dụng 5.900,59 ha; Đất rừng phòng hộ 2.588,41 ha; Đất rừng sản xuất 1.048,30 ha.

STT Trạng thái rừng Đặc trưng thảm phủ Diện tích (ha) Tỷ lệ (%)
1 IA Đất trống trảng cỏ nhiệt đới 54,16 0,57
2 IB Trảng cây bụi, cây gỗ rải rác 105,80 1,11
3 IC Cây gỗ tái sinh rải rác xen cây bụi 123,50 1,29
4 IIA Rừng phục hồi sau nương rẫy, gỗ nhỏ 3.272,50 34,31
5 IIB Rừng thứ sinh phục hồi nghèo 3.369,20 35,33
6 IIIA1 Rừng thứ sinh cây gỗ trung bình/nghèo kiệt 2.612,14 27,39
Tổng 9.537,30 100,00 
                       CƠ CẤU TRẠNG THÁI RỪNG XÃ THẦN SA
+-----------------------------------------------------------------------------------+
| Trạng thái IIB   [=====================================] 35.33% (3.369,2 ha)     |
| Trạng thái IIA   [===================================]   34.31% (3.272,5 ha)     |
| Trạng thái IIIA1 [============================]          27.39% (2.612,1 ha)     |
| Trạng thái IC    [=]                                      1.29% (  123,5 ha)     |
| Trạng thái IB    [=]                                      1.11% (  105,8 ha)     |
| Trạng thái IA    [>]                                      0.57% (   54,2 ha)     |
+-----------------------------------------------------------------------------------+

2. Đặc điểm Cấu trúc Trắc lượng Tầng Cây Gỗ

Kết quả đo đếm chi tiết tại 6 ô tiêu chuẩn đại diện (3 OTC trạng thái IIB và 3 OTC trạng thái IIIA1) phản ánh sự suy giảm trữ lượng rõ rệt:

Trạng thái Ô tiêu chuẩn (OTC) Vị trí thôn Mật độ ($N$, 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}$)
IIB OTC 1 Xuyên Sơn 495 14,21 11,5 5,16 26,65
IIB OTC 2 Xuyên Sơn 450 15,35 12,2 6,85 32,84
IIB OTC 1 Ngọc Sơn 2 760 13,80 10,8 7,32 36,45
IIIA1 OTC 2 Ngọc Sơn 1 530 16,42 13,8 7,26 39,06
IIIA1 OTC 1 Ngọc Sơn 2 475 17,10 14,5 7,85 45,20
IIIA1 OTC 2 Ngọc Sơn 2 815 16,85 14,1 8,32 55,12

Nhận xét cấu trúc: Tiết diện ngang tầng cây gỗ trạng thái IIB dao động từ $5,16 - 7,32\text{ m}^2/\text{ha}$ với trữ lượng rất thấp ($26,65 - 36,45\text{ m}^3/\text{ha}$). Trạng thái IIIA1 có trữ lượng dao động từ $39,06 - 55,12\text{ m}^3/\text{ha}$. Các giá trị này chỉ bằng 25-35% so với rừng nguyên sinh núi đá vôi chuẩn ($180 - 220\text{ m}^3/\text{ha}$), chứng minh tình trạng suy thoái cấu trúc lâm phần nghiêm trọng.

3. Phân tích Tổ thành Loài Ưu thế và Chỉ số IVI

  • Trạng thái IIB (OTC 1 - Xuyên Sơn): Có 14 loài cây gỗ, trong đó 10 loài có $IVI > 5%$. Công thức tổ thành sinh thái: $$\text{Tổ thành} = 12,86\text{Kh} + 10,49\text{Bu} + 8,98\text{Tra} + 7,81\text{Da} + 7,81\text{Th} + 7,77\text{Tr} + 7,44\text{Ken} + 7,44\text{Su} + 6,40\text{Ke} + 5,59\text{Sa}$$ (Ghi chú: Kh: Kháo, Bu: Bứa, Tra: Trám, Da: Dâu da xoan, Th: Thẩu tấu, Tr: Trẩu, Ken: Kẹn, Su: Sui, Ke: Kè đuôi dông, Sa: Sảng).
  • Trạng thái IIIA1 (OTC 2 - Ngọc Sơn 1): Có 16 loài cây gỗ, 11 loài đạt $IVI > 5%$. Loài Sấu (Dracontomelon duperreanum) chiếm ưu thế tuyệt đối ($IVI = 13,00%$), tiếp theo là Vàng anh (Saraca dives, $IVI = 10,88%$) và Dẻ (Castanopsis, $IVI = 8,44%$).
  • Trạng thái IIIA1 (OTC 2 - Ngọc Sơn 2): Mật độ đạt 815 cây/ha, tổ thành ưu thế gồm: Sấu ($9,60%$), Sau sau ($8,45%$), Sảng ($8,43%$), Bồ đề ($8,06%$), Chẹo ($7,96%$), Kháo ($7,62%$), Lim xẹt ($7,46%$).
  • Hiện trạng vắng bóng cây gỗ quý: Các loài đặc trưng đỉnh cao của núi đá vôi như Nghiến (Burretiodendron hsienmu), Lát hoa (Chukrasia tabularis), Trai lý (Garcinia fagraeoides) có giá trị $IVI$ tụt dốc dưới $4-7%$ hoặc hoàn toàn biến mất ở tầng tán ưu thế trên $20\text{ m}$, chỉ còn lại các cá thể non hoặc cụt ngọn.

4. Định lượng Đa dạng Sinh thái và Hệ số Tương đồng

Trạng thái Vị trí OTC Số lượng loài ($S$) Chỉ số Shannon-Wiener ($H'$) Chỉ số Simpson ($Cd$)
IIB OTC 1 (Xuyên Sơn) 14 2,578 0,079
IIB OTC 2 (Xuyên Sơn) 16 2,649 0,011
IIB OTC 1 (Ngọc Sơn 2) 18 2,806 0,063
IIIA1 OTC 2 (Ngọc Sơn 1) 16 2,696 0,071
IIIA1 OTC 1 (Ngọc Sơn 2) 16 2,642 0,077
IIIA1 OTC 2 (Ngọc Sơn 2) 15 2,602 0,082

Hệ số tương đồng Sorensen ($SI$) giữa các ô tiêu chuẩn dao động từ 0,12 đến 0,80 (đặc biệt giữa OTC 1 IIB và OTC 1 IIIA1 đạt $SI = 0,80$), minh chứng các trạng thái rừng đều có nguồn gốc phân hóa từ cùng một thảm thực vật nguyên sinh núi đá vôi nhưng bị thoái hóa ở các mức độ can thiệp khác nhau.

5. Kết quả Đánh giá Áp lực Nhân sinh và Nguy cơ Suy thoái Rừng (PRA)

Khảo sát 21 phiếu điều tra bán cấu trúc tại các cụm dân cư bản địa cho thấy bức tranh áp lực toàn diện:

                      MỨC ĐỘ ẢNH HƯỞNG CỦA CÁC TÁC ĐỘNG NHÂN SINH
+-----------------------------------------------------------------------------------+
| 1. Khai thác gỗ thương phẩm/làm nhà [====================================] 100.0%  |
| 2. Khai thác củi đốt hàng ngày      [====================================] 100.0%  |
| 3. Khai thác lâm sản ngoài gỗ (LSNG)[==============================>     ]  85.7%  |
| 4. Săn bắt động vật hoang dã        [=========================>          ]  71.4%  |
| 5. Chăn thả gia súc tự do (Dê/Bò)   [====================>               ]  57.1%  |
| 6. Xâm lấn đất rừng làm nương rẫy   [==========>                         ]  28.6%  |
| 7. Cháy rừng do bất cẩn mang lửa    [===>                                ]   9.5%  |
+-----------------------------------------------------------------------------------+
  • Mục đích khai thác gỗ: 100% (21/21 hộ) đốn gỗ dựng nhà sàn truyền thống và lấy củi đun nấu sinh hoạt hàng ngày; 9,52% khai thác bán lấy tiền mặt.
  • Đối tượng loài bị khai thác kiệt quệ nhất: Nghiến (100% số phiếu ghi nhận bị đốn hạ), Lát hoa (80,95%), Kháo (76,19%), Trai lý (66,67%), Xoan nhừ (52,38%), Trám (42,86%).
  • Nguy cơ từ khai khoáng: Tình trạng khai thác vàng sa khoáng trái phép tại các lưu vực sông suối gây bồi lắng dòng chảy, ô nhiễm nguồn nước ngầm và sạt lở chân núi đá vôi.

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

Đổi mới về Phương pháp luận (Technical Innovations)

  1. Chuẩn hóa kỹ thuật Ô tiêu chuẩn Vệ tinh thích ứng địa hình hiểm trở: Thay vì áp dụng cứng nhắc ô tiêu chuẩn vuông $2.000\text{ m}^2$ dễ gặp đứt gãy do vách đá karst dựng đứng, đề tài cải tiến mô hình 4 ô vệ tinh ($10\text{ m} \times 50\text{ m}$) của Zoehrer (1980), cho phép dịch chuyển góc linh hoạt (Đông Nam, Tây Nam) khi gặp trở ngại tự nhiên mà vẫn đảm bảo độ tin cậy thống kê đại diện cho lâm phần.
  2. Tích hợp Ma trận Sinh thái - Xã hội học (Eco-Sociological Matrix): Kết nối trực tiếp giữa các thông số suy giảm cấu trúc trắc lượng ($IVI$, $M$, $G$) với mức độ khai thác của từng nhóm cộng đồng dân cư thông qua điều tra PRA định lượng, tạo cơ sở thực chứng vững chắc cho việc thiết kế chính sách REDD+.

So sánh Hiệu quả với các Công trình Tiền nhiệm

Tiêu chí khoa học Đánh giá Lâm sinh Cổ điển (Thái Văn Trừng, 1978) Nghiên cứu Kiểm kê Rừng Cấp tỉnh (Sở NN&PTNT, 2008) Nghiên cứu Tích hợp tại Xã Thần Sa (Đề tài)
Tiêu chí phân loại Hình thái ngoại mạo sinh thái Bản đồ hiện trạng 3 loại rừng Định lượng $IVI$ kết hợp $H', Cd, SI$ chi tiết
Hiệu chỉnh độ dốc Ước tính trực giác Không hiệu chỉnh cục bộ 100% cạnh ô được giải tích qua $L/\cos(\alpha)$
Độ chính xác đo diện tích Sai số biên $15 - 20%$ Sai số vĩ mô $\pm 10%$ Giảm sai số trắc lượng xuống $< 3,5%$ nhờ GPS/GIS
Định vị nguyên nhân Không phân tích nhân khẩu Báo cáo hành chính chung Định danh chính xác 6 loài gỗ quý bị tận diệt

Đóng góp cho Ngành và Địa phương

  • Cung cấp bộ cơ sở dữ liệu định lượng gốc (Baseline Data) gồm 6 ô tiêu chuẩn mẫu với tọa độ GPS chuẩn xác, phục vụ công tác theo dõi diễn biến rừng định kỳ của Ban Quản lý KBTTN Thần Sa – Phượng Hoàng.
  • Lượng hóa chính xác mức độ suy thoái trữ lượng ($26,65 - 55,12\text{ m}^3/\text{ha}$ so với chuẩn $180\text{ m}^3/\text{ha}$), chứng minh lâm phần đang ở giai đoạn suy thoái thứ sinh nghiêm trọng, đòi hỏi phải chuyển trạng thái từ bảo tồn thụ động sang khoanh nuôi xúc tiến tái sinh tích cực có làm giàu rừng.

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

Kịch bản Ứng dụng Thực tế (Real-World Use Cases)

                    CƠ CHẾ PHỐI HỢP QUẢN LÝ VÀ PHÁT TRIỂN BỀN VỮNG
+------------------------------------------------------------------------------------+
|  1. PHÂN VÙNG QUẢN LÝ BẢO VỆ NGHIÊM NGẶT (5.900,59 ha Rừng Đặc Dụng)               |
|  - Cắm mốc ranh giới số hóa GPS; tuần tra kiểm soát chốt chặn cửa rừng             |
|  - Cấm tuyệt đối cưa máy cầm tay và khai thác gỗ Nghiến, Lát hoa, Trai lý          |
+------------------------------------------------------------------------------------+
                                         |
                                         v
+------------------------------------------------------------------------------------+
|  2. CHƯƠNG TRÌNH ĐỒNG QUẢN LÝ RỪNG & CHI TRẢ DVMTR (REDD+ / PES)                  |
|  - Ký hợp đồng khoán bảo vệ rừng trực tiếp cho 21 cụm hộ dân bản địa               |
|  - Chi trả tiền dịch vụ môi trường rừng gắn với tỷ lệ giữ vững trữ lượng G và M    |
+------------------------------------------------------------------------------------+
                                         |
                                         v
+------------------------------------------------------------------------------------+
|  3. CHUYỂN ĐỔI KỸ THUẬT & CẢI THIỆN SINH KẾ NÔNG LÂM KẾT HỢP                      |
|  - Hỗ trợ xây dựng hầm Biogas / Bếp đun cải tiến tiết kiệm 60-70% củi đốt          |
|  - Trồng làm giàu rừng bằng cây bản địa đa mục đích (Dược liệu dưới tán, Trám, Dẻ) |
+------------------------------------------------------------------------------------+

Chiến lược Triển khai và Lộ trình Thực hiện

  1. Giai đoạn 1 (Năm 1 - 2: Khẩn cấp & Ngăn chặn):
    • Ban hành quy ước, hương ước thôn bản về bảo vệ rừng đặc dụng tại các xóm Xuyên Sơn, Ngọc Sơn 1, Ngọc Sơn 2.
    • Kiểm kê và niêm phong các cơ sở cưa xẻ gỗ trái phép; giải tỏa dứt điểm các điểm khai thác vàng sa khoáng gây hại hệ sinh thái suối.
  2. Giai đoạn 2 (Năm 3 - 4: Phục hồi Cấu trúc & Tái sinh):
    • Triển khai khoanh nuôi xúc tiến tái sinh tự nhiên kết hợp trồng bổ sung 200 ha cây gỗ bản địa ưu thế sinh thái ($IVI > 8%$) như Kháo, Trám, Vàng anh, Sau sau tại trạng thái IIB.
    • Thay thế toàn bộ đàn dê chăn thả tự do trong rừng đặc dụng sang mô hình nuôi nhốt bán thâm canh tại vùng đệm.
  3. Giai đoạn 3 (Năm 5 trở đi: Bền vững & Tiếp cận REDD+):
    • Vận hành cơ chế đo đạc, báo cáo và thẩm định (MRV) trữ lượng carbon rừng dựa trên mạng lưới OTC định vị; đấu nối nguồn tài chính carbon quốc tế.

Phân tích Hiệu quả Kinh tế - Môi trường (Cost-Benefit Analysis)

  • Đầu tư công nghệ giảm tiêu thụ củi: Xây dựng 100 hầm biogas quy mô hộ gia đình và cấp phát 300 bếp đun cải tiến giúp giảm thiểu $1.200\text{ tấn củi khô/năm}$, bảo tồn tương đương $35 - 40\text{ ha}$ rừng phục hồi khỏi nguy cơ chặt cụt tán.
  • Tiềm năng hấp thụ Carbon: Khôi phục trữ lượng gỗ từ $35\text{ m}^3/\text{ha}$ lên $90\text{ m}^3/\text{ha}$ trong chu kỳ 15 năm tại $3.369,2\text{ ha}$ trạng thái IIB tương đương hấp thụ thêm khoảng $180.000\text{ tấn } \text{CO}_2$, mở ra nguồn thu hàng tỷ đồng từ tín chỉ carbon rừng.

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

Hạn chế Kỹ thuật và Nguồn lực

  • Quy mô mẫu điều tra dã ngoại: Do địa hình núi đá chia cắt hiểm trở và thời gian thực hiện đề tài giới hạn (4 tháng: 02/2016 – 05/2016), số lượng ô tiêu chuẩn mới dừng lại ở 6 OTC ($12.000\text{ m}^2$), tập trung chủ yếu ở hai trạng thái IIB và IIIA1.
  • Mẫu khảo sát xã hội học PRA: Phạm vi điều tra 21 hộ gia đình tuy mang tính đại diện điển hình nhưng chưa bao phủ toàn diện 100% các hộ dân thuộc 9 xóm bản trên toàn xã Thần Sa.
  • Thiếu hụt dữ liệu biến động chu kỳ dài: Nghiên cứu mang tính chất mặt cắt thời gian (Cross-sectional study), chưa có dữ liệu đo lặp hàng năm (Longitudinal data) trên ô tiêu chuẩn định vị vĩnh viễn (Permanent Sampling Plots).

Hướng phát triển Đề tài

  • Ứng dụng Công nghệ Viễn thám LiDAR và Drone (UAV): Sử dụng máy bay không người lái tích hợp cảm biến LiDAR quét mật độ điểm cao để tái lập mô hình chiều cao tán rừng (Canopy Height Model - CHM) trên toàn bộ vách núi đá vôi dốc đứng mà con người không thể tiếp cận.
  • Thiết lập Mạng lưới Ô tiêu chuẩn Cố định Vĩnh viễn (PSP): Đóng mốc bê tông kiên cố cho 20 ô tiêu chuẩn định vị GPS trên các đai cao khác nhau để theo dõi tốc độ tăng trưởng đường kính ($D_{1.3}$), sinh khối và trữ lượng carbon phục vụ dự án REDD+ cấp tỉnh.
  • Nghiên cứu Kỹ thuật Tái sinh Hạt giống Cây Gỗ Quý: Thử nghiệm nhân giống sinh dưỡng và kỹ thuật gieo ươm các loài cây có nguy cơ tuyệt chủng cục bộ như Nghiến, Trai lý trên giá thể đất karst nhân tạo.

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

+------------------------------------------------------------------------------------+
|                         MA TRẬN GIÁ TRỊ VÀ ĐỐI TƯỢNG HƯỞNG LỢI                     |
+------------------------------------------------------------------------------------+
|  1. Sinh viên & Giảng viên ngành Lâm nghiệp / Sinh thái                            |
|     -> Bộ dữ liệu mẫu thực tế về phương pháp OTC Zoehrer và giải thuật IVI, H', Cd |
|     -> Tài liệu tham khảo chuẩn mực cho đồ án trắc lượng rừng trên núi đá vôi      |
+------------------------------------------------------------------------------------+
|  2. Ban Quản lý KBTTN Thần Sa - Phượng Hoàng & Kiểm lâm                            |
|     -> Bản đồ và tọa độ GPS chính xác của các vùng rừng suy thoái trọng điểm       |
|     -> Căn cứ pháp lý và khoa học để xác lập danh mục tuần tra bảo vệ ưu tiên      |
+------------------------------------------------------------------------------------+
|  3. Chính quyền Xã Thần Sa & Cơ quan Quản lý Quy hoạch                             |
|     -> Cơ sở dữ liệu cơ cấu 6 trạng thái đất lâm nghiệp (9.537,30 ha)              |
|     -> Luận cứ phân bổ nguồn vốn các chương trình mục tiêu quốc gia và sinh kế     |
+------------------------------------------------------------------------------------+
|  4. Cộng đồng Dân cư Bản địa (Tày, Dao, H'Mông, Nùng)                              |
|     -> Cơ hội tiếp cận các dự án hỗ trợ hầm Biogas, bếp cải tiến, cây giống        |
|     -> Hưởng lợi lâu dài từ nguồn tiền chi trả dịch vụ môi trường rừng (PES)       |
+------------------------------------------------------------------------------------+
  • Sinh viên & Giới học thuật: Tiếp cận phương pháp luận định lượng sinh thái rừng chuẩn tắc; làm chủ quy trình kết hợp đo đạc dã ngoại, hiệu chỉnh góc dốc địa hình và lập trình thống kê sinh thái học.
  • Cơ quan Quản lý Rừng Đặc dụng: Nắm bắt chính xác tọa độ và danh sách loài bị suy thoái mạnh (Nghiến $100%$, Lát hoa $80,95%$) để bố trí các trạm kiểm soát lưu động tại cửa rừng.
  • Cộng đồng Địa phương: Chuyển đổi mô hình sinh kế từ khai thác rừng tự nhiên sang đồng quản lý và nhận khoán bảo vệ rừng, cải thiện thu nhập bình quân từ dịch vụ môi trường sinh thái.

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

1. Yêu cầu kỹ thuật cốt lõi để thiết lập ô tiêu chuẩn trên núi đá vôi hiểm trở là gì?

Cần trang bị máy đo độ dốc Suunto Clinometer để đo góc nghiêng sườn dốc $\alpha$ và áp dụng công thức $L' = L / \cos(\alpha)$ nhằm hiệu chỉnh cự ly thực tế trên thực địa, đảm bảo diện tích chiếu ngang của ô đạt chuẩn $2.000\text{ m}^2$ hoặc các ô vệ tinh $500\text{ m}^2$ ($10\text{ m} \times 50\text{ m}$). Đồng thời, máy định vị GPS phải được cài đặt hệ quy chiếu tọa độ VN-2000 kinh tuyến trục $106^\circ 30'$ tỉnh Thái Nguyên.

2. Làm thế nào để giải quyết tình trạng cưa máy cầm tay khai thác gỗ Nghiến lén lút trong KBT?

Giải pháp tích hợp gồm 3 tầng: (1) Cắm mốc ranh giới số hóa và lắp đặt các trạm camera bẫy ảnh/cảm biến âm thanh nhận diện tiếng cưa máy tại các hẻm núi hiểm trở; (2) Ký cam kết kiểm soát chuỗi cung ứng với các thợ mộc địa phương và cấm lưu thông gỗ thớt nghiến; (3) Hỗ trợ xi măng, ngói fibro và vật liệu thay thế để người dân bỏ tập quán làm nhà sàn bằng gỗ rừng tự nhiên.

3. Tại sao chỉ số Đa dạng sinh học Shannon-Wiener ($H'$) ở trạng thái IIB lại cao hơn IIIA1?

Ở trạng thái rừng phục hồi IIB (ví dụ OTC 1 Ngọc Sơn 2 đạt $H' = 2,806$), tầng tán chính bị phá vỡ tạo điều kiện cho nhiều loài cây tiên phong ưa sáng (Kháo, Bứa, Dâu da xoan, Sảng, Trẩu) cùng đồng thời xâm nhập và cạnh tranh với số lượng cá thể phân bố tương đối đồng đều. Trong khi đó, ở trạng thái IIIA1, một số loài cây tầng ưu thế còn sót lại như Sấu ($IVI = 13,00%$) hoặc Bồ đề chiếm tỷ trọng áp đảo, làm chỉ số ưu thế Simpson $Cd$ tăng lên và kéo giảm độ đồng đều sinh thái $H'$.

4. Chi phí đầu tư cho hệ thống bếp cải tiến và hầm biogas tại xã Thần Sa được thu hồi như thế nào?

Mỗi hầm biogas composite quy mô hộ gia đình có chi phí khoảng 10-12 triệu VNĐ, bếp đun cải tiến khoảng 300.000 - 500.000 VNĐ. Lợi ích thu hồi được tính bằng: (1) Tiết kiệm 2-3 giờ lao động lấy củi mỗi ngày/hộ (tương đương giá trị ngày công nông nhàn); (2) Giảm chi phí phân bón hóa học nhờ tận dụng phụ phẩm bùn thải biogas cho cây trồng; (3) Giảm thiểu phát thải khí nhà kính và bảo vệ nguyên vẹn tán rừng phục hồi.

5. Dữ liệu từ khóa luận này có thể tích hợp vào Hệ thống Đo đạc - Báo cáo - Thẩm định (MRV) của REDD+ không?

Hoàn toàn có thể. Toàn bộ thông số trắc lượng cây gỗ ($D_{1.3}, H_{vn}, G, M$) kết hợp với danh mục loài và tọa độ GPS của 6 ô tiêu chuẩn cung cấp bộ dữ liệu sinh khối trên mặt đất (Aboveground Biomass - AGB) gốc theo tiêu chuẩn IPCC Tier 2, phục vụ tính toán đường phát thải tham chiếu rừng (FREL/FRL) cho tỉnh Thái Nguyên.

Kết luận

Khóa luận tốt nghiệp "Đánh giá các nguy cơ gây suy thoái rừng tại xã Thần Sa thuộc 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" của sinh viên Nông Văn Tuần đã giải quyết trọn vẹn các mục tiêu khoa học và thực tiễn đề ra. Bằng việc kết hợp chặt chẽ giữa kỹ thuật trắc lượng rừng trên ô tiêu chuẩn hiệu chỉnh góc dốc địa hình karst, bản đồ không gian MapInfo và điều tra xã hội học nông thôn PRA, nghiên cứu đã:

  1. Lượng hóa chi tiết hiện trạng tài nguyên: Xác định diện tích phân bố 6 trạng thái rừng trên $9.537,30\text{ ha}$ đất lâm nghiệp, trong đó rừng nghèo và phục hồi chiếm ưu thế tuyệt đối (IIA: $34,31%$, IIB: $35,33%$, IIIA1: $27,39%$).
  2. Minh chứng khoa học về sự suy thoái cấu trúc: Trữ lượng gỗ tại các trạng thái IIB ($26,65 - 36,45\text{ m}^3/\text{ha}$) và IIIA1 ($39,06 - 55,12\text{ m}^3/\text{ha}$) cùng tiết diện ngang ($5,16 - 8,32\text{ m}^2/\text{ha}$) đang ở mức rất thấp; các loài cây gỗ đặc hữu có giá trị cao (Nghiến, Lát hoa, Trai lý) bị khai thác cạn kiệt khỏi tầng ưu thế sinh thái.
  3. Định danh chính xác chuỗi nguyên nhân nhân sinh: Khai thác gỗ làm nhà (100%), lấy củi đun nấu (100%), khai thác lâm sản ngoài gỗ (85,71%) và khai thác khoáng sản sa khoáng là những tác nhân trực tiếp phá hủy tính ổn định của hệ sinh thái núi đá vôi.
  4. Đề xuất hệ thống giải pháp toàn diện: Kết hợp bảo vệ nghiêm ngặt vùng lõi $5.900,59\text{ ha}$ với khoanh nuôi xúc tiến tái sinh, chuyển giao kỹ thuật năng lượng thay thế (biogas, bếp cải tiến) và thể chế hóa mô hình đồng quản lý rừng có sự tham gia của cộng đồng bản địa.

Đây là công trình nghiên cứu ứng dụng có giá trị học thuật cao, cung cấp cơ sở dữ liệu thực chứng cho các nhà quản lý rừng đặc dụng, nhà nghiên cứu sinh thái học và chính quyền địa phương trong việc hoạch định chiến lược bảo tồn đa dạng sinh học và phát triển sinh kế bền vững tại vùng Đông Bắc Việt Nam.