Tổng quan về luận án

Nghiên cứu động lực học trầm tích tại các vùng đất ngập nước ven biển, đặc biệt là hệ sinh thái rừng ngập mặn (RNM), đóng vai trò cốt lõi trong bảo vệ bờ biển, thích ứng biến đổi khí hậu và duy trì cân bằng sinh thái địa mạo. Luận án tiến sĩ ngành Vật lý địa cầu (Mã số: 62440111) của nghiên cứu sinh Lê Nguyễn Hoa Tiên, dưới sự hướng dẫn khoa học của PGS. Võ Lương Hồng Phước tại Trường Đại học Khoa học Tự nhiên – Đại học Quốc gia TP. Hồ Chí Minh (2023), mang tiêu đề "Nghiên cứu động lực học trầm tích trong vùng rừng ngập mặn", đại diện cho một bước tiến thực nghiệm và mô hình hóa toán học chuyên sâu trong hải dương học vật lý cận duyên.

                                  ĐỘNG LỰC HỌC SÓNG & TRẦM TÍCH RNM
      [Động lực học sóng biển]                                        [Động lực học trầm tích lơ lửng]
  - Phổ sóng (Blackman-Tukey, FFT)                                - Thí nghiệm cột chìm lắng (Settling column)
  - Hệ số suy giảm năng lượng r                                   - Xác định vận tốc lắng w_s (flocs kết dính)
  - Mô phỏng tiêu tán bằng WAPROMAN                               - Mô hình 1D SSC không dừng có tái lơ lửng
                               [ĐỐI SÁNH ĐỊA HÌNH ĐỘNG LỰC HỌC]
                    - Cần Giờ (Nàng Hai): Cửa sông Đồng Tranh, xói lở cục bộ
                    - Cù Lao Dung (Sóc Trăng): Cửa sông Hậu (Định An - Trần Đề), bồi tụ

Thực tế nghiên cứu tại Việt Nam cho thấy phần lớn các công trình trước đây chỉ tập trung khảo sát định tính hoặc đánh giá thực nghiệm suy giảm sóng tại các bãi bồi ven biển miền Bắc (Thái Bình, Nam Định, Hải Phòng). Miền Nam Việt Nam – nơi tập trung hai hệ thống rừng ngập mặn tiêu biểu là Khu dự trữ sinh quyển Cần Giờ (tiểu khu IVA) và Cù Lao Dung - Sóc Trăng thuộc Đồng bằng sông Cửu Long (tiểu khu IVB) – đang chứng kiến sự phân hóa địa mạo phức tạp: "Mặc dù cả hai đều thuộc khu vực ven biển Nam Bộ, nhưng quá trình xói lở và bồi tụ bờ biển là khác nhau: Cù Lao Dung có xu hướng bồi lắng trong khi khu vực Nàng Hai thuộc RNM Cần Giờ bị xói lở." Luận án giải quyết triệt để research gap về cơ chế tương tác vi mô và vĩ mô giữa sóng gió, dòng triều bán nhật triều không đều và trầm tích kết dính (cohesive sediments), khắc phục sự thiếu hụt các tham số vật lý then chốt như vận tốc lắng đọng ($w_s$) và hệ số khuếch tán thẳng đứng ($K_z$) trong các mô hình vận chuyển bùn cát hiện hữu.

Câu hỏi nghiên cứu và giả thuyết khoa học được chuẩn hóa:

  • RQ1: Năng lượng và phổ sóng suy giảm như thế nào qua các cấu trúc rễ - thân cây ngập mặn đặc trưng (Rhizophora apiculata, Avicennia alba, Sonneratia alba) trong các pha triều và mùa gió đối nghịch?
  • H1: Sự suy giảm độ cao sóng có nghĩa ($H_s$) phụ thuộc phi tuyến vào độ sâu cột nước ngập, mật độ rễ thở/rễ chống và hiện tượng sóng bể vỡ cục bộ.
  • RQ2: Vận tốc lắng đọng ($w_s$) của trầm tích kết dính trong RNM biến thiên ra sao theo hàm lượng trầm tích lơ lửng (SSC) và độ mặn?
  • H2: Trầm tích lơ lửng trong RNM tồn tại chủ yếu ở dạng khối kết dính (flocs), có quy luật chìm lắng phi tuyến chuyển tiếp qua 3 pha: chìm lắng tự do, chìm lắng kết tụ và chìm lắng cản trở.
  • RQ3: Quy luật phân bố SSC theo độ sâu được thiết lập như thế nào khi xét đồng thời quá trình khuếch tán rối, lắng đọng và tái lơ lửng từ đáy?
  • H3: Ứng suất đáy do sóng kết hợp dòng triều đóng vai trò chi phối cơ chế tái lơ lửng (resuspension), quyết định sự biến thiên nồng độ SSC trong tầng sát đáy.

