Tổng quan nghiên cứu

Tỉnh Bến Tre sở hữu diện tích tự nhiên khoảng 2.360 km² với hơn 65 km bờ biển tiếp giáp Biển Đông, nằm ở hạ lưu của 4 phân lưu sông Cửu Long gồm sông Tiền, Ba Lai, Hàm Luông và Cổ Chiên. Địa hình toàn tỉnh tương đối bằng phẳng và rất thấp trũng, trong đó cao độ dưới 1,0 m chiếm tới hơn 66% tổng diện tích tự nhiên, khiến địa phương này trở thành một trong những điểm nóng dễ bị tổn thương nhất do biến đổi khí hậu và nước biển dâng tại vùng Đồng bằng sông Cửu Long. Theo các số liệu quan trắc thực tế tại các trạm hải văn như An Thuận, Bình Đại, Bến Trại và Mỹ Hóa, tốc độ dâng lên của mực nước biển trung bình năm dao động từ 0,454 cm đến 1,076 cm/năm, làm gia tăng cường độ xâm nhập mặn sâu vào nội đồng và gây ngập lụt cục bộ.

Vấn đề cốt lõi đặt ra là sự suy giảm nghiêm trọng diện tích đất canh tác nông nghiệp truyền thống, đặc biệt là đất chuyên trồng lúa nước và cây ăn trái đặc sản, do tác động kép từ triều cường và nhiễm mặn. Mục tiêu trọng tâm của nghiên cứu là ứng dụng tích hợp công nghệ Hệ thống Thông tin Địa lý (GIS) và Viễn thám quang học đa thời gian để xây dựng cơ sở dữ liệu hiện trạng sử dụng đất qua các mốc năm 1995, 2005 và 2014. Qua đó, nghiên cứu tiến hành đánh giá định lượng mức độ ảnh hưởng của các kịch bản nước biển dâng do Bộ Tài nguyên và Môi trường ban hành đối với cơ cấu cây trồng của 8 đơn vị hành chính cấp huyện trong toàn tỉnh.

Nghiên cứu mang ý nghĩa khoa học và thực tiễn sâu sắc khi cung cấp bộ công cụ không gian số hóa chuẩn xác, định vị cụ thể các khu vực xung yếu. Theo kịch bản nước biển dâng 12 cm vào năm 2020, toàn tỉnh có khoảng 272,09 km² (chiếm 12,24% tổng diện tích) bị ngập úng, tạo tiền đề vững chắc cho công tác quy hoạch chuyển đổi mùa vụ và quản lý bền vững tài nguyên đất đai.

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 viễn thám đa phổ (Multispectral Remote Sensing) và nguyên lý phản xạ phổ của các đối tượng vỏ bọc mặt đất. Mỗi đối tượng tự nhiên như thực vật, thủy vực, đất trống và công trình nhân tạo đều sở hữu đặc trưng phản xạ quang phổ riêng biệt trên các dải sóng ánh sáng nhìn thấy (Visible), hồng ngoại gần (NIR) và hồng ngoại sóng ngắn (SWIR). Đặc biệt, hệ thống cảm biến quang học Landsat Thematic Mapper (TM), Enhanced Thematic Mapper Plus (ETM+) và Operational Land Imager (OLI) cung cấp các kênh phổ quan trọng phục vụ phân tách sinh khối thực vật và ranh giới mặt nước với độ phân giải không gian 30 m.

Bên cạnh đó, nghiên cứu kết hợp lý thuyết phân tích không gian trong GIS để mô hình hóa sự biến động lớp phủ đất (LULC - Land Use / Land Cover) và khung đánh giá tính dễ bị tổn thương theo kịch bản biến đổi khí hậu toàn cầu của Ban Liên chính phủ về Biến đổi Khí hậu (IPCC) cùng các kịch bản phát thải B1 (thấp), B2 (trung bình) và A1FI (cao). Các khái niệm then chốt được vận dụng bao gồm: thuật toán phân loại có kiểm định Maximum Likelihood (phân loại xác suất cực đại), ma trận nhầm lẫn (Confusion Matrix), hệ số tương đồng Kappa và thuật toán xử lý mất dữ liệu sọc Gapfill cho cảm biến ảnh vệ tinh.

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

