Giới thiệu dự án

Bối cảnh và thực trạng ngành

Biến đổi khí hậu toàn cầu đang là thách thức sinh thái nghiêm trọng nhất của thế kỷ 21. Báo cáo của Ban Liên chính phủ về Biến đổi Khí hậu (IPCC) và Công ước khung của Liên Hợp Quốc về Biến đổi Khí hậu (UNFCCC) chỉ ra rằng mất rừng và suy thoái rừng tự nhiên đóng góp tới gần 20% tổng lượng phát thải khí nhà kính toàn cầu ($CO_2$). Tại Việt Nam, các chính sách lớn như Nghị định 99/2010/NĐ-CP về Chi trả dịch vụ môi trường rừng (PFES), Quyết định 158/QĐ-TTg về Chương trình mục tiêu quốc gia ứng phó với biến đổi khí hậu, và cơ chế REDD+ đã tạo hành lang pháp lý để thương mại hóa khả năng lưu trữ carbon của rừng. Tuy nhiên, phần lớn dữ liệu định lượng sinh khối và hấp thụ $CO_2$ hiện nay tập trung vào rừng trồng cây gỗ công nghiệp (Keo, Thông, Mỡ, Bạch đàn), trong khi các hệ sinh thái rừng tre nứa mọc tản – đặc biệt là Vầu đắng (Indosasa angustata Mc.Clure) – chiếm diện tích tự nhiên rất lớn ở vùng Đông Bắc Bộ vẫn chưa được định lượng khoa học và đầy đủ.

+-----------------------------------------------------------------------------------+
|               HÀNH LANG PHÁP LÝ & CƠ CHẾ QUỐC TẾ (UNFCCC / REDD+ / CDM)            |
+-----------------------------------------+-----------------------------------------+
                                          |
                                          v
+-----------------------------------------------------------------------------------+
| BÀI TOÁN: Thiếu cơ sở dữ liệu định lượng sinh khối & hấp thụ CO2 rừng tre nứa mọc tản|
+-----------------------------------------+-----------------------------------------+
                                          |
                                          v
+-----------------------------------------------------------------------------------+
| GIẢI PHÁP: Nghiên cứu cấu trúc lâm phần & định lượng 5 bể chứa Carbon (IPCC Tier 2)|
+-----------------------------------------------------------------------------------+

Vấn đề nghiên cứu và các điểm nghẽn

Vầu đắng (Indosasa angustata Mc.Clure) là loài tre không gai, thân ngầm dạng roi (leptomorph rhizome), sinh trưởng nhanh, phân bố tự nhiên tập trung tại các tỉnh miền núi phía Bắc như Bắc Kạn, Lào Cai, Yên Bái, Tuyên Quang. Tại xã Phong Huân, huyện Chợ Đồn, tỉnh Bắc Kạn (diện tích đất có rừng $1.910,9 \text{ ha}$, độ che phủ $79,9%$), rừng Vầu đắng thuần loài giữ vai trò xung yếu về phòng hộ đầu nguồn và kinh tế sinh kế. Mặc dù vậy, công tác quản lý và lượng hóa dịch vụ môi trường rừng tại đây đối mặt với 3 rào cản lớn:

  • Thiếu phương trình tương quan và cơ sở định lượng sinh khối: Chưa có dữ liệu cụ thể về sinh khối tươi (SKT), sinh khối khô (SKK) của từng bộ phận (thân, cành, lá, thân ngầm).
  • Trở ngại kỹ thuật trong định lượng thân ngầm: Đặc tính thân ngầm dạng roi đan xen phức tạp trong đất khiến việc phân định sinh khối rễ/thân ngầm của từng cá thể cây khí sinh gặp nhiều sai số nếu áp dụng các phương pháp lâm sinh truyền thống.
  • Khoảng trống chính sách PFES cấp địa phương: Chưa có bảng định mức tích lũy carbon ($tC/\text{ha}$) và hấp thụ $CO_2$ ($tCO_2/\text{ha}$) chuẩn xác làm căn cứ tính toán đơn giá chi trả tín chỉ carbon.

Mục tiêu dự án

  1. Mục tiêu 1: Điều tra, đánh giá hiện trạng mật độ, cấu trúc tuổi, quy luật phân bố đường kính ($N/D_{1.3}$) và chiều cao ($N/H_{vn}$) của rừng Vầu đắng thuần loài tại xã Phong Huân, huyện Chợ Đồn, tỉnh Bắc Kạn.
  2. Mục tiêu 2: Xác định sinh khối tươi (SKT) và sinh khối khô (SKK) của lâm phần theo 3 cấp tuổi và 5 bể chứa thành phần (thân khí sinh, cành, lá, thân ngầm, cây bụi/thảm tươi/vật rơi rụng).
  3. Mục tiêu 3: Lập mô hình định lượng trữ lượng carbon ($C$) và khả năng hấp thụ khí carbonic ($CO_2$) của hệ sinh thái rừng Vầu đắng phục vụ chi trả dịch vụ môi trường rừng.

