Tổng quan nghiên cứu

Khu vực phía Nam đồng bằng sông Hồng, đặc biệt là tỉnh Nam Định, đóng vai trò trọng điểm về an ninh lương thực của vùng Bắc Bộ với diện tích đất nông nghiệp chiếm tới 65% tổng diện tích tự nhiên. Hơn 87,6% dân số địa phương sinh sống dựa vào sản xuất nông nghiệp, đóng góp khoảng 40% vào tổng sản phẩm trên địa bàn (GDP) toàn tỉnh. Tuy nhiên, dải bờ biển dài 72 km cùng hệ thống sông ngòi phân bố dày đặc đang khiến địa phương này trở thành tâm điểm dễ bị tổn thương nhất trước các tác động của biến đổi khí hậu. Theo số liệu quan trắc dài hạn, mực nước biển tại Nam Định dâng trung bình 2,15 mm mỗi năm, đồng thời đường bờ biển bị xói lở và lấn sâu vào đất liền trung bình khoảng 10 m mỗi năm. Hiện tượng biến đổi dòng chảy thượng nguồn kết hợp triều cường cực đoan đã đẩy nhanh quá trình xâm nhập mặn trên hai trục thủy văn huyết mạch là sông Đáy và sông Ninh Cơ, khiến khoảng 10% đến 20% diện tích đất nông nghiệp vụ xuân thường xuyên đối mặt với nguy cơ hạn hán và thiếu nước tưới trầm trọng.

Nghiên cứu được thực hiện nhằm mục tiêu đánh giá định lượng hiện trạng xâm nhập mặn, phân tích quy luật diễn biến theo thời gian và không gian, đồng thời ứng dụng công nghệ mô hình số thủy động lực để dự tính xu thế nhiễm mặn trong tương lai theo các kịch bản biến đổi khí hậu giai đoạn đến năm 2100. Phạm vi không gian nghiên cứu tập trung vào lưu vực hạ lưu sông Đáy và sông Ninh Cơ thuộc địa phận tỉnh Nam Định. Công trình có ý nghĩa khoa học và thực tiễn sâu sắc khi cung cấp hệ thống cơ sở dữ liệu số hóa tin cậy về ranh giới mặn, hỗ trợ đắc lực cho các nhà quản lý trong việc hoạch định chính sách thích ứng, chuyển đổi cơ cấu cây trồng và bảo vệ nguồn tài nguyên đất, nước phục vụ phát triển nông nghiệp bền vữ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 được xây dựng dựa trên nền tảng thủy động lực học dòng chảy một chiều và lý thuyết lan truyền vật chất trong môi trường nước hở. Cốt lõi của khung lý thuyết bao gồm hệ phương trình Saint-Venant mô tả dòng chảy không ổn định trong mạng lưới kênh sông, kết hợp giữa phương trình liên tục (bảo toàn khối lượng) và phương trình động lượng (bảo toàn động lượng có xét đến gia tốc trọng trường 9,81 m/s² cùng hệ số nhám lòng dẫn Manning n). Quá trình di chuyển và phân tán của muối trong dòng chảy được mô hình hóa thông qua phương trình truyền tải - khuếch tán một chiều Fick, phản ánh đồng thời hai cơ chế: vận chuyển bình lưu theo hướng dòng chảy chính và phân tán khuếch tán do gradien chênh lệch nồng độ chất hòa tan.

Về mặt tiêu chuẩn đánh giá nông nghiệp, luận văn áp dụng khung phân loại chất lượng nước tưới của Tổ chức Nông Lương Liên Hợp Quốc (FAO). Theo đó, nước có nồng độ muối dưới 0,5‰ được xếp vào nhóm không mặn (dùng cho sinh hoạt và tưới tiêu tối ưu), từ 0,5‰ đến 1,5‰ là nước mặn ít (vẫn có thể tưới cho đa số cây trồng thông thường), từ 1,5‰ đến 7‰ là mặn vừa, và trên 7‰ là nước mặn đến rất mặn gây độc hại cho thực vật. Trong điều kiện canh tác lúa nước tại đồng bằng sông Hồng, ngưỡng nồng độ 1,5‰ được chọn làm giới hạn an toàn tối đa cho việc lấy nước tưới; khi độ mặn vượt quá 1,5‰, áp suất thẩm thấu tăng cao sẽ ức chế rễ hấp thu dưỡng chất, phá hủy cấu trúc keo đất do thừa ion Na+, khiến cây trồng suy thoái và giảm năng suất nghiêm trọng. Ngoài ra, nghiên cứu cũng tham chiếu quy chuẩn của Trung tâm Phòng tránh và Giảm nhẹ Thiên tai với ngưỡng mặn 4‰ để nhận diện vùng ngập mặn cực đoan.

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

