Tổng quan về luận án
Biến động địa hình bờ biển, đặc biệt là hiện tượng xói lở bờ biển quy mô lớn, đang trở thành một thách thức địa mạo - môi trường toàn cầu trong bối cảnh biến đổi khí hậu và nước biển dâng. Tại đồng bằng châu thổ sông Mê Công, dải bờ biển Sóc Trăng - Cà Mau (STCM) vốn là khu vực bồi tích lấn biển điển hình với tốc độ bồi tụ lịch sử trung bình đạt 16 m/năm, riêng khu vực Mũi Cà Mau từng đạt 26 m/năm theo tính toán của Xue và cộng sự (2010). Tuy nhiên, những năm gần đây, xu thế bồi - xói đã đảo chiều phức tạp; hiện tượng xói lở nghiêm trọng đã phá hủy hàng loạt công trình đê kè, rừng ngập mặn và cơ sở hạ tầng ven biển suốt chiều dài bán đảo Cà Mau, đe dọa trực tiếp đến bãi Khai Long, Gành Hào, Rạch Gốc và Đất Mũi.
Nghiên cứu của NCS. Lưu Thành Trung (2017) với đề tài "Nghiên cứu đánh giá biến động đường bờ phục vụ quản lý bờ biển Sóc Trăng - Cà Mau" thuộc chuyên ngành Quản lý tài nguyên và môi trường (Mã số: 62 85 01 01) tại Trường Đại học Khoa học Tự nhiên, Đại học Quốc gia Hà Nội, dưới sự hướng dẫn khoa học của PGS.TS. Vũ Văn Phái và TS. Nguyễn Đắc Đồng, đã giải quyết trực tiếp khoảng trống nghiên cứu (research gap) cốt lõi: Sự thiếu hụt một khung tích hợp liên ngành giữa động lực học địa mạo (morphodynamics), địa tin học viễn thám đa thời gian, mô hình hóa thủy thạch động lực 2 chiều và đánh giá tính dễ bị tổn thương để thiết lập các công cụ quản lý bờ biển theo luật định.
Luận án xác lập hai câu hỏi nghiên cứu và hai giả thuyết khoa học tương ứng:
- Câu hỏi nghiên cứu 1 (RQ1): Hiện trạng, xu thế dài hạn (1965–2013) và quy luật biến động không gian - thời gian của địa hình đường bờ và đáy biển vùng STCM diễn ra như thế nào dưới tác động tổng hợp của các nhân tố tự nhiên và nhân sinh?
- Câu hỏi nghiên cứu 2 (RQ2): Làm thế nào để tích hợp đánh giá tính dễ bị tổn thương hình thái (CVI) và xác lập hành lang cảnh báo nguy cơ xói lở - bồi tụ (EPAW) nhằm phân vùng định hướng quản lý bờ biển hiệu quả theo Luật Tài nguyên, Môi trường biển và hải đảo năm 2015?
- Giả thuyết khoa học 1 (H1): Dưới tác động tổng hợp của chế độ động lực biển, sự suy giảm trầm tích lưu vực sông Mê Công và các hoạt động nhân sinh trong bối cảnh biến đổi khí hậu, đường bờ biển STCM giai đoạn 1965–2013 đã chuyển hóa cấu trúc hình thái từ trạng thái ưu thế bồi tụ sang xu thế gia tăng xói lở cục bộ và phân hóa sâu sắc theo các tiểu vùng địa mạo.
- Giả thuyết khoa học 2 (H2): Việc tích hợp phân tích xu thế biến động đường bờ lịch sử qua chỉ thị thực vật - viễn thám với chỉ số tổn thương CVI và mô hình tính toán vùng xói lở tiềm năng cho phép thiết lập hành lang cảnh báo nguy cơ tai biến và phân vùng chức năng quản lý vùng bờ có cơ sở khoa học định lượng vững chắc.
Phạm vi nghiên cứu bao quát dải bờ biển dài hàng trăm kilômét kéo dài từ cửa Định An (tỉnh Sóc Trăng) đến Tiểu Dừa (tỉnh Cà Mau), thuộc địa giới hành chính của 12 đơn vị cấp huyện (Cù Lao Dung, Trần Đề, Vĩnh Châu, TP. Bạc Liêu, Hòa Bình, Đông Hải, Đầm Dơi, Năm Căn, Ngọc Hiển, Phú Tân, Trần Văn Thời, U Minh) thuộc 3 tỉnh Tây Nam Bộ. Phạm vi không gian mở rộng từ đất liền ra phía biển đến đẳng sâu 20 m nước, khung thời gian khảo sát động thái kéo dài 48 năm (1965–2013) cùng các chuỗi số liệu mô phỏng thủy văn - động lực chi tiết năm 2009 và 2012.
Literature Review và Positioning
Nghiên cứu tiến hóa và biến động địa hình đới bờ biển là chủ đề trọng tâm của địa mạo học động lực (process geomorphology). Từ nền tảng lý thuyết tiến hóa địa hình trong hệ mở của Gilbert (1896), Strahler (1952), Hack (1960) và lý thuyết hệ thống địa mạo của Chorley (1962), bờ biển được nhìn nhận là một hệ thống cân bằng động nhạy cảm. Zencovich (1962) và Bird (2008) đã hệ thống hóa các quá trình hình thành, biến đổi bờ biển dưới tác động tương hỗ giữa sóng, triều, dòng chảy và trầm tích.
Trong y văn quốc tế tồn tại hai luồng quan điểm đối thoại sâu sắc về động lực biến đổi bờ biển bùn sét nhiệt đới có rừng ngập mặn:
- Luồng quan điểm địa chất - sinh thái học truyền thống (Woodroffe, 1992, 2002; Pedro Walfir et al., 2006): Rừng ngập mặn đóng vai trò như "bẫy trầm tích" và "kẻ kiến tạo bờ biển", làm giảm năng lượng sóng, giữ bùn cát mịn và thúc đẩy quá trình bồi tụ lấn biển ra phía trước.
- Luồng quan điểm động lực hình thái hiện đại (Slaymaker, 2009; Stephenson & Brander, 2003; Wilson et al., 2012): Quần xã rừng ngập mặn chỉ đóng vai trò là "chỉ thị thích ứng bị động" trên nền địa hình; chính sự biến động của nguồn cung trầm tích thượng nguồn kết hợp với biến đổi dòng dọc bờ và thay đổi sử dụng đất nhân sinh mới là tác nhân quyết định "số phận của vật liệu trầm tích" và vị trí thực tế của đường bờ.
Tại Việt Nam, các nghiên cứu của Nguyễn Địch Dỹ (1995, 2010), Vũ Trường Sơn (2004), Nguyễn Huy Dũng và Ngô Quang Toàn (2004) đã làm rõ tiến hóa trầm tích Đệ Tứ và nguồn gốc đồng bằng châu thổ sông Cửu Long sau pha biển tiến Flandrian cực đại cách đây 6.000 năm. Các công trình của Nguyễn Văn Cư và Phạm Huy Tiến (2003), Nguyễn Chu Hồi (1997), Đào Đình Bắc (2002), Vũ Văn Phái (2012, 2013) đã bước đầu tiếp cận xói lở bờ biển dưới góc độ tai biến địa mạo và điều tra cơ bản. Tuy nhiên, phần lớn các công trình trước đây chỉ dừng lại ở việc mô tả hiện trạng xói lở cục bộ tại miền Trung hoặc phân tích định tính tại các cửa sông lớn, chưa có công trình nào tích hợp toàn diện chuỗi viễn thám đa thời gian gần nửa thế kỷ với mô hình hóa số trị 2D đáy biển và chỉ số CVI cho toàn bộ hệ thống bờ biển bùn sét Tây Nam Bộ.
So sánh với các nghiên cứu quốc tế, luận án đặt vị trí nghiên cứu tương đồng với các công trình kinh điển:
- Nghiên cứu của Wilson và cộng sự (2012) tại vùng bờ biển rừng ngập mặn cửa sông Amazon (Brazil) giai đoạn 1996–2008, nơi các tác giả tích hợp Landsat và giải đoán chỉ thị thực vật trong môi trường bùn sét nhiệt đới bở rời.
- Khung phân tích quản lý nguy cơ tai biến đới bờ Địa Trung Hải của Sanò và cộng sự (2010) và nghiên cứu xác lập chỉ giới lùi bờ biển (setback lines) của Cambers (1998) tại vùng biển Caribe.
- Hệ thống quản lý đường bờ pháp lý và chỉ giới hạn chế phát triển quốc gia của Hà Lan thiết lập từ năm 1990 dựa trên chuỗi quan trắc liên tục từ 1963.
Luận án của Lưu Thành Trung định vị một bước tiến đột phá khi cụ thể hóa lý thuyết địa mạo học ứng dụng vào bối cảnh vùng ven biển đồng bằng bùn sét cận xích đạo chịu tác động của chế độ bán nhật triều không đều và nhật triều, làm cầu nối trực tiếp giữa khoa học trái đất và thực thi chính sách quản lý tài nguyên.
Đó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 học thuyết cân bằng động (Dynamic Equilibrium Theory) của Hack (1960) và Strahler (1952) bằng cách bổ sung cơ chế phản hồi phi tuyến của hệ thống địa mạo - sinh vật - trầm tích trong điều kiện chịu áp lực kép: suy thoái nguồn cấp phù sa sông xuyên biên giới kết hợp với biến đổi khí hậu toàn cầu.
Tác giả luận chứng rằng: Trầm tích ven biển không chỉ di chuyển thụ động dưới tác động của trường sóng và dòng chảy, mà quần xã rừng ngập mặn (Avicennia - Cây Mắm, Rhizophora - Cây Đước) hoạt động như một hệ chỉ thị địa mạo nhạy cảm. Sự phân tầng sinh thái, mật độ rễ khí sinh và chiều cao giảm dần từ đất liền ra biển phản ánh pha bồi tụ (aggradation), trong khi sự xuất hiện của vách xói lở đứng, bờ răng cưa và cây ngập mặn cổ thụ bị bật gốc là bằng chứng vật lý của pha bóc mòn xâm thực (degradation).
[Nguồn cung trầm tích sông Mê Công & Dòng dọc bờ]
│
▼
┌─────────────────────────────┐
│ CÂN BẰNG ĐỘNG ĐỚI BỜ │
└──────────────┬──────────────┘
│
┌────────────┴────────────┐
▼ ▼
[PHA BỒI TỤ] [PHA XÓI LỞ]
• Mắm tiên phong • Bờ dạng răng cưa
• Rừng hạ thấp ra biển • Cây cổ thụ bật gốc
• Đáy biển bồi lắng • Sụt lở vách triều
│ │
└────────────┬────────────┘
│
▼
[CHỈ SỐ DỄ BỊ TỔN THƯƠNG ĐƯỜNG BỜ (CVI)]
│
▼
[HÀNH LANG CẢNH BÁO NGUY CƠ (EPAW) & PHÂN VÙNG QUẢN LÝ]
Luận án tạo ra sự dịch chuyển nhận thức luận (paradigm shift): Thay vì xem xói lở bờ biển là một "sự cố kỹ thuật" cần ngăn chặn bằng các công trình đê kè cứng tốn kém, nghiên cứu khẳng định xói lở - bồi tụ là hai mặt đối lập thống nhất của một quá trình tự nhiên trong hệ bờ mở; giải pháp quản lý bền vững phải dựa trên nguyên tắc "biển tự xây bờ" và bảo toàn các dịch vụ địa mạo (geomorphic services).
Khung phân tích độc đáo
Khung phân tích của luận án là sự tích hợp đa tầng giữa 3 lý thuyết nền tảng:
- Thuyết Hệ thống Tổng quát (General Systems Theory - Bertalanffy, 1956): Xem đới bờ là hệ mở tương tác vật chất - năng lượng giữa biển, lục địa và sinh quyển.
- Lý thuyết Địa mạo Môi trường & Động lực hình thái (Environmental Geomorphology - Panizza, 1996; Boak & Turner, 2005): Phân tích tương tác giữa động lực thủy thạch và hình thái đường bờ.
- Lý thuyết Đánh giá Tổn thương & Quản lý tổng hợp vùng bờ (ICZM & Vulnerability Theory - Thieler & Hammar, 2000; Pendleton et al., 2010): Lượng hóa rủi ro tự nhiên thành các chỉ số không gian phục vụ quy hoạch.
Điều kiện biên (boundary conditions) của khung phân tích được xác lập nghiêm ngặt: Áp dụng tối ưu cho các đới bờ biển tích tụ bùn sét mịn nhiệt đới, có thảm thực vật ngập mặn phát triển, độ dốc đáy biển thoải (< 1°), biên độ triều lớn (từ 1,0 m đến trên 3,5 m) và chịu sự kiểm soát trầm tích mạnh mẽ từ mạng lưới sông ngòi châu thổ lớn.
Phương pháp nghiên cứu tiên tiến
Thiết kế nghiên cứu
Luận án vận hành trên lập trường nhận thức luận Hiện thực Phê phán (Critical Realism) kết hợp Thực chứng (Positivism), sử dụng phương pháp luận Phân tích hệ thống làm sợi chỉ đỏ. Cấu trúc nghiên cứu đa cấp độ (multi-level design) kết hợp đồng thời phương pháp phân tích vĩ mô (chuỗi dữ liệu viễn thám 48 năm), trung mô (mô hình hóa số trị 2 chiều mùa gió) và vi mô (khảo sát thực địa chi tiết theo mặt cắt và phỏng vấn xã hội học môi trường).
┌────────────────────────────────────────┐
│ BƯỚC 1: THU THẬP & CHUẨN HÓA DỮ LIỆU │
└───────────────────┬────────────────────┘
│
┌─────────────────────────────┼─────────────────────────────┐
▼ ▼ ▼
[Viễn thám & Bản đồ] [Thực địa & GPS] [Số liệu Hải văn/Độ sâu]
• Bản đồ UTM 1965 (1:50k) • 03 đợt (2012, 2013, 2015) • Độ sâu 0-20m Viện ĐC&ĐVB
• Landsat 5 TM (1990) • Lộ trình ven bờ + vỏ lãi • Khí tượng thủy văn trạm QG
• Landsat 8 OLI (2013) • 45 chỉ thị địa mạo • Nồng độ trầm tích lơ lửng
• SPOT 2 (1990), SPOT 5 (2013) • Phỏng vấn dân cư, kiểm lâm
│ │ │
└─────────────────────────────┼─────────────────────────────┘
│
┌───────────────────┴────────────────────┐
│ BƯỚC 2: PHÂN TÍCH & MÔ PHỎNG │
└───────────────────┬────────────────────┘
│
┌─────────────────────────────┴─────────────────────────────┐
▼ ▼
[Động lực Đường bờ DSAS] [Mô phỏng Thủy thạch MIKE-21]
• Tính tốc độ EPR & NSM • Mô đun HD (Nước nông 2D)
• Đối chiếu 48 năm biến động • Mô đun SW (Trường sóng cực đại/TB)
• Bóc tách bồi - xói chi tiết • Mô đun MT (Bồi xói đáy & bùn cát)
│ │
└─────────────────────────────┬─────────────────────────────┘
│
┌───────────────────┴────────────────────┐
│ BƯỚC 3: TỔNG HỢP & ỨNG DỤNG QUẢN LÝ │
└───────────────────┬────────────────────┘
│
┌─────────────────────────────┴─────────────────────────────┐
▼ ▼
[Chỉ số Tổn thương CVI] [Hành lang Cảnh báo EPAW]
• Tích hợp 6 tham số hình thái • Tính toán 3 xã Nam Mũi Cà Mau
• Phân tích đường cong lũy tích • Thiết lập phân vùng quản lý bờ
Quy trình nghiên cứu rigorous
Quy trình thu thập và xử lý dữ liệu trải qua các bước kiểm định nghiêm ngặt:
- Dữ liệu viễn thám và bản đồ: Sử dụng bản đồ địa hình UTM tỷ lệ 1:50.000 năm 1965; các cảnh ảnh vệ tinh Landsat 5 TM (tháng 8/1990) và Landsat 8 OLI (tháng 8/2013) với mã hàng/cột chuẩn xác (Path 125–126, Row 53–54) nhằm đảm bảo đồng nhất chu kỳ mùa và hạn chế tối đa độ che phủ của mây; ảnh SPOT 2 (1990) và SPOT 5 (2013) độ phân giải cao kết hợp kênh toàn sắc Panchromatic và đa phổ MultiSpectral (P+XS).
- Khảo sát thực địa: Tác giả tiến hành 3 đợt khảo sát thực địa thực tế vào các năm 2012, 2013 và 2015. Các lộ trình được thực hiện bằng ô tô, xe máy, đi bộ dọc đê biển và đặc biệt sử dụng thuyền vỏ lãi tiếp cận từ phía biển đối với các vùng sình lầy hiểm trở tại Vườn Quốc gia Mũi Cà Mau, bãi Gành Hào, Rạch Gốc, Rạch Tàu.
- Tam giác hóa phương pháp (Methodological Triangulation): Kết hợp đồng thời 4 trụ cột: (1) Đo vẽ bản đồ viễn thám số DSAS; (2) Đo đạc định vị GPS thực địa; (3) Mô hình hóa động lực học MIKE-21; (4) Khảo sát chỉ thị sinh học thực vật và cấu trúc trầm tích lỗ khoan bãi triều.
Data và phân tích
Phân tích biến động đường bờ được thực hiện thông qua công cụ chuyên dụng Digital Shoreline Analysis System (DSAS v4.2) do Cục Khảo sát Địa chất Hoa Kỳ (USGS) phát triển chạy trên môi trường ArcGIS 10. Dữ liệu vector đường bờ được hiệu chỉnh hình học về hệ quy chiếu chuẩn Quốc gia VN-2000 lưới chiếu UTM. Hệ thống các mặt cắt vuông góc (transects) với đường cơ sở (baseline) được thiết lập tự động để tính toán hai chỉ số cốt lõi: Tốc độ dịch chuyển điểm cuối (End Point Rate - EPR, đơn vị m/năm) và Khoảng cách biến động thực tế (Net Shoreline Movement - NSM, đơn vị mét).
Mô phỏng thủy động lực học và vận chuyển bùn cát đáy biển sử dụng bộ phần mềm MIKE-21 của Viện Thủy lực Đan Mạch (DHI) gồm 3 mô đun liên hoàn:
- Mô đun Thủy động lực học nước nông 2 chiều (MIKE-21 HD): Giải hệ phương trình liên tục và Navier-Stokes 2D tích phân theo độ sâu với tham số Coriolis ($f = 2\Omega\sin\varphi$) và tenxơ ứng suất bức xạ sóng ($S_{xx}, S_{xy}, S_{yy}$).
- Mô đun Sóng phổ (MIKE-21 SW): Tính toán mômen bậc 0 ($m_0$) và bậc 1 ($m_1$) của phổ sóng tác động, giải bài toán truyền sóng cực đại ngày 24/5/2012 và trường sóng trung bình các mùa ngày 04/7/2012, 15/12/2012.
- Mô đun Vận chuyển bùn cát hạt mịn (MIKE-21 MT): Tính toán biến động nồng độ trầm tích lơ lửng và cân bằng lắng đọng ($D$) - xói mòn ($E$) tại lớp đáy theo bước thời gian $\Delta t$ trên lưới tính toán độ sâu đáy biển năm 2009 và 2012.
Đánh giá tính dễ bị tổn thương bờ biển áp dụng công thức bán định lượng của Robert & Erika (2000) và Thieler & Hammar (2000):
$$CVI = \sqrt{\frac{X_1 \cdot X_2 \cdot X_3 \cdot X_4 \cdot X_5 \cdot X_6}{6}}$$
Trong đó $X_1$ đến $X_6$ là 6 biến số hình thái - động lực gồm: Đặc điểm địa mạo đới bờ ($X_1$), Độ dốc địa hình bờ ($X_2$), Tốc độ dâng mực nước biển ($X_3$), Tốc độ biến động bờ trung bình nhiều năm ($X_4$), Biên độ triều ($X_5$), và Độ cao sóng trung bình ($X_6$). Mỗi biến số được chuẩn hóa thành 5 cấp rủi ro (từ 1: Rất thấp đến 5: Rất cao), từ đó xây dựng biểu đồ đường cong lũy tích để phân cấp mức độ tổn thương thực tế.
Phát hiện đột phá và implications
Những phát hiện then chốt
Thứ nhất, sự đảo chiều lịch sử từ ưu thế bồi tích sang gia tăng xói lở trên toàn tuyến bờ biển STCM. Phân tích chuỗi dữ liệu 48 năm (1965–2013) chỉ ra rằng bức tranh tiến hóa đới bờ có sự phân hóa không gian cực kỳ gay gắt:
- Giai đoạn 1965–1990: Toàn dải bờ biển STCM chiếm ưu thế bồi tụ vượt trội, đặc biệt tại bán đảo Cà Mau với tốc độ tiến ra biển hàng chục mét mỗi năm.
- Giai đoạn 1990–2013: Hiện tượng xói lở bùng phát dữ dội. Nhiều đoạn bờ bồi tụ trước đây đã chuyển hóa hoàn toàn sang pha xói lở nghiêm trọng. Cán cân bồi - xói tại nhiều tiểu vùng rơi vào trạng thái âm, đường bờ bị giật lùi sâu vào đất liền hàng trăm mét, phá hủy trực tiếp các đai rừng phòng hộ nguyên sinh.
Thứ hai, phát hiện mối liên hệ chặt chẽ giữa trường thủy thạch động lực học đáy biển và hiện tượng xói lở bờ ngầm. Kết quả mô phỏng MIKE-21 năm 2012 chứng minh sự biến động địa hình đáy theo chiều thẳng đứng: Trong 4 tháng cuối năm 2012, các trường sóng cực đại trong mùa gió Đông Bắc tạo nên các dòng chảy dọc bờ có vận tốc lớn, cuốn trôi bùn sét đáy ở dải độ sâu từ 0 đến 5 m nước, làm gia tăng độ dốc sườn bờ ngầm, tạo điều kiện cho năng lượng sóng truyền thẳng vào chân đê và chân rừng phòng hộ gây sạt lở hàm ếch.
Thứ ba, sự phân hóa sâu sắc của chỉ số tổn thương CVI dọc theo bờ biển 3 tỉnh. Biểu đồ đường cong lũy tích CVI xác định các khu vực có mức độ tổn thương rất cao (Very High Vulnerability) tập trung tại các đoạn bờ cửa sông Gành Hào, Nhà Mát (Bạc Liêu), Tân Ân, Rạch Gốc, Tam Giang Đông (Cà Mau) và Vĩnh Tân (Sóc Trăng) - nơi ghi nhận sự kết hợp bất lợi giữa nền địa chất bùn sét bở rời, biên độ triều cao trên 3 m và tốc độ giật lùi bờ biển lớn hơn 15 m/năm.
Thứ tư, nhận diện nguyên nhân cốt lõi gây gia tăng xói lở bờ biển là tác động liên hợp giữa sụt giảm trầm tích thượng nguồn sông Mê Công và suy thoái rừng ngập mặn do chuyển đổi nuôi tôm. Dữ liệu đối sánh từ công trình của Goichot & Bravard (2013) cho thấy các bậc thang thủy điện thượng nguồn và hoạt động khai thác cát lòng sông tại hạ lưu Mê Công đã giữ lại phần lớn lượng trầm tích lơ lửng, làm giảm nguồn vật liệu bùn cát cung cấp cho dòng ven biển Tây Nam, kết hợp với việc chặt phá rừng ngập mặn làm đầm nuôi thủy sản nội đồng trong giai đoạn 1990–2005 đã triệt tiêu hoàn toàn khả năng tiêu tán năng lượng sóng tự nhiên.
Implications đa chiều
- Về mặt lý luận: Khẳng định tính tất yếu của việc đưa nghiên cứu địa mạo động lực và địa mạo sinh vật vào quy hoạch không gian biển. Luận án mở rộng phương pháp luận đánh giá biến động bờ biển bùn sét nhiệt đới thông qua tích hợp chuỗi viễn thám đa thời gian với các chỉ thị thực vật.
- Về mặt phương pháp: Cung cấp bộ quy trình chuẩn hóa kết hợp giữa GIS, DSAS, mô hình số trị MIKE-21 và CVI có thể chuyển giao và nhân rộng cho các vùng bờ biển châu thổ khác tại Việt Nam như châu thổ sông Hồng hoặc đồng bằng ven biển Bắc Trung Bộ.
- Về mặt thực tiễn và chính sách: Cung cấp luận cứ khoa học trực tiếp phục vụ triển khai Điều 23 và Điều 79 Luật Tài nguyên, Môi trường biển và hải đảo năm 2015. Kết quả phân vùng quản lý và bản đồ CVI là công cụ hỗ trợ ra quyết định phân bổ ngân sách đầu tư công, chuyển dịch từ đê kè cứng sang các giải pháp công trình mềm bảo vệ bờ (như kè ngầm phá sóng giảm sóng bằng hàng cọc tre hai hàng kết hợp đá hộc rọ hàn, phục hồi rừng ngập mặn).
Limitations và Future Research
Luận án thẳng thắn thừa nhận các giới hạn khoa học nội tại:
- Độ phân giải thời gian của tư liệu viễn thám: Khoảng cách giữa các mốc ảnh (1965–1990–2013) còn khá dài (trên 20 năm mỗi giai đoạn), chưa phản ánh hết được các biến động mang tính chu kỳ mùa hoặc các sự kiện thời tiết cực đoan bất thường (như bão số 5 Linda năm 1997).
- Số liệu trầm tích và độ sâu đo đạc trực tiếp: Độ sâu đáy biển và số liệu bùn cát lơ lửng ngoài khơi độ sâu 10–20 m còn phụ thuộc vào các đề tài đo vẽ trước đó; mạng lưới trạm quan trắc hải văn cố định liên tục tại vùng biển Tây Nam còn thưa thớt.
- Mô hình tính toán 2 chiều (2D): Chưa phản ánh đầy đủ hoàn lưu thẳng đứng 3 chiều (3D) tại các cửa sông lớn (Định An, Trần Đề, Bồ Đề, Cửa Lớn) nơi có sự phân tầng mật độ nước mặn - nước ngọt rõ rệt.
Chương trình nghiên cứu tiếp nối trong tương lai cần tập trung vào:
- Ứng dụng ảnh vệ tinh radar khẩu độ tổng hợp (Sentinel-1 SAR, COSMO-SkyMed) và ảnh quang học phân giải cao (PlanetScope, WorldView) theo chu kỳ tuần/tháng để theo dõi động thái bờ biển thời gian thực.
- Xây dựng mô hình coupled 3 chiều hoàn lưu - sóng - vận chuyển trầm tích (Delft3D 3D hoặc MIKE-3) có tính đến ảnh hưởng của biến động dòng chảy sông Mê Công theo mùa lũ - cạn.
- Nghiên cứu cơ chế tiêu tán năng lượng sóng qua các kiểu cấu trúc rễ rừng ngập mặn (Avicennia alba, Rhizophora apiculata) dưới các kịch bản nước biển dâng RCP 4.5 và RCP 8.5 của Bộ Tài nguyên và Môi trường.
- Đánh giá kinh tế - sinh thái lượng hóa toàn diện giá trị của các dịch vụ địa mạo bờ biển phục vụ chi trả dịch vụ môi trường rừng ngập mặn.
Tác động và ảnh hưởng
Luận án mang lại tác động học thuật và thực tiễn to lớn:
- Tác động học thuật: Đóng góp một công trình nghiên cứu hệ thống, toàn diện đầu tiên về địa mạo động lực dải bờ biển Sóc Trăng - Cà Mau; ước tính tiềm năng trích dẫn học thuật cao trong các nghiên cứu về địa mạo ven biển, hải dương học và biến đổi khí hậu tại Đông Nam Á.
- Tác động kinh tế - xã hội: Giúp tiết kiệm hàng trăm tỷ đồng ngân sách nhà nước thông qua việc định hướng thay thế các công trình kè bê tông kiên cố không phù hợp (vốn dễ bị sụp đổ do xói ngầm chân kè tại Hòn Đá Bạc, Nhà Mát) bằng các giải pháp công trình sinh thái mềm thân thiện môi trường.
- Tác động chính sách: Cung cấp cơ sở khoa học trực tiếp cho UBND 3 tỉnh Sóc Trăng, Bạc Liêu, Cà Mau và Tổng cục Biển và Hải đảo Việt Nam trong việc phê duyệt quy hoạch hành lang bảo vệ bờ biển, quy hoạch sử dụng không gian biển và thiết lập các khu bảo tồn đất ngập nước ven biển.
Đối tượng hưởng lợi
- Nghiên cứu sinh và Giới học thuật: Tiếp cận cơ sở dữ liệu số hóa phong phú về biến động đường bờ 48 năm, khung phương pháp luận kết hợp DSAS - MIKE-21 - CVI và hệ thống 45 chỉ thị địa mạo bờ biển bùn sét.
- Cơ quan Quản lý Nhà nước (Bộ TN&MT, Tổng cục Biển & Hải đảo, Sở TN&MT địa phương): Nhận được bản đồ phân vùng chức năng bờ biển, bản đồ chỉ số tổn thương CVI và tọa độ ranh giới hành lang cảnh báo nguy cơ xói lở - bồi tụ để tích hợp vào đồ án quy hoạch tổng thể đới bờ.
- Doanh nghiệp Xây dựng Thủy lợi & Năng lượng Tái tạo: Các chủ đầu tư dự án điện gió ngoài khơi (như khu vực Nhà Mát, Bạc Liêu, Vĩnh Châu) có số liệu nền tảng về trường sóng, phân bố trầm tích đáy và động thái bồi xói để tính toán an toàn móng tuabin gió.
- Cộng đồng dân cư ven biển: Được bảo vệ an toàn tính mạng và sinh kế, giảm thiểu rủi ro mất nhà cửa, mất đầm tôm và đất sản xuất nông nghiệp nhờ hệ thống cảnh báo sớm và phân vùng sử dụng đất hợp lý.
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 lý thuyết độc đáo nhất là việc xác lập cơ chế phản hồi động lực giữa địa mạo - sinh vật - trầm tích trong hệ thống bờ biển bùn sét nhiệt đới, mở rộng trực tiếp Học thuyết Hệ thống Địa mạo của Chorley (1962) và Học thuyết Cân bằng Động của Hack (1960). Luận án chứng minh rằng trong môi trường đới bờ châu thổ, thảm thực vật rừng ngập mặn không chỉ là lớp phủ thụ động mà là một tác nhân cấu trúc phản ánh trạng thái cân bằng động của hệ thống: Quá trình bồi tụ gắn liền với thế hệ rừng mắm tiên phong lấn biển, trong khi quá trình xói lở biểu hiện qua sự sụp đổ vách triều và đảo lộn cấu trúc sinh thái, đóng vai trò là chỉ thị hình thái định lượng cho sự biến đổi quỹ trầm tích khu vực.
2. Điểm mới về mặt phương pháp so với các nghiên cứu trước đây?
So sánh với các nghiên cứu của Nguyễn Thanh Ngà (1995), Lê Xuân Hồng (1997) hay Nguyễn Văn Cư (2001, 2005) - vốn chủ yếu dựa trên khảo sát trắc địa mặt cắt truyền thống hoặc phân tích định tính tại các điểm nóng xói lở riêng lẻ, luận án đã thực hiện một bước đột phá về công nghệ tích hợp:
- Chuẩn hóa chuỗi dữ liệu không gian viễn thám kéo dài 48 năm (1965–2013) bằng công cụ DSAS của USGS để tính toán tốc độ EPR và NSM trên hàng nghìn mặt cắt tự động.
- Liên kết kết quả viễn thám với mô hình số trị thủy thạch động lực MIKE-21 (HD, SW, MT) để giải thích cơ chế bồi xói đáy biển 2 chiều theo các trường sóng thực tế.
- Tích hợp kết quả vào chỉ số CVI và mô hình tính toán độ rộng vùng nhạy cảm xói lở EPAW theo chuẩn mực quốc tế của bang Queensland (Úc) và Liên minh Châu Âu.
3. Phát hiện bất ngờ nhất từ dữ liệu thực nghiệm là gì?
Phát hiện bất ngờ và mang tính nghịch lý nhất là hiện tượng xói lở nghiêm trọng diễn ra khốc liệt tại chính những khu vực vốn có truyền thống bồi tụ mạnh nhất trong lịch sử như bãi Khai Long, Mũi Bà Quan và khu vực phía Tây Rạch Gốc thuộc Mũi Cà Mau. Dữ liệu thực địa và mô hình hóa chỉ ra rằng mặc dù phù sa sông Mê Công vẫn tiếp tục đổ ra biển, nhưng sự thay đổi cán cân trầm tích lơ lửng kết hợp với sự gia tăng cục bộ năng lượng sóng gió mùa Đông Bắc khúc xạ quanh bán đảo Cà Mau đã làm xói sâu sườn bờ ngầm, biến bãi bồi nông thành các vách xói lở dạng răng cưa sâu hoắm, cuốn trôi hàng loạt cây đước, cây mắm cổ thụ có tuổi thọ hàng chục năm.
4. Luận án có cung cấp quy trình tái lập nghiên cứu (Replication Protocol) rõ ràng không?
Có. Luận án cung cấp một quy trình nghiên cứu 3 bước hoàn chỉnh và minh bạch:
- Các thông số ảnh vệ tinh Landsat (chu kỳ, mã cảnh 125053, 125054, 126053, 126054, thời điểm tháng 8) và SPOT được mô tả chính xác.
- Quy trình tạo Geodatabase, thiết lập Baseline, khoảng cách mặt cắt trong phần mềm DSAS được quy chuẩn hóa chi tiết.
- Hệ phương trình vi phân đạo hàm riêng, điều kiện biên, hệ số nhám Manning, độ nhớt xoáy và các tenxơ ứng suất bức xạ trong mô hình MIKE-21 được thiết lập đầy đủ, cho phép các nhà khoa học khác tái lập và kiểm chứng mô hình trên bất kỳ đoạn bờ biển nào.
5. Khung chương trình nghiên cứu 10 năm được vạch ra như thế nào?
Khung nghiên cứu 10 năm được vạch ra tập trung vào 4 trọng tâm:
- Xây dựng trạm quan trắc địa mạo - hải văn tự động ven bờ tích hợp camera độ phân giải cao và thiết bị đo sâu hồi âm đa chùm (Multi-beam).
- Phát triển mô hình dự báo xói lở bờ biển thời gian thực (Real-time Coastal Erosion Forecasting System) ứng dụng trí tuệ nhân tạo (AI) và dữ liệu vệ tinh Radar.
- Nghiên cứu sâu cơ chế biến động đường bờ dưới tác động của các kịch bản xây dựng chuỗi đập thủy điện dòng chính sông Mê Công tại Lào và Campuchia.
- Thử nghiệm và đánh giá hiệu quả dài hạn của các mô hình kè rọ đá giảm sóng kết hợp trồng rừng ngập mặn phục hồi bờ biển bền vững tại các điểm nóng xói lở Tây Nam Bộ.
Kết luận
Luận án "Nghiên cứu đánh giá biến động đường bờ phục vụ quản lý bờ biển Sóc Trăng - Cà Mau" của NCS. Lưu Thành Trung là một công trình khoa học công phu, nghiêm túc, có tính mới và giá trị ứng dụng thực tiễn vượt trội. Luận án đúc kết 6 đóng góp cốt lõi:
- Xác lập bức tranh toàn diện, định lượng về hiện trạng và xu thế biến động đường bờ dải ven biển Sóc Trăng - Cà Mau trong suốt chiều dài 48 năm (1965–2013), chứng minh sự đảo chiều lịch sử từ ưu thế bồi tụ sang gia tăng xói lở cục bộ sâu sắc.
- Vận dụng thành công hệ thống 45 chỉ thị địa mạo - sinh vật kết hợp viễn thám đa thời gian để nhận diện chính xác ranh giới bờ biển bùn sét ngập mặn.
- Ứng dụng thành công mô hình số trị MIKE-21 mô phỏng trường sóng cực đại, dòng chảy ven bờ và động lực bồi xói đáy biển theo chiều thẳng đứng giai đoạn 2009–2012.
- Xây dựng bản đồ phân cấp tính dễ bị tổn thương bờ biển (CVI) tích hợp 6 tham số hình thái - động lực cho toàn bộ bờ biển 3 tỉnh Sóc Trăng, Bạc Liêu, Cà Mau.
- Thiết lập thành công hành lang cảnh báo nguy cơ xói lở - bồi tụ (EPAW) tại 3 xã trọng điểm ven biển phía Nam Mũi Cà Mau (Nguyễn Huân, Tam Giang Đông, Tam Giang Tây).
- Đề xuất hệ thống giải pháp định hướng quản lý đới bờ theo tiểu vùng địa mạo, gắn kết chặt chẽ giữa khoa học địa mạo ứng dụng với việc triển khai Luật Tài nguyên, Môi trường biển và hải đảo năm 2015.
Công trình không chỉ mở ra hướng nghiên cứu mới về Địa mạo Môi trường và Địa lý Biển tại Việt Nam mà còn là tài liệu tham khảo mẫu mực cho công tác quy hoạch, quản lý tổng hợp và bảo vệ bền vững đới bờ biển quốc gia trước thách thức của biến đổi khí hậu toàn cầu.