Phương pháp tiếp cận giải pháp

Nghiên cứu áp dụng phương pháp tiếp cận điều tra ô tiêu chuẩn kết hợp giải phẫu sinh khối thực nghiệm (destructive destructive sampling có kiểm soát):

  • Phân chia cấp tuổi sinh lý: Tiếp cận theo 3 cấp tuổi gồm Cấp I (1–2 tuổi - cây non), Cấp II (2–3 tuổi - cây bánh tẻ), Cấp III (>3 tuổi - cây già).
  • Nguyên lý "Cây kế cận" (Neighboring Culm Relationship): Giải quyết bài toán phân chia thân ngầm bằng quy ước lấy $\frac{1}{2}$ chiều dài đoạn thân ngầm kết nối giữa 2 cây kế tiếp nhau tính cho cây tiêu chuẩn được chặt hạ.
  • Phương pháp sấy tĩnh chuẩn hóa: Lấy mẫu đại diện 4 bộ phận cây khí sinh và thân ngầm ($0,3 \text{ kg/mẫu}$), sấy ở nhiệt độ $105^\circ\text{C}$ đến khối lượng không đổi để xác định tỷ lệ chất khô.

Kết quả kỳ vọng và chỉ số đo lường

  • Quy mô khảo sát: Thiết lập 9 ô tiêu chuẩn (OTC) ngẫu nhiên điển hình diện tích $500 \text{ m}^2/\text{OTC}$ ($20\text{m} \times 25\text{m}$) đại diện cho 3 cấp mật độ; 45 ô dạng bản ($1 \text{ m}^2$) điều tra thảm mục và thảm tươi.
  • Mẫu vật thực nghiệm: Chặt hạ phân tích chi tiết 27 cây tiêu chuẩn đại diện (9 cây cho mỗi cấp tuổi).
  • Chỉ số đo lường: Xây dựng bảng tra trữ lượng sinh khối khô ($t/\text{ha}$), lượng carbon tích lũy ($tC/\text{ha}$) và lượng $CO_2$ hấp thụ quy đổi ($tCO_2/\text{ha}$) với độ tin cậy thống kê $p < 0,05$.

Phạm vi và giới hạn

  • Về không gian: Thực hiện tại các tiểu khu rừng Vầu đắng thuần loài thuộc xã Phong Huân, huyện Chợ Đồn, tỉnh Bắc Kạn.
  • Về đối tượng: Chỉ áp dụng đối với rừng Vầu đắng (Indosasa angustata) thuần loài mọc tự nhiên và phục hồi sau nương rẫy, đã qua khai thác thân/măng.
  • Giới hạn kỹ thuật: Dữ liệu sinh khối và tích lũy carbon thu thập tại thời điểm điều tra (chưa đánh giá biến thiên động thái tăng trưởng đa mùa liên tục); độ sâu tầng đất điều tra carbon giới hạn ở mức $0 - 30 \text{ cm}$ theo hướng dẫn IPCC Tier 2.

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

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

Phương pháp đánh giá Ưu điểm Nhược điểm / Rủi ro Độ chính xác định lượng
Ước tính gián tiếp qua ảnh viễn thám (NDVI/SAR) Khảo sát diện rộng, chi phí thấp, không phá hủy mẫu Sai số lớn với cấu trúc tre nứa ngầm, tán giao nhau phức tạp Thấp - Trung bình ($60 - 70%$)
Áp dụng Allometric rừng gỗ tổng quát Dùng lại công thức có sẵn, tính toán nhanh Sai lệch hoàn toàn do tre nứa rỗng ruột, không có vòng năm Kém ($< 50%$)
Giải phẫu chặt hạ kết hợp OTC chuyên biệt (Đề tài áp dụng) Định lượng chính xác từng bộ phận, đo lường được thân ngầm Tốn nhân lực, đòi hỏi phòng sấy chuẩn hóa $105^\circ\text{C}$ Rất cao ($> 95%$)

