Tổng quan nghiên cứu

Tỉnh Điện Biên sở hữu tổng diện tích tự nhiên là 956.290,37 ha (tương đương 9.562,90 km2), chiếm 16,96% diện tích vùng Tây Bắc và 2,89% diện tích cả nước. Đây là địa bàn xung yếu về quốc phòng - an ninh và đầu nguồn của ba hệ thống sông lớn gồm sông Đà, sông Mã và sông Mê Kông. Địa hình toàn tỉnh bị chia cắt hiểm trở với hơn 70% quỹ đất có độ dốc trên 25 độ, khiến tài nguyên rừng giữ vai trò sống còn trong việc phòng hộ đầu nguồn, điều tiết nguồn nước và chống xói mòn rửa trôi. Tuy nhiên, áp lực gia tăng dân số với quy mô 527.000 người (mật độ bình quân 55,1 người/km2) cùng tập quán canh tác nương rẫy đã khiến thảm phủ rừng biến động mạnh mẽ, phần lớn diện tích chỉ còn là rừng non tái sinh với trữ lượng bình quân khiêm tốn từ 80 đến 245 m3/ha.

Mục tiêu trọng tâm của nghiên cứu là đánh giá định lượng mức độ và xác định các nguyên nhân gây biến động lớp phủ rừng tỉnh Điện Biên trong giai đoạn 12 năm (2002 - 2014). Đề tài ứng dụng công nghệ viễn thám hiện đại với phương pháp phân loại hướng đối tượng (OBIA - Object-Based Image Analysis) kết hợp hệ thống thông tin địa lý (GIS) nhằm thành lập hệ thống bản đồ hiện trạng và bản đồ biến động rừng ở tỷ lệ 1:100.000 theo hệ tọa độ chuẩn quốc gia VN-2000.

Ý nghĩa khoa học và thực tiễn của công trình thể hiện ở việc nâng cao độ chính xác trong giải đoán ảnh viễn thám đa thời gian, khắc phục hoàn toàn hiện tượng lẫn phổ của các phương pháp phân loại điểm ảnh truyền thống. Kết quả nghiên cứu cung cấp cơ sở dữ liệu số hóa tin cậy, hỗ trợ đắc lực cho các nhà hoạch định chính sách địa phương trong công tác kiểm kê, quy hoạch và phát triển lâm nghiệp 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 nền tảng lý luận về hệ thái sinh rừng và phân loại thảm thực vật nhiệt đới. Cơ sở pháp lý chính dựa trên Luật Bảo vệ và Phát triển rừng năm 2004 cùng Thông tư số 34/2009/TT-BNNPTNT, trong đó xác định ba tiêu chí chuẩn hóa về rừng: độ tàn che của tán cây từ 0,1 trở lên, diện tích liền khoảnh tối thiểu từ 0,5 ha và chiều cao vút ngọn của cây gỗ từ 5,0 mét trở lên.

Về mặt học thuật, đề tài kế thừa hệ thống phân loại 14 kiểu thảm thực vật sinh thái phát sinh của Thái Văn Trừng và mô hình 4 cấp độ trạng thái rừng theo hiện trạng của Loeschau (từ đất trống đồi trọc loại I đến rừng giàu phục hồi hoàn toàn loại IV). Đồng thời, nghiên cứu vận dụng quy luật phản xạ phổ của diệp lục thực vật ở dải sóng cận hồng ngoại (NIR) và vùng ánh sáng nhìn thấy (Red) thông qua chỉ số khác biệt thực vật NDVI (Normalized Difference Vegetation Index) để phân tách trạng thái sinh trưởng của lớp phủ. Khái niệm biến động rừng được tiếp cận toàn diện theo quan điểm sinh thái cảnh quan, đặt mối quan hệ tương tác giữa các quá trình diễn thế tự nhiên và hoạt động kinh tế - xã hội của con người trong một không gian mở.

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

