Giới thiệu dự án

Hệ sinh thái rừng tự nhiên (HST RTN) đóng vai trò sống còn trong việc bảo vệ môi trường, duy trì vòng tuần hoàn sinh địa hóa và bảo tồn đa dạng sinh học. Theo Báo cáo Đánh giá Hệ sinh thái Thiên niên kỷ (Millennium Ecosystem Assessment - MEA, 2005), suy thoái rừng tại các khu vực nhiệt đới làm suy giảm nghiêm trọng các dịch vụ điều tiết nguồn nước, phòng chống xói mòn và lưu trữ carbon. Tại Việt Nam, áp lực chuyển đổi mục đích sử dụng đất sang sản xuất nông nghiệp và trồng cây công nghiệp đã khiến nhiều diện tích rừng tự nhiên bị chia cắt và chuyển dịch thành trạng thái rừng nghèo kiệt.

Tại xã Rã Bản, huyện Chợ Đồn, tỉnh Bắc Kạn – khu vực miền núi có địa hình dốc từ 20° đến 25° và độ cao trên 600m – tài nguyên rừng tự nhiên đối mặt với áp lực khai thác lâm sản ngoài gỗ (LSNG) thiếu kiểm soát và hiện tượng thiên tai cực đoan (lũ quét năm 2006 gây vùi lấp và mất trắng 27,75 ha đất nông nghiệp; sạt lở đất nghiêm trọng năm 2014). Việc thiếu hụt dữ liệu định lượng về các dịch vụ hệ sinh thái gây khó khăn cho việc hoạch định chính sách Chi trả Dịch vụ Môi trường Rừng (PFES) và quản lý tài nguyên bền vững tại địa phương.

