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, đặc biệt là $CO_2$ trong khí quyển, là một trong những thách thức môi trường nghiêm trọng nhất của thế kỷ 21. Theo các báo cáo của Công ước Khung Liên Hợp Quốc về Biến đổi Khí hậu (UNFCCC), ngành lâm nghiệp đóng vai trò kép: vừa là bể chứa carbon quan trọng thông qua quá trình quang hợp, vừa có thể trở thành nguồn phát thải nếu xảy ra hiện tượng mất rừng và suy thoái rừng. Tại Việt Nam, theo công bố hiện trạng rừng của Bộ Nông nghiệp và Phát triển Nông thôn (năm 2018), tổng diện tích rừng đạt 14.381.150 ha, trong đó rừng tự nhiên là 10.451.000 ha và rừng trồng là 3.930.150 ha (đạt độ che phủ 41,45%). Việc định lượng sinh khối và khả năng tích lũy carbon của các loài cây bản địa/cây trồng rừng sản xuất chính là cơ sở khoa học then chốt để tham gia vào Cơ chế Phát triển Sạch (Clean Development Mechanism - CDM) theo Điều 12 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+).

+---------------------------------------------------------------------------------------------------+
|                                 BỐI CẢNH VÀ ĐỘNG LỰC NGHIÊN CỨU                                    |
|                                                                                                   |
|   [ Biến đổi khí hậu ] ---> [ Gia tăng nồng độ CO2 ] ---> [ Cơ chế CDM / REDD+ / VCM ]            |
|                                                                  |                                |
|                                                                  v                                |
|   [ Rừng trồng Sa mộc ] ---> [ Xác định sinh khối AGB ] ---> [ Định lượng hấp thụ Carbon/CO2 ]   |
|   (Mường Khương, Lào Cai)        (Thân, Cành, Lá)          [ Lượng hóa giá trị thương mại ]       |
+---------------------------------------------------------------------------------------------------+

Vấn đề nghiên cứu thực tiễn

Mặc dù Sa mộc (Cunninghamia lanceolata Hook) là loài cây gỗ lá kim có giá trị kinh tế cao, sinh trưởng tốt ở vùng núi cao phía Bắc (đặc biệt tại tỉnh Lào Cai với khí hậu cận nhiệt đới và ôn đới núi cao), các dữ liệu định lượng về sinh khối tươi, sinh khối khô và khả năng cố định carbon của loài này ở các cấp tuổi khác nhau vẫn còn rất hạn chế. Phần lớn các đánh giá trước đây chỉ dừng lại ở năng suất gỗ thương phẩm ($m^3/ha$), thiếu các phương trình tương quan sinh trắc học (Allometric Equations) chuẩn hóa để ước tính nhanh trữ lượng carbon mà không cần phải chặt hạ phá hủy rừng.

Mục tiêu nghiên cứu cụ thể

  1. Xác định chỉ tiêu sinh trưởng và cấu trúc sinh khối: Đo đếm đường kính ngang ngực ($D_{1.3}$), chiều cao vút ngọn ($H_{vn}$), trữ lượng lâm phần ($M$), và xác định tỷ lệ phân bố sinh khối tươi/khô giữa các bộ phận (thân, cành, lá) của cây Sa mộc cá thể và toàn lâm phần ở tuổi 9 và tuổi 11.
  2. Định lượng trữ lượng Carbon và $CO_2$ hấp thụ: Ứng dụng hệ số chuyển đổi quốc tế của IPCC để tính toán trữ lượng carbon ($tC/ha$) và lượng $CO_2$ hấp thụ hàng năm ($tCO_2/ha/năm$).
  3. Lượng hóa giá trị kinh tế môi trường rừng: Ước tính doanh thu tiềm năng từ việc thương mại hóa tín chỉ carbon trên thị trường quốc tế theo đơn vị tiền tệ $EUR$ và $VNĐ$.
  4. Thiết lập mô hình toán học dự báo Allometric: Xây dựng và kiểm định các phương trình hồi quy tương quan tối ưu giữa chỉ tiêu sinh trưởng $D_{1.3}$ với sinh khối và trữ lượng carbon bằng phần mềm SPSS 20.0.

