Giới thiệu dự án

Biến đổi khí hậu toàn cầu do sự gia tăng nồng độ khí nhà kính (GHG) – đặc biệt là $\text{CO}_2$ với mức phát thải vượt ngưỡng 36 tỷ tấn/năm trên toàn cầu – đang thúc đẩy nhu cầu cấp bách về các giải pháp lưu trữ carbon tự nhiên. Theo Công ước Khung của Liên Hợp Quốc về Biến đổi Khí hậu (UNFCCC) và Cơ chế Phát triển Sạch (Clean Development Mechanism - CDM) theo Nghị định thư Kyoto, các hệ sinh thái rừng đóng vai trò là "bể chứa" hấp thụ khoảng 100 tỷ tấn $\text{CO}_2$ mỗi năm thông qua quá trình quang hợp. Tại Việt Nam, Thông mã vĩ (Pinus massoniana Lamb) là loài cây gỗ tiên phong được trồng rộng rãi tại vùng Đông Bắc (Lạng Sơn, Quảng Ninh, Cao Bằng) nhằm phục hồi đất trống đồi núi trọc và khai thác nhựa, gỗ. Tuy nhiên, việc lượng hóa giá trị sinh thái và khả năng thương mại hóa tín chỉ carbon của loài cây này vẫn còn thiếu các mô hình toán học thực nghiệm chính xác theo từng cấp tuổi và điều kiện lập địa cụ thể.

                    HỆ SINH THÁI RỪNG THÔNG MÃ VĨ (TUỔI 9)
   SINH KHỐI TRÊN MẶT ĐẤT                              SINH KHỐI DƯỚI MẶT ĐẤT
   (83,73% Sinh khối khô)                              (16,27% Sinh khối khô)
 Thân     Cành       Lá       Cây bụi/Thảm mục              Hệ rễ cây
(60,05%) (15,65%)  (8,03%)    (1,25% / 0,33%)               (> 2mm)

Vấn đề thực tiễn và mục tiêu dự án

Rừng trồng Thông mã vĩ tại xã Hữu Khánh, huyện Lộc Bình, tỉnh Lạng Sơn chủ yếu được quản lý theo phương thức truyền thống lấy gỗ và nhựa, chưa có cơ sở dữ liệu định lượng về sinh khối khô ($W_k$) và trữ lượng carbon ($W_C$). Thiếu hụt các phương trình tương quan sinh khối (allometric equations) địa phương dẫn đến sai số ước tính trữ lượng carbon lên tới 25-30% khi áp dụng các mô hình tổng quát.

Đề tài đặt ra 4 mục tiêu cụ thể:

  1. Xác định chỉ tiêu sinh trưởng và cấu trúc thể tích của 27 cây tiêu chuẩn đại diện cho 3 cấp kính tại lâm phần Thông mã vĩ thuần loài 9 tuổi.
  2. Định lượng sinh khối tươi ($W_t$) và sinh khối khô ($W_k$) chi tiết theo 4 bộ phận cây cá lẻ (thân, cành, lá, rễ) và các tầng thảm mục, cây bụi thảm tươi.
  3. Xây dựng hệ thống mô hình hồi quy tương quan tối ưu giữa sinh khối, carbon tích lũy với biến tổ hợp điều tra ($D_{1.3}^2 \cdot H_{vn}$).
  4. Lượng hóa tổng trữ lượng tích lũy carbon ($W_C$) của lâm phần trên 3 dạng địa hình (chân đồi, sườn đồi, đỉnh đồi) phục vụ mục tiêu chi trả dịch vụ môi trường rừng (PFES) và thị trường carbon.

Phạm vi và giới hạn nghiên cứu

  • Đối tượng: Rừng trồng Thông mã vĩ thuần loài 9 tuổi, mật độ trung bình 2.104 cây/ha, chưa qua tỉa thưa.
  • Không gian: 9 ô tiêu chuẩn ($500\text{ m}^2/\text{ô}$) bố trí tại 3 vị trí địa hình (chân, sườn, đỉnh) trên nền đất đỏ vàng Feralit phát triển trên đá sét (pH 4,5 - 5,5), độ dốc $>25^\circ$.
  • Giới hạn kỹ thuật: Nghiên cứu tập trung vào sinh khối sống (trên và dưới mặt đất) cùng thảm mục, chưa bao gồm phân tích hàm lượng carbon hữu cơ trong đất (SOC) ở các tầng sâu $>30\text{ cm}$.

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

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

