Tổng quan nghiên cứu

Việt Nam sở hữu đường bờ biển dài hơn 3.000 km cùng hệ sinh thái rừng ngập mặn phong phú với hơn 50 loài thực vật đặc trưng. Trong bối cảnh biến đổi khí hậu toàn cầu ngày càng gay gắt, thảm rừng ngập mặn tại các cửa sông ven biển đóng vai trò đặc biệt quan trọng như một lá chắn phòng hộ tự nhiên và là bể hấp thụ carbon khổng lồ. Vườn quốc gia Xuân Thủy thuộc tỉnh Nam Định, nằm ở cửa Ba Lạt của hạ lưu sông Hồng, được UNESCO công nhận là khu Ramsar thứ 50 của thế giới từ tháng 01/1989 và được Thủ tướng Chính phủ ký Quyết định số 01/2003/QĐ-TTg ngày 02/01/2003 nâng hạng thành vườn quốc gia.

Tuy nhiên, công tác kiểm kê và giám sát tài nguyên rừng ngập mặn hiện nay chủ yếu dựa vào ảnh viễn thám quang học và đo đạc thủ công ngoài thực địa. Dữ liệu quang học thụ động thường bị giới hạn bởi mây mù và chỉ phản ánh thông tin bề mặt hai chiều, không cung cấp được trữ lượng sinh khối ba chiều của thảm cây. Để giải quyết rào cản này, đề tài nghiên cứu tập trung ứng dụng ảnh vệ tinh viễn thám Radar khẩu độ tổng hợp băng C từ bộ cảm ENVISAT ASAR kết hợp với ảnh quang học SPOT năm 2010 và dữ liệu điều tra thực địa nhằm xác định chính xác sinh khối rừng ngập mặn ven biển tỉnh Thái Bình và Nam Định.