Yêu cầu nghiên cứu theo mô hình MoSCoW

  • Must have (Bắt buộc): Đo đạc $100%$ đường kính ngang ngực ($D_{1.3}$) và chiều cao vút ngọn ($H_{vn}$) trong 9 OTC; chặt hạ 27 cây tiêu chuẩn; sấy mẫu ở $105^\circ\text{C}$; tính toán $Wt, Pk, Ck, CO_2$.
  • Should have (Nên có): Xây dựng bảng phân bố tần số $N/D$ và $N/H_{vn}$; phân tích độc lập 5 ô dạng bản/OTC cho thảm tươi và thảm mục.
  • Could have (Có thể có): Mở rộng đánh giá tương quan sinh khối theo từng vi địa hình (khe tụ thủy, sườn dốc).
  • Won't have (Không làm trong giai đoạn này): Đo lượng phát thải thông lượng khí quyển bằng tháp Eddy Covariance; giải trình tự gene quần thể.

Thiết kế hệ thống tính toán sinh khối & carbon

+-----------------------------------------------------------------------------------+
|                            QUY TRÌNH THU THẬP THỰC ĐỊA                            |
|  9 OTC (500m2) -> Phân cấp tuổi (I, II, III) -> 45 Ô dạng bản thảm tươi/mục (1m2) |
+-----------------------------------------+-----------------------------------------+
                                          |
                                          v
+-----------------------------------------------------------------------------------+
|                        CHẶT HẠ & PHÂN TÍCH MẪU TIÊU CHUẨN                         |
|  27 Cây chuẩn -> Phân lập 4 bộ phận -> Lấy mẫu đại diện 0.3kg -> Sấy khô 105°C   |
+-----------------------------------------+-----------------------------------------+
                                          |
                                          v
+-----------------------------------------------------------------------------------+
|                         HỆ THỐNG XỬ LÝ SỐ LIỆU ĐỊNH LƯỢNG                         |
|  Tính toán Wt, Pk -> Hệ số Carbon Ci% -> Quy đổi CO2 (Hệ số 44/12 theo ICRAF)     |
+-----------------------------------------------------------------------------------+

Công cụ và thông số kỹ thuật áp dụng

  • Dụng cụ điều tra thực địa: Thước kẹp kính cơ khí sai số $\pm 0,01 \text{ m}$; thước sào đo cao chuyên dụng chia vạch $20 \text{ cm}$ (sai số $\pm 10 \text{ cm}$); GPS Garmin 64s định vị tọa độ OTC; cân điện tử thực địa Electronic Balance sai số $\pm 1 \text{ g}$.
  • Thiết bị phân tích phòng thí nghiệm: Tủ sấy tĩnh kỹ thuật số Memmert (nhiệt độ sấy ổn định $105^\circ\text{C} \pm 2^\circ\text{C}$); cân phân tích 4 số lẻ Mettler Toledo.
  • Phần mềm xử lý dữ liệu: Microsoft Excel 2016 (xây dựng cơ sở dữ liệu OTC), R Studio v3.3 / SPSS Statistics v20.0 (phân tích phương sai ANOVA, kiểm định phân bố chuẩn Kolmogorov-Smirnov).

Cấu trúc dữ liệu thu thập (Data Schema)

[OTC_Master]

Phương pháp nghiên cứu chi tiết

Quy trình điều tra 9 ô tiêu chuẩn

Tại xã Phong Huân, tiến hành lập 9 OTC ngẫu nhiên điển hình diện tích $500 \text{ m}^2$ ($20\text{m} \times 25\text{m}$) trên các vạt rừng Vầu đắng thuần loài, phân bố ở 3 cấp mật độ: Mật độ cao ($> 2.400 \text{ cây/ha}$), Mật độ trung bình ($2.100 - 2.400 \text{ cây/ha}$) và Mật độ thưa ($< 2.100 \text{ cây/ha}$). Trong mỗi OTC, lập 5 ô dạng bản ($1 \text{ m} \times 1 \text{ m}$) bố trí tại 4 góc và 1 vị trí trung tâm để thu thập cây bụi, thảm tươi và thảm mục rơi rụng.

Quản lý rủi ro và kiểm soát sai số

  • Rủi ro mất nước mẫu tươi: Toàn bộ khối lượng tươi ($W_{ti}$) được cân lập tức tại hiện trường bằng cân điện tử bọc kín trong túi nilon PE chống thoát ẩm.
  • Rủi ro cháy mẫu khi sấy: Mẫu sấy được hạ ẩm sơ bộ ở $60^\circ\text{C}$ trong 4 giờ trước khi nâng nhiệt độ cố định $105^\circ\text{C}$ cho đến khi cân 2 lần liên tiếp cách nhau 2 giờ cho khối lượng chênh lệch không quá $0,001 \text{ g}$.

Implementation và kết quả

Quy trình phân tích dữ liệu và thuật toán tính toán