Hiện nay, việc đánh giá sinh khối và carbon rừng dựa trên ba hướng tiếp cận chính với ưu nhược điểm rõ rệt:

Phương pháp Ưu điểm Nhược điểm Chi phí & Khả thi
Đo đếm chặt hạ trực tiếp (Chặt cây mẫu) Độ chính xác cao tuyệt đối ($\pm 1-2%$), phân tích chi tiết từng bộ phận ($W_t, W_k$). Phá hủy mẫu, tốn nhân lực và thời gian thực địa. Chi phí cao, là tiêu chuẩn vàng để lập mô hình.
Hệ số mở rộng sinh khối (BEF / IPCC Tier 1) Áp dụng nhanh từ số liệu thể tích gỗ thương phẩm. Sai số lớn ($20 - 40%$), không phản ánh được sai khác lập địa. Chi phí thấp, độ tin cậy hạn chế cho dự án tín chỉ.
Viễn thám kết hợp LiDAR / UAV (Tier 3) Bao phủ diện tích lớn, giám sát định kỳ liên tục. Cần dữ liệu mặt đất (ground-truth) chuẩn xác để hiệu chuẩn. Đầu tư công nghệ cao, phụ thuộc mô hình thực địa.

Nghiên cứu này lựa chọn phương pháp chặt hạ cây tiêu chuẩn phân tầng theo cấp kính kết hợp phân tích phòng thí nghiệm nhằm tạo ra bộ phương trình allometric chuẩn xác làm "ground-truth" cho khu vực Lạng Sơn.

          QUY TRÌNH THU THẬP VÀ XỬ LÝ SỐ LIỆU SINH KHỐI
          
  [5 Ô thứ cấp 25m2 & 1m2]       [Chặt hạ 27 cây tiêu chuẩn]
                                 [Xác định tỷ lệ Wk & Model hóa]

Thiết kế hệ thống mô hình toán học

Hệ thống xử lý thống kê sử dụng phần mềm IBM SPSS Statistics v20.0 và Microsoft Excel 2013 để đánh giá 5 dạng phương trình hồi quy phi tuyến và tuyến tính:

  • Tuyến tính: $Y = a + bX$
  • Đa thức bậc 2: $Y = aX^2 + bX + c$
  • Logarit: $Y = a + b \cdot \ln(X)$
  • Hàm mũ: $Y = a \cdot e^{bX}$
  • Hàm lũy thừa (Power Model - Tối ưu): $Y = a \cdot X^b$ (với $X = D_{1.3}^2 \cdot H_{vn}$)
import numpy as np
from scipy.optimize import curve_fit

def allometric_power_model(D, H, a, b):
    """
    Tính toán sinh khối khô (Wk) dựa trên đường kính D1.3 (m) và chiều cao Hvn (m).
    Mô hình: Wk = a * (D^2 * H)^b
    """
    D2H = (D ** 2) * H
    return a * (D2H ** b)

# Tham số đã được kiểm định từ 27 cây mẫu lâm phần Thông mã vĩ tuổi 9
params = {
    'stem': {'a': 38.606, 'b': 1.120, 'R2': 0.981},
    'branch': {'a': 230.715, 'b': 0.925, 'R2': 0.954},
    'root': {'a': 175.221, 'b': 1.085, 'R2': 0.982},
    'total': {'a': 22.945, 'b': 1.103, 'R2': 0.983}
}

