Giới thiệu dự án

Sự biến đổi khí hậu (BĐKH) toàn cầu do gia tăng nồng độ khí nhà kính (KNK/GHG - Greenhouse Gas) trong khí quyển đang là một trong những thách thức sinh thái và kinh tế lớn nhất của thế kỷ 21. Theo Báo cáo Đánh giá của Ban Liên chính phủ về Biến đổi Khí hậu (IPCC - Intergovernmental Panel on Climate Change), các hoạt động chuyển đổi mục đích sử dụng đất và suy thoái rừng nhiệt đới phát thải khoảng 1,5 tỷ tấn carbon mỗi năm, chiếm gần 20% tổng lượng phát thải khí nhà kính toàn cầu – con số vượt qua cả tổng phát thải từ ngành giao thông vận tải thế giới. Đứng trước sức ép này, các cơ chế tài chính carbon quốc tế như Cơ chế Phát triển Sạch (CDM - Clean Development Mechanism) theo Nghị định thư Kyoto và chương trình Giảm phát thải từ mất rừng và suy thoái rừng (REDD+ / Reducing Emissions from Deforestation and Forest Degradation) đã xác lập vai trò sống còn của việc định lượng, giám sát và thương mại hóa trữ lượng carbon rừng.

Tại Việt Nam, phần lớn các chương trình kiểm kê carbon rừng quốc gia và địa phương chỉ tập trung vào sinh khối trên mặt đất của tầng cây gỗ sống (AGB - Above-Ground Biomass), trong khi tầng vật rơi rụng (litterfall - bao gồm lá già, cành khô, hoa, quả và thảm mục đang phân giải) thường bị bỏ qua hoặc ước lượng sơ sài. Trên thực tế, vật rơi rụng là một bể chứa carbon (carbon pool) trung gian vô cùng năng động, đóng vai trò bản lề trong chu trình dinh dưỡng khoáng và lưu giữ carbon ngắn hạn của hệ sinh thái rừng nhiệt đới. Việc thiếu dữ liệu định lượng cụ thể ở trạng thái rừng phục hồi sau nương rẫy hoặc khai thác (trạng thái IIA theo phân loại lâm nghiệp Việt Nam) dẫn đến sự sai lệch từ 12% đến 18% trong việc xây dựng đường carbon cơ sở (carbon baseline) và thẩm định tín chỉ carbon.

Đề tài "Nghiên cứu khả năng tích lũy $\text{CO}_2$ của vật rơi rụng dưới tán rừng phục hồi (IIA) tại Xã La Bằng, Huyện Đại Từ, Tỉnh Thái Nguyên" được thực hiện nhằm giải quyết khoảng trống học thuật và thực tiễn này với các mục tiêu cụ thể:

  1. Xác định đặc điểm cấu trúc lâm phần, mật độ tầng cây cao và tầng cây tái sinh của trạng thái rừng phục hồi IIA tại xã La Bằng.
  2. Định lượng sinh khối tươi và sinh khối khô của tầng vật rơi rụng, phân định chi tiết theo hai hợp phần: cành rụng và lá/hoa/quả rụng.
  3. Xác định tổng lượng carbon ($C$) tích lũy và quy đổi lượng $\text{CO}2$ hấp thụ tương đương ($\text{CO}{2\text{-eq}}$) trên cơ sở hướng dẫn của IPCC.
  4. Xây dựng mô hình ước tính giá trị thương mại carbon từ bể vật rơi rụng theo khung giá thị trường carbon tự nguyện quốc tế.

Nghiên cứu áp dụng phương pháp điều tra thực địa kết hợp thu mẫu ô dạng bản (nested subplots) và phân tích giải tích sấy khô tiêu chuẩn trong phòng thí nghiệm. Phạm vi nghiên cứu tập trung vào 12 ô tiêu chuẩn ($2.500\text{ m}^2/\text{ô}$) trên diện tích 422 ha rừng sản xuất và 1.096 ha rừng đặc dụng thuộc vùng đệm Vườn Quốc gia Tam Đảo tại xã La Bằng, trong khung thời gian từ tháng 01/2015 đến tháng 05/2015.

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

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