Nguồn dữ liệu đầu vào bao gồm 3 bộ ảnh vệ tinh Landsat đa thời gian chụp khu vực Bến Tre (Path 125, Row 53) vào các năm 1995 (Landsat 5 TM), năm 2005 (Landsat 7 ETM+) và năm 2014 (Landsat 8 OLI) với tỷ lệ mây dưới 1%. Dữ liệu không gian được đối chiếu với bản đồ địa hình tỷ lệ 1:50.000, bản đồ hiện trạng sử dụng đất năm 2005 và số liệu niên giám thống kê giai đoạn 1995 - 2012 do Sở Tài nguyên và Môi trường tỉnh Bến Tre cung cấp.

Về quy mô mẫu và phương pháp chọn mẫu, nghiên cứu đã thiết lập mạng lưới với hơn 8 điểm khống chế mặt đất (Ground Control Points - GCPs) phân bố đồng đều trên các cù lao Bảo, cù lao Minh và cù lao An Hóa nhằm nắn chỉnh hình học, đảm bảo sai số toàn phương trung bình (RMS Error) luôn nhỏ hơn 1 pixel (tương đương dưới 30 m). Đồng thời, hàng trăm điểm mẫu khảo sát thực địa bằng thiết bị định vị GPS cầm tay đã được thu thập đại diện cho 7 loại hình sử dụng đất chính: đất lúa, đất cây lâu năm, đất nuôi trồng thủy sản, đất lâm nghiệp, đất đô thị - công trình, đất chưa sử dụng và mặt nước sông ngòi.

Quy trình phân tích được thực hiện tuần tự: tiền xử lý ảnh và chạy thuật toán Gapfill trên nền tảng ENVI để loại bỏ lỗi sọc cảm biến SLC-off của ảnh năm 2005; phân loại ảnh tự động bằng thuật toán Maximum Likelihood nhờ khả năng tính toán hàm mật độ xác suất chuẩn xác cho từng nhóm pixel; đánh giá độ chính xác phân loại thông qua ma trận sai số; chuyển đổi dữ liệu raster sang định dạng vector shapefile và tiến hành phân tích không gian giao cắt (Spatial Overlay) trên phần mềm ArcMap nhằm xác lập biến động cơ cấu sử dụng đất trong giai đoạn 19 năm (1995 - 2014).

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

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

Thứ nhất, diện tích đất trồng lúa nước tại tỉnh Bến Tre ghi nhận mức độ sụt giảm chưa từng có. Theo kết quả giải đoán ảnh vệ tinh và thống kê địa phương, diện tích đất chuyên trồng lúa từ mức 140.857 ha vào năm 2005 đã sụt giảm mạnh mẽ chỉ còn 22.887 ha vào năm 2014, tương ứng mức giảm kỷ lục hơn 83,7% chỉ trong vòng một thập kỷ. Sự suy giảm này tập trung rõ rệt nhất tại các vùng trũng giáp ranh thuộc huyện Ba Tri, Giồng Trôm và Thạnh Phú.

Thứ hai, diện tích đất nuôi trồng thủy sản ghi nhận xu hướng tăng trưởng bùng nổ, mở rộng từ hơn 8.000 ha vào năm 1995 lên xấp xỉ 50.000 ha vào năm 2014, tăng gấp hơn 6 lần. Hoạt động chuyển đổi mặt nước phục vụ nuôi tôm sú, tôm thẻ chân trắng và nghêu sò diễn ra mạnh mẽ tại 3 huyện ven biển gồm Ba Tri, Bình Đại và Thạnh Phú, phản ánh sự thích ứng tự phát của cộng đồng cư dân trước tình trạng nguồn nước ngọt bị suy giảm.

