Tổng quan luận án
Luận án tiến sĩ chuyên ngành Hải dương học (mã số: 62440228) với đề tài "Nghiên cứu mô phỏng chế độ thủy văn và hoàn lưu Vịnh Bắc Bộ" do nghiên cứu sinh Hà Thanh Hương thực hiện tại Trường Đại học Khoa học Tự nhiên – Đại học Quốc gia Hà Nội, dưới sự hướng dẫn khoa học của GS. Đinh Văn Mạnh, bảo vệ năm 2017.
Tính cấp thiết và khoảng trống nghiên cứu
Vịnh Bắc Bộ là vùng biển nửa khép kín ở Tây Bắc Biển Đông, có vị trí địa lý trong phạm vi tọa độ từ 17°00'N đến 21°30'N và từ 105°40'E đến 110°00'E, tổng diện tích khoảng 126.250 km². Vịnh có địa hình lòng chảo nông tương đối phức tạp với hơn 3.000 đảo lớn nhỏ, độ sâu trung bình chưa tới 60 m và không vượt quá 100 m. Vịnh thông với Biển Đông qua hai cửa: eo biển Quỳnh Châu ở phía Đông Bắc (rộng 35,2 km) và cửa chính phía Nam mở ra Biển Đông (nối từ đảo Cồn Cỏ, tỉnh Quảng Trị tới mũi Oanh Ca, đảo Hải Nam, Trung Quốc với bề rộng khoảng 200 km).
Về mặt khí tượng và thủy động lực, Vịnh Bắc Bộ chịu sự chi phối mạnh mẽ của chế độ gió mùa nhiệt đới (gió mùa Đông Bắc lạnh và khô/ẩm vào mùa đông; gió mùa Tây Nam nóng ẩm vào mùa hè), chế độ nhật triều và nhật triều không đều với biên độ lớn, cùng nguồn nước ngọt bổ sung theo mùa từ các hệ thống sông chính (sông Hồng, sông Thái Bình, sông Mã, sông Cả). Sự tương tác phức hợp giữa biển, khí quyển và lục địa tạo nên tính chất biến động mạnh mẽ của trường nhiệt độ, độ muối và hệ thống dòng chảy.
Mặc dù các đợt khảo sát đơn lẻ và các chương trình hợp tác quốc tế (Việt – Trung từ thập niên 1960, Việt – Xô giai đoạn 1960–1995) cùng các đề tài nghiên cứu biển trong nước đã thu thập được khối lượng số liệu thủy văn đáng kể, nhưng tại thời điểm nghiên cứu vẫn chưa có công trình chuyên khảo nào công bố một cách đầy đủ và đồng bộ về chế độ thủy văn và cấu trúc thủy động lực 3 chiều (3D) của toàn vịnh. Các nghiên cứu trước đây chủ yếu dừng lại ở các báo cáo chuyên đề về thủy triều, dòng triều hoặc phân bố nhiệt độ bề mặt; các cơ sở dữ liệu toàn cầu như World Ocean Atlas (WOA) chỉ cung cấp độ phân giải 1°, chưa phản ánh chi tiết các biến động quy mô nhỏ hơn 1/12°. Đặc biệt, cấu trúc hoàn lưu mùa hè trong vịnh vẫn tồn tại nhiều mâu thuẫn giữa các quan điểm (tồn tại xoáy thuận hay xoáy nghịch). Do đó, việc xây dựng một quy trình kết hợp phân tích xử lý số liệu đo đạc thực tế với mô hình số trị 3 chiều để tái hiện quy luật biến động của trường nhiệt – muối và dòng chảy là yêu cầu cấp thiết.
Mục tiêu và nội dung nghiên cứu
Mục tiêu nghiên cứu tổng quát:
Nghiên cứu phát triển và ứng dụng quy trình kết hợp mô đun (modul) phân tích số liệu và mô hình thủy nhiệt động lực 3 chiều (3D) GHER, nhằm mô phỏng quy luật hình thành và biến động các cấu trúc thủy văn (trường nhiệt độ, độ muối) và hoàn lưu nước trong Vịnh Bắc Bộ; xác lập các đặc trưng chế độ phục vụ công tác giám sát và dự báo môi trường biển.
Các nội dung nghiên cứu cụ thể:
- Nghiên cứu và hoàn thiện kỹ thuật phân tích số liệu vật lý hải dương Vịnh Bắc Bộ nhằm xác định đặc trưng chế độ các trường nhiệt độ, độ muối, dòng chảy theo quy mô mùa.
- Phát triển và hoàn thiện mô hình 3 chiều thủy nhiệt động lực đáp ứng yêu cầu mô phỏng các đặc trưng chế độ thủy văn có tính đến vai trò của các tác động quy mô vừa (triều, sóng, lưu lượng sông).
- Phát triển và ứng dụng quy trình kết hợp kỹ thuật phân tích và mô phỏng các đặc trưng chế độ thủy văn và dòng chảy Vịnh Bắc Bộ, làm rõ cấu trúc phân bố không gian – thời gian và hoàn thiện đầu vào cho hệ thống mô hình dự báo hải dương học nghiệp vụ.
Đối tượng và phạm vi nghiên cứu
- Đối tượng nghiên cứu: Cấu trúc 3 chiều của các trường nhiệt độ nước biển, độ muối và trường vận tốc dòng chảy (hoàn lưu) trong Vịnh Bắc Bộ.
- Phạm vi không gian: Toàn bộ thủy vực Vịnh Bắc Bộ và vùng biển kế cận (trích xuất tính toán trong phạm vi tọa độ từ 16°00'N đến 22°00'N và từ 105°00'E đến 110°00'E).
- Phạm vi thời gian: Quy mô biến động trung bình tháng nhiều năm và biến động theo mùa (mùa đông, mùa hè, các giai đoạn chuyển tiếp); khai thác chuỗi số liệu khảo sát lịch sử và cập nhật từ năm 1961 đến 2013.
Tổng quan tài liệu và vị trí của luận án
Luận án đã tổng hợp hệ thống các nghiên cứu về trường thủy văn và hoàn lưu ở Biển Đông và Vịnh Bắc Bộ qua các giai đoạn phát triển:
Các nghiên cứu về Biển Đông
Các công bố đầu tiên về hoàn lưu mặt Biển Đông xuất hiện trong Atlas của Hải quân Hoa Kỳ (1945), công trình của Dale (1956) dựa trên sự trôi dạt của tàu thuyền, và công trình tổng hợp của Wyrtki K. (1961) thiết lập bức tranh biến động mùa của dòng chảy mặt quy mô toàn biển. Xu và cộng sự (1982) phân tích chuỗi dữ liệu 1921–1970, chỉ ra cấu trúc xoáy thuận mùa đông (xoáy Aw gần đảo Luzon, Bw ở Đông Nam bán đảo Đông Dương) và xoáy nghịch mùa hè (xoáy Bs ở Đông Nam bán đảo Đông Dương, xoáy thuận Cs ở phía Đông Việt Nam). Các kết quả này được tái khẳng định qua nghiên cứu trường dòng địa chuyển của Yu và Liu (1993), Rong (1994), Zhou và cộng sự (1995).
Từ thập niên 1980–1990, các nghiên cứu ứng dụng mô hình số trị phát triển mạnh:
- Mô hình 2 chiều: Wang (1985), Zeng và cộng sự (1989, 1992), Su và Liu (1992), Li và cộng sự (1992a, 1992b, 1994), Liu và Su (1993) mô phỏng dòng chảy trung bình tầng mặt và trung bình độ sâu.
- Mô hình 3 chiều: Pohlmann (1987) và Mao và cộng sự (1992) ứng dụng mô hình nghiêng áp của Backhaus (1985); Zhang và cộng sự (1994, 1995), Shaw và Chao (1994), Chao và cộng sự (1996), Ly và Luong (1997, 1999), Chu và cộng sự (1999a, 1999b).
- Ứng dụng ảnh vệ tinh và đồng hóa dữ liệu: Chen (1983), Soong và cộng sự (1995), Shaw và cộng sự (1999), Wu và cộng sự (1999), Wang X. và cộng sự (2000), Morimoto và cộng sự (2000), Qi Quan và cộng sự (2016) sử dụng dữ liệu độ cao vệ tinh TOPEX/Poseidon (T/P) và Geosat.
- Tại Việt Nam: Báo cáo đề tài 48B.01 do Võ Văn Lành chủ nhiệm (2001) xây dựng bản đồ dòng chảy tầng mặt, 50 m, 100 m trên lưới 1°; Đinh Văn Ưu (1983, 1997) nghiên cứu hàm phân bố vạn năng và biến đổi mùa; đề tài KHCN 06-02 (2000) do Đinh Văn Ưu chủ trì ứng dụng mô hình 3D GHER với lưới phân giải 0,25°; đề tài KT03.03 (1991–1995) của Nguyễn Ngọc Thụy về thủy triều và năng lượng triều; đề tài KC09.24 (2005) xuất bản tập bản đồ hoàn lưu Biển Đông tầng mặt, 50 m và 150 m.
Các nghiên cứu tại khu vực Vịnh Bắc Bộ
Các số liệu điều tra tổng hợp Việt – Trung (1960–1964) và Việt – Xô (đến 1995) ban đầu ghi nhận dòng chảy tầng mặt trong vịnh chủ yếu chảy theo hướng Nam vào mùa đông và bảo toàn hướng trong mùa hè. Nghiên cứu thủy triều của Cao và Fang (1990), Liu và Zhang (1997), Hu và cộng sự (2003), Nguyễn Nguyệt Minh (2014), Phạm Văn Huấn (2010 – đề tài QG08.11) xác định sóng nhật triều $O_1$ và $K_1$ giữ vai trò thống trị, bán nhật triều $M_2$ là thứ yếu; dòng dư do triều đạt độ lớn khoảng 5 cm/s (Li, 1987; Cao và Fang, 1990).
Về hoàn lưu gió mùa trong vịnh, các nghiên cứu của Shaw và Chao (1994), Wu và cộng sự (1998, 2008), Xue và cộng sự (2004), Gan và cộng sự (2006), Chern và cộng sự (2010), Gao và cộng sự (2013), Yang Ding (2013), Trần Văn Chung (2015) đều thống nhất vào mùa đông trường gió Đông Bắc tạo dòng chảy hướng Nam ở bờ Tây và dòng hướng Bắc ở bờ Đông, cấu thành một xoáy thuận lớn trên toàn vịnh.
Tuy nhiên, cấu trúc hoàn lưu mùa hè tồn tại hai nhóm quan điểm trái ngược:
- Quan điểm tồn tại xoáy nghịch mùa hè: Sun và cộng sự (2001), Đinh Văn Mạnh và Yanagi (2000) xác định dòng dư tầng trên hướng Đông/Đông Bắc; Wu và cộng sự (2008) cho rằng xoáy nghịch hình thành do dòng chảy mạnh từ eo Quỳnh Châu; Gao và cộng sự (2013) dùng mô hình 3D POM chỉ ra xoáy nghịch ở phía Nam vịnh; Đinh Văn Ưu và cộng sự (2004 – đề tài KC09.17) mô phỏng lưới 1/12° ghi nhận xoáy nghịch.
- Quan điểm tồn tại xoáy thuận mùa hè: Xu và cộng sự (1980), Zhuang và cộng sự (1981), Bao và cộng sự (2005), Xia và cộng sự (1997) cho rằng vịnh duy trì cấu trúc xoáy thuận quanh năm do dòng mật độ; Zu (2005) và Yang Ding (2013) sử dụng mô hình FVCOM kết hợp quan trắc 1988–1989 chỉ ra sự hiện diện của xoáy thuận từ vĩ tuyến 19°N lên phía Bắc vịnh.
Khoảng trống luận án giải quyết
Luận án tập trung giải quyết sự thiếu nhất quán trong bức tranh hoàn lưu mùa hè và cấu trúc phân tầng nhiệt – muối 3D bằng cách thiết lập một công cụ xử lý số liệu chuyên biệt (LOWESS) để xây dựng trường nhiệt – muối ban đầu có độ tin cậy cao, kết hợp vận hành mô hình động lực học 3D GHER theo chu trình tương tác thuận nghịch có kiểm soát sai số.
Cơ sở lý thuyết và phương pháp nghiên cứu
Cơ sở dữ liệu sử dụng
Nghiên cứu tích hợp các nguồn số liệu lịch sử và quan trắc hiện đại bao gồm:
- Số liệu nhiệt độ bề mặt biển (SST): Trích xuất từ cơ sở dữ liệu MODAS (Modular Ocean Data Assimilation System) của Phòng Nghiên cứu Hải quân Hoa Kỳ, độ phân giải không gian 1/8° x 1/8°, tính trung bình tháng nhiều năm trong giai đoạn 1993–2013.
- Số liệu Atlas đại dương thế giới: Dữ liệu nhiệt độ, độ muối từ WOA (2001, 2005, 2009) độ phân giải 1° và phần mềm đồng hóa DIVA.
- Số liệu khảo sát Việt – Xô (1989–1995): 7.160 điểm đo theo tầng sâu, phân bố tập trung ở vùng biển phía Đông Vịnh Bắc Bộ.
- Số liệu đo đạc CTD của Viện Nghiên cứu Hải sản (2006–2013): 30.796 điểm đo với độ phân giải thẳng đứng 1 m/tầng, tập trung dày đặc ở vùng biển ven bờ phía Tây vịnh.
| Nguồn số liệu |
Thời gian |
Quy mô / Số lượng điểm đo |
Đặc điểm phân bố |
| MODAS (SST) |
1993–2013 |
Lưới 1/8° x 1/8°, trung bình tháng |
Toàn bộ bề mặt Vịnh Bắc Bộ |
| WOA & DIVA |
2001, 2005, 2009 |
Lưới 1° x 1° |
Toàn Biển Đông và Vịnh Bắc Bộ |
| Khảo sát Việt – Xô |
1989–1995 |
7.160 điểm đo tầng sâu |
Phía Đông vịnh (tầng đo 10–15 m) |
| CTD (Viện NCHS) |
2006–2013 |
30.796 điểm đo thẳng đứng |
Phía Tây vịnh (bước đo 1 m độ sâu) |
Phương pháp phân tích số liệu và kỹ thuật LOWESS
Luận án đánh giá các phương pháp biểu diễn đường phân bố thẳng đứng (profile) truyền thống:
- Phương pháp hàm vạn năng: Dựa trên giả định tựa đồng dạng $\theta = \frac{T_0 - T_z}{T_0 - T_H} = f(\eta)$, bộc lộ sai số lớn khi lớp đồng nhất nhiệt $h$ thay đổi liên tục trong vùng nước nông.
- Phương pháp hàm Spline: Bắt buộc quy đổi số liệu về các tầng chuẩn cố định, làm mất mát thông tin thực đo chi tiết và gây sai số tại vùng có gradien ngang lớn.
- Phương pháp đa thức bậc cao (bậc 5, bậc 6): Dễ phát sinh dao động không kiểm soát ở hai đầu mút (biên trên và biên đáy).
Để khắc phục, tác giả áp dụng phương pháp bình phương tối thiểu có trọng số địa phương (LOWESS - Locally Weighted Scatter plot Smooth, Cleveland 1979) lập trình trên phần mềm MATLAB. Phương pháp cho phép xây dựng mặt cong phân bố nhiệt độ và độ muối theo độ sâu cho từng nút lưới, liên kết trực tiếp với trường biến thiên bề mặt từ MODAS.
Trọng số $w_{kl}$ của một điểm $p_{kl}(x_{kl}, y_{kl})$ đối với điểm tính toán $p_{ij}(x_{ij}, y_{ij})$ trong lân cận chọn trước được xác định qua hàm tam thức bậc ba (tricube function):
$$w_{kl} = T\left(\frac{|p_{ij} - p_{kl}|}{d_{ij}}\right)$$
Trong đó $d_{ij}$ là khoảng cách từ điểm xa nhất trong lân cận đến $p_{ij}$, và hàm trọng số $T(u)$ định nghĩa:
$$T(u) = \begin{cases} (1 - |u|^3)^3 & \text{khi } |u| \le 1 \ 0 & \text{khi } |u| > 1 \end{cases}$$
Mô hình thủy động lực học 3 chiều GHER
Mô hình GHER (Geo-Hydrodynamic Environment Research) giải hệ phương trình nguyên thủy 3 chiều thủy nhiệt động lực học trong hệ tọa độ cong $\sigma$ bám đáy. Mô hình tính toán tiến triển của các biến trạng thái: vận tốc dòng chảy 3 chiều $(u, v, w)$, nhiệt độ $(T)$, độ muối $(S)$ và dao động mực nước mặt biển $(\zeta)$ bằng phương pháp thể tích hữu hạn (Finite Volume Method), tính đến ứng suất gió bề mặt, ứng suất đáy, trao đổi nhiệt – ẩm và tác động của triều.
Quy trình liên kết thuận nghịch
Hệ thống tính toán được tổ chức theo quy trình tương tác hai chiều khép kín:
- Quá trình thuận: Modul LOWESS tiếp nhận cơ sở dữ liệu đa nguồn để tính toán, tái lập trường phân bố 3D nhiệt – muối ban đầu cung cấp cho mô hình GHER.
- Quá trình nghịch: Kết quả mô phỏng trung gian từ mô hình GHER được so sánh trực tiếp với trường phân tích ban đầu. Nếu sai số vượt quá ngưỡng tiêu chuẩn:
- Nhiệt độ: chênh lệch $> 0,5^\circ\text{C}$
- Độ muối: chênh lệch $> 0,2‰$
- Vận tốc dòng chảy: chênh lệch $> 5\text{ cm/s}$
thì các giá trị tính toán sẽ được bổ sung, trộn ngược vào cơ sở dữ liệu để chạy lại modul LOWESS, tái thiết lập trường ban đầu cho đến khi đạt độ hội tụ ổn định.
Nội dung chính theo từng chương
Chương 1: Tổng quan các nghiên cứu về các trường thủy văn và hoàn lưu khu vực Vịnh Bắc Bộ và Biển Đông
Chương 1 hệ thống hóa toàn diện các công trình nghiên cứu trong nước và quốc tế về hoàn lưu và đặc trưng nhiệt – muối tại Biển Đông và Vịnh Bắc Bộ. Tác giả phân tích chi tiết các kết quả khảo sát từ giai đoạn 1945 đến 2016, phân định rõ các trường phái quan điểm về cấu trúc dòng chảy mùa hè trong vịnh. Nội dung chương khẳng định việc thiếu vắng một trường nhiệt – muối ban đầu có độ phân giải cao và thiếu sự phối hợp giữa số liệu thực đo với mô hình động lực học là nguyên nhân chính dẫn đến các kết luận không đồng nhất giữa các tác giả trước đây. Từ đó, chương đề xuất sơ đồ quy trình liên kết phân tích – mô phỏng thuận nghịch làm cơ sở phương pháp luận cho toàn bộ luận án.
Chương 2: Modul phân tích số liệu và mô hình thủy động lực 3 chiều GHER
Chương 2 trình bày chi tiết về cấu trúc cơ sở dữ liệu và thuật toán xử lý. Tác giả phân tích các tập dữ liệu CTD (30.796 điểm), số liệu Việt – Xô (7.160 điểm) và MODAS. Chương đi sâu so sánh định lượng giữa các phương pháp xây dựng profile nhiệt – muối (hàm vạn năng, Spline, đa thức bậc 5–6) và chứng minh ưu thế của kỹ thuật LOWESS trong việc tái lập mặt cong 3D.
Dưới đây là bảng trích xuất số liệu phân bố suy giảm nhiệt độ ($\Delta T$) và gia tăng độ muối ($\Delta S$) theo độ sâu tại các khu vực đại diện được tính toán trong chương:
| Tầng sâu (m) |
Suy giảm $\Delta T$ tháng 1 ($21^\circ\text{N} \pm 0,5^\circ$) |
Phân bố nhiệt độ $T$ ($^\circ\text{C}$) |
Gia tăng $\Delta S$ tháng 1 ($20^\circ\text{N} \pm 0,5^\circ$) |
Phân bố độ muối $S$ (‰) |
| 0 |
0,00 |
20,40 |
0,00 |
32,20 |
| 10 |
-0,06 |
20,34 |
+0,05 |
32,25 |
| 20 |
-0,18 |
20,22 |
+0,12 |
32,32 |
| 30 |
-0,45 |
19,95 |
+0,25 |
32,45 |
| 40 |
-1,10 |
19,30 |
+0,42 |
32,62 |
| 50 |
-1,85 |
18,55 |
+0,60 |
32,80 |
Kết quả tính toán sai số trong chương cho thấy trường nhiệt – muối tái lập bằng LOWESS bám sát số liệu thực đo tại các trạm kiểm chứng, loại bỏ hiện tượng dao động biên của đa thức bậc cao và khắc phục tình trạng thiếu hụt dữ liệu ở các vùng biển sâu hơn 50 m. Chương cũng thiết lập cấu trúc toán học của mô hình GHER, lưới tính theo phương đứng với hệ tọa độ $\sigma$, các điều kiện biên hở, thông lượng gió và thông lượng sông.
Chương 3: Ứng dụng mô hình 3D nghiên cứu biến động các trường nhiệt - muối và hoàn lưu Vịnh Bắc Bộ
Chương 3 triển khai ứng dụng hệ thống mô hình đã thiết lập để tính toán chi tiết cấu trúc hoàn lưu và trường nhiệt – muối theo 12 tháng trong năm, tập trung phân tích các tháng điển hình (tháng 1, 2, 5, 7, 8, 9, 10) tại các tầng mặt, 10 m, 30 m và 50 m.
Về cấu trúc hoàn lưu:
- Thời kỳ mùa đông (tháng 1, tháng 2): Toàn vịnh hình thành một hoàn lưu xoáy thuận quy mô lớn bao trùm. Dòng chảy tầng mặt dọc bờ biển Việt Nam chảy mạnh về hướng Nam với vận tốc đạt từ 25–40 cm/s dưới tác dụng của gió mùa Đông Bắc. Dòng chảy này được bù trừ bởi dòng chảy hướng Bắc ở phía bờ Tây đảo Hải Nam, tạo nên cấu trúc xoáy thuận khép kín. Nước từ Biển Đông chảy vào vịnh qua eo biển Quỳnh Châu với lưu tốc đáng kể.
- Thời kỳ mùa hè (tháng 7, tháng 8): Bức tranh hoàn lưu có sự phân hóa không gian phức tạp. Ở khu vực phía Bắc vịnh (từ khoảng $19^\circ\text{N}$ trở lên phía Bắc), mô hình ghi nhận sự tồn tại của xoáy thuận cục bộ chịu tác động kết hợp của dòng chảy mật độ và địa hình. Ở khu vực phía Nam vịnh, hình thành cấu trúc xoáy nghịch mở ra phía Biển Đông, phù hợp với trường ứng suất gió mùa Tây Nam và kết quả mô phỏng của Gao và cộng sự (2013).
- Các giai đoạn chuyển tiếp (tháng 5, tháng 9, tháng 10): Hoàn lưu chuyển biến dần trạng thái, xuất hiện các tâm xoáy quy mô nhỏ cục bộ dọc theo sườn dốc địa hình giữa vịnh.
Về cấu trúc nhiệt độ và độ muối:
- Trường nhiệt độ: Vào mùa đông, nước biển bị xáo trộn đối lưu mạnh do gió mùa, tầng đồng nhất nhiệt phát triển sâu xuống 30–40 m, chênh lệch nhiệt độ tầng mặt và đáy nhỏ ($< 1,5^\circ\text{C}$). Nhiệt độ nước tăng dần từ Tây Bắc sang Đông Nam vịnh ($17^\circ\text{C}$ lên $22^\circ\text{C}$). Vào mùa hè, hiện tượng phân tầng nhiệt xuất hiện rõ rệt, lớp nêm nhiệt hình thành mạnh ở độ sâu từ 15 m đến 35 m; nhiệt độ nước tầng mặt đạt $29–30,5^\circ\text{C}$, trong khi tầng đáy ($> 50\text{ m}$) duy trì khối nước lạnh cục bộ $22–24^\circ\text{C}$.
- Trường độ muối: Phân bố độ muối phản ánh rõ ảnh hưởng của dòng xả từ các cửa sông lớn (sông Hồng, sông Mã) vào mùa kiệt (độ muối vùng cửa sông 28–30‰) và mùa lũ (độ muối giảm xuống dưới 20–25‰ ở dải ven bờ Tây). Nước mặn từ Biển Đông ($> 33,5‰$) xâm nhập sâu vào đáy vịnh qua cửa phía Nam và eo Quỳnh Châu.
Các mặt cắt dọc kinh tuyến (107°E, 108°E) và vĩ tuyến (18°N, 19°N, 20°N) trong chương thể hiện rõ ranh giới các khối nước và quá trình xâm nhập mặn đáy theo mùa.
Kết quả và những đóng góp mới
Những đóng góp mới về mặt khoa học và lý luận
- Hoàn thiện phương pháp luận phân tích số liệu: Đã nghiên cứu, phát triển và ứng dụng thành công kỹ thuật bình phương tối thiểu có trọng số địa phương (LOWESS) kết hợp với thuật toán tam thức bậc ba để tích hợp các tập dữ liệu đa nguồn (viễn thám MODAS, số liệu đo sâu CTD độ phân giải cao và số liệu khảo sát lịch sử Việt – Xô). Kỹ thuật này cho phép tái lập mặt cong phân bố 3D nhiệt độ và độ muối liên tục theo không gian – thời gian, hạn chế sai số do tính không đồng bộ của số liệu đo đạc.
- Lý giải bản chất hoàn lưu mùa hè Vịnh Bắc Bộ: Kết quả mô phỏng 3 chiều đã giải thích được tính hai mặt của hoàn lưu mùa hè mà các nghiên cứu trước đây tranh luận: khu vực phía Bắc vịnh tồn tại cấu trúc xoáy thuận do vai trò chi phối của dòng mật độ và phân bố địa hình nông, trong khi khu vực phía Nam vịnh phát triển cấu trúc xoáy nghịch dưới tác động của trường gió mùa Tây Nam.
- Làm rõ bức tranh phân tầng thủy văn 3D: Luận án đã mô tả chi tiết quy luật biến đổi theo mùa của cấu trúc phân tầng nhiệt – muối, xác định độ sâu lớp nêm nhiệt mùa hè (15–35 m), vùng nước lạnh đáy ở trung tâm vịnh và mức độ xáo trộn đối lưu sâu vào mùa đông.
Đóng góp thực tiễn và giải pháp đề xuất
- Thiết lập thành công quy trình kết hợp hai chiều thuận nghịch giữa modul phân tích số liệu và mô hình số trị 3D GHER, tạo ra công cụ tạo trường ban đầu có độ tin cậy cao phục vụ cho hệ thống mô hình dự báo hải văn và môi trường biển nghiệp vụ.
- Cung cấp bộ bản đồ chi tiết về dòng chảy, nhiệt độ và độ muối trung bình 12 tháng tại các tầng chuẩn (0 m, 10 m, 30 m, 50 m) làm tài liệu tham khảo chuẩn xác cho các bài toán tính toán lan truyền chất ô nhiễm, quy hoạch khai thác hải sản và an ninh quốc phòng trên Vịnh Bắc Bộ.
Hạn chế và hướng nghiên cứu tiếp
Hạn chế của luận án
- Cơ sở dữ liệu thực đo tầng sâu phân bố không đồng đều: số liệu CTD tập trung chủ yếu ở dải ven bờ Tây, trong khi số liệu phía Đông chủ yếu dựa vào các chuyến khảo sát cũ (Việt – Xô trước năm 1995), làm hạn chế mức độ kiểm chứng tức thời ở một số khu vực nước sâu giáp đảo Hải Nam.
- Mô hình chưa tích hợp liên kết thời gian thực với mô hình sóng biển và mô hình khí tượng phân giải cao trong điều kiện có bão mạnh hoặc áp thấp nhiệt đới dị thường; việc tham số hóa dòng chảy sông chủ yếu dựa trên lưu lượng trung bình tháng.
Hướng nghiên cứu tiếp theo
- Nâng cao độ phân giải không gian của mô hình 3D GHER xuống quy mô dưới 1/24° kinh vĩ, lồng ghép chi tiết các vùng cửa sông ven bờ.
- Mở rộng kết nối quy trình phân tích – mô phỏng thuận nghịch này với các mô hình lan truyền sinh thái, vận chuyển bùn cát và ô nhiễm dầu phục vụ công tác cảnh báo sự cố môi trường biển.
Giá trị tham khảo
Luận án là tài liệu chuyên khảo có giá trị cao đối với nhiều nhóm đối tượng:
- Nghiên cứu sinh, học viên cao học và giảng viên: Chuyên ngành Hải dương học, Khí tượng học, Thủy văn học và Mô hình hóa môi trường có thể tham khảo thuật toán ứng dụng LOWESS trong xử lý số liệu hải văn phân tán và phương pháp cài đặt hệ phương trình nguyên thủy trong mô hình GHER (Chương 2).
- Cơ quan dự báo khí tượng thủy văn và môi trường: Trung tâm Dự báo Khí tượng Thủy văn Quốc gia, Viện Hải dương học, Viện Nghiên cứu Biển và Hải đảo có thể sử dụng quy trình tính toán thuận nghịch và bộ trường chế độ nhiệt – muối 12 tháng để làm trường ban đầu/trường biên cho các mô hình dự báo nghiệp vụ.
- Ngành thủy sản và kinh tế biển: Viện Nghiên cứu Hải sản, các cơ quan quản lý khai thác thủy sản có thể khai thác các dữ liệu về cấu trúc khối nước lạnh tầng đáy, lớp nêm nhiệt mùa hè và vùng phân bố dòng chảy cục bộ (Chương 3) để xác định ngư trường và mùa vụ đánh bắt.
Câu hỏi thường gặp
1. Điểm khác biệt cơ bản giữa kỹ thuật LOWESS và phương pháp hàm Spline hay đa thức bậc cao trong xử lý số liệu nhiệt – muối là gì?
Hàm Spline đòi hỏi quy đổi dữ liệu về tầng chuẩn cố định làm mất chi tiết đo đạc 1 m của CTD, còn đa thức bậc cao dễ bị dao động lớn ở hai đầu mút (mặt biển và đáy). Kỹ thuật LOWESS sử dụng lân cận cục bộ có trọng số tam thức bậc ba, cho phép xây dựng mặt cong biến thiên liên tục trực tiếp từ mọi điểm đo thực tế không đồng bộ, gắn liền với trường nhiệt độ bề mặt MODAS mà không làm biến dạng cấu trúc phân tầng.
2. Luận án giải quyết mâu thuẫn giữa hai quan điểm xoáy thuận và xoáy nghịch mùa hè ở Vịnh Bắc Bộ như thế nào?
Luận án chỉ ra rằng hoàn lưu mùa hè không tồn tại đơn nhất một xoáy toàn vịnh mà phân hóa thành hai cấu trúc rõ rệt: nửa phía Bắc vịnh (từ vĩ tuyến 19°N trở lên) duy trì xoáy thuận do dòng mật độ kết hợp địa hình chi phối; trong khi nửa phía Nam vịnh hình thành xoáy nghịch mở rộng ra Biển Đông dưới tác động của trường gió mùa Tây Nam.
3. Ngưỡng sai số nào được sử dụng để kiểm soát quy trình lặp thuận nghịch giữa modul LOWESS và mô hình GHER?
Hệ thống sử dụng các ngưỡng sai số giới hạn: độ lệch nhiệt độ $\le 0,5^\circ\text{C}$, độ lệch độ muối $\le 0,2‰$ và độ lệch vận tốc dòng chảy $\le 5\text{ cm/s}$. Khi kết quả tính toán trung gian từ GHER vượt các ngưỡng này so với trường ban đầu, dữ liệu sẽ được phản hồi ngược để tái phân tích LOWESS.
4. Khối lượng số liệu thực nghiệm được xử lý trong luận án gồm những nguồn nào?
Tổng khối lượng số liệu bao gồm: chuỗi ảnh viễn thám SST MODAS trung bình tháng giai đoạn 1993–2013 (độ phân giải 1/8°); 7.160 điểm đo sâu từ các chuyến khảo sát hợp tác Việt – Xô (1989–1995); 30.796 điểm đo CTD thẳng đứng của Viện Nghiên cứu Hải sản (2006–2013) cùng các tập dữ liệu WOA và DIVA.
5. Cấu trúc nhiệt độ tầng đáy Vịnh Bắc Bộ vào mùa hè có đặc điểm gì nổi bật?
Vào mùa hè, dù nhiệt độ tầng mặt tăng cao ($29–30,5^\circ\text{C}$), ở các vùng trũng sâu trung tâm vịnh ($> 50\text{ m}$) vẫn tồn tại khối nước lạnh tầng đáy với nhiệt độ duy trì từ $22–24^\circ\text{C}$, được phân cách với tầng mặt bởi lớp nêm nhiệt sắc nét nằm ở độ sâu khoảng 15–35 m.
Kết luận
Luận án tiến sĩ của tác giả Hà Thanh Hương đã phát triển và thực nghiệm thành công quy trình công nghệ kết hợp modul phân tích số liệu LOWESS với mô hình thủy nhiệt động lực 3 chiều GHER cho vùng biển Vịnh Bắc Bộ. Công trình đã tái hiện và làm sáng tỏ quy luật biến động theo mùa của các trường nhiệt độ, độ muối và cấu trúc hoàn lưu 3D, đặc biệt là bản chất phân hóa hoàn lưu trong mùa hè và cấu trúc phân tầng nhiệt – muối 12 tháng trong năm. Các kết quả đạt được cung cấp cơ sở khoa học và bộ công cụ đầu vào tin cậy cho các bài toán nghiên cứu hải văn, giám sát môi trường và dự báo nghiệp vụ biển tại Việt Nam.