Phương pháp xác định trữ lượng carbon hệ sinh thái rừng hiện nay tồn tại nhiều cách tiếp cận với các ưu và nhược điểm khác nhau:

Tiêu chí so sánh Phương pháp chặt hạ toàn phần (Destructive Sampling) Phương pháp Allometric không đo thảm mục Phương pháp ô dạng bản phân tích thảm mục (Đề tài áp dụng)
Độ chính xác bể thảm mục Rất cao ($> 95%$) Bị bỏ qua ($0%$) Rất cao ($90 - 95%$)
Tác động môi trường Phá hủy hoàn toàn cấu trúc lâm phần Không phá hủy Thu mẫu tối thiểu ($1\text{ m}^2/\text{ôDB}$), không ảnh hưởng rừng
Chi phí & Nhân lực Cực kỳ tốn kém, phức tạp Chi phí thấp nhưng sai số cao Chi phí tối ưu, dễ nhân rộng theo chuẩn MRV
Tính khả thi cho REDD+ Không khả thi trên diện rộng Thiếu bể chứa carbon theo chuẩn IPCC Tier 2 Tương thích hoàn toàn với tiêu chuẩn IPCC và CDM

Nhiều công trình trước đây như nghiên cứu của Bảo Huy (2009) tại Tây Nguyên chỉ giới hạn ở việc lập phương trình sinh khối cây gỗ sống; nghiên cứu của Ngô Đình Quế (2006) tập trung vào rừng trồng thuần loài (Acacia, Pinus); hoặc Vũ Tấn Phương (2006) khảo sát tầng cây bụi thảm tươi. Do đó, nghiên cứu này thu hẹp khoảng trống về cấu trúc thảm mục trong rừng tự nhiên nhiệt đới hỗn loài phục hồi.

Yêu cầu kỹ thuật và phân loại MoSCoW cho phương pháp nghiên cứu:

  • Must have (Bắt buộc): Lập ô tiêu chuẩn định vị GPS; thu thập toàn bộ vật rơi rụng ($D < 5\text{ cm}, L < 50\text{ cm}$) trong 5 ô dạng bản/ô tiêu chuẩn; sấy khô ở $75 - 80^\circ\text{C}$ đến khối lượng không đổi; áp dụng hệ số chuyển đổi carbon $CF = 0.5$.
  • Should have (Nên có): Phân chia chi tiết thành phần cành rụng độc lập với lá/hoa/quả; phân tích biến động sinh khối theo 3 vị trí địa hình (chân núi, sườn núi, đỉnh núi).
  • Could have (Có thể mở rộng): Kết nối dữ liệu với phần mềm chuyên dụng như CO2FIX V3.1 để mô phỏng động thái tích lũy dài hạn.
  • Won't have (Không thực hiện): Không đo lường bể carbon rễ ngầm dưới đất (BGB) và hàm lượng carbon hữu cơ trong đất sâu $> 30\text{ cm}$.

Thiết kế hệ thống

Hệ thống quy trình lấy mẫu, xử lý dữ liệu và tính toán carbon được thiết kế khép kín từ thực địa đến phòng thí nghiệm:

Technology Stack và Thiết bị tiêu chuẩn:

  • Thiết bị đo đạc thực địa: Thước dây sợi thủy tinh 50m, thước kính lùi đo đường kính thân ($D_{1.3}$), máy đo cao Blume-Leiss / Suunto Clinometer, la bàn địa chất Harbin DQL-8, định vị cầm tay GPS Garmin eTrex 20x.
  • Thiết bị phân tích phòng thí nghiệm: Tủ sấy đối lưu cưỡng bức Memmert UN55 (kiểm soát nhiệt độ $\pm 0.5^\circ\text{C}$), cân phân tích điện tử 3 số lẻ Sartorius ($0,001\text{g}$), cân kỹ thuật Ohaus ($0,01\text{g}$), bình hút ẩm hạt silica gel.
  • Công cụ tính toán & Xử lý số liệu: Microsoft Excel 2016 (Data Modeling & Statistical Analysis Toolpak), kiểm tra đối chiếu bằng thư viện Python pandas và scipy.stats (v3.8+), hỗ trợ kiểm chứng lý thuyết qua phần mềm CO2FIX V3.1.