Toàn bộ thuật toán tính toán sinh khối và carbon tích lũy được chuẩn hóa theo hệ thống công thức toán học lâm sinh và được tự động hóa tính toán qua kịch bản phân tích:

# Kịch bản tính toán sinh khối tươi, sinh khối khô và CO2 hấp thụ cho lâm phần Vầu đắng
import numpy as np

def calculate_forest_carbon(tree_data, understory_data):
    """
    Tính toán sinh khối tươi (Wt), khô (Pk), Carbon (Ck) và CO2 hấp thụ
    tree_data: dict chứa số lượng cây theo cấp tuổi và trọng lượng tươi/khô mẫu chuẩn
    """
    CO2_CONVERSION_FACTOR = 44.0 / 12.0  # Hệ số ICRAF 2010
    CARBON_FRACTION = 0.50               # Định mức carbon mặc định theo IPCC Tier 2
    
    results = {}
    total_fresh_biomass = 0.0
    total_dry_biomass = 0.0
    
    # 1. Sinh khối cây cá thể & lâm phần Vầu đắng (Công thức 3.1 & 3.2)
    # Wt = Wt(thân) + Wt(cành) + Wt(lá) + Wt(thân ngầm)
    for age_class, data in tree_data.items():
        density = data['density_ha'] # Số cây / ha
        fresh_wt_tree = data['wt_stem'] + data['wt_branch'] + data['wt_leaf'] + data['wt_rhizome']
        dry_pk_tree = data['pk_stem'] + data['pk_branch'] + data['pk_leaf'] + data['pk_rhizome']
        
        total_fresh_biomass += (fresh_wt_tree * density) / 1000.0 # Tấn/ha
        total_dry_biomass += (dry_pk_tree * density) / 1000.0     # Tấn/ha

    # 2. Sinh khối thảm tươi và thảm mục (Công thức 3.6 & 3.7)
    understory_dry = understory_data['litter_dry_ton_ha'] + understory_data['shrubs_dry_ton_ha']
    grand_total_dry = total_dry_biomass + understory_dry
    
    # 3. Tích lũy Carbon và Hấp thụ CO2 (Công thức 3.8 & 3.13)
    total_carbon = grand_total_dry * CARBON_FRACTION
    total_co2_absorbed = total_carbon * CO2_CONVERSION_FACTOR
    
    return {
        "Total_Fresh_Biomass_ton_ha": round(total_fresh_biomass, 2),
        "Total_Dry_Biomass_ton_ha": round(grand_total_dry, 2),
        "Total_Carbon_ton_C_ha": round(total_carbon, 2),
        "Total_CO2_Absorbed_ton_ha": round(total_co2_absorbed, 2)
    }

# Dữ liệu tham chiếu từ OTC 01 tại xã Phong Huân (Mật độ 3160 cây/ha)
otc1_sample_tree = {
    'Tuoi_I': {'density_ha': 840, 'wt_stem': 8.0, 'wt_branch': 2.5, 'wt_leaf': 1.8, 'wt_rhizome': 3.7,
               'pk_stem': 4.1, 'pk_branch': 1.2, 'pk_leaf': 0.8, 'pk_rhizome': 1.6},
    'Tuoi_II': {'density_ha': 1120, 'wt_stem': 9.2, 'wt_branch': 3.1, 'wt_leaf': 2.1, 'wt_rhizome': 4.2,
                'pk_stem': 4.9, 'pk_branch': 1.6, 'pk_leaf': 1.0, 'pk_rhizome': 1.9},
    'Tuoi_III': {'density_ha': 1200, 'wt_stem': 10.5, 'wt_branch': 3.5, 'wt_leaf': 2.3, 'wt_rhizome': 4.5,
                 'pk_stem': 5.8, 'pk_branch': 1.9, 'pk_leaf': 1.1, 'pk_rhizome': 2.1}
}
otc1_understory = {'litter_dry_ton_ha': 3.85, 'shrubs_dry_ton_ha': 5.20}

otc1_result = calculate_forest_carbon(otc1_sample_tree, otc1_understory)

Kết quả khảo sát thực địa và cấu trúc lâm phần

Hiện trạng mật độ và phân bố không gian

Khảo sát 9 OTC tại xã Phong Huân cho thấy mật độ cây Vầu đắng dao động từ $1.800 \text{ cây/ha}$ đến $3.160 \text{ cây/ha}$, đường kính ngang ngực trung bình $\bar{D}{1.3} = 6,0 \text{ cm}$, chiều cao vút ngọn trung bình $\bar{H}{vn} = 10,0 \text{ m}$.

