Tổng quan về luận án

Vùng biển ven bờ giữ vị thế chiến lược đặc biệt trong cơ cấu phát triển kinh tế - xã hội toàn cầu, là nơi tập trung khoảng 60% dân số thế giới, đóng góp 25% tổng năng suất sơ cấp toàn cầu và cung cấp tới 90% sản lượng khai thác thủy hải sản tự nhiên. Tại Việt Nam, với hơn 3.260 km bờ biển và không gian kinh tế biển rộng lớn, Nghị quyết số 36-NQ/TW về Chiến lược phát triển bền vững kinh tế biển Việt Nam đến năm 2030, tầm nhìn đến năm 2045 đã xác định nhiệm vụ bảo vệ môi trường biển là nội dung xuyên suốt. Vùng biển ven bờ Nam Trung Bộ và Nam Bộ (đặc biệt là dải ven biển từ Khánh Hòa đến Ninh Thuận) là trung tâm kinh tế biển tổng hợp về hàng hải, khai khoáng, khai thác và nuôi trồng thủy sản. Tuy nhiên, áp lực từ phát triển đô thị, xả thải công nghiệp, nước thải nuôi trồng thủy sản và dòng chảy lục địa đã làm gia tăng nguy cơ suy thoái môi trường cục bộ, điển hình là hiện tượng phú dưỡng và tảo nở hoa độc hại từng gây thiệt hại nghiêm trọng (như sự cố tảo nở hoa tại vịnh Vân Phong cuối năm 2016).

Khoảng trống nghiên cứu (research gap) cốt lõi xuất phát từ sự hạn chế của mạng lưới quan trắc truyền thống tại Việt Nam: mật độ trạm mỏng (chỉ khoảng 0,03 điểm/km bờ biển) với tần suất thấp (2–3 đợt/năm đối với cấp quốc gia và hàng tháng đối với cấp tỉnh), không thể bao quát quy mô không gian lớn cũng như không bắt kịp diễn biến nhanh của các sự cố môi trường biển. Mặt khác, các nghiên cứu viễn thám màu đại dương trước đây ở Việt Nam chủ yếu dựa trên tư liệu MODIS độ phân giải không gian trung bình (1 km) hoặc MERIS đã dừng hoạt động từ 2012; việc ứng dụng thế hệ vệ tinh mới Sentinel-3 với bộ cảm OLCI (Ocean and Land Colour Instrument, độ phân giải không gian 300 m) chưa từng được thực hiện đồng bộ. Hơn thế, việc xử lý tín hiệu quang học tại vùng nước ven bờ có đặc tính quang học phức tạp (nước Case-2) dưới ảnh hưởng của khí hậu nhiệt đới gió mùa có độ che phủ mây dày đặc đặt ra thách thức khoa học lớn chưa có lời giải toàn diện.

Luận án tiến sĩ "Nghiên cứu giám sát chất lượng môi trường nước biển ven bờ phía Nam bằng công nghệ viễn thám" của nghiên cứu sinh Nguyễn An Bình (chuyên ngành Bản đồ, viễn thám và hệ thông tin địa lý, Mã số: 9 44 02 14, Học viện Khoa học và Công nghệ – Viện Hàn lâm Khoa học và Công nghệ Việt Nam, dưới sự hướng dẫn của PGS. Phạm Việt Hòa và TS. Hồ Đình Duẩn) giải quyết trực diện bài toán này. Nghiên cứu xác lập các câu hỏi và giả thuyết khoa học trọng tâm:

  1. Câu hỏi nghiên cứu 1 (RQ1): Mô hình hiệu chỉnh khí quyển nào tối ưu hóa việc bóc tách phổ phản xạ rời mặt nước ($p_w$) cho dữ liệu Sentinel-3 OLCI tại vùng biển ven bờ phức tạp phía Nam Việt Nam?
  2. Câu hỏi nghiên cứu 2 (RQ2): Thuật toán quang sinh học biển nào cho phép định lượng chính xác hàm lượng diệp lục tố a (chlorophyll-a - Chl-a) trong các phân vùng quang học nước biển ven bờ Case-2?
  3. Câu hỏi nghiên cứu 3 (RQ3): Làm thế nào để đồng bộ hóa nguồn dữ liệu vệ tinh đa cảm biến và giải quyết triệt để khiếm khuyết mất dữ liệu do mây che phủ trong chế độ gió mùa nhiệt đới?

Hệ thống giả thuyết tương ứng:

  • Giả thuyết 1 (H1): Thuật toán hiệu chỉnh khí quyển dựa trên mạng nơ-ron C2RCC (Case 2 Regional Coast Colour) vượt trội hơn phương pháp DSF (Dark Spectrum Fitting) trong việc khử nhiễu sol khí (aerosol) và hiệu ứng liền kề tại vùng nước Case-2 nhiệt đới.
  • Giả thuyết 2 (H2): Thuật toán bán thực nghiệm đa kênh dải chuyển tiếp xanh lam - xanh lục (blue-to-green) chuyên biệt cho nước ven bờ (như OC5v7) mang lại độ chính xác cao hơn thuật toán OC4ME truyền thống khi so sánh với số liệu thực địa ($p < 0,01$).
  • Giả thuyết 3 (H3): Tích hợp dữ liệu cặp vệ tinh Sentinel-3A/3B kết hợp thuật toán tái cấu trúc dữ liệu không gian DINEOF (Data Interpolating Empirical Orthogonal Functions) cho phép tạo lập chuỗi dữ liệu giám sát Chl-a liên tục hàng ngày với độ tin cậy thống kê cao.

