Tổng quan nghiên cứu

Trong giai đoạn 1943 - 1990, diện tích rừng tự nhiên tại Việt Nam bị suy giảm nghiêm trọng khoảng 5 triệu ha, gây ảnh hưởng tiêu cực đến chu trình sinh thái và làm gia tăng phát thải khí nhà kính toàn cầu. Bước vào giai đoạn phục hồi 1990 - 2005, tổng diện tích rừng cả nước đã tăng lên xấp xỉ 12,6 triệu ha với độ che phủ đạt 38%. Cùng với xu thế toàn cầu sau sự ra đời của Nghị định thư Kyoto (1997), ngành lâm nghiệp chuyển đổi mạnh mẽ từ tư duy khai thác gỗ truyền thống sang công nhận và thương mại hóa các giá trị dịch vụ môi trường rừng, tiêu biểu là khả năng hấp thụ và lưu giữ cacbon (quy đổi thành tín chỉ CO2).

Tại huyện Hoàng Su Phì, tỉnh Hà Giang, loài Thông ba lá (Pinus kesiya) được đưa vào trồng tập trung từ những năm 1990 với quy mô đạt khoảng 2.400 ha, chiếm 8,6% tổng diện tích rừng toàn huyện. Mặc dù đóng vai trò quan trọng trong việc phòng hộ đầu nguồn sông Chảy và phát triển kinh tế địa phương, các dữ liệu khoa học định lượng về sinh khối và bể chứa cacbon của loài cây này tại vùng Đông Bắc vẫn còn là một khoảng trống lớn.

Luận văn thạc sĩ khoa học lâm nghiệp của tác giả Nguyễn Thanh Hải do Giáo sư, Tiến sĩ Khoa học Đỗ Đình Sâm hướng dẫn (thực hiện năm 2011 tại Trường Đại học Lâm nghiệp) đã tập trung giải quyết vấn đề cốt lõi này. Mục tiêu cụ thể của nghiên cứu là xác định sinh khối, trữ lượng cacbon ở các cấp tuổi từ 5 đến 26, xây dựng hệ thống mô hình toán học dự đoán sinh khối qua đường kính ngang ngực, đồng thời lượng giá kinh tế giá trị hấp thụ cacbon của rừng Thông ba lá. Với tổng trữ lượng cacbon hệ sinh thái đạt từ 66,5 đến 144,3 tấn C/ha, công trình cung cấp cơ sở dữ liệu quan trọng cho các chương trình chi trả dịch vụ môi trường rừng (PES), cơ chế REDD+ và công tác kiểm kê khí nhà kính quốc gia.

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

Khung lý thuyết áp dụng

Khung lý thuyết của nghiên cứu được phát triển dựa trên nguyên lý chu trình sinh địa hóa cacbon, lý thuyết năng suất sinh học sơ cấp và hướng dẫn kiểm kê khí nhà kính của Ban Liên chính phủ về Biến đổi Khí hậu (IPCC 2003, 2006). Theo đó, tổng trữ lượng cacbon của hệ sinh thái rừng được cấu thành từ 4 bể chứa chính: sinh khối cây sống trên mặt đất và dưới mặt đất (AGB và BGB), thảm mục/vật rơi rụng (DWB), thảm tươi cây bụi, và lượng cacbon hữu cơ trong đất (SOC).

Các khái niệm then chốt được vận dụng gồm: Đường kính ngang ngực (DBH đo tại vị trí 1,3 m), Hệ số mở rộng sinh khối (BEF), Tăng trưởng trung bình hàng năm (MAI), và Hệ số quy đổi phân tử từ cacbon nguyên tố sang CO2 tương đương (hệ số 3,67). Mô hình nghiên cứu kế thừa các phương trình alometric dạng hàm lũy thừa phi tuyến $Y = a \cdot X^b$. Điểm đặc thù của luận văn là xác định hàm lượng cacbon thực tế qua phân tích hóa nghiệm cho từng bộ phận cây, thay vì áp dụng cứng nhắc hệ số mặc định 0,45 hoặc 0,5 của IPCC, giúp phản ánh chân thực đặc tính sinh thái của Thông ba lá tại vùng núi cao phía Bắc.

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