Về quy mô, nghiên cứu tích hợp chuỗi số liệu đo đạc hải văn – khí tượng đồng bộ qua nhiều chiến dịch thực địa: Cần Giờ (tháng 02/2012 - gió mùa Đông Bắc, tháng 06/2014 - gió mùa Tây Nam) và Cù Lao Dung (tháng 09-10/2014 - gió mùa Tây Nam, tháng 03/2015 - gió mùa Đông Bắc), kết hợp phân tích chuỗi thời gian khí tượng - mực nước 20 năm (1987-2006) tại trạm kiểm chứng Vũng Tàu và trạm Trần Đề (2018).

Literature Review và Positioning

Khung lý thuyết về động lực học trầm tích vùng triều ngập mặn được định hình qua nhiều trường phái học thuật quốc tế. Woodroffe (1992) đã phân loại 5 kiểu môi trường RNM dựa trên tương quan lực chi phối: sông chiếm ưu thế (river-dominated), triều chiếm ưu thế (tide-dominated), sóng chiếm ưu thế (wave-dominated), hỗn hợp sông - sóng và thung lũng đá ngập chìm. Trong khi đó, Van Rijn (1993) chuẩn hóa ba cơ chế vận chuyển trầm tích: tải đáy (bed load), tải lơ lửng (suspended load) và tải rửa trôi (wash load). Đối với môi trường trầm tích hạt mịn kết dính, Krone (1962), Partheniades (1965) và Mehta (1989) đã đặt nền móng cho động lực học bùn sét thông qua tương tác điện hóa bề mặt, tạo ra các chu trình kết tụ (aggregation), phá vỡ (breakup), lắng đọng (deposition) và xói mòn (erosion).

Trong y văn tồn tại hai luồng quan điểm tranh luận sâu sắc:

  1. Trường phái thủy triều chiếm ưu thế (Tide-dominated trapping): Wolanski và cộng sự (1992, 1995) qua các nghiên cứu tại rạch Middle và rạch Cocoa (Townsville, Australia) khẳng định dòng triều và tính chất bất đối xứng của dòng triều là nhân tố quyết định đưa bùn cát vào rừng; hiện tượng nước đứng (slack water) kết hợp bẫy cơ học của thảm rễ làm lắng đọng 80% trầm tích lơ lửng, gia tăng bồi tụ 0,1 cm/năm.
  2. Trường phái sóng chiếm ưu thế và tái lơ lửng (Wave-induced resuspension): Sanford (1994) tại Upper Chesapeake Bay (Mỹ), Christie và cs (1999) tại cửa sông Humber (Anh), và Green và cs (1997) tại cảng Manukau (New Zealand) chỉ ra rằng sóng gió chu kỳ ngắn dù biên độ nhỏ (<0,2 m) vẫn tạo ứng suất tiếp đáy $\tau_b$ vượt ngưỡng xói tới hạn ($\tau_{cr}$), khiến sự tái lơ lửng do sóng cao gấp 3-5 lần so với tác động thuần túy của triều.

So sánh với các nghiên cứu quốc tế tiêu biểu:

  • Tại bờ biển Andaman, tỉnh Trang, Thái Lan (Bao, 2011; Mazda et al., 2006), năng lượng sóng gió ($T < 10$ s) giảm 63% - 72%, đạt cực đại suy giảm 79% trên cự ly 141 m trong điều kiện bão.
  • Tại Vịnh Leizhou, Trung Quốc (Chen et al., 2018), RNM tiêu tán năng lượng sóng vượt trội so với bãi triều trần nhờ cấu trúc tán và rễ ngập triều, đặc biệt dưới tác động sóng bão cực đoan.

