Tổng quan nghiên cứu

Biến đổi khí hậu toàn cầu đang diễn ra gay gắt với nồng độ khí nhà kính gia tăng, đòi hỏi các giải pháp khoa học nhằm nâng cao năng lực hấp thụ carbon từ hệ sinh thái rừng. Trong cơ chế Giảm phát thải từ mất rừng và suy thoái rừng (REDD+) cũng như dịch vụ chi trả môi trường rừng (PES), việc lượng hóa chính xác trữ lượng carbon là điều kiện tiên quyết để thương mại hóa tín chỉ carbon. Tại Việt Nam, cây Luồng (Dendrocalamus barbatus) là loài lâm sản ngoài gỗ quan trọng với diện tích hơn 10.000 ha tại riêng huyện Lang Chánh (Thanh Hóa) và hàng nghìn hecta tại tỉnh Hòa Bình. Tuy nhiên, các tính toán bể chứa carbon trước đây chủ yếu áp dụng hệ số quy đổi mặc định 0,50 (50% sinh khối khô) của Ủy ban Liên chính phủ về Biến đổi Khí hậu (IPCC) vốn dành cho cây thân gỗ, dẫn đến độ lệch số liệu đáng kể khi áp dụng cho loài tre trúc sinh trưởng nhanh.

Nghiên cứu được thực hiện tại Lâm trường Lương Sơn (Hòa Bình, diện tích quản lý 2.467,6 ha) và Công ty Lâm nghiệp Lang Chánh (Thanh Hóa) trong năm 2015 nhằm giải quyết nút thắt khoa học này. Mục tiêu cụ thể là xác định đặc điểm sinh trưởng, cấu trúc sinh khối tươi, sinh khối khô và phân tích chính xác hàm lượng carbon thực tế tích lũy trong từng bộ phận (thân khí sinh, cành, lá, thân ngầm, rễ) của cây Luồng qua 4 cấp tuổi (tuổi 1, 2, 3 và từ 4 trở lên). Kết quả cung cấp cơ sở dữ liệu thực nghiệm chuẩn xác, hỗ trợ nâng cao độ tin cậy của các báo cáo kiểm kê phát thải quốc gia và tối ưu hóa giá trị kinh tế của rừng Luồng thông qua thị trường tín chỉ carbon tự nguyện.

Cơ sở lý thuyết và phương pháp nghiên cứu

Khung lý thuyết áp dụng

Nghiên cứu vận dụng lý thuyết tích lũy sinh khối và chu chuyển carbon trong hệ sinh thái rừng nhiệt đới của Mc Kenzie (2001) kết hợp hướng dẫn kiểm kê khí nhà kính của IPCC (2003, 2006). Hệ thống lý thuyết xác định 5 bể chứa carbon cốt lõi, trong đó sinh khối sống trên mặt đất và dưới mặt đất của thực vật giữ vai trò hấp thụ $CO_2$ chủ động. Bên cạnh đó, đề tài kế thừa quy luật phân bố sinh khối thực vật rừng theo các mô hình phân tích sinh thái học của Cannell (1982) và Newbould (1967).

Ba khái niệm học thuật trung tâm được chuẩn hóa xuyên suốt gồm:

  1. Sinh khối tươi (SKT): Tổng khối lượng vật chất hữu cơ tươi của từng bộ phận cây tại thời điểm chặt hạ ngoài hiện trường.
  2. Sinh khối khô (SKK): Khối lượng chất khô tuyệt đối sau khi sấy mẫu ở nhiệt độ $105^\circ\text{C}$ đến khối lượng không đổi, biểu thị lượng chất khô hữu cơ tích lũy thực tế.
  3. Hệ số tích lũy carbon ($K%$): Tỷ lệ phần trăm carbon nguyên chất có trong một đơn vị khối lượng sinh khối khô, được xác lập qua phân tích đốt cháy buồng kín thay cho việc ước tính cố định 0,50.

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

