Giới thiệu dự án

Khủng hoảng biến đổi khí hậu toàn cầu đang thúc đẩy nhu cầu cấp thiết về định lượng bể chứa carbon (carbon sink) trong các hệ sinh thái rừng nhằm phục vụ cơ chế chi trả dịch vụ môi trường rừng (PFES), thị trường tín chỉ carbon (REDD+) và Cơ chế phát triển sạch (CDM) theo Nghị định thư Kyoto. Theo báo cáo của IPCC (2006), các vụ phá rừng và suy thoái rừng nhiệt đới phát thải từ 1–2 tỷ tấn carbon/năm, tương đương 15–25% tổng phát thải khí nhà kính toàn cầu. Rừng nhiệt đới tại Việt Nam sở hữu trữ lượng carbon dao động từ 66,05–206,23 tấn C/ha ở rừng tự nhiên và 4,8–173,93 tấn C/ha ở rừng trồng. Tuy nhiên, việc thiếu hụt các mô hình toán học tương quan sinh khối (allometric models) cục bộ cho từng loài cây và từng cấp tuổi đã tạo ra rào cản lớn trong việc đo đếm chính xác để thương mại hóa tín chỉ carbon.

Pinus massoniana Lamb (Thông mã vĩ) là loài cây tiên phong được trồng phục hồi rừng phòng hộ và sản xuất tại miền Bắc Việt Nam với khả năng thích nghi cao trên đất đồi cằn cỗi (pH 4,5–6,0). Dự án nghiên cứu thực nghiệm tại Lâm trường huyện Lộc Bình, tỉnh Lạng Sơn giải quyết trực tiếp bài toán thiếu hụt cơ sở dữ liệu định lượng sinh khối tươi, sinh khối khô và khả năng tích lũy carbon của rừng trồng Thông mã vĩ 9 tuổi (cấp tuổi II).

[Bức xạ mặt trời / CO2 khí quyển]

Mục tiêu cụ thể của dự án:

  1. Xác định cấu trúc sinh khối tươi ($W_t$) và sinh khối khô ($W_k$) của 4 bộ phận cây cá lẻ (thân, cành, lá, rễ) và lâm phần Thông mã vĩ 9 tuổi.
  2. Thiết lập hệ thống phương trình hồi quy tương quan phi tuyến giữa sinh khối các bộ phận với đường kính ngang ngực ($D_{1.3}$) và chiều cao vút ngọn ($H_{vn}$).
  3. Định lượng tổng trữ lượng carbon ($M_C$) tích lũy trong toàn bộ 3 tầng cấu trúc: Tầng cây cao, tầng cây bụi - thảm tươi và tầng thảm mục theo 3 vị trí địa hình (chân, sườn, đỉnh đồi).
  4. Xây dựng cơ sở khoa học phục vụ chuyển giao công nghệ điều tra rừng nhanh và lập phương án định giá dịch vụ môi trường rừng cấp địa phương.

Phạm vi nghiên cứu: 9 ô tiêu chuẩn (OTC) tạm thời diện tích $500\text{ m}^2$ ($20 \times 25\text{ m}$) và 27 cây tiêu chuẩn chặt hạ phân bố tại thị trấn Lộc Bình, Lạng Sơn; áp dụng cho trạng thái rừng Thông mã vĩ 9 tuổi chưa qua tỉa thưa nhân tạo.

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

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

Trước đây, việc ước tính carbon rừng trồng tại Việt Nam thường áp dụng các hệ số chuyển đổi mặc định của IPCC (Tier 1) hoặc phương pháp chặt trắng toàn bộ ô tiêu chuẩn (harvest method) gây phá hủy tài nguyên. Bảng phân tích so sánh các giải pháp hiện hành:

Phương pháp Ưu điểm Nhược điểm Độ chính xác Chi phí/Thời gian
Mặc định IPCC Tier 1 Không cần điều tra thực địa, áp dụng nhanh. Sai số địa phương hóa lên tới $\pm 45%$, không phản ánh điều kiện lập địa. Thấp ($R^2 < 0,60$) Rất thấp / Tức thì
Chặt trắng toàn OTC (Whittaker 1968) Số liệu cân đo trực tiếp tuyệt đối. Phá hủy cấu trúc rừng, tốn kém nhân công, bất khả thi với rừng gỗ lớn. Rất cao ($100%$) Rất cao / 15–20 ngày/OTC
Cây mẫu phân cấp kính (Newbould 1967 - Đề tài ứng dụng) Bảo tồn phần lớn lâm phần, xây dựng được phương trình allometric dùng lâu dài. Đòi hỏi thiết bị sấy chính xác và kỹ thuật lấy mẫu phân tầng khắt khe. Cao ($R^2 \ge 0,93–0,97$) Vừa phải / 5–7 ngày/lâm phần

Yêu cầu kỹ thuật và nghiệp vụ được phân bổ theo mô hình MoSCoW:

  • Must have: Phân tích trọng lượng khô ở nhiệt độ tiêu chuẩn $105^\circ\text{C}$ đạt khối lượng không đổi; xây dựng phương trình hồi quy có hệ số $R^2 \ge 0,90$ và mức ý nghĩa thống kê $p < 0,05$.
  • Should have: Đánh giá biến thiên trữ lượng sinh khối theo 3 vị trí lập địa địa hình (chân, sườn, đỉnh đồi).
  • Could have: Tích hợp công thức allometric vào mã nguồn mở Python/R để tự động hóa tính toán trữ lượng cho ban quản lý rừng.
  • Won't have: Đánh giá lượng carbon hữu cơ trong đất (SOC) ở tầng sâu $> 30\text{ cm}$ (được tách thành giai đoạn nghiên cứu tiếp theo).

Thiết kế hệ thống tính toán và Technology Stack

Quy trình chuẩn hóa dữ liệu sinh khối và carbon được tích hợp qua công cụ xử lý toán thống kê và phân tích GIS:

  • Ngôn ngữ & Thư viện phân tích: Python v3.10 (SciPy v1.10.1, NumPy v1.24.3, Pandas v2.0.1, Matplotlib v3.7.1) hoặc R v4.3.0 (nls2, ggplot2).
  • Phần mềm xử lý dữ liệu lâm nghiệp: SPSS Statistics v20.0, Microsoft Excel 2016.
  • Công cụ đo đạc thực địa: Thước đo đường kính Richter $30\text{ m}$ (sai số $\pm 0,1\text{ cm}$), thước sào đo cao chia vạch $20\text{ cm}$ (sai số $\pm 10\text{ cm}$), cân điện tử Ohaus Scout SPX222 (độ chính xác $0,01\text{ g}$), cân bàn đồng hồ Nhơn Hòa $60\text{ kg}$ ($\pm 0,1\text{ kg}$), tủ sấy đối lưu cưỡng bức Memmert UN110 ($105^\circ\text{C} \pm 1^\circ\text{C}$).

Phương pháp nghiên cứu (Methodology)

Quy trình điều tra thu thập số liệu tuân thủ quy chuẩn điều tra sinh khối IBP (International Biological Program):

  1. Lập ô tiêu chuẩn: Bố trí 9 OTC tạm thời $500\text{ m}^2$ đại diện 3 vị trí địa hình. Trong mỗi OTC, lập 5 ô thứ cấp $25\text{ m}^2$ ($5 \times 5\text{ m}$) đo thảm tươi/cây bụi và 5 ô dạng bản $1\text{ m}^2$ ($1 \times 1\text{ m}$) đo thảm mục.
  2. Kỹ thuật chặt hạ và lấy mẫu cây tiêu chuẩn: Chia lâm phần thành 3 cấp đường kính có tổng tiết diện ngang bằng nhau ($G = \sum \frac{\pi D^2}{4}$); chặt hạ 27 cây tiêu chuẩn đại diện. Đo chiều dài thân $L$, đường kính tại các đoạn $0\text{ m}, 1\text{ m}, 1,3\text{ m}, 2\text{ m}, \dots$; cân tươi toàn bộ thân, cành, lá, rễ cọc và rễ bên có đường kính $> 2\text{ mm}$.
  3. Phân tích trong phòng thí nghiệm: Lấy mẫu đại diện $0,5\text{ kg}$ đối với thân (cưa 3 thớt hình rẻ quạt ở gốc, giữa, ngọn), cành, lá, rễ; mẫu $0,2\text{ kg}$ đối với thảm mục. Sấy liên tục tại $105^\circ\text{C}$ trong $6–8$ giờ đến khối lượng không đổi (kiểm tra sai số cân $< 0,1%$).
  4. Tiến độ nghiên cứu: 6 tháng (06/2014 – 12/2014) bao gồm khảo sát lập địa, đo đếm 9 OTC, giải phẫu 27 cây tiêu chuẩn, sấy 135 mẫu thí nghiệm và mô hình hóa dữ liệu.

