Tổng quan nghiên cứu

Rừng ngập mặn tại Việt Nam đóng vai trò là vành đai xanh phòng hộ xung yếu trước biến đổi khí hậu, trải dài dọc theo bờ biển hơn 3000 km với hơn 50 loài thực vật đặc trưng. Tuy nhiên, tình trạng suy thoái tài nguyên rừng và áp lực chuyển đổi mục đích sử dụng đất ven biển đang đặt ra thách thức cấp bách trong công tác giám sát định lượng carbon sinh khối. Các phương pháp viễn thám quang học truyền thống thường bị hạn chế nghiêm trọng bởi tầng mây dày đặc và chỉ phân biệt được ranh giới diện tích bề mặt mà không phản ánh chính xác cấu trúc ba chiều của lâm phần.

Nghiên cứu tập trung giải quyết bài toán định lượng sinh khối trên mặt đất thông qua công nghệ viễn thám RADAR độ mở tổng hợp chủ động, khắc phục triệt để rào cản thời tiết và mây mù. Mục tiêu then chốt của luận văn là nghiên cứu đặc tính tán xạ, phân cực của sóng viễn thám băng C từ vệ tinh ENVISAT ASAR nhằm xây dựng mô hình toán học ước tính sinh khối khô cho thảm rừng ngập mặn ven biển hai tỉnh Thái Bình và Nam Định thuộc vùng Đồng bằng sông Hồng.

Phạm vi không gian bao quát dải ven biển châu thổ sông Hồng với chu kỳ dữ liệu ảnh vệ tinh cập nhật 35 ngày ở độ cao quỹ đạo khoảng 800 km, đối chiếu trực tiếp cùng số liệu đo đạc mặt đất. Kết quả nghiên cứu mang lại ý nghĩa then chốt cho tiến trình kiểm kê rừng quốc gia, lượng hóa bể hấp thụ khí nhà kính và cung cấp tham số phục vụ cơ chế chi trả dịch vụ môi trường rừng cũng như các chương trình giảm phát thải khí nhà kính toàn cầu.

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, lý thuyết tương tác sóng điện từ viễn thám chủ động và mô hình ước lượng sinh khối phi tham số. Sinh khối rừng trên mặt đất được định nghĩa là toàn bộ chất hữu cơ sống của thảm thực vật bao gồm thân, cành, lá trên một đơn vị diện tích, tính theo đơn vị tấn/ha trọng lượng khô. Cấu trúc không gian lâm phần được xác định qua đường kính thân cây tại vị trí 1,3 mét, chiều cao vút ngọn và mật độ cá thể.

Hệ thống lý thuyết viễn thám RADAR sóng cực ngắn băng C với bước sóng 5,7 cm phản xạ chủ yếu tại tầng tán lá và cành nhỏ của thảm phủ thực vật. Cơ chế tán xạ ngược là kết quả tổng hợp của tán xạ bề mặt, tán xạ khối và đa tán xạ thể tích giữa tán cây và nền bùn ngập nước. Cường độ tán xạ ngược phụ thuộc vào hằng số điện môi của thực vật và góc quét quan sát dao động từ 14 độ đến 45 độ. Các nghiên cứu quốc tế chỉ ra rằng tín hiệu viễn thám RADAR nhạy cảm với dải sinh khối rừng từ 50 tấn/ha đến 200 tấn/ha trước khi đạt ngưỡng bão hòa khoảng 250 tấn/ha.

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

Nguồn dữ liệu đầu vào tích hợp đa nguồn bao gồm ảnh viễn thám ENVISAT ASAR phân giải không gian 30 mét, dữ liệu quang học SPOT và hệ thống tọa độ định vị toàn cầu GPS thu thập trong giai đoạn thực địa năm 2010 đến năm 2011. Phương pháp chọn mẫu ngẫu nhiên phân tầng được triển khai với 30 ô tiêu chuẩn kích thước 100 mét vuông tại các kiểu rừng thuần loài bần chua và trang, kèm theo 3 ô mẫu độc lập phục vụ công tác kiểm chứng mô hình.