Nguồn dữ liệu ngoại nghiệp được thu thập thông qua hệ thống 6 ô tiêu chuẩn (OTC) tạm thời có kích thước $25\text{m} \times 40\text{m}$ (diện tích $1.000,\text{m}^2/\text{OTC}$), bố trí đại diện cho 3 cấp lập địa tốt, trung bình và xấu tại Lương Sơn và Lang Chánh. Phương pháp chọn mẫu điển hình được áp dụng để hạ 24 cây mẫu tiêu chuẩn (12 cây tại mỗi vùng nghiên cứu), phân tầng đều đặn cho 4 cấp tuổi sinh trưởng từ tuổi 1 đến tuổi 4 trở lên nhằm bao quát quy luật tích lũy vật chất hữu cơ theo tuổi thọ lâm phần.

Để khảo sát sinh khối rễ dưới mặt đất, đề tài thiết lập 4 ô dạng bản (ODB) kích thước $0,5\text{m} \times 0,5\text{m}$ xuyên tâm bụi Luồng, tiến hành đào sâu $30,\text{cm}$ theo hướng đồng mức và vuông góc đường đồng mức. Phân tích nội nghiệp được thực hiện tại Phòng thí nghiệm Đất và Môi trường thuộc Viện Nghiên cứu Sinh thái và Môi trường rừng. Các mẫu vật sau khi sấy khô kiệt được nghiền mịn thành bột, cân chính xác $50,\text{mg}$ mẫu và phân tích trên hệ thống máy buồng đốt kín Multi N/C 3000 ở nhiệt độ cao có chất xúc tác $CeO_2$. Khí $CO_2$ giải phóng được định lượng bằng đầu đo quang phổ hồng ngoại không phân tán (NDIR) để xác định hàm lượng carbon tuyệt đối. Dữ liệu được xử lý bằng kiểm định Chi-bình phương ($\chi^2$) để đánh giá phẩm chất lâm phần và phân tích phương sai một nhân tố (One-way ANOVA) nhằm so sánh độ sai lệch giữa hệ số carbon thực nghiệm và hệ số lý thuyết IPCC. Toàn bộ chu trình khảo sát và phân tích phòng thí nghiệm được hoàn thành đồng bộ trong năm 2015.

Kết quả nghiên cứu và thảo luận

Những phát hiện chính

Kết quả đo đếm lâm học ghi nhận sự khác biệt rõ nét về sinh trưởng giữa hai vùng sinh thái. Tại Lang Chánh (Thanh Hóa), đường kính ngang ngực ($D_{1.3}$) trung bình đạt $7,66,\text{cm}$ và chiều cao vút ngọn ($H_{vn}$) đạt $10,17,\text{m}$, vượt trội hơn so với Lương Sơn (Hòa Bình) có $D_{1.3}$ trung bình $5,80,\text{cm}$ và $H_{vn}$ trung bình $9,08,\text{m}$ (đường kính cao hơn $32,07%$ và chiều cao lớn hơn $12,00%$). Cây đạt kích thước cực đại tại ô tiêu chuẩn 1 ở Lang Chánh tại tuổi 4 với $D_{1.3} = 8,57,\text{cm}$ và $H_{vn} = 11,93,\text{m}$.

Về phẩm chất lâm phần, cấu trúc chất lượng cây tại Lương Sơn có tính đồng nhất cao giữa các ô tiêu chuẩn với tỷ lệ cây sinh trưởng tốt chiếm từ $42,74%$ đến $46,62%$ (trung bình $44,37%$), kiểm định $\chi^2_n = 1,44 < \chi^2_{0.05} = 9,49$. Ngược lại, tại Lang Chánh xuất hiện sự phân hóa mạnh mẽ theo vi địa hình ($\chi^2_n = 17,21 > 9,49$), trong đó tỷ lệ cây phẩm chất xấu ở ô tiêu chuẩn 2 lên tới $37,44%$ do hiện tượng cụt ngọn và sâu bệnh.