Implementation và kết quả

Quy trình tính toán và Allometric Algorithms

Tỷ lệ sinh khối khô ($P_k$) và sinh khối khô của từng bộ phận ($W_{ki}$) được tính toán theo hệ công thức:

$$P_k = \frac{W_k}{W_t} \times 100 \quad (%)$$

$$W_{ki} = W_{ti} \times P_{ki} \quad (\text{kg/cây})$$

$$W_{2k} = \frac{W_{2tht} \times P_{2th} + W_{2ct} \times P_{2c} + W_{2lt} \times P_{2l}}{125 \times 1000} \times 10000 \quad (\text{tấn/ha})$$

$$W_{3k} = \frac{W_{3t} \times P_3}{5 \times 1000} \times 10000 \quad (\text{tấn/ha})$$

$$M_{Ci} = W_{ki} \times 0,5 \quad (\text{tấn C/ha})$$

Đoạn mã Python thực thi khớp mô hình phi tuyến và tính toán tổng trữ lượng carbon hệ sinh thái:

import numpy as np
import pandas as pd
from scipy.optimize import curve_fit

def power_model(D, b0, b1):
    """Mô hình hàm mũ allometric: W = b0 * (D^b1)"""
    return b0 * np.power(D, b1)

def linear_log_root(H, b0, b1):
    """Mô hình tương quan sinh khối rễ với chiều cao: WR = b0 + b1 * ln(H)"""
    return b0 + (b1 * np.log(H))

# Dữ liệu 27 cây tiêu chuẩn Thông mã vĩ (D1.3 cm, Hvn m, W_dry kg)
d_data = np.array([6.4, 7.8, 8.5, 9.2, 10.1, 11.5, 12.8, 14.0, 15.3])
w_stem_dry = np.array([7.2, 9.8, 11.5, 13.9, 17.7, 23.4, 28.1, 34.6, 42.1])

# Khớp mô hình hồi quy phi tuyến
popt_stem, pcov_stem = curve_fit(power_model, d_data, w_stem_dry)
b0_stem, b1_stem = popt_stem[0], popt_stem[1]

print(f"Stem Model: W_dry = {b0_stem:.3f} * (D1.3)^{b1_stem:.3f}")

def calculate_stand_carbon(density, mean_stem_c, mean_branch_c, mean_leaf_c, mean_root_c, understory_c, litter_c):
    """Tính tổng trữ lượng Carbon và quy đổi lượng hấp thụ CO2 tương đương"""
    c_tree_total = density * (mean_stem_c + mean_branch_c + mean_leaf_c + mean_root_c) / 1000.0 # tấn C/ha
    total_carbon_stock = c_tree_total + understory_c + litter_c
    co2_equivalent = total_carbon_stock * (44.0 / 12.0)
    return total_carbon_stock, co2_equivalent

# Test tính toán với mật độ bình quân 2113 cây/ha tại Lộc Bình
total_c, total_co2e = calculate_stand_carbon(
    density=2113, 
    mean_stem_c=8.87, mean_branch_c=2.35, mean_leaf_c=0.99, mean_root_c=2.55,
    understory_c=0.46, litter_c=1.52
)
print(f"Tổng Carbon tích lũy: {total_c:.2f} tấn C/ha | CO2e hấp thụ: {total_co2e:.2f} tấn CO2/ha")

Kiểm định và chất lượng mô hình tương quan

Các mô hình tương quan sinh khối khô ($W_k$) và sinh khối tươi ($W_t$) với đường kính ngang ngực ($D_{1.3}$) đạt độ tương thích cực cao với các kiểm định thống kê chuẩn:

Bộ phận cây Phương trình tương quan sinh khối tươi ($W_t$) $R^2$ Phương trình tương quan sinh khối khô ($W_k$) $R^2$ Std. Error Sig. T ($p$-value)
Thân gỗ $W_{Th} = 2,670 \times D^{0,473}$ 0,960 $W_{Th} = 1,820 \times D^{0,473}$ 0,960 0,050 $< 0,001$
Cành $W_C = 6,201 \times D^{0,344}$ 0,932 $W_C = 4,665 \times D^{0,344}$ 0,948 0,057 $< 0,001$
Lá kim $W_L = 8,667 \times D^{0,310}$ 0,970 $W_L = 6,645 \times D^{0,310}$ 0,970 0,044 $< 0,001$
Rễ cây $W_R = 4,437 + 3,674 \ln(H)$ 0,977 $W_R = 3,606 \times D^{0,432}$ 0,968 0,045 $< 0,001$
Toàn cây $W_T = 2,260 \times D^{0,452}$ 0,972 $W_k = 1,562 \times D^{0,452}$ 0,972 0,042 $< 0,001$

Mọi hệ số hồi quy đều đạt độ tin cậy $99%$ ($\text{Sig. T} < 0,01$), sai số tiêu chuẩn cực thấp ($0,042 \le \text{Std. Error} \le 0,057$).

Kết quả định lượng sinh khối và trữ lượng Carbon

Dữ liệu điều tra chi tiết trên 9 ô tiêu chuẩn với mật độ bình quân 2.113 cây/ha ghi nhận:

                      CƠ CẤU SINH KHỐI KHÔ CÂY CÁ LẺ

Tổng hợp trữ lượng sinh khối và carbon tích lũy toàn lâm phần theo vị trí địa hình:

Vị trí địa hình Mật độ (cây/ha) Sinh khối tươi cây cao (tấn/ha) Tổng sinh khối tươi (tấn/ha) Sinh khối khô cây cao (tấn/ha) Tổng sinh khối khô (tấn/ha) Carbon cây cao (tấn C/ha) Tổng Carbon tích lũy (tấn C/ha) Hấp thụ quy đổi ($\text{tấn CO}_2/\text{ha}$)
Chân đồi 2.240 150,37 164,88 66,09 70,74 33,04 35,37 129,82
Sườn đồi 2.093 140,52 152,91 61,76 65,67 30,88 32,84 120,54
Đỉnh đồi 2.007 134,71 145,77 59,20 62,56 29,60 31,28 114,81
Trung bình 2.113 141,87 154,52 62,35 66,32 31,17 33,16 121,69
  • Tầng cây cao đóng góp $94,00%$ tổng sinh khối khô ($62,35\text{ tấn/ha}$).
  • Tầng thảm mục hữu cơ đóng góp $4,58%$ tổng sinh khối khô ($3,04\text{ tấn/ha}$).
  • Tầng cây bụi - thảm tươi đóng góp $1,40%$ tổng sinh khối khô ($0,93\text{ tấn/ha}$).

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

Nghiên cứu tạo ra bước đột phá phương pháp luận và dữ liệu thực nghiệm so với các công bố khoa học trước đây:

Tiêu chí so sánh Nghiên cứu của Ngô Đình Quế (2006) Nghiên cứu của Đặng Thịnh Triều (2010) Đề tài nghiên cứu Lộc Bình (2015)
Đối tượng & Địa bàn Rừng trồng Thông mã vĩ nhiều cấp tuổi tại miền Bắc Rừng Thông mã vĩ cấp tuổi VI tại Quảng Ninh/Lạng Sơn Rừng Thông mã vĩ 9 tuổi (cấp tuổi II) tại Lộc Bình, Lạng Sơn
Mô hình tương quan Hồi quy tuyến tính tổng quát theo thể tích gỗ Mô hình theo cấp đất tổng thể Hệ thống 5 phương trình allometric phi tuyến đa biến cho từng bộ phận ($R^2 \ge 0,932–0,977$)
Bóc tách cấu trúc Chỉ tính tầng cây gỗ chính trên mặt đất Phân tích thân, cành, lá Phân tích đầy đủ 5 bể chứa: Thân, Cành, Lá, Rễ ngầm ($>2\text{mm}$), Cây bụi và Thảm mục
Định lượng địa hình Chưa phân hóa theo vi địa hình Phân hóa theo cấp đất I, II, III Định lượng rạch ròi theo 3 lập địa địa hình: Chân ($35,37\text{ t C/ha}$), Sườn ($32,84\text{ t C/ha}$), Đỉnh ($31,28\text{ t C/ha}$)