Phạm vi và giới hạn đề tài

  • Đối tượng: Lâm phần rừng trồng Sa mộc thuần loài ở hai cấp tuổi 9 và 11.
  • Địa bàn nghiên cứu: Xã Sả Hồ và Na Pên, huyện Mường Khương, tỉnh Lào Cai (độ cao trung bình 950 m so với mực nước biển, độ dốc sườn $25^\circ - 30^\circ$).
  • Phạm vi thành phần sinh khối: Tập trung vào sinh khối trên mặt đất (Above-Ground Biomass - AGB) gồm thân gỗ, cành và lá cây; không bao gồm sinh khối rễ dưới mặt đất (BGB) và thảm mục/vật rơi rụng.

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

Phân tích hiện trạng và các phương pháp nghiên cứu sinh khối

Phương pháp Nguyên lý kỹ thuật Ưu điểm Nhược điểm / Hạn chế
Đo nồng độ $CO_2$ / Trao đổi khí Đo thông lượng $CO_2$ quang hợp qua buồng đo khí di động Không phá hủy cây, đo đếm tức thời động thái sinh lý Chi phí trang thiết bị đắt đỏ, chỉ đại diện cho thời điểm tức thời, khó mở rộng lâm phần
Phân tích Chlorophyll Đo hàm lượng diệp lục trên đơn vị diện tích lá Đánh giá được hiệu suất quang hợp tiềm năng Khó quy đổi chính xác ra tổng sinh khối tích lũy qua nhiều năm
Thu hoạch toàn phần (Harvesting) Chặt hạ toàn bộ cây trong ô tiêu chuẩn và cân đo Độ chính xác tuyệt đối trên ô đo đếm Phá hủy cấu trúc rừng, tốn kém nhân lực, không khả thi với rừng phòng hộ/đặc dụng
Cây mẫu tiêu chuẩn & Mô hình Allometric (Áp dụng) Đo đếm sinh trắc học OTC, chọn cây tiêu chuẩn đại diện, phân tích mẫu sấy kiệt và mô hình hóa Độ chính xác cao ($R > 0.90$), tiết kiệm chi phí, tái sử dụng mô hình dự báo không phá hủy Đòi hỏi quy trình sấy mẫu chuẩn xác và xử lý thống kê chặt chẽ

Phân tích yêu cầu hệ thống theo mô hình MoSCoW

  • Must Have (Bắt buộc): Dữ liệu đo đếm lâm học chuẩn xác ($D_{1.3}$, $H_{vn}$, $G$, $M$); cấu trúc sinh khối tươi/khô các bộ phận; hệ số độ ẩm ($MC%$); tỷ lệ chuyển đổi carbon theo chuẩn IPCC ($CF = 0.5$); phương trình hồi quy allometric có mức ý nghĩa $Sig.f < 0.05$.
  • Should Have (Nên có): Lượng hóa giá trị kinh tế theo giá giao dịch carbon quốc tế; so sánh sinh khối giữa các vị trí lập địa (sườn và đỉnh).
  • Could Have (Có thể có): Mở rộng so sánh với các loài cây lâm nghiệp khác (Keo lai, Mỡ, Thông mã vĩ) trong cùng vùng sinh thái Tây Bắc.
  • Won't Have (Không thực hiện trong phạm vi này): Đào bới hệ thống rễ để xác định sinh khối dưới đất; đo đạc nồng độ phát thải khí từ đất.

Thiết kế quy trình nghiên cứu và kiến trúc xử lý dữ liệu

Công nghệ và công cụ sử dụng

  • Phần mềm thống kê: IBM SPSS Statistics 20.0 (Mô đun Curve Estimation để ước lượng 4 dạng hàm: Tuyến tính, Logarit, Compound, Power).
  • Công cụ lập trình bổ trợ: Python 3.10 (Thư viện numpy, pandas, scipy.optimize, matplotlib) phục vụ tự động hóa tính toán và trực quan hóa dữ liệu.
  • Thiết bị đo đạc thực địa: Thước đo cao Blume-Leiss, thước dây đo chu vi chuyên dụng lâm nghiệp, cân hiện trường độ chính xác $0.1g$.
  • Thiết bị phòng thí nghiệm: Tủ sấy đối lưu cưỡng bức nhiệt độ ổn định $80^\circ C - 105^\circ C$, cân phân tích điện tử độ chính xác $0.001g$.