Nghiên cứu được triển khai trên phạm vi dải rừng ngập mặn ven biển châu thổ sông Hồng trong giai đoạn năm 2010 đến năm 2011. Kết quả nghiên cứu mang ý nghĩa thực tiễn to lớn khi định lượng chính xác các cấp trữ lượng sinh khối khô dưới 150 tấn/ha, giúp tiết kiệm hơn 70% thời gian và chi phí khảo sát truyền thống, đồng thời cung cấp cơ sở dữ liệu nền tảng cho việc chi trả dịch vụ môi trường rừng và tham gia thị trường tín chỉ carbon quốc tế.

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 tương tác sóng điện từ vi sóng Radar khẩu độ tổng hợp (SAR) và các quy luật cấu trúc sinh thái học rừng ngập mặn.

  • Lý thuyết tán xạ sóng Radar chủ động: Bộ cảm SAR sử dụng bước sóng vi ba băng C (khoảng 5,6 cm đến 5,7 cm) có khả năng xuyên qua mây mù, mưa gió và thấu qua các khoảng hở của tầng tán thực vật. Tín hiệu thu nhận là hệ số tán xạ ngược biểu thị bằng decibel (dB), phản ánh mức độ tán xạ khối từ thân cành, tán lá và tán xạ góc đôi giữa thân cây với bề mặt bùn ngập triều.
  • Lý thuyết phân cực vi sóng: Sóng Radar phát và thu ở các dạng phân cực ngang và phân cực đứng (HH, VV, HV, VH). Tỷ số giữa các kênh phân cực (HH/VV) có khả năng triệt tiêu các tác động bất lợi do độ ẩm nền và độ gồ ghề của mặt đất bùn lầy, làm nổi bật phản xạ thuần túy từ sinh khối cây rừng.
  • Cấu trúc sinh thái rừng ngập mặn: Quần xã thực vật rừng được phân tích qua cấu trúc thẳng đứng (tầng tán chính, tầng dưới tán), cấu trúc nằm ngang (mật độ cây, đường kính thân D1.3, chiều cao vút ngọn Hvn) và cấu trúc tuổi theo từng giai đoạn phát triển 5 năm.
  • Mô hình quan hệ allometric thực nghiệm: Thiết lập mối quan hệ hàm số giữa các chỉ tiêu đo đếm hình thái cây cá thể với tổng sinh khối tươi và sinh khối khô của lâm phần.

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

  • Nguồn dữ liệu sử dụng: Dữ liệu ảnh vệ tinh viễn thám Radar ENVISAT ASAR năm 2010 chế độ phân cực luân phiên băng C (độ phân giải 30 m) và ảnh vệ tinh quang học đa phổ SPOT năm 2010 (độ phân giải 10 m).
  • Cỡ mẫu và phương pháp chọn mẫu: Đề tài tiến hành thu thập số liệu trên 13 ô tiêu chuẩn thực địa điển hình với kích thước mỗi ô là 100 m2 (10 m x 10 m) phân bố tại các khu vực rừng Vẹt và Trang thuộc Vườn quốc gia Xuân Thủy. Phương pháp chọn mẫu phân tầng ngẫu nhiên có đại diện được áp dụng dựa trên mật độ và chiều cao cây. Tập dữ liệu gồm 10 ô tiêu chuẩn được dùng để thiết lập hàm tương quan hồi quy và 3 ô tiêu chuẩn độc lập (ô số 5, 14, 15) được sử dụng để kiểm chứng độ chính xác mô hình.
  • Phương pháp phân tích và lý do lựa chọn:
    1. Phương pháp phân loại ảnh quang học có kiểm định: Ứng dụng trên ảnh SPOT để bóc tách chính xác đường ranh giới lớp phủ rừng ngập mặn, làm mặt nạ không gian cắt dữ liệu ảnh Radar.
    2. Tiền xử lý ảnh SAR chuyên sâu: Thực hiện chuyển đổi biên độ sang năng lượng, định chuẩn bức xạ, nắn chỉnh hình học tọa độ khớp với ảnh SPOT và lọc nhiễu hạt bằng thuật toán lọc Frost. Bộ lọc Frost được lựa chọn vì ứng dụng mô hình nhiễu thích ứng cục bộ, giúp loại bỏ tối đa nhiễu đóm hạt mà vẫn bảo tồn nguyên vẹn các đường ranh giới và chi tiết cấu trúc rừng.
    3. Phương pháp phân tích hồi quy thống kê: Xây dựng phương trình tương quan giữa hệ số tán xạ ngược (dB) tại các phân cực HH, VV, tỷ số HH/VV với sinh khối khô thực địa (tấn/ha). Việc lựa chọn hồi quy tỷ số phân cực HH/VV giúp triệt tiêu biến động của mực nước thủy triều tại vùng ven biển, tạo ra mô hình ước tính sinh khối có độ tin cậy cao nhất.
  • Timeline nghiên cứu: Toàn bộ quá trình thu thập số liệu thực địa, tiền xử lý ảnh và tính toán mô hình được hoàn thành trong giai đoạn từ tháng 03/2011 đến tháng 12/2011.

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

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

Nghiên cứu đã đưa ra các kết quả định lượng cụ thể về khả năng ứng dụng các kênh phân cực Radar trong việc ước tính sinh khối rừng ngập mặn:

  • Tương quan đơn phân cực VV đạt mức rất thấp: Mối quan hệ giữa hệ số tán xạ ngược phân cực VV và sinh khối khô chỉ đạt hệ số xác định R2 = 0,32. Giá trị này quá thấp, không đủ tiêu chuẩn khoa học để xây dựng mô hình toán học tính toán sinh khối cho khu vực.
  • Tương quan phân cực HH có sự cải thiện nhưng sai số phân loại lớn: Phương trình hồi quy phân cực HH đạt tương quan tuyến tính khá tốt, tuy nhiên khi áp dụng vào phân loại không gian, diện tích nhóm sinh khối thấp bị phóng đại quá mức. Cụ thể, diện tích nhóm sinh khối từ 0 đến 50 tấn/ha bị ước tính lên tới 956,96 ha.
  • Tỷ số hai phân cực HH/VV cho độ tin cậy vượt trội: Tỷ số phân cực HH/VV thể hiện mối tương quan chặt chẽ nhất với sinh khối khô thực địa theo phương trình hồi quy Y = -387,4 * (HH/VV). Mô hình này đạt độ chính xác kiểm chứng 100% trên cả 3 ô mẫu kiểm tra độc lập (ô số 5, 14 và 15).
  • Phân bố diện tích theo các cấp sinh khối thực tế: Khi sử dụng ảnh tỷ số HH/VV, toàn bộ diện tích rừng ngập mặn được phân định thành 3 cấp rõ rệt: Cấp 1 (0 đến 50 tấn/ha) có diện tích 629,56 ha; Cấp 2 (50 đến 100 tấn/ha) và Cấp 3 (100 đến 150 tấn/ha) chiếm phần diện tích còn lại. So với phân cực HH, diện tích nhóm 1 giảm từ 956,96 ha xuống 629,56 ha (chênh lệch hơn 51,9%), phản ánh đúng thực trạng phân bố mật độ cây rừng ngoài thực địa.

