Tổng quan về luận án

Sự biến đổi khí hậu toàn cầu đang đặt tài nguyên nước trước nguy cơ suy thoái và cạn kiệt nghiêm trọng, đòi hỏi các chiến lược quản lý, khai thác và sử dụng bền vững. Trong chu trình thủy văn, bốc thoát hơi nước (Evapotranspiration - ET) là thành phần quan trọng thứ hai chỉ sau lượng mưa, đóng vai trò then chốt trong cân bằng nước và trao đổi năng lượng giữa bề mặt đất và khí quyển. Nông nghiệp là ngành tiêu thụ nguồn nước ngọt lớn nhất, do đó việc ước tính chính xác nhu cầu nước của cây trồng thông qua chỉ số bốc thoát hơi nước thực tế ($ET_a$) là tiền đề cốt lõi để nâng cao hiệu quả tưới tiêu và ứng phó hạn hán.

Tại Việt Nam, mạng lưới quan trắc khí tượng thủy văn mặt đất còn rất thưa thớt, trung bình mỗi tỉnh chỉ có trên dưới 10 trạm đo đạc thủ công với chi phí vận hành cao nhưng độ bao phủ không gian hạn chế. Đặc biệt tại khu vực Tây Bắc — nơi có địa hình núi cao chia cắt hiểm trở, chênh cao lớn với dãy Hoàng Liên Sơn dài 180 km, rộng 30 km, dãy Sông Mã dài 500 km, lưu vực sông Đà phức tạp và độ che phủ rừng đạt 44,7% — các trạm mặt đất hoàn toàn không thể phản ánh tính dị thường của các tiểu vùng khí hậu. Dữ liệu đo đạc trực tiếp (như ống Piche hay Thủy tiêu kế - Lysimeter) chỉ mang tính đại diện cục bộ, dẫn đến khoảng trống nghiên cứu lớn (research gap) trong việc giám sát $ET_a$ liên tục trên diện rộng với độ phân giải cao.

Luận án "Nghiên cứu xây dựng mô hình giám sát sự bốc - thoát hơi nước của lớp phủ khu vực Tây Bắc Việt Nam từ dữ liệu ảnh vệ tinh" (Chuyên ngành: Kỹ thuật Trắc địa - Bản đồ / Viễn thám và GIS, Trường Đại học Mỏ - Địa chất) giải quyết triệt để bài toán này thông qua hai câu hỏi nghiên cứu và giả thuyết khoa học:

  • Câu hỏi nghiên cứu 1 (RQ1): Làm thế nào để chiết xuất chính xác các tham số vật lý bề mặt phục vụ ước tính $ET_a$ tại địa hình phức tạp Tây Bắc hoàn toàn từ tư liệu viễn thám quang học/nhiệt mà không phụ thuộc vào chuỗi đo khí tượng mặt đất thời gian thực?
  • Câu hỏi nghiên cứu 2 (RQ2): Cặp hệ số tuyến tính thực nghiệm ($a, b$) trong mô hình bức xạ Priestley - Taylor được tối ưu hóa như thế nào để phản ánh trung thực đặc trưng lớp phủ đa dạng tại tỉnh Hòa Bình?
  • Giả thuyết nghiên cứu 1 (H1): Các tham số bức xạ ròng trung bình ngày ($R_{nd}$), nhiệt ẩn hóa hơi ($\lambda$), hằng số Psychrometric ($\gamma$) và độ dốc áp suất hơi bão hòa ($\Delta$) tính toán trực tiếp từ ảnh Landsat 8 kết hợp Mô hình số độ cao (DEM) có độ chính xác tương đương và có thể thay thế hoàn toàn chuỗi số liệu khí tượng trạm bề mặt.
  • Giả thuyết nghiên cứu 2 (H2): Mô hình lai ghép giữa thuật toán Cân bằng năng lượng bề mặt (SEBAL) và mô hình Priestley - Taylor với hệ số ($a, b$) hiệu chỉnh thực nghiệm trên nền tảng Google Earth Engine cho phép định lượng hóa $ET_a$ tức thời và ngày với sai số tuyệt đối dưới 10% so với trạm thực nghiệm.

Phạm vi nghiên cứu tập trung thực nghiệm tại tỉnh Hòa Bình (Tây Bắc Việt Nam), khai thác chuỗi dữ liệu khí tượng giai đoạn 2015–2021 và dữ liệu ảnh vệ tinh Landsat 8 tại các thời điểm điển hình (ngày 01/7/2015, 04/6/2017 và 18/8/2021). Ý nghĩa khoa học đột phá của nghiên cứu là thiết lập cơ sở tự động hóa giám sát tài nguyên nước địa không gian, cung cấp công cụ chuyển giao cho quản lý tưới tiêu, dự báo hạn hán và phòng cháy rừng vùng cao.

Literature Review và Positioning

Nghiên cứu về bốc thoát hơi nước bề mặt đã trải qua sự tiến hóa mạnh mẽ từ các phương pháp đo cơ học, mô hình dựa trên trạm khí tượng đến các thuật toán viễn thám vệ tinh đa phổ:

flowchart LR
    A["Mô hình vật lý cổ điển<br/>Penman 1948<br/>Monteith 1965"] --> B["Tiêu chuẩn hóa quốc tế<br/>FAO-56 Penman-Monteith<br/>Allen et al. 1998"]
    B --> C["Mô hình đơn giản hóa bức xạ<br/>Makkink 1957<br/>Priestley-Taylor 1972"]
    C --> D["Thuật toán cân bằng năng lượng viễn thám<br/>SEBAL, SEBS, METRIC, S-SEBI<br/>Bastiaanssen 1998, Su 2002"]
    D --> E["Mô hình tích hợp viễn thám lai<br/>Luận án: SEBAL + DEM + Priestley-Taylor<br/>Google Earth Engine"]

Tổng hợp các dòng nghiên cứu quốc tế và tranh luận học thuật

  1. Dòng mô hình kết hợp truyền thống: Dựa trên phương trình Penman (1948) và Monteith (1965), Tổ chức Lương thực và Nông nghiệp Liên Hiệp Quốc (FAO) đã chuẩn hóa phương pháp FAO-56 Penman - Monteith (Allen et al., 1998) làm tiêu chuẩn vàng để tính bốc thoát hơi nước tham chiếu ($ET_0$). Mặc dù có độ tin cậy vượt trội, phương pháp này đòi hỏi tập dữ liệu khí hậu dày đặc (nhiệt độ, ẩm độ, bức xạ, tốc độ gió ở độ cao 2m), khiến việc ngoại suy không gian tại các vùng núi cao gặp sai số nội suy rất lớn.
  2. Dòng mô hình dựa trên bức xạ (Radiation-based models): Nhằm giảm thiểu tham số đầu vào, Makkink (1957) và Priestley - Taylor (1972) đã đề xuất các mô hình tính $ET_a$ dựa trên bức xạ mặt trời ($R_s$) hoặc bức xạ ròng ($R_n$). Priestley - Taylor giả định rằng trong điều kiện bề mặt bão hòa ẩm, dòng nhiệt đối lưu khí quyển đạt trạng thái cân bằng và lượng bốc thoát hơi nước chủ yếu do bức xạ chi phối qua phương trình: $$ET_a = a \frac{\Delta}{\Delta + \gamma} \frac{R_n}{\lambda} + b$$ Tuy nhiên, các cuộc tranh luận thực nghiệm quốc tế cho thấy cặp hệ số ($a, b$) không phải là hằng số toàn cầu. Priestley - Taylor (1972) ban đầu đề xuất $a = 0,61; b = 0,12$; kiểm định tại Thụy Sĩ (1984) cho kết quả $a = 0,90; b = 0$; trong khi khảo sát tại Đài Loan (2005) ghi nhận $a = 1,00; b = 0$. Xu và Singh (2001) khi so sánh 4 mô hình bức xạ với chậu đo bốc hơi Pan đã kết luận rằng mô hình Makkink và Priestley - Taylor đạt độ tin cậy cao nhất nhưng bắt buộc phải hiệu chuẩn hệ số theo từng vùng khí hậu vi mô.
  3. Dòng thuật toán cân bằng năng lượng bề mặt từ viễn thám: Sự ra đời của các mô hình SEBAL (Bastiaanssen et al., 1998), S-SEBI (Roerink et al., 2000), SEBS (Su, 2002) và METRIC (Allen et al., 2007) đã mở ra kỷ nguyên mới. SEBAL giải phương trình cân bằng năng lượng: $$\lambda LE = R_n - G - H$$ thông qua việc tự động hiệu chuẩn nội bộ dựa trên việc chọn điểm ảnh "nóng" (hot pixel - nơi $ET = 0$) và điểm ảnh "lạnh" (cold pixel - nơi $ET$ đạt cực đại). Tuy nhiên, các công trình quốc tế như Santos et al. (2005) tại Tây Ban Nha, Gowda et al. (2005) và Tasumi et al. (2005) tại vùng bán khô hạn Hoa Kỳ đã chỉ ra nhược điểm chí mạng của SEBAL và METRIC: việc xác định pixel nóng/lạnh mang tính chủ quan cao và mô hình giả định bề mặt tương đối bằng phẳng, gây sai số nghiêm trọng khi áp dụng cho địa hình núi chia cắt có hiệu ứng che khuất và biến thiên nhiệt độ theo độ cao.

So sánh với các nghiên cứu trong nước và định vị khoảng trống