Nghiên cứu áp dụng quy trình điều tra lâm học kết hợp thực nghiệm phân tích sinh khối theo tiêu chuẩn phối hợp của Trung tâm Hợp tác Quốc tế Nhật Bản và Viện Khoa học Lâm nghiệp (FORDA & JICA, 2005). Cỡ mẫu được thiết lập gồm 27 ô tiêu chuẩn lâm phần (kích thước 20m x 20m, diện tích 400 m2/ô) được bố trí ngẫu nhiên đại diện cho 4 cấp tuổi: 5 tuổi (mật độ 1.600 cây/ha), 12 tuổi (1.100 cây/ha), 17 tuổi (700 cây/ha) và 26 tuổi (500 cây/ha).

Phương pháp chọn mẫu cây tiêu chuẩn dựa trên chỉ số đường kính bình quân theo tiết diện ngang (Dg) của từng ô đo đếm nhằm đảm bảo tính đại diện thống kê cao nhất. Cây mẫu được chặt giải tích, cân trọng lượng tươi ngoài hiện trường và tách riêng 5 bộ phận: thân, vỏ, cành, lá, rễ. Mẫu sinh khối được sấy khô kiệt ở nhiệt độ 75 - 85 độ C trong 6 - 8 giờ tại phòng thí nghiệm để xác định tỷ lệ chất khô.

Sinh khối thảm mục và tầng cây bụi thảm tươi được thu thập trong các ô dạng bản 1 m2. Mẫu đất được lấy ở hai tầng sâu 0 - 10 cm và 10 - 30 cm bằng ống dung trọng 50 cm3. Hàm lượng cacbon trong mẫu đất và thực vật được xác định bằng phương pháp oxy hóa ướt Walkley - Black. Phương pháp này được lựa chọn vì tính chính xác cao, loại bỏ sai số do tạp chất vô cơ và đạt chuẩn hóa quốc tế. Toàn bộ dữ liệu được phân tích hồi quy tương quan trên phần mềm SPSS và Microsoft Excel, kiểm tra độ tin cậy tham số bằng kiểm định Student (t) và kiểm định Fisher (F) ở mức ý nghĩa $P < 0,05$. Quá trình thu thập và xử lý số liệu hoàn thành trọn vẹn trong năm 2011.

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

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

Nghiên cứu mang lại 4 phát hiện thực nghiệm mang tính định lượng cao về cấu trúc sinh khối và khả năng tích lũy cacbon của rừng trồng Thông ba lá:

Thứ nhất, tổng sinh khối khô cá thể tăng trưởng vượt bậc theo cấp tuổi, đạt 31,01 kg/cây ở tuổi 5; 90,79 kg/cây ở tuổi 12; 227,94 kg/cây ở tuổi 17 và đạt 645,38 kg/cây ở tuổi 26. Tốc độ tăng trưởng sinh khối bình quân tăng từ 6,20 kg/cây/năm lên 24,82 kg/cây/năm. Tỷ lệ sinh khối thân chiếm ưu thế lớn nhất và tăng dần từ 39% ở tuổi 5 lên 68% tổng sinh khối ở tuổi 26. Tỷ lệ cành giảm dần từ 28,4% (tuổi 5) xuống 9,5% (tuổi 26). Tỷ lệ sinh khối rễ duy trì ổn định trong khoảng 14,5% - 19,5%, tiệm cận tỷ lệ chuẩn 4:1 giữa phần trên mặt đất và dưới mặt đất theo khuyến cáo của IPCC.