Methodology

  • Phương pháp luận: Áp dụng Hướng dẫn thực hành tốt của IPCC (2006 IPCC Guidelines for National Greenhouse Gas Inventories - Volume 4: AFOLU) cho cấp độ Tier 2.
  • Kế hoạch triển khai (Project Timeline):
    • Tháng 01/2015: Thu thập bản đồ địa hình, khảo sát ngoại nghiệp, phân định ranh giới tiểu khu và trạng thái rừng IIA tại La Bằng.
    • Tháng 02/2015 – Tháng 03/2015: Thiết lập 12 ô tiêu chuẩn ngoại nghiệp ($50\text{m} \times 50\text{m}$), đo đếm $D_{1.3}$, $H_{vn}$, $D_t$, điều tra cây tái sinh và thu gom vật rơi rụng trong 60 ô dạng bản ($1\text{m} \times 1\text{m}$).
    • Tháng 04/2015: Xử lý giải tích phòng thí nghiệm (băm mẫu, cân khối lượng tươi, vận hành tủ sấy $75 - 80^\circ\text{C}$, cân lặp lại ở các mốc thời gian 2h, 4h, 6h, 8h).
    • Tháng 05/2015: Tổng hợp cơ sở dữ liệu, phân tích phương sai, tính toán hệ số quy đổi $\text{CO}_2$, định giá thương mại và hoàn thiện báo cáo khoa học.
  • Đánh giá rủi ro & Biện pháp giảm thiểu:
    • Rủi ro mất độ ẩm tự nhiên khi vận chuyển: Cân khối lượng tươi ngay tại hiện trường bằng cân kỹ thuật di động; đóng gói kín trong túi zip polyetylen dán nhãn niêm phong.
    • Rủi ro sai số sấy chưa kiệt: Thiết lập quy tắc kiểm tra khối lượng sau mỗi 2 giờ sấy; chỉ công nhận khối lượng khô khi 3 lần cân liên tiếp có độ lệch nhỏ hơn $0,01\text{g}$.

Implementation và kết quả

Development process

Quá trình định lượng cấu trúc lâm phần và vật rơi rụng được mô hình hóa qua các công thức toán học và thuật toán giải tích:

  1. Chỉ số giá trị quan trọng (Importance Value Index - IVI): $$IVI% = \frac{N% + G%}{2}$$ Trong đó:

    • $N%$: Tỷ lệ phần trăm số cá thể của một loài so với tổng số cây điều tra ($N% = \frac{n_i}{N} \times 100$).
    • $G%$: Tỷ lệ phần trăm tiết diện ngang thân cây ở vị trí $1,3\text{ m}$ của loài so với tổng tiết diện ngang ($G% = \frac{g_i}{G} \times 100$, với $g_i = \sum \frac{\pi \cdot D_{1.3}^2}{4}$).

  2. Nội suy sinh khối khô và quy đổi trữ lượng $\text{CO}_2$:

    • Khối lượng khô vật rơi rụng ($DW_{\text{ha}}$, tấn/ha): $$DW_{\text{ha}} = \left( \frac{1}{k} \sum_{j=1}^{k} dw_{\text{odb}, j} \right) \times 10.000 \times 10^{-3}$$
    • Lượng Carbon tích lũy ($C_{\text{stock}}$, tấn/ha): $$C_{\text{stock}} = DW_{\text{ha}} \times CF \quad (\text{với } CF = 0,5)$$
    • Lượng $\text{CO}2$ hấp thụ tương đương ($\text{CO}{2\text{-eq}}$, tấn/ha): $$\text{CO}{2\text{-eq}} = C{\text{stock}} \times \frac{\text{MW}{\text{CO}2}}{\text{AW}{\text{C}}} = C{\text{stock}} \times \frac{44}{12}$$