Tại Việt Nam, các nghiên cứu viễn thám trước đây chủ yếu sử dụng ảnh MODIS độ phân giải thấp (250m - 1km) hoặc chỉ áp dụng đơn lẻ trên địa hình đồng bằng:

  • Lê Minh, Lương Chính Kế, Trần Ngọc Tưởng et al. (2013, 2014) sử dụng ảnh MODIS tính bức xạ ròng và $ET$ theo Makkink/Priestley-Taylor với sai số từ 6% đến 10% ở miền Bắc. Trần Thị Tâm et al. (2016) áp dụng thuật toán MODIS 16 cho Bắc Trung Bộ với sai lệch so với Penman-Monteith từ -53,49 đến 230,10 mm/năm.
  • Nguyễn Văn Hùng (2017) ứng dụng S-SEBI và Landsat 7 cho lưu vực Sông Cầu với sai số 11% - 18%, nhưng chỉ thử nghiệm các cặp hệ số giả định mà chưa xây dựng cơ sở thực nghiệm xác định hệ số ($a, b$) tối ưu.
  • Phạm Thị Thanh Ngà et al. (2020) dùng Himawari-8 ước tính bức xạ mặt trời tại Hà Nội và Đà Nẵng với sai số 1,41 - 2,14 $\text{MJ/m}^2/\text{ngày}$, tương đương sai số ET khoảng 0,6 - 0,9 mm/ngày.
  • Trần Thị Hồng Ngọc và Mark Honti (2017) khảo sát vùng Tứ giác Long Xuyên (An Giang) giai đoạn 2010–2015 cho thấy ET dao động 3,73 - 4,83 mm/ngày (1.761 mm/năm) nhưng hoàn toàn dựa vào trạm đo mặt đất.

Khoảng trống nghiên cứu được xác định: Chưa có bất kỳ công trình nào tích hợp mô hình viễn thám SEBAL để chiết xuất các tham số vi khí hậu phân giải cao (30m) từ Landsat 8 kết hợp mô hình số độ cao DEM và ghép nối trực tiếp vào mô hình năng lượng Priestley - Taylor với bộ hệ số ($a, b$) được hiệu chuẩn hóa riêng biệt cho vùng núi Tây Bắc Việt Nam.

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

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

Luận án đã mở rộng và làm sâu sắc thêm lý thuyết Cân bằng năng lượng bề mặt (Surface Energy Balance) của Bastiaanssen et al. (1998) và lý thuyết Bốc thoát hơi nước bức xạ của Priestley - Taylor (1972) trong điều kiện địa hình phức tạp:

Thành phần lý thuyết Lý thuyết nguyên bản Mở rộng của Luận án Bằng chứng thực nghiệm
Cơ chế xác định tham số khí tượng Phụ thuộc trạm đo mặt đất hoặc nội suy không gian thô. Tham số hóa phân tán 100% từ kênh phổ quang học và nhiệt Landsat 8 kết hợp DEM. Chiết xuất thành công $R_{nd}, T_s, \lambda, \gamma, \Delta$ ở độ phân giải 30m.
Mô hình hóa độ cao địa hình Coi bề mặt đẳng thế hoặc hiệu chỉnh tuyến tính đơn giản. Tích hợp DEM hiệu chỉnh áp suất khí quyển $P$, nhiệt độ bề mặt $T_s$ và hằng số $\gamma$. Triệt tiêu sai số do hiệu ứng sườn phơi và độ cao địa hình Tây Bắc.
Hệ số kinh nghiệm ($a, b$) Priestley - Taylor Cố định toàn cầu ($a=0,61; b=0,12$ hoặc $a=1,26$). Tái cấu trúc cặp hệ số tuyến tính động theo chuỗi thực nghiệm khí tượng 2015–2021. Đảm bảo độ tương quan cao và sai số tính toán dưới 10% tại Hòa Bình.

Mô hình lý thuyết đề xuất xác lập mối liên hệ trực tiếp giữa bức xạ ròng trung bình ngày ($R_{nd}$) hấp thụ bởi thảm thực vật với lượng bốc thoát hơi nước thực tế $ET_a$ mà không cần giải bài toán lặp phức tạp để tìm dòng nhiệt cảm ứng ($H$) vốn chứa đựng nhiều rủi ro sai số chủ quan khi chọn pixel nóng/lạnh.

graph TD
    subgraph INPUT["Dữ liệu đầu vào vệ tinh & Địa hình"]
        L8["Ảnh Landsat 8 (Quang học & Nhiệt TIRS)"]
        DEM["Mô hình số độ cao DEM (SRTM/ASTER)"]
        MET["Chuỗi số liệu khí tượng thủy văn 2015-2021"]
    end

    subgraph SEBAL_PARAM["Chiết xuất tham số vật lý viễn thám (SEBAL Algorithm)"]
        ALB["Suất phân sai bề mặt (α)"]
        TS["Nhiệt độ bề mặt đất (Ts)"]
        VI["Chỉ số thực vật (NDVI, LAI, SAVI)"]
        EMISS["Độ phát xạ bề mặt (εo)"]
        RN["Bức xạ ròng mặt trời (Rn & Rnd)"]
    end

    subgraph THERMO_PARAM["Tham số nhiệt động học bề mặt"]
        LAMBDA["Nhiệt ẩn hóa hơi (λ)"]
        GAMMA["Hằng số Psychrometric hiệu chỉnh DEM (γ)"]
        DELTA["Độ dốc áp suất hơi bão hòa (Δ)"]
        COEFF["Hiệu chuẩn thực nghiệm hệ số (a, b)"]
    end

    subgraph OUTPUT["Mô hình tích hợp Priestley - Taylor trên GEE"]
        ETA["Ước tính & Lập bản đồ phân bố bốc thoát hơi nước thực tế (ETa)"]
    end

    L8 --> ALB & TS & VI & EMISS
    ALB & TS & VI & EMISS --> RN
    DEM --> GAMMA
    TS --> DELTA & LAMBDA
    MET --> COEFF
    RN & LAMBDA & GAMMA & DELTA & COEFF --> ETA

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