Thứ ba, mô hình hóa kịch bản ngập lụt cho thấy mối đe dọa trực tiếp của mực nước biển dâng. Khi mực nước biển dâng 12 cm theo kịch bản A1FI vào năm 2020, tổng diện tích ngập lụt ước tính là 272,09 km² (chiếm 12,24% diện tích toàn tỉnh). Trong đó, huyện Ba Tri chịu tổn thương nặng nề nhất do có cao trình phổ biến từ 0,5 m đến 0,75 m, kế đến là huyện Thạnh Phú với hệ thống đê bao chịu áp lực triều cường từ cửa sông Cổ Chiên và Hàm Luông.

Thứ tư, độ chính xác tổng thể của các bản đồ phân loại hiện trạng sử dụng đất từ ảnh Landsat qua các năm đều đạt trên 85%, với hệ số Kappa vượt ngưỡng 0,82, khẳng định tính khả thi và độ tin cậy vượt trội của chu trình xử lý ảnh viễn thám kết hợp kiểm chứng thực địa.

Thảo luận kết quả

Nguyên nhân cốt lõi dẫn đến sự chuyển dịch cơ cấu sử dụng đất mạnh mẽ xuất phát từ tác động kép của thủy văn tự nhiên và yếu tố kinh tế xã hội. Về mặt tự nhiên, địa hình toàn tỉnh nghiêng dần ra biển với hơn 46% diện tích nằm ở dải cao độ 0,75 - 1,0 m. Khi mực nước biển dâng kết hợp biên độ bán nhật triều Biển Đông cao, ranh mặn 4‰ xâm nhập sâu hơn 40 - 50 km vào các nhánh sông chính, biến các cánh đồng lúa nước ngọt truyền thống thành vùng đất nhiễm mặn không thể tiếp tục canh tác. Về mặt xã hội, nông dân chủ động đào ao chuyển dịch sang mô hình nuôi tôm nước lợ để tìm kiếm lợi nhuận kinh tế cao hơn.

Các kết quả nghiên cứu được trực quan hóa sinh động thông qua biểu đồ hình cột so sánh biến động diện tích qua 3 thời kỳ 1995, 2005, 2014 và bảng ma trận chuyển dịch đất đai (Transition Matrix) trong GIS. Bản đồ không gian thể hiện rõ nét ranh giới các mảng màu biến đổi: sắc xanh của đất lúa lùi dần về phía Tây Bắc (huyện Chợ Lách, Châu Thành), nhường chỗ cho sắc tím của đầm nuôi thủy sản lan rộng từ duyên hải vào trung tâm tỉnh. Phát hiện này hoàn toàn tương đồng với các nghiên cứu quan trắc đất ngập nước tại các tỉnh ven biển lân cận như Trà Vinh, Sóc Trăng và Cà Mau, chứng minh tính quy luật của quá trình chuyển dịch sinh kế thích ứng mặn tại đồng bằng hạ lưu sông Mekong.

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

Thứ nhất, khẩn trương gia cố, nâng cấp và xây mới đồng bộ hệ thống đê bao ngăn triều, cống đập điều tiết xâm nhập mặn và giữ ngọt tại các vùng trũng xung yếu thuộc huyện Ba Tri, Giồng Trôm và Bình Đại. Mục tiêu cụ thể là bảo vệ an toàn cho hơn 20.000 ha vùng chuyên canh cây ăn trái đặc sản và hoa màu giá trị cao trước năm 2030, do Sở Nông nghiệp và Phát triển nông thôn phối hợp với Công ty TNHH MTV Khai thác Công trình Thủy lợi tỉnh Bến Tre chủ trì thực hiện.

Thứ hai, đẩy mạnh tái cơ cấu ngành nông nghiệp theo hướng thích ứng thông minh với biến đổi khí hậu (Climate-Smart Agriculture). Triển khai nhân rộng mô hình luân canh lúa - tôm, nuôi thủy sản sinh thái và phát triển các giống cây trồng chịu hạn mặn như dừa công nghiệp, bưởi da xanh ghép gốc chống chịu mặn tại 4 huyện cù lao Minh giai đoạn 2025 - 2035. Chủ thể thực hiện là Sở Nông nghiệp và Phát triển nông thôn liên kết với Viện Cây ăn quả Miền Nam và Hội Nông dân các cấp, nhằm bảo đảm sinh kế bền vững cho hơn 150.000 hộ gia đình nông thôn.

