Tổng quan về luận án

Biến đổi khí hậu (BĐKH) toàn cầu đang làm thay đổi căn bản quy luật thủy văn và gia tăng tần suất các hiện tượng thời tiết cực đoan. Tại Việt Nam, một trong những quốc gia chịu tổn thất nặng nề nhất, tác động của BĐKH thể hiện rõ rệt qua sự tàn phá của thiên tai: "chỉ tính trong 15 năm trở lại đây, các loại thiên tai như: bão, lũ, lũ quét, sạt lở đất, úng ngập, hạn hán, xâm nhập mặn và các thiên tai khác đã làm chết và mất tích hơn 10.711 người, thiệt hại tài sản ước tính chiếm khoảng 1,5% GDP/năm" (Ủy ban Quốc gia về phòng chống thiên tai, 2017). Trong bức tranh chung đó, dải đất Nam Trung Bộ—trải dài từ Đà Nẵng đến Bình Thuận—là khu vực xung yếu với sườn Đông Trường Sơn dốc đứng, sông suối ngắn và hẹp, thường xuyên hứng chịu các trận lũ lịch sử như đợt mưa lũ tháng 11/1999, năm 2007 (với lượng mưa 600–800 mm) và bão số 12 (Damrey) năm 2017 gây ngập sâu 0,5–2,0 m trên diện rộng.

Khoảng trống nghiên cứu then chốt (Research Gap) được chỉ rõ từ thực tiễn: "hầu hết các kết quả nghiên cứu tính toán mưa, lũ có xét đến biến đổi khí hậu đã công bố trước đây, đều được lấy trung bình hóa từ kết quả của các mô hình khí hậu toàn cầu với lưới tính toán theo phạm vi quốc gia hoặc chỉ một số ít các vùng được tính chi tiết... gây ra hệ quả là sai số khá lớn khi tính toán chi tiết cho các vùng khác nhau". Việc thiếu vắng dữ liệu chi tiết hóa cấp trạm khiến các kỹ sư thủy lợi, giao thông gặp bế tắc khi tính toán khẩu độ công trình và kiểm tra an toàn hồ đập cho các lưu vực vừa và nhỏ.

Luận án tiến sĩ kỹ thuật chuyên ngành Thủy văn học (Mã số: 62-44-02-24) của nghiên cứu sinh Lê Thị Hải Yến, dưới sự hướng dẫn khoa học của PGS. TS Ngô Lê Long và TS Trần Thanh Tùng tại Trường Đại học Thủy Lợi (2018), đã giải quyết triệt để khoảng trống này. Nghiên cứu tập trung giải quyết ba câu hỏi trọng tâm:

  1. RQ1: Xu thế biến đổi lượng mưa ngày cực trị ($R_{1\text{day max}}$) tại Nam Trung Bộ diễn biến như thế nào dưới các kịch bản phát thải khí nhà kính mới của IPCC?
  2. RQ2: Làm thế nào để loại bỏ sai số hệ thống (bias) của các mô hình hoàn lưu chung toàn cầu (GCMs) nhằm chi tiết hóa dữ liệu lượng mưa cực trị về từng vị trí trạm đo mặt đất cụ thể?
  3. RQ3: Mức độ gia tăng của đỉnh lũ thiết kế ($Q_{\max}$) tại các lưu vực vừa, nhỏ và lớn dưới tác động của BĐKH là bao nhiêu, và phân vùng nguy cơ biến động dòng chảy lũ được xác lập ra sao?

Luận án kiểm định hai giả thuyết khoa học:

  • H1: Phương pháp chi tiết hóa thống kê (Statistical Downscaling) kết hợp hiệu chỉnh sai số theo hàm phân phối xác suất thực đo cho phép hạ quy mô không gian từ tổ hợp 11 mô hình GCMs về 93 trạm khí tượng với độ tin cậy vượt trội so với phương pháp nội suy trung bình hóa thô.
  • H2: Dưới các kịch bản RCP4.5 và RCP8.5 trong hai giai đoạn tương lai (2040–2069 và 2070–2099), cường độ mưa 1 ngày lớn nhất và đỉnh lũ thiết kế tại Nam Trung Bộ sẽ gia tăng bất đối xứng theo không gian, đòi hỏi phải tái lập bản đồ phân vùng biến động phục vụ tiêu chuẩn thiết kế công trình đa ngành.

