Tổng quan về luận án

Nghiên cứu của nghiên cứu sinh Nguyễn Thị Tình, dưới sự hướng dẫn khoa học của GS. Bảo Huy tại Trường Đại học Tây Nguyên (2021) 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), đại diện cho một công trình tiên phong trong lĩnh vực sinh trắc học rừng và lượng hóa carbon lâm nghiệp. Rừng khộp (Tropical Deciduous Dry Dipterocarp Forest - DF) là hệ sinh thái rừng khô rụng lá nhiệt đới đặc thù chỉ có ở châu Á, phân bố trải dài từ Ấn Độ, Myanmar, Thái Lan, Lào, Campuchia đến Việt Nam (Maury-Lechon & Curtet, 1998; Rundel et al., 2017). Tại Việt Nam, hệ sinh thái này chiếm diện tích khoảng 650.000 ha, tập trung chủ yếu tại vùng Tây Nguyên và Đông Nam Bộ (Huy et al., 2018). Mặc dù đóng vai trò then chốt trong dịch vụ hấp thụ $CO_2$ nhằm giảm nhẹ biến đổi khí hậu toàn cầu theo khung hành động UN-REDD+ và Công ước Khung của Liên Hợp Quốc về Biến đổi Khí hậu (UNFCCC), rừng khộp Việt Nam trong nhiều thập kỷ qua vẫn thiếu vắng một hệ thống mô hình sinh trắc chuẩn tắc, đa cấp độ phân loại và được thẩm định sai số nghiêm ngặt.

Khoảng trống nghiên cứu (research gap) cốt lõi xuất phát từ thực tế: các mô hình sinh khối toàn cầu dạng tổng quát (pantropical allometric equations) như của Brown (1997) hay Chave et al. (2005, 2014) thường bộc lộ độ lệch hệ thống rất lớn khi áp dụng vào các hệ sinh thái rừng khô có tính chất lập địa khắc nghiệt; nghiên cứu của Jara et al. (2015) và Chave et al. (2014) đã chứng minh rằng các phương trình toàn cầu có thể gây sai số cục bộ lên đến 400%. Tại khu vực Đông Nam Á, các công bố hiếm hoi của Basuki et al. (2009) tại Indonesia hay Hanpattanakit et al. (2016) tại Thái Lan chỉ mới dừng lại ở quy mô mẫu hạn chế, trong khi công trình sơ khởi của Huy et al. (2016c) và Kralicek et al. (2017) tại Việt Nam với 222 cây mẫu chặt hạ vẫn bộc lộ hạn chế do cỡ mẫu ở các cấp kính lớn còn thiếu và chưa tích hợp đồng thời các yếu tố sinh thái lập địa cũng như tính tương quan đa thành phần của cấu trúc cây gỗ.

Nghiên cứu tập trung giải quyết 3 câu hỏi nghiên cứu (Research Questions - RQ) và 3 giả thuyết tương ứng (Hypotheses - H):

  • RQ1: Phương pháp toán học nào (hóa tuyến tính logarit qua OLS hay phi tuyến có trọng số qua Maximum Likelihood) mô phỏng chính xác nhất quy luật tương quan sinh khối cây rừng khộp mà không làm lệch phương sai?
    Giả thuyết H1: Mô hình phi tuyến có trọng số ước lượng bằng Maximum Likelihood kiểm soát hiệu quả hiện tượng phân hóa phương sai (heteroscedasticity) vượt trội hơn phương pháp log-log OLS truyền thống.
  • RQ2: Việc thiết lập đồng thời hệ thống phương trình các bộ phận (thân, cành, lá, vỏ) và tổng sinh khối trên mặt đất (AGB) theo phương pháp Seemingly Unrelated Regression (SUR) có nâng cao độ tin cậy và đảm bảo tính cộng gộp sinh học hơn các mô hình độc lập?
    Giả thuyết H2: Phương pháp SUR khắc phục triệt để sai số tương quan chéo giữa các cơ quan thực vật, đảm bảo tổng sinh khối thành phần bằng chính xác $AGB$.
  • RQ3: Các nhân tố sinh thái môi trường (lượng mưa $P$, độ cao, độ dốc) và cấu trúc lâm phần (tiết diện ngang $BA$, mật độ $N$) đóng vai trò như thế nào trong việc hiệu chỉnh độ chính xác của hàm sinh khối?
    Giả thuyết H3: Mô hình tích hợp cấu trúc hàm điều chỉnh sinh thái ($AGB = \text{AVERAGE} \times \text{MODIFIER}$) làm giảm đáng kể sai số ước tính sinh khối ở cấp độ lâm phần.

