Giới thiệu dự án

Bối cảnh và thực trạng ngành

Biến đổi khí hậu toàn cầu đang diễn ra với tốc độ phức tạp, biểu hiện qua sự gia tăng nhiệt độ khí quyển, tan băng và dâng cao của mực nước biển. Theo báo cáo từ Ủy ban Liên chính phủ về Biến đổi Khí hậu (IPCC), phát thải khí nhà kính (KNK)—đặc biệt là $\text{CO}_2$—chiếm tỷ trọng lớn nhất trong các tác nhân gây hiệu ứng nhà kính. Trong đó, mất rừng và suy thoái rừng nhiệt đới đóng góp xấp xỉ 17–20% tổng lượng phát thải toàn cầu.

Trước thực trạng này, các công cụ tài chính và cơ chế quốc tế như Cơ chế Phát triển Sạch (CDM/AR-CDM) theo Nghị định thư Kyoto và Sáng kiến Giảm phát thải từ Mất rừng và Suy thoái rừng (REDD+) ra đời nhằm biến bể chứa carbon sinh thái thành nguồn tài chính bền vững. Tại Việt Nam, chính sách Chi trả Dịch vụ Môi trường Rừng (PFES) đang được luật hóa sâu rộng, đòi hỏi việc lượng hóa chính xác trữ lượng carbon cơ sở (carbon baseline) để phục vụ chứng nhận và giao dịch tín chỉ carbon.

+-------------------------------------------------------------------------------+
|                           BỂ CHỨA CARBON RỪNG TỰ NHIÊN                       |
+------------------------------------+------------------------------------------+
|  Sinh khối trên mặt đất (AGB)      |  Sinh khối dưới mặt đất & Đất (BGB & SOC)|
|  - Tầng cây gỗ sống (DBH >= 5cm)   |  - Hệ rễ cây gỗ                          |
|  - Tầng cây bụi, thảm tươi         |  - Hữu cơ trong đất                      |
|  - Thảm mục & vật rơi rụng         |                                          |
+------------------------------------+------------------------------------------+

Vấn đề nghiên cứu cụ thể

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) là khu vực có tính đa dạng sinh học cao với địa hình núi đá vôi karst đặc thù. Tuy nhiên, tại xã Nghinh Tường—vùng lõi khu bảo tồn:

  • Chưa có dữ liệu định lượng chuẩn xác về sinh khối trên mặt đất (Aboveground Biomass - AGB) và trữ lượng carbon tích lũy cho từng trạng thái rừng phục hồi (IIB) và rừng nghèo (IIIA1).
  • Thiếu các phương trình tương quan sinh khối và dữ liệu cấu trúc tầng tán phục vụ trực tiếp cho việc xây dựng kịch bản phát thải cơ sở (FREL/FRL) cấp xã/tiểu khu.
  • Áp lực từ sinh kế người dân địa phương dẫn tới nguy cơ suy thoái rừng nếu không có cơ chế bù đắp tài chính carbon minh bạch.

Mục tiêu dự án

  1. Xác định diện tích, các chỉ số cấu trúc lâm phần và Chỉ số Giá trị Quan trọng (Importance Value Index - IVI) của trạng thái rừng IIB và IIIA1 tại xã Nghinh Tường.
  2. Lượng hóa sinh khối khô trên mặt đất của 3 thành phần chính: tầng cây gỗ ($D_{1.3} \ge 5\text{ cm}$), tầng cây bụi – thảm tươi, và tầng thảm mục hữu cơ.
  3. Xác định hệ số tích lũy carbon thực tế và tính toán tổng trữ lượng carbon tích lũy ($tC$) cũng như lượng $\text{CO}_2$ tương đương cô lập được trên toàn bộ diện tích hai trạng thái rừng.
  4. Đánh giá nguy cơ suy giảm trữ lượng carbon và đề xuất giải pháp kỹ thuật quản trị rừng bền vững gắn liền với thị trường carbon.

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

