Tổng quan về luận án

Biến đổi khí hậu toàn cầu và sự gia tăng nồng độ $\text{CO}_2$ trong khí quyển từ mức 280 ppm thời kỳ tiền công nghiệp lên vượt ngưỡng 390 ppm đã đặt hệ sinh thái rừng vào vị trí trung tâm của các cam kết quốc tế về giảm phát thải theo Công ước Khung của Liên Hợp Quốc về Biến đổi Khí hậu (UNFCCC, 1992) và Nghị định thư Kyoto (1997). Rừng nhiệt đới lưu giữ khoảng 46% trữ lượng carbon trên mặt đất và 11,6% carbon dưới mặt đất toàn cầu (Brown et al., 1989; Soepadmo, 1993). Tuy nhiên, tại các vùng sinh thái đặc thù như vùng khô hạn và bán khô hạn Nam Trung Bộ Việt Nam, dữ liệu định lượng về sinh khối và trữ lượng carbon còn thiếu hụt trầm trọng. Luận án tiến sĩ khoa học lâm nghiệp của tác giả Trịnh Minh Hoàng với đề tài "Nghiên cứu khả năng tích tụ carbon của rừng tự nhiên ở tỉnh Ninh Thuận" (Chuyên ngành Lâm sinh, Mã số 62 62 02 05; Người hướng dẫn: PGS. Viên Ngọc Nam, PGS. Phạm Văn Hiền; Trường Đại học Nông Lâm TP. Hồ Chí Minh, 2016) là công trình tiên phong giải quyết khoảng trống học thuật này.

Khoảng trống nghiên cứu (research gap) cốt lõi được xác định là: các hướng dẫn cấp quốc tế của IPCC (2003, 2006) và AFOLU chỉ cung cấp các hệ số chuyển đổi mặc định mang tính khái quát cao (Tier 1/Tier 2), gây sai số tích lũy lên tới 20–40% khi áp dụng cho các kiểu rừng nhiệt đới khô hạn có cấu trúc phân tầng phức tạp (Brown, 2002; Chave et al., 2005). Tại Ninh Thuận, diện tích rừng tự nhiên đạt 133.101 ha, trong đó hai kiểu rừng chính là rừng kín thường xanh hơi khô nhiệt đới ($\text{R}{\text{kx}}$) và rừng thưa nửa thường xanh hơi khô nhiệt đới ($\text{R}{\text{tr}}$) chiếm tới 94.635 ha (tập trung tại Bác Ái 42.080 ha và Ninh Sơn 26.000 ha), nhưng hoàn toàn chưa có hệ thống phương trình sinh khối sinh trắc học cục bộ (allometric biomass equations), hệ số mở rộng sinh khối tầng thứ cấp ($\text{BEF}_i$), hay tỷ lệ phân chia sinh khối ($\text{R}_i$).

Nghiên cứu vận hành trên hệ thống câu hỏi và giả thuyết khoa học chặt chẽ:

  • Câu hỏi nghiên cứu 1 (RQ1): Quy luật cấu trúc lâm phần và tổ thành loài ưu thế ảnh hưởng như thế nào đến sự tích lũy sinh khối cá thể và quần thụ tại $\text{R}{\text{kx}}$ và $\text{R}{\text{tr}}$?
  • Câu hỏi nghiên cứu 2 (RQ2): Dạng phương trình toán học phi tuyến nào (lũy thừa tương quan hay Korsun–Strand) mô tả chính xác nhất mối quan hệ giữa sinh khối các bộ phận (thân, cành, lá) với đường kính ngang ngực ($D$), chiều cao vút ngọn ($H$), tiết diện ngang ($G$) và trữ lượng ($M$)?
  • Câu hỏi nghiên cứu 3 (RQ3): Trữ lượng carbon tích tụ trên mặt đất biến động như thế nào qua các trạng thái phục hồi sinh thái (IIB, IIIA1, IIIA2, IIIA3) tại Vườn Quốc gia Phước Bình?
  • Giả thuyết 1 (H1): Phương trình allometric phi tuyến đa biến kết hợp biến đường kính ($D$) và chiều cao ($H$) không nhất thiết cho độ chính xác vượt trội hơn phương trình đơn biến ($D$) khi xét đến sai số đo đạc tán rừng nhiệt đới ngoài thực địa.
  • Giả thuyết 2 (H2): Hệ số điều chỉnh sinh khối ($\text{BEF}_i$) và hệ số chuyển đổi ($\text{BCEF}_i$) không phải là hằng số cố định mà biến thiên nghịch biến phi tuyến theo cấp kính ($D$) và thể tích thân cây ($V$).