Thứ hai, hàm lượng cacbon trong các bộ phận thực vật dao động từ 36,90% đến 43,36%, cao nhất ở thân cây (43,36%), tiếp theo là rễ (41,00%), lá (40,61%), cành (40,21%), thảm mục (38,42%) và thấp nhất ở vỏ cây (36,90%). Trữ lượng cacbon cá thể tích lũy tăng từ 12,64 kg C/cây ở tuổi 5 lên 259,17 kg C/cây ở tuổi 26, trong đó cacbon tích lũy trong thân chiếm 41% đến 69%.

Thứ ba, tổng trữ lượng cacbon của lâm phần dao động từ 66,5 tấn C/ha đến 144,3 tấn C/ha. Tỷ trọng cacbon trên mặt đất mở rộng theo tuổi từ 38,4% lên 80,7% (tương đương 25,5 đến 116,4 tấn C/ha). Ngược lại, tỷ lệ cacbon dưới mặt đất giảm mạnh từ 160,66% xuống còn 23,92% so với phần trên mặt đất. Trữ lượng thảm mục dao động 1,69 - 2,98 tấn/ha và cây bụi thảm tươi đạt 4,93 - 14,28 tấn/ha.

Thứ tư, phương trình tương quan lũy thừa $Y = a \cdot X^b$ giữa sinh khối khô, trữ lượng cacbon với đường kính ngang ngực (DBH) thể hiện độ tương thích rất cao với hệ số tương quan $r \ge 0,99$ và mức ý nghĩa thống kê $Sig F < 0,001$.

Thảo luận kết quả

Khi so sánh với các nghiên cứu lâm học tại Tây Nguyên, sinh khối cá thể Thông ba lá tại Hoàng Su Phì đạt mức tương đương với rừng trồng ở cấp đất I tại Lâm Đồng, Gia Lai và Kon Tum (ở tuổi 5 - 25 dao động 30,4 - 500,1 kg/cây). Tuy nhiên, tổng trữ lượng cacbon trên mỗi hecta của lâm phần lại chỉ tương đương cấp đất IV và V (43,89 - 144,31 tấn C/ha). Nguyên nhân của sự chênh lệch này là do mật độ lâm phần tại Hoàng Su Phì giảm nhanh theo thời gian (từ 1.600 cây/ha xuống 500 cây/ha ở tuổi 26) cùng với địa hình đồi núi dốc trên 25 độ khiến đất bị rửa trôi, làm giảm trữ lượng cacbon tích lũy trong tầng đất 0 - 30 cm (dao động 27,8 - 41,0 tấn C/ha).

Để minh họa trực quan, toàn bộ cơ cấu sinh khối và tỷ lệ tích lũy cacbon của các bộ phận có thể được biểu diễn trực quan qua biểu đồ tròn phân đoạn tỷ lệ phần trăm và bảng so sánh đa biến theo cấp tuổi. Mối quan hệ giữa đường kính ngang ngực và sinh khối được mô hình hóa tối ưu qua đồ thị phân tán XY dạng đường cong hồi quy lũy thừa, cho phép các cơ quan kiểm lâm và chủ rừng tra cứu nhanh trữ lượng mà không phải giải tích phá hủy cây rừng.

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

Nhằm nâng cao hiệu quả hấp thụ cacbon và tối đa hóa giá trị kinh tế của rừng trồng Thông ba lá, luận văn đề xuất 4 nhóm giải pháp cụ thể:

Thứ nhất, tối ưu hóa mật độ lâm phần và kỹ thuật tỉa thưa. Ủy ban nhân dân huyện Hoàng Su Phì và Ban Quản lý rừng phòng hộ cần chỉ đạo duy trì mật độ ban đầu từ 1.600 đến 2.000 cây/ha ở giai đoạn 1 - 5 năm đầu, thực hiện tỉa thưa nuôi dưỡng ở giai đoạn 10 - 15 tuổi để giữ mật độ 700 - 800 cây/ha, đảm bảo mục tiêu nâng tổng trữ lượng cacbon lâm phần đạt trên 160 tấn C/ha ở tuổi 20. Quy trình tỉa thưa cần thực hiện định kỳ 5 năm một lần.