Khung lý thuyết của luận án được xây dựng dựa trên Lý thuyết Tỷ lệ Sinh trắc học (Allometric Scaling Theory - West, Brown & Enquist, 1997, 1999), Lý thuyết Ước lượng Hệ phương trình Tương quan Dường như Không liên quan (SUR - Zellner, 1962; Parresol, 2001) và Nguyên lý Điều chỉnh Sinh thái học Quần xã (Community Functional Ecology). Về quy mô và phạm vi, nghiên cứu thu thập dữ liệu trên 28 ô tiêu chuẩn 0,25 ha ($50 \times 50\text{ m}$) tại Tây Nguyên (CH) và 1 ô định vị 1,0 ha ($100 \times 100\text{ m}$) tại Đông Nam Bộ (SE), khảo sát trên các quần hợp ưu thế gồm 4 loài chủ chốt: 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.), kết hợp tập dữ liệu giải phẫu chặt hạ 222 cây mẫu kế thừa từ chương trình UN-REDD cùng các cây bổ sung ở các cấp kính lớn ($D \ge 5\text{ cm}$ đến $D_{\max} \approx 50\text{ cm}$).

Literature Review và Positioning

Nghiên cứu sinh khối và định lượng carbon rừng nhiệt đới trong 3 thập kỷ qua được phân hóa thành ba dòng học thuật chính:

  1. Dòng mô hình sinh khối toàn cầu (Pantropical Global Models): Khởi xướng bởi Brown et al. (1989, 1997), Ketterings et al. (2001), Chave et al. (2005, 2014). Trường phái này xây dựng các phương trình đa loài dựa trên tập hợp dữ liệu lớn từ nhiều châu lục với các biến số độc lập tiêu chuẩn gồm đườ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$). Dù có tính tổng quát cao, dòng nghiên cứu này vấp phải sự chỉ trích gay gắt từ Picard et al. (2015) và Jara et al. (2015) do tính chất che lấp các đặc trưng hình thái học đặc thù của từng hệ sinh thái cục bộ.
  2. Dòng mô hình sinh thái chuyên biệt theo vùng (Regional & Forest-Type Specific Models): Được thúc đẩy bởi Basuki et al. (2009) tại Indonesia, Hanpattanakit et al. (2016) tại Thái Lan, Cairns et al. (2003) tại Mexico, và Bảo Huy et al. (2012, 2016a,b,c) tại Việt Nam. Các tác giả khẳng định hệ sinh thái rừng khô rụng lá họ Dầu có tỷ lệ phân chia sinh khối cành - thân - vỏ hoàn toàn dị biệt so với rừng mưa thường xanh nhiệt đới ẩm, do cơ chế thích nghi giải phẫu với chu kỳ cháy rừng và khô hạn 3-4 tháng kéo dài hàng năm (Erskine & Huy, 2003).
  3. Dòng mô hình nâng cao tích hợp phụ thuộc tương quan và biến số lâm phần: Được phát triển bởi Parresol (2001), Temesgen et al. (2007, 2015), Poudel & Temesgen (2016), và Kralicek et al. (2017). Hướng tiếp cận này xem xét cơ thể cây gỗ như một hệ thống cấu trúc hữu cơ có tương quan nội tại giữa các bộ phận, yêu cầu mô hình hóa đồng thời và lồng ghép các yếu tố hoàn cảnh rừng.

Trong y văn tồn tại hai cuộc tranh luận học thuật lớn. Thứ nhất là tranh luận về phương pháp tuyến tính hóa logarit ($ \ln(AGB) = \ln(a) + b\ln(D) $) sử dụng bình phương tối thiểu (OLS) có nhân tử hiệu chỉnh Baskerville (CF) so với phương pháp hồi quy phi tuyến trực tiếp (Weighted Non-Linear Regression). Các tác giả truyền thống (Brown, 1997; Chave et al., 2005) ủng hộ biến đổi logarit vì tính đơn giản, trong khi Picard et al. (2012) và Bảo Huy (2017a,b) chỉ ra rằng phép biến đổi logarit làm biến dạng cấu trúc sai số thực tế và đánh giá sai phương sai ở các cây có kích thước lớn. Thứ hai là tranh luận về tính độc lập hay phụ thuộc của các phương trình thành phần sinh khối: việc lập mô hình riêng lẻ từng bộ phận thân ($Bst$), cành ($Bbr$), lá ($Ble$), vỏ ($Bba$) dẫn đến hiện tượng tổng sinh khối các bộ phận không bằng sinh khối toàn cây ($AGB \ne Bst + Bbr + Ble + Bba$), tạo ra độ chệch nghiêm trọng trong các hệ thống Đo đạc, Báo cáo và Thẩm định (MRV) của REDD+.

