Tổng quan nghiên cứu

Vịnh Tiên Yên thuộc tỉnh Quảng Ninh là vùng biển ven bờ có vị thế địa lý đặc biệt với diện tích phần nước ngập thường xuyên khoảng 400 km2, chiều dài 57 km và chiều rộng trung bình 9 km. Khu vực này sở hữu hơn 20.000 ha bãi triều cùng hệ thống sinh thái phong phú gồm 9.903 ha rừng ngập mặn và hàng trăm hecta cỏ biển. Tuy nhiên, áp lực phát triển kinh tế xã hội từ 134,7 nghìn dân cư trên địa bàn 3 huyện Tiên Yên, Đầm Hà, Hải Hà cùng hoạt động nuôi trồng thủy sản trên 2.900,7 ha đầm nuôi đã tạo ra những tác động tiêu cực đến môi trường nước biển ven bờ.

Vấn đề cốt lõi được đặt ra là nhu cầu giám sát chất lượng môi trường nước thường xuyên, trong đó độ đục là thông số vật lý then chốt phản ánh sự vận chuyển trầm tích và ô nhiễm môi trường. Các phương pháp trắc địa truyền thống bằng thuyền lấy mẫu tốn kém kinh phí, tiêu tốn nhiều thời gian và không thể bao quát toàn bộ 400 km2 diện tích mặt nước trong cùng một thời điểm.

Mục tiêu của nghiên cứu là xây dựng quy trình tích hợp công nghệ viễn thám đa thời gian từ vệ tinh MODIS và hệ thống thông tin địa lý GIS để tính toán định lượng độ đục, thành lập bản đồ phân bố không gian và phân tích quy luật biến động độ đục nước vịnh Tiên Yên trong các tháng 7 giai đoạn 2006 - 2011.

Phạm vi nghiên cứu tập trung vào toàn bộ vùng nước mặt vịnh Tiên Yên và các cửa sông chính như sông Tiên Yên, sông Ba Chẽ, sông Hà Cối và Ka Long. Kết quả nghiên cứu mang lại ý nghĩa thực tiễn to lớn khi giúp giảm thiểu trên 40% chi phí khảo sát thực địa định kỳ, nâng cao hiệu quả giám sát không gian lên 100% diện tích mặt vịnh và cung cấp cơ sở khoa học chính xác cho công tác quản lý tài nguyên môi trường ven biể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 dựa trên sự kết hợp chặt chẽ giữa quang học thủy văn và địa thống kê không gian:

  1. Lý thuyết quang học môi trường nước ven bờ: Nước biển ven bờ thuộc nhóm thủy vực phức tạp loại 2 (Case 2 waters), nơi tính chất quang học bị chi phối bởi các thành phần vật chất lơ lửng, chất hữu cơ hòa tan và diệp lục tố. Phổ phản xạ bức xạ điện từ của nước mặt ở dải sóng ánh sáng đỏ (kênh 1 của MODIS, bước sóng 0,620 - 0,670 µm) có mối tương quan thuận rõ rệt với hàm lượng tổng chất rắn lơ lửng (TSS) và độ đục của nước.

  2. Lý thuyết mô hình địa thống kê không gian (Geostatistics): Kỹ thuật sử dụng toán đề bán phần (semi-variogram) để định lượng cấu trúc không gian và tính biến thiên dị hướng của các thông số môi trường qua ba đại lượng đặc trưng: range (khoảng tương quan), sill (độ biến thiên cực đại) và nugget effect (hiệu ứng khối mẫu).

Các khái niệm cốt lõi bao gồm: Độ đục quang học (Nephelometric Turbidity Units - NTU), phổ phản xạ bề mặt sau hiệu chỉnh khí quyển, phương pháp nội suy Kriging thông thường (Ordinary Kriging - OK) và mô hình hồi quy viễn thám kết hợp địa thống kê đa biến Co-Kriging (CK).

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

