Tổng quan nghiên cứu

Mô phỏng khí hậu khu vực đóng vai trò then chốt trong việc cung cấp thông tin dự báo độ phân giải cao phục vụ quy hoạch kinh tế xã hội và ứng phó biến đổi khí hậu tại Việt Nam. Nghiên cứu thực hiện đánh giá toàn diện khả năng mô phỏng hai yếu tố cốt lõi là nhiệt độ không khí bề mặt mực 2 mét và lượng mưa trong chuỗi thời gian liên tục 22 năm, từ năm 1986 đến năm 2007. Vấn đề khoa học trọng tâm xuất phát từ việc các mô hình hoàn lưu chung toàn cầu thường bỏ sót các quá trình quy mô địa phương do giới hạn lưới thô, trong khi mô hình thủy tĩnh truyền thống bộc lộ nhiều sai số tại khu vực nhiệt đới gió mùa phức tạp.

Mục tiêu cụ thể của công trình là so sánh định lượng kỹ năng mô phỏng giữa mô hình khí hậu khu vực phi thủy tĩnh NHRCM và mô hình thủy tĩnh RegCM phiên bản 4.2 trên lưới tính có độ phân giải ngang 20 km với 40 mực thẳng đứng. Phạm vi không gian bao phủ toàn bộ lãnh thổ Việt Nam và Biển Đông, tâm miền tính đặt tại 110 độ Đông và 15 độ Bắc, tích hợp dữ liệu kiểm chứng từ 58 trạm quan trắc synop cùng bộ dữ liệu lưới Aphrodite độ phân giải 0.25 độ kinh vĩ.

Ý nghĩa thực tiễn của nghiên cứu được xác lập thông qua việc chứng minh mô hình phi thủy tĩnh NHRCM giúp giảm thiểu sai số trung bình nhiệt độ xuống dưới mức 1 độ C đến 2 độ C, đồng thời hạn chế mức thiếu hụt lượng mưa trung bình năm từ 800 mm ở mô hình RegCM xuống chỉ còn khoảng 150 mm so với quan trắc thực tế, tạo tiền đề vững chắc cho việc ứng dụng công nghệ mô hình phi thủy tĩnh trong bài toán dự tính khí hậu tương lai tại nước ta.

Cơ sở lý thuyết và phương pháp nghiên cứu

Khung lý thuyết áp dụng

Nghiên cứu ứng dụng lý thuyết động lực học khí quyển phi thủy tĩnh của mô hình NHRCM kết hợp với lý thuyết cân bằng thủy tĩnh truyền thống của mô hình RegCM nhằm làm rõ sự khác biệt trong việc giải đoán các hệ thống thời tiết nhiệt đới. Trong mô hình NHRCM, các phương trình chuyển động không áp dụng xấp xỉ thủy tĩnh, cho phép mô phỏng chính xác các hoàn lưu thẳng đứng quy mô nhỏ và mây đối lưu sâu. Hệ thống vật lý của NHRCM tích hợp sơ đồ tham số hóa đối lưu Kain-Fritsch kết hợp mô hình sinh quyển đơn giản SiB để mô tả tương tác mặt đệm đất liền và khí quyển, đồng thời áp dụng kỹ thuật ghép phổ Spectral Coupling tại các biên ngoài.

Ngược lại, mô hình RegCM phiên bản 4 vận hành trên hệ tọa độ áp suất chuẩn hóa sigma, sử dụng sơ đồ tham số hóa đối lưu hỗn hợp Grell-Emanuel, trong đó sơ đồ Grell được kích hoạt trên đất liền và sơ đồ Emanuel áp dụng trên đại dương theo khuyến nghị từ Trung tâm Vật lý Lý thuyết Quốc tế. Ba khái niệm nền tảng được vận dụng xuyên suốt bao gồm phương pháp hạ quy mô động lực, sai số hệ thống trong mô phỏng khí hậu dài hạn, và phân bố xác suất Gamma đối với chuỗi số liệu lượng mưa ngày cực đoan.

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

Nguồn dữ liệu đầu vào làm điều kiện biên và điều kiện ban đầu cho cả hai mô hình là bộ số liệu tái phân tích khí quyển toàn cầu JRA-25 của Cơ quan Khí tượng Nhật Bản với độ phân giải phổ T106 gồm 40 mực thẳng đứng, cập nhật định kỳ 6 giờ một lần, kết hợp trường nhiệt độ mặt nước biển phân tích MGDSST. Dữ liệu thực nghiệm kiểm chứng sử dụng chuỗi quan trắc liên tục 22 năm giai đoạn 1986-2007 từ 58 trạm synop đại diện phân bố trên 4 khu vực địa lý theo vĩ độ gồm Bắc Bộ, Bắc Trung Bộ, Nam Trung Bộ và Nam Bộ. Phương pháp chọn mẫu theo vĩ độ được thiết kế có chủ đích nhằm cô lập và đánh giá tác động của hiệu ứng phi thủy tĩnh khi tiến dần về khu vực cận xích đạo.