def calculate_carbon_stock(dw_total_ton_per_ha, cf=0.5):
    """Quy đổi sinh khối khô kiệt sang trữ lượng Carbon và CO2 tương đương"""
    carbon_ton = dw_total_ton_per_ha * cf
    co2_equivalent = carbon_ton * (44.0 / 12.0)
    return carbon_ton, co2_equivalent

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

  1. Thiết lập ô tiêu chuẩn (OTC): Bố trí 9 OTC tạm thời kích thước $20\text{ m} \times 25\text{ m}$ ($500\text{ m}^2$) phân bổ đều ở 3 vị trí chân, sườn, đỉnh.
  2. Điều tra tầng cây cao: Đo $D_{1.3}$ bằng thước dây chuyên dụng, $H_{vn}$ và chiều cao dưới cành ($H_{dc}$) bằng sào đo cao $\pm 10\text{ cm}$, đo đường kính tán ($D_t$) theo 2 hướng Đông - Tây và Nam - Bắc.
  3. Điều tra cây bụi thảm tươi và thảm mục: Lập 5 ô thứ cấp $25\text{ m}^2$ ($5\text{ m} \times 5\text{ m}$) tại 4 góc và tâm OTC; 5 ô dạng bản $1\text{ m}^2$ ($1\text{ m} \times 1\text{ m}$) ở tâm ô thứ cấp để thu thập thảm mục.
  4. Phương pháp sấy khô kiệt: Mẫu đại diện ($0,5\text{ kg}$ mỗi bộ phận) được sấy trong tủ sấy chuyên dụng ở nhiệt độ $105^\circ\text{C}$ trong $10 - 12\text{ giờ}$ cho đến khi khối lượng không đổi giữa 4 lần cân liên tiếp trên cân phân tích điện tử độ chính xác $0,1%$.
  5. Tính toán chuyển đổi Carbon: Áp dụng hệ số IPCC (2003) chuẩn $C = 0,5 \times DW_{\text{trên mặt đất}}$.

Implementation và kết quả

Quá trình phân tích thực nghiệm cây tiêu chuẩn

Điều tra thực địa 27 cây tiêu chuẩn chặt hạ phân theo 3 cấp kính (Cấp I: $D_{1.3} < 8\text{ cm}$; Cấp II: $8\text{ cm} \le D_{1.3} \le 10\text{ cm}$; Cấp III: $D_{1.3} > 10\text{ cm}$) thu được các thông số cơ bản:

Cấp kính $D_{1.3}$ trung bình (cm) $H_{vn}$ trung bình (m) Tiết diện ngang $G$ ($\text{m}^2$) Thể tích $V$ ($\text{m}^3$/cây)
Cấp I 7,60 7,20 0,0046 0,0166
Cấp II 9,70 8,05 0,0075 0,0306
Cấp III 12,48 9,04 0,0124 0,0566
Bình quân lâm phần 9,93 8,10 0,0082 0,0346

Hệ thống mô hình tương quan sinh khối Allometric

Kết quả hồi quy phi tuyến giữa sinh khối tươi ($W_t$), sinh khối khô ($W_k$) của các bộ phận với biến tổ hợp $D_{1.3}^2 \cdot H_{vn}$ ($D$ tính bằng mét, $H$ tính bằng mét) cho độ tương quan đặc biệt cao:

Bộ phận cây Phương trình sinh khối tươi ($W_t$, kg) $R^2$ Phương trình sinh khối khô ($W_k$, kg) $R^2$ Sai tiêu chuẩn (Std.E) Mức ý nghĩa ($p$-value)
Thân $W_t = 15,235 \cdot (D^2 \cdot H)^{1,120}$ 0,981 $W_k = 38,606 \cdot (D^2 \cdot H)^{1,120}$ 0,981 0,084 $< 0,001$
Cành $W_t = 103,849 \cdot (D^2 \cdot H)^{0,925}$ 0,954 $W_k = 230,715 \cdot (D^2 \cdot H)^{0,925}$ 0,954 0,126 $< 0,001$
Lá $\ln(W_t) = \ln(287,974) + 0,178(D^2 \cdot H)$ 0,974 $\ln(W_k) = \ln(287,973) + 0,435(D^2 \cdot H)$ 0,974 0,095 $< 0,001$
Rễ $W_t = 67,684 \cdot (D^2 \cdot H)^{1,085}$ 0,982 $W_k = 175,221 \cdot (D^2 \cdot H)^{1,085}$ 0,982 0,080 $< 0,001$
Toàn cây $W_t = 9,043 \cdot (D^2 \cdot H)^{1,102}$ 0,983 $W_k = 22,945 \cdot (D^2 \cdot H)^{1,103}$ 0,983 0,077 $< 0,001$

