Tổng quan nghiên cứu

Tỉnh Vĩnh Phúc có tổng diện tích tự nhiên đạt 123.054,4 ha với cơ cấu hành chính gồm 9 huyện, thành phố và dân số đạt hơn 1.054,5 ngàn người vào năm 2015. Theo Quyết định số 667/QĐ-CT của Ủy ban nhân dân tỉnh Vĩnh Phúc về việc phê duyệt kết quả kiểm kê rừng, toàn tỉnh sở hữu 33.928,7 ha đất lâm nghiệp, chiếm 27,5% diện tích tự nhiên. Trong cơ cấu đất lâm nghiệp, diện tích đất rừng phòng hộ chiếm 4.170,8 ha, giữ vai trò đặc biệt xung yếu đối với việc bảo vệ nguồn nước, điều hòa khí hậu và chống xói mòn đất cho các lưu vực sông lớn như sông Hồng (đoạn qua tỉnh dài 28 km), sông Lô (dài 31 km) và sông Phó Đáy (dài 36 km).

Vấn đề nghiên cứu trọng tâm xuất phát từ áp lực gia tăng dân số và tốc độ công nghiệp hóa nhanh chóng của tỉnh, khiến cảnh quan và lớp phủ rừng phòng hộ đầu nguồn chịu nhiều tác động phức tạp. Trước đây, việc giám sát tài nguyên rừng địa phương chủ yếu dựa trên các phương pháp đo vẽ thủ công truyền thống, tiêu tốn nhiều thời gian, công sức và thiếu tính đồng bộ theo chuỗi thời gian. Mục tiêu cụ thể của luận văn là ứng dụng công nghệ viễn thám kết hợp hệ thống thông tin địa lý (GIS) nhằm định lượng và xây dựng bản đồ đánh giá biến động diện tích rừng phòng hộ tỉnh Vĩnh Phúc trong suốt 27 năm (1990 - 2017) qua 7 mốc thời gian chính: 1990, 1995, 2000, 2005, 2010, 2015 và 2017.

Ý nghĩa của nghiên cứu thể hiện qua việc cung cấp hệ thống số liệu khoa học xác thực, làm cơ sở vững chắc cho các nhà hoạch định chính sách tối ưu hóa giải pháp quản lý lưu vực bền vững, nâng cao chất lượng môi trường sinh thái cho hơn 1.247 loài thực vật bậc cao và 1.188 loài động vật rừng sinh sống trên địa bàn.

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 trên nền tảng lý thuyết sinh thái học rừng kinh điển của Morozov và Tcachenco, khẳng định rừng là một thể thống nhất hữu cơ giữa thực vật, động vật, vi sinh vật và môi trường lập địa. Hệ thống phân loại rừng của luận văn tuân thủ nghiêm ngặt Thông tư số 34/2009/TT-BNNPTNT của Bộ Nông nghiệp và Phát triển nông thôn, quy định tiêu chuẩn xác định rừng với độ tàn che của tán cây thân gỗ từ 0,1 trở lên và diện tích liền khoảnh tối thiểu 0,5 ha. Đồng thời, nghiên cứu kết hợp quan điểm phân loại trạng thái rừng của Loeschau với 4 cấp độ diễn thế và quan điểm sinh thái phát sinh quần thể thực vật gồm 14 kiểu rừng của Thái Văn Trừng.

