Tổng quan nghiên cứu

Biến đổi khí hậu toàn cầu đang diễn ra với tốc độ nhanh chóng, trong đó hiện tượng dâng cao của mực nước biển là một trong những hệ quả nghiêm trọng nhất, đe dọa trực tiếp đến đời sống kinh tế - xã hội của các quốc gia ven biển. Theo các báo cáo quốc tế, mực nước biển trung bình toàn cầu đã tăng khoảng 1,7 ± 0,2 mm/năm trong giai đoạn 1901-2010 với tổng mức dâng tích lũy đạt 0,19 ± 0,2 m, và tốc độ này đã tăng vọt lên mức 3,2 ± 0,4 mm/năm trong giai đoạn quan trắc vệ tinh 1993-2012. Đối với Việt Nam, quốc gia có hơn 3.260 km đường bờ biển trải dài qua 28 tỉnh, thành phố, số liệu viễn thám ghi nhận tốc độ dâng tại Biển Đông lên tới 4,7 mm/năm giai đoạn 1993-2009 và 4,0 mm/năm giai đoạn 1993-2013, cao hơn đáng kể so với mức trung bình toàn cầu.

Vấn đề nghiên cứu trọng tâm của luận văn là khắc phục các hạn chế của những phương pháp ngoại suy thống kê đơn thuần trước đây, vốn chưa phản ánh đầy đủ bản chất vật lý và động lực học hải dương học khu vực. Mục tiêu cụ thể của công trình là đánh giá toàn diện xu thế biến đổi mực nước biển trong quá khứ tại vùng biển Việt Nam từ chuỗi số liệu thực đo và vệ tinh, đồng thời ứng dụng kết quả mô phỏng từ 21 mô hình hoàn lưu chung khí quyển - đại dương (AOGCMs) thuộc Dự án Đối chứng Mô hình Kết hợp giai đoạn 5 (CMIP5) để dự tính chi tiết mực nước biển dâng cho khu vực Biển Đông đến năm 2100 theo hai kịch bản nồng độ khí nhà kính đại diện RCP4.5 và RCP8.5. Phạm vi không gian bao trùm toàn bộ Biển Đông và dải ven biển Việt Nam với mốc thời gian nền chuẩn 1986-2005. Kết quả nghiên cứu mang ý nghĩa thực tiễn sâu sắc, cung cấp cơ sở dữ liệu định lượng tin cậy giúp các nhà quản lý xây dựng chiến lược ứng phó, quy hoạch hạ tầng ven bờ và bảo vệ sinh kế bền vững cho hàng chục triệu người dân.

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

Khung lý thuyết áp dụng

Luận văn được xây dựng dựa trên khung lý thuyết đánh giá biến đổi khí hậu của Ban Liên chính phủ về Biến đổi Khí hậu (IPCC) trong Báo cáo Đánh giá lần thứ năm (AR5). Nghiên cứu tích hợp hệ thống kịch bản đường nồng độ đại diện (RCP), tập trung vào hai kịch bản chính gồm RCP4.5 với mức bức xạ cưỡng bức đạt 4,5 W/m² (tương đương nồng độ khí nhà kính ổn định khoảng 650 ppm CO₂ tương đương vào năm 2100) và RCP8.5 với mức bức xạ cưỡng bức 8,5 W/m² (kịch bản phát thải cao với nồng độ vượt 1.370 ppm CO₂ tương đương vào năm 2100).

Mô hình nghiên cứu phân rã toàn diện 5 thành phần vật lý chính đóng góp vào mực nước dâng: hiệu ứng giãn nở nhiệt của nước biển (thermosteric), sự tan chảy của các dòng sông băng trên núi cao (glaciers), cân bằng khối lượng bề mặt và động lực học của hai dải băng lớn Greenland và Nam Cực, sự thay đổi trữ lượng nước trên lục địa do khai thác và tích trữ, cùng hiệu ứng điều chỉnh đẳng tĩnh băng (Glacial Isostatic Adjustment - GIA) thông qua mô hình động học vỏ Trái Đất ICE-5G VM2. Bên cạnh đó, các khái niệm then chốt như độ cao mặt biển động lực học, mực chuẩn Geoid và dao động khí hậu quy mô thập kỷ (như hiện tượng El Niño - Dao động Nam ENSO với chu kỳ 2-7 năm và Dao động thập kỷ Thái Bình Dương PDO) được áp dụng đồng bộ để làm rõ các biến thiên mực nước cục bộ.

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

