Tổng quan nghiên cứu

Biến đổi khí hậu đang tạo ra những áp lực chưa từng có lên hệ sinh thái rừng nhiệt đới tại Việt Nam. Theo thống kê của Cục Kiểm lâm, chỉ riêng trong năm 2007, cả nước đã xảy ra 773 vụ cháy rừng, tăng 46,4% so với năm 2006, làm thiệt hại 4.249,12 ha rừng. Sang năm 2008, toàn quốc tiếp tục ghi nhận 146 vụ cháy với diện tích rừng bị tàn phá là 693 ha. Báo cáo của Ủy ban Liên chính phủ về Biến đổi Khí hậu (IPCC) chỉ ra rằng nhiệt độ bề mặt Trái Đất trong giai đoạn 1996–2005 đã tăng 0,74°C, kéo theo tần suất nắng nóng kéo dài và sự thất thường của lượng mưa. Thực trạng này đặt ra yêu cầu cấp bách phải đánh giá định lượng mối quan hệ giữa các nhân tố thời tiết cực đoan và nguy cơ hỏa hoạn tại các vùng sinh thái trọng điểm.

Luận văn thạc sĩ tập trung giải quyết vấn đề cốt lõi: Xác định quy luật biến động của các yếu tố khí tượng và định lượng mức độ ảnh hưởng của biến đổi khí hậu đến nguy cơ cháy rừng. Mục tiêu cụ thể của công trình gồm:

  • Thử nghiệm phương pháp đánh giá định lượng tác động của biến đổi khí hậu đến nguy cơ cháy rừng thông qua chỉ số khí tượng tổng hợp.
  • Đánh giá mức độ biến động của nguy cơ cháy rừng tại hai vùng sinh thái trọng điểm là Tây Bắc (đại diện bởi trạm Sơn La và Lai Châu) và Bắc Trung Bộ (đại diện bởi trạm Thanh Hóa và Vinh).
  • Dự báo số ngày có nguy cơ cháy ở mức nguy hiểm và cực kỳ nguy hiểm theo các mốc thời gian 2020, 2050, 2070 và 2100 ứng với các kịch bản phát thải khí nhà kính.
  • Đề xuất các giải pháp kỹ thuật và quản lý nhằm giảm thiểu thiệt hại tài nguyên rừng.

Nghiên cứu có ý nghĩa thực tiễn to lớn trong việc nâng cao độ chính xác của công tác dự báo cháy rừng, giúp các cơ quan quản lý chuyển dịch từ thế bị động ứng phó sang chủ động phòng ngừa, hướng tới mục tiêu giảm thiểu từ 20% đến 30% diện tích rừng bị thiệt hại hàng năm do cháy rừng gây ra.

Cơ sở lý thuyết và phương pháp nghiên cứu

Khung lý thuyết áp dụng

Nghiên cứu được xây dựng dựa trên nền tảng lý thuyết dự báo cháy rừng bằng chỉ tiêu khí tượng tổng hợp của V. Nesterov và hệ thống kịch bản biến đổi khí hậu toàn cầu của IPCC:

  • Chỉ tiêu khí tượng tổng hợp P của V. Nesterov: Phương pháp đánh giá nguy cơ cháy dựa trên tích số giữa nhiệt độ không khí lúc 13 giờ ($t_{13}$) và độ chênh lệch bão hòa độ ẩm không khí ($d_{13}$), tích lũy qua các ngày không mưa hoặc có lượng mưa dưới ngưỡng quy định. Tại Việt Nam, công thức được chuẩn hóa với ngưỡng mưa ngày điều chỉnh là 5 mm (thay vì 3 mm như nguyên bản tại Nga). Nguy cơ cháy được phân định thành 5 cấp độ rõ rệt: Cấp I (< 1.000, ít khả năng cháy), Cấp II (1.000 - 2.500, nguy cơ thấp), Cấp III (2.501 - 5.000, nguy cơ trung bình), Cấp IV (5.001 - 10.000, nguy cơ cao/nguy hiểm) và Cấp V (> 10.000, nguy cơ cực kỳ nguy hiểm).
  • Mô hình kịch bản biến đổi khí hậu MAGICC/SCENGEN: Công cụ tổ hợp mô hình chu trình khí quyển và hoàn lưu chung (GCM) do Bộ Khoa học, Công nghệ và Môi trường cung cấp, cho phép mô phỏng sự gia tăng nhiệt độ và biến động lượng mưa tại Việt Nam theo các họ kịch bản phát thải cao, trung bình và thấp (A1, A2, B1, B2) từ năm 2020 đến năm 2100.
  • Các khái niệm then chốt: Độ thiếu hụt bão hòa độ ẩm không khí ($d_{13}$ tính bằng mb), chỉ số khô hạn liên tục $H$ của Phạm Ngọc Hưng, độ ẩm vật liệu cháy và chu kỳ mùa cháy đặc thù (mùa cháy Bắc Trung Bộ từ tháng 3 đến tháng 8; mùa cháy Tây Bắc từ tháng 10 đến tháng 3 năm sau).

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

