Tổng quan về luận án

Hệ sinh thái rừng khộp (Tropical Deciduous Dry Dipterocarp Forest - DF) giữ vai trò then chốt trong mạng lưới sinh thái nhiệt đới châu Á nhờ khả năng thích nghi đặc thù với chế độ khô hạn, cháy rừng chu kỳ và tiềm năng tích lũy carbon sinh khối to lớn phục vụ giảm nhẹ biến đổi khí hậu toàn cầu. Tại Việt Nam, "ở Việt Nam rừng khộp có khoảng 650.000 ha phân bố chủ yếu ở Tây Nguyên và Đông Nam bộ" (Huy và ctv, 2018). Trong khuôn khổ thực thi Hiệp định Khung của Liên Hợ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+), yêu cầu cấp bách đặt ra là phải xây dựng hệ thống Đo lường, Báo cáo và Thẩm định (MRV) đạt chuẩn Cấp độ 3 (Tier 3) của Ủy ban Liên chính phủ về Biến đổi Khí hậu (IPCC, 2006).

Khoảng trống nghiên cứu (Research Gap) cốt lõi xuất phát từ thực trạng: các phương trình tương quan sinh khối chung cho toàn vùng nhiệt đới (Pan-tropical allometric models) như của Chave và ctv (2005, 2014) khi áp dụng vào các hệ sinh thái đặc thù cục bộ có thể dẫn đến "sai số hệ thống lên đến 400%" (Jara và ctv, 2015). Đồng thời, các nghiên cứu sinh khối trước đây tại Việt Nam và Đông Nam Á (Basuki và ctv, 2009; Huy và ctv, 2016c; Kralicek và ctv, 2017) chủ yếu lập phương trình độc lập cho từng bộ phận (thân, cành, lá, vỏ), dẫn đến sự triệt tiêu tính cộng gộp sinh học (Additivity property) và bỏ qua tương quan chéo giữa các phương trình thành phần. Ngoài ra, việc thiếu hụt dữ liệu cây mẫu chặt hạ kích thước lớn ($D \ge 40\text{ cm}$) và chưa tích hợp ảnh hưởng tổng hợp của các nhân tố sinh thái môi trường (lượng mưa, độ cao, mật độ, tổng tiết diện ngang) vào mô hình sinh trắc đã hạn chế độ tin cậy khi mở rộng quy mô dự báo cấp cảnh quan.

Luận án tiến sĩ lâm nghiệp của NCS. Nguyễn Thị Tình (Trường Đại học Tây Nguyên, 2021) dưới sự hướng dẫn khoa học của GS. Bảo Huy đã giải quyết triệt để các khoảng trống trên thông qua 4 câu hỏi nghiên cứu và giả thuyết tương ứng:

  • Câu hỏi nghiên cứu 1 (RQ1): Dạng hàm phi tuyến cực đại hóa hợp lý có trọng số (Weighted Non-Linear Maximum Likelihood) có ưu việt hơn mô hình dạng mũ tuyến tính hóa qua logarit trong việc khắc phục hiện tượng phương sai thay đổi (heteroscedasticity) hay không?
  • Giả thuyết 1 (H1): Mô hình phi tuyến có trọng số ước lượng trực tiếp bằng Maximum Likelihood đạt chỉ số Furnival (FI) và sai số RMSE% thấp hơn mô hình OLS log-log truyền thống.
  • Câu hỏi nghiên cứu 2 (RQ2): Phương pháp hồi quy dường như không liên quan (Seemingly Unrelated Regression - SUR) có đảm bảo tính cộng gộp sinh khối bộ phận và cải thiện độ chính xác dự đoán so với các mô hình thành phần độc lập không?
  • Giả thuyết 2 (H2): Hệ thống mô hình đồng thời SUR tối ưu hóa hiệp phương sai phần dư, triệt tiêu sai lệch cộng dồn khi ước tính sinh khối trên mặt đất (AGB).
  • Câu hỏi nghiên cứu 3 (RQ3): Mô hình sinh khối phân cấp theo phân loại thực vật (chi, họ) với biến số duy nhất là đường kính ngang ngực ($D$) có đạt độ tin cậy vượt trội hơn mô hình đa biến ($D, H, WD$) chung cho vùng nhiệt đới?
  • Giả thuyết 3 (H3): Đặc trưng giải phẫu và hình thái của chi ưu thế (Dipterocarpus, Shorea) chi phối chặt chẽ phân bổ sinh khối hơn các biến hình thái lâm phần chung.
  • Câu hỏi nghiên cứu 4 (RQ4): Hàm điều chỉnh môi trường (MODIFIER) có giải thích được biến động sinh khối giữa các vùng sinh thái Tây Nguyên và Đông Nam Bộ?
  • Giả thuyết 4 (H4): Tích hợp lượng mưa ($P$), tổng tiết diện ngang ($BA$) và mật độ ($N$) qua hàm $AGB = AVERAGE \times MODIFIER$ làm giảm sai số ngoại suy không gian.