Khung lý thuyết của nghiên cứu dựa trên lý thuyết truyền xạ bức xạ quang học hải dương, phân loại quang sinh học nước biển Case-1/Case-2 và toán học hàm trực giao thực nghiệm EOF. Về phạm vi, luận án tập trung vào dải ven biển từ Khánh Hòa đến Ninh Thuận (bao gồm vịnh Vân Phong, vịnh Nha Phu, vịnh Nha Trang, vịnh Cam Ranh và vịnh Phan Rang) trong giai đoạn 2018–2019, thiết lập một quy trình chuẩn hóa từ xử lý ảnh vệ tinh đến thành lập bản đồ phân vùng dinh dưỡng phục vụ giám sát môi trường biển.

Literature Review và Positioning

Lịch sử phát triển của viễn thám màu đại dương (ocean color remote sensing) trải qua hơn 4 thập kỷ kể từ khi vệ tinh CZCS (Coastal Zone Color Scanner, 1978–1986) chứng minh tính khả thi của việc ước tính sắc tố thực vật phù du từ vũ trụ (Gordon, 1978; Gordon et al., 1980). Dòng nghiên cứu này phát triển mạnh mẽ qua các thế hệ cảm biến SeaWiFS (1997), MODIS Aqua/Terra (2002), MERIS (2002–2012), VIIRS (2011) và đỉnh cao hiện nay là Sentinel-3 OLCI (Sentinel-3A phóng năm 2016, Sentinel-3B phóng năm 2018).

Trong y văn quốc tế, tồn tại những tranh luận lý thuyết và phương pháp luận sâu sắc giữa hai trường phái chính:

  1. Trường phái thuật toán tỷ số kênh truyền thống (Standard Empirical Ocean Color Algorithms): Đại diện bởi O'Reilly et al. (1998, 2019) và nhóm nghiên cứu NASA Ocean Biology Processing Group (OBPG), ủng hộ mô hình đa thức bậc 4 dựa trên tỷ số cực đại của các kênh xanh lam/xanh lục (như OC2, OC3, OC4, OC4ME, OC6). Nhóm này lập luận rằng các thuật toán thực nghiệm được hiệu chuẩn từ cơ sở dữ liệu toàn cầu SeaBASS (Werdell et al., 2003, 2019) có tính tổng quát hóa cao và độ ổn định thống kê tốt trên các vùng đại dương mở (nước Case-1).
  2. Trường phái mô hình bán giải tích và trí tuệ nhân tạo chuyên biệt cho nước đục ven bờ (Semi-analytical & Machine Learning Algorithms for Case-2 Waters): Đại diện bởi Doerffer & Schiller (2007) với mô hình mạng nơ-ron C2RCC, Gohin et al. (2002) với thuật toán OC5, Mishra & Mishra (2012) với chỉ số diệp lục chuẩn hóa NDCI (Normalized Difference Chlorophyll Index), và Hu et al. (2012) với chỉ số màu CI (Color Index). Các tác giả này chỉ ra rằng tại vùng nước Case-2 (chiếm 2% diện tích đại dương nhưng tập trung phần lớn hoạt động nhân sinh), sự hiện diện độc lập của vật chất lơ lửng vô cơ (TSM/TSS) và chất hữu cơ hòa tan có màu (CDOM) làm vô hiệu hóa mối tương quan đơn nhất giữa Chl-a và phổ phản xạ blue-to-green, đòi hỏi các thuật toán phải bóc tách riêng biệt các thành phần quang học nội tại (IOP).

Về hiệu chỉnh khí quyển, mâu thuẫn học thuật nảy sinh giữa phương pháp trừ phổ tối đa kênh DSF trong Acolite (Vanhellemont, 2019) vốn giả định tín hiệu cận hồng ngoại/hồng ngoại sóng ngắn (NIR/SWIR) hoàn toàn bằng không trên nước trong, và phương pháp C2RCC sử dụng mạng nơ-ron lan truyền ngược đa lớp được huấn luyện trên mô sở dữ liệu mô phỏng truyền xạ bức xạ Hydrolight khổng lồ để giải quyết đồng thời tán xạ sol khí và phản xạ rời mặt nước.

Luận án định vị chính xác khoảng trống học thuật: Hầu hết các thuật toán chuẩn quốc tế được xây dựng từ tập dữ liệu ôn đới hoặc đại dương hở. Khi áp dụng vào vùng biển nhiệt đới gió mùa ven bờ Nam Trung Bộ Việt Nam – nơi có sự phân tầng thủy văn mạnh mẽ, hiện tượng nước trồi mùa hè đặc thù, cùng lượng trầm tích lục địa cao – các mô hình này bộc lộ sai số lớn do hiệu ứng mờ liền kề từ đất liền và nhiễu do mây cục bộ.

