Tổng quan về luận án

Bối cảnh biến đổi khí hậu toàn cầu đang làm thay đổi căn bản tần suất, cường độ và quy luật diễn biến của các cực đoan khí tượng thủy văn, đặc biệt là các tổ hợp thiên tai liên hoàn vùng duyên hải. Tại dải ven biển miền Trung Việt Nam, các hiện tượng thời tiết cực đoan như bão mạnh kết hợp mưa lớn cục bộ và lũ lụt sau bão đã gây ra những tổn thất thảm khốc về người, tài sản và cơ sở hạ tầng. Luận án tiến sĩ chuyên ngành Biến đổi khí hậu (Mã số: 9440221) của nghiên cứu sinh Trần Thanh Thủy, dưới sự hướng dẫn khoa học của GS. Trần Thục tại Viện Khoa học Khí tượng Thủy văn và Biến đổi khí hậu (Bộ Tài nguyên và Môi trường, 2021) mang tên: "Nghiên cứu đánh giá rủi ro đa thiên tai đối với các tỉnh ven biển Trung Trung Bộ". Đây là công trình học thuật tiên phong tại Việt Nam chuyển dịch mô hình nghiên cứu từ đánh giá rủi ro đơn thiên tai cục bộ sang đánh giá rủi ro đa thiên tai (Multi-Hazard Risk Assessment) có tính đến tương tác chuỗi, tác động đồng thời và nối tiếp trong bối cảnh biến đổi khí hậu.

Khoảng trống nghiên cứu (Research gap) cốt lõi được xác định rõ ràng: Phần lớn các nghiên cứu truyền thống tại Việt Nam (như Nguyễn Mai Đăng, 2010; Nguyễn Xuân Hiển, 2013; Cấn Thu Văn và Nguyễn Thanh Sơn, 2016) chủ yếu tiếp cận theo hướng đơn thiên tai (single-hazard approach), đánh giá độc lập từng hiểm họa như bão, ngập lụt, hoặc nước dâng mà bỏ qua hiệu ứng cộng hưởng, khuếch đại phi tuyến tính khi các hiểm họa xuất hiện đồng thời hoặc nối tiếp nhau (cascading and compound hazards). Đồng thời, các nghiên cứu quốc tế về đa thiên tai (như Liu et al., 2015; Gallina et al., 2016) chưa được bản địa hóa và hoàn thiện cấu trúc toán học để tích hợp đầy đủ các hàm phân bố xác suất đồng thời (Copula đa biến) với bộ chỉ thị tổn thương đa chiều (kinh tế, xã hội, môi trường, hạ tầng) ở quy mô cấp quận/huyện ven biển.

Luận án tập trung giải quyết 3 câu hỏi nghiên cứu (Research Questions - RQ) và kiểm chứng 3 luận điểm/giả thuyết khoa học (Hypotheses - H):

  1. RQ1: Phương pháp luận và mô hình toán học nào có thể đánh giá định lượng xác thực rủi ro đa thiên tai xảy ra đồng thời hoặc nối tiếp ở dải ven biển?
    • H1: Rủi ro đa thiên tai được xác định chuẩn xác thông qua tổ hợp phi tuyến giữa xác suất đồng thời của đa hiểm họa (tính bằng hàm Copula 3 biến) và sự gia tăng tính dễ bị tổn thương tổng hợp trước các hiểm họa liên hoàn.
  2. RQ2: Mức độ rủi ro tại vùng ven biển Trung Trung Bộ gia tăng định lượng ra sao khi bão, mưa trong bão và mưa sau bão xảy ra đồng thời hoặc nối tiếp so với từng đơn thiên tai riêng lẻ?
    • H2: Sự tương tác đồng thời và nối tiếp giữa gió mạnh trong bão và các đợt mưa lớn làm gia tăng đáng kể mức độ rủi ro tổng hợp so với tổng đại số của các rủi ro đơn lẻ.
  3. RQ3: Biến đổi khí hậu theo các kịch bản phát thải tác động như thế nào đến sự gia tăng đa hiểm họa và rủi ro đa thiên tai trong tương lai?
    • H3: Biến đổi khí hậu theo kịch bản nồng độ khí nhà kính (RCP4.5 và RCP8.5) làm gia tăng rõ rệt lượng mưa cực đoan một ngày lớn nhất và tần suất xuất hiện các tổ hợp đa hiểm họa cực hạn.