Kết quả sinh khối và trữ lượng Carbon lâm phần

  1. Cấu trúc sinh khối cây cá lẻ:
    • Sinh khối tươi bình quân: $62,99\text{ kg/cây}$ (Thân $58,97%$, Rễ $16,74%$, Cành $15,87%$, Lá $8,42%$).
    • Sinh khối khô bình quân: $26,96\text{ kg/cây}$ (Thân $16,09\text{ kg} = 60,05%$, Rễ $4,40\text{ kg} = 16,27%$, Cành $4,26\text{ kg} = 15,65%$, Lá $2,20\text{ kg} = 8,03%$). Tỷ lệ sinh khối khô/tươi đạt $42,8%$.
  2. Sinh khối lâm phần theo lập địa (Tấn/ha):
    • Chân đồi: Tổng sinh khối tươi đạt cao nhất với $139,40\text{ tấn/ha}$ (Tầng cây cao $136,89\text{ t/ha}$; Cây bụi $2,11\text{ t/ha}$; Thảm mục $0,40\text{ t/ha}$).
    • Sườn đồi: Đạt $132,88\text{ tấn/ha}$ (Tầng cây cao $130,59\text{ t/ha}$; Cây bụi $1,82\text{ t/ha}$; Thảm mục $0,47\text{ t/ha}$).
    • Đỉnh đồi: Đạt $131,76\text{ tấn/ha}$ (Tầng cây cao $130,17\text{ t/ha}$; Cây bụi $1,13\text{ t/ha}$; Thảm mục $0,46\text{ t/ha}$).
    • Bình quân chung: Đạt $134,68\text{ tấn/ha}$ tươi, tương ứng sinh khối khô tầng cây gỗ là $56,73\text{ tấn/ha}$.
  3. Trữ lượng tích lũy Carbon ($W_C$):
    • Trữ lượng carbon tầng cây gỗ bình quân đạt $28,37\text{ tấn C/ha}$ (tương đương hấp thụ cố định $104,02\text{ tấn CO}_2\text{/ha}$).
    • Phân bố carbon theo bộ phận: Thân ($17,04\text{ tấn C/ha}$), Rễ ($4,62\text{ tấn C/ha}$), Cành ($4,44\text{ tấn C/ha}$), Lá ($2,28\text{ tấn C/ha}$).
PHÂN BỐ TRỮ LƯỢNG CARBON THEO BỘ PHẬN TẦNG CÂY GỖ (TỔNG: 28,37 TẤN C/HA)

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

  1. Xây dựng bộ phương trình Allometric đặc thù cho cấp tuổi 9 tại vùng Đông Bắc: Thay vì áp dụng phương trình chung của IPCC vốn có biên độ sai lệch lớn, mô hình lũy thừa $W_k = 22,945 \cdot (D^2 \cdot H)^{1,103}$ đạt $R^2 = 0,983$, giúp giảm sai số ước tính sinh khối địa phương từ $>20%$ xuống dưới $3,5%$.
  2. Lượng hóa thực nghiệm hệ rễ cây ngầm ($>2\text{ mm}$): Trực tiếp đào và cân đo hệ rễ của 27 cây mẫu, chứng minh rễ chiếm tới $16,27%$ tổng sinh khối khô ($9,23\text{ tấn khô/ha}$). Con số này bác bỏ nhận định trước đây cho rằng rễ thông dưới 10 tuổi chỉ chiếm $<10%$.
  3. So sánh với các công trình nghiên cứu tương đồng:
    • So với nghiên cứu của Đặng Thịnh Triều (2008) trên Thông mã vĩ: Tỷ lệ phân bố thân ($60,05%$ so với $62,85%$) và cành ($15,65%$ so với $18,35%$) có tính đồng nhất cao, phản ánh chính xác quy luật sinh trưởng giai đoạn trước tỉa thưa.
    • So với rừng Thông ba lá (Pinus kesiya) của Lê Hồng Phúc (1996): Tỷ lệ sinh khối lá của Thông mã vĩ thấp hơn ($8,03%$ so với $11,2%$), cho thấy đặc tính rụng lá tập trung và hiệu suất tích lũy sinh khối thân cao hơn.
  4. Cơ sở định giá dịch vụ hấp thụ carbon: Cung cấp định lượng $104,02\text{ tấn CO}_2\text{/ha}$ cho rừng 9 tuổi làm căn cứ khoa học tính toán đơn giá chi trả tín chỉ carbon lâm nghiệp.

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