So sánh với các nghiên cứu quốc tế điển hình:

  • Tại Vịnh Chesapeake (Hoa Kỳ): Nghiên cứu của Werdell et al. (2009) và Son & Wang (2012) cho thấy thuật toán chuẩn của MODIS đánh giá vượt mức (overestimate) nồng độ Chl-a từ 2–5 lần do ảnh hưởng của CDOM từ cửa sông.
  • Tại Biển Baltic (Châu Âu): Các nghiên cứu của Attila et al. (2013) và Toming et al. (2016) trên dữ liệu Sentinel-2 MSI và Sentinel-3 OLCI khẳng định thuật toán C2RCC và NDCI ưu việt hơn các thuật toán NASA OCx trong môi trường nước có độ đục và hàm lượng hữu cơ cao.

Luận án của nghiên cứu sinh Nguyễn An Bình kế thừa các khung tiếp cận tiên tiến trên thế giới, nhưng tạo bước đột phá trong việc tối ưu hóa và đánh giá chéo đồng thời các mô hình C2RCC, DSF, OC5, OC4ME, NDCI, CI trên nguồn dữ liệu đa phổ Sentinel-3A/3B tại vùng biển Việt Nam.

Đóng góp lý thuyết và khung phân tích

Đóng góp cho lý thuyết

Nghiên cứu mang lại những đóng góp nền tảng cho lý thuyết quang học hải dương và viễn thám môi trường biển:

  1. Mở rộng lý thuyết phân loại quang học nước biển Morel & Prieur (1977) trong điều kiện nhiệt đới ven bờ: Luận án làm sáng tỏ ranh giới quang học giữa nước Case-1 và Case-2 tại vùng biển Nam Trung Bộ, chứng minh rằng sự biến thiên theo mùa của trường bức xạ quang sinh học chịu sự chi phối phức tạp bởi tương tác ba thành phần: sắc tố phytoplankton ($P$), chất hữu cơ hòa tan có màu CDOM ($Y$) và vật chất lơ lửng trầm tích ($S$).
  2. Phát triển mô hình truyền xạ bức xạ vi sai đỉnh khí quyển: Luận án kế thừa và chuẩn hóa công thức phân tách tổng hệ số phản xạ tại đỉnh khí quyển $p_t(\lambda)$ theo Gordon (1997) và Mobley (1994): $$p_t(\lambda) = p_a(\lambda) + p_r(\lambda) + p_{ra}(\lambda) + T(\lambda)p_g(\lambda) + t_0 \cdot t_v \cdot p_w(\lambda)$$ Trong đó bóc tách thành công thành phần phản xạ sol khí $p_a(\lambda)$, tán xạ phân tử Rayleigh $p_r(\lambda)$, tương tác Rayleigh - sol khí $p_{ra}(\lambda)$, phản xạ tia sáng mặt trời $p_g(\lambda)$ (sun glint), bọt trắng sóng biển (whitecaps) và hệ số lan truyền sóng khuếch tán $t_0, t_v$, từ đó thu hồi chính xác phổ phản xạ rời mặt nước $p_w(\lambda)$ – biến số biểu kiến (AOP) cốt lõi chứa đựng thông tin chất lượng nước.
  3. Định lập các mệnh đề khoa học (Propositions):
    • Mệnh đề 1: Trong điều kiện khí quyển ven biển nhiệt đới ẩm, mạng nơ-ron C2RCC cho phép xấp xỉ phi tuyến tính chính xác hơn các thông số quang học nội tại (IOP: hệ số hấp thụ $a(\lambda)$ và hệ số tán xạ ngược $b_b(\lambda)$) so với các phép ngoại suy hàm mũ aerosol truyền thống.
    • Mệnh đề 2: Việc kết hợp tỷ số đa phổ cực đại mở rộng (OC5v7) hiệu chỉnh được độ lệch do tán xạ của hạt vô cơ lơ lửng, tạo ra đường cong tương quan logarithmic tuyến tính chặt chẽ giữa tỷ số kênh và logarit nồng độ Chl-a.
graph TD
    A["Tín hiệu bức xạ đỉnh khí quyển (TOA) L_TOA"] --> B["Hiệu chỉnh khí quyển chuyên sâu"]
    B --> C["Thuật toán C2RCC (Neural Network)"]
    B --> D["Thuật toán DSF (Dark Spectrum Fitting)"]
    C --> E["Phổ phản xạ rời mặt nước p_w(λ) tại 413, 443, 490, 510, 560, 665 nm"]
    D --> E
    E --> F["Mô hình quang sinh học biển ước tính Chl-a"]
    F --> G["Thuật toán OC5v7 (Coastal Case-2)"]
    F --> H["Thuật toán OC4ME, NDCI, Color Index"]
    G --> I["Sản phẩm không gian Chl-a ban đầu"]
    H --> I
    I --> J["Xử lý mặt nạ mây và bóng mây"]
    J --> K["Mô hình tái cấu trúc dữ liệu không gian DINEOF"]
    K --> L["Chuỗi bản đồ Chl-a liên tục 300m hàng ngày"]
    L --> M["Phân vùng trạng thái dinh dưỡng (Trophic State Mapping)"]

Khung phân tích độc đáo

