to study on ng dng mo hinh toan nghien

Tạp chí Khoa học Cơng nghệ Biển T12 (2012) Số Tr - 11 ỨNG DỤNG MÔ HÌNH TỐN NGHIÊN CỨU VÙNG NƯỚC ĐỤC NHẤT Ở VÙNG CỬA SÔNG BẠCH ĐẰNG VŨ DUY VĨNH, TRẦN ĐỨC THẠNH Viện Tài ngun Mơi trường Biển Tóm tắt: Các nghiên cứu động lực tương tác sông biển - biển vùng cửa sông gần biết đến xuất vùng nước đục (Maximum Turbidity Zone MTZ) so với xung quanh Đây nơi giàu dinh dưỡng liên quan đến tài nguyên sinh vật phong phú, nơi có khả tích luỹ cao chất gây ô nhiễm Các kết nghiên cứu mơ hình Delf 3D cho phép phát vùng MTZ vùng cửa sông Bạch Đằng Chúng có quy mơ vị trí ln biến đổi theo vai trò ảnh hưởng dòng vật chất từ lục địa biển theo mùa dao động mực nước thuỷ triều MTZ có vị trí xa cách bờ khoảng 15km vào pha triều xuống mùa mưa Vào mùa khơ, chúng có qui mơ nhỏ vị trí xuất gần bờ I MỞ ĐẦU Hình Sơ đồ mơ tả khái niệm vùng nước đục cửa sông (MTZ) [2] Các vùng cửa sơng (VCS) chiếm chưa tới 10% diện tích bề mặt đại dương giới [14] lại có vai trị quan trọng chu trình vật chất toàn cầu Gần đây, nghiên cứu động lực tương tác lục địa - biển VCS phát triển mạnh mẽ, bật lên vấn đề xuất vùng nước đục (Maximum Turbidity Zone MTZ), (hình 1) Đây nơi hội tụ - tương tác khối nước sông nằm có hàm lượng chất lơ lửng cao nằm khối nước biển mặn có tỷ trọng lớn nằm phía hình thành nên nêm mặn [1, 5, 7, 14] MTZ khơng nơi lắng đọng, tích tụ vật chất từ lục địa đưa ra, giàu chất dinh dưỡng [15] mà nơi tích tụ hố chất độc hại kim loại nặng [3, 4, 18] Các chất gây ô nhiễm hấp thụ vào sinh vật phù du, qua chuỗi thức ăn ảnh hưởng đến toàn hệ sinh thái khu vực Mặt khác, MTZ nơi có nguồn thức ăn phong phú nên loài cá sinh vật khác tập trung đến sinh sống đẻ trứng [8, 10] Vì vậy, nghiên cứu tương tác sơng - biển thơng qua chếhình thành, xuất hiện, vị trí yếu tố mơi trường MTZ VCS có ý nghĩa khoa học thực tiễn VCS Bạch Đằng nằm Hải Phòng Quảng Ninh hàng năm nhận 10km3 nước từ lục địa dao động thuỷ triều lớn, cực đại 4,35m, chế độ sóng dịng chảy sơng biến đổi mạnh theo mùa Bài báo trình bày kết bước đầu việc ứng dụng mơ hình tốn học nghiên cứu vùng nước đục - hệ trình tương tác lục địa - biển VCS II TÀI LIỆU VÀ PHƯƠNG PHÁP Tài liệu Các tài liệu liệu gió, sóng, mực nước độ muối trạm hải văn Hòn Dáu từ năm 1970 đến 2007 Số liệu quan trắc, lưu lượng nước hàm lượng trầm tích lơ lửng sông Bạch Đằng, Cấm, Lạch Tray, Văn Úc Thái Bình thời gian gần Tài liệu địa hình từ đồ địa hình tỷ lệ 1:50000, hệ lưới chiếu VN2000 khu vực ven biển Hải Phòng Cơ sở liệu số điều hòa thủy triều FES2004 [16] sở liệu nhiệt muối World Ocean Atlas 2005 [9] Phương pháp Phương pháp thiết lập mơ hình tốn học tổng hợp (mơ hình Delf-3D) song - dịng chảy - thủy triều vận chuyển trầm tích lơ lửng cho khu vực ven bờ cửa Bạch Đằng Khu vực tính tốn độ sâu số hố từ đồ địa hình khu vực ven biển Hải Phòng từ mảnh ghép F48-83, F48-82, F48-94 F48-95 Đây đồ địa hình 1:50000 hệ tọa độ UTM - VN2000, Trung tâm Trắc địa Bản đồ, Bộ Tài nguyên Mơi trường xuất năm 2004 Ngồi ra, cịn có số tài liệu địa hình đo đạc bổ sung gần Kích thước phạm vi miền tính thể hình 2, với chiều đơng Bắc - Tây Nam khoảng 34km Tây Bắc - Đông Nam khoảng 65km Hệ thống lưới cong trực giao chọn sử dụng để tính tốn với lưới có kích thước từ 40,4 đến 415,1m tồn miền tính chia làm 257 × 460 điểm Chiều thẳng đứng từ mặt xuống đáy chia làm lớp độ sâu, lớp 20%H (H độ sâu điểm tính) Mơ hình thuỷ động lực áp dụng tính tốn cho trường hợp: mùa mưa từ 0h ngày 01/8/2007 đến 0h ngày 30/8/2007 mùa khô từ 0h ngày 1/3/2007 đến 0h ngày 30/3/2008 Bước thời gian tính tốn mơ hình thuỷ động lực 0,5 phút Các trình vật lý tính đến mơ hình bao gồm lực Coriolis, ma sát đáy, nhiệt độ nước độ muối, vận tốc hướng gió °4 ' °5 ' °0 ' Hải Phòng 20°50' 20°50' Cát Hải Cát Bà Đồ Sơn 20°40' 20°40' Độ sâu (m) 20°30' 20°30' °4 ' °5 ' °0 ' Hình Lưới độ sâu miền tính Điều kiện ban đầu mơ hình thủy động lực chọn giá trị cho mực nước, nhiệt độ độ muối Ở lần chạy lại sau, điều kiện ban đầu sử dụng restart file từ kết lần chạy trước Điều kiện biên mở Điều kiện nhiệt độ, độ muối dao động mực nước biên mở phía biển mơ hình thủy động lực sử dụng kết NESTHD mơ hình khác có phạm vi miền tính lớn bao bên ngồi [21] Tại biên mở sơng chính, điều kiện biên lưu lượng nước, hàm lượng trầm tích lơ lửng sử dụng từ kết khảo sát thủy văn nhóm CAMELIA mùa mưa 2008 mùa khơ năm 2009 cửa sơng Hải Phòng [12] Nhằm đánh giá biến đổi dòng chảy, độ muối từ lục địa phía biển, mặt cắt hướng Tây Bắc - Đông Nam gần trùng hướng chảy cửa sơng Bạch Đằng thiết lập (hình 3) Trong số này, Mặt cắt I gần trùng trục luồng Nam Triệu Mặt cắt IV khống chế khu vực Đông Bắc Tây Nam bán đảo Đồ Sơn Để kiểm chứng hiệu chỉnh kết tính tốn mơ hình, chúng tơi tiến hành so sánh kết tính tốn quan trắc vận tốc dòng chảy, mực nước, hàm 106°40' 106°50' 107°00' Cát Hải 107°10' 20°50' 20°50' lượng trầm tích lơ lửng khu vực nghiên cứu nhận thấy phù hợp tương đối số liệu mực nước tính tốn mực nước quan trắc Hòn Dấu tháng 8/2007 tháng 3/2008 Sau hiệu chỉnh, hệ số tương quan tương ứng 0,9716 0,9658 mùa mưa mùa khô Sai số chấp nhận điều kiện địa hình phức tạp biên độ dao động mực nước khu vực tính lớn Cát Bà Hải Phòng Mặt cắt I 20°40' 20°40' Đồ Sơn Mặt cắt II Mặt cắt III 20°30' 20°30' Mặt cắt IV Mặt cắt V 106°40' 106°50' 107°00' 107°10' Hình Vị trí mặt cắt tính tốn III KẾT QUẢ VÀ THẢO LUẬN Phân bố trầm tích lơ lửng vùng cửa sông Bạch Đằng Các kết nghiên cứu MTZ có sử dụng cơng cụ tốn học mơ hình [6, 11, 13, 17, 19, 20, 22] hoàn lưu mật độ (gravitational circulation) chênh lệch tỷ trọng nước nước mặn, lắng đọng trầm tích lở lửng tái lơ lửng bào mòn đáy tác động dòng triều ba nguyên nhân tạo thành MTZ Kết tính tốn mơ chúng tơi cho thấy, chế độ động lực, đặc điểm phân bố cấu trúc muối xuất MTZ VCS Bạch Đằng có liên quan đến phân bố trầm tích lơ lửng (TTLL) phụ thuộc vào dao động mực nước biến đổi mùa dòng vật chất từ lục địa đưa Trong pha triều xuống, dịng trầm tích từ sơng phát triển mạnh phía ngồi biển Vào mùa mưa, vùng nước có hàm lượng TTLL khoảng 40mg/l xuất cách cửa sông khoảng 20-23km (đối với tầng mặt) 15-18km tầng đáy Vùng nước có hàm lượng trầm tích lơ lửng cao 150mg/l tầng mặt trải rộng từ cửa sông xa khoảng 15-18km, tầng đáy rộng 3km Ở tầng đáy tồn vùng nhỏ phía ngồi, MTZ, có hàm lượng trầm tích lơ lửng tương đối cao so với xung quanh (hình 4) Vào mùa khô, hàm lượng lưu lượng bùn cát sông đưa giảm nên phạm vi phân bố TTLL gần cửa sông với hàm lượng nhỏ, 80mg/l Sự chênh lệch phạm vi hàm lượng TTLL tầng đáy tầng mặt tương tự mùa mưa Tại tầng đáy xuất số MTZ nhỏ với hàm lượng TTLL khoảng 70mg/l Tây Nam đảo Cát Hải Hình Phân bố trầm tích lơ lửng tầng đáy VCS Bạch Đằng vào mùa mưa (triều lên) Trong pha triều lên, khối nước có hàm lượng TTLL cao từ lục địa tiếp tục đưa khối nước biển với độ mặn cao di chuyển phía lục địa Sự tương tác hai khối nước khiến phạm vi ảnh hưởng khối nước có hàm lượng TTLL cao bị thu hẹp đáng kể so với nước ròng Do ảnh hưởng dịng triều có hướng Đơng Bắc, phần lượng trầm tích phía Tây Nam bán đảo Đồ Sơn vận chuyển lên số vùng MTZ nhỏ thấy xuất tầng đáy VCS Bạch Đằng Lưu lượng hàm lượng TTLL giảm đáng kể mùa khô pha triều không thấy dịng trầm tích bổ sung cho VCS Bạch Đằng từ phía Tây Nam bán đảo Đồ Sơn Ở thời điểm nước lớn, ảnh hưởng khối nước biển vào lục địa lớn Vào mùa mưa, vùng nước có hàm lượng TTLL cao từ sơng bị dồn sát phía Tây Nam đảo Cát Hải Một phần TTLL từ phía Tây Nam bán đảo Đồ Sơn đẩy lên phía Đơng Bắc (hình 5) Các nêm mặn lấn mạnh tương tác với dòng TTLL dòng chảy sơng đẩy dồn góp phần tạo nên số MTZ, nhận thấy qua phân bố TTLL tầng đáy Vào mùa khơ, dịng TTLL từ lục địa đưa khơng lớn vào pha triều lên, ảnh hưởng khối nước biển chiếm ưu tuyệt đối tương tác lục địa - biển VCS Các khối nước biển với hàm lượng TTLL nhỏ tiến sát vào cửa sông, đồng thời không thấy vận chuyển TTLL vào VCS Bạch Đằng từ Tây Nam Đồ Sơn không xuất MTZ mùa mưa Trong pha triều xuống vào mùa mưa, phạm vi vùng nước sơng có hàm lượng TTLL cao mở rộng phía biển Ở lớp nước tầng đáy xuất vùng MTZ TTLL từ VCS Bạch Đằng di chuyển phía đơng nhiều di chuyển xuống phía Tây Nam bán đảo Đồ Sơn nhỏ Trong mùa khô, tải lượng TTLL từ lục địa đưa nhỏ, nên dù dòng triều chảy xuống, khối nước có hàm lượng TTLL cao khơng mở rộng phía ngồi theo chiều nằm ngang TTLL VCS rõ xuất MTZ Hình5 Phân bố trầm tích lơ lửng tầng đáy VCS Bạch Đằng vào mùa mưa (nước lớn) Trắc diện trầm tích lơ lửng vùng đục VCS Bạch Đằng Các kết tính tốn cho thấy trắc diện TTLL theo chiều thẳng đứng VCS Bạch Đằng thể rõ tương tác sông - biển biến đổi theo pha dao động mực nước triều biến động mùa hàm lượng lượng TTLL 2.1 Biến động mùa mưa Trong kỳ triều kém, biên độ dao động triều nhỏ nên ảnh hưởng khối nước biển có hàm lượng TTLL nhỏ đến vùng cửa sơng khơng lớn Vào pha triều xuống kỳ triều kém, khối nước sơng có hàm lượng TTLL cao phát triển phía ngồi, dù khơng lớn pha triều lên Trên mặt cắt I II, khối nước có hàm lượng TTLL 100mg/l mở rộng phía ngồi khoảng 5-10km Đáng ý mặt cắt I thấy xuất vùng nhỏ cách cửa sông khoảng 10-12km phía ngồi có hàm lượng TTLL lớn, khoảng 170mg/l Đó MTZ Trên mặt cắt II thấy MTZ tương tự, phạm vi nhỏ Do lưu lượng TTLL từ cửa Lạch Tray không lớn nên mặt cắt III không thấy xuất MTZ Phân bố TTLL mặt cắt IV V khơng có đáng ý Trong pha triều lên kỳ triều kém, khối nước biển dồn vùng nước có hàm lượng TTLL cao phía cửa sơng đường đẳng hàm lượng TTLL phía tiếp giáp với khối nước biển dày Phân bố TTLL mặt cắt thể xu hướng chung giảm dần hàm lượng phía ngồi Riêng mặt cắt V phía ngồi cửa Văn Úc, có vùng hàm lượng TTLL cao nằm cách cửa sông khoảng 10-12km Đây MTZ khu vực Trong kỳ triều cường, tương tác sông-biển diễn mạnh mẽ hơn, tải lượng trầm tích lơ lửng từ sông biên độ dao động mực nước triều lớn Vào pha triều xuống, ảnh hưởng khối nước có hàm lượng TTLL cao lớn nhất, tương tác hai khối nước sông-biển thể trắc diện TTLL vị trí dao động cách bờ khoảng 1522km Trên mặt cắt thấy xuất vùng có hàm lượng cao khu vực khác Các trắc diện TTLL I II giống nhau, có vùng MTZ xuất cách bờ khoảng 8-15km với qui mơ khác Ở phía ngồi MTZ này, hàm lượng TTLL giảm mạnh theo hướng biển, lại đồng theo chiều sâu Tải lượng TTLL từ sông Lạch Tray nhỏ, nên mặt cắt III, MTZ xuất sát cửa sông Hàm lượng TTLL mặt cắt biến động giảm mạnh theo hướng biển biến động theo độ sâu Mặt cắt IV nằm xa cửa sông chịu tác động đồng thời khối nước từ biển từ phía đơng bắc tây nam bán đảo Đồ Sơn nên hàm lượng TTLL biến động phức tạp Trên mặt cắt xuất số MTZ tầng mặt tầng đáy hàm lượng TTLL nước nhỏ Có thể độ muối cao tầng đáy phía ngồi biển khiến cho dịng TTLL khơng lớn khó xâm nhập Tại mặt cắt V, ảnh hưởng mạnh dịng TTLL từ sơng Văn Úc nên xuất MTZ phía ngồi, cách bờ khoảng 10km Xung quanh MTZ này, hàm lượng TTLL biến đổi nhanh theo chiều ngang không thay đổi theo độ sâu Trong pha triều lên kỳ triều cường, tác động dòng triều khối nước từ biển hướng vào bờ, vùng nước có hàm lượng TTLL cao bị đẩy trở lại vùng cửa sông Tuy nhiên, tải lượng TTLL từ sông đưa lớn lên tồn vùng nước có hàm lượng TTLL cao phía ngồi cửa sơng Sự tương tác sơng-biển pha triều tạo thành vùng nước có hàm lượng TTLL cao so với khu vực xung quanh Tại mặt cắt I xuất vùng MTZ nhỏ sát phía cửa sơng, mặt cắt thấy phân bố đồng hàm lượng TTLL theo độ sâu Trong dọc theo mặt cắt hướng biển, hàm lượng TTLL giảm dần, vị trí cách bờ khoảng 15km cịn xuất thêm vùng có hàm lượng TTLL khoảng 90-110mg/l Trên mặt cắt II pha triều xuất MTZ phía ngồi cửa sơng, cách bờ khoảng 3-5km, với hàm lượng TTLL 140mg/l vùng khác có hàm lượng TTLL tương đối cao khoảng 80-110mg/l, cách bờ khoảng 12-15km Trên mặt cắt III phía ngồi cửa sơng Lạch Tray, ta thấy rõ MTZ cách bờ khoảng 12-15km hàm lượng TTLL xung quanh biến động mạnh theo chiều ngang (hình 6), biến động theo độ sâu Tại mặt cắt IV, MTZ bị đẩy dần khoảng 1.5-2km phía bờ từ vị trí pha triều xuống lưỡi nước mặn phía nguyên nhân khiến hàm lượng TTLL lớp nước phía lớn so với lớp nước tầng đáy Xu hướng biến động chung TTLL mặt cắt vị trí MTZ bị đẩy lùi phía cửa sơng, cách bờ khoảng 8-10km, hàm lượng TTLL quanh MTZ biến đổi đồng theo độ sâu ĐN TB Hình Trắc diện trầm tích lơ lửng mặt cắt III - mùa mưa Pha triều lên - nước lớn khoảng thời gian ảnh hưởng lớn biển VCS Phân tích mặt cắt TTLL vào thời gian cho thấy tương tác sông-biển với ưu khối nước từ biển vào, nên khối nước sơng có hàm lượng TTLL cao bị đẩy sát vào cửa sông Tuy biến đổi khác nhau, MTZ xuất mặt cắt II, III V Trên mặt cắt II III, xuất vùng có hàm lượng TTLL cao so với xung quanh, thấp hàm lượng TTLL phía sơng Diễn biến phân bố theo trắc diện TTLL pha triều xuống có xu ngược lại với triều lên Trong tương tác sông-biển triều xuống chịu ảnh hưởng khối nước có hàm lượng TTLL cao, tăng cường thêm dịng chảy sơng hướng với dịng triều Hàm lượng TTLL biến động mạnh theo chiều ngang theo độ sâu Các MTZ xuất trở lại có xu hướng dịch chuyển dần phía ngồi 2.2 Biến động mùa khô Xu hướng biến động TTLL theo trắc diện mùa khô tương tự mùa mưa Tuy nhiên, lưu lượng nước sông, tải lượng hàm lượng TTLL giảm rõ rệt so với mùa mưa có tác động lớn đến tương tác sơng-biển hình thành MTZ VCS Bạch Đằng Trong pha triều xuống, ảnh hưởng khối nước từ biển vào nhỏ tạo điều kiện cho khối nước sông có hàm lượng TTLL cao xa phía ngồi Lưu lượng hàm lượng TTLL chung nhỏ Tại mặt cắt I II, TTLL khối nước sơng tầng mặt di chuyển phía ngồi nhanh hơn, tầng đáy di chuyển chậm cản trở khối nước biển có tỷ trọng cao Kết phía ngồi vùng tiếp giáp khối nước sông-biển, hàm lượng TTLL lớp nước phía có giá trị cao lớp nước tầng đáy Cũng hai mặt cắt này, xuất MTZ nhỏ phía ngồi cửa sơng, cách bờ khoảng 3-5km, có hàm lượng TTLL cao (trên 50mg/l) Trên mặt cắt III, hàm lượng TTLL nhỏ xuất vùng MTZ cách bờ khoảng 5-6km TTLL mặt cắt cịn lại khơng thể rõ xu hướng biến động hàm lượng TTLL (hình 7) Các mặt cắt cịn lại có hàm lượng TTLL nước nhỏ rõ quy luật biến đổi TB ĐN Hình Trắc diện trầm tích lơ lửng mặt cắt III - mùa khô Trong pha triều lên, khối nước biển đẩy khối nước sông với hàm lượng TTLL cao phía cửa sơng Khác với trường hợp triều lên mùa mưa, hàm lượng TTLL biến động mạnh theo hướng phía biển theo độ sâu Mặc dù hàm lượng TTLL nhỏ mặt cắt xuất MTZ với hàm lượng TTLL lớn so với xung quanh Tương tự mùa mưa, kỳ triều cường mùa khô tương tác khối nước biển nước sơng có hàm lượng TTLL cao diễn cửa sơng bị đẩy sâu vào sông Mặc dù hàm lượng TTLL nhỏ, mặt cắt III V thấy xuất vùng nhỏ có hàm lượng TTLL lớn xung quanh IV KẾT LUẬN Một kết tương tác lục địa - biển VCS Bạch Đằng xuất MTZ phía ngồi số cửa sơng Quy mơ vị trí chúng ln biến đổi theo vai trò ảnh hưởng dòng vật chất từ lục địa dao động mực nước thuỷ triều Các MTZ có vị trí xa cách bờ khoảng 15km vào thời điểm nước rịng mùa mưa Vào mùa khơ chúng có qui mơ nhỏ vị trí xuất gần bờ Các MTZ di chuyển dọc theo mặt cắt hướng từ cửa sơng phía ngồi, nhiên, chúng tạo bồi lắng vị trí mặt cắt mơi trường trầm tích khác biệt nhiều so với xung quanh Các MTZ nơi tập trung chất dinh dưỡng chất nhiễm từ lục địa có ý nghĩa quan trọng khía cạnh môi trường sinh thái Đây vấn đề có ý nghĩa thực tiễn, cần mở rộng phạm vi cho tồn vùng ven bờ châu thổ sơng Hồng TÀI LIỆU THAM KHẢO Bowden, K.F., 1984 Turbulence and mixing in estuaries In: Kennedy, V.S (Ed.), The Estuary as a Filter Acad Press, Orlando, pp 15-26 Charlotte harbor national estuary program, 2009 CHNEP Citizens Advisory Committee Meeting, CAC 7-15-09, pg 151 http://www.chnep.org Cauwet, G., Mackenzie, F.T., 1993 Carbon inputs and distribution in estuaries of turbid rivers: the Yang Tze and Yellow rivers (China) Marine Chemistry 43, 235-246 Fishes T R., 1991 Phytoplankton nutrients and turbiding in the Chaspeak, Delaware and Husdson esturies, 16(1) 104-112 Gordeev, V.V., Shevchenko, V.P., 1995 Chemical composition of suspended sediments in the Lena river and its mixing zone In: Kassens, H., et al (Ed.), Berichte zur Polar- und Meeresforschung AWI Bremerhaven, pp 154-169 Hamblin, P.F., 1989 Observations and model of sediment transport near the turbidity maximum of the upper Saint Lawrence estuary Journal of Geophysical Research 94 (C10): 14419-14428 Jiufa, L., Chen, Z., 1998 Sediment resuspension and implictions for turbidity maximum in the Changjiang Estuary Marine Geology 148: 117-124 John L L, 1993 Estuarine front How important are they? Estuarine , 16(1): 1-11 Johnson D.R., H E Garcia, and T P Boyer, 2005 World Ocean Database 2005:Tutorial, NODC Internal Report 19, U.S Government Printing Office, Wash., D.C., 18 pp 10 Julian J D, Jean - Claude D, Grant R I, et al., 1989 Abudance of rainbow smelt in relation to the maximum turbidity zone and asscociate macroplanktonic fauna of the middle St Lawrence Estuary Estuarine , 12 (2): 66-81 11 Kistner, D.A., Pettigrew, N.R., 2001 A variable turbidity maximum in the Kennebec Estuary Estuaries 24: 680- 687 12 Lefebvre J.P , S Ouillon, Vu Duy Vinh, R Arfi, Do Trong Binh, J.Y Panché, X Mari, J.P Torréton, 2009 Hydro-sedimentary behaviour of the Cam and Bach Dang estuaries flowing into the Haiphong bay Physical Part of the HAIPHONG project (funded by EC2CO) 13 Lin, J., Kuo, A.Y., 2003 A model study of turbidity maxima in the York River Estuary, Virginia Estuaries 26: 1269-1280 14 Lisitsyn, A.P., 1995 The marginal filter of the ocean Oceanology 34 (5), 671-682 (engl transl.) 15 Lucotte, M., 1989 Organic carbon isotope ratios and implications for the maximum turbidity zone of the St Lawrence upper estuary Estuarine, Coastal and Shelf Science 29: 293-304 10 16 Lyard F., F Lefevre, T Letellier, and O Francis, 2006 Modelling the global ocean tides: modern insights from FES2004 Ocean Dynamics, 56: 394-415: 2006 17 Sheng Fangting, 1995 A new reconization of estuaries tubidity maximum in china Advance in Earth Sciences, 10(2): 210-213 18 Shi, Z., Kirby, R., 2003 Observations of fine suspended sediment processes in the turbidity maximum at the north passage of the Changjiang Estuary, China Journal of Coastal Research 19: 529-540 19 Simpson, J.E., 1997 Gravity Currents in the Environment and the Laboratory, second ed Cambridge University Press, Cambridge, 244 pp 20 Wang Hui, 1998 Several basic problem of marine ecosystem modeling Oceanologia et Limnologia Sinica, 29(2): 341-346 21 WL|Delft Hydraulics, 1999 Delft3D-FLOW User Manual Version 3.05, Delft3DWaq User Manual Version 3.01, Delft3D-Part User Manual Version 1.0 WL| Delft Hydraulics, Delft, Netherlands 22 Wolanski, E., Spagnol, S., 2003 Dynamics of the turbidity maximum in King Sound, tropical Western Australia Estuarine, Coastal and Shelf Science 56: 877-890 APLICATION NUMERICAL MODEL TO STUDY ON MAXIMUM TURBIDITY ZONES IN BACH DANG ESTUARY VU DUY VINH, TRAN DUC THANH Summary: Recent researches on dynamics of land- sea interaction in estuaries have discovered a zone with very high turbidity or maximum turbidity zone (MTZ) This zone is rich in the nutrients concerning abundance of marine biological resources, and high in potential of pollutant accumulation The studied using Delf-3D Model in this paper has discovered MTZ in Bach Dang Estuary Its scale and location are changed by river discharge in seasons, and tidal oscillation The furthest location of MTZ in this estuary is about 15km from shoreline during the ebb tide in rainy season In the dry season, it is smaller in scale and moves closed to the coast Ngày nhận bài: 22 - 12 - 2011 Người nhận xét: TS Trần Đình Lân 11