Đề tài áp dụng phương pháp điều tra lâm học kết hợp thực nghiệm phân tích phòng thí nghiệm:

  • Thiết lập ô tiêu chuẩn (OTC) lồng nhau: Ứng dụng kỹ thuật bố trí cụm ô tiêu chuẩn vệ tinh theo Zoehrer (1980) cải tiến, diện tích $2.000\text{ m}^2$/OTC, chia nhỏ thành các ô thứ cấp $20 \times 25\text{ m}$ và ô dạng bản $1 \times 1\text{ m}$ nhằm thu thập số liệu không phá hủy tầng cây cao và lấy mẫu phá hủy tầng dưới tán/thảm mục.
  • Phương trình tương quan Allometric: Ứng dụng mô hình toán học Brown et al. (1989) chuẩn hóa cho rừng nhiệt đới ẩm để tính toán sinh khối cây gỗ từ đường kính ngang ngực ($D_{1.3}$).
  • Quy đổi trữ lượng Carbon: Ứng dụng hệ số chuyển đổi mặc định $0{,}46$ theo khuyến nghị của IPCC (2003) và nhân tỷ lệ chuyển đổi phân tử $\frac{44}{12} \approx 3{,}67$ để quy đổi sang $\text{CO}_2$ hấp thụ.
                                                                                    [Quy đổi CO2e (x3.67)]

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

  • Hoàn thành bộ cơ sở dữ liệu định lượng về mật độ cây ($N/\text{ha}$), tiết diện ngang ($G/\text{ha}$), trữ lượng gỗ ($M/\text{ha}$) và chỉ số IVI cho 6 OTC đại diện.
  • Xác định chính xác mật độ tích lũy carbon ($tC/\text{ha}$) cho từng trạng thái rừng với độ sai số ước tính dưới $10%$.
  • Lập bảng tổng hợp trữ lượng carbon quy đổi trên tổng quy mô $4.264{,}8\text{ ha}$ rừng IIB và IIIA1 của xã Nghinh Tường.

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

  • Không gian: Trạng thái rừng phục hồi (IIB) và rừng nghèo (IIIA1) phân bố tại xã Nghinh Tường, thuộc Khu Bảo tồn Thiên nhiên Thần Sa – Phượng Hoàng.
  • Bể chứa carbon: Giới hạn ở 3 bể chứa carbon trên mặt đất (AGB), không thực hiện đào xới xác định bể rễ ngầm (BGB) và bể carbon hữu cơ trong đất (SOC) nhằm bảo toàn nghiêm ngặt hiện trạng đất rừng đặc dụng.

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

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

Tiêu chí Phương pháp Viễn thám thuần túy (Remote Sensing) Phương pháp Khai thác Chặt hạ (Destructive Sampling) Giải pháp Đề xuất: OTC Lồng nhau + Mô hình Tương quan
Độ chính xác Trung bình (sai số 20–35% do tán che phức tạp) Rất cao (sai số < 2%) Cao (sai số 5–8% theo chuẩn IPCC Tier 2)
Bảo tồn hệ sinh thái Hoàn toàn không xâm phạm Phá hủy nghiêm trọng quần xã sinh vật Không phá hủy tầng cây gỗ, chỉ lấy mẫu vi mô tầng mặt
Chi phí & Nhân lực Cần ảnh viễn thám độ phân giải cao đắt tiền Chi phí khai thác, vận chuyển và nhân công cực lớn Tối ưu, phù hợp điều kiện nghiên cứu địa phương
Tính khả thi Khó áp dụng độc lập cho rừng karst tầng tán dày Bị cấm trong phân khu bảo vệ nghiêm ngặt khu bảo tồn Khả thi cao, tuân thủ quy chế rừng đặc dụng