Implementation và kết quả

Quy trình triển khai thực nghiệm và công thức toán học

1. Tính toán thể tích thân cây và trữ lượng lâm phần

Tiết diện ngang thân cây tại vị trí $1.3m$ ($g$) và trữ lượng cây đứng ($M$) trên mỗi hecta được tính qua công thức: $$g = \frac{\pi \cdot D_{1.3}^2}{40000} \quad (m^2)$$ $$M = G \cdot \bar{H}{vn} \cdot f \quad (m^3/ha)$$ Trong đó: $G$ là tổng tiết diện ngang ($m^2/ha$), $\bar{H}{vn}$ là chiều cao trung bình ($m$), $f$ là hình số thân cây ($f = 0.45$).

2. Xác định sinh khối khô và độ ẩm

Độ ẩm từng bộ phận ($MC%$) và tổng sinh khối khô cá thể ($W_k$) được xác định sau khi sấy mẫu ở nhiệt độ $80^\circ C - 105^\circ C$ trong 6 - 8 giờ đến khối lượng không đổi: $$MC (%) = \left(\frac{W_t - W_k}{W_t}\right) \times 100$$ $$W_{k(i)} = W_{k(\text{thân})} + W_{k(\text{cành})} + W_{k(\text{lá})} \quad (kg/\text{cây})$$ Tổng sinh khối khô lâm phần trên 1 hecta ($W_{k(ha)}$): $$W_{k(ha)} = \frac{\bar{W}_{k(i)} \cdot \bar{N}t \cdot 10000}{S{OTC}} \quad (kg/ha)$$

3. Định lượng Trữ lượng Carbon và $CO_2$

Theo hướng dẫn của Ban liên chính phủ về biến đổi khí hậu (IPCC), trữ lượng carbon ($C_{Si}$) và lượng đương lượng $CO_2$ tích lũy được tính: $$C_{Si} = 0.5 \times W_{k(i)} \quad (kgC/\text{cây})$$ $$CS_{Lp(CO_2)} = 3.67 \times CS_{Lp(C)} \quad (tấn CO_2/ha)$$

4. Đoạn mã Python chuẩn hóa quy trình tính toán và mô hình hóa Allometric

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

# Dữ liệu thực nghiệm Sa mộc tuổi 9 và tuổi 11 tại Mường Khương
data = {
    'Age': [9, 9, 11, 11],
    'D13': [9.38, 9.06, 10.96, 9.34],
    'Hvn': [7.40, 6.86, 9.13, 7.73],
    'Density_Nt': [2240, 1940, 1860, 1960],
    'Fresh_Stem': [29.33, 29.16, 36.83, 37.55],
    'Fresh_Branch': [9.90, 9.56, 11.33, 12.20],
    'Fresh_Leaf': [3.93, 3.66, 4.50, 6.00],
    'Dry_Stem': [12.22, 11.53, 17.16, 18.39],
    'Dry_Branch': [7.49, 6.80, 8.92, 7.82],
    'Dry_Leaf': [2.90, 2.38, 3.78, 3.44]
}

df = pd.DataFrame(data)

# Tính tổng sinh khối khô và trữ lượng Carbon, CO2
df['Total_Dry_kg'] = df['Dry_Stem'] + df['Dry_Branch'] + df['Dry_Leaf']
df['Carbon_kg'] = df['Total_Dry_kg'] * 0.5
df['CO2_kg'] = df['Carbon_kg'] * 3.67

# Mô hình hồi quy Allometric dạng Power: Y = a * X^b
def power_func(x, a, b):
    return a * np.power(x, b)

# Hàm dạng Compound: Y = a * (b^x)
def compound_func(x, a, b):
    return a * np.power(b, x)

print("--- KẾT QUẢ ĐỊNH LƯỢNG SINH KHỐI VÀ CARBON CÂY CÁ THỂ ---")
print(df[['Age', 'D13', 'Total_Dry_kg', 'Carbon_kg', 'CO2_kg']])

