Giới thiệu dự án

Biến đổi khí hậu toàn cầu đang diễn biến phức tạp với sự gia tăng nồng độ khí nhà kính trong khí quyển, đặc biệt là carbon dioxide ($\text{CO}_2$). Theo các báo cáo khoa học quốc tế, nếu toàn bộ sinh khối rừng nhiệt đới bị suy thoái và đốt cháy trong vòng 50 năm tới, nồng độ $\text{CO}_2$ sẽ tăng gấp đôi, làm nhiệt độ toàn cầu tăng từ $2 - 5^\circ\text{C}$ và mực nước biển dâng từ $1 - 3\text{m}$. Tại Việt Nam, hệ thống hành lang pháp lý về Cơ chế phát triển sạch (CDM), Giảm phát thải từ mất rừng và suy thoái rừng (REDD+) cũng như dịch vụ môi trường rừng đã được hình thành qua Nghị định 48/2007/NĐ-CP, Quyết định 380/QĐ-TTg, Quyết định 158/QĐ-TTg và Nghị định 99/2010/NĐ-CP. Tuy nhiên, việc thực thi chi trả dịch vụ môi trường rừng (PFES) đối với khả năng hấp thụ và lưu giữ carbon vẫn gặp nhiều rào cản do thiếu cơ sở dữ liệu định lượng chính xác cho từng trạng thái rừng cụ thể, đặc biệt là nhóm lâm sản ngoài gỗ.

Vầu đắng (Indosasa angustata Mc. Clure) là loài tre mọc tản bản địa phân bố rộng rãi tại vùng Đông Bắc Bộ, có vai trò to lớn trong phòng hộ đầu nguồn, chống xói mòn và phát triển kinh tế địa phương. Mặc dù giá trị khai thác gỗ tre và măng tre đã được ghi nhận, tiềm năng cố định carbon của loài cây này vẫn chưa được lượng hóa một cách bài bản, dẫn đến việc định giá tài nguyên rừng bị thiếu hụt. Đề tài "Nghiên cứu lượng Carbon tích lũy của rừng Vầu đắng (Indosasa angustata Mc. Clure) tại xã Lam Vỹ, huyện Định Hóa, tỉnh Thái Nguyên" được triển khai nhằm giải quyết trực tiếp khoảng trống khoa học này.

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

  1. Xác định quy luật cấu trúc lâm phần rừng Vầu đắng thuần loài qua phân bố số cây theo đường kính ($N/D$), phân bố số cây theo chiều cao ($N/H$) và tương quan sinh trưởng ($H - D$).
  2. Xác định sinh khối tươi ($\text{SKT}$) và sinh khối khô ($\text{SKK}$) của lâm phần Vầu đắng theo 3 cấp tuổi, tầng cây bụi thảm tươi và tầng vật rơi rụng.
  3. Định lượng hàm lượng và trữ lượng carbon tích lũy trong các bể chứa sinh khối, cung cấp luận cứ khoa học cho chính sách chi trả dịch vụ hấp thụ carbon rừng.

Phạm vi nghiên cứu tập trung vào trạng thái rừng Vầu đắng thuần loài tự nhiên tại xã Lam Vỹ, huyện Định Hóa, tỉnh Thái Nguyên trong khoảng thời gian từ tháng 01/2014 đến tháng 05/2014. Giới hạn đề tài không phân tích bể carbon hữu cơ trong đất mà tập trung toàn diện vào sinh khối trên mặt đất và thân ngầm của tầng cây đứng, cây bụi thảm tươi và thảm mục.

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

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

Nghiên cứu sinh khối và carbon trên thế giới đã phát triển mạnh từ chương trình International Biology Program (1964) và Man and Biosphere (1971). Tại Việt Nam, phần lớn các công trình tập trung vào rừng trồng gỗ lớn như Thông mã vĩ, Keo tai tượng, Keo lai, Bạch đàn hoặc Mỡ. Nhóm tre trúc, điển hình là Vầu đắng, có chu kỳ sinh trưởng ngầm phức tạp, chưa được định lượng đầy đủ.

