Tổng quan nghiên cứu

Khu vực cửa sông ven biển là môi trường động lực phức tạp, nơi diễn ra quá trình tương tác mạnh mẽ giữa lục địa và đại dương. Cửa Đáy đóng vai trò là một trong những cửa sông quan trọng ở phía Bắc Việt Nam, tiếp nhận lượng bùn cát từ sông Đáy khoảng 34.000 tấn mỗi năm cùng lượng bùn cát dịch chuyển dọc bờ từ phía Bắc xuống lên tới 220.000 tấn mỗi năm. Sự biến động của tổng chất rắn lơ lửng (TSS) không chỉ chi phối trực tiếp đến chế độ bồi tụ, xói lở mà còn tác động sâu sắc đến hệ sinh thái và năng suất sinh học của hơn 6.058 ha mặt nước nuôi trồng thủy sản trong khu vực.

Tuy nhiên, việc quan trắc TSS bằng các phương pháp đo đạc truyền thống trên phạm vi bờ biển dài hơn 30 km thường tốn kém chi phí, đòi hỏi nhiều nhân lực và khó nắm bắt được bức tranh phân bố liên tục theo không gian và thời gian. Đề tài luận văn thạc sĩ khoa học chuyên ngành Địa chất Môi trường của tác giả Phạm Thị Nhung, dưới sự hướng dẫn của TS. Nguyễn Thị Thu Hà tại Trường Đại học Khoa học Tự nhiên - Đại học Quốc gia Hà Nội, đã tập trung giải quyết bài toán cấp thiết này.

Mục tiêu cụ thể của nghiên cứu bao gồm: xác định mối tương quan định lượng giữa hàm lượng TSS và phổ phản xạ kênh 4 (dải sóng đỏ) của vệ tinh Landsat 8; mô hình hóa sự phân bố không gian - thời gian của TSS vùng Cửa Đáy qua các mùa; và phân tích quá trình vận chuyển trầm tích từ sông Đáy ra biển Đông trong giai đoạn 2013 - 2015. Kết quả nghiên cứu mang lại ý nghĩa khoa học và thực tiễn cao, giúp giảm thiểu khoảng 50% đến 60% chi phí quan trắc thực địa, đồng thời cung cấp công cụ giám sát môi trường nước mặt với độ phân giải không gian 30 m chính xác và kịp thời.

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 truyền xạ quang học môi trường nước (Case-2 water optical theory), lý thuyết tích hợp công nghệ địa không gian (Geospatial Technology bao gồm Viễn thám và GIS), và lý thuyết địa thống kê không gian (Spatial Geostatistics). Trong môi trường nước cửa sông, các phân tử nước, chất hữu cơ hòa tan (CDOM), diệp lục (Chlorophyll-a) và các phần tử chất rắn lơ lửng hấp thụ cũng như tán xạ bức xạ mặt trời theo các quy luật quang phổ đặc thù.

Ba khái niệm trung tâm được chuẩn hóa xuyên suốt nghiên cứu bao gồm:

  1. Tổng chất rắn lơ lửng (TSS): Tập hợp các hạt vô cơ (ISS) và hữu cơ (OMC) không hòa tan có kích thước lớn hơn 0,45 micromet lơ lửng trong cột nước.
  2. Phổ phản xạ rời mặt nước: Tỷ số giữa độ chói rời mặt nước và độ chiếu sáng tới, biểu thị đặc trưng quang học riêng biệt của thủy vực tại từng bước sóng.
  3. Địa thống kê không gian (Ordinary Kriging): Phương pháp nội suy tuyến tính không thiên vị tối ưu hóa dựa trên hàm tương quan không gian bán phương sai (semi-variogram).

Mô hình nghiên cứu kế thừa các phát hiện kinh điển của Ritchie và cộng sự năm 1976 cùng các nghiên cứu ven biển châu thổ sông Hồng của Nguyễn Văn Thảo năm 2016, xác định rằng trong dải phổ ánh sáng nhìn thấy từ 600 nm đến 800 nm, đặc biệt tại bước sóng đỏ khoảng 655 nm, phổ phản xạ có mức độ nhạy cảm cao nhất và tỷ lệ thuận chặt chẽ với hàm lượng bùn cát lơ lửng.

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

