Tổng quan về luận án

Nghiên cứu sinh khối và định lượng trữ lượng carbon rừng đóng vai trò hạt nhân trong các nỗ lực toàn cầu nhằm giảm thiểu biến đổi khí hậu theo Công ước khung của Liên Hiệp Quốc về Biến đổi Khí hậu (UNFCCC) và cơ chế Chi trả Giảm phát thải từ Mất rừng và Suy thoái rừng (UN-REDD+). Trong bối cảnh đó, hệ sinh thái Rừng khộp (Tropical Deciduous Dry Dipterocarp Forest - DF) – kiểu rừng khô rụng lá nhiệt đới đặc thù chỉ phân bố tại châu Á và chiếm khoảng 650.000 ha tại vùng Tây Nguyên và Đông Nam Bộ Việt Nam – sở hữu tiềm năng cô lập carbon vượt trội nhưng lại chưa có một hệ thống biểu phương trình sinh khối tương thích, chuẩn hóa và được kiểm định sai số khách quan. Luận án tiến sĩ lâm nghiệp của nghiên cứu sinh Nguyễn Thị Tình với đề tài "Thiết lập và thẩm định chéo hệ thống mô hình ước tính sinh khối trên mặt đất cây rừng khộp ở Việt Nam" (chuyên ngành Lâm sinh, mã số 9 62 02 05, thực hiện tại Trường Đại học Tây Nguyên dưới sự hướng dẫn khoa học của GS. Bảo Huy) là công trình nghiên cứu tiên phong giải quyết trọn vẹn khoảng trống khoa học này.

Khoảng trống nghiên cứu (research gap) then chốt mà luận án nhận diện xuất phát từ việc các phương trình sinh khối toàn cầu dạng tổng quát (pan-tropical generic allometric equations) như nghiên cứu của Chave và cộng sự (2005, 2014) thường dẫn đến sai số hệ thống nghiêm trọng khi áp dụng cho các lâm phần địa phương. Cụ thể, trích dẫn trực tiếp từ luận án chỉ ra: "Phương trình sinh khối chung cho vùng nhiệt đới như vậy có thể dẫn đến sai số hệ thống lên đến 400%, do vậy các mô hình tại địa phương có thể là một thay thế tốt hơn và dự kiến sẽ cung cấp độ tin cậy cao hơn so với phương trình chung (Chave và ctv, 2014)". Đồng thời, các nghiên cứu trước đây tại Việt Nam và Đông Nam Á (Basuki và ctv, 2009; Huy và ctv, 2016c; Kralicek và ctv, 2017) hoặc có dung lượng mẫu cây chặt hạ hạn chế ở cấp chi/loài và cây kích thước lớn, hoặc chỉ thiết lập các phương trình độc lập rời rạc, vi phạm tính cộng gộp sinh học tự nhiên giữa các bộ phận cây.

Để giải quyết triệt để vấn đề trên, luận án xác lập hệ thống 4 câu hỏi nghiên cứu (Research Questions - RQ) và 4 giả thuyết khoa học (Hypotheses - H) tương ứng:

  • RQ1: Dạng mô hình phi tuyến tính có trọng số theo phương pháp Hợp lý cực đại (Weighted Nonlinear fit by Maximum Likelihood) có triệt tiêu được hiện tượng phương sai thay đổi (heteroscedasticity) và cho độ chính xác cao hơn dạng hàm Power tuyến tính hóa logarit hay không?
  • RQ2: Hệ thống phương trình hồi quy dường như phi tương quan (Seemingly Unrelated Regression - SUR) có đảm bảo tính cộng gộp sinh học và giảm thiểu sai số chéo giữa các bộ phận thân ($B_{st}$), cành ($B_{br}$), lá ($B_{le}$), vỏ ($B_{ba}$) và tổng sinh khối trên mặt đất ($AGB$) so với phương pháp thiết lập độc lập không?
  • RQ3: Việc xây dựng mô hình sinh khối theo phân cấp thực vật học (họ Dipterocarpaceae, chi Dipterocarpus, Shorea và 4 loài ưu thế) với biến số đo đếm đơn giản (đường kính ngang ngực $D$) có đạt độ tin cậy vượt trội hơn mô hình đa loài phức tạp đòi hỏi nhiều biến số ($D, H, WD$) hay không?
  • RQ4: Các nhân tố sinh thái môi trường (lượng mưa $P$, độ cao $Altitude$) và cấu trúc lâm phần (tổng tiết diện ngang $BA$, mật độ $N$) tương tác như thế nào lên tham số ước tính sinh khối thông qua hàm điều chỉnh (Modifier function)?