So sánh với hai nghiên cứu quốc tế điển hình:

  • Basuki et al. (2009) tại Indonesia: Phát triển mô hình $AGB$ cho các loài thuộc họ Dầu ở vùng thấp nhưng chỉ áp dụng hồi quy tuyến tính hóa logarit OLS và dừng lại ở cấp độ biến $D$, chưa giải quyết được tính cộng gộp thành phần và chưa kiểm định chéo K-fold.
  • Hanpattanakit et al. (2016) tại Thái Lan: Xây dựng phương trình cho 5 loài rừng khộp nhưng sử dụng cỡ mẫu nhỏ, mô hình thành phần tách rời và không kiểm soát sự tương quan môi trường lâm phần.

Luận án này định vị chính xác vào giao điểm của các khoảng trống học thuật trên: lần đầu tiên thiết lập hệ thống phương trình sinh khối đồng thời (SUR) phi tuyến có trọng số, trải rộng từ cấp loài, chi (Dipterocarpus, Shorea), họ (Dipterocarpaceae) đến cấp hỗn loài rừng khộp, kết hợp cấu trúc hàm điều chỉnh sinh thái môi trường ($P$, $BA$) và thẩm định chéo sai số toàn diện.

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

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

Luận án đã mở rộng và làm sâu sắc thêm Lý thuyết Phân bổ Sinh trắc học Thực vật (Plant Allometric Partitioning Theory) của Enquist & Niklas (2001) trong điều kiện lập địa khắc nghiệt của rừng khô rụng lá nhiệt đới:

  1. Khẳng định tính ưu việt của mô hình phân cấp phân loại học (Taxonomic Hierarchy Specificity): Nghiên cứu chứng minh rằng ở cấp độ chi thực vật (Dipterocarpus và Shorea), mối quan hệ sinh trắc giữa đường kính ngang ngực ($D$) và tổng sinh khối ($AGB$) đạt mức độ ổn định hình thái rất cao. Việc chỉ sử dụng một biến số duy nhất là $D$ ở cấp chi đem lại độ tin cậy và chỉ số thông tin tương đương hoặc vượt trội so với mô hình hỗn loài chung phải tích hợp cả 3 biến số phức tạp ($D$, $H$, $WD$).
  2. Minh chứng quy luật cộng gộp cấu trúc sinh khối cây rừng qua SUR: Đóng góp về mặt lý thuyết sinh thái - toán học khi chỉ ra rằng ma trận hiệp phương sai sai số giữa các cơ quan sinh dưỡng (thân, cành, lá, vỏ) chứa đựng các tương tác sinh lý nội tại; việc mô hình hóa đồng thời giải tỏa hoàn toàn độ chệch ước tính, xác lập một bước chuyển biến mô hình (paradigm shift) từ phương pháp hồi quy đơn biến rời rạc sang tiếp cận hệ thống động lực tích hợp.
  3. Mở rộng lý thuyết tương tác môi trường - sinh trắc (Ecological-Allometric Interaction Theory): Cung cấp bằng chứng thực nghiệm khẳng định sinh khối cây rừng cá thể không chỉ là hàm số của kích thước hình thái mà còn chịu sự điều tiết có ý nghĩa thống kê của lượng mưa năm ($P$) và mật độ không gian cạnh tranh qua tổng tiết diện ngang lâm phần ($BA$).

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

Khung phân tích của luận án tích hợp liên ngành giữa:

  • Lý thuyết Sinh trắc Rừng dạng Hàm Lũy thừa (Power Function Formulation): $y_i = a \cdot X_j^b$ hoặc $y_i = a \cdot (D^2 H WD)^b$.
  • Hệ thống Phương trình Hồi quy Dường như Không tương quan (SUR): Cho phép ước lượng đồng thời hệ thống gồm 5 phương trình ($Bst, Bbr, Ble, Bba, AGB$) dưới dạng phi tuyến có trọng số (Weighted Nonlinear SUR): $$\begin{cases} Bst = a_1 \cdot D^{b_1} \ Bbr = a_2 \cdot D^{b_2} \ Ble = a_3 \cdot D^{b_3} \ Bba = a_4 \cdot D^{b_4} \ AGB = Bst + Bbr + Ble + Bba \end{cases}$$
  • Mô hình Hiệu chỉnh Môi trường Lâm phần (Environmental Modifier Framework): $$AGB = \text{AVERAGE} \times \text{MODIFIER} = a \cdot D^b \cdot WD^d \times \exp\left[b_2(P - 1502) + b_3(BA - 12.5)\right]$$ Trong đó giá trị 1502 mm và 12,5 $\text{m}^2\text{/ha}$ là các giá trị trung bình chuẩn hóa của lượng mưa năm và tổng tiết diện ngang lâm phần nghiên cứu.