Khung lý thuyết của nghiên cứu được tích hợp chặt chẽ từ Khung quản lý rủi ro thiên tai của UNISDR (2009, 2017) và Báo cáo đặc biệt SREX của IPCC (2012, 2014), cấu trúc rủi ro thành hàm số của ba thành phần cơ bản: Hiểm họa ($Hazard - H$), Mức độ phơi bày ($Exposure - E$) và Tính dễ bị tổn thương ($Vulnerability - V$), trong đó $V$ được phân rã thành Độ nhạy cảm ($Sensitivity - S$) và Năng lực chống chịu/Nguồn lực ($Capacity - C$).

Phạm vi nghiên cứu bao quát toàn bộ 6 tỉnh/thành phố ven biển Trung Trung Bộ gồm: Quảng Bình, Quảng Trị, Thừa Thiên - Huế, Đà Nẵng, Quảng Nam và Quảng Ngãi. Chuỗi số liệu khí tượng thủy văn thực đo kéo dài 58-60 năm (1961-2020 với 91 cơn bão/áp thấp nhiệt đới; 1961-2018 cho quan trắc gió và mưa chi tiết), kết hợp với số liệu thống kê kinh tế - xã hội 2018-2019 tại tất cả các đơn vị hành chính cấp quận/huyện, tạo nên một cơ sở dữ liệu thực nghiệm quy mô và vững chắc.

Literature Review và Positioning

Tổng quan y văn quốc tế và trong nước chỉ ra sự chuyển dịch căn bản của hệ hình nghiên cứu rủi ro thiên tai qua các giai đoạn:

[Tiếp cận Đơn thiên tai (Single-hazard)] 
         │ 
         ▼ (Bỏ qua tương tác chuỗi & cộng hưởng)
[Mô hình Thiệt hại/Chỉ số Đa hiểm họa sơ khai (HAZUS-MH, CAPRA, Barua et al.)]
         │ 
         ▼ (Thiếu hàm xác suất đồng thời & TDBTT động)
[Tiếp cận Rủi ro Đa thiên tai Tích hợp (Luận án: Copula 3 biến + Ma trận TDBTT + IPCC Framework)]

Dòng nghiên cứu đầu tiên tập trung vào các công cụ định lượng thiệt hại đơn thiên tai quy mô lớn. Ban Quản lý Khẩn cấp Liên bang Hoa Kỳ (FEMA) phát triển công cụ HAZUS-MH (FEMA, 2010; Gallina et al., 2016), cho phép ước tính thiệt hại vật lý và kinh tế do bão, lũ lụt, động đất. Tuy nhiên, HAZUS-MH chỉ dừng lại ở việc so sánh thiệt hại song song giữa các đơn thiên tai, hoàn toàn không đánh giá được sự cộng hưởng rủi ro khi các hiểm họa kích hoạt lẫn nhau. Tương tự, dự án RiskScape tại Hawke's Bay, New Zealand (GNS Science & NIWA; Gallina et al., 2016) xây dựng đường cong thiệt hại vượt ngưỡng cho động đất, lũ lụt và bão biển, chỉ ra rằng với ngưỡng thiệt hại trên 10.000 USD/km², bão có tỷ lệ diện tích ảnh hưởng cao nhất (>80%), nhưng mô hình này chỉ dựa vào hàm thiệt hại quá khứ mà bỏ qua mức độ phơi bày và cấu trúc tính dễ bị tổn thương của xã hội.

Dòng nghiên cứu thứ hai hướng đến việc tích hợp đa hiểm họa thông qua phần mềm CAPRA (CEPREDENAC, IDB, World Bank; Gallina et al., 2016) tại Trung Mỹ. CAPRA phân tích hiểm họa chuỗi bằng cách tính tổng thiệt hại các hiểm họa thứ cấp (như sạt lở do bão, sóng ngập). Dù vậy, CAPRA đòi hỏi cơ sở dữ liệu thiệt hại chi tiết mà các nước đang phát triển không có, đồng thời chưa tích hợp kịch bản biến đổi khí hậu. Nghiên cứu của Rafiq et al. (2012) tại Pakistan và Barua et al. (2016) tại Bangladesh áp dụng phương pháp cho điểm chỉ số đa hiểm họa cấp huyện, nhưng lại mắc hạn chế khi cộng gộp cơ học các hiểm họa không thể xảy ra đồng thời (như hạn hán và lũ lụt) và bỏ qua thành phần tính dễ bị tổn thương.