Khung phân tích của luận án tích hợp liên ngành ba cấu phần lý thuyết:

  1. Lý thuyết Quang học Hải dương (Ocean Optics Theory): Liên kết các đặc tính quang học vốn có (Inherent Optical Properties - IOP: hệ số hấp thụ thực vật phù du $a_{ph}$, hấp thụ CDOM $a_{cdom}$, tán xạ ngược hạt lơ lửng $b_{bp}$) với các đặc tính quang học bên ngoài (Apparent Optical Properties - AOP: phổ phản xạ rời mặt nước vi sai $R_{rs}(\lambda)$).
  2. Lý thuyết Xử lý Tín hiệu Địa không gian và Tối ưu hóa Toàn cục: Tích hợp thuật toán DINEOF (Data Interpolating Empirical Orthogonal Functions) của Beckers & Rixen (2003) và Alvera-Azcárate et al. (2005). Thuật toán sử dụng kỹ thuật phân tích thành phần chính thực nghiệm trực giao (EOF) lặp để tái tạo ma trận không gian - thời gian bị khuyết thiếu dữ liệu do mây mà không cần tham số giả định ngoại suy nhân tạo.
  3. Lý thuyết Đánh giá Trạng thái Dinh dưỡng Thủy vực Ven bờ (Coastal Trophic State Theory): Chuẩn hóa thang phân loại Carlson và tiêu chuẩn đánh giá phú dưỡng quốc tế sang điều kiện nhiệt đới Việt Nam, thiết lập ngưỡng định lượng nồng độ Chl-a để phân chia 4 mức độ: Nghèo dinh dưỡng (Oligotrophic: Chl-a $< 0,5\ \mu\text{g/L}$), Trung dưỡng (Mesotrophic: $0,5 \le \text{Chl-a} < 2,0\ \mu\text{g/L}$), Phú dưỡng (Eutrophic: $2,0 \le \text{Chl-a} < 5,0\ \mu\text{g/L}$), và Siêu phú dưỡng (Hypertrophic: Chl-a $\ge 5,0\ \mu\text{g/L}$).

Điều kiện biên (boundary conditions) của khung phân tích: Áp dụng hiệu quả tối đa cho vùng nước ven bờ có độ sâu $> 5\text{ m}$ nhằm triệt tiêu phản xạ đáy (bottom reflectance), góc thiên đỉnh mặt trời $\theta_0 < 60^\circ$, và vùng không bị che phủ hoàn toàn bởi mây dày trong chuỗi thời gian phân tích DINEOF.

Phương pháp nghiên cứu tiên tiến

Thiết kế nghiên cứu

Luận án tuân thủ triệt để thế giới quan thực chứng (positivism paradigm) với cách tiếp cận định lượng nghiêm ngặt, kết hợp đồng bộ giữa nghiên cứu thực địa (in-situ oceanographic survey), vật lý truyền xạ bức xạ và phân tích dữ liệu lớn địa không gian. Thiết kế nghiên cứu đa tầng (multi-level design) được cấu trúc qua 4 cấp độ:

  • Cấp độ 1 (Sensor Level): Thu nhận bức xạ phổ toàn phần từ hai bộ cảm OLCI trên vệ tinh Sentinel-3A và Sentinel-3B thuộc chương trình Copernicus của Cơ quan Vũ trụ Châu Âu (ESA). Bộ cảm OLCI sở hữu 21 kênh phổ (trong đó 6 kênh quan trọng phục vụ nước biển: 413, 443, 490, 510, 560, 665 nm), tỷ số tín hiệu trên nhiễu (SNR) vượt trội (SNR $> 1000$ tại dải nhìn thấy), độ phân giải không gian 300 m và độ rộng dải quét 1.270 km.
  • Cấp độ 2 (Atmospheric & Bio-optical Processing Level): Kiểm định song song các giải thuật hiệu chỉnh khí quyển (C2RCC vs. DSF) và các thuật toán quang sinh học ước tính nồng độ Chl-a.
  • Cấp độ 3 (Spatio-temporal Reconstruction Level): Đồng bộ hóa và tái cấu trúc chuỗi thời gian không gian bằng mô hình toán học DINEOF nhằm triệt tiêu khiếm khuyết mất dữ liệu do mây nhiệt đới.
  • Cấp độ 4 (Validation & Application Level): Đánh giá độ chính xác độc lập với số liệu thực địa thu thập tại các vịnh ven biển Nam Trung Bộ và xây dựng bản đồ phân vùng dinh dưỡng theo tuần/tháng.

Quy trình nghiên cứu rigorous

Quy trình nghiên cứu được thực hiện qua các bước thực nghiệm chuẩn hóa:

flowchart LR
    subgraph S1["1. Thu thập & Tiền xử lý"]
        A1["Dữ liệu vệ tinh Sentinel-3A/3B OLCI"]
        A2["Số liệu thực địa in-situ Chl-a"]
        A3["Dữ liệu khí quyển ECMWF"]
    end

    subgraph S2["2. Hiệu chỉnh & Bóc tách"]
        B1["Hiệu chỉnh khí quyển C2RCC & DSF"]
        B2["Lọc mây, bóng mây, Sun glint, Whitecaps"]
        B3["Trích xuất phổ phản xạ R_rs(λ)"]
    end

    subgraph S3["3. Mô hình hóa Quang sinh học"]
        C1["Tính toán Chl-a (OC5v7, OC4ME, NDCI, CI)"]
        C2["Kiểm định chéo & Đánh giá sai số (RMSE, R², Taylor)"]
        C3["Lựa chọn thuật toán tối ưu (OC5v7 + C2RCC)"]
    end

    subgraph S4["4. Tái cấu trúc & Ứng dụng"]
        D1["Tích hợp dữ liệu cặp vệ tinh S-3A + S-3B"]
        D2["Tái cấu trúc không gian DINEOF"]
        D3["Bản đồ phân bố Chl-a & Phân vùng dinh dưỡng"]
    end

    S1 --> S2 --> S3 --> S4