Điều kiện biên (boundary conditions) của khung phân tích được xác định tường minh: áp dụng cho cây gỗ tự nhiên có $D \ge 5\text{ cm}$, sinh trưởng trong các kiểu rừng nhiệt đới khô rụng lá ưu thế họ Dầu ở dải cao độ 171 - 417 m, lượng mưa 1000 - 1600 mm/năm và có mùa khô sâu sắc kéo dài từ 3 đến 4 tháng.

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

Thiết kế nghiên cứu

Nghiên cứu theo đuổi trường phái thực chứng định lượng chặt chẽ (Empirical Positivism). Thiết kế nghiên cứu đa tầng (Multi-level Hierarchical Design) kết hợp giữa điều tra lâm phần thực địa không tàn phá và phương pháp rút mẫu chặt hạ giải phẫu (Destructive Sampling) toàn diện nhằm thu thập dữ liệu chân thực tuyệt đối (ground-truth data).

Quy mô mẫu và tiêu chí lựa chọn:

  • Khảo sát lâm phần trên 29 ô tiêu chuẩn định vị (28 ô tại Tây Nguyên và 1 ô tại Đông Nam Bộ).
  • Toàn bộ cây mẫu chặt hạ ($N = 222$ cây kế thừa từ chương trình UN-REDD của GS. Bảo Huy cùng các cây bổ sung) được phân bổ theo tỷ lệ phân bố mật độ theo cấp kính ($N/D$) của lâm phần thực tế, đảm bảo tính đại diện sinh học từ cấp kính nhỏ ($D \approx 5\text{ cm}$) đến cây cổ thụ ($D \approx 50\text{ cm}$).

Quy trình nghiên cứu rigorous

Quy trình thu thập và xử lý mẫu tuân thủ nghiêm ngặt chuẩn mực quốc tế của IPCC (2003, 2006) và Silva (2011):

  1. Phân đoạn sinh khối tươi tại hiện trường: Cây mẫu sau khi chặt hạ được tách thành 4 phân đoạn cơ thể độc lập: Thân gỗ ($Bst$), cành ($Bbr$), lá ($Ble$), và vỏ cây ($Bba$). Toàn bộ sinh khối tươi của từng phân đoạn được cân trực tiếp tại hiện trường bằng cân đồng hồ cơ học và cân điện tử có độ chính xác cao.
  2. Lấy mẫu phụ và sấy khô tuyệt đối: Lấy mẫu đại diện cho từng bộ phận (đĩa gỗ thân ở các vị trí gốc, ngang ngực, giữa thân, ngọn; mẫu cành các cấp; mẫu lá hỗn hợp; mẫu vỏ), niêm phong và vận chuyển về phòng thí nghiệm. Mẫu được sấy ở nhiệt độ $80^\circ\text{C} - 105^\circ\text{C}$ liên tục trong 48 giờ cho đến khi khối lượng không đổi để xác định tỷ lệ sinh khối khô/tươi ($DW/FW$).
  3. Quy đổi trữ lượng Carbon và $CO_2$ tương đương: Áp dụng hệ số chuyển đổi quốc tế chuẩn hóa: $$\text{Carbon} = 0,47 \times \text{Sinh khối khô} \quad (\text{IPCC, 2006; Bảo Huy, 2013})$$ $$\text{CO}_2\text{ tương đương} = 3,67 \times \text{Carbon}$$
  4. Tam giác đạc và kiểm soát độ tin cậy (Triangulation & Reliability): Kiểm tra chéo giữa dữ liệu giải phẫu thực tế, cơ sở dữ liệu khối lượng thể tích gỗ toàn cầu (Global Wood Density Database - Chave et al., 2009) và hệ thống đo đạc sinh thái môi trường WorldClim (Hijmans et al., 2005; Fischer et al., 2008).

Data và phân tích

Đặc trưng thống kê lâm phần nghiên cứu thể hiện sự biến động phong phú: mật độ lâm phần ($N$) dao động từ 228 đến 1291 cây/ha (đối với cây $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}$, độ cao tuyệt đối từ 171 đến 417 m (trung bình 246 m, độ lệch chuẩn $Sd = 64,6\text{ m}$), độ dốc từ 0 đến $8^\circ$ (trung bình $2,0^\circ$, $Sd = 2,2^\circ$).