Công trình của Lê Nguyễn Hoa Tiên định vị chính xác khoảng trống học thuật: Trong khi các mô hình toàn cầu như SWAN (sử dụng công thức Dalrymple et al., 1984), DELFT-3D hay FVCOM-ESSed (sử dụng tại cảng Darwin - Úc) gặp sai số lớn khi áp dụng vào vùng bãi triều có rừng ngập mặn Đông Nam Á do chưa phản ánh đúng đặc tính kết tụ hạt sét nhiệt đới và dạng hình học rễ cây (Sonneratia, Avicennia, Rhizophora), luận án đã kết nối hoàn chỉnh giữa dữ liệu đo đạc vi thủy động lực học, thực nghiệm chìm lắng chuyên biệt và mô phỏng 1D động lực trầm tích có xét đến nguồn tái lơ lửng từ đáy.

Đó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 các lý thuyết thủy động lực học công trình biển và địa động lực học trầm tích:

  1. Phát triển lý thuyết tiêu tán sóng phi tuyến của Massel et al. (1999) và Võ Lương Hồng Phước (2008): Bổ sung cơ chế tiêu tán kết hợp giữa lực cản hình học của thảm thực vật ($C_D$) và sự tiêu tán do sóng vỡ (wave breaking) khi độ sâu nước ($h$) biến thiên theo pha triều. Minh chứng thực nghiệm khẳng định: "Sự suy giảm sóng trong RNM chủ yếu là do sự tổn thất năng lượng do lực cản của thảm thực vật, vì sự suy giảm sóng do ma sát đáy và sự tiêu tán tính nhớt trên bãi triều thấp hơn đáng kể so với bên trong thảm thực vật ngập mặn."
  2. Xác lập mô hình lắng đọng trầm tích kết dính 3 pha theo lý thuyết Mehta (1989) và Winterwerp (2002): Luận án chứng minh trầm tích lơ lửng tại Nam Bộ tuân thủ nghiêm ngặt chu trình tạo bông (flocculation): "Đối với trầm tích trong RNM, 99% trầm tích lơ lửng sẽ bị dính kết tạo thành các khối kết dính... kích thước của flocs thay đổi từ 30 đến 300 $\mu$m. Kích thước trung bình của flocs là 100 $\mu$m tương ứng với vận tốc lắng đọng là 0,005 m/s, vận tốc lắng đọng này lớn gấp 100 lần vận tốc của từng hạt sét và bùn riêng cấu tạo nên các flocs."
  3. Hình thành khung mệnh đề lý thuyết động lực học trầm tích:
    • Mệnh đề 1: Tốc độ tiêu tán năng lượng sóng tỷ lệ nghịch với chiều sâu ngập nước nhưng tỷ lệ thuận với mật độ diện tích cản của thân và rễ thở.
    • Mệnh đề 2: Trong thảm rừng ngập mặn mật độ cao, hệ số khuếch tán rối thẳng đứng tầng đáy ($K_{z,bottom}$) suy giảm mạnh, tạo điều kiện giữ lại các hạt sét mịn ($<4$ $\mu$m) và bùn mịn ($4 - 63$ $\mu$m).
    • Mệnh đề 3: Tại ranh giới ngoài bãi triều tiếp giáp biển, tương tác phi tuyến sóng - dòng triều tạo ra vùng cực đại SSC do hiện tượng tái lơ lửng tức thời (instantaneous bottom resuspension).

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

Khung phân tích của luận án tích hợp đa tầng 3 trục lý thuyết: Lý thuyết lan truyền sóng nước nông phi tuyến, Động lực học trầm tích kết dính (Cohesive Sediment Dynamics) và Thủy khí động lực học môi trường cản xốp (Flow-vegetation porous media dynamics).

                      KHUNG PHÂN TÍCH TỔNG HỢP CỦA LUẬN ÁN

Phương trình chủ đạo mô phỏng phân bố nồng độ SSC ($C$) theo phương thẳng đứng ($z$) và thời gian ($t$) được thiết lập:

$$\frac{\partial C}{\partial t} = \frac{\partial}{\partial z}\left(K_z \frac{\partial C}{\partial z}\right) + w_s \frac{\partial C}{\partial z}$$

Điều kiện biên mặt ($z = h$):

$$K_z \frac{\partial C}{\partial z} + w_s C = 0$$

Điều kiện biên đáy ($z = 0$):

$$K_z \frac{\partial C}{\partial z} + w_s C = E - D$$