Khung phân tích của luận án tích hợp đồng thời 3 lý thuyết nền tảng:

  1. Lý thuyết truyền xạ quang phổ điện từ: Tính toán bức xạ ròng $R_n$ dựa trên phương trình cân bằng bức xạ sóng ngắn và sóng dài: $$R_n = R_{S\downarrow} - \alpha R_{S\downarrow} + R_{L\downarrow} - R_{L\uparrow} - (1 - \varepsilon_o)R_{L\downarrow}$$ trong đó suất phân sai bề mặt ($\alpha$) được tích hợp từ các kênh phổ quang học, nhiệt độ bề mặt ($T_s$) giải đoán từ kênh hồng ngoại nhiệt (TIRS), và độ phát xạ bề mặt ($\varepsilon_o$) được xác định thông qua hàm chỉ số thực vật NDVI và SAVI.
  2. Lý thuyết nhiệt động học khí quyển: Xác định độ dốc đường cong áp suất hơi bão hòa $\Delta = \frac{4098 \cdot [0,6108 \exp(\frac{17,27 T_s}{T_s + 237,3})]}{(T_s + 237,3)^2}$, nhiệt hóa hơi $\lambda = 2,501 - (0,002361 \cdot T_s)$, và hằng số $\gamma = 0,00163 \frac{P}{\lambda}$ với áp suất khí quyển $P$ được điều chỉnh theo độ cao DEM.
  3. Lý thuyết ước lượng thống kê thực nghiệm: Sử dụng chuỗi số liệu khí tượng đa năm tại Hòa Bình để giải bài toán tối ưu hóa tham số tuyến tính ($a, b$).

Điều kiện biên (Boundary conditions): Mô hình áp dụng tối ưu cho các cảnh ảnh Landsat 8 có độ che phủ mây dưới 10% trong mùa sinh trưởng thảm phủ vùng nhiệt đới gió mùa, trên nền địa hình có độ cao biến thiên từ 100m đến trên 2.000m.

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

Thiết kế nghiên cứu

  • Triết lý nghiên cứu (Research Philosophy): Nghiên cứu theo chủ nghĩa thực chứng định lượng (Positivism) kết hợp tiếp cận thực nghiệm mô hình hóa địa không gian (Spatial empirical modeling). Mọi quy luật trao đổi nhiệt ẩm bề mặt đều được định lượng hóa thông qua các phương trình vật lý truyền xạ và kiểm chứng bằng dữ liệu trạm quan trắc mặt đất khách quan.
  • Thiết kế tích hợp đa tầng (Multi-level design):
    • Tầng 1 (Vật lý viễn thám): Xử lý mức pixel (30m x 30m) để chiết xuất các đặc trưng phản xạ phổ, albedo, chỉ số thực vật ($NDVI, SAVI, LAI$) và nhiệt độ bề mặt ($T_s$).
    • Tầng 2 (Khí động và bức xạ): Tính toán các trường thông lượng năng lượng bức xạ ròng ($R_n, R_{nd}$) và các tham số khí vi mô ($\lambda, \gamma, \Delta$) tích hợp DEM.
    • Tầng 3 (Tối ưu hóa mô hình): Hiệu chuẩn thống kê chuỗi thời gian trắc nghiệm hệ số ($a, b$) và tổng hợp bản đồ $ET_a$ trên toàn bộ không gian tỉnh Hòa Bình.
flowchart TD
    subgraph S1["Giai đoạn 1: Tiền xử lý dữ liệu viễn thám & DEM"]
        D1["Ảnh Landsat 8 C1/C2 (30m)"] --> P1["Hiệu chỉnh bức xạ & Khí quyển FLAASH"]
        D2["DEM SRTM 30m"] --> P2["Tính toán độ dốc, hướng sườn & Áp suất P"]
    end

    subgraph S2["Giai đoạn 2: Chiết xuất tham số bề mặt (SEBAL)"]
        P1 --> F1["Suất phân sai α & Chỉ số NDVI, SAVI"]
        P1 --> F2["Nhiệt độ bề mặt Ts (Kênh nhiệt 10/11)"]
        F1 & F2 & P2 --> F3["Bức xạ ròng tức thời Rn & Ngày Rnd"]
        F2 & P2 --> F4["Nhiệt ẩn λ, Hằng số γ, Độ dốc Δ"]
    end

    subgraph S3["Giai đoạn 3: Hiệu chuẩn hệ số thực nghiệm (a, b)"]
        MET["Số liệu trạm khí tượng Hòa Bình 2015-2021"] --> CAL["Phân tích hồi quy & Thực nghiệm mô hình"]
        CAL --> COEFF["Bộ hệ số tối ưu (a, b) cho Tây Bắc"]
    end

    subgraph S4["Giai đoạn 4: Tự động hóa & Kiểm định trên GEE"]
        F3 & F4 & COEFF --> GEE["Script xử lý tự động Google Earth Engine"]
        GEE --> MAP["Bản đồ phân bố không gian ETa (mm/ngày)"]
        MAP --> VAL["So sánh kiểm chứng với số liệu đo trạm khí tượng thủy văn"]
    end