Thứ ba, bảo vệ nghiêm ngặt và phục hồi vành đai rừng phòng hộ ven biển tại 3 huyện Ba Tri, Bình Đại và Thạnh Phú với chỉ tiêu trồng mới tối thiểu 100 đến 150 ha rừng ngập mặn mỗi năm trong giai đoạn 2025 - 2030. Giải pháp này giúp gia tăng khả năng tiêu tán năng lượng sóng biển, chống xói lở bờ biển và bồi tụ bãi triều, do Chi cục Kiểm lâm tỉnh Bến Tre phối hợp cùng chính quyền các xã ven biển triển khai.

Thứ tư, xây dựng và chuẩn hóa hệ thống cơ sở dữ liệu số GIS - Viễn thám giám sát biến động tài nguyên đất tự động với chu kỳ cập nhật định kỳ 12 tháng/lần. Sở Tài nguyên và Môi trường tỉnh Bến Tre chịu trách nhiệm vận hành, tích hợp trực tiếp dữ liệu này vào Cổng thông tin quy hoạch không gian của tỉnh nhằm hỗ trợ quá trình ra quyết định cấp phép chuyển mục đích sử dụng đất đến năm 2030.

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

Thứ nhất, các nhà quản lý và hoạch định chính sách tại Sở Tài nguyên và Môi trường, Sở Nông nghiệp và Phát triển nông thôn tỉnh Bến Tre cũng như các tỉnh ven biển Đồng bằng sông Cửu Long. Nghiên cứu cung cấp bức tranh dữ liệu không gian định lượng giúp tối ưu hóa công tác lập quy hoạch sử dụng đất và kế hoạch hành động ứng phó biến đổi khí hậu tầm nhìn 2030 - 2050.

Thứ hai, các chuyên gia phân tích dữ liệu không gian, kỹ sư GIS và chuyên viên viễn thám môi trường. Luận văn là tài liệu tham khảo chi tiết về quy trình kỹ thuật tiền xử lý ảnh vệ tinh Landsat đa thời gian, phương pháp khắc phục lỗi sọc SLC-off bằng thuật toán Gapfill và kỹ thuật phân loại Maximum Likelihood trên phần mềm ENVI và ArcMap.

Thứ ba, các giảng viên, nhà nghiên cứu, học viên cao học và sinh viên thuộc các chuyên ngành Quản lý Môi trường, Biến đổi Khí hậu, Khoa học Đất, Địa lý Tài nguyên và Quản lý Đất đai. Tài liệu này cung cấp phương pháp luận chặt chẽ kết hợp giữa khoa học không gian và các kịch bản phát thải khí nhà kính toàn cầu.

Thứ tư, các tổ chức phi chính phủ (NGOs), tổ chức phát triển quốc tế và các doanh nghiệp đầu tư nông - thủy sản. Luận văn hỗ trợ đánh giá rủi ro thiên tai, hỗ trợ thiết kế các dự án viện trợ sinh kế dựa vào cộng đồng và định hướng vùng đầu tư nuôi trồng thủy hải sản bền vững tại 3 huyện ven biển Bến Tre.

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

Tại sao ảnh vệ tinh Landsat 7 năm 2005 gặp sự cố sọc đen và giải pháp khắc phục trong nghiên cứu là gì?

Từ ngày 31 tháng 3 năm 2003, bộ phận hiệu chỉnh đường quét (Scan Line Corrector - SLC) trên vệ tinh Landsat 7 bị hỏng khiến toàn bộ ảnh chụp sau đó bị mất khoảng 22% dữ liệu dưới dạng các vệt sọc. Nghiên cứu đã áp dụng thành công thuật toán Gapfill do NASA phát triển, lồng ghép ảnh đối chứng năm 2002 để khôi phục hoàn chỉnh dữ liệu phổ trước khi giải đoán.

Mực nước biển dâng 12 cm theo kịch bản năm 2020 ảnh hưởng cụ thể như thế nào đến diện tích tỉnh Bến Tre?