Nguồn dữ liệu của luận văn bao gồm 6 cảnh ảnh vệ tinh MODIS/Terra mức 1B độ phân giải không gian 250 m được thu thập vào tháng 7 hoặc tháng 8 các năm từ 2006 đến 2011 (cụ thể: 20/07/2006, 14/07/2007, 20/08/2008, 10/07/2009, 06/07/2010 và 03/07/2011) tải từ hệ thống NASA. Vệ tinh Terra bay qua vùng nghiên cứu vào khoảng 10:30 giờ Việt Nam, tạo điều kiện thuận lợi cho việc đối soát dữ liệu mặt đất.

Phương pháp khảo sát thực địa được thực hiện với cỡ mẫu 40 điểm đo đạc trực tiếp trên toàn vịnh Tiên Yên vào ngày 06/07/2010. Phương pháp chọn mẫu là phân bố ngẫu nhiên phân tầng theo mạng lưới lưới ô vuông đại diện cho các vùng sinh thái cửa sông, trung tâm vịnh và luồng lạch ngoài khơi. Độ đục tại 40 điểm được đo bằng máy đo chất lượng nước Horiba Water Quality Meter U-10 kết hợp máy định vị GPS cầm tay, thực hiện đồng thời trong khung thời gian trước và sau 6 giờ khi vệ tinh Terra chụp ảnh.

Quy trình xử lý ảnh được tiến hành trên phần mềm ENVI 4.8 và ArcGIS 10:

  • Hiệu chỉnh hình học: Sử dụng công cụ MODIS Conversion Toolkit để đưa ảnh về hệ tọa độ UTM WGS84 múi 48N bằng thuật toán đạc tam giác (Triangulation).
  • Hiệu chỉnh khí quyển: Áp dụng phương pháp loại bỏ điểm tối (Dark Object Subtraction - DOS) để triệt tiêu ảnh hưởng tán xạ khí dung và phân tử khí quyển, chiết xuất phổ phản xạ thực của bề mặt nước.
  • Mô hình hóa không gian: Sử dụng thuật toán Ordinary Kriging để nội suy giá trị độ đục trên toàn bộ mặt nước vịnh từ phương trình tương quan xây dựng được.

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

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

Quá trình phân tích số liệu thực nghiệm và giải đoán ảnh viễn thám đa thời gian đã chỉ ra các kết quả quan trọng sau:

Thứ nhất, thiết lập thành công phương trình tương quan định lượng giữa phổ phản xạ kênh 1 ảnh MODIS (dải sóng 0,620 - 0,670 µm) và giá trị độ đục thực tế. Dữ liệu thực nghiệm tại 40 điểm đo ngày 06/07/2010 cho thấy độ đục nước vịnh dao động từ 1,0 NTU đến 16,0 NTU, giá trị trung bình đạt 10,3 NTU, trung vị đạt 6,0 NTU và độ lệch chuẩn đạt khoảng 4,2 NTU. Mô hình hồi quy phi tuyến thể hiện sự đồng biến cao giữa hệ số phản xạ phổ và nồng độ chất rắn gây đục.

Thứ hai, xác định quy luật phân hóa không gian rõ rệt của độ đục nước vịnh Tiên Yên. Độ đục cao nhất tập trung tại các khu vực cửa sông Tiên Yên, cửa sông Ba Chẽ và cửa sông Hà Cối với giá trị đạt từ 12,0 NTU đến 16,0 NTU. Ngược lại, khu vực ngoài khơi và các luồng thông biển gần đảo Cái Bầu, đảo Cái Chiên, đảo Vĩnh Thực có độ đục thấp nhất, chỉ dao động từ 1,0 NTU đến 4,0 NTU.

Thứ ba, ghi nhận xu thế biến động độ đục theo mùa mưa lũ trong tháng 7 giai đoạn 2006 - 2011. Độ đục nước vịnh tăng mạnh vào các thời điểm có lượng mưa lớn. Năm 2011, với lượng mưa tháng 7 đạt 363,4 mm (và tháng 6 đạt 643,0 mm), diện tích vùng nước có độ đục trên 8,0 NTU mở rộng chiếm hơn 65% tổng diện tích mặt nước vịnh Tiên Yên, mở rộng phạm vi lan truyền ra phía các cửa luồng Cửa Mô và Cửa Đại.