Trong đó $E$ là thông lượng tái lơ lửng khi ứng suất đáy vượt ngưỡng ($\tau_b > \tau_{ce}$), $D$ là thông lượng lắng đọng khi $\tau_b < \tau_{cd}$. Khung phân tích xác lập rõ boundary conditions: không xét nguồn bùn cát nội lưu sâu từ thượng nguồn sông chuyển vận dạng trôi nổi tự do trong mùa lũ cực đại, tập trung vào phạm vi đới bãi triều ngập triều định kỳ.

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

Thiết kế nghiên cứu

Luận án vận dụng triết lý nghiên cứu thực chứng (positivism) kết hợp suy luận diễn dịch - thực nghiệm (deductive-empirical approach). Thiết kế nghiên cứu đa cấp độ (multi-level design) kết hợp đồng thời 4 trụ cột: quan trắc hải văn in-situ độ phân giải cao, phân tích quang phổ năng lượng sóng, thực nghiệm cột lắng thủy lực trong phòng thí nghiệm và mô phỏng số giải tích.

Hệ thống trạm đo được bố trí dọc theo các transect thẳng góc với đường bờ:

  • Khu vực Cần Giờ: Tuyến đo gồm 3 trạm đại diện: Trạm cửa sông ngoài khơi (ST0-CG), Trạm bãi triều trần (ST1-CG) và Trạm nằm sâu trong thảm rừng ngập mặn Nàng Hai (ST2-CG).
  • Khu vực Cù Lao Dung: Tuyến đo gồm 4 trạm phân cấp: Trạm biển ngoài khơi (ST0-CLD), Trạm bãi triều phía ngoài (ST1-CLD), Trạm rìa rừng ngập mặn (ST2-CLD) và Trạm nằm sâu trong lõi rừng ngập mặn (ST3-CLD).

Quy trình nghiên cứu rigorous

Thiết bị đo đạc và giao thức chuẩn hóa:

  • Đo trường sóng và dao động mực nước: Thiết bị đo áp suất sóng độ chính xác cao Wave Gauge kết hợp cảm biến áp suất tự ghi RBR-Virtuoso. Tần số lấy mẫu đạt 2 Hz - 4 Hz cho phép trích xuất đầy đủ các sóng gió chu kỳ ngắn ($T = 1 - 10$ s).
  • Đo trường dòng chảy và SSC: Cảm biến quang học OBS-3+ (Optical Backscatter Sensors) đo độ đục (NTU) được hiệu chuẩn trực tiếp với mẫu nước lọc qua màng lọc sợi thủy tinh Whatman GF/C (kích thước lỗ $1,2$ $\mu$m) để chuyển đổi chính xác sang hàm lượng trầm tích lơ lửng ($mg/l$). Vận tốc dòng chảy được đo bằng máy đo lưu tốc Doppler âm thanh (Acoustic Doppler Velocimeter - ADV) và ADCP SonTek 2.0 MHz.
  • Đo tham số môi trường nước: Thiết bị CTD (Conductivity-Temperature-Depth) đo độ dẫn điện, nhiệt độ và độ mặn liên tục trong cột nước.
  • Quan trắc biến động vi địa mạo: Kỹ thuật cọc đánh dấu (tracer sticks) kết hợp đo đạc mặt cắt địa hình đáy định kỳ.
                  QUY TRÌNH THỰC NGHIỆM CỘT CHÌM LẮNG TRẦM TÍCH

Thí nghiệm cột chìm lắng trầm tích (Settling Column Test) trong phòng thí nghiệm được tiến hành với cột trụ plexiglass chuyên dụng có chiều cao $1,8$ m, đường kính trong $0,15$ m, bố trí các cổng hút mẫu tại 3 tầng độ sâu khác nhau ($z = 0,3$ m, $0,8$ m, $1,3$ m). Mẫu trầm tích nguyên dạng lấy từ bãi triều và trong rừng tại Cần Giờ và Cù Lao Dung được khuấy trộn đều với nước mặn hiện trường, sau đó rút mẫu định kỳ tại các mốc thời gian: $0, 5, 10, 15, 20, 30, 45, 60, 90, 120, 180, 240, 360, 1440$ phút nhằm xác định chính xác đường cong bán thực nghiệm lắng đọng.