Khung lý thuyết của luận án tích hợp Lý thuyết Tương quan Sinh trưởng Allometric (Allometric Scaling Theory - WBE) và Khung hướng dẫn Thực hành Tốt về Sử dụng Đất, Thay đổi Sử dụng Đất và Lâm nghiệp của Ban Liên chính phủ về Biến đổi Khí hậu (IPCC GPG-LULUCF/AFOLU). Phạm vi không gian nghiên cứu bao quát 35 ô tiêu chuẩn định vị (kích thước 0,2–1,0 ha) và 88 cây mẫu đại diện được chặt hạ giải tích cân đo trực tiếp (47 cây $\text{R}{\text{kx}}$ và 41 cây $\text{R}{\text{tr}}$) tại Vườn Quốc gia Phước Bình. Đóng góp đột phá được lượng hóa của luận án là xác định chính xác tổng sinh khối trên mặt đất ($\text{AGB}$) và trữ lượng carbon tương ứng: $\text{R}{\text{kx}}$ đạt $87,5 \text{ tấn/ha}$ ($41,1 \text{ tấn C/ha}$); $\text{R}{\text{tr}}$ đạt $57,0 \text{ tấn/ha}$ ($26,8 \text{ tấn C/ha}$). Tổng lượng carbon lưu giữ trên mặt đất tại Vườn Quốc gia Phước Bình được định lượng đạt $114,5 \times 10^3 \text{ tấn C}$ đối với $\text{R}{\text{kx}}$ và $40,3 \times 10^3 \text{ tấn C}$ đối với $\text{R}{\text{tr}}$, cung cấp cơ sở dữ liệu Tier 3 chuẩn xác cho cơ chế Chi trả Dịch vụ Môi trường Rừng (PFES) và các sáng kiến REDD+.

Literature Review và Positioning

Tổng quan tài liệu quốc tế và trong nước về lượng hóa sinh khối và carbon rừng được phân hóa thành ba dòng nghiên cứu chính:

  1. Dòng phương pháp giải tích chặt hạ trực tiếp (Destructive Sampling): Bắt đầu từ các nghiên cứu nền tảng của Ebermeyer (1876), Satoo (1982) qua Chương trình Sinh học Quốc tế (IBP, 1960), đến Brown (1997, 2002). Đây là tiêu chuẩn vàng về độ chính xác nhưng bị giới hạn bởi chi phí, nhân lực và tính phá hủy sinh cảnh đặc dụng.
  2. Dòng mô hình hóa allometric đa biến: Ketterings et al. (2001) tại Indonesia, Chave et al. (2005) trên phạm vi liên lục địa, Zianis et al. (2005) với 607 mô hình tại châu Âu, và Návar (2009) tại Mexico. Trọng tâm của dòng này là tối ưu hóa hàm toán học ($82,6%$ dạng biến đổi logarit, $12,0%$ phi tuyến) dựa trên đường kính ($D$), chiều cao ($H$), mật độ gỗ ($\rho$) và hình số ($f$).
  3. Dòng chuyển đổi sinh khối qua hệ số $\text{BEF}/\text{BCEF}$ và tỷ lệ $\text{R}_i$: Jenkins et al. (2003) tại Hoa Kỳ, Lehtonen et al. (2004) tại Phần Lan, Jalkanen et al. (2005) và Repola (2009) tại Thụy Điển, Teobaldelli et al. (2009) ở quy mô toàn cầu.

Trong y văn tồn tại hai cuộc tranh luận học thuật sâu sắc:

  • Tranh luận 1: Số lượng biến dự đoán trong hàm allometric: Chave et al. (2005) và Návar et al. (2002) lập luận rằng tích hợp $D$, $H$ và $\rho$ vào mô hình ($B = a \cdot (D^2 H \rho)^b$) giảm thiểu căn bản sai số dư (residual error). Ngược lại, Ketterings et al. (2001) và Brown (2002) chỉ ra rằng việc đo đạc chiều cao vút ngọn ($H$) trong tán rừng nhiệt đới rậm rạp thường phát sinh sai số chủ quan từ 10–25%, biến $H$ thành nguồn gây nhiễu thay vì cải thiện độ chính xác. Do đó, hàm tương quan đơn biến lũy thừa dựa vào $D$ ($B = a \cdot D^b$) có tính khả thi và độ tin cậy thực tiễn cao hơn.
  • Tranh luận 2: Tính tĩnh hay động của hệ số $\text{BEF}$: Brown (1997) và các báo cáo ban đầu của IPCC thường gán $\text{BEF}$ là hằng số trung bình theo kiểu rừng. Tuy nhiên, Lehtonen et al. (2004) và Jenkins et al. (2003) phản biện rằng $\text{BEF}$ thay đổi mạnh theo cấp tuổi và kích thước thân cây do quá trình đào thải tự nhiên và biến đổi tỷ lệ phân bổ sinh khối thân/cành lá.

