Tổng quan nghiên cứu

Vịnh Tiên Yên, còn được gọi là vùng vịnh Tiên Yên - Hà Cối thuộc địa phận tỉnh Quảng Ninh, là một thủy vực ven bờ đặc biệt nằm trong tọa độ địa lý từ 21°15'00'' đến 21°30'00'' vĩ độ Bắc và 107°27'00'' đến 108°21'35'' kinh độ Đông. Vùng vịnh sở hữu chiều dài khoảng 57 km, chiều rộng trung bình 9 km và diện tích mặt nước ngập thường xuyên đạt xấp xỉ 400 km2. Nơi đây là trung tâm hội tụ của 13 kiểu đất ngập nước sinh thái với hơn 754 loài sinh vật biển, bao gồm 9.903 ha rừng ngập mặn cùng hệ sinh thái cỏ biển quý hiếm.

Tuy nhiên, áp lực kép từ các quá trình xói lở tự nhiên với tốc độ trung bình 5 mm/năm và hoạt động kinh tế nhân sinh như việc chuyển đổi 272,9 ha rừng ngập mặn sang đầm nuôi thủy sản giai đoạn 2000-2009 đã làm gia tăng đáng kể lượng trầm tích lơ lửng, đe dọa trực tiếp chất lượng nước vịnh. Độ đục nước biển là chỉ tiêu cốt lõi phản ánh mức độ tích tụ bùn cát và lan truyền ô nhiễm, song việc quan trắc truyền thống bằng tàu thuyền trên diện tích 400 km2 thường tiêu tốn kinh phí lớn và hạn chế về độ bao phủ không gian.

Nghiên cứu được xây dựng nhằm giải quyết bài toán cấp thiết trên thông qua việc tích hợp ảnh viễn thám vệ tinh đa thời gian và mô hình hệ thống thông tin địa lý GIS để lượng hóa, lập bản đồ phân bố độ đục nước vịnh trong tháng 7 giai đoạn 2006-2011. Kết quả nghiên cứu không chỉ thiết lập phương trình tương quan thực nghiệm có độ tin cậy cao mà còn chứng minh giải pháp viễn thám giúp giảm hơn 65% chi phí và rút ngắn 80% thời gian quan trắc định kỳ so với các phương pháp trắc địa thông thườ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 ứng dụng lý thuyết quang học thủy văn và mô hình phản xạ phổ điện từ của các hợp phần trong nước biển. Khi bức xạ mặt trời chiếu xuống mặt nước, các hạt vật chất lơ lửng và độ đục sẽ làm tăng đáng kể hệ số tán xạ ngược tại dải sóng ánh sáng nhìn thấy và cận hồng ngoại. Dựa trên nguyên lý này, mối quan hệ hàm số giữa hệ số phản xạ phổ kênh đỏ của cảm biến vệ tinh và nồng độ chất rắn lơ lửng được thiết lập làm nền tảng toán học.

Bên cạnh đó, lý thuyết địa thống kê không gian được áp dụng để tối ưu hóa việc mô hình hóa phân bố độ đục liên tục. Ba khái niệm chuyên ngành then chốt chi phối cấu trúc nghiên cứu bao gồm: độ đục quang học biểu thị mức độ tán xạ ánh sáng của nước, tổng chất rắn lơ lửng (TSS) đại diện cho khối lượng hạt trầm tích mịn, và phương pháp nội suy không gian Kriging dựa trên mô hình hàm bán phương sai (Semivariogram).

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

Nguồn dữ liệu đầu vào sử dụng chuỗi ảnh vệ tinh quang học MODIS độ phân giải không gian 250 m thu nhận vào tháng 7 và tháng 8 từ năm 2006 đến năm 2011, kết hợp số liệu đo đạc thực địa độ đục và tổng chất rắn lơ lửng tại 40 điểm khảo sát vào ngày 06/07/2010. Phương pháp chọn mẫu phân tầng ngẫu nhiên được triển khai dọc theo các mặt cắt từ cửa sông Tiên Yên, sông Ba Chẽ ra các luồng lạch sâu từ 2 m đến 30 m tại cửa biển.

