Nghiên cứu ứng dụng hệ thống mô hình thuỷ nhiệt động lực sinh thái biển phục vụ nghiên cứu và quản lý tài nguyên biển vùng ven bờ việt nam

BỘ KHOA HỌC VÀ CÔNG NGHỆ VIỆN KHOA HỌC VÀ CƠNG NGHỆ VIỆT NAM CHƯƠNG TRÌNH HỢP TÁC VỀ KHOA HỌC VÀ CÔNG NGHỆ THEO NGHỊ ĐỊNH THƯ VIỆT NAM- VƯƠNG QUỐC BỈ GIAI ĐOẠN 2010-2011 BÁO CÁO TỔNG HỢP KẾT QUẢ KHOA HỌC CÔNG NGHỆ NHIỆM VỤ NGHIÊN CỨU ỨNG DỤNG HỆ THỐNG MƠ HÌNH THỦY NHIỆT ĐỘNG LỰC-SINH THÁI BIỂN PHỤC VỤ NGHIÊN CỨU VÀ QUẢN LÝ TÀI NGUYÊN BIỂN VÙNG VEN BỜ VIỆT NAM Cơ quan chủ trì: Chủ trì nhiệm vụ: Viện Tài ngun Mơi trường biển (Viện Khoa học Công nghệ Việt Nam) Ths Vũ Duy Vĩnh Hải Phòng, 2013 BỘ KHOA HỌC VÀ CÔNG NGHỆ VIỆN KHOA HỌC VÀ CÔNG NGHỆ VIỆT NAM CHƯƠNG TRÌNH HỢP TÁC VỀ KHOA HỌC VÀ CƠNG NGHỆ THEO NGHỊ ĐỊNH THƯ VIỆT NAM- VƯƠNG QUỐC BỈ GIAI ĐOẠN 2010-2011 BÁO CÁO TỔNG HỢP KẾT QUẢ KHOA HỌC CÔNG NGHỆ NHIỆM VỤ NGHIÊN CỨU ỨNG DỤNG HỆ THỐNG MƠ HÌNH THỦY NHIỆT ĐỘNG LỰC-SINH THÁI BIỂN PHỤC VỤ NGHIÊN CỨU VÀ QUẢN LÝ TÀI NGUYÊN BIỂN VÙNG VEN BỜ VIỆT NAM Cơ quan chủ trì nhiệm vụ Chủ trì nhiệm vụ Vũ Duy Vĩnh Ban chủ nhiệm chương trình Bộ Khoa học Cơng nghệ Hải Phịng, 2013 Báo cáo tổng kết Nhiệm vụ: Nghiên cứu ứng dụng hệ thống mơ hình thủy nhiệt động lực-sinh thái biển phục vụ nghiên cứu quản lý tài nguyên biển vùng ven bờ Việt Nam DANH SÁCH THÀNH VIÊN THAM GIA THÀNH VIÊN VIỆT NAM TT Họ tên Ths Phạm Hải An KS Nguyễn Thị Kim Anh Đơn vị công tác Nhiệm vụ đề tài Viện Tài nguyên Môi trường biển Tham gia nt Thư ký HC nt Đồng Chủ trì TS Đỗ Trọng Bình CN Nguyễn Minh Hải nt Tham gia CN Lê Đức Cường nt Tham gia Ths Chu Thế Cường nt Tham gia Ths Nguyễn Thị Minh Huyền nt Tham gia KS Vũ Thị Lựu nt Tham gia Ths Nguyễn Trung Kiên Đại học KHTN Tham gia 10 Ths Dương Thanh Nghị Viện Tài nguyên Môi trường biển Tham gia 11 Ths Nguyễn Thị Thu nt Tham gia 12 Ths Cao Thị Thu Trang nt Tham gia 13 Ths Trần Anh Tú nt Thư ký KH 14 Ths Vũ Duy Vĩnh nt Chủ trì Đơn vị cơng tác Nhiệm vụ đề tài Tham gia Tham gia THÀNH VIÊN BỈ TT Họ tên TS Katrijn Baetens Ban quản lý mơ hình tốn học Biển bắc, Viện Khoa học Tự nhiên Hồng gia Vương quốc Bỉ TS Patrick Luyten nt ViÖn Tài nguyên Môi trờng biển 246 - Đà Nẵng, Ngô Quyền, Hải Phòng iii Bỏo cỏo tng kt Nhim vụ: Nghiên cứu ứng dụng hệ thống mơ hình thủy nhiệt động lực-sinh thái biển phục vụ nghiên cứu quản lý tài nguyên biển vùng ven bờ Việt Nam VIỆN KHOA HỌC VÀ CƠNG NGHỆ VIỆT NAM CỘNG HỒ XÃ HỘI CHỦ NGHĨA VIỆT NAM Độc lập - Tự - Hạnh phúc Viện Tài nguyên Môi trường biển Hải Phòng, ngày tháng năm 2013 BÁO CÁO THỐNG KÊ KẾT QUẢ THỰC HIỆN NHIỆM VỤ NGHỊ ĐỊNH THƯ I THÔNG TIN CHUNG Tên nhiệm vụ: “Nghiên cứu ứng dụng hệ thống mơ hình thủy nhiệt động lực sinh thái biển phục vụ nghiên cứu quản lý tài nguyên biển vùng ven bờ Việt Nam” Thuộc: Chương trình hợp tác khoa học công nghệ Cộng hòa XHCN Việt Nam Vương quốc Bỉ Nghị định thư Việt Nam - Bỉ giai đoạn 2010-2011 theo Nghị định thư Chủ nhiệm nhiệm vụ: Họ tên: Vũ Duy Vĩnh Ngày, tháng, năm sinh: 18/11/1978 Nam/ Nữ: Nam Học hàm, học vị: Thạc sỹ Chức danh khoa học: Nghiên cứu viên Điện thoại: Tổ chức: 0313760605; Nhà riêng: Mobile: 0912799629 Fax: 0313761521; E-mail: vinhvd@imer.ac.vn Tên tổ chức công tác: Viện Tài nguyên Môi trường biển Địa tổ chức: 246 phố Đà Nẵng, quận Ngô Quyền, thành phố Hải Phòng Địa nhà riêng: Cầu Đen, Núi Đối, Huyện Kiến Thuỵ, thành phố Hải Phòng Tổ chức chủ trì nhiệm vụ: Tên tổ chức chủ trì nhiệm vụ: Viện Tài ngun Mơi trường biển Điện thoại: 0313761523; Fax: 0313761521 E-mail: imervn@imer.ac.vn Website: http//www.imer.ac.vn Địa tổ chức: 246 phố Đà Nẵng, quận Ngô Quyền, thành phố Hải Phòng Họ tên thủ trưởng t chc: TRN C THNH Viện Tài nguyên Môi trờng biển 246 - Đà Nẵng, Ngô Quyền, Hải Phòng iv Báo cáo tổng kết Nhiệm vụ: Nghiên cứu ứng dụng hệ thống mơ hình thủy nhiệt động lực-sinh thái biển phục vụ nghiên cứu quản lý tài nguyên biển vùng ven bờ Việt Nam Số tài khoản: 3711, cấp NS Ngân hàng: Kho bạc Nhà nước Hải Phòng Tên quan chủ quản đề tài: Viện Khoa học Cơng nghệ Việt Nam II TÌNH HÌNH THỰC HIỆN Thời gian thực nhiệm vụ - Theo Hợp đồng ký kết: từ tháng năm 2010 đến tháng 12 năm 2011 - Thực tế thực hiện: từ tháng năm 2010 đến tháng năm 2012 Kinh phí sử dụng kinh phí: a) Tổng số kinh phí thực hiện: 650 triệu đồng, đó: + Kính phí hỗ trợ từ SNKH: 650 triệu đồng + Kinh phí từ nguồn khác: tr.đ + Tỷ lệ kinh phí thu hồi dự án (nếu có): b) Tình hình cấp sử dụng kinh phí từ nguồn SNKH (tính đến 31/12/2011) Thời gian (Tháng, năm) Kinh phí (Tr.đ) Thời gian (Tháng, năm) Kinh phí (Tr.đ) Ghi (Số đề nghị toán) 2010 900 2010 900 900 2011 750 2011 750 750 Số TT Theo kế hoạch Thực tế đạt c) Kết sử dụng kinh phí theo khoản chi: Đơn vị tính: Triệu đồng Số TT Nội dung khoản chi Theo kế hoạch Tổng SNKH Nguồn khác 600 600 600 73 73 80,9 80,9 487 487 487 487 Tổng SNKH Trả công lao động (khoa học, phổ thông) 600 Nguyên, vật liệu, lượng Thiết bị, máy móc Xây dựng, sửa chữa nhỏ Chi khác Tổng cộng Thực tế đạt Nguồn khác 490 490 650 650 Viện Tài nguyên Môi trờng biển 246 - Đà Nẵng, Ngô Quyền, Hải Phòng 482,1 482,1 650 650 v Báo cáo tổng kết Nhiệm vụ: Nghiên cứu ứng dụng hệ thống mơ hình thủy nhiệt động lực-sinh thái biển phục vụ nghiên cứu quản lý tài nguyên biển vùng ven bờ Việt Nam Các văn hành trình thực nhiệm vụ (Liệt kê định, văn quan quản lý từ công đoạn xác định nhiệm vụ, xét chọn, phê duyệt kinh phí, hợp đồng, điều chỉnh (thời gian, nội dung, kinh phí thực có); văn tổ chức chủ trì đề tài, dự án (đơn, kiến nghị điều chỉnh có) Số Số, thời gian ban hành văn TT Tên văn Ghi Số 61/2010/HĐ-NĐT Hợp đồng thực nhiệm vụ khoa tháng năm 2010 học công nghệ theo nghị định thư Bộ Khoa học Công nghệ, Viện Khoa học Công nghệ Việt Nam Thời gian gian thực nhiệm vụ từ 1/1/2010 đến 31/12/2011 Số 70/QĐ-TMB ngày Về việc thành lập BCN Chủ nhiệm 4/3/2010 Viện nhiệm vụ trưởng Viện Tài nguyên Môi trường biển Ban Chủ trì Nhiệm vụ Số 2681/BKHCNVề việc thay đổi Chủ nhiệm nhiệm Đồng ý thay XHTN ngày 01/11/2010 vụ HTQT theo Nghị định thư Việt-Bỉ đổi chủ trì Bộ Khoa học Nhiệm vụ Cơng nghệ Số 438/TMB ngày Về việc xin phép tổ chức hội thảo 19/12/2011 Viện trưởng Viện Tài nguyên Môi trường biển Về việc xin gia hạn thực nhiệm Xin gia hạn Số 436/TMB ngày vụ thực 23/11/2011 Viện trưởng Viện Tài nguyên Nhiệm vụ Môi trường biển đến tháng 6/2012 Số 3359/BKHCNVề việc gia hạn thời gian thực XHTN ngày 26/12/2011 nhiệm vụ nghị định thư với Bỉ Bộ Khoa hc v Cụng ngh Viện Tài nguyên Môi trờng biển 246 - Đà Nẵng, Ngô Quyền, Hải Phòng Xin phép tổ chức Hội thảo quốc tế đánh giá kết thực nhiệm vụ Đồng ý gia hạn thực Nhiệm vụ đến tháng 6/2012 vi Báo cáo tổng kết Nhiệm vụ: Nghiên cứu ứng dụng hệ thống mơ hình thủy nhiệt động lực-sinh thái biển phục vụ nghiên cứu quản lý tài nguyên biển vùng ven bờ Việt Nam Tổ chức phối hợp thực nhiệm vụ Tên tổ chức STT đăng ký theo Thuyết minh Trường Đại học Khoa học tự nhiên (Đại học Quốc gia Hà Nội) Tên tổ chức tham gia thực Trường Đại học Khoa học tự nhiên (Đại học Quốc gia Hà Nội) Các sở KH&CN, sở TN&MT, tỉnh Quảng Ninh, Hải Phịng, Thái Bình số tỉnh có liên quan Sản phẩm Nội dung Ghi chủ yếu đạt tham gia chủ yếu chú* Phối hợp tham gia thực hiện, tham gia hội thảo, đóng góp ý kiến việc xây dựng kịch mô Tham gia vào mạng lưới người sử dụng mơ hình Tham gia xây dựng hệ thống bó máy tính Tham gia hội thảo, hội nghị đóng góp ý kiến việc xây dựng kịch mô Tham gia vào mạng lưới người sử dụng mơ hình kết mơ hình Cá nhân tham gia thực nhiệm vụ STT Tên cá nhân đăng ký theo Thuyết minh Tên cá nhân tham gia thực Ths Phạm Hải An Ths Phạm Hải An TS Trương Văn Ths Chu Th Viện Tài nguyên Môi trờng biển 246 - Đà Nẵng, Ngô Quyền, Hải Phòng Ni dung tham gia Cài đặt thiết lập hệ thống bó máy tính Biến đổi khí Sản phẩm chủ yếu đạt Hệ thống bó máy tính cài đặt để chạy mơ hình Báo cáo Ghi chú* vii Báo cáo tổng kết Nhiệm vụ: Nghiên cứu ứng dụng hệ thống mô hình thủy nhiệt động lực-sinh thái biển phục vụ nghiên cứu quản lý tài nguyên biển vùng ven bờ Việt Nam Bốn* CN Lê Đức Cường TS Đỗ Trọng Bình 10 11 12 13 Cường CN Lê Đức Cường TS Đỗ Trọng Bình hậu tác động đến sinh thái Điều kiện tự nhiên, KTXH, MT khu vực nghiên cứu Mơ hình thủy động lực- sinh thái Thu thập, xử lý số liệu chuyên đề Báo cáo chuyên đề Số liệu phân tích xử lý Số liệu Các muối dinh KS Vũ Thị Lựu TS Lưu Văn Diệu* phân tích dưỡng xử lý CN Nguyễn Minh Thu thập, xử lý Số liệu CN Nguyễn Minh Hải số liệu phân tích Hải xử lý Sinh vật phù Số liệu Ths Nguyễn Ths Nguyễn Minh du phân tích Minh Huyền Huyền xử lý Hệ thống bó máy Cài đặt thiết tính PGS.TS Trần Đức Ths Nguyễn Trung lập hệ thống Kiên Thạnh* cài đặt để bó máy tính chạy mơ hình Số liệu Ths Nguyễn Thị Sinh vật phù Ths Nguyễn Thị du phân tích Thu Thu xử lý Ths Cao Thu Các muối dinh Số liệu Ths Cao Thu Trang dưỡng phân tích Trang xử lý Mơ hình thủy Báo cáo, động lực, thu tài liệu Ths Trần Anh Tú Ths Trần Anh Tú thập xử lý tài liệu, tổ chức khảo sát GS TS Đinh Văn Ths Dương Các muối dinh Số liệu TS Trần Đình Lân* KS Nguyễn Thị Kim Anh Viện Tài nguyên Môi trờng biển 246 - Đà Nẵng, Ngô Quyền, Hải Phòng viii Bỏo cỏo tng kết Nhiệm vụ: Nghiên cứu ứng dụng hệ thống mô hình thủy nhiệt động lực-sinh thái biển phục vụ nghiên cứu quản lý tài nguyên biển vùng ven bờ Việt Nam 14 Ưu* Thanh Nghị dưỡng Ths Vũ Duy Vĩnh Mơ hình thủy Ths Vũ Duy Vĩnh động lực- sinh thái phân tích xử lý Báo cáo chuyên đề * đăng ký tham gia thuyết minh lý bận nên khơng tham gia Tình hình hợp tác quốc tế Theo kế hoạch Số Thực tế đạt (Nội dung, thời gian, kinh phí, địa điểm, tên tổ chức hợp tác, số đoàn, số lượng người tham gia ) (Nội dung, thời gian, kinh phí, địa điểm, tên tổ chức hợp tác, số đoàn, số lượng người tham gia ) Đoàn vào hội thảo triển khai kế hoạch phân cơng nhiệm vụ Đồn gồm 02 người từ Ban quản lý mơ hình tốn học biển Bắc (Viện KHTN Hoàng gia vương quốc Bỉ): TS Patrick Luyten, TS Katrijn Baetens Đoàn vào hội thảo triển khai kế hoạch phân cơng nhiệm vụ Đồn gồm 02 người từ Ban quản lý mơ hình tốn học biển Bắc (Viện KHTN Hoàng gia vương quốc Bỉ): TS Patrick Luyten, TS Katrijn Baetens Đoàn vào tập huấn sử dụng mơ hình trao đổi khoa học Đồn gồm 02 người từ Ban quản lý mơ hình tốn học biển Bắc (Viện KHTN Hoàng gia vương quốc Bỉ): TS Patrick Luyten, TS Katrijn Baetens Đoàn vào tập huấn sử dụng mơ hình trao đổi khoa học Đồn gồm 02 người từ Ban quản lý mơ hình tốn học biển Bắc (Viện KHTN Hoàng gia vương quốc Bỉ): TS Patrick Luyten, TS Katrijn Baetens Đoàn vào tham gia hội thảo tổng kết nhiệm vụ, trao đổi khoa học bàn hướng hợp tác tới: Đoàn gồm 02 người từ Ban quản lý mơ hình tốn học biển Bắc (Viện KHTN Hoàng gia vương quốc Bỉ): TS Patrick Luyten, TS Katrijn Baetens Đoàn vào tham gia hội thảo tổng kết nhiệm vụ, trao đổi khoa học bàn hướng hợp tác tới: Đoàn gồm 02 người từ Ban quản lý mơ hình tốn học biển Bắc (Viện KHTN Hoàng gia vương quốc Bỉ): TS Patrick Luyten, TS Katrijn Baetens Đồn tập huấn sử dụng mơ Đồn tập huấn sử dụng mơ hình Ban quản lý mơ hình hình Ban quản lý mơ hình TT Viện Tài nguyên Môi trờng biển 246 - Đà Nẵng, Ngô Quyền, Hải Phòng Ghi chỳ* ix Báo cáo tổng kết Nhiệm vụ: Nghiên cứu ứng dụng hệ thống mơ hình thủy nhiệt động lực-sinh thái biển phục vụ nghiên cứu quản lý tài nguyên biển vùng ven bờ Việt Nam toán học biển Bắc (Viện KHTN Hoàng gia vương quốc Bỉ) Đoàn gồm 02 thành viên từ Viện Tài nguyên Môi trường biển: Vũ Duy Vĩnh; Chu Thế Cường toán học biển Bắc (Viện KHTN Hoàng gia vương quốc Bỉ) Đoàn gồm 02 thành viên từ Viện Tài nguyên Môi trường biển Vũ Duy Vĩnh; Chu Thế Cường Đoàn tham dự hội thảo trình bày kết tại Ban quản lý mơ hình tốn học biển Bắc Đồn gồm 02 thành viên từ Viện Tài nguyên Môi trường biển: Vũ Duy Vĩnh; Trần Anh Tú Đoàn tham dự hội thảo trình bày kết tại Ban quản lý mơ hình tốn học biển Bắc Đoàn gồm 02 thành viên từ Viện Tài nguyên Môi trường biển Vũ Duy Vĩnh; Trần Anh Tú Tình hình tổ chức hội thảo, hội nghị Số TT Theo kế hoạch (Nội dung, thời gian, kinh phí, địa điểm ) Thực tế đạt (Nội dung, thời gian, địa điểm ) kinh phí, Hội thảo triển khai kế hoạch phân công nhiệm vụ Tại Viện Tài nguyên Môi trường biển Hội thảo triển khai kế hoạch phân cơng nhiệm vụ (có tham gia 02 người phía đối tác Bỉ) Tại Viện Tài nguyên Môi trường biển Hội thảo tập huấn sử dụng mơ hình Hội thảo đánh giá kết kỳ tập huấn sử dụng mơ hình Với tham dự 19 học viên đến từ quan nghiên cứu, trường ĐH nước 02 giảng viên từ Bỉ Hội thảo tổng kết nhiệm vụ Hội thảo tổng kết nhiệm vụ có tham dự 02 chuyên gia Bỉ đại diện sử KHCN, TNMT địa phương: Hải Phòng, Quảng Ninh, Thái Bình, Nam Định Ghi chú* Tóm tắt nội dung, công việc chủ yếu (Nêu mục 15 thuyết minh, không bao gồm: Hội thảo khoa học, iu tra kho sỏt Viện Tài nguyên Môi trờng biển 246 - Đà Nẵng, Ngô Quyền, Hải Phòng x The influence of wind on the hydrodynamics of the Red River Delta, North Vietnam Fig Contour plots of surface salinity distributions at different river discharges, wind directions and wind velocities Low and high river discharges refer to the river discharges of the dry season and rainy season, the values are as listed in table The arrows show the wind direction a) a) b) Fig 11 Surface salinity during the rainy season at the balat transect during neap tide (a) and spring tide (b) Dashed lines indicate high tide, solid lines indicate low tide b) Fig 10 Salinity profile of the Balat transect during the rainy season with no surface wind (a) and a surface wind of 3.5 m s−1 (b) 3.3 Analysis of model resulst with a typhoon wind profile Discussion The model is validated through surface elevations measured at tidal station, so it is rather concise The fit has an RM SE value of 0.35 m for the rainy and 0.31 m for the dry season Though the model lacked skill in properly simulating the amplitudes, the phases were captured well The wind measurements were variable in time, but not in space The measurements at Hon Dau station were generalized for the whole domain Hon Dau station is lying lee-ward of land for both the dry and the rainy season and may have lower or higher (due to tunneling effects) wind velocities than the rest of the domain This may cause the discrepancy in measured and modeled sea elevation The measurements show a certain skewness of the amplitudes, which is not captured by the model This can be caused by grid resolution of the model, which did not capture many of the small islands at Halong Bay The accuracy of the model grid resolution was tested by implementing a nesting procedure around Hon Dau station The nested grid has a resolution of 0.0021 ◦ , this grid captures a lot more of the islands of Halong Bay in the North East of the domain This nesting procedure did however not change the results significantly The skewness can also be caused by sedimentation processes, less accurate tidal boundary conditions or some measurement anomalies of the tidal elevation The hypothesis of the less accurate boundary conditions is supported by the work of Fang et al (1999) Fang et al (1999) showed that for the region the resonance of the M2 and the O1 tide result in asymmetrical flows This is not picked up by the model boundary input Sedimentation and land flooding can not be studied with the current code In the near future it will be possible to study the effect of inundation in more detail if the inundation module will be coupled to COHERENS (Heredia et al 2011) The validation with remote satellite data did not add extra value due to the high resolution of the added region More data are needed to give an in depth validation of the model According to the model results, wind and wind direction have a considerable impact on the hydrodynamics of the region Salinity is lower in the rainy season than in the dry season due to higher river discharges in the rainy season The formation and evolution of river plumes has been examined numerically by e.g Chao and Boicourt (1986), Chao (1988); Kourafalou et al (1996) in the absence of tides and by Ruddick et al (1995) who simulated the tidally modulated plume of the river Rhine The general picture emerging from these studies is that the fresh water released at the river mouth, first expands seawards and then turns anticyclonically (i.e to the right looking seawards in the northern hemisphere) Before reaching the coast again, the plume water deflects in the cyclonic sense forming a buoyancy driven coastal jet The plume deflection K Baetens et al in the case study is not very clear because the anticyclonic movement is modulated by the tidal current The phenomenon is clearer in the dry season since the river discharge is lower and a more thin salinity belt is formed along the coast line The SW wind of the rainy season balances for the plume deflection resulting in a river plume nearly perpendicular to the coast line The results of the transect in Fig 10 show that river water and sea water mix better when a South Western wind is applied The vertical residual data could not clearly demonstrate that a SW wind is a favourable upwelling wind It is demonstrated that this wind balances for the anticyclonic force on the plume and the influence of the Ekman drift becomes visible (Fig 6) The wind enhances the relative importance of offshore dispersion of the salinity plume, this in contrast to a NE wind Fig shows that close to the coast, the wind direction has a strong influence on the salinity circulation of the region, this is supported by the shape of the plumes shown in Fig The vectors of the wind plot are making several turns of 360 ◦ during a short period, suggesting that the center of the typhoon passed the measurement stations.However, two station measurements are insufficient to recover a correct wind profile The profile of Hon Dau station was selected, since this station is closer to the coast Typhoons are often accompagnied by strong rainfall and a changed river discharge profile No measurements were available, therefor the river discharge profile of the rainy season was taken as a boundary condition for this run The results of the wind sensitivity study can be extrapolated to extreme values of wind velocities and rapidly changing wind direction: This stresses even more the high sensitivity of the studied region to wind velocity and direction Conclusion A 3D hydrodynamical model for the coastal zone of the Red River delta area in Vietnam is presented The several river mouths are responsible for a characteristic salinity distribution in the modeled region Wind and wind direction are important for the hydrodynamic profile of the modeled domain According to the model results, the surface current profile close to the coast is sensitive to wind velocity changes, at the open sea, it is determined mainly by tidal forces A SW wind spreads the coastal salinity belt to the open sea A NE wind presses the salinity belt against the coast The model demonstrated more horizontal and vertical mixing of salinity when a SW wind was applied The typhoon study is an extreme application of the results gained in the wind sensitivity study The influence of wind on the hydrodynamics of the Red River Delta, North Vietnam The developed model showed that it can provide insights in the circulation regime of the region, which is important for good management Close collaboration with real concerns and strategy planning as explained in Mai et al (2009) and more detailed input and validation data are necessary References Chao SY (1988) River-forced estuarine plumes J Phys Oceanogr 18:72–88 Chao SY, Boicourt WC (1986) Onset of estuarine plumes J Phys Oceanogr 16:2137–2149 Fang G, Kwok YK, Yu K, Zhu Y (1999) Numerical simulation of principal tidal constituents in the south china sea, gulf of tonkin and gulf of thailand Cont Shelf Res 19(7):845–869 Hedstrom GW (1979) Nonreflecting boundary conditions for nonlinear hyperbolic systems J Comput Phys 30:222–237 Hellerman S, Rosenstein M (1983) Normal monthly wind stress over the world ocean with error estimates J Phys Oceanogr 13:1093–1104 Heredia M, Rocabado I, Hyde P, Luyten P, van Holland G (2011) Implementation of Inundation Schemes into the COHERENS Mode Fifth International Conference on Advanced Computational Methods in Engineering, Liege, Belgium Hirsch C (1990) Numerical computation of internal and external flows, vol 2: Computational methods for inviscid and viscous flows New York: Wiley Kourafalou VH, Oey LY, Wang JD, Lee TN (1996) The fate of river discharge on the continental shelf - modeling the river plume and the inner shelf coastal current J Geophys Res 101:3415–3434 Luyten P (2011) COHERENS - A Coupled Hydrodynamical-Ecological Model for Regional and Shelf Seas: User Documentation Version V2.0 RBINS-MUMM, Royal Belgian Institute of Natural Sciences, Brussels, Belgium Mai CV, Stive MJF, Van Gelder PHAJM (2009) Coastal protection strategies for the red river delta J Coastal Res 25(1):105–116 McDougall TJ, Jackett DG D R and Wright, Feistel R (2003) Accurate and computationally efficient algorithms for potential temperature and density of seawater J Atmos Oceanic Technol 20:730–741 Phillips NA (1957) A coordinate system having some special advantages for numerical forecasting Journal Meteor 14:184–185 Rød LP, Cooper CK (1987) A study of various open boundary conditions for wind-forced barotropic numerical ocean models Amsterdam: Elsevier Ruddick KG, Deleersnijder E, Luyten PJ, Ozer J (1995) Haline stratification in the rhine-meuse freshwater plume: A three-dimensional model sensitivity analysis Cont Shelf Res 15:1597–1630 Zu T, Gan J, Erofeeva S (2008) Numerical study of the tide and tidal dynamics in the south china sea Deep-Sea Res 55:137–154 ẢNH HƯỞNG CỦA GIĨ BỀ MẶT ĐẾN PHÂN BỐ ĐỘ MẶN VÀ HỒN LƯU VÙNG VEN BỜ CHÂU THỔ SÔNG HỒNG VŨ DUY VĨNH(1), KATRIJN BAETENS(2), PATRICK LUYTEN(2), TRẦN ANH TÚ(1), NGUYỄN THỊ KIM ANH(1) (1) Viện Tài nguyên Môi trường biển, 246 Đà Nẵng Ngơ Quyền, Hải Phịng Management Unit of the North Sea Mathematical Models, Royal Belgian Institute of Natural Sciences, Gulledelle 100, 1200 Brussels Belgium (2) TÓM TẮT: Bài viết trình bày kết áp dụng mơ hình vật lý- thủy động lực cho vùng ven bờ châu thổ sơng Hồng Đây mơ hình phát triển với mã nguồn mở COHERENS V2.0- mơ hình chiều áp dùng cho vùng ven bờ thềm lục địa dựa phương pháp phần tử hữu hạn Các kết tính tốn với 10 kịch khác cho thấy vai trò điều kiện gió kết hợp với thủy triều tải lượng nước sơng đến phân bố độ mặn hồn lưu ven bờ khu vực Theo trường gió mùa khơ làm tăng cường vận tốc dịng chảy dư xuống phía tây nam, tăng xâm nhập mặn vào vùng ven bờ gradient độ mặn theo phương thẳng đứng Trong vào mùa mưa, trường gió làm tăng cường vận chuyển khối nước từ sông phía ngồi, tăng phạm vi ảnh hưởng khối nước sông lớp nước bề mặt ảnh hưởng nước biển tầng đáy vào vùng ven bờ I MỞ ĐẦU Vùng ven bờ châu thổ sông Hồng (CTSH) có vị trí quan trọng phát triển kinh tế xã hội Việt Nam Với vị trí thuận lợi mối quan hệ với tam giác phát triển Hải Phòng- Hà Nội –Quảng Ninh nguồn tài nguyên biển phong phú, khu vực trở thành nơi phát triển kinh tế động nước ta Tuy nhiên phát triển kinh tế - xã hội gây sức ép lớn môi trường tự nhiên khu vực [17] Nghiên cứu môi trường biển khu vực thực từ nhiều góc độ khác nhau, có phương pháp tiếp cận từ mơ hình tốn học để nghiên cứu trình động lực nhằm tăng cường hiểu biết trình mối liên hệ với điều kiện môi trường khác Nghiên cứu điều kiện động lực môi trường vùng ven bờ châu thổ sông Hồng nhận quan tâm tác giả nước quy mô nhỏ [20, 21, 22, 23, 24] quy mô lớn [13, 19, 25] Bài viết đưa số kết nghiên cứu quy mô vừa ảnh hưởng trường gió bề mặt đến phân bố độ mặn hoàn lưu vùng ven bờ châu thổ sông Hồng thông qua việc sử dụng mô hình COHERENS II TÀI LIỆU VÀ PHƯƠNG PHÁP Khu vực nghiên cứu nằm vùng biển ven bờ phía tây Vịnh Bắc Bộ, hệ tọa độ 19.321.01 độ vĩ bắc 105.6 – 107.71 độ kinh đông Vùng biển có đặc điểm thủy triều mang tính chất nhật triều với biên độ lớn Độ dốc đáy biển nhỏ độ sâu lớn khoảng 40m Khu vực chịu ảnh hưởng mạnh khối nước từ hệ thống sơng Hồng-Thái Bình đưa tải lượng nước phân phối không năm mà chủ yếu tập trung vào tháng mùa mưa Khu vực chịu chi phối hệ thống gió mùa đơng bắc mùa khơ gió mùa đơng nam mùa mưa 1 Mơ hình COHERENS COHERENS (COupled Hydrodynamical Ecological model for REgioNal Shelf) mơ hình thủy động lực chiều cho vùng biển thềm lục địa [11] Mơ hình sử dụng hệ tọa độ khác Tuy nhiên nghiên cứu lựa chọn hệ tọa độ cầu Hệ tọa độ theo phương thẳng đứng mơ hình xác định theo hệ tọa độ σ [14] với 20 lớp Độ sâu mơ hình xử lý từ đồ địa hình cho vùng ven bờ số liệu từ sở liệu địa hình GEBCO [10] phía ngồi đồ địa hình tỷ lệ 1: 50000 ven bờ Lưới tính mơ hình cho vùng ven bờ châu thổ sông Hồng lưới chữ nhật với độ phân giải 0.01 độ Mơ hình gồm biên mở phía biển biên sơng Hình Độ sâu vị trí biên mở, biên sơng mơ hình Các phương trình sử dụng cho mơ hình thủy động lực khu vực phương trình liên tục, phương trình động lượng (momentum equations) phương trình liên tục cho độ mặn Các phương trình xuất phát từ phương pháp xấp xỉ Boussinesq với giả thiết áp suất thủy tĩnh, thành phần bình lưu vận tốc quay trái đất bỏ qua Mối liên hệ mật độ nước, độ muối nhiệt độ xác định phương trình trạng thái tham khảo theo phương pháp McDougall [9] Đây phương pháp tác giả khác ứng dụng mơ hình COHERENS cho kết xác phương trình trạng thái quốc tế sử dụng số mơ hình khác UNESCO 1980 [18] có hiệu tính tốn cao Điều kiện nhiệt – muối mơ hình biên mở sông sử dụng số liệu quan trắc, biên mở biển biển sử dụng số liệu đặc trưng trung bình tháng sở liệu Word Ocean Atlas 2009 (WOA09) với phương pháp nội suy tuyến tính [26] Thuật tốn mơ hình số cho gradient áp suất thủy tĩnh (baroclinic pressure gradient) hệ tọa độ σ truyền thống (bậc hai) Sơ đồ khép kín rối k sử dụng cho mơ hình Sơ đồ RANS (Reynolds Averaged Navier-Stokes) mơ hình HR82 Hossain Rodi [7] Mơ hình COHERENS cho vùng ven bờ CTSH tính tới thuật tốn khơ/ướt (drying and wetting algorithm) dựa mơ hình Burchard Bolding [2] Mơ hình áp dụng yếu tố khơ tới điều kiện phương trình động lượng độ sâu cột nước giảm tới giá trị nhỏ xác định Các giá trị tham số mặc định mơ hình phương trình xem chi tiết tài liệu mơ hình COHERENS [11] Độ nhám đáy mơ hình áp dụng cho điều kiện địa phương khu vực nghiên cứu 0.0035m Phương pháp thiết lập điều kiện biên mở biển mơ hình thảo luận nghiên cứu Hedstrom [5] Hirsch [6] Việc áp dụng điều chỉnh phương pháp tiến hành nghiên cứu Roed Cooper [15] với điều kiện cụ thể địa phương việc giải phương trình vận chuyển Pha biên độ thành phần sóng triều điểm biên gần bờ có sử dụng kết phân tích số điều hòa thủy triều từ số liệu thực đo Tại điểm biên phía ngồi sử dụng số liệu pha biên độ thủy triều từ sở liệu Fes2004 [12] Việc nội suy tiến hành mơ hình AG (Andersen-Gill) Andersen [1] Đối với điều kiện biên sông, sử dụng giá trị lưu lượng nước trung bình (mùa mưa mùa khơ) cửa sơng khu vực nghiên cứu Cửa Ông, Cấm-Bạch Đằng, Lạch Tray, Văn Úc, Thái Bình, Trà Lý, Ba Lạt, Ninh Cơ Đáy (bảng 1) Các số liệu xử lý từ số liệu đo nhiệm vụ hợp tác theo nghị định thư Việt Nam- Vương quốc Bỉ Độ mặn sông sử dụng số liệu độ mặn trung bình thời gian quan trắc Bảng Lưu lượng nước trung bình (m3/s) cửa sơng khu vực nghiên cứu Sơng Mơng Dương (Cửa Ơng) Bạch Đằng-Cấm Lạch Tray Văn Úc Thái Bình Trà Lý Ba Lạt Ninh Cơ Đáy Lưu lượng nước trung bình mùa Mùa khô 120 252 55 217 96 121 352 88 360 Mùa mưa 300 1505 525 1261 560 665 1927 483 1969 Mơ hình chạy với thời gian tháng cho kịch bản, điều kiện ban đầu mơ hình sử dụng kết sau ngày chạy Bước thời gian tính tốn mơ hình 20 giây Điều kiện khí tượng dùng cho mơ hình số liệu gió thực đo trạm Hịn Dáu thời gian tính tốn Hiệu chỉnh kiểm chứng kết mơ hình Mơ hình hiệu chỉnh kiểm chứng mùa mưa mùa khô thông qua việc so sánh số liệu đo đạc mực nước Hòn Dáu, số liệu đo dòng chảy số điểm khu vực nghiên cứu kết tính tốn từ mơ hình So sánh kết sau lần hiệu chỉnh cuối cho thấy phù hợp tương đối tính tốn số liệu đo đạc (hình hình 3) modeled Mơ hình observed Quan trắc 1.5 zeta (m) 0.5 0 200 400 600 800 1000 1200 1400 -0.5 -1 -1.5 24 +/- 14 days -2 Start date: 08/01/2010 time (h) (a) Mơ hình modeled observed Quan trắc 1.5 zeta (m) 0.5 0 200 400 600 800 1000 1200 1400 -0.5 -1 -1.5 -2 +/- 14 days (b) time (h) Start date: 03/01/2010 End date: 04/30/2010 Hình Dao động mực nước tính mơ hình quan trắc Hòn Dáu (a- mùa mưa, b mùa khơ) quan trắc tính tốn (a) vận tốc (cm/s) 40 12 14 16 18 20 22 10 12 -40 12 14 16 18 20 22 18 19 20 10 12 -40 -80 -120 -120 40 60 quan trắc tính tốn (c) 10 11 12 13 14 15 16 17 18 19 20 21 22 -20 quan trắc tính tốn (d) 40 vận tốc (cm/s) vận tốc (cm/s) (b) 40 -80 20 quan trắc tính tốn 80 vận tốc (cm/s) 80 20 10 11 12 13 14 15 16 17 21 22 -20 -40 -40 Hình Số liệu đo đạc dịng chảy tính tốn mơ hình khu vực cửa Ba Lạt (a- thành phần vận tốc theo hướng u, mùa mưa; b- thành phần vận tốc theo hướng v, mùa mưa; c- thành phần vận tốc theo hướng u, mùa khô; d- thành phần vận tốc theo hướng v, mùa khô) Để đánh giá ảnh hưởng gió đến phân bố độ mặn hồn lưu khu vực nghiên cứu, 10 kịch tính tốn liên quan thiết lập Trong đó, kịch khơng tính tới điều kiện gió (một cho mùa mưa cho mùa khơ) Các kịch cịn lại thiết lập với điều kiện vận tốc gió biến đổi khác nhau: 3.0m/s, 3.5m/s, 4m/s, 4.5m/s mùa mùa mùa khô (bảng 2) Bảng Các kịch thiết lập với điều kiện gió khác STT Kịch 10 Mùa khơ Điều kiện gió Vận tốc Hướng - Mùa mưa 3.0 3.5 4.0 4.5 3.0 3.5 4.0 4.5 Mùa khô Mùa mưa NE NE NE NE SE SE SE SE III KẾT QUẢ VÀ THẢO LUẬN Các kết tính tốn phân tích mơ hình cho thấy trường dòng chảy dư vùng ven bờ châu thổ sơng Hồng có khác biệt mùa mưa mùa khơ (hình 4) (a) (b) Hình Dịng dư tầng mặt ven bờ châu thổ sơng Hồng (a- mùa khô; b- mùa mưa) Vào mùa khô, tải lượng nước từ sơng đưa nhỏ, dịng dư có hướng chủ yếu dọc bờ từ đơng bắc xuống tây nam với vận tốc khoảng 0.2-0.4m/s vùng ven bờ giảm dần tới giá trị nhỏ phía ngồi (hình 4a) Trong đó, vào mùa mưa tải lượng nước sông đưa lớn nên vận tốc dòng dư lớn với giá trị khoảng 0.3-0.5m/s, hướng dòng dư vùng ven bờ chủ yếu đơng bắc xuống tây nam (hình 4b) Các kết tính tốn cho thấy phân bố độ mặn khu vực nghiên cứu có khác biệt đáng kể mùa mưa mùa khô Vào mùa mưa ảnh hưởng lưu lượng nước từ sông đưa lớn nên vùng nước ven bờ có độ mặn nhỏ 15‰ (hình 4b), vào mùa khơ, độ mặn với giá trị khoảng 25‰ tiến sát vào ven bờ, đặc biệt độ mặn cao vùng ven bờ phía tây nam (hình 4a) (b) (a) Hình Phân bố độ mặn vùng ven bờ CTSH mùa mưa (a- khơng có gió, b- trường gió với vận tốc trung bình 3.5m/s) (a) (b) Hình Phân bố độ mặn vùng ven bờ CTSH mùa khơ (a- khơng có gió, b- trường gió với vận tốc trung bình 4.0m/s) Phân tích ảnh hưởng trường gió thơng qua kịch có gió khơng có gió tác động cho thấy vào mùa mưa ảnh hưởng trường gió, khối nước sơng phát triển mạnh phía ngồi hơn, với dịng chảy tức thời hướng từ phía cửa sơng phía ngồi khơi nhiều định hướng dọc bờ (hình 5) Vào mùa khơ, ảnh hưởng trường gió làm tăng cường vận tốc dòng chảy dọc bờ (do với hướng gió) Trong độ mặn vùng ven bờ tăng lên xâm nhập khối nước biển vào vùng ven bờ ảnh hưởng gió bề mặt (hình 6) Những nghiên cứu động thái phát triển khối nước sông thực mơ hình tốn Chao Boicourt [3], Chao [4]; Kourafalou nnk [8] vùng có biên độ triều nhỏ Ruddick nnk [16] vùng ảnh hưởng thủy triều mạnh cửa sông Rhin Các kết nghiên cứu khối nước sông đưa trước hết mở rộng phía biển sau dịch chuyển phía bên phải (ở vùng bắc bán cầu) Trước khối nước sông quay trở lại vùng ven bờ, khối nước chệch hướng tạo thành vệt nước sông ven bờ Tuy nhiên nghiên cứu này, vệt nước cách rõ rệt vào mùa mưa di chuyển khối nước bên phải số trường hợp cân với dịng triều (hình 5) Hiện tượng thể rõ rệt kịch mùa khơ (hình 6) tải lượng nước từ sơng đưa nhỏ vùng nước hòa trộn dải ven bờ nhỏ so với mùa mưa (a) (b) Hình Profiler độ mặn mặt cắt phía ngồi cửa Ba Lạt mùa khơ (a- khơng có gió; b- gió với vận tốc 4m/s) (a) (b) Hình Biến đổi độ mặn mặt cắt phía ngồi cửa Ba Lạt mùa khô (a-kỳ triều kém; b- kỳ triều cường; nét đứt – thời điểm nước lớn; nét liền -thời điểm nước rịng) Phân tích profile độ mặn mặt cắt phía ngồi cửa Ba Lạt (hình 10) cho thấy vào mùa khơ, ảnh hưởng gió làm vùng nước có độ mặn lớn tiến sâu vào gần bờ (khoảng 23km), biến đổi độ mặn theo chiều thẳng đứng (phân tầng) tăng lên đáng kể (hình 7b) Trong khơng có ảnh hưởng gió độ mặn đồng theo phương thẳng đứng (hình 7a) Cũng mắt cắt phía ngồi cửa Ba Lạt, kỳ triều mùa khô độ mặn vùng nước xáo trộn mạnh phía ngồi cửa Ba Lạt có phạm vi khoảng 15km vào kỳ triều cường vùng nước xáo trộn giới hạn khoảng 2-10km (hình 8) Chênh lệch độ mặn thời điểm nước lớn nước ròng vào kỳ triều cường khu vực lớn vào kỳ triều (b) (a) Hình Profile độ mặn mặt cắt phía ngồi cửa Ba Lạt mùa mưa (a- khơng có gió; b- gió với vận tốc 4m/s) (a) (b) Hình 10 Biến đổi độ mặn mặt cắt cửa Ba Lạt mùa mưa (a-kỳ triều kém; b-kỳ triều cường; nét đứt – thời điểm nước lớn; nét liền -thời điểm nước ròng) Vào mùa mưa mặt cắt phía ngồi cửa Ba Lạt, ảnh hưởng gió làm cho khối nước sơng phát triển mạnh phía ngồi tầng mặt khối nước biển với độ mặn lớn lại tiến sâu vào phía cửa sơng tầng đáy (hình 9b) Tại mặt cắt vào mùa mưa tác động gió khối nước sơng tồn chủ yếu lớp nước mặt (hình 9a) với độ dày khoảng 5m (với phạm vi khoảng 25km từ bờ ra), ảnh hưởng gió khối nước sơng phát triển phía ngồi lớp bề mặt với phạm vi khoảng 40-45km (hình 9) Tại mặt cắt mùa mưa, độ mặn biến đổi lớn theo pha dao động mực nước triều: vào kỳ triều khu vực xáo trộn mạnh vùng nước cách bờ 10-20km, vào kỳ triều cường vùng nước xáo trộn mạnh có phạm vi khoảng 5-25km cách bờ (hình 10) Chênh lệch độ mặn nước lớn nước ròng nhỏ vào kỳ triều tương đối lớn (có thể lên tới 10-15‰) kỳ triều cường IV KẾT LUẬN Điều kiện thủy động lực vùng ven bờ châu thổ sông Hồng chịu tác động tổng hợp yếu tố trường gió, dao động mực nước khối nước từ hệ thống sơng Hồng-Thái Bình đưa Dao động mực nước triều làm thay đổi phạm vi vùng xáo trộn mạnh khối nước sông nước biển: cách bờ 5-25km vào mùa mưa nhỏ cách bờ 15km vào mùa khô Vào mùa mưa, phạm vi vùng xáo trộn vào kỳ triều cường mở rộng kỳ triều Ngược lại, vào mùa khô kỳ triều phạm vi vùng xáo trộn rộng vào kỳ triều Trường gió vào mùa khơ làm tăng cường vận tốc dịng chảy dư vận chuyển khối nước dọc bờ từ đông bắc xuống tây nam Trong vào mùa mưa ảnh hưởng gió làm tăng vận tốc dòng chảy tức thời hướng từ bờ phía ngồi khơng làm ảnh hưởng nhiều đến tính chất dịng dư vùng ven bờ dòng chảy dư vào mùa chịu tác động lớn từ khối nước sông Vào mùa khô, ảnh hưởng gió làm cho phạm vi phân bố khối nước sông thu hẹp, khối nước biển tiến sát vào gần bờ Gió mùa khơ làm tăng gradient theo chiều thẳng đứng độ mặn vùng ven bờ Trong vào mùa mưa, gió làm tăng cường phát tán khối nước sơng phía ngồi, tăng phạm vi ảnh hưởng khối nước sông lớp nước bề mặt ảnh hưởng nước biển tầng đáy V LỜI CẢM ƠN Bãi báo này có sử dụng tư liệu Nhiệm vụ Hợp tác theo nghị định thư Khoa học Công nghệ Việt Nam Vương quốc Bỉ giai đoạn 2010-2011: “Nghiên cứu ứng dụng hệ thống mơ hình thủy nhiệt động lực-sinh thái biển phục vụ nghiên cứu quản lý tài nguyên biển vùng ven bờ Việt Nam, tác giả xin chân thành cảm ơn hỗ trợ quý báu Nhóm tác giả xin bày tỏ lòng biết ơn sâu sắc tới TS Trần Đình Lân PGS TS Trần Đức Thạnh, có góp ý khích lệ quan trọng để báo hoàn thành TÀI LIỆU THAM KHẢO Andersen O., 1995 Global ocean tides from ers1 and to-pex/poseidon altimetry J Geophys Res 100(25): 249-259 Burchard H, Bolding K., 2002 GETM- A general estuarine transport model Institute for Environment and Sustainability, Joint Research Centre, Ispra, Italy Chao S, Boicourt W., 1986 Onset of estuarine plumes J Phys Oceanogr 16:2137-2149 Chao S., 1988 River-forced estuarine plumes J Phys Oceanogr 18:72-88 Hedstrom G., 1979 Nonreecting boundary conditions for nonlinear hyperbolic systems J Comput Phys 30: 222-237 Hirsch C., 1990 Numerical computation of internal and external ows, vol 2: Computational methods for inviscid and viscous ows New York: Wiley Hossain M, Rodi W., 2002 A turbulence model for buoyant ows and its application to vertical buoyant jets, vol Oxford: Pergamon Press Kourafalou V, Oey L, Wang J, Lee T., 1996 The fate of river discharge on the continental shelf modeling the river plume and the inner shelf coastal current J Geophys Res 101: 3415-3434 McDougall T, Jackett D DR and Wright, Feistel R., 2003 Accurate and computationally eficient algorithms for potential temperature and density of seawater J Atmos Oceanic Technol 20:730-741 10 10 Merri T Jone, Pauline W., Raymond N Cramer, 2009 User Guide to the centernary edition of the GEBCO digital atlas and its datasets Natural environment research council 11 Luyten P., 2011 COHERENS - A Coupled Hydrodynamical-Ecological Model for Regional and Shelf Seas: User Documentation Version V2.0 RBINS-MUMM, Royal Belgian Institute of Natural Sciences, Brussels, Belgium 12 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 13 Manh D, Yanagi T, 2000 A study on the residual flow in the gulf of tongkin J Oceanogr 56: 59-68 14 Phillips N., 1957 A coordinate system having some special advantages for numerical forecasting Journal Meteor 14: 184-185 15 Roed L, Cooper C., 1987 A study of various open boundary conditions for wind-forced barotropic numerical ocean models Amsterdam: Elsevier 16 Ruddick K, Deleersnijder E, Luyten P, Ozer J., 1995 Haline stratification in the rhinemeuse freshwater plume: A three-dimensional model sensitivity analysis Cont Shelf Res 15: 1597-1630 17 Trần Đức Thạnh, Lưu Văn Diệu, 2006 Những vấn đề môi trường bật dải ven bờ tây vịnh Bắc Bộ Khoa học công nghệ biển Phụ trương (T.6) Tr 3-14 18 UNESCO, 1981 Background papers and supporting data on the international equation of state 1980 UNESCO Technical Papers in Marine Science 37, 144 pp 19 Dinh Van Uu, Ha Thanh Huong, 2008 Model for water circulation in tidal dominated estuarine regions Int J Water Res Environ Eng 23: 33-39 20 Van Maren D, Hoekstra A, Hoitink A, 2004 Tidal flow asymmetry in the diurnal regime: bed-load transport and morphologic changes around the red river delta Ocean Dynam 54: 424-434 21 Vũ Duy Vĩnh, 2006 Áp dụng mơ hình tốn mơ lan truyền dầu xảy cố tràn dầu vịnh Bái Tử Long Khoa học Công nghệ Biển Phụ trương (Tháng 6/2006) 22 Vũ Duy Vĩnh, 2007 Mơ hình tốn nghiên cứu thuỷ động lực chất lượng nước khu vực Vịnh Bái Tử Long Tài nguyên Môi trường biển, tập XII, Nxb KH&KT, Hà Nội 2007 Tr.93-116 23 Vũ Duy Vĩnh, Đỗ Đình Chiến, Trần Anh Tú, 2008 Mơ hình tốn chiều nghiên cứu chất lượng nước vịnh Hạ Long Tài nguyên Môi trường biển, tập XIII Nxb KH&KT, Hà Nội Tr 318-327 24 Vũ Duy Vĩnh, 2011 Một số kết bước đầu ứng dụng mơ hình sinh thái vùng biển Cát Bà-Hạ Long Tài nguyên Môi trường biển, tập XVI Nxb Khoa học Tự nhiên Công nghệ Hà Nội Tr.215-229 25 Vũ Duy Vĩnh, Nguyễn Đức Cự, Trần Đức Thạnh, 2011 Ảnh hưởng đập Hịa Bình đến phân bố trầm tích lơ lửng vùng ven bờ châu thổ sông Hồng Kỷ yếu Hội nghị Khoa học Công nghệ biển lần thứ 5, tập Địa chất, Địa lý- Địa vật lý Tr 465-475 26 World Ocean Atlas 2009 National Oceanographic Data Center 30-03-2010 http://www.nodc.noaa.gov/OC5/WOA09/pr_woa09.html Retrieved 19-5-2010 11 THE INFLUENCE OF SURFACE WIND ON THE SALINITY DISTRIBUTION AND CIRCULATION IN THE COASTAL WATERS OF THE RED RIVER DELTA, VIETNAM VU DUY VINH, KATRIJN BAETENS, PATRICK LUYTEN, TRẦN ANH TÚ, NGUYỄN THỊ KIM ANH A sufficiently well validated hydro-physics model for the coastal zone of the Red River delta area in Vietnam was presented The development code was COHERENS V2.0, a threedimensional hydrodynamic multi-purpose model for coastal and shelf seas based on the finite element method The results of the model with 10 scenarios simulations show the role of wind condition combined with tidal oscillation and river discharge on salinity distribution and circulation of the red river coastal area Wind NE in the dry season make increasing velocity of residual current, intensify salinity in the coastal zone and gradient of salinity on the vertical On the other hand, wind SE in the rainy season makes intensify fresh water to seaward, increasing range of river water on the surface and seawater on the bottom in the coastal zone 12