Tổng quan về luận án

Cháy rừng tại các hệ sinh thái rừng lá kim á nhiệt đới núi cao là một trong những thách thức sinh thái và quản trị tài nguyên nghiêm trọng nhất dưới tác động của biến đổi khí hậu toàn cầu. Nghiên cứu của NCS. Lê Văn Hương (2020) với đề tài "Nghiên cứu đề xuất giải pháp phòng cháy cho rừng Thông ba lá (Pinus kesiya) tại Vườn quốc gia Bidoup - Núi Bà, tỉnh Lâm Đồng" (chuyên ngành Quản lý tài nguyên rừng, mã số 96 20 211, Trường Đại học Lâm nghiệp) đại diện cho một công trình khoa học tiên phong, đặt nền móng định lượng cho công tác phòng cháy, chữa cháy rừng (PCCCR) tại vùng lõi Khu dự trữ sinh quyển thế giới Langbiang và Vườn di sản ASEAN.

                  +----------------------------------------------+
                  |  HỆ THỐNG DỰ BÁO & PHÒNG CHÁY TỔNG HỢP       |
                  |     RỪNG THÔNG BA LÁ (LÊ VĂN HƯƠNG, 2020)    |
                  +----------------------------------------------+
                                          |
          +-------------------------------+-------------------------------+
          |                                                               |
+---------------------+                                         +-------------------+
| TIỂU MÔ HÌNH NHIÊN LIỆU                                       | TIỂU MÔ HÌNH KHÍ TƯỢNG
| - Vật liệu khô: m1                                            | - Biến cơ bản: P, T, H
| - Vật liệu tươi: m2                                           | - Chỉ số hạn: LANG, DEMA,
| - Hệ số bắt cháy: K = m1/(a+b*m1)                             |   SELY, IVA, THORW
| - Phần trăm cháy: Pc                                          | - Chỉ số cháy: ANGS, SHAR,
| - Thời gian cháy: Tc                                          |   SUL, VIN, MNI
+---------------------+                                         +-------------------+
          |                                                               |
          +-------------------------------+-------------------------------+
                                          |
                                          v
                  +----------------------------------------------+
                  | TÍCH HỢP ĐA BIẾN (PCA, FA, MDSA, CCA)         |
                  | & PHÂN CẤP NGUY CƠ (MAHALANOBIS D2, FCF)      |
                  +----------------------------------------------+
                                          |
                                          v
                  +----------------------------------------------+
                  | GIẢI PHÁP CAN THIỆP SINH THÁI THỰC ĐỊA       |
                  | - Đốt chỉ định (Prescribed Burning)          |
                  | - Xử lý thảm mục cục bộ (Tải lượng >20 t/ha) |
                  | - Bản đồ phân vùng trường nhiệt viễn thám    |
                  +----------------------------------------------+

Về mặt bối cảnh khoa học, Vườn quốc gia Bidoup - Núi Bà (VQG BDNB) sở hữu diện tích tự nhiên 70.038 ha, trong đó diện tích rừng có nguy cơ cháy cao lên tới khoảng 30.000 ha (chiếm hơn 42% diện tích toàn vườn). Thống kê thực tế giai đoạn 2005–2017 ghi nhận 84 vụ cháy rừng gây thiệt hại 229,4 ha, tập trung chủ yếu vào trạng thái rừng Thông ba lá (Pinus kesiya Royle ex Gordon). Khoảng trống nghiên cứu (research gap) cốt lõi được nhận diện là sự đứt gãy giữa các mô hình dự báo khí tượng vĩ mô truyền thống (chỉ số Nesterov, Angstrom, Lang) với các đặc tính lý - hóa và động thái tích lũy vi mô của vật liệu cháy (VLC) dưới tán rừng thông tại các đai cao nhiệt đới gió mùa. Phần lớn các phương pháp trước đây chỉ dựa vào biến thời tiết đơn lẻ hoặc phân cấp định tính, thiếu vắng mô hình hóa toán học chuẩn xác về mối tương quan giữa tỷ lệ VLC khô ($m_1$), VLC tươi ($m_2$), tổng sinh khối VLC ($M$), phần trăm cháy hết ($P_c$) và thời gian cháy ($T_c$).

Luận án thiết lập 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: Cấu trúc lâm phần và đặc tính sinh thái của các trạng thái rừng Thông ba lá tại các đai cao VQG BDNB ảnh hưởng như thế nào đến sự hình thành tải lượng VLC?
    • H1: Sự tích lũy VLC và nguy cơ cháy biến thiên có quy luật theo cấp tuổi rừng trồng (cấp I đến cấp IV) và trạng thái rừng tự nhiên.
  • RQ2: Mối quan hệ tương quan toán học giữa khối lượng VLC khô ($m_1$), VLC tươi ($m_2$) và tổng sinh khối ($M$) quyết định khả năng bắt cháy như thế nào?
    • H2: Tồn tại một mô hình phi tuyến xác định hệ số bắt cháy $K = m_1 / M$ biểu diễn mức độ nhạy cảm nhiệt của thảm thực bì.
  • RQ3: Các tổ hợp chỉ số khô hạn và chỉ số nguy cơ cháy rừng đa biến có khả năng tối ưu hóa độ chính xác dự báo mùa cháy và cấp độ nguy hiểm của lửa rừng so với mô hình đơn biến hay không?
    • H3: Phân tích đa biến (PCA, FA, MDSA, FCF, khoảng cách Mahalanobis) cho phép phân tách chuẩn xác các cấp nguy cơ cháy từ dữ liệu chuỗi thời gian khí tượng.
  • RQ4: Giải pháp can thiệp kỹ thuật nào (đốt chỉ định, xử lý cơ giới, băng cản lửa) khả thi và an toàn sinh thái cho rừng Thông ba lá tại VQG BDNB?
    • H4: Kỹ thuật đốt chỉ định (prescribed burning) dựa trên định lượng ngưỡng VLC và hệ số $K$ có thể làm giảm triệt để nguy cơ cháy lớn mà không làm suy thoái thảm thực vật bản địa.