Phạm vi nghiên cứu bao quát 8 tỉnh/thành phố Nam Trung Bộ (Đà Nẵng, Quảng Nam, Quảng Ngãi, Bình Định, Phú Yên, Khánh Hòa, Ninh Thuận, Bình Thuận), khai thác chuỗi dữ liệu thực đo tại 93 trạm khí tượng thủy văn kết hợp tổ hợp 11 mô hình GCMs thuộc dự án CMIP5 dưới hai kịch bản nồng độ đại diện RCP4.5 và RCP8.5.

Literature Review và Positioning

Nghiên cứu về biến đổi của mưa lớn và dòng chảy lũ cực trị trên thế giới đã trải qua nhiều giai đoạn phát triển. Ở quy mô toàn cầu, Frich et al. (2002) và Alexander et al. (2006) khi phân tích chuỗi số liệu quan trắc thế kỷ 20 đã chứng minh xu thế tăng vượt trội của các chỉ số mưa cực trị tại nhiều khu vực thuộc Bắc bán cầu. Ở quy mô châu lục, Re & Barros (2009) tại Đông Nam Nam Mỹ và Penalba & Robledo (2009) tại lưu vực sông La Plata ghi nhận sự gia tăng mạnh mẽ của tần suất và cường độ mưa lớn theo mùa. Ngược lại, tại châu Âu, Moberg et al. (2006) và Dankers & Feyen (2008) sử dụng chuỗi liên kết mô hình khí hậu HIRHAM (độ phân giải 12 km) với mô hình thủy văn LISFLOOD đã chỉ ra rằng các trận lũ chu kỳ lặp lại 100 năm tại nhiều lưu vực sông Tây Âu có nguy cơ rút ngắn chu kỳ xuất hiện xuống còn 50 năm vào cuối thế kỷ 21 dưới kịch bản SRES A2.

Trong bức tranh học thuật quốc tế, việc đưa BĐKH vào tiêu chuẩn tính toán công trình được triển khai với nhiều mức độ khác nhau:

  • Trường hợp Na Uy (Norwegian Hydrological Guidelines): Quy định bắt buộc cập nhật tính toán an toàn đập sau mỗi 15–20 năm, trong đó lưu lượng đỉnh lũ thiết kế được định lượng gia tăng thêm từ 20% đến 40% tùy theo từng vùng thủy văn cụ thể.
  • Trường hợp Hàn Quốc (Hyun-Han Kwon et al., 2011): Ứng dụng chuỗi Markov Bayes kết hợp mô hình thời tiết ngẫu nhiên đa biến để mô phỏng trạng thái ẩm đất và tính toán lại lũ thiết kế cho đập Soyang trong bối cảnh BĐKH bất định.
  • Trường hợp Ấn Độ (IS: 5477 và IS: 11223-1985): Thiết lập khung tiêu chuẩn an toàn cho đập dâng và hồ chứa dựa trên các chỉ số mưa lớn nhất khả năng (PMP) và lũ lớn nhất khả năng (PMF).

Tại Việt Nam, các nghiên cứu kinh điển trước đây của Chow & Maidment (1988), Shaw (1964), Lê Đình Thành (1997), Phạm Ngọc Quý et al. (2005), cùng các quy phạm kỹ thuật như QP.TL C-6-77 và TCVN 9845:2013 chủ yếu dựa trên giả thiết chuỗi thủy văn đồng nhất và dừng (stationary assumption). Khi BĐKH được xem xét, các công trình của Lê Mạnh Hùng & Lê Thị Kim Cúc (2011), Vũ Thanh Tâm et al. (2013), Đoàn Thị Nội (2016), La Đức Dũng (2017) và Ngô Lê An et al. (2015) đã có những đóng góp ban đầu nhưng bộc lộ hạn chế lớn: chỉ sử dụng kịch bản phát thải cũ (SRES A1B, B2, A2) hoặc chỉ chạy đơn lẻ một vài mô hình GCM, dẫn đến việc thiếu đánh giá tính bất định giữa các mô hình và không thể chi tiết hóa xuống quy mô lưu vực vừa và nhỏ. Luận án của Lê Thị Hải Yến định vị chính xác ở điểm giao thoa giữa việc khắc phục tính bất định của tổ hợp đa mô hình CMIP5 và bài toán chuyển tiếp quy mô không gian (downscaling) phục vụ thiết kế kỹ thuật công trình.

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

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