Quy trình phân tích dữ liệu được thực hiện qua các bước chuẩn hóa nghiêm ngặt:

  • Xử lý tiền kỳ ảnh viễn thám trên phần mềm ENVI bao gồm hiệu chỉnh khí quyển và nắn chỉnh hình học tọa độ.
  • Thiết lập phương trình hồi quy tương quan tuyến tính giữa giá trị độ phản xạ phổ kênh 1 MODIS (dải sóng 620 - 670 nm) với số liệu độ đục đo thực tế.
  • Ứng dụng mô hình địa thống kê Ordinary Kriging (OK) và Co-Kriging (CK) trong phần mềm ArcGIS để nội suy và thành lập bản đồ phân bố độ đục trên toàn bộ mặt vịnh. Lý do lựa chọn mô hình Co-Kriging là khả năng kết hợp dữ liệu ảnh viễn thám làm biến phụ trợ liên tục, giúp tối thiểu hóa phương sai sai số dự báo tại các khu vực thiếu trạm đo trực tiếp.

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

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

Nghiên cứu đã xây dựng thành công hàm thực nghiệm phản ánh mối quan hệ chặt chẽ giữa độ phản xạ phổ kênh 1 ảnh MODIS và độ đục nước vịnh Tiên Yên với hệ số tương quan đạt mức cao, khẳng định kênh phổ đỏ có độ nhạy vượt trội trong việc phát hiện nồng độ hạt bùn cát mịn lơ lửng ven bờ.

Phân bố độ đục tại vịnh Tiên Yên thể hiện tính phân hóa không gian vô cùng rõ nét:

  • Vùng nước nông ven bờ tiếp giáp các cửa sông Tiên Yên, Hà Cối và Ba Chẽ có độ đục cao nhất, giá trị thường xuyên dao động từ 25 NTU đến trên 50 NTU.
  • Khu vực vịnh kín và các bãi triều cát bùn trung gian duy trì độ đục ở mức trung bình từ 12 NTU đến 24 NTU.
  • Vùng nước sâu ngoài cửa Mô và cửa Đại được che chắn bởi hệ thống đảo Cái Bầu, Vạn Vược có độ đục thấp nhất, giảm mạnh xuống dưới 8 NTU.

Xu thế biến động theo thời gian trong tháng 7 giai đoạn 2006-2011 cho thấy độ đục nước mặt có xu hướng tăng trung bình từ 12% đến 18% tại các khu vực tập trung hoạt động san lấp mặt bằng khu công nghiệp cảng biển Hải Hà và các xã ven biển có diện tích đầm tôm gia tăng. Đồng thời, việc ứng dụng thuật toán Co-Kriging kết hợp biến ảnh vệ tinh đã giúp giảm phương sai sai số nội suy hơn 25% so với phương pháp Ordinary Kriging thông thường.

Thảo luận kết quả

Sự gia tăng độ đục tại các khu vực cửa sông bắt nguồn từ đặc điểm thủy văn - hải văn đặc thù của lưu vực. Tổng lượng dòng chảy hàng năm từ sông Tiên Yên (1,1 km3) và sông Ba Chẽ (0,96 km3) đổ dồn vào mùa mưa lũ kết hợp lượng mưa trung bình năm đạt 2.537,6 mm cuốn theo khối lượng lớn vật liệu xói mòn từ vùng núi thượng nguồn. Thêm vào đó, chế độ nhật triều thuần nhất với biên độ triều cực đại lên tới 4,0 m tạo nên vận tốc dòng triều rút mạnh, đạt 74,2 cm/s tại cửa Mô và 49,3 cm/s tại cửa Đại, liên tục khuấy động tầng trầm tích đáy bùn sét Holocen dày từ 1 m đến 10 m.