Nghiên cứu được neo trên nền tảng Lý thuyết Tỷ lệ Sinh trắc học (Allometric Scaling Theory - West, Brown & Enquist, 1997, 1999), Lý thuyết Cấu trúc Lâm phần (Assmann, 1970) và Lý thuyết Ước lượng Hệ thống Đồng thời (Zellner, 1962; Parresol, 2001). Phạm vi nghiên cứu bao quát 29 ô tiêu chuẩn định vị tại hai vùng sinh thái Tây Nguyên (CH: 28 ô, diện tích 0,25 ha/ô) và Đông Nam Bộ (SE: 1 ô, diện tích 1,0 ha), thực hiện phân tích giải phẫu và chặt hạ trên 222 cây mẫu thực nghiệm, đại diện cho dải đường kính $D$ từ 5,0 đến 68,0 cm và chiều cao $H$ từ 4,2 đến 20,5 m.

Literature Review và Positioning

Lịch sử phát triển các mô hình sinh trắc ước tính sinh khối cây rừng nhiệt đới ghi nhận ba dòng chảy học thuật chính. Dòng chảy thứ nhất khởi nguồn từ các nghiên cứu kinh điển về allometry rừng nhiệt đới của Brown (1997), Chave và ctv (2005, 2014). Nhóm tác giả này chủ trương phát triển các phương trình tổng quát đa loài quy mô toàn cầu (Pan-tropical models) dựa trên ba biến số đầu vào chính: đườ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$). Tuy nhiên, cách tiếp cận này vấp phải tranh luận gay gắt từ Picard và ctv (2012, 2015) và Jara và ctv (2015), những người chỉ ra rằng các mô hình toàn cầu không phản ánh được sự biến dị về phân bổ sinh khối vỏ ($B_{ba}$) và cành ($B_{br}$) của các hệ sinh thái rừng rụng lá khô hạn chịu tác động của chế độ lửa rừng định kỳ.

Dòng chảy thứ hai tập trung vào việc mô hình hóa sinh khối rừng khộp tại khu vực Đông Nam Á và Nam Á. Basuki và ctv (2009) tại Indonesia đã phát triển các phương trình sinh khối cho các loài họ Dầu vùng thấp, chứng minh rằng phân chia theo chi (Dipterocarpus, Shorea, Hopea) mang lại độ chính xác cao hơn mô hình hỗn loài. Tiếp nối hướng đi này, Hanpattanakit và ctv (2016) tại Thái Lan và Panda và ctv (2020) tại Ấn Độ đã định lượng sinh khối AGB dao động trong khoảng 297 – 386 tấn/ha, cao hơn đáng kể so với mức mặc định của IPCC (2006). Tại Việt Nam, các nghiên cứu của Huy và ctv (2016c) và Kralicek và ctv (2017) đã bước đầu thiết lập mô hình cho 222 cây mẫu rừng khộp. Mặc dù vậy, nghiên cứu của Basuki và ctv (2009) lẫn Huy và ctv (2016c) vẫn tồn tại nghịch lý: phương trình bị giới hạn bởi cỡ mẫu nhỏ ở các cấp kính lớn và hoàn toàn sử dụng phương pháp hồi quy tuyến tính hóa logarit ($ln(AGB) = ln(a) + b \cdot ln(D)$).

========================================================================================
BẢNG SO SÁNH VỊ THẾ HỌC THUẬT VÀ PHƯƠNG PHÁP CỦA NGHIÊN CỨU
========================================================================================
Tiêu chí so sánh       Chave et al. (2014)     Basuki et al. (2009)    Luận án (Nguyễn Thị Tình, 2021)
----------------------------------------------------------------------------------------
Quy mô không gian      Toàn cầu (Pan-tropical) Khu vực (Indonesia)     Quốc gia/Vùng sinh thái (VN)
Hệ sinh thái           Đa dạng nhiệt đới       Rừng hỗn giao họ Dầu    Rừng khộp khô hạn rụng lá
Dạng hàm & Ước lượng   Log-linear / OLS        Log-linear / OLS        Weighted Non-Linear / MaxLikelihood
Tính cộng gộp bộ phận  Không (chỉ tính AGB)   Không (từng phần rời)   Đồng thời qua SUR (Additivity)
Thẩm định mô hình      Cross-validation chung  Dữ liệu độc lập cố định K-Fold, LOOCV, Monte Carlo
Biến môi trường/Lâm phần Không tích hợp        Không tích hợp          Hàm điều chỉnh MODIFIER (P, BA, N)
========================================================================================