Nghiên cứu áp dụng quy trình phân tích định lượng kết hợp mô hình hóa thống kê:

  • Cỡ mẫu và nguồn dữ liệu: Bộ dữ liệu gồm 7.665 ngày quan trắc khí tượng liên tục trong giai đoạn 21 năm (1980–2000) thu thập từ 4 trạm khí tượng quốc gia đại diện gồm Thanh Hóa, Vinh, Sơn La và Lai Châu. Dữ liệu tương lai kế thừa từ 6 kịch bản biến đổi khí hậu chi tiết hóa cho Việt Nam tại các mốc năm 2020, 2050, 2070 và 2100.
  • Phương pháp chọn mẫu: Phương pháp chọn mẫu phân tầng có chủ đích (purposive stratified sampling) được lựa chọn nhằm đại diện cho hai tiểu vùng sinh thái có tính chất khí hậu và chế độ gió mùa phân hóa sâu sắc (vùng chịu ảnh hưởng của gió Tây Nam khô nóng tại Bắc Trung Bộ và vùng địa hình núi cao khô hạn mùa đông tại Tây Bắc).
  • Phương pháp phân tích và lý do lựa chọn: Sử dụng toán thống kê đa biến trên phần mềm SPSS và Excel để thiết lập phương trình hồi quy tương quan giữa độ thiếu hụt bão hòa $d_{13}$ với nhiệt độ $t_{13}$ và lượng mưa $M$. Lý do lựa chọn phương pháp này là vì các kịch bản khí hậu tương lai chỉ cung cấp giá trị nhiệt độ và lượng mưa trung bình, không có sẵn dữ liệu độ ẩm chi tiết từng ngày. Việc xây dựng phương trình tương quan chặt ($R > 0,80$) cho phép nội suy chuẩn xác chỉ số $d_{13}$, từ đó tính toán chính xác chỉ số P và phân cấp nguy cơ cháy cho các thập kỷ tương lai.

Kết quả nghiên cứu và thảo luận

Những phát hiện chính