Yêu cầu hệ thống theo mô hình MoSCoW

  • Must Have: Đo đếm đầy đủ $D_{1.3} \ge 5\text{ cm}$, tính toán IVI từng loài; xác định khối lượng tươi/khô của tầng thảm tươi và thảm mục qua sấy $85^\circ\text{C}$ trong $48\text{ giờ}$; tính trữ lượng carbon trung bình trên đơn vị hecta.
  • Should Have: Ứng dụng phần mềm GIS (MapInfo) để bóc tách hiện trạng ranh giới 3 loại rừng và diện tích từng tiểu khu.
  • Could Have: Xây dựng phương trình allometric riêng cho các loài ưu thế địa phương (Táu, Sổ, Chẹo, Dẻ).
  • Won't Have: Đo đạc biến động dòng khí nhà kính thời gian thực qua tháp Eddy Covariance.

Thiết kế hệ thống kỹ thuật

Công cụ và Công nghệ thực hiện

  • Phần mềm xử lý không gian & số liệu: MapInfo Professional v15.0 (biên tập bản đồ địa hình và trạng thái rừng), Microsoft Excel 2016 / R Studio v4.2 (xử lý mô hình hồi quy và thống kê toán học).
  • Thiết bị ngoại nghiệp: Thước đo đường kính cây Richter, thước đo góc/cao bách phân Blume-Leiss, GPS Garmin 64s độ chính xác $\pm 3\text{ m}$.
  • Thiết bị nội nghiệp: Tủ sấy đối lưu cưỡng bức Memmert UN55 (Đức) dải nhiệt độ đến $300^\circ\text{C}$, Cân phân tích điện tử chuẩn Ohaus NV222 (độ chia $0{,}01\text{ g}$).

Phương pháp luận và Quy trình triển khai

  1. Thiết lập Ô mẫu vệ tinh Zoehrer (1980):
    • Mỗi OTC $2.000\text{ m}^2$ gồm 4 ô thứ cấp $20 \times 25\text{ m}$ theo 4 hướng Đông, Tây, Nam, Bắc.
    • Trong ô thứ cấp đầu tiên, thiết lập 01 ô cấp 1 ($10 \times 20\text{ m} = 200\text{ m}^2$) đo cây gỗ.
    • Bố trí 05 ô dạng bản ($1 \times 1\text{ m}$) tại 4 góc và tâm ô cấp 1 để thu thập toàn bộ sinh khối tươi của cây bụi, thảm tươi và thảm mục.
  2. Kế hoạch quản trị rủi ro:
    • Rủi ro địa hình karst dốc đứng, vách đá: Thay thế hướng đặt ô tiêu chuẩn thứ cấp bằng các hướng phụ (Đông Nam, Tây Nam) khi gặp trở ngại vật lý.
    • Rủi ro ẩm mốc mẫu sinh vật: Đóng gói kín trong túi giấy chuyên dụng, vận chuyển về phòng thí nghiệm trong vòng $12\text{ giờ}$ và đưa ngay vào chu trình sấy khử ẩm.

Implementation và kết quả

Quy trình phân tích số liệu và thuật toán lõi

Dữ liệu đo đạc thực địa được chuẩn hóa qua các công thức toán - sinh học lâm nghiệp chuyên sâu:

  1. Công thức Mật độ ($N$) và Tiết diện ngang ($G$): $$N = \frac{n}{S} \times 10.000 \quad (\text{cây/ha})$$ $$G = \sum_{i=1}^n \frac{\pi \times D_i^2}{40.000} \quad (\text{m}^2/\text{ha})$$

  2. Chỉ số Giá trị Quan trọng (Importance Value Index - IVI): $$IVI_i (%) = \frac{A_i% + D_i% + RF_i%}{3}$$ Trong đó: $A_i%$ là độ phong phú tương đối, $D_i%$ là độ ưu thế tương đối theo tiết diện ngang, $RF_i%$ là tần suất xuất hiện tương đối.

  3. Mô hình Sinh khối Cây gỗ Allometric (Brown et al., 1989): $$Y = \exp\left[-2{,}134 + 2{,}530 \times \ln(D_{1.3})\right] \quad (\text{kg/cây})$$

  4. Sinh khối khô Tầng dưới tán và Thảm mục ($DW$): $$DW = F_w \times \frac{D_{w_\text{sub}}}{F_{w_\text{sub}}} \quad (\text{kg/m}^2 \text{ hoặc tấn/ha})$$

  5. Trữ lượng Carbon và $\text{CO}_2$ quy đổi: $$W_{\text{Carbon}} = 0{,}46 \times DW_{\text{Tổng}} \quad (\text{tấn C/ha})$$ $$\text{CO}{2\text{e}} = W{\text{Carbon}} \times \frac{44}{12} \approx W_{\text{Carbon}} \times 3{,}667 \quad (\text{tấn }\text{CO}_2/\text{ha})$$