Khung lý thuyết của nghiên cứu tích hợp ba trụ cột: Động học nhiệt hóa học của phản ứng cháy (Combustion Kinetics của Belov, 1982; Kaplukov, 1982), Lý thuyết chu kỳ Cỏ - Lửa (Grass-Fire Cycle của D’Antonio & Vitousek, 1992) và Khung quản lý lửa rừng tích hợp / Nghịch lý của Lửa (Fire Paradox của Kraus, 2010; FAO, 2003). Đóng góp đột phá của luận án là việc xây dựng thành công phương trình hồi quy phi tuyến xác định hệ số bắt cháy $K = m_1 / (2,902691 + 0,46979 \cdot m_1)$ với hệ số tương quan $r = 0,94778$ ($p = 1,99 \times 10^{-15} \ll 0,05$), đồng thời chuẩn hóa ma trận phân loại nguy cơ cháy bằng hàm biệt định Fisher (FCF) trên tập dữ liệu thực địa 340 ô tiêu chuẩn trải rộng trên 4 khu vực đại diện (Đưng K’nớ, Đưng Iar Jiêng, Cổng Trời, Bidoup).

Literature Review và Positioning

Tổng quan y văn quốc tế và trong nước về sinh thái học lửa rừng và quản lý lâm phần Pinus kesiya phản ánh một tiến trình phát triển học thuật sâu rộng qua nhiều thập kỷ:

+-------------------------------------------------------------------------------+
|                             TIẾN TRÌNH Y VĂN LỬA RỪNG                         |
+-------------------------------------------------------------------------------+
|  1920 - 1950: Khởi xướng chỉ số khí hậu đơn biến                              |
|  - Lang (1920): Chỉ số P/T                                                    |
|  - De Martonne (1926): P/(T+10)                                               |
|  - Angstrom (1942): I = (R/20) + (27-T)/10                                    |
|  - Nesterov (1949): Chỉ số khô hạn thực nghiệm NI = Sum[T(T-D)]               |
+-------------------------------------------------------------------------------+
                                      |
                                      v
+-------------------------------------------------------------------------------+
|  1960 - 1990: Mô hình động học lan truyền & Động thái vật liệu rơi rụng       |
|  - Schnell (1962), Schmid (1964): Nguồn gốc phát sinh Thông ba lá & lửa rừng  |
|  - McArthur (1967): McArthur Forest Fire Danger Index (FFDI)                  |
|  - Rothermel (1972): Mô hình toán học lan truyền mặt cầu lửa                  |
|  - Belov (1982), Kaplukov (1982): Nhiệt hóa học VLC thảm mục (C:48.5%, H:6.1%)|
|  - Das & Ramakrishnan (1985): Rơi rụng Pinus kesiya 6.090-8.984 kg/ha/năm     |
+-------------------------------------------------------------------------------+
                                      |
                                      v
+-------------------------------------------------------------------------------+
|  1990 - 2015: Sinh thái học tương hỗ & Tích hợp viễn thám đa biến            |
|  - D'Antonio & Vitousek (1992): Thuyết Chu kỳ Cỏ - Lửa (Grass-Fire Cycle)    |
|  - Cheney (1996), Gould (2007), Anderson (2008): Cấu trúc tầng nhiên liệu     |
|  - Vương Văn Quỳnh (2005), Bế Minh Châu (2012): Độ ẩm VLC & Chỉ số P7i         |
|  - Kraus (2010), FAO (2003): Quản lý lửa tích hợp & Nghịch lý của Lửa         |
+-------------------------------------------------------------------------------+
                                      |
                                      v
+-------------------------------------------------------------------------------+
|  LUẬN ÁN LÊ VĂN HƯƠNG (2020): Định lượng vi mô K, m1, m2 tích hợp MDSA,       |
|  PCA, FCF & Quy chuẩn hóa Đốt chỉ định cho Pinus kesiya nhiệt đới núi cao    |
+-------------------------------------------------------------------------------+

Nghiên cứu vị trí của loài Pinus kesiya bắt đầu từ các khảo cứu thực vật học của Gaussen (1954), Styles & Burley (1972), Burley (1980, 1990), chỉ rõ biên độ phân bố từ vĩ độ 28°N (biên giới Ấn - Trung) trải dài xuống cực Nam tại vĩ độ 12°N thuộc cao nguyên Di Linh và Langbiang (Việt Nam). Về nguồn gốc sinh thái, tồn tại hai luồng quan điểm đối lập: giả thuyết của B. Schnell (1962) khẳng định cao nguyên Langbiang là cái nôi phát sinh nguyên thủy nhờ các vách đá cheo leo và tác động chọn lọc lặp đi lặp lại của lửa nhân sinh; ngược lại, M. Schmid (1964) đề xuất giả thuyết di lưu quần xã từ dãy Himalaya qua ngả Miến Điện, Vân Nam, men theo dãy Trường Sơn để định cư nhờ đặc tính chịu lửa vượt trội so với cây lá rộng.