Dòng nghiên cứu thứ ba – phát triển các phương pháp ma trận và xác suất hiện đại – được khởi xướng bởi Liu et al. (2015, 2017) và Gallina (2015). Liu et al. (2015) đề xuất khung phân tích 3 cấp (Định tính, Bán định lượng bằng ma trận trọng số 0-3, và Định lượng bằng Mạng Bayes). Gallina (2015) ứng dụng phân tích đa tiêu chuẩn (MCDA) để đánh giá rủi ro đa thiên tai vùng ven biển dưới tác động của biến đổi khí hậu. Tuy nhiên, Gallina (2015) coi các đối tượng phơi bày có trọng số bằng nhau và tính xác suất đa hiểm họa bằng tổ hợp biến cố độc lập, không phản ánh đúng cấu trúc phụ thuộc phi tuyến giữa các cực đoan thời tiết. Nhằm giải quyết hạn chế này, Xu et al. (2018) và các tác giả quốc tế (Kew et al., 2019) đã đưa lý thuyết Copula vào phân tích phụ thuộc đa biến của bão, triều cường và mưa lớn.

Tại Việt Nam, các nghiên cứu rủi ro ven biển (Trần Cảnh Dương, 2018; Nguyễn Kim Anh et al., 2019; WB5 Project, 2017; Viện KTTVBĐKH, 2020) chủ yếu phân vùng nguy cơ bão đơn lẻ hoặc xây dựng chỉ số tổn thương tĩnh. Chưa có công trình nào tích hợp mô hình Copula 3 biến thực nghiệm với ma trận gia tăng tính dễ bị tổn thương để định lượng rủi ro chuỗi "Gió bão – Mưa trong bão – Mưa sau bão". Luận án của Trần Thanh Thủy đã định vị chính xác khoảng trống này và giải quyết triệt để trên bình diện lý thuyết lẫn thực tiễn.

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

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

Luận án đã đóng góp mở rộng khung lý thuyết quản lý rủi ro thiên tai của UNISDR/IPCC qua việc cấu tạo lại bản chất toán học và thực nghiệm của hàm rủi ro đa thiên tai:

  1. Chuỗi hiểm họa liên hoàn (Cascading Hazards): Luận án chỉ ra rằng hiểm họa đa thiên tai không phải là tổng số rời rạc của các đơn hiểm họa mà là hàm phụ thuộc phi tuyến. Sự xuất hiện của gió mạnh trong bão cấp $\ge 8$ ($H_{GM}$) kết hợp với mưa lớn trong bão $\ge 100\text{ mm/ngày}$ ($H_{MTB}$) và mưa lớn sau bão $\ge 100\text{ mm/ngày}$ ($H_{MSB}$) tạo ra một trạng thái cực đoan mới về mặt thủy văn và động lực học ven bờ.
  2. Khái niệm gia tăng Tính dễ bị tổn thương động (Dynamic Vulnerability Amplification): Khung lý thuyết của luận án chứng minh rằng tính dễ bị tổn thương của hệ thống kinh tế - xã hội đối với đa hiểm họa ($V_{MH}$) luôn lớn hơn tính dễ bị tổn thương đối với từng đơn hiểm họa ($V_{SH}$). Khi bão đổ bộ làm sập đổ công trình và hư hại hệ thống thoát nước, độ nhạy cảm ($S$) của khu vực tăng đột biến trong khi nguồn lực ứng phó tại chỗ ($C$) bị suy giảm nghiêm trọng, khiến tác động của đợt mưa lớn nối tiếp sau bão bị khuếch đại gấp nhiều lần.
  3. Mô hình toán rủi ro đa thiên tai tổng quát: $$\text{Risk}{MH} = f(H{MH}, E_{MH}, V_{MH}) = H_{MH} \times E_{MH} \times V_{MH}$$ Trong đó $H_{MH}$ được xác định bằng xác suất vượt ngưỡng đồng thời 3 chiều thông qua hàm Copula; $E_{MH}$ là chỉ số phơi bày tổng hợp chuẩn hóa; và $V_{MH}$ là tính dễ bị tổn thương tương tác có xét ma trận điều chỉnh trọng số.
       ┌─────────────────────────────────────────────────────────┐
       │             RỦI RO ĐA THIÊN TAI (Risk_MH)               │
       └────────────────────────────┬────────────────────────────┘
                                    │
         ┌──────────────────────────┼──────────────────────────┐
         ▼                          ▼                          ▼
 ┌───────────────┐          ┌───────────────┐          ┌───────────────┐
 │ Đa hiểm họa   │          │ Mức phơi bày  │          │ TDBTT Đa TT   │
 │    (H_MH)     │          │    (E_MH)     │          │    (V_MH)     │
 └───────┬───────┘          └───────┬───────┘          └───────┬───────┘
         │                          │                          │
 ┌───────┴───────┐          ┌───────┴───────┐          ┌───────┴───────┐
 │Copula 3 biến: │          │Dân số, GRDP,  │          │V_MH = S_MH/   │
 │- Gió bão      │          │Đất nông nghiệp│          │       C_MH    │
 │- Mưa trong bão│          │Hạ tầng giao   │          │Ma trận tương  │
 │- Mưa sau bão  │          │thông ven biển │          │tác tổn thương │
 └───────────────┘          └───────────────┘          └───────────────┘

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