Giao thức thu mẫu và kiểm định thực địa: Các chiến dịch đo đạc thực địa được triển khai đồng bộ với thời gian vệ tinh bay qua (overpass time trong khoảng $\pm 2$ giờ) tại các điểm quan trắc trên các vịnh Vân Phong, Nha Phu, Nha Trang, Cam Ranh. Mẫu nước bề mặt ($0 - 0,5\text{ m}$) được thu thập bằng bình Nansen/Niskin chuẩn, lọc qua màng sợi thủy tinh Whatman GF/F ($0,7\ \mu\text{m}$), bảo quản đông sâu $-20^\circ\text{C}$ trong bóng tối và phân tích nồng độ Chl-a bằng phương pháp quang phổ huỳnh quang (spectrophotometric/fluorometric method) theo tiêu chuẩn ISO 10260.

Phương pháp tam giác đạc (Triangulation): Kết hợp đa nguồn dữ liệu (ảnh vệ tinh Sentinel-3A, Sentinel-3B, số liệu thực địa in-situ và dữ liệu khí quyển tái phân tích ECMWF), đa thuật toán hiệu chỉnh khí quyển, và đa mô hình quang học hải dương để đảm bảo độ tin cậy và giá trị nội suy của kết quả.

Data và phân tích

Dữ liệu ảnh vệ tinh được xử lý tự động hóa bằng ngôn ngữ Python kết hợp các công cụ chuyên dụng: phần mềm SNAP (Sentinel Application Platform) của ESA tích hợp module C2RCC, công cụ Acolite (Royal Belgian Institute of Natural Sciences) cho thuật toán DSF, và thư viện khoa học NumPy, SciPy, Matplotlib để lập trình mô hình giải thuật DINEOF.

Mô hình toán học DINEOF thiết lập trên thuật toán lặp phân tích kỳ dị (Singular Value Decomposition - SVD): $$X = U \cdot S \cdot V^T$$ Trong đó ma trận dữ liệu không gian - thời gian $X$ (kích thước $M \times N$) chứa các điểm ảnh bị mây che phủ được nội suy lặp qua các mode trực giao EOF tối ưu nhất, được kiểm soát tự động thông qua sai số kiểm định chéo điểm ẩn (cross-validation error) $RMSE_{\text{DINEOF}}$.

Đánh giá độ chính xác mô hình được thực hiện thông qua các chỉ số thống kê chuẩn mực: Hệ số xác định ($R^2$), Sai số căn bậc hai trung bình (Root Mean Square Error - $RMSE$), Sai số tuyệt đối trung bình (Mean Absolute Error - $MAE$), hệ số tương quan tuyến tính Pearson ($r$), và biểu đồ Taylor Diagram để tổng hợp độ lệch chuẩn và độ tương quan không gian. Mọi phân tích thống kê đều được kiểm định mức ý nghĩa với $p\text{-value} < 0,01$.

Phát hiện đột phá và implications

Những phát hiện then chốt

  1. Xác lập hiệu năng vượt trội của mô hình C2RCC trong hiệu chỉnh khí quyển vùng ven bờ: Kết quả thực nghiệm trên ảnh Sentinel-3 OLCI chỉ ra rằng thuật toán C2RCC bóc tách phổ phản xạ rời mặt nước ($p_w$) tại các kênh 413, 443, 490, 510, 560, 665 nm đạt độ tương thích cao, triệt tiêu hiện tượng giá trị phản xạ âm thường gặp ở các kênh bước sóng ngắn (413, 443 nm) của thuật toán DSF khi gặp nước đục ven bờ. Sai số hiệu chỉnh khí quyển của C2RCC thấp hơn đáng kể và duy trì tính nhất quán quang học giữa hai bộ cảm Sentinel-3A và 3B ($p < 0,01$).

  2. Mô hình OC5v7 đạt độ chính xác cao nhất trong ước tính nồng độ Chl-a vùng biển Case-2: Đánh giá chéo 4 thuật toán quang sinh học (OC5v7, OC4ME, NDCI, Color Index) trên nền ảnh hiệu chỉnh C2RCC cho thấy thuật toán OC5v7 chuyên biệt cho vùng nước ven bờ đạt độ chính xác cao nhất với $RMSE$ thấp nhất và hệ số tương quan $R^2 > 0,82$ so với dữ liệu thực địa. Thuật toán OC4ME tiêu chuẩn của đại dương mở bị sai lệch đáng kể (đánh giá cao hơn thực tế 30–60%) do không tách được tín hiệu tán xạ của vật chất lơ lửng vô cơ.