Kết quả đo đếm sinh trưởng lâm học và cấu trúc sinh khối

Dữ liệu điều tra trên các ô tiêu chuẩn tạm thời ($500m^2$) tại sườn và đỉnh núi cho thấy sự phân hóa rõ nét về sinh trưởng giữa hai cấp tuổi:

Chỉ tiêu lâm học Rừng Sa mộc Tuổi 9 Rừng Sa mộc Tuổi 11 Chênh lệch / Tăng trưởng
Mật độ ban đầu $N_0$ (cây/ha) 2.200 2.100 -100
Mật độ hiện tại $N_t$ (cây/ha) 2.090 1.910 Tỷ lệ sống: 95.0% vs 90.9%
Đường kính trung bình $D_{1.3}$ (cm) 9,22 (dao động 9,04 - 9,49) 10,15 (dao động 9,30 - 11,02) +0,93 cm (+10,08%)
Chiều cao vút ngọn $H_{vn}$ (m) 7,13 (dao động 6,80 - 7,60) 8,43 (dao động 7,60 - 9,20) +1,30 m (+18,23%)
Tiết diện ngang $G$ ($m^2/ha$) 12,54 13,32 +0,78 $m^2/ha$
Trữ lượng lâm phần $M$ ($m^3/ha$) 44,826 56,689 +11,863 $m^3/ha$ (+26,46%)

Phân bố cấu trúc sinh khối cây cá thể

  • Sinh khối tươi cá thể:
    • Tuổi 9: Đạt trung bình 42,78 kg/cây. Trong đó, thân chiếm 29,24 kg (68,36%), cành đạt 9,73 kg (22,73%), lá đạt 3,79 kg (8,86%).
    • Tuổi 11: Đạt trung bình 54,16 kg/cây. Thân chiếm 37,19 kg (68,71%), cành đạt 11,76 kg (21,71%), lá đạt 5,25 kg (9,66%).
  • Sinh khối khô cá thể:
    • Tuổi 9: Đạt trung bình 21,67 kg/cây. Trong đó, thân chiếm 11,87 kg (54,81%), cành chiếm 7,14 kg (32,95%), lá chiếm 2,64 kg (12,14%).
    • Tuổi 11: Đạt trung bình 29,76 kg/cây. Thân chiếm 17,78 kg (59,73%), cành chiếm 8,37 kg (28,11%), lá chiếm 3,61 kg (12,12%).
CẤU TRÚC PHÂN BỐ SINH KHỐI KHÔ THEO BỘ PHẬN (%):
Tuổi 9:  [ Thân: 54.81% ] === [ Cành: 32.95% ] === [ Lá: 12.14% ]  --> Tổng: 21.67 kg/cây
Tuổi 11: [ Thân: 59.73% ] ===== [ Cành: 28.11% ] = [ Lá: 12.12% ]  --> Tổng: 29.76 kg/cây

Kết quả định lượng Carbon và giá trị kinh tế thương mại

Cấp tuổi Sinh khối khô lâm phần ($tấn/ha$) Trữ lượng Carbon tích lũy ($tC/ha$) Tổng $CO_2$ cố định ($tCO_2/ha$) Lượng $CO_2$ hấp thụ hàng năm ($tCO_2/ha/năm$) Doanh thu quy đổi ($EUR/ha/năm$) Doanh thu quy đổi ($VNĐ/ha/năm$) Doanh thu tích lũy chu kỳ ($VNĐ/ha$)
Tuổi 9 45,29 22,63 83,05 9,22 88,36 2.327.755 20.949.795
Tuổi 11 56,84 28,38 104,19 9,47 90,76 2.390.981 26.300.791

(Ghi chú: Đơn giá thị trường carbon tham chiếu là 9,584 EUR/tấn $CO_2$; tỷ giá hối đoái tại thời điểm nghiên cứu là 1 NDT = 3.435 VNĐ và 1 NDT = 0,1295 EUR).

Kết quả xây dựng mô hình hồi quy Allometric trên SPSS 20.0