Phương pháp nghiên cứu Ưu điểm Nhược điểm Khả năng ứng dụng cho Vầu đắng
Ước tính gián tiếp qua thể tích gỗ ($V \times \text{BEF}$) Nhanh, chi phí thấp, không phá hủy lâm phần. Sai số lớn với họ Tre trúc do thân rỗng, không có vòng năm. Không phù hợp
Chặt hạ toàn diện (Destructive Method) Độ chính xác tuyệt đối, cân đo trực tiếp từng bộ phận. Tốn nhân lực, phá hủy một phần lâm phần. Phù hợp với cỡ mẫu chọn lọc
Chặt hạ cây tiêu chuẩn kết hợp giải tích thân ngầm Cân bằng giữa độ chính xác và bảo tồn lâm phần, đo được thân ngầm. Cần kỹ thuật phân cấp tuổi chính xác trên thực địa. Tối ưu nhất (Được lựa chọn) 
                       [Ma trận yêu cầu kỹ thuật MoSCoW]

Thiết kế hệ thống

Khung nghiên cứu định lượng carbon lâm phần được mô hình hóa theo cấu trúc module đa tầng:

Hệ thống công cụ và công nghệ sử dụng trong nghiên cứu:

  • Thước dây sợi thủy tinh bọc nhựa đo chu vi thân chính xác $\pm 0,1\text{cm}$.
  • Thước đo cao chuyên dụng Blumleiss kết hợp sào đo chia vạch $20\text{cm}$, sai số $\pm 10\text{cm}$.
  • Cân hiện trường điện tử định lượng sinh khối tươi tải trọng $50\text{kg}$ (độ nhạy $10\text{g}$) và cân phân tích phòng thí nghiệm Ohaus Pioneer (độ chính xác $0,001\text{g}$).
  • Tủ sấy đối lưu nhiệt Memmert Universal Oven UFE 500 duy trì dải nhiệt độ $70 - 105^\circ\text{C}$.
  • Phần mềm xử lý thống kê toán học: IBM SPSS Statistics 13.0 và Microsoft Excel 2013.

Methodology

Nghiên cứu áp dụng phương pháp tiếp cận điều tra ô tiêu chuẩn điển hình kết hợp phân tích phòng thí nghiệm theo tiêu chuẩn IPCC:

  • Bố trí thí nghiệm: Lập 09 Ô tiêu chuẩn (ÔTC) diện tích $2.500\text{ m}^2$ ($50\text{m} \times 50\text{m}$) trên 3 cấp mật độ (cao: $>5.000\text{ cây/ha}$; trung bình: $3.000 - 5.000\text{ cây/ha}$; thưa: $<3.000\text{ cây/ha}$).
  • Trong mỗi ÔTC lập 05 ô thứ cấp diện tích $25\text{ m}^2$ ($5\text{m} \times 5\text{m}$) tại 4 góc và tâm để thu thập cây bụi thảm tươi (tổng số: 45 ô thứ cấp).
  • Trong mỗi ô thứ cấp lập 01 ô dạng bản $1\text{ m}^2$ ($1\text{m} \times 1\text{m}$) để thu gom vật rơi rụng và thảm mục (tổng số: 45 ô dạng bản).
  • Phân cấp tuổi cá thể Vầu đắng: Do thân ngầm không tạo vòng năm, nghiên cứu kết hợp hình thái vỏ mo, màu sắc lóng thân và tri thức địa phương để phân chia 3 cấp tuổi:
    • Cấp tuổi I: Cây non ($1 - 2$ tuổi).
    • Cấp tuổi II: Cây bánh tẻ ($3 - 4$ tuổi).
    • Cấp tuổi III: Cây già ($> 4$ tuổi).
  • Nguyên lý cây kế cận cho thân ngầm: Đoạn thân ngầm nối giữa hai cây liên tiếp được chia đôi, một nửa được tính là sinh khối thân ngầm của cây tiêu chuẩn chặt hạ.