Theo tính toán từ mô hình số độ cao DEM và kịch bản phát thải cao A1FI, mức nước dâng 12 cm làm ngập lụt trực tiếp 272,09 km² đất đai, chiếm 12,24% diện tích tự nhiên của tỉnh. Khu vực chịu thiệt hại nặng nhất là huyện Ba Tri do có địa hình trũng thấp và huyện Thạnh Phú giáp hai cửa sông lớn.

Nguyên nhân cốt lõi khiến diện tích trồng lúa tại Bến Tre giảm mạnh hơn 80% từ năm 2005 đến 2014 là gì?

Sự sụt giảm từ 140.857 ha xuống 22.887 ha bắt nguồn từ tác động cộng hưởng của hiện tượng xâm nhập mặn mùa khô kéo dài khiến canh tác lúa không đạt hiệu quả, kết hợp với xu hướng người dân chủ động chuyển đổi sang mô hình nuôi tôm nước lợ với diện tích toàn tỉnh đạt gần 50.000 ha vào năm 2014.

Phương pháp phân loại Maximum Likelihood trên phần mềm ENVI đạt độ tin cậy ra sao?

Phương pháp này đạt độ chính xác tổng thể trên 85% và hệ số Kappa trên 0,82 cho cả 3 cảnh ảnh. Kết quả phân loại được kiểm chứng chặt chẽ thông qua việc đối chiếu với hơn 8 điểm khống chế mặt đất, bản đồ hiện trạng sử dụng đất thực tế và các điểm định vị GPS khảo sát thực địa.

Ba huyện nào của tỉnh Bến Tre chịu mức độ tổn thương nặng nề nhất do biến đổi khí hậu và nước biển dâng?

Ba huyện ven biển gồm Ba Tri, Bình Đại và Thạnh Phú chịu tổn thương cao nhất. Các địa phương này sở hữu đường bờ biển dài 65 km, cao trình mặt đất trung bình chỉ từ 0,5 m đến 1,0 m và chịu ảnh hưởng trực tiếp bởi chế độ bán nhật triều từ các cửa sông Hàm Luông, Cổ Chiên và Ba Lai.

Kết luận

  • Số hóa và thành lập thành công bản đồ hiện trạng sử dụng đất tỉnh Bến Tre qua 3 mốc thời gian 1995, 2005 và 2014 từ dữ liệu vệ tinh Landsat với độ chính xác tổng thể trên 85%.
  • Định lượng rõ nét sự dịch chuyển cơ cấu nông nghiệp: diện tích đất trồng lúa giảm mạnh từ 140.857 ha (năm 2005) xuống 22.887 ha (năm 2014), trong khi diện tích nuôi trồng thủy sản tăng vọt lên gần 50.000 ha.
  • Đánh giá nguy cơ tổn thương không gian theo kịch bản nước biển dâng: mức dâng 12 cm vào năm 2020 gây ngập 272,09 km² (12,24% diện tích toàn tỉnh), tập trung tại các vùng trũng huyện Ba Tri và Thạnh Phú.
  • Xây dựng hoàn chỉnh cơ sở dữ liệu GIS chuyên đề làm nền tảng khoa học phục vụ công tác quản lý tài nguyên và dự báo biến động đất đai cho 8 đơn vị hành chính.
  • Đề xuất hệ thống 4 nhóm giải pháp kỹ thuật công trình đê cống và chuyển đổi sinh kế thông minh thích ứng mặn giai đoạn 2025 - 2035.

Đóng góp lớn nhất của luận văn là chứng minh hiệu quả vượt trội của việc tích hợp GIS và công nghệ Viễn thám trong giám sát diễn biến môi trường tự nhiên với chi phí tối ưu. Trong giai đoạn tiếp theo từ năm 2025 đến năm 2030, các cơ quan quản lý cần tích hợp thêm ảnh vệ tinh độ phân giải cao như Sentinel-2 và dữ liệu radar ngập lụt nhằm hoàn thiện hệ thống cảnh báo sớm. Hãy chủ động ứng dụng công nghệ không gian số ngay hôm nay để bảo vệ nguồn tài nguyên đất đai và kiến tạo sinh kế bền vững cho vùng đồng bằng châu thổ trước thách thức biến đổi khí hậu toàn cầu.