Mã nguồn Python minh họa thuật toán xử lý dữ liệu và kiểm chứng tính toán:

import pandas as pd
import numpy as np

def calculate_forest_carbon(plot_data: pd.DataFrame) -> pd.DataFrame:
    """
    Hàm tính toán sinh khối khô, trữ lượng Carbon và lượng CO2 tương đương
    dựa trên số liệu sấy mẫu vật rơi rụng (Litterfall).
    """
    IPCC_CARBON_FRACTION = 0.50
    CO2_C_RATIO = 44.0 / 12.0
    CARBON_PRICE_USD_PER_TON = 5.0
    
    results = []
    for idx, row in plot_data.iterrows():
        dry_twig_ton_ha = row['dry_twig_otc'] * 4.0
        dry_leaf_ton_ha = row['dry_leaf_otc'] * 4.0
        total_dry_biomass = dry_twig_ton_ha + dry_leaf_ton_ha
        
        c_twig = dry_twig_ton_ha * IPCC_CARBON_FRACTION
        c_leaf = dry_leaf_ton_ha * IPCC_CARBON_FRACTION
        total_c = total_dry_biomass * IPCC_CARBON_FRACTION
        
        co2_equivalent = total_c * CO2_C_RATIO
        commercial_value = co2_equivalent * CARBON_PRICE_USD_PER_TON
        
        results.append({
            'plot_id': row['plot_id'],
            'position': row['position'],
            'dry_biomass_ton_ha': round(total_dry_biomass, 3),
            'c_accumulated_ton_ha': round(total_c, 4),
            'co2_equivalent_ton_ha': round(co2_equivalent, 4),
            'value_usd_ha': round(commercial_value, 2)
        })
    return pd.DataFrame(results)

# Dữ liệu kiểm chứng thực tế từ 12 OTC tại Xã La Bằng
raw_sample_data = pd.DataFrame([
    {'plot_id': 1, 'position': 'Chân', 'dry_twig_otc': 1.902, 'dry_leaf_otc': 1.888},
    {'plot_id': 2, 'position': 'Sườn', 'dry_twig_otc': 0.847, 'dry_leaf_otc': 1.375},
    {'plot_id': 3, 'position': 'Đỉnh', 'dry_twig_otc': 1.132, 'dry_leaf_otc': 1.169},
    {'plot_id': 4, 'position': 'Chân', 'dry_twig_otc': 1.297, 'dry_leaf_otc': 1.436},
    {'plot_id': 5, 'position': 'Sườn', 'dry_twig_otc': 0.544, 'dry_leaf_otc': 0.558},
    {'plot_id': 6, 'position': 'Đỉnh', 'dry_twig_otc': 1.183, 'dry_leaf_otc': 0.397},
    {'plot_id': 7, 'position': 'Chân', 'dry_twig_otc': 0.431, 'dry_leaf_otc': 0.696},
    {'plot_id': 8, 'position': 'Sườn', 'dry_twig_otc': 1.072, 'dry_leaf_otc': 0.552},
    {'plot_id': 9, 'position': 'Đỉnh', 'dry_twig_otc': 1.503, 'dry_leaf_otc': 0.647},
    {'plot_id': 10, 'position': 'Chân', 'dry_twig_otc': 2.250, 'dry_leaf_otc': 2.323},
    {'plot_id': 11, 'position': 'Sườn', 'dry_twig_otc': 1.027, 'dry_leaf_otc': 1.043},
    {'plot_id': 12, 'position': 'Đỉnh', 'dry_twig_otc': 0.681, 'dry_leaf_otc': 0.634}
])