Qua phân tích chuỗi dữ liệu thực nghiệm 21 năm và mô phỏng kịch bản tương lai, nghiên cứu đã rút ra 4 phát hiện quan trọng:

  1. Xu hướng tăng nhiệt độ và suy giảm lượng mưa mùa khô trong hai thập kỷ (1980–2000): Nhiệt độ trung bình lúc 13 giờ tại các trạm nghiên cứu đều có xu hướng tăng từ 0,1°C đến 0,3°C (tăng 0,37% đến 1,13%). Tại Sơn La, nhiệt độ trung bình lúc 13 giờ tăng 0,3°C (từ 25,1°C lên 25,4°C), trong khi lượng mưa các tháng mùa khô (tháng 10 đến tháng 12) giảm mạnh từ 10 đến 50 mm. Tại Vinh, tổng lượng mưa trung bình năm trong thập kỷ 1990 giảm 211,0 mm (tương đương giảm 9,6%) so với thập kỷ 1980, lượng mưa giảm sâu vào các tháng đầu mùa cháy (tháng 3 đến tháng 6).

  2. Chỉ số khí tượng tổng hợp P tăng trưởng đồng loạt trên quy mô vùng: Chỉ số nguy cơ cháy P trung bình năm của cả hai vùng sinh thái đều tăng từ 15% đến 17% qua hai thập kỷ. Khu vực Tây Bắc ghi nhận chỉ số P tăng từ 2.763,2 lên 3.232,4 (tăng 17,0%). Khu vực Bắc Trung Bộ tăng từ 2.639,4 lên 3.035,8 (tăng 15,0%). Tại Sơn La, chỉ số P các tháng đầu mùa cháy (tháng 9 đến tháng 11) tăng vọt tới 45,6%.

  3. Gia tăng mạnh mẽ số ngày có nguy cơ cháy cấp IV và cấp V: Tổng số ngày đạt cấp cực kỳ nguy hiểm (cấp IV và cấp V) tăng rõ rệt ở tất cả các trạm. Trạm Thanh Hóa tăng từ 61 ngày/năm (thập kỷ 1980) lên 78 ngày/năm (thập kỷ 1990), tương ứng mức tăng 27,87%. Trạm Sơn La tăng từ 58 ngày lên 81 ngày/năm, tương ứng mức tăng 39,66%. Trạm Lai Châu có số ngày cấp IV–V cao nhất toàn vùng với trung bình 85 ngày/năm.

  4. Dự báo nguy cơ cháy rừng leo thang theo các kịch bản 2020–2100: Đến năm 2100, nhiệt độ mùa cháy tại khu vực nghiên cứu được dự báo tăng từ 2,0°C đến 4,3°C. Kéo theo đó:

  • Tại Bắc Trung Bộ: Tổng số ngày cấp IV–V tăng từ 2% đến 55%, tổng chỉ tiêu P tăng 2% đến 47%. Đến năm 2100, số ngày nguy hiểm cháy tại Thanh Hóa đạt từ 60 đến 74 ngày (tăng tới 57,4% so với năm 2000), tại Vinh đạt từ 74 đến 89 ngày.
  • Tại Tây Bắc: Tổng số ngày cấp IV–V tăng từ 20% đến 70%, tổng chỉ tiêu P tăng 4% đến 45%. Đến năm 2100, số ngày nguy hiểm cháy tại Sơn La tăng lên 58 đến 71 ngày (tăng tới 97,2% so với mốc năm 2000), tại Lai Châu đạt từ 77 đến 86 ngày.

Thảo luận kết quả

Sự gia tăng đột biến của chỉ tiêu P và số ngày nguy hiểm cháy tại Bắc Trung Bộ bắt nguồn từ cơ chế hoạt động của gió mùa Tây Nam khô nóng (hiệu ứng phơn) kéo dài liên tục từ tháng 3 đến tháng 8, đẩy nhiệt độ không khí lên trên 38°C và hạ thấp độ ẩm tương đối xuống dưới 35%. Tại Tây Bắc, địa hình núi cao chia cắt kết hợp với mùa khô kéo dài từ tháng 10 đến tháng 3 năm sau làm cho độ ẩm vật liệu cháy dưới tán rừng giảm xuống mức bén lửa tối thiểu.

Về mặt trình bày dữ liệu, các phát hiện trên có thể được hệ thống hóa trực quan thông qua:

  • Biểu đồ đường đa trục (Line Chart) so sánh tương quan giữa đà tăng nhiệt độ lúc 13 giờ và sự suy giảm lượng mưa trung bình tháng giữa hai thập kỷ 1980 và 1990.
  • Biểu đồ thanh xếp chồng (Stacked Bar Chart) biểu thị sự dịch chuyển tỷ lệ phần trăm số ngày cháy phân theo 5 cấp độ (từ cấp I đến cấp V) cho 4 trạm quan trắc.
  • Bảng ma trận kịch bản (Scenario Heatmap Matrix) thể hiện số ngày cấp IV–V tương ứng với 3 mức phát thải (thấp, vừa, cao) tại các mốc 2020, 2050, 2070 và 2100.