Nguồn dữ liệu của luận văn được tích hợp từ chuỗi dữ liệu thực nghiệm mặt đất và tư liệu viễn thám đa thời gian:

  • Dữ liệu viễn thám: Sử dụng 9 cảnh ảnh vệ tinh Landsat 8 OLI (Path/Row: 126/46) thu nhận từ Cục Khảo sát Địa chất Hoa Kỳ giai đoạn 2013 - 2015 với độ che phủ mây dưới 20%, kết hợp dữ liệu ảnh Sentinel-2A đa phổ độ phân giải 10 m nhằm thử nghiệm nâng cao chi tiết không gian.
  • Cỡ mẫu và phương pháp chọn mẫu thực địa: Tổng cộng 24 mẫu nước mặt được thu thập qua 2 đợt khảo sát đại diện (12 mẫu ngày 30/12/2014 thuộc cao điểm mùa khô và 12 mẫu ngày 10/07/2015 thuộc mùa mưa lũ). Mẫu được lấy đồng thời tại thời điểm vệ tinh quét qua, phân bố theo tuyến mạng lưới mặt cắt ngang dọc từ bến phà Quỹ Nhất chạy dài ra cửa sông đến khu vực đảo Cồn Nổi ở độ sâu 10 cm cách mặt nước.

Quy trình phân tích phòng thí nghiệm tuân thủ tiêu chuẩn APHA 1998 bằng phương pháp sấy cân màng lọc sợi thủy tinh Whatman GF/F tại nhiệt độ 105 độ C và nung ở 500 độ C để bóc tách thành phần vô cơ và hữu cơ. Dữ liệu ảnh được hiệu chỉnh bức xạ tuyệt đối và hiệu chỉnh khí quyển theo thuật toán FLAASH trên phần mềm ENVI 5.2, giúp loại bỏ trên 90% nhiễu tán xạ khí quyển. Sau đó, mô hình hồi quy phi tuyến được kiểm định bằng phần mềm IBM SPSS Statistics 20 ở độ tin cậy 95%, và bản đồ phân bố trường TSS được nội suy tự động bằng mô hình Ordinary Kriging trong ArcGIS 10.2. Lý do lựa chọn Ordinary Kriging là khả năng tính toán trọng số không gian tối ưu, giảm thiểu phương sai sai số ước tính và phản ánh chính xác tính dị hướng của dòng thủy triều ven bờ.

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 chi tiết về đặc trưng phân bố vật chất lơ lửng tại khu vực Cửa Đáy:

Thứ nhất, phân tích mẫu thực địa cho thấy sự phân hóa sâu sắc về hàm lượng TSS giữa các mùa khí hậu. Trong đợt khảo sát mùa khô (tháng 12/2014), hàm lượng TSS dao động từ 9,1 mg/l đến 246,0 mg/l với giá trị trung bình là 68,8 mg/l và độ lệch chuẩn đạt 72,2. Ngược lại, vào mùa mưa lũ (tháng 07/2015), hàm lượng TSS tăng đột biến, dao động từ 12,2 mg/l đến 612,3 mg/l với giá trị trung bình đạt 174,6 mg/l và độ lệch chuẩn là 184,5. Như vậy, nồng độ chất rắn lơ lửng trung bình trong mùa mưa cao gấp 2,54 lần (tăng hơn 153,7%) so với mùa khô.

Thứ hai, nghiên cứu đã thiết lập thành công phương trình hồi quy hàm mũ mô tả mối quan hệ giữa hàm lượng TSS và phổ phản xạ kênh 4 (bước sóng đỏ 655 nm) của ảnh Landsat 8 sau khi xử lý qua mô đun FLAASH: Hàm lượng TSS tính toán theo công thức hàm mũ với hệ số tương quan xác định đạt mức cao (R bình phương = 0,79). Kết quả đo phổ thực địa tại cửa Ba Lạt lân cận cũng khẳng định độ tương thích giữa phổ đo ngoài trời và phổ ảnh vệ tinh đạt hệ số tương quan R = 0,86.