processed_results = calculate_forest_carbon(raw_sample_data)
print(f"Giá trị trung bình sinh khối khô: {processed_results['dry_biomass_ton_ha'].mean():.3f} tấn/ha")
print(f"Trữ lượng CO2 tương đương trung bình: {processed_results['co2_equivalent_ton_ha'].mean():.4f} tấn/ha")

Testing và validation

Quá trình kiểm định chất lượng số liệu được thực hiện nghiêm ngặt qua các bước:

  1. Kiểm tra độ ổn định sấy khô: Mẫu sấy đạt chuẩn khi độ biến thiên khối lượng qua các mốc 6h và 8h đạt $\Delta m \le 0,005\text{g}$, loại trừ hoàn toàn ảnh hưởng của độ ẩm dư thừa.
  2. Kiểm tra độ biến động theo địa hình:
    • Khu vực Chân núi (OTC 1, 4, 7, 10): Sinh khối khô trung bình đạt $12,222\text{ tấn/ha}$ (cao nhất do tích tụ cơ học từ sườn dốc và độ tàn che dày).
    • Khu vực Sườn núi (OTC 2, 5, 8, 11): Sinh khối khô trung bình đạt $7,018\text{ tấn/ha}$.
    • Khu vực Đỉnh núi (OTC 3, 6, 9, 12): Sinh khối khô trung bình đạt $7,346\text{ tấn/ha}$.

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

1. Đặc điểm cấu trúc lâm phần và tầng cây tái sinh rừng IIA

  • Mật độ tầng cây cao: Biến động từ $288$ đến $324\text{ cây/ha}$, trung bình toàn lâm phần đạt $316\text{ cây/ha}$. Số loài cây gỗ xuất hiện dao động từ $20$ đến $23\text{ loài/OTC}$ (trung bình $22\text{ loài/OTC}$).
  • Tổ thành loài ưu thế ($IVI \ge 5%$): Lâm phần không có loài độc tôn chiếm tuyệt đối mà đặc trưng bởi nhóm cây gỗ phục hồi nhanh: Kẹn (Castanopsis), Lim xẹt (Peltophorum dasyrhachis), Thành ngạnh (Cratoxylum formosum), Chẹo tía (Engelhardtia chrysolepis), Mán đỉa (Archidendron clypearia), Thẩu tấu (Aporosa dioica).
  • Mật độ tầng cây tái sinh: Dao động từ $4.000$ đến $5.520\text{ cây/ha}$, trung bình đạt $4.720\text{ cây/ha}$ với độ phong phú trung bình $19\text{ loài/OTC}$. Có 6 đến 10 loài tái sinh có triển vọng kinh tế và sinh thái cao (Thành ngạnh, Mỡ, Kháo, Kẹn, Lim xẹt, Chò chỉ), khẳng định tiềm năng diễn thế phục hồi mạnh mẽ thành rừng giàu trong tương lai.

2. Định lượng sinh khối vật rơi rụng

  • Sinh khối tươi: Trung bình đạt $15,500\text{ tấn/ha}$ (biến động từ $10,380$ đến $22,100\text{ tấn/ha}$). Trong đó:
    • Bộ phận cành rụng: $8,143\text{ tấn/ha}$ (chiếm $52,5%$).
    • Bộ phận lá, hoa, quả: $7,363\text{ tấn/ha}$ (chiếm $47,5%$).
  • Sinh khối khô: Trung bình đạt $8,863\text{ tấn/ha}$ (biến động từ $4,408$ đến $18,290\text{ tấn/ha}$). Trong đó:
    • Bộ phận cành rụng: $4,623\text{ tấn/ha}$ (chiếm $52,2%$).
    • Bộ phận lá, hoa, quả: $4,240\text{ tấn/ha}$ (chiếm $47,8%$).

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

  • Lượng Carbon ($C$) lưu trữ: Trung bình toàn lâm phần đạt $4,4315\text{ tấn C/ha}$ (Cành rụng tích lũy $2,3115\text{ tấn C/ha}$; Lá, hoa, quả tích lũy $2,1200\text{ tấn C/ha}$).
  • Lượng $\text{CO}_2$ hấp thụ tương đương: Đạt $16,2488\text{ tấn }\text{CO}_2\text{/ha}$ (Cành rụng đóng góp $8,4755\text{ tấn }\text{CO}_2\text{/ha}$; Lá, hoa, quả đóng góp $7,7733\text{ tấn }\text{CO}_2\text{/ha}$).