Về phương pháp phân tích, nghiên cứu lựa chọn chỉ số sai số trung bình ME thay vì chỉ số RMSE. Lý do khoa học là RMSE thường trộn lẫn moment bậc một và moment bậc hai, làm mất tính đặc trưng của sai số hệ thống vốn rất lớn trong các bài toán khí hậu dài hạn. Nhiệt độ không khí mực 2 mét được hiệu chỉnh theo độ cao địa hình thực tế với số gia tiêu chuẩn 0.65 độ C trên 100 mét độ cao. Đối với trường lượng mưa có phân bố dạng Gamma lệch phải, phương pháp đánh giá thông qua hàm logarit và tỷ số phần trăm được áp dụng trên biểu đồ histogram thang đo logarit với độ dài khoảng chia 5 mm cho mưa và 1 độ C cho nhiệt độ.

Kết quả nghiên cứu và thảo luận

Những phát hiện chính

Mô hình NHRCM thể hiện ưu thế vượt trội trong việc tái tạo trường nhiệt độ và lượng mưa so với mô hình RegCM trên toàn bộ lãnh thổ Việt Nam trong suốt giai đoạn 1986-2007:

  • Độ chính xác trường nhiệt độ: Mô hình RegCM liên tục mô phỏng nhiệt độ thiên thấp so với thực tế, độ lệch âm đạt từ 1 độ C đến 3 độ C trên phần lớn diện tích và lên tới 4 độ C tại các vùng núi cao Tây Bắc và Tây Nguyên. Ngược lại, mô hình NHRCM cho kết quả sát thực tế hơn với độ lệch trung bình chỉ khoảng 0.6 độ C trong các tháng mùa đông và dao động tối đa khoảng 2 độ C tại một số khu vực cục bộ.
  • Tái tạo lượng mưa trung bình năm: RegCM đánh giá thấp nghiêm trọng lượng mưa trên quy mô toàn quốc, mức thiếu hụt trung bình lên tới 800 mm một năm so với số liệu quan trắc trạm. Trong khi đó, NHRCM kiểm soát sai số rất tốt khi mức thiếu hụt chỉ dừng ở khoảng 150 mm một năm, tỷ lệ mô phỏng mưa của NHRCM đạt độ tương thích cao và vượt trội từ 20% đến 70% so với dữ liệu phân tích Aphrodite tại các vùng trọng điểm.
  • Khả năng định vị các tâm mưa lớn: NHRCM bắt trọn vẹn vị trí và cường độ các trung tâm mưa cực đại tại Bắc Quang, Mường Tè ở vùng Tây Bắc, Móng Cái ở Đông Bắc, Bảo Lộc thuộc Lâm Đồng và khu vực Cà Mau ở Nam Bộ. Ngược lại, trường mô phỏng của RegCM làm mờ nhạt hoàn toàn các tâm mưa địa hình này.
  • Mô phỏng chu trình mùa và cực trị ngày: Trong các tháng gió mùa mùa hè và mùa thu, NHRCM bám sát diễn biến biến trình năm của lượng mưa tại các khu vực phía Nam. Về phân bố tần suất ngày, NHRCM tái tạo chính xác tần suất mưa lớn trong ngưỡng từ 70 mm đến 120 mm một ngày tại Bắc Trung Bộ, vượt trội so với độ suy giảm nhanh của RegCM.

Thảo luận kết quả

Sự vượt trội của mô hình NHRCM bắt nguồn từ nền tảng động lực học phi thủy tĩnh, cho phép mô tả chân thực các hoàn lưu đối lưu nhiệt đứng và tương tác phức tạp giữa gió mùa Tây Nam với địa hình đón gió hiểm trở của dãy Trường Sơn và khối núi Hoàng Liên Sơn. Khi xem xét dữ liệu qua các biểu đồ phân tán dạng tụ điểm, các điểm biểu diễn của NHRCM phân bố tập trung sát đường chéo tối ưu 45 độ, trong khi các điểm của RegCM bị trượt sâu xuống phía dưới trục hoành, phản ánh xu hướng đánh giá thấp có hệ thống. Bảng số liệu thống kê sai số trung bình theo từng tháng cũng chỉ ra rằng RegCM gặp nhược điểm nghiêm trọng khi tiến dần về phía Nam Bộ và Nam Trung Bộ do các giả thiết thủy tĩnh không còn phù hợp ở vĩ độ thấp.