Khung lý thuyết viễn thám tập trung vào đặc tính phản xạ phổ của các đối tượng tự nhiên trong dải sóng nhìn thấy và cận hồng ngoại. Trọng tâm là việc áp dụng chỉ số thực vật chuẩn hóa NDVI (Normalized Difference Vegetation Index) được tính theo công thức tỷ lệ sai khác giữa kênh cận hồng ngoại và kênh đỏ. Chỉ số NDVI dao động trong dải từ -1 đến +1, trong đó các giá trị nhỏ hơn 0,1 tương ứng với vùng mặt nước hoặc đất trống, giá trị từ 0,299 đến 0,469 phản ánh rừng non chưa có trữ lượng, và giá trị từ 0,469 đến 1,0 tương ứng với trạng thái rừng tự nhiên phục hồi hoặc rừng giàu. Ba khái niệm nền tảng được vận hành xuyên suốt gồm: phân loại ảnh hướng đối tượng (OBIA), ma trận chuyển dịch không gian và quản lý bền vững lưu vực đầu nguồn.

       +-------------------------------------------------------------+
       |   Hệ thống dữ liệu ảnh vệ tinh đa thời gian (1990 - 2017)   |
       |   - Landsat 5 TM (1990, 1995, 2005, 2010)                   |
       |   - Landsat 7 ETM+ (2000), Sentinel-2A (2015), Landsat 8 (2017) |
       +------------------------------+------------------------------+
                                      |
                                      v
       +-------------------------------------------------------------+
       | Tiền xử lý & Chuẩn hóa dữ liệu (Hệ tọa độ không gian WGS84) |
       +------------------------------+------------------------------+
                                      |
                                      v
       +-------------------------------------------------------------+
       |   Phân loại hướng đối tượng (OBIA) trên phần mềm eCognition  |
       |   - Phân đoạn đa độ phân giải (Multiresolution Segmentation)|
       |   - Cây phân loại tự động kết hợp chỉ số phổ (NDVI, GLCM)   |
       +------------------------------+------------------------------+
                                      |
                                      v
       +-------------------------------------------------------------+
       |    Kiểm chứng thực địa bằng GPS & Đánh giá ma trận biến động |
       |    trên ArcGIS (4 trạng thái: RTN, RTG, DT, DKH)            |
       +-------------------------------------------------------------+

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

Nghiên cứu sử dụng nguồn dữ liệu viễn thám đa thời gian gồm 7 cảnh ảnh vệ tinh chất lượng cao được thu thập từ Cục Khảo sát Địa chất Hoa Kỳ (USGS) và Cơ quan Không gian Châu Âu (ESA). Cụ thể bao gồm ảnh Landsat 5 TM các năm 1990, 1995, 2005, 2010; Landsat 7 ETM+ năm 2000; Sentinel-2A năm 2015 (độ phân giải 20 m) và Landsat 8 OLI năm 2017 (độ phân giải 30 m). Toàn bộ dữ liệu bản đồ vector và raster được quy chuẩn thống nhất về hệ quy chiếu tọa độ quốc tế WGS84 bằng phần mềm ERDAS IMAGINE và ENVI.

Về phương pháp chọn mẫu, nghiên cứu thiết lập tối thiểu 20 điểm mẫu thực địa cho mỗi trạng thái rừng trên từng cảnh ảnh thông qua 3 đến 5 tuyến điều tra ngoại nghiệp bằng máy định vị vệ tinh GPS. Các điểm mẫu kiểm tra được bố trí nghiêm ngặt, cách ranh giới chuyển tiếp giữa các trạng thái tối thiểu 50 m nhằm loại bỏ sai số biên. Phương pháp phân loại hướng đối tượng (Object-Based Image Analysis - OBIA) trên phần mềm eCognition Developer 8.7 được lựa chọn thay thế cho phương pháp phân loại theo điểm ảnh (pixel) truyền thống. Lý do lựa chọn OBIA là khả năng tổng hợp đồng thời thông tin phổ, chỉ số cấu trúc bề mặt (GLCM Homogeneity, GLCM Contrast) và mối quan hệ không gian giữa các đối tượng, giúp triệt tiêu hiện tượng nhiễu hạt "muối tiêu" và tăng cường độ chính xác khi phân tách 4 lớp phủ chính: Rừng tự nhiên (RTN), Rừng trồng (RTG), Đất lâm nghiệp chưa có rừng (DT) và Đất khác (DKH).

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

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

Quá trình giải đoán ảnh vệ tinh đa thời gian đã xây dựng thành công 7 bộ bản đồ hiện trạng rừng phòng hộ và 6 bản đồ ma trận biến động không gian tại tỉnh Vĩnh Phúc. Độ chính xác tổng thể của phương pháp giải đoán ảnh qua các thời kỳ đều đạt trên 85%, hệ số Kappa đạt trên 0,82, khẳng định tính chuẩn xác và độ tin cậy cao của quy trình phân loại tự động.