Khung phân tích của luận án tích hợp ba trụ cột phương pháp luận tiên tiến:

  • Lý thuyết xác suất đa biến Copula: Sử dụng Copula 3 biến Archimedean (Gumbel-Hougaard, Clayton, Frank) để liên kết các phân bố biên cực trị đơn biến phi chuẩn của gió bão và lượng mưa, mô hình hóa cấu trúc phụ thuộc phức tạp mà các phân bố chuẩn truyền thống không thể mô tả.
  • Phương pháp trọng số Iyengar - Sudarshan (1982): Loại bỏ tính chủ quan trong việc gán trọng số cho các biến số kinh tế - xã hội, đảm bảo tính khách quan toán học dựa trên phương sai dữ liệu thực tế.
  • Ma trận tương tác tính dễ bị tổn thương bán định lượng: Kế thừa phương pháp ma trận cấp độ của Liu et al. (2015) để xác lập hệ số khuếch đại tổn thương tương hỗ giữa các loại hình thiên tai diễn ra liên tiếp.

Điều kiện biên (Boundary conditions): Khung phân tích áp dụng đặc hiệu cho vùng dải ven biển nhiệt đới gió mùa chịu tác động trực tiếp của xoáy thuận nhiệt đới kết hợp không khí lạnh và địa hình hẹp, dốc từ Tây sang Đông của dãy Trường Sơn.

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

Thiết kế nghiên cứu

Nghiên cứu đứng trên lập trường nhận thức luận thực chứng phản biện (Critical Realism) kết hợp phương pháp luận đa cấp (Multi-level mixed-method design):

  • Cấp độ vĩ mô (Macro-level): Mô hình hóa hoàn lưu khí quyển, thống kê cực trị khí hậu toàn vùng ven biển Trung Trung Bộ dựa trên chuỗi số liệu 1961-2020 và số liệu tái phân tích ERA5 của Trung tâm Dự báo Thời tiết Hạn vừa Châu Âu (ECMWF).
  • Cấp độ vi mô (Micro-level): Đánh giá chi tiết mức độ phơi bày ($E$) và tính dễ bị tổn thương ($V$) tại cấp đơn vị hành chính quận/huyện (gồm toàn bộ các huyện thuộc 6 tỉnh nghiên cứu) dựa trên dữ liệu điều tra dân số và niên giám thống kê.

Quy trình nghiên cứu rigorous

Quy trình đánh giá rủi ro đa thiên tai được chuẩn hóa qua 5 giai đoạn nghiêm ngặt:

[BƯỚC 1: Xử lý dữ liệu & Kiểm định phân bố biên đơn biến (Gió, Mưa trong bão, Mưa sau bão)]
                                    │
                                    ▼
[BƯỚC 2: Mô phỏng Xác suất Đồng thời Đa hiểm họa bằng Copula 3 biến (Kiểm định AIC, K-S)]
                                    │
                                    ▼
[BƯỚC 3: Xây dựng Bộ chỉ số Phơi bày (E) & Tổn thương (V = S/C) (Trọng số Iyengar-Sudarshan)]
                                    │
                                    ▼
[BƯỚC 4: Nội suy Không gian Kriging & Tính toán Ma trận Gia tăng TDBTT Đa thiên tai]
                                    │
                                    ▼
