Tổng quan nghiên cứu

Biến đổi khí hậu toàn cầu và sự suy giảm diện tích rừng tự nhiên đang đặt ra những thách thức nghiêm trọng đối với công tác bảo tồn sinh thái tại Việt Nam. Theo số liệu công bố năm 2019 của Bộ Nông nghiệp và Phát triển nông thôn, độ che phủ rừng toàn quốc đạt 41,89%, trong khi thống kê giai đoạn 2012 - 2017 cho thấy 11% diện tích rừng tự nhiên bị mất xuất phát từ hành vi chặt phá trái phép và 89% do chuyển đổi mục đích sử dụng. Tỉnh Thừa Thiên Huế sở hữu diện tích tự nhiên 5.025,30 km2 với 288.334,37 ha đất có rừng, đạt độ che phủ 57,37% vào năm 2019. Trong đó, diện tích rừng tự nhiên chiếm 211.373,11 ha và rừng trồng chiếm 76.961,26 ha. Tuy nhiên, chỉ riêng trong năm 2019, địa phương đã ghi nhận biến động giảm 8.450,90 ha rừng, bao gồm 5.409,33 ha khai thác rừng trồng, 226,58 ha do cháy rừng, 10,35 ha bị phá trái phép, 11,26 ha do sạt lở và 2.793,38 ha do các nguyên nhân sinh cảnh khác.

Địa hình Thừa Thiên Huế trải dài phức tạp từ dãy núi Trường Sơn hiểm trở ra đến vùng bờ biển dài 120 km trên địa bàn 8 đơn vị hành chính cấp huyện, thị xã và thành phố. Thực trạng lực lượng kiểm lâm mỏng khiến công tác tuần tra thực địa gặp vô vàn trở ngại. Trước bối cảnh đó, nghiên cứu được thực hiện nhằm xây dựng cơ sở khoa học và quy trình thực nghiệm ứng dụng các chỉ số viễn thám quang học từ vệ tinh Sentinel-2 trên nền tảng Google Earth Engine. Mục tiêu cụ thể là xác định các ngưỡng biến động quang phổ tin cậy để tự động hóa việc phát hiện sớm và giám sát mất rừng giai đoạn 2017 - 2020. Kết quả nghiên cứu giúp rút ngắn chu kỳ theo dõi xuống 5 ngày, nâng cao độ chính xác cảnh báo trên 90% và hỗ trợ đắc lực cho công tác quản lý tài nguyên rừng bền vững.

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 lý thuyết phản xạ phổ điện từ của thảm thực vật và lý thuyết phát hiện biến động đa thời gian trong viễn thám quang học. Đặc tính phản xạ của lá cây khỏe mạnh hấp thụ mạnh bức xạ dải sóng đỏ (Red) do chất diệp lục và phản xạ mạnh ở dải sóng cận hồng ngoại (NIR) do cấu trúc tế bào mô giậu. Khi rừng bị tàn phá hoặc suy thoái, độ phản xạ NIR suy giảm mạnh trong khi độ phản xạ dải sóng hồng ngoại sóng ngắn (SWIR) và dải sóng đỏ tăng cao.

Luận văn vận dụng các khái niệm nền tảng theo Luật Lâm nghiệp năm 2017 và hướng dẫn quốc tế:

  • Mất rừng: Sự suy giảm và chuyển trạng thái từ đất có rừng sang đất trống với độ tàn che giảm xuống dưới 0,1 và diện tích liên vùng suy giảm dưới 0,3 ha.
  • Chỉ số khác biệt thực vật chuẩn hóa (NDVI): Tỷ số giữa dải sóng NIR và Red, dao động từ -1,0 đến +1,0, trong đó thảm phủ rừng nhiệt đới thường đạt giá trị cao từ 0,60 đến 0,85.
  • Chỉ số hỏa hoạn và độ ẩm chuẩn hóa (NBR): Tỷ số giữa dải sóng NIR và SWIR, phản ánh sinh khối ẩm và mức độ tổn thương sau cháy rừng hoặc khai thác.
  • Chỉ số viễn thám tích hợp (IRSI): Chỉ số tổng hợp được tính bằng tổng của NDVI và NBR nhằm tối ưu hóa khả năng tách lọc biến động thảm phủ.

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

Nghiên cứu sử dụng nguồn dữ liệu ảnh vệ tinh Sentinel-2 đa phổ với độ phân giải không gian 10m và chu kỳ chụp lặp lại 5 ngày, được truy xuất và xử lý trực tiếp trên nền tảng điện toán đám mây Google Earth Engine (GEE). Dữ liệu đối chứng bao gồm bản đồ kiểm kê rừng năm 2016, cơ sở dữ liệu diễn biến rừng FRMS tỉnh Thừa Thiên Huế đến năm 2019 và số liệu quan trắc khí tượng lượng mưa giai đoạn 2002 - 2019.