Thứ hai, ứng dụng mô hình tương quan vào công tác quản lý tài nguyên. Chi cục Kiểm lâm tỉnh Hà Giang và Hạt Kiểm lâm Hoàng Su Phì cần tích hợp các phương trình hồi quy $Y = a \cdot X^b$ vào phần mềm kiểm kê rừng từ năm 2026, giúp giảm thiểu 80% chi phí điều tra ngoại nghiệp và rút ngắn 50% thời gian đo đếm hiện trường.

Thứ ba, xúc tiến dự án thương mại hóa tín chỉ cacbon rừng. Sở Nông nghiệp và Phát triển nông thôn tỉnh Hà Giang phối hợp với các tổ chức tư vấn quốc tế xây dựng hồ sơ tham gia thị trường cacbon tự nguyện (thị trường OTC/cơ chế REDD+) cho 2.400 ha rừng thông tại địa phương, tạo nguồn doanh thu bổ sung từ việc bán tín chỉ CO2 trong giai đoạn 2026 - 2030.

Thứ tư, tăng cường quản lý thảm mục và bảo vệ đất dốc. Ban Quản lý rừng phòng hộ phối hợp với chính quyền 25 xã và thị trấn tuyên truyền người dân chấm dứt việc thu gom thảm mục và đốt thực bì, giữ độ dày tầng vật rơi rụng từ 2 đến 3 tấn/ha nhằm bảo toàn 30 - 40 tấn C/ha trong tầng đất 0 - 30 cm, triển khai thường xuyên hàng năm trước mùa khô hanh.

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

Luận văn là nguồn tài liệu chuyên khảo giá trị cho 4 nhóm đối tượng chính:

  1. Cơ quan quản lý nhà nước về lâm nghiệp: Sở Nông nghiệp và PTNT, Chi cục Kiểm lâm các tỉnh miền núi phía Bắc sử dụng số liệu của 2.400 ha rừng thông tại Hoàng Su Phì để xây dựng chính sách chi trả dịch vụ môi trường rừng (PES) và định mức kỹ thuật lâm sinh địa phương.
  2. Ban Quản lý rừng và các công ty lâm nghiệp: Các đơn vị quản lý 15.502,73 ha rừng phòng hộ và 27.702,64 ha rừng sản xuất trên địa bàn áp dụng biểu mật độ từ 500 đến 1.600 cây/ha để lập kế hoạch tỉa thưa, nuôi dưỡng và phòng cháy chữa cháy rừng hiệu quả.
  3. Tổ chức phát triển dự án và chuyên gia biến đổi khí hậu: Các nhà phát triển dự án REDD+, CDM hoặc thị trường cacbon tự do sử dụng hệ số chuyển đổi sinh khối và hàm lượng cacbon thực nghiệm (36,9% - 43,36%) làm thông số đầu vào chính xác để xây dựng văn kiện dự án (PDD).
  4. Giảng viên, học viên sau đại học và sinh viên ngành Lâm học: Tham khảo khung phương pháp chặt giải tích cây mẫu, quy trình phân tích đất Walkley - Black và các thuật toán kiểm định hồi quy ANOVA trong SPSS phục vụ công tác đào tạo và nghiên cứu học thuật.

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

  1. Hàm lượng cacbon thực tế của Thông ba lá có điểm gì khác so với hướng dẫn mặc định của IPCC? Hàm lượng cacbon thực nghiệm trong các bộ phận Thông ba lá dao động từ 36,90% (vỏ) đến 43,36% (thân), thấp hơn đáng kể so với mức mặc định 50% do IPCC khuyến nghị. Việc sử dụng số liệu phân tích hóa nghiệm trực tiếp giúp loại bỏ sai số phóng đại khi tính toán tín chỉ CO2 cho rừng trồng địa phương.

  2. Trữ lượng cacbon của rừng Thông ba lá tại Hoàng Su Phì biến động ra sao theo độ tuổi? Tổng trữ lượng cacbon lâm phần tăng liên tục theo cấp tuổi, đạt từ 66,5 tấn C/ha ở tuổi 5 lên 144,3 tấn C/ha ở tuổi 26. Trong đó, tầng cây gỗ đứng trên mặt đất đóng vai trò chủ đạo, mở rộng tỷ trọng từ 38,4% lên 80,7% tổng trữ lượng của toàn bộ hệ sinh thái rừng.

  3. Tại sao tỷ lệ cacbon dưới mặt đất lại suy giảm tương đối so với trên mặt đất khi lâm phần già đi? Ở tuổi 5, tỷ lệ cacbon dưới mặt đất so với trên mặt đất đạt 160,66% do sinh khối đất chiếm ưu thế, nhưng con số này giảm xuống còn 23,92% ở tuổi 26. Nguyên nhân là do tầng thân gỗ tích lũy sinh khối vượt bậc (đạt 441,73 kg/cây), trong khi trữ lượng cacbon trong đất giữ mức ổn định từ 27,8 đến 41 tấn C/ha.

  4. Mô hình toán học nào ước tính sinh khối Thông ba lá chính xác nhất hiện nay? Phương trình dạng hàm lũy thừa $Y = a \cdot X^b$ (với $X$ là đường kính ngang ngực DBH) cho độ chính xác cao nhất với hệ số tương quan $r \ge 0,99$ và mức ý nghĩa thống kê $P < 0,05$. Mô hình này cho phép tính toán nhanh trữ lượng rừng mà không cần chặt hạ cây mẫu.

  5. Việc định lượng trữ lượng cacbon mang lại lợi ích kinh tế trực tiếp nào cho người trồng rừng? Với diện tích 2.400 ha rừng Thông ba lá tại Hoàng Su Phì, việc lượng giá bể chứa cacbon mở ra cơ hội thương mại hóa tín chỉ CO2 trên thị trường tự nguyện OTC. Nguồn thu từ dịch vụ môi trường rừng tạo thêm sinh kế bền vững cho người dân địa phương bên cạnh việc khai thác gỗ.