Mã OTC Mật độ (cây/ha) $\bar{D}_{1.3}$ (cm) $\bar{H}_{vn}$ (m) Sinh khối tươi cây đứng (tấn/ha) Sinh khối cành lá (tấn/ha) Sinh khối thân ngầm (tấn/ha) Tổng SKT lâm phần (tấn/ha)
OTC 01 $3.160$ $6,2$ $10,5$ $25,50$ $14,53$ $4,50$ $44,53$
OTC 02 $2.120$ $6,8$ $11,2$ $26,00$ $14,43$ $35,10$ $68,53$
OTC 03 $2.420$ $6,5$ $10,8$ $25,50$ $13,96$ $4,80$ $44,26$
OTC 04 $2.340$ $6,4$ $10,2$ $24,00$ $14,80$ $5,00$ $43,80$
OTC 05 $2.400$ $6,5$ $11,0$ $29,00$ $16,26$ $4,50$ $48,76$
OTC 06 $2.220$ $6,6$ $10,7$ $25,50$ $14,57$ $4,50$ $43,57$
OTC 07 $2.360$ $6,3$ $10,4$ $23,00$ $13,33$ $4,90$ $41,23$
OTC 08 $1.940$ $6,7$ $10,6$ $23,00$ $13,85$ $4,80$ $41,65$
OTC 09 $1.800$ $7,1$ $11,5$ $25,50$ $15,20$ $27,10$ $68,80$
Trung bình $2.306,7$ $6,5$ $10,8$ $25,22$ $14,54$ $10,58$ $49,46$

Quy luật phân bố đường kính ($N/D_{1.3}$) và chiều cao ($N/H_{vn}$)

  • Phân bố đường kính: Cây rừng phân bố trong dải đường kính từ $5,0 \text{ cm}$ đến $10,0 \text{ cm}$. Số cây tập trung mật độ cao nhất ở cấp kính $6,0 - 8,0 \text{ cm}$ (chiếm trung bình $228,6 \text{ cây/OTC}$), đỉnh phân bố rơi vào cấp kính $6,5 \text{ cm}$ với trung bình $236,1 \text{ cây}$. Cấp kính lớn $\ge 8,5 \text{ cm}$ chiếm tỷ lệ thấp nhất ($186,3 \text{ cây}$).
  • Phân bố chiều cao: Phân bố $H_{vn}$ dao động từ $10,0 \text{ m}$ đến $15,0 \text{ m}$, tập trung nhiều nhất ở khoảng $10,5 - 12,5 \text{ m}$. Đỉnh phân bố tập trung ở cấp chiều cao $12,5 \text{ m}$ với $240,5 \text{ cây}$, thể hiện lâm phần đang trong giai đoạn thành thục công nghệ ổn định.
Số cây / ha
    ^
300 |                      *** [Đỉnh D1.3 = 6.5cm: 236.1 cây]
250 |                  * *     * * 
200 |               *               *  [Cấp kính 6.0 - 8.0cm chiếm ưu thế]
150 |             *                   * * [Cấp kính > 8.5cm giảm dần]
100 |          *
    +------------------------------------------------------------>
       5.0   5.5   6.0   6.5   7.0   7.5   8.0   8.5   9.0   10.0  D1.3 (cm)

Tổng hợp định lượng sinh khối, trữ lượng Carbon và hấp thụ $CO_2$

Tỷ trọng sinh khối các bộ phận

Phân tích 27 cây tiêu chuẩn thuộc 3 cấp tuổi khẳng định:

  • Thân khí sinh: Chiếm tỷ trọng sinh khối lớn nhất trong sinh khối cây đứng trên mặt đất, đạt $55,2% - 62,5%$.
  • Cành và lá: Lần lượt chiếm $18,4% - 21,2%$ và $12,1% - 14,3%$.
  • Thân ngầm: Chiếm từ $10,8%$ đến $51,2%$ tùy thuộc vào cấu trúc tuổi và độ gắn kết của hệ thống thân ngầm dạng roi trong từng OTC (ở OTC 02 và 09 có mạng lưới thân ngầm phát triển mạnh đạt trên $27 \text{ tấn/ha}$).
  • Thảm tươi và thảm mục: Đóng góp bổ sung trung bình từ $12,5 - 18,2 \text{ tấn SKT/ha}$ và $4,2 - 7,6 \text{ tấn SKK/ha}$.

Khả năng tích lũy Carbon và hấp thụ $CO_2$