Kết quả này hoàn toàn nhất quán với công trình nghiên cứu của các tác giả quốc tế tại khu vực Đông Nam Australia khi mô phỏng thời tiết cháy rừng theo các kịch bản phát thải khí nhà kính tương lai. Tại Việt Nam, nghiên cứu đã tiên phong lượng hóa thành công tác động của biến đổi khí hậu đến hỏa hoạn rừng, cung cấp luận cứ khoa học tin cậy cho công tác lâm nghiệp.

Đề xuất và khuyến nghị

Để chủ động giảm thiểu rủi ro cháy rừng trước bối cảnh biến đổi khí hậu phức tạp, luận văn đề xuất 4 nhóm giải pháp trọng tâm:

  1. Hiện đại hóa hệ thống quan trắc và cảnh báo cháy sớm tự động:
  • Động từ hành động: Lắp đặt, tích hợp, tự động hóa mạng lưới trạm đo khí tượng chuyên dụng dưới tán rừng và camera giám sát nhiệt hồng ngoại.
  • Chủ thể thực hiện: Cục Kiểm lâm chủ trì, phối hợp với Sở Nông nghiệp và Phát triển Nông thôn các tỉnh Tây Bắc và Bắc Trung Bộ.
  • Target metric: Rút ngắn thời gian phát hiện điểm cháy xuống dưới 30 phút kể từ khi bén lửa; giảm 25% diện tích rừng bị cháy hàng năm.
  • Timeline: Giai đoạn 2026–2030.
  1. Thiết kế và quy hoạch hệ thống công trình phòng cháy lâm sinh:
  • Động từ hành động: Xây dựng, làm mới và duy tu hệ thống đường băng cản lửa xanh (sử dụng cây bản địa phòng cháy), đường băng cản lửa trắng và các bể chứa nước dự phòng tại các tiểu khu nguy cơ cao.
  • Chủ thể thực hiện: Các Ban Quản lý rừng phòng hộ, Vườn quốc gia và Công ty TNHH MTV Lâm nghiệp.
  • Target metric: Hoàn thành 100% mạng lưới băng cản lửa rộng 15–20 mét tại các khu vực rừng có số ngày cấp IV–V vượt ngưỡng 60 ngày/năm.
  • Timeline: Hoàn thành trước năm 2028.
  1. Ứng dụng mô hình dự báo tương quan khí tượng tổng hợp P vào phần mềm tác nghiệp:
  • Động từ hành động: Số hóa, chuyển giao thuật toán tương quan $d_{13}$ theo nhiệt độ và lượng mưa vào hệ thống thông tin địa lý (GIS) của kiểm lâm địa phương.
  • Chủ thể thực hiện: Chi cục Kiểm lâm các tỉnh Sơn La, Điện Biên, Lai Châu, Thanh Hóa, Nghệ An.
  • Target metric: Nâng độ chính xác của bản tin dự báo cấp cháy rừng hàng ngày lên trên 85%.
  • Timeline: Triển khai thử nghiệm hàng năm vào đầu mỗi mùa cháy (tháng 3 tại Bắc Trung Bộ và tháng 10 tại Tây Bắc).
  1. Kiện toàn năng lực chữa cháy tại chỗ và tập huấn cộng đồng vùng đệm:
  • Động từ hành động: Đào tạo, trang bị công cụ chữa cháy cơ giới nhẹ, diễn tập phương án phối hợp 4 tại chỗ cho lực lượng dân quân tự vệ và người dân nhận khoán bảo vệ rừng.
  • Chủ thể thực hiện: Ủy ban nhân dân cấp huyện, xã phối hợp với Hạt Kiểm lâm sở tại.
  • Target metric: 100% hộ dân sống gần rừng ký cam kết phòng cháy; 100% tổ đội quần chúng được huấn luyện kỹ thuật xử lý thực bì an toàn.
  • Timeline: Thực hiện định kỳ hàng quý trong giai đoạn 2026–2030.

Đối tượng nên tham khảo luận văn

