Tổng quan về luận án
Nghiên cứu về sự tiến hóa của các đồng bằng châu thổ ven biển trong kỷ Đệ tứ, đặc biệt là giai đoạn Holocen, đóng vai trò nền tảng trong khoa học Địa chất và Địa mạo ven bờ. Chiếm giữ vị trí địa kinh tế - sinh thái chiến lược, đồng bằng châu thổ sông Mã là hệ thống châu thổ lớn thứ ba tại Việt Nam (xếp sau châu thổ sông Cửu Long và châu thổ sông Hồng). Khu vực ven biển Thanh Hóa có mạng lưới thủy văn gồm 90 nhánh sông bắt nguồn từ Điện Biên, chảy qua Lào và đổ ra Biển Đông qua ba cửa sông chính: cửa Lèn, cửa Lạch Trường và cửa Hới (trong đó cửa Hới có quy mô động lực lớn nhất). Trầm tích Holocen phân bố rộng khắp lãnh thổ ven biển Thanh Hóa với bề dày trung bình từ 15 đến 20m; tại khu vực lân cận cửa Hới và cửa Lèn, chiều dày đạt từ 30 đến 32m, và đạt cực đại lên đến 41m tại cửa Lạch Trường.
Khoảng trống nghiên cứu (research gap) trọng tâm nằm ở sự thiếu hụt các mô hình địa tầng phân tập (sequence stratigraphy) phân giải cao tích hợp đa phương pháp (multi-proxy) để làm sáng tỏ cơ chế tiến hóa trầm tích chịu sự chi phối áp đảo của động lực sóng (wave-dominated delta) trong bối cảnh tương tác phức tạp giữa chuyển động kiến tạo cục bộ dọc hệ thống đứt gãy sông Mã và dao động mực nước biển chân tĩnh (eustatic sea-level change) sau băng hà cuối cùng. Các nghiên cứu trước đây tại khu vực chủ yếu dừng lại ở công tác đo vẽ bản đồ địa chất khu vực tỷ lệ 1:200.000 (Đặng Trần Quân, 1980; Lê Duy Bách, 1995) hoặc điều tra địa chất đô thị (Hoàng Văn Khổn, 1997), chưa giải quyết thấu đáo cấu trúc không gian - thời gian của các thung lũng cắt xẻ (incised valley) và quy luật phân dị tướng trầm tích chi tiết.
Hệ thống câu hỏi nghiên cứu và giả thuyết khoa học được thiết lập chặt chẽ:
- RQ1: Cấu trúc tướng và môi trường lắng đọng trầm tích Holocen khu vực ven biển Thanh Hóa được phân hóa như thế nào dưới tác động tương hỗ giữa động lực sông, sóng và thủy triều?
- Giả thuyết H1: Trầm tích Holocen khu vực ven biển Thanh Hóa bao gồm tổ hợp 23 tướng trầm tích đơn lẻ thuộc 4 nhóm tướng chính, phản ánh sự chuyển biến năng lượng thủy động lực từ lục địa sang vũng vịnh/estuary và tiến dần sang châu thổ - bồi tích aluvi hiện đại.
- RQ2: Quy luật dao động mực nước biển Flandrian và khung địa tầng phân tập chi phối không gian tích tụ (accommodation space) của vùng cửa sông ven biển Thanh Hóa diễn ra theo những pha nào?
- Giả thuyết H2: Quá trình biển tiến Flandrian khởi đầu từ cuối Pleistocen muộn lấp đầy các thung lũng cắt xẻ, đạt đỉnh ngập lụt cực đại (MFS) vào đầu Holocen giữa (+4m đến +5m), tạo lập hệ thống trầm tích biển tiến (TST) và kế tiếp là hệ thống trầm tích biển cao (HST) dạng nêm lấn (progradation).
- RQ3: Động lực sóng thể hiện vai trò thống trị như thế nào trong việc định hình hình thái học đới bờ và sự bất đối xứng của đồng bằng châu thổ sông Mã trong Holocen muộn?
- Giả thuyết H3: Sự chiếm ưu thế của trường sóng biển Đông và dòng vận chuyển trầm tích dọc bờ (longshore drift) đã tạo nên các phức hệ cồn cát ven biển, đầm lầy sau bar và làm biến tính hình thái cửa sông theo mô hình châu thổ sóng điển hình.
Khung lý thuyết của luận án được xây dựng dựa trên sự tích hợp giữa: Lý thuyết địa tầng phân tập hiện đại (Posamentier et al., 1988; Van Wagoner et al., 1990; Catuneanu, 2006, 2009), Định luật tương quan tướng Walther (1894), Mô hình phân loại động lực tam giác cửa sông - châu thổ (Galloway, 1975; Coleman & Wright, 1975; Boyd et al., 2006) và Lý thuyết chỉ số bất đối xứng châu thổ chịu tác động sóng (Bhattacharya & Giosan, 2003).
Phạm vi không gian nghiên cứu trải dài từ huyện Nga Sơn đến hết huyện Quảng Xương, vươn sâu vào đất liền từ 15 đến 20km, giới hạn trong tọa độ 105°45'02" đến 106°01'02" kinh độ Đông và 19°37'10" đến 19°59'53" vĩ độ Bắc. Cơ sở dữ liệu thực nghiệm bao gồm: 09 lỗ khoan máy lấy mẫu nguyên dạng (LKTH1 đến LKTH9) với 213m lõi khoan, 531 mẫu phân tích thành phần độ hạt, 25 mẫu lát mỏng thạch học bở rời, 186 mẫu bào tử - phấn hoa, 186 mẫu tảo Diatomeae, 140 mẫu Foraminifera, 10 mẫu hóa - lý môi trường, 10 mẫu khoáng vật sét XRD, 21 mẫu định tuổi tuyệt đối bằng carbon phóng xạ AMS $^{14}\text{C}$, 11,7 km tuyến địa chấn nông phân giải cao ven bờ, 1,3 km tuyến địa chấn khúc xạ và 107 tài liệu địa tầng lỗ khoan kế thừa.
Literature Review và Positioning
Lịch sử nghiên cứu trầm tích châu thổ trên thế giới ghi nhận những mốc phát triển mang tính cách mạng từ quan sát kinh điển của Herodotus (400 TCN) tại sông Nile, đến các công trình định lượng hóa tại châu thổ sông Mississippi của Barrell (1912), Fisk (1944, 1955, 1961), Coleman & Wright (1975). Galloway (1975) đã đặt nền móng cho việc phân loại châu thổ thành ba cực động lực: sông (river-dominated), sóng (wave-dominated) và triều (tide-dominated). Tiếp tục mở rộng khung lý thuyết này, Bhattacharya & Giosan (2003) đưa ra "Chỉ số bất đối xứng A" biểu thị tỷ số giữa lưu lượng vận chuyển trầm tích dọc bờ với lưu lượng nước sông để giải thích cấu trúc bất đối xứng của các châu thổ Danube hay Ebro. Ron Boyd et al. (2006) và Siddiqui et al. (2017) đã hoàn thiện biểu đồ tam giác phân định 6 kiểu môi trường lắng đọng ven biển: đồng bằng châu thổ, estuary do sóng thống trị, đầm lầy sau bar, đồng bằng cồn cát (strand plain), estuary do triều thống trị và đồng bằng triều (tidal flat).
ĐỘNG LỰC SÔNG (River)
/\
/ \
/ 1 \
/ \
/ 2 5 \
/ 3 4 \
/____________\
ĐỘNG LỰC SÓNG (Wave) ĐỘNG LỰC TRIỀU (Tide)
(Mô hình 6 môi trường trầm tích ven bờ theo Boyd et al., 2006: 1-Châu thổ sông; 2-Estuary do sóng; 3-Đầm lầy sau bar; 4-Đồng bằng cồn cát; 5-Estuary do triều; 6-Đồng bằng triều)
Tại Việt Nam, các nghiên cứu tiến hóa châu thổ Holocen tập trung sâu sắc tại hai bồn trũng lớn: châu thổ sông Hồng (Trần Nghi, 2000, 2017; Doãn Đình Lâm, 2004, 2005; Tanabe et al., 2004, 2006; Phùng Văn Phách et al., 2018; Nguyễn Tài Tuệ, 2019) và châu thổ sông Cửu Long (Lê Đức An, 1996; Nguyễn Địch Dỹ et al., 2010; Tạ Thị Kim Oanh, Nguyễn Văn Lập, 2004; Nguyễn Thị Thu Cúc, 2015; Vũ Văn Hà, 2015). Trong khi châu thổ sông Hồng và sông Cửu Long thể hiện sự chuyển dịch phức hợp giữa chế độ triều và sông, thì các đồng bằng duyên hải miền Trung, tiêu biểu là đồng bằng Thanh Hóa, lại phản ánh sắc nét vai trò kiểm soát của động lực sóng và đới bờ hẹp.
Trong y văn tồn tại hai cuộc tranh luận học thuật sâu sắc:
- Biên độ và thời điểm đỉnh biển tiến Flandrian: Một trường phái cho rằng mực nước biển dâng cao vượt mức hiện tại từ +5m đến +6m trong khoảng 6.000 - 5.000 năm Bp (Fairbridge; Pirazzoli, 1991; Tanabe et al., 2006; Trần Nghi, 2017; Nguyễn Biểu et al., 2005); trong khi Shepard và cộng sự cho rằng mực nước biển dâng sau băng hà không bao giờ vượt quá độ cao hiện đại, hoặc chỉ đạt đỉnh cục bộ ở mức +2m đến +3m (Lưu Tỳ, Nguyễn Thế Tiệp, 1985; Phùng Văn Phách, 2009).
- Cơ chế khống chế hình thái châu thổ sông Mã: Luận điểm truyền thống xem châu thổ sông Mã là dạng châu thổ sông thuần túy (river-dominated delta) được hình thành đơn thuần do bồi tích phù sa (Đặng Trần Quân, 1980); trường phái đối lập cho rằng đây là một hệ thống cửa sông đứt gãy bị tái biến tính mạnh mẽ bởi sóng biển và dòng ven bờ (Lưu Đức Dũng, 2017).
Luận án này định vị chính xác khoảng trống khoa học bằng cách so sánh đối chiếu trực tiếp với các mô hình quốc tế như châu thổ Danube (Romania/Ukraine - Giosan et al., 2005) và châu thổ Tiber/Ebro (Địa Trung Hải - Bhattacharya & Giosan, 2003), chứng minh rằng châu thổ sông Mã thời kỳ Holocen muộn đại diện cho mô hình châu thổ sóng bất đối xứng (asymmetric wave-dominated delta), trong đó các cồn cát ven biển (sand ridges/dunes) và đồng bằng lagoon đóng vai trò che chắn cơ học cho hệ thống bồi tích nội đồng.
Đó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 và kiểm chứng thực nghiệm mô hình địa tầng phân tập của Posamentier (1988), Van Wagoner (1990) và Catuneanu (2006, 2009) trong môi trường châu thổ nhiệt đới gió mùa có biên độ triều nhỏ và năng lượng sóng mạnh. Các đóng góp lý thuyết cốt lõi gồm:
+-------------------------------------------------------------------------+
| KHUNG TIẾN HÓA ĐỊA TẦNG PHÂN TẬP |
+-------------------------------------------------------------------------+
| Pha 4: HST Aluvi (Hiện đại) |
| -> Cát lòng sông, đê ven sông, đồng bằng ngập lụt |
+-------------------------------------------------------------------------+
| Pha 3: HST Châu thổ (Đầu Holocen giữa - Holocen muộn: ~6.0 - 1.0 ka Bp) |
| -> Chân châu thổ -> Tiền châu thổ -> Bãi triều/Cồn cát/Lagoon |
+-------------------------------------------------------------------------+
| === BỀ MẶT NGẬP LỤT CỰC ĐẠI (MFS) (~6.0 ka Bp: Đỉnh biển tiến +4 đến +5m)==|
+-------------------------------------------------------------------------+
| Pha 2: TST Estuary - Vũng vịnh (Giữa Holocen sớm - Đầu Holocen giữa) |
| -> Sét bột vũng vịnh, cồn cát triều, lạch triều |
+-------------------------------------------------------------------------+
| === BỀ MẶT BIỂN TIẾN (TS) / BẤT CHỈNH HỢP ĐÁY (SB) ======================|
+-------------------------------------------------------------------------+
| Pha 1: TST Bồi lấp thung lũng cắt xẻ (Cuối Pleistocen - Holocen sớm) |
| -> Cát sạn sỏi lòng sông cổ, sét bột đầm lầy lấp đầy TLCX |
+-------------------------------------------------------------------------+
- Mệnh đề P1: Trong giai đoạn biển tiến nhanh (tốc độ 9 - 12 mm/năm ở đầu Holocen sớm), không gian tích tụ ($A$) tăng nhanh hơn tốc độ cung cấp trầm tích ($S$) ($A > S$), tạo điều kiện cho hệ thống estuary - vũng vịnh lấn sâu vào đất liền, bồi lấp mặt bóc mòn của thung lũng cắt xẻ Pleistocen muộn.
- Mệnh đề P2: Bề mặt ngập lụt cực đại (MFS) được thiết lập tại thời điểm khoảng 6.000 năm Bp khi mực nước biển dâng cực đại (+4m đến +5m), đánh dấu bước ngoặt chuyển dịch từ hệ thống trầm tích biển tiến (TST) sang hệ thống trầm tích biển cao (HST).
- Mệnh đề P3: Giai đoạn biển lùi tương đối từ sau 6.000 năm Bp diễn ra theo dạng dao động con lắc tắt dần, tốc độ hạ thấp mực nước biển giảm dần theo hàm sin kết hợp tốc độ hạ lún kiến tạo khu vực (~2 mm/năm), khiến tốc độ bồi tụ vượt không gian tích tụ ($S > A$), thúc đẩy châu thổ tiến ra biển (delta progradation) dưới sự tái phân bố trầm tích của sóng.
Khung phân tích độc đáo
Khung phân tích của luận án kết hợp đa chiều giữa Thạch học trầm tích định lượng, Cổ sinh thái học (Biofacies), Địa hóa môi trường và Địa chấn - Địa tầng:
- Liên kết cấu trúc trầm tích theo Định luật Walther (1894): Tái lập trật tự không gian của 23 tướng trầm tích đơn lẻ trên các mặt cắt thẳng đứng qua 9 lỗ khoan máy đồng bộ.
- Định lượng môi trường qua chỉ số Kation trao đổi ($Kt$) và Sulfua sắt ($Fe^{2+}S/C{org}$): Ứng dụng công thức của Grim (1974) và hệ chỉ tiêu Trần Nghi (2002):
$$Kt = \frac{\text{Na}^+ + \text{K}^+}{\text{Ca}^{2+} + \text{Mg}^{2+}}$$
Xác định rõ ba ngưỡng môi trường: Lục địa ($Kt < 0,5$; $Fe^{2+}S/C{org} < 0,06$), Chuyển tiếp/Lợ ($0,5 \le Kt \le 1,0$; $0,06 \le Fe^{2+}S/C{org} \le 0,2$) và Biển mở ($Kt > 1,0$; $Fe^{2+}S/C{org} > 0,2$).
- Phân tích quần xã chỉ thị sinh thái: Phối hợp định danh giữa Diatomeae (tảo Silic nước mặn, lợ, ngọt bám đáy và trôi nổi) và Foraminifera (trùng lỗ đáy biển nông) cùng phức hệ Bào tử - phấn hoa đại diện cho thảm thực vật ngập mặn (Rhizophora, Sonneratia) và đầm lầy lục địa.
Phương pháp nghiên cứu tiên tiến
Thiết kế nghiên cứu
Luận án tuân thủ thế giới quan khoa học thực chứng luận (positivism) kết hợp trường phái hiện thực phản tư (critical realism). Thiết kế nghiên cứu tích hợp phương pháp hỗn hợp (mixed-methods) ở cấp độ cao, kết nối dữ liệu địa vật lý trường sâu (geophysics), hình thái trắc địa viễn thám không gian (remote sensing) và dữ liệu lõi khoan phân tích hóa - lý - sinh siêu vi trong phòng thí nghiệm.
+-------------------------------------------------------------------------+
| SƠ ĐỒ THIẾT KẾ NGHIÊN CỨU TÍCH HỢP |
+-------------------------------------------------------------------------+
| [Thực địa & Khoan sâu] -> 09 lỗ khoan (LKTH1-9, 213m lõi nguyên dạng) |
| -> Bóc tách bề mặt phong hóa Pleistocen muộn |
+-------------------------------------------------------------------------+
| [Địa vật lý & Viễn thám] -> Địa chấn nông ven bờ (11,7 km, sâu -10m nước)|
| -> Địa chấn khúc xạ ABEM TerralogPro (1,3 km) |
| -> Chuỗi ảnh Landsat (1975-2017) & Sentinel |
+-------------------------------------------------------------------------+
| [Phân tích Lab] -> 531 mẫu độ hạt (Folk, 1954: Md, So, Sk) |
| -> 25 lát mỏng thạch học (kính hiển vi phân cực)|
| -> 10 mẫu sét XRD (PANAlytical CuKα 45kV/40mA) |
| -> 10 mẫu hóa - lý (pH, Eh, Kt, Fe2+S/Corg) |
| -> 186 Diatomeae + 186 Bào tử-phấn + 140 Foram |
| -> 21 mẫu tuổi tuyệt đối AMS 14C (DirectAMS) |
+-------------------------------------------------------------------------+
| [Mô hình hóa] -> Khung Địa tầng phân tập (SB, TS, TST, MFS, HST)|
| -> Bản đồ Tướng đá - Cổ địa lý 4 giai đoạn |
+-------------------------------------------------------------------------+
Quy trình nghiên cứu rigorous
- Quy trình khoan và bảo quản lõi: Thực hiện 09 lỗ khoan máy (LKTH1 đến LKTH9), khoan qua toàn bộ tầng trầm tích Holocen và khoan cắm sâu vào bề mặt phong hóa loang lổ thuộc hệ tầng Vĩnh Phúc ($Q_1^{3}vp$) tuổi Pleistocen muộn. Lấy 100% mẫu lõi nguyên dạng với tổng chiều dài 213m.
- Quy trình địa vật lý địa chấn khúc xạ: Sử dụng máy ABEM TerralogPro (Thụy Điển) 24 kênh. Mỗi trạm máy bố trí 115m tuyến cáp vật lý với 7 điểm nổ chập ($DN_{-55m}$, $DN_{0m}$, $DN_{27,5m}$, $DN_{57,5m}$, $DN_{82,5m}$, $DN_{115m}$, $DN_{170m}$), thực hiện thu nhận 2 biểu đồ thời khoảng xuôi và 2 biểu đồ thời khoảng ngược để dựng mặt cắt ranh giới âm học đáy thung lũng cắt xẻ tại lỗ khoan LKTH9.
- Quy trình kiểm soát tính giá trị (validity) và độ tin cậy (reliability): Kiểm soát chéo tam giác (triangulation): Dữ liệu mặt cắt địa chấn phân giải cao tuyến NB-PX95 và tuyến DVL-21 được chuẩn hóa đối sánh trực tiếp với cột thạch học lỗ khoan LKTH9; kết quả phân tích vi cổ sinh được giám định độc lập trên hai hệ kính Carl Zeiss và Nikon Eclipse.
Data và phân tích
- Phân tích độ hạt: 531 mẫu trầm tích bở rời được xử lý bằng hệ thống rây cơ học và phương pháp ống hút Pipette, tính toán các tham số thống kê độ hạt theo Folk (1954): kích thước hạt trung vị ($Md$), độ chọn lọc ($So$), hệ số bất đối xứng ($Sk$).
- Nhiễu xạ tia X (XRD): 10 mẫu tổng và mẫu sét định hướng phân tích trên thiết bị nhiễu xạ tia X PANalytical, bức xạ CuKα, điện áp 45 kV, dòng điện 40 mA, bước nhảy quét $0,02^\circ$. Xác định các khoáng vật sét ưu thế: Kaolinit (chỉ thị phong hóa nóng ẩm lục địa), Smectit (môi trường kiềm, khí hậu khô nóng), Illit và Clorit.
- Định tuổi carbon phóng xạ AMS $^{14}\text{C}$: 21 mẫu hữu cơ, vỏ sò ốc và than bùn được phân tích bằng phương pháp khối phổ gia tốc (AMS) tại Phòng thí nghiệm DirectAMS (Hoa Kỳ). Dữ liệu tuổi được hiệu chỉnh chuẩn hóa (calibrated years BP) bằng phần mềm Calib Radiocarbon Calibration Program phiên bản Rev 8.0 (Stuiver & Reimer).
Phát hiện đột phá và implications
Những phát hiện then chốt
- Xác lập cấu trúc và quy mô thung lũng cắt xẻ (TLCX) Pleistocen muộn: Luận án đã làm sáng tỏ hình thái học của mạng lưới lòng sông cổ bị xâm thực sâu trong kỳ bằng hà cực đại (LGM - Last Glacial Maximum) khi mực nước biển hạ thấp xuống -100m đến -120m. Tại LKTH9, đáy thung lũng cắt xẻ nằm ở độ sâu -32m, lấp đầy bởi cuội sỏi, cát thô lòng sông thuộc hệ thống trầm tích biển tiến (TST), phủ bất chỉnh hợp rõ rệt lên bề mặt phong hóa loang lổ chứa kết vón laterit của hệ tầng Vĩnh Phúc.
- Khám phá và phân loại chi tiết 23 tướng trầm tích Holocen: Luận án đã thiết lập 23 tướng trầm tích đơn lẻ thuộc 4 nhóm tướng tiến hóa có quy luật:
- Nhóm tướng aluvi lấp đầy TLCX (3 tướng): Cát sạn sỏi lòng sông cổ ($a_1Q_2^{1-2}$); sét bột đồng bằng ngập lụt ($a_2Q_2^{1-2}$); sét bột đầm lầy cửa sông ven biển ($amQ_2^{1-2}$).
- Nhóm tướng estuary - vũng vịnh (5 tướng): Sét bột bãi triều vũng vịnh ($mQ_2^{1-2}$); cồn cát triều ($mQ_2^{1-2}$); sét bột nhánh triều và lạch triều ($tmQ_2^{1-2}$); sét, bột bãi triều ($tmQ_2^{1-2}$) và sét bột vũng vịnh ($mbQ_2^{1-2}$).
- Nhóm tướng châu thổ (10 tướng): Sét, sét bột chân châu thổ ($mQ_2^{2-3}$); bột cát tiền châu thổ ($maQ_2^{2-3}$); sét bột cát bãi triều ($tmQ_2^{2-3}$); cát bột đồng bằng cát ($mvQ_2^{2-3}$); cát cồn cát ($vQ_2^{2-3}$); sét bột trũng sau cồn cát ($bmQ_2^{2-3}$); sét bột đồng bằng lagoon ($lbQ_2^{2-3}$); sét bột bãi trên triều ($tm_1Q_2^{2-3}$); sét bột cát bãi gian triều ($tm_2Q_2^{2-3}$) và cát bột sét bãi dưới triều ($tm_3Q_2^{2-3}$).
- Nhóm tướng aluvi hiện đại (5 tướng): Cát lòng sông ($a_1Q_2^3$); cát bột đê ven sông ($a_2Q_2^3$); sét bột đồng bằng ngập lụt ($a_3Q_2^3$); cát bột bãi ven lòng ($a_4Q_2^3$) và cát bột bãi giữa lòng ($a_5Q_2^3$).
MẶT CẮT TIẾN HÓA TRẦM TÍCH ĐIỂN HÌNH TẠI LỖ KHOAN LKTH6
--------------------------------------------------------------------------
Độ sâu (m) Tướng trầm tích Môi trường & Hệ thống phân tập
--------------------------------------------------------------------------
0,0 - 4,5m Cát bột bãi ven lòng / Đê Aluvi hiện đại (HST)
4,5 - 9,5m Cát cồn cát / Đồng bằng cát Đồng bằng châu thổ sóng (HST)
9,5 - 12,2m Cát bột tiền châu thổ Tiền châu thổ (HST)
12,2 - 15,3m Sét bột chân châu thổ Chân châu thổ (HST)
--- 15,3m -- [ BỀ MẶT NGẬP LỤT CỰC ĐẠI - MFS (Tuổi ~6.000 cal Bp) ] ------
15,3 - 24,8m Sét bột vũng vịnh / Bãi triều Estuary - Vũng vịnh (TST)
24,8 - 28,6m Bùn đầm lầy xen mùn thực vật Cửa sông ven biển (TST)
28,6 - 32,0m Cát sạn sỏi lót đáy TLCX Aluvi lấp đầy TLCX (TST)
--- 32,0m -- [ BẤT CHỈNH HỢP RANH GIỚI TẬP SB (Bề mặt Vĩnh Phúc) ] --------
> 32,0m Sét loang lổ phong hóa Pleistocen muộn (Q1-3 vp)
- Làm sáng tỏ diễn biến dao động mực nước biển Flandrian tại Bắc Trung Bộ: Kết quả 21 mẫu định tuổi AMS $^{14}\text{C}$ xác thực rằng mực nước biển bắt đầu dâng tràn ngập thung lũng sông Mã vào khoảng 10.000 năm Bp với tốc độ 9 - 12 mm/năm. Mực nước biển đạt đỉnh cực đại ở độ cao +4m đến +5m vào khoảng 6.000 năm Bp (giai đoạn văn hóa Đa Bút). Sau đó, mực nước rút dần theo chu kỳ dao động sóng tắt dần: hạ xuống mức +1,5m đến +3m vào khoảng 4.000 - 3.500 năm Bp (giai đoạn văn hóa Hoa Lộc) và tiến dần về vị trí đường bờ hiện đại.
- Chứng minh sự thống trị của động lực sóng trong Holocen muộn: Luận án chỉ ra rằng sau đỉnh biển tiến, sự suy giảm không gian tích tụ cùng với trường sóng Đông Bắc mạnh mẽ trong vịnh Bắc Bộ đã biến cải toàn diện bờ biển Thanh Hóa. Trầm tích cát hạt trung - mịn bị sóng biển bồi tụ tạo thành các dải đê cát, cồn cát chắn bờ cao từ 3 - 5m đến 10m, hình thành đồng bằng lagoon và đồng bằng cồn cát dạng chuỗi (beach ridge strand plains) chạy song song với đường bờ cổ.
- Định lượng biến động đường bờ hiện đại (1965 - 2017): Dữ liệu viễn thám đa thời gian cho thấy sự phân hóa cực đoan: cửa Lèn và cửa Lạch Trường chịu sự xói lở dữ dội tại các mũi nhô và bồi tích nhanh tại các bãi triều thứ cấp; cửa Hới có tốc độ tiến bờ đạt từ 15 đến 25 m/năm ở các giai đoạn bồi tụ mạnh, nhưng hiện nay đang chuyển sang trạng thái xói lở bờ cục bộ do suy giảm phù sa thượng nguồn.
Implications đa chiều
- Đóng góp lý thuyết: Bổ sung luận cứ vững chắc cho mô hình tiến hóa châu thổ bất đối xứng do sóng thống trị tại khu vực Đông Nam Á, chứng minh vai trò chuyển tiếp từ bồn trũng estuary mở sang châu thổ nêm lấn trong bối cảnh biên độ triều nhỏ (<2m).
- Đổi mới phương pháp luận: Chuẩn hóa quy trình liên kết đa chỉ số: kết hợp địa chấn khúc xạ mặt đất để định vị thung lũng cắt xẻ, dùng chỉ số địa hóa $Kt$ và $Fe^{2+}S/C{org}$ để phân định ranh giới nước lợ - biển, và dùng Diatomeae để định lượng chính xác độ mặn cổ môi trường.
- Ứng dụng thực tiễn: Cung cấp bản đồ phân vùng tướng đá và địa kỹ thuật phục vụ quy hoạch tuyến đường bộ ven biển Bắc - Nam đoạn qua Thanh Hóa dài 107,5 km; xác định các tầng chứa nước nhạt trong trầm tích cát cồn cổ ($vQ_2^{2-3}$) và cảnh báo nguy cơ sụt lún công trình khi xây dựng trên tầng bùn sét vũng vịnh/lagoon ($lbQ_2^{2-3}$, $mbQ_2^{1-2}$) có độ chịu tải yếu.
- Khuyến nghị chính sách: Đề xuất giải pháp công trình và phi công trình bảo vệ bờ biển: nghiêm cấm khai thác cát trái phép tại các đê cát ven biển Sầm Sơn, Hoằng Hóa, Quảng Xương; phục hồi dải rừng ngập mặn trên bãi gian triều Nga Sơn - Nga Tân để giảm xung lực sóng và chủ động ứng phó với kịch bản nước biển dâng do biến đổi khí hậu.
Limitations và Future Research
- Các giới hạn học thuật:
- Độ phủ địa chấn nông ven biển: Tuyến địa chấn nông phân giải cao 11,7 km mới chỉ thực hiện ở độ sâu -10m nước ngoài khơi, chưa thể bao phủ toàn diện dải biến tính ven bờ nông (<5m) do hạn chế của mớn nước tàu khảo sát.
- Phương pháp định tuổi cồn cát gió: Các thể cồn cát eolian trên bề mặt chủ yếu được gián tiếp định tuổi qua tương quan địa tầng và khảo cổ học ($^{14}\text{C}$ mùn thực vật/than bùn lân cận), chưa áp dụng phương pháp định tuổi quang phát quang (OSL - Optically Stimulated Luminescence) trực tiếp trên hạt thạch anh.
- Mô hình hóa số trị thủy động lực (Hydrodynamic numerical modeling): Chưa tích hợp các phần mềm mô phỏng 3D dòng chảy - sóng - bùn cát (như Delft3D hoặc MIKE 21) để tái hiện định lượng trường sóng cổ và dòng vận chuyển dọc bờ theo từng mùa gió mùa.
- Chương trình nghiên cứu tương lai (5 hướng cụ thể):
- Ứng dụng kỹ thuật định tuổi OSL trên toàn bộ các chuỗi cồn cát từ Nga Sơn đến Tĩnh Gia để xác lập niên biểu chi tiết của các đợt gió mùa cổ (paleo-monsoon).
- Mô phỏng số trị 3D (Delft3D) tương tác thủy - hải văn cổ để lượng hóa chỉ số bất đối xứng $A$ của châu thổ sông Mã qua từng mốc thời gian 6.000, 4.000 và 2.000 năm Bp.
- Mở rộng mạng lưới khoan sâu và địa chấn phản xạ phân giải cao ra toàn bộ thềm lục địa vịnh Bắc Bộ để kết nối hoàn chỉnh hệ thống thung lũng cắt xẻ sông Mã với sông Hồng.
- Phân tích đồng vị bền $\delta^{13}\text{C}$, $\delta^{18}\text{O}$ và tỷ số $\text{C}/\text{N}$ trên các lõi khoan trầm tích để giải mã chi tiết các biến cố khí hậu cực đoan (sự kiện khô hạn 4.2 ka và 8.2 ka).
- Đánh giá tác động tích lũy của hệ thống các đập thủy điện bậc thang thượng nguồn sông Mã và sông Chu đến sự mất cân bằng cán cân bùn cát đới bờ hiện đại.
Tác động và ảnh hưởng
- Tác động học thuật: Luận án thiết lập một hệ quy chiếu trầm tích luận chuẩn mực cho vùng ven biển Bắc Trung Bộ, có tiềm năng trích dẫn cao trong các nghiên cứu quốc tế về tiến hóa châu thổ Đông Nam Á và biến đổi môi trường Đệ tứ.
- Chuyển dịch công nghiệp & Xây dựng hạ tầng: Dữ liệu tướng trầm tích cung cấp cơ sở địa chất công trình thiết yếu cho các tổ hợp đô thị ven biển, khu du lịch sinh thái FLC Sầm Sơn, khu kinh tế Nghi Sơn, hệ thống cảng Lễ Môn, Quảng Châu và các tuyến đê biển xung yếu.
- Tác động chính sách & Quản trị công: Cung cấp cơ sở khoa học chính thống cho UBND tỉnh Thanh Hóa và Bộ Tài nguyên & Môi trường trong việc lập "Quy hoạch không gian biển tỉnh Thanh Hóa thời kỳ 2021 - 2030, tầm nhìn đến năm 2050" và xây dựng hành lang bảo vệ bờ biển chống xâm nhập mặn.
- Lợi ích xã hội: Nâng cao khả năng thích ứng cho hơn 500.000 cư dân ven biển vùng hạ lưu sông Mã trước hiểm họa thiên tai bão lũ, triều cường, sạt lở bờ sông, bờ biển và suy thoái nguồn nước ngầm.
Đối tượng hưởng lợi
- Nghiên cứu sinh & Nhà khoa học trẻ: Tiếp cận một mô hình nghiên cứu mẫu mực về tích hợp địa tầng phân tập, thạch học định lượng và vi cổ sinh; nhận diện rõ các khoảng trống nghiên cứu về trầm tích Đệ tứ thềm lục địa.
- Giảng viên & Giáo sư chuyên ngành Địa chất - Môi trường: Tài liệu giảng dạy chuyên sâu, cung cấp các lát cắt thực tế và tập dữ liệu phân tích chuẩn xác phục vụ đào tạo sau đại học ngành Địa chất học, Trầm tích học và Địa mạo ven bờ.
- Doanh nghiệp tư vấn Địa kỹ thuật & Xây dựng: Nắm vững cấu trúc địa tầng dưới sâu (độ sâu tầng cát cuội sỏi chịu tải, phân bố các túi bùn sét yếu $lbQ_2^{2-3}$) để tối ưu hóa thiết kế nền móng công trình và định tuyến giao thông ven biển.
- Các nhà hoạch định chính sách & Quản lý đới bờ: Nhận được các khuyến nghị có cơ sở khoa học rõ ràng nhằm phân bổ không gian sử dụng đất hợp lý, kết hợp giữa bảo tồn hệ sinh thái bãi triều ngập mặn với phát triển kinh tế biển bền vững.
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 độc đáo nhất là việc xác lập và luận chứng thành công mô hình tiến hóa trầm tích của một châu thổ do sóng thống trị (wave-dominated delta) trong điều kiện bồn trũng ven biển hẹp chịu sự khống chế của đứt gãy tân kiến tạo sông Mã. Luận án đã mở rộng Lý thuyết địa tầng phân tập của Catuneanu (2006) và Mô hình phân loại 6 môi trường đới bờ của Boyd et al. (2006) bằng cách chứng minh rằng: Trong giai đoạn biển tiến (TST), năng lượng triều và vũng vịnh chiếm ưu thế lấp đầy thung lũng cắt xẻ; nhưng khi chuyển sang giai đoạn biển cao (HST), động lực sóng biển Đông trở thành nhân tố quyết định tái cấu trúc toàn bộ trầm tích tiền châu thổ thành các phức hệ cồn cát chắn bờ và đầm lầy sau cồn, tạo nên kiểu đồng bằng châu thổ bất đối xứng đặc trưng.
2. Sự đổi mới về phương pháp luận của luận án thể hiện qua việc so sánh với các nghiên cứu trước đây như thế nào?
So với các nghiên cứu trước đây tại đồng bằng Thanh Hóa (Đặng Trần Quân, 1980; Hoàng Văn Khổn, 1997; Nguyễn Ngọc Mên, 1988) vốn chủ yếu dựa vào mô tả thạch học định tính và đo vẽ địa mạo bề mặt:
- Luận án đã tiến hành chuẩn hóa địa tầng phân tập phân giải cao, sử dụng máy đo địa chấn khúc xạ hiện đại 24 kênh ABEM TerralogPro kết hợp 11,7 km địa chấn nông phản xạ phân giải cao để định vị chính xác bề mặt bóc mòn cơ sở thung lũng cắt xẻ (SB/TS).
- Luận án áp dụng phân tích đa chỉ số định lượng đồng bộ trên 9 lỗ khoan nguyên dạng (213m lõi) gồm: 531 mẫu độ hạt Folk, 21 tuổi tuyệt đối AMS $^{14}\text{C}$ (DirectAMS, Hoa Kỳ), 186 mẫu Diatomeae, 140 mẫu Foraminifera và các chỉ số địa hóa định lượng ($Kt$, $Fe^{2+}S/C{org}$), tạo nên độ tin cậy vượt trội so với các công trình đơn lẻ trước đó.
3. Phát hiện nào trong luận án gây bất ngờ nhất và có bằng chứng thực nghiệm nào chứng minh?
Phát hiện bất ngờ nhất là sự tồn tại của một thung lũng cắt xẻ (incised valley) cổ sâu tới -32m dưới mực nước biển tại khu vực LKTH9 (cửa Lèn - Hậu Lộc), được lấp đầy bởi tầng cát sạn sỏi hạt thô dày gần 4m có tuổi chuyển tiếp Pleistocen muộn - Holocen sớm. Bằng chứng thực nghiệm gồm:
- Kết quả đo địa chấn khúc xạ tuyến DVL-21 ghi nhận mặt phản xạ sóng âm biến đổi vận tốc đột ngột tại độ sâu -31,5m đến -32m.
- Mẫu thạch học lõi khoan LKTH9 tại độ sâu 31,9 - 32,0m phát hiện cuội sạn thạch anh mài tròn kém, gắn kết bở rời, kích thước hạt trung vị $Md > 2\text{mm}$, nằm trực tiếp trên nền sét loang lổ phong hóa chứa kết vón laterit của hệ tầng Vĩnh Phúc ($Q_1^3 vp$).
4. Luận án có cung cấp quy trình tái lập (replication protocol) chi tiết không?
Luận án cung cấp một quy trình thực nghiệm chi tiết và minh bạch hoàn toàn:
- Ghi rõ tọa độ không gian chính xác của 9 lỗ khoan máy và các tuyến địa chấn.
- Mô tả chi tiết quy trình chuẩn bị mẫu độ hạt theo phương pháp sàng rây và lắng đọng pipette; quy trình gia công 25 lát mỏng thạch học bở rời; thông số máy nhiễu xạ XRD (PANalytical, ống phóng CuKα, bước nhảy $0,02^\circ$).
- Công thức tính chỉ số kation trao đổi $Kt$ và tỷ số $Fe^{2+}S/C{org}$.
- Phần mềm hiệu chỉnh tuổi carbon phóng xạ (Calib Rev 8.0) và phòng thí nghiệm thực hiện (DirectAMS, USA), cho phép các nhà nghiên cứu độc lập tái lập hoàn toàn quy trình phân tích.
5. Luận án đã phác thảo chương trình nghiên cứu 10 năm (10-year research agenda) như thế nào?
Luận án định hướng chương trình nghiên cứu dài hạn tập trung vào:
- Xác lập mối liên kết địa tầng - trầm tích giữa châu thổ sông Mã và hệ thống bồn trũng Sông Hồng ngoài khơi vịnh Bắc Bộ.
- Ứng dụng công nghệ mô phỏng số 3D đánh giá tương tác sóng - dòng triều - bùn cát trong điều kiện kịch bản biến đổi khí hậu cực đoan.
- Ứng dụng kỹ thuật đồng vị bền để giải mã các biến cố khí hậu đột biến trong Holocen.
- Xây dựng bản đồ phân vùng biến động động lực hình thái đới bờ phục vụ bảo vệ hệ sinh thái rừng ngập mặn và quy hoạch hạ tầng cảng biển nước sâu.
Kết luận
- Xác lập khung phân loại 23 tướng trầm tích đơn lẻ thuộc 4 nhóm tướng: Luận án đã phân chia chi tiết và làm sáng tỏ đặc điểm thạch học, cổ sinh, hóa - lý môi trường của nhóm tướng aluvi lấp đầy thung lũng cắt xẻ (3 tướng), nhóm tướng estuary - vũng vịnh (5 tướng), nhóm tướng châu thổ (10 tướng) và nhóm tướng aluvi hiện đại (5 tướng).
- Khôi phục 4 giai đoạn tiến hóa trầm tích Holocen gắn liền với dao động mực nước biển: (1) Giai đoạn bồi lấp thung lũng cắt xẻ ($Q_2^1$); (2) Giai đoạn estuary - vũng vịnh ($Q_2^{1-2}$); (3) Giai đoạn châu thổ ($Q_2^{2-3}$); và (4) Giai đoạn bồi tích aluvi hiện đại ($Q_2^3$).
- Làm sáng tỏ quy luật dao động mực nước biển Flandrian tại khu vực ven biển Thanh Hóa: Biển tiến bắt đầu dâng ngập thung lũng cổ từ 10.000 năm Bp với tốc độ 9 - 12 mm/năm, đạt đỉnh ngập lụt cực đại (MFS) ở độ cao +4m đến +5m vào khoảng 6.000 năm Bp, sau đó hạ dần theo quy luật dao động tắt dần tạo nên các bậc thềm và đường bờ cổ (giai đoạn Đa Bút 6.000 năm Bp và Hoa Lộc 4.000 năm Bp).
- Chứng minh vai trò thống trị của động lực sóng trong Holocen muộn: Sự tương tác giữa dòng phù sa sông Mã và trường sóng Đông Bắc đã hình thành hệ thống đê cát, cồn cát chắn bờ, đồng bằng lagoon và cấu trúc châu thổ bất đối xứng.
- Mở ra 3 hướng nghiên cứu học thuật mới: Địa tầng phân tập phân giải cao vùng chuyển tiếp thềm lục địa; cổ khí hậu học phân tử Đệ tứ; và mô hình hóa động lực học trầm tích ven bờ dưới tác động của các hồ chứa thượng nguồn.
- Đóng góp giá trị thực tiễn to lớn: Cung cấp cơ sở khoa học chuẩn xác phục vụ quy hoạch lãnh thổ bền vững, định tuyến đường bộ ven biển 107,5 km, phát triển các cảng biển và đề xuất giải pháp công trình - sinh học ứng phó với xói lở bờ biển và biến đổi khí hậu toàn cầu.
