Tổng quan về luận án

Nghiên cứu giám sát động thái bốc - thoát hơi nước ($ET$ - Evapotranspiration) đóng vai trò then chốt trong việc phân tích chu trình thủy văn, quản lý an ninh nguồn nước, dự báo hạn hán và phòng chống cháy rừng trong bối cảnh biến đổi khí hậu toàn cầu. Tại các khu vực địa hình đồi núi phức tạp như vùng Tây Bắc Việt Nam—nơi có độ chia cắt địa hình sâu sắc, chênh cao lớn với dãy Hoàng Liên Sơn dài 180 km (đỉnh cao trên 2.000 m), dãy Sông Mã dài 500 km (đỉnh cao trên 1.800 m) xen kẽ các cao nguyên đá vôi Mộc Châu, Tà Phình, Nà Sản và độ che phủ rừng đạt 44,7%—việc định lượng chính xác lượng bốc thoát hơi nước thực tế ($ET_a$) bằng mạng lưới trạm quan trắc mặt đất truyền thống gặp phải những rào cản nghiêm trọng về chi phí, nhân lực và tính đại diện không gian.

Khoảng trống nghiên cứu (research gap) cốt lõi được xác định: Mạng lưới trạm khí tượng thủy văn tại Việt Nam nói chung và khu vực Tây Bắc nói riêng rất thưa thớt (trung bình dưới 10 trạm trên phạm vi một tỉnh), dữ liệu thu thập thủ công dạng điểm (point-scale) không thể phản ánh chính xác sự biến thiên liên tục của $ET_a$ trên các bề mặt lớp phủ không đồng nhất và địa hình chia cắt mạnh. Ngược lại, các mô hình viễn thám cân bằng năng lượng bề mặt kinh điển như SEBAL (Surface Energy Balance Algorithms for Land - Bastiaanssen et al., 1998) hoặc METRIC (Mapping ET with Internalized Calibration - Allen et al., 2007) khi áp dụng độc lập vào khu vực đồi núi nhiệt đới thường bộc lộ hạn chế lớn do việc xác định thủ công các điểm ảnh giới hạn "nóng" (pixel khô) và "lạnh" (pixel ướt) mang tính chủ quan cao, đồng thời việc tính toán thông lượng nhiệt cảm ứng ($H$) đòi hỏi cấu trúc tham số khí động học phức tạp dễ dẫn đến sai số lan truyề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" do nghiên cứu sinh Lê Hùng Chiến thực hiện (chuyên ngành Kỹ thuật Trắc địa - Bản đồ, mã số: 9.520503; người hướng dẫn khoa học: PGS.TS. Trần Xuân Trường và PGS.TS. Doãn Hà Phong, Trường Đại học Mỏ - Địa chất, 2022) đã giải quyết triệt để khoảng trống này thông qua hai câu hỏi nghiên cứu và hai giả thuyết khoa học:

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

Khung lý thuyết của công trình tích hợp lý thuyết cân bằng bức xạ bề mặt đất (Bastiaanssen et al., 1998), phương trình bức xạ cân bằng Priestley - Taylor (1972), tiêu chuẩn tham chiếu FAO-56 Penman - Monteith (Allen et al., 1998) và nhiệt động học khí quyển hiệu chỉnh theo mô hình độ cao số (SRTM DEM). Phạm vi thực nghiệm của luận án tập trung tại tỉnh Hòa Bình (tọa độ địa lý đặc trưng cho vùng Tây Bắc), sử dụng chuỗi dữ liệu khí tượng thủy văn giai đoạn 2015-2021 kết hợp các cảnh ảnh vệ tinh Landsat 8 vào ba mốc thời gian đại diện mùa vụ (01/07/2015, 04/06/2017 và 18/08/2021).

Literature Review và Positioning

Lịch sử nghiên cứu bốc thoát hơi nước bề mặt trải qua ba dòng tiếp cận chính với sự đóng góp của các học giả quốc tế và trong nước:

Dòng tiếp cận thứ nhất dựa trên các phép đo vi khí tượng và cân bằng thủy văn mặt đất. Allen et al. (1990, 1998) định nghĩa bốc hơi ($E$) là sự khuếch tán của phân tử nước từ đất, khối nước và bề mặt ẩm ướt vào khí quyển, trong khi thoát hơi ($T$) là phản ứng sinh lý của thực vật qua lỗ khí khổng để chống chịu sự khô hạn. Để chuẩn hóa nhu cầu bốc thoát hơi của khí quyển độc lập với loại cây trồng, Doorenbos và Pruitt (1975) cùng Allen et al. (1998) đã chuẩn hóa khái niệm bốc thoát hơi nước tham chiếu ($ET_0$) thông qua mô hình FAO-56 Penman - Monteith. Mặc dù các công cụ đo trực tiếp như ống đo bốc hơi Piche hay thủy tiêu kế (Lysimeter - Lê Anh Tuấn, 2009) cung cấp độ chính xác cao tại chỗ, chúng đòi hỏi kinh phí duy trì rất lớn và hoàn toàn bất khả thi trong việc mở rộng quy mô không gian cho các lưu vực sông miền núi.