Quy trình nghiên cứu rigorous

  1. Thu thập dữ liệu thứ cấp: Thu thập dữ liệu khí tượng đo đạc hàng ngày (nhiệt độ, độ ẩm, số giờ nắng, tốc độ gió, lượng bốc hơi Piche) tại các trạm khí tượng thủy văn tỉnh Hòa Bình giai đoạn 2015–2021. Dữ liệu viễn thám bao gồm các cảnh ảnh Landsat 8 (Path/Row: 127/045 và 127/046) thời điểm 01/7/2015, 04/6/2017 và 18/8/2021 do Cục Khảo sát Địa chất Hoa Kỳ (USGS) cung cấp.
  2. Quy trình tính toán bức xạ ròng $R_n$ từ ảnh Landsat 8:
    • Bức xạ sóng ngắn tới ($R_{S\downarrow}$): Tính toán từ hằng số mặt trời ($G_{sc} = 1367 \text{ W/m}^2$), góc cao mặt trời ($\theta$), khoảng cách tương đối Trái đất - Mặt trời ($d_r$), và độ truyền dẫn khí quyển sóng ngắn ($\tau_{sw}$).
    • Suất phân sai bề mặt ($\alpha$): Tính từ tổ hợp phản xạ dải phổ nhìn thấy và cận hồng ngoại có trọng số theo phương pháp Tasumi et al.
    • Bức xạ sóng dài tới và đi ($R_{L\downarrow}, R_{L\uparrow}$): Ứng dụng định luật Stefan - Boltzmann với độ phát xạ khí quyển ($\varepsilon_a$) và độ phát xạ bề mặt ($\varepsilon_o$).
  3. Triangulation và kiểm định:
    • Data triangulation: Đối chiếu số liệu bức xạ tính từ viễn thám với chuẩn FAO-56 và số liệu trạm quan trắc mặt đất.
    • Methodological triangulation: Thực nghiệm so sánh 4 mô hình năng lượng bức xạ độc lập: Priestley - Taylor, Makkink, Hargreaves & Samani, và Turc.

Data và phân tích

  • Đặc trưng khu vực thực nghiệm: Tỉnh Hòa Bình đại diện cho vùng cửa ngõ Tây Bắc, độ che phủ rừng cao (rừng phòng hộ, đặc dụng, sản xuất xen kẽ cây nông nghiệp), địa hình đồi núi dốc phân cắt mạnh tạo nên các lòng chảo và thung lũng có vi khí hậu nhiệt đới gió mùa ẩm khác biệt rõ rệt.
  • Nền tảng công nghệ xử lý: Toàn bộ thuật toán tích hợp SEBAL - Priestley - Taylor được lập trình tự động hóa bằng ngôn ngữ JavaScript trên nền tảng điện toán đám mây Google Earth Engine (GEE) kết hợp phần mềm chuyên dụng ArcGIS và ENVI, giúp xử lý dữ liệu không gian tức thời ở quy mô toàn vùng mà không đòi hỏi hạ tầng tính toán cục bộ đắt đỏ.

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

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

Nghiên cứu mang lại 4 phát hiện khoa học mang tính đột phá với các bằng chứng thực nghiệm rõ ràng:

  1. Khả năng thay thế hoàn toàn dữ liệu trạm mặt đất bằng tham số viễn thám kết hợp DEM: Các tham số then chốt gồm $R_{nd}, T_s, \lambda, \gamma, \Delta$ tính trực tiếp từ ảnh vệ tinh Landsat 8 và mô hình số độ cao DEM phản ánh chính xác các biến động nhiệt ẩm bề mặt theo độ cao và lớp phủ. Khi kiểm định bức xạ ròng $R_n$ tính từ viễn thám so với phương pháp chuẩn FAO-56, sai số tuyệt đối trung bình chỉ đạt mức xấp xỉ 6%, tương thích với các công bố quốc tế của Lương Chính Kế et al. (2014) và tốt hơn mức sai số 10% của ảnh MODIS (Lê Minh et al., 2013).

  2. Xác lập bộ hệ số tuyến tính thực nghiệm ($a, b$) đặc thù cho Tây Bắc: Qua phân tích chuỗi dữ liệu khí tượng 2015–2021 tại tỉnh Hòa Bình, luận án đã xác định được bộ hệ số thực nghiệm tối ưu cho mô hình Priestley - Taylor phù hợp với địa hình phân cắt mạnh và lớp phủ rừng hỗn giao. Kết quả này khắc phục hoàn toàn hiện tượng thiên lệch sai số khi áp dụng hệ số chuẩn quốc tế ($a = 1,26$ hoặc $a = 0,61; b = 0,12$) vốn chỉ phù hợp cho đồng bằng hoặc cây trồng thuần loài.

  3. Độ chính xác vượt trội của mô hình lai SEBAL - Priestley - Taylor: Ước tính $ET_a$ tại 3 mốc thời gian thực nghiệm (01/7/2015, 04/6/2017 và 18/8/2021) cho thấy sự tương thích rất cao với giá trị đo thực tế tại các trạm khí tượng thủy văn. Sai số tuyệt đối trung bình toàn vùng duy trì ở mức dưới 10%, vượt trội hơn đáng kể so với mô hình S-SEBI trên Landsat 7 của Nguyễn Văn Hùng (2017) tại Sông Cầu (sai số 11% - 18%).

  4. Phản ánh trung thực quy luật phân hóa không gian của bốc thoát hơi nước: Dữ liệu ảnh viễn thám phân giải 30m bóc tách rõ nét sự tương phản $ET_a$:

    • Khu vực rừng đặc dụng, rừng tự nhiên trên núi cao và mặt nước các hồ chứa (như lòng hồ thủy điện Hòa Bình) có giá trị $ET_a$ đạt cực đại (dao động từ 4,5 đến 5,8 mm/ngày trong mùa hè) do sinh khối lớn, tán lá rậm rạp và nguồn cung cấp nước dồi dào.
    • Ngược lại, các khu vực đất trống, đồi trọc, bề mặt đô thị hóa hoặc đất nông nghiệp sau thu hoạch có $ET_a$ thấp nhất (chỉ từ 1,2 đến 2,4 mm/ngày), chủ yếu chỉ có bốc hơi đất ($E$) mà thiếu thành phần thoát hơi thực vật ($T$).