Kịch bản thương mại hóa và ứng dụng thực tiễn

  • Thị trường Tín chỉ Carbon Rừng: Với diện tích rừng Thông mã vĩ tại Lộc Bình ước tính hàng nghìn hecta, tại mức giá thị trường carbon tự nguyện (VCM) dao động từ $5 - 10\text{ USD/tấn CO}_2$, mỗi hecta rừng Thông tuổi 9 tích lũy giá trị môi trường tương đương $520 - 1.040\text{ USD}$ bên cạnh nguồn thu từ nhựa và gỗ.
  • Tối ưu hóa quản trị lâm sinh: Dữ liệu cho thấy ở vị trí chân đồi sinh khối đạt $139,40\text{ tấn/ha}$, cao hơn sườn ($132,88\text{ t/ha}$) và đỉnh ($131,76\text{ t/ha}$). Do đó, các biện pháp bón phân bổ sung và tỉa thưa cần tập trung vào các đai sườn và đỉnh đồi để nâng cao đồng đều trữ lượng rừng.
                    LỘ TRÌNH TRIỂN KHAI ỨNG DỤNG

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

Hạn chế kỹ thuật

  • Đơn cấp tuổi: Nghiên cứu chỉ thực hiện trên lâm phần thuần loài tuổi 9, chưa bao quát toàn bộ chu kỳ kinh doanh (từ tuổi 1 đến tuổi 25-30).
  • Hệ số quy đổi carbon mặc định: Nghiên cứu áp dụng hệ số chuyển đổi carbon $C=0,5$ của IPCC (2003) thay vì phân tích đốt cháy trực tiếp mẫu bằng máy phân tích nguyên tố CHN Elemental Analyzer.
  • Carbon trong đất: Chưa tính toán lượng carbon hữu cơ tích lũy trong đất rừng (SOC), vốn có thể chiếm từ $30 - 40%$ tổng trữ lượng carbon của hệ sinh thái.

Hướng nghiên cứu tiếp theo

  • Xây dựng chuỗi số liệu thời gian (chronosequence) theo các cấp tuổi 5, 10, 15, 20 và 25 tuổi.
  • Kết hợp dữ liệu bay chụp UAV LiDAR mật độ cao để ngoại suy sinh khối quy mô cấp vùng huyện/tỉnh.
  • Đo đạc chính xác tỷ lệ phát thải khí nhà kính từ quá trình phân hủy thảm mục theo mùa.

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

  • Sinh viên & Học viên Lâm nghiệp/Môi trường: Tài liệu tham khảo chuẩn mực về phương pháp chọn cây tiêu chuẩn phân tầng cấp kính, kỹ thuật cân sấy mẫu và xây dựng mô hình allometric.
  • Kỹ sư lâm nghiệp & Đơn vị điều tra quy hoạch rừng: Sử dụng trực tiếp biểu phương trình $W_k = 22,945 \cdot (D^2 \cdot H)^{1,103}$ để tính nhanh sinh khối mà không cần chặt hạ cây.
  • Chủ rừng & Doanh nghiệp phát triển dự án Carbon: Có số liệu đầu vào chuẩn xác để lập hồ sơ thiết kế dự án (PDD) theo cơ chế CDM hoặc Tiêu chuẩn Carbon Đã Xác Minh (VCS).
  • Cơ quan quản lý Nhà nước (Chi cục Kiểm lâm, Sở NN&PTNT): Căn cứ khoa học phục vụ chính sách chi trả dịch vụ môi trường rừng và định giá tài nguyên rừng quốc gia.

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