So sánh với các nghiên cứu thủy văn quang học tại các vịnh ven bờ lân cận, kết quả nghiên cứu hoàn toàn tương đồng về quy luật phân dị độ đục từ cửa sông ra biển khơi. Các dữ liệu biến động không gian - thời gian này được mô hình hóa trực quan qua hệ thống biểu đồ phân tán hồi quy phổ, bản đồ phân bố phương sai sai số lưới ô vuông 250 m x 250 m, và bảng thống kê diện tích biến động nồng độ hạt lơ lửng, cung cấp bức tranh định lượng hoàn chỉnh cho công tác giám sát thủy vực.

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

Nhằm nâng cao hiệu quả quản trị môi trường nước và bảo vệ tính toàn vẹn sinh thái vịnh Tiên Yên, các giải pháp trọng tâm cần được thực thi đồng bộ:

  • Thiết lập hệ thống giám sát viễn thám tự động hóa: Chi cục Bảo vệ Môi trường tỉnh Quảng Ninh cần chủ trì triển khai phần mềm quan trắc trực tuyến tích hợp ảnh vệ tinh định kỳ 15 ngày/lần, hoàn thành khung vận hành trong giai đoạn 2024-2026 nhằm kiểm soát sớm các điểm nóng gia tăng độ đục bất thường.
  • Phục hồi vành đai rừng ngập mặn phòng hộ: Sở Nông nghiệp và Phát triển Nông thôn phối hợp với Ủy ban nhân dân huyện Tiên Yên thực hiện kế hoạch tái sinh tối thiểu 300 ha rừng ngập mặn tại các xã Đồng Rui và Hải Lạng trước năm 2027, đặt mục tiêu giảm 30% lượng phù sa xói mòn đổ ra vịnh.
  • Kiểm soát nghiêm ngặt nguồn xả thải nuôi trồng thủy sản: Bắt buộc 100% các cơ sở thâm canh tôm trên 650 ha mặt nước tại huyện Hải Hà và Đầm Hà phải xây dựng hồ lắng sinh học đạt tiêu chuẩn trước năm 2025, hướng đến giảm 40% lượng chất thải hữu cơ và hạt lơ lửng xả thẳng ra bãi triều.
  • Xây dựng công trình kè sinh thái và nạo vét luồng hàng hải: Ban Quản lý dự án giao thông thủy tỉnh thực hiện nạo vét duy tu 15 km luồng lạch tại cảng Mũi Chùa và Đầm Buôn, đồng thời gia cố kè mềm sinh thái tại các đoạn bờ xói lở 5 mm/năm thuộc xã Quảng Điền trong lộ trình 2025-2028.

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

Nội dung và kết quả thực nghiệm của luận văn mang lại giá trị thực tiễn và học thuật cho các nhóm đối tượng:

  • Cơ quan quản lý nhà nước về môi trường và tài nguyên biển: Sở Tài nguyên và Môi trường, Chi cục Biển và Hải đảo khai thác bản đồ độ đục và dữ liệu 13 kiểu đất ngập nước để hoạch định không gian biển, quy hoạch phân vùng bảo tồn và cấp phép dự án cảng biển.
  • Chuyên gia, nghiên cứu sinh và học viên ngành Địa tin học, Khoa học Môi trường: Tham khảo quy trình tiền xử lý ảnh quang học MODIS, thuật toán hồi quy kênh phổ đỏ và mô hình địa thống kê Co-Kriging ứng dụng cho các thủy vực nông ven bờ.
  • Doanh nghiệp và cộng đồng nuôi trồng thủy sản ven biển: Sử dụng thông tin phân bố độ đục và chu kỳ dòng chảy để lựa chọn tọa độ đặt bè nuôi cá, khoanh vùng bãi thả ngao giống nhằm tránh hiện tượng bồi lắng gây ngạt thủy sản.
  • Các tổ chức bảo tồn thiên nhiên và biến đổi khí hậu: Kế thừa số liệu điều tra về sự suy giảm 394,44 ha rừng ngập mặn và thảm cỏ biển để xây dựng các đề án tài trợ phục hồi sinh thái biển Đông Bắc Việt Nam.

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

Vì sao kênh 1 của ảnh vệ tinh MODIS lại tối ưu trong việc đánh giá độ đục nước biển?