Phát hiện này hoàn toàn tương đồng với các kết luận của các chuyên gia khí tượng quốc tế tại Nhật Bản và Cộng hòa Liên bang Đức khi so sánh các phiên bản mô hình khu vực độ phân giải cao. Tuy nhiên, một hạn chế chung được nhận diện là cả hai mô hình đều giảm đáng kể kỹ năng dự báo khi chạm ngưỡng mưa cực đoan trên 150 mm một ngày, đòi hỏi các nghiên cứu tiếp theo phải thu nhỏ bước lưới xuống dưới 5 km.

Đề xuất và khuyến nghị

Dựa trên kết quả thực nghiệm, bốn nhóm giải pháp hành động được đề xuất nhằm nâng cao năng lực ứng dụng mô hình số trị trong nghiên cứu và nghiệp vụ:

  1. Nâng cấp hạ tầng siêu máy tính và chuyển đổi mô hình: Bộ Tài nguyên và Môi trường cùng Trung tâm Dự báo Khí tượng Thủy văn Quốc gia cần ưu tiên đầu tư hệ thống tính toán hiệu năng cao để triển khai chạy nghiệp vụ mô hình phi thủy tĩnh NHRCM với độ phân giải siêu cao từ 2 km đến 5 km, mục tiêu hoàn thành lộ trình tích hợp vào hệ thống dự báo hạn mùa trước năm 2030.
  2. Chuẩn hóa quy trình hiệu chỉnh sai số hệ thống: Các viện nghiên cứu chuyên ngành cần áp dụng bắt buộc phương pháp hiệu chỉnh địa hình với hệ số suy giảm 0.65 độ C trên 100 mét kết hợp hàm phân bố xác suất Gamma cho trường mưa nhằm kéo giảm sai số nhiệt độ xuống dưới 0.5 độ C và sai số lượng mưa dưới 10% trong vòng 12 đến 24 tháng tới.
  3. Tối ưu hóa sơ đồ tham số hóa vật lý cho vùng nhiệt đới: Các trường đại học trọng điểm cần chủ trì các đề tài khoa học thử nghiệm phối hợp sơ đồ Kain-Fritsch cải tiến với mô hình bề mặt đất SiB cho các kịch bản thời tiết cực đoan giai đoạn 2025-2035, tập trung giải quyết hiện tượng mô phỏng thiếu mưa cục bộ từ tháng 4 đến tháng 8 tại Bắc Bộ.
  4. Mở rộng và đồng bộ hóa mạng lưới quan trắc tự động: Tổng cục Khí tượng Thủy văn cần tăng cường mật độ trạm đo mưa tự động từ khoảng cách 50 km hiện nay xuống dưới 15 km trên toàn quốc trước năm 2030, phục vụ việc kiểm định mô hình số trị và hiệu chuẩn các bộ dữ liệu dạng lưới.

Đối tượng nên tham khảo luận văn

Nội dung luận văn mang giá trị ứng dụng thực tiễn và học thuật sâu sắc cho bốn nhóm đối tượng chuyên môn:

  1. Chuyên gia dự báo khí tượng thủy văn: Khai thác cấu hình miền tính 20 km và quy trình thiết lập sơ đồ tham số hóa của mô hình NHRCM để nâng cao độ chính xác trong các bản tin cảnh báo mưa lớn diện rộng và dự báo khí hậu hạn mùa.
  2. Giảng viên, nghiên cứu sinh và học viên sau đại học: Sử dụng tài liệu như một khung tham chiếu chuẩn mực về phương pháp hạ quy mô động lực, kỹ thuật xử lý số liệu tái phân tích JRA-25 và phương pháp đánh giá sai số phân bố thống kê bậc cao ME.
  3. Cơ quan quản lý thiên tai và quy hoạch tài nguyên nước: Vận dụng bản đồ phân bố các tâm mưa lớn tại Bắc Quang, Bảo Lộc và dải ven biển miền Trung để tính toán lưu lượng dòng chảy, điều tiết hồ chứa thủy điện và thiết kế công trình đê điều thủy lợi an toàn.
  4. Cơ quan xây dựng chính sách biến đổi khí hậu: Căn cứ vào bức tranh mô phỏng 22 năm để cập nhật kịch bản biến đổi khí hậu quốc gia, đưa ra các biện pháp thích ứng phù hợp cho từng vùng sinh thái nông nghiệp tại Nam Trung Bộ và Nam Bộ.