Luận án định vị chính xác vị thế học thuật bằng việc đối chiếu với các công trình quốc tế điển hình trong khu vực Đông Nam Á:

  • So sánh với Ketterings et al. (2001) tại Sumatra (Indonesia) trên rừng thứ sinh hỗn giao (khối lượng thể tích gỗ $\rho = 0,60 \text{ g/cm}^3$, dao động $0,53–0,67 \text{ g/cm}^3$): Luận án không chỉ dừng lại ở việc ước lượng sinh khối thân cây mà phân giải chi tiết tới từng cấu phần lá, cành và thân riêng biệt thông qua hệ phương trình Korsun–Strand.
  • So sánh với Basuki et al. (2009) trên 122 cây họ Dầu (Dipterocarpaceae) ở Đông Kalimantan và Chaiyo et al. (2011) trên rừng rụng lá khô ở Ratchaburi (Thái Lan): Nghiên cứu tại Ninh Thuận khai thác đặc thù của kiểu rừng nhiệt đới bán khô hạn – nơi hình thái cây có tỷ lệ phân cành thấp và tán rộng – từ đó chứng minh sự vượt trội của hàm Korsun–Strand đối với sinh khối tán lá so với mô hình lũy thừa truyền thống của Ogawa et al. (1965).
  • So sánh với các nghiên cứu trong nước của Viên Ngọc Nam (1998) trên rừng ngập mặn Cần Giờ, Bảo Huy (2010, 2012) tại Tây Nguyên, và Vũ Tiến Hinh (2012): Luận án hoàn thiện khoảng trống cho hệ sinh thái rừng khô hạn ven biển Nam Trung Bộ, thiết lập bộ tham số thực nghiệm hoàn chỉnh cho cả cấp độ cây cá thể lẫn cấp độ quần thụ lâm phần.

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

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

Nghiên cứu mở rộng và kiểm chứng thực nghiệm Lý thuyết Tương quan Sinh trưởng Allometric (WBE Model của West, Brown, & Enquist, 1997) và Lý thuyết Mô hình Đường ống (Pipe Model Theory của Shinozaki et al., 1964) trong điều kiện sinh thái khắc nghiệt của rừng tự nhiên bán khô hạn:

  1. Mở rộng quy luật phân bổ sinh khối: Chứng minh rằng quy luật phân chia vật chất quang hợp giữa cơ quan nâng đỡ (thân gỗ) và cơ quan đồng hóa (lá, cành) có sự chuyển dịch phi tuyến tính theo đường kính $D$. Tỷ lệ sinh khối cành lá ($\text{R}{\text{CL}}$) và tổng sinh khối ($\text{R}{\text{To}}$) so với thân cây tăng dần theo cấp kính theo hàm logarit $\text{R}_i = a + b \cdot \ln(D)$ và hàm lũy thừa $\text{R}_i = a \cdot D^b$.
  2. Phá vỡ định kiến về hệ số $\text{BEF}$ tĩnh: Bác bỏ việc áp dụng $\text{BEF}$ hằng số của IPCC; chứng minh hệ số điều chỉnh sinh khối giảm đơn điệu theo đường kính cây cá thể theo mô hình $\text{BEF}_i = (a + b/D)^2$ và $\text{BEF}_i = a \cdot D^{-b}$.
  3. Định chế hóa mô hình Korsun–Strand trong sinh thái học: Xác lập vị trí của hàm phân bố Korsun–Strand như một quy luật lý thuyết mô tả chính xác sự sinh trưởng có giới hạn trần sinh lý của bộ phận cành và lá cây rừng nhiệt đới.
                  [ Khung Phân Tích Allometric Đa Cấp ]
[ Cấp Cây Cá Thể (n = 88) ]                               [ Cấp Quần Thụ (35 OTC) ]
       [ Trữ Lượng Sinh Khối (AGB) & Bể Carbon Rừng Tự Nhiên ]

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

Khung phân tích của luận án tích hợp đồng thời 3 hệ lý thuyết: (1) Lý thuyết cấu trúc lâm phần sinh thái học rừng nhiệt đới; (2) Lý thuyết toán thống kê sinh trắc học phi tuyến; (3) Khung kiểm kê carbon AFOLU/IPCC.

Tính mới trong tiếp cận phân tích: Luận án không gộp chung các cấu phần cây gỗ mà tách rời thành 4 hệ phương trình độc lập nhưng tương thích sinh học: hàm sinh khối thân ($\text{B}{\text{T}}$), hàm sinh khối cành ($\text{B}{\text{C}}$), hàm sinh khối lá ($\text{B}{\text{L}}$), và hàm tổng sinh khối trên mặt đất ($\text{B}{\text{To}}$). Đồng thời, luận án xây dựng cầu nối chuyển đổi trực tiếp giữa dữ liệu điều tra lâm nghiệp truyền thống (Tiết diện ngang $G$, Trữ lượng $M$) sang Trữ lượng Sinh khối ($\text{M}_{\text{B}}$) thông qua hàm phân tích:

$$\text{M}_{\text{B}} = f(G, M) = \text{Exp}\left(a + \frac{b}{G} + \frac{c}{G^2}\right) \quad \text{hoặc} \quad \text{Exp}\left(a + \frac{b}{M} + \frac{c}{M^2}\right)$$