Luận án tạo ra bước tiến quan trọng trong việc mở rộng lý thuyết tính toán thủy văn công trình truyền thống (Hydrological Engineering Design Theory) từ trạng thái chuỗi dừng sang trạng thái phi dừng (non-stationary paradigm) dưới tác động của BĐKH. Cụ thể:

  • Mở rộng lý thuyết chuyển đổi mưa – dòng chảy của Chow & Maidment (1988) và lý thuyết đường lũ đơn vị SCS thông qua việc tích hợp biến động phân bố xác suất lượng mưa cực trị tương lai vào hệ số dòng chảy và mô hình cường độ giới hạn.
  • Bổ sung cơ sở lý thuyết cho phương pháp hạ quy mô không gian thống kê dựa trên hàm truyền phi tham số (Quantile Mapping) của Gudmundsson et al. (2012), chuẩn hóa quá trình hiệu chỉnh sai số cho các chuỗi biến lượng mưa ngày cực đoan tại vùng địa hình đón gió phức tạp ven biển nhiệt đới gió mùa.
  • Xây dựng mô hình quan hệ toán học giữa lưu lượng đỉnh lũ tức thời ($Q_{\max}$) và lưu lượng trung bình ngày lớn nhất ($Q_{\text{ngày max}}$) tại các trạm thủy văn chính (Nông Sơn, Thành Mỹ, Củng Sơn, An Khê, Bình Tường), chứng minh mối tương quan phi tuyến chặt chẽ ($R^2 > 0,85$) được duy trì ngay cả khi quy mô mưa ngoại suy tăng cao.

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

Khung phân tích của luận án tích hợp liên hoàn ba trụ cột lý thuyết:

  1. Lý thuyết hoàn lưu khí quyển quy mô lớn (Large-scale Atmospheric Circulation): Khai thác cấu trúc động lực học từ tổ hợp 11 mô hình GCMs (Access1-3, CanESM2, CMCC-CMS, CNRM-CM5, HadGEM2-ES, NorESM1-M, MPI-ESM-LR, CCSM4, GFDL-CM3, NCAR-SST, HadGEM2-SST) theo kịch bản RCP4.5 và RCP8.5 do IPCC (Moss et al., 2010; Wayne, 2013) chuẩn hóa.
  2. Lý thuyết hiệu chỉnh phân phối thống kê (Statistical Bias Correction): Thiết lập thuật toán điều chỉnh sai số trung bình và độ lệch chuẩn giữa dữ liệu mô phỏng của 11 GCMs với chuỗi số liệu thực đo tại 93 trạm quan trắc địa phương.
  3. Lý thuyết mô hình hóa thủy văn bán phân bố (Semi-distributed Hydrological Modeling): Ứng dụng mô hình NAM (Nielsen & Hansen, 1973) kết hợp thuật toán truyền lũ đoạn sông (Muskingum/cân bằng lượng trữ) cho lưu vực lớn, song song với công thức Xokolopsky và phương pháp cường độ giới hạn cho các lưu vực vừa và nhỏ.

Điều kiện biên (Boundary Conditions) của khung phân tích được xác định rõ: áp dụng cho các lưu vực sông có độ dốc lớn, thời gian tập trung lũ nhanh ($< 24\text{ giờ}$), chịu sự chi phối trực tiếp của hoàn lưu gió mùa Đông Bắc kết hợp địa hình sườn Đông Trường Sơn và các xoáy thuận nhiệt đới.

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

Thiết kế nghiên cứu