Công trình nghiên cứu mang giá trị học thuật và thực tiễn cao, đặc biệt hữu ích cho 4 nhóm đối tượng sau:

  1. Cán bộ quản lý nhà nước ngành Kiểm lâm và Lâm nghiệp:
  • Lợi ích: Nắm bắt cơ sở khoa học về quy luật biến đổi của thời tiết cháy rừng theo từng vùng sinh thái.
  • Use case: Xây dựng phương án phòng cháy chữa cháy rừng (PCCCR) cấp tỉnh, phân bổ ngân sách mua sắm trang thiết bị và bố trí lực lượng thường trực theo các cấp dự báo P.
  1. Các nhà nghiên cứu, giảng viên và học viên cao học chuyên ngành Lâm nghiệp và Quản lý Môi trường:
  • Lợi ích: Tiếp cận phương pháp luận định lượng mới trong việc kết hợp chỉ tiêu Nesterov với mô hình kịch bản biến đổi khí hậu MAGICC/SCENGEN.
  • Use case: Kế thừa phương trình hồi quy tương quan đa biến ($d_{13}$, $t_{13}$, lượng mưa) để phát triển các mô hình dự báo nguy cơ cháy cho các vùng sinh thái khác như Tây Nguyên hay Đông Bắc Bộ.
  1. Ban Quản lý các Vườn Quốc gia, Khu Bảo tồn Thiên nhiên và Chủ rừng lớn:
  • Lợi ích: Hiểu rõ đặc điểm phân hóa mùa cháy và nguy cơ gia tăng số ngày nguy hiểm cấp IV–V tại từng tiểu vùng.
  • Use case: Xây dựng bản đồ phân vùng rủi ro cháy rừng theo tháng, xác định vị trí tối ưu để bố trí chòi canh lửa, đường băng cản lửa và lịch tuần tra thực địa.
  1. Các tổ chức phi chính phủ (NGO) và cơ quan hoạch định chính sách thích ứng khí hậu:
  • Lợi ích: Khai thác số liệu dự báo định lượng về nguy cơ cháy rừng giai đoạn 2020–2100 theo các kịch bản phát thải.
  • Use case: Thiết kế các dự án tài trợ phát triển sinh kế gắn liền với quản lý lửa rừng dựa vào cộng đồng và bảo tồn đa dạng sinh học trong bối cảnh biến đổi khí hậu toàn cầu.

Câu hỏi thường gặp

1. Chỉ tiêu khí tượng tổng hợp P của V. Nesterov được tính toán và áp dụng tại Việt Nam như thế nào?

Chỉ số P được tính bằng tổng tích số giữa nhiệt độ không khí lúc 13 giờ ($t_{13}$, tính bằng °C) và độ chênh lệch bão hòa độ ẩm không khí ($d_{13}$, tính bằng mb) qua các ngày liên tiếp không mưa. Tại Việt Nam, hệ số điều chỉnh lượng mưa $K = 0$ khi ngày có mưa từ 5 mm trở lên (thay vì 3 mm như ở Nga) và $K = 1$ khi lượng mưa dưới 5 mm. Khi giá trị P vượt ngưỡng 5.000, nguy cơ cháy bước vào cấp IV (nguy hiểm) và trên 10.000 là cấp V (cực kỳ nguy hiểm).

2. Vì sao mùa cháy rừng ở Tây Bắc và Bắc Trung Bộ lại diễn ra vào hai thời điểm hoàn toàn lệch nhau?

Mùa cháy rừng tại Bắc Trung Bộ kéo dài từ tháng 3 đến tháng 8 do chịu tác động trực tiếp của gió Tây Nam khô nóng vượt dãy Trường Sơn, khiến nhiệt độ tăng vọt và độ ẩm tụt giảm mạnh. Ngược lại, mùa cháy tại Tây Bắc diễn ra từ tháng 10 đến tháng 3 năm sau do đây là thời kỳ mùa khô sâu sắc của miền núi phía Bắc, lượng mưa tháng 12 tại Sơn La chỉ đạt khoảng 14,8 mm và thường xuyên chịu ảnh hưởng của gió mùa Đông Bắc hanh khô.

3. Biến đổi khí hậu làm gia tăng số ngày có nguy cơ cháy cấp IV và cấp V trong tương lai ra sao?