Điều kiện biên (Boundary conditions): Khung phân tích được xác lập giới hạn ứng dụng chuẩn xác cho các lâm phần rừng kín thường xanh hơi khô và rừng thưa nửa thường xanh nhiệt đới có đường kính $D \ge 5\text{ cm}$ (đối với cây cá thể) và độ che phủ tán $> 0,3$ (đối với quần thụ), phân bố ở đai cao dưới 1.000 m thuộc vùng khí hậu nhiệt đới gió mùa khô hạn Nam Trung Bộ.

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

Thiết kế nghiên cứu

  • Triết lý nghiên cứu (Research Philosophy): Thực chứng luận (Positivism) kết hợp Hiện thực phê phán (Critical Realism), dựa hoàn toàn trên dữ liệu cân đo thực nghiệm khách quan và quy nạp toán thống kê.
  • Thiết kế nghiên cứu: Thiết kế thứ bậc đa cấp (Multi-level hierarchical design): Cấp độ giải phẫu mô học (mẫu khô phòng thí nghiệm) $\rightarrow$ Cấp độ cây cá thể (88 cây giải tích) $\rightarrow$ Cấp độ ô tiêu chuẩn lâm phần (35 ô mẫu) $\rightarrow$ Cấp độ cảnh quan sinh thái toàn vùng (Vườn Quốc gia Phước Bình).
  • Cỡ mẫu và tiêu chí chọn mẫu:
    • Cỡ mẫu ô tiêu chuẩn: 35 ô tiêu chuẩn tạm thời đại diện có diện tích từ 0,2 ha ($20 \times 100\text{ m}$) đến 1,0 ha ($100 \times 100\text{ m}$), bao phủ đầy đủ các trạng thái phục hồi: IIB (rừng phục hồi sau nương rẫy), IIIA1 (rừng nghèo), IIIA2 (rừng trung bình), IIIA3 (rừng giàu).
    • Cỡ mẫu cây chặt hạ: 88 cây mẫu thuộc tầng ưu thế và đồng ưu thế sinh thái (47 cây thuộc $\text{R}{\text{kx}}$ và 41 cây thuộc $\text{R}{\text{tr}}$). Tiêu chí chọn cây: cây sinh trưởng bình thường, không cụt ngọn, không sâu bệnh, đại diện cho phân bố đường kính từ cấp nhỏ ($D < 10\text{ cm}$) đến cây gỗ lớn ($D > 50\text{ cm}$).

Quy trình nghiên cứu rigorous

Quy trình thu thập và xử lý số liệu tuân thủ nghiêm ngặt các giao thức quốc tế của FAO/IPCC:

  1. Đo đạc trắc lượng lâm phần: Xác định tọa độ GPS, độ dốc, hướng phơi, loại đất. Đo đếm toàn bộ cây gỗ có $D_{1.3} \ge 5,0\text{ cm}$, đo chiều cao vút ngọn ($H$), chiều cao dưới cành ($H_{\text{dc}}$), đường kính tán ($D_{\text{T}}$).
  2. Giao thức chặt hạ và phân đoạn sinh khối: Cây sau khi hạ được bóc tách riêng 4 bộ phận: thân gỗ, cành lớn/nhỏ, lá và hoa quả. Thân cây được cắt đoạn giải tích 1–2 m, cân đo khối lượng tươi ($\text{B}_{\text{tươi}}$) ngay tại hiện trường với độ chính xác $\pm 0,1\text{ kg}$. Lấy các thớt gỗ giải tích dày 3–5 cm tại gốc, $1,3\text{ m}$, giữa thân và ngọn để xác định tỷ trọng và độ ẩm.
  3. Quy trình sấy khô và phân tích carbon: Rút mẫu phụ (sub-sampling) 0,5–1,0 kg cho từng bộ phận (thân, cành, lá), sấy hồi nhiệt tại phòng thí nghiệm ở nhiệt độ $105^\circ\text{C}$ đến khi khối lượng không đổi để xác định tỷ lệ chất khô. Phân tích hàm lượng carbon hữu cơ ($\text{PC}%$) trong sinh khối khô bằng phương pháp oxy hóa hoặc sử dụng hệ số chuẩn hóa theo tiêu chuẩn IPCC.
  4. Tam giác đạc (Triangulation): Kiểm chứng chéo giữa số liệu cân đo giải tích trực tiếp, số liệu tính toán qua thể tích $\times$ khối lượng thể tích gỗ ($\text{V} \times \rho$), và mô hình tương quan allometric. Độ tin cậy của mô hình được bảo đảm với hệ số tương quan $R > 0,85$.