Luận án áp dụng triết lý thực chứng định lượng (Positivism), kết hợp chặt chẽ giữa thống kê khí hậu đa mô hình và mô phỏng toán thủy văn. Thiết kế nghiên cứu đa cấp (Multi-level Design) được cấu trúc qua 4 giai đoạn logic:

  • Cấp 1 (Global/Macro-level): Trích xuất trường dữ liệu lượng mưa ngày từ lưới tính toán của 11 mô hình GCMs thuộc CMIP5 giai đoạn quá khứ (thời kỳ nền) và tương lai (2006–2100).
  • Cấp 2 (Meso-level): Hạ quy mô không gian và hiệu chỉnh sai số thống kê (Quantile Mapping) về 93 trạm quan trắc mặt đất; áp dụng thuật toán nghịch đảo khoảng cách (Inverse Distance Weighting - IDW) để nội suy trường mưa không gian.
  • Cấp 3 (Micro/Basin-level): Xác định lượng mưa thiết kế 1 ngày lớn nhất ứng với các tần suất $P = 0,1%; 0,5%; 1%; 2%; 5%; 10%$ cho từng lưu vực.
  • Cấp 4 (Hydraulic/Impact-level): Mô phỏng dòng chảy lũ và tính toán đỉnh lũ thiết kế $Q_{\max}$ cho các lưu vực điển hình (Vu Gia – Thu Bồn, Kôn – Hà Thanh, Sông Ba) bằng mô hình NAM và phương pháp giải tích.

Quy trình nghiên cứu rigorous

Quy trình thu thập và xử lý số liệu tuân thủ các chuẩn mực khí tượng thủy văn khắt khe:

  • Mạng lưới trạm: Thu thập chuỗi số liệu quan trắc mưa ngày thực đo liên tục từ sau năm 1975 đến 2014 tại 93 trạm khí tượng trải rộng khắp 8 tỉnh Nam Trung Bộ (Đà Nẵng, Trà My, Kham Đức, Nông Sơn, Quy Nhơn, Hoài Nhơn, An Khê, Sơn Hòa, Củng Sơn, Nha Trang, Phan Thiết...).
  • Kiểm soát chất lượng: Kiểm định tính đồng nhất của chuỗi số liệu thực đo bằng phương pháp đường cong tích lũy kép (Double Mass Curve) và kiểm định phi tham số Mann-Kendall.
  • Triangulation: Luận án thực hiện đối sánh chéo giữa 11 mô hình GCMs độc lập nhằm giảm thiểu tính bất định của cấu trúc tham số mô hình đơn lẻ.
  • Đánh giá độ tin cậy mô hình thủy văn: Mô hình NAM được hiệu chỉnh và kiểm định trên chuỗi dòng chảy thực đo giai đoạn 1977–2010 tại 4 trạm khống chế trọng yếu (Nông Sơn, Thành Mỹ, Bình Tường, Củng Sơn), đạt hệ số Nash-Sutcliffe ($NSE$) từ 0,82 đến 0,91 và sai số tổng lượng dòng chảy ($W%$) nhỏ hơn $\pm 8%$.

Data và phân tích

  • Phần mềm và công cụ ứng dụng: Xử lý dữ liệu khí hậu bằng mã nguồn R và Python; mô hình hóa thủy lực - thủy văn bằng MIKE 11 (Module NAM); tích hợp không gian và xây dựng bản đồ phân vùng bằng ArcGIS.
  • Phân tích thống kê sai số: Tính toán độ lệch chuẩn ($\sigma$) và sai số trung bình tuyệt đối ($MAE$) giữa 11 mô hình GCMs với 93 trạm thực đo thời kỳ nền trước và sau khi hiệu chỉnh, đảm bảo sai số phân bố xác suất lượng mưa cực trị giảm dưới $5%$.
  • Phương trình hồi quy đỉnh lũ: Thiết lập các phương trình tương quan tuyến tính và phi tuyến giữa $Q_{\max}$ và $Q_{\text{ngày max}}$:
    • Trạm Nông Sơn (1977–2010): $Q_{\max} = 1,32 \cdot Q_{\text{ngày max}} + 245,6$ ($R^2 = 0,89$)
    • Trạm Thành Mỹ (1977–2010): $Q_{\max} = 1,41 \cdot Q_{\text{ngày max}} + 112,4$ ($R^2 = 0,92$)
    • Trạm Củng Sơn (1977–1997): $Q_{\max} = 1,28 \cdot Q_{\text{ngày max}} + 304,1$ ($R^2 = 0,87$)

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

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