Câu hỏi thường gặp

Mô hình NHRCM thể hiện ưu điểm gì vượt trội so với RegCM tại Việt Nam?
Mô hình NHRCM sử dụng hệ phương trình phi thủy tĩnh giúp tái tạo chính xác các chuyển động thẳng đứng và quá trình đối lưu nhiệt đới quy mô nhỏ. Nhờ đó, NHRCM giảm sai số thiếu hụt lượng mưa năm từ 800 mm ở RegCM xuống chỉ còn 150 mm và mô phỏng xuất sắc các tâm mưa lớn tại Mường Tè, Bắc Quang và Bảo Lộc.

Tại sao luận văn sử dụng sai số ME thay vì RMSE để đánh giá mô hình?
Trong dự báo khí hậu dài hạn, sai số hệ thống giữa mô hình và thực tế thường rất lớn. Chỉ số RMSE trộn lẫn cả sai số bậc một và phương sai bậc hai, làm sai lệch nhận định về độ lệch trung bình. Việc sử dụng ME giúp phản ánh trực quan và khách quan bản chất độ lệch hệ thống của từng mô hình.

Hạn chế lớn nhất của hai mô hình khi mô phỏng mưa cực đoan là gì?
Cả hai mô hình đều bộc lộ hạn chế rõ rệt khi nắm bắt các sự kiện mưa cực đoan có cường độ trên 150 mm một ngày. Nguyên nhân chủ yếu do độ phân giải ngang 20 km vẫn chưa đủ mịn để phân giải hoàn toàn các đám mây đối lưu bão tố cực mạnh mà không cần tham số hóa.

Hệ số hiệu chỉnh nhiệt độ theo độ cao 0.65 độ C trên 100 mét có vai trò như thế nào?
Địa hình mô hình số thường bị làm trơn so với địa hình thực tế của các trạm đo trên núi cao. Việc áp dụng gradient suy giảm nhiệt độ 0.65 độ C trên 100 mét giúp đưa giá trị mô phỏng về đúng độ cao hình thái trạm, loại bỏ sai số nhân tạo và phản ánh chính xác nền nhiệt độ bề mặt thực tế.

Bộ số liệu lưới Aphrodite có ưu điểm và nhược điểm gì trong nghiên cứu?
Aphrodite cung cấp bức tranh phân bố lượng mưa liên tục theo không gian với độ phân giải 0.25 độ kinh vĩ suốt 22 năm. Tuy nhiên, do mật độ trạm quan trắc trên bán đảo Đông Dương tương đối thưa khoảng 50 km, thuật ngữ phân tích của Aphrodite đã làm trơn số liệu, khiến lượng mưa cực đại thấp hơn thực tế tại trạm.

Kết luận

  • Mô hình phi thủy tĩnh NHRCM chứng minh tính ưu việt vượt bậc so với mô hình thủy tĩnh RegCM trong việc mô phỏng các trường nhiệt độ và lượng mưa trên khu vực Việt Nam suốt giai đoạn 1986-2007.
  • NHRCM kiểm soát xuất sắc sai số nhiệt độ với độ lệch chỉ từ 0.6 độ C đến 2 độ C và giảm thiểu độ hụt mưa năm xuống mức 150 mm so với mức sai số 800 mm của RegCM.
  • Các trung tâm mưa lớn có tính quy luật địa hình tại Bắc Quang, Mường Tè, Móng Cái, Bảo Lộc và Cà Mau được NHRCM tái lập rõ nét và chuẩn xác.
  • Luận văn khẳng định sự cần thiết phải chuyển dịch từ các mô hình thủy tĩnh truyền thống sang các hệ thống phi thủy tĩnh phân giải cao khi nghiên cứu khí hậu vùng nhiệt đới gió mùa.
  • Nghiên cứu mở ra lộ trình thử nghiệm tiếp theo với bước lưới từ 2 km đến 5 km giai đoạn 2026-2030, đóng góp nguồn tài liệu tham khảo giá trị cho cộng đồng khoa học khí tượng Việt Nam. Các cơ quan nghiên cứu và chuyên gia khí hậu hãy khai thác ngay cơ sở dữ liệu và cấu hình mô hình NHRCM này để tối ưu hóa công tác dự báo thời tiết hạn mùa và xây dựng kịch bản ứng phó biến đổi khí hậu trong tương lai.