Data và phân tích

  • Đặc trưng thống kê mẫu nghiên cứu: Tỷ trọng gỗ khô ($\rho$) của các loài cây ưu thế dao động từ $0,48\text{ g/cm}^3$ đến $0,72\text{ g/cm}^3$, trung bình đạt $0,57\text{ g/cm}^3$ (trùng khớp với giá trị trung bình cây gỗ nhiệt đới châu Á theo Reyes et al., 1992). Tỷ lệ hàm lượng carbon trong chất khô được chuẩn hóa ở mức $\text{PC} = 47,0%$ (theo hướng dẫn của IPCC, 2003, 2006).
  • Kỹ thuật phân tích phi tuyến: Ước lượng tham số bằng phương pháp hồi quy phi tuyến tính (Nonlinear Regression) sử dụng thuật toán lặp Gauss–Newton và Levenberg–Marquardt.
  • Kiểm định độ vững (Robustness Checks): So sánh và lựa chọn mô hình tối ưu dựa trên 4 chỉ số thống kê định lượng:
    1. Tổng sai lệch bình phương nhỏ nhất ($\text{SSR} = \sum (y_i - \hat{y}_i)^2 \rightarrow \text{Min}$).
    2. Hệ số xác định hiệu chỉnh ($R^2 \rightarrow \text{Max}$).
    3. Sai số chuẩn của ước lượng ($S_e$).
    4. Sai số tuyệt đối phần trăm trung bình ($\text{MAPE}% = \frac{1}{n}\sum \left|\frac{y_i - \hat{y}_i}{y_i}\right| \times 100% \rightarrow \text{Min}$).
  • So sánh đánh giá 6 phương pháp ước lượng sinh khối cây cá thể và 2 phương pháp ước lượng sinh khối quần thụ lâm phần nhằm loại trừ triệt để sai số phương pháp luận.

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

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

Thứ nhất, quy luật tương quan hình thái và sự thích hợp của mô hình allometric:

  • Đối với rừng kín thường xanh hơi khô ($\text{R}{\text{kx}}$), hàm lũy thừa biến đơn $D$ là mô hình tối ưu nhất để ước lượng tổng sinh khối trên mặt đất ($\text{B}{\text{To}}$) và sinh khối thân ($\text{B}_{\text{T}}$):

$$\text{B}_{\text{To}} = a \cdot D^b \quad (R^2 > 0,92; \text{MAPE} < 12,5%)$$

  • Tuy nhiên, đối với sinh khối cành ($\text{B}{\text{C}}$) và sinh khối lá ($\text{B}{\text{L}}$) của $\text{R}_{\text{kx}}$, hàm phân bố Korsun–Strand thể hiện sự vượt trội tuyệt đối về mặt thống kê:

$$\text{B}_i = \text{Exp}\left(a + \frac{b}{D} + \frac{c}{D^2}\right)$$

Hàm Korsun–Strand phản ánh đúng giới hạn sinh lý khi cây đạt cấp kính lớn thì tốc độ tích lũy sinh khối tán lá chậm lại.

  • Đối với rừng thưa nửa thường xanh hơi khô ($\text{R}{\text{tr}}$), hàm lũy thừa tương quan $\text{B}i = a \cdot D^b$ thích hợp cho toàn bộ các thành phần sinh khối ($\text{B}{\text{To}}, \text{B}{\text{T}}, \text{B}{\text{C}}, \text{B}{\text{L}}$) do cấu trúc tán đơn giản và mật độ thưa.

Thứ hai, phát hiện về động thái của hệ số $\text{BEF}_i$ và tỷ lệ $\text{R}_i$:

  • Trích dẫn trực tiếp từ dữ liệu luận án: "Những hệ số điều chỉnh sinh khối ($\text{BEF}i$) đối với những thành phần cây gỗ ở $\text{R}{\text{kx}}$ và $\text{R}_{\text{tr}}$ giảm dần theo sự gia tăng $D$ tương ứng với hàm $\text{BEF}_i = (a + b/D)^2$ và $\text{BEF}_i = a \cdot D^{-b}$".
  • Ngược lại, tỷ lệ tổng sinh khối và tỷ lệ cành lá so với thân cây tăng dần theo cấp kính: "Tỷ lệ tổng sinh khối trên mặt đất ($\text{R}{\text{To}}$) và tỷ lệ sinh khối cành và lá ($\text{R}{\text{CL}}$) so với sinh khối thân cây gỗ ở $\text{R}{\text{kx}}$ và $\text{R}{\text{tr}}$ gia tăng dần theo sự gia tăng $D$ tương ứng với hàm $\text{R}_i = a + b \cdot \ln(D)$ và $\text{R}_i = a \cdot D^b$".

Thứ ba, mô hình sinh khối cấp quần thụ lâm phần: Hàm Korsun–Strand với biến dự đoán là Tiết diện ngang ($G, \text{m}^2/\text{ha}$) hoặc Trữ lượng lâm phần ($M, \text{m}^3/\text{ha}$) là mô hình chính xác nhất để ước lượng trực tiếp sinh khối quần thụ mà không cần tổng hợp qua từng cây cá thể ($R^2 > 0,88$).