Thứ ba, trường phân bố không gian của TSS cho thấy quy luật suy giảm rõ rệt theo hướng từ sông ra biển khơi. Nồng độ TSS cao nhất tập trung tại lòng dẫn cửa sông và khu vực cồn cát ngập triều trước cửa Đáy (thường xuyên vượt ngưỡng 300 - 500 mg/l vào mùa mưa), sau đó giảm nhanh xuống dưới 30 mg/l khi ra đến đới biển khơi cách bờ từ 5 km đến 10 km.

Thảo luận kết quả

Sự biến động mạnh mẽ của hàm lượng TSS tại Cửa Đáy phản ánh sự tương tác phức tạp giữa chế độ thủy thạch động lực và các yếu tố khí tượng hải văn. Vào mùa lũ (từ tháng 6 đến tháng 10), lượng nước từ thượng nguồn sông Đáy đổ ra chiếm từ 75% đến 85% tổng lượng dòng chảy cả năm với lưu lượng trung bình khoảng 195 mét khối trên giây, mang theo khối lượng phù sa mịn khổng lồ kích thước trung bình 0,17 mm. Dưới tác động của sóng hướng Đông Nam và Nam có chiều cao sóng biển ven bờ đạt từ 4 m đến 5 m trong các đợt áp thấp nhiệt đới, trầm tích đáy bị khuấy đục mạnh mẽ tạo nên các vệt phù sa lơ lửng trải rộng ra phía Đông Nam.

Vào mùa khô (từ tháng 11 đến tháng 4 năm sau), dòng chảy sông suy giảm chỉ còn chiếm 7% đến 8% tổng lượng dòng chảy năm, tạo điều kiện cho triều cường và nêm mặn lấn sâu vào nội địa. Gió mùa Đông Bắc thịnh hành với tần suất 60% đến 70% tạo ra dòng chảy ven bờ theo hướng Bắc - Nam, kết hợp độ cao sóng trung bình từ 0,1 m đến 0,4 m, giới hạn phạm vi phát tán phù sa trong một dải hẹp sát bờ phía huyện Nghĩa Hưng.

Khi so sánh với các nghiên cứu trước đây của Nguyễn Quốc Phi năm 2014 và Nguyễn Văn Thảo năm 2016, mô hình hàm mũ của luận văn thể hiện tính ưu việt vượt trội so với mô hình hồi quy đa biến tuyến tính nhờ phản ánh đúng bản chất vật lý quang học nước khi nồng độ phù sa tăng cao. Dữ liệu này có thể được trình bày trực quan thông qua bản đồ viễn thám phân vùng nồng độ đa thời gian kết hợp biểu đồ biến thiên mặt cắt ngang thủy lực, giúp các nhà quản lý dễ dàng nhận diện ranh giới các vùng xói lở và bồi tụ cục bộ.

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

Dựa trên kết quả nghiên cứu thực nghiệm và mô hình hóa không gian, luận văn đưa ra 4 nhóm giải pháp chiến lược nhằm tối ưu hóa công tác quản lý môi trường nước và quy hoạch vùng bờ Cửa Đáy:

  1. Thiết lập hệ thống giám sát chất lượng nước mặt tự động bằng viễn thám định kỳ: Sở Tài nguyên và Môi trường tỉnh Ninh Bình cùng tỉnh Nam Định cần chủ trì triển khai hệ thống giải đoán ảnh vệ tinh Landsat 8 và Sentinel-2A với tần suất 5 đến 10 ngày một lần. Mục tiêu cụ thể là kiểm soát tự động biến động TSS trên diện tích hơn 200 km2 vùng cửa sông, giảm 50% kinh phí đo đạc thủ công trong giai đoạn 2026 - 2028.

  2. Tối ưu hóa quy hoạch nạo vét luồng hàng hải Cửa Đáy: Ban Quản lý Cảng Ninh Phúc và Cục Hàng hải Việt Nam cần sử dụng bản đồ phân bố bùn cát lắng đọng để xác định chính xác các điểm nghẽn sa bồi trước cửa sông. Mục tiêu đảm bảo cho tàu trọng tải từ 3.000 DWT đến 5.000 DWT lưu thông an toàn vào hệ thống cảng nội địa theo Nghị định 04/2011/NĐ-CP, đồng thời hạn chế tình trạng sạt lở đê kè ven biển giai đoạn 2027 - 2030.

  3. Quy hoạch vùng nuôi trồng thủy sản an toàn thích ứng biến đổi môi trường: Chi cục Thủy sản hai tỉnh Ninh Bình và Nam Định cần tích hợp dữ liệu biến động TSS vào quy hoạch vùng nuôi ngao và tôm nước lợ tại các xã Kim Trung, Kim Đông và Nam Điền. Mục tiêu duy trì ngưỡng TSS trong nước cấp dưới 100 mg/l cho 6.058 ha đầm nuôi, nâng sản lượng thủy sản toàn vùng vượt mức 45.000 tấn mỗi năm từ nay đến năm 2028.

  4. Mở rộng hệ thống đai rừng ngập mặn phòng hộ và kè chắn sóng sinh thái: Ủy ban nhân dân các huyện Kim Sơn và Nghĩa Hưng cần đẩy mạnh trồng mới khoảng 300 ha đến 500 ha rừng ngập mặn tại bãi bồi ngoài đê Bình Minh-3 và khu vực Cồn Nổi trước năm 2030. Giải pháp này giúp giữ lại 30% lượng bùn cát bồi tụ tự nhiên, bảo vệ vùng bờ trước các cơn bão có gió giật trên cấp 12 và hiện tượng nước biển dâng.

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