Dòng tiếp cận thứ hai là các mô hình ước tính dựa trên bức xạ mặt trời (radiation-based models). Makkink (1957) đề xuất công thức tính $ET_a$ qua bức xạ sóng ngắn đến ($R_s$) với hệ số $u = 0,61; v = 0,12$, sau đó được Hansen (1984) hiệu chỉnh tại Hà Lan ($u = 0,70; v = 0$) và Xue & Singh (1999) khảo sát ($u = 0,77; v = 0,22$). Priestley và Taylor (1972) đưa ra bước đột phá khi biểu diễn $ET_a$ như một hàm của bức xạ ròng ($R_n$), độ dốc áp suất hơi bão hòa ($\Delta$) và hằng số tâm lý ẩm ($\gamma$): $$ET_a = a \frac{\Delta}{\Delta + \gamma} \frac{R_n}{\lambda} + b$$ Hệ số $(a, b)$ ban đầu được Priestley - Taylor xác định là $(0,61; 0,12)$, được kiểm định tại Thụy Sĩ (1984) là $(0,90; 0)$ và tại Đài Loan (2005) là $(1,00; 0)$. Xu và Singh (2001) khi đánh giá 4 mô hình bức xạ so với chậu đo bốc hơi Pan đã kết luận rằng mô hình Makkink và Priestley - Taylor mang lại độ tin cậy cao nhất trong việc mô phỏng $ET_0$. Ngoài ra, Irmak et al. (2003) đề xuất mô hình thực nghiệm đơn giản dựa trên bức xạ ròng ngày ($R_{nd}$) và nhiệt độ không khí ($T$).

Dòng tiếp cận thứ ba khai thác ảnh viễn thám nhiệt và cân bằng năng lượng bề mặt (Surface Energy Balance). Điển hình là mô hình SEBAL (Bastiaanssen et al., 1998, 2000), SEBI (Choudhury & Menenti, 1993), S-SEBI (Roerink et al., 2000), SEBS (Su, 2002) và METRIC (Allen et al., 2007). Bastiaanssen et al. (1998) chứng minh SEBAL đạt độ chính xác 85% theo ngày và 95% theo mùa trên các lưu vực bằng phẳng. Tuy nhiên, Trezza (2002) ghi nhận sai số SEBAL dao động từ 2,7% đến 35% (trung bình 18,2%) khi áp dụng cho các vùng cây trồng phức tạp, và Bastiaanssen (1995) thừa nhận mô hình tồn tại nhiều hạn chế tại khu vực đồi núi hiểm trở.

So sánh với các nghiên cứu quốc tế đương đại:

  • Ghaderi et al. (2020) áp dụng SEBAL tại đồng bằng Ein Khosh (Iran), đạt $RMSE = 0,5\text{ mm/ngày}$, $MAPE = 2,9%$ và hệ số tương quan $R = 0,97$ khi đối chiếu với FAO-56 Penman - Monteith. Tuy nhiên, vùng nghiên cứu của Ghaderi là địa hình sa mạc đồng bằng đồng nhất.
  • Bruno et al. (2018) kết hợp Landsat 8 và trạm khí tượng tại Brazil cho thấy sai số ước tính ET khoảng $1,00\text{ mm/ngày}$, nhưng phụ thuộc chặt vào mật độ trạm mặt đất.
  • Madugundu et al. (2017) dùng METRIC trên cánh đồng cỏ linh lăng tại Ả Rập Xê Út đạt RMSE $0,13 - 4,15\text{ mm/ngày}$, đòi hỏi quy trình hiệu chuẩn nội bộ khắt khe.

Vị thế khoa học (positioning) của luận án: Luận án không áp dụng máy móc mô hình cân bằng năng lượng viễn thám nguyên bản (SEBAL/METRIC) vốn dễ sai lệch trong việc lặp tính thông lượng nhiệt cảm ứng $H$ ở địa hình đồi núi; thay vào đó, tác giả thực hiện một giải pháp đột phá: Khai thác module cân bằng bức xạ của SEBAL để tính toán năng lượng bức xạ ròng trung bình ngày ($R_{nd}$), sau đó tích hợp vào mô hình Priestley - Taylor đã được tối ưu hóa hệ số $(a, b)$ địa phương kết hợp chiết xuất nhiệt động học từ Landsat 8 và SRTM DEM.