Trong lĩnh vực động lực học VLC và đánh giá tính dễ cháy, các công trình của Belov (1982) và Kaplukov (1982) đã đặt nền móng nhiệt động lực học bằng cách chứng minh 1 kg VLC khô kiệt của rừng thông chứa 48,5% Carbon, 6,1% Hydro, 38,7% Oxy, sinh nhiệt lượng cháy hoàn toàn $Q_c = 18.800\text{ kJ/kg}$ theo công thức: $$Q = 3,8 \times 10^3 \cdot K \cdot C$$ (với $K$ là hệ số oxy và $C$ là hàm lượng carbon). D’Antonio & Vitousek (1992) nâng tầm thành lý thuyết hệ sinh thái với mô hình "Chu kỳ Cỏ - Lửa" (Grass-Fire Cycle), chỉ ra rằng sự xâm lấn của các loài cỏ ưa sáng làm gia tăng VLC khô thoáng khí, thúc đẩy lửa bùng phát thường xuyên, làm suy giảm cây gỗ và gây thất thoát từ 8 đến 50 triệu tấn Carbon vào khí quyển như trường hợp cỏ Bromus tectorum ở vùng Great Basin (Hoa Kỳ).

Đối với công tác dự báo khí tượng, y văn thế giới phát triển từ các mô hình bán thực nghiệm đơn biến như chỉ số Lang ($P/T$, 1920), De Martonne ($P/(T+10)$, 1926), Selyaninov ($P / (0,1 \sum T)$, 1928), Thornthwaite ($P/\text{PET}$, 1948), Angstrom ($I = (R/20) + (27-T)/10$, 1942), Nesterov ($NI = \sum T(T-D)$, 1949), cho đến các hệ thống phức hợp như McArthur Forest Fire Danger Index (FFDI) (McArthur, 1967; Dowdy et al., 2009), Canadian Fire Weather Index (FWI) (Van Wagner, 1987), và chỉ số Sharples ($F = (1/H) \cdot f(U)$, 2009). Tại Việt Nam, các công trình kế thừa tiêu biểu gồm chỉ số ngày khô hạn liên tục $H$ và chỉ số Nesterov hiệu chỉnh của Phạm Ngọc Hưng (1988, 1994), phương pháp dự báo theo độ ẩm tuyệt đối VLC của Bế Minh Châu (1999, 2002, 2012), chỉ số tổng hợp $P_{7i}$ của Vương Văn Quỳnh (2005, 2012), và mô hình hồi quy đa biến địa phương của Trần Xuân Hiền (2007).

So sánh quốc tế:

  1. So với nghiên cứu của Ashesh Kumar Das & P. Ramakrishnan (1985) tại Mawlai, bang Meghalaya, Đông Bắc Ấn Độ (độ cao 1.250 m), lượng rơi rụng hàng năm của rừng Thông ba lá 7, 15 và 22 tuổi đạt từ 6.090 đến 8.984 kg/ha/năm, trong đó lá kim chiếm 78–98%. Luận án của Lê Văn Hương xác nhận tính quy luật này tại VQG BDNB nhưng phát hiện thêm rằng tại các quần thụ già lâu năm không cháy, tổng tích lũy VLC tích tụ vượt ngưỡng 20 tấn/ha (20.000 kg/ha), tạo thành thảm nhiên liệu cực kỳ nguy hiểm.
  2. So với nghiên cứu đốt trước phòng cháy của Mc Callum (2006) tại căn cứ quân sự Avon Park Air Force Range (bang Florida, Hoa Kỳ) và các công trình của Cheney (1996), Gould et al. (2007) trên rừng bạch đàn và thông tại Australia, luận án đã mở rộng mô hình hóa mối quan hệ giữa thời gian cháy hết ($T_c$) và tỷ lệ cháy sạch ($P_c$) trong điều kiện địa hình chia cắt phức tạp của vùng núi cao nhiệt đới Tây Nguyên.

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

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

Luận án tạo ra bước tiến quan trọng trong sinh thái học lửa rừng bằng việc bổ sung và hoàn thiện ba luận điểm lý thuyết cốt lõi:

+---------------------------------------------------------------------------------+
|                       ĐÓNG GÓP LÝ THUYẾT & MÔ HÌNH HÓA                          |
+---------------------------------------------------------------------------------+
|  Mô hình Động học Khả năng Bắt cháy (Ignitability Kinetics Model):             |
|                                                                                 |
|                   m1                                                            |
|       K = ---------------------      (r = 0,94778; p = 1,99 x 10^-15)           |
|           2,902691 + 0,46979*m1                                                 |
|                                                                                 |
|  - Khắc phục sự phụ thuộc đơn thuần vào biến khí tượng vĩ mô.                   |
|  - Định lượng hóa ranh giới chuyển pha bắt lửa mặt đất sang cháy tán.           |
+---------------------------------------------------------------------------------+
                                      |
                                      v