Implementation và kết quả

Development process

Quá trình thu thập và phân tích dữ liệu thực địa diễn ra qua 4 giai đoạn chặt chẽ:

Tại 09 ÔTC, tiến hành đo đếm toàn bộ đường kính $D_{1.3}$ và chiều cao vút ngọn $H_{vn}$. Lựa chọn 27 cây tiêu chuẩn đại diện (3 cây/ÔTC $\times$ 9 ÔTC) phân bổ đều theo cấp kính và cấp tuổi để chặt hạ. Tách rời 4 bộ phận: Thân, cành, lá và thân ngầm.

Thu thập mẫu phân tích sinh khối khô: Mỗi cây lấy 04 mẫu đại diện gồm thân (gốc $0,0\text{m}$, giữa thân $1/2H$, ngọn $3/4H$), cành, lá và thân ngầm. Mẫu được sấy ở $105^\circ\text{C}$ đến khối lượng không đổi.

Hệ thống công thức toán học tính toán sinh khối và carbon:

# Script mô phỏng thuật toán xử lý sinh khối và trữ lượng Carbon
def calculate_forest_carbon(wt_fresh, sample_dry, sample_fresh, density_n, c_rate=0.5):
    """
    wt_fresh: Sinh khối tươi bộ phận cây mẫu (kg/cây)
    sample_dry: Khối lượng mẫu sau sấy kiệt (g)
    sample_fresh: Khối lượng mẫu tươi ban đầu (g)
    density_n: Mật độ lâm phần (cây/ha)
    c_rate: Tỷ lệ quy đổi Carbon mặc định theo IPCC/thực nghiệm
    """
    # 1. Tính sinh khối khô cây cá thể (Pki)
    p_ki = wt_fresh * (sample_dry / sample_fresh)  # kg/cây
    
    # 2. Quy đổi sinh khối khô trên hecta (Pk_ha)
    pk_ha = (p_ki * density_n) / 1000.0  # tấn/ha
    
    # 3. Tính lượng Carbon tích lũy (Cki)
    c_ki = pk_ha * c_rate  # tấn C/ha
    return p_ki, pk_ha, c_ki

Testing và validation

Số liệu điều tra được kiểm định phân phối chuẩn và phân tích tương quan hồi quy đơn biến $H = f(D)$ trên phần mềm SPSS 13.0.

Ô tiêu chuẩn Phương trình hồi quy sinh trưởng $H - D$ Hệ số tương quan ($R$) Sai số tiêu chuẩn ($S%$) Ý nghĩa thống kê ($F$)
OTC 1 $H_{vn} = 6,9486 + 0,8205 \times D_{1.3}$ $0,87$ $6,42%$ $p < 0,001$
OTC 2 $H_{vn} = 5,1240 + 1,0532 \times D_{1.3}$ $0,94$ $5,11%$ $p < 0,001$
OTC 3 $H_{vn} = 7,3112 + 0,7645 \times D_{1.3}$ $0,82$ $7,85%$ $p < 0,001$
OTC 4 $H_{vn} = 6,1050 + 0,9120 \times D_{1.3}$ $0,89$ $6,18%$ $p < 0,001$
OTC 5 $H_{vn} = 4,8920 + 1,1140 \times D_{1.3}$ $0,97$ $4,20%$ $p < 0,001$
OTC 6 $H_{vn} = 6,4500 + 0,8830 \times D_{1.3}$ $0,88$ $6,72%$ $p < 0,001$
OTC 7 $H_{vn} = 7,0120 + 0,7950 \times D_{1.3}$ $0,79$ $8,90%$ $p < 0,001$
OTC 8 $H_{vn} = 5,8760 + 0,9450 \times D_{1.3}$ $0,91$ $5,64%$ $p < 0,001$
OTC 9 $H_{vn} = 6,2340 + 0,8910 \times D_{1.3}$ $0,85$ $7,15%$ $p < 0,001$