Nguồn dữ liệu sơ cấp phục vụ đề tài bao gồm 5 cảnh ảnh vệ tinh Landsat độ phân giải không gian 30 mét: 2 cảnh ảnh Landsat 7 ETM+ chụp vào tháng 03/2002 và 3 cảnh ảnh Landsat 8 OLI/TIRS chụp vào tháng 03/2014. Toàn bộ ảnh được lựa chọn trong cùng mùa khô để đồng nhất điều kiện góc chiếu mặt trời và giảm thiểu tối đa độ che phủ của mây. Dữ liệu thứ cấp gồm bản đồ nền địa lý tỉnh Điện Biên tỷ lệ 1:100.000 (hệ tọa độ VN-2000), số liệu niên giám thống kê và báo cáo ngành lâm nghiệp địa phương.

Về phương pháp chọn mẫu, nghiên cứu áp dụng kỹ thuật lấy mẫu thực địa phân tầng ngẫu nhiên với cỡ mẫu tối thiểu 10 ô tiêu chuẩn định vị và bán định vị cho mỗi trạng thái rừng mục tiêu, kết hợp phỏng vấn sâu người dân địa phương để kiểm chứng lịch sử canh tác.

Lý do lựa chọn phương pháp phân loại hướng đối tượng (OBIA) trên phần mềm eCognition thay vì phương pháp dựa trên pixel truyền thống là nhằm khai thác đồng thời cả thông tin phổ, hình thái, cấu trúc và ngữ cảnh không gian. Thuật toán phân vùng đa phân giải (multi-resolution segmentation) tự động gộp các pixel liền kề có độ đồng nhất cao thành các đa giác đối tượng (polygons), hạn chế tối đa hiện tượng nhiễu "muối tiêu". Sau đó, thuật toán láng giềng gần nhất (Standard Nearest Neighbor) kết hợp khóa giải đoán NDVI được áp dụng để phân loại chi tiết các lớp phủ. Cuối cùng, công cụ phân tích không gian GIS với các phép toán chồng xếp bản đồ (raster and vector overlay) và ma trận chuyển đổi diện tích được sử dụng để bóc tách các vùng biến động trong suốt chu kỳ 2002 - 2014.

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

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

Quá trình xử lý ảnh vệ tinh và phân tích không gian GIS đã mang lại các phát hiện cốt lõi về thực trạng và động thái rừng Điện Biên:

  • Phân bố diện tích theo lập địa: Cơ cấu đất lâm nghiệp có rừng năm 2013 đạt 637.525,80 ha, chiếm tới 81,83% tổng diện tích đất nông nghiệp toàn tỉnh. Tài nguyên rừng tập trung chủ yếu trên hai nhóm đất then chốt: đất đỏ vàng với 485.559,62 ha (chiếm 50,78% diện tích tự nhiên) và đất mùn vàng đỏ trên núi cao với 424.086,44 ha (chiếm 44,35%). Nhóm đất phù sa ven sông suối chỉ chiếm 35.222,43 ha (3,68%) nhưng chịu áp lực chuyển đổi nông nghiệp mạnh nhất.
  • Biến động cơ cấu ba loại rừng: Trong tổng thể diện tích đất có rừng, rừng phòng hộ chiếm tỷ trọng lớn nhất với 305.421 ha (khoảng 47,9%), tiếp theo là rừng sản xuất với 213.211 ha (khoảng 33,4%) và rừng đặc dụng chiếm 118.893 ha (khoảng 18,7%). Giai đoạn 2002 - 2014 ghi nhận sự phân hóa rõ rệt giữa diện tích rừng tự nhiên suy giảm cục bộ và diện tích rừng trồng mới gia tăng tại các dải đồi thấp.
  • Thực trạng chất lượng thảm thực vật: Rừng tự nhiên nguyên sinh giàu gỗ quý (như lát hoa, pơ mu, nghiến, táu) chỉ còn tồn tại rải rác ở các đỉnh núi hiểm trở thuộc huyện Mường Nhé và khu vực biên giới. Phần lớn diện tích là rừng non tái sinh nghèo kiệt với trữ lượng thấp từ 80 đến 245 m3/ha. Quỹ đất chưa sử dụng còn rất lớn với 148.562,11 ha (chiếm 15,54% diện tích tự nhiên toàn tỉnh), chủ yếu là trảng cỏ và cây bụi trơ sỏi đá.
  • Độ chính xác phân loại vượt trội: Kết quả đánh giá ma trận sai số trên ảnh Landsat 8 năm 2014 bằng phương pháp OBIA đạt hệ số tin cậy Kappa vượt 0,85 và độ chính xác tổng thể trên 88%, cao hơn đáng kể so với mức xấp xỉ 73% của các phương pháp phân loại theo điểm ảnh trước đây.