Thảo luận kết quả

Nguyên nhân dẫn đến sự khác biệt rõ rệt giữa các kênh phân cực bắt nguồn từ cơ chế vật lý tương tác sóng vi ba trong môi trường ngập triều. Rừng ngập mặn ven biển đồng bằng sông Hồng có nền đất ngập nước định kỳ, độ ẩm bùn cao tạo nên hằng số điện môi lớn. Kênh đơn phân cực VV phát sóng đứng và thu sóng đứng nên dễ bị chi phối bởi độ nhám bề mặt nước và bùn lầy, gây nhiễu loạn tín hiệu tán xạ từ tán cây. Ngược lại, việc thiết lập tỷ số phân cực HH/VV giúp triệt tiêu hoàn toàn tín hiệu phản xạ nền không mong muốn, giữ lại năng lượng tán xạ khối thuần túy phát ra từ cành, thân và vòm lá vẹt.

Trong thực tế xử lý dữ liệu viễn thám, kết quả này được trình bày trực quan thông qua biểu đồ phân tán hồi quy đa biến và bản đồ phân vùng sinh khối định dạng raster. Dữ liệu diện tích được tổng hợp thành bảng thống kê so sánh đa tiêu chí, giúp các nhà quản lý dễ dàng nhận diện sự dịch chuyển trữ lượng giữa các phân cấp rừng.

Khi so sánh với các công bố khoa học quốc tế tại vùng Guiana thuộc Pháp hay miền Bắc Australia, các phát hiện này hoàn toàn tương thích về ngưỡng bão hòa của sóng Radar băng C ở mức xấp xỉ 80 đến 150 tấn/ha. Đúng như khẳng định khoa học: "rừng là nguồn lưu trữ và tiêu thụ lượng các bon trong tự nhiên", kết quả luận văn đã chứng minh viễn thám Radar băng C là công cụ tối ưu để giám sát trữ lượng rừng ngập mặn mà không bị phụ thuộc vào điều kiện thời tiết hay mây mù nhiệt đới.

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