Công trình nghiên cứu mang lại giá trị ứng dụng thiết thực cho 4 nhóm đối tượng chính:

  1. Cơ quan quản lý nhà nước về tài nguyên, môi trường và đê điều: Sở Tài nguyên và Môi trường, Chi cục Biển và Hải đảo, Ban Chỉ huy Phòng chống thiên tai các tỉnh ven biển đồng bằng Bắc Bộ có thể sử dụng phương pháp luận và quy trình xử lý ảnh trong luận văn để xây dựng bản đồ quan trắc chất lượng nước và đánh giá rủi ro bồi lấp luồng tàu.

  2. Ban quản lý cảng biển và các đơn vị giao thông vận tải đường thủy: Các đơn vị quản lý Cảng Ninh Phúc, Cảng Ninh Bình và doanh nghiệp khai thác luồng hàng hải áp dụng kết quả phân tích trường dòng chảy bùn cát để lựa chọn thời điểm và vị trí nạo vét luồng tối ưu, tiết kiệm hàng tỷ đồng kinh phí nạo vét định kỳ.

  3. Hợp tác xã và doanh nghiệp nuôi trồng thủy sản ven biển: Các hộ nuôi tôm sú, ngao giống tại Kim Sơn và Nghĩa Hưng có thể ứng dụng thông tin về chu kỳ phát tán TSS theo mùa để chủ động điều tiết lịch lấy nước phù sa giàu dinh dưỡng, phòng tránh hiện tượng phù sa quá đậm đặc làm tắc nghẽn hô hấp của thủy sản.

  4. Giảng viên, nghiên cứu sinh và sinh viên chuyên ngành Khoa học Trái đất: Các nhà nghiên cứu thuộc ngành Địa chất Môi trường, Hải dương học, Viễn thám và GIS sử dụng tài liệu như một nguồn tham khảo chuẩn mực về quy trình hiệu chỉnh khí quyển FLAASH, tích hợp thống kê địa lý Kriging và thuật toán bán thực nghiệm trong thủy vực nước nông ven bờ.

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

  1. Tại sao kênh phổ đỏ (kênh 4) của Landsat 8 lại tối ưu nhất để xác định hàm lượng TSS? Kênh 4 của Landsat 8 hoạt động ở dải bước sóng 640 nm đến 670 nm (trung tâm 655 nm). Trong dải sóng này, nước hấp thụ ánh sáng mạnh nhưng các hạt vô cơ và phù sa lơ lửng lại có độ tán xạ cực đại. Thực nghiệm tại Cửa Đáy cho thấy độ phản xạ kênh 4 tương quan chặt chẽ với TSS theo hàm mũ đạt hệ số xác định R bình phương bằng 0,79, vượt trội hoàn toàn so với các kênh xanh lam hay cận hồng ngoại.

  2. Hiệu chỉnh khí quyển FLAASH đóng vai trò quan trọng như thế nào trong nghiên cứu? Bức xạ do các hạt sol khí và phân tử không khí tán xạ có thể chiếm hơn 90% tổng tín hiệu mà cảm biến vệ tinh thu nhận được trên bề mặt nước. Mô đun FLAASH trên phần mềm ENVI 5.2 sử dụng mô hình truyền xạ MODTRAN để loại bỏ triệt để nhiễu khí quyển, giúp phổ phản xạ bề mặt đạt độ tương quan cao R = 0,86 so với dữ liệu đo phổ trực tiếp ngoài thực địa.

  3. Sự phân bố hàm lượng TSS giữa mùa mưa và mùa khô tại Cửa Đáy khác nhau như thế nào? Vào mùa mưa lũ (tháng 7), lượng phù sa từ sông Đáy đổ ra rất lớn khiến TSS trung bình đạt 174,6 mg/l (cực đại 612,3 mg/l), lan tỏa rộng ra biển theo hướng Đông Nam. Vào mùa khô (tháng 12), dòng chảy sông suy kiệt khiến TSS giảm xuống mức trung bình 68,8 mg/l (cực đại 246,0 mg/l) và bị ép sát vào dải ven bờ dưới tác động của dòng chảy gió mùa Đông Bắc.

  4. Dữ liệu vệ tinh Sentinel-2A mang lại ưu thế gì khi mở rộng nghiên cứu giám sát TSS? Sentinel-2A cung cấp độ phân giải không gian vượt trội 10 m ở các kênh ánh sáng nhìn thấy so với 30 m của Landsat 8, cùng chu kỳ lặp ngắn từ 5 ngày khi kết hợp Sentinel-2B. Điều này cho phép bóc tách chi tiết các dải phù sa hẹp trong các lạch triều tự nhiên và hệ thống kênh dẫn thủy sản tại Nghĩa Hưng và Kim Sơn với độ chính xác hình học cao hơn.

  5. Phương pháp địa thống kê Ordinary Kriging đem lại lợi ích gì cho bản đồ phân bố chất lượng nước? Ordinary Kriging cho phép mô hình hóa cấu trúc không gian của TSS thông qua hàm bán phương sai semi-variogram. Phương pháp này thực hiện ước tính tuyến tính không thiên vị tại các điểm không có mẫu đo, giảm thiểu tối đa phương sai sai số và tạo ra bề mặt bản đồ liên tục, mượt mà, phản ánh chân thực các ranh giới thủy động lực ven bờ.