Luận án mang lại 5 phát hiện mang tính đột phá cho ngành khí tượng thủy văn miền Trung:

  1. Đặc trưng hình thế synop gây mưa cực trị: Xác lập bằng chứng định lượng chứng minh "Trong 25 năm gần đây, khoảng 68% các đợt mưa lớn xảy ra trên khu vực Nam Trung Bộ là do XTNĐ hoạt động riêng lẻ (đơn thuần) hoặc là do XTNĐ kết hợp với các loại hình thế khác như KKL và dải HTNĐ gây ra". Mưa lũ tập trung cao độ vào tháng X và XI (chiếm 40–50% tổng lượng mưa năm và 50–70% số trận lũ lớn nhất năm).
  2. Phân hóa không gian rõ rệt của lượng mưa 1 ngày lớn nhất ($R_{1\text{day max}}$):
    • Tại lưu vực Vu Gia – Thu Bồn và Kôn – Hà Thanh, $R_{1\text{day max}}$ có xu thế tăng mạnh từ 15% đến 35% dưới kịch bản RCP8.5 giai đoạn 2070–2099 so với thời kỳ nền.
    • Ngược lại, khu vực phía Nam (Bình Định đến Bình Thuận) xuất hiện sự phân hóa nghịch chiều: các trạm ven biển và trung du (như An Khê, Ayun Pa, Sơn Hòa) có xu thế giảm nhẹ lượng mưa cực trị theo chuỗi quan trắc lịch sử, nhưng lại gia tăng đột biến tính cực đoan của các trận mưa cục bộ trong tương lai.
  3. Sự gia tăng định lượng của đỉnh lũ thiết kế ($Q_{\max}$):
    • Dưới kịch bản RCP4.5 (2040–2069), đỉnh lũ thiết kế tần suất $P = 1%$ tại Thành Mỹ và Nông Sơn tăng từ 12% đến 18%.
    • Dưới kịch bản RCP8.5 (2070–2099), đỉnh lũ thiết kế $Q_{\max}$ có thể gia tăng từ 25% đến 38%, làm thay đổi hoàn toàn cấp bảo đảm công trình theo các tiêu chuẩn hiện hành.
  4. Thành lập bản đồ phân vùng biến động dòng chảy lũ thiết kế: Lần đầu tiên số hóa và phân chia dải lãnh thổ Nam Trung Bộ thành 4 cấp độ rủi ro biến động lũ: Nguy cơ rất cao ($> 30%$), Nguy cơ cao ($20–30%$), Nguy cơ trung bình ($10–20%$), và Nguy cơ thấp ($< 10%$).
  5. Phát hiện tính phi đối xứng theo thời gian: Giai đoạn 2070–2099 ghi nhận mức độ gia tăng mưa lũ cực đoan cao gấp 1,8–2,3 lần so với giai đoạn giữa thế kỷ (2040–2069), phản ánh quán tính nhiệt và sự tích tụ nồng độ phát thải khí nhà kính toàn cầu.

Implications đa chiều

  • Về mặt lý thuyết: Bổ sung bằng chứng thực nghiệm khẳng định sự cần thiết phải thay thế giả định chuỗi dừng bằng mô hình xác suất thích ứng trong tính toán thủy văn cực trị vùng nhiệt đới gió mùa.
  • Về mặt phương pháp luận: Cung cấp quy trình chuẩn (Protocols) kết hợp downscaling thống kê 11 mô hình GCMs với mô hình thủy văn bán phân bố, có thể chuyển giao áp dụng cho các vùng kinh tế – sinh thái khác như Bắc Trung Bộ hay Tây Nguyên.
  • Về ứng dụng thực tiễn và chính sách:
    • Cung cấp cơ sở khoa học để Bộ Nông nghiệp & PTNT, Bộ Giao thông Vận tải sửa đổi, cập nhật tiêu chuẩn tính toán khẩu độ thoát nước cầu cống và dung tích phòng lũ của các hồ chứa thủy điện/thủy lợi (A Vương, Sông Tranh 2, Đắk Mi 4, Sông Ba Hạ).
    • Định hình lộ trình quản lý rủi ro thiên tai thích ứng cho chính quyền 8 tỉnh Nam Trung Bộ, tích hợp kịch bản lũ BĐKH vào quy hoạch không gian đô thị và đê kè ven sông.

Limitations và Future Research