Hệ số tương quan $R$ dao động từ $0,79$ đến $0,97$, chứng minh mối quan hệ giữa chiều cao và đường kính thân cây tuân theo quy luật sinh học chặt chẽ, tạo điều kiện thuận lợi cho việc mô hình hóa sinh khối gián tiếp.

Kết quả đạt được

1. Đặc điểm phân bố lâm phần ($N/D$ và $N/H$)

  • Phân bố $N/D$ có dạng 1 đỉnh lệch trái điển hình của rừng tự nhiên, tập trung cao nhất ở cấp kính $7,0\text{ cm}$ (trung bình $122,4\text{ cây/OTC}$), phân bố từ $4,0\text{ cm}$ ($36,7\text{ cây/OTC}$) đến $10,0\text{ cm}$ ($38,1\text{ cây/OTC}$).
  • Phân bố $N/H$ tập trung dày đặc ở dải chiều cao $12,5 - 14,5\text{ m}$ (trung bình đạt $85,4 - 99,2\text{ cây/OTC}$), thấp nhất ở cấp $10\text{ m}$ ($21,6\text{ cây/OTC}$) và $17\text{ m}$ ($17,4\text{ cây/OTC}$).

2. Tổng hợp sinh khối tươi và sinh khối khô

[Tổng sinh khối khô lâm phần: 62,50 tấn/ha]
  • Sinh khối tươi tầng cây đứng: Cấp tuổi II đạt cao nhất với trung bình $37,5\text{ tấn/ha}$ (thân chiếm $61,2%$, cành $16,0%$, thân ngầm $12,1%$, lá $10,7%$). Cấp tuổi III đạt $34,8\text{ tấn/ha}$ và Cấp tuổi I đạt $31,2\text{ tấn/ha}$.
  • Sinh khối khô tầng cây đứng: Cấp tuổi II đạt $25,0\text{ tấn/ha}$ (thân chiếm $66,1%$, cành $16,9%$, thân ngầm $9,1%$, lá $7,9%$). Cấp III đạt $22,7\text{ tấn/ha}$ và Cấp I đạt $14,8\text{ tấn/ha}$.
  • Cây bụi, thảm tươi: Sinh khối tươi trung bình $3,36\text{ tấn/ha}$ (thân/cành chiếm $72,0%$, lá chiếm $28,0%$); Sinh khối khô trung bình $1,90\text{ tấn/ha}$ (thân/cành $86,0%$, lá $14,0%$).
  • Vật rơi rụng (Thảm mục): Sinh khối tươi trung bình $14,10\text{ tấn/ha}$ (lá/hoa/quả chiếm $65,1%$, cành mục $34,9%$); Sinh khối khô trung bình $8,75\text{ tấn/ha}$ (lá/hoa/quả $70,69%$, cành $29,31%$).

3. Trữ lượng Carbon tích lũy

Cấp tuổi / Bể chứa Lượng Carbon Thân (tấn C/ha) Lượng Carbon Cành (tấn C/ha) Lượng Carbon Lá (tấn C/ha) Lượng Carbon Thân ngầm (tấn C/ha) Tổng trữ lượng Carbon (tấn C/ha)
Cấp tuổi I (1-2 tuổi) $5,94$ ($64,57%$) $2,16$ ($23,74%$) $0,47$ ($5,22%$) $0,59$ ($6,47%$) $9,16$
Cấp tuổi II (3-4 tuổi) $6,55$ ($57,38%$) $3,03$ ($25,68%$) $0,97$ ($8,24%$) $1,03$ ($8,69%$) $11,57$
Cấp tuổi III (>4 tuổi) $6,30$ ($63,11%$) $1,90$ ($18,14%$) $0,85$ ($8,51%$) $1,04$ ($10,24%$) $10,08$
Cây bụi thảm tươi - - - - $0,95$
Thảm mục rơi rụng - - - - $4,38$

