Tổng quan về luận án
Biến đổi khí hậu (BĐKH) toàn cầu với biểu hiện cốt lõi là sự gia tăng nhiệt độ bề mặt Trái Đất và sự xáo trộn chế độ thủy văn cực đoan đang trở thành động lực chính thúc đẩy tần suất cũng như mức độ khốc liệt của các vụ cháy rừng trên toàn thế giới. Việt Nam được Ủy ban Liên chính phủ về Biến đổi Khí hậu (IPCC) xếp vào nhóm 5 quốc gia chịu tổn thương nặng nề nhất bởi BĐKH. Trong cấu trúc tài nguyên quốc gia tính đến năm 2013, tổng diện tích rừng đạt khoảng 13,86 triệu ha, trong đó có hơn 6 triệu ha thuộc các hệ sinh thái rừng đặc biệt dễ bắt lửa như rừng tràm (Melaleuca cajuputi), rừng khộp (rừng thưa cây họ Dầu - Dipterocarpaceae), rừng thông (Pinus spp.), rừng bạch đàn (Eucalyptus spp.) và rừng tre nứa. Số liệu thực chứng từ Cục Kiểm lâm giai đoạn 2002 – 2011 ghi nhận 9.689 vụ cháy rừng trên quy mô cả nước, thiêu rụi 55.505 ha rừng tự nhiên và rừng trồng, gây thiệt hại nghiêm trọng về kinh tế sinh thái và phát thải lượng lớn khí nhà kính (CO₂, CO, NOₓ, SO₂) vào khí quyển.
Khoảng trống nghiên cứu (research gap) trọng tâm nằm ở sự thiếu vắng một công cụ định lượng tích hợp có khả năng mô hình hóa nguy cơ cháy rừng theo các kịch bản biến đổi khí hậu vĩ mô tại Việt Nam. Các phương pháp dự báo cháy rừng truyền thống (tiêu biểu là chỉ số Nesterov nguyên bản) bộc lộ sự thiếu tương thích khi áp dụng vào các vùng khí hậu phân mùa phức tạp, đồng thời chỉ dừng lại ở việc dự báo thời tiết tức thời (vài ngày) mà không đánh giá được xu thế dịch chuyển không gian - thời gian của nguy cơ cháy theo các kịch bản phát thải dài hạn.
Luận án của nghiên cứu sinh Lê Sỹ Doanh với đề tài "Nghiên cứu tác động của biến đổi khí hậu đến nguy cơ cháy rừng ở Việt Nam và đề xuất giải pháp ứng phó" (Chuyên ngành Lâm sinh, mã số 62.05, Đại học Lâm nghiệp, 2014) dưới sự hướng dẫn của GS. Vương Văn Quỳnh và PGS. Trần Quang Bảo đã giải quyết toàn diện khoảng trống này. Luận án đặt ra ba câu hỏi nghiên cứu cốt lõi:
- RQ1: Làm thế nào để thiết lập một chỉ số khí hậu đặc thù ($Q_i$) phản ánh chính xác nguy cơ tích lũy nhiệt ẩm của vật liệu cháy dưới tán rừng trong điều kiện khí hậu nhiệt đới gió mùa?
- RQ2: Quy luật biến đổi không gian và thời gian của số ngày có nguy cơ cháy rừng cao và rất cao ($S_{nc45}$) diễn ra như thế nào trên 8 vùng sinh thái của Việt Nam qua các mốc thời gian 2020, 2050 và 2090?
- RQ3: Những giải pháp lâm sinh và quản trị phòng cháy chữa cháy rừng (PCCCR) nào cần được thiết kế thích ứng có phân hóa theo từng vùng sinh thái trọng điểm?
Nghiên cứu vận hành trên ba giả thuyết khoa học:
- H1: Mức độ khô kiệt của vật liệu cháy không chỉ phụ thuộc vào thời tiết ngày hiện tại mà chịu tác động tích lũy phi tuyến tính của nhiệt độ và độ ẩm không khí các ngày liền trước.
- H2: Dưới kịch bản biến đổi khí hậu trung bình B2, số ngày có nguy cơ cháy cao ($S_{nc45}$) gia tăng phi tuyến tính và phân hóa sâu sắc giữa các vùng sinh thái, trong đó Tây Nguyên và Tây Nam Bộ là hai tâm điểm chịu tác động nặng nề nhất.
- H3: Nguy cơ cháy rừng tích hợp là hàm đa biến phụ thuộc đồng thời vào chỉ số khí hậu $Q_i$, cấu trúc thảm thực vật - vật liệu cháy ($E_{ct}, f_{ij}$) và can thiệp nhân sinh.
Khung lý thuyết của công trình kết hợp Lý thuyết Tam giác lửa (Fire Triangle Theory), Cơ chế cân bằng ẩm vật liệu cháy (Fuel Moisture Equilibrium) và Lý thuyết Sinh thái cảnh quan (Landscape Ecology). Về quy mô và tầm ảnh hưởng, nghiên cứu bao phủ toàn diện 8 vùng sinh thái lâm nghiệp quốc gia, khai thác dữ liệu từ 95 trạm khí tượng bề mặt và thiết lập hệ thống ô tiêu chuẩn thực địa đo lường khối lượng thảm khô ($M_{tk}$), thảm tươi ($M_{tt}$) và độ ẩm lúc 13 giờ ($W_{13}$).
Literature Review và Positioning
Nghiên cứu dự báo nguy cơ cháy rừng trên thế giới trải qua ba dòng nghiên cứu chính. Dòng thứ nhất tập trung vào các mô hình khí tượng thực nghiệm cổ điển. Nesterov (1949) phát triển chỉ số khí tượng tổng hợp $P = \sum (T \cdot d)$ (với $T$ là nhiệt độ lúc 13 giờ, $d$ là độ hụt bão hòa không khí). McArthur (1986) tại Úc thiết lập Chỉ số Nguy cơ Cháy rừng (FFDI) tích hợp nhiệt độ, độ ẩm tương đối, tốc độ gió và nhân tố khô hạn đất ($D$). Tại Trung Quốc, Yangmei xây dựng chỉ số bén lửa dựa trên phương trình hồi quy đa biến của bốc hơi, số giờ nắng và nhiệt ẩm lúc 14 giờ. Dòng nghiên cứu thứ hai ứng dụng công nghệ không gian địa lý và viễn thám: Hệ thống Thông tin Cháy rừng Châu Âu (EFFIS) cung cấp bản đồ dự báo nguy cơ cháy cho hơn 30 quốc gia; Hệ thống WFAS (Wildland Fire Assessment System) của Hoa Kỳ mô phỏng độ ẩm vật liệu cháy 7 ngày liên tiếp từ ảnh vệ tinh; các công trình của Chowdhury và Hassan (2011, 2013) tại Alberta (Canada) ứng dụng ảnh MODIS kết hợp chỉ số hạn hán NMDI và chỉ số ẩm TVWI để dự báo chính xác 91,63% đến 98,19% các điểm cháy rừng phương Bắc. Dòng nghiên cứu thứ ba phân tích tác động của BĐKH đến cháy rừng: Nicholls và cộng sự (2005) dự báo số ngày nguy cơ cháy cao ở Đông Nam Úc tăng 15–70% vào năm 2050; Flannigan và cộng sự (2011) mô hình hóa nguy cơ cháy toàn cầu tăng gấp 2–3 lần vào giai đoạn 2081–2090; Tian và cộng sự (2011) sử dụng mô hình PRECIS và chỉ số FWI dự báo mùa cháy tại Đông Bắc Trung Quốc kéo dài thêm 21–26 ngày vào năm 2080 theo kịch bản A2 và B2.
Tại Việt Nam, các công trình tiền đề đã xác lập nền móng quan trọng: Phạm Ngọc Hưng (1988) hiệu chỉnh chỉ số Nesterov với ngưỡng mưa triệt tiêu 5 mm và phát triển phương pháp số ngày khô hạn liên tục ($H$). Võ Đình Tiến (1995) thiết lập chỉ số dự báo tháng tại Bình Thuận tích hợp 6 yếu tố khí hậu và nhân sinh. Đề tài cấp nhà nước KC.24 (2002–2005) do Vương Văn Quỳnh, Trần Quang Bảo và cộng sự thực hiện đã tiên phong ứng dụng ảnh vệ tinh MODIS, GIS và cơ sở dữ liệu thời tiết của 95 trạm khí tượng để vận hành phần mềm cảnh báo cháy rừng tự động toàn quốc. Nguyễn Đăng Quế và Đặng Văn Thắng (2010), Bế Minh Châu (2011) bước đầu thăm dò tác động của BĐKH đến cháy rừng ở phạm vi cấp tỉnh (Sơn La) hoặc một số khu vực cục bộ.
Tranh luận học thuật cốt lõi nảy sinh giữa hai quan điểm: (1) Nhóm nghiên cứu cổ điển cho rằng có thể dùng chuỗi số liệu khí hậu quá khứ trung bình nhiều năm để xác định cấp cháy rừng cố định theo lịch mùa; (2) Nhóm quan điểm hiện đại (Bế Minh Châu 2001, Vương Văn Quỳnh 2012) chứng minh rằng tại các khu vực chịu ảnh hưởng của hoàn lưu gió mùa xen kẽ khối khí biển và lục địa, các chỉ số tĩnh $P$ hay $H$ bộc lộ sai số lớn, không phản ánh được sự biến đổi động của độ ẩm vật liệu cháy và hoàn toàn bất khả tri trước các kịch bản BĐKH tương lai.
Luận án của Lê Sỹ Doanh định vị chính xác tại giao điểm đột phá: chuyển đổi từ các mô hình dự báo tác nghiệp ngắn hạn sang mô hình hóa tác động BĐKH dài hạn toàn quốc. Công trình khắc phục triệt để hạn chế của các nghiên cứu trước bằng cách xây dựng một chỉ số khí hậu tích lũy mới ($Q_i$), thiết lập mối tương quan thực nghiệm với số ngày cháy nguy hiểm ($S_{nc45}$) và lần đầu tiên bản đồ hóa nguy cơ cháy rừng theo kịch bản B2 cho toàn bộ lãnh thổ Việt Nam tới năm 2090.
Đó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 Lý thuyết Tam giác lửa (Fire Triangle Theory của Brown 1979, Belop 1982, Chandler 1983) và lý thuyết bay hơi vật liệu cháy của Nesterov trong bối cảnh nhiệt đới ẩm. Luận án chỉ ra rằng cơ chế khô kiệt của vật liệu cháy dưới tán rừng nhiệt đới không chỉ chịu sự chi phối tức thời của thời tiết mà mang tính quán tính nhiệt ẩm (hygrothermal inertia).
Mô hình chỉ số khí hậu $Q_i$ được luận án thiết lập dưới dạng hàm trọng số phân rã theo thời gian:
$$Q_i = 0,1 \cdot \left(K_{i-2} \cdot T_{i-2} \cdot |R_{i-2}-100|^{0,8}\right) + 0,2 \cdot \left(K_{i-1} \cdot T_{i-1} \cdot |R_{i-1}-100|^{0,8}\right) + 0,7 \cdot \left(K_i \cdot T_i \cdot |R_i-100|^{0,8}\right)$$
Trong đó:
- $T_i, T_{i-1}, T_{i-2}$: Nhiệt độ không khí trung bình tại ngày/tháng thứ $i$, $i-1$, $i-2$ ($^\circ\text{C}$).
- $R_i, R_{i-1}, R_{i-2}$: Độ ẩm không khí tương đối tương ứng (%).
- $K_i, K_{i-1}, K_{i-2}$: Hệ số hiệu chỉnh theo lượng mưa ngày/tháng.
- Lũy thừa $0,8$ phản ánh mối quan hệ phi tuyến giữa độ hụt bão hòa không khí và tốc độ mất nước của thảm khô. Trọng số $0,1 - 0,2 - 0,7$ lượng hóa chính xác mức độ suy giảm ảnh hưởng của điều kiện khí hậu quá khứ đến độ ẩm vật liệu cháy hiện tại.
Mô hình lý thuyết thiết lập mối quan hệ tuyến tính xác thực giữa số ngày nguy cơ cháy cao/rất cao ($S_{nc45}$) với chỉ số $Q_i$:
$$S_{nc45} = 7,284 \cdot Q_i + 1,029 \quad (R^2 = 0,588)$$
Đây là bước chuyển biến luận giải thực nghiệm (empirical paradigm shift), cung cấp cầu nối toán học trực tiếp giữa dữ liệu kịch bản khí hậu vĩ mô và chỉ số nguy cơ cháy sinh thái vi mô.
+-----------------------------------------------------------------------------------+
| KHUNG PHÂN TÍCH TỔNG THỂ |
+-----------------------------------------------------------------------------------+
| Kịch bản BĐKH B2 (IPCC) -> Nhiệt độ (T), Lượng mưa (P), Độ ẩm (R) tại 95 trạm |
| Chỉ số khí hậu phản ánh nguy cơ cháy rừng (Qi) tích lũy |
| Mô hình thực nghiệm: Snc45 = 7,284 * Qi + 1,029 (R² = 0,588) |
| Phân cấp nguy cơ khí hậu 5 bậc (I: <10d -> V: >25d/tháng) |
| Tích hợp Ma trận Trọng số Trạng thái rừng (fij, Ect) |
| (Khối lượng thảm khô Mtk, thảm tươi Mtt, độ ẩm W13) |
| Phân vùng Nguy cơ Cháy rừng Tích hợp Không gian - Thời gian (2020, 2050, 2090) |
+-----------------------------------------------------------------------------------+
Khung phân tích độc đáo
Khung phân tích của luận án tích hợp đa chiều ba hệ thống lý thuyết:
- Lý thuyết Động lực học Khí hậu (Climate Dynamics): Tiếp cận kịch bản phát thải trung bình B2 của IPCC công bố năm 2009 cho Việt Nam, chi tiết hóa cho 8 vùng khí hậu.
- Lý thuyết Sinh thái Nhiệt ẩm Vật liệu cháy (Fuel Moisture Ecology): Xác định cấu trúc lý hóa thảm che phủ mặt đất thông qua chỉ số hiệu quả canh tác thảm khô/thảm tươi ($E_{ct}$) và chỉ số tích lũy vật liệu cháy theo trạng thái ($f_{ij}$).
- Lý thuyết Quản lý Rủi ro Thích ứng (Adaptive Risk Management): Tích hợp không gian địa lý qua ma trận ghép tầng giữa cấp nguy cơ thời tiết ($S_{nc45}$) và cấp dễ cháy của từng loại rừng.
Điều kiện biên (boundary conditions) của khung phân tích: Áp dụng cho các hệ sinh thái rừng đồi núi và ngập phèn tại Việt Nam chịu sự chi phối của khí hậu nhiệt đới gió mùa; kịch bản khí hậu nền tảng là B2 (phát thải trung bình phù hợp với đặc thù phát triển kinh tế - xã hội, dân số và công nghệ của Việt Nam); ngưỡng lượng mưa triệt tiêu nguy cơ cháy được hiệu chỉnh cục bộ theo từng vùng sinh thái.
Phương pháp nghiên cứu tiên tiến
Thiết kế nghiên cứu
Luận án vận dụng thế giới quan thực chứng luận (positivism) kết hợp cách tiếp cận hệ thống (systems approach), tiếp cận đa ngành (lâm sinh, khí tượng thủy văn, toán tin học, sinh thái cảnh quan) và tiếp cận nghiên cứu phát triển.
Thiết kế nghiên cứu đa tầng (multi-level design) bao gồm:
- Tầng vĩ mô: Phân tích chuỗi số liệu khí hậu quan trắc nhiều năm và dữ liệu dự báo BĐKH từ mô hình khí hậu khu vực độ phân giải cao CCAM tại 95 trạm khí tượng quốc gia đại diện cho 8 vùng sinh thái.
- Tầng trung mô: Không gian hóa các yếu tố nhiệt - ẩm - mưa và số ngày nguy cơ cháy cao $S_{nc45}$ trên nền tảng bản đồ địa hình, độ cao ($Dc$), vĩ độ ($Vd$), kinh độ ($Kd$) sử dụng công nghệ GIS.
- Tầng vi mô: Khảo sát thực địa trên hệ thống ô tiêu chuẩn định vị tại các trạng thái rừng trọng điểm để cân đo sinh khối thảm khô, thảm tươi và độ ẩm vật liệu cháy lúc 13 giờ ($W_{13}$).
Quy trình nghiên cứu rigorous
Quy trình thực hiện được chuẩn hóa nghiêm ngặt qua 4 bước:
- Bước 1: Chuẩn hóa cơ sở dữ liệu khí hậu và điều tra rừng: Thu thập dữ liệu khí tượng bề mặt (nhiệt độ trung bình, nhiệt độ tối cao, tối thấp, lượng mưa, độ ẩm không khí) qua các thập kỷ; điều tra cấu trúc tầng tán, thảm mục tại các ô tiêu chuẩn đại diện cho rừng tự nhiên lá rộng thường xanh, rừng rụng lá/khộp, rừng thông, rừng keo, rừng tràm.
- Bước 2: Phát triển và tối ưu hóa toán học chỉ số khí hậu $Q_i$: Thử nghiệm các tổ hợp tham số trọng số ($\alpha, \beta, \mu$), xác định hàm tối ưu liên hệ giữa các giá trị khí tượng lùi thời gian ($i, i-1, i-2$) với độ khô kiệt vật liệu cháy.
- Bước 3: Xây dựng mô hình hồi quy và kiểm định độ phù hợp: Thiết lập tương quan thực nghiệm giữa $Q_i$ và số ngày có nguy cơ cháy cao thực tế $S_{nc45}$; phân cấp 5 mức độ nguy hiểm:
- Cấp I (Ít khả năng cháy): $S_{nc45} < 10\text{ ngày/tháng}$
- Cấp II (Nguy cơ cháy thấp): $S_{nc45} = 10 - 15\text{ ngày/tháng}$
- Cấp III (Nguy cơ cháy trung bình): $S_{nc45} = 16 - 20\text{ ngày/tháng}$
- Cấp IV (Nguy cơ cháy cao): $S_{nc45} = 21 - 25\text{ ngày/tháng}$
- Cấp V (Nguy cơ cháy rất cao): $S_{nc45} > 25\text{ ngày/tháng}$
- Bước 4: Chiếu xạ kịch bản BĐKH B2 và tích hợp ma trận loại rừng: Đưa dữ liệu nhiệt ẩm dự báo các năm 2020, 2050, 2090 vào mô hình $Q_i \to S_{nc45}$, nhân chập không gian với bản đồ phân bố các trạng thái rừng để phân vùng trọng điểm cháy.
Tam giác đạc dữ liệu (data triangulation) được thực hiện thông qua việc đối chiếu giữa: (a) Số liệu quan trắc khí tượng bề mặt; (b) Dữ liệu đo đạc thực nghiệm độ ẩm thảm mục tại ô tiêu chuẩn; (c) Thống kê lịch sử 9.689 vụ cháy rừng thực tế của Cục Kiểm lâm giai đoạn 2002–2011; (d) Ảnh viễn thám MODIS phát hiện điểm nhiệt độ cao bất thường.
Data và phân tích
Cơ sở dữ liệu thống kê bao quát toàn diện các yếu tố:
- Chuỗi số liệu khí tượng quan trắc từ mạng lưới trạm Khí tượng Quốc gia phân bổ khắp 8 vùng sinh thái: Đông Bắc, Tây Bắc, Đồng bằng Bắc Bộ, Bắc Trung Bộ, Nam Trung Bộ, Tây Nguyên, Đông Nam Bộ, Tây Nam Bộ.
- Bộ kịch bản BĐKH trung bình B2 do Bộ Tài nguyên và Môi trường công bố năm 2009.
- Dữ liệu thực nghiệm thực địa: đo khối lượng thảm khô ($M_{tk}$, kg/ha), khối lượng thảm tươi ($M_{tt}$, kg/ha), độ ẩm thảm khô lúc 13 giờ ($W_{13}$, %).
- Phần mềm phân tích: Xử lý tương quan hồi quy và kiểm định thống kê toán học; tích hợp công nghệ GIS để nội suy không gian và hiển thị hóa bản đồ chuyên đề.
Phát hiện đột phá và implications
Những phát hiện then chốt
Luận án mang lại 4 phát hiện mang tính đột phá về quy luật cháy rừng dưới tác động của BĐKH tại Việt Nam:
-
Sự gia tăng có quy luật và không đồng tốc của số ngày có nguy cơ cháy rừng cao ($S_{nc45}$):
Tính trung bình toàn quốc vào thời điểm năm 2090 theo kịch bản B2, số ngày có nguy cơ cháy rừng cao trong năm đạt tới 84 ngày/năm với hệ số biến động giữa các tỉnh là 41%. Mức nguy cơ này tăng 6 ngày/năm so với mốc năm 2050 và tăng 12 ngày/năm so với mốc năm 2030, khẳng định tác động khuếch đại lũy tiến của quá trình ấm lên toàn cầu lên tính dễ cháy của rừng Việt Nam.
-
Sự phân hóa sâu sắc về không gian và xác lập hai tâm điểm cháy cực đoan:
Hệ số biến động số ngày nguy cơ cháy cao giữa các vùng sinh thái lên tới 36%. Trong đó:
- Tây Nam Bộ là khu vực có nguy cơ cháy rừng cao nhất cả nước với dự báo đạt tới 123 ngày/năm vào năm 2090.
- Tây Nguyên là tâm điểm cháy lớn thứ hai với 101 ngày/năm vào năm 2090 do hiện tượng hạn hán kéo dài và suy giảm lượng mưa mùa khô nghiêm trọng.
- Trái lại, các khu vực như Đông Bắc có số ngày nguy cơ cháy tăng chậm hơn do lượng mưa một số tháng mùa hè có xu hướng tăng nhẹ bù trừ nhiệt độ.
+-----------------------------------------------------------------------------------+
| DỰ BÁO BIẾN THIÊN SỐ NGÀY NGUY CƠ CHÁY CAO (Snc45) ĐẾN NĂM 2090 (KỊCH BẢN B2) |
+------------------------------------+---------------------+------------------------+
| Vùng Sinh Thái Lâm Nghiệp | Snc45 (ngày/năm) | Mức độ gia tăng rủi ro |
+------------------------------------+---------------------+------------------------+
| Tây Nam Bộ | 123 ngày | Cực kỳ nghiêm trọng |
| Tây Nguyên | 101 ngày | Cực kỳ nghiêm trọng |
| Nam Trung Bộ | 92 ngày | Nghiêm trọng |
| Bắc Trung Bộ | 86 ngày | Nghiêm trọng |
| Trung bình toàn quốc | 84 ngày | Gia tăng lũy tiến |
| Tây Bắc Bộ | 79 ngày | Trung bình - Cao |
| Đông Nam Bộ | 76 ngày | Trung bình - Cao |
| Đông Bắc Bộ | 62 ngày | Tăng cục bộ mùa khô |
| Đồng Bằng Bắc Bộ | 58 ngày | Thấp - Trung bình |
+------------------------------------+---------------------+------------------------+
-
Hiện tượng nghịch lý khí hậu - vật liệu cháy (Counter-intuitive Phenomenon):
Tại một số tiểu vùng sinh thái thuộc Đông Bắc và Bắc Trung Bộ, mặc dù tổng lượng mưa năm dự báo có xu hướng tăng 2–5%, nhưng nguy cơ cháy rừng không hề suy giảm mà phân cực gay gắt. Lượng mưa gia tăng chủ yếu tập trung vào các trận mưa bão cực đoan trong mùa mưa (gây lũ lụt), trong khi mùa khô lại bị kéo dài và nền nhiệt tăng vọt ($+1,6^\circ\text{C}$ đến $+2,6^\circ\text{C}$), khiến tốc độ bốc thoát hơi nước của thảm thực vật diễn ra cực nhanh, làm độ ẩm vật liệu cháy ($W_{13}$) tụt xuống dưới ngưỡng bén lửa ($< 15%$) sớm hơn từ 15–20 ngày so với quy luật tự nhiên trước đây.
-
Sự kéo dài và dịch chuyển mùa cháy rừng lịch sử:
Mùa cháy rừng truyền thống tại Tây Bắc (tháng 11 đến tháng 4 năm sau) và Bắc Trung Bộ có xu hướng bắt đầu sớm hơn 10–15 ngày và kết thúc muộn hơn. Các tháng 3 và tháng 4 ghi nhận sự tăng vọt về nhiệt độ kết hợp hiệu ứng gió Phơn Tây Nam khô nóng làm cho số ngày đạt cấp IV và cấp V chiếm tới 70–80% thời gian tháng.
Implications đa chiều
- Về mặt lý thuyết: Mở rộng khung lý thuyết sinh thái học lửa nhiệt đới, chứng minh rằng mô hình dự báo cháy phải tích hợp tính trễ nhiệt ẩm (thermal-moisture lag) và cấu trúc thảm thực vật thay vì chỉ dùng các biến số khí tượng đơn lẻ.
- Về phương pháp luận: Cung cấp quy trình 4 bước chuẩn hóa từ tính toán chỉ số $Q_i$, thiết lập hàm tương quan thực nghiệm $S_{nc45}$ đến chiết xuất bản đồ phân cấp nguy cơ, có khả năng chuyển giao và áp dụng cho các quốc gia trong khu vực Đông Nam Á có cùng hoàn cảnh gió mùa.
- Về ứng dụng thực tiễn và chính sách:
- Cung cấp cơ sở khoa học định lượng để Bộ Nông nghiệp và PTNT, Cục Kiểm lâm tái cơ cấu phân bổ ngân sách và trang thiết bị chữa cháy rừng trọng điểm theo không gian ưu tiên (dồn nguồn lực cho Tây Nam Bộ và Tây Nguyên).
- Tích hợp kết quả phân vùng nguy cơ cháy vào Quy hoạch Lâm nghiệp Quốc gia và Đề án thực hiện Quyết định số 158/2008/QĐ-TTg, Quyết định số 1474/QĐ-TTg của Thủ tướng Chính phủ về ứng phó với BĐKH.
Limitations và Future Research
Mặc dù đạt được những đóng góp học thuật và thực tiễn to lớn, luận án duy trì sự trung thực khoa học khi ghi nhận các giới hạn nghiên cứu:
- Giới hạn kịch bản khí hậu: Nghiên cứu tập trung chuyên sâu vào kịch bản biến đổi khí hậu trung bình B2 của Bộ Tài nguyên và Môi trường (2009). Chưa thực hiện mô phỏng so sánh đa kịch bản với kịch bản phát thải cao A2 hoặc các kịch bản nồng độ phát thải mới nhất theo cấu trúc RCP (RCP 4.5, RCP 8.5) của Báo cáo IPCC AR5.
- Độ phân giải không gian khí tượng: Số liệu khí hậu dựa trên 95 trạm quan trắc bề mặt, tại một số vùng địa hình chia cắt phức tạp (như núi cao hiểm trở Hoàng Liên Sơn hay Trường Sơn), mạng lưới trạm còn thưa, đòi hỏi các thuật toán nội suy không gian có thể tạo ra sai số cục bộ.
- Định lượng hóa biến số nhân sinh: Hoạt động nhân sinh (đốt nương làm rẫy, du lịch, săn bắt) mới được xem xét gián tiếp qua chỉ số tích hợp loại rừng và vùng kinh tế xã hội, chưa xây dựng được mô hình hành vi động của cộng đồng cư dân vùng đệm tiếp giáp rừng.
Chương trình nghiên cứu tương lai (Future Research Agenda):
- Ứng dụng các kịch bản phát thải mới nhất (CMIP6/SSPs và RCPs) để tái đánh giá độ nhạy của chỉ số $Q_i$.
- Tích hợp ảnh vệ tinh viễn thám độ phân giải không gian - thời gian cao (Sentinel-2, Landsat-8/9) để cập nhật liên tục sinh khối và độ ẩm thảm khô ($W_{13}$) theo thời gian thực.
- Thiết lập hệ thống WebGIS tự động kết nối mạng lưới trạm khí tượng tự động chuyên dùng ngành lâm nghiệp và cảm biến IoT tại các khu rừng trọng điểm cháy.
- Phát triển mô hình đa tác tử (Agent-Based Modeling) để lượng hóa tác động tương hỗ giữa BĐKH, thay đổi sinh kế cộng đồng và nguy cơ phát sinh nguồn lửa.
Tác động và ảnh hưởng
Luận án của Lê Sỹ Doanh xác lập dấu ấn học thuật và thực tiễn sâu rộng:
- Tác động học thuật: Trở thành tài liệu tham khảo nền tảng tại Trường Đại học Lâm nghiệp, Viện Khoa học Lâm nghiệp Việt Nam và các cơ sở đào tạo sau đại học chuyên ngành Lâm sinh, Quản lý Tài nguyên Rừng và Biến đổi Khí hậu. Phương pháp luận xây dựng chỉ số $Q_i$ mở ra hướng nghiên cứu mới cho hàng loạt luận văn, luận án kế thừa.
- Chuyển giao và tác động ngành: Cung cấp cơ sở khoa học trực tiếp cho Cục Kiểm lâm và Chi cục Kiểm lâm các tỉnh (tiêu biểu như Thanh Hóa, Phú Thọ, Lâm Đồng, Kiên Giang, Cà Mau) nâng cấp phần mềm cảnh báo sớm nguy cơ cháy rừng trên nền WebGIS.
- Định hình chính sách vĩ mô: Đóng góp luận cứ định lượng cho Kế hoạch hành động ứng phó với Biến đổi khí hậu ngành Nông nghiệp và PTNT giai đoạn 2011–2015, tầm nhìn 2050 (Quyết định số 543/QĐ-BNN-KHCN) và Đề án giảm phát thải khí nhà kính trong nông nghiệp, nông thôn (Quyết định số 3199/QĐ-BNN-KHCN).
- Lợi ích xã hội và môi trường: Giúp các địa phương chủ động xây dựng phương án PCCCR theo phương châm "4 tại chỗ", giảm thiểu diện tích rừng bị thiêu rụi, bảo vệ tính mạng và sinh kế của hàng triệu người dân sống phụ thuộc vào rừng, bảo tồn đa dạng sinh học và giảm hàng triệu tấn khí thải nhà kính ra môi trường.
Đối tượng hưởng lợi
+-----------------------------------------------------------------------------------+
| MA TRẬN CÁC ĐỐI TƯỢNG HƯỞNG LỢI |
+-----------------------------------+-----------------------------------------------+
| Nhóm đối tượng | Giá trị khoa học và ứng dụng thực tiễn |
+-----------------------------------+-----------------------------------------------+
| Nghiên cứu sinh & Nhà khoa học | Hệ thống công thức toán học Qi - Snc45; |
| (Lâm nghiệp, BĐKH, Địa lý) | Khung phương pháp đánh giá tác động 4 bước. |
+-----------------------------------+-----------------------------------------------+
| Nhà quản lý & Cơ quan hoạch định | Bản đồ phân vùng nguy cơ cháy rừng 2020-2090; |
| (Bộ NN&PTNT, Cục Kiểm lâm) | Cơ sở phân bổ kinh phí và lực lượng PCCCR. |
+-----------------------------------+-----------------------------------------------+
| Lực lượng Kiểm lâm địa phương | Dự báo chính xác cấp cháy theo trạng thái |
| (Vườn Quốc gia, Khu Bảo tồn) | rừng và mùa khô hạn cục bộ. |
+-----------------------------------+-----------------------------------------------+
| Chủ rừng, Doanh nghiệp Lâm nghiệp | Thiết kế băng cản lửa, chọn loài cây phòng |
| & Cộng đồng ven rừng | cháy sinh học thích ứng BĐKH bền vững. |
+-----------------------------------+-----------------------------------------------+
- Nghiên cứu sinh và Giảng viên: Khai thác khung phương pháp luận hệ thống và các hàm hồi quy thực nghiệm để mở rộng nghiên cứu sinh thái lửa và quản trị tài nguyên dưới tác động BĐKH.
- Lãnh đạo ngành Lâm nghiệp và Quản lý Nhà nước: Nắm bắt bức tranh toàn cảnh về sự dịch chuyển các điểm nóng cháy rừng trên toàn quốc đến năm 2090, từ đó ban hành các chính sách đầu tư trang thiết bị, quy hoạch vùng đệm và chiến lược thích ứng dài hạn.
- Cán bộ tác nghiệp tại hiện trường: Ứng dụng bảng phân cấp 5 bậc nguy cơ cháy rừng ($S_{nc45}$) kết hợp đặc điểm thảm mục để triển khai tuần tra, ứng trực chữa cháy chính xác vào các tháng cao điểm.
- Cộng đồng và Doanh nghiệp phát triển rừng: Nâng cao năng lực thích ứng thông qua việc áp dụng các giải pháp lâm sinh: trồng băng cây xanh phòng cháy, tỉa cành, xử lý thực bì có kiểm soát và chuyển đổi cơ cấu cây trồng chịu lửa cao.
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à nó 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 Chỉ số khí hậu phản ánh nguy cơ cháy rừng $Q_i$ tích hợp tính trễ nhiệt ẩm và phương trình tương quan thực nghiệm:
$$S_{nc45} = 7,284 \cdot Q_i + 1,029 \quad (R^2 = 0,588)$$
Công trình đã mở rộng Lý thuyết Tam giác lửa (Brown 1979) và Lý thuyết bốc hơi Nesterov (1949) bằng cách chứng minh rằng trong điều kiện nhiệt đới gió mùa, tính dễ cháy của thảm mục là một quá trình tích lũy phi tuyến tính với độ hụt bão hòa không khí dạng $|R-100|^{0,8}$ và suy giảm ảnh hưởng theo chuỗi thời gian lùi ($0,1 - 0,2 - 0,7$). Điều này giải quyết triệt để hạn chế của các chỉ số khí tượng tức thời vốn không có khả năng liên kết với các mô hình kịch bản BĐKH toàn cầu.
2. Điểm đổi mới phương pháp luận khi so sánh với ít nhất 2 nghiên cứu quốc tế?
- So với nghiên cứu của Nicholls và cộng sự (2005) tại Đông Nam Úc: Nghiên cứu của Nicholls sử dụng chỉ số FFDI đòi hỏi đầu vào chi tiết về tốc độ gió hàng ngày và nhân tố khô hạn đất $D$, vốn rất khó dự báo chính xác trong các kịch bản BĐKH dài hạn. Luận án của Lê Sỹ Doanh đã tinh gọn hóa vào các biến số nhiệt - ẩm - mưa cốt lõi nhưng tích lũy đa thời gian ($Q_i$), cho phép tích hợp trực tiếp vào kết quả của mô hình khí hậu vùng CCAM mà vẫn đảm bảo độ tin cậy thực nghiệm ($R^2 = 0,588$).
- So với nghiên cứu của Chowdhury và Hassan (2011, 2013) tại Canada: Nghiên cứu tại Alberta phụ thuộc chủ yếu vào dữ liệu viễn thám tức thời (MODIS NMDI, TVWI) để phát hiện trạng thái cháy hiện tại mà không dự báo được xu thế biến đổi khí hậu đến cuối thế kỷ. Luận án đã kết hợp hài hòa giữa kịch bản BĐKH vĩ mô (IPCC B2) với đo đạc thực địa thảm khô/thảm tươi ($M_{tk}, M_{tt}, W_{13}$) và không gian hóa bằng GIS trên quy mô toàn bộ 8 vùng sinh thái của một quốc gia.
3. Phát hiện bất ngờ nhất (surprising finding) được hỗ trợ bởi số liệu là gì?
Phát hiện bất ngờ nhất là nghịch lý phân hóa nguy cơ cháy tại Tây Nam Bộ và Đông Bắc Bộ:
Mặc dù Tây Nam Bộ là vùng đồng bằng trũng ngập nước và rừng tràm ngập phèn, nhưng dưới tác động của kịch bản B2, đây lại là khu vực có số ngày nguy cơ cháy cao lớn nhất cả nước (123 ngày/năm vào năm 2090), vượt xa vùng khô hạn truyền thống Nam Trung Bộ (92 ngày/năm). Ngược lại, vùng Đông Bắc có lượng mưa mùa hè tăng làm giảm số ngày nguy cơ cháy trung bình năm (62 ngày/năm), nhưng lại tạo ra sinh khối thảm khô cực lớn vào cuối mùa mưa, biến các đợt rét khô hạn vào mùa đông thành các sự kiện cháy bùng phát dữ dội không theo quy luật cũ.
4. Luận án có cung cấp quy trình nhân bản (replication protocol) không?
Có. Luận án cung cấp một quy trình 4 bước chuẩn hóa và định lượng hoàn toàn:
- Xác định ma trận chỉ số khí hậu $Q_i$ từ bộ dữ liệu nhiệt độ, độ ẩm và lượng mưa theo công thức trọng số lũy thừa.
- Tính toán số ngày nguy cơ cháy cao $S_{nc45}$ theo phương trình hồi quy tuyến tính chuẩn.
- Phân cấp 5 bậc nguy cơ thời tiết theo ngưỡng số ngày cháy trong tháng ($<10$, $10–15$, $16–20$, $21–25$, $>25$ ngày).
- Ghép tầng ma trận với bản đồ hiện trạng rừng và chỉ số vật liệu cháy ($f_{ij}, E_{ct}$) để xuất ra bản đồ nguy cơ cháy rừng tổng hợp trên phần mềm GIS.
5. Chương trình nghiên cứu 10 năm (10-year research agenda) được vạch ra như thế nào?
Chương trình 10 năm tập trung vào 4 trọng tâm:
- Nâng cấp mô hình với các kịch bản nồng độ phát thải mới (CMIP6/SSPs và RCP 4.5, 8.5).
- Tích hợp dữ liệu viễn thám vệ tinh siêu phổ và viễn thám radar để cập nhật bản đồ vật liệu cháy động với chu kỳ 5 ngày/lần.
- Nghiên cứu sâu về đặc tính bắt lửa hóa sinh của từng loài cây bản địa phục vụ thiết kế các đai rừng phòng hộ cản lửa sinh học thích ứng BĐKH.
- Xây dựng nền tảng cảnh báo cháy rừng trực tuyến cấp quốc gia tích hợp trí tuệ nhân tạo (AI) và mạng lưới cảm biến mặt đất (IoT).
Kết luận
Luận án Tiến sĩ của Lê Sỹ Doanh là công trình nghiên cứu khoa học tiên phong, toàn diện và có hệ thống đầu tiên tại Việt Nam về tác động của biến đổi khí hậu đến nguy cơ cháy rừng. Các đóng góp cốt lõi của công trình được tóm lược qua 5 điểm nhấn then chốt:
- Thiết lập thành công Chỉ số khí hậu phản ánh nguy cơ cháy rừng $Q_i$: Lần đầu tiên xây dựng công thức toán học lượng hóa chính xác ảnh hưởng tích lũy phi tuyến tính của nhiệt ẩm quá khứ đến nguy cơ cháy rừng nhiệt đới.
- Xác lập Phương trình hồi quy thực nghiệm $S_{nc45}$: Tìm ra quy luật tương quan chặt chẽ $S_{nc45} = 7,284 \cdot Q_i + 1,029$ ($R^2 = 0,588$), cho phép chuyển đổi trực tiếp dữ liệu kịch bản BĐKH thành số ngày nguy cơ cháy nguy hiểm.
- Hoàn thiện Phương pháp luận đánh giá tác động BĐKH 4 bước: Xây dựng quy trình tích hợp giữa mô hình khí hậu khu vực CCAM, công nghệ không gian địa lý GIS và cấu trúc vật liệu cháy thực địa ($M_{tk}, M_{tt}, W_{13}$).
- Bản đồ hóa định lượng nguy cơ cháy rừng toàn quốc đến năm 2090: Dự báo chính xác số ngày cháy cao trung bình toàn quốc đạt 84 ngày/năm vào năm 2090, đồng thời xác định Tây Nam Bộ (123 ngày/năm) và Tây Nguyên (101 ngày/năm) là hai tâm điểm cháy khốc liệt nhất.
- Hệ thống hóa các giải pháp ứng phó lâm sinh phân vùng: Đề xuất gói giải pháp kỹ thuật lâm sinh, quản lý vật liệu cháy và phòng cháy sinh học đặc thù cho từng vùng sinh thái nhằm chủ động thích ứng với BĐKH.
Công trình tạo tiền đề khai mở 3 hướng nghiên cứu học thuật mới: (1) Mô hình hóa sinh thái lửa thích ứng BĐKH đa kịch bản; (2) Ứng dụng công nghệ không gian địa lý viễn thám thời gian thực trong kiểm soát vật liệu cháy rừng; (3) Kinh tế sinh thái và quản trị rủi ro thiên tai rừng trong bối cảnh cam kết giảm phát thải ròng toàn cầu. Luận án khẳng định giá trị khoa học bền vững và là kim chỉ nam thực tiễn cho chiến lược bảo vệ tài nguyên rừng Việt Nam qua các thập kỷ tới.