Kỹ thuật phân tích dữ liệu nâng cao:

  • Xử lý hiện tượng phân hóa phương sai (Heteroscedasticity): Sử dụng hàm trọng số $w = 1/D^k$ hoặc $w = 1/Var(e)$ tích hợp trong mô hình phi tuyến ước lượng bằng phương pháp Hợp lý cực đại (Weighted Nonlinear fit by Maximum Likelihood) thực hiện trên ngôn ngữ R.
  • So sánh mô hình khác biệt biến phụ thuộc qua Chỉ số Furnival (Furnival’s Index - FI): $$FI = \left[\exp\left(\frac{1}{n}\sum \ln(y)\right)\right] \times RMSE_{\text{log}}$$ Mô hình có chỉ số $FI$ thấp hơn phản ánh chất lượng điều chỉnh dữ liệu vượt trội mà không bị thiên kiến bởi dạng biến đổi toán học.
  • Mô hình hồi quy hệ thống SUR trên phần mềm SAS: Sử dụng thủ tục PROC MODEL để ước lượng đồng thời hệ phương trình phi tuyến tính có trọng số, ràng buộc tham số tổng sinh khối.
  • Quy trình Thẩm định chéo (Cross-Validation Pipeline): Triển khai đồng thời Thẩm định chéo K-lần (K-Fold với $K = 5 - 10$), Thẩm định loại một (Leave-One-Out - LOOCV) và mô phỏng ngẫu nhiên lặp Monte Carlo ($R$ lần lặp với tỷ lệ 70% Training / 30% Validation). Đánh giá sai số qua ba chỉ số định lượng: Sai lệch phần trăm (Bias %), Sai số trung phương tương đối (RMSE %), và Sai số tuyệt đối trung bình phần trăm (MAPE %). Đồ thị Bland-Altman được áp dụng để kiểm định giới hạn tương đồng sai số ở độ tin cậy 95%.

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

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

  1. Phép ước lượng Phi tuyến có trọng số vượt trội hoàn toàn so với Log-Linear OLS: Kết quả tính toán chỉ số Furnival (FI) chứng minh phương pháp phi tuyến ước lượng bằng Maximum Likelihood có trọng số triệt tiêu hiện tượng phương sai thay đổi (heteroscedasticity) khi kích thước cây tăng lên, đạt độ chuẩn xác cao hơn hẳn phương pháp hồi quy tuyến tính hóa logarit truyền thống ở mọi cấp phân loại.
  2. Hệ thống mô hình SUR xác lập tính cộng gộp tuyệt đối và giảm thiểu sai số dự báo: Việc áp dụng phương pháp SUR để thiết lập đồng thời các hàm sinh khối thành phần ($Bst, Bbr, Ble, Bba$) và $AGB$ đã giải quyết triệt để sự mâu thuẫn giữa tổng các phương trình thành phần độc lập và phương trình toàn cây. Kết quả thẩm định chéo K-Fold cho thấy hệ thống SUR làm giảm đáng kể sai số toàn phần (MAPE% và RMSE%), bảo đảm tính vững chắc sinh học trong cấu trúc phân bổ vật chất.
  3. Mô hình cấp Chi thực vật ưu thế với một biến số ($D$) có độ chính xác vượt trội mô hình hỗn loài đa biến: Đối với hai chi ưu thế Dipterocarpus (Dầu) và Shorea (Cà chít), phương trình dạng $AGB = a \cdot D^b$ ước lượng theo SUR cho độ tin cậy cao hơn, sai số dự báo thấp hơn so với mô hình chung hỗn loài sử dụng 3 biến ($D, H, WD$). Điều này chứng minh tính đồng nhất cao về mặt hình thái học giải phẫu của các loài trong cùng một chi thuộc hệ sinh thái rừng khộp.
  4. Xác lập mô hình điều chỉnh tích hợp nhân tố lâm phần và lượng mưa: Phương trình tổng quát dạng: $$AGB = a \cdot D^b \cdot WD^d \times \exp\left[b_2(P - 1502) + b_3(BA - 12.5)\right]$$ được ước lượng thành công bằng phương pháp phi tuyến cố định có trọng số. Sự tham gia của lượng mưa ($P$) và tổng tiết diện ngang ($BA$) đóng vai trò biến số điều tiết (moderators), phản ánh chính xác sự biến thiên sinh thái giữa vùng khô hạn Đông Nam Bộ ($P \approx 1003\text{ mm}$) và vùng Tây Nguyên ($P \approx 1600\text{ mm}$).
  5. Bác bỏ tính ứng dụng trực tiếp của các mô hình Pantropical toàn cầu cho rừng khộp Việt Nam: Khi kiểm định chéo các mô hình toàn cầu (như Chave et al., 2005, 2014) trên tập dữ liệu cây mẫu rừng khộp Việt Nam, sai số hệ thống xuất hiện ở mức rất cao do cấu trúc phân cành thấp, thân cong queo và vỏ dày đặc trưng chịu lửa của cây họ Dầu khô rụng lá; trong khi đó, mô hình cấp chi của luận án có sự tương thích cao với các nghiên cứu vùng tại Indonesia (Basuki et al., 2009).