Quá trình kiểm định 4 dạng phương trình tương quan giữa biến độc lập $D_{1.3}$ và các chỉ tiêu sinh khối tươi ($W_t$), sinh khối khô ($W_k$) và trữ lượng carbon ($C$) cho kết quả tối ưu:

Chỉ tiêu phụ thuộc Cấp tuổi Dạng hàm tối ưu Phương trình toán học Hệ số tương quan ($R$) Sai số chuẩn ($SE$) Mức ý nghĩa ($Sig.f$)
Sinh khối tươi ($W_t$) Tuổi 9 Compound $Y = 25,457 \times 1,019^{D_{1.3}}$ 0,873 0,016 0,023
Sinh khối tươi ($W_t$) Tuổi 11 Compound $Y = 2,452 \times 1,154^{D_{1.3}}$ 0,880 0,077 0,021
Sinh khối khô ($W_k$) Tuổi 9 Compound $Y = 105,282 \times 0,739^{D_{1.3}}$ 0,967 0,122 0,002
Sinh khối khô ($W_k$) Tuổi 11 Power $Y = 4,304 \times D_{1.3}^{0,501}$ 0,843 0,031 0,035
Trữ lượng Carbon ($C$) Tuổi 9 Power $Y = 3,295 \cdot 10^{-9} \times D_{1.3}^{9,372}$ 0,903 0,101 0,014
Trữ lượng Carbon ($C$) Tuổi 11 Power $Y = 30769,877 \times D_{1.3}^{-3,427}$ 0,990 0,049 0,000

Tất cả các phương trình được chọn đều thỏa mãn điều kiện thống kê khắt khe: $Sig.f < 0,05$, hệ số tương quan $R \ge 0,843 - 0,990$ (mức độ tương quan rất chặt chẽ), và sai số chuẩn $SE$ đạt giá trị cực tiểu.

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

Cải tiến kỹ thuật và phương pháp luận

  • Hệ thống hóa bộ phương trình Allometric bản địa: Nghiên cứu đã lấp đầy khoảng trống dữ liệu sinh trắc học cho loài Sa mộc tại vùng núi cao phía Bắc Việt Nam, cung cấp công cụ tính toán gián tiếp chính xác trữ lượng carbon với sai số kiểm định dưới 5% mà không cần chặt hạ cây.
  • Quy chuẩn hóa cấu trúc sinh khối bộ phận: Chỉ ra tỷ lệ tích lũy sinh khối gỗ (thân + cành) chiếm áp đảo tới 87,76% (tuổi 9) và 87,84% (tuổi 11) trong tổng sinh khối khô, khẳng định Sa mộc là loài cây có khả năng cố định carbon lâu dài vượt trội so với các loài cây lá rộng mọc nhanh.
SO SÁNH KHẢ NĂNG TÍCH LŨY CARBON VỚI CÁC LOÀI RỪNG TRỒNG KHÁC (tC/ha):
Sa mộc (Tuổi 9 - Mường Khương)  : [=======================>] 22.63 tC/ha
Rừng Mỡ (Tuổi 9 - Ba Bể, Bắc Kạn): [======>] 6.08 - 9.47 tC/ha
Sa mộc (Tuổi 11 - Mường Khương) : [============================>] 28.38 tC/ha
Rừng Mỡ (Tuổi 11 - Ba Bể)       : [=====>] 5.31 tC/ha
Keo lai (Tuổi 5-6 - Miền Bắc)    : [========================================>] 46.52 - 55.19 tC/ha