Khung lý thuyết của luận án tích hợp Lý thuyết Tương quan Sinh trắc học (Allometric Scaling Theory), Kinh tế lượng hệ thống SUR (Parresol, 2001), và Hướng dẫn Thực hành Tốt của IPCC (2003, 2006) về Đo đạc, Báo cáo và Thẩm định (MRV Tier 3). Phạm vi nghiên cứu bao phủ 28 ô tiêu chuẩn định vị diện tích 0,25 ha ($50 \times 50\text{ m}$) tại Tây Nguyên (CH) và 01 ô tiêu chuẩn diện tích 1,0 ha ($100 \times 100\text{ m}$) tại Đông Nam Bộ (SE). Dữ liệu lâm phần ghi nhận mật độ lâm phần dao động từ 228 đến 1291 cây/ha (đối với cây có $D \ge 5\text{ cm}$), tổng tiết diện ngang $BA$ dao động từ 3,8 đến 23,4 $\text{m}^2/\text{ha}$, lượng mưa trung bình năm biến thiên từ 1003 đến 1600 mm/năm, và độ cao tuyệt đối từ 171 đến 417 m. Toàn bộ hệ thống mô hình được xây dựng từ tập dữ liệu cây mẫu chặt hạ giải phẫu sinh khối trực tiếp kế thừa và bổ sung từ chương trình UN-REDD Việt Nam (với hơn 222 cây mẫu), tạo ra bước đột phá định lượng cho công tác kiểm kê carbon rừng cấp quốc gia.

Literature Review và Positioning

Tổng quan tài liệu học thuật quốc tế về lượng hóa sinh khối và carbon rừng nhiệt đới ghi nhận sự phân nhánh rõ nét thành ba luồng nghiên cứu chủ đạo (major research streams):

Luồng nghiên cứu 1 - Các mô hình toàn cầu/vùng nhiệt đới tổng quát: Khởi xướng từ các công trình kinh điển của Brown và cộng sự (1989, 1997) và được phát triển mạnh mẽ bởi Chave và cộng sự (2005, 2014). Nhóm tác giả này sử dụng tập dữ liệu toàn cầu đa loài (pan-tropical multi-species datasets với 2410 cây mẫu) để đề xuất dạng hàm Power tích hợp ba biến số: đường kính ngang ngực ($D$), chiều cao vút ngọn ($H$), và khối lượng thể tích gỗ ($WD$). Mặc dù mang tính khái quát cao, các mô hình này bộc lộ nhược điểm lớn khi áp dụng vào các hệ sinh thái rụng lá khô đặc thù do hình thái tán và tỷ lệ phân chia sinh khối thân-cành-lá sai khác hoàn toàn so với rừng mưa thường xanh.

Luồng nghiên cứu 2 - Các mô hình khu vực cho rừng họ Dầu Đông Nam Á và Nam Á: Điển hình là các công trình của Basuki và cộng sự (2009) tại Đông Kalimantan (Indonesia), Hanpattanakit và cộng sự (2016) tại Thái Lan, Panda và cộng sự (2020) tại Ấn Độ, và Laumonier và cộng sự (2010) tại Sumatra. Các nghiên cứu này chỉ ra sinh khối trên mặt đất của rừng khộp dao động từ 297 đến 386 tấn/ha, cao hơn đáng kể so với giá trị mặc định của IPCC (2006). Basuki và cộng sự (2009) khẳng định rằng việc phân nhóm mô hình theo chi thực vật (Dipterocarpus, Shorea, Hopea) mang lại độ chính xác vượt bậc so với mô hình gộp chung mọi loài.

Luồng nghiên cứu 3 - Phát triển hệ thống phương trình cộng gộp và thẩm định sai số: Đại diện bởi Parresol (2001), Picard và cộng sự (2012, 2015), Poudel & Temesgen (2016), Kralicek và cộng sự (2017), và Huy và cộng sự (2016a,b,c, 2019). Luồng nghiên cứu này chỉ rõ rằng việc lập mô hình độc lập cho từng bộ phận dẫn đến tổng sinh khối các bộ phận không bằng sinh khối toàn cây, đồng thời thúc đẩy ứng dụng kỹ thuật hồi quy dường như phi tương quan (SUR) và phương pháp Thẩm định chéo (Cross-Validation) nhằm tránh hiện tượng quá khớp (overfitting).