Thứ tư, định lượng trữ lượng sinh khối và carbon theo trạng thái rừng:

  • Dữ liệu nguyên văn từ luận án: "Sinh khối và dự trữ carbon trên mặt đất đối với $\text{R}{\text{kx}}$ tương ứng là $87,5 \text{ tấn/ha}$ và $41,1 \text{ tấn/ha}$. Sinh khối và dự trữ carbon trên mặt đất đối với $\text{R}{\text{tr}}$ tương ứng là $57,0 \text{ tấn/ha}$ và $26,8 \text{ tấn/ha}$".
  • Sự phân hóa sinh khối và carbon thể hiện rõ nét qua các pha diễn thế phục hồi:
    • Kiểu $\text{R}_{\text{kx}}$: Trạng thái IIIA1 đạt $47,9 \text{ tấn/ha}$ ($22,5 \text{ tấn C/ha}$); Trạng thái IIIA2 đạt $85,3 \text{ tấn/ha}$ ($40,1 \text{ tấn C/ha}$); Trạng thái IIIA3 đạt $147,0 \text{ tấn/ha}$ ($69,1 \text{ tấn C/ha}$).
    • Kiểu $\text{R}_{\text{tr}}$: Trạng thái IIIA1 đạt $35,1 \text{ tấn/ha}$ ($16,5 \text{ tấn C/ha}$); Trạng thái IIIA2 đạt $92,0 \text{ tấn/ha}$ ($43,2 \text{ tấn C/ha}$).

Thứ năm, tổng trữ lượng bể carbon cảnh quan tại Vườn Quốc gia Phước Bình:

  • Trích dẫn trực tiếp từ luận án: "Tổng sinh khối và dự trữ carbon trên mặt đất đối với $\text{R}{\text{kx}}$ thuộc Vườn quốc gia Phước Bình tương ứng là $243,7 \times 10^3 \text{ tấn}$ và $114,5 \times 10^3 \text{ tấn}$, còn $\text{R}{\text{tr}}$ tương ứng là $85,5 \times 10^3 \text{ tấn}$ và $40,3 \times 10^3 \text{ tấn}$".
Kiểu Thảm Thực Vật / Trạng Thái Sinh Khối Trên Mặt Đất (AGB, tấn/ha) Trữ Lượng Carbon (C, tấn/ha) Tổng Sinh Khối VQG Phước Bình ($10^3$ tấn) Tổng Carbon VQG Phước Bình ($10^3$ tấn C) Mô Hình Allometric Tối Ưu Cho Tổng Sinh Khối Cá Thể
Rừng Kín Thường Xanh ($\text{R}_{\text{kx}}$) 87,5 (TB) 41,1 (TB) 243,7 114,5 $\text{B}_{\text{To}} = a \cdot D^b$
Rừng Thưa Nửa Thường Xanh ($\text{R}_{\text{tr}}$) 57,0 (TB) 26,8 (TB) 85,5 40,3 $\text{B}_{\text{To}} = a \cdot D^b$

Implications đa chiều

  • Đóng góp lý luận lâm học: Cung cấp minh chứng sinh thái định lượng về chu trình trao đổi vật chất và năng lượng của thảm thực vật nhiệt đới khô hạn Nam Trung Bộ, khẳng định vai trò của cấu trúc loài ưu thế trong việc duy trì ổn định bể chứa carbon.
  • Đổi mới phương pháp luận: Đơn giản hóa quy trình điều tra kiểm kê carbon bằng cách chỉ ra rằng việc sử dụng biến đơn đường kính $D$ kết hợp với hàm Korsun–Strand mang lại hiệu quả tương đương hoặc vượt trội so với các mô hình phức tạp đòi hỏi đo chiều cao $H$ hoặc bóc tách thể tích $V$.
  • Ứng dụng thực tiễn và chính sách:
    • Cung cấp bộ công cụ toán học và bảng tra hệ số chuyển đổi chuẩn xác để Sở Nông nghiệp & PTNT Ninh Thuận, Vườn Quốc gia Phước Bình và các công ty lâm nghiệp (như Công ty TNHH MTV Lâm nghiệp Tân Tiến) tính toán chính xác đơn giá chi trả dịch vụ môi trường rừng (PFES) theo Nghị định 99/2010/NĐ-CP.
    • Thiết lập mức phát thải tham chiếu rừng (FREL/FRL) cấp tỉnh phục vụ trực tiếp cho các dự án tín chỉ carbon thị trường tự nguyện (VCM) và REDD+.

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 mang tính khách quan:

  1. Giới hạn về bể carbon nghiên cứu: Đề tài tập trung chuyên sâu vào sinh khối trên mặt đất ($\text{AGB}$); chưa tiến hành đo đạc thực nghiệm bể sinh khối dưới mặt đất (hệ thống rễ cây - BGB), vật rơi rụng (litter), gỗ chết mục (deadwood) và bể carbon hữu cơ trong đất ($\text{SOC}$). Việc gán tỷ lệ chuyển đổi cho rễ cây vẫn phải dựa vào các hệ số mặc định của IPCC.
  2. Giới hạn về đại diện loài thu mẫu: Do quy định bảo tồn nghiêm ngặt của Luật Bảo vệ và Phát triển Rừng tại Vườn Quốc gia, mẫu chặt hạ giải tích 88 cây tập trung vào các loài cây ưu thế và đồng ưu thế sinh thái phổ biến; chưa bao phủ được toàn bộ các loài cây gỗ quý hiếm hoặc các loài rải rác thuộc tầng dưới tán.
  3. Giới hạn không gian và khí hậu: Nghiên cứu thực hiện tập trung tại Vườn Quốc gia Phước Bình (Ninh Thuận), mang đặc trưng vi khí hậu nhiệt đới gió mùa khô hạn cục bộ; việc ngoại suy cho các vùng sinh thái có lượng mưa cao hơn đòi hỏi phải tái chuẩn hóa tham số.