Implications đa chiều

  • Về mặt lý thuyết: Đặt nền móng khoa học cho việc phân loại sinh thái - sinh trắc học cây rừng nhiệt đới, khẳng định vai trò của cấp phân loại học chi (genus-level allometry) trong việc tối ưu hóa quy luật tương quan sinh trưởng.
  • Về phương pháp luận: Thiết lập quy chuẩn xử lý dữ liệu sinh trắc học hiện đại tại Việt Nam: kết hợp ngôn ngữ R và SAS để triển khai hệ thống hồi quy phi tuyến SUR có trọng số và quy trình thẩm định chéo K-Fold/Monte Carlo, thay thế hoàn toàn phương pháp chia tách dữ liệu kiểm chứng độc lập đơn lẻ truyền thống vốn gây lãng phí dữ liệu giải phẫu quý giá.
  • Về ứng dụng thực tiễn: Cung cấp trọn bộ công cụ bảng tra sinh khối và phần mềm tính toán trữ lượng carbon chính xác cho 650.000 ha rừng khộp, phục vụ trực tiếp cho các công ty lâm nghiệp, ban quản lý rừng đặc dụng/phòng hộ tại Đắk Lắk, Gia Lai, Đắk Nông, Bình Phước, Tây Ninh.
  • Về chính sách và kinh tế carbon: Cung cấp hệ số phát thải và hấp thụ carbon cấp quốc gia đạt chuẩn Tier 3 của IPCC, tạo cơ sở kỹ thuật vững chắc để Việt Nam tham gia thị trường tín chỉ carbon quốc tế, đàm phán chi trả tài chính dịch vụ môi trường rừng trong các thỏa thuận REDD+ và Thỏa thuận Paris (Điều 6).

Limitations và Future Research

Luận án thẳng thắn thừa nhận các giới hạn nghiên cứu sau:

  1. Giới hạn về dữ liệu cây mẫu kích thước cực đại ($D > 50\text{ cm}$): Mặc dù đã bổ sung mẫu, nhưng do phân bố mật độ tự nhiên theo cấp kính của rừng khộp có dạng giảm nhanh, số lượng cây mẫu giải phẫu ở các cấp kính lớn ($D > 50\text{ cm}$) vẫn chiếm tỷ lệ khiêm tốn, có thể tạo ra độ bất định nhất định khi ngoại suy cho các lâm phần nguyên sinh cổ thụ.
  2. Hạn chế về dữ liệu sinh khối dưới mặt đất (BGB - Hệ rễ): Nghiên cứu tập trung chủ yếu vào sinh khối trên mặt đất ($AGB$). Việc khai quật hệ rễ cây họ Dầu trên nền đất xám sỏi sạn xương xẩu cực kỳ tốn kém và phá hủy môi trường, do đó việc ước tính $BGB$ trong luận án vẫn phải kế thừa từ các công bố chuyên sâu của Kralicek et al. (2017) hoặc sử dụng hệ số chuyển đổi mặc định $R = 0,20 - 0,24$ (IPCC, 2006).
  3. Tính động thái theo mùa của sinh khối lá ($Ble$): Rừng khộp có đặc tính rụng lá hoàn toàn vào mùa khô; mô hình sinh khối lá của luận án phản ánh trạng thái tán lá phát triển tối đa trong mùa sinh trưởng (mùa mưa), chưa mô phỏng liên tục biến thiên rụng lá theo từng tháng trong năm.