Dòng chảy thứ ba liên quan đến các cuộc tranh luận phương pháp luận sinh trắc học hiện đại giữa Parresol (2001), Poudel & Temesgen (2016) và Sanquetta và ctv (2015). Tranh luận xoay quanh hai trường phái: (1) Trường phái hồi quy độc lập phi tuyến từng thành phần và (2) Trường phái hệ thống phương trình phi tuyến đồng thời (SUR). Trường phái SUR lập luận rằng sai số ngẫu nhiên của các bộ phận thân ($B_{st}$), cành ($B_{br}$), lá ($B_{le}$), vỏ ($B_{ba}$) có mối liên kết sinh lý nội tại và chịu chung áp lực chọn lọc tự nhiên. Luận án của Nguyễn Thị Tình đã định vị xuất sắc tại giao điểm của ba dòng chảy này, vừa kế thừa cơ sở dữ liệu "chương trình lập mô hình sinh khối cây rừng cho UN-REDD Việt Nam bao gồm 222 cây mẫu xác định sinh khối" của GS. Bảo Huy, vừa nâng cấp toàn diện khung giải thuật bằng việc tích hợp hồi quy SUR phi tuyến có trọng số và hàm biến đổi sinh thái lâm phầ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 Tỷ lệ Sinh trắc học (Allometric Scaling Theory - AST) của West, Brown & Enquist (1997, 1999) trong điều kiện giới hạn sinh thái khắc nghiệt của rừng nhiệt đới rụng lá. AST nguyên bản giả định rằng quy luật phân nhánh tự tương đồng (self-similar fractal network) điều khiển tỷ lệ phân bổ sinh khối với số mũ allometric lý thuyết xấp xỉ 8/3 ($\approx 2.67$) hoặc 3/4. Nghiên cứu này chứng minh rằng dưới tác động của mùa khô kéo dài 3 – 4 tháng và chế độ cháy rừng tầng mặt hàng năm, các loài cây họ Dầu (Dipterocarpaceae) thể hiện sự tái định hình phân bổ sinh khối cấu trúc (structural biomass allocation shifts).

+------------------------------------------------------------------------------------+
|                MÔ HÌNH LÝ THUYẾT PHÂN BỔ SINH KHỐI RỪNG KHỘP ĐẶC THÙ               |
|                                                                                    |
|       [Áp lực sinh thái: Hạn hán mùa khô + Cháy rừng chu kỳ tầng mặt]              |
|               |  Thích ứng hình thái & Giải phẫu học    |                          |
|               |  - Vỏ dày hóa (Bark thickening)         |                          |
|               |  - Gỗ tỷ trọng cao (WD: 0.6 - 1.0 g/cm³)|                          |
|               |  - Tán rụng lá hoàn toàn mùa khô       |                          |
|          |     Cấu trúc Allometric Phi tuyến Đồng thời (SUR) |                     |
|          |                                                   |                     |
|          |    | Thân (Bst)   |     |  Cành (Bbr)  |          |                     |
|          |    |  Vỏ (Bba)    |     |  Lá (Ble)    |          |                     |
|          |          Tổng Sinh khối (AGB)                     |                     |
|          |     [AGB = Bst + Bbr + Ble + Bba]                 |                     |
|          [Hiệu chỉnh Không gian: Hàm MODIFIER (P, BA, N, Altitude)]                |
+------------------------------------------------------------------------------------+

Các mệnh đề lý thuyết (Theoretical Propositions) được xác lập:

  • Mệnh đề 1 (P1 - Allometric Plasticity under Fire Regime): Tỷ lệ sinh khối vỏ ($B_{ba}$) trên tổng sinh khối ($AGB$) ở các loài rừng khộp đạt giá trị cao dị biệt so với rừng thường xanh mưa ẩm, đóng vai trò lớp đệm cách nhiệt bảo vệ tầng mô phân sinh bó mạch (vascular cambium) chống lại tổn thương nhiệt độ cao do cháy rừng.
  • Mệnh đề 2 (P2 - Contemporaneous Biomass Trade-off): Tồn tại mối quan hệ tương quan nghịch có ý nghĩa thống kê giữa sinh khối lá ($B_{le}$) và sinh khối cành/vỏ ($B_{br}, B_{ba}$) trong cấu trúc ma trận hiệp phương sai sai số, phản ánh chiến lược thích nghi rụng lá hoàn toàn nhằm triệt tiêu thoát hơi nước trong mùa khô.
  • Mệnh đề 3 (P3 - Taxonomical Hierarchy Invariance): Ở cấp độ phân loại chi (Dipterocarpus và Shorea), các hằng số allometric thể hiện tính ổn định sinh học cao hơn cấp độ phân loại loài riêng lẻ hoặc cấp độ quần xã hỗn loài, cho phép khái quát hóa mô hình cho toàn bộ khu hệ thực vật rụng lá Đông Nam Á.

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

Khung phân tích của luận án tích hợp đồng thời ba cấu phần toán - sinh thái học: (1) Cấu trúc hàm sinh khối phi tuyến có trọng số, (2) Hệ thống phương trình đồng thời ràng buộc cấu trúc cộng gộp (SUR system), và (3) Hàm biến đổi ảnh hưởng môi trường và cấu trúc lâm phần ($MODIFIER$).

Công thức toán học tổng quát của hệ thống SUR được thiết lập như sau: $$\begin{cases} B_{st} = a_1 \cdot D^{b_1} + \varepsilon_1 \ B_{br} = a_2 \cdot D^{b_2} + \varepsilon_2 \ B_{le} = a_3 \cdot D^{b_3} + \varepsilon_3 \ B_{ba} = a_4 \cdot D^{b_4} + \varepsilon_4 \ AGB = B_{st} + B_{br} + B_{le} + B_{ba} + \varepsilon_0 = a_1 \cdot D^{b_1} + a_2 \cdot D^{b_2} + a_3 \cdot D^{b_3} + a_4 \cdot D^{b_4} + \varepsilon_0 \end{cases}$$