+---------------------------------------------------------------------------------+
|  Mệnh đề Lý thuyết (Theoretical Propositions):                                  |
|  - P1: Trạng thái bắt cháy bão hòa tiệm cận phi tuyến theo khối lượng khô m1.   |
|  - P2: Tỷ lệ tươi m2 đóng vai trò bức cản nhiệt làm chậm quá trình nhiệt phân.  |
|  - P3: Tỷ lệ Pc và thời gian Tc phụ thuộc đồng thời vào ma trận không gian VLC. |
|  - P4: Chuyển dịch mô hình quản trị: "Triệt tiêu lửa tuyệt đối" sang           |
|       "Chủ động điều tiết tải lượng bằng Đốt chỉ định" (Kraus, 2010).          |
+---------------------------------------------------------------------------------+
  • Mở rộng lý thuyết Động học Cháy (Combustion Kinetics): Luận án chứng minh rằng tính bắt cháy của thảm thực bì không chỉ phụ thuộc vào độ ẩm tuyệt đối tức thời, mà được kiểm soát trực tiếp bởi tỷ lệ khối lượng vật liệu khô $m_1$ trên tổng khối lượng $M$ thông qua hệ số $K = m_1 / M$. Tác giả đã lượng hóa mối quan hệ phi tuyến giữa $K$ và $m_1$ bằng phương trình: $$K = \frac{m_1}{2,902691 + 0,46979 \cdot m_1}$$ với $r = 0,94778$ ($p = 1,99 \times 10^{-15} \ll 0,05$), chứng minh bản chất động học bão hòa tiệm cận của quá trình sấy khô tự nhiên trong thảm mục rừng thông.
  • Hệ thống hóa các mệnh đề lý thuyết (Theoretical Propositions):
    • Mệnh đề P1: Tính dễ cháy của quần xã Thông ba lá tỉ lệ thuận phi tuyến với khối lượng VLC khô ($m_1$) và bị ức chế bởi khối lượng VLC tươi ($m_2$).
    • Mệnh đề P2: Tỷ lệ phần trăm cháy hết ($P_c$) đạt ngưỡng bùng phát khi hệ số $K \ge 0,5$, tại đó năng lượng giải phóng từ $m_1$ vượt qua nhiệt dung riêng cần thiết để bốc thoát hơi nước của $m_2$.
    • Mệnh đề P3: Thời gian cháy hết ($T_c$) là hàm đa biến đồng quy phụ thuộc vào tổng sinh khối $M$, độ lèn chặt và tỷ lệ $m_1/m_2$.
    • Mệnh đề P4: Quản trị lửa rừng nhiệt đới đòi hỏi sự chuyển đổi mô hình (paradigm shift) từ "dập tắt lửa thụ động" sang "quản lý tải lượng nhiên liệu chủ động thông qua lửa chỉ định" phù hợp với luận thuyết Nghịch lý của Lửa (Fire Paradox).

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

Khung phân tích của luận án tích hợp đa tầng dữ liệu:

[DỮ LIỆU ĐỊA HÌNH & LÂM PHẦN] ---> [THÍ NGHIỆM ĐỐT HIỆN TRƯỜNG: m1, m2, M]
                                                   |
                                                   v
[DỮ LIỆU KHÍ TƯỢNG VĨ MÔ] --------> [TÍNH TOÁN CHỈ SỐ: LANG, SELY, MNI, FCF]
                                                   |
                                                   v
                              [PHÂN TÍCH ĐA BIẾN: PCA, FA, MDSA, CCA]
                                                   |
                                                   v
                              [QUY CHUẨN KỸ THUẬT: ĐỐT CHỈ ĐỊNH & DỰ BÁO]
  • Tích hợp liên ngành: Kết hợp Nhiệt động học quá trình cháy, Sinh thái học cảnh quan và Khí tượng sinh vật học (Gaussen-Walter, Thornthwaite).
  • Tiếp cận phân tích thứ bậc: Phân tích từ cấp độ vi mô (thành phần hóa học, độ ẩm bắt cháy của lá kim, cành nhánh $<6\text{ mm}$) đến cấp độ trung mô (cấu trúc lâm phần cấp tuổi I–IV, rừng tự nhiên tại 4 vùng sinh thái) và cấp độ vĩ mô (trường nhiệt, chỉ số khô hạn toàn cảnh quan VQG BDNB).
  • Điều kiện biên (Boundary Conditions): Áp dụng chuẩn xác cho các hệ sinh thái rừng Thông ba lá phân bố ở đai cao từ 1.300 m đến trên 2.000 m thuộc vùng Tây Nguyên và các sinh cảnh tương đồng ở Đông Nam Á, có chế độ khí hậu mùa khô rõ rệt kéo dài từ 4 đến 6 tháng.

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 (Positivism) kết hợp với chủ nghĩa hiện thực phản biện (Critical Realism) trong quản lý sinh thái. Thiết kế nghiên cứu kết hợp chặt chẽ giữa quan trắc thực nghiệm diện rộng, phân tích nhiệt hóa phòng thí nghiệm và mô hình hóa thống kê đa biến:

+-------------------------------------------------------------------------------+
|                       THIẾT KẾ & QUY MÔ NGHIÊN CỨU                            |
+-------------------------------------------------------------------------------+
|  Tổng số ô tiêu chuẩn (OTC): N = 340 ô nghiên cứu                             |
|  - Bố trí tại 4 phân vùng sinh thái độ cao:                                   |
|    1. Đưng K'nớ          2. Đưng Iar Jiêng                                    |
|    3. Cổng Trời          4. Bidoup                                            |
|  - Trạng thái lâm phần:                                                       |
|    + Rừng trồng Thông ba lá: Cấp tuổi I, II, III, IV                           |
|    + Rừng Thông ba lá tự nhiên (trên đất mùn xám, vỏ nứt dày >3cm)            |
+-------------------------------------------------------------------------------+

Quy trình nghiên cứu rigorous

  • Quy trình thu thập mẫu VLC: Bố trí các ô nghiên cứu dạng vòng tròn và ô bàn cờ đại diện. Tiến hành phân tách tại hiện trường thành phần VLC khô ($m_1$: cành khô, lá rụng, vỏ rụng, thảm mục) và VLC tươi ($m_2$: cỏ tranh, cỏ lúa, dương xỉ, cây bụi) theo tiêu chuẩn đường kính $<6\text{ mm}$ (VLC tinh) và $\ge 6\text{ mm}$ (VLC thô).
  • Tam giác hóa phương pháp (Methodological Triangulation):
    1. Thực nghiệm hiện trường: Cân khối lượng tươi, sấy khô kiệt ở 105°C để xác định độ ẩm tuyệt đối và sinh khối khô.
    2. Thực nghiệm buồng đốt: Đốt mẫu đo thời gian cháy ($T_c$), phần trăm cháy hết ($P_c$), tốc độ lan truyền theo các cấp độ ẩm.
    3. Khí tượng trạm & Viễn thám: Khai thác chuỗi số liệu khí hậu dài hạn (nhiệt độ $T$, độ ẩm tương đối $H$, lượng mưa $P$, tốc độ gió $v$, số giờ nắng $S$) tích hợp dữ liệu kênh nhiệt hồng ngoại vệ tinh Landsat TM/ETM+, MODIS (kênh 3.27 µm) và ASTER.
  • Độ tin cậy và giá trị đo lường: Các mô hình hồi quy và kiểm định thống kê đều đạt mức ý nghĩa cao ($p < 0,001$), sai số hệ thống được triệt tiêu thông qua chuẩn hóa thang đo đa biến.

Data và phân tích

Phân tích số liệu sử dụng hai gói phần mềm chuyên dụng SPSS và Statgraphics với các kỹ thuật cao cấp:

+-------------------------------------------------------------------------------+
|                   CÁC KỸ THUẬT PHÂN TÍCH THỐNG KÊ ĐA BIẾN                     |
+-------------------------------------------------------------------------------+
|  1. Phân tích đồ hình đa chiều (MDSA): Định vị cấu trúc không gian biến.      |
|  2. Phân tích thành phần chính (PCA) & Phân tích nhân tố (FA):                |
|     - Giảm chiều 5 biến khô hạn: LANG, DEMA, SELY, IVA, THORW.                |
|     - Giảm chiều 4 biến chỉ số cháy: ANGS, SHAR, SUL, VIN.                    |
|  3. Phân tích tương quan chuẩn tắc (CCA) & Hàm biệt định (CDF):               |
|     - Liên kết tập biến khí hậu độc lập với tập biến hành vi cháy.            |
|  4. Hàm phân loại Fisher (FCF) & Khoảng cách Mahalanobis (D2):                |
|     - Xây dựng ma trận phân định tự động 5 cấp nguy cơ cháy rừng.             |
+-------------------------------------------------------------------------------+

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

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

  1. Quy luật tích lũy sinh khối VLC theo cấu trúc lâm phần: Rừng trồng Thông ba lá cấp tuổi I và II có tỷ lệ cỏ tươi ($m_2$) lớn nhưng sinh khối tích lũy thấp; ngược lại, rừng cấp tuổi III, IV và rừng tự nhiên có lớp thảm mục và lá rụng khô ($m_1$) đan dày từ 2 đến 5 cm, tổng tải lượng nhiên liệu thường xuyên vượt mức 20 tấn/ha. Thảm cỏ dày rậm kết hợp với vật rụng lá kim (chiếm 78–98% vật rụng hàng năm) tạo thành khối nhiên liệu cực kỳ nguy hiểm trong mùa khô.
  2. Xác lập mô hình động học hệ số bắt cháy $K$: Mối quan hệ giữa hệ số $K$ và khối lượng $m_1$ là một hàm phi tuyến có hệ số tương quan $r = 0,94778$. Khi $m_1$ tăng, $K$ tăng tiệm cận về 1. Kết quả thực nghiệm phân cấp:
    • $K < 0,2$: Không có khả năng cháy.
    • $0,2 \le K < 0,3$: Ít có khả năng cháy.
    • $0,3 \le K < 0,5$: Có khả năng cháy.
    • $0,5 \le K \le 0,7$: Khả năng cháy cao.
    • $K > 0,7$: Khả năng cháy rất cao (cháy bùng phát dữ dội).
  3. Mô hình hóa thời gian cháy ($T_c$) và tỷ lệ cháy sạch ($P_c$): Luận án xây dựng các phương trình bề mặt hồi quy biểu diễn $T_c = f(m_1, m_2, M)$ và $P_c = f(m_1, m_2, K)$. Kết quả chứng minh $P_c$ đạt trên 85–95% khi $K > 0,5$ ngay cả khi độ ẩm không khí chưa đạt mức cực hạn, giải thích hiện tượng cháy lan nhanh bất thường trên các sườn dốc khu vực Cổng Trời và Bidoup.
  4. Phát hiện chiều không gian tối ưu trong dự báo khí tượng: Phân tích PCA và FA trên 5 chỉ số khô hạn (Lang, De Martonne, Selyaninov, Ivanov, Thornthwaite) và 4 chỉ số cháy (Angstrom, Sharples, Sullivan, Viney) đã rút gọn không gian biến thành 2 trục nhân tố chính giải thích trên 85% phương sai. Dự báo nguy cơ cháy rừng đạt độ chính xác phân loại cao nhất khi tích hợp khoảng cách đa biến Mahalanobis ($D^2$) và các hàm phân loại Fisher (FCF) thay vì dùng chỉ số Nesterov đơn lẻ.
  5. Đặc trưng nguồn nhiệt gây cháy: 100% các vụ cháy rừng tại VQG BDNB trong giai đoạn nghiên cứu đều có nguồn gốc từ con người (đốt dọn nương rẫy, xử lý thực bì sai quy trình, săn bắt, đốt đồng cỏ chăn thả và bất cẩn của khách du lịch). Không ghi nhận cháy do sét đánh hay tự bốc cháy sinh học.