Luận án thừa nhận 4 giới hạn nghiên cứu khách quan:

  1. Độ phân giải thời gian của dữ liệu đầu vào: Phương pháp chi tiết hóa hiện tại dựa trên số liệu mưa bước thời gian ngày (daily), chưa thể hiện chi tiết diễn biến mưa cường độ lớn thời đoạn ngắn (5 phút, 1 giờ, 3 giờ, 6 giờ)—yếu tố then chốt gây ra lũ quét cho các sườn dốc nhỏ hẹp.
  2. Tính bất biến của thảm phủ bề mặt: Mô hình thủy văn giả định điều kiện sử dụng đất, thảm phủ rừng và hạ tầng đô thị là cố định theo thời gian, chưa liên kết động với các kịch bản phát triển kinh tế – xã hội (SSPs).
  3. Hiện tượng tương tác vùng cửa sông: Nghiên cứu chưa tính toán chi tiết bài toán thủy lực ngập lụt vùng hạ du dưới tác động kép của lũ thượng lưu kết hợp triều cường và nước biển dâng cực đoan do bão.

Chương trình nghiên cứu tương lai (Future Research Agenda) đề xuất 4 định hướng:

  • Ứng dụng mô hình khí hậu khu vực phân giải siêu cao (Convection-permitting Regional Climate Models, độ phân giải $< 3\text{ km}$) để chi tiết hóa lượng mưa thời đoạn ngắn sub-daily (giờ/phút).
  • Tích hợp mô hình biến động sử dụng đất (Cellular Automata/CLUE-s) vào mô hình dòng chảy tràn bề mặt.
  • Thiết lập mô hình nối kết thủy văn - thủy lực 2D mô phỏng lan truyền lũ ngập lụt toàn diện hạ lưu có xét tương tác nước biển dâng.
  • Ứng dụng lý thuyết Copula đa chiều để mô hình hóa cấu trúc tương quan phụ thuộc phi tuyến giữa mưa đỉnh lũ và thời gian duy trì lũ.

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

  • Tác động học thuật: Công trình tạo tiền đề trích dẫn quan trọng cho các nghiên cứu về thủy văn phi dừng và an toàn công trình tại Đông Nam Á. Kết quả luận án được công bố trên các tạp chí khoa học chuyên ngành uy tín trong nước và kỷ yếu hội thảo quốc tế.
  • Chuyển đổi ngành công nghiệp xây dựng & thủy lợi: Cung cấp bảng tra thông số lượng mưa 1 ngày lớn nhất và lưu lượng đỉnh lũ thiết kế tương lai cho các đơn vị tư vấn thiết kế công trình thủy lợi (HEC), tư vấn giao thông (TEDI), giúp tối ưu hóa chi phí đầu tư và ngăn ngừa nguy cơ vỡ đập hoặc xói lở mố cầu.
  • Tác động chính sách và an sinh xã hội: Giúp giảm thiểu thiệt hại hàng nghìn tỷ đồng mỗi năm do thiên tai tại miền Trung thông qua việc nâng cao dung tích phòng lũ an toàn cho các công trình hạ tầng trọng yếu, bảo vệ sinh mạng và tài sản cho hơn 9 triệu người dân vùng Nam Trung Bộ.

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

  • Nghiên cứu sinh & Giới học thuật: Tiếp cận quy trình chi tiết hóa thống kê tổ hợp 11 mô hình GCMs và phương pháp luận tính toán thủy văn phi dừng.
  • Kỹ sư thiết kế công trình Thủy lợi - Thủy điện - Giao thông: Khai thác trực tiếp số liệu mưa 1 ngày lớn nhất và đỉnh lũ thiết kế có xét BĐKH để tính toán khẩu độ công trình chính xác, tránh hiện tượng thiếu hụt dung tích thoát lũ.
  • Cơ quan Quản lý Quy hoạch và Vận hành Hồ chứa (EVN, Cục Thủy lợi): Ứng dụng quy trình điều tiết lũ liên hồ chứa linh hoạt trên lưu vực Vu Gia – Thu Bồn và sông Ba.
  • Nhà hoạch định chính sách cấp tỉnh và quốc gia: Có căn cứ khoa học định lượng để xây dựng Kế hoạch hành động ứng phó với BĐKH giai đoạn 2021–2030 tầm nhìn 2050 cho 8 tỉnh 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ì? Trả lời: Luận án đã mở rộng lý thuyết xác suất thống kê thủy văn truyền thống sang môi trường biến đổi khí hậu phi dừng thông qua việc lượng hóa tác động của tổ hợp 11 mô hình khí hậu toàn cầu CMIP5 lên phân bố xác suất lượng mưa cực trị $R_{1\text{day max}}$ cấp trạm và đỉnh lũ thiết kế $Q_{\max}$, lấp đầy khoảng trống thiếu hụt dữ liệu quy mô nhỏ cho miền Trung Việt Nam.

  2. Cải tiến phương pháp luận nào tạo nên sự vượt trội so với các nghiên cứu trước đây? Trả lời: Khác với các nghiên cứu trước đây vốn chỉ sử dụng một vài kịch bản trung bình hóa thô (như kịch bản Bộ TN&MT 2012 hay mô hình đơn lẻ HadGEM2 trong Ngô Lê An et al., 2015), luận án đã thiết lập quy trình tích hợp 11 mô hình GCMs với phương pháp hiệu chỉnh sai số Quantile Mapping tại 93 trạm thực đo, kết hợp mô hình bán phân bố NAM và phân tích GIS không gian.

  3. Phát hiện nào gây bất ngờ nhất về mặt số liệu thực nghiệm? Trả lời: Sự phân hóa không đồng nhất của lượng mưa cực đoan: trong khi vùng núi phía Bắc (Quảng Nam, Đà Nẵng) có xu thế mưa 1 ngày lớn nhất tăng vọt đến $+35%$, thì nhiều trạm vùng cao nguyên và ven biển phía Nam (như An Khê, Ayun Pa, Mang Yang) chuỗi số liệu thực đo lịch sử lại thể hiện xu thế giảm (từ $-0,18$ đến $-1,39\text{ mm/năm}$), nhưng dưới tác động tương lai của RCP8.5 tính bất định và độ biến thiên cực trị lại tăng vọt.

  4. Luận án có cung cấp quy trình nhân bản (Replication Protocol) rõ ràng không? Trả lời: Hoàn toàn có. Quy trình 4 bước gồm: Trích xuất lưới GCM $\rightarrow$ Hiệu chỉnh sai số hàm phân phối xác suất trạm $\rightarrow$ Mô hình hóa mưa – dòng chảy $\rightarrow$ Phân vùng GIS được mô tả chi tiết với đầy đủ thông số tham số hóa của 93 trạm và 4 lưu vực khống chế.

  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: Nghiên cứu định hướng phát triển hệ thống cảnh báo sớm lũ quét – sạt lở đất thời gian thực tích hợp dự báo thời tiết phân giải siêu cao (Radar/WRF), kết hợp chuỗi dữ liệu mưa sub-hourly và mô hình ngập lụt số trị 2D liên kết lưu vực - cửa sông.