graph LR
    subgraph HCKQ["Hiệu chỉnh khí quyển"]
        A["C2RCC (Tối ưu)"] -->|Khử nhiễu sol khí & hấp thụ| C["Phổ phản xạ R_rs"]
        B["DSF (Nhiễu kênh ngắn)"] -.->|Sai số tại 413-443nm| C
    end
    subgraph QSH["Mô hình Quang sinh học"]
        C --> D["OC5v7 (R² > 0.82, RMSE tối ưu)"]
        C --> E["OC4ME (Overestimate 30-60%)"]
    end
    subgraph TK["Tái cấu trúc không gian"]
        D --> F["DINEOF + Cặp S-3A/S-3B"]
        F --> G["Phục hồi 75-85% điểm ảnh mất do mây"]
    end

  1. Bước đột phá trong xử lý mất dữ liệu do mây bằng mô hình DINEOF kết hợp cặp vệ tinh Sentinel-3A/3B: Lần đầu tiên tại Việt Nam, sự kết hợp đồng thời nguồn ảnh Sentinel-3A và Sentinel-3B (rút ngắn chu kỳ chụp lặp xuống 1–2 ngày) cùng thuật toán DINEOF đã phục hồi thành công từ 75% đến 85% diện tích mặt nước bị che phủ bởi mây và bóng mây trong các tháng mùa mưa bão (tháng 7, 8 và 10). Phân bố không gian của trường Chl-a sau tái cấu trúc hoàn toàn liền mạch, bảo toàn đặc tính động lực học của các xoáy nước và vệt nước trồi ven bờ.

  2. Bản đồ hóa quy luật phân bố Chl-a và trạng thái dinh dưỡng chi tiết theo tuần: Phân tích chuỗi dữ liệu tái cấu trúc giai đoạn tháng 7 – tháng 8 năm 2019 cho thấy hàm lượng Chl-a trung bình tuần biến động mạnh từ $0,2\ \mu\text{g/L}$ đến trên $8,5\ \mu\text{g/L}$. Các khu vực vịnh kín có mật độ nuôi trồng thủy sản lồng bè dày đặc (vịnh Cam Ranh, vịnh Nha Phu, đầm Nha Phu) và cửa sông thể hiện trạng thái phú dưỡng và siêu phú dưỡng thường trực (Chl-a $> 3,0\ \mu\text{g/L}$). Trong khi đó, vùng nước mở ngoài khơi duy trì trạng thái nghèo đến trung dưỡng (diện tích vùng nghèo dinh dưỡng $< 0,2%$ được tích hợp trong nhóm trung dưỡng chiếm hơn 80% tổng diện tích vùng nghiên cứu).

Implications đa chiều

  • Về lý thuyết: Cung cấp bằng chứng thực nghiệm khẳng định tính khả thi của việc mô hình hóa các quá trình sinh địa hóa đại dương ở quy mô khu vực nhiệt đới; hoàn thiện cơ sở khoa học cho chu trình carbon và đánh giá năng suất sơ cấp đại dương tại thềm lục địa Biển Đông.
  • Về phương pháp luận: Thiết lập một quy trình công nghệ hoàn chỉnh, tự động hóa từ xử lý ảnh viễn thám cấp độ thô (Level-1B) đến sản phẩm chuyên đề cấp cao (Level-4 Spatio-temporally reconstructed products), có thể chuyển giao và mở rộng áp dụng cho toàn bộ 28 tỉnh thành ven biển Việt Nam.
  • Về ứng dụng thực tiễn: Cung cấp công cụ cảnh báo sớm hiện tượng phú dưỡng và nguy cơ bùng phát tảo độc nở hoa cho các vùng nuôi trồng thủy hải sản trọng điểm tại Nam Trung Bộ, giúp ngư dân và các cơ sở nuôi biển chủ động giảm thiểu rủi ro kinh tế.
  • Về chính sách và quản trị biển: Đóng góp cơ sở dữ liệu số hóa địa không gian chuẩn xác phục vụ Bộ Tài nguyên và Môi trường, Ủy ban Nhân dân các tỉnh Khánh Hòa, Ninh Thuận trong việc quy hoạch không gian biển (Marine Spatial Planning), kiểm soát nguồn thải lục địa theo tinh thần Luật Bảo vệ Môi trường 2020 và Nghị quyết 36-NQ/TW.

Limitations và Future Research

Luận án thừa nhận một cách khách quan các giới hạn khoa học nội tại:

  1. Hiệu ứng đáy nông và rạn san hô: Tại các dải nước cực nông ven bờ và quanh các đảo có rạn san hô (độ sâu $< 3\text{ m}$), bức xạ phản xạ từ đáy cát và san hô vẫn gây nhiễu nhất định lên việc tính toán $p_w(\lambda)$, đòi hỏi các mô hình bóc tách đáy sâu chuyên biệt trong tương lai.
  2. Mật độ mẫu thực địa kiểm định: Do điều kiện kinh phí và khó khăn trong công tác đi biển mùa sóng gió, số lượng điểm đo thực địa đồng thời với thời gian vệ tinh bay qua còn giới hạn ở các tháng đại diện, chưa bao phủ toàn bộ 12 tháng liên tục trong nhiều năm.
  3. Hiện tượng che phủ mây cực đoan: Mặc dù thuật toán DINEOF phục hồi xuất sắc các khoảng trống dữ liệu mây cục bộ, nhưng đối với các giai đoạn thời tiết bão nhiệt đới kéo dài nhiều ngày với mây bao phủ $100%$ toàn bộ khung cảnh quan sát, mô hình toán học gặp hạn chế trong việc dự báo xu thế ngắn hạn.