Nghiên cứu sử dụng nguồn dữ liệu đa tầng bao gồm: chuỗi số liệu đo cao vệ tinh AVISO giai đoạn 1993-2013 với độ phân giải không gian 0,25 độ kinh vĩ; chuỗi số liệu mực nước thực đo tại 17 trạm hải văn ven bờ Việt Nam với thời gian quan trắc liên tục từ 30 đến 50 năm; và tập dữ liệu mô phỏng từ 21 mô hình AOGCMs toàn cầu thuộc CMIP5. Cỡ mẫu nghiên cứu bao gồm 21 mô hình khí hậu kết hợp đại dương - khí quyển độc lập, cùng 17 trạm quan trắc thực địa ven bờ đại diện đầy đủ cho 3 miền khí hậu hải văn Bắc, Trung, Nam.

Phương pháp chọn mẫu được tiến hành theo tiêu chí lựa chọn toàn bộ các mô hình CMIP5 có đầy đủ biến số về nhiệt độ tầng nước, độ mặn và động lực học mặt biển, đồng thời chọn các trạm hải văn có chuỗi số liệu dài, độ tin cậy thống kê cao và liên tục. Luận văn lựa chọn phương pháp phân tích hiệu chỉnh sai số (bias correction), kỹ thuật chi tiết hóa không gian kết hợp nội suy khoảng cách nghịch đảo (IDW) và phân tích thống kê đa biến. Lý do lựa chọn tổ hợp phương pháp này là nhằm loại bỏ độ lệch có hệ thống của các mô hình toàn cầu ở vùng biển nông ven bờ, đồng thời bảo toàn tính chính xác của các quy luật vật lý khi chuyển đổi từ lưới mô hình toàn cầu sang quy mô chi tiết 0,5 x 0,5 độ tại Biển Đông. Timeline nghiên cứu được xác lập xuyên suốt từ thời kỳ nền 1986-2005 đến giai đoạn giữa thế kỷ (2046-2065) và cuối thế kỷ 21 (2081-2100).

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

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

Thứ nhất, phân tích chuỗi dữ liệu lịch sử cho thấy tốc độ dâng mực nước biển tại Biển Đông có sự gia tăng rõ rệt: số liệu quan trắc từ 17 trạm hải văn ven bờ ghi nhận mức tăng trung bình đạt 2,8 đến 3,0 mm/năm, trong khi dữ liệu vệ tinh giai đoạn 1993-2013 chỉ ra tốc độ dâng lên tới 4,0 ± 0,3 mm/năm, phản ánh sự gia tốc mạnh mẽ trong hai thập kỷ gần đây.

Thứ hai, kết quả dự tính từ tổ hợp 21 mô hình AOGCMs cho thấy vào giữa thế kỷ 21 (2046-2065), mực nước biển trung bình khu vực Biển Đông sẽ dâng từ 22 đến 28 cm theo kịch bản RCP4.5 và từ 25 đến 32 cm theo kịch bản RCP8.5 so với thời kỳ nền 1986-2005. Đến cuối thế kỷ 21 (2081-2100), mức dâng tổng cộng đạt 52 đến 62 cm (giá trị trung vị khoảng 57 cm) đối với kịch bản RCP4.5 và từ 73 đến 88 cm (giá trị trung vị khoảng 79 cm) đối với kịch bản RCP8.5.

