Giới thiệu dự án

Theo báo cáo của Ban Liên chính phủ về Biến đổi Khí hậu (IPCC), hệ sinh thái rừng toàn cầu hiện lưu giữ hơn 547 tỷ tấn các-bon (C), trong đó rừng mưa nhiệt đới chiếm đến 62,16% tổng lượng các-bon lục địa. Sự gia tăng nồng độ $\text{CO}_2$ trong khí quyển – nguyên nhân chính khiến nhiệt độ Trái Đất dự báo tăng từ 1,5°C đến 4,5°C vào năm 2050 – đòi hỏi những giải pháp định lượng trữ lượng sinh khối rừng chính xác nhằm cung cấp cơ sở khoa học cho cơ chế Chi trả dịch vụ môi trường rừng (PFES) và chương trình Giảm phát thải từ mất rừng và suy thoái rừng (REDD+).

Tại Việt Nam, Khu Bảo tồn Thiên nhiên (BTTN) Thần Sa – Phượng Hoàng (huyện Võ Nhai, tỉnh Thái Nguyên) là khu vực rừng trên núi đá vôi đặc dụng có tính đa dạng sinh học cao với diện tích đất lâm nghiệp lên đến 9.037,58 ha (chiếm 88,06% tổng diện tích tự nhiên của xã Thần Sa). Tuy nhiên, các dữ liệu định lượng cụ thể về sinh khối và trữ lượng tích lũy các-bon theo từng trạng thái rừng phục hồi (IIb) và rừng nghèo (IIIa1) trên nền địa hình Karst hiểm trở tại đây vẫn còn thiếu hụt trầm trọng.

                                SƠ ĐỒ PHÂN BỐ BỂ CHỨA CÁC-BON RỪNG
+-----------------------------------------------------------------------------------+
|                        TỔNG TRỮ LƯỢNG CÁC-BON TRÊN MẶT ĐẤT                        |
|                               (Above-Ground Carbon - AGC)                         |
+----------------------------------------+------------------------------------------+
                                         |
         +-------------------------------+-------------------------------+
         |                               |                               |
+--------v-------+              +--------v-------+              +--------v-------+
|  TẦNG CÂY GỖ   |              | CÂY BỤI / TÁI  |              |   THẢM MỤC /   |
| (D1.3 >= 5 cm) |              | SINH / THẢM    |              | VẬT RƠI RỤNG   |
| Tỷ trọng:      |              | TƯƠI (<5 cm)   |              | Tỷ trọng:      |
| 50.0% - 61.8%  |              | Tỷ trọng:      |              | 20.6% - 31.8%  |
|                |              | 16.8% - 17.6%  |              |                |
+----------------+              +----------------+              +----------------+

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

  1. Xác định đặc điểm cấu trúc lâm phần và tổ thành loài: Đo đếm chỉ số mức độ quan trọng (Importance Value Index - IVI), mật độ thân cây (cây/ha), tiết diện ngang ($G - \text{m}^2/\text{ha}$), đường kính ngang ngực ($D_{1.3}$) và chiều cao vút ngọn ($H_{vn}$) của các trạng thái rừng IIb và IIIa1.
  2. Định lượng sinh khối tươi và sinh khối khô (Dry Weight - DW): Phân tích chi tiết 3 bể chứa các-bon trên mặt đất: Tầng cây gỗ sống ($D_{1.3} \ge 5\text{ cm}$), tầng dưới tán (cây bụi, thảm tươi, cây tái sinh $D_{1.3} < 5\text{ cm}$) và tầng thảm mục (litter layer).
  3. Ước tính tổng trữ lượng tích lũy $\text{C}$ và $\text{CO}_2$ tương đương: Quy đổi sinh khối khô sang các-bon bằng hệ số mặc định IPCC ($CF = 0,46$), tích hợp dữ liệu GIS để lượng hóa trên toàn cảnh quan 5.981,34 ha rừng tự nhiên xã Thần Sa.