Tranh luận học thuật cốt lõi trong y văn tập trung vào hai quan điểm đối lập: (1) Tranh luận giữa tuyến tính hóa logarit truyền thống (Brown, 1997; Chave et al., 2005) và hồi quy phi tuyến trực tiếp có trọng số (Weighted Nonlinear Regression) (Davidian & Giltinan, 1995; Picard et al., 2012). Việc biến đổi logarit $\ln(y) = \ln(a) + b\ln(x)$ làm sai lệch cấu trúc sai số gốc và đòi hỏi hệ số hiệu chỉnh chuyển đổi ngược (Correction Factor - CF), trong khi hồi quy phi tuyến trực tiếp kết hợp trọng số $w$ tối ưu hóa hàm hợp lý cực đại xử lý triệt để hiện tượng phương sai dị sai. (2) Tranh luận về độ phức tạp của biến đầu vào: liệu một mô hình phức tạp 3 biến ($D, H, WD$) có thực sự tốt hơn mô hình thực vật học đơn biến ($D$) khi xét đến sai số đo đếm chiều cao ($H$) và sự bất định của mật độ gỗ ($WD$) trên hiện trường điều tra.

Luận án định vị chính xác vào giao điểm của các khoảng trống học thuật nêu trên: là công trình đầu tiên tại khu vực hạ lưu sông Mê Kông thiết lập hoàn chỉnh hệ thống mô hình sinh khối cộng gộp SUR theo cấu trúc phân loại thực vật 3 cấp (Họ - Chi - Loài), kết hợp đồng thời hàm điều chỉnh sinh thái môi trường và quy trình thẩm định chéo nghiêm ngặt K-Fold, Monte Carlo, LOOCV.

Đóng góp lý thuyết và khung phân tích

Đóng góp cho lý thuyết

Luận án đóng góp trực tiếp vào việc mở rộng và tinh chỉnh Lý thuyết Tương quan Sinh trắc học (Allometric Scaling Theory của West, Brown & Enquist, 1997) và Hệ thống Sinh khối Cộng gộp (Additive Biomass Systems của Parresol, 2001) trên các khía cạnh:

  1. Mở rộng lý thuyết phân chia sinh khối cơ quan thực vật: Chứng minh rằng sự phân bổ sinh khối giữa thân ($B_{st}$), cành ($B_{br}$), lá ($B_{le}$) và vỏ ($B_{ba}$) của các loài cây họ Dầu chịu sự kiểm soát chặt chẽ của cơ chế thích ứng với lửa rừng và khô hạn định kỳ. Vỏ cây chiếm tỷ lệ sinh khối đáng kể ($B_{ba}$) đóng vai trò cách nhiệt sinh học. Hệ thống mô hình SUR đã mô hình hóa đồng thời ma trận hiệp phương sai sai số giữa các cơ quan, phản ánh chính xác cấu trúc sinh lý học thực vật mà các mô hình độc lập bỏ qua.
  2. Xác lập mô hình tương tác môi trường - sinh trắc: Luận án phát triển khung lý thuyết điều chỉnh sinh thái: $$AGB = \text{AVERAGE} \times \text{MODIFIER}$$ Trong đó, thành phần MODIFIER tích hợp trực tiếp ảnh hưởng của lượng mưa ($P$), độ cao ($Altitude$), và diện tích tiết diện ngang ($BA$), chứng minh rằng tham số allometric không phải là hằng số cố định mà biến thiên có quy luật theo gradient sinh thái cảnh quan.
  3. Chuyển dịch hệ hình (Paradigm shift): Chuyển từ hệ hình ước tính sinh khối dựa trên các phương trình hồi quy đơn lẻ, phi cộng gộp sang hệ hình Hệ thống mô hình sinh khối đồng thời, ràng buộc sinh học và kiểm định chéo khách quan.

Khung phân tích độc đáo