Chương trình nghiên cứu tương lai (10-Year Research Agenda):

  • Tích hợp công nghệ viễn thám LiDAR và ảnh siêu phổ (Hyperspectral): Kết hợp mô hình SUR mặt đất với dữ liệu quét LiDAR không gian để tự động hóa kiểm kê carbon trên toàn bộ 650.000 ha rừng khộp.
  • Nghiên cứu giải phẫu rễ sâu bằng công nghệ khí nén Supersonic Air-Knife: Khai quật không tàn phá hệ rễ để hoàn thiện hệ phương trình SUR tích hợp cả 6 bộ phận ($Bst, Bbr, Ble, Bba, Bstump, Broot$).
  • Mạng lưới giám sát biến đổi khí hậu dài hạn (Long-term Ecological Research Plots - LTER): Theo dõi tác động của hiện tượng El Niño/Hạn hán cực đoan đến tốc độ tích lũy sinh khối và tỷ lệ tử vong của các loài cây họ Dầu ưu thế.

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

  • Tác động học thuật: Định hình tiêu chuẩn nghiên cứu sinh trắc học rừng tại Việt Nam và Đông Nam Á; ước tính tiềm năng tạo ra các trích dẫn quốc tế quan trọng trên các tạp chí chuyên ngành lâm nghiệp hàng đầu (Forest Ecology and Management, Annals of Forest Science).
  • Chuyển đổi ngành lâm nghiệp: Chuyển đổi phương thức kiểm kê rừng từ ước tính trữ lượng gỗ ($m^3$) đơn thuần sang định lượng chính xác sinh khối đa thành phần và trữ lượng carbon ($tC/ha$), thúc đẩy chuyển dịch cơ cấu kinh tế lâm nghiệp sang kinh tế dịch vụ môi trường.
  • Tác động chính sách quốc gia: Trực tiếp hỗ trợ Tổng cục Lâm nghiệp (Bộ Nông nghiệp & PTNT) hoàn thiện Báo cáo Định kỳ Quốc gia (National Communications) gửi UNFCCC, củng cố tính minh bạch của Hệ thống Đo đạc, Báo cáo và Thẩm định (MRV) quốc gia.
  • Lợi ích xã hội và cộng đồng: Tạo cơ sở pháp lý và định lượng khoa học để thực hiện chi trả tiền dịch vụ môi trường rừng từ thị trường carbon đến các cộng đồng nghèo, đồng bào dân tộc thiểu số tại Tây Nguyên tham gia nhận khoán 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 Lâm nghiệp: Tiếp cận phương pháp luận chuẩn tắc về mô hình hóa phi tuyến có trọng số, kỹ thuật SUR trên SAS và quy trình thẩm định chéo K-Fold/Monte Carlo trên R để áp dụng cho các hệ sinh thái rừng khác.
  • Các Ban Quản lý Rừng và Công ty Lâm nghiệp: Ứng dụng trực tiếp hệ thống phương trình theo loài và chi ưu thế để tính toán chính xác biến động carbon của từng lâm phần quản lý mà không cần chặt hạ cây.
  • Các Nhà phát triển Dự án Carbon Rừng (Carbon Project Developers): Rút ngắn thời gian, giảm thiểu chi phí đo đếm hiện trường và vượt qua các kỳ thẩm định khắt khe của các tổ chức chứng nhận tín chỉ carbon quốc tế (Verra VCS, Plan Vivo, Gold Standard).
  • Cơ quan Quản lý Nhà nước và Hoạch định Chính sách: Nắm giữ bộ công cụ định lượng khoa học để quy hoạch vùng bảo tồn rừng khộp, xây dựng chính sách hạn chế chuyển đổi rừng khộp nghèo kiệt sang trồng cây nông nghiệp ngắn ngày.

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à đã mở rộng Lý thuyết Sinh thái Sinh trắc học Phân cấp (Hierarchical Allometric Scaling Theory) bằng việc chứng minh: Trong hệ sinh thái rừng khô rụng lá họ Dầu, tính đồng nhất hình thái học ở cấp chi (Dipterocarpus và Shorea) mạnh đến mức mô hình sinh khối đơn biến ($D$) cấp chi đạt độ chính xác tương đương hoặc vượt trội mô hình hỗn loài 3 biến ($D, H, WD$). Đồng thời, hoàn thiện lý thuyết tương quan hệ thống thực vật thông qua việc ứng dụng SUR để giải quyết triệt để tính cộng gộp sinh thái của các cơ quan cây gỗ.

2. Đột phá về phương pháp luận của luận án khi so sánh với ít nhất 2 nghiên cứu quốc tế trước đó?
Trả lời: So với Basuki et al. (2009) (Indonesia) sử dụng hồi quy tuyến tính hóa logarit OLS đơn lẻ từng biến và Chave et al. (2014) (Toàn cầu) sử dụng phương trình tổng quát không xét tương quan bộ phận, luận án đã tạo bước đột phá phương pháp luận khi: (1) Ứng dụng kỹ thuật hồi quy phi tuyến có trọng số qua Maximum Likelihood trên R triệt tiêu hiện tượng heteroscedasticity; (2) Thiết lập đồng thời hệ phương trình thành phần qua SUR trên SAS đảm bảo $AGB = Bst + Bbr + Ble + Bba$; và (3) Chuẩn hóa quy trình thẩm định chéo K-Fold và Monte Carlo giúp tận dụng 100% dữ liệu giải phẫu mà vẫn cung cấp sai số khách quan toàn cầu.

