Tổng quan nghiên cứu
Rừng trồng đóng vai trò chiến lược trong việc cung ứng nguyên liệu gỗ công nghiệp toàn cầu, với tỷ lệ đóng góp tăng mạnh từ 5% vào năm 1960 lên 30% năm 2005 và ước tính đạt 75% vào năm 2050 theo các dự báo lâm nghiệp quốc tế. Tuy nhiên, biến đổi khí hậu toàn cầu đang đe dọa trực tiếp đến tính ổn định của hệ sinh thái rừng trồng; tại nhiều khu vực, các hiện tượng hạn hán kéo dài và nhiệt độ cực đoan có thể làm suy giảm tới 20% năng suất lâm phần. Tại Việt Nam, Keo tai tượng (Acacia mangium Willd.) là một trong những loài cây trồng rừng kinh tế chủ lực, cung cấp hàng triệu mét khối gỗ dăm và gỗ xẻ mỗi năm. Mặc dù vậy, công tác quản lý và dự báo sản lượng rừng trồng hiện nay chủ yếu vẫn dựa vào các mô hình kinh nghiệm truyền thống. Những mô hình tĩnh này bộc lộ nhược điểm lớn khi không thể phản ánh phản ứng sinh lý của cây rừng trước sự biến đổi của các yếu tố khí hậu, dinh dưỡng đất đai và các biện pháp kỹ thuật lâm sinh như mật độ hoặc tỉa thưa.
Nhằm giải quyết nút thắt này, luận văn thạc sĩ lâm học của tác giả Nguyễn Tuấn Linh dưới sự hướng dẫn khoa học của Tiến sĩ Phan Minh Sáng tại Trường Đại học Lâm nghiệp Việt Nam đã tiến hành đề tài mô hình hóa động thái năng suất, sinh khối và hấp thụ carbon của rừng Keo tai tượng bằng phần mềm động thái 3-PG. Mục tiêu trọng tâm của nghiên cứu là xác định năng suất, sinh khối và lượng carbon tích lũy thực tế của rừng Keo tai tượng ở các cấp tuổi từ 3 đến 17; đồng thời hiệu chỉnh, chuẩn hóa bộ tham số sinh lý - lập địa cho phần mềm 3-PG để kiểm tra khả năng ứng dụng công cụ này tại Việt Nam. Phạm vi nghiên cứu được thực hiện trên 66 lâm phần trải dài từ Bắc vào Nam qua 8 tỉnh trọng điểm gồm Tuyên Quang, Quảng Trị, Thừa Thiên Huế, Quảng Nam, Gia Lai, Lâm Đồng, Đồng Nai và Bình Dương trong giai đoạn 2009 - 2011. Nghiên cứu mang lại ý nghĩa học thuật và thực tiễn sâu sắc, cung cấp bộ công cụ mô phỏng đạt hệ số tương quan $R^2 = 0{,}8771$ giữa lý thuyết và thực nghiệm, mở ra hướng đi mới cho việc tối ưu hóa chu kỳ kinh doanh gỗ và định lượng giá trị dịch vụ môi trường rừng theo cơ chế chi trả carbon.
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 được xây dựng trên nền tảng lý thuyết sinh thái học hệ sinh thái, lý thuyết cân bằng carbon sinh học và nguyên lý mô hình hóa tăng trưởng rừng. Trong khoa học lâm nghiệp, các mô hình dự báo sinh trưởng được phân loại thành ba nhóm chính: mô hình thực nghiệm thống kê, mô hình động thái sinh lý và mô hình lai. Trong khi mô hình thực nghiệm phụ thuộc vào chuỗi số liệu quá khứ và không phản ánh được tác động của biến đổi môi trường, mô hình động thái 3-PG (Physiological Principles in Predicting Growth) do Landsberg và Waring phát triển từ năm 1997 hoạt động dựa trên các nguyên lý sinh lý học kết hợp với các tham số điều kiện tự nhiên.
Khung lý thuyết của mô hình 3-PG vận hành dựa trên cơ chế hấp thụ bức xạ quang hợp hoạt động (APAR) từ tổng lượng bức xạ mặt trời chiếu xuống lâm phần, sau đó chuyển hóa thành tổng năng suất sơ cấp (GPP) thông qua hiệu quả sử dụng ánh sáng của tán lá. Năng suất sơ cấp thô sau đó bị chiết giảm bởi các chỉ số giới hạn môi trường bao gồm: nhiệt độ bình quân tháng ($f_T$), độ hụt áp suất hơi bão hòa ($f_D$), lượng nước hữu hiệu trong đất ($f_\theta$), số ngày sương giá ($f_F$), độ phì đất ($f_N$) và độ tuổi lâm phần ($f_{\text{age}}$). Sản lượng quang hợp thuần (NPP) được ấn định tương đương tỷ lệ cố định khoảng 0,45 đến 0,47 so với GPP sau khi trừ tiêu hao hô hấp duy trì. Cuối cùng, NPP được phân bổ động thái đến các bể sinh khối lá, cành thân và rễ dựa trên điều kiện lập địa và độ tuổi. Lượng carbon tích lũy trong các bộ phận thực vật được quy đổi theo hướng dẫn tiêu chuẩn quốc tế của Ủy ban Liên chính phủ về Biến đổi Khí hậu (IPCC) với hệ số chuyển đổi là 0,5 lần sinh khối khô tuyệt đối.
Phương pháp nghiên cứu
Nghiên cứu sử dụng phương pháp kết hợp giữa kế thừa số liệu điều tra thứ cấp từ đề tài nghiên cứu cấp Bộ thuộc Viện Khoa học Lâm nghiệp Việt Nam giai đoạn 2010 - 2012 và thu thập dữ liệu sơ cấp ngoài hiện trường. Cỡ mẫu nghiên cứu bao gồm 66 lâm phần Keo tai tượng thuần loại từ 3 đến 17 tuổi, trong đó mỗi lâm phần được bố trí 3 ô tiêu chuẩn hình tròn ngẫu nhiên với dung lượng tối thiểu 30 cây mỗi ô nhằm đảm bảo tính đại diện thống kê.
Phương pháp chặt hạ cây mẫu giải tích được tiến hành trên 16 cây tiêu chuẩn đại diện cho các cấp đường kính từ 10 đến 30 cm để cân sinh khối tươi riêng biệt cho thân, cành, lá. Các mẫu thí nghiệm có khối lượng từ 50 đến 100 gam được sấy khô kiệt ở nhiệt độ 105 độ C đối với thân, cành và 75 độ C đối với lá đến khối lượng không đổi nhằm xác định hàm lượng ẩm và tỷ lệ sinh khối khô. Các phương trình tương quan sinh khối allometric được xây dựng bằng phương pháp hồi quy phi tuyến với đường kính ngang ngực ($D_{1.3}$). Dữ liệu khí tượng thu thập từ trạm khí tượng thủy văn khu vực bao gồm nhiệt độ, lượng mưa, bức xạ mặt trời, độ bốc thoát hơi nước. Năm nhóm đất chính theo phân loại FAO-UNESCO gồm Ferralic Acrisols, Dystric Cambisols, Ferric Acrisols, Xanthic Ferrasols và Haplic Acrisols được phân tích các chỉ tiêu hóa lý như pH, dung trọng, hàm lượng cát - sét, dung tích hấp thu cation (CEC), hàm lượng carbon tổng số và nitơ tổng số.
Kết quả nghiên cứu và thảo luận
Những phát hiện chính
Thứ nhất, nghiên cứu ghi nhận sự biến động rõ rệt về năng suất tăng trưởng bình quân hàng năm (MAI) của rừng Keo tai tượng giữa các vùng sinh thái và các cấp tuổi. Tại Tuyên Quang, MAI dao động ổn định từ 11,2 đến 19,6 $m^3$/ha/năm. Ngược lại, khu vực miền Trung và Nam Bộ có sự phân hóa mạnh mẽ: Thừa Thiên Huế dao động từ 5,5 đến 24,2 $m^3$/ha/năm, Lâm Đồng từ 13,9 đến 29,9 $m^3$/ha/năm, Đồng Nai từ 16,1 đến 32,4 $m^3$/ha/năm. Tỉnh Bình Dương ghi nhận mức tăng trưởng trung bình cao nhất đạt 30,0 $m^3$/ha/năm, trong khi Quảng Nam đạt mức thấp nhất với 10,2 $m^3$/ha/năm, thể hiện mức chênh lệch năng suất lên tới gần 300% giữa các lập địa tối ưu và lập địa nghèo kiệt.
Thứ hai, các phương trình allometric thực nghiệm về sinh khối bộ phận cây cá lẻ đạt độ tin cậy thống kê vượt trội. Cụ thể: sinh khối thân cây tuân theo phương trình $W_{\text{thân}} = 0{,}0573 \times D^{2{,}378}$ với sai số tuyệt đối trung bình chỉ 5,0%; sinh khối lá $W_{\text{lá}} = 0{,}0709 \times D^{1{,}468}$ với sai số 11,0%; sinh khối cành $W_{\text{cành}} = 0{,}0116 \times D^{2{,}117}$ với sai số 9,0% ($P < 0{,}001$). Tổng sinh khối trên mặt đất tích lũy cao nhất đạt trên 250 tấn khô/ha ở các lâm phần 16 - 17 tuổi, tương ứng với lượng hấp thụ carbon đạt trên 125 tấn C/ha.
Thứ ba, nghiên cứu đã chuẩn hóa thành công phương pháp gán tham số độ phì đất ($FR$) và lượng nước hữu hiệu tối đa ($ASW_{\text{max}}$) cho mô hình 3-PG tại Việt Nam. Giá trị $FR$ được lượng hóa lần lượt là 0,9 cho đất Dystric Cambisols và Xanthic Ferrasols; 0,7 cho Ferric Acrisols; và 0,3 cho Ferralic Acrisols và Haplic Acrisols. Giá trị $ASW_{\text{max}}$ tương ứng là 235 mm, 155 mm, 185 mm, 135 mm và 115 mm. Kết quả chạy mô phỏng trên phần mềm 3-PG đối chiếu với số liệu đo đếm thực tế đạt hệ số tương quan tuyến tính $R^2 = 0{,}8771$, khẳng định độ tương thích trên 70% các lâm phần nghiên cứu mà không cần qua quy trình hiệu chỉnh thử sai phức tạp.
Thảo luận kết quả
Sự sai khác về sinh trưởng giữa các vùng nghiên cứu bắt nguồn từ sự tương tác phức tạp giữa chế độ nhiệt - ẩm và đặc tính lý hóa của đất. Đất Dystric Cambisols và Xanthic Ferrasols tại Bình Dương hay Đồng Nai có tầng dày, độ phì cao và khả năng giữ nước hữu hiệu vượt trội ($ASW > 155\text{ mm}$), giúp cây duy trì chỉ số diện tích lá (LAI) cao trong suốt mùa khô kéo dài 4 - 6 tháng. Ngược lại, tại Quảng Nam và các vùng đất cát, đất xói mòn trơ sỏi đá có hàm lượng lân hữu dụng thấp (dưới 0,2 ppm), chỉ số hụt áp suất hơi bão hòa cao vào mùa khô làm đóng khí khổng, ức chế quang hợp và dẫn tới sinh trưởng suy giảm.
Dữ liệu so sánh giữa thực nghiệm và mô phỏng được thể hiện trực quan qua các biểu đồ phân bố sai số và biểu đồ phân tán hồi quy tuyến tính giữa MAI mô phỏng và MAI thực tế. So với các nghiên cứu của Landsberg tại Úc hay Almeida tại Nam Mỹ trên cây Bạch đàn, việc ứng dụng 3-PG cho loài Keo tai tượng tại Việt Nam cho thấy sai số mô phỏng trữ lượng được kéo giảm đáng kể khi thay thế bộ tham số mặc định của phần mềm bằng bộ tham số thực nghiệm được đề tài thiết lập. Sự phù hợp này chứng minh rằng việc tích hợp các hàm giới hạn sinh lý khí hậu vào dự báo sản lượng rừng tại Việt Nam là hoàn toàn khả thi và có độ tin cậy vượt trội so với các biểu sản lượng tĩnh truyền thống.
Đề xuất và khuyến nghị
Thứ nhất, chuyển giao và ứng dụng phần mềm 3-PG vào hệ thống quản lý của các Công ty Lâm nghiệp và Ban Quản lý rừng phòng hộ trên cả nước. Chủ rừng cần thiết lập cơ sở dữ liệu khí tượng và bản đồ đất chi tiết để vận hành mô phỏng, đặt mục tiêu nâng cao độ chính xác trong dự báo sản lượng khai thác lên trên 85% trong giai đoạn 2026 - 2030.
Thứ hai, áp dụng các biện pháp kỹ thuật thâm canh lâm sinh theo từng cấp lập địa được mô hình hóa. Tại các lập địa có độ phì thấp ($FR = 0{,}3$) như Haplic Acrisols, các đơn vị quản lý cần bổ sung phân bón lót chứa lân và duy trì thảm phủ hữu cơ nhằm tăng khả năng giữ ẩm, phấn đấu nâng năng suất rừng trồng tại các vùng này từ mức 10 - 12 $m^3$/ha/năm lên tối thiểu 18 $m^3$/ha/năm trong chu kỳ 7 năm.
Thứ ba, điều chỉnh mật độ trồng ban đầu và kỹ thuật tỉa thưa dựa trên dự báo động thái sinh khối của mô hình. Khuyến nghị các hộ gia đình và doanh nghiệp trồng rừng duy trì mật độ trồng ban đầu từ 1.330 đến 1.660 cây/ha và tiến hành tỉa thưa ở giai đoạn tuổi 4 - 5 để tập trung năng suất sinh khối vào thân cây gỗ lớn, hướng tới mục tiêu gia tăng giá trị kinh tế thêm 25 - 30% so với kinh doanh gỗ dăm nhỏ.
Thứ tư, xây dựng hồ sơ kỹ thuật carbon rừng cho các dự án trồng rừng Keo tai tượng nhằm tiếp cận thị trường tín chỉ carbon tự nguyện. Sở Nông nghiệp và Phát triển Nông thôn các tỉnh cần chủ trì phối hợp với các viện nghiên cứu sử dụng thuật toán 3-PG để lập báo cáo giám sát hấp thụ carbon định kỳ hàng năm, tạo nguồn thu bổ sung ước đạt 50 - 100 USD/ha/năm từ dịch vụ chi trả môi trường rừng.
Đối tượng nên tham khảo luận văn
Thứ nhất, các nhà nghiên cứu lâm nghiệp và học viên cao học chuyên ngành Lâm sinh, Quản lý tài nguyên rừng: Tài liệu cung cấp cơ sở phương pháp luận chi tiết, hệ thống phương trình sinh khối allometric và quy trình chuẩn hóa tham số sinh lý cho mô hình động thái 3-PG.
Thứ hai, các kỹ sư lâm nghiệp, nhà quy hoạch lâm nghiệp tại các Viện điều tra quy hoạch rừng và Chi cục Kiểm lâm: Cung cấp công cụ khoa học để dự báo trữ lượng gỗ, đánh giá độ thích nghi loài cây trồng theo từng tiểu vùng khí hậu và lập kế hoạch điều chế rừng dài hạn.
Thứ ba, lãnh đạo các Công ty Trách nhiệm Hữu hạn Một thành viên Lâm nghiệp và chủ rừng quy mô lớn: Luận văn hỗ trợ xây dựng kịch bản kinh doanh tối ưu, lựa chọn thời điểm luân kỳ khai thác hợp lý và đánh giá rủi ro hạn hán tác động đến năng suất lâm phần.
Thứ tư, các chuyên gia tư vấn phát triển dự án biến đổi khí hậu và thị trường carbon (REDD+): Nguồn dữ liệu thực nghiệm chuẩn xác về sinh khối khô và hệ số hấp thụ carbon trên 8 tỉnh thành giúp thiết lập đường phát thải tham chiếu và tính toán tín chỉ carbon chính xác cho rừng trồng thương mại.
Câu hỏi thường gặp
Mô hình động thái 3-PG có ưu điểm gì vượt trội so với các biểu sản lượng truyền thống?
Mô hình 3-PG tích hợp dữ liệu khí tượng tháng và đặc tính đất đai để mô phỏng quá trình sinh lý quang hợp và phân bổ sinh khối, giúp dự báo chính xác năng suất khi khí hậu biến động, trong khi biểu truyền thống chỉ phản ánh số liệu quá khứ tĩnh.
Hệ số chuyển đổi sinh khối sang lượng hấp thụ carbon của rừng Keo tai tượng là bao nhiêu?
Theo nghiên cứu và tiêu chuẩn của IPCC, hệ số chuyển đổi được xác định là 0,5 lần sinh khối khô tuyệt đối; nghĩa là cứ 1 tấn sinh khối khô trên mặt đất sẽ tích lũy cố định 0,5 tấn carbon nguyên tố tương đương 1,83 tấn CO2.
Tại sao năng suất rừng Keo tai tượng tại Quảng Nam thấp hơn đáng kể so với Đồng Nai và Bình Dương?
Nguyên nhân chủ yếu do chất lượng đất tại Quảng Nam nghèo dinh dưỡng ($FR = 0{,}3$), khả năng giữ nước hữu hiệu thấp ($ASW = 115\text{ mm}$), kết hợp thời tiết khắc nghiệt mùa khô làm cây bị hạn sinh lý, giảm năng suất quang hợp xuống 10,2 $m^3$/ha/năm so với 30 - 32 $m^3$/ha/năm ở Nam Bộ.
Làm thế nào để xác định tham số độ phì đất cho mô hình 3-PG mà không cần thử sai?
Luận văn đã lượng hóa độ phì theo 5 nhóm đất FAO-UNESCO: gán $FR = 0{,}9$ cho Dystric Cambisols và Xanthic Ferrasols, $FR = 0{,}7$ cho Ferric Acrisols, và $FR = 0{,}3$ cho Ferralic Acrisols, Haplic Acrisols, đạt độ chính xác mô phỏng trên 70%.
Mô hình 3-PG có thể áp dụng cho các loài cây rừng trồng khác tại Việt Nam không?
Hoàn toàn có thể áp dụng, với điều kiện phải hiệu chỉnh lại các tham số sinh lý đặc thù của loài như hệ số dập tắt ánh sáng, hiệu quả lượng tử ánh sáng của tán lá và phương trình phân bổ sinh khối cho từng loài cụ thể.
Kết luận
- Luận văn đã xác định toàn diện động thái năng suất và sinh khối của rừng Keo tai tượng từ 3 đến 17 tuổi trên 66 lâm phần tại 8 tỉnh thành, với năng suất trung bình dao động từ 10,2 đến 32,4 $m^3$/ha/năm.
- Xây dựng thành công hệ phương trình sinh khối tương quan allometric cho các bộ phận thân, cành, lá có độ chính xác cao với sai số dao động từ 5,0% đến 11,0%.
- Tiên phong trong việc hiệu chỉnh và gán giá trị tham số độ phì đất và độ ẩm hữu hiệu cho 5 nhóm đất theo chuẩn FAO-UNESCO, đạt hệ số tương quan mô phỏng $R^2 = 0{,}8771$.
- Cung cấp cơ sở khoa học định lượng phục vụ công tác quản lý rừng bền vững, chuyển đổi rừng gỗ nhỏ sang gỗ lớn và phát triển các dự án tín chỉ carbon lâm nghiệp.
- Mở ra giai đoạn tiếp theo cần tập trung số hóa công cụ 3-PG trên nền tảng web GIS để hỗ trợ người trồng rừng tra cứu năng suất trực tuyến từ năm 2026.
Các nhà quản lý lâm nghiệp, viện nghiên cứu và chủ rừng hãy chủ động ứng dụng mô hình 3-PG vào thực tiễn sản xuất để tối ưu hóa năng suất rừng trồng và gia tăng nguồn thu từ thị trường tín chỉ carbon ngay hôm nay!