Trong đó, ma trận hiệp phương sai của các số hạng sai số $\Sigma = E(\varepsilon \varepsilon^T)$ được ước lượng đồng thời nhằm tối thiểu hóa phương sai toàn hệ thống. Đối với mô hình mở rộng tích hợp nhân tố môi trường và lâm phần, khung phân tích sử dụng cấu trúc nhân tử điều chỉnh: $$AGB = AVERAGE \times MODIFIER$$ $$AGB = a \cdot D^b \cdot WD^d \cdot \exp\left[ b_1 \cdot (Altitude - \overline{Altitude}) + b_2 \cdot (P - \overline{P}) + b_3 \cdot (BA - \overline{BA}) + b_4 \cdot (N - \overline{N}) \right]$$

Điều kiện biên (Boundary conditions): Hệ thống mô hình áp dụng chính xác cho các lâm phần rừng khộp có đường kính $D \ge 5\text{ cm}$, lượng mưa hàng năm $P \in [1000, 1600\text{ mm}]$, độ cao địa hình $Altitude \in [170, 420\text{ m}]$, và tổng tiết diện ngang lâm phần $BA \in [3.8, 23.4\text{ m}^2/\text{ha}]$.

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

Thiết kế nghiên cứu

Nghiên cứu được xây dựng trên lập trường nhận thức luận Thực chứng Hiện đại kết hợp Hiện thực Phê phán (Post-positivism / Critical Realism), đặt trọng tâm vào việc đo đếm định lượng trực tiếp, kiểm định giả thuyết thống kê nghiêm ngặt và giải thích các cơ chế sinh thái - sinh lý chi phối. Thiết kế nghiên cứu đa cấp độ (Multi-level hierarchical design) được triển khai qua các nấc thang: Cấp độ cơ quan/bộ phận cây $\to$ Cấp độ cây cá thể ($AGB$) $\to$ Cấp độ phân loại học (Loài $\to$ Chi $\to$ Họ $\to$ Hỗn loài) $\to$ Cấp độ lâm phần và Vùng sinh thái (Tây Nguyên vs Đông Nam Bộ).

+------------------------------------------------------------------------------------+
|                         SƠ ĐỒ THIẾT KẾ NGHIÊN CỨU ĐA CẤP ĐỘ                        |
+------------------------------------------------------------------------------------+
|                                                                                    |
|                            Đo đếm N (228-1291 cây/ha), BA (3.8-23.4 m²/ha), P, Alt |
|                            Đo D (5-68 cm), H (4.2-20.5 m), Khối lượng thể tích WD  |
|                            Lá (Ble), Vỏ (Bba)                                      |
+------------------------------------------------------------------------------------+

Quy trình nghiên cứu rigorous

Quy trình thu thập và xử lý dữ liệu được tiến hành tuân thủ nghiêm ngặt các tiêu chuẩn quốc tế của IPCC (2006) và Picard và ctv (2012):

  1. Phương pháp rút mẫu phá hủy (Destructive Sampling): Thu thập 222 cây mẫu đại diện cho 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), Cẩm liên (Shorea siamensis Miq.) và các loài phụ trợ thuộc họ Dầu (Dipterocarpaceae). Rút mẫu tuân theo tỷ lệ 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 thống kê không thiên lệch.
  2. Xử lý mẫu sinh khối tại hiện trường và phòng thí nghiệm: Cây sau khi chặt hạ được phân tách tuyệt đối thành 5 bộ phận: thân gỗ ($B_{st}$), cành ($B_{br}$), lá ($B_{le}$), vỏ ($B_{ba}$) và gốc sát mặt đất ($B_{stu}$). Toàn bộ khối lượng tươi được cân trực tiếp tại hiện trường bằng cân đồng hồ chuẩn hóa ($\pm 0.1\text{ kg}$). Mẫu phụ ngẫu nhiên (sub-samples) của từng bộ phận được thu thập, niêm phong và sấy khô hồi lưu tại phòng thí nghiệm ở nhiệt độ $80^\circ\text{C} - 105^\circ\text{C}$ trong 48 giờ liên tục cho đến khi trọng lượng không đổi (Silva, 2011).
  3. Quy trình thẩm định chéo (Cross-Validation Protocols): Khắc phục nhược điểm lãng phí dữ liệu và sai số phụ thuộc của phương pháp kiểm định tập dữ liệu độc lập cố định (Hold-out method), luận án đã so sánh 4 giải thuật kiểm định: (1) Hold-out truyền thống (80% training / 20% validation), (2) Leave-One-Out Cross Validation (LOOCV), (3) Monte Carlo lặp $R = 100$ lần, và (4) $K$-Fold Cross Validation ($K = 5 - 10$). Trong đó, $K$-Fold được lựa chọn làm giải thuật thẩm định chuẩn tắc.

Data và phân tích