Giai đoạn biến động Xu hướng dịch chuyển lớp phủ chính Tác động diện tích tương đối Đánh giá độ chính xác phân loại
1990 - 1995 Rừng tự nhiên suy giảm sang đất trống Giảm khoảng 3,5% - 5,2% Độ chính xác đạt 85,3%
1995 - 2000 Rừng phục hồi bước đầu, đất trống thu hẹp Tăng nhẹ diện tích che phủ Độ chính xác đạt 86,1%
2000 - 2005 Rừng trồng mở rộng mạnh mẽ Tăng bình quân > 15% diện tích Độ chính xác đạt 87,4%
2005 - 2010 Ổn định cấu trúc rừng tự nhiên và rừng trồng Biến động diện tích cục bộ < 2% Độ chính xác đạt 88,0%
2010 - 2015 Đất lâm nghiệp chưa có rừng giảm mạnh Rừng trồng đạt bình quân 747,6 ha/năm Độ chính xác đạt 89,2%
2015 - 2017 Rừng phòng hộ duy trì ổn định bền vững Tổng diện tích đạt 4.170,8 ha Độ chính xác đạt 90,5%

Các phát hiện chính của luận văn bao gồm:

  • Giai đoạn 1990 - 1995: Diện tích đất trống chưa có rừng chiếm tỷ trọng lớn trong quy hoạch rừng phòng hộ do tình trạng khai thác gỗ củi tự do và tập quán du canh nương rẫy, khiến diện tích rừng tự nhiên sụt giảm cục bộ từ 3,5% đến 5,2%.
  • Giai đoạn 1995 - 2005: Chứng kiến sự đảo chiều mạnh mẽ khi diện tích rừng trồng phòng hộ tăng trưởng vượt bậc với tỷ lệ tăng bình quân trên 15% mỗi chu kỳ 5 năm, thu hẹp đáng kể các diện tích đất đồi núi trọc.
  • Giai đoạn 2005 - 2017: Diện tích đất có rừng tiếp tục được củng cố ổn định, đưa tổng diện tích đất rừng phòng hộ toàn tỉnh đạt mức 4.170,8 ha. Diện tích trồng rừng tập trung toàn tỉnh trong giai đoạn 2011 - 2015 đạt mức bình quân 747,6 ha/năm, song song với 149 ha/năm diện tích cây trồng phân tán quy đổi.
Diễn biến tỷ lệ lớp phủ rừng phòng hộ Vĩnh Phúc (1990 - 2017)
100% |============================================================|
 80% | [Đất trống: ~45%] -------------> [Đất trống: <15%]         |
 60% |                                                            |
 40% | [Rừng trồng: ~10%] ------------> [Rừng trồng: >40%]        |
 20% | [Rừng tự nhiên: ~45%] ---------> [Rừng tự nhiên: ~45%]     |
  0% +------------------------------------------------------------+
     Năm 1990                                             Năm 2017

Thảo luận kết quả

Nguyên nhân căn bản của xu hướng gia tăng diện tích rừng phòng hộ bắt nguồn từ việc thực thi hiệu quả các chương trình mục tiêu quốc gia về lâm nghiệp (Chương trình 327 và Dự án 661) kết hợp với chính sách giao đất khoán rừng theo Quyết định số 667/QĐ-CT của Ủy ban nhân dân tỉnh. Cơ chế này đã tạo động lực kinh tế trực tiếp cho người dân vùng đệm tham gia trồng rừng và bảo vệ lưu vực.

Khi so sánh với các công trình nghiên cứu của Viện Điều tra, Quy hoạch rừng (FIPI) và kết quả nghiên cứu lưu vực của một số tác giả trong nước cùng thời kỳ, xu thế biến động tại Vĩnh Phúc hoàn toàn tương thích với quy luật phục hồi rừng chung của vùng Đông Bắc Bộ. Việc sử dụng ma trận chuyển dịch không gian (cross-tabulation matrix) trong nghiên cứu đã làm rõ bản chất của việc chuyển hóa hàng ngàn hecta đất trống đồi trọc thành các thảm rừng trồng có giá trị phòng hộ cao, góp phần giảm thiểu lưu lượng đỉnh lũ tại lưu vực sông Phó Đáy (vốn có lưu lượng mùa lũ lên tới 833 m³/s) và bảo vệ an toàn cho các hồ thủy lợi lớn như hồ Thanh Lanh, hồ Xạ Hương và hồ Đại Lải.

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