pie title Tỷ lệ bốc thoát hơi nước theo trạng thái sinh trưởng của thảm phủ
    "Thoát hơi nước qua khí khổng lá (T)" : 90
    "Bốc hơi trực tiếp từ bề mặt đất (E)" : 10

[!NOTE] "Khi cây trồng còn nhỏ gần như 100% ET đến từ bốc hơi nước, trong khi ở thời điểm sinh trưởng phát triển mạnh độ che phủ kín diện tích đất thì hơn 90% ET đến từ thoát hơi nước" (Trích dẫn cơ sở sinh lý thảm phủ từ luận án).

Implications đa chiều

  • Đóng góp lý thuyết: Khẳng định tính đúng đắn của việc kết hợp mô hình vật lý cân bằng năng lượng viễn thám với mô hình giải tích bức xạ bán thực nghiệm, hình thành nhánh tiếp cận lai ghép (hybrid approach) mới trong thủy văn địa không gian.
  • Đổi mới phương pháp luận: Thiết lập quy trình tự động hóa khép kín trên Google Earth Engine, cho phép trích xuất $ET_a$ tự động theo chu kỳ lặp của vệ tinh Landsat (8-16 ngày khi kết hợp Landsat 8 và 9), giải phóng sự phụ thuộc vào các trạm đo mặt đất thủ công.
  • Ứng dụng thực tiễn và chính sách:
    • Quản lý tưới tiêu chính xác: Cung cấp bản đồ nhu cầu nước cây trồng đến từng ô đất 30m x 30m, giúp ngành nông nghiệp Hòa Bình tối ưu hóa lịch đóng mở cống tưới.
    • Cảnh báo hạn hán và cháy rừng: Chỉ số $ET_a$ suy giảm bất thường trên các trạng thái rừng phòng hộ là tín hiệu cảnh báo sớm nguy cơ thiếu ẩm đất và vật liệu cháy khô nổ, phục vụ đắc lực cho Chi cục Kiểm lâm và Ban Chỉ huy Phòng chống thiên tai.

Limitations và Future Research

Mặc dù đạt được những bước tiến đột phá, luận án thẳng thắn chỉ ra các giới hạn nội tại:

  1. Rào cản mây che phủ trong ảnh quang học: Vùng Tây Bắc có lượng mây bao phủ lớn trong mùa mưa ẩm, làm gián đoạn chuỗi ảnh quang học Landsat 8, gây khó khăn cho việc quan trắc $ET_a$ liên tục theo ngày.
  2. Độ phân giải thời gian của dữ liệu: Chu kỳ quay trở lại 16 ngày của Landsat 8 (hoặc 8 ngày khi kết hợp Landsat 9) chưa đáp ứng hoàn hảo các yêu cầu theo dõi bốc thoát hơi nước theo từng giờ trong ngày phục vụ điều tiết hồ chứa khẩn cấp.
  3. Phạm vi kiểm chứng thực địa: Nghiên cứu mới thực nghiệm chi tiết tại tỉnh Hòa Bình; các tiểu vùng Tây Bắc khác (như lòng chảo Điện Biên, cao nguyên Mộc Châu, vùng núi cao Lai Châu - Lào Cai) cần tiếp tục mở rộng mạng lưới đo kiểm chứng.