Nhằm nâng cao hiệu quả ứng dụng công nghệ viễn thám vi sóng trong giám sát tài nguyên và thích ứng với biến đổi khí hậu, luận văn đề xuất 4 giải pháp trọng tâm:

  1. Chuẩn hóa quy trình viễn thám Radar đa băng tần: Đề nghị Chi cục Kiểm lâm và Sở Nông nghiệp và Phát triển nông thôn các tỉnh Nam Định, Thái Bình tích hợp dữ liệu Radar băng C (Sentinel-1, ENVISAT) với Radar băng L (ALOS PALSAR) có bước sóng dài hơn 23 cm. Mục tiêu nâng ngưỡng đo lường sinh khối từ 150 tấn/ha lên trên 200 tấn/ha, đạt độ chính xác phân loại trên 90% trong giai đoạn 2026 - 2028.
  2. Mở rộng mạng lưới quan trắc cố định định vị vệ tinh: Nâng quy mô điều tra thực địa từ 13 ô mẫu lên 50 ô tiêu chuẩn cố định tại vùng lõi và vùng đệm Vườn quốc gia Xuân Thủy. Thực hiện đo đạc định kỳ 6 tháng/lần nhằm hiệu chỉnh liên tục các phương trình allometric, hoàn thành xây dựng mạng lưới trước quý IV năm 2027 do Ban quản lý Vườn quốc gia Xuân Thủy chủ trì.
  3. Thương mại hóa dữ liệu sinh khối phục vụ thị trường tín chỉ carbon: Ứng dụng bản đồ phân vùng sinh khối khô để tính toán chính xác trữ lượng phát thải hấp thụ CO2 tương đương theo chuẩn quốc tế. Mục tiêu gia tăng nguồn thu từ dịch vụ môi trường rừng thêm 20% đến 30% cho cộng đồng cư dân vùng đệm huyện Giao Thủy đến năm 2030, do Quỹ Bảo vệ và Phát triển rừng phối hợp với các tổ chức quốc tế thực hiện.
  4. Xây dựng hệ sinh thái WebGIS giám sát đới rừng phòng hộ ven biển: Thiết lập hệ thống cơ sở dữ liệu không gian thời gian thực phục vụ cảnh báo suy thoái rừng chắn sóng và sạt lở bãi bồi cửa Ba Lạt. Thời gian triển khai hoàn thiện trong 12 đến 18 tháng do Sở Tài nguyên và Môi trường tỉnh Nam Định chịu trách nhiệm vận hành.

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

  • Nhóm nghiên cứu và học viên cao học chuyên ngành Viễn thám, Hệ thống thông tin địa lý (GIS) và Quản lý tài nguyên môi trường: Cung cấp tài liệu tham khảo chuyên sâu về thuật toán xử lý ảnh Radar ASAR, kỹ thuật lọc nhiễu Frost và phương pháp xây dựng mô hình hồi quy vi sóng. Use case: Ứng dụng làm tài liệu giảng dạy và khung phương pháp luận cho các đề tài nghiên cứu định lượng tài nguyên rừng.
  • Ban quản lý các Vườn quốc gia và Khu bảo tồn đất ngập nước: Giúp nắm bắt công nghệ giám sát diện tích và chất lượng lâm phần từ xa với độ tin cậy cao. Use case: Xây dựng đề án bảo tồn đa dạng sinh học, theo dõi diễn thế rừng ngập mặn và bảo vệ bãi ăn của các loài chim di trú quý hiếm như Cò thìa.
  • Cơ quan quản lý nhà nước về Lâm nghiệp và Biến đổi khí hậu: Cung cấp công cụ định lượng nhanh trữ lượng gỗ và sinh khối thực vật phục vụ công tác quy hoạch đê điều, rừng phòng hộ. Use case: Thực hiện kiểm kê rừng định kỳ cấp tỉnh và lập báo cáo cam kết giảm phát thải khí nhà kính.
  • Các tổ chức phát triển quốc tế và doanh nghiệp tư vấn tín chỉ carbon: Cung cấp cơ sở khoa học độc lập để thẩm định dung lượng lưu trữ carbon của hệ sinh thái ven biển châu thổ sông Hồng. Use case: Xây dựng hồ sơ thiết kế dự án carbon rừng ngập mặn theo tiêu chuẩn thị trường tự nguyện.

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

Tại sao ảnh viễn thám Radar SAR lại vượt trội hơn ảnh quang học trong ước tính sinh khối rừng ngập mặn?

Ảnh quang học bị giới hạn bởi hiện tượng mây mù thường xuyên ở vùng biển nhiệt đới và chỉ ghi nhận phản xạ bề mặt lá hai chiều. Trong khi đó, sóng vi ba Radar là hệ thống viễn thám chủ động, có thể chụp ảnh xuyên mây trong mọi điều kiện thời tiết. Tín hiệu sóng băng C tương tác trực tiếp với các cấu trúc nhánh cây và vòm tán, phản ánh chính xác khối lượng sinh khối khô ba chiều của lâm phần.

Vì sao tỷ số phân cực HH/VV lại cho độ tin cậy cao hơn từng phân cực đơn lẻ?

Môi trường rừng ngập mặn có nền đất bùn lầy và thủy triều thay đổi liên tục, làm sai lệch năng lượng phản xạ ở các phân cực đơn HH hoặc VV. Phép chia tỷ số hai phân cực HH/VV giúp triệt tiêu các tác động nền gây nhiễu từ mặt nước và độ ẩm đất, làm nổi bật phản xạ tán xạ khối từ thân cành cây ngập mặn, từ đó tạo ra phương trình hồi quy có hệ số tương quan cao nhất.