Luận văn triển khai đồng bộ ba hướng tiếp cận cốt lõi: tiếp cận hệ thống, tiếp cận lịch sử và tiếp cận liên vùng - liên ngành. Phương pháp nghiên cứu kết hợp giữa khảo sát thực địa, phân tích hóa nghiệm phòng thí nghiệm và mô phỏng số học nâng cao:

  • Khảo sát thực địa và lấy mẫu: Cỡ mẫu khảo sát gồm 48 mẫu nước mặt được thu thập tại 16 vị trí trắc dọc chiến lược trên sông Đáy và sông Ninh Cơ vào tháng 3 năm 2016 (giai đoạn đỉnh điểm mùa kiệt). Phương pháp chọn mẫu tuân thủ quy chuẩn đo đạc thủy văn: tại mỗi điểm lấy mẫu ở 3 tầng độ sâu khác nhau gồm 0,2H, 0,5H và 0,8H (với H là độ sâu lòng sông) nhằm tính toán độ mặn trung bình mặt cắt ngang theo công thức bình quân gia quyền.
  • Phân tích hóa nghiệm: Độ mặn được xác định trong phòng thí nghiệm thông qua phương pháp chuẩn độ Bạc Nitrat (AgNO3) kinh điển của Knudsen với chất chỉ thị màu Kali Cromat (K2CrO4). Phương pháp này cho phép định lượng chính xác tổng hàm lượng ion halogen (chủ yếu là Cl-, chiếm 55,61% tổng lượng khoáng hòa tan) trong môi trường dung dịch có độ pH từ 7,0 đến 8,5, đảm bảo sai số ở mức tối thiểu.
  • Mô hình hóa số trị: Bộ phần mềm MIKE 11 do Viện Thủy lực Đan Mạch (DHI) phát triển được lựa chọn làm công cụ tính toán chính, tích hợp đồng thời module thủy lực HD và module truyền tải khuếch tán AD. Cơ sở dữ liệu đầu vào bao gồm 59 mặt cắt địa hình thực đo dọc hai triền sông, kết hợp chuỗi số liệu thủy văn lưu lượng thượng nguồn Q(t), mực nước triều hạ lưu H(t) và độ mặn biên được ghi nhận liên tục với tần suất 1 giờ/lần trong suốt tháng 1 năm 2010 do Sở Tài nguyên và Môi trường tỉnh Nam Định cung cấp. Mô hình sau đó được hiệu chỉnh và kiểm định độc lập với chuỗi số liệu thực đo tháng 3 năm 2016, đem lại độ tin cậy và mức độ khớp nối rất cao.

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

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