Kênh 1 của cảm biến MODIS hoạt động ở dải bước sóng đỏ 620 - 670 nm, đây là vùng quang phổ có độ nhạy quang học cao nhất với sự tán xạ ánh sáng của các hạt trầm tích mịn lơ lửng trong môi trường nước nông, giúp phản ánh chân thực nồng độ bùn cát mà ít bị hấp thụ bởi các phân tử nước tinh khiết.

Sai số phương sai giữa phương pháp Ordinary Kriging và Co-Kriging chênh lệch như thế nào?

Phương pháp Co-Kriging sử dụng thêm kênh phổ ảnh MODIS độ phân giải 250 m làm biến phụ trợ liên tục, giúp giảm hơn 25% phương sai sai số dự báo so với Ordinary Kriging chỉ dựa vào 40 điểm đo thực địa rời rạc, từ đó nâng cao tính chính xác của bản đồ nội suy.

Yếu tố tự nhiên chính nào làm gia tăng độ đục nước vịnh Tiên Yên vào mùa mưa?

Tổng lưu lượng nước lũ trên 2,06 km3 từ hệ thống sông Tiên Yên và sông Ba Chẽ kết hợp lượng mưa vượt 2.537,6 mm cuốn trôi đất đá xói mòn xuống hạ lưu, cùng với biên độ nhật triều cực đại 4,0 m khuấy đảo lớp trầm tích đáy bùn sét là các nguyên nhân chính gây đục nước.

Sự suy giảm diện tích rừng ngập mặn tác động thế nào đến độ trong của nước?

Việc mất đi 394,44 ha rừng ngập mặn làm suy giảm khả năng giữ hạt trầm tích của hệ thống rễ cây phòng hộ, khiến bùn cát dễ dàng bị sóng triều bào mòn và phân tán rộng ra các vụng vịnh, làm tăng độ đục cục bộ từ 12% đến 18%.

Quy trình công nghệ của luận văn có thể chuyển giao cho các thủy vực khác không?

Quy trình tích hợp viễn thám quang học và mô hình địa thống kê trong luận văn hoàn toàn có thể hiệu chỉnh hệ số tương quan để áp dụng hiệu quả cho các vùng đầm phá, cửa sông ven biển miền Trung và Đồng bằng sông Cửu Long với chi phí tiết kiệm.

Kết luận

  • Xây dựng thành công quy trình công nghệ kết hợp ảnh vệ tinh MODIS 250 m và mô hình địa thống kê Co-Kriging trong việc giám sát độ đục nước mặt quy mô lớn.
  • Xác lập phương trình hồi quy thực nghiệm tin cậy giữa độ phản xạ kênh phổ đỏ và độ đục thực địa tại 40 trạm khảo sát trên toàn vùng vịnh.
  • Hoàn thiện bộ bản đồ chuyên đề đánh giá biến động độ đục theo không gian và thời gian tháng 7 giai đoạn 2006-2011 trên diện tích 400 km2.
  • Chứng minh tính ưu việt của mô hình Co-Kriging trong việc giảm 25% sai số nội suy và tiết kiệm 65% ngân sách quan trắc định kỳ so với phương pháp cổ điển.
  • Cung cấp luận cứ khoa học vững chắc về mối liên hệ giữa thủy triều, lưu lượng sông ngòi và các áp lực nhân sinh đối với chất lượng môi trường nước biển.

Đóng góp then chốt của công trình là mở ra hướng tiếp cận định lượng hiện đại, hỗ trợ hiệu quả cho việc kiểm soát tài nguyên biển. Trong giai đoạn tiếp theo đến năm 2030, việc tích hợp nguồn ảnh có độ phân giải cao như Sentinel-2 và Landsat 8/9 sẽ tiếp tục hoàn thiện khung cảnh báo môi trường tự động. Các cơ quan chuyên trách cần nhanh chóng ứng dụng giải pháp này vào thực tiễn quản lý để hướng tới mục tiêu phát triển kinh tế biển bền vững.