Kịch bản xử lý dữ liệu (Python Script Reference)

import numpy as np
import pandas as pd

def calculate_forest_carbon(df_trees, fw_understory, ratio_understory, fw_litter, ratio_litter, area_ha=1.0):
    """
    Tính toán sinh khối, trữ lượng Carbon và CO2 tương đương cho lâm phần
    """
    # 1. Sinh khối cây gỗ theo Brown et al. (1989)
    df_trees['Biomass_kg'] = np.exp(-2.134 + 2.530 * np.log(df_trees['D13_cm']))
    agb_trees_ton_ha = (df_trees['Biomass_kg'].sum() / 1000.0) * (10000.0 / (area_ha * 10000.0))
    
    # 2. Sinh khối tầng dưới tán và thảm mục (tấn/ha)
    agb_understory_ton_ha = (fw_understory * ratio_understory) * 10.0  # kg/m2 -> ton/ha
    agb_litter_ton_ha = (fw_litter * ratio_litter) * 10.0
    
    # 3. Tổng sinh khối khô AGB
    total_agb = agb_trees_ton_ha + agb_understory_ton_ha + agb_litter_ton_ha
    
    # 4. Trữ lượng Carbon (IPCC 2003: hệ số 0.46)
    c_trees = agb_trees_ton_ha * 0.46
    c_understory = agb_understory_ton_ha * 0.46
    c_litter = agb_litter_ton_ha * 0.46
    total_carbon = total_agb * 0.46
    
    # 5. Quy đổi CO2 tương đương
    total_co2e = total_carbon * (44.0 / 12.0)
    
    return {
        'AGB_Trees_t_ha': round(agb_trees_ton_ha, 2),
        'AGB_Understory_t_ha': round(agb_understory_ton_ha, 2),
        'AGB_Litter_t_ha': round(agb_litter_ton_ha, 2),
        'Total_AGB_t_ha': round(total_agb, 2),
        'Carbon_Trees_tC_ha': round(c_trees, 2),
        'Total_Carbon_tC_ha': round(total_carbon, 2),
        'Total_CO2e_t_ha': round(total_co2e, 2)
    }

# Dữ liệu kiểm chứng mẫu cho Trạng thái IIIA1 - OTC 04
sample_trees = pd.DataFrame({'D13_cm': [12.5, 18.2, 8.4, 22.1, 15.6, 9.8, 14.2, 26.5]})
result = calculate_forest_carbon(sample_trees, fw_understory=1.32, ratio_understory=0.46, 
                                 fw_litter=1.45, ratio_litter=0.446, area_ha=0.08)

Kết quả nghiên cứu thực nghiệm

1. Đặc điểm cấu trúc và tổ thành loài lâm phần

  • Trạng thái rừng phục hồi (IIB): Mật độ dao động từ $315 - 435\text{ cây/ha}$, tiết diện ngang trung bình đạt $3{,}40 - 4{,}33\text{ m}^2/\text{ha}$, trữ lượng cây gỗ đứng đạt $17{,}04 - 23{,}71\text{ m}^3/\text{ha}$. Tổ thành loài chiếm ưu thế sinh thái ($IVI > 5%$) gồm: Táu ($IVI = 7{,}39%$), Sổ ($IVI = 4{,}56 - 9{,}61%$), Kháo, Vạng trứng.
  • Trạng thái rừng nghèo (IIIA1): Mật độ đạt từ $455 - 470\text{ cây/ha}$, tiết diện ngang tăng vượt bậc lên $8{,}51 - 9{,}68\text{ m}^2/\text{ha}$, trữ lượng gỗ từ $49{,}65 - 58{,}51\text{ m}^3/\text{ha}$. Các loài ưu thế tầng tán gồm: Chẹo ($IVI = 11{,}80%$), Dẻ ($IVI = 10{,}96%$), Kẹn, Đinh lăng rừng.