Kết quả phân tích từ chuỗi số liệu lịch sử, thực nghiệm đo đạc năm 2016 và mô phỏng MIKE 11 đã chỉ ra những đặc trưng cơ bản về quy luật xâm nhập mặn tại phía Nam đồng bằng sông Hồng:

  • Tính quy luật theo thời gian: Xâm nhập mặn biến thiên sâu sắc theo mùa và chu kỳ nhật triều. Độ mặn cao nhất trong năm luôn tập trung vào các tháng mùa kiệt (tháng 1, tháng 2 và tháng 3). Đỉnh mặn trên sông thường xuất hiện trễ hơn đỉnh triều ngoài cửa biển khoảng 1 giờ, trong khi chân mặn xuất hiện gần như đồng thời với chân triều.
  • Sự phân hóa không gian giữa hai hệ thống sông: Sông Ninh Cơ có mức độ và phạm vi nhiễm mặn khốc liệt hơn hẳn sông Đáy. Điển hình vào mùa kiệt năm 2007, ranh mặn 1‰ đã ăn sâu tới 46 km trên sông Ninh Cơ và 43 km trên sông Đáy; năm 2006, cự ly này lần lượt là 38,5 km (Ninh Cơ) và 26,5 km (Đáy).
  • Bằng chứng từ số liệu đo thực địa tháng 3 năm 2016: Tại sông Ninh Cơ, nồng độ mặn ghi nhận ở vị trí cách biển 1 km đạt tới 27,55‰, cách 10,7 km đạt 19,25‰, cách 19,8 km đạt 5,70‰ và chỉ giảm về mức an toàn 0,50‰ khi đến vị trí cách biển 27,3 km. Trên sông Đáy, nồng độ mặn tại điểm cách biển 2,8 km là 28,54‰, cách 11 km là 22,65‰, cách 19,6 km là 16,22‰ và giảm xuống 0,37‰ ở vị trí cách biển 25,6 km.
  • Xu thế gia tăng dưới tác động của biến đổi khí hậu: Khi chạy mô hình dự tính theo Kịch bản Biến đổi Khí hậu và Nước biển dâng năm 2016 của Bộ Tài nguyên và Môi trường, dưới kịch bản phát thải trung bình RCP 4.5 (nước biển dâng 13 cm vào năm 2030, 22 cm vào năm 2050 và 53 cm vào năm 2100) cũng như kịch bản cao RCP 8.5 (nước biển dâng 25 cm vào năm 2050 và 72 cm vào năm 2100), ranh giới mặn 1,5‰ trên cả hai dòng sông đều tịnh tiến sâu thêm vào nội đồng từ 3 đến 7 km. Đáng chú ý, nếu mực nước biển dâng 50 cm, ước tính có tới 26% diện tích tự nhiên của tỉnh Nam Định sẽ rơi vào vùng nguy cơ ngập úng và nhiễm mặn.

Thảo luận kết quả

Sự khác biệt về mức độ xâm nhập mặn giữa sông Ninh Cơ và sông Đáy bắt nguồn từ mối tương tác phức tạp giữa chế độ nhật triều đều có biên độ triều cực đại lên tới 3,31 m tại vịnh Bắc Bộ với đặc điểm hình thái đáy sông và lưu lượng nước ngọt thượng nguồn đổ về. Dòng chảy sông Đáy nhận thêm một phần nước hồi lưu và có hệ thống phân lưu khác biệt, trong khi sông Ninh Cơ có lòng dẫn dốc thoải, mặt cắt mở rộng dần ra phía cửa Lạch Giang, tạo điều kiện thuận lợi cho sóng triều truyền nhanh với vận tốc 15 đến 20 km/h, đẩy nêm mặn đi thẳng vào sâu nội địa.

Về phương diện biểu diễn trực quan, các kết quả mô phỏng đã được trực quan hóa thông qua hệ thống biểu đồ đường quá trình độ mặn theo thời gian thực và biểu đồ trắc dọc không gian biểu diễn sự suy giảm nồng độ muối từ cửa sông vào nội đồng. Các bảng thống kê cự ly xâm nhập mặn cực đại tương ứng với các ngưỡng 1‰, 1,5‰, 4‰ và 7‰ cho thấy rõ nguy cơ các công trình cống lấy nước tưới nông nghiệp tại các huyện Giao Thủy, Hải Hậu, Nghĩa Hưng sẽ bị vô hiệu hóa trong những đợt triều cường nếu không có giải pháp kiểm soát chủ động.

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