Cỡ mẫu nghiên cứu gồm 250 điểm giải đoán thực địa và mẫu khóa ảnh được phân bố đại diện trên toàn bộ các huyện miền núi và đồng bằng. Phương pháp chọn mẫu phân tầng ngẫu nhiên (Stratified Random Sampling) được áp dụng dựa trên các trạng thái: rừng tự nhiên giàu, trung bình, nghèo, rừng trồng keo/tràm và đất trống sau khai thác.

Lý do lựa chọn phương pháp phân tích sai khác ngưỡng chỉ số đa thời gian (Image Differencing) kết hợp đánh giá chéo (Cross-validation) và kiểm chuẩn chất lượng ảnh NIIRS là nhằm tự động hóa quá trình xử lý ảnh trên diện rộng, loại bỏ sai số do mây mù và tiết kiệm trên 80% thời gian khảo sát thực tế. Toàn bộ quá trình đánh giá được chia thành 3 kỳ quan sát: Kỳ 1 (tháng 4/2017 đến tháng 4/2018), Kỳ 2 (tháng 4/2018 đến tháng 4/2019) và Kỳ 3 (tháng 4/2019 đến tháng 5/2020), tập trung vào các tháng mùa khô từ tháng 2 đến tháng 6 khi tỷ lệ mây bình quân dưới 20%.

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

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

Quá trình phân tích thực nghiệm trên không gian địa bàn tỉnh Thừa Thiên Huế đã mang lại các phát hiện khoa học quan trọng:

Thứ nhất, nghiên cứu đã lượng hóa chính xác ngưỡng giá trị quang phổ đặc trưng cho hiện tượng mất rừng. Ngưỡng sụt giảm của chỉ số NDVI (dNDVI) được xác định ở mức lớn hơn 0,32 (tương đương mức suy giảm trên 42% giá trị sinh khối xanh). Ngưỡng sụt giảm của chỉ số NBR (dNBR) đạt mức lớn hơn 0,28 (tương đương mức suy giảm trên 38% độ ẩm và cấu trúc tán). Ngưỡng biến động của chỉ số tích hợp IRSI được thiết lập ở mức dIRSI lớn hơn 0,60.

Thứ hai, kết quả đánh giá chéo độ chính xác phát hiện điểm mất rừng cho thấy sự vượt trội rõ rệt của chỉ số IRSI. Độ chính xác tổng thể của chỉ số tích hợp IRSI đạt 92,5%, cao hơn 5,2% so với khi sử dụng đơn lẻ chỉ số NDVI (đạt 87,3%) và cao hơn 6,9% so với chỉ số NBR (đạt 85,6%).

Thứ ba, tổng diện tích mất rừng được bóc tách trong cả giai đoạn 2017 - 2020 biến động trên 24.500 ha (trung bình hơn 8.000 ha mỗi năm). Trong đó, diện tích mất rừng tập trung trên 65% ở đối tượng rừng trồng sản xuất đến kỳ thu hoạch luân canh tại các huyện Nam Đông, Phong Điền và Phú Lộc; khoảng 35% diện tích còn lại phân bố rải rác do cháy rừng, sạt lở tự nhiên và chuyển đổi mục đích tại các khu vực đồi núi hiểm trở.

Thảo luận kết quả

Khả năng nhận diện vượt trội của chỉ số IRSI bắt nguồn từ cơ chế kết hợp đồng thời độ phản xạ diệp lục ở bước sóng ánh sáng nhìn thấy và độ nhạy cảm cấu trúc ẩm ở bước sóng ngắn SWIR. Khi một diện tích rừng bị triệt hạ, tán cây xanh biến mất làm giảm mạnh NDVI, đồng thời bề mặt đất trống bị phơi lộ làm gia tăng mạnh bức xạ nhiệt và phản xạ dải sóng SWIR, khiến giá trị NBR rơi tự do. Việc kết hợp hai chỉ số này giúp triệt tiêu hiện tượng báo động giả do hiện tượng rụng lá theo mùa hoặc thay đổi độ ẩm mặt đất cục bộ.

So sánh với các công bố khoa học tương đồng, kết quả này hoàn toàn tương thích với mô hình phân tích BFAST của Michael Schultz (2016) trên dữ liệu nhiệt đới và báo cáo thực nghiệm năm 2020 của Bộ Nông nghiệp và Phát triển nông thôn tại Điện Biên (dNDVI lớn hơn 0,308 và dNBR lớn hơn 0,248). Sự hiệu chỉnh ngưỡng dNDVI lên mức 0,32 và dNBR lên 0,28 phản ánh chính xác đặc thù sinh khối rừng mưa nhiệt đới thường xanh của vùng duyên hải Bắc Trung Bộ.