Thảo luận kết quả

Nguyên nhân chính dẫn đến biến động rừng tại Điện Biên là sự tổng hòa của các yếu tố nhân sinh và điều kiện tự nhiên khắc nghiệt. Sự gia tăng dân số cùng tỷ lệ hộ nghèo cao tại các huyện vùng sâu như Mường Nhé, Tuần Giáo, Mường Chà dẫn đến nhu cầu mở rộng đất nương rẫy và khai thác gỗ củi làm chất đốt. Thêm vào đó, việc phát triển các công trình hạ tầng giao thông và xây dựng các bậc thang thủy điện dọc theo lưu vực sông Đà (diện tích lưu vực nội tỉnh hơn 5.000 km2), sông Mã và sông Nậm Rốm đã chuyển đổi trực tiếp nhiều diện tích thảm thực vật ven sông suối thành đất công trình và lòng hồ ngập nước.

Về phương thức biểu diễn trực quan, các dữ liệu biến động không gian của nghiên cứu được mô hình hóa tối ưu thông qua ma trận chuyển dịch hai chiều (ma trận biến động lớp phủ), cho phép đối chiếu chính xác diện tích chuyển hóa qua lại giữa từng cặp trạng thái (như từ rừng giàu sang rừng nghèo, từ rừng tự nhiên sang đất nương rẫy, hoặc từ đất trống sang rừng phục hồi). Kết hợp với đó là hệ thống biểu đồ cột trực quan so sánh diện tích lớp phủ năm 2002 và năm 2014 cùng bản đồ biến động số hóa phân tầng không gian.

So sánh với các nghiên cứu quốc tế của M. Hansen (2013) hay nghiên cứu tại Việt Nam của Nguyễn Trường Sơn (2009), kết quả này hoàn toàn tương đồng về xu thế: các khu vực rừng có độ dốc cao và xa nguồn nước thường bảo tồn tốt hơn, trong khi vùng đệm tiếp giáp sông suối trong phạm vi bán kính 200 đến 500 mét có tỷ lệ suy thoái cao gấp 3 lần do khả năng tiếp cận khai thác của con người quá dễ dàng.

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

Dựa trên các phân tích không gian và nguyên nhân thực tế, luận văn đưa ra 4 nhóm giải pháp cụ thể:

  1. Thiết lập hệ thống giám sát rừng tự động bằng công nghệ viễn thám và GIS: Chi cục Kiểm lâm tỉnh Điện Biên cần chủ trì triển khai quy trình giải đoán ảnh vệ tinh đa thời gian định kỳ 6 tháng/lần, tích hợp dữ liệu GPS tuần tra thực địa để cập nhật biến động lớp phủ rừng. Mục tiêu phát hiện sớm và xử lý 100% các vụ xâm hại, phá rừng trái phép ngay trong quý đầu tiên phát sinh.
  2. Khoanh nuôi xúc tiến tái sinh trên nhóm đất mùn vàng đỏ núi cao: Sở Nông nghiệp và Phát triển nông thôn cần tập trung nguồn vốn quy hoạch bảo vệ nghiêm ngặt 424.086,44 ha đất rừng trên độ cao trên 900 mét. Đặt mục tiêu nâng độ che phủ rừng toàn tỉnh thêm từ 3% đến 5% trong giai đoạn 2015 - 2025, đặc biệt ưu tiên phục hồi các dải rừng phòng hộ xung yếu đầu nguồn sông Đà và sông Mã.
  3. Mở rộng chính sách chi trả dịch vụ môi trường rừng (PES): Ủy ban nhân dân tỉnh Điện Biên cần hoàn thiện cơ chế giao khoán bảo vệ 305.421 ha rừng phòng hộ và 118.893 ha rừng đặc dụng cho các cộng đồng thôn bản. Mục tiêu nâng mức thu nhập từ nghề rừng lên 20% đến 30% tổng thu nhập hộ gia đình vào năm 2025, tạo động lực kinh tế bền vững để đồng bào dân tộc thiểu số từ bỏ tập quán phá rừng làm nương.
  4. Chuyển đổi mô hình nông lâm kết hợp trên đất dốc: Chính quyền các huyện Mường Nhé, Tuần Giáo, Điện Biên Đông cần chỉ đạo chuyển đổi canh tác nương rẫy du canh trên 148.562,11 ha đất chưa sử dụng sang mô hình trồng cây ăn quả, cây công nghiệp dài ngày kết hợp cây lâm nghiệp bản địa (như trẩu, keo, cao su). Biện pháp này giúp cố định đất, chống xói mòn trên các sườn đồi dốc trên 25 độ trước năm 2030.

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