Đó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à hoàn thiện lý thuyết truyền nhiệt viễn thám và vi khí tượng bề mặt trong các khía cạnh:

  1. Tái cấu trúc và địa phương hóa phương trình bức xạ Priestley - Taylor (1972): Luận án chứng minh rằng các hệ số mặc định toàn cầu của Priestley - Taylor ($a = 1,26$ hoặc $a = 0,61; b = 0,12$) không phản ánh đúng động thái bốc thoát hơi nước trong điều kiện thảm thực vật đa tầng tán và độ ẩm gió mùa nhiệt đới Tây Bắc Việt Nam. Việc xác định cặp hệ số thực nghiệm $(a, b)$ thông qua chuỗi số liệu khí tượng thủy văn nhiều năm đã cung cấp một cơ sở lý thuyết chuẩn xác cho vùng sinh thái đồi núi Đông Nam Á.
  2. Xác lập tiên đề thay thế tham số khí tượng bằng viễn thám: Luận án chứng minh một cách chặt chẽ rằng bốn biến nhiệt động học then chốt—gồm năng lượng bức xạ ròng ngày ($R_{nd}$), nhiệt hóa hơi ($\lambda$), hằng số tâm lý ẩm ($\gamma$) và độ dốc áp suất hơi bão hòa ($\Delta$)—có thể được giải tích hóa hoàn toàn từ kênh phổ phản xạ (Bands 2-7), kênh hồng ngoại nhiệt (Band 10/11) của Landsat 8 và độ cao tuyệt đối ($z$) từ SRTM DEM mà không làm suy giảm độ chính xác so với phương trình tiêu chuẩn FAO-56.

Mô hình lý thuyết được thiết lập qua hai mệnh đề cốt lõi:

  • Mệnh đề 1 (P1): Giá trị bức xạ ròng tức thời ($R_{ni}$) và trung bình ngày ($R_{nd}$) trích xuất từ dữ liệu viễn thám quang học/nhiệt phản ánh chính xác cấu trúc cân bằng bức xạ bề mặt theo định luật Stefan-Boltzmann và định luật bảo toàn năng lượng.
  • Mệnh đề 2 (P2): Sự suy giảm của áp suất khí quyển $P(z)$ theo độ cao địa hình chi phối trực tiếp hằng số tâm lý ẩm $\gamma(z)$, tạo ra mối quan hệ phi tuyến giữa $ET_a$ và mô hình số độ cao DEM.

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

Khung phân tích của luận án tích hợp đồng thời ba cấu trúc lý thuyết: (1) Mô hình cân bằng bức xạ SEBAL, (2) Mô hình bức xạ tương đối Priestley - Taylor, và (3) Mô hình nhiệt động học độ cao khí quyển.

                  +-----------------------------------+
                  |      Ảnh vệ tinh Landsat 8        |
                  | (Bands 2-7, Band 10) + SRTM DEM   |
                  +-----------------+-----------------+
                                    |
          +-------------------------+-------------------------+
          |                                                   |
          v                                                   v
+-----------------------+                           +-------------------+
|  Mô hình bức xạ SEBAL |                           | Tham số nhiệt động|
|  - Albedo (α), NDVI,  |                           | học bề mặt & DEM  |
|    SAVI, LAI          |                           | - Ts (Band 10)    |
|  - ε0, εNB, RL↑, RL↓  |                           | - λ = f(Ts)       |
|  - Rni -> Rnd_VT      |                           | - γ = f(DEM, Ts)  |
+-----------+-----------+                           | - Δ = f(Ts)       |
            |                                       +---------+---------+
            |                                                 |
            +-----------------------+-------------------------+
                                    |
                                    v
                     +-----------------------------+
                     |  Mô hình Priestley - Taylor |
                     |   Hiệu chuẩn hệ số (a, b)   |
                     |         khu vực Tây Bắc     |
                     +--------------+--------------+
                                    |
                                    v
                     +-----------------------------+
                     | Bản đồ phân bố không gian   |
                     | lượng bốc thoát hơi nước    |
                     | thực tế ETa (mm/ngày)       |
                     +-----------------------------+