Thứ ba, nghiên cứu xác định cơ cấu đóng góp của các thành phần vật lý: giãn nở nhiệt đại dương là yếu tố chi phối lớn nhất, chiếm từ 38% đến 45% tổng mức dâng; băng tan trên đỉnh núi đóng góp khoảng 25% đến 28%; băng tan từ dải băng Greenland và Nam Cực đóng góp 18% đến 22%; phần còn lại khoảng 5% đến 10% bắt nguồn từ thay đổi trữ lượng nước lục địa và hiệu ứng điều chỉnh đẳng tĩnh địa chất GIA.

Thứ tư, phân bố mực nước biển dâng thể hiện tính dị thường không gian rõ rệt: khu vực biển ven bờ Trung Bộ và Nam Bộ có xu thế dâng cao hơn vùng vịnh Bắc Bộ từ 3 đến 7 cm, do đặc điểm dòng chảy hoàn lưu mùa hè và tương tác địa hình thềm lục địa sâu hơn 1.000 m.

Thảo luận kết quả

Nguyên nhân chính dẫn đến sự chênh lệch mực nước dâng giữa Biển Đông và trung bình toàn cầu là do sự kết hợp của hiệu ứng nhiệt động lực học đại dương khu vực và các dao động khí hậu quy mô lớn. Hoàn lưu gió mùa Đông Bắc và Tây Nam xen kẽ hàng năm, cùng ảnh hưởng từ các pha nóng El Niño và pha lạnh La Niña (như đợt La Niña mạnh năm 2011), đã làm tái phân bố khối nước tầng mặt và gây biến thiên mực nước liên năm từ 5 đến 15 mm. Khi so sánh với Kịch bản Biến đổi khí hậu năm 2012 của Bộ Tài nguyên và Môi trường (vốn dựa trên kịch bản phát thải SRES cũ và phương pháp thống kê đơn giản), kết quả dự tính mới từ 21 mô hình CMIP5 có độ tương quan cao (hệ số r > 0,85) nhưng cung cấp khoảng tin cậy hẹp hơn và phản ánh chân thực hơn sự phân hóa không gian.

Các kết quả tính toán trong luận văn được trực quan hóa rất hiệu quả thông qua biểu đồ biến trình thời gian (time series plot) thể hiện đường trung bình đa mô hình kèm dải bóng mờ biểu thị khoảng tin cậy 90% từ năm 2006 đến 2100. Đồng thời, các bản đồ trường đẳng trị không gian (spatial contour maps) lưới 0,5 x 0,5 độ kết hợp bảng số liệu phân rã chi tiết 5 thành phần vật lý theo từng mốc thời gian 10 năm đã giúp làm rõ bức tranh ngập lụt tiềm tàng, tạo tiền đề vững chắc cho việc xây dựng bản đồ nguy cơ ngập úng ven bờ.

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

Dựa trên kết quả dự tính khoa học, luận văn đề xuất 4 nhóm giải pháp chiến lược nhằm giảm thiểu rủi ro và tăng cường khả năng thích ứng cho dải ven biển Việt Nam:

  1. Nâng cấp và kiên cố hóa hệ thống đê biển, kè chống sạt lở: Tiến hành rà soát, nâng cao trình thiết kế cho hơn 3.260 km đê biển trên toàn quốc nhằm đảm bảo khả năng chống chịu mức nước biển dâng 30 cm vào năm 2030 và 80 cm vào năm 2050 kết hợp triều cường cực đoan; giao Bộ Nông nghiệp và Phát triển nông thôn chủ trì phối hợp với Ủy ban nhân dân 28 tỉnh, thành phố ven biển triển khai trong giai đoạn 2025-2035.

  2. Điều chỉnh quy hoạch cốt xây dựng và hạ tầng đô thị ven biển: Bắt buộc áp dụng cốt nền xây dựng mới tăng thêm từ 0,5 đến 0,8 m đối với các khu đô thị mới, khu công nghiệp và hạ tầng giao thông huyết mạch tại các địa bàn trũng thấp như Hải Phòng, Đà Nẵng, Cần Thơ và Thành phố Hồ Chí Minh; giao Bộ Xây dựng ban hành quy chuẩn kỹ thuật trước năm 2026 và các địa phương hoàn thành rà soát quy hoạch tổng thể trước năm 2028.

  3. Phục hồi và mở rộng hành lang rừng ngập mặn phòng hộ: Triển khai trồng mới và khoanh nuôi tái sinh tối thiểu 50.000 ha rừng ngập mặn tại các cửa sông và dải ven biển Đồng bằng sông Cửu Long nhằm giảm từ 30% đến 50% năng lượng sóng biển tàn phá đê kè khi xảy ra bão mạnh; giao Bộ Tài nguyên và Môi trường phối hợp cùng chính quyền địa phương thực hiện trong lộ trình 2025-2030.

  4. Mở rộng mạng lưới quan trắc hải văn tự động và nâng cấp hệ thống cảnh báo sớm: Lắp đặt thêm ít nhất 20 trạm đo triều tự động truyền dữ liệu vệ tinh thời gian thực và tích hợp mô hình dự báo ngập lụt độ phân giải cao vào hệ thống cảnh báo quốc gia, phấn đấu nâng độ chính xác dự báo ngập lụt ven bờ lên trên 90% vào năm 2028; giao Tổng cục Khí tượng Thủy văn trực tiếp vận hành.

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