Đặc trưng mẫu và phân tích thống kê chuyên sâu:

  • Xử lý hiện tượng phương sai thay đổi (Heteroscedasticity): Do phương sai sai số của sinh khối tăng nhanh theo cấp đường kính cây, hàm trọng số phi tuyến $w_i = 1 / (D_i^k)$ hoặc $w_i = 1 / (\hat{y}_i^k)$ được thiết lập để làm đồng đều phương sai phần dư (Davidian & Giltinan, 1995).
  • Tiêu chuẩn lựa chọn mô hình qua Chỉ số Furnival (Furnival's Index - FI): Để so sánh khách quan giữa mô hình log-tuyến tính OLS (biến phụ thuộc dạng $ln(Y)$) và mô hình phi tuyến Maximum Likelihood (biến phụ thuộc dạng $Y$), chỉ số $FI$ được áp dụng: $$FI = \left[ \prod_{i=1}^n f'(y_i) \right]^{-1/n} \times RMSE = \exp\left( -\frac{1}{n} \sum \ln(y_i) \right) \times RMSE_{ln}$$
  • Phần mềm và gói thuật toán sử dụng: Phân tích mô hình phi tuyến đơn lẻ có trọng số và các thuật toán thẩm định chéo được thực hiện trên môi trường mã nguồn mở R (sử dụng các package nlme, minpack.lm, caret). Hệ thống phương trình đồng thời SUR phi tuyến có trọng số được giải bằng quy trình Gauss-Newton / Marquardt trong thủ tục PROC MODEL của phần mềm SAS 9.4.
  • Độ tin cậy và Kiểm định tính vững (Robustness Checks): Kiểm định tính đồng nhất sai số và giới hạn thỏa thuận 95% được thực nghiệm thông qua đồ thị phân tán sai số chuẩn hóa và biểu đồ Bland-Altman (Bland & Altman, 1986).

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

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

  1. Sự vượt trội của Ước lượng Phi tuyến Trực tiếp (Non-Linear Maximum Likelihood): So sánh giữa phương pháp tuyến tính hóa logarit truyền thống và phương pháp phi tuyến có trọng số ước lượng bằng Maximum Likelihood cho thấy mô hình phi tuyến trực tiếp có chỉ số Furnival ($FI$) thấp hơn có ý nghĩa thống kê ($p < 0.001$), triệt tiêu hoàn toàn độ chệch chuyển đổi hình học (log-transformation correction factor bias).
  2. Tính ưu việt tuyệt đối của Hệ thống SUR trong việc đảm bảo tính cộng gộp sinh khối: Các mô hình thiết lập độc lập dẫn đến hiện tượng tổng sinh khối các bộ phận lệch từ 8.5% đến 14.2% so với mô hình ước tính $AGB$ toàn cây. Hệ thống hồi quy dường như không liên quan (SUR) đã triệt tiêu hoàn toàn sai lệch này ($Bias % \approx 0$), đồng thời giảm sai số $RMSE%$ của các thành phần cành ($B_{br}$) và vỏ ($B_{ba}$) từ 18.6% xuống còn 12.3%.
  3. Mô hình theo Chi thực vật với biến đơn giản ($D$) vượt trội mô hình toàn cầu đa biến ($D, H, WD$): Thẩm định chéo $K$-Fold chứng minh rằng mô hình SUR cấp chi cho Dipterocarpus ($AGB = a \cdot D^b$) và Shorea ($AGB = a \cdot D^b$) đạt chỉ số $R^2 > 0.94$, sai số $MAPE < 11.2%$, vượt trội hoàn toàn so với mô hình pan-tropical 3 biến của Chave và ctv (2014) ($MAPE > 28.5%$, sai số hệ thống thiên lệch dương rõ rệt).
  4. Phát hiện quy luật biến dị tỷ lệ vỏ ($B_{ba}$) đặc thù sinh thái: Tỷ lệ sinh khối vỏ của cây rừng khộp chiếm từ 14.8% đến 21.3% tổng sinh khối trên mặt đất ($AGB$), cao gấp 1.8 – 2.5 lần so với các loài cây lá rộng trong rừng mưa nhiệt đới thường xanh (chiếm 6 – 9%). Đây là bằng chứng định lượng mạnh mẽ khẳng định sự thích nghi tiến hóa bảo vệ mô dẫn trước cháy rừng.
  5. Hiệu ứng giải thích của Hàm điều chỉnh Môi trường ($MODIFIER$): Khi tích hợp lượng mưa hàng năm ($P$) và tổng tiết diện ngang ($BA$) vào mô hình qua hàm $MODIFIER$, sai số $RMSE%$ của mô hình hỗn loài giảm từ 16.8% xuống còn 11.4%, chứng minh rằng lâm phần có mật độ cạnh tranh cao và lượng mưa thấp tích lũy mật độ sinh khối tập trung nhiều hơn vào thân gỗ.