Dựa trên kết quả dự tính diễn biến mặn và đánh giá rủi ro nông nghiệp, luận văn đề xuất 4 nhóm giải pháp chiến lược mang tính hành động cao:

  • Thiết lập mạng lưới quan trắc và cảnh báo xâm nhập mặn tự động: Sở Tài nguyên và Môi trường tỉnh Nam Định cần chủ trì triển khai lắp đặt hệ thống cảm biến đo mặn tự động theo công nghệ SCADA tại 100% các cống đầu mối và trạm bơm chính dọc sông Đáy và sông Ninh Cơ trong giai đoạn 2018 đến 2025. Mục tiêu hướng tới việc cập nhật dữ liệu độ mặn theo thời gian thực lên cổng thông tin điện tử, giúp các hợp tác xã nông nghiệp chủ động đóng mở cống lấy nước ngọt khi độ mặn dưới 1,5‰.
  • Nâng cấp và kiên cố hóa hệ thống công trình thủy lợi ngăn mặn: Sở Nông nghiệp và Phát triển nông thôn cần ưu tiên phân bổ ngân sách trong kế hoạch trung hạn 2021-2030 để nâng cấp các tuyến đê sông, đê biển chịu ăn mòn muối, đồng thời xây dựng mới các cống ngăn mặn tự động theo kết cấu đập trụ hiện đại (tương tự mô hình đập Ba Lai). Giải pháp này nhằm mục tiêu kiểm soát triệt để ranh mặn 1,5‰, bảo vệ an toàn nguồn nước tưới cho hơn 70.000 ha đất trồng lúa.
  • Chuyển đổi cơ cấu giống và lịch thời vụ thích ứng mặn: Ủy ban nhân dân các huyện ven biển gồm Giao Thủy, Hải Hậu và Nghĩa Hưng cần phối hợp với các viện nghiên cứu nông nghiệp chuyển đổi từ 15% đến 20% diện tích đất trũng nhiễm mặn thường xuyên sang mô hình luân canh lúa - tôm hoặc trồng các giống lúa chịu mặn (ngưỡng chịu đựng từ 2‰ đến 3‰). Lộ trình thực hiện hoàn thành trước năm 2025 nhằm đảm bảo sinh kế bền vững cho người dân.
  • Điều tiết liên hồ chứa thủy điện thượng nguồn sông Hồng: Đề nghị Bộ Công Thương và Bộ Nông nghiệp và Phát triển nông thôn phối hợp chỉ đạo xả nước gia tăng từ các hồ chứa thủy điện lớn (Hòa Bình, Tuyên Quang, Thác Bà) trong các đợt cao điểm mùa kiệt từ tháng 1 đến tháng 3 hàng năm. Việc duy trì lưu lượng xả đáy tối thiểu đủ lớn sẽ tạo dòng chảy đẩy nêm mặn ra xa cửa biển từ 5 đến 10 km, tạo điều kiện thuận lợi cho hạ lưu lấy nước phục vụ gieo cấy vụ đông xuân.

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

Công trình nghiên cứu mang giá trị chuyên khảo sâu sắc, phù hợp cho 4 nhóm đối tượng chính:

  • Cơ quan quản lý nhà nước và hoạch định chính sách: 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ại Nam Định và các tỉnh ven biển đồng bằng Bắc Bộ có thể khai thác các kịch bản dự tính để xây dựng Kế hoạch hành động ứng phó với biến đổi khí hậu và quy hoạch vùng sản xuất nông nghiệp đến năm 2030, tầm nhìn 2050.
  • Kỹ sư công trình và chuyên gia quy hoạch thủy lợi: Đội ngũ tư vấn thiết kế hạ tầng thủy lợi có thể sử dụng các thông số thủy lực, mặt cắt sông và phân bố ranh mặn trong luận văn làm cơ sở dữ liệu đầu vào để tính toán khẩu độ cống, vị trí đặt đập ngăn mặn và tuyến đê bao.
  • Cộng đồng học thuật và nghiên cứu sinh: Giảng viên, nhà nghiên cứu và học viên cao học thuộc các ngành Biến đổi khí hậu, Thủy văn học, Khoa học Môi trường và Quản lý Đất đai có thể tham khảo quy trình tích hợp mô hình MIKE 11 và phương pháp chuẩn độ hóa nghiệm Knudsen làm tài liệu giảng dạy và nghiên cứu mở rộng.
  • Doanh nghiệp nông nghiệp và hợp tác xã ven biển: Các tổ chức kinh tế nông nghiệp có thể ứng dụng các cảnh báo về thời điểm đỉnh mặn và bản đồ ranh mặn để chủ động lập kế hoạch tích trữ nước ngọt nội đồng, bảo vệ mùa màng và chuyển đổi mô hình kinh doanh thủy sản thích ứng.

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

  • Ngưỡng độ mặn nào là an toàn cho nước tưới nông nghiệp tại đồng bằng sông Hồng? Theo tiêu chuẩn phân loại của FAO và điều kiện canh tác thực tế tại Nam Định, nguồn nước có độ mặn dưới 1,5‰ được coi là an toàn cho việc tưới dưỡng lúa và hoa màu. Khi độ mặn vượt quá 1,5‰, nước bắt đầu gây ức chế sinh trưởng cây trồng; nếu độ mặn đạt ngưỡng 4‰ trở lên, cây lúa sẽ bị cháy lá, héo khô và chết hàng loạt do sốc thẩm thấu.

  • Tại sao sông Ninh Cơ lại ghi nhận hiện tượng xâm nhập mặn sâu và khốc liệt hơn sông Đáy? Hiện tượng này bắt nguồn từ sự kết hợp giữa địa mạo lòng dẫn và động lực dòng chảy. Sông Ninh Cơ có lòng dẫn thẳng, dốc thoải và khẩu độ mở rộng về phía cửa Lạch Giang, tạo điều kiện cho sóng triều truyền vào với tốc độ cao từ 15 đến 20 km/h mà ít bị tiêu hao năng lượng. Trong khi đó, dòng chảy sông Đáy chịu sự phân lưu phức tạp và có các khúc uốn giúp giảm bớt cường độ xâm nhập của triều mặn.

  • Mô hình MIKE 11 có những ưu thế vượt trội nào trong tính toán lan truyền mặn? MIKE 11 tích hợp đồng bộ module thủy lực một chiều HD và module truyền tải khuếch tán AD, cho phép giải hệ phương trình Saint-Venant và định luật Fick với tốc độ tính toán nhanh và độ ổn định cao. Mô hình dễ dàng liên kết với các hệ thống GIS, xử lý linh hoạt mạng lưới sông phức hợp và phản ánh chính xác tác động đồng thời của lưu lượng thượng nguồn và triều cường hạ lưu.

  • Dưới kịch bản biến đổi khí hậu RCP 4.5, xâm nhập mặn tại Nam Định sẽ biến đổi ra sao đến năm 2050 và 2100? Theo kịch bản RCP 4.5, mực nước biển dâng thêm 22 cm vào năm 2050 và 53 cm vào năm 2100. Dưới tác động này, ranh mặn 1,5‰ dự tính sẽ tịnh tiến sâu thêm vào nội đồng từ 3 đến 7 km so với hiện trạng năm 2010. Đặc biệt, nếu mực nước biển dâng 50 cm, khoảng 26% diện tích toàn tỉnh Nam Định sẽ đối diện với nguy cơ ngập úng và nhiễm mặn.

  • Phương pháp Knudsen xác định độ mặn trong phòng thí nghiệm dựa trên nguyên lý gì? Phương pháp Knudsen dựa trên định luật hằng số tỉ lệ giữa các thành phần ion chính trong nước biển, trong đó ion Cl- chiếm tỷ trọng lớn nhất (khoảng 55,61%). Bằng cách chuẩn độ dung dịch mẫu với Bạc Nitrat (AgNO3) và chỉ thị màu Kali Cromat (K2CrO4), lượng ion Clo được xác định chính xác, từ đó tính toán ra độ Clo và quy đổi trực tiếp sang tổng độ mặn với độ tin cậy tuyệt đối trong dải pH trung tính.