+-------------------------------------------------------------------------------+
|            TIÊU CHÍ ĐÁNH GIÁ KHẢ NĂNG CHÁY CỦA VLC THEO HỆ SỐ K              |
+-------------------+-----------------------------------+-----------------------+
|  Giá trị hệ số K  |  Khả năng cháy thảm thực bì       |  Hành vi ngọn lửa      |
+-------------------+-----------------------------------+-----------------------+
|  K < 0,20         |  Không cháy                       |  Tắt ngấm tự nhiên    |
|  0,20 <= K < 0,30 |  Ít có khả năng cháy              |  Cháy âm ỉ, cục bộ    |
|  0,30 <= K < 0,50 |  Có khả năng cháy                 |  Cháy lan mặt đất     |
|  0,50 <= K <= 0,70|  Khả năng cháy cao                |  Cháy mạnh, nhiệt cao |
|  K > 0,70         |  Khả năng cháy rất cao            |  Bùng phát lên tán    |
+-------------------+-----------------------------------+-----------------------+

Implications đa chiều

  • Về mặt lý luận học thuật: Đóng góp luận cứ khoa học chuẩn xác làm phong phú thêm chuyên ngành Quản lý tài nguyên rừng và Sinh thái học lửa rừng nhiệt đới; cung cấp công cụ toán học định lượng thay thế phương pháp chuyên gia định tính.
  • Về mặt phương pháp luận: Thiết lập quy trình mẫu kết hợp giữa điều tra lâm học truyền thống, thực nghiệm nhiệt động học và phân tích thống kê đa biến (PCA, FCF, Mahalanobis) có thể nhân rộng cho các Vườn quốc gia và Khu bảo tồn khác trên toàn quốc.
  • Về ứng dụng thực tiễn:
    • Đề xuất lịch thời vụ mùa cháy chính xác cho VQG BDNB kéo dài từ tháng 12 đến hết tháng 4 năm sau (cao điểm cực hạn vào tháng 3 và 4 khi thời tiết chịu ảnh hưởng của áp thấp nóng phía Tây).
    • Xác định đối tượng và thời điểm áp dụng giải pháp đốt chỉ định (prescribed burning): thực hiện vào đầu mùa khô (tháng 11 - 12), chọn ngày có gió nhẹ ($v < 5\text{ km/h}$), độ ẩm không khí $H > 60%$, đốt ngược chiều gió hoặc đốt theo băng cục bộ để triệt tiêu lớp VLC tích lũy trước khi bước vào cao điểm khô hạn.
  • Về chính sách và quản trị: Cung cấp cơ sở khoa học để Ban quản lý VQG BDNB và Chi cục Kiểm lâm tỉnh Lâm Đồng tái cấu trúc phân vùng kiểm soát cháy, điều chỉnh biểu phân cấp dự báo cháy rừng phù hợp với đặc thù sinh thái vùng cao nguyên.

Limitations và Future Research

  1. Giới hạn không gian và địa hình: Nghiên cứu tập trung chủ yếu tại 4 phân khu trọng điểm của VQG BDNB; tính chất biến thiên vi khí hậu tại các hẻm vực sâu và sườn đón gió đông chưa được bao phủ hoàn toàn bởi các trạm quan trắc tự động.
  2. Giới hạn về mô phỏng 3D hành vi lửa: Luận án tập trung chủ yếu vào động học cháy bề mặt và lớp thảm mục; sự liên kết động lực học giữa cháy mặt đất và hiện tượng cháy cầu ngọn (crown fire) dưới tác động của các cơn lốc nhiệt cục bộ cần được mô phỏng sâu hơn bằng mô hình thủy động lực học tính toán (CFD).
  3. Ảnh hưởng dài hạn của Đốt chỉ định: Nghiên cứu mới dừng lại ở việc đánh giá hiệu quả giảm tải VLC tức thời; tác động sinh thái dài hạn của kỹ thuật đốt chỉ định đối với chu trình dinh dưỡng đất ($N, P$ bốc hơi), hệ vi sinh vật đất và cấu trúc tái sinh tự nhiên của cây họ Dẻ (Fagaceae), Trâm (Myrtaceae) dưới tán thông cần được theo dõi qua các ô định vị lâu dài (Permanent Sample Plots).