4. Giá trị thương mại Carbon ước tính

  • Với đơn giá tham chiếu chuẩn trong các chương trình lâm nghiệp tự nguyện là $5,00\text{ USD/tấn }\text{CO}_2$:
    • Giá trị kinh tế từ bể vật rơi rụng đạt $81,244\text{ USD/ha}$ (tương đương khoảng $1.750.000\text{ VNĐ/ha}$).
    • Tổng tiềm năng giảm phát thải của bể vật rơi rụng trên toàn bộ diện tích rừng sản xuất của xã La Bằng ($422\text{ ha}$) đạt $6.857\text{ tấn }\text{CO}_2$, tương đương giá trị tài sản carbon là $34.285\text{ USD}$.
    • Tính trên tổng quy mô $1.518\text{ ha}$ đất có rừng toàn xã, trữ lượng tích lũy đạt $24.665,68\text{ tấn }\text{CO}_2$, mang lại giá trị tiềm năng xấp xỉ $123.328\text{ USD}$.

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

Nghiên cứu mang lại nhiều đổi mới về phương pháp luận và cơ sở dữ liệu thực nghiệm:

  1. Phân tích tách biệt cấu phần cành và lá/hoa/quả: Thay vì đo gộp thảm mục như các nghiên cứu thông thường, việc phân chia cho thấy cành rơi rụng đóng góp tới $52,2%$ sinh khối khô thực tế. Cành cây có tỷ lệ hóa gỗ và lignin cao hơn, giúp carbon được cố định lâu hơn trong hệ sinh thái so với lá phân hủy nhanh.
  2. Khẳng định tầm quan trọng của bể vật rơi rụng trong hạch toán REDD+: Cung cấp bằng chứng thực nghiệm rằng tầng thảm mục lưu giữ tới $16,25\text{ tấn }\text{CO}_2\text{/ha}$, chứng minh việc bỏ qua bể chứa này sẽ làm suy giảm nghiêm trọng giá trị thẩm định tín chỉ carbon của rừng.
So sánh sinh khối thảm mục giữa các nghiên cứu quốc tế và Việt Nam:
+------------------------------------------------------+---------------------+
| Nghiên cứu / Hệ sinh thái                            | Năng suất rơi rụng  |
+------------------------------------------------------+---------------------+
| La Bằng, Thái Nguyên (Rừng phục hồi IIA - Đề tài)   | 8,863 tấn khô/ha    |
| Daniel & Adam (1984, New Zealand - Rừng lá rộng)     | 6,180 kg/ha/năm     |
| Levett et al. (1985, Rừng Sồi Nothofagus truncata)   | 7,252 kg/ha/năm     |
| Xiaoniu Xu et al. (2004, Okinawa - Rừng á nhiệt đới) | 7,558 kg/ha/năm     |
| Weaver et al. (1986, Puerto Rico - Rừng lùn núi cao) | 3,100 kg/ha/năm     |
+------------------------------------------------------+---------------------+

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

Kịch bản ứng dụng thực tế

  • Xây dựng Hồ sơ Thiết kế Dự án (PDD - Project Design Document) theo tiêu chuẩn VCS/CCB: Tích hợp số liệu thảm mục vào đường carbon cơ sở cho các dự án REDD+ và phục hồi rừng bền vững tại khu vực miền núi phía Bắc.
  • Chi trả Dịch vụ Môi trường Rừng (PFES / DVMTR): Cung cấp cơ sở định lượng khoa học để các nhà máy thủy điện, khu công nghiệp và doanh nghiệp phát thải lớn chi trả tiền hấp thụ và lưu giữ carbon rừng cho các cộng đồng địa phương tại Đại Từ, Thái Nguyên.