Định hướng nghiên cứu tương lai (Future Research Agenda):

  • Tích hợp kỹ thuật kết hợp dữ liệu đa nguồn (data fusion) giữa ảnh quang học Landsat 8/9, Sentinel-2 với ảnh radar khẩu độ tổng hợp (SAR) Sentinel-1 để xuyên mây trong mùa mưa (tương tự định hướng của Singh et al., 2020).
  • Ứng dụng dữ liệu viễn thám địa tĩnh nhiệt độ cao (như vệ tinh Himawari-8/9) để ngoại suy $ET_a$ theo từng giờ phục vụ mô hình hóa thủy văn thời gian thực.
  • Ứng dụng các thuật toán Trí tuệ nhân tạo (Machine Learning/Deep Learning) để tự động hóa việc nhận diện và hiệu chỉnh hệ số vi khí hậu theo từng vi sinh cảnh phức tạp.

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

  • Ảnh hưởng học thuật: Luận án mở ra hướng nghiên cứu mới trong việc tham số hóa vi khí hậu viễn thám tại các vùng núi nhiệt đới, dự kiến đóng góp vào các diễn đàn khoa học uy tín trong nước và quốc tế về Khoa học Trái đất và Tài nguyên nước.
  • Tác động ngành nông - lâm nghiệp: Cung cấp giải pháp định lượng bốc thoát hơi nước chính xác cho 44,7% diện tích rừng Tây Bắc và hàng trăm nghìn hecta hoa màu lưu vực sông Đà, giúp cắt giảm 15-20% lượng nước tưới tiêu lãng phí nhờ điều tiết theo nhu cầu thực tế.
  • Chuyển đổi số trong quản trị tài nguyên: Hỗ trợ Sở Nông nghiệp & PTNT, Sở Tài nguyên & Môi trường các tỉnh Tây Bắc chuyển từ hình thức kiểm tra thủ công sang giám sát số hóa trực quan trên bản đồ WebGIS/GEE với chi phí vận hành giảm trên 60%.
  • Ý nghĩa quốc tế: Phương pháp luận tích hợp SEBAL - Priestley - Taylor có khả năng chuyển giao và áp dụng rộng rãi cho các quốc gia có địa hình đồi núi tương đồng trong khu vực Đông Nam Á và lưu vực sông Mê Kông.

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

  • Nghiên cứu sinh & Nhà khoa học trẻ: Tiếp cận quy trình giải đoán tham số nhiệt viễn thám và thuật toán cân bằng bức xạ trên nền tảng đám mây GEE để kế thừa phát triển các đề tài thủy văn chuyên sâu.
  • Giảng viên & Chuyên gia Viễn thám/GIS: Nguồn tài liệu tham khảo chuẩn mực về phương pháp kết hợp mô hình viễn thám với mô hình thực nghiệm địa phương hóa.
  • Cơ quan Quản lý Thủy lợi & Tài nguyên Nước: Tiếp nhận công cụ giám sát lượng bốc thoát hơi nước diện rộng, hỗ trợ tính toán cân bằng nước lưu vực sông và quy hoạch hồ chứa thủy điện/thủy lợi.
  • Chi cục Kiểm lâm & Phòng chống cháy rừng: Công cụ theo dõi độ ẩm thảm phủ vi mô từ xa để khoanh vùng các điểm nóng có nguy cơ cháy rừng cao trong mùa khô nóng Tây Bắc.

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?

Đóng góp lý thuyết độc đáo nhất là việc thiết lập thành công cầu nối giải tích giữa Thuật toán cân bằng năng lượng bề mặt viễn thám (SEBAL) và Mô hình bức xạ Priestley - Taylor (1972) trên nền tảng mô hình số độ cao DEM. Luận án đã chứng minh về mặt khoa học rằng toàn bộ 5 tham số cốt lõi phục vụ tính toán $ET_a$ ($R_{nd}, T_s, \lambda, \gamma, \Delta$) hoàn toàn có thể trích xuất ở mức độ phân giải 30m từ ảnh Landsat 8 mà không cần dựa vào chuỗi đo nhiệt ẩm mặt đất. Công trình đã mở rộng lý thuyết Priestley - Taylor từ một mô hình điểm (point-based model) dựa vào trạm thời tiết thành một mô hình phân tán không gian (spatially-distributed model) có khả năng tự thích ứng với địa hình núi cao nhiệt đới.

2. Tính đổi mới trong phương pháp luận so với các nghiên cứu trước đây thể hiện như thế nào?

So với nghiên cứu của Nguyễn Văn Hùng (2017) trên lưu vực Sông Cầu (dùng S-SEBI và chỉ thử nghiệm hệ số giả định với sai số 11% - 18%) và nghiên cứu của Trần Thị Tâm et al. (2016) trên ảnh MODIS (độ phân giải thấp 1km, độ lệch lớn tới 230 mm/năm):

  • Luận án đã khai thác triệt để độ phân giải không gian 30m của Landsat 8 kết hợp kênh hồng ngoại nhiệt TIRS.
  • Phương pháp luận tiến hành giải mã bức xạ ròng theo chuẩn cân bằng năng lượng thực, sau đó chuẩn hóa và kiểm chứng với tiêu chuẩn FAO-56.
  • Đổi mới căn bản nằm ở việc xác định thực nghiệm cặp hệ số ($a, b$) từ chuỗi dữ liệu đa năm (2015–2021) tại Hòa Bình, đưa sai số toàn diện của mô hình xuống dưới 10%.

3. Phát hiện bất ngờ nhất từ thực nghiệm dữ liệu là gì?