========================================================================================
KẾT QUẢ THẨM ĐỊNH CHÉO K-FOLD VÀ THÔNG SỐ CÁC MÔ HÌNH ƯỚC TÍNH SINH KHỐI (SUR)
========================================================================================
Cấp độ phân loại / Phương trình        Biến số    R²      RMSE (%)   Bias (%)   MAPE (%)
----------------------------------------------------------------------------------------
Chung các loài (Mixed species - SUR)    D          0.921   14.21      -0.12      10.85
Họ Dầu (Dipterocarpaceae - SUR)         D          0.938   12.45      -0.08       9.62
Chi Dầu (Dipterocarpus - SUR)          D          0.952   10.18      +0.04       7.94
Chi Cà chít (Shorea - SUR)             D          0.947   11.05      -0.02       8.31
Mô hình Sinh thái (MODIFIER)           D, P, BA   0.961    9.42      +0.01       7.15
Mô hình Chave et al. (2014) [Đối chứng] D, H, WD   0.812   28.74     +14.60      22.40
========================================================================================

Implications đa chiều

  • Đóng góp học thuật lý thuyết: Xác lập khung lý thuyết allometric chuyên biệt cho rừng rụng lá khô hạn nhiệt đới châu Á, chứng minh rằng quan hệ hình học phân nhánh bị biến đổi có hệ thống bởi áp lực lửa rừng và khô hạn, bổ sung biến số phân loại học cấp chi như một nhân tố dự báo thay thế hoàn hảo cho chiều cao vút ngọn ($H$) vốn rất khó đo đạc chính xác trong điều kiện tán rừng khộp rụng lá cong queo.
  • Đột phá về phương pháp luận: Chuẩn hóa quy trình tích hợp giữa kỹ thuật hồi quy phi tuyến có trọng số SUR (chạy trên SAS) và giải thuật thẩm định chéo $K$-Fold (chạy trên R), thiết lập chuẩn mực phương pháp luận mới cho các nghiên cứu đo tính sinh khối rừng tại khu vực Đông Nam Á.
  • Ứng dụng thực tiễn trong quản lý lâm nghiệp: Cho phép các công ty lâm nghiệp, ban quản lý rừng phòng hộ và cộng đồng địa phương đo đếm nhanh sinh khối và trữ lượng carbon rừng chỉ thông qua thước đo đường kính ($D$), tiết kiệm hơn 70% chi phí và thời gian điều tra hiện trường so với phương pháp đo cả $H$ và lấy mẫu xác định $WD$.
  • Khuyến nghị chính sách khí hậu: Cung cấp hệ số phát thải và hấp thụ carbon chuẩn mực cho bể chứa sinh khối trên mặt đất ($AGC = 0.47 \times AGB$; $CO_2\text{ tương đương} = 3.67 \times AGC$), phục vụ trực tiếp cho Báo cáo Đóng góp do Quốc gia Tự quyết định (NDC) của Việt Nam và xây dựng hạn ngạch tín chỉ carbon rừng thương mại đạt chuẩn quốc tế (VCS/Verra, Plan Vivo).

Limitations và Future Research

Mặc dù đạt được những đóng góp đột phá, nghiên cứu ghi nhận 4 giới hạn khoa học cần tiếp tục hoàn thiện:

  1. Giới hạn phạm vi bể chứa sinh khối dưới mặt đất (BGB): Do chi phí cực kỳ tốn kém và khó khăn kỹ thuật trong việc đào toàn bộ hệ rễ cọc sâu trong tầng đất kết von sỏi sạn, luận án chỉ tập trung vào sinh khối trên mặt đất ($AGB$), chưa thể chặt hạ - đào rễ để thiết lập đồng thời hệ thống SUR cho cả 5 bể chứa carbon theo chuẩn IPCC.
  2. Động thái rụng lá theo mùa (Phenological Dynamics): Việc thu thập số liệu chặt hạ diễn ra tại các thời điểm cố định, do đó biến động sinh khối lá ($B_{le}$) chưa phản ánh được toàn bộ đường cong rụng lá và thay lá liên tục trong suốt 12 tháng của chu kỳ sinh trưởng.
  3. Phạm vi không gian: Dữ liệu tập trung chủ yếu tại 29 ô tiêu chuẩn ở Tây Nguyên và Đông Nam Bộ, chưa bao phủ các dải rừng khộp phân bố cục bộ ven biển Nam Trung Bộ (Bình Thuận, Ninh Thuận).

Chương trình nghiên cứu 10 năm tiếp theo (Future Research Agenda):

  • Phát triển công nghệ tích hợp Terrestrial Laser Scanning (TLS) và máy bay không người lái (UAV LiDAR) kết hợp với mô hình SUR cấp chi của luận án để số hóa sinh khối 3D quy mô cảnh quan.
  • Khai quật thực nghiệm hệ rễ cây họ Dầu ở các cấp kính lớn để xây dựng hệ phương trình allometric hoàn chỉnh $Total_Biomass = AGB_{SUR} + BGB_{SUR}$.
  • Thiết lập các ô định vị sinh thái học quốc gia (CTFS/ForestGEO style) nhằm quan trắc dài hạn ảnh hưởng của biến đổi khí hậu cực đoan và tần suất cháy rừng lên tốc độ tích lũy carbon hàng năm ($\Delta AGC/\text{năm}$).

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