Lộ trình triển khai hệ thống MRV (Đo đạc - Báo cáo - Thẩm tra) Carbon rừng:
[Giai đoạn 1: Q1] Thiết lập mạng lưới OTC chuẩn & Thu mẫu vật rơi rụng 
        ↓
[Giai đoạn 2: Q2] Phân tích giải tích sấy khô tại Lab & Tính toán Carbon Pool
        ↓
[Giai đoạn 3: Q3] Tích hợp số liệu vào GIS/Web-GIS Quản lý Rừng cộng đồng
        ↓
[Giai đoạn 4: Q4] Thẩm định độc lập cấp chứng chỉ & Thương mại hóa tín chỉ

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

Hạn chế kỹ thuật

  • Thời gian nghiên cứu tập trung từ tháng 01 đến tháng 05/2015, chưa bao quát trọn vẹn chu kỳ 12 tháng để ghi nhận đỉnh rụng lá vào mùa thu - đông khi có bão hoặc sương muối kéo dài.
  • Chưa đo lường trực tiếp tốc độ phân hủy ($k$-constant) của thảm mục trong từng tiểu vùng vi khí hậu.

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

  • Lắp đặt bẫy vật rơi rụng cố định (litter traps) thu mẫu liên tục 2 tuần/lần trong 3-5 năm để lập mô hình chu chuyển vật chất và dinh dưỡng khoáng (N, P, K, Ca, Mg) trả lại cho đất.
  • Ứng dụng dữ liệu viễn thám vệ tinh đa phổ (Sentinel-2, Landsat-8) kết hợp thuật toán Random Forest/Machine Learning để ước lượng sinh khối thảm mục trên toàn bộ lưu vực dãy núi Tam Đảo.

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

  • Cộng đồng dân cư địa phương & Chủ rừng (Xã La Bằng): Tạo cơ sở pháp lý và kinh tế để nhận tiền chi trả dịch vụ môi trường rừng với mức đóng góp tiềm năng $81,24\text{ USD/ha}$, nâng cao sinh kế gắn liền với bảo vệ rừng.
  • Cơ quan Quản lý Lâm nghiệp & Ban Quản lý Vườn Quốc gia Tam Đảo: Sở hữu bộ dữ liệu thực chứng phục vụ công tác theo dõi diễn biến rừng, phòng cháy chữa cháy rừng (đánh giá vật liệu cháy từ thảm mục) và quy hoạch bảo tồn đa dạng sinh học.
  • Các nhà phát triển dự án Carbon & Doanh nghiệp: Có nguồn số liệu tin cậy để hoàn thiện hồ sơ đăng ký tín chỉ carbon quốc tế, rút ngắn thời gian thẩm định PDD từ 6-8 tháng.
  • Sinh viên & Giới nghiên cứu Lâm sinh: Cung cấp quy trình thực nghiệm chuẩn hóa về thu mẫu, xử lý giải tích và công thức quy đổi carbon áp dụng cho rừng nhiệt đới hỗn loài.

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

1. Yêu cầu kỹ thuật cốt lõi để lập ô dạng bản (ÔDB) thu thập vật rơi rụng là gì?

Trên mỗi ô tiêu chuẩn ($50\text{m} \times 50\text{m}$), bố trí 5 ô dạng bản kích thước $1\text{m} \times 1\text{m}$ ($1\text{m}^2$) tại 4 góc và tâm ô theo hai đường chéo. Toàn bộ cành có đường kính $D < 5\text{ cm}$, chiều dài $L < 50\text{ cm}$, lá, hoa, quả rụng phải được thu gom sạch sẽ, cạo sạch đất đá dính kèm trước khi cân khối lượng tươi.

2. Vì sao áp dụng hệ số chuyển đổi Carbon $CF = 0,5$ và tỷ lệ $44/12$?