Kết luận

  • Xác định sinh khối khô cá thể Thông ba lá tăng trưởng mạnh từ 31,01 kg/cây (tuổi 5) lên 645,38 kg/cây (tuổi 26) với hàm lượng cacbon trong thân đạt 43,36%.
  • Định lượng tổng trữ lượng cacbon của lâm phần đạt từ 66,5 đến 144,3 tấn C/ha, trong đó tỷ trọng tích lũy trên mặt đất chiếm từ 38,4% đến 80,7%.
  • Thiết lập thành công hệ thống phương trình tương quan lũy thừa giữa đường kính ngang ngực (DBH) và sinh khối với độ tin cậy vượt trội ($r \ge 0,99$).
  • Đánh giá tiềm năng kinh tế to lớn của 2.400 ha rừng thông tại Hoàng Su Phì khi tham gia thị trường chuyển nhượng tín chỉ giảm phát thải CO2.
  • Cung cấp bộ dữ liệu thực nghiệm chuẩn xác phục vụ trực tiếp cho chiến lược quản lý rừng bền vững và kiểm kê khí nhà kính quốc gia.

Đóng góp lớn nhất của luận văn là thiết lập hệ thống cơ sở khoa học định lượng toàn diện về sinh khối và bể chứa cacbon của Thông ba lá tại vùng núi cao phía Bắc Việt Nam. Trong lộ trình 2026 - 2030, các địa phương cần đẩy mạnh số hóa mô hình tính toán và tích cực xúc tiến thị trường dịch vụ môi trường rừng. Các tổ chức lâm nghiệp và nhà nghiên cứu quan tâm hãy liên hệ các cơ quan chuyên môn để tiếp cận toàn văn công trình và ứng dụng bộ công cụ này vào thực tiễn quản lý rừng ngay hôm nay.