Về cấu trúc sinh khối tươi (trừ rễ), tổng khối lượng trung bình mỗi cây tại Lang Chánh đạt $22,68,\text{kg}$, cao hơn mức $19,95,\text{kg}$ ở Lương Sơn. Tỷ trọng phân bố sinh khối tươi tập trung cao nhất ở thân khí sinh (Lương Sơn: $11,28,\text{kg/cây}$ chiếm $56,51%$; Lang Chánh: $13,66,\text{kg/cây}$ chiếm $60,24%$), tiếp theo là cành ($17,68% - 20,68%$), thân ngầm ($12,91% - 15,03%$) và thấp nhất ở lá ($7,06% - 9,90%$). Khi chuyển sang sinh khối khô, tỷ lệ chất khô trong thân khí sinh và thân ngầm tăng lên đáng kể do hàm lượng xơ sợi lignin cao, trong khi hàm lượng carbon thực nghiệm trong các bộ phận dao động từ $45,20%$ đến $49,15%$, không đạt mức cố định $50,00%$ như quy ước truyền thống.

Thảo luận kết quả

Sự chênh lệch về sinh trưởng và khả năng tích lũy sinh khối giữa hai địa bàn xuất phát từ điều kiện lập địa và cường độ xói mòn đất. Huyện Lang Chánh có tầng đất phong hóa feralit dày, lượng mưa dồi dào khoảng $2.500,\text{mm/năm}$, tạo điều kiện tối ưu cho hệ thân ngầm phát triển chồi măng. Tuy nhiên, trên những sườn đồi có độ dốc từ $20^\circ$ đến $50^\circ$, hiện tượng rửa trôi bề mặt diễn ra mạnh mẽ làm suy giảm độ phì, dẫn đến việc gia tăng tỷ lệ cây xấu nếu không có biện pháp bảo vệ đất. Tại cả hai khu vực, mức sinh khối cá thể nhìn chung vẫn thấp hơn tiềm năng di truyền của loài do tập quán canh tác chỉ chú trọng thu hoạch măng vụ xuân mà bỏ qua các biện pháp lâm sinh thâm canh như bón phân hay xới xáo cải tạo đất.

Khi so sánh với các công trình nghiên cứu trên rừng tự nhiên của GS. Võ Đại Hải (2012) và rừng lá rộng của Bảo Huy (2009), hàm lượng carbon trong cây Luồng phân bố không đều giữa các cơ quan. Việc áp dụng hệ số mặc định 0,50 của IPCC làm sai lệch tổng trữ lượng carbon của lâm phần rừng Luồng khoảng từ $3,5%$ đến $7,8%$ tùy thuộc vào cơ cấu tuổi và tỷ trọng sinh khối thân - ngầm. Trong báo cáo khoa học, dữ liệu này có thể được trực quan hóa hiệu quả thông qua biểu đồ cột đa trục mô tả mối quan hệ giữa sinh khối khô từng bộ phận theo 4 cấp tuổi và đồ thị hồi quy tương quan phi tuyến giữa đường kính bụi, chiều cao với lượng hấp thụ carbon quy đổi.