Mô hình tính toán sinh khối trong luận văn áp dụng tối ưu cho ngưỡng trữ lượng nào?

Nghiên cứu chứng minh rằng bộ cảm Radar ENVISAT ASAR băng C đạt độ chính xác rất cao cho các quần thể rừng ngập mặn có sinh khối khô dưới 150 tấn/ha (tương đương với các thảm rừng Vẹt, Trang tái sinh hoặc trồng phòng hộ tại miền Bắc). Đối với các khu rừng cổ thụ có sinh khối vượt quá 150 tấn/ha đến 200 tấn/ha, tín hiệu băng C sẽ xuất hiện hiện tượng bão hòa, đòi hỏi phải phối hợp thêm bước sóng băng L.

Làm thế nào để loại bỏ triệt để hiện tượng nhiễu hạt trên ảnh viễn thám Radar?

Nghiên cứu sử dụng thuật toán lọc Frost dựa trên tỷ số giữa phương sai và giá trị trung bình cục bộ của các điểm ảnh trong cửa sổ trượt. Thuật toán này hoạt động như một bộ lọc thông thấp tại các vùng đồng nhất để khử các đốm sáng tối do sai pha sóng, đồng thời kích hoạt cơ chế bảo tồn biên tại các khu vực ranh giới giữa rừng ngập mặn và mặt nước biển.

Kết quả sinh khối từ luận văn có thể chuyển đổi thành chỉ số hấp thụ khí nhà kính CO2 như thế nào?

Sinh khối khô là tham số đầu vào trực tiếp để tính lượng carbon tích lũy. Theo quy chuẩn lâm nghiệp, lượng carbon nguyên chất chiếm khoảng 47% đến 50% tổng trọng lượng sinh khối khô. Bằng cách nhân diện tích 629,56 ha của cấp sinh khối ưu thế với hệ số chuyển đổi tương đương sang phân tử CO2 (nhân với hệ số 3,67), các nhà quản lý có thể xác định ngay tổng lượng phát thải được cô lập tại khu vực.

Kết luận

  • Khẳng định tính khả thi và độ tin cậy vượt trội của việc ứng dụng ảnh vệ tinh viễn thám Radar ENVISAT ASAR băng C trong việc xác định sinh khối rừng ngập mặn ven bờ châu thổ sông Hồng.
  • Xác lập mô hình hồi quy dựa trên tỷ số hai phân cực HH/VV (Y = -387,4 * HH/VV) với độ chính xác kiểm chứng thực địa đạt mức tuyệt đối trên các ô mẫu kiểm tra.
  • Thành lập thành công bản đồ phân vùng trữ lượng sinh khối rừng ngập mặn tỉnh Nam Định và Thái Bình, phân định rõ 3 cấp sinh khối với tổng diện tích nhóm ưu thế từ 0 đến 50 tấn/ha đạt 629,56 ha.
  • Đóng góp quy trình tiền xử lý ảnh Radar chuyên biệt kết hợp bộ lọc nhiễu thích ứng Frost và phân loại mặt nạ thảm phủ từ ảnh quang học SPOT 10 m.
  • Cung cấp nguồn tài liệu học thuật nền tảng phục vụ chuyển giao công nghệ giám sát tài nguyên rừng và chi trả dịch vụ môi trường rừng trong kỷ nguyên số.

Luận văn đã giải quyết thành công bài toán định lượng sinh khối rừng ngập mặn vùng cửa sông ven biển mà không phụ thuộc vào điều kiện thời tiết, tạo tiền đề khoa học vững chắc cho việc quản lý rừng bền vững. Trong giai đoạn 2 - 3 năm tới, hướng nghiên cứu cần tiếp tục mở rộng sang việc tích hợp dữ liệu Radar đa tần số băng L và băng P nhằm giám sát toàn diện các hệ sinh thái rừng nhiệt đới trên toàn quốc. Hãy tải và tham khảo ngay toàn văn công trình nghiên cứu xuất sắc này để làm chủ công nghệ viễn thám Radar trong phân tích môi trường hiện đại!