Kết luận

  1. Luận án đã giải quyết thành công bài toán chi tiết hóa thông tin khí hậu toàn cầu từ 11 mô hình GCMs về 93 trạm quan trắc mặt đất, cung cấp chuỗi dữ liệu lượng mưa 1 ngày lớn nhất có xét BĐKH với độ chính xác cao.
  2. Xác lập bức tranh định lượng về sự gia tăng của đỉnh lũ thiết kế ($Q_{\max}$) tại các lưu vực trọng điểm Nam Trung Bộ với mức tăng từ 12% đến 38% tùy theo kịch bản (RCP4.5, RCP8.5) và giai đoạn (2040–2069, 2070–2099).
  3. Hoàn thành hệ thống bản đồ phân vùng biến động dòng chảy lũ thiết kế tỷ lệ chi tiết cho toàn bộ 8 tỉnh thành phố vùng Nam Trung Bộ.
  4. Mở ra 3 hướng nghiên cứu học thuật mới: mô hình hóa mưa thời đoạn ngắn sub-daily, tương tác ngập lụt lũ – triều – nước biển dâng cửa sông, và thủy văn phi dừng đa biến.
  5. Cung cấp luận cứ khoa học vững chắc phục vụ rà soát, điều chỉnh các quy chuẩn, tiêu chuẩn kỹ thuật thiết kế công trình thủy lợi, giao thông, hạ tầng đô thị ứng phó với BĐKH tại Việt Nam.