2. Định lượng sinh khối trên mặt đất (AGB)

Trạng thái Ô Tiêu Chuẩn Cây gỗ ($t/\text{ha}$) Cây bụi - Thảm tươi ($t/\text{ha}$) Thảm mục ($t/\text{ha}$) Tổng Sinh khối ($t/\text{ha}$)
IIB (Phục hồi) OTC 01 (Nà Hấu) 24,21 10,75 14,90 49,86
OTC 02 (Thượng Lương) 18,50 10,60 15,46 44,56
OTC 03 (Nà Hấu) 25,50 9,58 15,03 50,11
Trung bình IIB 22,74 10,31 15,13 48,18
IIIA1 (Nghèo) OTC 04 (Thượng Lương) 58,68 13,20 14,06 85,94
OTC 05 (Hạ Lương) 63,04 13,80 15,65 92,49
OTC 06 (Thượng Lương) 53,94 15,50 12,24 81,68
Trung bình IIIA1 58,55 14,17 13,98 86,70

3. Trữ lượng Carbon tích lũy và quy đổi $\text{CO}_2$

+-------------------------------------------------------------------------------+
|             CƠ CẤU BỂ CHỨA CARBON GIỮA HAI TRẠNG THÁI RỪNG                    |
+-------------------------------------------------------------------------------+
| Trạng thái IIB:                                                               |
| [████████████████████ 47.2%] Cây gỗ                                           |
| [███████████ 21.4%] Thảm tươi                                                 |
| [███████████████ 31.4%] Thảm mục                                              |
|                                                                               |
| Trạng thái IIIA1:                                                             |
| [████████████████████████████████ 67.5%] Cây gỗ                              |
| [███████ 16.3%] Thảm tươi                                                     |
| [████████ 16.2%] Thảm mục                                                     |
+-------------------------------------------------------------------------------+
  • Trạng thái IIB:
    • Mật độ carbon trung bình: $\mathbf{22{,}16\text{ tấn C/ha}}$ (tương đương hấp thụ $\mathbf{81{,}33\text{ tấn CO}_{2\text{e}}/\text{ha}}$).
    • Tỷ trọng phân bổ: Tầng cây gỗ chiếm $47{,}2%$, thảm mục chiếm $31{,}4%$, thảm tươi chiếm $21{,}4%$.
  • Trạng thái IIIA1:
    • Mật độ carbon trung bình: $\mathbf{39{,}88\text{ tấn C/ha}}$ (tương đương hấp thụ $\mathbf{146{,}36\text{ tấn CO}_{2\text{e}}/\text{ha}}$).
    • Tỷ trọng phân bổ: Tầng cây gỗ chiếm ưu thế tuyệt đối với $67{,}5%$, thảm tươi chiếm $16{,}3%$, thảm mục chiếm $16{,}2%$.

4. Tổng trữ lượng Carbon trên quy mô toàn xã Nghinh Tường

Trạng thái Rừng Diện tích ($ha$) Mật độ Carbon ($tC/ha$) Tổng Carbon ($tC$) Quy đổi $\text{CO}_{2\text{e}}$ Hấp thụ ($t\text{CO}_2$)
Trạng thái IIB 2.743,9 22,16 60.804,62 223.153,00
Trạng thái IIIA1 1.520,9 39,88 60.653,49 222.598,30
Tổng cộng 4.264,8 — 121.458,11 445.751,30

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