Định hướng nghiên cứu tương lai (Future Research Agenda):

  • Hướng 1: Phát triển phương pháp đào mẫu đất và siêu âm rễ cây không phá hủy để xây dựng hệ phương trình allometric chuyên biệt cho sinh khối dưới mặt đất ($\text{BGB}$) và đánh giá động thái bể carbon đất ($\text{SOC}$) theo độ sâu $0–100\text{ cm}$.
  • Hướng 2: Ứng dụng công nghệ viễn thám hiện đại (LiDAR quét tia laser mặt đất và UAV gắn cảm biến đa phổ) tích hợp với các phương trình allometric của luận án nhằm lập bản đồ phân bố trữ lượng carbon rừng theo thời gian thực (real-time carbon mapping).
  • Hướng 3: Thiết lập hệ thống ô định vị sinh thái lâu dài (Permanent Sample Plots - PSP) để theo dõi tốc độ tăng trưởng sinh khối hàng năm ($\Delta \text{AGB}$) và khả năng hấp thụ $\text{CO}_2$ thuần (Net Primary Productivity - NPP) dưới tác động của biến đổi khí hậu.

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

  • Ảnh hưởng học thuật: Công trình tạo tiền lệ nghiên cứu sinh thái định lượng cho kiểu rừng khô hạn nhiệt đới tại Việt Nam; cung cấp nguồn dữ liệu tham chiếu giá trị cao cho các nghiên cứu tổng hợp (meta-analysis) quốc tế về sinh khối rừng nhiệt đới tại khu vực Đông Nam Á.
  • Chuyển đổi quản trị ngành lâm nghiệp: Cung cấp công cụ kỹ thuật số hóa dữ liệu lâm nghiệp, giúp các ban quản lý rừng đặc dụng và phòng điều tra quy hoạch rừng chuyển đổi từ phương pháp ước tính trữ lượng gỗ truyền thống sang kiểm kê tích hợp trữ lượng sinh khối - carbon.
  • Tác động chính sách công: Cung cấp luận cứ khoa học chuẩn xác để Ủy ban Nhân dân tỉnh Ninh Thuận và Bộ Nông nghiệp & PTNT phê duyệt biểu giá dịch vụ môi trường rừng, phân bổ công bằng quỹ bảo vệ và phát triển rừng cho cộng đồng nhận khoán bảo vệ rừng.
  • Lợi ích môi trường và xã hội: Góp phần bảo tồn diện tích 94.635 ha rừng tự nhiên tại Bác Ái và Ninh Sơn; nâng cao sinh kế bền vững cho đồng bào dân tộc thiểu số địa phương thông qua nguồn thu từ chi trả dịch vụ môi trường rừng và tiềm năng thương mại hóa tín chỉ 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 phương pháp luận thu thập mẫu phá hủy chuẩn tắc kết hợp mô hình hóa hồi quy phi tuyến Korsun–Strand; kế thừa các phương trình tương quan để mở rộng nghiên cứu sinh thái học.
  • Các Nhà Khoa học Lâm nghiệp & Sinh thái: Sở hữu bộ thông số thực nghiệm đầy đủ về cấu trúc loài ($\text{IVI}%$) và tương quan hình thái cây rừng khô hạn ven biển Nam Trung Bộ.
  • Ban Quản lý Rừng Đặc dụng & Doanh nghiệp Lâm nghiệp: (Điển hình như VQG Phước Bình, Công ty Lâm nghiệp Tân Tiến) Áp dụng trực tiếp bảng tính sinh khối dựa trên đường kính $D$ và tiết diện ngang $G$ để kiểm kê tài nguyên rừng nhanh chóng, chính xác, tiết kiệm kinh phí.
  • Các Nhà Hoạch định Chính sách & Cơ quan Quản lý: (Sở NN&PTNT Ninh Thuận, Cục Lâm nghiệp) Sử dụng các hằng số tích lũy carbon ($41,1 \text{ tấn C/ha}$ và $26,8 \text{ tấn C/ha}$) để lập báo cáo kiểm kê khí nhà kính cấp tỉnh và xây dựng đề án tín chỉ carbon tham gia thị trường quốc tế.

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 lý thuyết độc đáo nhất là việc chứng minh sự suy giảm phi tuyến của hệ số điều chỉnh sinh khối ($\text{BEF}_i$) theo đường kính cây cá thể ($D$) và sự gia tăng của tỷ lệ cành lá theo hàm logarit ($\text{R}_i = a + b \cdot \ln(D)$). Luận án đã mở rộng Lý thuyết Tương quan Sinh trưởng Allometric (WBE Model) trong sinh thái học thực vật, bác bỏ giả định hệ số $\text{BEF}$ hằng số của các khung hướng dẫn IPCC mặc định đối với rừng khô hạn.

