Bài viết này ứng dụng mô hình ROMS để mô phỏng trường nhiệt mặt biển khu vực miền Trung Việt Nam thông qua sử dụng các sơ đồ khép kín rối khác nhau. Các sơ đồ khép kín rối Brunt-Vaisala và KPP được ứng dụng với cùng điều kiện mô phỏng nhằm đưa ra những so sánh về tính thích hợp và hiệu quả giữa hai phương pháp tại khu vực biển miền Trung Việt Nam.
Tạp chí Khoa học Cơng nghệ Biển; Tập 16, Số 4; 2016: 364-372 DOI: 10.15625/1859-3097/16/4/7212 http://www.vjs.ac.vn/index.php/jmst ĐÁNH GIÁ ẢNH HƯỞNG CỦA CÁC PHƯƠNG PHÁP KHÉP KÍN RỐI TRONG MƠ PHỎNG TRƯỜNG NHIỆT ĐỘ NƯỚC MẶT BIỂN KHU VỰC BIỂN MIỀN TRUNG VIỆT NAM Vũ Thị Vui*, Nguyễn Kim Cương, Nguyễn Minh Huấn Khoa Khí tượng Thủy văn Hải dương học, Trường Đại học Khoa học Tự nhiên, ĐHQGHN * E-mail: vuivt89@gmail.com Ngày nhận bài: 25-12-2015 TÓM TẮT: Bài báo ứng dụng mơ hình ROMS để mơ trường nhiệt mặt biển khu vực miền Trung Việt Nam thơng qua sử dụng sơ đồ khép kín rối khác Các sơ đồ khép kín rối Brunt-Vaisala KPP ứng dụng với điều kiện mô nhằm đưa so sánh tính thích hợp hiệu hai phương pháp khu vực biển miền Trung Việt Nam Các kết mô hai sơ đồ so sánh với tập đồ “điều kiện tự nhiên môi trường vùng biển Việt Nam kế cận” (năm 2009) so sánh với Những gợi ý việc sử dụng sơ đồ rút Kết nghiên cứu ứng dụng việc sử dụng lựa chọn sơ đồ khép kín rối để mơ trường nhiệt cho khu vực cụ thể Đây vấn đề tương đối cần thiết sử dụng mơ hình số Từ khóa: ROMS, KPP, Brunt-Vaisala, SST, sơ đồ khép kín rối MỞ ĐẦU Trong lĩnh vực nghiên cứu biển, mơ hình số trị phương pháp nghiên cứu ưu việt, vừa tiết kiệm vừa hiệu quả, cho phép mô trạng thái mơi trường biển sóng, dòng chảy, nhiệt-muối, … quy mơ khu vực biển tồn cầu Để khép kín hệ phương trình, mơ hình số trị 3D cần sử dụng phương pháp khép kín rối, việc đánh giá tính hiệu thích hợp phương pháp khép kín rối khu vực biển thử nghiệm quan trọng sử dụng mơ hình số trị chiều hải dương học Bài báo sử dụng hệ thống mô hình đại dương quy mơ khu vực ROMS (Regional Ocean Modeling System) mô trường nhiệt lớp nước mặt biển khu vực biển miền Trung Việt Nam sử dụng hai phương án khép kín rối: Brunt-Vaisala KPP, nhằm so sánh 364 đánh giá khác biệt, tính hiệu sử dụng phương án tính tốn mơ PHƯƠNG PHÁP TÍNH TỐN Giới thiệu tần số Brunt-Vaisala sơ đồ khép kín rối KPP Tần số Brunt-Vaisala Tần số Brunt-Vaisala gọi tần số ổn định tần số phân tầng Tần số đánh giá độ ổn định lớp nước biển, biến động lực học chất lỏng phân tầng Trong dạng đơn giản nhất, tần số coi tần số dịch chuyển theo phương thẳng đứng phần tử chất lỏng Nước biển đại dương thường phân bố tương đối ổn định theo phương thẳng đứng Tuy nhiên, thực tế tác động khác nhau, xảy tượng nước có mật độ thấp lại nằm nơi có độ sâu lớn Đánh giá ảnh hưởng phương pháp … lớp nước có mật độ cao, theo quy luật vật lý thể qua định luật độ Archimedes, xảy tượng dịch chuyển lên phía khối nước có mật độ thấp đồng thời nước có mật độ cao hơn chìm xuống Khi phân tầng mật độ nước ổn định, khối nước dịch chuyển có xu quay vị trí ban đầu Trong q trình dịch chuyển, qn tính, thể tích nước vượt q vị trí trở lại, từ làm xuất dao động quán tính Tần số dao động qn tính tần số Brunt-Vaisala Để xác định tần số BruntVaisala, sử dụng cơng thức: N2 - gE Trong đó: N- tần số Brunt-Vaisala; g- gia tốc trọng trường; E- độ ổn định nước biển [1] Sơ đồ khép kín rối KPP Sơ đồ KPP (K-Profile Parameterization: tham số hóa cấu trúc phân bố lượng rối K) xây dựng cho lớp biên bề mặt kết hợp với tham số đơn giản để xử lý trình xáo trộn bên đại dương Các sơ đồ KPP đánh giá cao ứng dụng vùng đại dương nước sâu Phương pháp sơ đồ KPP kết nối tham số rối phương thẳng đứng riêng biệt lớp biên mặt lớp nước bên Cơng thức tính dựa thuyết tương tự áp dụng cột nước độ sâu lớp biên tính tốn Tham số sau kết hợp bên để xác định ranh giới địa phương, sóng nội hiệu ứng xáo trộn khuếch tán kép (khuếch tán kép sinh từ tốc độ khuếch tán phân tử khác biệt nhiệt độ độ muối) [2] Giới thiệu hệ thống mơ hình ROMS Hệ thống mơ hình đại dương khu vực (ROMS) hệ thống mơ hình số trị có hệ tọa độ thích ứng địa hình, mặt thống, ba chiều, với mục đích tính tốn hồn lưu, sinh thái chu trình sinh-địa-hóa học khu vực biển ven khác Hệ thống mơ hình ROMS phát triển thành hệ thống mơ hình biển thực đa mục đích, có khả thể nhiều quy mơ khơng gian/thời gian nhiều vấn đề đa dạng [3] Bài báo sử dụng phiên ROMS tổ chức IRD ROMS_AGRIF, hỗ trợ công cụ ROMSTOOLS [4], tập trung khai thác ứng dụng tính tốn trường nhiệt cho khu vực biển miền Trung Việt Nam với hai phương pháp khép kín rối: theo tần số Brunt-Vaisala theo sơ đồ khép kín rối KPP Các phương trình mơ hình ROMS viết tọa độ Đề theo phương ngang tọa độ Sigma theo phương thẳng đứng, độ sâu phi thứ nguyên tính theo cơng thức: z h Trong đó: z tọa độ theo phương thẳng đứng hệ tọa độ Đề các; ζ- mực mặt thoáng; h- độ sâu đáy biển Với cách chuyển đổi -1≤ σ ≤ 0, σ = tương ứng với mặt thoáng z = ζ, σ = -1 tương ứng với đáy biển z = -h ROMS giải hệ phương trình Navier-Stokes trung bình Reynolds sử dụng xấp xỉ Boussinesq xấp xỉ thủy tĩnh Hệ phương trình thủy động lực chủ đạo - dạng thơng lượng, hệ tọa độ ngang Đề sigma theo độ sâu, viết sau: HZ u uH Z u vHZ u H Z u H p u fH Z v Z HZ g (u ' w ' ) t x y 0 x x H Z HZ v t uH Z v x vH Z v y H Z v 0 fH Z u H Z p y HZ g y (v ' w ' v ) H Z p g Hz 0 0 Với phương trình liên tục: 365 Vũ Thị Vui, Nguyễn Kim Cương, … t H Zu x H Z v y H Z 0 phương trình vận chuyển đại lượng vơ hướng: HZC t uH Z C x vH Z C y H Z C Trong đó: u, v Ω thành phần vận tốc dòng chảy theo phương ngang (x y) theo phương thẳng đứng (theo hệ tọa độ sigma, σ) tương ứng; ζ độ cao nước dâng mặt thống trung bình; Hz hệ số tỉ lệ theo chiều thẳng đứng f tham số Coriolis Dấu gạch ngang biến thể giá trị trung bình theo thời gian dấu nháy (’) thể nhiễu động rối Áp suất p; ρ ρ0 mật độ tổng cộng mật độ chuẩn; g gia tốc trọng trường; υ hệ số nhớt phân tử C yếu tố vật lý yếu tố mơi trường (ví dụ nhiệt độ, độ muối, hàm lượng trầm tích lơ lửng, chất vẩn, thực vật nổi, động vật nổi, nitơ vơ hòa tan…); Csource nguồn sinh/mất thành phần yếu tố C Cuối cùng, hàm: f (T , S , p ) (C ' w ' C ) C source H Z Với T- nhiệt độ, S- độ muối, p- áp suất, phương trình trạng thái nước biển Hệ phương trình khép kín theo nhiều phương án Trong nghiên cứu này, để so sánh sử dụng phương án: theo tần số Brunt-Vaisala sơ đồ khép kín rối KPP Khi sử dụng tần số Brunt-Vaisala, hệ số xáo trộn tính theo giá trị tương ứng tần số Brunt-Vaisala Còn phương án sử dụng sơ đồ khép kín rối KPP, với giả thiết trình xáo trộn diễn cột nước, hệ số xáo trộn tính theo tham số K theo phương thẳng đứng hình vẽ 1c Khơng gian tính rời rạc hóa theo phương ngang lưới tính Arakawa-C theo phương thẳng đứng lưới sigma Mơ hình sử dụng phương pháp sai phân trung tâm Hệ số xáo trộn Lớp xáo trộn mặt Lớp xáo trộn Lớp xáo trộn đáy a) b) c) Hình a) Lưới Arakawa-C theo phương ngang; b) Lưới theo phương thẳng đứng mơ hình ROMS [3]; c) Các hệ số xáo trộn sử dụng lớp xáo trộn [5] Thiết lập phương án tính tốn Các bước chuẩn bị cho việc thiết lập mơ hình Mơ hình thiết lập với điều kiện 366 dựa bước chuẩn bị số liệu đầu vào nhờ công cụ ROOMSTOOLS [5] Các file số liệu sử dụng làm đầu vào đầu mơ hình ROMS lưu theo dạng netCDF (*.nc) Đánh giá ảnh hưởng phương pháp … Miền tính cho khu vực biển miền Trung Việt Nam nghiên cứu thể hình vẽ có đặc điểm sau: Vĩ độ: từ 90N đến 190N; Kinh độ: từ 1050E đến 1180E Lưới tính có độ phân giải khơng gian theo phương ngang 0,250 chia thành 10 tầng sigma theo độ sâu Kết tính tốn từ mơ hình trường giá trị trung bình nhiều năm: độ cao mặt mực, nhiệt độ, độ muối, vận tốc dòng chảy theo hướng, … Bài báo quan tâm đưa đánh giá kết mô trường nhiệt độ nước biển theo phương ngang tầng nước theo phương thẳng đứng (từ mặt biển xuống độ sâu 700 m) mặt cắt từ (140N, 1090E30’) đến (140N, 1180E), thể hình ảnh kết Địa hình khu vực nghiên cứu: Thiết lập từ nguồn số liệu ETOTO2 - số liệu địa hình tồn cầu với độ phân giải phút Việc thiết lập trường tác động cho mơ hình: Sử dụng số liệu tồn cầu trung bình tháng thơng lượng khí tượng bề mặt biển từ số liệu COADS05 với độ phân giải 0,50 Để thiết lập điều kiện biên cho mơ hình, giá trị trung bình tháng thời gian dài yếu tố hải văn cho tồn khu vực tính với độ phân giải độ kinh vĩ theo sở liệu WOA2009, riêng điều kiện ban đầu trường dòng chảy giá trị cho tồn khu vực QĐ Hoàng Sa Mặt cắt quan tâm QĐ Trường Sa Hình Địa hình khu vực nghiên cứu mặt cắt quan tâm từ (140N, 1090E30’) đến (140N, 1180E) Các thơng số thiết lập cho mơ hình theo hai sơ đồ khép kín rối Thơng qua việc khai báo tham số phù hợp file cppdefs.h mơ hình, mơ hình thực thi với hai phương pháp khép kín rối: Theo tần số Brunt-Vaisala theo sơ đồ khép kín rối KPP Hệ số ma sát đáy: 0,0003 m/s [5]; Sử dụng điều kiện trượt tự biên cứng Bước thời gian: 1.200 s; Bước thời gian bước thời gian tối ưu cho phép chạy mơ hình theo phương pháp tần số Brunt-Vaisala Phương pháp khép kín rối theo tần số Brunt-Vaisala đòi hỏi bước thời gian nhỏ so với sơ đồ khép kín rối KPP, để so sánh kết từ phương pháp, tác giả sử dụng bước thời gian hai phương án Mật độ nước: 1.025 kg/m3 [5]; KẾT QUẢ NGHIÊN CỨU Mơ hình thiết lập tính tốn thời gian năm từ thời điểm ban đầu với tham số thiết lập sau: 367 Vũ Thị Vui, Nguyễn Kim Cương, … QĐ Hoàng Sa QĐ Hoàng Sa QĐ Trường Sa QĐ Trường Sa a) Bản đồ phân bố trường nhiệt độ trung bình tầng mặt mùa đông [6] d) Bản đồ phân bố trường nhiệt độ trung bình tầng mặt mùa hè [6] 19 19 18 18 17 17 QĐ Hoàng Sa QĐ Hoàng Sa 31 16 16 30 15 29 28 14 15 14 27 13 13 26 12 25 24 11 11 23 QĐ Trường Sa 10 12 QĐ Trường Sa 10 9 106 108 110 112 114 116 118 106 b) Sơ đồ phân bố trường nhiệt độ trung bình nước biển tầng mặt mùa đông theo phương pháp khép kín rối Brunt-Vaisala 108 110 112 114 116 118 e) Sơ đồ phân bố trường nhiệt độ trung bình nước biển tầng mặt mùa hè theo phương pháp khép kín rối Brunt-Vaisala 19 19 18 18 17 17 31 QĐ Hoàng Sa 16 QĐ Hoàng Sa 16 30 15 29 28 14 15 14 27 13 13 26 12 25 24 11 QĐ Trường Sa 10 12 11 QĐ Trường Sa 23 10 9 106 108 110 112 114 116 118 c) Sơ đồ phân bố trường nhiệt độ trung bình nước biển tầng mặt mùa đơng theo phương pháp khép kín rối KPP 106 108 110 112 114 116 Hình Trường nhiệt độ nước mặt biển trung bình (0C) khu vực miền Trung Việt Nam 368 118 f) Sơ đồ phân bố trường nhiệt độ trung bình nước biển tầng mặt mùa hè theo phương pháp khép kín rối KPP Đánh giá ảnh hưởng phương pháp … a) Bản đồ phân bố trường nhiệt độ nước biển trung bình mặt cắt quan tâm mùa đông [6] d) Bản đồ phân bố trường nhiệt độ nước biển trung bình mặt cắt quan tâm mùa hè [6] b) Sơ đồ phân bố trường nhiệt độ nước biển trung bình mặt cắt quan tâm mùa đơng theo phương pháp khép kín rối Brunt-Vaisala e) Sơ đồ phân bố trường nhiệt độ nước biển trung bình mặt cắt quan tâm mùa hè theo phương pháp khép kín rối Brunt-Vaisala c) Sơ đồ phân bố trường nhiệt độ nước biển trung bình mặt cắt quan tâm mùa đông theo phương pháp khép kín rối KPP f) Sơ đồ phân bố trường nhiệt độ nước biển trung bình mặt cắt quan tâm mùa hè theo phương pháp khép kín rối KPP Hình Trường nhiệt độ nước biển trung bình (0C) mặt cắt quan tâm 369 Vũ Thị Vui, Nguyễn Kim Cương, … Hình trình bày kết mơ hai mơ hình khép kín rối phân bố số liệu nhiệt độ từ ATLAS “Điều kiện tự nhiên môi trường vùng biển Việt Nam kế cận” [6] Vào mùa đông (tháng - hình 3a, 3b, 3c), nhiệt độ khu vực biển miền Trung Việt Nam dao động khoảng từ 220C đến 280C, với nhiệt độ thấp khu vực cửa vịnh Bắc Bộ, nhiệt độ cao phía đơng nam khu vực nghiên cứu Trường nhiệt thể rõ xâm nhập lưỡi nước lạnh vào mùa đơng ảnh hưởng gió mùa Đơng Bắc Các kết mô với hai phương pháp khép kín rối cho kết tương đối tương đồng với trường nhiệt từ số liệu thực đo Kết mơ theo phương pháp khép kín rối KPP có xu hướng thể vùng nhiệt độ thấp ven bờ Ninh ThuậnBình Thuận rõ rệt so với phương pháp khép kín rối theo tần số Brunt-Vaisala nhiệt độ diện tích khu vực nước trồi phù hợp hơn, xuất rõ rệt khu vực nước từ bờ khơi từ tâm trồi (tại vĩ độ 140N), kéo theo nhiệt độ nước giảm Hình thể kết tính mặt cắt quan tâm - khu vực nghiên cứu, từ (140N, 1090E 30’) đến (140N, 1180E) Các kết cho thấy, việc sử dụng phương pháp khép kín rối (hình 4b, 4c, 4e, 4f) đưa trường nhiệt độ mặt cắt quan tâm tương đồng với giá trị thực đo (hình 4a, 4d), tương đồng so sánh với Vào mùa đông, trường nhiệt độ dao động khoảng từ 9-260C từ độ sâu 500 m lên mặt Còn vào mùa hè, từ mặt xuống đến độ sâu 500 m, nhiệt độ mặt cắt giảm dần từ 290C đến 90C Tuy nhiên, tầng nước sát mặt, kết tính theo phương pháp khép kín rối KPP thể độ dày lớp tựa đồng có xu hướng lớn phương pháp khép kín rối tần số BruntVaisala, vào mùa đơng, kết tính gần với giá trị thực đo hơn, lớp tựa đồng có độ dày khoảng 50 m Còn mùa hè, độ dày lớp tựa đồng theo thực đo khoảng 30 m, kết từ phương pháp tần số Brunt-Vaisala lại gần với giá trị thực đo Tại tầng nước sâu hơn, khác biệt phương pháp không đáng kể Vào mùa hè (tháng - hình 3d, 3e, 3f), nhiệt độ nước mặt biển khu vực miền Trung Việt Nam dao động khoảng từ 240C đến 300C Nhiệt độ thấp khu vực tâm nước trồi ven bờ Ninh Thuận-Bình Thuận nhiệt độ cao đạt 300C khu vực khơi khu vực nghiên cứu Kết từ phương pháp khép kín rối theo tần số BruntVaisala cho thấy tương đồng với thực đo hơn, 19 19 18 18 17 17 QĐ Hoàng Sa 16 QĐ Hoàng Sa 16 1.6 1.4 1.2 15 15 14 14 13 13 12 12 11 11 0.8 0.6 0.4 0.2 -0.2 -0.4 -0.6 -0.8 -1 QĐ Trường Sa 10 -1.2 -1.4 QĐTrường Sa 10 -1.6 106 108 110 112 114 116 118 a) 106 108 110 112 114 116 118 b) Hình Độ lệch nhiệt độ nước mặt biển trung bình (0C) phương pháp khép kín rối Brunt-Vaisala KPP khu vực miền Trung Việt Nam mùa đông (a) mùa hè (b) (hình ảnh kết nhiệt độ theo phương pháp tần số Brunt-Vaisala trừ kết nhiệt độ theo phương pháp KPP) Để so sánh kết tính tốn sử dụng phương pháp khép kín rối, chênh lệch 370 giá trị nhiệt độ mặt biển (SST) xử lý hình Hai phương pháp cho kết Đánh giá ảnh hưởng phương pháp … tương đối tương đồng thể xu phân bố trường nhiệt số liệu thực đo Vào mùa đơng, khu vực có trường nhiệt độ thấp phía bắc khu vực nghiên cứu, phương pháp khép kín rối KPP cho giá trị trường nhiệt độ nước bề mặt biển lớn so với giá trị trường nhiệt mô phương pháp khép kín rối KPP (hình 5a) Xu hồn tồn ngược lại vào mùa hè (hình 5b) Các kết mơ trường nhiệt sử dụng sơ đồ khép kín rối BruntVaisala cho kết nhiệt độ thấp vùng lạnh khu vực lưỡi lạnh vào mùa đông khu vực nước trồi vào mùa hè Đây khu vực có vận tốc thẳng đứng đối lưu nước trồi lớn Các khu vực ngồi khơi xa, nơi có sai khác nhiệt độ không lớn, hai phương pháp cho kết mơ có khác biệt khơng đáng kể cho thấy, mô cho khu vực cụ thể, cần lựa chọn ứng dụng sơ đồ khép kín rối phù hợp sơ đồ ảnh hưởng đáng kể tới kết mô Khi mô trường nhiệt độ nước biển tầng mặt khu vực Nam Trung Bộ, nên sử dụng phương pháp tần số Brunt-Vaisala Việc lựa chọn sơ đồ khép kín rối đòi hỏi kinh nghiệm hiểu biết định khu vực nghiên cứu Tuy nhiên, với khu vực ngồi khơi, nơi có phân bố nhiệt độ tương đối đồng ứng dụng hai phương pháp khép kín rối KẾT LUẬN Đinh Văn Ưu, Nguyễn Minh Huấn, 2003 Vật lý biển Nxb ĐHQG Hà Nội Tr 29-31 Website: https://www.myroms.org/wiki/Vertical_Mixing_Parameterizations Haidvogel, D B., Arango, H., Budgell, W P., Cornuelle, B D., Curchitser, E., Di Lorenzo, E., Fennel, K., Geyer, W R., Hermann, A J., Lanerolle, L., Levin, J McWilliams, J C., Miller, A J., Moore, A M., Powell, T M., Shchepetkin, A F., Sherwood, C R., Signell, R P., Warner, J C., and Wilkin, J., 2008 Ocean forecasting in terrain-following coordinates: Formulation and skill assessment of the Regional Ocean Modeling System Journal of Computational Physics, 227(7): 3595-3624 Penven, P., Cambon, G., Tan, T., Marchesiello, P., and Debreu, L., 2010 ROMS AGRIF/ROMSTOOLS user’s guide Institut de Recherche pour le Developpement, Dunkerque Gildas Cambon, Elodie Gutknecht, Marine Herrmann, 2012 The document of Ocean Modeling ROMS AGRIF/ROMS TOOLS, Summer school in Hanoi, Vietnam 2012, Vietnam ATLAS “Điều kiện tự nhiên môi trường vùng biển Việt Nam kế cận”, 2009 Nxb Khoa học tự nhiên Công nghệ Tr 68, 72, 75 Bài báo trình bày mô phân bố trường nhiệt độ lớp nước mặt biển khu vực biển miền Trung Việt Nam vào mùa đơng mùa hè mơ hình ROMS ứng dụng hai phương pháp khép kín rối khác nhau: Brunt-Vaisala KPP Các kết trường nhiệt mô tương đối tương đồng với với đồ phân bố chế độ trường nhiệt mặt biển, khác biệt hai phương pháp khơng lớn khu vực ngồi khơi Việc mô trường nhiệt lớp nước mặt biển sử dụng tần số Brunt-Vaisala cho kết nhiệt độ thấp trường nhiệt mô sơ đồ KPP vùng lạnh khu vực lưỡi lạnh phía bắc khu vực tính vào mùa đơng khu vực nước trồi vào mùa hè Điều chứng tỏ vận tốc thẳng đứng tăng lên phương pháp tần số Brunt-Vaisala cho giá trị nhiệt độ nước bề mặt biển nhỏ phù hợp với số liệu ATLAS so với phương pháp KPP Theo phương thẳng đứng, đánh giá định tính cho thấy độ dày lớp tựa đồng tính theo phương pháp tần số Brunt-Vaisala nhỏ tính theo phương pháp KPP, dẫn đến việc mô theo phương pháp KPP tin cậy mùa đông ngược lại vào mùa hè Tuy nhiên, đánh giá định lượng tầng nước sâu, số liệu nhiệt độ nước biển thực đo có chưa cho phép đánh giá cách tin cậy kết tính theo hai phương pháp Vấn đề cần tiếp tục nghiên cứu tương lai Các kết có Lời cảm ơn: Nghiên cứu tài trợ Trường Đại học Khoa học tự nhiên đề tài mã số TN.16.22 Các tác giả xin cảm ơn tài trợ TÀI LIỆU THAM KHẢO 371 Vũ Thị Vui, Nguyễn Kim Cương, … ASSESSMENT OF THE EFFECT OF DIFFERENT TURBULENT CLOSURE SCHEMES ON SIMULATING THE SEA SURFACE TEMPERATURE IN THE CENTRAL REGION OF VIETNAM’S EAST SEA Vu Thi Vui, Nguyen Kim Cuong, Nguyen Minh Huan Faculty of Hydro-Meteorology and Oceanography, Hanoi University of Science-VNU ABSTRACT: In this paper, the ROMS model was used to simulate the sea surface temperature (SST) in the central region of Vietnam’s East Sea Two turbulent closure schemes: Brunt-Vaisala and KPP were applied with the same conditions in order to compare the differences and the ability for using them The simulated SST fields were compared and in good agreement with the mechanism maps of SST The suggestions for using these schemes were analyzed This result raises a question and answers how to choose a better closure scheme to apply in simulating SST The solved problem is very necessary for applying numerical models to simulate SST field Keywords: ROMS, KPP, Brunt-Vaisala, closure scheme 372 ... Hình Trường nhiệt độ nước mặt biển trung bình (0C) khu vực miền Trung Việt Nam 368 118 f) Sơ đồ phân bố trường nhiệt độ trung bình nước biển tầng mặt mùa hè theo phương pháp khép kín rối KPP Đánh. .. đơng, khu vực có trường nhiệt độ thấp phía bắc khu vực nghiên cứu, phương pháp khép kín rối KPP cho giá trị trường nhiệt độ nước bề mặt biển lớn so với giá trị trường nhiệt mơ phương pháp khép kín. .. 3f), nhiệt độ nước mặt biển khu vực miền Trung Việt Nam dao động khoảng từ 240C đến 300C Nhiệt độ thấp khu vực tâm nước trồi ven bờ Ninh Thuận-Bình Thuận nhiệt độ cao đạt 300C khu vực khơi khu vực