3. Phát hiện nào trong nghiên cứu mang tính bất ngờ nhất và được chứng minh bằng dữ liệu ra sao?
Trả lời: Phát hiện bất ngờ nhất là việc đưa thêm biến chiều cao cây ($H$) và khối lượng thể tích gỗ ($WD$) vào mô hình sinh khối của chi Dipterocarpus không làm tăng đáng kể độ chính xác của mô hình so với việc chỉ sử dụng duy nhất biến đường kính ($D$). Dữ liệu thẩm định chéo K-Fold cho thấy chỉ số sai số MAPE% và RMSE% của mô hình một biến $AGB = a \cdot D^b$ cho chi Dipterocarpus sai khác không có ý nghĩa thống kê so với mô hình đa biến phức tạp, do đặc tính thân cây rừng khộp có quan hệ hình học $H/D$ cực kỳ ổn định dưới áp lực chọn lọc của gió lốc và lửa rừng hàng năm.

4. Nghiên cứu có cung cấp giao thức tái lập (Replication Protocol) rõ ràng không?
Trả lời: Luận án cung cấp giao thức tái lập chi tiết ở mức độ tuyệt đối: từ quy chuẩn lập 29 ô tiêu chuẩn định vị ($50 \times 50\text{ m}$ và $100 \times 100\text{ m}$), phương pháp rút mẫu chặt hạ theo phân bố cấp kính $N/D$, quy trình sấy mẫu chuẩn $80 - 105^\circ\text{C}$ trong 48 giờ, đến toàn bộ cấu trúc hàm toán học, cú pháp mã lệnh ước lượng SUR trên SAS (PROC MODEL), mã phân tích thẩm định chéo K-Fold trên R và các hệ số tham số $a, b, c, d$ cùng giá trị chuẩn hóa môi trường ($P = 1502\text{ mm}, BA = 12,5\text{ m}^2\text{/ha}$).

5. Lộ trình nghiên cứu 10 năm tới được phác thảo như thế nào?
Trả lời: Lộ trình 10 năm tập trung vào 3 hướng mũi nhọn: (1) Mở rộng hệ thống mô hình SUR cho sinh khối dưới mặt đất ($BGB$) và carbon hữu cơ trong đất ($SOC$) bằng phương pháp đào rễ cơ giới hóa; (2) Tích hợp phương trình sinh khối mặt đất với dữ liệu quét viễn thám LiDAR đa thời gian và máy bay không người lái (UAV) để lập bản đồ biến động carbon thời gian thực; và (3) Xây dựng mạng lưới ô mẫu định vị quan trắc động thái chu chuyển carbon rừng khộp dài hạn dưới tác động của biến đổi khí hậu và tần suất cháy rừng.

Kết luận

Luận án tiến sĩ của Nguyễn Thị Tình đã giải quyết trọn vẹn và xuất sắc mục tiêu khoa học đề ra thông qua 6 đóng góp cốt lõi:

  1. Xác lập hệ thống mô hình sinh khối chuẩn tắc: Lần đầu tiên xây dựng hoàn chỉnh hệ thống mô hình sinh khối trên mặt đất ($AGB$) và 4 bộ phận ($Bst, Bbr, Ble, Bba$) cho rừng khộp Việt Nam từ cấp loài, chi, họ đến hỗn loài.
  2. Khẳng định tính ưu việt của phương pháp SUR phi tuyến có trọng số: Giải quyết triệt để tính cộng gộp sinh học và kiểm soát hiện tượng phân hóa phương sai, khắc phục hoàn toàn các sai số thiên kiến của phương pháp hồi quy logarit truyền thống.
  3. Chứng minh hiệu quả vượt trội của mô hình cấp Chi thực vật: Khẳng định mô hình một biến $D$ cho hai chi ưu thế Dipterocarpus và Shorea có độ tin cậy và tính ứng dụng thực tiễn cao nhất.
  4. Tích hợp thành công hàm điều chỉnh sinh thái lâm phần: Định lượng hóa rõ ràng ảnh hưởng của lượng mưa năm ($P$) và tiết diện ngang ($BA$) lên biến thiên sinh khối cây cá thể ($AGB = \text{AVERAGE} \times \text{MODIFIER}$).
  5. Thiết lập quy chuẩn thẩm định chéo nghiêm ngặt: Cung cấp bộ chỉ số sai số khách quan qua K-Fold, LOOCV và Monte Carlo, bác bỏ tính thích dụng trực tiếp của các mô hình Pantropical toàn cầu cho rừng khộp Việt Nam.
  6. Cung cấp cơ sở định lượng đạt chuẩn IPCC Tier 3: Mở ra kỷ nguyên chuyển đổi số trong quản lý tài nguyên rừng và tạo lập đòn bẩy kỹ thuật vững chắc để Việt Nam thương mại hóa tín chỉ carbon rừng khộp trên thị trường quốc tế.