Theo các kịch bản phát thải khí nhà kính, đến năm 2100, nhiệt độ mùa cháy tăng từ 2,0°C đến 4,3°C khiến độ ẩm vật liệu cháy suy giảm nhanh hơn. Tính toán cho thấy số ngày cấp IV–V tại Bắc Trung Bộ sẽ tăng từ 2% đến 55% (Vinh đạt tới 89 ngày/năm), trong khi tại Tây Bắc con số này tăng từ 20% đến 70% (Sơn La tăng tới 97,2% so với năm 2000, đạt 71 ngày nguy hiểm), làm kéo dài đáng kể khung thời gian trực cháy.

4. Phương pháp xác định độ thiếu hụt bão hòa độ ẩm không khí ($d_{13}$) trong luận văn có điểm gì nổi bật?

Do các kịch bản biến đổi khí hậu tương lai không cung cấp dữ liệu độ ẩm không khí chi tiết hàng ngày, luận văn đã thiết lập thành công phương trình tương quan hồi quy giữa $d_{13}$ với nhiệt độ 13 giờ ($t_{13}$) và lượng mưa ($M$). Với hệ số tương quan đạt từ 0,81 đến 0,84, phương pháp này cho phép nội suy tin cậy giá trị $d_{13}$, giải quyết triệt để bài toán thiếu dữ liệu độ ẩm trong dự báo dài hạn.

5. Ngoài các yếu tố khí tượng, những nhân tố nào cần được tích hợp để nâng cao độ chính xác dự báo cháy thực địa?

Bên cạnh các chỉ tiêu thời tiết (nhiệt độ, lượng mưa, độ ẩm), mô hình thực địa cần tích hợp thêm: cấu trúc và tải lượng vật liệu cháy dưới tán rừng, đặc điểm địa hình (độ dốc, hướng phơi sườn núi) và các nhân tố kinh tế - xã hội như tập quán đốt nương làm rẫy của người dân địa phương để tạo nên bản đồ dự báo nguy cơ cháy tích hợp hoàn chỉnh.

Kết luận

  • Quy luật biến đổi khí hậu rõ nét: Nhiệt độ trung bình lúc 13 giờ tại hai vùng nghiên cứu đã tăng 0,1°C đến 0,3°C trong giai đoạn 1980–2000, lượng mưa mùa khô giảm từ 10 đến 50 mm/tháng.
  • Chỉ tiêu nguy cơ cháy leo thang: Chỉ số P trung bình năm tăng 15% đến 17% trên toàn vùng; số ngày đạt cấp cực kỳ nguy hiểm (cấp IV–V) tăng 27,87% tại Thanh Hóa và 39,66% tại Sơn La.
  • Tương quan toán học chặt chẽ: Thiết lập thành công mô hình hồi quy tuyến tính xác định độ thiếu hụt bão hòa $d_{13}$ theo nhiệt độ và lượng mưa với hệ số tương quan cao ($R = 0,81 - 0,84$).
  • Dự báo nguy cơ tương lai nghiêm trọng: Đến năm 2100, số ngày nguy cơ cháy cấp IV–V được dự báo tăng tới 55% tại Bắc Trung Bộ và tăng tới 70% tại Tây Bắc theo các kịch bản phát thải khí nhà kính.
  • Hệ thống giải pháp đồng bộ: Đề xuất hoàn chỉnh 4 nhóm giải pháp kỹ thuật lâm sinh, quan trắc tự động, chuyển giao công nghệ dự báo và nâng cao năng lực cộng đồng tại chỗ.

Đóng góp lớn nhất của luận văn là đã đặt nền móng nghiên cứu định lượ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, chuyển hóa các kịch bản phát thải trừu tượng thành các chỉ số quản lý lửa rừng cụ thể.

Trong giai đoạn tiếp theo (2026–2030), hướng nghiên cứu cần mở rộng tích hợp công nghệ viễn thám, ảnh vệ tinh thời gian thực và mô hình hóa không gian GIS để xây dựng hệ thống phần mềm cảnh báo cháy rừng tự động trên toàn quốc. Các nhà quản lý, nhà nghiên cứu và chủ rừng hãy chủ động áp dụng các phương pháp luận và dữ liệu dự báo từ luận văn này vào việc xây dựng kế hoạch quản lý lửa rừng bền vững ngay hôm nay.