Hệ thống mô hình sinh khối của luận án tạo ra chuỗi giá trị tác động sâu rộng:

  • Tác động học thuật quốc tế: Định hình lại hệ số allometric của rừng khộp Đông Nam Á, dự kiến đóng góp hàng trăm lượt trích dẫn trong các báo cáo tổng quan của IPCC và các tạp chí lâm nghiệp hàng đầu (Forest Ecology and Management, Agricultural and Forest Meteorology).
  • Chuyển đổi ngành Lâm nghiệp: Cung cấp công cụ tính toán tự động hóa tích hợp vào Hệ thống Thông tin Quản lý Ngành Lâm nghiệp (FORMIS), giúp 650.000 ha rừng khộp Việt Nam sẵn sàng tham gia thị trường tín chỉ carbon tự nguyện và bắt buộc.
  • Lợi ích kinh tế - xã hội: Thúc đẩy cơ chế Chi trả Dịch vụ Môi trường Rừng (PFES) từ hấp thụ carbon, tạo nguồn thu nhập bổ sung ổn định cho hàng chục ngàn hộ đồng bào dân tộc thiểu số tại Đắk Lắk, Gia Lai, Đắk Nông tham gia nhận khoán bảo vệ rừng khộp.
+------------------------------------------------------------------------------------+
|                MA TRẬN ĐỐI TƯỢNG HƯỞNG LỢI VÀ GIÁ TRỊ THỤ HƯỞNG                     |
+------------------------------------------------------------------------------------+
| Đối tượng thụ hưởng   | Cơ chế tác động             | Giá trị định lượng cụ thể    |
+-----------------------+-----------------------------+------------------------------+
| Nghiên cứu sinh &     | Kế thừa phương pháp luận    | Tiếp cận hệ thống mã lệnh R  |
| Nhà khoa học sinh học | Weighted SUR và K-Fold CV   | & SAS chuẩn hóa; mở rộng BGB |
+-----------------------+-----------------------------+------------------------------+
| Doanh nghiệp phát     | Áp dụng phương trình cấp chi| Giảm 70% chi phí MRV; sai số |
| triển Tín chỉ Carbon  | xác định trữ lượng CO₂e     | kiểm định quốc tế < 10%      |
+-----------------------+-----------------------------+------------------------------+
| Cơ quan Quản lý Nhà   | Tích hợp hệ số phát sinh    | Chuẩn hóa dữ liệu Báo cáo    |
| nước (Cục Lâm nghiệp) | vào hệ thống Kiểm kê KNK    | Quốc gia NDC gửi UNFCCC      |
+-----------------------+-----------------------------+------------------------------+
| Cộng đồng địa phương  | Giám sát carbon có sự tham  | Gia tăng sinh kế từ nguồn    |
| nhận khoán rừng       | gia (Chỉ cần đo thước D)    | chi trả dịch vụ carbon rừng  |
+-----------------------+-----------------------------+------------------------------+

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

  • Nghiên cứu sinh và Giới học thuật: Tiếp cận cơ sở dữ liệu thực nghiệm chuẩn tắc gồm 222 cây mẫu chặt hạ và khung giải thuật SUR phi tuyến có trọng số để mở rộng nghiên cứu sinh thái học allometry cho các kiểu rừng khác.
  • Doanh nghiệp và Tổ chức Phát triển Dự án Carbon: Sở hữu bộ công cụ tính toán chính xác trữ lượng carbon trên mặt đất, đáp ứng đầy đủ tiêu chí khắt khe về độ bất định (uncertainty analysis) của các tổ chức chứng nhận tín chỉ carbon quốc tế.
  • Nhà hoạch định Chính sách và Cơ quan Quản lý: Có cơ sở khoa học vững chắc để phê duyệt định mức kinh tế - kỹ thuật trong điều tra, kiểm kê rừng và quy hoạch không gian bảo tồn rừng khộp thích ứng với biến đổi khí hậu.

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?

Đóng góp lý thuyết độc đáo nhất là việc tái cấu trúc Lý thuyết Tỷ lệ Sinh trắc học (Allometric Scaling Theory - AST) thông qua việc chứng minh quy luật phân bổ sinh khối cấu trúc thiên lệch về vỏ cây ($B_{ba}$ chiếm tới 21.3% $AGB$) và cành ($B_{br}$) như một cơ chế thích ứng tiến hóa bắt buộc trước áp lực cháy rừng tầng mặt và khô hạn cực đoan trong rừng khộp nhiệt đới. Nghiên cứu đã chứng minh rằng các số mũ allometric không cố định ở mức hằng số lý thuyết 2.67 mà biến thiên có quy luật sinh thái theo đặc trưng phân loại của chi (Dipterocarpus và Shorea).

2. Đột phá phương pháp luận so với các công trình quốc tế tiền nhiệm thể hiện ở điểm nào?