Nội dung và phương pháp luận của luận văn mang lại giá trị thực tiễn cho 4 nhóm đối tượng chính:

  • Cơ quan quản lý nhà nước về lâm nghiệp và Chi cục Kiểm lâm: Hưởng lợi từ quy trình chuẩn hóa bản đồ hiện trạng tỷ lệ 1:100.000 và cơ sở dữ liệu số hóa VN-2000, phục vụ trực tiếp công tác kiểm kê tài nguyên rừng định kỳ 5 năm và lập quy hoạch 3 loại rừng cấp tỉnh.
  • Chuyên gia và kỹ sư công nghệ địa không gian (GIS/Remote Sensing): Tiếp cận quy trình chi tiết về kỹ thuật phân đoạn đa phân giải (multi-resolution segmentation) và xây dựng bộ khóa giải đoán NDVI trên phần mềm eCognition, làm tài liệu tham khảo để xử lý ảnh vệ tinh quang học độ phân giải trung bình và cao.
  • Giảng viên, nghiên cứu sinh và học viên cao học chuyên ngành Quản lý Tài nguyên - Môi trường: Khai thác khung lý thuyết hệ sinh thái rừng, phương pháp chọn mẫu ô tiêu chuẩn thực địa và kỹ thuật phân tích ma trận biến động để triển khai các đề tài nghiên cứu liên quan đến suy thoái đất và thảm phủ.
  • Các tổ chức phi chính phủ và ban quản lý dự án bảo tồn thiên nhiên: Sử dụng số liệu kiểm kê 740 loài thực vật thân gỗ và 59 loài động vật quý hiếm làm căn cứ xây dựng hồ sơ bảo tồn đa dạng sinh học, kêu gọi tài trợ cho các dự án chi trả dịch vụ carbon rừng (REDD+) tại khu vực Tây Bắc.

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

1. Phương pháp phân loại hướng đối tượng (OBIA) có ưu điểm gì vượt trội so với phân loại điểm ảnh truyền thống?
Phân loại OBIA không phân tích riêng lẻ từng pixel mà gộp các pixel đồng nhất thành các đối tượng không gian (polygons) dựa trên thuật toán đa phân giải. Phương pháp này tận dụng đồng thời giá trị phản xạ phổ, hình thái, cấu trúc và ngữ cảnh lân cận, giúp loại bỏ hoàn toàn hiện tượng nhiễu "muối tiêu" và nâng hệ số chính xác Kappa từ 0,73 lên trên 0,85.

2. Nguồn dữ liệu ảnh vệ tinh Landsat được ứng dụng trong nghiên cứu gồm những loại nào?
Nghiên cứu sử dụng 5 cảnh ảnh vệ tinh Landsat độ phân giải không gian 30 mét, bao gồm 2 cảnh ảnh Landsat 7 ETM+ chụp vào tháng 03/2002 và 3 cảnh ảnh Landsat 8 OLI/TIRS chụp vào tháng 03/2014. Tất cả các cảnh ảnh đều được chụp vào thời điểm mùa khô ít mây, đảm bảo tính đồng nhất cao khi đối chiếu chuỗi thời gian 12 năm.