Thứ tư, phương pháp địa thống kê Ordinary Kriging (OK) và Co-Kriging (CK) chứng minh độ tin cậy vượt trội trong việc thành lập bản đồ chuyên đề, với phương sai sai số ước tính duy trì ở mức tối thiểu, phản ánh chính xác ranh giới các khối nước có độ đục khác nhau.

Thảo luận kết quả

Cơ chế phân bố độ đục tại vịnh Tiên Yên chịu sự chi phối tổng hợp của các yếu tố động lực tự nhiên và hoạt động nhân sinh. Về tự nhiên, đặc thù địa hình lưu vực sông ngắn và dốc (sông Tiên Yên có độ dốc 28,1%, sông Hà Cối 22,5%) kết hợp chế độ nhật triều biên độ lớn (lên tới 4,0 m trong tháng 7) tạo nên động năng dòng triều rút mạnh trên 1,0 m/s tại các cửa sông, cuốn theo lượng lớn trầm tích bùn sét tầng mặt kích thước mịn. Vận tốc dòng chảy tại các cửa thoát triều như Cửa Mô đạt 74,2 cm/s và Cửa Đại đạt 49,3 cm/s làm gia tăng sự xáo trộn cơ học của lớp trầm tích đáy.

Về hoạt động nhân sinh, quá trình chuyển đổi 394,44 ha rừng ngập mặn giai đoạn 2000 - 2009 (trong đó 272,9 ha biến thành ao nuôi tôm) đã làm suy giảm nghiêm trọng khả năng giữ trầm tích tự nhiên của dải ven bờ. Hoạt động san lấp mặt bằng xây dựng khu công nghiệp cảng biển Hải Hà, cảng Đầm Buôn và việc xả thải trực tiếp từ 2.874 hộ nuôi thủy sản đã gia tăng tải lượng chất lơ lửng đổ ra vịnh.

Các chuỗi sơ đồ chuyên đề thành lập trong luận văn thể hiện trực quan diễn biến độ đục qua từng năm. Biểu đồ đường biến thiên phản ánh sự tương thích hoàn toàn giữa đỉnh độ đục viễn thám và đỉnh lượng mưa quan trắc tại trạm Quảng Hà. Khi so sánh với các nghiên cứu vùng biển hở ven bờ miền Trung, môi trường nước vịnh Tiên Yên có tính chất kín gió nhờ các đảo chắn ngoài khơi nên sự lan truyền trầm tích chủ yếu diễn ra theo hành lang luồng lạch sông - triều, giúp khẳng định tính ưu việt của ảnh viễn thám MODIS kênh 1 độ phân giải 250 m trong quan trắc các thủy vực ven biển nhiệt đới.

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

Nhằm nâng cao hiệu quả giám sát chất lượng nước và bảo vệ hệ sinh thái biển vịnh Tiên Yên, các giải pháp cụ thể được đề xuất như sau:

Thứ nhất, tích hợp quy trình giám sát viễn thám MODIS vào mạng lưới quan trắc môi trường định kỳ của tỉnh Quảng Ninh.

  • Hành động: Xây dựng cổng dữ liệu tự động tiếp nhận ảnh vệ tinh và chiết xuất bản đồ độ đục định kỳ 15 ngày/lần.
  • Mục tiêu: Giảm thiểu 40% chi phí vận hành quan trắc thực địa và nâng cao tần suất cập nhật dữ liệu chất lượng nước.
  • Thời gian thực hiện: Giai đoạn 2024 - 2026.
  • Chủ thể thực hiện: Sở Tài nguyên và Môi trường tỉnh Quảng Ninh phối hợp cùng Trung tâm Viễn thám Quốc gia.

Thứ hai, phục hồi và bảo tồn hệ sinh thái rừng ngập mặn và thảm cỏ biển ven vịnh.

  • Hành động: Trồng mới và phục hồi 200 ha rừng ngập mặn tại xã Đồng Rui, Đông Hải và bảo vệ nghiêm ngặt 150 ha cỏ biển tại vụng Hà Cối, 80 ha tại vụng Đầm Hà.
  • Mục tiêu: Giảm 20% lượng trầm tích xói lở trôi dạt vào lòng vịnh trong mùa mưa bão.
  • Thời gian thực hiện: Giai đoạn 2024 - 2028.
  • Chủ thể thực hiện: Chi cục Kiểm lâm Quảng Ninh, UBND các huyện Tiên Yên, Đầm Hà và Hải Hà.