[BƯỚC 5: Tích hợp Bản đồ Rủi ro Đa thiên tai & Đánh giá Tác động Kịch bản BĐKH (RCP4.5/8.5)]
  • Xác lập ngưỡng thiên tai:
    • Gió mạnh trong bão: Vận tốc gió $V_{max} \ge \text{Cấp 8}$ (tương đương $\ge 17.2\text{ m/s}$).
    • Mưa lớn trong bão ($MTB$) và mưa lớn sau bão ($MSB$): Lượng mưa ngày lớn nhất $Rx1day \ge 50\text{ mm/ngày}$ và ngưỡng cực đoan $\ge 100\text{ mm/ngày}$.
  • Kiểm định phân bố xác suất biên: Dữ liệu gió và mưa tại các trạm quan trắc (Ba Đồn, Đồng Hới, Tuyên Hóa, Khe Sanh, Cồn Cỏ, Đông Hà, A Lưới, Nam Đông, Huế, Đà Nẵng, Tam Kỳ, Trà My, Ba Tơ, Quảng Ngãi, Lý Sơn) được kiểm định qua các hàm phân bố Pearson III, GEV, Gumbel, Weibull. Kiểm định mức độ phù hợp bằng tiêu chuẩn khoảng cách Kolmogorov - Smirnov ($D_{max}$) và Tiêu chuẩn thông tin Akaike (AIC).
  • Thiết lập Copula 3 chiều: Xác suất xuất hiện đồng thời của đa hiểm họa được tính toán: $$P(H_{GM} \ge h_1, H_{MTB} \ge h_2, H_{MSB} \ge h_3) = C(u_1, u_2, u_3)$$ Trong đó $u_1, u_2, u_3$ là các hàm phân bố tích lũy biên. Hàm Copula tối ưu được chọn dựa trên giá trị AIC nhỏ nhất và kiểm định K-S nhiều chiều.

Data và phân tích

  • Chuẩn hóa dữ liệu: Sử dụng phương pháp chuẩn hóa Min-Max đưa các biến số về thang đo không thứ nguyên $[0, 1]$:
    • Đối với chỉ số đồng biến với tổn thương (chỉ số nhạy cảm): $$x_{ij}^* = \frac{x_{ij} - \min(x_j)}{\max(x_j) - \min(x_j)}$$
    • Đối với chỉ số nghịch biến với tổn thương (chỉ số nguồn lực/năng lực thích ứng): $$x_{ij}^* = \frac{\max(x_j) - x_{ij}}{\max(x_j) - \min(x_j)}$$
  • Xác định trọng số khách quan: Phương pháp Iyengar và Sudarshan (1982) giả định trọng số $w_j$ tỉ lệ nghịch với độ lệch chuẩn của chỉ số thành phần đã chuẩn hóa, nhằm đảm bảo không một chỉ số đơn lẻ nào chi phối toàn bộ kết quả phân tích: $$w_j = \frac{c}{\sqrt{\text{Var}(x_j^)}}, \quad \text{với } c = \left[\sum_{j=1}^K \frac{1}{\sqrt{\text{Var}(x_j^)}}\right]^{-1}$$
  • Nội suy không gian: Ứng dụng kỹ thuật Geostatistics Kriging tối ưu trên phần mềm ArcGIS để xây dựng lưới bản đồ liên tục về phân bố xác suất hiểm họa, mức độ phơi bày và rủi ro trên toàn bộ dải duyên hải Trung Trung Bộ.

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

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

  1. Hiệu ứng khuếch đại phi đại số của rủi ro đa thiên tai: Luận án cung cấp bằng chứng thực nghiệm rõ nét rằng chỉ số rủi ro đa thiên tai ($RRĐTT$) tại dải ven biển Trung Trung Bộ cao hơn đáng kể (từ $1.5$ đến $2.8$ lần) so với rủi ro của từng đơn thiên tai riêng lẻ. Khu vực chịu rủi ro đa thiên tai cao nhất tập trung dọc dải đồng bằng ven biển từ Thừa Thiên - Huế, Đà Nẵng, Quảng Nam đến Quảng Ngãi – nơi có sự trùng hợp về không gian giữa vùng đổ bộ của bão mạnh cấp 11-12, lượng mưa tích lũy trong và sau bão vượt $500\text{ mm}$, cùng mật độ dân số và tài sản phơi bày cao.
  2. Gia tăng cực đoan mưa trong bão do biến đổi khí hậu: Lượng mưa một ngày lớn nhất thực đo trong bão giai đoạn 1961-2018 đạt kỷ lục trên $747\text{ mm/ngày}$ được ghi nhận tại trạm Trà Bồng (Quảng Ngãi) trong cơn bão Ketsana năm 2009. Xu thế biến đổi lượng mưa $Rx1day$ trong 58 năm qua tăng mạnh tại hầu hết các trạm, điển hình như Ba Đồn tăng $36.2%/58\text{ năm}$, Trà My tăng $31.4%/58\text{ năm}$, Đồng Hới tăng $30.8%/58\text{ năm}$, Quảng Ngãi tăng $30.1%/58\text{ năm}$.
  3. Hiện tượng nghịch lý nhiệt độ tại Thừa Thiên - Huế: Trong khi toàn dải ven biển Trung Trung Bộ ghi nhận mức tăng nhiệt độ trung bình năm từ $+0.44^\circ\text{C}/61\text{ năm}$ (Đà Nẵng) đến $+1.08^\circ\text{C}/61\text{ năm}$ (Ba Đồn), trạm khí tượng Huế là trạm duy nhất trên cả nước ghi nhận xu thế giảm nhiệt độ trung bình năm ($-0.27^\circ\text{C}/61\text{ năm}$). Luận án đã giải thích hiện tượng này từ góc độ vi khí hậu: sự gia tăng của lượng mưa năm (Huế tăng $25.0%/58\text{ năm}$, A Lưới tăng $29.3%$) và số đợt mưa kéo dài kết hợp với tần suất không khí lạnh tràn về sớm làm gia tăng độ che phủ mây và hiệu ứng làm mát bốc hơi bề mặt.
  4. Dự báo gia tăng rủi ro tương lai theo kịch bản RCP: Theo kịch bản RCP4.5 và RCP8.5, lượng mưa 1 ngày lớn nhất trung bình năm ($Rx1day$) vào giữa và cuối thế kỷ 21 có xu thế tăng lên tới $70%$, tập trung trọng điểm từ Thừa Thiên - Huế đến Quảng Nam. Cường độ bão mạnh và rất mạnh tăng từ $2 \div 11%$, lượng mưa trong bán kính $100\text{ km}$ từ tâm bão tăng khoảng $20%$, đẩy xác suất xuất hiện các đợt đa thiên tai cực hạn từ chu kỳ lặp 100 năm xuống còn 25-30 năm.