Chương trình nghiên cứu tương lai (4 hướng trọng tâm):

  • Xây dựng mạng lưới cảm biến IoT không dây kết hợp camera AI cảnh báo sớm trường nhiệt tại các tiểu khu xung yếu.
  • Ứng dụng ảnh viễn thám siêu phổ (Hyperspectral) và dữ liệu LiDAR để lập bản đồ 3D độ dày thảm mục và phân bố không gian sinh khối VLC theo thời gian thực.
  • Nghiên cứu sâu về động thái bốc thoát khí nhà kính ($CO_2, CH_4, N_2O$) từ các vụ cháy rừng thông nhằm phục vụ thị trường tín chỉ carbon rừng.
  • Đánh giá khả năng chống chịu và phục hồi sau cháy của các loài thực vật đặc hữu trong Vùng lõi Khu dự trữ sinh quyển Langbiang.

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

  • Tác động học thuật: Định hình chuẩn mực nghiên cứu định lượng về VLC tại Việt Nam, đóng vai trò tài liệu tham khảo thiết yếu cho các chương trình đào tạo sau đại học chuyên ngành Lâm học và Quản lý tài nguyên rừng.
  • Chuyển đổi ngành Lâm nghiệp: Thay đổi căn bản phương thức tác nghiệp PCCCR tại Lâm Đồng từ ứng phó thụ động sang điều tiết tải lượng nhiên liệu chủ động, giúp giảm thiểu chi phí huy động lực lượng dập lửa hàng năm.
  • Bảo tồn giá trị toàn cầu: Bảo vệ an toàn sinh thái cho 70.038 ha vùng lõi Khu dự trữ sinh quyển thế giới Langbiang và Vườn di sản ASEAN, ngăn ngừa nguy cơ suy thoái đa dạng sinh học của các quần xã thực vật hạt trần và động vật quý hiếm.
  • Hiệu quả kinh tế - xã hội: Giảm thiểu thiệt hại tài nguyên rừng gỗ, ổn định nguồn nước đầu nguồn cho hệ thống lưu vực sông Đồng Nai, đồng thời bảo đảm an toàn cho các cộng đồng dân cư sinh sống xen kẽ ven rừng.

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

+-------------------+-----------------------------------------------------------+
|  NHÓM ĐỐI TƯỢNG   |  GIÁ TRỊ VÀ LỢI ÍCH TIẾP NHẬN TỪ CÔNG TRÌNH               |
+-------------------+-----------------------------------+-----------------------+
|  Nghiên cứu sinh  |  Khung phân tích thống kê đa biến (PCA, FA, MDSA, FCF);   |
|  & Học viên SĐH   |  Phương pháp luận thực nghiệm đo đạc VLC chuẩn xác.       |
+-------------------+-----------------------------------------------------------+
|  Các nhà Khoa học |  Cơ sở dữ liệu định lượng về sinh thái lửa Pinus kesiya;  |
|  & Giảng viên     |  Mô hình toán học phi tuyến về hệ số bắt cháy K.          |
+-------------------+-----------------------------------------------------------+
|  Ban quản lý VQG, |  Quy trình kỹ thuật đốt chỉ định an toàn;                 |
|  Khu bảo tồn      |  Bản đồ phân vùng nguy cơ cháy theo cấp tuổi lâm phần.    |
+-------------------+-----------------------------------------------------------+
|  Cơ quan Quản lý  |  Luận cứ xây dựng Tiêu chuẩn Quốc gia (TCVN) về PCCCR;    |
|  (Cục Lâm nghiệp) |  Công cụ tối ưu hóa ngân sách và nhân lực mùa khô hạn.    |
+-------------------+-----------------------------------------------------------+

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 thiết lập và kiểm định thành công mô hình toán học thực nghiệm phi tuyến biểu diễn hệ số bắt cháy $K$ theo khối lượng VLC khô $m_1$: $$K = \frac{m_1}{2,902691 + 0,46979 \cdot m_1} \quad (r = 0,94778; ; p = 1,99 \times 10^{-15})$$ Công trình này mở rộng trực tiếp Lý thuyết Động học Quá trình Cháy (Combustion Kinetics của Kaplukov & Belov, 1982) và Lý thuyết Chu kỳ Cỏ - Lửa (D’Antonio & Vitousek, 1992) bằng cách định lượng hóa chính xác tỷ lệ pha trộn giữa vật liệu chết ($m_1$) và vật liệu sống ($m_2$), khắc phục hoàn toàn hạn chế của các chỉ số khí hậu thuần túy vốn không phản ánh được cấu trúc thực bì tại chỗ.

2. Điểm cải tiến phương pháp luận so với các nghiên cứu trước đây là gì?