Hệ số $CF = 0,5$ là thông số mặc định được Ủy ban Liên chính phủ về Biến đổi Khí hậu (IPCC) công nhận toàn cầu đối với sinh khối khô thực vật nhiệt đới (hàm lượng carbon hữu cơ chiếm xấp xỉ $50%$ khối lượng khô cấu tạo từ cellulose, hemicellulose và lignin). Tỷ lệ $44/12$ là tỷ số khối lượng phân tử giữa khí $\text{CO}_2$ ($\text{MW} = 44\text{ g/mol}$) và nguyên tử Carbon ($\text{AW} = 12\text{ g/mol}$).

3. Dữ liệu này có thể tích hợp trực tiếp vào phần mềm CO2FIX V3.1 không?

Hoàn toàn khả thi. Dữ liệu sinh khối khô cành và lá ($4,623\text{ tấn/ha}$ và $4,240\text{ tấn/ha}$) cùng thông số lâm phần ($316\text{ cây/ha}$) là các biến số đầu vào chuẩn (Input Parameters) cho mô-đun Litter & Soil của phần mềm CO2FIX để mô phỏng dòng luân chuyển carbon qua các năm.

4. Vị trí địa hình ảnh hưởng như thế nào đến sự phân bố sinh khối thảm mục?

Sinh khối khô vật rơi rụng tập trung cao nhất ở khu vực chân núi ($12,222\text{ tấn/ha}$), cao hơn đáng kể so với đỉnh núi ($7,346\text{ tấn/ha}$) và sườn núi ($7,018\text{ tấn/ha}$). Nguyên nhân do gió, trọng lực và dòng chảy bề mặt cuốn trôi vật liệu nhẹ từ sườn dốc xuống chân núi, kết hợp với mật độ che phủ tầng cây cao tại chân núi dày hơn.

5. Chi phí đo kiểm và thời gian hoàn vốn đầu tư dự án tín chỉ Carbon từ rừng IIA?

Chi phí đo đạc thực địa và giải tích phòng thí nghiệm ước tính khoảng $3 - 5\text{ USD/ha}$. Với giá trị thương mại đạt $81,24\text{ USD/ha}$ chỉ tính riêng cho tầng vật rơi rụng (chưa tính tầng cây đứng), dự án có tỷ suất sinh lời cao và thời gian hoàn vốn ngay sau chu kỳ thẩm định tín chỉ đầu tiên (khoảng 12 - 18 tháng).

Kết luận

Nghiên cứu đã xác lập thành công cơ sở dữ liệu định lượng chính xác về cấu trúc lâm phần và năng lực hấp thụ $\text{CO}_2$ của tầng vật rơi rụng dưới tán rừng phục hồi trạng thái IIA tại xã La Bằng, huyện Đại Từ, tỉnh Thái Nguyên. Với mật độ tầng cây cao $316\text{ cây/ha}$ và mật độ tái sinh mạnh mẽ $4.720\text{ cây/ha}$, hệ sinh thái rừng IIA tại khu vực thể hiện khả năng phục hồi sinh thái vượt trội. Tầng vật rơi rụng tích lũy trung bình $8,863\text{ tấn/ha}$ sinh khối khô, lưu giữ $4,4315\text{ tấn C/ha}$ và hấp thụ tương đương $16,2488\text{ tấn }\text{CO}_2\text{/ha}$, mang lại giá trị kinh tế ước tính $81,244\text{ USD/ha}$.

Kết quả này khẳng định tầng thảm mục là một bể chứa carbon quan trọng không thể tách rời trong các chương trình chi trả dịch vụ môi trường rừng, dự án REDD+ và chiến lược phát triển thị trường tín chỉ carbon quốc gia. Các tổ chức lâm nghiệp, doanh nghiệp và nhà nghiên cứu quan tâm có thể ứng dụng trực tiếp quy trình thực nghiệm và bộ tham số này để chuẩn hóa công tác kiểm kê khí nhà kính và thúc đẩy phát triển lâm nghiệp bền vững thích ứng với biến đổi khí hậu toàn cầu.