Dữ liệu nghiên cứu được biểu diễn trực quan qua biểu đồ biến thiên 12 tháng của năm 2019, chỉ ra giá trị NDVI đạt đỉnh trung bình 0,78 vào tháng 5 và sụt giảm tự nhiên xuống 0,55 vào mùa mưa lũ tháng 10 và 11. Bản đồ phân bố không gian thể hiện rõ các điểm nóng mất rừng tập trung tại lưu vực sông Hương, sông Ô Lâu và khu vực rừng giáp ranh huyện A Lưới, tạo tiền đề vững chắc cho việc trực quan hóa dữ liệu trên các bảng biểu thống kê chuyên ngành.

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

Dựa trên các kết luận thực nghiệm, nghiên cứu đề xuất 4 nhóm giải pháp chiến lược nhằm hiện đại hóa công tác quản lý bảo vệ rừng:

  • Phát triển và tích hợp module phát hiện mất rừng tự động: Xây dựng công cụ phần mềm chuyên dụng vận hành trên nền tảng Google Earth Engine kết hợp phần mềm mã nguồn mở QGIS cho Chi cục Kiểm lâm Thừa Thiên Huế. Mục tiêu là tự động quét ảnh Sentinel-2 định kỳ và giảm thời gian phát hiện điểm nghi vấn mất rừng từ 30 ngày xuống dưới 5 ngày, hoàn thành triển khai trong giai đoạn 2021 - 2022 do Chi cục Kiểm lâm phối hợp với Sở Nông nghiệp và Phát triển nông thôn chủ trì.
  • Chuẩn hóa quy trình ứng dụng chỉ số tích hợp IRSI: Ban hành hướng dẫn kỹ thuật thống nhất việc sử dụng ngưỡng dIRSI lớn hơn 0,60 trong công tác rà soát biến động rừng định kỳ hàng tháng. Phấn đấu nâng tỷ lệ xác thực chính xác các điểm biến động ngoài thực địa đạt trên 95%, áp dụng ngay từ quý I/2021 bởi Hạt Kiểm lâm 8 huyện, thị xã và thành phố.
  • Nâng cao năng lực chuyên môn cho lực lượng kiểm lâm cơ sở: Tổ chức các chương trình tập huấn chuyên sâu về xử lý ảnh viễn thám, sử dụng thiết bị GPS cầm tay và ứng dụng phần mềm FRMS Mobile cho 100% cán bộ kiểm lâm địa bàn (với hơn 150 biên chế toàn tỉnh). Lộ trình thực hiện trong vòng 12 tháng dưới sự tài trợ của các chương trình phát triển lâm nghiệp bền vững.
  • Hoàn thiện hạ tầng cơ sở dữ liệu không gian rừng: Đầu tư hệ thống máy chủ và liên kết chia sẻ dữ liệu ảnh vệ tinh độ phân giải 10m giữa các Ban Quản lý rừng phòng hộ, Vườn Quốc gia Bạch Mã và các Khu bảo tồn thiên nhiên. Mục tiêu hướng tới giảm thiểu 30% diện tích rừng tự nhiên bị xâm hại trái phép đến năm 2025 theo chỉ đạo của Ủy ban nhân dân tỉnh Thừa Thiên Huế.

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

Công trình nghiên cứu mang tính ứng dụng thực tiễn cao và là tài liệu tham khảo giá trị cho 4 nhóm đối tượng chính:

  • Cán bộ quản lý và lực lượng kiểm lâm địa phương: Nắm vững quy trình xử lý ảnh tự động và thuật toán chỉ số IRSI để khoanh vùng chính xác điểm mất rừng, giúp tiết kiệm 60% công sức tuần tra phân tán trên tổng diện tích 288.334,37 ha rừng của tỉnh.
  • Các nhà nghiên cứu và học viên cao học ngành Quản lý tài nguyên, Lâm nghiệp và GIS: Khai thác khung phương pháp luận xử lý dữ liệu lớn đa thời gian trên nền tảng điện toán đám mây GEE, phương pháp xác định ngưỡng biến động và kỹ thuật đánh giá chéo mẫu ảnh viễn thám độ phân giải 10m.
  • Các tổ chức phi chính phủ và cơ quan tài trợ bảo tồn quốc tế: Sử dụng cơ sở dữ liệu biến động 3 kỳ (2017 - 2020) và bản đồ phân bố mất rừng để hoạch định các chương trình tài trợ phục hồi hệ sinh thái, giảm phát thải khí nhà kính REDD+ tại hành lang đa dạng sinh học dãy Trường Sơn.
  • Doanh nghiệp lâm nghiệp và chủ rừng trồng quy mô lớn: Ứng dụng công nghệ giám sát viễn thám để quản lý chu kỳ khai thác và trồng lại rừng, hỗ trợ duy trị chứng chỉ quản lý rừng bền vững FSC trên hơn 7.000 ha rừng sản xuất.

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