Phạm vi & Giới hạn: Đề tài tập trung đo đếm 6 ô tiêu chuẩn (OTC) cố định dạng cụm vệ tinh đại diện trên các vị trí chân, sườn và đỉnh núi đá vôi tại các xóm Xuyên Sơn, Ngọc Sơn 1, Ngọc Sơn 2; không nghiên cứu bể chứa các-bon dưới đất (rễ cây) và các-bon hữu cơ trong đất (SOC).

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á Ưu điểm Nhược điểm Độ chính xác Khả thi thực địa
Chặt hạ toàn phần (Destructive Sampling) Đo đếm trực tiếp, chính xác 100% sinh khối thực tế Phá hủy cấu trúc rừng, bị nghiêm cấm trong Khu BTTN, tốn kém chi phí Rất cao ($\pm 1%$) Không khả thi tại vùng lõi bảo tồn
Viễn thám thuần túy (Optical Remote Sensing) Phủ diện tích lớn, theo dõi biến động nhanh qua chuỗi thời gian Bị bão hòa chỉ số NDVI ở tán rừng rậm, sai số cao ở địa hình núi đá Trung bình ($\pm 25% - 35%$) Khả thi nhưng độ tin cậy thấp
Ô tiêu chuẩn lồng nhau + Allometric Modeling (Giải pháp đề xuất) Bảo tồn nguyên vẹn hệ sinh thái, đo đếm đa tầng, kết hợp mô hình allometric chuẩn quốc tế Đòi hỏi nhân lực điều tra thực địa trên vách đá dốc, sấy mẫu phòng thí nghiệm Cao ($\pm 5% - 8%$) Hoàn toàn khả thi và tối ưu

Yêu cầu kỹ thuật theo mô hình MoSCoW:

  • Must have: Hệ thống ô tiêu chuẩn lồng nhau ($2.000\text{ m}^2$), phương trình tương quan sinh khối $Y = \exp(-2.134 + 2.530 \times \ln(D_{1.3}))$, tủ sấy mẫu chuẩn $85^\circ\text{C}$ trong 48 giờ, phần mềm GIS MapInfo Pro v15.0.
  • Should have: Phân tích chỉ số ưu thế sinh thái IVI $\ge 5%$ cho từng loài cây gỗ.
  • Could have: Tích hợp công cụ đánh giá nhanh trữ lượng các-bon RacSA thuộc bộ công cụ TULSEA (ICRAF).
  • Won't have (lần này): Quét LiDAR 3D tán rừng mặt đất do hạn chế kinh phí.

Thiết kế hệ thống điều tra và xử lý dữ liệu

Hệ thống điều tra sinh khối rừng được chuẩn hóa theo mô hình ô tiêu chuẩn vệ tinh 4 ô thứ cấp của Zoehrer (1980) có cải tiến:

Technology Stack & Tiêu chuẩn kỹ thuật

  • Nền tảng GIS & Bản đồ: MapInfo Professional v15.0, hệ quy chiếu VN-2000 kinh tuyến trục $105^\circ 45'$.
  • Mô hình toán sinh thái: Allometric Equations theo Brown et al. (1989) & Ketterings et al. (2001).
  • Hệ số phát thải / lưu giữ quốc tế: IPCC Good Practice Guidance for LULUCF (2003) với $C = \text{Biomass} \times 0,46$.
  • Thiết bị đo đếm thực địa: Thước đo đường kính Richter 5m, thước đo độ cao Blume-Leiss, cân điện tử chính xác $0,01\text{g}$, tủ sấy Memmert Universal Oven UN55.

Implementation và kết quả

Development Process & Thuật toán tính toán

Quy trình tính toán sinh khối và trữ lượng các-bon được tự động hóa thông qua pipeline phân tích toán sinh thái. Thuật toán xử lý số liệu ô tiêu chuẩn được thiết lập như sau:

import numpy as np
import pandas as pd