Tổng trữ lượng carbon tích lũy của tầng cây Vầu đắng đạt đỉnh ở cấp tuổi II ($11,57\text{ tấn C/ha}$), khẳng định giai đoạn $3 - 4$ tuổi là giai đoạn thành thục công nghệ và hấp thụ carbon mạnh mẽ nhất của lâm phần.

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

  1. Phương pháp luận đo đếm thân ngầm sáng tạo: Ứng dụng thành công nguyên lý "Cây kế cận" (Neighboring Plant Approach) giải quyết triệt để bài toán tách đoạn thân ngầm tre tản vốn là điểm nghẽn kỹ thuật trong điều tra sinh khối lâm nghiệp.
  2. Hệ thống phân cấp tuổi thực nghiệm: Chuẩn hóa phương pháp phân cấp 3 giai đoạn tuổi cho Vầu đắng dựa trên biến đổi hình thái bề mặt lóng thân và độ hóa gỗ, đạt độ tin cậy tương đương với các loài cây có vòng năm.
  3. Bộ dữ liệu gốc về sinh khối tre tản tại Đông Bắc Bộ: Cung cấp định lượng chi tiết cho 4 bộ phận riêng biệt (thân, cành, lá, thân ngầm) thay vì chỉ ước tính phần khí sinh như các nghiên cứu trước đây.
Tiêu chí so sánh Đề tài Vầu đắng (I. angustata) Rừng trồng Keo tai tượng (Nguyễn Duy Kiên) Rừng trồng Bạch đàn Uro (Nguyễn Văn Tấn)
Đối tượng Tre mọc tản tự nhiên Cây gỗ lá rộng mọc nhanh Cây gỗ ngoại lai mọc nhanh
Trữ lượng C cây đứng $9,16 - 11,57\text{ tấn C/ha}$ $55,80 - 74,95\text{ tấn C/ha}$ $25,32 - 37,17\text{ tấn C/ha}$
Chu kỳ tích lũy tối ưu $3 - 4\text{ năm}$ (Rất ngắn) $7 - 10\text{ năm}$ $5 - 6\text{ năm}$
Tỷ lệ sinh khối thân ngầm/rễ $9,1% - 12,5%$ $16,4%$ $16,7% - 25,4%$
Tính bền vững khai thác Khai thác tỉa thưa hàng năm, rừng không bị mất tán Khai thác trắng theo luân kỳ Khai thác trắng theo luân kỳ

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

Kịch bản thương mại hóa tín chỉ Carbon (Carbon Offset)