Đề xuất và khuyến nghị

  1. Chuẩn hóa hệ số phát thải và chuyển đổi carbon tre nứa: Đề nghị Cục Lâm nghiệp và các cơ quan thẩm định tín chỉ REDD+ cập nhật dải hệ số tích lũy carbon thực nghiệm theo từng bộ phận ($45,20% - 49,15%$) thay vì áp dụng hệ số cứng 0,50. Giải pháp này cần hoàn thành trong vòng 12 đến 24 tháng nhằm giảm thiểu sai số kiểm kê khí nhà kính quốc gia xuống dưới $2,0%$.
  2. Áp dụng gói kỹ thuật lâm sinh thâm canh phục tráng rừng Luồng: Các Công ty Lâm nghiệp tại Hòa Bình và Thanh Hóa cần tổ chức bón phân NPK cân đối với liều lượng $0,5 - 1,0,\text{kg/bụi/năm}$, kết hợp xới đất sâu $15 - 20,\text{cm}$ và vun gốc định kỳ vào đầu mùa mưa (tháng 3 đến tháng 4). Mục tiêu kỹ thuật là nâng tốc độ tăng trưởng đường kính lên trên $8,0,\text{cm}$ và tăng tổng sinh khối lâm phần thêm $20% - 25%$ trong giai đoạn 2026 - 2028.
  3. Điều chỉnh quy chế khai thác và tối ưu hóa cơ cấu cấp tuổi: Chủ rừng cần duy trì tỷ lệ cơ cấu cây trong bụi theo tỷ lệ tuổi 1, 2, 3, 4 ở mức cân bằng 1:1:1:1, khống chế cường độ khai thác hàng năm không quá $30%$ tổng số thân cây trong bụi. Biện pháp này giúp nâng tỷ lệ cây chất lượng tốt (loại A) đạt trên $60%$ trên toàn lâm phần trước năm 2029.
  4. Xây dựng dự án thí điểm chi trả dịch vụ môi trường rừng gắn với tín chỉ carbon: Ủy ban nhân dân hai tỉnh Thanh Hóa và Hòa Bình cần chủ trì liên kết giữa các hợp tác xã lâm nghiệp với các tổ chức quốc tế để xây dựng hồ sơ thương mại hóa tín chỉ carbon từ 10.000 ha rừng Luồng, phấn đấu hiện thực hóa nguồn thu bổ sung đạt khoảng 5 đến 10 USD/tấn $CO_2$ tương đương trong giai đoạn 2026 - 2030.

Đối tượng nên tham khảo luận văn

  • Cơ quan quản lý nhà nước và nhà hoạch định chính sách lâm nghiệp: Cung cấp bộ thông số đầu vào chính xác để xây dựng khung giá chi trả dịch vụ môi trường rừng (PES) và định mức kỹ thuật quốc gia về đo đếm carbon rừng tre nứa trong khuôn khổ Công ước khung Liên Hợp Quốc về Biến đổi Khí hậu.
  • Ban quản lý rừng phòng hộ, công ty lâm nghiệp và hợp tác xã trồng Luồng: Giúp cán bộ kỹ thuật nắm bắt quy luật sinh khối theo cấp tuổi để thiết lập quy trình tỉa thưa, điều tiết mật độ cây đứng và kéo dài chu kỳ tích lũy vật chất hữu cơ của lâm phần.
  • Giảng viên, nghiên cứu sinh và học viên cao học ngành Lâm học, Quản lý tài nguyên rừng: Cung cấp một tài liệu tham khảo thực nghiệm chuẩn mực về quy trình giải phẫu mẫu sinh khối tre nứa, kỹ thuật vận hành máy quang phổ buồng đốt kín Multi N/C 3000 và phương pháp xử lý thống kê sinh học.
  • Doanh nghiệp tư vấn và tổ chức phát triển dự án carbon tự nguyện: Hỗ trợ các chuyên gia tính toán đường phát thải cơ sở (carbon baseline) khoa học, loại bỏ rủi ro bị bác bỏ tín chỉ do sử dụng sai hệ số quy đổi trong các tiêu chuẩn quốc tế như VCS hay Gold Standard.

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

Tại sao không nên dùng hệ số quy đổi 0,50 của IPCC cho rừng Luồng?

Hệ số 0,50 được IPCC tính toán trung bình cho các loài cây thân gỗ có tỷ lệ hóa gỗ cao. Kết quả phân tích nhiệt phân thực nghiệm cho thấy hàm lượng carbon trong các bộ phận cây Luồng chỉ dao động trong khoảng $45,20% - 49,15%$. Áp dụng hệ số 0,50 sẽ gây hiện tượng phóng đại trữ lượng carbon thực tế từ $3,5%$ đến $7,8%$, làm suy giảm tính chính xác của hồ sơ tín chỉ carbon.

Bộ phận nào của cây Luồng tích lũy khối lượng carbon cao nhất?