def calculate_forest_carbon(tree_data, understorey_fw, understorey_sub_dw, understorey_sub_fw,
                            litter_fw, litter_sub_dw, litter_sub_fw, plot_area_ha=0.2):
    """
    Tính toán sinh khối khô (DW) và trữ lượng Các-bon (tấn C/ha)
    tree_data: DataFrame chứa ['species', 'dbh_cm', 'h_m']
    """
    # 1. Sinh khối cây gỗ sống (kg/cây) theo Brown et al. (1989) cho rừng mưa nhiệt đới
    # ln(Y) = -2.134 + 2.530 * ln(D) <=> Y = exp(-2.134 + 2.530 * ln(D))
    tree_biomass_kg = np.exp(-2.134 + 2.530 * np.log(tree_data['dbh_cm']))
    total_tree_biomass_ton_ha = (tree_biomass_kg.sum() / 1000.0) / plot_area_ha

    # 2. Sinh khối tầng dưới tán (cây bụi, thảm tươi) (tấn/ha)
    # DW = (Sub_DW / Sub_FW) * Total_FW
    understorey_dw_g_m2 = (understorey_sub_dw / understorey_sub_fw) * understorey_fw
    understorey_biomass_ton_ha = understorey_dw_g_m2 * (10000.0 / (1000.0 * 1000.0))

    # 3. Sinh khối tầng thảm mục (tấn/ha)
    litter_dw_g_m2 = (litter_sub_dw / litter_sub_fw) * litter_fw
    litter_biomass_ton_ha = litter_dw_g_m2 * (10000.0 / (1000.0 * 1000.0))

    # 4. Quy đổi sang hàm lượng Carbon (IPCC default factor = 0.46)
    IPCC_FACTOR = 0.46
    c_tree = total_tree_biomass_ton_ha * IPCC_FACTOR
    c_understorey = understorey_biomass_ton_ha * IPCC_FACTOR
    c_litter = litter_biomass_ton_ha * IPCC_FACTOR
    c_total_ha = c_tree + c_understorey + c_litter

    return {
        "Tree_C_ton_ha": round(c_tree, 4),
        "Understorey_C_ton_ha": round(c_understorey, 4),
        "Litter_C_ton_ha": round(c_litter, 4),
        "Total_AGC_ton_ha": round(c_total_ha, 4),
        "Total_CO2e_ton_ha": round(c_total_ha * (44.0 / 12.0), 4)
    }

          CÔNG THỨC XÁC ĐỊNH CHỈ SỐ MỨC ĐỘ QUAN TRỌNG (IVI %)
          
             IVI(%) = [ A_i(%) + F_i(%) + D_i(%) ] / 3   (hoặc tổng tương đối)
  Trong đó:
  - A_i (Độ phong phú tương đối): A_i = (N_i / N_total) * 100
  - F_i (Tần suất xuất hiện tương đối): F_i = (n_i / n_total_plots) * 100
  - D_i (Độ ưu thế tương đối theo tiết diện ngang): D_i = (G_i / G_total) * 100

Testing và Validation thực địa

Dữ liệu điều tra từ 6 ô tiêu chuẩn đại diện tại Khu BTTN Thần Sa – Phượng Hoàng cho thấy các đặc trưng cấu trúc lâm phần:

Mã OTC Trạng thái rừng Vị trí địa hình Độ dốc (độ) Độ tàn che (%) Mật độ (cây/ha) $\bar{D}_{1.3}$ (cm) $\bar{H}_{vn}$ (m) Tiết diện ngang $G$ ($\text{m}^2/\text{ha}$) Trữ lượng gỗ ($m^3/\text{ha}$)
OTC 1 IIb Đỉnh núi $25^\circ$ 75% 495 10,46 11,15 5,16 26,65
OTC 2 IIb Đỉnh núi $22^\circ$ 70% 450 11,19 12,03 5,43 31,07
OTC 3 IIb Chân núi $45^\circ$ 75% 760 10,89 11,46 7,32 36,45
OTC 4 IIIa1 Chân núi $40^\circ$ 85% 530 11,88 12,87 7,26 39,06
OTC 5 IIIa1 Sườn núi $35^\circ$ 85% 475 12,38 13,75 6,27 42,85
OTC 6 IIIa1 Sườn núi $30^\circ$ 80% 815 11,47 14,45 8,32 55,12