Phương trình tổng quát của khung phân tích được phát biểu: $$ET_a = a \left[ \frac{\Delta(T_s)}{\Delta(T_s) + \gamma(z)} \right] \frac{R_{nd_VT}}{\lambda(T_s)} + b$$ Trong đó:

  • $T_s$: Nhiệt độ bức xạ bề mặt đất trích xuất từ kênh nhiệt Landsat 8 (K).
  • $z$: Độ cao tuyệt đối chiết xuất từ mô hình số độ cao SRTM DEM (m).
  • $R_{nd_VT}$: Năng lượng bức xạ ròng trung bình ngày tính từ viễn thám ($\text{MJ/m}^2\text{/ngày}$).
  • $\lambda(T_s)$: Nhiệt ẩn hóa hơi ($\text{MJ/kg}$).
  • $\gamma(z)$: Hằng số tâm lý ẩm hiệu chỉnh theo áp suất không khí tại độ cao $z$ ($\text{kPa/}^\circ\text{C}$).
  • $\Delta(T_s)$: Độ dốc đường cong áp suất hơi bão hòa tại nhiệt độ bề mặt ($\text{kPa/}^\circ\text{C}$).
  • $(a, b)$: Cặp hệ số thực nghiệm tối ưu cho vùng Tây Bắc.

Điều kiện biên (boundary conditions): Khung phân tích áp dụng tối ưu cho các cảnh ảnh có tỷ lệ mây dưới 10-20%, thảm phủ nhiệt đới đa dạng (rừng tự nhiên, rừng trồng, đất nông nghiệp, mặt nước, đô thị) trong phạm vi độ cao từ 50 m đến 2.500 m so với mực nước biển.

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

Thiết kế nghiên cứu

Nghiên cứu tuân thủ thế giới quan thực chứng (positivism) và phương pháp định lượng vật lý - thực nghiệm nghiêm ngặt. Thiết kế nghiên cứu đa cấp độ (multi-level design) kết hợp xử lý dữ liệu không gian vĩ mô ở độ phân giải 30m trên toàn tỉnh Hòa Bình và đối soát vi mô tại các trạm khí tượng bề mặt.

Không gian thực nghiệm bao gồm toàn bộ diện tích tự nhiên tỉnh Hòa Bình với mạng lưới 5 trạm khí tượng thủy văn đại diện cho các tiểu vùng địa hình: Trạm Hòa Bình, Chi Nê, Kim Bôi, Lạc Sơn và Mai Châu. Chuỗi số liệu khí tượng quan trắc ngày được thu thập liên tục trong 6 năm (2015-2021), kết hợp các cảnh ảnh vệ tinh Landsat 8 OLI/TIRS chất lượng cao (WGS-84, múi chiếu UTM 48N / múi 6 độ) thu nhận tại các thời điểm điển hình: 01/07/2015, 04/06/2017 và 18/08/2021.

Quy trình nghiên cứu rigorous

Quy trình xử lý dữ liệu viễn thám và tính toán tham số được chuẩn hóa qua các bước giải tích:

  1. Tiền xử lý và Hiệu chỉnh Khí quyển: Chuyển đổi giá trị số điểm ảnh ($DN$) sang năng lượng bức xạ phổ tại đỉnh khí quyển ($L_\lambda$), sau đó áp dụng mô hình truyền xạ FLAASH (Fast Line-of-sight Atmospheric Analysis of Spectral Hypercubes) để triệt tiêu ảnh hưởng của tán xạ và hấp thụ khí quyển, thu được hệ số phản xạ bề mặt ($\rho_\lambda$) cho các kênh quang học.
  2. Chiết xuất các chỉ số thực vật và tham số bề mặt:
    • Chỉ số thực vật khác biệt chuẩn hóa: $NDVI = \frac{\rho_{NIR} - \rho_{Red}}{\rho_{NIR} + \rho_{Red}}$
    • Chỉ số thực vật điều chỉnh ảnh hưởng của đất: $SAVI = \frac{(1 + L)(\rho_{NIR} - \rho_{Red})}{\rho_{NIR} + \rho_{Red} + L}$ (với $L = 0,5$)
    • Chỉ số diện tích lá ($LAI$) và Suất phân sai bề mặt ($\alpha$ - Albedo): $$\alpha = \sum (w_\lambda \cdot \rho_\lambda)$$
    • Hệ số phát xạ nhiệt bề mặt dải hẹp ($\varepsilon_{NB}$) và dải rộng ($\varepsilon_0$) được xác định từ hàm phi tuyến của $NDVI$ và $LAI$.
    • Nhiệt độ bề mặt ($T_s$) được tính từ kênh hồng ngoại nhiệt (Band 10) thông qua việc nghịch đảo hàm bức xạ Planck có hiệu chỉnh phát xạ $\varepsilon_{NB}$.
  3. Mô hình hóa bức xạ ròng ($R_{ni}$ và $R_{nd}$): Phương trình cân bằng bức xạ tức thời tại thời điểm chụp ảnh: $$R_{ni} = (1 - \alpha) R_{S\downarrow} + R_{L\downarrow} - R_{L\uparrow} - (1 - \varepsilon_0) R_{L\downarrow}$$ Trong đó: Bức xạ sóng ngắn tới ($R_{S\downarrow}$) tính qua hằng số mặt trời ($G_{sc} = 1367\text{ W/m}^2$), góc thiên đỉnh mặt trời ($\theta$) và độ truyền dẫn khí quyển ($\tau_{sw}$). Bức xạ sóng dài tới ($R_{L\downarrow}$) và phát xạ sóng dài đi ($R_{L\uparrow}$) được xác định theo định luật Stefan-Boltzmann với nhiệt độ bề mặt $T_s$ và nhiệt độ không khí $T_a$. Chuyển đổi bức xạ ròng tức thời ($R_{ni}$) sang trung bình ngày ($R_{nd_VT}$): $$R_{nd_VT} = \frac{24}{\pi \sin\left(\frac{\pi \cdot t_d}{N}\right)} R_{ni} \cdot \frac{\tau_{sw_d}}{\tau_{sw}}$$
  4. Chiết xuất tham số nhiệt động học kết hợp DEM:
    • Hằng số nhiệt ẩn hóa hơi: $\lambda = 2,501 - (2,361 \times 10^{-3}) T_s$ ($\text{MJ/kg}$)
    • Áp suất khí quyển theo độ cao $z$ từ SRTM DEM: $$P = 101,3 \left( \frac{293 - 0,0065 z}{293} \right)^{5,26} \text{ (kPa)}$$
    • Hằng số tâm lý ẩm: $\gamma = 0,665 \times 10^{-3} P = 0,665 \times 10^{-3} \left[ 101,3 \left( \frac{293 - 0,0065 z}{293} \right)^{5,26} \right]$
    • Độ dốc đường cong áp suất hơi bão hòa: $$\Delta = \frac{4098 \left[ 0,6108 \exp\left( \frac{17,27 T_s}{T_s + 237,3} \right) \right]}{(T_s + 237,3)^2}$$
  5. Tam giác đạc dữ liệu (Triangulation) và Kiểm định: So sánh kết quả $R_{nd_VT}$ với $R_{nd_FAO}$ tính theo tiêu chuẩn FAO-56 từ số liệu quan trắc trạm. So sánh lượng bốc thoát hơi nước mô hình $ET_{a_VT}$ với số liệu thực tế đo bằng chậu Pan/Piche tại các trạm khí tượng ($ET_{a_Đo}$).