So với phương pháp dự báo đơn biến của Phạm Ngọc Hưng (1988), Cooper (1991) hay Bế Minh Châu (1999), luận án đã thực hiện bước nhảy vọt về phương pháp luận:

  • Tích hợp đồng thời 5 chỉ số khô hạn cổ điển và 4 chỉ số nguy cơ cháy quốc tế qua các kỹ thuật giảm chiều dữ liệu hiện đại (PCA, FA, MDSA).
  • Sử dụng hàm phân loại Fisher (FCF) và khoảng cách đa biến Mahalanobis ($D^2$) để xây dựng ma trận phân định nguy cơ cháy tự động, giúp triệt tiêu hiện tượng báo động giả và tăng độ tin cậy dự báo lên mức vượt trội so với việc tra cứu bảng thủ công.

3. Phát hiện bất ngờ nhất từ dữ liệu thực nghiệm là gì?

Phát hiện bất ngờ nhất là tính chất bão hòa nhiệt của thảm mục rừng thông già: khi lớp VLC khô $m_1$ đạt ngưỡng tích lũy trên 12–15 tấn/ha, hệ số bắt cháy $K$ tiệm cận ngưỡng nguy hiểm ($>0,7$) và tốc độ lan truyền ngọn lửa không còn phụ thuộc tuyến tính vào độ ẩm tương đối của không khí ($H$). Tại thời điểm này, chỉ cần một mồi lửa nhỏ từ hoạt động nhân sinh cũng đủ kích hoạt đám cháy bùng phát với nhiệt lượng cực đại thiêu hủy toàn bộ thảm thực bì mà các biện pháp chữa cháy thủ công không thể tiếp cận dập tắt.

4. Luận án có cung cấp quy trình tái lặp (Replication Protocol) rõ ràng không?

Có. Luận án cung cấp chi tiết quy trình 5 bước chuẩn hóa có khả năng tái lập hoàn toàn:

  1. Phương pháp bố trí 340 ô tiêu chuẩn theo đai cao và cấp tuổi.
  2. Quy chuẩn lấy mẫu, phân loại $m_1, m_2$ (ngưỡng đường kính 6 mm) và sấy khô tuyệt đối ở 105°C.
  3. Bộ phương trình tương quan hồi quy tường minh kèm tham số $b_0, b_i$.
  4. Bộ ma trận hệ số cho hàm phân loại Fisher (FCF).
  5. Quy trình kỹ thuật đốt chỉ định chi tiết (thời điểm, vận tốc gió, dải phân cách an toàn).

5. Chương trình nghị sự nghiên cứu 10 năm (2020–2030) được định hình ra sao?

Chương trình nghiên cứu 10 năm tiếp theo tập trung vào ba trụ cột:

  • Giai đoạn 2020–2023: Số hóa toàn diện bản đồ nguy cơ cháy dựa trên GIS - Viễn thám trường nhiệt MODIS/Landsat cho toàn bộ lâm phần Lâm Đồng.
  • Giai đoạn 2024–2027: Hoàn thiện quy chuẩn kỹ thuật quốc gia về Đốt chỉ định trong rừng thông bản địa và thiết lập mạng lưới trạm quan trắc tự động IoT.
  • Giai đoạn 2028–2030: Tích hợp mô hình lan truyền lửa 3D với dữ liệu LiDAR viễn thám và đánh giá cân bằng phát thải carbon rừng trong khuôn khổ cam kết Net-Zero của quốc gia.

Kết luận

Công trình nghiên cứu tiến sĩ của Lê Văn Hương (2020) đã giải quyết trọn vẹn và xuất sắc mục tiêu đề ra, mang lại 5 đóng góp cốt lõi:

  1. Định lượng hóa toàn diện cấu trúc VLC: Xác định chính xác quy luật phân bố sinh khối thảm mục và thảm tươi trên 4 cấp tuổi rừng trồng và rừng tự nhiên tại VQG Bidoup - Núi Bà, chỉ rõ ngưỡng tích lũy nhiên liệu nguy hiểm ($>20\text{ tấn/ha}$).
  2. Phát triển mô hình toán học động học cháy đột phá: Xây dựng thành công phương trình tương quan phi tuyến giữa hệ số bắt cháy $K$ và khối lượng VLC khô $m_1$ ($r = 0,94778; p < 0,001$), tạo cơ sở định lượng vững chắc cho sinh thái học lửa rừng.
  3. Cách tân hệ thống dự báo bằng thống kê đa biến: Ứng dụng thành công PCA, FA, MDSA, FCF và khoảng cách Mahalanobis trên chuỗi dữ liệu khí tượng và các chỉ số khô hạn quốc tế, nâng cao độ chính xác phân cấp nguy cơ cháy.
  4. Chuẩn hóa giải pháp đốt chỉ định (Prescribed Burning): Đưa ra luận cứ khoa học và quy trình kỹ thuật can thiệp thực địa an toàn, tạo bước chuyển dịch mô hình từ chữa cháy thụ động sang quản trị nhiên liệu chủ động.
  5. Bảo tồn bền vững Khu dự trữ sinh quyển thế giới: Cung cấp bộ giải pháp kỹ thuật khả thi bảo vệ an toàn diện tích 70.038 ha của VQG BDNB, mở ra 3 hướng nghiên cứu mới về mô phỏng 3D hành vi lửa, viễn thám siêu phổ sinh khối và cân bằng carbon sau cháy cho khoa học lâm nghiệp Việt Nam và khu vực.