┌───────────────────────────────┬───────────────────────────────┬───────────────────────────────┐
│ Cực đoan Lượng mưa (1961-2018)│ Hiện tượng Khí hậu Dị thường  │ Dự báo Tương lai (Kịch bản)   │
├───────────────────────────────┼───────────────────────────────┼───────────────────────────────┤
│ • Trà Bồng (2009): 747 mm/ngày│ • Huế: Nhiệt độ giảm -0.27°C  │ • Rx1day tăng tới 70%         │
│ • Ba Đồn: Tăng +36.2%/58 năm  │   (Duy nhất toàn quốc)        │ • Cường độ bão tăng +2÷11%    │
│ • Trà My: Tăng +31.4%/58 năm  │ • Tăng lượng mưa năm +29.3%   │ • Mưa bán kính 100km tăng +20%│
│ • Đồng Hới: Tăng +30.8%/58 năm│   kết hợp mây dày & lạnh ẩm   │ • Chu kỳ lặp 100 năm -> 25 năm│
└───────────────────────────────┴───────────────────────────────┴───────────────────────────────┘

Implications đa chiều

  • Về mặt lý thuyết: Cung cấp cơ sở học thuật hoàn chỉnh để tích hợp khái niệm "rủi ro tương tác chuỗi" vào các chương trình nghiên cứu biến đổi khí hậu và quản lý thiên tai quốc gia.
  • Về mặt phương pháp luận: Thiết lập quy trình 5 bước chuẩn hóa, sử dụng Copula và ma trận trọng số khách quan có thể nhân rộng chuyển giao cho tất cả 28 tỉnh/thành ven biển Việt Nam cũng như các vùng duyên hải Đông Nam Á.
  • Về mặt quản lý thực tiễn và chính sách:
    • Chuyển đổi căn bản quy hoạch đê điều, thoát lũ đô thị tại Đà Nẵng, Tam Kỳ, Đồng Hới: Tiêu chuẩn thiết kế công trình ven biển không thể chỉ căn cứ vào cấp gió bão đơn thuần mà bắt buộc phải tính đến kịch bản cộng hưởng "Gió bão cấp 12 + Triều dâng + Mưa sau bão $>300\text{ mm}$".
    • Tích hợp bản đồ phân vùng rủi ro đa thiên tai vào Quy hoạch tỉnh thời kỳ 2021-2030, tầm nhìn 2050; ưu tiên nguồn lực đầu tư cho các huyện có tính dễ bị tổn thương cao như các huyện miền núi nghèo (Nam Giang, Tây Giang, Nam Trà My - nơi có tỷ lệ nước sạch thấp và kết nối giao thông dễ chia cắt) và các vùng bãi ngang ven biển.