Data và phân tích

Toàn bộ thuật toán được lập trình tự động hóa trên nền tảng Google Earth Engine (GEE API Javascript), cho phép người dùng định nghĩa khu vực quan tâm (ROI), lọc ảnh Landsat 8 theo tỷ lệ mây, tính toán chuỗi chỉ số và trích xuất bản đồ phân bố $ET_a$ dạng raster (GeoTIFF) và bảng thống kê giá trị điểm.

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

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

Kết quả thực nghiệm của luận án mang lại 4 phát hiện khoa học mang tính đột phá:

  1. Khả năng thay thế hoàn hảo của bức xạ ròng viễn thám ($R_{nd_VT}$): Kết quả so sánh giữa bức xạ ròng trung bình ngày tính từ ảnh Landsat 8 ($R_{nd_VT}$) và tính từ số liệu khí tượng trạm theo chuẩn FAO-56 ($R_{nd_FAO}$) tại 5 trạm quan trắc trên địa bàn tỉnh Hòa Bình trong 3 thời điểm nghiên cứu (01/07/2015, 04/06/2017, 18/08/2021) cho thấy sự tương đồng rất cao. Chênh lệch tuyệt đối giữa $R_{nd_VT}$ và $R_{nd_FAO}$ tại hầu hết các trạm chỉ dao động trong khoảng $0,05 - 0,42\text{ MJ/m}^2\text{/ngày}$ (sai số tương đối dưới 3,5%), chứng minh thuật toán cân bằng bức xạ SEBAL hoàn toàn đủ độ tin cậy để chiết xuất năng lượng mặt trời hấp thụ tại mặt đất mà không cần trạm đo bức xạ chuyên dụng.

  2. Xác lập thành công cặp hệ số thực nghiệm $(a, b)$ tối ưu cho vùng Tây Bắc: Qua phân tích hồi quy tương quan đa biến từ chuỗi số liệu khí tượng thủy văn 2015-2021, luận án đã xác định được bộ hệ số thực nghiệm cho mô hình Priestley - Taylor phù hợp đặc thù khí hậu nhiệt đới gió mùa đồi núi Tây Bắc. Khác biệt rõ rệt với các hệ số kinh điển ở vùng ôn đới ($a = 0,61; b = 0,12$), bộ hệ số mới phản ánh chính xác tác động của độ ẩm cao và bức xạ nhiệt đới, loại bỏ hoàn toàn hiện tượng ước tính thiếu (underestimation) vào mùa hè và ước tính thừa (overestimation) vào mùa đông.

  3. Độ chính xác vượt trội của mô hình lai ghép $ET_{a_VT}$ so với dữ liệu đo thực tế: Khi đối chiếu lượng bốc thoát hơi nước thực tế $ET_{a_VT}$ (tính toán hoàn toàn từ ảnh Landsat 8 và SRTM DEM) với giá trị thực đo $ET_{a_Đo}$ tại các trạm khí tượng Hòa Bình, Chi Nê, Kim Bôi, Lạc Sơn, Mai Châu:

    • Thời điểm ngày 01/07/2015: Giá trị $ET_{a_VT}$ dao động từ $3,82 - 4,45\text{ mm/ngày}$, độ lệch so với trạm đo $ET_{a_Đo}$ ($3,90 - 4,50\text{ mm/ngày}$) chỉ từ $0,05 - 0,18\text{ mm/ngày}$ (sai số dưới 4,2%).
    • Thời điểm ngày 04/06/2017: $ET_{a_VT}$ đạt $3,65 - 4,28\text{ mm/ngày}$, chênh lệch so với trạm đo không quá $0,21\text{ mm/ngày}$.
    • Thời điểm ngày 18/08/2021: $ET_{a_VT}$ đạt $3,40 - 4,12\text{ mm/ngày}$, tương thích tuyệt vời với thực tế quan trắc mặt đất.

  4. Đặc trưng phân hóa không gian của $ET_a$ theo cấu trúc lớp phủ và địa hình: Bản đồ $ET_a$ thành lập cho thấy sự phân hóa không gian sâu sắc: Lượng bốc thoát hơi nước đạt giá trị cực đại ($4,5 - 5,8\text{ mm/ngày}$) tại các vùng rừng đặc dụng, rừng phòng hộ đầu nguồn có tầng tán dày và độ ẩm đất cao (khu bảo tồn Ngọc Sơn - Ngổ Luông, Thượng Tiến, Hang Kia - Bà Cò); đạt mức trung bình ($3,2 - 4,2\text{ mm/ngày}$) trên đất trồng cây nông nghiệp, cây ăn quả và lúa nước lưu vực sông Đà; và giảm mạnh ($1,5 - 2,5\text{ mm/ngày}$) tại các diện tích đất trống, đồi trọc, vùng đô thị hóa cao hoặc các đỉnh núi đá vôi trơ sỏi đá.