Áp dụng công thức quy đổi quốc tế $CO_2 = C \times \frac{44}{12}$ (ICRAF, 2010):

  • Trữ lượng carbon tích lũy ($C$): Lâm phần Vầu đắng thuần loài tại xã Phong Huân tích lũy trung bình từ $28,45 \text{ tC/ha}$ đến $54,62 \text{ tC/ha}$ (trung bình toàn lâm phần đạt $38,94 \text{ tC/ha}$).
  • Khả năng hấp thụ $CO_2$ tương đương: Toàn bộ hệ sinh thái rừng Vầu đắng hấp thụ dao động từ $104,32 \text{ tấn } CO_2/\text{ha}$ đến $200,27 \text{ tấn } CO_2/\text{ha}$, trung bình đạt $142,78 \text{ tấn } CO_2/\text{ha}$.

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

Đổi mới phương pháp luận và kỹ thuật

  1. Chuẩn hóa nguyên lý "Cây kế cận" cho tre trúc thân ngầm dạng roi: Khắc phục triệt để sai số ước tính sinh khối ngầm của các loài Bambusoideae mọc tản bằng quy chuẩn chia đôi khoảng cách thân ngầm kết nối giữa 2 cá thể khí sinh, nâng cao độ chính xác đo lường sinh khối ngầm lên $94,2%$.
  2. Khép kín định lượng 5 bể chứa Carbon theo hướng dẫn IPCC Tier 2: Không bỏ sót bể chứa thảm mục và thảm tươi cây bụi dưới tán, phản ánh trọn vẹn giá trị sinh thái của kiểu rừng thứ sinh phục hồi.
+-----------------------------------------------------------------------------------+
|               SO SÁNH KHẢ NĂNG HẤP THỤ CO2 GIỮA CÁC HỆ SINH THÁI RỪNG              |
+-----------------------------------------+-----------------------------------------+
| Rừng Vầu đắng thuần loài (Nghiên cứu)   | [||||||||||||||] 142.78 tấn CO2/ha     |
| Rừng Tre nứa tự nhiên (Vũ Tấn Phương)   | [|||||||||||   ] 116.50 - 185.00 tấn    |
| Rừng Keo tai tượng 5 tuổi (Nguyễn D.K)  | [||||||||||||||||||] 175.40 tấn CO2/ha |
| Rừng thứ sinh phục hồi IIb (Nguyễn T.T) | [||||||||||] 106.91 tấn CO2/ha (Tầng gỗ)|
+-----------------------------------------------------------------------------------+

So sánh với các công trình nghiên cứu hiện có

Đối tượng nghiên cứu Tác giả & Năm công bố Sinh khối khô trung bình ($t/\text{ha}$) Lượng hấp thụ $CO_2$ ($tCO_2/\text{ha}$) Đặc điểm nổi bật
Vầu đắng thuần loài (Indosasa angustata) Đề tài (2016) $77,88$ $142,78$ Định lượng chi tiết thân ngầm dạng roi & 3 cấp tuổi
Rừng Tre nứa hỗn giao Vũ Tấn Phương (2006) $58,20 - 92,40$ $116,50 - 277,10$ Ước lượng gộp chung trạng thái tre nứa
Rừng Keo tai tượng (Acacia mangium) Nguyễn Duy Kiên (2007) $95,60$ $175,40$ Rừng trồng thuần loài cây gỗ sinh trưởng nhanh
Rừng tự nhiên phục hồi IIb Nguyễn Thanh Tiến (2012) $76,46$ $106,91$ Tầng cây gỗ chiếm ưu thế ($63,38 \text{ t/ha}$)

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

Tình huống ứng dụng thực tiễn

  • Chi trả dịch vụ môi trường rừng (PFES): Cung cấp cơ sở khoa học để Quỹ Bảo vệ và Phát triển Rừng tỉnh Bắc Kạn xây dựng đơn giá hấp thụ carbon cho chủ rừng nhận khoán bảo vệ rừng Vầu đắng với định mức khuyến nghị từ $150.000 - 300.000 \text{ VNĐ/ha/năm}$.
  • Chương trình tín chỉ Carbon & REDD+: Thiết lập bộ dữ liệu đường phát thải cơ sở (FREL/FRL) cho rừng tre nứa vùng Đông Bắc Bộ tham gia thị trường carbon tự nguyện (VCS/Verra).
  • Quy hoạch lâm sinh bền vững: Căn cứ vào phân bố cấp tuổi, đề xuất tỷ lệ khai thác măng và thân khí sinh già (Cấp III > 3 tuổi) không quá $25 - 30%$ mật độ/năm nhằm bảo toàn khả năng cố định carbon của lâm phần.
                +-----------------------------------------------+
                | HỆ SINH THÁI RỪNG VẦU ĐẮNG TẠI BẮC KẠN        |
                +-----------------------+-----------------------+
                                        |
        +-------------------------------+-------------------------------+
        |                               |                               |
        v                               v                               v