Luận văn thạc sĩ là tài liệu tham khảo học thuật và ứng dụng thực tiễn chuyên sâu, đặc biệt hữu ích cho 4 nhóm đối tượng sau:

  1. Giảng viên, nghiên cứu sinh và học viên cao học chuyên ngành Khí tượng, Thủy văn, Hải dương học và Biến đổi khí hậu: Sử dụng luận văn như tài liệu tham khảo chuẩn mực về phương pháp xử lý dữ liệu viễn thám AVISO, kỹ thuật khai thác tổ hợp 21 mô hình khí hậu CMIP5 và mô hình hóa phân rã các hợp phần vật lý hải dương.

  2. Cơ quan quản lý nhà nước và nhà hoạch định chính sách môi trường: Cán bộ chuyên trách tại Bộ Tài nguyên và Môi trường, Bộ Nông nghiệp và Phát triển nông thôn cùng Sở Tài nguyên và Môi trường của 28 tỉnh ven biển có thể ứng dụng trực tiếp các số liệu dự tính 52-88 cm mực nước dâng để cập nhật Kịch bản Biến đổi khí hậu địa phương giai đoạn 2025-2050.

  3. Kỹ sư tư vấn thiết kế công trình thủy lợi, cảng biển và hạ tầng giao thông ven bờ: Khai thác số liệu về biên độ dâng cực đại và phân bố không gian để tính toán tải trọng sóng, cao trình đỉnh đê, độ sâu luồng hàng hải và tuổi thọ công trình với vòng đời dự án từ 50 đến 100 năm.

  4. Các tổ chức phi chính phủ quốc tế và đơn vị tài trợ phát triển bền vững: Sử dụng kết quả đánh giá nguy cơ ngập lụt để xây dựng các đề án hỗ trợ tài chính, sinh kế thích ứng biến đổi khí hậu cho hơn 20 triệu người dân tại các vùng nhạy cảm như Đồng bằng sông Cửu Long và Duyên hải miền Trung.

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

  1. Tại sao tốc độ dâng mực nước biển tại Biển Đông lại cao hơn trung bình toàn cầu? Số liệu vệ tinh giai đoạn 1993-2013 cho thấy tốc độ dâng tại Biển Đông đạt 4,0 mm/năm, cao hơn mức 3,2 mm/năm của toàn cầu. Nguyên nhân là do Biển Đông là biển nửa kín nằm trong khu vực nhiệt đới gió mùa, chịu ảnh hưởng mạnh mẽ từ hiện tượng giãn nở nhiệt tầng nước mặt và sự thay đổi hoàn lưu gió mùa kết hợp với các chu kỳ dao động ENSO kéo dài từ 2 đến 7 năm.

  2. Việc sử dụng tổ hợp 21 mô hình AOGCMs mang lại ưu thế gì so với một mô hình đơn lẻ? Mỗi mô hình khí hậu có độ phân giải và cách tham số hóa vật lý riêng biệt nên dễ dẫn đến sai số cục bộ. Việc tích hợp tổ hợp 21 mô hình thuộc CMIP5 giúp triệt tiêu các sai số ngẫu nhiên, mở rộng không gian khảo sát và cung cấp dải dự tính với khoảng tin cậy 90%, nâng cao tính chính xác và độ tin cậy khoa học của kết quả.

  3. Kịch bản RCP4.5 và RCP8.5 dự tính mức dâng khác nhau như thế nào vào cuối thế kỷ 21? Vào giai đoạn 2081-2100, theo kịch bản phát thải trung bình thấp RCP4.5 (bức xạ cưỡng bức 4,5 W/m²), mực nước biển khu vực Biển Đông dự tính dâng từ 52 đến 62 cm. Trong khi đó, theo kịch bản phát thải cao RCP8.5 (bức xạ cưỡng bức 8,5 W/m²), mức dâng tăng vọt lên 73 đến 88 cm so với thời kỳ nền 1986-2005.

  4. Hợp phần vật lý nào đóng góp nhiều nhất vào mực nước biển dâng tại Việt Nam? Quá trình giãn nở nhiệt của nước biển do đại dương hấp thụ nhiệt lượng từ khí quyển là thành phần đóng góp lớn nhất, chiếm từ 38% đến 45% tổng lượng dâng. Tiếp theo là lượng nước tan từ các sông băng trên núi cao chiếm khoảng 25% đến 28%, và sự suy giảm thể tích băng ở Greenland và Nam Cực chiếm khoảng 18% đến 22%.

  5. Số liệu dự tính trong luận văn có tác động như thế nào đến kế hoạch phát triển Đồng bằng sông Cửu Long? Với địa hình có hơn 60% diện tích thấp dưới 1,0 m so với mực nước biển, mức dâng dự tính 73 đến 88 cm theo kịch bản RCP8.5 sẽ khiến diện tích ngập úng tăng thêm từ 20% đến 35%, đe dọa trực tiếp vựa lúa quốc gia. Số liệu này là căn cứ tối quan trọng để chuyển đổi mô hình nông nghiệp thích ứng và xây dựng hệ thống đê bao khép kín.