Implications đa chiều

  • Ý nghĩa Lý luận (Theoretical Advances): Giải quyết dứt điểm nghịch lý giữa tính phức tạp của phương trình truyền nhiệt vi mô và tính khả thi trong mô hình hóa vĩ mô. Kết hợp thành công nguyên lý cân bằng năng lượng viễn thám với mô hình giải tích Priestley - Taylor, mở ra hướng tiếp cận mới trong mô hình hóa thủy văn viễn thám tại các vùng sinh thái đồi núi nhiệt đới.
  • Đột phá Phương pháp luận (Methodological Innovations): Xây dựng quy trình khép kín chiết xuất 100% tham số động học ($R_{nd}, \lambda, \gamma, \Delta$) từ dữ liệu viễn thám mở (Landsat 8 và SRTM DEM). Việc tích hợp thuật toán trên Google Earth Engine loại bỏ rào cản tính toán phần cứng cục bộ, cho phép xử lý dữ liệu không gian lớn trong vài giây.
  • Ứng dụng Thực tiễn (Practical Applications): Cung cấp công cụ phục vụ quy hoạch tưới tiêu tiết kiệm nước cho cây trồng nông nghiệp (cam Cao Phong, mía tím Tân Lạc), đánh giá nhu cầu nước sinh thái lưu vực hồ chứa thủy điện Hòa Bình và hệ thống sông Đà.
  • Hàm ý Chính sách (Policy Recommendations): Cung cấp cơ sở khoa học cho Bộ Nông nghiệp & Phát triển Nông thôn và Sở Tài nguyên & Môi trường các tỉnh Tây Bắc trong việc xây dựng bản đồ cảnh báo sớm hạn hán nông nghiệp và bản đồ nguy cơ cháy rừng theo thời gian thực dựa trên chỉ số thiếu hụt độ ẩm bốc thoát hơi nước.

Limitations và Future Research