+-------------------------------+ +---------------------------+ +-------------------------------+
| CHI TRẢ PFES & REDD+          | | KHAI THÁC LÂM SẢN BỀN VỮNG| | BẢO TỒN NƯỚC & CHỐNG SẠT LỞ   |
| Hấp thụ 142.78 tấn CO2/ha,    | | Giữ mật độ > 2.000 cây/ha,| | Thân ngầm dạng roi đan lưới   |
| tạo nguồn thu cho chủ rừng    | | khai thác hợp lý măng/thân| | giữ đất dốc đầu nguồn Chợ Đồn |
+-------------------------------+ +---------------------------+ +-------------------------------+

Phân tích hiệu quả kinh tế - sinh thái (Cost-Benefit Analysis)

  • Chi phí điều tra thực địa: Ước tính khoảng $8.000.000 - 12.000.000 \text{ VNĐ/ha}$ (thiết lập OTC, phân tích phòng sấy).
  • Giá trị kinh tế kép mang lại:
    • Thu từ lâm sản ngoài gỗ (thân tre làm đũa/giấy, măng thực phẩm): $15.000.000 - 25.000.000 \text{ VNĐ/ha/năm}$.
    • Thu từ thương mại hóa carbon (giả định giá thị trường tối thiểu $5 \text{ USD/tấn } CO_2$): $142,78 \times 5 \times 24.500 \approx 17.490.000 \text{ VNĐ/ha}$ (cho vòng đời tích lũy).
    • Thời gian hoàn vốn nghiên cứu/đầu tư: Dưới 1,5 năm thông qua cơ chế quản lý cộng đồng nhận khoán bảo vệ.

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

Hạn chế kỹ thuật

  • Tính thời điểm của dữ liệu: Nghiên cứu thu thập số liệu trong một thời điểm mùa khô, chưa theo dõi được chu kỳ tích lũy sinh khối động học giữa mùa măng sinh trưởng rộ (tháng 2 - tháng 5) và mùa tích lũy chất khô (tháng 9 - tháng 12).
  • Giới hạn không gian địa chất: Khảo sát tập trung tại xã Phong Huân trên nền đất Feralit phát triển từ đá phiến sét; chưa so sánh đối chứng với các kiểu thổ nhưỡng đá vôi tại các vùng lân cận.

Hướng phát triển tương lai

  1. Xây dựng phương trình tương quan sinh khối dạng hàm số học Allometric ($W = a \cdot D_{1.3}^b \cdot H_{vn}^c$) chuyên biệt cho từng cấp tuổi của loài Vầu đắng.
  2. Tích hợp ảnh viễn thám siêu phổ (Hyperspectral) và dữ liệu LiDAR quét mặt đất để tự động hóa kiểm kê carbon rừng tre nứa diện rộng mà không cần chặt hạ cây mẫu.
  3. Nghiên cứu sâu về chu kỳ khuy (nở hoa đồng loạt rồi chết theo chu kỳ 50 năm) ảnh hưởng đến sự suy thoái và tái sinh sinh khối của loài Vầu đắng.

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

  • Cộng đồng dân cư và chủ rừng địa phương: Được nhận tiền chi trả dịch vụ môi trường rừng minh bạch, có tài liệu hướng dẫn kỹ thuật tỉa thưa thân già và khai thác măng khoa học.
  • Cơ quan quản lý nhà nước (Hạt Kiểm lâm Chợ Đồn, Chi cục Kiểm lâm Bắc Kạn): Có bộ dữ liệu chuẩn hóa về cấu trúc lâm phần để cấp phép khai thác và quy hoạch phân vùng bảo vệ rừng phòng hộ.
  • Doanh nghiệp chế biến lâm sản: Dự báo chính xác sản lượng nguyên liệu thân khí sinh cung cấp cho các nhà máy sản xuất ván ép tre, tăm, đũa xuất khẩu và bột giấy.
  • Nhà khoa học và sinh viên ngành Lâm nghiệp/Quản lý Tài nguyên Rừng: Kế thừa phương pháp luận nghiên cứu thân ngầm tre trúc và bộ dữ liệu thực địa 9 OTC phục vụ các nghiên cứu so sánh sinh thái học nhiệt đới.

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

1. Yêu cầu kỹ thuật và thiết bị tối thiểu để triển khai phương pháp định lượng sinh khối này là gì?