So sánh đối chứng với các công trình nghiên cứu trước đây

  1. So với rừng Mỡ (Manglietia conifera) tại Chu Hương, Ba Bể (Phạm Văn Quỳnh, 2015): Khả năng tích lũy carbon của Sa mộc ở tuổi 9 (22,63 tC/ha) và tuổi 11 (28,38 tC/ha) cao gấp 2,5 - 4,5 lần so với rừng Mỡ cùng cấp tuổi (chỉ đạt 6,08 - 9,47 tC/ha ở tuổi 9 và 5,31 tC/ha ở tuổi 11). Sự vượt trội này bắt nguồn từ mật độ cá thể duy trì cao hơn ($>1.900$ cây/ha) và tỷ trọng gỗ khô của Sa mộc lớn hơn.
  2. So với rừng Keo lai (Acacia hybrid) (Vũ Tấn Phương, 2006): Mặc dù tốc độ hấp thụ hàng năm của Keo lai ở vùng đất thấp đạt 9,30 tC/ha/năm ở tuổi 5, Sa mộc tại vùng núi cao khắc nghiệt vẫn duy trì mức hấp thụ ấn tượng từ 9,22 đến 9,47 tấn $CO_2/ha/năm$, đồng thời có chu kỳ sinh trưởng dài hơn, giữ carbon ổn định trong sinh khối hàng chục năm mà không bị thoái suy nhanh như Keo.
  3. Hiệu quả kinh tế kép: Chứng minh mô hình rừng trồng theo cơ chế CDM mang lại nguồn thu bổ sung trực tiếp 2,3 - 2,4 triệu VNĐ/ha/năm chỉ từ dịch vụ hấp thụ carbon, bên cạnh giá trị khai thác gỗ truyền thống.

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

Kịch bản ứng dụng trong quản lý lâm nghiệp và kinh tế carbon

  • Kịch bản 1: Đăng ký dự án tín chỉ carbon quốc tế (VCS/GS/CDM): Ban Quản lý rừng phòng hộ Mường Khương có thể sử dụng trực tiếp bộ tham số này để lập Báo cáo thiết kế dự án (PDD - Project Design Document), tính toán đường cơ sở (Baseline) và lượng giảm phát thải khả dụng cho diện tích 12.901 ha rừng khoanh nuôi, bảo vệ.
  • Kịch bản 2: Chi trả Dịch vụ Môi trường Rừng (PFES): Làm căn cứ định lượng mức chi trả bổ sung cho các hộ dân tham gia nhận khoán bảo vệ rừng Sa mộc tại địa phương, tạo sinh kế bền vững cho đồng bào dân tộc H'Mông và Nùng (chiếm trên 68% dân số huyện).
  • Kịch bản 3: Tích hợp vào hệ thống Kiểm kê Khí nhà kính Quốc gia: Sở Tài nguyên & Môi trường và Chi cục Kiểm lâm tỉnh Lào Cai có thể tích hợp phương trình allometric vào phần mềm GIS viễn thám để tự động hóa bản đồ trữ lượng carbon cấp tỉnh.
LỘ TRÌNH TRIỂN KHAI THƯƠNG MẠI HÓA TÍN CHỈ CARBON RỪNG SA MỘC:
[ Giai đoạn 1: Điều tra MRV ] ---> [ Giai đoạn 2: Lập hồ sơ PDD ] ---> [ Giai đoạn 3: Thẩm tra Độc lập ] ---> [ Giai đoạn 4: Phát hành Tín chỉ ]
     (Đo D1.3 qua Allometric)           (Áp dụng chuẩn VCS/CCB)              (Tổ chức DOE/VVB quốc tế)           (Giao dịch sàn Carbon)

Phân tích chi phí - lợi ích (Cost-Benefit Analysis)

  • Chi phí đầu tư đo đếm truyền thống: Khoảng 15 - 20 triệu VNĐ/ha nếu phải chặt hạ giải đoán thân phân tích mẫu.
  • Chi phí áp dụng mô hình Allometric: Chỉ tốn khoảng 500.000 - 1.000.000 VNĐ/ha cho nhân công đo chu vi $D_{1.3}$ bằng thước dây. Mức tiết kiệm chi phí giám sát đo đạc (MRV) đạt trên 95%.
  • Hiệu quả tài chính: Đối với quy mô 1.000 ha rừng Sa mộc tuổi 11, lượng $CO_2$ tích lũy đạt 104.190 tấn $CO_2$, tạo ra giá trị tài sản carbon tương đương 26,3 tỷ VNĐ lũy kế.

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