3. Đặc điểm địa hình và thổ nhưỡng nào chi phối mạnh nhất đến phân bố rừng Điện Biên?
Điện Biên có hơn 70% diện tích đất tự nhiên sở hữu độ dốc trên 25 độ với địa hình núi cao chia cắt mạnh. Rừng phân bố chủ yếu trên nhóm đất đỏ vàng chiếm 485.559,62 ha (50,78%) và đất mùn vàng đỏ núi cao chiếm 424.086,44 ha (44,35%), tạo nên các vành đai thực vật phân tầng rõ nét theo đai cao địa hình.

4. Luận văn đã áp dụng những tiêu chuẩn pháp lý nào để định danh và phân loại trạng thái rừng?
Nghiên cứu tuân thủ nghiêm ngặt Luật Bảo vệ và Phát triển rừng năm 2004 cùng Thông tư số 34/2009/TT-BNNPTNT của Bộ Nông nghiệp và Phát triển nông thôn. Các tiêu chí pháp lý bao gồm: độ tàn che của tán rừng đạt từ 0,1 trở lên, diện tích khoảnh rừng tập trung từ 0,5 ha trở lên và chiều cao vút ngọn cây gỗ đạt tối thiểu 5,0 mét.

5. Ma trận biến động lớp phủ trong GIS giúp ích gì cho công tác quản lý tài nguyên địa phương?
Ma trận biến động là bảng thống kê hai chiều đối chiếu diện tích của từng trạng thái lớp phủ giữa năm 2002 và năm 2014. Công cụ này cho phép nhà quản lý nhìn thấy rõ bản chất của sự thay đổi (diện tích rừng giàu chuyển sang nghèo, rừng chuyển sang đất nương rẫy hoặc đất trống phục hồi thành rừng), từ đó định vị chính xác điểm nóng cần can thiệp.

Kết luận

  • Đề tài đã lượng hóa chi tiết diễn biến tài nguyên rừng trên toàn bộ 956.290,37 ha diện tích tự nhiên tỉnh Điện Biên xuyên suốt chu kỳ 12 năm (2002 - 2014).
  • Ứng dụng thành công phương pháp phân loại hướng đối tượng (OBIA) trên ảnh Landsat 7 và Landsat 8, đạt độ chính xác Kappa trên 0,85 và tối ưu hóa quy trình chiết tách thông tin lớp phủ.
  • Xác lập bức tranh phân bố thực trạng của 637.525,80 ha đất có rừng, chứng minh vai trò trọng yếu của nhóm đất đỏ vàng và đất mùn núi cao đối với chức năng phòng hộ đầu nguồn.
  • Chỉ rõ các nguyên nhân nhân sinh chính làm suy giảm trữ lượng rừng xuống mức 80 đến 245 m3/ha, bao gồm áp lực dân số, du canh nương rẫy và phát triển thủy điện.
  • Đề xuất đồng bộ 4 nhóm giải pháp gắn kết giữa công nghệ giám sát viễn thám hiện đại với chính sách chi trả dịch vụ môi trường rừng và chuyển đổi sinh kế bền vững.

Đóng góp lớn nhất của luận văn là đã chuẩn hóa một quy trình công nghệ kết hợp giữa OBIA và GIS có chi phí thấp, độ chính xác cao và khả năng chuyển giao thực tế cho ngành lâm nghiệp vùng cao. Kế hoạch tiếp theo trong giai đoạn 2015 - 2025 là tích hợp các nguồn ảnh vệ tinh độ phân giải siêu cao (như Sentinel, SPOT) để tự động hóa hoàn toàn hệ thống cảnh báo mất rừng theo thời gian thực. Các cơ quan quản lý và nhà nghiên cứu quan tâm có thể khai thác bộ dữ liệu số hóa và bản đồ chuyên đề tỷ lệ 1:100.000 để ứng dụng trực tiếp vào công tác quy hoạch lâm nghiệp tại địa phương.