Dựa trên kết quả phân tích chuỗi biến động không gian 27 năm, luận văn đề xuất 4 nhóm giải pháp quản lý bền vững lưu vực rừng phòng hộ tỉnh Vĩnh Phúc:

  • Giải pháp 1 (Kỹ thuật lâm sinh): Tiến hành khoanh nuôi xúc tiến tái sinh tự nhiên kết hợp trồng bổ sung các loài cây bản địa có bộ rễ sâu, tán rộng như kim giao, vù hương, lim xanh và re gừng trên 4.170,8 ha đất rừng phòng hộ. Mục tiêu cụ thể là nâng cao độ tàn che của toàn bộ diện tích rừng phòng hộ lên mức trên 0,7 trong giai đoạn 2020 - 2030, do Chi cục Kiểm lâm tỉnh Vĩnh Phúc chủ trì thực hiện.
  • Giải pháp 2 (Ứng dụng công nghệ giám sát số): Xây dựng hệ thống cơ sở dữ liệu GIS tích hợp ảnh viễn thám Sentinel-2 và Landsat 8 với chu kỳ cập nhật định kỳ từ 15 đến 30 ngày/lần. Hệ thống nhằm tự động phát hiện sớm các điểm mất rừng, cháy rừng và sạt lở lưu vực, hoàn thành triển khai trước năm 2025 dưới sự phối hợp của Sở Nông nghiệp và Phát triển nông thôn cùng Sở Tài nguyên và Môi trường.
  • Giải pháp 3 (Chính sách an sinh và chi trả dịch vụ môi trường): Hoàn thiện chính sách chi trả dịch vụ môi trường rừng (DVMTR) từ nguồn thu của các cơ sở sản xuất nước sinh hoạt, du lịch sinh thái và thủy lợi. Nâng mức hỗ trợ khoán bảo vệ rừng lên khoảng 400.000 - 500.000 đồng/ha/năm nhằm tạo sinh kế ổn định cho người dân vùng đệm, đặc biệt là đồng bào dân tộc thiểu số như dân tộc Sán Dìu (chiếm 3,93% dân số toàn tỉnh), do Ủy ban nhân dân tỉnh Vĩnh Phúc phê duyệt.
  • Giải pháp 4 (Quản lý lưu vực và phòng cháy chữa cháy rừng): Thiết kế và tu bổ hệ thống đường băng cản lửa rộng từ 10 đến 15 m kết hợp trạm quan trắc thủy văn tự động tại các khu vực xung yếu thuộc dãy Tam Đảo và vùng thượng nguồn sông Phó Đáy, đảm bảo giảm thiểu trên 90% nguy cơ cháy rừng trong mùa khô hanh từ tháng 10 đến tháng 3 năm sau.

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

Nội dung luận văn mang lại giá trị thực tiễn và học thuật sâu sắc cho 4 nhóm đối tượng trọng tâm sau:

  • Cơ quan quản lý nhà nước và Chi cục Kiểm lâm: Cung cấp bộ công cụ bản đồ số và cơ sở dữ liệu biến động rừng 27 năm làm căn cứ trực tiếp để rà soát, điều chỉnh quy hoạch 3 loại rừng và xây dựng phương án kiểm kê rừng định kỳ 5 năm cấp tỉnh.
  • Nhà nghiên cứu, giảng viên và học viên chuyên ngành Lâm nghiệp, Quản lý tài nguyên, GIS: Cung cấp tài liệu tham khảo chuẩn mực về phương pháp phân loại ảnh hướng đối tượng (OBIA) trên eCognition và kỹ thuật xử lý ảnh viễn thám đa thời gian (Landsat, Sentinel) trong nghiên cứu tài nguyên thiên nhiên.
  • Ban quản lý rừng phòng hộ, rừng đặc dụng và Vườn quốc gia Tam Đảo: Ứng dụng quy trình giám sát lớp phủ thực vật để kiểm soát ranh giới đất rừng, bảo tồn 83 loài thực vật nguy cấp và 64 loài động vật quý hiếm ghi trong Sách Đỏ Việt Nam.
  • Các tổ chức phát triển bền vững và doanh nghiệp tham gia thị trường tín chỉ carbon (REDD+): Sử dụng chuỗi số liệu viễn thám liên tục làm đường phát thải cơ sở (FREL/FRL), phục vụ việc đo đạc, báo cáo và thẩm định (MRV) trữ lượng sinh khối hấp thụ carbon của các thảm rừng đầu nguồn.

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