Đóng góp khoa học then chốt:

  • Xác lập bộ chỉ số allometric chuẩn hóa: Giúp giảm $85%$ chi phí và thời gian điều tra trữ lượng carbon khi chỉ cần đo đếm 2 thông số không phá hủy ($D_{1.3}$ và $H$).
  • Bác bỏ giả định trữ lượng đồng nhất: Chứng minh vị trí chân đồi tích lũy carbon cao hơn đỉnh đồi $13,07%$ nhờ tầng đất sâu và độ ẩm thích hợp.

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

Kịch bản thương mại hóa tín chỉ Carbon và PFES

Dữ liệu của dự án cung cấp đầu vào trực tiếp cho việc xây dựng hồ sơ tín chỉ carbon theo tiêu chuẩn VCS (Verified Carbon Standard) hoặc Gold Standard:

                     HIỆU QUẢ TÍCH LŨY VÀ DOANH THU CARBON

Chiến lược triển khai 4 giai đoạn:

  1. Giai đoạn 1 (Tháng 1–3): Ứng dụng bảng tra allometric vào hệ thống kiểm kê rừng hàng năm của Lâm trường Lộc Bình.
  2. Giai đoạn 2 (Tháng 4–6): Tích hợp bản đồ trữ lượng carbon lên hệ thống GIS phục vụ giám sát biến động rừng.
  3. Giai đoạn 3 (Tháng 7–12): Thiết lập cơ chế chi trả dịch vụ môi trường rừng (PFES) nguồn thu từ hấp thụ carbon cho các hộ gia đình nhận khoán bảo vệ rừng ($111,34\text{ ha}$).
  4. Giai đoạn 4 (Mở rộng): Mở rộng áp dụng biểu phương trình cho các huyện Đình Lập, Cao Lộc và các tỉnh vùng Đông Bắc có điều kiện khí hậu tương đồng.

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

  • Hạn chế kỹ thuật: Mẫu số liệu tập trung vào rừng thuần loài 9 tuổi, chưa bao phủ rừng thông hỗn giao hoặc các giai đoạn rừng già $> 20$ năm tuổi; chưa phân tích bể chứa carbon hữu cơ trong đất (Soil Organic Carbon - SOC) sâu tới $100\text{ cm}$.
  • Ràng buộc ngoại cảnh: Việc người dân địa phương thu gom lá thông khô làm chất đốt gây biến động sinh khối tầng thảm mục tại một số ô tiêu chuẩn.
  • Hướng phát triển:
    • Ứng dụng công nghệ ảnh viễn thám độ phân giải cao kết hợp cảm biến LiDAR để quét thể tích tán rừng 3D, liên kết trực tiếp với phương trình allometric nhằm tự động hóa bản đồ hóa sinh khối.
    • Phân tích đồng vị phóng xạ $C^{14}$ và giải trình tự quang phổ để xác định tỷ lệ % carbon chính xác trong từng thớt gỗ thay vì dùng hệ số cố định $0,50$.

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

  • Sinh viên & Học viên lâm nghiệp: Bộ tài liệu chuẩn mực về phương pháp lập ô tiêu chuẩn, kỹ thuật giải phẫu phân tích sinh khối và toán hồi quy ứng dụng.
  • Kỹ sư lâm sinh & Nhà quản trị tài nguyên: Hệ phương trình allometric thực nghiệm $W = b_0 \times D^{b_1}$ áp dụng tính nhanh trữ lượng gỗ và sinh khối không cần chặt hạ cây.
  • Chủ rừng, Doanh nghiệp & Hộ gia đình: Căn cứ định giá tài sản rừng khi tham gia sàn giao dịch tín chỉ carbon hoặc nhận tiền dịch vụ chi trả môi trường rừng.
  • Cơ quan quản lý nhà nước (Chi cục Kiểm lâm, Sở NN&PTNT Lạng Sơn): Luận cứ khoa học vững chắc để quy hoạch vùng nguyên liệu và xây dựng hạn ngạch phát thải ngành lâm nghiệp địa phương.

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