Thứ ba, kiểm soát nghiêm ngặt nguồn thải từ hoạt động nuôi trồng thủy sản ven bờ.

  • Hành động: Bắt buộc 100% các trang trại nuôi tôm công nghiệp trên diện tích 2.900,7 ha phải đầu tư ao lắng lọc xử lý chất thải hữu cơ và bùn đáy trước khi xả ra môi trường biển.
  • Mục tiêu: Cắt giảm ít nhất 30% tải lượng chất rắn lơ lửng và chất hữu cơ thải vào các bãi triều gian triều.
  • Thời gian thực hiện: Hoàn thành trước quý IV năm 2025.
  • Chủ thể thực hiện: Phòng Nông nghiệp và Phát triển Nông thôn các huyện ven vịnh phối hợp với Công an Môi trường.

Thứ tư, áp dụng biện pháp giảm thiểu bùn cát tại các dự án nạo vét và san lấp ven biển.

  • Hành động: Lắp đặt màn chắn bùn lơ lửng (silt curtains) tại khu vực thi công nạo vét luồng tàu cảng Mũi Chùa, Cảng Đầm Buôn và khu công nghiệp Hải Hà.
  • Mục tiêu: Khống chế độ đục vùng lân cận công trình xây dựng không vượt quá ngưỡng 15 NTU.
  • Thời gian thực hiện: Áp dụng liên tục trong suốt tiến độ dự án từ 2024 đến 2030.
  • Chủ thể thực hiện: Ban Quản lý Khu kinh tế Quảng Ninh và các nhà thầu thi công hạ tầng cảng biển.

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

Nội dung và phương pháp luận của luận văn mang lại giá trị thực tiễn cao cho 4 nhóm đối tượng chính:

  1. Cơ quan quản lý nhà nước về tài nguyên và môi trường:
  • Lợi ích: Sở Tài nguyên và Môi trường, Chi cục Biển và Hải đảo tỉnh Quảng Ninh có thể ứng dụng trực tiếp thuật toán và quy trình GIS để thiết lập hệ thống cảnh báo sớm suy thoái môi trường nước mặt.
  • Use case: Tích hợp bản đồ độ đục vào quy hoạch không gian biển và đánh giá tác động môi trường của các dự án công nghiệp cảng biển.
  1. Các nhà nghiên cứu và viện chuyên ngành:
  • Lợi ích: Cung cấp cơ sở khoa học chi tiết về việc kết hợp dữ liệu vệ tinh đa thời gian MODIS với thuật toán địa thống kê Kriging trong điều kiện thủy vực phức tạp ven bờ Việt Nam.
  • Use case: Tham khảo phương pháp hiệu chỉnh khí quyển DOS và mô hình hóa lan truyền trầm tích cho các vùng vịnh khác như vịnh Cam Ranh hay đầm phá Tam Giang - Cầu Hai.
  1. Doanh nghiệp và hợp tác xã nuôi trồng thủy sản:
  • Lợi ích: Hiểu rõ quy luật phân bố dòng chảy và độ đục theo mùa triều để điều chỉnh lịch lấy nước, thả giống và thu hoạch an toàn.
  • Use case: Lựa chọn vị trí đặt lồng bè nuôi cá, bãi nuôi ngao tại các khu vực có độ đục ổn định dưới 6,0 NTU, tránh thiệt hại kinh tế do hiện tượng bồi lắng bùn.
  1. Giảng viên và học viên cao học, sinh viên đại học:
  • Lợi ích: Nguồn tài liệu học thuật hoàn chỉnh phục vụ nghiên cứu và giảng dạy các môn học chuyên đề Viễn thám ứng dụng, Hệ thống thông tin địa lý GIS và Quản lý môi trường biển.
  • Use case: Sử dụng bộ dữ liệu thực địa 40 điểm đo và chuỗi ảnh MODIS làm bài tập thực hành xử lý không gian trong ENVI và ArcGIS.

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