+-----------------------------------------------------------------------------------+
|                            MỤC TIÊU VÀ PHẠM VI NGHIÊN CỨU                         |
+-----------------------------------------------------------------------------------+
| 1. Định lượng trữ lượng gỗ, mật độ lâm phần và trạng thái sinh khối (AGB).        |
| 2. Đo lường khả năng lưu trữ carbon hữu cơ (C-pool) theo tiêu chuẩn IPCC.         |
| 3. Xác định đa dạng loài thực vật thân gỗ qua chỉ số Shannon-Wiener (H').         |
| 4. Thống kê, phân loại các nhóm dịch vụ cung ứng (LSNG, dược liệu, thực phẩm).     |
| 5. Đề xuất khung giải pháp quản trị rừng dựa vào cộng đồng (CBFM).                |
+-----------------------------------------------------------------------------------+

Dự án áp dụng phương pháp tiếp cận liên ngành kết hợp điều tra sinh thái thực địa theo ô tiêu chuẩn (OTC) và đánh giá nhanh nông thôn có sự tham gia của người dân (Participatory Rapid Appraisal - PRA). Kết quả nghiên cứu xác lập bộ chỉ số cơ sở về carbon, đa dạng sinh học và giá trị tài nguyên, hỗ trợ chính quyền địa phương triển khai chiến lược bảo tồn đa dạng sinh học gắn liền với sinh kế cộng đồng.

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

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

Phương pháp quản lý rừng truyền thống tại cấp xã chủ yếu dựa trên báo cáo hành chính và thống kê diện tích che phủ thô, bỏ qua việc định lượng các bể chứa carbon và dịch vụ sinh thái.

Tiêu chí phân tích Điều tra kiểm kê truyền thống Mô hình viễn thám GIS (InVEST) Phương pháp điều tra đa dịch vụ (Đề tài)
Độ chính xác mặt đất Trung bình (chỉ đo $D_{1.3}$ cây gỗ) Thấp đến Trung bình (thiếu ground-truth) Rất cao (OTC $2500m^2$ + ô dạng bản $1m^2$)
Bể carbon đo đếm Duy nhất tầng cây gỗ đứng Ước tính gián tiếp qua thảm phủ Đo toàn diện (Cây gỗ + Thảm tươi/Cây bụi + Thảm mục)
Đánh giá xã hội Không tích hợp Không tích hợp Tích hợp sâu (Phỏng vấn PRA, bảng hỏi bán cấu trúc)
Chi phí & Thiết bị Thấp Rất cao (bản quyền phần mềm) Tối ưu, khả thi chuyển giao kỹ thuật 
+-----------------------------------------------------------------------------------+
|                  KHUNG YÊU CẦU ĐỊNH LƯỢNG HỆ SINH THÁI (MoSCoW)                   |
+-----------------------------------------------------------------------------------+
| [Must Have]    : Đo đường kính D1.3 (cm), chiều cao H (m), sinh khối tươi/khô DW. |
| [Should Have]  : Chỉ số Shannon-Wiener (H'), trữ lượng carbon theo hướng dẫn IPCC.|
| [Could Have]   : Ma trận tương quan giữa suy giảm thảm phủ và tần suất thiên tai.  |
| [Won't Have]   : Lượng giá kinh tế toàn phần bằng phương pháp CVM phức tạp.       |
+-----------------------------------------------------------------------------------+

Thiết kế hệ thống

Quy trình thu thập và xử lý dữ liệu được thiết kế thành một chuỗi pipeline chuẩn hóa từ thực địa đến mô hình tính toán sinh khối và carbon tích lũy:

[Khảo sát thực địa] ---> [Thiết lập 6 OTC (50m x 50m)] ---> [Đo đếm D1.3 >= 5cm]
                                  |
                                  v
[Lập 5 Ô dạng bản (1m x 1m)] ---> [Thu thập Thảm tươi + Thảm mục] ---> [Sấy 80°C / 48h]
                                  |
                                  v
[Tính toán Sinh khối AGB]   ---> [Áp dụng Allometric Models & Hệ số Carbon IPCC 0.46]
                                  |
                                  v
[Phân tích Đa dạng sinh học] ---> [Chỉ số Shannon-Wiener (H') & Báo cáo Dịch vụ HST]

Technology Stack và các công cụ đo đếm:

  • Thiết bị đo đạc: Thước kính chuyên dụng, thước dây lâm nghiệp, GPS Garmin eTrex 20x, khung chuẩn $1\text{m} \times 1\text{m}$, cân điện tử độ chính xác $\pm 0.01\text{g}$, tủ sấy tĩnh nhiệt phân tích sinh khối.
  • Tiêu chuẩn lâm học: Thông tư 34/2009/TT-BNNPTNT về tiêu chí phân loại trạng thái rừng.
  • Công cụ tính toán & xử lý số liệu: Microsoft Excel 2016, Hệ thống phương trình Allometric Ketterings et al. (2001), Hướng dẫn kiểm kê khí nhà kính IPCC (2003).

Methodology

Dự án áp dụng mô hình điều tra lâm học kết hợp sinh thái học quần xã và xã hội học môi trường trong vòng 6 tháng:

+-----------------------------------------------------------------------------------+
|                          KẾ HOẠCH VÀ MILESTONES DỰ ÁN                             |
+-----------------------------------------------------------------------------------+
| Phase 1: Tổng quan tài liệu, thiết kế bảng hỏi và tuyến điều tra (Tháng 1-2).     |
| Phase 2: Khảo sát thực địa, lập 6 OTC tại thôn Pắc Giả và Khuổi Giả (Tháng 3-4). |
| Phase 3: Xử lý mẫu trong phòng thí nghiệm (sấy khô, xác định tỷ lệ DW) (Tháng 5). |
| Phase 4: Phân tích thống kê, tổng hợp trữ lượng carbon và viết báo cáo (Tháng 6).  |
+-----------------------------------------------------------------------------------+

Kiểm soát rủi ro và chất lượng dữ liệu:

  • Xử lý ranh giới OTC: Cây nằm trên ranh giới ô đo đếm được áp dụng quy tắc biên: chỉ tính các cây nằm trên cạnh trước và cạnh bên phải theo hướng tiến của tuyến khảo sát.
  • Kiểm soát độ ẩm: Mẫu thảm mục và thảm tươi ($200\text{g}$ mẫu tươi đại diện) được sấy chuẩn hóa ở $80^\circ\text{C}$ liên tục trong $48\text{ giờ}$ cho đến khi đạt khối lượng không đổi để triệt tiêu sai số bốc hơi nước.

Implementation và kết quả

Development process

  1. Phương pháp bố trí ô tiêu chuẩn: Thiết lập 6 OTC diện tích $2500\text{ m}^2$ ($50\text{m} \times 50\text{m}$) ngẫu nhiên có hệ thống đại diện cho các trạng thái rừng tự nhiên nghèo tại 2 thôn Pắc Giả (PG1, PG2, PG3) và Khuổi Giả (KG1, KG2, KG3).
  2. Mô hình toán học và thuật toán xử lý dữ liệu:
    • Mật độ cây rừng ($N/\text{ha}$): $$N/\text{ha} = \frac{n}{S} \times 10000$$ (Trong đó: $n$ là tổng số cây đo đếm có $D_{1.3} \ge 5\text{cm}$, $S = 2500\text{ m}^2$ là diện tích OTC).
    • Sinh khối cây gỗ trên mặt đất (AGB - Above Ground Biomass): Sử dụng phương trình sinh thái nhiệt đới của Ketterings et al. (2001): $$Y = 0.11 \times \rho \times (D_{1.3})^{2.62}$$ (Trong đó: $Y$ là sinh khối khô $(\text{kg/cây})$, $D_{1.3}$ là đường kính ngang ngực $(\text{cm})$, $\rho$ là dung trọng gỗ trung bình $= 0.62\text{ g/cm}^3$).
    • Trữ lượng carbon tích lũy ($W_{\text{carbon}}$): Áp dụng hệ số carbon mặc định theo IPCC (2003): $$W_{\text{carbon}} = DW \times 0.46$$
    • Chỉ số đa dạng sinh học Shannon-Wiener ($H'$): $$H' = -\sum_{i=1}^{s} \left( \frac{n_i}{N} \right) \ln \left( \frac{n_i}{N} \right)$$

Đoạn mã Python triển khai tính toán sinh khối và chỉ số sinh thái:

import numpy as np
import pandas as pd

def calculate_forest_metrics(tree_data, understory_dw_per_ha, litter_dw_per_ha):
    """
    Tính toán sinh khối, carbon và chỉ số Shannon-Wiener cho ô tiêu chuẩn
    """
    # 1. Sinh khối cây gỗ (Ketterings et al. 2001), rho = 0.62
    tree_data['AGB_kg'] = 0.11 * 0.62 * (tree_data['dbh_cm'] ** 2.62)
    agb_tonnes_ha = (tree_data['AGB_kg'].sum() / 1000.0) * (10000.0 / 2500.0)
    
    # 2. Tổng sinh khối khô (tonnes/ha)
    total_biomass = agb_tonnes_ha + understory_dw_per_ha + litter_dw_per_ha
    
    # 3. Trữ lượng Carbon (IPCC Tier 2 Factor = 0.46)
    carbon_tree = agb_tonnes_ha * 0.46
    carbon_understory = understory_dw_per_ha * 0.46
    carbon_litter = litter_dw_per_ha * 0.46
    total_carbon = total_biomass * 0.46
    
    # 4. Chỉ số đa dạng sinh học Shannon-Wiener (H')
    species_counts = tree_data['species'].value_counts().values
    total_individuals = len(tree_data)
    proportions = species_counts / total_individuals
    shannon_h = -np.sum(proportions * np.log(proportions))
    
    return {
        "AGB_Tree_t_ha": round(agb_tonnes_ha, 2),
        "Total_Biomass_t_ha": round(total_biomass, 2),
        "Total_Carbon_t_ha": round(total_carbon, 2),
        "Tree_Density_N_ha": int(total_individuals * (10000 / 2500)),
        "Shannon_Index_H": round(shannon_h, 5)
    }

Testing và validation

Số liệu thực nghiệm qua 6 ô tiêu chuẩn được kiểm chuẩn độc lập, phân tích đối chiếu giữa các bể chứa:

+-----------------------------------------------------------------------------------+
|                        TỔNG HỢP KẾT QUẢ ĐO ĐẾM THỰC ĐỊA                           |
+-----------------------------------------------------------------------------------+
| Chỉ tiêu đo đếm                 | Giá trị Min   | Giá trị Max   | Trung bình      |
+---------------------------------+---------------+---------------+-----------------+
| Trữ lượng gỗ lâm phần (m3/ha)   | 11.22         | 14.51         | 12.95 (Nghèo)   |
| Mật độ cây gỗ (cây/ha)          | 108           | 132           | 118             |
| Sinh khối tầng cây gỗ (tấn/ha)  | 10.27         | 14.34         | 12.31           |
| Sinh khối tầng thảm mục (tấn/ha)| 2.69          | 3.30          | 2.99            |
| Sinh khối cây bụi/tươi (tấn/ha) | 1.72          | 2.28          | 1.98            |
| Tổng trữ lượng Carbon (tấn C/ha)| 6.75          | 8.96          | 7.95            |
| Chỉ số Shannon-Wiener (H')      | 2.51815       | 3.04439       | 2.79340         |
+-----------------------------------------------------------------------------------+

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

  1. Dịch vụ cung ứng (Provisioning Services):
    • Cây dược liệu: Xác định 17 loài bản địa thuộc nhiều họ thực vật (Sa nhân, Bách bộ, Hoàng đằng già, Chân chim, Thảo đậu khấu, Sim, Thồm lồm, Ráy đứng, Bạc thau...).
    • Cây thực phẩm: Khai thác 16 loài bản địa có giá trị thương phẩm cao (Rau bò khai, Rau sắng, Rau dớn, Măng rừng, Quả trám đen, Trám trắng, Giảo cổ lam...).
    • Động vật rừng: Ghi nhận 6 loài thường gặp (Gà rừng, sóc, dúi, chồn, rắn, cầy), trong đó nhóm rắn và cầy nằm trong danh mục cấm săn bắt trái phép.
  2. Dịch vụ điều tiết và hỗ trợ (Regulating & Supporting Services):
    • Điều tiết nước: 7 công trình nước sạch tự chảy phục vụ 372 hộ dân (đạt tỷ lệ 98,1%); cung cấp nước tưới ổn định cho 145,19 ha đất trồng trọt (54 ha tại Pắc Giả, 12 ha tại Khuổi Giả).
    • Phòng hộ đất dốc: Hệ thống thảm mục dày ($2,69 - 3,30\text{ tấn/ha}$) đóng vai trò giảm thiểu dòng chảy mặt, bảo vệ lớp đất canh tác trong lưu vực.

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

  • Tích hợp đa tầng bể carbon: Thay vì chỉ ước lượng tầng cây gỗ như các nghiên cứu lâm nghiệp truyền thống, đề tài giải mã định lượng đồng thời cả 3 tầng sinh khối (cây gỗ, cây bụi/thảm tươi, thảm mục), chứng minh tầng thảm mục đóng góp tới $15% - 18%$ tổng trữ lượng carbon của rừng nghèo nhiệt đới.
  • Tiếp cận sinh thái học cộng đồng: Kết hợp công cụ PRA với giải pháp lâm sinh, cung cấp danh mục 33 loài thực vật hữu ích (dược liệu và thực phẩm) làm căn cứ xây dựng mô hình nông lâm kết hợp dưới tán rừng.
  • Bằng chứng thực tiễn về giá trị rừng phục hồi: Chứng minh dù trữ lượng gỗ thấp ($11,22 - 14,51\text{ m}^3/\text{ha}$), trạng thái rừng nghèo tại Bắc Kạn vẫn duy trì chỉ số đa dạng loài ở mức trung bình ($H' = 2,793$) và lưu giữ trung bình $7,95\text{ tấn C/ha}$.
+-----------------------------------------------------------------------------------+
|                  SO SÁNH ĐÓNG GÓP VỚI CÁC CÔNG TRÌNH CÙNG LĨNH VỰC                |
+-----------------------------------------------------------------------------------+
| 1. Bùi Thị Minh Nguyệt & Trần Quang Bảo (2013): Nghiên cứu VQG Ba Vì chỉ tập      |
|    trung vào giá trị du lịch sinh thái (97% doanh thu từ thuê môi trường rừng).   |
| 2. Đinh Đức Trường (2009): Lượng giá đất ngập nước Ba Lạt (196 tỷ đồng) bằng kinh |
|    tế lượng, thiếu dữ liệu thực địa chi tiết về cấu trúc loài lâm sinh.           |
| 3. Đề tài tại Rã Bản: Cung cấp đầy đủ tham số vi mô về mật độ, loài chỉ thị,     |
|    sinh khối từng tầng và bộ giải pháp phục hồi cộng đồng.                        |
+-----------------------------------------------------------------------------------+

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

[Xã Rã Bản / Huyện Chợ Đồn]

Kịch bản triển khai thực tế

  1. Phát triển kinh tế dưới tán rừng: Quy hoạch vùng thu hái bền vững cho các loài có giá trị kinh tế cao như Rau bò khai, Rau sắng, Sa nhân tím tại thôn Pắc Giả và Khuổi Giả, nâng cao thu nhập cho 379 hộ dân.
  2. Chi trả dịch vụ môi trường rừng (PFES): Áp dụng dữ liệu trữ lượng carbon ($6,75 - 8,96\text{ tấn C/ha}$) và chức năng bảo vệ nguồn nước ngọt (cấp nước 98,1% dân cư) làm căn cứ phân bổ nguồn thu PFES từ các nhà máy thủy điện và cơ sở cấp nước hạ lưu.
  3. Phòng chống thiên tai lưu vực: Tăng cường mật độ cây đứng từ 118 cây/ha lên mức tiêu chuẩn bền vững (>300 cây/ha) thông qua khoanh nuôi tái sinh tự nhiên có xúc tiến, hạn chế nguy cơ trượt lở đất và lũ quét mùa mưa.

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

  • Hạn chế kỹ thuật: Chưa thực hiện khoan lấy mẫu phân tích carbon trong đất (Soil Organic Carbon - SOC) và rễ cây ngầm (Below Ground Biomass - BGB) do điều kiện thiết bị phân tích dã ngoại có hạn.
  • Giới hạn không gian: Phạm vi thực nghiệm tập trung vào 6 ô tiêu chuẩn trên trạng thái rừng nghèo kiệt, chưa mở rộng sang các trạng thái rừng trung bình và rừng giàu trên núi đá vôi.
  • Hướng phát triển:
    1. Tích hợp dữ liệu viễn thám độ phân giải cao (Sentinel-2, LiDAR) để xây dựng bản đồ số hóa không gian phân bố carbon toàn huyện Chợ Đồn.
    2. Ứng dụng mô hình SWAT (Soil and Water Assessment Tool) để lượng hóa dòng chảy mặt và hệ số giữ nước của lưu vực sông suối tại xã Rã Bản.

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

  • Sinh viên & Học viên Lâm nghiệp/Môi trường: Tiếp cận phương pháp thiết lập ô tiêu chuẩn $2500\text{ m}^2$, kỹ thuật đo đếm lâm học và code mẫu tính toán chỉ số $H'$, Allometric equations.
  • Cán bộ Kiểm lâm & Ban Quản lý Rừng phòng hộ: Sở hữu bộ số liệu định lượng về trữ lượng gỗ ($12,95\text{ m}^3/\text{ha}$) và mật độ loài để lập phương án quản lý rừng bền vững cấp xã.
  • Chính quyền địa phương & Doanh nghiệp Carbon: Có dữ liệu nền tảng xác thực (MRV - Measurement, Reporting, Verification) phục vụ các dự án tín chỉ carbon rừng và chính sách chi trả PFES.
  • Cộng đồng dân cư địa phương: Hưởng lợi trực tiếp từ việc bảo vệ nguồn nước sinh hoạt ổn định và khai thác bền vững 33 loài dược liệu - thực phẩm bản địa.

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

1. Phương pháp tính sinh khối của Ketterings et al. (2001) có ưu điểm gì so với các mô hình khác?

Phương pháp Ketterings et al. (2001) được tối ưu hóa cho rừng nhiệt đới thứ sinh với hệ số tương quan phụ thuộc trực tiếp vào đường kính ngang ngực ($D_{1.3}$) và dung trọng gỗ ($\rho = 0.62\text{ g/cm}^3$). Mô hình này loại bỏ sự phụ thuộc vào chiều cao cây ($H$) vốn có sai số đo đạc rất lớn ở địa hình núi dốc rừng rậm nhiệt đới.

2. Tại sao trữ lượng carbon của rừng tự nhiên Rã Bản chỉ đạt 6.75 - 8.96 tấn C/ha?

Khu vực nghiên cứu có 100% diện tích khảo sát là rừng nghèo kiệt (trữ lượng gỗ đứng bình quân chỉ $12,95\text{ m}^3/\text{ha}$, mật độ cây lớn chỉ 118 cây/ha). Đây là giai đoạn đầu của quá trình phục hồi sau nương rẫy, sinh khối chủ yếu tập trung ở tầng cây non, thảm tươi và thảm mục.

3. Tầng thảm mục có vai trò như thế nào trong cấu trúc sinh khối tại khu vực?

Thảm mục tại Rã Bản đạt khối lượng khô từ $2,69$ đến $3,30\text{ tấn/ha}$, tương đương lượng tích lũy carbon từ $1,24$ đến $1,52\text{ tấn C/ha}$ (chiếm gần $18%$ tổng carbon mặt đất). Tầng này đóng vai trò quan trọng trong việc thấm giữ nước mưa, chống xói mòn bề mặt trên sườn dốc $20^\circ - 25^\circ$.

4. Hệ số carbon 0.46 của IPCC được áp dụng như thế nào?

Theo Hướng dẫn Đánh giá Khí nhà kính Quốc gia của IPCC (2003) dành cho rừng nhiệt đới, hàm lượng carbon nguyên tố trung bình chiếm $46%$ tổng sinh khối khô hữu cơ ($DW$). Trữ lượng carbon được tính bằng công thức: $W_{\text{carbon}} = DW \times 0.46$.

5. Giải pháp nào khả thi nhất để phục hồi rừng tự nhiên tại xã Rã Bản?

Giải pháp tối ưu là khoanh nuôi xúc tiến tái sinh tự nhiên kết hợp làm giàu rừng bằng cây bản địa đa mục đích (Trám, Sa nhân, Bò khai), kết hợp giao khoán bảo vệ rừng gắn với quy ước thôn bản (hương ước) của cộng đồng Pắc Giả và Khuổi Giả.

Kết luận

Đề tài "Nghiên cứu xác định các loại dịch vụ hệ sinh thái rừng tự nhiên tại xã Rã Bản, huyện Chợ Đồn, tỉnh Bắc Kạn" đã định lượng thành công cấu trúc sinh thái và giá trị dịch vụ của trạng thái rừng tự nhiên nghèo nhiệt đới:

  • Xác định mật độ cây rừng trung bình đạt 118 cây/ha, trữ lượng bình quân 12,95 m³/ha.
  • Đo lường tổng sinh khối khô mặt đất đạt 17,28 tấn/ha, tương ứng trữ lượng carbon trung bình 7,95 tấn C/ha.
  • Đánh giá chỉ số đa dạng loài Shannon-Wiener đạt mức trung bình khá ($H' = 2,793$), ghi nhận 33 loài LSNG làm thực phẩm và dược liệu giá trị cao.
  • Chứng minh vai trò điều tiết trọng yếu của hệ sinh thái trong việc cung cấp nước sạch cho 98,1% số hộ dân và bảo vệ an toàn lưu vực canh tác nông nghiệp trước nguy cơ lũ quét, sạt lở đất.

Bộ dữ liệu thực nghiệm là cơ sở khoa học tin cậy hỗ trợ các nhà quản lý lâm nghiệp, chính quyền địa phương và các nhà nghiên cứu xây dựng chính sách lâm phận bền vững, phát triển sinh kế dưới tán rừng và tích hợp cơ chế chi trả dịch vụ môi trường rừng hiệu quả tại vùng cao Đông Bắc Việt Nam.