2. Đổi mới phương pháp luận nổi bật nhất so với các nghiên cứu trước đây là gì?
Trả lời: So với các mô hình quốc tế như Chave et al. (2005) vốn phụ thuộc vào biến chiều cao ($H$) và khối lượng thể tích gỗ ($\rho$), luận án đã chứng minh rằng tại hiện trường rừng tự nhiên nhiệt đới gió mùa, mô hình đơn biến đường kính ($D$) kết hợp hàm phân bố Korsun–Strand cho các cấu phần tán lá mang lại độ chính xác cao hơn, đồng thời triệt tiêu hoàn toàn sai số hệ thống do đo đạc chiều cao tán rừng rậm rạp gây ra.

3. Phát hiện bất ngờ nhất và có ý nghĩa thống kê cao nhất từ dữ liệu thực nghiệm là gì?
Trả lời: Phát hiện bất ngờ nhất là dạng hàm tối ưu cho sinh khối tán lá ($\text{B}{\text{L}}$) và cành ($\text{B}{\text{C}}$) ở rừng kín thường xanh không phải là hàm lũy thừa truyền thống mà là hàm Korsun–Strand ($\text{Exp}(a + b/D + c/D^2)$). Điều này phản ánh chính xác quy luật sinh lý học: khi đường kính thân cây tăng lên mức cổ thụ, sinh khối cành lá đạt trạng thái bão hòa do giới hạn cơ học và sinh lý vận chuyển nước của mạch gỗ.

4. Luận án có cung cấp giao thức nhân rộng (Replication Protocol) rõ ràng không?
Trả lời: Có. Luận án cung cấp hệ thống giao thức chi tiết từng bước: từ kỹ thuật bố trí ô tiêu chuẩn 0,2–1,0 ha, kỹ thuật chặt hạ phân đoạn thân–cành–lá, phương pháp cắt thớt giải tích 3–5 cm, quy trình sấy mẫu chuẩn $105^\circ\text{C}$, đến các thuật toán hồi quy phi tuyến tính và bộ tham số thực nghiệm cụ thể cho từng phương trình.

5. Chương trình nghị sự nghiên cứu 10 năm tới được phác thảo như thế nào?
Trả lời: Luận án định hình lộ trình nghiên cứu 10 năm tập trung vào: (1) Số hóa và tích hợp phương trình allometric lên cảm biến viễn thám LiDAR/UAV; (2) Đo đạc thực nghiệm bể carbon dưới mặt đất ($\text{BGB}$) và carbon hữu cơ trong đất ($\text{SOC}$); (3) Xây dựng mô hình cân bằng carbon cảnh quan dưới các kịch bản biến đổi khí hậu và can thiệp lâm sinh tỉa thưa.

Kết luận

  1. Xây dựng thành công hệ thống phương trình allometric phi tuyến đơn biến và đa biến tối ưu để ước lượng sinh khối trên mặt đất cho từng bộ phận (thân, cành, lá) và toàn cây đối với hai kiểu rừng tự nhiên đặc trưng ($\text{R}{\text{kx}}$ và $\text{R}{\text{tr}}$) tại Ninh Thuận.
  2. Phát hiện và mô hình hóa thành công quy luật suy giảm phi tuyến của hệ số điều chỉnh sinh khối ($\text{BEF}_i$) theo đường kính cây cá thể theo dạng hàm $\text{BEF}_i = (a + b/D)^2$ và $\text{BEF}_i = a \cdot D^{-b}$.
  3. Xác lập mô hình tương quan Korsun–Strand cấp quần thụ lâm phần cho phép dự báo trực tiếp sinh khối và trữ lượng carbon từ các chỉ số điều tra truyền thống (Tiết diện ngang $G$ và Trữ lượng $M$).
  4. Định lượng chính xác trữ lượng sinh khối và carbon trung bình: Rừng kín thường xanh hơi khô đạt $87,5 \text{ tấn sinh khối/ha}$ ($41,1 \text{ tấn C/ha}$); Rừng thưa nửa thường xanh hơi khô đạt $57,0 \text{ tấn sinh khối/ha}$ ($26,8 \text{ tấn C/ha}$).
  5. Tổng kết trữ lượng carbon lưu giữ trên mặt đất tại Vườn Quốc gia Phước Bình đạt $114,5 \times 10^3 \text{ tấn C}$ ($\text{R}{\text{kx}}$) và $40,3 \times 10^3 \text{ tấn C}$ ($\text{R}{\text{tr}}$), xác lập bộ thông số Tier 3 vững chắc cho các cơ chế chi trả dịch vụ môi trường rừng (PFES) và kiểm kê khí nhà kính quốc gia.
  6. Mở ra 3 hướng nghiên cứu sinh thái - viễn thám mới, đóng góp di sản học thuật thực nghiệm có giá trị lâu dài cho ngành lâm học và quản lý tài nguyên thiên nhiên Việt Nam.