Khung phân tích của luận án tích hợp đồng bộ ba trụ cột lý thuyết và phương pháp luận:

  • Trụ cột 1 - Ràng buộc cộng gộp SUR: Thiết lập hệ phương trình đồng thời dạng: $$\begin{cases} B_{st} = f_1(X) + \varepsilon_1 \ B_{br} = f_2(X) + \varepsilon_2 \ B_{le} = f_3(X) + \varepsilon_3 \ B_{ba} = f_4(X) + \varepsilon_4 \ AGB = f_1(X) + f_2(X) + f_3(X) + f_4(X) + \varepsilon_5 \end{cases}$$ Với cấu trúc ma trận hiệp phương sai sai số $\boldsymbol{\Sigma}$ được ước lượng đồng thời, đảm bảo $\sum B_i = AGB$ trên mọi kích thước cây.
  • Trụ cột 2 - Chuẩn hóa sai số qua Chỉ số Furnival (Furnival's Index - FI): Sử dụng chỉ số FI (Furnival, 1961; Jayaraman, 1999) để so sánh khách quan giữa phương trình logarit tuyến tính và phi tuyến trực tiếp: $$FI = \left[ \prod_{i=1}^n \frac{1}{y_i} \right]^{-1/n} \times RMSE_{\ln}$$
  • Trụ cột 3 - Điều kiện biên (Boundary conditions): Khung phân tích được xác lập giới hạn cho hệ sinh thái rừng nhiệt đới rụng lá khô ưu thế họ Dầu, cấp kính $D \ge 5\text{ cm}$, sinh trưởng trên đất phong hóa từ đá mẹ bazan/xương xẩu, chịu áp lực cháy rừng mùa khô và ngập úng mùa mưa.

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

Thiết kế nghiên cứu

Luận án tuân thủ thế giới quan thực chứng hậu hiện đại (Post-positivism / Critical Realism), kết hợp giữa định lượng giải phẫu thực nghiệm và tối ưu hóa toán học thống kê. Thiết kế nghiên cứu đa tầng bậc (Multi-level design) được cấu trúc chặt chẽ:

  • Cấp độ lâm phần/vùng sinh thái: 28 ô mẫu ($50 \times 50\text{ m}$) tại Tây Nguyên và 01 ô mẫu ($100 \times 100\text{ m}$) tại Đông Nam Bộ.
  • Cấp độ phân loại thực vật: Toàn bộ loài gộp chung (Mixed species) $\rightarrow$ Họ Dipterocarpaceae $\rightarrow$ Chi Dipterocarpus và Chi Shorea $\rightarrow$ 4 loài ưu thế: Dầu đồng (Dipterocarpus tuberculatus Roxb.), Dầu trà beng (Dipterocarpus obtusifolius Teijsm.), Cà chít (Shorea obtusa Wall. Ex Blume), và Cẩm liên (Shorea siamensis Miq.).
  • Cấp độ cơ quan sinh khối: Phân tách thành 4 bộ phận cấu thành ($B_{st}, B_{br}, B_{le}, B_{ba}$) và tổng thể $AGB$.

Quy trình nghiên cứu rigorous

Quy trình thu thập dữ liệu giải phẫu sinh khối cây mẫu tuân thủ nghiêm ngặt các tiêu chuẩn quốc tế:

  1. Rút mẫu chặt hạ (Destructive sampling): Cây mẫu được chọn đại diện theo phân bố số cây theo cấp kính ($N/D$) của lâm phần tự nhiên nhằm đảm bảo tính đại diện cho toàn bộ dải đường kính.
  2. Xử lý mẫu sinh khối tươi và khô: Toàn bộ thân, cành, lá, vỏ của từng cây mẫu được cân tươi trực tiếp tại hiện trường bằng cân điện tử chuyên dụng. Các mẫu phụ đại diện cho từng bộ phận được thu thập, niêm phong và sấy tại phòng thí nghiệm ở nhiệt độ $80^\circ\text{C} - 105^\circ\text{C}$ cho đến khi khối lượng không đổi (thời gian sấy chuẩn hóa 48 giờ) theo giao thức của Silva (2011).
  3. Quy đổi trữ lượng Carbon và $CO_2$ tương đương: Luận án trích dẫn và áp dụng chính xác công thức định lượng tiêu chuẩn quốc tế:

    "Lượng carbon tích lũy và $CO_2$ tương đương cây rừng hấp thụ được tính toán thông qua hệ số chuyển đổi như sau: Carbon = 0,47 × Sinh khối (IPCC, 2006, Bảo Huy, 2013); $CO_2$ tương đương = 3,67 × Carbon"

  4. Tam giác đạc và kiểm soát độ tin cậy: Sử dụng 4 hệ số đánh giá sai số độc lập: Phần trăm sai lệch (Bias %), Sai số trung phương trung bình phần trăm (RMSE %), Sai số tuyệt đối trung bình phần trăm (MAPE %), và Hệ số xác định hiệu chỉnh ($R^2_{adj}$).

Data và phân tích

Đặc tính mẫu thống kê thể hiện độ biến động tự nhiên sâu sắc của hệ sinh thái rừng khộp:

  • Mật độ lâm phần ($N$): 228 – 1291 cây/ha (trung bình 685 cây/ha).
  • Tổng tiết diện ngang ($BA$): 3,8 – 23,4 $\text{m}^2/\text{ha}$.
  • Lượng mưa trung bình năm ($P$): 1003 – 1600 mm.
  • Nhiệt độ trung bình năm ($T$): $25,3^\circ\text{C} - 25,5^\circ\text{C}$.
  • Độ dốc địa hình: $0^\circ - 8^\circ$ (trung bình $2,0^\circ$).

Kỹ thuật phân tích dữ liệu nâng cao được thực thi trên hai nền tảng phần mềm thống kê chuyên sâu:

  • Ngôn ngữ R (các gói nlme, minpack.lm): Ước lượng mô hình phi tuyến có trọng số theo phương pháp Hợp lý cực đại (Weighted Nonlinear fit by Maximum Likelihood). Trọng số $w = 1/D^k$ hoặc $w = 1/(D^2H \cdot WD)^k$ được tích hợp để vô hiệu hóa hiện tượng dị sai (heteroscedasticity).
  • Phần mềm SAS (thủ tục PROC MODEL): Thiết lập hệ thống mô hình đồng thời SUR dạng phi tuyến có trọng số (Weighted Nonlinear SUR).
  • Thẩm định chéo đa giao thức (Cross-Validation Protocols): So sánh 4 phương pháp kiểm định sai số: Dữ liệu độc lập (tách 80% training / 20% validation), Thẩm định chéo Leave-One-Out (LOOCV), Thẩm định chéo K-Fold ($K = 5 - 10$), và Thẩm định chéo lặp lại Monte Carlo ($R = 500$ lần lặp). Kết quả chứng minh K-Fold và Monte Carlo cung cấp ước lượng sai số khách quan nhất mà không làm tổn thất dung lượng mẫu quý giá.

Phát hiện đột phá và implications

Những phát hiện then chốt

  1. Hồi quy phi tuyến Maximum Likelihood vượt trội hồi quy tuyến tính Logarit: Đánh giá qua chỉ số Furnival (FI) khẳng định dạng hàm phi tuyến ước lượng bằng Maximum Likelihood có trọng số đạt giá trị FI thấp hơn rõ rệt so với mô hình log-log biến đổi tuyến tính, loại bỏ hoàn toàn sai số chệch khi chuyển đổi ngược hàm số mũ.
  2. Hệ thống SUR loại trừ hoàn toàn sai số phi cộng gộp: Phương pháp SUR thiết lập thành công hệ thống phương trình các bộ phận ($B_{st}, B_{br}, B_{le}, B_{ba}$) và $AGB$. Đồ thị sai số Bland-Altman với độ tin cậy 95% chứng minh mô hình SUR duy trì tính cộng hưởng tuyệt đối và thu hẹp dải tin cậy sai số so với mô hình độc lập.
  3. Phát hiện nghịch lý - Đột phá về cấp độ phân loại thực vật: Luận án khám phá ra rằng mô hình sinh khối thiết lập riêng cho chi ưu thế (Dipterocarpus và Shorea) chỉ sử dụng duy nhất một biến số đường kính ngang ngực ($D$) đạt độ chính xác và độ ổn định thống kê (RMSE% và MAPE% thấp hơn) tương đương hoặc vượt trội so với mô hình chung đa loài đòi hỏi phối hợp cả ba biến phức tạp ($D, H, WD$).
  4. Tác động định lượng của nhân tố môi trường và lâm phần: Khi đưa các biến sinh thái vào mô hình thông qua hàm điều chỉnh MODIFIER, mô hình dạng: $$AGB = a \times D^b \times WD^d \times \exp\left(b_1 \cdot Altitude + b_2(P - 1502) + b_3(BA - 12,4)\right)$$ cho thấy sinh khối cá thể có tương quan nghịch có ý nghĩa thống kê ($p < 0,05$) với mật độ/tiết diện ngang lâm phần (do cạnh tranh không gian dinh dưỡng) và biến thiên rõ rệt theo lượng mưa mùa khô.
  5. Tính chuyển giao quốc tế của mô hình cấp chi: So sánh chéo mô hình cấp chi của luận án với mô hình của Basuki và cộng sự (2009) tại Indonesia cho thấy không có sự sai biệt thống kê đáng kể, mở ra khả năng áp dụng thống nhất mô hình cấp chi cho toàn bộ sinh cảnh rừng khộp Đông Nam Á.

Implications đa chiều

  • Về mặt lý thuyết: Củng cố sâu sắc các nguyên lý phân bổ năng lượng và sinh khối trong điều kiện lập địa khắc nghiệt; làm rõ mối quan hệ giữa đặc tính giải phẫu gỗ, độ dày vỏ chịu lửa và khả năng tích lũy carbon của thực vật thân gỗ họ Dầu.
  • Về mặt phương pháp luận: Chuẩn hóa quy trình 6 bước từ thu thập mẫu chặt hạ, xử lý dị sai, ước lượng phi tuyến SUR đến thẩm định chéo K-Fold/Monte Carlo. Đây là cẩm nang phương pháp luận mẫu mực cho các nghiên cứu sinh trắc học rừng nhiệt đới.
  • Về mặt thực tiễn và quản lý: Việc chứng minh mô hình cấp chi chỉ cần biến số đơn giản $D$ giúp đơn giản hóa triệt để công tác điều tra kiểm kê carbon tại hiện trường. Cán bộ kiểm lâm, ban quản lý rừng phòng hộ/đặc dụng và cộng đồng địa phương có thể dễ dàng đo đếm bằng thước dây/thước kẹp kính mà không cần các thiết bị đo chiều cao đắt tiền hay phân tích mật độ gỗ phức tạp.
  • Về mặt chính sách: Cung cấp trực tiếp bộ hệ số phát thải (Emission Factors - Tier 3) chuẩn xác phục vụ Hệ thống Đo đạc, Báo cáo và Thẩm định (MRV) của Chương trình Quốc gia UN-REDD+ Việt Nam, tạo cơ sở pháp lý và kỹ thuật vững chắc để thương mại hóa tín chỉ carbon trên thị trường quốc tế theo tiêu chuẩn Verified Carbon Standard (VCS) và ART-TREES.

Limitations và Future Research

Luận án thẳng thắn chỉ ra các giới hạn nghiên cứu (limitations) và điều kiện biên:

  • Giới hạn về bể chứa sinh khối dưới mặt đất (BGB): Do chi phí khai quật hệ rễ cây gỗ lớn cực kỳ tốn kém và phá hủy sinh thái nghiêm trọng, luận án tập trung sâu vào sinh khối trên mặt đất ($AGB$), trong khi dữ liệu sinh khối dưới mặt đất rễ cây ($BGB$) kế thừa từ các công bố chuyên biệt của Kralicek và cộng sự (2017) và hệ số chuyển đổi mặc định $R$ của IPCC (2006) ($BGB = R \times AGB$, với $R \approx 0,20 - 0,24$).
  • Không gian thu thập mẫu: Dữ liệu tập trung tại hai vùng sinh thái trọng điểm là Tây Nguyên và Đông Nam Bộ, chưa mở rộng ra các dải rừng khộp rải rác tại vùng duyên hải Nam Trung Bộ.
  • Biến động thời gian của vật liệu rơi rụng: Do đặc thù cháy rừng hàng năm thiêu rụi tầng thảm mục, nghiên cứu chưa mô hình hóa chu kỳ phát thải carbon từ vật liệu cháy theo thời gian thực.

Chương trình nghiên cứu tương lai (Future research agenda) định hình 4 hướng đột phá:

  1. Mở rộng nghiên cứu giải phẫu hệ rễ kết hợp kỹ thuật quét radar xuyên đất (Ground Penetrating Radar - GPR) để lập hệ thống SUR tích hợp cả $AGB$ và $BGB$.
  2. Ứng dụng công nghệ viễn thám LiDAR tầng quét dải hẹp kết hợp máy bay không người lái (UAV-LiDAR) để tích hợp mô hình allometric cấp chi vào bản đồ hóa carbon 3D độ phân giải siêu cao.
  3. Nghiên cứu ảnh hưởng của biến đổi khí hậu dài hạn và biến động chu kỳ cháy rừng đến tốc độ tích lũy sinh khối của các loài cây tái sinh chồi gốc.
  4. Mở rộng thử nghiệm mô hình cấp chi cho các quần xã rừng khộp tại Lào và Campuchia nhằm hoàn thiện biểu sinh khối chung cho toàn tiểu vùng sông Mê Kông.

Tác động và ảnh hưởng

  • Ảnh hưởng học thuật: Tạo lập chuẩn mực phương pháp luận mới trong phân tích sinh trắc rừng tại Việt Nam; cung cấp nguồn dữ liệu thực nghiệm chuẩn hóa cho các mạng lưới nghiên cứu carbon quốc tế. Ước tính công trình sẽ thu hút hàng trăm lượt trích dẫn trong các nghiên cứu lâm nghiệp nhiệt đới, sinh thái học phục hồi và chính sách REDD+.
  • Chuyển đổi ngành lâm nghiệp: Thay đổi căn bản phương pháp điều tra trữ lượng rừng truyền thống (vốn chỉ đo thể tích gỗ thương phẩm $V$) sang điều tra đa mục tiêu tích hợp sinh khối, cấu trúc tán và trữ lượng carbon tương đương $CO_2\text{e}$.
  • Tác động chính sách và kinh tế - xã hội: Trực tiếp đóng góp vào việc thực thi Đóng góp do Quốc gia Tự quyết định (NDC) của Việt Nam cam kết tại COP26. Cung cấp công cụ tính toán minh bạch để chi trả dịch vụ môi trường rừng (PES) từ nguồn thu tín chỉ carbon cho các cộng đồng đồng bào dân tộc thiểu số tại Tây Nguyên tham gia quản lý bảo vệ rừng khộp.

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

  • Nghiên cứu sinh và Giới học thuật: Tiếp cận một khung phương pháp luận hoàn chỉnh về hồi quy phi tuyến SUR, kỹ thuật xử lý heteroscedasticity và thuật toán thẩm định chéo K-Fold/Monte Carlo trên phần mềm R và SAS.
  • Chuyên gia phát triển dự án Carbon rừng (Carbon Project Developers): Sở hữu trọn bộ phương trình sinh khối Tier 3 được thẩm định nghiêm ngặt, đáp ứng toàn diện các tiêu chí khắt khe của UNFCCC, Gold Standard và Verra/VCS.
  • Cơ quan Quản lý Nhà nước (Tổng cục Lâm nghiệp, Sở NN&PTNT các tỉnh Tây Nguyên và Đông Nam Bộ): Có được bộ công cụ chuẩn xác để giám sát diễn biến tài nguyên rừng, định giá dịch vụ hấp thụ carbon của hệ sinh thái rừng khộp.
  • Cộng đồng địa phương và Chủ rừng: Hưởng lợi kinh tế trực tiếp từ cơ chế chia sẻ lợi ích carbon rừng thông qua quy trình đo đếm sinh khối đơn giản hóa, khả thi hóa việc tham gia giám sát carbon có sự tham gia của người dân (Participatory Carbon Monitoring - PCM).

Câu hỏi chuyên sâu

1. Đóng góp lý thuyết độc đáo nhất của luận án là gì và đã mở rộng lý thuyết nào?
Trả lời: Đóng góp độc đáo nhất là việc mở rộng Lý thuyết Hệ thống Sinh khối Cộng gộp (Parresol, 2001) và Thuyết Tương quan Sinh thái cảnh quan thông qua việc thiết lập đồng thời mô hình SUR phi tuyến có trọng số kết hợp hàm điều chỉnh sinh thái môi trường $AGB = \text{AVERAGE} \times \text{MODIFIER}$. Mô hình này giải quyết triệt để tính bất đẳng thức sinh học giữa các bộ phận thân-cành-lá-vỏ và nội suy hóa tác động của các gradient khí hậu/lập địa vào hệ số phương trình.

2. Đột phá phương pháp luận của luận án khi so sánh với các công bố quốc tế trước đây là gì?
Trả lời: So với mô hình đa loài toàn cầu của Chave và cộng sự (2005, 2014) vốn mắc sai số hệ thống lên tới 400% tại rừng khô rụng lá, và mô hình độc lập của Basuki và cộng sự (2009), luận án đã kết hợp đột phá giữa kỹ thuật hồi quy phi tuyến Maximum Likelihood trực tiếp (tránh sai số chuyển đổi ngược logarit) với quy trình kiểm định chéo khách quan K-Fold và Monte Carlo ($R = 500$). Điều này cho phép tối đa hóa việc sử dụng 222 cây mẫu chặt hạ mà vẫn cung cấp sai số độc lập không chệch.

3. Phát hiện bất ngờ và phản trực giác nhất (Counter-intuitive finding) từ dữ liệu thực nghiệm là gì?
Trả lời: Phát hiện phản trực giác nhất là: Mô hình cấp chi (Dipterocarpus, Shorea) chỉ sử dụng duy nhất một biến số đầu vào là đường kính ngang ngực ($D$) lại đạt độ chính xác tương đương hoặc vượt trội (về RMSE% và MAPE%) so với mô hình chung đa loài tích hợp cả ba biến số ($D, H, WD$). Về mặt sinh thái, sự đồng nhất về cấu trúc giải phẫu và hình thái tán ở cấp chi trong rừng khộp đã triệt tiêu sự cần thiết của biến chiều cao ($H$) và khối lượng thể tích gỗ ($WD$) – hai biến số vốn chứa đựng nhiều sai số đo đạc tại hiện trường.

4. Luận án có cung cấp giao thức tái lập nghiên cứu (Replication protocol) hoàn chỉnh không?
Trả lời: Có. Luận án chi tiết hóa toàn bộ giao thức thực địa từ quy cách ô tiêu chuẩn (0,25 ha và 1,0 ha), phương pháp giải phẫu 4 bộ phận cây mẫu, quy chuẩn sấy mẫu ở $80^\circ\text{C}-105^\circ\text{C}$ trong 48 giờ, hệ số chuyển đổi chuẩn Carbon = $0,47 \times AGB$, $CO_2\text{e} = 3,67 \times \text{Carbon}$, đến mã lệnh phân tích hồi quy SUR trên SAS (PROC MODEL) và ước lượng Maximum Likelihood trên R.

5. Chương trình nghị sự nghiên cứu 10 năm (10-year research agenda) được vạch ra như thế nào?
Trả lời: Định hướng 10 năm tới tập trung vào: (1) Tích hợp hoàn chỉnh bể carbon rễ ($BGB$) và carbon hữu cơ trong đất ($SOC$) vào hệ thống SUR toàn diện 5 bể chứa; (2) Tự động hóa kiểm kê carbon quy mô cảnh quan bằng cách liên kết hệ phương trình allometric cấp chi với dữ liệu không gian từ cảm biến LiDAR viễn thám vệ tinh/UAV; (3) Thiết lập mạng lưới ô định vị sinh thái lâu dài để lượng hóa động thái tích lũy carbon dưới tác động của biến đổi khí hậu và chế độ lửa rừng.

Kết luận

  1. Luận án đã giải quyết thành công bài toán khoa học cấp bách về định lượng sinh khối và carbon cho 650.000 ha rừng khộp tại Việt Nam, đáp ứng chuẩn mực quốc tế cao nhất (Tier 3 của IPCC).
  2. Khẳng định tính ưu việt tuyệt đối của phương pháp hồi quy phi tuyến có trọng số theo Hợp lý cực đại (Weighted Nonlinear Maximum Likelihood) so với mô hình tuyến tính hóa logarit truyền thống thông qua chuẩn định lượng Chỉ số Furnival (FI).
  3. Thiết lập thành công hệ thống phương trình sinh khối đồng thời SUR cho 4 bộ phận ($B_{st}, B_{br}, B_{le}, B_{ba}$) và toàn cây ($AGB$), đảm bảo tính cộng gộp sinh học tuyệt đối và triệt tiêu sai số tương quan chéo giữa các cơ quan cây rừng.
  4. Chứng minh đột phá thực tiễn: mô hình sinh khối cấp chi (Dipterocarpus và Shorea) đơn biến đường kính ($D$) mang lại độ tin cậy vượt trội, loại bỏ sự phụ thuộc phức tạp vào các biến số chiều cao ($H$) và khối lượng thể tích gỗ ($WD$), tạo điều kiện tối đa cho việc ứng dụng tại hiện trường điều tra lâm nghiệp.
  5. Mô hình hóa thành công hàm điều chỉnh sinh thái MODIFIER ($Altitude, P, BA$), mở ra kỷ nguyên mới trong việc tích hợp các yếu tố sinh thái cảnh quan vào phương trình sinh thái sinh trắc rừng.
  6. Cung cấp bộ công cụ đo tính carbon và $CO_2$ tương đương chuẩn xác, phục vụ trực tiếp cho Chiến lược Quốc gia về Biến đổi Khí hậu, Chương trình UN-REDD+ và thị trường tín chỉ carbon rừng bền vững tại Việt Nam và khu vực Đông Nam Á.