Phân tích tổ thành loài cây gỗ:

  • Trạng thái IIb (Rừng phục hồi sau nương rẫy): Tổ thành loài chiếm ưu thế gồm Kháo (Actinodaphne sp., IVI = 12,86%), Vàng anh (Saraca dives, IVI = 11,74%), Bứa (Garcinia oblongifolia, IVI = 10,49%), Bồ đề (Styrax tonkinensis, IVI = 9,19%), Trám (Canarium album, IVI = 9,48%).
  • Trạng thái IIIa1 (Rừng nghèo): Xuất hiện nhiều loài gỗ lớn bản địa có giá trị bảo tồn như Sấu (Dracontomelon duperreanum, IVI = 12,99%), Vàng anh (IVI = 10,87%), Dẻ (Castanopsis sp., IVI = 8,43%), Lim xẹt (Peltophorum dasyrrhachis, IVI = 7,45%), Nghiến (Burretiodendron hsienmu), Lát hoa (Chukrasia tabularis).

Kết quả tích lũy sinh khối và các-bon đạt được

Tổng hợp dữ liệu phân tích phòng thí nghiệm xác định sinh khối khô và trữ lượng các-bon chi tiết theo từng thành phần:

              TRỮ LƯỢNG CÁC-BON THEO CÁC BỂ CHỨA (TẤN C/HA)
  35 +-----------------------------------------------------------------------+
     |                                                          [33.56]      |
  30 |                                                           ######      |
     |                                      [27.81]              ######      |
  25 |                                       ######              ######      |
     |                                       ######              ##++##      |
  20 |                                       ##++##              ##++##      |
     |                                       ##++##              ##**##      |
  15 |                                       ##**##              ##**##      |
     |                                       ##**##              ##**##      |
  10 |                                       ##**##              ##**##      |
     |                                       ##**##              ##**##      |
   5 |                                       ##**##              ##**##      |
     |                                       ##**##              ##**##      |
   0 +-----------------------------------------------------------------------+
                                           Trạng thái IIb      Trạng thái IIIa1
     [Chú giải: (**): Cây gỗ sống | (++) : Thảm mục | (##) : Cây bụi thảm tươi]
Hợp phần tích lũy Trạng thái IIb (TB $\pm$ SD) Tỷ lệ IIb (%) Trạng thái IIIa1 (TB $\pm$ SD) Tỷ lệ IIIa1 (%) Mức tăng trưởng IIIa1 vs IIb
Cây gỗ sống ($D_{1.3} \ge 5\text{cm}$) $14,30 \pm 0,81\text{ tấn C/ha}$ 51,40% $20,73 \pm 1,62\text{ tấn C/ha}$ 61,77% +44,97%
Cây bụi, thảm tươi, tái sinh $4,67 \pm 0,31\text{ tấn C/ha}$ 16,79% $5,90 \pm 0,75\text{ tấn C/ha}$ 17,58% +26,34%
Tầng thảm mục, rơi rụng $8,85 \pm 0,13\text{ tấn C/ha}$ 31,81% $6,93 \pm 0,98\text{ tấn C/ha}$ 20,65% -21,69%
Tổng tích lũy AGC/ha $27,81 \pm 0,95\text{ tấn C/ha}$ 100,0% $33,56 \pm 1,44\text{ tấn C/ha}$ 100,0% +20,68%
Quy đổi đương lượng $\text{CO}_2$ $101,97\text{ tấn }\text{CO}_2/\text{ha}$ - $123,06\text{ tấn }\text{CO}_2/\text{ha}$ - +20,68%

Tổng trữ lượng các-bon toàn cảnh quan xã Thần Sa:

  • Diện tích trạng thái IIb: $3.369,20\text{ ha} \times 27,81\text{ tấn C/ha} = \mathbf{93.697,45\text{ tấn C}}$ (tương đương $343.870\text{ tấn }\text{CO}_2$).
  • Diện tích trạng thái IIIa1: $2.612,14\text{ ha} \times 33,56\text{ tấn C/ha} = \mathbf{87.663,42\text{ tấn C}}$ (tương đương $321.725\text{ tấn }\text{CO}_2$).
  • Tổng trữ lượng 2 trạng thái: $\mathbf{181.360,87\text{ tấn C}}$, cố định tổng cộng $665.595\text{ tấn }\text{CO}_2$ tương đương.

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

  1. Phương pháp điều tra ô lồng vệ tinh cải tiến trên núi đá vôi: Nghiên cứu áp dụng thành công thiết kế ô tiêu chuẩn vệ tinh của Zoehrer (1980) chia nhỏ thành 4 ô thứ cấp $20\text{m} \times 25\text{m}$ định hướng linh hoạt (Tây, Bắc, Đông, Nam hoặc Tây Nam/Đông Nam khi vướng vách đá). Thiết kế này giảm thiểu $35%$ thời gian di chuyển nguy hiểm trên địa hình Karst hiểm trở so với ô vuông truyền thống $50\text{m} \times 40\text{m}$.
  2. Định lượng toàn diện 3 bể chứa các-bon trên mặt đất: Thay vì chỉ tính tầng cây gỗ như các nghiên cứu lâm sinh thông thường, đề tài chứng minh tầng thảm mục và cây bụi thảm tươi chiếm đến $48,60%$ tổng lượng $\text{C}$ ở trạng thái IIb và $38,23%$ ở trạng thái IIIa1. Việc bỏ qua hai bể chứa này sẽ gây đánh giá thấp nghiêm trọng tiềm năng tín chỉ các-bon của rừng phục hồi.
  3. So sánh tương quan với các công bố khoa học:
Công trình nghiên cứu Đối tượng & Địa bàn Trữ lượng Các-bon (tấn C/ha) Đương lượng $\text{CO}_2$ (tấn/ha) Ghi chú so sánh
Chu Văn Tâm (Nghiên cứu này, 2016) Rừng tự nhiên IIb & IIIa1, Núi đá Thần Sa 27,81 - 33,56 101,97 - 123,06 Đo chi tiết 3 bể chứa AGC trên nền Karst
Vũ Tấn Phương (2006) Rừng nghèo & Rừng trung bình (toàn quốc) 35,00 - 65,00 128,33 - 238,55 Rừng trên đất đồi có tầng đất dày hơn núi đá
Nguyễn Văn Dũng (2005) Rừng Thông mã vĩ & Keo lá tràm (Hà Nội) 62,50 - 122,00 229,38 - 447,74 Rừng trồng thuần loài mật độ cao
Brown et al. (1980 / 1986) Rừng nhiệt đới Đông Nam Á 25,00 - 215,00 91,75 - 789,05 Khoảng biến động cảnh quan vĩ mô

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

Kịch bản thương mại hóa Tín chỉ Các-bon và Chi trả Dịch vụ Môi trường Rừng (PFES)

Với tổng trữ lượng $665.595\text{ tấn }\text{CO}_2\text{e}$ được lưu giữ ổn định trên $5.981,34\text{ ha}$, khu vực Thần Sa có tiềm năng kinh tế - môi trường rất lớn khi tham gia thị trường các-bon tự nguyện (Voluntary Carbon Market - VCM) theo tiêu chuẩn VCS (Verified Carbon Standard) hoặc Plan Vivo:

                              MÔ HÌNH DÒNG TIỀN DỰ ÁN REDD+
+-----------------------------------------------------------------------------------+
| BẢO VỆ 5.981,34 HA RỪNG XÃ THẦN SA TRÁNH MẤT RỪNG (Baseline Deforestation: 1.2%/năm)|
+----------------------------------------+------------------------------------------+
                                         |
                                         v
+-----------------------------------------------------------------------------------+
| Lượng giảm phát thải tạo tín chỉ: ~7.987 tấn CO2e/năm                             |
| Giá bán tín chỉ tự nguyện dự kiến: 8.00 USD / tấn CO2e                           |
| Doanh thu hằng năm: 63.896 USD (~1,59 tỷ VNĐ / năm)                               |
+----------------------------------------+------------------------------------------+
                                         |
         +-------------------------------+-------------------------------+
         |                               |                               |
+--------v-------+              +--------v-------+              +--------v-------+
| HỖ TRỢ SINH KẾ |              | QUỸ TUẦN TRA   |              | GIÁM SÁT MRV & |
| 562 HỘ DÂN     |              | BẢO VỆ RỪNG    |              | ĐO ĐẾM ĐỊNH KỲ |
| (60% ~ 954 tr) |              | (25% ~ 397 tr) |              | (15% ~ 238 tr) |
+----------------+              +----------------+              +----------------+

Roadmap Triển khai Giám sát MRV (Measurement, Reporting, Verification)

  • Giai đoạn 1 (Tháng 1 - 3): Số hóa ranh giới tiểu khu và lô rừng trên hệ thống GIS Web, thiết lập mạng lưới 30 ô tiêu chuẩn định vị GPS cố định.
  • Giai đoạn 2 (Tháng 4 - 8): Tập huấn cho Hạt Kiểm lâm Võ Nhai và người dân địa phương quy trình thu thập mẫu thảm tươi, thảm mục và đo đếm $D_{1.3}$ phi phá hủy.
  • Giai đoạn 3 (Tháng 9 - 12): Tích hợp công nghệ viễn thám Sentinel-2 để theo dõi biến động chỉ số lá rậm rạp (Leaf Area Index - LAI) và cập nhật đường phát thải cơ sở (FREL).

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

Hạn chế kỹ thuật

  • Cỡ mẫu điều tra còn khiêm tốn: Do điều kiện địa hình núi đá dốc đứng và kinh phí khóa luận giới hạn, số lượng ô tiêu chuẩn mới dừng lại ở 6 OTC ($1,2\text{ ha}$ diện tích thực nghiệm).
  • Chưa đo đếm bể chứa dưới mặt đất: Bể sinh khối rễ cây ($BGB$) và các-bon hữu cơ trong đất đá vôi ($SOC$) chưa được lấy mẫu trực tiếp do hạn chế trang thiết bị khoan đất tầng sâu.

Hướng nâng cấp đề xuất

  • Ứng dụng Drone/UAV LiDAR: Quét đám mây điểm 3D cấu trúc tán rừng trên núi đá để trích xuất trực tiếp chiều cao cây ($H$) và thể tích thân cây cho toàn bộ lưu vực mà không cần leo trèo thực địa.
  • Xây dựng phương trình sinh khối bản địa: Thực hiện giải tích thân cây rụng tự nhiên của các loài ưu thế bản địa (Kháo, Vàng anh, Bứa, Sấu) để tối ưu hóa hệ số allometric thay vì dùng công thức chung của Brown et al. (1989).

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

  • Sinh viên & Giảng viên Lâm nghiệp: Tiếp cận bộ tài liệu hướng dẫn chuẩn về phương pháp lập ô tiêu chuẩn vệ tinh Zoehrer, kỹ thuật sấy mẫu phân tích sinh khối và công thức tính toán chỉ số IVI.
  • Kỹ sư GIS & Quản lý Tài nguyên Rừng: Ứng dụng mã nguồn phân tích sinh khối tự động và quy trình tích hợp lớp bản đồ hiện trạng rừng MapInfo phục vụ kiểm kê khí nhà kính.
  • Chủ rừng & Ban Quản lý Khu BTTN: Nắm giữ số liệu khoa học xác thực để lập hồ sơ phê duyệt chi trả Dịch vụ Môi trường Rừng (PFES) và thu hút các quỹ bảo tồn quốc tế.
  • Cộng đồng cư dân bản địa (Tày, Dao, Mông, Nùng): Tăng thu nhập ổn định thông qua các hợp đồng khoán bảo vệ rừng và chia sẻ lợi ích từ thị trường các-bon.

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

1. Yêu cầu thiết bị kỹ thuật tối thiểu để triển khai đo đếm sinh khối rừng là gì?

Cần chuẩn bị thước đo đường kính ngang ngực Richter 5m, thước đo độ cao Blume-Leiss (hoặc súng bắn góc Suunto Hypsometer), thiết bị định vị GPS cầm tay Garmin 64s, khung gỗ $1\text{m} \times 1\text{m}$ lấy mẫu thảm mục, cân điện tử hiện số độ nhạy $0,01\text{g}$, túi bảo quản mẫu bằng giấy kraft và tủ sấy đối lưu nhiệt độ $85^\circ\text{C}$.

2. Tại sao tỷ trọng các-bon thảm mục ở trạng thái IIb (31,81%) lại cao hơn trạng thái IIIa1 (20,65%)?

Ở trạng thái IIb (rừng phục hồi), tổ thành loài gồm nhiều cây tiên phong ưa sáng lá rộng có chu kỳ rụng lá nhanh và thảm tươi dễ gãy đổ, kết hợp độ tàn che thấp hơn ($70% - 75%$) làm tăng quá trình rụng cành nhánh. Ngược lại, ở trạng thái IIIa1, sinh khối tập trung chủ yếu vào tầng cây gỗ lớn ($61,77%$), tầng tán rậm rạp ($80% - 85%$) làm ẩm độ nền rừng cao, thúc đẩy vi sinh vật phân giải thảm mục nhanh hơn.

3. Phương trình sinh khối của Brown et al. (1989) có đảm bảo độ tin cậy cho rừng núi đá vôi?

Phương trình $Y = \exp(-2.134 + 2.530 \times \ln(D))$ được IPCC khuyến nghị áp dụng cho rừng mưa nhiệt đới ẩm có lượng mưa từ $1.500 - 4.000\text{ mm/năm}$. Mặc dù có sai số nhất định do đặc thù cây trên núi đá sinh trưởng chậm, đây là phương trình phi phá hủy chuẩn mực nhất được quốc tế công nhận khi chưa có phương trình chuyên biệt cho cây gỗ núi đá vôi Việt Nam.

4. Chi phí duy trì hệ thống giám sát trữ lượng các-bon định kỳ là bao nhiêu?

Chi phí đo đếm nhắc lại 5 năm/lần cho mạng lưới ô tiêu chuẩn cố định ước tính khoảng $30.000 - 50.000\text{ VNĐ/ha}$, chỉ chiếm chưa đầy $3% - 5%$ doanh thu thu được từ việc bán tín chỉ các-bon hoặc chi trả dịch vụ môi trường rừng.

5. Khả năng tích hợp dữ liệu nghiên cứu vào hệ thống báo cáo phát thải quốc gia?

Toàn bộ dữ liệu mật độ, sinh khối khô và hệ số chuyển đổi $CF=0,46$ hoàn toàn tương thích với hướng dẫn IPCC Tier 2, sẵn sàng tích hợp trực tiếp vào Báo cáo Cập nhật Hai năm một lần (BUR) và Đóng góp do quốc gia tự quyết định (NDC) của Việt Nam gửi Công ước Khung Liên Hợp Quốc về Biến đổi Khí hậu (UNFCCC).

Kết luận

Công trình nghiên cứu đã lượng hóa thành công năng lực tích lũy các-bon phần trên mặt đất của hai trạng thái rừng tự nhiên đặc trưng tại xã Thần Sa, Khu BTTN Thần Sa – Phượng Hoàng. Kết quả khẳng định:

  • Trạng thái IIb lưu giữ trung bình $27,81\text{ tấn C/ha}$ ($101,97\text{ tấn }\text{CO}_2\text{e/ha}$), tổng trữ lượng toàn xã đạt $93.697,45\text{ tấn C}$.
  • Trạng thái IIIa1 lưu giữ trung bình $33,56\text{ tấn C/ha}$ ($123,06\text{ tấn }\text{CO}_2\text{e/ha}$), tổng trữ lượng toàn xã đạt $87.663,42\text{ tấn C}$.
  • Tổng khối lượng các-bon cố định trên $5.981,34\text{ ha}$ rừng tự nhiên đạt $181.360,87\text{ tấn C}$ (tương đương $665.595\text{ tấn }\text{CO}_2$).

Đây là nguồn tư liệu khoa học then chốt chứng minh giá trị sinh thái to lớn của rừng đặc dụng núi đá vôi Thái Nguyên, tạo tiền đề vững chắc để địa phương thu hút các nguồn tài chính khí hậu quốc tế, triển khai dự án REDD+ và nâng cao sinh kế bền vững cho đồng bào các dân tộc thiểu số tại khu vực vùng đệm.