Với giá thị trường carbon tự nguyện dao động từ $5 - 10\text{ USD/tấn }\text{CO}_2\text{ tương đương}$ ($\text{tCO}_2\text{e}$), khả năng cố định carbon của $1\text{ ha}$ rừng Vầu đắng thuần loài (quy đổi theo tỷ số phân tử $\text{CO}_2/\text{C} = 44/12 \approx 3,67$) đem lại nguồn thu bổ sung đáng kể:

  • Tổng trữ lượng carbon sinh khối trung bình: $\approx 36,14\text{ tấn C/ha}$ (tương đương $\approx 132,63\text{ tCO}_2\text{e/ha}$).
  • Doanh thu dự kiến từ dịch vụ hấp thụ carbon: $660 - 1.320\text{ USD/ha/chu kỳ}$, tạo sinh kế bền vững cho đồng bào dân tộc thiểu số tại huyện Định Hóa.
                  [Lộ trình triển khai cơ chế PFES Carbon]
  Tháng 1-3          Tháng 4-6          Tháng 7-9          Tháng 10-12

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

  • Hạn chế kỹ thuật: Chưa phân tích hàm lượng carbon hữu cơ trong đất tầng sâu ($0 - 30\text{cm}$ và $30 - 60\text{cm}$); số lượng cây chặt hạ thực nghiệm ($27\text{ cây}$) còn chịu giới hạn bởi kinh phí bảo tồn hiện trường.
  • Hướng phát triển:
    1. Xây dựng phương trình tương quan Allometric đa biến biểu diễn sinh khối theo cả $D_{1.3}$ và chiều cao dưới cành $H_{dc}$.
    2. Ứng dụng ảnh viễn thám độ phân giải cao kết hợp cảm biến LiDAR để lập bản đồ sinh khối rừng Vầu tự động trên toàn bộ tỉnh Thái Nguyên.

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

  • Sinh viên & Học viên cao học: Tiếp cận phương pháp luận điều tra lâm học thực địa chuẩn mực, thuật toán phân chia cấp tuổi và quy trình xử lý thống kê mẫu sinh khối rỗng ruột.
  • Cán bộ quản lý lâm nghiệp & Kiểm lâm: Sử dụng bảng tra phân bố $N/D$, $N/H$ và trữ lượng carbon làm định mức phục vụ công tác quy hoạch rừng và giám định chi trả dịch vụ môi trường rừng.
  • Cộng đồng dân cư sống gần rừng: Nâng cao nhận thức bảo vệ tài nguyên rừng Vầu bản địa; tiếp cận nguồn vốn chi trả carbon để tăng thu nhập song song với khai thác măng và thân tre hợp lý.

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

  1. Yêu cầu kỹ thuật cốt lõi để xác định sinh khối rừng Vầu là gì? Cần tuân thủ phân cấp 3 cấp tuổi, lập ô tiêu chuẩn đại diện tối thiểu $2.500\text{ m}^2$, sấy mẫu các bộ phận ở nhiệt độ ổn định $105^\circ\text{C}$ cho đến khi cân 2 lần liên tiếp không đổi khối lượng.
  2. Rừng Vầu đắng có khả năng tích lũy carbon cao hơn rừng gỗ lớn không? Trữ lượng tức thời trên một hecta thấp hơn cây gỗ lớn lâu năm, nhưng tốc độ luân chuyển sinh khối và khả năng hấp thụ hàng năm của rừng Vầu cao hơn nhờ khả năng tái sinh măng chồi liên tục mà không cần trồng lại sau khai thác.
  3. Tại sao cấp tuổi II (3-4 tuổi) lại có trữ lượng carbon cao nhất? Ở giai đoạn $3 - 4$ tuổi, tế bào vách sợi tre đã hóa gỗ hoàn toàn, mật độ thành phần cơ lý đạt cực đại trong khi mật độ cây cá thể của cấp tuổi này trong lâm phần tự nhiên chiếm tỷ trọng lớn nhất ($37,5\text{ tấn tươi/ha}$).
  4. Việc khai thác măng có làm suy giảm trữ lượng carbon của rừng không? Khai thác tối đa $50%$ lượng măng nhú đầu mùa không làm suy giảm sinh khối rừng, vì trong tự nhiên khoảng $50%$ số măng non sẽ tự đào thải do cạnh tranh dinh dưỡng dưới tán.
  5. Dữ liệu này có thể tích hợp vào hệ thống MRV quốc gia không? Hoàn toàn khả thi. Bộ dữ liệu cung cấp hệ số phát thải và bể hấp thụ cấp địa phương (Tier 2/Tier 3) phục vụ kiểm kê khí nhà kính ngành LULUCF theo hướng dẫn của IPCC.

Kết luận

Công trình nghiên cứu đã lượng hóa thành công các quy luật kết cấu lâm học, đặc điểm tích lũy sinh khối và trữ lượng carbon của rừng Vầu đắng (Indosasa angustata Mc. Clure) tại huyện Định Hóa, tỉnh Thái Nguyên. Kết quả khẳng định rừng Vầu đắng đạt đỉnh tích lũy carbon ở cấp tuổi II ($11,57\text{ tấn C/ha}$), tổng sinh khối khô lâm phần đạt $62,50\text{ tấn/ha}$. Đây là cơ sở khoa học trực tiếp phục vụ công tác định giá rừng, thực thi Nghị định 99/2010/NĐ-CP và mở ra tiềm năng thương mại hóa tín chỉ carbon rừng tre nứa tại Việt Nam.