1. Dữ liệu đầu vào tối thiểu để áp dụng mô hình tính nhanh sinh khối là gì?

Người điều tra chỉ cần đo 2 chỉ tiêu lâm sinh cơ bản không phá hủy ngoài thực địa: Đường kính ngang ngực $D_{1.3}$ (đổi ra đơn vị mét) và Chiều cao vút ngọn $H_{vn}$ (mét), sau đó đưa vào phương trình $W_k = 22,945 \cdot (D^2 \cdot H)^{1,103}$ để tính tổng sinh khối khô của cây ($kg$).

2. Tại sao mô hình chọn biến tổ hợp $D_{1.3}^2 \cdot H_{vn}$ thay vì chỉ dùng riêng $D_{1.3}$?

Biến tổ hợp $D^2 \cdot H$ tỷ lệ thuận trực tiếp với thể tích hình học của thân cây hình nón/trụ. Kết quả kiểm định cho thấy $D^2 \cdot H$ đạt hệ số $R^2 = 0,983$, cao hơn đáng kể so với mô hình đơn biến chỉ dùng $D_{1.3}$ ($R^2 \approx 0,91 - 0,93$), giúp giảm thiểu tối đa sai số do phân cành và biến động chiều cao theo điều kiện lập địa.

3. Tỷ lệ phần trăm sinh khối khô so với sinh khối tươi ở Thông mã vĩ là bao nhiêu?

Tỷ lệ sinh khối khô kiệt toàn cây bình quân đạt $42,8%$ ($26,96\text{ kg khô} / 62,99\text{ kg tươi}$). Trong đó: Thân cây đạt $43,6%$, Rễ đạt $41,5%$, Cành đạt $42,2%$ và Lá đạt $40,7%$.

4. Khả năng hấp thụ $\text{CO}_2$ của 1 ha rừng Thông mã vĩ tuổi 9 là bao nhiêu?

Tầng cây gỗ tích lũy trung bình $28,37\text{ tấn C/ha}$. Chuyển đổi theo hệ số phân tử $\text{CO}_2 / \text{C} = 44 / 12 \approx 3,667$, lượng $\text{CO}_2$ được cố định tương đương $104,02\text{ tấn CO}_2\text{/ha}$.

5. Tại sao có sự chênh lệch sinh khối giữa các vị trí chân đồi, sườn đồi và đỉnh đồi?

Vị trí chân đồi có tầng đất dày hơn, độ ẩm cao và lượng dinh dưỡng tích tụ do quá trình rửa trôi từ đỉnh dốc, giúp cây sinh trưởng tốt hơn ($139,40\text{ tấn tươi/ha}$). Ngược lại, đỉnh đồi bị xói mòn mạnh, đất mỏng và khô hạn hơn nên sinh khối thấp nhất ($131,76\text{ tấn tươi/ha}$).

Kết luận

Công trình nghiên cứu đã hoàn thành toàn diện việc định lượng cấu trúc sinh khối và khả năng tích lũy carbon của rừng trồng Thông mã vĩ (Pinus massoniana Lamb) thuần loài 9 tuổi tại xã Hữu Khánh, huyện Lộc Bình, tỉnh Lạng Sơn. Các kết quả then chốt bao gồm:

  • Thiết lập hệ thống 10 phương trình tương quan sinh khối allometric với độ tin cậy vượt trội ($R^2 > 0,95; p < 0,001$).
  • Xác định cấu trúc sinh khối khô cây cá thể tập trung chủ yếu ở Thân ($60,05%$) và Hệ rễ ($16,27%$).
  • Định lượng trữ lượng carbon tầng cây gỗ đạt $28,37\text{ tấn C/ha}$ ($104,02\text{ tấn CO}_2\text{/ha}$), khẳng định giá trị dịch vụ hệ sinh thái to lớn của loài cây này.

Nghiên cứu cung cấp cơ sở khoa học và thực tiễn vững chắc cho công tác quản lý rừng bền vững, quy hoạch kinh doanh gỗ lớn kết hợp bán tín chỉ carbon, góp phần thực hiện mục tiêu giảm phát thải ròng (Net Zero) của quốc gia.