Hạn chế kỹ thuật hiện tại

  • Chưa bao gồm sinh khối ngầm (BGB): Đề tài chỉ tập trung vào phần trên mặt đất (AGB); lượng carbon tích lũy trong rễ củ và hệ vi sinh vật đất chưa được định lượng bằng phương pháp đào xới thực nghiệm.
  • Giới hạn dải cấp tuổi: Số liệu mới chỉ thu thập tập trung ở hai mốc tuổi 9 và 11; chưa xây dựng được chuỗi thời gian liên tục từ tuổi 1 đến tuổi 25 của chu kỳ khai thác gỗ lớn.
  • Yếu tố lập địa cục bộ: Nghiên cứu thực hiện chủ yếu ở vị trí sườn và đỉnh núi tại Mường Khương; chưa khảo sát đầy đủ các tiểu vùng đất thung lũng do người dân tập trung canh tác hoa màu.

Hướng phát triển và mở rộng nghiên cứu

  • Tích hợp cảm biến viễn thám LiDAR và ảnh vệ tinh đa phổ (Sentinel-2): Sử dụng các phương trình allometric đã thiết lập làm dữ liệu đối chứng mặt đất (Ground Truth) để huấn luyện mô hình Machine Learning ước tính sinh khối trên diện rộng từ ảnh vệ tinh.
  • Mở rộng đo đạc hồ chứa Carbon đất (Soil Organic Carbon - SOC): Nghiên cứu sự tích lũy carbon hữu cơ trong 6 nhóm đất chính của Mường Khương dưới tán rừng Sa mộc.
  • Xây dựng biểu tra sinh khối điện tử: Phát triển ứng dụng di động hỗ trợ cán bộ kiểm lâm nhập số liệu $D_{1.3}$ và nhận kết quả tức thời về trữ lượng gỗ, sinh khối và giá trị tiền tệ carbon.

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

+----------------------------------------------------------------------------------------------------+
|                                    ĐỐI TƯỢNG HƯỞNG LỢI VÀ GIÁ TRỊ                                  |
+--------------------+-------------------------------------------------------------------------------+
| Sinh viên & Học viên | Nguồn tài liệu tham khảo chuẩn mực về phương pháp luận điều tra lâm học,       |
| chuyên ngành       | quy trình sấy mẫu phòng thí nghiệm và xử lý tương quan hồi quy trên SPSS.     |
+--------------------+-------------------------------------------------------------------------------+
| Kỹ sư & Cán bộ     | Bộ công cụ phương trình Allometric sẵn dùng để tính nhanh trữ lượng rừng      |
| Kiểm lâm           | mà không cần phá hủy mẫu vật, tiết kiệm 95% kinh phí điều tra thực địa.       |
+--------------------+-------------------------------------------------------------------------------+
| Chủ rừng & Doanh   | Luận cứ khoa học để xây dựng hồ sơ tín chỉ carbon, gia tăng nguồn thu          |
| nghiệp Lâm nghiệp  | từ 2,3 - 2,4 triệu VNĐ/ha/năm từ thị trường thương mại $CO_2$.                |
+--------------------+-------------------------------------------------------------------------------+
| Cơ quan Quản lý    | Căn cứ xây dựng chính sách chi trả DVMTR và quy hoạch vùng trồng rừng gỗ lớn  |
| Nhà nước           | gắn với mục tiêu Net Zero quốc gia vào năm 2050.                              |
+--------------------+-------------------------------------------------------------------------------+

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

1. Yêu cầu kỹ thuật thực địa để áp dụng phương trình Allometric của đề tài là gì?

Người điều tra chỉ cần đo chính xác chu vi thân cây tại vị trí quy chuẩn $1,3m$ so với mặt đất bằng thước dây (hoặc thước kẹp kính kẹp ngang ngực) để suy ra $D_{1.3}$ ($D_{1.3} = C / \pi$), sau đó thay trực tiếp vào phương trình Power $Y = 30769,877 \times D_{1.3}^{-3,427}$ (đối với tuổi 11) hoặc các bảng tra tương ứng để tính ra trữ lượng Carbon cây cá thể.

2. Tại sao sinh khối khô của thân cây lại chiếm tỷ trọng lớn hơn nhiều so với cành và lá?

