Mục tiêu của nghiên cứu này là đánh giá tác động của NBD đến sự phân phối lại nguồn bùn cát tại các khu vực khác nhau tại vùng cửa sông Soài Rạp. Phương pháp mô hình toán Telemac-Sisyphe-Tomawac, trong đó vận chuyển bùn cát hỗn hợp (mix-sediment) được sử dụng trong nghiên cứu này. Mời các bạn cùng tham khảo!
BÀI BÁO KHOA HỌC TÁC ĐỘNG CỦA NƯỚC BIỂN DÂNG ĐẾN VẬN CHUYỂN VÀ PHÂN PHỐI NGUỒN BÙN CÁT TẠI CỬA SƠNG SỒI RẠP Lê Ngọc Anh1,2, Hồng Trung Thống2 Tóm tắt: Q trình biến đổi hình thái cửa sơng thông tin cần thiết cho việc hoạch định chiến lược quản lý tài nguyên nước Trong điều kiện nước biển dâng (NBD), chế độ thủy động lực vận chuyển bùn cát vùng cửa sông thay đổi theo, hệ làm thay đổi lại nguồn bùn cát từ thượng lưu đến cửa biển Mục tiêu nghiên cứu đánh giá tác động NBD đến phân phối lại nguồn bùn cát khu vực khác vùng cửa sơng Sồi Rạp Phương pháp mơ hình tốn Telemac-Sisyphe-Tomawac, vận chuyển bùn cát hỗn hợp (mix-sediment) sử dụng nghiên cứu Các kịch NBD+0.0 m (kịch nền), NBD+0.5 m, NBD+1.0 m đưa để dự báo trình vận chuyển bùn cát Kết mô sau năm cho thấy rằng, so với kịch nền, KV1, tổng lượng bùn cát giảm -579.103 m3 (12,8%) NBD+0.5 m, -913.103 m3 (-35%) NBD+1.0 m; KV2, thay đổi không đáng kể; KV3, tổng lượng bùn cát tăng lên 3,5.106 m3 (10%) NBD+0.5 m, 6,9.103 m6 (20%) NBD+1.0 m Từ khố: Telemac, Tomawac, Sisyphe, Cửa sơng Sồi Rạp, vận chuyển bùn cát, hình thái lịng dẫn GIỚI THIỆU * Nước biển dâng thách thức lớn khu vực ven biển vùng cửa sông tương lai, mà nhiều chứng gia tăng mực nước biển đại dương giới ngày thể rõ Mực nước biển tồn cầu tăng trung bình 1,7 mm/năm kỷ 20 (IPCC, 2007) Tại Nha Trang (biển Đông – Việt Nam), từ năm 1976 đến năm 2008, cho thấy xu biến đổi mực nước biển với chu kì 5,7 năm; từ 1976 đến 1992 mực nước (theo xu thế) giảm từ 1993 đến 2008 mực nước tăng (Nguyễn Tác An, et al., 2015) Tại khu vực cửa biển Sồi Rạp, biến trình mực nước trung bình 18 năm Vũng Tàu cho thấy gia tăng mực nước biển, mực nước lớn trung bình 18 năm thời kỳ 19902007 cao mực nước lớn trung bình 18 năm thời kỳ 1982-1999 46,7 mm, trung bình năm gia tăng 6,2 mm (Hội Đập lớn & Phát triển Nguồn nước Việt Nam, 2012) Nhiều nghiên NCS Trường ĐH Bách Khoa Tp.HCM Trường Đại học Tài Nguyên & Môi trường TP.HCM 120 cứu tác động NBD lên đặc tính triều, chế độ thủy lực, xâm nhập mặn biến đổi hình thái sơng/biển, thay đổi đường bờ Tại cửa sơng Sồi Rạp, NBD làm thay đổi đặc tính triều, mực nước đáy đỉnh triều tăng, gia tăng chân triều lớn so với đỉnh triều (Lưu Xuân Lộc, et al., 2015); gây biến đổi vận tốc dòng chảy mùa lũ kiệt vùng cửa sông Soài Rạp (Lê Ngọc Anh, et al., 2017) NBD làm cho hình thái đường bờ biển, cửa sơng thay đổi theo chiều hướng bất lợi: xói lở mạnh hơn, hình dạng đường bờ thay đổi theo hướng xâm thực (Hồng Văn Hn, et al.) Cửa sơng Sồi Rạp nằm vùng hạ lưu hệ thống sông Đồng Nai – Sài Gịn thuộc hệ thống sơng Đồng Nai Cửa Sồi Rạp điểm cuối hệ thống sông Đồng Nai có diện tích lưu vực 40.700 km2 NBD gây thay đổi chế độ thủy động lực, vận chuyển bùn cát hình thái lịng dẫn vùng cửa sơng Sồi Rạp, ảnh hưởng đến giao thơng thủy gây tình trạng sạt lở bờ phức tạp năm gần Nghiên cứu trình biến đổi hình thái cửa sơng KHOA HỌC KỸ THUẬT THỦY LỢI VÀ MÔI TRƯỜNG - SỐ 74 (6/2021) điều kiện NBD thông tin cần thiết cho việc hoạch định chiến lược quản lý tài nguyên nước tương lai PHƯƠNG PHÁP NGHIÊN CỨU Trong nghiên cứu chúng tơi sử dụng mơ hình tốn Telemac-2D – Tomawac – Sisyphe kết nối với để mơ q trình thủy động lực đồng thời với q trình biến đổi hình thái (Hình 2) Trong mơ hình Telemac-2D mơ q trình thủy động lực đóng vai trị trung tâm kết nối với mơ hình khác Tomawac mơ sóng có xét đến tác động sóng lên dịng ngược lại, đồng thời cung cấp thơng số cho việc tính tốn vận chuyển bùn cát Sisyphe mơ q trình vận chuyển bùn cát, cung cấp thông tin thay đổi đáy cập nhật cho mơ hình khác bước thời gian mô Cơ sở lý thuyết hệ thống mơ hình Telemac-2D – Sisyphe - Tomawac trình bày chi tiết (ATA Riadh, et al., 2014; TASSI Pablo, et al., 2014; Fouquet Thierry, 2015) Trong nghiên cứu này, thay đổi lại cách ứng xử vận chuyển bùn cát dựa đặc tính bùn cát hỗn hợp thay đổi theo khơng gian thông qua việc thiết lập biểu thức quan hệ (ES1, ES2) = f(f1, f2, ES, n) trình bày chi tiết báo (Lê Ngọc Anh, et al., 2020) Cấu trúc file liệu để đưa vào để xác định thành thần ES1 (chiều dày lớp cát) ES2 (chiều dày lớp bùn) layer dựa vào chiều dày layer1, layer2 (ES_1, ES_2), tỷ lệ khối lượng bùn (f2-1, f2-2) minh họa Hình 1-a THIẾT LẬP MƠ HÌNH Các bước thiết lập mơ hình, hiệu chỉnh kiểm định mơ hình, thơng số cài đặt trình bày chi tiết báo (Lê Ngọc Anh, et al., 2020) (a) Hình (a) Cấu trúc file liệu áp đặt bùn cát hỗn hợp theo không gian (Point: số thứ tự điểm miền tính; f2-1, f2-1 tỷ lệ khối lượng bùn layer 1, 2; C1, C2 hàm lượng bùn đất ứng với layer 1, (kg/m3); (b) phạm vi khu vực KHOA HỌC KỸ THUẬT THỦY LỢI VÀ MÔI TRƯỜNG - SỐ 74 (6/2021) 121 hCk hUCk hVCk C C h s k h s k E k D k t x y x x y y E k Dk Zref s C eqk C refk (1) (2) Hình Sơ đồ chung kết nối mơ hình TELEMAC; (U,V vận tốc theo phương; H: chiều sâu nước; Fx, Fy thành phần lực ứng suất sóng gây theo phương), k: tên class (k =1 cát, k =2 bùn) KẾT QUẢ NGHIÊN CỨU Trong nghiên cứu này, biến đổi hình thái cửa sơng Sồi Rạp mô dựa kịch NBD (+0,5 m, +1,0 m) kịch NBD +0,0 m với thời kỳ mô kéo dài năm NBD gây thay đổi chế độ triều vùng cửa sơng dẫn đến thay đổi q trình tích lũy bùn cát Tại khu vực KV1, dựa đồ thị Hình 5, nhận định rằng, NBD gây thay đổi lượng cát tích lũy đây, mức độ thay đổi phụ thuộc vào thời kỳ triều theo mùa Vào mùa kiệt, thay đổi lượng bùn cát khu vực không đáng kể, đường lũy tích ứng với kịch NBD tiệm cận Tuy nhiên, vào mùa lũ, ảnh hưởng NBD lớn so với mùa kiệt; NBD có xu hướng làm giảm lượng bùn cát bồi lắng gia tăng xói kỳ triều cường, tăng lượng bùn cát bồi 122 kỳ triều Kết tổng hợp lượng bùn cát tích lũy ngày triều điển hình trình bày bảng cho thấy, vào mùa kiệt, kỳ triều cường, NBD làm giảm lượng bùn cát tích lũy đây, tổng lượng bùn cát chu kỳ triều 2.900 m3 (NBD+0,0 m), 1.700 m3 (NBD+0,5 m), 1.200 m3 (NBD+1,0 m); kỳ triều thay đổi không đáng kể Vào mùa lũ, NBD làm gia tăng lượng bùn cát dịch chuyển phía cửa sơng làm gia tăng xói mạnh kỳ triều cường; tổng lượng bùn cát chu kỳ triều –97.200 m3 (NBD +0,0 m), –102.300 m3 (NBD+0,5 m), – 109.200 m (NBD+1,0 m) Đối với khu vực KV2, dựa đồ thị Hình 5-a thấy rằng, NBD làm giảm lượng bùn cát bồi lắng mùa kiệt tăng lượng bùn cát bồi lắng mùa lũ so với kịch Vào mùa kiệt, KHOA HỌC KỸ THUẬT THỦY LỢI VÀ MÔI TRƯỜNG - SỐ 74 (6/2021) NBD làm giảm lượng bùn cát bồi lắng đây, kỳ triều cường, tổng lượng bùn cát chu kỳ triều 36.150 m3 (NBD+0,0 m), 120.950 m3 (NBD+0,5 m) 7.200 m3 (NBD+1,0 m); kỳ kém, tổng lượng bùn cát chu kỳ triều 510 m3 (NBD + 0,0 m), 180 m3 (NBD + 0,5 m) –190 m3 (NBD+1,0 m) Vào mùa lũ, NBD làm tăng lượng bùn cát bồi lắng đây, kỳ triều kém, lượng bùn cát bồi có thay đổi không đáng kể so với kịch nền; kỳ triều kém, NBD làm gia tăng lượng bùn cát bồi đây, tổng lượng bùn cát chu kỳ triều 157.400 m3 (NBD+0,0 m), 170.700 m3 (NBD + 0,5 m) 190.200 m3 (NBD + m) Đối với KV3, đồ thị Hình 5-b cho thấy rằng, NBD làm gia tăng lượng bùn cát lắng đọng so với kịch Vào mùa kiệt, kỳ triều cường, xét chu kỳ triều, tổng lượng bùn cát bồi 114.000 m3 (NBD + 0,0 m), 148.000 m3 (NBD + 0,5 m) 183.000 m3 (NBD + 1,0 m); kỳ triều kém, tổng lượng bùn cát bồi –48.000 m3 (NBD + 0,0 m), –37.900 m3 (NBD + 0,5 m), –26.500 m3 (NBD + 1,0 m) Vào mùa lũ, kỳ triều cường, xét chu kỳ triều, tổng lượng bùn cát bồi 357.000 m3 (NBD + 0,0 m), 388.700 m3 (NBD + 0,5 m) 440.800 m3 (NBD + 1,0 m); kỳ triều kém, tổng lượng bùn cát bồi 448.000 m3 (NBD +0,0 m), 469.800 m3 (NBD + 0,5 m) 485.500 m3 (NBD + 1,0 m) Tổng lượng cát tích lũy năm khu vực khác thể đồ thị Hình 5-d nhận định rằng, khu vực KV1, NBD làm giảm lượng bùn cát bồi lắng đây, mùa lũ ảnh hưởng lớn so với mùa kiệt Tổng lượng bùn cát lắng đọng năm 2,6 tr.m3 (NBD+0,0 m), 2,0 tr.m3 (NBD+0,5 m) 1,7 tr.m3 (NBD+1,0 m) Đối với khu vực KV2, NBD làm giảm lượng bùn cát lắng đọng mùa kiệt lại làm gia tăng lượng bùn cát lắng đọng mùa lũ, gia tăng bù đắp thiếu hụt mùa kiệt nên năm tổng lượng bùn cát không thay đổi nhiều điều kiện NBD Tổng lượng bùn cát lắng đọng năm 10,8 tr.m3 (NBD+0,0 m), 10,8 tr.m3 (NBD+0,5 m) 10,6 tr.m3 (NBD+1,0 m) Đối với khu vực KV3, NBD làm gia tăng lượng bùn cát tích lũy năm, vào mùa lũ gia tăng lớn so với mùa kiệt Tổng lượng bùn cát lắng đọng năm 32,0 tr.m3 (NBD+0,0 m), 34,8 tr.m3 (NBD+0,5 m) 37,4 tr.m3 (NBD+1,0 m) Các khu vực ngồi biển Đơng, NBD có xu hướng làm giảm lượng bùn cát; riêng khu vực BĐ2 (vùng biển bên ngồi cửa sơng Tiền), NBD làm gia tăng thêm lượng bùn cấp Tổng lượng bùn cát nhận khu vực biển Đông (BD1+BD2+BD3) gần cửa sông 335,2 tr.m3 (NBD+0,0 m), 325,4 tr.m3 (NBD+0,5 m) 313,4 tr.m3 (NBD+1,0 m) Sự phân phối lại tổng lượng bùn cát điều kiện nước biển dâng làm thay đổi đáy cửa sơng Sồi Rạp Diễn biến tích lũy bùn cát khu vực diễn biến theo thời gian trình bày bảng Theo đó, thấy KV1, chịu tác động mạnh NBD, suy giảm mạnh lượng bùn cát điều kiện NBD kết gia tăng lượng vận chuyển bùn cát phía cửa sơng triều rút xảy mạnh vào mùa lũ Hình thể thay đổi đáy tương ứng với kịch NBD+0,5 m NBD+1,0 m so với kịch +0,0 m sau năm mô Kết cho thấy, NBD gây bồi khu vực cửa sông khu vực nước nơng gây xói mạnh cho khu vực thượng lưu KẾT LUẬN VÀ THẢO LUẬN Sóng triều xem sóng dài với vận tốc truyền sóng tính , điều kiện NBD, chiều sâu nước h tăng lên làm cho vận tốc truyền sóng tăng, dẫn đến pha số sóng tăng lên Trong nghiên cứu này, xử dụng phương pháp phân tích điều hịa vị trí cửa sơng Sồi Rạp thành phần sóng (4 sóng bán nhật: M2, S2, N2, K2 ; sóng tồn nhật: O1, P1, Q1, K1; sóng nước nơng: M4, M6, MS4) Kết cho thấy, NBD làm tăng biên độ hầu hết phân triều Đối với thành phần sóng bán nhật, sóng M2 (sóng mặt trăng chính) có biên độ tăng lớn +2,49 cm (NBD + 0,5 m), +4,43 cm (NBD + m) sớm pha 7,5 phút (NBD + 0,5 m), 14,1 phút (NBD + 1,0 m) so với kịch (NBD + 0,0 m) KHOA HỌC KỸ THUẬT THỦY LỢI VÀ MÔI TRƯỜNG - SỐ 74 (6/2021) 123 NBD làm tăng biên độ triều, gia tăng động lực triều vận chuyển bùn cát biển Đông bên hệ thống sông Kết KV3, tổng lượng bùn cát tăng lên trái ngược với khu vực phía thượng lưu (KV1) tổng lượng bùn cát có xu hướng giảm KV3 nguồn bùn cát tăng lên (xem Hình 3) Nguồn tài ngun cát sơng có giá trị cao xây dựng (cát xây) thường phân bố KV1 có nguy giảm đi; nguồn cát KV3 tăng ên nguồn cát chủ yếu cát san lấp có giá trị thấp Do đó, đặt vấn đề cần phải xây dựng chiến lược khai thác tài nguyên hợp lý tương lai Trong điều kiện NBD, nguồn bùn cát có xu hướng dịch chuyển mạnh lên phía Bắc làm cho tổng lượng bùn cát khu vực BD3 giảm khu vực BD2 tăng, đặc biệt xảy mạnh thời kỳ gió mùa Tây Nam với thời điểm mùa lũ mà lượng phù sa từ lưu vực sông đổ biển nhiều % % Hình Mức độ ảnh hưởng NBD đến nguồn bùn cát khu vực (tính % KB0) Hình Biến đổi đáy so sánh kịch NBD+0,5 m, NBD+1,0 m, so với KB +0,0 m Bảng Tổng lượng bùn cát thống kê chu kỳ triều đại diện Đvt: 1.000 m3 Mùa Mùa kiệt 124 Kỳ triều Triều cường Triều Kịch +0,0 +0,5 +1,0 +0,0 +0,5 +1,0 Bồi 68,9 71,5 75,8 48,6 46,9 45,6 KV1 Xói -66,0 -69,8 -74,6 -32,5 -31,4 -29,5 Tổng 2,9 1,7 1,2 16,1 15,5 16,1 Bồi 74,4 72,3 69,3 3,9 3,3 2,7 KV2 Xói -38,3 -51,4 -62,1 -3,4 -3,1 -2,9 Tổng 36,2 20,9 7,2 0,5 0,2 -0,2 Bồi 535,4 532,4 518,7 31,5 27,0 21,7 KV3 Xói -421,3 -384,3 -335,7 -79,5 -64,9 -48,2 Tổng 114,1 148,1 183,0 -48,0 -37,9 -26,5 KHOA HỌC KỸ THUẬT THỦY LỢI VÀ MÔI TRƯỜNG - SỐ 74 (6/2021) Mùa Mùa lũ Kỳ triều Triều cường Triều Kịch +0,0 +0,5 +1,0 +0,0 +0,5 +1,0 Bồi 78,5 79,1 79,3 139,4 148,3 156,1 KV1 Xói -175,7 -181,4 -188,9 -85,1 -86,7 -87,9 Tổng -97,2 -102,3 -109,6 54,3 61,6 68,2 Bồi 111,5 127,4 138,9 231,8 241,9 260,6 KV2 Xói -65,6 -78,8 -89,6 -74,3 -71,2 -70,3 Tổng 45,9 48,6 49,3 157,5 170,7 190,3 Bồi 604,5 637,6 680,3 546,9 551,8 544,9 KV3 Xói -246,7 -248,8 -239,4 -98,8 -81,9 -59,3 Tổng 357,7 388,8 440,9 448,1 469,9 485,6 Bảng Tổng lượng bùn cát tháng khu vực ứng với trường hợp NBD Đvt: 1.000 m3 Tháng 10 11 12 Năm +0,0 m 967 911 844 704 556 161 138 -367 -646 -487 -133 -46 2.602 KV1 +0,5 m 940 898 834 702 550 123 14 -487 -701 -497 -242 -113 2.023 +1,0 m 923 886 835 721 567 134 -41 -536 -749 -545 -341 -165 1.689 +0,0 m 908 451 459 786 1.593 1.558 1.188 40 807 1.872 834 317 10.813 KV2 +0,5 m 780 321 324 612 1.330 2.053 1.425 312 876 1.898 740 210 10.880 (a) KV1 KHOA HỌC KỸ THUẬT THỦY LỢI VÀ MÔI TRƯỜNG - SỐ 74 (6/2021) +1,0 m 633 224 168 462 1.014 2.382 1.800 579 989 1.833 548 -21 10.611 +0,0 m -1.442 -58 465 369 1.009 2.573 6.574 5.346 5.176 8.215 4.592 2.092 34.911 KV3 +0,5 m -725 142 616 370 1.063 2.185 6.798 5.676 5.464 8.874 5.165 2.828 38.456 +1,0 m 57 404 811 385 1.062 1.769 6.764 5.865 5.808 9.763 5.742 3.449 41.879 (b) KV2 125 (d) (c) KV3 Hình Tổng lượng cát tích lũy ngày triều điển hình (a) KV1, (b) KV2, (c) KV3 (+): bồi, (-): xói Ký hiệu số 1: diễn biến lượng bùn cát thời điểm, 2: lượng bùn cát tích lũy theo thời gian (d) Tổng lượng cát tích lũy TÀI LIỆU THAM KHẢO Hoàng Văn Huân, Trần Thị Xuân Mỹ 'Tác động trình nước biển dâng vùng cửa sông, ven biển đồng Nam Bộ định hướng hành động ứng phó', Tuyển tập KHCN 50 năm XD&PT Hội Đập lớn & Phát triển Nguồn nước Việt Nam (2012) Quy hoạch thủy lợi Đồng Sơng Cửu Long điều kiện biến đổi khí hậu nước biển dâng http://www.vncold.vn Available at: http://www.vncold.vn/Web/Content.aspx?distid=2931 (Accessed: 15/02/2012 2012) Lê Ngọc Anh, Hồng Trung Thống, Châu Ngơ Bảo (2020) 'Thiết lập mơ hình vận chuyển bùn cát hỗn hợp theo khơng gian mơ diễn biến hình thái lịng dẫn cửa sơng Sồi Rạp', Khoa học kỹ thuật Thủy lợi Môi trường, 69, pp 175-181 Lê Ngọc Anh, Lộc Lưu Xuân (2017) 'Tác động nước biển dâng đến chế độ dòng chảy cửa sơng Sồi Rạp', Tài ngun Nước, Số chun đề, pp 30-38 Lưu Xuân Lộc, Nguyễn Thống, et al (2015) 'Tác động nước biển dâng đến chế độ thủy triều sơng Sài Gịn - Đồng Nai', Nơng nghiệp Phát triển Nông thôn, pp 60-65 Nguyễn Tác An, Nguyễn Kỳ Phùng 'Biến đổi khí hậu giải ven bờ biển tỉnh Khánh Hịa, tiếp cận thích ứng ứng phó' Tập 1: Khí tượng - Khí hậu, Khí tượng Nơng nghiệp Biến đổi Khí hậu Hội thảo khoa học Quốc gia Khí tượng Thuỷ văn, Mơi trường Biến đổi Khí hậu, 231-237 ATA Riadh, Goeury Cédric, Hervouet Jean Michel (2014), Telemac2d v7.0 User's Manual, EDF-R&D Fouquet Thierry (2015), Tomawac release 7.1 Operating manual, EDF R&D IPCC (2007), Climate Change 2007: The Physical Science Basis Contribution of Working Group I to the Fourth Assessment Report of the Intergovernmental Panel on Climate Change Cambridge University Press TASSI Pablo, GOUTAL Nicole, BODEL Charles (2014), Sisyphe v6.3 User's Manual, EDF-R&D 126 KHOA HỌC KỸ THUẬT THỦY LỢI VÀ MÔI TRƯỜNG - SỐ 74 (6/2021) Abstract: IMPACTS OF SEA LEVEL RISE ON TRANSPORTATION OF SEDIMENT REALLOCATION AMOUNT OF SEDIMENT AT SOAI RAP ESTUARY The morphology changing process of an estuary is essential information for planning strategies of water resources management In the context of sea level rise (SLR), the hydrodynamic mode and transportation of sediment at the estuaries also changed, as a result it leads to a change in the source of sediment from upstream to the sea The objective of this study is to evaluate the impact of SLR on redistribution of sediment resources in different areas at the Soai Rap estuary Telemac-SisypheTomawac mathematical modeling method, in which transportation of sediment with mixed sediment characteristics (mix-sediment) was used in this study The scenarios of SLR + 0.0 m (base scenario), SLR + 0.5 m, SLR + 1.0 m were proposed to forecast the process of sediment transportation The results of simulation after one year show that, compared to the base scenario, for zone KV1, the total amount of sediment decreases by –579.103 m3 (12.8%) SLR +0.5 m, –913.103 m3 (-35%) SLR +1.0 m; for zone KV2, the change is not significant; for KV3, the total amount of sediment increases by 3.5*106 m3 (10%) SLR + 0.5 m, 6.9*103 m6 (20%) SLR + 1.0 m Keywords: Telemac, Tomawac, Sisyphe, Soai Rap estuary, sediment transport, morphology Ngày nhận bài: 21/5/2021 Ngày chấp nhận đăng: 18/6/2021 KHOA HỌC KỸ THUẬT THỦY LỢI VÀ MÔI TRƯỜNG - SỐ 74 (6/2021) 127 ... (2017) 'Tác động nước biển dâng đến chế độ dịng chảy cửa sơng Soài Rạp' , Tài nguyên Nước, Số chuyên đề, pp 30-38 Lưu Xuân Lộc, Nguyễn Thống, et al (2015) 'Tác động nước biển dâng đến chế độ thủy... lượng bùn cát có xu hướng giảm KV3 nguồn bùn cát tăng lên (xem Hình 3) Nguồn tài ngun cát sơng có giá trị cao xây dựng (cát xây) thường phân bố KV1 có nguy giảm đi; nguồn cát KV3 tăng ên nguồn cát. .. Trần Thị Xuân Mỹ 'Tác động trình nước biển dâng vùng cửa sông, ven biển đồng Nam Bộ định hướng hành động ứng phó', Tuyển tập KHCN 50 năm XD&PT Hội Đập lớn & Phát triển Nguồn nước Việt Nam (2012)