Quy trình xử lý dữ liệu bao gồm hiệu chỉnh bức xạ, lọc nhiễu đốm Speckle bằng thuật toán thích nghi, chuyển đổi hệ số tán xạ ngược sang đơn vị decibel và nắn chỉnh hình học theo lưới chiếu bản đồ phẳng UTM. Phương pháp phân tích thống kê hồi quy đơn biến và đa biến được lựa chọn nhằm thiết lập tương quan phi tuyến giữa hệ số tán xạ ngược các kênh phân cực HH, VV và tỷ số phân cực HH/VV với trữ lượng sinh khối khô thực địa. Phương pháp này tối ưu hóa độ chính xác nhờ việc lượng hóa trực tiếp mối quan hệ tán xạ thể tích của cấu trúc tán cây ven biển.

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

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

Phân tích định lượng dữ liệu phản xạ cho thấy hệ số tán xạ ngược của kênh phân cực ngang HH và tỷ số hai kênh phân cực HH/VV có mối tương quan rất chặt chẽ với sinh khối khô của thảm thực vật ngập mặn. Cụ thể, khi sinh khối rừng tăng từ mức thấp 16 tấn/ha ở các dải cây tái sinh chiều cao dưới 2 mét lên mức 85 tấn/ha ở các quần xã cây trưởng thành cao trên 5 mét, giá trị tán xạ ngược kênh HH tăng tuyến tính trong khoảng từ âm 14 decibel đến âm 7 decibel.

Mô hình hồi quy phi tuyến dạng hàm mũ thực nghiệm giữa sinh khối khô và tỷ số phân cực HH/VV đạt hệ số xác định R bình phương ở mức 0,78, cao hơn 25% so với việc chỉ sử dụng kênh phân cực đơn lẻ VV. Nghiên cứu đã xác lập bản đồ phân bố trữ lượng sinh khối tại khu vực cửa sông Ba Lạt và các bãi bồi ven biển Kim Sơn, Giao Thủy với mức độ phân hóa rõ rệt: nhóm sinh khối thấp dưới 20 tấn/ha chiếm 42% diện tích thảm ngập mặn mới bồi tụ, nhóm sinh khối trung bình từ 20 đến 60 tấn/ha chiếm 38% và nhóm sinh khối cao trên 60 tấn/ha chiếm 20% tập trung tại các khu bảo tồn được quản lý nghiêm ngặt.

Thảo luận kết quả

Mối liên hệ tỷ lệ thuận giữa hệ số tán xạ ngược và sinh khối khô bắt nguồn từ sự gia tăng của diện tích bề mặt tán lá và mạng lưới cành cây, làm tăng mạnh cơ chế tán xạ khối khi bước sóng băng C 5,7 cm tương tác trực tiếp với tầng vòm cây. Tại các vùng có độ ẩm nền cao hoặc ngập triều, tỷ số phân cực HH/VV đã loại bỏ hiệu quả ảnh hưởng phản xạ gương từ mặt nước, giúp giá trị hồi quy phản ánh trung thực sinh khối thực vật thay vì độ ẩm nền đất.

Dữ liệu nghiên cứu được biểu diễn trực quan thông qua đồ thị tán xạ tương quan hình chuông và hệ thống bản đồ phân lớp raster với thang màu phân cấp mật độ trữ lượng 5 bậc từ 0 đến 100 tấn/ha. Khi so sánh với kết quả kiểm chứng tại 3 ô mẫu độc lập ngoài thực địa, sai số tuyệt đối trung bình RMSE ghi nhận dưới 12,5%, khẳng định tính khả thi vượt trội so với các mô hình viễn thám quang học đơn thuần. Hiện tượng bão hòa tín hiệu tán xạ chỉ bắt đầu xuất hiện cục bộ khi mật độ sinh khối vượt ngưỡng 80 tấn/ha, phù hợp với các quy luật vật lý đã được kiểm chứng tại vùng rừng ngập mặn nhiệt đới Đông Nam Á.

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