Kết luận

  • Luận văn đã xác lập bức tranh tổng thể, chi tiết về hiện trạng xâm nhập mặn tại phía Nam đồng bằng sông Hồng, chứng minh quy luật mặn đạt đỉnh vào các tháng 1, 2, 3 mùa kiệt và phân hóa rõ rệt giữa hai dòng sông Đáy và Ninh Cơ.
  • Ứng dụng thành công bộ mô hình số MIKE 11 kết hợp dữ liệu khảo sát thực địa 48 mẫu nước, xây dựng công cụ mô phỏng đáng tin cậy phản ánh mối tương quan giữa thủy động lực dòng chảy và lan truyền nồng độ muối.
  • Dự tính định lượng sự gia tăng phạm vi nhiễm mặn theo các kịch bản phát thải RCP 4.5 và RCP 8.5 đến năm 2100, cung cấp căn cứ khoa học cho thấy ranh mặn 1,5‰ có thể xâm nhập sâu thêm tới 7 km vào nội đồng.
  • Đề xuất hệ thống giải pháp toàn diện kết hợp công trình ngăn mặn giữ ngọt với các biện pháp phi công trình như quan trắc SCADA thời gian thực, điều tiết xả nước hồ chứa thủy điện và chuyển đổi mô hình kinh tế nông nghiệp thích ứng.
  • Đóng góp bộ cơ sở dữ liệu số hóa quan trọng cho tỉnh Nam Định và vùng đồng bằng Bắc Bộ; lộ trình nghiên cứu tiếp theo cần mở rộng tích hợp mô hình thủy văn hai chiều (MIKE 21) và đánh giá chi tiết quá trình mặn hóa tầng chứa nước ngầm nhằm bảo vệ toàn diện an ninh nguồn nước ven biển.