Limitations và Future Research

  1. Hạn chế về dữ liệu chuỗi thời gian thiệt hại: Cơ sở dữ liệu thiệt hại thiên tai lịch sử tại cấp xã/phường chưa được phân tách chi tiết theo từng cấu phần nguyên nhân (do gió, do ngập hay do sạt lở), buộc nghiên cứu phải sử dụng phương pháp bán định lượng thông qua bộ chỉ số kinh tế - xã hội thay vì xây dựng các hàm tổn thất vật lý trực tiếp (damage functions).
  2. Độ phân giải mô hình không gian: Việc nội suy Kriging từ mạng lưới trạm quan trắc mặt đất (15 trạm chính) dù được hỗ trợ bởi số liệu tái phân tích ERA5 nhưng vẫn còn hạn chế nhất định ở các vùng địa hình chia cắt phức tạp sâu trong dãy Trường Sơn.
  3. Định hướng nghiên cứu tiếp theo:
    • Ứng dụng mô hình khí hậu động lực khu vực phân giải siêu cao (dưới $3\text{ km}$) kết hợp trí tuệ nhân tạo (Machine Learning/Deep Learning) để dự báo cực ngắn rủi ro đa thiên tai thời gian thực.
    • Tích hợp thêm các hiểm họa liên quan như trượt lở đất, lũ quét và xâm nhập mặn vào hàm Copula đa chiều ($N \ge 4$ biến).
    • Đánh giá tính dễ bị tổn thương động theo thời gian thực (real-time dynamic vulnerability) dựa trên dữ liệu lớn viễn thông và hành vi di chuyển của dân cư.

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

  • Tác động học thuật: Luận án mở ra một hướng nghiên cứu mũi nhọn về rủi ro đa thiên tai tại Việt Nam, đóng góp trực tiếp vào hệ thống giáo trình sau đại học của Viện Khoa học Khí tượng Thủy văn và Biến đổi khí hậu, cung cấp tài liệu tham khảo cốt lõi cho các nghiên cứu sinh chuyên ngành Khí tượng, Thủy văn, Biến đổi khí hậu và Quản lý Môi trường.
  • Chuyển đổi quy hoạch ngành: Kết quả nghiên cứu là đầu vào khoa học quan trọng cho Đề tài khoa học công nghệ cấp quốc gia (Mã số: KC.24/16-20) thuộc Chương trình KC.08/16-20 của Bộ Khoa học và Công nghệ, phục vụ xây dựng công cụ hỗ trợ ra quyết định ứng phó đa thiên tai tại miền Trung.
  • Tác động chính sách cấp tỉnh và quốc gia: Cung cấp cơ sở thực chứng để Ban Chỉ đạo Quốc gia về Phòng, chống thiên tai và UBND các tỉnh Quảng Bình, Quảng Trị, Thừa Thiên - Huế, Đà Nẵng, Quảng Nam, Quảng Ngãi rà soát kịch bản sơ tán dân, quy hoạch các điểm tránh trú bão an toàn (hiện có từ 20 đến 145 điểm/huyện) và định hình chiến lược thích ứng biến đổi khí hậu dài hạn.

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

  • Nghiên cứu sinh và Giới học thuật: Tiếp cận phương pháp luận Copula 3 biến và ma trận tổn thương tương tác, giải quyết khoảng trống học thuật trong phân tích hiểm họa chuỗi.
  • Các nhà hoạch định chính sách và Quản lý đô thị ven biển: Sở Kế hoạch & Đầu tư, Sở Xây dựng, Sở Nông nghiệp & MPTNT các tỉnh miền Trung có được bản đồ số hóa phân vùng rủi ro đa thiên tai chi tiết đến cấp huyện phục vụ phân bổ ngân sách công.
  • Cộng đồng dân cư ven biển: Giảm thiểu thiệt hại về tính mạng và sinh kế thông qua các phương án cảnh báo sớm đa thiên tai chính xác hơn, khắc phục tình trạng bị động khi bão tan nhưng lũ lớn và sạt lở lại bất ngờ ập đến.

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ì?
    Đóng góp đột phá nhất là việc mở rộng khung lý thuyết rủi ro IPCC/UNISDR bằng cách tích hợp cấu trúc phụ thuộc phi tuyến tính của hàm Copula 3 biến với ma trận tương tác tính dễ bị tổn thương động, định lượng hóa thành công sự khuếch đại rủi ro phi tuyến của chuỗi đa thiên tai "Bão - Mưa trong bão - Mưa sau bão".
  2. Cải tiến phương pháp luận so với các nghiên cứu quốc tế trước đây?
    So với mô hình HAZUS-MH (FEMA) chỉ đánh giá các đơn thiên tai song song độc lập và nghiên cứu của Gallina (2015) tính xác suất đa hiểm họa theo biến cố độc lập, luận án đã ứng dụng thành công Copula Archimedean đa chiều để mô hình hóa chính xác hàm phân bố xác suất đồng thời thực nghiệm của các cực đoan khí hậu có tương quan phi chuẩn.
  3. Phát hiện bất ngờ nhất từ số liệu thực nghiệm là gì?
    Đó là nghịch lý giảm nhiệt độ trung bình năm tại trạm khí tượng Huế ($-0.27^\circ\text{C}/61\text{ năm}$) ngược chiều hoàn toàn với xu thế nóng lên toàn cầu và xu thế tăng nhiệt của toàn miền Trung ($+0.44 \div +1.08^\circ\text{C}$), do tác động của sự gia tăng lượng mưa ($+25.0%$) và tần suất mây mù, ẩm độ cao kéo dài.
  4. Quy trình nghiên cứu có khả năng nhân rộng (Replication protocol) không?
    Hoàn toàn khả thi. Toàn bộ quy trình 5 bước – từ chuẩn hóa Min-Max, tính trọng số Iyengar-Sudarshan, kiểm định AIC/K-S cho Copula đến nội suy Kriging trên GIS – đã được xây dựng thành một khung thuật toán chuẩn mực, sẵn sàng áp dụng cho bất kỳ vùng lãnh thổ ven biển nào có chuỗi số liệu quan trắc.
  5. Chương trình nghiên cứu 10 năm tiếp theo được định hình ra sao?
    Trọng tâm 10 năm tới là phát triển hệ thống giám sát và cảnh báo sớm rủi ro đa thiên tai thời gian thực (Real-time Multi-hazard Early Warning System - MHEWS) dựa trên trí tuệ nhân tạo, đồng hóa dữ liệu viễn thám độ phân giải cao và mô hình ngập lụt 2 chiều (2D Hydrodynamic modeling).