Thân khí sinh chiếm tỷ trọng carbon vượt trội nhất với khối lượng tươi trung bình từ $11,28,\text{kg}$ đến $13,66,\text{kg/cây}$ (chiếm $56,51% - 60,24%$ tổng sinh khối tươi của cây). Sau đó đến cành chiếm khoảng $17,68% - 20,68%$ và hệ thân ngầm chiếm $12,91% - 15,03%$.

Vì sao sinh trưởng của Luồng tại Thanh Hóa lại vượt trội hơn ở Hòa Bình?

Khu vực Lang Chánh (Thanh Hóa) có điều kiện khí hậu nhiệt đới nóng ẩm tối ưu với tổng tích ôn năm khoảng $7.000^\circ\text{C}$ và lượng mưa đạt $2.500,\text{mm/năm}$, kết hợp nền đất phong hóa sâu trên đá trầm tích, giúp cây Luồng đạt đường kính trung bình $7,66,\text{cm}$, lớn hơn đáng kể so với mức $5,80,\text{cm}$ tại Lương Sơn (Hòa Bình).

Độ tuổi nào của cây Luồng có khả năng tích lũy sinh khối tốt nhất?

Sinh khối của cây Luồng tăng dần theo cấp tuổi và đạt giá trị tối đa ở tuổi 4 trở lên với sinh khối tươi cá thể đạt từ $22,3,\text{kg}$ đến $29,5,\text{kg/cây}$. Do đó, việc duy trì tỷ lệ cây tuổi 3 và tuổi 4 hợp lý trong bụi là yếu tố then chốt để tối đa hóa trữ lượng carbon rừng.

Phương pháp đốt mẫu kín bằng máy Multi N/C 3000 có ưu thế gì so với phương pháp cổ điển?

Phương pháp đốt mẫu buồng kín với xúc tác $CeO_2$ ở nhiệt độ cao chuyển hóa triệt để toàn bộ carbon hữu cơ thành khí $CO_2$ và đo đạc trực tiếp qua cảm biến hồng ngoại NDIR. Phương pháp này loại bỏ hoàn toàn sai số do oxy hóa không triệt để của phương pháp Walkley-Black ướt truyền thống, mang lại độ chính xác đạt ngưỡng 0,1%.

Kết luận

  • Đề tài đã hoàn thành đo đếm thực nghiệm và xác định rõ kích thước sinh trưởng trung bình của Luồng tại Lang Chánh ($D_{1.3} = 7,66,\text{cm}$; $H_{vn} = 10,17,\text{m}$) vượt trội hơn khu vực Lương Sơn ($D_{1.3} = 5,80,\text{cm}$; $H_{vn} = 9,08,\text{m}$).
  • Xác lập được cấu trúc sinh khối tươi cá thể (trừ rễ) với thân khí sinh chiếm tỷ trọng lớn nhất ($56,51% - 60,24%$), cành ($17,68% - 20,68%$), thân ngầm ($12,91% - 15,03%$) và lá ($7,06% - 9,90%$).
  • Chứng minh hàm lượng carbon thực nghiệm trong các bộ phận cây Luồng dao động trong dải $45,20% - 49,15%$, bác bỏ tính chính xác tuyệt đối của việc áp đặt hệ số mặc định 0,50 cho loài tre trúc.
  • Đóng góp luận cứ khoa học thực tiễn vững chắc cho việc hiệu chỉnh công thức tính toán trữ lượng carbon rừng tre nứa trong khuôn khổ các chương trình REDD+ và dịch vụ môi trường rừng tại Việt Nam.
  • Đề xuất lộ trình 2026 - 2030 tập trung chuẩn hóa hướng dẫn kỹ thuật cấp bộ và nhân rộng mô hình thâm canh thương mại hóa tín chỉ carbon rừng Luồng tại các tỉnh miền núi phía Bắc.

Hãy chủ động áp dụng hệ số chuyển đổi carbon chuyên biệt từ công trình nghiên cứu này để xây dựng các đề án phát triển rừng bền vững và gia tăng giá trị kinh tế xanh từ tài nguyên rừng Luồng Việt Nam ngay hôm nay!