1. Nghiên cứu sử dụng những nguồn tư liệu ảnh viễn thám nào để đánh giá biến động rừng phòng hộ? Nghiên cứu tích hợp chuỗi 7 cảnh ảnh vệ tinh quang học đa thời gian gồm Landsat 5 TM (các năm 1990, 1995, 2005, 2010), Landsat 7 ETM+ (năm 2000), Sentinel-2A (năm 2015, độ phân giải không gian 20 m) và Landsat 8 OLI (năm 2017, độ phân giải 30 m) trên nền hệ tọa độ WGS84.

2. Chỉ số thực vật NDVI đóng vai trò gì trong quá trình phân loại trạng thái rừng? Chỉ số NDVI được tính toán từ kênh phổ đỏ và cận hồng ngoại nhằm định lượng mật độ sinh khối và độ tươi tốt của thảm thực vật. Với các ngưỡng phân loại cụ thể (NDVI < 0,1 là đất trống/mặt nước, 0,299 - 0,469 là rừng non, và > 0,469 là rừng tự nhiên/rừng giàu), NDVI giúp thuật toán phân tách chính xác các kiểu lớp phủ.

3. Tại sao phương pháp phân loại hướng đối tượng (OBIA) bằng eCognition lại được ưu tiên sử dụng? Khác với phương pháp phân loại theo từng pixel dễ gây hiện tượng "muối tiêu", OBIA trên eCognition gom cụm các điểm ảnh liền kề có cùng đặc tính phổ thành các đối tượng (objects). Phương pháp này kết hợp thông tin phổ, cấu trúc bề mặt (GLCM) và hình học, nâng độ chính xác giải đoán lên trên 85%.

4. Yếu tố nào đóng vai trò then chốt giúp gia tăng diện tích rừng phòng hộ tại Vĩnh Phúc qua 27 năm? Sự gia tăng diện tích rừng trồng phòng hộ bắt nguồn từ việc thực hiện đồng bộ các chương trình trồng rừng quốc gia kết hợp chính sách giao đất khoán rừng theo Quyết định số 667/QĐ-CT của Ủy ban nhân dân tỉnh, giúp duy trì diện tích trồng rừng tập trung đạt bình quân 747,6 ha/năm.

5. Kết quả nghiên cứu có ý nghĩa thực tiễn như thế nào đối với công tác bảo vệ lưu vực sông tại địa phương? Bản đồ biến động và ma trận chuyển dịch giúp xác định chính xác các điểm nóng suy thoái đất tại vùng đầu nguồn sông Phó Đáy và dãy Tam Đảo. Điều này hỗ trợ các cấp quản lý xây dựng giải pháp công trình và phi công trình nhằm giảm thiểu nguy cơ lũ quét, sạt lở và xói mòn đất cho các hồ chứa lớn.

Kết luận

  • Luận văn đã ứng dụng thành công công nghệ viễn thám đa thời gian và hệ thống thông tin địa lý (GIS) để số hóa toàn diện diễn biến rừng phòng hộ tỉnh Vĩnh Phúc trong chuỗi thời gian 27 năm (1990 - 2017).
  • Quy trình phân loại hướng đối tượng (OBIA) trên phần mềm eCognition đạt độ chính xác cao trên 85% và hệ số Kappa trên 0,82, chuẩn hóa bộ khóa phân loại 4 trạng thái lớp phủ chính.
  • Phân tích định lượng khẳng định sự phục hồi tích cực của diện tích rừng trồng phòng hộ, đưa tổng diện tích đất rừng phòng hộ đạt mức 4.170,8 ha, đóng góp trực tiếp vào an ninh nguồn nước và ổn định môi trường lưu vực.
  • Đề tài đóng góp cơ sở lý luận và thực tiễn vững chắc thông qua 4 nhóm giải pháp kỹ thuật lâm sinh, công nghệ giám sát số, cơ chế chính sách dịch vụ môi trường và quản lý lưu vực xung yếu đến năm 2030.
  • Các nhà quản lý, nghiên cứu sinh và chủ rừng quan tâm đến giải pháp quản lý lưu vực bền vững có thể khai thác và ứng dụng ngay quy trình công nghệ này để nâng cao hiệu quả bảo vệ tài nguyên rừng tại địa phương.