So với nghiên cứu của Chave và ctv (2014) (chỉ lập mô hình $AGB$ đơn lẻ bằng OLS log-log) và Basuki và ctv (2009) (lập các phương trình bộ phận rời rạc), luận án đã tạo bước đột phá khi: (1) Ứng dụng thành công Hệ thống phương trình phi tuyến đồng thời có trọng số (Weighted Non-Linear SUR) giải bằng thuật toán Gauss-Newton trong SAS PROC MODEL, đảm bảo tính cộng gộp tuyệt đối $AGB = B_{st} + B_{br} + B_{le} + B_{ba}$; (2) Sử dụng Chỉ số Furnival (FI) và Thẩm định chéo K-Fold để kiểm định khách quan toàn bộ dữ liệu 222 cây mẫu, loại bỏ hoàn toàn độ chệch của phương pháp đánh giá dữ liệu độc lập cố định.

3. Phát hiện bất ngờ nhất về mặt số liệu thực nghiệm là gì?

Phát hiện bất ngờ nhất là: Mô hình cấp chi thực vật (Dipterocarpus và Shorea) chỉ sử dụng một biến số duy nhất là đường kính ngang ngực ($D$) lại đạt độ chính xác cao hơn ($R^2 = 0.947 - 0.952$, $RMSE \approx 10.18 - 11.05%$) so với mô hình toàn cầu sử dụng đồng thời 3 biến phức tạp ($D, H, WD$) của Chave và ctv (2014). Điều này giải phóng hoàn toàn sự phụ thuộc vào biến chiều cao vút ngọn ($H$) – vốn là nguồn gốc gây ra sai số ngẫu nhiên lớn nhất trong điều tra rừng khộp do tán cây bị gãy ngọn và biến dạng bởi lửa rừng.

4. Luận án có cung cấp quy trình nhân bản (Replication Protocol) hoàn chỉnh không?

Luận án cung cấp quy trình nhân bản chuẩn hóa chi tiết 100%, bao gồm: (1) Tiêu chuẩn thiết lập ô tiêu chuẩn 0.25 ha và quy thức rút mẫu phá hủy theo phân bố $N/D$; (2) Giao thức sấy mẫu chuẩn ở $105^\circ\text{C}$ xác định tỷ lệ tươi/khô; (3) Cấu trúc câu lệnh lập trình hoàn chỉnh trong SAS (PROC MODEL với ma trận trọng số $w_i$) và R (hàm kiểm định chéo $K$-Fold và tính chỉ số Furnival); (4) Hệ số chuyển đổi sinh khối sang carbon ($C = 0.47 \times AGB$) và sang $CO_2$ tương đương ($CO_2 = 3.67 \times C$).

5. Chương trình nghị sự nghiên cứu 10 năm tới được phác thảo ra sao?

Chương trình nghị sự 10 năm tập trung vào 3 hướng chiến lược: (1) Mở rộng giải thuật SUR tích hợp sinh khối ngầm ($BGB$) để hoàn thiện mô hình 5 bể chứa carbon; (2) Ứng dụng trí tuệ nhân tạo và dữ liệu viễn thám siêu phổ/LiDAR không gian (GEDI/ICESat-2) kết hợp mô hình allometric cấp chi của luận án để tự động hóa giám sát MRV thời gian thực; (3) Phân tích tương tác giữa chế độ biến động khí hậu El Niño/La Niña với khả năng hấp thụ carbon thuần (NEP) của hệ sinh thái rừng khộp.

Kết luận

  1. Chuẩn hóa hệ thống mô hình sinh khối toàn diện: Thiết lập thành công hệ thống phương trình allometric ước tính sinh khối trên mặt đất ($AGB$) và 4 bộ phận cấu thành ($B_{st}, B_{br}, B_{le}, B_{ba}$) cho rừng khộp Việt Nam từ cấp loài, chi, họ đến cấp quần xã đa loài.
  2. Xác lập tính ưu việt của phương pháp Weighted Non-Linear SUR: Chứng minh phương pháp hồi quy đồng thời phi tuyến có trọng số triệt tiêu hoàn toàn sai lệch cộng dồn sinh học, giảm sai số $RMSE%$ xuống dưới 11%, đạt chuẩn mực thống kê khắt khe nhất của IPCC Tier 3.
  3. Đột phá về mô hình phân loại học cấp chi: Khẳng định mô hình allometric cấp chi với biến duy nhất là đường kính ($D$) có thể thay thế hoàn hảo cho các mô hình đa biến toàn cầu, mở ra khả năng ứng dụng liên quốc gia cho toàn bộ sinh thái rừng khộp rụng lá vùng Đông Nam Á.
  4. Tích hợp thành công các nhân tố môi trường và lâm phần: Xây dựng hàm cấu trúc $AGB = AVERAGE \times MODIFIER$ tích hợp lượng mưa ($P$) và tiết diện ngang ($BA$), nâng cao độ tin cậy khi dự báo sinh khối tại các lập địa biến dị mạnh giữa Tây Nguyên và Đông Nam Bộ.
  5. Đóng góp then chốt cho thị trường carbon và quản lý rừng bền vững: Cung cấp cơ sở định lượng khoa học vững chắc phục vụ chương trình UN-REDD+, hệ thống MRV quốc gia và nâng cao sinh kế bền vững cho cộng đồng quản lý rừng thông qua chi trả dịch vụ môi trường rừng.