Điểm mới khoa học và kỹ thuật

  1. Chuẩn hóa quy trình quan trắc carbon trên địa hình karst: Khắc phục nhược điểm định vị sai số lớn trên núi đá vôi bằng cách kết hợp trắc địa vệ tinh phân giải cao và phương pháp bố trí ô mẫu lồng đa tầng theo Zoehrer biến tính.
  2. Lượng hóa chi tiết thành phần thảm mục và thảm tươi: Không bỏ qua bể chứa thứ cấp (vốn chiếm tới $52{,}8%$ trữ lượng carbon ở rừng phục hồi IIB), tạo nên bộ tham số đầy đủ cho mô hình IPCC Tier 2 tại miền núi phía Bắc.
  3. Cơ sở dữ liệu thực chứng cho phân khu bảo tồn: Cung cấp bộ thông số cấu trúc lâm phần ($D_{1.3}, H_{vn}, G, M, IVI$) cho các loài cây bản địa có giá trị bảo tồn cao.

So sánh với các nghiên cứu trước đây

Nghiên cứu & Địa bàn Đối tượng nghiên cứu Mật độ Carbon Cây gỗ ($tC/\text{ha}$) Tổng Carbon AGB ($tC/\text{ha}$) Ghi chú & Khác biệt
Dương Văn Dũng (Nghiên cứu này - Võ Nhai) Rừng IIB (Phục hồi)
Rừng IIIA1 (Nghèo)		 10,46
26,93		 22,16
39,88		 Rừng núi đá vôi, phân tích đủ 3 bể thành phần mặt đất
Nghiên cứu tại Tam Đảo (Vĩnh Phúc) Nông lâm kết hợp (Keo/Thông) 7,50 – 12,30 16,40 – 21,80 Rừng trồng thuần/hỗn loài đồi đất, sinh khối thảm mục thấp
Nghiên cứu tại Thừa Thiên Huế Trảng cỏ, cây bụi 0,80 – 2,10 7,30 – 11,20 Chưa hình thành cấu trúc tầng tán cây gỗ
Nghiên cứu toàn quốc (Viện ĐTQH Rừng) Rừng nghèo tự nhiên 25,00 – 35,00 45,00 – 60,00 Trung bình đa dạng lập địa, bao gồm cả rừng đất ẩm 
Hiệu quả tối ưu hóa:
- Giảm thiểu 100% việc chặt hạ cây gỗ nhờ mô hình tương quan sinh khối Allometric.
- Giảm 45% thời gian thu thập mẫu ngoại nghiệp thông qua kỹ thuật lấy mẫu lồng 200m2.

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

Kịch bản thương mại hóa Tín chỉ Carbon (VCS/REDD+)