Thứ nhất, Bộ Nông nghiệp và Phát triển nông thôn cùng Chi cục Kiểm lâm các tỉnh ven biển cần ban hành quy chuẩn kỹ thuật ứng dụng ảnh vệ tinh RADAR băng C và băng L vào quy trình giám sát tài nguyên rừng thường niên từ năm 2025. Việc này hướng tới mục tiêu tự động hóa 100% bản đồ biến động sinh khối và giảm 60% chi phí điều tra đo đếm thủ công trên các vùng đầm lầy hiểm trở.

Thứ hai, Sở Nông nghiệp và Phát triển nông thôn hai tỉnh Thái Bình và Nam Định cần thành lập trạm quan trắc viễn thám chuyên dụng, tích hợp thuật toán tỷ số phân cực đa kênh để cập nhật dữ liệu sinh khối định kỳ 3 tháng một lần. Mục tiêu then chốt là kiểm soát chặt chẽ 12.000 hecta rừng ngập mặn hiện có và phát hiện sớm các vụ phá rừng trái phép phục vụ nuôi trồng thủy hải sản.

Thứ ba, Ban quản lý các Vườn quốc gia và Khu bảo tồn thiên nhiên Đất ngập nước cần triển khai số hóa hệ thống 50 ô tiêu chuẩn định vị GPS cố định từ quý 3 năm 2024 để cung cấp chuỗi số liệu đối soát mặt đất liên tục cho các thế hệ vệ tinh quan sát Trái Đất mới.

Thứ tư, các viện nghiên cứu và cơ sở đào tạo lâm nghiệp cần mở rộng việc tích hợp dữ liệu RADAR giao thoa và thuật toán học máy nhằm nâng cao giới hạn phát hiện sinh khối vượt ngưỡng bão hòa 150 tấn/ha trong giai đoạn 2025 đến 2030, phục vụ việc xây dựng hạn ngạch tín chỉ carbon rừng đạt tiêu chuẩn quốc tế.

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

Luận văn là tài liệu chuyên khảo giá trị dành cho các nhà hoạch định chính sách và cơ quan quản lý nhà nước về lâm nghiệp, giúp xây dựng phương án định giá dịch vụ hệ sinh thái và quản lý hơn 200.000 hecta rừng ngập mặn trên phạm vi toàn quốc dựa trên nền tảng dữ liệu số chính xác.

Nghiên cứu là cẩm nang hữu ích cho các chuyên gia phân tích hệ thống thông tin địa lý và viễn thám, cung cấp bộ quy trình tiền xử lý ảnh RADAR, hiệu chỉnh góc tới và kỹ thuật lọc nhiễu sóng với độ chuẩn xác cao, giúp rút ngắn 40% thời gian xử lý dữ liệu viễn thám chuyên sâu.

Các tổ chức phát triển quốc tế và doanh nghiệp tham gia thị trường carbon có thể sử dụng phương pháp tính toán trong luận văn để xây dựng hệ thống đo đạc, báo cáo và thẩm định lượng phát thải phục vụ phát hành tín chỉ carbon rừng theo tiêu chuẩn toàn cầu.

Giảng viên, nghiên cứu sinh và học viên cao học thuộc các ngành Địa lý, Lâm nghiệp, Quản lý tài nguyên môi trường và Công nghệ không gian sẽ tìm thấy khung lý thuyết bài bản, nguồn tài liệu tham khảo phong phú cùng các công thức hồi quy thực nghiệm ứng dụng trực tiếp trong giảng dạy và nghiên cứu mở rộng.

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

Viễn thám RADAR có ưu điểm gì vượt trội so với ảnh quang học trong đánh giá sinh khối rừng ngập mặn? Công nghệ RADAR sử dụng sóng siêu cao tần chủ động có khả năng xuyên qua tầng mây dày, sương mù và hoạt động liên tục cả ngày lẫn đêm. Băng sóng C 5,7 cm tương tác trực tiếp với thể tích tán lá và cấu trúc cành cây, mang lại thông tin cấu trúc ba chiều thay vì chỉ ghi nhận màu sắc bề mặt như ảnh quang học.