Ảnh vệ tinh Sentinel-2 có những ưu thế gì nổi bật trong giám sát tài nguyên rừng tại Thừa Thiên Huế? Sentinel-2 cung cấp nguồn tư liệu đa phổ hoàn toàn miễn phí với độ phân giải không gian cao 10m và chu kỳ chụp lặp lại 5 ngày. Tần suất quan sát dày đặc này giúp các nhà nghiên cứu dễ dàng chọn lọc các cảnh ảnh không mây trong các tháng mùa khô (tháng 2 đến tháng 6), đảm bảo độ tin cậy phát hiện biến động rừng đạt trên 90%.

Chỉ số tích hợp IRSI vượt trội hơn các chỉ số truyền thống như NDVI và NBR ở điểm nào? Chỉ số IRSI là sự kết hợp cộng gộp giữa NDVI (nhạy với sắc tố diệp lục tươi) và NBR (nhạy với độ ẩm và sinh khối khô). Sự kết hợp này triệt tiêu các tín hiệu nhiễu do thay đổi độ ẩm đất hoặc rụng lá sinh lý, nâng độ chính xác nhận diện mất rừng lên 92,5%, cao hơn chỉ số NDVI đơn lẻ 5,2%.

Làm thế nào để loại bỏ sai số do mây mù che phủ khi phân tích viễn thám tại miền Trung? Nghiên cứu ứng dụng thuật toán lọc mây tự động trên Google Earth Engine kết hợp thang đo chất lượng ảnh chuẩn quốc tế NIIRS. Đồng thời, tác giả tập trung xử lý dữ liệu vào khung thời gian có lượng mưa bình quân thấp nhất trong năm (dưới 100 mm/tháng), giúp tỷ lệ mây bình quân trên các cảnh ảnh thu nhận luôn ở mức dưới 20%.

Luận văn định nghĩa tiêu chí mất rừng như thế nào để phục vụ bóc tách dữ liệu vệ tinh? Theo Luật Lâm nghiệp 2017 và FAO, mất rừng trong nghiên cứu được xác định là sự chuyển đổi từ trạng thái đất có rừng sang đất trống với diện tích suy giảm từ 0,3 ha trở lên và độ tàn che tụt xuống dưới 0,1. Về mặt quang phổ, hiện tượng này tương ứng với mức giảm dNDVI lớn hơn 0,32.

Quy trình phát hiện mất rừng tự động được ứng dụng vào thực tế kiểm lâm theo các bước nào? Quy trình gồm 4 bước liên hoàn: tải và lọc ảnh Sentinel-2 trên GEE, tính toán chỉ số chênh lệch delta đa thời kỳ, lọc điểm mất rừng theo ngưỡng biến động tối ưu và giao cắt dữ liệu với bản đồ quy hoạch 3 loại rừng để loại bỏ các cảnh báo nhầm ngoài lâm phận. Toàn bộ chu trình chỉ mất dưới 24 giờ xử lý dữ liệu.

Kết luận

  • Luận văn đã xây dựng thành công quy trình ứng dụng ảnh vệ tinh Sentinel-2 độ phân giải 10m trên nền tảng Google Earth Engine để giám sát biến động tài nguyên rừng định kỳ 5 ngày.
  • Nghiên cứu chứng minh chỉ số tích hợp IRSI đạt độ chính xác cao nhất (92,5%), vượt trội hoàn toàn so với việc sử dụng riêng rẽ NDVI (87,3%) và NBR (85,6%).
  • Bóc tách và định lượng chính xác diện tích mất rừng trên 288.334,37 ha đất lâm nghiệp tỉnh Thừa Thiên Huế qua 3 kỳ giám sát liên tục từ năm 2017 đến năm 2020.
  • Cung cấp luận cứ khoa học thực tiễn giúp ngành lâm nghiệp địa phương cắt giảm 60% thời gian và chi phí tuần tra thực địa nhờ tự động hóa cảnh báo sớm.
  • Đề xuất lộ trình tích hợp thuật toán viễn thám vào phần mềm FRMS giai đoạn 2021 - 2025, tạo bước đột phá trong chuyển đổi số quản lý tài nguyên và bảo tồn đa dạng sinh học.

Công trình là tài liệu tham khảo thiết thực cho các nhà quản lý lâm nghiệp, cơ quan bảo tồn và giới nghiên cứu. Hãy tải toàn văn luận văn thạc sĩ để khai thác chi tiết thuật toán xử lý ảnh, bộ mã nguồn phân tích viễn thám và dữ liệu không gian thực nghiệm giá trị.