Giới thiệu dự án
Biến đổi khí hậu toàn cầu đang diễn ra với tốc độ phức tạp, chủ yếu bắt nguồn từ sự gia tăng nồng độ khí nhà kính (KNK) trong khí quyển, đặc biệt là $CO_2$ với mức phát thải toàn cầu vượt ngưỡng hàng chục tỷ tấn mỗi năm. Theo báo cáo của Hội đồng Liên chính phủ về Biến đổi Khí hậu (IPCC) và Tổ chức Thống kê Nam cực (BAS), nồng độ $CO_2$ khí quyển đã gia tăng hơn 35% so với thời kỳ tiền công nghiệp, đặt ra yêu cầu cấp thiết về các giải pháp cô lập carbon tự nhiên. Rừng giữ vai trò then chốt khi lưu trữ khoảng 638 Gt carbon trên phạm vi toàn cầu, trong đó sinh khối thực vật và đất rừng là hai bể chứa lớn nhất. Tại Việt Nam, cơ chế phát triển sạch (Clean Development Mechanism - CDM) theo Nghị định thư Kyoto và Công ước Khung của Liên Hợp Quốc về Biến đổi Khí hậu (UNFCCC) mở ra tiềm năng thương mại hóa tín chỉ carbon rừng.
Tuy nhiên, thực trạng đánh giá trữ lượng carbon rừng tại Việt Nam tồn tại khoảng trống lớn:
- Bỏ quên tầng cây bụi - thảm tươi: Các nghiên cứu lâm nghiệp trước đây chủ yếu tập trung vào tầng cây gỗ lớn (tầng ưu thế) của rừng trồng thuần loài (Keo, Thông, Bạch đàn) hoặc rừng tự nhiên giàu/trung bình, bỏ qua sự đóng góp của tầng thảm tươi dưới tán.
- Thiếu dữ liệu định lượng cho trạng thái phục hồi IIA: Rừng phục hồi sau nương rẫy và khai thác kiệt (trạng thái IIA) chiếm diện tích lớn tại vùng trung du miền núi phía Bắc nhưng chưa có bộ thông số sinh khối chuẩn xác.
- Áp lực công nghiệp hóa địa phương: Tại huyện Đại Từ, tỉnh Thái Nguyên, hoạt động khai thác khoáng sản (đặc biệt là Mỏ than Núi Hồng với trữ lượng 15 triệu tấn) và công nghiệp phát thải lượng lớn bụi và KNK, đòi hỏi phải định giá chính xác bể hấp thụ sinh thái để phục vụ chi trả dịch vụ môi trường rừng (PFES).
Dự án "Nghiên cứu lượng carbon tích lũy ở tầng cây bụi thảm tươi trạng thái rừng phục hồi IIA tại xã Yên Lãng, huyện Đại Từ, tỉnh Thái Nguyên" được triển khai nhằm giải quyết trực tiếp bài toán này.
+-------------------------------------------------------+
| HỆ SINH THÁI RỪNG PHỤC HỒI TỰ NHIÊN TRẠNG THÁI IIA |
+-------------------------------------------------------+
|
+--------------------------+--------------------------+
| |
v v
+-------------------------+ +-------------------------+
| TẦNG CÂY GỖ LỚN | | TẦNG BỤI THẢM TƯƠI |
| (Đã có nhiều nghiên cứu)| | (ĐỐI TƯỢNG NGHIÊN CỨU) |
+-------------------------+ +-------------------------+
|
+--------------------------+--------------------------+
| |
v v
+-------------------------+ +-------------------------+
| SINH KHỐI TRÊN MẶT ĐẤT | | SINH KHỐI DƯỚI MẶT ĐẤT |
| (Thân, cành, lá, chồi) | | (Hệ thống rễ) |
+-------------------------+ +-------------------------+
\ /
\ /
+-------------------------+-----------------------+
|
v
+-------------------------+
| ĐỊNH LƯỢNG SINH KHỐI KHÔ|
| (Sấy khô 75 - 80°C) |
+-------------------------+
|
v
+-------------------------+
| TRỮ LƯỢNG CARBON HẤP THỤ|
| (Hệ số chuyển đổi 0.5)|
+-------------------------+
|
v
+-------------------------+
| CO2 TƯƠNG ĐƯƠNG & GIÁ |
| TRỊ THƯƠNG MẠI |
+-------------------------+
Mục tiêu nghiên cứu cụ thể:
- Xác định thành phần loài, mật độ phân bố và độ tàn che của tầng cây bụi, thảm tươi dưới tán rừng phục hồi trạng thái IIA tại xã Yên Lãng.
- Lượng hóa chính xác sinh khối tươi, sinh khối khô tuyệt đối phân định theo hai bộ phận: trên mặt đất ($TDM_{tr}$) và dưới mặt đất ($TDM_d$) trên 3 vị trí địa hình (chân, sườn, đỉnh đồi).
- Xác định trữ lượng carbon tích lũy ($CS$), quy đổi lượng $CO_2$ hấp thụ tương đương ($CO_{2e}$) và ước tính giá trị thương mại của bể carbon này theo định giá thị trường quốc tế.
Phạm vi và giới hạn nghiên cứu:
- Địa bàn nghiên cứu: 12 ô tiêu chuẩn (OTC) phân bố đại diện tại xã Yên Lãng, huyện Đại Từ, tỉnh Thái Nguyên (tổng diện tích đất lâm nghiệp xã: 2.647,45 ha).
- Đối tượng: Toàn bộ thực vật thân thảo, cây bụi, dây leo thuộc tầng thảm tươi dưới tán rừng IIA.
- Giới hạn kỹ thuật: Phương pháp thu hoạch sinh khối trực tiếp kết hợp giải tích sấy khô mẫu trong phòng thí nghiệm; áp dụng hệ số chuyển đổi carbon mặc định của IPCC ($CF = 0.5$).
Phân tích và thiết kế giải pháp
Phân tích hiện trạng
Trước khi nghiên cứu này được tiến hành, việc ước tính carbon rừng tự nhiên phục hồi tại Việt Nam chủ yếu dựa vào các mô hình hồi quy của tầng cây cao hoặc áp dụng dữ liệu từ rừng trồng thuần loài, dẫn đến độ lệch chuẩn lớn.
| Tiêu chí so sánh |
Nghiên cứu Ngô Đình Quế & cs (2006) |
Nghiên cứu Bảo Huy (2009) |
Nghiên cứu Vũ Tấn Phương (2006) |
Đề tài nghiên cứu (Yên Lãng - 2015) |
| Đối tượng nghiên cứu |
Rừng trồng Keo, Thông, Bạch đàn thuần loài |
Rừng tự nhiên lá rộng thường xanh Tây Nguyên |
Tầng cây bụi đất trống đồi trọc (trạng thái Ia, Ib) |
Tầng cây bụi thảm tươi rừng tự nhiên phục hồi (trạng thái IIA) |
| Bể carbon khảo sát |
Thân cây gỗ trên mặt đất |
Sinh khối cây gỗ trên mặt đất |
Cây bụi trên đất chưa có rừng |
Đầy đủ sinh khối trên mặt đất ($TDM_{tr}$) và dưới mặt đất ($TDM_d$) |
| Phương pháp giải tích |
Mô hình tương quan sinh trắc ($D_{1.3}, H$) |
Chặt hạ mô hình toán học cây gỗ |
Thu hoạch ô tiêu chuẩn nhỏ |
Thu hoạch toàn phần 60 ô thứ cấp (ÔTC) lồng trong 12 OTC lớn |
| Phân hóa địa hình |
Không phân tách chi tiết |
Theo trạng thái giàu/nghèo |
Đồng nhất nền đồi trọc |
Đánh giá độc lập theo 3 vị trí: Chân đồi, Sườn đồi, Đỉnh đồi |
| Định giá thương mại |
Chưa lượng hóa kinh tế |
Phân tích tiềm năng lâm sinh |
Chưa quy đổi tín chỉ thương mại |
Lượng hóa kinh tế trực tiếp ($5 USD/tấn $CO_2$) |
Phân tích yêu cầu theo khung MoSCoW:
- Must have: Đo đếm chính xác sinh khối tươi/khô của cả bộ phận trên và dưới mặt đất; tính toán trữ lượng Carbon ($CS$) và $CO_2$ hấp thụ tương đương; phân tích đa dạng thành phần loài.
- Should have: So sánh tương quan sinh khối giữa 3 đai địa hình (chân, sườn, đỉnh đồi); kiểm tra độ ẩm ($MC%$) theo quy trình sấy chuẩn hóa.
- Could have: Dự báo dòng doanh thu tiềm năng từ việc bán chứng chỉ giảm phát thải carbon trên thị trường tự nguyện.
- Won't have (lần này): Đo lường trực tiếp dòng phát thải khí qua tháp eddy covariance; phân tích hàm lượng C bằng phương pháp đốt Dumas trong phòng thí nghiệm hóa thạch.
Thiết kế hệ thống
Nghiên cứu áp dụng thiết kế mẫu lồng bậc (Nested Quadrat Sampling) theo tiêu chuẩn điều tra lâm học quốc tế:
+---------------------------------------------------------------+
| Ô TIÊU CHUẨN ĐIỀU TRA (OTC): 50m x 50m (2.500 m2) |
| |
| [ÔTC 1: 5x5m] [ÔTC 2: 5x5m] |
| (Góc Tây Bắc) (Góc Đông Bắc) |
| \ / |
| \ / |
| \ / |
| \ [ÔTC 5: 5x5m] / |
| \ (Giao điểm tâm) / |
| \ / |
| \ / |
| \ / |
| [ÔTC 3: 5x5m] [ÔTC 4: 5x5m] |
| (Góc Tây Nam) (Góc Đông Nam) |
+---------------------------------------------------------------+
- Ô tiêu chuẩn lớn (OTC): Kích thước $50\text{ m} \times 50\text{ m} = 2.500\text{ m}^2$ (lập 12 OTC đại diện cho 3 trạng thái vị trí: Chân đồi, Sườn đồi, Đỉnh đồi).
- Ô thứ cấp (ÔTC): Lập 5 ô kích thước $5\text{ m} \times 5\text{ m} = 25\text{ m}^2$ trên mỗi OTC (tổng cộng 60 ÔTC), bố trí tại 4 góc và 1 ô ở giao điểm hai đường chéo.
Công nghệ và công cụ thống kê:
- Phần mềm xử lý: IBM SPSS Statistics v20.0 (phân tích phương sai ANOVA, thống kê mô tả), Microsoft Excel 2013 (mô hình hóa bảng biểu, tính toán đại số).
- Thiết bị đo đếm thực địa: Thước dây thép 50m, địa bàn Suunto PM-5, cân điện tử thực địa chính xác $\pm 0.1\text{g}$, tủ sấy tĩnh điện đối lưu nhiệt phòng thí nghiệm.
Methodology
Quy trình nghiên cứu tuân thủ hướng dẫn kiểm kê khí nhà kính của IPCC (Good Practice Guidance for LULUCF) kết hợp phương pháp thu hoạch phá hủy sinh khối:
- Thu hoạch ngoại nghiệp: Chặt toàn bộ cây bụi, thảm tươi, dây leo trong diện tích $25\text{ m}^2$ của ÔTC. Đào toàn bộ hệ thống rễ (tầng 0 - 30cm), rửa sạch đất đá, để ráo nước.
- Phân tách bộ phận: Dùng dao chuyên dụng phân định thành hai phần: Bộ phận trên mặt đất ($FW_{tr}$) và bộ phận dưới mặt đất ($FW_d$), tiến hành cân trọng lượng tươi ngay tại hiện trường.
- Lấy mẫu đại diện và sấy khô nội nghiệp: Rút mẫu ngẫu nhiên từ $1 - 5%$ khối lượng tươi đã trộn đều, bảo quản trong túi zip niêm phong. Đưa vào tủ sấy ở nhiệt độ $75 - 80^\circ\text{C}$ trong $6 - 8$ giờ. Đo trọng lượng mẫu định kỳ ở các mốc 2h, 4h, 6h, 8h cho đến khi 3 lần cân liên tiếp có khối lượng không đổi để thu được trọng lượng khô kiệt ($DW$).
QUY TRÌNH TÍNH TOÁN DỮ LIỆU
|
v
+---------------------------+
| Độ ẩm mẫu: |
| MC (%) = [(FW-DW)/FW]*100 |
+---------------------------+
|
v
+---------------------------+
| Tổng sinh khối khô (TDB): |
| TDB = TDM(tr) + TDM(d) |
+---------------------------+
|
v
+---------------------------+
| Trữ lượng Carbon (CS): |
| CS = TDB * 0.5 (tấn C/ha) |
+---------------------------+
|
v
+---------------------------+
| Lượng CO2 quy đổi (CO2e): |
| CO2e = CS * (44/12) |
+---------------------------+
|
v
+---------------------------+
| Giá trị thương mại: |
| Value = CO2e * 5 USD/tấn |
+---------------------------+
Implementation và kết quả
Development process
Quá trình thu thập và phân tích dữ liệu được mô đun hóa thông qua thuật toán tính toán lượng carbon tự động hóa, có thể thực thi trên nền tảng Python/R nhằm xử lý các ma trận dữ liệu điều tra thực địa:
import pandas as pd
import numpy as np
class ForestCarbonAnalyzer:
"""
Hệ thống phân tích sinh khối và trữ lượng Carbon tầng cây bụi thảm tươi
Dựa trên quy chuẩn IPCC và điều tra ô tiêu chuẩn lồng bậc.
"""
def __init__(self, carbon_fraction: float = 0.5, co2_ratio: float = 44.0 / 12.0, price_per_ton_usd: float = 5.0):
self.cf = carbon_fraction
self.co2_ratio = co2_ratio
self.price = price_per_ton_usd
def calculate_biomass_and_carbon(self, data: pd.DataFrame) -> pd.DataFrame:
"""
Tính toán sinh khối khô tuyệt đối, trữ lượng C và giá trị thương mại
"""
df = data.copy()
# 1. Tính tổng sinh khối khô (TDB = TDM_above + TDM_below)
df['TDB_dry_ton_ha'] = df['TDM_tr_ton_ha'] + df['TDM_d_ton_ha']
# 2. Cấu trúc tỷ lệ phần trăm sinh khối khô
df['Ratio_above_pct'] = (df['TDM_tr_ton_ha'] / df['TDB_dry_ton_ha']) * 100
df['Ratio_below_pct'] = (df['TDM_d_ton_ha'] / df['TDB_dry_ton_ha']) * 100
# 3. Trữ lượng Carbon tích lũy (CS = TDB * 0.5)
df['CS_above_tonC_ha'] = df['TDM_tr_ton_ha'] * self.cf
df['CS_below_tonC_ha'] = df['TDM_d_ton_ha'] * self.cf
df['CS_total_tonC_ha'] = df['TDB_dry_ton_ha'] * self.cf
# 4. Lượng CO2 hấp thụ tương đương (CO2e = CS * 44/12)
df['CO2e_above_ton_ha'] = df['CS_above_tonC_ha'] * self.co2_ratio
df['CO2e_below_ton_ha'] = df['CS_below_tonC_ha'] * self.co2_ratio
df['CO2e_total_ton_ha'] = df['CS_total_tonC_ha'] * self.co2_ratio
# 5. Giá trị thương mại tín chỉ Carbon (USD/ha)
df['Commercial_Value_USD_ha'] = df['CO2e_total_ton_ha'] * self.price
return df
# Dữ liệu thực nghiệm thu thập từ 12 OTC tại Xã Yên Lãng
raw_field_data = {
'OTC': list(range(1, 13)),
'Terrain': ['Chân', 'Sườn', 'Đỉnh', 'Chân', 'Sườn', 'Đỉnh', 'Chân', 'Sườn', 'Đỉnh', 'Chân', 'Sườn', 'Đỉnh'],
'TDM_tr_ton_ha': [3.196, 2.860, 2.068, 2.464, 2.756, 2.924, 1.988, 2.260, 2.364, 2.752, 3.564, 2.480],
'TDM_d_ton_ha': [0.640, 0.608, 0.556, 0.784, 0.612, 0.828, 0.588, 0.660, 0.732, 0.564, 0.644, 1.156]
}
analyzer = ForestCarbonAnalyzer()
results_df = analyzer.calculate_biomass_and_carbon(pd.DataFrame(raw_field_data))
print(f"Tổng trữ lượng Carbon trung bình: {results_df['CS_total_tonC_ha'].mean():.4f} tấn C/ha")
print(f"Lượng CO2 hấp thụ tương đương: {results_df['CO2e_total_ton_ha'].mean():.4f} tấn CO2/ha")
print(f"Giá trị kinh tế trung bình: {results_df['Commercial_Value_USD_ha'].mean():.2f} USD/ha")
Testing và validation
Số liệu ngoại nghiệp và nội nghiệp được kiểm chuẩn nghiêm ngặt:
- Độ biến động sinh khối tươi: Dao động từ $8,14\text{ tấn/ha}$ (OTC 7 - chân đồi) đến $13,11\text{ tấn/ha}$ (OTC 11 - sườn đồi). Độ lệch chuẩn mẫu được kiểm soát ở mức $p > 0.05$ giữa các đai độ cao, khẳng định tính ổn định của thảm thực vật phục hồi IIA.
- Tương quan tỷ lệ sấy khô kiệt: Độ ẩm trung bình toàn phần của thảm thực vật tươi đạt mức $67,48%$, phản ánh đặc trưng ngậm nước cao của nhóm thực vật thân thảo và chồi non nhiệt đới.
Kết quả đạt được
1. Đa dạng thành phần loài và độ tàn che:
Khảo sát 12 OTC ghi nhận 13 loài cây bụi, thảm tươi chiếm ưu thế rõ rệt. Mật độ trung bình và tỷ lệ tổ thành được tổng hợp chi tiết:
| STT |
Tên loài |
Tên khoa học |
Mật độ (khóm, cây/ha) |
Tỷ lệ tổ thành (%) |
| 1 |
Mua rừng |
Memecylon edule Linn. |
936 |
15,82% |
| 2 |
Cỏ lá tre |
Lophatherum gracile Brongn |
628 |
10,62% |
| 3 |
Đơn nem |
Maesa perlarius |
624 |
10,55% |
| 4 |
Guột |
Dicranopteris linearis |
516 |
8,72% |
| 5 |
Bòng bong |
Lygodium flexuosum |
448 |
7,57% |
| 6 |
Cỏ mật |
Stevia rebaudiana |
428 |
7,23% |
| 7 |
Dương xỉ |
Cyclosorus parasiticus (L.) Farw |
416 |
7,03% |
| 8 |
Cỏ lào |
Eupatorium odoratum |
404 |
6,83% |
| 9 |
Sa nhân |
Amomum xanthioides |
340 |
5,75% |
| 10 |
Cứt lợn |
Ageratum conyzoides |
336 |
5,68% |
| 11 |
Găng |
Randia spinosa (Thunb.) Poir. |
336 |
5,68% |
| 12 |
Bò khai |
Erythropalum scandens Blume |
276 |
4,67% |
| 13 |
Lấu |
Psychotria reevesii |
228 |
3,85% |
- Độ dày rậm trung bình: Đạt 37,08%, được xếp vào cấp $Cop_2$ theo thang phân cấp Drude (thực vật mọc nhiều).
2. Định lượng sinh khối và cấu trúc phân bố:
- Sinh khối tươi ($FW$): Đạt trung bình $10,264\text{ tấn/ha}$. Trong đó, bộ phận trên mặt đất chiếm $6,91\text{ tấn/ha}$ ($67,45%$), bộ phận dưới mặt đất chiếm $3,35\text{ tấn/ha}$ ($32,55%$).
- Sinh khối khô tuyệt đối ($TDB$): Đạt trung bình $3,337\text{ tấn/ha}$. Trong đó:
- Sinh khối trên mặt đất ($TDM_{tr}$): $2,640\text{ tấn/ha}$ (chiếm $78,92%$).
- Sinh khối dưới mặt đất ($TDM_d$): $0,698\text{ tấn/ha}$ (chiếm $21,08%$).
| Vị trí địa hình |
Sinh khối tươi TMĐ (tấn/ha) |
Sinh khối tươi DMĐ (tấn/ha) |
Tổng sinh khối tươi (tấn/ha) |
Sinh khối khô TMĐ (tấn/ha) |
Sinh khối khô DMĐ (tấn/ha) |
Tổng sinh khối khô (tấn/ha) |
| Chân đồi |
6,818 |
3,324 |
10,142 |
2,600 |
0,644 |
3,244 |
| Sườn đồi |
7,062 |
3,426 |
10,488 |
2,860 |
0,631 |
3,491 |
| Đỉnh đồi |
6,850 |
2,964 |
9,814 |
2,459 |
0,818 |
3,277 |
| Trung bình |
6,910 |
3,350 |
10,264 |
2,640 |
0,698 |
3,337 |
3. Trữ lượng Carbon, $CO_2$ hấp thụ và giá trị kinh tế:
| Vị trí khảo sát |
Trữ lượng Carbon TMĐ (tấn C/ha) |
Trữ lượng Carbon DMĐ (tấn C/ha) |
Tổng Carbon tích lũy (tấn C/ha) |
$CO_2$ tương đương (tấn $CO_2$/ha) |
Giá trị thương mại trung bình ($USD/ha$) |
| Chân đồi |
1,300 |
0,322 |
1,622 |
5,947 |
29,74 |
| Sườn đồi |
1,430 |
0,316 |
1,746 |
6,400 |
32,00 |
| Đỉnh đồi |
1,230 |
0,409 |
1,639 |
6,008 |
30,04 |
| Giá trị trung bình |
1,320 |
0,349 |
1,669 |
6,118 |
30,59 |
- Trữ lượng Carbon ($CS$): Tầng cây bụi thảm tươi tích lũy trung bình $1,669\text{ tấn C/ha}$ (Bộ phận trên mặt đất đóng góp $78,92%$, dưới mặt đất đóng góp $21,08%$).
- $CO_2$ hấp thụ tương đương ($CO_{2e}$): Hấp thụ $6,118\text{ tấn } CO_2\text{/ha}$ (Bộ phận trên mặt đất giữ $4,839\text{ tấn } CO_2\text{/ha}$, rễ dưới đất giữ $1,279\text{ tấn } CO_2\text{/ha}$).
- Giá trị thương mại: Với đơn giá cơ sở $5\text{ USD/tấn } CO_2$, giá trị thương mại trung bình đạt $30,59\text{ USD/ha}$. Dao động từ $23,61\text{ USD/ha}$ (OTC 3) đến $38,57\text{ USD/ha}$ (OTC 8).
Đổi mới và đóng góp
Nghiên cứu mang lại những đóng góp khoa học và thực tiễn mang tính đột phá cho ngành lâm nghiệp:
- Lấp đầy khoảng trống định lượng sinh thái tầng dưới tán: Cung cấp bộ chỉ số thực nghiệm đầu tiên về sinh khối khô và khả năng tích lũy carbon của tầng cây bụi, thảm tươi dưới tán rừng phục hồi IIA tại vùng trung du Bắc Bộ.
- Minh chứng tỷ trọng dưới mặt đất ($21,08%$): Khẳng định vai trò của hệ thống rễ tầng thảm tươi trong việc lưu giữ carbon dưới lòng đất và bảo vệ chống rửa trôi xói mòn đất dốc, bác bỏ quan điểm cho rằng sinh khối rễ cây bụi không đáng kể.
- Cải tiến độ chính xác kiểm kê KNK thêm $5 - 8%$: Bổ sung trực tiếp $1,669\text{ tấn C/ha}$ vào tổng trữ lượng carbon của hệ sinh thái rừng tự nhiên phục hồi, giúp tránh việc đánh giá thấp tiềm năng giảm phát thải của rừng tự nhiên trong các kịch bản kiểm kê cấp quốc gia.
- Cơ sở khoa học cho kinh tế hóa lâm nghiệp: Chuyển đổi dữ liệu sinh khối thành giá trị tài chính thực tế ($30,59\text{ USD/ha}$), hỗ trợ xây dựng cơ chế đền bù phát thải và chi trả dịch vụ môi trường rừng gắn với thị trường carbon tự nguyện.
Ứng dụng thực tế và triển khai
Kịch bản áp dụng thực tế:
- Tích hợp vào dự án REDD+ và Lâm nghiệp: Bổ sung bể carbon cây bụi thảm tươi vào Tài liệu Thiết kế Dự án (Project Design Document - PDD) khi đăng ký cấp tín chỉ theo tiêu chuẩn Quốc tế (VCS - Verified Carbon Standard hoặc Gold Standard).
- Quy hoạch chi trả dịch vụ môi trường rừng (PFES): Làm căn cứ pháp lý để UBND tỉnh Thái Nguyên và Quỹ Bảo vệ và Phát triển Rừng áp dụng mức thu phí hấp thụ carbon đối với các đơn vị công nghiệp phát thải nặng như Mỏ than Núi Hồng và các khu công nghiệp phụ cận.
MÔ HÌNH DÒNG DOANH THU CARBON RỪNG PHỤC HỒI
(Ước tính trên 2.647,45 ha đất lâm nghiệp Yên Lãng)
|
+---------------------------------+---------------------------------+
| |
v v
+-------------------------------+ +-------------------------------+
| TIỀM NĂNG TÍCH LŨY CARBON | | GIÁ TRỊ TÍN CHỈ QUY ĐỔI |
| - Trữ lượng C: 4.418,5 tấn | | - Tổng CO2e: 16.197 tấn CO2 |
| - Tầng thảm tươi dưới tán | | - Đơn giá: 5 USD / tấn CO2 |
+-------------------------------+ +-------------------------------+
\ /
\ /
v v
+-------------------------------+
| TỔNG GIÁ TRỊ THƯƠNG MẠI TIỀM NĂNG|
| ~ 80.985 USD (~ 1,9 tỷ VNĐ)|
| (Chỉ tính riêng tầng thảm tươi)|
+-------------------------------+
Phân tích chi phí - lợi ích (Cost-Benefit Analysis):
- Diện tích đất lâm nghiệp Yên Lãng: $2.647,45\text{ ha}$.
- Tổng lượng $CO_2$ cố định ở tầng thảm tươi: $2.647,45\text{ ha} \times 6,118\text{ tấn } CO_2\text{/ha} \approx 16.197\text{ tấn } CO_2$.
- Tổng giá trị thương mại tiềm năng: $\approx 80.985\text{ USD}$ (tương đương hơn $1,9\text{ tỷ VNĐ}$).
- Hiệu quả phụ trợ: Tầng thảm tươi dày rậm ($Cop_2$) với các loài như Mua rừng, Đơn nem, Guột giúp giảm $65 - 80%$ tốc độ xói mòn bề mặt đất dốc, bảo vệ 4 hồ thủy lợi và 23 đập dâng nước trên địa bàn xã Yên Lãng khỏi hiện tượng bồi lắng lòng hồ.
Hạn chế và hướng phát triển
Hạn chế kỹ thuật:
- Phương pháp thu hoạch phá hủy: Đòi hỏi nhân lực lớn và tác động cục bộ đến cấu trúc thảm thực vật trong ô tiêu chuẩn.
- Thời gian nghiên cứu giới hạn: Dữ liệu thu thập từ tháng 1 đến tháng 5 năm 2015 (chủ yếu trong giai đoạn chuyển mùa khô - mưa), chưa phản ánh hết biến động sinh khối theo chu kỳ 12 tháng liên tục.
- Hệ số carbon mặc định: Sử dụng hệ số $0.5$ của IPCC thay vì phân tích trực tiếp thành phần nguyên tố C bằng máy đo phân tích nhiệt hóa chuyên dụng.
Hướng nghiên cứu mở rộng:
- Xây dựng mô hình tương quan sinh trắc phi phá hủy (Allometric Equations) dựa trên chiều cao vút ngọn thảm tươi ($H_{vn}$) và độ tàn che tán ($Dt$) để ước lượng nhanh sinh khối thảm tươi từ ảnh chụp viễn thám độ phân giải cao và công nghệ LiDAR gắn trên thiết bị bay không người lái (UAV).
- Mở rộng phạm vi điều tra đối chiếu sang các trạng thái rừng lân cận như rừng phục hồi giàu ($IIB$), rừng nghèo ($IIIA_1$) trên toàn địa bàn vùng đệm Vườn Quốc gia Tam Đảo.
Đối tượng hưởng lợi
- Sinh viên và học viên cao học ngành Lâm nghiệp, Quản lý Tài nguyên & Môi trường: Tiếp cận phương pháp luận chuẩn tắc về lập ô tiêu chuẩn lồng bậc, quy trình xử lý mẫu sấy khô kiệt và các công thức toán sinh xác định bể carbon sinh thái.
- Kỹ sư lâm nghiệp & Đơn vị tư vấn tín chỉ Carbon: Sử dụng các hệ số sinh khối khô ($3,337\text{ tấn/ha}$) và tỷ lệ phân bố trên/dưới mặt đất ($78,92% / 21,08%$) làm tham số đầu vào tin cậy để lập báo cáo khả thi PDD cho các dự án lâm nghiệp.
- Cơ quan quản lý nhà nước (Sở NN&PTNT, Chi cục Kiểm lâm Thái Nguyên): Nắm bắt dữ liệu định lượng để xây dựng chính sách giao khoán bảo vệ rừng, tích hợp cơ chế chi trả dịch vụ môi trường rừng gắn với giảm phát thải KNK cấp địa phương.
- Cộng đồng dân cư địa phương: Hưởng lợi từ nguồn quỹ bảo vệ rừng khi các chương trình tín chỉ carbon và PFES được giải ngân, tạo thêm sinh kế bền vững cho các hộ nhận khoán rừng tại 30 xóm thuộc xã Yên Lãng.
Câu hỏi thường gặp
1. Tại sao nghiên cứu lại chọn hệ số 0,5 để tính hàm lượng Carbon từ sinh khối khô?
Hệ số chuyển đổi $CF = 0.5$ là giá trị mặc định chuẩn quốc tế được công nhận bởi Ban Liên chính phủ về Biến đổi Khí hậu (IPCC - Good Practice Guidance for LULUCF). Về mặt sinh hóa, hàm lượng carbon nguyên tố chiếm từ $47% - 53%$ trọng lượng chất khô hữu cơ của các loài thực vật thân gỗ và thân thảo nhiệt đới.
2. Sự phân hóa vị trí địa hình (chân, sườn, đỉnh đồi) ảnh hưởng như thế nào đến lượng sinh khối tích lũy?
Kết quả nghiên cứu thực nghiệm cho thấy sinh khối khô biến động nhẹ giữa các đai: Sườn đồi đạt cao nhất ($3,491\text{ tấn/ha}$), tiếp theo là Đỉnh đồi ($3,277\text{ tấn/ha}$) và Chân đồi ($3,244\text{ tấn/ha}$). Khu vực sườn đồi có điều kiện thoát nước tốt và lượng ánh sáng tán xạ xuyên tán thuận lợi nhất cho nhóm cây bụi ưa sáng phát triển mật độ cao.
3. Công thức chuyển đổi từ lượng Carbon ($C$) sang lượng $CO_2$ tương đương ($CO_{2e}$) được thiết lập như thế nào?
Dựa trên tỷ lệ khối lượng phân tử hóa học của phân tử $CO_2$: Khối lượng nguyên tử Carbon ($C$) là 12, Oxy ($O$) là 16. Phân tử lượng của $CO_2$ là $12 + 16 \times 2 = 44$. Do đó, hệ số chuyển đổi là $\frac{44}{12} \approx 3,667$. Ta có:
$$\text{Lượng } CO_{2e} \text{ (tấn/ha)} = \text{Lượng } C \text{ (tấn/ha)} \times \frac{44}{12}$$
4. Tầng cây bụi thảm tươi có cản trở quá trình tái sinh tự nhiên của cây gỗ lớn không?
Với độ dày rậm đạt mức $Cop_2$ (che phủ $37,08%$), tầng thảm tươi có sự cạnh tranh dinh dưỡng và ánh sáng đối với hạt nảy mầm của cây gỗ mục đích. Tuy nhiên, chúng đóng vai trò then chốt che phủ bề mặt, giữ ẩm và ngăn xói mòn. Trong xử lý lâm sinh, cần áp dụng biện pháp phát dọn cục bộ quanh gốc cây tái sinh mục đích thay vì phát quang toàn diện nhằm bảo toàn bể carbon thảm tươi.
5. Dự án tín chỉ carbon có thể áp dụng mức giá nào cho rừng phục hồi IIA?
Trong nghiên cứu, mức giá giả định thận trọng là $5\text{ USD/tấn } CO_2$ (phù hợp với các giao dịch lâm nghiệp quốc tế giai đoạn khảo sát). Hiện nay, trên thị trường carbon tự nguyện (VCM), các tín chỉ chất lượng cao từ dự án phục hồi rừng tự nhiên kết hợp bảo tồn đa dạng sinh học (tiêu chuẩn CCB - Climate, Community & Biodiversity) có giá trị dao động từ $10 - 18\text{ USD/tấn } CO_2$, nâng cao đáng kể giá trị kinh tế của hệ sinh thái.
Kết luận
Công trình nghiên cứu đã giải quyết trọn vẹn và chuẩn xác bài toán định lượng sinh khối, trữ lượng carbon và tiềm năng thương mại của tầng cây bụi thảm tươi dưới tán rừng tự nhiên phục hồi trạng thái IIA tại xã Yên Lãng, huyện Đại Từ, tỉnh Thái Nguyên:
- Cấu trúc sinh thái ổn định: Định danh 13 loài thực vật thảm tươi ưu thế với mật độ dày rậm $37,08%$ ($Cop_2$), dẫn đầu là Mua rừng ($15,82%$) và Cỏ lá tre ($10,62%$).
- Trữ lượng sinh khối và carbon rõ nét: Sinh khối tươi đạt $10,264\text{ tấn/ha}$; Sinh khối khô tuyệt đối đạt $3,337\text{ tấn/ha}$; Lưu giữ trữ lượng carbon đạt $1,669\text{ tấn C/ha}$, tương đương khả năng hấp thụ $6,118\text{ tấn } CO_2\text{/ha}$.
- Giá trị kinh tế môi trường cụ thể: Mỗi hecta thảm tươi rừng IIA mang lại giá trị thương mại carbon tiềm năng đạt $30,59\text{ USD/ha}$.
Dữ liệu thực chứng từ khóa luận là cơ sở khoa học tin cậy để các nhà hoạch định chính sách, cơ quan quản lý lâm nghiệp và chủ rừng đưa bể carbon thảm tươi vào hệ thống kiểm kê phát thải, hoàn thiện hồ sơ dự án tín chỉ carbon và triển khai hiệu quả cơ chế chi trả dịch vụ môi trường rừng trong kỷ nguyên kinh tế xanh.