Chương trình nghiên cứu tương lai (Future Research Agenda):

  • Hướng 1: Tích hợp dữ liệu siêu phổ (Hyperspectral Ocean Color) từ các thế hệ vệ tinh mới như NASA PACE (Ocean Color Instrument - OCI) để phân giải chi tiết các nhóm loài thực vật phù du (phytoplankton functional types - PFTs) gây hại.
  • Hướng 2: Ứng dụng các kiến trúc mạng học sâu tiên tiến (như Convolutional Auto-Encoders DINACE, Physics-Informed Neural Networks - PINN) để nâng cao độ chính xác tái cấu trúc không gian - thời gian.
  • Hướng 3: Mở rộng nghiên cứu ước tính đồng thời các biến số sinh địa hóa khác bao gồm tổng carbon hữu cơ hòa tan (DOC), carbon hữu cơ dạng hạt (POC), vật chất lơ lửng (TSS) và nhiệt độ bề mặt biển (SST) để mô hình hóa toàn diện chu trình năng suất sơ cấp ven biển Việt Nam.

Tác động và ảnh hưởng

Nghiên cứu tạo ra tác động học thuật và thực tiễn sâu rộng:

  • Ảnh hưởng học thuật (Academic Impact): Đặt nền móng tiên phong cho việc khai thác dữ liệu vệ tinh Sentinel-3 OLCI trong hải dương học viễn thám tại Việt Nam; tạo tiền đề cho hàng loạt công bố quốc tế uy tín trong danh mục ISI/Scopus và dự kiến thu hút trích dẫn cao trong cộng đồng nghiên cứu màu đại dương Đông Nam Á.
  • Chuyển đổi ngành kinh tế biển: Giúp chuyển đổi phương thức giám sát môi trường nuôi biển truyền thống (lấy mẫu thủ công, tốn kém hàng tỷ đồng mỗi năm) sang phương thức số hóa tự động tầm xa diện rộng, giảm $80%$ chi phí quan trắc định kỳ mà vẫn nâng cao độ phân giải thời gian từ hàng tháng lên hàng ngày.
  • Hỗ trợ chính sách quốc gia: Cung cấp công cụ kiểm chứng khoa học phục vụ quản lý tổng hợp đới bờ (ICM), bảo vệ các khu bảo tồn biển (như Vịnh Nha Trang, Núi Chúa - Ninh Thuận), và cung cấp dữ liệu pháp lý - kỹ thuật khẳng định chủ quyền, quyền chủ quyền và quyền tài phán của Việt Nam trên vùng đặc quyền kinh tế Biển Đông.
  • Đóng góp quốc tế: Kết quả nghiên cứu đóng góp bộ dữ liệu thực địa nhiệt đới quý giá vào mạng lưới màu đại dương toàn cầu IOCCG (International Ocean Colour Coordinating Group), hỗ trợ hiệu chỉnh chéo các thuật toán vệ tinh quốc tế tại khu vực biển nhiệt đới Đông Nam Á.

Đối tượng hưởng lợi

  • Nghiên cứu sinh và học viên cao học: Tiếp cận quy trình phương pháp luận chuẩn mực, các script xử lý dữ liệu viễn thám mở (SNAP/Python/DINEOF), và các hướng nghiên cứu mới về viễn thám môi trường biển Case-2.
  • Các nhà khoa học và giảng viên đại học: Có được bộ tư liệu tham khảo học thuật chuyên sâu về quang học hải dương, phục vụ giảng dạy các chuyên ngành Bản đồ, Viễn thám, GIS, Hải dương học và Khoa học Môi trường.
  • Cơ quan quản lý nhà nước (Bộ Tài nguyên và Môi trường, Sở TN&MT các tỉnh ven biển): Sở hữu giải pháp công nghệ hiện đại để nâng cấp Hệ thống Quan trắc Môi trường Biển Quốc gia, thiết lập cơ chế giám sát và cảnh báo ô nhiễm môi trường nước biển theo thời gian thực.
  • Doanh nghiệp và cộng đồng nuôi trồng thủy sản: Được hưởng lợi từ các bản đồ cảnh báo chất lượng nước và nguy cơ phú dưỡng/tảo độc, giúp bảo vệ an toàn sản xuất lồng bè quy mô hàng nghìn tỷ đồng tại khu vực Nam Trung Bộ.