Mặc dù đạt được những kết quả đột phá, luận án thẳng thắn chỉ ra các giới hạn nghiên cứu:

  1. Độ phân giải thời gian của ảnh Landsat 8: Chu kỳ quay vòng 16 ngày của vệ tinh Landsat 8 cùng với xác suất mây mù cao tại vùng núi Tây Bắc trong mùa mưa gây khó khăn cho việc giám sát liên tục theo chuỗi ngày kế tiếp.
  2. Ảnh hưởng của góc nghiêng sườn dốc và bóng địa hình: Mô hình cân bằng bức xạ SEBAL chưa tích hợp sâu mô hình hiệu chỉnh bức xạ địa hình phức tạp (topographic solar radiation modeling) trên các sườn núi dốc đứng hướng khuất nắng, có thể dẫn đến sai số cục bộ tại các vách núi hiểm trở.
  3. Số lượng trạm kiểm định thực nghiệm: Mạng lưới 5 trạm khí tượng tại Hòa Bình tuy phân bố đều nhưng chủ yếu đặt tại các thung lũng hoặc vùng đồi thấp ven sông, thiếu vắng các trạm quan trắc vi khí tượng trên các đỉnh núi cao trên 1.500 m để kiểm chứng cực hạn.

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

  • Tích hợp dữ liệu đa vệ tinh (Multi-sensor Data Fusion): Kết hợp Landsat 8/9, Sentinel-2 (quang học 10-20m) với Sentinel-3 và MODIS/Himawari-8 (độ phân giải thời gian 10 phút đến 1 ngày) để xây dựng chuỗi ảnh $ET_a$ độ phân giải không - thời gian cao không gián đoạn.
  • Ứng dụng Trí tuệ nhân tạo (AI/Machine Learning): Sử dụng các mô hình học sâu (Deep Learning - CNN, LSTM) để tái tạo các vùng ảnh bị che phủ bởi mây (cloud gap-filling).
  • Mở rộng mô hình hóa 3D: Bổ sung thuật toán hiệu chỉnh bức xạ trực tiếp - khuếch tán theo độ dốc và hướng dốc địa hình 3D cho toàn bộ vòng cung Tây Bắc.

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

  • Tác động Học thuật: Cung cấp tài liệu tham khảo nền tảng cho các viện nghiên cứu, trường đại học trong ngành Trắc địa - Bản đồ, Địa tin học, Khí tượng Thủy văn và Quản lý Tài nguyên Nước. Dự báo công trình sẽ thu hút trích dẫn khoa học cao trong các nghiên cứu ứng dụng GEE và viễn thám nhiệt tại khu vực Đông Nam Á.
  • Chuyển đổi Ngành & R&D: Thúc đẩy chuyển đổi số trong công tác điều hành thủy lợi, tối ưu hóa vận hành điều tiết các hồ chứa thủy điện lớn (Hòa Bình, Sơn La, Lai Châu) thông qua việc định lượng lượng nước thất thoát vào khí quyển từ lưu vực.
  • Tác động Xã hội & Môi trường: Giảm thiểu thiệt hại kinh tế nông nghiệp do hạn hán cục bộ, nâng cao năng lực ứng phó biến đổi khí hậu cho cộng đồng dân tộc thiểu số vùng cao Tây Bắc.

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

  • Nghiên cứu sinh và Học viên cao học: Tiếp cận phương pháp luận lai ghép mô hình viễn thám - khí tượng và kỹ thuật lập trình giải tích không gian trên nền tảng Google Earth Engine.
  • Các nhà khoa học cao cấp: Kế thừa bộ tham số thực nghiệm $(a, b)$ và mô hình toán nhiệt động học hiệu chỉnh DEM để phát triển các mô hình hoàn lưu khí quyển - thủy văn quy mô khu vực.
  • Kỹ sư Thủy lợi và R&D Nông nghiệp: Sử dụng bản đồ phân bố $ET_a$ để thiết kế định mức tưới chính xác cho từng loại cây trồng theo từng tiểu vùng khí hậu.
  • Các nhà hoạch định chính sách: Sử dụng hệ thống bản đồ giám sát tự động để ra quyết định phân bổ ngân sách phòng chống hạn hán và phòng cháy chữa cháy rừng kịp thời.

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 độc đáo nhất là việc tái cấu trúc và địa phương hóa mô hình bức xạ Priestley - Taylor (1972) bằng cách tích hợp trực tiếp năng lượng bức xạ ròng viễn thám ($R_{nd_VT}$ từ SEBAL) và các tham số nhiệt động học ($\lambda, \gamma, \Delta$) phụ thuộc độ cao DEM, thay thế hoàn toàn nhu cầu đo đạc khí tượng mặt đất tức thời cho khu vực đồi núi nhiệt đới.