Kết luận

  • Luận văn đã xác lập thành công mô hình toán học hàm mũ liên kết giữa phổ phản xạ kênh 4 ảnh vệ tinh Landsat 8 với hàm lượng TSS vùng Cửa Đáy với độ tin cậy cao (R bình phương = 0,79).
  • Định lượng rõ nét biên độ dao động nồng độ TSS theo mùa, từ trung bình 68,8 mg/l trong mùa khô lên tới 174,6 mg/l trong mùa mưa, với giá trị cao nhất đạt 612,3 mg/l tại cửa sông.
  • Làm sáng tỏ cơ chế vận chuyển phù sa chịu sự chi phối của dòng chảy sông Đáy (34.000 tấn/năm) vào mùa mưa và dòng dọc bờ Bắc - Nam (220.000 tấn/năm) vào mùa khô.
  • Chứng minh tính khả thi vượt trội của việc tích hợp công nghệ viễn thám quang học độ phân giải cao và mô hình nội suy địa thống kê Ordinary Kriging trong giám sát môi trường nước mặt ven biển.
  • Đề xuất hệ thống giải pháp liên ngành phục vụ hiệu quả công tác nạo vét luồng tàu 3.000 - 5.000 DWT tại Cảng Ninh Phúc và bảo vệ hệ sinh thái 6.058 ha nuôi trồng thủy sản.

Trong giai đoạn 2026 - 2030, việc mở rộng tích hợp nguồn ảnh Sentinel-2 và xây dựng nền tảng viễn thám trực tuyến thời gian thực sẽ là bước đi tất yếu. Các nhà khoa học, cơ quan quản lý và doanh nghiệp ven biển hãy chủ động ứng dụng công nghệ địa không gian tiên tiến này để thúc đẩy phát triển kinh tế biển bền vững và bảo vệ tài nguyên môi trường hiệu quả.