Câu hỏi chuyên sâu

  1. Đóng góp lý thuyết độc đáo nhất của luận án là gì và đã mở rộng lý thuyết nào? Trả lời: Đóng góp lý thuyết độc đáo nhất là việc mở rộng và kiểm chứng thành công lý thuyết phân loại quang học nước biển Case-2 của Morel & Prieur (1977) kết hợp với mô hình mạng nơ-ron C2RCC cho vùng biển ven bờ nhiệt đới gió mùa Việt Nam. Luận án đã chứng minh được mối tương quan phi tuyến giữa các đặc tính quang học nội tại (IOP) và biểu kiến (AOP) trong điều kiện nồng độ CDOM và trầm tích lơ lửng biến động mạnh, bẻ gãy giả định tuyến tính đơn giản của các mô hình đại dương mở truyền thống.

  2. Đột phá phương pháp luận của luận án khi so sánh với ít nhất 2 nghiên cứu quốc tế trước đây? Trả lời: So với nghiên cứu của Gohin et al. (2002) tại vùng biển châu Âu (chỉ dùng đơn vệ tinh SeaWiFS độ phân giải 1 km) và nghiên cứu của Son & Wang (2012) tại Vịnh Chesapeake (dùng MODIS với thuật toán khí quyển tiêu chuẩn bị bão hòa kênh đỏ), đột phá của luận án nằm ở giải pháp tích hợp chéo cặp vệ tinh Sentinel-3A/3B (độ phân giải cao 300 m) với thuật toán hiệu chỉnh khí quyển C2RCC và mô hình toán học tái cấu trúc không gian DINEOF. Giải pháp này giải quyết triệt để bài toán mất dữ liệu do mây nhiệt đới – điều mà các nghiên cứu trước đây chưa thực hiện được đồng bộ.

  3. Phát hiện nào gây bất ngờ nhất và bằng chứng số liệu chứng minh? Trả lời: Phát hiện bất ngờ nhất là tại vùng nước vịnh kín ven bờ Nam Trung Bộ, thuật toán tiêu chuẩn OC4ME của ESA/NASA đánh giá vượt mức (overestimate) nồng độ Chl-a từ $35%$ đến hơn $58%$ so với giá trị đo thực địa in-situ. Ngược lại, thuật toán OC5v7 khi kết hợp với C2RCC cho kết quả bám sát thực tế với sai số $RMSE < 0,28\ \mu\text{g/L}$ ($p < 0,01$), chứng minh rằng việc áp dụng trực tiếp các sản phẩm Chl-a tiêu chuẩn toàn cầu chưa qua hiệu chuẩn khu vực vào vùng biển Việt Nam sẽ dẫn đến những nhận định sai lệch nghiêm trọng về hiện trạng phú dưỡng.

  4. Luận án có cung cấp giao thức tái lặp (replication protocol) đầy đủ không? Trả lời: Luận án xây dựng quy trình thực nghiệm minh bạch, chi tiết và có khả năng tái lặp hoàn toàn: từ cấu hình thông số kỹ thuật bộ cảm OLCI, các bước thiết lập tham số trong SNAP/Acolite, danh mục mask phân loại (land, cloud, sun glint, whitecaps), ma trận trọng số trong mô hình C2RCC, đến thuật toán lập trình DINEOF trong môi trường Python. Mọi nhà nghiên cứu đều có thể sử dụng quy trình này để tái lập kết quả trên các vùng biển khác.

  5. Chương trình nghiên cứu 10 năm tiếp theo được vạch ra như thế nào? Trả lời: Luận án định hình lộ trình nghiên cứu 10 năm (2023–2033) tập trung vào 3 trụ cột: (1) Xây dựng Hệ sinh thái Dữ liệu lớn Viễn thám Môi trường Biển Việt Nam (Vietnam Ocean Big Data Platform) tích hợp đa thế hệ vệ tinh Sentinel-3C/3D, PACE OCI và VNREDSat; (2) Phát triển các mô hình học sâu Physics-Informed AI để dự báo trước 72 giờ các đợt bùng phát tảo nở hoa độc hại; (3) Định lượng hóa chu trình hấp thụ carbon xanh (blue carbon) và năng suất sơ cấp đại dương toàn vùng biển Đông Việt Nam phục vụ mục tiêu Net-Zero 2050.

Kết luận

Công trình luận án tiến sĩ của tác giả Nguyễn An Bình đại diện cho một bước tiến quan trọng trong ứng dụng công nghệ địa không gian và viễn thám màu đại dương tại Việt Nam với 5 đóng góp cốt lõi:

  1. Xác lập cơ sở khoa học và quy trình công nghệ chuẩn hóa trong việc bóc tách thông số chất lượng nước biển từ dữ liệu vệ tinh thế hệ mới Sentinel-3 OLCI độ phân giải cao 300 m.
  2. Chứng minh bằng thực nghiệm và kiểm định chéo chuẩn xác sự vượt trội của mô hình hiệu chỉnh khí quyển C2RCC và thuật toán quang sinh học OC5v7 cho vùng biển ven bờ Case-2 nhiệt đới.
  3. Đột phá giải quyết triệt để điểm nghẽn mây che phủ thông qua mô hình tái cấu trúc dữ liệu không gian DINEOF kết hợp cặp vệ tinh Sentinel-3A và Sentinel-3B, hiện thực hóa khả năng giám sát chất lượng nước hàng ngày.
  4. Cung cấp bộ bản đồ chuyên đề định lượng về biến động Chl-a và phân vùng mức độ phú dưỡng tại dải ven biển Nam Trung Bộ, phát hiện các điểm nóng ô nhiễm cục bộ tại các vịnh nuôi trồng thủy sản trọng điểm.
  5. Mở ra hướng nghiên cứu liên ngành mới, liên kết giữa viễn thám quang sinh học hải dương, trí tuệ nhân tạo địa không gian và quản trị môi trường biển bền vững.

Di sản khoa học của luận án không chỉ đóng góp trực tiếp vào kho tàng tri thức hải dương học nước nhà, mà còn là công cụ đắc lực hỗ trợ các cơ quan quản lý trong công cuộc bảo vệ tài nguyên môi trường biển, phục vụ phát triển kinh tế biển bền vững và khẳng định vững chắc chủ quyền biển đảo thiêng liêng của Tổ quốc.