Phổ năng lượng sóng $S(f)$ được tính toán và kiểm chứng chéo (cross-validation) bằng hai phương pháp toán học độc lập:

  1. Phương pháp Blackman-Tukey: Tính hàm tự tương quan $R(\tau)$ qua phép dịch bước $m$ ($m = 32, 48, 64, 80, 96, 112, 128$), sau đó nhân cửa sổ làm mịn Hamming/Tukey và biến đổi Fourier.
  2. Phương pháp Fast Fourier Transform (FFT): Chia nhỏ chuỗi thời gian thành các đoạn độ dài $N = 512, 1024$ điểm với độ chồng lấp $50%$ để giảm thiểu phương sai phổ.

Data và phân tích

Phân tích cấp phối hạt theo thang Wentworth (1922) trên máy phân tích kích thước hạt laser (Laser Diffraction Particle Size Analyzer) ghi nhận:

  • Mẫu trầm tích Cần Giờ: Tỷ lệ bùn ($4 - 63$ $\mu$m) chiếm $58,4% - 64,2%$, sét ($<4$ $\mu$m) chiếm $31,5% - 38,1%$, cát mịn ($>63$ $\mu$m) chiếm $<3,5%$.
  • Mẫu trầm tích Cù Lao Dung: Trầm tích đáy vùng bãi triều ngoài có đường kính trung vị $D_{50} = 15 - 28$ $\mu$m, sâu trong rừng $D_{50} = 6,5 - 11,2$ $\mu$m; tại các luồng cửa sông Định An trầm tích thô hơn ($D_{50} = 112 - 310$ $\mu$m) do cát từ sông Hậu đổ ra.

Dữ liệu mô phỏng sóng được tích hợp vào mô hình giải tích WAPROMAN (Wave Propagation in Mangrove Forests). Các tham số thực vật được số hóa chính xác:

  • Loài Avicennia alba (Mắm trắng): Mật độ rễ thở $N_r = 180 - 320$ rễ/m², đường kính trung bình $d_r = 0,008 - 0,012$ m, chiều cao rễ $h_r = 0,15 - 0,25$ m.
  • Loài Rhizophora apiculata (Đước đôi): Mật độ rễ chống $N_r = 45 - 85$ rễ/m², đường kính $d_r = 0,035 - 0,065$ m, chiều cao tương tác lên đến $1,2$ m.

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

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

      PHỔ NĂNG LƯỢNG SÓNG TỪ NGOÀI KHƠI VÀO SÂU TRONG RỪNG NGẬP MẶN
  Mật độ phổ S(f)
  1. Hiệu ứng suy giảm sóng phi đối xứng và tiêu tán chọn lọc theo tần số:
    • Dữ liệu đo đạc tại Cù Lao Dung qua hai mùa gió chứng minh năng lượng sóng giảm mạnh từ $65,8%$ đến $84,2%$ khi lan truyền qua dải rừng ngập mặn rộng $200 - 300$ m. Hệ số suy giảm sóng thực nghiệm ($r$) trong mùa gió Tây Nam dao động từ $0,0035$ đến $0,0058$ m⁻¹, trong khi mùa gió Đông Bắc đạt $0,0042 - 0,0071$ m⁻¹.
    • Phân tích phổ năng lượng sóng khẳng định rừng ngập mặn hoạt động như một bộ lọc thông thấp (low-pass filter): Các đỉnh sóng gió tần số cao ($f > 0,2$ Hz, chu kỳ $T < 5$ s) bị triệt tiêu gần như hoàn toàn ($>90%$), trong khi các thành phần sóng chu kỳ dài hoặc dao động triều duy trì khả năng truyền sâu.
  2. Quy luật vận tốc lắng đọng ($w_s$) phi tuyến của trầm tích kết dính:
    • Từ thí nghiệm cột lắng, phương trình bán thực nghiệm xác định vận tốc lắng đọng theo nồng độ SSC ($C$) được thiết lập:

$$w_s = a \cdot C^m \quad (\text{khi } C < C_{crit}) \quad \text{và} \quad w_s = a \cdot C^m \left(1 - b \cdot C\right)^n \quad (\text{khi } C \ge C_{crit})$$

  • Tại RNM Cần Giờ: Bộ tham số hiệu chỉnh đạt $a = 0,0012$, $m = 0,85$, $b = 0,00045$, $n = 1,5$, cho vận tốc lắng đọng cực đại $w_{s,max} = 0,48 \times 10^{-3} - 1,06 \times 10^{-3}$ m/s ở nồng độ bùn cát tối ưu $C \approx 1,2 - 1,8$ g/l.
  • Tại RNM Cù Lao Dung: Do độ mặn thấp hơn và tỷ lệ khoáng sét kaolinite/illite khác biệt, $w_{s,max}$ đạt $0,35 \times 10^{-3} - 0,82 \times 10^{-3}$ m/s. Khi nồng độ $C > 2,5$ g/l, hiện tượng chìm lắng cản trở xuất hiện khiến $w_s$ sụt giảm nhanh chóng.
  1. Cơ chế tái lơ lửng trầm tích tầng đáy do tương tác sóng - dòng triều:
    • Tại trạm bãi triều Cần Giờ (ST1-CG), SSC không suy giảm đơn điệu theo triều rút mà ghi nhận các xung nồng độ cực đại ($SSC > 350$ mg/l) xuất hiện đồng thời với thời điểm độ cao sóng $H_s$ tăng cục bộ kết hợp dòng triều rút mạnh. Điều này bác bỏ giả định truyền thống cho rằng trầm tích chỉ lắng đọng tĩnh trong chu kỳ triều.
  2. Sự tương phản động lực học bồi - xói giữa hai hệ sinh thái:
    • Tại Cù Lao Dung, lưu lượng trầm tích qua cửa Định An chiếm ưu thế ($60% - 70%$ tổng lượng bùn cát sông Hậu) kết hợp biên độ triều $1,5 - 2,0$ m tạo điều kiện cho thảm thực vật Avicennia alba tiên phong mở rộng bãi bồi ra biển với tốc độ bồi lắng trung bình $1,2 - 2,5$ cm/năm.
    • Tại Cần Giờ (vùng Nàng Hai), sự thiếu hụt nguồn phù sa bổ sung từ sông Đồng Tranh cùng với sự biến dạng của trường sóng khúc xạ qua luồng tàu Soài Rạp đã biến khu vực này thành điểm nóng xói lở bờ biển.

Implications đa chiều

  • Về mặt học thuật và lý thuyết: Luận án cung cấp bộ số liệu thực nghiệm chuẩn mực về các hệ số vật lý địa cầu ($w_s, K_z, r, C_D$) đặc thù cho vùng ven biển nhiệt đới gió mùa Đông Nam Á, tạo cơ sở hiệu chỉnh các mô hình thủy động lực học 2D/3D (như MIKE 21, DELFT-3D, TELEMAC).
  • Về mặt kỹ thuật công trình: Cung cấp công thức định lượng chiều rộng đai rừng an toàn tối thiểu ($B_{min} \ge 250 - 350$ m) phục vụ thiết kế đê biển sinh thái kết hợp thảm cây ngập mặn giảm sóng thay thế kè cứng bê tông.
  • Về mặt quản trị tài nguyên và chính sách: Cung cấp luận cứ khoa học để Sở Nông nghiệp & PTNT TP.HCM và Sóc Trăng quy hoạch phân vùng bảo tồn lõi rừng ngập mặn, thiết lập các giải pháp công trình mỏ hàn mềm giữ bùn giảm sóng tại các vị trí xói lở trọng điểm.

Limitations và Future Research

Nghiên cứu duy trì góc nhìn khoa học khách quan khi thừa nhận các giới hạn nội tại:

  1. Giới hạn không gian mô hình 1D: Mô hình tính toán phân bố SSC theo phương thẳng đứng (1D vertical) chưa mô phỏng trọn vẹn hoàn lưu ngang 2D và các dòng chảy xoáy cục bộ phức tạp quanh từng gốc cây đơn lẻ trong các lạch triều uốn khúc.
  2. Độ phân giải thời gian của các đợt đo đạc bão: Chuỗi số liệu đo đạc hiện trường tuy bao quát mùa gió Đông Bắc và Tây Nam nhưng chưa đo đạc trực tiếp trong tâm các cơn bão nhiệt đới cực đoan (cấp 11-12) do rủi ro an toàn thiết bị đo in-situ.
  3. Giả định về tính đồng nhất của hạt flocs: Mô hình hóa vận tốc lắng $w_s$ giả định kích thước hạt flocs trung bình đại diện, chưa tích hợp mô hình biến thiên động kích thước hạt bông bùn theo gradient ứng suất cắt rối ($\tau_{turb}$).