Với tổng lượng hấp thụ tích lũy trên mặt đất đạt $\mathbf{445.751{,}30\text{ tấn CO}_{2\text{e}}$, xã Nghinh Tường sở hữu tiềm năng rất lớn để tham gia thị trường carbon tự nguyện (VCM - Voluntary Carbon Market) theo tiêu chuẩn Tiêu chuẩn Carbon Đảm bảo (VCS - Verified Carbon Standard) của Verra hoặc chương trình REDD+ cấp tỉnh.

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

  • Doanh thu tiềm năng từ Tín chỉ Carbon: Với mức giá carbon trung bình thận trọng từ $5 - 10\text{ USD/tấn CO}_{2\text{e}}$, tổng giá trị bể carbon mặt đất hiện hữu ước đạt 2,22 - 4,45 triệu USD.
  • Chính sách Chi trả Dịch vụ Môi trường Rừng (PFES): Dữ liệu là căn cứ pháp lý để Ban Quản lý Khu bảo tồn Thần Sa – Phượng Hoàng và UBND xã Nghinh Tường phân bổ nguồn thu PFES từ các cơ sở thủy điện và nhà máy nước ngầm về các tổ đội nhận khoán bảo vệ rừng tại 12 xóm bản.
+-------------------------------------------------------------------------------+
|                    LỘ TRÌNH TRIỂN KHAI THỰC HÓA GIÁ TRỊ CARBON                |
+-------------------------------------------------------------------------------+
| [Giai đoạn 1: Số hóa & Hồ sơ]                                                |
|   - Tích hợp dữ liệu OTC vào WebGIS Khu bảo tồn Thần Sa                      |
|   - Xây dựng đường phát thải cơ sở (FREL)                                     |
|                                                                               |
| [Giai đoạn 2: Quản lý & Chi trả PFES]                                        |
|   - Giao khoán quản lý bảo vệ 4.264,8 ha cho cộng đồng thôn xóm               |
|   - Giải ngân chi trả dịch vụ môi trường rừng gắn với tỷ lệ che phủ          |
|                                                                               |
| [Giai đoạn 3: Đăng ký Tín chỉ Carbon Quốc tế]                                |
|   - Đánh giá thẩm định bên thứ ba (Validation / Verification body)            |
|   - Phát hành tín chỉ carbon thương mại trên sàn giao dịch quốc tế            |
+-------------------------------------------------------------------------------+

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

Hạn chế kỹ thuật

  1. Thiếu bể chứa ngầm: Chưa tiến hành đo đạc trực tiếp Bể sinh khối rễ dưới mặt đất (BGB) và Carbon hữu cơ trong đất (Soil Organic Carbon - SOC) do yêu cầu nghiêm ngặt về bảo tồn đất rừng đặc dụng.
  2. Dung lượng mẫu: Quy mô 6 ô tiêu chuẩn phân bố tại 3 xóm (Nà Hấu, Thượng Lương, Hạ Lương) phản ánh mẫu đại diện tốt cho xã nhưng cần tăng dày lưới ô mẫu khi mở rộng quy mô toàn lâm phận Khu bảo tồn ($19.914\text{ ha}$).
  3. Mô hình sinh khối chung: Đang sử dụng phương trình của Brown (1989) thiết kế chung cho rừng nhiệt đới, chưa xây dựng được phương trình allometric bản địa đặc thù cho các loài cây karst như Táu, Chẹo.

Định hướng mở rộng

  • Ứng dụng công nghệ viễn thám hiện đại: Tích hợp dữ liệu LiDAR quét tán rừng kết hợp ảnh vệ tinh đa phổ Sentinel-2 / Landsat 9 để xây dựng bản đồ số hóa sinh khối 3D độ phân giải $10\text{ m}$.
  • Quan trắc động thái sinh trưởng: Thiết lập hệ thống ô mẫu định vị cố định (Permanent Sampling Plots - PSP) để theo dõi tốc độ tăng trưởng sinh khối ($tC/\text{ha}/\text{năm}$) hàng năm.

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

+-------------------+-----------------------------------------------------------+
| Nhóm đối tượng    | Giá trị và Lợi ích cụ thể                                 |
+-------------------+-----------------------------------------------------------+
| Sinh viên & Học   | Nắm vững phương pháp luận đo đếm sinh khối rừng, kỹ thuật |
| viên Lâm nghiệp   | lập ô tiêu chuẩn Zoehrer và thuật toán allometric chuẩn.  |
+-------------------+-----------------------------------------------------------+
| Ban Quản lý Rừng  | Có bộ dữ liệu số hóa chính xác để xây dựng phương án quản |
| & Kiểm lâm        | lý rừng bền vững và bảo vệ nghiêm ngặt các diện tích IIIA1|
+-------------------+-----------------------------------------------------------+
| Doanh nghiệp      | Căn cứ khoa học minh bạch để thẩm định dự án tín chỉ     |
| Đầu tư Carbon     | carbon (Carbon Offsetting) theo chuẩn quốc tế.            |
+-------------------+-----------------------------------------------------------+
| Cộng đồng cư dân  | Nâng cao thu nhập bền vững thông qua tiền khoán bảo vệ    |
| bản địa           | rừng từ quỹ dịch vụ môi trường rừng và dự án carbon.     |
+-------------------+-----------------------------------------------------------+

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

1. Yêu cầu kỹ thuật tối thiểu để tiến hành điều tra trữ lượng carbon rừng là gì?

Cần trang bị tối thiểu thiết bị định vị GPS chuyên dụng, thước đo chu vi/đường kính thân cây, thước đo cao Blume-Leiss, dụng cụ khoanh mẫu dạng bản $1\text{ m}^2$, cân điện tử chính xác $0{,}01\text{ g}$ và tủ sấy sấy mẫu đạt chuẩn $85 - 105^\circ\text{C}$. Nhân sự cần được đào tạo quy chuẩn phân định loài cây và kỹ thuật lập ô tiêu chuẩn.

2. Tại sao tầng thảm mục ở trạng thái rừng IIB lại chiếm tỷ lệ carbon cao (31,4%)?

Ở trạng thái rừng phục hồi IIB, tầng cây gỗ chưa khép tán hoàn toàn và đang trong giai đoạn diễn thế mạnh, mật độ cây nhỏ nhiều dẫn đến lượng rơi rụng vật lý (cành khô, lá rụng) tích tụ trên bề mặt núi đá vôi phân hủy chậm trong mùa khô, tạo nên bể thảm mục lớn tạm thời trước khi khoáng hóa hoàn toàn thành mùn đất.

3. Phương trình allometric của Brown (1989) có độ tin cậy như thế nào đối với rừng núi đá vôi?

Mô hình Brown et al. (1989) ($Y = \exp[-2{,}134 + 2{,}530 \times \ln(D)]$) là phương trình thực nghiệm quốc tế được IPCC khuyến nghị sử dụng cho cấp độ Tier 1 và Tier 2 trong điều kiện thiếu phương trình bản địa, với hệ số tương quan $R^2 > 0{,}88$ đối với rừng mưa nhiệt đới ẩm hỗn loài lá rộng thường xanh.

4. Chi phí để chuyển hóa kết quả nghiên cứu thành tín chỉ carbon thương mại là bao nhiêu?

Chi phí phát triển dự án carbon (Project Design Document - PDD) bao gồm: khảo sát mở rộng lưới ô tiêu chuẩn, thẩm định độc lập bởi đơn vị được chỉ định (VVB), và phí đăng ký tài khoản trên sàn giao dịch Verra/Gold Standard thường dao động từ $40.000 - 80.000\text{ USD}$ cho toàn bộ vòng đời thẩm định ban đầu.

5. Làm thế nào để bảo toàn trữ lượng carbon rừng trước áp lực sinh kế của người dân địa phương?

Cần kết hợp hài hòa giữa giải pháp hành chính (tuần tra kiểm lâm liên ngành) và kinh tế (chia sẻ lợi ích tài chính từ tiền dịch vụ môi trường rừng PFES, hỗ trợ phát triển mô hình nông lâm kết hợp dưới tán rừng tại vùng đệm).

Kết luận

Đề tài tốt nghiệp của tác giả Dương Văn Dũng đã hoàn thành toàn diện các mục tiêu nghiên cứu đề ra:

  • Thực chứng khoa học vững chắc: Đã lượng hóa chính xác sinh khối và trữ lượng carbon tích lũy trên mặt đất của trạng thái rừng IIB ($22{,}16\text{ tấn C/ha}$) và trạng thái rừng IIIA1 ($39{,}88\text{ tấn C/ha}$) tại xã Nghinh Tường, Khu Bảo tồn Thiên nhiên Thần Sa – Phượng Hoàng.
  • Giá trị kinh tế và môi trường: Toàn bộ diện tích $4.264{,}8\text{ ha}$ của hai trạng thái rừng đang lưu giữ $121.458{,}11\text{ tấn Carbon}$, tương đương khả năng cô lập $445.751{,}30\text{ tấn CO}_{2\text{e}}$.
  • Tính ứng dụng cao: Kết quả nghiên cứu là luận cứ kỹ thuật quan trọng phục vụ trực tiếp cho chính sách Chi trả Dịch vụ Môi trường Rừng (PFES) tại tỉnh Thái Nguyên, đồng thời mở ra tiềm năng thương mại hóa tín chỉ carbon rừng theo cơ chế REDD+ và thị trường carbon quốc tế.