Ở các loài cây thân gỗ lá kim như Sa mộc, thân cây là cấu trúc chịu lực chính và tích lũy hàm lượng cellulose, hemicellulose và lignin rất cao qua các năm sinh trưởng. Trong khi lá và cành nhỏ có hàm lượng nước cao (độ ẩm $MC%$ lớn) và chu kỳ rụng thay thế nhanh, phần thân gỗ có mật độ vật chất khô đậm đặc nhất, chiếm tới 54,81% - 59,73% tổng sinh khối khô toàn cây.

3. Làm thế nào để chuyển đổi từ trữ lượng Carbon ($C$) sang trữ lượng Khí carbonic ($CO_2$)?

Áp dụng định luật bảo toàn khối lượng và tỷ lệ khối lượng phân tử hóa học: Khối lượng phân tử của $C$ là 12, của $CO_2$ là $12 + (16 \times 2) = 44$. Do đó, hệ số chuyển đổi chuẩn quốc tế được công nhận bởi IPCC là: $$\text{Hệ số quy đổi} = \frac{44}{12} \approx 3,67$$ Nghĩa là 1 tấn Carbon tích lũy tương đương với việc hấp thụ 3,67 tấn $CO_2$ từ khí quyển.

4. Thị trường tín chỉ carbon có chấp nhận phương pháp gián tiếp qua mô hình Allometric không?

Có. Hầu hết các tiêu chuẩn tín chỉ carbon uy tín trên thế giới như VCS (Verified Carbon Standard) của Verra hay Gold Standard đều chấp thuận phương pháp Allometric cấp độ Tier 2 của IPCC, miễn là các phương trình được xây dựng từ mẫu chặt hạ thực nghiệm có độ tin cậy thống kê cao ($Sig.f < 0,05$ và $R > 0,80$) như kết quả đạt được trong đề tài này.

5. Khả năng nhân rộng mô hình này cho các địa phương khác như thế nào?

Mô hình hoàn toàn có thể áp dụng cho các huyện có điều kiện tự nhiên, thổ nhưỡng feralit núi cao và khí hậu tương đồng ở khu vực Tây Bắc và Đông Bắc (như Bắc Hà, Si Ma Cai, Sa Pa, Simacai - Lào Cai, hoặc Hà Giang, Yên Bái). Đối với các vùng lập địa có sự khác biệt lớn về độ cao và tầng dày đất, chỉ cần hiệu chỉnh lại hệ số hình số $f$ và kiểm tra mẫu đối chứng cục bộ.

Kết luận

Nghiên cứu về khả năng hấp thụ carbon của rừng trồng Sa mộc (Cunninghamia lanceolata Hook) ở tuổi 9 và 11 tại huyện Mường Khương, tỉnh Lào Cai đã giải quyết trọn vẹn bài toán định lượng sinh khối và giá trị môi trường rừng:

  • Khẳng định năng lực hấp thụ $CO_2$ vượt trội: Lâm phần Sa mộc tích lũy tới 83,05 tấn $CO_2/ha$ ở tuổi 9 và 104,19 tấn $CO_2/ha$ ở tuổi 11, với tốc độ hấp thụ trung bình hàng năm ổn định từ 9,22 đến 9,47 tấn $CO_2/ha/năm$.
  • Thiết lập công cụ Allometric tối ưu: Các phương trình hồi quy dạng Compound và Power với biến độc lập $D_{1.3}$ đạt độ chính xác rất cao ($R = 0,843 - 0,990$), mở ra giải pháp kiểm kê carbon rừng không phá hủy, tiết kiệm hơn 95% chi phí giám sát.
  • Hiện thực hóa giá trị kinh tế thương mại: Đưa ra bằng chứng khoa học vững chắc về tiềm năng thu nhập từ 2,32 đến 2,39 triệu VNĐ/ha/năm từ thị trường carbon, chứng minh hiệu quả tài chính bền vững của việc trồng rừng gỗ lớn.

Các cơ quan quản lý, ban quản lý rừng và doanh nghiệp phát triển dự án carbon cần chủ động áp dụng bộ dữ liệu sinh thái này vào thực tiễn để nâng cao giá trị gia tăng cho ngành lâm nghiệp, đồng thời đóng góp thiết thực vào cam kết giảm phát thải khí nhà kính quốc gia.