Chương trình nghiên cứu tiếp nối (Future Agenda):

  • Phát triển mô hình ghép nối 3D Hydro-Sediment-Vegetation tích hợp trực tiếp phương trình động lực học phân rã - kết tụ bông bùn (Flocculation Population Balance Equation).
  • Triển khai hệ thống trạm quan trắc hải văn ngầm tự động thời gian thực (Real-time oceanographic seabed nodes) đo đạc liên tục qua nhiều chu kỳ năm El Niño/La Niña.
  • Đánh giá vai trò của chất hữu cơ hòa tan (Extracellular Polymeric Substances - EPS) do thảm vi sinh vật đáy tiết ra trong việc gia tăng ứng suất xói tới hạn của bùn kết dính.

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

                               TÁC ĐỘNG ĐA PHƯƠNG DIỆN

Công trình của Lê Nguyễn Hoa Tiên tạo ra tầm ảnh hưởng rộng lớn:

  • Tác động học thuật: Định hình tiêu chuẩn nghiên cứu động lực học trầm tích ven biển tại Việt Nam, thu hẹp khoảng cách với các nhóm nghiên cứu hàng đầu thế giới tại Hà Lan (TU Delft), Úc (AIMS) và Mỹ (USGS).
  • Ứng dụng thực tiễn: Các kết quả tính toán suy giảm sóng bằng WAPROMAN và bộ tham số lắng đọng trầm tích được chuyển giao trực tiếp cho các đề án phục hồi rừng ngập mặn phòng hộ ven biển Đồng bằng sông Cửu Long, phục vụ chương trình mục tiêu quốc gia ứng phó biến đổi khí hậu.
  • Bảo tồn sinh thái: Khẳng định giá trị dịch vụ hệ sinh thái của RNM Cần Giờ (75.740 ha) và Cù Lao Dung trong việc bẫy bùn cát, hấp thụ carbon xanh và ổn định dải ven bờ trước hiểm họa nước biển dâng.

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

  • Nghiên cứu sinh và học viên cao học: Tiếp cận quy trình phương pháp luận mẫu mực kết hợp đo đạc vi thủy động lực, xử lý phổ Fourier/Blackman-Tukey và thí nghiệm cột lắng bùn sét.
  • Các nhà khoa học Vật lý địa cầu / Hải dương học: Sở hữu bộ cơ sở dữ liệu thực nghiệm chuẩn xác về trường sóng, dòng chảy và SSC tại cửa sông Đồng Tranh và sông Hậu.
  • Kỹ sư công trình thủy lợi - công trình thủy: Sử dụng công thức thực nghiệm suy giảm sóng và hệ số cản $C_D$ để tính toán tải trọng sóng tác dụng lên đê biển sau thảm rừng.
  • Nhà hoạch định chính sách môi trường: Có căn cứ khoa học định lượng để phân bổ ngân sách bảo vệ rừng ngập mặn dựa trên hiệu quả giảm sóng và bồi tụ thực tế.

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?

Luận án đã mở rộng Lý thuyết tiêu tán năng lượng sóng trong thảm thực vật của Dalrymple et al. (1984) và Massel et al. (1999) bằng việc tích hợp hiệu ứng tương tác phi tuyến giữa hình học rễ cây phức tạp (Avicennia/Rhizophora) và sự tiêu tán do sóng bể vỡ trên nền đáy dốc biến đổi trong mô hình WAPROMAN. Đồng thời, luận án đã cụ thể hóa Lý thuyết chìm lắng trầm tích kết dính của Mehta (1989) bằng việc thiết lập chính xác các hệ số thực nghiệm ($a, b, m, n$) mô tả trọn vẹn cả 3 pha chìm lắng (tự do, tạo bông và cản trở) cho bùn sét ven biển Nam Bộ.

2. Đột phá phương pháp luận của luận án khi so sánh với các nghiên cứu trước đây?

So với công trình của Bao (2011) tại Thái Lan chỉ đo đạc phân tích thống kê sóng đơn thuần hoặc công trình của Mazda et al. (2006) tại Nhật Bản chỉ khảo sát trong kênh dẫn hẹp, luận án đã thực hiện bước đột phá: Kết hợp chuỗi đo đạc in-situ đa trạm đồng bộ ngoài khơi - bãi triều - trong rừng với thí nghiệm cột chìm lắng cao $1,8$ m và mô hình hóa 1D SSC có xét đến điều kiện biên tái lơ lửng đáy theo thời gian thực ($\partial C/\partial t$).