Kênh 1 ảnh vệ tinh MODIS có ưu điểm gì trong việc xác định độ đục nước biển? Kênh 1 của cảm biến MODIS hoạt động ở dải sóng 0,620 - 0,670 µm với độ phân giải không gian 250 m, là dải sóng ánh sáng đỏ có độ nhạy tối ưu nhất với hàm lượng chất rắn lơ lửng trong nước. Bước sóng này ít bị hấp thụ bởi các sắc tố diệp lục hòa tan so với dải sóng xanh, giúp phản ánh trung thực mức độ tán xạ của hạt trầm tích khoáng.

Phương pháp hiệu chỉnh khí quyển Dark Object Subtraction (DOS) đóng vai trò gì? Khí quyển và các hạt khí dung (aerosols) chiếm tới hơn 90% tổng tín hiệu quang học mà cảm biến vệ tinh thu nhận được trên bề mặt nước biển. Kỹ thuật DOS giúp loại bỏ phần bức xạ do khí quyển tán xạ ngược, khôi phục lại giá trị phổ phản xạ thực tế phát ra từ cột nước, từ đó bảo đảm độ chính xác của phương trình tính toán độ đục.

Vì sao độ đục nước vịnh Tiên Yên lại tăng cao đột biến trong tháng 7? Tháng 7 là thời điểm đỉnh điểm của mùa mưa bão tại Quảng Ninh với lượng mưa trung bình đạt trên 363 mm kết hợp chế độ triều cường cực đại lên tới 4,0 m. Lượng nước lũ đổ về từ 5 hệ thống sông chính cuốn theo một lượng lớn phù sa, đồng thời động năng dòng triều rút mạnh khuấy đảo lớp trầm tích bùn đáy khiến độ đục toàn vịnh tăng cao.

Mô hình địa thống kê Ordinary Kriging mang lại lợi ích gì so với phương pháp nội suy truyền thống? Khác với phương pháp nghịch đảo khoảng cách (IDW) chỉ tính toán dựa trên khoảng cách hình học, Ordinary Kriging sử dụng hàm biến sai bán phần (semi-variogram) để phân tích cấu trúc tương quan không gian thực tế và tính dị hướng của dòng nước. Điều này giúp ước tính chính xác giá trị tại các điểm không đo đạc và cung cấp bản đồ phương sai sai số tin cậy.

Dữ liệu ảnh MODIS độ phân giải 250 m có đáp ứng tốt cho các vịnh ven bờ diện tích nhỏ không? Đối với vịnh Tiên Yên có diện tích ngập nước 400 km2 và chiều dài 57 km, độ phân giải 250 m của MODIS kênh 1 cung cấp hàng nghìn điểm ảnh phủ kín mặt nước, hoàn toàn đáp ứng tốt việc phân tích biến động không gian tổng thể. Tuy nhiên, ở các lạch triều hẹp dưới 200 m, cần tích hợp thêm ảnh vệ tinh Landsat hoặc Sentinel để tăng độ nét chi tiết.

Kết luận

  • Luận văn đã chứng minh tính hiệu quả vượt trội của việc tích hợp công nghệ viễn thám MODIS độ phân giải 250 m và công cụ địa thống kê GIS trong quan trắc định lượng độ đục nước biển ven bờ.
  • Thiết lập thành công phương trình toán học chuyển đổi phổ phản xạ kênh 1 sang độ đục dựa trên 40 điểm đo thực nghiệm ngày 06/07/2010 với độ đục trung bình đạt 10,3 NTU.
  • Làm sáng tỏ quy luật phân bố độ đục không gian và chu kỳ biến động tháng 7 các năm từ 2006 đến 2011, khẳng định sự tác động đồng thời của dòng chảy sông - triều và sự suy giảm 394,44 ha rừng ngập mặn.
  • Cung cấp cơ sở khoa học thực tiễn giúp giảm thiểu hơn 40% chi phí quan trắc môi trường truyền thống cho các cơ quan quản lý tại tỉnh Quảng Ninh giai đoạn 2024 - 2030.
  • Các nhà quản lý và đơn vị nghiên cứu quan tâm có thể áp dụng ngay quy trình kỹ thuật này để tối ưu hóa công tác giám sát tài nguyên biển và bảo tồn đa dạng sinh học cho các vùng vịnh ven bờ trên toàn quốc.