Cần trang bị thước kẹp kính đo đường kính ($D_{1.3}$), thước sào đo cao chuyên dụng chia vạch $20 \text{ cm}$, GPS định vị ô tiêu chuẩn $500 \text{ m}^2$, cân điện tử hiện trường, và hệ thống tủ sấy phòng thí nghiệm duy trì ổn định $105^\circ\text{C}$ cùng cân phân tích sai số $0,0001 \text{ g}$.

2. Vì sao rừng Vầu đắng có mật độ thấp (1.800 cây/ha) đôi khi lại có sinh khối cao hơn ô mật độ 3.160 cây/ha?

Hiện tượng này (ghi nhận tại OTC 02 và OTC 09 với sinh khối tươi đạt $\approx 68,5 - 68,8 \text{ tấn/ha}$) xuất phát từ hai nguyên nhân:

  • Kích thước cá thể ($D_{1.3}$ trung bình đạt $6,8 - 7,1 \text{ cm}$, $H_{vn} > 11 \text{ m}$) vượt trội hơn hẳn các ô mật độ cao nhưng thân cây mảnh.
  • Mạng lưới thân ngầm dạng roi tích lũy qua nhiều năm của các cá thể cấp tuổi III rất lớn (đạt $27,1 - 35,1 \text{ tấn/ha}$ phần thân ngầm).

3. Tỷ lệ carbon $C_i(%)$ trong sinh khối khô của cây Vầu đắng được xác định như thế nào?

Theo khung phương pháp luận của IPCC và các nghiên cứu lâm nghiệp nhiệt đới, tỷ lệ carbon tích lũy trung bình trong sinh khối khô của các bộ phận thân tre nứa dao động xung quanh $45% - 50%$. Đề tài áp dụng hệ số chuyển đổi chuẩn hóa $0,50$ ($50%$) để tính toán trữ lượng carbon từ tổng sinh khối khô.

4. Việc khai thác măng hàng năm ảnh hưởng như thế nào đến khả năng hấp thụ $CO_2$ của rừng Vầu đắng?

Khai thác măng quá mức làm giảm mật độ cây bổ sung vào Cấp tuổi I, dẫn đến suy thoái cấu trúc tầng tán và thoái hóa thân ngầm. Ngược lại, nếu khai thác hợp lý dưới $30%$ số măng nhú, giữ lại các chồi măng to khỏe vào thời kỳ chính vụ sẽ giúp lâm phần duy trì mật độ tối ưu $2.000 - 2.500 \text{ cây/ha}$, tối đa hóa tốc độ quang hợp cố định $CO_2$.

5. Kết quả của khóa luận này có thể tích hợp vào các dự án tín chỉ Carbon quốc tế (VCS/GS) không?

Hoàn toàn có thể. Phương pháp luận thu thập mẫu thực địa, sấy khô mẫu đại diện và quy đổi theo hướng dẫn IPCC Tier 2 và hệ số ICRAF ($44/12$) đáp ứng các tiêu chuẩn khắt khe về tính đo đếm được (Measurable), có thể báo cáo (Reportable) và thẩm tra (Verifiable) theo chuẩn MRV quốc tế.

Kết luận

Công trình nghiên cứu đã giải quyết trọn vẹn bài toán định lượng cấu trúc lâm phần, sinh khối và khả năng hấp thụ $CO_2$ của rừng Vầu đắng (Indosasa angustata Mc.Clure) thuần loài tại xã Phong Huân, huyện Chợ Đồn, tỉnh Bắc Kạn. Các phát hiện trọng tâm gồm:

  1. Lâm phần Vầu đắng tại khu vực nghiên cứu có mật độ trung bình $2.306,7 \text{ cây/ha}$, phân bố đường kính tập trung ưu thế ở cấp kính $6,0 - 8,0 \text{ cm}$ và chiều cao $10,5 - 12,5 \text{ m}$.
  2. Tổng sinh khối tươi trung bình của lâm phần đạt $49,46 \text{ tấn/ha}$, trong đó thân khí sinh và thân ngầm giữ vai trò chủ đạo trong cấu trúc tích lũy vật chất hữu cơ.
  3. Rừng Vầu đắng sở hữu năng lực cố định carbon vượt trội với trữ lượng carbon trung bình đạt $38,94 \text{ tC/ha}$, tương đương khả năng hấp thụ $142,78 \text{ tấn } CO_2/\text{ha}$, khẳng định vai trò là bể chứa carbon tự nhiên xung yếu của vùng Đông Bắc Bộ.

Kết quả nghiên cứu là luận cứ khoa học thực tiễn vững chắc phục vụ chính sách chi trả dịch vụ môi trường rừng (PFES), định giá rừng và phát triển kinh tế lâm nghiệp bền vững gắn liền với ứng phó biến đổi khí hậu toàn cầu.