2. Đột phá phương pháp luận của luận án khi so sánh với các công trình quốc tế? Khác với mô hình SEBAL nguyên bản (Bastiaanssen et al., 1998) hoặc METRIC (Allen et al., 2007) đòi hỏi quy trình chọn điểm ảnh khô/ướt thủ công dễ sai lệch ở địa hình phức tạp, và khác với nghiên cứu của Bruno et al. (2018) hay Ghaderi et al. (2020) phụ thuộc vào trạm khí tượng mặt đất, phương pháp luận của luận án loại bỏ hoàn toàn bước giải lặp thông lượng nhiệt cảm ứng $H$ bằng cách chuyển hóa trực tiếp năng lượng bức xạ ròng sang $ET_a$ qua phương trình Priestley - Taylor hiệu chuẩn, tự động hóa 100% trên Google Earth Engine.

3. Phát hiện bất ngờ và có ý nghĩa thống kê rõ nét nhất từ dữ liệu thực nghiệm? Phát hiện bất ngờ nhất là giá trị bức xạ ròng ngày $R_{nd_VT}$ trích xuất từ vệ tinh Landsat 8 có sự tương thích gần như tuyệt đối với bức xạ tính toán theo tiêu chuẩn FAO-56 tại các trạm mặt đất, với sai số chênh lệch chỉ từ $0,05 - 0,42\text{ MJ/m}^2\text{/ngày}$, dẫn đến sai số ước tính $ET_a$ toàn mạng lưới trạm luôn duy trì ở mức dưới $0,21\text{ mm/ngày}$ (sai số tương đối $< 4,5%$).

4. Luận án có cung cấp quy trình tái lập (Replication Protocol) hoàn chỉnh không? Có. Toàn bộ quy trình giải tích từ tính chuyển $DN \to L_\lambda \to \rho_\lambda$, tính toán các chỉ số $\alpha, NDVI, SAVI, LAI, \varepsilon_0, T_s$, mô hình hóa $R_{ni} \to R_{nd_VT}$, tích hợp độ cao SRTM DEM để tính $\gamma(z), \lambda(T_s), \Delta(T_s)$ và mã nguồn trên Google Earth Engine đều được công bố chi tiết với các sơ đồ khối giải thuật rõ ràng.

5. Định hướng nghiên cứu 10 năm tiếp theo từ nền tảng luận án? Xây dựng hệ thống giám sát chu trình nước tích hợp đa tầng viễn thám (Multi-scale Earth Observation Digital Twin) kết hợp trí tuệ nhân tạo (AI/Deep Learning) để dự báo nhu cầu tưới tiêu và rủi ro sinh thái thời gian thực trên toàn bộ lưu vực sông Hồng - sông Đà.

Kết luận

Công trình luận án tiến sĩ của tác giả Lê Hùng Chiến đã tạo nên những đóng góp khoa học và thực tiễn xuất sắc:

  1. Hoàn thiện khung lý thuyết: Thiết lập thành công mô hình toán học tích hợp giữa mô hình cân bằng bức xạ viễn thám SEBAL và mô hình bức xạ Priestley - Taylor hiệu chỉnh theo mô hình số độ cao DEM cho vùng nhiệt đới đồi núi.
  2. Xác lập hệ số thực nghiệm địa phương: Đề xuất chính xác bộ hệ số $(a, b)$ tối ưu cho phương trình Priestley - Taylor đặc trưng cho điều kiện khí hậu, địa hình và thảm phủ tỉnh Hòa Bình nói riêng và khu vực Tây Bắc Việt Nam nói chung.
  3. Chứng minh tính khả thi thay thế dữ liệu mặt đất: Khẳng định các đại lượng nhiệt động học và bức xạ trích xuất từ dữ liệu ảnh Landsat 8 và SRTM DEM hoàn toàn có thể thay thế chuỗi số liệu đo đạc tại các trạm khí tượng thủy văn mặt đất với sai số $< 4,5%$.
  4. Tự động hóa công nghệ trên nền tảng điện toán đám mây: Phát triển chương trình giám sát bốc thoát hơi nước tự động trên nền Google Earth Engine, mở ra khả năng trích xuất dữ liệu không - thời gian quy mô lớn với chi phí vận hành tối thiểu.
  5. Thành lập bản đồ phân bố chuyên đề độ chính xác cao: Xây dựng thành công hệ thống bản đồ phân bố lượng bốc thoát hơi nước thực tế $ET_a$ tỉnh Hòa Bình tại các thời điểm điển hình, phản ánh trung thực mối tương quan sinh thái giữa $ET_a$ với từng trạng thái rừng và thảm phủ thực vật.
  6. Mở ra các hướng nghiên cứu liên ngành mới: Đặt nền móng vững chắc cho việc ứng dụng công nghệ viễn thám tầm cao trong việc giám sát biến động tài nguyên nước, cảnh báo hạn hán và quản lý bền vững các hệ sinh thái rừng nhiệt đới trước tác động của biến đổi khí hậu toàn cầu.