Kết luận

  • Luận văn đã xác lập cơ sở khoa học định lượng vững chắc khi dự tính mực nước biển dâng tại khu vực Biển Đông vào cuối thế kỷ 21 sẽ đạt từ 52 đến 62 cm theo kịch bản RCP4.5 và từ 73 đến 88 cm theo kịch bản RCP8.5 so với thời kỳ nền 1986-2005.
  • Đột phá phương pháp luận của công trình nằm ở việc kết hợp đồng bộ 21 mô hình AOGCMs thuộc CMIP5 với chuỗi số liệu thực đo tại 17 trạm hải văn và vệ tinh đo cao trong hơn 20 năm, giải quyết triệt để bài toán phân rã 5 hợp phần vật lý hải dương học.
  • Đóng góp chính của nghiên cứu là cung cấp bức tranh chi tiết hóa không gian về độ dâng mực nước biển cho toàn dải ven bờ Việt Nam, chỉ rõ vùng biển Trung Bộ và Nam Bộ chịu mức dâng cao hơn Bắc Bộ từ 3 đến 7 cm.
  • Lộ trình nghiên cứu tiếp theo trong giai đoạn 2025-2030 cần tích hợp dữ liệu dâng mực nước biển này vào các mô hình thủy lực 2D/3D độ phân giải cao tỉ lệ 1:10.000 để mô phỏng chính xác phạm vi ngập úng và xâm nhập mặn tại từng xã ven biển.
  • Kêu gọi các bộ, ngành và chính quyền 28 tỉnh, thành ven biển cần khẩn trương lồng ghép ngay các chỉ số dự tính của luận văn vào quy hoạch tổng thể kinh tế - xã hội, bảo vệ an toàn sinh kế cho hàng chục triệu người dân trước thách thức biến đổi khí hậu.