3. Phát hiện bất ngờ và phản trực giác nhất từ số liệu thực nghiệm là gì?

Phát hiện phản trực giác nhất là: Hàm lượng trầm tích lơ lửng (SSC) trong thảm rừng ngập mặn Cần Giờ không suy giảm tuyến tính khi triều rút mà ghi nhận hiện tượng tái lơ lửng tăng vọt cục bộ. Nguyên nhân do sự kết hợp đồng thời giữa sóng gió chu kỳ ngắn xâm nhập tầng nước nông và ứng suất cắt đáy gia tăng đột biến khi dòng triều đổi hướng, chứng minh thảm rừng không chỉ là bể lắng thụ động mà là một hệ động lực biến tính liên tục.

4. Giao thức lặp lại nghiên cứu (Replication Protocol) có được cung cấp đầy đủ không?

Có. Luận án mô tả chi tiết quy trình chế tạo và vận hành cột chìm lắng thủy lực, định dạng cấu trúc tệp dữ liệu đầu vào, các bước tiền xử lý lọc nhiễu tín hiệu OBS bằng mẫu cân sấy GF/C, thuật toán tính phổ Blackman-Tukey với các bậc dịch bước $m$ và mã nguồn các phương trình sai phân hữu hạn giải bài toán 1D khuếch tán - lắng đọng.

5. Luận án đã phác thảo chương trình nghiên cứu 10 năm tiếp theo như thế nào?

Luận án định hình lộ trình nghiên cứu 10 năm tập trung vào: Mở rộng mô hình WAPROMAN thành phần mềm mã nguồn mở 2D/3D mô phỏng lan truyền sóng và biến động địa mạo vùng cửa sông ngập mặn; thiết lập mạng lưới quan trắc vi địa mạo tự động tại toàn bộ các cửa sông Đồng bằng sông Cửu Long; và định lượng hóa mối liên hệ giữa động lực bùn cát với khả năng lưu trữ carbon hữu cơ (Blue Carbon sequestration) trong trầm tích rừng ngập mặn.

Kết luận

Công trình luận án tiến sĩ của Lê Nguyễn Hoa Tiên đã hoàn thành xuất sắc các mục tiêu khoa học đề ra, mang lại những giá trị cốt lõi:

  1. Xác lập định lượng hệ số suy giảm sóng ($r$): Chứng minh năng lượng sóng suy giảm từ $65,8%$ đến $84,2%$ qua $200 - 300$ m rừng ngập mặn, làm sáng tỏ vai trò lọc thông thấp triệt tiêu sóng gió tần số cao của thảm cây nhiệt đới.
  2. Khám phá bản chất động lực học hạt bùn kết dính: Xác định $99%$ trầm tích lơ lửng trong RNM tồn tại dưới dạng flocs kích thước $30 - 300$ $\mu$m với vận tốc lắng đọng thực nghiệm $w_s = 0,35 \times 10^{-3} - 1,06 \times 10^{-3}$ m/s, lớn gấp 100 lần hạt sét đơn lẻ.
  3. Hoàn thiện mô hình toán học WAPROMAN và mô hình 1D SSC: Giải thành công bài toán phân bố nồng độ trầm tích lơ lửng theo độ sâu có xét đến quá trình khuếch tán rối ($K_z$) và tái lơ lửng đáy ($E$).
  4. Giải mã tương phản địa mạo Cần Giờ - Cù Lao Dung: Cung cấp bức tranh toàn cảnh về cơ chế xói lở tại cửa sông Đồng Tranh và bồi tụ mở rộng bãi bồi tại cửa sông Hậu dưới tác động tổng hợp của sóng, triều và nguồn bùn cát thượng nguồn.
  5. Mở ra 3 hướng nghiên cứu mới: Động lực học bông bùn trong dòng rối thảm thực vật, mô hình hóa số trị 3D tương tác sóng - dòng - bùn cát - rễ cây, và tích hợp động lực trầm tích vào lượng giá sinh thái carbon xanh.
  6. Giá trị di sản khoa học: Đóng góp một hệ thống dữ liệu thực nghiệm và lý thuyết chuẩn mực cho ngành Vật lý địa cầu và Hải dương học Việt Nam, khẳng định vai trò bảo vệ bờ biển không thể thay thế của hệ sinh thái rừng ngập mặn trong kỷ nguyên biến đổi khí hậu toàn cầu.