Phát hiện bất ngờ và có ý nghĩa thực tiễn sâu sắc là: Mặc dù khu vực Tây Bắc có nền nhiệt độ trung bình cao hơn vùng Đông Bắc từ 2 - $3^\circ\text{C}$ và số giờ nắng lớn, nhưng lượng bốc thoát hơi nước thực tế ($ET_a$) tại các vùng đồi núi dốc trơ sỏi đá lại thấp hơn rất nhiều so với dự báo từ các mô hình khí tượng thuần túy. Lý do là độ ẩm đất cạn kiệt đã khống chế bốc hơi đất ($E \approx 0$), trong khi thảm thực vật bị suy thoái không thể thực hiện thoát hơi ($T \approx 0$). Ngược lại, tại các thảm rừng phòng hộ có độ che phủ cao (đạt trung bình 44,7% toàn vùng), quá trình thoát hơi nước qua khí khổng thực vật chiếm tới hơn 90% tổng lượng $ET$, đóng vai trò máy điều hòa sinh thái giữ ẩm cho toàn bộ lưu vực.

4. Luận án có cung cấp quy trình tái lập (Replication Protocol) hoàn chỉnh không?

Hoàn toàn có. Luận án đã chuẩn hóa một quy trình công nghệ 4 bước chặt chẽ:

  • Bước 1: Thu nhận ảnh Landsat 8 C1/C2 và DEM SRTM 30m tự động trên Google Earth Engine.
  • Bước 2: Chạy mã nguồn thuật toán SEBAL chiết xuất các layer raster ($\alpha, NDVI, SAVI, T_s, R_n, R_{nd}, \lambda, \gamma, \Delta$).
  • Bước 3: Nạp bộ hệ số thực nghiệm ($a, b$) đã tối ưu hóa cho điều kiện vùng nghiên cứu.
  • Bước 4: Tính toán tự động raster $ET_a$ (mm/ngày) và xuất báo cáo không gian. Toàn bộ quy trình có thể tái lập 100% cho bất kỳ tỉnh nào thuộc khu vực miền núi phía Bắc Việt Nam.

5. Lộ trình nghiên cứu 10 năm (10-year Research Agenda) được phác thảo ra sao?

  • Giai đoạn 1 (2026–2028): Mở rộng mạng lưới hiệu chuẩn hệ số ($a, b$) cho toàn bộ 6 tỉnh Tây Bắc (Hòa Bình, Sơn La, Điện Biên, Lai Châu, Lào Cai, Yên Bái); hoàn thiện module hợp nhất ảnh Landsat 8/9 và Sentinel-2.
  • Giai đoạn 2 (2029–2032): Tích hợp ảnh radar viễn thám Sentinel-1/SAR để giải quyết triệt để bài toán mù mây trong mùa mưa bão; xây dựng hệ thống giám sát hạn hán nông - lâm nghiệp tự động cập nhật hàng tuần.
  • Giai đoạn 3 (2033–2036): Phát triển nền tảng Trí tuệ nhân tạo địa không gian (Geo-AI) dự báo nhu cầu tiêu thụ nước cây trồng và cảnh báo cháy rừng sớm theo thời gian thực (real-time) dưới các kịch bản biến đổi khí hậu cực đoan.

Kết luận

Luận án "Nghiên cứu xây dựng mô hình giám sát sự bốc - thoát hơi nước của lớp phủ khu vực Tây Bắc Việt Nam từ dữ liệu ảnh vệ tinh" là một công trình khoa học công phu, nghiêm túc và mang tính thời sự cao. Nghiên cứu đã hoàn thành xuất sắc các mục tiêu đặt ra với các đóng góp cốt lõi:

  1. Chuẩn hóa cơ sở khoa học: Chứng minh tính khả thi tuyệt đối của việc sử dụng dữ liệu viễn thám Landsat 8 kết hợp mô hình số độ cao DEM để chiết xuất trọn bộ các tham số nhiệt vi khí hậu ($R_{nd}, T_s, \lambda, \gamma, \Delta$) phục vụ tính toán $ET_a$ vùng núi cao.
  2. Đột phá về tham số địa phương hóa: Xác định thành công cặp hệ số thực nghiệm ($a, b$) của mô hình Priestley - Taylor từ chuỗi số liệu khí tượng thủy văn 2015–2021, phản ánh trung thực đặc tính vi khí hậu và lớp phủ tỉnh Hòa Bình.
  3. Thiết lập mô hình lai ưu việt: Đề xuất thành công mô hình tích hợp SEBAL - Priestley - Taylor, khống chế sai số tính toán $ET_a$ dưới 10% so với số liệu quan trắc thực địa.
  4. Tự động hóa công nghệ cao: Xây dựng hoàn chỉnh quy trình và chương trình giám sát bốc thoát hơi nước tự động trên nền tảng điện toán đám mây Google Earth Engine.
  5. Giá trị ứng dụng thực tiễn to lớn: Cung cấp công cụ số hóa phục vụ đắc lực cho công tác quản lý tài nguyên nước, điều tiết thủy lợi, dự báo hạn hán và phòng chống cháy rừng tại khu vực Tây Bắc Việt Nam.

Công trình không chỉ tạo tiền đề vững chắc cho việc hiện đại hóa công tác giám sát tài nguyên nước quốc gia mà còn mở ra hướng tiếp cận học thuật tiên tiến cho cộng đồng nghiên cứu Trắc địa - Viễn thám và Thủy văn trong kỷ nguyên số.