Tại sao tỷ số phân cực HH/VV lại cho kết quả ước tính sinh khối chính xác hơn phân cực đơn kênh? Kênh phân cực đơn lẻ thường bị nhiễu bởi độ gồ ghề bề mặt bùn lầy và mực nước thủy triều lên xuống. Tỷ số phân cực HH/VV giúp triệt tiêu phần lớn ảnh hưởng phản xạ từ mặt đất ẩm ướt, làm nổi bật tín hiệu tán xạ thể tích thuần túy từ tán rừng, giúp tăng hệ số tương quan thêm khoảng 25%.

Ngưỡng bão hòa của sóng viễn thám RADAR băng C đối với sinh khối rừng ngập mặn là bao nhiêu? Theo các tính toán thực nghiệm, sóng RADAR băng C bắt đầu bão hòa tín hiệu khi mật độ sinh khối khô chạm mức 80 đến 100 tấn/ha. Đối với các khu rừng nguyên sinh có trữ lượng vượt quá ngưỡng 200 tấn/ha, nghiên cứu khuyến nghị kết hợp thêm dữ liệu sóng băng L hoặc băng P có bước sóng dài hơn.

Quy mô diện tích tối thiểu của ô mẫu thực địa để khớp với ảnh ENVISAT ASAR là bao nhiêu? Ảnh vệ tinh ENVISAT ASAR ở chế độ Image Mode có độ phân giải không gian là 30 mét x 30 mét, tương đương kích thước một điểm ảnh khoảng 900 mét vuông. Do đó, các ô mẫu thực địa cần đạt diện tích tối thiểu từ 100 mét vuông đến 400 mét vuông tại các vị trí thảm rừng đồng nhất để hạn chế sai số vị trí điểm ảnh.

Dữ liệu của luận văn có thể áp dụng trực tiếp cho các vùng rừng ngập mặn khác tại Việt Nam không? Mô hình toán học hoàn toàn có thể chuyển giao cho các dải rừng ngập mặn ven biển Bắc Bộ và Bắc Trung Bộ có cấu trúc tương đồng. Tuy nhiên, khi áp dụng cho các khu rừng đước đôi phát triển cao lớn tại Cà Mau với sinh khối vượt trên 150 tấn/ha, các tham số hồi quy cần được hiệu chỉnh lại theo mẫu đo đạc thực địa tại chỗ.

Kết luận

Luận văn đã giải quyết xuất sắc nhiệm vụ khoa học và thực tiễn trong việc định lượng cấu trúc rừng ngập mặn ven biển Đồng bằng sông Hồng bằng công nghệ viễn thám vi sóng tiên tiến:

  • Xác lập hoàn chỉnh cơ sở vật lý tương tác sóng viễn thám băng C trên quần xã thực vật ngập mặn đặc thù miền Bắc.
  • Chứng minh tính ưu việt của chỉ số tỷ số phân cực HH/VV trong việc nâng cao độ chính xác mô hình hóa sinh khối rừng.
  • Xây dựng thành công bản đồ phân bố trữ lượng sinh khối khô cho vùng bờ biển Thái Bình và Nam Định với sai số dưới 12,5%.
  • Đề xuất quy trình công nghệ chuẩn hóa từ tiền xử lý ảnh viễn thám đến thiết lập hàm tương quan phi tuyến phục vụ ngành lâm nghiệp.
  • Cung cấp luận cứ khoa học vững chắc phục vụ chương trình kiểm kê carbon và chi trả dịch vụ môi trường rừng quốc gia.

Trong giai đoạn 6 tháng tới, các đơn vị quản lý cần số hóa toàn diện quy trình xử lý dữ liệu và tiếp tục tích hợp thêm dữ liệu vệ tinh đa bước sóng để mở rộng phạm vi ứng dụng trên quy mô toàn quốc. Hãy chủ động liên hệ các chuyên gia viễn thám địa lý và triển khai giải pháp công nghệ không gian ngay hôm nay nhằm tối ưu hóa công tác bảo tồn tài nguyên thiên nhiên và phát triển bền vững dải ven biển Việt Nam.