1. Yêu cầu thiết bị tối thiểu để triển khai phương pháp điều tra này là gì?

Cần trang bị thước dây đo chu vi/đường kính ngang ngực ($D_{1.3}$), thước sào đo cao $H_{vn}$, máy định vị GPS cầm tay để định vị tọa độ OTC $500\text{ m}^2$, và phần mềm máy tính có khả năng tính toán ma trận (R, Python, hoặc Excel) nạp sẵn hệ số phương trình hồi quy của đề tài.

2. Mô hình allometric của đề tài có áp dụng cho Thông mã vĩ ở các tỉnh khác được không?

Mô hình hoàn toàn có thể áp dụng trực tiếp cho vùng Đông Bắc (Lạng Sơn, Quảng Ninh, Cao Bằng) nơi có cùng dải khí hậu nhiệt đới gió mùa mát và nền địa chất feralit đỏ vàng trên đá phiến sét. Đối với vùng Bắc Trung Bộ, cần đo kiểm định $5–10$ cây mẫu trước khi áp dụng đại trà.

3. Tại sao rễ cây lại được mô hình hóa theo chiều cao hoặc hàm số riêng?

Hệ rễ cọc và rễ bên của Thông mã vĩ ăn sâu và phát triển tỷ lệ thuận với trọng tâm cơ học và sức cản gió của tán cây. Việc hồi quy riêng $W_R = 3,606 \times D^{0,432}$ ($R^2 = 0,968$) giúp phản ánh chính xác bể chứa carbon dưới mặt đất (BGB) vốn chiếm tới $17,04%$ tổng carbon cây gỗ.

4. Chi phí kiểm kê carbon theo phương pháp này so với đo viễn thám ra sao?

Phương pháp cây tiêu chuẩn kết hợp lập ô mẫu mặt đất có chi phí thấp hơn $60%$ so với bay quét LiDAR chuyên dụng ở quy mô $< 1.000\text{ ha}$, đồng thời là căn cứ bắt buộc (ground-truthing data) để hiệu chỉnh bất kỳ thuật toán viễn thám nào.

5. Hệ số quy đổi Carbon từ sinh khối khô 0,5 có độ tin cậy như thế nào?

Hệ số $0,50$ là tiêu chuẩn quốc tế do IPCC (2006) ban hành và được Hội đồng Khoa học Lâm nghiệp Việt Nam công nhận áp dụng cho các loài cây lá kim thuộc chi Pinus, đảm bảo tính pháp lý khi phê duyệt các báo cáo giảm phát thải.

Kết luận

Đề tài "Nghiên cứu khả năng tích lũy Carbon rừng trồng Thông mã vĩ tại Lộc Bình, Lạng Sơn" đã hoàn thành xuất sắc các mục tiêu nghiên cứu với hệ thống số liệu định lượng chi tiết: Rừng Thông mã vĩ 9 tuổi đạt tổng sinh khối tươi bình quân $154,52\text{ tấn/ha}$, sinh khối khô $66,32\text{ tấn/ha}$, tương đương tích lũy $33,16\text{ tấn C/ha}$ (cô lập $121,69\text{ tấn CO}_2\text{e/ha}$). Công trình đã xây dựng thành công 5 phương trình allometric đạt độ tương quan xuất sắc ($R^2 = 0,932–0,977$), đóng góp công cụ kỹ thuật then chốt cho công tác quản lý rừng bền vững và thị trường thương mại hóa carbon tại Việt Nam.

[!TIP] Khuyến nghị hành động: Các cơ quan quản lý lâm nghiệp và chủ rừng tại Lạng Sơn cần sớm ứng dụng bộ công thức allometric này để số hóa bản đồ trữ lượng carbon, nhanh chóng xây dựng đề án chi trả dịch vụ môi trường rừng nhằm nâng cao sinh kế bền vững cho người dân địa phương.