Kết luận

Công trình nghiên cứu tiến sĩ của Trần Thanh Thủy đã hoàn thành xuất sắc các mục tiêu khoa học đề ra, tạo nên những đóng góp nền tảng mang tính bước ngoặt:

  1. Xác lập hoàn thiện hệ phương pháp luận đánh giá định lượng rủi ro đa thiên tai xảy ra đồng thời hoặc nối tiếp đặc thù cho dải ven biển nhiệt đới.
  2. Ứng dụng thành công lý thuyết Copula 3 biến trong việc xác định xác suất xuất hiện đồng thời của tổ hợp cực đoan gió bão, mưa trong bão và mưa sau bão.
  3. Xây dựng bộ chỉ số phơi bày và tính dễ bị tổn thương đa chiều (kinh tế, xã hội, môi trường, cơ sở hạ tầng) cấp quận/huyện bằng phương pháp trọng số khách quan Iyengar - Sudarshan.
  4. Cung cấp bằng chứng thực nghiệm khẳng định rủi ro đa thiên tai làm tăng mức độ tổn thất gấp nhiều lần so với rủi ro đơn thiên tai tại 6 tỉnh ven biển Trung Trung Bộ.
  5. Lượng hóa xu thế biến đổi của mưa cực đoan trong quá khứ ($Rx1day$ tăng trên $30%$) và dự tính tương lai theo kịch bản RCP ($Rx1day$ tăng tới $70%$).
  6. Mở ra 3 nhánh nghiên cứu học thuật mới: Mô hình hóa hiểm họa chuỗi thời gian thực, TDBTT động đô thị ven biển, và Tích hợp rủi ro đa thiên tai vào quy hoạch không gian thích ứng biến đổi khí hậu bền vững.